A discrete spectral method for the chemical master equation

As an equivalent formulation of the Markov-assumption of stochastic processes, the master equation of chemical reactions is an accurate description of general systems in chemistry. ForD reacting species this is a differential-difference equation in D dimensions, exactly soluble for very simple systems only. We present and analyze a novel solution strategy based upon a Galerkin spectral method with an inherent natural adaptivity and a very favorable choice of basis functions. The method is demonstrated by the numerical solution of two model problems followed by two more realistic systems taken from molecular biology. It is shown that the method remains effective and accurate, providing a viable alternative to other solution methods when the dimensionality is not too high.

[1]  E. Dynkin.  Functionals of Markov processes , 1965 .

[2]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[3]  F. Olver.  Asymptotics and Special Functions , 1974 .

[4]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[5]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[6]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[7]  D. Sherrington.  Stochastic Processes in Physics and Chemistry , 1983 .

[8]  Gene H. Golub,et al.  Matrix computations , 1983 .

[9]  Douglas Quadling,et al.  Ramanujan's Notebooks , 1986 .

[10]  M. Haviv.  Aggregation/disaggregation methods for computing the stationary distribution of a Markov chain , 1987 .

[11]  Peter Deuflhard,et al.  Computational treatment of polyreaction kinetics by orthogonal polynomials of a discrete variable , 1989, IMPACT Comput. Sci. Eng..

[12]  W. J. Anderson.  Continuous-Time Markov Chains , 1991 .

[13]  MIcHAEL WULKoW,et al.  Adaptive treatment of polyreactions in weighted sequence spaces , 1992, IMPACT Comput. Sci. Eng..

[14]  D. Gillespie.  A rigorous derivation of the chemical master equation , 1992 .

[15]  P. Guptasarma,et al.  Does replication-induced transcription regulate synthesis of the myriad low copy number proteins of Escherichia coli? , 1995, BioEssays : news and reviews in molecular, cellular and developmental biology.

[16]  M. Wulkow,et al.  The simulation of molecular weight distributions in polyreaction kinetics by discrete Galerkin methods , 1996 .

[17]  Rene F. Swarttouw,et al.  The Askey-scheme of hypergeometric orthogonal polynomials and its q-analogue Report Fac , 1996, math/9602214.

[18]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[19]  B. Berndt.  Ramanujan’s Notebooks: Part V , 1997 .

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

[21]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[22]  Jie Shen,et al.  Stable and Efficient Spectral Methods in Unbounded Domains Using Laguerre Functions , 2000, SIAM J. Numer. Anal..

[23]  T. Mark Dunster.  Uniform Asymptotic Expansions for Charlier Polynomials , 2001, J. Approx. Theory.

[24]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[25]  J. Elf,et al.  Problems of high dimension in molecular biology , 2003 .

[26]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[27]  Hermann G. Matthies,et al.  Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations , 2005 .

[28]  Yang Cao,et al.  Multiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems , 2005 .

[29]  Stefan Engblom.  Gaussian quadratures with respect to discrete measures , 2006 .

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

[31]  Per Lötstedt,et al.  Conservative solution of the Fokker–Planck equation for stochastic chemical reactions , 2006 .

[32]  Per Lötstedt,et al.  Dimensional Reduction of the Fokker-Planck Equation for Stochastic Chemical Reactions , 2006, Multiscale Model. Simul..

[33]  Adam P Arkin,et al.  Deviant effects in molecular reaction pathways , 2006, Nature Biotechnology.

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

[35]  Peter Deuflhard,et al.  Adaptive Discrete Galerkin Methods Applied to the Chemical Master Equation , 2008, SIAM J. Sci. Comput..

[36]  Sigal Gottlieb,et al.  Spectral Methods , 2019, Numerical Methods for Diffusion Phenomena in Building Physics.

[37]  Per Lötstedt,et al.  Fokker–Planck approximation of the master equation in molecular biology , 2009 .