Soliton-like Waves in a Two-Dimensional Recurrent Spiking Neural Network with Weighted Spike-Timing-Dependent Plasticity

기술 요약: 가중치 기반 스파이크 타이밍 의존 가소성(WSTDP)을 갖춘 2차원 재귀적 스파이킹 신경망에서의 솔리톤 유사 파동

문제 제기
신경 활동의 이동 파동(traveling waves)은 피질에서 도처에 존재하지만, 재귀적 회로 내에서 이러한 파동이 어떻게 발생하고 스스로 유지되는지에 대한 메커니즘은 여전히 완전히 이해되지 않은 상태로 남아 있다. 신경장(neural field) 모델들은 "멕시칸 모자(Mexican-hat)" 형태의 연결 프로파일(단거리 흥분과 장거리 억제의 균형)을 바탕으로 이동 파동과 고립된 범프(bumps)를 예측해 왔으나, 이러한 모델들은 일반적으로 개별 스파이크와 시냅스 가소성을 추상화하여 생략한다. 반대로, 스파이크 타이밍 의존 가소성(STDP)이 수용장(receptive fields)을 형성하고 경로를 강화하는 것으로 알려져 있음에도 불구하고, STDP를 갖춘 재귀적 흥분/억제(E/I) 네트워크가 균일한 상태로부터 안정적이고 자가 전파되는 파동 패킷(솔리톤 유사 구조)을 자율적으로 생성할 수 있는지, 그리고 이러한 구조가 소멸(annihilation)과 같은 특정한 충돌 역학을 보이는지는 불분명하다.

방법론
저자들은 200 × 100 쌍의 흥분성(E) 및 억제성(I) 뉴런이 이산 격자(discrete grid) 상에 배치된 2차원 재귀적 네트워크를 구성하는 최소한의 생물학적 타당성을 갖춘 계산 모델을 구축하였다. 모델은 이산 시간 단계(discrete time steps)에서 작동하며 동기식 업데이트를 수행한다.

• 뉴런 모델: 각 뉴런은 전시냅스 이웃들로부터 가중치가 부여된 스파이크를 통합한다. 핵심 특징은 시냅스 입력의 **분할 정규화(divisive normalization)**로, 총 흥분성 구동력은 흥분성 시냅스 가중치의 합($S_E$)에 의해 정규화된다. 이 정규화는 이진 접합부(binary junctions)에서의 국소적 감쇄를 이용한 생물학적으로 타당한 수상돌기 구현을 제안한다. 발화 임계값($\theta$)은 항상성 적응을 거친다: "정지 상태(stall states)"를 방지하기 위해 침묵 중에는 최소값($\theta_{min}$)을 향해 감소하며, 스파이크 발생 시에는 스파이크 빈도 적응을 구현하기 위해 최대값($\theta_{max}$)을 향해 증가한다.
• 연결성: 네트워크는 기하학적 비대칭성을 특징으로 한다: 흥분-흥분($E \to E$) 연결은 흥분-억제($E \to I$) 및 억제-흥분($I \to E$) 연결보다 더 큰 반경($r_{EE}=9$)을 갖는다 ($r_{EI}=r_{IE}=5$). $I \to I$ 연결은 존재하지 않는다. 초기 가중치는 균일하며, 억제성 입력의 수가 훨씬 적음을 보상하기 위해 억제성 가중치를 흥분성보다 훨씬 높게 설정하였다.
• 학습 규칙: 시냅스 가중치는 곱셈적이고 시간적 비대칭성을 가진 **가중치 기반 STDP(WSTDP)**를 통해 진화한다. 이 규칙은 가중치가 하드 클리핑(hard clipping) 없이 경계 내에 유지되도록 보장하며, 장기 증강(LTP)과 장기 억제(LTD) 이벤트 비율에 의해 결정되는 고정점으로 수렴한다.
• 자극: 네트워크는 고정된 주기로 반복적인 국소 자극(흥분성 뉴런의 디스크 영역이 발화하도록 강제함)을 받는다. 실험은 단일 소스 자극과 다양한 위상 및 주파수를 가진 이중 소스 자극을 포함하여 파동 상호작용을 연구하기 위해 수행되었다.

주요 기여 및 결과

1. 소산 솔리톤(Dissipative Solitons)의 출현: 공간적으로 균일한 연결성으로부터 시작하여, 반복적인 국소 자극은 안정적이고 자가 전파되는 파동 패킷의 자발적 형성을 유도한다. 이 패킷들은 일정한 공간적 프로파일(밝은 흥분성 링과 그 뒤를 따르는 억제성 궤적)을 유지하며 일정한 속도로 전파된다. 결정적으로, 일단 형성되면 이 파동들은 외부 자극이 제거된 후에도 자율적으로 계속 전파되는데, 이는 이들이 자가 유지 구조임을 입증한다.
2. 기하학적 비대칭성과 WSTDP의 역할: 이 파동의 출현에는 두 가지 조건이 필요하다:
   • 기하학적 비대칭성: 기능적인 멕시칸 모자 프로파일을 만들기 위해 흥분성 연결 반경이 억제성 반경보다 커야 한다.
   • WSTDP: 학습 규칙은 전파 방향을 시냅스 가중치 지형에 "각인"한다. 전파 방향의 시냅스는 순 LTP를 겪는 반면, 역방향의 시냅스는 순 LTD를 겪는다. 이는 역방향 전파를 방지하고 파동을 안정화하는 방향성 편향을 생성한다.
3. 충돌 역학 및 소멸: 두 파동 패킷이 정면으로 충돌할 때, 이들은 서로 통과하는 대신 소멸한다. 이러한 행동은 소산 솔리톤의 특징이다. 소멸은 한 파동의 억제성 궤적이 다른 파동의 흥분성 전면을 억제함으로써 발생한다.
4. 반-지속적 공간 경계(Semi-Persistent Spatial Boundaries): 이중 소스 실험에서 충돌 지점은 시냅스 가중치에 인코딩된 "반-지속적" 경계를 형성한다.
   • 위상 인코딩: 소스들이 동위상(in phase)이면 경계는 기하학적 중간 지점에 형성된다. 위상 차이가 있으면 경계는 뒤처지는 소스 쪽으로 이동한다.
   • 주파수 경쟁: 소스들이 서로 다른 주파수를 가지면, 빠른 소스가 느린 소스의 시냅스 가중치를 점진적으로 덮어쓰며 경계가 표류(drift)하게 된다. 충분히 큰 주파수 우위는 완전한 영토 탈취로 이어진다.
   • 견고성: 이러한 경계는 한 소스가 비활성화된 후에도 지속되지만, 그 안정성은 적응형 임계값 매개변수에 의해 조절된다.
5. 강한 초기 억제의 필요성: 모델은 강한 초기 억제 가중치가 기하학적 필수 요소임을 보여준다. 이 가중치를 줄이면 국소적 솔리톤 대신 전역적 동기화 진동(간질 유사 상태)이 발생하며, 이는 국소적 파동 전면을 가두는 데 있어 억제의 결정적 역할을 강조한다.

의의 및 주장
본 논문은 단순한 국소 가소성 규칙으로부터 뇌 피질의 이동 파동, 활동 영역 경계 설정, 그리고 공간적 기억의 출현을 연구하기 위한 최소한의 계산 프레임워크를 제공한다고 주장한다.

• 출현 메커니즘: 본 연구는 (WSTDP 및 항상성 임계값 적응과 같은) 단순한 국로 규칙으로부터 안정적인 이동 파동 및 충돌 소멸과 같은 복잡한 시공간 역학이 (사전 설계된 공간 구조 없이도) 어떻게 나타날 수 있는지를 시사한다.
• 공간 기억: 저자들은 이러한 네트워크의 시냅스 가중치 지형이 일종의 공간 기억을 구현한다고 제안한다. 네트워크는 과거 입력의 상대적 타이밍(위상)과 비율(주파수)을 소멸 경계의 위치와 안정성을 통해 "기억"한다.
• 방향성 안정성: WSTDP 메커니즘은 전파를 가능하게 할 뿐만 아니라, 불안정한 역방향 모드를 억제함으로써 전파를 능동적으로 안정화하여 폭주하는 진동을 방지한다.
• 생물학적 타당성: 이 모델은 국소적인 스파이크 기반 가소성 틀 안에서 피질 회로가 어떻게 별개의 활동 영역을 분리하고 유지하며, 이동 파동이 경쟁적인 입력에 의해 어떻게 조절될 수 있는지에 대한 잠재적 설명을 제공한다.

저자들은 모델이 최소한의 형태(이산 시간, 이진 스파이크)이지만, 솔리톤 출현의 견고함은 이러한 현상이 특정 매개변수 선택의 산물이 아니라 재귀적 E/I 네트워크의 근본적인 특성임을 시사한다고 언급한다. 향후 연구에서는 비동기적 전파 및 구조화된 시냅스 구심성(afferents)에 대한 탐구를 제안한다.