Precise Genomic Riboregulator Control of Metabolic Flux in Microbial Systems.

Engineered microbes can be used for producing value-added chemicals from renewable feedstocks, relieving the dependency on nonrenewable resources such as petroleum. These microbes often are composed of synthetic metabolic pathways; however, one major problem in establishing a synthetic pathway is the challenge of precisely controlling competing metabolic routes, some of which could be crucial for fitness and survival. While traditional gene deletion and/or coarse overexpression approaches do not provide precise regulation, cis-repressors (CRs) are RNA-based regulatory elements that can control the production levels of a particular protein in a tunable manner. Here, we describe a protocol for a generally applicable fluorescence-activated cell sorting technique used to isolate eight subpopulations of CRs from a semidegenerate library in Escherichia coli, followed by deep sequencing that permitted the identification of 15 individual CRs with a broad range of protein production profiles. Using these new CRs, we demonstrated a change in production levels of a fluorescent reporter by over two orders of magnitude and further showed that these CRs are easily ported from E. coli to Pseudomonas putida. We next used four CRs to tune the production of the enzyme PpsA, involved in pyruvate to phosphoenolpyruvate (PEP) conversion, to alter the pool of PEP that feeds into the shikimate pathway. In an engineered P. putida strain, where carbon flux in the shikimate pathway is diverted to the synthesis of the commodity chemical cis,cis-muconate, we found that tuning PpsA translation levels increased the overall titer of muconate. Therefore, CRs provide an approach to precisely tune protein levels in metabolic pathways and will be an important tool for other metabolic engineering efforts.

Immobilization of Cyanines in DNA Produces Systematic Increases in Fluorescence Intensity.

Cyanines are useful fluorophores for a myriad of biological labeling applications, but their interactions with biomolecules are unpredictable. Cyanine fluorescence intensity can be highly variable due to complex photoisomerization kinetics, which are exceedingly sensitive to the surrounding environment. This introduces large errors in Förster resonance energy transfer (FRET)-based experiments where fluorescence intensity is the output parameter. However, this environmental sensitivity is a strength from a biological sensing point of view if specific relationships between biomolecular structure and cyanine photophysics can be identified. We describe a set of DNA structures that modulate cyanine fluorescence intensity through the insertion of adenine or thymine bases. These structures simultaneously provide photophysical predictability and tunability. We characterize these structures using steady-state fluorescence measurements, fluorescence correlation spectroscopy (FCS), and time-resolved photoluminescence (TRPL). We find that the photoisomerization rate decreases over an order of magnitude across the adenine series, which is consistent with increasing immobilization of the cyanine moiety by the surrounding DNA structure.

Chelated Magnesium Logic Gate Regulates Riboswitch Pseudoknot Formation.

Magnesium plays a critical role in the structure, dynamics, and function of RNA. The precise microscopic effect of chelated magnesium on RNA structure is yet to be explored. Magnesium is known to act through its diffuse cloud around RNA, through the outer sphere (water-mediated), inner sphere, and often chelated ion-mediated interactions. A mechanism is proposed for the role of experimentally discovered site-specific chelated magnesium ions on the conformational dynamics of SAM-I riboswitch aptamers in bacteria. This mechanism is observed with atomistic simulations performed in a physiological mixed salt environment at a high temperature. The simulations were validated with phosphorothioate interference mapping experiments that help to identify crucial inner-sphere Mg2+ sites prescribing an appropriate initial distribution of inner- and outer-sphere magnesium ions to maintain a physiological ion concentration of monovalent and divalent salts. A concerted role of two chelated magnesium ions is newly discovered since the presence of both supports the formation of the pseudoknot. This constitutes a logical AND gate. The absence of any of these magnesium ions instigates the dissociation of long-range pseudoknot interaction exposing the inner core of the RNA. A base triple is the epicenter of the magnesium chelation effect. It allosterically controls RNA pseudoknot by bolstering the direct effect of magnesium chelation in protecting the functional fold of RNA to control ON and OFF transcription switching.

Zinc-finger protein CNBP alters the 3-D structure of lncRNA Braveheart in solution.

Long non-coding RNAs (lncRNAs) constitute a significant fraction of the transcriptome, playing important roles in development and disease. However, our understanding of structure-function relationships for this emerging class of RNAs has been limited to secondary structures. Here, we report the 3-D atomistic structural study of epigenetic lncRNA, Braveheart (Bvht), and its complex with CNBP (Cellular Nucleic acid Binding Protein). Using small angle X-ray scattering (SAXS), we elucidate the ensemble of Bvht RNA conformations in solution, revealing that Bvht lncRNA has a well-defined, albeit flexible 3-D structure that is remodeled upon CNBP binding. Our study suggests that CNBP binding requires multiple domains of Bvht and the RHT/AGIL RNA motif. We show that RHT/AGIL, previously shown to interact with CNBP, contains a highly flexible loop surrounded by more ordered helices. As one of the largest RNA-only 3-D studies, the work lays the foundation for future structural studies of lncRNA-protein complexes.

Sensor-Enabled Alleviation of Product Inhibition in Chorismate Pyruvate-Lyase.

Product inhibition is a frequent bottleneck in industrial enzymes, and testing mutations to alleviate product inhibition via traditional methods remains challenging as many variants need to be tested against multiple substrate and product concentrations. Further, traditional screening methods are conducted in vitro, and resulting enzyme variants may perform differently in vivo in the context of whole-cell metabolism and regulation. In this study, we address these two problems by establishing a high-throughput screening method to alleviate product inhibition in an industrially relevant enzyme, chorismate pyruvate-lyase (UbiC). First, we engineered a highly specific, genetically encoded biosensor for 4-hydroxybenzoate (4HB) in an industrially relevant host, Pseudomonas putida KT2440. We subsequently applied the biosensor to detect the activity of a heterologously expressed UbiC that converts chorismate into 4HB and pyruvate. By using benzoate as a product surrogate that inhibits UbiC without activating the biosensor, we were able to efficiently create and screen a diversified library for UbiC variants with reduced product inhibition. Introduction of the improved UbiC enzyme variant into an experimental production strain for the industrial precursor cis,cis-muconic acid (muconate), enabled a >2-fold yield improvement for glucose to muconate conversion when the new UbiC variant was expressed from a plasmid and a 60% yield increase when the same UbiC variant was genomically integrated into the strain. Overall, this work demonstrates that by coupling a library of enzyme variants to whole-cell catalysis and biosensing, variants with reduced product inhibition can be identified, and that this improved enzyme can result in increased titers of a downstream molecule of interest.

Magnesium controls aptamer-expression platform switching in the SAM-I riboswitch.

Investigations of most riboswitches remain confined to the ligand-binding aptamer domain. However, during the riboswitch mediated transcription regulation process, the aptamer domain and the expression platform compete for a shared strand. If the expression platform dominates, an anti-terminator helix is formed, and the transcription process is active (ON state). When the aptamer dominates, transcription is terminated (OFF state). Here, we use an expression platform switching experimental assay and structure-based electrostatic simulations to investigate this ON-OFF transition of the full length SAM-I riboswitch and its magnesium concentration dependence. Interestingly, we find the ratio of the OFF population to the ON population to vary non-monotonically as magnesium concentration increases. Upon addition of magnesium, the aptamer domain pre-organizes, populating the OFF state, but only up to an intermediate magnesium concentration level. Higher magnesium concentration preferentially stabilizes the anti-terminator helix, populating the ON state, relatively destabilizing the OFF state. Magnesium mediated aptamer-expression platform domain closure explains this relative destabilization of the OFF state at higher magnesium concentration. Our study reveals the functional potential of magnesium in controlling transcription of its downstream genes and underscores the importance of a narrow concentration regime near the physiological magnesium concentration ranges, striking a balance between the OFF and ON states in bacterial gene regulation.

Functional and Structural Analysis of a Highly-Expressed Yersinia pestis Small RNA following Infection of Cultured Macrophages.

Non-coding small RNAs (sRNAs) are found in practically all bacterial genomes and play important roles in regulating gene expression to impact bacterial metabolism, growth, and virulence. We performed transcriptomics analysis to identify sRNAs that are differentially expressed in Yersinia pestis that invaded the human macrophage cell line THP-1, compared to pathogens that remained extracellular in the presence of host. Using ultra high-throughput sequencing, we identified 37 novel and 143 previously known sRNAs in Y. pestis. In particular, the sRNA Ysr170 was highly expressed in intracellular Yersinia and exhibited a log2 fold change ~3.6 higher levels compared to extracellular bacteria. We found that knock-down of Ysr170 expression attenuated infection efficiency in cell culture and growth rate in response to different stressors. In addition, we applied selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis to determine the secondary structure of Ysr170 and observed structural changes resulting from interactions with the aminoglycoside antibiotic gentamycin and the RNA chaperone Hfq. Interestingly, gentamicin stabilized helix 4 of Ysr170, which structurally resembles the native gentamicin 16S ribosomal binding site. Finally, we modeled the tertiary structure of Ysr170 binding to gentamycin using RNA motif modeling. Integration of these experimental and structural methods can provide further insight into the design of small molecules that can inhibit function of sRNAs required for pathogen virulence.

COOLAIR Antisense RNAs Form Evolutionarily Conserved Elaborate Secondary Structures.

There is considerable debate about the functionality of long non-coding RNAs (lncRNAs). Lack of sequence conservation has been used to argue against functional relevance. We investigated antisense lncRNAs, called COOLAIR, at the A. thaliana FLC locus and experimentally determined their secondary structure. The major COOLAIR variants are highly structured, organized by exon. The distally polyadenylated transcript has a complex multi-domain structure, altered by a single non-coding SNP defining a functionally distinct A. thaliana FLC haplotype. The A. thaliana COOLAIR secondary structure was used to predict COOLAIR exons in evolutionarily divergent Brassicaceae species. These predictions were validated through chemical probing and cloning. Despite the relatively low nucleotide sequence identity, the structures, including multi-helix junctions, show remarkable evolutionary conservation. In a number of places, the structure is conserved through covariation of a non-contiguous DNA sequence. This structural conservation supports a functional role for COOLAIR transcripts rather than, or in addition to, antisense transcription.

Tunable Riboregulator Switches for Post-transcriptional Control of Gene Expression.

Until recently, engineering strategies for altering gene expression have focused on transcription control using strong inducible promoters or one of several methods to knock down wasteful genes. Recently, synthetic riboregulators have been developed for translational regulation of gene expression. Here, we report a new modular synthetic riboregulator class that has the potential to finely tune protein expression and independently control the concentration of each enzyme in an engineered metabolic pathway. This development is important because the most straightforward approach to altering the flux through a particular metabolic step is to increase or decrease the concentration of the enzyme. Our design includes a cis-repressor at the 5' end of the mRNA that forms a stem-loop helix, occluding the ribosomal binding sequence and blocking translation. A trans-expressed activating-RNA frees the ribosomal-binding sequence, which turns on translation. The overall architecture of the riboregulators is designed using Watson-Crick base-pairing stability. We describe here a cis-repressor that can completely shut off translation of antibiotic-resistance reporters and a trans-activator that restores translation. We have established that it is possible to use these riboregulators to achieve translational control of gene expression over a wide dynamic range. We have also found that a targeting sequence can be modified to develop riboregulators that can, in principle, independently regulate translation of many genes. In a selection experiment, we demonstrated that by subtly altering the sequence of the trans-activator it is possible to alter the ratio of the repressed and activated states and to achieve intermediate translational control.

Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT.

Integration and calibration of molecular dynamics simulations with experimental data remain a challenging endeavor. We have developed a novel method to integrate chemical probing experiments with molecular simulations of RNA molecules by using a native structure-based model. Selective 2'-hydroxyl acylation by primer extension (SHAPE) characterizes the mobility of each residue in the RNA. Our method, SHAPE-FIT, automatically optimizes the potential parameters of the force field according to measured reactivities from SHAPE. The optimized parameter set allows simulations of dynamics highly consistent with SHAPE probing experiments. Such atomistic simulations, thoroughly grounded in experiment, can open a new window on RNA structure-function relations.

Tackling structures of long noncoding RNAs.

RNAs are important catalytic machines and regulators at every level of gene expression. A new class of RNAs has emerged called long non-coding RNAs, providing new insights into evolution, development and disease. Long non-coding RNAs (lncRNAs) predominantly found in higher eukaryotes, have been implicated in the regulation of transcription factors, chromatin-remodeling, hormone receptors and many other processes. The structural versatility of RNA allows it to perform various functions, ranging from precise protein recognition to catalysis and metabolite sensing. While major housekeeping RNA molecules have long been the focus of structural studies, lncRNAs remain the least characterized class, both structurally and functionally. Here, we review common methodologies used to tackle RNA structure, emphasizing their potential application to lncRNAs. When considering the complexity of lncRNAs and lack of knowledge of their structure, chemical probing appears to be an indispensable tool, with few restrictions in terms of size, quantity and heterogeneity of the RNA molecule. Probing is not constrained to in vitro analysis and can be adapted to high-throughput sequencing platforms. Significant efforts have been applied to develop new in vivo chemical probing reagents, new library construction protocols for sequencing platforms and improved RNA prediction software based on the experimental evidence.

3S: shotgun secondary structure determination of long non-coding RNAs.

Long non-coding RNAs (lncRNAs) have emerged as an important class of RNAs playing key roles in development, disease and epigenetics. Knowledge of lncRNA structure may be critical in understanding function for many lncRNA systems. Due to the enormous number of possible folds for these sequences, secondary structure determination presents a significant challenge, both experimentally and computationally. Here, we present a new strategy capable of determining the RNA secondary structure in the wet lab without significant reliance on computational predictions. First, we chemically probe the entire lncRNA. Next, using a shotgun approach, we divide the RNA into overlapping fragments and probe these fragments. We then compare probing profiles of fragments with the profiles of the full RNA and identify similarities. Sequence regions with profiles that are similar in the fragment and full-length transcript possess only base pairing partners within the fragment. Thus, by experimentally folding smaller and smaller fragments of the full RNA and probing these chemically, we are able to isolate modular sub-domains, dramatically reducing the number of possible folds. The method also eliminates the possibility of pseudoknots within a modular sub-domain. The 3S technique is ideally suited for lncRNAs because it is designed for long RNA sequences. The 3S-determined secondary structure of a specific lncRNA in one species (e.g., human) enables searches for instances of the same lncRNA in other species.

Rise of the RNA machines: exploring the structure of long non-coding RNAs.

Novel, profound and unexpected roles of long non-coding RNAs (lncRNAs) are emerging in critical aspects of gene regulation. Thousands of lncRNAs have been recently discovered in a wide range of mammalian systems, related to development, epigenetics, cancer, brain function and hereditary disease. The structural biology of these lncRNAs presents a brave new RNA world, which may contain a diverse zoo of new architectures and mechanisms. While structural studies of lncRNAs are in their infancy, we describe existing structural data for lncRNAs, as well as crystallographic studies of other RNA machines and their implications for lncRNAs. We also discuss the importance of dynamics in RNA machine mechanism. Determining commonalities between lncRNA systems will help elucidate the evolution and mechanistic role of lncRNAs in disease, creating a structural framework necessary to pursue lncRNA-based therapeutics.

The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch.

Riboswitch operation involves the complex interplay between the aptamer domain and the expression platform. During transcription, these two domains compete against each other for shared sequence. In this study, we explore the cooperative effects of ligand binding and Magnesium interactions in the SAM-I riboswitch in the context of aptamer collapse and anti-terminator formation. Overall, our studies show the apo-aptamer acts as (i) a pre-organized aptamer competent to bind ligand and undergo structural collapse and (ii) a conformation that is more accessible to anti-terminator formation. We show that both Mg(2+) ions and SAM are required for a collapse transition to occur. We then use competition between the aptamer and expression platform for shared sequence to characterize the stability of the collapsed aptamer. We find that SAM and Mg(2+) interactions in the aptamer are highly cooperative in maintaining switch polarity (i.e. aptamer 'off-state' versus anti-terminator 'on-state'). We further show that the aptamer off-state is preferentially stabilized by Mg(2+) and similar divalent ions. Furthermore, the functional switching assay was used to select for phosphorothioate interference, and identifies potential magnesium chelation sites while characterizing their coordinated role with SAM in aptamer stabilization. In addition, we find that Mg(2+) interactions with the apo-aptamer are required for the full formation of the anti-terminator structure, and that higher concentrations of Mg(2+) (>4 mM) shift the equilibrium toward the anti-terminator on-state even in the presence of SAM.

Sizing up long non-coding RNAs: do lncRNAs have secondary and tertiary structure?

Long noncoding RNAs (lncRNAs) play a key role in many important areas of epigenetics, stem cell biology, cancer, signaling and brain function. This emerging class of RNAs constitutes a large fraction of the transcriptome, with thousands of new lncRNAs reported each year. The molecular mechanisms of these RNAs are not well understood. Currently, very little structural data exist. We review the available lncRNA sequence and secondary structure data. Since almost no tertiary information is available for lncRNAs, we review crystallographic structures for other RNA systems and discuss the possibilities for lncRNAs in the context of existing constraints.

Magnesium fluctuations modulate RNA dynamics in the SAM-I riboswitch.

Experiments demonstrate that Mg(2+) is crucial for structure and function of RNA systems, yet the detailed molecular mechanism of Mg(2+) action on RNA is not well understood. We investigate the interplay between RNA and Mg(2+) at atomic resolution through ten 2-µs explicit solvent molecular dynamics simulations of the SAM-I riboswitch with varying ion concentrations. The structure, including three stemloops, is very stable on this time scale. Simulations reveal that outer-sphere coordinated Mg(2+) ions fluctuate on the same time scale as the RNA, and that their dynamics couple. Locally, Mg(2+) association affects RNA conformation through tertiary bridging interactions; globally, increasing Mg(2+) concentration slows RNA fluctuations. Outer-sphere Mg(2+) ions responsible for these effects account for 80% of Mg(2+) in our simulations. These ions are transiently bound to the RNA, maintaining interactions, but shuttled from site to site. Outer-sphere Mg(2+) are separated from the RNA by a single hydration shell, occupying a thin layer 3-5 Å from the RNA. Distribution functions reveal that outer-sphere Mg(2+) are positioned by electronegative atoms, hydration layers, and a preference for the major groove. Diffusion analysis suggests transient outer-sphere Mg(2+) dynamics are glassy. Since outer-sphere Mg(2+) ions account for most of the Mg(2+) in our simulations, these ions may change the paradigm of Mg(2+)-RNA interactions. Rather than a few inner-sphere ions anchoring the RNA structure surrounded by a continuum of diffuse ions, we observe a layer of outer-sphere coordinated Mg(2+) that is transiently bound but strongly coupled to the RNA.

Structural architecture of the human long non-coding RNA, steroid receptor RNA activator.

While functional roles of several long non-coding RNAs (lncRNAs) have been determined, the molecular mechanisms are not well understood. Here, we report the first experimentally derived secondary structure of a human lncRNA, the steroid receptor RNA activator (SRA), 0.87 kB in size. The SRA RNA is a non-coding RNA that coactivates several human sex hormone receptors and is strongly associated with breast cancer. Coding isoforms of SRA are also expressed to produce proteins, making the SRA gene a unique bifunctional system. Our experimental findings (SHAPE, in-line, DMS and RNase V1 probing) reveal that this lncRNA has a complex structural organization, consisting of four domains, with a variety of secondary structure elements. We examine the coevolution of the SRA gene at the RNA structure and protein structure levels using comparative sequence analysis across vertebrates. Rapid evolutionary stabilization of RNA structure, combined with frame-disrupting mutations in conserved regions, suggests that evolutionary pressure preserves the RNA structural core rather than its translational product. We perform similar experiments on alternatively spliced SRA isoforms to assess their structural features.

Excited states of ribosome translocation revealed through integrative molecular modeling.

The dynamic nature of biomolecules leads to significant challenges when characterizing the structural properties associated with function. While X-ray crystallography and imaging techniques (such as cryo-electron microscopy) can reveal the structural details of stable molecular complexes, strategies must be developed to characterize configurations that exhibit only marginal stability (such as intermediates) or configurations that do not correspond to minima on the energy landscape (such as transition-state ensembles). Here, we present a methodology (MDfit) that utilizes molecular dynamics simulations to generate configurations of excited states that are consistent with available biophysical and biochemical measurements. To demonstrate the approach, we present a sequence of configurations that are suggested to be associated with transfer RNA (tRNA) movement through the ribosome (translocation). The models were constructed by combining information from X-ray crystallography, cryo-electron microscopy, and biochemical data. These models provide a structural framework for translocation that may be further investigated experimentally and theoretically to determine the precise energetic character of each configuration and the transition dynamics between them.

Tertiary contacts control switching of the SAM-I riboswitch.

Riboswitches are non-coding RNAs that control gene expression by sensing small molecules through changes in secondary structure. While secondary structure and ligand interactions are thought to control switching, the exact mechanism of control is unknown. Using a novel two-piece assay that competes the anti-terminator against the aptamer, we directly monitor the process of switching. We find that the stabilization of key tertiary contacts controls both aptamer domain collapse and the switching of the SAM-I riboswitch from the aptamer to the expression platform conformation. Our experiments demonstrate that SAM binding induces structural alterations that indirectly stabilize the aptamer domain, preventing switching toward the expression platform conformer. These results, combined with a variety of structural probing experiments performed in this study, show that the collapse and stabilization of the aptamer domain are cooperative, relying on the sum of key tertiary contacts and the bimodal stability of the kink-turn motif for function. Here, ligand binding serves to shift the equilibrium of aptamer domain structures from a more open toward a more stable collapsed form by stabilizing tertiary interactions. Our data show that the thermodynamic landscape for riboswitch operation is finely balanced to allow large conformational rearrangements to be controlled by small molecule interactions.

Free state conformational sampling of the SAM-I riboswitch aptamer domain.

Riboswitches are highly structured elements residing in the 5' untranslated region of messenger RNAs that specifically bind cellular metabolites to alter gene expression. While there are many structures of ligand-bound riboswitches that reveal details of bimolecular recognition, their unliganded structures remain poorly characterized. Characterizing the molecular details of the unliganded state is crucial for understanding the riboswitch's mechanism of action because it is this state that actively interrogates the cellular environment and helps direct the regulatory outcome. To develop a detailed description of the ligand-free form of an S-adenosylmethionine binding riboswitch at the local and global levels, we have employed a series of biochemical, biophysical, and computational methods. Our data reveal that the ligand binding domain adopts an ensemble of states that minimizes the energy barrier between the free and bound states to establish an efficient decision making branchpoint in the regulatory process.