Structural and evolutionary bioinformatics of the SPOUT superfamily of methyltransferases

Structural insight into the functional mechanism of Nep1/Emg1 N1-specific pseudouridine methyltransferase in ribosome biogenesis

Nucleolar Essential Protein 1 (Nep1) is required for small subunit (SSU) ribosomal RNA (rRNA) maturation and is mutated in Bowen–Conradi Syndrome. Although yeast (Saccharomyces cerevisiae) Nep1 interacts with a consensus sequence found in three regions of SSU rRNA, the molecular details of the interaction are unknown. Nep1 is a SPOUT RNA methyltransferase, and can catalyze methylation at the N1 of pseudouridine. Nep1 is also involved in assembly of Rps19, an SSU ribosomal protein. Mutations in Nep1 that result in decreased methyl donor binding do not result in lethality, suggesting that enzymatic activity may not be required for function, and RNA binding may play a more important role. To study these interactions, the crystal structures of the scNep1 dimer and its complexes with RNA were determined. The results demonstrate that Nep1 recognizes its RNA site via base-specific interactions and stabilizes a stem-loop in the bound RNA. Furthermore, the RNA structure observed contradicts the predicted structures of the Nep1-binding sites within mature rRNA, suggesting that the Nep1 changes rRNA structure upon binding. Finally, a uridine base is bound in the active site of Nep1, positioned for a methyltransfer at the C5 position, supporting its role as an N1-specific pseudouridine methyltransferase.

Structural basis for the methylation of A1408 in 16S rRNA by a panaminoglycoside resistance methyltransferase NpmA from a clinical isolate and analysis of the NpmA interactions with the 30S ribosomal subunit

NpmA, a methyltransferase that confers resistance to aminoglycosides was identified in an Escherichia coli clinical isolate. It belongs to the kanamycin–apramycin methyltransferase (Kam) family and specifically methylates the 16S rRNA at the N1 position of A1408. We determined the structures of apo-NpmA and its complexes with S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy) at 2.4, 2.7 and 1.68 Å, respectively. We generated a number of NpmA variants with alanine substitutions and studied their ability to bind the cofactor, to methylate A1408 in the 30S subunit, and to confer resistance to kanamycin in vivo. Residues D30, W107 and W197 were found to be essential. We have also analyzed the interactions between NpmA and the 30S subunit by footprinting experiments and computational docking. Helices 24, 42 and 44 were found to be the main NpmA-binding site. Both experimental and theoretical analyses suggest that NpmA flips out the target nucleotide A1408 to carry out the methylation. NpmA is plasmid-encoded and can be transferred between pathogenic bacteria; therefore it poses a threat to the successful use of aminoglycosides in clinical practice. The results presented here will assist in the development of specific NpmA inhibitors that could restore the potential of aminoglycoside antibiotics.

YibK is the 2'-O-methyltransferase TrmL that modifies the wobble nucleotide in Escherichia coli tRNA(Leu) isoacceptors

Transfer RNAs are the most densely modified nucleic acid molecules in living cells. In Escherichia coli, more than 30 nucleoside modifications have been characterized, ranging from methylations and pseudouridylations to more complex additions that require multiple enzymatic steps. Most of the modifying enzymes have been identified, although a few notable exceptions include the 2'-O-methyltransferase(s) that methylate the ribose at the nucleotide 34 wobble position in the two leucyl isoacceptors tRNA(Leu)(CmAA) and tRNA(Leu)(cmnm5UmAA). Here, we have used a comparative genomics approach to uncover candidate E. coli genes for the missing enzyme(s). Transfer RNAs from null mutants for candidate genes were analyzed by mass spectrometry and revealed that inactivation of yibK leads to loss of 2'-O-methylation at position 34 in both tRNA(Leu)(CmAA) and tRNA(Leu)(cmnm5UmAA). Loss of YibK methylation reduces the efficiency of codon-wobble base interaction, as demonstrated in an amber suppressor supP system. Inactivation of yibK had no detectable effect on steady-state growth rate, although a distinct disadvantage was noted in multiple-round, mixed-population growth experiments, suggesting that the ability to recover from the stationary phase was impaired. Methylation is restored in vivo by complementing with a recombinant copy of yibK. Despite being one of the smallest characterized α/β knot proteins, YibK independently catalyzes the methyl transfer from S-adenosyl-L-methionine to the 2'-OH of the wobble nucleotide; YibK recognition of this target requires a pyridine at position 34 and N⁶-(isopentenyl)-2-methylthioadenosine at position 37. YibK is one of the last remaining E. coli tRNA modification enzymes to be identified and is now renamed TrmL.

The Bowen-Conradi syndrome protein Nep1 (Emg1) has a dual role in eukaryotic ribosome biogenesis, as an essential assembly factor and in the methylation of Ψ1191 in yeast 18S rRNA

The Nep1 (Emg1) SPOUT-class methyltransferase is an essential ribosome assembly factor and the human Bowen–Conradi syndrome (BCS) is caused by a specific Nep1^D86G mutation. We recently showed in vitro that Methanocaldococcus jannaschii Nep1 is a sequence-specific pseudouridine-N1-methyltransferase. Here, we show that in yeast the in vivo target site for Nep1-catalyzed methylation is located within loop 35 of the 18S rRNA that contains the unique hypermodification of U1191 to 1-methyl-3-(3-amino-3-carboxypropyl)-pseudouri-dine (m1acp3Ψ). Specific ¹⁴C-methionine labelling of 18S rRNA in yeast mutants showed that Nep1 is not required for acp-modification but suggested a function in Ψ1191 methylation. ESI MS analysis of acp-modified Ψ-nucleosides in a Δnep1-mutant showed that Nep1 catalyzes the Ψ1191 methylation in vivo. Remarkably, the restored growth of a nep1-1^ts mutant upon addition of S-adenosylmethionine was even observed after preventing U1191 methylation in a Δsnr35 mutant. This strongly suggests a dual Nep1 function, as Ψ1191-methyltransferase and ribosome assembly factor. Interestingly, the Nep1 methyltransferase activity is not affected upon introduction of the BCS mutation. Instead, the mutated protein shows enhanced dimerization propensity and increased affinity for its RNA-target in vitro. Furthermore, the BCS mutation prevents nucleolar accumulation of Nep1, which could be the reason for reduced growth in yeast and the Bowen-Conradi syndrome.

The protein composition of mitotic chromosomes determined using multiclassifier combinatorial proteomics

Despite many decades of study, mitotic chromosome structure and composition remain poorly characterized. Here, we have integrated quantitative proteomics with bioinformatic analysis to generate a series of independent classifiers that describe the ∼4,000 proteins identified in isolated mitotic chromosomes. Integrating these classifiers by machine learning uncovers functional relationships between protein complexes in the context of intact chromosomes and reveals which of the ∼560 uncharacterized proteins identified here merits further study. Indeed, of 34 GFP-tagged predicted chromosomal proteins, 30 were chromosomal, including 13 with centromere-association. Of 16 GFP-tagged predicted nonchromosomal proteins, 14 were confirmed to be nonchromosomal. An unbiased analysis of the whole chromosome proteome from genetic knockouts of kinetochore protein Ska3/Rama1 revealed that the APC/C and RanBP2/RanGAP1 complexes depend on the Ska complex for stable association with chromosomes. Our integrated analysis predicts that up to 97 new centromere-associated proteins remain to be discovered in our data set.

New archaeal methyltransferases forming 1-methyladenosine or 1-methyladenosine and 1-methylguanosine at position 9 of tRNA

Two archaeal tRNA methyltransferases belonging to the SPOUT superfamily and displaying unexpected activities are identified. These enzymes are orthologous to the yeast Trm10p methyltransferase, which catalyses the formation of 1-methylguanosine at position 9 of tRNA. In contrast, the Trm10p orthologue from the crenarchaeon Sulfolobus acidocaldarius forms 1-methyladenosine at the same position. Even more surprisingly, the Trm10p orthologue from the euryarchaeon Thermococcus kodakaraensis methylates the N¹-atom of either adenosine or guanosine at position 9 in different tRNAs. This is to our knowledge the first example of a tRNA methyltransferase with a broadened nucleoside recognition capability. The evolution of tRNA methyltransferases methylating the N¹ atom of a purine residue is discussed.

Megadalton complexes in the chloroplast stroma of Arabidopsis thaliana characterized by size exclusion chromatography, mass spectrometry, and hierarchical clustering

To characterize MDa-sized macromolecular chloroplast stroma protein assemblies and to extend coverage of the chloroplast stroma proteome, we fractionated soluble chloroplast stroma in the non-denatured state by size exclusion chromatography with a size separation range up to approximately 5 MDa. To maximize protein complex stability and resolution of megadalton complexes, ionic strength and composition were optimized. Subsequent high accuracy tandem mass spectrometry analysis (LTQ-Orbitrap) identified 1081 proteins across the complete native mass range. Protein complexes and assembly states above 0.8 MDa were resolved using hierarchical clustering, and protein heat maps were generated from normalized protein spectral counts for each of the size exclusion chromatography fractions; this complemented previous analysis of stromal complexes up to 0.8 MDa (Peltier, J. B., Cai, Y., Sun, Q., Zabrouskov, V., Giacomelli, L., Rudella, A., Ytterberg, A. J., Rutschow, H., and van Wijk, K. J. (2006) The oligomeric stromal proteome of Arabidopsis thaliana chloroplasts. Mol. Cell. Proteomics 5, 114-133). This combined experimental and bioinformatics analyses resolved chloroplast ribosomes in different assembly and functional states (e.g. 30, 50, and 70 S), which enabled the identification of plastid homologues of prokaryotic ribosome assembly factors as well as proteins involved in co-translational modifications, targeting, and folding. The roles of these ribosome-associating proteins will be discussed. Known RNA splice factors (e.g. CAF1/WTF1/RNC1) as well as uncharacterized proteins with RNA-binding domains (pentatricopeptide repeat, RNA recognition motif, and chloroplast ribosome maturation), RNases, and DEAD box helicases were found in various sized complexes. Chloroplast DNA (>3 MDa) was found in association with the complete heteromeric plastid-encoded DNA polymerase complex, and a dozen other DNA-binding proteins, e.g. DNA gyrase, topoisomerase, and various DNA repair enzymes. The heteromeric >or=5-MDa pyruvate dehydrogenase complex and the 0.8-1-MDa acetyl-CoA carboxylase complex associated with uncharacterized biotin carboxyl carrier domain proteins constitute the entry point to fatty acid metabolism in leaves; we suggest that their large size relates to the need for metabolic channeling. Protein annotations and identification data are available through the Plant Proteomics Database, and mass spectrometry data are available through Proteomics Identifications database.

Experimental detection of knotted conformations in denatured proteins

Structures that contain a knot formed by the path of the polypeptide backbone represent some of the most complex topologies observed in proteins. How or why these topological knots arise remains unclear. By developing a method to experimentally trap and detect knots in nonnative polypeptide chains, we find that two knotted methyltransferases, YibK and YbeA, can exist in a trefoil-knot conformation even in their chemically unfolded states. The unique denatured-state topology of these molecules explains their ability to efficiently fold to their native knotted structures in vitro and offers insights into the potential role of knots in proteins. Furthermore, the high prevalence of the denatured-state knots identified here suggests that they are either difficult to untie or that threading of any untied molecules is rapid and spontaneous. The occurrence of such knotted topologies in unfolded polypeptide chains raises the possibility that they could play an important, and as yet unexplored, role in folding and misfolding processes in vivo.

Structural basis for the methylation of G1405 in 16S rRNA by aminoglycoside resistance methyltransferase Sgm from an antibiotic producer: a diversity of active sites in m7G methyltransferases

Sgm (Sisomicin-gentamicin methyltransferase) from antibiotic-producing bacterium Micromonospora zionensis is an enzyme that confers resistance to aminoglycosides like gentamicin and sisomicin by specifically methylating G1405 in bacterial 16S rRNA. Sgm belongs to the aminoglycoside resistance methyltransferase (Arm) family of enzymes that have been recently found to spread by horizontal gene transfer among disease-causing bacteria. Structural characterization of Arm enzymes is the key to understand their mechanism of action and to develop inhibitors that would block their activity. Here we report the structure of Sgm in complex with cofactors S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy) at 2.0 and 2.1 Å resolution, respectively, and results of mutagenesis and rRNA footprinting, and protein-substrate docking. We propose the mechanism of methylation of G1405 by Sgm and compare it with other m⁷G methyltransferases, revealing a surprising diversity of active sites and binding modes for the same basic reaction of RNA modification. This analysis can serve as a stepping stone towards developing drugs that would specifically block the activity of Arm methyltransferases and thereby re-sensitize pathogenic bacteria to aminoglycoside antibiotics.

Flexible recognition of the tRNA G18 methylation target site by TrmH methyltransferase through first binding and induced fit processes

Transfer RNA (Gm18) methyltransferase (TrmH) catalyzes methyl transfer from S-adenosyl-l-methionine to a conserved G18 in tRNA. We investigated the recognition mechanism of Thermus thermophilus TrmH for its guanosine target. Thirteen yeast tRNA(Phe) mutant transcripts were prepared in which the modification site and/or other nucleotides in the D-loop were substituted by dG, inosine, or other nucleotides. We then conducted methyl transfer kinetic studies, gel shift assays, and inhibition experiments using these tRNA variants. Sites of methylation were confirmed with RNA sequencing or primer extension. Although the G18G19 sequence is not essential for methylation by TrmH, disruption of G18G19 severely reduces the efficiency of methyl transfer. There is strict recognition of guanosine by TrmH, in that methylation occurs at the adjacent G19 when the G18 is replaced by dG or adenosine. The fact that TrmH methylates guanosine in D-loops from 4 to 12 nucleotides in length suggests that selection of the position of guanosine within the D-loop is relatively flexible. Our studies also demonstrate that the oxygen 6 atom of the guanine base is a positive determinant for TrmH recognition. The recognition process of TrmH for substrate is inducible and product-inhibited, in that tRNAs containing Gm18 are excluded by TrmH. In contrast, substitution of G18 with dG18 results in the formation of a more stable TrmH-tRNA complex. To address the mechanism, we performed the stopped-flow pre-steady state kinetic analysis. The result clearly showed that the binding of TrmH to tRNA is composed of at least three steps, the first bi-molecular binding and the subsequent two uni-molecular induced-fit processes.

The ribosome assembly factor Nep1 responsible for Bowen-Conradi syndrome is a pseudouridine-N1-specific methyltransferase

Nep1 (Emg1) is a highly conserved nucleolar protein with an essential function in ribosome biogenesis. A mutation in the human Nep1 homolog causes Bowen-Conradi syndrome-a severe developmental disorder. Structures of Nep1 revealed a dimer with a fold similar to the SPOUT-class of RNA-methyltransferases suggesting that Nep1 acts as a methyltransferase in ribosome biogenesis. The target for this putative methyltransferase activity has not been identified yet. We characterized the RNA-binding specificity of Methanocaldococcus jannaschii Nep1 by fluorescence- and NMR-spectroscopy as well as by yeast three-hybrid screening. Nep1 binds with high affinity to short RNA oligonucleotides corresponding to nt 910-921 of M. jannaschii 16S rRNA through a highly conserved basic surface cleft along the dimer interface. Nep1 only methylates RNAs containing a pseudouridine at a position corresponding to a previously identified hypermodified N1-methyl-N3-(3-amino-3-carboxypropyl) pseudouridine (m1acp3-Psi) in eukaryotic 18S rRNAs. Analysis of the methylated nucleoside by MALDI-mass spectrometry, HPLC and NMR shows that the methyl group is transferred to the N1 of the pseudouridine. Thus, Nep1 is the first identified example of an N1-specific pseudouridine methyltransferase. This enzymatic activity is also conserved in human Nep1 suggesting that Nep1 is the methyltransferase in the biosynthesis of m1acp3-Psi in eukaryotic 18S rRNAs.

Insights into the structure, function and evolution of the radical-SAM 23S rRNA methyltransferase Cfr that confers antibiotic resistance in bacteria

The Cfr methyltransferase confers combined resistance to five classes of antibiotics that bind to the peptidyl tranferase center of bacterial ribosomes by catalyzing methylation of the C-8 position of 23S rRNA nucleotide A2503. The same nucleotide is targeted by the housekeeping methyltransferase RlmN that methylates the C-2 position. Database searches with the Cfr sequence have revealed a large group of closely related sequences from all domains of life that contain the conserved CX(3)CX(2)C motif characteristic of radical S-adenosyl-l-methionine (SAM) enzymes. Phylogenetic analysis of the Cfr/RlmN family suggests that the RlmN subfamily is likely the ancestral form, whereas the Cfr subfamily arose via duplication and horizontal gene transfer. A structural model of Cfr has been calculated and used as a guide for alanine mutagenesis studies that corroborate the model-based predictions of a 4Fe-4S cluster, a SAM molecule coordinated to the iron-sulfur cluster (SAM1) and a SAM molecule that is the putative methyl group donor (SAM2). All mutations at predicted functional sites affect Cfr activity significantly as assayed by antibiotic susceptibility testing and primer extension analysis. The investigation has identified essential amino acids and Cfr variants with altered reaction mechanisms and represents a first step towards understanding the structural basis of Cfr activity.

Backbone resonance assignments of the 48 kDa dimeric putative 18S rRNA-methyltransferase Nep1 from Methanocaldococcus jannaschii

Nep1 from Methanocaldococcus jannaschii is a 48 kDa dimeric protein belonging to the SPOUT-class of S-adenosylmethionine dependent RNA-methyltransferases and acting as a ribosome assembly factor. Mutations in the human homolog are the cause of Bowen-Conradi syndrome. We report here 1H, 15N and 13C chemical shift assignments for the backbone of the protein in its apo state.

Bioinformatic Identification of Novel Methyltransferases

Methylation of DNA, protein, and even RNA species are integral processes in epigenesis. Enzymes that catalyze these reactions using the donor S-adenosylmethionine fall into several structurally distinct classes. The members in each class share sequence similarity that can be used to identify additional methyltransferases. Here, we characterize these classes and in silico approaches to infer protein function. Computational methods such as hidden Markov model profiling and the Multiple Motif Scanning program can be used to analyze known methyltransferases and relay information into the prediction of new ones. In some cases, the substrate of methylation can be inferred from hidden Markov model sequence similarity networks. Functional identification of these candidate species is much more difficult; we discuss one biochemical approach.

Trm13p, the tRNA:Xm4 modification enzyme from Saccharomyces cerevisiae is a member of the Rossmann-fold MTase superfamily: prediction of structure and active site

2'-O-ribose methylation is one of the most common posttranscriptional modifications in RNA. Methylations at different positions are introduced by enzymes from at least two unrelated superfamilies. Recently, a new family of eukaryotic RNA methyltransferases (MTases) has been identified, and its representative from yeast (Yol125w, renamed as Trm13p) has been shown to 2'-O-methylate position 4 of tRNA. Trm13 is conserved in Eukaryota, but exhibits no sequence similarity to other known MTases. Here, I present the results of bioinformatics analysis which suggest that Trm13 is a strongly diverged member of the Rossmann-fold MTase (RFM) superfamily, and therefore is evolutionarily related to 2'-O-MTases such as Trm7 and fibrillarin. However, the character of conserved residues in the predicted active site of the Trm13 family suggests it may use a different mechanism of ribose methylation than its relatives. A molecular model of the Trm13p structure has been constructed and evaluated for potential accuracy using model quality assessment methods. The predicted structure will facilitate experimental analyses of the Trm13p mechanism of action.

Structure of the thiostrepton resistance methyltransferase.S-adenosyl-L-methionine complex and its interaction with ribosomal RNA

The x-ray crystal structure of the thiostrepton resistance RNA methyltransferase (Tsr).S-adenosyl-L-methionine (AdoMet) complex was determined at 2.45-A resolution. Tsr is definitively confirmed as a Class IV methyltransferase of the SpoU family with an N-terminal "L30-like" putative target recognition domain. The structure and our in vitro analysis of the interaction of Tsr with its target domain from 23 S ribosomal RNA (rRNA) demonstrate that the active biological unit is a Tsr homodimer. In vitro methylation assays show that Tsr activity is optimal against a 29-nucleotide hairpin rRNA though the full 58-nucleotide L11-binding domain and intact 23 S rRNA are also effective substrates. Molecular docking experiments predict that Tsr.rRNA binding is dictated entirely by the sequence and structure of the rRNA hairpin containing the A1067 target nucleotide and is most likely driven primarily by large complementary electrostatic surfaces. One L30-like domain is predicted to bind the target loop and the other is near an internal loop more distant from the target site where a nucleotide change (U1061 to A) also decreases methylation by Tsr. Furthermore, a predicted interaction with this internal loop by Tsr amino acid Phe-88 was confirmed by mutagenesis and RNA binding experiments. We therefore propose that Tsr achieves its absolute target specificity using the N-terminal domains of each monomer in combination to recognize the two distinct structural elements of the target rRNA hairpin such that both Tsr subunits contribute directly to the positioning of the target nucleotide on the enzyme.

Mutation of a gene essential for ribosome biogenesis, EMG1, causes Bowen-Conradi syndrome

Bowen-Conradi syndrome (BCS) is an autosomal-recessive disorder characterized by severely impaired prenatal and postnatal growth, profound psychomotor retardation, and death in early childhood. Nearly all reported BCS cases have been among Hutterites, with an estimated birth prevalence of 1/355. We previously localized the BCS gene to a 1.9 Mbp interval on human chromosome 12p13.3. The 59 genes in this interval were ranked as candidates for BCS, and 35 of these, including all of the best candidates, were sequenced. We identified variant NM_006331.6:c.400A-->G, p.D86G in the 18S ribosome assembly protein EMG1 as the probable cause of BCS. This mutation segregated with disease, was not found in 414 non-Hutterite alleles, and altered a highly conserved aspartic acid (D) residue. A structural model of human EMG1 suggested that the D86 residue formed a salt bridge with arginine 84 that would be disrupted by the glycine (G) substitution. EMG1 mRNA was detected in all human adult and fetal tissues tested. In BCS patient fibroblasts, EMG1 mRNA levels did not differ from those of normal cells, but EMG1 protein was dramatically reduced in comparison to that of normal controls. In mammalian cells, overexpression of EMG1 harboring the D86G mutation decreased the level of soluble EMG1 protein, and in yeast two-hybrid analysis, the D86G substitution increased interaction between EMG1 subunits. These findings suggested that the D-to-G mutation caused aggregation of EMG1, thereby reducing the level of the protein and causing BCS.

YgdE is the 2'-O-ribose methyltransferase RlmM specific for nucleotide C2498 in bacterial 23S rRNA

The rRNAs of Escherichia coli contain four 2'-O-methylated nucleotides. Similar to other bacterial species and in contrast with Archaea and Eukaryota, the E. coli rRNA modifications are catalysed by specific methyltransferases that find their nucleotide targets without being guided by small complementary RNAs. We show here that the ygdE gene encodes the methyltransferase that catalyses 2'-O-methylation at nucleotide C2498 in the peptidyl transferase loop of E. coli 23S rRNA. Analyses of rRNAs using MALDI mass spectrometry showed that inactivation of the ygdE gene leads to loss of methylation at nucleotide C2498. The loss of ygdE function causes a slight reduction in bacterial fitness. Methylation at C2498 was restored by complementing the knock-out strain with a recombinant copy of ygdE. The recombinant YgdE methyltransferase modifies C2498 in naked 23S rRNA, but not in assembled 50S subunits or ribosomes. Nucleotide C2498 is situated within a highly conserved and heavily modified rRNA sequence, and YgdE's activity is influenced by other modification enzymes that target this region. Phylogenetically, YgdE is placed in the cluster of orthologous groups COG2933 together with S-adenosylmethionine-dependent, Rossmann-fold methyltransferases such as the archaeal and eukaryotic RNA-guided fibrillarins. The ygdE gene has been redesignated rlmM for rRNA large subunit methyltransferase M.

RNomics and Modomics in the halophilic archaea Haloferax volcanii: identification of RNA modification genes

Background

Naturally occurring RNAs contain numerous enzymatically altered nucleosides. Differences in RNA populations (RNomics) and pattern of RNA modifications (Modomics) depends on the organism analyzed and are two of the criteria that distinguish the three kingdoms of life. If the genomic sequences of the RNA molecules can be derived from whole genome sequence information, the modification profile cannot and requires or direct sequencing of the RNAs or predictive methods base on the presence or absence of the modifications genes.

Results

By employing a comparative genomics approach, we predicted almost all of the genes coding for the t+rRNA modification enzymes in the mesophilic moderate halophile Haloferax volcanii. These encode both guide RNAs and enzymes. Some are orthologous to previously identified genes in Archaea, Bacteria or in Saccharomyces cerevisiae, but several are original predictions.

Conclusion

The number of modifications in t+rRNAs in the halophilic archaeon is surprisingly low when compared with other Archaea or Bacteria, particularly the hyperthermophilic organisms. This may result from the specific lifestyle of halophiles that require high intracellular salt concentration for survival. This salt content could allow RNA to maintain its functional structural integrity with fewer modifications. We predict that the few modifications present must be particularly important for decoding, accuracy of translation or are modifications that cannot be functionally replaced by the electrostatic interactions provided by the surrounding salt-ions. This analysis also guides future experimental validation work aiming to complete the understanding of the function of RNA modifications in Archaeal translation.

YbeA is the m3Psi methyltransferase RlmH that targets nucleotide 1915 in 23S rRNA

Pseudouridines in the stable RNAs of Bacteria are seldom subjected to further modification. There are 11 pseudouridine (Psi) sites in Escherichia coli rRNA, and further modification is found only at Psi1915 in 23S rRNA, where the N-3 position of the base becomes methylated. Here, we report the identity of the E. coli methyltransferase that specifically catalyzes methyl group addition to form m(3)Psi1915. Analyses of E. coli rRNAs using MALDI mass spectrometry showed that inactivation of the ybeA gene leads to loss of methylation at nucleotide Psi1915. Methylation is restored by complementing the knockout strain with a plasmid-encoded copy of ybeA. Homologs of the ybeA gene, and thus presumably the ensuing methylation at nucleotide m(3)Psi1915, are present in most bacterial lineages but are essentially absent in the Archaea and Eukaryota. Loss of ybeA function in E. coli causes a slight slowing of the growth rate. Phylogenetically, ybeA and its homologs are grouped with other putative S-adenosylmethionine-dependent, SPOUT methyltransferase genes in the Cluster of Orthologous Genes COG1576; ybeA is the first member to be functionally characterized. The YbeA methyltransferase is active as a homodimer and docks comfortably into the ribosomal A site without encroaching into the P site. YbeA makes extensive interface contacts with both the 30S and 50S subunits to align its active site cofactor adjacent to nucleotide Psi1915. Methylation by YbeA (redesignated RlmH for rRNA large subunit methyltransferase H) possibly functions as a stamp of approval signifying that the 50S subunit has engaged in translational initiation.

Identification of pseudouridine methyltransferase in Escherichia coli

In ribosomal RNA, modified nucleosides are found in functionally important regions, but their function is obscure. Stem-loop 69 of Escherichia coli 23S rRNA contains three modified nucleosides: pseudouridines at positions 1911 and 1917, and N3 methyl-pseudouridine (m(3)Psi) at position 1915. The gene for pseudouridine methyltransferase was previously not known. We identified E. coli protein YbeA as the methyltransferase methylating Psi1915 in 23S rRNA. The E. coli ybeA gene deletion strain lacks the N3 methylation at position 1915 of 23S rRNA as revealed by primer extension and nucleoside analysis by HPLC. Methylation at position 1915 is restored in the ybeA deletion strain when recombinant YbeA protein is expressed from a plasmid. In addition, we show that purified YbeA protein is able to methylate pseudouridine in vitro using 70S ribosomes but not 50S subunits from the ybeA deletion strain as substrate. Pseudouridine is the preferred substrate as revealed by the inability of YbeA to methylate uridine at position 1915. This shows that YbeA is acting at the final stage during ribosome assembly, probably during translation initiation. Hereby, we propose to rename the YbeA protein to RlmH according to uniform nomenclature of RNA methyltransferases. RlmH belongs to the SPOUT superfamily of methyltransferases. RlmH was found to be well conserved in bacteria, and the gene is present in plant and in several archaeal genomes. RlmH is the first pseudouridine specific methyltransferase identified so far and is likely to be the only one existing in bacteria, as m(3)Psi1915 is the only methylated pseudouridine in bacteria described to date.

YccW is the m5C methyltransferase specific for 23S rRNA nucleotide 1962

Methylation at the 5-position of cytosine [m(5)C (5-methylcytidine)] occurs at three RNA nucleotides in Escherichia coli. All these modifications are at highly conserved nucleotides in the rRNAs, and each is catalyzed by its own m(5)C methyltransferase enzyme. Two of the enzymes, RsmB and RsmF, are already known and methylate 16S rRNA at nucleotides C967 and C1407, respectively. Here, we report the identity of the third E. coli m(5)C methyltransferase. Analysis of rRNAs by matrix-assisted laser desorption/ionization mass spectrometry showed that inactivation of the yccW gene leads to loss of m(5)C methylation at nucleotide 1962 in E. coli 23S rRNA. This methylation is restored by complementing the knockout strain with a plasmid-encoded copy of the yccW gene. Purified recombinant YccW protein retains its specificity for C1962 in vitro and methylates naked 23S rRNA isolated from the yccW knockout strain. However, YccW does not methylate assembled 50S subunits, and this is somewhat surprising as the published crystal structures show nucleotide C1962 to be fully accessible at the subunit interface. YccW-directed methylation at nucleotide C1962 is conserved in bacteria, and loss of this methylation in E. coli marginally reduces its growth rate. YccW had previously eluded identification because it displays only limited sequence similarity to the m(5)C methyltransferases RsmB and RsmF and is in fact more similar to known m(5)U (5-methyluridine) RNA methyltransferases. In keeping with the previously proposed nomenclature system for bacterial rRNA methyltransferases, yccW is now designated as the rRNA large subunit methyltransferase gene rlmI.

Stabilization of tRNA (mG37) methyltransferase [TrmD] from Aquifex aeolicus by an intersubunit disulfide bond formation

Recombinant Aquifex aeolicus TrmD protein has a Cys20-Cys20 disulfide bond between its two subunits. This was demonstrated by SDS-polyacrylamide gel analysis of wild-type enzyme and C20S mutant protein (in which the Cys20 residue is substituted by serine), in the absence or presence of various concentrations of dithiothreitol. Analytical gel-filtration chromatography revealed that the C20S mutant protein forms a dimer structure even though it is missing the disulfide bond. Western blotting analysis suggests that the Cys20-Cys20 disulfide bond is formed in native TrmD protein in living A. aeolicus cells. Incubation at 85 degrees C for 20 min caused the precipitation of more than half of the C20S protein, while more than 70% of the wild-type enzyme was soluble at that temperature. This assay clearly demonstrates that the disulfide bond enhances the protein stability at 85 degrees C. A kinetic assay showed that the methyl-transfer activity of the C20S mutant protein was slightly less than that of the wild-type enzyme at 70 degrees C. Comparison of the CD-spectra of wild-type and C20S proteins reveals that some of the alpha-helices in the C20S mutant protein are less tightly packed than those of the wild-type enzyme at 70 degrees C.

Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses

For a very long time, Type II restriction enzymes (REases) have been a paradigm of ORFans: proteins with no detectable similarity to each other and to any other protein in the database, despite common cellular and biochemical function. Crystallographic analyses published until January 2008 provided high-resolution structures for only 28 of 1637 Type II REase sequences available in the Restriction Enzyme database (REBASE). Among these structures, all but two possess catalytic domains with the common PD-(D/E)XK nuclease fold. Two structures are unrelated to the others: R.BfiI exhibits the phospholipase D (PLD) fold, while R.PabI has a new fold termed 'half-pipe'. Thus far, bioinformatic studies supported by site-directed mutagenesis have extended the number of tentatively assigned REase folds to five (now including also GIY-YIG and HNH folds identified earlier in homing endonucleases) and provided structural predictions for dozens of REase sequences without experimentally solved structures. Here, we present a comprehensive study of all Type II REase sequences available in REBASE together with their homologs detectable in the nonredundant and environmental samples databases at the NCBI. We present the summary and critical evaluation of structural assignments and predictions reported earlier, new classification of all REase sequences into families, domain architecture analysis and new predictions of three-dimensional folds. Among 289 experimentally characterized (not putative) Type II REases, whose apparently full-length sequences are available in REBASE, we assign 199 (69%) to contain the PD-(D/E)XK domain. The HNH domain is the second most common, with 24 (8%) members. When putative REases are taken into account, the fraction of PD-(D/E)XK and HNH folds changes to 48% and 30%, respectively. Fifty-six characterized (and 521 predicted) REases remain unassigned to any of the five REase folds identified so far, and may exhibit new architectures. These enzymes are proposed as the most interesting targets for structure determination by high-resolution experimental methods. Our analysis provides the first comprehensive map of sequence-structure relationships among Type II REases and will help to focus the efforts of structural and functional genomics of this large and biotechnologically important class of enzymes.

The crystal structure of Nep1 reveals an extended SPOUT-class methyltransferase fold and a pre-organized SAM-binding site

Ribosome biogenesis in eukaryotes requires the participation of a large number of ribosome assembly factors. The highly conserved eukaryotic nucleolar protein Nep1 has an essential but unknown function in 18S rRNA processing and ribosome biogenesis. In Saccharomyces cerevisiae the malfunction of a temperature-sensitive Nep1 protein (nep1-1(ts)) was suppressed by the addition of S-adenosylmethionine (SAM). This suggests the participation of Nep1 in a methyltransferase reaction during ribosome biogenesis. In addition, yeast Nep1 binds to a 6-nt RNA-binding motif also found in 18S rRNA and facilitates the incorporation of ribosomal protein Rps19 during the formation of pre-ribosomes. Here, we present the X-ray structure of the Nep1 homolog from the archaebacterium Methanocaldococcus jannaschii in its free form (2.2 A resolution) and bound to the S-adenosylmethionine analog S-adenosylhomocysteine (SAH, 2.15 A resolution) and the antibiotic and general methyltransferase inhibitor sinefungin (2.25 A resolution). The structure reveals a fold which is very similar to the conserved core fold of the SPOUT-class methyltransferases but contains a novel extension of this common core fold. SAH and sinefungin bind to Nep1 at a preformed binding site that is topologically equivalent to the cofactor-binding site in other SPOUT-class methyltransferases. Therefore, our structures together with previous genetic data suggest that Nep1 is a genuine rRNA methyltransferase.

The yeast ribosome synthesis factor Emg1 is a novel member of the superfamily of alpha/beta knot fold methyltransferases

Emg1 was previously shown to be required for maturation of the 18S rRNA and biogenesis of the 40S ribosomal subunit. Here we report the determination of the crystal structure of Emg1 at 2 Å resolution in complex with the methyl donor, S-adenosyl-methionine (SAM). This structure identifies Emg1 as a novel member of the alpha/beta knot fold methyltransferase (SPOUT) superfamily. In addition to the conserved SPOUT core, Emg1 has two unique domains that form an extended surface, which we predict to be involved in binding of RNA substrates. A point mutation within a basic patch on this surface almost completely abolished RNA binding in vitro. Three point mutations designed to disrupt the interaction of Emg1 with SAM each caused>100-fold reduction in SAM binding in vitro. Expression of only Emg1 with these mutations could support growth and apparently normal ribosome biogenesis in strains genetically depleted of Emg1. We conclude that the catalytic activity of Emg1 is not essential and that the presence of the protein is both necessary and sufficient for ribosome biogenesis.

Crystal structure and mutational study of a unique SpoU family archaeal methylase that forms 2'-O-methylcytidine at position 56 of tRNA

The conserved cytidine residue at position 56 of tRNA contributes to the maintenance of the L-shaped tertiary structure. aTrm56 catalyzes the 2'-O-methylation of the cytidine residue in archaeal tRNA, using S-adenosyl-L-methionine. Based on the amino acid sequence, aTrm56 is the most distant member of the SpoU family. Here, we determined the crystal structure of Pyrococcus horikoshii aTrm56 complexed with S-adenosyl-L-methionine at 2.48 A resolution. aTrm56 consists of the SPOUT domain, which contains the characteristic deep trefoil knot, and a unique C-terminal beta-hairpin. aTrm56 forms a dimer. The S-adenosyl-L-methionine binding and dimerization of aTrm56 were similar to those of the other SpoU members. A structure-based sequence alignment revealed that aTrm56 conserves only motif II, among the four signature motifs. However, an essential Arg16 residue is located at a novel position within motif I. Biochemical assays showed that aTrm56 prefers the L-shaped tRNA to the lambda form as its substrate.