[논문 리뷰] Topologically-constrained fluctuations and thermodynamics regulate nonequilibrium response

Gabriela Fernandes Martins, Jordan M. Horowitz|PubMed|2023. 05. 30.

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

한 줄 요약

이 논문은 상태공간의 위상적 제약 조건과 열역학적 구동력의 조합을 통해 마코프 시스템에서 비평형 반응의 기본 한계를 도출한다. 반응성의 관측 가능도—예를 들어 수용체 결합—는 상태공간의 위상적 및 열역학적 요소에 의해 강화된 힐 계수를 통해 최대화될 수 있으며, 상태공간의 구조 외에는 동역학적 속도 상수의 세부 정보가 필요하지 않다.

ABSTRACT

The limits on a system's response to external perturbations inform our understanding of how physical properties can be shaped by microscopic characteristics. Here, we derive constraints on the steady-state nonequilibrium response of physical observables in terms of the topology of the microscopic state space and the strength of thermodynamic driving. Notably, evaluation of these limits requires no kinetic information beyond the state-space structure. When applied to models of receptor binding, we find that sensitivity is bounded by the steepness of a Hill function with a Hill coefficient enhanced by the chemical driving beyond the structural equilibrium limit.

연구 동기 및 목표

• 비평형 시스템에서 물리적 관측 가능도의 정 steady-state 반응에 대한 보편적 한계를 규명하는 것.
• 기존 이론에서 미세 상태 간의 위상적 상관관계를 간과한 반응 한계 이론의 격차를 메우는 것.
• 상태공간의 위상적 특성과 열역학적 구동력을 동시에 통합하여 반응 민감도를 제약하는 프레임워크를 개발하는 것.
• 이러한 한계가 동역학적 속도 상수의 값에 의존하지 않고, 상태공간의 연결성과 순환력에만 의존함을 보여주는 것.

제안 방법

• 열역학적 일관성을 확보하기 위해 방향성 있는 전이와 그 역전이 쌍을 갖는 그래프 위의 연속시간 마코프 점프 과정으로 시스템을 모델링하는 것.
• 기저 그래프 이론적 방법을 사용하여 동역학적 속도 상수에 대한 정상분포 및 그 도함수를 계산하며, 키르히호프의 매트릭스-트리 정리를 활용하는 것.
• 폐쇄된 순환 경로에서의 전진 및 역방향 속도 곱의 로그 비율을 통해 순환 친화도를 정의함으로써 열역학적 구동력을 도입하는 것.
• 정상분포에 대한 제약 최적화를 통해 반응 한계를 도출하며, 이는 관측 가능도의 값과 상태공간의 위상적 분할에 의존하는 것.
• 반응 기여를 분리하고 극값 구성을 식별하기 위해 상태를 위상적으로 일관된 방식으로 부분집합으로 나누는 것.
• 수용체 결합 모델에 이 형식을 적용하여, 최대 민감도가 열역학적 구동력으로 인해 증가된 위상적 요소를 반영한 힐 함수로 제약됨을 보여주는 것.

Figure 1: Representative example: (a) Five-state house graph with perturbed edge $e_{12}$ highlighted in pink. (b) Cycles associated to the thermodynamic forces $F_{C}$ . Only cycles through the perturbed edge $e_{12}$ constrain the response. (c) Topologically-consistent splittings ${\mathcal{V}}^{1

실험 결과

연구 질문

• RQ1 상태공간의 위상적 특성이 비평형 정상 상태에서 물리적 관측 가능도의 최대 가능한 반응을 어떻게 제약하는가?
• RQ2 열역학적 구동력은 평형 한계를 초월해 시스템의 민감도를 얼마나 향상시킬 수 있으며, 이는 어떻게 정량화되는가?
• RQ3 개별 동역학적 속도 상수의 지식 없이도 상태공간의 구조만으로도 반응 한계를 도출할 수 있는가?
• RQ4 다른 미세 상태에서의 반응 상관관계가 시스템 민감도의 기본 한계를 설정하는 데 어떤 역할을 하는가?

주요 결과

• 관측 가능도의 최대 반응은 개별 속도 상수의 값과 무관하게 관측 가능도의 범위와 상태공간의 위상적 구조에 의존하는 함수로 제약된다.
• 수용체 결합 모델에서의 민감도 한계는 열역학적 구동력으로 인해 구조적 평형 한계를 초월해 증가된 힐 계수를 갖는 힐 함수와 동일한 형태이다.
• 최적의 반응은 정상분포가 두 상태에 집중될 때 달성되며, 최대 민감도는 관측 가능도의 평균에서의 편차가 특정 값일 때 발생한다.
• 도달 가능한 최대 민감도는 관측 가능도의 최소 및 최대 편차의 곱에 의해 제약되며, 가능한 값의 영역이 이러한 극단의 비율에 대해 임계 비율을 허용할 경우 한계가 달성된다.
• 주어진 열역학적 구동력 하에서, 상태공간 위상적 구조가 상태를 균형 있게 분할할 수 있을 경우 시스템의 민감도가 최대화된다.
• 수용체 결합 모델을 이용한 수치적 검증을 통해 유도된 한계가 특정 속도 구성 조건 하에서 날카롭고 물리적으로 실현 가능함을 확인하였다.

Figure 2: Response ratio $\mathcal{R}_{B}$ as a function of the maximum cycle force ${\mathcal{F}}_{\rm max}$ around cycles containing edge $e_{12}$ of the house graph (Fig. 1 ) for 15000 rate matrices with the logarithmic asymmetric $\ln(W_{ij}/W_{ji})$ and symmetric $\ln(W_{ij}W_{ji})$ parts sampl