Abstract
Platelets have been seen traditionally as fragments of blood mediating coagulation. However, evidence during malaria infection suggests that platelets also act against merozoites, an infectious form of malaria in the bloodstream, and megakaryocytes can release giant platelets with a larger volume than normal platelets. We propose a mathematical model to study the interaction between red blood cells, merozoites, and platelets during malaria infection. We analyzed two cases of the interaction of platelets with malaria infection. In the first one, we considered the isolated action of normal platelets and, in the second one, the joint antiparasitic action of both normal and giant platelets. Numerical simulations were performed to evaluate the stability of the equilibrium points of the system of equations. The model showed that the isolated antiparasitic action of normal platelets corroborates malaria infection control. However, the system can converge to a presence-merozoite equilibrium point, or an oscillatory behavior may appear. The joint antiparasitic action of both normal and giant platelets eliminated the oscillatory behavior and drove the dynamics to converge to lower parasitic concentration than the case of isolated action of normal platelets. Moreover, the joint antiparasitic action of platelets proved more easily capable of eliminating the infection.
1 Introduction
In tropical and subtropical regions, diseases transmitted by mosquitoes represent a major public health problem, and malaria is one of them. In 2019, there were 409 thousand deaths by malaria reported by the World Health Organization (WHO). Malaria is caused by protozoan belonging to the genus Plasmodium, which is transmitted by the female mosquito Anopheles, being among the diseases that kill the most in the world [1].
During Anopheles bite, sporozoites are introduced into the human bloodstream and captured by hepatocytes. After a time in the liver, the sporozoites differentiate into merozoites, which is the parasite form that penetrates erythrocytes and starts the pathogenic phase of malaria infection [5,9]. Inside this infected erythrocyte, after a replication time, more merozoites are released to the bloodstream, and they can infect other naive erythrocytes. This cycle of erythrocyte infection repeats at regular intervals, which is different for each species. The periodic release of merozoites causes typical symptoms like chills, sweating, and fever [6].
Platelets have been traditionally seen as fragments of blood mediating coagulation. However, evidence suggests that platelets are critical in immunology, particularly innate immunity [2,16,19]. They preferentially bind to the infected erythrocytes, and this junction is associated with the death of the parasite [12]. This binding involves interactions between the platelets CD36 molecules and specifically expressed molecules in infected erythrocytes [11]. After binding to infected erythrocytes by merozoites, the platelets are specifically activated and release α-granules. The α-granule protein, mainly responsible for platelets’ cytocidal action, is the platelet factor 4, denoted by PF4 [13].
Mc Morran et al. (2012) found that both PF4 and Duffy antigen (expressed in erythrocytes) are required to kill the parasite by the action of platelets [13]. Recent studies have shown that PF4 is specifically internalized by the parasite and kills inside the cell. The Duffy antigen is equally essential, and it was discovered that in erythrocytes that lack Duffy antigen, the parasites were not affected by platelets [11].
Lacerda (2010) suggested that during malaria infection, the megakaryocytes can release giant platelets into the bloodstream to compensate for the low absolute number of peripheral platelets. These giant platelets appear to have a shorter time of maturation in detriment of the normal platelets [8]. Thus, platelets, directly on the merozoites or indirectly due to the platelets’ adhesion in the infected cells, may help to control the infection of erythrocyte by merozoites. Here we propose a mathematical model composed of ordinary differential equations to investigate the interaction between platelets and the course of malaria infection. We emphasize that other control mechanisms will not be considered; only the investigation of platelets will be done as an auxiliary response in the controlling of erythrocyte merozoite infection.
The paper is structured as follows. In Section 2, a mathematical model is presented to describe platelets’ action on malaria infection. In Section 3, the model is analyzed in two parts: (i) the model without antiparasitic action of normal and giant platelets, (ii) the effect of antiparasitic action of both normal and giant platelets against the parasite. Section 4 presents a discussion. The conclusion is given in Section 5.
2 Model formulation
We assume that uninfected red blood cells, whose concentration is denoted by H, are maintained at a constant level (homeostasis) by a constant production rate K1 and death µH. These cells, interacting with free merozoites, which are denoted by X, are infected at a constant infection rate β.
The infected red blood cells, denoted by HI, have an additional mortality rate µI due to infection. Each infected red blood cell releases α new merozoites in the bloodstream during its average lifetime. The merozoites have a death rate µX, and δ1 represents the number of released merozoites that invade uninfected red blood cells.
To describe platelets’ production, we assume a compartment of megakaryocyte precursor cells, denoted by Z. We admitted that megakaryocyte precursor cells are produced at a constant rate K2 and have a death rate µZ. After a time, they differentiate into megakaryocytes, which concentration is denoted by M1. This megakaryocyte have death rate µ1 and, after a ε1−1 time, releases δ2 normal platelets, which concentration is denoted by P1, and death rate .
We consider that, in the presence of infected red blood cells in the bloodstream, with a rate σ2, the megakaryocyte precursor cell can also be differentiated into a different type of megakaryocyte, denoted by M2, and produce giant platelets, denoted by P2. After time (less than ), the megakaryocyte M2 releases δ3 giant platelets with a larger volume compared to normal platelets. These giant platelets appear during malaria infection to assist the normal platelets, and they have a death rate .
We also assume that some megakaryocytes can be recruited to produce giant platelets during malaria infection. In this way, we consider that ε2 represents the recruitment rate of megakaryocyte M1 to produce giant platelets.
We consider that both normal and giant platelets can act against infected red blood cells to describe platelets’ immune action. The direct action of normal platelets (giant platelets) in merozoites and infected red blood cells is given by the per-capita rate of antiparasitic action γX and γa , respectively. The parameter δ4 (δ5) represents the number of ordinary platelets (giant platelets) bind to the infected red blood cell, and δ6 (δ7) represents the number of common platelets (giant platelets) bind to the merozoite.
Based on the above assumptions, we have the following system of differential equations: Table 1 summarizes the variables, and Table 2 summarizes the parameters assumed in the model described by the system of equation (1).
3 Analysis
We performed the steady-state analysis in order to evaluate the platelet antiparasitic action during malaria infection. This section is divided into two main parts. In the first part, we study the case where there is no antiparasitic activity of normal or giant platelets. Hence, there are only normal platelets circulating in the bloodstream without interfering with merozoite infection. That is, malaria infection interacts only with erythrocytes. This part is denoted by “Non-interacting platelets and malaria infection”.
In the second part, we analyze both normal and giant platelet antiparasitic activity, how they interact with malaria infection. This part is denoted by “Platelets antiparasitic action “. To investigate the individual contribution of normal platelets and the joint action on the infection control, we analyze the two cases: (i) normal platelet acts against merozoites and infected red blood cells, and (ii) joint action normal and giant platelets.
3.1 Non-interacting platelets and malaria infection
Here, platelets do not have antiparasitic action. Therefore, we have two decoupled systems, a system that models the dynamics of normal platelet production and another related to the interaction of merozoites with the target cells.
For the system that models the dynamics of platelet production, we have In this case, we have only one equilibrium point given by , where , and .
By the change of variables , and , we can prove that T is globally asymptotically stable. For more details, see Appendix A.
The system that models the interaction between merozoites and red blood cells without antiparasitic action of platelets is given by The merozoites-free equilibrium (MFE), , is given by W 0 = (K1/µH, 0, 0). The characteristic polynomial corresponding to the Jacobian matrix evaluated at W 0 has one negative root −µH and two roots given by the solutions of the second-order polynomial where the coefficients are, and R0 is given by, Applying the Routh-Hurwitz criteria [15], MFE is locally asymptotically stable if R0 < 1 and unstable if R0 > 1. Therefore, the local stability of this point is stated in Theorem 1.
The merozoites-free equilibrium point W 0 is locally asymptotically stable if R0 < 1, and unstable if R0 > 1, where R0 is given by (5).
Proof: The polynomial q(λ), which is given by equation (4), satisfies all Routh-Hurwitz conditions if R0 < 1, that is, q1 > 0 and q0 > 0. □
With respect to R0, it can be rewritten as , where Notice that [(β K1/µH)α/µX] represents the number of new merozoites produced by one merozoite introduced into the bloodstream (containing a susceptible erythrocyte population of size K1/µH) free of this parasite. The second term, , is the number of merozoites consumed by the infected erythrocyte during the cell infection process.
Therefore, R0 is the net number of merozoites produced by invading merozoites in the early stage of malaria infection.
The merozoites-presence equilibrium (MPE), which is , has the coordinates given by Therefore, MPE exists only if R0 > 1 due to the biological condition of existence. The characteristic polynomial corresponding to the Jacobian matrix evaluated at W ∗ is given by where The local stability of W ∗ is stated in Theorem 2.
The merozoites-presence equilibrium point W ∗ is locally asymptotically stable if R0 > 1, and unstable if R0 < 1.
Proof. See Appendix B.
To obtain global stability of W 0 and W ∗, we used the direct method of Lyapunov. The proofs were adapted from Gómez (2016) [4].
If R0 < 1 then merozoite-free equilibrium W 0 is globally asymptotically stable.
Proof. See Appendix C.
The merozoite-presence equilibrium W ∗ is globally asymptotically stable if R0 > 1 and .
Proof. See Appendix D.
We showed that MPE is globally asymptotically stable only if R0 > 1 and . The additional restriction, , appears because the model considered the cell invasion by the parasite necessary to occur the erythrocyte infection. If we had considered only the contact between the parasite and the erythrocyte to occur the infection, there was no term δ1β HX. Thus, we would have a classic model of Mathematical Epidemiology, and W ∗ would be globally asymptotically stable, without additional restriction, only R0 > 1.
3.2 Platelets antiparasitic action
Here, we consider the antiparasitic action of both normal and giant platelets. The merozoites-free steady state is given by , where the nonnull coordinates values are the same as appeared in W 0 and .
The characteristic polynomial corresponding to the Jacobian matrix evaluated at S0 has six negative roots: and −µH; and two roots are given by the solutions of second-order polynomial w(λ) = λ2 + w1λ + w0, where wherein Rc is given by and .
Applying the Routh-Hurwitz criteria, S0 is locally asymptotically stable if Rc < 1, and unstable if Rc > 1.
The merozoites-free equilibrium point S0 is locally asymptotically stable if Rc > 1, and unstable if Rc < 1.
Proof. The proof is analogue of the proof in Theorem 1. So, as w1 and w0 are positive values, the equilibrium point S0 is locally asymptotically stable if Rc < 1 and unstable if Rc > 1.□
Let us interpret the parameter Rc given by (9). Considering , where and are defined by The parameter has the following interpretation. Consider that one merozoite during its lifespan , which suffers a decrease due to normal platelets’ action, infects with a rate β one susceptible red blood cell in the bloodstream (containing a susceptible erythrocyte population of size K1/µH). This infected erythrocyte has a probability of surviving during its lifespan. If this occurs, this cell will release to the bloodstream α new merozoites. To the parameter , we have the following interpretation, δ1 merozoites during their lifespan at a rate β infect one red blood cell in the bloodstream (containing a susceptible erythrocyte population of size K1/µH). This infected erythrocyte has a probability of suffering antiparasitic action.
As a result, the threshold Rc is the net number of merozoites produced by one invading merozoite on the occurrence of antiparasitic action of normal platelets.
Observe that, if , then we have Rc < 0. Thus, the action of normal platelets on infected red blood cells is much more efficient than their action on the merozoites, as it is able to make the threshold Rc negative. The negative value of Rc shows that the infection would already be controlled because no new merozoite would be produced from the primary infection. It is also important to emphasize that if all Rc factors are positives, then Rc is less than R0. It shows that the antiparasitic action of normal platelets decreases the net number of merozoites produced by one invading merozoite.
The merozoites-presence equilibrium S∗ = (H∗, X ∗, HI∗, Z∗, M1∗, M2∗, P1∗, P2∗) has the coordinates given by where X ∗ and are solution of the system of equations Due to the system (11) be strongly nonlinear, we did not obtain all coordinates of S∗. So, in the next subsections, we analyze the existence and stability of S∗ restricted to two different cases: (i) Antiparasitic action of normal platelets without giant platelets and (ii) Joint antiparasitic action of both normal and giant platelets.
3.2.1 Antiparasitic action of normal platelets without giant platelets
In this case, only normal platelets have antiparasitic action, and the giant platelets are not being produced. Simplifying the system of equations (11), we obtain S∗ = (H∗, X ∗, HI∗, Z∗, M1∗, 0, P1∗, 0), wherein H∗, Z∗, M1∗, and P1∗ are given by (10) and HI∗ is given by Substituting all equilibrium values of (10) and (12) in the first equation of (11), we obtain an equation with only the variable X ∗. Simplifying this equation, we obtain p(X ∗) = f (X ∗) −g(X ∗) = 0, where f (X ∗) and g(X ∗) are fourth-degree polynomials. The expressions are presented in Appendix E. By the Descartes’ rule of sign, we found that if 0 < Rc < 1, then there are zero or two positive roots (biologically viable) for X ∗. If Rc > 1 there is a unique positive root (biologically viable) for X ∗.
We analyze two cases to evaluate the dependence of the parameters γa and γX on the existence of positive equilibrium points for X ∗ with Rc ∈ [0, 1]. Firstly, we analyze the influence of normal platelets being added to infected red blood cells, and after that, the influence of normal platelets being added into merozoites. The analysis of both cases will be through bifurcation diagrams, where we fixed all parameters (Table 2), except β, γa, and γX, in order to explore the influence of γa and γX for different values of the infection parameter β.
Following the same idea, we analyze the influence of γX in the merozoites-presence equilibrium point. We found a similar behavior to the diagrams varying γa so that the diagrams will not be shown here.
Figure 1 shows that increasing the value of β increases the concentration of merozoite, which was expected. In addition, the higher β, the higher the value of γa needed to make the value of Rc smaller than one, that is, eliminating the parasite concentration. Here, we found a forward bifurcation, that is, when Rc < 1 there only is MFE point S0, which is LAS and if Rc > 1 a LAS MPE point S∗ appears and the MFE point S0 becomes unstable.
Appendix F illustrates a case with low platelet production. In that case, a backward bifurcation could appear. It is possible to have two MPE points when Rc is less than one. Also, depending on the β and γ values, oscillatory behaviors may arise.
3.2.2 Joint antiparasitic action of both normal and giant platelets
Here, we analyze the joint action of normal and giant platelets on combating merozoites infection. To evaluate the antiparasitic action of giant platelets, we take four cases for the action of giant platelets on infected red blood cells and merozoites so that the values of and increase: (i) and , (ii and , (iii) and , and (iv) and .
In order to investigate the effect of giant platelets to assist antiparasitic action of normal platelets, we consider the case of low production rate for megakaryocyte precursor cells; that is, we have a low platelet production. This case is presented in Appendix F.
The first case that will be portrayed is when only normal platelets are acting against merozoite infection. We chose a case from item (c) of Figure 3, with γa = 7.0 × 10−9, which has only unstable equilibrium points.
In Figure 2b, for a small value (blue color), the giant platelet action contributes to eliminating the oscillatory behavior of the dynamics of the system, stabilizing their values. As the antiparasitic action of giant platelets on infected red cells and merozoites increases, parasite concentration decreases, so in item (d), we have a convergence close to zero for X.
Thus, even in the case of low production of megakaryocyte precursor cells, with the antiparasitic action of giant platelets, it is possible to decrease the parasite concentration or even control it.
Therefore, giant platelets assisting the existing normal platelet antiparasitic action assist malaria infection control.
4 Discussion
We proposed a mathematical model to evaluate the antiparasitic action of both normal and giant platelets against merozoites and infected red blood cells during malaria infection. Our model showed that merozoites’ concentration when they suffered antiparasitic action was lower than the case without the antiparasitic response of platelets.
Without normal and giant platelet antiparasitic action, the model had two equilibrium points, the merozoite-free equilibrium S0, which was globally asymptotically stable if the threshold R0 is less than one, and unstable if R0 is greater than one. The merozoites-presence equilibrium S∗ was locally asymptotically stable if R0 > 1 and unstable if R0 < 1. The threshold R0 represents the net amount of new merozoites that one merozoite can produce. Therefore, if R0 < 1 even in a situation of colossal parasite inoculation, the host could fight against the merozoite and eliminate the infection. However, if R0 > 1, the introduction of only a single merozoite led to the persistence of parasitic concentration in the host.
When we included to the model the antiparasitic action of platelets, we obtained a new threshold, Rc, which was related to R0, with Rc ≤ R0. The threshold Rc represents the net number of new merozoites that one merozoite can produce under the effect of normal platelet antiparasitic action. Also, the antiparasitic action of normal platelets contributed to decreasing R0. Thus, in a case with R0 > 1, we can have a Rc < 1, leading to population dynamics converge to S0, but without antiparasitic action of platelets would converge to S∗. Moreover, in the threshold Rc, there were no giant platelet antiparasitic action parameters. Our model considered that the giant platelets were produced in the bone marrow only after the red blood cell infection by merozoites. In the existence of infected red blood cells, the megakaryocyte also produces this platelet type. Therefore, the antiparasitic action of normal platelets occurred as a primary response against the infection. The antiparasitic action of giant platelets occurred as a secondary response to assist the normal platelet antiparasitic action.
The antiparasitic action of normal platelets on infected erythrocytes, by means of the parameter γa proved to be able to control infection by itself, as long as γa is greater than the threshold Rc will be negative. Hence, no infection would persist even if we introduce a huge merozoite concentration in the bloodstream. However, this effect was not possible using only the antiparasitic action of normal platelets on merozoites, utilizing the parameter γX.
Considering a low platelet production, which may lead to a case of thrombocytopenia, a change of a forward to backward bifurcation as we increase the value of β occurred. In this case, we obtained the possibility of having up to three equilibrium points, S0, and two positive equilibrium points.
Depending on the parameters γa and γX, there was a possibility of appearing a limit cycle, which led to an oscillatory behavior for the population dynamics. The introduction of the giant platelet action to the case of a limit cycle could eliminate this oscillatory behavior. The higher the giant platelet’s action, the lower the concentration of the parasite.
When we considered a normal platelet production, the system (1) had only two equilibrium points. In that case, we had a classical behavior for the existence and stability of the equilibrium points; that is, we obtained only one positive equilibrium point, which was locally asymptotically stable if Rc is higher than one and unstable if Rc is less than one.
5 Conclusion
In conclusion, through a simple mathematical model, we determined the antiparasitic effect that normal and giant platelets can have on infected red blood cells and circulating merozoites. We emphasize that our goal with this model was not to have a quantitative framework but to present a theoretical approach to the dynamic behavior between antiparasitic action of platelets and the course of malaria infection to elucidate the role of platelets during malaria infection. In the absence of platelet antiparasitic activity, high levels of parasitic concentration were achieved, as expected. The antiparasitic action of normal platelets contributed to decreased parasite concentration levels in the bloodstream or even eliminated the parasite. However, the low production of platelets led to the occurrence of backward bifurcation in this model. It implied the existence of a merozoite-presence equilibrium point even when the threshold parameter for the net number of merozoites produced was below the unit, which showed that the malaria infection could persist or be controlled in a case of low platelet count. Also, in the case of low production of platelet, a limit cycle may appear. The appearance of giant platelets occurred to assist the antiparasitic action of normal platelets to eliminate the limit cycle and contributed to the control of the parasite. Models that take into account the immune system and antimalarial drugs merit further investigations.
Data Availability
The data used are available online.
Conflict of interest
The authors declare that they have no conflict of interest.
A Global stability of the equilibrium point T̅
Using the change of variables and , the system (2) results in a homogeneous linear system. The solution is given by where , and are eigenvectors associated with the eigenvalues λ1 = −(µZ + σ1), λ2 = −(µ1 + ε1), and λ3 = −µP1, respectively.
Notice that when t → ∞, the solution (13) tends to the null vector, so we concluded that . Therefore, the system converges to for any initial condition, resulting in a globally asymptotically stability.
B Local stability of the point W∗
To show the local stability of the point W ∗, we proved that the coefficients of the characteristic polynomial s(λ), given by equation (8), satisfy the Routh-Hurwitz conditions: a2 > 0, a1 > 0, a0 > 0 and a2a1 > a0. It’s easy to verify that a2, a1, and a0 are positives when R0 > 1. To prove the condition a2a1 > a0, we consider a2a1 − a0, and this term proved to be positive. rewriting, we found that The two last terms are canceled, so there are only positive terms. Thus, we concluded that a2a1 > a0.
C Proof of global stability of W0
Consider f a function defined by f : Ω1 → ℝ, such that and As W 0 = (k1/µH, 0, 0), observe that f (W 0) = 0 and f > 0 in Ω1 − {W 0}. In addition, Replacing the expressions for and we find that Therefore, in Ω1 −{W 0} if and only if R0 ≤ 1. It is concluded that merozoites-free equilibrium point W 0 is globally asymptotically stable.
D Proof of global stability of W∗
Consider the function g defined by g :, where Therefore, the derivative function is given by Simplifying the terms of the above equation, we have Being the arithmetic mean greater than or equal to the geometric mean, that is, we have and .
Thus, if we have ġ < 0 because Hence, we have ġ < 0 in . In addition, g(W ∗) = 0 and g > 0 in . Therefore, the merozoites-presence equilibrium, W ∗, is globally asymptotically stable.
E Polynomials f and g shown in section 3.2.1
F A case of normal platelet antiparasitic action with a low production rate for megakaryocyte precursor cells
Here, the megakaryocyte precursor cells have a low production rate; we consider K2 = 3.0 × 103. Through bifurcation diagrams varying the parameters γa and γX for different values of β, we analyze the antiparasitic action of normal platelets considering that there is a low platelet production.
we obtain that the graphs only have a reverse forward bifurcation for low β values, and the merozoitespresence equilibrium point is stable for 0 < Rc < 1.
In Figure 3, we obtain that the graphs only have a reverse forward bifurcation for low β values, and the merozoites-presence equilibrium point is stable for 0 < Rc < 1. However, as we increase β and high γa values an unstable region appears to the merozoites-presence equilibrium point. As we increase the value of β, the unstable region grows. Also, two positive equilibrium points may appear. In this case, the bifurcation diagrams have reverse backward bifurcation, which is given by Figures 4 to 5.
From Figures 3 and 4, we can see that as we increase the value of β, we pass from the existence of two equilibrium points to the presence of three equilibrium points. Depending on the value of γa, it is possible to have the following situations:
- there are two equilibrium points: a stable (non-zero equilibrium point) and the other is unstable (merozoites-free equilibrium point),
- there are only two unstable equilibrium points, according to the item (d) of Figure 3 (merozoites-free equilibrium point and a non-zero equilibrium point),
- there are three equilibrium points, one stable and two unstable, according to the item (b) of Figure 7 or two stable and one unstable equilibrium points, according to the item (b) of Figure 4.
In the case of three equilibrium point existence, we found that the stable region decreases when we increase β. Lastly, we can see that only to low values of γa, the merozoites-presence equilibrium point is stable; otherwise, it is unstable.
When γa surpasses a critical value , the merozoites-presence equilibrium point becomes unstable, and a cycle limit appears, which can be numerically assessed through the eigenvalues of the associated Jacobian matrix. We found a pair of complex eigenvalue with zero real part. For another value of γa, in this case, we denote by , we found a limit cycle, and when γa surpasses , the merozoites-presence equilibrium point becomes stable again.
As presented, a Hopf bifurcation appears in the diagrams since it presents a complex eigenvalue pair whose real part is null. To illustrate this situation, we chose item (c) of Figure 3. Figure 8 illustrates the value of (a) the real part and (b) imaginary part, which are assumed in the complex eigenvalue pair as we increase the value of γa.
We found an eigenvalue with null real part for and after for , showing Hopf bifurcations when varying γa. In the other bifurcation diagrams, it is also possible to find two Hopf bifurcations; they arise precisely at the point where the change of stability occurs.
Acknowledgements
The authors are grateful to Coordination of Superior Level Staff Improvement (CAPES) - Finance Code 001 by the financial support.
Footnotes
Tel.: +55-19-3521-6031, E-mail: hyunyang{at}ime.unicamp.br