중심 대식세포가 조절하는 적혈구 생산: 하이브리드 모델을 통한 erythroblastic island 분석

초록

본 논문은 골수 내 적혈구 전구세포와 중앙 대식세포가 형성하는 ‘erythroblastic island’를 다중스케일 하이브리드 모델로 구현한다. 세포 내 Erk‑Fas 단백질 경쟁을 ODE로, 세포 간 물리·화학 상호작용을 개별 에이전트 기반의 이산 모델로, 성장인자 확산을 PDE로 결합한다. 모델 시뮬레이션과 수학적 안정성 분석을 통해, 대식세포가 생존 시토카인을 분비할 경우 섬의 구조적 안정성이 크게 향상되고, 외부 피드백에 대한 반응성이 높아짐을 확인한다.

상세 분석

이 연구는 세 가지 스케일을 통합한 하이브리드 프레임워크를 제시한다. 첫 번째는 세포 내 수준으로, Erk와 Fas라는 두 핵심 단백질의 상호 억제 관계를 비선형 ODE(식 2.1‑2.2)로 기술한다. Erk는 에리트로포이에틴(Epo)과 기타 성장인자(GF)의 활성화와 자체 활성화(β·k)를 통해 증식 신호를 전달하고, Fas는 주변 적혈구(reticulocyte)에서 방출되는 Fas‑리간드(FL)에 의해 활성화된다. 두 단백질의 농도는 각각 0‒1 사이의 정규화된 변수(E, F)로 표현되며, 상호 억제(b·F, c·E)와 선형 분해(a·E, d·F)를 포함한다. 파라미터 α, β, γ 등은 외부 신호 강도에 따라 조정 가능해, 세포가 자기 재생(self‑renewal), 분화(differentiation), 혹은 세포자멸(apoptosis) 중 하나의 운명을 선택하도록 만든다.

두 번째 스케일은 개별 세포를 입자처럼 다루는 이산 모델이다. 각 세포는 위치와 반경을 갖고, 주변 세포와 물리적 충돌·접촉을 통해 신호 교환을 수행한다. 세포 분열 시에는 무작위로 결정되는 세포주기 길이와 유사분열면의 방향성을 도입해, 실제 골수 조직에서 관찰되는 변동성을 재현한다.

세 번째 스케일은 성장인자와 사이토카인의 확산·소멸을 기술하는 PDE(확산‑반응 방정식)이다. 특히 대식세포가 분비하는 ‘생존 사이토카인’(예: IL‑6, TGF‑β 등)의 농도장은 공간적으로 계산되며, 이는 주변 전구세포의 α와 γ 파라미터에 직접적인 영향을 미친다.

모델 결합은 다음과 같이 이루어진다. 세포가 경험하는 외부 신호(Epo, GF, FL)는 PDE 해에서 추출된 국소 농도값을 통해 ODE 파라미터에 매핑된다. 반대로, 세포가 분비하거나 소모하는 물질은 PDE의 소스/소멸 항에 반영된다. 이렇게 순환적인 피드백 루프를 통해 전구세포 집단의 동역학이 전반적인 섬 구조와 연동된다.

수학적 분석에서는 고정점 존재와 안정성을 검증한다. Erk‑Fas 시스템은 파라미터 조합에 따라 세 개의 고정점(고Erk/저Fas, 저Erk/고Fas, 중간 상태)을 가질 수 있으며, 각각이 세포 운명에 대응한다. Jacobian 행렬의 고유값을 통해 고정점의 지역 안정성을 평가했으며, 특히 대식세포가 제공하는 고농도 사이토카인 영역에서는 고Erk 고정점이 넓은 매개변수 구역에서 안정적으로 유지됨을 확인했다.

시뮬레이션 결과는 두 가지 시나리오를 비교한다. (1) 대식세포가 없는 섬: 초기 전구세포는 빠르게 분화·사멸하고, 섬 구조가 파괴되어 세포 수가 급감한다. (2) 대식세포가 중심에 존재하고 생존 사이토카인을 지속적으로 방출하는 경우: 전구세포는 지속적인 자기 재생을 유지하며, 섬이 장기간 안정적으로 존재한다. 또한, 외부 Epo 농도 변화에 대한 섬의 반응 속도와 진폭이 대식세포가 있을 때 현저히 개선된다. 이는 대식세포가 ‘핵심 조절자’로 작용해 섬의 내재적 안정성을 부여하고, 전신 피드백(혈액 내 적혈구 수에 따른 Epo 조절)과의 연계를 강화한다는 생물학적 의미를 시사한다.

이러한 결과는 기존의 비공간적, 결정론적 혈액세포 모델이 놓치고 있던 미세환경의 역할을 정량적으로 입증한다. 특히, 대식세포가 제공하는 국소 성장인자와 세포 간 물리적 접촉이 전구세포의 운명 결정에 미치는 영향을 동시에 고려함으로써, 골수 내 ‘니치(niche)’ 개념을 수학적으로 구체화하였다.