Estimation of kinetic and thermodynamic ligand-binding parameters using computational strategies.

Kinetic and thermodynamic ligand-protein binding parameters are gaining growing importance as key information to consider in drug discovery. The determination of the molecular structures, using particularly x-ray and NMR techniques, is crucial for understanding how a ligand recognizes its target in the final binding complex. However, for a better understanding of the recognition processes, experimental studies of ligand-protein interactions are needed. Even though several techniques can be used to investigate both thermodynamic and kinetic profiles for a ligand-protein complex, these procedures are very often laborious, time consuming and expensive. In the last 10 years, computational approaches have enormous potential in providing insights into each of the above effects and in parsing their contributions to the changes in both kinetic and thermodynamic binding parameters. The main purpose of this review is to summarize the state of the art of computational strategies for estimating the kinetic and thermodynamic parameters of a ligand-protein binding.

[1]  H. Eyring The Activated Complex in Chemical Reactions , 1935 .

[2]  M. Perutz,et al.  Structure of Hæmoglobin: A Three-Dimensional Fourier Synthesis at 5.5-Å. Resolution, Obtained by X-Ray Analysis , 1960, Nature.

[3]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[4]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[5]  K. Laidler,et al.  Development of transition-state theory , 1983 .

[6]  W. Gilbert Origin of life: The RNA world , 1986, Nature.

[7]  P. Wolynes,et al.  Intermediates and barrier crossing in a random energy model , 1989 .

[8]  R. Wade,et al.  New hydrogen-bond potentials for use in determining energetically favorable binding sites on molecules of known structure. , 1989, Journal of medicinal chemistry.

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

[10]  PatrickY.-S. Lam,et al.  Rational design of potent, bioavailable, nonpeptide cyclic ureas as HIV protease inhibitors. , 1994, Science.

[11]  J. Dunitz The entropic cost of bound water in crystals and biomolecules. , 1994, Science.

[12]  P. Tavan,et al.  Ligand Binding: Molecular Mechanics Calculation of the Streptavidin-Biotin Rupture Force , 1996, Science.

[13]  J. Ladbury Just add water! The effect of water on the specificity of protein-ligand binding sites and its potential application to drug design. , 1996, Chemistry & biology.

[14]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[15]  J. Mccammon,et al.  Electrostatic Influence on the Kinetics of Ligand Binding to Acetylcholinesterase , 1997, The Journal of Biological Chemistry.

[16]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[17]  I. Kuntz,et al.  The maximal affinity of ligands. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[18]  K. Dill,et al.  The Strength of Hydrogen Bonds in Liquid Water and Around Nonpolar Solutes , 2000 .

[19]  R. Nussinov,et al.  Folding and binding cascades: Dynamic landscapes and population shifts , 2008, Protein science : a publication of the Protein Society.

[20]  G. Gilli,et al.  Towards an unified hydrogen-bond theory , 2000 .

[21]  R. Norel,et al.  Electrostatic aspects of protein-protein interactions. , 2000, Current opinion in structural biology.

[22]  R. Wade,et al.  How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms. , 2000, Journal of molecular biology.

[23]  U. Gether Uncovering molecular mechanisms involved in activation of G protein-coupled receptors. , 2000, Endocrine reviews.

[24]  K. Schulten,et al.  Steered molecular dynamics and mechanical functions of proteins. , 2001, Current opinion in structural biology.

[25]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[26]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[27]  V. Hruby Designing peptide receptor agonists and antagonists , 2002, Nature Reviews Drug Discovery.

[28]  Early mechanistic events in biotin dissociation from streptavidin , 2002, Nature Structural Biology.

[29]  L. Tong,et al.  Inhibition of p38 MAP kinase by utilizing a novel allosteric binding site , 2002, Nature Structural Biology.

[30]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[31]  Marcel L. Verdonk,et al.  The consequences of translational and rotational entropy lost by small molecules on binding to proteins , 2002, J. Comput. Aided Mol. Des..

[32]  Gerhard Klebe,et al.  Utilising structural knowledge in drug design strategies: applications using Relibase. , 2003, Journal of molecular biology.

[33]  Xiaomin Luo,et al.  Steered molecular dynamics simulation on the binding of NNRTI to HIV-1 RT. , 2003, Biophysical journal.

[34]  Richard H. Henchman,et al.  The dynamics of ligand barrier crossing inside the acetylcholinesterase gorge. , 2003, Biophysical journal.

[35]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[36]  M. Karplus,et al.  How Enzymes Work: Analysis by Modern Rate Theory and Computer Simulations , 2004, Science.

[37]  Johannes C. Hermann,et al.  A combined QM/MM approach to protein--ligand interactions: polarization effects of the HIV-1 protease on selected high affinity inhibitors. , 2004, Journal of medicinal chemistry.

[38]  J. Bajorath,et al.  Docking and scoring in virtual screening for drug discovery: methods and applications , 2004, Nature Reviews Drug Discovery.

[39]  A. Laio,et al.  Flexible docking in solution using metadynamics. , 2005, Journal of the American Chemical Society.

[40]  Jeremy C. Smith,et al.  Protein/ligand binding free energies calculated with quantum mechanics/molecular mechanics. , 2005, The journal of physical chemistry. B.

[41]  A. Laio,et al.  Assessing the accuracy of metadynamics. , 2005, The journal of physical chemistry. B.

[42]  G. Hummer,et al.  Reaction coordinates and rates from transition paths. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[43]  J. Gready,et al.  Combining docking and molecular dynamic simulations in drug design , 2006, Medicinal research reviews.

[44]  G. Klebe Virtual ligand screening: strategies, perspectives and limitations , 2006, Drug Discovery Today.

[45]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[46]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

[47]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[48]  John P. Overington,et al.  How many drug targets are there? , 2006, Nature Reviews Drug Discovery.

[49]  A. Wand,et al.  Characterization of the fast dynamics of protein amino acid side chains using NMR relaxation in solution. , 2006, Chemical reviews.

[50]  R. Copeland,et al.  Drug–target residence time and its implications for lead optimization , 2007, Nature Reviews Drug Discovery.

[51]  P. V. Hippel From "Simple" DNA-Protein Interactions to the Macromolecular Machines of Gene Expression , 2007 .

[52]  A. Laio,et al.  A bias-exchange approach to protein folding. , 2007, The journal of physical chemistry. B.

[53]  S. Moro,et al.  The ATP‐Binding Site of Protein Kinase CK2 Holds a Positive Electrostatic Area and Conserved Water Molecules , 2007, Chembiochem : a European journal of chemical biology.

[54]  A. Wand,et al.  Conformational entropy in molecular recognition by proteins , 2007, Nature.

[55]  R. Friesner,et al.  Structural analysis of protein dynamics by MD simulations and NMR spin‐relaxation , 2008, Proteins.

[56]  M. Peräkylä Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations , 2009, European Biophysics Journal.

[57]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[58]  P. Ball Water as an active constituent in cell biology. , 2008, Chemical reviews.

[59]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[60]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[61]  B. Brooks,et al.  Multiscale methods for macromolecular simulations. , 2008, Current opinion in structural biology.

[62]  Lina M. Nilsson,et al.  Beyond induced-fit receptor-ligand interactions: structural changes that can significantly extend bond lifetimes. , 2008, Structure.

[63]  J. Šponer,et al.  Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM. , 2009, Methods.

[64]  T. Oas,et al.  Conformational selection or induced fit: A flux description of reaction mechanism , 2009, Proceedings of the National Academy of Sciences.

[65]  A. Laio,et al.  Substrate binding mechanism of HIV-1 protease from explicit-solvent atomistic simulations. , 2009, Journal of the American Chemical Society.

[66]  Kateri H. DuBay,et al.  Calculation of proteins' total side-chain torsional entropy and its influence on protein-ligand interactions. , 2009, Journal of molecular biology.

[67]  N. Vermeulen,et al.  The role of water molecules in computational drug design. , 2010, Current topics in medicinal chemistry.

[68]  P. Tonge,et al.  Drug-target residence time: critical information for lead optimization. , 2010, Current opinion in chemical biology.

[69]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[70]  C. Higgs,et al.  Hydration Site Thermodynamics Explain SARs for Triazolylpurines Analogues Binding to the A2A Receptor. , 2010, ACS medicinal chemistry letters.

[71]  M. Trevethick,et al.  Ligand binding assays at equilibrium: validation and interpretation , 2010, British journal of pharmacology.

[72]  A. Cavalli,et al.  Molecular basis of cyclooxygenase enzymes (COXs) selective inhibition , 2010, Proceedings of the National Academy of Sciences.

[73]  Albert C. Pan,et al.  Pathway and mechanism of drug binding to G-protein-coupled receptors , 2011, Proceedings of the National Academy of Sciences.

[74]  F. J. Luque,et al.  Shielded hydrogen bonds as structural determinants of binding kinetics: application in drug design. , 2011, Journal of the American Chemical Society.

[75]  Angelo D. Favia,et al.  Substrate Binding Process and Mechanistic Functioning of Type 1 11β-Hydroxysteroid Dehydrogenase from Enhanced Sampling Methods , 2011, PloS one.

[76]  Claudio N. Cavasotto,et al.  Quantum mechanical binding free energy calculation for phosphopeptide inhibitors of the Lck SH2 domain , 2011, J. Comput. Chem..

[77]  M. Holloway,et al.  Thermodynamics of ligand binding and efficiency. , 2011, ACS medicinal chemistry letters.

[78]  Claudio N. Cavasotto,et al.  Computational and experimental studies of the interaction between phospho-peptides and the C-terminal domain of BRCA1 , 2011, J. Comput. Aided Mol. Des..

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

[80]  F. Diederich,et al.  Aromatic rings in chemical and biological recognition: energetics and structures. , 2011, Angewandte Chemie.

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

[82]  J. Kästner Umbrella sampling , 2011 .

[83]  M. Parrinello,et al.  Metadynamics with Adaptive Gaussians. , 2012, Journal of chemical theory and computation.

[84]  F. Noé,et al.  Calculating kinetics and pathways of protein-ligand association. , 2012, European journal of cell biology.

[85]  Gregory A Ross,et al.  Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites , 2012, PloS one.

[86]  Levi C. T. Pierce,et al.  Routine Access to Millisecond Time Scale Events with Accelerated Molecular Dynamics , 2012, Journal of chemical theory and computation.

[87]  Samuel Genheden,et al.  A semiempirical approach to ligand‐binding affinities: Dependence on the Hamiltonian and corrections , 2012, J. Comput. Chem..

[88]  W. L. Jorgensen,et al.  Methyl effects on protein-ligand binding. , 2012, Journal of medicinal chemistry.

[89]  T. Woolf,et al.  Towards the prediction of order parameters from molecular dynamics simulations in proteins. , 2012, The Journal of chemical physics.

[90]  K. D. Collins Why continuum electrostatics theories cannot explain biological structure, polyelectrolytes or ionic strength effects in ion-protein interactions. , 2012, Biophysical chemistry.

[91]  P. Seeman,et al.  Clozapine, atypical antipsychotics, and the benefits of fast-off D2 dopamine receptor antagonism , 2012, Naunyn-Schmiedeberg's Archives of Pharmacology.

[92]  E. Di Cera,et al.  Conformational selection or induced fit? A critical appraisal of the kinetic mechanism. , 2012, Biochemistry.

[93]  M. Parrinello,et al.  Sampling protein motion and solvent effect during ligand binding , 2012, Proceedings of the National Academy of Sciences.

[94]  J. Mccammon,et al.  Population Based Reweighting of Scaled Molecular Dynamics , 2013, The journal of physical chemistry. B.

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

[96]  Jianhan Chen,et al.  Efficiency of Adaptive Temperature-Based Replica Exchange for Sampling Large-Scale Protein Conformational Transitions. , 2013, Journal of chemical theory and computation.

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

[98]  Albert C. Pan,et al.  Molecular determinants of drug-receptor binding kinetics. , 2013, Drug discovery today.

[99]  Francesca Deflorian,et al.  High end GPCR design: crafted ligand design and druggability analysis using protein structure, lipophilic hotspots and explicit water networks , 2013, In Silico Pharmacology.

[100]  M. Parrinello,et al.  From metadynamics to dynamics. , 2013, Physical review letters.

[101]  Lin Li,et al.  On the Dielectric “Constant” of Proteins: Smooth Dielectric Function for Macromolecular Modeling and Its Implementation in DelPhi , 2013, Journal of chemical theory and computation.

[102]  G. Dahl,et al.  Pharmacokinetics and the drug-target residence time concept. , 2013, Drug discovery today.

[103]  G. Whitesides,et al.  Water networks contribute to enthalpy/entropy compensation in protein-ligand binding. , 2013, Journal of the American Chemical Society.

[104]  Jonathan W. Essex,et al.  Water Network Perturbation in Ligand Binding: Adenosine A2A Antagonists as a Case Study , 2013, J. Chem. Inf. Model..

[105]  Heike Schönherr,et al.  Profound methyl effects in drug discovery and a call for new C-H methylation reactions. , 2013, Angewandte Chemie.

[106]  D. Shaw Millisecond-Long Molecular Dynamics Simulations of Proteins on a Special-Purpose Machine , 2013 .

[107]  U. Ryde A fundamental view of enthalpy–entropy compensation , 2014 .

[108]  R. Altman,et al.  Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways , 2013, Nature chemistry.

[109]  E. Di Cera,et al.  Essential role of conformational selection in ligand binding. , 2014, Biophysical chemistry.

[110]  Andrea Cavalli,et al.  Steered Molecular Dynamics Simulations for Studying Protein-Ligand Interaction in Cyclin-Dependent Kinase 5 , 2014, J. Chem. Inf. Model..

[111]  Shuo Gu,et al.  Quantitatively Characterizing the Ligand Binding Mechanisms of Choline Binding Protein Using Markov State Model Analysis , 2014, PLoS Comput. Biol..

[112]  Michele Parrinello,et al.  Assessing the Reliability of the Dynamics Reconstructed from Metadynamics. , 2014, Journal of chemical theory and computation.

[113]  Stefano Moro,et al.  Supervised Molecular Dynamics (SuMD) as a Helpful Tool To Depict GPCR-Ligand Recognition Pathway in a Nanosecond Time Scale , 2014, J. Chem. Inf. Model..

[114]  Frank Noé,et al.  Markov state models of biomolecular conformational dynamics. , 2014, Current opinion in structural biology.

[115]  Stefano Moro,et al.  Bridging Molecular Docking to Membrane Molecular Dynamics To Investigate GPCR-Ligand Recognition: The Human A2A Adenosine Receptor as a Key Study , 2014, J. Chem. Inf. Model..

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

[117]  J. Mccammon,et al.  Accelerated molecular dynamics simulations of ligand binding to a muscarinic G-protein-coupled receptor , 2015, Quarterly Reviews of Biophysics.

[118]  Jonathan S. Mason,et al.  Decoding the Role of Water Dynamics in Ligand-Protein Unbinding: CRF1R as a Test Case , 2015, J. Chem. Inf. Model..

[119]  F. Totzke,et al.  Optimization of potent DFG-in inhibitors of platelet derived growth factor receptorβ (PDGF-Rβ) guided by water thermodynamics. , 2015, Journal of medicinal chemistry.

[120]  S. Moro,et al.  Understanding allosteric interactions in G protein-coupled receptors using Supervised Molecular Dynamics: A prototype study analysing the human A3 adenosine receptor positive allosteric modulator LUF6000. , 2015, Bioorganic & medicinal chemistry.

[121]  Luca Mollica,et al.  Kinetics of protein-ligand unbinding via smoothed potential molecular dynamics simulations , 2015, Scientific Reports.

[122]  A. Cavalli,et al.  The ligand binding mechanism to purine nucleoside phosphorylase elucidated via molecular dynamics and machine learning , 2015, Nature Communications.

[123]  F. Noé,et al.  Protein conformational plasticity and complex ligand-binding kinetics explored by atomistic simulations and Markov models , 2015, Nature Communications.

[124]  J. Ulander,et al.  Ligand Binding Thermodynamics in Drug Discovery: Still a Hot Tip? , 2015, Journal of medicinal chemistry.

[125]  A. Minarini,et al.  An optimized polyamine moiety boosts the potency of human type II topoisomerase poisons as quantified by comparative analysis centered on the clinical candidate F14512. , 2015, Chemical communications.

[126]  M. Parrinello,et al.  Kinetics of protein–ligand unbinding: Predicting pathways, rates, and rate-limiting steps , 2015, Proceedings of the National Academy of Sciences.

[127]  R. Leurs,et al.  Ligand Residence Time at G-protein–Coupled Receptors—Why We Should Take Our Time To Study It , 2015, Molecular Pharmacology.

[128]  S. Moro,et al.  Exploring the recognition pathway at the human A2A adenosine receptor of the endogenous agonist adenosine using supervised molecular dynamics simulations , 2015 .

[129]  Jürgen Bajorath,et al.  Predicting bioactive conformations and binding modes of macrocycles , 2016, Journal of Computer-Aided Molecular Design.

[130]  L. Pardo,et al.  The pathway of ligand entry from the membrane bilayer to a lipid G protein-coupled receptor , 2016, Scientific Reports.

[131]  Albert C. Pan,et al.  Molecular Basis of Ligand Dissociation from the Adenosine A2A Receptor , 2016, Molecular Pharmacology.

[132]  Michael S Bodnarchuk,et al.  Water, water, everywhere… It's time to stop and think. , 2016, Drug discovery today.

[133]  B. Berne,et al.  Spectral gap optimization of order parameters for sampling complex molecular systems , 2015, Proceedings of the National Academy of Sciences.

[134]  I. Dale,et al.  Can residence time offer a useful strategy to target agonist drugs for sustained GPCR responses? , 2016, Drug discovery today.

[135]  J. Izaguirre,et al.  Accelerated molecular dynamics simulations of the octopamine receptor using GPUs: discovery of an alternate agonist‐binding position , 2016, Proteins.

[136]  Hakan Gunaydin,et al.  Stacking with No Planarity? , 2016, ACS medicinal chemistry letters.

[137]  C. Tautermann Impact, determination and prediction of drug-receptor residence times for GPCRs. , 2016, Current opinion in pharmacology.

[138]  A. Cavalli,et al.  Molecular Dynamics Simulations and Kinetic Measurements to Estimate and Predict Protein-Ligand Residence Times. , 2016, Journal of medicinal chemistry.

[139]  B. Berne,et al.  How wet should be the reaction coordinate for ligand unbinding? , 2016, The Journal of chemical physics.

[140]  R. Copeland The drug–target residence time model: a 10-year retrospective , 2015, Nature Reviews Drug Discovery.

[141]  Stefano Moro,et al.  Deciphering the Complexity of Ligand-Protein Recognition Pathways Using Supervised Molecular Dynamics (SuMD) Simulations , 2016, J. Chem. Inf. Model..

[142]  Nargis Tabassum,et al.  Investigation of α-conotoxin unbinding using umbrella sampling. , 2016, Bioorganic & medicinal chemistry letters.

[143]  Huanxiang Liu,et al.  Computational study on the binding and unbinding mechanism of HCV NS5B with the inhibitor GS-461203 and substrate using conventional and steered molecular dynamics simulations , 2016 .