Регулярные и хаотические колебания в модели астроцита с регуляцией динамики кальциевого выброса Текст научной статьи по специальности «Математика»

УДК 621.391.01

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

В. В. Матросов1, Н. В. Болдырева1, В. Б. Казанцев1,2

1 Нижегородский государственный университет им. Н.И. Лобачевского 2 Институт прикладной физики РАН

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

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

Введение

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

• При срабатывании синаптической связи между нейронами происходит выброс во внеклеточное пространство нейроактивных химических веществ, нейропе-редатчиков. Эти вещества детектируются соседним нейроном, что приводит к

возникновению ответного сигнала. Однако часть нейропередатчика диффундирует от синаптического контакта и может быть захвачена соседними астро-цитами. Это приводит к каскаду молекулярных преобразований, результатом которых является генерации кальциевых импульсов. С точки зрения нелинейной динамики, этот процесс подобен классической схеме автогенерации колебаний [3].

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

• Астроциты связаны между собой посредством так называемых щелевых контактов, динамика которых описывается линейной диффузионной связью, образуя таким образом в мозгу сети астроцитарных клеток. Через контакты, в частности, могут проходить молекулы инозитол-трифосфата ( ИТФ). Как было показано в недавних модельных работах, эффект диффузионной связи может привести к спонтанной генерации кальциевых колебаний в сети астроцитов, возникновению мультистабильности и сложным (хаотическим) режимам [7-9]. В распределенных моделях кальциевой сигнализации внутриклеточная диффузия ИТФ приводит к возникновению сложных колебательных режимов, включая режимы самоорганизованной критичности [10].

• Растущий интерес исследователей к астроцитам связан также и с тем, что эти клетки контактируют как с нейронами, так и с клетками кровеносных сосудов, что делает их удобной «мишенью» для доставки веществ, воздействующих на мозг, и лекарственных средств [2].

Экспериментальные исследования астроцитов показывают, что динамика ИТФ во многом является определяющей в процессе возникновения кальциевых колебаний, поскольку ИТФ фактически «запускает» кальциевый генератор, обеспечивая выход ионов кальция из внутриклеточных хранилищ. ИТФ вовлечен в множество сложных молекулярных каскадов, часть из которых не учитывается в традиционных упрощенных моделей астроцитов. Учет динамики ИТФ приводит к трехкомпонент-ным моделям, где динамические режимы могут включать возникновение сложных хаотических колебаний (см., например, [6]). В одной из недавних моделей астроци-тарной динамики, предложенной де Питта с соавторами [5], были учтены процессы нелинейной релаксации ИТФ за счет активности протеинкиназы С (ПКС). Такая модификация позволила получить кальциевые сигналы более сложной и реалистичной формы. Однако, динамические механизмы возникновения таких колебаний описаны не были.

В данной работе мы исследуем модель де Питта [5] кальциевого генератора с регуляцией кальциевого выброса с помощью бифуркационного анализа и показываем, что учетом такой регуляции приводит к усложнению бифуркационных механизмов генерации колебаний, возникновению мультистабильности и хаотических колебаний даже в одиночной клетке. Мы также исследуем основные динамические характеристики колебаний, изменения их формы с помощью анализа фазовых проекций и построения одномерных бифуркационных диаграмм.

1. Модель

Рассмотрим модель химической сигнализации астроцита с учетом нелинейной релаксации ИТФ, предложенную в работе [5]:

61 = ° + 0+1н(С2,Кь) - °шн(°4,Кв)н{1,Кш) -

dC ~dt

— = pa(QNmlh3 + Ql)(Cq - (1 - PA)C) - OpH(C2,KP), (1)

бК = - К),

где I - концентрация ИТФ, С - внутриклеточная концентрация кальция, К - доля инактивированных ИТФ-зависимых каналов эндоплазматического ретикулума (ЭР), (0<К<1), Н(х,К,п) = хп/(хп+Кп) - функция Хилла (х - переменная, п и К -параметры), т= Н(I, 6{)Н(С, 4), = °2(Я2 + С), Q2 = 6^(1 + )/(1 + 6з), = Я2/(Я2 + С). Константы и параметры, входящие в выражения для правых частей модели (1), определяются кинетикой биохимических преобразований, и установлены в соответствии с экспериментальными данными исследования культуры астроцитов. Значения констант и диапазоны изменения параметров для модели (1) приведены в Таблице согласно данным статьи [5].

Таблица

Параметр Описание Значения Размерность

Ор Скорость продуцирования ИТФ

при активации фосфолипазы С в, PLCP 0

Оs Скорость продуцирования ИТФ при

активации фосфолипазы С S, PLCS >1

ks Константа подавления фосфолипазы С S 1.0 цМ

Озк Максимальная скорость деградации ИТФ

за счет 3-киназы >1

Q6p Максимальная скорость деградации ИТФ

за счет дефосфолирирования 0.01 s-1

Аффинность Ca2+ к PLCS 0.5 цМ

Kd Аффинность Ca2+ к 3-киназе 0.5 цМ

Кзк Аффинность ИТФ к 3-киназе 1.0 цМ

Pa Соотношение между объемами ЭР и цитозоли 0.5 -

Qn Скорость выброса Ca2+ через ИТФ зависимые

рецепторы ЭР 6.0 s-1

Ql Скорость утечки Ca2+ через мембрану ЭР 0.05 s-1

Co Полная концентрация Ca2+ в клетке 7.2 цМ

Ор Максимальная скорость закачки кальция в ЭР

активными насосами 1 цМ-1s-1

Kp Аффинность Ca2+ к активным насосам 0.1 цМ

di Константа диссоциации ИТФ 0.1 цМ

d2 Инактивационная константа диссоциации Ca2+ 1.05 цМ

ds Константа диссоциации ИТФ 0.1 цМ

d6 Активационная константа диссоциации Ca2+ 0.08 цМ

O2 Скорость связывания инактивационных гейтов

ИТФ-зависимых рецепторов с Ca2+ 0.2 s-1

Щх,К,п)

Отметим, что процессами обмена кальция с внешней средой в модели (1) пренебрегается. Виды функции Хилла для модельных значений параметров представлены рис. 1. Основное отличие модели (1) от предшественников (например, также трехкомпонент-ной модели Уллаха [11]) заключается в учете молекулярных механизмов нелинейной деградации ИТФ, опосредованной активностью различных ферментов.

0.8

0.6

0.4

0.2

0

0

0.5 1.0 1.5

х

Рис. 1. Вид функции Хилла при значениях парамет- Некоторые результаты исследования моров: п = 1, К = 1 (линия 1), п = 2, К = 0.5 (ли- дели (1) приведены в [5]. Эти результа-

ыля 7); п 4, К °.5 (линия 3), п 2, К ты подтверждают адекватность предло-

сывают базовые динамические режимы в виде однопараметрических зависимостей. В данной работе мы проводим расширенный бифуркационный анализ динамики модели (1) в плоскости двух параметров. Он основан на методах теории колебаний и теории бифуркаций с использованием программного комплекса ДНС [12]. Далее остановимся на анализе полученных результатов.

1.1. Анализ равновесных состояний. В фазовом пространстве системы (1) в зависимости от значений параметров может существовать либо одно - 0\ (Оз), либо три - О1, О2, Оз состояния равновесия (рис. 2, a). Координаты устойчивых состояний равновесия определяют стационарные концентрации ИТФ, кальция и долю открытых каналов ЭР. На рис. 2, б представлено разбиение плоскости параметров (Озк, Ой) модели (1) на области с различным числом и устойчивостью состояний равновесия. Здесь область существования трех состояний равновесия ограничена штрихпунктирными линиями 1 и 2. При выходе из этой области через линию 1 сливаются и исчезают состояния равновесия Оз и О2, на линии 2 - состояния равновесия О1 и О2. Пунктирные линии 3 соответствуют смене устойчивости состояний равновесия, которая происходит в результате бифуркации Андронова-Хопфа. Точками N1, N2, N3, N4 отмечены значения параметров, где первая ляпуновская величина Ь обращается в ноль. На участках линии 3, расположенных ниже точки N3 и выше N4, смена устойчивости состояний равновесия происходит жестко, поскольку здесь Ь>0. На остальных участках линии 3 смена устойчивости состояний равновесия происходит мягко (Ь<0). В точках N1 и N2 состояние равновесия О1 имеет два нулевых собственных числа, в этих точках линия 2 и линия 3 соприкасаются.

На рис. 2, б светло-серым цветом отмечены области Б^Б^, Б^, где модель (1) имеет одно устойчивое состояние равновесия, темно-серым цветом - область Б3 с двумя устойчивым состояниями равновесия. В области Во устойчивые состояниями равновесия в фазовом пространстве модели (1) отсутствуют. В областях БА и Бв состояния равновесия О 1 и Оз устойчивы глобально. В области Вз в фазовом пространстве системы (1) существуют два устойчивых состояния равновесия О , Оз, бассейны притяжения которых разделяет устойчивое инвариантное многообразие седла О2 (рис. 2, a). Координаты С *, I * устойчивых состояний равновесия О 1 и Оз таковы, что они определяют высокую и низкую концентрацию кальция и ИТФ, соответственно. Таким образом, равновесные значения концентраций при значени-

(линия 4)

женной математической модели и опи-

Рис. 2. К анализу состояний равновесия модели (1): a - качественный фазовый портрет при значениях параметров из области В3, б - структура плоскости параметров (Озк, Ой) при = 0.01

ях параметров из области Бз определяются начальными условиями. Отметим также, что при вариациях параметров концентрации кальция и ИТФ могут изменяться по разному. На рис. 3 приведены однопараметрические диаграммы, иллюстрирующие некоторые возможные сценарии эволюции концентраций кальция, ИТФ и доли инактивированных ИТФ-зависимых каналов ЭР при вариациях параметров Озк и О§. На диаграммах устойчивые и неустойчивые ветви отмечены сплошными и пунктирными линиями соответственно, точки К и К2 отражают смену устойчивости режимов стационарной концентрации, серые области характеризуют диапазоны изменения концентраций кальция, ИТФ и доли инактивированных каналов в автоколебательных режимах. Остановимся на анализе представленных диаграмм более подробно.

При построении диаграмм, представленных на рис. 3, б,г,е, начальные значения параметров располагались в области ВА (Озк = 0.08, Ой = 0.06). Далее, не меняя Озк, значения параметров плавно перемещались в область Бв, проходя при этом области Бз и Бс. При достижении параметром Ой значений 0.0001, осуществлялся обратный ход, из области Бв в область Б

В области БА реализуется режим высокой концентрации кальция и ИТФ, при низкой доли инактивированных каналов ЭР. С уменьшением Ой концентрации падают, доля инактивированных каналов растет, в точке К устойчивый режим стационарной концентрации замещается автоколебательным. В автоколебательном режиме концентрации кальция, ИТФ и доля инактивированных каналов ЭР периодически меняются, диапазоны изменения переменных на диаграммах отмечены серым цветом. Видно, что концентрация кальция и доля инактивированных каналов изменяются в широком диапазоне, а концентрация ИТФ практически не осциллирует. При убывании Ой амплитуда колебаний концентрации кальция уменьшается, при Ой = 0.01298 автоколебательный режим разрушается, все фазовые траектории устремляются к устойчивому состоянию равновесия Оз, определяющему режим низкой равновесной концентрации кальция и ИТФ и высокой доли инактивирован-ных каналов. Дальнейшее уменьшение Ой влечет снижение концентрации кальция и ИТФ, но эти изменения плавные и небольшие. Доля инактивированных каналов практически не меняется. При изменении значений параметра Ой в обратную сто-

О 0.01 0.02 0.03 0.04 0.05 05 0 0.01 0.02 0.03 0.04 0.05 05

д е

Рис. 3. Однопараметрические диаграммы, отражающие эволюции концентраций кальция, ИТФ и переменной Н в зависимости от Оь при О3к = 0.08 (а,в,д) и О3к = 10 (б,г,е), = 0.01

рону (при увеличении °§) координаты и состояния равновесия °з растут, что характеризует плавное увеличение концентраций кальция, ИТФ. При достижении = 0.048273 происходит резкое увеличение концентраций, которые при дальнейшем росте повышаются, но уже плавно. Доля инактивированных каналов при прохождении резко снижается. В фазовом пространстве резкое увеличение концентраций связано с исчезновением состояния равновесия °з, и переходом фазовых состояний к устойчивому состоянию равновесия °1.

Заметим, что в точке К1 теоретически смена устойчивости происходит мягко1, то есть она должна сопровождаться рождением автоколебаний малой амплитуды. В действительности же, из рис. 3, a следует, что при прохождении точки К1 концен-

1 Смена устойчивости состояния равновесия О1 происходит в результате бифуркации Андронова-Хопфа, при этом первая ляпуновская величина Ь < 0.

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

Теперь рассмотрим эволюцию концентраций кальция и ИТФ при вариациях Ой в случае Озк=10 и □бз=0.01, когда в фазовом пространстве модели (1) существует одно состояние равновесия. Эти изменения характеризуют диаграммы на рис. 3, в,г. Здесь параметр Ой сначала уменьшался от 0.6 до 0.0001, потом увеличивался от 0.0001 до 0.6. В результате наблюдались следующие явления: при уменьшении Ой до точки К1 концентрации кальция и ИТФ практически не меняются, в точке К1 наблюдается жесткий переход на автоколебательный режим, амплитуда которого при убывании Ой уменьшается до нуля (в точке К2), далее в районе Ой = 0.5 наблюдается резкое снижение концентраций. При обратном ходе наблюдаемые ранее явления сохраняются вплоть до точки К1. После точки К1 существует небольшой интервал, где модель (1) демонстрирует бистабильное поведение. При рассматриваемых значениях параметров интервал бистабильности настолько мал, что на представленной диаграмме он не виден. Из анализа представленных диаграмм и сравнения их с диаграммами на рис. 3, а,б следует, что при увеличении Озк сужаются диапазоны изменения концентраций кальция и ИТФ (как в стационарных, так и автоколебательных режимах), уменьшается зона гистерезиса, при этом интервалы скачкообразного возрастания сохраняются.

1.2. Анализ автоколебательных режимов. Область существования автоколебательных режимов приведена на рис. 4, она выделена серым цветом. Границами области служат бифуркационные кривые предельных циклов модели (1), в частности, касательной бифуркации (линия 4), бифуркации удвоения периода (линия 6), петли сепаратрис седла О2 (линия 7), многообходных петель сепаратрис седла О2 и удвоения периода циклов большой кратности (линия 9), бифуркации Андронова-Хопфа (пунктирные линии 3). Кроме линий, ограничивающих область существования автоколебаний, на рис. 4 приведены бифуркационные кривые предельных циклов (линии 5, 8), не приводящие к исчезновению автоколебаний, и состояний равновесия (линии 1-3). Из описания границ области существования автоколебательных режимов следует, что в фазовом пространстве модели (1) устойчивые предельные циклы появляются в результате бифуркации Андронова-Хопфа, касательной бифуркации, а также из го-моклинической траектории седла О2. Остановимся на анализе этих бифуркации более подробно.

В результате бифуркации Андро-нова-Хопфа теряет устойчивость состояние равновесия О1, соответствующее высокой концентрации кальция и ИТФ. На рис. 4 этой бифуркации отвечают пунктирные линии 3. При пересечении участков этой линии, расположенных выше точки N2, и между точками Nз и N4 в фазовом пространстве модели (1) в

0.001 0.01 0.1 1.0 О

Рис. 4. К анализу состояний атоколебательных режимов модели (1): двухпараметрическая бифуркационная диаграмма в плоскости параметров (Озк, Ой) при Й5Р = 0.01

окрестности состояния равновесия 0\ появляется предельный цикл Поскольку первая ляпуновская величина на рассматриваемых участках линии 3 отрицательна, то в момент рождения цикл Ь2 устойчив, а его амплитуда мала. Напомним, что смена устойчивости состояния равновесия О1 при фиксированных значениях Озк на участке линии 3, примыкающем к точке N2, происходит с ростом Об, а на участке между точками N3 и N4 при убывании Об. На участке линии 3 выше точки N4 и ниже точки N3 смена устойчивости О1 происходит за счет стягивания в точку седлово-го предельного цикла, который появляется при пересечении линии 4 слева направо, в результате бифуркации двукратного предельного цикла. Обратим внимание на то, что линия 4 проходит левее линии 3, то есть в модели (1) с ростом параметра Озк автоколебания всегда появляются жестко. В серой области, расположенной между линиями 4 и 3, устойчивый предельный цикл Ь1 существует совместно с устойчивыми состояниями равновесия, амплитуда цикла Ь1 изначально не мала, при увеличении Озк она растет. Появление в фазовом пространстве модели (1) устойчивого предельного цикла Ь1 до смены устойчивости состояния равновесия О1 объясняет наблюдаемую картину жесткого перехода на автоколебательный режим при увеличении Озк (уменьшении Об) даже в случае мягкой смены устойчивости О1.

Рассмотрим два сценария возбуждения автоколебаний в модели (1) при выходе из области В^ с уменьшением Об при фиксированных значениях Озк: когда смена устойчивости О1 происходит жестко (при Озк = 1) и мягко (при Озк = 0.1) (рис. 5).

Случай Озк = 1, рис. 5, а. В области ВА система (1) глобально асимптотически устойчива, единственным аттрактором модели является состояние равновесия О1, в области единственным аттрактором модели является цикл Ь^ Двигаясь с убыванием Об из области ВА в область ^1, последовательно пересекаем линии 4,3 (рис. 5, а): на линии 4 в результате касательной бифуркации возникают устойчивый Ь1 и седловой Г - предельные циклы, при уменьшении Об циклы Ь и Г расходятся, при этом амплитуда цикла Г убывает до нуля. Цикл Г стягивается в точку О1 на линии 3, состояние равновесия О1 теряет устойчивость жестко (Ь > 0), фазовые траектории из его окрестности устремляются к циклу Ь1 большой амплитуды. При дальнейшем уменьшении Об амплитуда цикла Ь1 убывает. При обратном ходе (при увеличении Об) амплитуда цикла Ь1 растет, при Об = 0.1784 цикл Ь1 сталкивается с неустойчивым циклом Г и исчезает. Фазовые траектории возвращаются к устойчи-

Рис. 5. Диаграмма эволюции максимумов концентрации кальция в автоколебательном режиме при изменении Об при О3к = 1 (а), Озк = 0.1 (б), = 0.01

вому состоянию равновесия 0\. Таким образом, модель (1) имеет зону бистабильно-го поведения, обусловливающая гистерезисные явления на границе существования автоколебательных режимов.

Случай Озк = 0.1, рис. 5, б. Двигаясь с убыванием Об из области ВА в область ^2, последовательно пересекаем линии 1,4,3,5: на линии 1 в фазовом пространстве модели (1) появляется устойчивое состояние равновесия Оз, которое не вносит изменений в поведение модели; на линии 4 в результате касательной бифуркации возникают устойчивый Ь\ и седловой Г\ предельные циклы, при уменьшении Об циклы Ь\ и Г расходятся; на линии 3 состояние равновесия О1 теряет устойчивость, в его окрестности возникает устойчивый предельный цикл £2 малой амплитуды, при убывании Об амплитуда цикла £2 растет, однако этот рост ограничен временем жизни цикла, которое пропорционально расстоянию между линиями 3 и 5; на линии 5 цикла £2 сталкивается с циклом Г1 и исчезает, в результате фазовые траектории из окрестности цикла £2 сравнительно малой амплитуды устремляются к циклу £1 большой амплитуды. Так как смена устойчивости состояния равновесия О1, рождение и исчезновения предельного цикла £1 происходит на малом интервале изменения параметров, то картина возникновения автоколебаний при мягкой смене устойчивости О1, аналогична картине, наблюдаемой при жесткой смене устойчивости О1.

Теперь проследим эволюцию автоколебательного режима при выходе из областей и ^2 существования автоколебательнных режимов. В процессе численных экспериментов нами выявлено большое разнообразие сценариев разрушения автоколебаний. Рассмотрим некоторые из них.

При выходе из области через участок линии 3, примыкающий к точке N2, амплитуда автоколебаний уменьшается до нуля - в фазовом пространстве модели (1) устойчивый предельный цикл £1 стягивается в точку О1, состояние равновесия О1 ставится устойчивым, и система переходит в режим стационарный концентрации. При обратном движении по параметру через кривую 3 состояние равновесия О1 теряет устойчивость, в его окрестности мягко рождается устойчивый предельный цикл £1. Таким образом, в окрестности кривой 3, примыкающей к точке N2, гистерезис отсутствует.

Рис. 6 отражает изменения концентрации кальция и межспайкового интервала в автоколебательном режиме при вариациях Озк. Из анализа представленных диа-

Рис. 6. Эволюция максимальных значений концентраций кальция (а) и межспайкового интервала (б) в автоколебательном режиме при вариациях О3к при =0.01, 0б=0.08

грамм следует, что с ростом Озк в автоколебательном режиме концентрация кальция монотонно убывает до некоторого стационарного значения, межспайковый интервал сначала возрастает, при Озк=2.0975 принимает значение ЛЬ = 30.12, далее убывает до Лt = 9.12 при Озк=9.5575.

В диапазоне 0зк € [0.16; 1.95] нижней границей существования автоколебательных режимов служит бифуркационная кривая петли сепаратрис (линия 7). В этом интервале при переходе из области в область Вс с убыванием Об сначала пересекаем линию 1 потом линию 7. При пересечении линии 1 в фазовом пространстве появляется устойчивое состояние равновесия Оз, глобальная устойчивость предельного цикла Ь\ нарушается. При пересечении линии 7 предельный цикл Ь\ исчезает, влипая в петлю сепаратрис, фазовые траектории из окрестности предельного цикла устремляются к устойчивому состоянию равновесия Оз. Поскольку цикл Ь\ исчезает в петлю сепаратрис, то период колебаний при приближении к линии 7 стремится к бесконечности, при этом амплитуда колебаний по переменным С и I остается практически неизменной. Изменение концентраций кальция и межспайко-вого интервала в автоколебательном режиме с убыванием Об иллюстрирует рис. 7.

Здесь цикл исчезает при О^ = 0.011316, где состояние равновесия Оз имеет следующие координаты Ц = 0.04491, = 0.05207, = 0.95275. Если, находясь в области Вс начать обратное движение в область увеличивая Об, то восстановление автоколебательного режима произойдет при пересечении линии 1. На рис. 7 автоколебательный режима восстанавливается при Об = 0.0384, где состояние равновесия О3 имеет координаты Ц = 0.036, С3 = 0.050, Н\ = 0.954.

В интервале Озк € [1.95, 4.345] динамика модели (1) богата и разнообразна. Здесь в зависимости от значений параметров можно наблюдать регулярные колебания с простой (рис. 8, a, б) и сложной структурой (рис. 8, в, г) (пачечную активность с различным числом колебаний в пачке), а также разнообразные хаотические колебания (рис. 8, д, е). При выходе из области существования автоколебательных режимов в рассматриваемом интервале Озк можно наблюдать хаотизацию колебаний, мягкие и жесткие переходы между хаотическими и регулярными движениями различной сложности, формирование пачечной активности и ее разрушение. Упомянутые перестройки динамических режимов происходят как при вариациях Об, так и при изменении Озк. При этом наблюдаемые картины бифуркационных переходов при увеличении О£ или Озк не совпадают с картинами, полученными при обратном

Рис. 7. Эволюция концентрации кальция (максимальных значений переменной С в автоколебательном режиме) (а) и межспайкового интервала (б), при уменьшении Об для Ош = 0.5, =0.01

Рис. 8. (1,С)-проекции и (£, С)-осциллограммы аттракторов модели (1) при Озк = 2.1, О8 = = 0.0375 (а), Озк = 4.4, О8 = 0.059 (б), Озк = 2.15, О8 = 0.037925 (в), Озк = 2.1, О8 = 0.0373 (г), Озк = 2.14, О8 = 0.037895 (д), Озк = 3.85, О8 = 0.048 (е)

Рис. 9. Эволюция характеристик (концентраций кальция С и межспайкового интервала АЬ) автоколебательного режима при уменьшении О8 в случае Озк = 2.13 и увеличении Озк в случае О8 = 0.03768, йбР = 0.01

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

Выводы

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

Проанализированы изменения концентрации кальция и ИТФ при вариациях параметров модели. Установлено, что в области малых значениях Озк существует диапазон по параметру Об, где система демонстрирует мультистабильное поведение, при котором результат установления высокой или низкой концентрации кальция и ИТФ определяется начальным состоянием (начальными условиями) или, например, внешним импульсом, перебрасывающим систему в область притяжения того или иного состояния. Динамика в области малых значениях Озк при вариациях Об характеризуется гистерезисными явлениями. В области больших значений Озк ги-стерезисные явления при вариациях Об отсутствуют, но при этом явление скачкообразного изменения концентраций сохраняется. Отметим, что возникновение биста-бильности является следствием учета в модели (1) нелинейной кинетики релаксации

ИТФ. Стационарные режимы, по сути, определяются балансом между генерацией ИТФ за счет PLC-ô и нелинейной релаксацией ИТФ.

Далее, были исследованы бифуркационные механизмы возникновения автоколебаний. Установлено, что автоколебания могут появляться в результате бифуркаций Андронова-Хопфа, двукратного предельного цикла (касательная бифуркация) и петли сепаратрис седла. Показано, что ростом параметра O3K автоколебания всегда появляются жестко. Выявлена возможность существования сложных периодических и хаотических колебаний, которые определяют пачечные (спайк-берстовые) колебания внутриклеточного кальция. Изучены бифуркационные механизмы хаотизации регулярных колебаний.

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

Работа поддержана грантом РФФИ соглашение №13-02-01223-14 от 12.03.2014.

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

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

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

3. Гордлеева С.Ю., Матросов В.В., Казанцев В.Б. Кальциевые колебания в аст-роцитах. Часть 1. Астроцит как генератор кальциевых колебаний // Известия вузов. Прикладная нелинейная динамика. 2012. Т. 20, № 3. C. 29.

4. Gordleeva S.Yu., Stasenko S.V., Semyanov A.V., Dityatev A.E., Kazantsev V.B. Bidirectional astrocytic regulation of neuronal activity within a network // Frontiers in computational neuroscience. 2012. Vol. 6. P.92.

5. De Pitta' M., Goldberg M., Volman V., Berry H., Ben-Jacob E. Glutamate regulation of calcium and IP3 oscillating and pulsating dynamics in astrocytes // Journal of Biological Physics. 2009. 35:383-411.

6. Houart G., Dupont G., Goldbeter A. Bursting, chaos and birhythmicity originating from self-modulation of the inositol 1,4,5-trisphosphate signal in a model for intracellular Ca2+ oscillations // Bulletin of Mathematical Biology. 1999. Vol. 61. P. 507.

7. Казанцев В.Б., Воробьев А.В. Осцилляторная неустойчивость и спонтанные подпороговые колебания в сети диффузионно связанных кальциевых осцилляторов // Известия вузов. Прикладная нелинейная динамика. 2009. T. 17, № 2. C. 123.

8. Матросов В.В., Гордлеева С.Ю., Казанцев В.Б. Кальциевые колебания в аст-роцитах. Часть 2. Астроцит как генератор кальциевых колебаний // Известия вузов. Прикладная нелинейная динамика. 2012. Т. 20, № 3. C. 40.

9. Matrosov V.V., Kazantsev V.B. Bifurcation mechanisms of regular and chaotic network signaling in brain astrocytes // Chaos. 2011. Vol. 21, № 2. P. 023103.

10. Wu Y.W., TangX., Arizono M., Bannai H., Shih P.Y., Dembitskaya Y., Kazantsev V.,

Tanaka M., Itohara S., Mikoshiba K., Semyanov A. Spatiotemporal calcium dynamics in single astrocytes and its modulation by neuronal activity // Philosophical Transactions of the Royal Society B: Biological Sciences. 2014. 55(2):119-129.

11. 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.

12. Матросов В.В. Динамика нелинейных систем. Программный комплекс для исследования нелинейных динамических систем с непрерывным временем. Н. Новгород: ННГУ, 2002.

13. Болдырева Н.В., Матросов В.В. Анализ динамики модифицированной модели астроцита // Труды 17-й научной конференции по радиофизике, посвященной 100-летию со дня рождения В.С. Троицкого / Под ред. C.M. Грача, А.В. Якимова. Н. Новгород: ННГУ, 2013. С. 89.

References

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

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

3. Gordleeva S.Yu., Matrosov V.V., Kazantsev V.B. // Izvestiya VUZ. Applied Nonlinear Dynamics. 2012. Vol.20, №3. Ss.29-39. (In Russian).

4. Gordleeva S.Yu., Stasenko S.K, Semyanov A.V., Dityatev A.E., Kazantsev V.B. Bidirectional astrocytic regulation of neuronal activity within a network // Frontiers in computational neuroscience. 2012. Vol.6.

5. De Pitta' M., Goldberg M., Volman V., Berry H., Ben-Jacob E. Glutamate regulation of calcium and IP3 oscillating and pulsating dynamics in astrocytes // Journal of Biological Physics. 2009. 35:383-411.

6. Houart G., Dupont G., Goldbeter A. Bursting, chaos and birhythmicity originating from self-modulation of the inositol 1,4,5-trisphosphate signal in a model for intracellular Ca2+ oscillations // Bulletin of Mathematical Biology. 1999. Vol. 61. P. 507.

7. Kazantsev V.B., Vorobiev A.V. // Izvestiya VUZ. Applied Nonlinear Dynamics. 2009. Vol.17, №2. Ss.123-137. (In Russian).

8. Matrosov V.V., Gordleeva S.Yu., Kazantsev V.B. // Izvestiya VUZ. Applied Nonlinear Dynamics. 2012. Vol.20, №3. Ss.40-50. (In Russian).

9. Matrosov V.V., Kazantsev V.B. Bifurcation mechanisms of regular and chaotic network signaling in brain astrocytes // Chaos. 2011. Vol. 21, № 2. P. 023103.

10. Wu Y.W., Tang X., Arizono M., Bannai H., Shih P.Y., Dembitskaya Y., Kazantsev V., Tanaka M., Itohara S., Mikoshiba K., Semyanov A. Spatiotemporal calcium dynamics in single astrocytes and its modulation by neuronal activity // Cell Calcium. 2014. Vol. 55(2). P. 119-129.

11. 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.

12. Matrosov V.V. Dinamika nelineynyh sistem. Programmnyi complex dlya issledo-vanya nelineynyh dynamicheskih sistem s nepreryvnym vremenem. N.Novgorod: NNGU, 2002. 54 s. (In Russian).

13. Boldyreva N.V., Matrosov V.V. // Proceedings of the XVII-th Scientific Conference on Radiophysics dedicated to the 100th Anniversary of V.S. Troitskiy's Birth (Nizhni Novgorod, May 13-17, 2013). /Ed. by S.M. Grach, A.V. Yakimov. Nizhni Novgorod: UNN, 2013. P. 89-90.

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

REGULAR AND CHAOTIC OSCILLATIONS IN ASTROCYTE MODEL WITH REGULATION OF CALCIUM RELEASE KINETICS

V. V. Matrosov1, N. V. Boldyreva1, V.B. Kazantsev1'2

xNizhny Novgorod State University 2IPF RAN, Nizhny Novgorod

The dynamics of an astrocyte model is investigated. The astrocytes represent a type of glial cells regulating oscillations of major signaling cells, e.g. neurons. Subserved by complex molecular mechanisms the astrocytes generate calcium auto-oscillations which, in turn, are associated with the release of neuroactive chemicals into extracellular space. At variance with classical astrocyte models the three-component model considered takes into account a regulation of calcium release due to nonlinear dynamics of inositol-1,4,5 trisphosphate (IP3). Stability of steady states, generation of regular and chaotic oscillations were analysed using bifurcation analysis. In particular, it was found that nonlinear IP3 relaxation dynamics can lead to appearance of complex modes in astrocyte dynamics including multistability and pulse oscillations of complex shape.

Keywords: Astrocyte, calcium, generation, bifurcation, chaos.

Матросов Валерий Владимирович - родился в 1960 году, окончил факультет вычислительной математики и кибернетики Горьковского госуниверситета им. Н.И. Лобачевского в 1982 году по специальности прикладная математика. Защитил диссертации на соискание ученой степени кандидата физико-математических наук (1994), доктора физико-математических наук (2007), профессор. С 1999 года работает на радиофизическом факультете ННГУ, в настоящее время декан радиофизического факультета и заведующий кафедрой теории колебаний и автоматического регулирования. Область научных интересов - динамика нелинейных систем, динамический хаос, синхронизация и управление хаосом, математическое моделирование. Имеет более 100 публикаций в отечественных и зарубежных изданиях, трех монографий и трех учебных пособий (в соавторстве).

603950, Нижний Новгород, пр. Гагарина, 23

Нижегородский государственный университет им. Н.И. Лобачевского E-mail: matrosov@rf.unn.ru

Болдырева Наталья Викторовна - родилась в 1990 году, окончила магистратуру радиофизического факультета Нижегородского госуниверситета им. Н.И. Лобачевского по программе «Нелинейные колебания и волны» (2013). В настоящее время - аспирант кафедры теории колебаний и автоматического регулирования радиофизического факультета ННГУ. Область научных интересов - моделирование динамики нелинейных систем. Имеет 5 публикаций.

603950, Нижний Новгород, пр. Гагарина, 23

Нижегородский государственный университет им. Н.И. Лобачевского E-mail: B0-209@yandex.ru

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

603950, Нижний Новгород, пр. Гагарина, 23

Нижегородский государственный университет им. Н.И. Лобачевского 603950 Нижний Новгород, ул. Ульянова 46 Институт прикладной физики РАН E-mail: kazantsev@neuron.nnov.ru