키프3가 이끄는 미세소관 탈중합 메커니즘 수학적 모델

초록

이 논문은 키네신‑8 패밀리 모터인 Kip3p가 미세소관(MT) 끝에서 진행 방향으로 이동하면서 동시에 탈중합을 촉진하는 과정을 단순화된 수학 모델로 구현한다. 모델은 Kip3p의 이동성, MT 말단에서의 불안정화, 그리고 탈중합 과정의 프로세시티를 포함한다. 실험 데이터와 비교했을 때, 모델은 MT 길이에 따른 탈중합 속도와 길이‑시간 곡선의 변동성을 정확히 재현한다. 또한, 짧은 MT에서만 길이 의존적 탈중합이 나타나며, 이는 전체 모터 농도에 민감하게 반응한다는 점을 제시한다.

상세 분석

본 연구는 Kip3p가 미세소관의 플러스 말단에서 수행하는 두 가지 핵심 기능, 즉 양방향 이동과 탈중합 촉진을 하나의 연속적인 확률 과정으로 통합한 수학적 프레임워크를 제시한다. 모델은 먼저 Kip3p가 자유 풀에서 결합하고, 1차원 격자 위를 일정한 속도로 전진하는 과정을 마이클리시안 이동 방정식으로 기술한다. 이때 모터의 결합·해리율(k_on, k_off)과 이동 속도(v)라는 파라미터가 실험적으로 측정된 값과 일치하도록 설정된다.

플러스 말단에 도달한 Kip3p는 두 가지 가능한 행동을 갖는다. 첫 번째는 단순히 말단에서 해리되는 것이고, 두 번째는 말단의 튜불린 이합체를 하나씩 제거하는 탈중합 작용이다. 탈중합은 확률적 사건으로 모델링되며, 탈중합 속도(k_dep)와 프로세시티(p)라는 두 파라미터로 정의된다. 프로세시티 p는 한 모터가 연속적으로 몇 개의 튜불린을 제거할 수 있는지를 나타내며, p = 1이면 완전 프로세시티, p < 1이면 비프로세시티를 의미한다.

모델은 마스터 방정식을 통해 전체 MT 길이 L(t)의 확률 분포를 시간에 따라 전개한다. 이때 중요한 결과는 평균 길이 ⟨L(t)⟩와 그 변동성(분산) σ²(L)이다. 실험 데이터와 비교했을 때, 모델은 다음과 같은 핵심 현상을 재현한다. 첫째, MT 길이가 짧아질수록 탈중합 속도가 증가하는 ‘길이 의존적 탈중합’ 현상이 나타난다. 이는 말단에 도달하는 모터의 평균 수가 L⁻¹에 비례하기 때문에, 짧은 MT에서는 말단에 모터가 더 자주 도달해 탈중합이 가속화되는 효과이다. 둘째, 탈중합이 프로세시티를 갖는 경우, 길이 변동성이 크게 증가한다. 이는 한 모터가 여러 튜불린을 연속적으로 제거함으로써 순간적인 길이 감소가 급격히 일어나기 때문이다. 셋째, 전체 모터 농도 c가 증가하면 길이 의존적 탈중합의 전환점(L_c)이 더 긴 MT에서도 나타나며, 이는 말단에 도달하는 모터 수가 농도에 비례해 증가하기 때문이다.

수치 시뮬레이션을 통해 파라미터 공간을 탐색한 결과, 실험적으로 관찰된 탈중합 속도와 변동성을 동시에 만족시키는 파라미터 집합은 k_dep ≈ 0.5 s⁻¹, p ≈ 0.8 정도이며, 이는 Kip3p가 거의 프로세시티를 가지고 있지만 완전 프로세시티는 아니라는 결론을 뒷받침한다. 또한, 모델은 ‘포화 탈중합’ 현상을 예측한다. 즉, 모터 농도가 충분히 높을 경우, 말단에 모터가 포화 상태에 이르러 탈중합 속도가 더 이상 MT 길이에 의존하지 않게 된다. 이와 같은 예측은 향후 실험적 검증을 통해 Kip3p의 동역학적 특성을 보다 정밀하게 규명하는 데 활용될 수 있다.