Semi-Quantitative Abstraction and Analysis of Chemical Reaction Networks

Analysis of large continuous-time stochastic systems is a computationally intensive task. In this work we focus on population models arising from chemical reaction networks (CRNs), which play a fundamental role in analysis and design of biochemical systems. Many relevant CRNs are particularly challenging for existing techniques due to complex dynamics including stochasticity, stiffness or multimodal population distributions. We propose a novel approach allowing not only to predict, but also to explain both the transient and steady-state behaviour. It focuses on qualitative description of the behaviour and aims at quantitative precision only in orders of magnitude. Firstly, we abstract the CRN into a compact model preserving rough timing information, distinguishing only signifcinatly different populations, but capturing relevant sequences of behaviour. Secondly, we approximately analyse the most probable temporal behaviours of the model through most probable transitions. As demonstrated on complex CRNs from literature, our approach reproduces the known results, but in contrast to the state-of-the-art methods, it runs with virtually no computational cost and thus offers unprecedented~scalability.

[1]  Peter Buchholz Exact Performance Equivalence: An Equivalence Relation for Stochastic Automata , 1999, Theor. Comput. Sci..

[2]  Joost-Pieter Katoen,et al.  Three-Valued Abstraction for Continuous-Time Markov Chains , 2007, CAV.

[3]  Thomas A. Henzinger,et al.  Model Checking Gene Regulatory Networks , 2015, TACAS.

[4]  Fabian J Theis,et al.  Method of conditional moments (MCM) for the Chemical Master Equation , 2013, Journal of Mathematical Biology.

[5]  Marta Z. Kwiatkowska,et al.  Probabilistic model checking of complex biological pathways , 2008, Theor. Comput. Sci..

[6]  Per Lötstedt,et al.  Adaptive solution of the master equation in low dimensions , 2009 .

[7]  Luca Bortolussi,et al.  Fluid Model Checking , 2012, CONCUR.

[8]  Luca Cardelli,et al.  Maximal aggregation of polynomial dynamical systems , 2017, Proceedings of the National Academy of Sciences.

[9]  Heinz Koeppl,et al.  Jump-Diffusion Approximation of Stochastic Reaction Dynamics: Error Bounds and Algorithms , 2014, Multiscale Model. Simul..

[10]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[11]  Yiannis Kaznessis,et al.  Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions. , 2005, The Journal of chemical physics.

[12]  David Eisenstat,et al.  The computational power of population protocols , 2006, Distributed Computing.

[13]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[14]  Christel Baier,et al.  Principles of model checking , 2008 .

[15]  Luca Cardelli,et al.  Design and analysis of DNA strand displacement devices using probabilistic model checking , 2012, Journal of The Royal Society Interface.

[16]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[17]  Thomas A. Henzinger,et al.  Fast Adaptive Uniformization of the Chemical Master Equation , 2009 .

[18]  G. Seelig,et al.  DNA as a universal substrate for chemical kinetics , 2010, Proceedings of the National Academy of Sciences.

[19]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[20]  Alessandro Abate,et al.  Adaptive and Sequential Gridding Procedures for the Abstraction and Verification of Stochastic Processes , 2013, SIAM J. Appl. Dyn. Syst..

[21]  Lubos Brim,et al.  Adaptive Aggregation of Markov Chains: Quantitative Analysis of Chemical Reaction Networks , 2015, CAV.

[22]  Thomas A. Henzinger,et al.  Sliding Window Abstraction for Infinite Markov Chains , 2009, CAV.

[23]  S. Ethier,et al.  Markov Processes: Characterization and Convergence , 2005 .

[24]  Robert H. Singer,et al.  Transcription of functionally related constitutive genes is not coordinated , 2010, Nature Structural &Molecular Biology.

[25]  François Laviolette,et al.  Approximate Analysis of Probabilistic Processes: Logic, Simulation and Games , 2008, 2008 Fifth International Conference on Quantitative Evaluation of Systems.

[26]  Kim G. Larsen,et al.  Bisimulation through Probabilistic Testing , 1991, Inf. Comput..

[27]  Zhen Zhang,et al.  Utilizing stochastic model checking to analyze genetic circuits , 2012, 2012 IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology (CIBCB).

[28]  Richard M. Karp,et al.  Parallel Program Schemata , 1969, J. Comput. Syst. Sci..

[29]  Jingwei Zhang,et al.  Adaptive aggregation method for the chemical master equation , 2008, 2008 8th IEEE International Conference on BioInformatics and BioEngineering.

[30]  Tadao Murata,et al.  Petri nets: Properties, analysis and applications , 1989, Proc. IEEE.

[31]  S. Bhat,et al.  Modeling and analysis of mass-action kinetics , 2009, IEEE Control Systems.

[32]  Luca Cardelli,et al.  Two-domain DNA strand displacement , 2010, Mathematical Structures in Computer Science.

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

[34]  Thomas A. Henzinger,et al.  Fast Adaptive Uniformization of the Chemical Master Equation , 2009, 2009 International Workshop on High Performance Computational Systems Biology.

[35]  J. Goutsias Quasiequilibrium approximation of fast reaction kinetics in stochastic biochemical systems. , 2005, The Journal of chemical physics.

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

[37]  E. Cox,et al.  Real-Time Kinetics of Gene Activity in Individual Bacteria , 2005, Cell.

[38]  Joost-Pieter Katoen,et al.  Robust PCTL model checking , 2012, HSCC '12.

[39]  Luca Cardelli,et al.  A Stochastic Hybrid Approximation for Chemical Kinetics Based on the Linear Noise Approximation , 2016, CMSB.

[40]  Thomas A. Henzinger,et al.  Hybrid numerical solution of the chemical master equation , 2010, CMSB '10.

[41]  Ankit Gupta,et al.  Adaptive hybrid simulations for multiscale stochastic reaction networks. , 2014, The Journal of chemical physics.

[42]  L. You,et al.  Stochastic vs. deterministic modeling of intracellular viral kinetics. , 2002, Journal of theoretical biology.

[43]  Alessandro Abate,et al.  Precise Approximations of the Probability Distribution of a Markov Process in Time: An Application to Probabilistic Invariance , 2014, TACAS.

[44]  Chris Thachuk,et al.  Probabilistic Model Checking for Biology , 2014, Software Systems Safety.

[45]  Jan Kvret'insk'y,et al.  Semi-Quantitative Abstraction and Analysis of Chemical Reaction Networks , 2019, CAV.

[46]  Joost-Pieter Katoen,et al.  Approximate Model Checking of Stochastic Hybrid Systems , 2010, Eur. J. Control.