Riboswitches: small-molecule recognition by gene regulatory RNAs

rsRNASP: A residue-separation-based statistical potential for RNA 3D structure evaluation

Knowledge-based statistical potentials have been shown to be rather effective in protein 3-dimensional (3D) structure evaluation and prediction. Recently, several statistical potentials have been developed for RNA 3D structure evaluation, while their performances are either still at a low level for the test datasets from structure prediction models or dependent on the "black-box" process through neural networks. In this work, we have developed an all-atom distance-dependent statistical potential based on residue separation for RNA 3D structure evaluation, namely rsRNASP, which is composed of short- and long-ranged potentials distinguished by residue separation. The extensive examinations against available RNA test datasets show that rsRNASP has apparently higher performance than the existing statistical potentials for the realistic test datasets with large RNAs from structure prediction models, including the newly released RNA-Puzzles dataset, and is comparable to the existing top statistical potentials for the test datasets with small RNAs or near-native decoys. In addition, rsRNASP is superior to RNA3DCNN, a recently developed scoring function through 3D convolutional neural networks. rsRNASP and the relevant databases are available to the public.

Large Phenotypic Enhancement of Structured Random RNA Pools

Laboratory evolution of functional RNAs has applications in many areas of chemical and synthetic biology. In vitro selections critically depend on the presence of functional molecules, such as aptamers and ribozymes, in the starting sequence pools. For selection of novel functions the pools are typically transcribed from random-sequence DNA templates, yielding a highly diverse set of RNAs that contain a multitude of folds and biochemical activities. The phenotypic potential, the frequency of functional RNAs, is very low, requiring large complexity of starting pools, surpassing 10¹⁵ different sequences, to identify highly active isolates. Furthermore, the majority of random sequences is not structured and has a high propensity for aggregation; the in vitro selection process thus involves not just enrichment of functional RNAs, but also their purification from aggregation-prone "free-riders". We reasoned that purification of the nonaggregating, monomeric subpopulation of a random-sequence RNA pool will yield pools of folded, functional RNAs. We performed six rounds of selection for monomeric sequences and show that the enriched population is compactly folded. In vitro selections originating from various mixtures of the compact pool and a fully random pool showed that sequences from the compact pool always dominate the population once a biochemical activity is detectable. A head-to-head competition of the two pools starting from a low (5 × 10¹²) sequence diversity revealed that the phenotypic potential of the compact pool is about 1000-times higher than the fully random pool. A selection for folded and monomeric RNA pools thus greatly increases the frequency of functional RNAs from that seen in random-sequence pools, providing a facile experimental approach to isolation of highly active functional RNAs from low-diversity populations.

A Modified In Vitro Transcription Approach to Improve RNA Synthesis and Ribozyme Cleavage Efficiency

RNA elements such as catalytic RNA, riboswitch, microRNA, and long non-coding RNA perform a major role in cellular processes. A complete understanding of cellular processes is impossible without knowing the structure-function relationship of participating RNA molecules that ultimately requires large quantities of pure RNAs. Thus, structural/functional analyses of emerging RNAs necessitate revised protocols for improved RNA quantity and quality. Here we present a modified in vitro transcription protocol to enhance ribozyme cleaving efficiency and RNA yield by working on two variables, i.e., incubation temperature and limiting GTPs. Following an improved RNA synthesis, the target RNA is purified from transcription mixture components through denaturing size-exclusion chromatography. The protocol confirms that cyclic elevated incubation temperatures during transcription and increased concentrations of GTPs improve the production rate of RNA. Our modified in vitro transcription method improves the ribozyme cleaving efficiency and targets RNA yield by four- to fivefold that can benefit almost any RNA-related study from protein-RNA interaction analysis to crystallography.

Tracking RNA with light: selection, structure, and design of fluorescence turn-on RNA aptamers

Fluorescence turn-on aptamers, in vitro evolved RNA molecules that bind conditional fluorophores and activate their fluorescence, have emerged as RNA counterparts of the fluorescent proteins. Turn-on aptamers have been selected to bind diverse fluorophores, and they achieve varying degrees of specificity and affinity. These RNA-fluorophore complexes, many of which exceed the brightness of green fluorescent protein and their variants, can be used as tags for visualizing RNA localization and transport in live cells. Structure determination of several fluorescent RNAs revealed that they have diverse, unrelated overall architectures. As most of these RNAs activate the fluorescence of their ligands by restraining their photoexcited states into a planar conformation, their fluorophore binding sites have in common a planar arrangement of several nucleobases, most commonly a G-quartet. Nonetheless, each turn-on aptamer has developed idiosyncratic structural solutions to achieve specificity and efficient fluorescence turn-on. The combined structural diversity of fluorophores and turn-on RNA aptamers has already produced combinations that cover the visual spectrum. Further molecular evolution and structure-guided engineering is likely to produce fluorescent tags custom-tailored to specific applications.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Towards an understanding of RNA structural modalities: a riboswitch case study

A riboswitch is a type of RNA molecule that regulates important biological functions by changing structure, typically under ligand-binding. We assess the extent that these ligand-bound structural alternatives are present in the Boltzmann sample, a standard RNA secondary structure prediction method, for three riboswitch test cases. We use the cluster analysis tool RNAStructProfiling to characterize the different modalities present among the suboptimal structures sampled. We compare these modalities to the putative base pairing models obtained from independent experiments using NMR or fluorescence spectroscopy. We find, somewhat unexpectedly, that profiling the Boltzmann sample captures evidence of ligand-bound conformations for two of three riboswitches studied. Moreover, this agreement between predicted modalities and experimental models is consistent with the classification of riboswitches into thermodynamic versus kinetic regulatory mechanisms. Our results support cluster analysis of Boltzmann samples by RNAStructProfiling as a possible basis for de novo identification of thermodynamic riboswitches, while highlighting the challenges for kinetic ones.

Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χ_OL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Structural Studies of the 3',3'-cGAMP Riboswitch Induced by Cognate and Noncognate Ligands Using Molecular Dynamics Simulation

Riboswtich RNAs can control gene expression through the structural change induced by the corresponding small-molecule ligands. Molecular dynamics simulations and free energy calculations on the aptamer domain of the 3′,3′-cGAMP riboswitch in the ligand-free, cognate-bound and noncognate-bound states were performed to investigate the structural features of the 3′,3′-cGAMP riboswitch induced by the 3′,3′-cGAMP ligand and the specificity of ligand recognition. The results revealed that the aptamer of the 3′,3′-cGAMP riboswitch in the ligand-free state has a smaller binding pocket and a relatively compact structure versus that in the 3′,3′-cGAMP-bound state. The binding of the 3′,3′-cGAMP molecule to the 3′,3′-cGAMP riboswitch induces the rotation of P1 helix through the allosteric communication from the binding sites pocket containing the J1/2, J1/3 and J2/3 junction to the P1 helix. Simultaneously, these simulations also revealed that the preferential binding of the 3′,3′-cGAMP riboswitch to its cognate ligand, 3′,3′-cGAMP, over its noncognate ligand, c-di-GMP and c-di-AMP. The J1/2 junction in the 3′,3′-cGAMP riboswitch contributing to the specificity of ligand recognition have also been found.

Molecular dynamics simulation on the allosteric analysis of the c-di-GMP class I riboswitch induced by ligand binding

Riboswitches are RNA molecules that regulate gene expression using conformation change, affected by binding of small molecule ligands. 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, we use dynamics simulations on a series of models to characterize the ligand-free and ligand-bound aptamer domain of the c-di-GMP class I (GEMM-I) riboswitch. The results revealed that the ligand-free aptamer has a stable state with a folded P2 and P3 helix, an unfolded P1 helix and open binding pocket. The first Mg ions binding to the aptamer is structurally favorable for the successive c-di-GMP binding. The P1 helix forms when c-di-GMP is successive bound. Three key junctions J1/2, J2/3 and J1/3 in the GEMM-I riboswitch contributing to the formation of P1 helix have been found. The binding of the c-di-GMP ligand to the GEMM-I riboswitch induces the riboswitch's regulation through the direct allosteric communication network in GEMM-I riboswitch from the c-di-GMP binding sites in the J1/2 and J1/3 junctions to the P1 helix, the indirect ones from those in the J2/3 and P2 communicating to P1 helix via the J1/2 and J1/3 media.

Structure of the 30 kDa HIV-1 RNA Dimerization Signal by a Hybrid Cryo-EM, NMR, and Molecular Dynamics Approach

Cryoelectron microscopy (cryo-EM) and nuclear magnetic resonance (NMR) spectroscopy are routinely used to determine structures of macromolecules with molecular weights over 65 and under 25 kDa, respectively. We combined these techniques to study a 30 kDa HIV-1 dimer initiation site RNA ([DIS]2; 47 nt/strand). A 9 Å cryo-EM map clearly shows major groove features of the double helix and a right-handed superhelical twist. Simulated cryo-EM maps generated from time-averaged molecular dynamics trajectories (10 ns) exhibited levels of detail similar to those in the experimental maps, suggesting internal structural flexibility limits the cryo-EM resolution. Simultaneous inclusion of the cryo-EM map and 2H-edited NMR-derived distance restraints during structure refinement generates a structure consistent with both datasets and supporting a flipped-out base within a conserved purine-rich bulge. Our findings demonstrate the power of combining global and local structural information from these techniques for structure determination of modest-sized RNAs.

Ligand recognition and helical stacking formation are intimately linked in the SAM-I riboswitch regulatory mechanism

Riboswitches are noncoding mRNA elements that control gene expression by altering their structure upon metabolite binding. Although riboswitch crystal structures provide detailed information about RNA-ligand interactions, little knowledge has been gathered to understand how riboswitches modulate gene expression. Here, we study the molecular recognition mechanism of the S-adenosylmethionine SAM-I riboswitch by characterizing the formation of a helical stacking interaction involving the ligand-binding process. We show that ligand binding is intimately linked to the formation of the helical stacking, which is dependent on the presence of three conserved purine residues that are flanked by stacked helices. We also find that these residues are important for the formation of a crucial long-range base pair formed upon SAM binding. Together, our results lend strong support to a critical role for helical stacking in the folding pathway and suggest a particularly important function in the formation of the long-range base pair.

Covalent attachment and Pro-Pro endopeptidase (PPEP-1)-mediated release of Clostridium difficile cell surface proteins involved in adhesion

In the past decade, Clostridium difficile has emerged as an important gut pathogen. This anaerobic, Gram-positive bacterium is the main cause of infectious nosocomial diarrhea. Whereas much is known about the mechanism through which the C. difficile toxins cause diarrhea, relatively little is known about the dynamics of adhesion and motility, which is mediated by cell surface proteins. This review will discuss the recent advances in our understanding of the sortase-mediated covalent attachment of cell surface (adhesion) proteins to the peptidoglycan layer of C. difficile and their release through the action of a highly specific secreted metalloprotease (Pro-Pro endopeptidase 1, PPEP-1). Specific emphasis will be on a model in which PPEP-1 and its substrates control the switch from a sessile to motile phenotype in C. difficile, and how this is regulated by the cyclic dinucleotide c-di-GMP (3'-5' cyclic dimeric guanosine monophosphate).

Role of lysine binding residues in the global folding of the lysC riboswitch

Riboswitches regulate gene expression by rearranging their structure upon metabolite binding. The lysine-sensing lysC riboswitch is a rare example of an RNA aptamer organized around a 5-way helical junction in which ligand binding is performed exclusively through nucleotides located at the junction core. We have probed whether the nucleotides involved in ligand binding play any role in the global folding of the riboswitch. As predicted, our findings indicate that ligand-binding residues are critical for the lysine-dependent gene regulation mechanism. We also find that these residues are not important for the establishment of key magnesium-dependent tertiary interactions, suggesting that folding and ligand recognition are uncoupled in this riboswitch for the formation of specific interactions. However, FRET assays show that lysine binding results in an additional conformational change, indicating that lysine binding may also participate in a specific folding transition. Thus, in contrast to helical junctions being primary determinants in ribozymes and rRNA folding, we speculate that the helical junction of the lysine-sensing lysC riboswitch is not employed as structural a scaffold to direct global folding, but rather has a different role in establishing RNA-ligand interactions required for riboswitch regulation. Our work suggests that helical junctions may adopt different functions such as the coordination of global architecture or the formation of specific ligand binding site.

Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.

The dual aptamer approach: rational design of a high-affinity FAD aptamer

A design strategy for high-affinity aptamers of complex biomolecules is presented. We developed an RNA with FAD-binding properties by combining known ATP- and FMN-aptamers. Cooperative binding of FAD was shown by SPR spectroscopy and fluorescence assays. The strategy should be transferable to several other biomolecules.

Riboswitches: From living biosensors to novel targets of antibiotics

Riboswitches are generally located in 5'-UTR region of mRNAs and specifically bind small ligands. Following ligand binding, gene expression is controlled mostly by transcription termination, translation inhibition or mRNA degradation processes. More than 30 classes of known riboswitches have been identified by now. Most riboswitches consist of an aptamer domain and an expression platform. The aptamer domain of each class of riboswitch is a conserved structure and stabilizes specific structures of the expression platforms through binding to specific compounds. In this review, we are highlighting most aspects of riboswitch research including the novel riboswitch discoveries, routine methods for discovering and investigating riboswitches along with newly discovered classes and mechanistic principles of riboswitch-mediated gene expression control. Moreover, we will give an overview about the potential of riboswitches as therapeutic targets for antibiotic design and also their utilization as biosensors for molecular detection.

Quantitative and predictive model of kinetic regulation by E. coli TPP riboswitches

Riboswitches are non-coding elements upstream or downstream of mRNAs that, upon binding of a specific ligand, regulate transcription and/or translation initiation in bacteria, or alternative splicing in plants and fungi. We have studied thiamine pyrophosphate (TPP) riboswitches regulating translation of thiM operon and transcription and translation of thiC operon in E. coli, and that of THIC in the plant A. thaliana. For all, we ascertained an induced-fit mechanism involving initial binding of the TPP followed by a conformational change leading to a higher-affinity complex. The experimental values obtained for all kinetic and thermodynamic parameters of TPP binding imply that the regulation by A. thaliana riboswitch is governed by mass-action law, whereas it is of kinetic nature for the two bacterial riboswitches. Kinetic regulation requires that the RNA polymerase pauses after synthesis of each riboswitch aptamer to leave time for TPP binding, but only when its concentration is sufficient. A quantitative model of regulation highlighted how the pausing time has to be linked to the kinetic rates of initial TPP binding to obtain an ON/OFF switch in the correct concentration range of TPP. We verified the existence of these pauses and the model prediction on their duration. Our analysis also led to quantitative estimates of the respective efficiency of kinetic and thermodynamic regulations, which shows that kinetically regulated riboswitches react more sharply to concentration variation of their ligand than thermodynamically regulated riboswitches. This rationalizes the interest of kinetic regulation and confirms empirical observations that were obtained by numerical simulations.

It's a loop world - single strands in RNA as structural and functional elements

Unpaired regions in RNA molecules - loops - are centrally involved in defining the characteristic three-dimensional (3D) architecture of RNAs and are of high interest in RNA engineering and design. Loops adopt diverse, but specific conformations stabilised by complex tertiary structural interactions that provide structural flexibility to RNA structures that would otherwise not be possible if they only consisted of the rigid A-helical shapes usually formed by canonical base pairing. By participating in sequence-non-local contacts, they furthermore contribute to stabilising the overall fold of RNA molecules. Interactions between RNAs and other nucleic acids, proteins, or small molecules are also generally mediated by RNA loop structures. Therefore, the function of an RNA molecule is generally dependent on its loops. Examples include intermolecular interactions between RNAs as part of the microRNA processing pathways, ribozymatic activity, or riboswitch-ligand interactions. Bioinformatics approaches have been successfully applied to the identification of novel RNA structural motifs including loops, local and global RNA 3D structure prediction, and structural and conformational analysis of RNAs and have contributed to a better understanding of the sequence-structure-function relationships in RNA loops.

Prediction of hydrogen and carbon chemical shifts from RNA using database mining and support vector regression

The Biological Magnetic Resonance Data Bank (BMRB) contains NMR chemical shift depositions for over 200 RNAs and RNA-containing complexes. We have analyzed the (1)H NMR and (13)C chemical shifts reported for non-exchangeable protons of 187 of these RNAs. Software was developed that downloads BMRB datasets and corresponding PDB structure files, and then generates residue-specific attributes based on the calculated secondary structure. Attributes represent properties present in each sequential stretch of five adjacent residues and include variables such as nucleotide type, base-pair presence and type, and tetraloop types. Attributes and (1)H and (13)C NMR chemical shifts of the central nucleotide are then used as input to train a predictive model using support vector regression. These models can then be used to predict shifts for new sequences. The new software tools, available as stand-alone scripts or integrated into the NMR visualization and analysis program NMRViewJ, should facilitate NMR assignment and/or validation of RNA (1)H and (13)C chemical shifts. In addition, our findings enabled the re-calibration a ring-current shift model using published NMR chemical shifts and high-resolution X-ray structural data as guides.

Current perspectives on RNA secondary structure probing

The range of roles played by structured RNAs in biological systems is vast. At the same time as we are learning more about the importance of RNA structure, recent advances in reagents, methods and technology mean that RNA secondary structural probing has become faster and more accurate. As a result, the capabilities of laboratories that already perform this type of structural analysis have increased greatly, and it has also become more widely accessible. The present review summarizes established and recently developed techniques. The information we can derive from secondary structural analysis is assessed, together with the areas in which we are likely to see exciting developments in the near future.

Computational prediction of riboswitches

Riboswitches present a ubiquitous genetic regulatory mechanism for prokaryotes and have been found in HIV1, fungi, plants, and even H. sapiens. We present an overview of approaches to predict riboswitch aptamers and, more generally, RNA conformational switches.

TEMPO-derived spin labels linked to the nucleobases adenine and cytosine for probing local structural perturbations in DNA by EPR spectroscopy

Three 2´-deoxynucleosides containing semi-flexible spin labels, namely (T)A, (U)A and (U)C, were prepared and incorporated into deoxyoligonucleotides using the phosphoramidite method. All three nucleosides contain 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) connected to the exocyclic amino group; (T)A directly and (U)A as well as (U)C through a urea linkage. (T)A and (U)C showed a minor destabilization of a DNA duplex, as registered by a small decrease in the melting temperature, while (U)A destabilized the duplex by more than 10 °C. Circular dichroism (CD) measurements indicated that all three labels were accommodated in B-DNA duplex. The mobility of the spin label (T)A varied with different base-pairing partners in duplex DNA, with the (T)A•T pair being the least mobile. Furthermore, (T)A showed decreased mobility under acidic conditions for the sequences (T)A•C and (T)A•G, to the extent that the EPR spectrum of the latter became nearly superimposable to that of (T)A•T. The reduced mobility of the (T)A•C and (T)A•G mismatches at pH 5 is consistent with the formation of (T)AH(+)•C and (T)AH(+)•G, in which protonation of N1 of A allows the formation of an additional hydrogen bond to N3 of C and N7 of G, respectively, with G in a syn-conformation. The urea-based spin labels (U)A and (U)C were more mobile than (T)A, but still showed a minor variation in their EPR spectra when paired with A, G, C or T in a DNA duplex. (U)A and (U)C had similar mobility order for the different base pairs, with the lowest mobility when paired with C and the highest when paired with T.

Site-specific fluorescence dynamics in an RNA 'thermometer' reveals the role of ribosome binding in its temperature-sensitive switch function

RNA thermometers control the translation of several heat shock and virulence genes by their temperature-sensitive structural transitions. Changes in the structure and dynamics of MiniROSE RNA, which regulates translation in the temperature range of 20–45°C, were studied by site specifically replacing seven adenine residues with the fluorescent analog, 2-aminopurine (2-AP), one at a time. Dynamic fluorescence observables of 2-AP-labeled RNAs were compared in their free versus ribosome-bound states for the first time. Noticeably, position dependence of fluorescence observables, which was prominent at 20°C, was persistent even at 45ºC, suggesting the persistence of structural integrity up to 45ºC. Interestingly, position-dependent dispersion of fluorescence lifetime and quenching constant at 45°C was ablated in ribosome-bound state, when compared to those at 20°C, underscoring loss of structural integrity at 45°C, in ribosome-bound RNA. Significant increase in the value of mean lifetime for 2-AP corresponding to Shine–Dalgarno sequences, when the temperature was raised from 20 to 45°C, to values seen in the presence of urea at 45°C was a strong indicator of melting of the 3D structure of MiniROSE RNA at 45°C, only when it was ribosome bound. Taken all together, we propose a model where we invoke that ribosome binding of the RNA thermometer critically regulates temperature sensing functions in MiniROSE RNA.

The importance of helix P1 stability for structural pre-organization and ligand binding affinity of the adenine riboswitch aptamer domain

We report here an in-depth characterization of the aptamer domain of the transcriptional adenine-sensing riboswitch (pbuE) by NMR and fluorescence spectroscopy. By NMR studies, the structure of two aptamer sequences with different lengths of the helix P1, the central element involved in riboswitch conformational switching, was characterized. Hydrogen-bond interactions could be mapped at nucleotide resolution providing information about secondary and tertiary structure, structure homogeneity and dynamics. Our study reveals that the elongation of helix P1 has pronounced effects not only on the local but on the global structure of the apo aptamer domain. The structural differences induced by stabilizing helix P1 were found to be linked to changes of the ligand binding affinity as revealed from analysis of kinetic and thermodynamic data obtained from stopped-flow fluorescence studies. The results provide new insight into the sequence-dependent fine tuning of the structure and function of purine-sensing 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.

Structural determinants for ligand capture by a class II preQ1 riboswitch

Prequeuosine (preQ1) riboswitches are RNA regulatory elements located in the 5' UTR of genes involved in the biosynthesis and transport of preQ1, a precursor of the modified base queuosine universally found in four tRNAs. The preQ1 class II (preQ1-II) riboswitch regulates preQ1 biosynthesis at the translational level. We present the solution NMR structure and conformational dynamics of the 59 nucleotide Streptococcus pneumoniae preQ1-II riboswitch bound to preQ1. Unlike in the preQ1 class I (preQ1-I) riboswitch, divalent cations are required for high-affinity binding. The solution structure is an unusual H-type pseudoknot featuring a P4 hairpin embedded in loop 3, which forms a three-way junction with the other two stems. (13)C relaxation and residual dipolar coupling experiments revealed interhelical flexibility of P4. We found that the P4 helix and flanking adenine residues play crucial and unexpected roles in controlling pseudoknot formation and, in turn, sequestering the Shine-Dalgarno sequence. Aided by divalent cations, P4 is poised to act as a "screw cap" on preQ1 recognition to block ligand exit and stabilize the binding pocket. Comparison of preQ1-I and preQ1-II riboswitch structures reveals that whereas both form H-type pseudoknots and recognize preQ1 using one A, C, or U nucleotide from each of three loops, these nucleotides interact with preQ1 differently, with preQ1 inserting into different grooves. Our studies show that the preQ1-II riboswitch uses an unusual mechanism to harness exquisite control over queuosine metabolism.

Sequence-dependent folding landscapes of adenine riboswitch aptamers

Expression of a large fraction of genes in bacteria is controlled by riboswitches, which are found in the untranslated region of mRNA. Structurally riboswitches have a conserved aptamer domain to which a metabolite binds, resulting in a conformational change in the downstream expression platform. Prediction of the functions of riboswitches requires a quantitative description of the folding landscape so that the barriers and time scales for the conformational change in the switching region in the aptamer can be estimated. Using a combination of all atom molecular dynamics (MD) and coarse-grained model simulations we studied the response of adenine (A) binding add and pbuE A-riboswitches to mechanical force. The two riboswitches contain a structurally similar three-way junction formed by three paired helices, P1, P2, and P3, but carry out different functions. Using pulling simulations, with structures generated in MD simulations, we show that after P1 rips the dominant unfolding pathway in the add A-riboswitch is the rupture of P2 followed by unraveling of P3. In the pbuE A-riboswitch, after P1 unfolds P3 ruptures ahead of P2. The order of unfolding of the helices, which is in accord with single molecule pulling experiments, is determined by the relative stabilities of the individual helices. Our results show that the stability of isolated helices determines the order of assembly and response to force in these non-coding regions. We use the simulated free energy profile for the pbuE A-riboswitch to estimate the time scale for allosteric switching, which shows that this riboswitch is under kinetic control lending additional support to the conclusion based on single molecule pulling experiments. A consequence of the stability hypothesis is that a single point mutation (U28C) in the P2 helix of the add A-riboswitch, which increases the stability of P2, would make the folding landscapes of the two riboswitches similar. This prediction can be tested in single molecule pulling experiments.

Co-crystal structure of a T-box riboswitch stem I domain in complex with its cognate tRNA

In Gram-positive bacteria, T-box riboswitches regulate the expression of aminoacyl-tRNA synthetases and other proteins in response to fluctuating transfer RNA aminoacylation levels under various nutritional states. T-boxes reside in the 5'-untranslated regions of the messenger RNAs they regulate, and consist of two conserved domains. Stem I contains the specifier trinucleotide that base pairs with the anticodon of cognate tRNA. 3' to stem I is the antiterminator domain, which base pairs with the tRNA acceptor end and evaluates its aminoacylation state. Despite high phylogenetic conservation and widespread occurrence in pathogens, the structural basis of tRNA recognition by this riboswitch remains ill defined. Here we demonstrate that the ~100-nucleotide T-box stem I is necessary and sufficient for specific, high-affinity (dissociation constant (Kd) ~150 nM) tRNA binding, and report the structure of Oceanobacillus iheyensis glyQ stem I in complex with its cognate tRNA at 3.2 Å resolution. Stem I recognizes the overall architecture of tRNA in addition to its anticodon, something accomplished by large ribonucleoproteins such as the ribosome, or proteins such as aminoacyl-tRNA synthetases, but is unprecedented for a compact mRNA domain. The C-shaped stem I cradles the L-shaped tRNA, forming an extended (1,604 Å(2)) intermolecular interface. In addition to the specifier-anticodon interaction, two interdigitated T-loops near the apex of stem I stack on the tRNA elbow in a manner analogous to those of the J11/12-J12/11 motif of RNase P and the L1 stalk of the ribosomal E-site. Because these ribonucleoproteins and T-boxes are unrelated, this strategy to recognize a universal tRNA feature probably evolved convergently. Mutually induced fit of stem I and the tRNA exploiting the intrinsic flexibility of tRNA and its conserved post-transcriptional modifications results in high shape complementarity, which in addition to providing specificity and affinity, globally organizes the T-box to orchestrate tRNA-dependent transcription regulation.

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.

Structure of an A-form RNA duplex obtained by degradation of 6S RNA in a crystallization droplet

In the course of a crystallographic study of a 132 nt variant of Aquifex aeolicus 6S RNA, a crystal structure of an A-form RNA duplex containing 12 base pairs was solved at a resolution of 2.6 Å. In fact, the RNA duplex is part of the 6S RNA and was obtained by accidental but precise degradation of the 6S RNA in a crystallization droplet. 6S RNA degradation was confirmed by microscopic observation of crystals and gel electrophoresis of crystallization droplets. The RNA oligomers obtained form regular A-form duplexes containing three GoU wobble-type base pairs, one of which engages in intermolecular contacts through a ribose-zipper motif at the crystal-packing interface.

Quality control certification of RNA aptamer-based detection

Aptamers are single-stranded DNA or RNA molecules with a defined tertiary structure for molecular recognition. Numerous RNA aptamers with excellent binding affinity and specificity have been reported; they constitute an attractive reservoir of molecular recognition elements for biosensor development. However, RNA is relatively unstable owing to spontaneous hydrolysis and nuclease degradation. Thus, RNA aptamer-based biosensors are prone to producing false-positive signals. Here, we present an RNA aptamer biosensor design strategy that utilises an internal control to distinguish target binding from false-positive signals. The sequence of a chosen RNA aptamer is expanded so that it can form three consecutive short RNA-DNA duplexes with 1) a quencher-labelled DNA strand (Q(1)DNA), 2) a dual-fluorophore-labelled DNA strand (F(1)DNAF(2)) and 3) another quencher-labelled DNA strand (Q(2)DNA). The addition of a target releases Q(2)DNA from the duplex assembly, and produces the expected positive signal from F(2). However, the authenticity of target recognition is validated only if no signal is generated from F(1). We have successfully engineered two fluorescent reporters by using an RNA aptamer that binds thrombin and one that binds theophylline. Both reporters show the expected binding affinity and specificity, and are capable of reporting system malfunction when treated with nucleases and chemical denaturants. This strategy provides a simple and reliable way to ensure high-quality detection when RNA aptamers are employed as molecular-recognition elements.

Metabolite recognition principles and molecular mechanisms underlying riboswitch function

Riboswitches are mRNA elements capable of modulating gene expression in response to specific binding by cellular metabolites. Riboswitches exert their function through the interplay of alternative ligand-free and ligand-bound conformations of the metabolite-sensing domain, which in turn modulate the formation of adjacent gene expression controlling elements. X-ray crystallography and NMR spectroscopy have determined three-dimensional structures of virtually all the major riboswitch classes in the ligand-bound state and, for several riboswitches, in the ligand-free state. The resulting spatial topologies have demonstrated the wide diversity of riboswitch folds and revealed structural principles for specific recognition by cognate metabolites. The available three-dimensional information, supplemented by structure-guided biophysical and biochemical experimentation, has led to an improved understanding of how riboswitches fold, what RNA conformations are required for ligand recognition, and how ligand binding can be transduced into gene expression modulation. These studies have greatly facilitated the dissection of molecular mechanisms underlying riboswitch action and should in turn guide the anticipated development of tools for manipulating gene regulatory circuits.

Database proton NMR chemical shifts for RNA signal assignment and validation

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the (1)H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson-Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 4(3) possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA (1)H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.

Resistance to environmental stress requires the RNA chaperones CspC and CspE

Stress response is essential for adaptation and for survival during environmental changes. A major factor in these responses is RpoS (σS), the master regulator of stationary phase and of the general stress response in Escherichia coli. RpoS is regulated by a complex network at several levels - transcription, translation and proteolysis. Previous studies indicated that rpoS transcripts are stabilized by overexpression of the cold shock proteins CspC and CspE. Here we demonstrate the importance of this transcript stabilization in the regulatory networks governing σS activity. We show that upon entry into stationary-phase rpoS transcripts are stabilized and this stabilization is necessary for the increased activity of σS. The increase in rpoS transcript stability requires at least one of the cold shock proteins CspC and CspE. We also show that the cellular concentration of CspC - but not CspE - increases concurrently with the increase in rpoS transcript stability, probably accounting for this stabilization. These data expand previous data showing that upon heat shock there is a reduction in CspC levels, coupled to a reduced half-life of heat shock gene transcripts. Taken together, it appears that CspC levels modulate transcript stability upon exposure to environmental stress while CspE acts as a 'housekeeping RNA chaperone' under general stress conditions.

Allosteric tertiary interactions preorganize the c-di-GMP riboswitch and accelerate ligand binding

Cyclic diguanylate (c-di-GMP) is a bacterial second messenger important for physiologic adaptation and virulence. Class-I c-di-GMP riboswitches are phylogenetically widespread and thought to mediate pleiotropic genetic responses to the second messenger. Previous studies suggest that the RNA aptamer domain switches from an extended free state to a compact, c-di-GMP-bound conformation in which two helical stacks dock side-by-side. Single molecule fluorescence resonance energy transfer (smFRET) experiments now reveal that the free RNA exists in four distinct populations that differ in dynamics in the extended and docked conformations. In the presence of c-di-GMP and Mg(2+), a stably docked population (>30 min) becomes predominant. smFRET mutant analysis demonstrates that tertiary interactions distal to the c-di-GMP binding site strongly modulate the RNA population structure, even in the absence of c-di-GMP. These allosteric interactions accelerate ligand recognition by preorganizing the RNA, favoring rapid c-di-GMP binding.

Synthesis of (6-(13)C)pyrimidine nucleotides as spin-labels for RNA dynamics

We present a (13)C-based isotope labeling protocol for RNA. Using (6-(13)C)pyrimidine phosphoramidite building blocks, site-specific labels can be incorporated into a target RNA via chemical oligonucleotide solid-phase synthesis. This labeling scheme is particularly useful for studying milli- to microsecond dynamics via NMR spectroscopy, as an isolated spin system is a crucial prerequisite to apply Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion type experiments. We demonstrate the applicability for the characterization and detection of functional dynamics on various time scales by incorporating the (6-(13)C)uridine and -cytidine labels into biologically relevant RNAs. The refolding kinetics of a bistable terminator antiterminator segment involved in the gene regulation process controlled by the preQ(1) riboswitch class I was investigated. Using (13)C CPMG relaxation dispersion NMR spectroscopy, the milli- to microsecond dynamics of the HIV-1 transactivation response element RNA and the Varkud satellite stem loop V motif was addressed.

Analysis of lysine recognition and specificity of the Bacillus subtilis L box riboswitch

The ever-changing environment of a bacterial cell requires sophisticated mechanisms to adjust gene expression in response to changes in nutrient availability. L box riboswitch RNAs regulate gene expression in response to cellular lysine (lys) concentrations in the absence of additional regulatory factors. In Bacillus subtilis, binding of lysine (lys) to the L box RNA causes premature transcription termination in the leader region upstream of the lysC coding sequence. To date, little is known about the specific RNA-lys interactions required for transcription termination. In this study, we characterize features of the B. subtilis lysC leader RNA responsible for lys specificity, and structural elements of the lys molecule required for recognition. The wild-type lysC leader RNA can recognize and discriminate between lys and lys analogs. We identified leader RNA variants with mutations in the lys-binding pocket that exhibit changes in the specificity of ligand recognition. These data demonstrate that lysC leader RNA specificity is the result of recognition of ligand features through a series of distinct interactions between lys and nucleotides that comprise the lys-binding pocket, and provide insight into the molecular mechanisms employed by L box riboswitch RNAs to bind and recognize lys.

Folding of the SAM-I riboswitch: a tale with a twist

Riboswitches are ligand-dependent RNA genetic regulators that control gene expression by altering their structures. The elucidation of riboswitch conformational changes before and after ligand recognition is crucial to understand how riboswitches can achieve high ligand binding affinity and discrimination against cellular analogs. The detailed characterization of riboswitch folding pathways suggest that they may use their intrinsic conformational dynamics to sample a large array of structures, some of which being nearly identical to ligand-bound molecules. Some of these structural conformers can be "captured" upon ligand binding, which is crucial for the outcome of gene regulation. Recent studies about the SAM-I riboswitch have revealed unexpected and previously unknown RNA folding mechanisms. For instance, the observed helical twist of the P1 stem upon ligand binding to the SAM-I aptamer adds a new element in the repertoire of RNA strategies for recognition of small metabolites. From an RNA folding perspective, these findings also strongly indicate that the SAM-I riboswitch could achieve ligand recognition by using an optimized combination of conformational capture and induced-fit approaches, a feature that may be shared by other RNA regulatory sequences.

Use of a coenzyme by the glmS ribozyme-riboswitch suggests primordial expansion of RNA chemistry by small molecules

The glmS ribozyme-riboswitch is the first known example of a naturally occurring catalytic RNA that employs a small molecule as a coenzyme. Binding of glucosamine-6-phosphate (GlcN6P) activates self-cleavage of the bacterial ribozyme, which is part of the mRNA encoding the metabolic enzyme GlcN6P-synthetase. Cleavage leads to negative feedback regulation. GlcN6P binds in the active site of the ribozyme, where its amine could function as a general acid and electrostatic catalyst. The ribozyme is pre-folded but inactive in the absence of GlcN6P, demonstrating it has evolved strict dependence on the exogenous small molecule. The ribozyme showcases the ability of RNA to co-opt non-covalently bound small molecules to expand its chemical repertoire. Analogue studies demonstrate that some molecules other than GlcN6P, such as l-serine (but not d-serine), can function as weak activators. This suggests how coenzyme use by RNA world ribozymes may have led to evolution of proteins. Primordial cofactor-dependent ribozymes may have evolved to bind their cofactors covalently. If amino acids were used as cofactors, this could have driven the evolution of RNA aminoacylation. The ability to make covalently bound peptide coenzymes may have further increased the fitness of such primordial ribozymes, providing a selective pressure for the invention of translation.

Prospects for riboswitch discovery and analysis

An expanding number of metabolite-binding riboswitch classes are being discovered in the noncoding portions of bacterial genomes. Findings over the last decade indicate that bacteria commonly use these RNA genetic elements as regulators of metabolic pathways and as mediators of changes in cell physiology. Some riboswitches are surprisingly complex, and they rival protein factors in their structural and functional sophistication. Each new riboswitch discovery expands our knowledge of the biochemical capabilities of RNA, and some give rise to new questions that require additional study to be addressed. Some of the greatest prospects for riboswitch research and some of the more interesting mysteries are discussed in this review.

Influence of ground-state structure and Mg2+ binding on folding kinetics of the guanine-sensing riboswitch aptamer domain

Riboswitch RNAs fold into complex tertiary structures upon binding to their cognate ligand. Ligand recognition is accomplished by key residues in the binding pocket. In addition, it often crucially depends on the stability of peripheral structural elements. The ligand-bound complex of the guanine-sensing riboswitch from Bacillus subtilis, for example, is stabilized by extensive interactions between apical loop regions of the aptamer domain. Previously, we have shown that destabilization of this tertiary loop–loop interaction abrogates ligand binding of the G37A/C61U-mutant aptamer domain (Gsw^loop) in the absence of Mg²⁺. However, if Mg²⁺ is available, ligand-binding capability is restored by a population shift of the ground-state RNA ensemble toward RNA conformations with pre-formed loop–loop interactions. Here, we characterize the striking influence of long-range tertiary structure on RNA folding kinetics and on ligand-bound complex structure, both by X-ray crystallography and time-resolved NMR. The X-ray structure of the ligand-bound complex reveals that the global architecture is almost identical to the wild-type aptamer domain. The population of ligand-binding competent conformations in the ground-state ensemble of Gsw^loop is tunable through variation of the Mg²⁺ concentration. We quantitatively describe the influence of distinct Mg²⁺ concentrations on ligand-induced folding trajectories both by equilibrium and time-resolved NMR spectroscopy at single-residue resolution.

Long-range pseudoknot interactions dictate the regulatory response in the tetrahydrofolate riboswitch

Tetrahydrofolate (THF), a biologically active form of the vitamin folate (B(9)), is an essential cofactor in one-carbon transfer reactions. In bacteria, expression of folate-related genes is controlled by feedback modulation in response to specific binding of THF and related compounds to a riboswitch. Here, we present the X-ray structures of the THF-sensing domain from the Eubacterium siraeum riboswitch in the ligand-bound and unbound states. The structure reveals an "inverted" three-way junctional architecture, most unusual for riboswitches, with the junction located far from the regulatory helix P1 and not directly participating in helix P1 formation. Instead, the three-way junction, stabilized by binding to the ligand, aligns the riboswitch stems for long-range tertiary pseudoknot interactions that contribute to the organization of helix P1 and therefore stipulate the regulatory response of the riboswitch. The pterin moiety of the ligand docks in a semiopen pocket adjacent to the junction, where it forms specific hydrogen bonds with two moderately conserved pyrimidines. The aminobenzoate moiety stacks on a guanine base, whereas the glutamate moiety does not appear to make strong interactions with the RNA. In contrast to other riboswitches, these findings demonstrate that the THF riboswitch uses a limited number of available determinants for ligand recognition. Given that modern antibiotics target folate metabolism, the THF riboswitch structure provides insights on mechanistic aspects of riboswitch function and may help in manipulating THF levels in pathogenic bacteria.

Riboswitches: structures and mechanisms

A critical feature of the hypothesized RNA world would have been the ability to control chemical processes in response to environmental cues. Riboswitches present themselves as viable candidates for a sophisticated mechanism of regulatory control in RNA-based life. These regulatory elements in the modern world are most commonly found in the 5'-untranslated regions of bacterial mRNAs, directly interacting with metabolites as a means of regulating expression of the coding region via a secondary structural switch. In this review, we focus on recent insights into how these RNAs fold into complex architectures capable of both recognizing a specific small molecule compound and exerting regulatory control over downstream sequences, with an emphasis on transcriptional regulation.

Molecular insights into the ligand-controlled organization of the SAM-I riboswitch

S-adenosylmethionine (SAM) riboswitches are widespread in bacteria, and up to five different SAM riboswitch families have been reported, highlighting the relevance of SAM regulation. On the basis of crystallographic and biochemical data, it has been postulated, but never demonstrated, that ligand recognition by SAM riboswitches involves key conformational changes in the RNA architecture. We show here that the aptamer follows a two-step hierarchical folding selectively induced by metal ions and ligand binding, each of them leading to the formation of one of the two helical stacks observed in the crystal structure. Moreover, we find that the anti-antiterminator P1 stem is rotated along its helical axis upon ligand binding, a mechanistic feature that could be common to other riboswitches. We also show that the nonconserved P4 helical domain is used as an auxiliary element to enhance the ligand-binding affinity. This work provides the first comprehensive characterization, to our knowledge, of a ligand-controlled riboswitch folding pathway.