The experimental and computational studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across variants of concern(VOCs), linking the effect of this mutation with the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies characterized the closed and open states of the B.1.1.7 (B.1.1.7) and B.1.351 (Beta) spike variants allowing for a more detailed atomistic characterization of the conformational landscapes and functional changes. In this study, we employed coarse-grained simulations of the SARS-CoV-2 spike variant trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration analysis of the conformational ensembles with collective dynamics and residue-based mutational scanning of protein stability, we determine protein stability hotspots and identify potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. Through mutational scanning of protein stability changes we quantify mutational adaptability of the S-G614, S-B.1.1.7 and S-B.1.351 variants in different functional forms. Using this analysis, we found a significant conformational and mutational plasticity of the open states for all studied variants. The results of this study suggest that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites, the inter-protomer hinges of functional motions and motions of the receptor-binding domain required for binding of the host cell receptor. The proposed mechanism of mutation-induced energetic frustration may result in the greater adaptability and the emergence of multiple conformational substates in the open form. This study also suggested functional relationships between mutation-induced modulation of protein dynamics, local frustration and allosteric regulation of the SARS-CoV-2 spike protein.
Functional Significance of Hierarchical Tiers in Carbonmonoxy Myoglobin: Conformational Substates and Transitions Studied by Conformational Flooding Simulations