Abstract
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.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
