Realistic boundary conditions for stochastic simulations of reaction-diffusion processes

Abstract

Many cellular and subcellular biological processes can be described in terms of diffusing and chemically reacting species (e.g. enzymes). Such reaction-diffusion processes can be mathematically modelled using either deterministic partial-differential equations or stochastic simulation algorithms. The latter provide a more detailed and precise picture, and several stochastic simulation algorithms have been proposed in recent years. Such models typically give the same description of the reaction-diffusion processes far from the boundary of the simulated domain, but the behaviour close to a reactive boundary (e.g. a membrane with receptors) is unfortunately model-dependent. In this paper, we study four different approaches to stochastic modelling of reaction-diffusion problems and show the correct choice of the boundary condition for each model. The reactive boundary is treated as partially reflective, which means that some molecules hitting the boundary are adsorbed (e.g. bound to the receptor) and some molecules are reflected. The probability that the molecule is adsorbed rather than reflected depends on the reactivity of the boundary (e.g. on the rate constant of the adsorbing chemical reaction and on the number of available receptors), and on the stochastic model used. This dependence is derived for each model.

24 pages, submitted to Physical Biology

References in corpus (4)

• What is the time scale of random sequential adsorption?
• Spatially Distributed Stochastic Systems: equation-free and equation-assisted preconditioned computation
• On chemisorption of polymers to solid surfaces
• Partially Reflected Diffusion

Cited by in corpus (42)

• Grain-Boundary Kinetics: A Unified Approach
• Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions
• A practical guide to stochastic simulations of reaction-diffusion processes
• On the Reaction Diffusion Master Equation in the Microscopic Limit
• Mesoscopic and continuum modelling of angiogenesis
• Simulating with AcCoRD: Actor-Based Communication via Reaction-Diffusion
• The snapping out Brownian motion
• From Molecular Dynamics to Brownian Dynamics
• Irreversibility in dynamical phases and transitions
• Dynamical large deviations of reflected diffusions
• Adaptive two-regime method: application to front propagation
• Efficient Reactive Brownian Dynamics
• A second order scheme for a Robin boundary condition in random walk algorithms
• What is the time scale of random sequential adsorption?
• The narrow capture problem: an encounter-based approach to partially reactive targets
• Modelling of interacting dimer adsorption
• Properties of random sequential adsorption of generalized dimers
• Extreme hitting probabilities for diffusion
• Flow Characteristics in a Crowded Transport Model
• Subdiffusion--reaction process with $A⟶B$ reactions versus subdiffusion--reaction process with $A+B⟶B$ reactions
• Motion of a Brownian molecule in the presence of reactive boundaries
• Exponential and Weibull models for spherical and spherical-shell diffusion-controlled release systems with semi-absorbing boundaries
• Quasi-Monte Carlo methods applied to tau-leaping in stochastic biological systems
• Close encounters of the sticky kind: Brownian motion at absorbing boundaries
• Supervised parallel-in-time algorithm for long-time Lagrangian simulations of stochastic dynamics: Application to hydrodynamics
• Random walks and moving boundaries: Estimating the penetration of diffusants into dense rubbers
• Accuracy of position determination in Ca${}^{2+}$ signaling
• Uniformisation techniques for stochastic simulation of chemical reaction networks
• Positional information readout in $C{a}^{2+}$ signaling
• Modeling and Simulation of Molecular Communication Systems with a Reversible Adsorption Receiver
• Efficient parameter sensitivity computation for spatially-extended reaction networks
• Reactions, Diffusion and Volume Exclusion in a Heterogeneous System of Interacting Particles
• An equivalence framework for an age-structured multi-stage representation of the cell cycle
• Simplified models of diffusion in radially-symmetric geometries
• Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures
• Three-dimensional cytoplasmic calcium propagation with boundaries
• Membrane Fusion-Based Transmitter Design for Molecular Communication Systems
• Opinion Dynamics with Continuous Age Structure
• Adsorption and permeation events in molecular diffusion
• Reverse-time analysis uncovers universality classes in directional biological dynamics
• Encounter-based approach to target search problems: a review
• Molecular Communication with a Reversible Adsorption Receiver