단일 세포와 3D 배양에서 HIF1α 동역학 모델링

초록

본 연구는 저산소 충격을 받은 단일 세포와 3차원 구형체(spheroid)에서 HIF‑1α 단백질의 일시적 변화를 수학적으로 모델링한다. 피드백 억제 회로인 HIF‑1α와 PHD(프로릴 하이드록실라제) 사이의 상호작용을 미분방정식으로 기술하고, 파라미터 변이를 통해 세포 간 이질성을 반영한다. PHD2가 주요 조절자임을 확인하고, 반복 저산소 펄스에 대한 예측을 제시한다. 또한 구형체 내부의 산소 프로파일을 실험적으로 측정해 모델에 통합함으로써, HIF‑1α와 그 전사 표적의 비단조적 응답을 설명한다.

상세 분석

이 논문은 저산소 환경에서 세포 내 신호 전달이 시간적 패턴을 통해 세포 운명을 결정한다는 최신 개념을 수학적 모델링으로 구체화한다. 핵심은 HIF‑1α와 산소 감지 효소인 PHD 사이의 부정적 피드백 루프를 2차 미분방정식 형태로 구현한 점이다. 모델은 HIF‑1α의 합성, PHD에 의한 프로리놀화 및 그에 따른 프로테아좀 매개 분해를 각각 독립적인 속도 상수(k_syn, k_oh, k_deg)로 분리하고, PHD 활성이 산소 농도와 HIF‑1α 농도에 비례적으로 증가한다는 가정을 둔다. 이를 통해 저산소 충격 직후 HIF‑1α가 급격히 축적되었다가 PHD가 활성화되면서 급격히 감소하는 전형적인 ‘피크‑감쇠’ 형태를 재현한다.

단일 세포 데이터의 이질성을 설명하기 위해 파라미터 분포를 도입했는데, 이는 베이즈 추정 혹은 최대우도 추정으로 각 세포별 최적 파라미터를 도출함으로써 가능했다. 특히 PHD2의 발현 수준과 효소 활성도가 다른 세포군에서 HIF‑1α 피크 높이와 지속 시간이 크게 달라짐을 확인하였다. 이는 PHD2가 HIF‑1α 동역학의 주된 조절자임을 실험적으로 뒷받침한다.

모델을 3D 구형체에 확장할 때는 산소 확산‑소비 방정식을 풀어 구형체 중심부와 외부의 산소 농도 구배를 얻었다. 이 산소 프로파일을 시간에 따라 변하는 입력값으로 사용해, 구형체 내부 세포들의 HIF‑1α 축적이 단순히 평균 산소 농도에 비례하지 않고, 급격한 저산소 파동에 민감하게 반응한다는 비단조적 현상을 설명한다. 결과적으로 HIF‑1α의 순간적 피크는 전사 표적인 VEGF, GLUT1 등과 시간 차이를 두고 발현되며, 이는 저산소 펄스의 빈도와 지속 시간에 따라 서로 다른 유전자 프로그램을 유도한다는 중요한 생물학적 함의를 제공한다.

이 모델은 실험적 검증을 위해 반복 저산소 펄스를 적용한 후 HIF‑1α와 PHD2의 동시 측정, 그리고 PHD2 억제제(예: DMOG) 처리 실험을 제안한다. 또한, 파라미터 민감도 분석을 통해 k_oh(산소 의존성 하이드록실화 속도)와 k_deg(분해 속도)의 변화가 피크 형태와 지속 시간에 미치는 영향을 정량화함으로써, 치료적 표적으로서 PHD2 조절의 가능성을 제시한다.