혈액 생성 과정의 수학적 분석과 만성 골수성 백혈병 적용

초록

본 논문은 골수 내 혈액세포 생산을 연령 구조 편미분 방정식 두 개로 모델링하고, 이를 연령에 대해 적분해 분포 지연을 포함한 두 비선형 미분방정식으로 변환한다. 라우프노프 함수 구축을 통해 유일한 평형점이 존재할 경우 영평형이 전역적으로 안정함을 증명하고, 비영평형(생물학적으로 의미 있는 평형)은 파라미터 변화에 따라 Hopf 분岐를 일으켜 주기적 진동을 발생시킬 수 있음을 보인다. 수치 시뮬레이션으로 이론적 결과를 검증한다.

상세 분석

이 연구는 혈액세포의 생성·분화 과정을 수학적으로 정량화하기 위해 연령 구조 모델을 채택한다. 구체적으로, 골수 내 미성숙 세포와 성숙 세포 각각에 대해 연령 a와 시간 t에 대한 편미분 방정식(∂/∂t + ∂/∂a) 형태를 설정하고, 세포주기의 평균 지속시간을 분포 지연(kernel) 형태로 도입한다. 연령에 대해 전체 세포수를 적분하면 두 개의 비선형 미분방정식이 도출되며, 여기서 비선형성은 세포 간의 피드백(예: 성장인자 억제)으로부터 발생한다.

주요 수학적 기법은 라우프노프 함수(Lyapunov functional) 구성이다. 저자들은 시스템이 영평형(모든 세포가 사멸)만을 갖는 경우, 적절히 정의된 라우프노프 함수를 통해 dV/dt ≤ 0 를 증명하고, V가 영점에서만 최소가 되므로 전역적 안정성을 확보한다. 이는 실제 생리학적 상황에서 골수 기능이 완전히 정지될 경우(예: 방사선 치료 후) 세포 수가 회복되지 않음을 의미한다.

비영평형에 대해서는 고정점 해를 구하고, 선형화 후 특성 방정식에 분포 지연 항이 포함된 형태가 나타난다. 지연 커널의 평균값과 분산이 증가하면 특성 방정식의 근이 복소수 쌍으로 이동하여 실수부가 양수가 되는 Hopf 임계점에 도달한다. 저자들은 임계 파라미터 값을 분석하고, Hopf 분岐가 초임계 영역에서 안정적인 주기해(limit cycle)를 생성함을 보인다. 이는 만성 골수성 백혈병(CML) 환자에서 관찰되는 장기적인 혈액 수치 진동을 설명하는 메커니즘으로 해석될 수 있다.

수치 실험에서는 실제 임상에서 보고된 세포주기(약 24~48시간)와 말초 혈액 진동 주기(수주에서 수개월)를 매개변수로 설정한다. 시뮬레이션 결과, 짧은 세포주기 지연에도 불구하고 분포 지연의 폭이 충분히 크면 시스템이 Hopf 분岐를 겪어 장기적인 진동을 보인다. 이는 “짧은 내부 주기 → 긴 외부 진동” 현상을 수학적으로 정당화한다. 또한, 파라미터 민감도 분석을 통해 성장인자 억제 강도와 세포 사멸률이 진동 발생에 결정적인 역할을 함을 확인한다.

이 논문은 연령 구조 PDE → 분포 지연 ODE 변환, 라우프노프 안정성 증명, Hopf 분岐 분석이라는 세 단계의 수학적 절차를 통해 혈액 생성 메커니즘을 체계적으로 탐구한다. 특히, 분포 지연을 명시적으로 모델에 포함시킨 점이 기존의 단일 지연 모델보다 생리학적 현실성을 크게 향상시킨다. 향후 연구에서는 다중 세포계(예: 골수 미성숙 세포, 혈관 내피 세포)와 공간 확산 효과를 추가하여 보다 복합적인 혈액 질환 모델링이 가능할 것으로 기대된다.