An evolutionary approach uncovers a diverse response of tRNA 2-thiolation to elevated temperatures in yeast

Chemical modifications of transfer RNA (tRNA) molecules are evolutionarily well conserved and critical for translation and tRNA structure. Little is known how these nucleoside modifications respond to physiological stress. Using mass spectrometry and complementary methods, we defined tRNA modification levels in six yeast species in response to elevated temperatures. We show that 2-thiolation of uridine at position 34 (s²U₃₄) is impaired at temperatures exceeding 30°C in the commonly used Saccharomyces cerevisiae laboratory strains S288C and W303, and in Saccharomyces bayanus. Upon stress relief, thiolation levels recover and we find no evidence that modified tRNA or s²U₃₄ nucleosides are actively removed. Our results suggest that loss of 2-thiolation follows accumulation of newly synthesized tRNA that lack s²U₃₄ modification due to temperature sensitivity of the URM1 pathway in S. cerevisiae and S. bayanus. Furthermore, our analysis of the tRNA modification pattern in selected yeast species revealed two alternative phenotypes. Most strains moderately increase their tRNA modification levels in response to heat, possibly constituting a common adaptation to high temperatures. However, an overall reduction of nucleoside modifications was observed exclusively in S288C. This surprising finding emphasizes the importance of studies that utilize the power of evolutionary biology, and highlights the need for future systematic studies on tRNA modifications in additional model organisms.     

Complete modification maps for the cytosolic small and large subunit rRNAs of Euglena gracilis: functional and evolutionary implications of contrasting patterns between the two rRNA components

In the protist Euglena gracilis, the cytosolic small subunit (SSU) rRNA is a single, covalently continuous species typical of most eukaryotes; in contrast, the large subunit (LSU) rRNA is naturally fragmented, comprising 14 separate RNA molecules instead of the bipartite (28S + 5.8S) eukaryotic LSU rRNA typically seen. We present extensively revised secondary structure models of the E. gracilis SSU and LSU rRNAs and have mapped the positions of all of the modified nucleosides in these rRNAs (88 in SSU rRNA and 262 in LSU rRNA, with only 3 LSU rRNA modifications incompletely characterized). The relative proportions of ribose-methylated nucleosides and pseudouridine (∼ 60% and ∼ 35%, respectively) are closely similar in the two rRNAs; however, whereas the Euglena SSU rRNA has about the same absolute number of modifications as its human counterpart, the Euglena LSU rRNA has twice as many modifications as the corresponding human LSU rRNA. The increased levels of rRNA fragmentation and modification in E. gracilis LSU rRNA are correlated with a 3-fold increase in the level of mispairing in helical regions compared to the human LSU rRNA. In contrast, no comparable increase in mispairing is seen in helical regions of the SSU rRNA compared to its homologs in other eukaryotes. In view of the reported effects of both ribose-methylated nucleoside and pseudouridine residues on RNA structure, these correlations lead us to suggest that increased modification in the LSU rRNA may play a role in stabilizing a ‘looser’ structure promoted by elevated helical mispairing and a high degree of fragmentation.     

The open reading frame TTC1157 of Thermus thermophilus HB27 encodes the methyltransferase forming N²-methylguanosine at position 6 in tRNA

N(2)-methylguanosine (m(2)G) is found at position 6 in the acceptor stem of Thermus thermophilus tRNA(Phe). In this article, we describe the cloning, expression, and characterization of the T. thermophilus HB27 methyltransferase (MTase) encoded by the TTC1157 open reading frame that catalyzes the formation of this modified nucleoside. S-adenosyl-L-methionine is used as donor of the methyl group. The enzyme behaves as a monomer in solution. It contains an N-terminal THUMP domain predicted to bind RNA and contains a C-terminal Rossmann-fold methyltransferase (RFM) domain predicted to be responsible for catalysis. We propose to rename the TTC1157 gene trmN and the corresponding protein TrmN, according to the bacterial nomenclature of tRNA methyltransferases. Inactivation of the trmN gene in the T. thermophilus HB27 chromosome led to a total absence of m(2)G in tRNA but did not affect cell growth or the formation of other modified nucleosides in tRNA(Phe). Archaeal homologs of TrmN were identified and characterized. These proteins catalyze the same reaction as TrmN from T. thermophilus. Individual THUMP and RFM domains of PF1002 from Pyrococcus furiosus were produced. These separate domains were inactive and did not bind tRNA, reinforcing the idea that the THUMP domain acts in concert with the catalytic domain to target a particular position of the tRNA molecule.     

Synthesis of RNA containing O-beta-D-ribofuranosyl-(1'-2')-adenosine-5'-phosphate and 1-methyladenosine, minor components of tRNA

tRNA is best known for its function as amino acid carrier in the translation process, using the anticodon loop in the recognition process with mRNA. However, the impact of tRNA on cell function is much wider, and mutations in tRNA can lead to a broad range of diseases. Although the cloverleaf structure of tRNA is well-known based on X-ray-diffraction studies, little is known about the dynamics of this fold, the way structural dynamics of tRNA is influenced by the modified nucleotides present in tRNA, and their influence on the recognition of tRNA by synthetases, ribosomes, and other biomolecules. One of the reasons for this is the lack of good synthetic methods to incorporate modified nucleotides in tRNA so that larger amounts become available for NMR studies. Except of 2'-O-methylated nucleosides, only one other sugar-modified nucleoside is present in tRNA, i.e., 2'-O-beta-D-ribofuranosyl nucleosides. The T loop of tRNA often contains charged modified nucleosides, of which 1-methyladenosine and phosphorylated disaccharide nucleosides are striking examples. A protecting-group strategy was developed to introduce 1-methyladenosine and 5'-O-phosphorylated 2'-O-(beta-D-ribofuranosyl)-beta-D-ribofuranosyladenine in the same RNA fragment. The phosphorylation of the disaccharide nucleoside was performed after the assembly of the RNA on solid support. The modified RNA was characterized by mass-spectrometry analysis from the RNase T1 digestion fragments. The successful synthesis of this T loop of the tRNA of Schizosaccharomyces pombe initiator tRNA(Met) will be followed by its structural analysis by NMR and by studies on the influence of these modified nucleotides on dynamic interactions within the complete tRNA.     

A fungal anticodon nuclease ribotoxin exploits a secondary cleavage site to evade tRNA repair

PaOrf2 and γ-toxin subunits of Pichia acaciae toxin (PaT) and Kluyveromyces lactis zymocin are tRNA anticodon nucleases. These secreted ribotoxins are assimilated by Saccharomyces cerevisiae, wherein they arrest growth by depleting specific tRNAs. Toxicity can be recapitulated by induced intracellular expression of PaOrf2 or γ-toxin in S. cerevisiae. Mutational analysis of γ-toxin has identified amino acids required for ribotoxicity in vivo and RNA transesterification in vitro. Here, we report that PaOrf2 residues Glu9 and His287 (putative counterparts of γ-toxin Glu9 and His209) are essential for toxicity. Our results suggest a similar basis for RNA transesterification by PaOrf2 and γ-toxin, despite their dissimilar primary structures and distinctive tRNA target specificities. PaOrf2 makes two sequential incisions in tRNA, the first of which occurs 3' from the mcm(5)s(2)U wobble nucleoside and depends on mcm(5). A second incision two nucleotides upstream results in the net excision of a di-nucleotide. Expression of phage and plant tRNA repair systems can relieve PaOrf2 toxicity when tRNA cleavage is restricted to the secondary site in elp3 cells that lack the mcm(5) wobble U modification. Whereas the endogenous yeast tRNA ligase Trl1 can heal tRNA halves produced by PaOrf2 cleavage in elp3 cells, its RNA sealing activity is inadequate to complete the repair. Compatible sealing activity can be provided in trans by plant tRNA ligase. The damage-rescuing ability of tRNA repair systems is lost when PaOrf2 can break tRNA at both sites. These results highlight the logic of a two-incision mechanism of tRNA anticodon damage that evades productive repair by tRNA ligases.     

The Pivotal Role of the Mitochondrial Amidoxime Reducing Component 2 in Protecting Human Cells against Apoptotic Effects of the Base Analog N6-Hydroxylaminopurine

N-Hydroxylated nucleobases and nucleosides as N-hydroxylaminopurine (HAP) or N-hydroxyadenosine (HAPR) may be generated endogenously in the course of cell metabolism by cytochrome P450, by oxidative stress or by a deviating nucleotide biosynthesis. These compounds have shown to be toxic and mutagenic for procaryotic and eucaryotic cells. For DNA replication fidelity it is therefore of great importance that organisms exhibit effective mechanisms to remove such non-canonical base analogs from DNA precursor pools. In vitro, the molybdoenzymes mitochondrial amidoxime reducing component 1 and 2 (mARC1 and mARC2) have shown to be capable of reducing N-hydroxylated base analogs and nucleoside analogs to the corresponding canonical nucleobases and nucleosides upon reconstitution with the electron transport proteins cytochrome b₅ and NADH-cytochrome b₅ reductase. By RNAi-mediated down-regulation of mARC in human cell lines the mARC-dependent N-reductive detoxication of HAP in cell metabolism could be demonstrated. For HAPR, on the other hand, the reduction to adenosine seems to be of less significance in the detoxication pathway of human cells as HAPR is primarily metabolized to inosine by direct dehydroxylamination catalyzed by adenosine deaminase. Furthermore, the effect of mARC knockdown on sensitivity of human cells to HAP was examined by flow cytometric quantification of apoptotic cell death and detection of poly (ADP-ribose) polymerase (PARP) cleavage. mARC2 was shown to protect HeLa cells against the apoptotic effects of the base analog, whereas the involvement of mARC1 in reductive detoxication of HAP does not seem to be pivotal.     

Familial dysautonomia (FD) patients have reduced levels of the modified wobble nucleoside mcmsU in tRNA

Familial dysautonomia (FD) is a recessive neurodegenerative genetic disease. FD is caused by a mutation in the IKBKAP gene resulting in a splicing defect and reduced levels of full length IKAP protein. IKAP homologues can be found in all eukaryotes and are part of a conserved six subunit protein complex, Elongator complex. Inactivation of any Elongator subunit gene in multicellular organisms cause a wide range of phenotypes, suggesting that Elongator has a pivotal role in several cellular processes. In yeast, there is convincing evidence that the main role of Elongator complex is in formation of modified wobble uridine nucleosides in tRNA and that their absence will influence translational efficiency. To date, no study has explored the possibility that FD patients display defects in formation of modified wobble uridine nucleosides as a consequence of reduced IKAP levels. In this study, we show that brain tissue and fibroblast cell lines from FD patients have reduced levels of the wobble uridine nucleoside 5-methoxycarbonylmethyl-2-thiouridine (mcm⁵s²U). Our findings indicate that FD could be caused by inefficient translation due to lower levels of wobble uridine nucleosides.     

Dynamics of Recognition between tRNA and elongation factor Tu

Elongation factor Tu (EF-Tu) binds to all standard aminoacyl transfer RNAs (aa-tRNAs) and transports them to the ribosome while protecting the ester linkage between the tRNA and its cognate amino acid. We use molecular dynamics simulations to investigate the dynamics of the EF-Tu·guanosine 5′-triphosphate·aa-tRNA^Cys complex and the roles played by Mg²⁺ ions and modified nucleosides on the free energy of protein·RNA binding. Individual modified nucleosides have pronounced effects on the structural dynamics of tRNA and the EF-Tu·Cys-tRNA^Cys interface. Combined energetic and evolutionary analyses identify the coevolution of residues in EF-Tu and aa-tRNAs at the binding interface. Highly conserved EF-Tu residues are responsible for both attracting aa-tRNAs as well as providing nearby nonbonded repulsive energies that help fine-tune molecular attraction at the binding interface. In addition to the 3′ CCA end, highly conserved tRNA nucleotides G1, G52, G53, and U54 contribute significantly to EF-Tu binding energies. Modification of U54 to thymine affects the structure of the tRNA common loop resulting in a change in binding interface contacts. In addition, other nucleotides, conserved within certain tRNA specificities, may be responsible for tuning aa-tRNA binding to EF-Tu. The trend in EF-Tu·Cys-tRNA^Cys binding energies observed as the result of mutating the tRNA agrees with experimental observation. We also predict variations in binding free energies upon misacylation of tRNA^Cys with d-cysteine or O-phosphoserine and upon changing the protonation state of l-cysteine. Principal components analysis in each case reveals changes in the communication network across the protein·tRNA interface and is the basis for the entropy calculations.     

Role of forefinger and thumb loops in production of Ψ54 and Ψ55 in tRNAs by archaeal Pus10

Pseudouridines (Ψ) are found in structurally and functionally important regions of RNAs. Six families of Ψ synthases, TruA, TruB, TruD, RsuA, RluA, and Pus10 have been identified. Pus10 is present in Archaea and Eukarya. While most archaeal Pus10 produce both tRNA Ψ54 and Ψ55, some produce only Ψ55. Interestingly, human PUS10 has been implicated in apoptosis and Crohn’s and Celiac diseases. Homology models of archaeal Pus10 proteins based on the crystal structure of human PUS10 reveal that there are subtle structural differences in all of these Pus10 proteins. These observations suggest that structural changes in homologous proteins may lead to loss, gain, or change of their functions, warranting the need to study the structure-function relationship of these proteins. Using comparison of structural models and a series of mutations, we identified forefinger loop (reminiscent of that of RluA) and an Arg and a Tyr residue of archaeal Pus10 as critical determinants for its Ψ54, but not for its Ψ55 activity. We also found that a Leu residue, in addition to the catalytic Asp, is essential for both activities. Since forefinger loop is needed for both rRNA and tRNA Ψ synthase activities of RluA, but only for tRNA Ψ54 activity of Pus10, archaeal Pus10 proteins must use a different mechanism of recognition for Ψ55 activity. We propose that archaeal Pus10 uses two distinct mechanisms for substrate uridine recognition and binding. However, since we did not observe any mutation that affected only Ψ55 activity, both mechanisms for archaeal Pus10 activities must share some common features.     

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³Ψ) at position 1915. The gene for pseudouridine methyltransferase was previously not known. We identified E. coli protein YbeA as the methyltransferase methylating Ψ1915 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³Ψ1915 is the only methylated pseudouridine in bacteria described to date.     

Solution Structure of an Archaeal RNase P Binary Protein Complex: Formation of the 30-kDa Complex between Pyrococcus furiosus RPP21 and RPP29 Is Accompanied by Coupled Protein Folding and Highlights Critical Features for Protein-Protein and Protein-RNA Interactions

Ribonuclease P (RNase P) is a ribonucleoprotein (RNP) enzyme that catalyzes the Mg²⁺-dependent 5′ maturation of precursor tRNAs. In all domains of life, it is a ribozyme: the RNase P RNA (RPR) component has been demonstrated to be responsible for catalysis. However, the number of RNase P protein subunits (RPPs) varies from 1 in bacteria to 9 or 10 in eukarya. The archaeal RPR is associated with at least 4 RPPs, which function in pairs (RPP21-RPP29 and RPP30-POP5). We used solution NMR spectroscopy to determine the three-dimensional structure of the protein-protein complex comprising Pyrococcus furiosus RPP21 and RPP29. We found that the protein-protein interaction is characterized by coupled folding of secondary structural elements that participate in interface formation. In addition to detailing the intermolecular contacts that stabilize this 30-kDa binary complex, the structure identifies surfaces rich in conserved basic residues likely vital for recognition of the RPR and/or precursor tRNA. Furthermore, enzymatic footprinting experiments allowed us to localize the RPP21-RPP29 complex to the specificity domain of the RPR. These findings provide valuable new insights into mechanisms of RNP assembly and serve as important steps towards a three-dimensional model of this ancient RNP enzyme.     

The Escherichia coli RlmN methyltransferase is a dual-specificity enzyme that modifies both rRNA and tRNA and controls translational accuracy

Modifying RNA enzymes are highly specific for substrate-rRNA or tRNA-and the target position. In Escherichia coli, there are very few multisite acting enzymes, and only one rRNA/tRNA dual-specificity enzyme, pseudouridine synthase RluA, has been identified to date. Among the tRNA-modifying enzymes, the methyltransferase responsible for the m(2)A synthesis at purine 37 in a tRNA set still remains unknown. m(2)A is also present at position 2503 in the peptidyl transferase center of 23S RNA, where it is introduced by RlmN, a radical S-adenosyl-L-methionine (SAM) enzyme. Here, we show that E. coli RlmN is a dual-specificity enzyme that catalyzes methylation of both rRNA and tRNA. The ΔrlmN mutant lacks m(2)A in both RNA types, whereas the expression of recombinant RlmN from a plasmid introduced into this mutant restores tRNA modification. Moreover, RlmN performs m(2)A(37) synthesis in vitro using a tRNA chimera as a substrate. This chimera has also proved useful to characterize some tRNA identity determinants for RlmN and other tRNA modification enzymes. Our data suggest that RlmN works in a late step during tRNA maturation by recognizing a precise 3D structure of tRNA. RlmN inactivation increases the misreading of a UAG stop codon. Since loss of m(2)A(37) from tRNA is expected to produce a hyperaccurate phenotype, we believe that the error-prone phenotype exhibited by the ΔrlmN mutant is due to loss of m(2)A from 23S rRNA and, accordingly, that the m(2)A2503 modification plays a crucial role in the proofreading step occurring at the peptidyl transferase center.     

Crystal structure of Thermus thermophilus tRNA m1A58 methyltransferase and biophysical characterization of its interaction with tRNA

Methyltransferases from the m¹A₅₈ tRNA methyltransferase (TrmI) family catalyze the S-adenosyl-l-methionine-dependent N₁-methylation of tRNA adenosine 58. The crystal structure of Thermus thermophilus TrmI, in complex with S-adenosyl-l-homocysteine, was determined at 1.7 Å resolution. This structure is closely related to that of Mycobacterium tuberculosis TrmI, and their comparison enabled us to enlighten two grooves in the TrmI structure that are large enough and electrostatically compatible to accommodate one tRNA per face of TrmI tetramer. We have then conducted a biophysical study based on electrospray ionization mass spectrometry, site-directed mutagenesis, and molecular docking. First, we confirmed the tetrameric oligomerization state of TrmI, and we showed that this protein remains tetrameric upon tRNA binding, with formation of complexes involving one to two molecules of tRNA per TrmI tetramer. Second, three key residues for the methylation reaction were identified: the universally conserved D170 and two conserved aromatic residues Y78 and Y194. We then used molecular docking to position a N₉-methyladenine in the active site of TrmI. The N₉-methyladenine snugly fits into the catalytic cleft, where the side chain of D170 acts as a bidentate ligand binding the amino moiety of S-adenosyl-l-methionine and the exocyclic amino group of the adenosine. Y194 interacts with the N₉-methyladenine ring, whereas Y78 can stabilize the sugar ring. From our results, we propose that the conserved residues that form the catalytic cavity (D170, Y78, and Y194) are essential for fashioning an optimized shape of the catalytic pocket.     

Complete analysis of tRNA-modified nucleosides by high-performance liquid chromatography: the 29 modified nucleosides of Salmonella typhimurium and Escherichia coli tRNA

A high-performance liquid chromatography (HPLC) method has been developed to quantify the major and modified nucleoside composition of total, unfractionated transfer RNA. The method is rapid and sensitive and offers a high degree of chromatographic resolution suitable for quantifying both stable and unstable modified nucleosides. It is nondestructive and allows the recovery of nucleosides for further characterization. We apply the method in the analysis of the 29 modified nucleosides in tRNA from Salmonella typhimurium (and Escherichia coli) and show it to be useful in examining changes in the modified nucleoside content of tRNA. Such changes may be important in regulation.     

A genome-wide screen identifies genes required for formation of the wobble nucleoside 5-methoxycarbonylmethyl-2-thiouridine in Saccharomyces cerevisiae

We recently showed that the gamma-subunit of Kluyveromyces lactis killer toxin (gamma-toxin) is a tRNA endonuclease that cleaves tRNA(mcm5s2UUC Glu), tRNA(mcm5s2UUU Lys), and tRNA(mcm5s2UUG Gln) 3' of the wobble nucleoside 5-methoxycarbonylmethyl-2-thiouridine (mcm(5)s(2)U). The 5-methoxycarbonylmethyl (mcm(5)) side chain was important for efficient cleavage by gamma-toxin, and defects in mcm(5) side-chain synthesis correlated with resistance to gamma-toxin. Based on this correlation, a genome-wide screen was performed to identify gene products involved in the formation of the mcm(5) side chain. From a collection of 4826 homozygous diploid Saccharomyces cerevisiae strains, each with one nonessential gene deleted, 63 mutants resistant to Kluyveromyces lactis killer toxin were identified. Among these, eight were earlier identified to have a defect in formation of the mcm(5) side chain. Analysis of the remaining mutants and other known gamma-toxin resistant mutants revealed that sit4, kti14, and KTI5 mutants also have a defect in the formation of mcm(5). A mutant lacking two of the Sit4-associated proteins, Sap185 and Sap190, displays the same modification defect as a sit4-null mutant. Interestingly, several mutants were found to be defective in the synthesis of the 2-thio (s(2)) group of the mcm(5)s(2)U nucleoside. In addition to earlier described mutants, formation of the s(2) group was also abolished in urm1, uba4, and ncs2 mutants and decreased in the yor251c mutant. Like the absence of the mcm(5) side chain, the lack of the s(2) group renders tRNA(mcm5s2UUC Glu) less sensitive to gamma-toxin, reinforcing the importance of the wobble nucleoside mcm(5)s(2)U for tRNA cleavage by gamma-toxin.