title: Simulation of biochemical reactions with time-dependent rates by the rejection-based algorithm.
abstract: We address the problem of simulating biochemical reaction networks with time-dependent rates and propose a new algorithm based on our rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)]. The computation for selecting next reaction firings by our time-dependent RSSA (tRSSA) is computationally efficient. Furthermore, the generated trajectory is exact by exploiting the rejection-based mechanism. We benchmark tRSSA on different biological systems with varying forms of reaction rates to demonstrate its applicability and efficiency. We reveal that for nontrivial cases, the selection of reaction firings in existing algorithms introduces approximations because the integration of reaction rates is very computationally demanding and simplifying assumptions are introduced. The selection of the next reaction firing by our approach is easier while preserving the exactness.
authors.affiliation: The Microsoft Research - University of Trento Centre for Computational and Systems Biology, Piazza Manifattura 1, Rovereto 38068, Italy.
authors.forename: Vo Hong
authors.initials: VH
authors.lastname: Thanh
authors.affiliation: The Microsoft Research - University of Trento Centre for Computational and Systems Biology, Piazza Manifattura 1, Rovereto 38068, Italy.
authors.forename: Corrado
authors.initials: C
authors.lastname: Priami
mesh_terms: D000465:Algorithms; D018562:Disease Transmission, Infectious; D058872:Epidemics; D005786:Gene Expression Regulation; D007700:Kinetics; D008954:Models, Biological; D013269:Stochastic Processes; D014158:Transcription, Genetic