Two classes of quasi-steady-state model reductions for stochastic kinetics
Abstract
The quasi-steady-state approximation (QSSA) is a model reduction technique used to remove highly reactive species from deterministic models of reaction mechanisms. In many reaction networks the highly reactive intermediates (QSSA species) have populations small enough to require a stochastic representation. In this work we apply singular perturbation analysis to remove the QSSA species from the chemical master equation for two classes of problems. The first class occurs in reaction networks where all the species have small populations and the QSSA species sample zero the majority of the time. The perturbation analysis provides a reduced master equation in which the highly reactive species can sample only zero, and are effectively removed from the model. The reduced master equation can be sampled with the Gillespie algorithm. This first stochastic QSSA reduction is applied to several example reaction mechanisms (including Michaelis-Menten kinetics) [Biochem. Z. 49, 333 (1913)]. A general framework for applying the first QSSA reduction technique to new reaction mechanisms is derived. The second class of QSSA model reductions is derived for reaction networks where non-QSSA species have large populations and QSSA species numbers are small and stochastic. We derive this second QSSA reduction from a combination of singular perturbation analysis and the Omega expansion. In some cases the reduced mechanisms and reaction rates from these two stochastic QSSA models and the classical deterministic QSSA reduction are equivalent; however, this is not usually the case.
Additional Information
©2007 American Institute of Physics (Received 7 February 2007; accepted 2 July 2007; published online 6 September 2007) Funding was provided by NIH grant 5R21AI071197-02 and NSF DDDAS grant CNS-0540147. One of the authors (E.L.H.) gratefully acknowledges support from the National Institutes of Health under Ruth L. Kirschstein National Research Service Award 5F32CA120055. The authors would like to thank Dr. J.W. Eaton for help preparing this paper. All simulations were performed using Octave (http://www.octave.org). Octave is freely distributed under the terms of the GNU General Public License.Files
Name | Size | Download all |
---|---|---|
md5:0c221f830c1a9aeffd3607c652d3eae9
|
277.2 kB | Preview Download |
Additional details
- Eprint ID
- 8707
- Resolver ID
- CaltechAUTHORS:MASjcp07
- Created
-
2007-09-07Created from EPrint's datestamp field
- Updated
-
2021-11-08Created from EPrint's last_modified field