Evaluating Free Energies of Binding and Conservation of Crystallographic Waters Using SZMAP

The SZMAP method computes binding free energies and the corresponding thermodynamic components for water molecules in the binding site of a protein structure [ SZMAP, 1.0.0 ; OpenEye Scientific Software Inc. : Santa Fe, NM, USA , 2011 ]. In this work, the ability of SZMAP to predict water structure and thermodynamic stability is examined for the X-ray crystal structures of a series of protein-ligand complexes. SZMAP results correlate with higher-level replica exchange thermodynamic integration double decoupling calculations of the absolute free energy of bound waters in the test set complexes. In addition, SZMAP calculations show good agreement with experimental data in terms of water conservation (across multiple crystal structures) and B-factors over a subset of the test set. In particular, the SZMAP neutral entropy difference term calculated at crystallographic water positions within each of the complex structures correlates well with whether that crystallographic water is conserved or displaceable. Furthermore, the calculated entropy of the water probe relative to the continuum shows a significant degree of correlation with the B-factors associated with the oxygen atoms of the water molecules. Taken together, these results indicate that SZMAP is capable of quantitatively predicting water positions and their energetics and is potentially a useful tool for determining which waters to attempt to displace, maintain, or build in through water-mediated interactions when evolving a lead series during a drug discovery program.

[1]  Gerhard Klebe,et al.  Ligand binding stepwise disrupts water network in thrombin: enthalpic and entropic changes reveal classical hydrophobic effect. , 2012, Journal of medicinal chemistry.

[2]  Charles L Brooks,et al.  Recent advances in implicit solvent-based methods for biomolecular simulations. , 2008, Current opinion in structural biology.

[3]  Y. Lindqvist,et al.  High‐resolution structures of scytalone dehydratase‐inhibitor complexes crystallized at physiological pH , 1999, Proteins.

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

[5]  Matthias Rarey,et al.  Evidence of Water Molecules - A Statistical Evaluation of Water Molecules Based on Electron Density , 2015, J. Chem. Inf. Model..

[6]  J M Bland,et al.  Statistical methods for assessing agreement between two methods of clinical measurement , 1986 .

[7]  Nathan A. Baker,et al.  Improving implicit solvent simulations: a Poisson-centric view. , 2005, Current opinion in structural biology.

[8]  Klaus R. Liedl,et al.  A GRID-Derived Water Network Stabilizes Molecular Dynamics Computer Simulations of a Protease , 2011, J. Chem. Inf. Model..

[9]  Glen E Kellogg,et al.  The Importance of Being Exhaustive. Optimization of Bridging Structural Water Molecules and Water Networks in Models of Biological Systems , 2004, Chemistry & biodiversity.

[10]  L D Jennings,et al.  A new potent inhibitor of fungal melanin biosynthesis identified through combinatorial chemistry. , 1999, Bioorganic & medicinal chemistry letters.

[11]  George M. Whitesides,et al.  Mechanism of the hydrophobic effect in the biomolecular recognition of arylsulfonamides by carbonic anhydrase , 2011, Proceedings of the National Academy of Sciences.

[12]  Andrea Bortolato,et al.  New insights from structural biology into the druggability of G protein-coupled receptors. , 2012, Trends in pharmacological sciences.

[13]  Felice C Lightstone,et al.  Approaches to efficiently estimate solvation and explicit water energetics in ligand binding: the use of WaterMap , 2013, Expert opinion on drug discovery.

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

[15]  Chao-Yie Yang,et al.  Binding free energy contributions of interfacial waters in HIV-1 protease/inhibitor complexes. , 2006, Journal of the American Chemical Society.

[16]  J. Andrew Grant,et al.  A smooth permittivity function for Poisson–Boltzmann solvation methods , 2001, J. Comput. Chem..

[17]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[18]  F. Allen.  The Cambridge Structural Database: a quarter of a million crystal structures and rising. , 2002, Acta crystallographica. Section B, Structural science.

[19]  Fumio Hirata,et al.  Analysis of biomolecular solvation sites by 3D-RISM theory. , 2013, The journal of physical chemistry. B.

[20]  Z. Xiang,et al.  On the role of the crystal environment in determining protein side-chain conformations. , 2002, Journal of molecular biology.

[21]  G. Basarab,et al.  Structure-based design of potent inhibitors of scytalone dehydratase: displacement of a water molecule from the active site. , 1998, Biochemistry.

[22]  Thomas A. Halgren,et al.  Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular. interactions , 1996, J. Comput. Chem..

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

[24]  G. Wagner,et al.  Detection of long-lived bound water molecules in complexes of human dihydrofolate reductase with methotrexate and NADPH. , 1995, Journal of molecular biology.

[25]  P. Andrew Karplus,et al.  Ordered water in macromolecular structure , 1994 .

[26]  M Karplus,et al.  Estimation of uncertainties in X‐ray refinement results by use of perturbed structures , 1987, Proteins.

[27]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

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

[29]  Vassilios Myrianthopoulos,et al.  An inhibitor-driven study for enhancing the selectivity of indirubin derivatives towards leishmanial Glycogen Synthase Kinase-3 over leishmanial cdc2-related protein kinase 3 , 2013, Parasites & Vectors.

[30]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[31]  Ricardo L. Mancera,et al.  WaterScore: a novel method for distinguishing between bound and displaceable water molecules in the crystal structure of the binding site of protein-ligand complexes , 2003, Journal of molecular modeling.

[32]  Franca Fraternali,et al.  Design and application of implicit solvent models in biomolecular simulations , 2014, Current opinion in structural biology.

[33]  Pietro Cozzini,et al.  Robust classification of "relevant" water molecules in putative protein binding sites. , 2008, Journal of medicinal chemistry.

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

[35]  David J. Huggins,et al.  Benchmarking the thermodynamic analysis of water molecules around a model beta sheet , 2012, J. Comput. Chem..

[36]  David Chandler,et al.  Density functional theory of nonuniform polyatomic systems. I. General formulation , 1986 .

[37]  James M. Hogle,et al.  Functional group placement in protein binding sites: a comparison of GRID and MCSS , 2001, J. Comput. Aided Mol. Des..

[38]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

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

[40]  A Wlodawer,et al.  Inhibitors of HIV-1 protease: a major success of structure-assisted drug design. , 1998, Annual review of biophysics and biomolecular structure.

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

[42]  Robert Abel,et al.  Computational methods for high resolution prediction and refinement of protein structures. , 2013, Current opinion in structural biology.

[43]  Ariel Fernández,et al.  Epistructural tension promotes protein associations. , 2012, Physical review letters.

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

[45]  Glen E. Kellogg,et al.  Hydrophobicity: is LogPo/w more than the sum of its parts? , 2000 .

[46]  David J Huggins,et al.  Application of inhomogeneous fluid solvation theory to model the distribution and thermodynamics of water molecules around biomolecules. , 2012, Physical chemistry chemical physics : PCCP.

[47]  Joanna Trylska,et al.  Thermodynamic linkage between the binding of protons and inhibitors to HIV‐1 protease , 2008, Protein science : a publication of the Protein Society.

[48]  T. C. Bruice,et al.  Role of a critical water in scytalone dehydratase-catalyzed reaction. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

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

[50]  J. Trylska,et al.  Continuum molecular electrostatics, salt effects, and counterion binding—A review of the Poisson–Boltzmann theory and its modifications , 2008, Biopolymers.

[51]  Martin Smiesko,et al.  AcquaAlta: A Directional Approach to the Solvation of Ligand-Protein Complexes , 2011, J. Chem. Inf. Model..

[52]  Supot Hannongbua,et al.  Accurate prediction of protonation state as a prerequisite for reliable MM‐PB(GB)SA binding free energy calculations of HIV‐1 protease inhibitors , 2008, J. Comput. Chem..

[53]  C. Spearman.  The proof and measurement of association between two things. , 2015, International journal of epidemiology.

[54]  P. Labute.  proteins STRUCTURE O FUNCTION O BIOINFORMATICS Protonate3D: Assignment of ionization , 2013 .

[55]  Carlos Simmerling,et al.  Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber. , 2010, Journal of chemical theory and computation.

[56]  Kengo Kinoshita,et al.  Blind prediction of interfacial water positions in CAPRI , 2014, Proteins.

[57]  S. Parthasarathy,et al.  Analysis of temperature factor distribution in high‐resolution protein structures , 1997, Protein science : a publication of the Protein Society.

[58]  Pietro Cozzini,et al.  Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules. , 2006, Journal of molecular biology.

[59]  Charles J. Eyermann,et al.  NMR and X-ray Evidence That the HIV Protease Catalytic Aspartyl Groups Are Protonated in the Complex Formed by the Protease and a Non-Peptide Cyclic Urea-Based Inhibitor , 1994 .

[60]  Stuart L. Schreiber,et al.  Ligand design by a combinatorial approach based on modeling and experiment: application to HLA-DR4 , 2007, J. Comput. Aided Mol. Des..

[61]  W. Punch,et al.  Predicting conserved water-mediated and polar ligand interactions in proteins using a K-nearest-neighbors genetic algorithm. , 1997, Journal of molecular biology.

[62]  A. Fedorov,et al.  Comparison of experimental and computational functional group mapping of an RNase A structure: implications for computer-aided drug design. , 1996, Protein engineering.

[63]  Glen Eugene Kellogg,et al.  Web application for studying the free energy of binding and protonation states of protein–ligand complexes based on HINT , 2009, J. Comput. Aided Mol. Des..

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

[65]  Woody Sherman,et al.  Contributions of water transfer energy to protein‐ligand association and dissociation barriers: Watermap analysis of a series of p38α MAP kinase inhibitors , 2013, Proteins.

[66]  M. Karplus,et al.  Functionality maps of binding sites: A multiple copy simultaneous search method , 1991, Proteins.

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

[68]  Guanglei Cui,et al.  SPAM: A Simple Approach for Profiling Bound Water Molecules. , 2013, Journal of chemical theory and computation.