Next Article in Journal
Fundamental Characteristics of the Product-Operated Metric Spaces
Next Article in Special Issue
Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method
Previous Article in Journal
Neutrosophic Sub-Group over t-Norm and t-Co-Norm
Previous Article in Special Issue
Study of Rotavirus Mathematical Model Using Stochastic and Piecewise Fractional Differential Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis in a Mathematical Model for Allergic Reactions

1
Department of Mathematics and Informatics, University Politehnica of Bucharest, 060042 Bucharest, Romania
2
Department of Mathematical Methods and Models, University Politehnica of Bucharest, 060042 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(2), 102; https://doi.org/10.3390/axioms13020102
Submission received: 27 December 2023 / Revised: 25 January 2024 / Accepted: 28 January 2024 / Published: 31 January 2024

Abstract

:
We present a mathematical model that captures the dynamics of the immune system during allergic reactions. Using delay differential equations, we depict the evolution of T cells, APCs, and IL6, considering cell migration between various body compartments. The biological discussions and interpretations within the article revolve around drug desensitization, highlighting one potential application of the model. We conduct stability analysis on certain equilibrium points, demonstrating stability in some cases and only partial stability in others. Numerical simulations validate the theoretical findings.

1. Introduction

An allergic response stems from an exaggerated reaction of the immune system. When the body encounters an allergen, it triggers the production of a substantial amount of IgE (immunoglobulin E) antibodies, responsible for instigating allergic reactions. Allergy symptoms can vary widely, from itching, rash, swelling, and cramps to the severe and potentially life-threatening anaphylactic shock.
Among the myriad potential allergens, drug-induced allergic reactions pose significant challenges. As remarked in [1], the widespread use of chemotherapy drugs has escalated hypersensitivity reactions (HSRs) globally. These specific allergic responses can endanger lives, restricting the use of primary medications and jeopardizing both patient survival and quality of life. Re-exposure to allergic-induced medications has tragically resulted in fatalities. A study of such drug-induced allergies was carried out in [2]. In this paper, we focus on modeling general allergic reactions.
The immune system comprises a complex network of biological elements and processes crucial for defending the body against diseases. Central to this system are white blood cells (lymphocytes), categorized into myeloid cells and lymphoid cells. Myeloid cells encompass various types, such as neutrophils, basophils, dendritic cells, and macrophages. Lymphoid cells consist of T cells, B cells, and natural killer cells.
T cells can be further categorized into helper T cells, memory T cells, cytotoxic T cells, and regulatory T cells. Regulatory T cells (Treg) play a pivotal role in modulating immune system activation and preventing autoimmune responses.
Antigen-presenting cells (APCs) aid in activating T cells. They are classified as professional APCs (including dendritic cells, macrophages, and B cells) or non-professional APCs (found in various nucleated cell types within the body). Upon contact with an allergen, APCs capture the allergen and present it to naive T helper cells.
The immune system relies significantly on cytokines, small proteins that coordinate immune responses. Interleukins, a type of cytokine, facilitate communication among white blood cells. These cytokines, produced predominantly by white blood cells, regulate the differentiation of other cells and can either stimulate or inhibit their functions. Notably, the same cytokine may be produced by various cell types.
Immunoglobulins, or antibodies, are generated by plasma cells derived from mature B cells. There are five primary antibody classes—IgA, IgD, IgE, IgG, and IgM—each with distinct structures, functions, and responses to antigens. Elevated levels of IgE are associated with allergic reactions.
Maintaining immunological balance is crucial in managing allergies, which involves a well-balanced Type 1 T helper (Th1) and Type 2 T helper (Th2) response. T helper cells, major cytokine producers, influence this balance. Certain Th2-type cytokines correlate with IgE production, suggesting that an imbalance toward Th2 may lead to allergies [3].
Mature T cells exhibit autocrine action, regulating their population and suppressing other cells via their cytokines. Th1 and Th2 cells suppress each other, while Treg cells suppress both Th1 and Th2 cells without being suppressed themselves.
Interleukin-6 (IL-6) is a cytokine known for promoting Th2 cell differentiation and inhibiting Th1 and Treg growth and differentiation in response to allergens [4]. It can be produced by both professional antigen-presenting cells (APCs) and T helper cells [4,5]. The authors of [6] were the first to demonstrate the influence of IL-6 on tilting the Th1/Th2 balance toward Th2.
Drug desensitization, a process for inducing temporary immunological tolerance in individuals hypersensitive to certain drugs, has proven safe and effective in improving clinical outcomes. This method involves gradually administering increasing doses of the drug [7,8]. Rapid IgE desensitization inhibits IL-6 [8], and studies [9] showed a 75% reduction in IL-6 production post-desensitization.
According to [10], the body consists of two compartments: the central compartment comprising the blood and well-perfused organs, and the peripheral compartment encompassing poorly perfused tissues and organs.
As argued in the paper by Fishman and Segel [11], effective immunotherapy for allergies necessitates considering at least two compartments within the peripheral immune system: one centered on mucosal lymph nodes (referred to as compartment- μ ) and another on non-mucosal lymphatic tissues draining the immunotherapeutic injection sites (compartment- η in [11]).
Based on this analysis, our model introduces three compartments: the central compartment and the peripheral compartment divided into mucosal and non-mucosal compartments.
Referring to [11], T cells initially recruited into the η -compartment will transit through lymphatic tissues in the μ -compartment and vice versa.
The mathematical modeling is performed through delayed differential equations (DDEs). Through these types of equations, we can take into account the cell migration time or the time required for the production of antibodies (which are not instantaneous). Thus, we obtain a more accurate description of the phenomena we want to model.
The paper’s structure is as follows: Section 2 introduces and elaborates on the mathematical model; Section 3 explores the relevant solutions’ properties; Section 4 investigates the biologically significant equilibria and their stability; Section 5 presents numerical simulations with biological interpretations; Section 6 lists the parameter values used in the simulations; and Section 7 concludes the article.

2. The Model Establishment

The model we propose comprises 10 delay differential equations, building upon the models presented in [12,13,14,15]. A notable difference lies in our consideration of three body compartments and the introduction of three delays to accurately capture the biological timeframe.
The initial four equations depict the dynamics within the non-mucosal compartment ( η -compartment) following allergen entry, describing the concentrations of naive T cells, Th1, Th2, and T r e g cells denoted as N η , T 1 η , T 2 η , and T r η , respectively.
Subsequently, Equations (5)–(7) detail the behavior within the mucosal compartment ( μ -compartment) concerning the concentrations of Th1, Th2, and T r e g cells ( T 1 μ , T 2 μ , T r μ ).
Equations (8) and (9) elucidate a two-stage maturation process of APCs post-allergen contact, where A 1 signifies the concentration of naive APCs and A 2 represents mature APCs. The final equation characterizes the production of IL-6, denoted by I.
Therefore, the model is as follows:
N ˙ η = α γ 3 N η N η A 2 ( t τ 2 ) T 1 η 1 + μ 2 T 2 η ϕ N η A 2 ( t τ 2 ) T 2 η κ N η A 2 ( t τ 2 ) T r η T ˙ 1 η = ( 1 + θ ) T 1 η + θ T 1 μ ( t τ 3 ) + + v N η A 2 ( t τ 2 ) ( 1 + μ r T r η ) T 1 η 1 + μ 2 T 2 η η 1 I T 1 η 1 + I T ˙ 2 η = ( 1 + θ ) T 2 η + θ T 2 μ ( t τ 3 ) + + ϕ v N η A 2 ( t τ 2 ) ( 1 + μ r T r η ) T 2 η 1 + μ 1 T 1 η 1 + μ 2 T 2 η + η 2 I T 2 η 1 + I T ˙ r η = ( 1 + θ ) T r η + θ T r μ ( t τ 3 ) + κ v N η A 2 ( t τ 2 ) T r η η r I T r η 1 + I T ˙ 1 μ = ( 1 + θ ) T 1 μ + θ T 1 η ( t τ 3 ) T ˙ 2 μ = ( 1 + θ ) T 2 μ + θ T 2 η ( t τ 3 ) T ˙ r μ = ( 1 + θ ) T r μ + θ T r η ( t τ 3 ) A ˙ 1 = λ γ 1 A 1 β Λ A 1 A ˙ 2 = β Λ A 1 γ 2 A 2 μ A 2 T r η I ˙ = γ 4 I + k 1 [ A 2 ( t τ 1 ) + N η ( t τ 1 ) + T 1 η ( t τ 1 ) + + T 2 η ( t τ 1 ) + T r η ( t τ 1 ) ]
The initial equation represents the variation in the concentration of naive T cells characterized by a constant production rate α and a constant death rate γ 3 . The last three terms signify the differentiation of naive cells into T h 1 , T h 2 , and T r e g cells upon encountering mature APCs ( A 2 ).
The subsequent three equations illustrate the evolution of the concentrations of T h 1 , T h 2 , and T r e g in the η -compartment. Following [14], these concentrations are related to naive cells, the presented allergen concentration, and their respective cytokines.
In each equation, the first two terms account for migration between the μ - and η -compartments, cell degradation, and the third term indicates the differentiation of naive cells into T h 1 , T h 2 , and T r e g . Notably, for T h 1 , differentiation is suppressed by T r e g and T h 2 cells, while for T h 2 , suppression occurs from T r e g and T h 1 cells. The last term in each equation represents the effect of I L 6 during chemotherapy with desensitization dose Λ . Additionally, the equations for T h 1 and T r e g incorporate inhibition rates η 1 and η r , while the equation for T h 2 includes a stimulation rate η 2 [6,16,17].
Parameters like v determine the proliferation rate of differentiated T cells from a single naive cell, while ϕ and κ account for autocrine action differences among mature T helper cell subsets. The suppression strengths of T h 1 , T h 2 , and T r e g are controlled by parameters μ 1 , μ 2 , and μ r .
Equations (5)–(7)model concentrations of T h 1 , T h 2 , and T r e g cells in the μ -compartment, considering migration to and from the η -compartment with migration probability ( θ ).
Equations (8) and (9) portray the progression of naive and mature APCs post-allergen contact. The chemotherapeutic drug-induced allergen Λ remains constant. According to [13], T r e g cells play a regulatory role in limiting APC maturation.
Naive APCs exhibit a constant production rate λ and a constant death rate γ 1 . The third term in Equation (9) represents APC activation by the antigen.
Equation (10) illustrates the production of IL-6, where mature APCs and mature T-cells, following [4], contribute to its production. The first term denotes the death rate of I L 6 cells, while the second term depicts I cell production by mature APCs, Naive T cells, Th1 cells, Th2 cells, and Treg cells.
To ensure biological accuracy, three delays are introduced as follows:
  • τ 1 is the time necessary for the production of IL-6;
  • τ 2 is the allergen’s propagation time from the central compartment to the peripheral compartment [18];
  • τ 2 is the transit time due to the inter-compartmental migration of Th1, Th2 and Treg cells between the mucosal and non-mucosal compartment. The delay is prescribed by the probability of migration ( k m ):
    τ 3 = 1 k m
A comprehensive list detailing the model’s parameters and their descriptions is available in Section 6.

Introducing New Notations for State Variables

Initially, the model is introduced using common notations frequently found in the literature, aiming for ease of comprehension from a biological perspective. To facilitate the mathematical analysis of the delay differential equation (DDE) system, subsequently, the following notations are introduced:
  • x 1 = concentration of naive T-cells ( N η ).
  • x 2 = concentration of Th1 cells in η -compartment.
  • x 3 = concentration of Th2 cells in η -compartment.
  • x 4 = concentration of T r e g cells in η -compartment.
  • x 5 = concentration of Th1 cells in μ -compartment.
  • x 6 = concentration of Th2 cells in μ -compartment.
  • x 7 = concentration of T r e g cells in μ -compartment.
  • x 8 = concentration of naive APCs.
  • x 9 = concentration of mature APCs.
  • x 10 = concentration of IL-6.
Thus, (1) becomes:
x ˙ 1 = α γ 3 x 1 x 1 x 9 τ 2 x 2 1 + μ 2 x 3 ϕ x 1 x 9 τ 2 x 3 κ x 1 x 9 τ 2 x 4 x ˙ 2 = ( 1 + θ ) x 2 + θ x 5 τ 3 + v x 1 x 9 τ 2 x 2 ( 1 + μ r x 4 ) ( 1 + μ 2 x 3 ) η 1 x 10 x 2 1 + x 10 x ˙ 3 = ( 1 + θ ) x 3 + θ x 6 τ 3 + ϕ v x 1 x 9 τ 2 x 3 ( 1 + μ 2 x 3 ) ( 1 + μ r x 4 ) ( 1 + μ 1 x 2 + μ 2 x 3 ) + η 2 x 10 x 3 1 + x 10 x ˙ 4 = ( 1 + θ ) x 4 + θ x 7 τ 3 + κ v x 1 x 9 τ 2 x 4 η r x 10 x 4 1 + x 10 x ˙ 5 = ( 1 + θ ) x 5 + θ x 2 τ 3 x ˙ 6 = ( 1 + θ ) x 6 + θ x 3 τ 3 x ˙ 7 = ( 1 + θ ) x 7 + θ x 4 τ 3 x ˙ 8 = λ γ 1 x 8 β Λ x 8 x ˙ 9 = β Λ x 8 γ 2 x 9 μ x 9 x 4 x ˙ 10 = γ 4 x 10 + k 1 ( x 9 τ 1 + x 1 τ 1 + x 2 τ 1 + x 3 τ 1 + x 4 τ 1 )

3. General Properties of the Solutions

Because the model describes a biological mechanism, a critical step is to demonstrate the positivity of the solutions.
Define τ = m a x { τ 1 , τ 2 , τ 3 } and let P C ( [ τ , 0 ] , R 10 ) denote the space of piecewise continuous functions defined on [ τ , 0 ] with values in R n . The norm in P C ( [ τ , 0 ] , R 10 ) will be defined by
φ τ = s u p { φ ( t ) 2 ; t [ τ , 0 ] }
with · 2 the Euclidean norm in R n . For (2), consider the initial data
x ( s ) = φ ( s ) , s [ τ , 0 ]
Proposition 1. 
If the initial data φ P C ( [ τ , 0 ] , R 10 ) satisfy φ j ( s ) > 0 for s [ τ , 0 ] , j = 1 , , 10 , then the solution to the Cauchy problem (2) + (3) will fulfill x j ( t ) 0 for all t within the domain of existence.
Proof. 
Suppose x j ( t ) > 0 t [ τ , t 0 ) . Then, if x 1 ( t 0 ) = 0 , it follows that x 1 ( t 0 ) = α > 0 so x 1 will increase for t > t 0 . If x 2 ( t 0 ) = 0 , then x 2 ( t 0 ) = θ x 5 ( t τ 3 ) > 0 so x 2 increases for t > t 0 . The same reasoning applies to x 3 , , x 10 . □
From now on, the initial data for (2) will be supposed positive.
Proposition 2. 
x 1 , x 8 , x 9 are bounded on the whole interval of existence.
Proof. 
From (2), it follows that
x 1 ˙ ( t ) = α γ 3 x 1 ( t ) x 1 ( t ) p ( t ) .
with p ( t ) 0 for positive initial data. Then,
x 1 ( t ) = x 1 ( 0 ) e γ 3 t 0 t p ( s ) d s + α 0 t e γ 3 s e 0 s p ( r ) d r d s e γ 3 t e 0 t p ( s ) d s .
The following estimation for the second term holds:
α 0 t e γ 3 s e 0 s p ( r ) d r d s e γ 3 t e 0 t p ( s ) d s α 0 t e γ 3 s e 0 t p ( r ) d r d s e γ 3 t e 0 t p ( s ) d s = α γ 3 ( 1 e γ 3 t ) α γ 3 , t 0 .
It follows that | x 1 ( t ) | M 1 for some positive M 1 .
It is elementary to see that
x 8 ( t ) = x 8 ( 0 ) λ γ 1 + β Λ e ( γ 1 + β Λ ) t + λ γ 1 + β Λ ,
so | x 8 ( t ) | M 8 , t [ τ , ) for some M 8 > 0 .
Concerning x 9 , the variation of constants formula gives
x 9 ( t ) = x 9 ( 0 ) e γ 2 t μ 0 t x 4 ( s ) d s + β Λ 0 t x 8 ( s ) e γ 2 s + μ 0 s x 4 ( r ) d r d s e γ 2 t μ 0 t x 4 ( s ) d s .
The first term is obviously bounded for t 0 . For the second one, we have the following estimations:
β Λ | 0 t x 8 τ 1 ( s ) e γ 2 s + 0 s x 4 ( r ) d r d s | e γ 2 t 0 t x 4 ( s ) d s β Λ 0 t | x 8 ( s ) | e γ 2 s + 0 s x 4 ( r ) d r d s e γ 2 t 0 t x 4 ( s ) d s β Λ M 8 0 t e γ 2 s + 0 s x 4 ( r ) d r d s e γ 2 t e 0 t x 4 ( s ) d s = = β Λ M 8 e γ 2 t 1 γ 2 e γ 2 t β Λ M 8 γ 2 ,
so | x 9 ( t ) | M 9 t [ τ , ) with M 9 > 0 . □
Proposition 3. 
The solution of (2) exists on [ τ , ) .
Proof. 
The proposition follows from a slight generalization of Theorem 1.2. in [19]. Namely, the condition imposed to the right-hand side of the equation should hold only for the solutions of the equation instead of all functions from P C [ τ , 0 ] ) . So we must show that, with φ = ( φ 1 , , φ 10 ) , a solution of (2) and f = ( f 1 , , f 10 ) , the right-hand side of (2),
| f j ( φ ) | h ( φ τ ) , r 0 1 h ( r ) d r = , r 0 > 0 .
We will show that there exist constants C 1 , C 2 so that
| f j ( φ ) | C 1 + C 2 φ τ , j = 1 , , 10 .   | f 1 ( φ ) | | α | + γ 3 | φ 1 ( t ) | + M 1 M 9 | φ 2 ( t ) | + ϕ M 1 M 9 | φ 3 ( t ) | + κ M 1 M 9 | φ 4 ( t ) |   | α | + ( γ 3 + M 1 M 9 + ϕ M 1 M 9 + κ M 1 M 9 ) φ τ ,   | f 2 ( φ ) | ( 1 + θ ) | φ 2 ( t ) | + θ | φ 5 ( t τ 3 ) | + ( v M 1 M 9 + η 1 ) | φ 2 ( t ) |   ( 1 + 2 θ + v M 1 M 9 + η 1 ) φ τ ,   | f 3 ( φ ) | ( 1 + θ ) | φ 3 ( t ) | + θ | φ 6 ( t τ 3 ) | + ( ϕ v M 1 M 9 + η 2 ) | φ 3 ( t ) |   + ( 1 + 2 θ + ϕ v M 1 M 9 + η 2 ) φ τ ,   | f 4 ( φ ) | ( 1 + θ ) | φ 4 ( t ) | + θ | φ 7 ( t τ 3 ) | + ( κ v M 1 M 9 + η r ) | φ 4 ( t ) |   ( 1 + 2 θ + κ v M 1 M 9 + η r ) φ τ ,   | f 5 ( φ ) | ( 1 + θ ) | φ 5 ( t ) | + θ | φ 2 ( t τ 3 ) | ( 1 + 2 θ ) φ τ ,   | f 6 ( φ ) | ( 1 + 2 θ ) φ τ , | f 7 ( φ ) | ( 1 + 2 θ ) φ τ ,   | f 8 ( φ ) | | λ | + ( γ 1 + β Λ ) φ τ ,   | f 9 ( φ ) | ( β Λ + γ 2 + μ M 9 ) φ τ ,   | f 10 ( φ ) | ( γ 4 + 5 k 1 ) φ τ .
The proposition is proved. □

4. Equilibria and Stability Analysis

The system (2) possesses two equilibrium points that are particularly significant for our study:
E 1 = ( x 1 * , 0 , 0 , 0 , 0 , 0 , 0 , 0 , x 8 * , x 9 * , x 10 , 10 * )
E 2 = ( x ^ 1 , x ^ 2 , 0 , 0 , x ^ 5 , 0 , 0 , x ^ 8 , x ^ 9 , x ^ 10 )
Both equilibrium points hold positive implications for patients. In the scenario of E 1 , there exists no immune response to the allergen, while for E 2 , the immune response lacks association with allergies, evidenced by x 2 > x 3 ( T h 1 > T h 2 ). The stability of these equilibrium points is pivotal for patient well-being.
Subsequently, we analyze the stability properties of these equilibrium points.

4.1. Equilibrium Point E 1

Let E 1 = ( x 1 * , 0 , 0 , 0 , 0 , 0 , 0 , 0 , x 8 * , x 9 * , x 10 , 10 * ) with
x 1 * = α , x 8 * = λ γ 1 + β Λ , x 9 * = β Λ λ γ 2 ( γ 1 + β Λ ) , x 10 * = k 1 x 9 * γ 4
The matrix A, comprising partial derivatives concerning the undelayed variables and computed at E 1 , exhibits the following non-zero elements:
a 11 = γ 3 , a 12 = x 1 * x 9 * , a 13 = ϕ x 1 * x 9 * , a 14 = k x 1 * x 9 * a 22 = ( 1 + θ ) + v x 1 * x 9 * η 1 x 10 * 1 + x 10 * a 33 = ( 1 + θ ) + ϕ x 1 * x 9 * + η 2 x 10 * 1 + x 10 * a 44 = ( 1 + θ ) + k v x 1 * x 9 * η r x 10 * 1 + x 10 * a 55 = a 66 = a 77 = ( 1 + θ ) , a 88 = γ 1 β Λ a 94 = μ x 9 * , a 98 = β Λ , a 99 = γ 2 , a 10 , 10 = γ 4 .
The matrix A τ 1 of partial derivatives concerning the τ 1 -delayed variables contains the only non-zero elements.
b 10 , 1 = b 10 , 2 = b 10 , 3 = b 10 , 4 = b 10 , 9 = k 1 & b 10 , 10 = γ 4 .
The matrix A τ 2 pertaining to partial derivatives concerning the τ 2 -delayed variables is zero. Lastly, the matrix A τ 3 , representing partial derivatives concerning the τ 3 -delayed variables, contains the nonzero elements:
d 25 = d 36 = d 47 = d 52 = d 63 = d 74 = θ .
The characteristic equation for E 1 is
d e t ( λ I 10 A 1 3 A τ j e λ τ j ) = 0 ( λ a 10 , 10 ) ( λ a 99 ) ( λ a 88 ) ( λ a 11 ) d 1 ( λ ) = 0
d 1 ( λ ) = λ a 22 0 0 θ e λ τ 3 0 0 0 λ a 33 0 0 θ e λ τ 3 0 0 0 λ a 44 0 0 θ e λ τ 3 θ e λ τ 3 0 0 λ a 55 0 0 0 θ e λ τ 3 0 0 λ a 66 0 0 0 θ e λ τ 3 0 0 λ a 77
To address the challenges associated with analyzing the latter part of the characteristic equation, we will employ an alternative approach for the case where d 1 = 0 . It is noteworthy that the equation d 1 ( λ ) = 0 represents the characteristic equation of a linear system:
x ˙ = A 1 x + B 1 x τ 3
where A 1 is diagonal, and B 1 follows the Meltzer structure (see [20]). One can then apply the criterion from [20] as presented in Proposition 2.3: If the matrix A 1 + B 1 is Hurwitz, and if the delay τ 3 is sufficiently small, the system (4) is exponentially stable. Consequently, it follows that if
a 10 , 10 < 0 , a 99 < 0 , a 88 < 0 , a 11 < 0 ,
and this is indeed the case, then the equilibrium point E 1 is asymptotically stable.
Below, sufficient conditions ensuring that A 1 + B 1 is Hurwitz are provided. It is worth noting that the polynomial derived from d 1 when τ 3 = 0 should be Hurwitz, meaning the roots of the corresponding equation must have negative real parts.
Now, when τ 3 = 0 , the equation d 1 ( λ ) = 0 transforms into:
( λ a 33 ) ( λ a 44 ) ( λ a 66 ) ( λ a 77 ) θ 2 [ ( λ a 33 ) ( λ a 66 ) + ( λ a 44 ) ( λ a 77 ) ] + θ 4 = 0
To address Equation (6), we introduce
u 1 = ( λ a 33 ) ( λ a 66 ) , u 2 = ( λ a 44 ) ( λ a 77 )
so, with s = θ 2 , Equation (6) becomes
s 2 s ( u 1 + u 2 ) + u 1 u 2 = 0
with solutions s 1 = u 1 , s 2 = u 2 .
Consequently, the roots of the following equations must have real parts smaller than zero:
λ 2 ( a 33 + a 66 ) λ + a 33 a 66 θ 2 = 0 , λ 2 ( a 44 + a 77 ) λ + a 44 a 77 θ 2 = 0 .
The following necessary and sufficient conditions ensure that Equation (8) possess roots with negative real parts:
a 33 + a 66 < 0 , a 33 a 66 θ 2 > 0 , a 44 + a 77 < 0 , a 44 a 77 θ 2 > 0
The following theorem results.
Theorem 1. 
Assuming conditions (5) and (9) are satisfied, then, for
τ 3 < 1 e | a j j | , j = 2 , , 7
the equilibrium point E 1 is asymptotically stable.
Proof. 
The conclusion is reached by applying Proposition 2.3. from [20] and the preceding reasoning. □

4.2. Equilibrium Point E 2

Let us now consider the equilibrium point
E 2 = ( x ^ 1 , x ^ 2 , 0 , 0 , x ^ 5 , 0 , 0 , x ^ 8 , x ^ 9 , x ^ 10 ) .
The values for E 2 are derived from the following equations:
x ^ 8 = λ γ 1 + β Λ x ^ 9 = λ β Λ γ 2 ( γ 1 + β Λ ) x ^ 2 = α γ 3 x ^ 1 x ^ 1 x ^ 9 x ^ 10 = k 1 γ 4 ( x ^ 9 + x ^ 1 + α γ 3 x ^ 1 x ^ 1 x ^ 9 ) x ^ 5 = θ x ^ 2 1 + θ
Substitute these values into the second Equation of (2) and simplify x ^ 2 . Let c 1 = 1 + 2 θ 1 + θ . The subsequent equation yields x ^ 1 :
k 1 v x ^ 9 2 x 3 + ( v γ 4 x ^ 9 2 + k 1 v x ^ 9 3 k 1 v γ 3 x ^ 9 c k 1 x ^ 9 η 1 k 1 x ^ 9 ) x 2 + + ( k 1 v α x ^ 9 c γ 4 x ^ 9 c k 1 x ^ 9 2 + c k 1 γ 3 η 1 k 1 x ^ 9 2 + η 1 k 1 γ 3 ) x c k 1 α η 1 k 1 α = 0
(10) can be expressed as:
a x 3 + b x 2 + c x + d = 0 .
Calculate, as previously done, the partial derivatives in E 2 , and denote
b 2 = ϕ v x ^ 1 x ^ 9 1 + μ 1 x ^ 2 + η 2 x ^ 10 1 + x ^ 10 ,
b 3 = κ v x ^ 1 x ^ 9 η r x ^ 10 1 + x ^ 10 .
Revisit now the following definitions from [21]. Consider the system
x ˙ = X ( t , x t ) , x t 0 = φ
Assume that τ > 0 is a given real number, x t : [ τ , 0 ] R n is defined by x t ( s ) = x ( t + s ) , with the norm | | x ( t ) | | 2 = x 1 2 ( t ) + + x n 2 ( t ) . Introduce the partition
x = ( z , u ) T (T denotes transposition), where z R m , u R n m ( 1 m < n ) .
X : R × C R n is assumed to be continuous and X maps every R × (bounded set) into a bounded set, in the domain D defined by
t 0 , | | z | | τ h , | | u | | τ <
The solutions of (11) are assumed to be unique and uncontinuable: the solutions are defined for every t t 0 where | | z ( t ; t 0 , φ ) | | 2 h . If we suppose that (12) holds for all t t 0 , then the solution x ( t ; t 0 , φ ) is defined for all t t 0 .
Taking into account the above partitions, the system (11) can be represented as:
z ˙ ( t ) = Z ( t , y t , z t ) , u ˙ ( t ) = U ( t , y t , z t )
Suppose that
Z ( t , 0 , 0 ) = 0 , U ( t , 0 , 0 ) = 0 t t 0
Let x ( t ; t 0 , φ ) denote a solution of the system (11) with initial condition x ( t 0 ; φ ) .The same notation will be applied for the solutions of each of the two equations in (14).
Definition 1 
([21,22,23]). The equilibrium point x = 0 in system (11), (13) is called the following:
1. 
z-stable, if for every t 0 0 and every ϵ > 0 , there exists a δ ( ϵ ; t 0 ) > 0 such that | | φ | | τ < δ implies | | z ( t ; t 0 ; φ ) | | 2 < ϵ for all t t 0 . It is called uniformly z-stable if δ does not depend on t 0 .
2. 
Asymptotically z-stable if it is z-stable and for every t 0 0 there exists a Δ ( t 0 ) > 0 such that for every solution x ( t ; t 0 ; φ ) of system (11), (13) that satisfies | | φ | | τ < Δ ( t 0 ) , the following holds true:
lim t | | z ( t ; t 0 ; φ ) | | = 0
3. 
Uniformly asymptotically z-stable, if it is uniformly z-stable with respect to t 0 in terms of point 1, and one can find Δ > 0 such that relation (15) is met uniformly with respect to ( t 0 , φ ) from the domain t 0 0 , | | φ | | τ < Δ (for any numbers η > 0 , t 0 0 one can find the number T = T ( η ) > 0 such that | | z ( t ; t 0 , φ ) | | < η for all t t 0 + T , | | φ | | τ < Δ )
The following theorem will be used to prove the stability of E 2 with respect to some of its variables.
Theorem 2 
([21], Th.5.2.1). Suppose there exists a function V : R × C R of class C 1 and continuously increasing functions.
a , b : R + R + with a ( 0 ) = b ( 0 ) = 0 so that, in the domain (12),
V ( t , φ ) ) a ( | | φ v ( 0 ) ) | | ) , V ˙ ( t , x t ) 0 .
where V ˙ denotes the derivative along system (11), then the equilibrium point x = 0 of system (11) is z-stable.
If, in addition,
V ( t , φ ) ) b ( | | φ | | τ ) ,
then the z-stability is uniform.
Theorem 3. 
Suppose b 2 < 1 , b 3 < 1 and β Λ 2 < γ 2 . Then the equilibrium point E 2 exhibits partial stability concerning the variables ( x 1 , x 3 , x 4 , x 6 , x 7 , x 8 , x 9 ) .
Proof. 
Shift the equilibrium to zero by utilizing z i = x i x ^ i , where i = 1 , , 10 (assuming x ^ 3 = 0 , x ^ 4 = 0 , x ^ 6 = 0 , x ^ 7 = 0 ). Introduce z = ( z 1 , , z 10 ) , and define
V ( z ) = i = 3 4 z i 2 2 + θ 2 t τ 3 t z i ( s ) 2 d s + i = 6 7 z i 2 2 + θ 2 t τ 3 t z i ( s ) 2 d s
+ z 1 2 2 + z 8 2 2 + z 9 2 2
Denote as y = ( z 1 , z 3 , z 4 , z 6 , z 7 , x 8 , x 9 ) = ( x 1 , x 3 , x 4 , x 6 , x 7 , x 8 , x 9 ) Additionally, it should be noted that the initial condition to establish a (partial) Lyapunov–Krasovskii functional is confirmed by V:
V ( z ) v ( | | y | | 2 ) .
Also,
V ( z ) M | | z | | τ 2
for some M 1 > 0 . The derivative of V along the adjusted system (11) is expressed as:
d V d t = z 3 f 3 + z 4 f 4 + z 6 f 6 + z 7 f 7 + θ 2 ( z 3 2 z 3 τ 3 2 + z 4 2 z 4 τ 3 2 + z 6 2 z 6 τ 3 2 + z 7 2 z 7 τ 3 2 ) + z 1 f 1 + z 8 f 8 + z 9 f 9 = ( 1 b 2 ) z 3 2 ( 1 b 3 ) z 4 2 ( 1 + θ ) z 6 2 ( 1 + θ ) z 7 2 θ 2 ( z 3 2 2 z 3 z 6 τ 3 + z 6 τ 3 2 ) ( 1 + θ ) z 5 2 θ 2 ( z 4 2 2 z 4 z 7 τ 3 + z 7 τ 3 2 ) θ 2 ( z 6 2 2 z 6 z 3 τ 3 + z 3 τ 3 2 ) θ 2 ( z 7 2 2 z 7 z 4 τ 3 + z 4 τ 3 2 ) + z 3 2 ϕ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( 1 + μ 2 z 3 ) ( 1 + μ r z 4 ) ( 1 + μ 2 z 3 + μ 1 ( z 2 + x ^ 2 ) ) ϕ v x ^ 1 x ^ 9 1 + μ x ^ 2 + η 2 z 10 ( 1 + z 10 + x ^ 10 ) ( 1 + x ^ 10 ) + z 4 2 κ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) η r z 10 + x ^ 10 1 + z 10 + x ^ 10 κ v x ^ 1 x ^ 9 + η r x ^ 10 ( 1 + x ^ 10 ) γ 3 z 1 2 ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( z 2 + x ^ 2 ) z 1 1 + μ 2 z 3 ϕ ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) z 3 z 1 κ ( z 1 + x ^ 1 * ) ( z 9 τ 2 + x ^ 9 ) z 4 z 1 ( γ 1 + β Λ ) z 8 2 + β Λ z 8 z 9 γ 2 z 9 2 μ ( z 9 + x ^ 9 ) z 4 z 9
By removing some negative terms and using a b < a 2 2 + b 2 2 , we obtain
d V d t < ( 1 b 2 ) z 3 2 ( 1 b 3 ) z 4 2 ( 1 + θ ) z 6 2 ( 1 + θ ) z 7 2 θ 2 ( z 3 2 2 z 3 z 6 τ 3 + z 6 τ 3 2 ) ( 1 + θ ) z 5 2 θ 2 ( z 4 2 2 z 4 z 7 τ 3 + z 7 τ 3 2 ) θ 2 ( z 6 2 2 z 6 z 3 τ 3 + z 3 τ 3 2 ) θ 2 ( z 7 2 2 z 7 z 4 τ 3 + z 4 τ 3 2 ) + z 3 2 ϕ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( 1 + μ 2 z 3 ) ( 1 + μ r z 4 ) ( 1 + μ 2 z 3 + μ 1 ( z 2 + x ^ 2 ) ) ϕ v x ^ 1 x ^ 9 1 + μ x ^ 2 + η 2 z 10 ( 1 + z 10 + x ^ 10 ) ( 1 + x ^ 10 ) + z 4 2 κ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) η r z 10 + x ^ 10 1 + z 10 + x ^ 10 κ v x ^ 1 x ^ 9 + η r x ^ 10 ( 1 + x ^ 10 ) γ 3 z 1 2 ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( z 2 + x ^ 2 ) z 1 1 + μ 2 z 3 ϕ ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) z 3 z 1 κ ( z 1 + x ^ 1 * ) ( z 9 τ 2 + x ^ 9 ) z 4 z 1 ( γ 1 + β Λ ) z 8 2 + β Λ z 8 z 9 γ 2 z 9 2 μ ( z 9 + x ^ 9 ) z 4 z 9 = ( 1 b 2 ) z 3 2 ( 1 b 3 ) z 4 2 ( 1 + θ ) z 6 2 ( 1 + θ ) z 7 2 θ 2 ( z 3 2 2 z 3 z 6 τ 3 + z 6 τ 3 2 ) ( 1 + θ ) z 5 2 θ 2 ( z 4 2 2 z 4 z 7 τ 3 + z 7 τ 3 2 ) θ 2 ( z 6 2 2 z 6 z 3 τ 3 + z 3 τ 3 2 ) θ 2 ( z 7 2 2 z 7 z 4 τ 3 + z 4 τ 3 2 ) + z 3 2 ϕ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( 1 + μ 2 z 3 ) ( 1 + μ r z 4 ) ( 1 + μ 2 z 3 + μ 1 ( z 2 + x ^ 2 ) ) ϕ v x ^ 1 x ^ 9 1 + μ x ^ 2 + η 2 z 10 ( 1 + z 10 + x ^ 10 ) ( 1 + x ^ 10 ) + z 4 2 κ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) η r z 10 + x ^ 10 1 + z 10 + x ^ 10 κ v x ^ 1 x ^ 9 + η r x ^ 10 ( 1 + x ^ 10 ) γ 3 z 1 2 ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( z 2 + x ^ 2 ) z 1 1 + μ 2 z 3 ϕ ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) z 3 z 1 κ ( z 1 + x ^ 1 * ) ( z 9 τ 2 + x ^ 9 ) z 4 z 1 ( γ 1 + β Λ ) z 8 2 + β Λ 2 z 8 2 + β Λ 2 z 9 2 γ 2 z 9 2 μ ( z 9 + x ^ 9 ) z 4 z 9 < ( 1 b 2 ) z 3 2 ( 1 b 3 ) z 4 2 ( 1 + θ ) z 6 2 ( 1 + θ ) z 7 2 θ 2 ( z 3 2 2 z 3 z 6 τ 3 + z 6 τ 3 2 ) ( 1 + θ ) z 5 2 θ 2 ( z 4 2 2 z 4 z 7 τ 3 + z 7 τ 3 2 ) θ 2 ( z 6 2 2 z 6 z 3 τ 3 + z 3 τ 3 2 ) θ 2 ( z 7 2 2 z 7 z 4 τ 3 + z 4 τ 3 2 ) + z 3 2 ϕ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) ( 1 + μ 2 z 3 ) ( 1 + μ r z 4 ) ( 1 + μ 2 z 3 + μ 1 ( z 2 + x ^ 2 ) ) ϕ v x ^ 1 x ^ 9 1 + μ x ^ 2 + η 2 z 10 ( 1 + z 10 + x ^ 10 ) ( 1 + x ^ 10 ) + z 4 2 κ v ( z 1 + x ^ 1 ) ( z 9 τ 2 + x ^ 9 ) η r z 10 + x ^ 10 1 + z 10 + x ^ 10 κ v x ^ 1 x ^ 9 + η r x ^ 10 ( 1 + x ^ 10 ) γ 3 z 1 2 ( γ 1 + β Λ 2 ) z 8 2 + ( β Λ 2 γ 2 ) z 9 2 = ω ( | | z t | | τ 2 ) + F ( z t )
with ω being strictly positively defined. It is important to note that, according to Propositions 1 and 2, the terms multiplying x 3 2 and x 4 2 in F are bounded, and therefore,
| F ( z t ) | C | | z t | | τ 3
for some positive C. Therefore, the derivative of V along system (3) is strictly negatively defined for z t < 1 , establishing uniform partial stability as mentioned (see [21,22,23,24,25,26]). □
Remark 1. 
The same reasoning can be applied to demonstrate partial stability for the equilibrium points where x 2 = x 3 = x 5 = x 6 = 0 and x 2 = x 4 = x 5 = x 7 = 0 .
Remark 2. 
The partial asymptotic stability extends to x 8 through direct computation of the solution’s behavior.

5. Numerical Simulations

The numerical simulations were obtained with the values for the parameters given in Section 6. These were taken from the existing literature.
When the stability for E 1 holds, we will have a successful therapy, and it shows that small quantities of allergens do not cause harm. Fortunately, using the above mentioned parameter values, we can show that the equilibrium point E 1 is stable (see Figure 1).
When the stability for E 2 holds, we will have successful immunotherapy of allergies induced by chemotherapy because the Th1 cell population will dominate Th2 cell population (see Figure 2 and Figure 3).

6. List of Parameters

The numerical values for the system’s parameters were taken from the relevant literature (see Table 1).

7. Conclusions

We proposed a mathematical model which describes the evolution of T helper cells, APCs and IL-6 in the case of drug allergies with exposure to drug desensitization during chemotherapy. The novelty of the model consists in a more detailed approach of the biological mechanism. We considered a partition of the body into a central compartment and a peripheral one, where the latter is comprised of the mucosal and the non-mucosal compartments.
The separation of the mucosal and non-mucosal compartments is biologically accurate. It influences the evolution of the cell population, as the migration impacts the number of cells remaining in the non-mucosal compartment.
Another important improvement from the existing literature is the presence of time delays. These help to better capture the correct timeframe of the biological processes. We considered three delays which model the time necessary for the production and migration of different cell populations.
From a mathematical point of view, we studied the existence, boundedness and positivity of solutions. The positivity of the solutions was an essential study to undertake, as we are modeling cell populations. We also gave parameter conditions for stability in the case of some equilibrium points. The stability of the equilibrium points in question illustrate the benefits of drug desensitization. The numerical simulations we performed validate the theoretical results, basically showing that a small amount of allergen does not cause any harm in therapy.
Simply put, this paper introduces a delay differential equations model that describes the cellular changes that occur in the action of the immune system, in the presence of an allergen. The biological insights provided by the model can aid in the desensitization process, a process often necessary in the case of strong allergic reactions or allergies induced by drug therapy. We improved the existing mathematical models on the evolution of the immune system during allergies by offering a more detailed perspective of the biological workings (we introduced time delays and considered that the body is partitioned into two relevant compartments). The stability study performed, along with the numerical simulations and biological interpretations presented, show how one can obtain scenarios for treatment without allergic reactions. The model can be used to find the optimal allergen dose for desensitization and to adjust it accordingly throughout the desensitization process.
Although our focus was on drug-induced allergies, the model is of a versatile nature and can be adapted to describe other types of allergies.

Author Contributions

Validation, I.B.; Formal analysis, A.H.; Writing—original draft, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data is included in the manucripts.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Castells, M. Drug hypersensitivity and anaphylaxis in cancer and chronic inflammatory diseases, the role of desensitizations. Front. Immunol. 2017, 8, 1472. [Google Scholar] [CrossRef]
  2. Abdullah, R.; Badralexi, I.; Halanay, A. Stability Analysis in a New Model for Desenzitization of Allergic Reactions Induced by Chemotherapy of Chronic Lymphocytic Leukemia. Mathemaics 2023, 11, 3225. [Google Scholar] [CrossRef]
  3. Berger, A. Th1 and Th2 responses: What are they? BMJ 2000, 321, 424. [Google Scholar] [CrossRef]
  4. Gubernatorova, E.; Gorshkova, E.; Namakanova, O.; Zvartsev, R.; Hidalgo, J.; Drutskaya, M.; Tumanov, A.; Nedospasov, S. Non-redundant functions of IL-6 produced by macrophages and dendritic cells in allergic airway inflammation. Front. Immunol. 2018, 9, 2718. [Google Scholar] [CrossRef] [PubMed]
  5. Brahmakshatriya, V.; Kuang, Y.; Devarajan, P.; Xia, J.; Zhang, W.; Vong, A.M.; Swain, S.L. IL-6 Production by TLR-Activated APC Broadly Enhances Aged Cognate CD4 Helper and B Cell Antibody Responses In Vivo. J. Immunol. 2017, 198, 2819–2833. [Google Scholar] [CrossRef]
  6. Dienz, O.; Rincon, M. The effects of IL-6 on CD4 T cell responses. Clin. Immunol. 2009, 130, 27–33. [Google Scholar] [CrossRef] [PubMed]
  7. Hong, D. Desensitization for allergic reactions to chemotherapy. Yonsei Med. J. 2019, 60, 119. [Google Scholar] [CrossRef]
  8. Sancho-Serra, M.D.; Simarro, M.; Castells, M. Rapid ige desensitization is antigen specific and impairs early and late mast cell responses targeting fceri internalization. Eur. J. Immunol. 2011, 41, 1004–1013. [Google Scholar] [CrossRef] [PubMed]
  9. Etemadifar, M.; Salari, M.; Saeri, M.; Sigari, A.A.; Ebrahimi, S. Rituximab induced cytokine release syndrome in an MS patient: A case report. Clin. Case Rep. 2021, 9, e04407. [Google Scholar] [CrossRef] [PubMed]
  10. Chen, W. Pharmacokinetic and Pharmacodynamic Modeling of Erythropoiesis. Ph.D Thesis, Drexel University, Philadelphia, PA, USA, 2010; pp. 27–30. [Google Scholar]
  11. Segel, L.A.; Fishman, M.A. Modeling immunotherapy for allergy. Bull. Math. Biol. 1996, 58, 1099–1121. [Google Scholar]
  12. Dambeck, H.; Metzner, G.; Behn, U. Modeling Th1-Th2 regulation, allergy, and hyposensitization. In Dynamical Modeling in Biotechnology; Lectures Presented at the EU Advanced Workshop; World Scientific: Singapore, 2001; pp. 227–243. [Google Scholar]
  13. Fouchet, D.; Regoes, R. A Population Dynamics Analysis of the Interaction between Adaptive Regulatory T Cells and Antigen Presenting Cells. PLoS ONE 2008, 3, e2306. [Google Scholar] [CrossRef]
  14. Gross, F.; Metzner, G.; Behn, U. Mathematical modeling of allergy and specific immunotherapy: Th1, Th2 and Treg interactions. J. Theor. Biol. 2011, 269, 70–78. [Google Scholar] [CrossRef] [PubMed]
  15. Vogel, R.; Behn, U. Th1–Th2 Regulation and Allergy: Bifurcation Analysis of the Non-Autonomous System. Math. Model. Biol. Syst. 2008, 2, 145–155. [Google Scholar]
  16. Diehl, S.; Anguita, J.; Hoffmeyer, A.; Zapton, T.; Ihle, J.N.; Fikrig, E.; Rinco, M. Inhibition of Th1 Differentiation by IL-6 is Mediated by SOCS1. Immunity 2000, 13, 805–815. [Google Scholar] [CrossRef]
  17. Friedman, A.; Liao, K.L.; Bai, X.F. Mathematical modeling of interleukin-27 induction of anti-tumor t cells response. PLoS ONE 2014, 9, e91844. [Google Scholar]
  18. Wu, G. Calculation of steady-state distribution delay between central and peripheral compartments in two-compartment models with infusion regimen. Eur. J. Drug Metab. Pharmacokinet. 2002, 27, 259–264. [Google Scholar] [CrossRef]
  19. Kharitonov, V.L. Time-Delay Systems: Lyapunov Functionals and Matrices; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  20. Domoshnitsky, A.; Friedman, E. A positivity-based approach to delay-dependent stability of systems with large time-varying delays. Syst. Control Lett. 2016, 97, 139–148. [Google Scholar] [CrossRef]
  21. Rumyantsev, V.V.; Vorotnikov, V.I. Foundations of Partial Stability and Control; Publisher: Nijnii Tagil, Russia, 2014. (In Russian) [Google Scholar]
  22. Aleksandrov, A.; Aleksandrova, E.; Zhabko, A.; Chen, Y. Partial stability analysis of some classes of nonlinear systems. Acta Math. Sci. 2017, 37, 329–341. [Google Scholar] [CrossRef]
  23. Corduneanu, C. On partial stability for delay systems. Ann. Pol. Math. 1975, 29, 357–362. [Google Scholar] [CrossRef]
  24. Halanay, A. Differential Equation: Stability, Oscillations, Time-Lags; Academic Press: New York, NY, USA, 1966. [Google Scholar]
  25. Hatvani, L. On partial asymptotic stability and instability, I (Autonomous systems). Acta Sci. Math. 1983, 45, 219–231. [Google Scholar]
  26. Vorotnikov, V.I. Partial Stability and Control: The State-of-the-Art and Development Prospects. Ph.D. Thesis, Ural State Technical University, Nizhni Tagil, Russia, 2003. [Google Scholar]
  27. Kim, P.; Lee, P.; Levy, D. A theory of immunodominance and adaptive regulation. Bull. Math. Biol. 2010, 73, 1645–1665. [Google Scholar] [CrossRef] [PubMed]
  28. Nazari, F.; Pearson, A.T.; NoÈ, J.E.; Jackson, T.L. A mathematical model for IL-6-mediated, stem cell driven tumor growth and targeted treatmen. PLoS Comput. Biol. 2018, 14, e1005920. [Google Scholar] [CrossRef] [PubMed]
  29. Silva, J.; Morais, R.; Rodrigues da Silva, I.; Adimy, M.; Mancera, P. A Mathematical model for treatment of papillary thyroid cancer using the allee effect. J. Biol. Syst. 2020, 28, 701–718. [Google Scholar] [CrossRef]
Figure 1. For a small disturbance in the initial conditions near E 1 , the system exhibits stability. Here, E 1 = [ 0.6667 , 0 , 0 , 0 , 0 , 0 , 0 , 0.3750 , 1.8750 , 0.3367 ] .
Figure 1. For a small disturbance in the initial conditions near E 1 , the system exhibits stability. Here, E 1 = [ 0.6667 , 0 , 0 , 0 , 0 , 0 , 0 , 0.3750 , 1.8750 , 0.3367 ] .
Axioms 13 00102 g001
Figure 2. The stable variables in E 2 as proved in Theorem 3.
Figure 2. The stable variables in E 2 as proved in Theorem 3.
Axioms 13 00102 g002
Figure 3. The remaining variables in E 2 . E 2 = [ 0.6623 , 1.0549 × 10−4, 0 , 0 , 9.5901 × 10−6, 0 , 0 , 0.3750 , 1.8750 , 0.3361 ] .
Figure 3. The remaining variables in E 2 . E 2 = [ 0.6623 , 1.0549 × 10−4, 0 , 0 , 9.5901 × 10−6, 0 , 0 , 0.3750 , 1.8750 , 0.3361 ] .
Axioms 13 00102 g003
Table 1. Parameter values.
Table 1. Parameter values.
The production rate of naive cells [27] α 0.02
The proliferation rate of stimulated T cells [14]v1
Cytokines produced by the immune system [14]c 10 4
The strength of the suppression rate of Th1 by Th2 [14] μ 2 0.1
The strength of suppression of Th2 by Th1 [14] μ 1 0.2
The strength of the suppression rate by Treg [14] μ r 0.25
The differences in the autocrine action of the three subsets in Equation (3) [14] ϕ 0.05
The differences in the autocrine action of the three subsets in Equation (4) [14] κ 0.1
The death rate of naive APCs [27] γ 1 0.3
The death rate of mature APCs [27] γ 2 0.1
The death rate of naive T cells [27] γ 3 0.03
The natural decay of IL-6 [9] γ 4 0.4152
The inhibition rate of Th1 cells by IL-6 [28] η 1 0.60
The stimulation rate of Th2 cells by IL-6 [28] η 2 0.007
The inhibition rate of Treg cells by IL-6 [28] η r 0.39
The amount of the injected drug dose during desensitization [13] Λ 1
Delay to production of IL-6 [28] τ 1 0.25
The second time delay [18] τ 2 0.0794
The third time delay [11] τ 3 0.08
The inter compartment migration rate [14] θ 0.1
The production of IL-6 by other cells [29] k 1 0.055
The birth rate of naive APCs [27] λ 0.3
The rate of APC activation by the antigen [13] β 0.5
The rate of APC inhibition by regulatory T cells [13] μ 10 2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abdullah, R.; Badralexi, I.; Halanay, A. Stability Analysis in a Mathematical Model for Allergic Reactions. Axioms 2024, 13, 102. https://doi.org/10.3390/axioms13020102

AMA Style

Abdullah R, Badralexi I, Halanay A. Stability Analysis in a Mathematical Model for Allergic Reactions. Axioms. 2024; 13(2):102. https://doi.org/10.3390/axioms13020102

Chicago/Turabian Style

Abdullah, Rawan, Irina Badralexi, and Andrei Halanay. 2024. "Stability Analysis in a Mathematical Model for Allergic Reactions" Axioms 13, no. 2: 102. https://doi.org/10.3390/axioms13020102

APA Style

Abdullah, R., Badralexi, I., & Halanay, A. (2024). Stability Analysis in a Mathematical Model for Allergic Reactions. Axioms, 13(2), 102. https://doi.org/10.3390/axioms13020102

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop