[논문 리뷰] Ising Models of Cooperativity in Muscle Contraction

Elaheh Saadat, Matthieu Caruel외 7명|arXiv (Cornell University)|2026. 03. 04.

Cardiomyopathy and Myosin Studies인용 수 0

한 줄 요약

이 논문은 근육 수축 중 얇은 필라멘트 활성화의 협동성을 설명하기 위해 1차원 아이징 모델을 개발하고, Ca2+ 농도와 마이오신 모터 힘을 Hill 유사 응답과 연결하며 OM(Omecamtiv Mecarbil)을 포함한 실험과의 예측을 검증합니다.

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)의 영향을 평가한다.
• 기존 다상 모델과의 정량적 비교를 제공하고 상관길이 예측을 도출한다.
• 실험적 force-pCa 데이터에 대해 모델을 보정하고 해석 가능한 매개변수(J 및 [Ca2+]50)를 추출한다.

제안 방법

• 얇은 필라멘트를 최근접 이웃 결합 J와 Ca2+ 농도를 포착하는 외부장 h를 갖는 N개의 이진 상태 단위(스핀 s_i = ±1)의 1D 체인으로 모델링한다. H1 = -sum_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]}를 얻는다.
• 25°C에서 OM 없이 [Ca2+]50 및 n_H를 추출하도록 모델을 실험적 force-pCa 데이터에 보정한다; J ≈ 0.57 (n_H ≈ 3.14)로 보고, 25°C에서 OM 없이도 12–35°C 전반에서 유사한 적합을 얻는다.
• 상관길이 xi = [log((λ_+)/(λ_-))]^{-1}를 c=1에서 xi = [log((n_H+1)/(n_H-1))]^{-1}로 분석한다.
• Rice 등의 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+ 농도와 attached-motor 힘이 얇은 필라멘트 전체의 협동 활성화를 어떻게 조절하는가?
• RQ2 단층 Ising(두 상태) 모델이 실험적으로 관찰된 force-pCa 곡선과 온도/OM 의존성을 재현하는가?
• RQ3 Hill 계수 n_H와 Ising 결합 J 간의 관계는 무엇이며 이것이 상관길이와 어떤 관련이 있는가?
• RQ4 모델에 따르면 Omecamtiv Mecarbil(OM)이 협동성 및 힘 생성에 어떤 영향을 미치는가?
• RQ5 단층 Ising 모델이 Rice 등 4상 모델과 force-pCa 데이터를 설명하는 데 얼마나 비교되는가?

주요 결과

• 단층 Ising 모델은 f = 1/2{1 + (c-1)/sqrt[(c-1)^2 + 4 n_H^{-2} c]} 및 bijection n_H = exp(2J)로 force-pCa 관계를 재현한다.
• 보정된 매개변수는 25°C에서 OM 없이 [Ca2+]50 ≈ 10^(-6.5) M 및 n_H ≈ 3.14 (J ≈ 0.57)로 도출되며, OM 없이도 12–35°C에서 유사한 적합을 보인다.
• Hill 계수 n_H는 모터 힘 F0가 증가함에 따라 증가하여 협동성이 증가하는 것을 보여준다; OM을 적용하면 n_H가 1 미만으로 떨어져 협조성(협동성) 감소(J<0)를 나타낸다.
• 상관길이 xi는 F0가 증가함에 따라 커지며 온도가 달라짐에 따라 대략 2에서 6~7 단위까지 확장되어 고힘에서 더 긴 거리의 협동 확산을 시사한다.
• 모델의 force-pCa 예측은 실험 데이터 및 OM 효과와 일치하여 모터 힘에 의해 활성화가 확산되는 협동 메커니즘을 뒷받침한다.
• Rice 등 4상 모델에 비해 이 2상 모델은 force-pCa 데이터에 대해 과도한 매개변수 구성 없이도 가능하며, 고 Ca2+와 저 Ca2+ 데이터를 하나의 매개변수 집합으로 잘 맞춘다는 점에서 비교적 단순하다.