심방세동 파동역학의 에너지 지형 분석: 쌍쌍 최대엔트로피 모델 적용

초록

본 연구는 인간 심방의 2차원 전기 시뮬레이션 데이터를 대상으로 쌍쌍 최대엔트로피 모델(MEM)을 적용해 에너지 지형을 구축하였다. 정상 동리듬, 단일 안정 로터, 로터와 파동 파괴가 결합된 경우는 높은 정확도와 신뢰도(>0.9)로 모델링되었으나, 다중 파동렛 현상은 설명하지 못했다. 에너지 최소점(흡인점)과 장벽의 분포를 통해 각 파동 패턴의 안정성 및 전이 특성을 정량화했으며, ‘드라이버’(동방 혹은 로터)의 존재가 전반적인 파동 동역학을 좌우한다는 결론을 제시한다.

상세 분석

본 논문은 복잡계 물리학에서 차용한 쌍쌍 최대엔트로피 모델을 심장 전기 활동에 적용한 최초 사례 중 하나로, 몇 가지 중요한 기술적·이론적 통찰을 제공한다. 첫째, 저자는 512×512 격자의 2‑D 인간 심방 모델을 Courtemanche‑Ramirez‑Nattel 전기생리학적 방정식으로 구현하고, 네 가지 전형적인 파동 상태(정상 동리듬, 단일 안정 로터, 로터‑파동 파괴, 다중 파동렛)를 이온 전류 파라미터 변조를 통해 유도하였다. 여기서 각 셀의 전위 변화를 1 ms 간격으로 샘플링하고, 9개의 대표 전극 위치를 이진화하여 N=9 차원의 상태 공간을 정의하였다.

둘째, MEM의 핵심은 관측된 평균 활성도와 쌍별 상관을 보존하면서 엔트로피를 최대화하는 확률 분포 p(V) = e^{−E(V)}/Z 를 추정하는 것이다. 저자는 Ising 형태의 해밀토니안을 E = −∑i h_i σ_i − ½∑{i≠j} J_{ij} σ_i σ_j 로 설정하고, 반복 스케일링(α=0.1)으로 h_i와 J_{ij}를 학습하였다. 모델 적합도는 Kullback‑Leibler 발산 기반 정확도와 엔트로피 차이를 이용해 평가했으며, 정상 동리듬·단일 로터·로터‑파동 파괴에서는 정확도 0.90 이상, 신뢰도 0.99 이상을 기록했다. 반면 다중 파동렛에서는 정확도가 0.51에 머물러, 쌍쌍 상관만으로는 복잡한 다중 파동 상호작용을 포착하지 못함을 보여준다.

셋째, 학습된 해밀토니안을 바탕으로 전역 에너지 지형을 구성하고, 로컬 최소점(흡인점)과 그 주변의 베이시스 크기, 최소점 간 에너지 장벽을 정량화하였다. 정상 동리듬과 단일 안정 로터는 소수의 깊은 최소점과 넓은 베이시스를 가지며, 높은 장벽으로 인해 상태 전이가 드물어 ‘드라이버’가 시스템을 강하게 고정시킨다. 로터‑파동 파괴는 최소점 수가 증가하고 장벽이 낮아져, 파동 파괴와 재생이 빈번히 일어나며 불안정성을 나타낸다. 다중 파동렛은 수백 개에 달하는 얕은 최소점과 낮은 장벽을 보이며, 전이 경로가 촘촘히 얽혀 있어 전반적인 동역학이 고도로 혼돈적임을 시각화한다.

넷째, 저자는 이러한 결과를 뇌 과학에서 활용된 에너지 지형 개념과 직접 연결한다. 뇌의 휴지기 활동을 설명하는 MEM과 유사하게, 심장 전기 활동에서도 쌍쌍 상관이 주요 구조적·기능적 연결을 반영한다는 점을 강조한다. 다만, 현재 모델은 쌍쌍 상호작용에 국한되므로 고차 상관이나 비선형 전도성(예: 이방성, 섬유 구조) 등을 반영하지 못한다는 한계가 명시된다. 향후 고차 MEM, 비정상 분포, 실제 임상 전극 데이터 적용 등을 통해 모델의 일반화와 예측력을 향상시킬 여지가 있다.

마지막으로, 이 연구는 전산심장학과 통계물리학을 융합한 방법론적 프레임워크를 제시함으로써, 심방세동의 메커니즘을 정량적으로 해석하고, 드라이버(동방·로터) 표적 치료 전략을 설계하는 데 이론적 기반을 제공한다는 점에서 큰 의의를 가진다.