Next Article in Journal
A Generalization of Szász–Mirakyan Operators Based on α Non-Negative Parameter
Next Article in Special Issue
Stability and Sensitivity Analysis of the COVID-19 Spread with Comorbid Diseases
Previous Article in Journal
A Constructive Method of Solving Local Boundary Value Problems for Nonlinear Systems with Perturbations and Control Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Study of Three Mathematical Models for the Interaction between the Human Immune System and a Virus

Department of Applied Mathematics, University of Craiova, Al. I. Cuza, 13, 200585 Craiova, Romania
Submission received: 10 July 2022 / Revised: 28 July 2022 / Accepted: 31 July 2022 / Published: 3 August 2022
(This article belongs to the Special Issue Mathematical Modeling of the Infectious Diseases and Their Controls)

Abstract

:
In this paper, we will consider three deterministic models for the study of the interaction between the human immune system and a virus: the logistic model, the Gompertz model, and the generalized logistic model (or Richards model). A qualitative analysis of these three models based on dynamical systems theory will be performed by studying the local behavior of the equilibrium points and obtaining the local dynamics properties from the linear stability point of view. Additionally, we will compare these models in order to understand which is more appropriate to model the interaction between the human immune system and a virus. Some natural medical interpretations will be obtained, which are available for all three models and can be useful to the medical community.

1. Introduction

According to the basic immunology treaties around the medical scientific world, from the structural point of view, the human immune system is formed by organs, cells, and molecules [1]. The organs of the human immune system are tonsils and adenoids, the thymus, lymph nodes, spleen, payer’s patches, appendix, lymphatic vessels, and bone marrow. The cells of the immune system are lymphocytes (T-lymphocytes, B-lymphocytes, plasma cells, and natural killer lymphocytes), monocytes, macrophages, and granulocytes (neutrophils, eosinophils, and basophils). The molecules of the immune system are the antibodies, complements, cytokines, interleukines, and interferons [1,2]. Immunity is the resistance of the host body to pathogen agents and their toxic effects. Every human body has a nonspecific response (innate immunity) and a specific response (acquired immunity) [1,2,3]. The innate immunity is based on mechanisms already existing before the pathogen agent infects the host and is the first line of the body’s defense, but it has no memory for subsequent exposure and is based on nonspecific mechanisms. Instead, the adaptive acquired immunity develops after the entrance of a pathogenic agent (virus, microbe, bacteria, parasites, or fungi) into the host and comes into action after innate immunity fails to get rid of these invaders. This acquired immunity has the memory to fight with pathogens in subsequent exposures and uses the specific cells (i.e., T cells (or cell-mediated) and B cells (or antibody-mediated)) [1,2,3].
Recently, a lot of modern epidemiological models, including mathematical models of the COVID-19 epidemic and other models for the spread of diseases such as Ebola, tuberculosis, or influenza, were studied by many experts in the dynamics of infectious diseases in [4,5,6,7].
In this paper, we consider the interaction between the human body’s immune system and a pathogenic invader, and we try to find an appropriate deterministic mathematical model in order to study the fight of the immune system with the invader [8,9]. These kinds of first-order systems of ordinary differential equations represent a class of Kolmogorov systems, and they are used very often in order to obtain mathematical models for the dynamics of populations, prey–predator ecological models, or other models for the spread of diseases or for the growth of tumors [10,11,12,13,14].
We will denote by x 1 ( t ) the level in time of the immune system (antibody level) and by v ( t ) the level of the virus. More exactly, v ( t ) represents the number of virus cells which exist in a body at time t. As soon as it notices the presence of a virus in the body, the immune system will fight with the virus (e.g., the present coronavirus, COVID-19) in order to stop its multiplication and eliminate it. Consider the time as being continuous, where t 0 . Thus, x ˙ 1 ( t ) and v ˙ ( t ) are the rates of change of these two quantities in a short unit of time, where x ˙ 1 ( t ) = d x 1 d t and v ˙ ( t ) = d v d t . Following the ideas of G. Moza from [15], this model is based on the following three hypotheses:
Hypothesis 1.
In the absence of a virus, the quantity of antibodies x 1 ( t ) can be present in the human body up to a threshold value K 1 . This hypothesis is based on the fact that the human body may have an innate immunity and also humoral and cellular immunity after prior possible contact with the virus. Thus, in the first stage, before a present contact with the virus v ( t ) , we consider that the evolution law of x 1 ( t ) is the following:
x ˙ 1 = a 1 x 1 b 1 x 1 2
with a 1 > 0 and b 1 > 0 . Taking into account that the solution of the prevoius equation with the initial condition x 1 ( 0 ) = x 10 is
x 1 ( t ) = K 1 1 + K 1 x 10 1 e a 1 t
then one can observe that x 1 ( t ) K 1 for t , where K 1 = a 1 b 1 is the threshold value.
If the term b 1 x 1 2 is missing, then x 1 ( t ) increases exponentially when t , because in this case ( b 1 = 0 ), the general solution of the equation in x 1 will be x 1 ( t ) = a 1 x 10 e t a 1 . Otherwise, if b 1 > 0 , then we find that the maximum threshold value of the antibody level x 1 is a 1 b 1 = K 1 .
Hypothesis 2.
Normally, in an healthy body without autoimmune diseases, the antibodies of the immune system do not attack other normal cells of the body. Thus, they can be destroyed only due to viruses, and as such, a term in the form c 1 x 1 v should be added to each equation in x 1 ( t ) , where v ( t ) is the number of virus cells which exist in the human body at time t.
Hypothesis 3.
In the absence of the immune system of the body, the virus would multiply indefinitely and exponentially, with v ( t ) satisfying the law v ˙ = p 2 v . However, of course, in the presence of the immune system, the number of virus cells will decrease, and consequently, we must add the terms p 1 x 1 v to the evolution law of v ( t ) .
If we denote the variable v by x 2 , then we can conclude that these hypotheses lead us to the following two-dimensional first-order differential system with five parameters, given by
x ˙ 1 = a 1 x 1 b 1 x 1 2 c 1 x 1 x 2 x ˙ 2 = p 2 x 2 p 1 x 1 x 2
where a 1 > 0 , b 1 > 0 , c 1 > 0 , p 1 > 0 , and p 2 > 0 .
Obviously, the model has medical relevance when x i 0 , i = 1 , 2 . Therefore, the solutions of the system lie in the set
Σ + 0 = ( x 1 , x 2 ) R 2 x i 0 , i = 1 , 2 .
Moreover, the lines { x i = 0 } are invariant manifolds with respect to the flow of the system (i.e., any orbit starting from a point which belongs to Σ + = ( x 1 , x 2 ) R 2 x i > 0 , i = 1 , 2 remains in Σ + ). Therefore, the orbits cannot cross any of these two invariant lines, and then the study of the system where it has medical relevance is well-defined in the sense that an orbit starting from a zone with medical relevance does not enter a zone with medical irrelevance, and vice versa.
This deterministic mathematical model of the interaction between the immune system and a pathogenic invader is based on the logistic differential equation
x ˙ = a x b x 2
or in the form
x ˙ = a x 1 x K , where K = a b
Next, using the same hypothesis, we can obtain another two models based on the Gompertz equation
x ˙ = a x ln K ln x
and generalized logistic equation (Richards equation):
x ˙ = a x 1 x K ν , where ν > 0 .
The graphs of the solutions of these ordinary differential equations belong to the so-called S-shaped curves (see Figure 1).
It can be remarked that only the logistic curve is symmetric with respect to the inflection point of 0.5 . The Gompertz curve has a lower inflection point, and the generalized logistic curve can have an inflection point over 0.5 or up to 0.5 , depending on the parameter ν .
The two-dimensional deterministic model represented by Equation (1) was developed using the logistic model in the second section, following the ideas from the three hypotheses of G. Moza in [15]. For the first time, in the third and the fourth sections, we will use the Gompertz model and the generalized logistic model (or Richards model) in order to construct a deterministic model for the interaction between the human immune system and a virus. Then, these models will be studied from the dynamical systems theory point of view, including the behavior and dynamics properties around the equilibrium points. The generalizations for higher dimensions of these three models will be introduced and studied in Section 5. For all three models, no matter the number of variables, the existence of transcritical bifurcations is obtained, but no Hopf-type bifurcations can occur. Therefore, there are no limit cycles. Finally, a section with relevant conclusions from the epidemiological and medical points of view will be presented.

2. The Logistic Model

The logistic function was introduced by Belgian mathematician Pierre-François Verhulst in a series of three papers [16,17,18] between the years 1838 and 1847. This model was created by Verhulst as a model of population growth by adjusting the exponential growth model through the introduction of a upper threshold called the carrying capacity.

2.1. A Short Presentation

The logistic equation is the following first-order differential equation:
x ˙ = a x b x 2
or equivalently
x ˙ = a x 1 x K
where K = a b is the carrying capacit and a is the inherent growth rate of the modeled population ( a > 0 , b > 0 ).
The solution of this differential equation with the initial condition x 0 = x ( 0 ) is the logistic function
x ( t ) = K 1 + K x 0 1 e a t
There are two equilibrium points for the dynamical system given by the logistic equation, namely x 1 = 0 and x 2 = K .
Since F ( x ) = a x 1 x K and d F d x ( x ) = a 1 2 x K , it follows that d F d x ( 0 ) = a > 0 and d F d x ( K ) = a < 0 (i.e., x 1 = 0 is a repeller, and x 2 = K is an attractor).

2.2. The Interaction between the Immune System and a Virus in a Logistic Model

We consider the following system of two differential equations:
x ˙ 1 = a 1 x 1 b 1 x 1 2 c 1 x 1 x 2 x ˙ 2 = p 2 x 2 p 1 x 1 x 2 ,
where x 1 > 0 represents the antibody level of the immune system and x 2 = v 0 represents the concentration of the virus in the body. Here, K 1 = a 1 b 1 is the carrying capacity of the antibodies, a 1 is his inherent growth rate, p 2 is the inherent growth rate of the virus, and c 1 , p 1 represent the interaction between the antibodies and virus. All parameters a i , b i , c i , and p i are strictly positive.
The Jacobi matrix at an equilibrium point ( x 1 , x 2 ) is
A = a 1 2 b 1 x 1 c 1 x 2 c 1 x 1 p 1 x 2 p 2 p 1 x 1
In order to find the equilibria, by analyzing the system
x 1 ( a 1 b 1 x 1 c 1 x 2 ) = 0 x 2 ( p 2 p 1 x 1 ) = 0
we obtain the following three equlibria: E 0 ( 0 , 0 ) with eigenvalues λ 1 = a 1 and λ 2 = p 2 , E 1 a 1 b 1 , 0 with eigenvalues λ 1 = a 1 and λ 2 = p 2 p 1 a 1 b 1 , and E 2 p 2 p 1 , 1 c 1 a 1 b 1 p 2 p 1 with eigenvalues
λ 1 , 2 = 1 2 p 1 b 1 p 2 ± b 1 2 p 2 2 + 4 p 2 p 1 2 a 1 b 1 p 2 p 1 .
Therefore, we have the following results:
Theorem 1. (a) 
The trivial equilibria E 0 ( 0 , 0 ) is always a repeller.
(b) 
The equilibria E 1 is an attractor if and only if p 2 p 1 < K 1 = a 1 b 1 .
Otherwise, E 1 is a saddle point.
(c) 
The equilibria E 2 is a saddle point whenever it lies on Σ + .
Moreover, E 2 Σ + if and only if p 2 p 1 < K 1 (i.e., E 1 is an attractor).
(d) 
The system does not undergo a Hopf bifurcation at E 2 on Σ + .
(e) 
E 2 collides with E 1 on the line S :   p 2 = p 1 K 1 .
(f) 
The equilibria E 2 bifurcates from the equilibria E 1 along the line
S = { ( p 1 , p 2 ) | p 2 = p 1 K 1 } by a transcritical bifurcation.

3. The Gompertz Model

Taking into account that when using the logistic equation, there are a lot of mathematical models for describing the behavior in time of the epidemiological or ecological systems, we will further introduce a new mathematical model for studying the interactions between different kinds of species using the Gompertz model. This model was introduced by British mathematician Benjamin Gompertz in 1825 in order to obtain a law of mortality as well as a demographic model (see [19]). One hundred forty years later in 1964, A.K. Laird used the Gompertz curve to fit data on the growth of tumors [20].

3.1. A Short Presentation

The Gompertz equation is the following first-order differential equation:
x ˙ = a x ln K x
or equivalently
x ˙ = a x ln K ln x
where K = a b is the carrying capacity and a is the inherent growth rate of the modeled population ( a > 0 , K > 0 ).
The solution of this differential equation with the initial condition x 0 = x ( 0 ) is the Gompertz function
x ( t ) = K e ln K x 0 e a t
There is a single equilibrium point for the dynamical system given by the Gompertz equation, namely x 1 = K .
Since F ( x ) = a x ln K ln x and d F d x ( x ) = a ln K x a , it follows that d F d x ( K ) = a < 0 (i.e., x 1 = K is an attractor).

3.2. The Interaction between the Immune System and a Virus in the Gompertz Model

We consider the following system of two differential equations:
x ˙ 1 = a 1 x 1 ln K 1 x 1 c 1 x 1 x 2 x ˙ 2 = p 2 x 2 p 1 x 1 x 2
where x 1 > 0 represents the antibody level of the immune system and x 2 = v 0 is the concentration of the virus in the body. Here, K 1 = a 1 b 1 is the carrying capacity of the antibodies, a 1 is the inherent growth rate, p 2 is the inherent growth rate of the virus, and c 1 and p 1 are the interaction’s coefficients between the antibodies and virus, respectively. All parameters a i , b i , c i , and p i are strictly positive.
The Jacobi matrix at an equilibrium point ( x 1 , x 2 ) is
A = a 1 ln K 1 x 1 a 1 c 1 x 2 c 1 x 1 p 1 x 2 p 2 p 1 x 1
In order to find the equilibrium points, we analyze the system
x 1 ( a 1 ln K 1 x 1 c 1 x 2 ) = 0 x 2 ( p 2 p 1 x 1 ) = 0
We obtain two equlibria: E 1 ( K 1 , 0 ) with eigenvalues λ 1 = a 1 and λ 2 = p 2 p 1 K 1 and E 2 p 2 p 1 , a 1 c 1 ln K 1 p 1 p 2 , which exists if and only if K 1 p 1 p 2 > 1 (i.e., p 2 p 1 < K 1 ).
Otherwise, if p 2 p 1 K 1 , the equilibrium E 2 does not exist (i.e., it does not belong to Σ + ).
For E 2 , we obtain the eigenvalues λ 1 , 2 = 1 2 a 1 ± 1 2 a 1 2 + 4 p 2 a 1 ln K 1 p 1 p 2 .
Since λ 1 + λ 2 = a 1 < 0 and λ 1 λ 2 = p 2 a 1 ln K 1 p 1 p 2 < 0 , it is found that E 2 is a saddle point.
In conclusion, we have the following results:
Theorem 2. (a) 
The equilibrium E 1 is an attractor if and only if p 2 p 1 < K 1 . Otherwise, E 1 is a saddle point.
(b) 
The equilibrium E 2 is a saddle point whenever it lies on Σ + . Moreover, E 2 Σ + if and only if p 2 p 1 < K 1 (i.e., E 1 is an attractor).
(c) 
The system does not undergo a Hopf bifurcation at E 2 on Σ + .
(d) 
E 2 collides with E 1 on the line S :   p 2 = p 1 K 1 .
(e) 
The equilibrium E 2 bifurcates from the equilibrium E 1 along the line
S = { ( p 1 , p 2 ) | p 2 = p 1 K 1 } by a transcritical bifurcation.

4. The Generalized Logistic Model

The generalized logistic function known as Richards’s function is a generalization of the logistic function, with more flexible S-shaped graph curves and an inflection point which depends on a parameter. This function is the solution of the Richards differential equation (RDE), otherwise known as the generalized logistic differential equation, and was used for modeling many growth phenomena from oncology and epidemiology. F. J. Richards proposed this growth function for the first time in 1959 in [21].

4.1. A Short Presentation

The generalized logistic equation (or Richards’s differential equation (RDE)) is the following first-order differential equation:
x ˙ = a x 1 x K ν
where K = a b is the carrying capacity, a is the inherent growth rate of the modeled population, and ν > 0 ( a > 0 , b > 0 ).
The solution of this differential equation with the initial condition x 0 = x ( 0 ) is the generalized logistic function or Richards function:
x ( t ) = K 1 + K x 0 ν 1 e a ν t 1 ν
There are two equilibrium points for the dynamical system given by the generalized logistic equation, namely x 1 = 0 and x 2 = K .
Since F ( x ) = a x 1 x K ν and d F d x ( x ) = a 1 2 x K ν , it follows that d F d x ( 0 ) = a > 0 and d F d x ( K ) = a < 0 (i.e., x 1 = 0 is a repeller, and x 2 = K is an attractor).
Let us remark that for ν = 1 , the classical logistic equation (Equation (3)) is just a particular case of the Richards differential equation (Equation (10)). Moreover, the Gompertz equation (Equation (7)) can be obtained from the generalized logistic equation (Equation (10)) when ν 0 , provided that a = O ( 1 ν ) . In fact, for ν > 0 , being small enough, we can take a = r ν , with r > 0 , and then the Richards equation x ˙ = a x 1 x K ν becomes
x ˙ = r ν x 1 x K ν = r x 1 e ν ln x K ν = r x ln ( K x )
if we take into account that 1 e ν Y ν = ν Y ν = Y for ν > 0 at a very small value. Here, the terms of order greater than two from e ν Y = n = 0 + ( ν Y ) n n ! are neglected because v is very small.

4.2. The Interaction between the Immune System and a Virus through the Generalized Logistic Model

We consider the following system of two differential equations:
x ˙ 1 = a 1 x 1 1 x 1 K 1 ν c 1 x 1 x 2 x ˙ 2 = p 2 x 2 p 1 x 1 x 2 ,
where x 1 > 0 represents the antibody level of the immune system and x 2 = v 0 is the concentration of the virus in the body. Here, K 1 = a 1 b 1 is the carrying capacity of the antibodies, a 1 is their inherent growth rate, p 2 is the inherent growth rate of the virus, and c 1 and p 1 are the interaction’s coefficients between the antibodies and the virus, respectively. All parameters a i , b i , c i , p i , and ν ( 0 , 1 ) are strictly positive.
The Jacobi matrix at an equilibrium point ( x 1 , x 2 ) is
A = a 1 a 1 ( 1 + ν ) x 1 K 1 ν c 1 x 2 c 1 x 1 p 1 x 2 p 2 p 1 x 1 .
In order to find the equilibrium points, by analyzing the system
x 1 a 1 1 x 1 K 1 ν c 1 x 2 = 0 x 2 ( p 2 p 1 x 1 ) = 0
we obtain the following three equilibria: E 0 ( 0 , 0 ) with eigenvalues λ 1 = a 1 and λ 2 = p 2 , E 1 K 1 , 0 with eigenvalues λ 1 = a 1 ν and λ 2 = p 2 p 1 K 1 , and E 2 p 2 p 1 , a 1 c 1 1 p 2 p 1 K 1 ν with eigenvalues λ 1 , 2 = 1 2 a 1 ν p 2 p 1 K 1 ν ± 1 2 a 1 2 ν 2 p 2 p 1 K 1 2 ν + 4 p 2 a 1 4 p 2 a 1 p 2 p 1 K 1 ν .
For E 2 , since λ 1 + λ 2 = a 1 ν e ν ln p 2 p 1 K 1 < 0 and λ 1 λ 2 = p 2 a 1 1 p 2 p 1 K 1 ν < 0 , it follow that E 2 is a saddle point whenever it belongs to Σ + .
Then, we have the following results:
Theorem 3. (a) 
The trivial equilibrium E 0 ( 0 , 0 ) is always a repeller.
(b) 
The equilibrium E 1 is an attractor if and only if p 2 p 1 < K 1 . Otherwise, E 1 is a saddle point.
(c) 
The equilibrium E 2 is a saddle point whenever it lies on Σ + . Moreover, E 2 Σ + if and only if p 2 p 1 < K 1 (i.e., E 1 is an attractor).
(d) 
The system does not undergo a Hopf bifurcation at E 2 on Σ + .
(e) 
E 2 collides with E 1 on the line S :   p 2 = p 1 K 1 .
(f) 
The equilibrium E 2 bifurcates from the equilibrium E 1 along the line
S = { ( p 1 , p 2 ) | p 2 = p 1 K 1 } by a transcritical bifurcation.

5. Generalizations for n 3

If we consider a mathematical model with interaction between two, three, or more types of human body immunities (or antibodies) and a virus, then we can similarly consider a system with three, four, or more differential equations using the logistic model, Gompertz model, or generalized logistic model.
For the three-dimensional model built by the logistic model, complete results can be found in the paper by G. Moza [15]. For the extended four-dimensional model, the results can be found in [22]. For a three-dimensional system, there are at most seven equilibrium points, among which only one can be an attractor. The rest of the equilibria are saddle points (if they exist), with the exception of the trivial equilibrium (the origin O), which is a repeller. More precisely, the equilibrium ( K 1 , K 2 , 0 ) is an attractor if p 3 < p 1 K 1 + p 2 K 2 , where K 1 and K 2 are the threshold values of the antibodies [15]. Similar results were obtained in [22] for the four-dimensional model, but the number of discovered equilibrium points was larger, having up to 15 equilibria. A generalization for higher dimensions can be discussed, but the number of equilibria becomes huge. For example, for a six-dimensional model, it will obtain up to 63 points of equilibrium [22].
However, if we consider the Gompertz model with two, three, or more types of human body immunities (or antibodies) and a virus, then the number of equilibria remains constant, and the obtained results do not depend of the number of variables of the system from the qualitative point of view.
For n = 3 , we consider the following system of three differential equations:
x ˙ 1 = a 1 x 1 ln K 1 x 1 c 1 x 1 x 3 x ˙ 2 = a 2 x 2 ln K 2 x 2 c 2 x 2 x 3 x ˙ 3 = p 3 x 3 p 1 x 1 x 3 p 2 x 2 x 3 ,
where x 1 , x 2 > 0 and x 3 = v 0 .
The Jacobi matrix at an equilibrium point ( x 1 , x 2 , x 3 ) is
A = a 1 ln K 1 x 1 a 1 c 1 x 3 0 c 1 x 1 0 a 2 ln K 2 x 2 a 2 c 2 x 3 c 2 x 2 p 1 x 3 p 2 x 3 p 3 p 1 x 1 p 2 x 2
In order to find the equilibrium points, by analyzing the system
x 1 ( a 1 ln K 1 x 1 c 1 x 3 ) = 0 x 2 ( a 2 ln K 2 x 2 c 2 x 3 ) = 0 x 3 ( p 3 p 1 x 1 p 2 x 2 ) = 0
we obtain the following two equilibria: E 1 ( K 1 , K 2 , 0 ) with eigenvalues λ 1 = a 1 , λ 2 = a 2 and λ 3 = p 3 p 1 K 1 p 2 K 2 and E 2 ( x 1 * , x 2 * , x 3 * ) , which exists if p 3 p 1 K 1 p 2 K 2 < 0 (i.e., E 1 is an attractor), where x 1 * , x 2 * , x 3 * are strictly positive real numbers which satisfy a i c i ln K i x i * = x 3 * or x i * = K i e x 3 * c i a i ( i = 1 , 2 ) and p 1 x 1 * + p 2 x 2 * = p 3 .
The characteristic polynomial at E 2 is
P ( λ ) = λ 3 + a 1 + a 2 λ 2 + a 1 a 2 x 3 * c 2 x 2 * p 2 + p 1 c 1 x 1 * λ x 3 * a 1 c 2 x 2 * p 2 + p 1 c 1 x 1 * a 2 .
Taking into account that α 2 = a 1 + a 2 > 0 , α 1 = a 1 a 2 x 3 * c 2 x 2 * p 2 + p 1 c 1 x 1 * , α 0 = x 3 * a 1 c 2 x 2 * p 2 + v p 1 c 1 x 1 * a 2 < 0 , and α 2 α 1 α 0 = a 1 a 2 ( a 1 + a 2 ) x 3 * p 1 c 1 x 1 * a 1 + c 2 x 2 * p 2 a 2 , following the Hurwitz criterion, the result is that E 2 is not an attractor.
Moreover, because we have λ 1 + λ 2 + λ 3 = ( a 1 + a 2 ) < 0 , λ 1 λ 2 + λ 2 λ 3 + λ 3 λ 1 = a 1 a 2 x 3 * c 2 x 2 * p 2 + p 1 c 1 x 1 * , λ 1 λ 2 λ 3 = x 3 * a 1 c 2 x 2 * p 2 + p 1 c 1 x 1 * a 2 > 0 , and λ 1 2 + λ 2 2 + λ 3 2 = ( a 1 + a 2 ) 2 2 a 1 a 2 x 3 * c 2 x 2 * p 2 + p 1 c 1 x 1 * = a 1 2 + a 2 2 + 2 c 2 x 2 * p 2 x 3 * + 2 p 1 x 3 * c 1 x 1 * > 0 , it follows that E 2 is a saddle point whenever it exists.
Next, we obtain the following results:
Theorem 4. (a) 
The equilibrium E 1 ( K 1 , K 2 , 0 ) is an attractor if and only if
p 3 p 1 K 1 p 2 K 2 < 0 . Otherwise, E 1 ( K 1 , K 2 , 0 ) is a saddle point.
(b) 
The equilibrium E 2 is a saddle point whenever it lies on Σ + . Moreover, E 2 Σ + if and only if p 3 p 1 K 1 p 2 K 2 < 0 (i.e., E 1 is an attractor).
(c) 
The system does not undergo a fold-Hopf bifurcation at E 2 .
(d) 
E 2 collides with E 1 on the surface S :   p 1 x 1 + p 2 x 2 = p 3 .
(e) 
The equilibrium E 2 bifurcates from E 1 along the surface
S = { ( p 1 , p 2 , p 3 ) | p 3 = p 1 K 1 + p 2 K 2 }
by a transcritical bifurcation.
For n = 4 , we consider the following system of four differential equations:
x ˙ 1 = a 1 x 1 ln K 1 x 1 c 1 x 1 x 4 x ˙ 2 = a 2 x 2 ln K 2 x 2 c 2 x 2 x 4 x ˙ 3 = a 3 x 3 ln K 3 x 3 c 3 x 3 x 4 x ˙ 4 = p 4 x 4 p 1 x 1 x 4 p 2 x 2 x 4 p 3 x 3 x 4 ,
where x 1 , x 2 , x 3 > 0 and x 4 = v 0 .
The Jacobi matrix at an equilibrium point ( x 1 , x 2 , x 3 , x 4 ) is
a 1 ln K 1 x 1 a 1 c 1 x 4 0 0 c 1 x 1 0 a 2 ln K 2 x 2 a 2 c 2 x 4 0 c 2 x 2 0 0 a 3 ln K 3 x 3 a 3 c 3 x 4 c 3 x 3 p 1 x 4 p 2 x 4 p 3 x 4 p 4 p 1 x 1 p 2 x 2 p 3 x 3
In order to find the equilibrium points, by analyzing the system
x 1 ( a 1 ln K 1 x 1 c 1 x 4 ) = 0 x 2 ( a 2 ln K 2 x 2 c 2 x 4 ) = 0 x 3 ( a 3 ln K 3 x 3 c 3 x 4 ) = 0 x 4 p 4 p 1 x 1 p 2 x 2 p 3 x 3 = 0
we obtain the following two equilibria: E 1 ( K 1 , K 2 , K 3 , 0 ) with eigenvalues λ 1 = a 1 , λ 2 = a 2 , λ 3 = a 3 , and λ 4 = p 4 p 1 K 1 p 2 K 2 p 3 K 3 and E 2 ( x 1 * , x 2 * , x 3 * , x 4 * ) , which exists if and only if p 4 p 1 K 1 p 2 K 2 p 3 K 3 < 0 (i.e., E 1 is an attractor), where x 1 * , x 2 * , x 3 * , and x 4 * are strictly positive real numbers which satisfy a i c i ln K i x i * = x 4 * or x i * = K i e x 4 * c i a i ( i = 1 , 2 , 3 ) and p 1 x 1 * + p 2 x 2 * + p 3 x 3 * = p 4 .
The characteristic polynomial at E 2 is P ( λ ) = λ 4 + α 3 λ 3 + α 2 λ 2 + α 1 λ + α 0 , where α 3 = a 1 + a 2 + a 3 , α 2 = a 1 a 2 + a 2 a 3 + a 3 a 1 x 4 * x 1 * c 1 p 1 + x 2 * c 2 p 2 + x 3 * c 3 p 3 , α 1 = a 1 a 2 a 3 x 4 * x 1 * c 1 p 1 a 2 + a 3 + x 2 * c 2 p 2 a 1 + a 3 + x 3 * c 3 p 3 a 1 + a 2 , and α 0 = x 4 * x 1 * c 1 p 1 a 2 a 3 + x 2 * c 2 p 2 a 1 a 3 + x 3 * c 3 p 3 a 1 a 2 = det A .
Following the relations of Viète for P ( λ ) , we have λ 1 λ 2 λ 3 λ 4 = α 0 < 0 and λ 1 + λ 2 + λ 3 + λ 4 = α 3 < 0 . Then, it follows that E 2 is a saddle point whenever it exists.
Moreover, the characteristic polynomial at E 2 has at least two real eigenvalues with different signs. Therefore, we have the following results:
Theorem 5. (a) 
The equilibrium E 1 ( K 1 , K 2 , K 3 , 0 ) is an attractor if and only if p 4 p 1 K 1 p 2 K 2 p 3 K 3 < 0 . Otherwise, E 1 ( K 1 , K 2 , K 3 , 0 ) is a saddle point.
(b) 
The equilibrium E 2 is a saddle point whenever it lies on Σ + . Moreover, E 2 Σ + if and only if p 4 p 1 K 1 p 2 K 2 p 3 K 3 < 0 (i.e., E 1 is an attractor).
(c) 
The system does not undergo a Hopf-Hopf bifurcation at E 2 .
(d) 
E 2 collides with E 1 on the hypersurface S :   p 4 = p 1 x 1 + p 2 x 2 + p 3 x 3 .
(e) 
The equilibrium E 2 bifurcates from E 1 along the hypersurface
S = { ( p 1 , p 2 , p 3 , p 4 ) | p 4 = p 1 x 1 + p 2 x 2 + p 3 x 3 }
by a transcritical bifurcation.
Finally, of course, we can write generalized results available for any n 3 :
Theorem 6. (a) 
E 1 ( K 1 , , K n 1 , 0 ) is an attractor if and only if p n < i = 1 n 1 p i K i . Otherwise, E 1 is a saddle point.
(b) 
The equilibrium E 2 is a saddle point whenever it lies on Σ + (i.e., E 1 is an attractor).
Let us remark that the coordinates of the second equilibrium E 2 ( x 1 * , , x n 1 * , x n * ) are strictly positive real numbers x i which satisfy x i * = K i e x n * c i a i ( i = 1 , n 1 ¯ ) and p 1 x 1 * + + p n 1 x n 1 * = p n . Moreover, x i * < K i for all i = 1 , n 1 ¯ .
Obviously, if we want to use the generalized logistic model (Richards model) for the interaction between two, three, or more types of human body immunities (or antibodies) and a virus, then the results will be similar, but the computations are very complicated. This study can be conducted in a future work.

6. Conclusions

Taking into account that our body has at least three types of immunity (innate immunity, humoral immunity and cellular immunity), and these three types of immunity fight jointly with the pathogenic invader (viruses, bacteria, parasites, or fungi) in order to stop its multiplication and eliminate it, in this work, I presented a study about the interactions between the antibodies, cells, and organs of the immune system, as well as other complex human body mechanisms (such as the complement system of the body), and a pathogenic agent, such as COVID-2019. A mathematical approach based on a first-order system of differential equations was used for modeling the interactions, as well as tools from dynamical systems theory for analysis of the models, following the ideas from [15,22]. This study reveals the importance of antibodies in the fight against viruses. Several conclusions relevant to the medical world arise from our study, which are as follows:
1. If the immune system is sufficiently weak when the virus starts to proliferate (e.g., in the neighborhood of the trivial equilibria and the origin O)—that is, the level of antibodies is very small, but the virus concentration is also very small—then the virus has a strong chance to win. This occurs in the neighborhood of the origin O, because O is a repeller (for the logistic and generalized logistic models), and then any orbit γ ( t ) starting at a point x 0 Σ + close to O will depart from O for a large t value, meaning that v ( t ) may progress to infinity when t is large. Therefore, a serious deficiency in the level of immunity in the early stages of virus proliferation may lead to the virus’s victory.
2. If the immunity system is within its normal level (at the carrying capacity) from the first moment it discovers the virus, and if the immune system is in a healthy condition to kill the virus at a high rate ( K 1 > p 2 p 1 ), then the immune system has the best chance to win the battle with the virus. See the case of the first equilibrium point E 1 ( K 1 , 0 ) when it is an attractor. Indeed, any orbit γ ( t ) starting at a point x 0 Σ + close to E 1 ( K 1 , 0 ) will converge to E 1 for large t values; that is, v ( t ) trends toward zero for large t values.
3. If the levels of immunity becomes considerably smaller than their normal concentrations ( x 1 < K 1 ) at a moment during the battle with the virus, then the virus may win even though the immune system kills the virus at a rate higher than the rate of virus proliferation ( K 1 > p 2 p 1 ). This is the case of the second equilibrium point E 2 ( x 1 , v ) , x 1 = p 2 p 1 , v > 0 , which is a saddle point whenever it exists. More precisely, any orbit γ ( t ) starting at a point x 0 Σ + close to E 2 , x 0 W E 2 s will depart from E 2 for a large t value, which means v ( t ) may progress to infinity when t is large. A stable limit cycle around E 2 cannot arise through a Hopf bifurcation since all eigenvalues are real.
4. If the immunity system is within its normal level (at the carrying capacity) from the first moment when it discovers the virus, but the immune system is not able to eliminate the virus at a high rate ( K 1 < p 2 p 1 ), then the immune system can lose the battle with the virus. This is exactly the situation when the first equilibrium point E 1 is a saddle point, and then any orbit γ ( t ) starting at a point x 0 Σ + close to E 1 , x 0 W E 1 s will depart from E 1 for large t values, that means v ( t ) may progress to infinity when t is large. A stable limit cycle around E 1 cannot arise through a Hopf bifurcation since all eigenvalues are real.
5. If the first equilibrium point E 1 ( K 1 , 0 ) (for which the level of the immunity is at its maximum at the thresholds K 1 , but the virus concentration is very, very low) is an attractor, then we can conclude that the immune system eliminates the virus for any orbit starting from a point belonging to a small enough neighborhood of E 1 . However, the virus has another opportunity in the neighborhood of the second equilibrium point E 2 , which is a saddle point whenever it exists. What is more, this second equilibrium exists if and only if the first equilibrium is an attractor.
Although only the logistic model of the three models has symmetry with respect to the inflection point, these conclusions are available no matter what model is used: logistic, Gompertz, or the generalized logistic model (Richards model). For all models, regardless of the dimensions, the existence of transcritical bifurcations is obtained, but no Hopf-type bifurcations can occur. Therefore, there are no limit cycles.
The results we obtained in this work are natural and tell the medical community to work harder on methods for strengthening the immune system to win battles with pathogenic viruses.

Funding

This research was funded by the University of Craiova, Romania.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was partially supported by the Horizon2020-2017-RISE-777911 project.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Available online: https://www.coursehero.com/ (accessed on 1 February 2021).
  2. Hoffman, W.; Lakkis, F.G.; Chalasani, G.B. Cells, antibodies and more. Clin. J. Am. Soc. Nephrol. 2016, 11, 137–154. [Google Scholar] [CrossRef] [PubMed]
  3. Chaplin, D.D. Overview of the immune response. J. Allergy Clin. Immunol. 2010, 125 (Suppl. 2), 3–23. [Google Scholar] [CrossRef] [PubMed]
  4. Trejos, D.Y.; Valverde, J.C.; Venturino, E. Dynamics of infectious diseases: A review of the main biological aspects and their mathematical translation. Appl. Math. Nonlinear Sci. 2021, 7, 1–26. [Google Scholar] [CrossRef]
  5. Kumar, A.; Prakash, A.; Baskonus, H.M. The epidemic COVID-19 model via Caputo–Fabrizio fractional operator. Waves Random Complex Media 2022. [Google Scholar] [CrossRef]
  6. Gao, W.; Baskonus, H.M. Deeper investigation of modified epidemiological computer virus model containing the Caputo operator. Chaos Solitons Fractals 2022, 158, 112050. [Google Scholar] [CrossRef]
  7. Dubey, R.S.; Goswami, P.; Baskonus, H.M.; Tailor Gomati, A. On the existence and uniqueness analysis of fractional blood glucose-insulin minimal model. Int. J. Model. Simul. Sci. 2022. [Google Scholar] [CrossRef]
  8. Brauer, F.; Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  9. Freedman, H.I. Deterministic Mathematical Models in Population Biology; Marcel Dekker: New York, NY, USA, 1980. [Google Scholar]
  10. Bacaër, N. Mathématiques et épidémies; Cassini: Paris, France, 2021; 320p. (In French)
  11. Kuznetsov, Y.A. Elements of Applied Bifurcation Theory, 2nd ed.; Applied Mathematical Sciences; Springer: New York, NY, USA, 1998; Volume 112. [Google Scholar]
  12. Perko, L. Differential Equations and Dynamical Systems, 3rd ed.; Springer: New York, NY, USA, 2001. [Google Scholar]
  13. Moza, G.; Lazureanu, C.; Munteanu, F.; Sterbeti, C.; Florea, A. Bifurcation diagrams in a class of Kolmogorov systems. Nonlinear Anal. RWA 2020, 56, 103154. [Google Scholar]
  14. Moza, G.; Lazureanu, C.; Munteanu, F.; Sterbeti, C.; Florea, A. Analysis of a class of Kolmogorov systems. Nonlinear Anal. RWA 2021, 57, 103202. [Google Scholar]
  15. Moza, G. Studying the Interactions between the Immune System and a Virus like COVID-19, Politehnica University of Timişoara: Timişoara, Romania, 2021; preprint.
  16. Verhulst, P.F. Notice sur la loi que la population poursuit dans son accroissement. Corresp. Math. Phys. 1838, 10, 113–121. [Google Scholar]
  17. Verhulst, P.F. Recherches mathématiques sur la loi d’accroissement de la population. Nouv. Mém. l’Acad. R. Sci. Belles-Lettres Brux. 1845, 18, 8. [Google Scholar]
  18. Verhulst, P.F. Deuxième mémoire sur la loi d’accroissement de la population. Mém. l’Acad. R. Sci. Lettres Beaux-Arts Belg. 1847, 20, 1–32. [Google Scholar]
  19. Gompertz, B. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. Philos. Trans. R. Soc. Lond. 1825, 115, 513–585. [Google Scholar]
  20. Laird, A.K. Dynamics of Tumor Growth. Br. J. Cancer 1964, 13, 490–502. [Google Scholar] [CrossRef]
  21. Richards, F.J. A Flexible Growth Function for Empirical Use. J. Exp. Bot. 1959, 10, 290–300. [Google Scholar] [CrossRef]
  22. Munteanu, F. A 4-Dimensional Mathematical Model for Interaction between the Human Immune System and a Virus. Preprints 2022, 2022070282. [Google Scholar] [CrossRef]
Figure 1. Graphs for the logistic model y = 1 1 + e x (red), for the Gompertz model y = e e x (blue), for the generalized logistic model with ν = 2 , y = 1 1 + e 2 x 1 / 2 (orange), and with ν = 1 2 , y = 1 1 + e ( 1 / 2 ) x 2 (green).
Figure 1. Graphs for the logistic model y = 1 1 + e x (red), for the Gompertz model y = e e x (blue), for the generalized logistic model with ν = 2 , y = 1 1 + e 2 x 1 / 2 (orange), and with ν = 1 2 , y = 1 1 + e ( 1 / 2 ) x 2 (green).
Symmetry 14 01594 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Munteanu, F. A Comparative Study of Three Mathematical Models for the Interaction between the Human Immune System and a Virus. Symmetry 2022, 14, 1594. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14081594

AMA Style

Munteanu F. A Comparative Study of Three Mathematical Models for the Interaction between the Human Immune System and a Virus. Symmetry. 2022; 14(8):1594. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14081594

Chicago/Turabian Style

Munteanu, Florian. 2022. "A Comparative Study of Three Mathematical Models for the Interaction between the Human Immune System and a Virus" Symmetry 14, no. 8: 1594. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14081594

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