Structural domains of transfer RNA molecules

Synthesis and antiproliferative evaluation of novel 3,5,8-trisubstituted coumarins against breast cancer

Aim: This study focused on designing and synthesizing novel derivatives of 3,5,8-trisubstituted coumarin. Results: The synthesized compounds, particularly compound 5, exhibited significant cytotoxic effects on MCF-7 cells, surpassing staurosporine, and reduced toxicity toward MCF-10A cells, highlighting potential pharmacological advantages. Further, compound 5 altered the cell cycle and significantly increased apoptosis in MCF-7 cells, involving both early (41.7-fold) and late stages (33-fold), while moderately affecting necrotic signaling. The antitumor activity was linked to a notable reduction (4.78-fold) in topoisomerase IIβ expression. Molecular modeling indicated compound 5's strong affinity for EGFR, human EGF2 and topoisomerase II proteins. Conclusion: These findings highlight compound 5 as a multifaceted antitumor agent for breast cancer.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Synthesis, NMR kinetics and dynamic structure of a 17-mer heptaloop RNA hairpin carrying a 3-N-methyluridine nucleotide residue in the loop region

A 17-mer RNA hairpin (5'GGGAGUXAGCGGCUCCC3') carrying 3-N-methyluridine (m3U) at position X (m3U7-RNA), designed to represent the anticodon stem-loop (ACSL) region of tRNAs to study an open loop state (O-state), was synthesized, purified by HPLC, and characterized by MALDI-ToF_MS and NMR methods. 1H-NMR data revealed primary (P-state in 56.1%), secondary (S-state in 43.9%) and tertiary (∼5-6%) ACSL conformations. Exchange rate constant (kex) for interconversion between P and S states is 112 sec-1 (<Δω ∼454 rad/sec), confirming a slow exchange regime between the two states. Forward (kPS) and backward (kSP) rate constants are 49.166 sec-1 and 62.792 sec-1, respectively, leading to a longer life-time (20.339 msec) for the P-state and a shorter life-time (15.926 msec) for the S-state. In accordance with conformational populations determined by 1H-NMR, dynamics of the P/S/tertiary states of m3U7-RNA and its wild-type counterpart (wt-RNA) were studied by three independent MD production simulations. Cluster analysis revealed that wt-RNA reflects the structural characteristics of the ACSL region of tRNAs. The P-state of m3U7-RNA was found to be structurally similar to wt-RNA but lacks an intraloop H-bond between m3U7 and C10 (U33 and nt36 in tRNAs). In the S-state of m3U7-RNA, m3U7 flips out of the loop region. O-state loop conformations of m3U7-RNA were also clustered (4.8%), wherein the loop nucleotides m3U7.A8.G9.C10.G11 stack one after another. We propose that the O-state of m3U7-RNA is the most suitable conformation that makes the loop accessible for complementary nucleotides and for non-enzymatic primordial replication of small circular RNAs.Communicated by Ramaswamy H. Sarma.

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

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

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

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

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

Structural and dynamic effects of pseudouridine modifications on noncanonical interactions in RNA

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.

General Strategies for RNA X-ray Crystallography

An extremely small proportion of the X-ray crystal structures deposited in the Protein Data Bank are of RNA or RNA-protein complexes. This is due to three main obstacles to the successful determination of RNA structure: (1) low yields of pure, properly folded RNA; (2) difficulty creating crystal contacts due to low sequence diversity; and (3) limited methods for phasing. Various approaches have been developed to address these obstacles, such as native RNA purification, engineered crystallization modules, and incorporation of proteins to assist in phasing. In this review, we will discuss these strategies and provide examples of how they are used in practice.

Structural basis for sequence-independent substrate selection by eukaryotic wobble base tRNA deaminase ADAT2/3

The essential deamination of adenosine A34 to inosine at the wobble base is the individual tRNA modification with the greatest effects on mRNA decoding, empowering a single tRNA to translate three different codons. To date, many aspects of how eukaryotic deaminases specifically select their multiple substrates remain unclear. Here, using cryo-EM, we present the structure of a eukaryotic ADAT2/3 deaminase bound to a full-length tRNA, revealing that the enzyme distorts the anticodon loop, but in contrast to the bacterial enzymes, selects its substrate via sequence-independent contacts of eukaryote-acquired flexible or intrinsically unfolded motifs distal from the conserved catalytic core. A gating mechanism for substrate entry to the active site is identified. Our multi-step tRNA recognition model yields insights into how RNA editing by A34 deamination evolved, shaped the genetic code, and directly impacts the eukaryotic proteome.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.

The bacterial yjdF riboswitch regulates translation through its tRNA-like fold

Unlike most riboswitches, which have one cognate effector, the bacterial yjdF riboswitch binds to diverse azaaromatic compounds, only a subset of which cause it to activate translation. We examined the yjdF aptamer domain by small-angle X-ray scattering, and found that in the presence of activating ligands, the RNA adopts an overall shape similar to that of tRNA. Sequence analyses suggested that the yjdF aptamer is a homolog of tRNA^Lys, and that two of the conserved loops of the riboswitch are equivalent to the D- and T-loops of tRNA, associating to form an elbow-like tertiary interaction. Chemical probing indicated that this association is promoted by activating ligands such as chelerythrine and harmine. In its native mRNA context, activator ligands stabilize the tRNA-like fold of the yjdF aptamer, outcompeting the attenuated state in which its T-loop base-pairs to the Shine-Dalgarno element of the mRNA. Moreover, we demonstrate that the liganded aptamer itself activates translation, as authentic tRNAs, when grafted into mRNA, can potently activate translation. Taken together, our data demonstrate the ability of tRNA to function as a small-molecule responsive cis regulatory element.

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

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

Genetic code degeneracy is established by the decoding center of the ribosome

The degeneracy of the genetic code confers a wide array of properties to coding sequences. Yet, its origin is still unclear. A structural analysis has shown that the stability of the Watson–Crick base pair at the second position of the anticodon–codon interaction is a critical parameter controlling the extent of non-specific pairings accepted at the third position by the ribosome, a flexibility at the root of degeneracy. Based on recent cryo-EM analyses, the present work shows that residue A1493 of the decoding center provides a significant contribution to the stability of this base pair, revealing that the ribosome is directly involved in the establishment of degeneracy. Building on existing evolutionary models, we show the evidence that the early appearance of A1493 and A1492 established the basis of degeneracy when an elementary kinetic scheme of translation was prevailing. Logical considerations on the expansion of this kinetic scheme indicate that the acquisition of the peptidyl transferase center was the next major evolutionary step, while the induced-fit mechanism, that enables a sharp selection of the tRNAs, necessarily arose later when G530 was acquired by the decoding center.

Cooperativity and Interdependency between RNA Structure and RNA-RNA Interactions

Complex RNA–RNA interactions are increasingly known to play key roles in numerous biological processes from gene expression control to ribonucleoprotein granule formation. By contrast, the nature of these interactions and characteristics of their interfaces, especially those that involve partially or wholly structured RNAs, remain elusive. Herein, we discuss different modalities of RNA–RNA interactions with an emphasis on those that depend on secondary, tertiary, or quaternary structure. We dissect recently structurally elucidated RNA–RNA complexes including RNA triplexes, riboswitches, ribozymes, and reverse transcription complexes. These analyses highlight a reciprocal relationship that intimately links RNA structure formation with RNA–RNA interactions. The interactions not only shape and sculpt RNA structures but also are enabled and modulated by the structures they create. Understanding this two-way relationship between RNA structure and interactions provides mechanistic insights into the expanding repertoire of noncoding RNA functions, and may inform the design of novel therapeutics that target RNA structures or interactions.

Unique anticodon loop conformation with the flipped-out wobble nucleotide in the crystal structure of unbound tRNAVal

During protein synthesis on ribosome, tRNA recognizes its cognate codon of mRNA through base-pairing with the anticodon. The 5'-end nucleotide of the anticodon is capable of wobble base-pairing, offering a molecular basis for codon degeneracy. The wobble nucleotide is often targeted for post-transcriptional modification, which affects the specificity and fidelity of the decoding process. Flipping-out of a wobble nucleotide in the anticodon loop has been proposed to be necessary for modifying enzymes to access the target nucleotide, which has been captured in selective structures of protein-bound complexes. Meanwhile, all other structures of free or ribosome-bound tRNA display anticodon bases arranged in stacked conformation. We report the X-ray crystal structure of unbound tRNA^Val1 to a 2.04 Å resolution showing two different conformational states of wobble uridine in the anticodon loop, one stacked on the neighboring base and the other swiveled out toward solvent. In addition, the structure reveals a rare magnesium ion coordination to the nitrogen atom of a nucleobase, which has been sampled very rarely among known structures of nucleic acids.

Structural model of the M7G46 Methyltransferase TrmB in complex with tRNA

TrmB belongs to the class I S-adenosylmethionine (SAM)-dependent methyltransferases (MTases) and introduces a methyl group to guanine at position 7 (m⁷G) in tRNA. In tRNAs m⁷G is most frequently found at position 46 in the variable loop and forms a tertiary base pair with C13 and U22, introducing a positive charge at G46. The TrmB/Trm8 enzyme family is structurally diverse, as TrmB proteins exist in a monomeric, homodimeric, and heterodimeric form. So far, the exact enzymatic mechanism, as well as the tRNA-TrmB crystal structure is not known. Here we present the first crystal structures of B. subtilis TrmB in complex with SAM and SAH. The crystal structures of TrmB apo and in complex with SAM and SAH have been determined by X-ray crystallography to 1.9 Å (apo), 2.5 Å (SAM), and 3.1 Å (SAH). The obtained crystal structures revealed Tyr193 to be important during SAM binding and MTase activity. Applying fluorescence polarization, the dissociation constant K_d of TrmB and tRNA^Phe was determined to be 0.12 µM ± 0.002 µM. Luminescence-based methyltransferase activity assays revealed cooperative effects during TrmB catalysis with half-of-the-site reactivity at physiological SAM concentrations. Structural data retrieved from small-angle x-ray scattering (SAXS), mass-spectrometry of cross-linked complexes, and molecular docking experiments led to the determination of the TrmB-tRNA^Phe complex structure.

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.

Post-Transcriptional Modifications of Conserved Nucleotides in the T-Loop of tRNA: A Tale of Functional Convergent Evolution

The high conservation of nucleotides of the T-loop, including their chemical identity, are hallmarks of tRNAs from organisms belonging to the three Domains of Life. These structural characteristics allow the T-loop to adopt a peculiar intraloop conformation able to interact specifically with other conserved residues of the D-loop, which ultimately folds the mature tRNA in a unique functional canonical L-shaped architecture. Paradoxically, despite the high conservation of modified nucleotides in the T-loop, enzymes catalyzing their formation depend mostly on the considered organism, attesting for an independent but convergent evolution of the post-transcriptional modification processes. The driving force behind this is the preservation of a native conformation of the tRNA elbow that underlies the various interactions of tRNA molecules with different cellular components.

Mitochondrial tRNA Mutations Associated With Essential Hypertension: From Molecular Genetics to Function

Essential hypertension (EH) is one of the most common cardiovascular diseases worldwide, entailing a high level of morbidity. EH is a multifactorial disease influenced by both genetic and environmental factors, including mitochondrial DNA (mtDNA) genotype. Previous studies identified mtDNA mutations that are associated with maternally inherited hypertension, including tRNA^Ile m.4263A>G, m.4291T>C, m.4295A>G, tRNA^Met m.4435A>G, tRNA^Ala m.5655A>G, and tRNA^Met/tRNA^Gln m.4401A>G, et al. These mtDNA mutations alter tRNA structure, thereby leading to metabolic disorders. Metabolic defects associated with mitochondrial tRNAs affect protein synthesis, cause oxidative phosphorylation defects, reduced ATP synthesis, and increase production of reactive oxygen species. In this review we discuss known mutations of tRNA genes encoded by mtDNA and the potential mechanisms by which these mutations may contribute to hypertension.

A stress-induced tyrosine-tRNA depletion response mediates codon-based translational repression and growth suppression

Eukaryotic transfer RNAs can become selectively fragmented upon various stresses, generating tRNA-derived small RNA fragments. Such fragmentation has been reported to impact a small fraction of the tRNA pool and thus presumed to not directly impact translation. We report that oxidative stress can rapidly generate tyrosine-tRNAGUA fragments in human cells-causing significant depletion of the precursor tRNA. Tyrosine-tRNAGUA depletion impaired translation of growth and metabolic genes enriched in cognate tyrosine codons. Depletion of tyrosine tRNAGUA or its translationally regulated targets USP3 and SCD repressed proliferation-revealing a dedicated tRNA-regulated growth-suppressive pathway for oxidative stress response. Tyrosine fragments are generated in a DIS3L2 exoribonuclease-dependent manner and inhibit hnRNPA1-mediated transcript destabilization. Moreover, tyrosine fragmentation is conserved in C. elegans. Thus, tRNA fragmentation can coordinately generate trans-acting small RNAs and functionally deplete a tRNA. Our findings reveal the existence of an underlying adaptive codon-based regulatory response inherent to the genetic code.

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature T_m and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

Characterization of UVA-Induced Alterations to Transfer RNA Sequences

Ultraviolet radiation (UVR) adversely affects the integrity of DNA, RNA, and their nucleoside modifications. By employing liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based RNA modification mapping approaches, we identified the transfer RNA (tRNA) regions most vulnerable to photooxidation. Photooxidative damage to the anticodon and variable loop regions was consistently observed in both modified and unmodified sequences of tRNA upon UVA (λ 370 nm) exposure. The extent of oxidative damage measured in terms of oxidized guanosine, however, was higher in unmodified RNA compared to its modified version, suggesting an auxiliary role for nucleoside modifications. The type of oxidation product formed in the anticodon stem-loop region varied with the modification type, status, and whether the tRNA was inside or outside the cell during exposure. Oligonucleotide-based characterization of tRNA following UVA exposure also revealed the presence of novel photoproducts and stable intermediates not observed by nucleoside analysis alone. This approach provides sequence-specific information revealing potential hotspots for UVA-induced damage in tRNAs.

Unusual tertiary pairs in eukaryotic tRNAAla

tRNA molecules have well-defined sequence conservations that reflect the conserved tertiary pairs maintaining their architecture and functions during the translation processes. An analysis of aligned tRNA sequences present in the GtRNAdb database (the Lowe Laboratory, University of California, Santa Cruz) led to surprising conservations on some cytosolic tRNAs specific for alanine compared to other tRNA species, including tRNAs specific for glycine. First, besides the well-known G3oU70 base pair in the amino acid stem, there is the frequent occurrence of a second wobble pair at G30oU40, a pair generally observed as a Watson-Crick pair throughout phylogeny. Second, the tertiary pair R15/Y48 occurs as a purine-purine R15/A48 pair. Finally, the conserved T54/A58 pair maintaining the fold of the T-loop is observed as a purine-purine A54/A58 pair. The R15/A48 and A54/A58 pairs always occur together. The G30oU40 pair occurs alone or together with these other two pairs. The pairing variations are observed to a variable extent depending on phylogeny. Among eukaryotes, insects display all variations simultaneously, whereas mammals present either the G30oU40 pair or both R15/A48 and A54/A58. tRNAs with the anticodon 34A(I)GC36 are the most prone to display all those pair variations in mammals and insects. tRNAs with anticodon Y34GC36 have preferentially G30oU40 only. These unusual pairs are not observed in bacterial, nor archaeal, tRNAs, probably because of the avoidance of A34-containing anticodons in four-codon boxes. Among eukaryotes, these unusual pairing features were not observed in fungi and nematodes. These unusual structural features may affect, besides aminoacylation, transcription rates (e.g., 54/58) or ribosomal translocation (30/40).

A Structural Basis for Restricted Codon Recognition Mediated by 2-thiocytidine in tRNA Containing a Wobble Position Inosine

Three of six arginine codons (CGU, CGC, and CGA) are decoded by two Escherichia coli tRNA^Arg isoacceptors. The anticodon stem and loop (ASL) domains of tRNA^Arg1 and tRNA^Arg2 both contain inosine and 2-methyladenosine modifications at positions 34 (I₃₄) and 37 (m²A₃₇). tRNA^Arg1 is also modified from cytidine to 2-thiocytidine at position 32 (s²C₃₂). The s²C₃₂ modification is known to negate wobble codon recognition of the rare CGA codon by an unknown mechanism, while still allowing decoding of CGU and CGC. Substitution of s²C₃₂ for C₃₂ in the Saccharomyces cerevisiae tRNA^Ile_IAU anticodon stem and loop domain (ASL) negates wobble decoding of its synonymous A-ending codon, suggesting that this function of s²C at position 32 is a generalizable property. X-ray crystal structures of variously modified ASL^Arg1_ICG and ASL^Arg2_ICG constructs bound to cognate and wobble codons on the ribosome revealed the disruption of a C_32-A₃₈ cross-loop interaction but failed to fully explain the means by which s²C₃₂ restricts I₃₄ wobbling. Computational studies revealed that the adoption of a spatially broad inosine-adenosine base pair at the wobble position of the codon cannot be maintained simultaneously with the canonical ASL U-turn motif. C_32-A₃₈ cross-loop interactions are required for stability of the anticodon/codon interaction in the ribosomal A-site.

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.

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.

The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA

The factors and mechanisms that govern tRNA stability in bacteria are not well understood. Here, we investigated the influence of posttranscriptional modification of bacterial tRNAs (tRNA modification) on tRNA stability. We focused on ThiI-generated 4-thiouridine (s⁴U), a modification found in bacterial and archaeal tRNAs. Comprehensive quantification of Vibrio cholerae tRNAs revealed that the abundance of some tRNAs is decreased in a ΔthiI strain in a stationary phase-specific manner. Multiple mechanisms, including rapid degradation of a subset of hypomodified tRNAs, account for the reduced abundance of tRNAs in the absence of thiI Through transposon insertion sequencing, we identified additional tRNA modifications that promote tRNA stability and bacterial viability. Genetic analysis of suppressor mutants as well as biochemical analyses revealed that rapid degradation of hypomodified tRNA is mediated by the RNA degradosome. Elongation factor Tu seems to compete with the RNA degradosome, protecting aminoacyl tRNAs from decay. Together, our observations describe a previously unrecognized bacterial tRNA quality control system in which hypomodification sensitizes tRNAs to decay mediated by the RNA degradosome.

7-Methylguanosine Modifications in Transfer RNA (tRNA)

More than 90 different modified nucleosides have been identified in tRNA. Among the tRNA modifications, the 7-methylguanosine (m⁷G) modification is found widely in eubacteria, eukaryotes, and a few archaea. In most cases, the m⁷G modification occurs at position 46 in the variable region and is a product of tRNA (m⁷G46) methyltransferase. The m⁷G46 modification forms a tertiary base pair with C13-G22, and stabilizes the tRNA structure. A reaction mechanism for eubacterial tRNA m⁷G methyltransferase has been proposed based on the results of biochemical, bioinformatic, and structural studies. However, an experimentally determined mechanism of methyl-transfer remains to be ascertained. The physiological functions of m⁷G46 in tRNA have started to be determined over the past decade. For example, tRNA m⁷G46 or tRNA (m⁷G46) methyltransferase controls the amount of other tRNA modifications in thermophilic bacteria, contributes to the pathogenic infectivity, and is also associated with several diseases. In this review, information of tRNA m⁷G modifications and tRNA m⁷G methyltransferases is summarized and the differences in reaction mechanism between tRNA m⁷G methyltransferase and rRNA or mRNA m⁷G methylation enzyme are discussed.

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.

Intrinsic Properties of tRNA Molecules as Deciphered via Bayesian Network and Distribution Divergence Analysis

The identity/recognition of tRNAs, in the context of aminoacyl tRNA synthetases (and other molecules), is a complex phenomenon that has major implications ranging from the origins and evolution of translation machinery and genetic code to the evolution and speciation of tRNAs themselves to human mitochondrial diseases to artificial genetic code engineering. Deciphering it via laboratory experiments, however, is difficult and necessarily time- and resource-consuming. In this study, we propose a mathematically rigorous two-pronged in silico approach to identifying and classifying tRNA positions important for tRNA identity/recognition, rooted in machine learning and information-theoretic methodology. We apply Bayesian Network modeling to elucidate the structure of intra-tRNA-molecule relationships, and distribution divergence analysis to identify meaningful inter-molecule differences between various tRNA subclasses. We illustrate the complementary application of these two approaches using tRNA examples across the three domains of life, and identify and discuss important (informative) positions therein. In summary, we deliver to the tRNA research community a novel, comprehensive methodology for identifying the specific elements of interest in various tRNA molecules, which can be followed up by the corresponding experimental work and/or high-resolution position-specific statistical analyses.

Rooted tRNAomes and evolution of the genetic code

We advocate for a tRNA- rather than an mRNA-centric model for evolution of the genetic code. The mechanism for evolution of cloverleaf tRNA provides a root sequence for radiation of tRNAs and suggests a simplified understanding of code evolution. To analyze code sectoring, rooted tRNAomes were compared for several archaeal and one bacterial species. Rooting of tRNAome trees reveals conserved structures, indicating how the code was shaped during evolution and suggesting a model for evolution of a LUCA tRNAome tree. We propose the polyglycine hypothesis that the initial product of the genetic code may have been short chain polyglycine to stabilize protocells. In order to describe how anticodons were allotted in evolution, the sectoring-degeneracy hypothesis is proposed. Based on sectoring, a simple stepwise model is developed, in which the code sectors from a 1→4→8→∼16 letter code. At initial stages of code evolution, we posit strong positive selection for wobble base ambiguity, supporting convergence to 4-codon sectors and ∼16 letters. In a later stage, ∼5–6 letters, including stops, were added through innovating at the anticodon wobble position. In archaea and bacteria, tRNA wobble adenine is negatively selected, shrinking the maximum size of the primordial genetic code to 48 anticodons. Because 64 codons are recognized in mRNA, tRNA-mRNA coevolution requires tRNA wobble position ambiguity leading to degeneracy of the code.

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.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into ¹H,³¹P-correlation experiments as well as a BEST-HNP experiment exploiting ^3hJ_N,P rather than ^2hJ_H,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

Intermolecular base stacking mediates RNA-RNA interaction in a crystal structure of the RNA chaperone Hfq

The RNA-chaperone Hfq catalyses the annealing of bacterial small RNAs (sRNAs) with target mRNAs to regulate gene expression in response to environmental stimuli. Hfq acts on a diverse set of sRNA-mRNA pairs using a variety of different molecular mechanisms. Here, we present an unusual crystal structure showing two Hfq-RNA complexes interacting via their bound RNA molecules. The structure contains two Hfq₆:A₁₈ RNA assemblies positioned face-to-face, with the RNA molecules turned towards each other and connected via interdigitating base stacking interactions at the center. Biochemical data further confirm the observed interaction, and indicate that RNA-mediated contacts occur between Hfq-RNA complexes with various (ARN)_X motif containing RNA sequences in vitro, including the stress response regulator OxyS and its target, fhlA. A systematic computational survey also shows that phylogenetically conserved (ARN)_X motifs are present in a subset of sRNAs, some of which share similar modular architectures. We hypothesise that Hfq can co-opt RNA-RNA base stacking, an unanticipated structural trick, to promote the interaction of (ARN)_X motif containing sRNAs with target mRNAs on a “speed-dating” fashion, thereby supporting their regulatory function.

tRNA structure and evolution and standardization to the three nucleotide genetic code

Cloverleaf tRNA with a 75 nucleotide (nt) core is posited to have evolved from ligation of three 31 nt minihelices followed by symmetric internal deletions of 9 nt within ligated acceptor stems. Statistical tests strongly support the model. Although the tRNA anticodon loop and T loop are homologs, their U-turns have been treated as distinct motifs. An appropriate comparison, however, shows that intercalation of D loop G19 between T loop bases 4 and 5 causes elevation of T loop base 5 and flipping of T loop bases 6 and 7 out of the 7 nt loop. In the anticodon loop, by contrast, loop bases 3–7 stack tightly to form a stiff connection to mRNA. Furthermore, we identify ancient repeat sequences of 3 (GCG), 5 (UAGCC) and 17 nt (∼CCGGGUUCAAAACCCGG) that comprise 75 out of 75 nts of the tRNA cloverleaf core. To present a sufficiently stiff 3-nt anticodon, a 7-nt anticodon loop was necessary with a U-turn between loop positions 2 and 3. Cloverleaf tRNA, therefore, was a radical evolutionary innovation essential for the 3-nt code. Conservation of GCG and UAGCC repeat sequences indicates that cloverleaf tRNA is at the interface between a strange RNA repeat world and the first evolution of molecules that fold to assume biologic functions. We posit that cloverleaf tRNA was the molecular archetype around which translation systems evolved.

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.

Revisiting GNRA and UNCG folds: U-turns versus Z-turns in RNA hairpin loops

When thinking about RNA three-dimensional structures, coming across GNRA and UNCG tetraloops is perceived as a boon since their folds have been extensively described. Nevertheless, analyzing loop conformations within RNA and RNP structures led us to uncover several instances of GNRA and UNCG loops that do not fold as expected. We noticed that when a GNRA does not assume its "natural" fold, it adopts the one we typically associate with a UNCG sequence. The same folding interconversion may occur for loops with UNCG sequences, for instance within tRNA anticodon loops. Hence, we show that some structured tetranucleotide sequences starting with G or U can adopt either of these folds. The underlying structural basis that defines these two fold types is the mutually exclusive stacking of a backbone oxygen on either the first (in GNRA) or the last nucleobase (in UNCG), generating an oxygen-π contact. We thereby propose to refrain from using sequences to distinguish between loop conformations. Instead, we suggest using descriptors such as U-turn (for "GNRA-type" folds) and a newly described Z-turn (for "UNCG-type" folds). Because tetraloops adopt for the largest part only two (inter)convertible turns, we are better able to interpret from a structural perspective loop interchangeability occurring in ribosomes and viral RNA. In this respect, we propose a general view on the inclination for a given sequence to adopt (or not) a specific fold. We also suggest how long-noncoding RNAs may adopt discrete but transient structures, which are therefore hard to predict.

RNA-based recognition and targeting: sowing the seeds of specificity

RNA is involved in the regulation of multiple cellular processes, often by forming sequence-specific base pairs with cellular RNA or DNA targets that must be identified among the large number of nucleic acids in a cell. Several RNA-based regulatory systems in eukaryotes, bacteria and archaea, including microRNAs (miRNAs), small interfering RNAs (siRNAs), CRISPR RNAs (crRNAs) and small RNAs (sRNAs) that are dependent on the RNA chaperone protein Hfq, achieve specificity using similar strategies. Central to their function is the presentation of short 'seed sequences' within a ribonucleoprotein complex to facilitate the search for and recognition of targets.

Stacking interactions involving non-Watson–Crick basepairs: dispersion corrected density functional theory studies

Stacking interactions between a non Watson–Crick G:A S:HT basepair and C:G basepair is predicted in terms of roll, twist and slide basepair step parameters using DFT-D augmented with coarse-grain energy penalty for sugar–phosphate backbone.

Identification of a GUAAY Pentaloop Sequence Involved in a Novel RNA Loop-Helix Interaction

Large RNAs often utilize GNRA tetraloops as structural elements to stabilize the overall tertiary fold. These tetraloop-receptor (TR) interactions have a conserved geometry in which the tetraloop docks into the receptor at an angle of ~15° from the helix containing the receptor. Here, we show that the conserved GUAAY pentaloop found in domain III of group IIB1 introns participates in a novel class of RNA tertiary interaction with a geometry and mode of binding that are significantly different from that found in GNRA TR interactions. This pentaloop is highly conserved within the IIB1 class and interacts with the minor groove of the catalytic domain V. The base planes of the loop and receptor nucleotides are not coplanar and greatly deviate from standard A-minor motifs. The helical axis of the GUAAY stem loop diverges ~70° from the angle of insertion found in a typical GNRA TR interaction. Therefore, the loop architecture and insertion orientation are distinctive, with in vitro splicing data indicating that a GNRA tetraloop is incompatible at this position. The GUAAY pentaloop-receptor motif is also found in the structure of the eukaryotic thiamine pyrophosphate riboswitch in the context of a hexanucleotide loop sequence. We therefore propose, based on phylogenetic, structural, and biochemical data, that the GUAAY pentaloop-receptor interaction represents a novel structural motif that is present in multiple structured RNAs.

Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding

The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg²⁺ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl₂ were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.

Identification of determinants for tRNA substrate recognition by Escherichia coli C/U34 2'-O-methyltransferase

Post-transcriptional modifications bring chemical diversity to tRNAs, especially at positions 34 and 37 of the anticodon stem-loop (ASL). TrmL is the prokaryotic methyltransferase that catalyzes the transfer of the methyl group from S-adenosyl-L-methionine to the wobble base of tRNA^Leu_CAA and tRNA^Leu_UAA isoacceptors. This Cm34/Um34 modification affects codon-anticodon interactions and is essential for translational fidelity. TrmL-catalyzed 2′-O-methylation requires its homodimerization; however, understanding of the tRNA recognition mechanism by TrmL remains elusive. In the current study, by measuring tRNA methylation by TrmL and performing kinetic analysis of tRNA mutants, we found that TrmL exhibits a fine-tuned tRNA substrate recognition mechanism. Anticodon stem-loop minihelices with an extension of 2 base pairs are the minimal substrate for EcTrmL methylation. A35 is a key residue for TrmL recognition, while A36-A37-A38 are important either via direct interaction with TrmL or due to the necessity for prior isopentenylation (i⁶) at A37. In addition, TrmL only methylates pyrimidines but not purine residues at the wobble position, and the 2′-O-methylation relies on prior N⁶-isopentenyladenosine modification at position 37.

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.

Folding of RNA tertiary structure: Linkages between backbone phosphates, ions, and water

The functional forms of many RNAs have compact architectures. The placement of phosphates within such structures must be influenced not only by the strong electrostatic repulsion between phosphates, but also by networks of interactions between phosphates, water, and mobile ions. This review first explores what has been learned of the basic thermodynamic constraints on these arrangements from studies of hydration and ions in simple DNA molecules, and then gives an overview of what is known about ion and water interactions with RNA structures. A brief survey of RNA crystal structures identifies several interesting architectures in which closely spaced phosphates share hydration shells or phosphates are buried in environments that provide intramolecular hydrogen bonds or site-bound cations. Formation of these structures must require strong coupling between the uptake of ions and release of water.

tRNA-targeting ribonucleases: molecular mechanisms and insights into their physiological roles

Most bacteria produce antibacterial proteins known as bacteriocins, which aid bacterial defence systems to provide a physiological advantage. To date, many kinds of bacteriocins have been characterized. Colicin has long been known as a plasmidborne bacteriocin that kills other Escherichia coli cells lacking the same plasmid. To defeat other cells, colicins exert specific activities such as ion-channel, DNase, and RNase activity. Colicin E5 and colicin D impair protein synthesis in sensitive E. coli cells; however, their physiological targets have not long been identified. This review describes our finding that colicins E5 and D are novel RNases targeting specific E. coli tRNAs and elucidates their enzymatic properties based on biochemical analyses and X-ray crystal structures. Moreover, tRNA cleavage mediates bacteriostasis, which depends on trans-translation. Based on these results and others, cell growth regulation depending on tRNA cleavage is also discussed.

The role of sequence in altering the unfolding pathway of an RNA pseudoknot: a steered molecular dynamics study

This work highlights a sequence dependent unfolding pathway of an RNA pseudoknot under force-induced pulling conditions.

An extended dsRBD is required for post-transcriptional modification in human tRNAs

In tRNA, dihydrouridine is a conserved modified base generated by the post-transcriptional reduction of uridine. Formation of dihydrouridine 20, located in the D-loop, is catalyzed by dihydrouridine synthase 2 (Dus2). Human Dus2 (HsDus2) expression is upregulated in lung cancers, offering a growth advantage throughout its ability to interact with components of the translation apparatus and inhibit apoptosis. Here, we report the crystal structure of the individual domains of HsDus2 and their functional characterization. HsDus2 is organized into three major modules. The N-terminal catalytic domain contains the flavin cofactor involved in the reduction of uridine. The second module is the conserved α-helical domain known as the tRNA binding domain in HsDus2 homologues. It is connected via a flexible linker to an unusual extended version of a dsRNA binding domain (dsRBD). Enzymatic assays and yeast complementation showed that the catalytic domain binds selectively NADPH but cannot reduce uridine in the absence of the dsRBD. While in Dus enzymes from bacteria, plants and fungi, tRNA binding is essentially achieved by the α-helical domain, we showed that in HsDus2 this function is carried out by the dsRBD. This is the first reported case of a tRNA-modifying enzyme carrying a dsRBD used to bind tRNAs.