Functional specialization of domains tandemly duplicated within 16S rRNA methyltransferase RsmC

RNA methyltransferases (MTases) are important players in the biogenesis and regulation of the ribosome, the cellular machine for protein synthesis. RsmC is a MTase that catalyzes the transfer of a methyl group from S-adenosyl-l-methionine (SAM) to G1207 of 16S rRNA. Mutations of G1207 have dominant lethal phenotypes in Escherichia coli, underscoring the significance of this modified nucleotide for ribosome function. Here we report the crystal structure of E. coli RsmC refined to 2.1 A resolution, which reveals two homologous domains tandemly duplicated within a single polypeptide. We characterized the function of the individual domains and identified key residues involved in binding of rRNA and SAM, and in catalysis. We also discovered that one of the domains is important for the folding of the other. Domain duplication and subfunctionalization by complementary degeneration of redundant functions (in particular substrate binding versus catalysis) has been reported for many enzymes, including those involved in RNA metabolism. Thus, RsmC can be regarded as a model system for functional streamlining of domains accompanied by the development of dependencies concerning folding and stability.     

RNA:(guanine-N2) methyltransferases RsmC/RsmD and their homologs revisited – bioinformatic analysis and prediction of the active site based on the uncharacterized Mj0882 protein structure

Background  

Escherichia coli guanine-N2 (m²G) methyltransferases (MTases) RsmC and RsmD modify nucleosides G1207 and G966 of 16S rRNA. They possess a common MTase domain in the C-terminus and a variable region in the N-terminus. Their C-terminal domain is related to the YbiN family of hypothetical MTases, but nothing is known about the structure or function of the N-terminal domain.     

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.     

Sequence-structure-function relationships of Tgs1, the yeast snRNA/snoRNA cap hypermethylase

The Saccharomyces cerevisiae Tgs1 methyltransferase (MTase) is responsible for conversion of the m⁷G caps of snRNAs and snoRNAs to a 2,2,7- trimethylguanosine structure. To learn more about the evolutionary origin of Tgs1 and to identify structural features required for its activity, we performed a structure–function study. By using sequence comparison and phylogenetic analysis, we found that Tgs1 shows strongest similarity to Mj0882, a protein related to a family comprised of bacterial rRNA:m²G MTases RsmC and RsmD. The structural information of Mj0882 was used to build a homology model of Tgs1p which allowed us to predict the range of the minimal globular MTase domain and the localization of other residues that may be important for enzyme function. To further characterize functional domains of Tgs1, mutants were constructed and tested for their effects on cell viability, subcellular localization and binding to the small nuclear ribonucleoproteins (snRNPs) and small nucleolar RNPs (snoRNPs). We found that the N-terminal domain of the hypermethylase is dispensable for binding to the common snRNPs and snoRNPs proteins but essential for correct nucleolar localization. Site- directed mutagenesis of Tgs1 allowed also the identification of the residues likely to be involved in the formation of the m7G-binding site and the catalytic center.     

Sequence analysis and structure prediction of aminoglycoside-resistance 16S rRNA:m7G methyltransferases

Methylation of G1405 within bacterial 16S ribosomal RNA results in high-level resistance to specific combinations of aminoglycoside antibiotics. Only a few closely related methyltransferases (MTases), which carry out the respective modification (here dubbed "Agr", for aminoglycoside resistance), are known. It is not clear, whether they are related to "typical" S-adenosylmethionine (AdoMet)-dependent MTases or not. Demydchuk et al., 1998 proposed that the cofactor-binding region is localized at the C-terminus of Agr MTases, which implies an interesting case of sequence permutation. Since the Agr MTases lack significant sequence similarity to other proteins, we tested that hypothesis using more sensitive sequence/structure threading approach. Structure prediction confirmed the presence of a putative AdoMet-binding site in these proteins, albeit at a distinct location, resembling that of "typical", non-permuted MTases. Additionally, a small alpha-helical domain dissimilar to other proteins in the database was identified in the N-terminal region of Agr MTases. Comparison of a three-dimensional model of the Agr family member with a recently solved structure of reovirus mRNA capping MTase suggests that the mechanism of guanine-N7 methylation in rRNA and mRNA may be different.     

Methyltransferase that modifies guanine 966 of the 16 S rRNA: functional identification and tertiary structure

N(2)-Methylguanine 966 is located in the loop of Escherichia coli 16 S rRNA helix 31, forming a part of the P-site tRNA-binding pocket. We found yhhF to be a gene encoding for m(2)G966 specific 16 S rRNA methyltransferase. Disruption of the yhhF gene by kanamycin resistance marker leads to a loss of modification at G966. The modification could be rescued by expression of recombinant protein from the plasmid carrying the yhhF gene. Moreover, purified m(2)G966 methyltransferase, in the presence of S-adenosylomethionine (AdoMet), is able to methylate 30 S ribosomal subunits that were purified from yhhF knock-out strain in vitro. The methylation is specific for G966 base of the 16 S rRNA. The m(2)G966 methyltransferase was crystallized, and its structure has been determined and refined to 2.05A(.) The structure closely resembles RsmC rRNA methyltransferase, specific for m(2)G1207 of the 16 S rRNA. Structural comparisons and analysis of the enzyme active site suggest modes for binding AdoMet and rRNA to m(2)G966 methyltransferase. Based on the experimental data and current nomenclature the protein expressed from the yhhF gene was renamed to RsmD. A model for interaction of RsmD with ribosome has been proposed.     

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.     

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.     

In silico identification, structure prediction and phylogenetic analysis of the 2'-O-ribose (cap 1) methyltransferase domain in the large structural protein of ssRNA negative-strand viruses

The Escherichia coli RrmJ gene product has recently been shown to be the 23S rRNA:U2552 specific 2'-O-ribose methyltransferase (MTase) (RrmJ). Its structure has been solved and refined to 1.5 A resolution, demonstrating conservation of the three-dimensional fold and key catalytic side chains with the vaccinia virus VP39 protein, which functions as an mRNA 5'm(7)G-cap-N-specific 2'-O-ribose MTase. Using the amino acid sequence of RrmJ as an initial probe in an iterative search of sequence databases, we identified a homologous domain in the sequence of the L protein of non-segmented, negative-sense, single-stranded RNA viruses. The plausibility of the prediction was confirmed by homology modeling and checking whether important residues at substrate/ligand-binding sites were conserved. The predicted structural compatibility and the conservation of the active site between the novel putative MTase domain and genuine 2'-O-ribose MTases, together with the available results of biochemical studies, strongly suggest that this domain is a 5'm(7)G-cap-N-specific 2'-O-ribose MTase (i.e. the cap 1 MTase). Evolutionary relationships between these proteins are also discussed.     

Structural insights into the function of 23S rRNA methyltransferase RlmG (m²G1835) from Escherichia coli

RlmG is a specific AdoMet-dependent methyltransferase (MTase) responsible for N2-methylation of G1835 in 23S rRNA of Escherichia coli. Methylation of m2G1835 specifically enhances association of ribosomal subunits and provides a significant advantage for bacteria in osmotic and oxidative stress. Here, the crystal structure of RlmG in complex with AdoMet and its structure in solution were determined. The structure of RlmG is similar to that of the MTase RsmC, consisting of two homologous domains: the N-terminal domain (NTD) in the recognition and binding of the substrate, and the C-terminal domain (CTD) in AdoMet-binding and the catalytic process. However, there are distinct positively charged protuberances and a distribution of conserved residues contributing to the charged surface patch, especially in the NTD of RlmG for direct binding of protein-free rRNA. The RNA-binding properties of the NTD and CTD characterized by both gel electrophoresis mobility shift assays and isothermal titration calorimetry showed that NTD could bind RNA independently and RNA binding was achieved by the NTD, accomplished by a coordinating role of the CTD. The model of the RlmG-AdoMet-RNA complex suggested that RlmG may unfold its substrate RNA in the positively charged cleft between the NTD and CTD, and then G1835 disengages from its Watson-Crick pairing with C1905 and flips out to insert into the active site. Our structure and biochemical studies provide novel insights into the catalytic mechanism of G1835 methylation.     

Sequence-structure-function studies of tRNA:m5C methyltransferase Trm4p and its relationship to DNA:m5C and RNA:m5U methyltransferases

Three types of methyltransferases (MTases) generate 5-methylpyrimidine in nucleic acids, forming m5U in RNA, m5C in RNA and m5C in DNA. The DNA:m5C MTases have been extensively studied by crystallographic, biophysical, biochemical and computational methods. On the other hand, the sequence-structure-function relationships of RNA:m5C MTases remain obscure, as do the potential evolutionary relationships between the three types of 5-methylpyrimidine-generating enzymes. Sequence analyses and homology modeling of the yeast tRNA:m5C MTase Trm4p (also called Ncl1p) provided a structural and evolutionary platform for identification of catalytic residues and modeling of the architecture of the RNA:m5C MTase active site. The analysis led to the identification of two invariant residues that are important for Trm4p activity in addition to the conserved Cys residues in motif IV and motif VI that were previously found to be critical. The newly identified residues include a Lys residue in motif I and an Asp in motif IV. A conserved Gln found in motif X was found to be dispensable for MTase activity. Locations of essential residues in the model of Trm4p are in very good agreement with the X-ray structure of an RNA:m5C MTase homolog PH1374. Theoretical and experimental analyses revealed that RNA:m5C MTases share a number of features with either RNA:m5U MTases or DNA:m5C MTases, which suggested a tentative phylogenetic model of relationships between these three classes of 5-methylpyrimidine MTases. We infer that RNA:m5C MTases evolved from RNA:m5U MTases by acquiring an additional Cys residue in motif IV, which was adapted to function as the nucleophilic catalyst only later in DNA:m5C MTases, accompanied by loss of the original Cys from motif VI, transfer of a conserved carboxylate from motif IV to motif VI and sequence permutation.     

Sequence-structure-function relationships of a tRNA (m7G46) methyltransferase studied by homology modeling and site-directed mutagenesis

The Escherichia coli TrmB protein and its Saccharomyces cerevisiae ortholog Trm8p catalyze the S-adenosyl-L-methionine-dependent formation of 7-methylguanosine at position 46 (m7G46) in tRNA. To learn more about the sequence-structure-function relationships of these enzymes we carried out a thorough bioinformatics analysis of the tRNA:m7G methyltransferase (MTase) family to predict sequence regions and individual amino acid residues that may be important for the interactions between the MTase and the tRNA substrate, in particular the target guanosine 46. We used site-directed mutagenesis to construct a series of alanine substitutions and tested the activity of the mutants to elucidate the catalytic and tRNA-recognition mechanism of TrmB. The functional analysis of the mutants, together with the homology model of the TrmB structure and the results of the phylogenetic analysis, revealed the crucial residues for the formation of the substrate-binding site and the catalytic center in tRNA:m7G MTases.     

Identification and preliminary characterization of an O6-methylguanine DNA repair methyltransferase in the yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae contains a DNA repair methyltransferase (MTase) that repairs O6-methylguanine. Methyl groups are irreversibly transferred from O6-methylguanine in DNA to a 25-kilodalton protein in S. cerevisiae cell extracts, and methyl transfer is accompanied by the formation of S-methylcysteine. The yeast MTase is expressed at approximately 150 molecules/cell in exponentially growing yeast cultures but is not detectable in stationary phase cells. Unlike mammalian and bacterial MTases, the yeast MTase is very temperature-sensitive, having a half-life of about 4 min at 37 degrees C, which may explain why others have failed to detect it. Like other DNA repair MTases, the S. cerevisiae MTase repairs O6-methylguanine more efficiently in double-stranded DNA than in single-stranded DNA. Synthesis of the yeast DNA MTase is apparently not inducible by sublethal exposures to alkylating agent, but rather MTase activity is depleted in cells exposed to low doses of alkylating agent. Judging from its molecular weight and substrate specificity, the yeast DNA MTase is more closely related to mammalian MTases than to Escherichia coli MTases.     

Cloning and analysis of the HaeIII and HaeII methyltransferase genes

The HaeIII methyltransferase (MTase) gene from Haemophilus aegyptius (recognition sequence: 5'-GGCC-3') was cloned into Escherichia coli in the plasmid vector pBR322. The gene was isolated on a single EcoRI fragment and on a single HindIII fragment. Clones carrying additional adjacent fragments were found to code also for the HaeII restriction endonuclease and HaeII modification MTase (recognition sequence: 5'-PuGCGCPy-3'). The sequence of the HaeIII modification gene was determined. The inferred amino acid sequence of the protein was found to share extensive similarity with other sequenced m5C-MTases. The central 'non-conserved' region of the M.HaeIII MTase, thought to form the nucleotide sequence-specificity domain, is almost identical to that of the M.BsuRI, M.BspRI and M.NgoPII MTases, which also recognize the sequence 5'-GGCC-3'.     

Refolding of a fully functional flavivirus methyltransferase revealed that S‐adenosyl methionine but not S‐adenosyl homocysteine is copurified with flavivirus methyltransferase

Methylation of flavivirus RNA is vital for its stability and translation in the infected host cell. This methylation is mediated by the flavivirus methyltransferase (MTase), which methylates the N7 and 2′‐O positions of the viral RNA cap by using S‐adenosyl‐l ‐methionine (SAM) as a methyl donor. In this report, we demonstrate that SAM, in contrast to the reaction by‐product S‐adenosyl‐l ‐homocysteine, which was assumed previously, is copurified with the Dengue (DNV) and West Nile virus MTases produced in Escherichia coli (E. coli). This endogenous SAM can be removed by denaturation and refolding of the MTase protein. The refolded MTase of DNV serotype 3 (DNV3) displays methylation activity comparable to native enzyme, and its crystal structure at 2.1 Å is almost identical to that of native MTase. We characterized the binding of Sinefungin (SIN), a previously described SAM‐analog inhibitor of MTase function, to the native and refolded DNV3 MTase by isothermal titration calorimetry, and found that SIN binds to refolded MTase with more than 16 times the affinity of SIN binding to the MTase purified natively. Moreover, we show that SAM is also copurified with other flavivirus MTases, indicating that purification by refolding may be a generally applicable tool for studying flavivirus MTase inhibition.     

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

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

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

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

Ribosomal RNA guanine-(N2)-methyltransferases and their targets

Five nearly universal methylated guanine-(N2) residues are present in bacterial rRNA in the ribosome. To date four out of five ribosomal RNA guanine-(N2)-methyltransferases are described. RsmC(YjjT) methylates G1207 of the 16S rRNA. RlmG(YgjO) and RlmL(YcbY) are responsible for the 23S rRNA m²G1835 and m²G2445 formation, correspondingly. RsmD(YhhF) is necessary for methylation of G966 residue of 16S rRNA. Structure of Escherichia coli RsmD(YhhF) methyltransferase and the structure of the Methanococcus jannaschii RsmC ortholog were determined. All ribosomal guanine-(N2)-methyltransferases have similar AdoMet-binding sites. In relation to the ribosomal substrate recognition, two enzymes that recognize assembled subunits are relatively small single domain proteins and two enzymes that recognize naked rRNA are larger proteins containing separate methyltransferase- and RNA-binding domains. The model for recognition of specific target nucleotide is proposed. The hypothetical role of the m²G residues in rRNA is discussed.