활성 유체 모델로 보는 종양 침윤의 손가락 불안정성

초록

본 논문은 세포가 자체적으로 힘을 발생시키는 활성 조직과, 성장하는 종양이 주변 수동 조직에 침투하면서 나타나는 손가락 모양 돌출 현상을 최소 4개의 물리적 파라미터(마찰, 활성 구동, 성장률, 표면 장력)만으로 설명하는 연속체 유체 모델을 제시한다. 2차원 원형 종양을 대상으로 선형 안정성 분석을 수행해 활동성(α)과 성장 속도(k)가 손가락 불안정성을 촉진한다는 결론을 도출하고, 전이 시작 시점의 종양 크기와 발생 가능한 손가락 수에 대한 정량적 예측식을 제시한다.

상세 분석

이 연구는 기존의 기계적 불안정성 모델이 간과해 온 ‘활성 힘’이라는 핵심 요소를 명시적으로 포함한다는 점에서 혁신적이다. 저자들은 조직을 두 종류의 유체—활성 성장 유체와 수동 비성장 유체—로 구분하고, 각각에 대해 압력·속도장을 다음과 같이 기술한다: ∇p = –β v + α v|v|, ∇·v = k (활성 조직)와 ∇p₀ = –β₀ v₀, ∇·v₀ = 0 (수동 조직). 여기서 β는 기판에 대한 마찰계수, α는 세포가 기판에 가하는 지속적인 추진력, k는 체적 성장률, γ는 인터페이스 표면 장력이다. α와 β가 같은 차원(힘·부피⁻¹)으로 정의된 점은 기존 Toner‑Tu 모델과 차별화되며, 속도 크기에 독립적인 활성 구동을 가정함으로써 분석을 단순화한다.

파라미터 추정은 실제 실험 데이터(세포 분열 시간, 세포‑기판 결합 수, 개별 세포가 발생시키는 힘 등)를 기반으로 하며, k≈10⁻⁴ s⁻¹, β≈10¹⁵–10¹⁶ Pa·s·m⁻², α≈10⁸–10¹⁰ Pa·m⁻¹, γ≈10⁻³–10⁻² Pa·m 로 설정한다. 이 값들을 이용해 특성 길이 ℓ = (2γβ/k)¹⁄³ ≈ 10 µm, 기준 활동 α* = βℓk ≈ 10⁷ Pa·m⁻¹ 등을 정의하고 무차원화한다.

선형 안정성 분석에서는 원형 인터페이스에 작은 주기적 변동 r(θ)=r₀+δrₙe^{inθ}를 가정하고, 각 모드 n에 대한 성장률 σₙ를 도출한다. σₙ는 α, φ=β₀/β, γ, k, r₀ 등에 의존하며, α=0(수동 경우)에서는 전형적인 비활성 유체의 라플라스-스털링 불안정성 식을 재현한다. 반면 α>0이면 σₙ에 추가적인 양의 항이 들어가, 특히 중간 규모(수십 마이크로미터) 종양에서 n≈5–10 정도의 모드가 가장 크게 성장한다는 결과가 나온다. 이는 실험적으로 관찰되는 ‘손가락’ 수와 일치한다. 또한, 표면 장력 γ이 클수록 불안정성 임계 반경 r_c가 증가하고, φ가 크면(즉, 주변 조직이 더 점성이면) 불안정성이 억제된다.

저자들은 이론적 결과를 이용해 전이가 시작되는 종양 반경 r_c ≈ ℓ·(α/α*)^{1/2}·(γ/γ₀)^{1/3} 형태의 스케일링 법칙을 제시하고, 손가락 수 N≈2πr_c/λ_max (λ_max는 가장 불안정한 파장의 2π/n) 로 추정한다. 파라미터 범위 내에서 예측된 r_c는 30–80 µm, N은 4–12개 정도이며, 이는 실제 종양 조직에서 보고된 값과 좋은 일치를 보인다.

이 모델의 강점은 (1) 최소 파라미터로 복잡한 현상을 포착, (2) 분석적 해를 제공해 파라미터 의존성을 명확히 파악, (3) 성장·활동·마찰·표면 장력 간의 상호작용을 정량화한다는 점이다. 한계로는 2D 평면 가정, 파라미터가 시간·공간에 균일하다는 가정, 그리고 세포 간 기계적 이질성(예: 리더 세포와 추종 세포의 차이)을 무시한다는 점이 있다. 향후 3D 확장, 비균일 α(r,t) 도입, 그리고 실험적 마이크로플루이드 플랫폼과의 정량적 검증이 필요하다.