Cloning of a mammalian DNA methyltransferase

Cloning and sequencing of cDNA clones has shown that mammalian DNA (cytosine-5)-methyltransferase comprises a 1000-amino acid (aa) N-terminal region of unknown function and a 570-aa C-terminal region that is clearly related to bacterial type-II cytosine restriction methyltransferases. These findings indicate that the mammalian enzyme contains at least two structural domains and suggest a common evolutionary origin for mammalian and prokaryotic DNA (cytosine-5)-methyltransferases.     

Two substrates are better than one: dual specificities for Dnmt2 methyltransferases

Dnmt2 enzymes have been widely conserved during evolution and contain all of the signature motifs of DNA (cytosine-5)-methyltransferases; however, the DNA methyltransferase activity of these proteins is comparatively weak and their biochemical and functional properties remain enigmatic. Recent evidence now shows that Dnmt2 has a novel tRNA methyltransferase activity, raising the possibility that the biological roles of these proteins might be broader than previously thought. This finding has important implications for understanding the evolutionary relationships among these enzymes.     

On the evolutionary origin of eukaryotic DNA methyltransferases and Dnmt2

The Dnmt2 enzymes show strong amino acid sequence similarity with eukaryotic and prokaryotic DNA-(cytosine C5)-methyltransferases. Yet, Dnmt2 enzymes from several species were shown to methylate tRNA-Asp and had been proposed that eukaryotic DNA methyltransferases evolved from a Dnmt2-like tRNA methyltransferase ancestor [Goll et al., 2006, Science, 311, 395-8]. It was the aim of this study to investigate if this hypothesis could be supported by evidence from sequence alignments. We present phylogenetic analyses based on sequence alignments of the methyltransferase catalytic domains of more than 2300 eukaryotic and prokaryotic DNA-(cytosine C5)-methyltransferases and analyzed the distribution of DNA methyltransferases in eukaryotic species. The Dnmt2 homologues were reliably identified by an additional conserved CFT motif next to motif IX. All DNA methyltransferases and Dnmt2 enzymes were clearly separated from other RNA-(cytosine-C5)-methyltransferases. Our sequence alignments and phylogenetic analyses indicate that the last universal eukaryotic ancestor contained at least one member of the Dnmt1, Dnmt2 and Dnmt3 families of enzymes and additional RNA methyltransferases. The similarity of Dnmt2 enzymes with DNA methyltransferases and absence of similarity with RNA methyltransferases combined with their strong RNA methylation activity suggest that the ancestor of Dnmt2 was a DNA methyltransferase and an early Dnmt2 enzyme changed its substrate preference to tRNA. There is no phylogenetic evidence that Dnmt2 was the precursor of eukaryotic Dnmts. Most likely, the eukaryotic Dnmt1 and Dnmt3 families of DNA methyltransferases had an independent origin in the prokaryotic DNA methyltransferase sequence space.     

Substitutions of a cysteine conserved among DNA cytosine methylases result in a variety of phenotypes.

The proposed mechanism for DNA (cytosine-5)-methyltransferases envisions a key role for a cysteine residue. It is expected to form a covalent link with carbon 6 of the target cytosine, activating the normally inactive carbon 5 for methyl transfer. There is a single conserved cysteine among all DNA (cytosine-5)-methyltransferases making it the candidate nucleophile. We have changed this cysteine to other amino acids for the EcoRII methylase; which methylates the second cytosine in the sequence 5'-CCWGG-3'. Mutants were tested for their methyl transferring ability and for their ability to form covalent complexes with DNA. The latter property was tested indirectly with the use of a genetic assay involving sensitivity of cells to 5-azacytidine. Replacement of the conserved cysteine with glycine, valine, tryptophan or serine led to an apparent loss of methyl transferring ability. Interestingly, cells carrying the mutant with serine did show sensitivity to 5-azacytidine, suggesting the ability to link to DNA. Unexpectedly, substitution of the cysteine with glycine results in the inhibition of cell growth and the mutant allele can be maintained in the cells only when it is poorly expressed. These results suggest that the conserved cysteine in the EcoRII methylase is essential for methylase action and it may play more than one role in it.     

Analysis of Xenopus dsRNA adenosine deaminase cDNAs reveals similarities to DNA methyltransferases.

We isolated two similar, but distinct, cDNA classes that encode Xenopus double-stranded RNA (dsRNA) adenosine deaminase. The longest, full-length open reading frame (ORF) predicts a 1,270-amino acid protein of 138,754 Da that is similar in size and about 50% identical to proteins encoded by mammalian cDNAs, yet larger than the 120-kDa protein purified from Xenopus eggs. Alignments of the Xenopus and mammalian ORFs show N-terminal heterogeneity, three conserved dsRNA binding motifs (dsRBMs), and strongly conserved carboxyl termini. Consistent with the observation of two cDNA classes, northern analyses of Xenopus oocyte poly A+ RNA show at least three mRNA species. Multiple nuclear polyadenylation hexamers and putative cytoplasmic polyadenylation elements were found in the 3'' UTRs of cDNAs corresponding to the largest mRNA. In vitro translation experiments show that the cDNAs encode active deaminases and that the entire N-terminus and first dsRBM are dispensable for deaminase activity. Importantly, an analysis of the C-termini of five known dsRNA adenosine deaminases, and two putative deaminases, reveals motifs that are strikingly similar to the conserved motifs of the DNA-(adenine-N6alpha)-aminomethyltransferases and the DNA-(cytosine-5)-methyltransferases.     

The origins and early evolution of DNA mismatch repair genes--multiple horizontal gene transfers and co-evolution

To understand the evolutionary process of the DNA mismatch repair system, we conducted systematic phylogenetic analysis of its key components, the bacterial MutS and MutL genes and their eukaryotic homologs. Based on genome-wide homolog searches, we identified three new MutS subfamilies (MutS3-5) in addition to the previously studied MutS1 and MutS2 subfamilies. Detailed evolutionary analysis strongly suggests that frequent ancient horizontal gene transfer (HGT) occurred with both MutS and MutL genes from bacteria to eukaryotes and/or archaea. Our results further imply that the origins of mismatch repair system in eukaryotes and archaea are largely attributed to ancient HGT from bacteria instead of vertical evolution. Specifically, the eukaryotic MutS and MutL homologs likely originated from endosymbiotic ancestors of mitochondria or chloroplasts, indicating that not only archaea, but also bacteria are important sources of eukaryotic DNA metabolic genes. The archaeal MutS1 and MutL homologs were also acquired from bacteria simultaneously through HGT. Moreover, the distribution and evolution profiles of the MutS1 and MutL genes suggest that they have undergone long-term coevolution. Our work presents an overall portrait of the evolution of these important genes in DNA metabolism and also provides further understanding about the early evolution of cellular organisms.     

Kinetic and catalytic mechanism of HhaI methyltransferase

Kinetic and catalytic properties of the DNA (cytosine-5)-methyltransferase HhaI are described. With poly(dG-dC) as substrate, the reaction proceeds by an equilibrium (or processive) ordered Bi-Bi mechanism in which DNA binds to the enzyme first, followed by S-adenosylmethionine (AdoMet). After methyl transfer, S-adenosylhomocysteine (AdoHcy) dissociates followed by methylated DNA. AdoHcy is a potent competitive inhibitor with respect to AdoMet (Ki = 2.0 microM) and its generation during reactions results in non-linear kinetics. AdoMet and AdoHcy significantly interact with only the substrate enzyme-DNA complex; they do not bind to free enzyme and bind poorly to the methylated enzyme-DNA complex. In the absence of AdoMet, HhaI methylase catalyzes exchange of the 5-H of substrate cytosines for protons of water at about 7-fold the rate of methylation. The 5-H exchange reaction is inhibited by AdoMet or AdoHcy. In the enzyme-DNA-AdoHcy complex, AdoHcy also suppresses dissociation of DNA and reassociation of the enzyme with other substrate sequences. Our studies reveal that the catalytic mechanism of DNA (cytosine-5)-methyltransferases involves attack of the C6 of substrate cytosines by an enzyme nucleophile and formation of a transient covalent adduct. Based on precedents of other enzymes which catalyze similar reactions and the susceptibility of HhaI to inactivation by N-ethylmaleimide, we propose that the sulfhydryl group of a cysteine residue is the nucleophilic catalyst. Furthermore, we propose that Cys-81 is the active-site catalyst in HhaI. This residue is found in a Pro-Cys doublet which is conserved in all DNA (cytosine-5)-methyltransferases whose sequences have been determined to date and is found in related enzymes. Finally, we discuss the possibility that covalent adducts between C6 of pyrimidines and nucleophiles of proteins may be important general components of protein-nucleic acid interactions.     

The eukaryotic DNMT2 genes encode a new class of cytosine-5 DNA methyltransferases

DNMT2 is a subgroup of the eukaryotic cytosine-5 DNA methyltransferase gene family. Unlike the other family members, proteins encoded by DNMT2 genes were not known before to possess DNA methyltransferase activities. Most recently, we have shown that the genome of Drosophila S2 cells stably expressing an exogenous Drosophila dDNMT2 cDNA became anomalously methylated at the 5'-positions of cytosines (Reddy, M. N., Tang, L. Y., Lee, T. L., and Shen, C.-K. J. (2003) Oncogene, in press). We present evidence here that the genomes of transgenic flies overexpressing the dDnmt2 protein also became hypermethylated at specific regions. Furthermore, transient transfection studies in combination with sodium bisulfite sequencing demonstrated that dDnmt2 as well as its mouse ortholog, mDnmt2, are capable of methylating a cotransfected plasmid DNA. These data provide solid evidence that the fly and mouse DNMT2 gene products are genuine cytosine-5 DNA methyltransferases.     

Exposition of a family of RNA m(5)C methyltransferases from searching genomic and proteomic sequences.

The Escherichia coli fmu gene product has recently been determined to be the 16S rRNA m(5)C 967 methyltransferase. As such, Fmu represents the first protein identified as an S -adenosyl-L-methionine (AdoMet)- dependent RNA m(5)C methyltransferase whose amino acid sequence is known. Using the amino acid sequence of Fmu as an initial probe in an iterative search of completed DNA sequence databases, 27 homologous ORF products were identified as probable RNA m(5)C methyltransferases. Further analysis of sequences in undeposited genomic sequencing data and EST databases yielded more than 30 additional homologs. These putative RNA m(5)C methyltransferases are grouped into eight subfamilies, some of which are predicted to consist of direct genetic counterparts, or orthologs. The enzymes proposed to be RNA m(5)C methyltransferases have sequence motifs closely related to signature sequences found in the well-studied DNA m(5)C methyltransferases and other AdoMet-dependent methyltransferases. Structure-function correlates in the known AdoMet methyltransferases support the assignment of this family as RNA m(5)C methyltransferases.     

Molecular phylogenetics of DNA 5mC-methyltransferases.

We have identified a total of 88 members of the DNA-(cytosine-5) methyltransferase (5mC MTase) family whose sequences have been deposited in the databases. The results of a comparison of these sequences is presented in the form of an alignment-based phylogenetic tree and sequence logos for 10 conserved motifs. Phylogenetic analysis showed that members of the family aggregate into subfamilies which are usually consistent with their target specificity. However, it was also shown that similar target specificity does not necessarily imply close homology of the catalytic domain of MTases, which strongly supports the hypothesis that target recognition evolved independently of catalytic properties. This analysis also indicate that the 5mC MTase was present in the cenancestor (last common ancestor) of eubacteria, archaebacteria, and eukaryotes. The phylogeny of the 5mC MTases catalytic domain provides the basis for establishing the patterns of evolutionary change that characterize this family of proteins with conserved structural core and variable and mobile modules not directly involved in formation of the active site.     

A mutant HpaII methyltransferase functions as a mutator enzyme.

DNA (cytosine-5)-methyltransferases can cause deamination of cytosine when the cofactor S-adenosylmethionine (AdoMet) is limiting and thus function as sequence-specific C-->U mutator enzymes. Here we explored whether mutations causing inactivation of the cofactor binding activity of the HpaII methyltransferase, thus mimicking conditions of limiting AdoMet concentration, could convert a DNA methyltransferase to a C-->U mutator enzyme. We created two mutator enzymes from the HpaII methyltransferase (F38S and G40D) which both showed enhanced cytosine deamination activities in vitro and in vivo. Interestingly, the G:U mispairs generated by these enzymes were not repaired completely in bacteria equipped with uracil-DNA glycosylase-initiated repair machinery, giving rise to a potent mutator phenotype. This is the first report showing the creation of mutator enzymes from a DNA methyltransferase and the demonstration of their mutagenicity in living cells.     

Transposition burst of mariner-like elements in the sequenced genome of Rhodnius prolixus

Transposable elements (TEs) are widespread in insect's genomes. However, there are wide differences in the proportion of the total DNA content occupied by these repetitive sequences in different species. We have analyzed the TEs present in R. prolixus (vector of the Chagas disease) and showed that 3.0% of this genome is occupied by Class II TEs, belonging mainly to the Tc1-mariner superfamily (1.65%) and MITEs (1.84%). Interestingly, most of this genomic content is due to the expansion of two subfamilies belonging to: irritans himar, a well characterized subfamily of mariners, and prolixus1, one of the two novel subfamilies here described. The high amount of sequences in these subfamilies suggests that bursts of transposition occurred during the life cycle of this family. In an attempt to characterize these elements, we performed an in silico analysis of the sequences corresponding to the DDD/E domain of the transposase gene. We performed an evolutionary analysis including network and Bayesian coalescent-based methods in order to infer the dynamics of the amplification, as well as to estimate the time of the bursts identified in these subfamilies. Given our data, we hypothesized that the TE expansions occurred around the time of speciation of R. prolixus around 1.4 mya. This suggestion lays on the “Transposon Model” of TE evolution, in which the members of a TE population that are replicative active are present at multiple loci in the genome, but their replicative potential varies, and of the “Life Cycle Model” that states that when present-day TEs have been involved in amplification bursts, they share an ancestral copy that dates back to this initial amplification.     

Biochemical and mutational analysis of EcoRII functional domains reveals evolutionary links between restriction enzymes

The archetypal Type IIE restriction endonuclease EcoRII is a dimer that has a modular structure. DNA binding studies indicate that the isolated C-terminal domain dimer has an interface that binds a single cognate DNA molecule whereas the N-terminal domain is a monomer that also binds a single copy of cognate DNA. Hence, the full-length EcoRII contains three putative DNA binding interfaces: one at the C-terminal domain dimer and two at each of the N-terminal domains. Mutational analysis indicates that the C-terminal domain shares conserved active site architecture and DNA binding elements with the tetrameric restriction enzyme NgoMIV. Data provided here suggest possible evolutionary relationships between different subfamilies of restriction enzymes.     

Methylation inhibitors can increase the rate of cytosine deamination by (cytosine-5)-DNA methyltransferase.

The target cytosines of (cytosine-5)-DNA methyltransferases in prokaryotic and eukaryotic DNA show increased rates of C-->T transition mutations compared to non-target cytosines. These mutations are induced either by the spontaneous deamination of 5-mC-->T generating inefficiently repaired G:T rather than G:U mismatches, or by the enzyme-induced C-->U deamination which occurs under conditions of reduced levels of S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy). We tested whether various inhibitors of (cytosine-5)-DNA methyltransferases analogous to AdoMet and AdoHcy would affect the rate of enzyme-induced deamination of the target cytosine by M.HpaII and M.SssI. Interestingly, we found two compounds, sinefungin and 5'-amino-5'-deoxyadenosine, that increased the rate of deamination 10(3)-fold in the presence and 10(4)-fold in the absence of AdoMet and AdoHcy. We have therefore identified the first mutagenic compounds specific for the target sites of (cytosine-5)-DNA methyltransferases. A number of analogs of AdoMet and AdoHcy have been considered as possible antiviral, anticancer, antifungal and antiparasitic agents. Our findings show that chemotherapeutic agents with affinities to the cofactor binding pocket of (cytosine-5)-DNA methyltransferase should be tested for their potential mutagenic effects.     

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

MT-A70 is the S-adenosylmethionine-binding subunit of human mRNA:m(6)A methyl-transferase (MTase), an enzyme that sequence-specifically methylates adenines in pre-mRNAs. The physiological importance yet limited understanding of MT-A70 and its apparent lack of similarity to other known RNA MTases combined to make this protein an attractive target for bioinformatic analysis. The sequence of MT-A70 was subjected to extensive in silico analysis to identify orthologous and paralogous polypeptides. This analysis revealed that the MT-A70 family comprises four subfamilies with varying degrees of interrelatedness. One subfamily is a small group of bacterial DNA:m(6)A MTases. The other three subfamilies are paralogous eukaryotic lineages, two of which have not been associated with MTase activity but include proteins having substantial regulatory effects. Multiple sequence alignments and structure prediction for members of all four subfamilies indicated a high probability that a consensus MTase fold domain is present. Significantly, this consensus fold shows the permuted topology characteristic of the b class of MTases, which to date has only been known to include DNA MTases.