Conformational Preferences of Modified Nucleoside 5-Taurinomethyluridine, τm(5)U Occur at 'wobble' 34th Position in the Anticodon Loop of tRNA

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

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

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

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

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

Epigenetic and epitranscriptomic regulation of viral replication

Eukaryotic gene expression is regulated not only by genomic enhancers and promoters, but also by covalent modifications added to both chromatin and RNAs. Whereas cellular gene expression may be either enhanced or inhibited by specific epigenetic modifications deposited on histones (in particular, histone H3), these epigenetic modifications can also repress viral gene expression, potentially functioning as a potent antiviral innate immune response in DNA virus-infected cells. However, viruses have evolved countermeasures that prevent the epigenetic silencing of their genes during lytic replication, and they can also take advantage of epigenetic silencing to establish latent infections. By contrast, the various covalent modifications added to RNAs, termed epitranscriptomic modifications, can positively regulate mRNA translation and/or stability, and both DNA and RNA viruses have evolved to utilize epitranscriptomic modifications as a means to maximize viral gene expression. As a consequence, both chromatin and RNA modifications could serve as novel targets for the development of antivirals. In this Review, we discuss how host epigenetic and epitranscriptomic processes regulate viral gene expression at the levels of chromatin and RNA function, respectively, and explore how viruses modify, avoid or utilize these processes in order to regulate viral gene expression.

The chemical diversity of RNA modifications

Nucleic acid modifications in DNA and RNA ubiquitously exist among all the three kingdoms of life. This trait significantly broadens the genome diversity and works as an important means of gene transcription regulation. Although mammalian systems have limited types of DNA modifications, over 150 different RNA modification types have been identified, with a wide variety of chemical diversities. Most modifications occur on transfer RNA and ribosomal RNA, however many of the modifications also occur on other types of RNA species including mammalian mRNA and small nuclear RNA, where they are essential for many biological roles, including developmental processes and stem cell differentiation. These post-transcriptional modifications are enzymatically installed and removed in a site-specific manner by writer and eraser proteins respectively, while reader proteins can interpret modifications and transduce the signal for downstream functions. Dysregulation of mRNA modifications manifests as disease states, including multiple types of human cancer. In this review, we will introduce the chemical features and biological functions of these modifications in the coding and non-coding RNA species.

Structural significance of hypermodified nucleoside 5-carboxymethylaminomethyluridine (cmnm5U) from 'wobble' (34th) position of mitochondrial tRNAs: Molecular modeling and Markov state model studies

A quantum chemical semi-empirical RM1 approach was used to deduce the structural role of hypermodified nucleoside 5-carboxymethylaminomethyluridine 5'-monophosphate (pcmnm⁵U) from 'wobble' (34th) position of mitochondrial tRNAs. The energetically preferred pcmnm⁵U₍₃₄₎ adopted a 'skew' conformation for C5-substituted side chain (-CH₂-NH₂⁺-CH₂-COO^-) moiety that orient towards the 5'-ribose-phosphate backbone, which support 'anti' orientation of glycosyl (χ₃₄) torsion angle. Preferred conformation of pcmnm⁵U₍₃₄₎ was stabilized by O(4) … HC(10), O1P⋯HN(11), O(15) … HN(11), O(15) … HC(10), O4' … HC(6) and O(2) … HC2' hydrogen bonding interactions. The high flexibility of side chain moiety displayed different structural properties for pcmnm⁵U₍₃₄₎. Three different conformations of pcmnm⁵U₍₃₄₎ were observed in molecular dynamics simulations and Markov state model studies. The unmodified uracil revealed 'syn' and 'anti' orientations for glycosyl (χ₃₄) torsion angle that substantiate the role of "-CH₂-NH₂⁺-CH₂-COO^-" moiety in maintaining the 'anti' orientation of pcmnm⁵U₍₃₄₎. The preferred conformation of pcmnm⁵U₍₃₄₎ helps to recognize Guanosine more proficiently than Adenosine from the third position of codons. The role of pcmnm⁵U₍₃₄₎ in tRNA biogenesis paves the way to understand its structural significance in usual mitochondrial metabolism and respiration.

A molecular mechanics study on GA codon box translation

The GA codon box incorporates the two-fold degeneracy of aspartic acid and of glutamic acid. Using the molecular mechanics approach of the AMBER suite, the four codons of the GA box are paired via H-bonding with two aspartic acid anticodons and two glutamic acid anticodons to yield 8 cognate and 11 non-cognate codon-anticodon duplexes. In addition four select non-cognate duplexes between the GA box codons and three alanine anticodons are also studied. These 23 duplexes display a variety of base-pairing possibilities at the wobble position. Cognate duplexes are differentiated from non-cognate duplexes on the grounds of structure and stability (chiefly the former). The results are in line with Crick's wobble hypothesis, and corroborate the observed reading properties of the aspartic acid anticodons GUC and QUC and of the glutamic acid anticodons CUC and S^mnUC.

Preferences of AAA/AAG codon recognition by modified nucleosides, τm5s2U34 and t6A37 present in tRNALys

Deficiency of 5-taurinomethyl-2-thiouridine, τm⁵s²U at the 34th 'wobble' position in tRNA^Lys causes MERRF (Myoclonic Epilepsy with Ragged Red Fibers), a neuromuscular disease. This modified nucleoside of mt tRNA^Lys, recognizes AAA/AAG codons during protein biosynthesis process. Its preference to identify cognate codons has not been studied at the atomic level. Hence, multiple MD simulations of various molecular models of anticodon stem loop (ASL) of mt tRNA^Lys in presence and absence of τm⁵s²U₃₄ and N⁶-threonylcarbamoyl adenosine (t⁶A₃₇) along with AAA and AAG codons have been accomplished. Additional four MD simulations of multiple ASL mt tRNA^Lys models in the context of ribosomal A-site residues have also been performed to investigate the role of A-site in recognition of AAA/AAG codons. MD simulation results show that, ASL models in presence of τm⁵s²U₃₄ and t⁶A₃₇ with codons AAA/AAG are more stable than the ASL lacking these modified bases. MD trajectories suggest that τm⁵s²U recognizes the codons initially by 'wobble' hydrogen bonding interactions, and then tRNA^Lys might leave the explicit codon by a novel 'single' hydrogen bonding interaction in order to run the protein biosynthesis process smoothly. We propose this model as the 'Foot-Step Model' for codon recognition, in which the single hydrogen bond plays a crucial role. MD simulation results suggest that, tRNA^Lys with τm⁵s²U and t⁶A recognizes AAA codon more preferably than AAG. Thus, these results reveal the consequences of τm⁵s²U and t⁶A in recognition of AAA/AAG codons in mitochondrial disease, MERRF.

Conformational preferences and structural analysis of hypermodified nucleoside, peroxywybutosine (o2yW) found at 37th position in anticodon loop of tRNAPhe and its role in modulating UUC codon-anticodon interactions

Hypermodified bases present at 3'-adjacent (37^th) position in anticodon loop of tRNA^Phe are well known for their contribution in modulating codon-anticodon interactions. Peroxywybutosine (o2yW), a wyosine family member, is one of such tricyclic modified bases observed at the 37^th position in tRNA^Phe. Conformational preferences and three-dimensional structural analysis of peroxywybutosine have not been investigated in detail at atomic level. Hence, in the present study quantum chemical semi-empirical RM1 and multiple molecular dynamics (MD) simulations have been used to study structural significance of peroxywybutosine in tRNA^Phe. Full geometry optimizations over the peroxywybutosine base have also been performed using ab-initio HF-SCF (6-31G**), DFT (B3LYP/6-31G**) and semi-empirical PM6 method to compare the salient properties. RM1 predicted most stable structure shows that the amino-carboxy-propyl side chain of o2yW remains 'distal' to the five membered imidazole ring of tricyclic guanosine. MD simulation trajectory of the isolated peroxy base showed restricted periodical fluctuations of peroxywybutosine side chain which might be helpful to maintain proper anticodon loop structure and mRNA reading frame during protein biosynthesis process. Another comparative MD simulation study of the anticodon stem loop with codon UUC showed various properties, which justify the functional implications of peroxywybutosine at 37^th position along with other modified bases present in ASL of tRNA^Phe. Thus, this study presents an atomic view into the structural properties of peroxywybutosine, which can be useful to determine its role in the anticodon stem loop in context of codon-anticodon interactions and frame shift mutations.

Idiosyncratic recognition of UUG/UUA codons by modified nucleoside 5-taurinomethyluridine, τm5U present at 'wobble' position in anticodon loop of tRNALeu: A molecular modeling approach

Lack of naturally occurring modified nucleoside 5-taurinomethyluridine (τm5U) at the 'wobble' 34th position in tRNALeu causes mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes (MELAS). The τm5U34 specifically recognizes UUG and UUA codons. Structural consequences of τm5U34 to read cognate codons have not been studied so far in detail at the atomic level. Hence, 50ns multiple molecular dynamics (MD) simulations of various anticodon stem loop (ASL) models of tRNALeu in presence and absence of τm5U34 along with UUG and UUA codons were performed to explore the dynamic behaviour of τm5U34 during codon recognition process. The MD simulation results revealed that τm5U34 recognizes G/A ending codons by 'wobble' as well as a novel 'single' hydrogen bonding interactions. RMSD and RMSF values indicate the comparative stability of the ASL models containing τm5U34 modification over the other models, lacking τm5U34. Another MD simulation study of 55S mammalian mitochondrial rRNA with tRNALeu showed crucial interactions between the A-site residues, A918, A919, G256 and codon-anticodon bases. Thus, these results could improve our understanding about the decoding efficiency of human mt tRNALeu with τm5U34 to recognize UUG and UUA codons.

Chemical and Conformational Diversity of Modified Nucleosides Affects tRNA Structure and Function

RNAs are central to all gene expression through the control of protein synthesis. Four major nucleosides, adenosine, guanosine, cytidine and uridine, compose RNAs and provide sequence variation, but are limited in contributions to structural variation as well as distinct chemical properties. The ability of RNAs to play multiple roles in cellular metabolism is made possible by extensive variation in length, conformational dynamics, and the over 100 post-transcriptional modifications. There are several reviews of the biochemical pathways leading to RNA modification, but the physicochemical nature of modified nucleosides and how they facilitate RNA function is of keen interest, particularly with regard to the contributions of modified nucleosides. Transfer RNAs (tRNAs) are the most extensively modified RNAs. The diversity of modifications provide versatility to the chemical and structural environments. The added chemistry, conformation and dynamics of modified nucleosides occurring at the termini of stems in tRNA's cloverleaf secondary structure affect the global three-dimensional conformation, produce unique recognition determinants for macromolecules to recognize tRNAs, and affect the accurate and efficient decoding ability of tRNAs. This review will discuss the impact of specific chemical moieties on the structure, stability, electrochemical properties, and function of tRNAs.

The influence of hypermodified nucleosides lysidine and t(6)A to recognize the AUA codon instead of AUG: a molecular dynamics simulation study

Hypermodified nucleosides lysidine (L) and N(6)-threonylcarbamoyladenosine (t(6)A) influence codon-anticodon interactions during the protein biosynthesis process. Lysidine prevents the misrecognition of the AUG codon as isoleucine and that of AUA as methionine. The structural significance of these modified bases has not been studied in detail at the atomic level. Hence, in the present study we performed multiple molecular dynamics (MD) simulations of anticodon stem loop (ASL) of tRNA(Ile) in the presence and absence of modified bases 'L' and 't(6)A' at the 34th and 37th positions respectively along with trinucleotide 'AUA' and 'AUG' codons. Hydrogen bonding interactions formed by the tautomeric form of lysidine may assist in reading the third base adenine of the 'AUA' codon, unlike the guanine of the 'AUG' codon. Such interactions might be useful to restrict codon specificity to recognize isoleucine tRNA instead of methionine tRNA. The t(6)A side chain interacts with the purine ring of the first codon nucleotide adenine, which might provide base stacking interactions and could be responsible for restricting extended codon-anticodon recognition. We found that ASL tRNA(Ile) in the absence of modifications at the 34th and 37th positions cannot establish proper hydrogen bonding interactions to recognize the isoleucine codon 'AUA' and subsequently disturbs the anticodon loop structure. The binding free energy calculations revealed that tRNA(Ile) ASL with modified nucleosides prefers the codon AUA over AUG. Thus, these findings might be useful to understand the role of modified bases L and t(6)A to recognize the AUA codon instead of AUG.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.