세포질 간헐적 운송 정량화

초록

이 논문은 바이러스와 같은 입자가 세포질에서 자유 확산과 미세소관을 따라 이동하는 유도 운동을 번갈아 가는 간헐적 동역학을 보이는 현상을 수학적으로 모델링한다. 평균 결합 시간과 유도 속도를 구해 방사형 및 원통형 세포 형태에 대한 효과적인 드리프트를 도출하고, 이를 이용해 작은 핵공을 찾는 확률과 평균 도달 시간을 추정한다.

상세 분석

본 연구는 세포질 내 입자 이동을 두 가지 기본 모드, 즉 순수 브라운 운동과 미세소관(MT) 상의 유도 운동으로 구분하고, 이들 사이의 전환을 확률적 스위칭 과정으로 모델링한다. 입자의 위치 x(t)는 확률 미분 방정식 d x = √(2D) dw (자유 확산) 혹은 d x = V dt (MT 결합) 로 기술되며, 여기서 D는 확산계수, V는 유도 속도이다. 핵심은 입자가 MT에 처음 결합하기까지의 평균 첫 통과 시간(Mean First Passage Time, MFPT)을 구하는 것으로, 이를 통해 ‘전방 결합률(forward binding rate)’을 해석적 형태로 얻는다.

방사형(2‑차원)와 원통형(3‑차원) 두 가지 기하학적 경우를 각각 분석한다. 방사형 경우, 세포를 외부 반경 R과 핵 반경 δ를 갖는 원환형 영역으로 가정하고, N개의 MT가 방사형으로 균등 배치된다. 작은 각도 Θ = 2π/N을 이용해 입자가 한 MT에 도달하기까지의 평균 확산 시간 τ(r₀)≈r₀² Θ²/(12D)와 평균 결합 반경 (\bar r(r₀)≈r₀(1+Θ²/12))를 도출한다. 이후, 유도 구간 동안 평균 이동 거리 d_m = V t_m (t_m은 결합 지속 시간)와 결합 후 방출 반경 r_f = (\bar r - d_m)를 정의하고, 전체 평균 소요 시간 τ(r₀)+t_m을 등가화함으로써 효과적인 드리프트 b(r₀)= (d_m − r₀ Θ²/12)/(t_m + r₀² Θ²/(12D)) 를 얻는다. 이 드리프트는 확산-유도 혼합 과정을 단일 확산 방정식 d x = b(x)dt+√(2D)dw 로 치환한다.

원통형 경우, MT를 원통 축에 평행한 얇은 원통으로 모델링하고, 횡단면에서의 첫 결합 시간 τ≈R² ln(1/ε)/(2ND) (ε는 MT 반경 대비 세포 반경 비) 를 구한다. 여기서도 동일하게 b = d_m/t_m + τ 로 정의하여, MT 수 N이 증가할수록 유도 속도가 선형적으로 증가함을 확인한다.

이러한 유도 드리프트를 이용해 Fokker‑Planck 방정식에 소멸 항 k(x) (예: 프로테아제에 의한 분해) 을 포함시켜, 작은 흡수 구멍(핵공)까지 도달할 확률 P_N과 평균 도달 시간 τ_N을 작은 구멍 이론(small hole theory)과 라플라스 근사법을 통해 구한다. 핵공 면적 비 ε≪1인 경우, P_N≈