Simulations of Phospholipids Using a Coarse Grain Model

The Martini coarse-grained force field

The Martini force field is a coarse-grained force field suited for molecular dynamics simulations of biomolecular systems. The force field has been parameterized in a systematic way, based on the reproduction of partitioning free energies between polar and apolar phases of a large number of chemical compounds. In this chapter the methodology underlying the force field is presented together with details of its parameterization and limitations. Then currently available topologies are described with a short overview of the key elements of their parameterization. These include the new polarizable Martini water model. A set of three selected ongoing studies using the Martini force field is presented. Finally the latest lines of development are discussed.

Accelerating membrane insertion of peripheral proteins with a novel membrane mimetic model

Characterizing atomic details of membrane binding of peripheral membrane proteins by molecular dynamics (MD) has been significantly hindered by the slow dynamics of membrane reorganization associated with the phenomena. To expedite lateral diffusion of lipid molecules without sacrificing the atomic details of such interactions, we have developed a novel membrane representation, to our knowledge, termed the highly mobile membrane-mimetic (HMMM) model to study binding and insertion of various molecular species into the membrane. The HMMM model takes advantage of an organic solvent layer to represent the hydrophobic core of the membrane and short-tailed phospholipids for the headgroup region. We demonstrate that using these components, bilayer structures are formed spontaneously and rapidly, regardless of the initial position and orientation of the lipids. In the HMMM membrane, lipid molecules exhibit one to two orders of magnitude enhancement in lateral diffusion. At the same time, the membrane atomic density profile of the headgroup region produced by the HMMM model is essentially identical to those obtained for full-membrane models, indicating the faithful representation of the membrane surface by the model. We demonstrate the efficiency of the model in capturing spontaneous binding and insertion of peripheral proteins by using the membrane anchor (γ-carboxyglutamic-acid-rich domain; GLA domain) of human coagulation factor VII as a test model. Achieving full insertion of the GLA domain consistently in 10 independent unbiased simulations within short simulation times clearly indicates the robustness of the HMMM model in capturing membrane association of peripheral proteins very efficiently and reproducibly. The HMMM model will provide significant improvements to the current all-atom models by accelerating lipid dynamics to examine protein-membrane interactions more efficiently.

A Simulation Study of the Self-Assembly of Coarse-Grained Skin Lipids

Computer simulations are an attractive means by which to probe the self-assembly and molecular level organization of lipids in biological membranes. In this work, we study a simple skin lipid system to demonstrate the ability of the coarse-grained models used for fatty acids, cholesterol, and water to self-assemble, thus validating the models for use in further studies of the complex lipid mixtures found in the outermost layer of the skin. Specifically, the ability of the models to predict the correct self-assembled structures from molecular dynamics simulations is compared against those seen experimentally and from all-atom simulations of preassembled bilayers. The nature of the molecular interactions and their roles in the self-assembly process is elucidated and heuristics for self-assembly established. Additionally, the coarse-grained models have been used to characterize the effect of varying cholesterol composition on bilayer properties and the mechanism of bilayer destabilization by short and long chain fatty acids in the presence of cholesterol.

Pore-spanning lipid membrane under indentation by a probe tip: a molecular dynamics simulation study

We study the indentation of a free-standing lipid membrane suspended over a nanopore on a hydrophobic substrate by means of molecular dynamics simulations. We find that in the course of indentation the membrane bends at the point of contact and the fringes of the membrane glide downward intermittently along the pore edges and stop gliding when the fringes reach the edge bottoms. The bending continues afterward, and the large strain eventually induces a phase transition in the membrane, transformed from a bilayered structure to an interdigitated structure. The membrane is finally ruptured when the indentation goes deep enough. Several local physical quantities in the pore regions are calculated, which include the tilt angle of lipid molecules, the nematic order, the included angle, and the distance between neighboring lipids. The variations of these quantities reveal many detailed, not-yet-specified local structural transitions of lipid molecules under indentation. The force-indentation curve is also studied and discussed. The results make a connection between the microscopic structure and the macroscopic properties and provide deep insight into the understanding of the stability of a lipid membrane spanning over nanopore.

A Solvent-Free Coarse Grain Model for Crystalline and Amorphous Cellulose Fibrils

Understanding biomass structure and dynamics on a range of time and length scales is important for the development of cellulosic biofuels. Here, to enable length and time scale extension, we develop a coarse grain (CG) model for molecular dynamics (MD) simulations of cellulose. For this purpose, we use distribution functions from fully atomistic MD simulations as target observables. A single bead per monomer level coarse graining is found to be sufficient to successfully reproduce structural features of crystalline cellulose. Without the use of constraints the CG crystalline fibril is found to remain stable over the maximum simulation length explored in this study (>1 μs). We also extend the CG representation to model fully amorphous cellulose fibrils. This is done by using an atomistic MD simulation of fully solvated individual cellulose chains as a target for developing the corresponding fully amorphous CG force field. Fibril structures with different degrees of crystallinity are obtained using force fields derived using a parameter coupling the crystalline and amorphous potentials. The method provides an accurate and constraint-free approach to derive CG models for cellulose with a wide range of crystallinity, suitable for incorporation into large-scale models of lignocellulosic biomass.

CELLmicrocosmos 2.2 MembraneEditor: a modular interactive shape-based software approach to solve heterogeneous membrane packing problems

New perspectives have been developed to understand the processes of modeling heterogeneous membranes. These are crucial steps prior to applying advanced techniques like molecular dynamic simulations of whole membrane systems. Lipid, protein, and membrane packing problems are addressed based on biochemical properties in combination with computational optimization techniques. The CELLmicrocosmos 2.2 MembraneEditor (CmME) is introduced as an appropriate framework to handle such problems by offering diverse algorithmic approaches. Its algorithm plug-in-interface enables modelers to generate problem-specific algorithms. Good solutions concerning runtime and lipid density are realized by focusing on the outer shapes of the PDB-based molecules. Application cases are presented like the publication-based modeling of inner and outer mitochondrial membrane-fragments, semiautomatic incorporation of proteins, and the assembly of rafts. Concerning geometrical aspects of the lipids, the achieved results are consistent with experimental observations related to lipid densities and distributions. Finally, two membranes simulated with GROMACS are analyzed and compared: the first is generated with conventional scripting techniques, the second with the CmME Distributor algorithm. The examples prove that CmME is a valuable and versatile tool for a broad set of applications in analysis and visualization of biomembranes.

Mimicking coarse-grained simulations without coarse-graining: enhanced sampling by damping short-range interactions

The damped-short-range-interaction (DSRI) method is proposed to mimic coarse-grained simulations by propagating an atomistic scale system on a smoothed potential energy surface. The DSRI method has the benefit of enhanced sampling provided by a typical coarse-grained simulation without the need to perform coarse-graining. Our method was used to simulate liquid water, alanine dipeptide folding, and the self-assembly of dimyristoylphosphatidylcholine lipid. In each case, our method appreciably accelerated the dynamics without significantly changing the free energy surface. Additional insights from DSRI simulations and the promise of coupling our DSRI method with Hamiltonian replica-exchange molecular dynamics are discussed.

Polarizable water model for the coarse-grained MARTINI force field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models.

Many biomolecular processes involve charged species moving between regions of high polarity, such as the water phase, and regions of lower polarity, such as the lipid membrane. Due to the change in electrostatic screening between these two environments, the strength of the interactions between the moving charge and the surrounding molecules also changes. This has important consequences for the way biological activity is controlled. To help understand the forces driving the movement of biomolecules, we developed a computational model which is capable of describing these processes at near-atomic detail. To do so efficiently, we use a coarse-grained description of the molecules, in which some of the atomistic detail is averaged out. To capture the inhomogeneous nature of the dielectric response, we re-introduce some detail in the water model; the new model effectively mimics the orientational polarizability of real water molecules, and screens electrostatic interactions realistically. This enables the study of a number of important biological processes that were hitherto considered challenging for coarse-grained models, such as the permeation of ions across a lipid membrane and the rupture of membranes due to an electrostatic field, at relatively low computational cost.

Zwitterionic lipid assemblies: molecular dynamics studies of monolayers, bilayers, and vesicles using a new coarse grain force field

A new coarse-grained (CG) intermolecular force field is presented for a series of zwitterionic lipids. The model is an extension of our previous work on nonionic surfactants and is designed to reproduce experimental surface/interfacial properties as well as distribution functions from all-atom molecular dynamics (MD) simulations. Using simple functional forms, the force field parameters are optimized for multiple lipid molecules, simultaneously. The resulting CG lipid bilayers have reasonable molecular areas, chain order parameters, and elastic properties. The computed surface pressure vs area (pi-A) curve for a dipalmitoyl phosphatidylcholine (DPPC) monolayer demonstrates a significant improvement over the previous CG models. The DPPC monolayer has a longer persistence length than a polyethyleneglycol (PEG) lipid monolayer, exhibiting a long-lived curved monolayer surface under negative tension. The bud ejected from an oversaturated DPPC monolayer has a large bicelle-like structure, which is different from the micellar bud formed from an oversaturated PEG lipid monolayer. We have successfully observed vesicle formation during CG-MD simulations, starting from an aggregate of dimyristoyl phosphatidylcholine (DMPC) molecules. Depending on the aggregate size, the lipid assembly spontaneously transforms into a closed vesicle or a bicelle. None of the various intermediate structures between these extremes seem to be stable. An attempt to observe fusion of two vesicles through the application of an external adhesion force was not successful. The present CG force field also supports stable multilamellar DMPC vesicles.

Coarse-grained molecular dynamics of tetrameric transmembrane peptide bundles within a lipid bilayer

The conformations of model transmembrane peptides are studied to understand the structural and dynamical aspects of tetrameric bundles using a series of coarse grain (CG) molecular dynamics (MD) simulations since membrane proteins play a crucial role in cell function. In this work, two different amphipathic models have been constructed using similar hydrophobic/hydrophilic characteristics with two structurally distinct morphologies to evaluate the effect of roughness and hydrophilic topology on the structure of tetrameric bundles, one class that forms an ion-channel and one class that does not. Free energy calculations of typical amphipathic peptide topologies show that using a relatively smooth surface morphology allows for a stable conformation of the tetramer bundle in a diamond formation. However, the model with side chains attached to the core in order to roughen the surface has a stable square tetramer bundle which is consistent with experimental data and all-atom (AA) MD simulations. Comparisons of the CG simulations with AA MD simulations are in reasonable agreement with the formation of tetrameric homo-oligomers, partitioning within the lipid bilayer and tilt angle with respect to the bilayer normal. We concluded that a square or diamond shape tetrameric homo-oligomers could be stabilized by rational design of the peptide morphology and topology of the surface, thus allowing us to tune the permeability of the bundle or channel.

Coarse-grained molecular dynamics simulation of the aggregation properties of multiheaded cationic surfactants in water

The aggregation property of multiheaded surfactants has been investigated by constant pressure molecular dynamics (MD) simulation in aqueous medium. The model multiheaded surfactants contain more than one headgroup (x = 2, 3, and 4) for a single tail group. This increases the hydrophilic charge progressively over the hydrophobic tail which has dramatic consequences in the aggregation behavior. In particular, we have looked at the change in the aggregation property such as critical micellar concentration (cmc), aggregation number, and size of the micelles for the multiheaded surfactants in water. We find with increasing number of headgroups of the multiheaded surfactants that the cmc values increase and the aggregation numbers as well as the size of the micelles decrease. These trends are in agreement with the experimental findings as reported earlier with x = 1, 2, and 3. We also predict the aggregation properties of multiheaded surfactant with four headgroups (x = 4) for which no experimental studies exist yet.

Simulations of membrane tubulation by lattices of amphiphysin N-BAR domains

Membrane compartments of manifold shapes are found in cells, often sculpted by cellular proteins. In particular, proteins of the BAR domain superfamily participate in membrane-sculpting processes in vivo and reshape also in vitro low-curvature membrane liposomes into high-curvature tubes and vesicles. Here we show by means of computer simulations totaling over 1 millisecond, how lattices involving parallel rows of amphiphysin N-BAR domains sculpt flat membranes into tubes. A highly detailed, dynamic picture of the 100-microsecond formation of membrane tubes by lattices of N-BAR domains is obtained. Lattice types inducing a wide range of membrane curvatures, with radii approximately 15-100 nm, are explored. The results suggest that multiple lattice types are viable for efficient membrane bending.

A systematic procedure to build a relaxed dense-phase atomistic representation of a complex amorphous polymer using a coarse-grained modeling approach

A systematic procedure has been developed to construct a relaxed dense-phase atomistic structure of a complex amorphous polymer. The numerical procedure consists of (1) coarse graining the atomistic model of the polymer into a mesoscopic model based on an iterative algorithm for potential inversion from distribution functions of the atomistic model, (2) relaxation of the coarse grained chain using a molecular dynamics scheme, and (3) recovery of the atomistic structure by reverse mapping based on the superposition of atomistic counterparts on the corresponding coarse grained coordinates. These methods are demonstrated by their application to construct a relaxed, dense-phase model of poly(DTB succinate), which is an amorphous tyrosine-derived biodegradable polymer that is being developed for biomedical applications. Both static and dynamic properties from the coarse-grained and atomistic simulations are analyzed and compared. The coarse-grained model, which contains the essential features of the DTB succinate structure, successfully described both local and global structural properties of the atomistic chain. The effective speedup compared to the corresponding atomistic simulation is substantially above 10(2), thus enabling simulation times to reach well into the characteristic experimental regime. The computational approach for reversibly bridging between coarse-grained and atomistic models provides an efficient method to produce relaxed dense-phase all-atom molecular models of complex amorphous polymers that can subsequently be used to study and predict the atomistic-level behavior of the polymer under different environmental conditions in order to optimally design polymers for targeted applications.

Maturation of high-density lipoproteins

Human high-density lipoproteins (HDLs) are involved in the transport of cholesterol. The mechanism by which HDL assembles and functions is not well understood owing to a lack of structural information on circulating spherical HDL. Here, we report a series of molecular dynamics simulations that describe the maturation of discoidal HDL into spherical HDL upon incorporation of cholesterol ester as well as the resulting atomic level structure of a mature circulating spherical HDL particle. Sixty cholesterol ester molecules were added in a stepwise fashion to a discoidal HDL particle containing two apolipoproteins wrapped around a 160 dipalmitoylphosphatidylcholine lipid bilayer. The resulting matured particle, captured in a coarse-grained description, was then described in a consistent all-atom representation and analysed in chemical detail. The simulations show that maturation results from the formation of a highly dynamic hydrophobic core comprised of cholesterol ester surrounded by phospholipid and protein; the two apolipoprotein strands remain in a belt-like conformation as seen in the discoidal HDL particle, but with flexible N- and C-terminal helices and a central region stabilized by salt bridges. In the otherwise flexible lipoproteins, a less mobile central region provides an ideal location to bind lecithin cholesterol acyltransferase, the key enzyme that converts cholesterol to cholesterol ester during HDL maturation.

Stability and rupture of archaebacterial cell membrane: a model study

It is known that the thermoacidophilic archaebacterium Sulfolobus acidocaldarius can grow in hot springs at 65-80 degrees C and live in acidic environments (pH 2-3); however, the origin of its unusual thermal stability remains unclear. In this work, using a vesicle as a model, we study the thermal stability and rupture of archaebacterial cell membrane. We perform a simulation investigation of the structure-property relationship of monolayer membrane formed by bolaform lipids and compare it with that of bilayer membrane formed by monopolar lipids. The origin of the unusually thermal stability of archaebacterial cell and the mechanism for its rupture are presented in molecular details.

Visualization of complex processes in lipid systems using computer simulations and molecular graphics

Computer simulation has become an increasingly popular tool in the study of lipid membranes, complementing experimental techniques by providing information on structure and dynamics at high spatial and temporal resolution. Molecular visualization is the most powerful way to represent the results of molecular simulations, and can be used to illustrate complex transformations of lipid aggregates more easily and more effectively than written text. In this chapter, we review some basic aspects of simulation methodologies commonly employed in the study of lipid membranes and we describe a few examples of complex phenomena that have been recently investigated using molecular simulations. We then explain how molecular visualization provides added value to computational work in the field of biological membranes, and we conclude by listing a few molecular graphics packages widely used in scientific publications.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.

The multiscale coarse-graining method. I. A rigorous bridge between atomistic and coarse-grained models

Coarse-grained (CG) models provide a computationally efficient method for rapidly investigating the long time- and length-scale processes that play a critical role in many important biological and soft matter processes. Recently, Izvekov and Voth introduced a new multiscale coarse-graining (MS-CG) method [J. Phys. Chem. B 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] for determining the effective interactions between CG sites using information from simulations of atomically detailed models. The present work develops a formal statistical mechanical framework for the MS-CG method and demonstrates that the variational principle underlying the method may, in principle, be employed to determine the many-body potential of mean force (PMF) that governs the equilibrium distribution of positions of the CG sites for the MS-CG models. A CG model that employs such a PMF as a "potential energy function" will generate an equilibrium probability distribution of CG sites that is consistent with the atomically detailed model from which the PMF is derived. Consequently, the MS-CG method provides a formal multiscale bridge rigorously connecting the equilibrium ensembles generated with atomistic and CG models. The variational principle also suggests a class of practical algorithms for calculating approximations to this many-body PMF that are optimal. These algorithms use computer simulation data from the atomically detailed model. Finally, important generalizations of the MS-CG method are introduced for treating systems with rigid intramolecular constraints and for developing CG models whose equilibrium momentum distribution is consistent with that of an atomically detailed model.

Soft patchy nanoparticles from solution-phase self-assembly of binary diblock copolymers

The fabrication of highly ordered, defect-free nanostructures is a key challenge in nanotechnology. Bottom-up fabrication approaches require nanobuilding blocks of precisely defined size and shape. In this work we propose a simple approach to obtain one type of building block--soft patchy nanoparticles--suggested by a series of coarse grain molecular dynamics simulations. A binary mixture of two different diblock copolymers with a common hydrophobic block but sufficiently dissimilar hydrophilic blocks reliably self-assembles into a "patchy" spherical micelle in water, with phase separation of the two hydrophilic blocks on the surface of the micelle core. Subsequent crosslinking of the core to solidify the patchy sphere geometry should allow further hierarchical assembly. Altering the hydrophilic versus hydrophobic composition of each polymer yields a change of morphology from "patchy spheres" to "patchy cylinders". Furthermore, by controlling the interaction strength of the blocks with solvent, the patches can be selectively placed either on the outer surface or inside the core of the micelle. The number and size of the patches are found to be largely controlled by the composition of the binary copolymer mixture.

Comparison of protein force fields for molecular dynamics simulations

In the context of molecular dynamics simulations of proteins, the term "force field" refers to the combination of a mathematical formula and associated parameters that are used to describe the energy of the protein as a function of its atomic coordinates. In this review, we describe the functional forms and parameterization protocols of the widely used biomolecular force fields Amber, CHARMM, GROMOS, and OPLS-AA. We also summarize the ability of various readily available noncommercial molecular dynamics packages to perform simulations using these force fields, as well as to use modern methods for the generation of constant-temperature, constant-pressure ensembles and to treat long-range interactions. Finally, we finish with a discussion of the ability of these force fields to support the modeling of proteins in conjunction with nucleic acids, lipids, carbohydrates, and/or small molecules.

Computer modeling in biotechnology: a partner in development

Computational modeling can be a useful partner in biotechnology, in particular, in nanodevice engineering. Such modeling guides development through nanoscale views of biomolecules and devices not available through experimental imaging methods. We illustrate the role of computational modeling, mainly of molecular dynamics, through four case studies: development of silicon bionanodevices for single molecule electrical recording, development of carbon nano-tube-biomolecular systems as in vivo sensors, development of lipoprotein nanodiscs for assays of single membrane proteins, and engineering of oxygen tolerance into the enzyme hydrogenase for photosynthetic hydrogen gas production. The four case studies show how molecular dynamics approaches were adapted to the specific technical uses through (i) multi-scale extensions, (ii) fast quantum chemical force field evaluation, (iii) coarse graining, and (iv) novel sampling methods. The adapted molecular dynamics simulations provided key information on device behavior and revealed development opportunities, arguing that the "computational microscope" is an indispensable nanoengineering tool.

Large multiprotein structures modeling and simulation: the need for mesoscopic models

Recent observational techniques based upon confocal microscopy make it possible to observe cells at a scale that has never been probed before: the mesoscopic scale. In the eukaryotic cell nucleus, many objects demonstrating phenomena occurring at this scale, such as nuclear bodies, are current subjects of investigations. But from a modeling perspective, this scale has not been widely explored, and hence there is a lack of suitable models for such studies. By reviewing higher and lower scale modeling techniques, we analyze their relevance in the context of mesoscale phenomena. We emphasize important characteristics that should be included in a mesoscopic model: an explicit continuous three-dimensional space with discrete simplified molecules that still have the characteristics of steric volume exclusion and realistic distant interaction forces. Then we present 3DSPI, a model dedicated to studies of nuclear bodies based on a simple formalism inspired from molecular dynamics and coarse-grained models: particles interacting through a potential energy function and driven by an overdamped Langevin equation. Finally, we present the features expected to be included in the model, pointing out the difficulties that might arise.

Coarse-grained free energy functions for studying protein conformational changes: a double-well network model

In this work, a double-well network model (DWNM) is presented for generating a coarse-grained free energy function that can be used to study the transition between reference conformational states of a protein molecule. Compared to earlier work that uses a single, multidimensional double-well potential to connect two conformational states, the DWNM uses a set of interconnected double-well potentials for this purpose. The DWNM free energy function has multiple intermediate states and saddle points, and is hence a "rough" free energy landscape. In this implementation of the DWNM, the free energy function is reduced to an elastic-network model representation near the two reference states. The effects of free energy function roughness on the reaction pathways of protein conformational change is demonstrated by applying the DWNM to the conformational changes of two protein systems: the coil-to-helix transition of the DB-loop in G-actin and the open-to-closed transition of adenylate kinase. In both systems, the rough free energy function of the DWNM leads to the identification of distinct minimum free energy paths connecting two conformational states. These results indicate that while the elastic-network model captures the low-frequency vibrational motions of a protein, the roughness in the free energy function introduced by the DWNM can be used to characterize the transition mechanism between protein conformations.

Comparison of atomic-level and coarse-grained models for liquid hydrocarbons from molecular dynamics configurational entropy estimates

Molecular liquids can be modeled at different levels of spatial resolution. In atomic-level (AL) models, all (heavy) atoms can be explicitly simulated. In coarse-grained (CG) models, particles (beads) that represent groups of covalently bound atoms are used as elementary units. Ideally, a CG model should reproduce the thermodynamic and structural properties of the corresponding AL model after mapping to the lower-resolution scale. In the present work, two such models are investigated: (i) the classical GROMOS atomic-level model; (ii) a CG model recently proposed by Marrink et al., which maps approximately four non-hydrogen atoms to one bead [J. Phys. Chem. B 2004, 108, 750]. The study is restricted to n-alkanes whose aliphatic fragments are abundantly found in lipids of biological interest. Additionally, cis-9-octadecene is included, as a template chain of the lipid dioleoylphosphatidylcholine (DOPC). The two representations of molecules in the liquid phase are compared in terms of average molecular structures, extent of configurational space sampled, and single-molecule entropies. An approximate method is used to estimate the rotational contributions to the absolute configurational entropy. Good correspondence between the AL and CG representations is found. The loss in configurational entropy due to the reduction in degrees of freedom upon coarse-graining of the model is estimated.

Coarse-grained peptide modeling using a systematic multiscale approach

A systematic new approach to derive multiscale coarse-grained (MS-CG) models has been recently developed. The approach employs information from atomistically detailed simulations to derive CG forces and associated effective potentials. In this work, the MS-CG methodology is extended to study two peptides representing distinct structural motifs, alpha-helical polyalanine and the beta-hairpin V(5)PGV(5). These studies represent the first known application of this approach to peptide systems. Good agreement between the MS-CG and atomistic models is achieved for several structural properties including radial distribution functions, root mean-square deviation, and radius of gyration. The new MS-CG models are able to preserve the native states of these peptides within approximately 1 A backbone root mean-square deviation during CG simulations. The MS-CG approach, as with most coarse-grained models, has the potential to increase the length and timescales accessible to molecular simulations. However, it is also able to maintain a clear connection to the underlying atomistic-scale interactions.

Stability and dynamics of virus capsids described by coarse-grained modeling

We report a study of the structural dynamics of viral capsids, simulated on a microsecond timescale, by employing a coarse-graining molecular dynamics method. The method was calibrated against an all-atom simulation of one complete virus. Among the studied capsids, some collapsed rapidly, while others were found to be stable. Interlocking between coat proteins is found to be a key factor determining the stability of the capsids.

Coarse-grained simulation of amphiphilic self-assembly

The authors present a computer simulation study of amphiphilic self-assembly performed using a computationally efficient single-site model based on Gay-Berne [J. Chem. Phys. 74, 3316 (1981)] and Lennard-Jones particles. Molecular dynamics simulations of these systems show that free self-assembly of micellar, bilayer, and inverse micelle arrangements can be readily achieved for a single model parametrization. This self-assembly is predominantly driven by the anisotropy of the amphiphile-solvent interaction, amphiphile-amphiphile dispersive interactions being found to be of secondary importance. While amphiphile concentration is the main determinant of phase stability, molecular parameters such as head group size and interaction strength also have measurable affects on system properties.