Semi-bottom-up coarse graining of water based on microscopic simulations

The generalized dissipative particle dynamics (DPD) equation derived from the generalized Langevin equation under Markovian approximations is used to simulate coarse-grained (CG) water cells. The mean force and the friction coefficients in the radial and transverse directions needed for DPD equation are obtained directly from the all atomistic molecular dynamics (AAMD) simulations. But the dissipative friction forces are overestimated in the Markovian approximation, which results in wrong dynamic properties for the CG water in the DPD simulations. To account for the non-Markovian dynamics, a rescaling factor is introduced to the friction coefficients. The value of the factor is estimated by matching the diffusivity of water. With this semi-bottom-up mapping method, the radial distribution function, the diffusion constant, and the viscosity of the coarse-grained water system computed with DPD simulations are all in good agreement with AAMD results. It bridges the microscopic level and mesoscopic level with consistent length and time scales.     

On Coarse-Graining by the Inverse Monte Carlo Method: Dissipative Particle Dynamics Simulations Made to a Precise Tool in Soft Matter Modeling

We present a promising coarse-graining strategy for linking micro- and mesoscales of soft matter systems. The approach is based on effective pairwise interaction potentials obtained from detailed atomistic molecular dynamics (MD) simulations, which are then used in coarse-grained dissipative particle dynamics (DPD) simulations. Here, the effective potentials were obtained by applying the inverse Monte Carlo method [Lyubartsev and Laaksonen, Phys. Rev. E. 52, 3730 (1995)] on a chosen subset of degrees of freedom described in terms of radial distribution functions. In our first application of the method, the effective potentials were used in DPD simulations of aqueous NaCl solutions. With the same computational effort we were able to simulate systems of one order of magnitude larger than the MD simulations. The results from the MD and DPD simulations are in excellent agreement.     

A new perspective on the coarse-grained dynamics of fluids

A computational methodology is presented that is designed to model, at a coarse-grained level, the mesoscale dynamics of fluids and potentially other forms of soft matter. Within a molecular dynamics simulation, "ghost" particles of a specific size, corresponding to the fundamental length-scale of coarse-graining, are used as micro-probes designed to respond to local mesoscale fluid flows and stress gradients. A subsequent coarse-grained model is then developed that incorporates both the coarse-grained mesoscale dynamics and isothermal compressibility of the original microscopic system. The method is applied to water and methanol. A contrast with dissipative particle dynamics (DPD) is also presented.     

The role of entanglements on the stability of microphase separated diblock copolymers in shear flow

Various possible orientations of lamellar structures of diblock copolymers under shear are investigated with respect to their stability. A Brownian dynamics model is put forward in which each diblock is modeled as a dumbbell. The blobs in each dumbbell are held together by a finite extendable nonlinear elastic (FENE) potential and interact with all surrounding blobs by soft dissipative particle dynamics (DPD) potentials. In addition to this, the blobs have the possibility to entangle with each other. The corresponding interactions lead to large viscosities, which, however, quickly diminish with increasing shear rate. This fact turns out to be crucial for the stabilization of structures with the lamellae parallel to the velocity-vorticity plane. As a second result it is found that asymmetry in the entanglement interactions stimulates the actual reorientation into this state.     

Dissipative Particle Dynamics Simulations of Polymer Brushes: Comparison with Molecular Dynamics Simulations

Summary: The structure of polymer brushes is investigated by dissipative particle dynamics (DPD) simulations that include explicit solvent particles. With an appropriate choice of the DPD interaction parameters , we obtain good agreement with previous molecular dynamics (MD) results where the good solvent behavior has been modeled by an effective Lennard–Jones potential. The present results confirm that DPD simulation techniques can be applied for large length scale simulations of polymer brushes. A relation between the different length scales and is established.

Polymer brush at a solid–liquid interface.     

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 dynamics simulations of polymer melts. I. Building potential of mean force for polyethylene and cis-polybutadiene

We present the derivation of coarse-grained force fields for two types of polymers, polyethylene (PE), and cis-polybutadiene (cis-PB), using the concept of potential of mean force. Coarse-grained force fields were obtained from microscopic simulations for several coarse-graining levels, i.e., different number of monomers lambda per mesoscopic unit called "bead." These force fields are then used in dissipative particle dynamics (DPD) simulations to study structural and dynamical properties of polymer melts of PE and cis-PB. The radial distribution functions g(R), the end-to-end distance R0, the end-to-end vector relaxation time tau, and the chain center of mass self-diffusion D(CM), are computed for different chain lengths at different coarse-graining factor lambda. Scaling laws typical of the Rouse regime are obtained for both polymers for chain lengths ranging from 6 to 50 beads. It is found that the end-to-end distance R0 obtained from DPD simulations agree well with values obtained from both microscopic simulations and experiments. The dependence of the friction coefficient used in DPD simulations versus the coarse-graining level is discussed in view of the overall scaling of the dynamical properties.     

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.     

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.     

Polymer-clay nanocomposites: a multiscale molecular modeling approach

A hierarchical procedure bridging the gap between atomistic and mesoscopic simulation for polymer-clay nanocomposite (PCN) design is presented. The dissipative particle dynamics (DPD) is adopted as the mesoscopic simulation technique, and the interaction parameters of the mesoscopic model are estimated by mapping the corresponding energy values obtained from atomistic molecular dynamics (MD) simulations. The predicted structure of the nylon 6 PCN system considered is in excellent agreement with previous experimental and atomistic simulation results.     

Peptide chain dynamics in light and heavy water: zooming in on internal friction

Frictional effects due to the chain itself, rather than the solvent, may have a significant effect on protein dynamics. Experimentally, such "internal friction" has been investigated by studying folding or binding kinetics at varying solvent viscosity; however, the molecular origin of these effects is hard to pinpoint. We consider the kinetics of disordered glycine-serine and α-helix forming alanine peptides and a coarse-grained protein folding model in explicit-solvent molecular dynamics simulations. By varying the solvent mass over more than two orders of magnitude, we alter only the solvent viscosity and not the folding free energy. Folding dynamics at the near-vanishing solvent viscosities accessible by this approach suggests that solvent and internal friction effects are intrinsically entangled. This finding is rationalized by calculation of the polymer end-to-end distance dynamics from a Rouse model that includes internal friction. An analysis of the friction profile along different reaction coordinates, extracted from the simulation data, demonstrates that internal as well as solvent friction varies substantially along the folding pathways and furthermore suggests a connection between friction and the formation of hydrogen bonds upon folding.     

Statistics of time-limited ensembles of bent DNA conformations

The paper considers statistical properties of ensembles of chain conformations obtained by short-time Brownian dynamics (BD) of a coarse-grained DNA model in order to find out if the conditions necessary for accurate evaluation of the polymer elasticity are attainable in atom-level molecular dynamics (MD) simulations. To measure the bending persistence length (PL) with a 10% error using data accumulated in a single trajectory of a double helix of 15 base pairs, dynamics should be continued for a few microseconds. However, these estimates should be scaled down by about 2 orders of magnitude because the bending dynamics of short double helices in MD features much smaller relaxation times. As a result, good qualitative agreement with the worm-like chain (WLC) theory is reached in MD after tens of nanoseconds. The presently accessible durations of MD trajectories provide reasonably accurate evaluation of DNA elasticity and allow modeling of its mesoscopic properties. The surprisingly fast bending dynamics of short double helices in MD suggests that the microscopic mechanisms of DNA flexibility differ from a simple harmonic model.     

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.     

Interfacial friction between semiflexible polymers and crystalline surfaces

The results obtained from molecular dynamics simulations of the friction at an interface between polymer melts and weakly attractive crystalline surfaces are reported. We consider a coarse-grained bead-spring model of linear chains with adjustable intrinsic stiffness. The structure and relaxation dynamics of polymer chains near interfaces are quantified by the radius of gyration and decay of the time autocorrelation function of the first normal mode. We found that the friction coefficient at small slip velocities exhibits a distinct maximum which appears due to shear-induced alignment of semiflexible chain segments in contact with solid walls. At large slip velocities, the friction coefficient is independent of the chain stiffness. The data for the friction coefficient and shear viscosity are used to elucidate main trends in the nonlinear shear rate dependence of the slip length. The influence of chain stiffness on the relationship between the friction coefficient and the structure factor in the first fluid layer is discussed.     

Elastic Network Models are Robust to Variations in Formalism

Understanding the functions of biomolecules requires insight not only from structures, but from dynamics as well. Often, the most interesting processes occur on time scales too slow for exploration by conventional molecular dynamics (MD) simulations. For this reason, alternative computational methods such as elastic network models (ENMs) have become increasingly popular. These simple, coarse-grained models represent molecules as beads connected by harmonic springs; the system's motions are solved analytically by normal mode analysis. In the past few years, many different formalisms for performing ENM calculations have emerged, and several have been optimized using all-atom MD simulations. In contrast to other studies, we have compared the various formalisms in a systematic, quantitative way. In this study, we optimize many ENM functional forms using a uniform dataset containing only long (> 1 μs) all-atom MD simulations. Our results show that all models once optimized produce spring constants for immediate neighboring residues that are orders of magnitude stiffer than more distal contacts. In addition, the statistical significance of ENM performance varied with model resolution. We also show that fitting long trajectories does not improve ENM performance due to a problem inherent in all network models tested: they underestimate the relative importance of the most concerted motions. Finally, we characterize ENMs' resilience by tessellating the parameter space to show that broad ranges of parameters produce similar quality predictions. Taken together our data reveals that choice of spring function and parameters are not vital to performance of a network model and that simple parameters can by derived "by hand" when no data is available for fitting, thus illustrating the robustness of these models.     

Mesoscale model of polymer melt structure: self-consistent mapping of molecular correlations to coarse-grained potentials

Development and application of coarse-graining methods to condensed phases of macromolecules is an active area of research. Multiscale modeling of polymeric systems using coarse-graining methods presents unique challenges. Here we apply a coarse-graining method that self-consistently maps structural correlations from detailed molecular dynamics (MD) simulations of alkane oligomers onto coarse-grained potentials using a combination of MD and inverse Monte Carlo methods. Once derived, the coarse-grained potentials allow computationally efficient sampling of ensemble of conformations of significantly longer polyethylene chains. Conformational properties derived from coarse-grained simulations are in excellent agreement with experiments. The level of coarse graining provides a control over the balance of computational efficiency and retention of chemical identity of the underlying polymeric system. Challenges to extension and application of this and similar structure-based coarse-graining methods to model dynamics and phase behavior in polymeric systems are briefly discussed.     

Semiflexible amphiphilic polymers: cylindrical-shaped, collagenlike, and toroidal structures

A coarse-grained model is used to study the conformational properties of semiflexible polymers with amphiphilic monomer units containing both hydrophilic and hydrophobic interaction sites. The hydrophobically driven conformational transitions are studied using molecular dynamics simulations for the chains of varying stiffness, as characterized by intrinsic Kuhn segment lengths that vary over a decade. It is shown that the energy of hydrophobic attraction required for the realization of the coil-to-globule transition increases with increasing chain stiffness. For rather stiff backbone, the coil-to-globule transition corresponds to a first order phase transition. We find that depending on the chain stiffness, a variety of thermodynamically stable anisometric chain morphologies are possible in a solvent selectively poor for hydrophobic sites of amphiphilic monomer units. For flexible chains, the amphiphilic polymer forms a cylindrical globule having blob structure with nearly spherical blobs. With increasing stiffness, the number of blobs composing the globule decreases and the shape of blobs transforms into elongated cylinder. Further increase in stiffness leads to compaction of macromolecules into a collagenlike structure when the chain folds itself several times and different strands wind round each other. In this state, the collagenlike structures coexist with toroidal globules, both conformations having approximately equal energies.     

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.     

Temperature‐dependent changes in the hydrogen bonded hard segment network and microphase morphology in a model polyurethane: Experimental and simulation studies

Hydrogen bonding between hard segments has a critical effect on the morphology and properties of polyurethanes. Influence of temperature on hydrogen bonded urethane network and melting behavior of a model semicrystalline segmented polyurethane was investigated by experiments and simulations. Polyurethane was synthesized by the stoichiometric reaction between p‐phenylene diisocyanate and poly(tetramethylene oxide) (PTMO) with a molecular weight of 1000 g/mol. Simulations were carried out using dissipative particle dynamics (DPD) and molecular dynamics (MD) approaches. Experimental melting behavior obtained by various techniques was compared with simulations. DPD simulations showed a room temperature microphase morphology consisting of a three‐dimensional hydrogen‐bonded urethane hard segment network in a continuous and amorphous PTMO matrix. The first‐order melting transitions of crystalline urethane hard segments observed during the continuous isobaric heating in DPD and MD simulations (340–360 K) were in reasonably good agreement with those observed experimentally, such as AFM (320–340 K), WAXS (330–360 K), and FTIR (320–350 K) measurements. Quantitative verification of the melting of urethane hard segments was demonstrated by sharp discontinuities in energy versus temperature plots obtained by MD simulations due to substantial decrease in the number of hydrogen bonds above 340 K. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 182–192