SP20 Phosphorylation Reaction Catalyzed by Protein Kinase A: QM/MM Calculations Based on Recently Determined Crystallographic Structures

The cAMP-dependent protein kinase (PKA) has been the most studied human protein kinase ever. Very recently, new X-ray crystallographic structures in which the SP20 substrate has been trapped in the ternary complex with PKA before and after the phosphoryl transfer have provided a few tentative snapshots of the evolution of the enzyme system along the catalytic reaction. In the present paper, we have studied the dissociative and the associative mechanisms for the phosphorylation reaction of the SP20 substrate catalyzed by PKA by means of MP2/aug-cc-pVTZ/CHARMM//B3LYP/6-31+G(d)/CHARMM electronic structure calculations using a complete solvated model of the PKAc-ATPMg2–SP20 system. Our results demonstrate that the dissociative mechanism (involving two consecutive steps: phosphoryl transfer and back protonation of the phosphorylated substrate) is clearly more favorable than the associative one. A comparison of Kemptide with SP20 shows that the catalytic mechanism is not substrate-dependent. However, the produc...

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

[2]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[3]  Àngels González-Lafont,et al.  A QM/MM study of the phosphoryl transfer to the Kemptide substrate catalyzed by protein kinase A. The effect of the phosphorylation state of the protein on the mechanism. , 2011, Physical chemistry chemical physics : PCCP.

[4]  J. Zheng,et al.  Structure of a peptide inhibitor bound to the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase. , 1991, Science.

[5]  Àngels González-Lafont,et al.  A QM/MM study of Kemptide phosphorylation catalyzed by protein kinase A. The role of Asp166 as a general acid/base catalyst. , 2015, Physical chemistry chemical physics : PCCP.

[6]  Susan S. Taylor,et al.  2.2 A refined crystal structure of the catalytic subunit of cAMP-dependent protein kinase complexed with MnATP and a peptide inhibitor. , 1993, Acta crystallographica. Section D, Biological crystallography.

[7]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[8]  Walter Thiel,et al.  Linear scaling geometry optimisation and transition state search in hybrid delocalised internal coordinates , 2000 .

[9]  Susan S. Taylor,et al.  Low- and room-temperature X-ray structures of protein kinase A ternary complexes shed new light on its activity. , 2012, Acta crystallographica. Section D, Biological crystallography.

[10]  Ajit Banerjee,et al.  Search for stationary points on surfaces , 1985 .

[11]  P. Langan,et al.  Phosphoryl Transfer Reaction Snapshots in Crystals , 2015, The Journal of Biological Chemistry.

[12]  Susan S. Taylor,et al.  Structure and Allostery of the PKA RIIβ Tetrameric Holoenzyme , 2012, Science.

[13]  W Smith,et al.  DL_POLY_2.0: a general-purpose parallel molecular dynamics simulation package. , 1996, Journal of molecular graphics.

[14]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[15]  J. C. Slater.  A Simplification of the Hartree-Fock Method , 1951 .

[16]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

[17]  Susan S. Taylor,et al.  cAMP‐dependent protein kinase: Crystallographic insights into substrate recognition and phosphotransfer , 1994, Protein science : a publication of the Protein Society.

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

[19]  M. Field,et al.  Insights into the phosphoryl-transfer mechanism of cAMP-dependent protein kinase from quantum chemical calculations and molecular dynamics simulations. , 2004, Journal of the American Chemical Society.

[20]  D. M. Jacobsen,et al.  Briefly bound to activate: transient binding of a second catalytic magnesium activates the structure and dynamics of CDK2 kinase for catalysis. , 2011, Structure.

[21]  J. Stephen Binkley,et al.  Theoretical models incorporating electron correlation , 2009 .

[22]  J Andrew McCammon,et al.  How does activation loop phosphorylation modulate catalytic activity in the cAMP‐dependent protein kinase: A theoretical study , 2006, Protein science : a publication of the Protein Society.

[23]  J. Adams,et al.  Participation of ADP dissociation in the rate-determining step in cAMP-dependent protein kinase. , 1997, Biochemistry.

[24]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[25]  Michael J. Frisch,et al.  MP2 energy evaluation by direct methods , 1988 .

[26]  Susan S. Taylor,et al.  PKA: lessons learned after twenty years. , 2013, Biochimica et biophysica acta.

[27]  J. Nielsen,et al.  Predicting the open conformations of protein kinases using molecular dynamics simulations. , 2012, Biopolymers.

[28]  Àngels González-Lafont,et al.  Comparative study of the prereactive protein kinase A Michaelis complex with Kemptide substrate , 2007, J. Comput. Aided Mol. Des..

[29]  Susan S. Taylor,et al.  Insights into the Phosphoryl Transfer Catalyzed by cAMP-Dependent Protein Kinase: An X-ray Crystallographic Study of Complexes with Various Metals and Peptide Substrate SP20 , 2013, Biochemistry.

[30]  Alessandro Cembran,et al.  Allostery and binding cooperativity of the catalytic subunit of protein kinase A by NMR spectroscopy and molecular dynamics simulations. , 2012, Advances in protein chemistry and structural biology.

[31]  M Karplus,et al.  Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison , 1988, Proteins.

[32]  J. Baker.  An algorithm for the location of transition states , 1986 .

[33]  Marat Valiev,et al.  Phosphorylation reaction in cAPK protein kinase-free energy quantum mechanical/molecular mechanics simulations. , 2007 .

[34]  S. Taylor,et al.  Dynamics of cAMP-dependent protein kinase. , 2001, Chemical reviews.

[35]  C. Gibbs,et al.  Rational scanning mutagenesis of a protein kinase identifies functional regions involved in catalysis and substrate interactions. , 1991, The Journal of biological chemistry.

[36]  Arieh Warshel,et al.  Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces. , 2013, The journal of physical chemistry. B.

[37]  Douglas M. Jacobsen,et al.  Price To Be Paid for Two-Metal Catalysis: Magnesium Ions That Accelerate Chemistry Unavoidably Limit Product Release from a Protein Kinase , 2012, Journal of the American Chemical Society.

[38]  A. Warshel,et al.  Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping. , 2013, The journal of physical chemistry. B.

[39]  R. Huber,et al.  Phosphotransferase and substrate binding mechanism of the cAMP‐dependent protein kinase catalytic subunit from porcine heart as deduced from the 2.0 A structure of the complex with Mn2+ adenylyl imidodiphosphate and inhibitor peptide PKI(5‐24). , 1993, The EMBO journal.

[40]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[41]  J. Mccammon,et al.  How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study. , 2005, Journal of the American Chemical Society.

[42]  Hua Guo,et al.  Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies. , 2011, The journal of physical chemistry. B.

[43]  Susan S. Taylor,et al.  Metal-Free cAMP-Dependent Protein Kinase Can Catalyze Phosphoryl Transfer , 2014, Biochemistry.

[44]  Paul Sherwood,et al.  Zeolite structure and reactivity by combined quantum-chemical- classical calculations , 1999 .

[45]  J. Zheng,et al.  Crystal structure of the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase. , 1991, Science.

[46]  Susan S. Taylor,et al.  Molecular Features of Product Release for the PKA Catalytic Cycle , 2014, Biochemistry.

[47]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[48]  M. Vincent,et al.  Computer simulation of zeolite structure and reactivity using embedded cluster methods , 1997 .

[49]  Àngels González-Lafont,et al.  A QM/MM study of the associative mechanism for the phosphorylation reaction catalyzed by protein kinase A and its D166A mutant , 2014, Journal of Computer-Aided Molecular Design.

[50]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

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

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

[53]  S. Taylor,et al.  Serine-53 at the tip of the glycine-rich loop of cAMP-dependent protein kinase: role in catalysis, P-site specificity, and interaction with inhibitors. , 2000, Biochemistry.

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

[55]  S. C. Rogers,et al.  QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis , 2003 .

[56]  Susan S. Taylor,et al.  Phosphoryl transfer by protein kinase A is captured in a crystal lattice. , 2013, Journal of the American Chemical Society.

[57]  Susan S. Taylor,et al.  Crystal structure of a transition state mimic of the catalytic subunit of cAMP-dependent protein kinase , 2002, Nature Structural Biology.

[58]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[59]  Tao Yu,et al.  Dynamics connect substrate recognition to catalysis in protein kinase A. , 2010, Nature chemical biology.

[60]  Jiali Gao,et al.  Dynamically committed, uncommitted, and quenched states encoded in protein kinase A revealed by NMR spectroscopy , 2011, Proceedings of the National Academy of Sciences.

[61]  S. Taylor,et al.  Identification of a partially rate-determining step in the catalytic mechanism of cAMP-dependent protein kinase: a transient kinetic study using stopped-flow fluorescence spectroscopy. , 1997, Biochemistry.

[62]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[63]  Ian H. Williams,et al.  Transition-state structural refinement with GRACE and CHARMM: Flexible QM/MM modelling for lactate dehydrogenase , 1999 .

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

[65]  Àngels González-Lafont,et al.  Influence of the enzyme phosphorylation state and the substrate on PKA enzyme dynamics. , 2012, Biophysical chemistry.

[66]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[67]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[68]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[69]  Àngels González-Lafont,et al.  Benchmark calculations on models of the phosphoryl transfer reaction catalyzed by protein kinase A , 2009 .

[70]  M. Valiev,et al.  The role of the putative catalytic base in the phosphoryl transfer reaction in a protein kinase: first-principles calculations. , 2003, Journal of the American Chemical Society.

[71]  J. Adams,et al.  Kinetic and catalytic mechanisms of protein kinases. , 2001, Chemical reviews.