[論文レビュー] Topologically-constrained fluctuations and thermodynamics regulate nonequilibrium response

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

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.

研究の動機と目的

• 非平衡系における物理的可視量の定常応答の普遍的限界を特定すること。
• 既存理論が微視的状態間の位相的相関を無視する応答限界のギャップを埋めること。
• 状態空間の位相構造と熱力学的駆動を両方組み込んだ応答感度の制約フレームワークを構築すること。
• これらの限界がキネティクスの速度定数に依存せず、状態空間の接続性とサイクル駆動力にのみ依存することを示すこと。

提案手法

• 熱力学的整合性を満たすために有向遷移とその逆ペアを持つグラフ上の連続時間マコフ跳躍過程としてシステムをモデル化する。
• キルヒホッフの行列-木定理を活用して、状態空間のグラフ理論的手法を用いて定常分布およびその速度定数に関する微分を計算する。
• サイクルのアフィニティ（力）を用いて熱力学的駆動力を定義し、状態空間グラフ内の閉ループ回りの前向き・逆向き速度定数の積の対数比として計算する。
• 定常分布に対する制約付き最適化を用いて応答限界を導出し、その限界が可視量の値と状態空間の位相的分割に依存することを示す。
• 応答寄与を分離し極値構成を特定するために、状態を位相的に整合性のある部分集合に分割する手法を導入する。
• 受容体結合モデルにこの形式を適用し、最大感度が熱力学的駆動によって強化されたヒル関数に相当するヒル係数を持つこと、かつその上限に達することを示す。

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 異なる微視的状態における応答の相関関係は、系の感度に根本的な限界を設ける役割を果たすか？

主な発見

• 可視量の最大応答は、可視量の範囲と状態空間の位相的構造に依存する関数によって制限され、個々の速度定数に依存しない。
• 受容体結合モデルにおける感度の上限は、構造的平衡限界を超えて熱力学的駆動によって強化されたヒル係数を持つヒル関数に等しい。
• 最適な応答は、定常分布が2つの状態に集中する場合に達成され、最大感度は可視量の平均からのずれが特定の値にあるときに現れる。
• 到達可能な最大感度は、可視量の最小および最大の平均からのずれの積によって制限され、可能な値の領域がこれらの極端な値の比に臨界的比を許す場合に限界が達成される。
• 与えられた熱力学的駆動力のもとで、系の感度は状態空間の位相構造が応答への摂動に対する増幅を可能にするバランスの取れた状態分割を許す場合に最大化される。
• 受容体結合モデルを用いた数値的検証により、導出された限界がタイトであり、特定の速度定数設定下で物理的に達成可能であることが確認された。

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