Аннотация результатов, полученных в 2018 году
Экспериментальные методы глобального скрининга и визуализации процессов в органах и тканях, действующих на различных пространственных и временных масштабах, используемые в современной иммунологии, привели к формированию беспрецедентной по степени разрешения и детализации системы представлений о структурной организации и регуляторных сетях имеющих место в функционирующей иммунной системе. Для понимания системной динамики процессов и рационального управления компонентами иммунной системы, необходима разработка эффективной методологии математического описания иерархически организованных, пространственно структурированных и нелинейно регулируемых физических, биохимических и физиологических процессов на основе многомасштабных и гибридных моделей с высоким уровнем детализации описания. Развитие методологии многомасштабного и гибридного моделирования иммунных процессов осуществляется применительно к исследованию инфекции вирусами иммунодефицита человека (ВИЧ-1) и экспериментальной инфекции мышей вирусами лимфоцитарного хориоменингита (ВЛХМ). В первый год реализации проекта нами был проведен критический обзор современных исследований в области гибридного моделирования биомедицинских процессов. Рассмотрены дискретно-непрерывные модели с нелинейными зависимостями скоростей процессов от концентрации участвующих компонент и многоуровневым описанием регуляции на внутриклеточном-, клеточно-популяционном-, тканевом- и организменном масштабах. Проанализированы существующие подходы к интегративному математическому моделированию организации и структуры иммунных процессов в лимфатической системе, в частности, лимфатических узлах (ЛУ), на основе непрерывно-детерминированного, дискретно-стохастического и гибридного описаний. Разработана детерминированная математическая модель развития ВИЧ инфекции с детальным описанием динамики иммунного ответа и процесса заражения клеток. Модель сформулирована в виде системы интегро-дифференциальных уравнений с запаздывающим аргументом. Стохастическая версия модели построена на основе ветвящегося процесса с взаимодействием частиц и учитывает стадии созревания клеток, вирионов и стадию размножения лимфоцитов-киллеров. Доказана глобальная разрешимость задачи Коши и положительная инвариантность системы. Получены условия асимптотической устойчивости тривиального положения равновесия, соответствующего отсутствию инфекции в организме человека. Численно исследованы возможности искоренения ВИЧ инфекции в течение фиксированного периода времени. Получена оценка вероятности попадания переменных стохастической модели в состояние, соответствующее тривиальному положению равновесия детерминированной модели. Нами разработано семейство вложенных математических моделей описывающих популяционную динамику вирусов и клеток иммунной системы на основе систем из двух нелинейных ОДУ, дифференциальных уравнений с запаздывающим аргументом, реакционно-диффузионных уравнений в одномерном пространственном рассмотрении и реакционно-диффузионных уравнений с запаздывающим временным аргументом. Для бистабильных систем дифференциальных уравнений с запаздыванием, моделирующих развитие вирусных инфекций и противовирусного иммунного ответа, реализован эффективный подход к численному построению оптимальных возмущений стационарных состояний с высокой вирусной нагрузкой, переводящих систему в состояние с низкой вирусной нагрузкой. Показана предпочтительность построения оптимальных возмущений с использованием в качестве локальных норм - нормы Соболева для оценки амплификации возмущения с целью разработки адекватных режимов терапии. Для систем уравнений нелинейных реакционно-диффузионных моделей с запаздыванием с заданными начальными и граничными условиями были разработаны эффективные конечно-разностные численные методы, использующие неявную по времени аппроксимация нестационарных членов и диффузионных слагаемых, трехточечное представление системы разностных уравнений, метод прогонки и два различных алгоритма аппроксимации запаздывающих переменных. Разработана методология качественного анализа двухкомпонентных реакционно-диффузионных систем с запаздыванием путем редукции их к одному уравнению того же класса. С помощью численного моделирования показано существование квазиволновых режимов пространственно-временной динамики, не обладающих регулярной периодической структурой и демонстрирующих сложные апериодические осцилляции. Разработана, откалибрована и валидирована математическая модель движения клеток иммунной системы в лимфатическом узле в двумерной пространственной постановке. Модель сформулирована виде системы стохастических дифференциальных уравнений движения в рамках ньютоновской механики для клеток, с параметрическим описанием сил межклеточного взаимодействия, диссипативных сил вязкости и сил, стохастически задающих активную подвижность клеток. Оценены характеристики потенциальных функций, определяющих взаимодействие и подвижность клеток. Проведен сравнительный анализа методов численного интегрирования уравнений модели с помощью симплектических и несимплектических разностных схем. Выполнена численная реализация гибридной модели пространственно-временной динамики противовирусного иммунного ответа в лимфатическом узле в трехмерной постановке. В модели используется случайно-дискретное описание популяций клеток, детерминированное описание пространственных полей вирусов и гуморальных факторов в Т-клеточной области ЛУ, наличие сети фибробластных ретикулярных клеток и изменение пространственной структуры ЛУ в процессе развития фиброза в ходе ВИЧ инфекции. Для эффективной реализации модели на графических процессорах, с использованием архитектуры параллельных вычислений CUDA, была переработана её алгоритмическая структура, что позволило реализовать параллельное выполнение ряда операций и минимизировать использование динамически выделяемой памяти. Разработана математическая модель внутриклеточной регуляции активности генов связанных с активацией синтеза интерферона при ВИЧ инфекции. Модель сформулирована в виде системы нелинейных дифференциальных уравнений с запаздывающими аргументом. Проведена калибровка модели для её интеграции в гибридную модель противовирусного иммунного ответа в лимфатическом узле. С помощью метода цифровой квантификации численности клеток по экспериментальным данным транскриптомного анализа селезенки у мышей построены кинетические кривые в фазовом пространстве системы из 125 переменных, описывающие численности клеточных популяций для острой и хронической инфекции ВЛХМ, с целью их последующего использования в задаче многомасштабного моделирования. Были выявлены субпопуляции моноцитов, макрофагов и CD8+ Т-клеток, согласованно меняющихся с экспрессией генов, контролирующих развитие инфекции ВЛХМ по хроническому или острому варианту. Нами построена уточненная математическая модель внутриклеточной репликации ВИЧ с большой размерностью пространства состояний (24 компоненты), характеризующих биохимические реакции всего цикла репликации от момента заражения клетки до секреции инфекционных вирусов. Модель описывает стадии репликации ВИЧ, на которые действуют ингибиторы обратной транскриптазы, интегразы, протеазы. Был выполнен начальный цикл исследований по оптимизации доз лекарственных препаратов при их комбинированном применении с целью уменьшения суммарной продукции вирусов клеткой при ограничении на суммарную дозу комбинации из трех лекарств. Исследована возможность расширения модели путем описания процессов метаболизма ингибитора обратной транскриптазы (AZT) с целью учета клеточной токсичности данного препарата. С помощью стохастической версии модели на основе алгоритма Гиллеспи были оценены характеристики трубки решений, внутри которой находится 95% случайных реализаций стохастической модели. Проведено исследование по построению редуцированной версии модели внутриклеточной репликации ВИЧ. Стадии репликации, избыточные с точки зрения конкретной задачи моделирования, параметризовывались введением нелинейных зависимостей с запаздыванием. Разработана математическая модель динамики ВИЧ инфекции, описывающая каскадный процесс клонального развития Т-клеточного иммунного ответа и комбинированной терапии ВИЧ-инфекции (сочетание антиретровирусной терапии и блокады рецептора PD-L1) с учетом динамики регуляторных Т-лимфоцитов. В модели описывается, что ингибиторы обратной транскрипции блокируют способность вирусных частиц эффективно инфицировать клетки, а ингибиторы протеазы приводят к образованию неинфекционных вирусных частиц. Путем применения информационно-теоретических критериев оценивания сложности моделей для всего множества возможных наборов параметров, определяющих фенотип развития инфекции и эффект блокировки рецептора, выбраны наиболее правдоподобные гипотезы (сочетания параметров) о механизмах, опосредующих эффект применения блокады рецептора PD-L1. С помощью модели предсказан эффект применения блокады рецептора для пяти ВИЧ-инфицированных пациентов для двух механизмов реализации эффекта, связанных с изменением (1) определённых кинетических параметров или (2) доли клеток, отвечающих на стимуляцию. Для решения задачи моделирования процессов переноса в лимфатической систем были выполнены исследования по разработке модели и численных методов расчета конвективного течения и диффузии лимфы для трехмерной модели лимфатического узла. Развиваемый подход связан с применением теории потенциалов к решению задачи Неймана для уравнения Лапласа для сферической области. Применение теории потенциалов позволяет перейти от уравнений в частных производных в трехмерной области описывающих течение жидкости в пористой среде, к интегральным уравнениям на границе области. Разработана математическая модель для решения задачи стационарного течения лимфы (вязкое течение) в пористой среде, описываемого законом Дарси-Бринкмана для упрощенной геометрической модели лимфатического узла в виде двух вложенных сфер. Пространство между сферами соответствует субкапсулярному синусу лимфатического узла, а внутренняя сфера – Т-клеточной области ЛУ. Данные области имеют различные характеристики гидравлической проводимости. В составной расчетной области было численно промоделировано вязкое течение и фильтрация лимфы. Разработан программный комплекс на языке FORTRAN 90 для решения поставленной задачи методом граничных конечных элементов.

Аннотация результатов, полученных в 2019 году
Современный этап математического моделирования иммунных процессов определяется необходимостью развития широкого спектра математических моделей различной природы, глубиной детализации переменных состояния системы, рассмотрением нескольких уровней и масштабов регуляции и описанием пространственно-временной динамики иммунных реакций. В отчетном году нами были разработаны, численно реализованы и исследованы математические модели клеточного, клеточно-популяционного, тканевого и системного уровней представления иммунных процессов. Был выполнен систематический обзор существующих стохастических и детерминистических моделей внутриклеточных процессов регуляции деления, дифференцировки, функционального подавления, гибели и пространственной подвижности иммунных клеток под действием цитокинов, хемокинов и вирусных компонент. На основе реакционно-диффузионных уравнений с запаздыванием по времени (РДУЗ) была разработана математическая модель вирусной инфекции в двумерной пространственной постановке. Исследованы условия формирования различных режимов пространственно-временной динамики в зависимости от геометрии области и характера начальных условий, например, переход от решений с радиальной симметрией в нерегулярные осциллирующие решения. Построена одномерная (в пространстве генетических признаков) модель вирусной инфекции на основе реакционно-диффузионного интегро-дифференциального уравнения, описывающая генетическую структуру квазивидов. Изучено влияние нелокальных взаимодействий вирусных геномов между собой при инфицировании клеток-мишеней, с кросс-реактивными иммунными клетками и генотип-зависимой резистентностью к антивирусной терапии на диверсификацию и генетическую изменчивость вирусов. Была разработана модель для исследования генетической изменчивости вирусов под действием иммунного ответа в двумерном пространстве признаков, т.е. рассмотрение генотипического и физического пространств. Установлено, что нелокальная зависимость иммунного ответа от пространственной плотности вирусной популяции приводит к качественному изменению режима динамики инфекции. Разработана многомасштабная модель пространственно-временной динамики распространения ВИЧ-инфекции в лимфоидных тканях, интегрирующая 1) двумерную модель движения клеток в рамках ньютоновской механики изотропных частиц, 2) стохастическую модель инфицирования клеток свободными вирионами и путем передачи вирусных геномов от клетки к клетке, 3) детерминистическую модель внутриклеточной репликации ВИЧ в продуктивно инфицированных клетках, 4) стохастическую модель последовательной активации противовирусного ответа системы молекул интерферона I типа, 5) двумерную реакционно-диффузионную модель распространения вирионов ВИЧ и молекул интерферона I типа в лимфоидных органах, 6) стохастическую модель гибели клеток и поступления новых клеток в лимфоидный орган. Модель позволяет предсказать кинетику суммарного числа вирионов ВИЧ в лимфоидном органе после поступления в него продуктивно инфицированной антигенпрезентирующей клетки при разных иммунопатофизиологических условиях. Выполнен переход к трехмерной постановке в задаче моделирования движения иммунных клеток в лимфоидных органах на основе второго закона Ньютона. Ранее разработанная двумерная модель программно реализована для случая трехмерного пространства. Разработан алгоритм, с использованием технологий OpenMP и CUDA, агентного моделирования клеточной динамики в лимфатическом узле в трехмерной пространственной постановке с учетом сети фибробластных ретикулярных клеток (ФРК) для описания миграции Т-лимфоцитов, межклеточных взаимодействий, пролиферации клеток, хемотаксиса и перераспределения веществ в межклеточном пространстве с течением времени. Исследовано влияние целостности сети ФРК на кинетику миграции иммунных клеток в процессе развития иммунного ответа. Были найдены с помощью алгоритмов, разрабатываемых в группе Ю.М. Нечепуренко, области значений параметров, для которых имеет место свойство мульти-стабильности (показано сосуществование до 4-х положений равновесия) для модели экспериментальной инфекции вирусами лимфоцитарного хориоменингита (ВЛХМ), определены условия их устойчивости и режимы переходов между ними на основе оптимальных возмущений. Для математической модели противовирусного иммунного ответа Марчука-Петрова исследована задача гарантированного вычисления всех стационарных состояний. В контексте задачи малопараметрическго моделирования внутриклеточной репликации ВИЧ для использования в гибридных моделях были рассмотрены существующие подходы к уменьшению размерности моделей и получены «оптимальные» по сложности и точности модели. Для моделирования процессов иммунной регуляции определяющих развитие острого или хронического вариантов инфекции ВЛХМ выявлены различия в динамике ответа системы интерферона для данных фенотипов инфекции. Были экспериментально (группа проф. А. Майерханса) получены данные численностей субпопуляций дендритных клеток (XCR1-cDC) и антиген-индуцированных CD8+ Т-клеток, находящихся в состоянии истощения в селезенке у мышей в динамике. Были подготовлены для моделирования процессов контроля хронической инфекции, выявленных на примере экспериментальной инфекции ВЛХМ, клинические данные для когорт ВИЧ-инфицированных пациентов с различным фенотипом инфекции, включая уровни вирусной нагрузки, концентрацию CD4+ T клеток и относительные численности субпопуляций дендритных клеток (pDC, cCD1, cCD2) среди клеток экспрессирующих маркер CD45. Разработана уточненная модель хронической фазы динамики ВИЧ инфекции, с помощью которой, впервые для ВИЧ инфекции получены практически значимые предсказания терапевтического эффекта применения препарата (анти-PD-L1), восстанавливающего функциональную активность как Т-, так и В-лимфоцитов. Показано, как количественный эффект зависит от фенотипа хронической фазы инфекции пациентов с различными характеристиками вирусной нагрузки и иммунного статуса. Был разработан гибридный подход (совместно с др. И. Сазоновым, Университет г. Суонси, Великобритания) к моделированию стохастической динамики вирусных инфекций на основе марковских процессов при учете запаздывания процесса продукции вирусов и при наличии переменных разного масштаба (численности популяций), позволяющий использовать стохастический и детерминистический режимы вычислений компонент траекторий в зависимости от численностей переменных модели вирусной инфекции. Построено и исследовано три варианта детерминированной модели, описывающей динамику ВИЧ инфекции в начальный период времени после инфицирования здорового человека. Модели различаются по набору переменных, степени усложнения схемы возникновения инфицированных клеток разных типов и вирусных частиц, а также схемы формирования клеточного иммунного ответа. Получены условия асимптотической устойчивости и неустойчивости тривиального положения равновесия, соответствующего отсутствию инфекции в организме человека. Построена двухкомпартментная модель ВИЧ инфекции, которая содержит уравнения для переменных, отражающих численности компонентов инфекционного процесса в лимфатических узлах (ЛУ) и находящихся в процессе переноса между лимфоузлами. Развитие ВИЧ инфекции в отдельно взятом ЛУ описывается высокоразмерной системой дифференциальных и интегральных уравнений, притоки клеток и вирусных частиц в ЛУ задаются с помощью запаздывающих переменных. Численности клеток и вирусных частиц, находящихся в стадии переноса между лимфоузлами, описываются интегральными уравнениями типа свертки или эквивалентными дифференциальными уравнениями с запаздыванием специального вида. Изучены условия на параметры модели, допускающие возможность искоренения ВИЧ инфекции за некоторый начальный промежуток времени после инфицирования здорового человека за счет эффективной работы специфического клеточного иммунного ответа. Разработана стохастическая версия многокомпартментной модели ВИЧ инфекции. Модель допускает описание продолжительности переходов клеток и вирусных частиц между компартментами в форме достаточно произвольных функций, удовлетворяющих ряду естественных ограничений. Для численного решения разработаны подход на основе метода Монте-Карло, с применением алгоритмов, сохраняющих предысторию развития случайного процесса, включая алгоритмы работы со списками, отражающими точечные распределения, которые являются стохастическими аналогами интегральных уравнений типа свертки. Показано, что вероятность искоренения инфекции после инфицирования здорового человека зависит от начального количества вирусных частиц, параметров, отражающих работу клеточного иммунного ответа в каждом лимфоузле и определяющих перенос клеток и вирусных частиц между лимфоузлами. Для решения задачи трехмерной фильтрации лимфы в лимфоузле была построена, обоснована и численно реализована модель на основе интегральных представлений для поля скоростей и давлений (совместная работа с д.ф.-м.н. А.В. Сетухой). Доказана эквивалентность предложенной интегральной и исходной дифференциальной постановок (закон Дарси-Бринкмана), в том числе, для решения более общей задачи фильтрации вязкой жидкости в ограниченной области с конечным числом источников (стоков). Доказано, что поля скоростей и давления могут быть представлены в виде соответствующих поверхностных потенциалов, удовлетворяющих построенной системе граничных интегральных уравнений. Для задачи сопряжения течения в неоднородных областях предложено соответствующее интегральное представление. Выполнена квантификация 3-х мерных морфометрических параметров и топологических характеристик сети кондуитов лимфатических узлов мышей (на базе Института иммунобиологии госпиталя Ст. Галлен, Швейцария, проф. Л. Буркхард). Структурные параметры включают гистограммы распределений длин сосудов, степеней вершин и углов в точках ветвления сосудов. Полученные данные позволяют перейти к разработке физиологически содержательной многомасштабной модели ЛУ описывающей количественный баланс жидких фаз в лимфоидных органах. Разработана многомасштабная ВИЧ инфекции, интегрирующая динамику иммунных процессов на двух физиологически различных уровнях рассмотрения – в лимфоидных тканях и крови – с учетом стохастической кинетики оборачиваемости лимфоцитов. Результаты исследований по проекту были опубликованы в нескольких профильных журналах Q1 и отражены в пресс-релизах центральных информационных агенств: https://nauka.tass.ru/nauka/7109351; https://ria.ru/20191113/1560890750.html https://www.gazeta.ru/science/news/2019/11/13/n_13692140.shtml; https://www.kommersant.ru/doc/4157546; https://indicator.ru/mathematics/matematicheskaya-model-rasschitala-chto-12-tysyach-podvizhnyh-limfocitov-pobedyat-vich-19-06-2019.htm https://ria.ru/20190619/1555691649.html http://www.ras.ru/news/shownews.aspx?id=3d3b667e-5a2a-4eaf-ac01-0a2390fbe0a3&print=1 https://www.gazeta.ru/science/news/2019/06/19/n_13110013.shtmlhttps://indicator.ru/news/2019/06/19/matematicheskaya-model-rasschitala-chto-12-tysyach-podvizhnyh-limfocitov-pobedyat-vich/

Аннотация результатов, полученных в 2020 году
Исследования по проекту в 2020 году были направлены на изучение широкого спектра вопросов – от анализа данных по динамике вирусных инфекций до построения, численной реализации и исследования различных классов математических моделей в рамках детерминистического, стохастического и гибридного подходов, в том числе – многомасштабных пространственно распределённых описаний иммунных процессов. В результате проведенных исследований были (1) сформулированы новые положения о возможных механизмах полной элиминации ВИЧ-1 инфекции из организма человека, (2) идентифицированы перспективные мишени для антиретровирусной терапии, (3) показано, что реакция системы интерферона I типа на ранних стадиях цитопатических вирусных инфекций играет ключевую роль в предотвращения тяжелых форм течения заболеваний. Разработан вычислительный подход к моделированию и анализу лимфатической системы человека, основанный на теории графов. В его основе лежат два разных метода построения графов лимфатической системы: на основе анатомических данных и на основе набора физиологических правил сетевой организации системы. Преобразование простого графа построенного по анатомическим данным в ориентированный граф реализовано путем количественной оценки направления потоков лимфы в сети лимфатических сосудов с использованием уравнения Пуазейля и закона сохранения потоков для точек ветвления сосудов. Получены фундаментальные топологические характеристики для графовых моделей лимфатической системы человека. Для проведения исследований по расчетам пространственных полей сигнальных молекулярных факторов в лимфатическом узле был разработан алгоритм построения трехмерной полномасштабной геометрической модели структурной организации лимфатического узла. Алгоритм позволяет осуществить интеграцию сети кровеносных сосудов с компонентами лимфатического узла - сетью фибробластных ретикулярных клеток, В-клеточными фолликулами, субкапсулярным синусом. В основе алгоритма лежим биофизический механизм роста кровеносного русла, обусловленный гипоксией, феноменологически представленной градиентом уровня кислорода. Разработана упрощенная модель дренажной функции лимфатического узла, характеризующаяся хорошим согласием свойств решений с экспериментальными данными. При прохождении через лимфоузел часть лимфы всасывается в кровь. Построена модель всасывания лимфы в кровеносную систему в лимфатическом узле. Предполагается, что кровеносные сосуды равномерно распределены в области всасывания. Сток лимфы в кровь описывается уравнением Старлинга для неизвестных скорости и давления лимфы и известного среднего давления в кровеносных сосудах лимфатического узла. В целом, данная проблема решена как задача фильтрации жидкости в однородной пористой среде с гладкой границей, на которой заданы неоднородные граничные условия. Построен алгоритм моделирования движения клеток в трехмерной пространственной постановке на основе второго закона Ньютона. Впервые для моделирования активной пространственной миграции клеток в лимфатическом узле применены методы пространственной статистики, в частности - распределение Мизеса-Фишера на сфере в 3D. Сила, определяющая активную подвижность клеток, определяется как реализация коррелированных случайных блужданий на основе Марковского процесса с дискретным временем. При этом использование распределений Мизеса-Фишера и Гауссова для определения направления и скачка величины силы подвижности обеспечивает обратную корреляцию угла поворота траектории движения клетки и амплитуды силы активной подвижности. Реализован прототип среды моделирования иммунных процессов с использованием программных средств разработки Qt Creator и языка Qt. В частности, разработан программный пакет Composite Mathematical Immunology, предоставляющий пользователю возможность решения следующих проблемно-ориентированных задач: развёртывать интерактивную модель графов лимфатической системы произвольной топологии, представимых в виде матриц связности и выполнять с ними различные действия; проводить моделирование клеточной динамики в двумерной постановке с возможностью гибко модифицировать свойства элементов расчётной области; выполнять детерминированный и стохастический расчёт модели внутриклеточной репликации ВИЧ-1 с возможностью модифицировать параметры модели в рамках физиологически допустимых значений. Для многостадийной модели ключевых процессов размножения ВИЧ в зараженной Т-клетке путем анализа чувствительности функционала, характеризующего суммарную продукцию вирусов зараженной клеткой выделены процессы, наиболее сильно влияющие на производство вирусов, и проведено их ранжированные по силе воздействия на суммарный выход вирусных частиц из клетки. Для некоторых из этих процессов пока не существует эффективных лекарственных ингибиторов. Новыми, дополнительными мишенями для антивирусной терапии могут быть (1) этап сборки ВИЧ-1 на внутренней стороне клеточной мембраны, зависящий от уровня белка Gag; (2) блок нелинейной регуляции репликации ВИЧ белками Tat и Rev на транскрипционном и пост-транскрипционном уровнях. Получены зависимости вероятности продуктивного инфицирования клеток от числа связанных вирионов и рассчитаны гистограммы распределений числа интегрированных в хромосомы клетки провирусных ДНК в диапазоне от 1 до 100 заражающих клетку вирусных частиц. На основе разработанной немарковской стохастической модели и ее детерминированного аналога найдены соотношения между параметрами инфекционного процесса, при которых вероятность искоренения ВИЧ-1 инфекции в течение фиксированного периода времени после инфицирования равна единице (случай R0 < 1, где R0 – базовое репродуктивное число). Для случая R0 > 1 построена оценка вероятности искоренения ВИЧ-1 инфекции в течение заданного периода времени после инфицирования в зависимости от начального количества вирусных частиц и параметров, отражающих эффективность иммунного ответа. Оценены значения параметров модели, при которых результаты вычислительного эксперимента воспроизводят известные эмпирические данные по динамике ВИЧ-1 инфекции в организме человека в локальной фазе инфекции. Для указанных параметров модели вероятность искоренения ВИЧ-1 инфекции при инфицировании человека небольшим количеством вирусных частиц в течение первых 8-12 суток после инфицирования равна единице. Результаты вычислений указывают на существенное влияние специфического иммунного ответа на искоренение ВИЧ-1 инфекции в организме при малом начальном количестве вирусных частиц (5-10 вирионов), участвующих в инфицировании отдельного лимфоузла. Изучены закономерности пространственно-временной динамики вирусной инфекции в организме с учетом распределения вирусных квазивидов по генотипу и с учетом распределения инфекции по ткани или органу. Была разработана математическая модель этого процесса с распределенными параметрами, основанная на интегро-дифференциальных уравнениях в частных производных. Показано, что распространение инфекции подобно распространению реакционно-диффузионной волны, и что в частном случае двумерная задача может быть сведена к одномерной задаче при помощи разделения переменных. Рассматривалась одномерная задача в зависимости от собственного значения, как приближение 2D проблемы. В общем случае разработанный переход позволяет приближенно исследовать решение двумерной задачи на основе решения одномерной задачи. Были определены зависимости вирулентности вируса (скорость волны) и уровня вирусной нагрузки (амплитуда волны) от ширины распределения генотипов, силы иммунного ответа, уровня иммунитета и исходной вирусной нагрузки (исходного состояния) и показано, что увеличение ширины распределения генотипов увеличивает вирулентность вируса и вирусную нагрузку. Путем системного междисциплинарного анализа множества разнообразных данных о развитии и структуре популяции ВИЧ в организме человека до и в ходе антиретровирусной терапии сформулирована новая концепция функционального лечения ВИЧ инфекции (элиминация резервуара латентно-инфицированных клеток). Она основана на использовании физиологического механизма вымывания части активированных субпопуляций (в том числе, латентно-инфицированных) CD4+ Т-клеток более зрелого фенотипа за счет поступления менее дифференцированных клеток в иммунную систему и их конкуренции за факторы выживания. Предложено, за счет применения препаратов дополнительной поликлональной активации, интенсифицировать процесс обновления Т-клеточной популяции (т.е., вымывания субпопуляций CD4+ T-клеток памяти, включая латентно-инфицированные CD4+ Т-лимфоциты). Данная концепция лечения ВИЧ-1 инфекции («rinse and replace»), открывает возможность реализации нового мультимодального подхода к элиминации всех латентно-инфицированных клеток, путем применения индукторов активации антиген-презентирующих клеток в сочетании с антиретровирусной терапией. На основе анализа соотношений между кинетикой репликации цитопатических вирусов, местом локализации инфекционного процесса, реакцией врожденной системы иммунитет (синтез интерферона первого типа, ИФН-1) и тяжестью инфекционного заболевания сформулирована гипотеза о том, что защита от цитопатической вирусной инфекции, развивающейся в периферических органах и тканях, обеспечиваемая реакцией ИФН-1 чрезвычайно уязвима по отношению к росту скорости репликации вирусов. Предложено, что ведение ИФН-1 зараженным людям на ранних стадиях цитопатических вирусных инфекций (ВИЧ-1, SARS-CoV-2) необходимо для предотвращения развития тяжелых форм течения периферических заболеваний. Выпущен сборник трудов под редакцией основных участников проекта: Bocharov, G., Ludewig, B., Meyerhans, A., Volpert, V., eds. (2020). Mathematical Modeling of the Immune System in Homeostasis, Infection and Disease. Lausanne: Frontiers Media SA. doi: 10.3389/978-2-88963-461-3 (276 pages). ISSN 1664-8714, ISBN 978-2-88963-461-3 Результаты исследований по проекту были опубликованы в нескольких профильных журналах первого квартиля (Q1) и, некоторые из них, отражены в следующих публикациях СМИ: https://nauka.tass.ru/nauka/8912173 https://www.gazeta.ru/science/news/2020/07/08/n_14644921.shtml https://indicator.ru/medicine/novyi-podkhod-lecheniyu-vich-09-07-2020.htm https://iz.ru/1032939/2020-07-08/rossiiskie-uchenye-nashli-novyi-podkhod-klecheniiu-vich