생체계의 동적 전이와 정적 무질서의 큰 변동성

초록

이 논문은 흥분성 세포와 수동성 세포가 혼합된 격자 모델에서, 세포 간 결합 강도가 증가할 때 자발적인 진동이 어떻게 나타나는지를 연구한다. 정적 무질서(수동 세포의 공간적 분포)가 복제마다 다른 전이 임계점을 만들며, 시스템 크기가 커질수록 이러한 변동이 로그 스케일로 확대된다. 결과는 큰 기관일수록 지속적 활동이 시작되는 시점에 더 큰 변동성을 보일 수 있음을 시사한다.

상세 분석

논문은 먼저 흥분성 세포와 수동성 세포를 각각 이산적인 미분 방정식으로 기술하고, 흥분성 세포는 FitzHugh‑Nagumo 형태, 수동성 세포는 선형 회복 모델을 사용한다는 점을 명시한다. 각 흥분성 세포는 네 개의 최근접 이웃과 전기적 커플링을 갖으며, 동시에 주변에 존재하는 수동성 세포와도 연결된다. 수동성 세포의 개수는 격자 각 점마다 확률적으로 할당되며, 이는 ‘quenched disorder’라 불리는 고정된 무작위 분포를 만든다. 저자는 이러한 무질서가 복제마다 다른 국소 평균 수동 세포 밀도를 초래함으로써, 전이 임계점인 결합 강도 Gc가 변동한다는 사실을 실험적으로 확인한다.

시뮬레이션 결과는 두 가지 뚜렷한 동역학적 상태를 보여준다. 첫 번째는 모든 세포가 정지 상태에 머무는 ‘quiescent’ 단계이며, 두 번째는 전반적인 동기화 진동이 지속되는 ‘oscillatory’ 단계이다. 전이 구간에서는 일부 복제가 진동을 보이는 반면, 다른 복제는 여전히 정지 상태에 머무르는 현상이 관찰된다. 이러한 현상은 국소적인 수동 세포 밀도가 평균보다 높을 경우, 해당 영역에서 흥분성 세포가 더 쉽게 탈동기화되고, 결국 전체 네트워크가 진동을 시작하게 만든다.

저자는 이 변동성을 정량화하기 위해 ‘effective local density’ 라는 변수를 도입하고, 이를 평균 수동 세포 밀도 ρ와 시스템 크기 N에 대한 로그 함수와 결합한 스케일링 관계를 제시한다. 구체적으로, 전이 임계점 Gc(N)≈G0+α·log(N)/ρ 형태로 근사할 수 있음을 보이며, 여기서 G0와 α는 모델 파라미터에 따라 결정되는 상수이다. 이 식은 큰 시스템일수록 작은 변화에도 전이 임계점이 크게 이동한다는 것을 의미한다.

또한, 저자는 이러한 스케일링이 실제 생리학적 현상, 예를 들어 임신 후기 자궁의 수축 시작 시점과 연관될 수 있음을 논의한다. 자궁은 수천에서 수만 개의 근육 세포와 다양한 비근육성 세포가 얽혀 있는 큰 조직이며, 수동성 세포(예: 섬유아세포)의 분포가 불균일할 경우, 전반적인 수축 활동이 시작되는 시점에 큰 개체 간 변동성을 보일 수 있다. 따라서 이 모델은 조직 규모와 세포 구성의 무질서가 기능적 전이 현상에 미치는 영향을 정량적으로 설명하는 틀을 제공한다.

결론적으로, 논문은 정적 무질서가 동적 전이의 임계점에 로그 스케일로 영향을 미치며, 이는 큰 생물학적 시스템에서 변동성을 증폭시킨다는 중요한 통찰을 제공한다. 이러한 결과는 신경망, 심장 조직, 그리고 자궁과 같은 다양한 흥분성-수동성 혼합 시스템에 일반화될 가능성을 시사한다.