key: cord-0727315-th7ke15k
authors: Woodcock, Clayton B; Horton, John R; Zhang, Xing; Blumenthal, Robert M; Cheng, Xiaodong
title: Beta class amino methyltransferases from bacteria to humans: evolution and structural consequences
date: 2020-05-26
journal: Nucleic Acids Res
DOI: 10.1093/nar/gkaa446
sha: 36d2fb747537d7f9411d10b6166fb4b4ec50ed47
doc_id: 727315
cord_uid: th7ke15k

S-adenosyl-l-methionine dependent methyltransferases catalyze methyl transfers onto a wide variety of target molecules, including DNA and RNA. We discuss a family of methyltransferases, those that act on the amino groups of adenine or cytosine in DNA, have conserved motifs in a particular order in their amino acid sequence, and are referred to as class beta MTases. Members of this class include M.EcoGII and M.EcoP15I from Escherichia coli, Caulobacter crescentus cell cycle–regulated DNA methyltransferase (CcrM), the MTA1-MTA9 complex from the ciliate Oxytricha, and the mammalian MettL3-MettL14 complex. These methyltransferases all generate N6-methyladenine in DNA, with some members having activity on single-stranded DNA as well as RNA. The beta class of methyltransferases has a unique multimeric feature, forming either homo- or hetero-dimers, allowing the enzyme to use division of labor between two subunits in terms of substrate recognition and methylation. We suggest that M.EcoGII may represent an ancestral form of these enzymes, as its activity is independent of the nucleic acid type (RNA or DNA), its strandedness (single or double), and its sequence (aside from the target adenine).

Early analysis, of the amino acid sequences of 13 bacterial DNA methyltransferases (MTases) generating 5methylcytosine (5mC), revealed a set of ten conserved blocks of amino acid residues (1). These conserved motifs, numbered I to X from amino to carboxyl end, were found to have a constant linear order, simplifying their identification in protein sequences (particularly for the shorter or less-conserved motifs), though one alternative permutation of 5mC MTase motif order was later found (2, 3) . This aided in the discovery of mammalian 5mC-generating DNA methyltransferases Dnmt1 (4) and Dnmt3 (5) . The structural analysis of HhaI MTase, from the bacterium Haemophilus haemolyticus, allowed the functions of conserved sequence motifs to be inferred, particularly those responsible for SAM binding (motif I consensus: FxGxG) and for the chemistry of methyltransfer reaction on an inert carbon (motif IV consensus: PC), along with a varied target recognition domain (TRD) for specific binding of the substrate DNA sequence (6, 7) .

Not surprisingly, the same set of sequence motifs also occur in mammalian Dnmt2 (8) (9) (10) , a tRNA 5mC MTase (11) . The conservation in Dnmt sequence and structure reflects the conserved nature of SAM binding, which occurs so as to optimize the catalysis of methyl transfer onto cytosine-C5 in nucleic acids. Before and after the discovery that Dnmt2 homologs are actually tRNA MTases, a number of studies on the activity of Dnmt2 on DNA yielded kinetic rate constants varying between zero and very low, though the human enzyme can methylate cytidine in a short segment of single-stranded DNA ligated into a tRNA molecule (12) .

In addition to cytosine-C5, methylated bases on DNA include N4-methylcytosine (N4mC) and N6-methyladenine (N6mA), both of which are modified on their exocyclic primary amines. Bacteriophages, particularly lytic ones, modify their DNA in many additional ways, including other methylations (13) , but our focus here is on methylations that are consistent with normal cell physiology. Recognition of the base to be methylated, and access to it by catalytic residues, occurs after the base is swung out of the DNA duplex into a typical concave binding pocket, in a process known as 'base flipping' (14) (15) (16) . Interestingly, the DNA MTases that generate N4mC and N6mA do not group separately from one another on the basis of sequence (17) . This highlights the constraints of the chemistry of amino methylation as being evolutionarily dominant, instead of recognition of cytosine versus adenine in the active-site binding pocket. In fact, there is evidence of independent derivations of N4mC MTases from N6mA MTases (18) , and at least one MTase that generates N4mC can generate N6mA on a DNA substrate in which a C-to-A substitution is made at the target base (19) .

A multiple sequence alignment of 42 known DNA amino MTases revealed the existence of (at least) three classes of amino MTases, differing from one another by circular permutation in their order of motifs important for three essential functions: binding the methyl donor SAM, binding substrate DNA and catalyzing the chemical reaction between the donor and substrate (20) (Figure 1 ). Briefly, class ␣ is arranged in the order (N-to-C termini): motif I (SAM-binding), TRD (substrate recognition and binding), and motif IV (methylation catalysis). Class ␤ is arranged in the order: motif IV-TRD-motif I. Class ␥ is arranged in the order: motif I-motif IV-TRD. The class ␥ amino MTases, exemplified by M.TaqI (21, 22) , have the motif order comparable to that of the 5mC MTases (23) . The residues of motif IV vary with the target atom--the PC (proline-cysteine) motif in 5mC MTases is responsible for ring carbon-C5 methylation of cytosine, with the Cys sulfhydryl acting as a nucleophile to attack the cytosine ring and initiate the reaction (24) (25) (26) (27) . In contrast, motif IV for MTases catalyzing exocyclic amino methylation to generate N6mA or N4mC have the four-residue consensus (D/N/S)-P-P-(Y/F/W) [(Asp/Asn/Ser)-Pro-Pro-(Tyr/Phe/Trp)]; where the last bulky hydrophobic side chain stacks against the target base, while the first one abstracts a proton from the target amino group (22, 28) .

As was illustrated by Dnmt2, the apparent sequence and structural similarities do not reveal with certainty whether an enzyme acts on a particular methylation substrate (DNA, RNA or protein). Other examples include Escherichia coli HemK and human HemK2, which were thought to be DNA N6mA MTases (29, 30) , yet were found to be protein MTases active on glutamine (31) (32) (33) (34) (35) (36) . The human HemK2 is also active as a histone lysine MTase (37) . The common feature of the potential substrates is the amino group (NH 2 ) of protein glutamine or lysine (in the deprotonated state) or DNA adenine, at which all three methylation reactions are catalyzed by a NPPY motif (a subset of the motif IV possibilities).

While the MTase families differ in motif order, their structures are remarkably similar, comprising a sevenstranded ␤ sheet (1-to-7) with a central topological switchpoint between strands ␤1 and ␤4, and a characteristic reversed ␤ hairpin (␤6 and ␤7) at one end of the sheet next to strand ␤5 (Figure 2A and B). This topology allows for circular permutation, where the same structure simply has the break between its amino beginning and carboxyl end at different points (38) and, in fact, it is possible to circularly permute a MTase in the laboratory with retention of function (3). The two most significant motifs -motif I (FxGxG) for binding SAM and motif IV (DPPY) for binding substrate Ade -juxtapose the target N6 atom of Ade in line with the methyl group and sulfur atom of SAM for catalysis, and are positioned at the carboxyl ends of the two parallel neighboring strands ␤1 and ␤4 ( Figure 2C ). Based on the chemical similarity of the SAM-adenosyl and DNAadenosyl moieties, it makes sense that the two structural elements responsible for the binding (␤1-␤2 and ␤4-␤5) do, in fact, have striking similarity (20) . It was suggested that the original MTases arose after tandem gene duplication (39) converted a SAM-binding domain into a protein that bound two molecules of SAM (20) , and there is some evidence to support that hypothesis (40) .

Based on the proposal that the ancestral tandem duplication had two copies of motif I (represented as I-I), one evolutionary model is that the two halves diverged with one SAM-binding motif I changing to an adenine-binding motif IV ( Figure 1A ). If this change involved the amino-proximal motif I, it would yield the ancestor to Class ␤ MTases (IV-I), while conversion of the carboxyl-proximal SAM binding site would yield the ancestor to Class ␣ and ␥ MTases (I-IV). An alternative would be that just one of the two changed, with the other motif order occurring subsequently via a true circular permutation event (38,41) (not shown). Presumably, the ancestral MTase(s) as described here would methylate free adenine or a free adenine nucleotide, perhaps contributing to pre-protein metabolism (42, 43) .

This ancestral nucleotide MTase could then have gained the ability to modify DNA or RNA adenines following an additional fusion of the target (nucleic acid substrate) recognition domain (TRD, Figure 1B ). The three classes are roughly the same size in number of known members ( Figure 1C ). While the majority of 5mC and N4mC DNA MTases each fall into single classes defined by motif order, the N6mA MTases are fairly evenly distributed among the three classes (␣, ␤ and ␥ ) (20) . This could also be explained by the tandem duplication model in which the ancestral nucleic acid MTase(s) generated N6mA. In addition to considering MTase evolution, we also note here that the ␤ arrangement of motifs has intrinsic structural consequences that, we suggest, impose a unique requirement for this family of MTases to dimerize in order to function. We begin by comparing two MTases that play similar biochemical and physiological roles, but one of them is a ␤ MTase and the other is not.

The great majority of bacterial and archaeal DNA MTases are associated with restriction-modification (RM) systems, where the MTase protects a cell's own DNA from digestion by the paired (cognate) restriction endonuclease (44, 45) . RM systems are important for defense against bacteriophage predation (46, 47) , although they play other roles as well (48, 49) . Bacterial 'orphan' MTases are so named as, unlike most bacterial DNA MTases, they are not paired with a restriction endonuclease as part of a RM system (50) . Orphan MTases are sometimes involved in chromosome replication, DNA repair, and epigenetic gene regulation (51) . Examples of such regulatory orphan MTases include the DNA adenine MTase (Dam) in Escherichia coli (Gammaproteobacteria) and cell cycleregulated DNA MTase (CcrM) in Caulobacter crescentus (Alphaproteobacteria) which are, respectively, responsible for maintenance of daughter strand adenine methylation at GATC or GAnTC sequences (n = any nucleotide), immediately after their replication (52, 53) . We will first compare these two 'orphan' enzymes, and then consider ␤ MTases associated with RM systems, and finally examine the mammalian ␤ MTases.

Although the substrate sequences of Dam (GATC) and CcrM (GAnTC) are both palindromic, both enzymes act on hemi-methylated DNA substrates ( Figure 2G -H) -that is, sequences that already contain N6mA in one strand. Both enzymes can act on unmethylated sites, but most often encounter hemimethylated sites resulting from replication of fully-methylated sites, and both enzymes show a preference for hemimethylated over unmethylated sites (54, 55) .

Besides the aforementioned difference in the order of motifs (Dam is a member of class ␣ and CcrM is a class ␤ MTase), there are three major differences between the two enzymes. First, though both enzymes use an arginine to interact with a 5 guanine ( Figure 2G and H) and both Arg-Gua interactions are sequence-discriminatory contacts, these contacts are made to different strands relative to the substrate Ade. In Dam, the recognized Gua is on the opposite strand, and 2-base pairs away, from the target Ade; whereas in CcrM, the recognized Gua is on the same strand as, and adjacent to, the target Ade. Second, except for the flipped-out target adenine, the Dam-bound DNA conformation has intact intra-helical paired bases, whereas CcrM pulls the two DNA strands apart, creating a bubble comprising four enzyme-recognized, unpaired bases ( Figure 2I ). These two features of CcrM, strand separation and base recognition on the same strand that contains the target Ade, allow CcrM (but not Dam) to be active on both double-stranded (ds) and single-stranded (ss) DNA, as well as on mismatches within or immediately outside of the recognition sequence (56) . Neither enzyme is active on ssRNA. Lastly, Dam is a monomer under most conditions (57, 58) , while CcrM is a homodimer ( Figure 3A ) (59) , reflecting a broader phenomenon that we will discuss next.

The homodimeric subunit structure of CcrM has been observed previously in other MTases of structurallycharacterized class ␤ MTases. M.PvuII of Proteus vulgaris is (at this writing) the only structurally-characterized N4mC MTase (60) , and yielded the first structure for a MTase in the ␤ family. It methylates the central cytosine in its symmetric recognition sequence 5 -CAGCTG-3 . The M.PvuII structure, consistent with the circular permutation model, shares a common fold with MTases of other families, while having the major functional regions (particularly motifs I and IV) permuted into distinct linear order. The major unexpected finding was that M.PvuII forms a homodimer ( Figure 3B ). This dimeric feature was subsequently observed by others in other ␤ MTase structures. M.RsrI, from Rhodobacter sphaeroides (61,62) ( Figure 3C ), methylates the internal adenine of the palindromic DNA sequence GAATTC, and its dimeric structure was consistent with earlier biochemical evidence for this MTase (63, 64) . M1.MboII, from Moraxella bovis ( Figure 3D ), methylates the 3 adenine of an asymmetric sequence 5 -GAAGA-3 , and the homodimeric structure again supported by biochemical results (65) . Dimer formation was also observed in (72)). However, none of these class ␤ MTases were structurally characterized in complex with DNA, and their postulated TRDs were disordered in the absence of bound DNA, so the significance of the dimeric character was not grasped until the structure determination of M.EcoP15I.

M.EcoP15I is an archetype of the Type III RM family, and consists of two methylation (M) and 0-2 restriction (R) subunits (73) , resulting in M 2 , M 2 R 1 or M 2 R 2 complexes. The M.EcoP15I M-subunit functions as a dimer (74, 75) but recognizes an asymmetric sequence, 5 -CAGCAG-3 and methylates the 3 internal adenine of the A-containing strand of the dsDNA (the complementary strand is 5 -CTGCTG-3 ). Other Type III MTases have also been shown biochemically to form dimers in solution (76, 77) .

The M.EcoP15I structure is the first one for a class-␤ amino MTase bound to its substrate DNA, and it suggests a division of labor between two M subunits in terms of DNA recognition (one M subunit) and catalysis of methyltransfer (the other M subunit) (78) ( Figure 3G ). The DNArecognizing subunit of the homodimer provides side chains, such as Arg305 and Arg350, for recognition of guanines on the non-target (non-A-containing) strand ( Figure 3H ). These Arg-Gua interactions with the non-target strand are reminiscent of Escherichia coli Dam ( Figure 2G ), and (like Dam) limit M.EcoP15I to methylating dsDNA and not ss-DNA. In contrast, while differing roles for each subunit are also seen for CcrM, each protomer contacts primarily a different DNA strand. Specifically, one CcrM subunit binds the target strand, recognizes the target sequence and catalyzes methyl transfer, while the second molecule simply binds the non-target strand (in the case of a dsDNA substrate) ( Figure 3A ) (59) . This division of labor allows CcrM to methylate both ds and ssDNA.

Looking at these dimeric class ␤ MTase structures together, whether in the presence or absence of DNA, they share a striking feature--the catalytic domain of one subunit and the DNA binding domain of the second subunit are arranged so as to face each other, forming one integral binding surface appropriate for one DNA duplex (dashed circles in Figure 3A , C and G). In contrast, the catalytic site and DNA binding of the same subunit face away from one another, such that they cannot cooperate to accomplish the DNA recognition and methylation in one binding event. This is distinct from the association of MTases with other proteins in order to modulate their activity (see Discussion); and also implies that--while there may be kinetic or other reasons for which dimerization may be advantageous to other MTases (examples include mammalian Dnm3L-3a-3a-3L tetramer (79-81))--for the ␤ MTases dimerization is essential for activity. 

M.EcoGII is encoded in the genome of the pathogenic strain E. coli O104:H4 C227-11 (82) , which was responsible for a severe outbreak of hemorrhagic uremia in Europe (83) . The gene appears to reside within an integrated prophage genome, but is not expressed during normal bacterial growth. When the gene encoding M.EcoGII is expressed in vivo -using a high copy plasmid vector and a methylationdeficient E. coli host--extensive in vivo adenine methylation activity is revealed. It appears that M.EcoGII methylates adenine residues in any sequence context, in DNA and RNA, double-stranded and single stranded, and even adenine bases in either strand of a DNA/RNA-hybrid oligonucleotide duplex (84). There is not yet a M.EcoGII-DNA complex structure, though modeling based on a homolog (using threading with the M.EcoP15I-DNA structure) predicts a dimeric M.EcoGII complexed with one duplex (Figure 4A) . Strikingly, from Figure 4A it appears to be impossible for M.EcoGII to function as a monomer. If it does function as a dimer, which seems likely, it is surprising that there appears to be no quadratic rate dependence of activity on the enzyme concentration [see figure 4 in (84)]. In contrast, such rate dependence is clear for two other ␤ MTases, M.KpnI and M.RsrI [see figure 4 in (64) and (70)]. One possible explanation of this apparent paradox is that M.EcoGII has a particularly high dimerization constant. In fact, the known homodimer CcrM also shows a linear rate dependence on the enzyme concentration [see figure 7 in (85)].

Even if we assume that M.EcoGII functions as a dimer, however, it is unclear how the 'division of labor' between the two M.EcoGII subunits is accomplished, given its lack of sequence specificity, and the fact that recognition and methylation occur to the same nucleotide. Considering the aforementioned gene duplication model of a tandemlyduplicated SAM-binding protein having one SAM pocket diverge to bind adenine, M.EcoGII might be an example of such a primordial adenine MTase, lacking a well-defined and functional TRD.

However, while it is difficult to identify TRDs from sequence analysis alone (86) (87) (88) , there is no evidence that the region in which the TRD is expected to occur in M.EcoGII and its orthologs is unusual (e.g., smaller than in other MTases; not shown). Thus it is also possible that M.EcoGII began as a standard N6mA MTase and lost specificity, under selection for its role in the physiology of the phage that carries it. In this regard, it is noteworthy that some class-␤ RM MTases tested under certain solvent conditions can in fact methylate the DNA strand of a DNA/RNA hybrid duplex (89) .

In contrast to metazoa, DNA N6mA is abundant in various unicellular eukaryotes, including ciliates, the green alga Chlamydomonas, and early-diverging fungi (90) (91) (92) . A recent study identified four ciliate proteins--two class-␤-like MTases (MTA1 and MTA9) and two homeoboxlike DNA binding proteins--as being necessary for deposition of N6mA in the Oxytricha genome (93) . Instead of a gene-targeting approach, the authors used classic biochemistry by identifying candidate proteins, in nuclear extracts, that co-purified with DNA MTase activity. A conserved catalytic motif IV (DPPW) is preserved in MTA1, but not MTA9 (which has an eroded motif NALW) (Figure 4B) . The four-subunit ciliate MTase complex preferentially methylates ApT dinucleotides in dsDNA -in agreement with earlier observations that in the unicellular eukaryotes Saccharomyces cerevesiae, Chlamydomonas reinhardii and Tetrahymena thermophila, N6mA is enriched in ApT dinucleotides (92) within nucleosome linker regions near promoters (94) . This ApT sequence is symmetrical, like the CpG methylation target in mammals, suggesting a similar mechanism for maintaining methylation following DNA replication. Indeed, the MTase activity was even higher on hemimethylated than on unmethylated dsDNA, though it showed no activity on ssDNA or RNA (93) . Homology modeling suggests a heterodimeric association of MTA1 and MTA9 ( Figure 4C ), in analogy to mammalian MettL3-MettL14, which we will discuss next.

Detection of N6mA in mammalian DNA was reported only recently, but the role of MettL3 and MettL14 in generating N6mA in RNA has been known for some time. In mammalian (HeLa) cells, the mRNA N6mA MTase activity on the degenerate consensus sequence RRACH (R = purine, H is not a G) (95) requires at least two separate subunits, MT-A and MT-B (96) . MT-A is itself a multimeric protein, that contains a 70-kDa MT-A70 subunit (now known as MettL3, containing 580-residues and with a predicted molecular weight of 65 kDa). The amino acid sequence of MT-A70/MettL3 was noted to include conserved motifs (e.g. motif IV: DPPW) of the bacterial DNA MTases but not of known RNA MTases (97) . MettL3 and MettL14 form a heterodimeric complex (98) , and have been studied widely for their role in generating N6mA in RNA ((99,100) and references therein). Like the ciliate MTA1 and MTA9 (preceding section), mammalian MettL3 is the catalytically active subunit, while MettL14 has an eroded motif IV (EPPL) ( Figure 4B ). An earlier study suggested that MettL14 catalyzed RNA methylation (98) . Unlike the ciliate complex, which requires two additional DNA binding subunits, MettL3 itself contains two tandem CCCH-type zinc fingers within its polypeptide, located N-terminal to its MTase domain, that are necessary for RNA binding and thus enzymatic activity (101) . Four research groups have independently determined the structures of the heterodimeric complex of the MTase domains of MettL3 and MettL14 ( (102) (103) (104) and Structural Genomics Consortium) ( Figure  4D ). However, these partial complexes, lacking the CCCH- zinc fingers, possess no enzymatic activity and there is no RNA substrate-bound complex structure currently available for MettL3-MettL14.

Low levels of N6mA have been reported in DNA from mouse (105, 106) , human (30, 107) , and human malignant brain tumor glioblastoma (108) , though other studies have failed to detect N6mA in mammalian genomes (109) (110) (111) . Two reports present evidence that N6mA in mammalian DNA can result from incorporation of RNA-derived ribo-N6mA that was converted to deoxy-N6mA via the nucleotide salvage pathway (112, 113) . Interestingly, the three class ␤ MTases MettL3, MettL4, and MettL14, from mouse and human, were recently assessed for the ability to generate N6mA on DNA (114, 115) . These three proteins were previously considered to be analogs of the MT-A70 subunit of the human mRNA N6mA MTase (116) .

Prompted by the observations noted above, revealing that characterized class ␤ MTases are active for DNA N6mA methylation, with some members having activity on ssDNA (CcrM) as well as ssRNA (M.EcoGII), we investigated whether the MettL3-MettL14 heterodimer also possesses methyl transfer activity onto DNA adenine. With synthetic substrate oligonucleotides, MettL3-MettL14 shows >10fold stronger catalytic efficiency of methylation on ssDNA than the corresponding ssRNA under the same conditions (114) . Furthermore, MettL3-MettL14 is active on an unpaired region in overall context of dsDNA. This appears to be consistent with the requirement of MettL3 methylation activity for DNA repair (117) .

The third member of MT-A70 class ␤ MTase in human and mouse is MettL4 (116) . The C. elegans homolog was the first of the MettL4 family members to be considered as a candidate DNA N6mA MTase, but there was no enzymatic evidence demonstrating that C. elegans MettL4 has intrinsic DNA N6mA activity (118) . Murine MettL4 was reported to be responsible for N6mA deposition in genic elements associated with transcriptional silencing (106) . Curiously, recombinant human MettL4 expressed in HEK293T (human embryonic kidney) cells has in vitro enzymatic activity on mitochondrial DNA (115) , whereas recombinant human MettL4 purified from E. coli has RNA MTase activity (119) . The former study showed that MettL4 localizes within mitochondria in tested tissues (115) , and exam-ination of 23 high-confidence N6mA-enriched sites of mitochondrial DNA revealed an ApT containing consensus sequence--CTTATC (in the main text) or CT(C/A)ATC (in the figure S2E of (115))--agrees at least partially with an earlier study that N6mA sites across the mitochondrial genome were generally at an ApT dinucleotide (107) . In contrast, the latter study found mainly nuclear localization of MettL4 expressed from a lentiviral vector, and failed to identify appreciable levels of N6mA in mitochondrial DNA, but did find N6mA in a small nuclear RNA (the spliceosome-associated U2 snRNA) (119). The differences among these various studies illustrates the complex nature of DNA versus RNA adenine methylation in mammalian genomes (genomic vs. mitochondrial). The accumulation of N6mA in DNA and/or RNA might also reflect diverse cellular and mitochondrial stress responses in different cell lines under different laboratory conditions. Homology modeling suggested a dimeric MettL4, but was indeterminate as to whether it is a homodimer or heterodimer (the latter of which would require a yet-to-be-identified binding partner) ( Figure 4E) . A proper biochemical approach, such as the one used to identify the ciliate complex (93), will be needed to identify necessary components of the catalytically active MettL4 MTase complex.

Here, we build on earlier suggestions that an ancestral protein that bound a single molecule of SAM evolved into an adenine-binding protein; undergoing first tandem gene duplication, and then divergence of the two SAM-binding pockets ( Figure 1A ). M.EcoGII might reflect properties of the ancestral enzymes, as it acts on any form of nucleic acid (ss or ds, RNA or DNA), in any sequence context, though many nucleic acids-modifying enzymes are still able to modify both DNA and RNA ((120) and references therein). To become a sequence-specific DNA/RNA adenine MTase, an additional fusion (or selection acting on an exposed surface loop) would have brought in a target recognition domain (TRD) ( Figure 1B) . Distinguishing between gain of specificity from a promiscuous ancestral MTase, and loss of specificity from a more recent sequence-specific MTase, is not straightforward. This question has been addressed in some cases, for example among metabolic enzymes (121) , where increased specificity is seen in essential or higher-flux enzymes, or via reconstruction of ancestral enzymes from phylogenetic analyses, e.g. (122) (123) (124) (125) (126) . It might be revealing to carry out such analysis with M.EcoGII, or even to determine if it can methylate various adenine derivatives.

The unique dimeric feature of class ␤ MTases, either in homo-or hetero-form, allows the enzyme to use 'division of labor' between two subunits in terms of DNA recognition and methylation. This division could well have resulted from the development or importation of a TRD in an orientation opposite to that of the methylation domain of the same polypeptide, thus requiring the two domains from each subunit to face each other to form one integral functional surface. In contrast, the monomeric enzymes in class ␣ (E. coli Dam), class ␥ (M.TaqI) and 5mC class (M.HhaI) have TRDs in the same orientation as the methylation domain ( Figure 5 ). The opposing TRD orientation in the class ␤ MTases may help explain why they are primarily nucleic acid enzymes. Examining ∼30 human methyltransferaselike proteins (MettL1 to 27), all--except the three class ␤ MTases discussed here--have the seven-stranded catalytic domain arising from motif order of I-to-IV, and many of them evolved to be non-nucleic acid enzymes (127) . We wish to clarify a distinction between the 'division of labor' in homodimeric MTases, and non-␤ MTases that participate in protein-protein interactions for other purposes. For example, Trm112, named for its role in tRNA methylation, is a relatively small protein conserved in all three domains of life (128) . In S. cerevisiae, Trm112 interacts with and activates at least four MTases (Bud23, Trm9, Trm11 and Mtq2--all are non-␤ MTases) that target different components of the translation machinery (rRNA, tRNA and release factors) (129) .

To summarize, the class ␤ MTases, defined by their unique order of conserved motifs, may have diverged from an ancestral tandemly-duplicated protein that bound two molecules of SAM, differing from the other MTase classes either by which SAM-binding pocket changed to binding adenine, or via circular permutation, with M.EcoGII possibly representing an ancestral ␤ MTase. While these elements are somewhat speculative, there are clear consequences of the ␤ MTase structure on their function--specifically, that the ␤ MTases are likely to be catalytically active only as dimers, while other MTase classes may or may not form dimers, and that this is an intrinsic property distinct from other protein-protein interactions in which MTases may engage.

Predictive motifs derived from cytosine methyltransferases

Cloning of the BssHII restriction-modification system in Escherichia coli: BssHII methyltransferase contains circularly permuted cytosine-5 methyltransferase motifs

Circularly permuted variants of two CG-specific prokaryotic DNA methyltransferases

Cloning and sequencing of a cDNA encoding DNA methyltransferase of mouse cells. The carboxyl-terminal domain of the mammalian enzymes is related to bacterial restriction methyltransferases

Cloning and characterization of a family of novel mammalian DNA (cytosine-5) methyltransferases

Crystal structure of the HhaI DNA methyltransferase complexed with S-adenosyl-L-methionine

HhaI methyltransferase flips its target base out of the DNA helix

A candidate mammalian DNA methyltransferase related to pmt1p of fission yeast

Dnmt2 is not required for de novo and maintenance methylation of viral DNA in embryonic stem cells

Structure of human DNMT2, an enigmatic DNA methyltransferase homolog that displays denaturant-resistant binding to DNA

Methylation of tRNAAsp by the DNA methyltransferase homolog Dnmt2

The RNA methyltransferase Dnmt2 methylates DNA in the structural context of a tRNA

Biosynthesis and function of modified bases in bacteria and their viruses

Base flipping

AdoMet-dependent methylation, DNA methyltransferases and base flipping

DNA base Flipping: A general mechanism for writing, reading, and erasing DNA modifications

Comparison of protein structures reveals monophyletic origin of the AdoMet-dependent methyltransferase family and mechanistic convergence rather than recent differentiation of N4-cytosine and N6-adenine DNA methylation

Molecular evolution of DNA-(cytosine-N4) methyltransferases: evidence for their polyphyletic origin

The cytosine N4-methyltransferase M.PvuII also modifies adenine residues

Structure-guided analysis reveals nine sequence motifs conserved among DNA amino-methyltransferases, and suggests a catalytic mechanism for these enzymes

Three-dimensional structure of the adenine-specific DNA methyltransferase M.Taq I in complex with the cofactor S-adenosylmethionine

Structure of the N6-adenine DNA methyltransferase M.TaqI in complex with DNA and a cofactor analog

Universal catalytic domain structure of AdoMet-dependent methyltransferases

Kinetic and catalytic mechanism of HhaI methyltransferase

Direct identification of the active-site nucleophile in a DNA (cytosine-5)-methyltransferase

Mechanism of DNA methylation: the double role of DNA as a substrate and as a cofactor

Human DNMT1 transition state structure

Theoretical study of the catalytic mechanism of DNA-(N4-cytosine)-methyltransferase from the bacterium Proteus vulgaris

Is the HemK family of putative S-adenosylmethionine-dependent methyltransferases a 'missing' zeta subfamily of adenine methyltransferases? A hypothesis

N(6)-methyladenine DNA modification in the human genome

HemK, a class of protein methyl transferase with similarity to DNA methyl transferases, methylates polypeptide chain release factors, and hemK knockout induces defects in translational termination

The hemK gene in Escherichia coli encodes the N(5)-glutamine methyltransferase that modifies peptide release factors

Molecular basis for bacterial class I release factor methylation by PrmC

HemK2 protein, encoded on human chromosome 21, methylates translation termination factor eRF1

Human HemK2/KMT9/N6AMT1 is an active protein methyltransferase, but does not act on DNA in vitro, in the presence of Trm112

Structural insight into human N6AMT1-Trm112 complex functioning as a protein methyltransferase

KMT9 monomethylates histone H4 lysine 12 and controls proliferation of prostate cancer cells

Circular permutations in the molecular evolution of DNA methyltransferases

Evolution of type II DNA methyltransferases. A gene duplication model

The PvuII DNA (cytosine-N4)-methyltransferase comprises two trypsin-defined domains, each of which binds a molecule of S-adenosyl-L-methionine

Sequence permutations in the molecular evolution of DNA methyltransferases

Coenzymes as fossils of an earlier metabolic state

Modified nucleotides may have enhanced early RNA catalysis

REBASE-a database for DNA restriction and modification: enzymes, genes and genomes

Bacterial DNA methylation and methylomes

Highlights of the DNA cutters: a short history of the restriction enzymes

Effects of mutations in phage restriction sites during escape from restriction-modification

Diverse functions of restriction-modification systems in addition to cellular defense. Microbiol

Regulation of genetic flux between bacteria by restriction-modification systems

Type II restriction endonucleases-a historical perspective and more

Roles of DNA adenine methylation in host-pathogen interactions: mismatch repair, transcriptional regulation, and more

Timing and targeting: the biological functions of Dam methylation in E. coli

A cell cycle-regulated adenine DNA methyltransferase from Caulobacter crescentus processively methylates GANTC sites on hemimethylated DNA

Specificity versus Processivity in the Sequential Modification of DNA: A Study of DNA Adenine Methyltransferase

The Caulobacter crescentus DNA-(adenine-N6)-methyltransferase CcrM methylates DNA in a distributive manner

Caulobacter crescentus cell cycle-regulated DNA methyltransferase uses a novel mechanism for substrate recognition

Transition from nonspecific to specific DNA interactions along the substrate-recognition pathway of dam methyltransferase

Structure and substrate recognition of the Escherichia coli DNA adenine methyltransferase

The cell cycle-regulated DNA adenine methyltransferase CcrM opens a bubble at its DNA recognition site

Structure of PvuII DNA-(cytosine N4) methyltransferase, an example of domain permutation and protein fold assignment

Structure of RsrI methyltransferase, a member of the N6-adenine beta class of DNA methyltransferases

Structures of liganded and unliganded RsrI N6-adenine DNA methyltransferase: a distinct orientation for active cofactor binding

Purification and characterization of the M.RsrI DNA methyltransferase from Escherichia coli

Dimerization of the bacterial RsrI N6-adenine DNA methyltransferase

Crystal structure of MboIIA methyltransferase

Crystal structure of a putative DNA methylase TTHA0409 from Thermus thermophilus HB8

Biochemical and structural characterization of a DNA N6-adenine methyltransferase from Helicobacter pylori

DNA-(N4-cytosine)-methyltransferase from Bacillus amyloliquefaciens: kinetic and substrate binding properties

Characterization of the LlaCI methyltransferase from Lactococcus lactis subsp. cremoris W15 provides new insights into the biology of type II restriction-modification systems

Kinetic and catalytic properties of dimeric KpnI DNA methyltransferase

Functional analysis of amino acid residues at the dimerisation interface of KpnI DNA methyltransferase

Kinetic and catalytic properties of M.HpyAXVII, a phase-variable DNA methyltransferase from Helicobacter pylori

Type III restriction-modification enzymes: a historical perspective

DNA recognition by the EcoP15I and EcoPI modification methyltransferases

Asymmetric DNA methylation by dimeric EcoP15I DNA methyltransferase

Functional analysis of an acid adaptive DNA adenine methyltransferase from Helicobacter pylori 26695

Type III restriction endonucleases are heterotrimeric: comprising one helicase-nuclease subunit and a dimeric methyltransferase that binds only one specific DNA

Structural basis of asymmetric DNA methylation and ATP-triggered long-range diffusion by EcoP15I

Structure of Dnmt3a bound to Dnmt3L suggests a model for de novo DNA methylation

DNMT3L facilitates DNA methylation partly by maintaining DNMT3A stability in mouse embryonic stem cells

Dissect the DNMT3A-and DNMT3B-mediated DNA Co-methylation through a covalent complex approach

Genome-wide mapping of methylated adenine residues in pathogenic Escherichia coli using single-molecule real-time sequencing

Genomic epidemiology of the Escherichia coli O104:H4 outbreaks in Europe

The non-specific adenine DNA methyltransferase M.EcoGII

Identification of the active oligomeric state of an essential adenine DNA methyltransferase from Caulobacter crescentus

Exact size and organization of DNA target-recognizing domains of multispecific DNA-(cytosine-C5)-methyltransferases

High plasticity of multispecific DNA methyltransferases in the region carrying DNA target recognizing enzyme modules

DNA binding properties in vivo and target recognition domain sequence alignment analyses of wild-type and mutant RsrI [N6-adenine] DNA methyltransferases

Relaxed specificity of prokaryotic DNA methyltransferases results in DNA site-specific modification of RNA/DNA heteroduplexes

Macronuclear DNA of the hypotrichous ciliate Oxytricha fallax

Comparative study of DNA methylation in three unicellular eucaryotes

Widespread adenine N6-methylation of active genes in fungi

Identification of a DNA N6-Adenine methyltransferase complex and its impact on chromatin organization

Deoxyribonucleic acid methylation and chromatin organization in Tetrahymena thermophila

Comparison of methylated sequences in messenger RNA and heterogeneous nuclear RNA from mouse L cells

Characterization and partial purification of mRNA N6-adenosine methyltransferase from HeLa cell nuclei. Internal mRNA methylation requires a multisubunit complex

Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferase

A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation

RNA modifications modulate gene expression during development

N (6)-methyladenosine of chromosome-associated regulatory RNA regulates chromatin state and transcription

Solution structure of the RNA recognition domain

N(6)-methyladenosine methyltransferase

Structural basis of N(6)-adenosine methylation by the METTL3-METTL14 complex

Structural basis for cooperative function of Mettl3 and Mettl14 methyltransferases

Structural insights into the molecular mechanism of the m(6)A writer complex. eLife

DNA methylation on N(6)-adenine in mammalian embryonic stem cells

An adversarial DNA N(6)-Methyladenine-Sensor network preserves polycomb silencing

Single-nucleotide-resolution sequencing of human N6-methyldeoxyadenosine reveals strand-asymmetric clusters associated with SSBP1 on the mitochondrial genome

N(6)-methyladenine DNA Modification in Glioblastoma

Undetectable levels of N6-methyl adenine in mouse DNA: Cloning and analysis of PRED28, a gene coding for a putative mammalian DNA adenine methyltransferase

Quantitative LC-MS provides no evidence for m(6) dA or m(4) dC in the genome of mouse embryonic stem cells and tissues

Metabolically generated stable Isotope-Labeled deoxynucleoside code for tracing DNA N(6)-methyladenine in human cells

The origin of genomic N(6)-methyl-deoxyadenosine in mammalian cells

N(6)-methyladenine is incorporated into mammalian genome by DNA polymerase

Human MettL3-MettL14 complex is a sequence-specific DNA adenine methyltransferase active on single-strand and unpaired DNA in vitro

N(6)-deoxyadenosine methylation in mammalian mitochondrial DNA

Structure prediction and phylogenetic analysis of a functionally diverse family of proteins homologous to the MT-A70 subunit of the human mRNA:m(6)A methyltransferase

RNA m(6)A methylation regulates the ultraviolet-induced DNA damage response

METTL4 is an snRNA m(6)Am methyltransferase that regulates RNA splicing

DNA and RNA Modification Enzymes: Structure, Mechanism, Function and Evolution

Network context and selection in the evolution to enzyme specificity

Resurrection of ancestral effector caspases identifies novel networks for evolution of substrate specificity

Contextualizing risk and protection: perceptions of Bedouin mothers from unrecognized villages in the Naqab

Reconstructed ancestral enzymes reveal that negative selection drove the evolution of substrate specificity in ADP-dependent kinases

Ancestral protein reconstruction yields insights into adaptive evolution of binding specificity in Solute-Binding proteins

Evolution of minimal specificity and promiscuity in steroid hormone receptors

Protein lysine methylation by seven-beta-strand methyltransferases

Evolutionary insights into Trm112-methyltransferase holoenzymes involved in translation between archaea and eukaryotes

Trm112, a protein activator of methyltransferases modifying actors of the eukaryotic translational apparatus

We thank members of the Cheng laboratory for discussion. Authors contributions: The thoughts discussed here in this survey were accumulated over the last two decades during the studies of M.PvuII, Dam, M.EcoGII, CcrM, Trm112-HemK2 and MettL3-14. R.M.B. and X.C. prepared the manuscript during COVID-19 shelter-in-place requirements. All were involved in discussion. 

Conflict of interest statement. None declared.