Chemical Master Equation Empirical Moment Closure

The numerical solution of the Chemical Master Equation (CME) governing gene regulatory networks and cell signaling processes remains a challenging task due to its complexity, exponentially growing with the number of species involved. When considering separated representations of the probability distribution function within the Proper Generalized Decomposition-PGD-frame-work the complexity of the CME grows only linearly with the number of state space dimensions. In order to speed up calculations moment-based descriptions are usually preferred, however these descriptions involve the necessity of using closure relations whose impact on the calculated solution is most of time unpredictable. In this work we propose an empirical closure, fitted from the solution of the chemical master equation, the last solved within the PGD framework.

[1]  Farren J. Isaacs,et al.  Computational studies of gene regulatory networks: in numero molecular biology , 2001, Nature Reviews Genetics.

[2]  Ankit Gupta,et al.  A Scalable Computational Framework for Establishing Long-Term Behavior of Stochastic Reaction Networks , 2013, PLoS Comput. Biol..

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

[4]  Wolf Verena,et al.  Reconstruction of Multimodal Distributions for Hybrid Moment-based Chemical Kinetics , 2014 .

[5]  C. J. Burden,et al.  A solver for the stochastic master equation applied to gene regulatory networks , 2007 .

[6]  Kwang-Hyun Cho,et al.  Modelling the dynamics of signalling pathways. , 2008, Essays in biochemistry.

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

[8]  Francisco Chinesta,et al.  Consistent closure schemes for statistical models of anisotropic fluids , 2008 .

[9]  Elías Cueto,et al.  Reduction of the chemical master equation for gene regulatory networks using proper generalized decompositions , 2012, International journal for numerical methods in biomedical engineering.

[10]  Francisco Chinesta,et al.  A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids , 2006 .

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

[12]  Francisco Chinesta,et al.  A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modelling of complex fluids - Part II: Transient simulation using space-time separated representations , 2007 .

[13]  Morgan Magnin,et al.  Kinetic Theory Modeling and Efficient Numerical Simulation of Gene Regulatory Networks Based on Qualitative Descriptions , 2015, Entropy.

[14]  D. Wilkinson,et al.  Moment closure approximations for stochastic kinetic models with rational rate laws. , 2011, Mathematical biosciences.

[15]  D. Gillespie Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[16]  C. Gillespie Moment-closure approximations for mass-action models. , 2009, IET systems biology.