Актуальность темы
Начиная с основополагающей работы А. Н. Колмогорова — И. Г. Петровского — Н. С. Пискунова, опубликованной в 1937 году, проблемам распространения самоподдерживающихся волн в активных средах было посвящено множество трудов. После публикации в 1951 году уравнений Ходжкин и Хаксли, описывающих распространение потенциала по мембране аксона кальмара, появились модели целого ряда активных сред. К ним относятся модели, описывающие протекание химических колебательных реакций в гомогенных и гетерогенных системах, в том числе реакций на каталитических поверхностях, модели динамики популяций, кинетические модели процессов свёртывания крови, модели фракционирования ионов и соединений различной хиральности на поверхности океана, и ряд других моделей, описывающих неравновесные явления в физике, химии и биологии.
Особое место в современной динамике активных сред занимает электрофизиология сердечной ткани. Последние тридцать лет в этой области доминирует гипотеза, согласно которой желудочковая фибрилляция, являющаяся причиной четверти всех внезапных сердечных смертей, возникает, как правило, в результате образования в стенках желудочков незатухающего самоподдерживающегося источника спиральных волн.
Основные свойства макроскопически протяженных автоволн, наблюдаемых во всех вышеперечисленных системах, хорошо известны. Это независимость стационарной скорости распространения такого рода автоволн от способа их инициации, взаимная аннигиляция автоволн при столкновении, способность к дифракции. Работы последних лет в этой области показали возможность существования, в возбудимых системах реакционно-диффузионного типа, автоволн в виде локализованных стационарных и нестационарных импульсов. На принципиальную возможность существования локализованных в пространстве, движущихся автоволновых объектов впервые указал А. Н. Заикин в 1993 году. Возник вопрос, является ли существование пространственно локализованных автоволновых решений исключительным случаем, свойственным только модели, использованной в оригинальной работе А. Н. Заикина, или же это достаточно общее свойство реакционно-диффузионных систем?
В работах последних лет А. С. Михайлова, X. Свинни, М. А. Цыганова, О. А. Морнева, Г. Р. Иваницкого, А. И. Лобанова, Т. К. Старожиловой и соавторов было показано, что пространственно локализованные объекты способны существовать в достаточно широком классе возбудимых систем. Можно сказать, что они свойственны основным базовым моделям возбудимых сред. Были тщательно исследованы механизмы зарождения, распространения и гибели пространственно локализованных объектов в однородных активных средах.
В самое последнее время заметно вырос научный интерес к вопросам управления динамикой возбудимых сред извне. Управление автоволновой динамикой активных сред представляет особый интерес в случае биофизически значимых слабовозбудимых сред, в которых возможно существование пространственно локализованных импульсов. Особый класс внешних управляющих воздействий представляют локальные, по сути точечные, неоднородности, внесение которых в активную среду способно существенно изменить целостную картину макроскопической автоволновой динамики. В данной работе систематическим образом исследуются механизмы трансформации автоволн в локально неоднородных активных средах.
Цель работы. Основной целью данной работы является исследование нелинейных эффектов, свойственных процессам распространения пространственно локализованных автоволн в активных средах при наличии в них неоднородностей. Центральными представлялись вопросы, каковы критические параметры локально невозбудимых участков, при которых существование автоволновых процессов возможно в принципе? Каковы основные качественно различные типы трансформаций автоволновых режимов при распространении автоволн в локально неоднородных средах?
В целях изучения проблем управления макроскопической автоволновой динамикой был поставлен вопрос, возможно ли с помощью создания локальных неоднородностей в активной среде обеспечить условия для асимметричной селекции правых и левых форм, в частности правои левовращающихся спиралей в активных средах?
Методы исследования. В работе использованы аналитические и численные математические методы. Компьютерное моделирование базовых моделей активных сред проводилось методами численного интегрирования систем дифференциальных уравнений. Программы написаны автором данной работы при помощи средства для разработки программного обеспечения Borland С++ Builder версии 6.0. Все графические пакеты отображения функций созданы автором.
Научная новизна работы. За исключением нескольких ранних работ 50-х — 60-х годов, выполненных Винером, Моу, Кринским и Заикиным, проблемы автоволновой динамики исследовались, как правило, для однородных активных сред. В связи с заметным расширением арсенала вычислительных компьютерных средств открылись новые, существенно более широкие возможности для исследования механизмов взаимодействия и трансформации автоволн. Следует надеяться, что проведенное в настоящей диссертации исследование распространения автоволн в локально неоднородных средах, хотя оно и основано, в значительной степени, на вычислительных экспериментах, представляет, тем не менее, достаточно общий научный интерес.
В частности, удалось обнаружить и исследовать своеобразный автоволновой туннельный эффект. Найдены условия трансформации автоволновых режимов в локализованные движущиеся возбуждения и спиральные волны. Обнаружены и исследованы взаимодействия частицеподобных автоволн. Для одной из основных базовых моделей возбудимой среды, учитывающей диффузию ингибитора, удалось обнаружить широкий спектр новых решений — стационарных и нестационарных структур, и построить соответствующие пространственно-временные картины. Очень важным является тот факт, что основные результаты, представленные в данной диссертации, были получены для одной из наиболее широко используемых базовых моделей возбудимой среды (модели Фитц-Хью — Нагумо). Это даёт основание утверждать, что полученные в работе выводы являются достаточно общими, применимыми к основным из известных активных сред.
Практическое значение работы. Нарушения процессов распространения автоволн в таких активных проводящих средах как сердечная, мышечная и нервная ткани могут быть губительны для организма. Углубление теоретического понимания механизмов распространения автоволн в локально неоднородных средах открывает новые перспективные пути для решения практических задач управления автоволновыми процессами в миокарде и других биологически значимых активных системах.
Структура работы. Диссертация посвящена исследованию процессов трансформации автоволн в неоднородных активных средах на примерах базовых моделей возбудимых сред. Работа состоит из пяти глав. В Главе I, посвящённой обзору литературы, рассматриваются существующие на сегодняшний день базовые модели разных типов активных сред — бистабильных, возбудимых и возбудимых с активным восстановлением, и обсуждаются их основные характеристики и свойства, а также рассматривается практическая актуальность решения соответствующих задач, как-то проблема фибрилляции или свёртывания крови. Во второй части первой главы даётся обзор публикаций последних десяти лет, в которых сообщается об открытии в разных моделях нетипичных по представлениям прошлых лет свойств автоволновых решений. К таким свойствам относится способность автоволн не аннигилировать при столкновении и возможность существования локализованных пространственных возбуждений. Обосновывается актуальность поставленных в работе проблем.
В Главе П рассмотрено распространение автоволн в одномерных локально неоднородных активных средах. Представлены результаты исследования распространения автоволн в неоднородных одномерных бистабильных и возбудимых средах. Рассмотрены среды, в которых неоднородность представлена в виде одиночного барьера, последовательности барьеров типа потенциальной решётки, а также асимметричного барьера. Найдены условия подбарьерного прохождения автоволн. Для возбудимых сред найдены условия возникновения источника периодической последовательности импульсов на неоднородном участке с повышенной возбудимостью.
В Главе Ш рассмотрено распространение автоволн в неоднородных двумерных средах. Найдены условия существования локализованных движущихся и покоящихся возбуждений. При исследовании взаимодействий локализованных возбуждений найдены режимы упругого отражения, а также режимы отражения под углом. При взаимодействии движущегося возбуждения со стационарным найдены условия трансформации локализованного возбуждения в спиральную волну. Рассмотрено распространение автоволнового фронта в среде с случайно разбросанными невозбудимыми зонами, размер которых в несколько раз меньше длины волны, и найден режим хаотического блуждания автоволнового фронта в подобным образом фрагментированной области. Также показано, что комбинированное препятствие, состоящее из невозбудимых зон и зон, позволяющих туннелирование, после прохождения через него плоского автоволнового фронта может стать источником спиральных волн.
В Главе IV исследуются механизмы двумерного структурообразования в слабовозбудимых локально неоднородных средах. Найдены режимы упругого столкновения частицеподобных импульсов, и режимы деления импульсов в рассмотренных возбудимых средах. Для двумерных активных сред показана возможность образования как стационарных, так и нестационарных динамических структур. Указаны пути дальнейшего исследования механизмов структурообразования в средах с локальными неоднородностями.
Глава V посвящена изучению возможных автоволновых механизмов в формировании режимов с различной динамической. хиральностью. В рамках исследованной базовой модели возбудимой среды продемонстрировано несколько способов трансформации плоского автоволнового фронта как в левую, так и в правую спиральную волну. Намечены дальнейшие пути исследования данной проблемы в связи с возможностью создания неоднородных поверхностей-фильтров, способных производить селекцию асимметричных форм. В свете полученных в Главе V результатов обсуждается общебиологическая значимость проблемы хиральности.
В Заключении обсуждаются полученные в работе результаты, их достоверность и возможная биофизическая значимость. Выражаются благодарности коллегам и соавторам.
В Выводах сформулированы основные выводы настоящей диссертационной работы.
В Приложении представлен анализ выбранной для изучения в данной работе модели активной среды, и показано, почему модель является базовой для изучения общих вопросов распространения и трансформации локализованных автоволн, содержащих неоднородности.
Основные выводы диссертации состоят в следующем:
1. Обнаружен эффект автоволнового туннелирования сквозь локальные препятствия в возбудимых средах. Установлено существование критического значения ширины симметричной невозбудимой зоны, при котором автоволна способна преодолеть препятствие и продолжить распространение. Обнаружено, что локальные асимметричные барьеры способны обеспечивать одностороннюю автоволновую проводимость в активной среде.
2. Обнаружен класс стационарных неоднородностей возбудимых сред, на которых может иметь место трансформация стационарных волн во вторичные генераторы асинхронных импульсов. Исследованы механизмы управления частотой вторичного источника автоволн и направлением генерируемых импульсов.
3. Исследованы свойства взаимодействия локализованных движущихся и покоящихся возбуждений в слабовозбудимых средах. Найдены условия, при которых при столкновении движущихся частицеподобных автоволн в активной среде могут иметь место эффекты их аннигиляции, отражения в перпендикулярном направлении и рассеяния под углом. При взаимодействии движущегося локализованного возбуждения с покоящимся обнаружен эффект трансформации квазичастиц в макроскопическую спиральную волну.
4. Установлен спектр частицеподобных пространственно-временных автоволновых решений модельной системы возбудимой среды. Исследованы режимы упругого взаимодействия частицеподобных импульсов и режим «эха», т. е. отражения импульсов от непроницаемых стенок. Выявлен механизм деления частицеподобных импульсов, лежащий в основе процессов установления хаотической импульсной активности в неравновесных активных системах. Обнаружены режимы «дышащих» пятен, рассеяния частицеподобных автоволн и их захвата препятствиями. Показана возможность управления пространственными макроскопическими структурами в активных средах, в том числе их подавлением за счёт введения локальных неоднородностей.
5. Установлено, что даже весьма кратковременные воздействия на систему, создающие градиент по степени ее возбудимости, могут трансформировать автоволновой фронт в правоили левовращающуюся спираль в зависимости от направления градиента. Показана возможность трансформации движущихся локализованных возбуждений в спиральные волны разной хиральности в результате столкновения со стационарным возмущением.
Автор выражает благодарность научным руководителям, Георгию Теодоровичу Гурия и Всеволоду Александровичу Твердислову, за проявленное внимание и благожелательность. Автор благодарит сотрудников кафедры биофизики физического факультета МГУ и сотрудников лаборатории криобиофизики Гематологического Научного Центра РАМН за плодотворные дискуссии, помощь и содействие.
ЗАКЛЮЧЕНИЕ
.
В данной работе мы показали, что свойства распространения автоволн в активных средах с локальными неоднородностями могут существенным образом отличаться от свойств распространения автоволн в однородных активных средах. Обнаруженный нами туннельный эффект для автоволн в активных средах свидетельствует о существовании критического значения локально невозбудимой зоны, при котором возможно автоволновое подбарьерное прохождение. Основным назначением автоволн в нервной, сердечной и мышечной тканях является передача импульса возбуждения, таким образом способность автоволны туннелировать через поражённые участки проводящей живой ткани представляется важным свойством автоволнового распространения импульса. Также нами было показано, что локальное асимметричное препятствие может создавать условия для односторонней проводимости импульсов, что представляет интерес в связи с организацией однонаправленных информационных потоков по нервным каналам.
Возникновение периодического источника импульсов в среде с локальной зоной повышенной возбудимости представляет собой яркий пример того, как можно существенно повлиять на развитие автоволнового режима точечным изменением свойств среды. Управление частотой возникающего источника может быть осуществлено при помощи создания локального участка пониженной возбудимости. В предельных случаях, когда зона пониженной возбудимости является непроходимой, мы можем получить источник, генерирующий импульсы только в одну сторону.
В двумерных средах наличие локальной неоднородности в среде может привести к уничтожению всего автоволнового фронта. Разрыв автоволнового фронта на препятствии в слабовозбудимых средах может стать критическим для существования автоволнового режима как такового. Оказалось, что в слабовозбудимых средах также возможно существование локализованных движущихся возбуждений, или автоволн-квазичастиц, которые получаются при прохождении плоского автоволнового фронта через щель в непроходимом барьере. Получившиеся локализованные возбуждения, или автоволны-квазичастицы, при взаимодействии друг с другом могут аннигилировать, переизлучаться в перпендикулярном направлении или расходится под углом. Кроме того, в слабовозбудимых средах могут существовать неподвижные локализованные возмущения, которые получаются при искусственном возмущении некоторого конечного участка среды. При взаимодействии движущихся и стационарных локализованных автоволн-квазичастиц был обнаружен эффект захвата движущейся автоволновой квазичастицы покоящейся с последующей трансформацией движущейся частицы в спиральную автоволну.
Явление трансформации плоского автоволнового фронта в спиральные волны было продемонстрировано на единичном препятствии, комбинирующем зоны полной непроходимости с зонами туннелирования. С другой стороны, было обнаружено, что локально непроходимые участки, размеры которых в несколько раз меньше длины волны, могут вызвать циркуляцию автоволнового возбуждения после прохождения плоского автоволнового фронта через определённым образом фрагментированную область.
При учёте диффузии ингибитора для базовой модели возбудимой среды было обнаружено разнообразие пространственно временных структур. При последовательном повышении соотношений коэффициентов диффузии активатора и ингибитора в одномерных средах были обнаружены режимы упругого столкновения импульсов, деления импульсов, режим «дышащих» пятен и стационарных структур. В двумерных средах обнаруженные пространственно-временные структуры соответствовали динамическому нестационарному режиму с делением, отражением и аннигиляцией пятен, или же формированию стационарных структур. Было показано, что стационарные структуры существенным образом зависят от начальных условий. Проблема формирования пространственно-временных структур в возбудимых средах, в которых диффузия ингибитора превышает диффузию активатора, при наличии локальных неоднородностей остаётся на сегодняшний день открытой. Были намечены пути дальнейшего исследования, и, в частности, показано, что внесение локальной неоднородности может быть критическим для формирования пространственной стационарной структуры в некоторой зоне исследуемой области. Каковы размеры критической области для динамически меняющихся структур, и насколько существенно влияние локальной неоднородности в этом случае — эти вопросы пока остаются открытыми.
Была затронута тема образования правых и левых спиральных волн в активных средах в связи с возможностью создания каталитических поверхностей-фильтров, осуществляющих отбор хиральных макроскопических структур в активных средах. Проблема хиральности в биологии вообще, и проблема образования и отбора хиральных молекулярных соединений в частности, весьма актуальна, однако в рамках настоящей работы мы не ставили целью найти связь между макроскопическим уровнем автоволновых процессов и молекулярным уровнем хиральных соединений. Было лишь продемонстрировано несколько способов образования правых и левых спиралей при помощи воздействия на систему путём кратковременного наложения внешнего поля и с помощью внесения локальной неоднородности в активную среду.
Таким образом, проведённое в данной работе исследование показывает, что в активных средах может иметь место широкий класс явлений трансформации автоволновых режимов на локальных неоднородностях. Трансформация автоволновых режимов может приводить к существенным качественным изменениям процессов распространения автоволн во всей среде. С точки зрения возможности управления автоволновыми процессами путём одиночного точечного воздействия на систему, изучение распространения автоволновых режимов в средах с локальными неоднородностями представляет огромный научный интерес. Автор надеется, что результаты, полученные в данной работе, смогут найти дальнейшее применение в электрофизиологии при изучении распространения волн возбуждения по нервной, мышечной и сердечной тканям, а также при моделировании широкого круга автоволновых процессов на разных уровнях биологической организации.