Pseudoknot preorganization of the preQ1 class I riboswitch

Flipping the script: Understanding riboswitches from an alternative perspective

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Advances in RNA Labeling with Trifluoromethyl Groups

Fluorine labeling of ribonucleic acids (RNA) in conjunction with 19 F NMR spectroscopy has emerged as a powerful strategy for spectroscopic analysis of RNA structure and dynamics, and RNA-ligand interactions. This study presents the first syntheses of 2'-OCF3 guanosine and uridine phosphoramidites, their incorporation into oligoribonucleotides by solid-phase synthesis and a comprehensive study of their properties. NMR spectroscopic analysis showed that the 2'-OCF3 modification is associated with preferential C2'-endo conformation of the U and G ribose in single-stranded RNA. When paired to the complementary strand, slight destabilization of the duplex caused by the modification was revealed by UV melting curve analysis. Moreover, the power of the 2'-OCF3 label for NMR spectroscopy is demonstrated by dissecting RNA pseudoknot folding and its binding to a small molecule. Furthermore, the 2'-OCF3 modification has potential for applications in therapeutic oligonucleotides. To this end, three 2'-OCF3 modified siRNAs were tested in silencing of the BASP1 gene which indicated enhanced performance for one of them. Importantly, together with earlier work, the present study completes the set of 2'-OCF3 nucleoside phosphoramidites to all four standard nucleobases (A, U, C, G) and hence enables applications that utilize the favorable properties of the 2'-OCF3 group without any restrictions in placing the modification into the RNA target sequence.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.

A chemical probe based on the PreQ1 metabolite enables transcriptome-wide mapping of binding sites

The role of metabolite-responsive riboswitches in regulating gene expression in bacteria is well known and makes them useful systems for the study of RNA-small molecule interactions. Here, we study the PreQ₁ riboswitch system, assessing sixteen diverse PreQ₁-derived probes for their ability to selectively modify the class-I PreQ₁ riboswitch aptamer covalently. For the most active probe (11), a diazirine-based photocrosslinking analog of PreQ₁, X-ray crystallography and gel-based competition assays demonstrated the mode of binding of the ligand to the aptamer, and functional assays demonstrated that the probe retains activity against the full riboswitch. Transcriptome-wide mapping using Chem-CLIP revealed a highly selective interaction between the bacterial aptamer and the probe. In addition, a small number of RNA targets in endogenous human transcripts were found to bind specifically to 11, providing evidence for candidate PreQ₁ aptamers in human RNA. This work demonstrates a stark influence of linker chemistry and structure on the ability of molecules to crosslink RNA, reveals that the PreQ₁ aptamer/ligand pair are broadly useful for chemical biology applications, and provides insights into how PreQ₁, which is similar in structure to guanine, interacts with human RNAs.

Microsecond Folding of preQ1 Riboswitch and Its Biological Significance Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

Riboswitches are regulatory elements of bacterial mRNA which function with conformational switching upon binding of specific cellular metabolites. In particular, transcriptional riboswitches regulate gene expression kinetically through the conformational change of the aptamer domain. In this study, we investigate the conformational dynamics and ligand binding mechanisms of the aptamer domain of a transcriptional prequeuosine (preQ₁) riboswitch from Bacillus subtilis using two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) with microsecond time resolution. The obtained time-resolved single-molecule data indicate that the aptamer domain undergoes folding/unfolding including three forms, which are attributed to hairpin (O), pseudoknot-like (pF), and H-type pseudoknot (fF) structures. It is found that a cofactor, Mg²⁺, binds only to the fF form with the conformational selection mechanism. In contrast, it is indicated that the ligand, preQ₁, binds to the O form with the induced-fit mechanism and significantly accelerates the microsecond O → pF folding process. It is also shown that the binding with preQ₁ substantially stabilizes the fF form that is generated from the pF form with a long time constant (>10 ms). Combining these results with the results of a former smFRET study on the slower time scale, we obtain an overall picture of the folding/unfolding dynamics of the aptamer domain as well as its energy landscape. On the basis of the picture obtained, we discuss the significance of the microsecond folding/unfolding of the aptamer domain for biological function of the riboswitch and propose the molecular mechanism of the gene expression controlled by the structural dynamics of the aptamer domain.

2'-O-Trifluoromethylated RNA - a powerful modification for RNA chemistry and NMR spectroscopy

New RNA modifications are needed to advance our toolbox for targeted manipulation of RNA. In particular, the development of high-performance reporter groups facilitating spectroscopic analysis of RNA structure and dynamics, and of RNA-ligand interactions has attracted considerable interest. To this end, fluorine labeling in conjunction with 19F-NMR spectroscopy has emerged as a powerful strategy. Appropriate probes for RNA previously focused on single fluorine atoms attached to the 5-position of pyrimidine nucleobases or at the ribose 2'-position. To increase NMR sensitivity, trifluoromethyl labeling approaches have been developed, with the ribose 2'-SCF3 modification being the most prominent one. A major drawback of the 2'-SCF3 group, however, is its strong impact on RNA base pairing stability. Interestingly, RNA containing the structurally related 2'-OCF3 modification has not yet been reported. Therefore, we set out to overcome the synthetic challenges toward 2'-OCF3 labeled RNA and to investigate the impact of this modification. We present the syntheses of 2'-OCF3 adenosine and cytidine phosphoramidites and their incorporation into oligoribonucleotides by solid-phase synthesis. Importantly, it turns out that the 2'-OCF3 group has only a slight destabilizing effect when located in double helical regions which is consistent with the preferential C3'-endo conformation of the 2'-OCF3 ribose as reflected in the 3 J (H1'-H2') coupling constants. Furthermore, we demonstrate the exceptionally high sensitivity of the new label in 19F-NMR analysis of RNA structure equilibria and of RNA-small molecule interactions. The study is complemented by a crystal structure at 0.9 Å resolution of a 27 nt hairpin RNA containing a single 2'-OCF3 group that well integrates into the minor groove. The new label carries high potential to outcompete currently applied fluorine labels for nucleic acid NMR spectroscopy because of its significantly advanced performance.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ₁ riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ₁ riboswitch aptamers without or with a ligand bound. We show that the preQ₁ riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3'-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ₁ binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

High-affinity recognition of specific tRNAs by an mRNA anticodon-binding groove

T-box riboswitches are modular bacterial noncoding RNAs that sense and regulate amino acid availability through direct interactions with tRNAs. Between the 5' anticodon-binding stem I domain and the 3' amino acid sensing domains of most T-boxes lies the stem II domain of unknown structure and function. Here, we report a 2.8-Å cocrystal structure of the Nocardia farcinica ileS T-box in complex with its cognate tRNA^Ile. The structure reveals a perpendicularly arranged ultrashort stem I containing a K-turn and an elongated stem II bearing an S-turn. Both stems rest against a compact pseudoknot, dock via an extended ribose zipper and jointly create a binding groove specific to the anticodon of its cognate tRNA. Contrary to proposed distal contacts to the tRNA elbow region, stem II locally reinforces the codon-anticodon interactions between stem I and tRNA, achieving low-nanomolar affinity. This study illustrates how mRNA junctions can create specific binding sites for interacting RNAs of prescribed sequence and structure.

Ligand Modulates Cross-Coupling between Riboswitch Folding and Transcriptional Pausing

Numerous classes of riboswitches have been found to regulate bacterial gene expression in response to physiological cues, offering new paths to antibacterial drugs. As common studies of isolated riboswitches lack the functional context of the transcription machinery, we here combine single-molecule, biochemical, and simulation approaches to investigate the coupling between co-transcriptional folding of the pseudoknot-structured preQ₁ riboswitch and RNA polymerase (RNAP) pausing. We show that pausing at a site immediately downstream of the riboswitch requires a ligand-free pseudoknot in the nascent RNA, a precisely spaced sequence resembling the pause consensus, and electrostatic and steric interactions with the RNAP exit channel. While interactions with RNAP stabilize the native fold of the riboswitch, binding of the ligand signals RNAP release from the pause. Our results demonstrate that the nascent riboswitch and its ligand actively modulate the function of RNAP and vice versa, a paradigm likely to apply to other cellular RNA transcripts.

High-Yield Spin Labeling of Long RNAs for Electron Paramagnetic Resonance Spectroscopy

Site-directed spin labeling is a powerful tool for investigating the conformation and dynamics of biomacromolecules such as RNA. Here we introduce a spin labeling strategy based on click chemistry in solution that, in combination with enzymatic ligation, allows highly efficient labeling of complex and long RNAs with short reaction times and suppressed RNA degradation. With this approach, a 34-nucleotide aptamer domain of the preQ1 riboswitch and an 81-nucleotide TPP riboswitch aptamer could be labeled with two labels in several positions. We then show that conformations of the preQ1 aptamer and its dynamics can be monitored in the absence and presence of Mg²⁺ and a preQ1 ligand by continuous wave electron paramagnetic resonance spectroscopy at room temperature and pulsed electron-electron double resonance spectroscopy (PELDOR or DEER) in the frozen state.

Structural Insight into the Unbound State of the DNA Analogue of the PreQ1 Riboswitch: A Thermodynamic Approach

The preQ₁ riboswitch aptamer domain is very dynamic in its unbound state with the ability to form multiple structures: a hairpin, kissing hairpins, and pseudoknot-like structure. The aim of this study is to determine whether the DNA analogue (PreQ₁) is able to form structures similar to that of the reported RNA aptamer. Using a thermodynamic approach, we report on structural determination using differential scanning calorimetry under different salt conditions. Further analysis of the primary sequence allowed us to design modified molecules to determine what potential structures are forming in this single-stranded DNA analogue. We found, in a 16 mM Na⁺ solution, PreQ₁ has three transitions with T_M values of 14.8, 19.4, and 26.2 °C and a total ΔH of -44.7 kcal/mol. With the increase in salt concentration to 116 mM, there are T_M values of 22.3, 28.7, and 38.9 °C and a ΔH of -69.1 kcal/mol, while at 216 mM, the three transitions have T_M values of 24.4, 31.6, and 42.9 °C with a total ΔH of -71.5 kcal/mol. Therefore, the increase in enthalpy is due to the formation of additional base-pair stacks. The modified molecules, which would inhibit pseudoknot formation, kissing hairpins, and internal loop interactions, were fully characterized and compared to the native DNA analogue. The analysis of the enthalpy and differential binding of counterions allows us to conclude this single-stranded DNA analogue under physiological conditions is not forming a pseudoknot-like structure. Instead, two potential structures, Compact-Hairpin and Kissing-Complex, are more likely and could be in equilibrium.

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

Allosteric pathways in tetrahydrofolate sensing riboswitch with dynamics correlation network

Riboswitches are cis-acting genetic control elements. Due to their fundamental importance in bacteria gene regulation, they have been proposed as antibacterial drug targets. Tetrahydrofolate (THF) is an essential cofactor of one-carbon transfer reactions and downregulates the expression of downstream genes. However, information on how to transfer from the binding site of THF to the expression platform is still unavailable. Herein, a nucleotide/nucleotide dynamics correlation network based on an all-atom molecular dynamic simulation was used to reveal the regulation mechanism of a THF-sensing riboswitch. Shortest pathway analysis based on the network illustrates that there is an allosteric pathway through the P2 helix to the pseudoknot, then to the P1 helix in the THF-riboswitch. Thus the hypothesis of "THF-binding induced allosteric switching" was proposed and evaluated using THF and pseudoknot weakened experiments. Furthermore, a possible allosteric pathway of C30-C31-G33-A34-G35-G36-G37-A38-G48-G47-U46-A90-U91-C92-G93-C94-G95-C96 was identified and confirmed through the perturbation of the network. The proposed allosteric mechanism and the underlying allosteric pathway provide fundamental insights for the regulation of THF sensing riboswitches.

The synthesis of 15N(7)-Hoogsteen face-labeled adenosine phosphoramidite for solid-phase RNA synthesis

ABSTRACT

We have developed an efficient route for the synthesis of 15N(7)-labeled adenosine as phosphoramidite building block for site- and atom-specific incorporation into RNA by automated solid-phase synthesis. Such labeled RNA is required for the evaluation of selected non-canonical base pair interactions in folded RNA using NMR spectroscopic methods.

GRAPHICAL ABSTRACT

Modulation of Conformational Equilibria in the S-Adenosylmethionine (SAM) II Riboswitch by SAM, Mg(2+), and Trimethylamine N-Oxide

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.

Mn(2+)-sensing mechanisms of yybP-ykoY orphan riboswitches

Gene regulation in cis by riboswitches is prevalent in bacteria. The yybP-ykoY riboswitch family is quite widespread, yet its ligand and function remained unknown. Here, we characterize the Lactococcus lactis yybP-ykoY orphan riboswitch as a Mn(2+)-dependent transcription-ON riboswitch, with a ∼30-40 μM affinity for Mn(2+). We further determined its crystal structure at 2.7 Å to elucidate the metal sensing mechanism. The riboswitch resembles a hairpin, with two coaxially stacked helices tethered by a four-way junction and a tertiary docking interface. The Mn(2+)-sensing region, strategically located at the highly conserved docking interface, has two metal binding sites. Whereas one site tolerates the binding of either Mg(2+) or Mn(2+), the other site strongly prefers Mn(2+) due to a direct contact from the N7 of an invariable adenosine. Mutagenesis and a Mn(2+)-free E. coli yybP-ykoY structure further reveal that Mn(2+) binding is coupled with stabilization of the Mn(2+)-sensing region and the aptamer domain.

Conformational Rearrangements of Individual Nucleotides during RNA-Ligand Binding Are Rate-Differentiated

A pronounced rate differentiation has been found for conformational rearrangements of individual nucleobases that occur during ligand recognition of the preQ1 class-I riboswitch aptamer from Thermoanaerobacter tengcongensis. Rate measurements rely on the 2ApFold approach by analyzing the fluorescence response of riboswitch variants, each with a single, strategically positioned 2-aminopurine nucleobase substitution. Observed rate discrimination between the fastest and the slowest conformational adaption is 22-fold, with the largest rate observed for the rearrangement of a nucleoside directly at the binding site and the smallest rate observed for the 3'-unpaired nucleoside that stacks onto the pseudo-knot-closing Watson-Crick base pair. Our findings provide novel insights into how compact, prefolded RNAs that follow the induced-fit recognition mechanism adapt local structural elements in response to ligand binding on a rather broad time scale and how this process culminates in a structural signal that is responsible for efficient downregulation of ribosomal translation.

Mg(2+) shifts ligand-mediated folding of a riboswitch from induced-fit to conformational selection

Bacterial riboswitches couple small-molecule ligand binding to RNA conformational changes that widely regulate gene expression, rendering them potential targets for antibiotic intervention. Despite structural insights, the ligand-mediated folding mechanisms of riboswitches are still poorly understood. Using single-molecule fluorescence resonance energy transfer (smFRET), we have investigated the folding mechanism of an H-type pseudoknotted preQ1 riboswitch in dependence of Mg(2+) and three ligands of distinct affinities. We show that, in the absence of Mg(2+), both weakly and strongly bound ligands promote pseudoknot docking through an induced-fit mechanism. By contrast, addition of as low as 10 μM Mg(2+) generally shifts docking toward conformational selection by stabilizing a folded-like conformation prior to ligand binding. Supporting evidence from transition-state analysis further highlights the particular importance of stacking interactions during induced-fit and of specific hydrogen bonds during conformational selection. Our mechanistic dissection provides unprecedented insights into the intricate synergy between ligand- and Mg(2+)-mediated RNA folding.

The "Speedy" Synthesis of Atom-Specific (15)N Imino/Amido-Labeled RNA

Although numerous reports on the synthesis of atom-specific (15)N-labeled nucleosides exist, fast and facile access to the corresponding phosphoramidites for RNA solid-phase synthesis is still lacking. This situation represents a severe bottleneck for NMR spectroscopic investigations on functional RNAs. Here, we present optimized procedures to speed up the synthesis of (15)N(1) adenosine and (15)N(1) guanosine amidites, which are the much needed counterparts of the more straightforward-to-achieve (15)N(3) uridine and (15)N(3) cytidine amidites in order to tap full potential of (1)H/(15)N/(15)N-COSY experiments for directly monitoring individual Watson-Crick base pairs in RNA. Demonstrated for two preQ1 riboswitch systems, we exemplify a versatile concept for individual base-pair labeling in the analysis of conformationally flexible RNAs when competing structures and conformational dynamics are encountered.

Native purification and labeling of RNA for single molecule fluorescence studies

The recent discovery that non-coding RNAs are considerably more abundant and serve a much wider range of critical cellular functions than recognized over previous decades of research into molecular biology has sparked a renewed interest in the study of structure-function relationships of RNA. To perform their functions in the cell, RNAs must dominantly adopt their native conformations, avoiding deep, non-productive kinetic traps that may exist along a frustrated (rugged) folding free energy landscape. Intracellularly, RNAs are synthesized by RNA polymerase and fold co-transcriptionally starting from the 5' end, sometimes with the aid of protein chaperones. By contrast, in the laboratory RNAs are commonly generated by in vitro transcription or chemical synthesis, followed by purification in a manner that includes the use of high concentrations of urea, heat and UV light (for detection), resulting in the denaturation and subsequent refolding of the entire RNA. Recent studies into the nature of heterogeneous RNA populations resulting from this process have underscored the need for non-denaturing (native) purification methods that maintain the co-transcriptional fold of an RNA. Here, we present protocols for the native purification of an RNA after its in vitro transcription and for fluorophore and biotin labeling methods designed to preserve its native conformation for use in single molecule fluorescence resonance energy transfer (smFRET) inquiries into its structure and function. Finally, we present methods for taking smFRET data and for analyzing them, as well as a description of plausible overall preparation schemes for the plethora of non-coding RNAs.

Expanding the Scope of 2'-SCF3 Modified RNA

The 2′-trifluoromethylthio (2′-SCF₃) modification endows ribonucleic acids with exceptional properties and has attracted considerable interest as a reporter group for NMR spectroscopic applications. However, only modified pyrimidine nucleosides have been generated so far. Here, the syntheses of 2′-SCF₃ adenosine and guanosine phosphoramidites of which the latter was obtained in highly efficient manner by an unconventional Boc-protecting group strategy, are reported. RNA solid-phase synthesis provided site-specifically 2′-SCF₃-modified oligoribonucleotides that were investigated intensively. Their excellent behavior in ¹⁹F NMR spectroscopic probing of RNA ligand binding was exemplified for a noncovalent small molecule–RNA interaction. Moreover, comparably to the 2′-SCF₃ pyrimidine nucleosides, the purine counterparts were also found to cause a significant thermodynamic destabilization when located in double helical regions. This property was considered beneficial for siRNA design under the aspect to minimize off-target effects and their performance in silencing of the BASP1 gene was demonstrated.

Free energy landscape and multiple folding pathways of an H-type RNA pseudoknot

How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.

Characterizing excited conformational states of RNA by NMR spectroscopy

Conformational dynamics is a hallmark of diverse non-coding RNA functions. During these functional processes, RNA molecules almost ubiquitously undergo conformational transitions that are tuned to meet distinct structural and kinetic requirements for proper function. A complete mechanistic understanding of RNA function requires comprehensive structural and dynamic knowledge of these complex transitions, which often involve alternative higher-energy conformational states that pose a major challenge for high-resolution structural study by conventional methods. In this review, we describe recent progress in RNA NMR that has started to unveil detailed structural, thermodynamic and kinetic insights into some of these excited conformational states of RNA and their functional roles in biology.

Polyfluorinated groups in medicinal chemistry

Introduction of novel and diverse functional groups in drug discovery is always seen with hesitancy until good activity and low toxicity characteristics are proven. The introduction of fluorine in drug-like compounds is now a well-accepted strategy in medicinal chemistry. However, polyfluoroalkyl groups, with the exception of trifluoromethyl substituents, are not well explored yet. Our aim is to show to the readers how polyfluorinated groups can be beneficial to the properties of pharmaceutically active compounds by highlighting the structure–activity relationship (SAR) studies that led to the selection of polyfluorinated moieties as key structural features. Despite the fact that the use of higher polyfluoroalkyl/aryl moieties is still in its infancy, we believe that they will soon acquire the same importance of their lower parents.

(Dis)similar Analogues of Riboswitch Metabolites as Antibacterial Lead Compounds

The rise of antimicrobial resistance in human pathogenic bacteria has increased the necessity for the discovery of novel, yet unexplored antibacterial drug targets. Riboswitches, which are embedded in untranslated regions of bacterial messenger RNA (mRNA), represent such an interesting target structure. These RNA elements regulate gene expression upon binding to natural metabolites, second messengers, and inorganic ions, such as fluoride with high affinity and in a highly discriminative manner. Recently, efforts have been directed toward the identification of artificial riboswitch activators by establishing high-throughput screening assays, fragment-based screening, and structure-guided ligand design approaches. Emphasis in this review is placed on the special requirements and synthesis of new potential antibiotic drugs that target riboswitches in which dissimilarity is an important aspect in the design of potential lead compounds.

Ligand-detected relaxation dispersion NMR spectroscopy: dynamics of preQ1-RNA binding

An NMR-based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand-detected Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments is described. A (15)N-modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand-receptor interaction: the off rate and the population of the small molecule-receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch's function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand-detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.

Syntheses of (15)N-labeled pre-queuosine nucleobase derivatives

Pre-queuosine or queuine (preQ1) is a guanine derivative that is involved in the biosynthetic pathway of the hypermodified tRNA nucleoside queuosine (Que). The core structure of preQ1 is represented by 7-(aminomethyl)-7-deazaguanine (preQ1 base). Here, we report the synthesis of three preQ1 base derivatives with complementary (15)N-labeling patterns, utilizing [(15)N]-KCN, [(15)N]-phthalimide, and [(15)N3]-guanidine as cost-affordable (15)N sources. Such derivatives are required to explore the binding process of the preQ1 base to RNA targets using advanced NMR spectroscopic methods. PreQ1 base specifically binds to bacterial mRNA domains and thereby regulates genes that are required for queuosine biosynthesis.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

Structure and function of pseudoknots involved in gene expression control

Natural RNA molecules can have a high degree of structural complexity but even the most complexly folded RNAs are assembled from simple structural building blocks. Among the simplest RNA elements are double-stranded helices that participate in the formation of different folding topologies and constitute the major fraction of RNA structures. One common folding motif of RNA is a pseudoknot, defined as a bipartite helical structure formed by base-pairing of the apical loop in the stem-loop structure with an outside sequence. Pseudoknots constitute integral parts of the RNA structures essential for various cellular activities. Among many functions of pseudoknotted RNAs is feedback regulation of gene expression, carried out through specific recognition of various molecules. Pseudoknotted RNAs autoregulate ribosomal and phage protein genes in response to downstream encoded proteins, while many metabolic and transport genes are controlled by cellular metabolites interacting with pseudoknotted RNA elements from the riboswitch family. Modulation of some genes also depends on metabolite-induced messenger RNA (mRNA) cleavage performed by pseudoknotted ribozymes. Several regulatory pseudoknots have been characterized biochemically and structurally in great detail. These studies have demonstrated a plethora of pseudoknot-based folds and have begun uncovering diverse molecular principles of the ligand-dependent gene expression control. The pseudoknot-mediated mechanisms of gene control and many unexpected and interesting features of the regulatory pseudoknots have significantly advanced our understanding of the genetic circuits and laid the foundation for modulation of their outcomes.

A neutral pH thermal hydrolysis method for quantification of structured RNAs

Riboswitch aptamers adopt diverse and complex tertiary structural folds that contain both single-stranded and double-stranded regions. We observe that this high degree of secondary structure leads to an appreciable hypochromicity that is not accounted for in the standard method to calculate extinction coefficients using nearest-neighbor effects, which results in a systematic underestimation of RNA concentrations. Here we present a practical method for quantifying riboswitch RNAs using thermal hydrolysis to generate the corresponding pool of mononucleotides, for which precise extinction coefficients have been measured. Thermal hydrolysis can be performed at neutral pH without reaction quenching, avoids the use of nucleases or expensive fluorescent dyes, and does not require generation of calibration curves. The accuracy of this method for determining RNA concentrations has been validated using quantitative (31)P-NMR calibrated to an external standard. We expect that this simple procedure will be generally useful for the accurate quantification of any sequence-defined RNA sample, which is often a critical parameter for in vitro binding and kinetic assays.

Structure and function of preQ1 riboswitches

PreQ₁ riboswitches help regulate the biosynthesis and transport of preQ₁ (7-aminomethyl-7-deazaguanine), a precursor of the hypermodified guanine nucleotide queuosine (Q), in a number of Firmicutes, Proteobacteria, and Fusobacteria. Queuosine is almost universally found at the wobble position of the anticodon in asparaginyl, tyrosyl, histidyl and aspartyl tRNAs, where it contributes to translational fidelity. Two classes of preQ₁ riboswitches have been identified (preQ₁-I and preQ₁-II), and structures of examples from both classes have been determined. Both classes form H-type pseudoknots upon preQ₁ binding, each of which has distinct unusual features and modes of preQ₁ recognition. These features include an unusually long loop 2 in preQ₁-I pseudoknots and an embedded hairpin in loop 3 in preQ₁-II pseudoknots. PreQ₁-I riboswitches are also notable for their unusually small aptamer domain, which has been extensively investigated by NMR, X-ray crystallography, FRET, and other biophysical methods. Here we review the discovery, structural biology, ligand specificity, cation interactions, folding, dynamics, and applications to biotechnology of preQ₁ riboswitches. This article is part of a Special Issue entitled: Riboswitches.

Global analysis of riboswitches by small-angle X-ray scattering and calorimetry

Riboswitches are phylogenetically widespread non-coding mRNA domains that directly bind cellular metabolites and regulate transcription, translation, RNA stability or splicing via alternative RNA structures modulated by ligand binding. The details of ligand recognition by many riboswitches have been elucidated using X-ray crystallography and NMR. However, the global dynamics of riboswitch-ligand interactions and their thermodynamic driving forces are less understood. By compiling the work of many laboratories investigating riboswitches using small-angle X-ray scattering (SAXS) and isothermal titration calorimetry (ITC), we uncover general trends and common themes. There is a pressing need for community-wide consensus experimental conditions to allow results of riboswitch studies to be compared rigorously. Nonetheless, our meta-analysis reveals considerable diversity in the extent to which ligand binding reorganizes global riboswitch structures. It also demonstrates a wide spectrum of enthalpy-entropy compensation regimes across riboswitches that bind a diverse set of ligands, giving rise to a relatively narrow range of physiologically relevant free energies and ligand affinities. From the strongly entropy-driven binding of glycine to the predominantly enthalpy-driven binding of c-di-GMP to their respective riboswitches, these distinct thermodynamic signatures reflect the versatile strategies employed by RNA to adapt to the chemical natures of diverse ligands. Riboswitches have evolved to use a combination of long-range tertiary interactions, conformational selection, and induced fit to work with distinct ligand structure, charge, and solvation properties. This article is part of a Special Issue entitled: Riboswitches.

Insights into ligand binding to PreQ1 Riboswitch Aptamer from molecular dynamics simulations

Riboswitches play roles in transcriptional or translational regulation through specific ligand binding of their aptamer domains. Although a number of ligand-bound aptamer complex structures have been solved, it is important to know ligand-free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, preQ₁ riboswitch aptamer domain from Bacillus subtilis is studied by overall 1.5 μs all-atom molecular dynamics simulations We found that the ligand-free aptamer has a stable state with a folded P1-L3 and open binding pocket. The latter forms a cytosine-rich pool in which the nucleotide C19 oscillates between close and open positions, making it a potential conformation for preQ₁ entrance. The dynamic picture further suggests that the specific recognition of preQ₁ by the aptamer domain is not only facilitated by the key nucleotide C19 but also aided and enhanced by other cytosines around the binding pocket. These results should help to understand the details of preQ₁ binding.

Riboswitch structure and dynamics by smFRET microscopy

Riboswitches are structured noncoding RNA elements that control the expression of their embedding messenger RNAs by sensing the intracellular concentration of diverse metabolites. As the name suggests, riboswitches are dynamic in nature so that studying their inherent conformational dynamics and ligand-mediated folding is important for understanding their mechanism of action. Single-molecule fluorescence energy transfer (smFRET) microscopy is a powerful and versatile technique for studying the folding pathways and intra- and intermolecular dynamics of biological macromolecules, especially RNA. The ability of smFRET to monitor intramolecular distances and their temporal evolution make it a particularly insightful tool for probing the structure and dynamics of riboswitches. Here, we detail the general steps for using prism-based total internal reflection fluorescence microscopy for smFRET studies of the structure, dynamics, and ligand-binding mechanisms of riboswitches.

A novel paramagnetic relaxation enhancement tag for nucleic acids: a tool to study structure and dynamics of RNA

In this work, we present a novel 2,2,6,6-tetramethylpiperidine 1-oxyl (TEMPO) radical phosphoramidite building block, which can be attached to the 5'-terminus of nucleic acids. To investigate the paramagnetic relaxation enhancement (PRE) emanating from this radical center, we incorporated the TEMPO label into various types of RNAs. We measured proton PREs for selectively (13)C-isotope labeled nucleotides to derive long-range distance restraints in a short 15 nucleotide stem-loop model system, underscoring the potential of the 5'-TEMPO tag to determine long-range distance restraints for solution structure determination. We subsequently applied the distance-dependent relaxation enhancement induced by the nitroxide radical to discern two folding states in a bistable RNA. Finally, we investigated the fast conformational sampling of the HIV-1 TAR RNA, a paradigm for structural flexibility in nucleic acids. With PRE NMR in combination with molecular dynamics simulations, the structural plasticity of this RNA was analyzed in the absence and presence of the ligand L-argininamide.

A personal perspective on chemistry-driven RNA research

In this mini review, we discuss how our understanding of ribonucleic acid (RNA) properties becomes significantly deepened when a broad range of modern chemical and biophysical methods is applied. We span our perspective from RNA solid-phase synthesis and site-specific labeling to single-molecule fluorescence-resonance-energy-transfer imaging and NMR spectroscopy approaches to explore the dynamics of RNA over a broad timescale. We then move on to Fourier-transform-ion-cyclotron-resonance mass spectrometry (FT-ICR-MS) as a powerful technique for RNA sequencing and modification analysis. The novel methodological developments are discussed for selected biological systems that include the thiamine-pyrophosphate riboswitch, HIV and ribosomal A-site RNA, and transfer RNA.

Single transcriptional and translational preQ1 riboswitches adopt similar pre-folded ensembles that follow distinct folding pathways into the same ligand-bound structure

Riboswitches are structural elements in the 5′ untranslated regions of many bacterial messenger RNAs that regulate gene expression in response to changing metabolite concentrations by inhibition of either transcription or translation initiation. The preQ₁ (7-aminomethyl-7-deazaguanine) riboswitch family comprises some of the smallest metabolite sensing RNAs found in nature. Once ligand-bound, the transcriptional Bacillus subtilis and translational Thermoanaerobacter tengcongensis preQ₁ riboswitch aptamers are structurally similar RNA pseudoknots; yet, prior structural studies have characterized their ligand-free conformations as largely unfolded and folded, respectively. In contrast, through single molecule observation, we now show that, at near-physiological Mg²⁺ concentration and pH, both ligand-free aptamers adopt similar pre-folded state ensembles that differ in their ligand-mediated folding. Structure-based Gō-model simulations of the two aptamers suggest that the ligand binds late (Bacillus subtilis) and early (Thermoanaerobacter tengcongensis) relative to pseudoknot folding, leading to the proposal that the principal distinction between the two riboswitches lies in their relative tendencies to fold via mechanisms of conformational selection and induced fit, respectively. These mechanistic insights are put to the test by rationally designing a single nucleotide swap distal from the ligand binding pocket that we find to predictably control the aptamers′ pre-folded states and their ligand binding affinities.

Tuning a riboswitch response through structural extension of a pseudoknot

Structural and dynamic features of RNA folding landscapes represent critical aspects of RNA function in the cell and are particularly central to riboswitch-mediated control of gene expression. Here, using single-molecule fluorescence energy transfer imaging, we explore the folding dynamics of the preQ1 class II riboswitch, an upstream mRNA element that regulates downstream encoded modification enzymes of queuosine biosynthesis. For reasons that are not presently understood, the classical pseudoknot fold of this system harbors an extra stem-loop structure within its 3'-terminal region immediately upstream of the Shine-Dalgarno sequence that contributes to formation of the ligand-bound state. By imaging ligand-dependent preQ1 riboswitch folding from multiple structural perspectives, we reveal that the extra stem-loop strongly influences pseudoknot dynamics in a manner that decreases its propensity to spontaneously fold and increases its responsiveness to ligand binding. We conclude that the extra stem-loop sensitizes this RNA to broaden the dynamic range of the ON/OFF regulatory switch.

Chemical biology 2012: from drug targets to biological systems and back

Multiple sites sharing a common target: This year's EMBO conference on chemical biology encouraged over 340 researchers to come to Heidelberg, Germany, and discuss the use of diverse chemical strategies and tools to investigate biological questions and better understand cellular processes.

2'-SCF3 uridine-a powerful label for probing structure and function of RNA by 19F NMR spectroscopy

Fluorishing: the Togni reagent allows efficient synthetic access to fluorine-labeled RNA molecules. These are in turn highly useful for NMR spectroscopic analyses of secondary and tertiary structures, RNA-protein interactions, and functionality of riboswitch modules.