{'fraction_non_alphanumeric': 0.032432432432432434, 'fraction_numerical': 0.02702702702702703, 'mean_word_length': 5.222689075630252, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Understanding the R882H mutation effects of DNA methyltransferase DNMT3A: a combination of molecular dynamics simulations and QM/MM calculations

DNA (cytosine-5)-methyltransferase 3A (DNMT3A), a key enzyme for de novo epigenetic methylation in human beings, was reported to undergo an R882H mutation in approximately 25% of M4/M5 subtype acute myeloid leukemia (AML) patients. In this work, a combination of classical molecular dynamics (MD) simulations and QM/MM calculation methods was utilized to reveal the molecular mechanism behind the activity attenuation caused by R882H mutation. We found that R882H mutation induces a “folded” conformation in the methyl donor S-adenosylmethionine (SAM) through different types of hydrogen bond formation at the terminal carbonyl oxygen atom and the hydroxyl O3′ atom of the ribose ring on SAM, with Arg891 as a mediator. Energetically, both the pre-reaction state (PRS) and transition state (TS) were stabilized in the R882H mutant. However, the energy barrier of the rate-determining step from the PRS to the TS was calculated to be roughly 1.0 kcal mol−1 larger in the R882H mutant than the WT. Also, a dynamic transformation occurred along the helix where R882H was located, tending to manifest in a quasi-“Newton's cradle” manner from the mutational site to the active site residues of DNMT3A. Our computational results provided molecular insights into the pathogenic R882H mutation and advanced the understanding of its mechanism.