A comparison of bimolecular reaction models for stochastic reaction-diffusion systems

Area of Science:

• Biophysics
• Chemical Kinetics
• Computational Biology

Background:

• Stochastic reaction-diffusion models are crucial for understanding molecular interactions in biological systems.
• Two prominent models, Smoluchowski and Doi, describe molecular motion and reactions.
• These models differ in how they represent chemical reaction events.

Purpose of the Study:

• To rigorously investigate the relationship between the Smoluchowski and Doi reaction-diffusion models.
• To analyze the convergence of the Doi model to the Smoluchowski model under specific conditions.
• To quantify differences in model predictions using biologically relevant parameters.

Main Methods:

• Analytical proof for model convergence.
• Numerical simulations of stochastic reaction-diffusion processes.
• Analysis of reaction time statistics and model solution differences.

Main Results:

• Proved the Doi model solution converges to the Smoluchowski model solution as the reaction rate (λ) approaches infinity, with a rigorous error bound.
• Numerical simulations confirmed convergence for a protein-DNA binding site scenario.
• Observed that differences between model solutions increase as the reaction radius decreases for large λ.

Conclusions:

• The Doi model provides a rigorous approximation of the Smoluchowski model in the high reaction rate limit.
• Understanding the relationship between these models is vital for accurate modeling of molecular encounters in biology.
• Model discrepancies highlight the importance of reaction radius and rate in diffusion-limited reactions.