[论文解读] Coupling particle-based reaction-diffusion simulations with reservoirs mediated by reaction-diffusion PDEs

Margarita Kostré, Christof Schütte+2 more|arXiv (Cornell University)|May 29, 2020

Block Copolymer Self-Assembly参考文献 53被引用 14

一句话总结

本文提出了一种在粒子基反应-扩散（PBRD）模拟与代表时变、空间非均匀池的宏观反应-扩散PDE之间实现数学一致耦合的方法。通过推导开放系统及二级反应的平均场理论，作者实现了对开放生化系统中跨尺度正确通量交换的精确、可扩展模拟。

ABSTRACT

Open biochemical systems of interacting molecules are ubiquitous in life-related processes. However, established computational methodologies, like molecular dynamics, are still mostly constrained to closed systems and timescales too small to be relevant for life processes. Alternatively, particle-based reaction-diffusion models are currently the most accurate and computationally feasible approach at these scales. Their efficiency lies in modeling entire molecules as particles that can diffuse and interact with each other. In this work, we develop modeling and numerical schemes for particle-based reaction-diffusion in an open setting, where the reservoirs are mediated by reaction-diffusion PDEs. We derive two important theoretical results. The first one is the mean-field for open systems of diffusing particles; the second one is the mean-field for a particle-based reaction-diffusion system with second-order reactions. We employ these two results to develop a numerical scheme that consistently couples particle-based reaction-diffusion processes with reaction-diffusion PDEs. This allows modeling open biochemical systems in contact with reservoirs that are time-dependent and spatially inhomogeneous, as in many relevant real-world applications.

研究动机与目标

• 开发一种将粒子基反应-扩散（PBRD）模拟与代表池的宏观反应-扩散PDE之间进行一致数学耦合的框架。
• 解决在粒子数不守恒的开放系统中建模的挑战，需采用广义系综方法。
• 建立微观反应速率与宏观PDE参数之间的正确关系，特别是针对二级反应。
• 确保PBRD模拟与PDE池之间的数值一致性，实现复杂生化系统准确的多尺度建模。

提出的方法

• 推导扩散粒子开放系统的平均场理论，以建模与池之间的粒子交换。
• 为具有二级反应的PBRD系统开发平均场公式，以关联微观与宏观速率常数。
• 利用推导出的平均场关系，构建PBRD模拟与PDE池之间的一致耦合方案。
• 使用Crank-Nicolson有限差分格式求解PDE，确保空间和时间上的二阶精度。
• 采用Strang分裂法处理PDE中的反应与扩散项，保持数值稳定性和准确性。
• 通过验证PDE池在PBRD模拟中正确调节粒子流入与流出，确保热力学一致性。

实验结果

研究问题

• RQ1 如何将粒子基反应-扩散模拟与代表时变、空间依赖池的宏观反应-扩散PDE实现一致耦合？
• RQ2 PBRD中微观反应速率与PDE中宏观速率常数之间的正确平均场关系是什么，特别是针对二级反应？
• RQ3 如何在广义系综框架下正确建模开放PBRD系统中的粒子数涨落？
• RQ4 何种数值格式可确保在多尺度下PBRD模拟与PDE池耦合时的一致性与稳定性？
• RQ5 包含空间非均匀且时变池如何影响PBRD模拟的精度与动力学特性？

主要发现

• 作者推导了扩散粒子开放系统的平均场理论，实现了在广义系综中对与池之间粒子交换的一致建模。
• 他们建立了PBRD中微观二级反应速率与宏观PDE速率常数之间的正确平均场关系，表明扩散系数不参与速率关系——这是对先前假设的关键修正。
• 所提出的耦合方案确保了PBRD模拟与PDE池之间的数学一致性，支持对具有空间和时间变化边界条件的开放生化系统进行精确模拟。
• 数值验证表明，Crank-Nicolson方法结合Strang分裂在空间和时间上均达到二阶精度，确保了PDE解的稳定与可靠。
• 该方法可实现对复杂、真实生物情景的模拟，如具有动态细胞外环境的细胞信号传导，其中池浓度随空间和时间变化。
• 该实现已在GitHub上以MIT许可证公开，便于复现，并可集成到多尺度建模工作流中。