Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations

Water is the natural medium of molecules in the cell and plays an important role in protein structure, function and interaction with small molecule ligands. However, the widely used molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method for binding energy calculation does not explicitly take account of water molecules that mediate key protein-ligand interactions. We have developed a protocol to include water molecules that mediate ligand-protein interactions as part of the protein structure in calculation of MM/PBSA binding energies (a method we refer to as water-MM/PBSA) for a series of JNK3 kinase inhibitors. Improved correlation between water-MM/PBSA binding energies and experimental IC50 values was obtained compared to that obtained from classical MM/PBSA binding energy. This improved correlation was further validated using sets of neuraminidase and avidin inhibitors. The observed improvement, however, appears to be limited to systems in which there are water-mediated ligand-protein hydrogen bond interactions. We conclude that the water-MM/PBSA method performs better than classical MM/PBSA in predicting binding affinities when water molecules play a direct role in mediating ligand-protein hydrogen bond interactions.

[1]  Paul Beroza,et al.  Improved Ligand Binding Energies Derived from Molecular Dynamics: Replicate Sampling Enhances the Search of Conformational Space , 2013, J. Chem. Inf. Model..

[2]  S Cusack,et al.  Influenza B virus neuraminidase can synthesize its own inhibitor. , 1993, Structure.

[3]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[4]  R. Neitz,et al.  Design and synthesis of a novel, orally active, brain penetrant, tri-substituted thiophene based JNK inhibitor. , 2011, Bioorganic & medicinal chemistry letters.

[5]  Ray Luo,et al.  Assessment of linear finite‐difference Poisson–Boltzmann solvers , 2010, J. Comput. Chem..

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

[7]  S. Hannongbua,et al.  Bridge water mediates nevirapine binding to wild type and Y181C HIV-1 reverse transcriptase--evidence from molecular dynamics simulations and MM-PBSA calculations. , 2009, Journal of molecular graphics & modelling.

[8]  L. Kuhn,et al.  Cluster analysis of consensus water sites in thrombin and trypsin shows conservation between serine proteases and contributions to ligand specificity , 1998, Protein science : a publication of the Protein Society.

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

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

[11]  Samuel Genheden,et al.  Comparison of the Efficiency of the LIE and MM/GBSA Methods to Calculate Ligand-Binding Energies. , 2011, Journal of chemical theory and computation.

[12]  Brian W Matthews,et al.  A review about nothing: Are apolar cavities in proteins really empty? , 2009, Protein science : a publication of the Protein Society.

[13]  R. Kurumbail,et al.  Discovery of N-substituted pyridinones as potent and selective inhibitors of p38 kinase. , 2009, Bioorganic & medicinal chemistry letters.

[14]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[15]  Paul W Smith,et al.  Drug design against a shifting target: a structural basis for resistance to inhibitors in a variant of influenza virus neuraminidase. , 1998, Structure.

[16]  A. Ferrari,et al.  Validation of an automated procedure for the prediction of relative free energies of binding on a set of aldose reductase inhibitors. , 2007, Bioorganic & medicinal chemistry.

[17]  J. Andrew McCammon,et al.  MM-PBSA Captures Key Role of Intercalating Water Molecules at a Protein−Protein Interface , 2009, Journal of chemical theory and computation.

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

[19]  Julien Michel,et al.  Effects of Water Placement on Predictions of Binding Affinities for p38α MAP Kinase Inhibitors. , 2010, Journal of chemical theory and computation.

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