고열치료와 단백질 분해 억제제 병용의 최적 전략 모델링

초록

본 연구는 기존 열충격 반응(HSR) 모델을 확장해 확률적 요소를 도입하고, 열충격 단백질(HSP)에 의한 열내성(thermotolerance)을 정량화한다. 또한, 임상에서 사용되는 프로테아좀 억제제 보르테조밉과 고열치료를 병합한 다중 치료 전략의 효능을 시뮬레이션으로 평가하여 최적의 투여 시점과 강도를 제시한다.

상세 분석

이 논문은 암 치료에 고열치료를 보조적으로 활용할 때 발생하는 열내성 현상을 수학적으로 규명하고, 이를 바탕으로 보르테조밉과의 병용 전략을 최적화한다는 두 가지 목표를 갖는다. 먼저 기존의 결정론적 HSR 모델을 stochastic differential equation 혹은 Gillespie 알고리즘 기반의 마코프 과정으로 확장함으로써, 세포 내 HSP 발현 및 분해 과정의 본질적인 변동성을 반영한다. 모델은 mRNA 전사, HSP 단백질 합성, 그리고 HSP와 손상된 단백질 간의 결합/해리 반응을 포함하며, 각 반응률은 온도 의존적인 Arrhenius 형태로 정의된다.

핵심적인 개념인 ‘HSP‑유도 열내성’은 고열 후 HSP 농도가 일정 수준 이상 유지될 때, 추가적인 열 스트레스에 대한 세포의 사멸률이 감소하는 현상으로 정량화된다. 이를 위해 저자들은 열내성 지속시간(thermotolerance window)과 강도(intensity)를 각각 HSP 농도 임계값 초과 시간과 사멸률 감소 비율로 정의하고, 시뮬레이션을 통해 온도 프로파일(예: 42 °C 60 min)과 회복 시간(0–24 h) 사이의 관계를 매핑하였다. 결과는 고열 후 4–8 h 사이에 열내성이 최고조에 이르며, 12 h 이후 급격히 소멸한다는 것을 보여준다.

다음으로 보르테조밉의 약동학·약력학(PK/PD) 모델을 HSR 모델에 통합한다. 보르테조밉은 프로테아좀 억제로 인해 비정상 단백질이 축적되고, 이는 HSP의 결합 파트너를 증가시켜 열내성을 간접적으로 강화한다. 저자들은 보르테조밉 농도 C(t)에 대해 Michaelis‑Menten 형태의 억제 효과와, 억제된 프로테아좀에 의해 발생하는 추가적인 단백질 손상을 ΔP(t)로 모델링하였다. 이때 ΔP(t)는 HSP와 경쟁적으로 결합하여 사멸률을 상승시키는 역할을 한다.

최적화 문제는 ‘고열 → 보르테조밉’ 혹은 ‘보르테조밉 → 고열’ 순서, 각 치료 간 간격, 그리고 각각의 강도를 변수로 두고, 전체 세포 사멸률(또는 종양 부피 감소)을 목표 함수로 설정한다. 다목적 유전 알고리즘을 이용한 탐색 결과, 고열 후 2 h 이내에 보르테조밉을 투여하는 전략이 열내성에 의한 보호 효과를 최소화하면서, 두 치료의 시너지 효과를 극대화한다는 것이 도출되었다. 또한, 고열 강도를 41.5 °C에서 43 °C 사이로 조절하고, 보르테조밉 용량을 0.5–1 µM 범위로 유지할 때 가장 높은 치료 효능을 보였다.

이러한 모델링은 실험적 데이터(예: HSP70 Western blot, 세포 생존률 측정)와의 정량적 일치를 보이며, 열내성의 시간‑강도 곡선을 실험적으로 검증할 수 있는 프레임워크를 제공한다. 한계점으로는 조직 수준의 혈류 및 열 전달 효과를 단순화했으며, 종양 이질성 및 면역 반응을 포함하지 않은 점을 들 수 있다. 향후 모델에 다중 세포 유형, 혈관 네트워크, 그리고 면역 세포 활성화를 통합한다면, 임상 전이 단계에서 보다 정교한 치료 스케줄링이 가능할 것이다.