Solution conformations of two naturally occurring RNA nucleosides: 3-methyluridine and 3-methylpseudouridine

The Ptch/SPOUT1 methyltransferase deposits an m3U modification on 28S rRNA for normal ribosomal function in flies and humans

The ribosomal RNA (rRNA) is one of the most heavily modified RNA species in nature. Although we have advanced knowledge of the sites, functions, and the enzymology of many of the rRNA modifications from all kingdoms of life, we lack basic understanding of many of those that are not universally present. A single N3 modified uridine base (m3U) was identified to be present on the 28S rRNA from humans and frogs but absent in bacteria or yeast. Here, we show that the equivalent m3U is present in Drosophila and that the Ptch/CG12128 enzyme and its human homolog SPOUT1 are both necessary and sufficient for carrying out the modification. The Ptch-modified U is at a functional center of the large ribosomal subunit, and, consistently, ptch-mutant cells suffer loss of ribosomal functions. SPOUT1, suggested to be the most druggable RNA methyltransferases in humans, represents a unique target where ribosomal functions could be specifically compromised in cancer cells.

Pseudouridine-modified RNA probe for label-free electrochemical detection of nucleic acids on 2D MoS2 nanosheets

RNA modification, particularly pseudouridine (Ψ), has played an important role in the development of the mRNA-based COVID-19 vaccine. This is because Ψ enhances RNA stability against nuclease activity and decreases the anti-RNA immune response. Ψ also provides structural flexibility to RNA by enhancing base stacking compared with canonical nucleobases. In this report, we demonstrate the first application of pseudouridine-modified RNA as a probe (Ψ-RNA) for label-free nucleic acid biosensing. It is known that MoS2 has a differential affinity for nucleic acids, which may be translated into a unique electronic signal. Herein, the Ψ-RNA probe interacts with the pristine MoS2 surface and causes a change in interfacial electrochemical charge transfer in the MoS2 nanosheets. Compared with an unmodified RNA probe, Ψ-RNA exhibited faster adsorption and higher affinity for MoS2. Moreover, Ψ-RNA could bind to complementary RNA and DNA targets with almost equal affinity when engaged with the MoS2 surface. Ψ-RNA maintained robust interactions with the MoS2 surface following the hybridization event, perhaps through its extra amino group. The detection sensitivity of the Ψ-RNA/MoS2 platform was as low as 500 attomoles, while the results also indicate that the probe can distinguish between complementary targets, single mismatches, and non-complementary nucleic acid sequences with statistical significance. This proof-of-concept study shows that the Ψ-RNA probe may solve numerous problems of adsorption-based biosensing platforms due to its stability and structural flexibility.

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.

The solid-state and solution conformations of 2'-deoxy-2'-fluoro-4'-substituted nucleosides

The conformational features of drugs are important with respect to their biological activity. In this report, we confirmed the solid-state conformation of 1-(4'-azido-2'-deoxy-2'-fluoro-β-d-arabinofuranosyl) cytosine using single-crystal X-ray crystallography and the conformations of three 2'-deoxy-2'-fluoro-4'-substituted nucleosides in solution using Altona-Haasnoot's equations and Nuclear Overhauser effect spectroscopy (NOESY). Furthermore, we compared the preferred solid-state and solution conformation of these nucleosides with thermodynamics cycles to obtain more evidence of their conformations. The results showed 1-(4'-azido-2'-deoxy-2'-fluoro-β-d-arabinofuranosyl) cytosine was south type conformation (C-3'-exo) in solid-state and three 4'-substituted nucleosides were north type conformations (C-3'-endo) in solution. The north type conformations in solution indicated these compounds were steady to acidic and enzymatic N-glycolysis.

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.

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.

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.

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.

2'-Methylseleno-modified oligoribonucleotides for X-ray crystallography synthesized by the ACE RNA solid-phase approach

Site-specifically modified 2'-methylseleno RNA represents a valuable derivative for phasing of X-ray crystallographic data. Several successful applications in three-dimensional structure determination of nucleic acids, such as the Diels-Alder ribozyme, have relied on this modification. Here, we introduce synthetic routes to 2'-methylseleno phosphoramidite building blocks of all four standard nucleosides, adenosine, cytidine, guanosine and uridine, that are tailored for 2'-O-bis(acetoxyethoxy)methyl (ACE) RNA solid-phase synthesis. We additionally report on their incorporation into oligoribonucleotides including deprotection and purification. The methodological expansion of 2'-methylseleno labeling via ACE RNA chemistry is a major step to make Se-RNA generally accessible and to receive broad dissemination of the Se-approach for crystallographic studies on RNA. Thus far, preparation of 2'-methylseleno-modified oligoribonucleotides has been restricted to the 2'-O-[(triisopropylsilyl)oxy]methyl (TOM) and 2'-O-tert-butyldimethylsilyl (TBDMS) RNA synthesis methods.

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

A series of 3-substituted uridine and pseudouridine derivatives, based on the naturally occurring 3-(3-amino-3-carboxypropyl) modification, were synthesized. Their aqueous solution conformations were determined by using circular dichroism and NMR spectroscopy. Functional group composition and chain length were shown to have only a subtle influence on the distribution of syn/anti conformations of the modified nucleosides. The dominating factor appears to be the glycosidic linkage (C- vs. N-glycoside) in determining the nucleoside conformation.