⚡ Neuron Hodgkina-Huxleya
Model neuronu oparty na przewodnościach (Nobel 1963). Podaj prąd i obserwuj, jak zmienne bramkowania Na⁺ i K⁺ — m, h, n — napędzają potencjał czynnościowy, okres refrakcji i wyładowania powtarzalne.
O modelu neuronu Hodgkina-Huxleya
Ta symulacja całkuje klasyczne równania Hodgkina-Huxleya, opisujące sposób, w jaki błona nerwowa generuje potencjał czynnościowy. Pojedyncze równanie pojemnościowe, C dV/dt = I − INa − IK − IL, jest sprzężone z trzema zmiennymi bramkowania m, h i n. Prąd Na⁺ skaluje się jak gNa m³h, a prąd K⁺ jak gK n⁴, przy czym zależne od napięcia stałe szybkości są tu rozwiązywane metodą Eulera w przód z krokiem czasowym 0,01 ms.
Panel Bodziec ustawia wstrzykiwany prąd Iext (0–50 µA/cm²), moment rozpoczęcia impulsu, czas trwania impulsu oraz okno symulacji, a presety takie jak Pojedynczy impuls, Powtarzalny czy Podprogowy wczytują typowe reżimy pracy. Wykres przedstawia napięcie błonowe nad krzywymi bramkowania m, h i n, a odczyty pokazują napięcie szczytowe, liczbę impulsów i częstość wyładowań. Hodgkin i Huxley otrzymali za tę pracę Nagrodę Nobla w 1963 roku; pozostaje ona fundamentem neuronauki obliczeniowej.
Najczęściej zadawane pytania
Co opisuje model Hodgkina-Huxleya?
Jest to model oparty na przewodnościach, opisujący sposób, w jaki błona neuronu wytwarza potencjał czynnościowy, czyli impuls nerwowy. Traktuje błonę jako kondensator połączony równolegle z bramkowanymi napięciowo kanałami sodowymi i potasowymi oraz biernym wyciekiem, i śledzi, jak napięcie błonowe zmienia się w czasie milisekund po wstrzyknięciu prądu.
Czym są zmienne bramkowania m, h i n?
Są to bezwymiarowe prawdopodobieństwa z przedziału od 0 do 1, opisujące stan kanału. Zmienna m to aktywacja sodowa, a h to inaktywacja sodowa, więc przewodność Na⁺ podlega m³h; n to aktywacja potasowa, więc przewodność K⁺ podlega n⁴. Każda z nich dąży do zależnego od napięcia stanu ustalonego z własnym tempem.
Co robią elementy sterujące na tej stronie?
Suwak wstrzykiwanego prądu ustawia siłę bodźca w µA/cm², natomiast suwaki momentu rozpoczęcia i czasu trwania impulsu określają, kiedy i jak długo prąd jest podawany. Suwak okna symulacji ustala całkowity pokazywany czas, a przyciski presetów przechodzą do typowych reżimów wyładowań. Każda zmiana natychmiast uruchamia ponowne obliczenia.
Dlaczego impuls pojawia się dopiero powyżej progu?
Poniżej około 6,5 µA/cm² niewielki napływ Na⁺ jest przeważany przez prądy K⁺ i wyciekowe, więc błona ładuje się nieznacznie, a następnie biernie wraca do spoczynku. Powyżej progu dodatnie sprzężenie zwrotne aktywacji sodowej wyprzedza prądy wypływające, wyzwalając skokowy, typu wszystko-albo-nic, wzrost potencjału czynnościowego.
Co powoduje okres refrakcji?
Po impulsie bramka inaktywacji sodowej h jest zamknięta, a bramka potasowa n jest wciąż otwarta, co utrzymuje błonę w stanie hiperpolaryzacji. W tej bezwzględnej fazie refrakcji drugi impuls nie może się pojawić, niezależnie od wielkości prądu, ponieważ kanały sodowe są tymczasowo niedostępne, dopóki h się nie zregeneruje.
Dlaczego napięcie przekracza około +50 mV?
Gdy kanały sodowe się otwierają, potencjał błonowy jest kierowany w stronę potencjału równowagi sodu ENa, ustawionego tutaj na +50 mV. Szybka przewodność sodowa m³h ciągnie napięcie w górę od poziomu spoczynkowego −65 mV, wytwarzając ostry, dodatni przerzut, zanim potas repolaryzuje komórkę.
Jak model wytwarza wyładowania powtarzalne?
Jeśli wstrzykiwany prąd jest podtrzymywany i wystarczająco silny, błona jest wielokrotnie ponownie popychana powyżej progu po każdym impulsie. Gdy bramka h regeneruje się, a błona ponownie się depolaryzuje, pojawia się kolejny potencjał czynnościowy, dając serię impulsów, których częstość rośnie wraz z prądem bodźca.
Czy ta symulacja jest fizycznie dokładna?
Wykorzystuje oryginalne parametry olbrzymiego aksonu kalmara z 1952 roku w temperaturze 6,3 °C, w tym gNa = 120, gK = 36 i gL = 0,3 mS/cm², wraz z opublikowanymi funkcjami szybkości α i β. Dynamika jest wierna, choć integracja metodą Eulera w przód i stała temperatura sprawiają, że jest to model dydaktyczny, a nie dopasowanie do konkretnego neuronu ssaczego.
Jaka metoda numeryczna rozwiązuje te równania?
Cztery sprzężone równania różniczkowe dla V, m, h i n są rozwiązywane metodą Eulera w przód z niewielkim krokiem czasowym 0,01 ms. W każdym kroku prądy jonowe są obliczane na podstawie bieżącego napięcia i bramek, wyznaczane są pochodne, a wszystkie cztery zmienne są aktualizowane przed kolejnym krokiem.
Dlaczego model ten ma znaczenie wykraczające poza neuronaukę?
Struktura Hodgkina-Huxleya jest wzorcem dla niemal wszystkich biofizycznych modeli neuronów i stanowi podstawę symulacji komórek pobudliwych, w tym mięśnia sercowego. Jest szeroko stosowana do badania wpływu leków na kanały jonowe, padaczki i arytmii, i pozostaje standardowym przykładem dydaktycznym nieliniowej dynamiki i układów pobudliwych.
Model neuronu oparty na przewodnościach (Nobel 1963). Podaj prąd i obserwuj, jak zmienne bramkowania Na⁺ i K⁺ — m, h, n — napędzają potencjał czynnościowy, okres refrakcji i wyładowania powtarzalne.
2D · HTML5 Canvas 2D · cel 60 FPS · działa w całości po stronie klienta, bez instalacji