Reaction rates for mesoscopic reaction-diffusion kinetics

Stefan Hellander,Andreas Hellander,Linda Petzold

DOI: https://doi.org/10.1103/PhysRevE.91.023312

2015-01-29

Abstract:The mesoscopic reaction-diffusion master equation (RDME) is a popular modeling framework, frequently applied to stochastic reaction-diffusion kinetics in systems biology. The RDME is derived from assumptions about the underlying physical properties of the system, and it may produce unphysical results for models where those assumptions fail. In that case, other more comprehensive models are better suited, such as hard-sphere Brownian dynamics (BD). Although the RDME is a model in its own right, and not inferred from any specific microscale model, it proves useful to attempt to approximate a microscale model by a specific choice of mesoscopic reaction rates. In this paper we derive mesoscopic reaction rates by matching certain statistics of the RDME solution to statistics of the solution of a widely used microscopic BD model: the Smoluchowski model with a mixed boundary condition at the reaction radius of two molecules. We also establish fundamental limits for the range of mesh resolutions for which this approach yields accurate results, and show both theoretically and in numerical examples that as we approach the lower fundamental limit, the mesoscopic dynamics approach the microscopic dynamics.

Numerical Analysis,Subcellular Processes

What problem does this paper attempt to address?

The key problem that this paper attempts to solve is how to accurately simulate reaction - diffusion dynamics at the mesoscopic scale, especially how to determine appropriate mesoscopic reaction rates to match the statistical properties of the microscopic model. Specifically, the authors adjust the reaction rates in the mesoscopic reaction - diffusion master equation (RDME) to make it more accurately reflect the reaction behavior between molecules in the Smoluchowski model. The following are the core problems and solutions in this paper:

### 1. **Research Background and Motivation**
The mesoscopic reaction - diffusion master equation (RDME) is a commonly used modeling framework in computational systems biology to describe stochastic reaction - diffusion processes under spatial heterogeneity. However, when the physical properties of the system do not conform to the assumptions of RDME, RDME may produce unphysical results. In this case, a more detailed microscopic model (such as hard - sphere Brownian dynamics BD) is more appropriate. Nevertheless, in order to make RDME better approximate the microscopic model, it is crucial to select appropriate mesoscopic reaction rates.

### 2. **Core Problem**
The main objective of this paper is to derive scale - dependent mesoscopic reaction rates by matching certain statistical properties of the RDME solution with those of the widely used microscopic BD model (Smoluchowski model). In addition, the author also explores the effective range of this matching method at different grid resolutions and shows that as the grid size approaches the lower limit, the mesoscopic dynamics gradually approximates the microscopic dynamics.

### 3. **Specific Solutions**
- **Derivation of Mesoscopic Reaction Rates**: By matching the average binding time, the steady - state level of the Smoluchowski model and the RDME model, and ensuring that the mesoscopic dissociation rate does not exceed the intrinsic dissociation rate, the authors derived mesoscopic reaction rates applicable to reversible reactions.

- **Determination of the Lower Limit of Grid Size**: The author established the basic limit of grid size, that is, there is a non - trivial lower limit. Above this lower limit, an appropriate mesoscopic reaction rate can be selected so that the average binding time of the RDME and the Smoluchowski model match. For irreversible reactions, this lower limit is \( h^*_{\infty}\approx3.2\sigma\) (3D) and \( h^*_{\infty}\approx5.1\sigma\) (2D).

- **Error Estimation**: The author also provided error estimates to evaluate the required grid resolution to achieve a specific precision of the recombination time distribution. The results show that as the mesoscopic grid size approaches the lower limit, the mesoscopic average recombination time gradually approaches the microscopic average recombination time.

### 4. **Theoretical and Numerical Verification**
Through theoretical analysis and numerical examples, the author proved that as the grid size approaches the lower limit, the mesoscopic dynamics can indeed approximate the microscopic dynamics. At the same time, they also showed that when the grid size is below the lower limit, the accuracy of the results will decrease.

### 5. **Conclusion**
By deriving mesoscopic reaction rates and determining the lower limit of grid size, this paper provides an improved method that enables RDME to accurately simulate microscopic reaction - diffusion dynamics in a wider range of scales. This not only improves the applicability and accuracy of RDME but also provides an important theoretical basis for future research.

### Related Formulas
- Boundary Conditions of the Smoluchowski Model:
\[
K\frac{\partial p}{\partial n}\bigg|_{|r| = \sigma}=k_rp(r,t|r_n,t_n)
\]
where \(K=\begin{cases}
4\pi\sigma^2D&(3D)\\
2\pi\sigma D&(2D)
\end{cases}\)

- Mesoscopic Reaction Rate:
\[
\rho^{(d)}(k_r,h)=\frac{k_r}{h^d\left(1+\frac{k_r}{D}G^{(d)}(h,\sigma)\right)}
\]
where \(G^{(d)}(h,\sigma)=\begin{cases}
...