Dynamics of immune response and drug resistance in malaria infection
- David Gurarie^{1, 2}Email author and
- F Ellis McKenzie^{2}
https://doi.org/10.1186/1475-2875-5-86
© Gurarie and McKenzie; licensee BioMed Central Ltd. 2006
Received: 27 July 2006
Accepted: 11 October 2006
Published: 11 October 2006
Abstract
Background
Malaria parasites that concurrently infect a host compete on the basis of their intrinsic growth rates and by stimulating cross-reactive immune responses that inhibit each others' growth. If the phenotypes also show different drug sensitivities ('sensitive' vs. 'resistant' strains), drug treatment can change their joint dynamics and the long-term outcome of the infection: most obviously, persistent drug pressure can permit the more resistant, but otherwise competitively-inferior, strains to dominate.
Methods
Here a mathematical model is developed to analyse how these and more subtle effects of antimalarial drug use are modulated by immune response, repeated re-inoculation of parasites, drug pharmacokinetic parameters, dose and treatment frequency.
Results
The model quantifies possible effects of single and multiple (periodic) treatment on the outcome of parasite competition. In the absence of further inoculation, the dosage and/or treatment frequency required for complete clearance can be estimated. With persistent superinfection, time-average parasite densities can be derived in terms of the basic immune-regulating parameters, the drug efficacy and treatment regimen.
Conclusion
The functional relations in the model are applicable to a wide range of conditions and transmission environments, allowing predictions to be made on both the individual and the community levels, and, in particular, transitions from drug-sensitive to drug-resistant parasite dominance to be projected on both levels.
Keywords
Background
Since the 1920s it has been clear that interactions of immune-system and drug dynamics are critical to the elimination or persistence of drug-resistant malaria parasites in individual patients, but details and dynamic principles remain to be established [1, 2]. Later, in the 1980s, it was found that an isolate from a malaria infection is typically a mixture of parasite genotypes, some of which differ in antigenic profile, growth rate and drug response, and show complex population dynamics in mixed cultures [3, 4]. Here a mathematical model is developed to investigate the dynamics of parasite phenotypes in a malaria-infected host, with respect to critical interactions between their immune-mediated competition, relative drug sensitivities and persistent superinfection. The results characterize conditions under which a malaria infection dominated by relatively drug-sensitive parasites transitions to one dominated by relatively drug-resistant parasites.
Antimalarial drugs taken prophylactically or early in infection limit the development of immune responses protective against clinical attacks in subsequent infections [5, 6]. In endemic regions, antimalarial drug treatment is usually meant to cure a symptomatic malaria infection, thus, usually, a host with incomplete clinical immunity. Hosts are repeatedly re-infected and superinfected with mixtures of parasites. Within an infection, in the absence of drug, the dynamics of the competing parasites are determined in large part by their differences and similarities in stimulating and succumbing to the host's immune responses. The frequency with which drug clears parasites, including drug-resistant parasites, increases with host age, a surrogate for exposure and acquired immunity [7, 8]. Stronger, more diverse antibody responses are associated with greater success in antimalarial drug treatment [9, 10]. Thus, for instance, all else equal, both a higher prevalence of drug-resistant parasites and a higher frequency of drug treatment would be expected in younger age groups.
Understanding the aggregate effects of immunity can help to guide the introduction of new antimalarial drugs, explain discrepancies between in-vitro and in-vivo tests in drug-resistance surveys [11], and, in general, inform population-based intervention strategies. Indeed, these community patterns summarize very broad effects of immunity across many infections, so their interpretation and application can be greatly improved by a more detailed understanding of immune-drug-parasite interactions within individual infections.
Most studies of intra-host parasite-drug dynamics focus on the clearing effect of drug (the key medical objective) or on optimal strategies for treatment (dosage, frequency). These issues become less relevant under persistent superinfection in endemic areas, however. The focus here is more on longer-term outcomes of treatment (or prophylaxis), and their effects on competition among parasites. For instance, when joined to a proper model of parasite transmission, the crude estimates derived in relating the mean parasitaemia of phenotypes in a 'typical' host to drug dose and frequency should allow community-level analysis of the spread and prevalence of drug -resistant parasites.
The model developed here builds on earlier models of the within-host dynamics of malaria [12, 13] to address critical interactions between immune response and drug resistance. In the baseline model, the relative cross-reactivities, immune efficiencies and growth rates of competing parasites determine their dynamics and joint equilibria. Further developments examine the effects of repeated re-inoculations of parasites, then incorporate differences in drug sensitivity between the competing parasites, with the costs of resistance summarized by differences in intrinsic growth rate, to examine how dose, decay rate and frequency of drug administration affect competitive dynamics.
Methods
Model of intra-host competition mediated by immunity
The basic model of intra-host competition of multiple parasite phenotypes is based on the models of cross-reactive 'stimulating-clearing' immunity formulated in previous papers [12, 13]. Briefly, it describes interactions in which parasite phenotypes stimulate the production of non-specific and specific immune effectors, which clear (kill) them at rates that depend on the strength, specificity and quantity of effector.
Mathematically the model is a coupled system of differential equations for parasite densities x, y, and immune-effector densities: non-specific I and specific J,K for x,y-species respectively (see [12, 13]),
$\begin{array}{l}\dot{x}=\left[{a}_{X}-f\left(\frac{x}{{x}_{p}}+\frac{y}{{y}_{p}}\right)-{c}_{n}I-{c}_{X}\left(J+\eta K\right)\right]x+{S}_{X};\\ \dot{y}=\left[{a}_{Y}-f\left(\frac{x}{{x}_{p}}+\frac{y}{{y}_{p}}\right)-{c}_{n}I-{c}_{Y}\left(\xi J+K\right)\right]y+{S}_{Y};\left(1\right)\\ \dot{I}={\sigma}_{n}\left(x+y\right)-{\mu}_{n}I;\dot{J}={\sigma}_{X}x-{\mu}_{s}J;\dot{K}={\sigma}_{Y}y-{\mu}_{s}K;\end{array}$
In the follow up analysis of drug treatment x will represent a resistant strain, while y -a sensitive one. This is similar to the basic model of [12] but augmented with two additional features:
(i) a 'fever regulator' f(z) = F_{0}φ_{ q }(z), described by a sigmoid function φ_{ q }(z) = z^{ q }/(1 + z^{ q }), switched 'on' and 'off by the combined parasite density (over their pyrogenic thresholds x_{ p }, y_{ p }), with maximal pyrogenic removal rate F_{0}.
(ii) super-infection (inoculation) sources S_{ X };S_{ Y }, which can be either stationary (convenient for qualitative analysis and numeric simulation), or stochastic with a prescribed 'spacing distribution' between inoculations.
The fever plays a limited role in the dynamics of (1). At the initial stage of infection it serves to arrest the parasite growth about the pyrogenic threshold, while immune effectors build up, but then becomes negligible below that level. Therefore it is omitted in the analyses of subfebrile equilibria, but is included in the dynamic simulations.
The immune-stimulation terms (in the I,J,K – equations) are taken here as simple linear functions, so that the effectors are stimulated in proportion to parasite densities. Other parameters in system (1) are similar to [12, 13]:
(i) Parasite growth rates: a_{ X };a_{ Y }
(ii) Stimulation coefficients, for nonspecific and specific effectors: σ_{ n };σ_{ X };σ_{ Y }
(iii) Clearing coefficients, for nonspecific and specific effectors: c_{ n };c_{ X };c_{ Y }
(iv) Decay/deactivation rates, for nonspecific and specific effectors: μ_{ n };μ_{ s }
(v) Cross-reactivities, in terms of the relative clearing rates of the x-effector on y and the y-effector on x: 0 ≤ ξ, η < 1
System (1) can be rescaled to a nondimensional form which yields the essential (non-dimensional) parameters. The most important of those are cross-reactivities and the immune efficiencies (non-specific e_{ n }and specific e_{ x }, e_{ y }) introduced in [13]. The immune efficiencies combine the appropriate 'stimulation' and 'clearing' coefficients against the product of 'immune decay' and 'parasite growth' into a single nondimensional parameter.
System (1) can be considered a 'multi-prey, multi-predator' model of parasites and immune effectors, or alternatively as 'immune-mediated' competition between parasites x and y. The outcome of such competition depends on phenotypic traits with respect to the host immune response, as expressed through the essential, non-dimensional parameters above (cross-reactivities and immune efficiencies), as well as relative parasite growth rate α = a_{ Y }/a_{ X }. Note, for instance, that different parasite species may differ in replication rates (a_{ X }≠ a_{ Y }) and show low cross-reactivity, while different "strains" of a species (e.g. drug-sensitive or drug-resistant/tolerant) may be more strongly cross-reactive
Basic pharmaco-dynamics of drug treatment
The next development is a simple pharmaco-dynamic model for a single parasite phenotype. The degree to which drug treatment reduces the growth rate of a given parasite depends on both the pharmacokinetic characteristics of the drug and the sensitivity of the parasite to the drug. Key assumptions are that the phenotypic cost of drug resistance is a lower intrinsic growth rate in the absence of drug, and that the drug-induced reduction can be approximated from the basic drug parameters and the frequency of (periodic) treatment. The effect of 'heavy treatment' may be to remove the more sensitive strain completely, or, for repeated treatment, with superinfection, to tip the effective relative growth rate α in favour of x (below its critical levels α_{1,2}, explained in the following sections), releasing the more resistant strain from competitive constraints.
The pharmacodynamic model adopts the basic premises of [14]:
(i) after intake, drug concentration decays from its initial dose d_{0} at an exponential rate d(t) = d_{0}e^{-βt}(for simplicity, blood and tissue decay are not distinguished);
(ii) the drug removes parasites through a 'clearing function' Bφ_{ p }(d/d_{ S }), where B denotes the maximal clearing rate, d_{ S }the parasite 'sensitivity threshold' (often denoted d_{50}, as it clears parasite at 50% maximal rate), and φ_{ p }(z) – a sigmoid function φ = φ_{ p }(z) = z^{ p }/(1 + z^{ p }) with a suitable Hill exponent p.
A large p in φ means that the clearing rate remains relatively steady ≈ B over a wide range of drug concentrations, provided that d stays above d_{50}. For the numeric simulations below, the exponent p = 3, close to its estimated value for mefloquine [15].
The resulting dynamic treatment model [14] with constant parasite growth and a source (inoculation) S, takes the form
where D(t) = d(t)/d_{ S }is the dimensionless drug concentration, and b = B/a the relative parasite clearing rate. Two special cases are:
(i) D(t) = ${D}_{0}{e}^{-\beta \left(t-{t}_{0}\right)}$ – for a single dose administered at t_{0}, written in terms of a dimensionless initial dose D_{0} = d_{0}/d_{ S }≫ 1;
(ii) a periodic function (equation 5; presented below) for multiple applications with period T.
Equation (2) contains the essential pharmaco-kinetics of drug (the relative parasite clearing rate b and 'clearing pattern' φ), as well as its pharmaco-dynamics (initial dose and the frequency of treatment encoded in D(t)).
Drug treatment will enter system (1) or its rescaled version (3) as additional attrition terms: b_{X/Y}φ(D_{X/Y}(t)) in the x, y equations. Two phenotypes are assumed to have different drug-sensitivities, with x the more tolerant/resistant and y the more sensitive strain. Hence the different clearing coefficients, b_{ X }<b_{ Y }(to indicate higher maximal removal rate of y), and different drug-thresholds, expressed through their relative initial dosages, D_{ X }<D_{ Y }(y being cleared at lower concentrations of drug, compared to x). Thus strain y has a higher natural growth rate, α = a_{ Y }/a_{ x }> 1, and is competitively superior in the absence of drug, but, as shown below, drug treatment can offset this advantage.
Results
Analysis of stationary equilibria
The analysis at this stage is meant to outline the long-term effect of parasite-immune interactions. In particular, it attempts to predict the resulting parasite densities (long term 'outcome of competition') through the basic host parameters (growth rates, immune efficiencies, cross-reactivities), and does not yet involve drug treatment.
The rescaled system (1) with drug treatment represented by time-dependent concentration D = D(t) is given by a coupled differential system
$\begin{array}{l}\dot{x}={a}_{X}\left[1-f\left(x+\frac{y}{{y}_{p}}\right)-{b}_{X}\phi \left(D\right)-{e}_{n}I-{e}_{X}\left(J+\eta K\right)\right]x+{S}_{1};\\ \dot{y}={a}_{X}\left[\alpha -f\left(x+\frac{y}{{y}_{p}}\right)-{b}_{Y}\phi \left(D\right)-{e}_{n}I-{e}_{Y}\left(\xi J+K\right)\right]y+{S}_{2};\left(3\right)\\ \dot{I}={\mu}_{n}\left[\left(x+y\right)-I\right];\dot{J}={\mu}_{s}\left(x-J\right);\dot{K}={\mu}_{s}\left(y-K\right);\end{array}$
It depends on several dimensionless parameters: fever efficiency ${f}_{0}=\frac{{F}_{0}}{{a}_{X}}$, cross-reactivities 0 ≤ ξ, η < 1, and immune efficiencies: ${e}_{n}=\frac{{c}_{n}{\sigma}_{n}}{{a}_{X}{\mu}_{n}}{x}_{p}$, ${e}_{X}=\frac{{c}_{X}{\sigma}_{X}}{{a}_{X}{\mu}_{s}}{x}_{p},{e}_{Y}=\frac{{c}_{Y}{\sigma}_{Y}}{{a}_{X}{\mu}_{s}}{x}_{p}$. To facilitate analysis we drop the fever term, assuming 'subfebrile equilibria' x*, y* ≪ 1 that result from (relatively) high immune efficiencies e_{ n };e_{X,Y}≫ 1. The equilibrium equations for (3) without treatment or fever, are then reduced to the classical Volterra-Lotka system
$\begin{array}{c}F\left(x,y\right)=\left[1-\left({e}_{X}+{e}_{n}\right)x-\left(\eta {e}_{X}+{e}_{n}\right)y\right]x+{S}_{X}=0\\ G\left(x,y\right)=\left[\alpha -\left(\xi {e}_{Y}+{e}_{n}\right)x-\left({e}_{Y}+{e}_{n}\right)y\right]y+{S}_{Y}=0\end{array}\left(4\right)$
${\alpha}_{1}=\frac{{e}_{Y}+{e}_{n}}{\eta {e}_{X}+{e}_{n}}>{\alpha}_{2}=\frac{\xi {e}_{Y}+{e}_{n}}{{e}_{X}+{e}_{n}}\left(5\right)$
The three fitness regions include:
$\begin{array}{l}\alpha >{\alpha}_{1}\text{-}\text{'}\text{y-domination}\text{'}\text{;}\\ {\alpha}_{\text{1}}>\alpha >{\alpha}_{2}\text{-coexistence;}\left(\text{6}\right)\\ {\alpha}_{2}>\alpha \text{-'x-domination'}\end{array}$
Here 'domination' in the absence of inoculation sources, means complete removal of the competitor. The familiar phase-plots of such competition are illustrated in Figure 1.
The presence of stationary sources of infection S_{ X }; S_{ Y }> 0 (no matter how small) will shift all three types of equilibria to the upper-right quadrant into a stable 'coexistence state', defined by two asymptotic hyperbolae (solid-x and dashed-y, Figure 1). Now 'domination' does not mean complete removal, but a 'dominant density', e.g. y*/x* ≫ 1. The use of α as the 'control parameter' is relevant for qualitative analysis of drug intervention (below). Indeed, given two parasites with different drug-sensitivities, drug treatment effectively lowers their growth rates, and thus tips the outcome of competition in favour of the resistant strain.
This qualitative analysis makes two simplifying assumptions – a reduced 2D model in place of the full 5D (3), and the omission of fever. The latter has minor consequences for high immune efficiencies. The dimensional reduction (5D to 2D) maintains the equilibrium values x*, y*, but changes their stability types somewhat, from '2D stable nodes' (real negative eigenvalues of the Jacobian matrix), to '5D spiral sinks' (complex negative eigenvalues), (see [13] for details).
Dynamic relaxation
Basic parameters for dynamic model (3)
Fever: Hill exponents, clearing efficiency and pyrogenic parasitaemia levels | Immune efficiencies | Cross-reactivities | Immune decay rate/half-life |
---|---|---|---|
q = l2; f_{0} = F_{0}/a_{ X }= 1.2 x_{ p }= y_{ p }= 1 | e_{ n }= 3 e_{ X }= e_{ Y }= 10 | ξ = .4 η = .3 | μ_{ n }= .45(≈ 1.5 days) μ_{ s }= .014(≈ 50 days) |
The role of cross-reactivities in the equilibrium analysis above is to narrow (ξ, η ≈ 0) or widen (ξ, η ≫ 0) the α - range of coexistence. Intermediate values are chosen to indicate a relative genetic (antigenic) proximity of the competing strains.
Under this choice of parameters the predicted thresholds for the relative growth rate α (5)-(6) become α_{1} = 2.17; α_{2} = .54. Also introduced is a lower cut-off for parasite densities, x_{ c }= 10^{-11} (in dimensionless units), assuming a pyrogenic threshold x_{0} = 10^{4}/μl, corresponding to a density of 10^{-7}/μl (i.e. fewer than one parasite for the entire blood volume of an adult). If either of the densities x,y falls below x_{ c }it is set to zero – i.e. complete clearing.
The inoculation sources are considered either stationary, of strength EIR × "injected density s_{0}", or 'random,' with exponentially distributed waiting times and mean spacing EIR ('entomological inoculation rate,' i.e. the frequency of infectious bites by mosquitoes). The 'injected density' is estimated in terms of the parasites (primary merozoites) released from the liver, assuming a mean of 10 parasites (sporozoites) per mosquito inoculum, with each sporozoite developing into 30,000 primary merozoites, thus a total 300,000 per 4.5 liters, or 0.67/ml. This gives a dimensionless value, s_{0} = 6.7·10^{-5}. The infection source is distributed among the two phenotypes in different proportions meant to reflect parasite prevalences in the community: p_{ X }+ p_{ Y }= 1. Dynamic simulations with stochastic S use relatively low EIR ε = .3/day, due to computational constraints, but the use of equivalent stationary (mean) sources allows sampling of a broad range of EIR. Figure 1(b) demonstrates the effect of stationary EIR on x,y-equilibria. It remains marginal for a wide range of EIR, but high values O(10 – 100) can bring about a significant shift.
For coexistence, in Figure 2(b) the equilibrium values show little change, but the amplitude and phases of 'damped oscillations' are shifted during the initial stage. Overall the dynamic simulations confirm the earlier conclusion that inoculation (stationary or random) will change the outcome of competition to a coexistence pattern, though the two equilibria may differ by orders of magnitude.
Single drug treatment
The simple growth-treatment model (2) for a single parasite without immune response allows exact analytic solution, expressed through the multiplier function $X\left(t\right)=a\left[t-\frac{b}{\beta p}\mathrm{ln}\left(\frac{1+{D}_{0}^{p}}{1+{D}_{0}^{p}{e}^{-p\beta t}}\right)\right]$,
$x\left(t\right)={e}^{-X\left(t\right)}\left[{x}_{0}+{\displaystyle {\int}_{0}^{t}{e}^{X\left(t-\tau \right)}S\left(\tau \right)d\tau}\right]\left(7\right)$
This form allows one to estimate the drug efficacy (maximal burden reduction) and the time required to reach it. Assuming zero inoculation (S = 0) in (2) we get
$\begin{array}{l}{t}_{\mathrm{min}}=\frac{1}{\beta}ln\left[{\left(b-1\right)}^{1/p}{D}_{0}\right];\hfill \\ {X}_{\mathrm{min}}=\frac{a}{\beta p}\left\{\mathrm{ln}\left[\frac{{b}^{b}}{{\left(b-1\right)}^{b-1}}\right]+p\mathrm{ln}{D}_{0}-b\mathrm{ln}\left(1+{D}_{0}^{p}\right)\right\}\hfill \end{array}\left(8\right)$
This method also works in the case of zero parasite growth (a = 0), when equation (2) is solved for initial value x_{0}. The corresponding solution x(t) (7) would then stabilize at its minimal value, ${x}_{\mathrm{min}}/{x}_{0}={\left(1+{D}_{0}^{p}\right)}^{-b/\beta p}$. So the 'clearing dose' to bring parasitaemia down from x_{0} below x_{ c }is estimated by
${D}_{0}\ge {\left[{\left({x}_{0}/{x}_{c}\right)}^{\beta p/b}-1\right]}^{1/p}\left(9\right)$
While immune regulation complicates the dynamics of drug treatment beyond the simple model (2), it still affords some clues for possible treatment strategies. Indeed, sustained immune levels (I*; J* > 0) in the 'mixed' system (3) or a simpler 'single-strain' version would effectively lower the parasite growth rate from its natural value a to a' = a(1 - c_{ n }I* - c_{ X }J*). This would automatically raise the clearing drug-efficiency for such an 'immune-competent' host from its 'natural' value b to b' = B/a', hence to higher clearing in a shorter time, as illustrated in Figure 3 (lower left panel). This result supports earlier findings [9, 10] on the positive impact of acquired immunity on antimalarial therapies.
Thus the above has shown how a high single-treatment dose can lead to the dominance of an otherwise less-fit drug-resistant/tolerant phenotype, in the absence of inoculation. Clearly, the introduction of a y-inoculum following such treatment would restore 'y-domination,' as predicted in Figure 4(a). Successful malaria prophylaxis or treatment in endemically-exposed hosts typically requires repeated doses.
Repeated periodic treatment
Now the effect of multiple, periodically-spaced drug doses at time intervals T is explored. The drug concentration and its 'clearing function' become periodic as well, namely,
$D\left(t\right)={D}_{0}\frac{{e}^{-\beta \text{Mod}\left(t,T\right)}}{1-{e}^{-\beta T}}\left(10\right)$
where Mod(t,T) designates a periodic linear function t on interval [0,T]. Turning to the clearing function φ(D(t)), observe that a high Hill exponent of φ(z) allows φ(D(t)) to be approximated by a periodic step-function, taking value 1 on interval [0, T_{ M }] and 0 on the complimentary range [T_{ M },T], with T_{ M }≈ ln D_{0}/β. Hence as the first-order approximation the periodic φ(D(t)) can be replaced by its mean value,
$\overline{\phi}=\frac{1}{T}{\displaystyle \underset{0}{\overset{T}{\int}}\phi}\left({D}_{0}{e}^{-\beta t}\right)dt\approx \frac{\mathrm{ln}{D}_{0}}{\beta T}\left(11\right)$
Therefore approximate formula (11) is applied to get crude estimates of the treatment frequency for the establishment of the resistant strain. Assuming for simplicity b_{ X }= 0 (fully-resistant strain), the equation is obtained for two 'critical' treatment frequencies/periods, T_{ C }(for coexistence) and T_{ R }(for 'x-domination'),
${T}_{C}=\frac{b\mathrm{ln}{D}_{0}}{\beta \left(\alpha -{\alpha}_{1}\right)}>{T}_{R}=\frac{b\mathrm{ln}{D}_{0}}{\beta \left(\alpha -{\alpha}_{2}\right)};\left(12\right)$
in terms of pharmacokinetic parameters β and $b=\frac{B}{{a}_{X}}$, initial dose D_{0}, and relative growth rate $\alpha =\frac{{a}_{Y}}{{a}_{X}}$. As in section 1 the sensitive strain y is assumed competitively superior in the absence of drug, i.e. α > α_{1}, and we expect the following long-term outcomes:
(i) treatment frequency T > T_{ C }will maintain y-domination;
(ii) T_{ C }> T > T_{ R }will bring about coexistence, with alternating x,y phases;
(iii) T <T_{ R }will lead to dominance by the resistant x-strain.
To test these predictions the y-dominant case of Figure 2(a) is subjected to random inoculation with EIR = 0.3/day, distributed among 2 strains in proportions p_{ X }= .6; p_{ Y }= .4 (60% of inoculations are x-type, 40% are y). Periodic treatment (10) is then applied to such a host, and a range of values of period T sampled, using the pharmacokinetic parameters of the two strains as in the previous section, b_{ X }= 1.5;b_{ Y }= 3; d_{ x }= 2d_{ Y }.
Figure 5 shows the resulting outcomes, ranging from y-dominance at period T = 120 days, to alternating x or y dominance (T = 45 days; 'periodic coexistence'), to x-dominance at T = 20. The inoculations, though barely visible in plots, still play a role here. Indeed, without them the process would terminate (at T = 120) after the first two cycles. Namely, the first cycle would drop 'treated y' to near-extinction, after which it would rebound, at a lower level, and drive x to extinction (Figure 5, left), but a second treatment would y altogether. These numeric results confirm the earlier qualitative predictions on the 'dominance – coexistence' transition, but the 'critical' periods T_{ C };T_{ R }differ from the simple estimates (12).
The above shows that persistent treatment creates a new 'effective environment' that could significantly alter the time-averaged distribution of phenotype densities within a host. Provided that such 'over-treated' hosts make up a sizable fraction of the population and that transmission to mosquitoes is in some way proportional to parasite densities, this could further promote the spread of drug resistance through the community.
Stationary treatment
The above has demonstrated some critical effects of periodic drug treatment and superinfection on immune-mediated competition. Many applications require a more 'exact' account of the relationships among all these factors, however. Specifically, this section examines the effects of (i) drug efficiencies b_{ X }≪ b_{ Y }, (ii) treatment intensity, as measured by the (period average) factor $\overline{\phi}$ (11), and (iii) superinfection sources {S_{ X }, S_{ Y }}, through analysis of an (approximate) stationary system obtained by averaging the sources and the treatment regimen. The resulting stationary model allows all three factors to be incorporated in a simple and efficient way.
The 'averaging' proceeds by replacing the stochastic inoculation source by its steady (mean) value (S_{ X }, S_{ Y }) = (p_{ X }, p_{ Y }) $\overline{S}$, where $\overline{S}$ is the product of EIR and 'mean inocula' (merozoite release s_{0}), and (p_{ X }, p_{ Y }) – the relative fractions of the two phenotypes. By the same pattern we replace the 'periodic drug function' φ [D(t)] with its mean value (treatment intensity) (11), $\overline{\phi}$ = $\overline{\phi}$(T,D_{0}); 0 ≤ $\overline{\phi}$ ≤ 1, considered as a function of the initial dose D_{ Y }= D_{0} and treatment period T. The resulting 'average' equilibrium system is similar to (4),
$\begin{array}{l}x\left[\left(\overline{\phi}-\frac{1}{{b}_{X}}\right)+\frac{1}{{b}_{X}}\left(\frac{x}{{m}_{1}}+\frac{y}{{n}_{1}}\right)\right]={S}_{1}/{b}_{X};\\ y\left[\left(\overline{\phi}-\frac{1}{{b}_{Y}}\right)+\frac{1}{{b}_{Y}}\left(\frac{x}{{m}_{2}}+\frac{y}{{n}_{2}}\right)\right]={S}_{2}/{b}_{Y}\left(13\right)\end{array}$
But somewhat different rescaling and notations are used for the 'x,y'-intercepts (designated here by {m_{ i }, n_{ i }}), and for the stationary (mean) sources (S_{1},S_{2}) = (αp_{ X }, p_{ Y })S of relative strength S = EIR × s_{0}/a_{ Y }, where α = a_{ Y }/a_{ x }– the above 'fitness' parameter. All three factors (intensity $\overline{\phi}$, clearing rates b_{ x }≪ b_{ Y }, and sources {S_{ i }}) enter (13) explicitly. The resulting 'stationary equilibria' of (13) are expected to approximate the stable periodic (or quasi-periodic/stochastic) equilibria of the original system.
The exact solution of the 4-th order algebraic system (13) is given by (grossly cumbersome) functions: x = x* ($\overline{\phi}$, S,b_{ X },b_{ Y }) and y = y* (φ,S,b_{ X },b_{ Y }). (Note that these can be brought into a more manageable 'analytic form' under simplifying assumptions, e.g. a 'highly sensitive' y-strain.) These can easily be manipulated, however (for numeric and graphic purposes), by any symbolic algebra package (here Mathematica 5). Having four independent parameters in x*, y* some of them can be fixed (e.g. clearing rates b_{ Y };b_{ X }– for the drug-sensitive/resistant phenotypes), to examine the effect of the treatment intensity 0 <$\overline{\phi}$ < 1, and/or the source strength S.
This suggests that treatment intensity is more crucial for the onset of resistance than inoculation frequency. The increased source by itself (at all treatment levels) would raise equilibria but diminish their relative difference: x*/y* ≈ 1 at high S.
Discussion
The above has examined how immune-mediated competition between parasites is perturbed by a drug to which the competing parasites are differentially resistant, and how drug dose, drug timing, and inoculations of new parasites affect these interactions. For clarity, the examples focused on the situation in which, in the absence of drug, the drug-sensitive strain is numerically dominant in immune-mediated competition, but the model spans the full range of possibilities. That is, while the outcome of immune-mediated competition is either that one phenotype dominates, or the phenotypes coexist, further possibilities arise when phenotypes of differing drug sensitivity are subjected to treatment and superinfection. Figure 7 illustrates possible outcomes for a range of sensitivities, drug efficacies, treatment intensities, and superinfection frequencies, and shows the potential for comparing model predictions across the range of possibilities that may arise in empirical studies.
Actual frequencies of parasite inoculation vary with the size of the vector mosquito population, and, like the composition of inocula, with the prevalence and characteristics of parasites and immunity in the human population. The baseline results reflect an inoculum that initiates an infection, i.e. after any previous infections have cleared. The parasites in superinfecting inocula are in general greatly outnumbered by those in an ongoing infection, particularly when the ongoing infection is at or near its peak. In the model, parasite inoculation enters either as a random (stochastic) source, or as its stationary (mean) value. The latter allowed exploration of a broad range of source strengths (EIR): overall, increased strength drives both equilibrium densities (or their 'period means') up, but increased intensity and efficacy of treatment eventually tips the outcome from the strain that is more 'fit' in terms of a host immune response to the one that is more 'fit' in terms of a drug.
What is transmitted from an infected human to a mosquito, and onward, poses another complex set of questions, but our results here seem to support the hypothesis that immune-mediated interactions can shape the spread of drug resistance, even if the phenotypic traits are not linked genetically [19, 20]. The fitness cost of drug resistance, considered here simply in terms of replication rate, is likely to be multifactorial in a population of hosts heterogeneous with respect to infection histories and immune profiles [21]. It would be interesting, in future work, to explicitly consider the genetics of parasite drug resistance, with respect to both origin and spread, and the common use of anti-malarial drugs to "treat" non-malarial fevers or other symptoms.
Antigenic variation can also play a role in parasite competition. From the standpoint of this model the only potentially relevant effect is a change of immune stimulation/clearing, since parasites' intrinsic replication or drug sensitivities should not be affected. A simple way to accommodate antigenic variation (in lieu of more complicated 'multi-strain/multi-clone' models [22, 23]) is to make the specific immune clearing function to change in time c(t). A drop in c can be thought of as resulting from a new variant of 'low cross-reactivity' to prior antibodies, taking over and growing into a dominant strain. As the clone keeps replicating, specific antibodies develop, and function c(t) gradually 'relaxes' to its normal (relatively high) value. Such a 'random drop + relaxation' form of c(t) was proposed in [13].
Time dependent c(t) would make the model non-stationary, the same way as the 'variable inoculation source' or 'non-stationary treatment'. In each case the strategy here was to 'average variable coefficients' over time. Applying the same 'averaging methodology' to variable clearing rates would replace them with somewhat lower 'mean values' {$\overline{c}$_{ X }; $\overline{c}$_{ Y }}, and the new 'effective' c would enter formulae and analyses. Indeed, lower {e_{ X };e_{ Y }} (proportional to {$\overline{c}$_{ X };$\overline{c}$_{ Y }}) would change (decrease) 'fitness thresholds' {α_{1},α_{2}}, and drive up equilibrium levels of x,y (as both are cleared less 'efficiently' on average).
Thus the implication of antigenic variation is that degrees of cross-reactivity between parasites may fluctuate during the course of an infection, and its net effect would be to (effectively) lower the immune efficiencies, and change the related equilibria, the fitness ranges of 'domination and coexistence' etc. The major points of interest here have to do with the relationship between the parasite entities, however, which would shift if antigenic variation differs between them in some systematic, biased way, which would then be expressed in terms of relative competitive advantage. Thus most of our results would maintain their qualitative form, but some quantitative changes on the predicted outcomes of treatment would be expected (Figures 6, 7). It might prove interesting, in future work, to explicitly consider the extent to which antigenic variation gives rise to fevers which give rise to drug-taking, for instance.
Apparently no previous work has examined these critical within-host interactions between parasite, immune-response and drug dynamics in a malaria infection, but several models have touched on important aspects of these analyses or on closely-related issues. Davis and Martin [24] compared several simple descriptive models of parasite clearance dynamics during curative drug treatment. Based on the results of pharmacokinetic-pharmacodynamic models, Hoshen et al. [25] argued that well-timed follow-up doses might eliminate even resistant phenotypes, and Simpson et al. [15] that, to prevent resistance, larger doses should be deployed as standard at the first introduction. While these models did not consider host immune responses or parasite replication, Austin et al. [14] noted that "most drugs act best against replicating pathogens in combination with effective immunological responses," joined a pharmacokinetic model to a simple model of pathogen-immune dynamics [26], and derived drug doses and frequencies necessary to reduce the average lifespan of infected RBCs below a critical threshold. Their model did not consider drug- or immune-mediated competition between strains, however.
Gatton et al. [27] modeled the risk of a drug-resistant mutant arising during an infection, in terms of the rate at which a single parasite genotype switches between antigenic variants, the response rate of the corresponding specific antibody, and a simple time-of-treatment model. They found mutants most likely to arise in hosts lacking any such specific antibodies, in hosts treated before antibody response to a switch of variants, and, in recent work [28], in hosts taking drugs with long half-lives or in subcurative doses. Hastings [29] recently developed parasite-population-genetic models that represent competing forces behind shifts in equilibrium gene frequencies in terms of within-host averages and proportions; though these models do not yet encompass considerations of immune responses or other dynamic factors, he emphasized that "intense competition between separate malaria clones co-infecting the same human can generate complex dynamics," and that "the dynamics underlying the evolution of antimalarial resistance may therefore be much more complex than previously realized." This appears irrefutable, hence the model above was developed to address these "complex dynamics."
Declarations
Acknowledgements
The research was supported by the Fogarty International Center, U.S. National Institutes of Health. The authors acknowledge critical input from anonymous referees and stimulating discussions with W.P. O'Meara, B.C. Sorkin, C.H. King, A.L. Graham.
Authors’ Affiliations
References
- Yorke W, Macfie JWS: Observations on malaria made during treatment of general paralysis. Trans Roy Soc Trop Med Hyg. 1924, 18: 13-33. 10.1016/S0035-9203(24)90664-X.View ArticleGoogle Scholar
- Wernsdorfer W, (ed): Drug-resistant malaria. UNDP/World Bank/WHO Geneva. 1982Google Scholar
- Thaithong S, Beale GH, Fenton B, McBride J, Rosario V, Walker A, Walliker D: Clonal diversity in a single isolate of the malaria parasite Plasmodium falciparum. Trans Roy Soc Trop Med Hyg. 1984, 78: 242-245. 10.1016/0035-9203(84)90287-6.View ArticlePubMedGoogle Scholar
- Willet GP, Milhous WK, Gerena L, Oduola AMJ: Mixed population dynamics in human malaria parasite cultures. Trans Roy Soc Trap Med Hyg. 1991, 85: 33-34. 10.1016/0035-9203(91)90142-L.View ArticleGoogle Scholar
- James SP: Some general results of a study of induced malaria in England. Trans Roy Soc Trop Med Hyg. 1931, 27: 477-538. 10.1016/S0035-9203(31)90068-0.View ArticleGoogle Scholar
- Greenwood BM: Malaria chemopropylaxis in endemic regions. Malaria, waiting for the vaccine. Edited by: Targett GAT. 1991, John Wiley, Chichester UK, 83-102.Google Scholar
- Dorsey G, Kamya MR, Ndeezi G, Babirye JN, Phares CR, Olsen JE, Katabira ET, Rosenthal PJ: Predictors of chloroquine treatment failure in children and adults with falciparum malaria in Kampala, Uganda. Am J Trap Med Hyg. 2000, 62: 686-692.Google Scholar
- Djimde AA, Doumbo OK, Traore O, Guindo AB, Kayentao K, Diourte Y, Niare-Doumbo S, Coulibaly D, Kone AK, Cissoko Y, Tekete M, Fofana B, Dicko A, Diallo DA, Wellems TE, Kwiatkowski D, Plowe CV: Clearance of drug-resistant parasites as a model for protective immunity in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2003, 69: 558-563.PubMedGoogle Scholar
- Mayxay M, Chotivanich K, Pukrittayakamee S, Newton P, Looareesuwan S, White NJ: Contribution of humoral immunity to the therapeutic response in falciparum malaria. Am J Trop Med Hyg. 2001, 65: 918-923.PubMedGoogle Scholar
- Mawili-Mboumba DP, Borrmann S, Cavanagh DR, McBride JS, Matsiegui PB, Missinou MA, Kremsner PG, Ntoumi F: Antibody responses to Plasmodium falciparum merozoite surface protein-1 and efficacy of amodiaquine in Gabonese children with P. falciparum malaria. J Inf Dis. 2003, 187: 1137-1141. 10.1086/368414.View ArticleGoogle Scholar
- Ringwald P, Basco LK: Comparison of in vivo and in vitro tests of resistance in patients treated with chloroquine in Yaounde, Cameroon. Bull World Health Organ. 1999, 77: 34-43.PubMed CentralPubMedGoogle Scholar
- Mason DP, McKenzie FE: Blood-stage dynamics and clinical implications of mixed Plasmodium vivax-Plasmodium falciparum infections. Am J Trop Med Hyg. 1999, 61: 367-374.PubMed CentralPubMedGoogle Scholar
- Gurarie D, Zimmerman PA, King CH: Dynamic regulation of single- and mixed-species malaria infection: Insights to specific and non-specific mechanisms of control. J Theor Biol. 2006, 240: 185-199. 10.1016/j.jtbi.2005.09.015.View ArticlePubMedGoogle Scholar
- Austin DJ, White NJ, Anderson RM: The dynamics of drug action on the within-host population growth of infectious agents: melding pharmacokinetics with pathogen population dynamics. J Theor Biol. 1998, 194: 313-339. 10.1006/jtbi.1997.0438.View ArticlePubMedGoogle Scholar
- Simpson JA, Watkins ER, Price RN, Aarons L, Kyle DE, White NJ: Mefloquine pharmacokinetic-pharmacodynamic models: implications for dosing and resistance. Antimicrob Agents Chemother. 2000, 44: 3414-3424. 10.1128/AAC.44.12.3414-3424.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Rouzine IM, McKenzie FE: Link between immune response and parasite synchronization in malaria. Proc Nat Acad Sci USA. 2003, 100: 3473-3478. 10.1073/pnas.262796299.PubMed CentralView ArticlePubMedGoogle Scholar
- McQueen PG, McKenzie FE: Age-structured red blood cell susceptibility and the dynamics of malaria infections. Proc Nat Acad Sci USA. 2004, 101: 9161-9166. 10.1073/pnas.0308256101.PubMed CentralView ArticlePubMedGoogle Scholar
- Kwiatkowski D: Febrile temperatures can synchronize the growth of Plasmodium falciparum in vitro. J Exp Med. 1989, 169: 357-361. 10.1084/jem.169.1.357.View ArticlePubMedGoogle Scholar
- McKenzie FE, Wong RC, Bossert WH: Discrete-event models of mixed- phenotype Plasmodium falciparum malaria. Simulation. 1999, 73: 213-217.PubMed CentralView ArticlePubMedGoogle Scholar
- McKenzie FE, Ferreira MU, Baird JK, Snounou G, Bossert WH: Meiotic recombination, cross-reactivity and persistence in Plasmodium falciparum. Evolution. 2001, 55: 1299-1307. 10.1554/0014-3820(2001)055[1299:MRCRAP]2.0.CO;2.PubMed CentralView ArticlePubMedGoogle Scholar
- Babiker HA, Satti G, Ferguson H, Bayoumi R, Walliker D: Drug resistant Plasmodium falciparum in an area of seasonal transmission. Acta Trop. 2005, 94: 260-268.View ArticlePubMedGoogle Scholar
- Molineaux L, Diebner HH, Eichner M, Collins WE, Jeffery GM, Dietz K: Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitology. 2001, 122: 379-391. 10.1017/S0031182001007533.View ArticlePubMedGoogle Scholar
- Paget-McNicol S, Gatton M, Hastings I, Saul A: The Plasmodium falciparum var gene switching rate, switching mechanism and patterns of parasite recrudescence described by mathematical modelling. Parasitology. 2002, 124: 225-235. 10.1017/S0031182001001160.View ArticlePubMedGoogle Scholar
- Davis TM, Martin RB: Clearance of young parasite forms following treatment of falciparum malaria in humans: comparison of three simple mathematical models. Epidemiol Infect. 1997, 119: 61-69. 10.1017/S0950268897007693.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoshen MB, Stein WD, Ginsburg H: Modelling the chloroquine chemotherapy of falciparum malaria: the value of spacing a split dose. Parasitology. 1998, 116: 407-416. 10.1017/S0031182098002480.View ArticlePubMedGoogle Scholar
- Hetzel C, Anderson RM: The within-host cellular dynamics of bloodstage malaria: theoretical and experimental studies. Parasitology. 1996, 113: 25-38.View ArticlePubMedGoogle Scholar
- Gatton ML, Hogarth W, Saul A: Time of treatment influences the appearance of drug-resistant parasites in Plasmodium falciparum infections. Parasitology. 2001, 123: 537-546. 10.1017/S0031182001008824.View ArticlePubMedGoogle Scholar
- Gatton ML, Martin LB, Cheng Q: Evolution of resistance to sulfadoxine- pyrimethamine in Plasmodium falciparum. Antimicrob Agents Chemother. 2004, 48: 2116-2123. 10.1128/AAC.48.6.2116-2123.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Hastings IM: Complex dynamics and stability of resistance to antimalarial drugs. Parasitology. 2006,Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.