시간에 따라 변하는 기질 강성에 대한 세포 반응의 히스테리시스

초록

본 연구는 수축성 액틴 번들 역학과 Rho‑신호전달 경로를 결합한 모델을 통해 기질 강성이 시간에 따라 변할 때 세포가 보여주는 이중안정성 및 히스테리시스 현상을 예측한다. 모델 분석 결과, 낮은 강도에서는 수축력이 형성되지 않고, 높은 강도에서는 지속적인 강한 수축 상태가 유지되는 두 개의 임계값이 존재한다. 또한 강성 변화 속도에 따라 세포의 트랙션 포스가 경로를 따라 다른 값을 갖는 속도 의존적 히스테리시스가 발생한다는 것을 제시한다.

상세 분석

이 논문은 세포가 외부 기질의 강성 변화를 감지하고 이에 따라 내부 수축력을 조절하는 메커니즘을 정량적으로 설명하기 위해, 두 개의 주요 서브시스템을 통합한 수학적 모델을 구축하였다. 첫 번째 서브시스템은 미세섬유와 근육성 액틴 번들의 기계적 거동을 묘사하며, 탄성 기질에 대한 접착점에서 발생하는 변형률과 그에 따른 내부 응력‑변형 관계를 선형 또는 비선형 스프링-댐퍼 형태로 구현한다. 두 번째 서브시스템은 Rho‑GTPase 신호전달 경로를 포함하는 화학적 네트워크로, 세포-기질 접촉에서 발생하는 기계적 자극이 Rho 활성화를 촉진하고, 이는 다시 미오신 라이트 체인 인산화와 액틴-미오신 수축을 강화하는 피드백 고리를 형성한다. 모델은 미분 방정식 형태로 기술되며, 주요 변수는 세포 내 응력(σ), Rho 활성도(R), 그리고 기질 변형률(ε)이다.

정상 상태 해석을 위해 σ와 R을 기질 강성(k_sub)의 함수로 두고, 고정점의 존재와 안정성을 파라미터 공간에서 탐색하였다. 결과적으로 두 개의 분기점, 즉 하한 임계값(k_low)와 상한 임계값(k_high)이 도출되었다. k_sub < k_low 인 경우, 피드백 고리가 충분히 활성화되지 않아 σ가 거의 0에 수렴하고, 세포는 비수축 상태를 유지한다. 반대로 k_sub > k_high 에서는 Rho 신호가 강하게 유지되어 σ가 높은 고정점으로 수렴하고, 세포는 지속적인 강한 수축 상태에 머문다. 이 구간 사이에서는 두 개의 안정 고정점이 공존하며, 초기 조건에 따라 저수축 또는 고수축 상태 중 하나를 선택한다.

시간에 따라 변하는 기질 강성(k_sub(t))을 도입하면, 시스템은 비정상적인 동적 궤적을 보이며 히스테리시스 루프를 형성한다. 특히, k_sub를 서서히 증가시킬 때는 고수축 고정점으로 전이되는 임계값이 k_high보다 낮게 나타나고, 반대로 서서히 감소시킬 때는 저수축 고정점으로 복귀하는 임계값이 k_low보다 높게 나타난다. 이러한 차이는 시스템이 경로 의존적인 전이 지연(lag)과 피드백 지연을 겪기 때문이며, 모델은 이 현상을 강성 변화 속도(dk/dt)에 대한 함수로 정량화한다. 따라서 빠른 강성 변동에서는 히스테리시스가 억제되고, 느린 변동에서는 명확한 히스테리시스 루프가 관찰된다.

이와 같은 결과는 실험적으로 보고된 세포의 기계적 기억 현상과 일치한다. 특히, 줄기세포나 섬유아세포가 경도 변화에 따라 분화 혹은 이동성을 전환하는 과정에서 관찰되는 이중안정성 및 히스테리시스는, 본 모델이 제시하는 Rho‑기계 피드백 고리와 기계‑화학 결합 메커니즘으로 설명될 수 있다.