Solution NMR of the magnetically active fluorine-19 (19F) isotope offers tremendous advantages in the study of RNA structure and conformational changes. However, exploiting these advantages requires the synthesis of fluorinated RNA-building blocks, and this comes with certain challenges in overall yield of nucleotides, purity, and cost. Here, we detail a scalable and reliable method to produce [1′,5-13C2, 5-19F, 6-2H]- and [1′,5-13C2, 5-19F, 6-2H, 1,3-15N2]-5-fluorouracil-5′-triphosphate. We further present the steps employed to characterize the labeled nucleotide before its incorporation into RNAs with relevant biological activities.

19 F- 13 C spy couple enables NMR visualization of long-lived signals of large RNA structures, folds, and drug binding pockets.

Conformational dynamics of RNA molecules play a critical role in governing their biological functions. Measurements of RNA dynamic behavior sheds important light on sites that interact with their binding partners or cellular stimulators. However, such measurements using solution-state NMR are difficult for large RNA molecules (>70 nt; nt=nucleotides) owing to severe spectral overlap, homonuclear 13 C scalar couplings, and line broadening. Herein, a strategic combination of solid-phase synthesis, site-specific isotopic labeled phosphoramidites, and enzymatic ligation is introduced. This approach allowed the position-specific insertion of isotopic probes into a 96 nt CCR5 RNA fragment. Accurate measurements of functional dynamics using the Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion (RD) experiments enabled extraction of the exchange rates and populations of this RNA. NMR chemical shift perturbation analysis of the RNA/microRNA-1224 complex indicated that A90-C1' of the pseudoknot exhibits similar changes in chemical shift observed in the excited state. This work demonstrates the general applicability of a NMR-labeling strategy to probe functional RNA structural dynamics.

Climate change and human globalization have spurred the rapid spread of mosquito-borne diseases to naïve populations. One such emerging virus of public health concern is chikungunya virus (CHIKV), a member of the Togaviridae family, genus Alphavirus CHIKV pathogenesis is predominately characterized by acute febrile symptoms and severe arthralgia, which can persist in the host long after viral clearance. CHIKV has also been implicated in cases of acute encephalomyelitis, and its vertical transmission has been reported. Currently, no FDA-approved treatments exist for this virus. Recoding elements help expand the coding capacity in many viruses and therefore represent potential therapeutic targets in antiviral treatments. Here, we report the molecular and structural characterization of two CHIKV translational recoding signals: a termination codon read-through (TCR) element located between the nonstructural protein 3 and 4 genes and a programmed -1 ribosomal frameshift (-1 PRF) signal located toward the 3' end of the CHIKV 6K gene. Using Dual-Luciferase and immunoblot assays in HEK293T and U87MG mammalian cell lines, we validated and genetically characterized efficient TCR and -1 PRF. Analyses of RNA chemical modification data with selective 2'-hydroxyl acylation and primer extension (SHAPE) assays revealed that CHIKV -1 PRF is stimulated by a tightly structured, triple-stem hairpin element, consistent with previous observations in alphaviruses, and that the TCR signal is composed of a single large multibulged hairpin element. These findings illuminate the roles of RNA structure in translational recoding and provide critical information relevant for design of live-attenuated vaccines against CHIKV and related viruses.

The alphaviruses Venezuelan equine encephalitis virus (VEEV), eastern equine encephalitis virus (EEEV), and western equine encephalitis virus (WEEV) are arthropod-borne positive-strand RNA viruses that are capable of causing acute and fatal encephalitis in many mammals, including humans. VEEV was weaponized during the Cold War and is recognized as a select agent. Currently, there are no FDA-approved vaccines or therapeutics for these viruses. The spread of VEEV and other members of this family due to climate change-mediated vector range expansion underscores the need for research aimed at developing medical countermeasures. These viruses utilize programmed -1 ribosomal frameshifting (-1 PRF) to synthesize the viral trans-frame (TF) protein, which has previously been shown to be important for neuropathogenesis in the related Sindbis virus. Here, the alphavirus -1 PRF signals were characterized, revealing novel -1 PRF stimulatory structures. -1 PRF attenuation mildly affected the kinetics of VEEV accumulation in cultured cells but strongly inhibited its pathogenesis in an aerosol infection mouse model. Importantly, the decreased viral titers in the brains of mice infected with the mutant virus suggest that the alphavirus TF protein is important for passage through the blood-brain barrier and/or for neuroinvasiveness. These findings suggest a novel approach to the development of safe and effective live attenuated vaccines directed against VEEV and perhaps other closely related -1 PRF-utilizing viruses.IMPORTANCE Venezuelan equine encephalitis virus (VEEV) is a select agent that has been weaponized. This arthropod-borne positive-strand RNA virus causes acute and fatal encephalitis in many mammals, including humans. There is no vaccine or other approved therapeutic. VEEV and related alphaviruses utilize programmed -1 ribosomal frameshifting (-1 PRF) to synthesize the viral trans-frame (TF) protein, which is important for neuropathogenesis. -1 PRF attenuation strongly inhibited VEEV pathogenesis in mice, and viral replication analyses suggest that the TF protein is critical for neurological disease. These findings suggest a new approach to the development of safe and effective live attenuated vaccines directed against VEEV and other related viruses.

Our 13C- and 1H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.

The dependence of the conformation of the S-adenosylmethionine (SAM) II riboswitch on the concentration of added Mg(2+) ions and SAM, individually and in mixtures, was monitored by circular dichroism (CD) spectroscopy and by measurement of the diffusion coefficient. The results are analyzed in the context of two complementary quantitative models, both of which are consistent with a single underlying physical model. Magnesium binding sites in the open state have an affinity on average higher than the affinity of those in the compact state, but formation of the compact state is accompanied by an increase in the number of binding sites. Consequently, at low Mg(2+) concentrations, Mg(2+) binds preferentially to the open state, favoring its formation, but at high concentrations, Mg(2+) binds preferentially to the compact state. The affinity of the riboswitch for SAM increases drastically with an increased level of binding of Mg(2+) to the compact pseudoknot conformation. The effect of increasing concentrations of trimethylamine N-oxide (TMAO), a well-studied molecular crowding agent, on the conformation of the riboswitch and its affinity for SAM were also monitored by CD spectroscopy and measurement of diffusion. In the absence of added Mg(2+), high concentrations of TMAO were found to induce a conformational change compatible with the formation of the pseudoknot form but have only a small effect on the affinity of the RNA for SAM.

Conformational equilibria are increasingly recognized as pivotal for biological function. Traditional structural analyses provide a static image of conformers in solution that sometimes present conflicting views. From (13) C and (1) H chemical exchange saturation transfer experiments, in concert with ligation and selective labeling strategies, we show that in the absence of metabolite, a Mg(2+) (0-0.5 mm)-bound apo SAM-II riboswitch RNA exists in a minor (≈10 %) partially closed state that rapidly exchanges with a predominantly (≈90 %) open form with a lifetime of ≈32 ms. The base and sugar (H6,C6, H1',C1') chemical shifts of C43 for the dominant conformer are similar to those of a free CMP, but those of the minor apo species are comparable to shifts of CMPs in helical RNA regions. Our results suggest that these transient, low populated states stabilized by Mg(2+) will likely enhance rapid ligand recognition and, we anticipate, will play potentially ubiquitous roles in RNA signaling.

RNAs are an important class of cellular regulatory elements, and they are well characterized by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy in their folded or bound states. However, the apo or unfolded states are more difficult to characterize by either method. Particularly, effective NMR spectroscopy studies of RNAs in the past were hampered by chemical shift overlap of resonances and associated rapid signal loss due to line broadening for RNAs larger than the median size found in the PDB (~25 nt); most functional riboswitches are bigger than this median size. Incorporation of selective site-specific (13)C/(15)N-labeled nucleotides into RNAs promises to overcome this NMR size limitation. Unlike previous isotopic enrichment methods such as phosphoramidite, de novo, uniform-labeling, and selective-biomass approaches, this newer chemical-enzymatic selective method presents a number of advantages for producing labeled nucleotides over these other methods. For example, total chemical synthesis of nucleotides, followed by solid-phase synthesis of RNA using phosphoramidite chemistry, while versatile in incorporating isotope labels into RNA at any desired position, faces problems of low yields (<10%) that drop precipitously for oligonucleotides larger than 50 nt. The alternative method of de novo pyrimidine biosynthesis of NTPs is also a robust technique, with modest yields of up to 45%, but it comes at the cost of using 16 enzymes, expensive substrates, and difficulty in making some needed labeling patterns such as selective labeling of the ribose C1' and C5' and the pyrimidine nucleobase C2, C4, C5, or C6. Biomass-produced, uniformly or selectively labeled NTPs offer a third method, but suffer from low overall yield per labeled input metabolite and isotopic scrambling with only modest suppression of (13)C-(13)C couplings. In contrast to these four methods, our current chemo-enzymatic approach overcomes most of these shortcomings and allows for the synthesis of gram quantities of nucleotides with >80% yields while using a limited number of enzymes, six at most. The unavailability of selectively labeled ribose and base precursors had prevented the effective use of this versatile method until now. Recently, we combined an improved organic synthetic approach that selectively places (13)C/(15)N labels in the pyrimidine nucleobase (either (15)N1, (15)N3, (13)C2, (13)C4, (13)C5, or (13)C6 or any combination) with a very efficient enzymatic method to couple ribose with uracil to produce previously unattainable labeling patterns (Alvarado et al., 2014). Herein we provide detailed steps of both our chemo-enzymatic synthesis of custom nucleotides and their incorporation into RNAs with sizes ranging from 29 to 155 nt and showcase the dramatic improvement in spectral quality of reduced crowding and narrow linewidths. Applications of this selective labeling technology should prove valuable in overcoming two major obstacles, chemical shift overlap of resonances and associated rapid signal loss due to line broadening, that have impeded studying the structure and dynamics of large RNAs such as full-length riboswitches larger than the ~25 nt median size of RNA NMR structures found in the PDB.

C-di-GMP has emerged as a ubiquitous second messenger, which regulates the transition between sessile and motile lifestyles and virulence factor expression in many pathogenic bacteria using both RNA riboswitches and protein effectors. We recently showed that two additional class I c-di-GMP riboswitch aptamers (Ct-E88 and Cb-17B) bind c-di-GMP with nanomolar affinity, and that Ct-E88 RNA binds 2'-F-c-di-GMP 422 times less tightly than class I Vc2 RNA. Based on sequence comparison, it was concluded that the global folds of Ct-E88 and Vc2 RNAs were similar and that differences in ligand binding were probably due to differences in binding site architectures. Herein, we utilized EMSA, aptamer sensing spinach modules, SAXS and 1D NMR titration to study the conformational transitions of Ct-E88. We conclude that whereas the global folds of the bound states of Vc2 and Ct-E88 RNAs are similar, the unbound states are different and this could explain differences in ligand affinities between these class I c-di-GMP riboswitches.

Characterization of the structure and dynamics of nucleic acids by NMR benefits significantly from position specifically labeled nucleotides. Here an E. coli strain deficient in the transketolase gene (tktA) and grown on glucose that is labeled at different carbon sites is shown to facilitate cost-effective and large scale production of useful nucleotides. These nucleotides are site specifically labeled in C1′ and C5′ with minimal scrambling within the ribose ring. To demonstrate the utility of this labeling approach, the new site-specific labeled and the uniformly labeled nucleotides were used to synthesize a 36-nt RNA containing the catalytically essential domain 5 (D5) of the brown algae group II intron self-splicing ribozyme. The D5 RNA was used in binding and relaxation studies probed by NMR spectroscopy. Key nucleotides in the D5 RNA that are implicated in binding Mg2+ ions are well resolved. As a result, spectra obtained using selectively labeled nucleotides have higher signal-to-noise ratio compared to those obtained using uniformly labeled nucleotides. Thus, compared to the uniformly 13C/15N-labeled nucleotides, these specifically labeled nucleotides eliminate the extensive 13C–13C coupling within the nitrogenous base and ribose ring, give rise to less crowded and more resolved NMR spectra, and accurate relaxation rates without the need for constant-time or band-selective decoupled NMR experiments. These position selective labeled nucleotides should, therefore, find wide use in NMR analysis of biologically interesting RNA molecules.

Riboswitches are a newly discovered large family of structured functional RNA elements that specifically bind small molecule targets out of a myriad of cellular metabolites to modulate gene expression. Structural studies of ligand-bound riboswitches by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy have provided insights into detailed RNA-ligand recognition and interactions. However, the structures of ligand-free riboswitches remain poorly characterized. In this study, we have used a variety of biochemical, biophysical and computational techniques including small-angle X-ray scattering and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of SAM-II riboswitch. Our data demonstrate that the RNA adopts multiple conformations along its folding pathway and suggest that the RNA undergoes marked conformational changes upon Mg2+ compaction and S-adenosylmethionine (SAM) metabolite binding. Further studies indicated that Mg2+ ion is not essential for the ligand binding but can stabilize the complex by facilitating loop/stem interactions. In the presence of millimolar concentration of Mg2+ ion, the RNA samples a more compact conformation. This conformation is near to, but distinct from, the native fold and competent to bind the metabolite. We conclude that the formation of various secondary and tertiary structural elements, including a pseudoknot, occur to sequester the putative Shine-Dalgarno sequence of the RNA only after metabolite binding.

Ribsowitches are a newly discovered large family of structured RNA elements that are typically located at 5’ untranslated regions of messenger RNAs. They are composed of two domains: aptamer domain and expression platform. The aptamer domain functions as a sensor to recognize and bind specifically cellular metabolites while expression platform modulates gene expression in response to the conformational changes in the aptamer by either inhibiting transcription elongation or blocking translation. More than twenty riboswitch subfamilies have been identified so far in a variety of organisms ranging from bacteria, fungi to plants, indicating that riboswitch is a widespread genetic control element. Structural studies of ligand-bound riboswitches by X-ray crystallography and NMR spectroscopy provide insight into detailed RNA/ligand recognition and interactions. The structure, however, of ligand-free riboswitches remains poorly characterized. For better understanding the mechanism of riboswitches’ functions and their transition from ligand-free and ligand-bound forms, it is critical to characterize the molecular details of the unliganded state. We have employed a variety of biochemical, biophysical, and computational techniques including SAXS and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of riboswitches. Our data reveal that ligand binding of the RNAs causes significant conformational change and that only after ligand-binding does the formation of various secondary and tertiary structural elements occur.

Escherichia coli ( E. coli ) is an ideal organism to tailor-make labeled nucleotides for biophysical studies of RNA. Recently, we showed that adding labeled formate enhanced the isotopic enrichment at protonated carbon sites in nucleotides. In this paper, we show that growth of a mutant E. coli strain DL323 (lacking succinate and malate dehydrogenases) on 13 C-2-glycerol and 13 C-1,3-glycerol enables selective labeling at many useful sites for RNA NMR spectroscopy. For DL323 E. coli grown in 13 C-2-glycerol without labeled formate, all the ribose carbon atoms are labeled except the C3′ and C5′ carbon positions. Consequently the C1′, C2′ and C4′ positions remain singlet. In addition, only the pyrimidine base C6 atoms are substantially labeled to ~96% whereas the C2 and C8 atoms of purine are labeled to ~5%. Supplementing the growth media with 13 C-formate increases the labeling at C8 to ~88%, but not C2. Not unexpectedly, addition of exogenous formate is unnecessary for attaining the high enrichment levels of ~88% for the C2 and C8 purine positions in a 13 C-1,3-glycerol based growth. Furthermore, the ribose ring is labeled in all but the C4′ carbon position, such that the C2′ and C3′ positions suffer from multiplet splitting but the C5′ position remains singlet and the C1′ position shows a small amount of residual C1′–C2′ coupling. As expected, all the protonated base atoms, except C6, are labeled to ~90%. In addition, labeling with 13 C-1,3-glycerol affords an isolated methylene ribose with high enrichment at the C5′ position (~90%) that makes it particularly attractive for NMR applications involving CH 2 -TROSY modules without the need for decoupling the C4′ carbon. To simulate the tumbling of large RNA molecules, perdeuterated glycerol was added to a mixture of the four nucleotides, and the methylene TROSY experiment recorded at various temperatures. Even under conditions of slow tumbling, all the expected carbon correlations were observed, which indicates this approach of using nucleotides obtained from DL323 E. coli will be applicable to high molecular weight RNA systems.