The Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614 Mutant: Conformational Plasticity and Frustration-Driven Allostery as Energetic Drivers of Highly Transmissible Spike Variant

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern(VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model which can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the inter-protomer hinge modulate the inter-domain interactions, global mobility change and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the inter-protomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

A multiscale computational study of the conformation of the full-length intrinsically disordered protein MeCP2

The malfunction of the Methyl CpG binding protein 2 (MeCP2) is associated to the Rett syndrome, one of the most common causes of cognitive impairment in females. MeCP2 is an intrinsically disordered protein (IDP), making its experimental characterization a challenge. There is currently no structure available for the full-length MeCP2 in any of the databases, and only the structure of its MBD domain has been solved. We used this structure to build a full-length model of MeCP2 by completing the rest of the protein via ab initio modelling. Using a combination of all-atom and coarse-grained simulations, we characterized its structure and dynamics as well as the conformational space sampled by the ID and TRD domains in the absence of the rest of the protein. The present work is the first computational study of the full-length protein. Two main conformations were sampled in the coarse-grained simulations: a globular structure similar to the one observed in the all-atom force field and a two-globule conformation. Our all-atom model is in good agreement with the available experimental data, predicting amino acid W104 to be buried, amino acids R111 and R133 to be solvent accessible, and having 4.1% of α-helix content, compared to the 4% found experimentally. Finally, we compared the model predicted by AlphaFold to our Modeller model. The model was not stable in water and underwent further folding. Together, these simulations provide a detailed (if perhaps incomplete) conformational ensemble of the full-length MeCP2, which is compatible with experimental data and can be the basis of further studies, e.g., on mutants of the protein or its interactions with its biological partners.

Actin crosslinker competition and sorting drive emergent GUV size-dependent actin network architecture

The proteins that make up the actin cytoskeleton can self-assemble into a variety of structures. In vitro experiments and coarse-grained simulations have shown that the actin crosslinking proteins α-actinin and fascin segregate into distinct domains in single actin bundles with a molecular size-dependent competition-based mechanism. Here, by encapsulating actin, α-actinin, and fascin in giant unilamellar vesicles (GUVs), we show that physical confinement can cause these proteins to form much more complex structures, including rings and asters at GUV peripheries and centers; the prevalence of different structures depends on GUV size. Strikingly, we found that α-actinin and fascin self-sort into separate domains in the aster structures with actin bundles whose apparent stiffness depends on the ratio of the relative concentrations of α-actinin and fascin. The observed boundary-imposed effect on protein sorting may be a general mechanism for creating emergent structures in biopolymer networks with multiple crosslinkers.