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

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

Shaping eukaryotic epigenetic systems by horizontal gene transfer

DNA methylation constitutes one of the pillars of epigenetics, relying on covalent bonds for addition and/or removal of chemically distinct marks within the major groove of the double helix. DNA methyltransferases, enzymes which introduce methyl marks, initially evolved in prokaryotes as components of restriction-modification systems protecting host genomes from bacteriophages and other invading foreign DNA. In early eukaryotic evolution, DNA methyltransferases were horizontally transferred from bacteria into eukaryotes several times and independently co-opted into epigenetic regulatory systems, primarily via establishing connections with the chromatin environment. While C5-methylcytosine is the cornerstone of plant and animal epigenetics and has been investigated in much detail, the epigenetic role of other methylated bases is less clear. The recent addition of N4-methylcytosine of bacterial origin as a metazoan DNA modification highlights the prerequisites for foreign gene co-option into the host regulatory networks, and challenges the existing paradigms concerning the origin and evolution of eukaryotic regulatory systems.

Bacterial DNA methyltransferase: A key to the epigenetic world with lessons learned from proteobacteria

Epigenetics modulates expression levels of various important genes in both prokaryotes and eukaryotes. These epigenetic traits are heritable without any change in genetic DNA sequences. DNA methylation is a universal mechanism of epigenetic regulation in all kingdoms of life. In bacteria, DNA methylation is the main form of epigenetic regulation and plays important roles in affecting clinically relevant phenotypes, such as virulence, host colonization, sporulation, biofilm formation et al. In this review, we survey bacterial epigenomic studies and focus on the recent developments in the structure, function, and mechanism of several highly conserved bacterial DNA methylases. These methyltransferases are relatively common in bacteria and participate in the regulation of gene expression and chromosomal DNA replication and repair control. Recent advances in sequencing techniques capable of detecting methylation signals have enabled the characterization of genome-wide epigenetic regulation. With their involvement in critical cellular processes, these highly conserved DNA methyltransferases may emerge as promising targets for developing novel epigenetic inhibitors for biomedical applications.

Beta class amino methyltransferases from bacteria to humans: evolution and structural consequences

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

The conserved aspartate in motif III of b family AdoMet-dependent DNA methyltransferase is important for methylation

S-adenosyl-L-methionine (AdoMet)-dependent methyltransferases (MTases) are involved in diverse cellular functions. These enzymes show little sequence conservation but have a conserved structural fold. The DNA MTases have characteristic motifs that are involved in AdoMet binding, DNA target recognition and catalysis. Motif III of these MTases have a highly conserved acidic residue, often an aspartate, whose functional significance is not clear. Here, we report a mutational study of the residue in the β family MTase of the Type III restriction-modification enzyme EcoP15I. Replacement of this residue by alanine affects its methylation activity. We propose that this residue contributes to the affinity of the enzyme for AdoMet. Analysis of the structures of DNA, RNA and protein MTases reveal that the acidic residue is conserved in all of them, and interacts with N6 of the adenine moiety of AdoMet. Interestingly, in the SET-domain protein lysine MTases, which have a fold different from other AdoMet-dependent MTases, N6 of the adenine moiety is hydrogen bonded to the main chain carbonyl group of the histidine residue of the highly conserved motif III. Our study reveals the evolutionary conservation of a carbonyl group in DNA, RNA and protein AdoMet-dependent MTases for specific interaction by hydrogen bond with AdoMet, despite the lack of overall sequence conservation.

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

The Caulobacter crescentus cell cycle-regulated DNA methyltransferase (CcrM) methylates the adenine of hemimethylated GANTC after replication. Here we present the structure of CcrM in complex with double-stranded DNA containing the recognition sequence. CcrM contains an N-terminal methyltransferase domain and a C-terminal nonspecific DNA-binding domain. CcrM is a dimer, with each monomer contacting primarily one DNA strand: the methyltransferase domain of one molecule binds the target strand, recognizes the target sequence, and catalyzes methyl transfer, while the C-terminal domain of the second molecule binds the non-target strand. The DNA contacts at the 5-base pair recognition site results in dramatic DNA distortions including bending, unwinding and base flipping. The two DNA strands are pulled apart, creating a bubble comprising four recognized base pairs. The five bases of the target strand are recognized meticulously by stacking contacts, van der Waals interactions and specific Watson–Crick polar hydrogen bonds to ensure high enzymatic specificity.

CcrM is a cell cycle-regulated DNA methyltransferase that methylates an adenine within a specific sequence following replication in the gram negative bacterium Caulobacter crescentus. Here the authors present a crystal structure of DNA-bound CcrM that reveals the molecular mechanism leading to sequence-specific methylation.

A model for the evolution of prokaryotic DNA restriction-modification systems based upon the structural malleability of Type I restriction-modification enzymes

Restriction Modification (RM) systems prevent the invasion of foreign genetic material into bacterial cells by restriction and protect the host's genetic material by methylation. They are therefore important in maintaining the integrity of the host genome. RM systems are currently classified into four types (I to IV) on the basis of differences in composition, target recognition, cofactors and the manner in which they cleave DNA. Comparing the structures of the different types, similarities can be observed suggesting an evolutionary link between these different types. This work describes the ‘deconstruction’ of a large Type I RM enzyme into forms structurally similar to smaller Type II RM enzymes in an effort to elucidate the pathway taken by Nature to form these different RM enzymes. Based upon the ability to engineer new enzymes from the Type I ‘scaffold’, an evolutionary pathway and the evolutionary pressures required to move along the pathway from Type I RM systems to Type II RM systems are proposed. Experiments to test the evolutionary model are discussed.

iDNA4mC: identifying DNA N4-methylcytosine sites based on nucleotide chemical properties

Motivation

DNA N4-methylcytosine (4mC) is an epigenetic modification. The knowledge about the distribution of 4mC is helpful for understanding its biological functions. Although experimental methods have been proposed to detect 4mC sites, they are expensive for performing genome-wide detections. Thus, it is necessary to develop computational methods for predicting 4mC sites.

Results

In this work, we developed iDNA4mC, the first webserver to identify 4mC sites, in which DNA sequences are encoded with both nucleotide chemical properties and nucleotide frequency. The predictive results of the rigorous jackknife test and cross species test demonstrated that the performance of iDNA4mC is quite promising and holds high potential to become a useful tool for identifying 4mC sites.

Availability and implementation

The user-friendly web-server, iDNA4mC, is freely accessible at http://lin.uestc.edu.cn/server/iDNA4mC.

Contact

chenweiimu@gmail.com or hlin@uestc.edu.cn.

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.

Structural insights into the RNA methyltransferase domain of METTL16

N⁶-methyladenosine (m⁶A) is an abundant modification in messenger RNA and noncoding RNAs that affects RNA metabolism. Methyltransferase-like protein 16 (METTL16) is a recently confirmed m⁶A RNA methyltransferase that methylates U6 spliceosomal RNA and interacts with the 3′-terminal RNA triple helix of MALAT1 (metastasis-associated lung adenocarcinoma transcript 1). Here, we present two X-ray crystal structures of the N-terminal methyltransferase domain (residues 1–291) of human METTL16 (METTL16_291): an apo structure at 1.9 Å resolution and a post-catalytic S-adenosylhomocysteine-bound complex at 2.1 Å resolution. The structures revealed a highly conserved Rossmann fold that is characteristic of Class I S-adenosylmethionine-dependent methyltransferases and a large, positively charged groove. This groove likely represents the RNA-binding site and it includes structural elements unique to METTL16. In-depth analysis of the active site led to a model of the methyl transfer reaction catalyzed by METTL16. In contrast to the major m⁶A methyltransferase heterodimer METTL3/METTL14, full-length METTL16 forms a homodimer and METTL16_291 exists as a monomer based on size-exclusion chromatography. A native gel-shift assay shows that METTL16 binds to the MALAT1 RNA triple helix, but monomeric METTL16_291 does not. Our results provide insights into the molecular structure of METTL16, which is distinct from METTL3/METTL14.

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

DNA N⁶-methyladenine modification plays an important role in regulating a variety of biological functions in bacteria. However, the mechanism of sequence-specific recognition in N⁶-methyladenine modification remains elusive. M1.HpyAVI, a DNA N⁶-adenine methyltransferase from Helicobacter pylori, shows more promiscuous substrate specificity than other enzymes. Here, we present the crystal structures of cofactor-free and AdoMet-bound structures of this enzyme, which were determined at resolutions of 3.0 Å and 3.1 Å, respectively. The core structure of M1.HpyAVI resembles the canonical AdoMet-dependent MTase fold, while the putative DNA binding regions considerably differ from those of the other MTases, which may account for the substrate promiscuity of this enzyme. Site-directed mutagenesis experiments identified residues D29 and E216 as crucial amino acids for cofactor binding and the methyl transfer activity of the enzyme, while P41, located in a highly flexible loop, playing a determinant role for substrate specificity. Taken together, our data revealed the structural basis underlying DNA N⁶-adenine methyltransferase substrate promiscuity.

Structural Basis of the Selectivity of GenN, an Aminoglycoside N-Methyltransferase Involved in Gentamicin Biosynthesis

Gentamicins are heavily methylated, clinically valuable pseudotrisaccharide antibiotics produced by Micromonospora echinospora. GenN has been characterized as an S-adenosyl-l-methionine-dependent methyltransferase with low sequence similarity to other enzymes. It is responsible for the 3″-N-methylation of 3″-dehydro-3″-amino-gentamicin A2, an essential modification of ring III in the biosynthetic pathway to the gentamicin C complex. Purified recombinant GenN also efficiently catalyzes 3″-N-methylation of related aminoglycosides kanamycin B and tobramycin, which both contain an additional hydroxymethyl group at the C5″ position in ring III. We have obtained eight cocrystal structures of GenN, at a resolution of 2.2 Å or better, including the binary complex of GenN and S-adenosyl-l-homocysteine (SAH) and the ternary complexes of GenN, SAH, and several aminoglycosides. The GenN structure reveals several features not observed in any other N-methyltransferase that fit it for its role in gentamicin biosynthesis. These include a novel N-terminal domain that might be involved in protein:protein interaction with upstream enzymes of the gentamicin X2 biosynthesis and two long loops that are involved in aminoglycoside substrate recognition. In addition, the analysis of structures of GenN in complex with different ligands, supported by the results of active site mutagenesis, has allowed us to propose a catalytic mechanism and has revealed the structural basis for the surprising ability of native GenN to act on these alternative substrates.

N6-Adenosine DNA Methyltransferase from H. pylori 98-10 Strain in Complex with DNA and AdoMet: Structural Insights from in Silico Studies

Helicobacter pylori is a primitive Gram-negative bacterium that resides in the acidic environment of the human gastrointestinal tract, and some strains of this bacterium cause gastric ulcers and cancer. DNA methyltransferases (MTases) are promising drug targets for the treatment of cancer and other diseases that are also caused by epigenetic alternations of the genome. The N6-adenine-specific DNA MTase from H. pylori (M. Hpy N6mA) catalyzes the transfer of a methyl group from the cofactor S-adenosyl-l-methionine (AdoMet) to the flipped adenine of the substrate DNA. In this work, we report the sequence analyses, three-dimensional structure modeling, and molecular dynamics simulations of M. Hpy N6mA, when complexed with AdoMet as well as DNA. We analyzed the protein-DNA interactions prominently established by the flipped cytosine and the interactions between protein cofactors in the active site. The comparable orientation of AdoMet in both systems confirms that AdoMet is in a catalytically competent orientation in the bimolecular system that is retained upon DNA binding in the termolecular system of M. Hpy N6mA. In both systems, AdoMet is stabilized in the binding pocket by hydrogen bonding (Thr84, Glu99, Asp122, and Phe123) as well as van der Waals (Ile100, Phe160, Arg104, and Cys76) interactions. We propose that the contacts made by flipped adenine DA6 with Asn138 (N6 and N1 atom of DA6) and Pro139 (N6) and π-stacking interactions with Phe141 and Phe219 play an important role in the methylation mechanism at the N6 position in our N6mA model. Specific recognition of DNA is mediated by residues 143-155, 183-189, 212-220, 280-293, and 308-325. These findings are further supported by alanine scanning mutagenesis studies. The conserved residues in distantly related sequences of the small domain are important in DNA binding. Results reported here elucidate the sequence, structure, and binding features necessary for the recognition between cofactor AdoMet and substrate DNA by the vital epigenetic enzyme, M. Hpy N6mA.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

DNA Methyltransferase Activity Assays: Advances and Challenges

DNA methyltransferases (MTases), a family of enzymes that catalyse the methylation of DNA, have a profound effect on gene regulation. A large body of evidence has indicated that DNA MTase is potentially a predictive biomarker closely associated with genetic disorders and genetic diseases like cancer. Given the attention bestowed onto DNA MTases in molecular biology and medicine, highly sensitive detection of DNA MTase activity is essential in determining gene regulation, epigenetic modification, clinical diagnosis and therapeutics. Conventional techniques such as isotope labelling are effective, but they often require laborious sample preparation, isotope labelling, sophisticated equipment and large amounts of DNA, rendering them unsuitable for uses at point-of-care. Simple, portable, highly sensitive and low-cost assays are urgently needed for DNA MTase activity screening. In most recent technological advances, many alternative DNA MTase activity assays such as fluorescent, electrochemical, colorimetric and chemiluminescent assays have been proposed. In addition, many of them are coupled with nanomaterials and/or enzymes to significantly enhance their sensitivity. Herein we review the progress in the development of DNA MTase activity assays with an emphasis on assay mechanism and performance with some discussion on challenges and perspectives. It is hoped that this article will provide a broad coverage of DNA MTase activity assays and their latest developments and open new perspectives toward the development of DNA MTase activity assays with much improved performance for uses in molecular biology and clinical practice.

DNA Base Flipping: A General Mechanism for Writing, Reading, and Erasing DNA Modifications

The modification of DNA bases is a classic hallmark of epigenetics. Four forms of modified cytosine-5-methylcytosine, 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine-have been discovered in eukaryotic DNA. In addition to cytosine carbon-5 modifications, cytosine and adenine methylated in the exocyclic amine-N4-methylcytosine and N6-methyladenine-are other modified DNA bases discovered even earlier. Each modified base can be considered a distinct epigenetic signal with broader biological implications beyond simple chemical changes. Since 1994, crystal structures of proteins and enzymes involved in writing, reading, and erasing modified bases have become available. Here, we present a structural synopsis of writers, readers, and erasers of the modified bases from prokaryotes and eukaryotes. Despite significant differences in structures and functions, they are remarkably similar regarding their engagement in flipping a target base/nucleotide within DNA for specific recognitions and/or reactions. We thus highlight base flipping as a common structural framework broadly applied by distinct classes of proteins and enzymes across phyla for epigenetic regulations of DNA.

Adenine methylation in eukaryotes: Apprehending the complex evolutionary history and functional potential of an epigenetic modification

While N⁶‐methyladenosine (m⁶A) is a well‐known epigenetic modification in bacterial DNA, it remained largely unstudied in eukaryotes. Recent studies have brought to fore its potential epigenetic role across diverse eukaryotes with biological consequences, which are distinct and possibly even opposite to the well‐studied 5‐methylcytosine mark. Adenine methyltransferases appear to have been independently acquired by eukaryotes on at least 13 occasions from prokaryotic restriction‐modification and counter‐restriction systems. On at least four to five instances, these methyltransferases were recruited as RNA methylases. Thus, m⁶A marks in eukaryotic DNA and RNA might be more widespread and diversified than previously believed. Several m⁶A‐binding protein domains from prokaryotes were also acquired by eukaryotes, facilitating prediction of potential readers for these marks. Further, multiple lineages of the AlkB family of dioxygenases have been recruited as m⁶A demethylases. Although members of the TET/JBP family of dioxygenases have also been suggested to be m⁶A demethylases, this proposal needs more careful evaluation.

Also watch the Video Abstract.

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

RNA/DNA hybrid duplexes regularly occur in nature, for example in transcriptional R loops. Their susceptibility to modification by DNA-specific or RNA-specific enzymes is, thus, a biologically relevant question, which, in addition, has possible biotechnological implications. In this study, we investigated the activity of four isospecific DNA methyltransferases (M.EcoVIII, M.LlaCI, M.HindIII, M.BstZ1II) toward an RNA/DNA duplex carrying one 5'-AAGCUU-3'/3'-TTCGAA-5' target sequence. The analyzed enzymes belong to the β-group of adenine N6-methyltransferases and recognize the palindromic DNA sequence 5'-AAGCTT-3'/3'-TTCGAA-5'. Under standard conditions, none of these isospecific enzymes could detectibly methylate the RNA/DNA duplex. However, the addition of agents that generally relax specificity, such as dimethyl sulfoxide (DMSO) and glycerol, resulted in substantial methylation of the RNA/DNA duplex by M.EcoVIII and M.LlaCI. Only the DNA strand of the RNA/DNA duplex was methylated. The same was not observed for M.HindIII or M.BstZ1II. This is, to our knowledge, the first report that demonstrates such activity by prokaryotic DNA methyltransferases. Possible applications of these findings in a laboratory practice are also discussed.

Naturally-occurring, dually-functional fusions between restriction endonucleases and regulatory proteins

Background

Restriction-modification (RM) systems appear to play key roles in modulating gene flow among bacteria and archaea. Because the restriction endonuclease (REase) is potentially lethal to unmethylated new host cells, regulation to ensure pre-expression of the protective DNA methyltransferase (MTase) is essential to the spread of RM genes. This is particularly true for Type IIP RM systems, in which the REase and MTase are separate, independently-active proteins. A substantial subset of Type IIP RM systems are controlled by an activator-repressor called C protein. In these systems, C controls the promoter for its own gene, and for the downstream REase gene that lacks its own promoter. Thus MTase is expressed immediately after the RM genes enter a new cell, while expression of REase is delayed until sufficient C protein accumulates. To study the variation in and evolution of this regulatory mechanism, we searched for RM systems closely related to the well-studied C protein-dependent PvuII RM system. Unexpectedly, among those found were several in which the C protein and REase genes were fused.

Results

The gene for CR.NsoJS138I fusion protein (nsoJS138ICR, from the bacterium Niabella soli) was cloned, and the fusion protein produced and partially purified. Western blots provided no evidence that, under the conditions tested, anything other than full-length fusion protein is produced. This protein had REase activity in vitro and, as expected from the sequence similarity, its specificity was indistinguishable from that for PvuII REase, though the optimal reaction conditions were different. Furthermore, the fusion was active as a C protein, as revealed by in vivo activation of a lacZ reporter fusion to the promoter region for the nsoJS138ICR gene.

Conclusions

Fusions between C proteins and REases have not previously been characterized, though other fusions have (such as between REases and MTases). These results reinforce the evidence for impressive modularity among RM system proteins, and raise important questions about the implications of the C-REase fusions on expression kinetics of these RM systems.

Recognition of guanosine by dissimilar tRNA methyltransferases

Guanosines are important for biological activities through their specific functional groups that are recognized for RNA or protein interactions. One example is recognition of N(1) of G37 in tRNA by S-adenosyl-methionine (AdoMet)-dependent tRNA methyltransferases to synthesize m(1)G37-tRNA, which is essential for translational fidelity in all biological domains. Synthesis of m(1)G37-tRNA is catalyzed by TrmD in bacteria and by Trm5 in eukarya and archaea, using unrelated and dissimilar structural folds. This raises the question of how dissimilar proteins recognize the same guanosine. Here we probe the mechanism of discrimination among functional groups of guanosine by TrmD and Trm5. Guanosine analogs were systematically introduced into tRNA through a combination of chemical and enzymatic synthesis. Single turnover kinetic assays and thermodynamic analysis of the effect of each analog on m(1)G37-tRNA synthesis reveal that TrmD and Trm5 discriminate functional groups differently. While both recognize N(1) and O(6) of G37, TrmD places a much stronger emphasis on these functional groups than Trm5. While the exocyclic 2-amino group of G37 is important for TrmD, it is dispensable for Trm5. In addition, while an adjacent G36 is obligatory for TrmD, it is nonessential for Trm5. These results depict a more rigid requirement of guanosine functional groups for TrmD than for Trm5. However, the sensitivity of both enzymes to analog substitutions, together with an experimental revelation of their low cellular concentrations relative to tRNA substrates, suggests a model in which these enzymes rapidly screen tRNA by direct recognition of G37 in order to monitor the global state of m(1)G37-tRNA.

Overcoming sequence misalignments with weighted structural superposition

An appropriate structural superposition identifies similarities and differences between homologous proteins that are not evident from sequence alignments alone. We have coupled our Gaussian-weighted RMSD (wRMSD) tool with a sequence aligner and seed extension (SE) algorithm to create a robust technique for overlaying structures and aligning sequences of homologous proteins (HwRMSD). HwRMSD overcomes errors in the initial sequence alignment that would normally propagate into a standard RMSD overlay. SE can generate a corrected sequence alignment from the improved structural superposition obtained by wRMSD. HwRMSD's robust performance and its superiority over standard RMSD are demonstrated over a range of homologous proteins. Its better overlay results in corrected sequence alignments with good agreement to HOMSTRAD. Finally, HwRMSD is compared to established structural alignment methods: FATCAT, secondary-structure matching, combinatorial extension, and Dalilite. Most methods are comparable at placing residue pairs within 2 Å, but HwRMSD places many more residue pairs within 1 Å, providing a clear advantage. Such high accuracy is essential in drug design, where small distances can have a large impact on computational predictions. This level of accuracy is also needed to correct sequence alignments in an automated fashion, especially for omics-scale analysis. HwRMSD can align homologs with low-sequence identity and large conformational differences, cases where both sequence-based and structural-based methods may fail. The HwRMSD pipeline overcomes the dependency of structural overlays on initial sequence pairing and removes the need to determine the best sequence-alignment method, substitution matrix, and gap parameters for each unique pair of homologs.

Substrate promiscuity in DNA methyltransferase M.PvuII. A mechanistic insight

M.PvuII is a DNA methyltransferase from the bacterium Proteus vulgaris that catalyzes methylation of cytosine at the N4 position. This enzyme also displays promiscuous activity catalyzing methylation of adenine at the N6 position. In this work we use QM/MM methods to investigate the reaction mechanism of this promiscuous activity. We found that N6 methylation in M.PvuII takes place by means of a stepwise mechanism in which deprotonation of the exocyclic amino group is followed by the methyl transfer. Deprotonation involves two residues of the active site, Ser53 and Asp96, while methylation takes place directly from the AdoMet cofactor to the target nitrogen atom. The same reaction mechanism was described for cytosine methylation in the same enzyme, while the reversal timing, that is methylation followed by deprotonation, has been described in M.TaqI, an enzyme that catalyzes the N6-adenine DNA methylation from Thermus aquaticus. These mechanistic findings can be useful to understand the evolutionary paths followed by N-methyltransferases.

DNA methyltransferases: mechanistic models derived from kinetic analysis

The sequence-specific transfer of methyl groups from donor S-adenosyl-L-methionine (AdoMet) to certain positions of DNA-adenine or -cytosine residues by DNA methyltransferases (MTases) is a major form of epigenetic modification. It is virtually ubiquitous, except for some notable exceptions. Site-specific methylation can be regarded as a means to increase DNA information capacity and is involved in a large spectrum of biological processes. The importance of these functions necessitates a deeper understanding of the enzymatic mechanism(s) of DNA methylation. DNA MTases fall into one of two general classes; viz. amino-MTases and [C5-cytosine]-MTases. Amino-MTases, common in prokaryotes and lower eukaryotes, catalyze methylation of the exocyclic amino group of adenine ([N6-adenine]-MTase) or cytosine ([N4-cytosine]-MTase). In contrast, [C5-cytosine]-MTases methylate the cyclic carbon-5 atom of cytosine. Characteristics of DNA MTases are highly variable, differing in their affinity to their substrates or reaction products, their kinetic parameters, or other characteristics (order of substrate binding, rate limiting step in the overall reaction). It is not possible to present a unifying account of the published kinetic analyses of DNA methylation because different authors have used different substrate DNAs and/or reaction conditions. Nevertheless, it would be useful to describe those kinetic data and the mechanistic models that have been derived from them. Thus, this review considers in turn studies carried out with the most consistently and extensively investigated [N6-adenine]-, [N4-cytosine]- and [C5-cytosine]-DNA MTases.

Methylation of class I translation termination factors: structural and functional aspects

During protein synthesis, release of polypeptide from the ribosome occurs when an in frame termination codon is encountered. Contrary to sense codons, which are decoded by tRNAs, stop codons present in the A-site are recognized by proteins named class I release factors, leading to the release of newly synthesized proteins. Structures of these factors bound to termination ribosomal complexes have recently been obtained, and lead to a better understanding of stop codon recognition and its coordination with peptidyl-tRNA hydrolysis in bacteria. Release factors contain a universally conserved GGQ motif which interacts with the peptidyl-transferase centre to allow peptide release. The Gln side chain from this motif is methylated, a feature conserved from bacteria to man, suggesting an important biological role. However, methylation is catalysed by completely unrelated enzymes. The function of this motif and its post-translational modification will be discussed in the context of recent structural and functional studies.

Cloning and characterization of the TneDI restriction: modification system of Thermotoga neapolitana

A putative Type II restriction-modification system of Thermotoga neapolitana, TneDI, was cloned into Escherichia coli XL1-Blue MRF' and characterized. Gene CTN_0339 specifies the endonuclease R.TneDI, while CTN_0340 encodes the cognate DNA methyltransferase M.TneDI. Both enzymes were purified simply by heating the cell lysates of E. coli followed by centrifugation. The enzymes were active over a broad range of temperatures, from 42°C to at least 77°C, with the highest activities observed at 77°C. R.TneDI cleaved at the center of the recognition sequence (CG↓CG) and generated blunt-end cuts. Overexpression of R.TneDI in BL21(DE3) was confirmed by both SDS-PAGE and Western blotting.

Structural and functional characterization of Rv2966c protein reveals an RsmD-like methyltransferase from Mycobacterium tuberculosis and the role of its N-terminal domain in target recognition

Nine of ten methylated nucleotides of Escherichia coli 16 S rRNA are conserved in Mycobacterium tuberculosis. All the 10 different methyltransferases are known in E. coli, whereas only TlyA and GidB have been identified in mycobacteria. Here we have identified Rv2966c of M. tuberculosis as an ortholog of RsmD protein of E. coli. We have shown that rv2966c can complement rsmD-deleted E. coli cells. Recombinant Rv2966c can use 30 S ribosomes purified from rsmD-deleted E. coli as substrate and methylate G966 of 16 S rRNA in vitro. Structure determination of the protein shows the protein to be a two-domain structure with a short hairpin domain at the N terminus and a C-terminal domain with the S-adenosylmethionine-MT-fold. We show that the N-terminal hairpin is a minimalist functional domain that helps Rv2966c in target recognition. Deletion of the N-terminal domain prevents binding to nucleic acid substrates, and the truncated protein fails to carry out the m(2)G966 methylation on 16 S rRNA. The N-terminal domain also binds DNA efficiently, a property that may be utilized under specific conditions of cellular growth.

Natural history of eukaryotic DNA methylation systems

Methylation of cytosines and adenines in DNA is a widespread epigenetic mark in both prokaryotes and eukaryotes. In eukaryotes, it has a profound influence on chromatin structure and dynamics. Recent advances in genomics and biochemistry have considerably elucidated the functions and provenance of these DNA modifications. DNA methylases appear to have emerged first in bacterial restriction-modification (R-M) systems from ancient RNA-modifying enzymes, in transitions that involved acquisition of novel catalytic residues and DNA-recognition features. DNA adenine methylases appear to have been acquired by ciliates, heterolobosean amoeboflagellates, and certain chlorophyte algae. Six distinct clades of cytosine methylases, including the DNMT1, DNMT2, and DNMT3 clades, were acquired by eukaryotes through independent lateral transfer of their precursors from bacteria or bacteriophages. In addition to these, multiple adenine and cytosine methylases were acquired by several families of eukaryotic transposons. In eukaryotes, the DNA-methylase module was often combined with distinct modified and unmodified peptide recognition domains and other modules mediating specialized interactions, for example, the RFD module of DNMT1 which contains a permuted Sm domain linked to a helix-turn-helix domain. In eukaryotes, the evolution of DNA methylases appears to have proceeded in parallel to the elaboration of histone-modifying enzymes and the RNAi system, with functions related to counter-viral and counter-transposon defense, and regulation of DNA repair and differential gene expression being their primary ancestral functions. Diverse DNA demethylation systems that utilize base-excision repair via DNA glycosylases and cytosine deaminases appear to have emerged in multiple eukaryotic lineages. Comparative genomics suggests that the link between cytosine methylation and DNA glycosylases probably emerged first in a novel R-M system in bacteria. Recent studies suggest that the 5mC is not a terminal DNA modification, with enzymes of the Tet/JBP family of 2-oxoglutarate- and iron-dependent dioxygenases further hydroxylating it to form 5-hydroxymethylcytosine (5hmC). These enzymes emerged first in bacteriophages and appear to have been transferred to eukaryotes on one or more occasions. Eukaryotes appear to have recruited three major types of DNA-binding domains (SRA/SAD, TAM/MBD, and CXXC) in discriminating DNA with methylated or unmethylated cytosines. Analysis of the domain architectures of these domains and the DNA methylases suggests that early in eukaryotic evolution they developed a close functional link with SET-domain methylases and Jumonji-related demethylases that operate on peptides in chromatin proteins. In several eukaryotes, other functional connections were elaborated in the form of various combinations between domains related to DNA methylation and those involved in ATP-dependent chromatin remodeling and RNAi. In certain eukaryotes, such as mammals and angiosperms, novel dependencies on the DNA methylation system emerged, which resulted in it affecting unexpected aspects of the biology of these organisms such as parent-offspring interactions. In genomic terms, this was reflected in the emergence of new proteins related to methylation, such as Stella. The well-developed methylation systems of certain heteroloboseans, stramenopiles, chlorophytes, and haptophyte indicate that these might be new model systems to explore the relevance of DNA modifications in eukaryotes.

Plasmid addiction systems: perspectives and applications in biotechnology

Biotechnical production processes often operate with plasmid-based expression systems in well-established prokaryotic and eukaryotic hosts such as Escherichia coli or Saccharomyces cerevisiae, respectively. Genetically engineered organisms produce important chemicals, biopolymers, biofuels and high-value proteins like insulin. In those bioprocesses plasmids in recombinant hosts have an essential impact on productivity. Plasmid-free cells lead to losses in the entire product recovery and decrease the profitability of the whole process. Use of antibiotics in industrial fermentations is not an applicable option to maintain plasmid stability. Especially in pharmaceutical or GMP-based fermentation processes, deployed antibiotics must be inactivated and removed. Several plasmid addiction systems (PAS) were described in the literature. However, not every system has reached a full applicable state. This review compares most known addiction systems and is focusing on biotechnical applications.

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

In this paper the reaction mechanism for methylation of cytosine at the exocyclic N4 position catalyzed by M.PvuII has been explored by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods. A reaction model was prepared by placing a single cytosine base in the active site of the enzyme. In this model the exocyclic amino group of the base establishes hydrogen bond interactions with the hydroxyl oxygen atom of Ser53 and the carbonyl oxygen atom of Pro54. The reaction mechanism involves a direct methyl transfer from AdoMet to the N4 atom and a proton transfer from this atom to Ser53, which in turn transfers a proton to Asp96. Different timings for the proton transfers and methylation steps have been explored at the AM1/MM and B3LYP/MM levels including localization and characterization of stationary structures. At our best estimate the reaction proceeds by means of a simultaneous but asynchronous proton transfer from Ser53 to Asp96 and from N4 of cytosine to Ser53 followed by a direct methyl transfer from AdoMet to the exocyclic N4 of cytosine.

Flavivirus RNA cap methyltransferase: structure, function, and inhibition

Many flaviviruses are significant human pathogens. The plus-strand RNA genome of a flavivirus contains a 5' terminal cap 1 structure (m(7)GpppAmG). The flavivirus encodes one methyltransferase (MTase), located at the N-terminal portion of the NS5 RNA-dependent RNA polymerase (RdRp). Here we review recent advances in our understanding of flaviviral capping machinery and the implications for drug development. The NS5 MTase catalyzes both guanine N7 and ribose 2'-OH methylations during viral cap formation. Representative flavivirus MTases, from dengue, yellow fever, and West Nile virus (WNV), sequentially generate GpppA → m(7)GpppA → m(7)GpppAm. Despite the existence of two distinct methylation activities, the crystal structures of flavivirus MTases showed a single binding site for S-adenosyl-L-methionine (SAM), the methyl donor. This finding indicates that the substrate GpppA-RNA must be repositioned to accept the N7 and 2'-O methyl groups from SAM during the sequential reactions. Further studies demonstrated that distinct RNA elements are required for the methylations of guanine N7 on the cap and of ribose 2'-OH on the first transcribed nucleotide. Mutant enzymes with different methylation defects can trans complement one another in vitro, demonstrating that separate molecules of the enzyme can independently catalyze the two cap methylations in vitro. In the context of the infectious virus, defects in both methylations, or a defect in the N7 methylation alone, are lethal to WNV. However, viruses defective solely in 2'-O methylation are attenuated and can protect mice from later wild-type WNV challenge. The results demonstrate that the N7 methylation activity is essential for the WNV life cycle and, thus, methyltransferase represents a novel and promising target for flavivirus therapy.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

Multiple Motif Scanning to identify methyltransferases from the yeast proteome

A new program (Multiple Motif Scanning) was developed to scan the Saccharomyces cerevisiae proteome for Class I S-adenosylmethionine-dependent methyltransferases. Conserved Motifs I, Post I, II, and III were identified and expanded in known methyltransferases by primary sequence and secondary structural analysis through hidden Markov model profiling of both a yeast reference database and a reference database of methyltransferases with solved three-dimensional structures. The roles of the conserved amino acids in the four motifs of the methyltransferase structure and function were then analyzed to expand the previously defined motifs. Fisher-based negative log statistical matrix sets were developed from the prevalence of amino acids in the motifs. Multiple Motif Scanning is able to scan the proteome and score different combinations of the top fitting sequences for each motif. In addition, the program takes into account the conserved number of amino acids between the motifs. The output of the program is a ranked list of proteins that can be used to identify new methyltransferases and to reevaluate the assignment of previously identified putative methyltransferases. The Multiple Motif Scanning program can be used to develop a putative list of enzymes for any type of protein that has one or more motifs conserved at variable spacings and is freely available (www.chem.ucla.edu/files/MotifSetup.Zip). Finally hidden Markov model profile clustering analysis was used to subgroup Class I methyltransferases into groups that reflect their methyl-accepting substrate specificity.

QM/MM hybrid calculation of biological macromolecules using a new interface program connecting QM and MM engines

An interface program connecting a quantum mechanics (QM) calculation engine, GAMESS, and a molecular mechanics (MM) calculation engine, AMBER, has been developed for QM/MM hybrid calculations. A protein-DNA complex is used as a test system to investigate the following two types of QM/MM schemes. In a 'subtractive' scheme, electrostatic interactions between QM/MM regions are truncated in QM calculations; in an 'additive' scheme, long-range electrostatic interactions within a cut-off distance from QM regions are introduced into one-electron integration terms of a QM Hamiltonian. In these calculations, 338 atoms are assigned as QM atoms using Hartree-Fock (HF)/density functional theory (DFT) hybrid all-electron calculations. By comparing the results of the additive and subtractive schemes, it is found that electronic structures are perturbed significantly by the introduction of MM partial charges surrounding QM regions, suggesting that biological processes occurring in functional sites are modulated by the surrounding structures. This also indicates that the effects of long-range electrostatic interactions involved in the QM Hamiltonian are crucial for accurate descriptions of electronic structures of biological macromolecules.

Tuning the relative affinities for activating and repressing operators of a temporally regulated restriction-modification system

Most type II restriction-modification (R-M) systems produce separate endonuclease (REase) and methyltransferase (MTase) proteins. After R-M genes enter a new cell, MTase activity must appear before REase or the host chromosome will be cleaved. Temporal control of these genes thus has life-or-death consequences. PvuII and some other R-M systems delay endonuclease expression by cotranscribing the REase gene with the upstream gene for an autogenous activator/repressor (C protein). C.PvuII was previously shown to have low levels early, but positive feedback later boosts transcription of the C and REase genes. The MTase is expressed without delay, and protects the host DNA. C.PvuII binds to two sites upstream of its gene: O_L, associated with activation, and O_R, associated with repression. Even when symmetry elements of each operator are made identical, C.PvuII binds preferentially to O_L. In this study, the intra-operator spacers are shown to modulate relative C.PvuII affinity. In light of a recently reported C.Esp1396I-DNA co-crystal structure, in vitro and in vivo effects of altering O_L and O_R spacers were determined. The results suggest that the GACTnnnAGTC consensus is the primary determinant of C.PvuII binding affinity, with intra-operator spacers playing a fine-tuning role that affects mobility of this R-M system.

Crystal structure of tRNA N2,N2-guanosine dimethyltransferase Trm1 from Pyrococcus horikoshii

Trm1 catalyzes a two-step reaction, leading to mono- and dimethylation of guanosine at position 26 in most eukaryotic and archaeal tRNAs. We report the crystal structures of Trm1 from Pyrococcus horikoshii liganded with S-adenosyl-l-methionine or S-adenosyl-l-homocysteine. The protein comprises N-terminal and C-terminal domains with class I methyltransferase and novel folds, respectively. The methyl moiety of S-adenosyl-l-methionine points toward the invariant Phe27 and Phe140 within a narrow pocket, where the target G26 might flip in. Mutagenesis of Phe27 or Phe140 to alanine abolished the enzyme activity, indicating their role in methylating G26. Structural analyses revealed that the movements of Phe140 and the loop preceding Phe27 may be involved in dissociation of the monomethylated tRNA*Trm1 complex prior to the second methylation. Moreover, the catalytic residues Asp138, Pro139, and Phe140 are in a different motif from that in DNA 6-methyladenosine methyltransferases, suggesting a different methyl transfer mechanism in the Trm1 family.

DNA base flipping by both members of the PspGI restriction-modification system

The PspGI restriction–modification system recognizes the sequence CCWGG. R.PspGI cuts DNA before the first C in the cognate sequence and M.PspGI is thought to methylate N4 of one of the cytosines in the sequence. M.PspGI enhances fluorescence of 2-aminopurine in DNA if it replaces the second C in the sequence, while R.PspGI enhances fluorescence when the fluorophore replaces adenine in the central base pair. This strongly suggests that the methyltransferase flips the second C in the recognition sequence, while the endonuclease flips both bases in the central base pair out of the duplex. M.PspGI is the first N4-cytosine MTase for which biochemical evidence for base flipping has been presented. It is also the first type IIP methyltransferase whose catalytic activity is strongly stimulated by divalent metal ions. However, divalent metal ions are not required for its base-flipping activity. In contrast, these ions are required for both base flipping and catalysis by the endonuclease. The two enzymes have similar temperature profiles for base flipping and optimal flipping occurs at temperatures substantially below the growth temperature of the source organism for PspGI and for the catalytic activity of endonuclease. We discuss the implications of these results for DNA binding by these enzymes and their evolutionary origin.

A comprehensive analysis of non-sequential alignments between all protein structures

Background

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.

Results

The study has been conducted on a non-redundant set of 8,865 protein structures aligned with the aid of the TOPOFIT method. It has been estimated that between 17.4% and 35.2% of all alignments are non-sequential depending on variations in the parameters. Analysis of the data revealed that non-sequential relations between proteins do occur systematically and in large quantities. Various sizes and numbers of non-sequential fragments have been observed with all possible complexities of fragment rearrangements found for alignments consisting of up to 12 fragments. It has been found that non-sequential alignments are not limited to proteins of any particular fold and are present in more than two hundred of them. Moreover, many of them are found between proteins with different fold assignments. It has been shown that protein structure symmetry does not explain non-sequential alignments. Therefore, compelling evidences have been provided that non-sequential alignments between proteins are systematic and widespread across the protein universe.

Conclusion

The phenomenon of the widespread occurrence of non-sequential alignments between proteins might represent a missing rule of protein structure organization. More detailed study of this phenomenon will enhance our understanding of protein stability, folding, and evolution.

Regulatory circuit based on autogenous activation-repression: roles of C-boxes and spacer sequences in control of the PvuII restriction-modification system

Type II restriction-modification (R-M) systems comprise a restriction endonuclease (REase) and a protective methyltransferase (MTase). After R-M genes enter a new cell, MTase must appear before REase or the chromosome will be cleaved. PvuII and some other R-M systems achieve this delay by cotranscribing the REase gene with the gene for an autogenous transcription activator (the controlling or 'C' protein C.PvuII). This study reveals, through in vivo titration, that C.PvuII is not only an activator but also a repressor for its own gene. In other systems, this type of circuit can result in oscillatory behavior. Despite the use of identical, symmetrical C protein-binding sequences (C-boxes) in the left and right operators, C.PvuII showed higher in vitro affinity for O(L) than for O(R), implicating the spacer sequences in this difference. Mutational analysis associated the repression with O(R), which overlaps the promoter -35 hexamer but is otherwise dispensable for activation. A nonrepressing mutant exhibited poor establishment in new cells. Comparing promoter-operator regions from PvuII and 29 R-M systems controlled by C proteins revealed that the most-highly conserved sequence is the tetranucleotide spacer separating O(L) from O(R). Any changes in that spacer reduced the stability of C.PvuII-operator complexes and abolished activation.

Convergent Evolution of Enzyme Active Sites Is not a Rare Phenomenon

Since convergent evolution of enzyme active sites was first identified in serine proteases, other individual instances of this phenomenon have been documented. However, a systematic analysis assessing the frequency of this phenomenon across enzyme space is still lacking. This work uses the Query3d structural comparison algorithm to integrate for the first time detailed knowledge about catalytic residues, available through the Catalytic Site Atlas (CSA), with the evolutionary information provided by the Structural Classification of Proteins (SCOP) database. This study considers two modes of convergent evolution: (i) mechanistic analogues which are enzymes that use the same mechanism to perform related, but possibly different, reactions (considered here as sharing the first three digits of the EC number); and (ii) transformational analogues which catalyse exactly the same reaction (identical EC numbers), but may use different mechanisms. Mechanistic analogues were identified in 15% (26 out of 169) of the three-digit EC groups considered, showing that this phenomenon is not rare. Furthermore 11 of these groups also contain transformational analogues. The catalytic triad is the most widespread active site; the results of the structural comparison show that this mechanism, or variations thereof, is present in 23 superfamilies. Transformational analogues were identified for 45 of the 951 four-digit EC numbers present within the CSA and about half of these were also mechanistic analogues exhibiting convergence of their active sites. This analysis has also been extended to the whole Protein Data Bank to provide a complete and manually curated list of the all the transformational analogues whose structure is classified in SCOP. The results of this work show that the phenomenon of convergent evolution is not rare, especially when considering large enzymatic families.

1 Protein Methyltransferases: Their Distribution Among the Five Structural Classes of AdoMet-Dependent Methyltransferases

S-adenosyl-l-methionine (AdoMet) dependent methyltransferases (MTases) are involved in biosynthesis, signal transduction, protein repair, chromatin regulation, and gene silencing. Five different structural folds (designated I through V) have been described that bind AdoMet and catalyze methyltransfer to diverse substrates, although the great majority of known MTases have the Class I fold. Even within a particular MTase class the amino-acid sequence similarity can be as low as 10%. Thus, the structural and catalytic requirements for methyltransfer from AdoMet appear to be remarkably flexible. MTases that act on protein substrates have been found to date among three of the five structural classes (I, the classical fold; III, the corrin MTase fold; and V, the SET fold). "There are many paths to the top of the mountain, but the view is always the same."-Chinese proverb The Columbia World of Quotations, New York, Columbia University Press, 1996.

A rapid and efficient method for cloning genes of type II restriction-modification systems by use of a killer plasmid

We present a method for cloning restriction-modification (R-M) systems that is based on the use of a lethal plasmid (pKILLER). The plasmid carries a functional gene for a restriction endonuclease having the same DNA specificity as the R-M system of interest. The first step is the standard preparation of a representative, plasmid-borne genomic library. Then this library is transformed with the killer plasmid. The only surviving bacteria are those which carry the gene specifying a protective DNA methyltransferase. Conceptually, this in vivo selection approach resembles earlier methods in which a plasmid library was selected in vitro by digestion with a suitable restriction endonuclease, but it is much more efficient than those methods. The new method was successfully used to clone two R-M systems, BstZ1II from Bacillus stearothermophilus 14P and Csp231I from Citrobacter sp. strain RFL231, both isospecific to the prototype HindIII R-M system.

Protein interfacial pocket engineering via coupled computational filtering and biological focusing criterion

To engineer bio-macromolecular systems, protein-substrate interactions and their configurations need to be understood, harnessed, and utilized. Due to the inherent large numbers of combinatorial configurations and conformational complexity, methods that rely on heuristics or stochastics, such as practical computational filtering (CF) or biological focusing (BF) criterions, when used alone rarely yield insights into these complexes or successes in (re)designing them. Here we use a coupled CF-BF criterion upon an amenable interfacial pocket (IP) of a protein scaffold complexed with its substrate to undergo residue replacement and R-group refinement (R4) to filter out energetically unfavorable residues and R-group conformations, and focus in on those that are evolutionarily favorable. We show that this coupled filtering and focusing can efficiently provide a putative engineered IP candidate and validate it computationally and empirically. The CF-BF criterion may permit holistic understanding of the nuances of existing protein IPs and their scaffolds and facilitate bioengineering efforts to alter substrate specificity. Such approach may contribute to accelerated elucidation of engineering principles of bio-macromolecular systems.

Topology of Type II REases revisited; structural classes and the common conserved core

Type II restriction endonucleases (REases) are deoxyribonucleases that cleave DNA sequences with remarkable specificity. Type II REases are highly divergent in sequence as well as in topology, i.e. the connectivity of secondary structure elements. A widely held assumption is that a structural core of five β-strands flanked by two α-helices is common to these enzymes. We introduce a systematic procedure to enumerate secondary structure elements in an unambiguous and reproducible way, and use it to analyze the currently available X-ray structures of Type II REases. Based on this analysis, we propose an alternative definition of the core, which we term the αβα-core. The αβα-core includes the most frequently observed secondary structure elements and is not a sandwich, as it consists of a five-strand β-sheet and two α-helices on the same face of the β-sheet. We use the αβα-core connectivity as a basis for grouping the Type II REases into distinct structural classes. In these new structural classes, the connectivity correlates with the angles between the secondary structure elements and with the cleavage patterns of the REases. We show that there exists a substructure of the αβα-core, namely a common conserved core, ccc, defined here as one α-helix and four β-strands common to all Type II REase of known structure.