Ribozyme-based insulator parts buffer synthetic circuits from genetic context

Synthetic genetic programs are built from circuits that integrate sensors and implement temporal control of gene expression. Transcriptional circuits are layered by using promoters to carry the signal between circuits. In other words, the output promoter of one circuit serves as the input promoter to the next. Thus, connecting circuits requires physically connecting a promoter to the next circuit. We show that the sequence at the junction between the input promoter and circuit can affect the input-output response (transfer function) of the circuit. A library of putative sequences that might reduce (or buffer) such context effects, which we refer to as 'insulator parts', is screened in Escherichia coli. We find that ribozymes that cleave the 5′ untranslated region (5′-UTR) of the mRNA are effective insulators. They generate quantitatively identical transfer functions, irrespective of the identity of the input promoter. When these insulators are used to join synthetic gene circuits, the behavior of layered circuits can be predicted using a mathematical model. The inclusion of insulators will be critical in reliably permuting circuits to build different programs.Life Technologies, Inc.United States. Defense Advanced Research Projects Agency (DARPA CLIO N66001-12-C-4018)United States. Office of Naval Research (N00014-10-1-0245)National Science Foundation (U.S.) (CCF-0943385)National Institutes of Health (U.S.) (AI067699)National Science Foundation (U.S.). Synthetic Biology Engineering Research Center (SynBERC, SA5284-11210

Global parameter identification of stochastic reaction networks from single trajectories

We consider the problem of inferring the unknown parameters of a stochastic biochemical network model from a single measured time-course of the concentration of some of the involved species. Such measurements are available, e.g., from live-cell fluorescence microscopy in image-based systems biology. In addition, fluctuation time-courses from, e.g., fluorescence correlation spectroscopy provide additional information about the system dynamics that can be used to more robustly infer parameters than when considering only mean concentrations. Estimating model parameters from a single experimental trajectory enables single-cell measurements and quantification of cell--cell variability. We propose a novel combination of an adaptive Monte Carlo sampler, called Gaussian Adaptation, and efficient exact stochastic simulation algorithms that allows parameter identification from single stochastic trajectories. We benchmark the proposed method on a linear and a non-linear reaction network at steady state and during transient phases. In addition, we demonstrate that the present method also provides an ellipsoidal volume estimate of the viable part of parameter space and is able to estimate the physical volume of the compartment in which the observed reactions take place.Comment: Article in print as a book chapter in Springer's "Advances in Systems Biology

Bridging Time Scales in Cellular Decision Making with a Stochastic Bistable Switch

Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms. Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.Comment: 14 pages, 4 figure

Complex and unexpected dynamics in simple genetic regulatory networks

Peer reviewedPublisher PD

Cell-Cycle Dependence of Transcription Dominates Noise in Gene Expression

The large variability in mRNA and protein levels found from both static and dynamic measurements in single cells has been largely attributed to random periods of transcription, often occurring in bursts. The cell cycle has a pronounced global role in affecting transcriptional and translational output, but how this influences transcriptional statistics from noisy promoters is unknown and generally ignored by current stochastic models. Here we show that variable transcription from the synthetic tetO promoter in S. cerevisiae is dominated by its dependence on the cell cycle. Real-time measurements of fluorescent protein at high expression levels indicate tetO promoters increase transcription rate ~2-fold in S/G2/M similar to constitutive genes. At low expression levels, where tetO promoters are thought to generate infrequent bursts of transcription, we observe random pulses of expression restricted to S/G2/M, which are correlated between homologous promoters present in the same cell. The analysis of static, single-cell mRNA measurements at different points along the cell cycle corroborates these findings. Our results demonstrate that highly variable mRNA distributions in yeast are not solely the result of randomly switching between periods of active and inactive gene expression, but instead largely driven by differences in transcriptional activity between G1 and S/G2/M.GM095733BBBE 103316MIT Startup Fun

Numerical Integration of the Master Equation in Some Models of Stochastic Epidemiology

The processes by which disease spreads in a population of individuals are inherently stochastic. The master equation has proven to be a useful tool for modeling such processes. Unfortunately, solving the master equation analytically is possible only in limited cases (e.g., when the model is linear), and thus numerical procedures or approximation methods must be employed. Available approximation methods, such as the system size expansion method of van Kampen, may fail to provide reliable solutions, whereas current numerical approaches can induce appreciable computational cost. In this paper, we propose a new numerical technique for solving the master equation. Our method is based on a more informative stochastic process than the population process commonly used in the literature. By exploiting the structure of the master equation governing this process, we develop a novel technique for calculating the exact solution of the master equation – up to a desired precision – in certain models of stochastic epidemiology. We demonstrate the potential of our method by solving the master equation associated with the stochastic SIR epidemic model. MATLAB software that implements the methods discussed in this paper is freely available as Supporting Information S1

Analysis of Stochastic Strategies in Bacterial Competence: A Master Equation Approach

Competence is a transiently differentiated state that certain bacterial cells reach when faced with a stressful environment. Entrance into competence can be attributed to the excitability of the dynamics governing the genetic circuit that regulates this cellular behavior. Like many biological behaviors, entrance into competence is a stochastic event. In this case cellular noise is responsible for driving the cell from a vegetative state into competence and back. In this work we present a novel numerical method for the analysis of stochastic biochemical events and use it to study the excitable dynamics responsible for competence in Bacillus subtilis. Starting with a Finite State Projection (FSP) solution of the chemical master equation (CME), we develop efficient numerical tools for accurately computing competence probability. Additionally, we propose a new approach for the sensitivity analysis of stochastic events and utilize it to elucidate the robustness properties of the competence regulatory genetic circuit. We also propose and implement a numerical method to calculate the expected time it takes a cell to return from competence. Although this study is focused on an example of cell-differentiation in Bacillus subtilis, our approach can be applied to a wide range of stochastic phenomena in biological systems

Linear mapping approximation of gene regulatory networks with stochastic dynamics

The intractability of most stochastic models of gene regulatory networks (GRNs) limits their utility. Here, the authors present a linear-mapping approximation mapping models onto simpler ones, giving approximate but accurate analytic or semi- analytic solutions for a wide range of model GRNs