상피세포 내 자유 라디칼과 항산화제 상호작용의 수학·컴퓨팅 모델링

초록

본 연구는 자유 라디칼과 항산화제의 동역학을 설명하기 위해 수학적 미분방정식 모델과 C++ 기반 컴퓨팅 시뮬레이션을 각각 상피세포 배양과 개인 수준에 적용하였다. 수학 모델은 2006년 Anne Hanneken 연구의 실험 데이터를 통해 검증되었으며, 컴퓨팅 모델은 노화 네트워크 이론을 도입해 개체 전체에서의 산화 스트레스와 항산화 방어 메커니즘을 재현한다. 두 모델 모두 기존 선형·서술적 접근의 한계를 극복하고 복합적인 생물학적 시스템을 정량적으로 분석하는 틀을 제공한다.

상세 분석

이 논문은 전통적인 생물학 서술 방식이 복잡한 산화·항산화 시스템을 충분히 포착하지 못한다는 점을 지적하고, 이를 보완하기 위한 두 가지 정량적 모델을 제시한다. 첫 번째는 연속 미분방정식으로 구성된 수학적 모델로, 자유 라디칼(R)과 항산화제(A) 사이의 생성·소멸, 결합 반응을 질량 보존 원칙에 따라 기술한다. 파라미터는 실험적으로 측정된 라디칼 생성률, 항산화제 재생 속도, 반응 상수 등을 사용하며, 초기 조건은 상피세포 배양 시점의 ROS 농도와 항산화 효소 활성을 반영한다. 모델은 0‑D(시간에만 의존) 형태이지만, 비선형 항을 포함해 피드백 억제와 포화 효과를 구현함으로써 실제 세포 내 동역학을 근사한다. 검증 단계에서는 Anne Hanneken 팀이 2006년 발표한 안구 상피세포의 ROS 측정 데이터를 활용해 시뮬레이션 결과와 실험값의 RMSE를 최소화하였다. 결과는 라디칼 피크 시점, 항산화제 고갈 및 회복 과정을 정밀히 재현했으며, 파라미터 민감도 분석을 통해 특정 효소(예: SOD, catalase)의 역할이 전체 시스템 안정성에 미치는 영향을 정량화했다.

두 번째는 C++로 구현된 에이전트 기반 컴퓨팅 모델이다. 여기서는 각 세포를 노드, 세포 간 신호 전달을 엣지로 보는 네트워크 구조를 채택하고, 노드마다 라디칼 및 항산화제 농도를 상태 변수로 둔다. 노드 간 상호작용은 ‘노화 네트워크 이론’에 기반해, 노드의 손상 정도가 증가할수록 연결 강도가 약화되고, 손상 신호가 주변 노드에 전파되는 확률이 높아지는 형태로 설계되었다. 이 모델은 개인 전체(예: 조직 수준)에서 산화 스트레스가 어떻게 전파되고, 항산화 방어가 지역적으로 차단되는지를 시뮬레이션한다. 파라미터는 연령에 따른 항산화 효소 감소율, 미토콘드리아 손상 축적, 세포 사멸 임계값 등을 포함한다. 시뮬레이션 결과는 노화 진행에 따라 ROS 파동이 조직 전반에 퍼지는 양상을 보여주며, 특정 부위(예: 망막)에서 항산화 보충이 전체 네트워크 안정성을 크게 향상시킬 수 있음을 시사한다.

두 모델 모두 선형 논리와 서술적 접근이 놓치는 비선형 피드백, 임계 현상, 공간적 전파 메커니즘을 포착한다는 점에서 학문적·실용적 의의를 가진다. 특히 수학 모델의 실험 검증과 컴퓨팅 모델의 가상 실험을 결합함으로써, 실험적 제한을 넘어선 가설 검증과 치료 전략 설계가 가능해진다. 향후 모델 확장은 다층 네트워크(세포‑조직‑기관)와 다변량 오믹스 데이터 통합을 통해 개인 맞춤형 항산화 치료 시뮬레이션으로 이어질 전망이다.