Purification, cloning and sequence analysis of RsrI DNA methyltransferase: lack of homology between two enzymes, RsrI and EcoRI, that methylate the same nucleotide in identical recognition sequences.

RsrI DNA methyltransferase (M-RsrI) from Rhodobacter sphaeroides has been purified to homogeneity, and its gene cloned and sequenced. This enzyme catalyzes methylation of the same central adenine residue in the duplex recognition sequence d(GAATTC) as does M-EcoRI. The reduced and denatured molecular weight of the RsrI methyltransferase (MTase) is 33,600 Da. A fragment of R. sphaeroides chromosomal DNA exhibited M.RsrI activity in E. coli and was used to sequence the rsrIM gene. The deduced amino acid sequence of M.RsrI shows partial homology to those of the type II adenine MTases HinfI and DpnA and N4-cytosine MTases BamHI and PvuII, and to the type III adenine MTases EcoP1 and EcoP15. In contrast to their corresponding isoschizomeric endonucleases, the deduced amino acid sequences of the RsrI and EcoRI MTases show very little homology. Either the EcoRI and RsrI restriction-modification systems assembled independently from closely related endonuclease and more distantly related MTase genes, or the MTase genes diverged more than their partner endonuclease genes. The rsrIM gene sequence has also been determined by Stephenson and Greene (Nucl. Acids Res. (1989) 17, this issue).     

How does a DNA interacting enzyme change its specificity during molecular evolution? A site-directed mutagenesis study at the DNA binding site of the DNA-(adenine-N6)-methyltransferase EcoRV

The EcoRV DNA-(adenine-N6)-methyltransferase (MTase) recognizes GATATC sequences and modifies the first adenine residue within this site. Parts of its DNA interface show high sequence homology to DNA MTases of the dam family which recognize and modify GATC sequences. A phylogenetic analysis of M.EcoRV and dam-MTases suggests that EcoRV arose in evolution from a primordial dam-MTase in agreement to the finding that M.EcoRV also methylates GATC sites albeit at a strongly reduced rate. GATCTC sites that deviate in only one position from the EcoRV sequence are preferred over general dam sites. We have investigated by site-directed mutagenesis the function of 17 conserved and nonconserved residues within three loops flanking the DNA binding cleft of M.EcoRV. M.EcoRV contacts the GATATC sequence with two highly cooperative recognition modules. The contacts to the GAT-part of the recognition sequence are formed by residues conserved between dam MTases and M.EcoRV. Mutations at these positions lead to an increase in the discrimination between GATATC and GATC substrates. Our data show that the change in sequence specificity from dam (GATC) to EcoRV (GATATC) was accompanied by the generation of a second recognition module that contacts the second half of the target sequence. The new DNA contacts are formed by residues from all three loops that are not conserved between M.EcoRV and dam MTases. Mutagenesis at important residues within this module leads to variants that show a decreased ability to recognize the TC-part of the GATATC sequence.     

DsaV methyltransferase and its isoschizomers contain a conserved segment that is similar to the segment in Hhai methyltransferase that is in contact with DNA bases.

The methyltransferase (MTase) in the DsaV restriction--modification system methylates within 5'-CCNGG sequences. We have cloned the gene for this MTase and determined its sequence. The predicted sequence of the MTase protein contains sequence motifs conserved among all cytosine-5 MTases and is most similar to other MTases that methylate CCNGG sequences, namely M.ScrFI and M.SsoII. All three MTases methylate the internal cytosine within their recognition sequence. The 'variable' region within the three enzymes that methylate CCNGG can be aligned with the sequences of two enzymes that methylate CCWGG sequences. Remarkably, two segments within this region contain significant similarity with the region of M.HhaI that is known to contact DNA bases. These alignments suggest that many cytosine-5 MTases are likely to interact with DNA using a similar structural framework.     

The DNA and S-adenosylmethionine-binding regions of EcoDam and related methyltransferases

Previous comparison of the amino acid sequences of the GATC-methylating Escherichia coli Dam methyltransferase (MTase) with those of other adenine MTases (M.EcoRV, M.DpnII and T4Dam) localized four conserved regions. Regions III and IV have similarities with many other MTases. The sequence DPPY (or NPPY) is always present in region IV. It was suggested to be the AdoMet binding site. Publication of the nucleotide and amino acid sequences of M.CviBIII, M.DpnA and MutH give further credence to this assignment: M.DpnA, which also methylates GATC, has strong similarities with regions III and IV; M.CviBIII, a cytosine methylase, has a characteristic NPPY sequence in region IV, and only limited resemblance in region III; MutH, the GATC-specific endonuclease in DNA mismatch repair, has significant similarities uniquely in region III. The presently available evidence suggests that region III is the GAT(C) binding site and region IV is the AdoMet binding site. This hypothesis is strengthened by recent genetic findings.     

DNA methyltransferases of the cyanobacterium Anabaena PCC 7120

From the characterization of enzyme activities and the analysis of genomic sequences, the complement of DNA methyltransferases (MTases) possessed by the cyanobacterium ANABAENA PCC 7120 has been deduced. ANABAENA has nine DNA MTases. Four are associated with Type II restriction enzymes (AVAI, AVAII, AVAIII and the newly recognized inactive AVAIV), and five are not. Of the latter, four may be classified as solitary MTases, those whose function lies outside of a restriction/modification system. The group is defined here based on biochemical and genetic characteristics. The four solitary MTases, DmtA/M.AVAVI, DmtB/M.AVAVII, DmtC/M. AVAVIII and DmtD/M.AVAIX, methylate at GATC, GGCC, CGATCG and rCCGGy, respectively. DmtB methylates cytosines at the N4 position, but its sequence is more similar to N6-adenine MTases than to cytosine-specific enzymes, indicating that it may have evolved from the former. The solitary MTases, appear to be of ancient origin within cyanobacteria, while the restriction MTases appear to have arrived by recent horizontal transfer as did five now inactive Type I restriction systems. One Mtase, M.AVAV, cannot reliably be classified as either a solitary or restriction MTase. It is structurally unusual and along with a few proteins of prokaryotic and eukaryotic origin defines a structural class of MTases distinct from all previously described.     

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

A genetic selection method, the P22 challenge-phage assay, was used to characterize DNA binding in vivo by the prokaryotic beta class [N:6-adenine] DNA methyltransferase M.RSR:I. M.RSR:I mutants with altered binding affinities in vivo were isolated. Unlike the wild-type enzyme, a catalytically compromised mutant, M.RSR:I (L72P), demonstrated site-specific DNA binding in vivo. The L72P mutation is located near the highly conserved catalytic motif IV, DPPY (residues 65-68). A double mutant, M.RSR:I (L72P/D173A), showed less binding in vivo than did M.RSR:I (L72P). Thus, introduction of the D173A mutation deleteriously affected DNA binding. D173 is located in the putative target recognition domain (TRD) of the enzyme. Sequence alignment analyses of several beta class MTases revealed a TRD sequence element that contains the D173 residue. Phylogenetic analysis suggested that divergence in the amino acid sequences of these methyltransferases correlated with differences in their DNA target recognition sequences. Furthermore, MTases of other classes (alpha and gamma) having the same DNA recognition sequence as the beta class MTases share related regions of amino acid sequences in their TRDs.     

DNA methyltransferase induced by PBCV-1 virus infection of a Chlorella-like green alga.

A DNA methyltransferase was isolated from a eucaryotic, Chlorella-like green alga infected with the virus PBCV-1. The enzyme recognized the sequence GATC and methylated deoxyadenosine solely in GATC sequences. Host DNA, which contains GATC sequences, but not PBCV-1 DNA, which contains GmATC sequences, was a good substrate for the enzyme in vitro. The DNA methyltransferase activity was first detected about 1 h after viral infection; PBCV-1 DNA synthesis and host DNA degradation also began at about this time. The appearance of the DNA methyltransferase activity required de novo protein synthesis, and the enzyme was probably virus encoded. Methylation of DNAs with the PBCV-1-induced methyltransferase conferred resistance of the DNAs to a PBCV-1-induced restriction endonuclease enzyme described previously (Y. Xia, D. E. Burbank, L. Uher, D. Rabussay, and J. L. Van Etten, Mol. Cell. Biol. 6:1430-1439). We propose that the PBCV-1-induced methyltransferase protects viral DNA from the PBCV-1-induced restriction endonuclease and is part of a virus-induced restriction and modification system in PBCV-1-infected Chlorella cells.     

Mycoplasma CG- and GATC-specific DNA methyltransferases selectively and efficiently methylate the host genome and alter the epigenetic landscape in human cells

Aberrant DNA methylation is frequently observed in disease, including many cancer types, yet the underlying mechanisms remain unclear. Because germline and somatic mutations in the genes that are responsible for DNA methylation are infrequent in malignancies, additional mechanisms must be considered. Mycoplasmas spp., including Mycoplasma hyorhinis, efficiently colonize human cells and may serve as a vehicle for delivery of enzymatically active microbial proteins into the intracellular milieu. Here, we performed, for the first time, genome-wide and individual gene mapping of methylation marks generated by the M. hyorhinis CG- and GATC-specific DNA cytosine methyltransferases (MTases) in human cells. Our results demonstrated that, upon expression in human cells, MTases readily translocated to the cell nucleus. In the nucleus, MTases selectively and efficiently methylated the host genome at the DNA sequence sites free from pre-existing endogenous methylation, including those in a variety of cancer-associated genes. We also established that mycoplasma is widespread in colorectal cancers, suggesting that either the infection contributed to malignancy onset or, alternatively, that tumors provide a favorable environment for mycoplasma growth. In the human genome, ∼11% of GATC sites overlap with CGs (e.g., CGAT^mCG); therefore, the methylated status of these sites can be perpetuated by human DNMT1. Based on these results, we now suggest that the GATC-specific methylation represents a novel type of infection-specific epigenetic mark that originates in human cells with a previous exposure to infection. Overall, our findings unveil an entirely new panorama of interactions between the human microbiome and epigenome with a potential impact in disease etiology.     

Modulation of Escherichia coli DNA methyltransferase activity by biologically derived GATC-flanking sequences

Escherichia coli DNA adenine methyltransferase (EcoDam) methylates the N-6 position of the adenine in the sequence 5'-GATC-3' and plays vital roles in gene regulation, mismatch repair, and DNA replication. It remains unclear how the small number of critical GATC sites involved in the regulation of replication and gene expression are differentially methylated, whereas the approximately 20,000 GATCs important for mismatch repair and dispersed throughout the genome are extensively methylated. Our prior work, limited to the pap regulon, showed that methylation efficiency is controlled by sequences immediately flanking the GATC sites. We extend these studies to include GATC sites involved in diverse gene regulatory and DNA replication pathways as well as sites previously shown to undergo differential in vivo methylation but whose function remains to be assigned. EcoDam shows no change in affinity with variations in flanking sequences derived from these sources, but methylation kinetics varied 12-fold. A-tracts immediately adjacent to the GATC site contribute significantly to these differences in methylation kinetics. Interestingly, only when the poly(A) is located 5' of the GATC are the changes in methylation kinetics revealed. Preferential methylation is obscured when two GATC sites are positioned on the same DNA molecule, unless both sites are surrounded by large amounts of nonspecific DNA. Thus, facilitated diffusion and sequences immediately flanking target sites contribute to higher order specificity for EcoDam; we suggest that the diverse biological roles of the enzyme are in part regulated by these two factors, which may be important for other enzymes that sequence-specifically modify DNA.     

Role of histidine residues in EcoP15I DNA methyltransferase activity as probed by chemical modification and site-directed mutagenesis

Towards understanding the catalytic mechanism of M.EcoP15I [EcoP15I MTase (DNA methyltransferase); an adenine methyltransferase], we investigated the role of histidine residues in catalysis. M.EcoP15I, when incubated with DEPC (diethyl pyrocarbonate), a histidine-specific reagent, shows a time- and concentration-dependent inactivation of methylation of DNA containing its recognition sequence of 5'-CAGCAG-3'. The loss of enzyme activity was accompanied by an increase in absorbance at 240 nm. A difference spectrum of modified versus native enzyme shows the formation of N-carbethoxyhistidine that is diminished by hydroxylamine. This, along with other experiments, strongly suggests that the inactivation of the enzyme by DEPC was specific for histidine residues. Substrate protection experiments show that pre-incubating the methylase with DNA was able to protect the enzyme from DEPC inactivation. Site-directed mutagenesis experiments in which the 15 histidine residues in the enzyme were replaced individually with alanine corroborated the chemical modification studies and established the importance of His-335 in the methylase activity. No gross structural differences were detected between the native and H335A mutant MTases, as evident from CD spectra, native PAGE pattern or on gel filtration chromatography. Replacement of histidine with alanine residue at position 335 results in a mutant enzyme that is catalytically inactive and binds to DNA more tightly than the wild-type enzyme. Thus we have shown in the present study, through a combination of chemical modification and site-directed mutagenesis experiments, that His-335 plays an essential role in DNA methylation catalysed by M.EcoP15I.     

Differential methylation kinetics of individual target site strands by T4Dam DNA methyltransferase

Prokaryote DNA methyltransferases (MTases) of the Dam family (including those of bacteriophages T2 and T4) catalyze methyl group transfer from S-adenosyl-L-methionine (AdoMet), producing S-adenosyl-L-homocysteine (AdoHcy) and methylated adenine residues in palindromic GATC sequences. Dam DNA MTases, as all site-specific enzymes interacting with polymeric DNA, require a mechanism of action that ensures a rapid search for specific targets for catalytic action, during both the initial and subsequent rounds of methylation. The results of pre-steady-state (reaction burst) and steady-state methylation analyses of individual targets permitted us to monitor the action of T4Dam, which has three degrees of freedom: sliding, reorientation and adaptation to the canonical GATC sequence. The salient results are as follows: (i) 40mer substrate duplexes containing two canonical GATC sites showed differential methylation of the potential targets, i.e., T4Dam exhibited a preference for one site/target, which may present the better 'kinetic trap' for the enzyme. (ii) Prior hemimethylation of the two sites made both targets equally capable of being methylated during the pre-steady-state reaction. (iii) Although capable of moving in either direction along double-stranded DNA, there are some restrictions on T4Dam reorientation/adaptation on 40mer duplexes.     

Cloning and expression of the BalI restriction-modification system.

BalI, a type II restriction-modification (R-M) system from the bacterium, Brevibacterium albidum, recognizes the DNA sequence 5'-TGGCCA-3'. We cloned the genes encoding the BalI restriction endonuclease and methyltransferase and expressed them in Escherichia coli. The two genes were aligned tail-to-tail and their termination codons overlapped. BalI restriction endonuclease and methyltransferase comprise 260 and 280 amino acids, respectively, and have molecular weights of 29 043 and 31 999 Da. The amino acid sequence of BalI methyltransferase is similar to that of other m6A MTases, although it has been categorized as a m5C methyltransferase. A high expression system for the BalI restriction endonuclease was constructed in E. coli for the production of large quantities of enzyme.     

The M.EcoRV DNA-(adenine N6)-methyltransferase uses DNA bending for recognition of an expanded EcoDam recognition site

The M.EcoRV DNA methyltransferase recognizes GATATC sites. It is related to EcoDam, which methylates GATC sites. The DNA binding domain of M.EcoRV is similar to that of EcoDam suggesting a similar mechanism of DNA recognition. We show that amino acid residue Lys11 of M.EcoRV is involved in recognition of Gua1 and Arg128 contacts the Gua in base pair 6. These residues correspond to Lys9 and Arg124 in EcoDam, which recognize the Gua residues in both strands of the Dam recognition sequence, indicating that M.EcoRV and EcoDam make similar contacts to outermost base pairs of their recognition sequences and M.EcoRV recognizes its target site as an expanded GATC site. In contrast to EcoDam, M.EcoRV considerably bends the DNA (59+/-4 degrees) suggesting indirect readout of the AT-rich inner sequence. Recognition of an expanded target site by DNA bending is a new principle for changing DNA recognition specificity of proteins during molecular evolution. R128A is inefficient in DNA bending and binding, whereas K11A bends DNA with relaxed sequence specificity. These results suggest a temporal order of the formation of protein-DNA contacts in which the Gua6-Arg128 contact forms early followed by DNA bending and, finally, the formation of the Lys11-Gua1 contact.     

Cloning and nucleotide sequence of the gene encoding the Ecal DNA methyltransferase.

The gene coding for the GGTNACC specific Ecal DNA methyltransferase (M.Ecal) has been cloned in E. coli from Enterobacter cloacae and its nucleotide sequence has been determined. The ecalM gene codes for a protein of 452 amino acids (Mr: 51,111). It was determined that M.Ecal is an adenine methyltransferase. M.Ecal shows limited amino acid sequence similarity to other adenine methyltransferases. A clone that expresses Ecal methyltransferase at high level was constructed.     

Homology modelling of the DNA 5mC methyltransferase M.BssHII. Is permutation of functional subdomains common to all subfamilies of DNA methyltransferases?

This work presents a full tertiary model of the M.BssHII methyltransferase (MTase) complexed with substrate DNA and cofactor S-adenosyl-L-methionine, built by homology modelling based on previously solved complete structures of DNA MTases M.HaeIII and M. HhaI. M.BssHII and the template proteins show high sequence similarity, which indicates that they are evolutionary related. However, they are topologically different: M.BssHII is a circularly permuted variant of template MTases (Xu et al. Nucleic Acids Res 1997;25:3991). The model developed in this work will be a good starting point and valuable help in designing mutagenesis experiments to better understand the biological function of methyltransferases and the process of domain swapping.     

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.     

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

KpnI DNA-(N6-adenine) methyltransferase (M.KpnI) recognises the sequence 5'-GGTACC-3' and transfers the methyl group from S-adenosyl-L-methionine (AdoMet) to the N6 position of the adenine residue in each strand. Earlier studies have shown that M.KpnI exists as a dimer in solution, unlike most other MTases. To address the importance of dimerisation for enzyme function, a three-dimensional model of M.KpnI was obtained based on protein fold-recognition analysis, using the crystal structures of M.RsrI and M.MboIIA as templates. Residues I146, I161 and Y167, the side chains of which are present in the putative dimerisation interface in the model, were targeted for site-directed mutagenesis. Methylation and in vitro restriction assays showed that the mutant MTases are catalytically inactive. Mutation at the I146 position resulted in complete disruption of the dimer. The replacement of I146 led to drastically reduced DNA and cofactor binding. Substitution of I161 resulted in weakening of the interaction between monomers, leading to both monomeric and dimeric species. Steady-state fluorescence measurements showed that the wild-type KpnI MTase induces structural distortion in bound DNA, while the mutant MTases do not. The results establish that monomeric MTase is catalytically inactive and that dimerisation is an essential event for M.KpnI to catalyse the methyl transfer reaction.     

Cloning of the Haemophilus influenzae Dam methyltransferase and analysis of its relationship to the Dam methyltransferase encoded by the HP1 phage

In this paper we report cloning and experimental characterization of the DNA adenine methyltransferase (dam) gene from Haemophilus influenzae and comparison of its product with the Dam protein from the lysogenic phage of H. influenzae, HP1. Molecular modeling of M.HinDam and M.HP1Dam was carried out, providing a framework for a comparative analysis of these enzymes and their close homologs in the structural context. Both proteins share the common fold and essential cofactor-binding and catalytic residues despite overall divergence. However, subtle but significant differences in the cofactor-binding pocket have been identified. Moreover, while M.HinDam seems to contact its target DNA sequence using a number of loops, most of them are missing from M.HP1Dam. Analysis of both MTases suggests that their catalytic activity was derived from a common ancestor, but similar sequence specificities arose by convergence.