Design Of Riboswitches Research Articles

Partial RNA design.

RNA design is a key technique to achieve new functionality in fields like synthetic biology or biotechnology. Computational tools could help to find such RNA sequences but they are often limited in their formulation of the search space. In this work, we propose partial RNA design, a novel RNA design paradigm that addresses the limitations of current RNA design formulations. Partial RNA design describes the problem of designing RNAs from arbitrary RNA sequences and structure motifs with multiple design goals. By separating the design space from the objectives, our formulation enables the design of RNAs with variable lengths and desired properties, while still allowing precise control over sequence and structure constraints at individual positions. Based on this formulation, we introduce a new algorithm, libLEARNA, capable of efficiently solving different constraint RNA design tasks. A comprehensive analysis of various problems, including a realistic riboswitch design task, reveals the outstanding performance of libLEARNA and its robustness. libLEARNA is open-source and publicly available at: https://github.com/automl/learna_tools.

Broad-Spectrum Detection of Tetracyclines by Riboswitch-Based Cell-Free Expression Biosensing.

We describe a sensitive and selective method for the determination of tetracycline content in foods using a riboswitch sensor. The sensor is based on a cell-free expression system that can be lyophilized to produce paper-based sensors or tube-based sensors for long-term storage. The riboswitch constructed using artificially screened tetracycline RNA aptamers was cloned into the pET-28a(+) vector of Escherichia coli TOP 10. The expression of the green fluorescent protein was positively correlated with the concentration of tetracyclines. The binding of tetracyclines to the aptamer domain results in a conformational change in the riboswitch secondary structure, resulting in the exposure of the ribosome binding site thereby promoting expression. The detection limits of the prepared sensor for the detection of tetracycline, oxytetracycline, chlortetracycline, and doxycycline were 0.47, 0.079, 0.084, and 0.43 μM, respectively. Moreover, the 1 μM tetracyclines allow for qualitative detection in milk samples by the naked eye. The work provides a proof-of-principle for riboswitch design to address global health and food safety.

Linking RNA Sequence, Structure, and Function on Massively Parallel High-Throughput Sequencers.

High-throughput sequencing methods have revolutionized our ability to catalog the diversity of RNAs and RNA-protein interactions that can exist in our cells. However, the relationship between RNA sequence, structure, and function is enormously complex, demonstrating the need for methods that can provide quantitative thermodynamic and kinetic measurements of macromolecular interaction with RNA, at a scale commensurate with the sequence diversity of RNA. Here, we discuss a class of methods that extend the core functionality of DNA sequencers to enable high-throughput measurements of RNA folding and RNA-protein interactions. Topics discussed include a description of the method and multiple applications to RNA-binding proteins, riboswitch design and engineering, and RNA tertiary structure energetics.

Tractable RNA\u2013ligand interaction kinetics

BackgroundThe binding of small ligands to RNA elements can cause substantial changes in the RNA structure. This constitutes an important, fast-acting mechanism of ligand-controlled transcriptional and translational gene regulation implemented by a wide variety of riboswitches. The associated refolding processes often cannot be explained by thermodynamic effects alone. Instead, they are governed by the kinetics of RNA folding. While the computational analysis of RNA folding can make use of well-established models of the thermodynamics of RNA structures formation, RNA–RNA interaction, and RNA–ligand interaction, kinetic effects pose fundamentally more challenging problems due to the enormous size of the conformation space. The analysis of the combined process of ligand binding and structure formation even for small RNAs is plagued by intractably large state spaces. Moreover, the interaction is concentration-dependent and thus is intrinsically non-linear. This precludes the direct transfer of the strategies previously used for the analysis of RNA folding kinetics.ResultsIn our novel, computationally tractable approach to RNA–ligand kinetics, we overcome the two main difficulties by applying a gradient-based coarse graining to RNA–ligand systems and solving the process in a pseudo-first order approximation. The latter is well-justified for the most common case of ligand excess in RNA–ligand systems. We present the approach rigorously and discuss the parametrization of the model based on empirical data. The method supports the kinetic study of RNA–ligand systems, in particular at different ligand concentrations. As an example, we apply our approach to analyze the concentration dependence of the ligand response of the rationally designed, artificial theophylline riboswitch RS3.ConclusionThis work demonstrates the tractability of the computational analysis of RNA–ligand interaction. Naturally, the model will profit as more accurate measurements of folding and binding parameters become available. Due to this work, computational analysis is available to support tasks like the design of riboswitches; our analysis of RS3 suggests strong co-transcriptional effects for this riboswitch.The method used in this study is available online, cf. Section “Availability of data and materials”.

StreAM-T_g: algorithms for analyzing coarse grained RNA dynamics based on Markov models of connectivity-graphs

BackgroundIn this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-T_g, a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA.ResultsWe benchmark StreAM-T_g (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-T_g on six long term molecular dynamics simulation of a synthetic tetracycline dependent riboswitch (500 ns) in combination with five different antibiotics.ConclusionsThe proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-T_g provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.

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.

Dynamic Energy Landscapes of Riboswitches Help Interpret Conformational Rearrangements and Function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.