Осцилляторная неустойчивость и спонтанные подпороговые колебания в сети диффузионно связанных кальциевых осцилляторов Текст научной статьи по специальности «Физика»

Изв. вузов «ПНД», т. 17, № 2, 2009 УДК 577.381, 530.182

ОСЦИЛЛЯТОРНАЯ НЕУСТОЙЧИВОСТЬ И СПОНТАННЫЕ ПОДПОРОГОВЫЕ КОЛЕБАНИЯ В СЕТИ ДИФФУЗИОННО СВЯЗАННЫХ КАЛЬЦИЕВЫХ ОСЦИЛЛЯТОРОВ

В.Б. Казанцев, А.В. Воробьёв

Работа посвящена исследованию динамики сети взаимодействующих астроцитов -глиальных клеток мозга, способных генерировать химические сигналы активности -кальциевые импульсы. Астроциты, подобно нервным клеткам (нейронам), формируют сети, где взаимодействие между клетками осуществляется через так называемые щелевые контакты (в англ. gap junctions). Такие контакты представляют собой специальные каналы, по которым химически активные вещества диффундируют между клетками. Как считается, кальциевые сигналы в астроцитах могут регулировать эффективность синаптической передачи в близлежащих нейронных клетках. В работе исследуются процессы возникновения колебательной активности в одномерной сети диффузионно связанных астроцитов. Динамика отдельной клетки описывается системой нелинейных дифференциальных уравнений 3-го порядка, полученной на основе детального описания кинетики биохимических преобразований в клетке (Де Янг и Кайзер, 1992; Ли и Ринцель, 1994; Уллах и др., 2006). Рассматриваемая модель сети при учете диффузионной связи представляет собой трехкомпонентную дискретную систему типа «реакция-диффузия» с одной диффундирующей переменной. Установлено, что существует критическое значение коэффициента диффузии, превышение которого приводит к появлению осцилляторной неустойчивости на частотах 0.1 Гц и генерации спонтанных квазисинусоидальных колебаний малой (0.05 мкмоль) амплитуды. Получен соответствующий спектр собственных значений, анализ которого показывает, что при дальнейшем увеличении коэффициента связи в сети возникает мультичастотный режим, что приводит, в частности, к модуляции амплитуды квазисинусоидальных колебаний и спонтанной генерации кальциевых импульсов.

Ключевые слова: Кальциевые сигналы, подпороговые колебания, диффузионная неустойчивость, астроцит.

Введение

Исследование динамических процессов возникновения структур или паттернов активности в нелинейных распределенных системах становится в последние годы одним из наиболее интенсивно развивающихся направлений нелинейной физики колебаний и волн. Ярким приложением этих исследований являются процессы генерации и преобразования сигналов в сетях мозга, лежащие в основе его функциональной активности и когнитивных функций [1-5]. Нервные клетки мозга способны

123

генерировать электрические импульсы возбуждения (потенциалы действия), которые передаются по сети за счет межнейронных связей. Эти процессы достаточно полно описаны с помощью динамических моделей различной степени сложности, от простейших триггерных сетей до сложных уравнений типа Ходжкина-Хаксли, детально учитывающих динамику ионных токов через клеточную мембрану нейрона [5]. Однако кроме нейронов в мозге существует большое количество глиальных клеток (астроциты, микроглия, миелиновые клетки и др.), которые обеспечивают поддержание жизнедеятельности нейронных систем. Интенсивные исследования физиологов последних десяти лет свидетельствуют о том, что кроме поддерживающей роли глиальные клетки, в частности астроциты, способны влиять на нейронную активность за счет синтеза нейроактивных веществ (нейропередатчиков, нейромедиаторов), выбрасываемых в межклеточное пространство и захватываемых рецепторами нейрональных клеток [6]. Кроме того, было показано, что астроциты также формируют сети и их химическая активность может существенно зависеть от коллективных эффектов за счет транспорта активных веществ между клетками. На уровне отдельных клеток астроциты, в отличие от нейронов, не обладают свойством электрической возбудимости - пороговой генерации электрических импульсов, поскольку не имеют потенциал-зависимых ионных каналов. Однако они способны генерировать импульсы внутриклеточной концентрации кальция (кальциевые импульсы или кальциевые колебания). Эти процессы могут быть описаны с помощью детального анализа кинетики биохимических преобразований внутри клетки [7]. Модель в этом случае представляет собой многомерную систему нелинейных дифференциальных уравнений, построенных аналогично формализму Ходжкина-Хаксли. Эти уравнения в некоторых случаях (в пренебрежении быстрыми процессами релаксации и при аппроксимации некоторых переменных их равновесными значениями) могут быть сведены к нелинейным системам второго [8] или третьего порядка [9]. Анализ этих моделей показывает, что при изменении контрольного параметра равновесная концентрация кальция становится неустойчивой и через бифуркацию Андронова-Хопфа происходит рождение устойчивого предельного цикла, описывающего колебания внутриклеточного кальция. В данной работе мы будем использовать трехкомпонентную модель клетки [9], переменными в которой являются: концентрация кальция в цитозоли, доля открытых кальциевых каналов эндоплазматического ретикулума (ЭР), являющегося источником кальция внутри клетки, и концентрация внутриклеточного посредника инозитол 1,4,5-трисфосфата (IP3). Взаимодействие между астроци-тами осуществляется посредством диффузии IP3 через щелевые контакты, проницаемость которых для IP3 существенно выше, чем для остальных молекул и ионов в клетке [10]. Одним из эффектов такого взаимодействия, исследованным ранее [9], является синхронизация (синфазная и антифазная) кальциевых колебаний в двух взаимодействующих клетках, подтвержденная экспериментально. Среди известных результатов по исследованию кальциевых сигналов в глиальных клетках и нейронах необходимо также упомянуть распространение кальциевых волн в диффузионных моделях IP3 чувствительных клеток [11], кальциевые колебания и волны, индуцированные глутаматом [9], нейрон-астроцитарное взаимодействие и формирование так называемого тройственного синапса [12,13] и др. [6]. Отметим, что процессы коллективной динамики взаимодействия астроцитов в больших сетях до настоящего времени практически не исследованы.

124

В данной работе рассмотрен эффект возникновения осцилляторной неустойчивости и генерации спонтанных квазисинусоидальных колебаний за счет диффузии IP3 в модели одномерной сети взаимодействующих астроцитов. Модель представляет собой дискретный аналог трехкомпонентной системы «реакция-диффузия» с одним диффундирующим компонентом (IP3). При отсутствии связи каждый из аст-роцитов находится в устойчивом состоянии равновесия и генерирует кальциевый импульс только при достаточно большом внешнем воздействии (например, в ответ на химическую стимуляцию глутаматом [9]). Мы покажем, что при введении диффузионной связи равновесная концентрация кальция становится неустойчивой при превышении критического значения коэффициента связи. Это соответствует переходу комплексной пары характеристических корней через мнимую ось и возникновению спонтанных квазисинусоидальных колебаний. Дальнейшее увеличение диффузии приводит к переходу в неустойчивую область целой части спектра и формированию многочастотного осцилляторного паттерна. Амплитуда кальциевых колебаний при этом растет и может достигнуть порога возбуждения, что приводит к спонтанной генерации кальциевых импульсов в отдельных клетках сети.

1. Модель

Рассмотрим модель одномерной сети локально связанных астроцитов, схематически представленную на рис. 1. Учет только локальных связей по внутриклеточной IP3-диффузии обоснован недавними экспериментальными исследованиями [14], которые показывают, что астроциты занимают «неперекрывающиеся территории» и, следовательно, эффективно взаимодействуют только с ближайшими соседями1.

Рис. 1. Одномерная сеть взаимодействующих кальциевых осцилляторов - астроцитов. PLC -phospholipase C; GJ - щелевой контакт; ER - эндоплазматический ретикулум. Стрелками показаны направления внутриклеточных токов активных веществ

1 Второй механизм взаимодействия астроцитов, допускающий нелокальные взаимодействия через внеклеточную диффузию АТФ [6], в данной работе не рассматривается.

125

Динамика сети из N элементов описывается следующей многомерной системой дифференциальных уравнений:

' dxj dt

dVj < dt

dzj

dt

JIP3 (xj, Vj) + к1Рз (xj-i - 2xj + xj+1 )> JCa(xj ,yj, zj),

HCa(xj, yj, zj ),

U = 1, 2,..., N.

(1)

Переменная xj = [IP3] описывает концентрацию посредника IP3, изменения которой определяются внутриклеточными преобразованиями Jip3(xj,Vj) и межклеточной диффузией с коэффициентом kip3. Переменная Vj = [Ca2+] описывает концентрацию кальция в цитозоли и определяется суммарным током кальция, Jca (xj, Vj, Zj). Переменная zj определяет долю открытых кальциево-зависимых каналов эндоплазмагического ретикулума с функцией инактивации Hca(xj ,Vj, Zj). Выражения для функций правых частей (1) имеют вид [7,9,12]

J1P3 (xj, Vj)

/ Vj + (1 - а)к4 V Vj + к4

-1 (xj - [IP3]*),

Tr

JCa(xj, Vj, zj) HCa (xj ,Vj ,zj )

JChannel JPump + JLeak + Jin Jout)

ZCO Zj

Tz '

(2)

Первое слагаемое в Jip3 описывает синтез IP3 за счет так называемого PLC-S механизма. Изменение внутриклеточной концентрации кальция JCa определяется током кальция через IPз-зависимые каналы мембраны ЭР (см. рис. 1)

•/Channel = Vixj Vj Р ■

((xj + di) (Vj + d5))3

активным током кальция, создаваемым АТФ-зависимыми насосами, закачивающими кальций обратно в ЭР,

JPump — v3

Vj2

к32 + Vj2

током утечки

/Leak = V2 (Со - (ci + 1) Vj)

и токами кальция через плазматическую мембрану астроцита

xj2

Jin = v5 + v6

Jout — k1Vj ■

к2 + x2 к2 xj

126

Функция инактивации кальциевых каналов ЭР zc деляются следующими выражениями:

и время релаксации tz опре-

z^ —

Q2

Q2 + yj 1

a2 (Q2 + Vj)!

Q2 — d2

Xj + d1 Xj + d3

Константы и параметры, входящие в выражения для функций правых частей, определяются кинетикой биохимических преобразований [7]. В соответствии с экспериментальными данными исследования культуры астроцитов [9], зафиксируем значения констант, перечисленные в таблице. Выберем в качестве контрольных параметры v4 и к\р3, определяющие динамику диффундирующей переменной. Предполагаем, что на границах сети (см. рис. 1) выполнены граничные условия «нулевого потока»:

Xo — X1, XN +1 — Xn ■

(3)

Таблица

Параметры модели

Параметр Значение

со 2.0 мкмоль

C1 0.185

vi 6 c-1

v2 0.11 c-1

v3 2.2 мкмоль-с-1

v5 0.025 мкмоль-с-1

ve 0.2 мкмоль-с-1

ki 0.5 c-1

k 1.0 мкмоль

k3 0.1 мкмоль

a3 0.14 мкмоль-1 - с-1

di 0.13 мкмоль

d2 1.049 мкмоль

d3 0.9434 мкмоль

d,5 0.082 мкмоль

a 0.8

Tr 7.143 c

[1P3]* 0.16 мкмоль

&4 1.1 мкмоль

Как было показано в предыдущих исследованиях [7,9], динамика изолированной клетки, kip3 — 0, сравнительно проста. При v4 < vAH система имеет единственное устойчивое состояние равновесия. При V4 — vAH в результате бифуркации Андронова-Хопфа мягко рождается устойчивый предельный цикл, соответствующий внутриклеточным колебаниям кальция.

т

Z

2. Возникновение осцилляторной неустойчивости в сети астроцитов

Рассмотрим теперь динамику системы при введении диффузионной связи, kiP3 > 0. Состояние равновесия каждого из элементов определяет пространственнооднородное состояние равновесия сети. Координаты этого состояния по соответствующим компонентам Xj — X*, Vj — y*, Zj — z* (j — 1, 2,..., N) определяются

следующим нелинейным уравнением:

v1

v4Tr (y* + (1 - a)&4) V* + k4

+ [IP3]* V*3a2d2 (Co - (C1 + 1)v*)

(V* + d5 )3

v4Tr (y* + (1 - a)&4) V* + k4

+ [IP3]* (y* + d2) + d3V* + did2

3

3

(4)

v3V*2 V*2 + Щ

+ v2 (C*

(Ci + 1)y*) — 0.

127

Численное решение уравнения (4) показывает, что в диапазоне изменения переменных (xj, yj, Zj > 0) и для выбранных значений констант (см. таблицу) оно имеет единственное решение.

2.1. Спектр собственных значений состояния равновесия. Для определения устойчивости пространственно-однородного состояния сети в многомерном фазовом пространстве линеаризуем систему (1) в окрестности этого состояния. Для малых возмущений вида

Xi - x* = Aieu, yi - y* = Bieu, Zi - z* = Cieu (5)

амплитуды Ai,Bi,Ci ^ 1 удовлетворяют следующей системе однородных линейных уравнений:

" Ai (JIP3x — 2kIp3 — + BiJIP3y + Ci JIP3z + Ai-1kIP3 + Ai+1kIP3 = °

AiJCax + Bi JCay - + CiJCaz = °

V ' (6)

AiHCax + BiHCay + Ci (HCaz - °

. A0 = A1, AN+1 = AN,

где J/P3xyz, JGaxyz, H'Gaxyz - частные производные функций правых частей (1) по соответствующим переменным. Исключая Bi и Ci из (6), легко получить

Ai

JIP3 x - 2kIP3 - ^ + JT1

TTf jf

H C axJC a z

JC ax(HC az

IP3

y (JCay ^)(HCaz ^) JCazHC

+

ay

(7)

+kIP3 (Ai-1 + Ai+1) = °-

Введем обозначение

5

kIP3

JIP3 x — 2kIP3 — ^ + JI

IP3 г

ff

H C ax J Caz

JC ax(HC az ^)

(JC a y ШC az JC az HC ay

1

Тогда (6) перепишется в виде

(8)

( (S + 1)A1 + A2 = 0,

< SAi + Ai-1 + Ai+1 = °, i = 2,3,..., N — 1, (9)

l (S + 1)an + an-1 = °-

Для нахождения спектра собственных значений системы (9) воспользуемся техникой, предложенной в [15] для цепочечных систем с локальной диффузионной связью. Характеристический определитель для (9) можно записать в виде

128

Qn (S)

S+1 1 0 0 • • • 0

1 S 1 0 • • • 0

0 0 0 1 S +

После преобразований Qn можно представить в виде

(10)

Qn(S) = (S + 1) ((S + 1)Tn-2(S) - Tn-s(S)) - (S + 1)Tn-s(S) + Tn-4(S), (11)

где Tn (S) - «симметричный» определитель

Tn (S)

S 1 0 •• •0

1 S 1 •• •0

0 0 0 S

для которого справедливо следующее рекуррентное уравнение:

Tn(S) = STn-i(S) - Tn-2(S).

Используя (13), можно записать

(12)

(13)

Qn (S ) = (S + 2)Tn-i(S ). (14)

Тогда корни характеристического уравнения Qn (S) = 0 будут определяться уравнениями

S = -2,Tn-i(S ) = 0. (15)

Согласно рекуррентному соотношению (13), определитель Tn (S) можно выразить через многочлены Чебышева второго рода, корни которых известны [16]. Для уравнения Tn-i(S) = 0 они имеют вид

S

2 cos

пк

N

к = 1,2,...,N - 1.

(16)

Для учёта корня S = -2 доопределим к значением к = N. Выражая S через параметры системы (1), получаем характеристическое уравнение в виде

X3 - l2(j[p3 x - kIP 3 (S + 2) + JC ay + HC az) +

+ X(JCayHCaz - JCazHCay + (JIP3x - к1Рз (S + 2))(JCay + HCaz) - JIP3y JCax)--(JIP3 y HCaxJC az-JIP3 y JC ax HC a z+(JIP3 x- к1Рз (S + 2))(JC ay HC az-JC az HC ay )) = 0

(17)

Кубическое уравнение (17) определяет 3N собственных значений пространственнооднородного состояния сети астроцитов.

129

2.2. Анализ спектра собственных значений. Исследуем расположение корней уравнения (17) на комплексной плоскости X в зависимости от контрольных параметров. Спектр состоит из N действительных отрицательных значений, не меняющих знак при изменении контрольных параметров V4,kip3, и 2N комплексносопряженных корней. При v4 < vAH и kIP3 < ЩРя (v4), где kj*P3 (v4) - некоторое критическое значение, весь спектр находится слева от мнимой оси и состояние равновесия является устойчивым (рис. 2, а). При значении коэффициента связи выше kj*p3 (V4) пара комплексно-сопряженных корней переходит через мнимую ось (рис. 2, б).

В многомерном фазовом пространстве происходит бифуркация Андронова-Хопфа и мягко рождается устойчивый предельный цикл2. Ему соответствует возникновение в сети малых квазисинусоидальных колебаний внутриклеточного кальция с частотой fi ~ Im Xi/2n. Дальнейшее увеличение параметра kip3 соответствует последовательному переходу комплексно-сопряженных пар в неустойчивую область (рис. 2, в, г). В сети формируется многочастностный режим с fi ~ Im Xi/2n, i = 1, 2,...,N. Зависимость действительной части собственных значений спектра от коэффициента связи представлена на рис. 3, а. Таким образом, в системе происходит каскад локальных бифуркаций, определяемых последовательностью бифуркационных параметров k\+, k22+, ...,k*j+, где M = M(v4) < N - максимальное число неустойчивых пар в спектре и, соответственно, частот, зависящее от параметра V4. Отметим, что речь идет о частотах собственных неустойчивых мод в линейном приближении, которое справедливо только при достаточно малых амплитудах колебаний сети. При сравнительно больших значениях kip3 в спектре происходит каскад обратных переходов k^~, k2-, ... ,k^-, при котором собственные значения

Рис. 2. Расположение комплексно-сопряженных корней спектра (17) пространственно-однородного состояния равновесия сети для N = 6 на комплексной плоскости. Параметры: а - к1Рз = 0.012 мкмоль, б - к1Рз = 0.023 мкмоль, в - к1Рз = 0.032 мкмоль, г - максимальный «порядок» неустойчивости M = 5, к1Рз = 0.8 мкмоль. Остальные параметры приведены в таблице

2 Сценарий мягкого рождения автоколебаний (первая ляпуновская величина отрицательна) будет проиллюстрирован ниже плавной монотонной зависимостью амплитуды от параметра.

130

Рис. 3. Зависимость действительных частей собственных значений спектра (17) от коэффициента связи к1Рз. Стрелками показаны соответствующие бифуркационные значения параметров. а - Переход собственных значений через ноль при сравнительно малых значениях к1Рз. б - Возврат собственных значений в устойчивую область при больших значениях коэффициента связи

возвращаются в устойчивую область и, соответственно, количество неустойчивых «мод» (колебаний на собственных частотах в линейном приближении) уменьшается, а при k > kM пространственно-однородное состояние сети вновь становится устойчивым (рис. 3, б).

3. Спонтанная колебательная активность сети астроцитов

3.1. Малые квазисинусоидальные колебания. Исследуем колебания сети астроцитов, возникновение которых при увеличении коэффициента связи предсказывает теоретический анализ, с помощью численного моделирования уравнений (1). При параметрах связи ниже первого критического значения к^+ весь спектр находится слева от мнимой оси и все возмущения затухают к устойчивому пространственнооднородному состоянию сети (рис. 4, а). При к\+ < kip3 < к22+ в системе существует единственный устойчивый предельный цикл, что соответствует квазисинусоидальным колебаниям элементов сети, амплитуда которых и частота (близкая к /1 ~ 0.1 Гц) не меняется во времени (рис. 4, б). Отметим, что амплитуда колебаний монотонно спадает при приближении к границам сети. В случае существования в спектре нескольких неустойчивых значений при kip3 > к\+ характер колебаний сети становится более сложным (рис. 4, в). Амплитуда колебаний становится медленно модулированной функцией времени. Этот эффект можно объяснить следующим образом. В линейном приближении несколько собственных значений в спектре определяют развитие неустойчивости по нескольким собственным направлениям (многообразиям - в многомерном фазовом пространстве) - «собственным модам». При достаточно малой амплитуде колебаний формирование колебательного паттерна можно трактовать как результат слабого нелинейного взаимодействия собственных мод, что ведет, в частности, к появлению комбинационных частот в спектре. Поскольку частоты неустойчивых мод достаточно близки между собой, наличие медленной огибающей соответствует сравнительно малой разностной частоте, А/ ^ /.

131

3.2. Генерация кальциевых импульсов. Появление амплитудной модуляции при развитии неустойчивости и взаимодействии нескольких неустойчивых мод является принципиально важным для появления кальциевых импульсов - выбросов сравнительно большой амплитуды. Эти импульсы возникают спонтанно в различных элементах сети (рис. 4, г), стимулируя передачу возбуждения к соседним элементам за счет диффузионной связи. Эффект связан, с одной стороны, с возбудимыми свойствами локального элемента сети - генерации кальциевого отклика при наложении достаточно сильного внешнего возмущения. С другой стороны, модуляция амплитуды приводит к повышению размаха колебаний в определенных элементах сети и достижению порога генерации кальциевых импульсов в определенные моменты времени. Таким образом, взаимодействие элементов сети за счет диффузионной связи и появление спонтанных подпороговых колебаний, фактически, обеспечивает снижение эффективного порога возбуждения и, следовательно, возможность межэлементной «передачи» возбуждения. Отметим, что особенностью генерации и передачи импульсов при наличии подпороговых колебаний является также определенный временной масштаб генерации импульсов (в англ. timing), поскольку возбуждение происходит только на максимумах квазисинусоидальных подпороговых возмущений. Этот временной масштаб составляет для сети астроцитов примерно 10 секунд. Обратный переход через последовательность бифуркационных значений ксоот-ветствует возвращению собственных значений в устойчивую область, уменьшению числа неустойчивых мод и приводит к снижению характерной амплитуды спонтанных колебаний. Как следствие, сеть возвращается в режим подпороговых колебаний с подавлением импульсной генерации.

Рис. 4. Характерный вид кальциевых колебаний в цепочке при различных значениях коэффициента диффузии к1Рз (мкмоль): 0.012 - затухающие колебания в окрестности устойчивого фокуса (а); 0.023 -квазисинусоидальные подпороговые колебания (б); 0.032 - амплитудная модуляция в многочастотном режиме (в); 0.8 - генерация кальциевых импульсов на пиках подпороговых колебаний (г). Остальные параметры приведены в таблице

132

3.3. Характеристики колебаний. Исследуем характеристики кальциевых сигналов, возникающих в сети астроцитов вследствие осцилляторной неустойчивости. Рассмотрим сначала представление колебаний в многомерном фазовом пространстве системы (1). При fcjp3 < к\+ состояние равновесия является устойчивым фокусом (рис. 5, а). При превышении критического значения, согласно проведенному теоретическому анализу, в фазовом пространстве системы (1) рождается устойчивый предельной цикл малой амплитуды, проекция которого на плоскость переменных (xi, yi) представлена на рис. 5, б. Предельный цикл располагается в окрестности интегрального многообразия состояния равновесия (плоскость в линейном приближении), которое при к\+ < kip3 < к^+ является седло-фокусом с размерностью устойчивого многообразия dim Wf_ = 3N — 2 и неустойчивого - dim WU = 2. Отметим, что в проекции на плоскости переменных других элементов сети амплитуда цикла будет спадать при приближении к границам. При переходе через следующее бифуркационное значение к\+ вторая пара корней переходит через мнимую ось и размерность неустойчивого многообразия возрастает, dim W2U = 4. Поскольку предельный цикл расположен вблизи многообразия W“, то этот переход будет соответствовать потере устойчивости предельного цикла и рождению в многомерном пространстве тора. Движения на этом торе отвечают режиму биений, проекция которых на плоскость переменных (xi,yi) представлена на рис. 5, в. Режим генерации кальциевых импульсов, возникающий при дальнейшем увеличении kip3 характеризует возбудимые свойства сети астроцитов (рис. 5, г). Кальциевым импульсам соответ-

Рис. 5. Проекции фазового пространства системы (1) на плоскость (x1,y1), отвечающие различным значениям коэффициента диффузионной связи к1Рз между элементами (мкмоль): 0.012 - устойчивый фокус (а); 0.023 - устойчивый предельный цикл (б); 0.032 - режим биений (в); 0.8 - режим импульсной генерации (г). Остальные параметры приведены в таблице

133

ствуют характерные выбросы траекторий в проекции на плоскость соответствующего элемента в область более высоких концентраций кальция. Динамика фазовых траекторий в этом случае напоминает спонтанную генерацию импульсов в возбудимых системах типа ФитцХью-Нагумо при воздействии внешних периодических или шумовых возмущений [17].

Наличие возбудимых свойств в генерации кальциевых сигналов в сети аст-роцитов характеризует рис. 6, а, иллюстрирующий зависимость средней амплитуды колебаний сети от коэффициента диффузии kip3. При увеличении kip3 выше первого бифуркационного значения происходит плавное увеличение средней амплитуды колебаний, что свидетельствует о мягком рождении предельного цикла в многомерном пространстве. При достижении порога генерации кальциевых импульсов амплитуда скачком возрастает и далее становится ограниченной, что является еще одной характеристикой сети астроцитов как возбудимой системы. Колебания на профиле амплитуды выше порога соответствуют различному числу импульсов, возникающих на фиксированном интервале времени. Другим словами, импульсы появляются не на каждом пике подпороговых колебаний (см. рис. 4, г), и характер их генерации является нерегулярным. Это свойство иллюстрирует бифуркационная диаграмма рис. 6, б, на котором представлена зависимость наблюдаемых величин межимпульсного интервала (в англ. interspike interval, ISI) для всех элементов сети от коэффициента диффузии kip3. Наряду со сравнительно малыми интервалами (один-два периода подпороговых колебаний), являющимися наиболее вероятными (наибольшая плотность точек на диаграмме), в элементах возможны сравнительно большие интервалы эволюции в подпороговом режиме без генерации.

к1Рз между элементами. Резкое увеличение амплитуды при достижении коэффициентом к1Рз значения 0.0346 мкмоль соответствует появлению кальциевых импульсов. б - Бифуркационная диаграмма наблюдаемых в N = 6 элементах сети интервалов следования кальциевых импульсов в зависимости от коэффициента связи

134

Заключение

В работе проведено исследование коллективной динамики сети диффузионносвязанных кальциевых генераторов, моделирующих активность глиальных клеток мозга - астроцитов. Отметим наиболее важные его результаты.

• Локальная диффузионная связь в сети взаимодействующих кальциевых осцилляторов (астроцитов) приводит к неустойчивости однородного распределения концентрации внутриклеточного кальция. Потеря устойчивости происходит через бифуркацию Андронова-Хопфа.

• Неустойчивость имеет осцилляторный характер и выражается в возникновении спонтанных квазисинусоидальных коллективных колебаний элементов.

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

• На пиках подпороговых колебаний возможна спонтанная генерация кальциевых импульсов сравнительно большой по отношению к подпороговым колебаниям амплитуды.

Возникновение подпороговых колебаний вследствие неустойчивости характеризует два важных свойства динамики сети астроцитов. С точки зрения элемента сети подпороговые колебания, возникающие за счет диффузии IP3, фактически снижают эффективный порог генерации импульсов (см. рис. 4, г). Более того, такая генерация вследствие коллективного взаимодействия обладает характерным временным масштабом (в англ. timing), определяемым периодом подпороговых колебаний. Наличие такого сравнительно медленного масштаба у глиальных сетей может оказывать существенное влияние на временные характеристики регуляции активности основных сигнальных клеток мозга - нейронов. Заметим также, что неустойчивость и подпороговые колебания в сети астроцитов (1) были получены в биологически релевантном диапазоне изменения контрольных параметров, выбранных согласно проведенным ранее экспериментальным исследованиям.

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

Работа поддержана грантами РФФИ (08-02-00724, 06-02-16137), грантом Президента Российской Федерации (МД-4602.2007.2), Фондом содействия отечественной науке, Программой Президиума РАН «Молекулярная и клеточная биология».

Библиографический список

1. Николлс Дж., Мартин Р, Валлас Б. Фукс П. От нейрона к мозгу. М.: Изд-во УРСС. 2003. 672 C.

2. Рубин А.Б. Биофизика: В 2 т. М.: Книжный дом «Университет» , 2000. 486 с.

135

3. Principles of Neural Science / Eds. E.R. Kandel, J.H. Schwartz, T.M. Jessell. Third Edition. Prentice-Hall Intern. Inc. 1991. 1135 p.

4. Scott A. Neuroscience: a mathematical premier. Springer-Verlag, Berlin, 2002.

5. Izhikevich E.M. Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting, MIT Press, 2007.

6. Verkhratsky A., Butt A. Glial Neurobiology, Wiley, 2007.

7. De Young G.W., Keizer J. A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration //Proc. Natl. Acad. Sci. USA. 1992. Vol. 89. P. 9895.

8. Li Y., Rinzel J. Equations for InsP3 receptor-mediated [Ca2+]i oscillations derived from a detailed kinetic model: a Hodgkin-Huxley-like formalism //J. Theor. Biol. 1994. Vol. 166. P. 461.

9. Ullah G., Jung P., Cornell-Bell A.H. Anti-phase calcium oscillations in astrocytes via inositol (1,4, 5)-trisphosphate regeneration //Cell Calcium. 2006. Vol. 39. P. 197.

10. Bennett M.V.L., Jorge E.C., Bukauskas F.F., Saez J.C. New roles for astrocytes: Gap junction hemichannels have something to communicate //Trends in Neurosciences. 2003. Vol. 26, № 11. P. 610.

11. Timofeeva Y., Coombes S.. Wave bifurcation and propagation failure in a model of Ca2+ release // J. Math. Biol. 2003. Vol. 47, № 3. P. 249.

12. Nadkarni S., Jung P Spontaneous oscillations of dressed neurons: a new mechanism for epilepsy? // Phys. Rev. Lett. 2003. Vol. 91. P. 268101.

13. Volman V., Ben-Jacob E., Levine H. The astrocyte as a gatekeeper of synaptic information transfer // Neural Comp. 2007. Vol. 19. P. 303.

14. Halassa M.M., Fellin T., Takano H., Dong J.-H., Haydon P.G. Synaptic islands defined by the territory of a single astrocyte // J. Neurosci. Vol. 27, № 24. P. 6473.

15. Nekorkin V.I., Velarde M.G. Synergetic Phenomena in Active Lattices. Berlin: Springer, 2002.

16. CRC Standard Mathematical Tables and Formulae / Ed. D. Zwillinger. Boca Raton, FL: CRC Press, 1995.

17. Makarov V.A., Nekorkin V.I., Velarde M.G. Spiking behavior in a noise-driven system combining oscillatory and excitatory properties // Phys. Rev. Lett. 2001. Vol. 86, № 15. P. 3431.

Поступила в редакцию 22.10.2008 После доработки 1.02.2009

OSCILLATORY INSTABILITY AND SPONTANEOUS SUBTHRESHOLD OSCILLATIONS IN A NETWORK OF DIFFUSIVELY COUPLED CALCIUM OSCILLATORS

V.B. Kazantsev, A.V. Vorobiev

The paper is devoted to the investigation of the dynamics of a network of interacting astrocytes. The astrocytes represent brain glial cells capable to generate chemical activity signals (calcium pulses). Similarly to nerve cells (neurons) the astrocytes form networks

136

of interacting units coupled by means of gap junctions. The junctions represent special protein channels providing the diffusion of chemically active species between neighboring cells. It is believed that calcium signals in astrocytes can regulate the efficiency of synaptic transmission in neighboring neuronal cells. In the present paper we investigate the processes of oscillatory activity generation in a one-dimensional network of coupled astrocytes. The dynamics of local cell is described by the third order nonlinear differential equation system that has been obtained from the detailed description of biochemical kinetics in the cell (de Young and Keizer, 1992; Li and Rinzel, 2003; Ullah, et al., 2006). The model accounting for the diffusive coupling represent a three-component reaction-diffusion network with single diffusing component. It is proven that there exists a critical value of diffusion coefficient above which the oscillatory instability at 0.1 Hz frequency develops and spontaneous low-amplitude quasisinusoidal oscillations (of 0.05 pM) appear. Corresponding eigenvalue spectrum is obtained and analyzed. It is found that further increase of the coupling coefficients leads to the appearance of multi-frequency mode with the modulation of the oscillation amplitude and spontaneous calcium pulse generation.

Keywords: Calcium signals, subthreshold oscillations, diffusive instability, astrocyte.

Казанцев Виктор Борисович - родился в Дзержинске, Горьковской области (1973), окончил радиофизический факультет ННГУ (1996). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1999), доктора физико-математических наук (2006). В настоящее время является заведующим лабораторией «Нелинейных процессов в живых системах» Института прикладной физики РАН и заведующим кафедрой нейродинамики и нейробиологии на биологическом факультете ННГУ им. Н.И. Лобачевского в Нижнем Новгороде. Область научных интересов - теория колебаний и волн, структуры и волны в многоэлементных сетевых системах, нейродинамика. Имеет более 70 научных публикаций.

E-mail: vkazan@neuron.appl.sci-nnov.ru Институт прикладной физики РАН 603950, Нижний Новгород, ул. Ульянова, 46

Воробьев Артем Викторович - родился в Горьком (1987). Студент магистратуры радиофизического факультета Нижегородского государственного университета им. Лобачевского по направлению «Радиофизические методы в нейробиологии». Область научных интересов - динамика нейроно-глиальных систем, вычислительная математика. Имеет 4 научные публикации.

E-mail: art.vorobiev@gmail.com

Нижегородский государственный университет им. Н.И. Лобачевского 603950, Нижний Новгород, пр. Гагарина, 23

137