[論文レビュー] Model comparison for single particle tracking in biological fluids

Martin Lysy, Natesh S. Pillai+5 more|arXiv (Cornell University)|Jul 22, 2014

Fractional Differential Equations Solutions参考文献 73被引用数 3

ひとこと要約

本研究では、ヒト肺粘液内の1-µmトレーサー粒子の運動を評価するためにベイジアンモデル比較を用い、分数 Browm運動（fBM）および一般化ランジュビン方程式（GLE）モデルを検討した。データは、小規模から中程度の記憶を持つモデルを除外し、長記憶スペクトルを有するfBMまたはGLEを支持しており、単純なブラウン運動を超える粘弾性モデルの必要性を示している。

ABSTRACT

State-of-the-art techniques in particle tracking microscopy provide high-resolution path trajectories for the dynamical analysis of diverse foreign particles in biological fluids. For particles on the order of one micron in diameter, these paths are not generally consistent with simple Brownian motion. Despite an abundance of data confirming these findings, stochastic modeling of the complex particle motion has received comparatively little attention. Yet, prediction beyond experimental time scales -- such as first-passage times through biological barriers -- is a fundamental scientific objective that requires a careful statistical approach. Even among existing candidate models, few attempts have been made at likelihood-based model comparisons and other quantitative evaluations. In this article, we employ a robust Bayesian methodology to address this gap. We analyze two competing models -- fractional Brownian motion (fBM) and a generalized Langevin equation (GLE) consistent with viscoelastic theory -- applied to an ensemble of 30 second paths of 1 micron diameter tracer particles in human lung mucus. We conclude that either fBM or a GLE with a long memory spectrum are favored by the data, whereas GLEs with small to moderate memory spectra are insufficient.

研究の動機と目的

• 複雑な生物学的媒体における粒子運動のための確率的モデル同士の厳密な統計的比較が不足しているという問題に取り組むこと。
• ヒト肺粘液内での1-µm粒子の観測軌道をよりよく説明するのは、分数 Browm運動（fBM）モデルか、一般化ランジュビン方程式（GLE）モデルかを評価すること。
• 具体的には、記憶スペクトル長が実験データに最もよく適合するモデル構造は何かを特定すること。
• 最も適切な確率的モデルを選択することで、生物学的バリアを通過するまでの初回通過時間などの長時間ダイナミクスを信頼性高く予測すること。
• 単粒子追跡におけるモデル評価のための堅牢な尤度ベースのフレームワークを確立し、現在の生物物理学的モデリング実務における空白を埋めること。

提案手法

• 競合するモデルの周辺尤度（証拠）を計算するためのベイジアン推論フレームワークを採用し、直接的なモデル比較を可能にした。
• 高解像度粒子追跡顕微鏡を用いて取得したヒト肺粘液内での1-µmトレーサー粒子の30秒間軌道にこの手法を適用した。
• 主な2つのモデルを評価した：ハーストパラメータHで特徴づけられる分数 Browm運動（fBM）、および異なる記憶カーネルスペクトルを有する一般化ランジュビン方程式（GLEs）。
• GLEフレームワークを用いて、特に長記憶と小〜中程度の記憶スペクトルを有するモデルを比較した。
• 相対的なモデル支持力度を定量化するためにベイズ因子を計算し、より高い周辺尤度を示すモデルを優遇した。
• 事後予測チェックと軌道アンサンブルにおける交差検証を通じて、モデルの適合度と予測性能を評価した。

実験結果

リサーチクエスチョン

• RQ1 ヒト肺粘液内での1-µm粒子運動を説明するにあたり、分数 Browm運動か一般化ランジュビン方程式のどちらの確率的モデルがより適切か。
• RQ2 GLEモデルにおける記憶スペクトル長が、実験的粒子軌道へのモデル適合度にどのように影響するか。
• RQ3 尤度ベースの比較に基づいて、fBMモデルとGLEモデルのどちらが統計的に優れているかの証拠はあるか。
• RQ4 GLEフレームワークにおける小〜中程度の記憶スペクトルを有するモデルは、粘弾性を持つ生物学的媒体内での粒子ダイナミクスを適切に記述できるか。
• RQ5 初回通過時間のような長時間ダイナミクスを予測するにあたり、データは単純なブラウン運動よりも粘弾性モデルの使用をどの程度支持しているか。

主な発見

• 分数 Browm運動（fBM）は、ヒト肺粘液内での1-µm粒子運動を説明する有効なモデルとして、データによって強く支持されている。
• 長記憶スペクトルを有する一般化ランジュビン方程式（GLEs）も支持されており、持続的な記憶を有する粘弾性効果がダイナミクスを記述するために不可欠であることが示された。
• 小〜中程度の記憶スペクトルを有するGLEモデルは統計的に除外され、データに観察された長距離時間相関を捉えられていないことが示された。
• ベイジアンモデル比較フレームワークは、最も妥当性の高い確率的モデルを的確に同定できており、複雑な生物学的系における尤度ベースの評価の重要性を示している。
• 結果から、粘弾性が粘液内粒子輸送において重要な役割を果たしており、長時間の挙動を正確に予測するには長距離時間相関を考慮したモデルが必要であることが示唆された。
• 30秒間の軌道アンサンブルは、標準的なブラウン運動とは整合しない持続的運動を示しており、より複雑な確率的モデルの導入が不可欠であることがわかった。