The Finite State Projection Approach for the Analysis of Stochastic Noise in Gene Networks

• Published in: IEEE transactions on automatic control 2008-01, Vol.53 (Special Issue), p.201-214
• Main Authors: Munsky, B., Khammash, M.
• Format: Article
• Language: English
• Subjects: Algorithm design and analysis; Algorithms; Biological system modeling; Chemical master equation; Chemicals; Computational biology; Computer simulation; Coronal mass ejection; Differential equations; gene regulatory network; Kinetic theory; Linear systems; Markov process; Mathematical analysis; Mathematical models; model reduction; Monte Carlo methods; Ordinary differential equations; Probability distribution; Projection; Reduction; stochastic process; Stochastic resonance; Studies; System theory; systems biology
• ISSN: 0018-9286; 1558-2523
• DOI: 10.1109/TAC.2007.911361

In order to capture important subcellular dynamics, researchers in computational biology have begun to turn to mesoscopic models in which molecular interactions at the gene level behave as discrete stochastic events. While the trajectories of such models cannot be described with deterministic expressions, the probability distributions of these trajectories can be described by the set of linear ordinary differential equations known as the chemical master equation (CME). Until recently, it has been believed that the CME could only be solved analytically in the most trivial of problems, and the CME has been analyzed almost exclusively with kinetic Monte Carlo (KMC) algorithms. However, concepts from linear systems theory have enabled the finite state projection (FSP) approach and have significantly enhanced our ability to solve the CME without resorting to KMC simulations. In this paper, we review the FSP approach and introduce a variety of systems-theory-based modifications and enhancements to the FSP algorithm. Notions such as observability, controllability, and minimal realizations enable large reductions and increase efficiency with little to no loss in accuracy. Model reduction techniques based upon linear perturbation theory allow for the systematic projection of multiple time-scale dynamics onto a slowly varying manifold of much smaller dimension. We also present a powerful new reduction approach, in which we perform computations on a small subset of configuration grid points and then interpolate to find the distribution on the full set. The power of the FSP and its various reduction approaches is illustrated on few important models of genetic regulatory networks.