호지킨-헉슬리 모형

계산물리 물리 수치해석 마지막 수정: 2026-08-06 04:12:41

1. 개요[편집]

호지킨-헉슬리 모형 Hodgkin–Huxley model
형태	4변수 비선형 ODE ($V, m, h, n$)
발표	1952년, J. Physiol. 5부작
실험 대상	대서양 오징어 거대축삭 (지름 0.5~1 mm)
수상	1963년 노벨생리의학상 (에클스와 공동)
수치적 성격	강성 ODE · 제2형 흥분성
축약형	피츠휴-나구모 모형

호지킨-헉슬리 모형(Hodgkin–Huxley model, HH 모형)은 신경세포 막을 전기회로로 놓고, 이온 채널의 열림 정도를 나타내는 게이팅 변수들과 막전위를 함께 푸는 4변수 비선형 상미분방정식계로 활동전위를 재현하는 모형이다. 앨런 호지킨과 앤드루 헉슬리가 1952년 오징어 거대축삭의 전압 고정(voltage clamp) 실험 데이터를 맞춰 만들었고, 이 공로로 1963년 노벨생리의학상을 받았다.¹

이 모형이 특별한 이유는 생물학에서 처음으로 “정량적으로 맞는” 미분방정식 모형이었다는 데 있다. 채널 단백질의 존재조차 확인되기 20년 전에, 호지킨과 헉슬리는 전류 데이터만 보고 “나트륨 전도도는 어떤 양의 세제곱에 비례하고, 칼륨은 네제곱에 비례한다”는 지수를 골라냈다. 그리고 그 지수는 나중에 밝혀진 채널의 서브유닛 구조와 놀랍도록 잘 맞았다. 계산 자체는 손으로 돌린 기계식 계산기로 했고, 한 번의 활동전위를 그리는 데 며칠이 걸렸다고 한다.²

2. 방정식[편집]

막을 축전기 하나와 세 개의 전도도-전지 가지가 병렬로 붙은 회로로 본다. 전류 보존을 쓰면

C_m \frac{dV}{dt} = -\bar g_{Na} m^3 h\,(V - E_{Na}) - \bar g_K n^4 (V - E_K) - g_L (V - E_L) + I

이고, 게이팅 변수 $x \in \{m, h, n\}$ 는 각각 1차 완화 방정식

\frac{dx}{dt} = \alpha_x(V)\,(1 - x) - \beta_x(V)\,x = \frac{x_\infty(V) - x}{\tau_x(V)}

를 따른다. 여기서 $x_\infty = \alpha_x/(\alpha_x + \beta_x)$ , $\tau_x = 1/(\alpha_x + \beta_x)$ 다. 표준 파라미터는 $C_m = 1\ \mu\mathrm{F/cm^2}$ , $\bar g_{Na} = 120$ , $\bar g_K = 36$ , $g_L = 0.3\ \mathrm{mS/cm^2}$ 이고, 역전위는 현대 관례(정지막전위 $\approx -65$ mV)로 $E_{Na} = +50$ , $E_K = -77$ , $E_L = -54.4$ mV 를 쓴다.³

속도 계수 $\alpha, \beta$ 는 이론에서 유도된 게 아니라 실험 곡선에 맞춘 경험식이다. 예를 들어

\alpha_m(V) = \frac{0.1\,(V + 40)}{1 - e^{-(V+40)/10}}, \qquad \beta_m(V) = 4\,e^{-(V+65)/18}

같은 식들이고, $V = -40$ mV 에서 분모가 0이 되는 제거 가능한 특이점이 있어 코드에서 그대로 계산하면 0/0 이 뜬다. HH를 처음 구현하는 사람이 거의 전부 한 번씩 밟는 지뢰이며, 해결책은 그 근방에서 급수 전개로 갈아타는 것이다.

3. 스파이크가 만들어지는 구조[편집]

HH 모형의 핵심은 시간척도 분리다. 세 게이트의 시상수가 다음처럼 크게 갈린다.

변수	역할	시상수 $\tau$ (탈분극 시)
$m$	Na $^+$ 활성화	약 0.1~0.5 ms — 매우 빠름
$h$	Na $^+$ 불활성화	약 1~8 ms — 느림
$n$	K $^+$ 활성화	약 1~6 ms — 느림

이 배치가 스파이크의 모든 국면을 설명한다. 막이 조금 탈분극되면 $m$ 이 먼저 열려 Na $^+$ 유입이 늘고, 그 유입이 다시 탈분극을 키운다. $m^3$ 의 세제곱 비선형성이 붙은 이 양의 되먹임이 상승기를 폭발적으로 만든다. 뒤늦게 $h$ 가 닫히면서 나트륨 전도도를 꺼버리고, 동시에 $n^4$ 로 켜진 칼륨 전류가 막을 $E_K$ 쪽으로 끌어내려 하강기와 과분극(afterhyperpolarization)을 만든다.

문턱(threshold) 은 이 경쟁의 균형점이다. HH에는 “이 전압을 넘으면 발화”라는 딱 떨어지는 상수가 실제로는 없다 — 자극의 세기뿐 아니라 지속시간과 그 시점의 $h$ 값에 따라 유효 문턱이 움직인다. 정확히는 위상공간에서 안정 정지점의 흡인영역 경계(분리선 비슷한 준안정 다양체)가 문턱 역할을 한다. 불응기도 같은 그림에서 나온다. 스파이크 직후에는 $h$ 가 아직 낮고 $n$ 이 아직 높아서, 아무리 세게 찔러도 Na $^+$ 되먹임이 점화되지 않는다(절대 불응기). $h$ 가 회복되는 동안에는 더 센 자극이 필요하다(상대 불응기).

4. 제2형 흥분성과 호프 분기[편집]

정상전류 $I$ 를 서서히 올리며 발화 주파수를 재면 f–I 곡선이 0에서 시작하지 않는다. 표준 파라미터에서 HH는 발화 개시 순간 이미 수십 Hz 수준의 유한한 주파수로 뛰어오른다. 이런 성질을 호지킨이 1948년에 제안한 분류에서 제2형(Type II) 흥분성이라 부르며, 그 수학적 원인은 정지점이 호프 분기를 통해 불안정해지기 때문이다. 야코비안의 복소 켤레 고유값 쌍이 허수축을 건너면 진동이 태어나는데, 이때 진동수는 횡단 시점의 허수부 $\omega$ 로 이미 결정돼 있으므로 “임의로 느린 발화”가 불가능하다.

호지킨-헉슬리 4변수 막방정식(V, m, h, n)을 RK4 Δt = 0.01 ms로 적분한다. I = 0에서 막전위는 −64.9964 mV로 멎고 그때 게이트값은 해석식 x∞ = α/(α+β)와 상대오차 1e−13 이내로 일치한다. 가운데 패널의 g_Na = 120 m³h와 g_K = 36 n⁴는 위 패널 V(t)와 같은 시간축이다 — I = 10 μA/cm², 6.3 °C에서 g_K 정점은 g_Na 정점보다 1.27 ms 늦고, 그 지연이 곧 상승-재분극-후과분극(−74.9 mV)이다. 아래 f-I 곡선은 41개 전류값을 각각 독립 채널로 계단 입력해 실측해 쌓은 것이고, 발화는 0 Hz에서 시작하지 않는다 — I* = 6.2595 μA/cm²에서 50 Hz 대로 불연속하게 켜진다(제2형 흥분성). 위 패널의 붉은 띠는 실측 절대 불응기 2.90 ms다. 온도 슬라이더는 φ = 3^((T−6.3)/10)로 게이트 속도상수만 바꾸며, 21.58 °C를 넘으면 I = 10에서 스파이크가 서지 못한다. 공간이 없는 단일 등전위 구획이라 전파하는 활동전위가 아니라 한 점의 막전위다.

게다가 표준 HH의 호프 분기는 아임계다. 그래서 분기점 바로 아래 구간에서 안정 정지점과 안정 극한 순환이 공존하는 쌍안정 구간이 생기고, 자극을 올릴 때와 내릴 때 발화 개시·정지 전류가 달라지는 이력현상이 나타난다. 반대로 안장-마디 분기(SNIC)를 거치는 뉴런은 제1형이라 f–I 곡선이 0부터 연속적으로 자란다 — 같은 “발화”라도 분기 종류가 다르면 정보를 싣는 방식이 다르다는 뜻이고, 이 구분은 분기 이론이 신경과학에 실질적으로 기여한 대표 사례다.

5. 수치적 성질과 확장[편집]

시간척도가 $\tau_m \sim 0.1$ ms 부터 막 시상수 $\sim 10$ ms 까지 두 자릿수 벌어져 있어, HH는 전형적인 강성 방정식이다. 명시적 룽게-쿠타법으로 풀면 안정성 한계가 가장 빠른 성분에 잡혀 $\Delta t \lesssim 0.01$ ms 급의 잘게 썬 시간간격을 요구한다. 스파이크가 없는 정지 구간에서도 같은 $\Delta t$ 를 유지해야 하니 낭비가 크고, 그래서 실무 코드는 두 가지 요령을 쓴다. 하나는 게이팅 방정식이 $V$ 를 고정하면 선형이라는 점을 이용한 지수 적분(exponential Euler): $x^{n+1} = x_\infty + (x^n - x_\infty)e^{-\Delta t/\tau}$ 로 게이트만 정확히 풀고 $V$ 는 크랭크-니콜슨류로 푼다. 다른 하나는 $\alpha, \beta$ 를 전압 테이블로 미리 계산해 두는 것으로, 지수함수 호출이 전체 시간의 대부분을 먹는 상황을 바로 해결한다.⁴

공간으로 확장하면 케이블 방정식이 붙는다. 반지름 $a$ , 축방향 저항률 $R_i$ 의 원통 축삭에서

\frac{a}{2R_i}\frac{\partial^2 V}{\partial x^2} = C_m \frac{\partial V}{\partial t} + I_{\mathrm{ion}}(V, m, h, n)

이 되고, 여기서 나오는 진행파 해가 실제 신경 전도속도를 준다. 수치적으로는 축삭·수상돌기를 등전위 구획으로 잘라 각 구획에 HH 방정식을 걸고 인접 구획을 축저항으로 잇는 구획 모형(compartmental model)으로 푼다. 사실상 유한체적법 1차원판이며, NEURON·Arbor 같은 전용 시뮬레이터가 이 이산화 위에 세워져 있다. 확산항이 붙으므로 계는 자연스럽게 흥분성 매질이 되고, 진행파·차단·나선파 같은 현상이 그대로 따라온다.

미시 방향 확장도 있다. $m^3h$ 라는 곱은 “네 개의 독립 게이트가 다 열려야 통한다”는 가정인데, 실제 채널은 이산 상태 사이를 뛰어다니는 마르코프 사슬이다. 이를 그대로 모형화하면 나트륨 채널은 8상태, 칼륨 채널은 5상태 사슬이 되고, 채널 수가 유한하므로 전도도에 잡음이 실린다(채널 잡음). 이 확률적 HH는 문턱 근처에서 자발 발화나 확률 공명 같은 결정론 모형에는 없는 현상을 만든다. 큰 채널 수 극한에서 화학 랑주뱅 근사를 쓰면 확률미분방정식 형태로 떨어진다.

반대 방향, 즉 축약의 대표가 피츠휴-나구모 모형이다. $m$ 이 워낙 빨라 준정상 상태 $m \approx m_\infty(V)$ 로 놓고, $h + n \approx \mathrm{const}$ 라는 실측 관계를 써서 느린 두 변수를 하나로 합치면 4변수가 2변수(빠른 전압류 + 느린 회복 변수)로 줄어든다. 위상평면에서 널클라인 두 개로 흥분성·불응기·극한 순환을 전부 그림으로 설명할 수 있게 되는 대신, 이온별 전류라는 물리적 해석은 잃는다. 정량 예측이 필요하면 HH, 메커니즘 직관이나 대규모 네트워크·반응속도론적 공간 모형이 필요하면 축약형 — 이 선택이 계산 신경과학의 상시 트레이드오프다.

6. 관련 문서[편집]

공동 수상자는 시냅스 전달을 연구한 존 에클스. 참고로 앤드루 헉슬리는 소설 멋진 신세계를 쓴 올더스 헉슬리의 이복동생이자, 다윈의 불도그 토머스 헨리 헉슬리의 손자다. 유전자 풀이 억울할 정도로 편중돼 있다. ↩
1952년 논문의 수치 적분은 브런즈비가(Brunsviga) 기계식 계산기로 손수 돌린 것이다. 오늘날 노트북에서는 같은 궤적이 밀리초 단위로 나온다. “일단 돌려” 라는 말을 하려면 최소한 그 정도 격세지감은 알고 있어야 한다. ↩
원 1952년 논문은 탈분극을 음수로 잡고 정지막전위를 0으로 두는 옛 부호 관례를 썼다. 그래서 원문에 실린 값들을 현대 교과서 값과 직접 비교하면 부호가 뒤집혀 보인다. 이 문서는 전부 현대 관례( $V$ = 실제 막전위, 정지 $\approx -65$ mV)로 통일했다. 옛 논문 그림을 그대로 재현하려다 파형이 위아래로 뒤집혀 나온다면 십중팔구 이 문제다. ↩
프로파일러를 켜 보면 HH 구현의 병목이 미분방정식이 아니라 exp() 인 경우가 흔하다. 그래서 NEURON을 비롯한 전용 시뮬레이터들은 예외 없이 전압 룩업 테이블을 들고 다닌다. 테이블 간격을 너무 키우면 문턱 근처 거동이 미묘하게 달라지므로, 여기서도 결국 격자 수렴 확인이 필요하다. ↩