사회아메바 군집 형성의 몬테카를로 시뮬레이션 모델

초록

본 논문은 사회아메바(Dictyostelium discoideum)의 군집 형성 과정을 이산형 몬테카를로 모델로 구현한다. 세포 간 cAMP 파동 전파, 초기 단계의 무작위 방출, 신호 전달 지연, 그리고 관성과 목적성을 갖는 무작위 보행을 모두 포함한다. 시뮬레이션 결과는 군집 시간의 평균과 분산이 일정한 평형값에 수렴함을 보여주며, 실험에서 관찰된 동시 방출, 파동 전파, 화학주성 등의 현상을 재현한다.

상세 분석

이 연구는 기존 연속 미분 방정식 기반 모델이 다루기 어려운 ‘임펄시브 방출’과 ‘파동 전파 지연’ 같은 비선형 현상을 이산형 셀룰러 오토마톤 형태로 구현한다는 점에서 혁신적이다. 신호 전파는 격자 상에서 15칸·초당(≈300 µm/min) 이동하는 파동으로 모델링되고, 거리와 시간에 따라 감쇠되는 식 s₍i,j₎ = k₁·rⁿ + w·dᵢⱼ을 사용한다. 여기서 r = 0.95는 PDE에 의한 감소, w = 0.038은 거리 가중치이며, n > 360이면 신호는 무시한다.

세포의 자율 방출은 초기 15~30분 간격에서 시작해, 나이(시뮬레이션 단계)와 함께 감소하는 함수 twl(n) = ⌈320·(k₂ − arctan(k₃·(n − k₄)/π)⌉ 로 구현된다. 기본 불응기 6분(120 단계) 이후 추가 불응기가 적용돼, 실제 실험에서 보고된 6분 주기의 동시 방출을 재현한다.

신호 전달(리레이)은 두 개의 임계값 thd₁ = 1.1과 thd₂ = 2.5×10⁻⁵을 도입해, 세포가 주변 신호의 절대값과 변화량을 동시에 만족할 때만 cAMP 펄스를 방출한다. 방출 전 5단계(15 s)의 지연을 두어 실제 세포의 반응 지연을 반영한다.

세포 이동 메커니즘은 cAMP 농도에 따라 두 가지 모드로 전환된다. thd₃ = 2.42 이하이면 무작위 보행을 수행하고, ½ 확률로 정지, 1/16 확률로 8방향 중 하나로 이동한다. 농도가 thd₃를 초과하면 ‘목적성 무작위 보행’으로 전환되며, 가장 강한 두 신호 방향 d₁, d₂의 내적을 이용해 이동 방향 d를 결정한다. 강도에 따라 m₁=6, m₂=4, m₃=2 단계 동안 해당 방향을 유지하고, 매 단계 주변 신호를 재평가해 방향을 수정하거나 감소시킨다. 이는 ‘관성’과 ‘스티어링’ 효과를 이산형 모델에 자연스럽게 삽입한 것이다.

시뮬레이션은 400×400 격자(40 µm×40 µm)에서 180개의 세포를 3 s 간격으로 진행한다. 결과는 군집 시간 평균과 분산이 일정한 평형값에 수렴함을 보여, 초기 무작위 방출이 충분히 누적될 경우 동시 방출과 파동 전파가 급격히 일어나 군집이 형성된다는 실험적 사실을 재현한다. 또한 파라미터 변동에 대한 민감도 분석을 통해, thd₁과 thd₃의 상대적 크기가 동시 방출 시점과 군집 속도에 결정적인 영향을 미침을 확인한다.

이 모델은 연속 모델이 놓치기 쉬운 ‘펄스형 신호’, ‘파동 전파 지연’, ‘불응기 가변성’ 등을 명시적으로 다루면서도, Monte Carlo 방식으로 대규모 시뮬레이션이 가능하도록 설계돼, 향후 다른 화학주성 시스템(예: 면역세포, 암 전이)에도 확장 가능성을 제시한다.