[論文レビュー] Maximal and minimal realizations of reaction kinetic systems: computation and properties

Gábor Szederkényi, Katalin M. Hangos+1 more|arXiv (Cornell University)|May 17, 2010

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

ひとこと要約

本稿では、与えられた複合体をもつ反応キネティクス系の最大および最小実現を計算するための混合整数線形計画法（MILP）に基づく手法を提案する。密度が最も高い実現の一意性を証明し、最小複雑性および完全に可逆な構造の計算を可能にする。主な貢献は、質量作用則に基づくCRN構造の最適化のための数値的に安定したフレームワークであり、欠損度解析および動的同等性に応用可能である。

ABSTRACT

This paper presents new results about the optimization based generation of chemical reaction networks (CRNs) of higher deficiency. Firstly, it is shown that the graph structure of the realization containing the maximal number of reactions is unique if the set of possible complexes is fixed. Secondly, a mixed integer programming based numerical procedure is given for computing a realization containing the minimal/maximal number of complexes. Moreover, the linear inequalities corresponding to full reversibility of the CRN realization are also described. The theoretical results are illustrated on meaningful examples.

研究の動機と目的

• 構造的に異なるが、同じ動的挙動を示す反応ネットワーク（CRN）を同定する長年の問題に取り組む。
• 複合体の集合が固定された場合に、反応数が最大の実現（密度が最も高い実現）の一意性を確立する。
• 最小の複雑度を持つ実現を計算的に効率的に生成するためのMILPフレームワークを開発する。
• 完全に可逆なCRN実現を計算可能にする線形制約を導出する。
• 欠損度解析および動的同等性の検証を含む、例題を用いて手法を示す。

提案手法

• 反応の存在を示すバイナリ変数と、速度定数を表す連続変数を用いて、実現問題を混合整数線形計画問題（MILP）として定式化する。
• 動的同等性を保証するため、化学 Stoichiometric 行列 Y と Kirchhoff 行列 Ak を用いる：M = Y · Ak。
• 質量作用則の制約と速度定数の非負性を確保する制約を課す。
• 反応の存在および可逆性をモデル化するため、ビッグM法を用いた論理的制約を適用する。
• 数値的安定性を確保し退化解を回避するため、パラメータ ε, ε₂, γ を導入する。
• Ak の非ゼロ要素数を最大化／最小化することで、それぞれ密度が最も高いおよび最も疎な実現を計算する。

実験結果

リサーチクエスチョン

• RQ1 複合体の集合が固定された場合、反応数が最大の実現は一意か？
• RQ2 与えられた動的モデルに対して、最小の複雑度を持つ反応ネットワークを効率的に計算できるか？
• RQ3 完全に可逆なCRN実現を保証するための必要十分な線形制約は何か？
• RQ4 MILPを用いて、特定の構造的性質を持つ動的同等実現をどのように計算できるか？
• RQ5 欠損度がゼロで可逆な実現の存在を検証し、その結果をシステムレベルの性質の推論に応用できるか？

主な発見

• 与えられた複合体の集合に対して、密度が最も高い実現（反応数が最大）は一意であり、標準的な基準構造を提供する。
• 3つの複雑度と3つの反応を持つ最小複雑度の実現が計算され、欠損度0を達成した。これにより、その場合の構造的唯一性が確認された。
• 9つの複雑度と8つの反応を持つ完全に可逆な実現が生成され、欠損度1の構造を示した。
• 3つの複雑度と3つの反応を持つ計算された実現は、欠損度ゼロ定理を満たし、複雑度安定性および散逸ハミルトニアン構造を示した。
• 元の不可逆ネットワーク（欠損度4）に対して、動的同等で、欠損度ゼロで可逆な実現が成功裏に同定された。
• 数値実験により、適切なパrameter化（ε=10⁻⁸, γ=0.01）を用いたMILPアプローチが、安定的かつ正確な実現をもたらすことが確認された。