Exact rule-based stochastic simulations for the system with unlimited number of molecular species

We introduce expandable partial propensity direct method (EPDM) - a new exact stochastic simulation algorithm suitable for systems involving many interacting molecular species. The algorithm is especially efficient for sparsely populated systems, where the number of species that may potentially be generated is much greater than the number of species actually present in the system at any given time. The number of operations per reaction scales linearly with the number of species, but only those which have one or more molecules. To achieve this kind of performance we are employing a data structure which allows to add and remove species and their interactions on the fly. When a new specie is added, its interactions with every other specie are generated dynamically by a set of user-defined rules. By removing the records involving the species with zero molecules, we keep the number of species as low as possible. This enables simulations of systems for which listing all species is not practical. The algorithm is based on partial propensities direct method (PDM) by Ramaswamy et al. for sampling trajectories of the chemical master equation.

[1]  D. Gillespie,et al.  Avoiding negative populations in explicit Poisson tau-leaping. , 2005, The Journal of chemical physics.

[2]  Mauro Mobilia,et al.  Stochastic dynamics of the prisoner's dilemma with cooperation facilitators. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  Sheng Wu,et al.  StochKit2: software for discrete stochastic simulation of biochemical systems with events , 2011, Bioinform..

[4]  A. Lees.  Modelling and Analysis , 2015 .

[5]  Linda R Petzold,et al.  Efficient step size selection for the tau-leaping simulation method. , 2006, The Journal of chemical physics.

[6]  Anne Auger,et al.  R-leaping: accelerating the stochastic simulation algorithm by reaction leaps. , 2006, The Journal of chemical physics.

[7]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[8]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[9]  G. Henkelman,et al.  Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table , 2001 .

[10]  Dranreb Earl Juanico,et al.  Critical network effect induces business oscillations in multi-level marketing systems , 2012, 1209.5576.

[11]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[12]  Jane Hillston,et al.  Fluid flow approximation of PEPA models , 2005, Second International Conference on the Quantitative Evaluation of Systems (QEST'05).

[13]  Sheldon M. Ross,et al.  Stochastic Processes , 2018, Gauge Integral Structures for Stochastic Calculus and Quantum Electrodynamics.

[14]  Dan ie l T. Gil lespie.  A rigorous derivation of the chemical master equation , 1992 .

[15]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[16]  K. Dill,et al.  How did prebiotic polymers become informational foldamers , 2016, 1604.08890.

[17]  Pejman Rohani,et al.  The Role of Environmental Transmission in Recurrent Avian Influenza Epidemics , 2009, PLoS Comput. Biol..

[18]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[19]  Ulf Dieckmann,et al.  Adaptive Evolution of Social Traits: Origin, Trajectories, and Correlations of Altruism and Mobility , 2004, The American Naturalist.

[20]  M. Keeling,et al.  Modeling Infectious Diseases in Humans and Animals , 2007 .

[21]  A Flahault,et al.  Understanding the dynamics of Ebola epidemics , 2006, Epidemiology and Infection.

[22]  T. Reichenbach,et al.  Mobility promotes and jeopardizes biodiversity in rock–paper–scissors games , 2007, Nature.

[23]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[24]  Ivo F Sbalzarini,et al.  A partial-propensity variant of the composition-rejection stochastic simulation algorithm for chemical reaction networks. , 2010, The Journal of chemical physics.

[25]  Ira B Schwartz,et al.  Recruitment dynamics in adaptive social networks , 2011, Journal of physics. A, Mathematical and theoretical.

[26]  Ivo F Sbalzarini,et al.  A new class of highly efficient exact stochastic simulation algorithms for chemical reaction networks. , 2009, The Journal of chemical physics.

[27]  M. A. Trump,et al.  The Dynamical Theory , 1999 .

[28]  Guy Karlebach,et al.  Modelling and analysis of gene regulatory networks , 2008, Nature Reviews Molecular Cell Biology.

[29]  Hong Li,et al.  Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. , 2004, The Journal of chemical physics.

[30]  Giuseppe Carbone,et al.  Model of human collective decision-making in complex environments , 2015, The European Physical Journal B.

[31]  Muruhan Rathinam,et al.  Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method , 2003 .

[32]  Spring Berman,et al.  Bio-Inspired Group Behaviors for the Deployment of a Swarm of Robots to Multiple Destinations , 2007, Proceedings 2007 IEEE International Conference on Robotics and Automation.

[33]  Spring Berman,et al.  Optimized Stochastic Policies for Task Allocation in Swarms of Robots , 2009, IEEE Transactions on Robotics.

[34]  U. Dieckmann,et al.  The Dynamical Theory of Coevolution : A Derivation from Stochastic Ecological Processes , 1996 .

[35]  Yi-fei Wang,et al.  Efficient binomial leap method for simulating chemical kinetics. , 2007, The Journal of chemical physics.

[36]  Xiaodong Cai,et al.  K-leap method for accelerating stochastic simulation of coupled chemical reactions. , 2007, The Journal of chemical physics.

[37]  Nagiza F. Samatova,et al.  The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior , 2006, Comput. Biol. Chem..

[38]  Roberto Barbuti,et al.  Tumour suppression by immune system through stochastic oscillations. , 2010, Journal of theoretical biology.

[39]  A J McKane,et al.  Predator-prey cycles from resonant amplification of demographic stochasticity. , 2005, Physical review letters.

[40]  彭新俊,et al.  L-leap: accelerating the stochastic simulation of chemically reacting systems , 2007 .

[41]  Darren J. Wilkinson.  Stochastic Modelling for Systems Biology , 2006 .