[論文レビュー] Ising Models of Cooperativity in Muscle Contraction

Elaheh Saadat, Matthieu Caruel+7 more|arXiv (Cornell University)|Mar 4, 2026

Cardiomyopathy and Myosin Studies被引用数 0

ひとこと要約

論文は、一重のイジングモデルを用いて筋収縮時の薄いフィラメント活性化の協同性を説明し、Ca2+濃度とミオシンモータ力をHill様応答に結び付け、OMを含む実験との予測検証を行う。

ABSTRACT

Regulation of contraction in striated muscle is controlled by a dual mechanism involving both thin filaments containing actin and thick filaments containing myosin. The thin filament is activated by calcium ions binding to troponin, leading to tropomyosin azimuthal displacement which allows the activation of a regulatory unit (composed of one troponin, one tropomyosin and seven actin monomers) that exposes the actin sites for interaction with the myosin motors. Motor attachment to actin contributes to spreading activation within and beyond a regulatory unit along the thin filament through a cooperative mechanism. We introduce a one-dimensional Ising model to elucidate the mechanism of cooperativity in thin filament activation in relation to the force generated by the attached myosin motor. The model characterizes thin filament activation and cooperativity using only two parameters: one related to calcium concentration and the other to the force exerted by the attached myosin motor, which is modulated by temperature. At any force, the model is able to determine the extent of actin-myosin interactions on a correlation length ranging from two to seven actin monomers in addition to the seven actin monomers of the regulatory unit. Our theoretical predictions are successfully tested on experimental data, and our tests also include the condition of hindered filament activation by the use of the specific drug Omecamtiv Mecarbil (OM). According to our model, the effect of OM results in an anti-cooperativity mechanism accounting for the experimental data.

研究の動機と目的

• Ca2+とモータ生成力に応答した薄いフィラメント活性化が協同性を示す仕組みを説明する。
• Hill係数と力-pCa関係へ対応する、最小限の機序的Isingモデルフレームワークを提供する。
• モータ力と温度が協同性に与える影響を定量化し、Omecamtiv Mecarbil（OM）の影響を評価する。
• 既存の多状態モデルとの定量的比較を提供し、相関長の予測を導出する。
• 実験の力-pCaデータとモデルを較正し、解釈可能なパラメータ（Jと[Ca2+]50）を抽出する。

提案手法

• 薄いフィラメントをN個の二状態ユニットの1次元チェーンとしてモデル化（スピン s_i = ±1）で近接相互作用JとCa2+濃度を捉える外場hを用い、H1 = -総和_i (J s_i s_{i+1} + h s_i)。
• Isingモデルを詳細平衡遷移率k_±^iを用いる二状態マルコフ過程へ写像し、k_+^i/k_-^i = exp{2[J(s_{i+1}+s_{i-1})+h]}を得る。
• hとJをCa2+濃度とHill係数に関連づけ、h = (1/2) log(c)および J = (1/2) log(n_H)。
• 伝播行列法を用いて力f = (1+<s>)/2を算出し、f = 1/2 {1 + (c-1)/sqrt[(c-1)^2 + 4 n_H^{-2} c]}を得る。
• 実験の力-pCaデータへモデルを較正して[Ca2+]50とn_Hを抽出；25°CでOMなしの場合、J ≈ 0.57（n_H ≈ 3.14）を報告。
• 相関長ξ = [log((λ_+)/(λ_-))]^{-1}をc=1でξ = [log((n_H+1)/(n_H-1))]^{-1}として分析。
• Rice et al.の4状態モデルと比較し、それを2パラメータ類似形 f = 1/2 {1 + (x-1)/sqrt[(x-1)^2 + 4 n_H^{-2} x]}へ再定式化。
• OMの影響を調べ、OMはn_Hを低下させ（しばしば1未満）、検証域では反協同性（J<0）を誘発する。

実験結果

リサーチクエスチョン

• RQ1 Ca2+濃度と付着モータ力が薄いフィラメントの協同性活性化をどのように調整するか？
• RQ2 単一層のIsing（二状態）モデルで実験的に観測される力-pCa曲線とその温度/OM依存性を再現できるか？
• RQ3 Hill係数n_HとIsing結合Jの関係はどのようで、相関長とどう結びつくか？
• RQ4 Omecamtiv Mecarbil（OM）はモデルにより協同性と力発生にどのような影響を与えるか？
• RQ5 単一層IsingモデルはRiceらの4状態モデルと比べて力-pCaデータを説明する際にどのような利点・制限があるか？

主な発見

• 単一層Isingモデルは力-pCa関係を f = 1/2{1 + (c-1)/sqrt[(c-1)^2 + 4 n_H^{-2} c]}で再現し、双射 n_H = exp(2J) を満たす。
• 較正されたパラメータは、25°CでOMなしの場合[Ca2+]50 ≈ 10^(-6.5) M、n_H ≈ 3.14（J ≈ 0.57）を与え、12–35°Cの範囲でも同様の適合を示す。
• Hill係数n_Hはモータ力F0の増加とともに増加し、協同性が高くなる；OMを用いるとn_Hは1未満に低下し、反協同性（J<0）を示す。
• 相関長ξはF0とともに増大し、温度が変化しても大体2～6–7単位程度の範囲で、力が大きいほど協同性の拡がりが長くなることを示す。
• モデルの力-pCa予測は実験データおよびOM効果と整合し、モータ力支援による活性化の拡散による協同メカニズムを支持する。
• Riceらの4状態モデルと比較すると、二状態モデルは力-pCaデータに対して過剰制約されず、同様の2パラメータ形へ再定式化可能である。一方、4状態モデルはより複雑で高Ca/低Caデータの単一パラメータセットで適合させるのが難しい。