비접착성 기판 위 활성 세포의 집합 역학

초록

본 연구는 세포 간 결합 친화성을 기반으로 한 셀룰러 오토마타 모델을 제시하여, 외부 신호 없이 무작위로 배치된 세포들이 어떻게 비접착성 표면에서 응집하고 컴팩트한 조직 구조를 형성하는지를 규명한다. 확산과 로밍(구조 재배열) 두 가지 과정을 통해 두 개의 시간 스케일이 나타나며, 군집 성장 속도, 면적의 비단조적 변화, 그리고 프랙탈 차원을 통한 거칠기 분석 결과가 실험과 일치함을 보여준다.

상세 분석

이 논문은 활성 세포 집합 현상을 설명하기 위해 2차원 정사각 격자(L × L) 위에 N₀개의 세포를 무작위로 배치하고, 각 세포를 하나의 클러스터로 간주하는 셀룰러 오토마타(CA) 프레임워크를 구축한다. 핵심 규칙은 두 가지 확률적 행동으로 나뉜다. 첫째, 클러스터 크기 n에 반비례하는 이동 확률 Pₘ(n)=μ/n(μ=1)으로, 클러스터가 주변 빈 칸으로 확산한다. 클러스터가 다른 클러스터와 접촉하면 즉시 결합(sticking probability=1)하여 새로운 클러스터를 형성한다. 둘째, 로밍 확률 Pᵣ(n)=β n(β는 세포 유형에 따라 조정)으로, 클러스터 내부에서 개별 세포가 주변 이웃 수를 최대화하도록 재배열한다. 재배열은 에너지 함수 E=−γ₁n₁−γ₂n₂−γ₃n₃(γ₁=3,γ₂=2,γ₃=1)로 정의된 결합 에너지 감소를 목표로 하며, Metropolis 기준(r<exp(−ΔE))을 적용해 확률적으로 허용한다.

시뮬레이션 결과는 순수 확산만을 고려했을 때 확산 제한 응집(DLA)과 유사한 가지형 구조가 형성되지만, 로밍 과정이 추가되면 시간이 지남에 따라 클러스터가 내부로 수축하고 최종적으로 원형에 가까운 고밀도 조직으로 변한다. 중요한 두 시간 스케일이 관찰된다. 첫 번째는 작은 클러스터가 빠르게 확산·충돌하여 집합하는 ‘확산 단계’이며, 두 번째는 클러스터가 임계 크기 n_c(실험값에 따라 10~30) 이상이 되면 이동이 거의 정지하고 로밍에 의해 내부 재배열이 진행되는 ‘압축 단계’이다.

동역학적 지표로는 (i) 클러스터 수 N(t)의 감소율, (ii) 전체 집합 면적 A(t)의 비단조적 변화(초기 증가 후 감소), (iii) 군집의 컴팩트성(평균 인접 이웃 수)와 결합 강도(시간에 따른 평균 E), (iv) 프랙탈 차원 D_f를 이용한 거칠기 분석이 있다. 특히 A(t)의 비단조성은 실험 보고와 일치하며, 로밍 파라미터 β가 클수록 더 빠른 컴팩트화와 낮은 D_f(≈1.2)로 나타나 정상 조직과 종양 조직의 구조적 차이를 정량화할 수 있음을 시사한다.

모델의 장점은 복잡한 화학·기계적 신호를 명시적으로 도입하지 않고도 세포 간 결합 친화성이라는 최소 조건만으로도 관찰된 현상을 재현한다는 점이다. 이는 Steinberg의 차등 부착 가설(Differential Adhesion Hypothesis)을 미시 수준에서 구현한 것으로, 세포 자체의 ‘결합 강도 최대화’ 원리가 조직 형성의 주요 구동력임을 뒷받침한다. 또한, 파라미터 μ와 β를 조절함으로써 다양한 세포 유형(예: 정상 vs. 암세포)의 집합 속도와 최종 형태를 예측할 수 있어, 조직 공학 및 암 연구에 실용적인 모델링 도구로 활용될 가능성이 크다.