[論文レビュー] Analysis and computation of some tumor growth models with nutrient: from cell density models to free boundary dynamics

Jian‐Guo Liu, Min Tang+2 more|arXiv (Cornell University)|Feb 2, 2018

Mathematical Biology Tumor Growth参考文献 14被引用数 4

ひとこと要約

本稿は栄養動態を伴う腫瘍成長モデルを分析し、圧力則 $ p(n) = n^\gamma $ に従う細胞密度モデルが $ \gamma \to \infty $ の極限において非圧縮状態としての自由境界 Hele-Shaw モデルにどのように対応するかを検討する。解析的推定を確立し、腫瘍フロントの伝播に関するベンチマーク解を導出し、保存則を満たし正値性を保つ数値スキームを用いて、細胞密度から自由境界ダイナミクスへの収束を検証する。

ABSTRACT

In this paper, we study the tumor growth equation along with various models for the nutrient component, including the \emph{in vitro} model and the \emph{in vivo} model. At the cell density level, the spatial availability of the tumor density $n$ is governed by the Darcy law via the pressure $p(n)=n^γ$. For finite $γ$, we prove some a priori estimates of the tumor growth model, such as boundedness of the nutrient density, and non-negativity and growth estimate of the tumor density. As $γ ightarrow \infty$, the cell density models formally converge to Hele-Shaw flow models, which determine the free boundary dynamics of the tumor tissue in the incompressible limit. We derive several analytical solutions to the Hele-Shaw flow models, which serve as benchmark solutions to the geometric motion of tumor front propagation. Finally, we apply a conservative and positivity preserving numerical scheme to the cell density models, with numerical results verifying the link between cell density models and the free boundary dynamical models.

研究の動機と目的

• 栄養依存的増殖を伴う腫瘍細胞密度モデルに対して、栄養の有界性、非負性、増殖制御のための厳密な解析的推定を確立すること。
• 圧力則 $ p(n) = n^\gamma $ を持つ細胞密度モデルが $ \gamma \to \infty $ の極限において、腫瘍組織の非圧縮状態を表す Hele-Shaw 型自由境界問題にどのように収束するかを形式化すること。
• 極限状態における Hele-Shaw モデルの解析的解を導出し、腫瘍フロント伝播のベンチマーク解としての役割を果たすこと。
• 細胞密度モデルから自由境界ダイナミクスへの遷移を正確に捉えることができる保存則を満たし、正値性を保つ数値スキームの開発と検証を行うこと。

提案手法

• 栄養消費が $ -\Delta c + \psi(n)c = 0 $ で記述される in vitro および in vivo の状況において、圧力 $ p(n) = n^\gamma $ に駆動される非線形拡散方程式により腫瘍細胞密度 $ n(x,t) $ をモデリングする。
• 漸近解析を用いて $ \gamma \to \infty $ の極限における Hele-Shaw 限界を導出し、$ p_\infty $ が $ D_\infty(t) $ 内で $ -\Delta p_\infty = G(c) $ を満たし、境界 $ \partial D_\infty(t) $ で $ p_\infty = 0 $ となる自由境界問題が得られることを示す。
• in vitro および in vivo モデルにおける球対称解を解析的に導出し、フロント伝播速度の正確なベンチマークを提供する。
• 有限の $ \gamma $ に対して細胞密度モデルを数値的に解くために、保存則を満たし正値性を保つ有限差分スキームを実装し、安定性と物理的整合性を確保する。
• 異なる $ \gamma $ におけるフロント伝播速度および空間的分布の比較を通じて、数値解が Hele-Shaw 限界に収束することを検証する。
• $ n $ および $ c $ に対してノイマンおよびディリクレ境界条件を用い、メッシュの細分化によりフロント追跡および栄養分布の精度を確保する。

実験結果

リサーチクエスチョン

• RQ1 栄養に影響を受けるモデルにおける圧力指数 $ \gamma $ が、腫瘍細胞密度の有界性、非負性、増殖推定にどのように影響するか。
• RQ2 圧力則 $ p(n) = n^\gamma $ を持つ細胞密度モデルと、$ \gamma \to \infty $ の極限における Hele-Shaw 自由境界モデルとの間の厳密な数学的関係は何か。
• RQ3 球対称な in vitro および in vivo の設定における Hele-Shaw モデルに解析的解は存在するか。また、それらはどのようにフロント伝播速度を予測するか。
• RQ4 保存則を満たし正値性を保つ数値スキームは、$ \gamma $ の増加に伴い、細胞密度ダイナミクスから自由境界運動への遷移をどの程度正確に再現できるか。

主な発見

• 有限の $ \gamma $ に対して、腫瘍細胞密度 $ n $ は非負であり、制御された割合で増加する。栄養濃度は一様に有界である。
• $ \gamma \to \infty $ の極限において、細胞密度モデルは形式的に Hele-Shaw 自由境界問題に収束し、腫瘍フロントは速度 $ v = -\nabla p_\infty $ で進化する。
• in vitro および in vivo モデルにおける球対称解の解析により、長時間にわたる極限において in vitro モデルのフロント伝播速度が in vivo モデルの2倍であることが確認された。
• $ \gamma = 80 $ の数値シミュレーションでは、細胞密度モデルと極限的な Hele-Shaw 解との間に良好な一致が得られ、フロントダイナミクスの収束が検証された。
• in vitro モデルではフロント伝播速度が $ R(t) \sim \sqrt{2t} $ と予測され、in vivo モデルでは $ R(t) \sim \sqrt{t} $ と予測され、数値結果と一致した。
• 数値スキームは正値性を正確に保持し質量保存を満たしており、高 $ \gamma $ に対しても安定かつ正確なフロント追跡が可能であった。ただし、$ \gamma $ が小さくかつ密度が1を超える場合、精度が低下した。