New Structural Insights into Translational Miscoding

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

The central role of transfer RNAs in mistranslation

Transfer RNAs (tRNA) are essential small non-coding RNAs that enable the translation of genomic information into proteins in all life forms. The principal function of tRNAs is to bring amino acid building blocks to the ribosomes for protein synthesis. In the ribosome, tRNAs interact with messenger RNA (mRNA) to mediate the incorporation of amino acids into a growing polypeptide chain following the rules of the genetic code. Accurate interpretation of the genetic code requires tRNAs to carry amino acids matching their anticodon identity and decode the correct codon on mRNAs. Errors in these steps cause the translation of codons with the wrong amino acids (mistranslation), compromising the accurate flow of information from DNA to proteins. Accumulation of mutant proteins due to mistranslation jeopardizes proteostasis and cellular viability. However, the concept of mistranslation is evolving, with increasing evidence indicating that mistranslation can be used as a mechanism for survival and acclimatization to environmental conditions. In this review, we discuss the central role of tRNAs in modulating translational fidelity through their dynamic and complex interplay with translation factors. We summarize recent discoveries of mistranslating tRNAs and describe the underlying molecular mechanisms and the specific conditions and environments that enable and promote mistranslation.

Tuning tRNAs for improved translation

Transfer RNAs have been extensively explored as the molecules that translate the genetic code into proteins. At this interface of genetics and biochemistry, tRNAs direct the efficiency of every major step of translation by interacting with a multitude of binding partners. However, due to the variability of tRNA sequences and the abundance of diverse post-transcriptional modifications, a guidebook linking tRNA sequences to specific translational outcomes has yet to be elucidated. Here, we review substantial efforts that have collectively uncovered tRNA engineering principles that can be used as a guide for the tuning of translation fidelity. These principles have allowed for the development of basic research, expansion of the genetic code with non-canonical amino acids, and tRNA therapeutics.

Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure. The infrastructure provides access to QM/MM models that may serve as the foundation for QM/MM-ΔMLP potentials, which supplement the semiempirical QM/MM model with a MLP correction trained to reproduce ab initio QM/MM energies and forces. Efficient optimization of minimum free energy pathways is enabled through a new surface-accelerated finite-temperature string method implemented in the FE-ToolKit package. Furthermore, we interfaced Sander with the i-PI software by implementing the socket communication protocol used in the i-PI client-server model. The new interface with i-PI allows for the treatment of nuclear quantum effects with semiempirical QM/MM-ΔMLP models. The modular interoperable software is demonstrated on proton transfer reactions in guanine-thymine mispairs in a B-form deoxyribonucleic acid helix. The current work represents a considerable advance in the development of modular software for performing free energy simulations of chemical reactions that are important in a wide range of applications.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

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.

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.

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA

Proton transfer across hydrogen bonds in DNA can produce non-canonical nucleobase dimers and is a possible source of single-point mutations when these forms mismatch under replication. Previous computational studies have revealed this process to be energetically feasible for the guanine-cytosine (GC) base pair, but the tautomeric product (G[Formula: see text]C[Formula: see text]) is short-lived. In this work we reveal, for the first time, the direct effect of the replisome enzymes on proton transfer, rectifying the shortcomings of existing models. Multi-scale quantum mechanical/molecular dynamics (QM/MM) simulations reveal the effect of the bacterial PcrA Helicase on the double proton transfer in the GC base pair. It is shown that the local protein environment drastically increases the activation and reaction energies for the double proton transfer, modifying the tautomeric equilibrium. We propose a regime in which the proton transfer is dominated by tunnelling, taking place instantaneously and without atomic rearrangement of the local environment. In this paradigm, we can reconcile the metastable nature of the tautomer and show that ensemble averaging methods obscure detail in the reaction profile. Our results highlight the importance of explicit environmental models and suggest that asparagine N624 serves a secondary function of reducing spontaneous mutations in PcrA Helicase.

The 3 31 Nucleotide Minihelix tRNA Evolution Theorem and the Origin of Life

There are no theorems (proven theories) in the biological sciences. We propose that the 3 31 nt minihelix tRNA evolution theorem be universally accepted as one. The 3 31 nt minihelix theorem completely describes the evolution of type I and type II tRNAs from ordered precursors (RNA repeats and inverted repeats). Despite the diversification of tRNAome sequences, statistical tests overwhelmingly support the theorem. Furthermore, the theorem relates the dominant pathway for the origin of life on Earth, specifically, how tRNAomes and the genetic code may have coevolved. Alternate models for tRNA evolution (i.e., 2 minihelix, convergent and accretion models) are falsified. In the context of the pre-life world, tRNA was a molecule that, via mutation, could modify anticodon sequences and teach itself to code. Based on the tRNA sequence, we relate the clearest history to date of the chemical evolution of life. From analysis of tRNA evolution, ribozyme-mediated RNA ligation was a primary driving force in the evolution of complexity during the pre-life-to-life transition. TRNA formed the core for the evolution of living systems on Earth.

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.

Translation dynamics in human cells visualized at high resolution reveal cancer drug action

Ribosomes catalyze protein synthesis by cycling through various functional states. These states have been extensively characterized in vitro, but their distribution in actively translating human cells remains elusive. We used a cryo-electron tomography-based approach and resolved ribosome structures inside human cells with high resolution. These structures revealed the distribution of functional states of the elongation cycle, a Z transfer RNA binding site, and the dynamics of ribosome expansion segments. Ribosome structures from cells treated with Homoharringtonine, a drug used against chronic myeloid leukemia, revealed how translation dynamics were altered in situ and resolve the small molecules within the active site of the ribosome. Thus, structural dynamics and drug effects can be assessed at high resolution within human cells.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

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.

The life and times of a tRNA

The study of eukaryotic tRNA processing has given rise to an explosion of new information and insights in the last several years. We now have unprecedented knowledge of each step in the tRNA processing pathway, revealing unexpected twists in biochemical pathways, multiple new connections with regulatory pathways, and numerous biological effects of defects in processing steps that have profound consequences throughout eukaryotes, leading to growth phenotypes in the yeast Saccharomyces cerevisiae and to neurological and other disorders in humans. This review highlights seminal new results within the pathways that comprise the life of a tRNA, from its birth after transcription until its death by decay. We focus on new findings and revelations in each step of the pathway including the end-processing and splicing steps, many of the numerous modifications throughout the main body and anticodon loop of tRNA that are so crucial for tRNA function, the intricate tRNA trafficking pathways, and the quality control decay pathways, as well as the biogenesis and biology of tRNA-derived fragments. We also describe the many interactions of these pathways with signaling and other pathways in the cell.

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.

“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.

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.

RNA Modifications in Pathogenic Bacteria: Impact on Host Adaptation and Virulence

RNA modifications are involved in numerous biological processes and are present in all RNA classes. These modifications can be constitutive or modulated in response to adaptive processes. RNA modifications play multiple functions since they can impact RNA base-pairings, recognition by proteins, decoding, as well as RNA structure and stability. However, their roles in stress, environmental adaptation and during infections caused by pathogenic bacteria have just started to be appreciated. With the development of modern technologies in mass spectrometry and deep sequencing, recent examples of modifications regulating host-pathogen interactions have been demonstrated. They show how RNA modifications can regulate immune responses, antibiotic resistance, expression of virulence genes, and bacterial persistence. Here, we illustrate some of these findings, and highlight the strategies used to characterize RNA modifications, and their potential for new therapeutic applications.

Quantum and classical effects in DNA point mutations: Watson-Crick tautomerism in AT and GC base pairs

Proton transfer along the hydrogen bonds of DNA can lead to the creation of short-lived, but biologically relevant point mutations that can further lead to gene mutation and, potentially, cancer. In this work, the energy landscape of the canonical A-T and G-C base pairs (standard, amino-keto) to tautomeric A*-T* and G*-C* (non-standard, imino-enol) Watson-Crick DNA base pairs is modelled with density functional theory and machine-learning nudge-elastic band methods. We calculate the energy barriers and tunnelling rates of hydrogen transfer between and within each base monomer (A, T, G and C). We show that the role of tunnelling in A-T tautomerisation is statistically unlikely due to the presence of a small reverse reaction barrier. On the contrary, the thermal populations of the G*-C* point mutation could be non-trivial and propagate through the replisome. For the direct intramolecular transfer, the reaction is hindered by a substantial energy barrier. However, our calculations indicate that tautomeric bases in their monomeric form have remarkably long lifetimes.

Evolution of the genetic code

Diverse models have been advanced for the evolution of the genetic code. Here, models for tRNA, aminoacyl-tRNA synthetase (aaRS) and genetic code evolution were combined with an understanding of EF-Tu suppression of tRNA 3rd anticodon position wobbling. The result is a highly detailed scheme that describes the placements of all amino acids in the standard genetic code. The model describes evolution of 6-, 4-, 3-, 2- and 1-codon sectors. Innovation in column 3 of the code is explained. Wobbling and code degeneracy are explained. Separate distribution of serine sectors between columns 2 and 4 of the code is described. We conclude that very little chaos contributed to evolution of the genetic code and that the pattern of evolution of aaRS enzymes describes a history of the evolution of the code. A model is proposed to describe the biological selection for the earliest evolution of the code and for protocell evolution.

Trends of mutation accumulation across global SARS-CoV-2 genomes: Implications for the evolution of the novel coronavirus

To understand SARS-CoV-2 microevolution, this study explored the genome-wide frequency, gene-wise distribution, and molecular nature of all point-mutations detected across its 71,703 RNA-genomes deposited in GISAID till 21 August 2020. Globally, nsp1/nsp2 and orf7a/orf3a were the most mutation-ridden non-structural and structural genes respectively. Phylogeny of 4618 spatiotemporally-representative genomes revealed that entities belonging to the early lineages are mostly spread over Asian countries, including India, whereas the recently-derived lineages are more globally distributed. Of the total 20,163 instances of polymorphism detected across global genomes, 12,594 and 7569 involved transitions and transversions, predominated by cytidine-to-uridine and guanosine-to-uridine conversions, respectively. Positive selection of nonsynonymous mutations (dN/dS >1) in most of the structural, but not the non-structural, genes indicated that SARS-CoV-2 has already harmonized its replication/transcription machineries with the host metabolism, while it is still redefining virulence/transmissibility strategies at the molecular level. Mechanistic bases and evolutionary/pathogenicity-related implications are discussed for the predominant mutation-types.

The nature of the purine at position 34 in tRNAs of 4-codon boxes is correlated with nucleotides at positions 32 and 38 to maintain decoding fidelity

Translation fidelity relies essentially on the ability of ribosomes to accurately recognize triplet interactions between codons on mRNAs and anticodons of tRNAs. To determine the codon-anticodon pairs that are efficiently accepted by the eukaryotic ribosome, we took advantage of the IRES from the intergenic region (IGR) of the Cricket Paralysis Virus. It contains an essential pseudoknot PKI that structurally and functionally mimics a codon-anticodon helix. We screened the entire set of 4096 possible combinations using ultrahigh-throughput screenings combining coupled transcription/translation and droplet-based microfluidics. Only 97 combinations are efficiently accepted and accommodated for translocation and further elongation: 38 combinations involve cognate recognition with Watson-Crick pairs and 59 involve near-cognate recognition pairs with at least one mismatch. More than half of the near-cognate combinations (36/59) contain a G at the first position of the anticodon (numbered 34 of tRNA). G34-containing tRNAs decoding 4-codon boxes are almost absent from eukaryotic genomes in contrast to bacterial genomes. We reconstructed these missing tRNAs and could demonstrate that these tRNAs are toxic to cells due to their miscoding capacity in eukaryotic translation systems. We also show that the nature of the purine at position 34 is correlated with the nucleotides present at 32 and 38.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Evolution of Life on Earth: tRNA, Aminoacyl-tRNA Synthetases and the Genetic Code

Life on Earth and the genetic code evolved around tRNA and the tRNA anticodon. We posit that the genetic code initially evolved to synthesize polyglycine as a cross-linking agent to stabilize protocells. We posit that the initial amino acids to enter the code occupied larger sectors of the code that were then invaded by incoming amino acids. Displacements of amino acids follow selection rules. The code sectored from a glycine code to a four amino acid code to an eight amino acid code to an ~16 amino acid code to the standard 20 amino acid code with stops. The proposed patterns of code sectoring are now most apparent from patterns of aminoacyl-tRNA synthetase evolution. The Elongation Factor-Tu GTPase anticodon-codon latch that checks the accuracy of translation appears to have evolved at about the eight amino acid to ~16 amino acid stage. Before evolution of the EF-Tu latch, we posit that both the 1st and 3rd anticodon positions were wobble positions. The genetic code evolved via tRNA charging errors and via enzymatic modifications of amino acids joined to tRNAs, followed by tRNA and aminoacyl-tRNA synthetase differentiation. Fidelity mechanisms froze the code by inhibiting further innovation.

Broad range of missense error frequencies in cellular proteins

Assessment of the fidelity of gene expression is crucial to understand cell homeostasis. Here we present a highly sensitive method for the systematic Quantification of Rare Amino acid Substitutions (QRAS) using absolute quantification by targeted mass spectrometry after chromatographic enrichment of peptides with missense amino acid substitutions. By analyzing incorporation of near- and non-cognate amino acids in a model protein EF-Tu, we show that most of missense errors are too rare to detect by conventional methods, such as DDA, and are estimated to be between <10⁻⁷–10^-5 by QRAS. We also observe error hotspots of up to 10⁻³ for some types of mismatches, including the G-U mismatch. The error frequency depends on the expression level of EF-Tu and, surprisingly, the amino acid position in the protein. QRAS is not restricted to any particular miscoding event, organism, strain or model protein and is a reliable tool to analyze very rare proteogenomic events.

Ribosome collisions alter frameshifting at translational reprogramming motifs in bacterial mRNAs

Ribosomes move along mRNAs in 3-nucleotide steps as they interpret codons that specify which amino acid is required at each position in the protein. There are multiple examples of genes with DNA sequences that do not match the produced proteins because ribosomes move to a new reading frame in the message before finishing translation (so-called frameshifting). This report shows that, when ribosomes stall at mRNA regions prone to cause frameshifting events, trailing ribosomes that collide with them can significantly change the outcome and potentially regulate protein production. This work highlights the principle that biological macromolecules do not function in isolation, and it provides an example of how physical interactions between neighboring complexes can be used to augment their performance.

Translational frameshifting involves the repositioning of ribosomes on their messages into decoding frames that differ from those dictated during initiation. Some messenger RNAs (mRNAs) contain motifs that promote deliberate frameshifting to regulate production of the encoded proteins. The mechanisms of frameshifting have been investigated in many systems, and the resulting models generally involve single ribosomes responding to stimulator sequences in their engaged mRNAs. We discovered that the abundance of ribosomes on messages containing the IS3, dnaX, and prfB frameshift motifs significantly influences the levels of frameshifting. We show that this phenomenon results from ribosome collisions that occur during translational stalling, which can alter frameshifting in both the stalled and trailing ribosomes. Bacteria missing ribosomal protein bL9 are known to exhibit a reduction in reading frame maintenance and to have a strong dependence on elongation factor P (EFP). We discovered that ribosomes lacking bL9 become compacted closer together during collisions and that the E-sites of the stalled ribosomes appear to become blocked, which suggests subsequent transpeptidation in transiently stalled ribosomes may become compromised in the absence of bL9. In addition, we determined that bL9 can suppress frameshifting of its host ribosome, likely by regulating E-site dynamics. These findings provide mechanistic insight into the behavior of colliding ribosomes during translation and suggest naturally occurring frameshift elements may be regulated by the abundance of ribosomes relative to an mRNA pool.

A tRNA-mimic Strategy to Explore the Role of G34 of tRNAGly in Translation and Codon Frameshifting

Decoding of the 61 sense codons of the genetic code requires a variable number of tRNAs that establish codon-anticodon interactions. Thanks to the wobble base pairing at the third codon position, less than 61 different tRNA isoacceptors are needed to decode the whole set of codons. On the tRNA, a subtle distribution of nucleoside modifications shapes the anticodon loop structure and participates to accurate decoding and reading frame maintenance. Interestingly, although the 61 anticodons should exist in tRNAs, a strict absence of some tRNAs decoders is found in several codon families. For instance, in Eukaryotes, G34-containing tRNAs translating 3-, 4- and 6-codon boxes are absent. This includes tRNA specific for Ala, Arg, Ile, Leu, Pro, Ser, Thr, and Val. tRNA^Gly is the only exception for which in the three kingdoms, a G34-containing tRNA exists to decode C3 and U3-ending codons. To understand why G34-tRNA^Gly exists, we analysed at the genome wide level the codon distribution in codon +1 relative to the four GGN Gly codons. When considering codon GGU, a bias was found towards an unusual high usage of codons starting with a G whatever the amino acid at +1 codon. It is expected that GGU codons are decoded by G34-containing tRNA^Gly, decoding also GGC codons. Translation studies revealed that the presence of a G at the first position of the downstream codon reduces the +1 frameshift by stabilizing the G34•U3 wobble interaction. This result partially explains why G34-containing tRNA^Gly exists in Eukaryotes whereas all the other G34-containing tRNAs for multiple codon boxes are absent.

Unique tRNA gene profile suggests paucity of nucleotide modifications in anticodons of a deep-sea symbiotic Spiroplasma

The position 34 of a tRNA is always modified for efficient recognition of codons and accurate integration of amino acids by the translation machinery. Here, we report genomics features of a deep-sea gut symbiotic Spiroplasma, which suggests that the organism does not require tRNA(34) anticodon modifications. In the genome, there is a novel set of tRNA genes composed of 32 species for recognition of the 20 amino acids. Among the anticodons of the tRNAs, we witnessed the presence of both U34- and C34-containing tRNAs required to decode NNR (A/G) 2:2 codons as countermeasure of probable loss of anticodon modification genes. In the tRNA fragments detected in the gut transcriptome, mismatches expected to be caused by some tRNA modifications were not shown in their alignments with the corresponding genes. However, the probable paucity of modified anticodons did not fundamentally change the codon usage pattern of the Spiroplasma. The tRNA gene profile that probably resulted from the paucity of tRNA(34) modifications was not observed in other symbionts and deep-sea bacteria, indicating that this phenomenon was an evolutionary dead-end. This study provides insights on co-evolution of translation machine and tRNA genes and steric constraints of codon-anticodon interactions in deep-sea extreme environment.

Understanding the Genetic Code

The universal triple-nucleotide genetic code is often viewed as a given, randomly selected through evolution. However, as summarized in this article, many observations and deductions within structural and thermodynamic frameworks help to explain the forces that must have shaped the code during the early evolution of life on Earth.

Accuracy of genetic code translation and its orthogonal corruption by aminoglycosides and Mg2+ ions

We studied the effects of aminoglycosides and changing Mg2+ ion concentration on the accuracy of initial codon selection by aminoacyl-tRNA in ternary complex with elongation factor Tu and GTP (T3) on mRNA programmed ribosomes. Aminoglycosides decrease the accuracy by changing the equilibrium constants of 'monitoring bases' A1492, A1493 and G530 in 16S rRNA in favor of their 'activated' state by large, aminoglycoside-specific factors, which are the same for cognate and near-cognate codons. Increasing Mg2+ concentration decreases the accuracy by slowing dissociation of T3 from its initial codon- and aminoglycoside-independent binding state on the ribosome. The distinct accuracy-corrupting mechanisms for aminoglycosides and Mg2+ ions prompted us to re-interpret previous biochemical experiments and functional implications of existing high resolution ribosome structures. We estimate the upper thermodynamic limit to the accuracy, the 'intrinsic selectivity' of the ribosome. We conclude that aminoglycosides do not alter the intrinsic selectivity but reduce the fraction of it that is expressed as the accuracy of initial selection. We suggest that induced fit increases the accuracy and speed of codon reading at unaltered intrinsic selectivity of the ribosome.

Importance of potassium ions for ribosome structure and function revealed by long-wavelength X-ray diffraction

The ribosome, the largest RNA-containing macromolecular machinery in cells, requires metal ions not only to maintain its three-dimensional fold but also to perform protein synthesis. Despite the vast biochemical data regarding the importance of metal ions for efficient protein synthesis and the increasing number of ribosome structures solved by X-ray crystallography or cryo-electron microscopy, the assignment of metal ions within the ribosome remains elusive due to methodological limitations. Here we present extensive experimental data on the potassium composition and environment in two structures of functional ribosome complexes obtained by measurement of the potassium anomalous signal at the K-edge, derived from long-wavelength X-ray diffraction data. We elucidate the role of potassium ions in protein synthesis at the three-dimensional level, most notably, in the environment of the ribosome functional decoding and peptidyl transferase centers. Our data expand the fundamental knowledge of the mechanism of ribosome function and structural integrity.

Optimal translational fidelity is critical for Salmonella virulence and host interactions

Translational fidelity is required for accurate flow of genetic information, but is frequently altered by genetic changes and environmental stresses. To date, little is known about how translational fidelity affects the virulence and host interactions of bacterial pathogens. Here we show that surprisingly, either decreasing or increasing translational fidelity impairs the interactions of the enteric pathogen Salmonella Typhimurium with host cells and its fitness in zebrafish. Host interactions are mediated by Salmonella pathogenicity island 1 (SPI-1). Our RNA sequencing and quantitative RT-PCR results demonstrate that SPI-1 genes are among the most down-regulated when translational fidelity is either increased or decreased. Further, this down-regulation of SPI-1 genes depends on the master regulator HilD, and altering translational fidelity destabilizes HilD protein via enhanced degradation by Lon protease. Our work thus reveals that optimal translational fidelity is pivotal for adaptation of Salmonella to the host environment, and provides important mechanistic insights into this process.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

A tRNA- and Anticodon-Centric View of the Evolution of Aminoacyl-tRNA Synthetases, tRNAomes, and the Genetic Code

Pathways of standard genetic code evolution remain conserved and apparent, particularly upon analysis of aminoacyl-tRNA synthetase (aaRS) lineages. Despite having incompatible active site folds, class I and class II aaRS are homologs by sequence. Specifically, structural class IA aaRS enzymes derive from class IIA aaRS enzymes by in-frame extension of the protein N-terminus and by an alternate fold nucleated by the N-terminal extension. The divergence of aaRS enzymes in the class I and class II clades was analyzed using the Phyre2 protein fold recognition server. The class I aaRS radiated from the class IA enzymes, and the class II aaRS radiated from the class IIA enzymes. The radiations of aaRS enzymes bolster the coevolution theory for evolution of the amino acids, tRNAomes, the genetic code, and aaRS enzymes and support a tRNA anticodon-centric perspective. We posit that second- and third-position tRNA anticodon sequence preference (C>(U~G)>A) powerfully selected the sectoring pathway for the code. GlyRS-IIA appears to have been the primordial aaRS from which all aaRS enzymes evolved, and glycine appears to have been the primordial amino acid around which the genetic code evolved.

Substrate-Induced Formation of Ribosomal Decoding Center for Accurate and Rapid Genetic Code Translation

Accurate translation of genetic information is crucial for synthesis of functional proteins in all organisms. We use recent experimental data to discuss how induced fit affects accuracy of initial codon selection on the ribosome by aminoacyl transfer RNA in ternary complex ( T₃) with elongation factor Tu (EF-Tu) and guanosine-5'-triphosphate (GTP). We define actual accuracy ([Formula: see text]) of a particular protein synthesis system as its current accuracy and the effective selectivity ([Formula: see text]) as [Formula: see text] in the limit of zero ribosomal binding affinity for T₃. Intrinsic selectivity ([Formula: see text]), defined as the upper thermodynamic limit of [Formula: see text], is determined by the free energy difference between near-cognate and cognate T₃ in the pre-GTP hydrolysis state on the ribosome. [Formula: see text] is much larger than [Formula: see text], suggesting the possibility of a considerable increase in [Formula: see text] and [Formula: see text] at negligible kinetic cost. Induced fit increases [Formula: see text] and [Formula: see text] without affecting [Formula: see text], and aminoglycoside antibiotics reduce [Formula: see text] and [Formula: see text] at unaltered [Formula: see text].

The multiple flavors of GoU pairs in RNA

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.

Identification of Lysine Misincorporation at Asparagine Position in Recombinant Insulin Analogs Produced in E. coli

PURPOSE

Identification of human insulin analogs' impurity with a mass shift +14 Da in comparison to a parent protein.

METHODS

The protein sequence variant was detected and identified with the application of peptide mapping, liquid chromatography, tandem mass spectrometric analysis, nuclear magnetic resonance spectroscopy (NMR) and Edman sequencing.

RESULTS

The misincorporated lysine (Lys) at asparagine (Asn) position A21 was detected in recombinant human insulin and its analogs.

CONCLUSIONS

Although there are three asparagine residues in the insulin derivative, the misincorporation of lysine occurred only at position A21. The process involves G/U or A/U wobble base pairing.

A viral RNA motif involved in signaling the initiation of translation on non-AUG codons

Noncanonical translation, and particularly initiation on non-AUG codons, are frequently used by viral and cellular mRNAs during virus infection and disease. The Sindbis virus (SINV) subgenomic mRNA (sgRNA) constitutes a unique model system to analyze the translation of a capped viral mRNA without the participation of several initiation factors. Moreover, sgRNA can initiate translation even when the AUG initiation codon is replaced by other codons. Using SINV replicons, we examined the efficacy of different codons in place of AUG to direct the synthesis of the SINV capsid protein. The substitution of AUG by CUG was particularly efficient in promoting the incorporation of leucine or methionine in similar percentages at the amino terminus of the capsid protein. Additionally, valine could initiate translation when the AUG is replaced by GUG. The ability of sgRNA to initiate translation on non-AUG codons was dependent on the integrity of a downstream stable hairpin (DSH) structure located in the coding region. The structural requirements of this hairpin to signal the initiation site on the sgRNA were examined in detail. Of interest, a virus bearing CUG in place of AUG in the sgRNA was able to infect cells and synthesize significant amounts of capsid protein. This virus infects the human haploid cell line HAP1 and the double knockout variant that lacks eIF2A and eIF2D. Collectively, these findings indicate that leucine-tRNA or valine-tRNA can participate in the initiation of translation of sgRNA by a mechanism dependent on the DSH. This mechanism does not involve the action of eIF2, eIF2A, or eIF2D.

Errors in translational decoding: tRNA wobbling or misincorporation?

As the central dogma of molecular biology, genetic information flows from DNA through transcription into RNA followed by translation of the message into protein by transfer RNAs (tRNAs). However, mRNA translation is not always perfect, and errors in the amino acid composition may occur. Mistranslation is generally well tolerated, but once it reaches superphysiological levels, it can give rise to a plethora of diseases. The key causes of mistranslation are errors in translational decoding of the codons in mRNA. Such errors mainly derive from tRNA misdecoding and misacylation, especially when certain codon-paired tRNA species are missing. Substantial progress has recently been made with respect to the mechanistic basis of erroneous mRNA decoding as well as the resulting consequences for physiology and pathology. Here, we aim to review this progress with emphasis on viral evolution and cancer development.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

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.

Translation of non-standard codon nucleotides reveals minimal requirements for codon-anticodon interactions

The precise interplay between the mRNA codon and the tRNA anticodon is crucial for ensuring efficient and accurate translation by the ribosome. The insertion of RNA nucleobase derivatives in the mRNA allowed us to modulate the stability of the codon-anticodon interaction in the decoding site of bacterial and eukaryotic ribosomes, allowing an in-depth analysis of codon recognition. We found the hydrogen bond between the N1 of purines and the N3 of pyrimidines to be sufficient for decoding of the first two codon nucleotides, whereas adequate stacking between the RNA bases is critical at the wobble position. Inosine, found in eukaryotic mRNAs, is an important example of destabilization of the codon-anticodon interaction. Whereas single inosines are efficiently translated, multiple inosines, e.g., in the serotonin receptor 5-HT2C mRNA, inhibit translation. Thus, our results indicate that despite the robustness of the decoding process, its tolerance toward the weakening of codon-anticodon interactions is limited.

Codon-specific effects of tRNA anticodon loop modifications on translational misreading errors in the yeast Saccharomyces cerevisiae

Protein synthesis requires both high speed and accuracy to ensure a healthy cellular environment. Estimates of errors during protein synthesis in Saccharomyces cerevisiae have varied from 10-3 to 10-4 errors per codon. Here, we show that errors made by ${\rm{tRNA}}^{\rm Glu}_{\rm UUC}$ in yeast can vary 100-fold, from 10-6 to 10-4 errors per codon. The most frequent errors require a G•U mismatch at the second position for the near cognate codon GGA (Gly). We also show, contrary to our previous results, that yeast tRNAs can make errors involving mismatches at the wobble position but with low efficiency. We have also assessed the effect on misreading frequency of post-transcriptional modifications of tRNAs, which are known to regulate cognate codon decoding in yeast. We tested the roles of mcm5s2U34 and t6A37 and show that their effects depend on details of the codon anticodon interaction including the position of the modification with respect to the base mismatch and the nature of that mismatch. Both mcm5 and s2 modification of wobble uridine strongly stabilizes G2•U35 mismatches when ${\rm{tRNA}}^{\rm Glu}_{\rm UUC}$ misreads the GGA Gly codon but has weaker effects on other mismatches. By contrast, t6A37 destabilizes U1•U36 mismatches when ${\rm{tRNA}}^{\rm Lys}_{\rm UUU}$ misreads UAA or UAG but stabilizes mismatches at the second and wobble positions.

Co-adaption of tRNA gene copy number and amino acid usage influences translation rates in three life domains

Although more and more entangled participants of translation process were realized, how they cooperate and co-determine the final translation efficiency still lacks details. Here, we reasoned that the basic translation components, tRNAs and amino acids should be consistent to maximize the efficiency and minimize the cost. We firstly revealed that 310 out of 410 investigated genomes of three domains had significant co-adaptions between the tRNA gene copy numbers and amino acid compositions, indicating that maximum efficiency constitutes ubiquitous selection pressure on protein translation. Furthermore, fast-growing and larger bacteria are found to have significantly better co-adaption and confirmed the effect of this pressure. Within organism, highly expressed proteins and those connected to acute responses have higher co-adaption intensity. Thus, the better co-adaption probably speeds up the growing of cells through accelerating the translation of special proteins. Experimentally, manipulating the tRNA gene copy number to optimize co-adaption between enhanced green fluorescent protein (EGFP) and tRNA gene set of Escherichia coli indeed lifted the translation rate (speed). Finally, as a newly confirmed translation rate regulating mechanism, the co-adaption reflecting translation rate not only deepens our understanding on translation process but also provides an easy and practicable method to improve protein translation rates and productivity.

Tautomeric G•U pairs within the molecular ribosomal grip and fidelity of decoding in bacteria

We report new crystallographic structures of Thermus thermophilus ribosomes complexed with long mRNAs and native Escherichia coli tRNAs. They complete the full set of combinations of Watson-Crick G•C and miscoding G•U pairs at the first two positions of the codon-anticodon duplex in ribosome functional complexes. Within the tight decoding center, miscoding G•U pairs occur, in all combinations, with a non-wobble geometry structurally indistinguishable from classical coding Watson-Crick pairs at the same first two positions. The contacts with the ribosomal grip surrounding the decoding center are all quasi-identical, except in the crowded environment of the amino group of a guanosine at the second position; in which case a G in the codons may be preferred. In vivo experimental data show that the translational errors due to miscoding by G•U pairs at the first two positions are the most frequently encountered ones, especially at the second position and with a G on the codon. Such preferred miscodings involve a switch from an A-U to a G•U pair in the tRNA/mRNA complex and very rarely from a G = C to a G•U pair. It is concluded that the frequencies of such occurrences are only weakly affected by the codon/anticodon structures but depend mainly on the stability and lifetime of the complex, the modifications present in the anticodon loop, especially those at positions 34 and 37, in addition to the relative concentration of cognate/near-cognate tRNA species present in the cellular tRNA pool.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.