Quantum mechanics study and Monte Carlo simulation on the hydrolytic deamination of 5-methylcytosine glycol

Abstract

The efficient formation of 5-methylcytosine glycol (mCg) and its facile deamination to thymine glycol (Tg) may account for the prevalent C → T transition mutation found at methylated CpG site (mCpG) in human p53 gene, a hallmark for many types of human tumors. In this work, the hydrolytic deamination of mCg was investigated at the MP2 and B3LYP levels of theory using the 6-311G(d,p) basis set. In the gas phase, three pathways were explored, paths A–C, and it indicates that the direct deamination of mCg with H₂O by either pathway is unlikely because of the high activation free energies involved in the rate-determining steps, the formation of the tetrahedral intermediate for paths A and B as well as the formation of the Tg tautomer for path C. In aqueous solution, the role of the water molecules in the deamination of mCg with H₂O was analyzed in two separate parts: the direct participation of one water molecule in the reaction pathway, called the water-assisted mechanism; and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for mCg and the cluster of two water molecules by quantum mechanical calculations in the gas phase. This indicates that the presence of the auxiliary water molecule significantly contributes to decreasing all the activation free energies. The bulk solution effect on the water-assisted mechanism was included by free energy perturbation implemented on Monte Carlo simulations, which is found to be substantial and decisive in the deamination mechanism of mCg. In this case, the water-assisted path A is the most plausible mechanism reported for the deamination of mCg, where the calculated activation free energy (22.6 kcal mol⁻¹ at B3LYP level of theory) agrees well with the experimentally determined activation free energy (24.8 kcal mol⁻¹). The main striking results of the present DFT computational study which is in agreement with previous experimental data is the higher rate of deamination displayed by mCg residues with respect to 5-methylcytosine (mC) bases, which supports that the deamination of mCg contributes significantly to the C → T transition mutation at mCpG dinucleotide site.