CC BY
14
КУТЫРКИНА М. А., МАМЕДОВА Т. Ф.
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ НЕЙРОННОЙ СЕТИ НА ОСНОВЕ МОДЕЛИ ХОДЖКИНА-ХАКСЛИ
Аннотация. Статья посвящена численному моделированию колебаний мембранного потенциала и воротных переменных в модели Ходжкина-Хаксли. Особое внимание уделяется исследованию изменения входного тока на нейронах, от которого зависит способность нейронов к обучению и запоминанию информации.
Ключевые слова: математическая модель Ходжкина-Хаксли, нейронная сеть, мембранная емкость, воротные переменные.
KUTYRKINA M. A., MAMEDOVA T. F.
NUMERICAL SIMULATION OF NEURAL NETWORK BASED ON HODGKIN-HUXLEY MODEL
Abstract. The article deals with numerical simulation of membrane potential fluctuations and gate variables in the Hodgkin-Huxley model. Particular attention is paid to the study of changes in the input current on neurons, which determines the ability of neurons to learn and remember information.
Keywords: Hodgkin-Huxley mathematical model, neural network, membrane capacitance, gate variables.
Актуальность проблемы. Современная неврология предъявляет высокие требования к квалификации врачей и методам лечения. Внедрение в неврологию математических методов, в частности методов математического моделирования, позволяет существенно расширить прогнозирование, диагностику и терапию заболеваний.
Нейрон — основной функциональный элемент нервной системы. Типичный нейрон состоит из трех частей: тела, дендритов и аксона. Дендрит принимает сигналы от другого нейрона, а аксон передаёт их другим нейронам. Место контакта аксона одного нейрона и дендрита другого нейрона называется синапсом. Сигналы, принимаемые с дендритов, суммируются в теле нейрона и если данная сумма превышает определённый порог, то возникает нервный импульс или спайк [1]. Кроме возбуждения нейроны могут и тормозить друг друга, уменьшая "сумму" внутри нейрона и препятствуя возникновению спайка.
Оболочка нейрона является хорошим изолятором, поэтому ионные составы вне и внутри нейрона различны. Заряженные частицы переносятся между внутренностью нейрона и окружающей его межклеточной жидкостью с помощью так называемых ионных насосов (существует, например, натрий-калиевый насос). Импульсы от возбуждающих и тормозящих
1
нейронов, приходящие к нашему нейрону, "включают" разные виды ионных насосов. Так регулируется электрическое поле нейрона. В зависимости от того, каким оно будет, импульс либо побежит дальше, либо нет.
Постановка задачи. На данный момент созданы различные как относительно простые модели, например, «Inregrate and Fire», в которой нейрон представляется в виде конденсатора и резистора, так и более детальные, биологически правдоподобные, модели, например, модель Ходжкина-Хаксли, которая гораздо сложнее как в вычислительном плане, так и в плане анализа её динамики, но она гораздо точнее описывает динамику мембранного потенциала нейрона [2].
Классическая модель Ходжкина-Хаксли описывается следующей системой дифференциальных уравнений [3]:
dV _
Gm ^ _ Iext + Isynif) — dm
- GK- Gl,
— =a(V)(1-m)-(3m(V)rn,
^ = ah(V)(l-h)-ph(V)h, dn
— = an(V)(1-n)-(3n(V)n,
где
am(V) =
0.1(25-V) exp(25-V)~
10
Pm(V) = 4exp(-^); ah(V) = 0.07exp (-^) ;
Ph(V) =
1
1 ^ f30 - ^ ' 1 +exp
an(V) =
0.01(10 — V)
77 (10-vy
1 +exp
/3n(V) = 0.125 exp
GNa = dNam3h(V - VNa); GK=gKn4(V-VK); GL=gL(V-VL).
Начальные условия для системы (1) имеют вид
7(0) = Vo;m(0) = то; h(0) = ho;n(0) = щ.
V.
1
Здесь V(t) - мембранный потенциал, m(t), h(t), n(t) - так называемые воротные переменные: т, п - активационные переменные для калиевых и натриевых каналов, h -инактивационная переменная натриевых каналов. Эти переменные задают изменения проводимости ионных каналов в зависимости от мембранного потенциала, Ст - мембранная емкость. Остальные параметры были получены экспериментальным путем и имеют следующие значения:
gNa = 120 мСм; ^ = 36 мСм,; gL = 0,3 мСм; VNa = 115 мВ; VK =-12 мВ; VL = 10,6 мВ.
Сумма синаптических токов рассчитывается по следующей формуле:
Ne Ni
hyn = Cmdexc[^i ^ l(t — tIi) — ^ ^^^^ — tin)]. k=1 I m=1 n
Здесь Ne - количество возбуждающих нейронов, Ni - количество тормозных нейронов, 1-(t — t[) и l(t — t%.) - индикаторные функции, принимающие значение 1, если их аргументы равны нулю, и 0 - в противном случае, t!k - время разряда 1-го спайка к-го нейрона, Ст — мембранная емкость. Константы, дехс = 1 и К = 4 вычислены эмпирически.
Численные расчеты получены методом Рунге-Кутты 4-го порядка при t Е [0, 300] и следующих параметрах: количество возбуждающих нейронов - 200, количество тормозящих нейронов - 50.
Данная модель была реализованы в среде разработки Google Colaboratory. При программировании использовались библиотеки NumPy и SciPy, а также библиотека Matplotlib, позволяющая визуализировать полученные данные.
Рассмотрим случай, когда значение входного тока lext принимает значение, равное
5мА.
На рис. 1 представлен график колебаний мембранного потенциала с течением времени, на рис. 2 - график изменения воротных переменных с течением времени, и на рис. 3 -зависимость мембранного потенциала V от времени.
Исходя из данных графиков, можно сделать вывод, что при значении входного тока lext = 5 мА колебания мембранного потенциала исчезают. То же происходит и с воротными переменными. Из рис. 3 видно, что траектория зависимости воротных переменных от мембранного потенциала V не является периодической.
Рис. 1. Зависимость мембранного потенциала V от времени по модели Ходжкина-Хаксли
при значении входного тока 1еХ{ = 5 мА.
Рис. 3. Зависимость воротных переменных от мембранного потенциала V при значении входного тока 1ех1 = 5 мА.
С другой стороны, если входной ток достигает значения 5 мА, то предельный цикл фактически вырождается. Это значит, что такое значение входного тока слишком мало, чтобы возбудить нейрон.
Увеличим значение входного тока до 10 мА. По графику колебаний мембранного потенциала, представленного на рис. 4, и графику воротных изменения воротных переменных (см. рис. 5), можно предположить, что процесс является периодическим. После 1-2 пиков потенциал стабилизируется. Аналогичным образом ведут себя воротные переменные (см. рис. 5). Можно предположить (см. рис. 6), что траектория зависимости воротных переменных от мембранного потенциала V представляет собой невырожденный предельный цикл.
Судя по результатам численного эксперимента, при значении входного тока в 10 мА имеем наиболее ярко выраженный предельный цикл. Данное значение является оптимальным для возбуждения нейрона и передачи сигналов от одного нейрона другим.
Мембранный потенциал
20
0
>
3 50 100 150 200 250 300
t (тэ)
Рис. 4. Зависимость мембранного потенциала от времени по модели Ходжкина-Хаксли при значении входного тока 1ех1: = 10 мА.
О 50 100 150 200 250 300
Мтэ)
Мембранный потенциал (V)
Рис. 6. Зависимость воротных переменных от мембранного потенциала V при значении входного тока 1еХ{ = 10 мА.
Пусть теперь значение входного тока достигает значения 1ех1 = 50 мА. Если мы посмотрим на графики колебаний мембранного потенциала (рис. 7) и изменения воротных переменных, то увидим, что после первоначального скачка амплитуда колебания мембранного потенциала в периодическом режиме снижается. Это же происходит и с воротными переменными см. (рис. 8). Из рис. 9 следует, что выход на предельный цикл теперь происходит медленнее (траектория делает большее число витков).
При таком значении тока для достижения предельного цикла нужно сделать большее число витков по фазовой траектории. Такое значение входного тока губительно для нейронов головного мозга.
Рис. 7. Зависимость мембранного потенциала от времени по модели Ходжкина-Хаксли при значении входного тока 1ех1: = 50 мА.
■ее -40 -20 О 20 40
Мембранный потенциал (V)
Рис. 9. Зависимость воротных переменных от мембранного потенциала V при значении входного тока lext = 50 мА.
Заключение. Таким образом, меняя значения входного тока в модели Ходжкина-Хаксли, можно сделать вывод, что наиболее оптимальным значением, при котором возможно возбуждение нейрона и передача сигналов другим нейронам, является значение входного тока, равное 10 мА.
СПИСОК ЛИТЕРАТУРЫ
1. Никифоров А. С., Гусев Е. И. Общая неврология: учебное пособие. - М.: ГЭОТАР-Медиа, 2007. - 720 с.
2. Шульговский В. В. Основы нейрофизиологии: учебное пособие для студентов вузов. - М.: Аспект Пресс, 2000. - 277 с.
3. Hodgkin A. L., Huxley A. F. Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo // Journal of Physiology. - 1952. - pp. 449-472.
