Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins

2003 ◽  
Vol 119 (11) ◽  
pp. 5740-5761 ◽  
Author(s):  
Michael R. Shirts ◽  
Jed W. Pitera ◽  
William C. Swope ◽  
Vijay S. Pande
Keyword(s):  
Free Energy ◽  
Amino Acid ◽  
Molecular Mechanics ◽  
Force Fields ◽  
Free Energy Calculations ◽  
Side Chain ◽  
Amino Acid Side Chain ◽  
Energy Calculations

scholarly journals Towards chemical accuracy for alchemical free energy calculations with hybrid physics-based machine learning / molecular mechanics potentials

2020 ◽  
Author(s):  
Dominic A. Rufa ◽  
Hannah E. Bruce Macdonald ◽  
Josh Fass ◽  
Marcus Wieder ◽  
Patrick B. Grinaway ◽  
...  
Keyword(s):  
Machine Learning ◽  
Free Energy ◽  
Molecular Mechanics ◽  
Force Fields ◽  
Free Energy Calculations ◽  
Binding Affinities ◽  
Free Energies ◽  
Energy Calculations ◽  
Alchemical Free Energy ◽  
Binding Free Energies

AbstractAlchemical free energy methods with molecular mechanics (MM) force fields are now widely used in the prioritization of small molecules for synthesis in structure-enabled drug discovery projects because of their ability to deliver 1–2 kcal mol−1 accuracy in well-behaved protein-ligand systems. Surpassing this accuracy limit would significantly reduce the number of compounds that must be synthesized to achieve desired potencies and selectivities in drug design campaigns. However, MM force fields pose a challenge to achieving higher accuracy due to their inability to capture the intricate atomic interactions of the physical systems they model. A major limitation is the accuracy with which ligand intramolecular energetics—especially torsions—can be modeled, as poor modeling of torsional profiles and coupling with other valence degrees of freedom can have a significant impact on binding free energies. Here, we demonstrate how a new generation of hybrid machine learning / molecular mechanics (ML/MM) potentials can deliver significant accuracy improvements in modeling protein-ligand binding affinities. Using a nonequilibrium perturbation approach, we can correct a standard, GPU-accelerated MM alchemical free energy calculation in a simple post-processing step to efficiently recover ML/MM free energies and deliver a significant accuracy improvement with small additional computational effort. To demonstrate the utility of ML/MM free energy calculations, we apply this approach to a benchmark system for predicting kinase:inhibitor binding affinities—a congeneric ligand series for non-receptor tyrosine kinase TYK2 (Tyk2)—wherein state-of-the-art MM free energy calculations (with OPLS2.1) achieve inaccuracies of 0.93±0.12 kcal mol−1 in predicting absolute binding free energies. Applying an ML/MM hybrid potential based on the ANI2x ML model and AMBER14SB/TIP3P with the OpenFF 1.0.0 (“Parsley”) small molecule force field as an MM model, we show that it is possible to significantly reduce the error in absolute binding free energies from 0.97 [95% CI: 0.68, 1.21] kcal mol−1 (MM) to 0.47 [95% CI: 0.31, 0.63] kcal mol−1 (ML/MM).

scholarly journals Challenges Encountered Applying Equilibrium and Non-Equilibrium Binding Free Energy Calculations

2020 ◽  
Author(s):  
Hannah Baumann ◽  
Vytautas Gapsys ◽  
Bert L. de Groot ◽  
David Mobley
Keyword(s):  
Free Energy ◽  
Binding Free Energy ◽  
Free Energy Calculations ◽  
Side Chain ◽  
Future Research ◽  
Energy Calculations ◽  
Binding Free Energy Calculations ◽  
Sampling Problems ◽  
Alchemical Free Energy ◽  
Non Equilibrium

<div>Binding free energy calculations have become increasingly valuable to drive decision making in drug discovery projects. </div><div>However, among other issues, inadequate sampling can reduce accuracy, limiting the value of the technique.</div><div>In this paper we apply absolute binding free energy calculations to ligands binding to T4 lysozyme L99A and HSP90 using equilibrium and non-equilibrium approaches. We highlight sampling problems encountered in these systems, such as slow side chain rearrangements and slow changes of water placement upon ligand binding. These same types of challenges are likely to show up in other protein-ligand systems as well and we propose some strategies to diagnose and test for such problems in alchemical free energy calculations. We also explore similarities and differences in how the equilibrium and the non-equilibrium approaches handle these problems. Our results show the large amount of work still to be done to make free energy calculations robust and reliable and provide insight for future research in this area. </div>

Insights to the Binding of a Selective Adenosine A3 Receptor Antagonist using Molecular Dynamics Simulations, Binding Free Energy Calculations and Mutagenesis

2019 ◽  
Author(s):  
Panagiotis Lagarias ◽  
Kerry Barkan ◽  
Eva Tzortzini ◽  
Eleni Vrontaki ◽  
Margarita Stampelou ◽  
...  
Keyword(s):  
Free Energy ◽  
Computational Model ◽  
Molecular Mechanics ◽  
Surface Area ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Free Energy Calculations ◽  
Binding Energies ◽  
Energy Calculations ◽  
Binding Free Energy Calculations

<p>Adenosine A<sub>3 </sub>receptor (A<sub>3</sub>R), is a promising drug target against cancer cell proliferation. Currently there is no experimentally determined structure of A<sub>3</sub>R. Here, we have investigate a computational model, previously applied successfully for agonists binding to A<sub>3</sub>R, using molecular dynamic (MD) simulations, Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) binding free energy calculations. Extensive computations were performed to explore the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18) to A<sub>3</sub>R. K18 is a new specific and competitive antagonist at the orthosteric binding site of A<sub>3</sub>R, discovered using virtual screening and characterized pharmacologically in our previous studies. The most plausible binding conformation for the dichlorophenyl group of K18 inside the A<sub>3</sub>R is oriented towards trans-membrane helices (TM) 5 and 6, according to the MM-PBSA and MM-GBSA binding free energy calculations, and by the previous results obtained by mutating residues of TM5, TM6 to alanine which reduce antagonist potency. The results from 14 site-directed mutagenesis experiments were interpreted using MD simulations and MM-GBSA calculations which show that the relative binding free energies of the mutant A<sub>3</sub>R - K18 complexes compare to the WT A<sub>3</sub>R are in agreement with the effect of the mutations, i.e. the reduction, maintenance or increase of antagonist potency. We show that when the residues V169<sup>5.30</sup>, M177<sup>5.38</sup>, I249<sup>6.54</sup> involved in direct interactions with K18 are mutated to alanine, the mutant A<sub>3</sub>R - K18 complexes reduce potency, increase the RMSD value of K18 inside the binding area and the MM-GBSA binding free energy compared to the WT A<sub>3</sub>R complex. Our computational model shows that other mutant A<sub>3</sub>R complexes with K18, including directly interacting residues, i.e. F168<sup>5.29</sup>A, L246<sup>6.51</sup>A, N250<sup>6.55</sup>A complexes with K18 are not stable. In these complexes of A<sub>3</sub>R mutated in directly interacting residues one or more of the interactions between K18 and these residues are lost. In agreement with the experiments, the computations show that, M174<sup>5.35</sup> a residue which does not make direct interactions with K18 is critical for K18 binding. A striking results is that the mutation of residue V169<sup>5.30</sup> to glutamic acid maintained antagonistic potency. This effect is in agreement with the binding free energy calculations and it is suggested that is due to K18 re-orientation but also to the plasticity of A<sub>3</sub>R binding area. The mutation of direct interacting L90<sup>3.32</sup> in the low region and the non-directly interacting L264<sup>7.35</sup> to alanine in the middle region increases the antagonistic potency, suggesting that chemical modifications of K18 can be applied to augment antagonistic potency. The calculated binding energies Δ<i>G</i><sub>eff</sub> values of K18 against mutant A<sub>3</sub>Rs displayed very good correlation with experimental potencies (pA<sub>2</sub> values). These results further approve the computational model for the description of K18 binding with critical residues of the orthosteric binding area which can have implications for the design of more effective antagonists based on the structure of K18.</p>

Sign in / Sign up

Export Citation Format

Share Document