통제와 무통제 확률적 노턴 사이먼 마사게 종양 성장 모델

초록

본 연구는 기존의 생물학적 근거를 가진 노턴‑사이먼‑마사게(NSM) 종양 성장 모델에 확률적 요소를 도입하여 무통제와 화학요법 통제 두 가지 상황을 분석한다. 무통제 모델에서는 종양 부피의 양성을 보장하는 파라미터 조건을 도출하고, 최대우도 추정과 혼합 효과 모델을 이용해 마우스 실험 데이터를 적합한다. 통제 모델에서는 약물 투여 효과를 포함시켜 개방형 및 폐쇄형 시뮬레이션을 통해 약물 파라미터와 확률적 변동성에 대한 민감도를 평가한다.

상세 분석

논문은 먼저 기존 NSM 모델의 결정론적 형태를 확률 미분 방정식으로 일반화한다. 종양 부피 V(t)를 상태 변수로 두고 성장률을 αV^β와 사망률을 γV 형태로 표현한다. 확률적 항은 Wiener 과정 W(t)와 스케일링 파라미터 σ를 곱한 형태로 추가되어, 환경 변동성이나 측정 오차를 내재화한다. 무통제 경우에는 dV = (αV^β – γV)dt + σV dW 로 기술되며, 여기서 V>0를 유지하기 위한 α,β,γ,σ의 관계를 이론적으로 증명한다. 특히, β가 1보다 작을 때 성장 항이 포화 효과를 나타내어 부피가 무한히 발산하지 않으며, σ가 너무 크면 음수 값이 발생할 위험이 있기에 σ^2 < 2γ 조건을 제시한다. 파라미터 추정은 개별 마우스 데이터를 기반으로 한 최대우도법을 사용하고, 인구 수준의 변동성을 반영하기 위해 비선형 혼합 효과 모델(NLME)을 구축한다. 이때 고정 효과는 α,β,γ와 같은 평균값을, 랜덤 효과는 개체별 편차와 σ를 포함한다. 추정 결과는 기존 문헌값과 일치하면서도 확률적 요소가 모델 적합도를 현저히 향상시킴을 보여준다. 통제 모델에서는 약물 농도 C(t)를 추가 변수로 두고, 약물 효과를 –kC(t)V 형태의 억제 항으로 도입한다. 약물 동역학은 일차적 일방향 흡수와 일차적 제거를 가정한 일차 미분식으로 기술되며, 제어 입력 u(t)로 약물 투여 스케줄을 정의한다. 폐쇄형 피드백 제어는 V(t)와 C(t) 상태에 기반한 선형 상태 피드백 u(t)=K1V+K2C 로 설계하고, 시뮬레이션을 통해 안정성 영역을 탐색한다. 결과는 약물 투여 강도가 σ에 비례해 민감하게 변함을 보여주며, 높은 σ(큰 환경 변동성)에서는 보다 보수적인 투여 전략이 필요함을 시사한다. 또한, 개방형 시나리오에서는 일정 간격 고정 투여가 평균 종양 부피 감소에 효과적이지만, 변동성이 큰 경우에는 과다 투여로 인한 부작용 위험이 증가한다는 점을 강조한다. 전반적으로 논문은 확률적 NSM 모델이 종양 성장과 치료 반응을 보다 현실적으로 묘사할 수 있음을 입증하고, 파라미터 조건과 제어 설계에 대한 이론적 가이드를 제공한다.