Hybrid particle‐field molecular dynamics simulations: Parallelization and benchmarks

The parallel implementation of a recently developed hybrid scheme for molecular dynamics (MD) simulations (Milano and Kawakatsu, J Chem Phys 2009, 130, 214106) where self‐consistent field theory (SCF) and particle models are combined is described. Because of the peculiar formulation of the hybrid method, considering single particles interacting with density fields, the most computationally expensive part of the hybrid particle‐field MD simulation can be efficiently parallelized using a straightforward particle decomposition algorithm. Benchmarks of simulations, including comparisons of serial MD and MD‐SCF program profiles, serial MD‐SCF and parallel MD‐SCF program profiles, and parallel benchmarks compared with efficient MD program GROMACS 4.5.4 are tested and reported. The results of benchmarks indicate that the proposed parallelization scheme is very efficient and opens the way to molecular simulations of large scale systems with reasonable computational costs. © 2012 Wiley Periodicals, Inc.     

Prediction of viscoelastic properties with coarse‐grained molecular dynamics and experimental validation for a benchmark polyurea system

To explore the relationship between microscopic structure and viscoelastic properties of polyurea, a coarse‐grained (CG) model is developed by a structure matching method and validated against experiments conducted on a controlled, benchmark material. Using the Green‐Kubo method, the relaxation function is computed from the autocorrelation of the stress tensor, sampled over equilibrium MD simulations, and mapped to a real time scale established by matching self‐diffusion rates of atomistic and CG models. Master curves computed from the predicted stress relaxation function are then compared with dynamic mechanical analysis experiments mapped to a wide frequency range by time–temperature superposition, as well as measurements of ultrasonic shear wave propagation. Computational simulations from monodisperse and polydisperse configurations, representative of the benchmark polyurea, show excellent agreement with the experimental measurements over a multidecade range of loading frequency. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 797–810     

IBIsCO: a molecular dynamics simulation package for coarse-grained simulation

IBIsCO is a parallel molecular dynamics simulation package developed specially for coarse-grained simulations with numerical potentials derived by the iterative Boltzmann inversion (IBI) method (Reith et al., J Comput Chem 2003, 24, 1624). In addition to common features of molecular dynamics programs, the techniques of dissipative particle dynamics (Groot and Warren, J Chem Phys 1997, 107, 4423) and Lowe-Andersen dynamics (Lowe, Europhys Lett 1999, 47, 145) are implemented, which can be used both as thermostats and as sources of friction to compensate the loss of degrees of freedom by coarse-graining. The reverse nonequilibrium molecular dynamics simulation method (Müller-Plathe, Phys Rev E 1999, 59, 4894) for the calculation of viscosities is also implemented. Details of the algorithms, functionalities, implementation, user interfaces, and file formats are described. The code is parallelized using PE_MPI on PowerPC architecture. The execution time scales satisfactorily with the number of processors.     

A coarse‐grained MARTINI‐like force field for DNA unzipping in nanopores

In nanopore force spectroscopy (NFS) a charged polymer is threaded through a channel of molecular dimensions. When an electric field is applied across the insulating membrane, the ionic current through the nanopore reports on polymer translocation, unzipping, dissociation, and so forth. We present a new model that can be applied in molecular dynamics simulations of NFS. Although simplified, it does reproduce experimental trends and all‐atom simulations. The scaled conductivities in bulk solution are consistent with experimental results for NaCl for a wide range of electrolyte concentrations and temperatures. The dependence of the ionic current through a nanopore on the applied voltage is symmetric and, in the voltage range used in experiments (up to 2 V), linear and in good agreement with experimental data. The thermal stability and geometry of DNA is well represented. The model was applied to simulations of DNA hairpin unzipping in nanopores. The results are in good agreement with all‐atom simulations: the scaled translocation times and unzipping sequence are similar. © 2015 Wiley Periodicals, Inc.     

Simulation of DNA electrophoresis in systems of large number of solvent particles by coarse-grained hybrid molecular dynamics approach

Simulation of DNA electrophoresis facilitates the design of DNA separation devices. Various methods have been explored for simulating DNA electrophoresis and other processes using implicit and explicit solvent models. Explicit solvent models are highly desired but their applications may be limited by high computing cost in simulating large number of solvent particles. In this work, a coarse-grained hybrid molecular dynamics (CGH-MD) approach was introduced for simulating DNA electrophoresis in explicit solvent of large number of solvent particles. CGH-MD was tested in the simulation of a polymer solution and computation of nonuniform charge distribution in a cylindrical nanotube, which shows good agreement with observations and those of more rigorous computational methods at a significantly lower computing cost than other explicit-solvent methods. CGH-MD was further applied to the simulation of DNA electrophoresis in a polymer solution and in a well-studied nanofluidic device. Simulation results are consistent with observations and reported simulation results, suggesting that CGH-MD is potentially useful for studying electrophoresis of macromolecules and assemblies in nanofluidic, microfluidic, and microstructure array systems that involve extremely large number of solvent particles, nonuniformly distributed electrostatic interactions, bound and sequestered water molecules.     

Coarse‐grained molecular dynamics simulations of protein–ligand binding

Coarse‐grained molecular dynamics (CGMD) simulations with the MARTINI force field were performed to reproduce the protein–ligand binding processes. We chose two protein–ligand systems, the levansucrase–sugar (glucose or sucrose), and LinB–1,2‐dichloroethane systems, as target systems that differ in terms of the size and shape of the ligand‐binding pocket and the physicochemical properties of the pocket and the ligand. Spatial distributions of the Coarse‐grained (CG) ligand molecules revealed potential ligand‐binding sites on the protein surfaces other than the real ligand‐binding sites. The ligands bound most strongly to the real ligand‐binding sites. The binding and unbinding rate constants obtained from the CGMD simulation of the levansucrase–sucrose system were approximately 10 times greater than the experimental values; this is mainly due to faster diffusion of the CG ligand in the CG water model. We could obtain dissociation constants close to the experimental values for both systems. Analysis of the ligand fluxes demonstrated that the CG ligand molecules entered the ligand‐binding pockets through specific pathways. The ligands tended to move through grooves on the protein surface. Thus, the CGMD simulations produced reasonable results for the two different systems overall and are useful for studying the protein–ligand binding processes. © 2014 Wiley Periodicals, Inc.     

Fast and accurate computation schemes for evaluating vibrational entropy of proteins

Standard normal mode analysis (NMA) method is able to calculate vibrational entropy of proteins, but it is computationally intensive, especially for large proteins. To evaluate vibrational entropy efficiently and accurately, we, here, propose computation schemes based on coarse-grained NMA methods. This can be achieved by rescaling coarse-grained results with a specific factor that is derived on the basis of the linear correlation of protein vibrational entropy between standard NMA and coarse-grained NMA. Our coarse-grained NMA computation schemes can repeat correctly and efficiently the results of standard NMA for large proteins.     

Sodium chloride's effect on self‐assembly of diphenylalanine bilayer

Understanding self‐assembling peptides becomes essential in nanotechnology, thereby providing a bottom‐up method for fabrication of nanostructures. Diphenylalanine constitutes an outstanding building block that can be assembled into various nanostructures, including two‐dimensional bilayers or nanotubes, exhibiting superb mechanical properties. It is known that the effect of the ions is critical in conformational and chemical interactions of bilayers or membranes. In this study, we analyzed the effect of sodium chloride on diphenylalanine bilayer using coarse‐grained molecular dynamics simulations, and calculated the bending Young's modulus and the torsional modulus by applying normal modal analysis using an elastic network model. The results showed that sodium chloride dramatically increases the assembling efficiency and stability, thereby promising to allow the precise design and control of the fabrication process and properties of bio‐inspired materials. © 2016 Wiley Periodicals, Inc.     

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. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Comparing normal modes across different models and scales: Hessian reduction versus coarse‐graining

Dimension reduction is often necessary when attempting to reach longer length and time scales in molecular simulations. It is realized by constraining degrees of freedom or by coarse‐graining the system. When evaluating the accuracy of a dimensional reduction, there is a practical challenge: the models yield vectors with different lengths, making a comparison by calculating their dot product impossible. This article investigates mapping procedures for normal mode analysis. We first review a horizontal mapping procedure for the reduced Hessian techniques, which projects out degrees of freedom. We then design a vertical mapping procedure for the “implosion” of the all‐atom (AA) Hessian to a coarse‐grained scale that is based upon vibrational subsystem analysis. This latter method derives both effective force constants and an effective kinetic tensor. Next, a series of metrics is presented for comparison across different scales, where special attention is given to proper mass‐weighting. The dimension‐dependent metrics, which require prior mapping for proper evaluation, are frequencies, overlap of normal mode vectors, probability similarity, Hessian similarity, collectivity of modes, and thermal fluctuations. The dimension‐independent metrics are shape derivatives, elastic modulus, vibrational free energy differences, heat capacity, and projection on a predefined basis set. The power of these metrics to distinguish between reasonable and unreasonable models is tested on a toy alpha helix system and a globular protein; both are represented at several scales: the AA scale, a Gō‐like model, a canonical elastic network model, and a network model with intentionally unphysical force constants. Published 2012 Wiley Periodicals, Inc.     

A new algorithm for construction of coarse‐grained sites of large biomolecules

The development of coarse‐grained (CG) models for large biomolecules remains a challenge in multiscale simulations, including a rigorous definition of CG representations for them. In this work, we proposed a new stepwise optimization imposed with the boundary‐constraint (SOBC) algorithm to construct the CG sites of large biomolecules, based on the s cheme of essential dynamics CG. By means of SOBC, we can rigorously derive the CG representations of biomolecules with less computational cost. The SOBC is particularly efficient for the CG definition of large systems with thousands of residues. The resulted CG sites can be parameterized as a CG model using the normal mode analysis based fluctuation matching method. Through normal mode analysis, the obtained modes of CG model can accurately reflect the functionally related slow motions of biomolecules. The SOBC algorithm can be used for the construction of CG sites of large biomolecules such as F‐actin and for the study of mechanical properties of biomaterials. © 2015 Wiley Periodicals, Inc.     

Understanding Creep Behavior of Semicrystalline Polymer via Coarse‐Grained Modeling

Understanding the deformational and failure behaviors of thermoplastic semicrystalline polymers is crucial due to the practical usages in various engineering applications. Taking isotactic polypropylene (iPP) as a semicrystalline polymer model system, atomistically informed coarse‐grained (CG) molecular dynamics (MD) simulations are employed to investigate the creep behavior of iPP. The simulations reveal that there exists a threshold stress of about 20.0 MPa, above which the maximum strain of iPP within the simulation time span increases dramatically. From the strain‐time analysis, it is observed that the iPP exhibits an initial elastic deformation stage and a subsequent plastic stage at lower stress levels, while a three‐stage creep behavior including a third fracture stage is observed at higher stress levels. Specifically, at lower stress levels, the bonded energy increases continuously as the chains stretch steadily, while the nonbonded energy shows an initial increase followed by a steady decrease due to the interchain sliding. At higher stress levels, both bonded and nonbonded energies change dramatically at the third stage, resulting from accelerated chain stretching, unfolding, sliding, and breaking. This study provides physical insight into the creep behavior of iPP at a fundamental molecular level and highlights the important role of microstructural evolution of chains in the deformation of semicrystalline polymer materials. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1779–1791     

Introduction of steered molecular dynamics into UNRES coarse‐grained simulations package

In this article, an implementation of steered molecular dynamics (SMD) in coarse‐grain UNited RESidue (UNRES) simulations package is presented. Two variants of SMD have been implemented: with a constant force and a constant velocity. The huge advantage of SMD implementation in the UNRES force field is that it allows to pull with the speed significantly lower than the accessible pulling speed in simulations with all‐atom representation of a system, with respect to a reasonable computational time. Therefore, obtaining pulling speed closer to those which appear in the atomic force spectroscopy is possible. The newly implemented method has been tested for behavior in a microcanonical run to verify the influence of introduction of artificial constrains on keeping total energy of the system. Moreover, as time dependent artificial force was introduced, the thermostat behavior was tested. The new method was also tested via unfolding of the Fn3 domain of human contactin 1 protein and the I27 titin domain. Obtained results were compared with Gø‐like force field, all‐atom force field, and experimental results. © 2017 Wiley Periodicals, Inc.     

Polarizable coarse‐grained models for molecular dynamics simulation of liquid cyclohexane

Force field parameters for polarizable coarse‐grained (CG) supra‐atomic models of liquid cyclohexane are proposed. Two different bead sizes were investigated, one representing two fine‐grained (FG) CH₂r united atoms of the cyclohexane ring, and one representing three FG CH₂r united atoms. Electronic polarizability is represented by a massless charge‐on‐spring particle connected to each CG bead. The model parameters were calibrated against the experimental density and heat of vaporization of liquid cyclohexane, and the free energy of cyclohexane hydration. Both models show good agreement with thermodynamic properties of cyclohexane, yet overestimate the self‐diffusion. The dielectric properties of the polarizable models agree very well with experiment. © 2015 Wiley Periodicals, Inc.     

Coarse-graining of protein structures for the normal mode studies

The coarse-grained structural model such as Gaussian network has played a vital role in the normal mode studies for understanding protein dynamics related to biological functions. However, for the large proteins, the Gaussian network model is computationally unfavorable for diagonalization of Hessian (stiffness) matrix for the normal mode studies. In this article, we provide the coarse-graining method, referred to as "dynamic model condensation," which enables the further coarse-graining of protein structures consisting of small number of residues. It is shown that the coarser-grained structures reconstructed by dynamic model condensation exhibit the dynamic characteristics, such as low-frequency normal modes, qualitatively comparable to original structures. This sheds light on that dynamic model condensation and may enable one to study the large protein dynamics for gaining insight into biological functions of proteins.     

A comparison between parallelization approaches in molecular dynamics simulations on GPUs

We test the relative performances of two different approaches to the computation of forces for molecular dynamics simulations on graphics processing units. A “vertex‐based” approach, where a computing thread is started per particle, is compared to an “edge‐based” approach, where a thread is started per each potentially non‐zero interaction. We find that the former is more efficient for systems with many simple interactions per particle while the latter is more efficient if the system has more complicated interactions or fewer of them. By comparing computation times on more and less recent graphics processing unit technology, we predict that, if the current trend of increasing the number of processing cores—as opposed to their computing power—remains, the “edge‐based” approach will gradually become the most efficient choice in an increasing number of cases. © 2014 Wiley Periodicals, Inc.     

Molecular simulation of structure formation and rheological properties of mixtures of telechelic and monofunctional associating polymer

We study the phase behavior, the characteristic times, and the rheological properties under the steady shear flow of the mixtures consisting of telechelic and monofunctional associating polymers by a coarse‐grained molecular dynamics simulation. The mixtures form the transient networks, the closely packed spherical micelles, and the wormlike micelles. We confirm the molecular origins of the several characteristic times of the mixtures. The dependencies of the characteristic times on the composition ratio between telechelic and monofunctional associating polymers show good agreement with reported experimental results. Under the steady shear flow, the mixtures show the shear thinning induced by the change of the spatial configuration of the micelles. The telechelic associating polymers especially play an important role in connecting the micelles at the shear thinning regime and enhance the steady shear viscosity. Furthermore, at the wormlike micellar region, the mixtures show the second shear thinning initiated by the transformation of the association conformation of the telechelic associating polymers.