针对多重时滞生化反应系统, 提出MD-Leaping模拟算法, 该算法能快速有效地刻画封闭系统中生化反应的进程和结果. 相较于D-Leaping算法, 针对单时滞反应建立以反应序列为主体的四元结构体, 该算法通过对多重时滞反应系统建立以不同时滞时间为主体的四元结构体, 使得在每一个跳跃时间步长内可以同时模拟多个生化反应, 提高了模拟速度. 该算法拓展D-Leaping算法, 并首次解决多重时滞生化反应系统的加速模拟问题. 通过对2个具体生化反应系统的模拟表明, 该算法可以广泛地应用于多重时滞生化反应系统, 并且与已有的多重时滞系统非跳跃方法(multi-delay stochastic simulation algorithm, MDSSA)相比, 其运行效率有显著提高.
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.
[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] 彭新俊, 王翼飞. 时滞化学反应系统中的D L算法[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.
