MODOMICS: a database of RNA modification pathways. 2021 update.

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5'-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

Mapping RNA Modifications Using Photo-Crosslinking-Assisted Modification Sequencing.

Epitranscriptomic RNA modifications function as an important layer of gene regulation that modulates the function of RNA transcripts. A key step in understanding how RNA modifications regulate biological processes is the mapping of their locations, which is most commonly done by RNA immunoprecipitation (RIP) using modification-specific antibodies. Here, we describe the use of a photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) method, in conjunction with RNA modification-specific antibodies, to map modification sites. First described as photo-crosslinking-assisted m6A sequencing (PA-m6A-seq), this method allows the mapping of RNA modifications at a higher resolution, with lower background than traditional RIP, and can be adapted to any RNA modification for which a specific antibody is available.

Naturally occurring modified ribonucleosides.

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

The ribonucleoside AICAr induces differentiation of myeloid leukemia by activating the ATR/Chk1 via pyrimidine depletion.

Metabolic pathways play important roles in proliferation and differentiation of malignant cells. 5-Aminoimidazole-4-carboxamide ribonucleoside (AICAr), a precursor in purine biosynthesis and a well-established activator of AMP-activated protein kinase (AMPK), induces widespread metabolic alterations and is commonly used for dissecting the role of metabolism in cancer. We have previously reported that AICAr promotes differentiation and inhibits proliferation of myeloid leukemia cells. Here, using metabolic assays, immunoblotting, flow cytometry analyses, and siRNA-mediated gene silencing in leukemia cell lines, we show that AICAr-mediated differentiation was independent of the known metabolic effects of AMPK, including glucose consumption, but instead depends on the activation of the DNA damage-associated enzyme checkpoint kinase 1 (Chk1) induced by pyrimidine depletion. LC/MS/MS metabolomics analysis revealed that AICAr increases orotate levels and decreases uridine monophosphate (UMP) levels, consistent with inhibition of UMP synthesis at a step downstream of dihydroorotate dehydrogenase (DHODH). AICAr and the DHODH inhibitor brequinar had similar effects on differentiation markers and S-phase arrest, and genetic or pharmacological Chk1 inactivation abrogated both of these effects. Our results delineate an AMPK-independent effect of AICAr on myeloid leukemia differentiation that involves perturbation of pyrimidine biosynthesis and activation of the DNA damage response network.

Simple In Vitro Assay To Evaluate the Incorporation Efficiency of Ribonucleotide Analog 5'-Triphosphates into RNA by Human Mitochondrial DNA-Dependent RNA Polymerase.

There is a growing body of evidence suggesting that some ribonucleoside/ribonucleotide analogs may be incorporated into mitochondrial RNA by human mitochondrial DNA-dependent RNA polymerase (POLRMT) and disrupt mitochondrial RNA synthesis. An assessment of the incorporation efficiency of a ribonucleotide analog 5'-triphosphate by POLRMT may be used to evaluate the potential mitochondrial toxicity of the analog early in the development process. In this report, we provide a simple method to prepare active recombinant POLRMT. A robust in vitro nonradioactive primer extension assay was developed to assay the incorporation efficiency of ribonucleotide analog 5'-triphosphates. Our results show that many ribonucleotide analogs, including some antiviral compounds currently in various preclinical or clinical development stages, can be incorporated into newly synthesized RNA by POLRMT and that the incorporation of some of them can lead to chain termination. The discrimination (D) values of ribonucleotide analog 5'-triphosphates over those of natural ribonucleotide triphosphates (rNTPs) were measured to evaluate the incorporation efficiency of the ribonucleotide analog 5'-triphosphates by POLRMT. The discrimination values of natural rNTPs under the condition of misincorporation by POLRMT were used as a reference to evaluate the potential mitochondrial toxicity of ribonucleotide analogs. We propose the following criteria for the potential mitochondrial toxicity of ribonucleotide analogs based on D values: a safe compound has a D value of >105; a potentially toxic compound has a D value of >104 but <105; and a toxic compound has a D value of <104 This report provides a simple screening method that should assist investigators in designing ribonucleoside-based drugs having lower mitochondrial toxicity.

Chemical suppressors of mlo-mediated powdery mildew resistance.

Loss-of-function of barley mildew locus o (Mlo) confers durable broad-spectrum penetration resistance to the barley powdery mildew pathogen, Blumeria graminis f. sp. hordei (Bgh). Given the importance of mlo mutants in agriculture, surprisingly few molecular components have been identified to be required for this type of resistance in barley. With the aim to identify novel cellular factors contributing to mlo-based resistance, we devised a pharmacological inhibitor screen. Of the 41 rationally chosen compounds tested, five caused a partial suppression of mlo resistance in barley, indicated by increased levels of Bgh host cell entry. These chemicals comprise brefeldin A (BFA), 2',3'-dideoxyadenosine (DDA), 2-deoxy-d-glucose, spermidine, and 1-aminobenzotriazole. Further inhibitor analysis corroborated a key role for both anterograde and retrograde endomembrane trafficking in mlo resistance. In addition, all four ribonucleosides, some ribonucleoside derivatives, two of the five nucleobases (guanine and uracil), some guanine derivatives as well as various polyamines partially suppress mlo resistance in barley via yet unknown mechanisms. Most of the chemicals identified to be effective in partially relieving mlo resistance in barley also to some extent compromised powdery mildew resistance in an Arabidopsis mlo2 mlo6 double mutant. In summary, our study identified novel suppressors of mlo resistance that may serve as valuable probes to unravel further the molecular processes underlying this unusual type of disease resistance.

Transcriptome-wide Identification of RNA-binding Protein Binding Sites Using Photoactivatable-Ribonucleoside-Enhanced Crosslinking Immunoprecipitation (PAR-CLIP).

RNA-binding proteins (RBPs) mediate important co- and post-transcriptional gene regulation by binding pre-mRNA in a sequence- and/or structure-specific manner. For a comprehensive understanding of RBP function, transcriptome-wide mapping of the RNA-binding sites is essential, and CLIP-seq methods have been developed to elucidate protein/RNA interactions at high resolution. CLIP-seq combines protein/RNA UV-crosslinking with immunoprecipitation (CLIP) followed by high-throughput sequencing of crosslinked RNA fragments. To overcome the limitations of low RNA-protein crosslinking efficiency in standard CLIP-seq, photoactivatable-ribonucleoside-enhanced CLIP (PAR-CLIP) has been developed. Here, living cells or whole organisms are fed photo-activatable nucleoside analogs that are incorporated into nascent RNA transcripts before UV treatment. This allows greater crosslinking efficiency at comparable radiation doses for enhanced RNA recovery and separation of crosslinked target RNA fragments from background RNA degradation products. Moreover, it facilitates the generation of specific UV-induced mutations that mark the crosslinking nucleotide and allow transcriptome-wide identification of RBP binding sites at single-nucleotide resolution. © by 2017 John Wiley & Sons, Inc.

A Platform for Discovery and Quantification of Modified Ribonucleosides in RNA: Application to Stress-Induced Reprogramming of tRNA Modifications.

Here we describe an analytical platform for systems-level quantitative analysis of modified ribonucleosides in any RNA species, with a focus on stress-induced reprogramming of tRNA as part of a system of translational control of cell stress response. This chapter emphasizes strategies and caveats for each of the seven steps of the platform workflow: (1) RNA isolation, (2) RNA purification, (3) RNA hydrolysis to individual ribonucleosides, (4) chromatographic resolution of ribonucleosides, (5) identification of the full set of modified ribonucleosides, (6) mass spectrometric quantification of ribonucleosides, (6) interrogation of ribonucleoside datasets, and (7) mapping the location of stress-sensitive modifications in individual tRNA molecules. We have focused on the critical determinants of analytical sensitivity, specificity, precision, and accuracy in an effort to ensure the most biologically meaningful data on mechanisms of translational control of cell stress response. The methods described here should find wide use in virtually any analysis involving RNA modifications.

Mycobacterial RNA isolation optimized for non-coding RNA: high fidelity isolation of 5S rRNA from Mycobacterium bovis BCG reveals novel post-transcriptional processing and a complete spectrum of modified ribonucleosides.

A major challenge in the study of mycobacterial RNA biology is the lack of a comprehensive RNA isolation method that overcomes the unusual cell wall to faithfully yield the full spectrum of non-coding RNA (ncRNA) species. Here, we describe a simple and robust procedure optimized for the isolation of total ncRNA, including 5S, 16S and 23S ribosomal RNA (rRNA) and tRNA, from mycobacteria, using Mycobacterium bovis BCG to illustrate the method. Based on a combination of mechanical disruption and liquid and solid-phase technologies, the method produces all major species of ncRNA in high yield and with high integrity, enabling direct chemical and sequence analysis of the ncRNA species. The reproducibility of the method with BCG was evident in bioanalyzer electrophoretic analysis of isolated RNA, which revealed quantitatively significant differences in the ncRNA profiles of exponentially growing and non-replicating hypoxic bacilli. The method also overcame an historical inconsistency in 5S rRNA isolation, with direct sequencing revealing a novel post-transcriptional processing of 5S rRNA to its functional form and with chemical analysis revealing seven post-transcriptional ribonucleoside modifications in the 5S rRNA. This optimized RNA isolation procedure thus provides a means to more rigorously explore the biology of ncRNA species in mycobacteria.

A pre-steady state and steady state kinetic analysis of the N-ribosyl hydrolase activity of hCD157.

hCD157 catalyzes the hydrolysis of nicotinamide riboside (NR) and nicotinic acid riboside (NAR). The release of nicotinamide or nicotinic acid from NR or NAR was confirmed by spectrophotometric, HPLC and NMR analyses. hCD157 is inactivated by a mechanism-based inhibitor, 2'-deoxy-2'-fluoro-nicotinamide arabinoside (fNR). Modification of the enzyme during the catalytic cycle by NR, NAR, or fNR increased the intrinsic protein fluorescence by approximately 50%. Pre-steady state and steady state data were used to derive a minimal kinetic scheme for the hydrolysis of NR. After initial complex formation a reversible step (360 and 30s(-1)) is followed by a slow irreversible step (0.1s(-1)) that defined the rate limiting step, or kcat. The calculated KMapp value for NR in the hydrolytic reaction is 6nM. The values of the kinetic constants suggest that one biological function of cell-surface hCD157 is to bind and slowly hydrolyze NR, possibly converting it to a ligand-activated receptor. Differences in substrate preference between hCD157 and hCD38 were rationalized through a comparison of the crystal structures of the two proteins. This comparison identified several residues in hCD157 (F108 and F173) that can potentially hinder the binding of dinucleotide substrates (NAD+).

Quantitative analysis of ribonucleoside modifications in tRNA by HPLC-coupled mass spectrometry.

Post-transcriptional modification of RNA is an important determinant of RNA quality control, translational efficiency, RNA-protein interactions and stress response. This is illustrated by the observation of toxicant-specific changes in the spectrum of tRNA modifications in a stress-response mechanism involving selective translation of codon-biased mRNA for crucial proteins. To facilitate systems-level studies of RNA modifications, we developed a liquid chromatography-mass spectrometry (LC-MS) technique for the quantitative analysis of modified ribonucleosides in tRNA. The protocol includes tRNA purification by HPLC, enzymatic hydrolysis, reversed-phase HPLC resolution of the ribonucleosides, and identification and quantification of individual ribonucleosides by LC-MS via dynamic multiple reaction monitoring (DMRM). In this approach, the relative proportions of modified ribonucleosides are quantified in several micrograms of tRNA in a 15-min LC-MS run. This protocol can be modified to analyze other types of RNA by modifying the steps for RNA purification as appropriate. By comparison, traditional methods for detecting modified ribonucleosides are labor- and time-intensive, they require larger RNA quantities, they are modification-specific or require radioactive labeling.

Developmental disruption of the serotonin system alters circadian rhythms.

Serotonin (5-HT) plays an important role in circadian rhythms, acting to modulate photic input to the mammalian clock, the suprachiasmatic nucleus (SCN), as well as playing a role in non-photic input. The transcription factor Pet-1 is an early developmental indicator of neurons that are destined for a 5-HTergic fate. Mice lacking the Pet-1 gene show a 70% loss of 5-HT immunopositive cell bodies in adult animals. 5-HT neurotoxic lesion studies using 5,7-dihydroxytryptamine (5,7-DHT) have highlighted species-specific differences in response to 5-HT depletion and studies using knockout mice lacking various 5-HT receptors have helped to elucidate the role of individual 5-HT receptors in mediating 5-HT's effects on circadian rhythms. Here we investigate the effects of a developmental disruption of the 5-HT system on the SCN and circadian wheel-running behavior. Immunohistochemical analysis confirmed depletion of 5-HT fiber innervation to the SCN as well as greatly reduced numbers of cell bodies in the raphe nuclei in Pet-1 knockout mice. These mice also display significantly longer free-running periods than wildtype or heterozygote counterparts. In light-dark cycles, knockouts showed a shift in peak wheel running behavior towards the late night as compared to wildtype and heterozygote animals. When kept in constant darkness for 70 days, wildtype animals showed decreases in free-running period over time while the period of knockout animals remained constant. Immunohistochemical analysis for neuropeptides within the SCN indicates that the behavioral changes observed in Pet-1 knockout mice were not due to gross changes in SCN structure. These results suggest that developmental loss of serotonergic input to the clock has long-term consequences for both circadian clock parameters and the temporal organization of activity.

Engineering ribonucleoside triphosphate specificity in a thymidylyltransferase.

Nature's glycosylation catalysts, glycosyltransferases, indirectly manipulate and control many important biological processes by transferring sugar nucleotide donors onto acceptors. Challenging chemical synthesis impedes synthetic access to sugar nucleotides and limits the study of many glycosyltransferases. Enzymatic access to sugar nucleotides is a rapidly expanding avenue of research, limited only by the substrate specificity of the enzyme. We have explored the promiscuous thymidylyltransferase from Streptococcus pneumoniae, Cps2L, and enhanced its uridylyltransferase and guanidyltransferase activities by active site engineering. Mutagenesis at position Q24 resulted in a variant with 10-, 3-, and 2-fold enhancement of UDP-glucosamine, UDP-mannose, and UDP- N-acetylglucosamine production, respectively. New catalytic activities were observed for the Cps2L variant over the wild-type enzyme, including the formation of GDP-mannose. The variant was evaluated as a catalyst for the formation of a series of dTDP- and UDP-furanoses and notably produced dTDP-Gal f in 90% yield and UDP-Ara f in 30% yield after 12 h. A series of 3- O-alkylglucose 1-phosphates were also evaluated as substrates, and notable conversions to UDP-3- O-methylglucose and UDP-3- O-dodecylglucose were achieved with the variant but not the wild-type enzyme. The Q24S variant also enhanced essentially all thymidylyltransferase activities relative to the wild-type enzyme. Comparison of active sites of uridylyltransferases and thymidylyltransferases with products bound indicate the Q24S variant to be a new approach in broadening nucleotidylyltransferase activity.

Histone hyperacetylation during SV40 transcription is regulated by p300 and RNA polymerase II translocation.

The effects of the RNA polymerase II (RNAPII) translocation inhibitors alpha amanitin and 5,6-dichloro-1-beta-D-ribobenzimidazole (DRB) and an siRNA targeting p300 on the presence of RNAPII, p300, hyperacetylated H4 and H3 and unmodified H4 and H3 in transcribing simian virus 40 (SV40) minichromosomes were determined. Following treatment with alpha amanitin we observed a profound reduction in the occupancy of the promoter by RNAPII, the loss of p300 from chromatin fragments containing RNAPII, and an increase in the amount of unmodified H4 and H3 associated with the RNAPII. Treatment with DRB had little effect on the presence of RNAPII or p300 but also resulted in a significant increase in the amount of unmodified H4 and H3 present in chromatin fragments associated with RNAPII. Following treatment with a p300 small interfering RNA (siRNA), we observed a significant decrease in late transcription and a corresponding reduction in the amounts of p300 and hyperacetylated histones associated with the transcribing SV40 minichromosomes. We conclude that in transcribing SV40 minichromosomes histone hyperacetylation and deacetylation is dependent upon the presence of p300 and an as yet unknown histone deacetylase associated with the RNAPII complex that occurs coordinately as the RNAPII complex moves through a nucleosome.

Elongation factor Ts can act as a steric chaperone by increasing the solubility of nucleotide binding-impaired elongation factor-Tu.

Several elongation factor (EF) Tu mutants (T25A, H22Y/T25S, D80N, D138N) that have impaired nucleotide binding show decreased solubility on overexpression in the E. coli cell, an indication that they do not fold correctly. Moreover, EF-Tu[T25A] and EF-Tu[D80N] were shown to inhibit cell growth on expression, an effect attributed to their sequestration of EF-Ts [Krab, I. M., and Parmeggiani, A. (1999) J. Biol. Chem. 274, 11132--11138; Krab, I. M., and Parmeggiani, A. (1999) Biochemistry 38, 13035--13041]. We present here results showing that the co-overexpression of EF-Ts at a 1:1 ratio dramatically improves the solubility of mutant EF-Tu, although in the case of EF-Tu[D138N]--which cannot at all bind the nucleotides available in the cell--this is a slow process. Moreover, with co-overexpression of EF-Ts, the mentioned growth inhibition is relieved. We conclude that for the formation of a correct EF-Tu structure the nucleotide plays an important role as a "folding nucleus", and also that in its absence EF-Ts can act as a folding template or steric chaperone for the correct folding of EF-Tu.

Excavating an active site: the nucleobase specificity of ribonuclease A.

Ribonuclease A (RNase A) catalyzes the cleavage of RNA after pyrimidine nucleotides. When bound in the active site, the base of a pyrimidine nucleotide forms hydrogen bonds with the side chain of Thr45. Here, the role of Thr45 was probed by using the wild-type enzyme, its T45G variant, X-ray diffraction analysis, and synthetic oligonucleotides as ligands and substrates. Catalytic specificity was determined with the fluorogenic substrate: 6-carboxyfluorescein approximately dArXdAdA approximately 6-carboxytetramethylrhodamine (6-FAM approximately dArXdAdA approximately 6-TAMRA), where X = C, U, A, or G. Wild-type RNase A cleaves 10(6)-fold faster when X = C than when X = A. Likewise, its affinity for the non-hydrolyzable oligonucleotide 6-FAM approximately d(CAA) is 50-fold greater than for 6-FAM approximately d(AAA). T45G RNase A cleaves 6-FAM approximately dArAdAdA approximately 6-TAMRA 10(2)-fold faster than does the wild-type enzyme. The structure of crystalline T45G RNase A, determined at 1.8-A resolution by X-ray diffraction analysis, does not reveal new potential interactions with a nucleobase. Indeed, the two enzymes have a similar affinity for 6-FAM approximately d(AAA). The importance of pentofuranosyl ring conformation to nucleotide specificity was probed with 6-FAM approximately d(AU(F)AA), where U(F) is 2'-deoxy-2'-fluorouridine. The conformation of the pentofuranosyl ring in dU(F) is known to be more similar to that in rU than dU. The affinity of wild-type RNase A for 6-FAM approximately d(AU(F)AA) is 50-fold lower than for 6-FAM approximately d(AUAA). This discrimination is lost in the T45G enzyme. Together, these data indicate that the side chain of Thr45 plays multiple roles-interacting favorably with pyrimidine nucleobases but unfavorably with purine nucleobases. Moreover, a ribose-like ring disfavors the interaction of Thr45 with a pyrimidine nucleobase, suggesting that Thr45 enhances catalysis by ground-state destabilization.

Site-specific mutagenesis of mistletoe lectin: the role of RIP activity in apoptosis.

Recombinant mistletoe lectin (rML) belongs to the class of type II ribosome-inactivating proteins (RIP) composed of a catalytically active A-chain with rRNA N-glycosidase activity and a B-chain with carbohydrate binding properties. To investigate the contribution of the enzymatic activity of the rML A-chain to the observed cytotoxic and apoptotic effects, an rMLA E166Q R169Q molecule was developed by means of site-specific mutagenesis. Following heterologous expression, the activity of mutant rMLA was measured in a cell-free assay for rRNA-N-glycosidase activity. Moreover, after generation of heterodimer, the activities of mutant rML E166Q R169Q and rML wild type were determined in a cytotoxicity and apoptosis assay. Although the reduction of activity as measured in the cell-free RIP assay was more pronounced (factor 237) than in both cellular assays (factors 20-22), the data clearly indicate a close correlation between cytotoxicity, apoptosis, and the enzymatic activity of the rML A-chain. Thus, RIP activity is an essential feature of rML and therefore a prerequisite for its biological function as an anticancer agent.