The present study proposes a stochastic simulation scheme to model reactive boundaries through a position jump process which can be readily implemented into the Inhomogeneous Stochastic Simulation Algorithm by modifying the propensity of the diffusive jump over the reactive boundary. As compared to the literature, the present approach does not require any correction factors for the propensity. Also, the current expression relaxes the constraint on the compartment size allowing the problem to be solved with a coarser grid and therefore saves considerable computational cost. The modified algorithm is then applied to simulate three reaction-diffusion systems with reactive boundaries.
References
[1]
Fedoroff, N. and Fontana, W. (2002) Small Numbers of Big Molecules. Science, 297, 1129-1131. https://doi.org/10.1126/science.1075988
[2]
Gillespie, D.T. (1977) Exact Stochastic Simulation of Coupled Chemical Reactions. Journal of Physical Chemistry, 81, 2340-2361. https://doi.org/10.1021/j100540a008
[3]
Andrews, S.S. and Bray, D. (2004) Stochastic Simulation of Chemical Reactions with Spatial Resolution and Single Molecule Detail. Physical Biology, 1, 137-151. https://doi.org/10.1088/1478-3967/1/3/001
[4]
Cao, Y., Gillespie, D.T. and Petzold, L.R. (2007) Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection. Journal of Chemical Physics, 126, Article ID: 224101. https://doi.org/10.1063/1.2745299
[5]
Smith, C.A. and Yates, C.A. (2018) Spatially Extended Hybrid Methods: A Review. Journal of the Royal Society Interface, 15, Article ID: 20170931. https://doi.org/10.1098/rsif.2017.0931
[6]
Gibson, M. and Bruck, J. (2000) Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Channels. Journal of Physical Chemistry, 104, 1876-1889. https://doi.org/10.1021/jp993732q
[7]
Elf, J. and Ehrenberg, M. (2004) Spontaneous Separation of Bi-Stable Biochemical Systems into Spatial Domains of Opposite Phases. Systems Biology, 1, 2230-2236. https://doi.org/10.1049/sb:20045021
[8]
Lampoudi, S., Gillespie, D.T. and Petzold, L.R. (2009) The Multinomial Simulation Algorithm for Discrete Stochastic Simulation of Reaction-Diffusion Systems. Journal of Chemical Physics, 130, Article ID: 094104. https://doi.org/10.1063/1.3074302
[9]
Drawert, B., Lawson, M.J., Petzold, L.R. and Kammash, M. (2010) The Diffusive Finite State Projection Algorithm for Efficient Simulation of the Stochastic Reaction-Diffusion Master Equation. Journal of Chemical Physics, 132, Article ID: 074101. https://doi.org/10.1063/1.3310809
[10]
Koh, W. and Blackwell, K.T. (2012) An Accelerated Algorithm for Discrete Stochastic Simulation of Reaction-Diffusion Systems Using Gradient-Based Diffusion and Tau-Leaping. Journal of Chemical Physics, 134, Article ID: 154103. https://doi.org/10.1063/1.3572335
[11]
Boccardo, G., Sokolov, I.M. and Paster, A. (2018) A Second Order Scheme for a Robin Boundary Condition in Random Walk Algorithms. Journal of Computational Physics, 374, 1152-1165. https://doi.org/10.1016/j.jcp.2018.08.009
[12]
Hoffmann, M. and Schawrz, U.S. (2014) Oscillations of Min-Proteins in Micropatterned Environments: A Three-Dimensional Particle-Based Stochastic Simulation Approach. Soft Matter, 10, 2388-2396. https://doi.org/10.1039/C3SM52251B
[13]
Erban, R., Chapman, S.J., Fisher, K., Kevrekidis, I. and Seymour, L.W. (2007) Dynamics of Polydisperse Irreversible Adsorption: A Pharmacological Example. Mathematical Models and Methods in Applied Sciences, 17, 759-781. https://doi.org/10.1142/S0218202507002091
[14]
Fisher, K., Stallwood, Y., Green, N., Ulbrich, K., Mautner, V. and Seymour, L.W. (2011) Polymer-Coated Adenovirus Permits Efficient Retargeting and Evades Neutralizing Antibodies. Gene Therapy, 8, 341-348. https://doi.org/10.1038/sj.gt.3301389
[15]
Erban, R. and Chapman, S.J. (2007) Reactive Boundary Conditions for Stochastic Simulations of Reaction-Diffusion Processes. Physical Biology, 4, 16-28. https://doi.org/10.1088/1478-3975/4/1/003
[16]
Singer, A., Schuss, Z., Lumsden, C., Osipov, A. and Holcman, D. (2008) Partially Reflected Diffusion. SIAM Journal on Applied Mathematics, 68, 844-868. https://doi.org/10.1137/060663258
[17]
Sayyidmousavi, A. and Rohlf, K. (2018) Reactive Multi-Particle Collision Dynamics with Reactive Boundary Conditions. Physical Biology, 15, Article ID: 046007. https://doi.org/10.1088/1478-3975/aabc35
[18]
Lier, A. and Marquez-Lago, T.T. (2011) Correction Factors for Boundary Diffusion in Reaction-Diffusion Master Equations. Journal of Chemical Physics, 135, Article ID: 134109. https://doi.org/10.1063/1.3634003
[19]
Isaacson, S.A. and Peskin, C.S. (2006) Incorporating Diffusion in Complex Geometries into Stochastic Chemical Kinetics Simulations. SIAM Journal on Scientific Computing, 28, 47-74. https://doi.org/10.1137/040605060
[20]
Fange, D. and Elf, J. (2006) Noise-Induced Min Phenotypes in E. coli. PLOS Computational Biology, 2, 637-648. https://doi.org/10.1371/journal.pcbi.0020080
[21]
Gillespie, D.T. (2001) Approximate Accelerated Stochastic Simulation of Chemically Reacting Systems. Journal of Chemical Physics, 155, 1716-1733. https://doi.org/10.1063/1.1378322
[22]
Chatterjee, A. and Vlachos, D.G. (2005) Binomial Distribution Based τ-Leap Accelerated Stochastic Simulation. Journal of Chemical Physics, 122, Article ID: 024112. https://doi.org/10.1063/1.1833357
[23]
Chiam, K.H., Tan, C.M., Bhargava, V. and Rajagopal, G. (2006) Hybrid Simulations of Stochastic Reaction-Diffusion Processes for Modeling Intracellular Signaling Pathways. Physical Review E, 74, Article ID: 051910. https://doi.org/10.1103/PhysRevE.74.051910
[24]
Rossinelli, D., Bayati, B. and Koumoutsakos, P. (2008) Accelerated Stochastic and Hybrid Methods for Spatial Simulations of Reaction-Diffusion Systems. Chemical Physics Letter, 451, 136-140. https://doi.org/10.1016/j.cplett.2007.11.055
[25]
Ferm, L., Hellander, A. and Lotstedt, P. (2010) An Adaptive Algorithm for Simulation of Stochastic Reaction-Diffusion Processes. Journal of Computational Physics, 229, 343-360. https://doi.org/10.1016/j.jcp.2009.09.030
[26]
Strehl, R. and Ilie, S. (2015) Hybrid Stochastic Simulation of Reaction-Diffusion Systems with Slow and Fast Dynamics. Journal of Chemical Physics, 143, Article ID: 234108. https://doi.org/10.1063/1.4937491
[27]
Padgett, J.M.A. and Ilie, S. (2016) An Adaptive Tau-Leaping Method for Stochastic Simulations of Reaction-Diffusion Systems. AIP Advances, 6, Article ID: 035217. https://doi.org/10.1063/1.4944952
[28]
Barrio, M., Burrage, K., Leier, A. and Tian, T. (2006) Oscillatory Regulation of Hes1: Discrete Stochastic Delay Modelling and Simulation. PLOS Computational Biology, 2, 1017-1030. https://doi.org/10.1371/journal.pcbi.0020117
[29]
Sayyidmousavi, A. and Ilie, S. (2017) An Efficient Hybrid Method for Stochastic Reaction-Diffusion Biochemical Systems with Delay. AIP Advances, 7, Article ID: 125305. https://doi.org/10.1063/1.5001760
[30]
Taylor, P.R., Baker, R.E. and Yates, C.A. (2014) Deriving Appropriate Boundary Conditions and Accelerating Position-Jump Simulations of Diffusion Using Non-Local Jumping. Physical Biology, 12, Article ID: 016006. https://doi.org/10.1088/1478-3975/12/1/016006
[31]
Smith, S., Cianci, C. and Grima, R. (2015) Analytical Approximations for Spatial Stochastic Gene Expression in Single Cells and Tissues. Journal of the Royal Society Interface, 13, Article ID: 20151051. https://doi.org/10.1098/rsif.2015.1051
