<title>Introduction</title> <p>The parasites that cause human malaria are inoculated by an <italic>Anopheles</italic> mosquito and initially multiply in the liver. After about a week, a primary wave of merozoite forms enters the bloodstream, invades RBCs (within seconds), and continues the asexual cycle of multiplication, developing into the schizont forms that burst and release more merozoites. The pathology of malaria is due to this asexual blood- stage cycle (reviewed in ##UREF##0##[1]##,##REF##11832955##[2]##). The sexual forms transmissible to mosquitoes (gametocytes) appear over time, but in much smaller numbers than the asexual forms. The two parasite species that cause the vast majority of human cases, <italic>Plasmodium vivax</italic> and <italic>Plasmodium falciparum</italic>, can induce severe anemia; with <italic>P. vivax</italic> especially, the anemia can appear far out of proportion to the percentage of RBCs infected ##REF##12875288##[3]##.</p> <p>Both the innate and adaptive arms of the human immune system mount responses to infections with both species (reviewed in ##REF##15950450##[4]##,##REF##10899022##[5]##). High fevers are a classic feature of infections. During <italic>P. vivax</italic> infections near-periodic episodes of fever (“paroxysms”) are associated with high levels of tumor necrosis factor α (TNF-α) and other cytokines associated with innate immunity ##REF##1565611##[6]##–##UREF##1##[8]##. Strong cytokine responses also occur in <italic>P. falciparum</italic> infection, though the timing of paroxysms tends to be irregular ##REF##11254594##[9]##–##REF##16203971##[11]##. These fever paroxysms are associated with the synchronized release of merozoites from bursting schizonts. This synchronization has been the subject of considerable experimental and theoretical work. It is possible that febrile temperatures induce synchronization by differentially influencing development rates of different parasite stages ##REF##15275344##[12]## and that immune responses ##REF##11085243##[13]##,##REF##12621156##[14]## as well as the host's melatonin release cycle ##REF##10878815##[15]##–##REF##12682086##[17]## contribute to this phenomenon. But it is not yet clear whether synchronization helps parasites, perhaps in the way a sudden overwhelming abundance of prey may overwhelm a predator's capacity, or hinders them. Malaria parasites have certainly evolved mechanisms of immune evasion; however, <italic>P. falciparum</italic> exhibits antigenic variation ##REF##6187885##[18]##, adheres to vascular endothelium in response to fever ##REF##12177447##[19]##, and produces prostaglandins which probably modulate host TNF-α production ##REF##9743538##[20]##. This species also manages to keep the membrane of infected RBCs deformable during its ring stage (RBCs in the first 24 hours after infection), apparently reducing exposure of ring-stage parasites to clearance by the spleen ##REF##17517609##[21]##. <italic>P. vivax</italic> can also evade spleen clearance (reviewed in ##REF##15582531##[22]##) and suppress immune responses directed against its liver stage ##REF##16793328##[23]##.</p> <p>Clinical investigations suggest that the malaria parasite and host immune response interact with the host erythropoietic system in a complex, dynamic manner (reviewed in ##UREF##2##[24]##). Increased production of TNF-α by the host apparently induces anemia ##REF##3179231##[25]##; experimental evidence suggests that hemozoin produced by <italic>P. falciparum</italic> suppresses RBC production ##REF##16804108##[26]##. Abnormalities observed in <italic>P. falciparum</italic>-infected patients include suppression of erythroid progenitor cells in bone marrow ##REF##8348780##[27]##, sequestration of parasites in the marrow ##REF##3304391##[28]##, decreased iron absorption and hemoglobin synthesis ##REF##982520##[29]##, and decreased RBC survival time ##REF##3322362##[30]##. Leucocyte infiltration of the marrow and erythroblast degradation by macrophages have been observed in <italic>P. vivax</italic> infections ##REF##2660903##[31]##. Phagocytosis of uninfected RBCs has been observed in vitro and in <italic>P. berghei</italic>-infected mice ##REF##374106##[32]##, and is suspected in human malaria ##REF##10466119##[33]##,##REF##12614579##[34]##. <italic>P. chabaudi</italic>-infected mice show enhanced erythropoiesis that can compensate for RBC loss ##REF##12614579##[34]##; in humans, elevated levels of erythropoietin are produced in response to <italic>P. falciparum</italic> infection ##REF##7485716##[35]##. Overall, the evidence suggests that RBC destruction or ineffective erythropoiesis may thwart erythropoietin-initiated processes that might otherwise compensate for RBC loss, although erythropoietin may have other protective effects ##REF##18263734##[36]##.</p> <p>Most fatal malaria infections are due to <italic>P. falciparum</italic>, which can induce cerebral complications as well as severe anemia. <italic>P. vivax</italic> infections are characterized by lower levels of parasitemia, and, though often debilitating, are rarely fatal. <italic>P. falciparum</italic> attacks RBCs of all ages, while <italic>P. vivax</italic> mainly attacks reticulocytes (RBCs&lt;1.5 days old, still showing remnants of nuclei) ##UREF##3##[37]##–##UREF##5##[39]##, and possibly RBCs up to two weeks old ##REF##10450440##[40]##. In a previous report, we argued that targeted depletion of the youngest RBCs makes <italic>P. vivax</italic> infection potentially much more dangerous than is commonly appreciated: unchecked growth of a <italic>P. vivax</italic> population would eventually prevent the replacement of older, uninfected RBCs as these senesce and are culled from the circulation ##REF##15178766##[41]##. Thus one might expect a strong immune response to <italic>P. vivax</italic>, despite its seemingly lower threat relative to <italic>P. falciparum</italic>. Furthermore, in that model—which did not incorporate an immune response—compensatory boosting of RBC production tended to increase parasitemia, while dyserythropoietic response (reducing RBC production below its basal rate) had the opposite effect.</p> <p>Here we consider compartmentalized ordinary differential equations (CODEs) representing <italic>P. vivax</italic> and <italic>P. falciparum</italic> infections which incorporate a quick-acting innate response and a longer-term acquired antibody response as well as a dynamic erythropoietic system. ##FIG##0##Figure 1## shows the basic scheme; the details are presented below in the <xref ref-type="sec" rid="s4">Model</xref> section. The innate response emulates aspects of the fever paroxysm response which is the hallmark symptom of malaria. We analyze how these components jointly affect parasite and RBC dynamics. We do not attempt to model the full complexity of the immune response in malaria infections: our aim is to assess potential trade-offs between host and parasite for given characteristics of immune and erythropoietic responses. However, we do consider several choices of targets for both the innate and antibody responses in the model. We assume that bursting schizonts activate the innate response, and that the stage that triggers the antibody response is the same as the stage it targets. (The intraerythrocyte stage terminology is explained in the <xref ref-type="sec" rid="s4">Model</xref> section.) We also consider different values for the time constant of decay for the effective clearance rate of parasites by the antibody response. This time constant is not necessarily the biochemical decay constant of actual antibodies as it also incorporates the possibility of a long lived population of B-cells producing the antibodies. The models also incorporate a nonzero standard deviation in the time parasites develop within RBCs. We show that host erythopoietic response can affect infection outcome even in the presence of sustained immunological action. We show how infection outcome varies with the life stages of parasite targeted by the model immune responses. Furthermore, we show that a tight synchronization of merozoite release does not necessarily help parasite populations evade host immune responses. Most of our simulated infections assume that the host has no pre-existing antibodies or memory B cells for the parasite (and thus is “antibody-naïve”), but for some examples we examine the effects of pre-existing antibodies.</p>
<title>Results</title> <title>Typical Simulations</title> <p>A measurable quantity that constrains the model parameters is the initial reproduction rate of a parasite, <bold>R<sub>0</sub></bold>, defined as the average number of descendants an individual parasite would have at the beginning of bloodstream infection in the absence of any host response to the parasite. As long as <italic>T<sub>D</sub></italic>
<sub>μ</sub> (duration of the merozoite stage)≪the duration of the infected RBC stage, then by probability arguments independent of model details ##REF##15178766##[41]##,
</p> <p>Here <italic>V</italic>
<sub>0</sub> is the initial density of the RBC population vulnerable to the parasite species, <italic>ζ</italic> is the binding affinity of merozoites to RBCs, and <italic>p</italic> is the average number of merozoites released per bursting schizont. (For <italic>P. falciparum</italic>, <italic>V</italic>
<sub>0</sub> = normal basal RBC count = 5×10<sup>6</sup> µl<sup>−1</sup>; for <italic>P. vivax</italic>, <italic>V</italic>
<sub>0</sub> = the basal reticulocyte count = 6.25×10<sup>4</sup> µl<sup>−1</sup>.) Experimental evidence suggests that <italic>T<sub>D</sub></italic>
<sub>μ</sub> is just minutes; we take it as 0.1 h ##REF##6771738##[42]##. Data from the time series of parasitemia in neurosyphilis patients infected with <italic>P. vivax</italic> for malariatherapy ##REF##12099421##[43]## suggest that for <italic>P. vivax</italic>, <bold>R<sub>0</sub></bold>∼15. Similar data for <italic>P. falciparum</italic> from neurosyphilis ##REF##10432044##[44]## and other patients ##REF##9210658##[45]##, as well as direct photographic evidence from studies of bursting schizonts ##REF##16169475##[46]## suggest <bold>R<sub>0</sub></bold> and p are both ≥15 for that species as well. Knowing <bold>R<sub>0</sub></bold>, <italic>T</italic>
<sub>Dμ</sub>, <italic>p</italic>, and <italic>V</italic>
<sub>0</sub> determines <italic>ζ</italic>. Because we wanted to simulate infections with parasite growth rates seen clinically, we used <italic>p</italic> = 16 with <bold>R<sub>0</sub></bold> = 15 for both species. The characteristic invasion time of an RBC by a merozoite in this model, (ξ<italic>V</italic>
<sub>0</sub>)<sup>−1</sup>, is 24 seconds. Equation 1 is consistent with our CODEs in that if the immune response is absent, a steady state with a nonzero parasite count exists only if <bold>R<sub>0</sub></bold>&gt;1 ##REF##15178766##[41]##.</p> <p>We simulated infections for up to one year after primary release of merozoites (when they enter the blood from the liver), unless the parasite population was cleared by the host first, or unless the RBC count of the host declined to less than 60% of the normal basal count of 5×10<sup>6</sup> µl<sup>−1</sup> (in which case it is assumed that the host died). For each type of parasite and immune response, we considered four model RBC source dynamics, described in detail in the <xref ref-type="sec" rid="s4">Model</xref> section: (1) RBC production <italic>increases</italic> to match the RBC loss due to infection (compensatory response), (2) the RBC production rate strongly <italic>declines</italic> in response to RBC loss, down to 65% of the basal rate (severe dyserythropoiesis), (3) a more moderate dyserythropoiesis occurs, with RBC production rate reduced down to 75% of the basal rate, and (4) RBC production remains fixed at the basal rate. We cannot emulate the full complexity of erythopoietic response in malaria, but these choices represent plausible modalities in that response. Altogether, for each of the two species of parasites we conducted about 3.8×10<sup>4</sup> simulations with varying values of the erythropoietic and immune parameters and targets for an antibody-naïve host. In addition, for model <italic>P. vivax</italic> infections, we considered about 8×10<sup>3</sup> simulations in which a pre-existing level of antibody from a previous infection was present; details are in the <xref ref-type="sec" rid="s4">Model</xref> section.</p> <p>
##FIG##1##Figure 2## shows results from typical simulations with an antibody-naïve host. The parasitemia and immune factor time series shown in panels A–D compare 2 model <italic>P. vivax</italic> infections over the first 52 days after the primary release of parasites from the liver into the bloodstream. For both simulations, the model parameters for the parasites are the same, the rate of RBC production is fixed, and the immune responses have the same targets, but the binding affinities to those targets differ between the 2 models. In ##FIG##1##Figure 2A##, note that the lower peak count of IBCs (infected RBCs) occurs in the model with the larger, more sustained innate-response-driven clearance rate, but it is in the other model, with the larger antibody-driven clearance rate, that the infection is completely cleared during this time span. (Note: clearance rate = binding affinity to target×nominal level of the component that attacks the target; see <xref ref-type="sec" rid="s4">Model</xref> section.) ##FIG##1##Figure 2E## plots peak IBC count (<italic>I</italic>
<sub>PK</sub>) against integrated IBC count (<italic>I</italic>
<sub>INT</sub>, the total number of IBCs produced during the infection) for all of the <italic>P. vivax</italic> models with a fixed rate of RBC production. <italic>I</italic>
<sub>PK</sub> versus <italic>I</italic>
<sub>INT</sub> plots for the other parasite species and erythropoietic responses are similar (see ##SUPPL##0##Figure S1## for further examples). To indicate severity of anemia, we use the RBC deficit, ΔRBC, defined as 5×10<sup>6</sup> µl<sup>−1</sup> – (minimum RBC count during infection). ##FIG##1##Figure 2F## plots ΔRBC against I<sub>INT</sub> for all of the <italic>P. vivax</italic> simulations for an antibody-naïve host with a fixed rate of RBC production (black points) and moderate dyserythropoeisis (gray points). Note that for the same <italic>I</italic>
<sub>INT</sub>, ΔRBC tends to be (but is not always) one-half to one order of magnitude higher (i.e., anemia is much more severe) for the models with moderate dyserythropoeisis. Results for <italic>P. falciparum</italic> infections and for <italic>P. vivax</italic> infections with pre-existing antibodies are similar (see ##SUPPL##0##Figure S1##).</p> <title>Effects of Variations in Model Parameters on Simulation Outcomes</title> <p>Peak and integrated parasitemia in individual simulations depend on details of the interactions within the CODE model, but the plots in ##FIG##2##Figure 3## indicate the trends in simulation outcomes with changes in the immune response parameters in antibody-naïve hosts. Note in panel A that an order of magnitude change in the trigger level for the innate response leads to an order of magnitude change in <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub> on average. <italic>I</italic>
<sub>PK</sub> and especially <italic>I</italic>
<sub>INT</sub> are sensitive to changes in the binding affinity of the antibody response to its target (panels B and C), on average, but less so to the decay time of the attacker component of the antibody response (panel D). Note, too, in comparing panels B and C, that in order to achieve the same degree of control of parasitemia, the binding affinity of an antibody response targeting only merozoites must be many orders of magnitude greater than one targeting an IBC stage. Because the maximum nominal concentration of the attacker component is fixed at 10 µl<sup>−1</sup> for all antibody models (see <xref ref-type="sec" rid="s4">Model</xref> section), the corresponding clearance rates must also differ by many orders of magnitudes for the same degree of parasitemia control. This might be expected since free merozoites either invade an RBC or else are cleared within a few minutes, and thus are exposed to antibodies for a much shorter time than are IBCs. Though the model antibody response against merozoites works by removing them from the system, the model is also plausible for antibodies which block merozoites from infecting RBCs, since, as discussed above, merozoites are so quickly removed from the blood.</p> <title>Antibodies But No Innate Response</title> <p>In this section we consider simulated infections in which the host has no innate immune response; in the next sections we show that, as a class, model infections in which the host has an innate response have qualitatively different behavior. This contrast might be expected based on ##FIG##2##Figure 3A##, which shows that the peak parasitemia <italic>I</italic>
<sub>PK</sub> for simulations without an innate response (but with an antibody response) is orders of magnitude higher, on average, than for simulations with an innate response even for the smallest innate-parasite binding affinity and largest threshold of innate activation considered. An examination of infection outcomes without an innate response distinguishes outcomes to which it is likely to be most and least crucial, and is not an academic exercise: similar outcomes might occur if the parasite were able to evade or suppress innate response such that it became ineffective in eliminating parasites.</p> <p>
##FIG##3##Figure 4## shows the values of several quantities averaged over all simulated infections for each type of model antibody and erythropoietic response that we considered. The host lacks an innate response, and is antibody naïve. As noted above, results for antibody responses which attack the merozoite stage cannot be directly compared with those which attack IBC stages since antibodies against merozoites must induce a much higher clearance rate to be similarly effective. However, note that for all types of antibody responses, the peak and integrated parasitemia tend to be higher with compensatory erythropoiesis than with dyserythropoiesis or with RBC production fixed at the basal rate, strikingly so for <italic>P. vivax</italic> infections. Peak parasitemia <italic>I</italic>
<sub>PK</sub> in model <italic>P. falciparum</italic> infections is ∼10×<italic>I</italic>
<sub>PK</sub> in <italic>P. vivax</italic> infections (##FIG##3##Figure 4A##), as would be expected since a much smaller subset of the RBC population can be infected by the latter species. However, ##FIG##3##Figure 4B## shows that I<sub>INT</sub> in <italic>P. vivax</italic> infections is of the same order of magnitude as in <italic>P. falciparum</italic> infections. (The same is also true for the ΔRBC; see F##SUPPL##1##igure S2A##.) Furthermore, for a given combination of antibody and erythropoietic responses, the fraction of simulated infections which ends in host death by catastrophic anemia is nearly the same for the two species (##FIG##3##Figure 4C##), except that the fraction of infections with host antibodies to merozoites ending in death is lower for <italic>P. vivax</italic> than for <italic>P. falciparum</italic>. (Incidentally, for infections with antibodies to merozoites of either species, this fraction is lower with compensatory erythropoiesis than dyserthropoiesis.) But ##FIG##3##Figure 4A and 4B## show a major difference between the infection dynamics produced by the two species: the ratio <italic>I</italic>
<sub>PK</sub>/<italic>I</italic>
<sub>INT</sub> in model <italic>P. falciparum</italic> infections is just under 1, while in <italic>P. vivax</italic> infections it is closer to 1/20. This is the case because, if host death occurs, it occurs within 2–5 weeks after primary merozoite release in <italic>P. falciparum</italic> infection, but requires 8–15 weeks in <italic>P. vivax</italic> infection (see ##SUPPL##2##Figure S3##). Thus, much of the integrated parasitemia in <italic>P. falciparum</italic> infection in this set of models is due to a final burst of parasite production before host death. Another pattern is apparent in ##FIG##3##Figure 4## among the antibody responses which target an IBC stage: the longer the duration of vulnerability, the better the response is at controlling parasitemia, reducing risk of host death, and inducing parasite clearance in infections with either species. For example, note in ##FIG##3##Figure 4D## that 100% of infections with an antibody response that attacked only the bursting schizont stage (with no innate response present) ended in host death, but only about 30% of those with antibodies that attacked infected RBCs of all stages died. Furthermore, the latter, broad antibody response resolved nearly 45% of infections by clearance of the parasite within one year (##FIG##3##Figure 4D##). For model infections in which parasites are cleared, clearance tends to occur within 2–6 weeks after primary release; see ##SUPPL##2##Figure S3##.</p> <p>In the absence of an innate response, infection outcomes were rather insensitive to differences in the standard deviation of the development time of the parasite within the IBC (<italic>σ</italic>
<sub>IBC</sub>) (see ##FIG##4##Figure 5## and ##SUPPL##1##Figure S2B##). However, a shorter σ<sub>IBC</sub> was associated with a slightly higher peak and integrated parasitemia in model <italic>P. falciparum</italic> infections (##FIG##4##Figure 5A and 5B##). Peak parasitemia in <italic>P. vivax</italic> infection tended to be higher for a shorter <italic>σ</italic>
<sub>IBC</sub> as well, except with antibodies that target the late trophozoite-schizont stage (##FIG##4##Figure 5A##).</p> <p>As mentioned above, we also did simulations that assumed a level of pre-existing antibodies in the host. For this set of simulations, we assumed that the species was <italic>P. vivax</italic>, and that the antibodies could attack any IBC stage. (For more details see <xref ref-type="sec" rid="s4">Model</xref> section below.) We found that if there was no innate response, the results were very similar to those shown in ##FIG##3##Figure 4## for this type of antibody response; see ##SUPPL##3##Figure S4##.</p> <title>Simulations with an Innate Response Present</title> <p>
##FIG##5##Figure 6## is the same as ##FIG##3##Figure 4##, in that it shows the values of several quantities averaged over all simulated infections for each type of model antibody and erythropoietic response, except that the host also has an innate response. (As noted above, the model innate responses can attack either the early trophozoite stage, the late trophozoite/schizont stage, or the entire intraerythrocytic stage.) Hosts are antibody naïve. Many differences between the two figures are apparent. (Note the one order of magnitude difference in the vertical scales between ##FIG##3##Figures 4## and ##FIG##5##6## in panels A and B.) First, parasitemia, on average, is much better controlled by the host if an innate response is present: <italic>I</italic>
<sub>PK</sub> (panel A) is reduced a factor of 100–1000 and <italic>I</italic>
<sub>INT</sub> (panel B) by a factor 5–10. The presence of an innate response generally diminishes the effects of erythropoietic response on I<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub>, notably of compensatory response on <italic>P. vivax</italic>. The striking exception is that here severe diseryrthropoiesis with no antibody response, or an antibody response only against bursting schizonts, decreases <italic>I</italic>
<sub>INT</sub> for both species. On the other hand, the fraction of hosts surviving infections is sensitive to the nature of the erythropoietic response: only with dyserythropoiesis did infections with model innate responses end in host death (panel C). (This sensitivity to erythropoietic response is mirrored in ΔRBC; see ##SUPPL##4##Figure S5A##.) Dyserythropoiesis with a broad antibody response slightly increased the fraction of infections cleared (##FIG##5##Figure 6D##), but, even with severe dyserythropoiesis, for all model antibody types, the fraction of infections ending in host death was much lower if the host mounted an innate as well as an antibody response. As with ##FIG##3##Figure 4##, for antibody responses that target an IBC stage, the longer the duration of the vulnerable stage, the lower the peak and integrated parasitemia, the lower the fraction of infections ending in host death, and the higher the fraction of infections in which the parasite is cleared.</p> <p>
##FIG##6##Figure 7## again shows values of several quantities averaged over all simulated infections (in antibody naïve hosts with innate responses) for each type of innate response and erythropoietic response. (The averaging process includes averaging over all model antibody responses for the given innate and erythropoietic response.) For all three innate responses studied, the risk of host death is highly sensitive to the nature of the erythropoietic response. (So is ΔRBC; see ##SUPPL##4##Figure S5B##.) If the early trophozoite IBC stage is targeted by the innate response, <italic>I</italic>
<sub>PK</sub> tends to be lower than for the other two types of innate targets (##FIG##6##Figure 7A##). However, <italic>I</italic>
<sub>INT</sub> is lower if the late trophozoite-schizont stage is targeted (##FIG##6##Figure 7B##). The early trophozoite is a slightly better target than the late trophozoite-schizont in terms of reducing the risk of host death (##FIG##6##Figure 7C##) but slightly worse in terms of enhancing clearance (##FIG##6##Figure 7D##). The broad response, covering both IBC stages, is best at enhancing clearance but, with severe dyserythropoiesis, worst at reducing the risk of host death. Because the late trophozoites and schizonts are believed to be the parasite stages most vulnerable to fever ##REF##2642531##[47]##, we will return to these points in the <xref ref-type="sec" rid="s3">Discussion</xref>.</p> <p>Another striking aspect of both ##FIG##5##Figures 6## and ##FIG##6##7## is that infections with <italic>P. vivax</italic> and <italic>P. falciparum</italic> have very similar outcomes, on average, when the host can mount an innate attack on the parasite. We found that the differences between the species in the values of <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub> averaged over a given immune and erythropoietic response were sometimes as small as 1%. One major difference between the species is that, with a dyserythropoeitic response, the fraction of infections ending in host death is higher for <italic>P. vivax</italic> infections than for <italic>P. falciparum</italic> infections (##FIG##5##Figures 6C## and ##FIG##6##7C##). This result holds true, on average, for each type of antibody response and innate response. Apparently, since <italic>P. vivax</italic> attacks the youngest RBCs, even a low-level parasitemia can deplete the number of RBCs that mature, which in turn can more often or more quickly lead to catastrophic anemia if the rate of RBC production is reduced. On the other hand, in <italic>P. falciparum</italic> infections, low-level depletion of the general RBC pool by the parasite may be less critical. (Even for infections with compensatory erythropoiesis or a fixed rate of RBC production, in which the outcome tends to be much less severe for the host, ΔRBC is higher for hosts infected with <italic>P. vivax</italic> than with <italic>P. falciparum</italic>, for all types of immune responses; see ##SUPPL##4##Figure S5##). With a host innate response, with either parasite species and each type of antibody and erythropoietic response, host deaths that ended infections occurred at least 20 weeks after primary merozoite release, on average, much longer than if no innate response were present; see ##SUPPL##5##Figure S6##. With infections resolved by parasite clearance, for each antibody and erythropoietic response the average time from primary release to clearance was 1–6 weeks; see ##SUPPL##5##Figure S6##).</p> <p>Among the simulations of <italic>P. vivax</italic> infections in hosts with pre-existing antibodies to IBCs of all stages, we found essentially the same trends as in ##FIG##6##Figure 7##: for all model erythropoietic responses, the innate response which attacks the late trophozoites and schizonts is more effective, on average, than the other two in reducing <italic>I</italic>
<sub>INT</sub> and ΔRBC, while dyserthyropoiesis greatly aggravates ΔRBC and increases the frequency of host death, even if it also slightly increases the frequency of parasite clearance (see ##SUPPL##6##Figure S7##).</p> <title>Consistency of Development Time and Infection Outcome with an Innate Response</title> <p>As shown in ##FIG##4##Figure 5##, if no host innate response is present, parasitemia is either insensitive to <italic>σ</italic>
<sub>IBC</sub> or higher for smaller <italic>σ</italic>
<sub>IBC</sub>, depending on the parasite species and antibody response. The situation is strikingly different if an innate response is present. ##FIG##7##Figure 8## shows the values of several quantities from simulations with an innate response averaged over the type of antibody response and <italic>σ</italic>
<sub>IBC</sub>, and ##FIG##8##Figure 9## shows the values of the same quantities averaged over the type of innate response and <italic>σ</italic>
<sub>IBC</sub>. (In both figures, results are for antibody naïve hosts.) Note in these two figures that <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub> decline with decreasing σ<sub>IBC</sub>, on average, for each type of antibody response, for each type of innate response, and for both parasite species. (The same is true for the RBC deficit ΔRBC; see ##SUPPL##7##Figure S8##.) Apparently, the faster rate of schizont bursting which occurs with smaller <italic>σ</italic>
<sub>IBC</sub> triggers the innate response more quickly and keeps parasitemia lower, on average. The fraction of infections ending in host death also declines with decreasing <italic>σ</italic>
<sub>IBC</sub>, though not necessarily monotonically (##FIG##7##Figures 8C## and ##FIG##8##9C##). (For simulations that end in host death, the time from primary merozoite release to host death tends to increase with smaller <italic>σ</italic>
<sub>IBC</sub>, on average, though this effect is not monotonic; see ##SUPPL##8##Figure S9##.) The antibody response that attacks all IBC stages and the innate response that attacks late trophozoites and schizonts have, on average, lower <italic>I</italic>
<sub>INT</sub> and ΔRBC than the other model antibody and innate responses, respectively, for all values of <italic>σ</italic>
<sub>IBC</sub> (##FIG##7##Figures 8B## and ##FIG##8##9B##, and ##SUPPL##7##Figure S8##). The fraction of infections ending in parasite clearance has a somewhat complicated dependence on <italic>σ</italic>
<sub>IBC</sub>. If the innate response can attack the entire IBC stage, then this fraction increases with decreasing <italic>σ</italic>
<sub>IBC</sub>, but is insensitive to the value of <italic>σ</italic>
<sub>IBC</sub> for the other two model innate responses (##FIG##8##Figure 9D##). Interestingly, for some types of model antibody response, simulations with an innate response never lead to parasite clearance unless <italic>σ</italic>
<sub>IBC</sub>&lt;4.8 hours (##FIG##7##Figure 8D##).</p> <p>For the set of simulations done for <italic>P. vivax</italic> infections assuming a level of pre-existing antibodies to the entire IBC stage, the variations of <italic>I</italic>
<sub>PK</sub>, <italic>I</italic>
<sub>INT</sub>, ΔRBC, fraction of infections ending in host death, and fraction of infections ending in clearance of parasite with σ<sub>IBC</sub> were similar to the trends observed in the antibody naïve hosts (see ##SUPPL##9##Figure S10##).</p>
<title>Discussion</title> <p>Our results provide some insights into real malaria infections and parasite biology, despite the vast simplifications involved in modeling immune and erythropoietic responses. First, anemia did not lead to death in any simulation unless there was either dyserythropoiesis or an absence of innate response. Apparently, for <italic>P. vivax</italic> or <italic>P. falciparum</italic> infection to produce catastrophic anemia in the presence of a fast-activating innate response, one of the following must happen: (a) dyserythropoiesis, (b) a weakening of the innate response, coupled with an adaptive response which fails to check parasite growth, or (c) sustained destruction and/or sequestering of uninfected RBCs (“bystander effect”). We did not consider a bystander effect directly, but note that dyserythropoiesis could amplify the resulting anemia. We cannot say anything directly about host death due to other causes, but we note that infections with high parasitemia (<italic>I</italic>
<sub>PK</sub>&gt;10<sup>4</sup> µl<sup>−1</sup> or <italic>I</italic>
<sub>INT</sub>&gt;10<sup>6</sup> µl<sup>−1</sup>) were the ones with no or very weak innate response (as illustrated in ##FIG##2##Figures 3## and ##FIG##3##4##). Although control of the magnitude of peak parasitemia was dominated by the presence or absence of an innate response, our results suggest that an antibody response which targets a significant duration of the intraerythrocyte stage will reduce anemia and integrated parasitemia and boost the likelihood that a host will clear the infection. Even an antibody response that targets only merozoites would reduce parasitemia if the reciprocal of the clearance rate is greater than the characteristic time for a merozoite to find and invade an RBC (as suggested by ##FIG##2##Figure 3D##).</p> <p>It is hardly surprising to find that the antibody response with the widest possible range of targets produces the lowest frequency of host death and the highest frequency of parasite clearance (##FIG##3##Figures 4##, ##FIG##5##6##, and ##FIG##7##8##). Our results emphasize the different durations over which stage-specific antigens are exposed to antibody attack—an antibody response that targets merozoites or bursting schizonts is much less effective, for instance—and the value of an array of antibody responses to an array of stage-specific antigens. With an immune response that includes innate as well as antibody components, rather than antibody alone, the frequency of host death is lower (compare ##FIG##5##Figures 6C## and ##FIG##3##4C##), and, as with real malaria infections, more infections persist for long periods at low parasitemia.</p> <p>It is the gametocytes that allow parasite genomes to be transported between hosts, as modified by meiotic recombination in the mosquito, and the conditions that regulate gametocyte production and infectivity remain unknown. However, if the production or infectivity of gametocytes is proportional to either peak or integrated parasitemia, it appears unlikely that inducing dyserythropoiesis would provide a transmission advantage to either <italic>P. vivax</italic> or <italic>P. falciparum</italic>. Our calculations suggest that a dyserythropoietic response to infection, compared to compensatory erythropoiesis or RBC production held at the basal rate, would suppress I<sub>PK</sub> and I<sub>INT</sub>, most dramatically in hosts lacking an innate response to <italic>P. vivax</italic> (##FIG##3##Figure 4##), but also for I<sub>INT</sub> for both species, in hosts with both immune-response components engaged (##FIG##5##Figure 6##). In the latter case, though dyserthropoiesis can slightly increase the frequency of parasite clearance, it can substantially increase the risk of severe anemia and host death. Our model does not consider the possibility that induction of dyserythropoiesis might be accompanied by reduced effectiveness of immune responses, however.</p> <p>It is possible that enhanced parasitemia enhances propagation of the parasite's genome, presumably by boosting gametocytemia, but it is also possible that enhanced parasitemia or gametocytemia triggers an immune response that limits transmission, and thus that the <italic>regulation</italic> of parasitemia or gametocytemia could benefit the parasite in an evolutionary sense. In vitro experiments have shown that typical febrile temperatures kill intraerythrocytic parasites in the last but not the first 24 hours of their development ##REF##2642531##[47]##. Interestingly, in our CODE models, I<sub>INT</sub> (and thus the total number of merozoites produced) is lower on average when the innate response attacks the late trophozoites and schizonts rather than the early trophozoites, a result that is a consequence of the model <italic>dynamics</italic>. If there are parasite adaptations that involve direct or indirect regulation of parasitemia that enhance gametocyte transmission, then the parasite may actually benefit by the vulnerability of the second half of the IBC stage to fever.</p> <p>When an innate response is present in the model, compensatory erythropoietic response leads to nearly the same <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub>, on average, as RBC production fixed at the basal rate. Our results show slightly lower anemia with a compensatory response than with a fixed RBC production rate, in the presence of innate immunity (##SUPPL##4##Figure S5##), but differences between the two in any outcome are remarkably few and small, which implies that, in the presence of an innate immune response, a compensatory erythropoietic response confers little or no advantage or disadvantage to either host or parasite. This marks a contrast to some of our results with an earlier model ##REF##15178766##[41]## in which, specifically to examine constraints on immune responses, we included only parasite and RBC dynamics. In that context, compensatory response boosted parasitemia and accelerated the onset of catastrophic anemia in a set of model <italic>P. vivax</italic> infections. Our conclusion here might change if there were long-term evasion or suppression of the innate response by the parasite—since, in the absence of an innate response, compensatory response boosts parasitemia (##FIG##3##Figure 4##)—or if removal of uninfected RBCs boosted RBC production well above rates that would compensate for infection-induced RBC loss ##REF##374106##[32]##. Overall, our results on the role of erythropoietic response point to larger questions of how trade-offs of compensatory and dyserythropoietic responses might balance, and how differences between individuals translate into differences in populations; the evolutionary aspects involve dynamics on multiple scales, to be incorporated in future work ##REF##15572065##[48]##.</p> <p>Our results showed lower <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub> for smaller values of the standard deviation in the IBC development time, <italic>σ</italic>
<sub>IBC</sub>, provided an innate response was present (see ##FIG##7##Figures 8## and ##FIG##8##9##). This was true for all innate and antibody targets and for both species, suggesting that a very tight synchronization of merozoite release by bursting schizonts would not benefit the parasite if parasitemia correlates with transmission probability. Furthermore, tight synchronization benefits the host overall, through enhanced clearing of infections, and reduced RBC deficits and frequencies of death. It has been suggested, however, that the release of prostaglandins or other factors by bursting schizonts during synchronization may suppress immunity ##REF##9743538##[20]##, an effect we did not consider. Furthermore, we assumed that the model immune responses were not overwhelmed by large increases in parasitemia, which should be true if the response is provided by fever directly, or by a large reservoir of antibody-producing cells. However, a limited pool of cells that target and kill parasites could perhaps be overwhelmed by a synchronized burst. And we did find that if the innate response is absent, parasite synchronization yields little benefit for the host. Again, these results indicate trade-offs that merit further investigation ##REF##12621156##[14]##.</p> <p>Finally, we address the relative abilities of <italic>P. vivax</italic> and <italic>P. falciparum</italic> to cause catastrophic illness. In our previous report—on a model which incorporated erythropoeitic but not immune responses—we inferred that without an effective immune response, <italic>P. vivax</italic> would be highly lethal, even for values of <bold>R<sub>0</sub></bold> just above 1, due to its depletion of reticulocytes and the consequent eventual lack of replacements for senescent RBCs ##REF##15178766##[41]##. Among the results of our new simulations, done with active immune responses, two trends stood out when comparing the outcome of simulated <italic>P. vivax</italic> infections with those caused by <italic>P. falciparum</italic>: (1) for the same host immune and erythropoietic responses, <italic>P. vivax</italic> produces ΔRBC values equal to or even greater than those of <italic>P. falciparum</italic>, and (2) under similar conditions of either dyserythropoiesis or an ineffective innate response, the fraction of simulations ending in catastrophic anemia was just as high or higher for <italic>P. vivax</italic> as for <italic>P. falciparum</italic>. In fact, as can be seen in ##FIG##5##Figures 6C## and ##FIG##6##7C##, no <italic>P. falciparum</italic> infections with both innate response and moderate dyserythropoiesis ended in host death, while a small fraction of <italic>P. vivax</italic> infections with those characteristics did. Though real <italic>P. vivax</italic> infections are rarely fatal, they are debilitating and trigger extremely strong immune responses (such as the paroxysm response described above). Our new results suggest that if these responses truly clear the parasite from the blood, their strength is not inappropriately aggressive, and may be a result of strong selective pressures from <italic>P. vivax</italic> malaria on human innate and adaptive immune systems. Indeed, recent reports of <italic>P. vivax</italic> malaria with severe organ complications ##REF##15705338##[49]##, splenic infarction ##REF##17374979##[50]##, and retinal hemorrhage ##REF##15099997##[51]## suggest that <italic>P. vivax</italic> can induce “malignant” illness in some individuals. There is also an increased risk of low birthweight in infants born to mothers with <italic>P. vivax</italic> malaria ##REF##10470698##[52]##.</p> <p>In both evolutionary and proximate terms, immune responses shaped by <italic>P. vivax</italic> infections are likely to be partly but not completely appropriate for <italic>P. falciparum</italic> infections, and vice versa ##REF##11121847##[53]##,##REF##16837717##[54]##. It may be that <italic>P. falciparum</italic> tends to cause more severe illness often than <italic>P. vivax</italic> because it is more effective at evading, suppressing or otherwise subverting the human immune responses on a population average. These is also the possibility that dyserythropoiesis is worse in <italic>P. falciparum</italic> infection. How these and other factors play off in determining the relative transmission success of the two species over multiple time scales is an intriguing topic for future investigation.</p>
<p>Conceived and designed the experiments: PGM FEM. Performed the experiments: PGM. Analyzed the data: PGM FEM. Contributed reagents/materials/analysis tools: PGM. Wrote the paper: PGM FEM. Performed numerical simulations and contributed the main analysis software used to analyze the simulation results: PGM.</p> <p>The two main agents of human malaria, <italic>Plasmodium vivax</italic> and <italic>Plasmodium falciparum</italic>, can induce severe anemia and provoke strong, complex immune reactions. Which dynamical behaviors of host immune and erythropoietic responses would foster control of infection, and which would lead to runaway parasitemia and/or severe anemia? To answer these questions, we developed differential equation models of interacting parasite and red blood cell (RBC) populations modulated by host immune and erythropoietic responses. The model immune responses incorporate both a rapidly responding innate component and a slower-responding, long-term antibody component, with several parasite developmental stages considered as targets for each type of immune response. We found that simulated infections with the highest parasitemia tended to be those with ineffective innate immunity even if antibodies were present. We also compared infections with dyserythropoiesis (reduced RBC production during infection) to those with compensatory erythropoiesis (boosted RBC production) or a fixed basal RBC production rate. Dyserythropoiesis tended to reduce parasitemia slightly but at a cost to the host of aggravating anemia. On the other hand, compensatory erythropoiesis tended to reduce the severity of anemia but with enhanced parasitemia if the innate response was ineffective. For both parasite species, sharp transitions between the schizont and the merozoite stages of development (i.e., with standard deviation in intra-RBC development time ≤2.4 h) were associated with lower parasitemia and less severe anemia. Thus tight synchronization in asexual parasite development might help control parasitemia. Finally, our simulations suggest that <italic>P. vivax</italic> can induce severe anemia as readily as <italic>P. falciparum</italic> for the same type of immune response, though <italic>P. vivax</italic> attacks a much smaller subset of RBCs. Since most <italic>P. vivax</italic> infections are nonlethal (if debilitating) clinically, this suggests that <italic>P. falciparum</italic> adaptations for countering or evading immune responses are more effective than those of <italic>P. vivax</italic>.</p> <title>Author Summary</title> <p>Of the four <italic>Plasmodium</italic> species that parasitize human erythrocytes and induce malaria, <italic>Plasmodium falciparum</italic> and <italic>Plasmodium vivax</italic> cause most of the public health burden. <italic>P. falciparum</italic> infection is typically characterized as “malignant” (due to severe, sometimes lethal consequences, particularly in immune-naïve individuals), and <italic>P. vivax</italic> malaria as (relatively) “benign.” Using the power of a Beowulf cluster, we tested hypotheses about host control of malaria by simulating ∼8.4×10<sup>4</sup> combinations of parasite species, host immune response, and erythropoietic response to infection. We tailored the models to specific details of the life cycles of the two species, which invade different subclasses of red blood cells. Our results challenge some standard assumptions. For example, we show that tight synchronization of the asexual reproduction of malaria parasites may actually benefit the host by reducing parasitemia. We also demonstrate that properties of host immunity or erythropoiesis that contribute to high parasitemia and severe anemia in <italic>P. falciparum</italic> malaria would do so in <italic>P. vivax</italic> infection as well, in line with recent reports indicating that <italic>P. vivax</italic> can indeed cause “malignant” illness in some patients. This suggests that <italic>P. falciparum</italic> is more effective overall at immune evasion or suppression than <italic>P. vivax</italic>.</p>
<title>Model</title> <title>General Considerations</title> <p>Our models are composed of compartmentalized ordinary differential equations (CODEs) and are closely similar to models in our previous work ##REF##15178766##[41]##,##REF##16837717##[54]## but include immune responses. CODEs are designed to simulate an ensemble of cells that, on average, take a certain duration <italic>T</italic>
<sub>D</sub> to pass through a development stage. The standard deviation in development time iswhere <italic>N</italic>
<sub>c</sub> is the number of compartments used to model the development stage ##UREF##6##[55]##,##UREF##7##[56]##. We stress that <italic>T</italic>
<sub>D</sub> and <italic>σ</italic> are the tangible quantities.</p> <p>
##FIG##0##Figure 1## is a schematic representation of the model. RBCs emerge from the bone marrow, and, if uninfected, circulate for 120 days ##UREF##8##[57]##. We assume that RBCs of all age cohorts are vulnerable to <italic>P. falciparum</italic>, but that <italic>P. vivax</italic> attacks only reticulocytes. The rate of RBC loss due to infection can feed back to the RBC source, to increase (compensatory response) or decrease (dyserythropoietic response) the RBC production rate. The immune components of our model are based on simple ideas from immunology, and resemble models of immune response in both viral infections ##REF##11602708##[58]##–##REF##12447982##[60]## and malaria ##REF##12614579##[34]##. An antigen signal over a threshold level triggers an initially self-amplifying (and eventually self-limiting) <italic>actuator</italic> phase, which in turns initiates either an <italic>attacker</italic> phase, which directly removes parasites, or a <italic>growth</italic> phase, which models the division of B or T cells, in which case the product of the growth phase initiates the attacker phase. The two-stage actuator-attacker response emulates a quick-acting but short-lived innate response, while the three-stage actuator-growth-attacker response emulates a slower-acting adaptive response, such as the production of antibodies against a target. As real antibody responses have control mechanisms that keep them from reaching destructive levels (reviewed in ##UREF##9##[61]##), in our three-stage model immune response the attacker phase is also self-limiting.</p> <p>We now introduce the equations of motion for the system, starting with those for the parasite population. Then we consider those for the RBC population and the immune factors.</p> <title>Parasite Population Dynamics</title> <p>Let <italic>T</italic>
<sub>DI</sub> be the average duration of the infected RBC stage, with standard deviation <italic>σ</italic>
<sub>IBC</sub>. Then <italic>N</italic>
<sub>cI</sub> = <italic>T</italic>
<sub>DI</sub>
<sup>2</sup>/<italic>σ</italic>
<sub>IBC</sub>. is the number of compartments needed to described the IBC development. Let <italic>I<sub>n</sub></italic> be number of infected RBCs per µl in compartment <italic>n</italic>. Then
</p> <p>Here κ<sub>I</sub> = <italic>N</italic>
<sub>cI</sub>/<italic>T</italic>
<sub>DI</sub>, <italic>μ</italic> is the merozoite density, <italic>ζ</italic> is the binding affinity between the merozoites and vulnerable RBCs, and <italic>V</italic> is the total density of vulnerable RBCs for a given species: for <italic>P. vivax</italic>, <italic>V</italic> = uninfected reticulocyte count, and for <italic>P. falciparum</italic>, <italic>V</italic> = <italic>E</italic>
<sub>T</sub>, the total uninfected RBC count. The infected RBCs (IBCs) could themselves be vulnerable to several immune system components: <italic>Att<sub>m</sub></italic> refers to the density of the attacker phases of component <italic>m</italic>, and <italic>ξ<sub>m</sub></italic>
<sub>,<italic>n</italic></sub> is the binding affinity between the attacker phase of immune component m and the parasite in compartment <italic>n</italic>. We allow <italic>ξ</italic> to have dependence on <italic>n</italic> in order to emulate an age-specific immune attack on IBCs. This is valid within the CODEs formalism provided one realizes that if <italic>N</italic>
<sub>sub</sub> is the number of compartments associated with a subinterval of the IBC stage, the standard deviation in the time any one IBC spends in the subinterval is ∼<italic>T</italic>
<sub>DI</sub>
<italic>N</italic>
<sub>sub</sub>
<sup>1/2</sup>/<italic>N</italic>
<sub>cI</sub>. We took <italic>T</italic>
<sub>DI</sub> = 48 h for both <italic>P. vivax</italic> and <italic>P. falciparum</italic>
##UREF##0##[1]##. For each immune and erythropoietic response, we evaluated the model's behavior for several values of <italic>σ</italic>
<sub>IBC</sub> (given in ##TAB##0##Table 1##). The discussion after Equation 1 gives the details on how the value of <italic>ζ</italic> is set.</p> <p>We use just one compartment for the merozoite stage, due to its short survival time in the blood, <italic>T</italic>
<sub>Dμ</sub> = <italic>σ</italic>
<sub>μ</sub> = 0.1 h ##REF##6771738##[42]##:
</p> <p>Here, <italic>p</italic> is the number of merozoites released per bursting IBC after development and asexual division are completed within an IBC; the choices for <italic>p</italic> are discussed below. Parameter <italic>ξ<sub>m</sub></italic>
<sub>,<italic>μ</italic></sub> is the binding affinity of immune response <italic>Att<sub>m</sub></italic> to merozoites. <italic>L</italic>(<italic>t</italic>) is the infusion of primary merozoites of the given species from the liver into blood, which is believed to happen quickly and to involve 10<sup>4</sup>–10<sup>5</sup> merozoites. For our simulations, the initial time (<italic>t</italic> = 0) corresponds to the release of the first parasite from the liver. For simplicity, we took <italic>L</italic>(<italic>t</italic>) = 0.002 (µl h)<sup>−1</sup>, for the first hour, (corresponding to a release of 10<sup>4</sup> primary merozoites into a blood volume of 5×10<sup>6</sup> µl) and then zero afterwards. We redid some simulations with 10<sup>5</sup> primary merozoites released, but this had little effect on the outcome (results not shown).</p> <p>In <italic>P. falciparum</italic> malaria, the IBCs sequester onto blood vessel walls during the second phase (last 24 hours) of parasite intra-RBC development. Sequestered parasites can readily release large numbers of merozoites that attack neighboring RBCs, however ##REF##16169475##[46]##. Thus, in our model, we can consider the IBC compartments to include sequestered as well as freely-circulating IBCs.</p> <title>RBC Population Dynamics</title> <p>The RBC development chain is divided into three parts. We first consider RBCs attacked by <italic>P. falciparum</italic>; the CODEs for reticulocytes areFor mature RBCs:For senescent RBCs:(We separate the senescent stage from the mature stage to allow for possible modeling of <italic>P. malariae</italic> infections, in which mainly senescent RBCs are attacked.) Here, <italic>κ</italic>
<sub>R</sub> = <italic>N</italic>
<sub>c<italic>R</italic></sub>/<italic>T</italic>
<sub>DR</sub>, <italic>κ</italic>
<sub>M</sub> = <italic>N</italic>
<sub>c<italic>M</italic></sub>/<italic>T</italic>
<sub>DM</sub>, and <italic>κ</italic>
<sub>S</sub> = <italic>N</italic>
<sub>c<italic>S</italic></sub>/<italic>T</italic>
<sub>DS</sub>, where the <italic>T</italic>
<sub>DR</sub>, <italic>T</italic>
<sub>DM</sub>, and <italic>T</italic>
<sub>DS</sub> are the durations of the respective RBC stages, and <italic>N</italic>
<sub>c<italic>R</italic></sub>, <italic>N</italic>
<sub>c<italic>M</italic></sub>, and <italic>N</italic>
<sub>c<italic>S</italic></sub> are the respective numbers of compartments used for each stage (as set by Equation 2 above). ES(<italic>t</italic>) is the rate of new RBC production at time <italic>t</italic>. Based on physiologically reasonable values ##UREF##8##[57]##, we took <italic>T</italic>
<sub>DR</sub> = 36 h with <italic>σ</italic> = 6 h, <italic>T</italic>
<sub>DM</sub> = 2,796 h with <italic>σ</italic> = 168 h, and <italic>T</italic>
<sub>DS</sub> = 48 h with σ = 12 h. For <italic>P. vivax</italic> infections, the model for RBC development is exactly the same, except that the mature and senescent stages are not attacked at all by merozoites. The total uninfected RBC count, <italic>E</italic>
<sub>T</sub>, is equal to the sum of all the <italic>R</italic>, <italic>M</italic>, and <italic>S</italic> compartments. If <italic>E</italic>
<sub>T</sub> drops below 3×10<sup>6</sup> RBCs µl<sup>−1</sup> (i.e., 0.6×the basal count of a typical healthy adult), then we assume that the host dies of catastrophic anemia. Studies of RBC or hemoglobin levels in patients with either <italic>P. falciparum</italic>
##REF##982520##[29]##,##REF##17060471##[62]## or <italic>P vivax</italic>
##REF##12875288##[3]## infections suggest that RBC counts can collapse to similarly low fractions of the basal count.</p> <title>Erythropoietic Response</title> <p>The dynamics for the marrow RBC source depend on the host response to losses of uninfected RBCs. Let Φ = ES<sub>0</sub>−d<italic>E</italic>
<sub>T</sub>/d<italic>t</italic>−<italic>δζμV</italic>, where ES<sub>0</sub> is the basal rate of RBC production (which maintains a healthy basal count of 5×10<sup>6</sup> µl<sup>−1</sup>), and <italic>ζ</italic>, <italic>μ</italic>, and <italic>V</italic> are as above in Equations 3 and 4. The factor <italic>δ</italic> is a simplistic way to account for the dyserythropoietic effects during infection. We model the dynamics of ES(<italic>t</italic>) with the following ODE:
</p> <p>ES<sub>MX</sub> is the maximum RBC production rate, and 1/<italic>λ</italic>
<sub>ES</sub> is a response time to changes in the rate of RBC loss. For a healthy compensatory response to RBC loss ##UREF##8##[57]##, ES<sub>MX</sub>/ES<sub>0</sub> = 5 and 1/<italic>λ</italic>
<sub>ES</sub> = 48 h. ES<sub>MN</sub> is a minimum production rate; with dyserythropoiesis, the RBC rate production would be driven to ES<sub>MN</sub>. As mentioned above, we considered four model RBC source dynamics, the parameters specified in ##TAB##1##Table 2##.</p> <title>Immune Response Dynamics</title> <p>As noted above, we consider two main types of immune actions: a fast-activating but short-acting innate response, and a slower-activating but long-acting antibody-like response. We consider the innate response first. Let <italic>Act</italic> be the actuator level, <italic>Att</italic> the attacker level, and <italic>A</italic> an antigen (targeted parasite stage) that triggers the response. Actuator production above its basal rate is triggered when the density of a target parasite stage <italic>A</italic> exceeds a detection threshold <italic>A</italic>
<sub>th</sub>:where Θ(<italic>x</italic>) = <italic>x</italic> if <italic>x</italic>&gt;0, zero otherwise, and
</p> <p>
<italic>S</italic>
<sub>Act,B</sub> is the basal actuator level of production, <italic>a</italic> the self-amplification, and <italic>Act</italic>
<sub>Mx</sub> the maximum allowed actuator level. Parameter <italic>S</italic>
<sub>Att,th</sub> is fixed by the other parameters so that production of the attacker does not begin unless the actuator production is above its basal level. Since schizont rupture is believed to trigger the short-term innate reaction associated with paroxysms in <italic>P. vivax</italic> infection ##REF##12689650##[7]##, we took the actuator phase to be sensitive to the count of IBCs in the last 4% of their development cycle (“bursting schizonts”). Clinical studies have shown rapid growth and quick decay of TNF-α levels in <italic>P. vivax</italic> infections ##REF##1565611##[6]##, on the order of 1–2 h. Fever spikes follow the TNF-α peaks and also decay quickly in such patients, on the order of 2–3 h. For the simulations reported here, we assumed 1/<italic>λ</italic>
<sub>Act</sub> = 1 h, <italic>Act</italic>
<sub>Mx</sub> = 2.5 µl<sup>−1</sup>, with <italic>S</italic>
<sub>Act,B</sub> set so that the basal actuator level is 0.1, <italic>a</italic> = 100, and 1/<italic>λ</italic>
<sub>Att</sub> = 2 h. Test simulations indicated that model results were rather insensitive to values of these parameters as long as <italic>λ</italic>
<sub>Act</sub> and <italic>λ</italic>
<sub>Att</sub> were ∼1 h<sup>−1</sup> (results not shown).</p> <p>Simulation outcomes were sensitive to the antigen detection threshold <italic>A</italic>
<sub>th</sub> (see Equation 9) and to <italic>ξ</italic>, the binding affinity to the target of the response (as defined in Equations 3 and 4). In our simulations, we considered three values of <italic>A</italic>
<sub>th</sub> (given in ##TAB##2##Table 3##). We also considered three targets for this response: (1) IBCs &lt;24 h old, (early trophozoite), (2) IBCs between 24 and 48 h old (late trophozoite and schizont), and (3) IBCs of any age. For each target and sensitivity, we considered three values of binding affinity <italic>ξ</italic> (given in ##TAB##2##Table 3##). In a sense, the product <italic>ξ</italic>
<italic>Att</italic> is the tangible quantity, as this is the actual clearance rate of the subpopulation targeted by this response.</p> <p>For the antibody-like response (Ab response), the actuator triggers an intermediate growth phase which imposes a delay on the production of the attacker; thus, we use a multicompartment model for it. Let <italic>G<sub>n</sub></italic> be the amount of intermediate phase in compartment <italic>n</italic>, and let <italic>T</italic>
<sub>DG</sub> be the duration of the growth phase with standard deviation <italic>S</italic>
<sub>DG</sub>:where
</p> <p>Parameters <italic>S</italic>
<sub>Act,B</sub>, <italic>a</italic>, <italic>Act</italic>
<sub>Mx</sub>, <italic>λ</italic>
<sub>Act</sub>, <italic>A</italic>
<sub>th</sub> and <italic>λ</italic>
<sub>Att</sub> have the same roles as for the innate response, <italic>N</italic>
<sub>c<italic>G</italic></sub> is set by Equation 2 above, κ<sub>G</sub> = <italic>NcG/T<sub>DG</sub></italic>, and <italic>Att</italic>
<sub>MX</sub> is the maximum allowed level of the attacker. S<sub>Gth</sub> is set so that there is no production of the growth phase unless the actuator is above its basal level. Note, too, that the definition for <italic>S</italic>
<sub>Act</sub> is slightly different from the one in equation 9, as we allow for a shutdown of actuator production once the attacker reaches its maximum level: all production of the attacker shuts down once the attacker level reaches <italic>Att</italic>
<sub>Mx</sub>. Results with mouse models suggest that the first specific activated B or T cells appear 4–7 days after inoculation with blood-stage parasites ##REF##17084842##[63]##,##REF##15987340##[64]##, so we take <italic>T</italic>
<sub>DG</sub> = 96 hours, with <italic>σ</italic> = 9.6 hours. As for the innate response, we assume 1/<italic>λ</italic>
<sub>Act</sub> = 1 h, <italic>Act</italic>
<sub>Mx</sub> = 2.5 µl<sup>−1</sup> and <italic>S</italic>
<sub>Act,B</sub> set so that the basal actuator level is 0.1. However, we took the self-amplification factor <italic>a</italic> to be much smaller, 0.01. We set <italic>Att</italic>
<sub>Mx</sub> = 10 µl<sup>−1</sup>. Again, the product <italic>ξAtt</italic> is the tangible quantity, not <italic>Att</italic> or <italic>ξ</italic> directly.</p> <p>For the Ab response, we assume that the targeted subpopulations of parasites also provide the trigger for the actuator. We considered five possible targets: (1) early trophozoite, (2) late trophozoite and schizont, (3) bursting schizonts (as defined above), (4) IBCs of any age, and (5) merozoites. Test simulations indicated that if merozoites were the target of the model antibody response, then the binding affinity <italic>ξ</italic> had to be many orders of magnitude larger and the trigger threshold <italic>A</italic>
<sub>th</sub> much smaller than if the other subpopulations were the target in order to achieve the same degree of control on parasitemia. (This is in line with other studies ##REF##11085243##[13]##,##REF##12614579##[34]##; see <xref ref-type="sec" rid="s3">Results</xref> section above.) The values used for <italic>A</italic>
<sub>th</sub> and <italic>ξ</italic> for each of the targeted parasite populations are given in ##TAB##3##Table 4##. For all choices of target population, <italic>A</italic>
<sub>th</sub> and <italic>ξ</italic>, we considered three values for the time decay constant of the attacker phase, <italic>λ</italic>
<sub>Att</sub> (also listed in ##TAB##3##Table 4##). For simplicity, we do not directly consider the effects of antigenic variation in <italic>P. falciparum</italic> infection. However, in order to consider the effects of an ineffective adaptive response, we did simulations with an innate response but no antibody response, and, symmetrically, simulations with an antibody response but no innate response.</p> <p>For most of the simulations, we started with an antibody-naïve host; i.e., the initial level of attacker component is zero. However, for all model innate and erythropoietic responses to <italic>P. vivax</italic> infection, we did addition simulations with pre-existing antibodies of the response that targets IBCs of any age. The parameters <italic>ξ</italic>, <italic>A</italic>
<sub>th</sub>, and <italic>λ</italic>
<sub>Att</sub> used for this set of simulations have the same values as those listed for the “All IBCs” model antibody response in ##TAB##3##Table 4##. As mentioned above, the initial antibody level is set at <italic>e</italic>
<sup>−1</sup>×<italic>Att</italic>
<sub>Mx</sub>, or ∼3.679 µl<sup>−1</sup>.</p> <title>Simulation Strategy</title> <p>
##TAB##0##Tables 1##–

##TAB##3##4## summarize the values used for all of the model parameters which we varied from simulation to simulation. A simulation was done for each possible choice of parameters in the tables, thus results of 38,080 simulations for antibody-naïve hosts are reported for each of the two <italic>Plasmodium</italic> species. In addition, we did a further 8,064 simulations of <italic>P. vivax</italic> infection with pre-existing antibodies. Model infections were simulated for up to one year after the primary release of parasites from the liver into the bloodstream, unless the host died from catastrophic anemia or the parasite was cleared from the host before a full year elapsed. Although the ODE system has hundreds to thousands of dependent compartmental variables to solve, it is not stiff in the parameter regimes of interest, and thus can be solved efficiently. Of course, we still only examined a small part of the parameter space of the model, and in many ways our model is a caricature of the real parasite-host system. Nonetheless, in sweeping across these many combinations of parameter values we identified several basic, broad trends with changes in the structures of the immune and erythropoietic responses.</p> <p>If at any point the merozoite count, the total infected RBC count, or the total uninfected RBC count dropped to &lt;1 in a total blood volume of 5×10<sup>6</sup> µl, the values of all compartments that contributed to that particular count were reset to zero. The CODE system was solved using the fifth-order Runge-Kutta-Fehlberg algorithm with adaptive stepsize control for time integration ##UREF##10##[65]##,##UREF##11##[66]##, so that the difference between the fourth- and fifth-order solutions for each component of the ODE systems was less than one part in 10<sup>6</sup>.</p> <title>Supporting Information</title>
<p>This study utilized the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health, Bethesda, Md. (<ext-link ext-link-type="uri" xlink:href="http://biowulf.nih.gov">http://biowulf.nih.gov</ext-link>).</p>
<fig id="pcbi-1000149-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g001</object-id><label>Figure 1</label><caption><title>Schematic of the model.</title><p>(A) The RBC development chain, (B) the parasite development chain, (C) the quick-acting immune response that emulates the innate response, (D) the longer-acting immune response that emulates an antibody response. <italic>T</italic>
<sub>D</sub> = average duration of a stage of development, and SD = standard deviation in the stage duration.</p></caption></fig> <fig id="pcbi-1000149-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g002</object-id><label>Figure 2</label><caption><title>Some typical simulation results.</title><p>Panels (A)–(D) show time series for several quantities for two model <italic>P. vivax</italic> infections with SD<italic><sub>I</sub></italic> = 2.4 h. <italic>T</italic>
<sub>PR</sub> is the time since primary release. For both, the source rate of RBC production is fixed, the innate response targets the last half of IBC development (late trophozoite-schizont stage), the host is antibody naïve, and the antibody response targets any part of the IBC stage. Black curve: <italic>A</italic>
<sub>th</sub> = 50 µl<sup>−1</sup>, <italic>ξ</italic> = 0.2 µl h<sup>−1</sup> for the innate response, and <italic>A</italic>
<sub>th</sub> = 0.5 µl<sup>−1</sup>, <italic>ξ</italic> = 0.008 µl h<sup>−1</sup>, 1/<italic>λ</italic>
<sub>Att</sub> = 8,760 h for the antibody response. Gray curve: <italic>A</italic>
<sub>th</sub> = 50 µl<sup>−1</sup>, <italic>ξ</italic> = 1.0 µl h<sup>−1</sup> for the innate response, and <italic>A</italic>
<sub>th</sub> = 0.5 µl<sup>−1</sup>, <italic>ξ</italic> = 0.0008 µl h<sup>−1</sup>, 1/<italic>λ</italic>
<sub>Att</sub> = 8,760 h for the antibody response. Where black and gray overlap, black is shown. (A) IBC count (inset shows the time series on a semilog plot), (B) merozoite count, (C) clearance rate of IBCs by the innate response, (D) clearance rate of IBCs by the antibody response, (E) the peak IBC count (<italic>I</italic>
<sub>PK</sub>) plotted against the integrated IBC count (<italic>I</italic>
<sub>INT</sub>) for all of the <italic>P. vivax</italic> simulations with fixed RBC production rate, (F) the RBC deficit (ΔRBC) plotted against the integrated IBC count (<italic>I</italic>
<sub>INT</sub>) for all of the <italic>P. vivax</italic> simulations with fixed RBC production rate (black) and all of the <italic>P. vivax</italic> simulations with moderate dyserythropoiesis (gray).</p></caption></fig> <fig id="pcbi-1000149-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g003</object-id><label>Figure 3</label><caption><title>Overall variations in parasitemia with changes in the model immune parameters.</title><p>The host was antibody-naïve in all simulations. (A) Plot of the peak (<italic>I</italic>
<sub>PK</sub>) and integrated (<italic>I</italic>
<sub>INT</sub>) parasitemia averaged over all simulations (results for <italic>P. vivax</italic> and <italic>P. falciparum</italic> combined) with the corresponding values of target binding affinity <italic>ξ</italic> and antigen detection threshold <italic>A</italic>
<sub>th</sub> for the innate responses. The dotted lines show the averages over simulations with no innate response. (B) Plot of <italic>I</italic>
<sub>PK</sub> and <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the corresponding values of target binding affinity <italic>ξ</italic> and threshold trigger sensitivity <italic>A</italic>
<sub>th</sub> for the antibody responses that do not target merozoites. (C) Plot of <italic>I</italic>
<sub>PK</sub> and integrated <italic>I<sub>INT</sub></italic> averaged over all simulations with the corresponding values of target binding affinity <italic>ξ</italic> for the antibody responses that target merozoites. (D) Plot of <italic>I</italic>
<sub>PK</sub> and integrated <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the corresponding values of attacker decay time <italic>T</italic>
<sub>Att</sub> for all antibody responses. In (B), (C), and (D), the dotted lines show the averages over all simulations with no antibody response. Note: In all panels, the vertical scale is logarithmic.</p></caption></fig> <fig id="pcbi-1000149-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g004</object-id><label>Figure 4</label><caption><title>Overall variation in outcome with different combinations of <italic>Plasmodium</italic> species, erythropoietic and antibody responses for model infections in antibody naïve hosts with no innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the given combination of species, erythropoietic response, and antibody target. (C) Fraction of all simulations in which the host died, and (D) fraction of simulations in which the host cleared the parasite for the given combination of species, erythropoietic response, and antibody target. (E) Color code for the data points and lines. Abbreviations for antibody responses: None: no antibodies present; IBC: any IBC stage; ET: early trophozoite; TS: late trophozoite-schizont stage; BS: bursting schizont; Mer: merozoite. Abbreviations for erythropoietic responses: SD: severe dyserythropoiesis; MD: moderate dyserythropoiesis; FX: RBC production rate fixed at basal rate; CO: compensatory erythropoiesis. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present. Note in Panel A that the vertical scale on the plots is different between the two species.</p></caption></fig> <fig id="pcbi-1000149-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g005</object-id><label>Figure 5</label><caption><title>Overall variation in outcomes with different combinations of <italic>Plasmodium</italic> species, antibody responses, and standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>) for model infections in antibody naïve hosts with no innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (C) Fraction of all simulations in which the host died for the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (D) Color code for the data points and lines. Abbreviations for antibody responses as in ##FIG##3##Figure 4##. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p></caption></fig> <fig id="pcbi-1000149-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g006</object-id><label>Figure 6</label><caption><title>Overall variation in outcome with different combinations of <italic>Plasmodium</italic> species, erythropoietic and antibody responses for model infections in antibody naïve hosts that mount an innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the given combination of species, erythropoietic response, and antibody target. (C) Fraction of all simulations in which the host died, and (D) fraction of simulations in which the host cleared the parasite for the given combination of species, erythropoietic response, and antibody target. (E) Color code for the data points and lines. Abbreviations as in ##FIG##3##Figure 4##. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present. (Note: the vertical scale for panel D is different from one in ##FIG##3##Figure 4D##.)</p></caption></fig> <fig id="pcbi-1000149-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g007</object-id><label>Figure 7</label><caption><title>Overall variation in outcome with different combinations of <italic>Plasmodium</italic> species, erythropoietic and innate responses for model infections in antibody naïve hosts that mount an innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the given combination of species, erythropoietic response, and innate immunity target. (C) Fraction of all simulations in which the host died, and (D) fraction of simulations in which the host cleared the parasite for the given combination of species, erythropoietic response, and innate target. (E) Color code for the data points and lines. Abbreviations for innate responses: IBC: any IBC stage; ET: early trophozoite; TS: late trophozoite-schizont stage. Abbreviations for erythropoietic responses as in ##FIG##3##Figure 4##. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p></caption></fig> <fig id="pcbi-1000149-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g008</object-id><label>Figure 8</label><caption><title>Overall variation in outcome with different combinations of <italic>Plasmodium</italic> species, standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>), and antibody responses for model infections in antibody naïve hosts that mount an innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) <italic>I</italic>
<sub>INT</sub> averaged over all simulations with the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (C) Fraction of all simulations in which the host died, and (D) fraction of simulations in which the host cleared the parasite for the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (E) Color code for the data points and lines. Abbreviations for antibody responses as in ##FIG##3##Figure 4##. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p></caption></fig> <fig id="pcbi-1000149-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.g009</object-id><label>Figure 9</label><caption><title>Overall variation in outcome with different combinations of <italic>Plasmodium</italic> species, standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>) and innate responses for model infections in antibody naïve hosts that mount an innate response.</title><p>(A) <italic>I</italic>
<sub>PK</sub> and (B) I<sub>INT</sub> averaged over all simulations with the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and innate immunity target. (C) Fraction of all simulations in which the host died, and (D) fraction of simulations in which the host cleared the parasite for the given combination of species, σ<sub>IBC</sub>, and innate immunity target. (E) Color code for the data points and lines. Abbreviations for innate responses as in ##FIG##6##Figure 7##. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p></caption></fig>
<table-wrap id="pcbi-1000149-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.t001</object-id><label>Table 1</label><caption><title>Values of Parameters for Model Parasites.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Species or Strain</td><td align="left" rowspan="1" colspan="1">
<italic>p</italic>
</td><td align="left" rowspan="1" colspan="1">R<sub>0</sub>
</td><td align="left" rowspan="1" colspan="1">Standard deviation <italic>σ</italic>
<sub>IBC</sub> in IBC duration (h)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>P vivax</italic>
</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">1.2, 2.4, 4.8, or 9.6</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>P falciparum</italic>
</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">1.2, 2.4, 4.8, or 9.6</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000149-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.t002</object-id><label>Table 2</label><caption><title>Values of Parameters for Model Erythropoietic Responses.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Name</td><td align="left" rowspan="1" colspan="1">1/<italic>λ</italic>
<sub>ES</sub> (h)</td><td align="left" rowspan="1" colspan="1">
<italic>δ</italic>
</td><td align="left" rowspan="1" colspan="1">ES<sub>MX</sub>/ES<sub>0</sub>
</td><td align="left" rowspan="1" colspan="1">ES<sub>MN</sub>/ES<sub>0</sub>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Fixed</td><td align="left" rowspan="1" colspan="1">Infinite</td><td align="left" rowspan="1" colspan="1">–</td><td align="left" rowspan="1" colspan="1">–</td><td align="left" rowspan="1" colspan="1">–</td></tr><tr><td align="left" rowspan="1" colspan="1">Compensatory</td><td align="left" rowspan="1" colspan="1">48</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.65</td></tr><tr><td align="left" rowspan="1" colspan="1">Moderate dyserythropoietic</td><td align="left" rowspan="1" colspan="1">48</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.75</td></tr><tr><td align="left" rowspan="1" colspan="1">Severe dyserythropoietic</td><td align="left" rowspan="1" colspan="1">48</td><td align="left" rowspan="1" colspan="1">8.5</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.65</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000149-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.t003</object-id><label>Table 3</label><caption><title>Parameter Values for Model Innate Responses.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Target</td><td align="left" rowspan="1" colspan="1">
<italic>A</italic>
<sub>th</sub> (µl<sup>−1</sup>)</td><td align="left" rowspan="1" colspan="1">
<italic>ξ</italic> (µl h<sup>−1</sup>)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">All IBCs</td><td align="left" rowspan="1" colspan="1">0.5, 5, or 50</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">Early trophozoites</td><td align="left" rowspan="1" colspan="1">0.5, 5, or 50</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">Late trophozoites and schizonts</td><td align="left" rowspan="1" colspan="1">0.5, 5, or 50</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">No target</td><td align="left" rowspan="1" colspan="1">–</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000149-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000149.t004</object-id><label>Table 4</label><caption><title>Values of Parameters for Model Antibody Responses.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Target</td><td align="left" rowspan="1" colspan="1">
<italic>A</italic>
<sub>th</sub> (µl<sup>−1</sup>)</td><td align="left" rowspan="1" colspan="1">
<italic>ξ</italic> (µl h<sup>−1</sup>)</td><td align="left" rowspan="1" colspan="1">1/<italic>λ</italic>
<sub>Att</sub> (h)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">All IBCs</td><td align="left" rowspan="1" colspan="1">0.5, 5</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td><td align="left" rowspan="1" colspan="1">160, 720, or 8760</td></tr><tr><td align="left" rowspan="1" colspan="1">Early trophozoites</td><td align="left" rowspan="1" colspan="1">0.5, 5</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td><td align="left" rowspan="1" colspan="1">160, 720, or 8760</td></tr><tr><td align="left" rowspan="1" colspan="1">Late trophozoites and schizonts</td><td align="left" rowspan="1" colspan="1">0.5, 5</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td><td align="left" rowspan="1" colspan="1">160, 720, or 8760</td></tr><tr><td align="left" rowspan="1" colspan="1">“Bursting” schizonts</td><td align="left" rowspan="1" colspan="1">0.5, 5</td><td align="left" rowspan="1" colspan="1">0.0008, 0.008, or 0.08</td><td align="left" rowspan="1" colspan="1">160, 720, or 8760</td></tr><tr><td align="left" rowspan="1" colspan="1">Merozoites</td><td align="left" rowspan="1" colspan="1">0.01</td><td align="left" rowspan="1" colspan="1">10, 100, 1000 or 10000</td><td align="left" rowspan="1" colspan="1">160, 720, or 8760</td></tr><tr><td align="left" rowspan="1" colspan="1">No target</td><td align="left" rowspan="1" colspan="1">–</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">–</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula><label>(1)</label></disp-formula> <disp-formula><label>(2)</label></disp-formula> <disp-formula><label>(3)</label></disp-formula> <disp-formula><label>(4)</label></disp-formula> <disp-formula><label>(5)</label></disp-formula> <disp-formula><label>(6)</label></disp-formula> <disp-formula><label>(7)</label></disp-formula> <disp-formula><label>(8)</label></disp-formula> <disp-formula><label>(9)</label></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(10)</label></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pcbi.1000149.s001"><label>Figure S1</label><caption><p>Plots illustrating results for certain subsets of simulated infections. (A) Peak parasitemia <italic>I</italic>
<sub>PK</sub> versus integrated parasitemia <italic>I</italic>
<sub>INT</sub> for all model <italic>P. falciparum</italic> infections with fixed RBC rate of production and antibody naïve host. (B) RBC deficit ΔRBC versus integrated parasitemia <italic>I</italic>
<sub>INT</sub> for all model <italic>P. falciparum</italic> infections with antibody naïve hosts and with fixed RBC rate of production (black points) and or with moderate dyserythropoiesis (gray points). (C) Peak parasitemia <italic>I</italic>
<sub>PK</sub> versus integrated parasitemia I<sub>INT</sub> for all model <italic>P. vivax</italic> infections with fixed RBC rate of production with pre-existing antibodies that attack all IBCs. (D) RBC deficit ΔRBC versus integrated parasitemia <italic>I</italic>
<sub>INT</sub> for all model <italic>P. vivax</italic> infections with pre-existing antibodies that attack all IBCs and with fixed RBC rate of production (black points) and or with moderate dyserythropoiesis (gray points).</p><p>(0.88 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s002"><label>Figure S2</label><caption><p>Overall variation in variation in anemia with different combinations of <italic>Plasmodium</italic> species, erythropoietic, antibody and intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>) for model infections in antibody naïve hosts with no innate response. (A) ΔRBC averaged over all simulations with the given combination of species, erythropoietic response, and antibody target. (B) ΔRBC averaged over all simulations with the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (C) Color code for the lines in panels (A) and (B). Abbreviations for antibody and erythropoietic responses as in ##FIG##3##Figure 4## of text. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p><p>(0.25 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s003"><label>Figure S3</label><caption><p>Variation in time until resolution of infection with different combinations of <italic>Plasmodium</italic> species, erythropoietic and antibody responses for infections in antibody naïve hosts without innate responses. (A) Time after primary release until death averaged over all those infections which ended in death of host by anemia, and (B) time from primary release until clearance of parasite from host averaged over all those infections in which the host cleared the parasite within one year for the given combination of species, erythropoietic response, and antibody target. (C) Color code for the data points and lines. Abbreviations as in ##FIG##3##Figure 4## in the main text. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present. Note: there was no clearance of infections in model host that either (i) lack an innate response and an antibody response or (ii) lack an innate response but had an antibody response to bursting schizonts.</p><p>(0.20 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s004"><label>Figure S4</label><caption><p>Overall variation in outcome for different erythropoietic responses to <italic>P. vivax</italic> infections in hosts with pre-existing antibodies to IBCs of any stage but with no innate immunity. (A) <italic>I</italic>
<sub>PK</sub>, (B) <italic>I</italic>
<sub>INT</sub> and (C) ΔRBC averaged over all simulations with the given response, (D) fraction of simulations with a given response in which the host died, (E) fraction of simulations with a given response in which the parasite was cleared. Abbreviations as in ##FIG##3##Figure 4## in the main text. Lines are just to guide the eye.</p><p>(0.20 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s005"><label>Figure S5</label><caption><p>Overall variation in variation in anemia with different combinations of <italic>Plasmodium</italic> species, erythropoietic, antibody and innate responses for model infections in antibody naïve hosts with an innate response. (A) ΔRBC averaged over all simulations with the given combination of species, erythropoietic response, and antibody target. (B) Color code for the lines in panel (A). Abbreviations for antibody responses as in ##FIG##3##Figure 4## of text. (C) ΔRBC averaged over all simulations with the given combination of species, erythropoietic response, and innate target. (D) Color code for the lines in panel (C). Abbreviations for innate responses as in ##FIG##6##Figure 7## of text. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p><p>(0.41 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s006"><label>Figure S6</label><caption><p>Variation in time until resolution of infection with different combinations of <italic>Plasmodium</italic> species, erythropoietic and antibody responses for model infections in antibody naïve hosts that mount an innate responses. (A) Time after primary release until death averaged over all those infections which ended in death of host by anemia, and (B) time from primary release until clearance of parasite from host averaged over all those infections in which the host cleared the parasite within one year for the given combination of species, erythropoietic response, and antibody target. (C) Color code for the data points and lines. Abbreviations as in ##FIG##3##Figure 4## in the main text. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present. Note: no host with <italic>P. vivax</italic> infection died unless there was dyserthropoiesis. No host with <italic>P. falciparum</italic> infection died unless there was severe dyserthropoiesis.</p><p>(0.23 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s007"><label>Figure S7</label><caption><p>Overall variation in outcome for different combinations of innate and erythropoietic responses to <italic>P. vivax</italic> infections in hosts with pre-existing antibodies to IBCs of any stage and that mount an innate immunity. (A) <italic>I</italic>
<sub>PK</sub>, (B) <italic>I</italic>
<sub>INT</sub> and (C) ΔRBC averaged over all simulations with the given combination of responses, (D) fraction of simulations with a given combination in which the host died, (E) fraction of simulations with a given combination in which the parasite was cleared. Abbreviations as in ##FIG##6##Figure 7## in text. Lines are just to guide the eye.</p><p>(0.27 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s008"><label>Figure S8</label><caption><p>Overall variation in variation in anemia with different combinations of <italic>Plasmodium</italic> species, standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>), antibody and innate responses for model infections in antibody naïve hosts with an innate response. (A) ΔRBC averaged over all simulations with the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (B) Color code for the lines in panel (A). Abbreviations for antibody responses as in ##FIG##3##Figure 4## of text. (C) ΔRBC averaged over all simulations with the given combination of species, σ<sub>IBC</sub>, and innate target. (D) Color code for the lines in panel (C). Abbreviations for innate responses as in ##FIG##6##Figure 7## of text. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present.</p><p>(0.30 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s009"><label>Figure S9</label><caption><p>Variation in time until resolution of infection with different combinations of <italic>Plasmodium</italic> species, standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>), and antibody responses for model infections in antibody naïve hosts that mount an innate responses. (A) Time after primary release until death averaged over all those infections which ended in death of host by anemia and (B) time from primary release until clearance of parasite from host averaged over all those infections in which the host cleared the parasite within one year for the given combination of species, <italic>σ</italic>
<sub>IBC</sub>, and antibody target. (C) Color code for the data points and lines. Abbreviations as in ##FIG##3##Figure 4## in the main text. Lines are just to guide the eye. If the data points and connecting lines for two or more antibody responses overlap in the plot, the lines are dashed to reveal all the responses present. Note: no host with either (i) antibodies to bursting schizonts, (ii) merozoites, or (iii) lacking antibodies could clear the parasite unless <italic>σ</italic>
<sub>IBC</sub>&lt;4.8 hours.</p><p>(0.27 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000149.s010"><label>Figure S10</label><caption><p>Overall variation in outcome for different combinations of standard deviations in the intraerythrocytic development time (<italic>σ</italic>
<sub>IBC</sub>) and innate responses to <italic>P. vivax</italic> infections in hosts with pre-existing antibodies to IBCs of all stages and that mount an innate immunity. (A) <italic>I</italic>
<sub>PK</sub>, (B) <italic>I</italic>
<sub>INT</sub> and (C) ΔRBC averaged over all simulations with the given combination of responses, (D) fraction of simulations with a given combination in which the host died, (E) fraction of simulations with a given combination in which the parasite was cleared. Abbreviations as in ##FIG##6##Figure 7## in text. Lines are just to guide the eye.</p><p>(0.28 MB TIF)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This research was supported by the Intramural Research Program of the National Institutes of Health through the Center for Information Technology and the John E. Fogarty International Center, US Department of Health and Human Services, Bethesda, Maryland.</p></fn></fn-group>
<graphic xlink:href="pcbi.1000149.g001"/> <graphic xlink:href="pcbi.1000149.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.g002"/> <graphic xlink:href="pcbi.1000149.g003"/> <graphic xlink:href="pcbi.1000149.g004"/> <graphic xlink:href="pcbi.1000149.g005"/> <graphic xlink:href="pcbi.1000149.g006"/> <graphic xlink:href="pcbi.1000149.g007"/> <graphic xlink:href="pcbi.1000149.g008"/> <graphic xlink:href="pcbi.1000149.g009"/> <graphic xlink:href="pcbi.1000149.e002.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e003.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000149-t001-1" xlink:href="pcbi.1000149.t001"/> <graphic xlink:href="pcbi.1000149.e004.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e005.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e006.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e007.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e008.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000149-t002-2" xlink:href="pcbi.1000149.t002"/> <graphic xlink:href="pcbi.1000149.e009.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e010.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000149-t003-3" xlink:href="pcbi.1000149.t003"/> <graphic xlink:href="pcbi.1000149.e011.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000149.e012.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000149-t004-4" xlink:href="pcbi.1000149.t004"/>
<media xlink:href="pcbi.1000149.s001.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s002.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s003.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s004.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s005.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s006.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s007.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s008.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s009.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000149.s010.tif"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["White", "Breman", "Braunwald", "Fauci", "Isslbachter", "Kasper", "Hauser", "Longo", "Jameson"], "given-names": ["NJ", "JG", "E", "AS", "JK", "DL", "SL", "DI", "JL"], "year": ["2001"], "article-title": ["Malaria and babesiosis."], "source": ["Harrison's Principles of Internal Medicine"], "publisher-loc": ["New York"], "publisher-name": ["McGraw-Hill"], "fpage": ["1203"], "lpage": ["1213"]}, {"label": ["8"], "element-citation": ["\n"], "surname": ["Zeyrek", "Kurcer", "Zeyrek", "Simsek"], "given-names": ["FY", "MA", "D", "Z"], "year": ["2000"], "article-title": ["Parasite density and serum cytokine levels in Plasmodium vivax malaria in Turkey."], "source": ["Parasite Immunol"], "volume": ["28"], "fpage": ["201"], "lpage": ["207"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Wickramasinghe", "Abdall"], "given-names": ["SN", "SH"], "year": ["2000"], "article-title": ["Blood and bone marrow changes in malaria."], "source": ["Ballieres Best Pract Res Clin Haematol"], "volume": ["13"], "fpage": ["277"], "lpage": ["299"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Craik"], "given-names": ["R"], "year": ["1920"], "article-title": ["A note on the erythrocytes in malaria."], "source": ["Lancet"], "volume": ["195"], "fpage": ["1110"]}, {"label": ["38"], "element-citation": ["\n"], "surname": ["Shushan", "Blitz", "Adams"], "given-names": ["M", "O", "CC"], "year": ["1937"], "article-title": ["The role of reticulocytes in malaria."], "source": ["J Lab Clin Med"], "volume": ["22"], "fpage": ["364"], "lpage": ["370"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Kitchen"], "given-names": ["SF"], "year": ["1938"], "article-title": ["The infection of reticulocytes by Plasmodium vivax."], "source": ["Am J Trop Med Hyg"], "volume": ["18"], "fpage": ["347"], "lpage": ["359"]}, {"label": ["55"], "element-citation": ["\n"], "surname": ["Lloyd"], "given-names": ["AL"], "year": ["2001"], "article-title": ["The dependence of viral parameter estimates on the assumed viral life cycle: limitations of studies of viral load data."], "source": ["Proc Biol Soc"], "volume": ["268"], "fpage": ["847"], "lpage": ["854"]}, {"label": ["56"], "element-citation": ["\n"], "surname": ["Lloyd"], "given-names": ["AL"], "year": ["2001"], "article-title": ["Destabilization of epidemic models with the inclusion of realistic distributions of infectious periods."], "source": ["Proc Roy London Soc Ser B"], "volume": ["268"], "fpage": ["985"], "lpage": ["993"]}, {"label": ["57"], "element-citation": ["\n"], "surname": ["Rapaport"], "given-names": ["SI"], "year": ["1987"], "source": ["Introduction to Hematology"], "publisher-loc": ["Philadelphia"], "publisher-name": ["J.B. Lippincott"]}, {"label": ["61"], "element-citation": ["\n"], "surname": ["Hjelm", "Carlsson", "Getahun", "Heyman"], "given-names": ["F", "F", "A", "B"], "year": ["2006"], "article-title": ["Antibody-mediated regulation of the immune response."], "source": ["Scandinavian J Immunol"], "volume": ["64"], "fpage": ["177"], "lpage": ["184"]}, {"label": ["65"], "element-citation": ["\n"], "surname": ["Cash", "Karp"], "given-names": ["JR", "AH"], "year": ["1990"], "article-title": ["A variable order Runge-Kutta method for initial-value problems with rapidly varying right-hand sides."], "source": ["ACM Trans Math Soft"], "volume": ["16"], "fpage": ["201"], "lpage": ["222"]}, {"label": ["66"], "element-citation": ["\n"], "surname": ["Press", "Teukolsky", "Vettering", "Flannery"], "given-names": ["WH", "SA", "WT", "BP"], "year": ["1992"], "source": ["Numerical Recipes in C"], "publisher-loc": ["Cambridge, UK"], "publisher-name": ["Cambridge University Press"]}]
PLoS Comput Biol. 2008 Aug 22; 4(8):e1000149