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     

Extension of CAVS coarse‐grained model to phospholipid membranes: The importance of electrostatics

It is evident from experiment that electrostatic potential (or dipole potential) is positive inside PC or PE lipid bilayers in the absence of ions. MARTINI coarse‐grained (CG) model, which has been widely used in simulating physical properties of lipid bilayers, fails to reproduce the positive value for the dipole potential in the membrane interior. Although the total dipole potential can be correctly described by the BMW/MARTINI model, the contribution from the ester dipoles, playing a nontrivial role in the electrostatic potential across lipid membranes, is neglected by this hybrid approach. In the ELBA CG model, the role of the ester dipoles is considered, but it is overweighed because various atomistic models have consistently shown that water is actually the leading contributor of dipole potential. Here, we present a CG approach by combining the BMW‐like water model (namely CAVS model) with the ELBA‐like lipid model proposed in this work. Our CG model was designed not only to correctly reproduce the positive values for the dipole potential inside PC and PE lipid bilayers but also to properly balance the individual contributions from the ester dipoles and water, surmounting the limitations of current CG models in the calculations of dipole potential. © 2017 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.     

Exploring the parameter space of the coarse‐grained UNRES force field by random search: Selecting a transferable medium‐resolution force field

We explored the energy‐parameter space of our coarse‐grained UNRES force field for large‐scale ab initio simulations of protein folding, to obtain good initial approximations for hierarchical optimization of the force field with new virtual‐bond‐angle bending and side‐chain‐rotamer potentials which we recently introduced to replace the statistical potentials. 100 sets of energy‐term weights were generated randomly, and good sets were selected by carrying out replica‐exchange molecular dynamics simulations of two peptides with a minimal α‐helical and a minimal β‐hairpin fold, respectively: the tryptophan cage (PDB code: 1L2Y) and tryptophan zipper (PDB code: 1LE1). Eight sets of parameters produced native‐like structures of these two peptides. These eight sets were tested on two larger proteins: the engrailed homeodomain (PDB code: 1ENH) and FBP WW domain (PDB code: 1E0L); two sets were found to produce native‐like conformations of these proteins. These two sets were tested further on a larger set of nine proteins with α or α + β structure and found to locate native‐like structures of most of them. These results demonstrate that, in addition to finding reasonable initial starting points for optimization, an extensive search of parameter space is a powerful method to produce a transferable force field. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

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.     

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.     

Implicit solvent coarse‐grained model of polyamidoamine dendrimers: Role of generation and pH

Highly branched polymers such as polyamidoamine (PAMAM) dendrimers are promising macromolecules in the realm of nanobiotechnology due to their high surface coverage of tunable functional groups. Modeling efforts of PAMAM can provide structural and morphological properties, but the inclusion of solvents and the exponential growth of atoms with generations make atomistic simulations computationally expensive. We apply an implicit solvent coarse‐grained model, called the Dry Martini force field, to PAMAM dendrimers. The reduced number of particles and the absence of a solvent allow the capture of longer spatiotemporal scales. This study characterizes PAMAM dendrimers of generations one through seven in acidic, neutral, and basic pH environments. Comparison with existing literature, both experimental and theoretical, is done using measurements of the radius of gyration, moment of inertia, radial distributions, and scaling exponents. Additionally, ion coordination distributions are studied to provide insight into the effects of interior and exterior protonation on counter ions. This model serves as a starting point for future designs of larger functionalized dendrimers. © 2015 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     

Melting thermodynamics of oligonucleic acids conjugated with relatively solvophobic linear polymers: A coarse‐grained molecular simulation study

Oligonucleic acids (ONAs), such as DNA and RNA, are used in various biotechnology and nanotechnology applications due to their ability to form a double helix that is stable at low temperature and melts at high temperatures. The melting temperature (T_m) of ONA duplexes can be tuned by the ONA composition, sequence, length and concentration, solvent quality, and salt concentration and by conjugation to other macromolecules. In this article, we use coarse‐grained (CG) molecular simulations to study ONAs conjugated with linear homopolymers that are relatively more solvophobic than the ONA. We study charged and stiff 8‐mer ONAs (e.g., DNA) and neutral and flexible 8‐mer ONAs (e.g., peptide nucleic acids or PNA), and vary the composition (or G‐C content) and sequence of ONA, conjugated homopolymer lengths and solvent quality for the polymer. For neutral and flexible ONAs, as the solvent quality worsens for the polymer, the ONA melting temperature increases from that of unconjugated ONA. The melting curves broaden with polymer length and worsening solvent quality, especially for ONAs with higher G‐C content. For charged and stiff ONAs, as the solvent quality worsens, the ONA melting temperature decreases compared to unconjugated ONA while the width of the melting curve remains the same. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1196–1208     

An anisotropic coarse‐grained model based on Gay–Berne and electric multipole potentials and its application to simulate a DMPC bilayer in an implicit solvent model

In this work, we aim at optimizing the performance of the anisotropic GBEMP model, which adopts a framework by combining a Gay–Berne (GB) anisotropic potential with an electric multipole (EMP) potential, in simulating a DMPC lipid bilayer in an implicit solvent model. First, the Gay–Berne parameters were initially obtained by fitting to atomistic profiles of van der Waals interactions between homodimers of molecular fragments while EMP parameters was directly derived from the expansion of point multipoles at predefined EMP sites. Second, the GB and EMP parameters for DMPC molecule were carefully optimized to be comparable to AMBER atomistic model in the calculations of the dipole moments of DMPC monomers adopting different conformations as well as the nonbonded interactions between two DMPC molecules adopting different conformations and separated at various distances. Finally, the GB parameters for DMPC were slightly adjusted in simulating a 72 DMPC bilayer system so that our GBEMP model would be able to reproduce a few important structural properties, namely, thickness (), area per lipid ( ) and volume per lipid ( ). Meanwhile, the atomistic and experimental results for electron density profiles and order parameters were reproduced reasonably well by the GBEMP model, demonstrating the promising feature of GBEMP model in modeling lipid systems. Finally, we have shown that current GBEMP model is more efficient by a factor of about 25 than AMBER atomistic point charge model. © 2015 Wiley Periodicals, Inc.     

Combining a polarizable force-field and a coarse-grained polarizable solvent model. II. Accounting for hydrophobic effects

A revised and improved version of our efficient polarizable force-field/coarse grained solvent combined approach (Masella, Borgis, and Cuniasse, J. Comput. Chem. 2008, 29, 1707) is described. The polarizable pseudo-particle solvent model represents the macroscopic solvent polarization by induced dipoles placed on mobile pseudo-particles. In this study, we propose a new formulation of the energy term handling the nonelectrostatic interactions among the pseudo-particles. This term is now able to reproduce the energetic and structural response of liquid water due to the presence of a hydrophobic spherical cavity. Accordingly, the parameters of the energy term handling the nonpolar solute/solvent interactions have been refined to reproduce the free-solvation energy of small solutes, based on a standard thermodynamic integration scheme. The reliability of this new approach has been checked for the properties of solvated methane and of the solvated methane dimer, as well as by performing 10 × 20 ns molecular dynamics (MD) trajectories for three solvated proteins. A long-time stability of the protein structures along the trajectories is observed. Moreover, our method still provides a measure of the protein solvation thermodynamic at the same accuracy as standard Poisson-Boltzman continuum methods. These results show the relevance of our approach and its applicability to massively coupled MD schemes to accurately and intensively explore solvated macromolecule potential energy surfaces.     

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.     

Molecule‐specific determination of atomic polarizabilities with the polarizable atomic multipole model

Recently, many polarizable force fields have been devised to describe induction effects between molecules. In popular polarizable models based on induced dipole moments, atomic polarizabilities are the essential parameters and should be derived carefully. Here, we present a parameterization scheme for atomic polarizabilities using a minimization target function containing both molecular and atomic information. The main idea is to adopt reference data only from quantum chemical calculations, to perform atomic polarizability parameterizations even when relevant experimental data are scarce as in the case of electronically excited molecules. Specifically, our scheme assigns the atomic polarizabilities of any given molecule in such a way that its molecular polarizability tensor is well reproduced. We show that our scheme successfully works for various molecules in mimicking dipole responses not only in ground states but also in valence excited states. The electrostatic potential around a molecule with an externally perturbing nearby charge also exhibits a near‐quantitative agreement with the reference data from quantum chemical calculations. The limitation of the model with isotropic atoms is also discussed to examine the scope of its applicability. © 2012 Wiley Periodicals, Inc.     

A fast method for large‐scale De Novo peptide and miniprotein structure prediction

Although peptides have many biological and biomedical implications, an accurate method predicting their equilibrium structural ensembles from amino acid sequences and suitable for large‐scale experiments is still missing. We introduce a new approach—PEP‐FOLD—to the de novo prediction of peptides and miniproteins. It first predicts, in the terms of a Hidden Markov Model‐derived structural alphabet, a limited number of local conformations at each position of the structure. It then performs their assembly using a greedy procedure driven by a coarse‐grained energy score. On a benchmark of 52 peptides with 9–23 amino acids, PEP‐FOLD generates lowest‐energy conformations within 2.8 and 2.3 Å Cα root‐mean‐square deviation from the full nuclear magnetic resonance structures (NMR) and the NMR rigid cores, respectively, outperforming previous approaches. For 13 miniproteins with 27–49 amino acids, PEP‐FOLD reaches an accuracy of 3.6 and 4.6 Å Cα root‐mean‐square deviation for the most‐native and lowest‐energy conformations, using the nonflexible regions identified by NMR. PEP‐FOLD simulations are fast—a few minutes only—opening therefore, the door to in silico large‐scale rational design of new bioactive peptides and miniproteins. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A systematic approach to calibrate a transferable polarizable force field parameter set for primary alcohols

In this work, parameters are optimized for a charge‐on‐spring based polarizable force field for linear alcohols. We show that parameter transferability can be obtained using a systematic approach in which the effects of parameter changes on physico‐chemical properties calculated from simulation are predicted. Our previously described QM/MM calculations are used to attribute condensed‐phase polarizabilities, and starting from the non‐polarizable GROMOS 53A5/53A6 parameter set, van der Waals and Coulomb interaction parameters are optimized to reproduce pure‐liquid (thermodynamic, dielectric, and transport) properties, as well as hydration free energies. For a large set of models, which were obtained by combining small perturbations of 10 distinct parameters, values for pure‐liquid properties of the series methanol to butanol were close to experiment. From this large set of models, we selected 34 models without special repulsive van der Waals parameters to distinguish between hydrogen‐bonding and non‐hydrogen‐bonding atom pairs, to make the force field simple and transparent. © 2017 Wiley Periodicals, Inc.     

Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short‐range penetration correction up to quadrupoles

We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short‐range charge penetration correction modifying the charge‐charge, charge‐dipole and charge‐quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry‐Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono‐ and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER‐HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration‐corrected polarizable force fields highlighting the mandatory need of non‐spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short‐range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio‐ or bioinorganic systems in periodic boundary conditions. © 2016 Wiley Periodicals, Inc.