노화가 나이브 T세포 수와 다양성에 미치는 영향에 대한 수학적 모델

초록

본 논문은 인간의 면역노화 현상을 설명하기 위해, 흉선에서의 신규 T세포 방출, 말초 증식 및 세포 사멸을 포함한 3가지 기본 메커니즘을 기반으로 한 미분방정식 모델을 구축한다. 시간에 따라 감소하는 흉선 방출률과 인구 의존적 사멸률을 고려해 전체 나이브 T세포 수와 클론 다양성을 정량화한다. 모델 분석 결과, 노화에 따른 다양성 감소는 주로 흉선 출력 감소와 작은 클론들의 소실에 기인함을 확인하였다.

상세 분석

이 연구는 면역노화의 핵심인 나이브 T세포 풀의 동태를 미시적 수준에서 기술한다. 총 세포 수 N(t)를 기술하는 기본 방정식 dN/dt = γ(t) + p N(t) – μ(N) N(t)에서, γ(t)=γ₀e^{–a t}는 흉선에서의 신규 세포 방출을, p는 말초에서의 자가 증식률을, μ(N)=μ₀+μ₁ N²/(N²+K²)는 IL‑7에 의한 포화 효과를 포함한 사멸률을 나타낸다. 파라미터 γ₀, a, p, μ₀, μ₁, K는 인간·쥐 실험 데이터를 통해 제한되었으며, 특히 γ₀·a^{–1}·K^{–1}와 μ₁/(p–μ₀)의 조합이 시스템의 장기 거동을 결정한다. μ₁/(p–μ₀)≥1이면 사멸이 증식을 억제해 안정적인 정착점 N_ss를 형성하고, γ₀·a^{–1}·K^{–1}가 작을 경우 증식‑주도(proliferation‑driven) 모드가, 크면 흉선‑주도(thymus‑driven) 모드가 나타난다.

다음으로, 클론별 세포 수 n_i(t)와 그 분포 P(n,t)를 도입해 다양성 지표 D_q = (∑_i n_i^q)^{1/(1–q)}(q≠1)와 같은 Rényi 엔트로피 기반 척도를 정의한다. 모델은 작은 클론(특히 n_i=1,2)의 소멸이 전체 다양성 감소에 비례함을 보여준다. 시뮬레이션 결과, 20‑30세 구간에서 평균 N≈10¹¹인 상태가 70‑80세에 약 50 % 감소하고, 동시에 D₁(Shannon)과 D₂(Simpson) 등 다양한 지표가 60‑80 % 수준으로 급격히 낮아진다. 이는 흉선 출력이 급격히 감소하면서 신규 클론이 거의 생성되지 못하고, 기존 작은 클론이 말초 증식에 의해 대형 클론으로 편향되는 현상을 수학적으로 재현한다.

또한, 샘플링 오류를 고려한 분석에서는 전체 풀에서 추출한 소규모 혈액 샘플이 실제 다양성을 크게 과소평가한다는 점을 강조한다. “보이지 않은 종” 문제를 해결하기 위해 Chao1, ACE 등 추정기를 적용했지만, 파라미터 불확실성 및 클론 크기 분포의 비대칭성 때문에 정확한 전체 다양성 복원은 여전히 어려움이 있다.

이러한 결과는 노화된 개인에서 백신 효능 저하와 감염 위험 증가를 설명하는 메커니즘을 제공한다. 특히, 작은 클론의 소실이 항원 특이적 반응 범위를 제한한다는 점은 면역 치료 설계 시 고려해야 할 중요한 요소이다. 모델은 파라미터 조정을 통해 IL‑7 보충, 흉선 재생 치료 등 잠재적 개입 전략의 정량적 효과를 예측할 수 있는 기반을 제공한다.