Chemical basis of glycine riboswitch cooperativity

Translation regulation by a guanidine-II riboswitch is highly tunable in sensitivity, dynamic range, and apparent cooperativity

Riboswitches function as important translational regulators in bacteria. Comprehensive mutational analysis of transcriptional riboswitches has been used to probe the energetic intricacies of interplay between the aptamer and expression platform, but translational riboswitches have been inaccessible to massively parallel techniques. The guanidine-II (gdm-II) riboswitch is an exclusively translational class. We have integrated RelE cleavage with next-generation sequencing to quantify ligand-dependent changes in translation initiation for all single and double mutations of the Pseudomonas aeruginosa gdm-II riboswitch, a total of more than 23,000 variants. This extensive mutational analysis is consistent with the prominent features of the bioinformatic consensus. These data indicate, unexpectedly, that direct sequestration of the Shine-Dalgarno sequence is dispensable for riboswitch function. Additionally, this comprehensive data set reveals important positions not identified in previous computational and crystallographic studies. Mutations in the variable linker region stabilize alternate conformations. The double mutant data reveal the functional importance of the previously modeled P0b helix formed by the 5' and 3' tails that serves as the basis for translational control. Additional mutations to GU wobble base pairs in both P1 and P2 reveal how the apparent cooperativity of the system involves an intricate network of communication between the two binding sites. This comprehensive examination of a translational riboswitch's expression platform illuminates how the riboswitch is precisely tuned and tunable with regard to ligand sensitivity, the amplitude of expression between ON and OFF states, and the cooperativity of ligand binding.

Architectures and complex functions of tandem riboswitches

Riboswitch architectures that involve the binding of a single ligand to a single RNA aptamer domain result in ordinary dose-response curves that require approximately a 100-fold change in ligand concentration to cover nearly the full dynamic range for gene regulation. However, by using multiple riboswitches or aptamer domains in tandem, these ligand-sensing structures can produce additional, complex gene control outcomes. In the current study, we have computationally searched for tandem riboswitch architectures in bacteria to provide a more complete understanding of the diverse biological and biochemical functions of gene control elements that are made exclusively of RNA. Numerous different arrangements of tandem homologous riboswitch architectures are exploited by bacteria to create more ‘digital’ gene control devices, which operate over a narrower ligand concentration range. Also, two heterologous riboswitch aptamers are sometimes employed to create two-input Boolean logic gates with various types of genetic outputs. These findings illustrate the sophisticated genetic decisions that can be made by using molecular sensors and switches based only on RNA.

Bioinformatics and Genomic Analyses of the Suitability of Eight Riboswitches for Antibacterial Drug Targets

Antibiotic resistance (AR) is an acute problem that results in prolonged and debilitating illnesses. AR mortality worldwide is growing and causes a pressing need to research novel mechanisms of action and untested target molecules. This article presents in silico analyses of eight bacterial riboswitches for their suitability for antibacterial drug targets. Most bacterial riboswitches are located in the 5′-untranslated region of messenger RNAs, act as allosteric cis-acting gene control elements, and have not been found in humans before. Sensing metabolites, the riboswitches regulate the synthesis of vital cellular metabolites in various pathogenic bacteria. The analyses performed in this article represent a complete and informative genome-wide bioinformatics analysis of the adequacy of eight riboswitches as antibacterial drug targets in different pathogenic bacteria based on four criteria. Due to the ability of the riboswitch to control biosynthetic pathways and transport proteins of essential metabolites and the presence/absence of alternative biosynthetic pathways, we classified them into four groups based on their suitability for use as antibacterial drug targets guided by our in silico analyses. We concluded that some of them are promising targets for antibacterial drug discovery, such as the PreQ1, MoCo RNA, cyclic-di-GMP I, and cyclic-di-GMP II riboswitches.

Development and characterization of a glycine biosensor system for fine-tuned metabolic regulation in Escherichia coli

Background

In vivo biosensors have a wide range of applications, ranging from the detection of metabolites to the regulation of metabolic networks, providing versatile tools for synthetic biology and metabolic engineering. However, in view of the vast array of metabolite molecules, the existing number and performance of biosensors is far from sufficient, limiting their potential applications in metabolic engineering. Therefore, we developed the synthetic glycine-ON and -OFF riboswitches for metabolic regulation and directed evolution of enzyme in Escherichia coli.

Results

The results showed that a synthetic glycine-OFF riboswitch (glyOFF6) and an increased-detection-range synthetic glycine-ON riboswitch (glyON14) were successfully screened from a library based on the Bacillus subtilis glycine riboswitch using fluorescence-activated cell sorting (FACS) and tetA-based dual genetic selection. The two synthetic glycine riboswitches were successfully used in tunable regulation of lactate synthesis, dynamic regulation of serine synthesis and directed evolution of alanine-glyoxylate aminotransferase in Escherichia coli, respectively. Mutants AGXT22 and AGXT26 of alanine-glyoxylate aminotransferase with an increase of 58% and 73% enzyme activity were obtained by using a high-throughput screening platform based on the synthetic glycine-OFF riboswitch, and successfully used to increase the 5-aminolevulinic acid yield of engineered Escherichia coli.

Conclusions

A synthetic glycine-OFF riboswitch and an increased-detection-range synthetic glycine-ON riboswitch were successfully designed and screened. The developed riboswitches showed broad application in tunable regulation, dynamic regulation and directed evolution of enzyme in E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01779-4.

Small RNA Profiling in Mycobacterium Provides Insights Into Stress Adaptability

Mycobacteria encounter a number of environmental changes during infection and respond using different mechanisms. Small RNA (sRNA) is a post-transcriptionally regulatory system for gene functions and has been investigated in many other bacteria. This study used Mycobacterium tuberculosis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) infection models and sequenced whole bacterial RNAs before and after host cell infection. A comparison of differentially expressed sRNAs using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and target prediction was carried out. Six pathogenically relevant stress conditions, growth rate, and morphology were used to screen and identify sRNAs. From these data, a subset of sRNAs was differentially expressed in multiple infection groups and stress conditions. Many were found associated with lipid metabolism. Among them, ncBCG427 was significantly downregulated when BCG entered into macrophages and was associated with increased biofilm formation. The reduction of virulence possibility depends on regulating lipid metabolism.

Regulation of Glycine Cleavage and Detoxification by a Highly Conserved Glycine Riboswitch in Burkholderia spp

The glycine riboswitch is a known regulatory element that is unique in having two aptamers that are joined by a linker region. In this study, we investigated a glycine riboswitch located in the 5' untranslated region of a glycine cleavage system homolog (gcvTHP) in Burkholderia spp. Structure prediction using the sequence generated a model with a glycine binding pocket composed of base-triple interactions (G62-A64-A86 and G65-U84-C85) that are supported by A/G minor interactions (A17-C60-G88 and G16-C61-G87, respectively) and two ribose-zipper motifs (C11-G12 interacting with A248-A247 and C153-U154 interacting with A79-A78) which had not been previously reported. The capacity of the riboswitch to bind to glycine was experimentally validated by native gel assays and the crucial role of interactions that make up the glycine binding pocket were proven by mutations of A17U and G16C which resulted in conformational differences that may lead to dysfunction. Using glycine supplemented minimal media, we were able to prove that the expression of the gcvTHP genes found downstream of the riboswitch responded to the glycine concentrations introduced thus confirming the role of this highly conserved Burkholderia riboswitch and its associated genes as a putative glycine detoxification system in Burkholderia spp.

A novel riboswitch classification based on imbalanced sequences achieved by machine learning

Riboswitch, a part of regulatory mRNA (50-250nt in length), has two main classes: aptamer and expression platform. One of the main challenges raised during the classification of riboswitch is imbalanced data. That is a circumstance in which the records of a sequences of one group are very small compared to the others. Such circumstances lead classifier to ignore minority group and emphasize on majority ones, which results in a skewed classification. We considered sixteen riboswitch families, to be in accord with recent riboswitch classification work, that contain imbalanced sequences. The sequences were split into training and test set using a newly developed pipeline. From 5460 k-mers (k value 1 to 6) produced, 156 features were calculated based on CfsSubsetEval and BestFirst function found in WEKA 3.8. Statistically tested result was significantly difference between balanced and imbalanced sequences (p < 0.05). Besides, each algorithm also showed a significant difference in sensitivity, specificity, accuracy, and macro F-score when used in both groups (p < 0.05). Several k-mers clustered from heat map were discovered to have biological functions and motifs at the different positions like interior loops, terminal loops and helices. They were validated to have a biological function and some are riboswitch motifs. The analysis has discovered the importance of solving the challenges of majority bias analysis and overfitting. Presented results were generalized evaluation of both balanced and imbalanced models, which implies their ability of classifying, to classify novel riboswitches. The Python source code is available at https://github.com/Seasonsling/riboswitch.

Structure and ligand binding of the ADP-binding domain of the NAD+ riboswitch

The nadA motif is the first known NAD+-dependent riboswitch, comprising two similar tandem bulged stem-loop structures. We have determined the structure of the 5' domain 1 of the riboswitch. It has three coaxial helical segments, separated by an ACANCCCC bulge and by an internal loop, with a tertiary contact between them that includes two C:G base pairs. We have determined the structure with a number of ligands related to NADH, but in each case only the ADP moiety is observed. The adenosine adopts an anti conformation, forms multiple hydrogen bonds across the width of the sugar edge of the penultimate C:G base pair of the helix preceding the bulge, and the observed contacts have been confirmed by mutagenesis and calorimetry. Two divalent metal ions play a key structural role at the narrow neck of the bulge. One makes direct bonding contacts to the diphosphate moiety, locking it into position. Thus the nucleobase, ribose, and phosphate groups of the ADP moiety are all specifically recognized by the RNA. The NAD+ riboswitch is modular. Domain 1 is an ADP binding domain that may be ancient and could potentially be used in combination with other ligand binding motifs such as CoA.

Structural Insights into RNA Dimerization: Motifs, Interfaces and Functions

In comparison with the pervasive use of protein dimers and multimers in all domains of life, functional RNA oligomers have so far rarely been observed in nature. Their diminished occurrence contrasts starkly with the robust intrinsic potential of RNA to multimerize through long-range base-pairing ("kissing") interactions, self-annealing of palindromic or complementary sequences, and stable tertiary contact motifs, such as the GNRA tetraloop-receptors. To explore the general mechanics of RNA dimerization, we performed a meta-analysis of a collection of exemplary RNA homodimer structures consisting of viral genomic elements, ribozymes, riboswitches, etc., encompassing both functional and fortuitous dimers. Globally, we found that domain-swapped dimers and antiparallel, head-to-tail arrangements are predominant architectural themes. Locally, we observed that the same structural motifs, interfaces and forces that enable tertiary RNA folding also drive their higher-order assemblies. These feature prominently long-range kissing loops, pseudoknots, reciprocal base intercalations and A-minor interactions. We postulate that the scarcity of functional RNA multimers and limited diversity in multimerization motifs may reflect evolutionary constraints imposed by host antiviral immune surveillance and stress sensing. A deepening mechanistic understanding of RNA multimerization is expected to facilitate investigations into RNA and RNP assemblies, condensates, and granules and enable their potential therapeutical targeting.

Former orphan riboswitches reveal unexplored areas of bacterial metabolism, signaling, and gene control processes

Comparative sequence analyses have been used to discover numerous classes of structured noncoding RNAs, some of which are riboswitches that specifically recognize small-molecule or elemental ion ligands and influence expression of adjacent downstream genes. Determining the correct identity of the ligand for a riboswitch candidate typically is aided by an understanding of the genes under its regulatory control. Riboswitches whose ligands were straightforward to identify have largely been associated with well-characterized metabolic pathways, such as coenzyme or amino acid biosynthesis. Riboswitch candidates whose ligands resist identification, collectively known as orphan riboswitches, are often associated with genes coding for proteins of unknown function, or genes for various proteins with no established link to one another. The cognate ligands for 16 former orphan riboswitch motifs have been identified to date. The successful pursuit of the ligands for these classes has provided insight into areas of biology that are not yet fully explored, such as ion homeostasis, signaling networks, and other previously underappreciated biochemical or physiological processes. Herein we discuss the strategies and methods used to match ligands with orphan riboswitch classes, and overview the lessons learned to inform and motivate ongoing efforts to identify ligands for the many remaining candidates.

The asymmetry and cooperativity of tandem glycine riboswitch aptamers

Glycine riboswitches utilize both single- and tandem-aptamer architectures. In the tandem system, the relative contribution of each aptamer toward gene regulation is not well understood. To dissect these contributions, the effects of 684 single mutants of a tandem ON switch from Bacillus subtilis were characterized for the wild-type construct and binding site mutations that selectively restrict ligand binding to either the first or second aptamer. Despite the structural symmetry of tandem aptamers, the response to these mutations was frequently asymmetrical. Mutations in the first aptamer often significantly weakened the K_1/2, while several mutations in the second aptamer improved the amplitude. These results demonstrate that this ON switch favors ligand binding to the first aptamer. This is in contrast to the tandem OFF switch variant from Vibrio cholerae, which was previously shown to have preferential binding to its second aptamer. A bioinformatic analysis of tandem glycine riboswitches revealed that the two binding pockets are differentially conserved between ON and OFF switches. Altogether, this indicates that tandem ON switch variants preferentially utilize binding to the first aptamer to promote helical switching, while OFF switch variants favor binding to the second aptamer. The data set also revealed a cooperative glycine response when both binding pockets were maximally stabilized with three GC base pairs. This indicates a cooperative response may sometimes be obfuscated by a difference in the affinities of the two aptamers. This conditional cooperativity provides an additional layer of tunability to tandem glycine riboswitches that adds to their versatility as genetic switches.

Classification of Riboswitch Families Using Block Location-Based Feature Extraction (BLBFE) Method

Purpose: Riboswitches are special non-coding sequences usually located in mRNAs' un-translated regions and regulate gene expression and consequently cellular function. Furthermore, their interaction with antibiotics has been recently implicated. This raises more interest in development of bioinformatics tools for riboswitch studies. Herein, we describe the development and employment of novel block location-based feature extraction (BLBFE) method for classification of riboswitches. Methods: We have already developed and reported a sequential block finding (SBF) algorithm which, without operating alignment methods, identifies family specific sequential blocks for riboswitch families. Herein, we employed this algorithm for 7 riboswitch families including lysine, cobalamin, glycine, SAM-alpha, SAM-IV, cyclic-di-GMP-I and SAH. Then the study was extended toward implementation of BLBFE method for feature extraction. The outcome features were applied in various classifiers including linear discriminant analysis (LDA), probabilistic neural network (PNN), decision tree and k-nearest neighbors (KNN) classifiers for classification of the riboswitch families. The performance of the classifiers was investigated according to performance measures such as correct classification rate (CCR), accuracy, sensitivity, specificity and f-score. Results: As a result, average CCR for classification of riboswitches was 87.87%. Furthermore, application of BLBFE method in 4 classifiers displayed average accuracies of 93.98% to 96.1%, average sensitivities of 76.76% to 83.61%, average specificities of 96.53% to 97.69% and average f-scores of 74.9% to 81.91%. Conclusion: Our results approved that the proposed method of feature extraction; i.e. BLBFE method; can be successfully used for classification and discrimination of the riboswitch families with high CCR, accuracy, sensitivity, specificity and f-score values.

Regulatory context drives conservation of glycine riboswitch aptamers

In comparison to protein coding sequences, the impact of mutation and natural selection on the sequence and function of non-coding (ncRNA) genes is not well understood. Many ncRNA genes are narrowly distributed to only a few organisms, and appear to be rapidly evolving. Compared to protein coding sequences, there are many challenges associated with assessment of ncRNAs that are not well addressed by conventional phylogenetic approaches, including: short sequence length, lack of primary sequence conservation, and the importance of secondary structure for biological function. Riboswitches are structured ncRNAs that directly interact with small molecules to regulate gene expression in bacteria. They typically consist of a ligand-binding domain (aptamer) whose folding changes drive changes in gene expression. The glycine riboswitch is among the most well-studied due to the widespread occurrence of a tandem aptamer arrangement (tandem), wherein two homologous aptamers interact with glycine and each other to regulate gene expression. However, a significant proportion of glycine riboswitches are comprised of single aptamers (singleton). Here we use graph clustering to circumvent the limitations of traditional phylogenetic analysis when studying the relationship between the tandem and singleton glycine aptamers. Graph clustering enables a broader range of pairwise comparison measures to be used to assess aptamer similarity. Using this approach, we show that one aptamer of the tandem glycine riboswitch pair is typically much more highly conserved, and that which aptamer is conserved depends on the regulated gene. Furthermore, our analysis also reveals that singleton aptamers are more similar to either the first or second tandem aptamer, again based on the regulated gene. Taken together, our findings suggest that tandem glycine riboswitches degrade into functional singletons, with the regulated gene(s) dictating which glycine-binding aptamer is conserved.

Characterization and Engineering of a Clostridium Glycine Riboswitch and Its Use To Control a Novel Metabolic Pathway for 5-Aminolevulinic Acid Production in Escherichia coli

A riboswitch, a regulatory RNA that controls gene expression by specifically binding a ligand, is an attractive genetic element for the control of conditional gene expression and metabolic pathways. In this study, we identified a glycine riboswitch located in the 5'-untranslated regions of a glycine:proton symporter gene in Clostridium pasteurianum. The glycine riboswitch is shown to contain two tandem aptamers and to function as an activator of expression of genes fused to its expression platform. Results of singlet aptamer experiments indicated that aptamer-2 has a much higher impact on regulating gene expression than aptamer-1. Further, we successfully obtained synthetic glycine-OFF riboswitches using a dual selection approach, and one of them repressed gene expression up to 10.2-fold with an improved dynamic range. The specific glycine-OFF riboswitch can function as an independent repressor in the presence of glycine, and its repression mechanism is inferred from predicted secondary structure. The selected glycine-OFF riboswitch was used to dynamically control the biosynthesis of 5-aminolevulinic acid (5-ALA) in Escherichia coli with an unnatural 5-ALA synthetic pathway, in which glycine plays a key role. It is demonstrated that the use of a synthetic Clostridium glycine-OFF riboswitch can lead to a significant increase (11%) of 5-ALA in E. coli harboring an unnatural biosynthetic pathway.

Ion-dependent mobility effects of the Fusobacterium nucleatum glycine riboswitch aptamer II via site-directed spin-labeling (SDSL) electron paramagnetic resonance (EPR)

Site-directed spin-labeling (SDSL) with continuous wave electron paramagnetic resonance (cw-EPR) spectroscopy was utilized to probe site-specific changes in backbone dynamics that accompany folding of the isolated 84 nucleotide aptamer II domain of the Fusobacterium nucleatum (FN) glycine riboswitch. Spin-labels were incorporated using splinted ligation strategies. Results show differential dynamics for spin-labels incorporated into the backbone at a base-paired and loop region. Additionally, the addition of a biologically relevant concentration of 5 mM  Mg²⁺, to an RNA solution with 100 mM K⁺, folds and compacts the structure, inferred by a reduction in spin-label mobility. Furthermore, when controlling for ionic strength, Mg²⁺ added to the RNA induces more folding/less flexibility at the two sites than RNA with K⁺ alone. Addition of glycine does not alter the dynamics of this singlet aptamer II, indicating that the full length riboswitch construct may be needed for glycine binding and induced conformational changes. This work adds to our growing understanding of how splinted-ligation SDSL can be utilized to interrogate differential dynamics in large dynamic RNAs, providing insights into how RNA folding and structure is differentially stabilized by monovalent versus divalent cations.

Classification of riboswitch sequences using k-mer frequencies

Riboswitches are non-coding RNAs that regulate gene expression by altering the structural conformation of mRNA transcripts. Their regulation mechanism might be exploited for interesting biomedical applications such as drug targets and biosensors. A major challenge consists in accurately identifying metabolite-binding RNA switches which are structurally complex and diverse. In this regard, we investigated the classification of 16 riboswitch families using supervised learning algorithms trained solely with sequence-based features. We generated a reduced feature set and proposed a visual representation to explore its components. We induced Support Vector Machine, Random Forest, Naive Bayes, J48, and HyperPipes classifiers with our proposed feature set and tested their performance over independent data. Our best multi-class classifier achieved F-measure values of 0.996 and 0.966 in the training and test phases, respectively, outperforming those of a previous approach. When compared against BLAST, our best classifiers yielded competitive results. This work shows that the classifiers trained with our sequence-based feature set accurately discriminate riboswitches.

In Vivo Behavior of the Tandem Glycine Riboswitch in Bacillus subtilis

In many bacterial species, the glycine riboswitch is composed of two homologous ligand-binding domains (aptamers) that each bind glycine and act together to regulate the expression of glycine metabolic and transport genes. While the structure and molecular dynamics of the tandem glycine riboswitch have been the subject of numerous in vitro studies, the in vivo behavior of the riboswitch remains largely uncharacterized. To examine the proposed models of tandem glycine riboswitch function in a biologically relevant context, we characterized the regulatory activity of mutations to the riboswitch structure in Bacillus subtilis using β-galactosidase assays. To assess the impact disruptions to riboswitch function have on cell fitness, we introduced these mutations into the native locus of the tandem glycine riboswitch within the B. subtilis genome. Our results indicate that glycine does not need to bind both aptamers for regulation in vivo and mutations perturbing riboswitch tertiary structure have the most severe effect on riboswitch function and gene expression. We also find that in B. subtilis, the glycine riboswitch-regulated gcvT operon is important for glycine detoxification.IMPORTANCE The glycine riboswitch is a unique cis-acting mRNA element that contains two tandem homologous glycine-binding domains that act on a single expression platform to regulate gene expression in response to glycine. While many in vitro experiments have characterized the tandem architecture of the glycine riboswitch, little work has investigated the behavior of this riboswitch in vivo In this study, we analyzed the proposed models of tandem glycine riboswitch regulation in the context of its native locus within the Bacillus subtilis genome and examined how disruptions to glycine riboswitch function impact organismal fitness. Our work offers new insights into riboswitch function in vivo and reinforces the potential of riboswitches as novel antimicrobial targets.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

Singlet glycine riboswitches bind ligand as well as tandem riboswitches

The glycine riboswitch often occurs in a tandem architecture, with two ligand-binding domains (aptamers) followed by a single expression platform. Based on previous observations, we hypothesized that "singlet" versions of the glycine riboswitch, which contain only one aptamer domain, are able to bind glycine if appropriate structural contacts are maintained. An initial alignment of 17 putative singlet riboswitches indicated that the single consensus aptamer domain is flanked by a conserved peripheral stem-loop structure. These singlets were sorted into two subtypes based on whether the active aptamer domain precedes or follows the peripheral stem-loop, and an example of each subtype of singlet riboswitch was characterized biochemically. The singlets possess glycine-binding affinities comparable to those of previously published tandem examples, and the conserved peripheral domains form A-minor interactions with the single aptamer domain that are necessary for ligand-binding activity. Analysis of sequenced genomes identified a significant number of singlet glycine riboswitches. Based on these observations, we propose an expanded model for glycine riboswitch gene control that includes singlet and tandem architectures.

Engineering and characterization of fluorogenic glycine riboswitches

A set of 12 fluorogenic glycine riboswitches with different thermodynamic and kinetic response properties was engineered. For the design of functional riboswitches, a three-part RNA approach was applied based on the idea of linking a RNA sensor, transmitter and actuator part together. For the RNA sensor and actuator part, we used the tandem glycine aptamer structure from Bacillus subtillis, and fluorogenic aptamer Spinach, respectively. To achieve optimal signal transduction from the sensor to the actuator, a riboswitch library with variable transmitter was screened with a microfluidic large-scale integration chip. This allowed us to establish the complete thermodynamic binding profiles of the riboswitch library. Glycine dissociation constants of the 12 strong fluorescence response riboswitches varied between 99.7 and 570 μM. Furthermore, the kinetic glycine binding (k^on), and dissociation (k^off) rates, and corresponding energy barriers of the 10 strongest fluorescence response riboswitches were determined with the same chip platform. k^on and k^off were in the order of 10⁻³s⁻¹ and 10⁻²s⁻¹, respectively. Conclusively, we demonstrate that systematic screening of synthetic and natural linked RNA parts with microfluidic chip technology is an effective approach to rapidly generate fluorogenic metabolite riboswitches with a broad range of biophysical response properties.

Characterization of a natural triple-tandem c-di-GMP riboswitch and application of the riboswitch-based dual-fluorescence reporter

c-di-GMP riboswitches are structured RNAs located in the 5'-untranslated regions (5'-UTRs) of mRNAs that regulate expression of downstream genes in response to changing concentrations of the second messenger c-di-GMP. We discovered three complete c-di-GMP riboswitches (Bc3, Bc4 and Bc5 RNA) with similar structures, which are arranged in tandem to constitute a triple-tandem (Bc3-5 RNA) riboswitch in the 5'-UTR of the cspABCDE mRNA in Bacillus thuringiensis subsp. chinensis CT-43. Our results showed that this natural triple-tandem riboswitch controlled the expression of the reporter gene more stringently and digitally than the double-tandem or single riboswitch. A sandwich-like dual-fluorescence reporter was further constructed by fusing the Bc3-5 RNA gene between the two fluorescence protein genes amcyan and turborfp. This reporter strain was found to exhibit detectable fluorescence color changes under bright field in response to intracellular c-di-GMP level altered by induced expression of diguanylate cyclase (DGC) PleD. Using this system, two putative membrane-bound DGCs from B. thuringiensis and Xanthomonas oryzae were verified to be functional by replacing pleD with the corresponding DGC genes. This report represented the first native triple-tandem riboswitch that was applied to serve as a riboswitch-based dual-fluorescence reporter for the efficient and convenient verification of putative DGC activity in vivo.

Fab Chaperone-Assisted RNA Crystallography (Fab CARC)

Recent discovery of structured RNAs such as ribozymes and riboswitches shows that there is still much to learn about the structure and function of RNAs. Knowledge learned can be employed in both biochemical research and clinical applications. X-ray crystallography gives unparalleled atomic-level structural detail from which functional inferences can be deduced. However, the difficulty in obtaining high-quality crystals and their phasing information make it a very challenging task. RNA crystallography is particularly arduous due to several factors such as RNA's paucity of surface chemical diversity, lability, repetitive anionic backbone, and flexibility, all of which are counterproductive to crystal packing. Here we describe Fab chaperone assisted RNA crystallography (CARC), a systematic technique to increase RNA crystallography success by facilitating crystal packing as well as expediting phase determination through molecular replacement of conserved Fab domains. Major steps described in this chapter include selection of a synthetic Fab library displayed on M13 phage against a structured RNA crystallization target, ELISA for initial choice of binding Fabs, Fab expression followed by protein A affinity then cation exchange chromatography purification, final choice of Fab by binding specificity and affinity as determined by a dot blot assay, and lastly gel filtration purification of a large quantity of chosen Fabs for crystallization.

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.

An Efficient Minimum Free Energy Structure-Based Search Method for Riboswitch Identification Based on Inverse RNA Folding

Riboswitches are RNA genetic control elements that were originally discovered in bacteria and provide a unique mechanism of gene regulation. They work without the participation of proteins and are believed to represent ancient regulatory systems in the evolutionary timescale. One of the biggest challenges in riboswitch research is to find additional eukaryotic riboswitches since more than 20 riboswitch classes have been found in prokaryotes but only one class has been found in eukaryotes. Moreover, this single known class of eukaryotic riboswitch, namely the TPP riboswitch class, has been found in bacteria, archaea, fungi and plants but not in animals. The few examples of eukaryotic riboswitches were identified using sequence-based bioinformatics search methods such as a combination of BLAST and pattern matching techniques that incorporate base-pairing considerations. None of these approaches perform energy minimization structure predictions. There is a clear motivation to develop new bioinformatics methods, aside of the ongoing advances in covariance models, that will sample the sequence search space more flexibly using structural guidance while retaining the computational efficiency of sequence-based methods. We present a new energy minimization approach that transforms structure-based search into a sequence-based search, thereby enabling the utilization of well established sequence-based search utilities such as BLAST and FASTA. The transformation to sequence space is obtained by using an extended inverse RNA folding problem solver with sequence and structure constraints, available within RNAfbinv. Examples in applying the new method are presented for the purine and preQ1 riboswitches. The method is described in detail along with its findings in prokaryotes. Potential uses in finding novel eukaryotic riboswitches and optimizing pre-designed synthetic riboswitches based on ligand simulations are discussed. The method components are freely available for use.

Cooperativity, allostery and synergism in ligand binding to riboswitches

Recent progress in identification and characterization of novel types of non-coding RNAs has proven that RNAs carry out a variety of cellular functions ranging from scaffolding to gene expression control. In both prokaryotic and eukaryotic cells, several classes of non-coding RNAs control expression of dozens of genes in response to specific cues. One of the most interesting and outstanding questions in the RNA field is whether regulatory RNAs are capable of employing basic biological concepts, such as allostery and cooperativity, previously attributed to the function of proteins. Aside from regulatory RNAs that form complementary base pairing with their nucleic acid targets, several RNA classes modulate gene expression via molecular mechanisms which can be paralleled to protein-mediated regulation. Among these RNAs are riboswitches, metabolite-sensing non-coding regulatory elements that adopt intrinsic three-dimensional structures and specifically bind various small molecule ligands. These characteristics of riboswitches make them well-suited for complex regulatory responses observed in allosteric and cooperative protein systems. Here we present an overview of the biochemical, genetic, and structural studies of riboswitches with a major focus on complex regulatory mechanisms and biological principles utilized by riboswitches for such genetic modulation.

Secondary structural entropy in RNA switch (Riboswitch) identification

BACKGROUND

RNA regulatory elements play a significant role in gene regulation. Riboswitches, a widespread group of regulatory RNAs, are vital components of many bacterial genomes. These regulatory elements generally function by forming a ligand-induced alternative fold that controls access to ribosome binding sites or other regulatory sites in RNA. Riboswitch-mediated mechanisms are ubiquitous across bacterial genomes. A typical class of riboswitch has its own unique structural and biological complexity, making de novo riboswitch identification a formidable task. Traditionally, riboswitches have been identified through comparative genomics based on sequence and structural homology. The limitations of structural-homology-based approaches, coupled with the assumption that there is a great diversity of undiscovered riboswitches, suggests the need for alternative methods for riboswitch identification, possibly based on features intrinsic to their structure. As of yet, no such reliable method has been proposed.

RESULTS

We used structural entropy of riboswitch sequences as a measure of their secondary structural dynamics. Entropy values of a diverse set of riboswitches were compared to that of their mutants, their dinucleotide shuffles, and their reverse complement sequences under different stochastic context-free grammar folding models. Significance of our results was evaluated by comparison to other approaches, such as the base-pairing entropy and energy landscapes dynamics. Classifiers based on structural entropy optimized via sequence and structural features were devised as riboswitch identifiers and tested on Bacillus subtilis, Escherichia coli, and Synechococcus elongatus as an exploration of structural entropy based approaches. The unusually long untranslated region of the cotH in Bacillus subtilis, as well as upstream regions of certain genes, such as the sucC genes were associated with significant structural entropy values in genome-wide examinations.

CONCLUSIONS

Various tests show that there is in fact a relationship between higher structural entropy and the potential for the RNA sequence to have alternative structures, within the limitations of our methodology. This relationship, though modest, is consistent across various tests. Understanding the behavior of structural entropy as a fairly new feature for RNA conformational dynamics, however, may require extensive exploratory investigation both across RNA sequences and folding models.

Ligand binding by the tandem glycine riboswitch depends on aptamer dimerization but not double ligand occupancy

The glycine riboswitch predominantly exists as a tandem structure, with two adjacent, homologous ligand-binding domains (aptamers), followed by a single expression platform. The recent identification of a leader helix, the inclusion of which eliminates cooperativity between the aptamers, has reopened the debate over the purpose of the tandem structure of the glycine riboswitch. An equilibrium dialysis-based assay was combined with binding-site mutations to monitor glycine binding in each ligand-binding site independently to understand the role of each aptamer in glycine binding and riboswitch tertiary interactions. A series of mutations disrupting the dimer interface was used to probe how dimerization impacts ligand binding by the tandem glycine riboswitch. While the wild-type tandem riboswitch binds two glycine equivalents, one for each aptamer, both individual aptamers are capable of binding glycine when the other aptamer is unoccupied. Intriguingly, glycine binding by aptamer-1 is more sensitive to dimerization than glycine binding by aptamer-2 in the context of the tandem riboswitch. However, monomeric aptamer-2 shows dramatically weakened glycine-binding affinity. In addition, dimerization of the two aptamers in trans is dependent on glycine binding in at least one aptamer. We propose a revised model for tandem riboswitch function that is consistent with these results, wherein ligand binding in aptamer-1 is linked to aptamer dimerization and stabilizes the P1 stem of aptamer-2, which controls the expression platform.

DNA-rescuable allosteric inhibition of aptamer II ligand affinity by aptamer I element in the shortened Vibrio cholerae glycine riboswitch

Glycine riboswitches contain two aptamers and turn on the expression of downstream genes in bacteria. Although full-length glycine riboswitches were shown to exhibit no glycine-binding cooperativity, the truncated glycine riboswitches were confirmed to bind two glycine molecules cooperatively. Thorough understanding of the ligand-binding cooperativity may shed light on the molecular basis of the cooperativity and help design novel intricate biosensing genetic circuits for application in synthetic biology. A previously proposed sequential model does not readily provide explanation for published data showing a deleterious mutation in the first aptamer inhibiting the glycine binding of the second one. Using the glycine riboswitch from Vibrio cholerae as a model system, we have identified a region in the first aptamer that modulates the second aptamer function especially in the shortened glycine riboswitch. Importantly, this modulation can be rescued by the addition of a complementary oligodeoxynucleotide, demonstrating the feasibility of developing this system into novel genetic circuits that sense both glycine and a DNA signal.

Specific RNA-binding antibodies with a four-amino-acid code

Numerous large non-coding RNAs are rapidly being discovered, and many of them have been shown to play vital roles in gene expression, gene regulation, and human diseases. Given their often structured nature, specific recognition with an antibody fragment becomes feasible and may help define the structure and function of these non-coding RNAs. As demonstrated for protein antigens, specific antibodies may aid in RNA crystal structure elucidation or the development of diagnostic tools and therapeutic drugs targeting disease-causing RNAs. Recent success and limitation of RNA antibody development has made it imperative to generate an effective antibody library specifically targeting RNA molecules. Adopting the reduced chemical diversity design and further restricting the interface diversity to tyrosines, serines, glycines, and arginines only, we have constructed a RNA-targeting Fab library. Phage display selection and downstream characterization showed that this library yielded high-affinity Fabs for all three RNA targets tested. Using a quantitative specificity assay, we found that these Fabs are highly specific, possibly due to the alternate codon design we used to avoid consecutive arginines in the Fab interface. In addition, the effectiveness of the minimal Fab library may challenge our view of the protein-RNA binding interface and provide a unique solution for future design of RNA-binding proteins.

Site-directed spin-labeling strategies and electron paramagnetic resonance spectroscopy for large riboswitches

Genetic regulation effected by RNA riboswitches is governed by ligand-induced structural reorganization with modulation of RNA conformation and dynamics. Characterization of the conformational states of riboswitches in the presence or absence of salts and ligands is important for understanding how interconversion of riboswitch RNA folding states influences function. The methodology of site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is suitable for such studies, wherein site-specific incorporation of a nitroxide radical spin probe allows for local dynamics and conformational changes to be investigated. This chapter reviews a strategy for SDSL-EPR studies of large riboswitches and uses the full length 232 nucleotide (nt) kink-turn motif-containing Vibrio cholerae (VC) glycine riboswitch as an example. Spin-labeling strategies and the challenges of incorporating spin labels into large riboswitches are reviewed and the approach to overcome these challenges is described. Results are subsequently presented illustrating changes in dynamics within the labeled region of the VC glycine riboswitch as observed using SDSL-EPR.

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.

Comparative genomics of metabolic capacities of regulons controlled by cis-regulatory RNA motifs in bacteria

Background

In silico comparative genomics approaches have been efficiently used for functional prediction and reconstruction of metabolic and regulatory networks. Riboswitches are metabolite-sensing structures often found in bacterial mRNA leaders controlling gene expression on transcriptional or translational levels.

An increasing number of riboswitches and other cis-regulatory RNAs have been recently classified into numerous RNA families in the Rfam database. High conservation of these RNA motifs provides a unique advantage for their genomic identification and comparative analysis.

Results

A comparative genomics approach implemented in the RegPredict tool was used for reconstruction and functional annotation of regulons controlled by RNAs from 43 Rfam families in diverse taxonomic groups of Bacteria. The inferred regulons include ~5200 cis-regulatory RNAs and more than 12000 target genes in 255 microbial genomes. All predicted RNA-regulated genes were classified into specific and overall functional categories. Analysis of taxonomic distribution of these categories allowed us to establish major functional preferences for each analyzed cis-regulatory RNA motif family. Overall, most RNA motif regulons showed predictable functional content in accordance with their experimentally established effector ligands. Our results suggest that some RNA motifs (including thiamin pyrophosphate and cobalamin riboswitches that control the cofactor metabolism) are widespread and likely originated from the last common ancestor of all bacteria. However, many more analyzed RNA motifs are restricted to a narrow taxonomic group of bacteria and likely represent more recent evolutionary innovations.

Conclusions

The reconstructed regulatory networks for major known RNA motifs substantially expand the existing knowledge of transcriptional regulation in bacteria. The inferred regulons can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. The obtained genome-wide collection of reference RNA motif regulons is available in the RegPrecise database (http://regprecise.lbl.gov/).

Glycols modulate terminator stem stability and ligand-dependency of a glycine riboswitch

The Bacillus subtilis glycine riboswitch comprises tandem glycine-binding aptamers and a putative terminator stem followed by the gcvT operon. Gene expression is regulated via the sensing of glycine. However, we found that the riboswitch behaves in a "glycine-independent" manner in the presence of polyethylene glycol (PEG) and ethylene glycol. The effect is related to the formation of a terminator stem within the expression platform under such conditions. The results revealed that increasing PEG stabilized the structure of the terminator stem. By contrast, the addition of ethylene glycol destabilized the terminator stem. PEG and ethylene glycol have opposite effects on transcription as well as on stable terminator stem formation. The glycine-independency of the riboswitch and the effects of such glycols might shed light on the evolution of riboswitches.

Refinement of the Listeria monocytogenes σB regulon through quantitative proteomic analysis

σ(B) is an alternative σ factor that regulates stress response and virulence genes in the foodborne pathogen Listeria monocytogenes. To gain further insight into σ(B)-dependent regulatory mechanisms in L. monocytogenes, we (i) performed quantitative proteomic comparisons between the L. monocytogenes parent strain 10403S and an isogenic ΔsigB mutant and (ii) conducted a meta-analysis of published microarray studies on the 10403S σ(B) regulon. A total of 134 genes were found to be significantly positively regulated by σ(B) at the transcriptomic level with >75 % of these genes preceded by putative σ(B)-dependent promoters; 21 of these 134 genes were also found to be positively regulated by σ(B) through proteomics. In addition, 15 proteins were only found to be positively regulated by σ(B) through proteomics analyses, including Lmo1349, a putative glycine cleavage system protein. The lmo1349 gene is preceded by a 5' UTR that functions as a glycine riboswitch, which suggests regulation of glycine metabolism by σ(B) in L. monocytogenes. Herein, we propose a model where σ(B) upregulates pathways that facilitate biosynthesis and uptake of glycine, which may then activate this riboswitch. Our data also (i) identified a number of σ(B)-dependent proteins that appear to be encoded by genes that are co-regulated by multiple transcriptional regulators, in particular PrfA, and (ii) found σ(B)-dependent genes and proteins to be overrepresented in the 'energy metabolism' role category, highlighting contributions of the σ(B) regulon to L. monocytogenes energy metabolism as well as a role of PrfA and σ(B) interaction in regulating aspects of energy metabolism in L. monocytogenes.

Metabolite recognition principles and molecular mechanisms underlying riboswitch function

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

A decade of riboswitches

Riboswitches were discovered in 2002 in bacteria as RNA-based intracellular sensors of vitamin derivatives. During the last decade, naturally occurring RNA sensor elements have been found to bind a range of small metabolites and ions and to exert regulatory control of transcription, translation, splicing, and RNA stability. Extensive biochemical, structural, and genetic studies have established the basic principles underpinning riboswitch function in all three kingdoms of life with implications for developing antibiotics, designing new molecular sensors, and integrating riboswitches into synthetic circuits.

Modulation of quaternary structure and enhancement of ligand binding by the K-turn of tandem glycine riboswitches

Most known glycine riboswitches have two homologous aptamer domains arranged in tandem and separated by a short linker. The two aptamers associate through reciprocal "quaternary" interactions that have been proposed to result in cooperative glycine binding. Recently, the interaptamer linker was found to form helix P0 with a previously unrecognized segment 5' to the first aptamer domain. P0 was shown to increase glycine affinity, abolish cooperativity, and conform to the K-turn motif consensus. We examine the global thermodynamic and structural role of P0 using isothermal titration calorimetry (ITC) and small-angle X-ray scattering (SAXS), respectively. To evaluate the generality of P0 function, we prepared glycine riboswitch constructs lacking and including P0 from Bacillus subtilis, Fusobacterium nucleatum, and Vibrio cholerae. We find that P0 indeed folds into a K-turn, supports partial pre-folding of all three glycine-free RNAs, and is required for ITC observation of glycine binding under physiologic Mg(2+) concentrations. Except for the unusually small riboswitch from F. nucleatum, the K-turn is needed for maximally compacting the glycine-bound states of the RNAs. Formation of a ribonucleoprotein complex between the B. subtilis or the F. nucleatum RNA constructs and the bacterial K-turn binding protein YbxF promotes additional folding of the free riboswitch, and enhances glycine binding. Consistent with the previously reported loss of cooperativity, P0-containing B. subtilis and V. cholerae tandem aptamers bound no more than one glycine molecule per riboswitch. Our results indicate that the P0 K-turn helps organize the quaternary structure of tandem glycine riboswitches, thereby facilitating ligand binding under physiologic conditions.

The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch

Riboswitch operation involves the complex interplay between the aptamer domain and the expression platform. During transcription, these two domains compete against each other for shared sequence. In this study, we explore the cooperative effects of ligand binding and Magnesium interactions in the SAM-I riboswitch in the context of aptamer collapse and anti-terminator formation. Overall, our studies show the apo-aptamer acts as (i) a pre-organized aptamer competent to bind ligand and undergo structural collapse and (ii) a conformation that is more accessible to anti-terminator formation. We show that both Mg(2+) ions and SAM are required for a collapse transition to occur. We then use competition between the aptamer and expression platform for shared sequence to characterize the stability of the collapsed aptamer. We find that SAM and Mg(2+) interactions in the aptamer are highly cooperative in maintaining switch polarity (i.e. aptamer 'off-state' versus anti-terminator 'on-state'). We further show that the aptamer off-state is preferentially stabilized by Mg(2+) and similar divalent ions. Furthermore, the functional switching assay was used to select for phosphorothioate interference, and identifies potential magnesium chelation sites while characterizing their coordinated role with SAM in aptamer stabilization. In addition, we find that Mg(2+) interactions with the apo-aptamer are required for the full formation of the anti-terminator structure, and that higher concentrations of Mg(2+) (>4 mM) shift the equilibrium toward the anti-terminator on-state even in the presence of SAM.

Challenges and opportunities for small molecule aptamer development

Aptamers are single-stranded oligonucleotides that bind to targets with high affinity and selectivity. Their use as molecular recognition elements has emerged as a viable approach for biosensing, diagnostics, and therapeutics. Despite this potential, relatively few aptamers exist that bind to small molecules. Small molecules are important targets for investigation due to their diverse biological functions as well as their clinical and commercial uses. Novel, effective molecular recognition probes for these compounds are therefore of great interest. This paper will highlight the technical challenges of aptamer development for small molecule targets, as well as the opportunities that exist for their application in biosensing and chemical biology.

Challenges and opportunities for small molecule aptamer development

Aptamers are single-stranded oligonucleotides that bind to targets with high affinity and selectivity. Their use as molecular recognition elements has emerged as a viable approach for biosensing, diagnostics, and therapeutics. Despite this potential, relatively few aptamers exist that bind to small molecules. Small molecules are important targets for investigation due to their diverse biological functions as well as their clinical and commercial uses. Novel, effective molecular recognition probes for these compounds are therefore of great interest. This paper will highlight the technical challenges of aptamer development for small molecule targets, as well as the opportunities that exist for their application in biosensing and chemical biology.

Molecular recognition and function of riboswitches

Regulatory mRNAs elements termed riboswitches respond to elevated concentrations of cellular metabolites by modulating expression of associated genes. Riboswitches attain their high metabolite selectivity by capitalizing on the intrinsic tertiary structures of their sensor domains. Over the years, riboswitch structure and folding have been amongst the most researched topics in the RNA field. Most recently, novel structures of single-ligand and cooperative double-ligand sensors have broadened our knowledge of architectural and molecular recognition principles exploited by riboswitches. The structural information has been complemented by extensive folding studies, which have provided several important clues on the formation of ligand-competent conformations and mechanisms of ligand discrimination. These studies have greatly improved our understanding of molecular events in riboswitch-mediated gene expression control and provided the molecular basis for intervention into riboswitch-controlled genetic circuits.

Identification of small molecule inhibitors of the HIV-1 nucleocapsid-stem-loop 3 RNA complex

Stem-loop 3 RNA (SL3) in ψ-RNA is a highly conserved motif in different strains of HIV-1 and serves as a principle determinant for viral packaging. Viral encapsulation is critical for viral replication, and disruption of the nucleocapsid-ψ-RNA complex interferes with viral replication. We have used SL3 RNA as a target for identification of small molecule inhibitors of the interactions of nucleocapsid protein (NCp7) and ψ-RNA. We report the use of computational and high-throughput screening approaches to identify 16 compounds that bind SL3 RNA with micromolar affinities. Among the identified ligands, two molecules, compounds 7 and 17, bind with higher affinity to SL3 RNA than to double- and single-stranded RNAs. Four of the 16 SL3 RNA ligands inhibit interactions between SL3 RNA and NCp7 with micromolar inhibition constants. In general, the identified SL3 ligands have simple molecular structures and low molecular weights and are, therefore, possible lead compounds for the development of ligands that target the elements of ψ-RNA of HIV-1 with high affinity and specificity.

An energetically beneficial leader-linker interaction abolishes ligand-binding cooperativity in glycine riboswitches

Comprised of two aptamers connected by a short nucleotide linker, the glycine riboswitch was the first example of naturally occurring RNA elements reported to bind small organic molecules cooperatively. Earlier works have shown binding of glycine to the second aptamer allows tertiary interactions to be made between the two aptamers, which facilitates binding of a separate glycine molecule to the first aptamer, leading to glycine-binding cooperativity. Prompted by a distinctive protection pattern in the linker region of a minimal glycine riboswitch construct, we have identified a highly conserved (>90%) leader-linker duplex involving leader nucleotides upstream of the previously reported consensus glycine riboswitch sequences. In >50% of the glycine riboswitches, the leader-linker interaction forms a kink-turn motif. Characterization of three glycine ribsowitches showed that the leader-linker interaction improved the glycine-binding affinities by 4.5- to 86-fold. In-line probing and native gel assays with two aptamers in trans suggested synergistic action between glycine-binding and interaptamer interaction during global folding of the glycine riboswitch. Mutational analysis showed that there appeared to be no ligand-binding cooperativity in the glycine riboswitch when the leader-linker interaction is present, and the previously measured cooperativity is simply an artifact of a truncated construct missing the leader sequence.

Automated RNA structure prediction uncovers a kink-turn linker in double glycine riboswitches

The tertiary structures of functional RNA molecules remain difficult to decipher. A new generation of automated RNA structure prediction methods may help address these challenges but have not yet been experimentally validated. Here we apply four prediction tools to a class of double glycine riboswitches that can bind two ligands cooperatively. A novel method (BPPalign), RMdetect, JAR3D, and Rosetta 3D modeling give consistent predictions for a new stem P0 and a kink-turn motif. These elements structure the linker between the RNAs' double aptamers. Chemical mapping on the Fusobacterium nucleatum riboswitch with N-methylisatoic anhydride, dimethyl sulfate and 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate probing, mutate-and-map studies, and mutation/rescue experiments all provide strong evidence for the structured linker. Under solution conditions that permit rigorous thermodynamic analysis, disrupting this helix-junction-helix structure gives 120- and 6-30-fold poorer dissociation constants for the RNA's two glycine-binding transitions, corresponding to an overall energetic impact of 4.3 ± 0.5 kcal/mol. Prior biochemical and crystallography studies did not include this critical element due to over-truncation of the RNA. We speculate that several further undiscovered elements are likely to exist in the flanking regions of this and other functional RNAs, and automated prediction tools can play a useful role in their detection and dissection.

Challenges and opportunities for small molecule aptamer development

Aptamers are single-stranded oligonucleotides that bind to targets with high affinity and selectivity. Their use as molecular recognition elements has emerged as a viable approach for biosensing, diagnostics, and therapeutics. Despite this potential, relatively few aptamers exist that bind to small molecules. Small molecules are important targets for investigation due to their diverse biological functions as well as their clinical and commercial uses. Novel, effective molecular recognition probes for these compounds are therefore of great interest. This paper will highlight the technical challenges of aptamer development for small molecule targets, as well as the opportunities that exist for their application in biosensing and chemical biology.

A two-dimensional mutate-and-map strategy for non-coding RNA structure

Non-coding RNAs fold into precise base-pairing patterns to carry out critical roles in genetic regulation and protein synthesis, but determining RNA structure remains difficult. Here, we show that coupling systematic mutagenesis with high-throughput chemical mapping enables accurate base-pair inference of domains from ribosomal RNA, ribozymes and riboswitches. For a six-RNA benchmark that has challenged previous chemical/computational methods, this 'mutate-and-map' strategy gives secondary structures that are in agreement with crystallography (helix error rates, 2%), including a blind test on a double-glycine riboswitch. Through modelling of partially ordered states, the method enables the first test of an interdomain helix-swap hypothesis for ligand-binding cooperativity in a glycine riboswitch. Finally, the data report on tertiary contacts within non-coding RNAs, and coupling to the Rosetta/FARFAR algorithm gives nucleotide-resolution three-dimensional models (helix root-mean-squared deviation, 5.7 Å) of an adenine riboswitch. These results establish a promising two-dimensional chemical strategy for inferring the secondary and tertiary structures that underlie non-coding RNA behaviour.

Understanding the errors of SHAPE-directed RNA structure modeling

Single-nucleotide-resolution chemical mapping for structured RNA is being rapidly advanced by new chemistries, faster readouts, and coupling to computational algorithms. Recent tests have shown that selective 2'-hydroxyl acylation by primer extension (SHAPE) can give near-zero error rates (0-2%) in modeling the helices of RNA secondary structure. Here, we benchmark the method using six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from Escherichia coli, the P4-P6 domain of the Tetrahymena group I ribozyme, and ligand-bound domains from riboswitches for adenine, cyclic di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs gave an overall false negative rate (FNR) of 17% and a false discovery rate (FDR) of 21%, with at least one helix prediction error in five of the six cases. Extensive variations of data processing, normalization, and modeling parameters did not significantly mitigate modeling errors. Only one varation, filtering out data collected with deoxyinosine triphosphate during primer extension, gave a modest improvement (FNR = 12%, and FDR = 14%). The residual structure modeling errors are explained by the insufficient information content of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping analysis. Beyond these benchmark cases, bootstrapping suggests a low level of confidence (<50%) in the majority of helices in a previously proposed SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA modeling is not always unambiguous, and helix-by-helix confidence estimates, as described herein, may be critical for interpreting results from this powerful methodology.