Search

The DNMT3A DNA methyltransferase is an important epigenetic enzyme that is frequently mutated in cancers, particularly in AML. The heterozygous R736H mutation in the FF-interface of the tetrameric enzyme is the second most frequently observed DNMT3A cancer mutation, but its pathogenic mechanism is unclear. We show here that R736H leads to a moderate reduction in catalytic activity of 20-40% depending on the substrate, but no changes in CpG specificity, flanking sequence preferences and subnuclear localization. Strikingly, R736H showed a very strong stimulation by DNMT3L and the R736H/DNMT3L complex was 3-fold more active than WT/DNMT3L. Similarly, formation of mixed WT/R736H FF-interfaces led to an increased activity. R736H/DNMT3L and mixed R736H/DNMT3A WT FF-interfaces were less stable than interfaces not involving R736H, suggesting that an increased flexibility of the mixed interfaces stimulates catalytic activity. Our data suggest that aberrant activity of DNMT3A R736H may lead to DNA hypermethylation in cancer cells which could cause changes in gene expression.

Somatic R882H DNMT3A mutations occur frequently in AML, but their pathogenic mechanism is unclear. As R882H mutations usually are heterozygous, wildtype (WT) and R882H subunits co-exist in affected cells. R882 is located in the RD interface of DNMT3A tetramers, which forms the DNA binding site. R882H causes strong changes in the flanking sequence preferences of DNMT3A. Here, we analyzed flanking sequence preferences for CGNNNN sites showing that most disfavored sites are methylated 4-5 fold slower by R882H than WT, while it methylates most preferred sites 2-fold faster. Overall, R882H was more active than WT at 13% and less active at 52% of all CGNNNN sites. We prepared mixed DNMT3A heterotetramers containing WT and R882H subunits and show that mixed complexes preferentially assemble with an R882H/R882H RD interface. Structural comparisons and MD simulations confirmed the conclusion that the R882H RD interface is more stable than that of WT, in part because H882 forms an inter-subunit contact in the RD interface, while R882 contacts the DNA. As the subunits at the RD interface contribute the two active centers to the DNMT3A tetramer, R882H characteristic flanking sequence preferences of DNMT3A were observed in mixed tetrameric complexes containing WT and R882H subunits, and they are not diluted by the "averaged" effects of mixed or WT interfaces. Hence, R882H has a dominant effect on the flanking sequence preferences and other catalytic properties of DNMT3A in samples containing WT and R882H subunits, which may explain its pathogenic effect in heterozygous state.

Somatic DNMT3A mutations at R882 are frequently observed in AML patients including the very abundant R882H, but also R882C, R882P and R882S. Using deep enzymology, we show here that DNMT3A-R882H has more than 70-fold altered flanking sequence preferences when compared with wildtype DNMT3A. The R882H flanking sequence preferences mainly differ on the 3′ side of the CpG site, where they resemble DNMT3B, while 5′ flanking sequence preferences resemble wildtype DNMT3A, indicating that R882H behaves like a DNMT3A/DNMT3B chimera. Investigation of the activity and flanking sequence preferences of other mutations of R882 revealed that they cause similar effects. Bioinformatic analyses of genomic methylation patterns focusing on flanking sequence effects after expression of wildtype DNMT3A and R882H in human cells revealed that genomic methylation patterns reflect the details of the altered flanking sequence preferences of R882H. Concordantly, R882H specific hypermethylation in AML patients was strongly correlated with the R882H flanking sequence preferences. R882H specific DNA hypermethylation events in AML patients were accompanied by R882H specific mis-regulation of several genes with strong cancer connection, which are potential downstream targets of R882H. In conclusion, our data provide novel and detailed mechanistic understanding of the pathogenic mechanism of the DNMT3A R882H somatic cancer mutation.

Recently, it has been discovered that different DNA-(cytosine C5)-methyltransferases including DNMT3A generate low levels of 3mC [Rosic et al. (2018), Nat. Genet., 50, 452–459]. This reaction resulted in the co-evolution of DNMTs and ALKB2 DNA repair enzymes, but its mechanism remained elusive. Here, we investigated the catalytic mechanism of DNMT3A for cytosine N3 methylation. We generated several DNMT3A variants with mutated catalytic residues and measured their activities in 5mC and 3mC generation by liquid chromatography linked to tandem mass spectrometry. Our data suggest that the methylation of N3 instead of C5 is caused by an inverted binding of the flipped cytosine target base into the active-site pocket of the DNA methyltransferase, which is partially compatible with the arrangement of catalytic amino acid residues. Given that all DNA-(cytosine C5)-methyltransferases have a common catalytic mechanism, it is likely that other enzymes of this class generate 3mC following the same mechanism.

DNMT3A/3L heterotetramers contain two active centers binding CpG sites at 12 bp distance, however their interaction with DNA not containing this feature is unclear. Using randomized substrates, we observed preferential co-methylation of CpG sites with 6, 9 and 12 bp spacing by DNMT3A and DNMT3A/3L. Co-methylation was favored by AT bases between the 12 bp spaced CpG sites consistent with their increased bending flexibility. SFM analyses of DNMT3A/3L complexes bound to CpG sites with 12 bp spacing revealed either single heterotetramers inducing 40° DNA bending as observed in the X-ray structure, or two heterotetramers bound side-by-side to the DNA yielding 80° bending. SFM data of DNMT3A/3L bound to CpG sites spaced by 6 and 9 bp revealed binding of two heterotetramers and 100° DNA bending. Modeling showed that for 6 bp distance between CpG sites, two DNMT3A/3L heterotetramers could bind side-by-side on the DNA similarly as for 12 bp distance, but with each CpG bound by a different heterotetramer. For 9 bp spacing our model invokes a tetramer swap of the bound DNA. These additional DNA interaction modes explain how DNMT3A and DNMT3A/3L overcome their structural preference for CpG sites with 12 bp spacing during the methylation of natural DNA., Graphical Abstract Graphical AbstractNGS methylation analysis with pool of randomized CpG double site substrates and SFM imaging of enzyme-DNA complexes with example substrates revealed novel interaction modes of DNMT3A and DNMT3A/3L with CpG double site substrates.

Methylation at the 5 position of cytosine in DNA (5meC), is a key epigenetic mark in eukaryotes. Once introduced, 5meC can be maintained through DNA replication due to the activity of “maintenance” DNA methyltransferases.. Despite their ancient origin, DNA methylation pathways differ widely across metazoans, such that 5meC is either confined to transcribed genes or lost altogether in several lineages. Here we use comparative epigenomics to investigate the evolution of DNA methylation. Although the model nematode C. elegans has lost DNA methylation, more basal nematodes retain cytosine DNA methylation, targeted to repeat loci. Unexpectedly, we find that DNA methylation coevolves with the DNA alkylation repair enzyme ALKB2 across eukaryotes. We further show that DNA methyltransferases introduce the toxic lesion 3meC into DNA both in vitro and in vivo. Alkylation damage is thus intrinsically associated with DNA methyltransferase activity, and this may promote the loss of DNA methylation in many species.

DNA interacting enzymes recognize their target sequences embedded in variable flanking sequence context. The influence of flanking sequences on enzymatic activities of DNA methyltransferases (DNMTs) can be systematically studied with “deep enzymology” approaches using pools of double-stranded DNA substrates, which contain target sites in random flanking sequence context. After incubation with DNMTs and bisulfite conversion, the methylation states and flanking sequences of individual DNA molecules are determined by NGS. Deep enzymology studies with different human and mouse DNMTs revealed strong influences of flanking sequences on their CpG and non-CpG methylation activity and the structures of DNMT-DNA complexes. Differences in flanking sequence preferences of DNMT3A and DNMT3B were shown to be related to the prominent role of DNMT3B in the methylation of human SATII repeat elements. Mutational studies in DNMT3B discovered alternative interaction networks between the enzyme and the DNA leading to a partial equalization of the effects of different flanking sequences. Structural studies in DNMT1 revealed striking correlations between enzymatic activities and flanking sequence dependent conformational changes upon DNA binding. Correlation of the biochemical data with cellular methylation patterns demonstrated that flanking sequence preferences are an important parameter that influences genomic DNA methylation patterns together with other mechanisms targeting DNMTs to genomic sites.

DNMT3A R882 mutations are frequently observed in AML including the abundant R882H and the rare R882C, R882P and R882S. Using deep enzymology we show here that the DNMT3A-R882H has more than 70-fold altered flanking sequence preferences when compared with wildtype DNMT3A. The R882H flanking sequence preferences mainly differ on the 3’ side of the CpG site, where they resemble DNMT3B, while 5’ flanking sequence preferences of R882H resemble wildtype DNMT3A, indicating that R882H behaves like a DNMT3A/DNMT3B chimera. Activities and flanking sequence preferences of R882C, R882P and R882S were determined as well. Genomic methylation patterns after expression of wildtype DNMT3A and R882H in human cells reflect the flanking sequence preferences. R882H specific hypermethylation in AML patients are correlated with R882H flanking sequence preferences. The hypermethylation events are accompanied by R882H specific misregulation of several genes with strong cancer connection in AML patients, which are potential downstream targets of R882H.

The PWWP domain of DNMT3 DNA methyltransferases binds to histone H3 tails containing methylated K36, and this activity is important for heterochromatic targeting. Here, we show that the PWWP domain of mouse DNMT3A binds to H3K36me2 and H3K36me3 with a slight preference for H3K36me2. PWWP domains have also been reported to bind to DNA, and the close proximity of H3K36 and nucleosomal DNA suggests a combined binding to H3K36me2/3 and DNA. We show here that the DNMT3A PWWP domain binds to DNA with a weak preference for AT-rich sequences and that the designed charge reversal R362E mutation disrupts DNA binding. The K295E mutation, as well as K295I recently identified in paraganglioma, a rare neuroendocrine neoplasm, disrupts both DNA and H3K36me2/3 binding, which is in agreement with the proximity of K295 to residues involved in K36me2/3 methyllysine binding. Nucleosome pulldown experiments show that DNA binding and H3K36me2/3 binding are important for the interaction of the DNMT3A PWWP domain with nucleosomes. Localization studies of transiently transfected fluorescently-tagged wild-type and PWWP-mutated full-length DNMT3A indicate that both interactions contribute to the subnuclear localization of DNMT3A in mouse cells. In summary, our data demonstrate that the combined binding of the DNMT3A PWWP domain to the H3 tail containing K36me2/3 and to the nucleosomal or linker DNA is important for its chromatin interaction and subnuclear targeting of DNMT3A in living cells.

The DNMT3A R882H mutation is frequently observed in acute myeloid leukemia (AML). It is located in the subunit and DNA binding interface of DNMT3A and has been reported to cause a reduction in activity and dominant negative effects. We investigated the mechanistic consequences of the R882H mutation on DNMT3A showing a roughly 40% reduction in overall DNA methylation activity. Biochemical assays demonstrated that R882H does not change DNA binding affinity, protein stability or subnuclear distribution of DNMT3A. Strikingly, DNA methylation experiments revealed pronounced changes in the flanking sequence preference of the DNMT3A-R882H mutant. Based on these results, different DNA substrates with selected flanking sequences were designed to be favored or disfavored by R882H. Kinetic analyses showed that the R882H favored substrate was methylated by R882H with 45% increased rate when compared with wildtype DNMT3A, while methylation of the disfavored substrate was reduced 7-fold. Our data expand the model of the potential carcinogenic effect of the R882H mutation by showing CpG site specific activity changes. This result suggests that R882 is involved in the indirect readout of flanking sequence preferences of DNMT3A and it may explain the particular enrichment of the R882H mutation in cancer patients by revealing mutation specific effects.

The Dnmt3a DNA methyltransferase has been shown to bind cooperatively to DNA and to form large multimeric protein/DNA fibers. However, it has also been reported to methylate DNA in a processive manner, a property that is incompatible with protein/DNA fiber formation. We show here that the DNA methylation rate of Dnmt3a increases more than linearly with increasing enzyme concentration on a long DNA substrate, but not on a short 30-mer oligonucleotide substrate. We also show that addition of a catalytically inactive Dnmt3a mutant, which carries an amino acid exchange in the catalytic center, increases the DNA methylation rate by wild type Dnmt3a on the long substrate but not on the short one. In agreement with this finding, preincubation experiments indicate that stable protein/DNA fibers are formed on the long, but not on the short substrate. In addition, methylation experiments with substrates containing one or two CpG sites did not provide evidence for a processive mechanism over a wide range of enzyme concentrations. These data clearly indicate that Dnmt3a binds to DNA in a cooperative reaction and that the formation of stable protein/DNA fibers increases the DNA methylation rate. Fiber formation occurs at low μm concentrations of Dnmt3a, which are in the range of Dnmt3a concentrations in the nucleus of embryonic stem cells. Understanding the mechanism of Dnmt3a is of vital importance because Dnmt3a is a hotspot of somatic cancer mutations one of which has been implicated in changing Dnmt3a processivity.