[論文レビュー] Hunting $\varepsilon$: The origin and validity of quasi-steady-state reductions in enzyme kinetics

Justin Eilertsen, Malgorzata Tyczynska+1 more|arXiv (Cornell University)|Sep 8, 2021

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

ひとこと要約

本稿は、幾何的特異摂動理論と中心多様体解析を用いて、分子間自己触媒的ゼノゲン活性化反応機構（IAZA）における準定常状態近似（QSSA）の有効性を規定する無次元パラメータ ε を同定する。ε は速い時定数と遅い時定数の比から生じ、QSS還元は片方の極限ではフェニケル理論、もう片方の極限では中心多様体理論によって正当化され、その有効性はパラメータ空間内の経路に依存する——特に動的トランスクリティカル分岐点付近では古典的理論が失敗する。

ABSTRACT

The estimation of the kinetic parameters requires the careful design of experiments under a constrained set of conditions. Many estimates reported in the literature incorporate protocols that leverage simplified mathematical models known as quasi-steady-state reductions. Such reductions often - but not always - emerge as the result of a singular perturbation scenario. However, the utilization of the singular perturbation reduction method requires knowledge of a dimensionless parameter, $\varepsilon$, that is proportional to the ratio of the reaction's fast and slow timescales. Using techniques from differential equations, Fenichel theory, and center manifold theory, we derive the appropriate $\varepsilon$ whose magnitude regulates the validity of the quasi-steady-state reduction employed in the reported experimental procedures for intermolecular autocatalytic zymogen activation reaction. Although the model equations are two-dimensional, the fast/slow dynamics are rich. The phase plane exhibits a dynamic transcritical bifurcation point in a particular singular limit. The existence of such a bifurcation is relevant, because the critical manifold losses normal hyperbolicity and classical Fenichel theory is inapplicable. Furthermore, we show that in some cases chemical reversibility can be interpreted dynamically as an imperfection, since the presence of reversibility can destroy the bifurcation structure present in the singular limit. We show that the reduction method by which QSS reductions are justified can depend on the path taken in parameter space. Specifically, we show that the standard quasi-steady-state reduction for this reaction is justifiable by center manifold theory in one limit, and via Fenichel theory in a different limit.

研究の動機と目的

• IAZA反応機構における準定常状態還元の有効性を規定する正しい無次元パラメータ ε を特定すること。
• QSS近似を用いた実験的手順が正確な動力学的パラメータ推定をもたらす条件の曖昧さを解消すること。
• 古典的フェニケル理論が正規双曲性の喪失により失敗する特異極限におけるIAZA機構の動的分岐構造を分析すること。
• 化学的可逆性や初期条件の変化が存在する場合に、QSS還元が数学的に正当化される条件を明確にすること。
• QSS還元の正当化がパラメータ空間内の経路に依存することを示すこと——極限を取る順序に依存する。

提案手法

• IAZA機構における速い時定数と遅い時定数の比として適切な ε を導出するため、幾何的特異摂動理論（GSPT）を適用する。
• 正規双曲性が保たれる条件下で、速い時定数が遅い時定数に比べて著しく短い極限において、フェニケル理論を用いてQSS還元を正当化する。
• 臨界多様体の正規双曲性が失われる小全酵素濃度極限（ET → 0）において、中心多様体理論を用いてQSS還元を正当化する。
• 漸近的解析により、線形化系の固有値比の平方根に比例する距離だけ分岐点から離れた領域に軌道が保たれることを示す。
• 位相平面のダイナミクスを分析し、特異極限において動的トランスクリティカル分岐が存在することを特定する——これは標準的GSPTの適用を妨げる。
• 化学的可逆性の影響を比較し、可逆性が不正則な摂動として作用し、不可逆極限に存在する分岐構造を破壊する可能性があることを示す。

実験結果

リサーチクエスチョン

• RQ1 分子間自己触媒的ゼノゲン活性化機構における準定常状態近似の有効性を規定する正しい無次元パラメータ ε は何か？
• RQ2 特異極限において動的トランスクリティカル分岐が存在する場合、なぜ標準的準定常状態還元が厳密に正当化されないのか？
• RQ3 QSS還元の正当化は、パラメータ空間内の経路にどのように依存するのか？特に、k2 が小さいか ET が小さいかといった異なる極限に近づく場合に。
• RQ4 化学的可逆性はどのようにして分岐構造を変化させる動的不正則性として作用するのか？
• RQ5 中心多様体理論とフェニケル理論は、同じQSS還元を同時に正当化できるか？もしそうなら、それぞれが成立する異なるパラメータ極限は何か？

主な発見

• パラメータ ε は速い時定数と遅い時定数の比として導出され、IAZA機構におけるQSS還元の有効性にはその小ささ（ε ≪ 1）が必要不可欠である。
• 通常の実験的初期条件から出発する軌道は、線形化系の固有値比の平方根に比例する距離だけ動的トランスクリティカル分岐点から離れて保たれる。
• 正規双曲性の喪失により、古典的フェニケル理論は分岐点では適用不可能であるが、分岐点を直接通過する「偽のカナード」が存在する。
• 標準的QSS還元は、小全酵素濃度極限（ET → 0）では中心多様体理論によって、小k2極限ではフェニケル理論によって正当化される。
• 化学的可逆性は動的不正則性として作用し、分岐構造を破壊し、QSS近似の有効性を変化させる可能性がある。
• QSS還元の正当化は経路依存的である——同じ還元がパラメータ空間内の異なる極限を取るに従い、異なる数学的理論によって正当化されることがある。