Andrzej J. Rzepiela

Martini Coarse-Grained Force Field: Extension to Carbohydrates.

We present an extension of the Martini coarse-grained force field to carbohydrates. The parametrization follows the same philosophy as was used previously for lipids and proteins, focusing on the reproduction of partitioning free energies of small compounds between polar and nonpolar phases. The carbohydrate building blocks considered are the monosaccharides glucose and fructose and the disaccharides sucrose, trehalose, maltose, cellobiose, nigerose, laminarabiose, kojibiose, and sophorose. Bonded parameters for these saccharides are optimized by comparison to conformations sampled with an atomistic force field, in particular with respect to the representation of the most populated rotameric state for the glycosidic bond. Application of the new coarse-grained carbohydrate model to the oligosaccharides amylose and Curdlan shows a preservation of the main structural properties with 3 orders of magnitude more efficient sampling than the atomistic counterpart. Finally, we investigate the cryo- and anhydro-protective effect of glucose and trehalose on a lipid bilayer and find a strong decrease of the melting temperature, in good agreement with both experimental findings and atomistic simulation studies.

Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes

Cell membranes are comprised of multicomponent lipid and protein mixtures that exhibit a complex partitioning behavior. Regions of structural and compositional heterogeneity play a major role in the sorting and self-assembly of proteins, and their clustering into higher-order oligomers. Here, we use computer simulations and optical microscopy to study the sorting of transmembrane helices into the liquid-disordered domains of phase-separated model membranes, irrespective of peptide–lipid hydrophobic mismatch. Free energy calculations show that the enthalpic contribution due to the packing of the lipids drives the lateral sorting of the helices. Hydrophobic mismatch regulates the clustering into either small dynamic or large static aggregates. These results reveal important molecular driving forces for the lateral organization and self-assembly of transmembrane helices in heterogeneous model membranes, with implications for the formation of functional protein complexes in real cells.

Reconstruction of atomistic details from coarse-grained structures.

We present an algorithm to reconstruct atomistic structures from their corresponding coarse‐grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the application of the reconstruction procedure to obtain low‐energy atomistic structural ensembles from their CG counterparts. We reconstructed individual molecules in vacuo (NCQ tripeptide, dipalmitoylphosphatidylcholine, and cholesterol), bulk water, and a WALP transmembrane peptide embedded in a solvated lipid bilayer. The first examples serve to optimize the parameters for the reconstruction procedure, whereas the latter examples illustrate the applicability to condensed‐phase biomolecular systems.

Membrane poration by antimicrobial peptides combining atomistic and coarse-grained descriptions

Antimicrobial peptides (AMPs) comprise a large family of peptides that include small cationic peptides, such as magainins, which permeabilize lipid membranes. Previous atomistic level simulations of magainin-H2 peptides show that they act by forming toroidal transmembrane pores. However, due to the atomistic level of description, these simulations were necessarily limited to small system sizes and sub-microsecond time scales. Here, we study the long-time relaxation properties of these pores by evolving the systems using a coarse-grain (CG) description. The disordered nature and the topology of the atomistic pores are maintained at the CG level. The peptides sample different orientations but at any given time, only a few peptides insert into the pore. Key states observed at the CG level are subsequently back-transformed to the atomistic level using a resolution-transformation protocol. The configurations sampled at the CG level are stable in the atomistic simulation. The effect of helicity on pore stability is investigated at the CG level and we find that partial helicity is required to form stable pores. We also show that the current CG scheme can be used to study spontaneous poration by magainin-H2 peptides. Overall, our simulations provide a multi-scale view of a fundamental biophysical membrane process involving a complex interplay between peptides and lipids.

Efficient Algorithms for Langevin and DPD Dynamics

In this article, we present several algorithms for stochastic dynamics, including Langevin dynamics and different variants of Dissipative Particle Dynamics (DPD), applicable to systems with or without constraints. The algorithms are based on the impulsive application of friction and noise, thus avoiding the computational complexity of algorithms that apply continuous friction and noise. Simulation results on thermostat strength and diffusion properties for ideal gas, coarse-grained (MARTINI) water, and constrained atomic (SPC/E) water systems are discussed. We show that the measured thermal relaxation rates agree well with theoretical predictions. The influence of various parameters on the diffusion coefficient is discussed.

Mixing Martinis: Hybrid Atomistic/Coarse-Grained Models for Protein Molecular Dynamics

In recent years, the development and deployment of coarse grained models for simulations of proteins has taken an enormous flight. The main reason for this is that such models provide significant alleviation of the time scale limits that otherwise restrict the use of molecular simulations for biological processes. Coarse graining allows assessment of processes that occur on the scale of microseconds and micrometers, rather than nanoseconds and nanometers, albeit with the obvious consequence that detail is lost. This loss of detail has proven acceptable in many cases, but poses problems for the assessment of mechanical features of proteins, especially where local dynamics is intimately linked with overall conformational changes.To bring back the detail, yet only where it is needed, we have developed an integrative approach, coupling a Martini Coarse Grained model to an atomistic description of part of the system. This method involves a novel treatment of the interaction of the all-atom parts with the surrounding coarse grained particles, using virtual sites, rather than specific cross interactions. The potential applications of the method are manifold and include high-throughput protein-ligand binding studies, adsorption and protein folding.