[논문 리뷰] van't Hoff-Arrhenius Analysis of Mesoscopic and Macroscopic Dynamics of Simple Biochemical Systems: Stochastic vs. Nonlinear Bistabilities

Yunxin Zhang, Hao Ge외 1명|arXiv (Cornell University)|2010. 11. 11.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 3인용 수 4

한 줄 요약

이 논문은 미세구조 생화학 시스템에 대한 반트호프-아레니우스 분석을 도입하여 시스템 크기 $V$를 역온도에 해당하는 것으로 간주함으로써 비선형 이중안정성(결정론적, 엔탈피 기반)과 소음 유도 스토케스틱 이중안정성(엔트로피 장벽 기반)을 구분한다. 주요 기여는 $\Delta\Phi^\ddag(V) \approx \Delta\phi_0^\ddag + \Delta\phi_1^\ddag/V$의 분해로, 비선형 이중안정성의 경우 $\Delta\phi_0^\ddag > 0$이고, 스토케스틱 이중안정성의 경우 $\Delta\phi_0^\ddag < 0$이지만 $\Delta\phi_1^\ddag > 0$이므로, 후자의 경우 전이 속도가 $V$ 증가에 따라 증가하는 이유를 설명한다.

ABSTRACT

Multistability of mesoscopic, driven biochemical reaction systems has implications to a wide range of cellular processes. Using several simple models, we show that one class of bistable chemical systems has a deterministic counterpart in the nonlinear dynamics based on the Law of Mass Action, while another class, widely known as noise-induced stochastic bistability, does not. Observing the system's volume ($V$) playing a similar role as the inverse temperature ($β$) in classical rate theory, an van't Hoff-Arrhenius like analysis is introduced. In one-dimensional systems, a transition rate between two states, represented in terms of a barrier in the landscape for the dynamics $Φ(x,V)$, $k\propto\exp\{-VΔΦ^‡(V)\}$, can be understood from a decomposition $ΔΦ^‡(V) \approxΔϕ_0^‡ Δϕ_1^‡/V$. Nonlinear bistability means $Δϕ_0^‡&gt;0$ while stochastic bistability has $Δϕ_0^‡&lt;0$ but $Δϕ_1^‡&gt;0$. Stochastic bistabilities can be viewed as remants (or "ghosts) of nonlinear bifurcations or extinction phenomenon, and $Δϕ_0^‡$ and $Δϕ_1^‡$ as "enthalpic" and "entropic" barriers to a transition.

연구 동기 및 목표

• 미세구조 생화학 시스템에서 결정론적 비선형 이중안정성과 소음 유도 스토케스틱 이중안정성을 구분하기.
• 전통적 전이 상태 이론에서의 역온도에 해당하는 시스템 크기 $V$를 제어 변수로 사용하여 이론적 프레임워크를 개발하기.
• 장애물 함수 $\Phi(x,V)$로 분해된 엔탈피($\phi_0$) 및 엔트로피($\phi_1/V$) 기여를 바탕으로 안정 상태 간 전이 속도를 특성화하기.
• 스토케스틱 이중안정성에서 전이 속도가 $V$ 증가에 따라 증가하는 이유를 설명하고, 이는 비선형 이중안정성에서 관찰되는 감소 경향과 정반대임을 밝히기.
• 대규모 $V$ 점근적 분석과 리아푸노프 함수를 통해 화학 마스터 방정식과 매크로스코픽 ODE 역학 간의 연결을 수립하기.

제안 방법

• 유한 부피 $V$를 갖는 반응 시스템에서의 미세구조 스토케스틱 역학을 모델링하기 위해 화학 마스터 방정식(CME)과 질스피 알고리즘을 사용한다.
• 정적 분포 $p_n^{st} \propto e^{-V\Phi(n,V)}$를 유도하며, 여기서 $\Phi(n,V) = \phi_0(x) + \phi_1(x)/V + \cdots$, $x = n/V$이다.
• 대규모 $V$ 점근적 전개를 적용하여 매크로스코픽 ODE 시스템 $dx/dt = F(x)$에 대한 리아푸노프 함수로 $\phi_0(x)$를 추출한다.
• 전이 장벽의 분해 $\Delta\Phi^{\ddag}(V) \approx \Delta\phi_0^{\ddag} + \Delta\phi_1^{\ddag}/V$를 이용하여 $\Delta\phi_0^{\ddag}$ 및 $\Delta\phi_1^{\ddag}$의 부호에 따라 이중안정성 유형을 분류한다.
• 크라머스 이론과 라플라스 방법을 적용하여 평균 첫 번째 통과 시간 $T_{x_1^* \to x_2^*} \propto e^{V[\Phi(x^{\ddag},V) - \Phi(x_1^*,V)]}$를 유도하며, 전이 속도와 장애물 함수 간의 연관성을 규명한다.
• 세부 균형과 광역정렬 집합 모델을 사용하여 $V\phi_0(x) + \phi_1(x)$의 볼록성을 검증함으로써, 닫힌 시스템에서 단일 최대값 균형이 보장됨을 입증한다.

실험 결과

연구 질문

• RQ1 시스템 크기 $V$는 스토케스틱 생화학 시스템에서 안정 상태 간 전이 속도에 어떻게 영향을 미치는가?
• RQ2 부피 의존성 측면에서 비선형 이중안정성(결정론적)과 소음 유도 스토케스틱 이중안정성은 어떻게 다를까?
• RQ3 반트호프-아레니우스 유사 분석을 $V^{-1}$을 역온도에 해당하는 것으로 간주함으로써 생화학 시스템에 적용할 수 있는가?
• RQ4 장애물 함수의 엔탈피 기여와 엔트로피 기여를 특성화하는 데 $\phi_0(x)$와 $\phi_1(x)$의 역할은 무엇인가?
• RQ5 왜 스토케스틱 이중안정성에서는 전이 속도가 $V$ 증가에 따라 증가하는 반면, 비선형 이중안정성에서는 감소하는가?

주요 결과

• 비선형 이중안정성은 $\Delta\phi_0^{\ddag} > 0$로 특징지어지며, 이는 엔탈피 장벽이 시스템 크기 증가에 따라 증가함을 의미하고, 이로 인해 $V$ 증가에 따라 전이 속도가 감소한다.
• 스토케스틱 이중안정성은 $\Delta\phi_0^{\ddag} < 0$이지만 $\Delta\phi_1^{\ddag} > 0$일 때 발생하며, 이는 증가하는 $V$에 따라 감소하는 엔트로피 장벽을 의미하여 전이 속도가 증가함을 설명한다.
• 화학 마스터 방정식의 정적 분포는 점점이 $p^{st}_n(V) \approx e^{-V\phi_0(x) - \phi_1(x)}$ 형태를 취하며, 여기서 $\phi_0(x)$는 매크로스코픽 ODE 시스템에 대한 리아푸노프 함수이다.
• 안정 상태 간 평균 첫 번째 통과 시간은 $T \propto e^{V[\Phi(x^{\ddag},V) - \Phi(x_1^*,V)]}$로 스케일링되며, 이는 장애물 함수 $\Phi(x,V)$가 전이 동역학을 결정하는 데 핵심적임을 확인한다.
• 세부 균형을 갖는 시스템에서는 평형 분포가 단일 최대값을 가지며, $V\phi_0(x) + \phi_1(x)$는 볼록하여 안정성 보장과 자발적 이중최대값 방지가 가능하다.
• 열린 시스템에서의 멸종의 유령은 $\phi_1(x)$만으로는 감지할 수 없으며, 경계에서의 점근적 전개가 균일하게 유효하지 않기 때문이다.