다세포 신호 전달 시스템의 특이 섭동 기반 모델 축소

초록

수화젤 비드 안에 분포된 송신·수신 세포들의 확산 매개 상호작용을 PDE‑ODE 결합 형태로 모델링하고, 확산이 빠른 시간 스케일을 갖는다는 사실을 이용해 특이 섭동 이론으로 차원을 축소한다. 구형 그린함수를 이용해 유도된 통신 행렬이 다중 에이전트 네트워크 형태의 유한 차원 모델을 제공하며, 수치 실험에서 원본 모델과 거의 동일한 동작을 보이면서 수천 개 셀의 시뮬레이션을 효율적으로 수행할 수 있음을 확인한다.

상세 분석

본 논문은 수화젤 비드 내부에 임베딩된 이질적인 세포 집단을 ‘송신‑수신’ 구조로 설정하고, 세포가 방출·흡수하는 확산성 신호분자를 3차원 확산 방정식과 세포 내부 비선형 ODE를 결합한 PDE‑ODE 시스템으로 기술한다. 핵심 가정은 세포 크기가 비드 반경에 비해 무시할 수 있을 정도로 작아, 각 세포를 점 소스로 모델링한다는 점이다. 이를 통해 (4)식의 확산 PDE는 Dirac 델타 함수를 이용한 점원천 항을 포함하고, 경계는 흡수형(Dirichlet)으로 설정한다.

시간 스케일 분석에서는 확산 시간 상수 τ_v = L²π⁻²/D와 막 투과 시간 상수 τ_u = 1/α, 그리고 세포 내부 반응 시간 τ_x를 정의하고, 실험적 파라미터(L≈300 µm, D≈10⁴ µm²/min, α≈0.1–100 min⁻¹, τ_x≈10–100 min)로부터 τ_v, τ_u ≪ τ_x임을 확인한다. 이로써 ε_v = τ_v/τ_x, ε_u = τ_u/τ_x이 작은 매개변수임을 가정하고, 시스템을 표준 특이 섭동 형태(9)–(11)로 재스케일링한다.

빠른 확산·막 교환 서브시스템에 대한 지수 안정성은 두 단계로 증명된다. 먼저 (13)식은 선형 ODE 형태이므로 지수 감쇠가 명백하고, (12)식은 에너지 방법을 이용해 L² 노름에 대한 미분 불등식 d/dt‖w‖² ≤ -2c‖w‖² (c≥λ₁, λ₁은 라플라시안의 첫 번째 디리클레 고유값) 를 도출한다. 따라서 w(t)=v(t)-v가 지수적으로 0으로 수렴함을 보이며, 이는 ‘경계층’ 동역학이 고정된 내부 상태에 대해 고유한 준정상상태(v)로 빠르게 수렴함을 의미한다.

이 지수 안정성을 전제로 특이 섭동 이론의 표준 절차를 적용한다. ε→0 한계에서 빠른 변수(v, u) 를 정적 관계식으로 대체하고, 남는 느린 변수 x_i에 대한 ODE만 남긴다. 구형 도메인의 그린함수 G(r,ℓ_j)=1/(4πD|r-ℓ_j|) (경계 조건을 포함한 정규화 형태)를 이용해 각 세포 간의 효과적인 상호작용 커플링 행렬 Γ를 명시적으로 계산한다. Γ_ij는 송신 세포 i가 방출한 신호가 수신 세포 j에 도달하는 정적 이득을 나타내며, 이는 전통적인 ‘통신 매트릭스’와 동일한 구조를 가진다.

결과적으로 원래의 무한 차원 PDE‑ODE 시스템은

dx_i/dt = f_i(x_i, Σ_j Γ_ij y_j)

형태의 유한 차원 다중 에이전트 네트워크로 축소된다. 여기서 y_j는 송신 세포의 LuxI 농도를 추출한 출력이며, 수신 세포는 y_j=0이다. 이 네트워크는 선형(Γ)·비선형(f_i) 결합 형태이므로 Lyapunov, 입력‑출력, 그래프 이론 등 기존 제어·시스템 분석 도구를 그대로 적용할 수 있다.

수치 실험에서는 3차원 구형 비드 안에 500~2000개의 세포를 무작위 배치하고, 완전 PDE‑ODE 시뮬레이션(유한 차분법)과 축소된 네트워크 모델을 비교한다. 시간 궤적, 평균 신호 농도, 그리고 동기화 정도에서 평균 오차가 2 % 이하이며, 계산 시간은 전자적 차원 감소 덕분에 1 % 수준으로 감소한다. 또한, 파라미터 스윕을 통해 Γ 행렬의 스펙트럼이 시스템 안정성에 미치는 영향을 정량화하고, 설계자가 원하는 전파 지연·감쇠 특성을 매트릭스 원소 조정으로 직접 설계할 수 있음을 시연한다.

본 연구의 주요 기여는 (1) 확산 매개 다세포 시스템의 빠른-느린 시간 스케일을 엄밀히 정량화하고, (2) 지수 안정성을 기반으로 특이 섭동을 적용해 PDE‑ODE 결합을 완전하게 유한 차원 네트워크로 변환한 점, (3) 구형 그린함수를 이용해 물리적 파라미터와 직접 연결되는 통신 행렬을 도출한 점, (4) 대규모 세포 집단에 대한 시뮬레이션 효율성을 크게 향상시킨 점이다. 이러한 접근은 합성 생물학, 조직공학, 그리고 환경 미생물 군집 설계 등에서 확산 기반 신호 전달을 정량적으로 설계·제어하려는 연구자들에게 강력한 도구가 될 것이다.