Identifying small-molecules binding sites in RNA conformational ensembles with SHAMAN.

The rational targeting of RNA with small molecules is hampered by our still limited understanding of RNA structural and dynamic properties. Most in silico tools for binding site identification rely on static structures and therefore cannot face the challenges posed by the dynamic nature of RNA molecules. Here, we present SHAMAN, a computational technique to identify potential small-molecule binding sites in RNA structural ensembles. SHAMAN enables exploring the conformational landscape of RNA with atomistic molecular dynamics simulations and at the same time identifying RNA pockets in an efficient way with the aid of probes and enhanced-sampling techniques. In our benchmark composed of large, structured riboswitches as well as small, flexible viral RNAs, SHAMAN successfully identifies all the experimentally resolved pockets and ranks them among the most favorite probe hotspots. Overall, SHAMAN sets a solid foundation for future drug design efforts targeting RNA with small molecules, effectively addressing the long-standing challenges in the field.

Structural Characterization of the Cotranscriptional Folding of the Thiamin Pyrophosphate Sensing thiC Riboswitch in Escherichia coli.

Riboswitches are RNA-regulating elements that mostly rely on structural changes to modulate gene expression at various levels. Recent studies have revealed that riboswitches may control several regulatory mechanisms cotranscriptionally, i.e., during the transcription elongation of the riboswitch or early in the coding region of the regulated gene. Here, we study the structure of the nascent thiamin pyrophosphate (TPP)-sensing thiC riboswitch in Escherichia coli by using biochemical and enzymatic conventional probing approaches. Our chemical (in-line and lead probing) and enzymatic (nucleases S1, A, T1, and RNase H) probing data provide a comprehensive model of how TPP binding modulates the structure of the thiC riboswitch. Furthermore, by using transcriptional roadblocks along the riboswitch sequence, we find that a certain portion of nascent RNA is needed to sense TPP that coincides with the formation of the P5 stem loop. Together, our data suggest that conventional techniques may readily be used to study cotranscriptional folding of nascent RNAs.

Partial RNA design.

MOTIVATION: RNA design is a key technique to achieve new functionality in fields like synthetic biology or biotechnology. Computational tools could help to find such RNA sequences but they are often limited in their formulation of the search space. RESULTS: In this work, we propose partial RNA design, a novel RNA design paradigm that addresses the limitations of current RNA design formulations. Partial RNA design describes the problem of designing RNAs from arbitrary RNA sequences and structure motifs with multiple design goals. By separating the design space from the objectives, our formulation enables the design of RNAs with variable lengths and desired properties, while still allowing precise control over sequence and structure constraints at individual positions. Based on this formulation, we introduce a new algorithm, libLEARNA, capable of efficiently solving different constraint RNA design tasks. A comprehensive analysis of various problems, including a realistic riboswitch design task, reveals the outstanding performance of libLEARNA and its robustness. AVAILABILITY AND IMPLEMENTATION: libLEARNA is open-source and publicly available at: https://github.com/automl/learna_tools.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration.

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single-molecule and bulk approaches, we discover how a single Mn2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the representative Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

Assessing the Quality of Cotranscriptional Folding Simulations.

Structural changes in RNAs are an important contributor to controlling gene expression not only at the posttranscriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5' region of the primary transcript regulates the downstream functional unit - usually an ORF - through premature termination of transcription. Not only such elements occur naturally, but they are also attractive devices in synthetic biology. The possibility to design such riboswitches or RNA thermometers is thus of considerable practical interest. Since these functional RNA elements act already during transcription, it is important to model and understand the dynamics of folding and, in particular, the formation of intermediate structures concurrently with transcription. Cotranscriptional folding simulations are therefore an important step to verify the functionality of design constructs before conducting expensive and labor-intensive wet lab experiments. For RNAs, full-fledged molecular dynamics simulations are far beyond practical reach because of both the size of the molecules and the timescales of interest. Even at the simplified level of secondary structures, further approximations are necessary. The BarMap approach is based on representing the secondary structure landscape for each individual transcription step by a coarse-grained representation that only retains a small set of low-energy local minima and the energy barriers between them. The folding dynamics between two transcriptional elongation steps is modeled as a Markov process on this representation. Maps between pairs of consecutive coarse-grained landscapes make it possible to follow the folding process as it changes in response to transcription elongation. In its original implementation, the BarMap software provides a general framework to investigate RNA folding dynamics on temporally changing landscapes. It is, however, difficult to use in particular for specific scenarios such as cotranscriptional folding. To overcome this limitation, we developed the user-friendly BarMap-QA pipeline described in detail in this contribution. It is illustrated here by an elaborate example that emphasizes the careful monitoring of several quality measures. Using an iterative workflow, a reliable and complete kinetics simulation of a synthetic, transcription-regulating riboswitch is obtained using minimal computational resources. All programs and scripts used in this contribution are free software and available for download as a source distribution for Linux® or as a platform-independent Docker® image including support for Apple macOS® and Microsoft Windows®.

Developing Complex RNA Design Applications in the Infrared Framework.

Applications in biotechnology and bio-medical research call for effective strategies to design novel RNAs with very specific properties. Such advanced design tasks require support by computational tools but at the same time put high demands on their flexibility and expressivity to model the application-specific requirements. To address such demands, we present the computational framework Infrared. It supports developing advanced customized design tools, which generate RNA sequences with specific properties, often in a few lines of Python code. This text guides the reader in tutorial format through the development of complex design applications. Thanks to the declarative, compositional approach of Infrared, we can describe this development as a step-by-step extension of an elementary design task. Thus, we start with generating sequences that are compatible with a single RNA structure and go all the way to RNA design targeting complex positive and negative design objectives with respect to single or even multiple target structures. Finally, we present a "real-world" application of computational design to create an RNA device for biotechnology: we use Infrared to generate design candidates of an artificial "AND" riboswitch, which activates gene expression in the simultaneous presence of two different small metabolites. In these applications, we exploit that the system can generate, in an efficient (fixed-parameter tractable) way, multiple diverse designs that satisfy a number of constraints and have high quality w.r.t. to an objective (by sampling from a Boltzmann distribution).

Switchable and orthogonal gene expression control inside artificial cells by synthetic riboswitches.

Here we report two novel synthetic riboswitches that respond to ASP2905 and theophylline and function in reconstituted cell-free protein synthesis (CFPS) system. We encapsulated the CFPS system as well as DNA-templated encoding reporter genes regulated by these orthogonal riboswitches inside liposomes, and achieved switchable and orthogonal control over gene expression by external stimulation with the cognate ligands.

Similarities and Differences in Ligand Binding to Protein and RNA Targets: The Case of Riboflavin.

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Fluorescent riboswitch-controlled biosensors for the genome scale analysis of metabolic pathways.

Fluorescent detection in cells has been tremendously developed over the years and now benefits from a large array of reporters that can provide sensitive and specific detection in real time. However, the intracellular monitoring of metabolite levels still poses great challenges due to the often complex nature of detected metabolites. Here, we provide a systematic analysis of thiamin pyrophosphate (TPP) metabolism in Escherichia coli by using a TPP-sensing riboswitch that controls the expression of the fluorescent gfp reporter. By comparing different combinations of reporter fusions and TPP-sensing riboswitches, we determine key elements that are associated with strong TPP-dependent sensing. Furthermore, by using the Keio collection as a proxy for growth conditions differing in TPP levels, we perform a high-throughput screen analysis using high-density solid agar plates. Our study reveals several genes whose deletion leads to increased or decreased TPP levels. The approach developed here could be applicable to other riboswitches and reporter genes, thus representing a framework onto which further development could lead to highly sophisticated detection platforms allowing metabolic screens and identification of orphan riboswitches.

Riboswitch and small RNAs modulate btuB translation initiation in Escherichia coli and trigger distinct mRNA regulatory mechanisms.

Small RNAs (sRNAs) and riboswitches represent distinct classes of RNA regulators that control gene expression upon sensing metabolic or environmental variations. While sRNAs and riboswitches regulate gene expression by affecting mRNA and protein levels, existing studies have been limited to the characterization of each regulatory system in isolation, suggesting that sRNAs and riboswitches target distinct mRNA populations. We report that the expression of btuB in Escherichia coli, which is regulated by an adenosylcobalamin (AdoCbl) riboswitch, is also controlled by the small RNAs OmrA and, to a lesser extent, OmrB. Strikingly, we find that the riboswitch and sRNAs reduce mRNA levels through distinct pathways. Our data show that while the riboswitch triggers Rho-dependent transcription termination, sRNAs rely on the degradosome to modulate mRNA levels. Importantly, OmrA pairs with the btuB mRNA through its central region, which is not conserved in OmrB, indicating that these two sRNAs may have specific targets in addition to their common regulon. In contrast to canonical sRNA regulation, we find that OmrA repression of btuB is lost using an mRNA binding-deficient Hfq variant. Together, our study demonstrates that riboswitch and sRNAs modulate btuB expression, providing an example of cis- and trans-acting RNA-based regulatory systems maintaining cellular homeostasis.

Characterization of the dual regulation by a c-di-GMP riboswitch Bc1 with a long expression platform from Bacillus thuringiensis.

A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.

Development and validation of a generic methyltransferase enzymatic assay based on an SAH riboswitch.

Methylation of proteins and nucleic acids plays a fundamental role in epigenetic regulation, and discovery of methyltransferase (MT) inhibitors is an area of intense activity. Because of the diversity of MTs and their products, assay methods that detect S-adenosylhomocysteine (SAH) - the invariant product of S-adenosylmethionine (SAM)-dependent methylation reactions - offer some advantages over methods that detect specific methylation events. However, direct, homogenous detection of SAH requires a reagent capable of discriminating between SAH and SAM, which differ by a single methyl group. Moreover, MTs are slow enzymes and many have submicromolar affinities for SAM; these properties translate to a need for detection of SAH at low nanomolar concentrations in the presence of excess SAM. To meet these needs, we leveraged the exquisite molecular recognition properties of a naturally occurring SAH-sensing RNA aptamer, or riboswitch. By splitting the riboswitch into two fragments, such that SAH binding induces assembly of a trimeric complex, we engineered sensors that transduce binding of SAH into positive fluorescence polarization (FP) and time resolved Förster resonance energy transfer (TR-FRET) signals. The split riboswitch configuration, called the AptaFluor™ SAH Methyltransferase Assay, allows robust detection of SAH (Z' > 0.7) at concentrations below 10 nM, with overnight signal stability in the presence of typical MT assay components. The AptaFluor assay tolerates diverse MT substrates, including histones, nucleosomes, DNA and RNA, and we demonstrated its utility as a robust, enzymatic assay method for several methyltransferases with SAM Km values < 1 µM. The assay was validated for HTS by performing a pilot screen of 1,280 compounds against the SARS-CoV-2 RNA capping enzyme, nsp14. By enabling direct, homogenous detection of SAH at low nanomolar concentrations, the AptaFluor assay provides a universal platform for screening and profiling MTs at physiologically relevant SAM concentrations.

Magnesium ions mitigate metastable states in the regulatory landscape of mRNA elements.

Residing in the 5' untranslated region of the mRNA, the 2'-deoxyguanosine (2'-dG) riboswitch mRNA element adopts an alternative structure upon binding of the 2'-dG molecule, which terminates transcription. RNA conformations are generally strongly affected by positively charged metal ions (especially Mg2+). We have quantitatively explored the combined effect of ligand (2'-dG) and Mg2+ binding on the energy landscape of the aptamer domain of the 2'-dG riboswitch with both explicit solvent all-atom molecular dynamics simulations (99 µsec aggregate sampling for the study) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments. We show that both ligand and Mg2+ are required for the stabilization of the aptamer domain; however, the two factors act with different modalities. The addition of Mg2+ remodels the energy landscape and reduces its frustration by the formation of additional contacts. In contrast, the binding of 2'-dG eliminates the metastable states by nucleating a compact core for the aptamer domain. Mg2+ ions and ligand binding are required to stabilize the least stable helix, P1 (which needs to unfold to activate the transcription platform), and the riboswitch core formed by the backbone of the P2 and P3 helices. Mg2+ and ligand also facilitate a more compact structure in the three-way junction region.

Riboswitch-controlled IL-12 gene therapy reduces hepatocellular cancer in mice.

Hepatocellular carcinoma (HCC) and solid cancers with liver metastases are indications with high unmet medical need. Interleukin-12 (IL-12) is a proinflammatory cytokine with substantial anti-tumor properties, but its therapeutic potential has not been realized due to severe toxicity. Here, we show that orthotopic liver tumors in mice can be treated by targeting hepatocytes via systemic delivery of adeno-associated virus (AAV) vectors carrying the murine IL-12 gene. Controlled cytokine production was achieved in vivo by using the tetracycline-inducible K19 riboswitch. AAV-mediated expression of IL-12 led to STAT4 phosphorylation, interferon-γ (IFNγ) production, infiltration of T cells and, ultimately, tumor regression. By detailed analyses of efficacy and tolerability in healthy and tumor-bearing animals, we could define a safe and efficacious vector dose. As a potential clinical candidate, we characterized vectors carrying the human IL-12 (huIL-12) gene. In mice, bioactive human IL-12 was expressed in a vector dose-dependent manner and could be induced by tetracycline, suggesting tissue-specific AAV vectors with riboswitch-controlled expression of highly potent proinflammatory cytokines as an attractive approach for vector-based cancer immunotherapy.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Riboswitch Distribution in the Human Gut Microbiome Reveals Common Metabolite Pathways.

Riboswitches are widely distributed, conserved RNAs which regulate metabolite levels in bacterial cells through direct, noncovalent binding of their cognate metabolite. Various riboswitch families are highly enriched in gut bacteria, suggestive of a symbiotic relationship between the host and bacteria. Previous studies of the distribution of riboswitches have examined bacterial taxa broadly. Thus, the distribution of riboswitches associated with bacteria inhabiting the intestines of healthy individuals is not well understood. To address these questions, we survey the gut microbiome for riboswitches by including an international database of prokaryotic genomes from the gut samples. Using Infernal, a program that uses RNA-specific sequence and structural features, we survey this data set using existing riboswitch models. We identify 22 classes of riboswitches with vitamin cofactors making up the majority of riboswitch-associated pathways. Our finding is reproducible in other representative databases from the oral as well as the marine microbiomes, underscoring the importance of thiamine pyrophosphate, cobalamin, and flavin mononucleotide in gene regulation. Interestingly, riboswitches do not vary significantly across microbiome representatives from around the world despite major taxonomic differences; this suggests an underlying conservation. Further studies elucidating the role of bacterial riboswitches in the host metabolome are needed to illuminate the consequences of our finding.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch.

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Direct fluoride monitoring using a fluorogenic RNA-based biosensor.

Fluorinated compounds, whether naturally occurring or from anthropogenic origin, have been extensively exploited in the last century. Degradation of these compounds by physical or biochemical processes is expected to result in the release of fluoride. Several fluoride detection mechanisms have been previously described. However, most of these methods are not compatible with high- and ultrahigh-throughput screening technologies, lack the ability to real-time monitor the increase of fluoride concentration in solution, or rely on costly reagents (such as cell-free expression systems). Our group recently developed "FluorMango" as the first completely RNA-based and direct fluoride-specific fluorogenic biosensor. To do so, we merged and engineered the Mango-III light-up RNA aptamer and the fluoride-specific aptamer derived from a riboswitch, crcB. In this chapter, we explain how this RNA-based biosensor can be produced in large scale before providing examples of how it can be used to quantitatively detect (end-point measurement) or monitor in real-time fluoride release in complex biological systems by translating it into measurable fluorescent signal.

Resolving the intricate binding of neomycin B to multiple binding motifs of a neomycin-sensing riboswitch aptamer by native top-down mass spectrometry and NMR spectroscopy.

Understanding small molecule binding to RNA can be complicated by an intricate interplay between binding stoichiometry, multiple binding motifs, different occupancies of different binding motifs, and changes in the structure of the RNA under study. Here, we use native top-down mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy to experimentally resolve these factors and gain a better understanding of the interactions between neomycin B and the 40 nt aptamer domain of a neomycin-sensing riboswitch engineered in yeast. Data from collisionally activated dissociation of the 1:1, 1:2 and 1:3 RNA-neomycin B complexes identified a third binding motif C of the riboswitch in addition to the two motifs A and B found in our previous study, and provided occupancies of the different binding motifs for each complex stoichiometry. Binding of a fourth neomycin B molecule was unspecific according to both MS and NMR data. Intriguingly, all major changes in the aptamer structure can be induced by the binding of the first neomycin B molecule regardless of whether it binds to motif A or B as evidenced by stoichiometry-resolved MS data together with titration data from 1H NMR spectroscopy in the imino proton region. Specific binding of the second and third neomycin B molecules further stabilizes the riboswitch aptamer, thereby allowing for a gradual response to increasing concentrations of neomycin B, which likely leads to a fine-tuning of the cellular regulatory mechanism.

The discovery of novel noncoding RNAs in 50 bacterial genomes.

Structured noncoding RNAs (ncRNAs) contribute to many important cellular processes involving chemical catalysis, molecular recognition and gene regulation. Few ncRNA classes are broadly distributed among organisms from all three domains of life, but the list of rarer classes that exhibit surprisingly diverse functions is growing. We previously developed a computational pipeline that enables the near-comprehensive identification of structured ncRNAs expressed from individual bacterial genomes. The regions between protein coding genes are first sorted based on length and the fraction of guanosine and cytidine nucleotides. Long, GC-rich intergenic regions are then examined for sequence and structural similarity to other bacterial genomes. Herein, we describe the implementation of this pipeline on 50 bacterial genomes from varied phyla. More than 4700 candidate intergenic regions with the desired characteristics were identified, which yielded 44 novel riboswitch candidates and numerous other putative ncRNA motifs. Although experimental validation studies have yet to be conducted, this rate of riboswitch candidate discovery is consistent with predictions that many hundreds of novel riboswitch classes remain to be discovered among the bacterial species whose genomes have already been sequenced. Thus, many thousands of additional novel ncRNA classes likely remain to be discovered in the bacterial domain of life.