全部 标题 作者
关键词 摘要

OALib Journal期刊
ISSN: 2333-9721
费用:99美元

查看量下载量

相关文章

更多...
PLOS ONE  2012 

Intrinsic Noise Analyzer: A Software Package for the Exploration of Stochastic Biochemical Kinetics Using the System Size Expansion

DOI: 10.1371/journal.pone.0038518

Full-Text   Cite this paper   Add to My Lib

Abstract:

The accepted stochastic descriptions of biochemical dynamics under well-mixed conditions are given by the Chemical Master Equation and the Stochastic Simulation Algorithm, which are equivalent. The latter is a Monte-Carlo method, which, despite enjoying broad availability in a large number of existing software packages, is computationally expensive due to the huge amounts of ensemble averaging required for obtaining accurate statistical information. The former is a set of coupled differential-difference equations for the probability of the system being in any one of the possible mesoscopic states; these equations are typically computationally intractable because of the inherently large state space. Here we introduce the software package intrinsic Noise Analyzer (iNA), which allows for systematic analysis of stochastic biochemical kinetics by means of van Kampen’s system size expansion of the Chemical Master Equation. iNA is platform independent and supports the popular SBML format natively. The present implementation is the first to adopt a complementary approach that combines state-of-the-art analysis tools using the computer algebra system Ginac with traditional methods of stochastic simulation. iNA integrates two approximation methods based on the system size expansion, the Linear Noise Approximation and effective mesoscopic rate equations, which to-date have not been available to non-expert users, into an easy-to-use graphical user interface. In particular, the present methods allow for quick approximate analysis of time-dependent mean concentrations, variances, covariances and correlations coefficients, which typically outperforms stochastic simulations. These analytical tools are complemented by automated multi-core stochastic simulations with direct statistical evaluation and visualization. We showcase iNA’s performance by using it to explore the stochastic properties of cooperative and non-cooperative enzyme kinetics and a gene network associated with circadian rhythms. The software iNA is freely available as executable binaries for Linux, MacOSX and Microsoft Windows, as well as the full source code under an open source license.

References

[1]  Gillespie D (2007) Stochastic simulation of chemical kinetics. Annual Review of Physical Chemistry 58: 35–55.
[2]  Grima R, Schnell S (2008) Modelling reaction kinetics inside cells. Essays in Biochemistry 45: 41.
[3]  Eldar A, Elowitz MB (2010) Functional roles for noise in genetic circuits. Nature 467: 167–173.
[4]  van Zon J, Ten Wolde P (2005) Greens-function reaction dynamics: A particle-based approach for simulating biochemical networks in time and space. The Journal of Chemical Physics 123: 234910.
[5]  Andrews S, Addy N, Brent R, Arkin A (2010) Detailed simulations of cell biology with smoldyn 2.1. PLoS computational biology 6: e1000705.
[6]  Kerr R, Bartol T, Kaminsky B, Dittrich M, Chang J, et al. (2008) Fast monte carlo simulation methods for biological reaction-diffusion systems in solution and on surfaces. SIAM Journal on Scientific Computing 30: 3126.
[7]  Ander M, Beltrao P, Di Ventura B, Ferkinghoff-Borg J, Foglierini M, et al. (2004) Smartcell, a framework to simulate cellular processes that combines stochastic approximation with diffusion and localisation: analysis of simple networks. Systems biology 1: 129.
[8]  Hattne J, Fange D, Elf J (2005) Stochastic reaction-diffusion simulation with MesoRD. Bioinformatics 21: 2923–2924.
[9]  van Kampen N (2007) Stochastic processes in physics and chemistry. North-Holland Personal Library, 3rd edition.
[10]  Gardiner C (2007) Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences. Springer, 3rd edition.
[11]  Gillespie D (1992) A rigorous derivation of the chemical master equation. Physica A: Statistical Mechanics and its Applications 188: 404–425.
[12]  Gillespie D (2009) A diffusional bimolecular propensity function. The Journal of Chemical Physics 131: 164109.
[13]  Munsky B, Khammash M (2006) The finite state projection algorithm for the solution of the chemical master equation. The Journal of Chemical Physics 124: 044104.
[14]  (2012) Cmepy v0.3 documentation. Accessed: 10 Feb. http://fcostin.github.com/cmepy/.
[15]  Gillespie C (2009) Moment-closure approximations for mass-action models. Systems Biology, IET 3: 52–58.
[16]  Hespanha J (2008) Moment closure for biochemical networks. In: Communications, Control and Signal Processing, 2008. ISCCSP 2008. 3rd International Symposium on. IEEE, 142–147.
[17]  Hucka M, Finney A, Bornstein B, Keating S, Shapiro B, et al. (2004) Evolving a lingua franca and associated software infrastructure for computational systems biology: the systems biology markup language (sbml) project. Systems Biology 1: 41–53.
[18]  Mendes P, Hoops S, Sahle S, Gauges R, Dada J, et al. (2009) Computational modeling of biochemical networks using copasi. Methods in Molecular Biology 500: 17–59.
[19]  Funahashi A, Matsuoka Y, Jouraku A, Morohashi M, Kikuchi N, et al. (2008) Celldesigner 3.5: a versatile modeling tool for biochemical networks. Proceedings of the IEEE 96: 1254–1265.
[20]  Nicolas R, Donizelli M, Le Novere N (2007) Sbmleditor: effective creation of models in the systems biology markup language (sbml). BMC Bioinformatics 8: 79.
[21]  Gillespie C, Wilkinson D, Proctor C, Shanley D, Boys R, et al. (2006) Tools for the sbml community. Bioinformatics 22: 628–629.
[22]  Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, et al. (2010) BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Systems Biology 4: 92.
[23]  Higham D (2008) Modeling and simulating chemical reactions. SIAM Review, Education Section 50: 347–368.
[24]  Turner T, Schnell S, Burrage K (2004) Stochastic approaches for modelling in vivo reactions. Computational Biology and Chemistry 28: 165–178.
[25]  Gillespie D (1977) Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry 81: 2340–2361.
[26]  (2012) The OpenMP API specification for parallel programming. Accessed: 10 Feb. http://openmp.org/wp/openmp-specificatio?ns/.
[27]  Cao Y, Li H, Petzold L (2004) Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. The Journal of Chemical Physics 121: 4059.
[28]  McQuarrie D (1967) Stochastic approach to chemical kinetics. Journal of Applied Probability 4: 413–478.
[29]  Elf J, Ehrenberg M (2003) Fast evaluation of uctuations in biochemical networks with the linear noise approximation. Genome Research 13: 2475–2484.
[30]  Hayot F, Jayaprakash C (2004) The linear noise approximation for molecular uctuations within cells. Physical Biology 1: 205.
[31]  Grima R (2010) An effective rate equation approach to reaction kinetics in small volumes: Theory and application to biochemical reactions in nonequilibrium steady-state conditions. The Journal of Chemical Physics 133: 035101.
[32]  Grima R (2009) Noise-induced breakdown of the Michaelis-Menten equation in steady-state conditions. Physical Review Letters 102: 218103.
[33]  Grima R (2009) Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments. BMC Systems Biology 3: 101.
[34]  Thomas P, Straube A, Grima R (2010) Stochastic theory of large-scale enzyme-reaction networks: Finite copy number corrections to rate equation models. The Journal of Chemical Physics 133: 195101.
[35]  Ramaswamy R, Gonzalez-Segredo N, Sbalzarini IF, Grima R (2012) Discreteness-induced concentration inversion in mesoscopic chemical systems. Nature Communications 3: 779.
[36]  Press W, Teukolsky S, Vetterling W, Flannery B (2007) Numerical recipes: the art of scientific computing. Cambridge University Press, 3rd edition.
[37]  Fersht A (1999) Structure and mechanism in protein science: a guide to enzyme catalysis and protein folding. WH Freeman.
[38]  Tyson J (2002) Biochemical oscillations. Computational Cell Biology. pp. 230–260.
[39]  Shou M, Mei Q, Ettore Jr M, Dai R, Baillie T, et al. (1999) Sigmoidal kinetic model for two co-operative substrate-binding sites in a cytochrome p450 3a4 active site: an example of the metabolism of diazepam and its derivatives. Biochemical Journal 340: 845.
[40]  Gonze D, Goldbeter A (2006) Circadian rhythms and molecular noise. Chaos 16: 26110.
[41]  De Ronde W, Daniels B, Mugler A, Sinitsyn N, Nemenman I (2009) Mesoscopic statistical properties of multistep enzyme-mediated reactions. Systems Biology, IET 3: 429–437.
[42]  Tam Z, Cai Y, Gunawan R (2010) Elucidating cytochrome c release from mitochondria: Insights from an in silico three-dimensional model. Biophysical Journal 99: 3155–3163.
[43]  Pike L (2009) The challenge of lipid rafts. Journal of Lipid Research 50: S323–S328.
[44]  Lodola A, Shore J, Parker D, Holbrook J (1978) Malate dehydrogenase of the cytosol. A kinetic investigation of the reaction mechanism and a comparison with lactate dehydrogenase. Biochemical Journal 175: 987.
[45]  Luby-Phelps K (1999) Cytoarchitecture and physical properties of cytoplasm: volume, viscosity, diffusion, intracellular surface area. International Review of Cytology 192: 189–221.
[46]  Hutzler J, Tracy T (2002) Atypical kinetic profiles in drug metabolism reactions. Drug metabolism and disposition 30: 355–362.
[47]  Bar-Even A, Noor E, Savir Y, Liebermeister W, Davidi D, et al. (2011) The moderately efficient enzyme: evolutionary and physico-chemical trends shaping enzyme parameters. Biochemistry 50: 4402.
[48]  McKane A, Nagy J, Newman T, Stefanini M (2007) Amplified biochemical oscillations in cellular systems. Journal of Statistical Physics 128: 165–191.
[49]  Davis K, Roussel M (2006) Optimal observability of sustained stochastic competitive inhibition oscillations at organellar volumes. FEBS Journal 273: 84–95.
[50]  Aronson B, Johnson K, Loros J, Dunlap J (1994) Negative feedback defining a circadian clock: autoregulation of the clock gene frequency. Science 263: 1578.
[51]  Goldbeter A (2002) Computational approaches to cellular rhythms. Nature 420: 238–245.
[52]  Sato T, Yamada R, Ukai H, Baggs J, Miraglia L, et al. (2006) Feedback repression is required for mammalian circadian clock function. Nature Genetics 38: 312–319.
[53]  Westermark P, Welsh D, Okamura H, Herzel H (2009) Quantification of circadian rhythms in single cells. PLoS Computational Biology 5: e1000580.
[54]  Thattai M, Van Oudenaarden A (2001) Intrinsic noise in gene regulatory networks. Proceedings of the National Academy of Sciences 98: 8614.
[55]  Shahrezaei V, Swain P (2008) The stochastic nature of biochemical networks. Current opinion in biotechnology 19: 369–374.
[56]  Luo C, Loros J, Dunlap J (1998) Nuclear localization is required for function of the essential clock protein FRQ. The EMBO Journal 17: 1228–1235.
[57]  Mettetal J, Muzzey D, Pedraza J, Ozbudak E, van Oudenaarden A (2006) Predicting stochastic gene expression dynamics in single cells. Proceedings of the National Academy of Sciences of the United States of America 103: 7304.
[58]  Ramaswamy R, Sbalzarini I (2011) Intrinsic noise alters the frequency spectrum of mesoscopic oscillatory chemical reaction systems. Scientific Reports 1.
[59]  Grima R, Thomas P, Straube A (2011) How accurate are the nonlinear chemical fokker-planck and chemical langevin equations? The Journal of Chemical Physics 135: 084103.
[60]  Welsh D, Yoo S, Liu A, Takahashi J, Kay S (2004) Bioluminescence imaging of individual fibroblasts reveals persistent, independently phased circadian rhythms of clock gene expression. Current Biology 14: 2289–2295.
[61]  Locke J, Millar A, Turner M (2005) Modelling genetic networks with noisy and varied experimental data: the circadian clock in arabidopsis thaliana. Journal of Theoretical Biology 234: 383–393.
[62]  Bornstein B, Keating S, Jouraku A, Hucka M (2008) Libsbml: an api library for sbml. Bioinformatics 24: 880.
[63]  Bauer C, Frink A, Kreckel R (2002) Introduction to the ginac framework for symbolic computation within the c++ programming language. Journal of Symbolic Computation 33: 1–12.
[64]  Machné R, Finney A, Müller S, Lu J, Widder S, et al. (2006) The sbml ode solver library: a native api for symbolic and fast numerical analysis of reaction networks. Bioinformatics 22: 1406.
[65]  Guennebaud G, Jacob B (2010) Eigen v3. http://eigen.tuxfamily.org.
[66]  Vallabhajosyula R, Chickarmane V, Sauro H (2006) Conservation analysis of large biochemical networks. Bioinformatics 22: 346–353.
[67]  Fehlberg E (1969) Low-order classical runge-kutta formulas with stepsize control and their application to some heat transfer problems. Technical Report 315, National Aeronautics and Space Administration.
[68]  Shampine LF (1982) Implementation of rosenbrock methods. ACM Transactions on Mathematical Software 8: 93–113.
[69]  Aho A, Lam M, Sethi R, Ullman J (2007) Compilers: principles, techniques, and tools. Pearson/Addison Wesley.
[70]  Van Kampen N (1976) The expansion of the master equation. Advances in Chemical Physics. pp. 245–309.
[71]  Rao C, Arkin A (2003) Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm. The Journal of Chemical Physics 118: 4999.
[72]  Adalsteinsson D, McMillen D, Elston T (2004) Biochemical network stochastic simulator (BioNetS): software for stochastic modeling of biochemical networks. BMC Bioinformatics 5: 24.
[73]  Mauch S (2009) Cain: Stochastic simulations for chemical kinetics. net.
[74]  Mauch S, Stalzer M (2011) Efficient formulations for exact stochastic simulation of chemical systems. Computational Biology and Bioinformatics, IEEE/ACM Transactions on 8: 27–35.
[75]  Caulfield E, Hellander A (2010) Cellmc – a multiplatform model compiler for the cell broadband engine and x86. Bioinformatics 26: 426–428.
[76]  Ramsey S, Orrell D, Bolouri H (2005) Dizzy: stochastic simulation of large-scale genetic regulatory networks. Journal of Bioinformatics and Computational Biology 3: 415–436.
[77]  Liu Q, Qi X, Fan S (2010) Simulating bioreaction processes based on simbiology. Computer Applications and Software 27: 212–214.
[78]  Sanft K, Wu S, Roh M, Fu J, Lim R, et al. (2011) Stochkit2: software for discrete stochastic simulation of biochemical systems with events. Bioinformatics 27: 2457–2458.
[79]  (2012) StochPy user guide v0.1. Accessed: 10 Feb. http://stompy.sourceforge.net/html/userg?uide.html.
[80]  Olivier B, Rohwer J, Hofmeyr J (2005) Modelling cellular systems with PySCeS. Bioinformatics 21: 560–561.

Full-Text

Contact Us

service@oalib.com

QQ:3279437679

WhatsApp +8615387084133