Speed up of dynamic observables in coarse-grained molecular-dynamics simulations of unentangled polymers

Coarse-grained models that preserve atomistic detail display faster dynamics than atomistic systems alone. We show that this " indirect speed up" is robust: coarse-grained dynamic observables computed with time scaled by a constant factor are in excellent agreement with their underlying atomistic counterparts. Borrowing from accelerated dynamics methods used in the field of rare events, we predict the scaling factor within 7%, based on reduced intermolecular attraction yielding faster neighbor cage escapes.     

Developing a coarse-grained force field for the diblock copolymer poly(styrene-b-butadiene) from atomistic simulation

We have developed a coarse-grained force field for the poly(styrene-b-butadiene) diblock copolymer. We describe the computational methods and discuss how they were applied to develop a coarse-grained force field for this diblock copolymer from the atomistic simulation. The new force field contains three different bonds, four angles, five dihedral angles, and three nonbonded terms. We successfully tested this coarse-grained model against the chain properties, including static and dynamic properties, derived from the atomistic simulation; the results suggest that the coarse-grained force field is an effective model.     

An automatic coarse-graining and fine-graining simulation method: application on polyethylene

Multiscale modeling of a polymeric system is a challenging task in polymer physics. Here we introduce a bottom-up and then top-down scheme for the simulation of polyethylene (PE). The coarse-grained numerical potential for PE is derived through an automatic updating program by mapping its radial distribution function (RDF) from the Lowe-Andersen temperature controlling (LA) simulation onto the one from detailed molecular dynamics (MD) simulation. This coarse-grained numerical potential can be applied in larger systems under the same thermodynamic conditions. We have tested the reliability of the derived potential in two ways. First, the blends of different linear low-density polyethylene (LLDPE) with high-density polyethylene (HDPE) have been simulated in LA with the coarse-grained numerical potentials and reasonable results are obtained. Moreover, Rouse scaling behavior is reproduced for monodispersed polymeric systems with different chain lengths. The atomistic details of the beads can be reintroduced into the coarse-grained HDPE and LLDPE/HDPE models, followed by a few MD runs to alleviate the local tension induced by this fine-graining procedure. The equilibrated large atomistic system can then be used for further studies.     

An extension of the quasicontinuum treatment of multiscale solid systems to nonzero temperature

Covering the solid lattice with a finite-element mesh produces a coarse-grained system of mesh nodes as pseudoatoms interacting through an effective potential energy that depends implicitly on the thermodynamic state. Use of the pseudoatomic Hamiltonian in a Monte Carlo simulation of the two-dimensional Lennard-Jones crystal yields equilibrium thermomechanical properties (e.g., isotropic stress) in excellent agreement with "exact" fully atomistic results.     

The multiscale coarse-graining method. VII. Free energy decomposition of coarse-grained effective potentials

The potential of mean force (PMF) with respect to coarse-grained (CG) coordinates is often calculated in order to study the molecular interactions in atomistic molecular dynamics (MD) simulations. The multiscale coarse-graining (MS-CG) approach enables the computation of the many-body PMF of an atomistic system in terms of the CG coordinates, which can be used to parameterize CG models based on all-atom configurations. We demonstrate here that the MS-CG method can also be used to analyze the CG interactions from atomistic MD trajectories via PMF calculations. In addition, MS-CG calculations at different temperatures are performed to decompose the PMF values into energetic and entropic contributions as a function of the CG coordinates, which provides more thermodynamic information regarding the atomistic system. Two numerical examples, liquid methanol and a dimyristoylphosphatidylcholine lipid bilayer, are presented. The results show that MS-CG can be used as an analysis tool, comparable to various free energy computation methods. The differences between the MS-CG approach and other PMF calculation methods, as well as the characteristics and advantages of MS-CG, are also discussed.     

Coarse-grained representation of beta-helical protein building blocks

A general strategy to develop coarse-grained models of beta-helical protein fragments is presented. The procedure has been applied to a building block formed by a two-turn repeat motif from E. coli galactoside acetyltransferase, which is able to provide a very stable self-assembled tubular nanoconstruct upon stacking of its replicas. For this purpose, first, we have developed a computational scheme to sample very efficiently the configurational space of the building block. This method, which is inspired by a strategy recently designed to study amorphous polymers and by an advanced Monte Carlo algorithm, provides a large ensemble of uncorrelated configurations at a very reasonable computational cost. The atomistic configurations provided by this method have been used to obtain a coarse-grained model that describes the amino acids with fewer particles than those required for full atomistic detail, i.e., two, three, or four depending on the chemical nature of the amino acid. Coarse-grained potentials have been developed considering the following types of interactions: (i) electrostatic and van der Waals interactions between residues i and i + n with n >/= 2; (ii) interactions between residues i and i + 1; and (c) intra-residue interactions. The reliability of the proposed model has been tested by comparing the atomistic and coarse-grained energies calculated for a large number of independent configurations of the beta-helical building block.     

Parameterization of a coarse-grained model of cholesterol with point-dipole electrostatics

We present a new coarse-grained (CG) model of cholesterol (CHOL) for the electrostatic-based ELBA force field. A distinguishing feature of our CHOL model is that the electrostatics is modeled by an explicit point dipole which interacts through an ideal vacuum permittivity. The CHOL model parameters were optimized in a systematic fashion, reproducing the electrostatic and nonpolar partitioning free energies of CHOL in lipid/water mixtures predicted by full-detailed atomistic molecular dynamics simulations. The CHOL model has been validated by comparison to structural, dynamic and thermodynamic properties with experimental and atomistic simulation reference data. The simulation of binary DPPC/cholesterol mixtures covering the relevant biological content of CHOL in mammalian membranes is shown to correctly predict the main lipid behavior as observed experimentally.     

A multiscale coarse-graining method for biomolecular systems

A new approach is presented for obtaining coarse-grained (CG) force fields from fully atomistic molecular dynamics (MD) trajectories. The method is demonstrated by applying it to derive a CG model for the dimyristoylphosphatidylcholine (DMPC) lipid bilayer. The coarse-graining of the interparticle force field is accomplished by an application of a force-matching procedure to the force data obtained from an explicit atomistic MD simulation of the biomolecular system of interest. Hence, the method is termed a "multiscale" CG (MS-CG) approach in which explicit atomistic-level forces are propagated upward in scale to the coarse-grained level. The CG sites in the lipid bilayer application were associated with the centers-of-mass of atomic groups because of the simplicity in the evaluation of the forces acting on them from the atomistic data. The resulting CG lipid bilayer model is shown to accurately reproduce the structural properties of the phospholipid bilayer.     

Dissipative particle dynamic simulation study of lipid membrane

The lipid membrane plays crucial roles in countless biologic processes, ranging from cell motility, endo- and exocytosis, and cell division to protein aggregation and trafficking. To gain a molecular insight in these biologic processes, the recently developed mesoscale simulation technique, dissipative particle dynamics (DPD) simulation, has become an invaluable tool. By providing a brief survey of existing atomistic and popular coarse-grained models used today in studying the dynamics (including vesicle formation and (protein-mediated) vesicle fusion) and phase behavior of lipid bilayers, this review illustrates how mesoscopic DPD models can be used to obtain a better understanding of these biologic processes currently inaccessible to atomistic and most coarse-grained models.     

Local stress and heat flux in atomistic systems involving three-body forces

Local densities of fundamental physical quantities, including stress and heat flux fields, are formulated for atomistic systems involving three-body forces. The obtained formulas are calculable within an atomistic simulation, in consistent with the conservation equations of thermodynamics of continuum, and can be applied to systems with general two- and three-body interaction forces. It is hoped that this work may correct some misuse of inappropriate formulas of stress and heat flux in the literature, may clarify the definition of site energy of many-body potentials, and may serve as an analytical link between an atomistic model and a continuum theory. Physical meanings of the obtained formulas, their relation with virial theorem and heat theorem, and the applicability are discussed.     

Mesoscale simulation of water

A coarse-grained model for water is developed which maps five water molecules onto one bead. The coarse-grained potential is derived by iteratively matching the radial distribution function of water in the coarse-grained and target models. It is shown that the coarse-grained model has the optimal balance between computational efficiency and accuracy in structural properties. The model has been used to calculate a number of static and dynamic properties, including the density, isobaric thermal expansivity, isothermal compressibility, surface tension, and diffusion coefficient of water. The effect of coarsening on these properties has been discussed. It is shown that while the present coarse-grained model well describes the structural properties of water, expectedly it has a much faster dynamics than the corresponding atomistic models.     

Monte Carlo simulations of polymer dynamics: Recent advances

A brief review is given of applications of Monte Carlo simulations to study the dynamical properties of coarse-grained models of polymer melts, emphasizing the crossover from the Rouse model toward reptation, and the glass transition. The extent to which Monte Carlo algorithms can mimic the actual chain dynamics is critically examined, and the need for the use of coarse-grained rather than fully atomistic models for such simulations is explained. It is shown that various lattice and continuum models yield qualitatively similar results, and the behavior agrees with the findings of corresponding molecular dynamics simulations and experiments, where available. It is argued that these simulations significantly enhance our understanding of the theoretical concepts on the dynamics of dense macromolecular systems. © 1997 John Wiley & Sons, Inc.