[論文レビュー] Generalized Mass-Action Systems and Positive Solutions of Polynomial Equations with Real and Symbolic Exponents

Stefan C. Müller, Georg Regensburger|arXiv (Cornell University)|Jun 25, 2014

Protein Structure and Dynamics参考文献 48被引用数 24

ひとこと要約

本稿は、実数および記号的反応速度論的指数を有する一般化質量作用則系へと化学反応ネットワーク理論（CRNT）を拡張し、化学Stoichiometricおよび反応速度論的指数の部分空間における符号ベクトル条件を用いて、正の定常状態の存在および一意性を証明する。これはバーチの定理を一般化し、オriented matroid理論と二項方程式を用いた多定常状態のためのアルゴリズム的基準を提供する。

ABSTRACT

Dynamical systems arising from chemical reaction networks with mass action kinetics are the subject of chemical reaction network theory (CRNT). In particular, this theory provides statements about uniqueness, existence, and stability of positive steady states for all rate constants and initial conditions. In terms of the corresponding polynomial equations, the results guarantee uniqueness and existence of positive solutions for all positive parameters. We address a recent extension of CRNT, called generalized mass-action systems, where reaction rates are allowed to be power-laws in the concentrations. In particular, the (real) kinetic orders can differ from the (integer) stoichiometric coefficients. As with mass-action kinetics, complex balancing equilibria are determined by the graph Laplacian of the underlying network and can be characterized by binomial equations and parametrized by monomials. In algebraic terms, we focus on a constructive characterization of positive solutions of polynomial equations with real and symbolic exponents. Uniqueness and existence for all rate constants and initial conditions additionally depend on sign vectors of the stoichiometric and kinetic-order subspaces. This leads to a generalization of Birch's theorem, which is robust with respect to certain perturbations in the exponents. In this context, we discuss the occurrence of multiple complex balancing equilibria. We illustrate our results by a running example and provide a MAPLE worksheet with implementations of all algorithmic methods.

研究の動機と目的

• 実数および記号的反応速度論的指数を有する系へと、従来の質量作用則の枠組みを超えて化学反応ネットワーク理論（CRNT）を一般化すること。
• すべての速度定数および初期条件に対して、一般化質量作用則系における正の定常状態の存在および一意性の条件を確立すること。
• 定Stoichiometricおよび反応速度論的指数部分空間の符号ベクトルの交差を用いて、多定常状態の性質を特徴づけること。
• コンピュータ代数を用いたアルゴリズム的手法を提供し、走査例を通じてその有効性を示すこと。
• 記号的指数を有する一般化多項式系へとバーチの定理を拡張し、指数の摂動に対しても解の一意性が保たれるよう保証すること。

提案手法

• 反応速度論的係数が実数値をとる反応速度論的複合体を用いて、一般化質量作用則系を形式化する。これはStoichiometric係数とは明確に区別される。
• 反応ネットワークのグラフラプラシアンから導かれる二項方程式を用いて、複雑バランスの均衡状態をモデル化する。
• Stoichiometric部分空間 $S$ および反応速度論的指数部分空間 $ ilde{S}$ の符号ベクトルを用いて、正の解の存在および一意性を決定する。
• オriented matroid理論を用いて、Stoichiometricおよび反応速度論的指数行列のチロトープを比較し、符号ベクトルの整合性を保証する。
• 行列の小行列および符号計算を用いて、多定常状態のための条件 $\sigma(S) \cap \sigma(\tilde{S}^\perp) \neq \{0\}$ を検証するアルゴリズムを構築する。
• すべての手法を再現可能かつ実用的応用のためのMAPLEワークシートに実装する。

実験結果

リサーチクエスチョン

• RQ1 反応速度論的指数および速度定数がどのような条件下にあっても、一般化質量作用則系がすべてのStoichiometric適合クラスにおいて一意な正の定常状態を有するのか。
• RQ2 実数および記号的指数を有する多項方程式の正の解の存在および一意性が、どのようにアルゴリズム的に決定できるか。
• RQ3 一般化質量作用則系が複数の複雑バランスの均衡状態を有する場合、どのような構造的条件がこれを可能にするか。
• RQ4 記号的指数を有する系に対して、一般化バーチの定理はどのように適用され、どのような摂動が解の一意性を保つのか。
• RQ5 符号ベクトルおよびオriented matroid構造は、一般化化学反応ネットワークにおける多定常状態の特徴づけにおいて、どのように機能するのか。

主な発見

• 走査例において、$a, b, c > 0$ かつ $a < c$ のとき、Stoichiometricおよび反応速度論的指数行列のチロトープが一致するため、すべての正のStoichiometric適合クラスに一意な複雑バランスの均衡状態が存在する。
• $a > c$ かつ $a, b, c > 0$ のとき、交差 $\sigma(S) \cap \sigma(\tilde{S}^\perp) \neq \{0\}$ が成立するため、複数の複雑バランスの均衡状態が存在する可能性がある。
• $\sigma(S) = \sigma(\tilde{S})$ および $(+,\dots,+) \in \sigma(S^\perp)$ である条件が満たされれば、すべての速度定数に対して正の定常状態の存在および一意性が保証される。
• 符号ベクトルの部分空間が整合的である限り、指数の小さな摂動に対しても一般化バーチの定理がロバストに成立する。
• MAPLEワークシートに実装されたアルゴリズムフレームワークにより、一般化多項式系の正則性および多定常状態の条件を体系的に検証可能である。
• 一般化質量作用則系における多定常状態は、Stoichiometricおよび反応速度論的指数部分空間の符号ベクトルが非自明に交差する場合に生じ、質量作用則系の結果を一般化する。