Synthesis and solution conformation studies of 3-substituted uridine and pseudouridine derivatives

acp³U: A Conserved RNA Modification with Lessons Yet to Unfold

RNA modifications are highly conserved across all domains of life, suggesting an early emergence and a fundamental role in cellular processes. The modification 3-(3-amino-3-carboxypropyl)uridine (acp³U) is found in tRNAs of eukaryotes and prokaryotes, and in the 16S rRNA of archaea. In eukaryotic rRNA, a complex modification containing the acp group, m1acp3Ψ is present at the analogous position. Although this modification was first identified in tRNA in 1969, only recently have the enzymes responsible for the synthesis of this modification on tRNA been identified. Despite its deep evolutionary conservation, the biological role of acp³U on tRNAs remains elusive. In Escherichia coli, it may contribute to genomic stability, while in human cells, loss of both tRNA acp³U-modifying enzymes impairs cell growth, though the underlying mechanisms are not yet understood. The conservation and multifunctionality of acp³U highlight the broader challenges of elucidating the roles of tRNA modifications in cellular homeostasis.

The modified RNA base acp3U is an attachment site for N-glycans in glycoRNA

GlycoRNA consists of RNAs modified with secretory N-glycans that are presented on the cell surface. Although previous work supported a covalent linkage between RNA and glycans, the direct chemical nature of the RNA-glycan connection was not described. Here, we develop a sensitive and scalable protocol to detect and characterize native glycoRNAs. Leveraging RNA-optimized periodate oxidation and aldehyde ligation (rPAL) and sequential window acquisition of all theoretical mass spectra (SWATH-MS), we identified the modified RNA base 3-(3-amino-3-carboxypropyl)uridine (acp3U) as a site of attachment of N-glycans in glycoRNA. rPAL offers sensitivity and robustness as an approach for characterizing direct glycan-RNA linkages occurring in cells, and its flexibility will enable further exploration of glycoRNA biology.

Systematic Engineering of Escherichia coli for Efficient Production of Pseudouridine from Glucose and Uracil

In this study, we proposed a biological approach to efficiently produce pseudouridine (Ψ) from glucose and uracil in vivo using engineered Escherichia coli. By screening host strains and core enzymes, E. coli MG1655 overexpressing Ψ monophosphate (ΨMP) glycosidase and ΨMP phosphatase was obtained, which displayed the highest Ψ concentration. Then, optimization of the RBS sequences, enhancement of ribose 5-phosphate supply in the cells, and overexpression of the membrane transport protein UraA were investigated. Finally, fed-batch fermentation of Ψ in a 5 L fermentor can reach 27.5 g/L with a yield of 89.2 mol % toward uracil and 25.6 mol % toward glucose within 48 h, both of which are the highest to date. In addition, the Ψ product with a high purity of 99.8% can be purified from the fermentation broth after crystallization. This work provides an efficient and environmentally friendly protocol for allowing for the possibility of Ψ bioproduction on an industrial scale.

Nanopore Direct RNA Sequencing for Modified Uridine Nucleotides Yields Signals Dependent on the Physical Properties of the Modified Base

Sequencing for RNA modifications with the nanopore direct RNA sequencing platform provides ionic current levels, helicase dwell times, and base call data that differentiate the modifications from the canonical form. Herein, model RNAs were synthesized with site-specific uridine (U) base modifications that enable the study of increasing an alkyl group size, halogen identity, or a change in base acidity to impact the nanopore data. The analysis concluded that increases in alkyl size trend with greater current blockage but a similar change in base-call error was not found. The addition of a halogen series to C5 of U revealed that the current levels recorded a trend with the water-octanol partition coefficient of the base, as well as the base call error. Studies with U modifications that are deprotonated (i.e., anionic) under the sequencing conditions gave broad current levels that influenced the base call error. Some modifications led to helicase dwell time changes. These insights provide design parameters for modification-specific chemical reagents that can shift nanopore signatures to minimize false positive reads, a known issue with this sequencing approach.

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA

Nucleophilic addition of bisulfite to pyrimidine bases has been known for a half century, and the reaction has been in use for at least a quarter of a century for identifying 5-methylcytidine in DNA. This account focuses on the chemistry of bisulfite with pseudouridine, an isomer of the RNA nucleoside uridine in which the uracil base is connected to C1' of ribose via C5 instead of N1. Pseudouridine, Ψ, is the most common nucleotide modification found in cellular RNA overall, in part due to its abundance in rRNAs and tRNAs. It has a stabilizing influence on RNA structure because N1 is now available for additional hydrogen bonding and because the heterocycle is slightly better at π stacking. The isomerization of U to Ψ in RNA strands is catalyzed by 13 different enzymes in humans and 11 in E. coli; some of these enzymes are implicated in disease states which is testament to the biological importance of pseudouridine in cells. Recently, pseudouridine came into the limelight as the key modification that, after N1 methylation, enables mRNA vaccines to be delivered efficiently into human tissue with minimal generation of a deleterious immunogenic response. Here we describe the bisulfite reaction with pseudouridine which gives rise to a chemical sequencing method to map the modified base in the epitranscriptome. Unlike the reaction with cytidine, the addition of bisulfite to Ψ leads irreversibly to form an adduct that is bypassed during cDNA synthesis by reverse transcriptases yielding a characteristic deletion signature. Although there were hints to the structure of the bisulfite adduct(s) 30 to 50 years ago, it took modern spectroscopic and computational methods to solve the mystery. Raman spectroscopy along with extensive NMR, ECD, and computational work led to the assignment of the major product as the (R) diastereomer of an oxygen adduct at C1' of a ring-opened pseudouridine. Mechanistically, this arose from a succession of conjugate addition, E2 elimination, and a [2,3] sigmatropic rearrangement, all of which are stereodefined reactions. A minor reaction with excess bisulfite led to the (S) isomer of a S-adducted SO3- group. Understanding structure and mechanism aided the design of a Ψ-specific sequencing reaction and guided attempts to improve the utility and specificity of the method. Separately, we have been investigating the use of nanopore direct RNA sequencing, a single-molecule method that directly analyzes RNA strands isolated from cells after end-ligation of adaptor sequences. By combining the electrical current and base-calling data from the nanopore with dwell-time analysis from the helicase employed to deliver RNA to the nanopore, we were able to map Ψ sites in nearly all sequence contexts. This analysis was employed to find Ψ residues in the SARS-CoV-2 vRNA, to analyze the sequence context effects of mRNA vaccine synthesis via in vitro transcription, and to evaluate the impact of stress on chemical modifications in the E. coli ribosome. Most recently, we found that bisulfite treatment of RNA leading to Ψ adducts could modulate the nanopore signal to help in mapping modifications of low occupancy.

A selective and atom-economic rearrangement of uridine by cascade biocatalysis for production of pseudouridine

As a crucial factor of their therapeutic efficacy, the currently marketed mRNA vaccines feature uniform substitution of uridine (U) by the corresponding C-nucleoside, pseudouridine (Ψ), in 1-N-methylated form. Synthetic supply of the mRNA building block (1-N-Me-Ψ-5'-triphosphate) involves expedient access to Ψ as the principal challenge. Here, we show selective and atom-economic 1N-5C rearrangement of β-D-ribosyl on uracil to obtain Ψ from unprotected U in quantitative yield. One-pot cascade transformation of U in four enzyme-catalyzed steps, via D-ribose (Rib)-1-phosphate, Rib-5-phosphate (Rib5P) and Ψ-5'-phosphate (ΨMP), gives Ψ. Coordinated function of the coupled enzymes in the overall rearrangement necessitates specific release of phosphate from the ΨMP, but not from the intermediary ribose phosphates. Discovery of Yjjg as ΨMP-specific phosphatase enables internally controlled regeneration of phosphate as catalytic reagent. With driving force provided from the net N-C rearrangement, the optimized U reaction yields a supersaturated product solution (∼250 g/L) from which the pure Ψ crystallizes (90% recovery). Scale up to 25 g isolated product at enzyme turnovers of ∼10⁵mol/mol demonstrates a robust process technology, promising for Ψ production. Our study identifies a multistep rearrangement reaction, realized by cascade biocatalysis, for C-nucleoside synthesis in high efficiency.

5-Ethynyluridine: A Bio-orthogonal Uridine Variant for mRNA-Based Therapies and Vaccines

The identification of pseudo- and N¹ -methylpseudo-uridine (Ψ and mΨ, respectively) as immunosilent uridine analogues has propelled the development of mRNA-based vaccines and therapeutics. Here, we have characterised another uridine analogue, 5-ethynyluridine (EU), which has an ethynyl moiety. We show that this uridine analogue does not cause immune activation in human macrophages, as it does not induce interleukin-6 secretion or expression of the inflammatory and antiviral genes MX1, PKR, and TAP2. Moreover, EU allows for prolonged expression, as shown with mRNA coding for yellow fluorescent protein (YFP). Side-by-side comparisons of EU with unmodified, Ψ, and mΨ revealed that EU-modified mRNA is expressed at lower levels, but confers similar stability and low immunogenicity to the other uridine analogues. Furthermore, structure analysis of modified mRNAs suggests that the observed phenotype is largely independent of RNA folding. Thus, EU is a potential candidate for RNA-based vaccines and therapeutics.

Nanopore sequencing for N1-methylpseudouridine in RNA reveals sequence-dependent discrimination of the modified nucleotide triphosphate during transcription

Direct RNA sequencing with a commercial nanopore platform was used to sequence RNA containing uridine (U), pseudouridine (Ψ) or N1-methylpseudouridine (m1Ψ) in >100 different 5-nucleotide contexts. The base calling data for Ψ or m1Ψ were similar but different from U allowing their detection. Understanding the nanopore signatures for Ψ and m1Ψ enabled a running start T7 RNA polymerase assay to study the selection of UTP versus ΨTP or m1ΨTP competing mixtures in all possible adjacent sequence contexts. A significant sequence context dependency was observed for T7 RNA polymerase with insertion yields for ΨTP versus UTP spanning a range of 20-65%, and m1ΨTP versus UTP producing variable yields that differ by 15-70%. Experiments with SP6 RNA polymerase, as well as chemically-modified triphosphates and DNA templates provide insight to explain the observations. The SP6 polymerase introduced m1ΨTP when competed with UTP with a smaller window of yields (15-30%) across all sequence contexts studied. These results may aid in future efforts that employ RNA polymerases to make therapeutic mRNAs with sub-stoichiometric amounts of m1Ψ.

Pseudouridine and N1-Methylpseudouridine Display pH-Independent Reaction Rates with Bisulfite Yielding Ribose Adducts

In RNA, pseudouridine (Ψ) and 5-methylcytidine (m5C) are located by their differential reactions with NaHSO3 at pH 5. The pyrimidines were allowed to react with NaHSO3, NaN3, NaCN, or NaSCN at pH 5 to find that NaHSO3 was unique in achieving quantitative yields. Pseudouridine reaction selectivity with NaHSO3 was found at pH 7 supported by the reaction rate constants. The Ψ derivative N1-methylpseudouridine found in mRNA vaccines reacts similarly with bisulfite to yield ribose adducts.

Stronger together for in-cell translation: natural and unnatural base modified mRNA

The preparation of highly modified mRNAs and visualization of their cellular distribution are challenging. We report in-cell application of in vitro transcribed mRNA containing natural base modifications and site-specifically introduced artificial nucleotides. Click chemistry on mRNA allows visualization in cells with excellent signal intensities. While non-specific introduction of reporter groups often leads to loss in mRNA functionality, we combined the benefits from site-specificity in the 3′-UTR incorporated unnatural nucleotides with the improved translation efficiency of the natural base modifications Ψ and 5mC. A series of experiments is described to observe, quantify and verify mRNA functionality. This approach represents a new way to visualize mRNA delivery into cells and monitor its spread on a cellular level and translation efficiency. We observed increased protein expression from this twofold chemically modified, artificial mRNA counterbalancing a reduced transfection rate. This synergetic effect can be exploited as a powerful tool for future research on mRNA therapeutics.

Introducing unnatural base modifications site-specifically into the 3′-UTR of an mRNA bearing natural base modifications allows efficient visualization in cells by click chemistry. An enhanced protein expression in cells is observed from this twofold modified mRNA.

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

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

Modifications in an Emergency: The Role of N1-Methylpseudouridine in COVID-19 Vaccines

The novel coronavirus SARS-CoV-2, the cause of the COVID-19 pandemic, has inspired one of the most efficient vaccine development campaigns in human history. A key aspect of COVID-19 mRNA vaccines is the use of the modified nucleobase N1-methylpseudouridine (m1Ψ) to increase their effectiveness. In this Outlook, we summarize the development and function of m1Ψ in synthetic mRNAs. By demystifying how a novel element within these medicines works, we aim to foster understanding and highlight future opportunities for chemical innovation.

Semi-enzymatic synthesis of pseudouridine

Modifications of RNA molecules have a significant effect on their structure and function. One of the most common modifications is the isomerization from uridine to pseudouridine. Despite its prevalence in natural RNA sequences, organic synthesis of pseudouridine has been challenging because of the stereochemistry requirement and the sensitivity of reaction steps to moisture. Herein, a semi-enzymatic synthetic route is developed for the synthesis of pseudouridine using adenosine 5'-monophosphate and uracil as the starting materials and a reverse reaction catalyzed by the pseudouridine monophosphate glycosidase. This synthetic route has only three steps and the overall yield of β-pseudouridine production was 68.4%.

A Small Molecule that Binds an RNA Repeat Expansion Stimulates Its Decay via the Exosome Complex

Many diseases are caused by toxic RNA repeats. Herein, we designed a lead small molecule that binds the structure of the r(CUG) repeat expansion [r(CUG)^exp] that causes myotonic dystrophy type 1 (DM1) and Fuchs endothelial corneal dystrophy (FECD) and rescues disease biology in patient-derived cells and in vivo. Interestingly, the compound's downstream effects are different in the two diseases, owing to the location of the repeat expansion. In DM1, r(CUG)^exp is harbored in the 3' untranslated region, and the compound has no effect on the mRNA's abundance. In FECD, however, r(CUG)^exp is located in an intron, and the small molecule facilitates excision of the intron, which is then degraded by the RNA exosome complex. Thus, structure-specific, RNA-targeting small molecules can act disease specifically to affect biology, either by disabling the gain-of-function mechanism (DM1) or by stimulating quality control pathways to rid a disease-affected cell of a toxic RNA (FECD).

Molecular Dynamics Simulation of the Conformational Preferences of Pseudouridine Derivatives: Improving the Distribution in the Glycosidic Torsion Space

There are only four derivatives of pseudouridine (Ψ) that are known to occur naturally in RNA as post-transcriptional modifications. We have studied the conformational consequences of pseudouridylation and further modifications using replica exchange molecular dynamics simulations at the nucleoside level, and the simulated conformational preferences were compared with the available experimental (NMR) data. We found that the existing AMBER FF99-derived parameters for these nucleosides did not reproduce the observed experimental features and while the recommended bsc0 correction could be combined with these parameters leading to an improvement in the description of sugar pucker distributions, the χ_OL3 correction could not be applied to these nucleosides as such because of base isomerization. On the other hand, the revised χ torsion parameters (χ_IDRP) for Ψ developed earlier by us (Deb, I., J. Comput. Chem., 2016, 37, 1576-1588) in combination with the AMBER provided parameters and the revised γ torsion parameters generated conformational distributions, which generally were in better agreement with the experimental data. A significant shift of the distribution of base orientation toward the syn conformation was observed with our revised parameter sets compared to the large excess of anti conformation predicted by the FF99 parameters. Overall, our observations indicated that our revised set of parameters (χ_IDRP) for Ψ were also able to generate conformational distributions for all of the derivatives of Ψ in better agreement with the experimental data.

Biogenesis and functions of aminocarboxypropyluridine in tRNA

Transfer (t)RNAs contain a wide variety of post-transcriptional modifications, which play critical roles in tRNA stability and functions. 3-(3-amino-3-carboxypropyl)uridine (acp3U) is a highly conserved modification found in variable- and D-loops of tRNAs. Biogenesis and functions of acp3U have not been extensively investigated. Using a reverse-genetic approach supported by comparative genomics, we find here that the Escherichia coli yfiP gene, which we rename tapT (tRNA aminocarboxypropyltransferase), is responsible for acp3U formation in tRNA. Recombinant TapT synthesizes acp3U at position 47 of tRNAs in the presence of S-adenosylmethionine. Biochemical experiments reveal that acp3U47 confers thermal stability on tRNA. Curiously, the ΔtapT strain exhibits genome instability under continuous heat stress. We also find that the human homologs of tapT, DTWD1 and DTWD2, are responsible for acp3U formation at positions 20 and 20a of tRNAs, respectively. Double knockout cells of DTWD1 and DTWD2 exhibit growth retardation, indicating that acp3U is physiologically important in mammals.

Total synthesis of pseudouridineviaHeck-typeC-glycosylation

Pseudouridine (1) was synthesized by functional group interconversions of the Heck adduct11from 2,4-dimethoxy-5-iodopyrimidine (8) and ribofuranoid glycal4.

Additive CHARMM force field for naturally occurring modified ribonucleotides

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

Time-dependent thermocontrol of the hydrophilic and lipophilic properties of DNA oligonucleotide prodrugs

This unit describes the preparation of alkylthioalkylated and formamidoalkylated alcohols, an amidoalkylated alcohol, a hydroxylalkylated phosphoramidate, and their phosphoramidothioate derivatives, all of which have been identified as heat-sensitive thiophosphate-protecting groups in the development of thermolytic immunostimulatory DNA prodrugs. The alcohols are converted to their deoxyribonucleoside phosphoramidite derivatives, which are then used in the preparation of thermosensitive dinucleoside phosphorothioates. The thiophosphate-protecting groups of these dinucleoside phosphorothioates presumably undergo thermolytic cyclodeesterification at elevated temperature under essentially neutral conditions to release the desired phosphorothioate diester function. On the basis of their thermolytic deprotection kinetics, one can identify those thiophosphate-protecting groups that (i) may be useful for thiophosphate protection of CpG motifs of immunostimulatory DNA oligonucleotides (CpG ODNs); (ii) are suitable for protection of phosphodiester functions flanking the CpG motifs; and (iii) offer adequate protection of terminal phosphodiester functions against ubiquitous extracellular and intracellular exonucleases that may be found in biological environments.

Reparameterization of RNA chi Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine

A reparameterization of the torsional parameters for the glycosidic dihedral angle, chi, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99chi. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99chi force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99chi force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) (1)H, steady-state 1D (1)H nuclear Overhauser effect (NOE), and transient 1D (1)H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2'-endo sugar puckering of the pyrimidines, while the AMBER99chi force field's predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310 degrees for the base orientation of purines. The AMBER99chi force field prefers anti conformations around 185 degrees , which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99chi force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.

Cytoplasmic delivery of liposomal contents mediated by an acid-labile cholesterol-vinyl ether-PEG conjugate

An acid-cleavable PEG lipid, 1'-(4'-cholesteryloxy-3'-butenyl)-omega-methoxy-polyethylene[112] glycolate (CVEP), has been developed that produces stable liposomes when dispersed as a minor component (0.5-5 mol %) in 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). Cleavage of CVEP at mildly acidic pHs results in dePEGylation of the latently fusogenic DOPE liposomes, thereby triggering the onset of content release. This paper describes the synthesis of CVEP via a six-step sequence starting from the readily available precursors 1,4-butanediol, cholesterol, and mPEG acid. The hydrolysis rates and release kinetics from CVEP/DOPE liposome dispersions as a function of CVEP loading, as well as the cryogenic transmission electron microscopy and pH-dependent monolayer properties of 9:91 CVEP/DOPE mixtures, also are reported. When folate receptor-positive KB cells were exposed to calcein-loaded 5:95 CVEP/DOPE liposomes containing 0.1 mol % folate-modified 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-polyethylene[76] glycolamide (folate-PEG-DSPE), delivery of the calcein cargo to the cytoplasm of the cells was observed as determined by fluorescence microscopy and flow cytometry. Fluorescence resonance energy transfer analysis of lipid mixing in these cells was consistent with membrane-membrane fusion between the liposome and endosomal membranes.

Synthesis and solution conformation studies of the modified nucleoside N(4),2'-O-dimethylcytidine (m(4)Cm) and its analogues

The dimethylated ribosomal nucleoside m(4)Cm and its monomethylated analogues Cm and m(4)C were synthesized. The conformations (syn vs anti) of the three modified nucleosides and cytidine were determined by CD and 1D NOE difference spectroscopy. The ribose sugar puckers were determined by the use of proton coupling constants. The position of modification (e.g., O vs N methylation) was found to have an effect on the sugar pucker of cytidine.