막전위 시뮬레이션 설명서

세포막을 사이에 두고 Na⁺는 세포 밖에, K⁺는 세포 안에 훨씬 많이 쌓여 있습니다(예: 세포외 Na⁺ 150mM vs 세포내 15mM, 세포내 K⁺ 150mM vs 세포외 5mM). 이 농도 차이는 저절로 생긴 게 아니라 Na⁺/K⁺ pump(나트륨-칼륨 펌프)가 ATP를 소모해가며 Na⁺ 3개를 퍼내고 K⁺ 2개를 퍼 넣는 능동수송을 끊임없이 반복한 결과입니다. 이렇게 만들어둔 농도 기울기 자체가 일종의 "저장된 에너지"이고, 이 에너지가 이온이 막을 통과해 확산하려는 힘(전기화학적 구동력)으로 바뀌면서 막전위가 만들어집니다.

이야기를 단순화해서, 막에 K⁺만 통과할 수 있는 통로가 있고 다른 이온은 전혀 못 지나간다고 해봅시다. K⁺는 농도가 높은 세포 안에서 밖으로 계속 확산하려 하겠지만, K⁺(양이온)가 빠져나갈수록 세포 안쪽이 점점 음전하를 띠면서 K⁺를 다시 안으로 끌어당기는 전기적 힘이 커집니다. 이 두 힘(농도차로 밀어내는 힘 vs 전하차로 끌어당기는 힘)이 정확히 균형을 이루는 지점의 전위를 그 이온의 평형전위(equilibrium potential)라 부르고, 네른스트 방정식으로 계산합니다:

이 시뮬레이션의 통계 타일에 뜨는 E_Na, E_K가 바로 이 값입니다. 기본값 기준으로 E_Na는 약 +60mV, E_K는 약 -90mV 근처가 나오는데, 이는 실제 뉴런에서 측정되는 값과 거의 일치합니다.

실제 세포막에는 Na⁺도, K⁺도 동시에 통과할 수 있는 통로가 있습니다. 이럴 때 막전위는 어느 한 이온의 평형전위로 고정되지 않고, 두 평형전위 사이 어딘가에서 결정됩니다. 정확히 어디쯤이냐를 계산하는 식이 골드만 방정식(2이온 단순화 버전)입니다:

네른스트 방정식과 모양이 비슷하지만, 각 이온의 농도에 막투과도(permeability, P)가 가중치로 곱해져 있다는 점이 다릅니다. 즉 골드만 방정식은 "E_Na와 E_K를 투과도 비율로 가중평균한 값"에 가깝습니다 — 정확히 산술평균은 아니지만(로그 안에서 가중되므로), 어떤 이온의 투과도가 크면 클수록 Vm이 그 이온의 평형전위 쪽으로 끌려간다는 직관은 그대로 성립합니다.

그래프의 가로축은 P_Na(K⁺ 투과도 P_K = 1을 기준으로 한 상대값)이고, 세로축은 그 조건에서 계산된 Vm입니다. 곡선을 왼쪽 끝(P_Na가 K⁺보다 훨씬 작을 때)부터 오른쪽 끝(P_Na가 K⁺보다 훨씬 클 때)까지 따라가 보면:

• 왼쪽 끝 — P_Na ≈ 0.001: 사실상 K⁺만 통과하는 상황과 같아서 Vm이 E_K에 거의 딱 붙습니다. 실제 휴지 상태의 뉴런이 이 부근에 있습니다(P_Na가 P_K의 약 1/20 수준) — 그래서 안정막전위(약 -70mV 근처)가 E_K(약 -90mV)에 가깝지만 완전히 같지는 않은 겁니다.
• 오른쪽 끝 — P_Na가 20 이상: 이번엔 반대로 Vm이 E_Na에 가까워집니다. 실제 활동전위의 정점 부근이 이 상황입니다 — 전압개폐성 Na⁺ 통로가 순간적으로 대량으로 열리면서 P_Na가 P_K보다 훨씬 커지고, Vm이 음수에서 양수로 확 뒤집힙니다(탈분극).

즉 Na⁺와 K⁺는 항상 서로 반대 방향으로 Vm을 끌어당기고 있고, 실제로 어느 쪽이 이기느냐는 그 순간 막이 어느 이온을 더 잘 통과시키느냐(투과도)에 달려 있습니다. "탈분극 — Na⁺ 채널 열림" 버튼을 눌러서 P_Na를 확 올려보면, 곡선을 따라 점이 오른쪽으로 이동하면서 Vm이 음수에서 양수로 뒤집히는 걸 바로 확인할 수 있습니다.

쌤 한마디: 재분극은 이 그래프의 반대 방향 이동으로 설명됩니다. Na⁺ 통로가 금방 불활성화(inactivation)되고 뒤이어 K⁺ 통로가 열리면서 P_Na가 다시 낮아지고 P_K가 일시적으로 더 커지면, 점이 다시 왼쪽 아래(E_K 근처, 심지어 E_K보다 더 아래인 과분극)로 돌아갑니다. 이 시뮬레이션은 P_K를 1로 고정해뒀지만, 개념은 "P_Na를 다시 낮춰보는 것"으로도 충분히 재현해볼 수 있습니다.

골드만 방정식 자체는 "이 순간의 조건이라면 막전위가 얼마다"를 계산하는 스냅샷 공식이라 시간 개념이 없습니다. 그런데 실제 세포막은 축전기(capacitor)처럼 전하를 저장하고 있어서, 투과도가 바뀐다고 막전위가 그 순간 바로 새 값으로 점프하지 않습니다 — 마치 콘덴서가 충전·방전되듯 지수적으로 새 목표값을 향해 서서히 따라갑니다. 이 관계를 식으로 쓰면:

τ(막 시간상수)는 막이 얼마나 빨리 반응하는지를 나타내는 값으로, 실제로는 막저항×막전기용량(τ = Rm × Cm)으로 정해집니다. τ가 작을수록(막이 얇고 저항이 낮을수록) 반응이 빠르고, τ가 클수록 느긋하게 변합니다.

재생 버튼을 누르면 "시간에 따른 막전위 변화" 그래프가 그려지기 시작합니다. 회색 점선은 지금 슬라이더 조건에서 계산된 목표값(슬라이더를 움직이면 즉시 계단처럼 뛰어오름/떨어짐)이고, 주황색 실선이 실제로 이 시뮬레이션이 τ를 반영해 계산한 막전위입니다. 재생 중에 "탈분극" 버튼을 눌러보면 목표값은 즉시 뛰지만 실선은 곡선을 그리며 뒤따라가는 게 보이고, 다시 "휴지 상태"를 누르면 원래대로 돌아오는 하강 곡선까지 볼 수 있습니다 — 두 버튼을 연달아 눌러보면 실제 활동전위와 닮은 파형(상승 → 짧은 유지 → 하강)이 그려집니다.

농도 슬라이더는 곡선 자체와 E_Na·E_K 점선을 통째로 움직입니다. 예를 들어 [K⁺] 세포외 농도를 올리면(고칼륨혈증과 비슷한 상황) E_K가 위로 올라가면서(음수값이 작아지면서) 휴지 상태 Vm도 함께 양의 방향으로 밀립니다 — 실제로 혈중 칼륨 농도가 비정상적으로 높아지면 심장 세포의 안정막전위가 불안정해져서 부정맥 위험이 커지는 것도 같은 원리입니다.

• 주황색 점은 지금 P_Na 슬라이더 값에서 실제로 계산된 (P_Na, Vm) 위치입니다.
• 가로축은 로그 스케일입니다 — 칸 하나를 지날 때마다 P_Na가 10배씩 커집니다. Na⁺ 투과도는 휴지 상태와 탈분극 상태 사이에 수백~수천 배 차이가 나기 때문에, 로그 스케일이 아니면 두 상태를 한 그래프에 같이 담을 수 없습니다.
• 회색 점선 두 개는 각각 E_Na, E_K입니다 — 곡선이 이 두 선 사이에서만 움직인다는 걸 눈으로 확인할 수 있습니다.

• 2이온 단순화 — 실제 세포막에는 Cl⁻ 등 다른 이온도 관여하지만, 이 시뮬레이션은 교육과정에서 가장 비중 있게 다루는 Na⁺·K⁺ 두 이온만으로 단순화했습니다.
• 일정한 전기장 가정(constant field assumption) — 골드만 방정식 자체가 막 안에서 전기장이 균일하다고 가정하고 유도된 근사식입니다. 실제로는 막 단백질 구조 때문에 완벽히 균일하지는 않습니다.
• 단순화된 시간 반응(τ 고정) — 재생 그래프는 막의 RC 반응(dVm/dt = (목표−Vm)/τ)까지는 반영하지만, τ 하나로 상승·하강을 모두 표현합니다. 실제 뉴런은 Na⁺ 통로가 열리는 속도·닫히는 속도(불활성화)·K⁺ 통로가 열리는 속도가 전부 다르고 심지어 전압에 따라서도 달라져서, 활동전위의 진짜 비대칭적인 파형(가파른 상승, 완만한 하강, 후과분극)을 정확히 재현하려면 호지킨-헉슬리 모델처럼 통로마다 별도의 시간 변수가 필요합니다.
• P_Na는 사용자가 직접 조절, 실제로는 전압에 따라 자동으로 열림 — 실제 뉴런은 Vm 자체가 어떤 전압 문턱을 넘으면 전압개폐성 Na⁺ 통로가 저절로 열립니다(양성 되먹임). 이 시뮬레이션은 그 자동 개폐 과정 대신 P_Na를 사람이 슬라이더·버튼으로 직접 조작하게 해서 "투과도가 바뀌면 막전위가 어떻게 반응하는가"라는 핵심 관계만 떼어서 보여줍니다.
• 온도 고정 — 사람 체온(37°C)으로 고정했습니다. 온도가 낮아지면 RT/F 값이 작아져서 같은 농도차에서도 평형전위의 절댓값이 살짝 줄어듭니다.