Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Rango | (−1,1) | (−∞,2) | (−2,0) | (0,1) | 1 | (−1,0) | (−1.0) | (0,1) | 1 | (−1,0) | (−1,0) | (−1,0) |
This paper is concerned with a nonlinear plate equation modeling the oscillations of suspension bridges. Under mixed boundary conditions consisting of simply supported and free boundary conditions, we obtain the global well-posedness of solutions in suitable function spaces. In addition, we use the perturbed energy method to prove the existence of a bounded absorbing set and establish a stabilizability estimate. Then, we derive the existence of a global attractor by verifying the asymptotic smoothness of the corresponding dissipative dynamical system.
Citation: Yang Liu. Global attractors for a nonlinear plate equation modeling the oscillations of suspension bridges[J]. Communications in Analysis and Mechanics, 2023, 15(3): 436-456. doi: 10.3934/cam.2023021
[1] | Prashant Tripathi, Vivek Kumar Gupta, Anurag Dixit, Raghvendra Kumar Mishra, Satpal Sharma . Development and characterization of low cost jute, bagasse and glass fiber reinforced advanced hybrid epoxy composites. AIMS Materials Science, 2018, 5(2): 320-337. doi: 10.3934/matersci.2018.2.320 |
[2] | Albert Uchenna Ude, Che Husna Azhari . Lateral crashworthiness response of bombyx mori fibre/glass–fibre/epoxy hybrid composite cylindrical tubes-experimental. AIMS Materials Science, 2019, 6(6): 1227-1239. doi: 10.3934/matersci.2019.6.1227 |
[3] | Venkatarajan Subbarayalu, Subbu Chinnaraman, Athijayamani Ayyanar, Jayaseelan Chinnapalanisamy . Mechanical properties of vinyl ester hybrid composite laminates reinforced with screw pine and glass fiber. AIMS Materials Science, 2024, 11(1): 114-128. doi: 10.3934/matersci.2024007 |
[4] | Kiyotaka Obunai, Daisuke Mikami, Tadao Fukuta, Koichi Ozaki . Microstructure and mechanical properties of newly developed SiC-C/C composites under atmospheric conditions. AIMS Materials Science, 2018, 5(3): 494-507. doi: 10.3934/matersci.2018.3.494 |
[5] | Sujith Bobba, Z. Leman, E.S. Zainudin, S.M. Sapuan . Hoop tensile strength behaviour between different thicknesses E-glass and S-glass FRP rings. AIMS Materials Science, 2019, 6(3): 315-327. doi: 10.3934/matersci.2019.3.315 |
[6] | M. Rajanish, N. V. Nanjundaradhya, Ramesh S. Sharma, H. K. Shivananda, Alok Hegde . Directional Interlaminar Shear Strength (ILSS) of nano-modified epoxy/unidirectional glass fibre composite. AIMS Materials Science, 2018, 5(4): 603-613. doi: 10.3934/matersci.2018.4.603 |
[7] | Iketut Suarsana, Igpagus Suryawan, NPG Suardana, Suprapta Winaya, Rudy Soenoko, Budiarsa Suyasa, Wijaya Sunu, Made Rasta . Flexural strength of hybrid composite resin epoxy reinforced stinging nettle fiber with silane chemical treatment. AIMS Materials Science, 2021, 8(2): 185-199. doi: 10.3934/matersci.2021013 |
[8] | Kator Jeff Jomboh, Adele Dzikwi Garkida, Emmanuel Majiyebo Alemaka, Mohammed Kabir Yakubu, Vershima Cephas Alkali, Wilson Uzochukwu Eze, Nuhu Lawal . Properties and applications of natural, synthetic and hybrid fiber reinforced polymer composite: A review. AIMS Materials Science, 2024, 11(4): 774-801. doi: 10.3934/matersci.2024038 |
[9] | Albert Uchenna Ude, Che Husna Azhari . Crashworthiness response of natural silk-fibre glass hybrid reinforced epoxy cylindrical composite tubes under quasi-static load. AIMS Materials Science, 2019, 6(5): 852-863. doi: 10.3934/matersci.2019.5.852 |
[10] | C. N. Panagopoulos, E. P. Georgiou, D.A. Lagaris, V. Antonakaki . The effect of nanocrystalline Ni-W coating on the tensile properties of copper. AIMS Materials Science, 2016, 3(2): 324-338. doi: 10.3934/matersci.2016.2.324 |
This paper is concerned with a nonlinear plate equation modeling the oscillations of suspension bridges. Under mixed boundary conditions consisting of simply supported and free boundary conditions, we obtain the global well-posedness of solutions in suitable function spaces. In addition, we use the perturbed energy method to prove the existence of a bounded absorbing set and establish a stabilizability estimate. Then, we derive the existence of a global attractor by verifying the asymptotic smoothness of the corresponding dissipative dynamical system.
Tuberculosis (TB) is an infectious disease caused by the Mycobacterium tuberculosis (Mtb). It is both preventable and generally curable. However, in 2022, TB was the second leading cause of death worldwide due to a single infectious agent, following coronavirus disease (COVID-19). Current statistics are encouraging, showing a decline in the incidence, prevalence, and mortality rates over the years. Between 2015 and 2022, the incidence rate decreased by 8.9%, and TB-related deaths fell by 19% during the same period [1]. Despite these improvements, the situation remains concerning, with approximately a quarter of the world's population estimated to be infected with Mtb, and over 10 million people falling ill with TB each year [1].
After an individual is infected with Mtb, the risk of developing TB disease is highest in the first 2 years (approximately 5%), after which it decreases. Without treatment, the mortality rate from TB is high (approximately 50%). With the treatments currently recommended by WHO (a 4- to 6-month course of anti-TB drugs), approximately 85% of people with TB can be cured [1]. The above allows us to assume that one in two sick people not treated with any medication overcomes the disease thanks to their immune system (a complex network of cells, tissues, organs, and substances they produce, which help the body fight infections and other diseases).
Defense mechanisms against substances that are considered harmful or foreign are divided into innate immune response and acquired immune response. The role of adaptive or acquired immunity in controlling Mtb infection is well-studied, and the protective role of T and B lymphocytes in anti-tuberculosis immunity is widely recognized [2]. However, following exposure to Mtb, it typically takes 4 to 6 weeks for a human host, and 2 to 3 weeks for mice, to develop antigen-specific T-cell responses in the lymph nodes [3]. Therefore, the non-specific innate immune response plays a pivotal role in protecting the host before the onset of adaptive immunity and even leads to early clearance of bacteria [3].
Innate immunity constitutes the first defense against harmful substances or antigens that try to enter the body, which also includes barriers such as the skin, mucous membranes, tears, and stomach acid [2,4]. In TB, the innate immune response is crucial in determining the course of infection. The processes of transmission, infection, and development of TB are well-defined [3].
Transmission of Mtb begins when an infected individual expels microdroplets containing the infectious bacillus into the environment, which are then inhaled by an uninfected individual. Infectious bacilli can persist for up to 3 hours in closed, non-ventilated environments, and those that overcome the intrinsic barriers of the innate immune response (skin, mucous membranes, tears) enter the respiratory tract, which serves as a means of transport, generally to the pulmonary alveoli [3,5]. It is important to note that there are many cells of the innate immune system that participate in the response against TB; however, the precise mechanisms of TB immunopathology have not yet been fully understood [3].
Once Mtb overcomes the first barriers, it encounters the alveolar epithelial cells (AEC), which are the first to interact with any pathogen that attempts to enter the respiratory tract. AECs are also known to be vital for host defenses against Mtb, and although knowledge about their immunopathological functions in TB has increased considerably, they are still not fully elucidated [3,6].
In the alveoli, bacteria are detected, recognized, phagocytosed, and processed by antigen-presenting cells (APCs). The first line of defense against Mtb are alveolar macrophages (AM). Through the expression of several cell surface receptors, the macrophage recognizes, binds, and internalizes foreign particles, including Mtb. These cells phagocytose the Mtb by placing it inside a vacuole called a phagosome, which fuses with a bag of enzymes called a lysosome, giving rise to the phagolysosome [4].
Through a complex process, the enzymes of the phagolysosome control or eliminate the bacteria and also generate a signaling cascade that involves the production of cytokines and chemokines. This stimulate the activation of the antimicrobial activities of the phagocytes and recruitment of other cells of the immune system to the site of infection. However, bacilli that are not destroyed manipulate phagosome maturation to prevent phagolysosome formation and evade their removal. Furthermore, the phagosome binds to other vesicles with nutrients necessary for its replication inside the AM, converting the infected macrophage into a niche for bacterial replication [7].
In an attempt to deprive Mtb of its chosen cellular niche, macrophages infected with attenuated strains of Mtb can undergo apoptosis in a manner involving tumor necrosis factor-alpha, TNF−α, (proinflammatory cytokine secreted in the immune system by monocytes and macrophages, T and B lymphocytes, NK cells, and polymorphonuclear leukocytes). This cell death process also results in reduced survival of mycobacteria and leads to the priming of Mtb-specific CD8 T-cells [8,9].
Macrophage apoptosis may contribute to the innate immune response against this intracellular infection by containing and limiting the growth of bacilli, as has been observed for other infectious agents, and may therefore be part of a successful host defense mechanism [10]. On the other hand, TNF induces mitochondrial reactive oxygen species (ROS) to cause necrosis of Mtb-infected macrophages. Necrosis of infected macrophages constitutes a critical pathogenic event in TB by releasing mycobacteria into the extracellular environment, permissive for Mtb population growth [11].
Neutrophils play a specific role of great clinical importance in TB, although controversial studies have established a close correlation between the development of TB and tissue infiltration by neutrophils. However, the functions of circulating neutrophils in the development of tuberculosis is largely unknown [11]. Through animal models and in vitro experiments with human cells, it has been suggested that endothelial cells (EC) also participate in the immunopathology of TB. In vitro studies have shown that human ECs succumb to Mtb infection, and can help in the dissemination of bacteria through the bloodstream and lymph [3]. Also, ECs produce ROS, nitrogen intermediates (RNI), and antimicrobial peptides such as defensins to control bacterial growth [3]. Similarly, it has been verified that other cells that also participate in the innate immune response are eosinophils, dendritic cells, inflammatory monocytes, natural killer cells, and innate lymphoid cells, among others [3,4].
As observed, Mtb interacts with various innate immune cells, however the immunopathology arising from these interactions is not fully determined. In this regard, further research is needed to understand the immunopathology of TB. It is well documented that macrophages exert an important role in innate immune response against Mtb but they are also a niche for Mtb replication and propagation [6,7,12,13]. However, new studies have revealed that other cells such as dendritic cells, neutrophils and monocytes also play a key role in the immunopathology of TB [14,15].
Mathematical models had been used to gain insights in within-host dynamics of TB, mainly, related to adaptive immunity [16,17,18,19]. Currently, it is known that the innate immune response is critical to the outcome of Mtb infection, since it can control bacterial progression or substantially influence the subsequent adaptive response of the host. In the literature, researchers have analyzed several aspects of the innate immune response against Mycobacterium tuberculosis throughout of mathematical modeling. In [20], the authors formulated and analyzed a model to describe the dynamics of innate immunology of macrophages against Mycobacterium tuberculosis. In [21], the authors formulated a system of ordinary differential equations that describes the interaction dynamics between Mtb, iron, lipids, and nitric oxide in the early phase of macrophage infection. They incorporated the bactericidal property of nitric oxide and the cellular regulation of iron and lipids to analyze the role they play in the outcome of infection. In [22], the authors developed a spatiotemporal model of the initial and innate response to tuberculosis. The model consists of coupled reaction-diffusion-advection equations governing the dynamics of macrophages (resting and infected), bacteria (extracellular and intracellular), and a chemokine released by the bacteria, each of which affects final granuloma size. The model describes the migration of uninfected macrophages to the site of infection (pulmonary alveoli) and their subsequent phagocytosis of bacteria. It aims to capture clearance by innate immunity or disease progression through granuloma growth. In [23], the early immune response to Mtb infection in the lungs was modeled by means of coupled reaction-diffusion-transport equations with chemotaxis. The results and conclusions are similar to those presented in [22]. Other computational approaches have also addressed the role of early innate immunity in TB [21,24,25,26].
In this work, we are interested in evaluating and comparing the response of macrophages with respect to the response of other cells of the innate immune system against TB. To this end, we formulate and analyze, at theoretical and numerical levels, a non-linear system of ordinary differential equations that describes the interaction between macrophages, the other cells of the innate immune system, and Mtb. We believe this is the first model that incorporates the dynamics of all innate immune cells against Mtb.
In order to characterize some aspects of TB immunology at the cellular level, in the context of the innate immune response, and focus on the role played by macrophage responses in the outcome of infection, in this section, we formulate a model that describes the dynamics of interaction between innate immune cells and Mycobacterium tuberculosis at the site of infection.
Since there are different innate immune cells that participate in the response against Mtb [3], we will group all the innate cells that act against TB into a single population that we will denote by C. We also consider the population of macrophages, M, which is divided into uninfected macrophages, MU (macrophages that have not had contact with Mtb or that phagocytize Mtb and eliminate it efficiently), and infected macrophages, MI (macrophages that phagocytose the bacteria but do not have the capacity to eliminate it, and the bacteria also reproduce inside it). Similarly, we consider the population of bacteria, B, which will be divided into two subpopulations: the population of internal bacteria, BI (bacteria that have been phagocytosed by innate cells), and the population of external bacteria, BE (bacteria that have not been phagocytosed or that have gone from internal to external due to the necrotic death of their host cells).
The cell population of the innate immune system at time t is represented by C(t), which is composed of the subgroups of alveolar epithelial cells, endothelial cells, macrophages, neutrophils, eosinophils, monocytes, dendritic cells, natural killer cells, invariant natural killer T lymphocytes and innate lymphoid cells. All of these cells play a role in the immunopathology of TB that favors the host and/or the bacteria. Unfortunately, for most of them, the functions they perform and the mechanisms they carry out are still a matter of research; available results come from individual studies, and the hypotheses established based on them must be empirically validated [3,4,15].
We assume that the site of infection are the pulmonary alveoli and that infectious bacilli have been transported there. At this site, resident cells recruited from the bone marrow via the thymus or lymph nodes, among others, have initiated an innate immune response. As mentioned in the previous section, the set of cells that participate in this response is quite broad. However, at this stage of the infection, none of them show a specific response to the antigen [3,4].
We suppose that innate cells are recruited at a constant rate, ΛC, and die at per capita natural death rate, μC. Since cells produced by the thymus and lymph differentiate into different types of innate immune cells, and we want to compare the macrophages response with respect to the response of the rest of the innate cells, we assume that a portion of innate cells differentiate into uninfected macrophages at a rate α.
Since a portion of the bacteria phagocytosed by macrophages is eliminated by them, and the other turns them into their reproduction niche, we assume that uninfected macrophages phagocytose Mtb at a rate ε, of which a fraction ν becomes infected. Other innate cells participate directly or indirectly in the death of infected macrophages [3,4]; for this reason, we assume that the innate cells eliminates the population of infected macrophages at a rate λ. Necrosis is classified by NCCD-201 as accidental cell death, and Mtb-infected macrophage necrosis occurs when the Mtb kills infected macrophages [27,28]. This process is a key part of the Mtb infection cycle, allowing the bacteria to spread and evade the host's immune response. In TB, the internal bacteria, in addition to reproducing within the macrophage, induce its death by necrosis, releasing necrotic tissue, and a certain amount of bacteria into the environment [3,28]. As a consequence, we assume that the natural death of infected macrophages is by necrosis, which occurs at a per capita rate μMI(BI), which is a function that depends on the internal bacteria. Consequently, the number of infected macrophages that die by necrosis is μMI(BI)MI, and in turn, these necrotic macrophages release on average r bacilli into the extracellular environment. This implies that the amount of external bacteria released into the environment would be proportional to the number of macrophages that die by necrosis, rμMI(BI). To date, the most appropriate functional response for μMI(BI) has not been empirically determined. However, there is experimental evidence that the amount of foreign bacteria produced by macrophage necrosis is proportional to the number of infected macrophages [29]. For this reason, in this model, we will assume that infected macrophages die by necrosis at a constant rate μMI. In addition, bacteria phagocytosed by macrophages are eliminated at a rate, κ.
Since there are other innate cells such as dendritic cells, neutrophils, or monocytes that phagocytose and eliminate Mtb [3,4], we assume that extracellular bacteria are eliminated by these innate cells at a rate η. We also assume that internal bacteria replicate inside an infected macrophage at a rate β, following a saturation kinetics with carrying capacity of the infected macrophages denoted by KI. Finally, μMU, μBE, and μBI denote the per capita natural death rates (reciprocal of the half-life) of uninfected macrophages, and intracellular and extracellular bacteria, respectively.
From the above assumptions, we obtain the following system of ordinary differential equations
dCdt=ΛC−αC−μCCdMUdt=αC−ενBEMU−μMUMUdMIdt=ενBEMU−λMIC−μMIMIdBEdt=rμMIMI−εκBEMU−ηBEC−μBEBEdBIdt=β(1−BIKI)MI−μBIBI. | (2.1) |
The set of biological interest is given by
Ω={(CMUMIBEBI)∈R5+:0≤C≤ΛCα+μC,0≤MU+MI≤αΛCμM(α+μC),0≤BE+BI≤ˉαμB}, | (2.2) |
where μM=min{μMU,μMI}, μB=min{μBE,μBI}, and ˉα=α(rμMI+β)ΛCμM(α+μC). The following lemma ensures that system (2.1) makes biological sense; that is, all solutions starting at Ω remain there for all t≥0.
Lemma 2.1. The set Ω defined in (2.2) is positively invariant for solutions of the system (2.1).
Proof. Let us start by verifying the existence of the solution to the initial value problem defined by the system (2.1) and x0=(C(0),MU(0),MI(0),BE(0),BI(0))∈Ω. Since the vector field defined by the right-hand side of (2.1) is C1(R5+), then the fundamental theorem of existence and uniqueness guarantees the existence of a solution x(t)∈R5+ in the interval [0,σ] [30]. Furthermore, if the compact set ˉΩ∈R5+ satisfies {y∈R5:y=f(t)for somet∈[0,σ)}, then by the extension theorem σ=∞ (Corollary 2, page 91, [30]). Now, suppose that x0∈Ω, then we have
0≤C(0)≤ΛCα+μC,0≤MU(0)+MI(0)≤αΛCμM(α+μC),0≤BE(0)+BI(0)≤ˉαμB. |
On the other hand, the solution of the first equation of the system (2.1) is given by
C(t)=ΛCα+μC+(C(0)−ΛCα+μC)e−(α+μC)t, |
since 0≤C(0)≤ΛC/(α+μC), then 0≤C(t)≤ΛC/(α+μC) for all t≥0. Adding the second and third equations of the system (2.1), we have
dMUdt+dMIdt=αC−μMUMU−λMIC−μMIMI, |
since μM≤μMU and μM≤μMI, then from the previous equation we obtain the following inequality
dMUdt+dMIdt≤αΛCα+μC−μM(MU+MI)−λMIC≤αΛCα+μC−μM(MU+MI), |
or equivalently
d(MU+MI)dt+μM(MU+MI)≤αΛCα+μC. |
The solution of the above inequality satisfies
MU(t)+MI(t)≤αΛCμU(α+μC)+(MU(0)+MI(0)−αΛCμM(α+μC))e−μMt, |
since 0≤MU(0)+MI(0)≤αΛCμM(α+μC), then 0≤MU(t)+MI(t)≤αΛCμM(α+μC) for all t≥0. Adding the fourth and fifth equations of (2.1) we have
dBEdt+dBIdt=rμMIMI+β(1−BIKI)MI−εκBEMU−ηBEC−μBEBE−μBIBI, |
since μB≤μBE, μB≤μBI, MI<αΛC/μM(α+μC), then we obtain the following inequality
d(BE+BI)dt+μB(BE+BI)<ˉα, |
where ˉα=α(rμMI+β)ΛCμM(α+μC). The solution of the above inequality satisfies
BE(t)+BI(t)<ˉαμB+(BE(0)+BI(0)−ˉαμB)e−μBt. |
Therefore, if B(0)≤ˉα/μB, then B(t)≤ˉα/μB for all t≥0.
Taking the left-hand side of (2.1) equal to the zero vector gives the following system of algebraic equations
ΛC−αC−μCC=0αC−ενBEMU−μMUMU=0ενBEMU−λMIC−μMIMI=0rμMIMI−εκBEMU−ηBEC−μBEBE=0β(1−BIKI)MI−μBIBI=0. | (3.1) |
From the first equation of (3.1), we obtain C=C∗, where
C∗=ΛCα+μC. | (3.2) |
Substituting (3.2) in the second equation of (3.1), we obtain
MU=αC∗ενBE+μMU. | (3.3) |
From the third equation of (3.1), we obtain
MI=ενBEMUλC∗+μMI=ενBEλC∗+μMIαC∗ενBE+μMU. | (3.4) |
Substituting (3.2), (3.3), and (3.4) in the fourth equation of (3.1), we obtain
rμMIενBEλC∗+μMIαC∗ενBE+μMU−εκBEαC∗ενBE+μMU−ηBEC∗−μBEBE=0. |
The solutions of the above equation are given by BE=0 and BE=B∗E, where
B∗E=[καC∗ν(ηC∗+μBE)+μMUεν](R0−1), | (3.5) |
and
R0=rμMIεναC∗(λC∗+μMI)[εκαC∗+μMU(ηC∗+μBE)]. | (3.6) |
If BE=0, we obtain the free-infection equilibrium
E0=(C∗,αC∗μMU,0,0,0). | (3.7) |
Now, substituting (3.5) in the Eqs (3.3) and (3.4) we obtain
M∗U=αC∗ενB∗E+μMUM∗I=ενB∗EλC∗+μMIαC∗ενB∗E+μMU. | (3.8) |
Substituting (3.8) in the fifth equation of (3.1), we obtain BI=B∗I, where
B∗I=βM∗IβKIM∗I+μBI. | (3.9) |
As a consequence, if R0>1, there exists a non-trivial equilibrium
E1=(C∗,M∗U,M∗I,B∗E,B∗I). | (3.10) |
The following theorem summarizes the results of the existence of equilibria.
Teorem 3.1. If R0≤1, then E0 defined in (3.7) is the only equilibrium in Ω. If R0>1, in addition to E0, there exists the infected equilibrium E1 defined in (3.10).
In this section, we will analyze the stability of equilibrium solutions. The linearization of the system (2.1) around an equilibrium solution ˉx is given by x=J(ˉx)x, where the Jacobian matrix J evaluated at x is given by
J(x)=(−(α+μC)0000α−(ενBE+μMU)0−ενMU0−λMIενBE−(λC+μMI)ενMU0−ηBE−εκBErμMI−(εκMU+ηC+μBE)000β(1−BIKI)0−(βMIKI+μBI)). | (3.11) |
From (3.11), we verify that the Jacobian matrix evaluated in trivial equilibrium E0 is given by
J(E0)=(−(α+μC)0000α−μMU0−εναC∗μMU000−(λC∗+μMI)εναC∗μMU000rμMI−(εκαC∗μMU+ηC∗+μBE)000β0−μBI.) | (3.12) |
The eigenvalues of J(E0) defined in (3.12) are given by −(α+μC), −μMU, −μBI, and the roots of the following quadratic equation
ς2+(λC∗+μMI+εκαC∗μMU+ηC∗+μBE)ς+1μMU(λC∗+μMI)(εκαC∗+μMU(ηC∗+μBE))(1−R0)=0. |
The above implies that E0 is locally asymptotically stable in Ω when R0<1. Similarly, it is verified that the eigenvalues of J(E1) are −(α+μC), −(βM∗I/KI+μBI) and the roots of the following polynomial
ς3+a2ς2+a1ς+a0, | (3.13) |
where
a2=εκM∗U+ηC∗+μBE+λC∗+μMI+ενB∗E+μMUa1=(εκM∗U+ηC∗+μBE+λC∗+μMI)(ενB∗E+μMU)−ενM∗UεκB∗Ea0=ενM∗U[rμMIενBE−εκB∗E(λC∗+μMI)]=ενM∗U(λC∗+μMI)(ηC∗+μBR)B∗EMU. |
Since the coefficients of the cubic polynomial (3.13) do not change sign, then from the rule of the signs of Descartes, it is verified that all its roots have a negative real part. Therefore, the equilibrium E1 is locally asymptotically stable in Ω. The following theorem summarizes the results of the existence of equilibria.
Teorem 3.2. If R0<1, then E0 defined in (3.7) is locally asymptotically stable in Ω. If R0>1, E0 is unstable, and E1 defined in (3.10) is locally asymptotically stable in Ω.
In the previous section, it was verified that the result of the infection depends on the basic radius of offspring R0 defined in (3.6). We observe that R0 depends on almost all the parameters of the model. In this sense, an interesting question is which parameters of the model (2.1) influence the most values of R0. To answer this question, it is necessary to perform a sensitivity analysis that determines the uncertainty of the model outputs based on the uncertainty of its inputs.
There are various theoretical and numerical methods to perform this type of analysis. In our case, we have an initial value problem defined by a system of nonlinear ordinary differential equations whose inputs are the model parameters and the initial conditions; and also, the behavior of the model outputs is conditioned by the value of R0. As such, good method to perform the analysis consists of using the normalized sensitivity index to determine the sensitivity of R0 with respect to each of the parameters that define it.
Let p be a parameter, and the normalized sensitivity index ΥR0p is given by the following partial derivative
ΥR0p=pR0∂R0∂p. | (4.1) |
From (4.1), we verify that
ΥR0ΛC=μMUμBEμMUμBE+(εκα+μMUη)C∗−λC∗λC∗+μMIΥR0α=1+αα+μC[λC∗λC∗+μMI−εκ(α+μC)C∗+μMUμBEμMUμBE+(εκα+μMUη)C∗]ΥR0μC=−μCα+μC[μMUμBEμMUμBE+(εκα+μMUη)C∗+λC∗λC∗+μMI]ΥR0ε=1−εκαC∗μMUμBE+(εκα+μMUη)C∗ΥR0ν=1ΥR0μMU=−μMUηC∗+μMUμBEμMUμBE+(εκα+μMUη)C∗ΥR0λ=−λC∗λC∗+μMIΥR0μMI=1−μMIλC∗+μMIΥR0r=1ΥR0κ=−εκαC∗μMUμBE+(εκα+μMUη)C∗ΥR0η=−μMUηC∗μMUμBE+(εκα+μMUη)C∗ΥR0μBE=−μMUμBEμMUμBE+(εκα+μMUη)C∗. | (4.2) |
From (4.2) we verify that sensitivity indices are within the ranges presented in Table 1.
Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Rango | (−1,1) | (−∞,2) | (−2,0) | (0,1) | 1 | (−1,0) | (−1.0) | (0,1) | 1 | (−1,0) | (−1,0) | (−1,0) |
From Table 1 we observe that ΥR0ΛC∈(−1,1), ΥR0α∈(−∞,2), ΥR0μC∈(−2,0), ΥR0ν=ΥR0r=1, {ΥR0μMU,ΥR0λ,ΥR0κ,ΥR0η,ΥR0μBE}∈(−1,0) and {ΥR0ε,ΥR0μMI}∈(0,1). The above indicates that the parameters μC, μMU, λ, κ, η, and μBE have a positive impact on infection control, in the sense that if the value of any of these parameters increases, then the value of R0 decreases. On the other hand, an increase in any of the parameters ν,r,μMI o ε leads to an increase in the value of R0.
Due to the form of the expressions that define ΥR0ΛC and ΥR0α, it was not possible to determine a smaller range; however, we can conclude that the parameters ΛC, α, and μC can have positive or negative effects on infection control. Note that, if we do not consider the three previous parameters, the parameters that most influence the outcome of the infection are r and ν.
Let us observe that the rate of differentiation of innate cells into macrophages, α, is the parameter that most influences the outcome of the infection (an increase of 5% in this parameter would generate an increase in R0 between 5.6% and 7.85%), followed by the parameters ν, r, ΛC, and μMI. It is also observed that when the phagocytic rate of macrophages has a high impact on the variation of R0, the rate of bacterial elimination by macrophages has a low impact, and vice versa. The incidence of the parameters μC and μBE on the variation of R0 is almost zero. Finally, the minimum value of μMU has a greater impact on the result of the infection than its maximum value.
Figures 1 and 2 were made with data from Table 2 and show graphs of the temporal behavior of the cell population. In Figure 1, it can be seen that in the first 10 hours, a large number of innate cells differentiate into macrophages, almost quintupling their initial population, all with the purpose of counteracting the spread of the infection. However, their response was insufficient because both the population of bacteria and infected macrophages increased considerably.
Parameter | Description | Values | References |
ΛC | Recruitment rate of innate cells | 50–120 cell/h | [31] |
α | Differentiation rate of macrophages | 0.02–0.06 h−1 | [32] |
ε | Phagocytosis rate of macrophages | 0.1–0.8 | [31] |
ν | Macrophage infection rate | 0.2–0.8 | [31] |
εν | Effective macrophage infection rate | 0.02–0.64 (h⋅ cell)−1 | [31] |
r | Nnumber of bact. released by necrotic MI | 1–50 cell | [31] |
κ | Bacterial elimination rate by macrophages | 0.001–0.01 | [31] |
εκ | Effective bacterial killing rate by macrophages | 0.0001–0.008 (h⋅ cell)−1 | [31] |
η | Bacterial elimination rate by innate cells | 0.003–0.03 (h⋅ cell)−1 | [31] |
β | BI replication rate within MI | 0.0032–0.032 h−1 | [31] |
λ | MI elimination rate by innate cells | 0.001–0.01 (h⋅ cell)−1 | [31] |
μC | Natural death rate of C | 0.0067–0.0083 h−1 | [31] |
μMU | Natural death rate of MU | 0.0032–0.0051 h−1 | [32] |
μMI | MI death rate due to necrosis | 0.0101–0.0139 h−1 | [31] |
μBE | Natural death rate of BE | 0.002–0.02 h−1 | [13] |
μBI | Natural death rate of BI | 0.002–0.03 h−1 | [13] |
KI | Carrying capacity of macrophages | 10–100 cell | [31] |
In contrast, in Figure 2, it can be seen that the initial population of macrophages declines very rapidly, and, although at the beginning the populations of external bacteria and infected macrophages grow rapidly, after 2 hours these decrease to less than 10 bacteria and 10 infected macrophages. Table 4 shows the values of sensitivity indices for the simulations presented in Figures 1 and 2.
Index | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Minimum | −0.99 | 1.57 | −0.24 | 0.83 | 1 | −0.83 | −0.99 | 0.99 | 1 | −0.17 | −0.83 | −2.95×10−4 |
Maximum | −0.99 | 1.12 | −0.12 | 0.24 | 1 | −0.24 | −0.99 | 0.99 | 1 | −0.76 | −0.24 | −9.17×10−5 |
Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Figure 1 | −0.89 | 1.74 | −0.11 | 0.95 | 1 | −0,95 | −0.90 | 0.90 | 1 | −0.05 | −0.94 | −7×10−3 |
Figure 2 | −0.99 | 0.89 | −0.12 | 0.01 | 1 | −0,01 | −0.99 | 0.99 | 1 | −0.99 | −0.01 | −1.2×10−4 |
It is important to highlight that although in Figure 1 it can be seen that at the beginning the population of uninfected macrophages grows, and in Figure 2 this population decreases from the beginning, the sensitivity analysis reveals that the innate immune response of uninfected macrophages was better in the second simulation. Indeed, the magnitude of ΥR0κ was larger in the second simulation, indicating that more bacteria were eliminated in the second simulation than in the first one.
In a similar way, indicating that more uninfected macrophages were killed in the first simulation than in the second one (Figure 2), and decreased steadily in the first simulation (Figure 1). Since the magnitude of ΥR0η was larger in the first simulation and the magnitude of ΥR0α was smaller in the second simulation, parameter analysis suggests that its incidence was significantly higher in the first one.
The infection and replication rates have a high impact on the outcome of the infection, but they have the same incidence in both cases. Although the recruitment of bone marrow cells is low in the first simulation and high in the second one, the level of incidence on the outcome of the infection is similar in both cases. The incidence level of both the death rate of innate cells that do not differentiate into infected macrophages as well as the death rate of external cells did not have a significant impact on bacterial progression.
Microbiological processes at any level, and in particular at the cellular level, are extremely complex and require contributions from different areas of knowledge to advance in their understanding. Our model is a very simple approximation to the interaction dynamics of innate cells (divided between macrophages and the rest of the innate cells without specifying) and Mycobacterium tuberculosis, at the beginning of the Mtb infection.
Qualitative analysis reveals a forward bifurcation in which the infection-free equilibrium E0 loses its stability when R0=1 and emerges a stable coexistence equilibrium E1. Note that the progression of the infection will depend on a condition very similar to that used in epidemiology for the spread of an infectious disease in a population of host individuals. The key quantity is the basic offspring proportion R0 defined in (3.6), which is interpreted as the number of secondary infections arising from a macrophage during its lifetime when all other innate cells are not infected. The fact that R0 depends on all the parameters of the model except for three parameters associated with internal bacteria dynamics suggests that both macrophages and the rest of innate cells must perform their effector functions appropriately to counteract Mtb pathogenesis. It was verified that the parameters associated with macrophages immunopathology (ε, ν, r, κ, μMU, μMI and α) significantly impact the interaction dynamics. The results above corroborates that macrophages play a fundamental role in the control or dissemination of the bacillus, during the first stage of infection. In addition, the analysis provides new insights that suggest focusing our attention on the rest of the innate cells that participate in the activation of the innate immune response against TB, since the parameters associated with these cells ΛC, α, μC, λ, and η are as relevant as those of the macrophages in the variation of R0.
The parameter α represents the differentiation rate of cells produced by the bone marrow into mature macrophages capable of performing their effector functions against any pathogen, in particular against Mtb. In reality, this parameter is an oversimplification homeostasis, a process with many unknown aspects at the biological level, but that begins with the production of stem cells in the bone marrow and involves the transition to other cells (such as the differentiation of promonocytic cells into mature monocytes); this process in turn includes self-regulation mechanisms, among other factors. Note that this parameter is associated with responses from macrophages and other innate cells. With the largest range among all the sensitivity indices (−∞,2), the results suggest that it can contribute to either reducing or increasing bacterial progression. Indeed, if ΥR0α∈(−∞,0) contributes to reducing the value of R0, which will decrease the number of secondary infections. On the other hand, if ΥR0α∈(0,2) increases secondary infections. Under the set of parameters used to perform the simulations in Figures 1 and 2, positive sensitivity indices were obtained for α. As a consequence, for both cases, the differentiation rate contributed to the progression of the infection. By carrying out a similar procedure, it is verified that the recruitment rate of innate cells ΛC can contribute to the decrease or increase of infection. In fact, in the numerical simulations of Figures 1 and 2, it positively affects the control of infection.
At the experimental level, there is a large body of literature supporting the theoretical results obtained here, which highlight that Mtb infection can be eliminated by the innate immune system before an adaptive immune response is initiated; this innate protection requires a variety of robust cell-autonomous responses from many different types of host immune cells [6].
For centuries, TB research has focused on the role of adaptive immunity in TB, leaving aside innate immunity. However, advances in recent decades have shown that the innate immune response determines the outcome of the infection. In this sense, scientific and academic production on the innate immunology of TB has increased, and generated new perspectives. However, there are still many gaps. In this article, we have formulated and analyzed a model whose results corroborate the importance of the innate response of macrophages in TB, and also suggest that other innate cells are also very important during primary infection.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
We want to thank anonymous referees for their valuable comments that helped us to improve the paper. Eduardo Ibarguen-Mondragón and Miller Cerón Gómez acknowledges the partial support from Universidad de Nariño. M. Victoria Otero-Espinar is partially supported by the Ministerio de Ciencia e Innovación, Agencia Estatal de Investigación (Spain), Grant PID2020-115155GB-I00 and the Consellería de Educación, Universidade e Formación Profesional (Xunta de Galicia), Grant ED431C 2023/31 with FEDER funds.
The authors declare there is no conflict of interest.
[1] | R. A. Adams, Sobolev Spaces, Academic Press, New York, 1975. |
[2] |
M. Al-Gwaiz, V. Benci, F. Gazzola, Bending and stretching energies in a rectangular plate modeling suspension bridges, Nonlinear Anal., 106 (2014), 18–34. https://doi.org/10.1016/j.na.2014.04.011 doi: 10.1016/j.na.2014.04.011
![]() |
[3] |
G. Arioli, F. Gazzola, A new mathematical explanation of what triggered the catastrophic torsional mode of the Tacoma Narrows Bridge, Appl. Math. Model., 39 (2015), 901–912. https://doi.org/10.1016/j.apm.2014.06.022 doi: 10.1016/j.apm.2014.06.022
![]() |
[4] |
U. Battisti, E. Berchio, A. Ferrero, F. Gazzola, Energy transfer between modes in a nonlinear beam equation, J. Math. Pures Appl., 108 (2017), 885–917. https://doi.org/10.1016/j.matpur.2017.05.010 doi: 10.1016/j.matpur.2017.05.010
![]() |
[5] |
E. Berchio, A. Ferrero, F. Gazzola, Structural instability of nonlinear plates modelling suspension bridges: mathematical answers to some long-standing questions, Nonlinear Anal. Real World Appl., 28 (2016), 91–125. https://doi.org/10.1016/j.nonrwa.2015.09.005 doi: 10.1016/j.nonrwa.2015.09.005
![]() |
[6] |
D. Bonheure, F. Gazzola, E. M. Dos Santos, Periodic solutions and torsional instability in a nonlinear nonlocal plate equation, SIAM J. Math. Anal., 51 (2019), 3052–3091. https://doi.org/10.1137/18M1221242 doi: 10.1137/18M1221242
![]() |
[7] | I. Chueshov, I. Lasiecka, Long-Time Behavior of Second Order Evolution Equations with Nonlinear Damping, Mem. Amer. Math. Soc., Providence, RI, 2008. https://doi.org/10.1090/memo/0912 |
[8] |
V. Ferreira Jr, F. Gazzola, E. M. dos Santos, Instability of modes in a partially hinged rectangular plate, J. Differential Equations, 261 (2016), 6302–6340. https://doi.org/10.1016/j.jde.2016.08.037 doi: 10.1016/j.jde.2016.08.037
![]() |
[9] |
A. Ferrero, F. Gazzola, A partially hinged rectangular plate as a model for suspension bridges, Discrete Contin. Dyn. Syst., 35 (2015), 5879–5908. https://doi.org/10.3934/dcds.2015.35.5879 doi: 10.3934/dcds.2015.35.5879
![]() |
[10] | F. Gazzola, Mathematical Models for Suspension Bridges: Nonlinear Structural Instability, in: MS & A. Modeling, Simulation and Applications, vol. 15, Springer, Cham, 2015. https://doi.org/10.1007/978-3-319-15434-3 |
[11] |
F. Gazzola, M. Squassina, Global solutions and finite time blow up for damped semilinear wave equations, Ann. Inst. H. Poincaré Anal. Non Linéaire, 23 (2006), 185–207. https://doi.org/10.1016/j.anihpc.2005.02.007 doi: 10.1016/j.anihpc.2005.02.007
![]() |
[12] |
J. B. Han, R. Z. Xu, C. Yang, Continuous dependence on initial data and high energy blowup time estimate for porous elastic system, Commun. Anal. Mech., 15 (2023), 214–244. https://doi.org/10.3934/cam.2023012 doi: 10.3934/cam.2023012
![]() |
[13] |
A. Haraux, E. Zuazua, Decay estimates for some semilinear damped hyperbolic problems, Arch. Rational Mech. Anal., 100 (1988), 191–206. https://doi.org/10.1007/BF00282203 doi: 10.1007/BF00282203
![]() |
[14] |
A. C. Lazer, P. J. McKenna, Large scale oscillatory behaviour in loaded asymmetric systems, Ann. Inst. H. Poincaré Anal. Non Linéaire, 4 (1987), 243–274. https://doi.org/10.1016/S0294-1449(16)30368-7 doi: 10.1016/S0294-1449(16)30368-7
![]() |
[15] |
A. C. Lazer, P. J. McKenna, Large-amplitude periodic oscillations in suspension bridges: some new connections with nonlinear analysis, SIAM Rev., 32 (1990), 537–578. https://doi.org/10.1137/1032120 doi: 10.1137/1032120
![]() |
[16] |
W. Lian, R. Z. Xu, Global well-posedness of nonlinear wave equation with weak and strong damping terms and logarithmic source term, Adv. Nonlinear Anal., 9 (2020), 613–632. https://doi.org/10.1515/anona-2020-0016 doi: 10.1515/anona-2020-0016
![]() |
[17] |
Y. Liu, B. Moon, V. D. Rǎdulescu, R. Z. Xu, C. Yang, Qualitative properties of solution to a viscoelastic Kirchhoff-like plate equation, J. Math. Phys., 64 (2023), 051511. https://doi.org/10.1063/5.0149240 doi: 10.1063/5.0149240
![]() |
[18] |
Y. Liu, J. Mu, Y. J. Jiao, A class of fourth order damped wave equations with arbitrary positive initial energy, Proc. Edinburgh Math. Soc., 62 (2019), 165–178. https://doi.org/10.1017/S0013091518000330 doi: 10.1017/S0013091518000330
![]() |
[19] |
Y. B. Luo, R. Z. Xu, C. Yang, Global well-posedness for a class of semilinear hyperbolic equations with singular potentials on manifolds with conical singularities, Calc. Var. Partial Differential Equations, 61 (2022), 210. https://doi.org/10.1007/s00526-022-02316-2 doi: 10.1007/s00526-022-02316-2
![]() |
[20] |
Q. F. Ma, S. H. Wang, C. K. Zhong, Necessary and sufficient conditions for the existence of global attractors for semigroups and applications, Indiana Univ. Math. J., 51 (2002), 1541–1559. https://doi.org/10.1512/iumj.2002.51.2255 doi: 10.1512/iumj.2002.51.2255
![]() |
[21] |
P. J. McKenna, W. Walter, Nonlinear oscillations in a suspension bridge, Arch. Ration. Mech. Anal., 98 (1987), 167–177. https://doi.org/10.1007/BF00251232 doi: 10.1007/BF00251232
![]() |
[22] |
P. J. McKenna, W. Walter, Travelling waves in a suspension bridge, SIAM J. Appl. Math., 50 (1990), 703–715. https://doi.org/10.1137/0150041 doi: 10.1137/0150041
![]() |
[23] |
J. Y. Park, J. R. Kang, Global attractors for the suspension bridge equations with nonlinear damping, Quart. Appl. Math., 69 (2011), 465–475. https://doi.org/10.1090/S0033-569X-2011-01259-1 doi: 10.1090/S0033-569X-2011-01259-1
![]() |
[24] | R. Scott, In the Wake of Tacoma: Suspension Bridges and the Quest for Aerodynamic Stability, ASCE Press, Reston, VA, 2001. https://doi.org/10.1061/9780784405420 |
[25] |
J. Simon, Compact sets in the space Lp(0,T;B), Ann. Mat. Pura Appl., 146 (1986), 65–96. https://doi.org/10.1007/BF01762360 doi: 10.1007/BF01762360
![]() |
[26] |
W. A. Strauss, On continuity of functions with values in various Banach spaces, Pacific J. Math., 19 (1966), 543–551. https://doi.org/10.2140/pjm.1966.19.543 doi: 10.2140/pjm.1966.19.543
![]() |
[27] | E. Ventsel, T. Krauthammer, Thin Plates and Shells: Theory, Analysis, and Applications, Marcel Dekker, New York, 2001. https://doi.org/10.1115/1.1483356 |
[28] |
Y. D. Wang, Finite time blow-up and global solutions for fourth order damped wave equations, J. Math. Anal. Appl., 418 (2014), 713–733. https://doi.org/10.1016/j.jmaa.2014.04.015 doi: 10.1016/j.jmaa.2014.04.015
![]() |
[29] |
X. C. Wang, R. Z. Xu, Global existence and finite time blowup for a nonlocal semilinear pseudoparabolic equation, Adv. Nonlinear Anal., 10 (2021), 261–288. https://doi.org/10.1515/anona-2020-0141 doi: 10.1515/anona-2020-0141
![]() |
[30] |
H. Y. Xu, Existence and blow-up of solutions for finitely degenerate semilinear parabolic equations with singular potentials, Commun. Anal. Mech., 15 (2023), 132–161. https://doi.org/10.3934/cam.2023008 doi: 10.3934/cam.2023008
![]() |
[31] |
R. Z. Xu, J. Su, Global existence and finite time blow-up for a class of semilinear pseudo-parabolic equations, J. Funct. Anal., 264 (2013), 2732–2763. https://doi.org/10.1016/j.jfa.2013.03.010 doi: 10.1016/j.jfa.2013.03.010
![]() |
[32] |
R. Z. Xu, X. C. Wang, Y. B. Yang, S. H. Chen, Global solutions and finite time blow-up for fourth order nonlinear damped wave equation, J. Math. Phys., 59 (2018), 061503. https://doi.org/10.1063/1.5006728 doi: 10.1063/1.5006728
![]() |
[33] |
R. Z. Xu, M. Y. Zhang, S. H. Chen, Y. B. Yang, J. H. Shen, The initial-boundary value problems for a class of six order nonlinear wave equation, Discrete Contin. Dyn. Syst., 37 (2017), 5631–5649. https://doi.org/10.3934/dcds.2017244 doi: 10.3934/dcds.2017244
![]() |
[34] |
Y. Yang, Z. B. Fang, On a strongly damped semilinear wave equation with time-varying source and singular dissipation, Adv. Nonlinear Anal., 12 (2022), 20220267. https://doi.org/10.1515/anona-2022-0267 doi: 10.1515/anona-2022-0267
![]() |
[35] |
W. H. Yang, J. Zhou, Global attractors of the degenerate fractional Kirchhoff wave equation with structural damping or strong damping, Adv. Nonlinear Anal., 11 (2022), 993–1029. https://doi.org/10.1515/anona-2022-0226 doi: 10.1515/anona-2022-0226
![]() |
[36] |
C. K. Zhong, Q. Z. Ma, C. Y. Sun, Existence of strong solutions and global attractors for the suspension bridge equations, Nonlinear Anal., 67 (2007), 442–454. https://doi.org/10.1016/j.na.2006.05.018 doi: 10.1016/j.na.2006.05.018
![]() |
1. | Mohammed Y. Abdellah, Mohamed K. Hassan, Ahmed F. Mohamed, Khalil Abdelrazek Khalil, A Novel and Highly Effective Natural Vibration Modal Analysis to Predict Nominal Strength of Open Hole Glass Fiber Reinforced Polymer Composites Structure, 2021, 13, 2073-4360, 1251, 10.3390/polym13081251 | |
2. | Rebai Billel, Contribution to study the effect of (Reuss, LRVE, Tamura) models on the axial and shear stress of sandwich FGM plate (Ti–6A1–4V/ZrO2) subjected on linear and nonlinear thermal loads, 2023, 10, 2372-0484, 26, 10.3934/matersci.2023002 | |
3. | Balaji Rajendran, K. Krishna Prabakar, Investigation on the damage mechanisms using acoustic emission method and damping characteristics of hybrid flax‐glass composites, 2024, 45, 0272-8397, 14521, 10.1002/pc.28779 | |
4. | Mario Ceddia, Giuseppe Solarino, Maria Tucci, Luciano Lamberti, Bartolomeo Trentadue, Stress Analysis of Tibial Bone Using Three Different Materials for Bone Fixation Plates, 2024, 8, 2504-477X, 334, 10.3390/jcs8090334 |
Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Rango | (−1,1) | (−∞,2) | (−2,0) | (0,1) | 1 | (−1,0) | (−1.0) | (0,1) | 1 | (−1,0) | (−1,0) | (−1,0) |
Parameter | Description | Values | References |
ΛC | Recruitment rate of innate cells | 50–120 cell/h | [31] |
α | Differentiation rate of macrophages | 0.02–0.06 h−1 | [32] |
ε | Phagocytosis rate of macrophages | 0.1–0.8 | [31] |
ν | Macrophage infection rate | 0.2–0.8 | [31] |
εν | Effective macrophage infection rate | 0.02–0.64 (h⋅ cell)−1 | [31] |
r | Nnumber of bact. released by necrotic MI | 1–50 cell | [31] |
κ | Bacterial elimination rate by macrophages | 0.001–0.01 | [31] |
εκ | Effective bacterial killing rate by macrophages | 0.0001–0.008 (h⋅ cell)−1 | [31] |
η | Bacterial elimination rate by innate cells | 0.003–0.03 (h⋅ cell)−1 | [31] |
β | BI replication rate within MI | 0.0032–0.032 h−1 | [31] |
λ | MI elimination rate by innate cells | 0.001–0.01 (h⋅ cell)−1 | [31] |
μC | Natural death rate of C | 0.0067–0.0083 h−1 | [31] |
μMU | Natural death rate of MU | 0.0032–0.0051 h−1 | [32] |
μMI | MI death rate due to necrosis | 0.0101–0.0139 h−1 | [31] |
μBE | Natural death rate of BE | 0.002–0.02 h−1 | [13] |
μBI | Natural death rate of BI | 0.002–0.03 h−1 | [13] |
KI | Carrying capacity of macrophages | 10–100 cell | [31] |
Index | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Minimum | −0.99 | 1.57 | −0.24 | 0.83 | 1 | −0.83 | −0.99 | 0.99 | 1 | −0.17 | −0.83 | −2.95×10−4 |
Maximum | −0.99 | 1.12 | −0.12 | 0.24 | 1 | −0.24 | −0.99 | 0.99 | 1 | −0.76 | −0.24 | −9.17×10−5 |
Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Rango | (−1,1) | (−∞,2) | (−2,0) | (0,1) | 1 | (−1,0) | (−1.0) | (0,1) | 1 | (−1,0) | (−1,0) | (−1,0) |
Parameter | Description | Values | References |
ΛC | Recruitment rate of innate cells | 50–120 cell/h | [31] |
α | Differentiation rate of macrophages | 0.02–0.06 h−1 | [32] |
ε | Phagocytosis rate of macrophages | 0.1–0.8 | [31] |
ν | Macrophage infection rate | 0.2–0.8 | [31] |
εν | Effective macrophage infection rate | 0.02–0.64 (h⋅ cell)−1 | [31] |
r | Nnumber of bact. released by necrotic MI | 1–50 cell | [31] |
κ | Bacterial elimination rate by macrophages | 0.001–0.01 | [31] |
εκ | Effective bacterial killing rate by macrophages | 0.0001–0.008 (h⋅ cell)−1 | [31] |
η | Bacterial elimination rate by innate cells | 0.003–0.03 (h⋅ cell)−1 | [31] |
β | BI replication rate within MI | 0.0032–0.032 h−1 | [31] |
λ | MI elimination rate by innate cells | 0.001–0.01 (h⋅ cell)−1 | [31] |
μC | Natural death rate of C | 0.0067–0.0083 h−1 | [31] |
μMU | Natural death rate of MU | 0.0032–0.0051 h−1 | [32] |
μMI | MI death rate due to necrosis | 0.0101–0.0139 h−1 | [31] |
μBE | Natural death rate of BE | 0.002–0.02 h−1 | [13] |
μBI | Natural death rate of BI | 0.002–0.03 h−1 | [13] |
KI | Carrying capacity of macrophages | 10–100 cell | [31] |
Index | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Minimum | −0.99 | 1.57 | −0.24 | 0.83 | 1 | −0.83 | −0.99 | 0.99 | 1 | −0.17 | −0.83 | −2.95×10−4 |
Maximum | −0.99 | 1.12 | −0.12 | 0.24 | 1 | −0.24 | −0.99 | 0.99 | 1 | −0.76 | −0.24 | −9.17×10−5 |
Índice | ΥR0ΛC | ΥR0α | ΥR0μC | ΥR0ε | ΥR0ν | ΥR0μMU | ΥR0λ | ΥR0μMI | ΥR0r | ΥR0κ | ΥR0η | ΥR0μBE |
Figure 1 | −0.89 | 1.74 | −0.11 | 0.95 | 1 | −0,95 | −0.90 | 0.90 | 1 | −0.05 | −0.94 | −7×10−3 |
Figure 2 | −0.99 | 0.89 | −0.12 | 0.01 | 1 | −0,01 | −0.99 | 0.99 | 1 | −0.99 | −0.01 | −1.2×10−4 |