[논문 리뷰] Molecular Systems with Infinite and Finite Degrees of Freedom. Part I: Multi-Scale Analysis

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

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.

연구 동기 및 목표

• 유한 상태 마코프 체인(예: 단백질 구조 변화)과 무한 차원의 탄생-죽음 과정(예: 대사물질 농도)을 조합한 결합된 확률적 시스템을 엄밀히 분석하는 것.
• 트로터 근사를 사용하여 이러한 하이브리드 시스템의 마스터 방정식에 대한 연속체 근사를 도출함으로써 이산적 동역학에서 포커-플랑크 유형의 동역학으로의 전이를 가능하게 하는 것.
• 유한 상태 전이가 농도 변화보다 훨씬 높은 주파수로 발생한다는 가정 하에 형식적인 아디아바틱 이론을 개발하는 것.
• 고전적 효소 동역학, 특히 힐 유형의 포화 동역학을 미세한 확률적 동역학에서 아디아바틱 제거를 통해 재도출하는 것.

제안 방법

• 시간 진동의 공동 확률 과정을 유한 차원 및 무한 차원 동역학으로 번갈아 가며 분해하기 위해 트로터 근사를 사용한다.
• 유한 상태 마코프 체인의 무한소 생성자 행렬의 특이한 구조를 다루기 위해 드라진 역행렬을 적용한다. 특히 전이 행렬이 가역적이지 않은 경우에 유용하다.
• 투영 연산자(예: $\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$ 로 명시적으로 계산되며, 이는 아디아바틱 투영을 정의한다.
• 느린 변수(대사물질 농도)에 대한 유도된 효과적 동역학은 미세한 확률적 가정에서 유도된 결정론적 힐 유형 동역학 $K([A]) = \frac{K_1 [A]}{1 + K_2 [A]}$ (여기서 $n=1$) 를 재현한다.
• 이 이론은 효소 동역학에 대한 엄밀한 확률적 기반을 제공하며, 포화 행동이 빠른 결합/해리 사건들이 시간 평균화됨으로써 발생함을 보여준다.