Clarifying binding difference of ATP and ADP to extracellular signal‐regulated kinase 2 by using molecular dynamics simulations

Extracellular signal‐regulated kinase 2 is a promising target for designs and development of anticancer drugs. Molecular dynamics simulations and molecular mechanics Poisson–Boltzmann method were applied to study binding difference of ADP and ATP to extracellular signal‐regulated kinase 2. The results prove that the binding ability of ATP to extracellular signal‐regulated kinase 2 is stronger than that of ADP. Principal component analysis performed by using molecular dynamics trajectories suggests that binding of ADP and ATP to extracellular signal‐regulated kinase 2 change motion directions of two helices α1 and α2. Residue‐based free energy decomposition method was adopted to calculate contributions of separate residues to associations of ADP and ATP with extracellular signal‐regulated kinase 2. The results show that ADP and ATP produce strong CH–π interactions with five residues Ile29, Val37, Ala50, Leu105, and Leu154. In addition, five hydrogen bonding interactions of ADP and ATP with residues Lys52, Gln103, Asp104, and Met106 also stabilize bindings of ADP and ATP to extracellular signal‐regulated kinase 2. Overall, the CH–π interactions of ATP with five residues Ile29, Val37, Ala50, Leu105, and Leu154 are stronger than ADP. This study is expected to contribute a significant theoretical hint for designs of anticancer drugs targeting extracellular signal‐regulated kinase 2.

[1]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[2]  F. Cunha,et al.  Peripheral analgesic blockade of hypernociception: activation of arginine/NO/cGMP/protein kinase G/ATP-sensitive K+ channel pathway. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[4]  Susan S. Taylor,et al.  Direct evidence that oncogenic tyrosine kinases and cyclic AMP-dependent protein kinase have homologous ATP-binding sites , 1984, Nature.

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

[6]  M. Meyerson,et al.  The T790M mutation in EGFR kinase causes drug resistance by increasing the affinity for ATP , 2008, Proceedings of the National Academy of Sciences.

[7]  Robert G. Weiss,et al.  Altered Creatine Kinase Adenosine Triphosphate Kinetics in Failing Hypertrophied Human Myocardium , 2006, Circulation.

[8]  Monique Laberge,et al.  Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics. , 2008, Biophysical journal.

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

[10]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[11]  P. Bharatam,et al.  Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity , 2015, Journal of biomolecular structure & dynamics.

[12]  Shaoyong Lu,et al.  Harnessing Allostery: A Novel Approach to Drug Discovery , 2014, Medicinal research reviews.

[13]  HuangDanzhi,et al.  Kinase selectivity potential for inhibitors targeting the ATP binding site , 2010 .

[14]  Ruibo Wu,et al.  Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study. , 2010, Journal of chemical theory and computation.

[15]  David L. Mobley,et al.  A Python tool to set up relative free energy calculations in GROMACS , 2015, Journal of Computer-Aided Molecular Design.

[16]  Bernard R Brooks,et al.  Calculations of Solvation Free Energy through Energy Reweighting from Molecular Mechanics to Quantum Mechanics. , 2016, Journal of chemical theory and computation.

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

[18]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[19]  P. Cohen.  Protein kinases — the major drug targets of the twenty-first century? , 2002, Nature reviews. Drug discovery.

[20]  Xingyu Wang,et al.  A Comparative Insight into Amprenavir Resistance of Mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 Protease Based on Thermodynamic Integration and MM-PBSA Methods , 2015, J. Chem. Inf. Model..

[21]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[22]  Dan Li,et al.  ADMET evaluation in drug discovery. 13. Development of in silico prediction models for P-glycoprotein substrates. , 2014, Molecular pharmaceutics.

[23]  John T. Powers,et al.  Targeting Bcr–Abl by combining allosteric with ATP-binding-site inhibitors , 2010, Nature.

[24]  Piotr Cieplak,et al.  Mechanism of influence of phosphorylation on serine 124 on a decrease of catalytic activity of human thymidylate synthase. , 2010, Bioorganic & medicinal chemistry.

[25]  Paul D. Martin,et al.  AZD6244 (ARRY-142886), a potent inhibitor of mitogen-activated protein kinase/extracellular signal-regulated kinase kinase 1/2 kinases: mechanism of action in vivo, pharmacokinetic/pharmacodynamic relationship, and potential for combination in preclinical models , 2007, Molecular Cancer Therapeutics.

[26]  Heather A. Carlson,et al.  Development of polyphosphate parameters for use with the AMBER force field , 2003, J. Comput. Chem..

[27]  Weiliang Zhu,et al.  A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation , 2015, Journal of biomolecular structure & dynamics.

[28]  Wei Wang,et al.  Revealing the binding modes and the unbinding of 14-3-3σ proteins and inhibitors by computational methods , 2015, Scientific Reports.

[29]  Xin Zhang,et al.  Molecular dynamics simulation reveals preorganization of the chloroplast FtsY towards complex formation induced by GTP binding. , 2011, Journal of structural biology.

[30]  Michael J. Bodkin,et al.  Accurate calculation of the absolute free energy of binding for drug molecules , 2015, Chemical science.

[31]  A. Bode,et al.  Caffeic Acid Directly Targets ERK1/2 to Attenuate Solar UV-Induced Skin Carcinogenesis , 2014, Cancer Prevention Research.

[32]  Xiao He,et al.  Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model. , 2013, Journal of chemical theory and computation.

[33]  S. Knapp,et al.  A unique inhibitor binding site in ERK1/2 is associated with slow binding kinetics , 2014, Nature chemical biology.

[34]  D. Sabatini,et al.  An ATP-competitive Mammalian Target of Rapamycin Inhibitor Reveals Rapamycin-resistant Functions of mTORC1* , 2009, Journal of Biological Chemistry.

[35]  Xin Zhang,et al.  Molecular dynamics simulation of SRP GTPases: Towards an understanding of the complex formation from equilibrium fluctuations , 2010, Proteins.

[36]  G. Labesse,et al.  Combining 'dry' co-crystallization and in situ diffraction to facilitate ligand screening by X-ray crystallography. , 2015, Acta crystallographica. Section D, Biological crystallography.

[37]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[38]  Weiliang Zhu,et al.  A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings , 2013, Journal of Computer-Aided Molecular Design.

[39]  Weiliang Zhu,et al.  Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation , 2015, Scientific Reports.

[40]  Zhichao Zhang,et al.  Mechanism of Mcl‐1 Conformational Regulation Upon Small Molecule Binding Revealed by Molecular Dynamic Simulation , 2016, Chemical biology & drug design.

[41]  Conformational selection or induced fit for Brinker and DNA recognition. , 2011, Physical chemistry chemical physics : PCCP.

[42]  E. Krebs,et al.  An adenosine 3',5'-monophosphate-dependant protein kinase from rabbit skeletal muscle. , 1968, The Journal of biological chemistry.

[43]  Alessandro Curioni,et al.  New Scoring Functions for Virtual Screening from Molecular Dynamics Simulations with a Quantum-Refined Force-Field (QRFF-MD). Application to Cyclin-Dependent Kinase 2 , 2006, J. Chem. Inf. Model..

[44]  A. Mildvan,et al.  ATP-binding site of adenylate kinase: mechanistic implications of its homology with ras-encoded p21, F1-ATPase, and other nucleotide-binding proteins. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[45]  W. Kolch,et al.  Mechanism of Suppression of the Raf/MEK/Extracellular Signal-Regulated Kinase Pathway by the Raf Kinase Inhibitor Protein , 2000, Molecular and Cellular Biology.

[46]  M. Lindvall,et al.  Ligand efficient tetrahydro-pyrazolopyridines as inhibitors of ERK2 kinase. , 2015, Bioorganic & medicinal chemistry letters.

[47]  Jiye Shi,et al.  Exploring transition pathway and free-energy profile of large-scale protein conformational change by combining normal mode analysis and umbrella sampling molecular dynamics. , 2014, The journal of physical chemistry. B.

[48]  J. Zou,et al.  Effect of mutation K85R on GSK-3beta: Molecular dynamics simulation. , 2008, Biochemical and biophysical research communications.

[49]  Robert I Cukier,et al.  Molecular dynamics of apo-adenylate kinase: a principal component analysis. , 2006, The journal of physical chemistry. B.

[50]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[51]  Samuel H. Wilson,et al.  Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects. , 2007, Journal of molecular biology.

[52]  Shao-Yong Lu,et al.  Molecular Modeling and Molecular Dynamics Simulation Studies of the GSK3β/ATP/Substrate Complex: Understanding the Unique P+4 Primed Phosphorylation Specificity for GSK3β Substrates , 2011, J. Chem. Inf. Model..

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

[54]  Ting Zhou,et al.  Kinase selectivity potential for inhibitors targeting the ATP binding site: a network analysis , 2010, Bioinform..

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

[56]  A. Mclachlan.  Gene duplications in the structural evolution of chymotrypsin. , 1979, Journal of molecular biology.

[57]  Ruth Nussinov,et al.  The Mechanism of ATP-Dependent Allosteric Protection of Akt Kinase Phosphorylation. , 2015, Structure.

[58]  P. Shapiro,et al.  Structure of extracellular signal-regulated kinase 2 in complex with ATP and ADP. , 2012, Acta crystallographica. Section F, Structural biology and crystallization communications.

[59]  L. Zawel,et al.  Dissecting Therapeutic Resistance to ERK Inhibition , 2016, Molecular Cancer Therapeutics.

[60]  Guohui Li,et al.  Coarse-Grained Modeling of Nucleic Acids Using Anisotropic Gay-Berne and Electric Multipole Potentials. , 2016, Journal of chemical theory and computation.

[61]  Jian Zhang,et al.  Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations. , 2005, Journal of the American Chemical Society.

[62]  K. Guan,et al.  Properties of MEKs, the kinases that phosphorylate and activate the extracellular signal-regulated kinases. , 1993, The Journal of biological chemistry.

[63]  Hualiang Jiang,et al.  Mechanically induced titin kinase activation studied by force-probe molecular dynamics simulations. , 2005, Biophysical journal.

[64]  T. Kinoshita,et al.  Identification of allosteric ERK2 inhibitors through in silico biased screening and competitive binding assay. , 2016, Bioorganic & medicinal chemistry letters.

[65]  Jiali Gao,et al.  Determination of the structure form of the fourth ligand of zinc in Acutolysin A using combined quantum mechanical and molecular mechanical simulation. , 2009, The journal of physical chemistry. B.

[66]  伊藤 善基.  Activation of mitogen-activated protein kinases/extracellular signal-regulated kinases in human hepatocellular carcinoma , 1998 .

[67]  Guodong Hu,et al.  Ligand selectivity of estrogen receptors by a molecular dynamics study. , 2014, European journal of medicinal chemistry.

[68]  Characterizing the Free‐Energy Landscape of MDM2 Protein–Ligand Interactions by Steered Molecular Dynamics Simulations , 2015, Chemical biology & drug design.

[69]  Ozlem Keskin,et al.  Ras Conformational Ensembles, Allostery, and Signaling. , 2016, Chemical reviews.

[70]  Mingjie Zhang,et al.  Phosphorylation of DCC by ERK2 is facilitated by direct docking of the receptor P1 domain to the kinase. , 2010, Structure.

[71]  Shulin Zhuang,et al.  Insights into unbinding mechanisms upon two mutations investigated by molecular dynamics study of GSK3β–axin complex: Role of packing hydrophobic residues , 2007, Proteins.

[72]  W. Rutter,et al.  Replacement of lysine residue 1030 in the putative ATP-binding region of the insulin receptor abolishes insulin- and antibody-stimulated glucose uptake and receptor kinase activity. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[73]  Nathanael Gray,et al.  Factors underlying sensitivity of cancers to small-molecule kinase inhibitors , 2009, Nature Reviews Drug Discovery.

[74]  R. Skeel,et al.  Langevin stabilization of molecular dynamics , 2001 .

[75]  D. Longo,et al.  Polo-like kinase 1-mediated phosphorylation of the GTP-binding protein Ran is important for bipolar spindle formation. , 2006, Biochemical and biophysical research communications.

[76]  P. Kollman,et al.  Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA. , 2001, Journal of the American Chemical Society.

[77]  P. Kollman,et al.  Computational study of protein specificity: The molecular basis of HIV-1 protease drug resistance , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[78]  A. R. Srinivasan,et al.  Quasi‐harmonic method for studying very low frequency modes in proteins , 1984, Biopolymers.

[79]  J. Boehm,et al.  p38 MAP kinases: key signalling molecules as therapeutic targets for inflammatory diseases , 2003, Nature Reviews Drug Discovery.

[80]  Karen Roberts,et al.  Structure-Guided Design of Highly Selective and Potent Covalent Inhibitors of ERK1/2. , 2015, Journal of medicinal chemistry.

[81]  Xin Zhang,et al.  Molecular dynamics simulations reveal structural coordination of Ffh‐FtsY heterodimer toward GTPase activation , 2011, Proteins.

[82]  Lisa L. Smith,et al.  AZD8055 is a potent, selective, and orally bioavailable ATP-competitive mammalian target of rapamycin kinase inhibitor with in vitro and in vivo antitumor activity. , 2010, Cancer research.

[83]  Michael Kracht,et al.  Targeting innate immunity protein kinase signalling in inflammation , 2009, Nature Reviews Drug Discovery.

[84]  K. Berg,et al.  Sustained EKR inhibition by EGFR targeting therapies is a predictive factor for synergistic cytotoxicity with PDT as neoadjuvant therapy , 2013 .

[85]  Tingjun Hou,et al.  Insight into Crizotinib Resistance Mechanisms Caused by Three Mutations in ALK Tyrosine Kinase using Free Energy Calculation Approaches , 2013, J. Chem. Inf. Model..

[86]  Shaoyong Lu,et al.  How calcium inhibits the magnesium‐dependent kinase gsk3β: A molecular simulation study , 2013, Proteins.