Naturally occurring modified ribonucleosides

Implications of RNA pseudouridylation for cancer biology and therapeutics: a narrative review

BACKGROUND

Pseudouridine (Ψ), a C5-glycoside isomer of uridine, stands as one of the most prevalent RNA modifications in all RNA types. Distinguishing from the C-N bond linking uridine to ribose, the link between Ψ and ribose is a C-C bond, endowing Ψ modified RNA distinct properties and functions in various biological processes. The conversion of uridine to Ψ is governed by pseudouridine synthases (PUSs). RNA pseudouridylation is implicated in cancer biology and therapeutics.

OBJECTIVES

In this review, we will summarize the methods for detecting Ψ, the process of Ψ generation, the impact of Ψ modification on RNA metabolism and gene expression, the roles of dysregulated Ψ and pseudouridine synthases in cancers, and the underlying mechanism.

METHODS

We conducted a comprehensive search of PubMed from its inception through February 2024. The search terms included "pseudouridine"; "pseudouridine synthase"; "PUS"; "dyskerin"; "cancer"; "tumor"; "carcinoma"; "malignancy"; "tumorigenesis"; "biomarker"; "prognosis" and "therapy". We included studies published in peer-reviewed journals that focused on Ψ detection, specific mechanisms involving Ψ and PUSs, and prognosis in cancer patients with high Ψ expression. We excluded studies lacking sufficient methodological details or appropriate controls.

RESULTS

Ψ has been recognized as a significant biomarker in cancer diagnosis and prognosis. Abnormal Ψ modifications mediated by various PUSs result in dysregulated RNA metabolism and impaired RNA function, promoting the development of various cancers. Overexpression of PUSs is common in cancer cells and predicts poor prognosis. PUSs inhibition arrests cell proliferation and enhances apoptosis in cancer cells, suggesting PUS-targeting cancer therapy may be a potential strategy in cancer treatment.

DISCUSSION

High Ψ levels in serum, urine, and saliva may suggest cancer, but do not specify the type, requiring additional lab markers and imaging for accurate diagnosis. Standardized detection methods are also crucial for reliable results. PUSs are linked to cancer, but more researches are needed to understand their mechanisms in different cancers. Anticancer treatments targeting PUSs are still under developed.

50 years in the making: acp3U, an amino-acid-containing nucleoside, links N-glycans and RNA in glycoRNA

In a recent publication in Cell, Xie et al.1 report a sensitive and scalable method for the detection and characterization of native glycoRNAs and identify acp3U, an abundant modified nucleoside discovered 50 years ago in tRNAPhe, as one of the primary attachment sites for N-glycans.

Combining Nanopore direct RNA sequencing with genetics and mass spectrometry for analysis of T-loop base modifications across 42 yeast tRNA isoacceptors

Transfer RNAs (tRNAs) contain dozens of chemical modifications. These modifications are critical for maintaining tRNA tertiary structure and optimizing protein synthesis. Here we advance the use of Nanopore direct RNA-sequencing (DRS) to investigate the synergy between modifications that are known to stabilize tRNA structure. We sequenced the 42 cytosolic tRNA isoacceptors from wild-type yeast and five tRNA-modifying enzyme knockout mutants. These data permitted comprehensive analysis of three neighboring and conserved modifications in T-loops: 5-methyluridine (m5U54), pseudouridine (Ψ55), and 1-methyladenosine (m1A58). Our results were validated using direct measurements of chemical modifications by mass spectrometry. We observed concerted T-loop modification circuits-the potent influence of Ψ55 for subsequent m1A58 modification on more tRNA isoacceptors than previously observed. Growing cells under nutrient depleted conditions also revealed a novel condition-specific increase in m1A58 modification on some tRNAs. A global and isoacceptor-specific classification strategy was developed to predict the status of T-loop modifications from a user-input tRNA DRS dataset, applicable to other conditions and tRNAs in other organisms. These advancements demonstrate how orthogonal technologies combined with genetics enable precise detection of modification landscapes of individual, full-length tRNAs, at transcriptome-scale.

Structural perspectives on adenosine to inosine RNA editing by ADARs

Adenosine deaminases acting on RNA (ADARs) are enzymes that catalyze the hydrolytic deamination of adenosine to inosine. The editing feature of ADARs has garnered much attention as a therapeutic tool to repurpose ADARs to correct disease-causing mutations at the mRNA level in a technique called site-directed RNA editing (SDRE). Administering a short guide RNA oligonucleotide that hybridizes to a mutant sequence forms the requisite dsRNA substrate, directing ADARs to edit the desired adenosine. However, much is still unknown about ADARs' selectivity and sequence-specific effects on editing. Atomic-resolution structures can help provide additional insight to ADARs' selectivity and lead to novel guide RNA designs. Indeed, recent structures of ADAR domains have expanded our understanding on RNA binding and the base-flipping catalytic mechanism. These efforts have enabled the rational design of improved ADAR guide strands and advanced the therapeutic potential of the SDRE approach. While no full-length structure of any ADAR is known, this review presents an exposition of the structural basis for function of the different ADAR domains, focusing on human ADAR2. Key insights are extrapolated to human ADAR1, which is of substantial interest because of its widespread expression in most human tissues.

Mycobacteria that cause tuberculosis have retained ancestrally acquired genes for the biosynthesis of chemically diverse terpene nucleosides

Mycobacterium tuberculosis (Mtb) releases the unusual terpene nucleoside 1-tuberculosinyladenosine (1-TbAd) to block lysosomal function and promote survival in human macrophages. Using conventional approaches, we found that genes Rv3377c and Rv3378c, but not Rv3376, were necessary for 1-TbAd biosynthesis. Here, we introduce linear models for mass spectrometry (limms) software as a next-generation lipidomics tool to study the essential functions of lipid biosynthetic enzymes on a whole-cell basis. Using limms, whole-cell lipid profiles deepened the phenotypic landscape of comparative mass spectrometry experiments and identified a large family of approximately 100 terpene nucleoside metabolites downstream of Rv3378c. We validated the identity of previously unknown adenine-, adenosine-, and lipid-modified tuberculosinol-containing molecules using synthetic chemistry and collisional mass spectrometry, including comprehensive profiling of bacterial lipids that fragment to adenine. We tracked terpene nucleoside genotypes and lipid phenotypes among Mycobacterium tuberculosis complex (MTC) species that did or did not evolve to productively infect either human or nonhuman mammals. Although 1-TbAd biosynthesis genes were thought to be restricted to the MTC, we identified the locus in unexpected species outside the MTC. Sequence analysis of the locus showed nucleotide usage characteristic of plasmids from plant-associated bacteria, clarifying the origin and timing of horizontal gene transfer to a pre-MTC progenitor. The data demonstrated correlation between high level terpene nucleoside biosynthesis and mycobacterial competence for human infection, and 2 mechanisms of 1-TbAd biosynthesis loss. Overall, the selective gain and evolutionary retention of tuberculosinyl metabolites in modern species that cause human TB suggest a role in human TB disease, and the newly discovered molecules represent candidate disease-specific biomarkers.

Chemo-enzymatic production of base-modified ATP analogues for polyadenylation of RNA

Base-modified adenosine-5'-triphosphate (ATP) analogues are highly sought after as building blocks for mRNAs and non-coding RNAs, for genetic code expansion or as inhibitors. Current synthetic strategies lack efficient and robust 5'-triphosphorylation of adenosine derivatives or rely on costly phosphorylation reagents. Here, we combine the efficient organic synthesis of base-modified AMP analogues with enzymatic phosphorylation by a promiscuous polyphosphate kinase 2 class III from an unclassified Erysipelotrichaceae bacterium (EbPPK2) to generate a panel of C2-, N6-, or C8-modified ATP analogues. These can be incorporated into RNA using template independent poly(A) polymerase. C2-halogenated ATP analogues were incorporated best, with incorporations of 300 to >1000 nucleotides forming hypermodified poly(A) tails.

Transfer RNA acetylation regulates in vivo mammalian stress signaling

Transfer RNA (tRNA) modifications are crucial for protein synthesis, but their position-specific physiological roles remain poorly understood. Here we investigate the impact of N4-acetylcytidine (ac 4 C), a highly conserved tRNA modification, using a Thumpd1 knockout mouse model. We find that loss of Thumpd1-dependent tRNA acetylation leads to reduced levels of tRNA Leu , increased ribosome stalling, and activation of eIF2α phosphorylation. Thumpd1 knockout mice exhibit growth defects and sterility. Remarkably, concurrent knockout of Thumpd1 and the stress-sensing kinase Gcn2 causes penetrant postnatal lethality, indicating a critical genetic interaction. Our findings demonstrate that a modification restricted to a single position within type II cytosolic tRNAs can regulate ribosome-mediated stress signaling in mammalian organisms, with implications for our understanding of translation control as well as therapeutic interventions.

ADATs: roles in tRNA editing and relevance to disease

Transfer RNAs (tRNAs) play central roles in protein biosynthesis. Post-transcriptional RNA modifications affect tRNA function and stability. Among these modifications, RNA editing is a widespread RNA modification in three domains of life. Proteins of the adenosine deaminase acting on tRNA (ADAT) family were discovered more than 20 years ago. They catalyze the deamination of adenosine to inosine (A-to-I) or cytidine to uridine (C-to-U) during tRNA maturation. The most studied example is the TadA- or ADAT2/3-mediated A-to-I conversion of the tRNA wobble position in the anticodon of prokaryotic or eukaryotic tRNAs, respectively. This review provides detailed information on A-to-I and C-to-U editing of tRNAs in different domains of life, presents recent new findings on ADATs for DNA editing, and finally comments on the association of mutations in the ADAT3 gene with intellectual disability.

A primer for junior trainees: Recognition of RNA modifications by RNA-binding proteins

The complexity of RNA cannot be fully expressed with the canonical A, C, G, and U alphabet. To date, over 170 distinct chemical modifications to RNA have been discovered in living systems. RNA modifications can profoundly impact the cellular outcomes of messenger RNAs (mRNAs), transfer and ribosomal RNAs, and noncoding RNAs. Additionally, aberrant RNA modifications are associated with human disease. RNA modifications are a rising topic within the fields of biochemistry and molecular biology. The role of RNA modifications in gene regulation, disease pathogenesis, and therapeutic applications increasingly captures the attention of the scientific community. This review aims to provide undergraduates, junior trainees, and educators with an appreciation for the significance of RNA modifications in eukaryotic organisms, alongside the skills required to identify and analyze fundamental RNA-protein interactions. The pumilio RNA-binding protein and YT521-B homology (YTH) family of modified RNA-binding proteins serve as examples to highlight the fundamental biochemical interactions that underlie the specific recognition of both unmodified and modified ribonucleotides, respectively. By instilling these foundational, textbook concepts through practical examples, this review contributes an analytical toolkit that facilitates engagement with RNA modifications research at large.

[Advances in mapping analysis of ribonucleic acid modifications through sequencing]

Over 170 chemical modifications have been discovered in various types of ribonucleic acids (RNAs), including messenger RNA (mRNA), ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA). These RNA modifications play crucial roles in a wide range of biological processes such as gene expression regulation, RNA stability maintenance, and protein translation. RNA modifications represent a new dimension of gene expression regulation known as the "epitranscriptome". The discovery of RNA modifications and the relevant writers, erasers, and readers provides an important basis for studies on the dynamic regulation and physiological functions of RNA modifications. Owing to the development of detection technologies for RNA modifications, studies on RNA epitranscriptomes have progressed to the single-base resolution, multilayer, and full-coverage stage. Transcriptome-wide methods help discover new RNA modification sites and are of great importance for elucidating the molecular regulatory mechanisms of epitranscriptomics, exploring the disease associations of RNA modifications, and understanding their clinical applications. The existing RNA modification sequencing technologies can be categorized according to the pretreatment approach and sequencing principle as direct high-throughput sequencing, antibody-enrichment sequencing, enzyme-assisted sequencing, chemical labeling-assisted sequencing, metabolic labeling sequencing, and nanopore sequencing technologies. These methods, as well as studies on the functions of RNA modifications, have greatly expanded our understanding of epitranscriptomics. In this review, we summarize the recent progress in RNA modification detection technologies, focusing on the basic principles, advantages, and limitations of different methods. Direct high-throughput sequencing methods do not require complex RNA pretreatment and allow for the mapping of RNA modifications using conventional RNA sequencing methods. However, only a few RNA modifications can be analyzed by high-throughput sequencing. Antibody enrichment followed by high-throughput sequencing has emerged as a crucial approach for mapping RNA modifications, significantly advancing the understanding of RNA modifications and their regulatory functions in different species. However, the resolution of antibody-enrichment sequencing is limited to approximately 100-200 bp. Although chemical crosslinking techniques can achieve single-base resolution, these methods are often complex, and the specificity of the antibodies used in these methods has raised concerns. In particular, the issue of off-target binding by the antibodies requires urgent attention. Enzyme-assisted sequencing has improved the accuracy of the localization analysis of RNA modifications and enables stoichiometric detection with single-base resolution. However, the enzymes used in this technique show poor reactivity, specificity, and sequence preference. Chemical labeling sequencing has become a widely used approach for profiling RNA modifications, particularly by altering reverse transcription (RT) signatures such as RT stops, misincorporations, and deletions. Chemical-assisted sequencing provides a sequence-independent RNA modification detection strategy that enables the localization of multiple RNA modifications. Additionally, when combined with the biotin-streptavidin affinity method, low-abundance RNA modifications can be enriched and detected. Nevertheless, the specificity of many chemical reactions remains problematic, and the development of specific reaction probes for particular modifications should continue in the future to achieve the precise localization of RNA modifications. As an indirect localization method, metabolic labeling sequencing specifically localizes the sites at which modifying enzymes act, which is of great significance in the study of RNA modification functions. However, this method is limited by the intracellular labeling of RNA and cannot be applied to biological samples such as clinical tissues and blood samples. Nanopore sequencing is a direct RNA-sequencing method that does not require RT or the polymerase chain reaction (PCR). However, challenges in analyzing the data obtained from nanopore sequencing, such as the high rate of false positives, must be resolved. Discussing sequencing analysis methods for various types of RNA modifications is instructive for the future development of novel RNA modification mapping technologies, and will aid studies on the functions of RNA modifications across the entire transcriptome.

The positional and numerical effect of N 6-methyladenosine in tracrRNA on the DNA cleavage activity of Cas9

Post-transcriptional modifications on the guide RNAs utilized in the Cas9 system may have the potential to impact the activity of Cas9. In this study, we synthesized a series of tracrRNAs containing N 6-methyadenosine (m6A), a prevalent post-transcriptional modification, at various positions. We evaluated the effect of these modifications on the DNA cleavage activity of Cas9. Our results show that multiple m6As in the anti-repeat region of tracrRNA reduce the DNA cleavage activity of Cas9. This suggests that the m6A-modified tracrRNA can be used for Cas9 only when the number and the position of the modified residue are properly chosen in tracrRNA.

A Visual Compendium of Principal Modifications within the Nucleic Acid Sugar Phosphate Backbone

Nucleic acid chemistry is a huge research area that has received new impetus due to the recent explosive success of oligonucleotide therapy. In order for an oligonucleotide to become clinically effective, its monomeric parts are subjected to modifications. Although a large number of redesigned natural nucleic acids have been proposed in recent years, the vast majority of them are combinations of simple modifications proposed over the past 50 years. This review is devoted to the main modifications of the sugar phosphate backbone of natural nucleic acids known to date. Here, we propose a systematization of existing knowledge about modifications of nucleic acid monomers and an acceptable classification from the point of view of chemical logic. The visual representation is intended to inspire researchers to create a new type of modification or an original combination of known modifications that will produce unique oligonucleotides with valuable characteristics.

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.

Comparative genomics of Deinococcus radiodurans: unveiling genetic discrepancies between ATCC 13939K and BAA-816 strains

The Deinococcus genus is renowned for its remarkable resilience against environmental stresses, including ionizing radiation, desiccation, and oxidative damage. This resilience is attributed to its sophisticated DNA repair mechanisms and robust defense systems, enabling it to recover from extensive damage and thrive under extreme conditions. Central to Deinococcus research, the D. radiodurans strains ATCC BAA-816 and ATCC 13939 facilitate extensive studies into this remarkably resilient genus. This study focused on delineating genetic discrepancies between these strains by sequencing our laboratory's ATCC 13939 specimen (ATCC 13939K) and juxtaposing it with ATCC BAA-816. We uncovered 436 DNA sequence differences within ATCC 13939K, including 100 single nucleotide variations, 278 insertions, and 58 deletions, which could induce frameshifts altering protein-coding genes. Gene annotation revisions accounting for gene fusions and the reconciliation of gene lengths uncovered novel protein-coding genes and refined the functional categorizations of established ones. Additionally, the analysis pointed out genome structural variations due to insertion sequence (IS) elements, underscoring the D. radiodurans genome's plasticity. Notably, ATCC 13939K exhibited a loss of six ISDra2 elements relative to BAA-816, restoring genes fragmented by ISDra2, such as those encoding for α/β hydrolase and serine protease, and revealing new open reading frames, including genes imperative for acetoin decomposition. This comparative genomic study offers vital insights into the metabolic capabilities and resilience strategies of D. radiodurans.

Identification of Amino Acids in Trm734 Required for 2'-O-Methylation of the tRNAPhe Wobble Residue

All organisms methylate their nucleic acids, and this methylation is critical for proper gene expression at both the transcriptional and translational levels. For proper translation in eukaryotes, 2'-O-methylation of C32 (Cm32) and G34 (Gm34) in the anticodon loop of tRNAPhe is critical, with defects in these modifications associated with human disease. In yeast, Cm32 is formed by an enzyme that consists of the methyltransferase Trm7 in complex with the auxiliary protein Trm732, and Gm34 is formed by an enzyme that consists of Trm7 in complex with Trm734. The role of Trm732 and Trm734 in tRNA modification is not fully understood, although previous studies have suggested that Trm734 is important for tRNA binding. In this report, we generated Trm734 variants and tested their ability to work with Trm7 to modify tRNAPhe. Using this approach, we identified several regions of amino acids that are important for Trm734 activity and/or stability. Based on the previously determined Trm7-Trm734 crystal structure, these crucial amino acids are near the active site of Trm7 and are not directly involved in Trm7-Trm734 protein-protein interactions. Immunoprecipitation experiments with these Trm734 variants and Trm7 confirm that these residues are not involved in Trm7-Trm734 binding. Further experiments should help determine if these residues are important for tRNA binding or have another role in the modification of the tRNA. Furthermore, our discovery of a nonfunctional, stable Trm734 variant will be useful in determining if the reported roles of Trm734 in other biological processes such as retromer processing and resistance to Ty1 transposition are due to tRNA modification defects or to other bona fide cellular roles of Trm734.

Functional redundancy in tRNA dihydrouridylation

Dihydrouridine (D) is a common modified base found predominantly in transfer RNA (tRNA). Despite its prevalence, the mechanisms underlying dihydrouridine biosynthesis, particularly in prokaryotes, have remained elusive. Here, we conducted a comprehensive investigation into D biosynthesis in Bacillus subtilis through a combination of genetic, biochemical, and epitranscriptomic approaches. Our findings reveal that B. subtilis relies on two FMN-dependent Dus-like flavoprotein homologs, namely DusB1 and DusB2, to introduce all D residues into its tRNAs. Notably, DusB1 exhibits multisite enzyme activity, enabling D formation at positions 17, 20, 20a and 47, while DusB2 specifically catalyzes D biosynthesis at positions 20 and 20a, showcasing a functional redundancy among modification enzymes. Extensive tRNA-wide D-mapping demonstrates that this functional redundancy impacts the majority of tRNAs, with DusB2 displaying a higher dihydrouridylation efficiency compared to DusB1. Interestingly, we found that BsDusB2 can function like a BsDusB1 when overexpressed in vivo and under increasing enzyme concentration in vitro. Furthermore, we establish the importance of the D modification for B. subtilis growth at suboptimal temperatures. Our study expands the understanding of D modifications in prokaryotes, highlighting the significance of functional redundancy in this process and its impact on bacterial growth and adaptation.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Targeted i6A-RNA degradation through sequential Fluorination-Azidation and Click reaction with imidazole-based probes

We report the use of trimethylsilyl azide and Selectfluor to implement a standard protocol targeted at the prenylated nucleic acid known as i6A-RNA. After optimizing the conditions, we applied this method to regulate a wide range of i6A-RNA species using synthetic imidazole-based probes (I-IV). We observed that prenylated nucleic acid plays a crucial role in the cell hemostasis in A549 cell lines.

A stochastic RNA editing process targets a limited number of sites in individual Drosophila glutamatergic motoneurons

RNA editing is a post-transcriptional source of protein diversity and occurs across the animal kingdom. Given the complete profile of mRNA targets and their editing rate in individual cells is unclear, we analyzed single cell RNA transcriptomes from Drosophila larval tonic and phasic glutamatergic motoneuron subtypes to determine the most highly edited targets and identify cell-type specific editing. From ∼15,000 genes encoded in the genome, 316 high confidence A-to-I canonical RNA edit sites were identified, with 102 causing missense amino acid changes in proteins regulating membrane excitability, synaptic transmission, and cellular function. Some sites showed 100% editing in single neurons as observed with mRNAs encoding mammalian AMPA receptors. However, most sites were edited at lower levels and generated variable expression of edited and unedited mRNAs within individual neurons. Together, these data provide insights into how the RNA editing landscape alters protein function to modulate the properties of two well-characterized neuronal populations in Drosophila .

Cryo-EM structures of the human Elongator complex at work

tRNA modifications affect ribosomal elongation speed and co-translational folding dynamics. The Elongator complex is responsible for introducing 5-carboxymethyl at wobble uridine bases (cm5U34) in eukaryotic tRNAs. However, the structure and function of human Elongator remain poorly understood. In this study, we present a series of cryo-EM structures of human ELP123 in complex with tRNA and cofactors at four different stages of the reaction. The structures at resolutions of up to 2.9 Å together with complementary functional analyses reveal the molecular mechanism of the modification reaction. Our results show that tRNA binding exposes a universally conserved uridine at position 33 (U33), which triggers acetyl-CoA hydrolysis. We identify a series of conserved residues that are crucial for the radical-based acetylation of U34 and profile the molecular effects of patient-derived mutations. Together, we provide the high-resolution view of human Elongator and reveal its detailed mechanism of action.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

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

Research progress of RNA pseudouridine modification in nervous system

Recent advances of pseudouridine (Ψ, 5-ribosyluracil) modification highlight its crucial role as a post-transcriptional regulator in gene expression and its impact on various RNA processes. Ψ synthase (PUS), a category of RNA-modifying enzymes, orchestrates the pseudouridylation reaction. It can specifically recognize conserved sequences or structural motifs within substrates, thereby regulating the biological function of various RNA molecules accurately. Our comprehensive review underscored the close association of PUS1, PUS3, PUS7, PUS10, and dyskerin PUS1 with various nervous system disorders, including neurodevelopmental disorders, nervous system tumors, mitochondrial myopathy, lactic acidosis and sideroblastic anaemia (MLASA) syndrome, peripheral nervous system disorders, and type II myotonic dystrophy. In light of these findings, this study elucidated how Ψ strengthened RNA structures and contributed to RNA function, thereby providing valuable insights into the intricate molecular mechanisms underlying nervous system diseases. However, the detailed effects and mechanisms of PUS on neuron remain elusive. This lack of mechanistic understanding poses a substantial obstacle to the development of therapeutic approaches for various neurological disorders based on Ψ modification.

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

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

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.

In silico λ-dynamics predicts protein binding specificities to modified RNAs

RNA modifications shape gene expression through a smorgasbord of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

MODOMICS: a database of RNA modifications and related information. 2023 update

The MODOMICS database was updated with recent data and now includes new data types related to RNA modifications. Changes to the database include an expanded modification catalog, encompassing both natural and synthetic residues identified in RNA structures. This addition aids in representing RNA sequences from the RCSB PDB database more effectively. To manage the increased number of modifications, adjustments to the nomenclature system were made. Updates in the RNA sequences section include the addition of new sequences and the reintroduction of sequence alignments for tRNAs and rRNAs. The protein section was updated and connected to structures from the RCSB PDB database and predictions by AlphaFold. MODOMICS now includes a data annotation system, with 'Evidence' and 'Estimated Reliability' features, offering clarity on data support and accuracy. This system is open to all MODOMICS entries, enhancing the accuracy of RNA modification data representation. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

Mitochondrial RNA maturation

The vast majority of oxygen-utilizing eukaryotes need to express their own mitochondrial genome, mtDNA, to survive. In comparison to size of their nuclear genome, mtDNA is minimal, even in the most exceptional examples. Having evolved from bacteria in an endosymbiotic event, it might be expected that the process of mtDNA expression would be relatively simple. The aim of this short review is to illustrate just how wrong this assumption is. The production of functional mitochondrial RNA across species evolved in many directions. Organelles use a dizzying array of RNA processing, modifying, editing, splicing and maturation events that largely require the import of nuclear-encoded proteins from the cytosol. These processes are sometimes driven by the unusual behaviour of the mitochondrial genome from which the RNA is originally transcribed, but in many examples the complex processes that are essential for the production of functional RNA in the organelle, are fascinating and bewildering.

Integrative Approach to Probe Alternative Redox Mechanisms in RNA Modifications

RNA modifications found in most RNAs, particularly in tRNAs and rRNAs, reveal an abundance of chemical alterations of nucleotides. Over 150 distinct RNA modifications are known, emphasizing a remarkable diversity of chemical moieties in RNA molecules. These modifications play pivotal roles in RNA maturation, structural integrity, and the fidelity and efficiency of translation processes. The catalysts responsible for these modifications are RNA-modifying enzymes that use a striking array of chemistries to directly influence the chemical landscape of RNA. This diversity is further underscored by instances where the same modification is introduced by distinct enzymes that use unique catalytic mechanisms and cofactors across different domains of life. This phenomenon of convergent evolution highlights the biological importance of RNA modification and the vast potential within the chemical repertoire for nucleotide alteration. While shared RNA modifications can hint at conserved enzymatic pathways, a major bottleneck is to identify alternative routes within species that possess a modified RNA but are devoid of known RNA-modifying enzymes. To address this challenge, a combination of bioinformatic and experimental strategies proves invaluable in pinpointing new genes responsible for RNA modifications. This integrative approach not only unveils new chemical insights but also serves as a wellspring of inspiration for biocatalytic applications and drug design. In this Account, we present how comparative genomics and genome mining, combined with biomimetic synthetic chemistry, biochemistry, and anaerobic crystallography, can be judiciously implemented to address unprecedented and alternative chemical mechanisms in the world of RNA modification. We illustrate these integrative methodologies through the study of tRNA and rRNA modifications, dihydrouridine, 5-methyluridine, queuosine, 8-methyladenosine, 5-carboxymethylamino-methyluridine, or 5-taurinomethyluridine, each dependent on a diverse array of redox chemistries, often involving organic compounds, organometallic complexes, and metal coenzymes. We explore how vast genome and tRNA databases empower comparative genomic analyses and enable the identification of novel genes that govern RNA modification. Subsequently, we describe how the isolation of a stable reaction intermediate can guide the synthesis of a biomimetic to unveil new enzymatic pathways. We then discuss the usefulness of a biochemical "shunt" strategy to study catalytic mechanisms and to directly visualize reactive intermediates bound within active sites. While we primarily focus on various RNA-modifying enzymes studied in our laboratory, with a particular emphasis on the discovery of a SAM-independent methylation mechanism, the strategies and rationale presented herein are broadly applicable for the identification of new enzymes and the elucidation of their intricate chemistries. This Account offers a comprehensive glimpse into the evolving landscape of RNA modification research and highlights the pivotal role of integrated approaches to identify novel enzymatic pathways.

The structural basis of mRNA recognition and binding by yeast pseudouridine synthase PUS1

The chemical modification of RNA bases represents a ubiquitous activity that spans all domains of life. Pseudouridylation is the most common RNA modification and is observed within tRNA, rRNA, ncRNA and mRNAs. Pseudouridine synthase or 'PUS' enzymes include those that rely on guide RNA molecules and others that function as 'stand-alone' enzymes. Among the latter, several have been shown to modify mRNA transcripts. Although recent studies have defined the structural requirements for RNA to act as a PUS target, the mechanisms by which PUS1 recognizes these target sequences in mRNA are not well understood. Here we describe the crystal structure of yeast PUS1 bound to an RNA target that we identified as being a hot spot for PUS1-interaction within a model mRNA at 2.4 Å resolution. The enzyme recognizes and binds both strands in a helical RNA duplex, and thus guides the RNA containing the target uridine to the active site for subsequent modification of the transcript. The study also allows us to show the divergence of related PUS1 enzymes and their corresponding RNA target specificities, and to speculate on the basis by which PUS1 binds and modifies mRNA or tRNA substrates.

TudS desulfidases recycle 4-thiouridine-5'-monophosphate at a catalytic [4Fe-4S] cluster

In all domains of life, transfer RNAs (tRNAs) contain post-transcriptionally sulfur-modified nucleosides such as 2- and 4-thiouridine. We have previously reported that a recombinant [4Fe-4S] cluster-containing bacterial desulfidase (TudS) from an uncultured bacterium catalyzes the desulfuration of 2- and 4-thiouracil via a [4Fe-5S] cluster intermediate. However, the in vivo function of TudS enzymes has remained unclear and direct evidence for substrate binding to the [4Fe-4S] cluster during catalysis was lacking. Here, we provide kinetic evidence that 4-thiouridine-5'-monophosphate rather than sulfurated tRNA, thiouracil, thiouridine or 4-thiouridine-5'-triphosphate is the preferred substrate of TudS. The occurrence of sulfur- and substrate-bound catalytic intermediates was uncovered from the observed switch of the S = 3/2 spin state of the catalytic [4Fe-4S] cluster to a S = 1/2 spin state upon substrate addition. We show that a putative gene product from Pseudomonas putida KT2440 acts as a TudS desulfidase in vivo and conclude that TudS-like enzymes are widespread desulfidases involved in recycling and detoxifying tRNA-derived 4-thiouridine monophosphate nucleosides for RNA synthesis.

Mitochondrial RNA m3C methyltransferase METTL8 relies on an isoform-specific N-terminal extension and modifies multiple heterogenous tRNAs

Methyltransferase-like 8 (METTL8) encodes a mitochondria-localized METTL8-Iso1 and a nucleolus-distributed METTL8-Iso4 isoform, which differ only in their N-terminal extension (N-extension), by mRNA alternative splicing. METTL8-Iso1 generates 3-methylcytidine at position 32 (m3C32) of mitochondrial tRNAThr and tRNASer(UCN). Whether METTL8-Iso4 is an active m3C32 methyltransferase and the role of the N-extension in mitochondrial tRNA m3C32 formation remain unclear. Here, we revealed that METTL8-Iso4 was inactive in m3C32 generation due to the lack of N-extension, which contains several absolutely conserved modification-critical residues; the counterparts were likewise essential in cytoplasmic m3C32 biogenesis by methyltransferase-like 2A (METTL2A) or budding yeasts tRNA N3-methylcytidine methyltransferase (Trm140), in vitro and in vivo. Cross-compartment/species tRNA modification assays unexpectedly found that METTL8-Iso1 efficiently introduced m3C32 to several cytoplasmic or even bacterial tRNAs in vitro. m3C32 did not influence tRNAThrN6-threonylcarbamoyladenosine (t6A) modification or aminoacylation. In addition to its interaction with mitochondrial seryl-tRNA synthetase (SARS2), we further discovered an interaction between mitochondrial threonyl-tRNA synthetase (TARS2) and METTL8-Iso1. METTL8-Iso1 substantially stimulated the aminoacylation activities of SARS2 and TARS2 in vitro, suggesting a functional connection between mitochondrial tRNA modification and charging. Altogether, our results deepen the mechanistic insights into mitochondrial m3C32 biogenesis and provide a valuable route to prepare cytoplasmic/bacterial tRNAs with only a m3C32 moiety, aiding in future efforts to investigate its effects on tRNA structure and function.

N4-acetylation of cytidine in mRNA plays essential roles in plants

The biological function of RNA can be modulated by base modifications. Here, we unveiled the occurrence of N4-acetylation of cytidine in plant RNA, including mRNA, by employing LC-MS/MS and acRIP-seq. We identified 325 acetylated transcripts from the leaves of 4-week-old Arabidopsis (Arabidopsis thaliana) plants and determined that 2 partially redundant N-ACETYLTRANSFERASEs FOR CYTIDINE IN RNA (ACYR1 and ACYR2), which are homologous to mammalian NAT10, are required for acetylating RNA in vivo. A double-null mutant was embryo lethal, while eliminating 3 of the 4 ACYR alleles led to defects in leaf development. These phenotypes could be traced back to the reduced acetylation and concomitant destabilization of the transcript of TOUGH, which is required for miRNA processing. These findings indicate that N4-acetylation of cytidine is a modulator of RNA function with a critical role in plant development and likely many other processes.

DNA Structure Design Is Improved Using an Artificially Expanded Alphabet of Base Pairs Including Loop and Mismatch Thermodynamic Parameters

We show that in silico design of DNA secondary structures is improved by extending the base pairing alphabet beyond A-T and G-C to include the pair between 2-amino-8-(1'-β-d-2'-deoxyribofuranosyl)-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one and 6-amino-3-(1'-β-d-2'-deoxyribofuranosyl)-5-nitro-(1H)-pyridin-2-one, abbreviated as P and Z. To obtain the thermodynamic parameters needed to include P-Z pairs in the designs, we performed 47 optical melting experiments and combined the results with previous work to fit free energy and enthalpy nearest neighbor folding parameters for P-Z pairs and G-Z wobble pairs. We find G-Z pairs have stability comparable to that of A-T pairs and should therefore be included as base pairs in structure prediction and design algorithms. Additionally, we extrapolated the set of loop, terminal mismatch, and dangling end parameters to include the P and Z nucleotides. These parameters were incorporated into the RNAstructure software package for secondary structure prediction and analysis. Using the RNAstructure Design program, we solved 99 of the 100 design problems posed by Eterna using the ACGT alphabet or supplementing it with P-Z pairs. Extending the alphabet reduced the propensity of sequences to fold into off-target structures, as evaluated by the normalized ensemble defect (NED). The NED values were improved relative to those from the Eterna example solutions in 91 of 99 cases in which Eterna-player solutions were provided. P-Z-containing designs had average NED values of 0.040, significantly below the 0.074 of standard-DNA-only designs, and inclusion of the P-Z pairs decreased the time needed to converge on a design. This work provides a sample pipeline for inclusion of any expanded alphabet nucleotides into prediction and design workflows.

Ghost authors revealed: The structure and function of human N6 -methyladenosine RNA methyltransferases

Despite the discovery of modified nucleic acids nearly 75 years ago, their biological functions are still being elucidated. N6 -methyladenosine (m6 A) is the most abundant modification in eukaryotic messenger RNA (mRNA) and has also been detected in non-coding RNAs, including long non-coding RNA, ribosomal RNA, and small nuclear RNA. In general, m6 A marks can alter RNA secondary structure and initiate unique RNA-protein interactions that can alter splicing, mRNA turnover, and translation, just to name a few. Although m6 A marks in human RNAs have been known to exist since 1974, the structures and functions of methyltransferases responsible for writing m6 A marks have been established only recently. Thus far, there are four confirmed human methyltransferases that catalyze the transfer of a methyl group from S-adenosylmethionine (SAM) to the N6 position of adenosine, producing m6 A: methyltransferase-like protein (METTL) 3/METTL14 complex, METTL16, METTL5, and zinc-finger CCHC-domain-containing protein 4. Though the methyltransferases have unique RNA targets, all human m6 A RNA methyltransferases contain a Rossmann fold with a conserved SAM-binding pocket, suggesting that they utilize a similar catalytic mechanism for methyl transfer. For each of the human m6 A RNA methyltransferases, we present the biological functions and links to human disease, RNA targets, catalytic and kinetic mechanisms, and macromolecular structures. We also discuss m6 A marks in human viruses and parasites, assigning m6 A marks in the transcriptome to specific methyltransferases, small molecules targeting m6 A methyltransferases, and the enzymes responsible for hypermodified m6 A marks and their biological functions in humans. Understanding m6 A methyltransferases is a critical steppingstone toward establishing the m6 A epitranscriptome and more broadly the RNome. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Widespread 8-oxoguanine modifications of miRNA seeds differentially regulate redox-dependent cancer development

Oxidative stress contributes to tumourigenesis by altering gene expression. One accompanying modification, 8-oxoguanine (o8G) can change RNA-RNA interactions via o8G•A base pairing, but its regulatory roles remain elusive. Here, on the basis of o8G-induced guanine-to-thymine (o8G > T) variations featured in sequencing, we discovered widespread position-specific o8Gs in tumour microRNAs, preferentially oxidized towards 5' end seed regions (positions 2-8) with clustered sequence patterns and clinically associated with patients in lower-grade gliomas and liver hepatocellular carcinoma. We validated that o8G at position 4 of miR-124 (4o8G-miR-124) and 4o8G-let-7 suppress lower-grade gliomas, whereas 3o8G-miR-122 and 4o8G-let-7 promote malignancy of liver hepatocellular carcinoma by redirecting the target transcriptome to oncogenic regulatory pathways. Stepwise oxidation from tumour-promoting 3o8G-miR-122 to tumour-suppressing 2,3o8G-miR-122 occurs and its specific modulation in mouse liver effectively attenuates diethylnitrosamine-induced hepatocarcinogenesis. These findings provide resources and insights into epitranscriptional o8G regulation of microRNA functions, reprogrammed by redox changes, implicating its control for cancer treatment.

Pyrimidine catabolism is required to prevent the accumulation of 5-methyluridine in RNA

5-Methylated cytosine is a frequent modification in eukaryotic RNA and DNA influencing mRNA stability and gene expression. Here we show that free 5-methylcytidine (5mC) and 5-methyl-2'-deoxycytidine are generated from nucleic acid turnover in Arabidopsis thaliana, and elucidate how these cytidines are degraded, which is unclear in eukaryotes. First CYTIDINE DEAMINASE produces 5-methyluridine (5mU) and thymidine which are subsequently hydrolyzed by NUCLEOSIDE HYDROLASE 1 (NSH1) to thymine and ribose or deoxyribose. Interestingly, far more thymine is generated from RNA than from DNA turnover, and most 5mU is directly released from RNA without a 5mC intermediate, since 5-methylated uridine (m5U) is an abundant RNA modification (m5U/U ∼1%) in Arabidopsis. We show that m5U is introduced mainly by tRNA-SPECIFIC METHYLTRANSFERASE 2A and 2B. Genetic disruption of 5mU degradation in the NSH1 mutant causes m5U to occur in mRNA and results in reduced seedling growth, which is aggravated by external 5mU supplementation, also leading to more m5U in all RNA species. Given the similarities between pyrimidine catabolism in plants, mammals and other eukaryotes, we hypothesize that the removal of 5mU is an important function of pyrimidine degradation in many organisms, which in plants serves to protect RNA from stochastic m5U modification.

Promoter-independent synthesis of chemically modified RNA by human DNA polymerase θ variants

Synthetic RNA oligonucleotides composed of canonical and modified ribonucleotides are highly effective for RNA antisense therapeutics and RNA-based genome engineering applications utilizing CRISPR-Cas9. Yet, synthesis of synthetic RNA using phosphoramidite chemistry is highly inefficient and expensive relative to DNA oligonucleotides, especially for relatively long RNA oligonucleotides. Thus, new biotechnologies are needed to significantly reduce costs, while increasing synthesis rates and yields of synthetic RNA. Here, we engineer human DNA polymerase theta (Polθ) variants and demonstrate their ability to synthesize long (95-200 nt) RNA oligonucleotides with canonical ribonucleotides and ribonucleotide analogs commonly used for stabilizing RNA for therapeutic and genome engineering applications. In contrast to natural promoter-dependent RNA polymerases, Polθ variants synthesize RNA by initiating from DNA or RNA primers, which enables the production of RNA without short abortive byproducts. Remarkably, Polθ variants show the lower capacity to misincorporate ribonucleotides compared to T7 RNA polymerase. Automation of this enzymatic RNA synthesis technology can potentially increase yields while reducing costs of synthetic RNA oligonucleotide production.

Direct sequencing of total Saccharomyces cerevisiae tRNAs by LC-MS/MS

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

The diverse structural modes of tRNA binding and recognition

tRNAs are short noncoding RNAs responsible for decoding mRNA codon triplets, delivering correct amino acids to the ribosome, and mediating polypeptide chain formation. Due to their key roles during translation, tRNAs have a highly conserved shape and large sets of tRNAs are present in all living organisms. Regardless of sequence variability, all tRNAs fold into a relatively rigid three-dimensional L-shaped structure. The conserved tertiary organization of canonical tRNA arises through the formation of two orthogonal helices, consisting of the acceptor and anticodon domains. Both elements fold independently to stabilize the overall structure of tRNAs through intramolecular interactions between the D- and T-arm. During tRNA maturation, different modifying enzymes posttranscriptionally attach chemical groups to specific nucleotides, which not only affect translation elongation rates but also restrict local folding processes and confer local flexibility when required. The characteristic structural features of tRNAs are also employed by various maturation factors and modification enzymes to assure the selection, recognition, and positioning of specific sites within the substrate tRNAs. The cellular functional repertoire of tRNAs continues to extend well beyond their role in translation, partly, due to the expanding pool of tRNA-derived fragments. Here, we aim to summarize the most recent developments in the field to understand how three-dimensional structure affects the canonical and noncanonical functions of tRNA.

Single-Base Resolution Detection of N6-Methyladenosine in RNA by Adenosine Deamination Sequencing

N⁶-Methyladenosine (m⁶A) is one of the most abundant and prevalent natural modifications occurring in diverse RNA species. m⁶A plays a wide range of roles in physiological and pathological processes. Revealing the functions of m⁶A relies on the faithful detection of individual m⁶A sites in RNA. However, developing a simple method for the single-base resolution detection of m⁶A is still a challenging task. Herein, we report an adenosine deamination sequencing (AD-seq) technique for the facile detection of m⁶A in RNA at single-base resolution. The AD-seq approach capitalizes on the selective deamination of adenosine, but not m⁶A, by the evolved tRNA adenosine deaminase (TadA) variant of TadA8e or the dimer protein of TadA-TadA8e. In AD-seq, adenosine is deaminated by TadA8e or TadA-TadA8e to form inosine, which pairs with cytidine and is read as guanosine in sequencing. m⁶A resists deamination due to the interference of the methyl group at the N6 position of adenosine. Thus, the m⁶A base pairs with thymine and is still read as adenosine in sequencing. The differential readouts from A and m⁶A in sequencing can achieve the single-base resolution detection of m⁶A in RNA. Application of the proposed AD-seq successfully identified individual m⁶A sites in Escherichia coli 23S rRNA. Taken together, the proposed AD-seq allows simple and cost-effective detection of m⁶A at single-base resolution in RNA, which provides a valuable tool to decipher the functions of m⁶A in RNA.

Methylated guanosine and uridine modifications in S. cerevisiae mRNAs modulate translation elongation

Chemical modifications to protein encoding messenger RNAs (mRNAs) influence their localization, translation, and stability within cells. Over 15 different types of mRNA modifications have been observed by sequencing and liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) approaches. While LC-MS/MS is arguably the most essential tool available for studying analogous protein post-translational modifications, the high-throughput discovery and quantitative characterization of mRNA modifications by LC-MS/MS has been hampered by the difficulty of obtaining sufficient quantities of pure mRNA and limited sensitivities for modified nucleosides. We have overcome these challenges by improving the mRNA purification and LC-MS/MS pipelines. The methodologies we developed result in no detectable non-coding RNA modifications signals in our purified mRNA samples, quantify 50 ribonucleosides in a single analysis, and provide the lowest limit of detection reported for ribonucleoside modification LC-MS/MS analyses. These advancements enabled the detection and quantification of 13 S. cerevisiae mRNA ribonucleoside modifications and reveal the presence of four new S. cerevisiae mRNA modifications at low to moderate levels (1-methyguanosine, N2-methylguanosine, N2,N2-dimethylguanosine, and 5-methyluridine). We identified four enzymes that incorporate these modifications into S. cerevisiae mRNAs (Trm10, Trm11, Trm1, and Trm2, respectively), though our results suggest that guanosine and uridine nucleobases are also non-enzymatically methylated at low levels. Regardless of whether they are incorporated in a programmed manner or as the result of RNA damage, we reasoned that the ribosome will encounter the modifications that we detect in cells. To evaluate this possibility, we used a reconstituted translation system to investigate the consequences of modifications on translation elongation. Our findings demonstrate that the introduction of 1-methyguanosine, N2-methylguanosine and 5-methyluridine into mRNA codons impedes amino acid addition in a position dependent manner. This work expands the repertoire of nucleoside modifications that the ribosome must decode in S. cerevisiae. Additionally, it highlights the challenge of predicting the effect of discrete modified mRNA sites on translation de novo because individual modifications influence translation differently depending on mRNA sequence context.

NMR of RNA - Structure and interactions

RNA was shown to have a more substantial role in the regulation of diverse cellular processes than anticipated until recently. Answers to questions what is the structure of specific RNAs, how structure changes to accommodate different functional roles, and how RNA senses other biomolecules and changes its fold upon interaction create a complete representation of RNA involved in cellular processes. Nuclear magnetic resonance (NMR) spectroscopy encompasses a collection of methods and approaches that offer insight into several structural aspects of RNAs. We review the most recent advances in the field of viral, long non-coding, regulatory, and four-stranded RNAs, with an emphasis on the detection of dynamic sub-states and in view of chemical modifications that expand RNA's function.

Identification of 2-Methylthio-methylenethio-N6 -(cis-4-hydroxyisopentenyl)adenosine (msms2 io6 A37 ) as a Novel Modification at Adenosine 37 of tRNAs from Salmonella typhimurium

Post-transcriptional modifications of tRNA nucleotide are important determinants in folding, structure and function. We have successfully identified and characterized a new modified base named 2-methylthio-methylenethio-N⁶ -(cis-4-hydroxyisopentenyl)adenosine, which is present at position 37 in some tRNAs. We also showed that this new modified adenosine is derived from the known 2-methylthio-methylenethio-N⁶ -(isopentenyl)adenosine nucleoside by a catalytic cycle of the tRNA-diiron monooxygenase, MiaE, present in Salmonella typhimurium.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Cytidine deaminases catalyze the conversion of N(S,O)4-substituted pyrimidine nucleosides

Cytidine deaminases (CDAs) catalyze the hydrolytic deamination of cytidine and 2'-deoxycytidine to uridine and 2'-deoxyuridine. Here, we report that prokaryotic homo-tetrameric CDAs catalyze the nucleophilic substitution at the fourth position of N⁴-acyl-cytidines, N⁴-alkyl-cytidines, and N⁴-alkyloxycarbonyl-cytidines, and S⁴-alkylthio-uridines and O⁴-alkyl-uridines, converting them to uridine and corresponding amide, amine, carbamate, thiol, or alcohol as leaving groups. The x-ray structure of a metagenomic CDA_F14 and the molecular modeling of the CDAs used in this study show a relationship between the bulkiness of a leaving group and the volume of the binding pocket, which is partly determined by the flexible β3α3 loop of CDAs. We propose that CDAs that are active toward a wide range of substrates participate in salvage and/or catabolism of variously modified pyrimidine nucleosides. This identified promiscuity of CDAs expands the knowledge about the cellular turnover of cytidine derivatives, including the pharmacokinetics of pyrimidine-based prodrugs.

m6A readers, writers, erasers, and the m6A epitranscriptome in breast cancer

Epitranscriptomic modification of RNA regulates human development, health, and disease. The true diversity of the transcriptome in breast cancer including chemical modification of transcribed RNA (epitranscriptomics) is not well understood due to limitations of technology and bioinformatic analysis. N-6-methyladenosine (m6A) is the most abundant epitranscriptomic modification of mRNA and regulates splicing, stability, translation, and intracellular localization of transcripts depending on m6A association with reader RNA-binding proteins. m6A methylation is catalyzed by the METTL3 complex and removed by specific m6A demethylase ALKBH5, with the role of FTO as an 'eraser' uncertain. In this review, we provide an overview of epitranscriptomics related to mRNA and focus on m6A in mRNA and its detection. We summarize current knowledge on altered levels of writers, readers, and erasers of m6A and their roles in breast cancer and their association with prognosis. We summarize studies identifying m6A peaks and sites in genes in breast cancer cells.

Nucleos'ID: A New Search Engine Enabling the Untargeted Identification of RNA Post-transcriptional Modifications from Tandem Mass Spectrometry Analyses of Nucleosides

As RNA post-transcriptional modifications are of growing interest, several methods were developed for their characterization. One of them established for their identification, at the nucleosidic level, is the hyphenation of separation methods, such as liquid chromatography or capillary electrophoresis, to tandem mass spectrometry. However, to our knowledge, no software is yet available for the untargeted identification of RNA post-transcriptional modifications from MS/MS data-dependent acquisitions. Thus, very long and tedious manual data interpretations are required. To meet the need of easier and faster data interpretation, a new user-friendly search engine, called Nucleos'ID, was developed for CE-MS/MS and LC-MS/MS users. Performances of this new software were evaluated on CE-MS/MS data from nucleoside analyses of already well-described Saccharomyces cerevisiae transfer RNA and Bos taurus total tRNA extract. All samples showed great true positive, true negative, and false discovery rates considering the database size containing all modified and unmodified nucleosides referenced in the literature. The true positive and true negative rates obtained were above 0.94, while the false discovery rates were between 0.09 and 0.17. To increase the level of sample complexity, untargeted identification of several RNA modifications from Pseudomonas aeruginosa 70S ribosome was achieved by the Nucleos'ID search following CE-MS/MS analysis.

Mitochondrial RNA methyltransferase TRMT61B is a new, potential biomarker and therapeutic target for highly aneuploid cancers

Despite being frequently observed in cancer cells, chromosomal instability (CIN) and its immediate consequence, aneuploidy, trigger adverse effects on cellular homeostasis that need to be overcome by anti-stress mechanisms. As such, these safeguard responses represent a tumor-specific Achilles heel, since CIN and aneuploidy are rarely observed in normal cells. Recent data have revealed that epitranscriptomic marks catalyzed by RNA-modifying enzymes change under various stress insults. However, whether aneuploidy is associated with such RNA modifying pathways remains to be determined. Through an in silico search for aneuploidy biomarkers in cancer cells, we found TRMT61B, a mitochondrial RNA methyltransferase enzyme, to be associated with high levels of aneuploidy. Accordingly, TRMT61B protein levels are increased in tumor cell lines with an imbalanced karyotype as well as in different tumor types when compared to control tissues. Interestingly, while TRMT61B depletion induces senescence in melanoma cell lines with low levels of aneuploidy, it leads to apoptosis in cells with high levels. The therapeutic potential of these results was further validated by targeting TRMT61B in transwell and xenografts assays. We show that TRM61B depletion reduces the expression of several mitochondrial encoded proteins and limits mitochondrial function. Taken together, these results identify a new biomarker of aneuploidy in cancer cells that could potentially be used to selectively target highly aneuploid tumors.

Data Analysis Pipeline for Detection and Quantification of Pseudouridine (ψ) in RNA by HydraPsiSeq

Pseudouridine, a modified RNA residue formed by the isomerization of its parental U nucleotide, is prevalent in a majority of cellular RNAs; its presence was reported in tRNA, rRNA, and sn/snoRNA as well as in mRNA/lncRNA. Multiple analytical deep sequencing-based approaches have been proposed for pseudouridine detection and quantification, among which the most popular relies on the use of soluble carbodiimide (termed CMCT). Recently, we developed an alternative protocol for pseudouridine mapping and quantification. The principle is based on protection of pseudouridine against random RNA cleavage by hydrazine/aniline treatment (HydraPsiSeq protocol). This "negative" detection mode requires higher sequencing depth and provides a precise quantification of the pseudouridine content. All "wet-lab" technical details of the HydraPsiSeq protocol have been described in recent publications. Here, we describe all bioinformatics analysis steps required for data processing from raw reads to the pseudouridylation profile of known or unknown RNA.

The role of post-transcriptional modifications during development

While the existence of post-transcriptional modifications of RNA nucleotides has been known for decades, in most RNA species the exact positions of these modifications and their physiological function have been elusive until recently. Technological advances, such as high-throughput next-generation sequencing (NGS) methods and nanopore-based mapping technologies, have made it possible to map the position of these modifications with single nucleotide accuracy, and genetic screens have uncovered the “writer”, “reader” and “eraser” proteins that help to install, interpret and remove such modifications, respectively. These discoveries led to intensive research programmes with the aim of uncovering the roles of these modifications during diverse biological processes. In this review, we assess novel discoveries related to the role of post-transcriptional modifications during animal development, highlighting how these discoveries can affect multiple aspects of development from fertilization to differentiation in many species.