Structure prediction and phylogenetic analysis of a functionally diverse family of proteins homologous to the MT-A70 subunit of the human mRNA:m(6)A methyltransferase

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

Is the HemK family of putative S-adenosylmethionine-dependent methyltransferases a "missing" zeta subfamily of adenine methyltransferases? A hypothesis

Previous comparative studies revealed close similarity among various groups of S-adenosyl-L-methionine (AdoMet)-dependent methyltransferases (MTases), indicating their common evolutionary origin. We present evidence for a remarkable similarity between the sequence and predicted structure of HemK (a widespread family of putative proteins encoded in genomes from bacteria to humans) and the catalytic domain of the gamma-subfamily of adenine-specific DNA MTases (N6mA MTases). We predict the structure and function of the putative catalytic domain of HemK proteins and speculate that the target-recognizing function may be conferred by the N-terminal variable region.     

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

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

Statistical correlation between protein secondary structure and messenger RNA stem-loop structure

A new integrated sequence-structure database, called IADE (Integrated ASTRAL-DSSP-EMBL), incorporating matching mRNA sequence, amino acid sequence, and protein secondary structural data, is constructed. It includes 648 protein domains. Based on the IADE database, we studied the relation between RNA stem-loop frequencies and protein secondary structure. It was found that the alpha-helices and beta-strands on proteins tend to be preferably "coded" by mRNA stem region, while the coils on proteins tend to be preferably "coded" by mRNA loop region. These tendencies are more obvious if we observe the structural words (SWs). An SW is defined by a four-amino-acid-fragment that shows the pronounced secondary structural (alpha-helix or beta-strand) propensity. It is demonstrated that the deduced correlation between protein and mRNA structure can hardly be explained as the stochastic fluctuation effect.     

Structure of human DNMT2, an enigmatic DNA methyltransferase homolog that displays denaturant-resistant binding to DNA

DNMT2 is a human protein that displays strong sequence similarities to DNA (cytosine-5)-methyltransferases (m(5)C MTases) of both prokaryotes and eukaryotes. DNMT2 contains all 10 sequence motifs that are conserved among m(5)C MTases, including the consensus S:-adenosyl-L-methionine-binding motifs and the active site ProCys dipeptide. DNMT2 has close homologs in plants, insects and Schizosaccharomyces pombe, but no related sequence can be found in the genomes of Saccharomyces cerevisiae or Caenorhabditis elegans. The crystal structure of a deletion mutant of DNMT2 complexed with S-adenosyl-L-homocysteine (AdoHcy) has been determined at 1.8 A resolution. The structure of the large domain that contains the sequence motifs involved in catalysis is remarkably similar to that of M.HHAI, a confirmed bacterial m(5)C MTase, and the smaller target recognition domains of DNMT2 and M.HHAI are also closely related in overall structure. The small domain of DNMT2 contains three short helices that are not present in M.HHAI. DNMT2 binds AdoHcy in the same conformation as confirmed m(5)C MTases and, while DNMT2 shares all sequence and structural features with m(5)C MTases, it has failed to demonstrate detectable transmethylase activity. We show here that homologs of DNMT2, which are present in some organisms that are not known to methylate their genomes, contain a specific target-recognizing sequence motif including an invariant CysPheThr tripeptide. DNMT2 binds DNA to form a denaturant-resistant complex in vitro. While the biological function of DNMT2 is not yet known, the strong binding to DNA suggests that DNMT2 may mark specific sequences in the genome by binding to DNA through the specific target-recognizing motif.     

mRNA:guanine-N7 cap methyltransferases: identification of novel members of the family, evolutionary analysis, homology modeling, and analysis of sequence-structure-function relationships

Background  

The 5'-terminal cap structure plays an important role in many aspects of mRNA metabolism. Capping enzymes encoded by viruses and pathogenic fungi are attractive targets for specific inhibitors. There is a large body of experimental data on viral and cellular methyltransferases (MTases) that carry out guanine-N7 (cap 0) methylation, including results of extensive mutagenesis. However, a crystal structure is not available and cap 0 MTases are too diverged from other MTases of known structure to allow straightforward homology-based interpretation of these data.     

Nucleotide sequence analysis of the gene specifying the bifunctional 6''-aminoglycoside acetyltransferase 2"-aminoglycoside phosphotransferase enzyme in Streptococcus faecalis and identification and cloning of gene regions specifying the two activities.

The gene specifying the bifunctional 6'-aminoglycoside acetyltransferase [AAC(6')] 2"-aminoglycoside phosphotransferase [APH(2")] enzyme from the Streptococcus faecalis plasmid pIP800 was cloned in Escherichia coli. A single protein with an apparent molecular weight of 56,000 was specified by this cloned determinant as detected in minicell experiments. Nucleotide sequence analysis revealed the presence of an open reading frame capable of specifying a protein of 479 amino acids and with a molecular weight of 56,850. The deduced amino acid sequence of the bifunctional AAC(6')-APH(2") gene product possessed two regions of homology with other sequenced resistance proteins. The N-terminal region contained a sequence that was homologous to the chloramphenicol acetyltransferase of Bacillus pumilus, and the C-terminal region contained a sequence homologous to the aminoglycoside phosphotransferase of Streptomyces fradiae. Subcloning experiments were performed with the AAC(6')-APH(2") resistance determinant, and it was possible to obtain gene segments independently specifying the acetyltransferase and phosphotransferase activities. These data suggest that the gene specifying the AAC(6')-APH(2") resistance enzyme arose as a result of a gene fusion.     

Cloning and analysis of the four genes coding for Bpu10I restriction-modification enzymes.

The Bpu 10I R-M system from Bacillus pumilus 10, which recognizes the asymmetric 5'-CCTNAGC sequence, has been cloned, sequenced and expressed in Escherichia coli . The system comprises four adjacent, similarly oriented genes encoding two m5C MTases and two subunits of Bpu 10I ENase (34.5 and 34 kDa). Both bpu10IR genes either in cis or trans are needed for the manifestation of R. Bpu 10I activity. Subunits of R. Bpu 10I, purified to apparent homogeneity, are both required for cleavage activity. This heterosubunit structure distinguishes the Bpu 10I restriction endonuclease from all other type II restriction enzymes described previously. The subunits reveal 25% amino acid identity. Significant similarity was also identified between a 43 amino acid region of R. Dde I and one of the regions of higher identity shared between the Bpu 10I subunits, a region that could possibly include the catalytic/Mg2+binding center. The similarity between Bpu 10I and Dde I MTases is not limited to the conserved motifs (CM) typical for m5C MTases. It extends into the variable region that lies between CMs VIII and IX. Duplication of a progenitor gene, encoding an enzyme recognizing a symmetric nucleotide sequence, followed by concerted divergent evolution, may provide a possible scenario leading to the emergence of the Bpu 10I ENase, which recognizes an overall asymmetric sequence and cleaves within it symmetrically.     

High plasticity of multispecific DNA methyltransferases in the region carrying DNA target recognizing enzyme modules.

Multispecific cytosine C5 DNA methyltransferases (MTases) methylate more than one specific DNA target. This is due to the presence of several target recognizing domains (TRDs) in these enzymes. Such TRDs form part of a variable centre in the MTase primary sequence, which separates conserved enzyme core sequences responsible for general steps in the methylation reaction. By deleting, rearranging and exchanging several TRDs of multispecific MTases, we demonstrate their modular character; they mediate target recognition independent of a particular TRD or core sequence context. We show also that multispecific MTases can accommodate inert material of non-MTase origin within their variable region without losing their activity. The remarkable plasticity with respect to the material that can be integrated into this region suggests that the enzyme core sequences preceding or following it form separable functional domains. In spite of the documented flexibility multispecific MTases could not be endowed with novel specificities by integration of putative TRDs of monospecific MTases, pointing to differences between multi- and monospecific MTases in the way their core and TRD sequences interact.     

A nuclear-encoded function essential for translation of the chloroplast psaB mRNA in chlamydomonas.

We report the analysis of a photosystem I-deficient mutant of Chlamydomonas reinhardtii, F15, that contains a mutation at the TAB1 (for translation of psaB mRNA) nuclear locus. Pulse labeling of chloroplast proteins revealed that the synthesis of the two photosystem I reaction center polypeptides PSAA and PSAB was undetectable in this mutant. The mRNA levels of these proteins were only moderately reduced, suggesting that the primary defect occurs at a step during or after translation. We constructed chimeric genes consisting of the psaA and psaB 5' untranslated region (5' UTR) fused to the aminoglycoside adenyltransferase (aadA) coding sequence, which confers spectinomycin resistance. Insertion of these genes into the chloroplast genome through biolistic transformation and analysis of their expression in the TAB1 mutant nuclear background revealed that the psaB (but not the psaA) 5' UTR is the target of the wild-type TAB1 function. This suggests that TAB1 is required for the initiation of psaB mRNA translation. The dependence of PSAA synthesis or accumulation on PSAB synthesis is strongly suggested by the identification of a suppressor mutation within the psaB 5' UTR. The suppressor specifically restores the synthesis of both proteins in the presence of the tab1-F15 mutation. The location of the suppressor mutation within a putative base-paired region near the psaB initiation codon suggests a role for TAB1 in the activation of translation of the psaB mRNA.     

Phylogenomic analysis of 16S rRNA:(guanine-N2) methyltransferases suggests new family members and reveals highly conserved motifs and a domain structure similar to other nucleic acid amino-methyltransferases.

The sequences of known Escherichia coli 16S rRNA:m2G1207 methyltransferase (MTase) RsmC and hypothetical 16S rRNA:m2G966 MTase encoded by the ygjo open reading frame were used to carry out a database search of other putative m2G-generating enzymes in finished and unfinished genomic sequences. Sequence comparison and phylogenetic analysis of 21 close homologs of RsmC and YgjO revealed the presence of the third paralogous lineage in E. coli and other gamma-Proteobacteria, which might correspond to the subfamily of MTases specific for G1516 in 16S rRNA. In addition, the comparative sequence analysis supported by sequence/structure threading suggests that rRNA:m2G MTases are very closely related to RNA and DNA:m6A MTases and that these two enzyme families share common architecture of the active site and presumably a similar mechanism of methyl group transfer onto the exocyclic amino group of their target bases.     

A novel cDNA from Drosophila encoding a protein with similarity to mammalian cysteine-rich secretory proteins,wasp venom antigen 5, and plant group 1 pathogenesis-related proteins

The CAP protein family is made up of a group of secreted proteins that share sequence similarity. Members of this family are found in animals, plants, and fungi, and their shared sequence similarity suggests that members share a common, but as yet unknown, molecular function. As a first step in defining the function of CAP family proteins, an 878 bp partial cDNA encoding a novel member of the CAP family was cloned by the polymerase chain reaction (PCR) from total RNA of adult Drosophila. The cDNA contained the complete coding sequence for a protein 256 amino acids in length, as well as the complete 3′ untranslated region (UTR) and a portion of the 5′ UTR. The protein, named Antigen 5-related (Agr), was most similar in sequence to antigen 5 (Ag5), a CAP family member found in social wasps and ants. The corresponding Agr RNA is about 1 kb in length and is present at all stages of development, with highest levels observed in adults. Agr RNA is transcribed from a single gene that is located within region 12F of the X chromosome. The identification of Agr in Drosophila expands the number of known CAP family members to well over four dozen. Further studies of Agr and the gene which encodes this protein using the Drosophila model system may help provide important insight into the molecular functioning of this little known, but increasingly significant protein family.     

In silico analysis of the tRNA:m1A58 methyltransferase family: homology-based fold prediction and identification of new members from Eubacteria and Archaea.

The amino acid sequences of Gcd10p and Gcd14p, the two subunits of the tRNA:(1-methyladenosine-58; m(1)A58) methyltransferase (MTase) of Saccharomyces cerevisiae, have been analyzed using iterative sequence database searches and fold recognition programs. The results suggest that the 'catalytic' Gcd14p and 'substrate binding' Gcd10p are related to each other and to a group of prokaryotic open reading frames, which were previously annotated as hypothetical protein isoaspartate MTases in sequence databases. It is predicted that the prokaryotic proteins are genuine tRNA:m(1)A MTases based on similarity of their predicted active site to the Gcd14p family. In addition to the MTase domain, an additional domain was identified in the N-terminus of all these proteins that may be involved in interaction with tRNA. These results suggest that the eukaryotic tRNA:m(1)A58 MTase is a product of gene duplication and divergent evolution of a possibly homodimeric prokaryotic enzyme.     

On the definition, nomenclature and classification of water channel proteins (aquaporins and relatives)

A water channel protein (WCP) or a water channel can be defined as a transmembrane protein that has a specific three-dimensional structure with a pore that provides a pathway for water permeation across biological membranes. The pore is formed by two highly conserved regions in the amino acid sequence, called NPA boxes (or motifs) with three amino acid residues (asparagine-proline-alanine, NPA) and several surrounding amino acids. The NPA boxes have been called the "signature" sequence of WCPs. WCPs are a family of proteins belonging to the Membrane Intrinsic Proteins (MIPs) superfamily. In addition, in the MIP superfamily (with more than 1000 members) there are also proteins with no channel activity. The WCP family include three subfamilies: aquaporins, aquaglyceroporins and S-aquaporins. (1) The aquaporins (AQPs) are water selective or specific water channels, also named by various authors as "orthodox", "ordinary", "conventional", "classical", "pure", "normal", or "sensu strictu" aquaporins); (2) The aquaglyceroporins are permeable to water, but also to other small uncharged molecules, in particular glycerol; this family includes the glycerol facilitators, abbreviated as GlpFs, from glycerol permease facilitators. The "signature" sequence for aquaglyceroporins is the aspartic acid residue (D) in the second NPA box. (3) The third subfamily of WCPs have little conserved amino acid sequences around the NPA boxes, unclassifiable to the first two subfamilies. I recommend to use always for this subfamily the name S-aquaporins. They are also named "superaquaporins", "aquaporins with unusual (or deviated) NPA boxes", "subcellular aquaporins", or "sip-like aquaporins". I also recommend to use always the spelling aquaporin (not aquaporine), and, for various AQPs, the abbreviation AQP followed immediately by the number, (e.g. AQP1), with no space or - which might create confusions with "minus".     

Molecular phylogenetics of the RrmJ/fibrillarin superfamily of ribose 2'-O-methyltransferases

Recent analyses identified a putative catalytic tetrad K-D-K-E common to several families of site-specific methyltransferases (MTases) that modify 2'-hydroxyl groups of ribose in mRNA, rRNA and tRNA (designated the RrmJ class after one of the structurally characterized members; 1eiz in Protein Data Bank) [Genome Biol. 2(9) (2001) 38]. Subsequently, three residues of the tetrad (K-D-K) were shown to be essential for catalysis in RrmJ [J. Biol. Chem. 277 (2002) 41978]. Here, we report identification of a similar conserved tetrad (K-D-K-H) in the family of snoRNA-guided ribose 2'-O-MTases related to fibrillarin (represented by the Mj0697 protein structure; 1fbn in PDB). The corresponding functional groups of putative catalytic tetrads of RrmJ and Mj0697 may be superimposed in space. However, one of the invariant residues (K(164) in RrmJ and K(179) in Mj0697) is observed in two distinct locations in the primary sequence, suggesting an interesting case of 'migration' of the conserved side chain within the framework of the active site. RrmJ and Mj0697 sequences were used as starting points to carry out comprehensive sequence database searches, resulting in identification of a similar conserved tetrad (and hence, prediction of a ribose 2'-O-specificity) in several families of putative MTases, including TlyA hemolysins, novel proteins from Trypanosoma, and large multidomain proteins from Flaviviriruses, Nidoviruses, and Alphaviruses. The results of our analysis of phylogenetic relationships in the RrmJ/fibrillarin superfamily provide insight into the evolution of site-specific and snoRNA-guided ribose 2'-O-MTases from a common ancestor.     

The amino acid sequence of the eukaryotic DNA [N6-adenine]methyltransferase, M.CviBIII, has regions of similarity with the prokaryotic isoschizomer M.TaqI and other DNA [N6-adenine] methyltransferases

The sequences of the genes coding for M.CviBIII (from virus NC-1A which infects a eukaryotic alga) [Narva et al., Nucleic Acids Res. 15 (1987) 9807-9823] and M.TaqI (from the bacterium Thermus aquaticus) [Slatko et al., Nucleic Acids Res. 15 (1987) 9781-9796] have been determined recently. Both enzymes methylate adenine in the sequence TCGA. We have compared the predicted amino acid sequences of these two methyltransferases (MTases), with each other and with ten other N6 A-MTases and find regions of similarity. M.CviBIII and M.TaqI were most closely related followed by M.PaeR7, whose recognition sequence (CTCGAG) contains the M.TaqI/M.CviBIII recognition sequence TCGA, and M.PstI, whose recognition sequence is CTGCAG. All of the N6-MTases contain the sequence Asp/Asn-Pro-Pro-Tyr (B-P-P-Y) referred to by Hattman et al. [J. Bacteriol. 164 (1985) 932-937] as region IV. The predicted secondary structure of this region forms a finger-like structure ('beta finger') containing a beta-pleated sheet (...XXXB), two beta-turns (P-P) followed by another beta-pleated sheet [Y/FXXX...].     

Bacteriophage Orphan DNA Methyltransferases: Insights from Their Bacterial Origin,Function, and Occurrence

Type II DNA methyltransferases (MTases) are enzymes found ubiquitously in the prokaryotic world, where they play important roles in several cellular processes, such as host protection and epigenetic regulation. Three classes of type II MTases have been identified thus far in bacteria which function in transferring a methyl group from S-adenosyl-l-methionine (SAM) to a target nucleotide base, forming N-6-methyladenine (class I), N-4-methylcytosine (class II), or C-5-methylcytosine (class III). Often, these MTases are associated with a cognate restriction endonuclease (REase) to form a restriction-modification (R-M) system protecting bacterial cells from invasion by foreign DNA. When MTases exist alone, which are then termed orphan MTases, they are believed to be mainly involved in regulatory activities in the bacterial cell. Genomes of various lytic and lysogenic phages have been shown to encode multi- and mono-specific orphan MTases that have the ability to confer protection from restriction endonucleases of their bacterial host(s). The ability of a phage to overcome R-M and other phage-targeting resistance systems can be detrimental to particular biotechnological processes such as dairy fermentations. Conversely, as phages may also be beneficial in certain areas such as phage therapy, phages with additional resistance to host defenses may prolong the effectiveness of the therapy. This minireview will focus on bacteriophage-encoded MTases, their prevalence and diversity, as well as their potential origin and function.