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
the level in time of the immune system (antibody level) and by
the level of the virus. More exactly,
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
. Thus,
and
are the rates of change of these two quantities in a short unit of time, where
and
. 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 can be present in the human body up to a threshold value . 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 , we consider that the evolution law of is the following:with and . Taking into account that the solution of the prevoius equation with the initial condition isthen one can observe that for , where is the threshold value. If the term is missing, then increases exponentially when , because in this case (), the general solution of the equation in will be . Otherwise, if , then we find that the maximum threshold value of the antibody level is .
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 should be added to each equation in , where 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 satisfying the law . However, of course, in the presence of the immune system, the number of virus cells will decrease, and consequently, we must add the terms to the evolution law of .
If we denote the variable
v by
, then we can conclude that these hypotheses lead us to the following two-dimensional first-order differential system with five parameters, given by
where
,
,
,
, and
.
Obviously, the model has medical relevance when
,
. Therefore, the solutions of the system lie in the set
Moreover, the lines are invariant manifolds with respect to the flow of the system (i.e., any orbit starting from a point which belongs to 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
or in the form
Next, using the same hypothesis, we can obtain another two models based on the Gompertz equation
and generalized logistic equation (Richards equation):
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 . The Gompertz curve has a lower inflection point, and the generalized logistic curve can have an inflection point over or up to , 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:
or equivalently
where
is the carrying capacit and
a is the inherent growth rate of the modeled population (
,
).
The solution of this differential equation with the initial condition
is the logistic function
There are two equilibrium points for the dynamical system given by the logistic equation, namely and .
Since and , it follows that and (i.e., is a repeller, and 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:
where
represents the antibody level of the immune system and
represents the concentration of the virus in the body. Here,
is the carrying capacity of the antibodies,
is his inherent growth rate,
is the inherent growth rate of the virus, and
,
represent the interaction between the antibodies and virus. All parameters
,
,
, and
are strictly positive.
The Jacobi matrix at an equilibrium point
is
In order to find the equilibria, by analyzing the system
we obtain the following three equlibria:
with eigenvalues
and
,
with eigenvalues
and
, and
with eigenvalues
Therefore, we have the following results:
- Theorem 1. (a)
The trivial equilibria is always a repeller.
- (b)
The equilibria is an attractor if and only if .
Otherwise, is a saddle point.
- (c)
The equilibria is a saddle point whenever it lies on .
Moreover, if and only if (i.e., is an attractor).
- (d)
The system does not undergo a Hopf bifurcation at on .
- (e)
collides with on the line .
- (f)
The equilibria bifurcates from the equilibria along the line
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:
or equivalently
where
is the carrying capacity and
a is the inherent growth rate of the modeled population (
,
).
The solution of this differential equation with the initial condition
is the Gompertz function
There is a single equilibrium point for the dynamical system given by the Gompertz equation, namely .
Since and , it follows that (i.e., 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:
where
represents the antibody level of the immune system and
is the concentration of the virus in the body. Here,
is the carrying capacity of the antibodies,
is the inherent growth rate,
is the inherent growth rate of the virus, and
and
are the interaction’s coefficients between the antibodies and virus, respectively. All parameters
,
,
, and
are strictly positive.
The Jacobi matrix at an equilibrium point
is
In order to find the equilibrium points, we analyze the system
We obtain two equlibria: with eigenvalues and and , which exists if and only if (i.e., ).
Otherwise, if , the equilibrium does not exist (i.e., it does not belong to ).
For , we obtain the eigenvalues .
Since and , it is found that is a saddle point.
In conclusion, we have the following results:
- Theorem 2. (a)
The equilibrium is an attractor if and only if . Otherwise, is a saddle point.
- (b)
The equilibrium is a saddle point whenever it lies on . Moreover, if and only if (i.e., is an attractor).
- (c)
The system does not undergo a Hopf bifurcation at on .
- (d)
collides with on the line .
- (e)
The equilibrium bifurcates from the equilibrium along the line
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:
where
is the carrying capacity,
a is the inherent growth rate of the modeled population, and
(
,
).
The solution of this differential equation with the initial condition
is the generalized logistic function or Richards function:
There are two equilibrium points for the dynamical system given by the generalized logistic equation, namely and .
Since and , it follows that and (i.e., is a repeller, and is an attractor).
Let us remark that for
, 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
, provided that
. In fact, for
, being small enough, we can take
, with
, and then the Richards equation
becomes
if we take into account that
for
at a very small value. Here, the terms of order greater than two from
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:
where
represents the antibody level of the immune system and
is the concentration of the virus in the body. Here,
is the carrying capacity of the antibodies,
is their inherent growth rate,
is the inherent growth rate of the virus, and
and
are the interaction’s coefficients between the antibodies and the virus, respectively. All parameters
,
,
,
, and
are strictly positive.
The Jacobi matrix at an equilibrium point
is
In order to find the equilibrium points, by analyzing the system
we obtain the following three equilibria:
with eigenvalues
and
,
with eigenvalues
and
, and
with eigenvalues
.
For , since and , it follow that is a saddle point whenever it belongs to .
Then, we have the following results:
- Theorem 3. (a)
The trivial equilibrium is always a repeller.
- (b)
The equilibrium is an attractor if and only if . Otherwise, is a saddle point.
- (c)
The equilibrium is a saddle point whenever it lies on . Moreover, if and only if (i.e., is an attractor).
- (d)
The system does not undergo a Hopf bifurcation at on .
- (e)
collides with on the line .
- (f)
The equilibrium bifurcates from the equilibrium along the line
by a transcritical bifurcation.
5. Generalizations for
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
is an attractor if
, where
and
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
, we consider the following system of three differential equations:
where
,
and
.
The Jacobi matrix at an equilibrium point
is
In order to find the equilibrium points, by analyzing the system
we obtain the following two equilibria:
with eigenvalues
,
and
and
, which exists if
(i.e.,
is an attractor), where
,
,
are strictly positive real numbers which satisfy
or
(
) and
.
The characteristic polynomial at
is
Taking into account that , , , and , following the Hurwitz criterion, the result is that is not an attractor.
Moreover, because we have , , , and , it follows that is a saddle point whenever it exists.
Next, we obtain the following results:
- Theorem 4. (a)
The equilibrium is an attractor if and only if
. Otherwise, is a saddle point.
- (b)
The equilibrium is a saddle point whenever it lies on . Moreover, if and only if (i.e., is an attractor).
- (c)
The system does not undergo a fold-Hopf bifurcation at .
- (d)
collides with on the surface .
- (e)
The equilibrium bifurcates from along the surfaceby a transcritical bifurcation.
For
, we consider the following system of four differential equations:
where
,
,
and
.
The Jacobi matrix at an equilibrium point
is
In order to find the equilibrium points, by analyzing the system
we obtain the following two equilibria:
with eigenvalues
,
,
, and
and
, which exists if and only if
(i.e.,
is an attractor), where
,
,
, and
are strictly positive real numbers which satisfy
or
(
) and
.
The characteristic polynomial at is , where , ,
, and
.
Following the relations of Viète for , we have and . Then, it follows that is a saddle point whenever it exists.
Moreover, the characteristic polynomial at has at least two real eigenvalues with different signs. Therefore, we have the following results:
- Theorem 5. (a)
The equilibrium is an attractor if and only if . Otherwise, is a saddle point.
- (b)
The equilibrium is a saddle point whenever it lies on . Moreover, if and only if (i.e., is an attractor).
- (c)
The system does not undergo a Hopf-Hopf bifurcation at .
- (d)
collides with on the hypersurface .
- (e)
The equilibrium bifurcates from along the hypersurfaceby a transcritical bifurcation.
Finally, of course, we can write generalized results available for any :
- Theorem 6. (a)
is an attractor if and only if . Otherwise, is a saddle point.
- (b)
The equilibrium is a saddle point whenever it lies on (i.e., is an attractor).
Let us remark that the coordinates of the second equilibrium are strictly positive real numbers which satisfy () and . Moreover, for all .
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 starting at a point close to O will depart from O for a large t value, meaning that 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 (), then the immune system has the best chance to win the battle with the virus. See the case of the first equilibrium point when it is an attractor. Indeed, any orbit starting at a point close to will converge to for large t values; that is, trends toward zero for large t values.
3. If the levels of immunity becomes considerably smaller than their normal concentrations () 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 (). This is the case of the second equilibrium point , , , which is a saddle point whenever it exists. More precisely, any orbit starting at a point close to , will depart from for a large t value, which means may progress to infinity when t is large. A stable limit cycle around 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 (), then the immune system can lose the battle with the virus. This is exactly the situation when the first equilibrium point is a saddle point, and then any orbit starting at a point close to , will depart from for large t values, that means may progress to infinity when t is large. A stable limit cycle around cannot arise through a Hopf bifurcation since all eigenvalues are real.
5. If the first equilibrium point (for which the level of the immunity is at its maximum at the thresholds , 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 . However, the virus has another opportunity in the neighborhood of the second equilibrium point , 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.