[論文レビュー] On the Mathematical Structure of Balanced Chemical Reaction Networks Governed by Mass Action Kinetics

Arjan van der Schaft, Shodhan Rao+1 more|arXiv (Cornell University)|Oct 27, 2011

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

ひとこと要約

本稿では、複素体グラフの対称的重み付きラプラシアン行列を用いて、バランス型化学反応ネットワークにおける質量作用速度論を幾何的・compactに定式化する。これにより、直接的な熱力学的解釈が可能となり、平衡状態の簡単な特徴付け、安定性の解析、シュール補行列を用いた体系的なモデル低減が可能となり、熱力学的整合性とネットワーク構造を保ったままとなる。

ABSTRACT

Motivated by recent progress on the interplay between graph theory, dynamics, and systems theory, we revisit the analysis of chemical reaction networks described by mass action kinetics. For reaction networks possessing a thermodynamic equilibrium we derive a compact formulation exhibiting at the same time the structure of the complex graph and the stoichiometry of the network, and which admits a direct thermodynamical interpretation. This formulation allows us to easily characterize the set of equilibria and their stability properties. Furthermore, we develop a framework for interconnection of chemical reaction networks. Finally we discuss how the established framework leads to a new approach for model reduction.

研究の動機と目的

• 反応ネットワークの複素体グラフ構造と化学计量を明示的に露わにする、化学反応ネットワークにおける質量作用速度論のcompactで幾何的な定式化を開発すること。
• 詳細なバランスとウェグシュタイナー条件を満たすバランス型ネットワークにおける平衡状態と安定性を、対称的重み付きラプラシアン行列を用いて特徴付けること。
• 境界種の化学ポテンシャルとして入力と出力を定義することにより、開いた反応ネットワークへのフレームワークの拡張を図ること。
• ラプラシアン行列のシュール補行列を用いた体系的モデル低減により、質量作用速度論と熱力学的平衡状態を保持したまま、低次元化されたネットワークを得ること。
• 欠損度に基づく手法（ゼロ欠損度や弱可逆性など）とは対照的に、統一的で洞察に基づいた代替手法を提供すること。

提案手法

• 反応ネットワークを、複素体を頂点とし、反応を有向辺とする有向複素体グラフとして表現する。
• 反応速度定数から導かれる対称的重み付きラプラシアン行列を用いて、動的挙動を定式化し、熱力学的整合性を保証する。
• 種の動的挙動を、複素体化学计量行列と、濃度の平衡値に対する対数比にラプラシアンを作用させた積として表現する。
• ラプラシアン行列のシュール補行列を用いて中間的複素体を削除することでモデル低減を行い、低次元化されたネットワークを導出する。
• 境界種と出力を化学ポテンシャルとして定義することで、開いたネットワークの接続を可能とし、非孤立系への拡張を図る。
• 複素体グラフのグラフ論的構造を活用してネットワークを分解し、安定性を解析する。

実験結果

リサーチクエスチョン

• RQ1 質量作用速度論に従うバランス型化学反応ネットワークの動的挙動を、ネットワーク構造と熱力学的性質を明確に露わにする、compactで幾何的な形でどのように表現できるか？
• RQ2 対称的重み付きラプラシアン行列は、このようなネットワークにおける平衡状態の特徴付けと、漸近的安定性の保証において果たす役割は何か？
• RQ3 入力と出力として境界種の化学ポテンシャルを定義することにより、フレームワークを流入・流出を伴う開いた反応ネットワークにどのように拡張できるか？
• RQ4 複素体グラフのラプラシアン行列を用いて、熱力学的平衡状態と質量作用速度論を保持したまま、体系的に低次元モデルを導出できるか？
• RQ5 この定式化と、欠損度理論やミヘリス・メンテン低減といった古典的手法との関係は何か？

主な発見

• バランス型化学反応ネットワークの動的挙動は、対称的重み付きラプラシアン行列を用いてcompactに表現可能であり、直接的な熱力学的解釈が可能となり、平衡状態と安定性解析が簡素化される。
• 平衡状態集合は、初期条件に対応する熱力学的平衡状態として一意に特徴付けられ、導出された定式化のもとで漸近的収束が保証される。
• ラプラシアン行列のシュール補行列によるモデル低減により得られる低次元ネットワークは、依然としてバランス型質量作用速度論系であり、平衡状態が保持される。
• 低次元化されたネットワークでは、中間的でない種の平衡濃度は同一であり、低次元化された反応定数は $\kappa = \frac{\kappa_X \kappa_Y}{\kappa_X + \kappa_Y}$ で与えられ、並列抵抗と類似した形を取る。
• 境界種の化学ポテンシャルとして入力と出力を定義することにより、開いたネットワークの接続が可能となり、ネットワーク接続の物理的理論が構築可能となる。
• この定式化は、コンSENSUSダイナミクスとの強い類似性を示しており、ネットワーキングシステムやシステム生物学への広範な応用可能性を示唆する。