Prospective strategies to delay the evolution of anti-malarial drug resistance: weighing the uncertainty
© Smith et al; licensee BioMed Central Ltd. 2010
Received: 19 February 2010
Accepted: 23 July 2010
Published: 23 July 2010
The evolution of drug resistance in malaria parasites highlights a need to identify and evaluate strategies that could extend the useful therapeutic life of anti-malarial drugs. Such strategies are deployed to best effect before resistance has emerged, under conditions of great uncertainty.
Here, the emergence and spread of resistance was modelled using a hybrid framework to evaluate prospective strategies, estimate the time to drug failure, and weigh uncertainty. The waiting time to appearance was estimated as the product of low mutation rates, drug pressure, and parasite population sizes during treatment. Stochastic persistence and the waiting time to establishment were simulated as an evolving branching process. The subsequent spread of resistance was simulated in simple epidemiological models.
Using this framework, the waiting time to the failure of artemisinin combination therapy (ACT) for malaria was estimated, and a policy of multiple first-line therapies (MFTs) was evaluated. The models quantify the effects of reducing drug pressure in delaying appearance, reducing the chances of establishment, and slowing spread. By using two first-line therapies in a population, it is possible to reduce drug pressure while still treating the full complement of cases.
At a global scale, because of uncertainty about the time to the emergence of ACT resistance, there was a strong case for MFTs to guard against early failure. Our study recommends developing operationally feasible strategies for implementing MFTs, such as distributing different ACTs at the clinic and for home-based care, or formulating different ACTs for children and adults.
Plasmodium falciparum, which causes malaria, is the most important parasite species that infects humans with approximately 2.37 billion people at risk [1, 2]. Prompt effective drug treatment can reduce the risk of mortality for those with clinical infections, and is a key component of malaria elimination and eradication plans both past and present . Diminished therapeutic efficacy due to the evolution of resistance to previous first-line drugs, however, contributed to the failure of initial eradication efforts and resulted in increases in infection and mortality . Drug stewardship, combination therapies, and other policies have been advocated to slow the evolution of resistance to anti-malarial drugs, however to be effective, a policy must delay emergence or at least slow the spread of resistant parasites. Combination of drugs, such as artemisinin combination therapy (ACT), delay emergence by eliminating resistant mutants except those that carry two different mutations. Since each drug can eliminate mutants that are resistant to the other component, mutations to both components must arise in the same parasite. Combination therapies thus delay the time until a viable resistant parasite appears . Alternative strategies call for reducing selection for the spread of resistance by using drugs more prudently or by rationing drugs . Such strategies assume that the rate of spread increases in proportion to the rate of drug use in a population. Multiple first-line therapies (MFTs) on the other hand, reduce pressure on each drug, defined as the proportion of clinical episodes that are treated, and slow spread while avoiding the ethical problem of leaving some patients untreated [7, 8].
De novo resistant mutants may fail to become established, as the parasite is present in only a few hosts, and thus its persistence is a matter of chance [10, 11] (Figure 1, middle). The fitness advantage for drug-resistant parasites under some level of drug pressure is partially offset by a fitness disadvantage, or a biological cost of resistance, if drugs are absent [12–15]. The biological cost will tend to be highest when resistance first appears and will reduce the probability of establishment because mutations that lower fitness tend to die out [12, 13]. Therefore, reducing drug pressure would also reduce the chances of establishment by increasing the effect of a high initial cost of resistance.
Thus, reducing drug pressure impacts the evolution of resistance in three ways: it delays its appearance, reduces the likelihood of its establishment, and slows its spread (Figure 1, bottom). However, the evolution of drug resistance involves complex nonlinear processes, and thus substantial uncertainty exists about the relative impact of deploying MFTs because of unknown mutation rates, poor information about drug pressure, and a poorly defined relationship between drug pressure and fitness. Weighing uncertainty requires a quantitative approach using mathematical models. In this paper, the effectiveness of deploying multiple ACTs to manage resistance was evaluated using a hybrid modelling approach that weighs all three potential benefits of deploying MFTs.
Most new first-line anti-malarial drugs as well as most anti-malarial drugs in the development pipeline contain an artemisinin-class drug. To extend their useful therapeutic lives, artemisinins are paired with a partner drug and given as combination therapies. While several different ACTs are already available on the market, artemisinin monotherapy is still readily available in many places. Cross-resistance among all the artemisinin derivatives is expected, and the evolution of resistance to a partner drug should be comparatively easy for a parasite with artemisinin resistance, so the loss of one ACT would probably lead to the rapid loss of all ACTs. Because new non-artemisinin anti-malarial drugs will likely not be available for a decade or more, preserving artemisinin is crucial. A significant threat to artemisinin comes from the evolution of resistance to the partner drugs, whose pharmacokinetic half-lives are generally significantly longer than the artemisin-class drugs they are paired with. This has led to concerns about selection for partner drug resistance after the artemisinin concentrations have waned [5, 16]. MFTs would provide one way for the partner drugs to protect each other and, by extension, artemisinin. To examine the full benefit of MFTs, a quantitative approach to assessing the waiting time to ACT failure was developed.
Appearance, defined as an infection with de novo mutants that are frequent enough within a single host to be transmitted to another host, requires a resistant mutant and treatment (Figure 1, top). The probability that at least one resistant mutant is present when chemotherapy begins is determined by the probability of spontaneously generating a resistant mutant per cell division, m, and the parasite population size at the time of treatment, P. The probability of appearance is thus approximately mP (assuming it is small, i.e. mP < 0.01).
After appearance, a parasite can only become established by spreading to other hosts. Parasite fitness is measured in terms of the replacement number, the expected number of human hosts that would be infected for each human host infected. The fitness of a resistant phenotype is determined by its biological cost and by the advantage conferred from drug pressure. If a parasite cannot replace itself, it will never become established, though parasite fitness may be too low for it to become established initially. The fitness of a resistant phenotype will tend to increase if they persist, as they can acquire compensatory mutations to overcome the biological cost of resistance. Both persistence and the acquisition of compensatory mutations are chance events, though the compensatory mutations increase the likelihood of persistence.
Establishment of resistance was considered in a population where malaria is endemic, and where on average parasites replace themselves approximately once. Parasites with de novo resistance are rare in a population when they first appear, by definition, and when replacement numbers are close to one, they frequently fail to establish.
The probability of establishment, Q, and the number of parasite generations required to establish, G, were estimated by simulating an evolving branching process [10, 11, 18]. A branching process does not keep track of when the progeny were produced, but rather the number of offspring in each generation that has elapsed (Figure 1, middle). The number of "offspring" (i.e., new human infections) was drawn from a negative binomial distribution to account for heterogeneous biting [19, 20]. The offspring evolved by inheriting a slightly different fitness than their parents. The branching process was iterated until no parasite offspring remained (a failure) or there were at least 100 offspring with replacement numbers greater than one (a success) (Additional file 1). It was repeated n times until there were s successes, and Q ≈ (s-1)/(n + s-1). Assuming that the time required to complete a parasite generation was T g (i.e. the average time from one human infection to another, at least one month and generally approximately two months), the waiting time to emerge was
Failure to establish has the effect of delaying emergence as much as requiring an additional mutation that occurs with frequency Q.
The spread of resistance was simulated using a simple epidemiological model that describes the dynamics of malaria. The models consider the fraction of humans in a population who are asymptomatically infected with drug-resistant or drug-sensitive parasites, x or w, respectively (Figure 1, bottom). The models also consider the incidence of clinical malaria in relation to infection status in a very general way (Additional file 1).
A therapy was considered to have "failed" when the frequency of resistance (i.e. x/(x + w)) exceeded a predefined threshold; the WHO threshold of 10% was adopted. The end of the branching process was the starting point for the epidemiological model, so in a population of H humans, the initial frequency of resistance was 100/H humans. The time to reach a frequency of 10% was longer in larger human populations where the initial frequency was lower (Figure 1, bottom).
Results and Discussion
Reducing parasite fitness
Reducing drug pressure and time to failure
Reducing drug pressure also slows the spread of resistance. The benefits of extending a policy are complicated because of a threshold on drug pressure, a tipping point that favors drug-sensitive or drug-resistant parasites (Additional file 1). If reducing drug pressure takes selection from one side of the tipping point to the other, then it would delay drug failure indefinitely (Figure 5b). The quantitative benefits of reducing drug pressure, however, are much lower if selection stays on the same side of the tipping point. If drug pressure were low enough that drug-sensitive parasites were favored, there would be no need to reduce drug pressure further. If drug pressure were sufficiently high that it favored drug-resistant parasites before and after adopting a policy, reducing drug pressure would approximately double the time to failure (Figure 5b).
Decisions to preserve the efficacy of anti-malarial drugs are often made over a finite time horizon; here, the goal was preserving drug efficacy over 20 years. The introduction of a finite time horizon for planning is similar to the analysis of tipping points. If the goal of a policy is to preserve a drug for at least 20 years, then the analysis must reduce drug pressure low enough to achieve this goal. If the waiting time to failure for a single ACT were longer than 20 years, then there would be no need to do anything. A maximum benefit of deploying MFTs is found when a single ACT would delay resistance by about a decade. The plot of years gained up to a maximum of 20 years plotted against drug pressure, has the shape of a shark fin (Figure 5b). The combined effects of delayed emergence and delayed spread all combine to extend the useful therapeutic life of a drug (Figure 5c).
The most important principle is that the benefits of implementing a policy in terms of a waiting time to failure are all proportional to the waiting time to failure in the absence of any policy. It follows that a policy to preserve a drug is worth implementing only if the waiting time to fail is in the right range; if the waiting time to fail is less than one year, then a policy must extend the useful life of a drug by a factor of 20 or more, a practical impossibility. If the drug already has a waiting time of more than 20 years, then a policy is not really necessary. This principle transcends the uncertainty and makes it possible to say something more concrete about which drugs are worth protecting with a policy.
Uncertainty and scale
The time to failure depends on the drug pressure, the mutation rates, the biological cost of resistance over time, the clinical incidence, and the generation time of the parasites. There is substantial uncertainty about all these parameters. To weigh the value of deploying MFTs, a probability distribution function was assigned to the parameters to compute times to failure (Additional file 1). The answers were generally uninformative.
Because there was so much uncertainty in all of the parameters, a slightly different question was asked: when would MFTs be most likely to provide a benefit? Asking the question in this way and doing a sensitivity analysis showed that the scale of the question was a dominant factor in the analysis; the more people with parasites who were treated, the shorter the time to appearance, but the longer the waiting time to reach a frequency of 10% because of a smaller initial frequency of resistance. The question of the appropriate scale came to dominate the analysis when it was noted that reducing the scale could increase the waiting time to appearance by millennia, depending on the size of the population at risk. Given the potential for ACT resistance to emerge anywhere and spread to nearly everywhere, as was the case with resistance to chloroquine and sulfadoxine-pyrimethamine , the global scale was considered to be appropriate.
Parasite population density at the time of treatment can vary by several orders of magnitude, which has a strong effect on the probability of appearance. Clinical malaria in an immunologically naïve adult can produce parasite densities that range up to 10 12 parasites, in extreme cases. In areas with highly endemic malaria and clinical immunity, adults tend to have less clinical malaria with lower parasite densities at the time of treatment. By comparison, children are often immunologically naïve, but they have less blood volume and a smaller number of parasites. Reasonably large parasite population sizes were considered (assuming 10 9 < P < 10 12 ).
The probability of appearance in an individual is extremely low, thus the waiting time to emergence somewhere in the world must account for the large populations at risk and global demand for drugs. There are approximately 450 million clinical malaria episodes per year, with a credible range of 350-550 million cases . If ACT were scaled up to meet global demand, there could be as many as 5 billion courses of ACT per decade. For a single combination therapy to delay the appearance of resistance by a decade, assuming that 1 billion people were treated over a decade, the probability a resistant mutant appears per treatment would have to be less than approximately 10 -18 .
Mutation rates in the dihydrofolate reductase gene of P. falciparum are estimated to be approximately 2.5 × 10-9 . If ACT resistance required two such mutations of similar frequency, the mutation rate required for resistance to an ACT would be closer to 10 -17 . By this arithmetic, the waiting time to appearance somewhere in the world is close to a year. For a reference set of parameters (Additional file 1), the estimate of the waiting time to failure, after scaling up ACT, was approximately eight years. There was, it should be recalled, very little confidence in this particular estimate.
One notable feature of this analysis was that it combined very small numbers (mutation rates) with very large numbers (parasite population densities, the global demand for anti-malarial drugs). A small change in one of the exponents could change the time to failure by a factor of 10 or more. The proportional changes that would be achieved by implementing MFTs were most important for policy when the waiting times were large enough to have a noticeable effect but not already sufficient to preserve ACTs for 20 years. Curiously, when the big numbers were combined, they roughly balanced, and the waiting time to appearance, at a point in time when ACTs are the de facto global first-line treatment for malaria, would be close to one year, precisely the value where MFTs would have maximum benefit. As an example, MFTs delayed the waiting time to failure from 8 to 17 years (Figure 2).
Treating fewer patients would likely reduce selection pressure and slow the spread of resistance, but there are few practical ways of reducing drug pressure without putting lives at risk. MFTs are one way: treating half the population with one ACT and the other half with a different ACT would reduce the fitness of resistant parasites approximately as well as treating half as many patients . A simple rule of thumb for deciding whether to choose MFTs is that if there is an advantage to cutting drug pressure by half or more, then there is an advantage to deploying two or more ACTs.
This analysis suggests that MFTs are deployed to best effect on different segments of the same population. One possibility would be to distribute one ACT for home-based care and use a different ACT in the clinic. Another promising strategy would partition the market by age: paediatric patients would be given one ACT, adults another. The paediatric and adult markets are already partitioned because adults and children require different formulations. Implementing MFTs would introduce operational challenges, but choosing different first-line ACT for paediatric and adult populations would work as a solution within a set of operational challenges that already exist.
MFTs are deployed to best effect before resistance emerges, but such decisions must be made in the face of massive uncertainty. Despite uncertainty, a strategic policy decision to protect artemisinin by promoting universal use of ACT has been proposed and partially implemented . If mutation rates for artemisinin resistance were known, the waiting time to the global emergence of ACT resistance could be calculated with increased confidence, and it would then be possible to evaluate whether existing policies were adequate. Recently, stable artemisinin resistance has been demonstrated in a murine model  and artemisinin-resistant infections were described along the Thai-Cambodia border  where artemisinin-class drugs have been used in Southeast Asia for twenty years . While, the rate that artemisinin resistance spontaneously appears is apparently low, its emergence as usage rates dramatically increased raises the question of if it is low enough to preserve ACT until a new class of anti-malarial drugs is available? The answer will not be known with confidence until after ACT resistance has fully emerged. Meanwhile, uncertainty will remain high.
Surprisingly, this analysis arrived at a robust conclusion despite the epidemiological and biological uncertainties. In small and isolated populations, combination therapies could have a huge benefit by delaying the time to appearance, so the benefit of MFTs would be comparatively small. In the global human population at risk of malaria, the enormous benefits of combination therapies shrank because ACT resistance could emerge anywhere and spread everywhere, and the relative benefits of delaying establishment and spread using MFTs were comparatively large.
Our study suggested that, with wide confidence intervals, the waiting time to the global emergence of resistance to ACT was approximately a decade. MFTs would provide some insurance against this calamity by delaying emergence further and then slowing the spread of resistance. The same logic would support a triple combination--three drugs given simultaneously, resistance to which requires three independent mutations. A triple combination would almost certainly delay the appearance of resistance by a decade or more if the genetic basis for resistance to all three components were mutually independent, but no such therapy is available today.
In retrospect, the results reflect that ACT resistance is unlikely to appear in any small population, but it will appear with virtual certainty in a large one, thus slowing spread has a comparatively large benefit at the global scale. Variability in drug pressure among populations would likely depend on access to health care, the level of clinical immunity, and socioeconomic status. The populations of the world can be considered a mosaic of sources and sinks for resistance, depending on whether resistance is favoured or not. In such a landscape, MFTs will shift the balance and increase the proportion of sinks that would naturally contain ACT resistance. Combination therapies are designed to make appearance a truly rare event, but if they don't work as well as intended, MFT strategies provide an insurance policy by delaying emergence further and by limiting spread after emergence. MFTs are thus a natural complement to combination therapies.
Combination therapies for malaria were adopted in the face of uncertainty. The basis for deploying artemisinin combination therapies comes mainly from a priori arguments about mutation rates and from experience with other drugs, such as the success of the combination drug sulphadoxine-pyrimethamine over pyrimethamine alone. The case for MFTs must be made on the same basis as the case for combination therapies but with a greater degree of uncertainty. Since the loss of artemisinin would likely result in the loss of all ACTs, existing ACTs must be used in the best possible way. Despite the uncertainty, the analysis makes a strong a priori case for a coordinated global policy, like the global subsidy for ACTs, to promote heterogeneous drug policies to ensure the long-term therapeutic efficacy of ACTs.
List of abbreviations
Artemisinin combination therapy
Multiple first-line therapies
This research was supported by Grant #4481, "Shaping the Global Subsidy for Anti-malarial Drugs" from the Bill & Melinda Gates Foundation to Ramanan Laxminarayan, Resources for the Future. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Preliminary findings were presented at a conference, "Anti-malarial Drug Strategies: Getting the Most from Anti-malarial Drugs," Prestana Kruger Lodge, Kruger National Park, South Africa, March 30-April 3, 2008. The authors gratefully acknowledge all the participants of that conference for their feedback and suggestions. The authors also thank Maciej Boni for his comments on earlier drafts.
- Guerra CA, Snow RW, Hay SI: Defining the global spatial limits of malaria transmission in 2005. Adv Parasitol. 2006, 62: 157-179. 10.1016/S0065-308X(05)62005-2.PubMed CentralView ArticlePubMed
- Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, Kabaria CW, Manh BH, Elyazar IR, Brooker S, Smith DL, Moyeed RA, Snow RW: A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6: e1000048.PubMed CentralView ArticlePubMed
- White NJ: The role of anti-malarial drugs in eliminating malaria. Malar J. 2008, 7 (Suppl 1): S8-10.1186/1475-2875-7-S1-S8.PubMed CentralView ArticlePubMed
- Snow RW, Trape JF, Marsh K: The past, present and future of childhood malaria mortality in Africa. Trends Parasitol. 2001, 17: 593-597. 10.1016/S1471-4922(01)02031-1.View ArticlePubMed
- White N: Antimalarial drug resistance and combination chemotherapy. Philos Trans R Soc Lond B Biol Sci. 1999, 354: 739-749. 10.1098/rstb.1999.0426.PubMed CentralView ArticlePubMed
- Foster KR, Grundmann H: Do we need to put society first? The potential for tragedy in antimicrobial resistance. PLoS Med. 2006, 3: e29-10.1371/journal.pmed.0030029.PubMed CentralView ArticlePubMed
- Bonhoeffer S, Lipsitch M, Levin BR: Evaluating treatment protocols to prevent antibiotic resistance. Proc Natl Acad Sci USA. 1997, 94: 12106-12111. 10.1073/pnas.94.22.12106.PubMed CentralView ArticlePubMed
- Boni MF, Smith DL, Laxminarayan R: Benefits of using multiple first-line therapies against malaria. Proc Natl Acad Sci USA. 2008, 105: 14216-14221. 10.1073/pnas.0804628105.PubMed CentralView ArticlePubMed
- Luria SE, Delbruck M: Mutations of bacteria from virus sensitivity to virus resistance. Genetics. 1943, 28: 491-511.PubMed CentralPubMed
- Antia R, Regoes RR, Koella JC, Bergstrom CT: The role of evolution in the emergence of infectious diseases. Nature. 2003, 426: 658-661. 10.1038/nature02104.View ArticlePubMed
- Handel A, Regoes RR, Antia R: The role of compensatory mutations in the emergence of drug resistance. PLoS Comput Biol. 2006, 2: e137-10.1371/journal.pcbi.0020137.PubMed CentralView ArticlePubMed
- Andersson DI: The biological cost of mutational antibiotic resistance: any practical conclusions?. Curr Opin Microbiol. 2006, 9: 461-465. 10.1016/j.mib.2006.07.002.View ArticlePubMed
- Andersson DI, Levin BR: The biological cost of antibiotic resistance. Curr Opin Microbiol. 1999, 2: 489-493. 10.1016/S1369-5274(99)00005-3.View ArticlePubMed
- Laufer MK, Thesing PC, Eddington ND, Masonga R, Dzinjalamala FK, Takala SL, Taylor TE, Plowe CV: Return of chloroquine antimalarial efficacy in Malawi. N Engl J Med. 2006, 355: 1959-1966. 10.1056/NEJMoa062032.View ArticlePubMed
- Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, Djimde AA, Kouriba B, Taylor TE, Plowe CV: Reemergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi. J Infect Dis. 2003, 187: 1870-1875. 10.1086/375419.View ArticlePubMed
- White NJ: Delaying antimalarial drug resistance with combination chemotherapy. Parassitologia. 1999, 41: 301-308.PubMed
- Guerra CA, Gikandi PW, Tatem AJ, Noor AM, Smith DL, Hay SI, Snow RW: The limits and intensity of Plasmodium falciparum transmission: implications for malaria control and elimination worldwide. PLoS Med. 2008, 5: e38-10.1371/journal.pmed.0050038.PubMed CentralView ArticlePubMed
- Becker N: Estimation for discrete time branching processes with application to epidemics. Biometrics. 1977, 33: 515-522. 10.2307/2529366.View ArticlePubMed
- Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.PubMed CentralView ArticlePubMed
- Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.PubMed CentralView ArticlePubMed
- Escalante AA, Smith DL, Kim Y: The dynamics of mutations associated with anti-malarial drug resistance in Plasmodium falciparum. Trends Parasitol. 2009, 25: 557-563. 10.1016/j.pt.2009.09.008.PubMed CentralView ArticlePubMed
- Hay SI, Okiro EA, Gething PW, Patil AP, Tatem AJ, Guerra CA, Snow RW: Estimating the global clinical burden of Plasmodium falciparum malaria in 2007. PLoS Med. 2010, 7: e1000290-10.1371/journal.pmed.1000290.PubMed CentralView ArticlePubMed
- Paget-McNicol S, Saul A: Mutation rates in the dihydrofolate reductase gene of Plasmodium falciparum. Parasitology. 2001, 122: 497-505. 10.1017/S0031182001007739.View ArticlePubMed
- Laxminarayan R, Gelband H: A global subsidy: key to affordable drugs for malaria?. Health Aff (Millwood). 2009, 28: 949-961. 10.1377/hlthaff.28.4.949.View Article
- Afonso A, Hunt P, Cheesman S, Alves AC, Cunha CV, do Rosario V, Cravo P: Malaria parasites can develop stable resistance to artemisinin but lack mutations in candidate genes atp6 (encoding the sarcoplasmic and endoplasmic reticulum Ca2+ ATPase), tctp, mdr1, and cg10. Antimicrob Agents Chemother. 2006, 50: 480-489. 10.1128/AAC.50.2.480-489.2006.PubMed CentralView ArticlePubMed
- Noedl H, Socheat D, Satimai W: Artemisinin-resistant malaria in Asia. N Engl J Med. 2009, 361: 540-541. 10.1056/NEJMc0900231.View ArticlePubMed
- White NJ: Qinghaosu (artemisinin): the price of success. Science. 2008, 320: 330-334. 10.1126/science.1155165.View ArticlePubMed
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.