양극성 미세소관 다발에서 모터‑매개 양방향 수송 모델링

초록

본 연구는 곰팡이 균사체의 중앙부에 존재하는 양극성 미세소관 다발을 배경으로, 키네신‑3와 다이네인의 협동에 의해 일어나는 조기 엔도좀의 양방향 이동을 비균등 입자 흐름(ASEP) 모델로 확장하였다. 입자는 각 미세소관에서 전·후방 전환률 Ω로 방향을 바꾸고, 마이너스 말단에서 장애물‑유도 전이율 q₁ 또는 말단‑유도 전이율 q₂로 다른 미세소관으로 이동한다. 수치 시뮬레이션과 평균장 해석을 통해 전체 밀도 Θ가 낮아도 플러스·마이너스 말단 근처에 충격(밀도 급변) 현상이 발생함을 확인했으며, q₁·q₂ 값에 따라 어느 모터가 주도적인지 전환되는 메커니즘을 제시한다.

상세 분석

이 논문은 기존에 단극성 미세소관(세포 말단 근처)에서의 양방향 운반을 설명하기 위해 사용된 ASEP(비대칭 단순 배제 과정) 모델을, 중앙부의 antipolar microtubule bundle(양극성 다발)로 일반화한다는 점에서 이론적 진보를 이룬다. 모델은 두 개의 반대 방향을 가진 미세소관을 평행하게 배치하고, 입자(엔도좀‑cargo)가 각 미세소관을 따라 일방향으로 이동한다. 이동 속도는 전방(kinesin‑3)과 후방(dynein) 모터에 의해 결정되며, 입자는 일정 확률 Ω로 같은 미세소관 내에서 전·후방을 전환한다. 이는 실제 세포 내에서 모터가 결합·해리하거나, ATP 농도 변화 등에 의해 전진·후진 전환이 일어나는 현상을 추상화한 것이다.

특히, 두 미세소관은 마이너스 말단에서 물리적으로 연결되어 있어, 입자는 그 지점에서 두 가지 전이 메커니즘을 통해 다른 미세소관으로 이동할 수 있다. 첫 번째는 장애물‑유도 전이율 q₁으로, 입자가 전방 이동 중에 다른 입자와 충돌하거나 ‘병목’에 직면했을 때 옆 미세소관으로 회피하는 과정을 모델링한다. 두 번째는 말단‑유도 전이율 q₂로, 입자가 마이너스 말단에 도달했을 때 자연스럽게 다른 미세소관으로 전이되는 현상을 반영한다. 이 두 전이율은 실제 세포에서 마이너스 말단에 존재하는 교차 결합 단백질(예: dynein‑anchor)이나, 미세소관 간의 물리적 간격에 따라 달라질 것으로 추정된다.

수치 시뮬레이션에서는 전체 입자 밀도 Θ를 0.1부터 0.9까지 변화시키며, 다양한 (Ω, q₁, q₂) 조합을 탐색한다. 결과는 평균장 해석과 일치하는데, 특히 Θ가 낮아도 플러스 말단(plus‑end)과 마이너스 말단 근처에 급격한 밀도 상승(‘shock’)이 나타난다. 이는 입자들이 말단에서 대기열을 형성하고, 전이율이 제한적일 경우 병목 현상이 확대되기 때문이다. Ω가 클수록 입자들이 방향을 자주 바꾸어 충돌 가능성이 감소해 shock의 폭이 얇아지지만, q₁·q₂가 작을 경우 여전히 강한 밀도 구배가 유지된다.

또한, q₁과 q₂의 비율에 따라 어느 모터가 ‘주도’하는지가 전환된다. q₂가 q₁보다 크게 설정되면 말단‑유도 전이가 우세해, 입자들이 마이너스 말단에서 다른 미세소관으로 쉽게 이동함으로써 dynein‑주도 후방 이동이 강화된다. 반대로 q₁이 우세하면 장애물 회피 전이가 빈번해 kinesin‑3에 의한 전방 이동이 주도적이다. 이러한 전이율 조절 메커니즘은 세포가 환경 변화(예: 영양 제한, 스트레스)에 따라 모터 사용 비율을 동적으로 조절할 수 있음을 시사한다.

마지막으로, 저자들은 모델이 실제 곰팡이 세포에서 관찰된 EE(early endosome) 이동 궤적과 정량적 속도 분포를 재현함을 보이며, 실험적 파라미터(예: 모터 발현량, 미세소관 길이)와 모델 파라미터 간의 매핑 가능성을 제시한다. 이는 향후 실시간 형광 현미경 데이터와 결합해, 모터‑미세소관 상호작용을 정밀하게 추정할 수 있는 기반을 제공한다.