뇌전증 발작 중 기능적 네트워크와 동기화 가능성의 동적 변화

초록

본 연구는 60명의 환자에서 기록된 100개의 국소 발작 전·중·후 EEG 데이터를 이용해, 시간에 따라 변화하는 기능적 뇌 네트워크의 구조적 지표(클러스터 계수, 평균 최단 경로)와 라플라시안 스펙트럼 기반 동기화 가능성(특징값 비율 S)을 분석하였다. 발작 진행 동안 네트워크는 무작위성에서 규칙성으로, 다시 무작위성으로 이동하는 ‘볼록형’ 변화를 보였으며, 동기화 가능성은 발작 중 가장 낮았다가 발작 종료 직전에 급격히 회복되는 패턴을 나타냈다. 이러한 결과는 발작 종료가 네트워크 재구성과 동기화 회복에 의해 자가조절되는 메커니즘일 수 있음을 시사한다.

상세 분석

본 논문은 복합 네트워크 이론을 인간 뇌 전기생리학에 적용한 선구적 연구로, 60명의 약물 저항성 뇌전증 환자에게서 수집된 100개의 국소 발작 EEG 데이터를 2.5 초 슬라이딩 윈도우로 나누어 기능적 연결 행렬 R을 구축하였다. 연결 강도는 정규화된 최대 지연 상관계수 ρmax ij 로 정의했으며, 각 윈도우마다 그래프가 연결성을 유지하도록 임계값 T(w)를 조정해 이진 인접 행렬 A(w)를 얻었다. 이렇게 얻은 시계열 네트워크에 대해 평균 최단 경로 L, 클러스터 계수 C, 에지 밀도 ε를 계산하고, 무작위 그래프(동일 차수 분포)와의 비율 C/Cr, L/Lr을 통해 구조적 변화를 정량화하였다.

동기화 가능성 평가는 라플라시안 L의 스펙트럼을 이용했으며, 최소 비영(λmin)과 최대 고유값(λmax) 사이의 비율 S = λmax/λmin을 사용하였다. S가 클수록 네트워크는 전역 동기화에 취약하고, 작을수록 동기화가 용이함을 의미한다. 결과는 세 가지 주요 특징을 보였다. 첫째, 발작 전후에 C와 L은 무작위 네트워크에 비해 약간 높은 값을 보이며, 발작 중반에 가장 큰 편차를 나타내어 네트워크가 보다 규칙적인(소월드) 형태로 전환됨을 보여준다. 둘째, ε는 발작 초반에 약 0.3 수준을 유지하다가 후반부에 서서히 증가하여, 연결 밀도가 발작 진행과 함께 강화됨을 시사한다. 셋째, S는 발작 초기에 감소했다가 중반에 최고점(동기화 저하)으로 상승하고, 발작 종료 직전 급격히 감소한다. 특히 S의 변동은 λmin의 동적 변화에 크게 좌우되며, λmin이 감소하면 네트워크는 일시적으로 ‘분절화’되어 동기화가 억제된다. λmax는 ε와 거의 일치하는 양상을 보여, 에지 밀도가 큰 고유값에 제한을 둔다는 이론적 기대와 일치한다.

이러한 구조‑동역학 관계는 발작이 시작될 때 뇌가 일시적으로 높은 동기화 상태를 피하기 위해 네트워크를 재구성하고, 발작 말기에 다시 동기화 가능성을 회복함으로써 자체 종료 메커니즘을 구현한다는 가설을 뒷받침한다. 또한, 클러스터 계수와 평균 최단 경로의 ‘볼록형’ 변화는 발작 전후에 무작위‑규칙‑무작위 전이 패턴을 일관되게 보이며, 이는 발작 부위의 해부학적 위치와 무관하게 보편적인 현상임을 확인한다.