Development of accurate coarse-grained force fields for weakly polar groups by an indirect parameterization strategy

2021 ◽  
Vol 23 (11) ◽  
pp. 6763-6774
Author(s):  
Junjie Song ◽  
Mingwei Wan ◽  
Ying Yang ◽  
Lianghui Gao ◽  
Weihai Fang

An indirect coarse-grained force field parameterization strategy for weakly polar groups.

2018 ◽  
Vol 39 (28) ◽  
pp. 2360-2370 ◽  
Author(s):  
Adam K. Sieradzan ◽  
Artur Giełdoń ◽  
Yanping Yin ◽  
Yi He ◽  
Harold A. Scheraga ◽  
...  

2020 ◽  
Vol 22 (48) ◽  
pp. 28325-28338
Author(s):  
Debdas Dhabal ◽  
Tanmoy Patra

By means of molecular simulation, the osmotic coefficient of aqueous solution of BMIMCl ionic liquid is calculated to compare with the experimental data and use that to optimize two popular force fields available in the literature for bulk ILs.


2021 ◽  
Author(s):  
Joe G Greener ◽  
David T Jones

AbstractFinding optimal parameters for force fields used in molecular simulation is a challenging and time-consuming task, partly due to the difficulty of tuning multiple parameters at once. Automatic differentiation presents a general solution: run a simulation, obtain gradients of a loss function with respect to all the parameters, and use these to improve the force field. This approach takes advantage of the deep learning revolution whilst retaining the interpretability and efficiency of existing force fields. We demonstrate that this is possible by parameterising a simple coarse-grained force field for proteins, based on training simulations of up to 2,000 steps learning to keep the native structure stable. The learned potential matches chemical knowledge and PDB data, can fold and reproduce the dynamics of small proteins, and shows ability in protein design and model scoring applications. Problems in applying differentiable molecular simulation to all-atom models of proteins are discussed along with possible solutions. The learned potential, simulation scripts and training code are made available at https://github.com/psipred/cgdms.


Author(s):  
Mingwei Wan ◽  
Junjie Song ◽  
Ying Yang ◽  
Lianghui Gao ◽  
Wei-hai Fang

Coarse-grained (CG) molecular dynamics are powerful tools to access mesoscopic phenomenon and simultaneously record microscopic details, but currently the CG force fields (FFs) are still limited by low parameterization efficiency...


2019 ◽  
Author(s):  
David Wang ◽  
Piotr E. Marszalek

AbstractMolecular mechanics force fields have been shown to differ in their predictions of processes such as protein folding. To test how force field differences affect predicted protein behavior, we created a mechanically perturbed model of the beta-stranded I91 titin domain based on atomic force spectroscopy data and examined its refolding behavior using six different force fields. To examine the transferability of the force field discrepancies identified by this model, we compared the results to equilibrium simulations of the weakly helical peptide Ac-(AAQAA)3-NH2. The total simulation time was 80 µs. From these simulations we found significant differences in I91 perturbation refolding ability between force fields. Concurrently, Ac-(AAQAA)3-NH2 equilibration experiments indicated that although force fields have similar overall helical frequencies, they can differ in helical lifetimes. The combination of these results suggests that differences in force field parameterization may allow a more direct transition between the beta and alpha regions of the Ramachandran plot thereby affecting both beta-strand refolding ability and helical lifetimes. Furthermore, the combination of results suggests that using mechanically perturbed models can provide a controlled method to gain more insight into how force fields affect protein behavior.


PLoS ONE ◽  
2021 ◽  
Vol 16 (9) ◽  
pp. e0256990
Author(s):  
Joe G. Greener ◽  
David T. Jones

Finding optimal parameters for force fields used in molecular simulation is a challenging and time-consuming task, partly due to the difficulty of tuning multiple parameters at once. Automatic differentiation presents a general solution: run a simulation, obtain gradients of a loss function with respect to all the parameters, and use these to improve the force field. This approach takes advantage of the deep learning revolution whilst retaining the interpretability and efficiency of existing force fields. We demonstrate that this is possible by parameterising a simple coarse-grained force field for proteins, based on training simulations of up to 2,000 steps learning to keep the native structure stable. The learned potential matches chemical knowledge and PDB data, can fold and reproduce the dynamics of small proteins, and shows ability in protein design and model scoring applications. Problems in applying differentiable molecular simulation to all-atom models of proteins are discussed along with possible solutions and the variety of available loss functions. The learned potential, simulation scripts and training code are made available at https://github.com/psipred/cgdms.


2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  

<div>Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.</div><div>In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol<sup>-1</sup> for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div>


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>


2019 ◽  
Author(s):  
Kateryna Goloviznina ◽  
José N. Canongia Lopes ◽  
Margarida Costa Gomes ◽  
Agilio Padua

A general, transferable polarisable force field for molecular simulation of ionic liquids and their mixtures with molecular compounds is developed. This polarisable model is derived from the widely used CL\&P fixed-charge force field that describes most families of ionic liquids, in a form compatible with OPLS-AA, one of the major force fields for organic compounds. Models for ionic liquids with fixed, integer ionic charges lead to pathologically slow dynamics, a problem that is corrected when polarisation effects are included explicitly. In the model proposed here, Drude induced dipoles are used with parameters determined from atomic polarisabilities. The CL\&P force field is modified upon inclusion of the Drude dipoles, to avoid double-counting of polarisation effects. This modification is based on first-principles calculations of the dispersion and induction contributions to the van der Waals interactions, using symmetry-adapted perturbation theory (SAPT) for a set of dimers composed of positive, negative and neutral fragments representative of a wide variety of ionic liquids. The fragment approach provides transferability, allowing the representation of a multitude of cation and anion families, including different functional groups, without need to re-parametrise. Because SAPT calculations are expensive an alternative predictive scheme was devised, requiring only molecular properties with a clear physical meaning, namely dipole moments and atomic polarisabilities. The new polarisable force field, CL\&Pol, describes a broad set set of ionic liquids and their mixtures with molecular compounds, and is validated by comparisons with experimental data on density, ion diffusion coefficients and viscosity. The approaches proposed here can also be applied to the conversion of other fixed-charged force fields into polarisable versions.<br>


Sign in / Sign up

Export Citation Format

Share Document