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

Periodische Lösungen in einem Tumor-Immun-Konkurrenzsystem mit Zeitverzögerung und Chemotherapieeffekten

Grundinformationen

  • Paper-ID: 2510.11135
  • Titel: Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects
  • Autoren: Pablo Amster, Andrés Rivera, John A. Arredondo
  • Klassifizierung: math.DS (Dynamische Systeme)
  • Veröffentlichungsdatum: 13. Oktober 2025
  • Paper-Link: https://arxiv.org/abs/2510.11135

Zusammenfassung

Dieser Artikel analysiert hauptsächlich die dynamischen Eigenschaften eines Differentialgleichungssystems mit Zeitverzögerung:

\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}$$ Dabei bezeichnen $T(t)$ und $E(t)$ die Konzentrationen von Tumorzellen bzw. Effektorzellen zum Zeitpunkt $t$; $f(t,T)$ stellt die relative Wachstumsrate der Tumorzellen dar, entsprechend einer verallgemeinerten logistischen Wachstumsfunktion, die periodische Chemotherapieeffekte beschreibt; $\tau_1$ ist die Reaktionsverzögerung des Immunsystems auf die Tumorinvasion, $\tau_2$ ist die Reaktionsverzögerung der Tumorzellen auf das Auftreten von Effektorzellen. ## Forschungshintergrund und Motivation ### Problemhintergrund 1. **Tumor-Immunsystem-Wechselwirkung**: Tumoren sind abnormales Gewebewachstum des Körpers; normalerweise kann das Immunsystem diese kontrollieren, aber die Tumorbildung ist oft mit Immundefekten verbunden 2. **Bedeutung von Zeitverzögerungseffekten**: Zeitverzögerungsphänomene sind in biologischen Systemen weit verbreitet; die Erkennung und Reaktion des Immunsystems auf Tumoren benötigt Zeit, und die Reaktion von Tumorzellen auf Immunangriffe ist ebenfalls verzögert 3. **Periodische Chemotherapieeffekte**: Chemotherapie wird typischerweise nach periodischen Dosierungsschemata verabreicht und erfordert die Berücksichtigung ihrer zeitvariablen Auswirkungen auf das Tumorwachstum ### Forschungsmotivation 1. **Einschränkungen bestehender Modelle**: Klassische Tumor-Immun-Modelle verwenden häufig logistisches Wachstum, aber experimentelle Daten zeigen, dass das Gompertz-Modell eine bessere Anpassung bietet 2. **Auswirkungen von Verzögerungsparametern**: Eine tiefere Analyse der Auswirkungen verschiedener Verzögerungsparameter auf die Systemdynamik ist erforderlich 3. **Existenz periodischer Lösungen**: Erforschung der Existenzbedingungen für periodische Lösungen des Systems unter periodischer Chemotherapiewirkung ## Kernbeiträge 1. **Etablierung eines verallgemeinerten Tumor-Immun-Verzögerungsmodells vom Richards-Typ** mit zwei unterschiedlichen Verzögerungsparametern $\tau_1$ und $\tau_2$ 2. **Vollständige Analyse der Existenz und Stabilität von Gleichgewichtspunkten** mit Klassifizierung unter verschiedenen Parameterbedingungen 3. **Bestimmung kritischer Verzögerungswerte für Stabilitätswechsel** mit expliziten Stabilitätsintervallen 4. **Beweis der Existenz periodischer Lösungen** durch Fortsetzungsmethoden zur Etablierung von Lösungszweigen vom autonomen zum nicht-autonomen System 5. **Angabe von Bedingungen für Hopf-Bifurkation** mit Analyse des Übergangs vom stabilen Gleichgewichtspunkt zu periodischen Lösungen ## Methodische Erläuterung ### Aufgabendefinition Untersuchung der dynamischen Eigenschaften eines Tumor-Immun-Konkurrenzsystems mit doppelter Zeitverzögerung, einschließlich: - Existenz- und Stabilitätsanalyse von Gleichgewichtspunkten - Auswirkungen von Verzögerungsparametern auf die Systemstabilität - Beweis der Existenz periodischer Lösungen - Hopf-Bifurkationsanalyse ### Modellarchitektur #### Grundlegendes Modellgerüst System (2) beschreibt die gegenseitige Wechselwirkung zwischen Tumorzellen T(t) und Effektorzellen E(t): **Tumorzellgleichung**: $$\dot{T}(t) = T(t)f(t,T(t)) - \gamma E(t)T(t)$$ - Erster Term: natürliches Wachstum von Tumorzellen - Zweiter Term: Tötung von Tumorzellen durch Effektorzellen **Effektorzellgleichung**: $$\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$: konstante Zufuhr von Effektorzellen - Zweiter Term: Aktivierung von Effektorzellen durch Tumorzellen (Holling-II-Funktionalantwort) - Dritter Term: Erschöpfung von Effektorzellen - Vierter Term: natürlicher Tod von Effektorzellen #### Wachstumsfunktionsdesign Verwendung einer Richards-Wachstumsfunktion: $$f(t,T) = r\left[\frac{(1-T^\beta)}{\beta} - b(t)\right]$$ Dabei gilt: - $\beta \in (0,1]$: Formparameter - $b(t)$: periodische Chemotherapiefunktion - Wenn $\beta=1$, degeneriert es zu logistischem Wachstum - Wenn $\beta \to 0$, nähert es sich Gompertz-Wachstum an ### Technische Innovationen 1. **Doppelter Verzögerungsmechanismus**: Unterscheidung zwischen Immunaktivierungsverzögerung $\tau_1$ und Immunerschöpfungsverzögerung $\tau_2$ für genauere Abbildung biologischer Prozesse 2. **Verallgemeinertes Wachstumsmodell**: Vereinheitlichung von logistischem und Gompertz-Wachstum durch Parameter $\beta$ 3. **Parameterskalierungstechnik**: Vereinfachung der Parameteranalyse durch Transformation (14) 4. **Stabilitätswechselkurven**: Angabe von Stabilitätsgrenzen in der $(\tau_1,\tau_2)$-Ebene ## Theoretische Analyse ### Gleichgewichtspunktanalyse #### Tumorfreier Gleichgewichtspunkt Für alle Parameterwerte existiert immer ein tumorfreier Gleichgewichtspunkt $(0,\sigma/\eta)$. **Stabilitätsbedingungen** (Theorem 10): Definiere $\Delta := \gamma\sigma - rb\eta$, dann gilt: - Wenn $\Delta < 0$: tumorfreier Gleichgewichtspunkt ist instabil - Wenn $\Delta = 0$: tumorfreier Gleichgewichtspunkt ist stabil - Wenn $\Delta > 0$: tumorfreier Gleichgewichtspunkt ist asymptotisch stabil #### Koexistenzgleichgewichtspunkt Wenn $a=0$, erfüllt der Koexistenzgleichgewichtspunkt $(T^*,E^*)$ die Gleichung (17): $$\mu(T^{*\beta} - b)T^* + T^{*\beta} = b - \sigma$$ **Existenzbedingungen** (Proposition 8): - Bedingung †: $\sigma < b$ und $\mu \leq 0$ - Bedingung ‡: $\sigma \leq b$ und $0 < \mu \leq \mu_c$ Dabei ist $\mu_c := \left(\frac{1}{b}\left(\frac{1-\beta}{1+\beta}\right)^{\beta-1}\right)^{1/\beta}$ ### Stabilitätsanalyse #### Fall kleiner Verzögerungen **Theorem 12**: Für hinreichend kleine $a,\tau_1,\tau_2 \geq 0$: 1. Unter Bedingung † ist $(T^*,E^*)$ lokal asymptotisch stabil 2. Unter Bedingung ‡ ist $(T^*,E^*)$ stabil, wenn zusätzliche Bedingungen erfüllt sind #### Allgemeiner Verzögerungsfall **Theorem 14**: Wenn $\tau_1 = \tau_2 = \tau$, existiert ein kritischer Wert $\tau_c$ derart, dass: - $\tau < \tau_c$: Gleichgewichtspunkt ist stabil - $\tau > \tau_c$: Gleichgewichtspunkt verliert Stabilität und erfährt Hopf-Bifurkation Der kritische Wert wird durch die folgende Gleichung bestimmt: $$G(y) = y^4 + \left[(\lambda_1^*)^2 + (\lambda_2^*)^2\right]y^2 + (\lambda_1^*\lambda_2^*)^2 - N^2 = 0$$ ### Existenz periodischer Lösungen **Theorem 16** (Fortsetzungssatz): Sei $(T^*,E^*)$ ein zulässiger Gleichgewichtspunkt für $b \equiv b_0$, dann existieren Konstanten $\phi^*,\tau_j^* > 0$ und $r^* > 0$, so dass wenn $\|\phi\|_\infty < \phi^*$, $\tau_1 < \tau_1^*$, $\tau_2 < \tau_2^*$, das System (2) mindestens eine $\omega$-periodische Lösung $(T,E)$ mit $\|(T-T^*,E-E^*)\|_\infty < r^*$ besitzt. **Beweisansatz**: 1. Verifikation, dass das linearisierte System keine nichttrivialen periodischen Lösungen hat (Lemma 4) 2. Anwendung des impliziten Funktionensatzes und des Fredholm-Alternativsatzes 3. Etablierung eines glatten Lösungszweigs vom autonomen verzögerungsfreien zum nicht-autonomen Verzögerungssystem ## Numerische Verifikation ### Parametereinstellung Der Artikel verwendet typische Parameterwerte aus der Literatur: - $g \approx 10^6$ (relativ groß) - Nach Skalierung $a \to a/g = O(10^{-6})$ (relativ klein) ### Stabilitätswechselkurven **Theorem 15** gibt die Stabilitätswechselkurven in der $(\tau_1,\tau_2)$-Ebene an: $$C_{s,k} = \{(\tau_1^{\pm}(y), \tau_2^{\pm}(y)) \in \mathbb{R}_+^2 : y \in I, s,k \in \mathbb{Z}\}$$ Dabei werden $\tau_1^{\pm}$ und $\tau_2^{\pm}$ durch die Gleichungen (35) und (36) gegeben. ## Verwandte Arbeiten ### Tumorwachstumsmodelle 1. **Klassische Modelle**: exponentielles Wachstum, logistisches Wachstum, Gompertz-Wachstum 2. **Experimentelle Verifikation**: Mehrere Studien zeigen, dass das Gompertz-Modell experimentelle Daten am besten anpasst 3. **Parameterkorrelation**: Die beiden Parameter des Gompertz-Modells zeigen signifikante Korrelation ### Tumor-Immun-Wechselwirkungsmodelle 1. **Kuznetsov et al. (1994)**: Etablierung des grundlegenden Tumor-Immun-Modellrahmens 2. **Li et al. (2021,2022)**: Untersuchung der Auswirkungen einfacher und doppelter Verzögerungen 3. **Beitrag dieses Artikels**: Erweiterung auf verallgemeinertes Richards-Wachstum und periodische Chemotherapieeffekte ### Theorie der Verzögerungsdifferentialgleichungen 1. **Stabilitätsanalyse**: Charakteristische Gleichungsmethode, Lyapunov-Funktionsmethode 2. **Bifurkationstheorie**: Hopf-Bifurkation, Stabilitätswechsel 3. **Periodische Lösungen**: Fortsetzungsmethode, Fixpunktsatz ## Schlussfolgerungen und Diskussion ### Hauptschlussfolgerungen 1. **Vollständige Gleichgewichtspunktklassifizierung**: Angabe aller möglichen Existenzbedingungen und Klassifizierung von Gleichgewichtspunkten 2. **Explizite Stabilitätsbedingungen**: Bereitstellung von expliziten Stabilitätskriterien und kritischen Verzögerungswerten 3. **Existenz periodischer Lösungen**: Beweis der Existenz periodischer Lösungen unter periodischer Chemotherapiewirkung 4. **Klare Bifurkationsphänomene**: Offenlegung des Übergangsmechanismus des Systems von Stabilität zu Oszillation ### Einschränkungen 1. **Parametervereinfachung**: Die Annahme $a=0$, obwohl biologisch begründet, schränkt die Allgemeinheit ein 2. **Unzureichende numerische Verifikation**: Theoretische Ergebnisse mangelt es an ausreichender numerischer Simulation 3. **Biologische Verifikation**: Modellvorhersagen benötigen mehr experimentelle Datenunterstützung ### Zukünftige Richtungen 1. **Dreidimensionale Erweiterung**: Berücksichtigung weiterer Zelltypen (z.B. regulatorische T-Zellen) 2. **Räumliche Heterogenität**: Einführung von Diffusionstermen zur Untersuchung räumlicher Muster 3. **Stochastische Effekte**: Berücksichtigung von Rauscheinflüssen auf die Systemdynamik 4. **Optimale Kontrolle**: Entwurf optimaler Chemotherapieschemata auf Grundlage des Modells ## Tiefgreifende Bewertung ### Stärken 1. **Theoretische Strenge**: Vollständige mathematische Analyse mit klaren Beweisen 2. **Angemessenes Modell**: Klarer biologischer Hintergrund und verständliche Parameterbedeutung 3. **Methodische Innovation**: Kombination von doppelter Verzögerung und verallgemeinerter Wachstumsfunktion 4. **Praktische Ergebnisse**: Theoretische Anleitung für die Gestaltung von Chemotherapieschemata ### Mängel 1. **Schwacher numerischer Teil**: Abschnitt 6 "Numerische Verifikation" ist als "in Bearbeitung" gekennzeichnet 2. **Fehlende Parameterschätzung**: Keine Methoden zur Parameterschätzung aus realen Daten angegeben 3. **Unzureichende Modellverifikation**: Fehlender Vergleich mit klinischen Daten ### Einfluss 1. **Theoretischer Beitrag**: Bereicherung der Anwendung von Verzögerungsdifferentialgleichungen in der Biomathematik 2. **Anwendungswert**: Neue Perspektiven für die mathematische Modellierung von Tumorbehandlung 3. **Methodologische Bedeutung**: Allgemeine Anwendbarkeit der Fortsetzungsmethode in der Analyse periodischer Lösungen ### Anwendungsszenarien 1. **Theoretische Forschung**: Verzögerungsdifferentialgleichungen, Theorie dynamischer Systeme 2. **Biomathematik**: Tumorwachstumsmodelle, Immunsystemmodellierung 3. **Medizinische Anwendung**: Optimierung von Chemotherapieschemata, Vorhersage von Behandlungsergebnissen ## Literaturverzeichnis Der Artikel zitiert 27 relevante Referenzen, hauptsächlich bestehend aus: - Klassischen Literaturstellen zu Tumorwachstumsmodellen - Stabilitätstheorie von Verzögerungsdifferentialgleichungen - Mathematischen Modellen der Tumor-Immun-Wechselwirkung - Bifurkationstheorie und Existenzsätzen für periodische Lösungen --- **Gesamtbewertung**: Dies ist ein theoretisch sehr starkes Papier der mathematischen Biologie, das wichtige Beiträge zur Theorie der Verzögerungsdifferentialgleichungen und zur Modellierung von Tumor-Immunsystem-Wechselwirkungen leistet. Die theoretische Analyse ist umfassend und die Methoden innovativ, aber der numerische Verifikationsteil bedarf noch der Verbesserung. Es hat große Bedeutung für das Verständnis der komplexen Dynamik von Tumor-Immunsystem-Wechselwirkungen.