Kinetics, in silico docking, molecular dynamics, and MM‐GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP‐binding site inhibitors: The role of water molecules examined

With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC50 > 50 μM), indirubin‐3′‐oxime (IC50 = 144 nM), KT5720 (Ki = 18.4 nM) and staurosporine (Ki = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP‐binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MM‐GBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc‐staurosporine complex has the greatest number of receptor‐ligand hydrogen bonds, while for the indirubin‐3′‐oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor‐enzyme contacts. The MM‐GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1‐water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schrödinger induced‐fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor‐ligand contacts. Proteins 2011. © 2010 Wiley‐Liss, Inc.

[1]  Hilla Peretz,et al.  Ju n 20 03 Schrödinger ’ s Cat : The rules of engagement , 2003 .

[2]  J. Dunitz.  The entropic cost of bound water in crystals and biomolecules. , 1994, Science.

[3]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[4]  Richard H. Henchman,et al.  Entropic cost of protein-ligand binding and its dependence on the entropy in solution. , 2009, The journal of physical chemistry. B.

[5]  Mark S. Gordon,et al.  Self‐consistent molecular orbital methods. XXIII. A polarization‐type basis set for second‐row elements , 1982 .

[6]  E. Krebs,et al.  [49a] Muscle phosphorylase b: x Glucose-1-phosphate+ Gn ⇄ Gn+x + x inorganic phosphate(where Gn designates glycogen containing n glucose residues) , 1962 .

[7]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[8]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[9]  Mindy I. Davis,et al.  A quantitative analysis of kinase inhibitor selectivity , 2008, Nature Biotechnology.

[10]  Annick Dejaegere,et al.  Protein–protein recognition and interaction hot spots in an antigen–antibody complex: Free energy decomposition identifies “efficient amino acids” , 2007, Proteins.

[11]  Michael K. Gilson,et al.  Theory of Free Energy and Entropy in Noncovalent Binding , 2009 .

[12]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[13]  G L Snyder,et al.  Indirubins inhibit glycogen synthase kinase-3 beta and CDK5/p25, two protein kinases involved in abnormal tau phosphorylation in Alzheimer's disease. A property common to most cyclin-dependent kinase inhibitors? , 2001, The Journal of biological chemistry.

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

[15]  Christodoulos A Floudas,et al.  Species specificity of the complement inhibitor compstatin investigated by all‐atom molecular dynamics simulations , 2010, Proteins.

[16]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

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

[18]  M. M. Bradford.  A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. , 1976, Analytical biochemistry.

[19]  Claudio N. Cavasotto and Narender Singh.  Docking and High Throughput Docking: Successes and the Challenge of Protein Flexibility , 2008 .

[20]  L. Johnson,et al.  Two structures of the catalytic domain of phosphorylase kinase: an active protein kinase complexed with substrate analogue and product. , 1995, Structure.

[21]  P. Cohen,et al.  The selectivity of protein kinase inhibitors: a further update. , 2007, The Biochemical journal.

[22]  Themis Lazaridis,et al.  Contributions to the binding free energy of ligands to avidin and streptavidin , 2002, Proteins.

[23]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[24]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

[25]  Bing Zhang,et al.  Significance of Water Molecules in the Inhibition of Cylin-Dependent Kinase 2 and 5 Complexes , 2007, J. Chem. Inf. Model..

[26]  G Eisenbrand,et al.  Inhibitor binding to active and inactive CDK2: the crystal structure of CDK2-cyclin A/indirubin-5-sulphonate. , 2001, Structure.

[27]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[28]  Alexander D. MacKerell,et al.  Computational evaluation of protein-small molecule binding. , 2009, Current opinion in structural biology.

[29]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[30]  Mark A Olson,et al.  Calculation of absolute protein-ligand binding affinity using path and endpoint approaches. , 2006, Biophysical journal.

[31]  Hans Briem,et al.  From the Insoluble Dye Indirubin towards Highly Active, Soluble CDK2‐Inhibitors , 2005, Chembiochem : a European journal of chemical biology.

[32]  V. Tan,et al.  Study of the inhibition of cyclin-dependent kinases with roscovitine and indirubin-3′-oxime from molecular dynamics simulations , 2006, Journal of molecular modeling.

[33]  V. Zoete,et al.  MM–GBSA binding free energy decomposition and T cell receptor engineering , 2010, Journal of molecular recognition : JMR.

[34]  Li-Jun Yang,et al.  Binding energy contributions of the conserved bridging water molecules in CDK2-inhibitor complexes: A combined QM/MM study , 2008 .

[35]  R. Friesner,et al.  Novel procedure for modeling ligand/receptor induced fit effects. , 2006, Journal of medicinal chemistry.

[36]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

[37]  Arieh Warshel,et al.  Absolute binding free energy calculations: On the accuracy of computational scoring of protein–ligand interactions , 2010, Proteins.

[38]  B. Honig,et al.  A hierarchical approach to all‐atom protein loop prediction , 2004, Proteins.

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

[40]  Richard A. Lewis,et al.  Lessons in molecular recognition: the effects of ligand and protein flexibility on molecular docking accuracy. , 2004, Journal of medicinal chemistry.

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

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

[43]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[44]  E F Garman,et al.  Expression, purification and crystallisation of phosphorylase kinase catalytic domain. , 1995, Journal of molecular biology.

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

[46]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[47]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[48]  C. Cori,et al.  THE ROLE OF ADENYLIC ACID IN THE ACTIVATION OF PHOSPHORYLASE. , 1964, Proceedings of the National Academy of Sciences of the United States of America.

[49]  N. Oikonomakos,et al.  Recognition of ribonuclease A by 3'-5'-pyrophosphate-linked dinucleotide inhibitors: a molecular dynamics/continuum electrostatics analysis. , 2007, Biophysical journal.

[50]  M. Bollen,et al.  An optimized assay of phosphorylase kinase in crude liver preparations. , 1991, The Biochemical journal.

[51]  L. L. Kastenschmidt,et al.  Subunit interactions and their relationship to the allosteric properties of rabbit skeletal muscle phosphorylase b. , 1968, Biochemistry.

[52]  L. Johnson,et al.  The Structure of Phosphorylase Kinase Holoenzyme at 9.9 Å Resolution and Location of the Catalytic Subunit and the Substrate Glycogen Phosphorylase , 2009, Structure.

[53]  Giulio Rastelli,et al.  Fast and accurate predictions of binding free energies using MM‐PBSA and MM‐GBSA , 2009, J. Comput. Chem..

[54]  M. Noble,et al.  Structures of P. falciparum PfPK5 test the CDK regulation paradigm and suggest mechanisms of small molecule inhibition. , 2003, Structure.

[55]  S. Zographos,et al.  Crystallographic and computational studies on 4‐phenyl‐N‐(β‐D‐glucopyranosyl)‐1H‐1,2,3‐triazole‐1‐acetamide, an inhibitor of glycogen phosphorylase: Comparison with α‐D‐glucose, N‐acetyl‐β‐D‐glucopyranosylamine and N‐benzoyl‐N′‐β‐D‐glucopyranosyl urea binding , 2007 .

[56]  Philip M. Dean,et al.  Hydration in drug design. 1. Multiple hydrogen-bonding features of water molecules in mediating protein-ligand interactions , 1995, J. Comput. Aided Mol. Des..

[57]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[58]  Doris Marko,et al.  Indirubin, the active constituent of a Chinese antileukaemia medicine, inhibits cyclin-dependent kinases , 1999, Nature Cell Biology.

[59]  Liisa Holm,et al.  Searching protein structure databases with DaliLite v.3 , 2008, Bioinform..

[60]  Andrea Musacchio,et al.  Mechanism of CDK5/p25 binding by CDK inhibitors. , 2005, Journal of medicinal chemistry.

[61]  P. Cohen,et al.  Specificity and mechanism of action of some commonly used protein kinase inhibitors. , 2000, The Biochemical journal.

[62]  L. Johnson,et al.  Catalytic mechanism of phosphorylase kinase probed by mutational studies. , 1999, Biochemistry.

[63]  L. Johnson,et al.  The crystal structure of a phosphorylase kinase peptide substrate complex: kinase substrate recognition , 1997, The EMBO journal.

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

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

[66]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

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

[68]  P. Greengard,et al.  Structural basis for the synthesis of indirubins as potent and selective inhibitors of glycogen synthase kinase-3 and cyclin-dependent kinases. , 2004, Journal of medicinal chemistry.

[69]  R. DeFronzo,et al.  Pathogenesis of NIDDM: A Balanced Overview , 1992, Diabetes Care.

[70]  Y. Kai,et al.  THE SYNTHESIS, ANTILEUKEMIC ACTIVITY, AND CRYSTAL STRUCTURES OF INDIRUBIN DERIVATIVES , 1996 .

[71]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[72]  G. Chang,et al.  An internal-coordinate Monte Carlo method for searching conformational space , 1989 .

[73]  M Levitt,et al.  Water: now you see it, now you don't. , 1993, Structure.

[74]  M Karplus,et al.  Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase. , 2001, Journal of molecular biology.

[75]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[76]  Gerhard Eisenbrand,et al.  Binding of the potential antitumour agent indirubin-5-sulphonate at the inhibitor site of rabbit muscle glycogen phosphorylase b. Comparison with ligand binding to pCDK2-cyclin A complex. , 2004, European journal of biochemistry.

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

[78]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[79]  Irwin D Kuntz,et al.  Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods. , 2003, Biopolymers.

[80]  PatrickY.-S. Lam,et al.  Rational design of potent, bioavailable, nonpeptide cyclic ureas as HIV protease inhibitors. , 1994, Science.

[81]  Philip M. Dean,et al.  Hydration in drug design. 3. Conserved water molecules at the ligand-binding sites of homologous proteins , 1995, J. Comput. Aided Mol. Des..