[論文レビュー] 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の間で数学的に整合性のある結合手法を提示する。平均場理論を用いて開放系および2次反応を扱い、正確でスケーラブルなスケール間のフラックス交換を可能にすることで、開放的な生化学的系のシミュレーションを実現する。

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を数学的に整合性を持って結合するためのフレームワークを構築すること。
• 粒子数が保存されない開放系をモデル化する課題に取り組み、グランドカノニカル集合のアプローチを必要とすること。
• 特に2次反応に対して、微視的反応速度と巨視的PDEパラメータの正しい関係を確立すること。
• PBRDシミュレーションとPDEリザボアの間で数値的一致性を確保し、複雑な生化学的系の正確なマルチスケールモデリングを可能にすること。

提案手法

• 拡散粒子系の開放系における平均場理論を導出し、リザボアとの粒子交換を記述する。
• 2次反応を含むPBRD系の平均場定式化を構築し、微視的および巨視的反応定数の関係を明示する。
• 導出された平均場関係を基に、PBRDシミュレーションとPDEリザボアの間で整合性のある結合スキームを構築する。
• PDEの解法にCrank-Nicolson有限差分法を適用し、空間的・時間的両方で2次精度を達成する。
• 反応項と拡散項を処理するためにStrang分割法を適用し、数値的安定性と精度を維持する。
• PDEリザボアがPBRDシミュレーションにおいて正しく粒子の流入・流出を媒介することを確認することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 時間的・空間的に変化するリザボアを表す巨視的反応拡散PDEと、粒子基盤反応拡散シミュレーションをどのように一貫して結合できるか？
• RQ2 特に2次反応において、PBRDの微視的反応速度とPDEの巨視的反応定数との正しい平均場的関係は何か？
• RQ3 グランドカノニカル枠組み内でのPBRD系における粒子数のゆらぎを、どのように適切にモデル化できるか？
• RQ4 PBRDシミュレーションとPDEリザボアを複数スケールにわたって結合する際、一貫性と安定性を保証する数値スキームは何か？
• RQ5 空間的に不均一で時間的に変化するリザボアを組み込むことで、PBRDシミュレーションの精度とダイナミクスにどのような影響が生じるか？

主な発見

• 著者らは、拡散粒子系の開放系における平均場理論を導出し、グランドカノニカル集合におけるリザボアとの粒子交換を一貫してモデル化可能であることを示した。
• PBRDにおける微視的2次反応速度とPDEの巨視的反応定数との正しい平均場的関係を確立した。特に、拡散係数が反応速度関係に寄与しないことが判明し、これまでは誤って考慮されてきた点が重要な是正となった。
• 提案された結合スキームにより、PBRDシミュレーションとPDEリザボアの間で数学的整合性が保たれ、空間的・時間的に変化する境界条件を有するオープン生化学的系の正確なシミュレーションが可能となった。
• 数値的妥当性検証により、Crank-Nicolson法にStrang分割法を組み合わせた手法が、空間的・時間的両方で2次精度を達成し、安定で信頼性の高いPDE解を得られることを確認した。
• 本手法により、リザボア濃度が空間的・時間的に変化する動的extracellular環境を伴う細胞シグナリングなどの複雑で現実的な生物学的シナリオのシミュレーションが可能となった。
• 実装はGitHub上でMITライセンスのもとで公開されており、再現性の確保およびマルチスケールモデリングワークフローへの統合を容易にしている。