Addressing Limitations with the MM-GB/SA Scoring Procedure using the WaterMap Method and Free Energy Perturbation Calculations

The MM-GB/SA scoring technique has become an important computational approach in drug design. We, and others, have demonstrated that for congeneric molecules the correlation with experimental data obtained with the physics-based scoring is usually superior to scoring functions from typical docking algorithms. Despite showing good accuracy when applied within a series, much work is necessary to improve the MM-GB/SA method in order to gain greater efficiency in drug design. Here, we investigate the poor estimation of protein desolvation provided by the GB/SA solvation model and the large dynamic range observed in the MM-GB/SA scoring compared to that of the experimental data. In the former, replacing the GB/SA protein desolvation in the MM-GB/SA method by the free energy associated with displacing binding site waters upon ligand binding estimated by WaterMap provides the best results when ranking congeneric series of factor Xa and cyclin-dependent kinase 2 (CDK2) inhibitors. However, the improvement is modest over results obtained with the MM-GB/SA and WaterMap methods individually, apparently due to the high correlation between the free energy liberation of the displaced solvent and the protein-ligand van der Waals interactions, which in turn may be interpretable as estimates of the hydrophobic effect and hydrophobic-like interactions, respectively. As for the large dynamic range, comparisons between MM-GB/SA and FEP calculations indicate that for the factor Xa test set this problem has its origin in the lack of shielding effects of protein--ligand electrostatic interactions; that overly favors ligands that engage in hydrogen bonds with the protein.

[1]  Yaxiong Sun,et al.  Improving Docking Accuracy through Molecular Mechanics Generalized Born Optimization and Scoring. , 2007, Journal of chemical theory and computation.

[2]  Thomas Simonson,et al.  Free Energy Simulations Come of Age: Protein—Ligand Recognition , 2002 .

[3]  Paul D Lyne,et al.  Accurate prediction of the relative potencies of members of a series of kinase inhibitors using molecular docking and MM-GBSA scoring. , 2006, Journal of medicinal chemistry.

[4]  Niu Huang,et al.  Physics-based methods for studying protein-ligand interactions. , 2007, Current opinion in drug discovery & development.

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

[6]  Jonathan W. Essex,et al.  An empirical boundary potential for water droplet simulations , 1995, J. Comput. Chem..

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

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

[9]  William L. Jorgensen,et al.  Improved convergence of binding affinities with free energy perturbation: Application to nonpeptide ligands with pp60src SH2 domain , 2001, J. Comput. Aided Mol. Des..

[10]  William L. Jorgensen,et al.  Free energy calculations: a breakthrough for modeling organic chemistry in solution , 1989 .

[11]  Matthew P Jacobson,et al.  Virtual Ligand Screening against Escherichia coli Dihydrofolate Reductase: Improving Docking Enrichment Using Physics-Based Methods , 2005, Journal of biomolecular screening.

[12]  Woody Sherman,et al.  High‐energy water sites determine peptide binding affinity and specificity of PDZ domains , 2009, Protein science : a publication of the Protein Society.

[13]  Michael H. Mazor,et al.  Hydration of cavities in proteins : a molecular dynamics approach , 1990 .

[14]  William L. Jorgensen,et al.  Free Energy Changes in Solution , 2002 .

[15]  Cristiano Ruch Werneck Guimarães,et al.  MM-GB/SA Rescoring of Docking Poses in Structure-Based Lead Optimization , 2008, J. Chem. Inf. Model..

[16]  M. Gilson,et al.  Ligand configurational entropy and protein binding , 2007, Proceedings of the National Academy of Sciences.

[17]  J. Andrew McCammon,et al.  Computation of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation , 1993 .

[18]  Cristiano Ruch Werneck Guimarães,et al.  Elucidation of fatty acid amide hydrolase inhibition by potent alpha-ketoheterocycle derivatives from Monte Carlo simulations. , 2005, Journal of the American Chemical Society.

[19]  P. Kollman,et al.  Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of molecular mechanics and continuum solvent models. , 2000, Journal of medicinal chemistry.

[20]  Hongming Wang,et al.  Virtual fragment screening: an exploration of various docking and scoring protocols for fragments using Glide , 2009, J. Comput. Aided Mol. Des..

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

[22]  J. A. McCammon,et al.  Dynamics and Design of Enzymes and Inhibitors. , 1986 .

[23]  Thomas Lengauer,et al.  A fast flexible docking method using an incremental construction algorithm. , 1996, Journal of molecular biology.

[24]  M. Jacobson,et al.  Molecular mechanics methods for predicting protein-ligand binding. , 2006, Physical chemistry chemical physics : PCCP.

[25]  William L. Jorgensen,et al.  Free Energies of Hydration from a Generalized Born Model and an All-Atom Force Field , 2004 .

[26]  Y. Martin,et al.  A general and fast scoring function for protein-ligand interactions: a simplified potential approach. , 1999, Journal of medicinal chemistry.

[27]  J Hermans,et al.  Hydrophilicity of cavities in proteins , 1996, Proteins.

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

[29]  Lauren Wickstrom,et al.  Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation. , 2007, The journal of physical chemistry. B.

[30]  R. Levy,et al.  Enthalpy−Entropy and Cavity Decomposition of Alkane Hydration Free Energies: Numerical Results and Implications for Theories of Hydrophobic Solvation , 2000 .

[31]  Peter Naur,et al.  The glutamate receptor GluR5 agonist (S)-2-amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]isoxazol-4-yl)propionic acid and the 8-methyl analogue: synthesis, molecular pharmacology, and biostructural characterization. , 2009, Journal of medicinal chemistry.

[32]  Julian Tirado-Rives,et al.  Contribution of conformer focusing to the uncertainty in predicting free energies for protein-ligand binding. , 2006, Journal of medicinal chemistry.

[33]  C. E. Peishoff,et al.  A critical assessment of docking programs and scoring functions. , 2006, Journal of medicinal chemistry.

[34]  Niu Huang,et al.  Physics-Based Scoring of Protein-Ligand Complexes: Enrichment of Known Inhibitors in Large-Scale Virtual Screening , 2006, J. Chem. Inf. Model..

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

[36]  J M Blaney,et al.  A geometric approach to macromolecule-ligand interactions. , 1982, Journal of molecular biology.

[37]  M. Orozco,et al.  How accurate can molecular dynamics/linear response and Poisson–Boltzmann/solvent accessible surface calculations be for predicting relative binding affinities? Acetylcholinesterase huprine inhibitors as a test case , 2001 .

[38]  A. Beatty,et al.  Design, synthesis, and activity of 2,6-diphenoxypyridine-derived factor Xa inhibitors. , 1999, Journal of medicinal chemistry.

[39]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[40]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[41]  Julien Michel,et al.  Prediction of the water content in protein binding sites. , 2009, The journal of physical chemistry. B.

[42]  Scott P. Brown,et al.  Healthy skepticism: assessing realistic model performance. , 2009, Drug discovery today.

[43]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[44]  Wei Xu,et al.  Design, synthesis, and activity of a novel series of factor Xa inhibitors: optimization of arylamidine groups. , 2002, Journal of medicinal chemistry.

[45]  D. Case,et al.  Generalized Born Models of Macromolecular Solvation Effects , 2001 .

[46]  Wilfred F. van Gunsteren,et al.  The importance of solute-solvent van der Waals interactions with interior atoms of biopolymers. , 2001 .

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

[48]  M. Murcko,et al.  Consensus scoring: A method for obtaining improved hit rates from docking databases of three-dimensional structures into proteins. , 1999, Journal of medicinal chemistry.

[49]  M Rarey,et al.  Detailed analysis of scoring functions for virtual screening. , 2001, Journal of medicinal chemistry.

[50]  P. A. Harris,et al.  Oxindole-based inhibitors of cyclin-dependent kinase 2 (CDK2): design, synthesis, enzymatic activities, and X-ray crystallographic analysis. , 2001, Journal of medicinal chemistry.

[51]  Caterina Barillari,et al.  Classification of water molecules in protein binding sites. , 2007, Journal of the American Chemical Society.

[52]  P. Charifson,et al.  Are free energy calculations useful in practice? A comparison with rapid scoring functions for the p38 MAP kinase protein system. , 2001, Journal of medicinal chemistry.

[53]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

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

[55]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[56]  N. Foloppe,et al.  Towards predictive ligand design with free-energy based computational methods? , 2006, Current medicinal chemistry.

[57]  David A Pearlman,et al.  Evaluating the molecular mechanics poisson-boltzmann surface area free energy method using a congeneric series of ligands to p38 MAP kinase. , 2005, Journal of medicinal chemistry.

[58]  Stefan Boresch,et al.  THE ROLE OF BONDED TERMS IN FREE ENERGY SIMULATIONS : 1. THEORETICAL ANALYSIS , 1999 .

[59]  W Patrick Walters,et al.  A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance , 2004, Proteins.

[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]  W F van Gunsteren,et al.  The importance of solute-solvent van der Waals interactions with interior atoms of biopolymers. , 2001, Journal of the American Chemical Society.