이온 채널

1. 개요[편집]

이온 채널 Ion channel
정체	막을 관통하는 이온 선택적 구멍 단백질
단일 채널 전도도	약 1~200 pS (K+ 채널 대체로 10~100 pS)
단일 채널 전류	0.1~10 pA — 사각파로 관측됨
수송률	107~108 ions/s (확산 한계 근처)
표준 모형	연속시간 마르코프 사슬 (Q-행렬)
측정법	패치 클램프 (네어·자크만, 1976)

이온 채널(ion channel)은 세포막을 관통해 특정 이온만 선택적으로 통과시키는 구멍을 가진 막단백질로, 열림과 닫힘 사이를 무작위로 오가는 두 상태 이상의 확률 스위치로 동작한다. 활동전위의 되먹임을 만드는 것도, 불응기를 만드는 것도 결국 이 스위치들의 집단 통계다.

거시적 관점 — 막전위와 게이팅 변수로 쓴 결정론적 방정식 — 은 호지킨-헉슬리 모형이 이미 다루고, 파형과 전파는 활동전위·케이블 방정식에 있다. 이 문서는 그 아래층, 채널 하나가 실제로 무엇을 하고 있으며 그것을 어떻게 수치적으로 다루는가를 다룬다. 요점 한 줄로 줄이면 이렇다: $g_K = \bar g_K n^4$ 의 $n^4$ 은 물리량이 아니라 어떤 마르코프 사슬의 정상해를 요약한 대수식이고, 채널 수가 유한하면 그 요약이 깨진다.

2. 단일 채널 기록 — 사각파가 보인다[편집]

1976년 네어(Neher)와 자크만(Sakmann)이 유리 미세피펫 끝을 세포막에 밀착시켜 기가옴 밀봉(giga-seal)을 만드는 데 성공하면서, 사람들은 처음으로 채널 한 개의 전류를 봤다. 결과는 충격적으로 단순했다. 부드러운 곡선이 아니라 거의 완벽한 사각파였다. 채널은 반쯤 열리지 않는다 — 닫혀 있거나(0 pA), 열려 있거나(예컨대 −2 pA), 둘뿐이고 그 사이 전이는 측정 대역폭보다 빠르다.¹

열려 있는 동안의 전류는 단순한 옴 법칙 꼴이다.

$$i = \gamma\,(V - E_{\rm rev})$$

$\gamma$ 가 단일 채널 전도도(수십 pS 규모)이고 $E_{\rm rev}$ 는 그 이온의 역전위다. 전압을 바꿔 가며 $i$ 를 재서 그린 $i$–$V$ 직선의 기울기가 $\gamma$, $x$ 절편이 $E_{\rm rev}$ 다. 즉 전도도는 게이팅과 무관한 상수이고, 게이팅은 오직 “열려 있는 시간 비율”만 바꾼다.

그래서 거시 전류는 곱으로 분해된다.

$$I = N\,P_o\,i$$

$N$ 은 패치 안의 채널 수, $P_o$ 는 열림 확률이다. 전압 고정 실험에서 같은 자극을 수백 번 반복해 단일 채널 기록을 평균 내면 매끄러운 거시 전류 파형이 그대로 복원된다 — 이것이 HH의 게이팅 변수에 미시적 의미를 부여한 결정적 실험이다. $n^4$ 은 결국 $P_o$ 의 근사식이었다.

$N$ 과 $i$ 를 따로 모르는 상태에서도 둘을 뽑아내는 표준 수법이 비정상 잡음 해석(nonstationary noise analysis)이다. 채널들이 독립이면 열린 개수가 이항분포이므로 분산과 평균 사이에

$$\sigma^2 = i\,I - \frac{I^2}{N}$$

라는 아래로 볼록한 포물선 관계가 성립한다. 반복 시행의 앙상블 분산을 평균 전류에 대해 찍고 이 포물선을 맞추면 초기 기울기에서 $i$, 곡률에서 $N$ 이 나온다. 채널 하나를 보지 않고도 채널 하나의 성질을 재는 이 트릭은 최소자승법 맞춤 문제 하나로 환원되며, 요즘도 표준 절차다.

3. 게이팅은 마르코프 사슬이다[편집]

단일 채널 기록에서 열림 지속시간의 히스토그램을 그리면 지수분포(혹은 지수의 합)가 나온다. 이것이 결정적 힌트다 — 채널의 다음 행동은 현재 상태에만 의존하고 얼마나 오래 그 상태에 있었는지에는 의존하지 않는다. 즉 게이팅은 이산 상태 위의 연속시간 마르코프 사슬(마르코프 사슬)이고, 상태 점유 확률 벡터 $\mathbf{p}(t)$ 는 마스터 방정식

$$\frac{d\mathbf{p}}{dt} = \mathbf{p}\,Q(V), \qquad \sum_j Q_{ij} = 0$$

를 따른다. 콜쿤과 호크스가 정리한 이 $Q$-행렬 형식이 채널 동역학의 공용어다. 전압을 고정하면 $Q$ 가 상수이므로 해는 행렬 지수 $\mathbf{p}(t) = \mathbf{p}(0)e^{Qt}$ 로 정확히 나온다.

HH의 곱 형태는 이 사슬의 특수한 경우로 정확히 환원된다. 칼륨 채널의 $n^4$ 은 네 개의 동일하고 독립인 서브유닛이 각각 확률 $n$ 으로 열린다는 뜻이고, 열린 서브유닛 개수를 상태로 삼으면 5상태 사슬이 된다.

$$C_0 \underset{\beta}{\overset{4\alpha}{\rightleftharpoons}} C_1 \underset{2\beta}{\overset{3\alpha}{\rightleftharpoons}} C_2 \underset{3\beta}{\overset{2\alpha}{\rightleftharpoons}} C_3 \underset{4\beta}{\overset{\alpha}{\rightleftharpoons}} O$$

앞의 계수 $4,3,2,1$ 은 “아직 닫힌 서브유닛이 몇 개 남았는가”라는 조합 인자다. 같은 논리로 $m^3h$ 는 $m$ 쪽 4상태와 $h$ 쪽 2상태의 곱, 즉 8상태 격자 사슬이 된다.

여기서 곱 형태의 한계가 드러난다. 곱은 활성화와 불활성화가 독립이라고 못 박는다. 실제 나트륨 채널의 불활성화는 활성화가 진행된 뒤라야 잘 일어나는 결합(coupled) 과정이라, 8상태 격자에서 몇몇 전이율이 0에 가깝거나 방향이 한쪽으로 쏠린 사슬이 실측에 훨씬 가깝다. 결과의 차이는 특히 불활성화로부터의 회복 시간 경과에서 나타나고, 그래서 회복이 곧 관심사인 심장 모형(부정맥 참고)에서는 마르코프 사슬 게이팅이 사실상 표준이다.

일반 사슬이 필요한 다른 이유도 있다.

• 상태 의존 약물 결합. 항부정맥제·국소마취제는 열린 상태나 불활성 상태에만 붙는다. “채널을 $x$ % 막는다”로는 표현할 수 없고, 결합·해리 상태를 사슬에 노드로 추가해야 사용 의존성(use-dependence)이 재현된다.
• 채널병증. 특정 변이가 특정 전이율만 바꾸는 경우, 사슬 위에서는 화살표 하나의 값이지만 곱 형태에서는 대응하는 파라미터가 아예 없다.
• 소프트웨어. NEURON의 NMODL이 KINETIC 블록을 제공하는 이유가 이것이고, 도구는 사슬을 받아 자동으로 야코비안까지 만들어 준다.

대가는 비용과 식별성이다. 상태가 늘면 전이율 파라미터가 제곱으로 늘어나는데 관측되는 것은 열림·닫힘 두 준위뿐이라, 서로 다른 사슬이 같은 데이터를 똑같이 잘 설명하는 일이 흔하다. 실질적으로 이것은 은닉 마르코프 모형 식별 문제이며, 사슬 위상을 데이터만으로 정하려는 시도는 대개 실패한다. 구조 정보와 실험 프로토콜(전압 계단, 회복 간격 쓸기)로 제약을 걸어야 한다.

4. 선택성 필터 — 어떻게 빠르면서 까다로운가[편집]

이온 채널의 역설은 선택성과 처리량을 동시에 가진다는 데 있다. K$^+$ 채널은 자기보다 작은 Na$^+$ 를 1000:1 이상으로 걸러내면서, 동시에 초당 $10^7$~$10^8$ 개를 통과시켜 자유 확산 한계에 근접한다. 보통 이런 조합은 불가능하다 — 까다로운 필터는 느리기 마련이다.

1998년 매키넌 그룹이 세균 K$^+$ 채널 KcsA의 원자 구조를 풀면서 답이 나왔다. 선택성 필터는 TVGYG라는 보존 서열의 주쇄 카보닐 산소들이 만든 좁은 통로이고, 이 산소들의 배치가 물 분자가 K$^+$ 를 둘러싼 수화껍질의 기하를 거의 그대로 흉내 낸다. K$^+$ 는 물을 벗은 대가를 필터가 대신 지불해 주므로 에너지 손해 없이 통과하고, 반지름이 작은 Na$^+$ 는 그 배치에 맞지 않아 탈수화 비용을 회수하지 못한다. 크기로 거르는 체가 아니라 배위 기하로 거르는 주형인 셈이다.

빠른 이유는 따로 있다. 필터 안에는 결합 자리가 네 개 있고, 두 개 이상이 동시에 점유되면 이온 사이의 정전 반발이 다음 이온을 밀어낸다 — 밀어내기 전도(knock-on)라 부르는 다중이온 메커니즘이다. 하나가 들어오면 반대쪽에서 하나가 튀어나가는 방식이라 결합의 깊이가 통과 속도를 죽이지 않는다.

이 층위는 상미분방정식이 다룰 수 있는 규모가 아니다. 전도 메커니즘 자체를 계산하려면 분자동역학이나 그 위의 자유에너지 계산으로 내려가야 하고, 반대로 세포·조직 시뮬레이션에서는 이 모든 것이 $\gamma$ 라는 상수 하나로 압축된다. 어느 층위에서 무엇을 상수로 취급할지 정하는 것이 다중 스케일 모형화의 전부라 해도 된다.

5. 채널 잡음과 스파이크 지터[편집]

$N$ 이 유한하다는 사실은 그냥 넘어갈 수 있는 디테일이 아니다. 열린 채널 수는 이항분포이므로 전도도의 상대 요동은 $1/\sqrt{N}$ 로 줄어들고, $N$ 은 막 면적에 비례한다. 즉 잡음은 면적의 제곱근에 반비례한다. 큰 세포체에서는 무시할 만하지만 가는 축삭이나 수상돌기 가시(spine)에서는 지배적이다.

구체적으로 무엇이 벌어지는가.

• 역치 근처의 자발 발화. 결정론 모형이라면 아무 일도 없을 정상 상태에서, 나트륨 채널이 우연히 여러 개 동시에 열려 되먹임이 점화되면 자극 없이 스파이크가 난다. 축삭이 가늘수록 빈발한다.
• 스파이크 시각 지터. 같은 자극을 반복해도 발화 시각이 흩어진다. 램프처럼 천천히 올라가는 자극일수록 역치 통과 시각이 잡음에 민감해 지터가 커지고, 계단처럼 급한 자극에서는 작다.
• 전도 실패와 축삭 지름의 하한. 축삭이 가늘어지면 잡음이 커져 엉뚱한 자리에서 스파이크가 나거나 진행 중인 스파이크가 죽는다. 이 때문에 무수축삭 지름에는 채널 잡음이 정하는 물리적 하한(0.1 µm 안팎)이 있다는 분석이 있고, 실측된 가장 가는 축삭들이 대략 그 근처에 있다.
• 확률 공명. 적당한 세기의 잡음이 역치 아래 신호의 검출을 오히려 개선한다. 잡음을 0으로 만드는 것이 언제나 이득은 아니라는, 이 바닥에서 가장 반직관적인 결과.

6. 확률적 게이팅을 수치적으로 돌리기[편집]

채널 잡음을 시뮬레이션에 넣는 방법은 정확도-비용 축 위에 정렬된다.

정확한 사건 구동. 각 채널의 상태를 정수로 들고, 다음 전이까지의 대기시간을 지수분포에서 뽑아 사건 단위로 전진한다(길레스피 알고리즘의 게이팅판). 채널이 많아지면 사건이 비례해서 늘어 비싸진다. 더 성가신 것은 전이율이 $V$ 를 통해 시간에 의존한다는 점이다. 표준 SSA의 대기시간 공식이 성립하지 않으므로, 스텝 안에서 $V$ 를 동결하거나 상한율을 잡아 기각 표집(thinning)을 해야 하고, 전자를 무비판적으로 쓰면 스파이크 상승기에서 편향이 생긴다.

이항 근사. 고정 스텝 $\Delta t$ 동안 각 전이의 발생 개수를 이항분포에서 한 번에 뽑는다. $\Delta t$ 안에 전이율이 거의 변하지 않는다는 가정이 필요하고, 여러 전이가 같은 상태의 채널을 놓고 경쟁할 때 개수가 음수가 되지 않도록 순서와 상한을 관리해야 한다. 비용이 채널 수와 무관해져 실용적이다.

랑주뱅 근사. $N$ 이 크면 이항 요동을 가우스로 바꿔 확률미분방정식으로 쓴다. 여기서 유명한 함정이 하나 있다 — 잡음을 게이팅 변수 $m,h,n$ 각각에 독립으로 붙이는 “서브유닛 기반” 구현은 틀린다. 실제 요동은 채널의 상태 점유수에 있고 서브유닛 조합을 거치며 상관을 가지므로, 잡음 항은 사슬의 상태별 전이에 붙여 확산 행렬을 그대로 써야 한다. 두 방식은 발화율 통계에서 눈에 띄게 갈리고, 이 구분을 흐린 옛 구현들이 한동안 문헌을 오염시켰다. 추가로 $[0,1]$ 이나 $[0,N]$ 경계에서 표본이 밖으로 나가므로 반사·절단 처리가 필요하다.²

결정론. $N \to \infty$ 극한. 채널 잡음이 관심사가 아니면 이것이 정답이고, 게이트는 지수 적분기 계열(러시-라슨)로 값싸게 정확히 밀 수 있다.

성능 얘기를 하나 덧붙이면, 상세 모형의 프로파일을 떠 보면 병목이 미분방정식이 아니라 전이율 계산의 exp() 인 경우가 흔하다. 전압 룩업 테이블을 미리 깔아 두는 것이 관행이고, 다만 테이블 간격이 성기면 활동전위 파형에 계단 오차가 실린다. 이것도 결국 메시 독립성 연구와 같은 문법의 수렴 검사 대상이다.³

7. 관련 문서[편집]

• 활동전위 · 호지킨-헉슬리 모형 · 불응기
• 케이블 방정식 · 구획 모형 · NEURON
• 부정맥 · 제세동 · 흥분성 매질
• 확률미분방정식 · 은닉 마르코프 모형 · 반응속도론
• 분자동역학 · 계산 신경과학 · 시냅스 가소성
• 지수 적분기 · 강성 방정식 · 최소자승법

8. Footnotes[편집]

1. 기가옴 밀봉이라는 이름은 피펫과 막 사이 누설 저항이 $10^9\ \Omega$ 을 넘는다는 뜻이다. pA 급 전류를 재려면 이 정도는 돼야 한다. 처음에는 밀봉이 잘 안 돼서 애를 먹다가 피펫을 깨끗이 하고 살짝 음압을 걸었더니 저항이 수백 배 뛰었다는 것이 이 분야의 창세기 설화다. 실험 팁 하나가 노벨상(1991)까지 간 흔치 않은 사례. ↩

2. “잡음을 어디에 붙이느냐”는 문제는 채널에만 있는 게 아니다. 미시 모형의 요동을 거시 변수로 옮길 때 변수 변환이 잡음의 상관 구조를 바꾼다는 것은 확률미분방정식 일반의 상식인데, 편리한 변수($m,h,n$)와 요동이 실제로 사는 변수(상태 점유수)가 다르면 이 상식을 밟게 된다. 결정론 극한에서는 두 구현이 같은 식으로 수렴하기 때문에 코드 리뷰로는 안 잡히고, 발화율 분포를 찍어 봐야 보인다. ↩

3. 채널이 유한하다는 사실이 반가운 순간도 있다. 결정론 모형은 초기조건이 완벽히 대칭이면 영원히 대칭인 답을 내놓는데, 실제 계는 잡음이 대칭을 깨 준다. 대칭 깨짐이 필요한 계산에서 난수를 일부러 넣는 대신 물리적으로 정당한 채널 잡음을 켜면 되는 셈이라, “잡음도 모형의 일부”라는 말이 여기서는 변명이 아니다. ↩