종양코드 성장의 다상 모델링과 정성 분석

초록

본 논문은 변형 가능한 다공성 매체 이론을 기반으로 종양코드가 혈관 근처에서 성장하는 과정을 거시적으로 모델링한다. 종양을 증식 세포, 세포외액, 세포외기질의 포화 혼합물로, 주변 정상 조직을 비증식 세포의 혼합물로 설정하고, 두 영역 간 상호작용을 핵심 메커니즘으로 포함한다. 세포외기질을 강체로 가정해 모델을 단순화하면 세포 부피비와 영양소 농도의 두 편미분 방정식이 도출되며, 영양소 농도에 따라 세포 증식·사멸이 결정된다. 2차원 수치 시뮬레이션과 정성적 수학적 분석을 통해 주요 현상을 검증한다.

상세 분석

이 연구는 다상 연속체 역학을 종양학에 적용한 점이 가장 큰 혁신이다. 종양 조직을 ‘포화 혼합물(saturated mixture)’이라 가정함으로써, 각 성분(증식 세포, 세포외액, 세포외기질)의 질량·운동량 보존식을 동시에 기술한다. 특히, 혈관으로부터 공급되는 영양소(산소·포도당 등)를 ‘영양소 농도 필드’로 도입해, 영양소 확산·소모와 세포 증식·사멸을 연계시킨다. 모델링 과정에서 다음과 같은 핵심 가정을 두었다. 첫째, 세포외기질(ECM)은 변형이 없고 재구성되지 않는 강체로 취급해 기계적 응답을 무시한다. 이는 ECM이 종양 성장 초기에 크게 변하지 않는다는 실험적 관찰에 근거한다. 둘째, 정상 조직은 비증식 세포만으로 구성된 포화 혼합물로, 영양소 소모율이 종양 조직보다 낮다고 가정한다. 셋째, 혈관은 영양소의 고정된 경계조건을 제공하는 ‘영양소 공급원’으로 모델링한다. 이러한 가정 하에 전체 시스템은 두 개의 비선형 편미분 방정식으로 축소된다. 첫 번째 방정식은 세포 부피비 φ(x,t)의 시간·공간 변화를 기술하며, φ는 영양소 농도 c에 의존하는 성장·사멸 함수 G(c)와 확산·대류 항을 포함한다. 두 번째 방정식은 영양소 농도 c(x,t)의 확산·소모를 기술하며, 소모율은 φ와 직접적으로 연결된다. 수학적으로는
∂tφ = ∇·(Dφ∇φ) + G(c)φ,
∂tc = ∇·(Dc∇c) – λφc,
와 같은 형태가 된다. 여기서 Dφ, Dc는 각각 세포와 영양소의 유효 확산계수, λ는 영양소 소모 계수이다.
논문은 또한 경계조건을 정교하게 설정한다. 혈관 측면에서는 c = c₀(고정 농도)와 φ = 0(세포가 아직 침투하지 않음)을 적용하고, 외부 조직 측면에서는 무플럭스(Neumann) 조건을 부여해 폐쇄된 시스템을 구현한다. 차원less화 과정을 통해 주요 무차원 파라미터(예: 영양소 공급 비율, 성장 억제 임계치)를 도출하고, 이 파라미터가 시스템의 정성적 거동에 미치는 영향을 탐구한다. 수학적 분석에서는 정적 해(steady state) 존재성, 안정성, 그리고 파라미터 변화에 따른 bifurcation 가능성을 검토한다. 특히, 영양소 농도가 임계값 c* 이하로 떨어지면 G(c) < 0이 되어 세포 부피비가 감소하는 ‘핵심 사멸 영역’이 형성되며, 이는 종양 코어의 괴사 현상을 재현한다. 수치 시뮬레이션에서는 2D 직사각형 도메인에 혈관을 한쪽 면에 배치하고, 초기 φ는 혈관 근처에서만 작은 값으로 설정한다. 시간 전진 시 φ는 혈관에서 멀어질수록 감소하고, c는 확산에 의해 균일해지면서 φ와의 피드백 루프가 형성된다. 결과는 종양 코어가 혈관에서 일정 거리 이상 떨어지면 사멸하고, 혈관 근처에 얇은 증식층이 유지되는 ‘코드 형태’를 성공적으로 재현한다. 이러한 정량·정성 결과는 기존의 단일-상 모델과 비교해 조직 이질성 및 영양소-세포 상호작용을 보다 현실적으로 포착한다는 점에서 학문적·임상적 의의를 가진다.