미세소관 위 키네신 집단 운동의 확률적 2차원 모델링

초록

본 연구는 키네신 모터 단백질이 미세소관을 따라 이동하는 과정을 2차원 확률적 분자동역학 시뮬레이션으로 구현한다. 손‑오버‑핸드 메커니즘을 재현하고, 모터의 평균 속도와 전류를 계산한다. 모터 밀도가 일정 임계값을 초과하면 교통 정체가 발생해 속도와 전류가 급격히 감소함을 보여준다.

상세 요약

이 논문은 세포 내 물질 수송을 담당하는 키네신의 집단적 움직임을 물리학적 모델링으로 정량화하려는 시도이다. 저자들은 2차원 격자 위에 미세소관을 선형 트랙으로 두고, 각 키네신을 두 개의 ‘발’(head)으로 구성된 입자로 표현한다. 발은 ATP 가수분해에 의해 발생하는 스프링‑같은 결합력과, 미세소관의 주기적 포텐셜(β‑튜불린 서열에 대응) 사이에서 확률적 전이(전이율은 Michaelis‑Menten 형태와 비슷한 함수)로 움직인다. 핵심은 ‘hand‑over‑hand’ 메커니즘을 구현하기 위해 두 발이 교대로 결합·해리하면서 앞발이 뒤발보다 약 8 nm 앞서 이동하도록 설계한 점이다.

시뮬레이션은 Langevin 방정식 기반의 stochastic differential equation을 시간 이산화하여 수행한다. 열잡음과 점성 저항을 포함해 실제 세포질 환경을 근사했으며, ATP 농도와 외부 부하(force)를 파라미터로 조절한다. 결과로 얻은 평균 속도 ⟨v⟩와 전류 J=ρ⟨v⟩(ρ는 모터 밀도)는 밀도 의존성을 보인다. 저밀도 구간에서는 ⟨v⟩가 거의 일정하고 J는 선형적으로 증가하지만, 임계 밀도(≈0.3 motors/µm) 이상에서는 뒤따르는 모터들이 앞 모터에 의해 물리적으로 방해받아 ‘교통 정체’가 형성된다. 이때 ⟨v⟩와 J는 급격히 감소하고, 시뮬레이션 영상에서는 고정된 군집이 장시간 지속되는 현상이 관찰된다.

이 모델은 기존 1차원 TASEP(총체적 비균등 입자 흐름) 모델과 달리 실제 2차원 공간에서의 회전·전위 변동을 포함함으로써, 키네신이 미세소관을 따라 비선형 경로를 탐색하거나, 다른 종류의 마이크로튜불린 결합 부위와 상호작용할 가능성을 탐색할 수 있다. 또한, 파라미터 스위핑을 통해 ATP 농도 감소, 외부 부하 증가, 혹은 변이된 결합 강도 등이 교통 정체에 미치는 영향을 정량적으로 예측한다. 이러한 결과는 세포 내 물질 수송 효율 저하와 관련된 신경퇴행성 질환 혹은 암세포의 과다 수송 현상을 이해하는 데 이론적 기반을 제공한다.