On Restraints in End‐Point Protein–Ligand Binding Free Energy Calculations

The impact of harmonic restraints on protein heavy atoms and ligand atoms on end‐point free energy calculations is systematically characterized for 54 protein–ligand complexes. We observe that stronger restraints reduce the equilibration time and statistical inefficiency, suppress conformational sampling, influence correlation with experiment, and monotonically decrease the estimated loss of entropy upon binding, leading to stronger estimated binding free energies in most systems. A statistical estimator that reweights for the biasing potential and includes data prior to the estimated equilibration time has the highest correlation with experiment. A spring constant of 20 cal mol−1 Å−2 maintains a near‐native energy landscape and suppresses artifactual energy minima while minimally limiting thermal fluctuations about the crystal structure. © 2019 Wiley Periodicals, Inc.

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

[2]  Ivet Bahar,et al.  ProDy: Protein Dynamics Inferred from Theory and Experiments , 2011, Bioinform..

[3]  Chris Oostenbrink,et al.  Improved ligand-protein binding affinity predictions using multiple binding modes. , 2010, Biophysical journal.

[4]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

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

[6]  C. David Andersson,et al.  Postprocessing of Docked Protein-Ligand Complexes Using Implicit Solvation Models , 2011, J. Chem. Inf. Model..

[7]  Wei Jiang,et al.  CHARMM-GUI Ligand Binder for Absolute Binding Free Energy Calculations and Its Application , 2013, J. Chem. Inf. Model..

[8]  Thomas J Lane,et al.  MDTraj: a modern, open library for the analysis of molecular dynamics trajectories , 2014, bioRxiv.

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

[10]  Alessandro Laio,et al.  Clustering by fast search and find of density peaks , 2014, Science.

[11]  Lingle Wang,et al.  On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities , 2012, Proceedings of the National Academy of Sciences.

[12]  P A Kollman,et al.  Free energy calculations on dimer stability of the HIV protease using molecular dynamics and a continuum solvent model. , 2000, Journal of molecular biology.

[13]  N. Skrynnikov,et al.  Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics , 2014, Protein science : a publication of the Protein Society.

[14]  John Z. H. Zhang,et al.  Interaction Entropy: A New Paradigm for Highly Efficient and Reliable Computation of Protein-Ligand Binding Free Energy. , 2016, Journal of the American Chemical Society.

[15]  David D. L Minh.  Alchemical Grid Dock (AlGDock): Binding Free Energy Calculations between Flexible Ligands and Rigid Receptors , 2020, J. Comput. Chem..

[16]  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.

[17]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[18]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

[19]  David L Mobley,et al.  Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site. , 2009, Journal of molecular biology.

[20]  A. Atilgan,et al.  Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. , 1997, Folding & design.

[21]  Kai Wang,et al.  Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics , 2013, Journal of Computer-Aided Molecular Design.

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

[23]  Arieh Warshel,et al.  A comprehensive examination of the contributions to the binding entropy of protein–ligand complexes , 2010, Proteins.

[24]  M. Kendall.  A NEW MEASURE OF RANK CORRELATION , 1938 .

[25]  Jacob Kongsted,et al.  An improved method to predict the entropy term with the MM/PBSA approach , 2009, J. Comput. Aided Mol. Des..

[26]  Stefano Piana,et al.  Refinement of protein structure homology models via long, all‐atom molecular dynamics simulations , 2012, Proteins.

[27]  Nathan A. Baker,et al.  PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations , 2004, Nucleic Acids Res..

[28]  J. Essex,et al.  Hit identification and binding mode predictions by rigorous free energy simulations. , 2008, Journal of medicinal chemistry.

[29]  David D L Minh,et al.  The entropic cost of protein-protein association: a case study on acetylcholinesterase binding to fasciculin-2. , 2005, Biophysical journal.

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

[31]  Wilfred F. van Gunsteren,et al.  Free energy calculations offer insights into the influence of receptor flexibility on ligand–receptor binding affinities , 2011, J. Comput. Aided Mol. Des..

[32]  David L Mobley,et al.  Perspective: Alchemical free energy calculations for drug discovery. , 2012, The Journal of chemical physics.

[33]  Ross C. Walker,et al.  An overview of the Amber biomolecular simulation package , 2013 .

[34]  B. Roux,et al.  Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis. , 2018, The journal of physical chemistry. B.

[35]  Michael J. Bodkin,et al.  Accurate calculation of the absolute free energy of binding for drug molecules , 2015, Chemical science.

[36]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[37]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[38]  C. Chipot,et al.  New Coarse Variables for the Accurate Determination of Standard Binding Free Energies. , 2017, Journal of chemical theory and computation.

[39]  Julien Michel,et al.  Reproducibility of Free Energy Calculations across Different Molecular Simulation Software Packages. , 2018, Journal of chemical theory and computation.

[40]  K. Dill,et al.  On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. , 2006, The Journal of chemical physics.

[41]  Charlotte M. Deane,et al.  Rapid, Precise, and Reproducible Prediction of Peptide-MHC Binding Affinities from Molecular Dynamics That Correlate Well with Experiment. , 2015, Journal of chemical theory and computation.

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

[43]  David L. Mobley,et al.  Let’s get honest about sampling , 2011, Journal of Computer-Aided Molecular Design.

[44]  Vijay S. Pande,et al.  OpenMM: A Hardware-Independent Framework for Molecular Simulations , 2010, Computing in Science & Engineering.

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

[46]  Robert H. Wood,et al.  Systematic errors in free energy perturbation calculations due to a finite sample of configuration space: sample-size hysteresis , 1991 .

[47]  David D. L. Minh,et al.  Simple Entropy Terms for End-Point Binding Free Energy Calculations. , 2018, Journal of chemical theory and computation.

[48]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[49]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[50]  Philip C Biggin,et al.  Predictions of Ligand Selectivity from Absolute Binding Free Energy Calculations , 2016, Journal of the American Chemical Society.

[51]  Michael K Gilson,et al.  Evaluating the Accuracy of the Quasiharmonic Approximation. , 2005, Journal of chemical theory and computation.

[52]  Chris Oostenbrink,et al.  Simulation of Reversible Protein–Protein Binding and Calculation of Binding Free Energies Using Perturbed Distance Restraints , 2017, Journal of chemical theory and computation.

[53]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[54]  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.

[55]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[56]  Edward D Harder,et al.  How To Deal with Multiple Binding Poses in Alchemical Relative Protein–Ligand Binding Free Energy Calculations , 2015, Journal of chemical theory and computation.

[57]  P. Coveney,et al.  Application of ESMACS binding free energy protocols to diverse datasets: Bromodomain-containing protein 4 , 2018, Scientific Reports.

[58]  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.

[59]  Benoît Roux,et al.  Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer. , 2012, Journal of chemical theory and computation.

[60]  Serdar Kuyucak,et al.  Accurate determination of the binding free energy for KcsA-charybdotoxin complex from the potential of mean force calculations with restraints. , 2011, Biophysical journal.

[61]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[62]  Benoît Roux,et al.  Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations. , 2010, Journal of chemical theory and computation.

[63]  Peter V. Coveney,et al.  Accurate Ensemble Molecular Dynamics Binding Free Energy Ranking of Multidrug-Resistant HIV-1 Proteases , 2010, J. Chem. Inf. Model..

[64]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

[65]  Peter V Coveney,et al.  Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study , 2016, Journal of chemical theory and computation.

[66]  U. Ryde,et al.  Ligand affinities predicted with the MM/PBSA method: dependence on the simulation method and the force field. , 2006, Journal of Medicinal Chemistry.

[67]  Benoît Roux,et al.  Computational analysis of the binding specificity of Gleevec to Abl, c-Kit, Lck, and c-Src tyrosine kinases. , 2013, Journal of the American Chemical Society.

[68]  Benoît Roux,et al.  Calculation of the standard binding free energy of sparsomycin to the ribosomal peptidyl‐transferase P‐site using molecular dynamics simulations with restraining potentials , 2010, Journal of molecular recognition : JMR.

[69]  A. Elcock,et al.  Accurate calculation of mutational effects on the thermodynamics of inhibitor binding to p38α MAP kinase: a combined computational and experimental study. , 2013, Journal of chemical theory and computation.

[70]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[71]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[72]  B. Roux,et al.  Explaining why Gleevec is a specific and potent inhibitor of Abl kinase , 2013, Proceedings of the National Academy of Sciences.

[73]  Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations. , 2019, Journal of the American Chemical Society.

[74]  Vijay S. Pande,et al.  OpenMM 7: Rapid development of high performance algorithms for molecular dynamics , 2016, bioRxiv.

[75]  William B. Smith,et al.  Selective inhibition of BET bromodomains , 2010, Nature.

[76]  Luigi Capoferri,et al.  Improving the iterative Linear Interaction Energy approach using automated recognition of configurational transitions , 2016, Journal of Molecular Modeling.

[77]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

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

[79]  Benoît Roux,et al.  Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations. , 2008, The Journal of chemical physics.

[80]  John D Chodera,et al.  A Simple Method for Automated Equilibration Detection in Molecular Simulations. , 2016, Journal of chemical theory and computation.

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

[82]  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.

[83]  F. Lightstone,et al.  A Comprehensive Docking and MM/GBSA Rescoring Study of Ligand Recognition upon Binding Antithrombin , 2017, Current topics in medicinal chemistry.

[84]  B. Roux,et al.  Absolute binding free energy calculations of sparsomycin analogs to the bacterial ribosome. , 2010, The journal of physical chemistry. B.