Time-dependent product-form Poisson distributions for reaction networks with higher order complexes

J Math Biol. 2020 May;80(6):1919-1951. doi: 10.1007/s00285-020-01485-y. Epub 2020 Mar 24.

Abstract

It is well known that stochastically modeled reaction networks that are complex balanced admit a stationary distribution that is a product of Poisson distributions. In this paper, we consider the following related question: supposing that the initial distribution of a stochastically modeled reaction network is a product of Poissons, under what conditions will the distribution remain a product of Poissons for all time? By drawing inspiration from Crispin Gardiner's "Poisson representation" for the solution to the chemical master equation, we provide a necessary and sufficient condition for such a product-form distribution to hold for all time. Interestingly, the condition is a dynamical "complex-balancing" for only those complexes that have multiplicity greater than or equal to two (i.e. the higher order complexes that yield non-linear terms to the dynamics). We term this new condition the "dynamical and restricted complex balance" condition (DR for short).

Keywords: Complex balancing; Deficiency; Poisson distribution; Reaction networks; Stochastic processes.

Publication types

  • Research Support, Non-U.S. Gov't
  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Gene Regulatory Networks
  • Kinetics
  • Linear Models
  • Markov Chains
  • Mathematical Concepts
  • Metabolic Networks and Pathways
  • Models, Biological*
  • Models, Chemical
  • Nonlinear Dynamics
  • Poisson Distribution
  • Signal Transduction
  • Stochastic Processes
  • Systems Biology / statistics & numerical data*