Диссертация (1143463), страница 22
Текст из файла (страница 22)
Обычно тета-гаммаколебания являются периодическими. Наконец, тета-ритм состоит из повторяющихся одиночных медленных спайков на почти спокойном фоне.Для медленной OLM-ячейки разница во временной динамике между тремя ритмами не столь выражены; здесь переключения между режимамиотражаются в разности сдвигов фаз между двумя клетками.Гамма-ритмЕсли значение проводимости равно нулю, то -клетки неоказывают ингибирующего действия на пирамидальную клетку; это соответствует удалению верхних соединительных дуг в схеме рис. 4.1.
Вэтой ситуации система совершает периодические колебания гамма-типа.Последовательные быстрые спайки происходят за определенный период,но точное повторение (другими словами, замыкание предельного циклав фазовом пространстве системы) требует нескольких спайков.В частности, малое значение может привести к асинхронностимежду медленными нейронами (рис. 4.2). Режим гамма сохраняется придостаточно низких ненулевых значениях , то есть при достаточнослабой связи между медленными и быстрой нейронами. Пример показанна рис.
4.3. При = 0.035 медленные нейроны колеблются в соотношении 1 : 1 (нижняя панель), но между максимумами спайков существуетсдвиг. В данном контексте естественно говорить о сдвиге фаз. Здесь фазазадавалась феноменологически: для каждого элемента приращение фазы равно 2 для каждого интервала между двумя последовательнымиположительными максимумами; между максимумами фаза была линейна интерполирована. Фазовый сдвиг между колебаниями, показаннымив нижней панели рисунке 4.3, близок к 2/3.Тета-гамма ритм155vP10-10100200t3004005000100200t300400500vL1,210-1Рис.
4.3:Гамма-ритм при =0.035, =0.01.Верхняя панель: динамика потенциала пирамидальной клетки. Нижняя панель: динамика потенциала OLM клеток;сплошная линия:1 ;пунктирная линия:1:1; сдвиг фаз близок к2 .OLM-клетки колеблятся в соотношении2/3В широком диапазоне значений выше 0.0362, соответствующий выбор начального условия приводит к другому ритмическому поведению: система проявляет тета-гамма колебания. На нижней панелирис. 4.4 показано, что тета-гамма режиму соответствует полная синхронизация между двумя медленными клетками.Тета ритмПри достаточно высоких значениях тета-гамма ритм заменяется тета-колебаниями.
Переключение между двумя режимами показанона рис. 4.5: динамика потенциала в начальном интервале времени имеетвид типичных колебаний тета-гамма, с пачками быстрых спайков, придостижении ≈ 200, динамика приобретает характерную форму тета156vP10-10100200t3004005000100200t300400500vL1,210-1Рис. 4.4:Тета-гамма ритм; =0.05.Верхняя панель: динамика потенциала пирамидальной клетки. Нижняя панель: динамика потенциала OLM клеток; сплошнаялиния:1 ;пунктирная линия:2 .ритма с равноудаленными спайками быстрой клетки. Надо отметить,что медленные нейроны, синфазно колеблясь на начальном этапе, “разъезжаются” при переключении в тета-гамма на : в тета-ритме клетки колеблются в противофазе.ГистерезисКак видно из предыдущего параграфа, в интервалах достаточнобольших значений , различные режимы сосуществуют как аттракторы систем (4.1, 4.2): при низких значениях из , в зависимости отначальных значений переменных, система может проявлять гамма-, атакже тета-гамма колебания, тогда как при средних и больших значениях существует гистерезис между тета-гамма- и чистые тета ритмы.157vP10-1-20100200t3004005000100200t300400500vL1,210-1-2Рис.
4.5:Распад тета-гамма ритма и включение тета-ритма. =2.274.Верхняяпанель: динамика потенциала пирамидальной клетки. Нижняя панель: динамика по клеток; сплошная линия: 1 ; пунктирная линия: 2 . Здесь синфазныеколебания клеток переключаются в антифазу сразу же после включения тетатенциаларежимаДля медленных -клеток, ни переход между гамма- и тета-гамма состояниями, ни дальнейший рост не влияет на период их колебаний. В противоположность, периодичность в паттернах быстрых клеток и изменяется довольно заметно. В таком контексте довольно удобноиспользовать характеристику, измеряемую как средняя величина интервала между спайками (ISI): среднее расстояние между положительнымимаксимумами колебаний потенциала клетки. На рис.
4.6 построена эта характеристика для колебаний мембранного потенциала в пирамидальнойклетки. Значение ISI увеличивается с 14.73 в гамма-ритме при = 0до 35.94 в тета-ритме при = 3. Аналогичная ситуация наблюдаетсяи для корзинчатой клетки.15830theta / gamma2030thetamagaminterspike intervaltheta100Рис. 4.6:theta/gamma20gamma1GLP2Длина интерспайкового интервала в зависимости от0.05GLP 0.1для различныхритмов. Сплошная линия: гамма-ритм. Точечная линия: тета-гамма ритм.
Пунктирная линия: тета-ритм. Вставка: переходы между режимами при низких значениях .При более детальном рассмотрении можно увидеть, что гамма- итета-ветви связаны; на самом деле они принадлежат к одному и тому женепрерывному семейству решений. Преобразование из гамма-тета-ритмапроисходит постепенно через промежуточные режимы, наблюдаемые вдиапазоне параметров 0.0724 < < 0.0734. Основные этапы этогопроцесса представлены на рис. 4.7, где показана динамика потенциалапирамидальной клетки.
Трансформации предшествует постепенная деформация гамма-паттерна: амплитуда двух из трех первоначально почтиравных пиков (см. рис. 4.3) уменьшается (рис. 4.7(а, b)). В результате периодический гамма-ритм приобретает характерную форму смешанногорежима колебаний [284] с чередующимися малыми подпороговыми и надпороговыми максимумами (порогом считается потенциал равный нулю).При = 0.0726 эти колебания приобретают форму сложнопериоиче159ского движения (с периодом 541.1), см. рис. 4.7(c).Небольшое увеличение приводит к частичному сглаживаниюподпороговых спайков (рис.
4.7(d)); следующее небольшое увеличениепараметра приводит к регулярному чередованию подпороговых и надпороговых спайков(рис. 4.7(е)). Фактически, это состояние уже оченьблизко к тета-ритму; дальнейшее увеличение приводит только к количественным изменениям: подпороговые колебания постепенно сглаживаются и исчезают, и устанавливается тета ритм. На фазовом портрете(см. рис. 4.8), подпороговые спайки представлены эпохи небольшими петлями траектории. В процессе сглаживания эти петли становятся меньше,превращаются в каспы (cusps) и исчезают, так что сохраняются толькобольшие петли (соответствующие одиночным пикам тета-ритма).Примечательно, что описанное преобразование изменяет разностьфаз между медленными нейронами: как уже упоминалось, в гамма-состоянии разница близка к 2/3, и близость к этому значению сохраняетсядо начала колебаний смешанного типа (mix-mode oscillations) на рис.
4.7(c). Как только появляется паттерн, похожий на тета-ритм (рис. 4.7(e)),разность фаз достигает значения, близкого к , и антифаза в колебаниях медленных клеток остается отличительной чертой тета состояния вовсем интервале существования ритма.Подобный процесс транформации происходит и на другой ветви решений, что соответствует состоянию тета-гамма. В этом состоянии пачки быстрых спайков разделяются очень малыми подпороговыми спайками(см.
рис. 4.4); последний соответствуют второстепенным петлям напроекциях соответствующих фазовых портретов. В процессе увеличенияпроводимости , форма притягивающей орбиты подвергается непрерывной деформации (аналогично показанному на рис. 4.8 для тета-ритма): второстепенный цикл постепенно уменьшается, сжимается и полно1602vp(a)0-202004006002004006002004006002004006002004006004006002vp(b)0-202vp(c)0-202vp(d)0-202vp(e)0-202vp(f)0-2Рис.
4.7:0200tТрансформация гамма в тета-ритм для мембранного потенциала пирамидальной клетки. (a) =0.0734;(f ) =0.07;(b) =0.0724;(c) =0.0726;(d) =0.0732;(e) =0.1000.стью исчезает. Этот процесс приводит к постепенному уменьшению иустранению одного из максимумов в временной структуре колебаний.Как видно на рис. 4.6, ветвь тета-гамма колебаний изолирована в161пространстве параметров. При более детальном раасмотрении его рождение при низких значениях начинается с очень короткого сегмента колебаний смешанного типа; они есть близко к значениям = 0.0362, ноуже в = 0.0363 устанавливается характерная форма периодическихколебаний тета-гамма.
На другом конце ветви в = 2.274 предельныйцикл, которому соответствует тета-гамма ритм, исчезает в результате бифуркации седло-узла.Влияние начальных условийДалее была исследована (численно) чувствительность отклика системы относительно изменения начальных условий. При низких (ниже0.036) значениях проводимости колебательное состояние гамма-типабыло достигнуто для всех исследованных начальных условий.
С другойстороны, при достаточно больших (выше 2.28) значениях система,похоже, обладает уникальным аттрактором- предельному циклу соответствующему тета-ритму. В интервале между этими значениями системаявляется бистабильной. Для дальнейшего исследования было взято значение = 0.8 из интервала, где различные колебания тета-типа сосуществуют с тета-гамма ритмом.Система была сканирована по сетке начальных условий для переменных и от -2 до 2: эти диапазоны соответствуют максимальномудиапазону колебаний гамма и тета-гамма для пирамидальных клеток.Исходный размер сетки был взят как 0.01 × 0.01 и вычисления проводились с переменным шагом по сетке. Оказалось, что выбор конечногоаттрактора зависит в основном от начального значения мембранного потенциала и мембранной переменной медленных клеток.Чтобы графически проиллюстрировать форму границы между бас162up10-1-2Рис.
4.8:-1vp01Эволюция аттрактора в фазовом пространстве (проекция для потенциалаи мембранной переменной пирамидальной клетки). Маленькие петли, соответствующие подпороговым колебаниям уменьшаются, певращаются в выступы и исчезают.Сплошная линия: =0.09, пунктирная линия: =0.2, точечная линия: =0.4.сейнами аттракторов в 12-мерном пространстве начальных значений, выберем две характеристики пересечения этой границы с координатнымиплоскостями. Первый график (левая панель рис. 4.9) является плоскостью на который все координаты, кроме (переменная потенциала пирамидальной ячейки) и 1 (мембранная переменная первой медленнойклетки) равны нулю.