소뇌 퍼킨지 세포의 토러스 카나드 현상

초록

이 논문은 퍼킨지 세포의 상세 모델을 단순화한 5차원 방정식에서, 전류 J의 조절에 따라 폭발적 발화(버스팅)에서 고속 연속 스파이킹으로 전이되는 과정에서 ‘토러스 카나드’라 불리는 새로운 동역학 현상이 나타남을 보여준다. 느린 M‑전류와 빠른 칼슘 전류가 결합된 슬로우‑패스트 구조 속에서, 제한 주기(리밋 사이클)의 안정·불안정 분기점을 따라 일시적으로 불안정한 리밋 사이클을 따라가는 현상이 관찰된다. 이는 기존의 고정점 기반 카나드와는 달리 리밋 사이클 기반 토러스 카나드로, 실제 생물학적 뉴런 모델에서도 재현된다.

상세 분석

본 연구는 559개의 구획과 12종 이온 전류를 포함하는 대규모 퍼킨지 세포 모델을, 전압 V와 네 개의 개폐 변수(n, h, c, M)로 구성된 5차원 축소 모델로 변환하였다. 여기서 M‑전류는 다른 전류에 비해 10배 이상 느린 시간 상수를 가지며, 슬로우‑패스트 분해에서 슬로우 변수로 작용한다. 전류 J를 외부 자극으로 간주하고, J 값을 감소시키면서 시스템의 동역학을 AUTO와 XPPAUT을 이용해 연속적으로 추적하였다. J가 -23 mV 정도일 때는 전통적인 fold‑of‑fixed‑points와 fold‑of‑limit‑cycles가 결합된 ‘fold/fold‑cycle bursting’이 나타난다. J가 J*≈‑32.938 mV에 접근하면, 빠른 스파이크가 시작된 후 M‑전류가 서서히 증가하여 제한 주기의 fold 지점에 도달한다. 이때 기대되는 것은 스파이크가 멈추고 고정점으로 돌아가는 것이지만, 실제로는 시스템이 일시적으로 불안정한(limit cycle) 분기(branch)를 따라 이동한다. 즉, 전압이 높은 리밋 사이클의 불안정한 부분을 따라가며 M‑전류는 감소하고, 다시 안정적인 리밋 사이클 혹은 고정점으로 회귀한다. 이러한 과정이 반복되면서 스파이크 진폭이 주기적으로 변조되는 AM(Amplitude Modulated) 사이클이 형성된다. J가 더 낮아지면 AM 사이클이 사라지고, 순수한 고속 스파이킹이 나타난다. 중요한 점은 이 전이 구간에서 발생하는 카나드가 고정점이 아닌 리밋 사이클의 안정·불안정 분기점 사이를 오간다는 점이다. 이는 전통적인 2차원 카나드(예: van der Pol)와는 근본적으로 다른 메커니즘이며, 토러스 분기와 결합된 ‘토러스 카나드’라는 새로운 용어를 제안한다. 토러스 카나드는 토러스(2차원 주기 궤도) 위에서 발생하는 분기와, 그 토러스가 불안정 리밋 사이클을 따라 흐르는 구간이 존재함을 의미한다. 또한, 토러스 카나드가 존재하는 파라미터 구간에서는 첫 번째 리아프노프 지수가 음수인 초임계 토러스 분기가 관찰되며, 이는 안정적인 토러스와 불안정한 리밋 사이클이 충돌해 새로운 안정적인 주기 궤도가 생성되는 과정과 일치한다. 이러한 동역학은 M‑전류와 칼슘 전류의 상호작용에 의해 조절되며, M‑전류가 충분히 강하면 버스팅이 유지되고, 약해지면 AM 스파이킹 혹은 순수 스파이킹으로 전이한다. 실험적으로는 퍼킨지 세포의 전류 클램프나 약물 차단을 통해 M‑전류를 조절함으로써 토러스 카나드 현상을 검증할 수 있다.