Dynamic energy landscapes of riboswitches help interpret conformational rearrangements and function

Riboswitch Distribution in the Human Gut Microbiome Reveals Common Metabolite Pathways

Riboswitches are widely distributed, conserved RNAs which regulate metabolite levels in bacterial cells through direct, noncovalent binding of their cognate metabolite. Various riboswitch families are highly enriched in gut bacteria, suggestive of a symbiotic relationship between the host and bacteria. Previous studies of the distribution of riboswitches have examined bacterial taxa broadly. Thus, the distribution of riboswitches associated with bacteria inhabiting the intestines of healthy individuals is not well understood. To address these questions, we survey the gut microbiome for riboswitches by including an international database of prokaryotic genomes from the gut samples. Using Infernal, a program that uses RNA-specific sequence and structural features, we survey this data set using existing riboswitch models. We identify 22 classes of riboswitches with vitamin cofactors making up the majority of riboswitch-associated pathways. Our finding is reproducible in other representative databases from the oral as well as the marine microbiomes, underscoring the importance of thiamine pyrophosphate, cobalamin, and flavin mononucleotide in gene regulation. Interestingly, riboswitches do not vary significantly across microbiome representatives from around the world despite major taxonomic differences; this suggests an underlying conservation. Further studies elucidating the role of bacterial riboswitches in the host metabolome are needed to illuminate the consequences of our finding.

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.

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

Out-of-Equilibrium Biophysical Chemistry: The Case for Multidimensional, Integrated Single-Molecule Approaches

Out-of-equilibrium processes are ubiquitous across living organisms and all structural hierarchies of life. At the molecular scale, out-of-equilibrium processes (for example, enzyme catalysis, gene regulation, and motor protein functions) cause biological macromolecules to sample an ensemble of conformations over a wide range of time scales. Quantifying and conceptualizing the structure–dynamics to function relationship is challenging because continuously evolving multidimensional energy landscapes are necessary to describe nonequilibrium biological processes in biological macromolecules. In this perspective, we explore the challenges associated with state-of-the-art experimental techniques to understanding biological macromolecular function. We argue that it is time to revisit how we probe and model functional out-of-equilibrium biomolecular dynamics. We suggest that developing integrated single-molecule multiparametric force–fluorescence instruments and using advanced molecular dynamics simulations to study out-of-equilibrium biomolecules will provide a path towards understanding the principles of and mechanisms behind the structure–dynamics to function paradigm in biological macromolecules.

RNA Modeling with the Computational Energy Landscape Framework

The recent advances in computational abilities, such as the enormous speed-ups provided by GPU computing, allow for large scale computational studies of RNA molecules at an atomic level of detail. As RNA molecules are known to adopt multiple conformations with comparable energies, but different two-dimensional structures, all-atom models are necessary to better describe the structural ensembles for RNA molecules. This point is important because different conformations can exhibit different functions, and their regulation or mis-regulation is linked to a number of diseases. Problematically, the energy barriers between different conformational ensembles are high, resulting in long time scales for interensemble transitions. The computational potential energy landscape framework was designed to overcome this problem of broken ergodicity by use of geometry optimization. Here, we describe the algorithms used in the energy landscape explorations with the OPTIM and PATHSAMPLE programs, and how they are used in biomolecular simulations. We present a recent case study of the 5'-hairpin of RNA 7SK to illustrate how the method can be applied to interpret experimental results, and to obtain a detailed description of molecular properties.

Beyond Plug and Pray: Context Sensitivity and in silico Design of Artificial Neomycin Riboswitches

Gene regulation in prokaryotes often depends on RNA elements such as riboswitches or RNA thermometers located in the 5' untranslated region of mRNA. Rearrangements of the RNA structure in response, e.g., to the binding of small molecules or ions control translational initiation or premature termination of transcription and thus mRNA expression. Such structural responses are amenable to computational modelling, making it possible to rationally design synthetic riboswitches for a given aptamer. Starting from an artificial aptamer, we construct the first synthetic transcriptional riboswitches that respond to the antibiotic neomycin. We show that the switching behaviour in vivo critically depends not only on the sequence of the riboswitch itself, but also on its sequence context. We therefore developed in silico methods to predict the impact of the context, making it possible to adapt the design and to rescue non-functional riboswitches. We furthermore analyse the influence of 5' hairpins with varying stability on neomycin riboswitch activity. Our data highlight the limitations of a simple plug-and-play approach in the design of complex genetic circuits and demonstrate that detailed computational models significantly simplify, improve, and automate the design of transcriptional circuits. Our design software is available under a free licence on GitHub (https://github.com/xileF1337/riboswitch_design).

Structural prediction of RNA switches using conditional base-pair probabilities

An RNA switch triggers biological functions by toggling between two conformations. RNA switches include bacterial riboswitches, where ligand binding can stabilize a bound structure. For RNAs with only one stable structure, structural prediction usually just requires a straightforward free energy minimization, but for an RNA switch, the prediction of a less stable alternative structure is often computationally costly and even problematic. The current sampling-clustering method predicts stable and alternative structures by partitioning structures sampled from the energy landscape into two clusters, but it is very time-consuming. Instead, we predict the alternative structure of an RNA switch from conditional probability calculations within the energy landscape. First, our method excludes base pairs related to the most stable structure in the energy landscape. Then, it detects stable stems (“seeds”) in the remaining landscape. Finally, it folds an alternative structure prediction around a seed. While having comparable riboswitch classification performance, the conditional-probability computations had fewer adjustable parameters, offered greater predictive flexibility, and were more than one thousand times faster than the sampling step alone in sampling-clustering predictions, the competing standard. Overall, the described approach helps traverse thermodynamically improbable energy landscapes to find biologically significant substructures and structures rapidly and effectively.

Towards an understanding of RNA structural modalities: a riboswitch case study

A riboswitch is a type of RNA molecule that regulates important biological functions by changing structure, typically under ligand-binding. We assess the extent that these ligand-bound structural alternatives are present in the Boltzmann sample, a standard RNA secondary structure prediction method, for three riboswitch test cases. We use the cluster analysis tool RNAStructProfiling to characterize the different modalities present among the suboptimal structures sampled. We compare these modalities to the putative base pairing models obtained from independent experiments using NMR or fluorescence spectroscopy. We find, somewhat unexpectedly, that profiling the Boltzmann sample captures evidence of ligand-bound conformations for two of three riboswitches studied. Moreover, this agreement between predicted modalities and experimental models is consistent with the classification of riboswitches into thermodynamic versus kinetic regulatory mechanisms. Our results support cluster analysis of Boltzmann samples by RNAStructProfiling as a possible basis for de novo identification of thermodynamic riboswitches, while highlighting the challenges for kinetic ones.

Development of an inducible anti-VEGF rAAV gene therapy strategy for the treatment of wet AMD

Vascular endothelial growth factor (VEGF) is a key mediator in the development and progression of choroidal neovascularization (CNV) in patients with wet age-related macular degeneration (AMD). As a consequence, current treatment strategies typically focus on the administration of anti-VEGF agents, such as Aflibercept (Eylea), that inhibit VEGF function. While this approach is largely successful at counteracting CNV progression, the treatment can require repetitive (i.e. monthly) intravitreal injections of the anti-VEGF agent throughout the patient’s lifetime, imposing a substantial financial and medical burden on the patient. Moreover, repetitive injection of anti-VEGF agents over a period of years may encourage progression of retinal and choroidal atrophy in patients with AMD, leading to a decrease in visual acuity. Herein, we have developed a single-injection recombinant adeno-associated virus (rAAV)-based gene therapy treatment for wet AMD that prevents CNV formation through inducible over-expression of Eylea. First, we demonstrate that by incorporating riboswitch elements into the rAAV expression cassette allows protein expression levels to be modulated in vivo through oral supplementation on an activating ligand (e.g. tetracycline). We subsequently utilized this technology to modulate the intraocular concentration of Eylea following rAAV delivery, leading to nearly complete (p = 0.0008) inhibition of clinically significant CNV lesions in an established mouse model of wet AMD. The results shown in this study pave the way for the development of a personalized gene therapy strategy for the treatment of wet AMD that is substantially less invasive and more clinically adaptable than the current treatment paradigm of repetitive bolus injections of anti-VEGF agents.

In silico design of ligand triggered RNA switches

This contribution sketches a work flow to design an RNA switch that is able to adapt two structural conformations in a ligand-dependent way. A well characterized RNA aptamer, i.e., knowing its K_d and adaptive structural features, is an essential ingredient of the described design process. We exemplify the principles using the well-known theophylline aptamer throughout this work. The aptamer in its ligand-binding competent structure represents one structural conformation of the switch while an alternative fold that disrupts the binding-competent structure forms the other conformation. To keep it simple we do not incorporate any regulatory mechanism to control transcription or translation. We elucidate a commonly used design process by explicitly dissecting and explaining the necessary steps in detail. We developed a novel objective function which specifies the mechanistics of this simple, ligand-triggered riboswitch and describe an extensive in silico analysis pipeline to evaluate important kinetic properties of the designed sequences. This protocol and the developed software can be easily extended or adapted to fit novel design scenarios and thus can serve as a template for future needs.

Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch

Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.

Synthetic Riboswitches: From Plug and Pray toward Plug and Play

In synthetic biology, metabolic engineering, and gene therapy, there is a strong demand for orthogonal or externally controlled regulation of gene expression. Here, RNA-based regulatory devices represent a promising emerging alternative to proteins, allowing a fast and direct control of gene expression, as no synthesis of regulatory proteins is required. Besides programmable ribozyme elements controlling mRNA stability, regulatory RNA structures in untranslated regions are highly interesting for engineering approaches. Riboswitches are especially well suited, as they show a modular composition of sensor and response elements, allowing a free combination of different modules in a plug-and-play-like mode. The sensor or aptamer domain specifically interacts with a trigger molecule as a ligand, modulating the activity of the adjacent response domain that controls the expression of the genes located downstream, in most cases at the level of transcription or translation. In this review, we discuss the recent advances and strategies for designing such synthetic riboswitches based on natural or artificial components and readout systems, from trial-and-error approaches to rational design strategies. As the past several years have shown dramatic development in this fascinating field of research, we can give only a limited overview of the basic riboswitch design principles that is far from complete, and we apologize for not being able to consider every successful and interesting approach described in the literature.

Changes in the Plasticity of HIV-1 Nef RNA during the Evolution of the North American Epidemic

Because of a high mutation rate, HIV exists as a viral swarm of many sequence variants evolving under various selective pressures from the human immune system. Although the Nef gene codes for the most immunogenic of HIV accessory proteins, which alone makes it of great interest to HIV research, it also encodes an RNA structure, whose contribution to HIV virulence has been largely unexplored. Nef RNA helps HIV escape RNA interference (RNAi) through nucleotide changes and alternative folding. This study examines Historic and Modern Datasets of patient HIV-1 Nef sequences during the evolution of the North American epidemic for local changes in RNA plasticity. By definition, RNA plasticity refers to an RNA molecule’s ability to take alternative folds (i.e., alternative conformations). Our most important finding is that an evolutionarily conserved region of the HIV-1 Nef gene, which we denote by R2, recently underwent a statistically significant increase in its RNA plasticity. Thus, our results indicate that Modern Nef R2 typically accommodates an alternative fold more readily than Historic Nef R2. Moreover, the increase in RNA plasticity resides mostly in synonymous nucleotide changes, which cannot be a response to selective pressures on the Nef protein. R2 may therefore be of interest in the development of antiviral RNAi therapies.

SwiSpot: modeling riboswitches by spotting out switching sequences

MOTIVATION

Riboswitches are cis-regulatory elements in mRNA, mostly found in Bacteria, which exhibit two main secondary structure conformations. Although one of them prevents the gene from being expressed, the other conformation allows its expression, and this switching process is typically driven by the presence of a specific ligand. Although there are a handful of known riboswitches, our knowledge in this field has been greatly limited due to our inability to identify their alternate structures from their sequences. Indeed, current methods are not able to predict the presence of the two functionally distinct conformations just from the knowledge of the plain RNA nucleotide sequence. Whether this would be possible, for which cases, and what prediction accuracy can be achieved, are currently open questions.

RESULTS

Here we show that the two alternate secondary structures of riboswitches can be accurately predicted once the 'switching sequence' of the riboswitch has been properly identified. The proposed SwiSpot approach is capable of identifying the switching sequence inside a putative, complete riboswitch sequence, on the basis of pairing behaviors, which are evaluated on proper sets of configurations. Moreover, it is able to model the switching behavior of riboswitches whose generated ensemble covers both alternate configurations. Beyond structural predictions, the approach can also be paired to homology-based riboswitch searches.

AVAILABILITY AND IMPLEMENTATION

SwiSpot software, along with the reference dataset files, is available at: http://www.iet.unipi.it/a.bechini/swispot/Supplementary information: Supplementary data are available at Bioinformatics online.

CONTACT

a.bechini@ing.unipi.it.

Mn(2+)-sensing mechanisms of yybP-ykoY orphan riboswitches

Gene regulation in cis by riboswitches is prevalent in bacteria. The yybP-ykoY riboswitch family is quite widespread, yet its ligand and function remained unknown. Here, we characterize the Lactococcus lactis yybP-ykoY orphan riboswitch as a Mn(2+)-dependent transcription-ON riboswitch, with a ∼30-40 μM affinity for Mn(2+). We further determined its crystal structure at 2.7 Å to elucidate the metal sensing mechanism. The riboswitch resembles a hairpin, with two coaxially stacked helices tethered by a four-way junction and a tertiary docking interface. The Mn(2+)-sensing region, strategically located at the highly conserved docking interface, has two metal binding sites. Whereas one site tolerates the binding of either Mg(2+) or Mn(2+), the other site strongly prefers Mn(2+) due to a direct contact from the N7 of an invariable adenosine. Mutagenesis and a Mn(2+)-free E. coli yybP-ykoY structure further reveal that Mn(2+) binding is coupled with stabilization of the Mn(2+)-sensing region and the aptamer domain.

Linking aptamer-ligand binding and expression platform folding in riboswitches: prospects for mechanistic modeling and design

The power of riboswitches in regulation of bacterial metabolism derives from coupling of two characteristics: recognition and folding. Riboswitches contain aptamers, which function as biosensors. Upon detection of the signaling molecule, the riboswitch transduces the signal into a genetic decision. The genetic decision is coupled to refolding of the expression platform, which is distinct from, although overlapping with, the aptamer. Early biophysical studies of riboswitches focused on recognition of the ligand by the aptamer-an important consideration for drug design. A mechanistic understanding of ligand-induced riboswitch RNA folding can further enhance riboswitch ligand design, and inform efforts to tune and engineer riboswitches with novel properties. X-ray structures of aptamer/ligand complexes point to mechanisms through which the ligand brings together distal strand segments to form a P1 helix. Transcriptional riboswitches must detect the ligand and form this P1 helix within the timescale of transcription. Depending on the cell's metabolic state and cellular environmental conditions, the folding and genetic outcome may therefore be affected by kinetics of ligand binding, RNA folding, and transcriptional pausing, among other factors. Although some studies of isolated riboswitch aptamers found homogeneous, prefolded conformations, experimental, and theoretical studies point to functional and structural heterogeneity for nascent transcripts. Recently it has been shown that some riboswitch segments, containing the aptamer and partial expression platforms, can form binding-competent conformers that incorporate an incomplete aptamer secondary structure. Consideration of the free energy landscape for riboswitch RNA folding suggests models for how these conformers may act as transition states-facilitating rapid, ligand-mediated aptamer folding.

Predicting Large RNA-Like Topologies by a Knowledge-Based Clustering Approach

An analysis and expansion of our resource for classifying, predicting, and designing RNA structures, RAG (RNA-As-Graphs), is presented, with the goal of understanding features of RNA-like and non-RNA-like motifs and exploiting this information for RNA design. RAG was first reported in 2004 for cataloging RNA secondary structure motifs using graph representations. In 2011, the RAG resource was updated with the increased availability of RNA structures and was improved by utilities for analyzing RNA structures, including substructuring and search tools. We also classified RNA structures as graphs up to 10 vertices (~200 nucleotides) into three classes: existing, RNA-like, and non-RNA-like using clustering approaches. Here, we focus on the tree graphs and evaluate the newly founded RNAs since 2011, which also support our refined predictions of RNA-like motifs. We expand the RAG resource for large tree graphs up to 13 vertices (~260 nucleotides), thereby cataloging more than 10 times as many secondary structures. We apply clustering algorithms based on features of RNA secondary structures translated from known tertiary structures to suggest which hypothetical large RNA motifs can be considered "RNA-like". The results by the PAM (Partitioning Around Medoids) approach, in particular, reveal good accuracy, with small error for the largest cases. The RAG update here up to 13 vertices offers a useful graph-based tool for exploring RNA motifs and suggesting large RNA motifs for design.

Secondary structural entropy in RNA switch (Riboswitch) identification

BACKGROUND

RNA regulatory elements play a significant role in gene regulation. Riboswitches, a widespread group of regulatory RNAs, are vital components of many bacterial genomes. These regulatory elements generally function by forming a ligand-induced alternative fold that controls access to ribosome binding sites or other regulatory sites in RNA. Riboswitch-mediated mechanisms are ubiquitous across bacterial genomes. A typical class of riboswitch has its own unique structural and biological complexity, making de novo riboswitch identification a formidable task. Traditionally, riboswitches have been identified through comparative genomics based on sequence and structural homology. The limitations of structural-homology-based approaches, coupled with the assumption that there is a great diversity of undiscovered riboswitches, suggests the need for alternative methods for riboswitch identification, possibly based on features intrinsic to their structure. As of yet, no such reliable method has been proposed.

RESULTS

We used structural entropy of riboswitch sequences as a measure of their secondary structural dynamics. Entropy values of a diverse set of riboswitches were compared to that of their mutants, their dinucleotide shuffles, and their reverse complement sequences under different stochastic context-free grammar folding models. Significance of our results was evaluated by comparison to other approaches, such as the base-pairing entropy and energy landscapes dynamics. Classifiers based on structural entropy optimized via sequence and structural features were devised as riboswitch identifiers and tested on Bacillus subtilis, Escherichia coli, and Synechococcus elongatus as an exploration of structural entropy based approaches. The unusually long untranslated region of the cotH in Bacillus subtilis, as well as upstream regions of certain genes, such as the sucC genes were associated with significant structural entropy values in genome-wide examinations.

CONCLUSIONS

Various tests show that there is in fact a relationship between higher structural entropy and the potential for the RNA sequence to have alternative structures, within the limitations of our methodology. This relationship, though modest, is consistent across various tests. Understanding the behavior of structural entropy as a fairly new feature for RNA conformational dynamics, however, may require extensive exploratory investigation both across RNA sequences and folding models.

Dynamic signal processing by ribozyme-mediated RNA circuits to control gene expression

Organisms have different circuitries that allow converting signal molecule levels to changes in gene expression. An important challenge in synthetic biology involves the de novo design of RNA modules enabling dynamic signal processing in live cells. This requires a scalable methodology for sensing, transmission, and actuation, which could be assembled into larger signaling networks. Here, we present a biochemical strategy to design RNA-mediated signal transduction cascades able to sense small molecules and small RNAs. We design switchable functional RNA domains by using strand-displacement techniques. We experimentally characterize the molecular mechanism underlying our synthetic RNA signaling cascades, show the ability to regulate gene expression with transduced RNA signals, and describe the signal processing response of our systems to periodic forcing in single live cells. The engineered systems integrate RNA–RNA interaction with available ribozyme and aptamer elements, providing new ways to engineer arbitrary complex gene circuits.

Using simulations and kinetic network models to reveal the dynamics and functions of riboswitches

Riboswitches, RNA elements found in the untranslated region, regulate gene expression by binding to target metaboloites with exquisite specificity. Binding of metabolites to the conserved aptamer domain allosterically alters the conformation in the downstream expression platform. The fate of gene expression is determined by the changes in the downstream RNA sequence. As the metabolite-dependent cotranscriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate both the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolite. Single molecule force experiments that decipher the free energy landscape of riboswitches from their mechanical responses, theoretical and computational studies have recently shed light on the distinct mechanism of folding dynamics in different classes of riboswitches. Here, we first discuss the dynamics of water around riboswitch, highlighting that water dynamics can enhance the fluctuation of nucleic acid structure. To go beyond native state fluctuations, we used the Self-Organized Polymer model to predict the dynamics of add adenine riboswitch under mechanical forces. In addition to quantitatively predicting the folding landscape of add-riboswitch, our simulations also explain the difference in the dynamics between pbuE adenine- and add adenine-riboswitches. In order to probe the function in vivo, we use the folding landscape to propose a system level kinetic network model to quantitatively predict how gene expression is regulated for riboswitches that are under kinetic control.

Design of transcription regulating riboswitches

In this chapter, we review both computational and experimental aspects of de novo RNA sequence design. We give an overview of currently available design software and their limitations, and discuss the necessary setup to experimentally validate proper function in vitro and in vivo. We focus on transcription-regulating riboswitches, a task that has just recently lead to first successful designs of such RNA elements.

RNA graph partitioning for the discovery of RNA modularity: a novel application of graph partition algorithm to biology

Graph representations have been widely used to analyze and design various economic, social, military, political, and biological networks. In systems biology, networks of cells and organs are useful for understanding disease and medical treatments and, in structural biology, structures of molecules can be described, including RNA structures. In our RNA-As-Graphs (RAG) framework, we represent RNA structures as tree graphs by translating unpaired regions into vertices and helices into edges. Here we explore the modularity of RNA structures by applying graph partitioning known in graph theory to divide an RNA graph into subgraphs. To our knowledge, this is the first application of graph partitioning to biology, and the results suggest a systematic approach for modular design in general. The graph partitioning algorithms utilize mathematical properties of the Laplacian eigenvector (µ2) corresponding to the second eigenvalues (λ2) associated with the topology matrix defining the graph: λ2 describes the overall topology, and the sum of µ2's components is zero. The three types of algorithms, termed median, sign, and gap cuts, divide a graph by determining nodes of cut by median, zero, and largest gap of µ2's components, respectively. We apply these algorithms to 45 graphs corresponding to all solved RNA structures up through 11 vertices (∼ 220 nucleotides). While we observe that the median cut divides a graph into two similar-sized subgraphs, the sign and gap cuts partition a graph into two topologically-distinct subgraphs. We find that the gap cut produces the best biologically-relevant partitioning for RNA because it divides RNAs at less stable connections while maintaining junctions intact. The iterative gap cuts suggest basic modules and assembly protocols to design large RNA structures. Our graph substructuring thus suggests a systematic approach to explore the modularity of biological networks. In our applications to RNA structures, subgraphs also suggest design strategies for novel RNA motifs.

A family of synthetic riboswitches adopts a kinetic trapping mechanism

Riboswitches are sequences of RNA that control gene expression via RNA–ligand interactions, without the need for accessory proteins. Riboswitches consist of an aptamer that recognizes the ligand and an expression platform that couples ligand binding to a change in gene expression. Using in vitro selection, it is possible to screen large (∼10¹³ members) libraries of RNA sequences to discover new aptamers. However, limitations in bacterial transformation efficiency make screening such large libraries for riboswitch function in intact cells impractical. Here we show that synthetic riboswitches function in an E. coli S30 extract in a manner similar to how they function in intact E. coli cells. We discovered that, although this family of riboswitches regulates the initiation of protein translation, the fate of whether an RNA message is translated is determined during transcription. Thus, ligand binding does not bias a population of rapidly equilibrating RNA structures, but rather, co-transcriptional ligand binding kinetically traps the RNA in a conformation that supports efficient translation. In addition to providing new insights into the mechanisms of action of a family of synthetic riboswitches, our experiments suggest that it may be possible to perform selections for novel synthetic riboswitches in an in vitro system.

Single-molecule studies of riboswitch folding

The folding dynamics of riboswitches are central to their ability to modulate gene expression in response to environmental cues. In most cases, a structural competition between the formation of a ligand-binding aptamer and an expression platform (or some other competing off-state) determines the regulatory outcome. Here, we review single-molecule studies of riboswitch folding and function, predominantly carried out using single-molecule FRET or optical trapping approaches. Recent results have supplied new insights into riboswitch folding energy landscapes, the mechanisms of ligand binding, the roles played by divalent ions, the applicability of hierarchical folding models, and kinetic vs. thermodynamic control schemes. We anticipate that future work, based on improved data sets and potentially combining multiple experimental techniques, will enable the development of more complete models for complex RNA folding processes. This article is part of a Special Issue entitled: Riboswitches.

Differences between cotranscriptional and free riboswitch folding

Riboswitches are part of noncoding regions of messenger RNA (mRNA) that act as RNA sensors regulating gene expression of the downstream gene. Typically, one out of two distinct conformations is formed depending on ligand binding when the transcript leaves RNA polymerase (RNAP). Elongation of the RNA chain by RNAP, folding and binding all occurs simultaneously and interdependently on the seconds' timescale. To investigate the effect of transcript elongation velocity on folding for the S-adenosylmethionine (SAM)-I and adenine riboswitches we employ two complementary coarse-grained in silico techniques. Native structure-based molecular dynamics simulations provide a 3D, atomically resolved model of folding with homogenous energetics. Energetically more detailed kinetic Monte Carlo simulations give access to longer timescale by describing folding on the secondary structure level and feature the incorporation of competing aptamer conformations and a ligand-binding model. Depending on the extrusion scenarios, we observe and quantify different pathways in structure formation with robust agreements between the two techniques. In these scenarios, free-folding riboswitches exhibit different folding characteristics compared with transcription-rate limited folding. The critical transcription rate distinguishing these cases is higher than physiologically relevant rates. This result suggests that in vivo folding of the analyzed SAM-I and adenine riboswitches is transcription-rate limited.

Kinetics of allosteric transitions in S-adenosylmethionine riboswitch are accurately predicted from the folding landscape

Riboswitches are RNA elements that allosterically regulate gene expression by binding cellular metabolites. The SAM-III riboswitch, one of several classes that binds S-adenosylmethionine (SAM), represses translation upon binding SAM (OFF state) by encrypting the ribosome binding sequence. We have carried out simulations of the RNA by applying mechanical force (f) to the ends of SAM-III, with and without SAM, to get quantitative insights into the f-dependent structural changes. Force-extension (z) curves (FECs) for the apo (ON) state, obtained in simulations in which f is increased at a constant loading rate, show three intermediates, with the first one being the rupture of SAM binding region, which is greatly stabilized in the OFF state. Force-dependent free energy profiles, G(z,f), as a function of z, obtained in equilibrium constant force simulations, reveal the intermediates observed in FECs. The predicted stability difference between the ON and OFF states using G(z,f) is in excellent agreement with experiments. Remarkably, using G(z,f)s and estimate of an effective diffusion constant at a single value of f allows us to predict the f-dependent transition rates using theory of first passage times for both the apo and holo states. To resolve the kinetics of assembly of SAM-III riboswitch in structural terms, we use force stretch-quench pulse sequences in which the force on RNA is maintained at a low (fq) value starting from a high value for a time period tq. Variation of tq over a wide range results in resolution of elusive states involved in the SAM binding pocket and leads to accurate determination of folding times down to fq = 0. Quantitative measure of the folding kinetics, obtained from the folding landscape, allows us to propose that, in contrast to riboswitches regulating transcription, SAM-III functions under thermodynamic control provided the basal concentration of SAM exceeds a small critical value. All of the predictions are amenable to tests in single molecule pulling experiments.

Network Theory Tools for RNA Modeling

An introduction into the usage of graph or network theory tools for the study of RNA molecules is presented. By using vertices and edges to define RNA secondary structures as tree and dual graphs, we can enumerate, predict, and design RNA topologies. Graph connectivity and associated Laplacian eigenvalues relate to biological properties of RNA and help understand RNA motifs as well as build, by computational design, various RNA target structures. Importantly, graph theoretical representations of RNAs reduce drastically the conformational space size and therefore simplify modeling and prediction tasks. Ongoing challenges remain regarding general RNA design, representation of RNA pseudoknots, and tertiary structure prediction. Thus, developments in network theory may help advance RNA biology.

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

Pseudoknot preorganization of the preQ1 class I riboswitch

To explore folding and ligand recognition of metabolite-responsive RNAs is of major importance to comprehend gene regulation by riboswitches. Here, we demonstrate, using NMR spectroscopy, that the free aptamer of a preQ(1) class I riboswitch preorganizes into a pseudoknot fold in a temperature- and Mg(2+)-dependent manner. The preformed pseudoknot represents a structure that is close to the ligand-bound state and that likely represents the conformation selected by the ligand. Importantly, a defined base pair mutation within the pseudoknot interaction stipulates whether, in the absence of ligand, dimer formation of the aptamer competes with intramolecular pseudoknot formation. This study pinpoints how RNA preorganization is a crucial determinant for the adaptive recognition process of RNA and ligand.