다중 가중치 보로노이로 보는 2차원 세포 조직 역학

초록

본 논문은 세포 크기가 서로 다른 경우에도 적용 가능한 가중치 거리 함수를 이용해 일반화된 보로노이 영역을 정의하고, 이 영역을 기반으로 세포 간 접착·반발 힘, 조직 가장자리의 돌출 운동력, 그리고 확률적 교란을 포함한 동역학 모델을 구축한다. 스토캐스틱 미분 방정식 시뮬레이션을 통해 실제 조직에서 관찰되는 구형·구면형 경계와 다양한 토폴로지를 재현함과 동시에, 단일 세포쌍의 경우 비선형 Ornstein‑Uhlenbeck 과정을 해석적으로 조사한다. 결과는 세포 크기 이질성과 람엘라(세포 돌출)의 확장 사이에 존재하는 트레이드오프를 제시한다.

상세 분석

이 연구는 기존 보로노이 모델이 평면 접촉면만을 가정한 한계를 극복하기 위해, 각 세포를 반경 r_i 를 갖는 원형 ‘세포체’로 보고, 거리 함수 d_i(x)=‖x−c_i‖/r_i 와 같은 곱셈 가중치를 도입한다. 이렇게 정의된 가중 거리로부터 유도된 일반화된 보로노이 영역 V_i={x|d_i(x)≤d_j(x) ∀j≠i}는 세포 간 겹침을 방지하면서도 경계가 구면형(또는 구면 조각)으로 형성되도록 한다. 수학적으로는 각 경계면 Σ_ij가 두 세포 중심 c_i, c_j 사이의 라플라시안 곡면으로 표현되며, 이는 곡률이 세포 반경 비율에 따라 조정된다.

세포 간 상호작용은 두 종류의 힘 밀도로 모델링된다. 첫째, 접착력은 경계면에 걸쳐 일정한 양의 음의 힘(σ_a)으로 구현되어, 세포가 일정 거리 이하로 접근하면 에너지 최소화를 유도한다. 둘째, 반발력은 짧은 거리에서 급격히 증가하는 양의 힘(σ_r)으로, 세포가 과도히 겹치지 않도록 제어한다. 이 두 힘은 각각 접촉면의 단위 법선 방향으로 작용하며, 전체 힘은 면적 적분을 통해 각 세포 중심에 대한 총 힘 F_i로 집계된다.

조직 가장자리에서는 세포가 외부 매트릭스와 상호작용한다는 가정 하에, 돌출(lamellipodia) 운동력을 t_i·n_i 형태로 적용한다. 여기서 t_i는 세포가 스스로 생성하는 추진력의 크기, n_i는 외부 법선 벡터이다. 이는 조직 성장·재생 상황을 재현하는 데 필수적이다.

동역학은 뉴턴식 방정식 m_i·ẍ_i = F_i + ξ_i(t) 로 기술되며, ξ_i(t) 는 백색 가우시안 잡음으로, 세포 이동의 비결정성을 반영한다. 따라서 시스템은 다차원 스토캐스틱 미분 방정식(SDE) 형태가 되며, 수치적 통합에는 Euler‑Maruyama 방법이 사용된다.

특히, 두 세포가 짧게 접촉하는 경우를 분석하여 비선형 Ornstein‑Uhlenbeck (OU) 과정을 도출한다. 여기서 거리 변수 d(t) 는 평균 복원력 k·(d_eq−d) 와 잡음 σ·η(t) 로 구성된 SDE: ẋ = −k·(x−x_eq)+σ·η(t) 로 근사된다. 이 해석은 접촉 지속 시간, 탈착 확률 등을 정량화하는 데 활용된다.

시뮬레이션 결과는 (1) 세포 크기 이질성이 클수록 경계가 더 복잡하고, (2) 람엘라 길이가 길어질수록 조직 전체의 연속성이 향상되지만, 과도한 이질성은 구조적 불안정을 초래한다는 두 가지 주요 현상을 보여준다. 이는 자연 조직이 세포 크기와 돌출 구조 사이에서 최적의 균형을 맞추려는 진화적 압력을 반영한다는 가설을 뒷받침한다.