Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy

Linker-Mediated Inactivation of the SAM-II Domain in the Tandem SAM-II/SAM-V Riboswitch

Tandem SAM-II/SAM-V riboswitch belongs to a class of riboswitches found in the marine bacterium 'Candidatus Pelagibacter ubique'. Previous studies have demonstrated that these riboswitches have the potential for digital modulation of gene expression at both the transcriptional and translational levels. In this study, we investigate the conformational changes in the tandem SAM-II/SAM-V riboswitch binding to S-adenosylmethionine (SAM) using selective 2'-hydroxyl acylation analyzed by the primer extension (SHAPE) assay, small-angle X-ray scattering (SAXS), and oligos depressing probing. Our findings reveal that the linker between SAM-II/SAM-V aptamers blocks the SAM response of the SAM-II domain. This result proposes a new mechanism for gene expression regulation, where the ligand-binding functions of tandem riboswitches can be selectively masked or released through a linker.

A sensitive and scalable fluorescence anisotropy single stranded RNA targeting approach for monitoring riboswitch conformational states

The capacity of riboswitches to undergo conformational changes in response to binding their native ligands is closely tied to their functional roles and is an attractive target for antimicrobial drug design. Here, we established a probe-based fluorescence anisotropy assay to monitor riboswitch conformational switching with high sensitivity and throughput. Using the Bacillus subtillis yitJ S-Box (SAM-I), Fusobacterium nucleatum impX RFN element of (FMN) and class-I cyclic-di-GMP from Vibrio cholerae riboswitches as model systems, we developed short fluorescent DNA probes that specifically recognize either ligand-free or -bound riboswitch conformational states. We showed that increasing concentrations of native ligands cause measurable and reproducible changes in fluorescence anisotropy that correlate with riboswitch conformational changes observed by native gel analysis. Furthermore, we applied our assay to several ligand analogues and confirmed that it can discriminate between ligands that bind, triggering the native conformational change, from those that bind without causing the conformational change. This new platform opens the possibility of high-throughput screening compound libraries to identify potential new antibiotics that specifically target functional conformational changes in riboswitches.

An optimized purification protocol for enzymatically synthesized S-adenosyl-L-methionine (SAM) for applications in solution state infrared spectroscopic studies

S-adenosyl-L-methionine (SAM) is an abundant biomolecule used by methyltransferases to regulate a wide range of essential cellular processes such as gene expression, cell signaling, protein functions, and metabolism. Despite considerable effort, there remain many specificity challenges associated with designing small molecule inhibitors for methyltransferases, most of which exhibit off-target effects. Interestingly, NMR evidence suggests that SAM undergoes conformeric exchange between several states when free in solution. Infrared spectroscopy can detect different conformers of molecules if present in appreciable populations. When SAM is noncovalently bound within enzyme active sites, the nature and the number of different conformations of the molecule are likely to be altered from when it is free in solution. If there are unique structures or different numbers of conformers between different methyltransferase active sites, solution-state information may provide promising structural leads to increase inhibitor specificity for a particular methyltransferase. Toward this goal, frequencies measured in SAM's infrared spectra must be assigned to the motions of specific atoms via isotope incorporation at discrete positions. The incorporation of isotopes into SAM's structure can be accomplished via an established enzymatic synthesis using isotopically labeled precursors. However, published protocols produced an intense and highly variable IR signal which overlapped with many of the signals from SAM rendering comparison between isotopes challenging. We observed this intense absorption to be from co-purifying salts and the SAM counterion, producing a strong, broad signal at 1100 cm-1. Here, we report a revised SAM purification protocol that mitigates the contaminating salts and present the first IR spectra of isotopically labeled CD3-SAM. These results provide a foundation for isotopic labeling experiments of SAM that will define which atoms participate in individual molecular vibrations, as a means to detect specific molecular conformations.

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.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

The Second Class of Tetrahydrofolate (THF-II) Riboswitches Recognizes the Tetrahydrofolic Acid Ligand via Local Conformation Changes

Riboswitches are regulatory noncoding RNAs found in bacteria, fungi and plants, that modulate gene expressions through structural changes in response to ligand binding. Understanding how ligands interact with riboswitches in solution can shed light on the molecular mechanisms of this ancient regulators. Previous studies showed that riboswitches undergo global conformation changes in response to ligand binding to relay information. Here, we report conformation switching models of the recently discovered tetrahydrofolic acid-responsive second class of tetrahydrofolate (THF-II) riboswitches in response to ligand binding. Using a combination of selective 2′-hydroxyl acylation, analyzed by primer extension (SHAPE) assay, 3D modeling and small-angle X-ray scattering (SAXS), we found that the ligand specifically recognizes and reshapes the THF-II riboswitch loop regions, but does not affect the stability of the P3 helix. Our results show that the THF-II riboswitch undergoes only local conformation changes in response to ligand binding, rearranging the Loop1-P3-Loop2 region and rotating Loop1 from a ~120° angle to a ~75° angle. This distinct conformation changes suggest a unique regulatory mechanism of the THF-II riboswitch, previously unseen in other riboswitches. Our findings may contribute to the fields of RNA sensors and drug design.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS)

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine-II riboswitch binds the small cation, guanidinium, and forms a kissing loop-loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg2+ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small-angle X-ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg2+ and ligand binding.

Incoherent dual regulation by a SAM-II riboswitch controlling translation at a distance

In Sinorhizobium meliloti, the methionine biosynthesis genes metA and metZ are preceded by S-adenosyl-L-methionine (SAM) riboswitches of the SAM-II class. Upon SAM binding, structural changes in the metZ riboswitch were predicted to cause transcriptional termination, generating the sRNA RZ. By contrast, the metA riboswitch was predicted to regulate translation from an AUG1 codon. However, downstream of the metA riboswitch, we found a putative Rho-independent terminator and an in-frame AUG2 codon, which may contribute to metA regulation. We validated the terminator between AUG1 and AUG2, which generates the sRNA RA1 that is processed to RA2. Under high SAM conditions, the activities of the metA and metZ promoters and the steady-state levels of the read-through metA and metZ mRNAs were decreased, while the levels of the RZ and RA2 sRNAs were increased. Under these conditions, the sRNAs and the mRNAs were stabilized. Reporter fusion experiments revealed that the Shine–Dalgarno (SD) sequence in the metA riboswitch is required for translation, which, however, starts 74 nucleotides downstream at AUG2, suggesting a novel translation initiation mechanism. Further, the reporter fusion data supported the following model of RNA-based regulation: Upon SAM binding by the riboswitch, the SD sequence is sequestered to downregulate metA translation, while the mRNA is stabilized. Thus, the SAM-II riboswitches fulfil incoherent, dual regulation, which probably serves to ensure basal metA and metZ mRNA levels under high SAM conditions. This probably helps to adapt to changing conditions and maintain SAM homoeostasis.

An integrative NMR-SAXS approach for structural determination of large RNAs defines the substrate-free state of a trans-cleaving Neurospora Varkud Satellite ribozyme

The divide-and-conquer strategy is commonly used for protein structure determination, but its applications to high-resolution structure determination of RNAs have been limited. Here, we introduce an integrative approach based on the divide-and-conquer strategy that was undertaken to determine the solution structure of an RNA model system, the Neurospora VS ribozyme. NMR and SAXS studies were conducted on a minimal trans VS ribozyme as well as several isolated subdomains. A multi-step procedure was used for structure determination that first involved pairing refined NMR structures with SAXS data to obtain structural subensembles of the various subdomains. These subdomain structures were then assembled to build a large set of structural models of the ribozyme, which was subsequently filtered using SAXS data. The resulting NMR-SAXS structural ensemble shares several similarities with the reported crystal structures of the VS ribozyme. However, a local structural difference is observed that affects the global fold by shifting the relative orientation of the two three-way junctions. Thus, this finding highlights a global conformational change associated with substrate binding in the VS ribozyme that is likely critical for its enzymatic activity. Structural studies of other large RNAs should benefit from similar integrative approaches that allow conformational sampling of assembled fragments.

Attenuator LRR - a regulatory tool for modulating gene expression in Gram-positive bacteria

With the rapid development of synthetic biology in recent years, particular attention has been paid to RNA devices, especially riboswitches, because of their significant and diverse regulatory roles in prokaryotic and eukaryotic cells. Due to the limited performance and context-dependence of riboswitches, only a few of them (such as theophylline, tetracycline and ciprofloxacin riboswitches) have been utilized as regulatory tools in biotechnology. In the present study, we demonstrated that a ribosome-dependent ribo-regulator, LRR, discovered in our previous work, exhibits an attractive regulatory performance. Specifically, it offers a 60-fold change in expression in the presence of retapamulin and a low level of leaky expression of about 1-2% without antibiotics. Moreover, LRR can be combined with different promoters and performs well in Bacillus thuringiensis, B. cereus, B. amyloliquefaciens, and B. subtilis. Additionally, LRR also functions in the Gram-negative bacterium Escherichia coli. Furthermore, we demonstrate its ability to control melanin metabolism in B. thuringiensis BMB171. Our results show that LRR can be applied to regulate gene expression, construct genetic circuits and tune metabolic pathways, and has great potential for many applications in synthetic biology.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Strategies for Covalent Labeling of Long RNAs

The introduction of chemical modifications into long RNA molecules at specific positions for visualization, biophysical investigations, diagnostic and therapeutic applications still remains challenging. In this review, we present recent approaches for covalent internal labeling of long RNAs. Topics included are the assembly of large modified RNAs via enzymatic ligation of short synthetic oligonucleotides and synthetic biology approaches preparing site-specifically modified RNAs via in vitro transcription using an expanded genetic alphabet. Moreover, recent approaches to employ deoxyribozymes (DNAzymes) and ribozymes for RNA labeling and RNA methyltransferase based labeling strategies are presented. We discuss the potentials and limits of the individual methods, their applicability for RNAs with several hundred to thousands of nucleotides in length and indicate future directions in the field.

Structures and target RNA preferences of the RNA-binding protein family of IGF2BPs: An overview

Insulin-like growth factor 2 mRNA-binding proteins (IMPs, IGF2BPs) act in mRNA transport and translational control but are oncofetal tumor marker proteins. The IMP protein family represents a number of bona fide multi-domain RNA-binding proteins with up to six RNA-binding domains, resulting in a high complexity of possible modes of interactions with target mRNAs. Their exact mechanism in stability control of oncogenic mRNAs is only partially understood. Our and other laboratories' recent work has significantly pushed the understanding of IMP protein specificities both toward RNA engagement and between each other from NMR and crystal structures serving the basis for systematic biochemical and functional investigations. We here summarize the known structural and biochemical information about IMP RNA-binding domains and their RNA preferences. The article also touches on the respective roles of RNA secondary and protein tertiary structures for specific RNA-protein complexes, including the limited knowledge about IMPs' protein-protein interactions, which are often RNA mediated.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Magnesium controls aptamer-expression platform switching in the SAM-I riboswitch

Investigations of most riboswitches remain confined to the ligand-binding aptamer domain. However, during the riboswitch mediated transcription regulation process, the aptamer domain and the expression platform compete for a shared strand. If the expression platform dominates, an anti-terminator helix is formed, and the transcription process is active (ON state). When the aptamer dominates, transcription is terminated (OFF state). Here, we use an expression platform switching experimental assay and structure-based electrostatic simulations to investigate this ON-OFF transition of the full length SAM-I riboswitch and its magnesium concentration dependence. Interestingly, we find the ratio of the OFF population to the ON population to vary non-monotonically as magnesium concentration increases. Upon addition of magnesium, the aptamer domain pre-organizes, populating the OFF state, but only up to an intermediate magnesium concentration level. Higher magnesium concentration preferentially stabilizes the anti-terminator helix, populating the ON state, relatively destabilizing the OFF state. Magnesium mediated aptamer-expression platform domain closure explains this relative destabilization of the OFF state at higher magnesium concentration. Our study reveals the functional potential of magnesium in controlling transcription of its downstream genes and underscores the importance of a narrow concentration regime near the physiological magnesium concentration ranges, striking a balance between the OFF and ON states in bacterial gene regulation.

Piperine Modulates Protein Mediated Toxicity in Fragile X-Associated Tremor/Ataxia Syndrome through Interacting Expanded CGG Repeat (r(CGG)exp) RNA

An expansion of CGG tandem repeats in the 5' untranslated region (5'-UTR) of fragile X mental retardation 1 (FMR1) gene causes fragile X-associated tremor/ataxia syndrome (FXTAS). The transcripts of these expanded repeats r(CGG)^exp either form RNA foci or undergo the repeat-associated non-ATG (RAN) translation that produces toxic homopolymeric proteins in neuronal cells. The discovery of small molecule modulators that possess a strong binding affinity and high selectivity to these toxic expanded repeats RNA could be a promising therapeutic approach to cure the expanded repeat-associated neurological diseases. Therefore, here we sought to test the therapeutic potential of a natural alkaloid, piperine, by assessing its ability to bind and neutralize the toxicity of r(CGG)^exp RNA motif. To accomplish this first, we have determined the affinity of piperine to r(CGG)^exp RNA using fluorescence-based binding assay and isothermal titration calorimetry assay. These assays showed that piperine forms a thermodynamically favorable interaction with r(CGG)^exp RNA with high selectivity to the G-rich RNA motif. Interaction of piperine with r(CGG)^exp motif was further validated using several biophysical techniques such as CD, CD melting, NMR spectroscopy, and gel retardation assay. Moreover, piperine was also found to be effective for improving the r(CGG)^exp associated splicing defects and RAN translation in a FXTAS cell model system. Our results effectively provided the evidence that piperine strongly interacts with r(CGG)^exp RNA and could be used as a suitable candidate for therapeutic development against FXTAS.

A Mn-sensing riboswitch activates expression of a Mn2+/Ca2+ ATPase transporter in Streptococcus

Maintaining manganese (Mn) homeostasis is important for the virulence of numerous bacteria. In the human respiratory pathogen Streptococcus pneumoniae, the Mn-specific importer PsaBCA, exporter MntE, and transcriptional regulator PsaR establish Mn homeostasis. In other bacteria, Mn homeostasis is controlled by yybP-ykoY family riboswitches. Here, we characterize a yybP-ykoY family riboswitch upstream of the mgtA gene encoding a PII-type ATPase in S. pneumoniae, suggested previously to function in Ca2+ efflux. We show that the mgtA riboswitch aptamer domain adopts a canonical yybP-ykoY structure containing a three-way junction that is compacted in the presence of Ca2+ or Mn2+ at a physiological Mg2+ concentration. Although Ca2+ binds to the RNA aptamer with higher affinity than Mn2+, in vitro activation of transcription read-through of mgtA by Mn2+ is much greater than by Ca2+. Consistent with this result, mgtA mRNA and protein levels increase ≈5-fold during cellular Mn stress, but only in genetic backgrounds of S. pneumoniae and Bacillus subtilis that exhibit Mn2+ sensitivity, revealing that this riboswitch functions as a failsafe 'on' signal to prevent Mn2+ toxicity in the presence of high cellular Mn2+. In addition, our results suggest that the S. pneumoniae yybP-ykoY riboswitch functions to regulate Ca2+ efflux under these conditions.

Systematic Evaluation of Genetic and Environmental Factors Affecting Performance of Translational Riboswitches

Since their discovery, riboswitches have been attractive tools for the user-controlled regulation of gene expression in bacterial systems. Riboswitches facilitate small molecule mediated fine-tuning of protein expression, making these tools of great use to the synthetic biology community. However, the use of riboswitches is often restricted due to context dependent performance and limited dynamic range. Here, we report the drastic improvement of a previously developed orthogonal riboswitch achieved through in vivo functional selection and optimization of flanking coding and noncoding sequences. The behavior of the derived riboswitches was mapped under a wide array of growth and induction conditions, using a structured Design of Experiments approach. This approach successfully improved the maximal protein expression levels 8.2-fold relative to the original riboswitches, and the dynamic range was improved to afford riboswitch dependent control of 80-fold. The optimized orthogonal riboswitch was then integrated downstream of four endogenous stress promoters, responsive to phosphate starvation, hyperosmotic stress, redox stress, and carbon starvation. These responsive stress promoter-riboswitch devices were demonstrated to allow for tuning of protein expression up to ∼650-fold in response to both environmental and cellular stress responses and riboswitch dependent attenuation. We envisage that these riboswitch stress responsive devices will be useful tools for the construction of advanced genetic circuits, bioprocessing, and protein expression.

Structure and ligand binding of the SAM-V riboswitch

SAM-V is one of the class of riboswitches that bind S-adenosylmethione, regulating gene expression by controlling translation. We have solved the crystal structure of the metY SAM-V riboswitch bound to its SAM ligand at 2.5 Å resolution. The RNA folds as an H-type pseudoknot, with a major-groove triple helix in which resides the SAM ligand binding site. The bound SAM adopts an elongated conformation aligned with the axis of the triple helix, and is held at either end by hydrogen bonding to the adenine and the amino acid moieties. The central sulfonium cation makes electrostatic interactions with an U:A.U base triple, so conferring specificity. We propose a model in which SAM binding leads to association of the triplex third strand that stabilizes a short helix and occludes the ribosome binding site. Thus the new structure explains both ligand specificity and the mechanism of genetic control.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Cooperation between Magnesium and Metabolite Controls Collapse of the SAM-I Riboswitch

The S-adenosylmethionine (SAM)-I riboswitch is a noncoding RNA that regulates the transcription termination process in response to metabolite (SAM) binding. The aptamer portion of the riboswitch may adopt an open or closed state depending on the presence of metabolite. Although the transition between the open and closed states is critical for the switching process, its atomistic details are not well understood. Using atomistic simulations, we calculate the effect of SAM and magnesium ions on the folding free energy landscape of the SAM-I riboswitch. These molecular simulation results are consistent with our previous wetlab experiments and aid in interpreting the SHAPE probing measurements. Here, molecular dynamics simulations explicitly identify target RNA motifs sensitive to magnesium ions and SAM. In the simulations, we observe that, whereas the metabolite mostly stabilizes the P1 and P3 helices, magnesium serves an important role in stabilizing a pseudoknot interaction between the P2 and P4 helices, even at high metabolite concentrations. The pseudoknot stabilization by magnesium, in combination with P1 stabilization by SAM, explains the requirement of both SAM and magnesium to form the fully collapsed metabolite-bound closed state of the SAM-I riboswitch. In the absence of SAM, frequent open-to-closed conformational transitions of the pseudoknot occur, akin to breathing. These pseudoknot fluctuations disrupt the binding site by facilitating fluctuations in the 5'-end of helix P1. Magnesium biases the landscape toward a collapsed state (preorganization) by coordinating pseudoknot and 5'-P1 fluctuations. The cooperation between SAM and magnesium in stabilizing important tertiary interactions elucidates their functional significance in transcription regulation.

An excited state underlies gene regulation of a transcriptional riboswitch

Riboswitches control gene expression through ligand-dependent structural rearrangements of the sensing aptamer domain. However, we found that the Bacillus cereus fluoride riboswitch aptamer adopts identical tertiary structures in solution with and without ligand. Using chemical-exchange saturation transfer (CEST) NMR spectroscopy, we revealed that the structured ligand-free aptamer transiently accesses a low-populated (∼1%) and short-lived (∼3 ms) excited conformational state that unravels a conserved 'linchpin' base pair to signal transcription termination. Upon fluoride binding, this highly localized, fleeting process is allosterically suppressed, which activates transcription. We demonstrated that this mechanism confers effective fluoride-dependent gene activation over a wide range of transcription rates, which is essential for robust toxicity responses across diverse cellular conditions. These results unveil a novel switching mechanism that employs ligand-dependent suppression of an aptamer excited state to coordinate regulatory conformational transitions rather than adopting distinct aptamer ground-state tertiary architectures, exemplifying a new mode of ligand-dependent RNA regulation.

Long-Range Interactions in Riboswitch Control of Gene Expression

Riboswitches are widespread RNA motifs that regulate gene expression in response to fluctuating metabolite concentrations. Known primarily from bacteria, riboswitches couple specific ligand binding and changes in RNA structure to mRNA expression in cis. Crystal structures of the ligand binding domains of most of the phylogenetically widespread classes of riboswitches, each specific to a particular metabolite or ion, are now available. Thus, the bound states-one end point-have been thoroughly characterized, but the unbound states have been more elusive. Consequently, it is less clear how the unbound, sensing riboswitch refolds into the ligand binding-induced output state. The ligand recognition mechanisms of riboswitches are diverse, but we find that they share a common structural strategy in positioning their binding sites at the point of the RNA three-dimensional fold where the residues farthest from one another in sequence meet. We review how riboswitch folds adhere to this fundamental strategy and propose future research directions for understanding and harnessing their ability to specifically control gene expression.

A magnesium-induced triplex pre-organizes the SAM-II riboswitch

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

In vivo, large-scale preparation of uniformly (15)N- and site-specifically (13)C-labeled homogeneous, recombinant RNA for NMR studies

Knowledge of how ribonucleic acid (RNA) structures fold to form intricate, three-dimensional structures has provided fundamental insights into understanding the biological functions of RNA. Nuclear magnetic resonance (NMR) spectroscopy is a particularly useful high-resolution technique to investigate the dynamic structure of RNA. Effective study of RNA by NMR requires enrichment with isotopes of (13)C or (15)N or both. Here, we present a method to produce milligram quantities of uniformly (15)N- and site-specifically (13)C-labeled RNAs using wild-type K12 and mutant tktA Escherichia coli in combination with a tRNA-scaffold approach. The method includes a double selection protocol to obtain an E. coli clone with consistently high expression of the recombinant tRNA-scaffold. We also present protocols for the purification of the tRNA-scaffold from a total cellular RNA extract and the excision of the RNA of interest from the tRNA-scaffold using DNAzymes. Finally, we showcase NMR applications to demonstrate the benefit of using in vivo site-specifically (13)C-labeled RNA.

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

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

SAXS studies of RNA: structures, dynamics, and interactions with partners

Small-angle X-ray scattering, SAXS, is a powerful and easily employed experimental technique that provides solution structures of macromolecules. The size and shape parameters derived from SAXS provide global structural information about these molecules in solution and essentially complement data acquired by other biophysical methods. As applied to protein systems, SAXS is a relatively mature technology: sophisticated tools exist to acquire and analyze data, and to create structural models that include dynamically flexible ensembles. Given the expanding appreciation of RNA's biological roles, there is a need to develop comparable tools to characterize solution structures of RNA, including its interactions with important biological partners. We review the progress toward achieving this goal, focusing on experimental and computational innovations. The use of multiphase modeling, absolute calibration and contrast variation methods, among others, provides new and often unique ways of visualizing this important biological molecule and its essential partners: ions, other RNAs, or proteins. WIREs RNA 2016, 7:512-526. doi: 10.1002/wrna.1349 For further resources related to this article, please visit the WIREs website.

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.

SAM-II Riboswitch Samples at least Two Conformations in Solution in the Absence of Ligand: Implications for Recognition

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

Chemo-enzymatic synthesis of site-specific isotopically labeled nucleotides for use in NMR resonance assignment, dynamics and structural characterizations

Stable isotope labeling is central to NMR studies of nucleic acids. Development of methods that incorporate labels at specific atomic positions within each nucleotide promises to expand the size range of RNAs that can be studied by NMR. Using recombinantly expressed enzymes and chemically synthesized ribose and nucleobase, we have developed an inexpensive, rapid chemo-enzymatic method to label ATP and GTP site specifically and in high yields of up to 90%. We incorporated these nucleotides into RNAs with sizes ranging from 27 to 59 nucleotides using in vitro transcription: A-Site (27 nt), the iron responsive elements (29 nt), a fluoride riboswitch from Bacillus anthracis (48 nt), and a frame-shifting element from a human corona virus (59 nt). Finally, we showcase the improvement in spectral quality arising from reduced crowding and narrowed linewidths, and accurate analysis of NMR relaxation dispersion (CPMG) and TROSY-based CEST experiments to measure μs-ms time scale motions, and an improved NOESY strategy for resonance assignment. Applications of this selective labeling technology promises to reduce difficulties associated with chemical shift overlap and rapid signal decay that have made it challenging to study the structure and dynamics of large RNAs beyond the 50 nt median size found in the PDB.

Structural and Dynamic Basis for Low-Affinity, High-Selectivity Binding of L-Glutamine by the Glutamine Riboswitch

Naturally occurring L-glutamine riboswitches occur in cyanobacteria and marine metagenomes, where they reside upstream of genes involved in nitrogen metabolism. By combining X-ray, NMR, and MD, we characterized an L-glutamine-dependent conformational transition in the Synechococcus elongatus glutamine riboswitch from tuning fork to L-shaped alignment of stem segments. This transition generates an open ligand-binding pocket with L-glutamine selectivity enforced by Mg(2+)-mediated intermolecular interactions. The transition also stabilizes the P1 helix through a long-range "linchpin" Watson-Crick G-C pair-capping interaction, while melting a short helix below P1 potentially capable of modulating downstream readout. NMR data establish that the ligand-free glutamine riboswitch in Mg(2+) solution exists in a slow equilibrium between flexible tuning fork and a minor conformation, similar, but not identical, to the L-shaped bound conformation. We propose that an open ligand-binding pocket combined with a high conformational penalty for forming the ligand-bound state provide mechanisms for reducing binding affinity while retaining high selectivity.

Stable isotope-labeled RNA phosphoramidites to facilitate dynamics by NMR

Given that Ribonucleic acids (RNAs) are a central hub of various cellular processes, methods to synthesize these RNAs for biophysical studies are much needed. Here, we showcase the applicability of 6-(13)C-pyrimidine phosphoramidites to introduce isolated (13)C-(1)H spin pairs into RNAs up to 40 nucleotides long. The method allows the incorporation of 6-(13)C-uridine and -cytidine residues at any desired position within a target RNA. By site-specific positioning of the (13)C-label using RNA solid phase synthesis, these stable isotope-labeling patterns are especially well suited to resolve resonance assignment ambiguities. Of even greater importance, the labeling pattern affords accurate quantification of important functional transitions of biologically relevant RNAs (e.g., riboswitch aptamer domains, viral RNAs, or ribozymes) in the μs- to ms time regime and beyond without complications of one bond carbon scalar couplings. We outline the chemical synthesis of the 6-(13)C-pyrimidine building blocks and their use in RNA solid phase synthesis and demonstrate their utility in Carr Purcell Meiboom Gill relaxation dispersion, ZZ exchange, and chemical exchange saturation transfer NMR experiments.

Methods for analysis of size-exclusion chromatography-small-angle X-ray scattering and reconstruction of protein scattering

Size-exclusion chromatography in line with small-angle X-ray scattering (SEC-SAXS) has emerged as an important method for investigation of heterogeneous and self-associating systems, but presents specific challenges for data processing including buffer subtraction and analysis of overlapping peaks. This paper presents novel methods based on singular value decomposition (SVD) and Guinier-optimized linear combination (LC) to facilitate analysis of SEC-SAXS data sets and high-quality reconstruction of protein scattering directly from peak regions. It is shown that Guinier-optimized buffer subtraction can reduce common subtraction artifacts and that Guinier-optimized linear combination of significant SVD basis components improves signal-to-noise and allows reconstruction of protein scattering, even in the absence of matching buffer regions. In test cases with conventional SAXS data sets for cytochrome c and SEC-SAXS data sets for the small GTPase Arf6 and the Arf GTPase exchange factors Grp1 and cytohesin-1, SVD-LC consistently provided higher quality reconstruction of protein scattering than either direct or Guinier-optimized buffer subtraction. These methods have been implemented in the context of a Python-extensible Mac OS X application known as Data Evaluation and Likelihood Analysis (DELA), which provides convenient tools for data-set selection, beam intensity normalization, SVD, and other relevant processing and analytical procedures, as well as automated Python scripts for common SAXS analyses and Guinier-optimized reconstruction of protein scattering.

Small-angle X-ray scattering: a bridge between RNA secondary structures and three-dimensional topological structures

Whereas the structures of small to medium-sized well folded RNA molecules often can be determined by either X-ray crystallography or NMR spectroscopy, obtaining structural information for large RNAs using experimental, computational, or combined approaches remains a major interest and challenge. RNA is very sensitive to small-angle X-ray scattering (SAXS) due to high electron density along phosphate-sugar backbones, whose scattering contribution dominates SAXS intensity. For this reason, SAXS is particularly useful in obtaining global RNA structural information that outlines backbone topologies and, therefore, molecular envelopes. Such information is extremely valuable in bridging the gap between the secondary structures and three-dimensional topological structures of RNA molecules, particularly those that have proven difficult to study using other structure-determination methods. Here we review published results of RNA topological structures derived from SAXS data or in combination with other experimental data, as well as details on RNA sample preparation for SAXS experiments.

Characterizing excited conformational states of RNA by NMR spectroscopy

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

Multiple conformational states of riboswitches fine-tune gene regulation

Riboswitches are structured regions of mRNAs that modulate gene expression in response to specific binding of low molecular-weight ligands. They function by induced transitions between different functional conformations. The standard model assumed that the two functional states, the ligand-bound and ligand-free state, populated only two stable conformations. Recent discoveries of multiple conformations for the apo-state and holo-state of riboswitches challenge this model. Moreover, it becomes evident that detected conformational heterogeneity--mostly in the apo-state--provides sensitivity to multiple environmental inputs for riboswitch-based gene-regulation.

Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT

Integration and calibration of molecular dynamics simulations with experimental data remain a challenging endeavor. We have developed a novel method to integrate chemical probing experiments with molecular simulations of RNA molecules by using a native structure-based model. Selective 2'-hydroxyl acylation by primer extension (SHAPE) characterizes the mobility of each residue in the RNA. Our method, SHAPE-FIT, automatically optimizes the potential parameters of the force field according to measured reactivities from SHAPE. The optimized parameter set allows simulations of dynamics highly consistent with SHAPE probing experiments. Such atomistic simulations, thoroughly grounded in experiment, can open a new window on RNA structure-function relations.

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

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

Structure and function of pseudoknots involved in gene expression control

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

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.

Common themes and differences in SAM recognition among SAM riboswitches

The recent discovery of short cis-acting RNA elements termed riboswitches has caused a paradigm shift in our understanding of genetic regulatory mechanisms. The three distinct superfamilies of S-adenosyl-l-methionine (SAM) riboswitches are the most commonly found riboswitch classes in nature. These RNAs represent three independent evolutionary solutions to achieve specific SAM recognition. This review summarizes research on 1) modes of gene regulatory mechanisms, 2) common themes and differences in ligand recognition, and 3) ligand-induced conformational dynamics among SAM riboswitch families. The body of work on the SAM riboswitch families constitutes a useful primer to the topic of gene regulatory RNAs as a whole. This article is part of a Special Issue entitled: Riboswitches.

Site-specific labeling of RNA at internal ribose hydroxyl groups: terbium-assisted deoxyribozymes at work

A general and efficient single-step method was established for site-specific post-transcriptional labeling of RNA. Using Tb(3+) as accelerating cofactor for deoxyribozymes, various labeled guanosines were site-specifically attached to 2'-OH groups of internal adenosines in in vitro transcribed RNA. The DNA-catalyzed 2',5'-phosphodiester bond formation proceeded efficiently with fluorescent, spin-labeled, biotinylated, or cross-linker-modified guanosine triphosphates. The sequence context of the labeling site was systematically analyzed by mutating the nucleotides flanking the targeted adenosine. Labeling of adenosines in a purine-rich environment showed the fastest reactions and highest yields. Overall, practically useful yields >70% were obtained for 13 out of 16 possible nucleotide (nt) combinations. Using this approach, we demonstrate preparative labeling under mild conditions for up to ~160-nt-long RNAs, including spliceosomal U6 small nuclear RNA and a cyclic-di-AMP binding riboswitch RNA.