Advances in binding free energies calculations: QM/MM-based free energy perturbation method for drug design.

Multiple approaches have been devised and evaluated to computationally estimate binding free energies. Results using a recently developed Quantum Mechanics (QM)/Molecular Mechanics (MM) based Free Energy Perturbation (FEP) method suggest that this method has the potential to provide the most accurate estimation of binding affinities to date. The method treats ligands/inhibitors using QM while using MM for the rest of the system. The method has been applied and validated for a structurally diverse set of fructose 1,6- bisphosphatase (FBPase) inhibitors suggesting that the approach has the potential to be used as an integral part of drug discovery for both lead identification lead optimization, where there is a structure available. In addition, this QM/MM-based FEP method was shown to accurately replicate the anomalous hydration behavior exhibited by simple amines and amides suggesting that the method may also prove useful in predicting physical properties of molecules. While the method is about 5-fold more computationally demanding than conventional FEP, it has the potential to be less demanding on the end user since it avoids development of MM force field parameters for novel ligands and thereby eliminates this time-consuming step that often contributes significantly to the inaccuracy of binding affinity predictions using conventional FEP methods. The QM/MM-based FEP method has been extensively tested with respect to important considerations such as the length of the simulation required to obtain satisfactory convergence in the calculated relative solvation and binding free energies for both small and large structural changes between ligands. Future automation of the method and parallelization of the code is expected to enhance the speed and increase its use for drug design and lead optimization.

[1]  U. Singh,et al.  Use of a QM/MM-based FEP method to evaluate the anomalous hydration behavior of simple alkyl amines and amides: application to the design of FBPase inhibitors for the treatment of type-2 diabetes. , 2011, Journal of the American Chemical Society.

[2]  D. Granner,et al.  The genes of hepatic glucose metabolism. , 1990, The Journal of biological chemistry.

[3]  Santosh A. Khedkar,et al.  Pharmacophore modeling in drug discovery and development: an overview. , 2007, Medicinal chemistry (Shariqah (United Arab Emirates)).

[4]  Yanli Wang,et al.  Structure-Based Virtual Screening for Drug Discovery: a Problem-Centric Review , 2012, The AAPS Journal.

[5]  Woody Sherman,et al.  Computer-Aided Drug Design of Falcipain Inhibitors: Virtual Screening, Structure-Activity Relationships, Hydration Site Thermodynamics, and Reactivity Analysis , 2012, J. Chem. Inf. Model..

[6]  M. Rami Reddy,et al.  Free energy calculations in rational drug design , 2001 .

[7]  Jan Hermans,et al.  Reversibility of Free Energy Simulations∶ Slow Growth May Have a Unique Advantage. (With a Note on Use of Ewald Summation) , 2002 .

[8]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[9]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[10]  Abby L. Parrill,et al.  Rational drug design : novel methodology and practical applications , 1999 .

[11]  E. Arnett,et al.  Thermodynamics of Ionization and Solution of Aliphatic Amines in Water , 2007 .

[12]  I. Weber,et al.  Molecular mechanics analysis of inhibitor binding to HIV-1 protease. , 1992, Protein engineering.

[13]  William L. Jorgensen,et al.  Relative partition coefficients for organic solutes from fluid simulations , 1990 .

[14]  Brian K. Shoichet,et al.  Virtual Screening in Drug Discovery , 2005 .

[15]  Sheng-Yong Yang,et al.  Pharmacophore modeling and applications in drug discovery: challenges and recent advances. , 2010, Drug discovery today.

[16]  R. Wolfenden,et al.  Influences of solvent water on protein folding: free energies of solvation of cis and trans peptides are nearly identical. , 1988, Biochemistry.

[17]  Christopher R. Corbeil,et al.  Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go , 2008, British journal of pharmacology.

[18]  M. Rami Reddy,et al.  Relative solvation free energies calculated using an ab initio QM/MM-based free energy perturbation method: dependence of results on simulation length , 2009, J. Comput. Aided Mol. Des..

[19]  J N Weinstein,et al.  Relative differences in the binding free energies of human immunodeficiency virus 1 protease inhibitors: a thermodynamic cycle-perturbation approach. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[21]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[22]  Kenneth M. Merz,et al.  Drug Design : Structure-and Ligand-Based Approaches , 2017 .

[23]  W. V. van Gunsteren,et al.  Amine hydration: a united-atom force-field solution. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[24]  Nagamani Sukumar,et al.  Current trends in virtual high throughput screening using ligand-based and structure-based methods. , 2011, Combinatorial chemistry & high throughput screening.

[25]  Tomasz Arodz,et al.  Computational methods in developing quantitative structure-activity relationships (QSAR): a review. , 2006, Combinatorial chemistry & high throughput screening.

[26]  M. Erion,et al.  Relative binding affinities of fructose-1,6-bisphosphatase inhibitors calculated using a quantum mechanics-based free energy perturbation method. , 2007, Journal of the American Chemical Society.

[27]  P Reddanna,et al.  Free energy calculations to estimate ligand-binding affinities in structure-based drug design. , 2014, Current pharmaceutical design.

[28]  W. Guida,et al.  Structure-based design of inhibitors of purine nucleoside phosphorylase. 1. 9-(arylmethyl) derivatives of 9-deazaguanine. , 1993, Journal of medicinal chemistry.

[29]  Jitender Verma,et al.  3D-QSAR in drug design--a review. , 2010, Current topics in medicinal chemistry.

[30]  W. L. Jorgensen,et al.  A priori calculations of pKa's for organic compounds in water. The pKa of ethane , 1987 .

[31]  S. Shaikh,et al.  From drug target to leads--sketching a physicochemical pathway for lead molecule design in silico. , 2007, Current pharmaceutical design.

[32]  M. R. Reddy,et al.  Use of quantum mechanics/molecular mechanics-based FEP method for calculating relative binding affinities of FBPase inhibitors for type-2 diabetes , 2012, Theoretical Chemistry Accounts.

[33]  Atul Agarwal,et al.  Free Energy Calculations: Use and Limitations in Predicting Ligand Binding Affinities , 2007 .

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

[35]  B. Honig,et al.  New Model for Calculation of Solvation Free Energies: Correction of Self-Consistent Reaction Field Continuum Dielectric Theory for Short-Range Hydrogen-Bonding Effects , 1996 .

[36]  J. Åqvist,et al.  Ligand binding affinities from MD simulations. , 2002, Accounts of chemical research.

[37]  R. Wolfenden,et al.  Interaction of the peptide bond with solvent water: a vapor phase analysis. , 1978, Biochemistry.

[38]  Matthew P. Repasky,et al.  Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. , 2006, Journal of medicinal chemistry.

[39]  J. Andrew McCammon,et al.  Free energy difference calculations by thermodynamic integration: Difficulties in obtaining a precise value , 1991 .

[40]  S E Ealick,et al.  Structure-based design of inhibitors of purine nucleoside phosphorylase. , 1995, Acta crystallographica. Section D, Biological crystallography.

[41]  M. Rami Reddy,et al.  Calculation of relative solvation free energy differences by thermodynamic perturbation method: Dependence of free energy results on simulation length , 1999, J. Comput. Chem..

[42]  M. Rami Reddy,et al.  Minimum MD simulation length required to achieve reliable results in free energy perturbation calculations: Case study of relative binding free energies of fructose‐1,6‐bisphosphatase inhibitors , 2011, J. Comput. Chem..

[43]  M. Rami Reddy,et al.  Ab initio quantum mechanics‐based free energy perturbation method for calculating relative solvation free energies , 2007, J. Comput. Chem..

[44]  U. Singh,et al.  Development of a quantum mechanics-based free-energy perturbation method: use in the calculation of relative solvation free energies. , 2004, Journal of the American Chemical Society.

[45]  William L. Jorgensen,et al.  OPLS ALL-ATOM MODEL FOR AMINES : RESOLUTION OF THE AMINE HYDRATION PROBLEM , 1999 .

[46]  W. Sherman,et al.  Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization , 2011, Proteins.

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

[48]  A Wlodawer,et al.  An approach to rapid estimation of relative binding affinities of enzyme inhibitors: application to peptidomimetic inhibitors of the human immunodeficiency virus type 1 protease. , 1996, Journal of medicinal chemistry.

[49]  Peter A. Kollman,et al.  The overlooked bond‐stretching contribution in free energy perturbation calculations , 1991 .

[50]  K. Dill,et al.  Binding of small-molecule ligands to proteins: "what you see" is not always "what you get". , 2009, Structure.

[51]  P. A. Bash,et al.  Calculation of the relative change in binding free energy of a protein-inhibitor complex. , 1987, Science.

[52]  M. Erion,et al.  Calculation of relative binding free energy differences for fructose 1,6-bisphosphatase inhibitors using the thermodynamic cycle perturbation approach. , 2001, Journal of the American Chemical Society.

[53]  Donald G. Truhlar,et al.  Continuum Solvation Models: Classical and Quantum Mechanical Implementations , 2007 .

[54]  M. Erion,et al.  Calculation of Relative Hydration Free Energy Differences for Heteroaromatic Compounds: Use in the Design of Adenosine Deaminase and Cytidine Deaminase Inhibitors , 1998 .

[55]  Tao Jiang,et al.  MB06322 (CS-917): A potent and selective inhibitor of fructose 1,6-bisphosphatase for controlling gluconeogenesis in type 2 diabetes , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Alexander Golbraikh,et al.  Predictive QSAR modeling workflow, model applicability domains, and virtual screening. , 2007, Current pharmaceutical design.

[57]  W. Lipscomb,et al.  Discovery of potent and specific fructose-1,6-bisphosphatase inhibitors and a series of orally-bioavailable phosphoramidase-sensitive prodrugs for the treatment of type 2 diabetes. , 2007, Journal of the American Chemical Society.

[58]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[59]  W. Lipscomb,et al.  Structure-guided design of AMP mimics that inhibit fructose-1,6-bisphosphatase with high affinity and specificity. , 2007, Journal of the American Chemical Society.

[60]  Xiaoqin Zou,et al.  Inclusion of Solvation and Entropy in the Knowledge-Based Scoring Function for Protein-Ligand Interactions , 2010, J. Chem. Inf. Model..

[61]  Charles L. Brooks,et al.  λ‐Dynamics free energy simulation methods , 2009, J. Comput. Chem..

[62]  Arieh Warshel,et al.  On possible pitfalls in ab initio quantum mechanics/molecular mechanics minimization approaches for studies of enzymatic reactions. , 2005, The journal of physical chemistry. B.

[63]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[64]  T. Hansson,et al.  Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations. , 1995, Protein engineering.

[65]  Peter A. Kollman,et al.  Solvation Free Energies of Amides and Amines: Disagreement between Free Energy Calculations and Experiment , 1995 .

[66]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[67]  V. N. Viswanadhan,et al.  Rapid estimation of relative binding affinities of enzyme inhibitors , 2002 .

[68]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[69]  Arieh Warshel,et al.  Towards accurate ab initio QM/MM calculations of free-energy profiles of enzymatic reactions. , 2006, The journal of physical chemistry. B.

[70]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[71]  R. Friesner,et al.  Novel procedure for modeling ligand/receptor induced fit effects. , 2006, Journal of medicinal chemistry.

[72]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[73]  U. Singh,et al.  Hydrophobic hydration: A free energy perturbation study , 1989 .

[74]  H. Kubinyi,et al.  3D QSAR in drug design. , 2002 .