두 종류 세포 분열 과정의 정확 해와 암 진행 모델링

초록

이 논문은 일방향 돌연변이를 포함한 두 종류(정상 A세포와 변이 B세포) 연속시간 출생-사망(branching) 과정을 정확히 풀어, 임의의 시점에서의 생성함수를 구하고, 다양한 대시간·소돌연변이 스케일링 한계에서의 거동을 분석한다. 특히 변이 세포 수 분포가 파워‑법 꼬리를 보이며 평균이 발산하는 비정상적 특성을 보이는 것을 밝혀, 암 치료 저항성 및 종양 진행의 수학적 이해에 기여한다.

상세 분석

본 연구는 Kendall이 제안한 두 종류 세포의 출생‑사망 과정(식 (2))을 일반화하여, A세포가 일정 비율(α₁=1)로 증식·소멸하고, 일정 비율(ν)로 불가역적으로 B세포로 변이되는 모델을 다룬다. λ₁=1−β₁−ν, λ₂=α₂−β₂ 로 정의된 ‘적합도’ 파라미터는 각각 A와 B 세포의 순증식률을 나타내며, 부호에 제한을 두지 않는다. 저자들은 후방 Kolmogorov 방정식(10a,10b)을 이용해 두 생성함수 A(x,y,t), B(x,y,t)를 도출하고, B에 대한 식은 독립적으로 해석 가능함을 보여준다(B식 (12)). A에 대한 방정식은 Riccati 형태가 되며, 이를 Sturm‑Liouville 형태로 변환해(14‑15) 일반 해를 hypergeometric 함수로 표현한다(식 (25)‑(27)). 핵심 파라미터 ω는 2차 방정식(16)으로부터 얻으며, C 상수는 초기조건(11a)으로부터 결정된다(28‑29). 결과적으로 A(x,y,t)는 복합적인 hypergeometric 함수와 그 파생함수들의 조합으로 정확히 기술된다.

특수 경우로 λ₁=0( A세포가 임계), λ₂=0( B세포가 임계) 혹은 양쪽이 모두 임계인 경우를 별도로 분석한다. 특히 양쪽이 임계인 경우(λ₁=λ₂=0)에는 Bessel 함수(Iₙ, Kₙ)를 이용한 간단한 형태(35‑36)로 축약된다. 이는 기존에 ‘해결 불가능’하다고 여겨졌던 모델이 실제로는 정확히 풀릴 수 있음을 보여준다.

대시간 스케일링 분석에서는 생존 확률 S(t)와 세포 수의 전형적인 성장 스케일 χ(t)를 정의하고, x→1 한계에서 생성함수가 라플라스 변환으로 전환됨을 이용해 확률밀도 p(ṁ)의 구조를 파악한다. 임계 경우에는 S_A(t)∼1/t, 전체 생존 확률 S(t)∼const·t⁻¹/⁄(α₂ν) 로서 B세포가 주된 생존자를 차지함을 확인한다. 초과임계(λ₂>0)에서는 최종 생존 확률 s_∞=λ₁/2+√(λ₁²/4+νλ₂α₂) 로 주어지며, 이는 B세포가 우세하게 되는 상황을 정량화한다. 또한, 작은 돌연변이율과 긴 시간의 동시극한에서 변이 세포 수 분포가 파워‑법 꼬리를 갖고 평균이 발산한다는 비정상적 현상을 발견한다. 이는 암 치료 저항성 세포가 희소하지만 장기적으로는 지배적인 역할을 할 수 있음을 수학적으로 뒷받침한다.

이러한 정확 해와 스케일링 결과는 실험 데이터(예: 유도형 라벨링을 통한 클론 크기 분포)와 직접 비교할 수 있는 기반을 제공한다. 또한, 변이가 세포 분열 시점에 일어나는 경우(식 (4))에 대한 해도 부록 A에 제시되어, 실제 생물학적 상황에 맞는 모델 선택이 가능하도록 한다. 전반적으로 복잡한 다형성 종양의 진화, 치료 저항성 발생, 전이 형성 등을 확률론적·분기 과정 관점에서 정밀하게 기술한 연구라 할 수 있다.