title: Modeling stochastic noise in gene regulatory systems.
abstract: The Master equation is considered the gold standard for modeling the stochastic mechanisms of gene regulation in molecular detail, but it is too complex to solve exactly in most cases, so approximation and simulation methods are essential. However, there is still a lack of consensus about the best way to carry these out. To help clarify the situation, we review Master equation models of gene regulation, theoretical approximations based on an expansion method due to N.G. van Kampen and R. Kubo, and simulation algorithms due to D.T. Gillespie and P. Langevin. Expansion of the Master equation shows that for systems with a single stable steady-state, the stochastic model reduces to a deterministic model in a first-order approximation. Additional theory, also due to van Kampen, describes the asymptotic behavior of multistable systems. To support and illustrate the theory and provide further insight into the complex behavior of multistable systems, we perform a detailed simulation study comparing the various approximation and simulation methods applied to synthetic gene regulatory systems with various qualitative characteristics. The simulation studies show that for large stochastic systems with a single steady-state, deterministic models are quite accurate, since the probability distribution of the solution has a single peak tracking the deterministic trajectory whose variance is inversely proportional to the system size. In multistable stochastic systems, large fluctuations can cause individual trajectories to escape from the domain of attraction of one steady-state and be attracted to another, so the system eventually reaches a multimodal probability distribution in which all stable steady-states are represented proportional to their relative stability. However, since the escape time scales exponentially with system size, this process can take a very long time in large systems.
authors: Computational Biology Lab, Bio-X Program, Stanford University, Stanford, CA 94305, USA.
authors: Arwen
authors: A
authors: Meister
authors: Computational Biology Lab, Bio-X Program, Stanford University, Stanford, CA 94305, USA.
authors: Chao
authors: C
authors: Du
authors: Computational Biology Lab, Bio-X Program, Stanford University, Stanford, CA 94305, USA.
authors: Ye Henry
authors: YH
authors: Li
authors: Computational Biology Lab, Bio-X Program, Stanford University, Stanford, CA 94305, USA.
authors: Wing Hung
authors: WH
authors: Wong
keywords: Gillespie algorithm; Langevin equation; Master equation; gene regulation; simulation; stochastic modeling
references: J Chem Phys. 2012 Dec 21;137(23):234109
references: 23267473
references: Biophys J. 2001 Dec;81(6):3116-36
references: 11720979
references: Nature. 2002 Nov 14;420(6912):224-30
references: 12432407
references: Nat Genet. 2002 May;31(1):69-73
references: 11967532
references: Cell. 2008 Oct 17;135(2):216-26
references: 18957198
references: Mol Cell. 2012 Feb 24;45(4):483-93
references: 22365828
references: PLoS One. 2012;7(2):e31407
references: 22363638
references: Nat Biotechnol. 2005 Mar;23(3):377-83
references: 15765094
references: J Mol Biol. 1985 Jan 20;181(2):211-30
references: 3157005
references: Proc Natl Acad Sci U S A. 1982 Feb;79(4):1129-33
references: 6461856
references: Mol Syst Biol. 2009;5:318
references: 19888213
references: Proc Natl Acad Sci U S A. 2002 Oct 1;99(20):12795-800
references: 12237400
references: Mol Syst Biol. 2007;3:78
references: 17299415
references: Nature. 2000 Jan 20;403(6767):339-42
references: 10659857
references: J Comput Sci Technol. 2010 Jan;25(1):154-168
references: 24999297
references: J Theor Biol. 2004 Dec 21;231(4):563-8
references: 15488533
references: Nature. 2000 Jan 20;403(6767):335-8
references: 10659856
references: J Chem Phys. 2006 Nov 28;125(20):204104
references: 17144687
references: BMC Syst Biol. 2010 Aug 09;4:108
references: 20696063
references: PLoS Biol. 2007 Jan;5(1):e8
references: 17214507
references: Science. 2002 Aug 16;297(5584):1183-6
references: 12183631
references: Genetics. 1998 Aug;149(4):1633-48
references: 9691025
references: Nature. 2003 Apr 10;422(6932):633-7
references: 12687005
references: Proc Natl Acad Sci U S A. 2001 Jul 17;98(15):8614-9
references: 11438714
references: Proc Natl Acad Sci U S A. 1997 Aug 19;94(17):9147-52
references: 9256450
references: Science. 2005 Mar 25;307(5717):1965-9
references: 15790857
references: Nature. 2002 Nov 14;420(6912):231-7
references: 12432408
references: J Chem Phys. 2008 Sep 7;129(9):095105
references: 19044893
references: Nature. 2004 Feb 19;427(6976):737-40
references: 14973486
references: Proc Natl Acad Sci U S A. 2010 Jul 27;107(30):13300-5
references: 20616054
references: Curr Opin Genet Dev. 2005 Apr;15(2):125-35
references: 15797195
references: Mol Cell. 2009 Sep 24;35(6):741-53
references: 19782025
references: Science. 2012 Apr 13;336(6078):183-7
references: 22499939
references: BMC Bioinformatics. 2009 Oct 19;10:343
references: 19840370
references: PLoS One. 2012;7(9):e44812
references: 22984566
references: Biol Direct. 2007 May 31;2:13
references: 17540019
references: Science. 2003 Jul 4;301(5629):102-5
references: 12843395
references: Nature. 2004 Jan 29;427(6973):415-8
references: 14749823
references: Phys Rev Lett. 2012 Dec 14;109(24):248107
references: 23368390
references: Curr Opin Genet Dev. 2005 Apr;15(2):116-24
references: 15797194
references: PLoS One. 2008;3(10 ):e3478
references: 18941526
references: Nat Rev Genet. 2005 Jun;6(6):451-64
references: 15883588
references: Science. 2005 Mar 25;307(5717):1962-5
references: 15790856
references: Biochem J. 1913 Oct;7(5):471-80
references: 16742267
references: Nat Genet. 2008 Dec;40(12):1493-8
references: 19029898