주기적 하중이 연골 병변 형성에 미치는 영향 모델링 및 시뮬레이션

초록

본 논문은 기존의 반응‑확산‑지연 모델을 확장하여, 주기적 압축 하중이 관절 연골의 세포 사멸과 염증 반응에 미치는 영향을 수학적으로 구현한다. 사이클릭 로딩을 손실 항으로 포함시켜 살아있는 연골세포의 감소를 모사하고, 에리트로포이에틴(EPO) 농도의 차이가 병변 억제에 어떻게 작용하는지를 시뮬레이션 결과로 제시한다.

상세 분석

이 연구는 Graham 등(2012)의 반응‑확산‑지연 모델을 기반으로, 주기적 압축 하중을 “사이클릭 싱크(term)” 형태의 손실 함수 Γ(ε,U,r) 로 구현한다는 점에서 독창적이다. Γ는 변형률 ε와 ECM 농도 U에 의존하며, 로그 함수를 이용해 세포 사멸률을 비선형적으로 표현한다. 이는 실험적으로 보고된 균일 변형에 따른 세포 사망률을 근사화한 것으로, 변형률이 40 % 이상일 때 급격히 증가하는 특성을 반영한다. 모델은 원형 대칭을 가정하고 1‑D 반경 좌표 r에만 의존하도록 축소했으며, 로딩 영역(r ≤ 0.25 cm)과 비로드 영역을 명확히 구분한다.

세포군은 건강한 세포 C, 카타볼릭 세포 S_T, EPOR‑활성 세포 S_A, 괴사 세포 D_N, 그리고 사멸 세포 D_A 로 구분된다. 각 군은 상호 전이율(α, β₁, β₂, γ, μ, ν 등)과 지연 시간(τ₁, τ₂)을 통해 연결되며, 특히 S_T → S_A 전이는 TNF‑α에 의해 촉진되고, S_A → C 전이는 EPO에 의해 억제된다. 화학물질은 ROS(R), DAMPs(M), TNF‑α(F), EPO(P) 네 종류이며, 각각 확산(D_R, D_M, D_F, D_P)과 자연 감쇠(δ_R 등) 및 세포에 의한 생산(σ_R, σ_M, σ_F, σ_P) 항을 가진다. ECM(U)은 TNF‑α에 의해 분해되며, 그 분해는 DAMPs 방출을 유도해 양성 피드백 루프를 형성한다.

파라미터 추정은 문헌값을 기반으로 하며, 특히 δ_U, δ_F, δ_P 등은 반감기 데이터를 이용해 일일 단위로 변환하였다. Γ 함수의 상수 K_U와 p₀는 실험적 셀 사망률 데이터를 그대로 적용하였다. 시뮬레이션은 ε = 0.3, 0.4, 0.6, 0.8(%)와 두 가지 EPO 생산률(σ_P = 4.2×10⁻⁵, 3.3×10⁻³) 조합으로 수행되었다. 결과는 변형률이 클수록 건강한 세포 C가 급격히 감소하고, 염증 영역( S_T, S_A )이 확대되는 것을 보여준다. 반면 높은 σ_P 값에서는 EPO가 충분히 축적되어 EPOR‑활성 세포가 건강 세포로 복귀하고, 염증 확산이 억제된다. 이는 EPO가 TNF‑α에 의해 제한되는 메커니즘을 수학적으로 재현한 것으로, 치료적 EPO 보강이 병변 진행을 늦출 가능성을 시사한다.

모델의 제한점으로는 1‑D 축소, 세포 이동성 무시, 그리고 Γ 함수가 시간에 따라 감소하지 않는 점이 있다. 또한, 실제 연골의 복합 비선형 기계적 거동(비탄성, 유동성)과 조직 수준의 스트레스‑스트레인 분포를 완전히 반영하지 못한다. 향후 연구에서는 2‑D/3‑D 유한요소 해석과 실험적 변형‑사망률 데이터를 결합해 Γ를 동적 함수로 확장하고, EPO‑TNF‑α 상호작용을 더 정밀히 모델링할 필요가 있다.