Catalytic Mechanism and Product Specificity of Rubisco Large Subunit Methyltransferase: QM/MM and MD Investigations
- 13 April 2007
- journal article
- Published by American Chemical Society (ACS) in Biochemistry
- Vol. 46 (18) , 5505-5514
- https://doi.org/10.1021/bi700119p
Abstract
Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT·Lys-NH3+·AdoMet), but is not present immediately after methyl transfer (LSMT·Lys-N(Me)H2+·AdoHcy). The water channel allows proton dissociation from both LSMT·AdoMet·Lys-NH3+ and LSMT·AdoMet·Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT·AdoMet·Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet → Lys-N(Me)H2+ + AdoHcy) is ΔG‡ = 22.8 ± 3.3 kcal/mol. This ΔG‡ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 × 10-5 s-1). The calculated ΔG‡ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H → AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 ± 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 × 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet → Lys-N(Me)3+ + AdoHcy) is associated with an allowed ΔG‡ of 25.9 ± 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.Keywords
This publication has 13 references indexed in Scilit:
- Computational Approaches: Reaction Trajectories, Structures, and Atomic Motions. Enzyme Reactions and ProficiencyChemical Reviews, 2006
- Catalytic Roles for Carbon-Oxygen Hydrogen Bonding in SET Domain Lysine MethyltransferasesJournal of Biological Chemistry, 2006
- Structural Insights of the Specificity and Catalysis of a Viral Histone H3 Lysine 27 MethyltransferaseJournal of Molecular Biology, 2006
- Structural and Sequence Motifs of Protein (Histone) Methylation EnzymesAnnual Review of Biophysics, 2005
- Structure and Catalytic Mechanism of a SET Domain Protein MethyltransferaseCell, 2002
- A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equationJournal of Computational Chemistry, 1991
- Solvent effects on protein motion and protein effects on solvent motionJournal of Molecular Biology, 1989
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics, 1977