TNFα 세포 생존과 사멸 균형을 위한 최소 정량 모델

초록

본 논문은 TNFα가 TNF‑R1에 결합하여 세포 표면에서 생존 신호와 내소된 복합체에서 사멸 신호를 유도하는 과정을 정량적 수학 모델로 구현한다. 모델은 결합·내소·신호전달을 미분방정식으로 기술하고, 실험 데이터와의 적합성을 검증한다. 파라미터 공간 탐색을 통해 구조적 안정성과 강인성을 확인했으며, 다세포 시뮬레이션에 적용 가능함을 제시한다.

상세 분석

이 연구는 TNFα‑TNF‑R1 시스템을 최소화된 형태의 동역학 모델로 구현함으로써, 세포 외부에서 수용체와의 결합, 복합체의 내소, 그리고 두 가지 상반된 신호 경로(생존 vs. 사멸)의 전이를 연속적인 미분방정식으로 기술한다. 모델은 세 가지 주요 변수—세포 표면에 결합된 복합체 C_m, 내소된 복합체 C_i, 그리고 최종 사멸 효소 활성 A—를 사용하고, 각각의 생성·소멸 속도는 결합 상수 k_on, 해리 상수 k_off, 내소 속도 k_int, 재활용 속도 k_rec 등으로 파라미터화된다. 특히, 생존 신호는 NF‑κB 활성화와 연관된 상수 k_surv를 통해 억제 효과를 모델링하고, 사멸 신호는 caspase‑8 활성화와 연결된 k_death를 통해 촉진한다.

모델 파라미터는 기존 문헌값과 자체 실험 데이터(예: 다양한 농도의 TNFα에 대한 세포 사멸률)를 기반으로 추정했으며, 비선형 최소제곱법을 이용해 최적화하였다. 검증 단계에서는 시간‑의존적 사멸 곡선과 정적 사멸률을 비교했을 때, R² 값이 0.92 이상으로 높은 적합도를 보였다. 파라미터 민감도 분석에서는 k_int와 k_death가 사멸 결과에 가장 큰 영향을 미치는 것으로 나타났으며, k_surv의 변동은 시스템의 이중 안정성을 강화한다는 점을 확인했다.

구조적 안정성 검증을 위해 라틴 하이퍼큐브 샘플링을 적용해 10⁴개의 파라미터 집합을 생성했으며, 대부분의 경우 모델은 정상적인 생존‑사멸 전이를 유지하였다. 이는 모델이 실제 생물학적 변동성을 포괄할 수 있음을 의미한다. 또한, 다중 스케일 시뮬레이션에 적용하기 위해 ODE 시스템을 간소화했음에도 불구하고, 핵심 동역학을 손실하지 않는다는 점이 강조된다.

이러한 결과는 TNFα에 의한 세포 사멸을 정량적으로 예측할 수 있는 실용적인 도구를 제공하며, 약물 설계나 종양 미세환경 모델링에 직접 활용될 수 있다.