Mathematics.Physics and Chemistry

MD-Leaping Algorithm for Multi-delayed Biochemical Reaction Systems

Expand
  • 1. College of Sciences, Shanghai University, Shanghai 200444, China;
    2. Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, Shanghai 200031, China

Received date: 2012-03-05

  Online published: 2013-02-28

Abstract

To simulate multi-delayed biochemical reaction systems, an MD-Leaping algorithm is proposed to efficiently analyze the process of biochemical reaction in closed systems. In contrast to the D-Leaping algorithm in which the quadruple structure is based on a single delay reaction, a new quadruple structure is built by different time-delay for the multi-delayed reactions. Several reactions can be simulated simultaneously in a time leap, leading to significant increase of the simulation speed. The proposed algorithm extends the D-Leaping method and speeds simulation of multi-delayed biochemical reaction system. Results of two specific biochemical reaction systems show effectiveness of the algorithm for multi-delayed biochemical reaction systems and the mprovement in speed compared to the multi-delay stochastic simulation algorithm (MDSSA) of multi-delayed systems.

Cite this article

ZHUANG Gang, YI Na, DA Liang, WANG Yi-fei . MD-Leaping Algorithm for Multi-delayed Biochemical Reaction Systems[J]. Journal of Shanghai University, 2013 , 19(1) : 75 -79 . DOI: 10.3969/j.issn.1007-2861.2013.01.015

References

 
[1] Gillespie D T. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions [J]. The Journal of Computational Physics,

1976, 22(4): 403-434.

[2] Gillespie D T. Exact stochastic simulation of coupled chemical reactions [J]. The Journal of Physical Chemistry, 1977, 81(25): 2340-2361.

[3] Barrio M, Burrage K, Leier A, et al. Oscillatory regulation of hes1: discrete stochastic delay modeling and simulation [J]. PLoS Computational Biology, 2006,

2(9): 1017-1030.

[4] Cai X D. Exact stochastic simulation of coupled chemical reactions with delays [J]. The Journal of Chemical Physics, 2006, 126(12): 124108.

[5] Anderson D F. A modified next reaction method for simulating chemical systems with time dependent propensities and delays [J]. The Journal of Chemical

Physics, 2007, 127(21): 214107.

[6] Gillespie D T. Approximate accelerated stochastic simulation of chemically reacting systems [J]. The Journal of Chemical Physics, 2001, 115(4): 1716-1733.

[7] Tian T H, Burrage K. Binomial leap methods for simulating stochastic chemical kinetics [J]. The Journal of Chemical Physics, 2004, 121(21): 10356-10364.

[8] Zhou W, Peng X J, Liu X, et al. Final all possible steps approach for accelerating stochastic simulation of coupled chemical reactions [J]. Applied Mathematics and

Mechanics, 2008, 29(3): 379-387.

[9] Leier A, Marquez-Lago T T, BURRAGE K. Generalized binomial tau-leap method for biochemical kinetics incorporating both delay and intrinsic noise [J]. The

Journal of Chemical Physics, 2008, 128(20): 205107.

[10] 彭新俊, 王翼飞. 时滞化学反应系统中的DL算法[J].上海大学学报:自然科学版, 2008, 14(4): 388-393.

[11] Bayati B, Chatelain P, Koumoutsakos P. DLeaping: accelerating stochastic simulation algorithms for reactions with delays [J]. Journal of Computational

Physics, 2009, 228(16): 5908-5916.

[12] Roussel M R, Zhu R. Validation of an algorithm for delay stochastic simulation of transcription and translation in prokaryotic gene expression [J]. Physical Biology,

2006, 3(4): 274-284.

[13] Ribeiro A S. Stochastic and delayed stochastic models of gene expression and regulation [J]. Mathematical Biosciences, 2010, 223(1): 1-11.
Outlines

/