면역 반응을 고려한 종양 성장 확률 모델

초록

본 논문은 로지스틱 성장식에 면역 억제 항목을 도입하고, 내부 성장 과정에서 발생하는 곱셈성 잡음과 외부 환경의 영향을 나타내는 가법성 잡음을 동시에 고려한 확률적 종양 성장 모델을 제시한다. 두 잡음은 유한한 상관시간을 갖는 색 잡음이며, 상호 교차 상관도 포함한다. 정적 확률분포 Ps를 유도하고, 면역 활성화 정도와 교차 상관 강도에 따른 Ps 의 피크 변화를 분석한다. 또한 평균 초과 시간(MFPT)을 계산해 면역 및 교차 상관이 종양 소멸에 미치는 영향을 정량화한다. 결과는 세 단계(정상, 성장, 소멸) 모델과 일치한다.

상세 분석

본 연구는 기존의 로지스틱 종양 성장 모델에 두 가지 중요한 확장성을 부여한다. 첫째, 사망률 항에 면역화(immunization)라는 결정적 항을 추가함으로써, 면역계가 종양 세포를 제거하는 효과를 수학적으로 구현한다. 이 항은 면역 활성화 비율 γ 으로 파라미터화되며, γ가 클수록 사망률이 증가해 종양 성장 억제가 강화된다. 둘째, 성장률(출산률)에는 내부 생물학적 변동을 반영한 곱셈성 색 잡음 ξ₁(t) 를, 외부 환경 변동을 반영한 가법성 색 잡음 ξ₂(t) 를 동시에 도입한다. 두 잡음은 각각 평균 0, 자기상관함수 ⟨ξ_i(t)ξ_i(t′)⟩= (D_i/τ_i) exp(−|t−t′|/τ_i) 를 갖고, 상관시간 τ_i 와 강도 D_i 가 정의된다. 특히, ξ₁과 ξ₂ 사이에 비대칭 교차 상관 ⟨ξ₁(t)ξ₂(t′)⟩= (λ/τ_c) exp(−|t−t′|/τ_c) 를 포함시켜, 내부 성장 과정과 외부 환경이 동시에 영향을 주는 현실적인 상황을 모델링한다.
이러한 설정을 Langevin 방정식 형태로 기술하고, 대응하는 Fokker‑Planck 방정식을 유도한다. 색 잡음의 유한 상관시간을 처리하기 위해 차원 축소와 적절한 마코프 근사를 적용하여, 효과적인 확산계수와 드리프트 항을 얻는다. 정적 해 Ps(x) 는 일반적인 베타‑가우시안 형태를 띠지만, 면역 파라미터 γ 와 교차 상관 강도 λ 에 따라 피크 위치와 폭이 크게 변한다. γ가 증가하면 Ps 의 피크가 낮은 종양 세포 수 영역으로 이동하면서, 종양 소멸 가능성이 높아진다. 반면 λ가 양의 값으로 커지면, 내부 성장 변동과 외부 환경 변동이 동조화되어 피크가 강화되고, 종양이 일정 수준 이상 성장할 확률이 증가한다.
MFPT 계산은 반사 경계(정상 상태)와 흡수 경계(소멸 상태)를 설정한 첫 통과 시간 문제로 전개된다. Kramers‑type 근사를 이용해 잠재적 장벽 높이와 폭을 추정하고, 면역 파라미터 γ와 교차 상관 λ 가 MFPT에 미치는 영향을 정량화한다. 결과적으로 γ가 클수록 장벽이 낮아져 MFPT가 급격히 감소하고, 이는 면역 치료가 종양을 빠르게 소멸시킬 수 있음을 시사한다. 반대로 λ가 양이면 장벽이 높아져 MFPT가 늘어나며, 이는 환경 요인과 내부 성장 변동이 동시에 강화될 때 종양이 더 오래 지속될 가능성을 의미한다.