에피제네틱 상태 상속이 약물내성 퍼시스터를 통해 고형암 저항을 유도하는 메커니즘

초록

본 연구는 에피제네틱 상태가 세포 분열 시 어떻게 전달되는지를 확률적 상속 커널로 모델링하고, 이를 기반으로 약물내성 퍼시스터(DTP)의 발생·진화 과정을 에이전트 기반 시뮬레이션으로 재현한다. NSCLC 모델에 적용한 결과, 약물 투여 초기의 가역적 DTP가 장기 노출 시 지속적 저항(DRC)으로 전이되는 시점을 정확히 예측했으며, 가상 환자군을 이용한 인‑실리코 임상시험에서 적응형 치료 전략이 재발을 현저히 지연시킴을 보였다.

상세 분석

본 논문은 기존의 ODE·PDE 기반 평균장 모델이 포착하지 못하는 ‘저항 연속체’를 미시적 수준에서 구현하기 위해, Lei(2020) 가 제시한 확률적 에피제네틱 상속 프레임워크를 차용하였다. 각 세포는 (x₁, x₂, x₃)∈ℝ³ 로 표현되는 연속형 상태벡터를 갖으며, x₁은 표적 경로 활성도(약물 민감도), x₂는 줄기성·휴면 정도, x₃는 우회 신호 활성화를 의미한다. 세포 주기는 G0 정지기와 고정된 준비기(τ)로 구분되고, G0 단계에서 재진입 비율 β, 사멸·분화 비율 κ, 준비기 중 사멸 비율 μ가 각각 상태벡터에 함수 형태로 정의된다. 특히 β₀(x)와 κ(x)는 x₂에 대해 비단조적(β₀는 중간 x₂에서 최대, κ는 x₂가 증가할수록 감소)인 Hill‑형식으로 설계돼, 줄기성 높은 세포는 증식이 억제되고, 중간 수준에서는 증식이 촉진되는 생물학적 논리를 반영한다.

에피제네틱 상속은 확률밀도 p(x,y) 로 기술되며, 이는 모세포 상태 y에서 딸세포 상태 x로 전이될 확률을 정의한다. p는 평균이 y이고 분산 σ²를 갖는 다변량 정규분포 형태를 취해, 세포 분열 시 ‘기억 손실’ 혹은 ‘재배열’ 정도를 조절한다. 이 커널 덕분에 모델은 DTP가 약물 억제 하에서 x₁이 감소하고 x₂가 상승하는 방향으로 drift하고, 장기 노출 시 x₃가 점차 증가해 영구 저항(DRC)으로 전이되는 과정을 자연스럽게 재현한다.

시뮬레이션은 에이전트 기반 알고리즘으로 구현돼 개별 세포의 계통추적이 가능하며, 고차원 IDE를 직접 수치해석하는 것보다 계산 효율성이 높다. 파라미터는 NSCLC PC9 세포주에 대한 erlotinib 실험 데이터를 기반으로 β₀, κ₀, a₁a₃, κ₁, nₖ 등 여러 생물학적 지표를 추정하였다. 결과적으로, 약물 투여 첫 2주 내에 전체 세포의 510%가 DTP로 전이되고, 이후 4~6주 사이에 x₃가 임계값을 초과하면서 DRC 비율이 급격히 상승하는 ‘재발 곡선’이 관찰되었다.

인‑실리코 임상시험에서는 가상 환자 200명을 무작위 배정해 표준 고정 용량 치료와, 종양 부피·DTP 비율에 기반한 적응형 용량 조절 전략을 비교하였다. 적응형 전략은 약물 강도를 단계적으로 낮추었다가 DTP 비율이 일정 수준을 초과하면 다시 강화하는 방식으로, 평균 progression‑free survival이 표준 치료 대비 약 45% 연장되는 효과를 보였다. 이는 에피제네틱 플라스틱성이 치료 반응에 미치는 영향을 정량화하고, 실제 임상에서 치료 스케줄을 최적화하는 데 유용한 정량적 근거를 제공한다.

한계점으로는 모델이 미세환경(예: CAF, TAM)에서 분비되는 사이토카인 등을 명시적으로 포함하지 않았으며, x₁·x₃ 간의 상호 억제 관계를 단순화한 점을 들 수 있다. 향후에는 싱글‑셀 전사체 데이터와 결합해 상속 커널을 데이터‑드리븐 방식으로 추정하고, 공간적 확산을 포함한 하이브리드 PDE‑ABM 모델로 확장하는 연구가 필요하다.