Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,593

СОБСТВЕННЫЙ ПОТЕНЦИАЛ ДЕЙСТВИЯ НЕЙРОНА С МОДЕЛЬНОЙ МЕМБРАНОЙ

Галанов Е.К. 1
1 Петербургский университет путей сообщения
Исследован механизм возникновения потенциала действия нейрона, находящегося в термодинамическом равновесии. Показано, что спусковым механизмом возникновения потенциала действия может быть фазовый переход липидной мембраны из жидкокристаллической в гель-фазу под действием электростатического поля, вызванного веерными колебаниями полярных «головок» фосфолипидов. Модельная мембрана взята в виде сферы, построенной из бислоя молекул фосфатидилхолина. Ориентация полярных «головок» молекул относительно направления углеводородной цепи определяет асимметрию ангармонического потенциала осциллятора веерного типа. Представлены волновые функции и энергии первых шести стационарных состояний изолированной молекулы. Рассчитано радиальное отклонение положительного иона диполя хср полярной «головки» фосфатидилхолина, которое определяет изменение направления вектора дипольного момента Δр полярной «головки» при изменении температуры. Показано, что при повышении температуры бислоя мембраны с 310 °К до 315 °К в пределах этого бислоя возникает дополнительное электростатическое поле Δ Е = 104 –107 В/м , которое вызывает фазовый переход из жидкокристаллической фазы в гель-фазу. Этот фазовый переход является спусковым механизмом изменения структуры бислоя, активации ионных каналов, включения ионных насосов, … и завершается возвращением мембраны в жидкокристаллическую фазу. Рассмотрены энергетические уровни и соответствующие энергетические зоны ангармонических осцилляторов веерного типа.
потенциал действия нейрона
оптические фононы мембраны
конформационные состояния фосфолипидов
фазовый переход липидной мембраны
1. Хьюбел Д. Глаз, мозг, зрение / Д. Хьюбел. – М.: Мир, 1990. – 241 с.
2. Альтман Я.А. Слуховая система / Альтман Я.А., И.А. Вартанян. – Л.: Наука, 1990. – 618 с.
3. Сотников О.О. Синцитиальная цитоплазматическая связь и слияние нейронов / О.О. Сотников. – СПб.: Наука, 2013. – 201 с.
4. Рубин А.Б. Биофизика / А.Б. Рубин. – М., Ижевск: Институт компьютерных технологий, 2013. – 381 с.
5. Капустин А.П. Экспериментальные исследования жидких кристаллов / А.П. Капустин. – М.: Мир, 1994. – 475 с.
6. Барфут Дж. Введение в физику сегнетоэлектрических явлений / Дж. Барфут. – М.: Мир, 1970. – 350 с.
7. Камзина Л.С. Влияние электрического поля на размытые фазовые переходы в кристаллах PbInNbO3 – PbTiO3 / Л.С. Камзина, С.П. Попов // Физика твёрдого тела. – 2014. – № 11. – С. 2142–2145.
8. Fridkin V., Ducbarme S. Ferroelectrricity at the nanoscale. Phys. Usp, 2014, v. 57, рр. 597–603.
9. Галанов Е.К. Оптические фононы модельной мембраны нейрона / Е.К. Галанов // Альманах современной науки и образования. – 2017. – № 1. – С. 19–22.
10. Фок В.А. Начала квантовой механики / В.А. Фок. – М.: Издательство ЛКИ, 2007. – 374 с.
11. Галанов Е.К. Потенциал покоя модельной мембраны нейрона / Е.К. Галанов // Альманах современной науки и образования. – 2017. – № 7. – С. 27–31.
12. Стокач Н.С. Электронно-колебательные спектры многоатомных молекул / Н.С. Cтокач, Г.А. Шигорян, Е.А. Гастолович. – М.: Наука, 1982. – 360 с.

Возникновение потенциала действия нейрона (состояние возбуждения нейрона) обычно трактуется как результат воздействия внешнего стимула. Частота следования потенциалов действия (ПД) первых (после рецепторов) нейронов имеет значение 30–5000 имп/с.

Однако, даже в отсутствии внешних стимулов нейроны ряда областей мозга проявляют спонтанную активность. В случае ганглиозных клеток (зрительная система) частота спонтанной активности этих клеток при отсутствии светового сигнала колеблется в пределах 2–20 имп/с. [1].

В нейронах кохлеарного ядра (слуховая система) спонтанная активность (при отсутствии звукового сигнала) имеет значение 2–30 имп/с. [2]. Спонтанная активность нейронов проявляется и когда они находятся в изолированном состоянии, её величина равна 1–2 имп/с. [3]. Во всех приведенных случаях система, содержащая нейрон, находится в состоянии термодинамического равновесия. Что является стимулом (спусковым механизмом) возникновения спонтанного ПД? Учитывая состояние термодинамического равновесия системы, ряд факторов, традиционно связываемых с возникновением ПД, необходимо исключить. Это активизация ионных каналов, включение (выключение) ионных насосов и др. [4].

Первоначальным стимулом (спусковым механизмом) активизации ПД может быть фазовый переход липидной мембраны при условии, что в состоянии термодинамического равновесия она находится в высокотемпературной (жидкокристаллической) фазе и под действием оптических фононов (веерные колебания полярных «головок») переходит в низкотемпературную (гель) фазу. Этот переход осуществляется в результате того, что оптические фононы веерных колебаний полярных «головок» липидов, структурированных в мембране, изменяют напряжённость электростатического поля мембраны.

Как известно [5–7], фазовый переход в неорганических и органических кристаллах из высокотемпературной фазы в низкотемпературную фазу (при постоянной температуре) совершается при приложении электрического поля, если кристалл находится при температуре выше точки фазового перехода на 1–5 °С в случае неорганических кристаллов [6–7] и на 2–10 °С в случае органических кристаллов [5]. Этот переход происходит в том числе, когда отмеченные кристаллы имеют наноразмерные толщины [8]. Напряжённость электрического поля, вызывающего этот фазовый переход, составляет 104–106 В/м. Фазовый переход можно считать спусковым механизмом (первой фазой) ПД. Последующая (во времени) структура ПД связана с активизацией ионных каналов, включением ионных насосов… и, конечно, с возвращением мембраны в термодинамическое равновесие, т.е. в жидкокристаллическую фазу.

В работе [9] рассмотрены оптические фононы модельной мембраны, представляющей собой сферу. Возьмём за основу эту же структуру мембраны нейрона, ориентация молекул в которой определяется двойным липидным слоем (≈70 % массы липидной мембраны).

Липидные молекулы могут находится в двух конформационных состояниях. Как и в [9], считаем, что эти два конформационных состояния определяются двумя основными электронными состояниями (i = 1,2) молекул Φi,j = φi (r, R0)∙ψi,j(Qj), (где r и R – совокупность координат электронов и ядер), вероятность перехода между которыми мала. В адиабатическом приближении энергии этих состояний равны Ei,j = Eiэл (r, R0) + ,Еiкj (Qj), где R0 – равновесная конфигурация ядер, Qj – обобщённые колебательные координаты молекул; ¦Еiэл (r, R0 )¦>> ¦EiкJ (Qj)¦, Eiкj – колебательная энергия.

1. Рассмотрим сначала колебания изолированной молекулы фосфолипида на примере фосфатидилхолина. Молекула находится в конформационном состоянии с более высокой энергией i = 1. На рис. 1 представлена структурная схема этой молекулы. Ориентацию (α0) полярной «головки» относительно гидрофобного «хвоста» фосфатидилхолина возьмём такую же, которую она имеет в мембране α0 ≈ 300. Из всех колебаний изолированной молекулы (3N-6) (N – число атомов в молекуле) рассмотрим веерные колебания полярной «головки» относительно гидрофобного «хвоста». Энергии квантов валентных и деформационных колебаний имеют значения 3500–400 см-1. Энергия квантов веерных колебаний больших фрагментов молекул равна 50–150 см-1. Возьмём значение кванта энергии веерного колебания полярной «головки» фосфатидилхолина равным 100 см-1.

gal1.tif

Рис. 1. Структура молекулы фосфатидилхолина; р – дипольный момент полярной «головки» молекулы

Структурное положение полярной «головки» в молекуле фосфатидилхолина (рис. 1) определяет асимметрию веерных колебаний (асимметрию потенциальной энергии U = f(α)). В работе [9] потенциальная энергия взята в форме U = [α4 + U0 exp[ – 2(α + 1)2]. Именно асимметрия веерных колебаний полярных «головок» должна обуславливать изменение ориентации (α) диполя «головок» при изменении температуры.

Поскольку энергия кванта веерных колебаний hν = 100 cm-1(hν(Дж)/hc) существенно меньше средней энергии теплового движения частиц ɛ = kT = 4,28·10-21 Дж (217 см-1) при Т = 315 °К, необходимо определить волновые функции и энергии колебательных уровней значительного числа возбуждённых состояний. Для их определения возьмём асимметричный ангармонический потенциал в более простой форме в сравнении с потенциалом работы [9], а именно:

U = (1/2·ε2 + β·ε3). (1)

Уравнение Шредингера для стационарных состояний такого ангармонического осциллятора имеет вид

--1/2·d2 ψ/dε2 + (1/2·ε2 + β·ε3)∙ψ = E∙ψ, (2)

где β – малый поправочный член, возьмём β = 0,03; ε = х/x0 (рис. 1), gal01.wmf

gal02.wmf – собственная частота гармонического осциллятора (β = 0);

mпр – приведенная масса полярной «головки»; x = d·sinα, d – плечо диполя полярной «головки».

Собственные функции ангармонического осциллятора (2) в первом приближении определяются выражением [10]:

ψn = ψn0 + β·ψ'n , (3)

gal03.wmf

gal04.wmf

где ψn0 – собственные функции гармонического осциллятора (β = 0).

gal05.wmf (4)

Собственные значения энергий стационарных состояний осциллятора (2) во втором приближении равны

En = En0 + β2 ·E'n = h·ν0 [(n + 1/2) – 15/4·β2∙(n2 + n + 11/30)],, (5)

где n = 0,1, 2,…

Волновые функции и энергии основного и первых возбуждённых колебательных состояний молекулы фосфатидилхолина определяются соотношениями

ψ1,0 = [0,75 – 0,0224 ε – 0,0075 ε3]·exp(–ε2/2),

ψ1,1 = [0,0543 + 1,06 ε – 0,122 ε2 – 0,0089 ε4]·exp(–ε2/2),…

. (6)

.

.

Ψ1,6 = [– 0,48 + 0,318ε + 2,52ε2 – 1,23ε3 – 1,68ε4 +0,75ε5 + 0,224ε6 –

– 0,107 ε7 + 0,00205ε9]∙exp(–ε2/2);

Е1,0 = 49,87 см-1; Е1,1 = 149,20см-1; Е1,2 = 248,19 см-1 ,…Е1,6 = 635,71 см-1.

Величина дипольного момента полярных головок фосфолипидов находится в пределах gal06.wmf Кл•м. Для фосфатидилхолина возьмём величину дипольного момента, принятую в работе [11] ¦p¦ = 0,6•10-9 м•1,6•10-19 Кл = 0,96•10-28 Кл•м. Асимметричность потенциальной энергии рассматриваемого осциллятора (1) обуславливает асимметричность волновых функций, изменение радиальной составляющей дипольного момента полярной «головки» и радиальной координаты (параллельной радиальному направлению мембраны) положительного иона диполя Δr = xср соsα0 (xср – среднее значение положения (рис. 1) иона диполя при изменении температуры):

gal07.wmf

gal08.wmf (7)

где Nj – заселённость j-го колебательного уровня осциллятора при заданной температуре.

Nj = N0·exp(Ej/kT). (8)

При температуре Т = 315 °К имеем: N1 = N0·0,633, N2 = N0·0,4016,... N6 = N0·0,067. Ограничение рассматриваемых энергетических состояний j ≤ 6 даёт погрешность при определении хср не более 15 % , хср = 0,00196 нм.

Асимметричность веерного колебания полярной «головки» фосфолипидов (хср ≠ 0) определяет величину изменения радиальной составляющей дипольного момента полярной «головки» (Δp = Δr∙q).

2. Рассмотрим ряд параметров системы взаимодействующих молекул фосфолипидов, образующих сферическую мембрану. Одним из таких параметров является напряжённость электростатического поля внутри бимолекулярного слоя и в пространстве, прилегающем к этому слою (рис. 2).

gal2.tif

Рис. 2. Напряжённость электростатического поля бимолекулярного слоя модельной мембраны

Напряжённость электростатического поля определим в пространстве, удалённом от полярных «головок» молекул не более чем на 3 нм. При нахождении величины напряжённости электростатического поля в каждой точке обозначенного пространства учитываем влияние ближайших (к данной точке) 1000 ионов. Благодаря экранирующему действию цитоплазмы сомы нейрона (ε = 80), погрешность, возникающая при таком ограничении, не превышает 10 %.

gal09.wmf (9)

В пределах пространства каждого иона, диаметр которых принят равным 0,1 нм, напряжённость поля не рассматривается. На рис. 2 представлена радиальная зависимость величины напряжённости электростатического поля; расстояние между соседними отрицательными ионами 0,4 нм, α0 = 300, плечо диполя d = 0,6 нм (рис. 1), заряды ионов q = ± 1,6·10-19 Кл.

При изменении температуры мембраны изменяется заселённость энергетических уровней осцилляторов веерного типа полярных «головок» фосфолипидов. Вследствие асимметричности этого колебания меняется радиальная составляющая дипольного момента (Δp = Δr∙q) полярных «головок» молекул и изменяется напряжённость электростатического поля в окружающем пространстве.

Возьмём температуру термостата, в котором находится нейрон (независимо от природы термостата) равной Т = 315 °К, температуру фазового перехода «жидкий кристалл – гель» равной 310 °К. При нагревании мембраны нейрона от 310 °К до 315 °К появляется дополнительное электростатическое поле ΔЕ. Средняя величина этого поля, определяемая изменением радиальной составляющей диполя полярной «головки» (Δр = xcp·cosα0 ·q), в пространстве бислоя мембраны составляет ≈1 % от величины электростатического поля, представленного на рис. 2, т.е. в пределах бислоя мембраны ΔЕ = 104 –107 В/м. Для жидкого кристалла (мембраны), находящегося при температуре 315 °К близкой к температуре фазового перехода (310 °К), напряжённости этого поля достаточно, чтобы перевести кристалл из высокотемпературной фазы в низкотемпературную. При этом фазовом переходе молекулы фосфолипидов переходят из одного конформационного состояния (Ф1,j) в другое (Ф2,j), изменяя свою структуру, а также геометрию электронного облака.

3. Взаимодействие диполей полярных «головок» липидов в мембране приводит к образованию из каждого энергетического уровня веерных осцилляторов энергетической зоны. Оператор диполь-дипольного взаимодействия имеет вид

gal10.wmf, (10)

где gal11.wmf – операторы дипольных моментов ангармонических осцилляторов n и m (m = n ± 1), gal12.wmf – радиус вектор, θ и φ – сферические углы, задающие положение осцилляторов в мембране.

Гамильтониан системы N осцилляторов мембраны с учётом диполь-дипольного взаимодейтвия ближайших соседей можно записать:

gal13.wmf (11)

где Hi,j – гамильтониан изолированной молекулы, находящейся в состоянии i (i = 1,2).

Учитывая трансляционную симметрию мембранной структуры, волновую функцию, соответствующую оператору Гамильтона (11), можно представить (в нулевом приближении) в виде

Ψijq(θ, φ) =N-1/2Σn Ψij(θ + nθo, φ + nφо) exp(iqn),(12)

где Ψij – волновые функции изолированных молекул, q = (2πn1/N1θo + 2πn2/N2φо) – волновой вектор, N1N2 = N; θo = 2π/N1 , φо = 2π/N2 (θo = φо = 2π α /2πRo ≈ 2∙10-5рад, R0 – радиус мембраны, α ≈ 0,4 нм – расстояние между соседними молекулами).

При малом перекрывании волновых функций соседних молекул получим выражение для собственных значений гамильтониана:

Еijq = Еij + 2•J•cos(q(θ, φ)), (13)

J = ∫Ψ*ij(θ, φ)∙Wn,n±1∙Ψij(θ, φ) dθ dφ

где Еi,j – собственные значения гамильтониана изолированных молекул фосфолипидов.

gal3.tif

Рис. 3. Зонная структура фононных состояний i = 1, j = 0,1,2…6 модельной мембраны, образованной молекулами фосфатидилхолина. Энергия колебательных состояний отсчитывается от основного электронного состояния, величина которого принята равной (–5 эв = 5000 см-1)

На рис. 3 представлена схема энергетических зон, образованных из колебательных уровней j = 0, 1, 2, … 6 полярных «головок» молекул фосфолипидов, находящихся в электронном состоянии i = 1. Энергия колебательных уровней отсчитывается от основного электронного состояния молекулы фосфолипида, величина которого принята равной (–5 эв = –5000 см-1) – типичная для органических молекул типа фосфолипидов [12].

Подобную схему энергетических зон Еi,j,q можно представить и для мембраны из молекул фосфолипидов, находящихся в конформационном состоянии i = 2 (низкотемпературная фаза).

Заключение

Фазовый переход мембраны из высокотемпературной жидкокристаллической фазы в низкотемпературную гель фазу под действием электростатического поля веерных колебаний полярных «головок» фосфолипидов должен сопровождаться квантовыми переходами между зонами, так как энергетические зоны этих двух фаз при одинаковой температуре должны иметь разную заселённость.


Библиографическая ссылка

Галанов Е.К. СОБСТВЕННЫЙ ПОТЕНЦИАЛ ДЕЙСТВИЯ НЕЙРОНА С МОДЕЛЬНОЙ МЕМБРАНОЙ // Международный журнал прикладных и фундаментальных исследований. – 2018. – № 5-2. – С. 312-317;
URL: https://applied-research.ru/ru/article/view?id=12261 (дата обращения: 19.04.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674