Flipping the script: Understanding riboswitches from an alternative perspective.

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

The RNA-Binding Domain of hnRNP U Extends beyond the RGG/RG Motifs.

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG motif within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV cross-linking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing 11 of these identified cross-link sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG motif is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG motif with RNA-binding properties comparable to those of the canonical RGG/RG motif. These RGG/RG motifs serve redundant functions, with neither serving as the primary RBD. While in isolation, each RGG/RG motif has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling the binding of most of the designed RNA set with low to midnanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlight the perils of a reductionist approach to defining biochemical activities in this system and pave the way for a detailed investigation of its RNA-binding specificity.

An Extended DNA Binding Domain of the Estrogen Receptor Alpha Directly Interacts with RNAs in Vitro.

Estrogen receptor alpha (ERα) is a ligand-responsive transcription factor critical for sex determination and development. Recent reports challenge the canonical view of ERα function by suggesting an activity beyond binding dsDNA at estrogen-responsive promotor elements: association with RNAs in vivo. Whether these interactions are direct or indirect remains unknown, which limits the ability to understand the extent, specificity, and biological role of ERα-RNA binding. Here we demonstrate that an extended DNA-binding domain of ERα directly binds a wide range of RNAs in vitro with structural specificity. ERα binds RNAs that adopt a range of hairpin-derived structures independent of sequence, while interacting poorly with single- and double-stranded RNA. RNA affinities are only 4-fold weaker than consensus dsDNA and significantly tighter than nonconsensus dsDNA sequences. Moreover, RNA binding is competitive with DNA binding. Together, these data show that ERα utilizes an extended DNA-binding domain to achieve a high-affinity/low-specificity mode for interacting with RNA.

The DNA-Binding High-Mobility Group Box Domain of Sox Family Proteins Directly Interacts with RNA In Vitro.

There is a growing body of evidence that a substantial number of protein domains identified as DNA-binding also interact with RNA to regulate biological processes. Several recent studies have revealed that the Sox2 transcription factor binds RNA through its high-mobility group box (HMGB) domain in vitro and in vivo. A high degree of conservation of this domain among members of the Sox family of transcription factors suggests that RNA-binding activity may be a general feature of these proteins. To address this hypothesis, we examined a subset of HMGB domains from human Sox family of proteins for their ability to bind both DNA and RNA in vitro. We observed selective, high-affinity interactions between Sox family HMGB domains and various model RNA elements, including a four-way junction RNA, a hairpin RNA with an internal bulge, G-quadruplex RNA, and a fragment of long noncoding RNA ES2, which is known to directly interact with Sox2. Importantly, the HMGB domains bind these RNA ligands significantly tighter than nonconsensus dsDNA and in some cases with affinities rivaling those of their consensus dsDNA sequences. These data suggest that RNA binding is a conserved feature of the Sox family of transcription factors with the potential to modulate unappreciated biological functions.

RNA-Puzzles Round IV: 3D structure predictions of four ribozymes and two aptamers.

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA 3D structure prediction. With agreement from crystallographers, the RNA structures are predicted by various groups before the publication of the crystal structures. We now report the prediction of 3D structures for six RNA sequences: four nucleolytic ribozymes and two riboswitches. Systematic protocols for comparing models and crystal structures are described and analyzed. In these six puzzles, we discuss (i) the comparison between the automated web servers and human experts; (ii) the prediction of coaxial stacking; (iii) the prediction of structural details and ligand binding; (iv) the development of novel prediction methods; and (v) the potential improvements to be made. We show that correct prediction of coaxial stacking and tertiary contacts is essential for the prediction of RNA architecture, while ligand binding modes can only be predicted with low resolution and simultaneous prediction of RNA structure with accurate ligand binding still remains out of reach. All the predicted models are available for the future development of force field parameters and the improvement of comparison and assessment tools.

SPRINT: a Cas13a-based platform for detection of small molecules.

Recent efforts in biological engineering have made detection of nucleic acids in samples more rapid, inexpensive and sensitive using CRISPR-based approaches. We expand one of these Cas13a-based methods to detect small molecules in a one-batch assay. Using SHERLOCK-based profiling of in vitrotranscription (SPRINT), in vitro transcribed RNA sequence-specifically triggers the RNase activity of Cas13a. This event activates its non-specific RNase activity, which enables cleavage of an RNA oligonucleotide labeled with a quencher/fluorophore pair and thereby de-quenches the fluorophore. This fluorogenic output can be measured to assess transcriptional output. The use of riboswitches or proteins to regulate transcription via specific effector molecules is leveraged as a coupled assay that transforms effector concentration into fluorescence intensity. In this way, we quantified eight different compounds, including cofactors, nucleotides, metabolites of amino acids, tetracycline and monatomic ions in samples. In this manner, hundreds of reactions can be easily quantified in a few hours. This increased throughput also enables detailed characterization of transcriptional regulators, synthetic compounds that inhibit transcription, or other coupled enzymatic reactions. These SPRINT reactions are easily adaptable to portable formats and could therefore be used for the detection of analytes in the field or at point-of-care situations.

hnRNPK recognition of the B motif of Xist and other biological RNAs.

Heterogeneous nuclear ribonuclear protein K (hnRNPK) is an abundant RNA-binding protein crucial for a wide variety of biological processes. While its binding preference for multi-cytosine-patch (C-patch) containing RNA is well documented, examination of binding to known cellular targets that contain C-patches reveals an unexpected breadth of binding affinities. Analysis of in-cell crosslinking data reinforces the notion that simple C-patch preference is not fully predictive of hnRNPK localization within transcripts. The individual RNA-binding domains of hnRNPK work together to interact with RNA tightly, with the KH3 domain being neither necessary nor sufficient for binding. Rather, the RG/RGG domain is implicated in providing essential contributions to RNA-binding, but not DNA-binding, affinity. hnRNPK is essential for X chromosome inactivation, where it interacts with Xist RNA specifically through the Xist B-repeat region. We use this interaction with an RNA motif derived from this B-repeat region to determine the RNA-structure dependence of C-patch recognition. While the location preferences of hnRNPK for C-patches are conformationally restricted within the hairpin, these structural constraints are relieved in the absence of RNA secondary structure. Together, these results illustrate how this multi-domain protein's ability to accommodate and yet discriminate between diverse cellular RNAs allows for its broad cellular functions.

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 glucocorticoid receptor DNA-binding domain recognizes RNA hairpin structures with high affinity.

The glucocorticoid receptor (GR) binds the noncoding RNA Gas5 via its DNA-binding domain (DBD) with functional implications in pro-apoptosis signaling. Here, we report a comprehensive in vitro binding study where we have determined that GR-DBD is a robust structure-specific RNA-binding domain. GR-DBD binds to a diverse range of RNA hairpin motifs, both synthetic and biologically derived, with apparent mid-nanomolar affinity while discriminating against uniform dsRNA. As opposed to dimeric recognition of dsDNA, GR-DBD binds to RNA as a monomer and confers high affinity primarily through electrostatic contacts. GR-DBD adopts a discrete RNA-bound state, as assessed by NMR, distinct from both free and DNA-bound. NMR and alanine mutagenesis suggest a heightened involvement of the C-terminal α-helix of the GR-DBD in RNA-binding. RNA competes for binding with dsDNA and occurs in a similar affinity range as dimer binding to the canonical DNA element. Given the prevalence of RNA hairpins within the transcriptome, our findings strongly suggest that many RNAs have potential to impact GR biology.

A multicolor riboswitch-based platform for imaging of RNA in live mammalian cells.

RNAs directly regulate a vast array of cellular processes, emphasizing the need for robust approaches to fluorescently label and track RNAs in living cells. Here, we develop an RNA imaging platform using the cobalamin riboswitch as an RNA tag and a series of probes containing cobalamin as a fluorescence quencher. This highly modular 'Riboglow' platform leverages different colored fluorescent dyes, linkers and riboswitch RNA tags to elicit fluorescence turn-on upon binding RNA. We demonstrate the ability of two different Riboglow probes to track mRNA and small noncoding RNA in live mammalian cells. A side-by-side comparison revealed that Riboglow outperformed the dye-binding aptamer Broccoli and performed on par with the gold standard RNA imaging system, the MS2-fluorescent protein system, while featuring a much smaller RNA tag. Together, the versatility of the Riboglow platform and ability to track diverse RNAs suggest broad applicability for a variety of imaging approaches.

A functional genetic screen reveals sequence preferences within a key tertiary interaction in cobalamin riboswitches required for ligand selectivity.

Riboswitches are structured mRNA sequences that regulate gene expression by directly binding intracellular metabolites. Generating the appropriate regulatory response requires the RNA rapidly and stably acquire higher-order structure to form the binding pocket, bind the appropriate effector molecule and undergo a structural transition to inform the expression machinery. These requirements place riboswitches under strong kinetic constraints, likely restricting the sequence space accessible by recurrent structural modules such as the kink turn and the T-loop. Class-II cobalamin riboswitches contain two T-loop modules: one directing global folding of the RNA and another buttressing the ligand binding pocket. While the T-loop module directing folding is highly conserved, the T-loop associated with binding is substantially less so, with no clear consensus sequence. To further understand the functional role of the binding-associated module, a functional genetic screen of a library of riboswitches with the T-loop and its interacting nucleotides was used to build an experimental phylogeny comprised of sequences that possess a wide range of cobalamin-dependent regulatory activity. Our results reveal conservation patterns of the T-loop and its interaction with the binding core that allow for rapid tertiary structure formation and demonstrate its importance for generating strong ligand-dependent repression of mRNA expression.

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.

RNA-Puzzles Round III: 3D RNA structure prediction of five riboswitches and one ribozyme.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5'-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson-Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Recurrent RNA motifs as scaffolds for genetically encodable small-molecule biosensors.

Allosteric RNA devices are increasingly being viewed as important tools capable of monitoring enzyme evolution, optimizing engineered metabolic pathways, facilitating gene discovery and regulators of nucleic acid-based therapeutics. A key bottleneck in the development of these platforms is the availability of small-molecule-binding RNA aptamers that robustly function in the cellular environment. Although aptamers can be raised against nearly any desired target through in vitro selection, many cannot easily be integrated into devices or do not reliably function in a cellular context. Here, we describe a new approach using secondary- and tertiary-structural scaffolds derived from biologically active riboswitches and small ribozymes. When applied to the neurotransmitter precursors 5-hydroxytryptophan and 3,4-dihydroxyphenylalanine, this approach yielded easily identifiable and characterizable aptamers predisposed for coupling to readout domains to allow engineering of nucleic acid-sensory devices that function in vitro and in the cellular context.

Intrinsically disordered RGG/RG domains mediate degenerate specificity in RNA binding.

RGG/RG domains are the second most common RNA binding domain in the human genome, yet their RNA-binding properties remain poorly understood. Here, we report a detailed analysis of the RNA binding characteristics of intrinsically disordered RGG/RG domains from Fused in Sarcoma (FUS), FMRP and hnRNPU. For FUS, previous studies defined RNA binding as mediated by its well-folded domains; however, we show that RGG/RG domains are the primary mediators of binding. RGG/RG domains coupled to adjacent folded domains can achieve affinities approaching that of full-length FUS. Analysis of RGG/RG domains from FUS, FMRP and hnRNPU against a spectrum of contrasting RNAs reveals that each display degenerate binding specificity, while still displaying different degrees of preference for RNA.

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.

B12 cofactors directly stabilize an mRNA regulatory switch.

Structures of riboswitch receptor domains bound to their effector have shown how messenger RNAs recognize diverse small molecules, but mechanistic details linking the structures to the regulation of gene expression remain elusive. To address this, here we solve crystal structures of two different classes of cobalamin (vitamin B(12))-binding riboswitches that include the structural switch of the downstream regulatory domain. These classes share a common cobalamin-binding core, but use distinct peripheral extensions to recognize different B(12) derivatives. In each case, recognition is accomplished through shape complementarity between the RNA and cobalamin, with relatively few hydrogen bonding interactions that typically govern RNA-small molecule recognition. We show that a composite cobalamin-RNA scaffold stabilizes an unusual long-range intramolecular kissing-loop interaction that controls mRNA expression. This is the first, to our knowledge, riboswitch crystal structure detailing how the receptor and regulatory domains communicate in a ligand-dependent fashion to regulate mRNA expression.

Structural basis for diversity in the SAM clan of riboswitches.

In bacteria, sulfur metabolism is regulated in part by seven known families of riboswitches that bind S-adenosyl-l-methionine (SAM). Direct binding of SAM to these mRNA regulatory elements governs a downstream secondary structural switch that communicates with the transcriptional and/or translational expression machinery. The most widely distributed SAM-binding riboswitches belong to the SAM clan, comprising three families that share a common SAM-binding core but differ radically in their peripheral architecture. Although the structure of the SAM-I member of this clan has been extensively studied, how the alternative peripheral architecture of the other families supports the common SAM-binding core remains unknown. We have therefore solved the X-ray structure of a member of the SAM-I/IV family containing the alternative "PK-2" subdomain shared with the SAM-IV family. This structure reveals that this subdomain forms extensive interactions with the helix housing the SAM-binding pocket, including a highly unusual mode of helix packing in which two helices pack in a perpendicular fashion. Biochemical and genetic analysis of this RNA reveals that SAM binding induces many of these interactions, including stabilization of a pseudoknot that is part of the regulatory switch. Despite strong structural similarity between the cores of SAM-I and SAM-I/IV members, a phylogenetic analysis of sequences does not indicate that they derive from a common ancestor.

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

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

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

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