인터페론과 리바비린 병용 치료를 통한 C형 간염 바이러스 동역학 모델링

초록

본 논문은 간세포, 감염세포, 감염성 바이러스, 비감염성 바이러스를 포함하는 4개의 연립 미분방정식으로 HCV 치료 과정을 모델링한다. 인터페론 단독 치료와 리바비린 병용 시 바이러스 부하의 2상·3상 감소를 설명하고, 약물 효율성 임계값을 도입해 효율이 이를 초과하면 바이러스가 소멸하고, 미달이면 낮은 새로운 평형 상태에 머무른다.

상세 분석

이 연구는 기존의 HCV 단일 약물 모델을 확장하여 인터페론(IFN)과 리바비린(RBV) 두 약물을 동시에 고려한 4차원 ODE 시스템을 제시한다. 변수는 T(정상 간세포), I(감염된 세포), V_I(감염성 바이러스), V_N(비감염성 바이러스)이며, 각각의 생산·소멸·감염 과정을 파라미터화한다. 인터페론은 감염세포의 사멸률을 증가시키고 바이러스 복제율을 감소시키는 효과를, 리바비린은 생산된 바이러스 중 일정 비율을 비감염성 형태(V_N)로 전환시키는 메커니즘을 수식에 반영한다.

모델은 다음과 같은 형태를 가진다:
dT/dt = s - d_T T - β T V_I,
dI/dt = β T V_I - δ I - ε_IFN I,
dV_I/dt = (1-ρ_RBV) p I - c V_I,
dV_N/dt = ρ_RBV p I - c V_N,
여기서 s는 간세포 생산, d_T는 자연 사멸, β는 감염률, δ는 감염세포 자연 사멸, ε_IFN은 인터페론에 의한 추가 사멸, p는 바이러스 생산, ρ_RBV는 리바비린에 의한 비감염성 전환 비율, c는 바이러스 청소율이다.

핵심적인 기여는 “약물 효율성 파라미터” η = (ε_IFN/δ)·(1-ρ_RBV) 를 정의하고, η_c 라는 임계값을 도출한 점이다. 선형 안정성 분석을 통해 η > η_c 일 때 감염성 바이러스 V_I와 감염세포 I가 모두 0으로 수렴하여 바이러스가 완전히 소멸함을 증명한다. 반대로 η < η_c 일 경우 새로운 비영(steady) 상태 (T*, I*, V_I*, V_N*)가 존재하며, 특히 V_I*는 치료 전 평형값보다 현저히 낮지만 완전 소멸은 아니다.

수치 시뮬레이션에서는 η가 임계값을 초과하면 바이러스 부하가 초기 급격한 1상 감소 후 완만한 2상 감소를 보이며, 경우에 따라 3상 감소(극히 낮은 잔류 바이러스)까지 관찰된다. 이는 임상에서 보고된 “biphasic” 및 “triphasic” 바이러스 감소 패턴을 정량적으로 재현한다. 또한, ρ_RBV가 증가할수록 η_c가 낮아져 리바비린의 보조효과가 인터페론 효율성을 크게 보강함을 확인한다.

이 모델은 파라미터 민감도 분석을 통해 β와 p가 치료 성공에 가장 큰 영향을 미치는 변수임을 밝히며, 환자별 파라미터 추정이 맞춤형 치료 전략 수립에 활용될 수 있음을 시사한다. 한계점으로는 면역 반응(특히 CD8⁺ T세포)의 동적 역할을 배제했으며, 간 조직 내 공간적 이질성을 고려하지 않은 점을 들 수 있다. 향후 연구에서는 면역학적 요소와 약물 동력학을 통합한 고차원 모델링이 필요하다.