[論文レビュー] On the steady states of weakly reversible chemical reaction networks

Jian Deng, Christopher R. Jones+2 more|arXiv (Cornell University)|Nov 10, 2011

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

ひとこと要約

本稿では、質量作用則に従う弱可逆な化学反応ネットワークが、各正の化学计量適合クラス（SCC）において少なくとも1つの正の平衡状態を有することを確立し、各SCCにおけるすべての平衡状態のインデックスの和が $(-1)^s$ に等しいというインデックス公式を提示する。ここで $s$ は化学计量部分空間の次元である。この結果は、このような系が常に平衡状態に到達するという実験的直観を裏付けるものである。

ABSTRACT

A natural condition on the structure of the underlying chemical reaction network, namely weak reversibility, is shown to guarantee the existence of an equilibrium (steady state) in each positive stoichiometric compatibility class for the associated mass-action system. Furthermore, an index formula is given for the set of equilibria in a given stoichiometric compatibility class.

研究の動機と目的

• 弱可逆な化学反応ネットワークが質量作用則に従う場合、各正の化学计量適合クラス（SCC）において常に少なくとも1つの正の定常状態を有するかどうかという長年の未解決問題を解消すること。
• 与えられたSCC内における平衡状態の集合に対して、トポロジー的インデックス公式を確立し、定常状態の数と安定性の種類をネットワークの構造的次元に関連付けること。
• 零の欠陊度および1の欠陊度を持つネットワークに関する先行結果を、より広い弱可逆ネットワークのクラスへ一般化し、構造的複雑さの変化に対しても平衡状態の存在が堅牢であることを確認すること。

提案手法

• 状態空間を分解し、濃度軌道のダイナミクスを解析するために、化学计量部分空間 $S$ 及びその直交補空間 $K^{ot}$ を用いる。
• 元のベクトル場を $K_2$ 上の新しいベクトル場に変換する微分同相写像 $\Psi: x_0 + S \to K_2$ を適用し、定常状態の構造とそのインデックスを保存する。
• 質量作用則のベクトル場表現を簡略化するために、縮小された構成行列 $\tilde{C}$ を用い、$K_2$ への射影に関して不変性を保証する。
• 次数論とホモトピー不変性を用いて、ベクトル場をより単純で同等のものに変形することで、与えられたSCC内におけるすべての定常状態のインデックスの和が $(-1)^s$ に等しいことを示す。ここで $s = \dim(S)$ である。
• ブロウアー次数およびホモトピーに関してインデックスの不変性といった代数的トポロジーの結果を応用し、インデックス公式を証明する。
• ネットワークの弱可逆性と速度定数の正の性質に依存し、ベクトル場が内向きを向く有界かつ凸な領域が存在することを保証し、少なくとも1つのゼロが存在することを保証する。

実験結果

リサーチクエスチョン

• RQ1 質量作用則に従う弱可逆な化学反応ネットワークは、各正の化学计量適合クラス（SCC）において常に少なくとも1つの正の定常状態を有するか？
• RQ2 このようなネットワークの与えられたSCC内におけるすべての定常状態のトポロジカルインデックスの和は何か？
• RQ3 インデックスの和は、特に化学计量部分空間の次元に関連したネットワークの構造的パrameterの関数として表現可能か？
• RQ4 反応ネットワークの構造、特に弱可逆性は、平衡状態の存在と分布にどのように影響するか？

主な発見

• 質量作用則に従う任意の弱可逆な化学反応ネットワークについて、空でない正の化学计量適合クラス（SCC）ごとに、少なくとも1つの正の定常状態が存在する。
• 与えられたSCCにおける定常状態の数は有限であり、それらのインデックスの和は $(-1)^s$ に等しい。ここで $s$ は化学计量部分空間の次元である。
• インデックスの和の公式 $\sum_{i=1}^{t} \text{ind}(x_i) = (-1)^s$ は、ネットワークが弱可逆である限り、特定の速度定数に依存せず、すべてのSCCに一様に成り立つ。
• 平衡状態の存在はトポロジカル不変性によって保証される：各SCCにおけるベクトル場は、次数が既知の場にホモトープであり、少なくとも1つのゼロが存在することを保証する。
• この結果は、実験的観察「通常の化学反応はしばしば定常状態に到達する」を裏付けるものであり、弱可逆系においてこの直観に厳密な基礎を与える。
• 証明は、状態空間の微分同相写像変換と、動的挙動の本質を保ちながら解析を簡略化する縮小された構成行列 $\tilde{C}$ の使用に依存する。