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
  • शीर्षक: समय-विलंब और कीमोथेरेपी प्रभावों के साथ ट्यूमर-प्रतिरक्षा प्रतिस्पर्धा प्रणाली में आवधिक समाधान
  • लेखक: पाब्लो एम्स्टर, एंड्रेस रिवेरा, जॉन ए. अरेडोंडो
  • वर्गीकरण: 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. **कुज़नेत्सोव आदि (1994)**: ट्यूमर-प्रतिरक्षा मॉडल ढांचा की स्थापना 2. **ली आदि (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 संबंधित संदर्भों का हवाला देता है, मुख्य रूप से शामिल हैं: - ट्यूमर वृद्धि मॉडल के शास्त्रीय साहित्य - विलंब अवकल समीकरण स्थिरता सिद्धांत - ट्यूमर-प्रतिरक्षा परस्पर क्रिया के गणितीय मॉडल - विभाजन सिद्धांत और आवधिक समाधान अस्तित्व प्रमेय --- **समग्र मूल्यांकन**: यह विलंब अवकल समीकरण और ट्यूमर-प्रतिरक्षा प्रणाली मॉडलिंग के क्षेत्र में महत्वपूर्ण योगदान देने वाला एक सैद्धांतिक रूप से मजबूत गणितीय जीव विज्ञान पेपर है। सैद्धांतिक विश्लेषण विस्तृत है, विधि नवीन है, लेकिन संख्यात्मक सत्यापन भाग में सुधार की आवश्यकता है। ट्यूमर-प्रतिरक्षा प्रणाली की जटिल गतिशीलता को समझने के लिए महत्वपूर्ण महत्व रखता है।