Water molecules in protein–ligand interfaces. Evaluation of software tools and SAR comparison

Targeting the interaction with or displacement of the ‘right’ water molecule can significantly increase inhibitor potency in structure-guided drug design. Multiple computational approaches exist to predict which waters should be targeted for displacement to achieve the largest gain in potency. However, the relative success of different methods remains underexplored. Here, we present a comparison of the ability of five water prediction programs (3D-RISM, SZMAP, WaterFLAP, WaterRank, and WaterMap) to predict crystallographic water locations, calculate their binding free energies, and to relate differences in these energies to observed changes in potency. The structural cohort included nine Bruton’s Tyrosine Kinase (BTK) structures, and nine bromodomain structures. Each program accurately predicted the locations of most crystallographic water molecules. However, the predicted binding free energies correlated poorly with the observed changes in inhibitor potency when solvent atoms were displaced by chemical changes in closely related compounds.Graphical abstract

[1]  T. Lazaridis.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 1. Theory , 1998 .

[2]  Alexander S. Bayden,et al.  The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery. , 2017, Journal of medicinal chemistry.

[3]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

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

[5]  D. J. Adams,et al.  Grand canonical ensemble Monte Carlo for a Lennard-Jones fluid , 1975 .

[6]  Michael K. Gilson,et al.  Erratum: “Grid inhomogeneous solvation theory: Hydration structure and thermodynamics of the miniature receptor cucurbit[7]uril” [J. Chem. Phys. 137, 044101 (2012)] , 2012 .

[7]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[8]  Demetri T. Moustakas,et al.  Evaluating Free Energies of Binding and Conservation of Crystallographic Waters Using SZMAP , 2015, J. Chem. Inf. Model..

[9]  W R Pitt,et al.  Modelling of solvent positions around polar groups in proteins. , 1991, Protein engineering.

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

[11]  M. Rarey,et al.  SIENA: Efficient Compilation of Selective Protein Binding Site Ensembles , 2016, J. Chem. Inf. Model..

[12]  P. Willett,et al.  SuperStar: improved knowledge-based interaction fields for protein binding sites. , 2001, Journal of molecular biology.

[13]  Matthias Rarey,et al.  From cheminformatics to structure-based design: Web services and desktop applications based on the NAOMI library. , 2017, Journal of biotechnology.

[14]  Gabriele Cruciani,et al.  A Common Reference Framework for Analyzing/Comparing Proteins and Ligands. Fingerprints for Ligands And Proteins (FLAP): Theory and Application , 2007, J. Chem. Inf. Model..

[15]  P. Sandy,et al.  Identification of a Benzoisoxazoloazepine Inhibitor (CPI-0610) of the Bromodomain and Extra-Terminal (BET) Family as a Candidate for Human Clinical Trials. , 2016, Journal of medicinal chemistry.

[16]  Adam R. Johnson,et al.  Battling Btk Mutants With Noncovalent Inhibitors That Overcome Cys481 and Thr474 Mutations. , 2016, ACS chemical biology.

[17]  Guanglei Cui,et al.  A Perspective on Water Site Prediction Methods for Structure Based Drug Design. , 2017, Current topics in medicinal chemistry.

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

[19]  K. Kosik,et al.  Structure-activity relationship study of 2,4-diaminothiazoles as Cdk5/p25 kinase inhibitors. , 2011, Bioorganic & medicinal chemistry letters.

[20]  S. Hymowitz,et al.  Potent and selective Bruton's tyrosine kinase inhibitors: discovery of GDC-0834. , 2015, Bioorganic & medicinal chemistry letters.

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

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

[23]  Brandon J. Bravo,et al.  Specific Btk inhibition suppresses B cell- and myeloid cell-mediated arthritis. , 2011, Nature chemical biology.

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

[25]  R. Knegtel,et al.  A Role for Hydration in Interleukin‐2 Inducible T Cell Kinase (Itk) Selectivity , 2011, Molecular informatics.

[26]  Robin Taylor,et al.  SuperStar: a knowledge-based approach for identifying interaction sites in proteins. , 1999, Journal of molecular biology.

[27]  James R Kiefer,et al.  Structural and thermodynamic characterization of the TYK2 and JAK3 kinase domains in complex with CP-690550 and CMP-6. , 2010, Journal of molecular biology.

[28]  Matthias Rarey,et al.  ASCONA: Rapid Detection and Alignment of Protein Binding Site Conformations , 2015, J. Chem. Inf. Model..

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

[30]  Heather A. Carlson,et al.  Predicting Displaceable Water Sites Using Mixed-Solvent Molecular Dynamics , 2018, J. Chem. Inf. Model..

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

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

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

[34]  Adam R. Johnson,et al.  Discovery of highly potent and selective Bruton's tyrosine kinase inhibitors: Pyridazinone analogs with improved metabolic stability. , 2016, Bioorganic & medicinal chemistry letters.

[35]  William R. Pitt,et al.  AQUARIUS2: Knowledge‐based modeling of solvent sites around proteins , 1993, J. Comput. Chem..

[36]  Andrew C Good,et al.  Diving into the Water: Inducible Binding Conformations for BRD4, TAF1(2), BRD9, and CECR2 Bromodomains. , 2016, Journal of medicinal chemistry.

[37]  Themis Lazaridis,et al.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 2. Applications to Simple Fluids , 1998 .

[38]  Michael K. Gilson,et al.  Thermodynamics of Water in an Enzyme Active Site: Grid-Based Hydration Analysis of Coagulation Factor Xa , 2014, Journal of chemical theory and computation.

[39]  Adam R. Johnson,et al.  Discovery of GDC-0853: A Potent, Selective, and Noncovalent Bruton's Tyrosine Kinase Inhibitor in Early Clinical Development. , 2018, Journal of medicinal chemistry.

[40]  Adam R. Johnson,et al.  Discovery of Potent and Selective Tricyclic Inhibitors of Bruton’s Tyrosine Kinase with Improved Druglike Properties , 2017, ACS medicinal chemistry letters.

[41]  Michael K Gilson,et al.  Grid inhomogeneous solvation theory: hydration structure and thermodynamics of the miniature receptor cucurbit[7]uril. , 2012, The Journal of chemical physics.

[42]  Jonathan W. Essex,et al.  Strategies to Calculate Water Binding Free Energies in Protein-Ligand Complexes , 2014, J. Chem. Inf. Model..

[43]  Denis Bucher,et al.  Shedding Light on Important Waters for Drug Design: Simulations versus Grid-Based Methods , 2018, J. Chem. Inf. Model..

[44]  Richard A. Friesner,et al.  Docking performance of the glide program as evaluated on the Astex and DUD datasets: a complete set of glide SP results and selected results for a new scoring function integrating WaterMap and glide , 2012, Journal of Computer-Aided Molecular Design.

[45]  Matthias Rarey,et al.  Estimating Electron Density Support for Individual Atoms and Molecular Fragments in X-ray Structures , 2017, J. Chem. Inf. Model..