[論文レビュー] On the quasi-steady-state approximation in an open Michaelis--Menten reaction mechanism

Justin Eilertsen, Marc R. Roussel+2 more|arXiv (Cornell University)|Jan 14, 2021

Nonlinear Dynamics and Pattern Formation参考文献 50被引用数 18

ひとこと要約

この論文は、基質供給があるオープンなマイケリス・メンテン機構における準定常状態近似（QSSA）を、特異摂動論および不変多様体論を用いて厳密に分析している。特徴的な遅い多様体を同定し、標準的な特異摂動仮定を越えてQSSAの有効性を確立した。臨界多様体が存在しない場合でも、遅い多様体が存在し、軌道を引きつけることが示され、オープンな生化学的系におけるQSSAの理論的基盤を、正確な時間スケール推定と高次補正を伴って拡張した。

ABSTRACT

The conditions for the validity of the standard quasi-steady-state approximation in the Michaelis--Menten mechanism in a closed reaction vessel have been well studied, but much less so the conditions for the validity of this approximation for the system with substrate inflow. We analyze quasi-steady-state scenarios for the open system attributable to singular perturbations, as well as less restrictive conditions. For both settings we obtain distinguished invariant slow manifolds and time scale estimates, and we highlight the special role of singular perturbation parameters in higher order approximations of slow manifolds. We close the paper with a discussion of distinguished invariant manifolds in the global phase portrait.

研究の動機と目的

• 閉系に限らない準定常状態近似（QSSA）の理論的裏付けを、基質供給があるオープンな反応機構にまで拡張すること。
• 標準的な特異摂動仮定が成り立たない場合の、オープンなマイケリス・メンテン系における不変遅い多様体の存在と構造を調査すること。
• 時間スケール推定と遅い多様体に対する高次補正を導出することで、QSSA低次元化の精度を定量化すること。
• Poincaré球を用いた無限遠における振る舞いの分析を含め、オープン系のグローバル位相図を解析し、平衡点の安定性を特徴づけること。

提案手法

• オープンなマイケリス・メンテン機構に特異摂動論を適用し、速い酵素-基質結合と遅い生成物形成に関連する小さなパラメータ ε を特定する。
• Tikhonov–Fenichelパラメータ値（TFPV）を用いて、非孤立な平衡点を形成するパラメータ領域を特定し、臨界多様体を構成する。
• 正規形理論および中心多様体還元を用いて、臨界多様体が存在しない場合の遅い多様体の高次近似を計算する。
• Poincaré球を用いてグローバルな力学を解析し、無限遠における振る舞いを評価し、Poincaré変換を用いて無限遠における平衡点を分類する。
• Butler-McGehee定理およびPoincaré-Bendixson理論を用いて、ω-極限集合および平衡点のグローバル吸引性を研究する。
• 明示的な正規形展開を3次および4次まで導出し、不変多様体の安定性および幾何学的構造を特徴づける。

実験結果

リサーチクエスチョン

• RQ1 標準的な特異摂動フレームワークが適用できない場合、基質供給があるオープンなマイケリス・メンテン機構における準定常状態近似（QSSA）が有効となる条件は何か？
• RQ2 臨界多様体が存在しない場合でも、特徴的な不変遅い多様体が依然として存在し、軌道を引きつけることができるか？
• RQ3 標準的なTikhonov-Fenichel還元が適用できない場合、遅い多様体を近似するための高次補正の役割は何か？
• RQ4 無限遠における振る舞いを含むオープン系のグローバル位相図は、QSSAの有効性および精度にどのように影響するか？
• RQ5 遅い多様体の正確な幾何学的・力学的構造は何か？また、その構造は平衡点の安定性およびω-極限集合とどのように関係するか？

主な発見

• 標準的な特異摂動フレームワークが成立しない場合でも、ある種の一般性条件が満たされていれば、オープンなマイケリス・メンテン系は特徴的な不変遅い多様体を有する。
• 臨界多様体が存在しない場合でも、遅い多様体は依然として存在し、近傍の軌道を引きつけるため、QSSAは古典的Tikhonov-Fenichel理論を越えて有効である可能性を示している。
• 原点における正規形展開を3次まで行うと、˙y3 = (k2eT − k0)y3³ + ... が得られ、原点の安定性は (k2eT − k0) の符号に依存する。k2eT = k0 のとき、退化した吸引ノードとなる。
• 無限遠では、P2にサドルノードが存在し、上半球面に反発ノード部分を有する。原点の不安定多様体に属する軌道のω-極限集合は、P1またはP3の反対点を含むが、反発ノードを除く。
• k2eT = k0 の場合、4次までの正規形は ˙y3 = −(k−1 + k2)k2eT/k1 y3⁴ + ... となる。これにより、P1が退化した吸引ノードであり、第1象限のすべての初期条件に対してグローバルに吸引されることを確認した。
• グローバル位相図は完全に特徴づけられた：k2eT < k0 のとき、第1象限のすべての軌道は原点P0に収束し、k2eT > k0 のときP1に収束する。後者の条件下ではP1がグローバルに吸引的である。