SnapShot: RNA Structure Probing Technologies.

Chemical probing coupled to high-throughput sequencing offers a flexible approach to uncover many aspects of RNA structure relevant to its cellular function. With a wide variety of chemical probes available that each report on different features of RNA molecules, a broad toolkit exists for investigating in vivo and in vitro RNA structure and interactions with other molecules.

Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates.

The series of RNA folding events that occur during transcription can critically influence cellular RNA function. Here, we present reconstructing RNA dynamics from data (R2D2), a method to uncover details of cotranscriptional RNA folding. We model the folding of the Escherichia coli signal recognition particle (SRP) RNA and show that it requires specific local structural fluctuations within a key hairpin to engender efficient cotranscriptional conformational rearrangement into the functional structure. All-atom molecular dynamics simulations suggest that this rearrangement proceeds through an internal toehold-mediated strand-displacement mechanism, which can be disrupted with a point mutation that limits local structural fluctuations and rescued with compensating mutations that restore these fluctuations. Moreover, a cotranscriptional folding intermediate could be cleaved in vitro by recombinant E. coli RNase P, suggesting potential cotranscriptional processing. These results from experiment-guided multi-scale modeling demonstrate that even an RNA with a simple functional structure can undergo complex folding and processing during synthesis.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch.

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

A transient intermediate RNA structure underlies the regulatory function of the E. coli thiB TPP translational riboswitch.

Riboswitches are cis-regulatory RNA elements that regulate gene expression in response to ligand binding through the coordinated action of a ligand-binding aptamer domain (AD) and a downstream expression platform (EP). Previous studies of transcriptional riboswitches have uncovered diverse examples that utilize structural intermediates that compete with the AD and EP folds to mediate the switching mechanism on the timescale of transcription. Here we investigate whether similar intermediates are important for riboswitches that control translation by studying the Escherichia coli thiB thiamin pyrophosphate (TPP) riboswitch. Using cellular gene expression assays, we first confirmed that the riboswitch acts at the level of translational regulation. Deletion mutagenesis showed the importance of the AD-EP linker sequence for riboswitch function. Sequence complementarity between the linker region and the AD P1 stem suggested the possibility of an intermediate nascent RNA structure called the antisequestering stem that could mediate the thiB switching mechanism. Experimentally informed secondary structure models of the thiB folding pathway generated from chemical probing of nascent thiB structures in stalled transcription elongation complexes confirmed the presence of the antisequestering stem, and showed it may form cotranscriptionally. Additional mutational analysis showed that mutations to the antisequestering stem break or bias thiB function according to whether the antisequestering stem or P1 is favored. This work provides an important example of intermediate structures that compete with AD and EP folds to implement riboswitch mechanisms.

Tuning strand displacement kinetics enables programmable ZTP riboswitch dynamic range in vivo.

A large body of work has shown that transcriptional riboswitches function through internal strand displacement mechanisms that guide the formation of alternative structures which drive regulatory outcomes. Here, we sought to investigate this phenomenon using the Clostridium beijerinckii pfl ZTP riboswitch as a model system. Using functional mutagenesis with Escherichia coli gene expression assays, we show that mutations designed to slow strand displacement of the expression platform enable precise tuning of riboswitch dynamic range (2.4-34-fold), depending on the type of kinetic barrier introduced, and the position of the barrier relative to the strand displacement nucleation site. We also show that expression platforms from a range of different Clostridium ZTP riboswitches contain sequences that impose these barriers to affect dynamic range in these different contexts. Finally, we use sequence design to flip the regulatory logic of the riboswitch to create a transcriptional OFF-switch, and show that the same barriers to strand displacement tune dynamic range in this synthetic context. Together, our findings further elucidate how strand displacement can be manipulated to alter the riboswitch decision landscape, suggesting that this could be a mechanism by which evolution tunes riboswitch sequence, and providing an approach to optimize synthetic riboswitches for biotechnology applications.

Programming cell-free biosensors with DNA strand displacement circuits.

Cell-free biosensors are powerful platforms for monitoring human and environmental health. Here, we expand their capabilities by interfacing them with toehold-mediated strand displacement circuits, a dynamic DNA nanotechnology that enables molecular computation through programmable interactions between nucleic acid strands. We develop design rules for interfacing a small molecule sensing platform called ROSALIND with toehold-mediated strand displacement to construct hybrid RNA-DNA circuits that allow fine-tuning of reaction kinetics. We use these design rules to build 12 different circuits that implement a range of logic functions (NOT, OR, AND, IMPLY, NOR, NIMPLY, NAND). Finally, we demonstrate a circuit that acts like an analog-to-digital converter to create a series of binary outputs that encode the concentration range of the molecule being detected. We believe this work establishes a pathway to create 'smart' diagnostics that use molecular computations to enhance the speed and utility of biosensors.

High-throughput determination of RNA structures.

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Cotranscriptional RNA strand exchange underlies the gene regulation mechanism in a purine-sensing transcriptional riboswitch.

RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.

Riboswitches have challenged our understanding of biological regulation for almost two decades. The ability of small molecules to bind to RNA and control gene expression offers another layer of regulation and the potential for direct action by compounds in the environment. While some riboswitches have been well studied, we lack a general understanding of how changes in RNA structure switch genetic expression from "On" to "Off". In this study, the authors propose an elegant "strand displacement" model to explain how the RNA structure shifts between "On" and "Off" states as the concentration of small molecule ligand changes. These observations help us to understand how riboswitches enable genetic decision-making. The data provide a possible general mechanism for understanding how the competition between different strand displacement outcomes can influence RNA folding. Understanding RNA folding pathways could advance the successful design of drugs that target RNA.

Zur and zinc increase expression of E. coli ribosomal protein L31 through RNA-mediated repression of the repressor L31p.

Bacteria can adapt in response to numerous stress conditions. One such stress condition is zinc depletion. The zinc-sensing transcription factor Zur regulates the way numerous bacterial species respond to severe changes in zinc availability. Under zinc sufficient conditions, Zn-loaded Zur (Zn2-Zur) is well-known to repress transcription of genes encoding zinc uptake transporters and paralogues of a few ribosomal proteins. Here, we report the discovery and mechanistic basis for the ability of Zur to up-regulate expression of the ribosomal protein L31 in response to zinc in E. coli. Through genetic mutations and reporter gene assays, we find that Zur achieves the up-regulation of L31 through a double repression cascade by which Zur first represses the transcription of L31p, a zinc-lacking paralogue of L31, which in turn represses the translation of L31. Mutational analyses show that translational repression by L31p requires an RNA hairpin structure within the l31 mRNA and involves the N-terminus of the L31p protein. This work uncovers a new genetic network that allows bacteria to respond to host-induced nutrient limiting conditions through a sophisticated ribosomal protein switching mechanism.

Dynamic RNA synthetic biology: new principles, practices and potential.

An increased appreciation of the role of RNA dynamics in governing RNA function is ushering in a new wave of dynamic RNA synthetic biology. Here, we review recent advances in engineering dynamic RNA systems across the molecular, circuit and cellular scales for important societal-scale applications in environmental and human health, and bioproduction. For each scale, we introduce the core concepts of dynamic RNA folding and function at that scale, and then discuss technologies incorporating these concepts, covering new approaches to engineering riboswitches, ribozymes, RNA origami, RNA strand displacement circuits, biomaterials, biomolecular condensates, extracellular vesicles and synthetic cells. Considering the dynamic nature of RNA within the engineering design process promises to spark the next wave of innovation that will expand the scope and impact of RNA biotechnologies.

Chemical roadblocking of DNA transcription for nascent RNA display.

Site-specific arrest of RNA polymerases (RNAPs) is fundamental to several technologies that assess RNA structure and function. Current in vitro transcription "roadblocking" approaches inhibit transcription elongation by blocking RNAP with a protein bound to the DNA template. One limitation of protein-mediated transcription roadblocking is that it requires inclusion of a protein factor extrinsic to the minimal in vitro transcription reaction. In this work, we developed a chemical approach for halting transcription by Escherichia coli RNAP. We first established a sequence-independent method for site-specific incorporation of chemical lesions into dsDNA templates by sequential PCR and translesion synthesis. We then show that interrupting the transcribed DNA strand with an internal desthiobiotin-triethylene glycol modification or 1,N6-etheno-2'-deoxyadenosine base efficiently and stably halts Escherichia coli RNAP transcription. By encoding an intrinsic stall site within the template DNA, our chemical transcription roadblocking approach enables display of nascent RNA molecules from RNAP in a minimal in vitro transcription reaction.

A ligand-gated strand displacement mechanism for ZTP riboswitch transcription control.

Cotranscriptional folding is an obligate step of RNA biogenesis that can guide RNA structure formation and function through transient intermediate folds. This process is particularly important for transcriptional riboswitches in which the formation of ligand-dependent structures during transcription regulates downstream gene expression. However, the intermediate structures that comprise cotranscriptional RNA folding pathways, and the mechanisms that enable transit between them, remain largely unknown. Here, we determine the series of cotranscriptional folds and rearrangements that mediate antitermination by the Clostridium beijerinckii pfl ZTP riboswitch in response to the purine biosynthetic intermediate ZMP. We uncover sequence and structural determinants that modulate an internal RNA strand displacement process and identify biases within natural ZTP riboswitch sequences that promote on-pathway folding. Our findings establish a mechanism for pfl riboswitch antitermination and suggest general strategies by which nascent RNA molecules navigate cotranscriptional folding pathways.

De novo-designed translation-repressing riboregulators for multi-input cellular logic.

Efforts to construct synthetic biological circuits with more complex functions have often been hindered by the idiosyncratic behavior, limited dynamic range and crosstalk of commonly utilized parts. Here, we employ de novo RNA design to develop two high-performance translational repressors with sensing and logic capabilities. These synthetic riboregulators, termed toehold repressors and three-way junction (3WJ) repressors, detect transcripts with nearly arbitrary sequences, repress gene expression by up to 300-fold and yield orthogonal sets of up to 15 devices. Automated forward engineering is used to improve toehold repressor dynamic range and SHAPE-Seq is applied to confirm the designed switching mechanism of 3WJ repressors in living cells. We integrate the modular repressors into biological circuits that execute universal NAND and NOR logic and evaluate the four-input expression NOT ((A1 AND A2) OR (B1 AND B2)) in Escherichia coli. These capabilities make toehold and 3WJ repressors valuable new tools for biotechnological applications.

DUETT quantitatively identifies known and novel events in nascent RNA structural dynamics from chemical probing data.

MOTIVATION: RNA molecules can undergo complex structural dynamics, especially during transcription, which influence their biological functions. Recently developed high-throughput chemical probing experiments that study RNA cotranscriptional folding generate nucleotide-resolution 'reactivities' for each length of a growing nascent RNA that reflect structural dynamics. However, the manual annotation and qualitative interpretation of reactivity across these large datasets can be nuanced, laborious, and difficult for new practitioners. We developed a quantitative and systematic approach to automatically detect RNA folding events from these datasets to reduce human bias/error, standardize event discovery and generate hypotheses about RNA folding trajectories for further analysis and experimental validation. RESULTS: Detection of Unknown Events with Tunable Thresholds (DUETT) identifies RNA structural transitions in cotranscriptional RNA chemical probing datasets. DUETT employs a feedback control-inspired method and a linear regression approach and relies on interpretable and independently tunable parameter thresholds to match qualitative user expectations with quantitatively identified folding events. We validate the approach by identifying known RNA structural transitions within the cotranscriptional folding pathways of the Escherichia coli signal recognition particle RNA and the Bacillus cereus crcB fluoride riboswitch. We identify previously overlooked features of these datasets such as heightened reactivity patterns in the signal recognition particle RNA about 12 nt lengths before base-pair rearrangement. We then apply a sensitivity analysis to identify tradeoffs when choosing parameter thresholds. Finally, we show that DUETT is tunable across a wide range of contexts, enabling flexible application to study broad classes of RNA folding mechanisms. AVAILABILITY AND IMPLEMENTATION: https://github.com/BagheriLab/DUETT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Probing of RNA structures in a positive sense RNA virus reveals selection pressures for structural elements.

In single stranded (+)-sense RNA viruses, RNA structural elements (SEs) play essential roles in the infection process from replication to encapsidation. Using selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) and covariation analysis, we explore the structural features of the third genome segment of cucumber mosaic virus (CMV), RNA3 (2216 nt), both in vitro and in plant cell lysates. Comparing SHAPE-Seq and covariation analysis results revealed multiple SEs in the coat protein open reading frame and 3' untranslated region. Four of these SEs were mutated and serially passaged in Nicotiana tabacum plants to identify biologically selected changes to the original mutated sequences. After passaging, loop mutants showed partial reversion to their wild-type sequence and SEs that were structurally disrupted by mutations were restored to wild-type-like structures via synonymous mutations in planta. These results support the existence and selection of virus open reading frame SEs in the host organism and provide a framework for further studies on the role of RNA structure in viral infection. Additionally, this work demonstrates the applicability of high-throughput chemical probing in plant cell lysates and presents a new method for calculating SHAPE reactivities from overlapping reverse transcriptase priming sites.

Achieving large dynamic range control of gene expression with a compact RNA transcription-translation regulator.

RNA transcriptional regulators are emerging as versatile components for genetic network construction. However, these regulators suffer from incomplete repression in their OFF state, making their dynamic range less than that of their protein counterparts. This incomplete repression causes expression leak, which impedes the construction of larger synthetic regulatory networks as leak propagation can interfere with desired network function. To address this, we demonstrate how naturally derived antisense RNA-mediated transcriptional regulators can be configured to regulate both transcription and translation in a single compact RNA mechanism that functions in Escherichia coli. Using in vivo gene expression assays, we show that a combination of transcriptional termination and ribosome binding site sequestration increases repression from 85% to 98%, or activation from 10-fold to over 900-fold, in response to cognate antisense RNAs. We also show that orthogonal repressive versions of this mechanism can be created through engineering minimal antisense RNAs. Finally, to demonstrate the utility of this mechanism, we use it to reduce network leak in an RNA-only cascade. We anticipate these regulators will find broad use as synthetic biology moves beyond parts engineering to the design and construction of more sophisticated regulatory networks.

Distributed biotin-streptavidin transcription roadblocks for mapping cotranscriptional RNA folding.

RNA folding during transcription directs an order of folding that can determine RNA structure and function. However, the experimental study of cotranscriptional RNA folding has been limited by the lack of easily approachable methods that can interrogate nascent RNA structure at nucleotide resolution. To address this, we previously developed cotranscriptional selective 2΄-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) to simultaneously probe all intermediate RNA transcripts during transcription by stalling elongation complexes at catalytically dead EcoRIE111Q roadblocks. While effective, the distribution of elongation complexes using EcoRIE111Q requires laborious PCR using many different oligonucleotides for each sequence analyzed. Here, we improve the broad applicability of cotranscriptional SHAPE-Seq by developing a sequence-independent biotin-streptavidin (SAv) roadblocking strategy that simplifies the preparation of roadblocking DNA templates. We first determine the properties of biotin-SAv roadblocks. We then show that randomly distributed biotin-SAv roadblocks can be used in cotranscriptional SHAPE-Seq experiments to identify the same RNA structural transitions related to a riboswitch decision-making process that we previously identified using EcoRIE111Q. Lastly, we find that EcoRIE111Q maps nascent RNA structure to specific transcript lengths more precisely than biotin-SAv and propose guidelines to leverage the complementary strengths of each transcription roadblock in cotranscriptional SHAPE-Seq.

Using in-cell SHAPE-Seq and simulations to probe structure-function design principles of RNA transcriptional regulators.

Antisense RNA-mediated transcriptional regulators are powerful tools for controlling gene expression and creating synthetic gene networks. RNA transcriptional repressors derived from natural mechanisms called attenuators are particularly versatile, though their mechanistic complexity has made them difficult to engineer. Here we identify a new structure-function design principle for attenuators that enables the forward engineering of new RNA transcriptional repressors. Using in-cell SHAPE-Seq to characterize the structures of attenuator variants within Escherichia coli, we show that attenuator hairpins that facilitate interaction with antisense RNAs require interior loops for proper function. Molecular dynamics simulations of these attenuator variants suggest these interior loops impart structural flexibility. We further observe hairpin flexibility in the cellular structures of natural RNA mechanisms that use antisense RNA interactions to repress translation, confirming earlier results from in vitro studies. Finally, we design new transcriptional attenuators in silico using an interior loop as a structural requirement and show that they function as desired in vivo. This work establishes interior loops as an important structural element for designing synthetic RNA gene regulators. We anticipate that the coupling of experimental measurement of cellular RNA structure and function with computational modeling will enable rapid discovery of structure-function design principles for a diverse array of natural and synthetic RNA regulators.

A flow cytometric approach to engineering Escherichia coli for improved eukaryotic protein glycosylation.

A synthetic pathway for production of the eukaryotic trimannosyl chitobiose glycan (mannose3-N-acetylglucosamine2, Man3GlcNAc2) and its transfer to specific asparagine residues in target proteins was previously engineered in Escherichia coli, providing this simple microbe with the ability to perform a complex post-translational protein modification. Here, we leveraged a flow cytometric fluorescence-based assay to improve Man3GlcNAc2 glycan biosynthesis in E. coli cells. Specifically, pathway improvements were identified, including reducing pathway enzyme expression levels and overexpressing nucleotide sugar biosynthesis genes, which enhanced production of lipid-linked Man3GlcNAc2 by nearly 50-fold to 13.9 µg/L. In turn, cells producing higher levels of the Man3GlcNAc2 substrate yielded up to 10 times more glycosylated acceptor protein (to ~ 14 mg/L) than their non-optimized counterparts. These results demonstrate the use of flow cytometry screening as a powerful tool for interrogating the surfaces of glyco-engineered bacteria and identifying meaningful improvements in glycan biosynthesis. We anticipate this approach will enable further optimization of bacterial glycan biosynthesis pathways using new strain engineering tools from metabolic engineering and synthetic biology.