Stabilization of the anticodon stem-loop of tRNALys,3 by an A+-C base-pair and by pseudouridine

Adjusting the Energy Profile for CH-O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

The molecular basis of tRNA selectivity by human pseudouridine synthase 3

Pseudouridine (Ψ), the isomer of uridine, is ubiquitously found in RNA, including tRNA, rRNA, and mRNA. Human pseudouridine synthase 3 (PUS3) catalyzes pseudouridylation of position 38/39 in tRNAs. However, the molecular mechanisms by which it recognizes its RNA targets and achieves site specificity remain elusive. Here, we determine single-particle cryo-EM structures of PUS3 in its apo form and bound to three tRNAs, showing how the symmetric PUS3 homodimer recognizes tRNAs and positions the target uridine next to its active site. Structure-guided and patient-derived mutations validate our structural findings in complementary biochemical assays. Furthermore, we deleted PUS1 and PUS3 in HEK293 cells and mapped transcriptome-wide Ψ sites by Pseudo-seq. Although PUS1-dependent sites were detectable in tRNA and mRNA, we found no evidence that human PUS3 modifies mRNAs. Our work provides the molecular basis for PUS3-mediated tRNA modification in humans and explains how its tRNA modification activity is linked to intellectual disabilities.

Decoding the 'Fifth' Nucleotide: Impact of RNA Pseudouridylation on Gene Expression and Human Disease

Cellular RNAs, both coding and noncoding are adorned by > 100 chemical modifications, which impact various facets of RNA metabolism and gene expression. Very often derailments in these modifications are associated with a plethora of human diseases. One of the most oldest of such modification is pseudouridylation of RNA, wherein uridine is converted to a pseudouridine (Ψ) via an isomerization reaction. When discovered, Ψ was referred to as the 'fifth nucleotide' and is chemically distinct from uridine and any other known nucleotides. Experimental evidence accumulated over the past six decades, coupled together with the recent technological advances in pseudouridine detection, suggest the presence of pseudouridine on messenger RNA, as well as on diverse classes of non-coding RNA in human cells. RNA pseudouridylation has widespread effects on cellular RNA metabolism and gene expression, primarily via stabilizing RNA conformations and destabilizing interactions with RNA-binding proteins. However, much remains to be understood about the RNA targets and their recognition by the pseudouridylation machinery, the regulation of RNA pseudouridylation, and its crosstalk with other RNA modifications and gene regulatory processes. In this review, we summarize the mechanism and molecular machinery involved in depositing pseudouridine on target RNAs, molecular functions of RNA pseudouridylation, tools to detect pseudouridines, the role of RNA pseudouridylation in human diseases like cancer, and finally, the potential of pseudouridine to serve as a biomarker and as an attractive therapeutic target.

Efficient Synthetic Access to Stable Isotope Labelled Pseudouridine Phosphoramidites for RNA NMR Spectroscopy

Here we report the efficient synthetic access to 13C/15N-labelled pseudouridine phosphoramidites, which were incorporated into a binary H/ACA box guide RNA/product complex comprising 77 nucleotides (nts) in total and into a 75 nt E. coli tRNAGly. The stable isotope (SI) labelled pseudouridines were produced via a highly efficient chemo-enzymatic synthesis. 13C/15N labelled uracils were produced via chemical synthesis and enzymatically converted to pseudouridine 5'-monophosphate (ΨMP) by using YeiN, a Ψ-5'-monophosphate C-glycosidase. Removal of the 5'-phosphate group yielded the desired pseudouridine nucleoside (Ψ), which was transformed into a phosphoramidite building suitable for RNA solid phase synthesis. A Ψ -building block carrying both a 13C and a 15N label was incorporated into a product RNA and the complex formation with a 63 nt H/ACA box RNA could be observed via NMR. Furthermore, the SI labelled pseudouridine building block was used to determine imino proton bulk water exchange rates of a 75 nt E. coli tRNAGly CCmnm5U, identifying the TΨC-loop 5-methyluridine as a modifier of the exchange rates. The efficient synthetic access to SI-labelled Ψ building blocks will allow the solution and solid-state NMR spectroscopic studies of Ψ containing RNAs and will facilitate the mass spectrometric analysis of Ψ-modified nucleic acids.

tRNA Modifications and Dysregulation: Implications for Brain Diseases

Transfer RNAs (tRNAs) are well-known for their essential function in protein synthesis. Recent research has revealed a diverse range of chemical modifications that tRNAs undergo, which are crucial for various cellular processes. These modifications are necessary for the precise and efficient translation of proteins and also play important roles in gene expression regulation and cellular stress response. This review examines the role of tRNA modifications and dysregulation in the pathophysiology of various brain diseases, including epilepsy, stroke, neurodevelopmental disorders, brain tumors, Alzheimer's disease, and Parkinson's disease. Through a comprehensive analysis of existing research, our study aims to elucidate the intricate relationship between tRNA dysregulation and brain diseases. This underscores the critical need for ongoing exploration in this field and provides valuable insights that could facilitate the development of innovative diagnostic tools and therapeutic approaches, ultimately improving outcomes for individuals grappling with complex neurological conditions.

Pseudouridine and N1-methylpseudouridine as potent nucleotide analogues for RNA therapy and vaccine development

Modified nucleosides are integral to modern drug development, serving as crucial building blocks for creating safer, more potent, and more precisely targeted therapeutic interventions. Nucleobase modifications often confer antiviral and anti-cancer activity as monomers. When incorporated into nucleic acid oligomers, they increase stability against degradation by enzymes, enhancing the drugs' lifespan within the body. Moreover, modification strategies can mitigate potential toxic effects and reduce immunogenicity, making drugs safer and better tolerated. Particularly, N1-methylpseudouridine modification improved the efficacy of the mRNA coding for spike protein of COVID-19. This became a crucial step for developing COVID-19 vaccine applied during the 2020 pandemic. This makes N1-methylpseudouridine, and its "parent" analogue pseudouridine, potent nucleotide analogues for future RNA therapy and vaccine development. This review focuses on the structure and properties of pseudouridine and N1-methylpseudouridine. RNA has a greater structural versatility, different conformation, and chemical reactivity than DNA. Watson-Crick pairing is not strictly followed by RNA that has more unusual base pairs and base-triplets. This requires detailed structural studies and structure-activity relationship analyses for RNA, also when modifications are incorporated. Recent successes in this direction are revised in this review. We describe recent successes with using pseudouridine and N1-methylpseudouridine in mRNA drug candidates. We also highlight remaining challenges that need to be solved to develop new mRNA vaccines and therapies.

Diversity in Biological Function and Mechanism of the tRNA Methyltransferase Trm10

Transfer ribonucleic acid (tRNA) is the most highly modified RNA species in the cell, and loss of tRNA modifications can lead to growth defects in yeast as well as metabolic, neurological, and mitochondrial disorders in humans. Significant progress has been made toward identifying the enzymes that are responsible for installing diverse modifications in tRNA, revealing a landscape of fascinating biological and mechanistic diversity that remains to be fully explored. Most early discoveries of tRNA modification enzymes were in model systems, where many enzymes were not strictly required for viability, an observation somewhat at odds with the extreme conservation of many of the same enzymes throughout multiple domains of life. Moreover, many tRNA modification enzymes act on more than one type of tRNA substrate, which is not necessarily surprising given the similar overall secondary and tertiary structures of tRNA, yet biochemical characterization has revealed interesting patterns of substrate specificity that can be challenging to rationalize on a molecular level. Questions about how many enzymes efficiently select a precise set of target tRNAs from among a structurally similar pool of molecules persist.The tRNA methyltransferase Trm10 provides an exciting paradigm to study the biological and mechanistic questions surrounding tRNA modifications. Even though the enzyme was originally characterized in Saccharomyces cerevisiae where its deletion causes no detectable phenotype under standard lab conditions, several more recently identified phenotypes provide insight into the requirement for this modification in the overall quality control of the tRNA pool. Studies of Trm10 in yeast also revealed another characteristic feature that has turned out to be a conserved feature of enzymes throughout the Trm10 family tree. We were initially surprised to see that purified S. cerevisiae Trm10 was capable of modifying tRNA substrates that were not detectably modified by the enzyme in vivo in yeast. This pattern has continued to emerge as we and others have studied Trm10 orthologs from Archaea and Eukarya, with enzymes exhibiting in vitro substrate specificities that can differ significantly from in vivo patterns of modification. While this feature complicates efforts to predict substrate specificities of Trm10 enzymes in the absence of appropriate genetic systems, it also provides an exciting opportunity for studying how enzyme activities can be regulated to achieve dynamic patterns of biological tRNA modification, which have been shown to be increasingly important for stress responses and human disease. Finally, the intriguing diversity in target nucleotide modification that has been revealed among Trm10 orthologs is distinctive among known tRNA modifying enzymes and necessitates unusual and likely novel catalytic strategies for methylation that are being revealed by biochemical and structural studies directed toward various family members. These efforts will no doubt yield more surprising discoveries in terms of tRNA modification enzymology.

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.

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.

The Repertoire of RNA Modifications Orchestrates a Plethora of Cellular Responses

Although a plethora of DNA modifications have been extensively investigated in the last decade, recent breakthroughs in molecular biology, including high throughput sequencing techniques, have enabled the identification of post-transcriptional marks that decorate RNAs; hence, epitranscriptomics has arisen. This recent scientific field aims to decode the regulatory layer of the transcriptome and set the ground for the detection of modifications in ribose nucleotides. Until now, more than 170 RNA modifications have been reported in diverse types of RNA that contribute to various biological processes, such as RNA biogenesis, stability, and transcriptional and translational accuracy. However, dysfunctions in the RNA-modifying enzymes that regulate their dynamic level can lead to human diseases and cancer. The present review aims to highlight the epitranscriptomic landscape in human RNAs and match the catalytic proteins with the deposition or deletion of a specific mark. In the current review, the most abundant RNA modifications, such as N6-methyladenosine (m⁶A), N5-methylcytosine (m⁵C), pseudouridine (Ψ) and inosine (I), are thoroughly described, their functional and regulatory roles are discussed and their contributions to cellular homeostasis are stated. Ultimately, the involvement of the RNA modifications and their writers, erasers, and readers in human diseases and cancer is also discussed.

Destabilization of mutated human PUS3 protein causes intellectual disability

Pseudouridine (Ψ) is an RNA base modification ubiquitously found in many types of RNAs. In humans, the isomerization of uridine is catalyzed by different stand-alone pseudouridine synthases (PUS). Genomic mutations in the human pseudouridine synthase 3 gene (PUS3) have been identified in patients with neurodevelopmental disorders. However, the underlying molecular mechanisms that cause the disease phenotypes remain elusive. Here, we utilize exome sequencing to identify genomic variants that lead to a homozygous amino acid substitution (p.[(Tyr71Cys)];[(Tyr71Cys)]) in human PUS3 of two affected individuals and a compound heterozygous substitution (p.[(Tyr71Cys)];[(Ile299Thr)]) in a third patient. We obtain wild-type and mutated full-length human recombinant PUS3 proteins and characterize the enzymatic activity in vitro. Unexpectedly, we find that the p.Tyr71Cys substitution neither affect tRNA binding nor pseudouridylation activity in vitro, but strongly impair the thermostability profile of PUS3, while the p.Ile299Thr mutation causes protein aggregation. Concomitantly, we observe that the PUS3 protein levels as well as the level of PUS3-dependent Ψ levels are strongly reduced in fibroblasts derived from all three patients. In summary, our results directly illustrate the link between the identified PUS3 variants and reduced Ψ levels in the patient cells, providing a molecular explanation for the observed clinical phenotypes.

Alternative RNA Conformations: Companion or Combatant

RNA molecules, in one form or another, are involved in almost all aspects of cell physiology, as well as in disease development. The diversity of the functional roles of RNA comes from its intrinsic ability to adopt complex secondary and tertiary structures, rivaling the diversity of proteins. The RNA molecules form dynamic ensembles of many interconverting conformations at a timescale of seconds, which is a key for understanding how they execute their cellular functions. Given the crucial role of RNAs in various cellular processes, we need to understand the RNA molecules from a structural perspective. Central to this review are studies aimed at revealing the regulatory role of conformational equilibria in RNA in humans to understand genetic diseases such as cancer and neurodegenerative diseases, as well as in pathogens such as bacteria and viruses so as to understand the progression of infectious diseases. Furthermore, we also summarize the prior studies on the use of RNA structures as platforms for the rational design of small molecules for therapeutic applications.

Role of main RNA modifications in cancer: N6-methyladenosine, 5-methylcytosine, and pseudouridine

Cancer is one of the major diseases threatening human life and health worldwide. Epigenetic modification refers to heritable changes in the genetic material without any changes in the nucleic acid sequence and results in heritable phenotypic changes. Epigenetic modifications regulate many biological processes, such as growth, aging, and various diseases, including cancer. With the advancement of next-generation sequencing technology, the role of RNA modifications in cancer progression has become increasingly prominent and is a hot spot in scientific research. This review studied several common RNA modifications, such as N⁶-methyladenosine, 5-methylcytosine, and pseudouridine. The deposition and roles of these modifications in coding and noncoding RNAs are summarized in detail. Based on the RNA modification background, this review summarized the expression, function, and underlying molecular mechanism of these modifications and their regulators in cancer and further discussed the role of some existing small-molecule inhibitors. More in-depth studies on RNA modification and cancer are needed to broaden the understanding of epigenetics and cancer diagnosis, treatment, and prognosis.

Data-informed reparameterization of modified RNA and the effect of explicit water models: application to pseudouridine and derivatives

Pseudouridine is one of the most abundant post-transcriptional modifications in RNA. We have previously shown that the FF99-derived parameters for pseudouridine and some of its naturally occurring derivatives in the AMBER distribution either alone or in combination with the revised γ torsion parameters (parmbsc0) failed to reproduce their conformational characteristics observed experimentally (Deb et al. in J Chem Inf Model 54:1129–1142, 2014; Deb et al. in J Comput Chem 37:1576–1588, 2016; Dutta et al. in J Chem Inf Model 60:4995–5002, 2020). However, the application of the recommended bsc0 correction did lead to an improvement in the description not only of the distribution in the γ torsional space but also of the sugar pucker distributions. In an earlier study, we examined the transferability of the revised glycosidic torsion parameters (χ_IDRP) for Ψ to its derivatives. We noticed that although these parameters in combination with the AMBER FF99-derived parameters and the revised γ torsional parameters resulted in conformational properties of these residues that were in better agreement with experimental observations, the sugar pucker distributions were still not reproduced accurately. Here we report a new set of partial atomic charges for pseudouridine, 1-methylpseudouridine, 3-methylpseudouridine and 2′-O-methylpseudouridine and a new set of glycosidic torsional parameters (χ_ND) based on chosen glycosidic torsional profiles that most closely corresponded to the NMR data for conformational propensities and studied their effect on the conformational distributions using REMD simulations at the individual nucleoside level. We have also studied the effect of the choice of water model on the conformational characteristics of these modified nucleosides. Our observations suggest that the current revised set of parameters and partial atomic charges describe the sugar pucker distributions for these residues more accurately and that the choice of a suitable water model is important for the accurate description of their conformational properties. We have further validated the revised sets of parameters by studying the effect of substitution of uridine with pseudouridine within single stranded RNA oligonucleotides on their conformational and hydration characteristics.

Increased Photostability of the Integral mRNA Vaccine Component N1 -Methylpseudouridine Compared to Uridine

N1 -Methylation of pseudouridine (m1 ψ) replaces uridine (Urd) in several therapeutics, including the Moderna and BioNTech-Pfizer COVID-19 vaccines. Importantly, however, it is currently unknown if exposure to electromagnetic radiation can affect the chemical integrity and intrinsic stability of m1 ψ. In this study, the photochemistry of m1 ψ is compared to that of uridine by using photoirradiation at 267 nm, steady-state spectroscopy, and quantum-chemical calculations. Furthermore, femtosecond transient absorption measurements are collected to delineate the electronic relaxation mechanisms for both nucleosides under physiologically relevant conditions. It is shown that m1 ψ exhibits a 12-fold longer 1 ππ* decay lifetime than uridine and a 5-fold higher fluorescence yield. Notably, however, the experimental results also demonstrate that most of the excited state population in both molecules decays back to the ground state in an ultrafast time scale and that m1 ψ is 6.7-fold more photostable than Urd following irradiation at 267 nm.

Characterizing RNA Pseudouridylation by Convolutional Neural Networks

Pseudouridine (Ψ) is the most prevalent post-transcriptional RNA modification and is widespread in small cellular RNAs and mRNAs. However, the functions, mechanisms, and precise distribution of Ψs (especially in mRNAs) still remain largely unclear. The landscape of Ψs across the transcriptome has not yet been fully delineated. Here, we present a highly effective model based on a convolutional neural network (CNN), called PseudoUridyLation Site Estimator (PULSE), to analyze large-scale profiling data of Ψ sites and characterize the contextual sequence features of pseudouridylation. PULSE, consisting of two alternatively-stacked convolution and pooling layers followed by a fully-connected neural network, can automatically learn the hidden patterns of pseudouridylation from the local sequence information. Extensive validation tests demonstrated that PULSE can outperform other state-of-the-art prediction methods and achieve high prediction accuracy, thus enabling us to further characterize the transcriptome-wide landscape of Ψ sites. We further showed that the prediction results derived from PULSE can provide novel insights into understanding the functional roles of pseudouridylation, such as the regulations of RNA secondary structure, codon usage, translation, and RNA stability, and the connection to single nucleotide variants. The source code and final model for PULSE are available at https://github.com/mlcb-thu/PULSE.

Pseudouridines in RNAs: switching atoms means shifting paradigms

The structure, stability, and function of various coding and noncoding RNAs are influenced by chemical modifications. Pseudouridine (Ψ) is one of the most abundant post‐transcriptional RNA base modifications and has been detected at individual positions in tRNAs, rRNAs, mRNAs, and snRNAs, which are referred to as Ψ‐sites. By allowing formation of additional bonds with neighboring atoms, Ψ strengthens RNA–RNA and RNA–protein interactions. Although many aspects of the underlying modification reactions remain unclear, the advent of new transcriptome‐wide methods to quantitatively detect Ψ‐sites has recently changed our perception of the functional roles and importance of Ψ. For instance, it is now clear that the occurrence of Ψs appears to be directly linked to the lifetime and the translation efficiency of a given mRNA molecule. Furthermore, the administration of Ψ‐containing RNAs reduces innate immune responses, which appears strikingly advantageous for the development of generations of mRNA‐based vaccines. In this review, we aim to comprehensively summarize recent discoveries that highlight the impact of Ψ on various types of RNAs and outline possible novel biomedical applications of Ψ.

Transfer RNAs: diversity in form and function

As the adaptor that decodes mRNA sequence into protein, the basic aspects of tRNA structure and function are central to all studies of biology. Yet the complexities of their properties and cellular roles go beyond the view of tRNAs as static participants in protein synthesis. Detailed analyses through more than 60 years of study have revealed tRNAs to be a fascinatingly diverse group of molecules in form and function, impacting cell biology, physiology, disease and synthetic biology. This review analyzes tRNA structure, biosynthesis and function, and includes topics that demonstrate their diversity and growing importance.

Revisiting tRNA chaperones: New players in an ancient game

tRNAs undergo an extensive maturation process including post-transcriptional modifications that influence secondary and tertiary interactions. Precursor and mature tRNAs lacking key modifications are often recognized as aberrant and subsequently targeted for decay, illustrating the importance of modifications in promoting structural integrity. tRNAs also rely on tRNA chaperones to promote the folding of misfolded substrates into functional conformations. The best characterized tRNA chaperone is the La protein, which interacts with nascent RNA polymerase III transcripts to promote folding and offers protection from exonucleases. More recently, certain tRNA modification enzymes have also been demonstrated to possess tRNA folding activity distinct from their catalytic activity, suggesting that they may act as tRNA chaperones. In this review, we will discuss pioneering studies relating post-transcriptional modification to tRNA stability and decay pathways, present recent advances into the mechanism by which the RNA chaperone La assists pre-tRNA maturation, and summarize emerging research directions aimed at characterizing modification enzymes as tRNA chaperones. Together, these findings shed light on the importance of tRNA folding and how tRNA chaperones, in particular, increase the fraction of nascent pre-tRNAs that adopt a folded, functional conformation.

Pseudouridines of tRNA Anticodon Stem-Loop Have Unexpected Role in Mutagenesis in Pseudomonas sp

Pseudouridines are known to be important for optimal translation. In this study we demonstrate an unexpected link between pseudouridylation of tRNA and mutation frequency in Pseudomonas species. We observed that the lack of pseudouridylation activity of pseudouridine synthases TruA or RluA elevates the mutation frequency in Pseudomonas putida 3 to 5-fold. The absence of TruA but not RluA elevates mutation frequency also in Pseudomonas aeruginosa. Based on the results of genetic studies and analysis of proteome data, the mutagenic effect of the pseudouridylation deficiency cannot be ascribed to the involvement of error-prone DNA polymerases or malfunctioning of DNA repair pathways. In addition, although the deficiency in TruA-dependent pseudouridylation made P. putida cells more sensitive to antimicrobial compounds that may cause oxidative stress and DNA damage, cultivation of bacteria in the presence of reactive oxygen species (ROS)-scavenging compounds did not eliminate the mutator phenotype. Thus, the elevated mutation frequency in the absence of tRNA pseudouridylation could be the result of a more specific response or, alternatively, of a cumulative effect of several small effects disturbing distinct cellular functions, which remain undetected when studied independently. This work suggests that pseudouridines link the translation machinery to mutation frequency.

Titration of Adenine in a GA Mismatch with Grand Canonical Simulations

Computational prediction of the pKa of ionizable groups remains a central challenge in biomolecular modeling. Although all-atom fixed-charge force fields could be accurate to describe the interaction network within the biomolecules, proper sampling techniques are required to obtain the thermodynamic information in the (de)protonation event. Sufficient sampling requires an ensemble of structures from simulations with proper treatments of the acid-base equilibria, and the grand canonical simulation technique could be used to model the growth/annihilation of hydrogen atoms by merging Hamiltonians of different protonation states into one simulation ensemble. The electrostatic feature of nucleotide systems is especially difficult to model, and the situation becomes more challenging when the ionizable site is highly perturbed. Although there are many successful predictions obtained from the grand canonical constant pH simulations, few reports focus on highly perturbed nucleotide systems with unconventional base-pair features. In this work, with the discrete constant pH method, we investigate the titration thermodynamics of an adenine in the catalytic triad in a 35-nucleotide single-stranded RNA hairpin, featuring an unconventional GA mismatch and a substantially shifted pKa value. Validation tests are performed with two system setups, both of which provide pKa predictions in good agreement with the experimental value. A single-configuration-based technique is used to calculate the pKa for comparison. The current success indicates the predictive power of the current nucleotide modeling framework.

An NMR-based approach reveals the core structure of the functional domain of SINEUP lncRNAs

Long non-coding RNAs (lncRNAs) are attracting widespread attention for their emerging regulatory, transcriptional, epigenetic, structural and various other functions. Comprehensive transcriptome analysis has revealed that retrotransposon elements (REs) are transcribed and enriched in lncRNA sequences. However, the functions of lncRNAs and the molecular roles of the embedded REs are largely unknown. The secondary and tertiary structures of lncRNAs and their embedded REs are likely to have essential functional roles, but experimental determination and reliable computational prediction of large RNA structures have been extremely challenging. We report here the nuclear magnetic resonance (NMR)-based secondary structure determination of the 167-nt inverted short interspersed nuclear element (SINE) B2, which is embedded in antisense Uchl1 lncRNA and upregulates the translation of sense Uchl1 mRNAs. By using NMR 'fingerprints' as a sensitive probe in the domain survey, we successfully divided the full-length inverted SINE B2 into minimal units made of two discrete structured domains and one dynamic domain without altering their original structures after careful boundary adjustments. This approach allowed us to identify a structured domain in nucleotides 31-119 of the inverted SINE B2. This approach will be applicable to determining the structures of other regulatory lncRNAs.

Regulation and Function of RNA Pseudouridylation in Human Cells

Recent advances in pseudouridine detection reveal a complex pseudouridine landscape that includes messenger RNA and diverse classes of noncoding RNA in human cells. The known molecular functions of pseudouridine, which include stabilizing RNA conformations and destabilizing interactions with varied RNA-binding proteins, suggest that RNA pseudouridylation could have widespread effects on RNA metabolism and gene expression. Here, we emphasize how much remains to be learned about the RNA targets of human pseudouridine synthases, their basis for recognizing distinct RNA sequences, and the mechanisms responsible for regulated RNA pseudouridylation. We also examine the roles of noncoding RNA pseudouridylation in splicing and translation and point out the potential effects of mRNA pseudouridylation on protein production, including in the context of therapeutic mRNAs.

Extracurricular Functions of tRNA Modifications in Microorganisms

Transfer RNAs (tRNAs) are essential adaptors that mediate translation of the genetic code. These molecules undergo a variety of post-transcriptional modifications, which expand their chemical reactivity while influencing their structure, stability, and functionality. Chemical modifications to tRNA ensure translational competency and promote cellular viability. Hence, the placement and prevalence of tRNA modifications affects the efficiency of aminoacyl tRNA synthetase (aaRS) reactions, interactions with the ribosome, and transient pairing with messenger RNA (mRNA). The synthesis and abundance of tRNA modifications respond directly and indirectly to a range of environmental and nutritional factors involved in the maintenance of metabolic homeostasis. The dynamic landscape of the tRNA epitranscriptome suggests a role for tRNA modifications as markers of cellular status and regulators of translational capacity. This review discusses the non-canonical roles that tRNA modifications play in central metabolic processes and how their levels are modulated in response to a range of cellular demands.

Adenine·cytosine substitutions are an alternative pathway of compensatory mutation in angiosperm ITS2

Compensatory mutations are crucial for functional RNA because they maintain RNA configuration and thus function. Compensatory mutation has traditionally been considered to be a two-step substitution through the GU-base-pair intermediate. We tested for an alternative AC-mediated compensatory mutation (ACCM). We investigated ACCMs by using a comprehensive sampling of ribosomal internal transcribed spacer 2 (ITS2) from 3934 angiosperm species in 80 genera and 55 families. We predicted ITS2 consensus secondary structures by using LocARNA for structure-based alignment and partitioning paired and unpaired regions. We examined and compared the substitution rates and frequencies among base pairs by using RNA-specific models. Base-pair states of ACCMs were mapped onto the inferred phylogenetic trees to infer their evolution. All types of compensatory mutations involving the AC intermediate were observed, but the most frequent substitutions were with AU or GC pairs, which are part of the AU-AC-GC pathway. Compared with the GU intermediate, AC had a lower frequency and higher mutability. Within the AU-AC-GC pathway, the AU-AC substitution rate was much slower than the AC-GC substitution rate. No consistently higher overall rate was identified for either pathway among all 80 sampled lineages, though compensatory mutations through the AC intermediate averaged about half that through the GU intermediate. These results demonstrate an alternative compensatory mutation between AU and GC that helps address the controversial inference of inferred simultaneous double substitutions.

Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair

Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5′-GΨC/3′-CAG and 5′-CΨG/3′-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5′-GΨC/3′-CAG, 5′-CΨG/3′-GAC, 5′-AΨU/3′-UAA and 5′-UΨA/3′-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.

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.

Importance of a tRNA anticodon loop modification and a conserved, noncanonical anticodon stem pairing in tRNACGGProfor decoding

Modification of anticodon nucleotides allows tRNAs to decode multiple codons, expanding the genetic code. Additionally, modifications located in the anticodon loop, outside the anticodon itself, stabilize tRNA–codon interactions, increasing decoding fidelity. Anticodon loop nucleotide 37 is 3′ to the anticodon and, in tRNA CGG Pro , is methylated at the N1 position in its nucleobase (m1G37). The m1G37 modification in tRNA CGG Pro stabilizes its interaction with the codon and maintains the mRNA frame. However, it is unclear how m1G37 affects binding at the decoding center to both cognate and +1 slippery codons. Here, we show that the tRNA CGG Pro m1G37 modification is important for the association step during binding to a cognate CCG codon. In contrast, m1G37 prevented association with a slippery CCC-U or +1 codon. Similar analyses of frameshift suppressor tRNASufA6, a tRNA CGG Pro derivative containing an extra nucleotide in its anticodon loop that undergoes +1 frameshifting, reveal that m1G37 destabilizes interactions with both the cognate CCG and slippery codons. One reason for this destabilization is the disruption of a conserved U32·A38 nucleotide pairing in the anticodon stem through insertion of G37.5. Restoring the tRNASufA6 U32·A37.5 pairing results in a high-affinity association on the slippery CCC-U codon. Further, an X-ray crystal structure of the 70S ribosome bound to tRNASufA6 U32·A37.5 at 3.6 Å resolution shows a reordering of the anticodon loop consistent with the findings from the high-affinity measurements. Our results reveal how the tRNA modification at nucleotide 37 stabilizes interactions with the mRNA codon to preserve the mRNA frame.

The emerging impact of tRNA modifications in the brain and nervous system

A remarkable number of neurodevelopmental disorders have been linked to defects in tRNA modifications. These discoveries place tRNA modifications in the spotlight as critical modulators of gene expression pathways that are required for proper organismal growth and development. Here, we discuss the emerging molecular and cellular functions of the diverse tRNA modifications linked to cognitive and neurological disorders. In particular, we describe how the structure and location of a tRNA modification influences tRNA folding, stability, and function. We then highlight how modifications in tRNA can impact multiple aspects of protein translation that are instrumental for maintaining proper cellular proteostasis. Importantly, we describe how perturbations in tRNA modification lead to a spectrum of deleterious biological outcomes that can disturb neurodevelopment and neurological function. Finally, we summarize the biological themes shared by the different tRNA modifications linked to cognitive disorders and offer insight into the future questions that remain to decipher the role of tRNA modifications. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Soller Matthias and Dr. Fray Rupert.

Contribution of protein Gar1 to the RNA-guided and RNA-independent rRNA:Ψ-synthase activities of the archaeal Cbf5 protein

Archaeal RNA:pseudouridine-synthase (PUS) Cbf5 in complex with proteins L7Ae, Nop10 and Gar1, and guide box H/ACA sRNAs forms ribonucleoprotein (RNP) catalysts that insure the conversion of uridines into pseudouridines (Ψs) in ribosomal RNAs (rRNAs). Nonetheless, in the absence of guide RNA, Cbf5 catalyzes the in vitro formation of Ψ2603 in Pyrococcus abyssi 23S rRNA and of Ψ55 in tRNAs. Using gene-disrupted strains of the hyperthermophilic archaeon Thermococcus kodakarensis, we studied the in vivo contribution of proteins Nop10 and Gar1 to the dual RNA guide-dependent and RNA-independent activities of Cbf5 on 23S rRNA. The single-null mutants of the cbf5, nop10, and gar1 genes are viable, but display a thermosensitive slow growth phenotype. We also generated a single-null mutant of the gene encoding Pus10, which has redundant activity with Cbf5 for in vitro formation of Ψ55 in tRNA. Analysis of the presence of Ψs within the rRNA peptidyl transferase center (PTC) of the mutants demonstrated that Cbf5 but not Pus10 is required for rRNA modification. Our data reveal that, in contrast to Nop10, Gar1 is crucial for in vivo and in vitro RNA guide-independent formation of Ψ2607 (Ψ2603 in P. abyssi) by Cbf5. Furthermore, our data indicate that pseudouridylation at orphan position 2589 (2585 in P. abyssi), for which no PUS or guide sRNA has been identified so far, relies on RNA- and Gar1-dependent activity of Cbf5.

Single-Molecule Mechanical Folding and Unfolding of RNA Hairpins: Effects of Single A-U to A·C Pair Substitutions and Single Proton Binding and Implications for mRNA Structure-Induced -1 Ribosomal Frameshifting

A wobble A·C pair can be protonated at near physiological pH to form a more stable wobble A⁺·C pair. Here, we constructed an RNA hairpin (rHP) and three mutants with one A-U base pair substituted with an A·C mismatch on the top (near the loop, U22C), middle (U25C), and bottom (U29C) positions of the stem, respectively. Our results on single-molecule mechanical (un)folding using optical tweezers reveal the destabilization effect of A-U to A·C pair substitution and protonation-dependent enhancement of mechanical stability facilitated through an increased folding rate, or decreased unfolding rate, or both. Our data show that protonation may occur rapidly upon the formation of an apparent mechanical folding transition state. Furthermore, we measured the bulk -1 ribosomal frameshifting efficiencies of the hairpins by a cell-free translation assay. For the mRNA hairpins studied, -1 frameshifting efficiency correlates with mechanical unfolding force at equilibrium and folding rate at around 15 pN. U29C has a frameshifting efficiency similar to that of rHP (∼2%). Accordingly, the bottom 2-4 base pairs of U29C may not form under a stretching force at pH 7.3, which is consistent with the fact that the bottom base pairs of the hairpins may be disrupted by ribosome at the slippery site. U22C and U25C have a similar frameshifting efficiency (∼1%), indicating that both unfolding and folding rates of an mRNA hairpin in a crowded environment may affect frameshifting. Our data indicate that mechanical (un)folding of RNA hairpins may mimic how mRNAs unfold and fold in the presence of translating ribosomes.

Detection of protonated non-Watson-Crick base pairs using electrospray ionization mass spectrometry

Many studies have shown that protonated nucleic acid base pairs are involved in a wide variety of nucleic acid structures. However, little information is available on relative stability of hemiprotonated self- and non-self-dimers at monomer level. We used electrospray ionization mass spectrometry (ESI-MS) to evaluate the relative stability under various concentrations of hydrogen ion. These enable conjecture of the formation of protonated non-Watson-Crick base pairs based on DNA and RNA base sequence. In the present study, we observed that ESI-MS peaks corresponded to respective self-dimers for all examined nucleosides except for adenosine. Peak heights depended on the concentration of hydrogen ion. The ESI-MS peak heights of the hemiprotonated cytidine dimers and the hemiprotonated thymidine dimer sharply increased with increased concentration of hydrogen ion, suggesting direct participation of hydrogen ion in dimer formations. In ESI-MS measurements of the solutions containing adenosine, cytidine, thymidine and guanosine, we observed protonated cytidine-guanosine dimer (CH+-G) and protonated cytidine-thymidine dimer (CH+-T) in addition to hemiprotonated cytidine-cytidine dimer (CH+-C) with following relative peak height, (CH+-C) > (CH+-G) ≈ (CH+-T) > (CH+-A). Additionally, in the ESI-MS measurements of solutions containing adenosine, thymidine and guanosine, we observed a considerable amount of protonated adenosine-guanosine (AH+-G) and protonated adenosine-thymidine (AH+-T).

A novel form of RNA double helix based on G·U and C·A+ wobble base pairing

Wobble base pairs are critical in various physiological functions and have been linked to local structural perturbations in double-helical structures of nucleic acids. We report a 1.38-Å resolution crystal structure of an antiparallel octadecamer RNA double helix in overall A conformation, which includes a unique, central stretch of six consecutive wobble base pairs (W helix) with two G·U and four rare C·A+ wobble pairs. Four adenines within the W helix are N1-protonated and wobble-base-paired with the opposing cytosine through two regular hydrogen bonds. Combined with the two G·U pairs, the C·A+ base pairs facilitate formation of a half turn of W-helical RNA flanked by six regular Watson-Crick base pairs in standard A conformation on either side. RNA melting experiments monitored by differential scanning calorimetry, UV and circular dichroism spectroscopy demonstrate that the RNA octadecamer undergoes a pH-induced structural transition which is consistent with the presence of a duplex with C·A+ base pairs at acidic pH. Our crystal structure provides a first glimpse of an RNA double helix based entirely on wobble base pairs with possible applications in RNA or DNA nanotechnology and pH biosensors.

Cooperativity between different tRNA modifications and their modification pathways

Ribonucleotide modifications perform a wide variety of roles in synthesis, turnover and functionality of tRNA molecules. The presence of particular chemical moieties can refine the internal interaction network within a tRNA molecule, influence its thermodynamic stability, contribute novel chemical properties and affect its decoding behavior during mRNA translation. As the lack of specific modifications in the anticodon stem and loop causes disrupted proteome homeostasis, diminished response to stress conditions, and the onset of human diseases, the underlying modification cascades have recently gained particular scientific and clinical interest. Nowadays, a complicated but conclusive image of the interconnectivity between different enzymatic modification cascades and their resulting tRNA modifications emerges. Here we summarize the current knowledge in the field, focusing on the known instances of cross talk among the enzymatic tRNA modification pathways and the consequences on the dynamic regulation of the tRNA modificome by various factors. This article is part of a Special Issue entitled: SI: Regulation of tRNA synthesis and modification in physiological conditions and disease edited by Dr. Boguta Magdalena.

The structure of an E. coli tRNAfMet A1-U72 variant shows an unusual conformation of the A1-U72 base pair

Translation initiation in eukaryotes and archaea involves a methionylated initiator tRNA delivered to the ribosome in a ternary complex with e/aIF2 and GTP. Eukaryotic and archaeal initiator tRNAs contain a highly conserved A₁-U₇₂ base pair at the top of the acceptor stem. The importance of this base pair to discriminate initiator tRNAs from elongator tRNAs has been established previously using genetics and biochemistry. However, no structural data illustrating how the A₁-U₇₂ base pair participates in the accurate selection of the initiator tRNAs by the translation initiation systems are available. Here, we describe the crystal structure of a mutant E. coli initiator tRNA_f^MetA₁-U₇₂, aminoacylated with methionine, in which the C₁:A₇₂ mismatch at the end of the tRNA acceptor stem has been changed to an A₁-U₇₂ base pair. Sequence alignments show that the mutant E. coli tRNA is a good mimic of archaeal initiator tRNAs. The crystal structure, determined at 2.8 Å resolution, shows that the A₁-U₇₂ pair adopts an unusual arrangement. A₁ is in a syn conformation and forms a single H-bond interaction with U₇₂ This interaction requires protonation of the N1 atom of A₁ Moreover, the 5' phosphoryl group folds back into the major groove of the acceptor stem and interacts with the N7 atom of G₂ A possible role of this unusual geometry of the A₁-U₇₂ pair in the recognition of the initiator tRNA by its partners during eukaryotic and archaeal translation initiation is discussed.

tRNA Modifications: Impact on Structure and Thermal Adaptation

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

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.

Pseudouridine Modification Inhibits Muscleblind-like 1 (MBNL1) Binding to CCUG Repeats and Minimally Structured RNA through Reduced RNA Flexibility

Myotonic dystrophy type 2 is a genetic neuromuscular disease caused by the expression of expanded CCUG repeat RNAs from the non-coding region of the CCHC-type zinc finger nucleic acid-binding protein (CNBP) gene. These CCUG repeats bind and sequester a family of RNA-binding proteins known as Muscleblind-like 1, 2, and 3 (MBNL1, MBNL2, and MBNL3), and sequestration plays a significant role in pathogenicity. MBNL proteins are alternative splicing regulators that bind to the consensus RNA sequence YGCY (Y = pyrimidine). This consensus sequence is found in the toxic RNAs (CCUG repeats) and in cellular RNA substrates that MBNL proteins have been shown to bind. Replacing the uridine in CCUG repeats with pseudouridine (Ψ) resulted in a modest reduction of MBNL1 binding. Interestingly, Ψ modification of a minimally structured RNA containing YGCY motifs resulted in more robust inhibition of MBNL1 binding. The different levels of inhibition between CCUG repeat and minimally structured RNA binding appear to be due to the ability to modify both pyrimidines in the YGCY motif, which is not possible in the CCUG repeats. Molecular dynamic studies of unmodified and pseudouridylated minimally structured RNAs suggest that reducing the flexibility of the minimally structured RNA leads to reduced binding by MBNL1.

TRUB1 is the predominant pseudouridine synthase acting on mammalian mRNA via a predictable and conserved code

Following synthesis, RNA can be modified with over 100 chemically distinct modifications, which can potentially regulate RNA expression post-transcriptionally. Pseudouridine (Ψ) was recently established to be widespread and dynamically regulated on yeast mRNA, but less is known about Ψ presence, regulation, and biogenesis in mammalian mRNA. Here, we sought to characterize the Ψ landscape on mammalian mRNA, to identify the main Ψ-synthases (PUSs) catalyzing Ψ formation, and to understand the factors governing their specificity toward selected targets. We first developed a framework allowing analysis, evaluation, and integration of Ψ mappings, which we applied to >2.5 billion reads from 30 human samples. These maps, complemented with genetic perturbations, allowed us to uncover TRUB1 and PUS7 as the two key PUSs acting on mammalian mRNA and to computationally model the sequence and structural elements governing the specificity of TRUB1, achieving near-perfect prediction of its substrates (AUC = 0.974). We then validated and extended these maps and the inferred specificity of TRUB1 using massively parallel reporter assays in which we monitored Ψ levels at thousands of synthetically designed sequence variants comprising either the sequences surrounding pseudouridylation targets or systematically designed mutants perturbing RNA sequence and structure. Our findings provide an extensive and high-quality characterization of the transcriptome-wide distribution of pseudouridine in human and the factors governing it and provide an important resource for the community, paving the path toward functional and mechanistic dissection of this emerging layer of post-transcriptional regulation.

Structural effects of modified ribonucleotides and magnesium in transfer RNAs

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

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.

Additive CHARMM force field for naturally occurring modified ribonucleotides

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Anions in Nucleic Acid Crystallography

Nucleic acid crystallization buffers contain a large variety of chemicals fitting specific needs. Among them, anions are often solely considered for pH-regulating purposes and as cationic co-salts while their ability to directly bind to nucleic acid structures is rarely taken into account. Here we review current knowledge related to the use of anions in crystallization buffers along with data on their biological prevalence. Chloride ions are frequently identified in crystal structures but display low cytosolic concentrations. Hence, they are thought to be distant from nucleic acid structures in the cell. Sulfate ions are also frequently identified in crystal structures but their localization in the cell remains elusive. Nevertheless, the characterization of the binding properties of these ions is essential for better interpreting the solvent structure in crystals and consequently, avoiding mislabeling of electron densities. Furthermore, understanding the binding properties of these anions should help to get clues related to their potential effects in crowded cellular environments.

Probing the influence of hypermodified residues within the tRNA3(Lys) anticodon stem loop interacting with the A-loop primer sequence from HIV-1

Replication of the HIV-1 virus requires reverse transcription of the viral RNA genome, a process that is specifically initiated by human tRNA3(Lys) packaged within the infectious virion. The primary binding site for the tRNA involves the 3' 18 nucleotides with an additional interaction between an adenine rich loop (A-loop) in the template and the anticodon stem-loop region of the tRNA3(Lys). The loop of the tRNA primer contains two hypermodified base residues and a pseudouridine that are required for a proper binding and activity. Here, we investigate the influence on the structure, dynamics and binding stability of the three modified residues (mnm(5)s(2)U34, t(6)A37 and Ψ39) using extensive molecular dynamics and Quantum Theory of Atoms in Molecules (QTAIM) analysis. Consistent with experiment, the results suggest that the three modified residues are required for faithful binding. Residues mnm(5)s(2)U34 and Ψ39 have a major influence in stabilizing the anticodon loop whereas mnm(5)s(2)U34 and t(6)A37 appear to stabilize the formation of the complex of tRNA3(Lys) with the HIV-1 A-loop.

Modified Nucleosides of Escherichia coli Ribosomal RNA

The modified nucleosides of RNA are chemically altered versions of the standard A, G, U, and C nucleosides. This review reviews the nature and location of the modified nucleosides of Escherichia coli rRNA, the enzymes that form them, and their known and/or putative functional role. There are seven Ψ (pseudouridines) synthases to make the 11 pseudouridines in rRNA. There is disparity in numbers because RluC and RluD each make 3 pseudouridines. Crystal structures have shown that the Ψ synthase domain is a conserved fold found only in all five families of Ψ synthases. The conversion of uridine to Ψ has no precedent in known metabolic reactions. Other enzymes are known to cleave the glycosyl bond but none carry out rotation of the base and rejoining to the ribose while still enzyme bound. Ten methyltransferases (MTs) are needed to make all the methylated nucleosides in 16S RNA, and 14 are needed for 23S RNA. Biochemical studies indicate that the modes of substrate recognition are idiosyncratic for each Ψ synthase since no common mode of recognition has been detected in studies of the seven synthases. Eight of the 24 expected MTs have been identified, and six crystal structures have been determined. Seven of the MTs and five of the structures are class I MTs with the appropriate protein fold plus unique appendages for the Ψ synthases. The remaining MT, RlmB, has the class IV trefoil knot fold.

Transfer RNA Modification

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

Accurate placement of substrate RNA by Gar1 in H/ACA RNA-guided pseudouridylation

H/ACA RNA-guided ribonucleoprotein particle (RNP), the most complicated RNA pseudouridylase so far known, uses H/ACA guide RNA for substrate capture and four proteins (Cbf5, Nop10, L7Ae and Gar1) for pseudouridylation. Although it was shown that Gar1 not only facilitates the product release, but also enhances the catalytic activity, the chemical role that Gar1 plays in this complicated machinery is largely unknown. Kinetics measurement on Pyrococcus furiosus RNPs at different temperatures making use of fluorescence anisotropy showed that Gar1 reduces the catalytic barrier through affecting the activation entropy instead of enthalpy. Site-directed mutagenesis combined with molecular dynamics simulations demonstrated that V149 in the thumb loop of Cbf5 is critical in placing the target uridine to the right position toward catalytic D85 of Cbf5. The enzyme elegantly aligns the position of uridine in the catalytic site with the help of Gar1. In addition, conversion of uridine to pseudouridine results in a rigid syn configuration of the target nucleotide in the active site and causes Gar1 to pull out the thumb. Both factors guarantee the efficient release of the product.

Evidence That Antibiotics Bind to Human Mitochondrial Ribosomal RNA Has Implications for Aminoglycoside Toxicity

Aminoglycosides are a well known antibiotic family used to treat bacterial infections in humans and animals, but which can be toxic. By binding to the decoding site of helix44 of the small subunit RNA of the bacterial ribosome, the aminoglycoside antibiotics inhibit protein synthesis, cause misreading, or obstruct peptidyl-tRNA translocation. Although aminoglycosides bind helix69 of the bacterial large subunit RNA as well, little is known about their interaction with the homologous human helix69. To probe the role this binding event plays in toxicity, changes to thermal stability, base stacking, and conformation upon aminoglycoside binding to the human cytoplasmic helix69 were compared with those of the human mitochondrial and Escherichia coli helix69. Surprisingly, binding of gentamicin and kanamycin A to the chemically synthesized terminal hairpins of the human cytoplasmic, human mitochondrial, and E. coli helix69 revealed similar dissociation constants (1.3-1.7 and 4.0-5.4 μM, respectively). In addition, aminoglycoside binding enhanced conformational stability of the human mitochondrial helix69 by increasing base stacking. Proton one-dimensional and two-dimensional NMR suggested significant and specific conformational changes of human mitochondrial and E. coli helix69 upon aminoglycoside binding, as compared with human cytoplasmic helix69. The conformational changes and similar aminoglycoside binding affinities observed for human mitochondrial helix69 and E. coli helix69, as well as the increase in structural stability shown for the former, suggest that this binding event is important to understanding aminoglycoside toxicity.

Controlling translation via modulation of tRNA levels

Transfer RNAs (tRNAs) are critical adaptor molecules that carry amino acids to a messenger RNA (mRNA) template during protein synthesis. Although tRNAs have commonly been viewed as abundant 'house-keeping' RNAs, it is becoming increasingly clear that tRNA expression is tightly regulated. Depending on a cell's proliferative status, the pool of active tRNAs is rapidly changed, enabling distinct translational programs to be expressed in differentiated versus proliferating cells. Here, I highlight several post-transcriptional regulatory mechanisms that allow the expression or functions of tRNAs to be altered. Modulating the modification status or structural stability of individual tRNAs can cause those specific tRNA transcripts to selectively accumulate or be degraded. Decay generally occurs via the rapid tRNA decay pathway or by the nuclear RNA surveillance machinery. In addition, the CCA-adding enzyme plays a critical role in determining the fate of a tRNA. The post-transcriptional addition of CCA to the 3' ends of stable tRNAs generates the amino acid attachment site, whereas addition of CCACCA to unstable tRNAs prevents aminoacylation and marks the tRNA for degradation. In response to various stresses, tRNAs can accumulate in the nucleus or be further cleaved into small RNAs, some of which inhibit translation. By implementing these various post-transcriptional control mechanisms, cells are able to fine-tune tRNA levels to regulate subsets of mRNAs as well as overall translation rates.