Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

Accurately predicting absolute binding free energies of protein-ligand complexes is important as a fundamental problem in both computational biophysics and pharmaceutical discovery. Calculating binding free energies for charged ligands is generally considered to be challenging because of the strong electrostatic interactions between the ligand and its environment in aqueous solution. In this work, we compare the performance of the potential of mean force (PMF) method and the double decoupling method (DDM) for computing absolute binding free energies for charged ligands. We first clarify an unresolved issue concerning the explicit use of the binding site volume to define the complexed state in DDM together with the use of harmonic restraints. We also provide an alternative derivation for the formula for absolute binding free energy using the PMF approach. We use these formulas to compute the binding free energy of charged ligands at an allosteric site of HIV-1 integrase, which has emerged in recent years as a promising target for developing antiviral therapy. As compared with the experimental results, the absolute binding free energies obtained by using the PMF approach show unsigned errors of 1.5-3.4 kcal mol-1, which are somewhat better than the results from DDM (unsigned errors of 1.6-4.3 kcal mol-1) using the same amount of CPU time. According to the DDM decomposition of the binding free energy, the ligand binding appears to be dominated by nonpolar interactions despite the presence of very large and favorable intermolecular ligand-receptor electrostatic interactions, which are almost completely cancelled out by the equally large free energy cost of desolvation of the charged moiety of the ligands in solution. We discuss the relative strengths of computing absolute binding free energies using the alchemical and physical pathway methods.

[1]  A. Engelman,et al.  The mechanism of H171T resistance reveals the importance of Nδ-protonated His171 for the binding of allosteric inhibitor BI-D to HIV-1 integrase , 2014, Retrovirology.

[2]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[3]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[4]  A. Engelman,et al.  Allosteric inhibition of HIV-1 integrase activity. , 2013, Current opinion in chemical biology.

[5]  Michael K. Gilson,et al.  Overcoming dissipation in the calculation of standard binding free energies by ligand extraction , 2013, J. Comput. Chem..

[6]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[7]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[8]  Richard H. Henchman,et al.  Standard Free Energy of Binding from a One-Dimensional Potential of Mean Force. , 2009, Journal of chemical theory and computation.

[9]  David L. Mobley,et al.  Predicting binding free energies: Frontiers and benchmarks , 2016, bioRxiv.

[10]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[11]  Jennifer L. Knight,et al.  Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. , 2015, Journal of the American Chemical Society.

[12]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[13]  Ronald M. Levy,et al.  On finite‐size effects in computer simulations using the Ewald potential , 1995 .

[14]  A. Engelman,et al.  Allosteric HIV‐1 integrase inhibitors promote aberrant protein multimerization by directly mediating inter‐subunit interactions: Structural and thermodynamic modeling studies , 2016, Protein science : a publication of the Protein Society.

[15]  Benoît Roux,et al.  Calculation of Standard Binding Free Energies:  Aromatic Molecules in the T4 Lysozyme L99A Mutant. , 2006, Journal of chemical theory and computation.

[16]  T. L. Hill.  Theory of Protein Solutions. I , 1955 .

[17]  Peng Zhang,et al.  Elucidating the energetics of entropically driven protein-ligand association: calculations of absolute binding free energy and entropy. , 2011, The journal of physical chemistry. B.

[18]  Piotr Cieplak,et al.  Free energy profile of RNA hairpins: a molecular dynamics simulation study. , 2010, Biophysical journal.

[19]  Mark A Olson,et al.  Calculation of absolute protein-ligand binding affinity using path and endpoint approaches. , 2006, Biophysical journal.

[20]  Christophe Chipot,et al.  Efficient determination of protein-protein standard binding free energies from first principles. , 2013, Journal of chemical theory and computation.

[21]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[22]  M. Parrinello,et al.  Funnel metadynamics as accurate binding free-energy method , 2013, Proceedings of the National Academy of Sciences.

[23]  Emilio Gallicchio,et al.  The Binding Energy Distribution Analysis Method (BEDAM) for the Estimation of Protein-Ligand Binding Affinities. , 2010, Journal of chemical theory and computation.

[24]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[25]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[26]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[27]  Chris Oostenbrink,et al.  Net charge changes in the calculation of relative ligand-binding free energies via classical atomistic molecular dynamics simulation , 2013, J. Comput. Chem..

[28]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[29]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[30]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Christophe Chipot,et al.  Standard binding free energies from computer simulations: What is the best strategy? , 2013, Journal of chemical theory and computation.

[32]  J A McCammon,et al.  Theoretical calculation of relative binding affinity in host-guest systems. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[33]  KumarShankar,et al.  The weighted histogram analysis method for free-energy calculations on biomolecules. I , 1992 .

[34]  Emilio Gallicchio,et al.  Recent theoretical and computational advances for modeling protein-ligand binding affinities. , 2011, Advances in protein chemistry and structural biology.

[35]  Robert D Skeel,et al.  Dual role of protein phosphorylation in DNA activator/coactivator binding. , 2011, Biophysical journal.

[36]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[37]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[38]  Benoît Roux,et al.  Computational Study of Gleevec and G6G Reveals Molecular Determinants of Kinase Inhibitor Selectivity , 2014, Journal of the American Chemical Society.

[39]  Stefano Alcaro,et al.  Ligand binding to telomeric G-quadruplex DNA investigated by funnel-metadynamics simulations , 2017, Proceedings of the National Academy of Sciences.

[40]  David L Mobley,et al.  Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: an accurate correction scheme for electrostatic finite-size effects. , 2013, The Journal of chemical physics.

[41]  Arthur J. Olson,et al.  Virtual screening of integrase inhibitors by large scale binding free energy calculations: the SAMPL4 challenge , 2014, Journal of Computer-Aided Molecular Design.

[42]  Bin W. Zhang,et al.  The Role of Interfacial Water in Protein-Ligand Binding: Insights from the Indirect Solvent Mediated Potential of Mean Force. , 2018, Journal of chemical theory and computation.

[43]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[44]  A. Engelman,et al.  The mechanism of H171T resistance reveals the importance of Nδ-protonated His171 for the binding of allosteric inhibitor BI-D to HIV-1 integrase , 2014, Retrovirology.

[45]  David D. L. Minh,et al.  Absolute Binding Free Energies between T4 Lysozyme and 141 Small Molecules: Calculations Based on Multiple Rigid Receptor Configurations. , 2017, Journal of chemical theory and computation.

[46]  Anthony K. Felts,et al.  Temperature weighted histogram analysis method, replica exchange, and transition paths. , 2005, The journal of physical chemistry. B.

[47]  Niel M. Henriksen,et al.  Attach-Pull-Release Calculations of Ligand Binding and Conformational Changes on the First BRD4 Bromodomain , 2017, Journal of chemical theory and computation.

[48]  Di Cui,et al.  Free energetics of rigid body association of ubiquitin binding domains: A biochemical model for binding mediated by hydrophobic interaction , 2014, Proteins.