The Time Dependent Propensity Function for Acceleration of Spatial Stochastic Simulation of Reaction-Diffusion...

Area of Science:

• Computational biology
• Biophysics
• Biochemical reaction modeling

Background:

• Spatial stochastic simulation is crucial for understanding biological processes at the molecular level.
• The inhomogeneous stochastic simulation algorithm (ISSA) is a standard method but struggles with high diffusion rates.
• Simulating frequent diffusion events in ISSA leads to high computational costs.

Purpose of the Study:

• To develop a more efficient algorithm for spatial stochastic simulation.
• To reduce the computational cost associated with frequent diffusion events in ISSA.
• To maintain accuracy while improving simulation speed.

Main Methods:

• Proposed using a time-dependent propensity function within the simulation algorithm.
• Replaced individual diffusion event simulations with larger time steps between reaction events.
• Validated the new algorithm against widely-used exact simulation methods.

Main Results:

• Achieved orders of magnitude efficiency gains compared to existing exact algorithms.
• Demonstrated excellent scalability with increasing grid resolution.
• Maintained a high level of accuracy in simulations.

Conclusions:

• The proposed method significantly enhances the efficiency of spatial stochastic simulations.
• This approach is particularly beneficial for systems with frequent diffusion relative to reaction events.
• The algorithm offers a computationally efficient and accurate alternative for complex biological modeling.