종양‑면역 상호작용의 공간적 동역학을 위한 동역학적 유도 모델

초록

본 논문은 종양 세포와 면역 세포 간의 상호작용을 미시적 충돌·전이 과정을 기반으로 기술한 두 종류의 동역학 방정식을 제시한다. 첫 번째는 면역 세포가 증식·사멸 없이 활성·비활성 상태를 전환하는 보존 모델이며, 두 번째는 면역 세포의 증식과 사멸을 포함한 비보존 모델이다. 각각의 미시 모델에서 적절한 스케일링을 적용해 선형 확산, 비선형 교차확산, 비선형 자체확산을 포함하는 거시 PDE 시스템을 유도한다. 거시 시스템의 균형점과 안정성을 동역학 시스템 이론으로 분석하고, 수치 실험을 통해 이론적 예측을 검증한다. 연구 결과는 공간적 상호작용이 종양‑면역 역학에 미치는 영향을 정량적으로 설명하며, 향후 치료 전략 설계에 활용될 수 있다.

상세 분석

이 연구는 세포 수준의 미시적 상호작용을 볼츠만형 적분 연산자로 모델링하고, 이를 통해 종양‑면역 시스템의 전이 방정식을 구축한다. 핵심 변수는 위치 x, 미세속도 v, 그리고 활성도 u(−1~1)이며, 종양 세포(i=1), 활성 면역 세포(i=2), 비활성 면역 세포(i=3) 각각에 대해 분포함수 f_i(t,x,v,u)를 정의한다. 충돌 연산자 Q_{ij}와 J_{ijk}는 각각 두 세포 간의 파괴·증식 상호작용과 세 종류 이상의 세포가 동시에 관여하는 복합 상호작용을 기술한다. 여기서 충돌 빈도 d_{ij}, μ_{ij}, α_{ijk}와 생성 밀도 ε_{ij}, ν_{ijk}는 생물학적 실험 데이터를 통해 추정 가능한 파라미터이다.

두 가지 시나리오를 구분한다. 보존 시나리오에서는 면역 세포가 활성↔비활성 전환만을 겪으며, μ_{ij}=0, 즉 증식·사멸 항이 없고, 전이 확률 p(u→u′)만 존재한다. 이 경우, 속도 전이 연산자 L는 런‑앤‑텀블 모델을 채택해 평균 자유 경로와 방향 전환 빈도 λ(x,v)를 도입한다. 스케일링 과정에서 작은 평균 자유 경로 ε→0와 시간 스케일 τ∼ε^2를 가정하면, 마크로스코픽 한계에서 종양·면역 세포 농도 n_i(t,x)는 선형 확산 항 D_iΔn_i와 함께 면역 세포에 대한 교차확산 항 ∇·(χ_{ij} n_i∇n_j) 형태를 갖는다. 교차확산 계수 χ_{ij}는 미시 충돌 커널의 비대칭성에 의해 유도되며, 면역 세포가 종양 농도 구배에 따라 이동하는 현상을 포착한다.

비보존(증식) 시나리오에서는 μ_{ij}>0 및 ε_{ij}가 존재해 면역 세포의 증식·사멸을 모델링한다. 이 경우, 마크로스코픽 한계에서 면역 세포에 대한 자체확산 항 ∇·(D_i(n_i)∇n_i)와 비선형 성장·소멸 항 r_i(n_1,n_2,n_3)도 동시에 나타난다. 특히, 자체확산 계수 D_i(n_i)=D_i^0+α_i n_i와 같은 형태가 도출되어, 세포 밀도가 높을수록 이동성이 증가하거나 감소하는 현상을 설명한다.

거시 PDE 시스템에 대해 균형점(steady states)을 구하고, 야코비 행렬을 이용해 선형 안정성을 분석한다. 종양 전용 균형점(종양만 존재)과 면역‑종양 공존 균형점이 존재함을 보이며, 파라미터(예: 면역 활성화율, 증식률, 교차확산 강도)에 따라 Hopf 및 Turing‑type bifurcation이 발생할 수 있음을 확인한다. 특히, 교차확산이 강할 경우 공간적 패턴(종양‑면역 구역화)이 나타나며, 자체확산이 비선형이면 패턴의 크기와 형태가 크게 변한다.

수치 실험에서는 1‑D 및 2‑D 도메인에 대해 유한 차분/유한 요소 방법을 적용하고, 초기 조건을 무작위 작은 교란으로 설정한다. 보존 모델에서는 면역 세포가 종양 주변으로 집결하면서 종양이 억제되는 동역학을, 증식 모델에서는 면역 세포가 과도하게 증식해 종양을 완전히 소멸시키거나, 반대로 면역 억제에 의해 종양이 폭증하는 두 가지 상이한 시나리오가 재현된다. 시뮬레이션 결과는 이론적 안정성 분석과 일치하며, 파라미터 스위칭을 통해 치료적 개입(예: 인터루킨 투여, 면역 세포 활성화) 효과를 정량적으로 평가할 수 있음을 시사한다.

전반적으로, 이 논문은 미시적 충돌·전이 메커니즘을 기반으로 한 동역학적 프레임워크를 제시함으로써, 종양‑면역 상호작용의 공간적 복잡성을 수학적으로 정형화하고, 교차확산·자체확산과 같은 비선형 확산 효과가 종양 진행 및 면역 억제에 미치는 역할을 명확히 밝힌다. 이러한 접근은 향후 맞춤형 면역 치료 전략 설계와 다중 스케일 모델링에 중요한 이론적 토대를 제공한다.