Differential Water Thermodynamics Determine PI3K-Beta/Delta Selectivity for Solvent-Exposed Ligand Modifications

Phosphoinositide 3-kinases (PI3Ks) are involved in important cellular functions and represent desirable targets for drug discovery efforts, especially related to oncology; however, the four PI3K subtypes (α, β, γ, and δ) have highly similar binding sites, making the design of selective inhibitors challenging. A series of inhibitors with selectivity toward the β subtype over δ resulted in compound 3(S), which has entered a phase I/Ib clinical trial for patients with advanced PTEN-deficient cancer. Interestingly, X-ray crystallography revealed that the modifications making inhibitor 3(S) and related compounds selective toward the β-isoform do not interact directly with either PI3Kβ or PI3Kδ, thereby confounding rationalization of the SAR. Here, we apply explicit solvent molecular dynamics and solvent thermodynamic analysis using WaterMap in an effort to understand the unusual affinity and selectivity trends. We find that differences in solvent energetics and water networks, which are modulated upon binding of different ligands, explain the experimental affinity and selectivity trends. This study highlights the critical role of water molecules in molecular recognition and the importance of considering water networks in drug discovery efforts to rationalize and improve selectivity.

[1]  J. Ptak,et al.  High Frequency of Mutations of the PIK3CA Gene in Human Cancers , 2004, Science.

[2]  Robert A Pearlstein,et al.  Estimation of Solvation Entropy and Enthalpy via Analysis of Water Oxygen-Hydrogen Correlations. , 2015, Journal of chemical theory and computation.

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

[4]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[5]  Alice T. Loo,et al.  PTEN-deficient cancers depend on PIK3CB , 2008, Proceedings of the National Academy of Sciences.

[6]  J. F. Mack,et al.  Rational Design, Synthesis, and SAR of a Novel Thiazolopyrimidinone Series of Selective PI3K-beta Inhibitors. , 2012, ACS medicinal chemistry letters.

[7]  M. Bodkin,et al.  Evaluation of water displacement energetics in protein binding sites with grid cell theory. , 2015, Physical chemistry chemical physics : PCCP.

[8]  Lei Shi,et al.  Molecular determinants of selectivity and efficacy at the dopamine D3 receptor. , 2012, Journal of medicinal chemistry.

[9]  Fumio Hirata,et al.  Solvation thermodynamics of protein studied by the 3D-RISM theory , 2004 .

[10]  Woody Sherman,et al.  Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments , 2013, Journal of Computer-Aided Molecular Design.

[11]  Shaoyong Lu,et al.  Crystal Structures of PI3Kα Complexed with PI103 and Its Derivatives: New Directions for Inhibitors Design. , 2014, ACS medicinal chemistry letters.

[12]  Mike C. Payne,et al.  Thermodynamic Properties of Water Molecules at a Protein–Protein Interaction Surface , 2011, Journal of chemical theory and computation.

[13]  Stefano Forli,et al.  A force field with discrete displaceable waters and desolvation entropy for hydrated ligand docking. , 2012, Journal of medicinal chemistry.

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

[15]  É. Hajós‐Korcsok,et al.  Spirocyclic sulfamides as β-secretase 1 (BACE-1) inhibitors for the treatment of Alzheimer's disease: utilization of structure based drug design, WaterMap, and CNS penetration studies to identify centrally efficacious inhibitors. , 2012, Journal of medicinal chemistry.

[16]  W. L. Jorgensen,et al.  Energetics of displacing water molecules from protein binding sites: consequences for ligand optimization. , 2009, Journal of the American Chemical Society.

[17]  Jean-Pierre Marquette,et al.  Discovery and optimization of pyrimidone indoline amide PI3Kβ inhibitors for the treatment of phosphatase and tensin homologue (PTEN)-deficient cancers. , 2014, Journal of medicinal chemistry.

[18]  R. Farid,et al.  The translocation kinetics of antibiotics through porin OmpC: Insights from structure‐based solvation mapping using WaterMap , 2013, Proteins.

[19]  I. Braña,et al.  Clinical development of phosphatidylinositol 3-kinase inhibitors for cancer treatment , 2012, BMC Medicine.

[20]  Paul Labute,et al.  A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies , 2014, Journal of chemical theory and computation.

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

[22]  P. Norman.  Selective PI3Kδ inhibitors, a review of the patent literature , 2011, Expert opinion on therapeutic patents.

[23]  G. Whitesides,et al.  Interactions between Hofmeister anions and the binding pocket of a protein. , 2015, Journal of the American Chemical Society.

[24]  Markus A. Lill,et al.  WATsite: Hydration site prediction program with PyMOL interface , 2014, J. Comput. Chem..

[25]  Maxim V Fedorov,et al.  Communication: Accurate hydration free energies at a wide range of temperatures from 3D-RISM. , 2015, The Journal of chemical physics.

[26]  Daniel Cappel,et al.  Novel Inverse Binding Mode of Indirubin Derivatives Yields Improved Selectivity for DYRK Kinases , 2012, ACS medicinal chemistry letters.

[27]  J. Nicolas,et al.  Discovery and optimization of new benzimidazole- and benzoxazole-pyrimidone selective PI3Kβ inhibitors for the treatment of phosphatase and TENsin homologue (PTEN)-deficient cancers. , 2012, Journal of medicinal chemistry.

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

[29]  G. Keating.  Idelalisib: a review of its use in chronic lymphocytic leukaemia and indolent non-Hodgkin’s lymphoma , 2015, Targeted Oncology.

[30]  Chao Xu,et al.  Chemical basis for the recognition of trimethyllysine by epigenetic reader proteins , 2015, Nature Communications.

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

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

[33]  J. F. Mack,et al.  Synthesis and structure-activity relationships of imidazo[1,2-a]pyrimidin-5(1H)-ones as a novel series of beta isoform selective phosphatidylinositol 3-kinase inhibitors. , 2012, Bioorganic & medicinal chemistry letters.

[34]  J. Nicolas,et al.  Preparation and optimization of new 4-(morpholin-4-yl)-(6-oxo-1,6-dihydropyrimidin-2-yl)amide derivatives as PI3Kβ inhibitors. , 2012, Bioorganic & medicinal chemistry letters.

[35]  Jinhe Kim,et al.  Discovery of new aminopyrimidine-based phosphoinositide 3-kinase beta (PI3Kβ) inhibitors with selectivity over PI3Kα. , 2011, Bioorganic & medicinal chemistry letters.

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

[37]  W. Sherman,et al.  Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. , 2010, Journal of chemical theory and computation.

[38]  Ariel Fernández,et al.  Dehydron: a structurally encoded signal for protein interaction. , 2003, Biophysical journal.

[39]  Ronald M. Levy,et al.  PrimeX and the Schrödinger computational chemistry suite of programs , 2012 .

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

[41]  G. Whitesides,et al.  The binding of benzoarylsulfonamide ligands to human carbonic anhydrase is insensitive to formal fluorination of the ligand. , 2013, Angewandte Chemie.

[42]  Robert Abel,et al.  Motifs for molecular recognition exploiting hydrophobic enclosure in protein–ligand binding , 2007, Proceedings of the National Academy of Sciences.

[43]  C. Higgs,et al.  Hydration Site Thermodynamics Explain SARs for Triazolylpurines Analogues Binding to the A2A Receptor. , 2010, ACS medicinal chemistry letters.

[44]  William A Weiss,et al.  A dual PI3 kinase/mTOR inhibitor reveals emergent efficacy in glioma. , 2006, Cancer cell.

[45]  I. Adcock,et al.  Therapeutic Potential of Phosphatidylinositol 3-Kinase Inhibitors in Inflammatory Respiratory Disease , 2007, Journal of Pharmacology and Experimental Therapeutics.

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

[47]  C. García-echeverría,et al.  Class IA phosphatidylinositol 3-kinase: from their biologic implication in human cancers to drug discovery , 2008, Expert opinion on therapeutic targets.

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

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

[50]  Woody Sherman,et al.  Contribution of Explicit Solvent Effects to the Binding Affinity of Small‐Molecule Inhibitors in Blood Coagulation Factor Serine Proteases , 2011, ChemMedChem.

[51]  J. Luengo,et al.  [3a,4]-Dihydropyrazolo[1,5a]pyrimidines: Novel, Potent, and Selective Phosphatidylinositol-3-kinase β Inhibitors. , 2013, ACS medicinal chemistry letters.

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

[53]  W. Sherman,et al.  Understanding Kinase Selectivity Through Energetic Analysis of Binding Site Waters , 2010, ChemMedChem.

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

[55]  Woody Sherman,et al.  New hypotheses about the structure–function of proprotein convertase subtilisin/kexin type 9: Analysis of the epidermal growth factor‐like repeat A docking site using WaterMap , 2010, Proteins.

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

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

[58]  J. Luengo,et al.  Synthesis and structure-activity relationships of 1,2,4-triazolo[1,5-a]pyrimidin-7(3H)-ones as novel series of potent β isoform selective phosphatidylinositol 3-kinase inhibitors. , 2012, Bioorganic & medicinal chemistry letters.

[59]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

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