Ця симуляція реалізує класичну модель Ходжкіна-Хаклі (ХХ) 1952 року — визначне досягнення обчислювальної нейронауки, яке принесло Алану Ходжкіну та Ендрю Хаклі Нобелівську премію з фізіології та медицини 1963 року. Модель описує, як ділянка мембрани велетенського аксону кальмара генерує потенціали дії через взаємодію трьох іонних струмів: швидкого вхідного натрієвого (Na+), що контролюється воротами активації m та інактивації h; повільнішого вихідного калієвого (K+), що контролюється воротою n; та пасивного омічного струму витоку. Ці струми об'єднані рівнянням Cm dV/dt = Iext − INa − IK − IL, яке чисельно інтегрується методом Рунге-Кутти четвертого порядку з кроком 0,02 мс. Перемикайтеся між режимами Напруга, Ворота та Фаза кнопками на панелі керування, а слайдери дозволяють досліджувати, як зовнішній струм, температура та провідності каналів впливають на збудливість нейрона.
Три режими відображення розкривають різні аспекти динаміки. Режим Напруга показує послідовність потенціалів дії у вікні 150 мс: підніміть Iext вище ~6 мкА/см² щоб побачити регулярні розряди, та спостерігайте, як частота зростає зі струмом (крива F-I). Режим Ворота одночасно відображає всі три змінні: m (золота) різко зростає під час деполяризації, h (рожева) спадає для рефрактерного періоду, а n (зелена) повільно зростає для реполяризації. Режим Фаза показує фазовий портрет V–n з нуль-клинами: бурштинова нуль-клина n (dn/dt = 0) та рожева нуль-клина V (dV/dt = 0) перетинаються у нерухомій точці стану спокою, а при повторних розрядах траєкторія обертається навколо неї як стійкий граничний цикл.
Що таке модель Ходжкіна-Хаклі?
Модель Ходжкіна-Хаклі (1952) — система чотирьох диференціальних рівнянь, що описує зміну потенціалу мембрани нейрона через швидкий натрієвий (Na+) та повільний калієвий (K+) струми разом із пасивним струмом витоку. За цю роботу автори отримали Нобелівську премію з фізіології та медицини 1963 року і донині залишаються золотим стандартом моделювання одиночного нейрона.
Що таке змінні воріт m, h і n?
m — ворота активації Na+ (швидко відкриваються при деполяризації), h — ворота інактивації Na+ (закриваються після спайку, зумовлюючи рефрактерний період), n — ворота активації K+ (відкриваються повільно для реполяризації). Кожна змінна підпорядковується кінетиці першого порядку з залежними від напруги константами швидкості, отриманими з оригінальних дослідів voltage-clamp на гігантському аксоні кальмара.
Як зовнішній струм I_ext генерує потенціали дії?
Нижче порогу реобази (~6 мкА/см²) мембрана повертається до стану спокою після підпорогового відхилення. Вище цього порогу деполяризація прискорює активацію Na+ швидше, ніж інактивація встигає закрити канали, генеруючи регенеративний спайк амплітудою ~100 мВ. Збільшення I_ext підвищує частоту розрядів — це крива F-I, ключова характеристика нейрональної збудливості.
Одразу після спайку ворота h (інактивація Na+) закриті, а ворота n (K+) максимально відкриті, тому другий спайк неможливий незалежно від величини стимулу. Цей абсолютний рефрактерний період тривалістю ~1–2 мс забезпечує однонаправлене поширення потенціалів дії вздовж аксонів.
Температура враховується через масштабування Q10: константи швидкості воріт множаться на 3^((T−6,3)/10), де 6,3 °C — початкова експериментальна температура. Вища температура пришвидшує воротну кінетику, скорочуючи тривалість спайку та збільшуючи максимальну частоту розрядів, але також підвищує поріг реобази. Вище ~31 °C модель може перейти в деполяризаційний блок.
Фазовий портрет V–n відображає траєкторію потенціалу мембрани проти змінної активації K+, відстежуючи граничний цикл під час повторних розрядів. Нуль-клина V (dV/dt=0) та нуль-клина n (dn/dt=0) перетинаються у нерухомій точці стану спокою. Коли I_ext перевищує поріг, нерухома точка стає нестійкою, і траєкторія обертається навколо неї як стійкий граничний цикл.
Провідність Na+ забезпечує вхідний регенеративний струм для підйому потенціалу дії. Зменшення gNa знижує максимальний вхідний струм. Нижче критичного значення деполяризація більше не може подолати вихідні струми K+ і витоку, перетворюючи збудливий стан на стійкий вузол — спайки не генеруються незалежно від I_ext.
Якщо струм наростає повільно замість ступінчастого включення, ворота h встигають інактивуватися при підпорогових напругах до досягнення порогу. Це прогресивне закриття Na+-каналів підвищує ефективний поріг, тому повільно зростаючий стимул може не генерувати спайк, хоча те саме кінцеве значення як ступінь викликало б розряд.
Симуляція використовує метод Рунге-Кутти 4-го порядку (RK4) з фіксованим кроком 0,02 мс, що значно менше найшвидших часових констант (~0,1 мс для m). RK4 забезпечує точність четвертого порядку та відмінну стабільність для жорстких рівнянь ХХ, даючи результати, що збігаються з опублікованими аналітичними еталонами.
Так. Вона є основою всіх кондуктансних моделей нейронів. Розширення включають багатокомпартментні моделі (дендрити + сома + аксон), кальцієві канали, M-струми та тисячі типів іонних каналів у симуляторах NEURON і Brian2. Великомасштабні проєкти мозку (Blue Brain, Human Connectome) зрештою спираються на рівняння типу ХХ для опису поведінки окремих нейронів.