다중 차선 ASEP 모델에서 양방향 운송과 펄싱 현상

초록

**
본 논문은 미세소관의 13개 원섬유를 차선으로 삼아, 양쪽 방향으로 움직이는 입자들의 집합적 이동을 ASEP(비대칭 단순 배제 과정) 형태로 모델링한다. 입자는 ‘당겨당김(tug‑of‑war)’에 의해 이동 방향을 바꾸고, 충돌 시 인접 차선으로 전이한다. 플러스‑엔드(plus‑end)에서 무유출 경계조건을 두고 입자 주입률을 변화시켰을 때, 팁에 축적되는 입자 수(팁 크기)가 주입률에 대해 비선형적으로 증가하고, 임계 주입률을 초과하면 시스템 전체가 주기적으로 충전·소진되는 펄싱 상태로 전이함을 보였다. 또한 실험에서 측정된 차선 전이율이 운송 효율을 최적화하는 값에 가깝다는 점을 제시한다.

**

상세 분석

**
본 연구는 미세소관(MT)의 구조적 특성을 반영해 13개의 원섬유를 ‘차선’으로 모델링하고, 각각의 차선에 입자를 배치한 다중 차선 ASEP 프레임워크를 구축하였다. 입자는 플러스‑타입(plus‑type)과 마이너스‑타입(minus‑type) 두 종류로 구분되며, 각각은 양극(plus‑end)과 음극(minus‑end) 방향으로 전진한다. 입자 이동은 전진(step forward)과 차선 전이(lane‑change) 두 가지 메커니즘으로 이루어지며, 전진은 전방에 빈 자리가 있을 때만 허용되는 배제 원리를 따른다. 차선 전이는 두 경우로 나뉜다. 첫째, 입자가 ‘막힘(blocked)’ 상태일 때는 인접 차선으로 전이할 확률 p_{l→k,±,b}가 적용되고, 둘째, ‘비막힘(unblocked)’ 상태일 때는 p_{l→k,±,u}가 적용된다. 또한, ‘당겨당김(tug‑of‑war)’ 이벤트를 통해 입자는 방향을 바꾸면서 차선을 전환할 수 있는데, 이 전이는 p_{l→k,+−}와 p_{l→k,−+}라는 두 개의 전이율로 기술된다.

경계조건은 한쪽 끝(plus‑end)에서 무유출(no‑flux) 조건을, 반대쪽 끝(minus‑end)에서는 입자 주입률 α_{l}^{+} (플러스‑타입) 혹은 α_{l}^{−} (마이너스‑타입)를 부여한다. 출구율 β_{l}^{±}는 각각의 차선에서 입자가 시스템을 떠나는 속도를 나타낸다. 논문은 특히 α_{l}^{−}=0, β_{l}^{+}=0인 경우, 즉 플러스‑타입 입자만이 minus‑end에서 주입되고 plus‑end에서 배출되지 않는 상황을 집중적으로 분석한다.

수학적 접근으로는 평균장(mean‑field) 근사를 이용해 연속적인 좌표 x=iδ (δ=1/N)로 변환하고, 입자 밀도 ρ_{l}(x)와 σ_{l}(x) (플러스·마이너스 입자 각각의 평균 점유율)를 정의한다. 이때, 차선 전이율이 대칭적이고 차선 간에 균일하면, 전체 시스템은 효과적으로 1차원 ASEP와 유사한 흐름 방정식을 만족한다. 그러나 차선 전이율이 비대칭이거나 ‘당겨당김’ 전이율이 크게 변할 경우, 입자들의 축적 현상이 급격히 강화된다.

주요 결과는 두 가지 비정상적 현상을 확인한 것이다. 첫째, 주입률 α^{+}가 낮을 때는 팁 크기 n_{tip}이 선형적으로 α^{+}에 비례한다는 기존 연구와 일치한다. 그러나 α^{+}가 일정 임계값 α_{c}에 접근하면, 팁에 축적되는 입자 수가 급격히 비선형적으로 증가한다. 이는 차선 전이와 방향 전이가 결합된 ‘트랩 메커니즘’에 의해 입자들이 플러스‑엔드 근처에 고정되기 때문이다. 수식적으로는 n_{tip}≈(1−x_{s})N≈α^{+}·p_{d}·exp(N·p_{u}/p_{−})/(N·p_{d}·p_{+})와 같은 형태로 표현되며, α^{+}→α_{c}에서 특이점이 발생한다.

둘째, α^{+}>α_{c} 구간에서는 시스템 전체가 주기적인 충전·소진(pulsing) 현상을 보인다. 이 펄싱 상태는 전체 입자 수가 시간에 따라 거의 정현파 형태로 변동하며, 팁에 축적된 입자들이 일정 시점에 대량 방출되고 다시 축적되는 과정을 반복한다. 시뮬레이션 결과는 이 현상이 차선 전이율 p_{±,u/b}와 ‘당겨당김’ 전이율 p_{+−}, p_{−+}에 민감하게 의존함을 보여준다. 특히, 실험에서 측정된 dynein의 차선 전이율(p_{lane}≈0.04 s⁻¹)과 kinesin의 차선 전이 억제(p_{lane}≈0) 사이의 차이가, 모델이 예측하는 최적 효율 영역에 거의 일치한다는 점은 생물학적 의미가 크다.

결론적으로, 다중 차선 ASEP 모델은 미세소관 기반 양방향 운송의 핵심 물리적 메커니즘을 정량적으로 포착한다. 차선 전이와 ‘당겨당김’ 전이율이 적절히 조절될 때, 입자 흐름은 효율적인 팁 축적과 최소 지연을 동시에 달성할 수 있으며, 과도한 주입은 펄싱 현상을 유발해 전체 운송 효율을 저하시킨다. 이러한 결과는 세포 내 소기관 운반 메커니즘을 이해하고, 인공 나노운송 시스템 설계에 중요한 설계 원칙을 제공한다.

**