만성 골수성 백혈병 모델의 안정성 분석: 간소화된 PDE 접근

초록

본 논문은 Roeder 등(2006)의 에이전트 기반 모델을 기반으로 한 혈액 생성 PDE 모델을 간소화하고, 그 간소화된 모델의 안정성을 수학적으로 분석한다. 간소화 과정과 파라미터 축소를 수치 시뮬레이션으로 검증한 뒤, 고유값 분석과 임계점(bifurcation) 연구를 통해 정상 혈구 생산과 CML(만성 골수성 백혈병) 진행 사이의 전이 메커니즘을 밝힌다.

상세 분석

본 연구는 먼저 Kim 등(2007)이 제시한 연속체 PDE 모델을 채택한다. 이 모델은 줄기세포(S), 전구세포(P), 분화세포(D) 세 군집을 연속적인 공간 변수와 시간 변수로 기술하며, 각 군집 간 전이율과 사멸율을 파라미터화한다. 원래 모델은 비선형 항이 복잡하고, 공간적 확산 항까지 포함해 8개 이상의 파라미터를 필요로 한다. 저자들은 생물학적 근거를 들어(예: 전구세포와 분화세포의 이동성 차이 무시, 확산 항의 규모가 작음) 이러한 비선형 항을 선형화하고, 확산 항을 제거함으로써 3개의 ODE 형태로 축소한다. 이 과정에서 핵심 동역학을 보존하는지 확인하기 위해 원 모델과 간소화 모델을 동일한 초기조건과 파라미터 집합으로 시뮬레이션하고, 시간 궤적과 고정점 위치를 비교하였다. 결과는 두 모델이 정성적으로 동일한 전이 현상을 보이며, 특히 정상 상태와 CML 상태 사이의 전이점이 거의 일치함을 보여준다.

수학적 분석 단계에서는 간소화된 시스템의 고정점을 찾고, 야코비안 행렬을 구성해 고유값을 계산한다. 정상 혈구 생산에 대응하는 고정점은 모든 고유값이 음의 실수부를 가져 안정적이며, CML에 대응하는 고정점은 하나 이상의 고유값이 양의 실수부를 가져 불안정함을 확인한다. 파라미터인 줄기세포 자기복제율 α와 전구세포 분화 억제율 β를 변동시키며, 두 파라미터가 특정 임계값을 초과할 때 고정점의 안정성이 전이(bifurcation)한다는 것을 발견한다. 특히, α가 증가하면 정상 고정점이 사라지고, β가 감소하면 CML 고정점이 안정화되는 서브크리티컬 피크스-노드 전이가 일어난다. 이러한 전이 현상은 실제 CML 환자에서 관찰되는 BCR‑ABL 유전자 발현 증가와 연관지어 해석될 수 있다.

또한, 저자들은 파라미터 민감도 분석을 수행해 α와 β가 모델 출력(예: 성숙 혈구 수)에 미치는 영향을 정량화한다. 결과는 α에 대한 민감도가 β보다 현저히 높으며, 이는 줄기세포 자기복제 조절이 치료 표적으로서의 중요성을 강조한다. 마지막으로, 간소화 모델이 원래 에이전트 기반 시뮬레이션과 거의 동일한 주기적 진동과 과도 현상을 재현함을 보여, 복잡한 미시적 상호작용을 거시적 연속체 모델로 효과적으로 대체할 수 있음을 입증한다.