The structure of a tetrahydrofolate-sensing riboswitch reveals two ligand binding sites in a single aptamer

Structure-based characterization and compound identification of the wild-type THF class-II riboswitch

Riboswitches are conserved regulatory RNA elements participating in various metabolic pathways. Recently, a novel RNA motif known as the folE RNA motif was discovered upstream of folE genes. It specifically senses tetrahydrofolate (THF) and is therefore termed THF-II riboswitch. To unravel the ligand recognition mechanism of this newly discovered riboswitch and decipher the underlying principles governing its tertiary folding, we determined both the free-form and bound-form THF-II riboswitch in the wild-type sequences. Combining structural information and isothermal titration calorimetry (ITC) binding assays on structure-based mutants, we successfully elucidated the significant long-range interactions governing the function of THF-II riboswitch and identified additional compounds, including alternative natural metabolites and potential lead compounds for drug discovery, that interact with THF-II riboswitch. Our structural research on the ligand recognition mechanism of the THF-II riboswitch not only paves the way for identification of compounds targeting riboswitches, but also facilitates the exploration of THF analogs in diverse biological contexts or for therapeutic applications.

A spermidine riboswitch class in bacteria exploits a close variant of an aptamer for the enzyme cofactor S-adenosylmethionine

Natural polyamines such as spermidine and spermine cations have characteristics that make them highly likely to be sensed by riboswitches, such as their general affinity to polyanionic RNA and their broad contributions to cell physiology. Despite previous claims that polyamine riboswitches exist, evidence of their biological functions has remained unconvincing. Here, we report that rare variants of bacterial S-adenosylmethionine-I (SAM-I) riboswitches reject SAM and have adapted to selectively sense spermidine. These spermidine-sensing riboswitch variants are associated with genes whose protein products are directly involved in the production of spermidine and other polyamines. Biochemical and genetic assays demonstrate that representatives of this riboswitch class robustly function as genetic "off" switches, wherein spermidine binding causes premature transcription termination to suppress the expression of polyamine biosynthetic genes. These findings confirm the existence of natural spermidine-sensing riboswitches in bacteria and expand the list of variant riboswitch classes that have adapted to bind different ligands.

Dynamical characterization and multiple unbinding paths of two PreQ1 ligands in one pocket

Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.

Riboswitches, from cognition to transformation

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

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.

Context-dependence of T-loop Mediated Long-range RNA Tertiary Interactions

The architecture and folding of complex RNAs is governed by a limited set of highly recurrent structural motifs that form long-range tertiary interactions. One of these motifs is the T-loop, which was first identified in tRNA but is broadly distributed across biological RNAs. While the T-loop has been examined in detail in different biological contexts, the various receptors that it interacts with are not as well defined. In this study, we use a cell-based genetic screen in concert with bioinformatic analysis to examine three different, but related, T-loop receptor motifs found in the flavin mononucleotide (FMN) and cobalamin (Cbl) riboswitches. As a host for different T-loop receptors, we employed the env8 class-II Cbl riboswitch, an RNA that uses two T-loop motifs for both folding and supporting the ligand binding pocket. A set of libraries was created in which select nucleotides that participate in the T-loop/T-loop receptor (TL/TLR) interaction were fully randomized. Library members were screened for their ability to support Cbl-dependent expression of a reporter gene. While T-loops appear to be variable in sequence, we find that the functional sequence space is more restricted in the Cbl riboswitch, suggesting that TL/TLR interactions are context dependent. Our data reveal clear sequence signatures for the different types of receptor motifs that align with phylogenic analysis of these motifs in the FMN and Cbl riboswitches. Finally, our data suggest the functional contribution of various nucleobase-mediated long-range interactions within the riboswitch subclass of TL/TLR interactions that are distinct from those found in other RNAs.

Crystal structures of the NAD+-II riboswitch reveal two distinct ligand-binding pockets

We present crystal structures of a new NAD+-binding riboswitch termed NAD+-II, bound to nicotinamide mononucleotide (NMN), nicotinamide adenine dinucleotide (NAD+) and nicotinamide riboside (NR). The RNA structure comprises a number of structural features including three helices, one of which forms a triple helix by interacting with an A5 strand in its minor-groove, and another formed from a long-range pseudoknot. The core of the structure (centrally located and coaxial with the triplex and the pseudoknot) includes two consecutive quadruple base interactions. Unusually the riboswitch binds two molecules of ligand, bound at distinct, non-overlapping sites in the RNA. Binding occurs primarily through the nicotinamide moiety of each ligand, held by specific hydrogen bonding and stacking interactions with the pyridyl ring. The mode of binding is the same for NMN, NR and the nicotinamide moiety of NAD+. In addition, when NAD+ is bound into one site it adopts an elongated conformation such that its diphosphate linker occupies a groove on the surface of the RNA, following which the adenine portion inserts into a pocket and makes specific hydrogen bonding interactions. Thus the NAD+-II riboswitch is distinct from the NAD+-I riboswitch in that it binds two molecules of ligand at separate sites, and that binding occurs principally through the nicotinamide moiety.

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

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

Structural insights into translation regulation by the THF-II riboswitch

In bacteria, expression of folate-related genes is controlled by the tetrahydrofolate (THF) riboswitch in response to specific binding of THF and its derivatives. Recently, a second class of THF riboswitches, named THF-II, was identified in Gram-negative bacteria, which exhibit distinct architecture from the previously characterized THF-I riboswitches found in Gram-positive bacteria. Here, we present the crystal structures of the ligand-bound THF-II riboswitch from Mesorhizobium loti. These structures exhibit a long rod-like fold stabilized by continuous base pair and base triplet stacking across two helices of P1 and P2 and their interconnecting ligand-bound binding pocket. The pterin moiety of the ligand docks into the binding pocket by forming hydrogen bonds with two highly conserved pyrimidines in J12 and J21, which resembles the hydrogen-bonding pattern at the ligand-binding site FAPK in the THF-I riboswitch. Using small-angle X-ray scattering and isothermal titration calorimetry, we further characterized the riboswitch in solution and reveal that Mg2+ is essential for pre-organization of the binding pocket for efficient ligand binding. RNase H cleavage assay indicates that ligand binding reduces accessibility of the ribosome binding site in the right arm of P1, thus down-regulating the expression of downstream genes. Together, these results provide mechanistic insights into translation regulation by the THF-II riboswitch.

Isothermal Titration Calorimetry Analysis of a Cooperative Riboswitch Using an Interdependent-Sites Binding Model

Isothermal titration calorimetry (ITC) is a powerful biophysical tool to characterize energetic profiles of biomacromolecular interactions without any alteration of the underlying chemical structures. In this protocol, we describe procedures for performing, analyzing, and interpreting ITC data obtained from a cooperative riboswitch-ligand interaction.

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.

Key players in regulatory RNA realm of bacteria

Precise regulation of gene expression is crucial for living cells to adapt for survival in diverse environmental conditions. Among the common cellular regulatory mechanisms, RNA-based regulators play a key role in all domains of life. Discovery of regulatory RNAs have made a paradigm shift in molecular biology as many regulatory functions of RNA have been identified beyond its canonical roles as messenger, ribosomal and transfer RNA. In the complex regulatory RNA network, riboswitches, small RNAs, and RNA thermometers can be identified as some of the key players. Herein, we review the discovery, mechanism, and potential therapeutic use of these classes of regulatory RNAs mainly found in bacteria. Being highly adaptive organisms that inhabit a broad range of ecological niches, bacteria have adopted tight and rapid-responding gene regulation mechanisms. This review aims to highlight how bacteria utilize versatile RNA structures and sequences to build a sophisticated gene regulation network.

•

The three major classes of prokaryotic ncRNAs and their characterized mechanisms of operation in gene regulation.

•

sRNAs emerging as major players in global gene regulatory networks.

•

Riboswitch mediated gene control mechanisms through on/off switches in response to ligand binding.

•

RNA thermo sensors for temperature-dependent gene expression.

•

Therapeutic importance of ncRNAs and computational approaches involved in the discovery of ncRNAs.

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

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

A small RNA that cooperatively senses two stacked metabolites in one pocket for gene control

Riboswitches are structured non-coding RNAs often located upstream of essential genes in bacterial messenger RNAs. Such RNAs regulate expression of downstream genes by recognizing a specific cellular effector. Although nearly 50 riboswitch classes are known, only a handful recognize multiple effectors. Here, we report the 2.60-Å resolution co-crystal structure of a class I type I preQ₁-sensing riboswitch that reveals two effectors stacked atop one another in a single binding pocket. These effectors bind with positive cooperativity in vitro and both molecules are necessary for gene regulation in bacterial cells. Stacked effector recognition appears to be a hallmark of the largest subgroup of preQ₁ riboswitches, including those from pathogens such as Neisseria gonorrhoeae. We postulate that binding to stacked effectors arose in the RNA World to closely position two substrates for RNA-mediated catalysis. These findings expand known effector recognition capabilities of riboswitches and have implications for antimicrobial development.

Riboswitches contain an aptamer domain that recognizes a metabolite and an expression platform that regulates gene expression. Here the authors report the crystal structure of a preQ₁-sensing riboswitch from Carnobacterium antarcticus that shows two metabolites in a single binding pocket.

Organic Bioelectronic Devices for Metabolite Sensing

Electrochemical detection of metabolites is essential for early diagnosis and continuous monitoring of a variety of health conditions. This review focuses on organic electronic material-based metabolite sensors and highlights their potential to tackle critical challenges associated with metabolite detection. We provide an overview of the distinct classes of organic electronic materials and biorecognition units used in metabolite sensors, explain the different detection strategies developed to date, and identify the advantages and drawbacks of each technology. We then benchmark state-of-the-art organic electronic metabolite sensors by categorizing them based on their application area (in vitro, body-interfaced, in vivo, and cell-interfaced). Finally, we share our perspective on using organic bioelectronic materials for metabolite sensing and address the current challenges for the devices and progress to come.

Riboswitch Mechanisms: New Tricks for an Old Dog

Discovered almost twenty years ago, riboswitches turned out to be one of the most common regulatory systems in bacteria, with representatives found in eukaryotes and archaea. Unlike many other regulatory elements, riboswitches are entirely composed of RNA and capable of modulating expression of genes by direct binding of small cellular molecules. While bacterial riboswitches had been initially thought to control production of enzymes and transporters associated with small organic molecules via feedback regulatory circuits, later findings identified riboswitches directing expression of a wide range of genes and responding to various classes of molecules, including ions, signaling molecules, and others. The 5'-untranslated mRNA regions host a vast majority of riboswitches, which modulate transcription or translation of downstream genes through conformational rearrangements in the ligand-sensing domains and adjacent expression-controlling platforms. Over years, the repertoire of regulatory mechanisms employed by riboswitches has greatly expanded; most recent studies have highlighted the importance of alternative mechanisms, such as RNA degradation, for the riboswitch-mediated genetic circuits. This review discusses the plethora of bacterial riboswitch mechanisms and illustrates how riboswitches utilize different features and approaches to elicit various regulatory responses.

Synthetic biology-driven microbial production of folates: Advances and perspectives

With the development and application of synthetic biology, significant progress has been made in the production of folate by microbial fermentation using cell factories, especially for using generally regarded as safe (GRAS) microorganism as production host. In this review, the physiological functions and applications of folates were firstly discussed. Second, the current advances of folate-producing GRAS strains development were summarized. Third, the applications of synthetic biology-based metabolic regulatory tools in GRAS strains were introduced, and the progress in the application of these tools for folate production were summarized. Finally, the challenges to folates efficient production and corresponding emerging strategies to overcome them by synthetic biology were discussed, including the construction of biosensors using tetrahydrofolate riboswitches to regulate metabolic pathways, adaptive evolution to overcome the flux limitations of the folate pathway. The combination of new strategies and tools of synthetic biology is expected to further improve the efficiency of microbial folate synthesis.

Tying the knot in the tetrahydrofolate (THF) riboswitch: A molecular basis for gene regulation

Effective gene regulation by the tetrahydrofolate riboswitch depends not only on ligand affinity but also on the kinetics of ligand association, which involves two cooperative binding sites. We have determined a 1.9-Å resolution crystal structure of the ligand-free THF riboswitch aptamer. The pseudoknot binding site 'unwinds' in the absence of ligand, whereby the adjacent helical domains (P1, P2, and P3) become disjointed, resulting in rotation and misalignment of the gene-regulatory P1 helix with respect to P3. In contrast, the second binding site at the three-way junction, which is the first to fold, is structurally conserved between apo and holo forms. This suggests a kinetic role for this site, in which binding of the first ligand molecule to the stably folded three-way junction promotes formation of the regulatory pseudoknot site and subsequent binding of the second molecule. As such, these findings provide a molecular basis for both conformational switching and kinetic control.

Cooperativity and Allostery in RNA Systems

Allostery is among the most basic biological principles employed by biological macromolecules to achieve a biologically active state in response to chemical cues. Although initially used to describe the impact of small molecules on the conformation and activity of protein enzymes, the definition of this term has been significantly broadened to describe long-range conformational change of macromolecules in response to small or large effectors. Such a broad definition could be applied to RNA molecules, which do not typically serve as protein-free cellular enzymes but fold and form macromolecular assemblies with the help of various ligand molecules, including ions and proteins. Ligand-induced allosteric changes in RNA molecules are often accompanied by cooperative interactions between RNA and its ligand, thus streamlining the folding and assembly pathways. This chapter provides an overview of the interplay between cooperativity and allostery in RNA systems and outlines methods to study these two biological principles.

Thermodynamic analysis of cooperative ligand binding by the ATP-binding DNA aptamer indicates a population-shift binding mechanism

The ATP-binding DNA aptamer is often used as a model system for developing new aptamer-based biosensor methods. This aptamer follows a structure-switching binding mechanism and is unusual in that it binds two copies of its ligand. We have used isothermal titration calorimetry methods to study the binding of ATP, ADP, AMP and adenosine to the ATP-binding aptamer. Using both individual and global fitting methods, we show that this aptamer follows a positive cooperative binding mechanism. We have determined the binding affinity and thermodynamics for both ligand-binding sites. By separating the ligand-binding sites by an additional four base pairs, we engineered a variant of this aptamer that binds two adenosine ligands in an independent manner. Together with NMR and thermal stability experiments, these data indicate that the ATP-binding DNA aptamer follows a population-shift binding mechanism that is the source of the positive binding cooperativity by the aptamer.

Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs)

Mechanical anisotropy is an essential property for many biomolecules to assume their structures, functions and applications, however, the mechanisms for their direction-dependent mechanical responses remain elusive. Herein, by using a single-molecule nanopore sensing technique, we explore the mechanisms of directional mechanical stability of the xrRNA1 RNA from ZIKA virus (ZIKV), which forms a complex ring-like architecture. We reveal extreme mechanical anisotropy in ZIKV xrRNA1 which highly depends on Mg²⁺ and the key tertiary interactions. The absence of Mg²⁺ and disruption of the key tertiary interactions strongly affect the structural integrity and attenuate mechanical anisotropy. The significance of ring structures in RNA mechanical anisotropy is further supported by steered molecular dynamics simulations in combination with force distribution analysis. We anticipate the ring structures can be used as key elements to build RNA-based nanostructures with controllable mechanical anisotropy for biomaterial and biomedical applications.

Small molecule targeting of RNA structures in neurological disorders

Aberrant RNA structure and function operate in neurological disease progression and severity. As RNA contributes to disease pathology in a complex fashion, that is, via various mechanisms, it has become an attractive therapeutic target for small molecules and oligonucleotides. In this review, we discuss the identification of RNA structures that cause or contribute to neurological diseases as well as recent progress toward the development of small molecules that target them, including small molecule modulators of pre-mRNA splicing and RNA repeat expansions that cause microsatellite disorders such as Huntington's disease and amyotrophic lateral sclerosis. The use of oligonucleotide-based modalities is also discussed. There are key differences between small molecule and oligonucleotide targeting of RNA. The former targets RNA structure, while the latter prefers unstructured regions. Thus, some targets will be preferentially targeted by oligonucleotides and others by small molecules.

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.

High Affinity Binding of N2-Modified Guanine Derivatives Significantly Disrupts the Ligand Binding Pocket of the Guanine Riboswitch

Riboswitches are important model systems for the development of approaches to search for RNA-targeting therapeutics. A principal challenge in finding compounds that target riboswitches is that the effector ligand is typically almost completely encapsulated by the RNA, which severely limits the chemical space that can be explored. Efforts to find compounds that bind the guanine/adenine class of riboswitches with a high affinity have in part focused on purines modified at the C6 and C2 positions. These studies have revealed compounds that have low to sub-micromolar affinity and, in a few cases, have antimicrobial activity. To further understand how these compounds interact with the guanine riboswitch, we have performed an integrated structural and functional analysis of representative guanine derivatives with modifications at the C8, C6 and C2 positions. Our data indicate that while modifications of guanine at the C6 position are generally unfavorable, modifications at the C8 and C2 positions yield compounds that rival guanine with respect to binding affinity. Surprisingly, C2-modified guanines such as N2-acetylguanine completely disrupt a key Watson–Crick pairing interaction between the ligand and RNA. These compounds, which also modulate transcriptional termination as efficiently as guanine, open up a significant new chemical space of guanine modifications in the search for antimicrobial agents that target purine riboswitches.

Structural basis for 2'-deoxyguanosine recognition by the 2'-dG-II class of riboswitches

A recent bioinformatic analysis of well-characterized classes of riboswitches uncovered subgroups unable to bind to the regulatory molecule of the parental class. Within the guanine/adenine class, seven groups of RNAs were identified that deviate from the consensus sequence at one or more of three positions directly involved purine nucleobase recognition, one of which was validated as a second class of 2'-deoxyguanosine riboswitch (called 2'-dG-II). To understand how 2'-dG-II riboswitches recognize their cognate ligand and how they differ from a previously identified class of 2'-deoxyguanosine binding riboswitches, we have solved the crystal structure of a 2'-dG-II aptamer domain bound to 2'-deoxyguanosine. This structure reveals a global architecture similar to other members of the purine riboswitch family, but contains key differences within the ligand binding core. Defining the 2'-dG-II riboswitches is a two-nucleotide insertion in the three-way junction that promotes novel base-base interactions. Unlike 2'-dG-I riboswitches, the 2'-dG-II class only requires local changes to the ligand binding pocket of the guanine/adenine class to achieve a change in ligand preference. Notably, members of the 2'-dG-II family have variable ability to discriminate between 2'-deoxyguanosine and riboguanosine, suggesting that a subset of 2'-dG-II riboswitches may bind either molecule to regulate gene expression.

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.

Mathematical approaches in estimating aptamer-target binding affinity

The performance of aptamers as versatile tools in numerous analytical applications is critically dependent on their high target binding specificity and selectivity. However, only the technical or methodological aspects of measuring aptamer-target binding affinities are focused, ignoring the equally important mathematical components that play pivotal roles in affinity measurements. In this study, we aim to provide a comprehensive review regarding the utilization of different mathematical models and equations, along with a detailed description of the computational steps involved in mathematically deriving the binding affinity of aptamers against their specific target molecules. Mathematical models ranging from one-site binding to multiple aptameric binding site-based models are explained in detail. Models applied in several different approaches of affinity measurements such as thermodynamics and kinetic analysis, including cooperativity and competitive-assay based mathematical models have been elaborately discussed. Mathematical models incorporating factors that could potentially affect affinity measurements are also further scrutinized.

Responsive self-assembly of tectoRNAs with loop-receptor interactions from the tetrahydrofolate (THF) riboswitch

Naturally occurring RNAs are known to exhibit a high degree of modularity, whereby specific structural modules (or motifs) can be mixed and matched to create new molecular architectures. The modular nature of RNA also affords researchers the ability to characterize individual structural elements in controlled synthetic contexts in order to gain new and critical insights into their particular structural features and overall performance. Here, we characterized the binding affinity of a unique loop–receptor interaction found in the tetrahydrofolate (THF) riboswitch using rationally designed self-assembling tectoRNAs. Our work suggests that the THF loop–receptor interaction has been fine-tuned for its particular role as a riboswitch component. We also demonstrate that the thermodynamic stability of this interaction can be modulated by the presence of folinic acid, which induces a local structural change at the level of the loop–receptor. This corroborates the existence of a THF binding site within this tertiary module and paves the way for its potential use as a THF responsive module for RNA nanotechnology and synthetic biology.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

Employing a ZTP Riboswitch to Detect Bacterial Folate Biosynthesis Inhibitors in a Small Molecule High-Throughput Screen

Various riboswitch classes are being discovered that precisely monitor the status of important biological processes, including metabolic pathway function, signaling for physiological adaptations, and responses to toxic agents. Biochemical components for some of these processes might make excellent targets for the development of novel antibacterial molecules, which can be broadly sought by using phenotypic drug discovery (PDD) methods. However, PDD data do not normally provide clues regarding the target for each hit compound. We have developed and validated a robust fluorescent reporter system based on a ZTP riboswitch that identifies numerous folate biosynthesis inhibitors with high sensitivity and precision. The utility of the riboswitch-based PDD strategy was evaluated using Escherichia coli bacteria by conducting a 128 310-compound high-throughput screen, which identified 78 sulfanilamide derivatives among the many initial hits. Similarly, representatives of other riboswitch classes could be employed to rapidly match antibacterial hits with the biological processes they target.

Biochemical validation of a second class of tetrahydrofolate riboswitches in bacteria

We previously reported a large collection of structured noncoding RNAs (ncRNAs) that includes many riboswitch candidates identified through comparative sequence analysis of bacterial intergenic regions. One of these candidates, initially named the "folE motif," adopts a simple architecture commonly found upstream of folE genes. FolE enzymes catalyze the first enzyme in the de novo folate biosynthesis pathway. Herein, we demonstrate that folE motif RNAs selectively bind the enzyme cofactor tetrahydrofolate (THF) and several of its close derivatives. These aptamers, commonly found in Gram-negative bacteria, are distinct from aptamers of the previous validated THF riboswitch class found in Gram-positive bacteria. Our findings indicate that folE motif RNAs are aptamer domains for a second THF riboswitch class, named THF-II. The biochemical validation of THF-II riboswitches further highlights the ability of bacteria to utilize diverse RNA structures to sense universal enzyme cofactors that are predicted to be of ancient origin.

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

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

Single-Molecule FRET Kinetics of the Mn2+ Riboswitch: Evidence for Allosteric Mg2+ Control of "Induced-Fit" vs "Conformational Selection" Folding Pathways

Gene expression in bacteria is often regulated dynamically by conformational changes in a riboswitch upon ligand binding, a detailed understanding of which is very much in its infancy. For example, the manganese riboswitch is a widespread RNA motif that conformationally responds in regulating bacterial gene expression to micromolar levels of its eponymous ligand, Mn²⁺, but the mechanistic pathways are poorly understood. In this work, we quantitatively explore the dynamic folding behavior of the manganese riboswitch by single-molecule fluorescence resonance energy transfer spectroscopy as a function of cation/ligand conditions. From the detailed analysis of the kinetics, the Mn²⁺ is shown to fold the riboswitch by a "bind-then-fold" (i.e., "induced-fit", IF) mechanism, whereby the ligand binds first and then promotes folding. On the other hand, the data also clearly reveal the presence of a folded yet ligand-free structure predominating due to the addition of physiological Mg²⁺ to a nonselective metal ion binding site. Of particular kinetic interest, such a Mg²⁺ "prefolded" conformation of the riboswitch is shown to exhibit a significantly increased affinity for Mn²⁺ and further stabilization by subsequent binding of the ligand, thereby promoting efficient riboswitch folding by a "fold-then-bind" (i.e., "conformational selection", CS) mechanism. Our results not only demonstrate Mg²⁺-controlled switching between IF and CS riboswitch folding pathways but also suggest a novel heterotropic allosteric control in the manganese riboswitch activity co-regulated by Mg²⁺ binding.

Small-Molecule-Binding Riboswitches

RNA is a versatile biomolecule capable of transferring information, taking on distinct three-dimensional shapes, and reacting to ambient conditions. RNA molecules utilize a wide range of mechanisms to control gene expression. An example of such regulation is riboswitches. Consisting exclusively of RNA, they are able to control important metabolic processes, thus providing an elegant and efficient RNA-only regulation system. Existing across all domains of life, riboswitches appear to represent one of the most highly conserved mechanisms for the regulation of a broad range of biochemical pathways. Through binding of a wide range of small-molecule ligands to their so-called aptamer domain, riboswitches undergo a conformational change in their downstream "expression platform." In consequence, the pattern of gene expression changes, which in turn results in increased or decreased protein production. Riboswitches unite the sensing and transduction of a signal that can directly be coupled to the metabolism of the cell; thus they constitute a very potent regulatory mechanism for many organisms. Highly specific RNA-binding domains not only occur in vivo but can also be evolved by means of the SELEX (systematic evolution of ligands by exponential enrichment) method, which allows in vitro selection of aptamers against almost any ligand. Coupling of these aptamers with an expression platform has led to the development of synthetic riboswitches, a highly active research field of great relevance and immense potential. The aim of this review is to summarize developments in the riboswitch field over the last decade and address key questions of recent research.

Riboswitch regulation in cyanobacteria is independent of their habitat adaptations

Cyanobacteria are one of the ancient bacterial species occupying a variety of habitats with diverse metabolic preferences. RNA regulators like riboswitches play significant role in controlling the gene expression in prokaryotes. The taxonomic distribution of riboswitches suggests that they might be one of the oldest mechanisms of gene control system. In this paper, we analyzed the distribution of different riboswitch families in various cyanobacterial genomes. It was observed that only four riboswitch classes were abundant in cyanobacteria, B12-element (Cob)/AdoCbl/AdoCbl-variant riboswitch being the most abundant. The analysis suggests that riboswitch mode of regulation is present in cyanobacterial species irrespective of their habitat types. A large number of unidentified genes regulated by riboswitches listed in this analysis indicate the wide range of targets for these riboswitch families. The analysis revealed a large number of genes regulated by riboswitches which may assist in elaborating the diversity among the cyanobacterial species.

Synthesis of All Possible Canonical (3'-5'-Linked) Cyclic Dinucleotides and Evaluation of Riboswitch Interactions and Immune-Stimulatory Effects

The cyclic dinucleotides (CDNs) c-di-GMP, c-di-AMP, and c-AMP-GMP are widely utilized as second messengers in bacteria, where they signal lifestyle changes such as motility and biofilm formation, cell wall and membrane homeostasis, virulence, and exo-electrogenesis. For all known bacterial CDNs, specific riboswitches have been identified that alter gene expression in response to the second messengers. In addition, bacterial CDNs trigger potent immune responses, making them attractive as adjuvants in immune therapies. Besides the three naturally occurring CDNs, seven further CDNs containing canonical 3'-5'-linkages are possible by combining the four natural ribonucleotides. Herein, we have synthesized all ten possible combinations of 3'-5'-linked CDNs. The binding affinity of novel CDNs and GEMM riboswitch variants was assessed utilizing a spinach aptamer fluorescence assay and in-line probing assays. The immune-stimulatory effect of CDNs was evaluated by induction of type I interferons (IFNs), and a novel CDN c-AMP-CMP was identified as a new immune-stimulatory agent.

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.

Allosteric pathways in tetrahydrofolate sensing riboswitch with dynamics correlation network

Riboswitches are cis-acting genetic control elements. Due to their fundamental importance in bacteria gene regulation, they have been proposed as antibacterial drug targets. Tetrahydrofolate (THF) is an essential cofactor of one-carbon transfer reactions and downregulates the expression of downstream genes. However, information on how to transfer from the binding site of THF to the expression platform is still unavailable. Herein, a nucleotide/nucleotide dynamics correlation network based on an all-atom molecular dynamic simulation was used to reveal the regulation mechanism of a THF-sensing riboswitch. Shortest pathway analysis based on the network illustrates that there is an allosteric pathway through the P2 helix to the pseudoknot, then to the P1 helix in the THF-riboswitch. Thus the hypothesis of "THF-binding induced allosteric switching" was proposed and evaluated using THF and pseudoknot weakened experiments. Furthermore, a possible allosteric pathway of C30-C31-G33-A34-G35-G36-G37-A38-G48-G47-U46-A90-U91-C92-G93-C94-G95-C96 was identified and confirmed through the perturbation of the network. The proposed allosteric mechanism and the underlying allosteric pathway provide fundamental insights for the regulation of THF sensing riboswitches.

Riboswitch diversity and distribution

Riboswitches are commonly used by bacteria to detect a variety of metabolites and ions to regulate gene expression. To date, nearly 40 different classes of riboswitches have been discovered, experimentally validated, and modeled at atomic resolution in complex with their cognate ligands. The research findings produced since the first riboswitch validation reports in 2002 reveal that these noncoding RNA domains exploit many different structural features to create binding pockets that are extremely selective for their target ligands. Some riboswitch classes are very common and are present in bacteria from nearly all lineages, whereas others are exceedingly rare and appear in only a few species whose DNA has been sequenced. Presented herein are the consensus sequences, structural models, and phylogenetic distributions for all validated riboswitch classes. Based on our findings, we predict that there are potentially many thousands of distinct bacterial riboswitch classes remaining to be discovered, but that the rarity of individual undiscovered classes will make it increasingly difficult to find additional examples of this RNA-based sensory and gene control mechanism.

The lost language of the RNA World

The possibility of an RNA World is based on the notion that life on Earth passed through a primitive phase without proteins, a time when all genomes and enzymes were composed of ribonucleic acids. Numerous apparent vestiges of this ancient RNA World remain today, including many nucleotide-derived coenzymes, self-processing ribozymes, metabolite-binding riboswitches, and even ribosomes. Many of the most common signaling molecules and second messengers used by modern organisms are also formed from RNA nucleotides or their precursors. For example, nucleotide derivatives such as cAMP, ppGpp, and ZTP, as well as the cyclic dinucleotides c-di-GMP and c-di-AMP, are intimately involved in signaling diverse physiological or metabolic changes in bacteria and other organisms. We describe the potential diversity of this "lost language" of the RNA World and speculate on whether additional components of this ancient communication machinery might remain hidden though still very much relevant to modern cells.

Naringenin-responsive riboswitch-based fluorescent biosensor module for Escherichia coli co-cultures

The ability to design and construct combinatorial synthetic metabolic pathways has far exceeded our capacity for efficient screening and selection of the resulting microbial strains. The need for high-throughput rapid screening techniques is of upmost importance for the future of synthetic biology and metabolic engineering. Here we describe the development of an RNA riboswitch-based biosensor module with dual fluorescent reporters, and demonstrate a high-throughput flow cytometry-based screening method for identification of naringenin over producing Escherichia coli strains in co-culture. Our efforts helped identify a number of key operating parameters that affect biosensor performance, including the selection of promoter and linker elements within the sensor-actuator domain, and the effect of host strain, fermentation time, and growth medium on sensor dynamic range. The resulting biosensor demonstrates a high correlation between specific fluorescence of the biosensor strain and naringenin titer produced by the second member of the synthetic co-culture system. This technique represents a novel application for synthetic microbial co-cultures and can be expanded from naringenin to any metabolite if a suitable riboswitch is identified. The co-culture technique presented here can be applied to a variety of target metabolites in combination with the SELEX approach for aptamer design. Due to the compartmentalization of the two genetic constructs responsible for production and detection into separate cells and application as independent modules of a synthetic microbial co-culture we have subsequently reduced the need for re-optimization of the producer module when the biosensor is replaced or removed. Biotechnol. Bioeng. 2017;114: 2235-2244. © 2017 Wiley Periodicals, Inc.

Salt-mediated two-site ligand binding by the cocaine-binding aptamer

Multisite ligand binding by proteins is commonly utilized in the regulation of biological systems and exploited in a range of biochemical technologies. Aptamers, although widely utilized in many rationally designed biochemical systems, are rarely capable of multisite ligand binding. The cocaine-binding aptamer is often used for studying and developing sensor and aptamer-based technologies. Here, we use isothermal titration calorimetry (ITC) and NMR spectroscopy to demonstrate that the cocaine-binding aptamer switches from one-site to two-site ligand binding, dependent on NaCl concentration. The high-affinity site functions at all buffer conditions studied, the low-affinity site only at low NaCl concentrations. ITC experiments show the two ligand-binding sites operate independently of one another with different affinities and enthalpies. NMR spectroscopy shows the second binding site is located in stem 2 near the three-way junction. This ability to control ligand binding at the second site by adjusting the concentration of NaCl is rare among aptamers and may prove a useful in biotechnology applications. This work also demonstrates that in vitro selected biomolecules can have functions as complex as those found in nature.

Cobalamin riboswitches exhibit a broad range of ability to discriminate between methylcobalamin and adenosylcobalamin

Riboswitches are a widely distributed class of regulatory RNAs in bacteria that modulate gene expression via small-molecule-induced conformational changes. Generally, these RNA elements are grouped into classes based upon conserved primary and secondary structure and their cognate effector molecule. Although this approach has been very successful in identifying new riboswitch families and defining their distributions, small sequence differences between structurally related RNAs can alter their ligand selectivity and regulatory behavior. Herein, we use a structure-based mutagenic approach to demonstrate that cobalamin riboswitches have a broad spectrum of preference for the two biological forms of cobalamin in vitro using isothermal titration calorimetry. This selectivity is primarily mediated by the interaction between a peripheral element of the RNA that forms a T-loop module and a subset of nucleotides in the cobalamin-binding pocket. Cell-based fluorescence reporter assays in Escherichia coli revealed that mutations that switch effector preference in vitro lead to differential regulatory responses in a biological context. These data demonstrate that a more comprehensive analysis of representative sequences of both previously and newly discovered classes of riboswitches might reveal subgroups of RNAs that respond to different effectors. Furthermore, this study demonstrates a second distinct means by which tertiary structural interactions in cobalamin riboswitches dictate ligand selectivity.

Long-Range Interactions in Riboswitch Control of Gene Expression

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

Characterization of Engineered PreQ1 Riboswitches for Inducible Gene Regulation in Mycobacteria

We report here the behavior of naturally occurring and rationally engineered preQ1 riboswitches and their application to inducible gene regulation in mycobacteria. Because mycobacteria lack preQ1 biosynthetic genes, we hypothesized that preQ1 could be used as an exogenous nonmetabolite ligand to control riboswitches in mycobacteria. Selected naturally occurring preQ1 riboswitches were assayed and successfully drove preQ1-dependent repression of a green fluorescent protein reporter in Mycobacterium smegmatis Using structure-based design, we engineered three preQ1 riboswitches from Thermoanaerobacter tencongensis, Bacillus subtilis, and Lactobacillus rhamnosus toward achieving higher response ratios and increased repression. Assuming a steady-state model, variants of the T. tencongensis riboswitch most closely followed the predicted trends. Unexpectedly, the preQ1 dose response was best described by a model with a second, independent preQ1 binding site. This behavior was general to the preQ1 riboswitch family, since the wild type and rationally designed mutants of riboswitches from all three bacteria behaved analogously. Across all variants, the response ratios, which describe expression in the absence versus the presence of preQ1, ranged from <2 to ∼10, but repression in all cases was incomplete up to 1 mM preQ1. By reducing the transcript expression level, we obtained a preQ1 riboswitch variant appropriate for inducible knockdown applications. We further showed that the preQ1 response is reversible, is titratable, and can be used to control protein expression in mycobacteria within infected macrophages. By engineering naturally occurring preQ1 riboswitches, we have not only extended the tools available for inducible gene regulation in mycobacteria but also uncovered new behavior of these riboswitches.IMPORTANCE Riboswitches are elements found in noncoding regions of mRNA that regulate gene expression, typically in response to an endogenous metabolite. Riboswitches have emerged as important tools for inducible gene expression in diverse organisms. We noted that mycobacteria lack the biosynthesis genes for preQ1, a ligand for riboswitches from diverse bacteria. Predicting that preQ1 is not present in mycobacteria, we showed that it controls optimized riboswitches appropriate for gene knockdown applications. Further, the riboswitch response is subject to a second independent preQ1 binding event that has not been previously documented. By engineering naturally occurring riboswitches, we have uncovered a new behavior, with implications for riboswitch function in its native context, and extended the tools available for inducible gene regulation in mycobacteria.

A human microRNA precursor binding to folic acid discovered by small RNA transcriptomic SELEX

RNA aptamers are structured motifs that bind to specific molecules. A growing number of RNAs bearing aptamer elements, whose functions are modulated by direct binding of metabolites, have been found in living cells. Recent studies have suggested that more small RNAs binding to metabolites likely exist and may be involved in diverse cellular processes. However, conventional methods are not necessarily suitable for the discovery of such RNA aptamer elements in small RNAs with lengths ranging from 50 to 200 nucleotides, due to the far more abundant tRNAs in this size range. Here, we describe a new in vitro selection method to uncover naturally occurring small RNAs capable of binding to a ligand of interest, referred to as small RNA transcriptomic SELEX (smaRt-SELEX). By means of this method, we identified a motif in human precursor microRNA 125a (hsa-pre-miR-125a) that interacts with folic acid. Mutation studies revealed that the terminal loop region of hsa-pre-miR-125a is important for this binding interaction. This method has potential for the discovery of new RNA aptamer elements or catalytic motifs in biological small RNA fractions.