Control mechanisms for stochastic biochemical systems via computation of reachable sets

Controlling the behaviour of cells by rationally guiding molecular processes is an overarching aim of much of synthetic biology. Molecular processes, however, are notoriously noisy and frequently non-linear. We present an approach to studying the impact of control measures on motifs of molecular interactions, that addresses the problems faced in biological systems: stochasticity, parameter uncertainty, and non-linearity. We show that our reachability analysis formalism can describe the potential behaviour of biological (naturally evolved as well as engineered) systems, and provides a set of bounds on their dynamics at the level of population statistics: for example, we can obtain the possible ranges of means and variances of mRNA and protein expression levels, even in the presence of uncertainty about model parameters.

[1]  Matthias Althoff,et al.  Reachable set computation for uncertain time-varying linear systems , 2011, HSCC '11.

[2]  Ian M. Mitchell,et al.  Level set methods for computing reachable sets of systems with differential algebraic equation dynamics , 2008, 2008 American Control Conference.

[3]  Michael P H Stumpf,et al.  Multivariate moment closure techniques for stochastic kinetic models. , 2015, The Journal of chemical physics.

[4]  M. Khammash,et al.  The finite state projection algorithm for the solution of the chemical master equation. , 2006, The Journal of chemical physics.

[5]  J. Aubin,et al.  Differential inclusions set-valued maps and viability theory , 1984 .

[6]  Gábor Szederkényi,et al.  Reachability of nonlinear fed‐batch fermentation processes , 2002 .

[7]  Xenofon D. Koutsoukos,et al.  Verification of Biochemical Processes Using Stochastic Hybrid Systems , 2007, 2007 IEEE 22nd International Symposium on Intelligent Control.

[8]  Komei Fukuda,et al.  From the zonotope construction to the Minkowski addition of convex polytopes , 2004, J. Symb. Comput..

[9]  Angelika Bayer,et al.  Ellipsoidal Calculus For Estimation And Control , 2016 .

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

[11]  Antonis Papachristodoulou,et al.  On validation and invalidation of biological models , 2009, BMC Bioinformatics.

[12]  C. Belta,et al.  Analysis of lactose metabolism in E. Coli using reachability analysis of hybrid systems. , 2007, IET systems biology.

[13]  Antoine Girard,et al.  Reachability of Uncertain Linear Systems Using Zonotopes , 2005, HSCC.

[14]  Michael P H Stumpf,et al.  Topological sensitivity analysis for systems biology , 2014, Proceedings of the National Academy of Sciences.

[15]  A. Papachristodoulou,et al.  Stochastic reachability analysis in complex biological networks , 2006, 2006 American Control Conference.

[16]  Michael P H Stumpf,et al.  A general moment expansion method for stochastic kinetic models. , 2013, The Journal of chemical physics.

[17]  Thao Dang,et al.  Analysis of parametric biological models with non-linear dynamics , 2012, HSB.

[18]  Marco Bozzano,et al.  Verifying Industrial Hybrid Systems with MathSAT , 2005, BMC@CAV.

[19]  Lubos Brim,et al.  Reachability in Biochemical Dynamical Systems by Quantitative Discrete Approximation , 2011, CompMod.

[20]  Matthias Althoff,et al.  Reachability Analysis of Nonlinear Differential-Algebraic Systems , 2014, IEEE Transactions on Automatic Control.

[21]  Michael P. H. Stumpf,et al.  Nonidentifiability of the Source of Intrinsic Noise in Gene Expression from Single-Burst Data , 2008, PLoS Comput. Biol..

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

[23]  Oded Maler,et al.  Computing reachable states for nonlinear biological models , 2009, Theor. Comput. Sci..

[24]  Paul D. W. Kirk,et al.  MEANS: python package for Moment Expansion Approximation, iNference and Simulation , 2016, Bioinform..

[25]  Mario di Bernardo,et al.  In-Vivo Real-Time Control of Protein Expression from Endogenous and Synthetic Gene Networks , 2014, PLoS Comput. Biol..

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

[27]  Bruce H. Krogh,et al.  Reachability Analysis of Large-Scale Affine Systems Using Low-Dimensional Polytopes , 2006, HSCC.

[28]  Sarah Filippi,et al.  Robustness of MEK-ERK Dynamics and Origins of Cell-to-Cell Variability in MAPK Signaling , 2016, Cell reports.

[29]  Francesca Parise,et al.  On the reachable set of the controlled gene expression system , 2014, 53rd IEEE Conference on Decision and Control.

[30]  B. Krogh,et al.  Reachability analysis of nonlinear systems using trajectory piecewise linearized models , 2006, 2006 American Control Conference.

[31]  Ian M. Mitchell,et al.  Reachability Analysis Using Polygonal Projections , 1999, HSCC.

[32]  Antoine Girard,et al.  Reachability Analysis of Nonlinear Systems Using Conservative Approximation , 2003, HSCC.

[33]  ICHAEL,et al.  Information Processing by Simple Molecular Motifs and Susceptibility to Noise , 2015, bioRxiv.

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

[35]  Kamil Erguler,et al.  Practical limits for reverse engineering of dynamical systems: a statistical analysis of sensitivity and parameter inferability in systems biology models. , 2011, Molecular bioSystems.

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