Цена доставки диссертации от 500 рублей 

Поиск:

Каталог / ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ / Математика

Математическое моделирование противовирусного Т-клеточного имунного ответа

Диссертация

Автор: Бочаров, Геннадий Алексеевич

Заглавие: Математическое моделирование противовирусного Т-клеточного имунного ответа

Справка об оригинале: Бочаров, Геннадий Алексеевич. Математическое моделирование противовирусного Т-клеточного имунного ответа : диссертация ... кандидата физико-математических наук : 01.01.11 Москва, 1984 149 c. : 61 85-1/2812

Физическое описание: 149 стр.

Выходные данные: Москва, 1984






Содержание:

ВВЕДШИЕ
ГЛАВА I» Аналитическое исследование математической модели противовирусного Т-клеточного иммунного ответа
11 Существование, единственность, неотрици-тельность и ограниченность решений модели
12 Исследование стационарного решения, описывающего состояние:-здорового организма
13 Исследование стационарных решений, описывающих хронические формы заболеваний
ГЛАВА 2 Численное моделирование противовирусного
Т-клеточного иммунного ответа
21 Алгоритм численного интегрирования дифференциальных уравнений с запаздывающим аргументом
22 Численное моделирование основных форм динамики вирусных заболеваний
23 Исследование защитно-повреждающего действия Т-киллеров
24 Исследование зависимости клинических характеристик вирусного заболевания от характеристик функциональной активности Т-хелперов
2*5 Численное моделирование реакции Т-системы иммунитета на увеличение количества вирусов в организме при хроническим течении вирусного заболевания
ГЛАВА 3 Численная имитация клинико-лабораторных данных по вирусному гепатиту В с помощью математической модели противовирусного Т-клеточного иммунного ответа Ю
3»1» Анализ клинико-лабораторных данных по трем формам вирусного гепатита В
3,2» Численная имитация клинико-лабораторных данных по вирусному гепатиту В

Введение:
Познание закономерностей функционирования иммунной системы и механизмов выздоровления при вирусных инфекциях является в настоящее время, в связи с заметным ростом числа хронических и гипертоксических форм вирусных заболеваний, одной из фундаментальных проблем современной иммунологии и медицины. Эффективным средством изучения этих закономерностей является математическое моделирование исследуемых процессов и систем, свидетельством чего может служить значительный рост интереса клиницистов и иммунологов к применению в медицинской практике разнообразных математических методов [7, 14, 19-20, 30-32, Зб].
В 1974 году в ВЦ СО АН СССР академиком Г.И.Марчуком в тесном сотрудничестве с академиком АМН СССР Р.В.Петровым и академиком АМН СССР Н.И.Нисевич были инициированы работы по математическому моделированию в иммунологии и медицине. В результате проведенных исследований к 1980 г. были созданы и внедрены в клиническую практику математические методы оценки тяжести и изучения динамики патологических процессов при вирусном гепатите В и пневмониях, сравнения эффективности различных методов лечения [7, 19-20, 30-32, 39]. Была построена и исследована простейшая математическая модель инфекционного заболевания, позволившая в рамках единой схемы объяснить общие закономерности формирования и развития инфекционных заболеваний [1-6, 28-29]. С 1980 г. исследования в этих направлениях продолжаются в Отделе вычислительной математики АН СССР.
К этому времени в клинической иммунологии были получены новые сведения о роли Т- и В-систем иммунитета при вирусных инфекциях, о кооперации клеток иммунной системы (в частности,были созданы теории "двойного сигнала" и "двойного распознавания" [33, 41, 65]), о механизмах развития вирусных заболеваний.(в частности, показано ведущее значение Т-системы иммунитета в освобождении организма от вирусо^. Вместе с тем выяснилось, что способность Т-лимфоцитов-киллеров разрушать зараженные вирусами клетки обуславливает не только защитный эффект, но и является причиной развития поражений различных органов и тканей, т.е. иммунопатологических реакций [47, 48]. Особенно ярко защитное и повреждающее действие Т-киллеров проявляется в случае вирусного гепатита В [9]. Значительная доля хронических и летальных исходов этого заболевания ставит перед врачами задачу выбора эффективных способов управления и коррекции Т-клеточного иммунного ответа. Клинические исследования в этом направлении активно ведутся на кафедре детских инфекционных болезней П МОЖЖ (зав. кафедрой академик АМН Н.И.Нисевич), сотрудники которой участвуют в постановке задачи регуляции иммунного ответа при инфекционных заболеваниях.
Неоднозначность результатов лечения, сложность и многоплановость реакции иммунной системы организма обусловили применение методов математического моделирования как одного из способов исследования процессов, связанных с функционированием Т-системы иммунитета и зараженного вирусами органа при вирусных заболеваниях.
Приведем краткий обзор результатов, полученных в области математического моделирования иммунологических процессов и инфекционных заболеваний к настоящему времени.
В современной медицине под инфекцией понимается совокупность физиологических и патологических процессов, возникающих и развивающихся в организме при внедрении в него патогенных микроорганизмов или вирусов, которые вызывают нарушение постоянства его внутренней среды (прежде всего антигенного36) и физиологических функций [51]. Одним из основных способов защиты организма от живых тел и веществ, несущих на себе признаки генетической чужеродности, является система иммунитета [41]. Суть реакции иммунной системы на появление генетически чужеродного материала (антигена), в том числе и возбудителей болезни, состоит в производстве специфических антител и клеток-киллеров способных нейтрализовать и уничтожить антигены [29].
Создание адекватных математических моделей инфекционных заболеваний связано, прежде всего, с построением и изучением моделей реакций иммунной системы организма на антигенное воздействие.
Как известно, реакции системы иммунитета осуществляемые лимфоцитами делятся на две основные группы: реакции клеточного типа (Т-система) и реакции гуморального типа (В-система). Большинство построенных математических моделей было посвящено моделированию гуморального иммунного ответа. В основу этих моделей положена клонально-селекционная теория Бернета [41].
Уравнения математических моделей строятся на основе соотношений баланса кавдого компонента процесса.
Наиболее полное описание гуморальной иммунной реакции было предложено в 1970 г. Белл ом, построившим математическую модель клональной селекции и антителопродукции [60-62], В свою модель прим, автора во] Белл включил следующие клеточные популяции: иммуноциты (клетки-мишени), пролиферирукхцие клетки, клетки памяти и плазматические клетки. В модели рассматривались также двухвалентные антитела и одновалентный неразмножающийся антиген. В работах [61-62] Белл ввел, в соответствии с экспериментальными данными, гетерогенность аффинитета иммунокомпетентных клеток и антител, исследовал роль поливалентности антигена, эффекты роста аффинитета в зависимости от дозы антигена, а также механизмы толерантности, вызываемой повторными инъекциями низких или высоких доз антигена. В модели Белла предполагалось, что стимуляция иммунокомпетентных клеток имеет место лишь тогда, когда доля их рецепторов связанных с антигеном лежит между двумя эмпирически определенными предельными значениями.
Представления, положенные в основу модели Белла были в дальнейшем развиты Бруни и соавт. в работах [66, 70-73]. Оставаясь в рамках теории клональной селекции, они ввели непрерывное распределение иммунокомпетентных клеток и антител по аффинитету К в некотором интервале (Кр 1^). Их модель представляет собой систему из пяти интегродифференциальных уравнений с распределенным параметром К и описывает динамику численности иммунокомпетентных клеток (В-лимфоцитов), плазмоклеток, антител, неразмножащегося антигена и иммунных комплексов. В модели предполагалось, что генерация иммунокомпетентных клеток стволовыми клетками происходит со скоростью, зависящей от аффинитета клеточных рецепторов. Скорость введения антигена и начальное распределение аффинитета являлись варьируемыми параметрами* Моделирование первичного гуморального ответа при умеренных дозах антигена дало хорошее согласие с экспериментальными данными.
Чтобы избавиться от неопределенности в числовых значениях ряда коэффициентов модели, была проведена идентификация параметров модели [66] на основе специально разработанной экспериментальной программы. Авторами было показано, что модель относится к классу систем с ограниченной производной [70], Исследуя поведение модели в пространстве состояний, были доказаны некоторые важные структурные свойства системы уравнений модели, в частности, глобальное существование и единственность решений для физически возможного диапазона входных воздействий и начальных условий, а также глобальная асимптотическая устойчивость единственного положения равновесия внутри физически допустимого множества состояний,
Молер и соавт. построили и исследовали сходную с моделью Бруни математическую модель гуморальной иммунной реакции, в которой предполагалось, что все антитела обладают одинаковым аффинитетом, а также модели с дискретным набором значений аффинитета [38]. Базовая модель представляет собой систему из пяти нелинейных дифференциальных уравнений описывающих динамику концентраций иммунокомпетентных клеток (В-лимфоцитов)# плазматических клеток, концентраций свободных активных центров антител, иммунных комплексов, свободного антигена. Уравнения состояния относятся либо к суммарной концентрации компонентов для данного среднего аффинитета К, либо к концентрации, соответствующей одному из многих значений К. Популяция свободного антигена является входным воздействием, на которое реагирует иммунная система. Молер и соавт. [38, 89], используя свою математическую модель для имитации опыта по введению мыши антигена в полном адъюванте Фрейнда, получили результаты согласующиеся с экспериментальными данными. В ходе математического анализа модели были доказаны теоремы о существовании и единственности решений, положительной инвариантности уравнений модели, а также выпуклость множества достижимых состояний для упрощенной системы уравнений модели. Показано, что модель представляет собой совокупность билинейных систем, поэтому для ее идентификации можно использовать функции Уолша. Это иллюстрируется на примере идентификации некоторых параметров модели. Обсуждаются вопросы управления моделью. Выделены две переменные аддитивного управления (скорость производства иммуноцитов из стволовых клеток и скорость введения антигена) и три - мультипликативного, являющиеся нелинейными функциями состояния системы. Молер считал, что скорость поступления иммунокомпетентных клеток из костного мозга является важнейшим фактором в иммунотерапии. Им было отмечено, что если в рассматриваемую модель гуморальной иммунной реакции ввести дополнительную компоненту учитывающую ухудшение состояния иммунной системы или органа, модель может быть использована для анализа эффективности такой терапии.
Была также предложена модификация базовой модели, которая позволила исследовать пространственную неоднородность иммунной системы, т.е. описать течение иммунного процесса в различных органах иммунной системы [88]. Эта модификация представляет собой совокупность базовых моделей, аппроксимирующих динамику процесса в крови, селезенке, лимфе и лимфатических узлах. Был проведен анализ чувствительности модели к изменениям ее параметров. Экспериментальная проверка модели дала хорошие результаты. Мо-лером отмечено, что непротиворечивые данные для идентификации ключевых параметров этой модели отсутствуют.
Заметим, что рассмотренные модели гуморальной иммуной реакции не описывали процессы кооперации Т- и В-систем иммунитета, которые в действительности имеют место при индукции гуморального иммунитета на большинство тимусзависимых антигенов.
В работе Молера [87, 90] рассмотрены модели Т- и В-системы иммунитета, описывающие кооперацию между ними при антителогене-зе. При моделировании иммунной реакции было получено близкое соответствие с некоторыми экспериментальными данными по первичной гуморальной реакции.
Представленные выше математические модели гуморальной иммунной реакции обладают одной общей особенностью (это подчеркивают и их авторы): модели дают качественное и количественное соответствие с экспериментальными данными при использовании умеренных доз антигена. Например, начальный вариант модели Белла [601 оказался неудовлетворительным при описании иммунного ответа на большие дозы антигена, так как привел к биологически невозможным количествам образованных в ходе реакции клеток. Это является следствием того предположения, что клоны иммунокомпе-тентных клеток реагирующие на конкретный антиген не зависят от остальной части иммунной системы и единственный регуляторный механизм - количество антигена. Открытие антигенных детерминант (идиотипов) у антител и рецепторов иммунокомпетентных клеток послужило основанием для создания Ерне сетевой теории иммунитета [82]. Согласно этой теории, когда антиген инициирует иммунную реакцию, то вновь созданные иммуноглобулины (циркулирующие антитела и рецепторы лимфоцитов) через свои идиотипы стимулируют вторую»антиидиотипическую реакцию и т.д. Наряду с такими прямыми стимулирующими воздействиями в теории иммунной сети постулируются и супрессорные обратные связи: каждое поколение антител и клеток стремится подавить предыдущее. Таким образом, предполагается существование в иммунной системе механизма регуляции действующего на идиотипы посредством активных центров других молекул антител или лимфоцитарных рецепторов. Идеи Ерне были развиты в дальнейшем в моделях Хоффмана [80-81]. Модели иммунной сети, представляя несомненный теоретический интерес, пока не являются конструктивной основой моделирования инфекционных заболеваний (в отличие от моделей основанных на селекции клонов антигеном).
Свойством любой инфекции является способность возбудителя заболевания к размножению. Поэтому первым шагом в переходе от математических моделей иммунных реакций на неразмножающийся антиген, ориентированных на конкретный эксперимент, к математическому моделированию инфекционных заболеваний является введение в уравнения моделей членов описывающих размножение антигена. К настоящему времени получен ряд результатов по моделированию иммунной реакции в случае самовоспроизводящегося антигена, т.е. реакции против бактерий, вирусов и опухолевых клеток. Общим для этих моделей является простота описания иммунной системы.
В работе [63] Белл предлагает для описания иммунного ответа на активные размножающиеся антигены простейшую из возможных моделей - систему из двух нелинейных- обыкновенных дифференциальных уравнений, описывающих динамику размножения антигена и антител. Позже в модель было добавлено уравнение динамики плазмокле-ток. На этих моделях была изучена зависимость решений от параметров, в частности, определены условия существования устойчивых стационарных состояний, предельных циклов, расходящихся колебаний [96, 97]. Решениям таких типов могут соответствовать состояние здорового организма, те или иные формы течения "болезней": хронические, затяжные, летальный исход, выздоровление. Хотя полной элиминации антигена в модели нет, под выздоровлением понимается падение концентрации антигена до некоторого минимума.
В работах Диброва Б.Ф. и соавт. [15-18] предложена математическая модель гуморальной иммунной реакции. Модель описывает динамику В-лимфоцитов, антител и антигена и учитывает запаздывание в появлении антител и клеток памяти. Показано, что ход иммунной реакции, т.е. амплитуда колебаний решений существенным образом зависит от продолжительности запаздывания: так, при малом запаздывании возможно стационарное существование антигена и антител, при его увеличении возможны циклическое течение процесса или неограниченное размножение антигена. При малом запаздывании велика вероятность рецидива в ходе заболевания, а его незначительное увеличение снижает указанную вероятность. Делается вывод, что изменение запаздывания может служить эффективным методом управления иммунными процессами. Определены области параметров, соответствующие различным режимам поведения: асимптотическому снижению количества антигена, стремлению его количества к постоянному значению, периодическому изменению количества антигена и неограниченному росту. В работе [18] рассматриваются стохастические аспекты моделирования иммунной реакции, связанные с вопросом о полном уничтожении антигена. Эта проблема вызвана тем, что в рамках детерминированной модели -в случае, если эффективная скорость размножения антигена в условиях здорового организма положительна, уничтожение антигена (даже асимптотическое) невозможно. В работе проведена оценка вероятности полного уничтожения антигена в тех случаях, когда согласно детерминистическим уравнениям модели количество антигена становится малым.
Для исследования закономерностей развития инфекционных заболеваний Г.И.Марчук в 1975 г. построил простейшую математическую модель инфекционного заболевания. В модели рассматриваются следующие процессы, из которых, в принципе, слагается взаимодействие организма с чужеродным размножающимся антигеном:
1) Размножение чужеродных клеток, обобщено названных антигеном или вирусом и нейтрализация их антителами. Понятие антитела также обобщено и под ним имеются в виду как иммуноглобулины, так и рецепторы Т-лимфоцитов.
2) Реакция иммунной системы организма на антигенное вторжение, характеризующаяся накоплением иммунокомпетентных антитело-образующих клеток, которые формируются через время Т . Стимулирующим сигналом считается комплекс антиген-антитело. Учитывается естественная гибель плазмоклеток и их генерация из стволовых клеток.
3) Производство антител, их расход на связывание антигена и естественный распад.
4) Вирус, проникая в организм, поражает какую-либо систему или орган. Поражение может достигать такой степени, при которой нарушается нормальная работа органа, что приводит к ослаблению функционирования иммунной системы. Для характеристики степени поражения органа-мишени вводится обобщенное понятие и масса пораженного органа". Она зависит от патогенности вирусов (различной для разных заболеваний) и скорости восстановления поражения. Снижение эффективности.иммунного ответа при поражении учитывается с помощью непрерывной невозрастащей функции (т) :
Введение в математическую модель уравнения для массы пораженного органа дает основания считать простейшую модель иммунологической реакции организма простейшей моделью инфекционного заболевания» В целом модель представляет собой систему из 4 нелинейных дифференциальных уравнений с запаздывающим аргументом. Переменные модели:- концентрация возбудителей — , концентрация антител - Р (-Ь) , концентрация плазмоклеток - С (•?) характеристика массы пораженного органа - т (-Ь).
Было выполнено математическое исследование модели [1-6, 28-29]. Для системы уравнений модели были доказаны теоремы о существовании, единственности и неотрицательности решений. Показано, что модель имеет два типа стационарных решений, проведено исследование их асимптотической устойчивости. Решения первого типа характеризуются нулевым стационарным количеством вирусов
0 и интерпретируются как состояние "здорового организма". Стационарные решения другого типа характеризующиеся тем, что сопэ-ЬЪ-О, интерпретируются как "хроническая" или латентная" форма болезни. Введено понятие иммунологического барьера организма. Численные эксперименты по моделированию заражения организма малой дозой вирусов показали, что существует четыре качественно отличающиеся друг от друга типа решений, которые мож но интерпретировать как формы протекания инфекционного заболевания: субклиническая, острая с выздоровлением, хроническая, летальный исход. Был математически исследован метод лече-зия хронических форм "обострением" и построена модель биостиму-щции, в которой один из антигенов является возбудителем хронической инфекции, а другой (биостимулятор) вводится в организм специально для нарушения устойчивости хронического процесса за счет перераспределения ресурсов иммунной системы. В ходе изучения влияния температурной реакции организма на динамику иммуного ответа показано, что повышение температуры снижает степень тяжести заболевания и уменьшает максимальное значение концентрации антигена, подавление температурной реакции может привести к переходу острой формы заболевания в затяжную или хроническую [I]. Получен ряд биологических следствий:
I* Возникновение той или иной формы болезни не зависит от дозы заражения, а определяется иммунологическим статусом организма по отношению к данному типу вирусов (набором параметров модели).
2. Возбудители хронических, инфекций обладают вялой динамикой.
3. Хронические формы болезни обусловлены недостаточно эффективной стимуляцией иммунной системы.
4. Возможно лечение хронических форм через повышение концентрации вирусов в организме, т.е. обострение болезни.
В рамках простейшей модели инфекционного заболевания был рассмотрен также вопрос лечения препаратом типа антибиотик ?3].
Таким образом, простейшая математическая модель инфекционного заболевания позволила выяснить ряд закономерностей развития реакций иммунной системы и протекания болезни и указать возможные пути их лечения, составив базу для дальнейших исследований в области математического моделирования инфекционных заболеваний.
Вместе с тем, простейшая математическая модель инфекционного заболевания не описывала ряда важных процессов, происходящих в организме при вирусном заболевании. Следующим шагом в развитии идей, положенных в основу простейшей модели инфекционного заболевания, стало построение в 1981 г. Г.И.Марчуком и Р.В.Петровым математической модели противовирусного иммунного ответа [33].
Перечислим основные положения незавершенной еще теории противовирусного иммунитета, которые не отражены в простейшей математической модели инфекционного заболевания, и которые обусловили необходимость создания математической модели противовирусного иммунного ответа [41, 471:
• иммунный ответ направлен не только против свободного вируса циркулирующего в крови или лимфе, но и против зараженных клеток хозяина;
• процесс выздоровления определяется скоростью и полнотой удаления из организма зараженных клеток;
• способность цитотоксических лимфоцитов разрушать инфицированные клетки своего организма не только обусловливает защитный эффект, но и является причиной развития поражений органов и тканей;
• процессы узнавания, включения эффекторных лимфоцитов в пролиферацию и цитотоксическое действие накопившихся Т-киллеров запускаются не самим вирусным антигеном, а комплексом: вирусный антиген - трансплантационный антиген главной системы гистосовме-стимости (МНС) хозяина;
• макрофаг является клеткой перерабатывающей антиген и формирующий на своей поверхности комплексные детерминанты из вирусных антигенов и антигенов МНС;
• для стимуляции В-лимфоцитов и предшественников Т-киллеров требуется двойной сигнал: от стимулированных вирусами макрофагов и Т-лимфоцитов-помощников.
Изложенные положения легли в основу математической модели противовирусного иммунного ответа» В настоящее время эта модель является предметом исследований, которые позволят, в конечном итоге, приступить на ее основе к моделированию конкретных вирусных заболеваний и выбору эффективных методов лечения [10, 68, 75],
В заключение обзора имеющихся результатов по математическому моделированию иммунологических процессов и инфекционных заболеваний отметим, что результаты полученные в других направлениях исследований в данной области моделирования, таких как применение теории управления и оптимизации, анализ специфических для теории систем свойств моделей, фильтрация и теория оценивания [Юб], не представлены главным образом потому, что содержание диссертационной работы не связано с этими вопросами. Отметим также, что к настоящему времени имеется несколько обзоров, посвященных различным аспектам математического моделирования иммунологических процессов [13, 23, 25, 34, 38]. Более полное и конкретное представление о математическом моделировании в иммунологии и медицине можно получить из монографии Г-И.Мар-чука "Математические модели в иммунологии" [29] и книги "Theoretical immunology" под редакцией Дж.Белла [102]. Следует выделить также сборник работ "Regulation of Immune Response Dynamics" под редакцией Ч.Делиси, и Ж.Эрно [100], посвященный описанию различных механизмов регуляции в иммунной системе и анализу их с помощью математических моделей.
В приведенном выше обзоре показано, что построение математической модели противовирусного иммунного ответа было обусловлено отсутствием адекватного описания в рамках имевшихся математических моделей иммунологических реакций и инфекционных заболеваний, отдельных существенных процессов противовирусного иммунитета. Отметим при этом, что в самое последнее время появились работы Белла и Перельсона [94, 95] посвященные изучению механизмов Т-клеточной цитотоксичности на уровне отдельной клетки.
Рассмотрим более подробно математическую модель противовирусного Т-клеточного иммунного ответа. Процессы протекающие в организме при попадании в него возбудителей вирусных инфекций связаны, превде всего, с развитием иммунного ответа на вирусные антигены и поверхностные антигены зараженных вирусами клеток, В математической модели противовирусного иммунного ответа рассматриваются два типа иммунного ответа: гуморальный с выработкой антител В-системой лимфоцитов и Т-клеточный с накоплением цитотоксических лимфоцитов киллеров. Антитела нейтрализуют свобода ©циркулирующие вирусные частицы, а Т-киллеры обнаруживают зараженные вирусами клетки и уничтожают их, причем общепринятым в настоящее время в теоретической иммунологии является тезис о ведущей эффекторной и регуляторной роли в противовирусном иммунитете популяции Т-лимфоцитов [33, 41, 48].
Математическая модель противовирусного Т-клеточного иммунного ответа была сформулирована Г.И.Марчуком и Р.В.Петровым в работе [33] как частный случай полной модели и в дальнейшем модифицирована в соответствии с их работой [35]. Уравнения модели описывают процессы заражения вирусами клеток чувствительной к вирусу ткани, воспроизводства последними полноценных вирусов, цитотоксическое действие вирусов в отношении зараженных клеток и реакцию Т-системы иммунитета на вирусные антигены. Соответственно в качестве переменных модели рассматриваются: Vj. {?) - количество свободноциркулирующих вирусов; CvW - доля зараженных клеток органа-мишени; m(i) - доля пораженных клеток органа-мишени; Му ("Ь) - количество стимулированных вирусами макрофагов; H(t) - количество Т-лимфоцитов-хелперов клеточного иммунитета; Е (t) - количество Т-лимфоцитов-киллеров. При построении уравнений модели были сделаны следующие предположения:
1. Количество макрофагов М в организме считается постоянным и достаточно большим для того, чтобы прирост стимулированных макрофагов М v был пропорционален количеству свободных вирусов
2. Все клетки одного клона, являющиеся результатом деления Т-хелперов или Т-киллеров, идентичны.
3. Для стимуляции Т-лимфоцита-помощника Н требуется один сигнал от Mv; для стимуляции Т-лимфоцита-киллера - двойной сигнал от Mv и Н.
4. Лимфоциты-помощники Н после осуществления акта взаимодействия с Т-киллерами свой жизненный цикл заканчивают, т.е. принята гипотеза "моногамии" Т-хелперов,
5. Т-лимфоциты-хелперы и Т-киллеры генерируются с некоторой постоянной скоростью в кроветворных органах.
6. В течение определенного времени зараженные клетки Cv выполняют свои нормальные функции. Их гибель обусловлена либо развитием необратимого вирусного поражения, либо цитотоксическим действием Т-киллеров.
7. Вирусы освободившиеся из разрушенных клеток, а также секретируемые в кровь и лимфу зараженными клетками, способны заражать имеющиеся здоровые клетки, число которых равно С (1-Сv- m), где С - обозначает общее число клеток органа мишени. Кроме того, свободноциркулирующие вирусы связываются антителами, уровень кос* торых V предполагается неизменным в течение всего процесса вирусного заболевания, поскольку в модели не рассматривается динамика гуморального иммунного ответа.
Уравнения модели являются следствием балансовых соотношении по отношению к основным компонентам участвующим в рассмотренных процессах в предположении, что "организм" описывается однородным замкнутым объемом, в котором все компоненты процесса равномерно перемешаны.
Система уравнений математической модели противовирусного Т-клеточного иммунного ответа имеет следующий вид:
S = p&A*4cv?-W'-^M-iciv, сИ-Ь
ЕЕ) iL=e oit Р j>esHMvlt-Te)H(i-te)E(t-g-M.HEj ,ile.v. gv=?C(l-Cv-m)VrgeCvE-6mCv H oit в v m v 0 m
В дальнейшем,для удобства переобозначим р8т на ?С на на , а величину Р +у„М + кб) через X.
Процесс заражения здорового организма в момент ?=0 может быть описан с помощью следующих начальных данных: оьуД мдо)=о, н(о)-н*
Е(0)«н*, СуЮ) = 0, т(0) = 0; (Б2)
МУ(5)ВД=0} зеГ-^,0); НД5)НС5)Е(8)=0, зеС-г^о).
В модели с помощью функции учитывается снижение эффективности функционирования Т-системы в процессе развития вирусного заболевания обусловленное поражением органа-мишени, причем принято, что:
КУ))= 1-т , 0< П04 1.
Учитывается запаздывание при появлении в результате деления клона новых клеток ( и "Се - величины запаздывания).
Целью диссертационной работы является исследование математической модели противовирусного Т-клеточного иммунного ответа и использование ее для анализа реальных клинико-лабораторных данных на примере вирусного гепатита В. В задачи исследования модели входят:
1) аналитическое исследование модели;
2) численное исследование модели с целью анализа:
• защитного и повреждающего действия Т-киллеров при вирусных заболеваниях,
• взаимосвязей между формами развития вирусных заболеваний и характеристиками Т-системы иммунитета, описываемыми в модели;
3) имитация и анализ с помощью модели клинико-лабораторных данных по вирусному гепатиту В.
Для решения указанных задач использовались методы теории дифференциальных уравнений с отклоняющим ся аргументом, теории возмущений, численные методы решения обыкновенных дифференциальных уравнений и интерполяции функций, имитационное моделирование и методы теории чувствительности,.
Настоящая диссертация является первой работой, в которой проведено исследование математической модели противовирусного Т-клеточного иммунного ответа и показана возможность ее применения для анализа реальных данных по конкретному заболеванию -вирусному гепатиту В.
Работа состоит из введения, трех глав и заключения.
В первой главе доказаны теоремы о существовании, единственности и неотрицательности решений модели и ограниченности отдельных компонент вектор-функции решения при неотрицательных начальных условиях; получены достаточные условия асимптотической устойчивости стационарного решения описывающего состояние "здорового организма"; при некоторых ограничениях проведено исследование устойчивости стационарных решений, описывающих хронические формы вирусных заболеваний.
Во второй главе проведено численное моделирование противовирусного Т-клеточного иммунного ответа. Рассмотрен и исследован алгоритм численного решения дифференциальных уравнений с запаздывающим аргументом. Показана возможность описания в рамках модели известных форм развития вирусных заболеваний; изучена связь за-цитного и повреждающего действия Т-киллеров с величиной коэффициентов модели характеризующих скорость стимуляции Т-киллеров и свойства вирусов; определена зависимость некоторых модельных характеристик вирусного заболевания (V^min7mmax,Еmax и др.) от коэффициентов модели, определяющих реактивные свойства Т-сис-темы; показано, что условием возникновения хронических форм заболеваний может являться как "вялая" динамика вирусов в организме, так и слабая реакция Т-системы иммунитета; показано, что эффективность применения метода лечения хронического заболевания обострением существенно зависит от иммунного статуса организма.
В третьей главе проведены численная имитация и анализ кли-нико-лабораторных данных по острой и хроническим формам гепатита В с помощью математической модели противовирусного Т-клеточ-ного иммунного ответа, обсуждаются функции чувствительности для решений модели имитирующих данные по хроническим формам гепатита В.
В заключении сформулированы основные результаты работы.
Перечислены названия институтов, семинаров и конференций, на которых обсуждались результаты работы. Приведен список работ, опубликованных по теме диссертации.