Accuracy Analysis of Hybrid Stochastic Simulation Algorithm on Linear Chain Reaction Systems

Noise in cellular systems is often modeled and simulated with Gillespie’s stochastic simulation algorithm (SSA), but the low efficiency of the SSA limits its application to large biochemical networks. To improve the efficiency of stochastic simulations, Haseltine and Rawlings (HR) proposed a hybrid...

Full description

Saved in:
Bibliographic Details
Published in:Bulletin of mathematical biology Vol. 81; no. 8; pp. 3024 - 3052
Main Authors: Chen, Minghan, Wang, Shuo, Cao, Yang
Format: Journal Article
Language:English
Published: New York Springer US 01.08.2019
Springer Nature B.V
Subjects:
Algorithms
Cell Biology
Chains
Computer simulation
Differential equations
Life Sciences
Mathematical and Computational Biology
Mathematical models
Mathematics
Mathematics and Statistics
Model accuracy
Molecular chains
Ordinary differential equations
Organic chemistry
Special Issue: Gillespie and His Algorithms
Stochastic models
SSA
Accuracy analysis
Next slow reaction firing time
Linear chain reaction systems
Hybrid stochastic algorithm
ISSN:0092-8240, 1522-9602, 1522-9602
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Noise in cellular systems is often modeled and simulated with Gillespie’s stochastic simulation algorithm (SSA), but the low efficiency of the SSA limits its application to large biochemical networks. To improve the efficiency of stochastic simulations, Haseltine and Rawlings (HR) proposed a hybrid algorithm, which combines ordinary differential equations for traditional deterministic models and the SSA for stochastic models. In this paper, accuracy of the HR hybrid method is studied based on a linear chain reaction system. Mathematical analysis and numerical results both show that the HR hybrid method is accurate if either the quantity of reactant molecules in fast reactions is above a certain threshold, or the reaction rates of fast reactions are much larger than those of slow reactions. This analysis also shows that the HR hybrid method approximates the chemical master equation well for a much greater region in system parameter space than the slow-scale SSA and the stochastic quasi-steady-state assumption methods.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
content type line 23
ISSN:0092-8240
1522-9602
1522-9602
DOI:10.1007/s11538-018-0461-z