[論文レビュー] Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class

Santiago Akle, Onkar Dalal+4 more|arXiv (Cornell University)|May 12, 2011

Microbial Metabolic Engineering and Bioproduction参考文献 11被引用数 4

ひとこと要約

本稿は、反応ネットワークの欠陥や速度定数にかかわらず、質量保存型化学反応ネットワークが単一の端末結合クラスを有する場合、質量作用速度則のもとで常に正の定常状態を有することを証明している。著者は、弱く再帰的かつ均一な系に対して収束する凸最適化定式化と固定点アルゴリズムを導入し、フラックスおよび濃度ベクトルの反復的精錬により、信頼性の高い定常状態の計算を可能にした。

ABSTRACT

We establish that mass conserving single terminal-linkage networks of chemical reactions admit positive steady states regardless of network deficiency and the choice of reaction rate constants. This result holds for closed systems without material exchange across the boundary, as well as for open systems with material exchange at rates that satisfy a simple sufficient and necessary condition. Our proof uses a fixed point of a novel convex optimization formulation to find the steady state behavior of chemical reaction networks that satisfy the law of mass-action kinetics. A fixed point iteration can be used to compute these steady states, and we show that it converges for weakly reversible homogeneous systems. We report the results of our algorithm on numerical experiments.

研究の動機と目的

• 質量保存型化学反応ネットワークにおける正の定常状態の存在について、必要十分条件を確立すること。
• 速度定数に依存しない、計算的に扱いやすい方法で正の定常状態を同定すること。
• 閉系からの結果を、外部物質交換を伴う開放系へ、交換レートに関する簡単な条件下で拡張すること。
• 弱く再帰的かつ均一なネットワークにおける定常状態を計算するための固定点アルゴリズムの設計および検証すること。
• アルゴリズムの収束特性および多様なネットワークトポロジーやサイズにおけるスケーラビリティを調査すること。

提案手法

• 新規の凸最適化フレームワークを用いて、フラックスおよび濃度ベクトル上の固定点反復として定常状態問題を定式化する。
• 対数形式での質量作用条件を定義する：$ Y^T \log(c^*\,\!)=\log(v^*) $、これにより種の濃度と反応フラックスを関連付ける。
• フラックスバランス条件 $ YA_k v^* = b $ と質量作用条件 $ \psi(c^*) = v^* $ を用いて、定常状態を特徴付ける。
• 質量保存および所望の端末結合クラス構造を満たす制約を課したErdős-Rényiグラフを用いてランダムネットワークを構築する。
• ストイキオメトリック行列 $ Y $ を正規化し、$ Y^T \mathbf{1} = \mathbf{1} $ とすることで、$ e^T Y A_k = 0 $ を満たし、質量保存を保証する。
• アルゴリズム1を実装し、$ v_k $ および $ c_k $ を反復的に更新し、$ \|YA_k v_k\|_\infty $ の不適合性を最小化することで収束まで進める。

実験結果

リサーチクエスチョン

• RQ1 質量保存型化学反応ネットワークが単一の端末結合クラスを有する場合、どのような条件下で正の定常状態が存在するか？
• RQ2 弱く再帰的かつ均一なネットワークに対して、固定点アルゴリズムが信頼性を持って正の定常状態を計算できるか？
• RQ3 アルゴリズムの収束特性は、ネットワークサイズやトポロジーにどのように依存するか？
• RQ4 複数の端末結合クラスを有するネットワークに対しても、正の定常状態の存在が拡張可能か？
• RQ5 外部物質交換レート $ b $ は、正の定常状態の存在および計算可能性にどのような役割を果たすか？

主な発見

• 欠陥や速度定数にかかわらず、すべての質量保存型で単一の端末結合クラスを有するネットワークは、正の定常状態を有する。
• 弱く再帰的かつ均一な系に対して、固定点アルゴリズムは単調に収束し、フラックスバランス条件の不適合性が各反復で減少する。
• 50種と500複合体を有するネットワークにおいて、ネットワークサイズが50倍に拡大されても反復回数は10%未満で増加する。
• 平均反復回数はネットワークサイズにかかわらず安定しており、大規模系においても強固なスケーラビリティを示している。
• 数値実験の結果、複数の端末結合クラスを有するネットワークに対しても正の定常状態が存在する可能性が示唆されるが、理論的証明は未解決のままである。
• 凸最適化定式化により、速度定数やネットワーク欠陥の事前知識がなくても、直接的に定常状態を計算可能である。