The wobble hypothesis revisited: uridine-5-oxyacetic acid is critical for reading of G-ending codons

tRNA modification profiling reveals epitranscriptome regulatory networks in Pseudomonas aeruginosa

Transfer RNA (tRNA) modifications have emerged as critical posttranscriptional regulators of gene expression affecting diverse biological and disease processes. While there is extensive knowledge about the enzymes installing the dozens of post-transcriptional tRNA modifications - the tRNA epitranscriptome - very little is known about how metabolic, signaling, and other networks integrate to regulate tRNA modification levels. Here we took a comprehensive first step at understanding epitranscriptome regulatory networks by developing a high-throughput tRNA isolation and mass spectrometry-based modification profiling platform and applying it to a Pseudomonas aeruginosa transposon insertion mutant library comprising 5,746 strains. Analysis of >200,000 tRNA modification data points validated the annotations of predicted tRNA modification genes, uncovered novel tRNA-modifying enzymes, and revealed tRNA modification regulatory networks in P. aeruginosa. Platform adaptation for RNA-seq library preparation would complement epitranscriptome studies, while application to human cell and mouse tissue demonstrates its utility for biomarker and drug discovery and 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.

Structures of the ribosome bound to EF-Tu-isoleucine tRNA elucidate the mechanism of AUG avoidance

The frequency of errors upon decoding of messenger RNA by the bacterial ribosome is low, with one misreading event per 1 × 104 codons. In the universal genetic code, the AUN codon box specifies two amino acids, isoleucine and methionine. In bacteria and archaea, decoding specificity of the AUA and AUG codons relies on the wobble avoidance strategy that requires modification of C34 in the anticodon loop of isoleucine transfer RNAIleCAU (tRNAIleCAU). Bacterial tRNAIleCAU with 2-lysylcytidine (lysidine) at the wobble position deciphers AUA while avoiding AUG. Here we report cryo-electron microscopy structures of the Escherichia coli 70S ribosome complexed with elongation factor thermo unstable (EF-Tu) and isoleucine-tRNAIleLAU in the process of decoding AUA and AUG. Lysidine in tRNAIleLAU excludes AUG by promoting the formation of an unusual Hoogsteen purine-pyrimidine nucleobase geometry at the third position of the codon, weakening the interactions with the mRNA and destabilizing the EF-Tu ternary complex. Our findings elucidate the molecular mechanism by which tRNAIleLAU specifically decodes AUA over AUG.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

A comprehensive analysis of translational misdecoding pattern and its implication on genetic code evolution

The universal genetic code is comprised of 61 sense codons, which are assigned to 20 canonical amino acids. However, the evolutionary basis for the highly conserved mapping between amino acids and their codons remains incompletely understood. A possible selective pressure of evolution would be minimization of deleterious effects caused by misdecoding. Here we comprehensively analyzed the misdecoding pattern of 61 codons against 19 noncognate amino acids where an arbitrary amino acid was omitted, and revealed the following two rules. (i) If the second codon base is U or C, misdecoding is frequently induced by mismatches at the first and/or third base, where any mismatches are widely tolerated; whereas misdecoding with the second-base mismatch is promoted by only U-G or C-A pair formation. (ii) If the second codon base is A or G, misdecoding is promoted by only G-U or U-G pair formation at the first or second position. In addition, evaluation of functional/structural diversities of amino acids revealed that less diverse amino acid sets are assigned at codons that induce more frequent misdecoding, and vice versa, so as to minimize deleterious effects of misdecoding in the modern genetic code.

Extensive breaking of genetic code degeneracy with non-canonical amino acids

Genetic code expansion (GCE) offers many exciting opportunities for the creation of synthetic organisms and for drug discovery methods that utilize in vitro translation. One type of GCE, sense codon reassignment (SCR), focuses on breaking the degeneracy of the 61 sense codons which encode for only 20 amino acids. SCR has great potential for genetic code expansion, but extensive SCR is limited by the post-transcriptional modifications on tRNAs and wobble reading of these tRNAs by the ribosome. To better understand codon-tRNA pairing, here we develop an assay to evaluate the ability of aminoacyl-tRNAs to compete with each other for a given codon. We then show that hyperaccurate ribosome mutants demonstrate reduced wobble reading, and when paired with unmodified tRNAs lead to extensive and predictable SCR. Together, we encode seven distinct amino acids across nine codons spanning just two codon boxes, thereby demonstrating that the genetic code hosts far more re-assignable space than previously expected, opening the door to extensive genetic code engineering.

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.

tRNA methylation resolves codon usage bias at the limit of cell viability

Codon usage of each genome is closely correlated with the abundance of tRNA isoacceptors. How codon usage bias is resolved by tRNA post-transcriptional modifications is largely unknown. Here we demonstrate that the N¹-methylation of guanosine at position 37 (m¹G37) on the 3′-side of the anticodon, while not directly responsible for reading of codons, is a neutralizer that resolves differential decoding of proline codons. A genome-wide suppressor screen of a non-viable Escherichia coli strain, lacking m¹G37, identifies proS suppressor mutations, indicating a coupling of methylation with tRNA prolyl-aminoacylation that sets the limit of cell viability. Using these suppressors, where prolyl-aminoacylation is decoupled from tRNA methylation, we show that m¹G37 neutralizes differential translation of proline codons by the major isoacceptor. Lack of m¹G37 inactivates this neutralization and exposes the need for a minor isoacceptor for cell viability. This work has medical implications for bacterial species that exclusively use the major isoacceptor for survival.

Masuda et al. show that loss of m¹G37 from the 3′ side of the tRNA anticodon renders a modified wobble nucleotide of the anticodon insufficient to decode a set of rare codons, providing a functional underpinning for the “modification circuit” between position 37 and the wobble position of the tRNA anticodon.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

“Superwobbling” and tRNA-34 Wobble and tRNA-37 Anticodon Loop Modifications in Evolution and Devolution of the Genetic Code

The genetic code evolved around the reading of the tRNA anticodon on the primitive ribosome, and tRNA-34 wobble and tRNA-37 modifications coevolved with the code. We posit that EF-Tu, the closing mechanism of the 30S ribosomal subunit, methylation of wobble U34 at the 5-carbon and suppression of wobbling at the tRNA-36 position were partly redundant and overlapping functions that coevolved to establish the code. The genetic code devolved in evolution of mitochondria to reduce the size of the tRNAome (all of the tRNAs of an organism or organelle). “Superwobbling” or four-way wobbling describes a major mechanism for shrinking the mitochondrial tRNAome. In superwobbling, unmodified wobble tRNA-U34 can recognize all four codon wobble bases (A, G, C and U), allowing a single unmodified tRNA-U34 to read a 4-codon box. During code evolution, to suppress superwobbling in 2-codon sectors, U34 modification by methylation at the 5-carbon position appears essential. As expected, at the base of code evolution, tRNA-37 modifications mostly related to the identity of the adjacent tRNA-36 base. TRNA-37 modifications help maintain the translation frame during elongation.

Implications of the Immune Polymorphisms of the Host and the Genetic Variability of SARS-CoV-2 in the Development of COVID-19

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for the current pandemic affecting almost all countries in the world. SARS-CoV-2 is the agent responsible for coronavirus disease 19 (COVID-19), which has claimed millions of lives around the world. In most patients, SARS-CoV-2 infection does not cause clinical signs. However, some infected people develop symptoms, which include loss of smell or taste, fever, dry cough, headache, severe pneumonia, as well as coagulation disorders. The aim of this work is to report genetic factors of SARS-CoV-2 and host-associated to severe COVID-19, placing special emphasis on the viral entry and molecules of the immune system involved with viral infection. Besides this, we analyze SARS-CoV-2 variants and their structural characteristics related to the binding to polymorphic angiotensin-converting enzyme type 2 (ACE2). Additionally, we also review other polymorphisms as well as some epigenetic factors involved in the immunopathogenesis of COVID-19. These factors and viral variability could explain the increment of infection rate and/or in the development of severe COVID-19.

Analysis and comparison of genetic variants and mutations of the novel coronavirus SARS-CoV-2

We present an analysis and comparison study of genetic variants and mutations of about 1200 genomes of SARS-CoV-2 virus sampled across the first seven months of 2020. The study includes 12 sets of about 100 genomes each collected between January and September. We analyzed the mutations, mutation frequency and count trends over time, and genomes trends over time from January through September. We show that certain mutations in the SARS-CoV-2 genome are not occurring randomly as it has been commonly believed. This finding is in agreement with other recently published research in this domain. Therefore, this validates other findings in this direction. This study includes approximately 1000 genomes and was able to identify over 35 different mutations most of which are common to almost all genomes groups. Some mutations' ratios (frequency percentage) fluctuate over time to adapt the virus to various environmental factors, climate, and populations. One of the interesting findings in this paper is that the coding region, at the nucleotide level for NSP13 protein is relatively conserved compared with other protein regions in the ORF1ab gene which makes this protein a good candidate for developing drug targets and treatment for the COVID-19 disease. Although this outcome was already reported by other researchers, we corroborated their result with our work in a different approach and another experimental setting with almost one thousand complete genome sequences. We presented and discussed all these results and findings with tables of results and illustrating figures.

Plastomes in the holoparasitic family Balanophoraceae: Extremely high AT content, severe gene content reduction, and two independent genetic code changes

The transition to a heterotrophic lifestyle in angiosperms is characterized by convergent evolutionary changes. Plastid genome remodeling includes dramatic functional and physical reductions with the highest degrees observed in fully heterotrophic plants. Genes related to photosynthesis are generally absent or pseudogenized, while a few genes related to other metabolic processes that take place within the plastid are almost invariably maintained. The family Balanophoraceae consists of root holoparasites that present reduced plastid genomes with an extraordinarily elevated AT content and the single genetic code change ever documented in land plant plastomes (the stop codon TAG now codes for tryptophan). Here, we studied the plastomes of Lophophytum leandri and Ombrophytum subterraneum (Balanophoraceae) that showed the remarkable absence of the gene trnE, a highly biased nucleotide composition, and an independent genetic code change (the standard stop codon TGA codes for tryptophan). This is the second genetic code change identified in land plant plastomes. Analysis of the transcriptome of Lophophytum indicated that the entire C5 pathway typical of plants is conserved despite the lack of trnE in its plastome. A hypothetical model of plastome evolution in the Balanophoraceae is presented.

Distinct evolutionary pathways for the synthesis and function of tRNA modifications

Transfer ribonucleicacids (RNAs) (tRNAs) are essential adaptor molecules for translation. The functions and stability of tRNAs are modulated by their post-transcriptional modifications (tRNA modifications). Each domain of life has a specific set of modifications that include ones shared in multiple domains and ones specific to a domain. In some cases, different tRNA modifications across domains have similar functions to each other. Recent studies uncovered that distinct enzymes synthesize the same modification in different organisms, suggesting that such modifications are acquired through independent evolution. In this short review, I outline the mechanisms by which various modifications contribute to tRNA function, including modulation of decoding and tRNA stability, using recent findings. I also focus on modifications that are synthesized by distinct biosynthetic pathways.

Construction of anti-codon table of the plant kingdom and evolution of tRNA selenocysteine (tRNASec)

Background

The tRNAs act as a bridge between the coding mRNA and incoming amino acids during protein translation. The anti-codon of tRNA recognizes the codon of the mRNA and deliver the amino acid into the protein translation chain. However, we did not know about the exact abundance of anti-codons in the genome and whether the frequency of abundance remains same across the plant lineage or not.

Results

Therefore, we analysed the tRNAnome of 128 plant species and reported an anti-codon table of the plant kingdom. We found that CAU anti-codon of tRNA^Met has highest (5.039%) whereas GCG anti-codon of tRNA^Arg has lowest (0.004%) abundance. However, when we compared the anti-codon frequencies according to the tRNA isotypes, we found tRNA^Leu (7.808%) has highest abundance followed by tRNA^Ser (7.668%) and tRNA^Gly (7.523%). Similarly, suppressor tRNA (0.036%) has lowest abundance followed by tRNA^Sec (0.066%) and tRNA^His (2.109). The genome of Ipomoea nil, Papaver somniferum, and Zea mays encoded the highest number of anti-codons (isoacceptor) at 59 each whereas the genome of Ostreococcus tauri was found to encode only 18 isoacceptors. The tRNA^Sec genes undergone losses more frequently than duplication and we found that tRNA^Sec showed anti-codon switch during the course of evolution.

Conclusion

The anti-codon table of the plant tRNA will enable us to understand the synonymous codon usage of the plant kingdom and can be very helpful to understand which codon is preferred over other during the translation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07216-3.

Identification and Characterization of Genes Required for 5-Hydroxyuridine Synthesis in Bacillus subtilis and Escherichia coli tRNA

In bacteria, tRNAs that decode 4-fold degenerate family codons and have uridine at position 34 of the anticodon are typically modified with either 5-methoxyuridine (mo⁵U) or 5-methoxycarbonylmethoxyuridine (mcmo⁵U). These modifications are critical for extended recognition of some codons at the wobble position. Whereas the alkylation steps of these modifications have been described, genes required for the hydroxylation of U34 to give 5-hydroxyuridine (ho⁵U) remain unknown. Here, a number of genes in Escherichia coli and Bacillus subtilis are identified that are required for wild-type (wt) levels of ho⁵U. The yrrMNO operon is identified in B. subtilis as important for the biosynthesis of ho⁵U. Both yrrN and yrrO are homologs to peptidase U32 family genes, which includes the rlhA gene required for ho⁵C synthesis in E. coli Deletion of either yrrN or yrrO, or both, gives a 50% reduction in mo⁵U tRNA levels. In E. coli, yegQ was found to be the only one of four peptidase U32 genes involved in ho⁵U synthesis. Interestingly, this mutant shows the same 50% reduction in (m)cmo⁵U as that observed for mo⁵U in the B. subtilis mutants. By analyzing the genomic context of yegQ homologs, the ferredoxin YfhL is shown to be required for ho⁵U synthesis in E. coli to the same extent as yegQ Additional genes required for Fe-S biosynthesis and biosynthesis of prephenate give the same 50% reduction in modification. Together, these data suggest that ho⁵U biosynthesis in bacteria is similar to that of ho⁵C, but additional genes and substrates are required for complete modification.IMPORTANCE Modified nucleotides in tRNA serve to optimize both its structure and function for accurate translation of the genetic code. The biosynthesis of these modifications has been fertile ground for uncovering unique biochemistry and metabolism in cells. In this work, genes that are required for a novel anaerobic hydroxylation of uridine at the wobble position of some tRNAs are identified in both Bacillus subtilis and Escherichia coli These genes code for Fe-S cluster proteins, and their deletion reduces the levels of the hydroxyuridine by 50% in both organisms. Additional genes required for Fe-S cluster and prephenate biosynthesis and a previously described ferredoxin gene all display a similar reduction in hydroxyuridine levels, suggesting that still other genes are required for the modification.

Bacterial wobble modifications of NNA-decoding tRNAs

Nucleotides of transfer RNAs (tRNAs) are highly modified, particularly at the anticodon. Bacterial tRNAs that read A-ending codons are especially notable. The U34 nucleotide canonically present in these tRNAs is modified by a wide range of complex chemical constituents. An additional two A-ending codons are not read by U34-containing tRNAs but are accommodated by either inosine or lysidine at the wobble position (I34 or L34). The structural basis for many N34 modifications in both tRNA aminoacylation and ribosome decoding has been elucidated, and evolutionary conservation of modifying enzymes is also becoming clearer. Here we present a brief review of the structure, function, and conservation of wobble modifications in tRNAs that translate A-ending codons. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1158-1166, 2019.

Dual pathways of tRNA hydroxylation ensure efficient translation by expanding decoding capability

In bacterial tRNAs, 5-carboxymethoxyuridine (cmo⁵U) and its derivatives at the first position of the anticodon facilitate non-Watson-Crick base pairing with guanosine and pyrimidines at the third positions of codons, thereby expanding decoding capabilities. However, their biogenesis and physiological roles remained to be investigated. Using reverse genetics and comparative genomics, we identify two factors responsible for 5-hydroxyuridine (ho⁵U) formation, which is the first step of the cmo⁵U synthesis: TrhP (formerly known as YegQ), a peptidase U32 family protein, is involved in prephenate-dependent ho⁵U formation; and TrhO (formerly known as YceA), a rhodanese family protein, catalyzes oxygen-dependent ho⁵U formation and bypasses cmo⁵U biogenesis in a subset of tRNAs under aerobic conditions. E. coli strains lacking both trhP and trhO exhibit a temperature-sensitive phenotype, and decode codons ending in G (GCG and UCG) less efficiently than the wild-type strain. These findings confirm that tRNA hydroxylation ensures efficient decoding during protein synthesis.

Wobbling Forth and Drifting Back: The Evolutionary History and Impact of Bacterial tRNA Modifications

Along with tRNAs, enzymes that modify anticodon bases are a key aspect of translation across the tree of life. tRNA modifications extend wobble pairing, allowing specific (“target”) tRNAs to recognize multiple codons and cover for other (“nontarget”) tRNAs, often improving translation efficiency and accuracy. However, the detailed evolutionary history and impact of tRNA modifying enzymes has not been analyzed. Using ancestral reconstruction of five tRNA modifications across 1093 bacteria, we show that most modifications were ancestral to eubacteria, but were repeatedly lost in many lineages. Most modification losses coincided with evolutionary shifts in nontarget tRNAs, often driven by increased bias in genomic GC and associated codon use, or by genome reduction. In turn, the loss of tRNA modifications stabilized otherwise highly dynamic tRNA gene repertoires. Our work thus traces the complex history of bacterial tRNA modifications, providing the first clear evidence for their role in the evolution of bacterial translation.

5-Oxyacetic Acid Modification Destabilizes Double Helical Stem Structures and Favors Anionic Watson-Crick like cmo5 U-G Base Pairs

Watson-Crick like G-U mismatches with tautomeric Genol or Uenol bases can evade fidelity checkpoints and thereby contribute to translational errors. The 5-oxyacetic acid uridine (cmo5 U) modification is a base modification at the wobble position on tRNAs and is presumed to expand the decoding capability of tRNA at this position by forming Watson-Crick like cmo5 Uenol -G mismatches. A detailed investigation on the influence of the cmo5 U modification on structural and dynamic features of RNA was carried out by using solution NMR spectroscopy and UV melting curve analysis. The introduction of a stable isotope labeled variant of the cmo5 U modifier allowed the application of relaxation dispersion NMR to probe the potentially formed Watson-Crick like cmo5 Uenol -G base pair. Surprisingly, we find that at neutral pH, the modification promotes transient formation of anionic Watson-Crick like cmo5 U- -G, and not enolic base pairs. Our results suggest that recoding is mediated by an anionic Watson-Crick like species, as well as bring an interesting aspect of naturally occurring RNA modifications into focus-the fine tuning of nucleobase properties leading to modulation of the RNA structural landscape by adoption of alternative base pairing patterns.

Transcriptional Regulation Buffers Gene Dosage Effects on a Highly Expressed Operon in Salmonella

A feature of bacterial chromosomes is that highly expressed essential genes are usually located close to the origin of replication. Because bacteria have overlapping cycles of replication, genes located close to the origin will often be present in multiple copies, and this is thought to be of selective benefit where high levels of expression support high growth rate. However, the magnitude of this selective effect and whether other forces could be at play are poorly understood. To study this, we translocated a highly expressed essential operon, tufB, to different locations and measured growth fitness. We found that transcriptional regulation buffered the effects of translocation and that even under conditions where growth rate was reduced, genetic changes that increased the expression of tufB were easily and rapidly selected. We conclude, at least for tufB, that forces other than gene dosage may be significant in selecting for chromosomal location.

Highly expressed genes are commonly located close to the origin of replication of bacterial chromosomes (OriC). This location skew is thought to reflect selective advantages associated with gene dosage effects during the replication cycle. The expression of constitutively expressed genes can vary up to fivefold based on chromosomal location, but it is not clear what level of variation would occur in naturally regulated operons. We tested the magnitude of the chromosome location effect using EF-Tu (tufA, tufB), an abundant protein whose cellular level correlates with, and limits, the maximum growth rate. We translocated the Salmonella tufB operon to four locations across the chromosome. The distance from OriC had only a small effect on growth rate, consistent with this operon having the natural ability to upregulate expression and compensate for reduced gene dosage. In contrast, when the total EF-Tu concentration was limiting for the growth rate (tufA deleted), we observed a strong gene dosage effect when tufB was located further from OriC. However, only a short period of experimental evolution was required before the bacteria adapted to this EF-Tu starvation situation by acquiring genetic changes that increased expression levels from the translocated tufB gene, restoring growth rates. Our findings demonstrate that, at least for the tufB operon, gene dosage is probably not the dominant force selecting for a chromosomal location close to OriC. We suggest that the colocation of highly expressed genes close to OriC might instead be selected because it enhances their coregulation during various growth states, with gene dosage being a secondary benefit.

Genetic code expansion via integration of redundant amino acid assignment by finely tuning tRNA pools

In all translation systems, the genetic code assigns codons to amino acids as building blocks of polypeptides, defining their chemical, structural and physiological properties. The canonical genetic code, however, utilizes only 20 proteinogenic amino acids redundantly encoded in 61 codons. In order to expand the building block repertoire, this redundancy was reduced by tuning composition of the transfer RNA (tRNA) mixture in vitro. Depletion of particular tRNAs from the total tRNA mixture or its reconstitution with in vitro-transcribed tRNA_SNNs (S = C or G, N = U, C, A or G) divided a codon box to encode two amino acids, expanding the repertoire to 23. The expanded genetic codes may benefit analysis of cellular regulatory pathways and drug screening.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Genome-Wide Analysis Reveals Ancestral Lack of Seventeen Different tRNAs and Clade-Specific Loss of tRNA-CNNs in Archaea

Transfer RNA (tRNA) is a category of RNAs that specifically decode messenger RNAs (mRNAs) into proteins by recognizing a set of 61 codons commonly adopted by different life domains. The composition and abundance of tRNAs play critical roles in shaping codon usage and pairing bias, which subsequently modulate mRNA translation efficiency and accuracy. Over the past few decades, effort has been concentrated on evaluating the specificity and redundancy of different tRNA families. However, the mechanism and processes underlying tRNA evolution have only rarely been investigated. In this study, by surveying tRNA genes in 167 completely sequenced genomes, we systematically investigated the composition and evolution of tRNAs in Archaea from a phylogenetic perspective. Our data revealed that archaeal genomes are compact in both tRNA types and copy number. Generally, no more than 44 different types of tRNA are present in archaeal genomes to decode the 61 canonical codons, and most of them have only one gene copy per genome. Among them, tRNA-Met was significantly overrepresented, with an average of three copies per genome. In contrast, the tRNA-UAU and 16 tRNAs with A-starting anticodons (tRNA-ANNs) were rarely detected in all archaeal genomes. The conspicuous absence of these tRNAs across the archaeal phylogeny suggests they might have not been evolved in the common ancestor of Archaea, rather than have lost independently from different clades. Furthermore, widespread absence of tRNA-CNNs in the Methanococcales and Methanobacteriales genomes indicates convergent loss of these tRNAs in the two clades. This clade-specific tRNA loss may be attributing to the reductive evolution of their genomes. Our data suggest that the current tRNA profiles in Archaea are contributed not only by the ancestral tRNA composition, but also by differential maintenance and loss of redundant tRNAs.

Genome recoding by tRNA modifications

RNA modifications are emerging as an additional regulatory layer on top of the primary RNA sequence. These modifications are particularly enriched in tRNAs where they can regulate not only global protein translation, but also protein translation at the codon level. Modifications located in or in the vicinity of tRNA anticodons are highly conserved in eukaryotes and have been identified as potential regulators of mRNA decoding. Recent studies have provided novel insights into how these modifications orchestrate the speed and fidelity of translation to ensure proper protein homeostasis. This review highlights the prominent modifications in the tRNA anticodon loop: queuosine, inosine, 5-methoxycarbonylmethyl-2-thiouridine, wybutosine, threonyl-carbamoyl-adenosine and 5-methylcytosine. We discuss the functional relevance of these modifications in protein translation and their emerging role in eukaryotic genome recoding during cellular adaptation and disease.

Celebrating wobble decoding: Half a century and still much is new

A simple post-transcriptional modification of tRNA, deamination of adenosine to inosine at the first, or wobble, position of the anticodon, inspired Francis Crick's Wobble Hypothesis 50 years ago. Many more naturally-occurring modifications have been elucidated and continue to be discovered. The post-transcriptional modifications of tRNA's anticodon domain are the most diverse and chemically complex of any RNA modifications. Their contribution with regards to chemistry, structure and dynamics reveal individual and combined effects on tRNA function in recognition of cognate and wobble codons. As forecast by the Modified Wobble Hypothesis 25 years ago, some individual modifications at tRNA's wobble position have evolved to restrict codon recognition whereas others expand the tRNA's ability to read as many as four synonymous codons. Here, we review tRNA wobble codon recognition using specific examples of simple and complex modification chemistries that alter tRNA function. Understanding natural modifications has inspired evolutionary insights and possible innovation in protein synthesis.

DAMBE6: New Tools for Microbial Genomics, Phylogenetics, and Molecular Evolution

DAMBE is a comprehensive software workbench for data analysis in molecular biology, phylogenetics, and evolution. Several important new functions have been added since version 5 of DAMBE: 1) comprehensive genomic profiling of translation initiation efficiency of different genes in different prokaryotic species, 2) a new index of translation elongation (I_TE) that takes into account both tRNA-mediated selection and background mutation on codon–anticodon adaptation, 3) a new and accurate phylogenetic approach based on pairwise alignment only, which is useful for highly divergent sequences from which a reliable multiple sequence alignment is difficult to obtain. Many other functions have been updated and improved including PWM for motif characterization, Gibbs sampler for de novo motif discovery, hidden Markov models for protein secondary structure prediction, self-organizing map for nonlinear clustering of transcriptomic data, comprehensive sequence alignment, and phylogenetic functions. DAMBE features a graphic, user-friendly and intuitive interface, and is freely available from http://dambe.bio.uottawa.ca.

The Importance of Being Modified: The Role of RNA Modifications in Translational Fidelity

The posttranscriptional modifications of tRNA's anticodon stem and loop (ASL) domain represent a third level, a third code, to the accuracy and efficiency of translating mRNA codons into the correct amino acid sequence of proteins. Modifications of tRNA's ASL domain are enzymatically synthesized and site specifically located at the anticodon wobble position-34 and 3'-adjacent to the anticodon at position-37. Degeneracy of the 64 Universal Genetic Codes and the limitation in the number of tRNA species require some tRNAs to decode more than one codon. The specific modification chemistries and their impact on the tRNA's ASL structure and dynamics enable one tRNA to decode cognate and "wobble codons" or to expand recognition to synonymous codons, all the while maintaining the translational reading frame. Some modified nucleosides' chemistries prestructure tRNA to read the two codons of a specific amino acid that shares a twofold degenerate codon box, and other chemistries allow a different tRNA to respond to all four codons of a fourfold degenerate codon box. Thus, tRNA ASL modifications are critical and mutations in genes for the modification enzymes and tRNA, the consequences of which is a lack of modification, lead to mistranslation and human disease. By optimizing tRNA anticodon chemistries, structure, and dynamics in all organisms, modifications ensure translational fidelity of mRNA transcripts.

An unmodified wobble uridine in tRNAs specific for Glutamine, Lysine, and Glutamic acid from Salmonella enterica Serovar Typhimurium results in nonviability-Due to increased missense errors?

In the wobble position of tRNAs specific for Gln, Lys, and Glu a universally conserved 5-methylene-2-thiouridine derivative (xm5s2U34, x denotes any of several chemical substituents and 34 denotes the wobble position) is present, which is 5-(carboxy)methylaminomethyl-2-thiouridine ((c)mnm5s2U34) in Bacteria and 5-methylcarboxymethyl-2-thiouridine (mcm5s2U34) in Eukarya. Here we show that mutants of the bacterium Salmonella enterica Serovar Typhimurium LT2 lacking either the s2- or the (c)mnm5-group of (c)mnm5s2U34 grow poorly especially at low temperature and do not grow at all at 15°C in both rich and glucose minimal media. A double mutant of S. enterica lacking both the s2- and the (c)mnm5-groups, and that thus has an unmodified uridine as wobble nucleoside, is nonviable at different temperatures. Overexpression of [Formula: see text] lacking either the s2- or the (c)mnm5-group and of [Formula: see text] lacking the s2-group exaggerated the reduced growth induced by the modification deficiency, whereas overexpression of [Formula: see text] lacking the mnm5-group did not. From these results we suggest that the primary function of cmnm5s2U34 in bacterial [Formula: see text] and mnm5s2U34 in [Formula: see text] is to prevent missense errors, but the mnm5-group of [Formula: see text] does not. However, other translational errors causing the growth defect cannot be excluded. These results are in contrast to what is found in yeast, since overexpression of the corresponding hypomodified yeast tRNAs instead counteracts the modification deficient induced phenotypes. Accordingly, it was suggested that the primary function of mcm5s2U34 in these yeast tRNAs is to improve cognate codon reading rather than prevents missense errors. Thus, although the xm5s2U34 derivatives are universally conserved, their major functional impact on bacterial and eukaryotic tRNAs may be different.

tRNA Modifications: Impact on Structure and Thermal Adaptation

Transfer RNAs (tRNAs) are central players in translation, functioning as adapter molecules between the informational level of nucleic acids and the functional level of proteins. They show a highly conserved secondary and tertiary structure and the highest density of post-transcriptional modifications among all RNAs. These modifications concentrate in two hotspots-the anticodon loop and the tRNA core region, where the D- and T-loop interact with each other, stabilizing the overall structure of the molecule. These modifications can cause large rearrangements as well as local fine-tuning in the 3D structure of a tRNA. The highly conserved tRNA shape is crucial for the interaction with a variety of proteins and other RNA molecules, but also needs a certain flexibility for a correct interplay. In this context, it was shown that tRNA modifications are important for temperature adaptation in thermophilic as well as psychrophilic organisms, as they modulate rigidity and flexibility of the transcripts, respectively. Here, we give an overview on the impact of modifications on tRNA structure and their importance in thermal adaptation.

tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence

Microbial pathogens adapt to the stress of infection by regulating transcription, translation and protein modification. We report that changes in gene expression in hypoxia-induced non-replicating persistence in mycobacteria—which models tuberculous granulomas—are partly determined by a mechanism of tRNA reprogramming and codon-biased translation. Mycobacterium bovis BCG responded to each stage of hypoxia and aerobic resuscitation by uniquely reprogramming 40 modified ribonucleosides in tRNA, which correlate with selective translation of mRNAs from families of codon-biased persistence genes. For example, early hypoxia increases wobble cmo⁵U in tRNA^Thr(UGU), which parallels translation of transcripts enriched in its cognate codon, ACG, including the DosR master regulator of hypoxic bacteriostasis. Codon re-engineering of dosR exaggerates hypoxia-induced changes in codon-biased DosR translation, with altered dosR expression revealing unanticipated effects on bacterial survival during hypoxia. These results reveal a coordinated system of tRNA modifications and translation of codon-biased transcripts that enhance expression of stress response proteins in mycobacteria.

Mycobacteria can adapt to the stress of human infection by entering a dormant state. Here the authors show that hypoxia-induced dormancy in M. bovis BCG involves the reprogramming of tRNA wobble modifications and copy numbers, coupled with biased use of synonymous codons in survival genes.

An integrated, structure- and energy-based view of the genetic code

The principles of mRNA decoding are conserved among all extant life forms. We present an integrative view of all the interaction networks between mRNA, tRNA and rRNA: the intrinsic stability of codon-anticodon duplex, the conformation of the anticodon hairpin, the presence of modified nucleotides, the occurrence of non-Watson-Crick pairs in the codon-anticodon helix and the interactions with bases of rRNA at the A-site decoding site. We derive a more information-rich, alternative representation of the genetic code, that is circular with an unsymmetrical distribution of codons leading to a clear segregation between GC-rich 4-codon boxes and AU-rich 2:2-codon and 3:1-codon boxes. All tRNA sequence variations can be visualized, within an internal structural and energy framework, for each organism, and each anticodon of the sense codons. The multiplicity and complexity of nucleotide modifications at positions 34 and 37 of the anticodon loop segregate meaningfully, and correlate well with the necessity to stabilize AU-rich codon-anticodon pairs and to avoid miscoding in split codon boxes. The evolution and expansion of the genetic code is viewed as being originally based on GC content with progressive introduction of A/U together with tRNA modifications. The representation we present should help the engineering of the genetic code to include non-natural amino acids.

Structural effects of modified ribonucleotides and magnesium in transfer RNAs

Modified nucleotides are ubiquitous and important to tRNA structure and function. To understand their effect on tRNA conformation, we performed a series of molecular dynamics simulations on yeast tRNA^Phe and tRNA_init, Escherichia coli tRNA_init and HIV tRNA^Lys. Simulations were performed with the wild type modified nucleotides, using the recently developed CHARMM compatible force field parameter set for modified nucleotides (J. Comput. Chem.2016, 37, 896), or with the corresponding unmodified nucleotides, and in the presence or absence of Mg²⁺. Results showed a stabilizing effect associated with the presence of the modifications and Mg²⁺ for some important positions, such as modified guanosine in position 37 and dihydrouridines in 16/17 including both structural properties and base interactions. Some other modifications were also found to make subtle contributions to the structural properties of local domains. While we were not able to investigate the effect of adenosine 37 in tRNA_init and limitations were observed in the conformation of E. coli tRNA_init, the presence of the modified nucleotides and of Mg²⁺ better maintained the structural features and base interactions of the tRNA systems than in their absence indicating the utility of incorporating the modified nucleotides in simulations of tRNA and other RNAs.

Biogenesis and growth phase-dependent alteration of 5-methoxycarbonylmethoxyuridine in tRNA anticodons

Post-transcriptional modifications at the anticodon first (wobble) position of tRNA play critical roles in precise decoding of genetic codes. 5-carboxymethoxyuridine (cmo⁵U) and its methyl ester derivative 5-methoxycarbonylmethoxyuridine (mcmo⁵U) are modified nucleosides found at the anticodon wobble position in several tRNAs from Gram-negative bacteria. cmo⁵U and mcmo⁵U facilitate non-Watson–Crick base pairing with guanosine and pyrimidines at the third positions of codons, thereby expanding decoding capabilities. By mass spectrometric analyses of individual tRNAs and a shotgun approach of total RNA from Escherichia coli, we identified mcmo⁵U as a major modification in tRNA^Ala1, tRNA^Ser1, tRNA^Pro3 and tRNA^Thr4; by contrast, cmo⁵U was present primarily in tRNA^Leu3 and tRNA^Val1. In addition, we discovered 5-methoxycarbonylmethoxy-2′-O-methyluridine (mcmo⁵Um) as a novel but minor modification in tRNA^Ser1. Terminal methylation frequency of mcmo⁵U in tRNA^Pro3 was low (≈30%) in the early log phase of cell growth, gradually increased as growth proceeded and reached nearly 100% in late log and stationary phases. We identified CmoM (previously known as SmtA), an AdoMet-dependent methyltransferase that methylates cmo⁵U to form mcmo⁵U. A luciferase reporter assay based on a +1 frameshift construct revealed that terminal methylation of mcmo⁵U contributes to the decoding ability of tRNA^Ala1.

Modification of the wobble uridine in bacterial and mitochondrial tRNAs reading NNA/NNG triplets of 2-codon boxes

Posttranscriptional modification of the uridine located at the wobble position (U34) of tRNAs is crucial for optimization of translation. Defects in the U34 modification of mitochondrial-tRNAs are associated with a group of rare diseases collectively characterized by the impairment of the oxidative phosphorylation system. Retrograde signaling pathways from mitochondria to nucleus are involved in the pathophysiology of these diseases. These pathways may be triggered by not only the disturbance of the mitochondrial (mt) translation caused by hypomodification of tRNAs, but also as a result of nonconventional roles of mt-tRNAs and mt-tRNA-modifying enzymes. The evolutionary conservation of these enzymes supports their importance for cell and organismal functions. Interestingly, bacterial and eukaryotic cells respond to stress by altering the expression or activity of these tRNA-modifying enzymes, which leads to changes in the modification status of tRNAs. This review summarizes recent findings about these enzymes and sets them within the previous data context.

Codon-Anticodon Recognition in the Bacillus subtilis glyQS T Box Riboswitch: RNA-DEPENDENT CODON SELECTION OUTSIDE THE RIBOSOME

Many amino acid-related genes in Gram-positive bacteria are regulated by the T box riboswitch. The leader RNA of genes in the T box family controls the expression of downstream genes by monitoring the aminoacylation status of the cognate tRNA. Previous studies identified a three-nucleotide codon, termed the "Specifier Sequence," in the riboswitch that corresponds to the amino acid identity of the downstream genes. Pairing of the Specifier Sequence with the anticodon of the cognate tRNA is the primary determinant of specific tRNA recognition. This interaction mimics codon-anticodon pairing in translation but occurs in the absence of the ribosome. The goal of the current study was to determine the effect of a full range of mismatches for comparison with codon recognition in translation. Mutations were individually introduced into the Specifier Sequence of the glyQS leader RNA and tRNA(Gly) anticodon to test the effect of all possible pairing combinations on tRNA binding affinity and antitermination efficiency. The functional role of the conserved purine 3' of the Specifier Sequence was also verifiedin this study. We found that substitutions at the Specifier Sequence resulted in reduced binding, the magnitude of which correlates well with the predicted stability of the RNA-RNA pairing. However, the tolerance for specific mismatches in antitermination was generally different from that during decoding, which reveals a unique tRNA recognition pattern in the T box antitermination system.

Determinants of the CmoB carboxymethyl transferase utilized for selective tRNA wobble modification

Enzyme-mediated modifications at the wobble position of tRNAs are essential for the translation of the genetic code. We report the genetic, biochemical and structural characterization of CmoB, the enzyme that recognizes the unique metabolite carboxy-S-adenosine-L-methionine (Cx-SAM) and catalyzes a carboxymethyl transfer reaction resulting in formation of 5-oxyacetyluridine at the wobble position of tRNAs. CmoB is distinctive in that it is the only known member of the SAM-dependent methyltransferase (SDMT) superfamily that utilizes a naturally occurring SAM analog as the alkyl donor to fulfill a biologically meaningful function. Biochemical and genetic studies define the in vitro and in vivo selectivity for Cx-SAM as alkyl donor over the vastly more abundant SAM. Complementary high-resolution structures of the apo- and Cx-SAM bound CmoB reveal the determinants responsible for this remarkable discrimination. Together, these studies provide mechanistic insight into the enzymatic and non-enzymatic feature of this alkyl transfer reaction which affords the broadened specificity required for tRNAs to recognize multiple synonymous codons.

Experimental strategies for functional annotation and metabolism discovery: targeted screening of solute binding proteins and unbiased panning of metabolomes

The rate at which genome sequencing data is accruing demands enhanced methods for functional annotation and metabolism discovery. Solute binding proteins (SBPs) facilitate the transport of the first reactant in a metabolic pathway, thereby constraining the regions of chemical space and the chemistries that must be considered for pathway reconstruction. We describe high-throughput protein production and differential scanning fluorimetry platforms, which enabled the screening of 158 SBPs against a 189 component library specifically tailored for this class of proteins. Like all screening efforts, this approach is limited by the practical constraints imposed by construction of the library, i.e., we can study only those metabolites that are known to exist and which can be made in sufficient quantities for experimentation. To move beyond these inherent limitations, we illustrate the promise of crystallographic- and mass spectrometric-based approaches for the unbiased use of entire metabolomes as screening libraries. Together, our approaches identified 40 new SBP ligands, generated experiment-based annotations for 2084 SBPs in 71 isofunctional clusters, and defined numerous metabolic pathways, including novel catabolic pathways for the utilization of ethanolamine as sole nitrogen source and the use of d-Ala-d-Ala as sole carbon source. These efforts begin to define an integrated strategy for realizing the full value of amassing genome sequence data.

A major controversy in codon-anticodon adaptation resolved by a new codon usage index

Two alternative hypotheses attribute different benefits to codon-anticodon adaptation. The first assumes that protein production is rate limited by both initiation and elongation and that codon-anticodon adaptation would result in higher elongation efficiency and more efficient and accurate protein production, especially for highly expressed genes. The second claims that protein production is rate limited only by initiation efficiency but that improved codon adaptation and, consequently, increased elongation efficiency have the benefit of increasing ribosomal availability for global translation. To test these hypotheses, a recent study engineered a synthetic library of 154 genes, all encoding the same protein but differing in degrees of codon adaptation, to quantify the effect of differential codon adaptation on protein production in Escherichia coli. The surprising conclusion that “codon bias did not correlate with gene expression” and that “translation initiation, not elongation, is rate-limiting for gene expression” contradicts the conclusion reached by many other empirical studies. In this paper, I resolve the contradiction by reanalyzing the data from the 154 sequences. I demonstrate that translation elongation accounts for about 17% of total variation in protein production and that the previous conclusion is due to the use of a codon adaptation index (CAI) that does not account for the mutation bias in characterizing codon adaptation. The effect of translation elongation becomes undetectable only when translation initiation is unrealistically slow. A new index of translation elongation I_TE is formulated to facilitate studies on the efficiency and evolution of the translation machinery.

Transfer RNA Modification: Presence, Synthesis, and Function

Transfer RNA (tRNA) from all organisms on this planet contains modified nucleosides, which are derivatives of the four major nucleosides. tRNA from Escherichia coli/Salmonella enterica serovar Typhimurium contains 33 different modified nucleosides, which are all, except one (Queuosine [Q]), synthesized on an oligonucleotide precursor, which by specific enzymes later matures into tRNA. The structural genes for these enzymes are found in mono- and polycistronic operons, the latter of which have a complex transcription and translation pattern. The synthesis of the tRNA-modifying enzymes is not regulated similarly, and it is not coordinated to that of their substrate, the tRNA. The synthesis of some of them (e.g., several methylated derivatives) is catalyzed by one enzyme, which is position and base specific, whereas synthesis of some has a very complex biosynthetic pathway involving several enzymes (e.g., 2-thiouridines, N  6-cyclicthreonyladenosine [ct6A], and Q). Several of the modified nucleosides are essential for viability (e.g., lysidin, ct6A, 1-methylguanosine), whereas the deficiency of others induces severe growth defects. However, some have no or only a small effect on growth at laboratory conditions. Modified nucleosides that are present in the anticodon loop or stem have a fundamental influence on the efficiency of charging the tRNA, reading cognate codons, and preventing missense and frameshift errors. Those that are present in the body of the tRNA primarily have a stabilizing effect on the tRNA. Thus, the ubiquitous presence of these modified nucleosides plays a pivotal role in the function of the tRNA by their influence on the stability and activity of the tRNA.

Recognition of Watson-Crick base pairs: constraints and limits due to geometric selection and tautomerism

The natural bases of nucleic acids have a strong preference for one tautomer form, guaranteeing fidelity in their hydrogen bonding potential. However, base pairs observed in recent crystal structures of polymerases and ribosomes are best explained by an alternative base tautomer, leading to the formation of base pairs with Watson-Crick-like geometries. These observations set limits to geometric selection in molecular recognition of complementary Watson-Crick pairs for fidelity in replication and translation processes.

The effect of mutation and selection on codon adaptation in Escherichia coli bacteriophage

Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.

Life without tRNAIle-lysidine synthetase: translation of the isoleucine codon AUA in Bacillus subtilis lacking the canonical tRNA2Ile

Translation of the isoleucine codon AUA in most prokaryotes requires a modified C (lysidine or agmatidine) at the wobble position of tRNA₂^Ile to base pair specifically with the A of the AUA codon but not with the G of AUG. Recently, a Bacillus subtilis strain was isolated in which the essential gene encoding tRNA^Ile-lysidine synthetase was deleted for the first time. In such a strain, C34 at the wobble position of tRNA₂^Ile is expected to remain unmodified and cells depend on a mutant suppressor tRNA derived from tRNA₁^Ile, in which G34 has been changed to U34. An important question, therefore, is how U34 base pairs with A without also base pairing with G. Here, we show (i) that unlike U34 at the wobble position of all B. subtilis tRNAs of known sequence, U34 in the mutant tRNA is not modified, and (ii) that the mutant tRNA binds strongly to the AUA codon on B. subtilis ribosomes but only weakly to AUG. These in vitro data explain why the suppressor strain displays only a low level of misreading AUG codons in vivo and, as shown here, grows at a rate comparable to that of the wild-type strain.

Changes in translational efficiency is a dominant regulatory mechanism in the environmental response of bacteria

To understand how cell physiological state affects mRNA translation, we used Shewanella oneidensis MR-1 grown under steady state conditions at either 20% or 8.5% O2. Using a combination of quantitative proteomics and RNA-Seq, we generated high-confidence data on >1000 mRNA and protein pairs. By using a steady state model, we found that differences in protein-mRNA ratios were primarily due to differences in the translational efficiency of specific genes. When oxygen levels were lowered, 28% of the proteins showed at least a 2-fold change in expression. Transcription levels were sp. significantly altered for 26% of the protein changes; translational efficiency was significantly altered for 46% and a combination of both was responsible for the remaining 28%. Changes in translational efficiency were significantly correlated with the codon usage pattern of the genes and measurable tRNA pools changed in response to altered O2 levels. Our results suggest that changes in the translational efficiency of proteins, in part due to altered tRNA pools, is a major determinant of regulated alterations in protein expression levels in bacteria.