Computational analysis of riboswitch-based regulation

Advances in computational analysis of riboswitches in the last decade have contributed greatly to our understanding of riboswitch regulatory roles and mechanisms. Riboswitches were originally discovered as part of the sequence analysis of the 5′-untranslated region of mRNAs in the hope of finding novel gene regulatory sites, and the existence of structural RNAs appeared to be a spurious phenomenon. As more riboswitches were discovered, they illustrated the diversity and adaptability of these RNA regulatory sequences. The fact that a chemically monotonous molecule like RNA can discern a wide range of substrates and exert a variety of regulatory mechanisms was subsequently demonstrated in diverse genomes and has hastened the development of sophisticated algorithms for their analysis and prediction. In this review, we focus on some of the computational tools for riboswitch detection and secondary structure prediction. The study of this simple yet efficient form of gene regulation promises to provide a more complete picture of a world that RNA once dominated and allows rational design of artificial riboswitches. This article is part of a Special Issue entitled: 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.     

The THI-box riboswitch, or how RNA binds thiamin pyrophosphate

Riboswitches are genetic control elements present mainly in the 5' untranslated regions of messenger RNAs that, upon binding of a small metabolite (like some vitamins, amino acids, and nucleobases), undergo conformational changes, affecting the expression of downstream genes. Structural studies of riboswitches are important for understanding how they recognize their ligands with high specificity and affinity. The thiamin pyrophosphate binding riboswitch (THI- box) is widely distributed in the three kingdoms of life and is involved in very distinct modes of gene regulation. Three recent THI-box structural analyses revealed how polyanionic RNA is able to bind a molecule with a negatively charged pyrophosphate group like thiamin pyrophosphate (TPP) and how it can discriminate between TPP and monophosphorylated analog molecules. These studies give insight into the genetic regulatory mechanisms in which the THI-box is involved.     

Three-way RNA junctions with remote tertiary contacts: A recurrent and highly versatile fold

Three-way junction RNAs adopt a recurrent Y shape when two of the helices form a coaxial stack and the third helix establishes one or more tertiary contacts several base pairs away from the junction. In this review, the structure, distribution, and functional relevance of these motifs are examined. Structurally, the folds exhibit conserved junction topologies, and the distal tertiary interactions play a crucial role in determining the final shape of the structures. The junctions and remote tertiary contacts behave as flexible hinge motifs that respond to changes in the other region, providing these folds with switching mechanisms that have been shown to be functionally useful in a variety of contexts. In addition, the juxtaposition of RNA domains at the junction and at the distal tertiary complexes enables the RNA helices to adopt unusual conformations that are frequently used by proteins, RNA molecules, and antibiotics as platforms for specific binding. As a consequence of these properties, Y-shaped junctions are widely distributed in all kingdoms of life, having been observed in small naked RNAs such as riboswitches and ribozymes or embedded in complex ribonucleoprotein systems like ribosomal RNAs, RNase P, or the signal recognition particle. In all cases, the folds were found to play an essential role for the functioning or assembly of the RNA or ribonucleoprotein systems that contain them.     

Discovery and characterization of a fourth class of guanidine riboswitches

Riboswitches are RNAs that specifically sense a small molecule and regulate genes accordingly. The recent discovery of guanidine-binding riboswitches revealed the biological significance of this compound, and uncovered genes related to its biology. For example, certain sugE genes encode guanidine exporters and are activated by the riboswitches to reduce toxic levels of guanidine in the cell. In order to study guanidine biology and riboswitches, we applied a bioinformatics strategy for discovering additional guanidine riboswitches by searching for new candidate motifs associated with sugE genes. Based on in vitro and in vivo experiments, we determined that one of our six best candidates is a new structural class of guanidine riboswitches. The expression of a genetic reporter was induced 80-fold in response to addition of 5 mM guanidine in Staphylococcus aureus. This new class, called the guanidine-IV riboswitch, reveals additional guanidine-associated protein domains that are extremely rarely or never associated with previously established guanidine riboswitches. Among these protein domains are two transporter families that are structurally distinct from SugE, and could represent novel types of guanidine exporters. These results establish a new metabolite-binding RNA, further validate a bioinformatics method for finding riboswitches and suggest substrate specificities for as-yet uncharacterized transporter proteins.     

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.     

Ribozymes, riboswitches and beyond: regulation of gene expression without proteins

Although various functions of RNA are carried out in conjunction with proteins, some catalytic RNAs, or ribozymes, which contribute to a range of cellular processes, require little or no assistance from proteins. Furthermore, the discovery of metabolite-sensing riboswitches and other types of RNA sensors has revealed RNA-based mechanisms that cells use to regulate gene expression in response to internal and external changes. Structural studies have shown how these RNAs can carry out a range of functions. In addition, the contribution of ribozymes and riboswitches to gene expression is being revealed as far more widespread than was previously appreciated. These findings have implications for understanding how cellular functions might have evolved from RNA-based origins.     

The aptamer core of SAM-IV riboswitches mimics the ligand-binding site of SAM-I riboswitches

A novel family of riboswitches, called SAM-IV, is the fourth distinct set of mRNA elements to be reported that regulate gene expression via direct sensing of S-adenosylmethionine (SAM or AdoMet). SAM-IV riboswitches share conserved nucleotide positions with the previously described SAM-I riboswitches, despite rearranged structures and nucleotide positions with family-specific nucleotide identities. Sequence analysis and molecular recognition experiments suggest that SAM-I and SAM-IV riboswitches share similar ligand binding sites, but have different scaffolds. Our findings support the view that RNA has considerable structural versatility and reveal that riboswitches exploit this potential to expand the scope of RNA in genetic regulation.     

Identification and characterization of small RNAs involved in RNA silencing

Double-stranded RNA (dsRNA) is a potent trigger of sequence-specific gene silencing mechanisms known as RNA silencing or RNA interference. The recognition of the target sequences is mediated by ribonucleoprotein complexes that contain 21- to 28-nucleotide (nt) guide RNAs derived from processing of the trigger dsRNA. Here, we review the experimental and bioinformatic approaches that were used to identify and characterize these small RNAs isolated from cells and tissues. The identification and characterization of small RNAs and their expression patterns is important for elucidating gene regulatory networks.     

Chemical and functional diversity of small molecule ligands for RNA

Functional RNAs such as ribosomal RNA and structured domains of mRNA are targets for small molecule ligands that can act as modulators of the RNA biological activity. Natural ligands for RNA display a bewildering structural and chemical complexity that has yet to be matched by synthetic RNA binders. Comparison of natural and artificial ligands for RNA may help to direct future approaches to design and synthesize potent novel scaffolds for specific recognition of RNA targets.     

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.     

Ancestral Roles of Small RNAs: An Ago-Centric Perspective

RNAi has existed at least since the divergence of prokaryotes and eukaryotes. This collection of pathways responds to a diversity of “abberant” RNAs and generally silences or eliminates genes sharing sequence content with the silencing trigger. In the canonical pathway, double-stranded RNAs are processed into small RNAs, which guide effector complexes to their targets by complementary base pairing. Many alternative routes from silencing trigger to small RNA are continuously being uncovered. Though the triggers of the pathway and the mechanisms of small RNA production are many, all RNAi-related mechanisms share Argonaute proteins as the heart of their effector complexes. These can act as self-contained silencing machines, binding directly to small RNAs, carrying out homology-based target recognition, and in some cases cleaving targets using an endogenous nuclease domain. Here, we discuss the diversity of Argonaute proteins from a structural and functional perspective.     

Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy

Riboswitches are a newly discovered large family of structured functional RNA elements that specifically bind small molecule targets out of a myriad of cellular metabolites to modulate gene expression. Structural studies of ligand-bound riboswitches by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy have provided insights into detailed RNA-ligand recognition and interactions. However, the structures of ligand-free riboswitches remain poorly characterized. In this study, we have used a variety of biochemical, biophysical and computational techniques including small-angle X-ray scattering and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of SAM-II riboswitch. Our data demonstrate that the RNA adopts multiple conformations along its folding pathway and suggest that the RNA undergoes marked conformational changes upon Mg(2+) compaction and S-adenosylmethionine (SAM) metabolite binding. Further studies indicated that Mg(2+) ion is not essential for the ligand binding but can stabilize the complex by facilitating loop/stem interactions. In the presence of millimolar concentration of Mg(2+) ion, the RNA samples a more compact conformation. This conformation is near to, but distinct from, the native fold and competent to bind the metabolite. We conclude that the formation of various secondary and tertiary structural elements, including a pseudoknot, occur to sequester the putative Shine-Dalgarno sequence of the RNA only after metabolite binding.     

An mRNA structure that controls gene expression by binding S-adenosylmethionine

Riboswitches are metabolite-binding RNA structures that serve as genetic control elements for certain messenger RNAs. These RNA switches have been identified in all three kingdoms of life and are typically responsible for the control of genes whose protein products are involved in the biosynthesis, transport or utilization of the target metabolite. Herein, we report that a highly conserved RNA domain found in bacteria serves as a riboswitch that responds to the coenzyme S-adenosylmethionine (SAM) with remarkably high affinity and specificity. SAM riboswitches undergo structural reorganization upon introduction of SAM, and these allosteric changes regulate the expression of 26 genes in Bacillus subtilis. This and related findings indicate that direct interaction between small metabolites and allosteric mRNAs is an important and widespread form of genetic regulation in bacteria.     

Structural RNA has lower folding energy than random RNA of the same dinucleotide frequency

We present results of computer experiments that indicate that several RNAs for which the native state (minimum free energy secondary structure) is functionally important (type III hammerhead ribozymes, signal recognition particle RNAs, U2 small nucleolar spliceosomal RNAs, certain riboswitches, etc.) all have lower folding energy than random RNAs of the same length and dinucleotide frequency. Additionally, we find that whole mRNA as well as 5'-UTR, 3'-UTR, and cds regions of mRNA have folding energies comparable to that of random RNA, although there may be a statistically insignificant trace signal in 3'-UTR and cds regions. Various authors have used nucleotide (approximate) pattern matching and the computation of minimum free energy as filters to detect potential RNAs in ESTs and genomes. We introduce a new concept of the asymptotic Z-score and describe a fast, whole-genome scanning algorithm to compute asymptotic minimum free energy Z-scores of moving-window contents. Asymptotic Z-score computations offer another filter, to be used along with nucleotide pattern matching and minimum free energy computations, to detect potential functional RNAs in ESTs and genomic regions.