Understanding the Catalytic Mechanism of Xanthosine Methyltransferase in Caffeine Biosynthesis from QM/MM Molecular Dynamics and Free Energy Simulations

S-Adenosyl-l-methionine (SAM) dependent xanthosine methyltransferase (XMT) is the key enzyme that catalyzes the first methyl transfer in the caffeine biosynthesis pathway to produce the intermediate 7-methylxanthosine (7mXR). Although XMT has been a subject of extensive discussions, the catalytic mechanism and nature of the substrate involved in the catalysis are still unclear. In this paper, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy (potential of mean force or PMF) simulations are undertaken to determine the catalytic mechanism of the XMT-catalyzed reaction. Both xanthosine and its monoanionic form with N3 deprotonated are used as the substrates for the methylation. It is found that while the methyl group can be transferred to the monoanionic form of xanthosine with a reasonable free energy barrier (about 17 kcal/mol), that is not the case for the neutral xanthosine. The results suggest that the substrate for the first methylation step in the caffeine biosynthesis pathway is likely to be the monoanionic form of xanthosine rather than the neutral form as widely adopted. This conclusion is supported by the pKa value on N3 of xanthosine both measured in aqueous phase and calculated in the enzymatic environment. The structural and dynamics information from both the X-ray structure and MD simulations is also consistent with the monoanionic xanthosine scenario. The implications of this conclusion for caffeine biosynthesis are discussed.

[1]  Jayarama,et al.  The coffee genome provides insight into the convergent evolution of caffeine biosynthesis , 2014, Science.

[2]  Hong Guo,et al.  Quantum Mechanical/Molecular Mechanical Study of Catalytic Mechanism and Role of Key Residues in Methylation Reactions Catalyzed by Dimethylxanthine Methyltransferase in Caffeine Biosynthesis , 2014, J. Chem. Inf. Model..

[3]  Guohui Li,et al.  QM/MM MD and free energy simulations of the methylation reactions catalyzed by protein arginine methyltransferase PRMT3 , 2013 .

[4]  Hong Guo,et al.  Understanding Product Specificity of Protein Lysine Methyltransferases from QM/MM Molecular Dynamics and Free Energy Simulations: The Effects of Mutation on SET7/9 beyond the Tyr/Phe Switch , 2012, J. Chem. Inf. Model..

[5]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[6]  Elizabeth E Howell,et al.  A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs. , 2011, Biochemistry.

[7]  Feng Chen,et al.  QM/MM free energy simulations of salicylic acid methyltransferase: effects of stabilization of TS-like structures on substrate specificity. , 2011, The journal of physical chemistry. B.

[8]  Y. Yoshizawa,et al.  Essential Region for 3-N Methylation in N-Methyltransferases Involved in Caffeine Biosynthesis , 2010, Zeitschrift fur Naturforschung. C, Journal of biosciences.

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

[10]  Yingkai Zhang,et al.  How do SET-domain protein lysine methyltransferases achieve the methylation state specificity? Revisited by Ab initio QM/MM molecular dynamics simulations. , 2008, Journal of the American Chemical Society.

[11]  A. Crozier,et al.  Caffeine and related purine alkaloids: biosynthesis, catabolism, function and genetic engineering. , 2008, Phytochemistry.

[12]  M. Jackson,et al.  A preference for edgewise interactions between aromatic rings and carboxylate anions: the biological relevance of anion-quadrupole interactions. , 2007, The journal of physical chemistry. B.

[13]  Hao-Bo Guo,et al.  Mechanism of histone methylation catalyzed by protein lysine methyltransferase SET7/9 and origin of product specificity , 2007, Proceedings of the National Academy of Sciences.

[14]  A. McCarthy,et al.  The Structure of Two N-Methyltransferases from the Caffeine Biosynthetic Pathway1[W][OA] , 2007, Plant Physiology.

[15]  Xiaodong Zhang,et al.  Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations. , 2007, Biochemistry.

[16]  H. Sano,et al.  Caffeine biosynthesis and adenine metabolism in transgenic Coffea canephora plants with reduced expression of N-methyltransferase genes. , 2006, Phytochemistry.

[17]  T. Fujimura,et al.  The first committed step reaction of caffeine biosynthesis: 7‐methylxanthosine synthase is closely homologous to caffeine synthases in coffee (Coffea arabica L.)1 , 2003, FEBS letters.

[18]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

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

[20]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[21]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[22]  M. Karplus,et al.  Active site dynamics in protein molecules: A stochastic boundary molecular‐dynamics approach , 1985, Biopolymers.

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

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

[25]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[26]  Efthimios Kaxiras,et al.  A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method , 2001 .

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