[論文レビュー] Molecular Systems with Infinite and Finite Degrees of Freedom. Part I: Multi-Scale Analysis

Luca Sbano, Markus Kirkilionis|ArXiv.org|Feb 28, 2008

Gene Regulatory Network Analysis参考文献 21被引用数 6

ひとこと要約

本稿は、有限および無限自由度を併せ持つ確率的分子系に対して、トロッター近似とドライジン逆行列を用いてマスター方程式の連続極限を導出し、厳密なマルチスケール漸近理論を構築する。高速に変化するマクロ分子（有限状態）が、ゆっくり変化する代謝物質濃度（無限自由度）を平均化するアディアバティック還元を確立し、古典的酵素動力学と整合する決定論的有効動力学を導出する。特に、ヒル型の飽和挙動を示す酵素動力学が、微視的確率的ダイナミクスから生じることを示す。

ABSTRACT

The paper analyses stochastic systems describing reacting molecular systems with a combination of two types of state spaces, a finite-dimensional, and an infinite dimenional part. As a typical situation consider the interaction of larger macro-molecules, finite and small in numbers per cell (like protein complexes), with smaller, very abundant molecules, for example metabolites. We study the construction of the continuum approximation of the associated Master Equation (ME) by using the Trotter approximation [27]. The continuum limit shows regimes where the finite degrees of freedom evolve faster than the infinite ones. Then we develop a rigourous asymptotic adiabatic theory upon the condition that the jump process arising from the finite degrees of freedom of the Markov Chain (MC, typically describing conformational changes of the macro-molecules) occurs with large frequency. In a second part of this work, the theory is applied to derive typical enzyme kinetics in an alternative way and interpretation within this framework.

研究の動機と目的

• 有限状態のマーカフ連鎖（例：タンパク質のコンformational 変化）と無限次元の誕生・死滅過程（例：代謝物質濃度）を組み合わせた結合確率的系を厳密に分析すること。
• トロッター近似を用いて、このようなハイブリッド系のマスター方程式の連続極限を導出し、離散的からフォッカー・プランク型ダイナミクスへの移行を可能にすること。
• 有限状態遷移が濃度変化に比べて高周波で発生するという仮定の下で、形式的なアディアバティック理論を構築すること。
• 特に、ヒル型の飽和動力学を含む古典的酵素動力学が、微視的確率的ダイナミクスからアディアバティック消去により再導出可能であることを示すこと。

提案手法

• 時間発展を交互に有限次元および無限次元のダイナミクスに分解するため、トロッター近似を用いる。
• 有限状態マーカフ連鎖における無限小生成行列の特異構造に対処するため、ドライジン逆行列を適用する。特に、遷移行列が非可逆である場合に有効である。
• 射影作用素（例：$\mathbf{\Pi}_\mu$）を用いて、確率分布を不変測度に沿った成分とその直交成分に分解する。
• 系の遅い変数を分離するために、周辺分布 $f(\mathbf{x}) = \mathbf{1}_\mu^T \rho(\mathbf{n})$ を導入する。
• 有限状態過程の不変測度に平行および直交する成分にベクトル場を分割することで、フォッカー・プランク方程式の弱形を導出する。
• アディアバティック極限において、特異な生成行列 ${{\mathcal{K}}^T}(\mathbf{x})$ の逆行列をドライジン逆行列 $({\mathcal{K}}^T_\mu)^D$ を用いて構成する。

実験結果

リサーチクエスチョン

• RQ1 有限および無限自由度を併せ持つハイブリッド確率的系が、連続極限においてどのように厳密に近似可能か。
• RQ2 有限状態遷移が濃度変化に比べて著しく高速に発生する場合に、有効マクロダイナミクスを導出するための数学的枠組みは何か。
• RQ3 ドライジン逆行列は、マルチスケールマーカフ過程における非可逆な無限小生成行列の解析をどのように支援するか。
• RQ4 ヒル型の飽和挙動を示すような古典的酵素動力学が、アディアバティック消去を経て微視的確率的モデルから再導出可能か。

主な発見

• トロッター近似により、結合確率的過程の時間離散化分解が有効に行え、連続極限の基礎が構築される。
• ドライジン逆行列 $({\mathcal{K}}^T_\mu)^D$ は、特異線形方程式 ${{\mathcal{K}}^T}(\mathbf{x}) \rho = \text{source}$ の解を提供し、アディアバティック極限の構築を可能にする。
• 確率分布 $\rho(\mathbf{x})$ は、${{\mathcal{K}}^T}(\mathbf{x})$ の核に属する成分とその直交成分に分解され、周辺分布 $f(\mathbf{x})$ が遅い時間スケールの進化を捉える。
• 有限状態マーカフ連鎖の不変測度 $\mu$ は明示的に $\mu = \left(\frac{k^{-}}{k^{-}+a k^{+}}, \frac{a k^{+}}{k^{-}+a k^{+}}\right)^T$ と計算され、アディアバティック射影を定義する。
• 遅い変数（代謝物質濃度）に対する導出された有効ダイナミクスは、$n=1$ のヒル型動力学の決定論的形 $K([A]) = \frac{K_1 [A]}{1 + K_2 [A]}$ を再現する。これは、微視的確率的仮定から導出可能である。
• 理論により、酵素動力学の厳密な確率的基盤が提供され、特に飽和挙動が、高速な結合・解離イベントの時間平均から生じることを示している。