2025-11-23T07:58:16.677890

Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects

Amster, Rivera, Arredondo
The main purpose of this paper is to analyze the dynamics of the system of time-delay differential equations (DDEs) \begin{equation*} \begin{split} \dot{T}(t)&=T(t) f(t,T(t))-γE(t)T(t),\\ \dot{E}(t)&=σ+ \frac{pE(t)T(t-τ_1)}{g+a T(t-τ_1)}-\frac{mE(t)T(t-τ_2)}{g+a T(t-τ_2)}-ηE(t), \end{split} \end{equation*} where $T=T(t)$ and $E=E(t)$ represent the concentrations of tumor and effector cells at the time $t$. The coefficients $σ$, $μ$, $γ$, and $η$ are all positive, and $f(t, T)$ represents the relative growth rate of tumor cells, corresponding to a generalized logistic growth function that describes periodic time chemotherapeutic effects. The parameter $τ_1 \in \mathbb{R}_{\ge 0}$ is the response time delay of the immune system (mediated by effector cells) to an invasion of tumor cells, while $τ_2 \in \mathbb{R}_{\ge 0}$ represents the time delay of tumor cells in response to the appearance of effector cells.
academic

Периодические решения в системе опухоль-иммунитет с временной задержкой и эффектами химиотерапии

Основная информация

  • ID статьи: 2510.11135
  • Название: Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects
  • Авторы: Pablo Amster, Andrés Rivera, John A. Arredondo
  • Классификация: math.DS (Динамические системы)
  • Дата публикации: 13 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2510.11135

Аннотация

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

\dot{T}(t) &= T(t)f(t,T(t)) - \gamma E(t)T(t) \\ \dot{E}(t) &= \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t) \end{align}$$ где $T(t)$ и $E(t)$ обозначают концентрации опухолевых клеток и эффекторных клеток в момент времени $t$ соответственно; $f(t,T)$ представляет относительную скорость роста опухолевых клеток, соответствующую обобщённой логистической функции роста, описывающей периодические эффекты химиотерапии; $\tau_1$ — временная задержка ответа иммунной системы на инвазию опухолевых клеток, $\tau_2$ — временная задержка ответа опухолевых клеток на появление эффекторных клеток. ## Исследовательский контекст и мотивация ### Предпосылки проблемы 1. **Взаимодействие опухоль-иммунная система**: Опухоль представляет собой аномальный рост тканей организма; в нормальных условиях иммунная система может её контролировать, однако образование опухоли часто связано с дефектами иммунной системы 2. **Значимость эффектов временной задержки**: Явления временной задержки широко распространены в биологических системах; распознавание опухоли иммунной системой и её ответ требуют времени, равно как и ответ опухоли на иммунную атаку 3. **Периодические эффекты химиотерапии**: Химиотерапия обычно применяется по периодической схеме введения, что требует учёта её временных эффектов на рост опухоли ### Исследовательская мотивация 1. **Ограничения существующих моделей**: Классические модели опухоль-иммунитет обычно используют логистический рост, однако экспериментальные данные показывают, что модель Гомпертца обеспечивает лучшее соответствие 2. **Влияние параметров задержки**: Необходимо глубокое исследование влияния различных параметров задержки на динамику системы 3. **Существование периодических решений**: Исследование условий существования периодических решений системы при периодическом воздействии химиотерапии ## Основные вклады 1. **Построена обобщённая модель опухоль-иммунитет с временной задержкой типа Ричардса**, содержащая два различных параметра задержки $\tau_1$ и $\tau_2$ 2. **Проведён полный анализ существования и устойчивости точек равновесия**, предоставлена классификация точек равновесия при различных условиях на параметры 3. **Определены критические значения задержки для переключения устойчивости**, предоставлены явные интервалы устойчивости 4. **Доказано существование периодических решений**, установлены ветви периодических решений от автономной системы к неавтономной системе методом продолжения 5. **Даны условия бифуркации Хопфа**, проанализирован переход системы от устойчивого равновесия к периодическим решениям ## Детальное описание методов ### Определение задачи Исследование динамических свойств системы опухоль-иммунитет с двойной временной задержкой, включая: - Анализ существования и устойчивости точек равновесия - Влияние параметров задержки на устойчивость системы - Доказательство существования периодических решений - Анализ бифуркации Хопфа ### Архитектура модели #### Основная структура модели Система (2) описывает взаимодействие опухолевых клеток T(t) и эффекторных клеток E(t): **Уравнение опухолевых клеток**: $$\dot{T}(t) = T(t)f(t,T(t)) - \gamma E(t)T(t)$$ - Первый член: естественный рост опухолевых клеток - Второй член: уничтожение опухолевых клеток эффекторными клетками **Уравнение эффекторных клеток**: $$\dot{E}(t) = \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t)$$ - $\sigma$: постоянный приток эффекторных клеток - Второй член: активация эффекторных клеток опухолевыми клетками (функциональный ответ типа Холлинга II) - Третий член: истощение эффекторных клеток - Четвёртый член: естественная смертность эффекторных клеток #### Конструкция функции роста Используется функция роста типа Ричардса: $$f(t,T) = r\left[\frac{(1-T^\beta)}{\beta} - b(t)\right]$$ где: - $\beta \in (0,1]$: параметр формы - $b(t)$: периодическая функция химиотерапии - При $\beta=1$ вырождается в логистический рост - При $\beta \to 0$ стремится к росту Гомпертца ### Технические инновации 1. **Двойной механизм задержки**: Различие между временной задержкой активации иммунитета $\tau_1$ и задержкой истощения иммунитета $\tau_2$ более точно отражает биологические процессы 2. **Обобщённая модель роста**: Параметр $\beta$ объединяет логистический и гомпертцевский режимы роста 3. **Техника масштабирования параметров**: Преобразование (14) упрощает анализ параметров 4. **Кривые переключения устойчивости**: Предоставлены границы устойчивости на плоскости $(\tau_1,\tau_2)$ ## Теоретический анализ ### Анализ точек равновесия #### Точка равновесия без опухоли Для всех значений параметров система всегда имеет точку равновесия без опухоли $(0,\sigma/\eta)$. **Условия устойчивости** (Теорема 10): Определим $\Delta := \gamma\sigma - rb\eta$, тогда: - Если $\Delta < 0$: точка равновесия без опухоли неустойчива - Если $\Delta = 0$: точка равновесия без опухоли устойчива - Если $\Delta > 0$: точка равновесия без опухоли асимптотически устойчива #### Точка сосуществования При $a=0$ точка сосуществования $(T^*,E^*)$ удовлетворяет уравнению (17): $$\mu(T^{*\beta} - b)T^* + T^{*\beta} = b - \sigma$$ **Условия существования** (Предложение 8): - Условие †: $\sigma < b$ и $\mu \leq 0$ - Условие ‡: $\sigma \leq b$ и $0 < \mu \leq \mu_c$ где $\mu_c := \left(\frac{1}{b}\left(\frac{1-\beta}{1+\beta}\right)^{\beta-1}\right)^{1/\beta}$ ### Анализ устойчивости #### Случай малых задержек **Теорема 12**: Для достаточно малых $a,\tau_1,\tau_2 \geq 0$: 1. При условии † точка $(T^*,E^*)$ локально асимптотически устойчива 2. При условии ‡, если выполнены дополнительные условия, то $(T^*,E^*)$ устойчива #### Случай общих задержек **Теорема 14**: Когда $\tau_1 = \tau_2 = \tau$, существует критическое значение $\tau_c$ такое, что: - $\tau < \tau_c$: точка равновесия устойчива - $\tau > \tau_c$: точка равновесия теряет устойчивость и происходит бифуркация Хопфа Критическое значение определяется уравнением: $$G(y) = y^4 + \left[(\lambda_1^*)^2 + (\lambda_2^*)^2\right]y^2 + (\lambda_1^*\lambda_2^*)^2 - N^2 = 0$$ ### Существование периодических решений **Теорема 16** (Теорема продолжения): Пусть $(T^*,E^*)$ — допустимая точка равновесия при $b \equiv b_0$. Тогда существуют константы $\phi^*,\tau_j^* > 0$ и $r^* > 0$ такие, что при $\|\phi\|_\infty < \phi^*$, $\tau_1 < \tau_1^*$, $\tau_2 < \tau_2^*$ система (2) имеет по крайней мере одно $\omega$-периодическое решение $(T,E)$, удовлетворяющее $\|(T-T^*,E-E^*)\|_\infty < r^*$. **Схема доказательства**: 1. Проверка отсутствия нетривиальных периодических решений линеаризованной системы (Лемма 4) 2. Применение теоремы об имплицитной функции и альтернативы Фредгольма 3. Построение гладкой ветви периодических решений от автономной системы без задержки к неавтономной системе с задержкой ## Численная верификация ### Установка параметров В работе используются типичные значения параметров из литературы: - $g \approx 10^6$ (относительно большое) - После масштабирования $a \to a/g = O(10^{-6})$ (относительно малое) ### Кривые переключения устойчивости **Теорема 15** даёт кривые переключения устойчивости на плоскости $(\tau_1,\tau_2)$: $$C_{s,k} = \{(\tau_1^{\pm}(y), \tau_2^{\pm}(y)) \in \mathbb{R}_+^2 : y \in I, s,k \in \mathbb{Z}\}$$ где $\tau_1^{\pm}$ и $\tau_2^{\pm}$ определяются уравнениями (35) и (36). ## Связанные работы ### Модели роста опухоли 1. **Классические модели**: Экспоненциальный рост, логистический рост, рост Гомпертца 2. **Экспериментальная верификация**: Множество исследований показывают, что модель Гомпертца лучше всего соответствует экспериментальным данным 3. **Корреляция параметров**: Два параметра модели Гомпертца демонстрируют значительную корреляцию ### Модели взаимодействия опухоль-иммунитет 1. **Kuznetsov и др. (1994)**: Построена базовая структура модели опухоль-иммунитет 2. **Li и др. (2021, 2022)**: Исследованы эффекты одиночной и двойной задержки 3. **Вклад данной работы**: Расширение на обобщённый рост типа Ричардса и периодические эффекты химиотерапии ### Теория дифференциальных уравнений с задержкой 1. **Анализ устойчивости**: Метод характеристического уравнения, метод функций Ляпунова 2. **Теория бифуркаций**: Бифуркация Хопфа, переключение устойчивости 3. **Периодические решения**: Метод продолжения, теория неподвижных точек ## Заключение и обсуждение ### Основные выводы 1. **Полная классификация точек равновесия**: Даны условия существования и классификация всех возможных точек равновесия 2. **Явные условия устойчивости**: Предоставлены явные критерии устойчивости и критические значения задержки 3. **Существование периодических решений**: Доказано существование периодических решений при периодическом воздействии химиотерапии 4. **Явления бифуркации**: Раскрыт механизм перехода системы от устойчивого состояния к колебаниям ### Ограничения 1. **Упрощение параметров**: Предположение $a=0$, хотя имеет биологическое обоснование, ограничивает общность 2. **Недостаточная численная верификация**: Теоретические результаты требуют более полного численного моделирования 3. **Отсутствие биологической верификации**: Предсказания модели требуют поддержки большего объёма экспериментальных данных ### Направления будущих исследований 1. **Трёхмерное расширение**: Рассмотрение дополнительных типов клеток (например, регуляторных T-клеток) 2. **Пространственная неоднородность**: Введение членов диффузии для исследования пространственных паттернов 3. **Стохастические эффекты**: Учёт влияния шума на динамику системы 4. **Оптимальное управление**: Разработка оптимальных схем химиотерапии на основе модели ## Глубокая оценка ### Преимущества 1. **Теоретическая строгость**: Полный математический анализ, ясное изложение доказательств 2. **Разумность модели**: Чёткий биологический контекст, ясное значение параметров 3. **Методологическая инновация**: Комбинация двойной задержки и обобщённой функции роста 4. **Практическая применимость**: Теоретическое обоснование для разработки схем химиотерапии ### Недостатки 1. **Слабая численная часть**: Раздел 6 "Численная верификация" отмечен как "в разработке" 2. **Отсутствие оценки параметров**: Не предоставлены методы оценки реальных параметров 3. **Недостаточная верификация модели**: Отсутствует сравнение с клиническими данными ### Влияние 1. **Теоретический вклад**: Обогащение приложений дифференциальных уравнений с задержкой в биоматематике 2. **Прикладная ценность**: Новые подходы к математическому моделированию лечения опухолей 3. **Методологическое значение**: Применение метода продолжения для анализа периодических решений имеет общий характер ### Области применения 1. **Теоретические исследования**: Дифференциальные уравнения с задержкой, теория динамических систем 2. **Биоматематика**: Модели роста опухолей, моделирование иммунной системы 3. **Медицинские приложения**: Оптимизация схем химиотерапии, прогнозирование эффективности лечения ## Библиография Статья цитирует 27 связанных работ, включая: - Классические работы по моделям роста опухолей - Теорию устойчивости дифференциальных уравнений с задержкой - Математические модели взаимодействия опухоль-иммунитет - Теоремы о бифуркациях и существовании периодических решений --- **Общая оценка**: Это высокотеоретическая работа по математической биологии, вносящая значительный вклад в область дифференциальных уравнений с задержкой и математическое моделирование системы опухоль-иммунитет. Теоретический анализ детален, методы инновационны, однако численная верификация требует дальнейшего развития. Работа имеет важное значение для понимания сложной динамики системы опухоль-иммунитет.