Computation of biochemical pathway fluctuations beyond the linear noise approximation using iNA

The linear noise approximation is commonly used to obtain intrinsic noise statistics for biochemical networks. These estimates are accurate for networks with large numbers of molecules. However it is well known that many biochemical networks are characterized by at least one species with a small number of molecules. We here describe version 0.3 of the software intrinsic Noise Analyzer (iNA) which allows for accurate computation of noise statistics over wide ranges of molecule numbers. This is achieved by calculating the next order corrections to the linear noise approximation's estimates of variance and covariance of concentration fluctuations. The efficiency of the methods is significantly improved by automated just-in-time compilation using the LLVM framework leading to a fluctuation analysis which typically outperforms that obtained by means of exact stochastic simulations. iNA is hence particularly well suited for the needs of the computational biology community.

[1]  Philipp Thomas,et al.  How accurate are the nonlinear chemical Fokker-Planck and chemical Langevin equations? , 2011, The Journal of chemical physics.

[2]  Johan Paulsson,et al.  Models of stochastic gene expression , 2005 .

[3]  Ramon Grima,et al.  A study of the accuracy of moment-closure approximations for stochastic chemical kinetics. , 2012, The Journal of chemical physics.

[4]  Ramon Grima,et al.  Discreteness-induced concentration inversion in mesoscopic chemical systems , 2012, Nature Communications.

[5]  Ertugrul M. Ozbudak,et al.  Regulation of noise in the expression of a single gene , 2002, Nature Genetics.

[6]  W. Ebeling.  Stochastic Processes in Physics and Chemistry , 1995 .

[7]  D. Frishman,et al.  Protein abundance profiling of the Escherichia coli cytosol , 2008, BMC Genomics.

[8]  Vahid Shahrezaei,et al.  Analytical distributions for stochastic gene expression , 2008, Proceedings of the National Academy of Sciences.

[9]  Van Kampen,et al.  The Expansion of the Master Equation , 2007 .

[10]  R. Grima,et al.  Intrinsic Noise Analyzer: A Software Package for the Exploration of Stochastic Biochemical Kinetics Using the System Size Expansion , 2012, PloS one.

[11]  Richard Kreckel,et al.  Introduction to the GiNaC Framework for Symbolic Computation within the C++ Programming Language , 2000, J. Symb. Comput..

[12]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[13]  Philipp Thomas,et al.  Stochastic theory of large-scale enzyme-reaction networks: finite copy number corrections to rate equation models. , 2010, The Journal of chemical physics.

[14]  Vikram S. Adve,et al.  LLVM: a compilation framework for lifelong program analysis & transformation , 2004, International Symposium on Code Generation and Optimization, 2004. CGO 2004..

[15]  L. Petzold.  Automatic Selection of Methods for Solving Stiff and Nonstiff Systems of Ordinary Differential Equations , 1983 .

[16]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[17]  Herbert M. Sauro,et al.  Conservation analysis of large biochemical networks , 2006, Bioinform..