Riboswitch Detection Using Profile Hidden Markov Models

Background  

Riboswitches are a type of noncoding RNA that regulate gene expression by switching from one structural conformation to another on ligand binding. The various classes of riboswitches discovered so far are differentiated by the ligand, which on binding induces a conformational switch. Every class of riboswitch is characterized by an aptamer domain, which provides the site for ligand binding, and an expression platform that undergoes conformational change on ligand binding. The sequence and structure of the aptamer domain is highly conserved in riboswitches belonging to the same class. We propose a method for fast and accurate identification of riboswitches using profile Hidden Markov Models (pHMM). Our method exploits the high degree of sequence conservation that characterizes the aptamer domain.     

Conformational changes in the expression domain of the Escherichia coli thiM riboswitch

The thiM riboswitch contains an aptamer domain that adaptively binds the coenzyme thiamine pyrophosphate (TPP). The binding of TPP to the aptamer domain induces structural rearrangements that are relayed to a second domain, the so-called expression domain, thereby interfering with gene expression. The recently solved crystal structures of the aptamer domains of the thiM riboswitches in complex with TPP revealed how TPP stabilizes secondary and tertiary structures in the RNA ligand complex. To understand the global modes of reorganization between the two domains upon metabolite binding the structure of the entire riboswitch in presence and absence of TPP needs to be determined. Here we report the secondary structure of the entire thiM riboswitch from Escherichia coli in its TPP-free form and its transition into the TPP-bound variant, thereby depicting domains of the riboswitch that serve as communication links between the aptamer and the expression domain. Furthermore, structural probing provides an explanation for the lack of genetic control exerted by a riboswitch variant with mutations in the expression domain that still binds TPP.     

Secondary structures and functional requirements for thiM riboswitches from Desulfovibrio vulgaris, Erwinia carotovora and Rhodobacter spheroides

Abstract Bacterial thiM riboswitches contain aptamer domains that bind the metabolite thiamine pyrophosphate (TPP). Binding of TPP to the aptamer domain induces structural rearrangements that are relayed to the expression domain, thereby interfering with gene expression. Here, we report identification of three putative thiM riboswitches from different bacteria and analysis of their secondary structures. Chemical probing revealed that the riboswitches share similar secondary structures in their aptamer domains that can communicate with the highly variant expression domains in a mechanism likely involving sequestration of the Shine-Dalgarno sequence. Remarkably, the aptamer domain of the thiM gene of Desulfovibrio vulgaris binds TPP with similar affinity and selectivity as that of Escherichia coli, although nucleotides previously shown to form direct contacts to the metabolite are mutated. We also designed small RNA hairpins for each riboswitch that bind the RNA only in the absence of the metabolite. Our study shows that aptamer domains in riboswitches with high similarity in their secondary structures can communicate with a broad variety of non-related expression domains by similar mechanisms.     

Rational design of SAM analogues targeting SAM-II riboswitch aptamer

Riboswitches are noncoding RNA elements embedded in 5′-untranslated region of many bacterial mRNAs regulating gene expression in response to essential metabolites. They are unique from other RNA targets because they have evolved to form specific structural receptors for the purpose of binding small molecular metabolites suggesting that structure-based rational drug design approach may be used in designing metabolite mimics targeting riboswitches. We have developed a fluorescence binding assay for SAM-II riboswitch aptamer and identified an S-adenosylmethionine (SAM) analogue that selectively binds to SAM-II riboswitch aptamer with comparable binding affinity to its native metabolite using structure-based design approach.     

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.     

Ligand binding and gene control characteristics of tandem riboswitches in Bacillus anthracis

Most riboswitches are composed of a single metabolite-binding aptamer and a single expression platform that function together to regulate genes in response to changing metabolite concentrations. In rare instances, two aptamers or sometimes two complete riboswitches reside adjacent to each other in untranslated regions (UTRs) of mRNAs. We have examined an example of a tandem riboswitch in the Gram-positive bacterium Bacillus anthracis that includes two complete riboswitches for thiamine pyrophosphate (TPP). Unlike other complex riboswitch systems described recently, tandem TPP riboswitches do not exhibit cooperative ligand binding and do not detect two different types of metabolites. In contrast, both riboswitches respond independently to TPP and are predicted to function in concert to mimic the more "digital" gene control outcome observed when two aptamers bind ligands cooperatively. Our findings further demonstrate that simple gene control elements made only of RNA can be assembled in different architectures to yield more complex gene control outcomes.     

Rational design of artificial riboswitches based on ligand-dependent modulation of internal ribosome entry in wheat germ extract and their applications as label-free biosensors

Riboswitches are RNA elements in mRNA that control gene expression in cis in response to their specific ligands. Because artificial riboswitches make it possible to regulate any gene with an arbitrary molecule, they are expected to function as biosensors, in which the output is easily detectable protein expression. I report herein a fully rational design strategy for artificially constructing novel riboswitches that work in a eukaryotic cell-free translation system (wheat germ extract). In these riboswitches, translation mediated by an internal ribosome entry site (IRES) is promoted only in the presence of a specific ligand (ON), while it is inhibited in the absence of the ligand (OFF). The first rationally designed riboswitch, which is regulated by theophylline, showed a high switching efficiency and dependency on theophylline. In addition, based on the design of the theophylline-dependent riboswitch, other three kinds of riboswitches controlled by FMN, tetracycline, and sulforhodamine B, were constructed only by calculating the ΔG value of one stem-loop structure. The rational design strategy described herein is therefore useful for easily producing various ligand-dependent riboswitches, which are available as biosensors for detecting their ligands.     

A structural intermediate pre-organizes the add adenine riboswitch for ligand recognition

Riboswitches are RNA sequences that regulate gene expression by undergoing structural changes upon the specific binding of cellular metabolites. Crystal structures of purine-sensing riboswitches have revealed an intricate network of interactions surrounding the ligand in the bound complex. The mechanistic details about how the aptamer folding pathway is involved in the formation of the metabolite binding site have been previously shown to be highly important for the riboswitch regulatory activity. Here, a combination of single-molecule FRET and SHAPE assays have been used to characterize the folding pathway of the adenine riboswitch from Vibrio vulnificus. Experimental evidences suggest a folding process characterized by the presence of a structural intermediate involved in ligand recognition. This intermediate state acts as an open conformation to ensure ligand accessibility to the aptamer and folds into a structure nearly identical to the ligand-bound complex through a series of structural changes. This study demonstrates that the add riboswitch relies on the folding of a structural intermediate that pre-organizes the aptamer global structure and the ligand binding site to allow efficient metabolite sensing and riboswitch genetic regulation.     

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.     

Riboswitch structure: an internal residue mimicking the purine ligand

The adenine and guanine riboswitches regulate gene expression in response to their purine ligand. X-ray structures of the aptamer moiety of these riboswitches are characterized by a compact fold in which the ligand forms a Watson–Crick base pair with residue 65. Phylogenetic analyses revealed a strict restriction at position 39 of the aptamer that prevents the G39–C65 and A39–U65 combinations, and mutational studies indicate that aptamers with these sequence combinations are impaired for ligand binding. In order to investigate the rationale for sequence conservation at residue 39, structural characterization of the U65C mutant from Bacillus subtilis pbuE adenine riboswitch aptamer was undertaken. NMR spectroscopy and X-ray crystallography studies demonstrate that the U65C mutant adopts a compact ligand-free structure, in which G39 occupies the ligand-binding site of purine riboswitch aptamers. These studies present a remarkable example of a mutant RNA aptamer that adopts a native-like fold by means of ligand mimicking and explain why this mutant is impaired for ligand binding. Furthermore, this work provides a specific insight into how the natural sequence has evolved through selection of nucleotide identities that contribute to formation of the ligand-bound state, but ensures that the ligand-free state remains in an active conformation.     

Folding of the SAM-I riboswitch: A tale with a twist

Riboswitches are ligand-dependent RNA genetic regulators that control gene expression by altering their structures. The elucidation of riboswitch conformational changes before and after ligand recognition is crucial to understand how riboswitches can achieve high ligand binding affinity and discrimination against cellular analogs. The detailed characterization of riboswitch folding pathways suggest that they may use their intrinsic conformational dynamics to sample a large array of structures, some of which being nearly identical to ligand-bound molecules. Some of these structural conformers can be "captured" upon ligand binding, which is crucial for the outcome of gene regulation. Recent studies about the SAM-I riboswitch have revealed unexpected and previously unknown RNA folding mechanisms. For instance, the observed helical twist of the P1 stem upon ligand binding to the SAM-I aptamer adds a new element in the repertoire of RNA strategies for recognition of small metabolites. From an RNA folding perspective, these findings also strongly indicate that the SAM-I riboswitch could achieve ligand recognition by using an optimized combination of conformational capture and induced-fit approaches, a feature that may be shared by other RNA regulatory sequences.     

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, B₁₂-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.     

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.     

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.