scholarly journals Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields

Molecules ◽  
2021 ◽  
Vol 26 (19) ◽  
pp. 6069
Author(s):  
Ilias Patmanidis ◽  
Riccardo Alessandri ◽  
Alex H. de de Vries ◽  
Siewert J. Marrink
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Force Fields ◽  
Supramolecular Assembly ◽  
Model Systems ◽  
Binding Modes ◽  
Free Energies ◽  
Care Needs ◽  
Strength Of Interaction ◽  
Dynamics Simulations

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly.

scholarly journals Binding thermodynamics of host-guest systems with SMIRNOFF99Frosst 1.0.5 from the Open Force Field Initiative

2019 ◽  
Author(s):  
David Slochower ◽  
Niel Henriksen ◽  
Lee-Ping Wang ◽  
John Chodera ◽  
David Mobley ◽  
...  
Keyword(s):  
Ligand Binding ◽  
Force Field ◽  
Small Molecule ◽  
Force Fields ◽  
Free Energy Calculations ◽  
Model Systems ◽  
Conformational Space ◽  
Free Energies ◽  
Binding Thermodynamics ◽  
Binding Free Energies

<div><div><div><p>Designing ligands that bind their target biomolecules with high affinity and specificity is a key step in small- molecule drug discovery, but accurately predicting protein-ligand binding free energies remains challenging. Key sources of errors in the calculations include inadequate sampling of conformational space, ambiguous protonation states, and errors in force fields. Noncovalent complexes between a host molecule with a binding cavity and a drug-like guest molecules have emerged as powerful model systems. As model systems, host-guest complexes reduce many of the errors in more complex protein-ligand binding systems, as their small size greatly facilitates conformational sampling, and one can choose systems that avoid ambiguities in protonation states. These features, combined with their ease of experimental characterization, make host-guest systems ideal model systems to test and ultimately optimize force fields in the context of binding thermodynamics calculations.</p><p><br></p><p>The Open Force Field Initiative aims to create a modern, open software infrastructure for automatically generating and assessing force fields using data sets. The first force field to arise out of this effort, named SMIRNOFF99Frosst, has approximately one tenth the number of parameters, in version 1.0.5, compared to typical general small molecule force fields, such as GAFF. Here, we evaluate the accuracy of this initial force field, using free energy calculations of 43 α and β-cyclodextrin host-guest pairs for which experimental thermodynamic data are available, and compare with matched calculations using two versions of GAFF. For all three force fields, we used TIP3P water and AM1-BCC charges. The calculations are performed using the attach-pull-release (APR) method as implemented in the open source package, pAPRika. For binding free energies, the root mean square error of the SMIRNOFF99Frosst calculations relative to experiment is 0.9 [0.7, 1.1] kcal/mol, while the corresponding results for GAFF 1.7 and GAFF 2.1 are 0.9 [0.7, 1.1] kcal/mol and 1.7 [1.5, 1.9] kcal/mol, respectively, with 95% confidence ranges in brackets. These results suggest that SMIRNOFF99Frosst performs competitively with existing small molecule force fields and is a parsimonious starting point for optimization.</p></div></div></div>

scholarly journals Non-equilibrium approach for binding free energies in cyclodextrins in SAMPL7: force fields and software

2020 ◽  
Author(s):  
Yuriy Khalak ◽  
Gary Tresadern ◽  
Bert L. de Groot ◽  
Vytautas Gapsys
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Force Fields ◽  
Simulation Software ◽  
Energy Estimates ◽  
Energy Calculation ◽  
Free Energies ◽  
Guest Molecules ◽  
Software Packages ◽  
Non Equilibrium

AbstractIn the current work we report on our participation in the SAMPL7 challenge calculating absolute free energies of the host–guest systems, where 2 guest molecules were probed against 9 hosts-cyclodextrin and its derivatives. Our submission was based on the non-equilibrium free energy calculation protocol utilizing an averaged consensus result from two force fields (GAFF and CGenFF). The submitted prediction achieved accuracy of $${1.38}\,\hbox {kcal}/\hbox {mol}$$ 1.38 kcal / mol in terms of the unsigned error averaged over the whole dataset. Subsequently, we further report on the underlying reasons for discrepancies between our calculations and another submission to the SAMPL7 challenge which employed a similar methodology, but disparate ligand and water force fields. As a result we have uncovered a number of issues in the dihedral parameter definition of the GAFF 2 force field. In addition, we identified particular cases in the molecular topologies where different software packages had a different interpretation of the same force field. This latter observation might be of particular relevance for systematic comparisons of molecular simulation software packages. The aforementioned factors have an influence on the final free energy estimates and need to be considered when performing alchemical calculations.

scholarly journals Binding thermodynamics of host-guest systems with SMIRNOFF99Frosst 1.0.5 from the Open Force Field Initiative

2019 ◽  
Author(s):  
David Slochower ◽  
Niel Henriksen ◽  
Lee-Ping Wang ◽  
John Chodera ◽  
David Mobley ◽  
...  
Keyword(s):  
Ligand Binding ◽  
Force Field ◽  
Small Molecule ◽  
Force Fields ◽  
Free Energy Calculations ◽  
Model Systems ◽  
Conformational Space ◽  
Free Energies ◽  
Binding Thermodynamics ◽  
Binding Free Energies

<div><div><div><p>Designing ligands that bind their target biomolecules with high affinity and specificity is a key step in small- molecule drug discovery, but accurately predicting protein-ligand binding free energies remains challenging. Key sources of errors in the calculations include inadequate sampling of conformational space, ambiguous protonation states, and errors in force fields. Noncovalent complexes between a host molecule with a binding cavity and a drug-like guest molecules have emerged as powerful model systems. As model systems, host-guest complexes reduce many of the errors in more complex protein-ligand binding systems, as their small size greatly facilitates conformational sampling, and one can choose systems that avoid ambiguities in protonation states. These features, combined with their ease of experimental characterization, make host-guest systems ideal model systems to test and ultimately optimize force fields in the context of binding thermodynamics calculations.</p><p><br></p><p>The Open Force Field Initiative aims to create a modern, open software infrastructure for automatically generating and assessing force fields using data sets. The first force field to arise out of this effort, named SMIRNOFF99Frosst, has approximately one tenth the number of parameters, in version 1.0.5, compared to typical general small molecule force fields, such as GAFF. Here, we evaluate the accuracy of this initial force field, using free energy calculations of 43 α and β-cyclodextrin host-guest pairs for which experimental thermodynamic data are available, and compare with matched calculations using two versions of GAFF. For all three force fields, we used TIP3P water and AM1-BCC charges. The calculations are performed using the attach-pull-release (APR) method as implemented in the open source package, pAPRika. For binding free energies, the root mean square error of the SMIRNOFF99Frosst calculations relative to experiment is 0.9 [0.7, 1.1] kcal/mol, while the corresponding results for GAFF 1.7 and GAFF 2.1 are 0.9 [0.7, 1.1] kcal/mol and 1.7 [1.5, 1.9] kcal/mol, respectively, with 95% confidence ranges in brackets. These results suggest that SMIRNOFF99Frosst performs competitively with existing small molecule force fields and is a parsimonious starting point for optimization.</p></div></div></div>

FOLDING FREE ENERGY LANDSCAPE OF THE DECAPEPTIDE CHIGNOLIN

2008 ◽  
Vol 22 (31) ◽  
pp. 3087-3098 ◽  
Author(s):  
XIANGHUA DOU ◽  
JIHUA WANG
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Force Fields ◽  
Energy Landscapes ◽  
Native Structure ◽  
Good Template ◽  
Free Energy Landscapes ◽  
Folded Structures ◽  
Dynamics Simulations ◽  
Folding Free Energy

Chignolin is an artificially designed ten-residue (GYDPETGTWG) folded peptide, which is the smallest protein and provides a good template for protein folding. In this work, we completed four explicit water molecular dynamics simulations of Chignolin folding using GROMOS and OPLS-AA force fields from extended initial states without any experiment informations. The four-folding free energy landscapes of the peptide has been drawn. The folded state of Chignolin has been successfully predicated based on the free energy landscapes. The four independent simulations gave similar results. (i) The four free energy landscapes have common characters. They are fairly smooth, barrierless, funnel-like and downhill without intermediate state, which consists with the experiment. (ii) The different extended initial structures converge at similar folded structures with the lowest free energy under GROMOS and OPLS-AA force fields. In the GROMOS force field, the backbone RMSD of the folded structures from the NMR native structure of Chignolin is only 0.114 nm, which is a stable structure in this force field. In the OPLS-AA force field, the similar results have been obtained. In addition, the smallest RMSD structure is in better agreement with the NMR native structure but unlikely stable in the force field.

Binding thermodynamics of host-guest systems with SMIRNOFF99Frosst 1.0.5 from the Open Force Field Initiative

2019 ◽  
Author(s):  
David Slochower ◽  
Niel Henriksen ◽  
Lee-Ping Wang ◽  
John Chodera ◽  
David Mobley ◽  
...  
Keyword(s):  
Ligand Binding ◽  
Force Field ◽  
Small Molecule ◽  
Force Fields ◽  
Free Energy Calculations ◽  
Model Systems ◽  
Conformational Space ◽  
Free Energies ◽  
Binding Thermodynamics ◽  
Binding Free Energies

<div><div><div><p>Designing ligands that bind their target biomolecules with high affinity and specificity is a key step in small- molecule drug discovery, but accurately predicting protein-ligand binding free energies remains challenging. Key sources of errors in the calculations include inadequate sampling of conformational space, ambiguous protonation states, and errors in force fields. Noncovalent complexes between a host molecule with a binding cavity and a drug-like guest molecules have emerged as powerful model systems. As model systems, host-guest complexes reduce many of the errors in more complex protein-ligand binding systems, as their small size greatly facilitates conformational sampling, and one can choose systems that avoid ambiguities in protonation states. These features, combined with their ease of experimental characterization, make host-guest systems ideal model systems to test and ultimately optimize force fields in the context of binding thermodynamics calculations.</p><p><br></p><p>The Open Force Field Initiative aims to create a modern, open software infrastructure for automatically generating and assessing force fields using data sets. The first force field to arise out of this effort, named SMIRNOFF99Frosst, has approximately one tenth the number of parameters, in version 1.0.5, compared to typical general small molecule force fields, such as GAFF. Here, we evaluate the accuracy of this initial force field, using free energy calculations of 43 α and β-cyclodextrin host-guest pairs for which experimental thermodynamic data are available, and compare with matched calculations using two versions of GAFF. For all three force fields, we used TIP3P water and AM1-BCC charges. The calculations are performed using the attach-pull-release (APR) method as implemented in the open source package, pAPRika. For binding free energies, the root mean square error of the SMIRNOFF99Frosst calculations relative to experiment is 0.9 [0.7, 1.1] kcal/mol, while the corresponding results for GAFF 1.7 and GAFF 2.1 are 0.9 [0.7, 1.1] kcal/mol and 1.7 [1.5, 1.9] kcal/mol, respectively, with 95% confidence ranges in brackets. These results suggest that SMIRNOFF99Frosst performs competitively with existing small molecule force fields and is a parsimonious starting point for optimization.</p></div></div></div>

Hydration Free Energies of Organic Molecules in the FreeSolv Database Calculated with Polarized Atom In Molecules Atomic Charges and the GAFF Force Field.

2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  
Keyword(s):  
Force Field ◽  
Functional Groups ◽  
Electrostatic Interactions ◽  
Organic Molecules ◽  
Force Fields ◽  
Chemical Elements ◽  
Atomic Charges ◽  
Free Energies ◽  
Molecular Systems ◽  
Solvent Model

<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>

scholarly journals Adsorption of Amino Acids on Graphene: Assessment of Current Force Fields

2019 ◽  
Author(s):  
Siva Dasetty ◽  
John K. Barrows ◽  
Sapna Sarupria
Keyword(s):  
Amino Acids ◽  
Amino Acid ◽  
Force Field ◽  
Force Fields ◽  
Principal Component ◽  
Free Energies ◽  
Adsorbed State ◽  
Free Energy Of Adsorption ◽  
Density Based Clustering ◽  
Structural Preferences

<div> <div> <div> <p>We compare the free energies of adsorption (∆Aads) and the structural preferences of amino acids obtained using the force fields — Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid–graphene interactions are favorable irrespective of the force field. While the magnitudes of ∆Aads differ between the force fields, the trends in the free energy of adsorption with amino acids are similar across the studied force fields. ∆Aads positively correlates with amino acid–graphene and negatively correlates with graphene–water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that while the thermodynamics of adsorption of proteins and peptides would be described consistently across different force fields, the structural preferences of peptides and proteins on graphene will be force field dependent. </p> </div> </div> </div>

scholarly journals GAFF/IPolQ-Mod+LJ-Fit: Optimized force field parameters for solvation free energy predictions

ADMET & DMPK ◽  
2020 ◽  
Author(s):  
Andreas Mecklenfeld ◽  
Gabriele Raabe
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Fluid Phase ◽  
Solvation Free Energy ◽  
Rational Drug Design ◽  
Phase Behaviour ◽  
Liquid Density ◽  
Free Energies ◽  
Implicit Representation ◽  
Amber Force Field

<p class="ADMETabstracttext">Rational drug design featuring explicit solubility considerations can greatly benefit from molecular dynamics simulations, as they allow for the prediction of the Gibbs free energy of solvation and thus relative solubilities. In our previous work (A. Mecklenfeld, G. Raabe. J. Chem. Theory Comput. <strong>13 </strong>no. 12 (2017) 6266–6274), we have compared solvation free energy results obtained with the General Amber Force Field (GAFF) and its default restrained electrostatic potential (RESP) partial charges to those obtained by modified implicitly polarized charges (IPolQ-Mod) for an implicit representation of impactful polarization effects. In this work, we have adapted Lennard-Jones parameters for GAFF atom types in combination with IPolQ-Mod to further improve the accuracies of solvation free energy and liquid density predictions. We thereby focus on prominent atom types in common drugs. For the refitting, 357 respectively 384 systems were considered for free energies and densities and validation was performed for 142 free energies and 100 densities of binary mixtures. By the in-depth comparison of simulation results for default GAFF, GAFF with IPolQ-Mod and our new set of parameters, which we label GAFF/IPolQ-Mod+LJ-Fit, we can clearly highlight the improvements of our new model for the description of both relative solubilities and fluid phase behaviour.</p>

scholarly journals W-RESP: Well-Restrained Electrostatic Potential Derived Charges. Revisiting the Charge Derivation Model

2020 ◽  
Author(s):  
Michal Janeček ◽  
Petra Kührová ◽  
Vojtěch Mlýnský ◽  
Michal Otyepka ◽  
Jiří Šponer ◽  
...  
Keyword(s):  
Force Field ◽  
Electrostatic Potential ◽  
Electrostatic Interactions ◽  
Force Fields ◽  
Amber Force Field ◽  
Heterocyclic Molecules ◽  
Partial Charges ◽  
Grid Points ◽  
Dynamics Simulations ◽  
Extra Point

ABSTRACTRepresentation of electrostatic interactions by a Coulombic pair-wise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules, and thus electrostatic interactions in the force field. In particular, we re-optimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme and included extra point charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model showed clear potential for improving the performance of the nucleic acid force fields, for which poor description of nonbonded interactions, such as underestimated base pairing, makes it difficult to describe the folding free energy landscape of small oligonucleotides.

Sign in / Sign up

Export Citation Format

Share Document