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

Основные компоненты моделей типа Ходжкина-Хаксли, которые представляют биофизические характеристики клеточных мембран . Липидный бислой представлен как емкость ( См ₎ . Потенциал-управляемые и ионные каналы утечки представлены нелинейной ( g _n ) и линейной ( g _L ) проводимостями соответственно. Электрохимические градиенты, управляющие потоком ионов, представлены батареями (Е), а ионные насосы и обменники представлены источниками тока ( I _p ).

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

Алан Ходжкин и Эндрю Хаксли описали эту модель в 1952 году, чтобы объяснить ионные механизмы, лежащие в основе инициирования и распространения потенциалов действия в гигантском аксоне кальмара . ^{[ 1 ]} За эту работу они получили Нобелевскую премию по физиологии и медицине 1963 года .

Основные компоненты

Типичная модель Ходжкина – Хаксли рассматривает каждый компонент возбудимой клетки как электрический элемент (как показано на рисунке). Липидный бислой представлен как емкость (C _m ). Потенциал-управляемые ионные каналы представлены электропроводностью ( g _n , где n — конкретный ионный канал), которая зависит как от напряжения, так и от времени. Каналы утечки представлены линейными проводимостями ( g _L ). Электрохимические градиенты, управляющие потоком ионов, представлены источниками напряжения ( En _{определяются} ), напряжения которых соотношением внутри- и внеклеточных концентраций интересующих видов ионов. Наконец, ионные насосы представлены источниками тока ( I _p ). ^{[ нужны разъяснения ]} Мембранный потенциал обозначается V _m .

Математически ток, текущий через липидный бислой, записывается как

I_{c}=C_{m}{\frac {{\mathrm {d} }V_{m}}{{\mathrm {d} }t}}

а ток через данный ионный канал является произведением проводимости этого канала и движущего потенциала для конкретного иона.

I_{i}={g_{i}}(V_{m}-V_{i})\;

где $V_{i}$ – реверсивный потенциал конкретного ионного канала. Таким образом, для клетки с натриевыми и калиевыми каналами полный ток через мембрану определяется выражением:

I=C_{m}{\frac {{\mathrm {d} }V_{m}}{{\mathrm {d} }t}}+g_{K}(V_{m}-V_{K})+g_{Na}(V_{m}-V_{Na})+g_{l}(V_{m}-V_{l})

где I - общий ток мембраны на единицу площади, C _m - емкость мембраны на единицу площади, g _K и g _Na - калиевая и натриевая проводимости на единицу площади соответственно, V _K и V _Na - потенциалы обращения калия и натрия. соответственно, а g _l и V _l — проводимость утечки на единицу площади и потенциал обращения утечки соответственно. Зависящими от времени элементами этого уравнения являются V _m , g _Na и g _K , где последние две проводимости явно зависят от мембранного напряжения ( V _m также ).

Характеристика ионного тока

В потенциалзависимых ионных каналах проводимость канала является функцией как времени, так и напряжения ( $g_{n}(t,V)$ на рисунке), находясь в каналах утечки, $g_{l}$ , это константа ( $g_{L}$ на рисунке). Ток, генерируемый ионными насосами, зависит от вида ионов, специфичного для этого насоса. В следующих разделах эти составы будут описаны более подробно.

Потенциал-управляемые ионные каналы

Используя серию экспериментов с фиксированием напряжения и варьируя внеклеточные концентрации натрия и калия, Ходжкин и Хаксли разработали модель, в которой свойства возбудимой клетки описываются набором четырех обыкновенных дифференциальных уравнений . ^{[ 1 ]} Вместе с упомянутым выше уравнением для полного тока это:

I=C_{m}{\frac {{\mathrm {d} }V_{m}}{{\mathrm {d} }t}}+{\bar {g}}_{\text{K}}n^{4}(V_{m}-V_{K})+{\bar {g}}_{\text{Na}}m^{3}h(V_{m}-V_{Na})+{\bar {g}}_{l}(V_{m}-V_{l}),

{\frac {dn}{dt}}=\alpha _{n}(V_{m})(1-n)-\beta _{n}(V_{m})n

{\frac {dm}{dt}}=\alpha _{m}(V_{m})(1-m)-\beta _{m}(V_{m})m

{\frac {dh}{dt}}=\alpha _{h}(V_{m})(1-h)-\beta _{h}(V_{m})h

где I - ток на единицу площади и $\alpha _{i}$ и $\beta _{i}$ — константы скорости i -го ионного канала, зависящие от напряжения, но не от времени. ${\bar {g}}_{n}$ – максимальное значение проводимости. n , m и h представляют собой безразмерные вероятности от 0 до 1, которые связаны с активацией субъединицы калиевого канала , активацией субъединицы натриевого канала и инактивацией субъединицы натриевого канала соответственно. Например, учитывая, что калиевые каналы в гигантском аксоне кальмара состоят из четырех субъединиц, все из которых должны находиться в открытом состоянии, чтобы канал мог проходить ионы калия, n необходимо возвести в четвертую степень. Для $p=(n,m,h)$ , $\alpha _{p}$ и $\beta _{p}$ принять форму

\alpha _{p}(V_{m})=p_{\infty }(V_{m})/\tau _{p}

\beta _{p}(V_{m})=(1-p_{\infty }(V_{m}))/\tau _{p}.

$p_{\infty }$ и $(1-p_{\infty })$ представляют собой стационарные значения активации и инактивации соответственно и обычно представляются уравнениями Больцмана как функции $V_{m}$ . В оригинальной статье Ходжкина и Хаксли ^{[ 1 ]} функции $\alpha$ и $\beta$ даны

{\begin{array}{lll}\alpha _{n}(V_{m})={\frac {0.01(10-V)}{\exp {\big (}{\frac {10-V}{10}}{\big )}-1}}&\alpha _{m}(V_{m})={\frac {0.1(25-V)}{\exp {\big (}{\frac {25-V}{10}}{\big )}-1}}&\alpha _{h}(V_{m})=0.07\exp {\bigg (}-{\frac {V}{20}}{\bigg )}\\\beta _{n}(V_{m})=0.125\exp {\bigg (}-{\frac {V}{80}}{\bigg )}&\beta _{m}(V_{m})=4\exp {\bigg (}-{\frac {V}{18}}{\bigg )}&\beta _{h}(V_{m})={\frac {1}{\exp {\big (}{\frac {30-V}{10}}{\big )}+1}}\end{array}}

где $V=V_{rest}-V_{m}$ обозначает отрицательную деполяризацию в мВ.

Во многих современных программах ^{[ 2 ]} Модели типа Ходжкина – Хаксли обобщают $\alpha$ и $\beta$ к

{\frac {A_{p}(V_{m}-B_{p})}{\exp {\big (}{\frac {V_{m}-B_{p}}{C_{p}}}{\big )}-D_{p}}}

Чтобы охарактеризовать потенциал-управляемые каналы, уравнения можно адаптировать к данным фиксации напряжения. О выводе уравнений Ходжкина – Хаксли при ограничении напряжения см. ^{[ 3 ]} Вкратце, когда мембранный потенциал поддерживается постоянным значением (т. е. с помощью фиксатора напряжения), для каждого значения мембранного потенциала нелинейные уравнения стробирования сводятся к уравнениям вида:

m(t)=m_{0}-[(m_{0}-m_{\infty })(1-e^{-t/\tau _{m}})]\,

h(t)=h_{0}-[(h_{0}-h_{\infty })(1-e^{-t/\tau _{h}})]\,

n(t)=n_{0}-[(n_{0}-n_{\infty })(1-e^{-t/\tau _{n}})]\,

Таким образом, для любого значения мембранного потенциала $V_{m}$ токи натрия и калия можно описать формулой

I_{\mathrm {Na} }(t)={\bar {g}}_{\mathrm {Na} }m(V_{m})^{3}h(V_{m})(V_{m}-E_{\mathrm {Na} }),

I_{\mathrm {K} }(t)={\bar {g}}_{\mathrm {K} }n(V_{m})^{4}(V_{m}-E_{\mathrm {K} }).

Чтобы прийти к полному решению для распространяющегося потенциала действия, необходимо записать текущий член I в левой части первого дифференциального уравнения через V , чтобы уравнение стало уравнением только для напряжения. Связь между I и V может быть получена из теории кабеля и определяется выражением

I={\frac {a}{2R}}{\frac {\partial ^{2}V}{\partial x^{2}}},

где а радиус аксона , R — удельное сопротивление аксоплазмы — , х — положение вдоль нервного волокна. Замена этого выражения вместо I преобразует исходную систему уравнений в систему дифференциальных уравнений в частных производных , поскольку напряжение становится функцией как x, так и t .

Алгоритм Левенберга – Марквардта часто используется для согласования этих уравнений с данными фиксации напряжения. ^{[ 4 ]}

Хотя в первоначальных экспериментах использовались только натриевые и калиевые каналы, модель Ходжкина-Хаксли также можно расширить для учета других видов ионных каналов .

Каналы утечки

Каналы утечки учитывают естественную проницаемость мембраны для ионов и принимают форму уравнения для потенциалзависимых каналов, где проводимость $g_{leak}$ является константой. Таким образом, ток утечки из-за пассивных ионных каналов утечки в формализме Ходжкина-Хаксли равен $I_{l}=g_{leak}(V-V_{leak})$ .

Насосы и теплообменники

Мембранный потенциал зависит от поддержания на ней градиентов концентрации ионов. Поддержание этих градиентов концентрации требует активного транспорта ионных частиц. натрий -калий и натрий-кальций Наиболее известны из них . Некоторые основные свойства Na/Ca-обменника уже хорошо известны: стехиометрия обмена составляет 3 Na. ⁺: 1 Калифорния ²⁺ а теплообменник электрогенен и чувствителен к напряжению. Обменник Na/K также был подробно описан с 3 Na ⁺: 2 К ⁺ стехиометрия. ^{[ 5 ]}^{[ 6 ]}

Математические свойства

Модель Ходжкина-Хаксли можно рассматривать как систему дифференциальных уравнений с четырьмя переменными состояния : $V_{m}(t),n(t),m(t)$ , и $h(t)$ , которые меняются со временем $t$ . Систему сложно изучать, поскольку она является нелинейной системой , не может быть решена аналитически и, следовательно, не имеет в замкнутой форме решения . Однако существует множество численных методов анализа системы. определенных свойств и общего поведения, таких как предельные циклы Можно доказать существование .

Центральный коллектор

Поскольку существует четыре переменных состояния, визуализация пути в фазовом пространстве может быть затруднена. Обычно выбираются две переменные: напряжение $V_{m}(t)$ и переменная гейтирования калия $n(t)$ , позволяющий визуализировать предельный цикл . Однако следует быть осторожным, поскольку это специальный метод визуализации четырехмерной системы. Это не доказывает существования предельного цикла.

Лучший прогноз можно построить на основе тщательного анализа якобиана системы , оцененного в точке равновесия . В частности, собственные значения якобиана указывают на существование центрального многообразия . Аналогично, собственные векторы якобиана показывают ориентацию центрального многообразия . Модель Ходжкина – Хаксли имеет два отрицательных собственных значения и два комплексных собственных значения со слегка положительными действительными частями. Собственные векторы, связанные с двумя отрицательными собственными значениями, уменьшатся до нуля с увеличением времени t . Остальные два комплексных собственных вектора определяют центральное многообразие. Другими словами, четырехмерная система коллапсирует на двухмерную плоскость. Любое решение, начинающееся с центрального многообразия, будет затухать по направлению к центральному многообразию. Более того, предельный цикл содержится в центральном многообразии.

Напряжение v ( t ) (в милливольтах) модели Ходжкина – Хаксли, построенное на графике за 50 миллисекунд. Инжектируемый ток варьируется от -5 наноампер до 12 наноампер. График проходит три стадии: стадию равновесия, стадию одного всплеска и стадию предельного цикла.

Бифуркации

Если подаваемый ток $I$ использовались в качестве параметра бифуркации , то модель Ходжкина–Хаксли претерпевает бифуркацию Хопфа . Как и в большинстве нейронных моделей, увеличение подаваемого тока увеличит скорость срабатывания нейрона. Одним из последствий бифуркации Хопфа является наличие минимальной скорострельности. Это означает, что либо нейрон вообще не срабатывает (что соответствует нулевой частоте), либо срабатывает с минимальной частотой. Из-за принципа «все или ничего не происходит » плавного увеличения амплитуды потенциала действия , а происходит резкий «скачок» амплитуды. Получающийся переход известен как утка .

Улучшения и альтернативные модели

Модель Ходжкина-Хаксли считается одним из величайших достижений биофизики XX века. Тем не менее, современные модели типа Ходжкина – Хаксли были расширены несколькими важными способами:

Дополнительные популяции ионных каналов были включены на основе экспериментальных данных.
Модель Ходжкина-Хаксли была модифицирована, чтобы включить теорию переходного состояния и создать термодинамические модели Ходжкина-Хаксли. ^{[ 7 ]}
Модели часто включают в себя очень сложную геометрию дендритов и аксонов , часто основанную на данных микроскопии.
Модели, основанные на проводимости, подобные модели Ходжкина-Хаксли, включают знания о типах клеток, определяемых транскриптомикой отдельных клеток. ^{[ 8 ]}
Стохастические модели поведения ионных каналов, приводящие к стохастическим гибридным системам. ^{[ 9 ]}
Модель Пуассона-Нернста-Планка (PNP) основана на приближении среднего поля ионных взаимодействий и континуальном описании концентрации и электростатического потенциала. ^{[ 10 ]}

Также было разработано несколько упрощенных нейронных моделей (таких как модель ФитцХью-Нагумо ), способствующих эффективному крупномасштабному моделированию групп нейронов, а также математическому пониманию динамики генерации потенциала действия.

См. также

Ссылки

^ Перейти обратно: ^а ^б ^с Ходжкин А.Л., Хаксли А.Ф. (август 1952 г.). «Количественное описание мембранного тока и его применение к проводимости и возбуждению нерва» . Журнал физиологии . 117 (4): 500–44. doi : 10.1113/jphysicalol.1952.sp004764 . ПМК 1392413 . ПМИД 12991237 .
^ Нельсон М.Э. (2005) Электрофизиологические модели в: База данных мозга: от данных к знаниям. (С. Козлоу и С. Субраманиам, ред.) Wiley, Нью-Йорк, стр. 285–301.
^ Грей диджей, Ву СМ (1997). Основы клеточной нейрофизиологии (3-е изд.). Кембридж, Массачусетс [ua]: MIT Press. ISBN 978-0-262-10053-3 .
^ Крапивин Владимир Федорович; Вароцос, Костас А.; Солдатов Владимир Юрьевич. (2015). Новые инструменты экоинформатики в науке об окружающей среде: приложения и принятие решений . Спрингер. стр. 37–38. ISBN 9783319139784 .
^ Раковски Р.Ф., Гадсби, округ Колумбия, Де Вир П. (май 1989 г.). «Стехиометрия и зависимость напряжения натриевого насоса в зажатом по напряжению, подвергнутом внутреннему диализу гигантском аксоне кальмара» . Журнал общей физиологии . 93 (5): 903–41. дои : 10.1085/jgp.93.5.903 . ПМК 2216238 . ПМИД 2544655 .
^ Хилле Б (2001). Ионные каналы возбудимых мембран (3-е изд.). Сандерленд, Массачусетс: Синауэр. ISBN 978-0-87893-321-1 .
^ Форрест, доктор медицины (май 2014 г.). «Может ли термодинамическая модель зависимой от напряжения проводимости экстраполироваться на температуру?» (PDF) . Вычисление . 2 (2): 47–60. дои : 10.3390/computation2020047 .
^ Нанди, Анирбан; Шартран, Томас; Ван Гейт, Вернер; Бучин, Анатолий; Яо, Цзычжэнь; Ли, Су Юн; Вэй, Ина; Калмбах, Брайан; Ли, Брайан; Лейн, Эд; Берг, Джим; Сюмбюль, Уйгар; Кох, Кристоф; Ташич, Босилька; Анастасиу, Костас А. (9 августа 2022 г.). «Модели одного нейрона, связывающие электрофизиологию, морфологию и транскриптомику типов клеток коры» . Отчеты по ячейкам . 40 (6): 111176. doi : 10.1016/j.celrep.2022.111176 . ISSN 2211-1247 . ПМЦ 9793758 . ПМИД 35947954 . S2CID 215790820 .
^ Пакдаман, К.; Тиуллен, М.; Вайнриб, Г. (2010). «Предельные теоремы о жидкости для стохастических гибридных систем с применением к моделям нейронов». Адв. Прил. Вероятно . 42 (3): 761–794. arXiv : 1001.2474 . Бибкод : 2010arXiv1001.2474P . дои : 10.1239/aap/1282924062 . S2CID 18894661 .
^ Чжэн, К.; Вэй, GW (май 2011 г.). «Модель Пуассона-Больцмана-Нернста-Планка» . Журнал химической физики . 134 (19): 194101. Бибкод : 2011JChPh.134s4101Z . дои : 10.1063/1.3581031 . ПМК 3122111 . ПМИД 21599038 .

Дальнейшее чтение

Ходжкин А.Л., Хаксли А.Ф. (апрель 1952 г.). «Токи, переносимые ионами натрия и калия через мембрану гигантского аксона Лолиго» . Журнал физиологии . 116 (4): 449–72. дои : 10.1113/jphysicalol.1952.sp004717 . ПМЦ 1392213 . ПМИД 14946713 .
Ходжкин А.Л., Хаксли А.Ф. (апрель 1952 г.). «Компоненты мембранной проводимости гигантского аксона Лолиго» . Журнал физиологии . 116 (4): 473–96. doi : 10.1113/jphysicalol.1952.sp004718 . ПМК 1392209 . ПМИД 14946714 .
Ходжкин А.Л., Хаксли А.Ф. (апрель 1952 г.). «Двойное влияние мембранного потенциала на натриевую проводимость в гигантском аксоне Лолиго» . Журнал физиологии . 116 (4): 497–506. doi : 10.1113/jphysicalol.1952.sp004719 . ПМЦ 1392212 . ПМИД 14946715 .
Ходжкин А.Л., Хаксли А.Ф. (август 1952 г.). «Количественное описание мембранного тока и его применение к проводимости и возбуждению нерва» . Журнал физиологии . 117 (4): 500–44. doi : 10.1113/jphysicalol.1952.sp004764 . ПМК 1392413 . ПМИД 12991237 .
Ходжкин А.Л., Хаксли А.Ф., Кац Б. (апрель 1952 г.). «Измерение вольт-амперных отношений в мембране гигантского аксона Лолиго» . Журнал физиологии . 116 (4): 424–48. дои : 10.1113/jphysicalol.1952.sp004716 . ПМЦ 1392219 . ПМИД 14946712 .

Внешние ссылки

Интерактивное моделирование модели HH на языке Javascript . Работает в любом браузере с поддержкой HTML5. Позволяет изменять параметры модели и подачу тока.
Интерактивный Java-апплет модели HH. Параметры модели могут быть изменены, а также возможны параметры возбуждения и построение фазовых пространств всех переменных.
Прямая ссылка на модель Ходжкина-Хаксли и описание в базе данных BioModels.
Нейронные импульсы: потенциал действия в действии Гаррет Неске, Демонстрационный проект Wolfram
Интерактивная модель Ходжкина – Хаксли Шимона Марома, Демонстрационный проект Вольфрама
ModelDB База данных исходного кода вычислительной нейронауки, содержащая 4 версии (в разных симуляторах) исходной модели Ходжкина-Хаксли и сотни моделей, которые применяют модель Ходжкина-Хаксли к другим каналам во многих типах электрически возбудимых клеток.
Несколько статей о стохастической версии модели и ее связи с исходной.

[HH-1] Перейти обратно: ^а ^б ^с Ходжкин А.Л., Хаксли А.Ф. (август 1952 г.). «Количественное описание мембранного тока и его применение к проводимости и возбуждению нерва» . Журнал физиологии . 117 (4): 500–44. doi : 10.1113/jphysicalol.1952.sp004764 . ПМК 1392413 . ПМИД 12991237 .

[2] Нельсон М.Э. (2005) Электрофизиологические модели в: База данных мозга: от данных к знаниям. (С. Козлоу и С. Субраманиам, ред.) Wiley, Нью-Йорк, стр. 285–301.

[JohnstonAndWu-3] Грей диджей, Ву СМ (1997). Основы клеточной нейрофизиологии (3-е изд.). Кембридж, Массачусетс [ua]: MIT Press. ISBN 978-0-262-10053-3 .

[4] Крапивин Владимир Федорович; Вароцос, Костас А.; Солдатов Владимир Юрьевич. (2015). Новые инструменты экоинформатики в науке об окружающей среде: приложения и принятие решений . Спрингер. стр. 37–38. ISBN 9783319139784 .

[Rakowski_1989-5] Раковски Р.Ф., Гадсби, округ Колумбия, Де Вир П. (май 1989 г.). «Стехиометрия и зависимость напряжения натриевого насоса в зажатом по напряжению, подвергнутом внутреннему диализу гигантском аксоне кальмара» . Журнал общей физиологии . 93 (5): 903–41. дои : 10.1085/jgp.93.5.903 . ПМК 2216238 . ПМИД 2544655 .

[Hille-6] Хилле Б (2001). Ионные каналы возбудимых мембран (3-е изд.). Сандерленд, Массачусетс: Синауэр. ISBN 978-0-87893-321-1 .

[7] Форрест, доктор медицины (май 2014 г.). «Может ли термодинамическая модель зависимой от напряжения проводимости экстраполироваться на температуру?» (PDF) . Вычисление . 2 (2): 47–60. дои : 10.3390/computation2020047 .

[8] Нанди, Анирбан; Шартран, Томас; Ван Гейт, Вернер; Бучин, Анатолий; Яо, Цзычжэнь; Ли, Су Юн; Вэй, Ина; Калмбах, Брайан; Ли, Брайан; Лейн, Эд; Берг, Джим; Сюмбюль, Уйгар; Кох, Кристоф; Ташич, Босилька; Анастасиу, Костас А. (9 августа 2022 г.). «Модели одного нейрона, связывающие электрофизиологию, морфологию и транскриптомику типов клеток коры» . Отчеты по ячейкам . 40 (6): 111176. doi : 10.1016/j.celrep.2022.111176 . ISSN 2211-1247 . ПМЦ 9793758 . ПМИД 35947954 . S2CID 215790820 .

[stochastic-9] Пакдаман, К.; Тиуллен, М.; Вайнриб, Г. (2010). «Предельные теоремы о жидкости для стохастических гибридных систем с применением к моделям нейронов». Адв. Прил. Вероятно . 42 (3): 761–794. arXiv : 1001.2474 . Бибкод : 2010arXiv1001.2474P . дои : 10.1239/aap/1282924062 . S2CID 18894661 .

[10] Чжэн, К.; Вэй, GW (май 2011 г.). «Модель Пуассона-Больцмана-Нернста-Планка» . Журнал химической физики . 134 (19): 194101. Бибкод : 2011JChPh.134s4101Z . дои : 10.1063/1.3581031 . ПМК 3122111 . ПМИД 21599038 .

[ 1 ]

[ 2 ]

[ 3 ]

[ 4 ]

[ 5 ]

[ 6 ]

[ 7 ]

[ 8 ]

[ 9 ]

[ 10 ]