종양 면역 시스템의 확률 모델 수학적 분석

초록

본 논문은 종양-면역 상호작용을 기술하기 위해 Wiener 과정 기반의 확률 모델을 도입하고, Lyapunov 지수를 활용한 평형점의 확률적 안정성을 분석한다. Kuznetsov‑Taylor형 모델과 Bell형 모델 두 가지를 제시하고, 2차 Euler 스키마를 이용해 수치 시뮬레이션을 수행하였다.

상세 분석

본 연구는 기존의 결정론적 종양‑면역 모델이 실제 생물학적 시스템에서 나타나는 불확실성을 충분히 반영하지 못한다는 점에 착안하여, 확률 미분 방정식(SDE) 형태로 모델을 재구성하였다. 특히, 외부 환경 변동이나 세포 간 상호작용의 미세한 변동을 Wiener 과정(표준 브라운 운동)으로 모델링함으로써, 잡음이 시스템의 안정성에 미치는 영향을 정량화하였다. 모델의 핵심은 두 개의 상태 변수, 즉 종양 세포 집단 (x(t))와 면역 세포 집단 (y(t))를 정의하고, 각각의 성장·소멸률에 비선형 상호작용 항과 확률적 교란 항을 추가한 형태이다.

안정성 분석에서는 평형점 ((x^,y^)) 주변의 선형화된 SDE를 고려하고, Ito 미분법을 적용하여 Lyapunov 지수를 도출하였다. Lyapunov 지수가 음수이면 해당 평형점이 확률적으로 안정함을 의미한다. 논문에서는 파라미터 공간을 체계적으로 탐색하여, 잡음 강도 (\sigma)와 면역 반응 파라미터 (c)가 Lyapunov 지수에 미치는 영향을 그래프와 수치표로 제시하였다. 흥미롭게도, 일정 수준 이상의 잡음은 종양 성장 억제에 기여하여, 결정론적 모델에서는 불가능한 “잡음에 의한 안정화” 현상을 확인하였다.

두 가지 구체적 모델인 Kuznetsov‑Taylor형과 Bell형은 각각 다른 면역 반응 메커니즘을 반영한다. Kuznetsov‑Taylor형은 면역 세포의 활성화가 종양 부피에 비선형적으로 의존하는 구조를 가지고 있으며, Bell형은 면역 세포가 종양 항원을 인식하는 확률적 전이 과정을 포함한다. 두 모델 모두 동일한 확률적 프레임워크 내에서 분석되었으며, Lyapunov 지수와 수치 시뮬레이션 결과가 일관되게 잡음이 시스템을 안정화시키는 방향으로 작용함을 보여준다.

수치 해법으로는 2차 Euler‑Maruyama 스키마를 채택하였다. 이는 1차 스키마에 비해 시간 전진 오차가 (O(\Delta t^{2}))로 감소하여, 잡음이 강하게 작용하는 상황에서도 정확한 궤적을 재현한다. 시뮬레이션은 다양한 초기 조건과 파라미터 조합에 대해 수행되었으며, 종양 세포 수가 일정 수준 이하로 감소하거나, 면역 세포가 지속적으로 활성화되는 두 종류의 안정된 궤적을 관찰하였다. 이러한 결과는 이론적 Lyapunov 분석과 수치적 증거가 서로 보강함을 시사한다.

본 논문의 주요 기여는 (1) 종양‑면역 상호작용에 확률적 요소를 체계적으로 도입한 수학적 모델을 제시하고, (2) Lyapunov 지수를 이용한 확률적 안정성 기준을 구축했으며, (3) 2차 Euler 스키마를 통한 고정밀 수치 시뮬레이션을 제공했다는 점이다. 특히, 잡음이 치료 전략 설계에 긍정적 역할을 할 수 있다는 통찰은 향후 면역 치료법 개발에 새로운 방향성을 제시한다.