[論文レビュー] Reachability in Biochemical Dynamical Systems by Quantitative Discrete Approximation

Luboš Brim, J. Fabriková+2 more|arXiv (Cornell University)|Jul 29, 2011

Gene Regulatory Network Analysis参考文献 20被引用数 8

ひとこと要約

本稿では、常微分方程式（ODE）で記述される生化学的動的システムに対して、パラメータ化された定量的離散近似（QDAA）手法を導入し、有界で有限状態の近似を用いて正確な到達可能性解析を可能にする。この手法は、局所的なベクトル場の挙動を組み込むことで、長方形抽象化における誤った遷移を低減し、従来の手法よりも到達可能な物質濃度の境界をより厳密に設定する。

ABSTRACT

In this paper a novel computational technique for finite discrete approximation of continuous dynamical systems suitable for a significant class of biochemical dynamical systems is introduced. The method is parameterized in order to affect the imposed level of approximation provided that with increasing parameter value the approximation converges to the original continuous system. By employing this approximation technique, we present algorithms solving the reachability problem for biochemical dynamical systems. The presented method and algorithms are evaluated on several exemplary biological models and on a real case study. This is a full version of the paper published in the proceedings of CompMod 2011.

研究の動機と目的

• 長方形抽象化における到達可能性解析の限界を解決すること。これは、長方形境界におけるベクトル場の過剰近似により、過剰な誤った挙動が生じるためである。
• パラメータを増加させることで元の連続的システムに収束する、パラメータ化された離散近似技術を開発すること。これにより、スケーラブルかつ正確な到達可能性計算が可能になる。
• 初期条件やパrameterの不確実性下での状態空間の到達可能性のグローバル解析を可能にする、数値シミュレーションの計算効率の良い代替手段を提供すること。
• 近似品質の尺度を用いて、モデル検査における誤った反例の程度を定量化・低減すること。
• 今後のシステム生物学用モデル検査ツールにこの手法を統合し、複雑で高次元の動的システムの解析を強化すること。

提案手法

• 連続的状態空間を長方形に分割することで有限離散近似を構築するが、境界内および境界間のベクトル場挙動に関する定量的測度を追加して強化する。
• 近似のレベルを制御するパラメータ化された精緻化プロセスを導入し、パrameterが増加するにつれて元のODEシステムに収束することを保証する。
• 到達可能な集合の精緻化に後向きシミュレーションを用い、パーティション密度の増加なしに物質濃度の区間境界の精度を向上させる。
• 長方形遷移系抽象化を適用するが、遷移ルールを局所的なベクトル場ダイナミクスに適合させることで、誤った遷移を低減する。
• アルゴリズム2を用いて、前方および後向きシミュレーションを統合し、状態変数の境界を厳密にすることで、到達可能な集合の過剰近似を計算する。
• プロトタイプソフトウェア実装により、小規模モデルおよび *E. coli* のアンモニウム輸送の実際の事例に対して手法を評価し、数値シミュレーションおよび長方形抽象化と結果を比較した。

実験結果

リサーチクエスチョン

• RQ1 生化学的ODEシステムの到達可能性解析において、誤った遷移を低減するパラメータ化された離散近似手法を開発できるか？
• RQ2 提案されたQDAA手法は、標準的な長方形抽象化と比較して、到達可能な集合推定の精度をどのように向上させるか？
• RQ3 パーティション密度の増加なしに、後向きシミュレーションが到達区間をどの程度精緻化し、過剰近似を低減できるか？
• RQ4 本手法は、NH₄⁺in や NH₃in などの物質濃度に対して、従来の抽象化手法よりも厳密な境界を提供できるか。また、計算的に実行可能か？
• RQ5 QDAAの近似品質は、選択されたパラメータ値とどの程度相関するか。また、真の連続的挙動に収束するか？

主な発見

• アンモニウム輸送モデルにおいて、QDAAはNH₄⁺inの上限を5.4×10⁻⁴と推定した。これは長方形抽象化の推定値5.4×10⁻⁴と同程度に近いが、過剰近似が低減されているため、QDAAの結果の方が正確である。
• QDAAはNH₃inに対して[10⁻⁸, 1.1×10⁻⁶]の到達可能な区間を達成し、初期濃度の境界と一致しており、シミュレーションデータと整合的である。
• NH₄⁺inに関しては、QDAAの結果[10⁻⁶, 10⁻⁵]は、長方形抽象化の推定値よりも1桁の精度で数値シミュレーションに近づいており、精度の向上が明確に示された。
• 本手法は、特に長方形境界における局所的なベクトル場ダイナミクスを組み込むことで、パーティションの細分化なしに誤った挙動を効果的に低減した。
• アルゴリズム2における後向きシミュレーションは、NH₃in や NH₄⁺in のような初期濃度が限定的な物質において、境界の精緻化を顕著に改善した。
• プロトタイプ実装により、QDAAが実際の生物学的モデルに効果的に適用可能であることが示された。形式的検証やシステム生物学的解析に適した信頼性の高い過剰近似を提供できる。