뇌종양 침윤 예측을 위한 확률적 다중스케일 모델

초록

본 논문은 개별 교세포의 미시적 행동을 기반으로, 뇌 조직의 이질성과 방사선 치료의 불확실성을 확률 과정으로 모델링한 다중스케일 확률적 반응‑확산‑택시 방정식 체계를 제시한다. 미시·중간·거시 단계의 연계 유도와 수치 실험을 통해 다양한 방사선 프로토콜의 종양 제어 확률(TCP), 정상 조직 합병증 확률(NTCP), 복합 성공 확률(UTCP)을 평가한다.

상세 분석

이 연구는 기존의 결정론적 종양 성장 모델이 갖는 한계를 극복하기 위해, 세포 수준의 수용체 결합·해리 동역학을 명시적으로 기술하고 이를 확률적 미분 방정식으로 전이시켰다. 특히, 세포가 조직 섬유와 상호작용하면서 발생하는 결합 변수 y₁, y₂를 질량 작용 법칙에 따라 서술하고, 방사선 피폭에 따른 살아남는 비율 S_Q, S_M을 선형‑이차(LQ) 모델 형태로 도입함으로써 치료 효과의 불확실성을 자연스럽게 포함시켰다.

중간 스케일에서는 속도‑점프 과정을 갖는 kinetic transport equation을 사용했으며, 전이율 λ(z)=λ₀−λ₁z가 세포의 내부 상태 z(=y*−y)와 연계되어 있다는 점이 핵심이다. 이는 세포가 과밀 지역을 회피하고 조직 섬유의 방향성에 따라 이동하는 ‘접촉 유도’ 현상을 수학적으로 구현한다. 전이 커널 K(x,v)는 DTI 기반의 섬유 방향 분포 q(x,θ)를 통해 환자 개별의 구조 정보를 반영한다는 점에서 임상 적용 가능성을 높인다.

또한, 증식 항 Pp를 수용체 매개 상호작용에 기반한 비선형 연산자로 정의하고, 증식률 μ(M,Q)가 종양 및 정상 조직 밀도에 의존하도록 설정함으로써 종양-미세환경의 피드백 루프를 포착한다. 방사선 효과는 Ornstein‑Uhlenbeck 형태의 확률 과정 ξ_t를 통해 시간에 따라 변동하는 선량 불확실성을 모델링한다. 이 stochastic term은 거시 방정식의 반응 항에 직접 삽입되어, 치료 결과의 확률 분포를 전파한다.

수학적 유도 단계에서는 확산‑수소동역학적 스케일링을 적용해 kinetic equation을 확률적 반응‑확산‑택시 PDE로 축소한다. 결과 방정식은 무작위 계수 D(x,ω), χ(x,ω) 등을 포함한 랜덤 PDE 형태이며, 이는 Monte‑Carlo 혹은 스펙트럼 방법으로 수치 해석이 가능하도록 설계되었다.

수치 실험에서는 하이퍼·하이퍼프랙션, 저용량 장기 방사선 등 다양한 프로토콜을 시뮬레이션하고, 각 프로토콜에 대한 TCP, NTCP, UTCP를 계산한다. 결과는 치료 계획의 강건성(robustness)을 정량화하고, 불확실성이 큰 경우 치료 성공 확률이 급격히 감소함을 보여준다. 특히, 조직 이질성이 높은 영역에서의 방사선 감수성 파라미터 ρ_M, ρ_Q 변동이 종양 재발 시점에 큰 영향을 미치는 것으로 나타났다.

이 논문의 주요 기여는 (1) 미시적 수용체 동역학을 기반으로 한 확률적 종양-치료 모델 구축, (2) kinetic theory of active particles를 이용한 다중스케일 연계, (3) 방사선 치료의 불확실성을 확률 과정으로 정량화하여 임상 지표(TCP, NTCP, UTCP)와 직접 연결한 점이다. 이러한 접근은 개인 맞춤형 방사선 치료 설계와 위험 기반 치료 최적화에 새로운 수학적 도구를 제공한다.