상처 치유 시 MAPK 파동을 설명하는 최소 확산‑대류 모델

초록

본 논문은 상피 상처 치유 과정에서 관찰되는 두 차례의 MAPK 파동을, ROS와 EGFR 리간드의 경쟁 및 세포 이동에 의한 기계적 스트레스를 포함한 최소 확산‑대류 모델로 재현한다. 모델은 MAPK 신호전달계의 이중안정성을 이용해 첫 번째 급성 파동과 두 번째 지속 파동을 각각 전파·반사와 느린 전진 형태로 설명하고, 수학적 해석을 통해 파동 속도와 존재 조건을 제시한다.

상세 분석

이 연구는 상피 세포층에서 손상 직후 발생하는 급격한 MAPK 활성화와, 상처가 닫히는 동안 지속되는 느린 MAPK 전파라는 두 가지 현상을 통합적으로 설명하려는 시도이다. 핵심 가정은 EGFR이 두 종류의 외부 자극, 즉 즉각적으로 생성되는 ROS와 지속적으로 분비되는 성장인자(EGF 등)와 경쟁적으로 결합한다는 점이다. ROS는 손상 부위에서 급격히 증가하여 EGFR을 일시적으로 활성화시키고, 이는 MAPK 카스케이드를 급속히 촉발한다. 이때 모델은 ROS의 확산과 소멸을 1차원 확산 방정식으로 기술하고, ROS‑EGFR 복합체 형성을 비선형 활성화 함수로 표현한다.

첫 번째 파동은 ROS 농도의 급격한 공간적 구배에 의해 전파되며, 경계에서 반사되는 특성을 보인다. 수학적으로는 반사 경계 조건을 갖는 확산‑대류 방정식의 해가 전파‑반사 파동으로 나타난다. 두 번째 파동은 세포 이동에 의해 발생하는 기계적 스트레스를 ROS와 유사한 신호원으로 가정한다. 세포가 상처를 메우면서 발생하는 장력은 EGFR에 추가적인 활성화를 제공하고, 이는 MAPK 신호를 지속적으로 유지한다. 여기서 모델은 세포 이동 속도를 대류 항으로, 스트레스‑유도 활성화를 또 다른 비선형 함수로 포함한다.

MAPK 카스케이드는 양의 피드백을 포함한 이중안정 시스템으로 모델링된다. 즉, 낮은 MAPK 수준에서는 활성화가 억제되고, 일정 임계값을 초과하면 높은 활성화 상태로 전이한다. 이 이중안정성은 파동 전파를 가능하게 하는 핵심 메커니즘이며, 수학적으로는 반응‑확산 방정식의 퍼텐셜 함수 형태로 나타난다. 저자들은 파라미터 공간을 탐색하여 두 파동이 동시에 존재할 수 있는 조건을 도출하고, 파동 속도는 활성화 함수의 기울기와 확산 계수, 대류 속도에 의해 결정된다는 것을 보였다.

또한, 저자들은 traveling wave 해를 가정하고 위상 평면 분석을 수행해 파동의 존재와 안정성을 증명한다. 특히, 첫 번째 파동은 높은 ROS 확산계수와 빠른 대류 속도에 의해 빠르게 소멸하고, 두 번째 파동은 낮은 확산계수와 지속적인 대류(세포 이동)로 인해 느리게 전파된다. 이러한 수학적 결과는 실험에서 관찰된 급성·지속 파동의 시간·공간적 특성과 일치한다.

결과적으로 이 모델은 복잡한 생물학적 네트워크를 최소한의 변수와 방정식으로 압축하면서도, ROS‑EGFR 경쟁, 기계적 스트레스, MAPK 이중안정성이라는 세 가지 핵심 메커니즘을 통해 상피 상처 치유 과정의 동역학을 정량적으로 설명한다.