Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2004, p. 2116-2123, Vol. 48, No. 6
0066-4804/04/$08.00+0 DOI: 10.1128/AAC.48.6.2116-2123.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
and Qin Cheng3
Malaria and Scabies Group, Queensland Institute of Medical Research, Royal Brisbane Hospital, Brisbane, Queensland 4029,1 Australian Centre for International and Tropical Health and Nutrition, University of Queensland, Queensland 4072,2 Department of Drug Resistance and Diagnostics, Australian Army Malaria Institute, Gallipoli Barracks, Enoggera, Queensland 4051, Australia3
Received 27 October 2003/ Returned for modification 22 December 2003/ Accepted 5 February 2004
|
|
|---|
|
|
|---|
Pyrimethamine and sulfadoxine act synergistically to inhibit two enzymes important in the parasite's folate biosynthetic pathway, dihydrofolate reductase (DHFR) and dihydropteroate synthetase (DHPS) (13). Point mutations in the DHFR and DHPS genes confer resistance to pyrimethamine and sulfadoxine, respectively, with decreasing in vitro Plasmodium falciparum susceptibility related to the number of mutations in each gene (6, 30, 34). The same mutations have been linked to treatment failure in the clinical setting (5, 7, 20, 33); the presence of mutations in DHFR appear to be more important in causing treatment failure than DHPS mutations (5). Although the molecular basis for SP resistance is understood, the factors promoting the development and transmission of these mutants are less clear.
It has been suggested that drug pharmacokinetics (12), overusage of drugs (28), cross-resistance between drugs (14), and inadequate treatment through inappropriate prescription or administration, noncompliance, or poor absorption (28, 36) contribute to the development of resistance. The timing of treatment relative to the initiation of an immune response in the patient has also been hypothesized as important in developing resistance (10), as have host immunity and transmission level (13, 38). Although individual factors involved in the evolution of drug resistance have been identified, the relative importance of these factors has not been reported in a quantitative format. This paper considers the influence of drug dosage and the timing of treatment on the rate of SP treatment failure predicted by a simulation model of P. falciparum infection. It also quantifies the relative importance of these factors in the development of SP resistance resulting in treatment failure. The selective pressure exerted by the long half-life of SP is also considered relative to its role in promoting the development and spread of resistance.
|
|
|---|
The probability of parasites' surviving at various time points after treatment was determined by combining information from SP isobolograms (32) and dose-response curves (39). Isobolograms for parasites containing triple mutations (3M) in both the DHFR and DHPS genes (3M/3M) and for parasites containing a triple mutation in DHFR and the wild-type DHPS gene (3M/WT) were digitized. These isobolograms contain the 50% inhibitory concentration (IC50) values for various SP concentrations in the presence of 45 nM folic acid. To estimate the probability of parasites' surviving, the dose-response curve for pyrimethamine (24) was used to estimate scaling factors, from which an IC50 in the isobolograms could be converted to an ICx (where x = 0,..., 100). The results from three dose-response curves, two from wild-type parasites (24, 39) and one from parasites containing a single mutation in DHFR (24), were averaged to achieve the scaling factors. It was assumed that these scaling factors were also applicable to sulfadoxine.
To calculate the survival probabilities for mutation combinations other than 3M/3M and 3M/WT, the following were assumed: a single mutation at codon 108 in the DHFR gene results in a
50-fold increase in the IC50 relative to the wild-type parasite (3, 5, 6, 25, 26); a double mutation in DHFR results in a
170-fold increase in the IC50 relative to the wild type (4, 6, 22, 26); a triple mutation in the DHFR gene results in a
700-fold increase in the IC50 compared to the wild type (3, 6); and single, double, and triple mutations in DHPS result in
15-,
45-, and
800-fold increases relative to the wild type in the IC50 for sulfadoxine, respectively (22, 34). These values were used to scale the isobolograms to represent other DHFR-DHPS mutation combinations.
Plasma drug concentrations were estimated by assuming that the average maximum concentration of pyrimethamine and sulfadoxine in the plasma of an adult after a standard three-tablet dose of Fansidar (75 mg of pyrimethamine and 1,500 mg of sulfadoxine) is 2.58 µM and 612 µM, respectively (35). Between 80 and 90% of pyrimethamine and 90 and 95% of sulfadoxine bind to plasma proteins (17), leaving approximately 0.52 and 60 µM unbound pyrimethamine and sulfadoxine, respectively, to act on the parasites. The elimination half-life of pyrimethamine and sulfadoxine was assumed to be 95.5 and 184 h, respectively (35). With this information, the drug concentration at various time points following treatment was determined. Three treatment regimes were simulated: the recommended three-tablet dose of Fansidar that reaches 100% of the expected plasma concentration, the recommended three-tablet dose of Fansidar that reaches 80% of the expected plasma concentration (e.g., poor absorption of drug), and a two-tablet dose of Fansidar. Figure 1 illustrates how the parasite genotype and drug kinetics were combined to estimate the probability of parasites' surviving drug treatment.
![]() View larger version (40K): [in a new window] |
FIG. 1. Parasite survival in the presence of various SP concentrations. (A) The wild-type (WT/WT) and (B) the 3M DHFR/3M DHPS parasite genotypes. Symbols on each of the surfaces indicate the SP concentrations that can be expected within a host 0, 12, 24, 36, and 48 days following treatment with a standard three-tablet dose of SP.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Pathology of simulated P. falciparum infections caused by the El Limon and Santee Cooper parasite strains treated at three different time points
|
![]() View larger version (23K): [in a new window] |
FIG. 2. Predicted treatment failure rates for various parasite genotypes when infections are treated with (A) three SP tablets, (B) three SP tablets which achieve 80% of the expected plasma concentration, and (C) two SP tablets. Genotypes not included in the graph had no simulated treatment failures. Simulations were conducted with treatment at 8 (solid bars), 10 (striped bars), or 12 days (open bars) after the start of the blood-stage infection.
|
|
View this table: [in a new window] |
TABLE 2. Prevalence of recrudescent infections with >1% of parasites carrying an additional mutation compared to the infecting parasite and with treatment administered 8 days after the start of asexual infection
|
![]() View larger version (23K): [in a new window] |
FIG. 3. Effect of residual drug on reinfection rates and development of clinical malaria. The percentage of reinfections that were successful by wild-type parasites (WT/WT) and parasites with a double mutation in either the DHFR (2M/WT) or DHPS (WT/2M) gene following SP treatment is indicated by the solid, dotted, and dashed lines, respectively. The bars indicate the mean time until the appearance of symptoms requiring treatment following reinfections occurring between 2 and 32 days post-SP treatment.
|
![]() View larger version (21K): [in a new window] |
FIG. 4. Illustrative changes in the prevalence of parasite genotypes caused by SP use in (A) low-transmission regions ( 1.8 infectious bites/year) and (B) high-transmission regions ( 18 infectious bites/year). In each region, the ratio of genotypes prior to the introduction of SP treatment was set to 9,900:30:30:10:10:10:7:3 for WT/WT, 1M/WT, WT/1M, 1M/1M, 2M/WT, WT/2M, 2M/1M, and 2M/2M parasites, respectively. The 1M/1M genotype was omitted because its prevalence curve was indistinguishable from the curve of the WT/1M genotype parasites.
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Schematic illustration of the evolution of SP resistance. The fate of parasites following treatment with SP is indicated for each genotype. At low levels of mutation, there is no mechanism for advancement of mutation level. However, the presence of residual drug concentrations provides a strong selection pressure against the indicated genotypes (*). At intermediate mutation levels, early suboptimal treatment provides an environment where parasites with additional mutations may be selected for within the host, leading to treatment failure. At high mutation levels, parasites are often able to survive therapeutic doses of SP, causing treatment failure.
|
|
|
|---|
The results presented need to be viewed relative to the assumptions of the model. The first and probably most important consideration is that the model mimics a P. falciparum infection in a malaria-naïve host such as a child, visitor, or immigrant from a nonmalarious area. Since semi-immune individuals typically respond better to treatment than hosts not previously exposed to malaria (37), the model output would be expected to overestimate the treatment failure rate for semi-immune individuals. The second consideration relates to modeling a single clone of parasites. In areas of high transmission, multiple infections are likely; the model ignores the effect of parasite recombination, which can act to reduce the probability that resistance will be retained by the parasite (11). Lastly, the determination of the speed with which resistance develops is restricted to limited situations in which everyone within a population of malaria-naïve hosts is treated when they become sick. This may occur in situations where a population is relocated into a malarious area and closely monitored. In other situations where only a proportion of the people are treated, the selection pressure on parasites would be expected to be less than that indicated. As such, the model predictions reported here represent a worse-case scenario for a naïve population.
Since the identification of molecular markers for SP resistance, numerous field studies have assessed the correlation between genetic mutations and treatment failure. A number of studies have attributed SP treatment failure to infection with parasites having at least a double mutation in DHFR coupled with a single mutation in the DHPS gene (1, 5, 7, 15, 20). The model results presented here are in agreement with these findings. A more detailed comparison of treatment failure rates for specific parasite genotypes indicated that the model predictions for mutation levels of 2M/1M or above encompass reported treatment failure rates at 3 days and 7 days posttreatment (4, 19). The simulation results also suggest that even within a naïve host, processes such as antigenic variation and the corresponding immune response that develops to it during the course of an infection may play a role in reducing the treatment failure rate. This type of interaction may explain the apparent discrepancy between the proportion of parasite isolates that are resistant in vitro and the much reduced incidence of in vivo treatment failure (1).
The model's prediction that treatment failure rates increase with suboptimal dosing agrees with field data from Kenyan children (29). The same study reported that treatment failure rates were higher in individuals who had been treated with SP in the 5 weeks prior to the current infection (29). Our predictions suggest that the increase in treatment failure observed may be due to the different length of protection afforded against reinfection for different parasite genotypes, so that only resistant parasites can reinfect individuals shortly after SP treatment.
The model output indicates that for the parasite characteristics considered, treatment of infections caused by parasites having at least two mutations in DHFR prior to the triggering of the specific immune response (day 8 in the model) resulted in 100% treatment failure, whereas treatment after the stimulation of the specific immune response (either day 10 or 12 in the model) resulted in lower failure rates. This result appears to contradict the general belief that treatment of an infection at a smaller parasite burden (earlier time) reduces the chances of treatment failure. However, closer examination of the factors influencing the survival of parasites following treatment indicates that the infectious load is not the only factor pertinent to determining treatment failure. This is particularly true for infections caused by parasites having at least two mutations in DHFR, for which even optimal doses of SP are often not sufficient to kill 100% of parasites. In these situations, the parasite burden is severely reduced following treatment but gradually increases again as the drug concentrations within the host wane. This high treatment failure rate can be improved by the development of a specific immune response to the variant antigens that can mop up the parasites not killed by the drug combination. Therefore, treating infections early so as to have a smaller biomass is important for infections with parasites having no or low levels of mutation but is predicted to lead to increased treatment failure in infections caused by highly resistant parasites.
For the parasite characteristics assumed in this model, the specific immunity for most of the expressed var genes has been triggered by day 12, coinciding with the onset of symptoms, while no or few antibodies have been triggered by day 8. However, these values are an approximation only and may vary for different parasites depending on the rate of var gene switching, the number of parasites expressing a variant required to trigger the specific immunity, and also the pyrogenic threshold (dictating the onset of symptoms).
In the majority of simulations in which recrudescence occurred, the model predicted that the recrudescent parasites exhibited the same genotype as those in the initial infection. However, parasites with additional mutations could become prominent under certain circumstances. These results agree with in vitro data showing no significant difference between the mean IC50 values of pyrimethamine and sulfadoxine for paired blood samples taken from patients prior to SP treatment and from recrudescent infections post-SP treatment (18). Epidemiological data collected from malaria-infected children also indicated that in the majority of SP treatment failures, the recrudescent parasites had the same DHFR-DHPS genotype as the parasites in the initial infection (23; A. Nzila, personal communication).
Although the model is able to predict the frequency of parasites with additional mutations becoming dominant within a host, it is not able to estimate the probability of these newly developed parasites being transmitted through the mosquito to a new human host and subsequently becoming successfully established in a community. Recent reports demonstrate that gene flow rather than new mutation is responsible for the high level of resistance mutations in the parasite's DHFR and DHPS in Africa (27) and DHFR in southeast Asia (21). This suggests that the successful establishment of these new mutated parasites in the community is a rare event.
The model predicted that the long half-life of SP provides strong selective pressure against parasites carrying wild-type DHFR, potentially resulting in the rapid decline of the wild-type DHFR genotype in high-transmission areas. Such a process may account for the rapid decrease (from 18 to 4% in 12 months) in the prevalence of the wild-type DHFR genotype in Tanzania (16). The model output also suggested that the infrequent observation of sulfadoxine-resistant but pyrimethamine-sensitive parasites in the field (18) may result from the greater selection pressure against wild-type DHFR/mutant DHPS parasites compared to mutant DHFR/wild-type DHPS parasites.
Although the long half-life of SP exerts a selection pressure favoring mutated parasites, it also has the potential to reduce the infection and transmission rate by providing some protection against reinfection. The actual magnitude of the potential decrease in transmission rate is related to the prevalence of parasites carrying mutations in DHFR and DHPS and the proportion of the population receiving treatment and therefore having residual drug in their blood. Such a decrease in transmission rate could partially offset the selection pressure caused by the long elimination half-life of SP. The current model, which was designed to mimic the in-host dynamics of P. falciparum infections in a naïve individual, requires further development to explore these types of interactions. Only then can the overall impact of SP treatment within a human population be assessed.
Hastings et al. (12) proposed that the evolution of drug resistance could be split into two distinct phases: the transition from wild-type parasites to slightly mutated forms that are less sensitive to the drug, but are still killed by therapeutic concentrations and the conversion from low to high levels of mutation and the emergence of clinical resistance to the drug. It was hypothesized that the level of drug usage in the population was a fundamental factor during the first stage, while the proportion of infections treated within the population was the primary factor influencing the second stage in the evolution of drug resistance (12). The results reported here support this general hypothesis, but suggest that there is an additional intermediate stage in which suboptimal treatment of low-grade infections is important.
The results presented indicate that the best protocol for slowing the development of drug resistance to SP is optimal treatment of individuals after the development of clinical symptoms and subsequent confirmation of malaria. The presumptive treatment of malaria by health workers or through self-medication has the potential to increase the speed with which resistance develops for two reasons. First, presumptive treatment of illnesses thought to be malaria creates a situation in which a larger proportion of the population has drug in their blood. This acts to increase the selection pressure against wild-type parasites, although it may also result in a reduction in transmission. Second, presumptive treatment of nonimmune individuals who are not ill from malaria but do carry low-grade resistant parasitemia due to a developing malaria infection results in treatment being administered prior to the triggering of any specific immune response. This has the potential to cause treatment failure, possibly selecting for a more resistant parasite population within the host.
The speed with which drug resistance develops is dictated by the prevalence of parasites carrying mutations in DHFR and/or DHPS when the drug is first introduced. Therefore, resistance to SP would be expected to develop more rapidly in areas where drugs targeting the same active site as either pyrimethamine or sulfadoxine have been or are being used (2, 14). Although these drugs may not necessarily be used to treat malaria, their use will ultimately result in malaria parasites' being exposed to subtherapeutic doses, providing an environment promoting the development of resistance to a component of the SP combination. Examples of this include the use of sulfa drugs to treat bacterial infections (31) or the use of trimethoprim, which cross-reacts with pyrimethamine, as a prophylactic treatment for opportunistic infections in human immunodeficiency virus-infected individuals (14). Since drug resistance is a problem for many infectious diseases, a combined effort from disease specialists is required to devise an overall treatment strategy which best suits the needs of an individual country or region. In this way, the effective life of many drugs used to treat multiple diseases may be extended.
The simulation model used in this analysis of SP resistance is equally applicable to the investigation of the development of drug resistance to other antimalarial combinations. A prime candidate for such analysis is the chlorproguanil-dapsone combination. With the pharmacokinetic data specific for this drug combination, it would be interesting to speculate on the factors promoting resistance and explore the likelihood of treatment failure with chlorproguanil-dapsone in regions already experiencing SP failure. Such analysis would provide useful information on the likely success and effective life of this new drug combination.
|
|
|---|
rk,i
1, depending on the drug concentration within the blood at the time of infection,
![]() | (1) |
![]() | (2) |
Ignoring parasite recombination and assuming that every person who becomes sick is treated and has an equal opportunity to transmit parasites to a susceptible mosquito, the change in the prevalence of genotypes can be tracked over time by using equation 3. The superscripts U and T represent the values of previously defined variables (S and r) in the untreated and treated populations, respectively. A time period representing 48 days is used in these calculations.
![]() | (3) |
![]() | (4) |
![]() | (5) |
![]() | (6) |
We thank Dennis Kyle for helpful discussions and valuable comments on the manuscript.
Present address: Malaria Vaccine Development Unit, NIAID, National Institutes of Health, Bethesda, Md. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»