On the reachable set of the controlled gene expression system

In this paper we investigate the reachable set of a standard stochastic model of controlled gene expression. Specifically, we explore what values of the protein mean and variance are achievable using the available external input, that is, the mRNA production rate. We proceed by constructing invariant sets in two-dimensional projections of the state space. We then use these sets to construct an outer approximation of the reachable region for the protein mean and variance. This can be computed solving a one-dimensional optimization problem and is tight enough to show that it is not possible, using such control input, to arbitrarily reduce the variance while maintaining a high mean. The lower bound on the variance derived with our approach turns out to be much higher than the one available in the literature for the case when both the mRNA production and degradation rate are controlled.

[1]  M. Bennett,et al.  A fast, robust, and tunable synthetic gene oscillator , 2008, Nature.

[2]  Nancy A. Lynch,et al.  Hybrid I/O automata , 1995, Inf. Comput..

[3]  Ian M. Mitchell,et al.  A Toolbox of Hamilton-Jacobi Solvers for Analysis of Nondeterministic Continuous and Hybrid Systems , 2005, HSCC.

[4]  Jean-Pierre Aubin,et al.  Characterization of Stochastic Viability of Any Nonsmooth Set Involving Its Generalized Contingent Curvature , 2003 .

[5]  F. Fages,et al.  Long-term model predictive control of gene expression at the population and single-cell levels , 2012, Proceedings of the National Academy of Sciences.

[6]  M. Elowitz,et al.  A synthetic oscillatory network of transcriptional regulators , 2000, Nature.

[7]  Mario di Bernardo,et al.  Analysis, design and implementation of a novel scheme for in-vivo control of synthetic gene regulatory networks , 2011, Autom..

[8]  Jesús Picó,et al.  Control of protein concentrations in heterogeneous cell populations , 2013, 2013 European Control Conference (ECC).

[9]  A. Kurzhanski,et al.  Ellipsoidal Calculus for Estimation and Control , 1996 .

[10]  M. Elowitz,et al.  Functional roles for noise in genetic circuits , 2010, Nature.

[11]  G. Vinnicombe,et al.  Fundamental limits on the suppression of molecular fluctuations , 2010, Nature.

[12]  Oded Maler,et al.  Systematic Simulation Using Sensitivity Analysis , 2007, HSCC.

[13]  John Lygeros,et al.  Grey-box techniques for the identification of a controlled gene expression model , 2014, 2014 European Control Conference (ECC).

[14]  J. Lygeros,et al.  Moment-based inference predicts bimodality in transient gene expression , 2012, Proceedings of the National Academy of Sciences.

[15]  J. Paulsson.  Summing up the noise in gene networks , 2004, Nature.

[16]  Jean-Pierre Aubin,et al.  Viability theory , 1991 .

[17]  Ahmad S. Khalil,et al.  Synthetic biology: applications come of age , 2010, Nature Reviews Genetics.

[18]  Jeffrey J. Tabor,et al.  Characterizing bacterial gene circuit dynamics with optically programmed gene expression signals , 2014, Nature Methods.

[19]  Glenn Vinnicombe,et al.  Single event molecular signalling for estimation and control , 2013, 2013 European Control Conference (ECC).

[20]  R. Durrett.  Probability: Theory and Examples , 1993 .

[21]  Franco Blanchini,et al.  Set invariance in control , 1999, Autom..

[22]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[23]  Corentin Briat,et al.  Computer control of gene expression: Robust setpoint tracking of protein mean and variance using integral feedback , 2012, 2012 IEEE 51st IEEE Conference on Decision and Control (CDC).

[24]  B. Krogh,et al.  Hyperplane method for reachable state estimation for linear time-invariant systems , 1991 .

[25]  Uwe Helmke,et al.  Global identifiability of a simple linear model for gene expression analysis , 2013, 52nd IEEE Conference on Decision and Control.

[26]  D. Pincus,et al.  In silico feedback for in vivo regulation of a gene expression circuit , 2011, Nature Biotechnology.

[27]  Kumpati S. Narendra,et al.  Reachable Sets for Linear Dynamical Systems , 1971, Inf. Control..

[28]  Brian Munsky,et al.  Listening to the noise: random fluctuations reveal gene network parameters , 2009, Molecular systems biology.

[29]  Stefan Engblom,et al.  Computing the moments of high dimensional solutions of the master equation , 2006, Appl. Math. Comput..