Numerical Methods for Stochastic Modeling of Genes and Proteins

Abstract: Stochastic models of biochemical reaction networks are used for understanding the properties of molecular regulatory circuits in living cells. The state of the cell is defined by the number of copies of each molecular species in the model. The chemical master equation (CME) governs the time evolution of the the probability density function of the often high-dimensional state space. The CME is approximated by a partial differential equation (PDE), the Fokker-Planck equation and solved numerically. Direct solution of the CME rapidly becomes computationally expensive for increasingly complex biological models, since the state space grows exponentially with the number of dimensions. Adaptive numerical methods can be applied in time and space in the PDE framework, and error estimates of the approximate solutions are derived. A method for splitting the CME operator in order to apply the PDE approximation in a subspace of the state space is also developed. The performance is compared to the most widely spread alternative computational method.