Исследование процесса распространения возбуждения в среде, имитирующей кластерную структуру миокарда тема диссертации и автореферата по ВАК РФ 03.00.13, кандидат биологических наук Коновалова, Ирина Николаевна

• Коновалова, Ирина Николаевна
• кандидат биологических наук
• 2000

Оглавление диссертации кандидат биологических наук Коновалова, Ирина Николаевна

ВВЕДЕНИЕ.

ГЛАВА 1. ЛИТЕРАТУРНЫЙ ОБЗОР

1.1. Строение и электрические характеристики миокарда.

1.2. Экспериментальные данные о фибрилляции сердечной мышцы

1.3. Обзор математических моделей возбудимых сред.

ГЛАВА 2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА РАСПРОСТРАНЕНИЯ ВОЗБУЖДЕНИЯ В СЕРДЕЧНОЙ ТКАНИ.

2.1. Назначение и основные особенности модели.

2.2. формальная модель кластера, сердечной ткани и процесса распространения возбуждения.

2.2.1. Кластер.

2.2.2. Геометрия проводящей среды.

2.2.3. Вероятности перехода между сврЩ&яшялш.

2.2.4. Моделирование процесса распространения возбуждения.

2.3. Количественное выражение уровня фибрилляции.

2.4. Пространственная неоднородность.

ГЛАВА 3. КОМПЬЮТЕРНАЯ ПРОГРАММА ДЛЯ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ ПРОЦЕССА РАСПРОСТРАНЕНИЯ ВОЗБУЖДЕНИЯ В СЕРДЕЧНОЙ ТКАНИ.

3.1. Принцип построения, возможности и особенности программы

3.1.1. Назначение программы.

3.1.2. Дополнения к математической модели.

3.2. Описание программы.

3.2.1. Раздел главного меню ФАЙЛ.

3.2.2. Раздел главного меню ПАРАМЕТРЫ.

3.2.3. Раздел главного меню ЭКСПЕРИМЕНТ.

ГЛАВА 4. НЕКОТОРЫЕ РЕЗУЛЬТАТЫ, ПОЛУЧЕННЫЕ С ПОМОЩЬЮ КОМПЬЮТЕРНОЙ МОДЕЛИ И ИХ ЭКСПЕРИМЕНТАЛЬНАЯ ПРОВЕРКА.

4.1. Экспериментальные данные о трансмембранных потенциалах кардиомиоцитов сердца в норме и при фибрилляции.

4.1.1. Общая характеристика опытов на интактном сердце.

4.1.2. Регистрация трансмембранных потенциалов кардиомиоцитов интактного сердца.

4.1.3. Электрическая активность клеток сердца морских свинок.

4.1.4. Электрическая активность кардиомиоцитов сократительного миокарда желудочков интактного сердца крысы.

4.2. Компьютерные эксперименты с использованием опции Наблюдение.

4.2.1. Нарушения в единичном кластере.

4.2.2. Нарушения в одном слое.

4.2.3. Нарушения в одной трубке.

4.2.4. Пример пространственной неоднородности сложной геометрической формы.

4.2.5. Наличие случайного фактора.

4.2.6. Минимальный объем среды, в котором возможна фибрилляция

4.2.7. Воздействие одиночного импульса.

4.3. Компьютерные эксперименты с использованием опции Гетерогенность.

4.3.1. Зависимость уровня фибрилляции от периода рефрактерности

4.3.2. Зависимость уровня фибрилляции от неоднородности по периоду рефрактерности.

4.3.3. Зависимость уровня фибрилляции от порога возбуждения.

4.3.4. Зависимость уровня фибрилляции от случайного фактора (от степени неоднородности по порогу возбуждения).

Введение диссертации (часть автореферата) на тему «Исследование процесса распространения возбуждения в среде, имитирующей кластерную структуру миокарда»

Актуальность.

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

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

Указанная модель была усовершенствована в работах И.М. Гельфанда и M.JI. Цетлина [25], а также В.И. Кринского [43]. В этих работах уже были рассмотрены процессы возникновения и прекращения фибрилляции.

В 70-х годах наблюдался большой поток работ по моделированию процессов распространения возбуждения в сердечной ткани (подробнее об этом в главе 1), однако, затем интерес к этому подходу заметно уменьшился.

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

1) Представления о физиологической природе указанных явлений, а также о механизме передачи возбуждения в миокарде претерпели существенные изменения. На передний план вышла и получила экспериментальные подтверждения гипотеза о том, что сердечная ткань представляет собой совокупность кластеров. Кластер, являясь группой кардиомиоцитов, выступает в процессе распространения возбуждения как единое целое, а ансамбль кластеров играет в этом процессе роль структурно и функционально образующего для сердечной ткани. Идея кластерной организации миокарда впервые наиболее четко была сформулирована еще в 1972 г. группой сотрудников профессора Г.И. Косицкого в дальнейшем активно разрабатывалась на кафедре нормальной физиологии РГМУ.

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

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

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

Итак, проблема изучения аритмий и фибрилляции сердца продолжает оставаться одной из основных как в теоретической, так и в клинической кардиологии и требует комплексного подхода и, как следствие, привлечения всех доступных современных методов исследования. Достижения последних лет в области математического моделирования, подкрепленные бурным развитием компьютерной техники, подталкивают обратиться к этим методам для продвижения в области изучения процессов распространения возбуждения в сердечной ткани. В частности, методы математического и компьютерного моделирования, которым уделялось недостаточно внимания, могли бы помочь охватить совокупность новых физиологических представлений и объяснить накопленные в последнее время экспериментальные факты. Такой новой физиологической теорией явилась кластерная гетерогенная концепция фибрилляции, которая в работе стала основанием и подверглась анализу с использованием новейших средств имитационного моделирования. Сказанное обосновывает актуальность работы, лежащей на стыке биологии, физиологии и математического моделирования.

Цель исследования.

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

Задачи исследования.

1) На базе электрофизиологических исследований, обосновывающих кластерную организацию миокарда, создать математическую модель сердечной ткани.

2) Реализовать ее в виде компьютерной имитационной модели.

3) Провести эксперименты с компьютерной моделью для выяснения ее работоспособности, диапазона сферы использования, а также для сопоставления результатов моделирования процессов с результатами реального эксперимента.

Научная новизна.

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

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

Введено понятие уровня фибрилляции, которое положено в основу количественного анализа процесса распространения возбуждения.

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

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

Практическая значимость.

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

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

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

Апробаиия диссертационного материала и внедрение результатов исследования в практику.

Основные положения диссертационной работы были доложены на: 4-ом международном симпозиуме по сравнительной электрокардиологии. ( Сыктывкар, 1997), Международной научно-технической и Российской научной школы (Москва-Сочи, 1998), Международная конференция "Компьютерная электрокардиография на рубеже столетий" (Москва, апрель 1999), XXVI Международном конгрессе по электрокардиологии (Сыктывкар, 1999), Международной научно-техническая конференция и российская научная школа молодых ученых и специалистов "Системные проблемы качества, математического моделирования и информационных технологий" (Москва-Сочи, сентябрь 1999).

Содержание работы было также доложено и обсуждено на совместной научно-практической конференции кафедры нормальной физиологии с курсом физиологии медико-биологического факультета и кафедры общей патологии РГМУ.

По теме диссертации имеются следующие публикации:

1. Коновалова И.Н. Механизмы гормонального контроля фибрилляции и дефибрилляции сердца. // В сб.: Интеграция механизмов регуляции висцеральных функций. Материалы симпозиума. - Краснодар, 1996, - С. 32. (Соавт.: Кобрин В.И., Маноах М., Коврижкин А.А, Порман Е.Е.)

2. Коновалова И.Н. Компьютерная модель кластерной структуры миокарда. // В сб.: Сравнительная электрокардиология. 4-й международный симпозиум по сравнительной электрокардиологии. -Сыктывкар, 1997. - С. 44. (Соавт.: Кобрин В.И.)

3. Коновалова И.Н. Имитационная модель распространения возбуждения в сердечной ткани. // В сб.: Системные проблемы надежности, математического моделирования и информационных технологий. Материалы Международной научно-технической и

Российской научной школы. - Москва-Сочи, 1998. - С. 93. (Соавт.: Кобрин В.И.)

4. Коновалова И.Н. Компьютерная модель для изучения распространения возбуждения в сердечной ткани. // Военная кибернетика. - № 3 (2). - 1998. - С. 23-25. (Соавт.: Кобрин В.И.)

5. Коновалова И.Н. Компьютерное моделирование явлений фибрилляции в сердечной ткани. // В сб.: Компьютерная электрокардиография на рубеже столетий. - Москва, 1999. - С. 283. (Соавт.: Кобрин В.И.)

6. Коновалова И.Н. Компьютерная модель для изучения распространения возбуждения в сердечной ткани. // Информатика -машиностроение. - Т. 1(23). - 1999. - С. 38-43. (Соавт.: Кобрин В.И.)

7. Konovalova I. The imitation model of stimulation spreading in the heart tissue. // Modern Electrocardiology. XXVI International Congress on Electrocardiology. - Syktyvkar, 1999. - P. 32. (Co-author: Kobrin V.)

8. Коновалова И.Н. Компьютерная модель как инструмент экспериментального исследования фибрилляции сердечной ткани. // В сб.: Системные проблемы надежности, математического моделирования и информационных технологий. - Москва-Сочи, 1999. - 79-80. (Соавт.: Кобрин В.И.)

Модель используется при проведении практических занятий на кафедре физиологии по модулю "сердечно-сосудистая система ГЕА им. Маймонида.

Объем и структура работы.

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

Заключение диссертации по теме «Физиология», Коновалова, Ирина Николаевна

Выводы

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

2. В среде Delphi 3.2 реализована компьютерная версия математической модели, позволяющая исследовать процессы распространения возбуждения в трехмерной среде, максимально приближенной к реальному миокарду.

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

4. Имитационная модель позволяет моделировать процесс нарушения проведения возбуждения как между одиночными кластерами, так и группами кластеров. Показано, что для возникновения процесса фибрилляции достаточно нарушения функции одного кластера.

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

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

Заключение

Цель работы, заключавшаяся в построении и реализации модели ткани миокарда, достигалась по следующей схеме: теоретические и экспериментальные предпосылки - математическая модель компьютерная программа - эксперименты.

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

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

Остановимся на тех предположениях относительно структурно-функциональных свойств сердечной ткани, которые непосредственно нашли отражение в математической модели. а) Пространственно-геометрическая организация.

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

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

Механизм возбуждения кластера подобен механизму возбуждения клетки: кластер переходит в возбужденное состояние, если его потенциал повышается до определенного порогового значения. Это повышение в реальной ткани либо спонтанно, либо вызвано внешним воздействием. В модели основным параметром, регулирующим переход в возбужденное состояние, является "уровень воздействия" со стороны окружающих кластеров, выраженный в условных единицах, а спонтанное возбуждение, если его желательно учитывать, может быть задано специальным "стохастическим" коэффициентом. Для того, чтобы хотя бы отчасти отразить "составное" устройство кластера и, следовательно, его не вполне совпадающее с отдельной клеткой поведение, в модель возбуждения кластера привносится еще один, дополнительный, случайный фактор.

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

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

В работе этому вопросу уделяется первостепенное значение, как определяющему фактору в явлении фибрилляции. Затрагиваются различные аспекты гетерогенности проводящей среды. Все они носят пространственный характер, в то время как неоднородность во времени не рассматривается. Структурная неоднородность, то есть неоднородность по взаимному расположению кластеров, также не рассматривается. К числу учитываемых признаков гетерогенности относится, прежде всего, индивидуальный характер изменения потенциала действия кластера, в частности, период рефрактерности и амплитуда. Это позволяет моделировать среды из кластеров, разнотипных по составу клеток. В модели допускается также различие в механизме возбуждения кластеров. Этот механизм моделируется, главным образом, двумя параметрами: уровнем возбуждения (модельный аналог порогового потенциала) и случайным фактором. Оба параметра могут быть выбраны индивидуально для всех кластеров, а их различия моделируют нарушения во взаимодействии. Гетерогенность во взаимодействии кластеров может также выражаться в индивидуальном характере реакции на влияние соседних кластеров. Наконец, в модели предусмотрены специальные коэффициенты, которые, оставляя неизменным качественный характер межкластерного взаимодействия, ослабляют или усиливают взаимное влияние для заданной пары соседних кластеров. д) Процесс распространения возбуждения.

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

Изложенный комплекс предположений реализуется в виде математической модели сердечной ткани как среды, проводящей импульсы возбуждения. Модель является дискретной, в частности, для описания процесса распространения возбуждения используется дискретное время. Кластер моделируется как объект, который в каждый момент времени может находиться в одном из конечного множества состояний. Среди возможных состояний различаются состояние покоя (диастолы), а также состояния, моделирующие фазы деполяризации и реполяризации. Количество состояний моделирует период рефрактерности для реального объекта и является регулируемым параметром модели, что позволяет с желаемой точностью приближать заданный (вообще говоря, непрерывный) характер изменения трансмембранного потенциала. Переходы между состояниями кластера обусловлены заданным (вероятностным) механизмом. В этом механизме выделяется следующий основной момент: кластер может перейти в состояние возбуждения (деполяризации) только из состояния покоя (или из близкого к нему состояния реполяризации) и только при условии, что уровень воздействия на него со стороны близлежащих кластеров достиг определенного, характерного для данного кластера, порогового значения. Действие указанного правила может быть "размыто", благодаря привнесению случайности в алгоритм перехода. В частности, кластер может возбуждаться, когда воздействие на него меньше порогового и, наоборот, оставаться в состоянии покоя, даже если порог превышен. Случайный фактор (также как и пороговый уровень возбуждения) является управляемым параметром модели, и может быть увеличен, либо уменьшен (устранен).

Размеры среды также являются параметрами модели.

С математической точки зрения модель представляет собой однородную конечную цепь Маркова, состояния которой - векторы с компонентами, означающими состояния кластеров. Хотя существует достаточно развитый аппарат для теоретического исследования таких объектов, в работе не делается попыток аналитического изучения модели. Формальная модель использована для разработки ее компьютерной реализации в виде программы, имитирующей процесс распространения возбуждения в ткани миокарда и предоставляющей достаточно широкие возможности для экспериментальной деятельности.

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

Еще одно важное назначение созданной модели заключается в том, что она дает инструмент для количественной оценки в экспериментальных исследованиях зависимости уровня фибрилляции от степени гетерогенности среды. В работе вводится математическое определение "уровня фибрилляции". Это определение базируется на следующих наглядных соображениях. Идеализированная проводящая среда, представленная моделью, в том случае, если она полностью однородна, обладает тем свойством, что в каждой своей точке воспроизводит (с некоторым запаздыванием, но без искажений) входной внешний импульс возбуждения. Всякое нарушение этого свойства (вызванное нарушением взаимодействия кластеров и, следовательно, неоднородностью среды) можно трактовать как "нарушение ритма", или "фибрилляцию". Вводя подходящим образом меру отклонения от "правильного поведения", и применяя усреднение по объему и по времени, можно получить количественное выражение уровня фибрилляции для произвольного фрагмента среды. Эта характеристика с математической точки зрения является функционалом предельного типа на траекториях упомянутой марковской цепи и зависит от весьма большого количества аргументов -параметров модели. В работе изучается зависимость уровня фибрилляции от четырех "усредненных" параметров: среднего периода рефрактерности, степени неоднородности по периоду рефрактерности, среднему порогу возбуждения и степени неоднородности по порогу возбуждения. Эти величины получают в работе точное определение. Программная реализация предоставляет возможность получать оценки уровня фибрилляции заданного фрагмента среды в любой области указанного четырехмерного пространства аргументов, а также автоматически строить графики зависимости уровня фибрилляции от любого из четырех перечисленных параметров при фиксированных остальных трех.

Работа содержит также экспериментальную часть.

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

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

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

Перспективы развития модели.

1) Более углубленное отражение свойств сердечной ткани. Существенные возможности для этого содержатся в модели уже в ее нынешнем виде. Прежде всего, в связи с этим отметим возможность моделировать произвольную (а не только гексагональную) геометрическую конфигурацию среды, в том числе не обладающую пространственной симметрией. Большую степень свободы предоставляет также выбор функции g, моделирующей трансмембранный потенциал (п.2.2.1), и коэффициентов С/, влияющих на межкластерное взаимодействие (п. 2.2.3). (В имеющейся программной реализации эти параметры зафиксированы; однако не представляет особого труда видоизменить программу так, чтобы сделать их выбор свободным.)

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

2) Усовершенствование программной реализации. Имеется в виду, прежде всего, применение более совершенной трехмерной графики, а также другое развитие интерфейса. Это способствовало бы большему удобству в использовании данного продукта и, следовательно, более широкому практическому внедрению.

3) Изучение и оптимизация с помощью модели воздействий, приводящих к прекращению фибрилляции. Такая постановка задачи представляется новой и перспективной. Кратко схема подобного исследования могла бы выглядеть следующим образом. Механизмы, которые реально используются в медицинской практике для компенсации нарушений ритма (фармакологические, электрофизические и пр.), должны быть переведены в формальной модели на язык абстрактных управлений, которые определенным образом влияют на параметры кластеров (период рефрактерности, порог возбуждения и т. д.). После этого ставится задача оптимального управления с целью нахождения стратегии, минимизирующей в кратчайший срок нарушения ритма. Решение задачи вновь формулируется в терминах реальных клинических мероприятий.

Список литературы диссертационного исследования кандидат биологических наук Коновалова, Ирина Николаевна, 2000 год

1. Аршавский Ю.И., Беркинблит М.Б., Ковалев С.А. Периодическая трансформация режима в одиночных нервных волокнах // Биофизика. -1962. - Т.7. - С. 449.

2. Аршавский Ю.И., Беркинблит М.Б., Ковалев С.А., Смолянинов В.В. // Сб. "Модели структурно-функциональной организации некоторых биологических систем". С. 112.

3. Аршавский Ю.И., Беркинблит М.Б., Ковалев С.А. Место возникновения трансформации ритма в нервном волокне с искусственно созданной неоднородностью // Биофизика. 1962. - Т. 7. - С. 619.

4. Аршавский Ю.И., Беркинблит М.Б., Ковалев С.А., Чайлахян Л.М. // Биофизика. 1964. - Т. 9. - С. 634.

5. Аршавский Ю.И., Беркинблит М.Б., Ковалев С.А., Чайлахян Л.М. // Биофизика. 1964. - Т. 9. - С. 365.

6. Асланиди О.В., Морнев О.А. // Биофизика. 1996. - Т. 41. - С.953.

7. Балаховский И.С. Некоторые режимы движения возбуждения в идеальной возбудимой ткани // Биофизика. 1965. - Т. 10. - С. 1063.7 8. Березовский В. А., Колотилов Н.Н. Биофизические характеристики тканей человека. Справочник. Киев, 1990.

8. Беркинблит М.Б., Божекова В.П., Бойцова Л.Ю. и др. Высокопроницаемые контактные мембраны. М.: Наука, 1981.

9. Беркинблит М.Б., Дунин-Барховский В.Л. // Биофизика. 1964. -Т. 14. - С. 324.

10. Беркинблит М.Б., Ковалев С.А., Смолянинов В.В., Чайлахян Л.М. // Сб. "Модели структурно-функциональной организации некоторых биологических систем". Институт биологической физики АН СССР. -М.: Наука. - 1966.

11. Беркинблит М.Б., Розенштраух Л.Б., Чайлахян JI.M. Электрофизиология миокарда. Руководство по кардиологии. Под ред. Чазова Е.И. М.: Медицина, 1982. - Т.1.

12. Беркинблит М.Б., Ковалев С.А., Смолянинов В.В., Чайлахян JI.M. Определение основных электрических характеристик миокарда желудочка лягушки // Биофизика. 1965. - T.l 1, № 5. - С. 861-867.

13. Беркинблит М.Б., Ковалев С.А., Смолянинов В.В., Чайлахян JI.M. Электрическая структура миокардиальной ткани // ДАН СССР. -1965. Т. 163. - С.741-744.

14. Беркинблит М.Б., Ковалев С.А., Чайлахян JI.M., Введенская Н.Б., Холопов А.В., Дудзльичус И. Исследования распространения возбуждения в волокнах Пуркинье сердца на математической модели // Биофизика. 1970. - Т.15. - С. 521.

15. Беркинблит М.Б., Ковалев С.А., Чайлахян JI.M., Введенская Н.Б., Холопов А.В., Гнеденко Л.С. // Биофизика. 1971. - Т. 16. - С. 103.

16. Бондарев А.П., Кобрин В.И., Павлов Б.Н. Усилитель для микроэлектродных исследований // Физиол. ж-л СССР. 1980. - Т.66. -№4. - С.597-599.

17. Бриккер В.Н. Материалы 1 Всероссийского съезда кардиологов. -Воронеж. 1968. - С. 197.

18. Быховский М.Л., Коротков А.Д., Русанов A.M., Ходоров Б.И. // Биофизика. 1966. - Т.П.

19. Букаускас Ф.Р. Межклеточное взаимодействие и его роль в процессах распространения возбуждения. Дис. канд. мед. наук. Каунас: мед. институт, 1975.

20. Вайдман С. Функциональные соединения между сердечными клетками. // Проблемы общей и клинической физиологии сердечнососудистой системы. Киев: Наукова думка. - 1976. - С. 43-49.

21. Вайсбейн С.Г. Неотложные состояния в клинике внутренних болезней. М.: Медицина, 1966.

22. Гельфанд И.М., Ковалев С.А., Чайлахян JI.M. Внутриклеточное раздражение различных отделов сердца лягушки. // ДАН СССР. 1963. -Т. 148(4). - С. 973-976.

23. Гельфанд И.М., Фомин С.В. Вариационное исчисление. М.: Наука, 1961.

24. Гельфанд И.М., Цетлин M.JI. О континуальных моделях управляющих систем. // ДАН СССР. 1966. - Т. 131(6). - С. 1242.

25. Глаголева В.В., Чечулин Ю.С. Ультраструктурная организация миокарда. Физиология сердца. Руководство по физиологии. JL: Наука, 1980.

26. Гофман Б., Крейнфилд П. Электрофизиология сердца. М.: Иностранная литература, 1962. 390 с.

27. Дощинин B.JI. Лечение нарушений сердечного ритма современными средствами. Автореферат дисс. к.м.н. Москва, 1970.

28. ЖаботинскийА.М. Концентрационные автоколебания. М.: Наука, 1974. 180 с.

29. Зыков B.C. Моделирование волновых процессов в возбудимых средах. М.: Наука, 1984. 165 с.

30. Катц Б. Нерв, мышца и синапс. М.: Мир, 1968.

31. Кобрин В.И., Алабовский В.В., Олейников О.Д., Алипов Н.Н. Электрическая активность сердца и сократимость миокарда при изменении внеклеточной концентрации натрия. // Физиол. жур. СССР им. И.М. Сеченова. 1988. - Т. 74, № 9. - С. 1257-1262.

32. Кобрин В.И. Асинхронизм возбуждения клеток миокарда в онто- и филогенезе // Дисс.на соиск.уч.ст.д.м.н., М.: 1989. 321с.

33. Кобрин В.И. Гетерогенность миокарда и аритмии сердца. // Успехи физиологических наук. 1993. - Т. 24, № 4.

34. Колмогоров А.Н., Петровский И.Г., Пискунов Н.С. Бюллетень МГУ, А, 1, 6, 1973.

35. Кринский В.И. // Биофизика. 1964. - Т.9, №2.

36. Кринский В.И. // ДАН ССР. 1964. - Т. 156, № 6.

37. Кринский В.И. Распространение возбуждения в неоднородной среде (режимы аналогичные фибриляции сердца). // Биофизика. 1966. -Т.Н.

38. Кринский В.И., Холопов А.В. "Эхо" возбудимой ткани. // Биофизика. -1967. Т. 12. - С. 524.

39. Кринский В.И., Холопов А.В., Проведение импульсов в возбудимой ткани с непрерывно распределенной рефрактерностью. // Биофизика. 1967. - Т. 12. - С. 669.

40. Кринский В.И., Холопов А.В. // Биофизика. -1967. Т. 12, № 4.

41. Кринский В.И., Холопов А.В. Тезисы докладов на Всесоюзном межвузовском симпозиуме по прикладной математике и кибернетике. -Горький. 1967. - С. 192.

42. Кринский В.И. Фибрилляция в возбудимых средах. // Сб. "Проблемы кибернетики". М. - 1968. - Т.20. - С. 59.

43. Кринский В.И., Фомин С.В., Холопов А.В. О критической массе фибриляции. // Биофизика. 1968. - Т. 12. - С. 908.

44. Кринский В.И. // Сб. "Мат. Модели биологических систем". -М.: Наука. -1971. С. 41-51.

45. Кринский В.И., Крейцер Г.П., Перцов А.М„ Кукушкин Н.И., Сергеев Ю.В. // Сб. "Модели структурно-функциональной организации биологических систем". Москва-Дубна. - 1972. - С. 70.

46. Кринский В.И., Фомин С.В., Холопов А.В. // Биофизика. 1967. -Т. 12, №5.

47. Кринский В.И. и др. Эволюция автоволновых вихрей. (Волны в сердце). М., 1986.

48. Кринский В.И. Возбудимые среды. Волны и автоколебательные процессы // Сб. "Математическая биофизика клетки". М.: Наука. - 1978. -С. 141-201.

49. Кукушкин Н.И., Саксон М.Е. // Биофизика. 1971. - Т. 16.1. С.135.

50. Лукашевич И.П. // Биофизика. 1963. - Т. 8. - С. 715.

51. Малосолов А.Н., Христалюбова Н.Б., Михайлова Н.Б. Электронная микроскопия сердечной мышцы. // ДАН СССР. Т. 119 (1). -Р. 168.

52. Морнев О.А., Асланиди О.В., Алиев P.P., Чайлахян JI.M. // ДАН. 1996. - Т. 347, № 1. - С. 123-125.

53. Морнев О.А., Асланиди О.В. Чайлахян JI.M. Солитонный режим в системе уравнений Фитцхью-Нагумо: динамика вращающейся спиральной волны. // ДАН. 1996. - Т. 353, № 5. - С. 682-686.

54. Модели структурно-функциональной организации некоторых биологических систем. // Сб. М.: Наука. - 1966.

55. Перцов А.М., Кринский В.И. Передача возбуждения в электрически связанных клетках: мембранный механизм задержек и сердечный аритмии. // Биофизика. 1973. - Т. 18, вып. 4.

56. Розенштраух Л.В. Распространение возбуждения по различным структурам сердца. Физиология сердца. Руководство по физиологии. Л.: Наука, 1980.

57. Розенштраух Л.В., Холопов А.В. Механизм возникновения нейрогенных аритмий (мат. модель, основанная на данных микроэлектродного исследования // Сб. "Ритм сердца в норме и паталогии". Вильнюс, 1970.

58. Розенштраух Л.В., Холопов А.В. Анализ механизма поддержания фибриляции. // Гл. информационные материалы. АН СССР. Т. 9-10. - С. 35.

59. Розенштраух Л.В., Холопов А.В., Юшманов А.В. // Биофизика. -1972. Т. 15. - С. 690.

60. Румянцев П.П. Кардиомиоциты в процессах репродукции, дифференцировки и регенрации. М.: Наука, 1982. - 288 с.

61. Смолянинов В.В. // Сб. "Модели структурно-функциональной организации некоторых биологических систем". М.: Наука, 1966. - 112с.

62. Смолянинов В.В. Теория возбудимых тканей. М., 1971.

63. Смолянинов В.В. Теория синцитиальных тканей. Пущино,1972.

64. Смолянинов В.В. Математические модели биологических тканей. М.: Наука, 1980. - 366 с.

65. Смолянинов В.В. // Биофизика. 1966. - Т. 11. - С. 337.

66. Смолянинов В.В. // Биофизика. 1968. - Т. 13. - С. 492.

67. Смолянинов В.В. // Биофизика. 1969. - Т. 14. - С. 336.

68. Телеснин В.Р. Изд. ВУЗов. Т. 6, № 3. 1963.

69. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. М. : Наука, 1967.

70. Ходоров Б.И., Биковский М.Л., Русанов A.M., Коротков А.Д. // Биофизика. 1966. - Т. 11. - С. 1049.

71. Ходоров Б.И. // Биофизика. 1967. - Т. 12. - С. 1050.

72. Ходоров Б.И. Проблема возбудимости. Л.: Медицина, 1969.

73. Ходоров Б.И., Тимин Е.М. // Биофизика. 1970. - Т. 15. С. 503.

74. Ходоров Б.И., Тимин Е.М. // Биофизика. 1971. - Т. 16. С. 494.

75. Ходоров Б.И., Позин Н.В., Шмелев Л.А. // Биофизика. 1971. Т. 16. С. 95.

76. Ходжкин А.К. Нервный импульс. М.: Мир, 1965.

77. Холопов А.В. // Биофизика. 1968. - Т. 13. С. 670.

78. Цукерман Б.М. Электроимпульсная терапия ритма сердца. Дисс. д.б.н, М., 1971.

79. Черненко Ю.П. Морфофункциональная организация межклеточных контактов. Киев: Здоров'я, 1982. - С. 9-24.

80. Шидловский В.А., Дауэ Н.С., Особенности гистоархитектоники венозного синуса и предсердий. // Сб. "Материалы по экспериментально-клинической электрокардиографии". -М., 1953. С. 5-17.

81. Ватт L., Dewey М.М., Berger W. Propagation of action potentials and structure of nexus in cardiac muscle. // J. Gen. Physiol. 1965. - V. 48. -P.797.

82. Barr L. Electrical transmission between the cells of vertebrate cardiac muscle. Comparative Physiology of the heart: Ourrent Trends. Ed. McCann F.V., Basel, Stuttgart, 1969. - P. 102-110.

83. Brooks C.M., Hoffinan B.F., Suckling E.E., Orias E. Excitability of the heart. Grune and Stratton, N. Y., London, 1955.

84. Coolly J., Dodge F. // Biophys, J. 1966. V. 6. - P. 583.

85. Dant J. The passive electrical properties of guinea-pig ventricular muscle as examined with a voltage-clamp techniques. // J. Physiol. (Engl.). -1982. V.330. - P.221-242.

86. F. Dodge. // Thesis of Rochfeller Institute. N.Y., 1963.

87. Draper M.H., Meja-Tu. Aconparisson of the conduction velocity in cardiac tissues of various mammals. M., 1959.

88. Fawcett D., McNutt N. The ultrastructure of cat myocardium. Ventricular papillary muscle. // J. Cell Biol. 1969. - V. 41, N 1. - P. 1-20.

89. De Felice J.L. Challice C.E. Anathomical and ultrastructural study of the electrophysiological atrioventricular node of the rabbit. // Circulat. Res.1969. V.24, N 3. - P. 457-474.

90. Fitz Hugh R. // Bull Math Biophys. 1955. - V. 17. - P. 257.

91. Fozzgard H.A. Membrane capacity of the cardiac Purkinje fiber. // J. Physiol. 1966, - V. 182. - P. 225-267.

92. B. Frankenheuser , A.F. Huxley. // J. Phisiol. 1964. V. 171.1. P.302.

93. Freygang W.H., Trautwein W. The structurial implications of the linear electrical properties of cardiac Purkinje strands. // J. Gen. Physiol.1970. V.55.-P. 524-547.

94. Garrey W.E. // Physiol. Rev. 1924. - V.4. - P. 215.

95. Garrey W.E. //Amer. J. Physiol. V. 33. - P. 397.

96. Harary J., Farley B. In vitro studies on single beating rat heart cells. Growth and organization. // Exptl. Cell Res. 1962. - V. 29, N 3. - P. 454-474.

97. Heathcote H.L. // Physical Chem. 1902. - V. 37. - P. 368.

98. Heathcote H.L. // J. Soc. Chem. I nol. 1907. - V. 26. - P. 899.

99. Hodghin A.L., A.F. Huxley. // J. Physiol. 1952. - V. 117.1. P.500.

100. Kamigama A. Matsuda F. Electrophysiological properties of the canine ventricular fiber. // Jap. J. Physiol. 1966. V. 16. - P. 407.

101. Kawamura K. Electron microscopic studies on the cardiac conduction system of the dog. // Jap. Circulation J. 1961. - V. 25. - P. 594597.

102. Kawamura K. Electron microscope study of the cardiac conduction system. // Jap. Circulation J. V. 30. - P. 149-151.

103. Kosek J., Marek M. // Phys. Rev. Lett. 1995. - V. 74. - P. 2134.

104. Legato M.J. Ultrastructure of the atrial, ventricular and Purkinje cell with special reference to the genesis to arythmia. // Circulation. 1973. - V. 7, N 1. - P. 178-189.

105. Lillie R.S. // Biol. Rev. 1936. - V. 16. - P. 216.

106. Lillie R.S. // J. Gen. Physiol. 1920. - V. 3. - P. 107.

107. Мое G.K., Rheinboedt W.C., Abiedskov T.A. // Amer. Heart. J. -1964. V. 67. P. 200.

108. Me Nutt N.S., Weinstein R.S.The ulrtastructure of the nexus. A corellated thin-section and freese-aleavy study. // J. Cell. Biol. 1970. - V. 47, N12. - P. 666-688.

109. Мое G.K., Mendes C. Modern concepts of cardiovascular disease (Basis of pharmacotherapy of card arch.), 31,739, 1962.

110. Мое G.K. // Arch. Int. Pharmacotym. 1962. V. 140. P. 183.

111. Мое G.K., Abiedskov T.A. // Amer. Heart. J. 1959. V. 54. P.59.

112. J.W. Moore. // Proc. IEEE. 1968. - V. 56. P. 895.

113. Muir A.R. An electron microscope study of the embryology of the intercalated disc in the heart of the rabbit. // J Biophys. Biochem. Cytol. -1957. V. 2.-P. 193-202.

114. D. Noble, R.W. Tsien J. // Physiol. 1969. V. 200. - P. 233.

115. D. Noble. // Physiol Rev. 46, 1, 1966.

116. D. Noble. //Physiol. 1962. V. 160. - P. 317.

117. Page E. Cat Heart muscle in vitro. The extracellular space. // J. Gen. Physiol. 1962. - V. 46. P. - 201-213.

118. Page E. Solomon A.K. Cat heart muscle in vitro. Cell volumes and intracellular concentrations in papillary muscle. // J. Gen. Physiol. 1960. -V.44. - P. 327-344.

119. Pappano A.I., Sperelakis N., Low К conductance and low resting potentials of isolated single cultured heart cells // Amer.J.Physiol.-1969-V.217, N4. P. 1076.

120. Sacomoto Y. Goto M., 1970, A study of the membrane constants in the dog myocardium. // Jap. J. Physiol. V. 20. - P. 30-41.

121. Sacomoto Y. Membrane characteristics of the canine papillary muscle fiber. // J. Gen. Physiol. V. 54. P. 765-781.

122. G.H. Sharp, R.W. Joynes. // Biophys. J. 1980. - V. 31, P. 403-424.

123. Sperelakis N. Propagation mechanisms in heart. // Annual Rev. Physiol. 1979. - V.41. - P. 441-457.

124. Szekeres L., Papp G.T. Experimental cardiac arhythmias and antiarrhythmic drugs. Academia Kiado, Budapest, 1971.

125. Taylor R.E. in Physical techniques in biological research, VG Academic Press, N.Y., 1963.

126. Tomita T. Membrane capacity and resistance of mammalian smooth muscle. // J. Theor. Biol. V. 12. P. 216-227.

127. Van der Pol. // A eta Med. Scand., Suppl. 1939. - V. 108. - P. 76.

128. Watanabe Y„ Dreifus L.S. // Amer. Heart J. 1968. - V. 76.1. P. 114.

129. West T.G., Landa J.F. // Amer. J. Physiol. 1962. - V. 202.1. P.232.

130. Weidmann S., 1966, The diffusion of radiopotassium across intercalated discs of mamalian cardiac muscle. // J. Physiol. V.187. P.323-342.

131. Weidmann S. Electrical constants of trabecular muscle from mammalian heart. // J. Physiol. 1970. - V. 210. - P. 1041-1054.

132. Weidmann S. Heart: Electrophysiology // Ann.Rev.Physiol. 1974. - V.36.-P.155-169.

133. Wiggers С.Т., Wegria R. // Amer. J. Physiol. 1939. - V. 128.1. P.500.

134. Wiggers C.T. // Amer. Heart J. 1940. - V. 20. - P. 399.

135. Woodbary J.H., Grill W.E. On the problem of impulse conduction in the atrium. // Su: Nervous inhibition. 1961. - P. 124.

136. Woodwary J.M. Cellnear electrophysiology of the heart. // Handbook of Physiology. Circulation. 1962. - V.l. - P. 237-286.

137. Yamagiwa K. Length-to-length transmission as the basis of conduction, and some related considerations about the lows of excitation. // Jap. J. Physiol. 1955. - V. 5. - P. 131-159.

138. Yamagiwa K. Active and resizing membrane resistance and safety factor co-operating for conducting. // Jap. J. Physiol. 1956. - V. 6. - P. 167174.

139. Yamagiwa K. Again on the local response as a small area activity. // Jap. J. Physiol. 1960. - V. 10. - P. 456-470.российсилО rsctycae*^annotc^i i- £ О!