Unraveling the Reaction Mechanism of Enzymatic C5-Cytosine Methylation of DNA. A Combined Molecular Dynamics and QM/MM Study of Wild Type and Gln119 Variant

详细信息    查看全文

• 作者：Juan Aranda ; Kirill Zinovjev ; Katarzyna Świderek ; Maite Roca ; Iñaki Tuñón
• 刊名：ACS Catalysis
• 出版年：2016
• 出版时间：May 6, 2016
• 年：2016
• 卷：6
• 期：5
• 页码：3262-3276
• 全文大小：642K
• 年卷期：0
• ISSN：2155-5435

文摘

M.HhaI is a DNA methyltransferase from Haemophilus hemolyticus that catalyzes the transfer of a methyl group from S-adenosyl-l-methionine (SAM) to the C₅ position of a cytosine. This enzyme is a paradigmatic model for C₅ DNA methyltransferases due to its major homology to mammalian enzymes and to the availability of high-resolution structures of the DNA–enzyme complex. In spite of the number of experimental and theoretical analyses carried out for this system, many mechanistic details remain unraveled. We have used full atomistic classical molecular dynamics simulations to explore the protein–SAM–DNA ternary complex, where the target cytosine base is flipped out into the active site for both the wild type and Glu119Gln mutant. The relaxed structure was used for a combined quantum mechanics/molecular mechanics exploration of the reaction mechanism using the string method. Exploration of multidimensional free energy surfaces allowed determining a complete picture of the reaction mechanism, in agreement with experimental observations. In our proposal Cys81 becomes unprotonated after formation of the ternary complex. This residue then attacks the C₆ position of the target cytosine, establishing a fast and reversible equilibrium. Methyl transfer from SAM to position C₅ takes place after the previous addition, and the transition state is stabilized by means of hydrogen bond interactions established between the cytosine base and Glu119. Although previous proposals suggest that a hydroxide anion could be the base abstracting the proton from the C₅ position, our calculations show that the free energy cost needed to have this anion in the active site is too high. Instead, a crystallographic water molecule can remove the excess proton once the cytosine has been activated by means of a proton transfer from Glu119 to the N₃ position. This observation explains the consequences of mutations of this residue on the pre-steady-state and steady-state rate constants. Our proposal then clarifies the role of Glu119 and identifies the nature of the base in charge of proton abstraction, two issues that have been the subject of a long debate in the literature.