생화학·기계 결합으로 본 세포 부착과 비균일 스트레스 섬유 수축 모델

초록

본 논문은 세포와 기질 사이의 부착 부위에서 발생하는 기계적 힘이 Rho 경로를 통해 생화학적 신호를 유도하고, 이 신호가 스트레스 섬유 내 마이오신 II 활성을 증가시켜 수축력을 강화한다는 양방향 피드백 메커니즘을 모델링한다. 저자들은 Rho‑경로 활성화를 반응‑확산 방정식으로 기술하고, 이를 비선형 점탄성 모델에 결합해 스트레스 섬유의 공간적 수축 변화를 예측한다. 모델은 실험에서 관찰된 활성화 및 변형의 강한 공간 구배를 재현하며, 세포 접착 성숙 과정의 물리·생화학적 연계를 정량적으로 설명한다.

상세 분석

이 연구는 세포 부착 부위에서 발생하는 기계‑생화학적 상호작용을 정량적으로 해석하기 위해 두 개의 주요 서브시스템을 결합한다. 첫 번째는 Rho‑GTPase‑ROCK‑MLC 경로를 기반으로 한 myosin II 활성화 메커니즘이며, 이는 반응‑확산 방정식으로 기술된다. 여기서 RhoA·GTP의 활성화는 외부 힘에 의해 촉진되고, 활성화된 RhoA는 ROCK을 활성화해 MLC(Myosin Light Chain)의 인산화를 촉진한다. 인산화된 MLC는 myosin II의 ATPase 활성을 증가시켜 섬유 내 수축력을 높인다. 반응‑확산 방정식은 활성화 물질의 생성·소멸 속도와 확산 계수를 포함해 공간적 비균일성을 자연스럽게 도입한다.

두 번째 서브시스템은 스트레스 섬유 자체를 비선형 점탄성(Viscoelastic) 모델로 묘사한다. 저자들은 섬유를 연속적인 1차원 체인으로 가정하고, 각 구간에 대해 탄성계수(E)와 점성계수(η)를 부여한다. myosin II에 의해 발생하는 등축성 수축력은 내부 응력 텐서에 추가되며, 이는 변형률-응력 관계식 σ = E·ε + η·∂ε/∂t + σ_contract 로 표현된다. 여기서 σ_contract는 myosin 활성도에 비례하는 수축 응력이다.

핵심은 두 서브시스템이 피드백 루프를 형성한다는 점이다. 외부 기계적 변형이 RhoA 활성화를 촉진하면, 반응‑확산 방정식에 의해 myosin 활성도가 증가하고, 이는 섬유의 수축력을 강화한다. 강화된 수축력은 다시 부착 부위에 더 큰 기계적 응력을 가해 RhoA를 추가로 활성화한다. 이러한 양성 피드백은 수치 해석을 통해 안정적인 고정점이 아닌, 공간적으로 비균일한 고활성 영역과 저활성 영역을 동시에 생성함을 보여준다.

모델 파라미터는 실험적 측정값(예: ROCK 활성화 속도, myosin ATPase 효율, 섬유 탄성계수 등)과 문헌값을 기반으로 설정했으며, 민감도 분석을 통해 확산 계수와 반응 속도 상수가 공간 구배 형성에 미치는 영향을 정량화했다. 특히 확산 계수가 작을수록 활성화 신호가 국소화되어 부착점 근처에 높은 수축력이 집중되고, 섬유 중앙부는 상대적으로 낮은 수축을 보인다. 이는 실제 현미경 영상에서 관찰되는 스트레스 섬유의 비대칭 수축 패턴과 일치한다.

이 모델의 강점은 복잡한 세포 내 신호망을 간소화하면서도, 기계적 변형과 화학적 신호가 동시에 고려된 통합 프레임워크를 제공한다는 점이다. 그러나 1차원 연속체 가정, 일정한 물성값 가정, 그리고 Rho 경로 외 다른 신호(예: Rac, Cdc42)의 무시 등은 모델의 적용 범위를 제한한다. 향후 2‑3차원 네트워크와 다중 신호 경로를 포함한 확장이 필요하다.