On the substrate specificity of DNA methyltransferases. adenine-N6 DNA methyltransferases also modify cytosine residues at position N4

Methylation of DNA is important in many organisms and essential in mammals. Nucleobases can be methylated at the adenine-N6, cytosine-N4, or cytosine-C5 atoms by specific DNA methyltransferases. We show here that the M.EcoRV, M.EcoRI, and Escherichia coli dam methyltransferases as well as the N- and C-terminal domains of the M. FokI enzyme, which were formerly all classified as adenine-N6 DNA methyltransferases, also methylate cytosine residues at position N4. Kinetic analyses demonstrate that the rate of methylation of cytosine residues by M.EcoRV and the M.FokI enzymes is reduced by only 1-2 orders of magnitude in relation to methylation of adenines. This result shows that although these enzymes methylate DNA in a sequence specific manner, they have a low substrate specificity with respect to the target base. This unexpected finding has implications on the mechanism of adenine-N6 DNA methyltransferases. Sequence comparisons suggest that adenine-N6 and cytosine-N4 methyltransferases have changed their reaction specificity at least twice during evolution, a model that becomes much more likely given the partial functional overlap of both enzyme types. In contrast, methylation of adenine residues by the cytosine-N4 methyltransferase M.BamHI was not detectable. On the basis of our results, we suggest that adenine-N6 and cytosine-N4 methyltransferases should be grouped into one enzyme family.     

Functional analysis of BamHI DNA cytosine-N4 methyltransferase

We show that the kinetic mechanism of the DNA (cytosine-N(4)-)-methyltransferase M.BamHI, which modifies the underlined cytosine (GGATCC), differs from cytosine C(5) methyltransferases, and is similar to that observed with adenine N(6) methyltransferases. This suggests that the obligate order of ternary complex assembly and disassembly depends on the type of methylation reaction. In contrast, the single-turnover rate of catalysis for M.BamHI (0.10s(-1)) is closer to the DNA (cytosine-C(5)-)-methyltransferases (0.14s(-1)) than the DNA (adenine-N(6)-)-methyltransferases (>200s(-1)). The nucleotide flipping transition dominates the single-turnover constant for adenine N(6) methyltransferases, and, since the disruption of the guanine-cytosine base-pair is essential for both types of cytosine DNA methyltransferases, this transition may be a common, rate-limiting step for methylation for these two enzyme subclasses. The similar overall rate of catalysis by M.BamHI and other DNA methyltransferases is consistent with a common rate-limiting catalytic step of product dissociation. Our analyses of M.BamHI provide functional insights into the relationship between the three different classes of DNA methyltransferases that complement both prior structural and evolutionary insights.     

Chemical display of thymine residues flipped out by DNA methyltransferases.

The DNA cytosine-C5 methyltransferase M. Hha I flips its target base out of the DNA helix during interaction with the substrate sequence GCGC. Binary and ternary complexes between M. Hha I and hemimethylated DNA duplexes were used to examine the suitability of four chemical methods to detect flipped-out bases in protein-DNA complexes. These methods probe the structural peculiarities of pyrimidine bases in DNA. We find that in cases when the target cytosine is replaced with thymine (GTGC), KMnO4proved an efficient probe for positive display of flipped-out thymines. The generality of this procedure was further verified by examining a DNA adenine-N6 methyltransferase, M. Taq I, in which case an enhanced reactivity of thymine replacing the target adenine (TCGT) in the recognition sequence TCGA was also observed. Our results support the proposed base-flipping mechanism for adenine methyltransferases, and offer a convenient laboratory tool for detection of flipped-out thymines in protein-DNA complexes.     

Symmetry elements in DNA structure important for recognition/methylation by DNA [amino]-methyltransferases

The phage T4Dam and EcoDam DNA-[adenine-N6] methyltransferases (MTases) methylate GATC palindromic sequences, while the BamHI DNA-[cytosine-N4] MTase methylates the GGATCC palindrome (which contains GATC) at the internal cytosine residue. We compared the ability of these enzymes to interact productively with defective duplexes in which individual elements were deleted on one chain. A sharp decrease in k_cat was observed for all three enzymes if a particular element of structural symmetry was disrupted. For the BamHI MTase, integrity of the ATCC was critical, while an intact GAT sequence was necessary for the activity of T4Dam, and an intact GA was necessary for EcoDam. Theoretical alignment of the region of best contacts between the protein and DNA showed that in the case of a palindromic interaction site, a zone covering the 5′-symmetric residues is located in the major groove versus a zone of contact covering the 3′-symmetric residues in the minor groove. Our data fit a simple rule of thumb that the most important contacts are aligned around the methylation target base: if the target base is in the 5′ half of the palindrome, the interaction between the enzyme and the DNA occurs mainly in the major groove; if it is in the 3′ half, the interaction occurs mainly in the minor groove.     

Identification of the binding site for the extrahelical target base in N6-adenine DNA methyltransferases by photo-cross-linking with duplex oligodeoxyribonucleotides containing 5-iodouracil at the target position.

DNA methyltransferases flip their target bases out of the DNA double helix for catalysis. Base flipping of C5-cytosine DNA methyltransferases was directly observed in the protein-DNA cocrystal structures of M.HhaI and M.HaeIII. Indirect structural evidence for base flipping of N6-adenine and N4-cytosine DNA methyltransferases was obtained by modeling DNA into the three-dimensional structures of M.TaqI and M.PvuII in complex with the cofactor. In addition, biochemical evidence of base flipping was reported for different N6-adenine DNA methyltransferases. As no protein-DNA cocrystal structure for the related N6-adenine and N4-cytosine DNA methyltransferases is available, we used light-induced photochemical cross-linking to identify the binding site of the extrahelical target bases. The N6-adenine DNA methyltransferases M.TaqI and M.CviBIII, which both methylate adenine within the double-stranded 5'-TCGA-3' DNA sequence, were photo-cross-linked to duplex oligodeoxyribonucleotides containing 5-iodouracil at the target position in 50-60% and almost quantitative yield, respectively. Proteolytic fragmentation of the M. CviBIII-DNA complex followed by Edman degradation and electrospray ionization mass spectrometry indicates photo-cross-linking to tyrosine 122. In addition, the mutant methyltransferases M. TaqI/Y108A and M.TaqI/F196A were photo-cross-linked with 6-fold and 2-fold reduced efficiency, respectively, which suggests that tyrosine 108 is the primary site of modification in M.TaqI. Our results indicate a close proximity between the extrahelical target base and tyrosine 122 in M.CviBIII or tyrosine 108 in M.TaqI. As both residues belong to the conserved motif IV ((N/D/S)(P/I)P(Y/F/W)) found in all N6-adenine and N4-cytosine DNA as well as in N6-adenine RNA methyltransferases, a similar spatial relationship between the target bases and the aromatic amino acid residue within motif IV is expected for all these methyltransferases.     

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

We have determined the structure of Pvu II methyltransferase (M. Pvu II) complexed with S -adenosyl-L-methionine (AdoMet) by multiwavelength anomalous diffraction, using a crystal of the selenomethionine-substituted protein. M. Pvu II catalyzes transfer of the methyl group from AdoMet to the exocyclic amino (N4) nitrogen of the central cytosine in its recognition sequence 5'-CAGCTG-3'. The protein is dominated by an open alpha/beta-sheet structure with a prominent V-shaped cleft: AdoMet and catalytic amino acids are located at the bottom of this cleft. The size and the basic nature of the cleft are consistent with duplex DNA binding. The target (methylatable) cytosine, if flipped out of the double helical DNA as seen for DNA methyltransferases that generate 5-methylcytosine, would fit into the concave active site next to the AdoMet. This M. Pvu IIalpha/beta-sheet structure is very similar to those of M. Hha I (a cytosine C5 methyltransferase) and M. Taq I (an adenine N6 methyltransferase), consistent with a model predicting that DNA methyltransferases share a common structural fold while having the major functional regions permuted into three distinct linear orders. The main feature of the common fold is a seven-stranded beta-sheet (6 7 5 4 1 2 3) formed by five parallel beta-strands and an antiparallel beta-hairpin. The beta-sheet is flanked by six parallel alpha-helices, three on each side. The AdoMet binding site is located at the C-terminal ends of strands beta1 and beta2 and the active site is at the C-terminal ends of strands beta4 and beta5 and the N-terminal end of strand beta7. The AdoMet-protein interactions are almost identical among M. Pvu II, M. Hha I and M. Taq I, as well as in an RNA methyltransferase and at least one small molecule methyltransferase. The structural similarity among the active sites of M. Pvu II, M. Taq I and M. Hha I reveals that catalytic amino acids essential for cytosine N4 and adenine N6 methylation coincide spatially with those for cytosine C5 methylation, suggesting a mechanism for amino methylation.     

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

DNA methylation is important in cellular, developmental and disease processes, as well as in bacterial restriction-modification systems. Methylation of DNA at the amino groups of cytosine and adenine is a common mode of protection against restriction endonucleases afforded by the bacterial methyltransferases. The first structure of an N:6-adenine methyltransferase belonging to the beta class of bacterial methyltransferases is described here. The structure of M. RSR:I from Rhodobacter sphaeroides, which methylates the second adenine of the GAATTC sequence, was determined to 1.75 A resolution using X-ray crystallography. Like other methyltransferases, the enzyme contains the methylase fold and has well-defined substrate binding pockets. The catalytic core most closely resembles the PVU:II methyltransferase, a cytosine amino methyltransferase of the same beta group. The larger nucleotide binding pocket observed in M. RSR:I is expected because it methylates adenine. However, the most striking difference between the RSR:I methyltransferase and the other bacterial enzymes is the structure of the putative DNA target recognition domain, which is formed in part by two helices on an extended arm of the protein on the face of the enzyme opposite the active site. This observation suggests that a dramatic conformational change or oligomerization may take place during DNA binding and methylation.     

M.HhaI binds tightly to substrates containing mismatches at the target base.

The (cytosine-5) DNA methyltransferase M.HhaI causes its target cytosine base to be flipped completely out of the DNA helix upon binding. We have investigated the effects of replacing the target cytosine by other, mismatched bases, including adenine, guanine, thymine and uracil. We find that M.HhaI binds more tightly to such mismatched substrates and can even transfer a methyl group to uracil if a G:U mismatch is present. Other mismatched substrates in which the orphan guanine is changed exhibit similar behavior. Overall, the affinity of DNA binding correlates inversely with the stability of the target base pair, while the nature of the target base appears irrelevant for complex formation. The presence of a cofactor analog. S-adenosyl-L-homocysteine, greatly enhances the selectivity of the methyltransferase for cytosine at the target site. We propose that the DNA methyltransferases have evolved from mismatch binding proteins and that base flipping was, and still is, a key element in many DNA-enzyme interactions.     

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

The 2.0 A crystal structure of the N6-adenine DNA methyltransferase M.TaqI in complex with specific DNA and a nonreactive cofactor analog reveals a previously unrecognized stabilization of the extrahelical target base. To catalyze the transfer of the methyl group from the cofactor S-adenosyl-l-methionine to the 6-amino group of adenine within the double-stranded DNA sequence 5'-TCGA-3', the target nucleoside is rotated out of the DNA helix. Stabilization of the extrahelical conformation is achieved by DNA compression perpendicular to the DNA helix axis at the target base pair position and relocation of the partner base thymine in an interstrand pi-stacked position, where it would sterically overlap with an innerhelical target adenine. The extrahelical target adenine is specifically recognized in the active site, and the 6-amino group of adenine donates two hydrogen bonds to Asn 105 and Pro 106, which both belong to the conserved catalytic motif IV of N6-adenine DNA methyltransferases. These hydrogen bonds appear to increase the partial negative charge of the N6 atom of adenine and activate it for direct nucleophilic attack on the methyl group of the cofactor.     

Sequence motifs characteristic of DNA[cytosine-N4]methyltransferases: similarity to adenine and cytosine-C5 DNA-methylases

The sequences coding for DNA[cytosine-N4]methyltransferases MvaI (from Micrococcus varians RFL19) and Cfr9I (from Citrobacter freundii RFL9) have been determined. The predicted methylases are proteins of 454 and 300 amino acids, respectively. Primary structure comparison of M.Cfr9I and another m4C-forming methylase, M.Pvu II, revealed extended regions of homology. The sequence comparison of the three DNA[cytosine-N4]-methylases using originally developed software revealed two conserved patterns, DPF-GSGT and TSPPY, which were found similar also to those of adenine and DNA[cytosine-C5]-methylases. These data provided a basis for global alignment and classification of DNA-methylase sequences. Structural considerations led us to suggest that the first region could be the binding site of AdoMet, while the second is thought to be directly involved in the modification of the exocyclic amino group.     

Changing the target base specificity of the EcoRV DNA methyltransferase by rational de novo protein-design

The EcoRV DNA-(adenine-N⁶)-methyltransferase (M.EcoRV) specifically modifies the first adenine residue within GATATC sequences. During catalysis, the enzyme flips its target base out of the DNA helix and binds it into a target base binding pocket which is formed in part by Lys16 and Tyr196. A cytosine residue is accepted by wild-type M.EcoRV as a substrate at a 31-fold reduced efficiency with respect to the k_cat/K_M values if it is located in a CT mismatch substrate (GCTATC/GATATC). Cytosine residues positioned in a CG base pair (GCTATC/GATAGC) are modified at much more reduced rates, because flipping out the target base is much more difficult in this case. We intended to change the target base specificity of M.EcoRV from adenine-N⁶ to cytosine-N⁴. To this end we generated, purified and characterized 15 variants of the enzyme, containing single, double and triple amino acid exchanges following different design approaches. One concept was to reduce the size of the target base binding pocket by site-directed mutagenesis. The K16R variant showed an altered specificity, with a 22-fold preference for cytosine as the target base in a mismatch substrate. This corresponds to a 680-fold change in specificity, which was accompanied by only a small loss in catalytic activity with the cytosine substrate. The K16R/Y196W variant no longer methylated adenine residues at all and its activity towards cytosine was reduced only 17-fold. Therefore, we have changed the target base specificity of M.EcoRV from adenine to cytosine by rational protein design. Because there are no natural paragons for the variants described here, a change of the target base specificity of a DNA interacting enzyme was possible by rational de novo design of its active site.     

Chemical mapping of cytosines enzymatically flipped out of the DNA helix

Haloacetaldehydes can be employed for probing unpaired DNA structures involving cytosine and adenine residues. Using an enzyme that was structurally proven to flip its target cytosine out of the DNA helix, the HhaI DNA methyltransferase (M.HhaI), we demonstrate the suitability of the chloroacetaldehyde modification for mapping extrahelical (flipped-out) cytosine bases in protein-DNA complexes. The generality of this method was verified with two other DNA cytosine-5 methyltransferases, M.AluI and M.SssI, as well as with two restriction endonucleases, R.Ecl18kI and R.PspGI, which represent a novel class of base-flipping enzymes. Our results thus offer a simple and convenient laboratory tool for detection and mapping of flipped-out cytosines in protein-DNA complexes.     

Probing the DNA interface of the EcoRV DNA-(adenine-N6)-methyltransferase by site-directed mutagenesis,fluorescence spectroscopy,and UV cross-linking

The EcoRV DNA-(adenine-N6)-methyltransferase recognizes GATATC sites and methylates the DNA as indicated. It is related to the large family of dam methyltransferases which modify GATC sites. We have studied the interaction of DNA with M.EcoRV and 12 M.EcoRV variants using oligonucleotides containing 2-aminopurine as a fluorescence probe in equilibrium and stopped-flow DNA binding studies and 5-iododeoxyuracil for UV cross-linking. M.EcoRV binds to DNA in a multistep binding reaction, including two different conformations of the specific enzyme-DNA complex, and induces a strong conformational change of the DNA at the fourth position of the recognition site. Mutations at residues forming contacts to the GAT part of the recognition site reduce the stability of both specific enzyme-DNA complexes. Two enzyme variants which fail to recognize the ATC part do not induce the deformation of the DNA which explains why they cannot interact properly with the recognition site. Other mutations at residues which interact with the ATC part selectively reduce the stability of the second enzyme-DNA complex. These results show that when approaching the DNA M.EcoRV first contacts the GAT part of the target site. Since the residues mediating these contacts are conserved among M.EcoRV and dam MTases, the kinetics of formation of the enzyme-DNA complex correspond to the evolutionary history of the protein. Whether the observation that evolutionarily conserved contacts are formed early during complex formation is a general rule for DNA interacting enzymes or proteins that change their specificity during evolution remains to be seen.     

Sequence, internal homology and high-level expression of the gene for a DNA-(cytosine N4)-methyltransferase, M.Pvu II.

The base sequence of the pvuIIM gene has been determined. This gene codes for a DNA-(cytosine N4)-methyltransferase, M.Pvu II. The base sequence contains a single large open reading frame that predicts a 38.3kDa polypeptide, consistent with experimental data. The pvuIIM gene contains some sequences common to DNA methyltransferases in general, but includes none of the sequences specifically conserved among DNA-(cytosine 5)-methyltransferases. The pvuIIM sequence also reveals an internal homology at the amino acid level, each half of which spans over 100 amino acids and is itself homologous to the sequences of some DNA-(adenine N6)-methyltransferases. A derivative of the pvuIIM plasmid was constructed to allow high-level production of M.Pvu II. Specifically, the composite Ptac promoter was inserted 5' to pvuIIM, intervening DNA was deleted, and the resulting construct was used to transform an mcrB laclq strain of Escherichia coli. When this transformant was induced with isopropyl-B-D-galactopyranoside (IPTG), growth rapidly ceased and M.Pvu II accumulated to the point of comprising over 10% of the total soluble protein.     

Isolation and Characterization of Site-Specific DNA-methyltransferases from Bacillus coagulans K

Two site-specific DNA methyltransferases, M.BcoKIA and M.BcoKIB, were isolated from the thermophilic strain Bacillus coagulans K. Each of the methylases protects the recognition site 5'-CTCTTC-3'/5'-GAAGAG-3' from cleavage with the cognate restriction endonuclease BcoKI. It is shown that M.BcoKIB is an N6-adenine specific methylase and M.BcoKIA is an N4-cytosine specific methylase. According to bisulfite mapping, M.BcoKIA methylates the first cytosine in the sequence 5'-CTCTTC-3'.     

Specificity of DNA binding and methylation by the M.FokI DNA methyltransferase

The M.FokI adenine-N(6) DNA methyltransferase recognizes the asymmetric DNA sequence GGATG/CATCC. It consists of two domains each containing all motifs characteristic for adenine-N(6) DNA methyltransferases. We have studied the specificity of DNA-methylation by both domains using 27 hemimethylated oligonucleotide substrates containing recognition sites which differ in one or two base pairs from GGATG or CATCC. The N-terminal domain of M.FokI interacts very specifically with GGATG-sequences, because only one of the altered sites is modified. In contrast, the C-terminal domain shows lower specificity. It prefers CATCC-sequences but only two of the 12 star sites (i.e. sites that differ in 1 bp from the recognition site) are not accepted and some star sites are modified with rates reduced only 2-3-fold. In addition, GGATGC- and CGATGC-sites are modified which differ at two positions from CATCC. DNA binding experiments show that the N-terminal domain preferentially binds to hemimethylated GGATG/C(m)ATCC sequences whereas the C-terminal domain binds to DNA with higher affinity but without specificity. Protein-protein interaction assays show that both domains of M.FokI are in contact with each other. However, several DNA-binding experiments demonstrate that DNA-binding of both domains is mutually exclusive in full-length M.FokI and both domains do not functionally influence each other. The implications of these results on the molecular evolution of type IIS restriction/modification systems are discussed.     

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.     

AdoMet-dependent methyl-transfer: Glu119 is essential for DNA C5-cytosine methyltransferase M.HhaI

The role of Glu119 in S-adenosyl-L-methionine-dependent DNA methyltransferase M.HhaI-catalyzed DNA methylation was studied. Glu119 belongs to the highly conserved Glu/Asn/Val motif found in all DNA C5-cytosine methyltransferases, and its importance for M.HhaI function remains untested. We show that formation of the covalent intermediate between Cys81 and the target cytosine requires Glu119, since conversion to Ala, Asp or Gln lowers the rate of methyl transfer 10(2)-10(6) fold. Further, unlike the wild-type M.HhaI, these mutants are not trapped by the substrate in which the target cytosine is replaced with the mechanism-based inhibitor 5-fluorocytosine. The DNA binding affinity for the Glu119Asp mutant is decreased 10(3)-fold. Thus, the ability of the enzyme to stabilize the extrahelical cytosine is coupled directly to tight DNA binding. The structures of the ternary protein/DNA/AdoHcy complexes for both the Glu119Ala and Glu119Gln mutants (2.70 A and 2.75 A, respectively) show that the flipped base is positioned nearly identically with that observed in the wild-type M.HhaI complex. A single water molecule in the Glu119Ala structure between Ala119 and the extrahelical cytosine N3 is lacking in the Glu119Gln and wild-type M.HhaI structures, and most likely accounts for this mutant's partial activity. Glu119 has essential roles in activating the target cytosine for nucleophilic attack and contributes to tight DNA binding.     

Novel procedure for the detection of 5-methylcytosine.

Bisulfite converts non-methylated cytosine in DNA to uracil leaving 5-methylcytosine unaltered. In this communication, we present a new approach omitting the conventional PCR amplification step. Bisulfite-converted methylated DNA is directly sequenced. The effectiveness of the new protocol is demonstrated by using it for the detection of 5-methylation of cytosine residues introduced by three different DNA methyltransferases (M.HaeIII, M.HpaII and M.HhaI). A simple experimental system useful to determine the sequence specificity of DNA methyltransferases is also presented.