Novel base-pairing interactions at the tRNA wobble position crucial for accurate reading of the genetic code

Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field

Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.

Engineering tRNAs for the Ribosomal Translation of Non-proteinogenic Monomers

Ribosome-dependent protein biosynthesis is an essential cellular process mediated by transfer RNAs (tRNAs). Generally, ribosomally synthesized proteins are limited to the 22 proteinogenic amino acids (pAAs: 20 l-α-amino acids present in the standard genetic code, selenocysteine, and pyrrolysine). However, engineering tRNAs for the ribosomal incorporation of non-proteinogenic monomers (npMs) as building blocks has led to the creation of unique polypeptides with broad applications in cellular biology, material science, spectroscopy, and pharmaceuticals. Ribosomal polymerization of these engineered polypeptides presents a variety of challenges for biochemists, as translation efficiency and fidelity is often insufficient when employing npMs. In this Review, we will focus on the methodologies for engineering tRNAs to overcome these issues and explore recent advances both in vitro and in vivo. These efforts include increasing orthogonality, recruiting essential translation factors, and creation of expanded genetic codes. After our review on the biochemical optimizations of tRNAs, we provide examples of their use in genetic code manipulation, with a focus on the in vitro discovery of bioactive macrocyclic peptides containing npMs. Finally, an analysis of the current state of tRNA engineering is presented, along with existing challenges and future perspectives for the field.

Structural insights into the decoding capability of isoleucine tRNAs with lysidine and agmatidine

The anticodon modifications of transfer RNAs (tRNAs) finetune the codon recognition on the ribosome for accurate translation. Bacteria and archaea utilize the modified cytidines, lysidine (L) and agmatidine (agm2C), respectively, in the anticodon of tRNAIle to decipher AUA codon. L and agm2C contain long side chains with polar termini, but their functions remain elusive. Here we report the cryogenic electron microscopy structures of tRNAsIle recognizing the AUA codon on the ribosome. Both modifications interact with the third adenine of the codon via a unique C-A geometry. The side chains extend toward 3' direction of the mRNA, and the polar termini form hydrogen bonds with 2'-OH of the residue 3'-adjacent to the AUA codon. Biochemical analyses demonstrated that AUA decoding is facilitated by the additional interaction between the polar termini of the modified cytidines and 2'-OH of the fourth mRNA residue. We also visualized cyclic N6-threonylcarbamoyladenosine (ct6A), another tRNA modification, and revealed a molecular basis how ct6A contributes to efficient decoding.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Inosine-Induced Base Pairing Diversity during Reverse Transcription

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

Anticodon stem-loop tRNA modifications influence codon decoding and frame maintenance during translation

RNAs are central to protein synthesis, with ribosomal RNA, transfer RNAs and messenger RNAs comprising the core components of the translation machinery. In addition to the four canonical bases (uracil, cytosine, adenine, and guanine) these RNAs contain an array of enzymatically incorporated chemical modifications. Transfer RNAs (tRNAs) are responsible for ferrying amino acids to the ribosome, and are among the most abundant and highly modified RNAs in the cell across all domains of life. On average, tRNA molecules contain 13 post-transcriptionally modified nucleosides that stabilize their structure and enhance function. There is an extensive chemical diversity of tRNA modifications, with over 90 distinct varieties of modifications reported within tRNA sequences. Some modifications are crucial for tRNAs to adopt their L-shaped tertiary structure, while others promote tRNA interactions with components of the protein synthesis machinery. In particular, modifications in the anticodon stem-loop (ASL), located near the site of tRNA:mRNA interaction, can play key roles in ensuring protein homeostasis and accurate translation. There is an abundance of evidence indicating the importance of ASL modifications for cellular health, and in vitro biochemical and biophysical studies suggest that individual ASL modifications can differentially influence discrete steps in the translation pathway. This review examines the molecular level consequences of tRNA ASL modifications in mRNA codon recognition and reading frame maintenance to ensure the rapid and accurate translation of proteins.

Contribution of tRNA sequence and modifications to the decoding preferences of E. coli and M. mycoides tRNAGlyUCC for synonymous glycine codons

tRNA superwobbling, used by certain bacteria and organelles, is an intriguing decoding concept in which a single tRNA isoacceptor is used to decode all synonymous codons of a four-fold degenerate codon box. While Escherichia coli relies on three tRNAGly isoacceptors to decode the four glycine codons (GGN), Mycoplasma mycoides requires only a single tRNAGly. Both organisms express tRNAGly with the anticodon UCC, which are remarkably similar in sequence but different in their decoding ability. By systematically introducing mutations and altering the number and type of tRNA modifications using chemically synthesized tRNAs, we elucidated the contribution of individual nucleotides and chemical groups to decoding by the E. coli and M. mycoides tRNAGly. The tRNA sequence was identified as the key factor for superwobbling, revealing the T-arm sequence as a novel pivotal element. In addition, the presence of tRNA modifications, although not essential for providing superwobbling, was shown to delicately fine-tune and balance the decoding of synonymous codons. This emphasizes that the tRNA sequence and its modifications together form an intricate system of high complexity that is indispensable for accurate and efficient decoding.

Structural aspects of RimP binding on small ribosomal subunit from Staphylococcus aureus

Ribosome biogenesis is an energy-intense multistep process where even minimal defects can cause severe phenotypes up to cell death. Ribosome assembly is facilitated by biogenesis factors such as ribosome assembly factors. These proteins facilitate the interaction of ribosomal proteins with rRNA and correct rRNA folding. One of these maturation factors is RimP which is required for efficient 16S rRNA processing and 30S ribosomal subunit assembly. Here, we describe the binding mode of Staphylococcus aureus RimP to the small ribosomal subunit and present a 4.2 Å resolution cryo-EM reconstruction of the 30S-RimP complex. Together with the solution structure of RimP solved by NMR spectroscopy and RimP-uS12 complex analysis by EPR, DEER, and SAXS approaches, we show the specificity of RimP binding to the 30S subunit from S. aureus. We believe the results presented in this work will contribute to the understanding of the RimP role in the ribosome assembly mechanism.

A connection between the ribosome and two S. pombe tRNA modification mutants subject to rapid tRNA decay

tRNA modifications are crucial in all organisms to ensure tRNA folding and stability, and accurate translation. In both the yeast Saccharomyces cerevisiae and the evolutionarily distant yeast Schizosaccharomyces pombe, mutants lacking certain tRNA body modifications (outside the anticodon loop) are temperature sensitive due to rapid tRNA decay (RTD) of a subset of hypomodified tRNAs. Here we show that for each of two S. pombe mutants subject to RTD, mutations in ribosomal protein genes suppress the temperature sensitivity without altering tRNA levels. Prior work showed that S. pombe trm8Δ mutants, lacking 7-methylguanosine, were temperature sensitive due to RTD, and that one class of suppressors had mutations in the general amino acid control (GAAC) pathway, which was activated concomitant with RTD, resulting in further tRNA loss. We now find that another class of S. pombe trm8Δ suppressors have mutations in rpl genes, encoding 60S subunit proteins, and that suppression occurs with minimal restoration of tRNA levels and reduced GAAC activation. Furthermore, trm8Δ suppression extends to other mutations in the large or small ribosomal subunit. We also find that S. pombe tan1Δ mutants, lacking 4-acetylcytidine, are temperature sensitive due to RTD, that one class of suppressors have rpl mutations, associated with minimal restoration of tRNA levels, and that suppression extends to other rpl and rps mutations. However, although S. pombe tan1Δ temperature sensitivity is associated with some GAAC activation, suppression by an rpl mutation only modestly inhibits GAAC activation. We propose a model in which ribosomal protein mutations result in reduced ribosome concentrations, leading to both reduced ribosome collisions and a reduced requirement for tRNA, with these effects having different relative importance in trm8Δ and tan1Δ mutants. This model is consistent with our results in S. cerevisiae trm8Δ trm4Δ mutants, known to undergo RTD, fueling speculation that this model applies across eukaryotes.

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

A connection between the ribosome and two S. pombe tRNA modification mutants subject to rapid tRNA decay

tRNA modifications are crucial in all organisms to ensure tRNA folding and stability, and accurate translation in the ribosome. In both the yeast Saccharomyces cerevisiae and the evolutionarily distant yeast Schizosaccharomyces pombe, mutants lacking certain tRNA body modifications (outside the anticodon loop) are temperature sensitive due to rapid tRNA decay (RTD) of a subset of hypomodified tRNAs. Here we show that for each of two S. pombe mutants subject to RTD, mutations in ribosomal protein genes suppress the temperature sensitivity without altering tRNA levels. Prior work showed that S. pombe trm8Δ mutants, lacking 7-methylguanosine, were temperature sensitive due to RTD and that one class of suppressors had mutations in the general amino acid control (GAAC) pathway, which was activated concomitant with RTD, resulting in further tRNA loss. We now find that another class of S. pombe trm8Δ suppressors have mutations in rpl genes, encoding 60S subunit proteins, and that suppression occurs with minimal restoration of tRNA levels and reduced GAAC activation. Furthermore, trm8Δ suppression extends to other mutations in the large or small ribosomal subunit. We also find that S. pombe tan1Δ mutants, lacking 4-acetylcytidine, are temperature sensitive due to RTD, that one class of suppressors have rpl mutations, associated with minimal restoration of tRNA levels, and that suppression extends to other rpl and rps mutations. However, although S. pombe tan1Δ temperature sensitivity is associated with some GAAC activation, suppression by an rpl mutation does not significantly inhibit GAAC activation. These results suggest that ribosomal protein mutations suppress the temperature sensitivity of S. pombe trm8Δ and tan1Δ mutants due to reduced ribosome concentrations, leading to both a reduced requirement for tRNA, and reduced ribosome collisions and GAAC activation. Results with S. cerevisiae trm8Δ trm4Δ mutants are consistent with this model, and fuel speculation that similar results will apply across eukaryotes.

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding

Accurate protein synthesis is determined by the two-subunit ribosome's capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site's physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

yaaJ, the tRNA-Specific Adenosine Deaminase, Is Dispensable in Bacillus subtilis

Post-transcriptional modifications of tRNA are crucial for their core function. The inosine (I; 6-deaminated adenosine) at the first position in the anticodon of tRNAArg(ICG) modulates the decoding capability and is generally considered essential for reading CGU, CGC, and CGA codons in eubacteria. We report here that the Bacillus subtilis yaaJ gene encodes tRNA-specific adenosine deaminase and is non-essential for viability. A β-galactosidase reporter assay revealed that the translational activity of CGN codons was not impaired in the yaaJ-deletion mutant. Furthermore, tRNAArg(CCG) responsible for decoding the CGG codon was dispensable, even in the presence or absence of yaaJ. These results strongly suggest that tRNAArg with either the anticodon ICG or ACG has an intrinsic ability to recognize all four CGN codons, providing a fundamental concept of non-canonical wobbling mediated by adenosine and inosine nucleotides in the anticodon. This is the first example of the four-way wobbling by inosine nucleotide in bacterial cells. On the other hand, the absence of inosine modification induced +1 frameshifting, especially at the CGA codon. Additionally, the yaaJ deletion affected growth and competency. Therefore, the inosine modification is beneficial for translational fidelity and proper growth-phase control, and that is why yaaJ has been actually conserved in B. subtilis.

Anionic G•U pairs in bacterial ribosomal rRNAs

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Decoding and Recoding of mRNA Sequences by the Ribosome

Faithful translation of messenger RNA (mRNA) into protein is essential to maintain protein homeostasis in the cell. Spontaneous translation errors are very rare due to stringent selection of cognate aminoacyl transfer RNAs (tRNAs) and the tight control of the mRNA reading frame by the ribosome. Recoding events, such as stop codon readthrough, frameshifting, and translational bypassing, reprogram the ribosome to make intentional mistakes and produce alternative proteins from the same mRNA. The hallmark of recoding is the change of ribosome dynamics. The signals for recoding are built into the mRNA, but their reading depends on the genetic makeup of the cell, resulting in cell-specific changes in expression programs. In this review, I discuss the mechanisms of canonical decoding and tRNA-mRNA translocation; describe alternative pathways leading to recoding; and identify the links among mRNA signals, ribosome dynamics, and recoding.

The thiolation of uridine 34 in tRNA, which controls protein translation, depends on a [4Fe-4S] cluster in the archaeum Methanococcus maripaludis

Thiolation of uridine 34 in the anticodon loop of several tRNAs is conserved in the three domains of life and guarantees fidelity of protein translation. U34-tRNA thiolation is catalyzed by a complex of two proteins in the eukaryotic cytosol (named Ctu1/Ctu2 in humans), but by a single NcsA enzyme in archaea. We report here spectroscopic and biochemical experiments showing that NcsA from Methanococcus maripaludis (MmNcsA) is a dimer that binds a [4Fe-4S] cluster, which is required for catalysis. Moreover, the crystal structure of MmNcsA at 2.8 Å resolution shows that the [4Fe-4S] cluster is coordinated by three conserved cysteines only, in each monomer. Extra electron density on the fourth nonprotein-bonded iron most likely locates the binding site for a hydrogenosulfide ligand, in agreement with the [4Fe-4S] cluster being used to bind and activate the sulfur atom of the sulfur donor. Comparison of the crystal structure of MmNcsA with the AlphaFold model of the human Ctu1/Ctu2 complex shows a very close superposition of the catalytic site residues, including the cysteines that coordinate the [4Fe-4S] cluster in MmNcsA. We thus propose that the same mechanism for U34-tRNA thiolation, mediated by a [4Fe-4S]-dependent enzyme, operates in archaea and eukaryotes.

Tautomerisation Mechanisms in the Adenine-Thymine Nucleobase Pair during DNA Strand Separation

The adenine-thymine tautomer (A*-T*) has previously been discounted as a spontaneous mutagenesis mechanism due to the energetic instability of the tautomeric configuration. We study the stability of A*-T* while the nucleobases undergo DNA strand separation. Our calculations indicate an increase in the stability of A*-T* as the DNA strands unzip and the hydrogen bonds between the bases stretch. Molecular Dynamics simulations reveal the time scales and dynamics of DNA strand separation and the statistical ensemble of opening angles present in a biological environment. Our results demonstrate that the unwinding of DNA, an inherently out-of-equilibrium process facilitated by helicase, will change the energy landscape of the adenine-thymine tautomerization reaction. We propose that DNA strand separation allows the stable tautomerization of adenine-thymine, providing a feasible pathway for genetic point mutations via proton transfer between the A-T bases.

Structural basis for reduced ribosomal A-site fidelity in response to P-site codon-anticodon mismatches

Rapid and accurate translation is essential in all organisms to produce properly folded and functional proteins. mRNA codons that define the protein coding sequences are decoded by tRNAs on the ribosome in the aminoacyl (A) binding site. The mRNA codon and the tRNA anticodon interaction is extensively monitored by the ribosome to ensure accuracy in tRNA selection. While other polymerases that synthesize DNA and RNA can correct for misincorporations, the ribosome is unable to correct mistakes. Instead, when a misincorporation occurs, the mismatched tRNA-mRNA pair moves to the peptidyl (P) site and from this location, causes a reduction in the fidelity at the A site, triggering post-peptidyl transfer quality control. This reduced fidelity allows for additional incorrect tRNAs to be accepted and for release factor 2 (RF2) to recognize sense codons, leading to hydrolysis of the aberrant peptide. Here, we present crystal structures of the ribosome containing a tRNA^Lys in the P site with a U•U mismatch with the mRNA codon. We find that when the mismatch occurs in the second position of the P-site codon-anticodon interaction, the first nucleotide of the A-site codon flips from the mRNA path to engage highly conserved 16S rRNA nucleotide A1493 in the decoding center. We propose that this mRNA nucleotide mispositioning leads to reduced fidelity at the A site. Further, this state may provide an opportunity for RF2 to initiate premature termination before erroneous nascent chains disrupt the cellular proteome.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Structural basis for reduced ribosomal A-site fidelity in response to P-site codon-anticodon mismatches

Rapid and accurate translation is essential in all organisms to produce properly folded and functional proteins. mRNA codons that define the protein coding sequences are decoded by tRNAs on the ribosome in the aminoacyl (A) binding site. The mRNA codon and the tRNA anticodon interaction is extensively monitored by the ribosome to ensure accuracy in tRNA selection. While other polymerases that synthesize DNA and RNA can correct for misincorporations, the ribosome is unable to correct mistakes. Instead, when a misincorporation occurs, the mismatched tRNA-mRNA pair moves to the peptidyl (P) site and from this location, causes a reduction in the fidelity at the A site, triggering post-peptidyl transfer quality control. This reduced fidelity allows for additional incorrect tRNAs to be accepted and for release factor 2 (RF2) to recognize sense codons, leading to hydrolysis of the aberrant peptide. Here, we present crystal structures of the ribosome containing a tRNA ^Lys in the P site with a U•U mismatch with the mRNA codon. We find that when the mismatch occurs in the second position of the P-site codon-anticodon interaction, the first nucleotide of the A-site codon flips from the mRNA path to engage highly conserved 16S rRNA nucleotide A1493 in the decoding center. We propose that this mRNA nucleotide mispositioning leads to reduced fidelity at the A site. Further, this state may provide an opportunity for RF2 to initiate premature termination before erroneous nascent chains disrupt the cellular proteome.

Expansion of the genetic code through reassignment of redundant sense codons using fully modified tRNA

Breaking codon degeneracy for the introduction of non-canonical amino acids offers many opportunities in synthetic biology. Yet, despite the existence of 64 codons, the code has only been expanded to 25 amino acids in vitro. A limiting factor could be the over-reliance on synthetic tRNAs which lack the post-transcriptional modifications that improve translational fidelity. To determine whether modified, wild-type tRNA could improve sense codon reassignment, we developed a new fluorous method for tRNA capture and applied it to the isolation of roughly half of the Escherichia coli tRNA isoacceptors. We then performed codon competition experiments between the five captured wild-type leucyl-tRNAs and their synthetic counterparts, revealing a strong preference for wild-type tRNA in an in vitro translation system. Finally, we compared the ability of wild-type and synthetic leucyl-tRNA to break the degeneracy of the leucine codon box, showing that only captured wild-type tRNAs are discriminated with enough fidelity to accurately split the leucine codon box for the encoding of three separate amino acids. Wild-type tRNAs are therefore enabling reagents for maximizing the reassignment potential of the genetic code.

Translation termination codons in protein synthesis and disease

Fidelity of protein synthesis, a process shaped by several mechanisms involving specialized ribosome regions and external factors, ensures the precise reading of sense as well as stop codons (UGA, UAG, UAA), which are usually localized at the 3' of mRNA and drive the release of the polypeptide chain. However, either natural (NTCs) or premature (PTCs) termination codons, the latter arising from nucleotide changes, can undergo a recoding process named ribosome or translational readthrough, which insert specific amino acids (NTCs) or subset(s) depending on the stop codon type (PTCs). This process is particularly relevant for nonsense mutations, a relatively frequent cause of genetic disorders, which impair gene expression at different levels by potentially leading to mRNA degradation and/or synthesis of truncated proteins. As a matter of fact, many efforts have been made to develop efficient and safe readthrough-inducing compounds, which have been challenged in several models of human disease to provide with a therapy. In this view, the dissection of the molecular determinants shaping the outcome of readthrough, namely nucleotide and protein contexts as well as their interplay and impact on protein structure/function, is crucial to identify responsive nonsense mutations resulting in functional full-length proteins. The interpretation of experimental and mechanistic findings is also important to define a possibly clear picture of potential readthrough-favorable features useful to achieve rescue profiles compatible with therapeutic thresholds typical of each targeted disorder, which is of primary importance for the potential translatability of readthrough into a personalized and mutation-specific, and thus patient-oriented, therapeutic strategy.

2-Selenouridine, a Modified Nucleoside of Bacterial tRNAs, Its Reactivity in the Presence of Oxidizing and Reducing Reagents

The 5-substituted 2-selenouridines are natural components of the bacterial tRNA epitranscriptome. Because selenium-containing biomolecules are redox-active entities, the oxidation susceptibility of 2-selenouridine (Se2U) was studied in the presence of hydrogen peroxide under various conditions and compared with previously reported data for 2-thiouridine (S2U). It was found that Se2U is more susceptible to oxidation and converted in the first step to the corresponding diselenide (Se2U)₂, an unstable intermediate that decomposes to uridine and selenium. The reversibility of the oxidized state of Se2U was demonstrated by the efficient reduction of (Se2U)₂ to Se2U in the presence of common reducing agents. Thus, the 2-selenouridine component of tRNA may have antioxidant potential in cells because of its ability to react with both cellular ROS components and reducing agents. Interestingly, in the course of the reactions studied, we found that (Se2U)₂ reacts with Se2U to form new ‘oligomeric nucleosides′ as linear and cyclic byproducts.

mRNA and tRNA modification states influence ribosome speed and frame maintenance during poly(lysine) peptide synthesis

Ribosome speed is dictated by multiple factors including substrate availability, cellular conditions, and product (peptide) formation. Translation slows during the synthesis of cationic peptide sequences, potentially influencing the expression of thousands of proteins. Available evidence suggests that ionic interactions between positively charged nascent peptides and the negatively charged ribosome exit tunnel impede translation. However, this hypothesis was difficult to test directly because of inability to decouple the contributions of amino acid charge from mRNA sequence and tRNA identity/abundance in cells. Furthermore, it is unclear if other components of the translation system central to ribosome function (e.g., RNA modification) influence the speed and accuracy of positively charged peptide synthesis. In this study, we used a fully reconstituted Escherichia coli translation system to evaluate the effects of peptide charge, mRNA sequence, and RNA modification status on the translation of lysine-rich peptides. Comparison of translation reactions on poly(lysine)-encoding mRNAs conducted with either Lys-tRNA^Lys or Val-tRNA^Lys reveals that that amino acid charge, while important, only partially accounts for slowed translation on these transcripts. We further find that in addition to peptide charge, mRNA sequence and both tRNA and mRNA modification status influence the rates of amino acid addition and the ribosome’s ability to maintain frame (instead of entering the −2, −1, and +1 frames) during poly(lysine) peptide synthesis. Our observations lead us to expand the model for explaining how the ribosome slows during poly(lysine) peptide synthesis and suggest that posttranscriptional RNA modifications can provide cells a mechanism to precisely control ribosome movements along an mRNA.

Structural basis for PoxtA-mediated resistance to phenicol and oxazolidinone antibiotics

PoxtA and OptrA are ATP binding cassette (ABC) proteins of the F subtype (ABCF). They confer resistance to oxazolidinone and phenicol antibiotics, such as linezolid and chloramphenicol, which stall translating ribosomes when certain amino acids are present at a defined position in the nascent polypeptide chain. These proteins are often encoded on mobile genetic elements, facilitating their rapid spread amongst Gram-positive bacteria, and are thought to confer resistance by binding to the ribosome and dislodging the bound antibiotic. However, the mechanistic basis of this resistance remains unclear. Here we refine the PoxtA spectrum of action, demonstrate alleviation of linezolid-induced context-dependent translational stalling, and present cryo-electron microscopy structures of PoxtA in complex with the Enterococcus faecalis 70S ribosome. PoxtA perturbs the CCA-end of the P-site tRNA, causing it to shift by ∼4 Å out of the ribosome, corresponding to a register shift of approximately one amino acid for an attached nascent polypeptide chain. We postulate that the perturbation of the P-site tRNA by PoxtA thereby alters the conformation of the attached nascent chain to disrupt the drug binding site.

PoxtA confers resistance to ribosome-targeting oxazolidinone (linezolid) and chloramphenicol antibiotics. Here, Crowe-McAuliffe et al. provide structural insights into how binding of PoxtA to the ribosome indirectly promotes drug dissociation.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Synthesis and structure elucidation of the human tRNA nucleoside mannosyl-queuosine

Queuosine (Q) is a structurally complex, non-canonical RNA nucleoside. It is present in many eukaryotic and bacterial species, where it is part of the anticodon loop of certain tRNAs. In higher vertebrates, including humans, two further modified queuosine-derivatives exist - galactosyl- (galQ) and mannosyl-queuosine (manQ). The function of these low abundant hypermodified RNA nucleosides remains unknown. While the structure of galQ was elucidated and confirmed by total synthesis, the reported structure of manQ still awaits confirmation. By combining total synthesis and LC-MS-co-injection experiments, together with a metabolic feeding study of labelled hexoses, we show here that the natural compound manQ isolated from mouse liver deviates from the literature-reported structure. Our data show that manQ features an α-allyl connectivity of its sugar moiety. The yet unidentified glycosylases that attach galactose and mannose to the Q-base therefore have a maximally different constitutional connectivity preference. Knowing the correct structure of manQ will now pave the way towards further elucidation of its biological function.

Stabilization of Ribosomal RNA of the Small Subunit by Spermidine in Staphylococcus aureus

Cryo-electron microscopy is now used as a method of choice in structural biology for studying protein synthesis, a process mediated by the ribosome machinery. In order to achieve high-resolution structures using this approach, one needs to obtain homogeneous and stable samples, which requires optimization of ribosome purification in a species-dependent manner. This is especially critical for the bacterial small ribosomal subunit that tends to be unstable in the absence of ligands. Here, we report a protocol for purification of stable 30 S from the Gram-positive bacterium Staphylococcus aureus and its cryo-EM structures: in presence of spermidine at a resolution ranging between 3.4 and 3.6 Å and in its absence at 5.3 Å. Using biochemical characterization and cryo-EM, we demonstrate the importance of spermidine for stabilization of the 30 S via preserving favorable conformation of the helix 44.

Crystal structure of human METTL6, the m3C methyltransferase

In tRNA, the epigenetic m³C modification at position 32 in the anticodon loop is highly conserved in eukaryotes, which maintains the folding and basepairing functions of the anticodon. However, the responsible enzymes METTL2 and METTL6 were identified only in recent years. The loss of human METTL6 (hMETTL6) affects the translational process and proteostasis in cells, while in mESCs cells, it leads to defective pluripotency potential. Despite its important functions, the catalytic mechanism of the C32 methylation by this enzyme is poorly understood. Here we present the 1.9 Å high-resolution crystal structure of hMETTL6 bound by SAH. The key residues interacting with the ligand were identified and their roles were confirmed by ITC. We generated a docking model for the hMETTL6-SAH-CMP ternary complex. Interestingly, the CMP molecule binds into a cavity in a positive patch with the base ring pointing to the inside, suggesting a flipped-base mechanism for methylation. We further generated a model for the quaternary complex with tRNA^Ser as a component, which reasonably explained the biochemical behaviors of hMETTL6. Taken together, our crystallographic and biochemical studies provide important insight into the molecular recognition mechanism by METTL6 and may aid in the METTL-based rational drug design in the future.

The uridine to pseudouridine modification at the wobble position of eukaryotic isoleucine tRNA species is unlikely to induce mistranslation

Replacement of a U in an RNA duplex with a pseudouridine (Ψ), in general, stabilize the duplex because of the stronger stacking interaction, even concerning the wobble pair with G. The tRNA species specific to the AUA isoleucine codon in many eukaryotes have a Ψ at the first position of the anticodon. This tRNA^Ile would cause mistranslation if it could recognize the AUG codon through formation of a Ψ-G base pair. Here, I propose rationales for the minimal promotive effect of the U to Ψ modification on the mistranslation of the AUG codon.

THUMPD3-TRMT112 is a m2G methyltransferase working on a broad range of tRNA substrates

Post-transcriptional modifications affect tRNA biology and are closely associated with human diseases. However, progress on the functional analysis of tRNA modifications in metazoans has been slow because of the difficulty in identifying modifying enzymes. For example, the biogenesis and function of the prevalent N2-methylguanosine (m2G) at the sixth position of tRNAs in eukaryotes has long remained enigmatic. Herein, using a reverse genetics approach coupled with RNA-mass spectrometry, we identified that THUMP domain-containing protein 3 (THUMPD3) is responsible for tRNA: m2G6 formation in human cells. However, THUMPD3 alone could not modify tRNAs. Instead, multifunctional methyltransferase subunit TRM112-like protein (TRMT112) interacts with THUMPD3 to activate its methyltransferase activity. In the in vitro enzymatic assay system, THUMPD3-TRMT112 could methylate all the 26 tested G6-containing human cytoplasmic tRNAs by recognizing the characteristic 3'-CCA of mature tRNAs. We also showed that m2G7 of tRNATrp was introduced by THUMPD3-TRMT112. Furthermore, THUMPD3 is widely expressed in mouse tissues, with an extremely high level in the testis. THUMPD3-knockout cells exhibited impaired global protein synthesis and reduced growth. Our data highlight the significance of the tRNA: m2G6/7 modification and pave a way for further studies of the role of m2G in sperm tRNA derived fragments.

Regulation of gene expression via translational buffering

Translation of an mRNA represents a critical step during the expression of protein-coding genes. As mechanisms governing post-transcriptional regulation of gene expression are progressively unveiled, it is becoming apparent that transcriptional programs are not fully reflected in the proteome. Herein, we highlight a previously underappreciated post-transcriptional mode of regulation of gene expression termed translational buffering. In principle, translational buffering opposes the impact of alterations in mRNA levels on the proteome. We further describe three types of translational buffering: compensation, which maintains protein levels e.g. across species or individuals; equilibration, which retains pathway stoichiometry; and offsetting, which acts as a reversible mechanism that maintains the levels of selected subsets of proteins constant despite genetic alteration and/or stress-induced changes in corresponding mRNA levels. While mechanisms underlying compensation and equilibration have been reviewed elsewhere, the principal focus of this review is on the less-well understood mechanism of translational offsetting. Finally, we discuss potential roles of translational buffering in homeostasis and disease.

Regulatory roles of RNA modifications in breast cancer

Collectively referred to as the epitranscriptome, RNA modifications play important roles in gene expression control regulating relevant cellular processes. In the last few decades, growing numbers of RNA modifications have been identified not only in abundant ribosomal (rRNA) and transfer RNA (tRNA) but also in messenger RNA (mRNA). In addition, many writers, erasers and readers that dynamically regulate the chemical marks have also been characterized. Correct deposition of RNA modifications is prerequisite for cellular homeostasis, and its alteration results in aberrant transcriptional programs that dictate human disease, including breast cancer, the most frequent female malignancy, and the leading cause of cancer-related death in women. In this review, we emphasize the major RNA modifications that are present in tRNA, rRNA and mRNA. We have categorized breast cancer-associated chemical marks and summarize their contribution to breast tumorigenesis. In addition, we describe less abundant tRNA modifications with related pathways implicated in breast cancer. Finally, we discuss current limitations and perspectives on epitranscriptomics for use in therapeutic strategies against breast and other cancers.

A Path to the Atomic-Resolution Structures of Prokaryotic and Eukaryotic Ribosomes

Resolving first crystal structures of prokaryotic and eukaryotic ribosomes by our group has been based on the knowledge accumulated over the decades of studies, starting with the first electron microscopy images of the ribosome obtained by J. Pallade in 1955. In 1983, A. Spirin, then a Director of the Protein Research Institute of the USSR Academy of Sciences, initiated the first study aimed at solving the structure of ribosomes using X-ray structural analysis. In 1999, our group in collaboration with H. Noller published the first crystal structure of entire bacterial ribosome in a complex with its major functional ligands, such as messenger RNA and three transport RNAs at the A, P, and E sites. In 2011, our laboratory published the first atomic-resolution structure of eukaryotic ribosome solved by the X-ray diffraction analysis that confirmed the conserved nature of the main ribosomal functional components, such as the decoding and peptidyl transferase centers, was confirmed, and eukaryote-specific elements of the ribosome were described. Using X-ray structural analysis, we investigated general principles of protein biosynthesis inhibition in eukaryotic ribosomes, along with the mechanisms of antibiotic resistance. Structural differences between bacterial and eukaryotic ribosomes that determine the differences in their inhibition were established. These and subsequent atomic-resolution structures of the functional ribosome demonstrated for the first time the details of binding of messenger and transport RNAs, which was the first step towards understanding how the ribosome structure ultimately determines its functions.

Graph Theoretical Methods and Workflows for Searching and Annotation of RNA Tertiary Base Motifs and Substructures

The increasing number and complexity of structures containing RNA chains in the Protein Data Bank (PDB) have led to the need for automated structure annotation methods to replace or complement expert visual curation. This is especially true when searching for tertiary base motifs and substructures. Such base arrangements and motifs have diverse roles that range from contributions to structural stability to more direct involvement in the molecule's functions, such as the sites for ligand binding and catalytic activity. We review the utility of computational approaches in annotating RNA tertiary base motifs in a dataset of PDB structures, particularly the use of graph theoretical algorithms that can search for such base motifs and annotate them or find and annotate clusters of hydrogen-bond-connected bases. We also demonstrate how such graph theoretical algorithms can be integrated into a workflow that allows for functional analysis and comparisons of base arrangements and sub-structures, such as those involved in ligand binding. The capacity to carry out such automatic curations has led to the discovery of novel motifs and can give new context to known motifs as well as enable the rapid compilation of RNA 3D motifs into a database.

Whole-Genome Transformation Promotes tRNA Anticodon Suppressor Mutations under Stress

tRNAs are encoded by a large gene family, usually with several isogenic tRNAs interacting with the same codon. Mutations in the anticodon region of other tRNAs can overcome specific tRNA deficiencies. Phylogenetic analysis suggests that such mutations have occurred in evolution, but the driving force is unclear. We show that in yeast suppressor mutations in other tRNAs are able to overcome deficiency of the essential TRT2-encoded tRNAThrCGU at high temperature (40°C). Surprisingly, these tRNA suppressor mutations were obtained after whole-genome transformation with DNA from thermotolerant Kluyveromyces marxianus or Ogataea polymorpha strains but from which the mutations did apparently not originate. We suggest that transient presence of donor DNA in the host facilitates proliferation at high temperature and thus increases the chances for occurrence of spontaneous mutations suppressing defective growth at high temperature. Whole-genome sequence analysis of three transformants revealed only four to five nonsynonymous mutations of which one causing TRT2 anticodon stem stabilization and two anticodon mutations in non-threonyl-tRNAs, tRNALysCUU and tRNAeMetCAU, were causative. Both anticodon mutations suppressed lethality of TRT2 deletion and apparently caused the respective tRNAs to become novel substrates for threonyl-tRNA synthetase. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) data could not detect any significant mistranslation, and reverse transcription-quantitative PCR results contradicted induction of the unfolded protein response. We suggest that stress conditions have been a driving force in evolution for the selection of anticodon-switching mutations in tRNAs as revealed by phylogenetic analysis.IMPORTANCE In this work, we have identified for the first time the causative elements in a eukaryotic organism introduced by applying whole-genome transformation and responsible for the selectable trait of interest, i.e., high temperature tolerance. Surprisingly, the whole-genome transformants contained just a few single nucleotide polymorphisms (SNPs), which were unrelated to the sequence of the donor DNA. In each of three independent transformants, we have identified a SNP in a tRNA, either stabilizing the essential tRNAThrCGU at high temperature or switching the anticodon of tRNALysCUU or tRNAeMetCAU into CGU, which is apparently enough for in vivo recognition by threonyl-tRNA synthetase. LC-MS/MS analysis indeed indicated absence of significant mistranslation. Phylogenetic analysis showed that similar mutations have occurred throughout evolution and we suggest that stress conditions may have been a driving force for their selection. The low number of SNPs introduced by whole-genome transformation may favor its application for improvement of industrial yeast strains.

An RNA-centric historical narrative around the Protein Data Bank

Some of the amazing contributions brought to the scientific community by the Protein Data Bank (PDB) are described. The focus is on nucleic acid structures with a bias toward RNA. The evolution and key roles in science of the PDB and other structural databases for nucleic acids illustrate how small initial ideas can become huge and indispensable resources with the unflinching willingness of scientists to cooperate globally. The progress in the understanding of the molecular interactions driving RNA architectures followed the rapid increase in RNA structures in the PDB. That increase was consecutive to improvements in chemical synthesis and purification of RNA molecules, as well as in biophysical methods for structure determination and computer technology. The RNA modeling efforts from the early beginnings are also described together with their links to the state of structural knowledge and technological development. Structures of RNA and of its assemblies are physical objects, which, together with genomic data, allow us to integrate present-day biological functions and the historical evolution in all living species on earth.

On the Nature of Nucleobase Stacking in RNA: A Comprehensive Survey of Its Structural Variability and a Systematic Classification of Associated Interactions

The astonishing diversity in folding patterns of RNA three-dimensional (3D) structures is crafted by myriads of noncovalent contacts, of which base pairing and stacking are the most prominent. A systematic and comprehensive classification and annotation of these interactions is necessary for a molecular-level understanding of their roles. However, unlike in the case of base pairing, where a widely accepted nomenclature and classification scheme exists in the public domain, currently available classification schemes for base-base stacking need major enhancements to comprehensively capture the necessary features underlying the rich stacking diversity in RNA. Here, we extend the previous stacking classification based on nucleobase interacting faces by introducing a structurally intuitive geometry-cum topology-based scheme. Specifically, a stack is first classified in terms of the geometry described by the relative orientation of the glycosidic bonds, which generates eight basic stacking geometric families for heterodimeric stacks and six of those for homodimeric stacks. Further annotation in terms of the identity of the bases and the region of involvement of purines (five-membered, six-membered, or both rings) leads to the enumeration of 384 distinct RNA base stacks. Based on our classification scheme, we present an algorithm for automated identification of stacks in RNA crystal structures and analyze the stacking context in selected RNA structures. Overall, the work described here is expected to greatly facilitate the structure-based RNA research.

Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations

Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.

Structural basis of early translocation events on the ribosome

Peptide-chain elongation during protein synthesis entails sequential aminoacyl-tRNA selection and translocation reactions that proceed rapidly (2-20 per second) and with a low error rate (around 10^-3 to 10^-5 at each step) over thousands of cycles¹. The cadence and fidelity of ribosome transit through mRNA templates in discrete codon increments is a paradigm for movement in biological systems that must hold for diverse mRNA and tRNA substrates across domains of life. Here we use single-molecule fluorescence methods to guide the capture of structures of early translocation events on the bacterial ribosome. Our findings reveal that the bacterial GTPase elongation factor G specifically engages spontaneously achieved ribosome conformations while in an active, GTP-bound conformation to unlock and initiate peptidyl-tRNA translocation. These findings suggest that processes intrinsic to the pre-translocation ribosome complex can regulate the rate of protein synthesis, and that energy expenditure is used later in the translocation mechanism than previously proposed.

Terpene Chain Length Affects the Base Pairing Discrimination of S-geranyl-2-thiouridine in RNA Duplex

Geranylation is a hydrophobic modification discovered in several bacteria tRNAs with the function of promoting codon bias during translation. However, why nature selects this C10-geranyl group remains a question. We conduct synthesis, UV-thermal denaturation, and molecular simulation studies in RNA duplexes and reveal possible reasons behind this natural selection. Among methyl-(C1), dimethylallyl-(C5), geranyl-(C10), and farnesyl-(C15) modified 2-thiouridines, only geranyl-group promotes U:G over U:A pair. Molecular simulation shows all the modified terpene groups point to the minor groove of RNA duplexes. The discrimination between U:G and U:A pairs derives from the difference in hydrogen bonding and interactions of the chain with the hydrophobic area in the minor groove. Geranyl group has perfect length to discriminate U:G and U:A pairs, whereas the others are either too long or too short to achieve the same behavior. This work indicates that geranyl group cannot be replaced by other terpene groups in promoting codon-specificity.

•

Geranyl group on ges²U has perfect length to enhance binding to G over A in RNA

•

Other ges²U analogues has no pairing preference between G or A

•

MD shows the molecular dynamic behavior of different terpene chains in RNA duplex

Molecular Insights into Determinants of Translational Readthrough and Implications for Nonsense Suppression Approaches

The fidelity of protein synthesis, a process shaped by several mechanisms involving specialized ribosome regions and external factors, ensures the precise reading of sense and stop codons. However, premature termination codons (PTCs) arising from mutations may, at low frequency, be misrecognized and result in PTC suppression, named ribosome readthrough, with production of full-length proteins through the insertion of a subset of amino acids. Since some drugs have been identified as readthrough inducers, this fidelity drawback has been explored as a therapeutic approach in several models of human diseases caused by nonsense mutations. Here, we focus on the mechanisms driving translation in normal and aberrant conditions, the potential fates of mRNA in the presence of a PTC, as well as on the results obtained in the research of efficient readthrough-inducing compounds. In particular, we describe the molecular determinants shaping the outcome of readthrough, namely the nucleotide and protein context, with the latter being pivotal to produce functional full-length proteins. Through the interpretation of experimental and mechanistic findings, mainly obtained in lysosomal and coagulation disorders, we also propose a scenario of potential readthrough-favorable features to achieve relevant rescue profiles, representing the main issue for the potential translatability of readthrough as a therapeutic strategy.

Deflating the RNA Mg2+ bubble. Stereochemistry to the rescue!

Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac⁴C nucleotides in thermophilic archaeal tRNAs, the presence of m⁷G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

CUBAP: an interactive web portal for analyzing codon usage biases across populations

Synonymous codon usage significantly impacts translational and transcriptional efficiency, gene expression, the secondary structure of both mRNA and proteins, and has been implicated in various diseases. However, population-specific differences in codon usage biases remain largely unexplored. Here, we present a web server, https://cubap.byu.edu, to facilitate analyses of codon usage biases across populations (CUBAP). Using the 1000 Genomes Project, we calculated and visually depict population-specific differences in codon frequencies, codon aversion, identical codon pairing, co-tRNA codon pairing, ramp sequences, and nucleotide composition in 17,634 genes. We found that codon pairing significantly differs between populations in 35.8% of genes, allowing us to successfully predict the place of origin for African and East Asian individuals with 98.8% and 100% accuracy, respectively. We also used CUBAP to identify a significant bias toward decreased CTG pairing in the immunity related GTPase M (IRGM) gene in East Asian and African populations, which may contribute to the decreased association of rs10065172 with Crohn's disease in those populations. CUBAP facilitates in-depth gene-specific and codon-specific visualization that will aid in analyzing candidate genes identified in genome-wide association studies, identifying functional implications of synonymous variants, predicting population-specific impacts of synonymous variants and categorizing genetic biases unique to certain populations.

Eukaryotic life without tQCUG: the role of Elongator-dependent tRNA modifications in Dictyostelium discoideum

In the Elongator-dependent modification pathway, chemical modifications are introduced at the wobble uridines at position 34 in transfer RNAs (tRNAs), which serve to optimize codon translation rates. Here, we show that this three-step modification pathway exists in Dictyostelium discoideum, model of the evolutionary superfamily Amoebozoa. Not only are previously established modifications observable by mass spectrometry in strains with the most conserved genes of each step deleted, but also additional modifications are detected, indicating a certain plasticity of the pathway in the amoeba. Unlike described for yeast, D. discoideum allows for an unconditional deletion of the single tQ^CUG gene, as long as the Elongator-dependent modification pathway is intact. In gene deletion strains of the modification pathway, protein amounts are significantly reduced as shown by flow cytometry and Western blotting, using strains expressing different glutamine leader constructs fused to GFP. Most dramatic are these effects, when the tQ^CUG gene is deleted, or Elp3, the catalytic component of the Elongator complex is missing. In addition, Elp3 is the most strongly conserved protein of the modification pathway, as our phylogenetic analysis reveals. The implications of this observation are discussed with respect to the evolutionary age of the components acting in the Elongator-dependent modification pathway.