Bottom-up coarse-graining of a simple graphene model: the blob picture

The coarse-graining of a simple all-atom 2D microscopic model of graphene, in terms of "blobs" described by center of mass variables, is presented. The equations of motion of the coarse-grained variables take the form of dissipative particle dynamics (DPD). The coarse-grained conservative forces and the friction of the DPD model are obtained via a bottom-up procedure from molecular dynamics (MD) simulations. The separation of timescales for blobs of 24 and 96 carbon atoms is sufficiently pronounced for the Markovian assumption, inherent to the DPD model, to provide satisfactory results. In particular, the MD velocity autocorrelation function of the blobs is well reproduced by the DPD model, provided that the effect of friction and noise is taken into account. However, DPD cross-correlations between neighbor blobs show appreciable discrepancies with respect to the MD results. Possible extensions to mend these discrepancies are briefly outlined.     

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.     

Coarse-grained rigid blob model for soft matter simulations

We have developed a coarse-grained multiscale molecular simulation method for soft matter systems that directly incorporates stereochemical information. We divide the material into disjoint groups of atoms or particles that move as separate rigid bodies; we call these groups "rigid blobs," hence the name coarse-grained rigid blob model. The method is enabled by the construction of transferable interblob potentials that approximate the net intermolecular interactions, as obtained from ab initio electronic structure calculations, other all-atom empirical potentials, experimental data, or any combination of the above. We utilize a multipolar expansion to obtain the interblob potential-energy functions. The series, which contains controllable approximations that allow us to estimate the errors, approaches the original intermolecular potential as the number of terms increases. Using a novel numerical algorithm, we can calculate the interblob potentials very efficiently in terms of a few interaction moment tensors. This reduces the labor well beyond what is required in standard molecular-dynamics calculations and allows large-scale simulations for temporal scales commensurate with characteristic times of nano- and mesoscale systems. A detailed derivation of the formulas is presented, followed by illustrative applications to several systems showing that the method can effectively capture realistic microscopic details and can easily extend to large-scale simulations.     

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.     

Efficient partitioning technique for computing the dynamics of intramolecular processes: radiationless transitions in pyrazine

An efficient QP partitioning algorithm to compute the eigenvalues, eigenvectors, and the dynamics of large molecular systems of a particular type is presented. Compared to straightforward diagonalization, the algorithm displays favorable scaling (proportional to N(T)2) as a function of N(T), the size of the Hamiltonian matrix. In addition, the algorithm is trivially parallelizable, necessitating no "cross-talk" between nodes, thus enjoying the full linear speedup of parallelization. Moreover, the method requires very modest storage space, even for extremely large matrices. The method has also been enhanced through the development of a coarse-grained approximation, enabling an increase of the basis set size to unprecedented levels (10(8)-10(10) in the current application). The QP algorithm is applied to the dynamics of electronic internal conversion in a 24 vibrational-mode model of pyrazine. A performance comparison with other dynamical methods is presented, along with results for the decay dynamics of pyrazine and a discussion of resonance line shapes.     

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.     

The central cell model: a mesoscopic hopping model for the study of the displacement autocorrelation function

On the mesoscale, the molecular motion in a microporous material can be represented as a sequence of hops between different pore locations and from one pore to the other. On the same scale, the memory effects in the motion of a tagged particle are embedded in the displacement autocorrelation function (DACF), the discrete counterpart of the velocity autocorrelation function (VACF). In this paper, a mesoscopic hopping model, based on a lattice-gas automata dynamics, is presented for the coarse-grained modeling of the DACF in a microporous material under conditions of thermodynamic equilibrium. In our model, that we will refer to as central cell model, the motion of one tagged particle is mimicked through probabilistic hops from one location to the other in a small lattice of cells where all the other particles are indistinguishable; the cells closest to the one containing the tagged particle are simulated explicitly in the canonical ensemble, whereas the border cells are treated as mean-field cells in the grand-canonical ensemble. In the present paper, numerical simulation of the central cell model are shown to provide the same results as a traditional lattice-gas simulation. Along with this a mean-field theory of self-diffusion which incorporates time correlations is discussed.     

Coarse-grained molecular dynamics modeling of associating fluids: thermodynamics, liquid structure, and dynamics in the limit of zero association strength

A continuous coarse-grained potential model for associating fluids, consisting of an off-center specific site bonded with a harmonic potential to a center particle, has been developed and used in canonical molecular dynamics simulations. The thermodynamic, structural, and dynamic properties of the limiting nonassociating reference coarse-grained fluid are investigated as functions of the mass distribution and bond strength between center and site particles. It is theoretically shown and confirmed by simulation that in this limit variations in these potential parameters do not alter the equation of state of the reference coarse-grained fluid but that they have profound influences on both the translational and the rotational dynamics. From the simulation results we arrive at some guidelines that should be kept in mind in the selection of appropriate values for the model parameters. This work provides the precursory knowledge for the study of coarse-grained associating fluids using the conventional molecular dynamics method.     

YUP: A Molecular Simulation Program for Coarse-Grained and Multi-Scaled Models

Coarse-grained models can be very different from all-atom models and are highly varied. Each class of model is assembled very differently and some models need customized versions of the standard molecular mechanics methods. The most flexible way to meet these diverse needs is to provide access to internal data structures and a programming language to manipulate these structures. We have created YUP, a general-purpose program for coarse-grained and multi-scaled models. YUP extends the Python programming language by adding new data types. We have then used the extended language to implement three classes of coarse-grained models. The coarse-grained RNA model type is an unusual non-linear polymer and the assembly was easily handled with a simple program. The molecular dynamics algorithm had to be extended for a coarse-grained DNA model so that it could detect a failure that is invisible to a standard implementation. A third model type took advantage of access to the force field to simulate the packing of DNA in viral capsid. We find that objects are easy to modify, extend and redeploy. Thus, new classes of coarse-grained models can be implemented easily.     

Coarse-grained molecular dynamics simulations of nanopatterning with multivalent inks

A coarse-grained molecular dynamics (MD) model is developed to study the multivalent, or multisite, binding of small functionalized dendrimer molecules to beta-cyclodextrin-terminated self-assembled monolayers, the so-called "molecular printboards" used to print "ink" molecules on surfaces with a high degree of positional control and specificity. Some current and future bionanotechnology applications are in the creation of nanoparticle assemblies, directed protein assembly, platforms for biosensing, and cell:surface attachment. The coarse-grained model allows us to probe up to microsecond timescales and model ink diffusion, crucial for the application of the printboard in, for example, medical diagnostics. Recent all-atom MD simulations identified and quantified the molecular strain limiting the stability of nanopatterns created with small dendrimer inks, and explained the different patterns obtained experimentally with different dendrimer inks. In the present work, the all-atom simulations are "scaled up" to longer timescales via coarse graining, without incurring significant additional computational expense, and, crucially, without significant loss in atom-scale detail, the coarse-grained MD simulations yielding properties similar to those obtained from the all-atom simulations. The anchoring of the ink molecules to the monolayer is of multivalent nature and the degree of multivalency shows a sharp dependence on temperature, control of temperature thus providing a further operational "switch" for directed molecular assembly. The computational protocol developed can, in principle, be extended to model any multivalent assembly, for example, virus-cell complexation.     

Coarse-graining kinetic networks in photosynthetic energy transport

Photosynthetic systems utilize hundreds of chlorophylls to collect sunlight and transport the energy to the reaction center with remarkably high quantum efficiency, however, the large size of the system together with the complex interactions among the components make it extremely challenging to understand the dynamics of light harvesting in large photosynthetic systems. To shed light on this problem, we present a structure-based theoretical framework that can be used to calculate transition rate matrix describing energy transport in photosynthetic systems and network clustering methods that provide simplified coarse-grained model revealing key structures guiding the light harvesting process. We constructed an effective model for energy transport in a Photosystem II supercomplex and applied several network clustering methods to generate coarse-grained kinetic cluster models for the system. Furthermore, we evaluated the performances of the network clustering methods, and show that a spectral clustering method and a minimum cut approach produce accurate coarse-grained models for the PSII-sc system. The results indicate that finding bottlenecks of energy transport is a crucial factor for reduced representations of photosynthetic light harvesting, and the overall work presented in this paper should provide a comprehensive theoretical framework to elucidate the dynamics of light harvesting in photosynthetic systems.     

A transferable coarse-grained model for hydrogen-bonding liquids

We present here a recent development of a generalized coarse-grained model for use in molecular simulations. In this model, interactions between coarse-grained particles consist of both van der Waals and explicit electrostatic components. As a result, the coarse-grained model offers the transferability that is lacked by most current effective-potential based approaches. The previous center-of-mass framework (P. A. Golubkov and P. Ren, J. Chem. Phys., 2006, 125, 64103) is generalized here to include arbitrary off-center interaction sites for both Gay-Berne and multipoles. The new model has been applied to molecular dynamic simulations of neat methanol liquid. By placing a single point multipole at the oxygen atom rather than at the center of mass of methanol, there is a significant improvement in the ability to capture hydrogen-bonding. The critical issue of transferability of the coarse-grained model is verified on methanol-water mixtures, using parameters derived from neat liquids without any modification. The mixture density and internal energy from coarse-grained molecular dynamics simulations show good agreement with experimental measurements, on a par with what has been obtained from more detailed atomic models. By mapping the dynamics trajectory from the coarse-grained simulation into the all-atom counterpart, we are able to investigate atomic-level structure and interaction. Atomic radial distribution functions of neat methanol, neat water and mixtures compare favorably to experimental measurements. Furthermore, hydrogen-bonded 6- and 7-molecule chains of water and methanol observed in the mixture are in agreement with previous atomic simulations.     

Mapping atomistic simulations to mesoscopic models: a systematic coarse-graining procedure for vinyl polymer chains

The paper introduces a systematic procedure to coarse-grain atomistic models of the largest family of synthetic polymers into a mesoscopic model that is able to keep detailed information about chain stereosequences. The mesoscopic model consists of sequences of superatoms centered on methylene carbons of two different types according to the kind of diad (m or r) they belong to. The corresponding force-field contains three different bonds, six angle and three nonbonded terms. Recently developed analytical potentials, based on sums of Gaussians for bond and angle terms of the mesoscale force field have been used. For the nonbonded part, numerical potentials optimized by pressure-corrected iterative Boltzmann inversion have been used. As test case we coarse-grained an atomistic all-atom model of atactic polystyrene. The proposed mesoscale model has been successfully tested against structural and dynamical properties for different chain lengths and opens the possibility of relaxing melts of high molecular weight vinyl polymers.     

Multigraining: an algorithm for simultaneous fine-grained and coarse-grained simulation of molecular systems

A method to combine fine-grained and coarse-grained simulations is presented. The coarse-grained particles are described as virtual particles defined by the underlying fine-grained particles are described as virtual particles defined by the underlying fine-grained particles. The contribution of the two grain levels to the interaction between particles is specified by a grain-level parameter lambda. Setting lambda = 0 results in a completely fine-grained simulation, whereas lambda = 1 yields a simulation governed by the coarse-grained potential energy surface with small contributions to keep the fine-grained covalently bound particles together. Simulations at different lambda values may be coupled using the replica-exchange molecular dynamics method to achieve enhanced sampling at the fine-grained level.     

Coarse-grained force field for the nucleosome from self-consistent multiscaling

A coarse-grained simulation model for the nucleosome is developed, using a methodology modified from previous work on the ribosome. Protein residues and DNA nucleotides are represented as beads, interacting through harmonic (for neighboring) or Morse (for nonbonded) potentials. Force-field parameters were estimated by Boltzmann inversion of the corresponding radial distribution functions obtained from a 5-ns all-atom molecular dynamics (MD) simulation, and were refined to produce agreement with the all-atom MD simulation. This self-consistent multiscale approach yields a coarse-grained model that is capable of reproducing equilibrium structural properties calculated from a 50-ns all-atom MD simulation. This coarse-grained model speeds up nucleosome simulations by a factor of 10(3) and is expected to be useful in examining biologically relevant dynamical nucleosome phenomena on the microsecond timescale and beyond.     

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.