소음이 돕는 종양‑면역 세포 상호작용

초록

본 논문은 종양 세포, 효소 세포, 종양 탐지 세포로 구성된 3상 모델에 자가상관 잡음과 교차상관 잡음을 도입하여, 잡음이 종양 성장 억제와 완전 소멸을 유도할 수 있음을 이론적으로 입증한다. 자동상관 시간 τ와 교차상관 강도 R의 변화에 따라 종양 휴면, 진행, 소멸 등 다양한 동역학적 상태가 나타난다.

상세 분석

본 연구는 기존의 결정론적 종양‑면역 모델이 안정적인 종양‑무료 상태를 제공하지 못한다는 한계를 인식하고, 모든 세포 종류에 내재된 내재적·외재적 잡음이 시스템의 동역학을 근본적으로 변화시킬 수 있음을 보여준다. 모델은 종양 세포 X, 효소 세포 Y, 종양 탐지 세포 Z의 세 가지 변수로 구성되며, 결정론적 부분은 로지스틱 성장(a·X‑b·X²)과 효소‑종양 상호작용(c·X·Y), 효소 세포의 탐지 세포에 의한 증식(e·Y·Z) 및 각각의 자연 사멸(ρ·Y, μ·Z)으로 기술된다. 이때 X·Y·Z 사이의 곱항이 존재함에도 불구하고, ρ와 μ가 양수이므로 효소와 탐지 세포가 지속적으로 공급되지 않으면 종양은 최종적으로 수용량 b‑1에 수렴한다.

잡음 도입은 각 변수에 대해 다중곱 형태의 곱셈 잡음 Ωij·ηj(t)를 추가함으로써 구현된다. 특히, 종양 세포와 탐지 세포 사이에 교차상관 잡음 R이 존재하도록 설정했으며, 이는 ηx와 ηz 사이의 백색 잡음 상관을 의미한다. 자가상관 잡음은 유한한 상관시간 τ를 갖는 색 잡음으로 모델링되어, τ→0일 때 백색 잡음으로 수렴한다. 잡음 강도 행렬 D는 D_x, D_y, D_z와 교차상관 강도 S, R, P를 포함한다.

색 잡음에 대한 Fokker‑Planck 방정식은 τ가 작은 경우에 한해 1차 τ 전개를 사용해 근사적으로 도출되었다. 연산자 L은 확산 항과 drift 항을 포함하며, τ와 D가 결합된 고차 항이 존재한다. 이로부터 평균값 ⟨x_j⟩의 시간 진화식(식 11)이 얻어지며, τ와 R, S, P가 시스템의 고정점 위치와 안정성을 결정한다는 점이 핵심이다.

수치 해석에서는 τ와 R을 주요 파라미터로 스위핑하여 세 가지 동역학적 시나리오를 확인한다. (1) 작은 τ와 약한 R에서는 종양이 기존의 결정론적 고정점인 b‑1에 머무르며, 잡음이 거의 영향을 미치지 않는다. (2) τ가 중간 수준이면서 R이 증가하면 종양 세포 수가 감소하고, 효소와 탐지 세포가 지속적으로 활성화되어 종양이 ‘휴면’ 상태에 머무른다. (3) τ가 충분히 크고 R이 강하면 효소와 탐지 세포의 생산이 잡음에 의해 촉진되어 종양 세포가 완전히 소멸한다. 특히, 교차상관 R이 양의 값일 때만 이러한 소멸 현상이 나타나며, R이 음이면 오히려 종양 성장 속도가 증가한다.

이러한 결과는 면역 편집(immune editing) 과정에서 발생하는 유전적 변이와 세포 간 신호 전달의 불확실성을 잡음으로 해석할 수 있음을 시사한다. 모델 파라미터는 실험적 데이터(예: 종양 성장률 a≈0.1 day⁻¹, 사멸률 ρ≈0.05 day⁻¹ 등)를 기반으로 설정되었으며, 잡음 강도 D_x≈10⁻³, D_y≈10⁻⁴, D_z≈10⁻⁴ 수준이 현실적인 변동성을 반영한다. 따라서, 이론적 분석은 실제 종양‑면역 상호작용에서 잡음이 치료적 면역 활성화에 기여할 가능성을 제시한다.