Coarse-Grained and Atomistic Simulations of the Salt-Stable Cowpea Chlorotic Mottle Virus (SS-CCMV) Subunit 26–49: β-Barrel Stability of the Hexamer and Pentamer Geometries

Dynamic stability of salt stable cowpea chlorotic mottle virus capsid protein dimers and pentamers of dimers

Intermediates of the self-assembly process of the salt stable cowpea chlorotic mottle virus (ss-CCMV) capsid can be modelled atomistically on realistic computational timescales either by studying oligomers in equilibrium or by focusing on their dissociation instead of their association. Our previous studies showed that among the three possible dimer interfaces in the icosahedral capsid, two are thermodynamically relevant for capsid formation. The aim of the current study is to evaluate the relative structural stabilities of the three different ss-CCMV dimers and to find and understand the conditions that lead to their dissociation. Long timescale molecular dynamics simulations at 300 K of the various dimers and of the pentamer of dimers underscore the importance of large contact surfaces on stabilizing the capsid subunits within an oligomer. Simulations in implicit solvent show that at higher temperature (350 K), the N-terminal tails of the protein units act as tethers, delaying dissociation for all but the most stable interface. The pentamer of dimers is also found to be stable on long timescales at 300 K, with an inherent flexibility of the outer protein chains.

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 with 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 modeling. 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 transcriptional repression domain (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 a 4.1% α-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.

Structural-kinetic-thermodynamic relationships identified from physics-based molecular simulation models

Coarse-grained molecular simulation models have provided immense, often general, insight into the complex behavior of condensed-phase systems but suffer from a lost connection to the true dynamical properties of the underlying system. In general, the physics that is built into a model shapes the free-energy landscape, restricting the attainable static and kinetic properties. In this work, we perform a detailed investigation into the property interrelationships resulting from these restrictions, for a representative system of the helix-coil transition. Inspired by high-throughput studies, we systematically vary force-field parameters and monitor their structural, kinetic, and thermodynamic properties. The focus of our investigation is a simple coarse-grained model, which accurately represents the underlying structural ensemble, i.e., effectively avoids sterically-forbidden configurations. As a result of this built-in physics, we observe a rather large restriction in the topology of the networks characterizing the simulation kinetics. When screening across force-field parameters, we find that structurally accurate models also best reproduce the kinetics, suggesting structural-kinetic relationships for these models. Additionally, an investigation into thermodynamic properties reveals a link between the cooperativity of the transition and the network topology at a single reference temperature.

The Role of Conformational Entropy in the Determination of Structural-Kinetic Relationships for Helix-Coil Transitions

Coarse-grained molecular simulation models can provide significant insight into the complex behavior of protein systems, but suffer from an inherently distorted description of dynamical properties. We recently demonstrated that, for a heptapeptide of alanine residues, the structural and kinetic properties of a simulation model are linked in a rather simple way, given a certain level of physics present in the model. In this work, we extend these findings to a longer peptide, for which the representation of configuration space in terms of a full enumeration of sequences of helical/coil states along the peptide backbone is impractical. We verify the structural-kinetic relationships by scanning the parameter space of a simple native-biased model and then employ a distinct transferable model to validate and generalize the conclusions. Our results further demonstrate the validity of the previous findings, while clarifying the role of conformational entropy in the determination of the structural-kinetic relationships. More specifically, while the global, long timescale kinetic properties of a particular class of models with varying energetic parameters but approximately fixed conformational entropy are determined by the overarching structural features of the ensemble, a shift in these kinetic observables occurs for models with a distinct representation of steric interactions. At the same time, the relationship between structure and more local, faster kinetic properties is not affected by varying the conformational entropy of the model.

Folding and insertion thermodynamics of the transmembrane WALP peptide

The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA)n (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of mean force characterizing the peptide's insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum-in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Testing the transferability of a coarse-grained model to intrinsically disordered proteins

The coarse-grained PLUM model is shown to capture structural and dimerization behaviour of the intrinsically disordered biomineralisation peptide n16N.

More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.

Mechanical and assembly units of viral capsids identified via quasi-rigid domain decomposition

Key steps in a viral life-cycle, such as self-assembly of a protective protein container or in some cases also subsequent maturation events, are governed by the interplay of physico-chemical mechanisms involving various spatial and temporal scales. These salient aspects of a viral life cycle are hence well described and rationalised from a mesoscopic perspective. Accordingly, various experimental and computational efforts have been directed towards identifying the fundamental building blocks that are instrumental for the mechanical response, or constitute the assembly units, of a few specific viral shells. Motivated by these earlier studies we introduce and apply a general and efficient computational scheme for identifying the stable domains of a given viral capsid. The method is based on elastic network models and quasi-rigid domain decomposition. It is first applied to a heterogeneous set of well-characterized viruses (CCMV, MS2, STNV, STMV) for which the known mechanical or assembly domains are correctly identified. The validated method is next applied to other viral particles such as L-A, Pariacoto and polyoma viruses, whose fundamental functional domains are still unknown or debated and for which we formulate verifiable predictions. The numerical code implementing the domain decomposition strategy is made freely available.

The genetic material of viruses is packaged inside capsids constituted from a few tens to thousands of proteins. The latter can organize in multimers that serve as fundamental blocks for the viral shell assembly or that control the capsid conformational transitions and response to mechanical stress. In this work, we introduce and apply a computational scheme that identifies the fundamental protein blocks from the structural fluctuations of the capsids in thermal equilibrium. These can be derived from phenomenological elastic network models with minimal computational expenditure. Accordingly, the basic functional protein units of a capsid can be obtained from the sole input of the capsid crystal structure. The method is applied to a heterogeneous set of viruses of various size and geometries. These include well-characterised instances for validation purposes, as well as debated ones for which predictions are formulated.

Coarse-grain modelling of protein-protein interactions

Here, we review recent advances towards the modelling of protein-protein interactions (PPI) at the coarse-grained (CG) level, a technique that is now widely used to understand protein affinity, aggregation and self-assembly behaviour. PPI models of soluble proteins and membrane proteins are separately described, but we note the parallel development that is present in both research fields with three important themes: firstly, combining CG modelling with knowledge-based approaches to predict and refine protein-protein complexes; secondly, using physics-based CG models for de novo prediction of protein-protein complexes; and thirdly modelling of large scale protein aggregates.

Optimization of an elastic network augmented coarse grained model to study CCMV capsid deformation

The major protective coat of most viruses is a highly symmetric protein capsid that forms spontaneously from many copies of identical proteins. Structural and mechanical properties of such capsids, as well as their self-assembly process, have been studied experimentally and theoretically, including modeling efforts by computer simulations on various scales. Atomistic models include specific details of local protein binding but are limited in system size and accessible time, while coarse grained (CG) models do get access to longer time and length scales but often lack the specific local interactions. Multi-scale models aim at bridging this gap by systematically connecting different levels of resolution. Here, a CG model for CCMV (Cowpea Chlorotic Mottle Virus), a virus with an icosahedral shell of 180 identical protein monomers, is developed, where parameters are derived from atomistic simulations of capsid protein dimers in aqueous solution. In particular, a new method is introduced to combine the MARTINI CG model with a supportive elastic network based on structural fluctuations of individual monomers. In the parametrization process, both network connectivity and strength are optimized. This elastic-network optimized CG model, which solely relies on atomistic data of small units (dimers), is able to correctly predict inter-protein conformational flexibility and properties of larger capsid fragments of 20 and more subunits. Furthermore, it is shown that this CG model reproduces experimental (Atomic Force Microscopy) indentation measurements of the entire viral capsid. Thus it is shown that one obvious goal for hierarchical modeling, namely predicting mechanical properties of larger protein complexes from models that are carefully parametrized on elastic properties of smaller units, is achievable.