Themes and variations in riboswitch structure and function

The complexity of gene expression control by non-coding RNA has been highlighted by the recent progress in the field of riboswitches. Discovered a decade ago, riboswitches represent a diverse group of non-coding mRNA regions that possess a unique ability to directly sense cellular metabolites and modulate gene expression through formation of alternative metabolite-free and metabolite-bound conformations. Such protein-free metabolite sensing domains utilize sophisticated three-dimensional folding of RNA molecules to discriminate between a cognate ligand from related compounds so that only the right ligand would trigger a genetic response. Given the variety of riboswitch ligands ranging from small cations to large coenzymes, riboswitches adopt a great diversity of structures. Although many riboswitches share structural principles to build metabolite-competent folds, form precise ligand-binding pockets, and communicate a ligand-binding event to downstream regulatory regions, virtually all riboswitch classes possess unique features for ligand recognition, even those tuned to recognize the same metabolites. Here we present an overview of the biochemical and structural research on riboswitches with a major focus on common principles and individual characteristics adopted by these regulatory RNA elements during evolution to specifically target small molecules and exert genetic responses. This article is part of a Special Issue entitled: Riboswitches.     

The Ligand-Free State of the TPP Riboswitch: A Partially Folded RNA Structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.     

Common themes and differences in SAM recognition among SAM riboswitches

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

Prospects for riboswitch discovery and analysis

An expanding number of metabolite-binding riboswitch classes are being discovered in the noncoding portions of bacterial genomes. Findings over the last decade indicate that bacteria commonly use these RNA genetic elements as regulators of metabolic pathways and as mediators of changes in cell physiology. Some riboswitches are surprisingly complex, and they rival protein factors in their structural and functional sophistication. Each new riboswitch discovery expands our knowledge of the biochemical capabilities of RNA, and some give rise to new questions that require additional study to be addressed. Some of the greatest prospects for riboswitch research and some of the more interesting mysteries are discussed in this review.     

Amino acid recognition and gene regulation by riboswitches

Riboswitches specifically control expression of genes predominantly involved in biosynthesis, catabolism and transport of various cellular metabolites in organisms from all three kingdoms of life. Among many classes of identified riboswitches, two riboswitches respond to amino acids lysine and glycine to date. Though these riboswitches recognize small compounds, they both belong to the largest riboswitches and have unique structural and functional characteristics. In this review, we attempt to characterize molecular recognition principles employed by amino acid-responsive riboswitches to selectively bind their cognate ligands and to effectively perform a gene regulation function. We summarize up-to-date biochemical and genetic data available for the lysine and glycine riboswitches and correlate these results with recent high-resolution structural information obtained for the lysine riboswitch. We also discuss the contribution of lysine riboswitches to antibiotic resistance and outline potential applications of riboswitches in biotechnology and medicine.     

Evolutionary Evidence for Alternative Structure in RNA Sequence Co-variation

Sequence conservation and co-variation of base pairs are hallmarks of structured RNAs. For certain RNAs (e.g. riboswitches), a single sequence must adopt at least two alternative secondary structures to effectively regulate the message. If alternative secondary structures are important to the function of an RNA, we expect to observe evolutionary co-variation supporting multiple conformations. We set out to characterize the evolutionary co-variation supporting alternative conformations in riboswitches to determine the extent to which alternative secondary structures are conserved. We found strong co-variation support for the terminator, P1, and anti-terminator stems in the purine riboswitch by extending alignments to include terminator sequences. When we performed Boltzmann suboptimal sampling on purine riboswitch sequences with terminators we found that these sequences appear to have evolved to favor specific alternative conformations. We extended our analysis of co-variation to classic alignments of group I/II introns, tRNA, and other classes of riboswitches. In a majority of these RNAs, we found evolutionary evidence for alternative conformations that are compatible with the Boltzmann suboptimal ensemble. Our analyses suggest that alternative conformations are selected for and thus likely play functional roles in even the most structured of RNAs.     

Coenzyme B12 riboswitches are widespread genetic control elements in prokaryotes

Recent studies have begun to reveal that numerous fundamental metabolic pathways in bacteria are regulated by riboswitches residing within certain messenger RNAs. These riboswitches selectively bind metabolites and modulate gene expression in response to changing ligand concentrations. Previously, we provided evidence that the btuB mRNAs of Escherichia coli and Salmonella typhimurium each carry a coenzyme B12-dependent riboswitch that causes repressed translation of the encoded cobalamin-transport protein at elevated coenzyme concentrations. Herein, we use a phylogenetic analysis to define a consensus sequence and secondary structure model for the ligand- binding domain of this riboswitch class. RNA structures that conform to this model are widespread in both Gram-positive and Gram-negative organisms. In addition, we find that the 5'-untranslated region (5'-UTR) of the cobalamin biosynthesis (cob) operon of S.typhimurium carries an RNA motif that matches this consensus sequence. Biochemical and genetic characterization of this motif confirms that the RNA directly binds coenzyme B12, and that it likely serves as a genetic control element for regulating expression of the 25-gene operon for cobalamin production in this pathogen.     

Riboswitches that sense S-adenosylmethionine and S-adenosylhomocysteine.

Numerous riboswitches have been discovered that specifically recognize metabolites and modulate gene expression. Each riboswitch class is defined either by the consensus sequence and structural features of its metabolite-binding aptamer domain, or by the distinct metabolite that the aptamer recognizes. Several distinct classes of riboswitches that respond to S-adenosylmethionine (SAM or AdoMet) have been discovered. Representatives of these classes have been shown to strongly discriminate against S-adenosylhomocystenine (SAH or AdoHcy), which is the metabolic byproduct produced when SAM is used as a cofactor for methylation reactions. However, a distinct class of riboswitches that selectively binds SAH, and strongly discriminates against SAM, also has been discovered. Herein we compare the features of SAM and SAH riboswitches, which help showcase the enormous structural diversity that RNA can harness to form precision genetic switches for compounds that are critical for fundamental metabolic processes.     

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.     

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or “azaaromatics,” trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.     

Riboswitch regulation mechanisms: RNA,metabolites and regulatory proteins

Riboswitches are RNA sensors that have been shown to modulate the expression of downstream genes by altering their structure upon metabolite binding. Riboswitches are unique among cellular regulators in that metabolite detection is strictly performed using RNA interactions with the sensed metabolite and in which no regulatory protein is needed to mediate the interaction. However, recent studies have shed light on riboswitch control mechanisms relying on protein regulators to harness metabolite binding for the mediation of gene expression, thereby increasing the range of cellular factors involved in riboswitch regulation. The interaction between riboswitches and proteins adds another level of evolutionary pressure as riboswitches must maintain key residues for metabolite detection, structural switching and protein binding sites. Here, we review regulatory mechanisms involving Escherichia coli riboswitches that have recently been shown to rely on regulatory proteins. We also discuss the implication of such protein-based riboswitch regulatory mechanisms for genetic regulation.     

Structural principles of nucleoside selectivity in a 2'-deoxyguanosine riboswitch

Purine riboswitches have an essential role in genetic regulation of bacterial metabolism. This family includes the 2'-deoxyguanosine (dG) riboswitch, which is involved in feedback control of deoxyguanosine biosynthesis. To understand the principles that define dG selectivity, we determined crystal structures of the natural Mesoplasma florum riboswitch bound to cognate dG as well as to noncognate guanosine, deoxyguanosine monophosphate and guanosine monophosphate. Comparison with related purine riboswitch structures reveals that the dG riboswitch achieves its specificity through modification of key interactions involving the nucleobase and rearrangement of the ligand-binding pocket to accommodate the additional sugar moiety. In addition, we observe new conformational changes beyond the junctional binding pocket extending as far as peripheral loop-loop interactions. It appears that re-engineering riboswitch scaffolds will require consideration of selectivity features dispersed throughout the riboswitch tertiary fold, and structure-guided drug design efforts targeted to junctional RNA scaffolds need to be addressed within such an expanded framework.     

Design Principles for Riboswitch Function

Scientific and technological advances that enable the tuning of integrated regulatory components to match network and system requirements are critical to reliably control the function of biological systems. RNA provides a promising building block for the construction of tunable regulatory components based on its rich regulatory capacity and our current understanding of the sequence–function relationship. One prominent example of RNA-based regulatory components is riboswitches, genetic elements that mediate ligand control of gene expression through diverse regulatory mechanisms. While characterization of natural and synthetic riboswitches has revealed that riboswitch function can be modulated through sequence alteration, no quantitative frameworks exist to investigate or guide riboswitch tuning. Here, we combined mathematical modeling and experimental approaches to investigate the relationship between riboswitch function and performance. Model results demonstrated that the competition between reversible and irreversible rate constants dictates performance for different regulatory mechanisms. We also found that practical system restrictions, such as an upper limit on ligand concentration, can significantly alter the requirements for riboswitch performance, necessitating alternative tuning strategies. Previous experimental data for natural and synthetic riboswitches as well as experiments conducted in this work support model predictions. From our results, we developed a set of general design principles for synthetic riboswitches. Our results also provide a foundation from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands.