부정맥

계산물리 물리 수치해석 마지막 수정: 2026-08-19 04:38:09

1. 개요[편집]

부정맥 Cardiac arrhythmia
기전 분류	자동능 이상 · 촉발 활동 · 재진입
핵심 양	파장 λ = 전도속도 × 불응기
불안정성	APD 회복곡선 기울기 > 1 → 교대맥
세포 모형	루오-루디 · 텐 투셰르 · 오하라-루디
조직 모형	단일영역(monodomain) 반응-확산
요구 해상도	격자 ~0.1 mm, 시간 ~0.005 ms

부정맥(arrhythmia)은 심장의 전기적 흥분이 동방결절이 지휘하는 정상 순서와 주기를 벗어난 상태를 통칭하며, 계산 심장학에서는 흥분성 매질 위의 반응-확산 방정식이 정상 표적파 대신 재진입 회로·나선·시공간 카오스로 넘어가는 분기 현상으로 다뤄진다. 임상 분류(심방세동, 심실빈맥, 방실 블록 …)는 겉으로 드러난 표현형이고, 모형 쪽에서 실제로 붙잡는 것은 그 아래의 동역학이다.

역할 분담을 먼저 정리해 두면 — 파면이 끊겨 감기는 2차원 구조(팁, 코어, 미앤더, 필라멘트)는 나선파, 불응기의 정의와 파장 공식의 유도는 불응기, 이온 전류의 미시 기전은 이온 채널, 충격으로 이것을 끄는 문제는 제세동에 있다. 이 문서는 부정맥이 성립하는 조건과, 그것을 세포에서 조직 규모까지 어떻게 계산하는가를 다룬다.

2. 세 가지 기전[편집]

교과서적 분류는 셋이다. 셋 다 모형에서 서로 다른 방식으로 구현된다.

1. 자동능 이상. 원래 자발 박동을 내지 않아야 할 심근이 스스로 발화하기 시작한다. 모형에서는 배경 전류나 과분극 활성화 전류( $I_f$ )의 크기를 건드리면 정지점이 불안정해지며 호프 분기로 극한 순환이 나온다 — 흥분성에서 진동성으로의 전이다.

2. 촉발 활동. 앞선 활동전위가 있어야만 나오는 이차 탈분극이 원인이다. 두 종류가 있다.

조기 후탈분극(EAD). 고평부가 비정상적으로 길어지면 L형 칼슘 전류가 불활성화에서 회복해 재활성화되고, 재분극 도중 전위가 다시 튀어 오른다. 재분극을 늦추는 요인(칼륨 전류 차단, 서맥, 저칼륨혈증)이 조건이며, 위상평면에서는 느린 다양체 위에서 안장-마디 분기가 사라졌다 생겼다 하는 구조로 읽힌다.
지연 후탈분극(DAD). 근소포체에 칼슘이 과부하되면 재분극이 끝난 뒤 자발적으로 칼슘이 방출되고, 나트륨-칼슘 교환체가 이를 안쪽 방향 전류로 바꿔 막을 밀어 올린다. 칼슘 순환을 명시적으로 들고 있는 모형에서만 재현되므로, 모형 선택이 곧 재현 가능한 현상 목록을 정한다는 것을 보여 주는 대표 사례다.

3. 재진입. 지속성 부정맥의 압도적 다수가 여기 속하고, 계산 쪽 관심도 거의 전부 여기 있다.

3. 재진입 — 파장이 회로보다 짧아야 한다[편집]

흥분파가 한 바퀴 돌아 출발점에 왔을 때 그 자리가 아직 불응이면 파는 죽고, 이미 회복했으면 다시 돈다. 그래서 결정적인 양이 파장이다.

\lambda = \theta \times \mathrm{ERP}

회로 길이가 $\lambda$ 보다 길어야 재진입이 유지된다. 유도와 임상적 함의(III군 약이 ERP를 늘려 $\lambda$ 를 키우는 것, I군 약이 $\theta$ 를 떨어뜨려 오히려 $\lambda$ 를 줄이는 역설)는 불응기에 정리돼 있다.

재진입이 시작되려면 조건이 하나 더 필요하다. 단방향 블록이다. 상대 불응기 한복판에 조기 자극이 들어가면 회복이 앞선 방향으로만 전도되고 아직 불응인 방향으로는 막히는데, 그 사이 뒤쪽이 회복되면 파가 되돌아 들어와 회로가 닫힌다. 조직이 이질적일수록(불응기가 자리마다 다를수록) 이 사고가 잘 난다. 정리하면 재진입의 성립 조건은 단방향 블록 + 느린 전도 + 회로 길이 > $\lambda$ 세 개다.

회로가 무엇으로 만들어지느냐에 따라 두 부류로 갈린다.

구분	회로를 정하는 것	흥분성 간극	페이싱으로 종료
해부학적 재진입	고정된 비흥분 장애물(흉터, 판막륜, 부전도로)	있음	대체로 가능
기능적 재진입	매질 자체의 동역학 (나선 코어)	거의 없음	어려움

흥분성 간극(excitable gap)은 파면의 머리와 앞선 파의 불응 꼬리 사이에 남은, 이미 회복된 구간이다. 이것이 있으면 외부에서 적절한 타이밍에 자극을 밀어 넣어 회로 안으로 파를 하나 더 집어넣을 수 있고, 그 파가 앞뒤로 충돌해 회로를 지운다. 항빈맥 페이싱(ATP)이 작동하는 물리적 이유가 정확히 이것이며, 간극이 사실상 없는 기능적 재진입에서는 통하지 않는다. 자세한 것은 제세동 참고.

기능적 재진입을 어떻게 볼 것인가는 이 분야의 오랜 논쟁거리였다. 1970년대의 선도 회로(leading circle) 개념은 “회로가 유지될 수 있는 가장 작은 크기로 저절로 줄어들고 그래서 간극이 없다”는 그림이었고, 이후 나선파 이론이 들어오면서 회전자(rotor) 관점 — 흥분하지 않는 코어를 파면이 감고 도는 위상 특이점 — 으로 대체됐다. 두 그림은 코어의 성격에서 결정적으로 다르다. 선도 회로의 중심은 계속 자극받지만 흥분하지 못하는 영역이고, 나선 코어는 곡률 때문에 파면이 도달하지 못하는 영역이다. 광학 매핑으로 위상 특이점이 화면 위를 기어다니는 것이 직접 관측되면서 회전자 쪽으로 무게가 실렸다.

4. 회복 곡선과 교대맥[편집]

부정맥이 어떻게 나빠지는가를 설명하는 표준 경로가 교대맥이다. 활동전위 지속시간이 직전 이완기 간격의 함수 $\mathrm{APD}_{n+1} = f(\mathrm{DI}_n)$ 라는 관계 아래서 고정 주기 자극은 1차원 사상의 반복이 되고, $\lvert f' \rvert > 1$ 에서 주기배가 분기가 일어나 길고 짧은 APD가 번갈아 나온다. 여기까지는 불응기에 있다.

조직으로 올라가면 새로운 것이 나온다.

동조 교대맥(concordant). 조직 전체가 같은 위상으로 길고-짧고를 반복한다. 그 자체로는 전도 블록을 만들지 않는다.
비동조 교대맥(discordant). 공간적으로 위상이 뒤집힌 영역이 생기고, 그 사이에 위상이 0인 마디선(nodal line)이 놓인다. 마디선 양쪽에서 APD 차이가 최대가 되므로 재분극 경사가 급격해지고, 짧은 쪽에서 온 파가 긴 쪽의 불응 꼬리에 걸려 국소 전도 블록이 난다. 파면이 끊기면 자유단이 생기고, 자유단은 나선이 된다.

비동조 교대맥이 생기는 데는 APD 회복곡선만으로는 부족하고 전도속도 회복곡선( $\theta$ 도 DI의 함수라는 것)과 전기긴장적 결합이 함께 작용한다. 실제로 회복 기울기 1이라는 기준은 필요조건도 충분조건도 아니라는 반례가 여럿 나와 있고, 대안 기전으로 칼슘 주도 교대맥이 있다. 근소포체의 칼슘 방출-재흡수 자체가 박동마다 교대하는 불안정성을 가지며, 이것이 나트륨-칼슘 교환체를 통해 전압에 되먹임되고 전압은 다시 칼슘 방출을 바꾼다. 전압과 칼슘이 양방향으로 결합된 2차원 사상이 되는 셈이라, 결합 부호에 따라 두 진동이 같은 위상으로 붙기도 하고 반대로 붙기도 한다.¹

5. 나선 붕괴와 세동[편집]

단일 안정 나선은 규칙적인 빈맥이다. 나선이 불안정해져 파면이 계속 새로 끊기면 파원 개수가 늘어나며 세동으로 간다. 이 그림의 원형은 놀랍게도 아주 이른 시기의 계산 실험이다 — 1964년 모(Moe) 그룹이 육각 격자 위 셀룰러 오토마타로 불응기를 자리마다 다르게 준 조직을 돌려, 파면이 계속 쪼개지며 여러 개의 잔물결(multiple wavelet)이 떠도는 상태를 재현했다. 당대 컴퓨터로 심방세동을 만들어 낸 이 모형이 다중 잔물결 가설의 근거가 됐고, 이후 수십 년의 표준 그림이었다.

현대적 관점은 여기에 두 가지를 얹는다. 하나는 잔물결이 무작위가 아니라 소수의 회전자가 만든 파면이 이질적 조직을 지나며 깨지는 것이라는 회전자 중심 해석이고, 다른 하나는 세동이 확률적 현상이 아니라 결정론적 반응-확산계의 시공간 카오스라는 것이다. 후자는 카오스 이론의 문법이 그대로 적용된다는 뜻이며, 실제로 세동 상태의 리아푸노프 지수나 상관 차원을 재는 연구가 있다.

조직 크기도 조건이다. 세동이 유지되려면 여러 개의 파가 동시에 존재할 만큼 조직이 커야 한다 — 임계 질량이다. $\lambda$ 가 조직 치수와 같은 자릿수인 작은 심장(예컨대 생쥐)에서는 세동이 잘 유지되지 않고 저절로 멈춘다. 사람 심실 정도라야 여러 파원이 공존한다. 계산 실험에서 도메인을 줄여 놓고 “이 모형은 세동을 안 만든다”고 결론 내리는 것이 흔한 실수인 이유다.

6. 단일세포 모형의 계보[편집]

조직 시뮬레이션의 반응항 $I_{\rm ion}$ 을 무엇으로 쓸 것인가가 절반이다. 계보는 대략 이렇게 흐른다.

노블(1962). 호지킨-헉슬리 모형을 푸르키녜 섬유에 맞춰 고친 최초의 심장 모형. 칼륨 전류를 두 개로 나눠 고평부를 만들어 냈다.
빌러-로이터(1977). 심실근 모형에 칼슘 전류를 명시적으로 넣었다. 상태변수 8개. 오랫동안 나선파 연구의 표준 반응항이었다.
루오-루디 I(1991)·II(1994). LR1은 기니피그 심실근 전류를 실측 데이터로 다시 맞춘 것이고, LR2가 결정적이다 — 세포 내 칼슘 농도와 근소포체 방출·재흡수를 동역학 변수로 들여와 DAD와 칼슘 교대맥을 모형이 스스로 만들어 낼 수 있게 했다. 이 문서 맥락에서 “루오-루디”라 하면 대개 이쪽을 가리킨다.
텐 투셰르-판필로프(2004, 2006). 사람 심실근. 상태변수 19개 규모로 조직 시뮬레이션에 쓸 만한 비용을 유지하면서 회복 곡선을 사람 데이터에 맞춘 것이 특징이라, 나선 붕괴 연구에서 가장 널리 쓰인다.
오하라-루디(2011). 사람 심실근, 상태변수 40여 개. 미제공 심장 100여 개의 실측을 기반으로 만들어 현재 사실상 표준이고, 약물의 심장 부작용을 다전류 차단으로 예측하는 규제용 파이프라인(CiPA)의 기반 모형이기도 하다.

반대 방향, 즉 일부러 줄인 모형도 한 축이다. 펜튼-카마(1998)는 전류 세 개·변수 세 개로, 미첼-셰퍼(2003)는 변수 두 개로 줄이되 회복 곡선을 목표 함수로 놓고 파라미터를 맞춘다. 조직 규모 파동 동역학은 대체로 회복 곡선과 파장이 결정하므로, 그 두 가지만 맞으면 상세 모형과 정성적으로 같은 나선·붕괴가 나온다. 비용 차이가 한두 자릿수라 대규모 파라미터 쓸기나 민감도 해석에는 이쪽이 정답이다.

선택 기준은 명확하다. 재현하려는 현상이 어느 변수에 사는지 보고 고른다. 나선 미앤더와 붕괴가 목적이면 축약 모형으로 충분하고, DAD·칼슘 교대맥·약물의 상태 의존 차단이 목적이면 그 변수를 실제로 들고 있는 상세 모형이 아니면 아예 재현이 불가능하다. 상세 모형을 쓴다고 답이 항상 나은 것도 아니다 — 파라미터가 40개 넘게 붙으면 그중 어느 것도 개별적으로 잘 식별되지 않는다.

7. 세포에서 조직으로[편집]

단일 세포 방정식에 확산을 붙이면 단일영역(monodomain) 모형이 된다.

\chi\left(C_m \frac{\partial V}{\partial t} + I_{\rm ion}(V,\mathbf{s})\right) = \nabla\cdot(\sigma\nabla V) + I_{\rm stim}, \qquad \frac{d\mathbf{s}}{dt} = \mathbf{f}(V,\mathbf{s})

$\chi$ 는 막의 표면적 대 부피비(심실근에서 1000~1500 cm $^{-1}$ 규모), $\sigma$ 는 유효 전도도 텐서다. 세포 안팎을 별도의 연속체로 분리하지 않고 하나로 뭉갠 것이며, 그래도 되는 이유와 그러면 안 되는 경우(충격 자극)는 제세동에서 이중영역 모형과 함께 다룬다.

실무에서 걸리는 지점들.

이방성. 심근 섬유 방향으로 전도가 2~3배 빠르고, 섬유 방향은 심실벽을 가로지르며 100° 넘게 회전한다. 즉 $\sigma$ 는 위치마다 주축이 도는 텐서이며, 이 회전이 3차원 필라멘트의 거동을 실제로 바꾼다. 등방 가정으로 얻은 결론을 심장에 옮기면 안 되는 이유.
해상도. 상승기의 파면 두께가 0.5 mm 수준이라 격자가 그보다 성기면 전도속도가 틀린다. 11개 코드가 같은 문제를 푼 국제 벤치마크(2011)에서, 성긴 격자(0.5 mm)에서는 코드 간 활성화 시각이 크게 갈렸고 0.1 mm 격자와 0.005 ms 스텝쯤 돼야 서로 수렴했다. 사람 심장 전체를 이 해상도로 채우면 격자점이 수천만~수억 개가 되며, 심장 시뮬레이션이 HPC 문제인 이유가 전부 여기 있다.²
전도속도를 맞추려고 $\sigma$ 를 손대는 유혹. 성긴 격자에서 느려진 전도속도를 $\sigma$ 를 키워 보정하는 것은 흔한 관행인데, 이렇게 맞춘 $\sigma$ 는 그 격자에만 유효하다. 격자를 세밀화하면 전도속도가 과대평가되고, 파장이 달라지므로 부정맥 유발 여부라는 결론 자체가 바뀐다. 메시 독립성 연구를 파형이 아니라 전도속도와 파장으로 판정해야 하는 이유다.
시간 적분. 반응항이 뻣뻣하고 확산항은 부드러우므로 연산자 분리가 표준이다. 확산은 암시적으로, 반응은 격자점별 독립 상미분방정식으로 푼다. 게이트 변수는 러시-라슨(지수 적분기의 심장판)으로 준해석적으로 밀면 강성 방정식의 부담이 크게 준다. 반응항이 격자점마다 완전히 독립이라 GPU 컴퓨팅과 궁합이 극단적으로 좋다는 것도 이 구조의 선물이다.
자극 프로토콜. 부정맥 유발 실험은 S1–S2 혹은 S1–S2–S3 조기 자극을 결합 간격을 쓸어가며 반복하는 것이라, 한 번의 “실험”이 수십~수백 회의 전체 시뮬레이션이다. 계산량이 폭발하는 지점이 방정식이 아니라 프로토콜이라는 것이 이 분야의 특징이다.
검증. 활성화 시각 지도, 전도속도, APD 지도, 심전도 모사 파형이 정량 비교 지표다. “나선이 예쁘게 나왔다”는 검증이 아니라는 것은 검증 및 확인의 일반 원칙 그대로다.³

8. 관련 문서[편집]

그래서 “회복 곡선 기울기를 1 아래로 누르면 세동을 막는다”는 한때의 약물 설계 지침은 지금은 반쪽짜리 취급을 받는다. 기울기를 눌렀는데 칼슘 쪽 불안정성이 남아 있으면 교대맥이 그대로 나오고, 반대로 기울기가 1을 넘는데도 전기긴장적 결합이 충분히 강하면 공간적으로 안 갈라진다. 저차원 지표 하나로 고차원 계의 안정성을 예측하려는 시도가 얼마나 자주 실패하는지에 대한 좋은 표본. ↩
이 벤치마크가 유명해진 진짜 이유는 “각 코드가 자기 격자에서 다 수렴했다고 주장했는데 서로 답이 달랐다”는 데 있다. 코드 내부의 수렴 검사는 그 코드의 이산화가 옳다는 증거가 되지 못한다. 서로 다른 구현이 같은 답에 모이는지를 봐야 한다는, 심장 분야에 국한되지 않는 교훈이다. ↩
심장 모형에는 자연스러운 검증 데이터가 있다는 점이 부럽기도 하다. 실제 심전도라는 것이 100년 넘게 축적된 비침습 관측이니까. 문제는 그 신호가 심장 전기활동을 몸통이라는 저역 필터에 통과시킨 결과라, 역방향으로 되돌리는 것이 전형적인 부적절 역문제라는 점이다. 정방향 모사가 훨씬 쉽다. ↩

부정맥