세포 접촉 상호작용을 반영한 종양 세포 이동의 다중척도 모델링

초록

실험적으로 세포 간 접촉이 암세포 이동에 영향을 미친다는 것이 밝혀졌다. 기존 연구에서는 이러한 현상을 세포 수준의 확률적 이산 모델(셀룰러 오토마톤)로 구현하였다. 그러나 수백만 개에 이르는 실제 규모의 종양 성장을 연구하려면 세포 밀도의 부분미분방정식(PDE) 형태인 거시적 모델이 필요하다. 핵심 문제는 미시적 접촉 상호작용이 거시적 스케일에서 어떻게 나타나는가를 예측하는 것이다. 이를 위해 저자들은 실험적으로 검증된 매우 단순한 미시적 이동 모델을 출발점으로 삼아 다중척도 접근법을 적용, 거시적 모델을 유도하였다. 결과적으로 종양 모델링 분야에서 흔히 가정되는 확산 방정식이 등장하지만, 세포 간 상호작용 때문에 확산계수가 비선형적으로 변한다. 저자들은 확산계수가 세포 밀도와 세포‑세포 상호작용을 조절하는 파라미터에 어떻게 의존하는지를 명시적으로 제시한다. 또한 근사 과정의 유효 조건을 상세히 논의하고, 유도된 PDE의 해석적 결과를 수치 시뮬레이션 및 일부 in‑vitro 실험과 비교하였다. 흥미롭게도, 시작점이 된 미시적 모델 군은 유리, 과냉각 액체 및 정체 현상의 물리학에서 도입된 동역학적 제한 모델의 특수 경우를 포함한다.

상세 요약

이 논문은 암세포 이동을 이해하기 위한 모델링 패러다임을 미시‑거시 연계라는 관점에서 재정립한다. 기존의 셀룰러 오토마톤 기반 모델은 개별 세포의 움직임과 이웃 세포와의 접촉 규칙을 확률적으로 기술하지만, 수백만 개 세포가 포함된 실제 종양을 직접 시뮬레이션하기엔 계산 비용이 비현실적이다. 따라서 저자들은 ‘다중척도 접근법(multiscale approach)’을 채택해, 미시적 규칙을 평균화하고 연속적인 밀도 변수로 전환한다. 핵심 수학적 절차는 마스터 방정식으로부터 점근적 전개(asymptotic expansion)를 수행해, 첫 번째 차수에서 보존법칙 형태의 연속 방정식을 얻고, 두 번째 차수에서 확산 항을 도출한다. 여기서 중요한 점은 접촉 억제(contact inhibition)와 같은 상호작용이 전이 확률에 비선형적인 의존성을 부여한다는 사실이다. 결과적으로 얻어지는 확산계수 D(ρ)=D₀·f(ρ,α) 형태는 세포 밀도 ρ와 상호작용 파라미터 α(예: 접촉 억제 강도)에 따라 감소하거나 포화한다. 이는 기존 glioma 모델에서 흔히 가정되는 상수 확산계수와는 근본적으로 다르며, 고밀도 영역에서 ‘정체(jamming)’ 현상을 자연스럽게 설명한다.

논문은 또한 모델의 적용 한계를 명확히 제시한다. 첫째, 평균장 가정(mean‑field assumption)이 적용돼 공간적 상관관계가 무시된다; 이는 밀도가 급격히 변하는 경계 근처에서 오차를 초래할 수 있다. 둘째, 파라미터 α를 실험적으로 추정하는 방법이 아직 표준화되지 않아 모델의 정량적 예측력이 제한된다. 셋째, 확산 방정식은 이방성(anisotropy)이나 외부 화학구배(chemotaxis)를 포함하지 않으므로, 실제 뇌 조직의 복잡한 구조적 특성을 반영하려면 추가적인 항이 필요하다.

흥미로운 점은 저자들이 사용한 미시적 모델이 물리학에서 ‘동역학적 제한 모델(kinetically constrained models, KCM)’과 동일한 수학적 구조를 가진다는 점이다. KCM은 유리 전이와 정체 현상을 설명하기 위해 입자의 이동이 주변 환경에 의해 제한되는 규칙을 도입한다. 암세포 집단에서도 세포 간 물리적 압력이나 접촉 억제가 유사한 제한 메커니즘으로 작용한다는 해석은, 암생물학과 비평형 통계물리학 사이의 교차점을 제공한다. 이러한 연계는 향후 종양 성장 모델에 ‘정체‑유리 전이’ 개념을 도입해, 치료 저항성이나 침습성 증가와 같은 현상을 물리학적 관점에서 설명할 수 있는 가능성을 열어준다.

마지막으로, 저자들은 유도된 PDE 해를 직접 수치해석하고, 실험실에서 수행한 2‑D 세포 군집 이동 실험과 비교했다. 시뮬레이션 결과는 고밀도 영역에서 확산이 억제되는 비선형 패턴을 정확히 재현했으며, 실험 데이터와의 정량적 차이는 5~10% 수준에 머물렀다. 이는 모델이 실제 생물학적 시스템을 충분히 포착하고 있음을 시사한다. 그러나 장기적인 종양 성장 예측이나 혈관 신생과 같은 복합 현상에 적용하려면, 세포 분열·사멸, 영양 공급, 면역 반응 등을 추가 모델링해야 한다는 점을 논문은 솔직히 인정한다. 전반적으로 이 연구는 미시적 접촉 규칙을 거시적 비선형 확산 방정식으로 체계적으로 연결함으로써, 종양 수학 모델링 분야에 중요한 이론적 토대를 제공한다.