scholarly journals Evaluating the use of absolute binding free energy in the fragment optimization process.

Author(s):  
Irfan Alibay ◽  
Aniket Mangakar ◽  
Daniel Seeliger ◽  
Philip Biggin

Key to the fragment optimization process is the need to accurately capture the changes in affinity that are associated with a given set of chemical modifications. Due to the weakly binding nature of fragments, this has proven to be a challenging task, despite recent advancements in leveraging experimental and computational methods. In this work, we evaluate the use of Absolute Binding Free Energy (ABFE) calculations in guiding fragment optimization decisions, retrospectively calculating binding free energies for 59 ligands across 4 fragment elaboration campaigns. We first demonstrate that ABFEs can be used to accurately rank fragment-sized binders with an overall Spearman’s r of 0.89 and a Kendall τ of 0.67, although often deviating from experiment in absolute free energy values with an RMSE of 2.75 kcal/mol. We then also show that in several cases, retrospective fragment optimization decisions can be supported by the ABFE calculations. Cases that were not supported were often limited by large uncertainties in the free energy estimates, however generally the right direction in ΔΔG is still observed. Comparing against cheaper endpoint methods, namely Nwat-MM/GBSA, we find that ABFEs offer better outcomes in ranking binders, improving correlation metrics, although a similar confidence in retrospective synthetic decisions is achieved. Our results indicate that ABFE calculations are currently at the level of accuracy that can be usefully employed to gauge which fragment elaborations are likely to offer the best gains in affinity.

2018 ◽  
Vol 20 (25) ◽  
pp. 17081-17092 ◽  
Author(s):  
Nanjie Deng ◽  
Di Cui ◽  
Bin W. Zhang ◽  
Junchao Xia ◽  
Jeffrey Cruz ◽  
...  

We compare the performance of the potential of mean force (PMF) method and double decoupling method (DDM) for computing absolute binding free energies for charged ligands.


Author(s):  
Germano Heinzelmann ◽  
Michael K. Gilson

AbstractAbsolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to fully automate the calculation of binding free energies for a protein with a series of ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost allows this procedure to be applied in a relatively high-throughput mode, and thus enables new applications in early-stage drug discovery.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Germano Heinzelmann ◽  
Michael K. Gilson

AbstractAbsolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.


2021 ◽  
Author(s):  
Yuriy Khalak ◽  
Gary Tresdern ◽  
Matteo Aldeghi ◽  
Hannah Magdalena Baumann ◽  
David L. Mobley ◽  
...  

The recent advances in relative protein-ligand binding free energy calculations have shown the value of alchemical methods in drug discovery. Accurately assessing absolute binding free energies, although highly desired, remains...


2021 ◽  
Author(s):  
Alexander Wade ◽  
Agastya Bhati ◽  
Shunzhou Wan ◽  
Peter Coveney

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat myriad diseases. In this work we examine the computation of alchemical relative binding free energies with an eye to assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2 and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2 and NAMD3. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between packages of 0.5 $kcal/mol$ The correlation between packages is very good with the lowest Spearman's, Pearson's and Kendall's tau correlation coefficient between two packages being 0.91, 0.89 and 0.74 respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.


Author(s):  
Mahdi Ghorbani ◽  
Phillip S. Hudson ◽  
Michael R. Jones ◽  
Félix Aviat ◽  
Rubén Meana-Pañeda ◽  
...  

AbstractIn this study, we report binding free energy calculations of various drugs-of-abuse to Cucurbit-[8]-uril as part of the SAMPL8 blind challenge. Force-field parameters were obtained from force-matching with different quantum mechanical levels of theory. The Replica Exchange Umbrella Sampling (REUS) approach was used with a cylindrical restraint to enhance the sampling of host–guest binding. Binding free energy was calculated by pulling the guest molecule from one side of the symmetric and cylindrical host, then into and through the host, and out the other side (bidirectional) as compared to pulling only to the bound pose inside the cylindrical host (unidirectional). The initial results with force-matched MP2 parameter set led to RMSE of 4.68 $${\text{kcal}}/{\text{mol}}$$ kcal / mol from experimental values. However, the follow-up study with CHARMM generalized force field parameters and force-matched PM6-D3H4 parameters resulted in RMSEs from experiment of $$2.65$$ 2.65 and $$1.72 {\text{kcal}}/{\text{mol}}$$ 1.72 kcal / mol , respectively, which demonstrates the potential of REUS for accurate binding free energy calculation given a more suitable description of energetics. Moreover, we compared the free energies for the so called bidirectional and unidirectional free energy approach and found that the binding free energies were highly similar. However, one issue in the bidirectional approach is the asymmetry of profile on the two sides of the host. This is mainly due to the insufficient sampling for these larger systems and can be avoided by longer sampling simulations. Overall REUS shows great promise for binding free energy calculations.


2021 ◽  
Vol 4 (1) ◽  
Author(s):  
Vytautas Gapsys ◽  
Ahmet Yildirim ◽  
Matteo Aldeghi ◽  
Yuriy Khalak ◽  
David van der Spoel ◽  
...  

AbstractThe accurate calculation of the binding free energy for arbitrary ligand–protein pairs is a considerable challenge in computer-aided drug discovery. Recently, it has been demonstrated that current state-of-the-art molecular dynamics (MD) based methods are capable of making highly accurate predictions. Conventional MD-based approaches rely on the first principles of statistical mechanics and assume equilibrium sampling of the phase space. In the current work we demonstrate that accurate absolute binding free energies (ABFE) can also be obtained via theoretically rigorous non-equilibrium approaches. Our investigation of ligands binding to bromodomains and T4 lysozyme reveals that both equilibrium and non-equilibrium approaches converge to the same results. The non-equilibrium approach achieves the same level of accuracy and convergence as an equilibrium free energy perturbation (FEP) method enhanced by Hamiltonian replica exchange. We also compare uni- and bi-directional non-equilibrium approaches and demonstrate that considering the work distributions from both forward and reverse directions provides substantial accuracy gains. In summary, non-equilibrium ABFE calculations are shown to yield reliable and well-converged estimates of protein–ligand binding affinity.


Author(s):  
David Huggins

<p>We present an approach to performing alchemical binding free energies which we term coupled topologies. Simultaneously coupling a molecule in the bound state while decoupling it in the unbound state allows us to calculate free energy changes where the system changes charge, without the need to correct for simulation artifacts. This solves a longstanding problem in computing free energy changes. The approach is applied to separated topology relative binding free energy calculations, but is appropriate for single topology calculations and dual topology calculations as well as absolute binding free energy calculations. We apply the method to small-molecule inhibitors of AmpC β-lactamase and show the coupled topologies approach yields results that are in excellent agreement with experiment and good agreement with a state-of-the-art separated topology approach. The promising results on this test case suggest that the coupled topologies approach will be a useful addition to the available arsenal of free-energy methods.</p>


2000 ◽  
Vol 47 (1) ◽  
pp. 1-9 ◽  
Author(s):  
W R Rudnicki ◽  
M Kurzepa ◽  
T Szczepanik ◽  
W Priebe ◽  
B Lesyng

A theoretical model for predicting the free energy of binding between anthracycline antibiotics and DNA was developed using the electron density functional (DFT) and molecular mechanics (MM) methods. Partial DFT-ESP charges were used in calculating the MM binding energies for complexes formed between anthracycline antibiotics and oligodeoxynucleotides. These energies were then compared with experimental binding free energies. The good correlation between the experimental and theoretical energies allowed us to propose a model for predicting the binding free energy for derivatives of anthracycline antibiotics and for quickly screening new anthracycline derivatives.


Author(s):  
Hafiz Saqib Ali ◽  
Arghya Chakravorty ◽  
Jas Kalayan ◽  
Samuel P. de Visser ◽  
Richard H. Henchman

AbstractFree energy drives a wide range of molecular processes such as solvation, binding, chemical reactions and conformational change. Given the central importance of binding, a wide range of methods exist to calculate it, whether based on scoring functions, machine-learning, classical or electronic structure methods, alchemy, or explicit evaluation of energy and entropy. Here we present a new energy–entropy (EE) method to calculate the host–guest binding free energy directly from molecular dynamics (MD) simulation. Entropy is evaluated using Multiscale Cell Correlation (MCC) which uses force and torque covariance and contacts at two different length scales. The method is tested on a series of seven host–guest complexes in the SAMPL8 (Statistical Assessment of the Modeling of Proteins and Ligands) “Drugs of Abuse” Blind Challenge. The EE-MCC binding free energies are found to agree with experiment with an average error of 0.9 kcal mol−1. MCC makes clear the origin of the entropy changes, showing that the large loss of positional, orientational, and to a lesser extent conformational entropy of each binding guest is compensated for by a gain in orientational entropy of water released to bulk, combined with smaller decreases in vibrational entropy of the host, guest and contacting water.


Sign in / Sign up

Export Citation Format

Share Document