Модель Ходжкина — Хаксли

Модель Ходжкина — Хаксли — это математическая модель, описывающая генерацию и распространение потенциала действия в гигантском аксоне кальмара. Разработанная Аланом Ходжкином и Эндрю Хаксли в 1952 году, она является фундаментальной основой современной нейробиологии и биофизики. Модель представляет собой систему нелинейных дифференциальных уравнений, которая количественно описывает, как изменения ионной проводимости мембраны приводят к возникновению электрического сигнала.

История

В 1930–1940-х годах физиологи активно изучали механизмы возбудимости нервных клеток. Ключевым объектом для исследований стал гигантский аксон кальмара (Doryteuthis pealeii), диаметр которого (до 1 мм) позволял вводить внутрь микроэлектроды и управлять ионным составом среды. В 1939 году Ходжкин и Хаксли впервые зарегистрировали потенциал действия внутри аксона, подтвердив, что его амплитуда превышает нулевой потенциал (явление «овершута»).

После Второй мировой войны, в 1947–1952 годах, Ходжкин и Хаксли совместно с Бернардом Кацем провели серию экспериментов по фиксации потенциала (voltage clamp). Этот метод, разработанный Кеннетом Коулом и Джорджем Мармонтом, позволял удерживать мембранный потенциал на заданном уровне и измерять проходящие через мембрану ионные токи. Ходжкин и Хаксли модифицировали установку для работы с аксоном кальмара и измерили зависимость натриевого и калиевого токов от времени и напряжения.

В 1952 году они опубликовали серию из пяти статей в Journal of Physiology, где представили математическую модель, описывающую свои экспериментальные данные. За эту работу в 1963 году Ходжкин, Хаксли и Джон Экклс (за исследования синаптической передачи) получили Нобелевскую премию по физиологии и медицине.

Основные принципы модели

Модель Ходжкина — Хаксли (HH-модель) описывает электрические свойства мембраны нейрона как эквивалентную электрическую цепь. Мембрана рассматривается как конденсатор ёмкостью C, а ионные каналы — как нелинейные проводимости (g), включённые последовательно с батареями (ионными равновесными потенциалами E).

Эквивалентная схема

В модели выделяются три основных ионных тока:

• Натриевый ток (Iₙₐ): переносится ионами Na⁺. Его проводимость gₙₐ зависит от напряжения и времени.
• Калиевый ток (Iₖ): переносится ионами K⁺. Его проводимость gₖ также зависит от напряжения и времени.
• Ток утечки (Iₗ): переносится в основном ионами Cl⁻ и другими ионами. Его проводимость gₗ считается постоянной.

Общий мембранный ток I равен сумме ёмкостного тока и ионных токов:

\[ I = C_m \frac{dV}{dt} + g_{Na}(V - E_{Na}) + g_K(V - E_K) + g_L(V - E_L) \]

где V — мембранный потенциал, C_m — удельная ёмкость мембраны, E_Na, E_K, E_L — равновесные потенциалы для натрия, калия и утечки соответственно.

Кинетика проводимостей

Ключевое открытие Ходжкина и Хаксли состояло в том, что ионные проводимости не являются постоянными, а изменяются во времени в зависимости от мембранного потенциала. Они описали это с помощью безразмерных переменных «ворот» (gating variables).

• Натриевая проводимость описывается тремя переменными активации m и одной переменной инактивации h:

\[ g_{Na} = \bar{g}_{Na} \cdot m^3 \cdot h \] где \(\bar{g}_{Na}\) — максимальная натриевая проводимость. Переменная m отвечает за активацию (открытие) канала, h — за инактивацию (закрытие). Степень 3 отражает, что для открытия канала должны одновременно активироваться три «ворота».

• Калиевая проводимость описывается четырьмя переменными активации n:

\[ g_K = \bar{g}_K \cdot n^4 \] где \(\bar{g}_K\) — максимальная калиевая проводимость. Степень 4 отражает кооперативность активации.

Каждая переменная ворот (m, h, n) подчиняется кинетическому уравнению первого порядка:

\[ \frac{dx}{dt} = \alpha_x(V)(1 - x) - \beta_x(V)x \]

где x — переменная (m, h или n), α_x(V) и β_x(V) — скоростные константы, зависящие от напряжения. Ходжкин и Хаксли эмпирически подобрали функции α и β, которые наилучшим образом описывали их экспериментальные данные при температуре 6,3 °C.

Уравнения модели

Полная система уравнений HH-модели для гигантского аксона кальмара выглядит следующим образом:

1. Уравнение мембранного потенциала:

\[ C_m \frac{dV}{dt} = I_{ext} - \bar{g}_{Na} m^3 h (V - E_{Na}) - \bar{g}_K n^4 (V - E_K) - g_L (V - E_L) \]

1. Кинетические уравнения:

\[ \frac{dm}{dt} = \alpha_m(V)(1 - m) - \beta_m(V)m \] \[ \frac{dh}{dt} = \alpha_h(V)(1 - h) - \beta_h(V)h \] \[ \frac{dn}{dt} = \alpha_n(V)(1 - n) - \beta_n(V)n \]

1. Функции α и β (для V в мВ, при 6,3 °C):

\[ \alpha_m(V) = \frac{0.1(V + 40)}{1 - \exp(-(V + 40)/10)} \] \[ \beta_m(V) = 4 \exp(-(V + 65)/18) \] \[ \alpha_h(V) = 0.07 \exp(-(V + 65)/20) \] \[ \beta_h(V) = \frac{1}{1 + \exp(-(V + 35)/10)} \] \[ \alpha_n(V) = \frac{0.01(V + 55)}{1 - \exp(-(V + 55)/10)} \] \[ \beta_n(V) = 0.125 \exp(-(V + 65)/80) \]

Параметры модели (типичные для аксона кальмара): C_m = 1 мкФ/см², \(\bar{g}_{Na}\) = 120 мСм/см², \(\bar{g}_K\) = 36 мСм/см², g_L = 0.3 мСм/см², E_Na = +50 мВ, E_K = -77 мВ, E_L = -54.4 мВ.

Биофизическая интерпретация

Модель Ходжкина — Хаксли позволила объяснить механизм генерации потенциала действия на молекулярном уровне задолго до того, как ионные каналы были идентифицированы структурно.

1. Покой: При потенциале покоя (около -65 мВ) калиевые каналы частично открыты, а натриевые каналы закрыты (инактивированы). Мембрана проницаема в основном для K⁺.

1. Деполяризация: Внешний стимул вызывает небольшое смещение потенциала в положительную сторону. Это приводит к открытию натриевых каналов (активация m). Входящий ток Na⁺ дополнительно деполяризует мембрану, вызывая лавинообразное открытие всё большего числа натриевых каналов (регенеративный процесс).

1. Пик потенциала действия: Потенциал мембраны быстро приближается к равновесному потенциалу натрия (+50 мВ). Однако натриевые каналы начинают инактивироваться (h уменьшается), а калиевые каналы — медленно открываться (n увеличивается). Выходящий ток K⁺ начинает реполяризовывать мембрану.

1. Реполяризация: Инактивация натриевых каналов и высокая калиевая проводимость возвращают потенциал к отрицательным значениям. Натриевые каналы остаются инактивированными в течение некоторого времени (рефрактерный период).

1. Следовая гиперполяризация: Из-за того, что калиевые каналы закрываются медленно, калиевая проводимость некоторое время остаётся выше уровня покоя, что приводит к небольшому снижению потенциала ниже уровня покоя.

Значение и применение

Фундаментальное значение

• Первая количественная модель возбудимости: HH-модель показала, что сложные биологические явления могут быть описаны точными математическими уравнениями.
• Подтверждение ионной теории: Модель доказала, что потенциал действия генерируется за счёт изменений ионной проводимости, а не «пробоя» мембраны.
• Основа нейробиологии: Все последующие модели нейронов (интегративно-разрядные, модели с двумя компартментами и т.д.) так или иначе основаны на принципах HH-модели.

Применение в науке

• Нейробиология: Модель используется для изучения механизмов возбудимости, синаптической передачи, пластичности.
• Кардиология: Модификации HH-модели применяются для описания потенциалов действия в кардиомиоцитах и изучения аритмий.
• Нейропротезирование: HH-модель лежит в основе алгоритмов для нейроинтерфейсов и искусственных нейронов.
• Вычислительная нейробиология: Модель является стандартным тестом для численных методов и симуляторов (NEURON, GENESIS).

Ограничения

• Температурная зависимость: Уравнения получены для фиксированной температуры (6,3 °C). Для других температур требуется введение поправочных коэффициентов Q10.
• Специфичность: Модель точно описывает только гигантский аксон кальмара. Для других нейронов (например, кортикальных пирамидных клеток) требуются другие наборы параметров и типов каналов.
• Отсутствие пространственной структуры: Базовая HH-модель описывает точечный нейрон (однородный потенциал). Для учёта распространения сигнала по аксону требуется добавление кабельных уравнений.
• Феноменологический характер: Переменные m, h, n являются эмпирическими подгоночными параметрами, а не прямыми отражениями конформационных изменений белков.

Интересные факты

• Ходжкин и Хаксли проводили эксперименты в Плимуте (Великобритания) на аксонах кальмаров, выловленных в Ла-Манше. Для введения микроэлектродов они использовали тонкие серебряные проволочки.
• Расчёты модели выполнялись вручную с помощью механического арифмометра. Каждый прогон модели занимал несколько недель.
• В 1960 году Ричард Фицхью упростил HH-модель до двух уравнений (модель Фицхью — Нагумо), сохранив её ключевые динамические свойства.
• Современные компьютерные симуляции HH-модели выполняются за миллисекунды, что позволяет моделировать сети из тысяч нейронов.

Источники

• Hodgkin, A. L., & Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology, 117(4), 500–544.
• Hodgkin, A. L., Huxley, A. F., & Katz, B. (1952). Measurement of current-voltage relations in the membrane of the giant axon of Loligo. The Journal of Physiology, 116(4), 424–448.
• Hille, B. (2001). Ion Channels of Excitable Membranes (3rd ed.). Sinauer Associates.
• Johnston, D., & Wu, S. M. (1995). Foundations of Cellular Neurophysiology. MIT Press.
• Koch, C. (1999). Biophysics of Computation: Information Processing in Single Neurons. Oxford University Press.