Agmenellum quadruplicatum M.AquI, a novel modification methylase

Circularly permuted variants of two CG-specific prokaryotic DNA methyltransferases

The highly similar prokaryotic DNA (cytosine-5) methyltransferases (C5-MTases) M.MpeI and M.SssI share the specificity of eukaryotic C5-MTases (5'-CG), and can be useful research tools in the study of eukaryotic DNA methylation and epigenetic regulation. In an effort to improve the stability and solubility of complementing fragments of the two MTases, genes encoding circularly permuted (CP) variants of M.MpeI and M.SssI were created, and cloned in a plasmid vector downstream of an arabinose-inducible promoter. MTase activity of the CP variants was tested by digestion of the plasmids with methylation-sensitive restriction enzymes. Eleven of the fourteen M.MpeI permutants and six of the seven M.SssI permutants had detectable MTase activity as indicated by the full or partial protection of the plasmid carrying the cpMTase gene. Permutants cp62M.MpeI and cp58M.SssI, in which the new N-termini are located between conserved motifs II and III, had by far the highest activity. The activity of cp62M.MpeI was comparable to the activity of wild-type M.MpeI. Based on the location of the split sites, the permutants possessing MTase activity can be classified in ten types. Although most permutation sites were designed to fall outside of conserved motifs, and the MTase activity of the permutants measured in cell extracts was in most cases substantially lower than that of the wild-type enzyme, the high proportion of circular permutation topologies compatible with MTase activity is remarkable, and is a new evidence for the structural plasticity of C5-MTases. A computer search of the REBASE database identified putative C5-MTases with CP arrangement. Interestingly, all natural circularly permuted C5-MTases appear to represent only one of the ten types of permutation topology created in this work.

Complementation between inactive fragments of SssI DNA methyltransferase

Background

Silencing mammalian genes by targeted DNA (cytosine-5) methylation of selected CG sites in the genome would be a powerful technique to analyze epigenomic information and to study the roles of DNA methylation in physiological and pathological states. A promising approach of targeted DNA methylation is based on the ability of split fragments of a monomeric DNA methyltransferase (C5-MTase) to associate and form active enzyme. A few C5-MTases of different specificities have been shown to possess the ability of fragment complementation, but a demonstration of this phenomenon for a C5-MTase, which has CG specificity and thus can be targeted to methylate any CG site, has been lacking. The purpose of this study was to test whether the CG-specific prokaryotic C5-MTase M.SssI shows the phenomenon of fragment complementation.

Results

We show that truncated inactive N-terminal fragments of M.SssI can assemble with truncated inactive C-terminal fragments to form active enzyme in vivo when produced in the same E. coli cell. Overlapping and non-overlapping fragments as well as fragments containing short appended foreign sequences had complementation capacity. In optimal combinations C-terminal fragments started between conserved motif VIII and the predicted target recognizing domain of M.SssI. DNA methyltransferase activity in crude extracts of cells with the best complementing fragment pairs was ~ 4 per cent of the activity of cells producing the full length enzyme. Fusions of two N-terminal and two C-terminal fragments to 21.6 kDa zinc finger domains only slightly reduced complementation ability of the fragments.

Conclusions

The CG-specific DNA methyltransferase M.SssI shows the phenomenon of fragment complementation in vivo in E. coli. Fusion of the split fragments to six unit zinc finger domains does not substantially interfere with the formation of active enzyme. These observations and the large number of complementing fragment combinations representing a wide range of MTase activity offer the possibility to develop M.SssI into a programmable DNA methyltransferase of high specificity.

Characterization of the large subunit of EcoHK31I methyltransferase by structural modeling and mutagenesis

M.EcoHK31I is a naturally occurring mC5-methyltransferase with a large alpha polypeptide and a small beta polypeptide. Polypeptide alpha contains conserved motifs I-VIII and X, and polypeptide beta contains motif IX. To understand how polypeptide alpha carries out its function, a molecular model of the large domain of polypeptide alpha was generated using M.HhaI and M.HaeIII as templates. The large domain is a mixed alpha/beta structure. Residues 15-19 in motif I (Phe-Naa-Gly-Naa) are conserved for cofactor binding. The key catalytic residue Cys-79 in motif IV is also conserved in comparison with other C-5 MTases. Comparing polypeptide alpha with M.HhaI and M.HaeIII revealed a unique region upstream of motif X. To understand the role of this region, 14 charged residues between R224 and E271 in the putative small domain were mutated. Activity assays indicated that most of these charges can be eliminated or changed conservatively. Among these charged residues, R224, E240, D245 and D251 may take part in proper interaction with DNA in the presence of polypeptide beta.

Functional studies of the small subunit of EcoHK31I DNA methyltransferase

EcoHK31I DNA methyltransferase recognizes the sequence 5'-YGGCCR-3' and adds a methyl group to the fifth position of the internal cytosine to protect the DNA from cleavage by its cognate endonuclease. M.EcoHK31I is composed of polypeptides alpha and beta. Polypeptide beta only contains the conserved IX motif of the C5-MTase family, and provides a unique example to show that this motif alone may be dislocated to another polypeptide. By electromobility shift assay, we found that the alpha/beta complex recognizes specific oligonucleotide substrates. Polypeptide alpha formed aggregates with DNA, while polypeptide beta alone did not bind DNA. Therefore, polypeptide beta assists in the proper binding of polypeptide alpha to DNA substrate. The complex of polypeptide alpha and a polypeptide beta variant with an N-terminal deletion of 41 amino acids showed a 16-fold reduction in methylation activity. Further deletion resulted in an inactive methyltransferase. The dissociation equilibrium constant (Kd) of the alpha/beta complex was 56.4 nM, while the Kd value for the alpha/deltaN46-polypeptide beta complex was increased approximately 95-fold, caused by a drastic decrease in dissociate rate constant (kd) and an increase in the association rate constant (ka). This indicates that the N-terminal region of polypeptide beta takes part in subunit interaction, while the C-terminal region is involved in DNA binding.

Protein fragment complementation in M.HhaI DNA methyltransferase

The 5mC DNA methyltransferase M.HhaI can be split into two individually inactive N- and C-terminal fragments that together can form an active enzyme in vivo capable of efficiently methylating DNA. This active fragment pair was identified by creating libraries of M.HhaI gene fragment pairs and then selecting for the pairs that code for an active 5mC methyltransferase. The site of bisection for successful protein fragment complementation in M.HhaI was in the variable region near the target recognition domain between motif VIII and TRD. This same region is the location of bifurcation in the naturally split 5mC methyltransferase M.AquI, the location for circular permutation in M.BssHII, and the location for previously engineered split versions of M.BspRI.

Structure of the Q237W mutant of HhaI DNA methyltransferase: an insight into protein-protein interactions

We have determined the structure of a mutant (Q237W) of HhaI DNA methyltransferase, complexed with the methyl-donor product AdoHcy. The Q237W mutant proteins were crystallized in the monoclinic space group C2 with two molecules in the crystallographic asymmetric unit. Protein-protein interface calculations in the crystal lattices suggest that the dimer interface has the specific characteristics for homodimer protein-protein interactions, while the two active sites are spatially independent on the outer surface of the dimer. The solution behavior suggests the formation of HhaI dimers as well. The same HhaI dimer interface is also observed in the previously characterized binary (M.HhaI-AdoMet) and ternary (M.HhaI-DNA-AdoHcy) complex structures, crystallized in different space groups. The dimer is characterized either by a non-crystallographic two-fold symmetry or a crystallographic symmetry. The dimer interface involves three segments: the amino-terminal residues 2-8, the carboxy-terminal residues 313-327, and the linker (amino acids 179-184) between the two functional domains--the catalytic methylation domain and the DNA target recognition domain. Both the amino- and carboxy-terminal segments are part of the methylation domain. We also examined protein-protein interactions of other structurally characterized DNA MTases, which are often found as a 2-fold related 'dimer' with the largest dimer interface area for the group-beta MTases. A possible evolutionary link between the Type I and Type II restriction-modification systems is discussed.

A nomenclature for restriction enzymes, DNA methyltransferases, homing endonucleases and their genes

A nomenclature is described for restriction endonucleases, DNA methyltransferases, homing endonucleases and related genes and gene products. It provides explicit categories for the many different Type II enzymes now identified and provides a system for naming the putative genes found by sequence analysis of microbial genomes.

The small subunit of M. AquI is responsible for sequence-specific DNA recognition and binding in the absence of the catalytic domain

AquI DNA methyltransferase (M. AquI) catalyzes the transfer of a methyl group from S-adenosyl-L-methionine to the C5 position of the outermost deoxycytidine base in the DNA sequence 5'-CCCGGG-3'. M. AquI is a heterodimer in which the polypeptide chain is separated at the junction between the two equivalent structural domains in the related enzyme M. HhaI. Recently, we reported the subcloning, overexpression, and purification of the subunits (alpha and beta) of M. AquI separately. Here we describe the DNA binding properties of M. AquI. The results presented here indicate that the beta subunit alone contains all of the information for sequence-specific DNA recognition and binding. The first step in the sequence-specific recognition of DNA by M. AquI involves the formation of binary complex with the target recognition domain in conjunction with conserved sequence motifs IX and X, found in all known C5 DNA methyltransferases, contained in the beta subunit. The alpha subunit enhances the binding of the beta subunit to DNA specifically and nonspecifically. It is likely that the addition of the alpha subunit to the beta subunit stabilizes the conformation of the beta subunit and thereby enhances its affinity for DNA indirectly. Addition of S-adenosyl-L-methionine and its analogues S-adenosyl-L-homocysteine and sinefungin enhances binding, but only in the presence of the alpha subunit. These compounds did not have any effect on DNA binding by the beta subunit alone. Using a 30-mer oligodeoxynucleotide substrate containing 5-fluorodeoxycytidine (5-FdC), it was found that the beta subunit alone did not form a covalent complex with its specific sequence in the absence or presence of S-adenosyl-L-methionine. However, the addition of the alpha subunit to the beta subunit led to the formation of a covalent complex with specific DNA sequence containing 5-FdC.

Sequence permutations in the molecular evolution of DNA methyltransferases

BACKGROUND

DNA methyltransferases (MTases), unlike MTases acting on other substrates, exhibit sequence permutation. Based on the sequential order of the cofactor-binding subdomain, the catalytic subdomain, and the target recognition domain (TRD), several classes of permutants have been proposed. The majority of known DNA MTases fall into the alpha, beta, and gamma classes. There is only one member of the zeta class known and no members of the delta and epsilon classes have been identified to date. Two mechanisms of permutation have been proposed: one involving gene duplication and in-frame fusion, and the other involving inter- and intragenic shuffling of gene segments.

RESULTS

Two novel cases of sequence permutation in DNA MTases implicated in restriction-modification systems have been identified, which suggest that members of the delta and zeta classes (M.MwoI and M.TvoORF1413P, respectively) evolved from beta-class MTases. This is the first identification of the delta-class MTase and the second known zeta-class MTase (the first zeta-class member among DNA:m4C and m6A-MTases).

CONCLUSIONS

Fragmentation of a DNA MTase gene may result from attack of nucleases, for instance when the RM system invades a new cell. Its reassembly into a functional form, the order of motifs notwithstanding, may be strongly selected for, if the cognate ENase gene remains active and poses a threat to the host's chromosome. The "cut-and-paste" mechanism is proposed for beta-delta permutation, which is non-circular and involves relocation of one segment of a gene. The circular beta-zeta permutation may be explained both by gene duplication or shuffling of gene fragments. These two mechanisms are not mutually exclusive and probably both played a role in the evolution of permuted DNA MTases.

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.

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.

Plant DNA methyltransferases

DNA methylation is an important modification of DNA that plays a role in genome management and in regulating gene expression during development. Methylation is carried out by DNA methyltransferases which catalyse the transfer of a methyl group to bases within the DNA helix. Plants have at least three classes of cytosine methyltransferase which differ in protein structure and function. The METI family, homologues of the mouse Dnmtl methyltransferase, most likely function as maintenance methyltransferases, but may also play a role in de novo methylation. The chromomethylases, which are unique to plants, may preferentially methylate DNA in heterochromatin; the remaining class, with similarity to Dnmt3 methyltransferases of mammals, are putative de novo methyltransferases. The various classes of methyltransferase may show differential activity on cytosines in different sequence contexts. Chromomethylases may preferentially methylate cytosines in CpNpG sequences while the Arabidopsis METI methyltransferase shows a preference for cytosines in CpG sequences. Additional proteins, for example DDM1, a member of the SNF2/SWI2 family of chromatin remodelling proteins, are also required for methylation of plant DNA.

Conserved plant genes with similarity to mammalian de novo DNA methyltransferases

DNA methylation plays a critical role in controlling states of gene activity in most eukaryotic organisms, and it is essential for proper growth and development. Patterns of methylation are established by de novo methyltransferases and maintained by maintenance methyltransferase activities. The Dnmt3 family of de novo DNA methyltransferases has recently been characterized in animals. Here we describe DNA methyltransferase genes from both Arabidopsis and maize that show a high level of sequence similarity to Dnmt3, suggesting that they encode plant de novo methyltransferases. Relative to all known eukaryotic methyltransferases, these plant proteins contain a novel arrangement of the motifs required for DNA methyltransferase catalytic activity. The N termini of these methyltransferases contain a series of ubiquitin-associated (UBA) domains. UBA domains are found in several ubiquitin pathway proteins and in DNA repair enzymes such as Rad23, and they may be involved in ubiquitin binding. The presence of UBA domains provides a possible link between DNA methylation and ubiquitin/proteasome pathways.

M.(phi)BssHII, a novel cytosine-C5-DNA-methyltransferase with target-recognizing domains at separated locations of the enzyme

In all cytosine-C5-DNA-methyltransferases (MTases) from prokaryotes and eukaryotes, remarkably conserved amino acid sequence elements responsible for general enzymatic functions are arranged in the same canonical order. In addition, one variable region, which includes the target-recognizing domain(s) (TRDs) characteristic for each enzyme, has been localized in one region between the same blocks of these conserved elements. This conservation in the order of conserved and variable sequences suggests stringent structural constraints in the primary structure to obtain the correct folding of the enzymes. Here we report the characterization of a new type of a multispecific MTase, M.(phiphi)BssHII, which is expressed as two isoforms. Isoform I is an entirely novel type of MTase which has, in addition to the TRDs at the conventional location, one TRD located at a non-canonical position at its N-terminus. Isoform II is represented by the same MTase, but without the N-terminal TRD. The N-terminal TRD provides HaeII methylation specificity to isoform I. The TRD is fully functional when engineered into either the conventional variable region of M.(phiphi)BssHII or the related monospecific M.phi3TII MTase. The implications of this structural plasticity with respect to the evolution of MTases are discussed.

The cyanobacterium Synechocystis sp. strain PCC 6803 expresses a DNA methyltransferase specific for the recognition sequence of the restriction endonuclease PvuI

By use of restriction endonucleases, the DNA of the cyanobacterium Synechocystis sp. strain PCC 6803 was analyzed for DNA-specific methylation. Three different recognition sites of methyltransferases, a dam-like site including N6-methyladenosine and two other sites with methylcytosine, were identified, whereas no activities of restriction endonucleases could be detected in this strain. slr0214, a Synechocystis gene encoding a putative methyltransferase that shows significant similarities to C5-methylcytosine-synthesizing enzymes, was amplified by PCR and cloned for further characterization. Mutations in slr0214 were generated by the insertion of an aphII gene cassette. Analyses of chromosomal DNAs of such mutants demonstrated that the methylation pattern was changed. The recognition sequence of the methyltransferase was identified as 5'-CGATCG-3', corresponding to the recognition sequence of PvuI. The specific methyltransferase activity was significantly reduced in protein extracts obtained from mutant cells. Mutation of slr0214 also led to changed growth characteristics of the cells compared to wild-type cells. These alterations led to the conclusion that the methyltransferase Slr0214 might play a regulatory role in Synechocystis. The Slr0214 protein was also overexpressed in Escherichia coli, and the purified protein demonstrated methyltransferase activity and specificity for PvuI recognition sequences in vitro. We propose the designation M.Ssp6803I [corrected] (Synechocystis methyltransferase I) for the slr0214-encoded enzyme.

A DNA methyltransferase homolog with a chromodomain exists in multiple polymorphic forms in Arabidopsis

Chromodomains are thought to mediate protein-protein interactions between chromatin components. We have detected a chromodomain embedded within the catalytic region of a predicted Arabidopsis DNA methyltransferase that is diverged from other eukaryotic enzymes. The 791 residue "chromomethylase" (CMT1) is encoded by a floral transcript that is spliced from 20 exons and is present at only approximately 1/10(-7) of total mRNA. Genomic sequencing reveals an ancient haplotype split at CMT1 between Col-0 + Metz and the other ecotypes examined. In the Col-0 + Metz haplotype, alternative mRNA processing at intron 13 truncates the coding region. In Ler, RLD, and No-0, similar truncation is caused by insertion of an intact retrotransposon, Evelknievel, which is present as a single copy in Ler and RLD and is currently methylated and inactive. Evelknievel is found at this site on a single branch that connects the Ler, RLD, and No-0 ecotypes but is absent from the genomes of all other ecotypes examined. A stop codon within exon 6 of the Metz ecotype confirms that CMT1 is nonessential. Nevertheless, comparison to CMT1 of Cardaminopsis arenosa, an outcrossing relative, indicates conservation for DNA methyltransferase function. We discuss how allelic diversity of CMT1 may reflect loosened selective constraints in a self-fertilizing species such as Arabidopsis thaliana.

Identification of a DNA cytosine methyltransferase gene in conjugative transposon Tn5252

The nucleotide sequence of the 3.5-kb right junction fragment of the streptococcal conjugative transposon Tn5252 was obtained. The DNA fragment was found to carry four putative genes one of which displayed a high degree of similarity to prokaryotic 5C-cytosine methyltransferases carrying multiple sequence specificities. No cognate endonuclease gene was detected in the sequenced DNA. Purified methylase polypeptide synthesized in a T7 promoter-controlled overexpression system was found to lack methylase activity while the cell extracts of host cells containing the recombinant plasmid carrying the methylase gene were active. In vivo mutations in the methylase gene did not seem to affect the transferability of the element.

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

BssHII restriction endonuclease cleaves 5'-GCGCGC-3' on double-stranded DNA between the first and second bases to generate a four base 5'overhang. BssHII restriction endonuclease was purified from the native Bacillus stearothermophilus H3 cells and its N-terminal amino acid sequence was determined. Degenerate PCR primers were used to amplify the first 20 codons of the BssHII restriction endonuclease gene. The BssHII restriction endonuclease gene (bssHIIR) and the cognate BssHII methyltransferase gene (bssHIIM) were cloned in Escherichia coli by amplification of Bacillus stearothermophilus genomic DNA using PCR and inverse PCR. BssHII methyltransferase (M.BssHII) contains all 10 conserved cytosine-5 methyltransferase motifs, but motifs IX and X precede motifs I-VIII. Thus, the conserved motifs of M. BssHII are circularly permuted relative to the motif organizations of other cytosine-5 methyltransferases. M.BssHII and the non-cognate multi-specific phiBssHII methyltransferase, M.phiBss HII [Schumann,J. et al . (1995) Gene, 157, 103-104] share 34% identity in amino acid sequences from motifs I-VIII, and 40% identity in motifs IX-X. A conserved arginine is located upstream of a TV dipeptide in the N-terminus of M.BssHII that may be responsible for the recognition of the guanine 5' of the target cytosine. The BssHII restriction endonuclease gene was expressed in E.coli via a T7 expression vector.

A prediction of the amino acids and structures involved in DNA recognition by type I DNA restriction and modification enzymes

The S subunits of type I DNA restriction/modification enzymes are responsible for recognising the DNA target sequence for the enzyme. They contain two domains of approximately 150 amino acids, each of which is responsible for recognising one half of the bipartite asymmetric target. In the absence of any known tertiary structure for type I enzymes or recognisable DNA recognition motifs in the highly variable amino acid sequences of the S subunits, it has previously not been possible to predict which amino acids are responsible for sequence recognition. Using a combination of sequence alignment and secondary structure prediction methods to analyse the sequences of S subunits, we predict that all of the 51 known target recognition domains (TRDs) have the same tertiary structure. Furthermore, this structure is similar to the structure of the TRD of the C5-cytosine methyltransferase, Hha I, which recognises its DNA target via interactions with two short polypeptide loops and a beta strand. Our results predict the location of these sequence recognition structures within the TRDs of all type I S subunits.

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.

Chemistry and biology of DNA methyltransferases

Recognition of a specific DNA sequence by a protein is probably the best example of macromolecular interactions leading to various events. It is a prerequisite to understanding the basis of protein-DNA interactions to obtain a better insight into fundamental processes such as transcription, replication, repair, and recombination. DNA methyltransferases with varying sequence specificities provide an excellent model system for understanding the molecular mechanism of specific DNA recognition. Sequence comparison of cloned genes, along with mutational analyses and recent crystallographic studies, have clearly defined the functions of various conserved motifs. These enzymes access their target base in an elegant manner by flipping it out of the DNA double helix. The drastic protein-induced DNA distortion, first reported for HhaI DNA methyltransferase, appears to be a common mechanism employed by various proteins that need to act on bases. A remarkable feature of the catalytic mechanism of DNA (cytosine-5) methyltransferases is the ability of these enzymes to induce deamination of the target cytosine in the absence of S-adenosyl-L-methionine or its analogs. The enzyme-catalyzed deamination reaction is postulated to be the major cause of mutational hotspots at CpG islands responsible for various human genetic disorders. Methylation of adenine residues in Escherichia coli is known to regulate various processes such as transcription, replication, repair, recombination, transposition, and phage packaging.

Cloning and analysis of the genes encoding the type IIS restriction-modification system HphI from Haemophilus parahaemolyticus

The genomic region encoding the type IIS restriction-modification (R-M) system HphI (enzymes recognizing the asymmetric sequence 5'-GGTGA-3'/5'-TCACC-3') from Haemophilus parahaemolyticus were cloned into Escherichia coli and sequenced. Sequence analysis of the R-M HphI system revealed three adjacent genes aligned in the same orientation: a cytosine 5 methyltransferase (gene hphIMC), an adenine N6 methyltransferase (hphIMA) and the HphI restriction endonuclease (gene hphIR). Either methyltransferase is capable of protecting plasmid DNA in vivo against the action of the cognate restriction endonuclease. hphIMA methylation renders plasmid DNA resistant to R.Hindill at overlapping sites, suggesting that the adenine methyltransferase modifies the 3'-terminal A residue on the GGTGA strand. Strong homology was found between the N-terminal part of the m6A methyltransferasease and an unidentified reading frame interrupted by an incomplete gaIE gene of Neisseria meningitidis. The HphI R-M genes are flanked by a copy of a 56 bp direct nucleotide repeat on each side. Similar sequences have also been identified in the non-coding regions of H.influenzae Rd DNA. Possible involvement of the repeat sequences in the mobility of the HphI R-M system is discussed.

Overproduction, purification and characterization of M.EcoHK31I, a bacterial methyltransferase with two polypeptides

The two overlapping genes coding for EcoHK31I methyltransferase have previously been cloned, sequenced and expressed [Lee, Kam and Shaw (1995) Nucleic Acids Res. 23, 103-108]. Here we describe protocols developed to purify polypeptides alpha and beta together or separately, to apparent homogeneity by various chromatographic media. M.EcoHK31I is a heterodimer with a native molecular mass of 61 kDa. Its specific activity towards non-methylated lambda DNA was 3.0 x 10(5) units per mg of protein. The respective denatured molecular masses of polypeptides alpha and beta were 38 and 23 kDa, and their pI values were 8.7 and 6.8. Initial rate kinetic parameters of the native enzyme were 2.0 nM, 0.58 microM and 3 min-1 for KmDNA, KmAdoMet and kcat. respectively, where AdoMet stands for S-adenosyl-L-methionine. Fully active enzyme was reconstituted by co-purifying the two separately synthesized polypeptides, and activity assays confirmed our previous finding that two polypeptides were needed to methylate substrate DNA.

A bacterial methyltransferase M.EcoHK311 requires two proteins for in vitro methylation

The genes encoding EcoHK311 restriction-modification (R-M) system were isolated from a clinically-isolated Escherichia coli strain HK31. The entire R-M system of EcoHK311 is located in a 2.1 kb fragment. R.EcoHK311 is an isoschizomer of Eael which recognizes and cleaves Y decreases GGCCR. M.EcoHK31l consists of two polypeptides alpha and beta with sizes 309 and 176 aa, respectively. Polypeptide beta is encoded within aa, alternative reading frame of polypeptide alpha. All the conserved motifs in mC5-MTases can be found in polypeptide alpha except motif IX which is present in polypeptide beta. Polypeptides alpha and beta were separately synthesized in a T7 promoter controlled over-expression system and in vitro methylation occurred only when the two extracts were mixed and thus confirms that two polypeptides are required for methylation.

M.phi 3TII: a new monospecific DNA (cytosine-C5) methyltransferase with pronounced amino acid sequence similarity to a family of adenine-N6-DNA-methyltransferases

The temperate B.subtilis phages phi 3T and rho 11s code, in addition to the multispecific DNA (cytosine-C5) methyltransferases (C5-MTases) M. phi 3TI and M. rho 11sI, which were previously characterized, for the identical monospecific C5-MTases M. phi 3TII and M. rho 11sII. These enzymes modify the C of TCGA sites, a novel target specificity among C5-MTases. The primary sequence of M. phi 3TII (326 amino acids) shows all conserved motifs typical of the building plan of C5-MTases. The degree of relatedness between M. phi 3TII and all other mono- or multispecific C5-MTases ranges from 30-40% amino acid identity. Particularly M. phi 3TII does not show pronounced similarity to M. phi 3TI indicating that both MTase genes were not generated from one another but were acquired independently by the phage. The amino terminal part of the M. phi 3TII (preceding the variable region 'V'), which predominantly constitutes the catalytic domain of the enzyme, exhibits pronounced sequence similarity to the amino termini of a family of A-N6-MTases, which--like M.TaqI--recognize the general sequence TNNA. This suggests that recently described similarities in the general three dimensional organization of C5- and A-N6-MTases imply divergent evolution of these enzymes originating from a common molecular ancestor.

The DNA recognition subunit of a DNA methyltransferase is predominantly a molten globule in the absence of DNA

Enzyme-catalysed DNA methylation provides an opportunity for the modulation of protein-DNA recognition in biological systems. Recently we have demonstrated that the smaller of the two subunits of the heterodimeric, cytosine-specific DNA methyltransferase, M. AquI, is largely responsible for sequence-specific DNA recognition. Here we present evidence from a series of NMR, fluorescence and circular dichroism spectroscopy experiments that the DNA binding subunit of M. AquI has the characteristics of a molten globule in the absence of the catalytic machinery. In this metastable state this subunit retains its ability to bind DNA in a sequence-specific manner. We believe this finding offers an insight into the structural flexibility which underpins the mechanism of action of these enzymes, and may provide a possible biological role for molten globules in protein function.

M.phi 3TII: a new monospecific DNA (cytosine-C5) methyltransferase with pronounced amino acid sequence similarity to a family of adenine-N6-DNA-methyltransferases

The temperate B.subtilis phages phi 3T and rho 11s code, in addition to the multispecific DNA (cytosine-C5) methyltransferases (C5-MTases) M.phi 3TI and M.rho 11sI, which were previously characterized, for the identical monospecific C5-MTases M.phi 3TII and M.rho 11sII. These enzymes modify the C to TCGA sites, a novel target specificity among C5-MTases. The primary sequence of M.phi 3TII (326 amino acids) shows all conserved motifs typical of the building plan of C5-MTases. The degree of relatedness between M.phi 3TII and all other mono- or multispecific C5-MTases ranges from 30-40% amino acid identity. Particularly M.phi 3TII does not show pronounced similarity to M.phi 3TI indicating that both MTase genes were not generated from one another but were acquired independently by the phage. The amino terminal part of the M.phi 3TII (preceding the variable region 'V'), which predominantly constitutes the catalytic domain of the enzyme, exhibits pronounced sequence similarity to the amino termini of a family of A-N6-MTases, which--like M.Taql--recognize the general sequence TNNA. This suggests that recently described similarities in the general three dimensional organization of C5- and A-N6-MTases imply divergent evolution of these enzymes originating from a common molecular ancestor.

The DNA (cytosine-5) methyltransferases

The m5C-MTases form a closely-knit family of enzymes in which common amino acid sequence motifs almost certainly translate into common structural and functional elements. These common elements are located predominantly in a single structural domain that performs the chemistry of the reaction. Sequence-specific DNA recognition is accomplished by a separate domain that contains recognition elements not seen in other structures. This, combined with the novel and unexpected mechanistic feature of trapping a base out of the DNA helix, makes the m5C-MTases an intriguing class of enzymes for further study. The reaction pathway has suddenly become more complicated because of the base-flipping and much remains to be learned about the DNA recognition elements in the family members for which structural information is not yet available.

Physical genome map of the unicellular cyanobacterium Synechococcus sp. strain PCC 7002

A physical restriction map of the genome of the cyanobacterium Synechococcus sp. strain PCC 7002 was assembled from AscI, NotI, SalI, and SfiI digests of intact genomic DNA separated on a contour-clamped homogeneous electric field pulsed-field gel electrophoresis system. An average genome size of 2.7 x 10(6) bp was calculated from 21 NotI, 37 SalI, or 27 SfiI fragments obtained by the digestions. The genomic map was assembled by using three different strategies: linking clone analysis, pulsed-field fragment hybridization, and individual clone hybridization to singly and doubly restriction-digested large DNA fragments. The relative positions of 21 genes or operons were determined, and these data suggest that the gene order is not highly conserved between Synechococcus sp. strain PCC 7002 and Anabaena sp. strain PCC 7120.

The M.AluI DNA-(cytosine C5)-methyltransferase has an unusually large, partially dispensable, variable region

The DNA methyltransferase of the AluI restriction-modification system, from Arthrobacter luteus, converts cytosine to 5-methylcytosine in the sequence AGCT. The gene for this methyltransferase, aluIM, was cloned into Escherichia coli and sequenced. A 525-codon open reading frame was found, consistent with deletion evidence, and the deduced amino acid sequence revealed all ten conserved regions common to 5-methylcytosine methyltransferases. The aluIM sequence predicts a protein of M(r) 59.0k, in agreement with the observed M(r), making M.AluI the largest known methyltransferase from a type II restriction-modification system. M.AluI also contains the largest known variable region of any monospecific DNA methyltransferase, larger than that of most multispecific methyltransferases. In other DNA methyltransferases the variable region has been implicated as the sequence-specific target recognition domain. An in-frame deletion that removes a third of this putative target-recognition region leaves the Alu I methyltransferase still fully active.

Cloning and sequence analysis of the genes coding for Eco57I type IV restriction-modification enzymes

A 6.3 kb fragment of E.coli RFL57 DNA coding for the type IV restriction-modification system Eco57I was cloned and expressed in E.coli RR1. A 5775 bp region of the cloned fragment was sequenced which contains three open reading frames (ORF). The methylase gene is 1623 bp long, corresponding to a protein of 543 amino acids (62 kDa); the endonuclease gene is 2991 bp in length (997 amino acids, 117 kDa). The two genes are transcribed convergently from different strands with their 3'-ends separated by 69 bp. The third short open reading frame (186 bp, 62 amino acids) has been identified, that precedes and overlaps by 7 nucleotides the ORF encoding the methylase. Comparison of the deduced Eco57I endonuclease and methylase amino acid sequences revealed three regions of significant similarity. Two of them resemble the conserved sequence motifs characteristic of the DNA[adenine-N6] methylases. The third one shares similarity with corresponding regions of the PaeR7I, TaqI, CviBIII, PstI, BamHI and HincII methylases. Homologs of this sequence are also found within the sequences of the PaeR7I, PstI and BamHI restriction endonucleases. This is the first example of a family of cognate restriction endonucleases and methylases sharing homologous regions. Analysis of the structural relationship suggests that the type IV enzymes represent an intermediate in the evolutionary pathway between the type III and type II enzymes.

Automated assembly of protein blocks for database searching

A system is described for finding and assembling the most highly conserved regions of related proteins for database searching. First, an automated version of Smith's algorithm for finding motifs is used for sensitive detection of multiple local alignments. Next, the local alignments are converted to blocks and the best set of non-overlapping blocks is determined. When the automated system was applied successively to all 437 groups of related proteins in the PROSITE catalog, 1764 blocks resulted; these could be used for very sensitive searches of sequence databases. Each block was calibrated by searching the SWISS-PROT database to obtain a measure of the chance distribution of matches, and the calibrated blocks were concatenated into a database that could itself be searched. Examples are provided in which distant relationships are detected either using a set of blocks to search a sequence database or using sequences to search the database of blocks. The practical use of the blocks database is demonstrated by detecting previously unknown relationships between oxidoreductases and by evaluating a proposed relationship between HIV Vif protein and thiol proteases.

Complementation by detached parts of GGCC-specific DNA methyltransferases

Individually inactive N- and C-terminal fragments of the m5C-methyltransferase M.BspRI can complement each other resulting in specific, in vivo methylation of the DNA. This was shown by cloning the coding regions for N- and C-terminal parts of the enzyme in compatible plasmids and co-transforming them into E.coli cells. The enzyme could be detached at several different sites, producing either non-overlapping or partially overlapping fragments capable of complementation. Reconstitution of the active methyltransferase from inactive fragments was demonstrated in vitro, as well. Another GGCC-specific methyltransferase, M.BsuRI, showed a similar complementation phenomenon. Moreover, interspecies complementation was observed between appropriate fragments of the two closely related enzymes M.BspRI and M.BsuRI. Fragments of structurally and functionally more different methyltransferases were unable to complement each other.

Organization of restriction-modification systems

The genes for over 100 restriction-modification systems have now been cloned, and approximately one-half have been sequenced. Despite their similar function, they are exceedingly heterogeneous. The heterogeneity is evident at three levels: in the gene arrangements; in the enzyme compositions; and in the protein sequences. This paper summarizes the main features of the R-M systems that have been cloned.

Stepwise cloning and molecular characterization of the HgiDI restriction-modification system from Herpetosiphon giganteus Hpa2

The restriction-modification system HgiDI from Herpetosiphon giganteus strain Hpa2 has been cloned in E. coli in a two-step procedure. Selection of the methyltransferase (M.HgiDI) gene in vitro was performed using the heterologous restriction endonuclease AhaII, an isoschizomer of Acyl and HgiDI (GRCGYC). Cloning of the complete HgiDI endonuclease (R.HgiDI) gene could only be achieved in recipient cells harbouring a recombinant plasmid, which was expressing the corresponding methyltransferase and could thereby prevent the host from self-destruction of its genetic material. The HgiDI restriction-modification system was sequenced and functionally correlated with two open reading frames of 309 (M) and 359 (R) codons. In homology studies M.HgiDI showed significant similarities to 20 other m5C-methyltransferases and turned out to be the most compact enzyme of this group described so far. Initial attempts for overexpression of M.HgiDI and partial purification of R.HgiDI have been successful.

Cloning and nucleotide sequence of the genes coding for the Sau96I restriction and modification enzymes

The genes coding for the GGNCC specific Sau96I restriction and modification enzymes were cloned and expressed in E. coli. The DNA sequence predicts a 430 amino acid protein (Mr: 49,252) for the methyltransferase and a 261 amino acid protein (Mr: 30,486) for the endonuclease. No protein sequence similarity was detected between the Sau96I methyltransferase and endonuclease. The methyltransferase contains the sequence elements characteristic for m5C-methyltransferases. In addition to this, M.Sau96I shows similarity, also in the variable region, with one m5C-methyltransferase (M.SinI) which has closely related recognition specificity (GGA/TCC). M.Sau96I methylates the internal cytosine within the GGNCC recognition sequence. The Sau96I endonuclease appears to act as a monomer.

Cloning and characterization of the HpaII methylase gene

The HpaII restriction-modification system from Haemophilus parainfluenzae recognizes the DNA sequence CCGG. The gene for the HpaII methylase has been cloned into E. coli and its nucleotide sequence has been determined. The DNA of the clones is fully protected against cleavage by the HpaII restriction enzyme in vitro, indicating that the methylase gene is active in E. coli. The clones were isolated in an McrA-strain of E. coli; attempts to isolate them in an McrA+ strain were unsuccessful. The clones do not express detectable HpaII restriction endonuclease activity, suggesting that either the endonuclease gene is not expressed well in E. coli, or that it is not present in its entirety in any of the clones that we have isolated. The derived amino acid sequence of the HpaII methylase shows overall similarity to other cytosine methylases. It bears a particularly close resemblance to the sequences of the HhaI, BsuFI and MspI methylases. When compared with three other methylases that recognize CCGG, the variable region of the HpaII methylase, which is believed to be responsible for sequence specific recognition, shows some similarity to the corresponding regions of the BsuFI and MspI methylases, but is rather dissimilar to that of the SPR methylase.