title: The validity of quasi-steady-state approximations in discrete stochastic simulations.
abstract: In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady-state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by nonelementary reaction-rate functions (e.g., Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the nonelementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of nonelementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the nonelementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in nonelementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the prefactor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when nonelementary reaction functions are obtained using the total QSSA. Our work provides an apparently novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.
authors.affiliation: Mathematical Biosciences Institute, The Ohio State University, Columbus, Ohio.
authors.forename: Jae Kyoung
authors.initials: JK
authors.lastname: Kim
authors.affiliation: Department of Mathematics, University of Houston, Houston, Texas; Department of Biology and Biochemistry, University of Houston, Houston, Texas. Electronic address: josic@math.uh.edu.
authors.forename: Krešimir
authors.initials: K
authors.lastname: Josić
authors.affiliation: Department of Biochemistry & Cell Biology, Rice University, Houston, Texas; Institute of Biosciences and Bioengineering, Rice University, Houston, Texas. Electronic address: matthew.bennett@rice.edu.
authors.forename: Matthew R
authors.initials: MR
authors.lastname: Bennett
mesh_terms: D007700:Kinetics; D053858:Metabolic Networks and Pathways; D008962:Models, Theoretical; D013269:Stochastic Processes
references.citation: Biophys J. 2001 Dec;81(6):3116-36
references.pmid: 11720979
references.citation: J Theor Biol. 2010 Jun 21;264(4):1177-89
references.pmid: 20346959
references.citation: J Chem Phys. 2012 Jul 28;137(4):044105
references.pmid: 22852595
references.citation: IET Syst Biol. 2011 Jan;5(1):58
references.pmid: 21261403
references.citation: Bull Math Biol. 2000 May;62(3):483-99
references.pmid: 10812718
references.citation: Biophys J. 2003 Mar;84(3):1606-15
references.pmid: 12609864
references.citation: J Chem Phys. 2005 Nov 15;123(19):194107
references.pmid: 16321076
references.citation: J Chem Phys. 2005 Jan 1;122(1):14116
references.pmid: 15638651
references.citation: Biophys J. 2014 May 6;106(9):2071-81
references.pmid: 24806939
references.citation: J Chem Phys. 2005 May 8;122(18):184102
references.pmid: 15918689
references.citation: J Biol Phys. 2002 Dec;28(4):637-53
references.pmid: 23345804
references.citation: Biophys J. 2008 Oct;95(8):3563-74
references.pmid: 18621809
references.citation: Bull Math Biol. 2003 Nov;65(6):1111-29
references.pmid: 14607291
references.citation: BMC Genomics. 2008 Feb 27;9:102
references.pmid: 18304323
references.citation: PLoS One. 2013 Jun 03;8(6):e65242
references.pmid: 23755198
references.citation: Annu Rev Phys Chem. 2007;58:35-55
references.pmid: 17037977
references.citation: Nature. 2003 Oct 16;425(6959):737-41
references.pmid: 14562106
references.citation: Proc Natl Acad Sci U S A. 2009 Jun 30;106(26):10546-51
references.pmid: 19525397
references.citation: J Theor Biol. 2003 Jan 21;220(2):261-9
references.pmid: 12468297
references.citation: PLoS Comput Biol. 2007 Mar 16;3(3):e45
references.pmid: 17367203
references.citation: J Chem Phys. 2011 Nov 14;135(18):181103
references.pmid: 22088045
references.citation: Genome Res. 2003 Nov;13(11):2475-84
references.pmid: 14597656
references.citation: Biophys J. 2007 May 15;92(10):3501-12
references.pmid: 17350994
references.citation: Methods Enzymol. 2011;487:171-215
references.pmid: 21187226
references.citation: J Theor Biol. 2011 Jun 7;278(1):87-106
references.pmid: 21377474
references.citation: Mol Syst Biol. 2012;8:630
references.pmid: 23212247
references.citation: J Chem Phys. 2008 Sep 7;129(9):095105
references.pmid: 19044893
references.citation: BMC Syst Biol. 2012 May 14;6:39
references.pmid: 22583770
references.citation: Biochem J. 1925;19(2):338-9
references.pmid: 16743508