Folding of the lysine riboswitch: importance of peripheral elements for transcriptional regulation

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.

Multidimensional Applications and Challenges of Riboswitches in Biosensing and Biotherapy

Riboswitches have received significant attention over the last two decades for their multiple functionalities and great potential for applications in various fields. This article highlights and reviews the recent advances in biosensing and biotherapy. These fields involve a wide range of applications, such as food safety detection, environmental monitoring, metabolic engineering, live cell imaging, wearable biosensors, antibacterial drug targets, and gene therapy. The discovery, origin, and optimization of riboswitches are summarized to help readers better understand their multidimensional applications. Finally, this review discusses the multidimensional challenges and development of riboswitches in order to further expand their potential for novel applications.

Riboswitches, from cognition to transformation

Riboswitches are functional RNA elements that regulate gene expression by directly detecting metabolites. Twenty years have passed since it was first discovered, researches on riboswitches are becoming increasingly standardized and refined, which could significantly promote people's cognition of RNA function as well. Here, we focus on some representative orphan riboswitches, enumerate the structural and functional transformation and artificial design of riboswitches including the coupling with ribozymes, hoping to attain a comprehensive understanding of riboswitch research.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

Cotranscriptional RNA strand exchange underlies the gene regulation mechanism in a purine-sensing transcriptional riboswitch

RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

A structural intermediate pre-organizes the add adenine riboswitch for ligand recognition

Riboswitches are RNA sequences that regulate gene expression by undergoing structural changes upon the specific binding of cellular metabolites. Crystal structures of purine-sensing riboswitches have revealed an intricate network of interactions surrounding the ligand in the bound complex. The mechanistic details about how the aptamer folding pathway is involved in the formation of the metabolite binding site have been previously shown to be highly important for the riboswitch regulatory activity. Here, a combination of single-molecule FRET and SHAPE assays have been used to characterize the folding pathway of the adenine riboswitch from Vibrio vulnificus. Experimental evidences suggest a folding process characterized by the presence of a structural intermediate involved in ligand recognition. This intermediate state acts as an open conformation to ensure ligand accessibility to the aptamer and folds into a structure nearly identical to the ligand-bound complex through a series of structural changes. This study demonstrates that the add riboswitch relies on the folding of a structural intermediate that pre-organizes the aptamer global structure and the ligand binding site to allow efficient metabolite sensing and riboswitch genetic regulation.

Requirements for efficient ligand-gated co-transcriptional switching in designed variants of the B. subtilis pbuE adenine-responsive riboswitch in E. coli

Riboswitches, generally located in the 5'-leader of bacterial mRNAs, direct expression via a small molecule-dependent structural switch that informs the transcriptional or translational machinery. While the structure and function of riboswitch effector-binding (aptamer) domains have been intensely studied, only recently have the requirements for efficient linkage between small molecule binding and the structural switch in the cellular and co-transcriptional context begun to be actively explored. To address this aspect of riboswitch function, we have performed a structure-guided mutagenic analysis of the B. subtilis pbuE adenine-responsive riboswitch, one of the simplest riboswitches that employs a strand displacement switching mechanism to regulate transcription. Using a cell-based fluorescent protein reporter assay to assess ligand-dependent regulatory activity in E. coli, these studies revealed previously unrecognized features of the riboswitch. Within the aptamer domain, local and long-range conformational dynamics influenced by sequences within helices have a significant effect upon efficient regulatory switching. Sequence features of the expression platform including the pre-aptamer leader sequence, a toehold helix and an RNA polymerase pause site all serve to promote strong ligand-dependent regulation. By optimizing these features, we were able to improve the performance of the B. subtilis pbuE riboswitch in E. coli from 5.6-fold induction of reporter gene expression by the wild type riboswitch to over 120-fold in the top performing designed variant. Together, these data point to sequence and structural features distributed throughout the riboswitch required to strike a balance between rates of ligand binding, transcription and secondary structural switching via a strand exchange mechanism and yield new insights into the design of artificial riboswitches.

High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Deep sea biology is known to thrive at pressures up to ≈1 kbar, which motivates fundamental biophysical studies of biomolecules under such extreme environments. In this work, the conformational equilibrium of the lysine riboswitch has been systematically investigated by single molecule FRET (smFRET) microscopy at pressures up to 1500 bar. The lysine riboswitch preferentially unfolds with increasing pressure, which signals an increase in free volume (ΔV0 > 0) upon folding of the biopolymer. Indeed, the effective lysine binding constant increases quasi-exponentially with pressure rise, which implies a significant weakening of the riboswitch-ligand interaction in a high-pressure environment. The effects of monovalent/divalent cations and osmolytes on folding are also explored to acquire additional insights into cellular mechanisms for adapting to high pressures. For example, we find that although Mg2+ greatly stabilizes folding of the lysine riboswitch (ΔΔG0 < 0), there is negligible impact on changes in free volume (ΔΔV0 ≈ 0) and thus any pressure induced denaturation effects. Conversely, osmolytes (commonly at high concentrations in deep sea marine species) such as the trimethylamine N-oxide (TMAO) significantly reduce free volumes (ΔΔV0 < 0) and thereby diminish pressure-induced denaturation. We speculate that, besides stabilizing RNA structure, enhanced levels of TMAO in cells might increase the dynamic range for competent riboswitch folding by suppressing the pressure-induced denaturation response. This in turn could offer biological advantage for vertical migration of deep-sea species, with impacts on food searching in a resource limited environment.

Unprecedented tunability of riboswitch structure and regulatory function by sub-millimolar variations in physiological Mg2

Riboswitches are cis-acting regulatory RNA biosensors that rival the efficiency of those found in proteins. At the heart of their regulatory function is the formation of a highly specific aptamer–ligand complex. Understanding how these RNAs recognize the ligand to regulate gene expression at physiological concentrations of Mg²⁺ ions and ligand is critical given their broad impact on bacterial gene expression and their potential as antibiotic targets. In this work, we used single-molecule FRET and biochemical techniques to demonstrate that Mg²⁺ ions act as fine-tuning elements of the amino acid-sensing lysC aptamer's ligand-free structure in the mesophile Bacillus subtilis. Mg²⁺ interactions with the aptamer produce encounter complexes with strikingly different sensitivities to the ligand in different, yet equally accessible, physiological ionic conditions. Our results demonstrate that the aptamer adapts its structure and folding landscape on a Mg²⁺-tunable scale to efficiently respond to changes in intracellular lysine of more than two orders of magnitude. The remarkable tunability of the lysC aptamer by sub-millimolar variations in the physiological concentration of Mg²⁺ ions suggests that some single-aptamer riboswitches have exploited the coupling of cellular levels of ligand and divalent metal ions to tightly control gene expression.

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.

Applicability of a computational design approach for synthetic riboswitches

Riboswitches have gained attention as tools for synthetic biology, since they enable researchers to reprogram cells to sense and respond to exogenous molecules. In vitro evolutionary approaches produced numerous RNA aptamers that bind such small ligands, but their conversion into functional riboswitches remains difficult. We previously developed a computational approach for the design of synthetic theophylline riboswitches based on secondary structure prediction. These riboswitches have been constructed to regulate ligand-dependent transcription termination in Escherichia coli. Here, we test the usability of this design strategy by applying the approach to tetracycline and streptomycin aptamers. The resulting tetracycline riboswitches exhibit robust regulatory properties in vivo. Tandem fusions of these riboswitches with theophylline riboswitches represent logic gates responding to two different input signals. In contrast, the conversion of the streptomycin aptamer into functional riboswitches appears to be difficult. Investigations of the underlying aptamer secondary structure revealed differences between in silico prediction and structure probing. We conclude that only aptamers adopting the minimal free energy (MFE) structure are suitable targets for construction of synthetic riboswitches with design approaches based on equilibrium thermodynamics of RNA structures. Further improvements in the design strategy are required to implement aptamer structures not corresponding to the calculated MFE state.

Transcriptional pausing at the translation start site operates as a critical checkpoint for riboswitch regulation

On the basis of nascent transcript sequencing, it has been postulated but never demonstrated that transcriptional pausing at translation start sites is important for gene regulation. Here we show that the Escherichia coli thiamin pyrophosphate (TPP) thiC riboswitch contains a regulatory pause site in the translation initiation region that acts as a checkpoint for thiC expression. By biochemically probing nascent transcription complexes halted at defined positions, we find a narrow transcriptional window for metabolite binding, in which the downstream boundary is delimited by the checkpoint. We show that transcription complexes at the regulatory pause site favour the formation of a riboswitch intramolecular lock that strongly prevents TPP binding. In contrast, cotranscriptional metabolite binding increases RNA polymerase pausing and induces Rho-dependent transcription termination at the checkpoint. Early transcriptional pausing may provide a general mechanism, whereby transient transcriptional windows directly coordinate the sensing of environmental cues and bacterial mRNA regulation.

Riboswitches are non-coding RNA elements that detect metabolites and control expression by regulating mRNA levels or translation. Here, the authors provide evidence that the E. coli thiC riboswitch has a pause site in the translation initiation region that acts as a checkpoint for thiC expression.

Protein-free ribosomal RNA folds to a near-native state in the presence of Mg2+

The assembled bacterial ribosome contains around 50 proteins and many counterions.

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.

Fluorescence-Based Strategies to Investigate the Structure and Dynamics of Aptamer-Ligand Complexes

In addition to the helical nature of double-stranded DNA and RNA, single-stranded oligonucleotides can arrange themselves into tridimensional structures containing loops, bulges, internal hairpins and many other motifs. This ability has been used for more than two decades to generate oligonucleotide sequences, so-called aptamers, that can recognize certain metabolites with high affinity and specificity. More recently, this library of artificially-generated nucleic acid aptamers has been expanded by the discovery that naturally occurring RNA sequences control bacterial gene expression in response to cellular concentration of a given metabolite. The application of fluorescence methods has been pivotal to characterize in detail the structure and dynamics of these aptamer-ligand complexes in solution. This is mostly due to the intrinsic high sensitivity of fluorescence methods and also to significant improvements in solid-phase synthesis, post-synthetic labeling strategies and optical instrumentation that took place during the last decade. In this work, we provide an overview of the most widely employed fluorescence methods to investigate aptamer structure and function by describing the use of aptamers labeled with a single dye in fluorescence quenching and anisotropy assays. The use of 2-aminopurine as a fluorescent analog of adenine to monitor local changes in structure and fluorescence resonance energy transfer (FRET) to follow long-range conformational changes is also covered in detail. The last part of the review is dedicated to the application of fluorescence techniques based on single-molecule microscopy, a technique that has revolutionized our understanding of nucleic acid structure and dynamics. We finally describe the advantages of monitoring ligand-binding and conformational changes, one molecule at a time, to decipher the complexity of regulatory aptamers and summarize the emerging folding and ligand-binding models arising from the application of these single-molecule FRET microscopy techniques.

Opportunities in the design and application of RNA for gene expression control

The past decade of synthetic biology research has witnessed numerous advances in the development of tools and frameworks for the design and characterization of biological systems. Researchers have focused on the use of RNA for gene expression control due to its versatility in sensing molecular ligands and the relative ease by which RNA can be modeled and designed compared to proteins. We review the recent progress in the field with respect to RNA-based genetic devices that are controlled through small molecule and protein interactions. We discuss new approaches for generating and characterizing these devices and their underlying components. We also highlight immediate challenges, future directions and recent applications of synthetic RNA devices in engineered biological systems.

Biophysical Approaches to Bacterial Gene Regulation by Riboswitches

The last decade has witnessed the discovery of a variety of non-coding RNA sequences that perform a broad range of crucial biological functions. Among these, the ability of certain RNA sequences, so-called riboswitches, has attracted considerable interest. Riboswitches control gene expression in response to the concentration of particular metabolites to which they bind without the need for any protein. These RNA switches not only need to adopt a very specific tridimensional structure to perform their function, but also their sequence has been evolutionary optimized to recognize a particular metabolite with high affinity and selectivity. Thus, riboswitches offer a unique opportunity to get fundamental insights into RNA plasticity and how folding dynamics and ligand recognition mechanisms have been efficiently merged to control gene regulation. Because riboswitch sequences have been mostly found in bacterial organisms controlling the expression of genes associated to the synthesis, degradation or transport of crucial metabolites for bacterial survival, they offer exciting new routes for antibiotic development in an era where bacterial resistance is more than ever challenging conventional drug discovery strategies. Here, we give an overview of the architecture, diversity and regulatory mechanisms employed by riboswitches with particular emphasis on the biophysical methods currently available to characterise their structure and functional dynamics.

A kissing loop is important for btuB riboswitch ligand sensing and regulatory control

RNA-based genetic regulation is exemplified by metabolite-binding riboswitches that modulate gene expression through conformational changes. Crystal structures show that the Escherichia coli btuB riboswitch contains a kissing loop interaction that is in close proximity to the bound ligand. To analyze the role of the kissing loop interaction in the riboswitch regulatory mechanism, we used RNase H cleavage assays to probe the structure of nascent riboswitch transcripts produced by the E. coli RNA polymerase. By monitoring the folding of the aptamer, kissing loop, and riboswitch expression platform, we established the conformation of each structural component in the absence or presence of bound adenosylcobalamin. We found that the kissing loop interaction is not essential for ligand binding. However, we showed that kissing loop formation improves ligand binding efficiency and is required to couple ligand binding to the riboswitch conformational changes involved in regulating gene expression. These results support a mechanism by which the btuB riboswitch modulates the formation of a tertiary structure to perform metabolite sensing and regulate gene expression.

Finding consensus stable local optimal structures for aligned RNA sequences and its application to discovering riboswitch elements

Many non-coding RNAs (ncRNAs) can fold into alternate native structures and perform different biological functions. The computational prediction of an ncRNA's alternate native structures can be conducted by analysing the ncRNA's energy landscape. Previously, we have developed a computational approach, RNASLOpt, to predict alternate native structures for a single RNA. In this paper, in order to improve the accuracy of the prediction, we incorporate structural conservation information among a family of related ncRNA sequences to the prediction. We propose a comparative approach, RNAConSLOpt, to produce all possible consensus SLOpt stack configurations that are conserved on the consensus energy landscape of a family of related ncRNAs. Benchmarking tests show that RNAConSLOpt can reduce the number of candidate structures compared with RNASLOpt, and can predict ncRNAs' alternate native structures accurately. Moreover, an application of the proposed pipeline to bacteria in Bacillus genus has discovered several novel riboswitch candidates.

Structure-guided mutational analysis of gene regulation by the Bacillus subtilis pbuE adenine-responsive riboswitch in a cellular context

Riboswitches are a broadly distributed form of RNA-based gene regulation in Bacteria and, more rarely, Archaea and Eukarya. Most often found in the 5'-leader sequence of bacterial mRNAs, they are generally composed of two functional domains: a receptor (aptamer) domain that binds an effector molecule and a regulatory domain (or expression platform) that instructs the expression machinery. One of the most studied riboswitches is the Bacillus subtilis adenine-responsive pbuE riboswitch, which regulates gene expression at the transcriptional level, up-regulating expression in response to increased intracellular effector concentrations. In this work, we analyzed sequence and structural elements that contribute to efficient ligand-dependent regulatory activity in a co-transcriptional and cellular context. Unexpectedly, we found that the P1 helix, which acts as the antitermination element of the switch in this RNA, supported ligand-dependent activation of a reporter gene over a broad spectrum of lengths from 3 to 10 bp. This same trend was also observed using a minimal in vitro single-turnover transcription assay, revealing that this behavior is intrinsic to the RNA sequence. We also found that the sequences at the distal tip of the terminator not directly involved in alternative secondary structure formation are highly important for efficient regulation. These data strongly support a model in which the switch is highly localized to the P1 helix adjacent to the ligand-binding pocket that likely presents a local kinetic block to invasion of the aptamer by the terminator.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

RNA structural analysis by evolving SHAPE chemistry

RNA is central to the flow of biological information. From transcription to splicing, RNA localization, translation, and decay, RNA is intimately involved in regulating every step of the gene expression program, and is thus essential for health and understanding disease. RNA has the unique ability to base-pair with itself and other nucleic acids to form complex structures. Hence the information content in RNA is not simply its linear sequence of bases, but is also encoded in complex folding of RNA molecules. A general chemical functionality that all RNAs have is a 2'-hydroxyl group in the ribose ring, and the reactivity of the 2'-hydroxyl in RNA is gated by local nucleotide flexibility. In other words, the 2'-hydroxyl is reactive at single-stranded and conformationally flexible positions but is unreactive at nucleotides constrained by base-pairing. Recent efforts have been focused on developing reagents that modify RNA as a function of RNA 2' hydroxyl group reactivity. Such RNA structure probing techniques can be read out by primer extension in experiments termed RNA SHAPE (selective 2'- hydroxyl acylation and primer extension). Herein, we describe the efforts devoted to the design and utilization of SHAPE probes for characterizing RNA structure. We also describe current technological advances that are being applied to utilize SHAPE chemistry with deep sequencing to probe many RNAs in parallel. The merging of chemistry with genomics is sure to open the door to genome-wide exploration of RNA structure and function.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Fluorescence tools to investigate riboswitch structural dynamics

Riboswitches are novel regulatory elements that respond to cellular metabolites to control gene expression. They are constituted of highly conserved domains that have evolved to recognize specific metabolites. Such domains, so-called aptamers, are folded into intricate structures to enable metabolite recognition. Over the years, the development of ensemble and single-molecule fluorescence techniques has allowed to probe most of the mechanistic aspects of aptamer folding and ligand binding. In this review, we summarize the current fluorescence toolkit available to study riboswitch structural dynamics. We fist describe those methods based on fluorescent nucleotide analogues, mostly 2-aminopurine (2AP), to investigate short-range conformational changes, including some key steady-state and time-resolved examples that exemplify the versatility of fluorescent analogues as structural probes. The study of long-range structural changes by Förster resonance energy transfer (FRET) is mostly discussed in the context of single-molecule studies, including some recent developments based on the combination of single-molecule FRET techniques with controlled chemical denaturation methods. This article is part of a Special Issue entitled: Riboswitches.

Comparative sequence and structure analysis reveals the conservation and diversity of nucleotide positions and their associated tertiary interactions in the riboswitches

The tertiary motifs in complex RNA molecules play vital roles to either stabilize the formation of RNA 3D structure or to provide important biological functionality to the molecule. In order to better understand the roles of these tertiary motifs in riboswitches, we examined 11 representative riboswitch PDB structures for potential agreement of both motif occurrences and conservations. A total of 61 unique tertiary interactions were found in the reference structures. In addition to the expected common A-minor motifs and base-triples mainly involved in linking distant regions the riboswitch structures three highly conserved variants of A-minor interactions called G-minors were found in the SAM-I and FMN riboswitches where they appear to be involved in the recognition of the respective ligand’s functional groups. From our structural survey as well as corresponding structure and sequence alignments, the agreement between motif occurrences and conservations are very prominent across the representative riboswitches. Our analysis provide evidence that some of these tertiary interactions are essential components to form the structure where their sequence positions are conserved despite a high degree of diversity in other parts of the respective riboswitches sequences. This is indicative of a vital role for these tertiary interactions in determining the specific biological function of riboswitch.

Modularity of select riboswitch expression platforms enables facile engineering of novel genetic regulatory devices

RNA-based biosensors and regulatory devices have received significant attention for their potential in a broad array of synthetic biology applications. One of the primary difficulties in engineering these molecules is the lack of facile methods to link sensory modules, or aptamers, to readout domains. Such efforts typically require extensive screening or selection of sequences that facilitate interdomain communication. Bacteria have evolved a widespread form of gene regulation known as riboswitches that perform this task with sufficient fidelity to control expression of biosynthetic and transport proteins essential for normal cellular homeostasis. In this work, we demonstrate that select riboswitch readout domains, called expression platforms, are modular in that they can host a variety of natural and synthetic aptamers to create novel chimeric RNAs that regulate transcription both in vitro and in vivo. Importantly, this technique does not require selection of device-specific "communication modules" required to transmit ligand binding to the regulatory domain, enabling rapid engineering of novel functional RNAs.

A new telomerase RNA element that is critical for telomere elongation

The stability of chromosome ends, the telomeres, is dependent on the ribonucleoprotein telomerase. In vitro, telomerase requires at least one RNA molecule and a reverse transcriptase-like protein. However, for telomere homeostasis in vivo, additional proteins are required. Telomerase RNAs of different species vary in size and sequence and only few features common to all telomerases are known. Here we show that stem-loop IVc of the Saccharomyces cerevisiae telomerase RNA contains a structural element that is required for telomerase function in vivo. Indeed, the distal portion of stem-loop IVc stimulates telomerase activity in vitro in a way that is independent of Est1 binding on more proximal portions of this stem-loop. Functional analyses of the RNA in vivo reveal that this distal element we call telomerase-stimulating structure (TeSS) must contain a bulged area in single stranded form and also show that Est1-dependent functions such as telomerase import or recruitment are not affected by TeSS. This study thus uncovers a new structural telomerase RNA element implicated in catalytic activity. Given previous evidence for TeSS elements in ciliate and mammalian RNAs, we speculate that this substructure is a conserved feature that is required for optimal telomerase holoenzyme function.

Abstract folding space analysis based on helices

RNA has many pivotal functions especially in the regulation of gene expression by ncRNAs. Identification of their structure is an important requirement for understanding their function. Structure prediction alone is often insufficient for this task, due to algorithmic problems, parameter inaccuracies, and biological peculiarities. Among the latter, there are base modifications, cotranscriptional folding leading to folding traps, and conformational switching as in the case of riboswitches. All these require more in-depth analysis of the folding space. The major drawback, which all methods have to cope with, is the exponential growth of the folding space. Therefore, methods are often limited in the sequence length they can analyze, or they make use of heuristics, sampling, or abstraction. Our approach adopts the abstraction strategy and remedies some problems of existing methods. We introduce a position-specific abstraction based on helices that we term helix index shapes, or hishapes for short. Utilizing a dynamic programming framework, we have implemented this abstraction in the program RNAHeliCes. Furthermore, we developed two hishape-based methods, one for energy barrier estimation, called HiPath, and one for abstract structure comparison, termed HiTed. We demonstrate the superior performance of HiPath compared to other existing methods and the competitive accuracy of HiTed. RNAHeliCes, together with HiPath and HiTed, are available for download at http://www.cyanolab.de/software/RNAHeliCes.htm.

Single-molecule studies of the lysine riboswitch reveal effector-dependent conformational dynamics of the aptamer domain

The lysine riboswitch is a cis-acting RNA genetic regulatory element found in the leader sequence of bacterial mRNAs coding for proteins related to biosynthesis or transport of lysine. Structural analysis of the lysine-binding aptamer domain of this RNA has revealed that it completely encapsulates the ligand and therefore must undergo a structural opening/closing upon interaction with lysine. In this work, single-molecule fluorescence resonance energy transfer (FRET) methods are used to monitor these ligand-induced structural transitions that are central to lysine riboswitch function. Specifically, a model FRET system has been developed for characterizing the lysine dissociation constant as well as the opening/closing rate constants for the Bacillus subtilis lysC aptamer domain. These techniques permit measurement of the dissociation constant (K(D)) for lysine binding of 1.7(5) mM and opening/closing rate constants of 1.4(3) s(-1) and 0.203(7) s(-1), respectively. These rates predict an apparent dissociation constant for lysine binding (K(D,apparent)) of 0.25(9) mM at near physiological ionic strength, which differs markedly from previous reports.

Basis for ligand discrimination between ON and OFF state riboswitch conformations: the case of the SAM-I riboswitch

Riboswitches are RNA elements that bind to effector ligands and control gene expression. Most consist of two domains. S-Adenosyl Methionine (SAM) binds the aptamer domain of the SAM-I riboswitch and induces conformational changes in the expression domain to form an intrinsic terminator (transcription OFF state). Without SAM the riboswitch forms the transcription ON state, allowing read-through transcription. The mechanistic link between the SAM/aptamer recognition event and subsequent secondary structure rearrangement by the riboswitch is unclear. We probed for those structural features of the Bacillus subtilis yitJ SAM-I riboswitch responsible for discrimination between the ON and OFF states by SAM. We designed SAM-I riboswitch RNA segments forming "hybrid" structures of the ON and OFF states. The choice of segment constrains the formation of a partial P1 helix, characteristic of the OFF state, together with a partial antiterminator (AT) helix, characteristic of the ON state. For most choices of P1 vs. AT helix lengths, SAM binds with micromolar affinity according to equilibrium dialysis. Mutational analysis and in-line probing confirm that the mode of SAM binding by hybrid structures is similar to that of the aptamer. Altogether, binding measurements and in-line probing are consistent with the hypothesis that when SAM is present, stacking interactions with the AT helix stabilize a partially formed P1 helix in the hybrids. Molecular modeling indicates that continuous stacking between the P1 and the AT helices is plausible with SAM bound. Our findings raise the possibility that conformational intermediates may play a role in ligand-induced aptamer folding.

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.

Conformational heterogeneity of the SAM-I riboswitch transcriptional ON state: a chaperone-like role for S-adenosyl methionine

Riboswitches are promising targets for the design of novel antibiotics and engineering of portable genetic regulatory elements. There is evidence that variability in riboswitch properties allows tuning of expression for genes involved in different stages of biosynthetic pathways by mechanisms that are not currently understood. Here, we explore the mechanism for tuning of S-adenosyl methionine (SAM)-I riboswitch folding. Most SAM-I riboswitches function at the transcriptional level by sensing the cognate ligand SAM. SAM-I riboswitches orchestrate the biosynthetic pathways of cysteine, methionine, SAM, and so forth. We use base-pair probability predictions to examine the secondary-structure folding landscape of several SAM-I riboswitch sequences. We predict different folding behaviors for different SAM-I riboswitch sequences. We identify several "decoy" base-pairing interactions involving 5' riboswitch residues that can compete with the formation of a P1 helix, a component of the ligand-bound "transcription OFF" state, in the absence of SAM. We hypothesize that blockage of these interactions through SAM contacts contributes to stabilization of the OFF state in the presence of ligand. We also probe folding patterns for a SAM-I riboswitch RNA using constructs with different 3' truncation points experimentally. Folding was monitored through fluorescence, susceptibility to base-catalyzed cleavage, nuclear magnetic resonance, and indirectly through SAM binding. We identify key decision windows at which SAM can affect the folding pathway towards the OFF state. The presence of decoy conformations and differential sensitivities to SAM at different transcript lengths is crucial for SAM-I riboswitches to modulate gene expression in the context of global cellular metabolism.

Structural architecture of an RNA that competitively inhibits RNase L

Activation of RNase L endonuclease activity is part of the mammalian innate immune response to viral infection. The poliovirus RNA genome contains a sequence in its protein-coding region that can act as a competitive inhibitor of RNase L. Mutation, sequence, and functional analysis of this competitive inhibitor RNA (ciRNA) revealed that its activity depends on specific sequences, showed that a loop-loop hairpin interaction forms in the ciRNA, and suggested the presence of a loop E motif. These features lead to the hypothesis that the ciRNA's function is conferred in part by a specific three-dimensional folded RNA architecture. By using a combination of biophysical, mutational, and functional studies, we have mapped features of the three-dimensional architecture of the ciRNA in its unbound form. We show that the loop-loop interaction forms in the free ciRNA and affects the overall structure, perhaps forming long-range tertiary interactions with the loop E motif. Local tight RNA-RNA backbone packing occurs in parts of the structure, but the fold appears to be less stable than many other tightly packed RNAs. This feature may allow the ciRNA to accommodate the translocation of ribosomes and polymerase across this multifunctional region of the viral RNA but also to function as an RNase L inhibitor.

Riboswitch structure in the ligand-free state

Molecular investigations of riboswitches bound to small-molecule effectors have produced a wealth of information on how these molecules achieve high affinity and specificity for a target ligand. X-ray crystal structures have been determined for the ligand-free state for representatives of the preQ₁-I, S-adenosylmethionine I, lysine, and glycine aptamer classes. These structures in conjunction with complimentary techniques, such as in-line probing, NMR spectroscopy, Förster resonance energy transfer, small-angle scattering, and computational simulations, have demonstrated that riboswitches adopt multiple conformations in the absence of ligand. Despite a number of investigations that support ligand-dependent folding, mounting evidence suggests that free-state riboswitches interact with their effectors in the sub-populations of largely prefolded states as embodied by the principle of conformational selection, which has been documented extensively for protein-mediated ligand interactions. Fundamental riboswitch investigations of the bound and free states have advanced our understanding of RNA folding, ligand recognition, and how these factors culminate in communication between an aptamer and its expression platform. An understanding of these topics is essential to comprehend riboswitch gene regulation at the molecular level, which has already provided a basis to understand the mechanism of action of natural antimicrobials.

Constitutive regulatory activity of an evolutionarily excluded riboswitch variant

The exquisite specificity of the adenine-responsive riboswitch toward its cognate metabolite has been shown to arise from the formation of a Watson-Crick interaction between the adenine ligand and residue U65. A recent crystal structure of a U65C adenine aptamer variant has provided a rationale for the phylogenetic conservation observed at position 39 for purine aptamers. The G39-C65 variant adopts a compact ligand-free structure in which G39 is accommodated by the ligand binding site and is base-paired to the cytosine at position 65. Here, we demonstrate using a combination of biochemical and biophysical techniques that the G39-C65 base pair not only severely impairs ligand binding but also disrupts the functioning of the riboswitch in vivo by constitutively activating gene expression. Folding studies using single-molecule FRET revealed that the G39-C65 variant displays a low level of dynamic heterogeneity, a feature reminiscent of ligand-bound wild-type complexes. A restricted conformational freedom together with an ability to significantly fold in monovalent ions are exclusive to the G39-C65 variant. This work provides a mechanistic framework to rationalize the evolutionary exclusion of certain nucleotide combinations in favor of sequences that preserve ligand binding and gene regulation functionalities.

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.