- Open Access
Population pharmacokinetics of mefloquine, piperaquine and artemether-lumefantrine in Cambodian and Tanzanian malaria patients
Malaria Journal volume 12, Article number: 235 (2013)
Inter-individual variability in plasma concentration-time profiles might contribute to differences in anti-malarial treatment response. This study investigated the pharmacokinetics of three different forms of artemisinin combination therapy (ACT) in Tanzania and Cambodia to quantify and identify potential sources of variability.
Drug concentrations were measured in 143 patients in Tanzania (artemether, dihydroartemisinin, lumefantrine and desbutyl-lumefantrine), and in 63 (artesunate, dihydroartemisinin and mefloquine) and 60 (dihydroartemisinin and piperaquine) patients in Cambodia. Inter- and intra-individual variabilities in the pharmacokinetic parameters were assessed and the contribution of demographic and other covariates was quantified using a nonlinear mixed-effects modelling approach (NONMEM®).
A one-compartment model with first-order absorption from the gastrointestinal tract fitted the data for all drugs except piperaquine (two-compartment). Inter-individual variability in concentration exposure was about 40% and 12% for mefloquine. From all the covariates tested, only body weight (for all antimalarials) and concomitant treatment (for artemether only) showed a significant influence on these drugs’ pharmacokinetic profiles. Artesunate and dihydroartemisinin could not be studied in the Cambodian patients due to insufficient data-points. Modeled lumefantrine kinetics showed that the target day 7 concentrations may not be achieved in a substantial proportion of patients.
The marked variability in the disposition of different forms of ACT remained largely unexplained by the available covariates. Dosing on body weight appears justified. The concomitance of unregulated drug use (residual levels found on admission) and sub-optimal exposure (variability) could generate low plasma levels that contribute to selecting for drug-resistant parasites.
Artemisinin-based combination therapy (ACT) is the current first-line treatment of malaria . At present, the forms of ACT recommended by the World Health Organization (WHO) contain artemether (AM) plus lumefantrine (LF), artesunate (AS) plus either amodiaquine (AQ), mefloquine (MQ), pyronaridine (PN) or sulphadoxine-pyrimethamine (SP), and dihydroartemisinin (DHA) plus piperaquine (PPQ) . While hundreds of thousand courses of ACT are deployed each year , there is a limited number of studies measuring levels of drug exposure and relating it to treatment efficacy and safety.
Acquiring this information is paramount in order to optimize treatment and, especially, prevent resistance which may result from inadequate dosing. One of the main questions is to know if giving the recommended dose produces the same level of exposure in all, or otherwise what proportion or categories of subjects, and under which circumstances, would be systematically over- or under-dosed.
In other words, one needs to know if dosing regimens are adequate or if there are systematic dosing errors in which populations, especially on account of inter-subject variability and special groups like children and pregnant women. There is evidence that SP was systematically under-dosed in children and that the lower drug levels have contributed to the emergence of parasite resistance to this drug [3, 4]. The situation is further complicated by the fact that the target doses and therapeutic windows have been established based mostly on data in adults, and assume all patients require the same level of exposure, while, for instance, the contribution of immunity to parasite clearance will change with age and malaria transmission.
In this respect, one will need to know how the pharmacokinetics contributes to efficacy or safety outcomes. Examples of proposed surrogate efficacy correlates are day 7 drug plasma concentrations for LF [5–7] and the time for drug plasma concentrations to fall below 500 μg/l (the minimal inhibitory concentration, MIC) for MQ .
Treating with a wrong dose may have both individual and general consequences. Over-exposure increases the risk of toxicity; under-dosing may lead to treatment failure, but also carries the risk of selecting for drug-resistant parasites, which can spread to the rest of the population [9–12].
Two clinical trials found that AM-LF was highly efficacious in Tanzania, but much less effective (71% cure rate) in Cambodia [13, 14]. The cure rate in Cambodia increased to 86.5% in the subsequent years when 250 ml milk and coconut biscuits were provided with each dose of the study medication to increase drug absorption . These findings raised the question of the factors that could have contributed to the lower efficacy of AM-LF in the Cambodian population. Parasite susceptibility is indeed a potential explanation [15–24], although no known molecular marker exists at the moment. Another possibility is that differences in drug levels induced by genetic or other factors could explain the difference in drug response between these two populations.
The objectives of this paper were to characterize the population pharmacokinetics of AS, DHA, MQ and PPQ in Cambodian patients and AM and LF in Tanzanian patients and to identify demographic and other factors that could explain variability in drug levels. In addition, day 7 concentrations have been shown to be a good surrogate marker of treatment success and model based-simulations of LF were performed to predict the proportion of patients with concentrations below the proposed day 7 cut-off values.
Study area, patients and data
Three studies have been conducted, one in Tanzania and two in Cambodia. The study profiles are described in Figure 1. The first study was performed during March to May 2008 at the Kibaoni Health Centre, Kilombero district, Morogoro region, Tanzania. A total of 1,672 patients with suspected malaria were screened by rapid diagnostic test (Paracheck, Orchid Biomedical Systems, India) and 389 (23%) had parasitologically-proven Plasmodium falciparum malaria (confirmed and quantified by microscopy). After giving their informed consent to participate, patients were included in the study if they did not present signs of complicated malaria or any other severe concomitant illness. Six AM-LF (Coartem®, Novartis Pharma, Switzerland) doses were administered at time 0, 8, 20, 32, 44 and 56 h according to body weight (see Table 1); patients were either admitted for three days or asked to come back to the health facility for each drug administration. Mothers of breastfed patients were encouraged to feed their children and patients admitted were provided with food. Patients who reported that they have not eaten within two hours to prior dose intake were instructed to eat as soon as possible. Patients were seen by the clinical officer on days 0, 1, 2, 3, 7, 14, 28, and 42. A blood sample for pharmacokinetic measurement was taken at baseline (pre-dose on day 0) and on days 1, 2, and 7 at pre-defined random times after the last dosing. At each visit, a filter paper sample and a thin and thick smear were taken, in addition to axillary temperature, respiratory rate and haemoglobin measurements (only on days 0, 28 and 42) as well as evaluation of symptoms (e.g. headache, vomiting, and diarrhea). The exact dose and time of last drug intake, body weight, height, age, sex, food intake and concomitant medications were carefully recorded. If patients suffered from concomitant illnesses they were provided with additional treatment (paracetamol, mebendazole, metronidazole, cloxacillin, amoxicillin).
The second study was conducted during October 2007 to February 2008 at the Phnom Dék Health Centre, Rovieng district, Preah Vihear province, Cambodia. Entry criteria and study procedures were identical as in the Tanzanian study with minor adaptations described below. In total, 234 suspected malaria cases were screened by microscopy, of whom 67 (29%) were found to be infected with P. falciparum and 74 (32%) with Plasmodium vivax (no mixed infections were detected). Pregnant or lactating women were excluded. Three doses of AS (Arsumax®, Sanofi-Aventis, France) and MQ (Eloquine®, Medochemie Ltd, Cyprus) were given according to body weight on three consecutive days (see Table 1). Patients were seen by the clinical officer on days 1, 2, 3, 7, 14, 21, 28, 35 and 42. Sampling for the pharmacokinetic study was done at pre-dose and approximately 1 h after first dose intake on day 0, and on days 1, 2, 7 and 14 at pre-defined random times after drug intake.
The third study was performed during July to October 2008 at Pramoy Health Centre, Veal Veng district, Pursat province, Cambodia. A total of 287 suspected malaria cases were screened by microscopy, of whom 82 (29%) were infected with P. falciparum and 50 (17%) with P. vivax (no mixed infections). Children younger than six years of age and pregnant or lactating women were excluded. Three doses of DHA-PPQ (Duo-Cotecxin®, Zhejiang Holley Nanhu Pharmaceutical Co., Ltd, China) were given according to age (as per national guidelines) on three consecutive days (see Table 1). Same follow-up and blood sampling as at Phnom Dék were performed.
Samples of 2 ml of venous blood were collected on an EDTA Vacutainer and kept on ice for no longer than 6 h after withdrawal and then aliquoted into whole blood, plasma and pellet and immediately stored in liquid nitrogen or a –80°C freezer. Plasma concentrations of 14 antimalarial drugs and their metabolites, i.e. AM, AS, DHA, amodiaquine, N-desethyl-amodiaquine, LF, desbutyl-lumefantrine (DLF), PPQ, PN, MQ, chloroquine, quinine and SP, were determined simultaneously using a liquid chromatography–tandem mass spectrometry method (LC–MS/MS) previously reported . The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effect variability, overall process efficiency, standard addition experiments as well as anti-malarials short- and long-term stability in plasma. The method is precise (inter-day coefficient of variation: 3.1–12.6%) and sensitive (lower limits of quantification 0.15–3.0 for basic/neutral anti-malarials and 0.75–5 ng/ml for artemisinin derivatives, respectively). The laboratory is part of the quality control system of the World-Wide Antimalarial Resistance Network (WWARN).
Model-based pharmacokinetic analysis
The pharmacokinetic analysis for each drug taken separately was performed using the NONMEM computer program  Version 6 (NM-TRAN version II). It uses mixed (fixed and random) effects regression to estimate population means and variances of the pharmacokinetic parameters and to identify factors that influence them.
One-, two- and three-compartment pharmacokinetic models with first-order absorption, with and without absorption lag times, were compared. Additional one or two-compartments were used for anti-malarials presenting metabolite concentrations (AM and LF). The final parameters estimated were systemic clearance (CL/F), inter-compartmental clearance (Q/F), central volume of distribution (V C /F), peripheral volume of distribution (V P /F) and absorption rate constant (ka). Since no intravenous drug concentration data were available, these parameters represent apparent values. Where available, metabolite data were included into the model and metabolism rate constant from drug compartment to metabolite compartment (k 23 ) and metabolite clearance (CL met ) were also estimated. Owing to identifiability problems, the volume of distribution of the metabolites (V M ) DLF and DHA were assumed to equal LF and AM V C , respectively. Analysis of baseline plasma samples (i.e. day 0 prior treatment) showed that some patients had non-zero concentration of the drug, probably resulting from the treatment of the previous malaria episode or intake of non-declared drugs [29, 30]. The observed baseline residual plasma concentrations were fitted by estimating a factor (F0) that provided an estimation of the residual doses from previous treatment. A schematic representation of the models is presented in Figure 2.
Exponential errors following a log-normal distribution were assumed for the description of inter-patient variability of the pharmacokinetic parameters and were of the form: where θ i is the individual pharmacokinetic parameter value in the jth individual, θ is the population parameter estimate, and η j is the random effect value, which is independently and normally distributed with a mean of 0 and variance ω2. Proportional and combined proportional-and-additive error models were compared to describe intra-patient (residual) variability for the mother compound, and if available for its metabolite using C pij (1 + ϵ1ij) + ϵ2ij or C mij (1 + ϵ3ij) + ϵ4ij, where C pij are the corresponding predicted ith drug plasma concentration and C mij are the predicted metabolite concentration for the jth individual, and ϵ1ij, ϵ2ij, ϵ3ij, ϵ4ij, are independent normally distributed residual error terms with a mean of zero and a variance of σ21, σ22, σ23, σ24.
Available covariates were body weight, height, age, sex, smoking status, pregnancy (for Tanzanians patients), and concomitant medications. Reported concomitant medications were coded as moderate to strong inhibitors or inducers of the cytochrome P450 isoenzymes (CYP) mostly involved in the metabolism of the anti-malarials (Table 2) [31–40]. This information was based on report of self-medication prior inclusion and prescription during the study.
The covariate analysis was performed using a stepwise insertion/deletion approach. Visual inspection of the correlation between post hoc individual estimates of the pharmacokinetic parameters and the available covariates was first conducted by graphical exploration. Potentially influential covariates were then incorporated sequentially into the pharmacokinetic model. The typical value of a given parameter θ (e.g., CL) was modeled to depend on the covariate (X) either linearly θ = θa × [1 + θ b × X], exponentially θ = θ a X or as an allometric power function , with θ a representing the population value of the pharmacokinetic parameter and θ b the contribution of the covariate X, centered on the mean value; categorical covariates were coded as 0 or 1. In the allometric power models, θ b was either estimated or fixed to literature values, i.e. 0.75 for CL and 1 for V C . At the end of the analysis, all patient characteristics showing an influence on the parameters were again confirmed by comparing the full model (with all factors included) to models from which each of the factors was removed sequentially.
Model selection and parameter estimation
NONMEM®  (version 6.0, NM-TRAN, version II) was used with the FOCE INTERACTION method to fit the data. The difference in the minimum objective function value (Δ OFV) provided by NONMEM®, (–2 log likelihood, approximate χ2 distribution) was used to discriminate between models using the likelihood ratio test. A model was considered superior to another nested model when the OF value was reduced by at least 3.84 points (p < 0.05). Covariate analysis comprised forward selection of influential factors followed by backward deletion. Covariates were retained in the final model at the statistical level of p < 0.01. Model assessment was based on diagnostic plots (goodness-of-fit plots) along with standard errors and correlation matrix of parameter estimates, size of residual errors and eta-shrinkage.
The stability and the performance of the final population pharmacokinetic model were validated by the bootstrap method. Two hundred data sets were reconstructed by re-sampling from the original data using the Perl-speaks-NONMEM (PsN) Toolkit Version 3.2.4 [42, 43]. The final population pharmacokinetic model was fitted repeatedly to the 200 bootstrapped samples and pharmacokinetic parameters were calculated for each dataset. The mean, standard error and 95% confidence interval of each parameter obtained from the bootstrap analysis were then compared to the corresponding parameters obtained with the original dataset. The statistical analysis was performed using PsN version 3.2.4 . The final model was also validated using visual predictive check (VPC) obtained by simulation of data for 1’000 individuals based on the final model and generating 2.5th, 50th and 97.5th percentiles. The observed concentrations were plotted against the 95% prediction interval (P.I.95%) of the simulated dataset at each time point and visually compared. Figures were generated using GraphPad Prism (Version 4.00 for Windows, GraphPad Software, San Diego California USA ).
Model-based simulations for lumefantrine
Concentration-time profiles of lumefantrine in 1’000 individuals receiving two different 6-dose regimens over 3 days (doses at 0, 8, 20, 32, 44 and 56 h) and 5 days (doses at 0, 8, 24, 48, 72, 96 h) were performed based on the final model including inter-patient variability. These simulations served to purpose of quantifying the percentage of patients at day 7 below the different cut-off thresholds of 50 ng/ml, 175 ng/ml, 280 ng/ml and 600 ng/ml associated with treatment outcome . In addition, the simulation-based predicted median time (95% P.I.), estimated from time of last dose to 168 h (day 7), at which patients would exhibit concentrations below the cut-off values of 50 ng/ml, 175 ng/ml and 280 ng/ml was derived.
Population pharmacokinetic analyses
Patients’ baseline characteristics are summarized in Table 3. The median (range) of samples available per subject was 3 (2–4) for LF, 3 (1–4) for DLF, 2 (1–3) for AM and 2 (1–3) for DHA in Tanzania, 5 (3–6) for MQ, 1 for AS and 1 for DHA in Phnom Dék and 5 (4–6) for PPQ and 1 for DHA in Pramoy. The number of measured samples per time point is presented in Additional file 1.
A one-compartment model with first-order absorption from the gastrointestinal tract adequately described the data; no improvement was obtained with a two-compartment model (difference in the objective function (Δ OFV) = 0). For this drug, no baseline residual plasma concentrations were measured. In addition to CL, the assignment of an inter-patient variability on ka (but on no other parameter) significantly improved the description of the data (Δ OFV = −63). Univariate analyses showed that body weight (Δ OFV ≤ −45.1), age (Δ OFV ≤ −49.0), height (Δ OFV ≤ −42.7) and sex (Δ OFV ≤ −13.3) significantly influenced CL. In multivariate analyses, only body weight remained significant since all other variables were correlated to body weight. Linear and allometric power functions described the effect of body weight on CL similarly well (no statistical difference was observed between the two models (Δ OFV = −2)); the latter was finally chosen based on goodness of fit plots. The exponent of the allometric power function was estimated to be 0.66 and finally fixed to the literature value (0.75), since statistically not different (Δ OFV = 1). Inhibitors of CYP2C9 and/or CYP3A4 significantly influenced CL as well (Δ OFV = −7), indicating a 70% decrease in CL in patients exposed to either a CYP2C9 or CYP3A4 inhibitor. Multivariate analysis showed an additive influence of body weight and CYP inhibitors on CL (Δ OFV = −61 relative to the model without covariates).
Metabolite concentrations (DHA) were included in the model using an additional compartment, assuming linear metabolism and elimination. The assignment of an inter-patient variability on the metabolism rate constant k 23 yielded a better fit of the data (Δ OFV = −81), while no improvement was observed when assigning variability to the metabolite clearance CL met (Δ OFV = 0). Finally, none of the available covariates significantly affected DHA pharmacokinetics. A proportional error model for drug and metabolite provided the best description of intra-patient variability. The parameter estimates for the final model and derived parameters are in Table 4. The concentration-time plots of AM and DHA in the 135 patients included in the analysis with average population predictions and 95% prediction intervals is presented in Figure 3.
A one-compartment model with first-order absorption from the gastrointestinal tract and linear metabolism into DLF described adequately the data; a two-compartment model for LF or for DLF did not improve the model fit (Δ OFV = 0). The average estimated residual dose from previous treatments (F 0 ) was 1.6 mg, which corresponds to 0.3–1.3% of the recommended LF first dose (120–480 mg). Adding an inter-patient variability on V C (Δ OFV = −75), k23 (Δ OFV = −199) and F0 (Δ OFV = −17) in addition to CL improved the description of the data, but no variability on the other parameters was significant. A proportional error model best described the residual intra-patient variability for LF and an additive one for DLF. Inclusion of age, height and body weight on both CL and V C improved the fit (Δ OFV ≤ −30). Since age, height and body weight were correlated, only body weight was retained for further testing. Linear and allometric power functions adequately described its influence on CL and V C equally well (Δ OFV ≤ −34); the latter was selected based on visual inspection of graphical analysis. The estimations of the exponents of the allometric power functions were 0.52 and 0.35 for CL and V C , respectively, and provided a better fit than the fixed literature values (Δ OFV ≤ −17). Sex, smoking status, pregnancy and concomitant medications did not affect CL or V C (Δ OFV ≥ −0.2). The parameter estimates for the final model and derived parameters are given in Table 4. Figure 4 shows the concentration-time plots of LF and DLF in the 143 patients included in the analysis with average population predictions and 95% prediction intervals.
A one-compartment model with first-order absorption from the gastrointestinal tract appropriately described the data, with no improvement using a two-compartment model (Δ OFV = 0). For this drug, the residual dose from previous treatments (F0) was estimated to be 33.1 mg, corresponding to 6.7–26.7% of an initial dose of 125–500 mg. A better fit was obtained by assigning an inter-patient variability on V C (Δ OFV = −172) and F0 (Δ OFV = −211) in addition to CL. The use of a proportional error model for the residual intra-patient variability fitted the data well. Again, inclusion of age, height, body weight and sex improved the fit in univariate analyses (Δ OFV ≤ −136). Plots of CL and V C as a function of body weight suggested that an allometric power function, with exponent fixed to literature values, should be preferred to a linear relationship. The addition of smoking status and concomitant medications on CL and V C did not improve the model significantly (Δ OFV ≥ −1.2). Multivariate analysis indicated that body weight remained the unique significant covariate on both CL and V C . The parameter estimates for the final model and derived parameters are given in Table 4. Figure 5 depicts the simulated plasma concentration-time plot of MQ in the 63 patients included in the analysis with average population predictions and 95% prediction intervals.
A two-compartment model with first-order absorption from the gastrointestinal tract described the data better than a one-compartment model (Δ OFV = −97), but no additional benefit was seen with a three-compartment model (Δ OFV = 0). The residual dose of PPQ was estimated to be 123 mg, which corresponds to 12.8–25.6% of an initial dose of 480–960 mg. Assigning an inter-patient variability on V C (Δ OFV = −129) and V P (Δ OFV = −17) in addition to CL improved the fit and the use of a proportional error model for the residual intra-patient variability fitted the data adequately. CL and V C were again influenced by body weight (Δ OFV = −14 and −27, respectively); the relationship was best described using an allometric power function with exponents fixed to the literature values, and was not statistically different from estimated values (Δ OFV ≥ −2). Addition of sex or smoking status as covariates of CL did not improve the model fit (Δ OFV ≥ −2). As the metabolizing CYP of PPQ are not known and few concomitant treatments were reported, this variable was not included in the model. The parameter estimates for the final model and derived parameters are given in Table 4. Figure 6 shows the simulated plasma concentration-time plot of PPQ in the 60 patients included in the analysis with average population predictions and 95% prediction intervals.
Concentration-time simulations of lumefantrine
The day 7 predicted median concentrations of lumefantrine after administration of a 6 dose-regimen over 3 days were 300.9 ng/ml (P.I.95% 12.2−2015.0 ng/ml). Considering the large inter-patient variability in LF kinetics, 11% of the patients would exhibit day 7 concentrations below the cut-off of 50 ng/ml, 33% below 175 ng/ml, 48% below 280 ng/ml and 71% below 600 ng/ml. Prolonging the time of drug administration over 5 days would provide median concentrations of 608.7 ng/ml (P.I.95% 69.5−3515 ng/ml), with 1%, 10%, 21% and 49% of patients with concentrations below the cut-off of 50 ng/ml, 175 ng/ml, 280 ng/ml and 600 ng/ml, respectively. In addition, simulations predicted that patients would exhibit concentrations below the cut-off values of 50 ng/ml, 175 ng/ml and 280 ng/ml in a median (95% P.I.) of 152 h (126.8−176.3 h), 142 h (108.7−175.3 h) 136 h (99.7−172.3 h), respectively after a standard dosing regimen of 6 doses over 3 days. Increasing the 6-dose regimen over 5 days would increase the median time to 160 h (147.2−173.4 h), to 156 h (137.3−174.3 h) and to 152 h (128.2−174.9 h) for the 3 proposed cut-off values, respectively (Figure 7).
This study describes the disposition of three widely used forms of ACT (namely; AM-LF, AS-MQ and DHA-PPQ) and explores factors potentially influencing the marked variability in drug exposure. The estimated values of clearance and volume of distribution for AM (including its metabolite DHA), LF (including its desbutyl metabolite DLF), MQ and PPQ are in line with previously published results (see Table 5 for review), so are the large inter-patient and marked intra-individual variability [8, 32, 46–55]. Below some of the key findings are discussed.
Interestingly and worryingly, more than half of the patients had residual concentrations above the lower limit of quantification of at least one antimalarial on admission (74.3% in Tanzania Kibaoni, 51.6% in Cambodia Phnom Dék and 68.9% in Cambodia Pramoy). Residual doses were low in Tanzania for LF, but much higher values were estimated in Cambodia for MQ and PPQ, with up to approximately one quarter of an initial dose already present at baseline. The levels reflect the different residence times of these drugs; the proportions show the extent of unregulated drug use and selective pressure going on in these countries [29, 30].
Determinants of exposure – implications for dosing
There was a clear correlation between CL, V C and body weight for all drugs, which accounted for about 10−30% of the inter-individual variability in these two parameters. Body weight was highly correlated with age, sex and height and remained the only significant parameter in the multivariate analyses. This result supports the use of antimalarial dosing regimens based on body weight, or age as a proxy for it. What the model cannot predict is whether an additional correction to dosing should be made for children (or naïve adults) on account of lack of immunity.
The scaling factor of 0.75 for CL and 1 for V C described the relationships with body weight adequately, with the exception of LF for which the usual allometric scaling function provided a worse description of the data. Although a controversy persists regarding the body weight-dependent allometric exponent in the literature, it is not clear whether a different scaling between children and adults should be expected for this specific drug, or whether some confounding factors (different compliance between adult and children, different food intake) could have contributed to this finding.
Interactions and metabolism
Most of the anti-malarial drugs are metabolized by CYPs and concomitant treatment with inhibitors or inducers of these enzymes might, therefore, influence their elimination. This study detected only an influence of CYP2C9 and CYP3A4 inhibitors on AM clearance, which was decreased by 70% in patients with concomitant treatment. The fact that very few co-medications were reported might explain the lack of interactions for LF, MQ or PPQ. For the latter, the metabolizing pathways are not known. Among other factors, genetic polymorphisms in the enzymes responsible for antimalarial drugs could represent another important source of variability. A population genetic- and pharmacokinetic-based analysis was conducted to address this issue and published elsewhere .
Other conditions potentially influencing exposure
Pregnancy is known to lower blood concentrations of AM and LF, thus putting pregnant women at risk of under-dosing . This study enrolled only three pregnant women in Tanzania, which prevented estimating the influence of pregnancy on LF or AM drug levels.
Food intake has been shown to affect strongly the bioavailability of LF, MQ and PPQ [56–59]. In practice, this is a source of systematic under-dosing; in a study of the adherence to treatment regimens in Tanzania, only 0.4% of patients were reported to take their antimalarials with food . While in the present study patients admitted to the health facility were provided with food, outpatients were advised to eat directly before or after supervised drug-intake but adherence to this recommendation was based on self-reporting only, which, against the background of the above-mentioned adherence study, made it unreliable to have food effect included in the analysis. The absence of food information represents, however, a clear limitation of this study.
Structural and variance model
These drugs are known to exhibit multi-compartmental pharmacokinetics (two- or three-compartment disposition model) that could not be well captured owing to the limited duration of sampling compared to other studies. Although our data could provide a good estimation of CL and variability, appropriate description of terminal elimination phase could not be done. In addition, due to the very sparse sampling design during the absorption phase, no estimation of different absorption model could be performed, neither could the variability in the absorption quantified for LF, MQ and PPQ, with the exception of AM that exhibited a large inter-patient variability in its absorption. This large variability could result from both inherent characteristics of the drug [46, 61], and practical issues with dosing using crushed AM-LF pills (the paediatric formulation was not available for this study) . Another limitation of this study is that no estimation could be made for AS and DHA for the AS-MQ and DHA-PPQ treatments.
Simulations for LF
While the notion of concentration-effect relationship for LF is generally accepted, there is yet no common understanding of what the therapeutic target concentration should be. The published day 7 LF concentrations associated with therapeutic response range from 175 ng/ml to 600 ng/ml . A recent, large pooled analysis of LF concentration-efficacy data confirmed that a strong association exists between low day 7 LF concentrations and an increased risk of recurrence until day 42, and until day 21 for new infection . However, the authors of the pooled analysis concluded that there is no clear cut-off value for the thresholds associated with risk of recrudescence or new infection, but that cut-offs can be defined based on achieving a proportion of the desirable effect. For example, in low transmission areas a cut-off of 125 ng/ml gave efficacy rates of 84% and 96% at 42 days, and in high transmission areas a cut-off of 50 ng/ml gave efficacy rates of 80% and 95% at 42 days. In the Tanzanian sample of this study, 35% of the patients had a concentration below the cut-off value of 175 ng/ml, but only one of the 7 patients who had recurrent parasitaemia (unadjusted failure rate) was in this group.
Owing to the important variability in LF pharmacokinetics, the simulations under the standard 6-dose over 3 days schedule shows that a substantial proportion of the patients would present concentrations below the various proposed therapeutic targets at day 7. The considerable inter-individual variability in LF plasma concentrations additionally suggests that in some patients plasma LF concentrations would fall below the proposed minimal concentrations between the fourth and seventh day (3rd and 4th cycles) after treatment (i.e. before all the parasite had been eradicated). Splitting the same recommended total dose over 5 days would greatly reduce the probability of exhibiting sub-therapeutic drug concentrations, as already shown by other studies [32, 63]. However, in practice, the potential increased exposure with this 5-day regimen may be impeded by the possible risk of lower adherence to the treatment. Very little evidence exists for the other compounds. For MQ, the time over the MIC seems an important component associated with treatment efficacy . Our results indicate that this drug exhibits the least variability in its disposition and it is thus not expected that differences in response would be strongly related to variable drug levels.
The opinions expressed in this paper are those of the authors and may not reflect those of their employing organizations. PO is a staff member of the WHO; the authors alone are responsible for the views expressed in this publication and they do not necessarily represent the decisions, policy or views of the WHO.
Difference in the minimum objective function value
Typical value of given parameter in the population
Artemisinin-based combination therapy
95% confidence interval
Cytochrom P450 isoenzyme
Residual amount of drug from the previous treatment
Metabolism rate constant from drug compartment to metabolite compartment
Absorption rate constant
Minimal inhibitory concentration
95% prediction interval
Central volume of distribution
Volume of distribution of the metabolites
Peripheral volume of distribution
Visual predictive check
Steady-state volume of distribution
World Health Organization
World-Wide Antimalarial Resistance Network
WHO: Guidelines for the treatment of malaria. 2010, Geneva: World Health Organization, Second Edition.
WHO: World Malaria Report 2011. 2011, Geneva: World Health Organization
Terlouw DJ, Nahlen BL, Courval JM, Kariuki SK, Rosenberg OS, Oloo AJ, Kolczak MS, Hawley WA, Lal AA, Kuile FO: Sulfadoxine-pyrimethamine in treatment of malaria in Western Kenya: increasing resistance and underdosing. Antimicrob Agents Chemother. 2003, 47: 2929-2932. 10.1128/AAC.47.9.2929-2932.2003.
Barnes KI, Watkins WM, White NJ: Antimalarial dosing regimens and drug resistance. Trends Parasitol. 2008, 24: 127-134. 10.1016/j.pt.2007.11.008.
Barnes KI, Little F, Smith PJ, Evans A, Watkins WM, White NJ: Sulfadoxine-pyrimethamine pharmacokinetics in malaria: Pediatric dosing implications. Clin Pharmacol Ther. 2006, 80: 582-596. 10.1016/j.clpt.2006.08.016.
WHO: Methods and techniques for assessing exposure to antimalarial drugs in clinical field studies. 2011, Geneva: World Health Organization
Stepniewska K, Workman L, Aweeka F, Borrmann S, Faucher J-F, Genton B, Lindegardh N, Martensson A, Mayxay M, Nosten F, Piola P, Ringwald P, White NJ, Barnes K: Defining adequate lumefantrine exposure in patients with uncomplicated malaria treated with artemether-lumefantrine. Poster LB-2212, presented at the 60th American Society of Tropical Medicine and Hygiene (ASTMH) Annual Meeting 3-8 December 2011, Philadelphia, PA, USA. 2011
Simpson JA, Price R, ter Kuile F, Teja-Isavatharm P, Nosten F, Chongsuphajaisiddhi T, Looareesuwan S, Aarons L, White NJ: Population pharmacokinetics of mefloquine in patients with acute falciparum malaria. Clin Pharmacol Ther. 1999, 66: 472-484. 10.1016/S0009-9236(99)70010-X.
Barnes KI, Lindegardh N, Ogundahunsi O, Olliaro P, Plowe CV, Randrianarivelojosia M, Gbotosho GO, Watkins WM, Sibley CH, White NJ: World antimalarial resistance network (WARN) IV: Clinical pharmacology. Malar J. 2007, 6: 122. 10.1186/1475-2875-6-122.
Sevene E, Gonzalez R, Menendez C: Current knowledge and challenges of antimalarial drugs for treatment and prevention in pregnancy. Expert Opin Pharmacother. 2010, 11: 1277-1293. 10.1517/14656561003733599.
White NJ, Stepniewska K, Barnes K, Price RN, Simpson J: Simplified antimalarial therapeutic monitoring: using the day-7 drug level?. Trends Parasitol. 2008, 24: 159-163. 10.1016/j.pt.2008.01.006.
Watkins WM, Mosobo M: Treatment of Plasmodium falciparum malaria with pyrimethamine-sulfadoxine: selective pressure for resistance is a function of long elimination half-life. Trans R Soc Trop Med Hyg. 1993, 87: 75-78. 10.1016/0035-9203(93)90431-O.
Denis MB, Tsuyuoka R, Lim P, Lindegardh N, Yi P, Top SN, Socheat D, Fandeur T, Annerberg A, Christophel EM, Ringwald P: Efficacy of artemether-lumefantrine for the treatment of uncomplicated falciparum malaria in northwest Cambodia. Trop Med Int Health. 2006, 11: 1800-1807. 10.1111/j.1365-3156.2006.01739.x.
Kabanywanyi AM, Mwita A, Sumari D, Mandike R, Mugittu K, Abdulla S: Efficacy and safety of artemisinin-based antimalarial in the treatment of uncomplicated malaria in children in southern Tanzania. Malar J. 2007, 6: 146. 10.1186/1475-2875-6-146.
Antimalarial drug resistance, Thai-Cambodian border. Wkly Epidemiol Rec. 2009, 84: 94-95.
Resistance to artemisinin derivatives along the Thai-Cambodian border. Wkly Epidemiol Rec. 2007, 82: 360.
Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, Ringwald P, Silamut K, Imwong M, Chotivanich K, Lim P, Herdman T, An SS, Yeung S, Singhasivanon P, Day NP, Lindegardh N, Socheat D, White NJ: Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009, 361: 455-467. 10.1056/NEJMoa0808859.
Burki T: Artemisinin resistance could endanger fight against malaria. Lancet Infect Dis. 2009, 9: 213.
Enserink M: Malaria. Signs of drug resistance rattle experts, trigger bold plan. Science. 2008, 322: 1776. 10.1126/science.322.5909.1776.
Duffy PE, Sibley H: Are we losing artemisinin combination therapy already?. Lancet. 2005, 366: 1908-1909. 10.1016/S0140-6736(05)67768-9.
Maude RJ, Pontavornpinyo W, Saralamba S, Aguas R, Yeung S, Dondorp AM, Day NPJ, White NJ, White LJ: The last man standing is the most resistant: eliminating artemisinin-resistant malaria in Cambodia. Malar J. 2009, 8: 31. 10.1186/1475-2875-8-31.
Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM, Artemisinin Resistance in Cambodia 1 (ARC1) Study Consortium: Evidence of artemisinin-resistant malaria in Western Cambodia. N Engl J Med. 2008, 359: 2619-2620. 10.1056/NEJMc0805011.
Samarasekera U: Countries race to contain resistance to key antimalarial. Lancet. 2009, 374: 277-280. 10.1016/S0140-6736(09)61349-0.
WHO: Drug resistance could set back malaria control success. 2009, World Health Organization,http://www.who.int/mediacentre/news/releases/2009/malaria_drug_resistance_20090225/en/,
Hodel EM, Zanolari B, Mercier T, Biollaz J, Keiser J, Olliaro P, Genton B, Decosterd LA: A single LC-tandem mass spectrometry method for the simultaneous determination of 14 antimalarial drugs and their metabolites in human plasma. J Chromatogr B Analyt Technol Biomed Life Sci. 2009, 877: 867-886. 10.1016/j.jchromb.2009.02.006.
Hodel EM, Ley SD, Qi W, Ariey F, Genton B, Beck HP: A microarray-based system for the simultaneous analysis of single nucleotide polymorphisms in human genes involved in the metabolism of anti-malarial drugs. Malar J. 2009, 8: 285. 10.1186/1475-2875-8-285.
Staehli Hodel EM, Csajka C, Ariey F, Guidi M, Kabanywanyi AM, Duong S, Decosterd LA, Olliaro P, Beck HP, Genton B: Effect of Single Nucleotide Polymorphisms in Cytochrome P450 Isoenzyme and N-Acetyltransferase 2 Genes on the Metabolism of Artemisinin-Based Combination Therapies in Malaria Patients from Cambodia and Tanzania. Antimicrob Agents Chemother. 2013, 57: 950-958. 10.1128/AAC.01700-12.
NONMEM Users Guides (1989-2006). Edited by: Beal SL, Sheiner LB, Boeckmann AJ. 2006, Ellicott City, Maryland, USA: Icon Development Solutions
Hodel EM, Genton B, Zanolari B, Mercier T, Duong S, Beck HP, Olliaro P, Decosterd LA, Ariey F: Residual antimalarial concentrations before treatment in patients with malaria from Cambodia: indication of drug pressure. J Infect Dis. 2010, 202: 1088-1094. 10.1086/655779.
Hodel EM, Kabanywanyi AM, Malila A, Zanolari B, Mercier T, Beck HP, Buclin T, Olliaro P, Decosterd LA, Genton B: Residual antimalarials in malaria patients from Tanzania–implications on drug efficacy assessment and spread of parasite resistance. PLoS One. 2009, 4: e8184. 10.1371/journal.pone.0008184.
Tarning J, Bergqvist Y, Day NP, Bergquist J, Arvidsson B, White NJ, Ashton M, Lindegardh N: Characterization of human urinary metabolites of the antimalarial piperaquine. Drug Metab Dispos. 2006, 34: 2011-2019. 10.1124/dmd.106.011494.
White NJ, van Vugt M, Ezzet F: Clinical pharmacokinetics and pharmacodynamics of artemether-lumefantrine. Clin Pharmacokinet. 1999, 37: 105-125. 10.2165/00003088-199937020-00002.
Navaratnam V, Mansor SM, Sit NW, Grace J, Li QG, Olliaro P: Pharmacokinetics of artemisinin-type compounds. Clin Pharmacokinet. 2000, 39: 255-270. 10.2165/00003088-200039040-00002.
Novartis Pharma Corporation: Coartem. Full prescribing information. 2009
Khoo S, Back D, Winstanley P: The potential for interactions between antimalarial and antiretroviral drugs. AIDS. 2005, 19: 995-1005. 10.1097/01.aids.0000174445.40379.e0.
Ilett KF, Ethell BT, Maggs JL, Davis TME, Batty KT, Burchell B, Binh TQ, Thu LTA, Hung NC, Pirmohamed M, Park BK, Edwards G: Glucuronidation of dihydroartemisinin in vivo and by human liver microsomes and expressed UDP-glucuronosyltransferases. Drug Metab Dispos. 2002, 30: 1005-1012. 10.1124/dmd.30.9.1005.
Fontaine F, de Sousa G, Burcham PC, Duchene P, Rahmani R: Role of cytochrome P450 3A in the metabolism of mefloquine in human and animal hepatocytes. Life Sci. 2000, 66: 2193-2212. 10.1016/S0024-3205(00)00546-4.
Ridtitid W, Wongnawa M, Mahatthanatrakul W, Chaipol P, Sunbhanich M: Effect of rifampin on plasma concentrations of mefloquine in healthy volunteers. J Pharm Pharmacol. 2000, 52: 1265-1269. 10.1211/0022357001777243.
Ridtitid W, Wongnawa M, Mahatthanatrakul W, Raungsri N, Sunbhanich M: Ketoconazole increases plasma concentrations of antimalarial mefloquine in healthy human volunteers. J Clin Pharm Ther. 2005, 30: 285-290. 10.1111/j.1365-2710.2005.00651.x.
van Agtmael MA, van der Graaf CAA, Dien TK, Koopmans RP, van Boxtel CJ: The contribution of the enzymes CYP2D6 and CYP2C19 in the demethylation of artemether in healthy subjects. Eur J Drug Metab Pharmacokinet. 1998, 23: 429-436. 10.1007/BF03192305.
Anderson BJ, Holford NHG: Mechanism-based concepts of size and maturity in pharmacokinetics. Annu Rev Pharmacol Toxicol. 2008, 48: 303-332. 10.1146/annurev.pharmtox.48.113006.094708.
Lindbom L, Pihlgren P, Jonsson EN: PsN-Toolkit--a collection of computer intensive statistical methods for non-linear mixed effect modeling using NONMEM. [Erratum appears in Comput Methods Programs Biomed. 2005 Dec;80(3):277 Note: Jonsson, Niclas [corrected to Jonsson, E Niclas]]. Comput Methods Programs Biomed. 2005, 79: 241-257. 10.1016/j.cmpb.2005.04.005.
Lindbom L, Ribbing J, Jonsson EN: Perl-speaks-NONMEM (PsN)–a Perl module for NONMEM related programming. Comput Methods Programs Biomed. 2004, 75: 85-94. 10.1016/j.cmpb.2003.11.003.
Hietala SF, Martensson A, Ngasala B, Dahlstrom S, Lindegardh N, Premji Z, Farnert A, Gil P, Bjorkman A, Ashton M: Population pharmacokinetics and pharmacodynamics of artemether and lumefantrine during combination treatment in children with uncomplicated falciparum malaria in Tanzania. Antimicrob Agents Chemother. 2010, 54: 4780-4788. 10.1128/AAC.00252-10.
Ezzet F, Mull R, Karbwang J: Population pharmacokinetics and therapeutic response of CGP56697 (artemether plus benflumetol) in malaria patients. Br J Clin Pharmacol. 1998, 46: 553-561.
Salman S, Page-Sharp M, Griffin S, Kose K, Siba PM, Ilett KF, Mueller I, Davis TM: Population pharmacokinetics of artemether, lumefantrine, and their respective metabolites in Papua New Guinean children with uncomplicated malaria. Antimicrob Agents Chemother. 2011, 55: 5306-5313. 10.1128/AAC.05136-11.
Tarning J, Kloprogge F, Piola P, Dhorda M, Muwanga S, Turyakira E, Nuengchamnong N, Nosten F, Day NP, White NJ, Guerin PJ, Lindegardh N: Population pharmacokinetics of artemether and dihydroartemisinin in pregnant women with uncomplicated Plasmodium falciparum malaria in Uganda. Malar J. 2012, 11: 293. 10.1186/1475-2875-11-293.
Tarning J, McGready R, Lindegardh N, Ashley EA, Pimanpanarak M, Kamanikom B, Annerberg A, Day NPJ, Stepniewska K, Singhasivanon P, White NJ, Nosten F: Population pharmacokinetics of lumefantrine in pregnant women treated with artemether-lumefantrine for uncomplicated Plasmodium falciparum malaria. Antimicrob Agents Chemother. 2009, 53: 3837-3846. 10.1128/AAC.00195-09.
Charles BG, Blomgren A, Nasveld PE, Kitchener SJ, Jensen A, Gregory RM, Robertson B, Harris IE, Reid MP, Edstein MD: Population pharmacokinetics of mefloquine in military personnel for prophylaxis against malaria infection during field deployment. Eur J Clin Pharmacol. 2007, 63: 271-278. 10.1007/s00228-006-0247-3.
Ashley EA, Stepniewska K, Lindegardh N, McGready R, Hutagalung R, Hae R, Singhasivanon P, White NJ, Nosten F: Population pharmacokinetic assessment of a new regimen of mefloquine used in combination treatment of uncomplicated falciparum malaria. Antimicrob Agents Chemother. 2006, 50: 2281-2285. 10.1128/AAC.00040-06.
Tarning J, Ashley EA, Lindegardh N, Stepniewska K, Phaiphun L, Day NPJ, McGready R, Ashton M, Nosten F, White NJ: Population pharmacokinetics of piperaquine after two different treatment regimens with dihydroartemisinin-piperaquine in patients with Plasmodium falciparum malaria in Thailand. Antimicrob Agents Chemother. 2008, 52: 1052-1061. 10.1128/AAC.00955-07.
Tarning J, Zongo I, Some FA, Rouamba N, Parikh S, Rosenthal PJ, Hanpithakpong W, Jongrak N, Day NP, White NJ, Nosten F, Ouedraogo JB, Lindegardh N: Population pharmacokinetics and pharmacodynamics of piperaquine in children with uncomplicated falciparum malaria. Clin Pharmacol Ther. 2012, 91: 497-505. 10.1038/clpt.2011.254.
Tarning J, Rijken MJ, McGready R, Phyo AP, Hanpithakpong W, Day NP, White NJ, Nosten F, Lindegardh N: Population pharmacokinetics of dihydroartemisinin and piperaquine in pregnant and nonpregnant women with uncomplicated malaria. Antimicrob Agents Chemother. 2012, 56: 1997-2007. 10.1128/AAC.05756-11.
Crevoisier C, Handschin J, Barre J, Roumenov D, Kleinbloesem C: Food increases the bioavailability of mefloquine. Eur J Clin Pharmacol. 1997, 53: 135-139. 10.1007/s002280050351.
Hoffmann-La Roche Ltd: Lariam. Full prescribing information. 2009
Price RN, Uhlemann AC, van Vugt M, Brockman A, Hutagalung R, Nair S, Nash D, Singhasivanon P, Anderson TJC, Krishna S, White NJ, Nosten F: Molecular and pharmacological determinants of the therapeutic response to artemether-lumefantrine in multidrug-resistant Plasmodium falciparum malaria. Clin Infect Dis. 2006, 42: 1570-1577. 10.1086/503423.
Sim IK, Davis TME, Ilett KF: Effects of a high-fat meal on the relative oral bioavailability of piperaquine. Antimicrob Agents Chemother. 2005, 49: 2407-2411. 10.1128/AAC.49.6.2407-2411.2005.
Kabanywanyi AM, Lengeler C, Kasim P, King'eng'ena S, Schlienger R, Mulure N, Genton B: Adherence to and acceptability of artemether-lumefantrine as first-line anti-malarial treatment: evidence from a rural community in Tanzania. Malar J. 2010, 9: 48. 10.1186/1475-2875-9-48.
Ali S, Najmi MH, Tarning J, Lindegardh N: Pharmacokinetics of artemether and dihydroartemisinin in healthy Pakistani male volunteers treated with artemether-lumefantrine. Malar J. 2010, 9: 275. 10.1186/1475-2875-9-275.
Abdulla S, Sagara I, Borrmann S, D'Alessandro U, Gonzalez R, Hamel M, Ogutu B, Martensson A, Lyimo J, Maiga H, Sasi P, Nahum A, Bassat Q, Juma E, Otieno L, Björkman A, Beck HP, Andriano K, Cousin M, Lefèvre G, Ubben D, Premji Z: Efficacy and safety of artemether-lumefantrine dispersible tablets compared with crushed commercial tablets in African infants and children with uncomplicated malaria: a randomised, single-blind, multicentre trial. Lancet. 2008, 372: 1819-1827. 10.1016/S0140-6736(08)61492-0.
Ezzet F, van Vugt M, Nosten F, Looareesuwan S, White NJ: Pharmacokinetics and pharmacodynamics of lumefantrine (benflumetol) in acute falciparum malaria. Antimicrob Agents Chemother. 2000, 44: 697-704. 10.1128/AAC.44.3.697-704.2000.
The authors would like to thank the patients who participated in the study. Furthermore, the authors would like to thank the staff from the National Center for Parasitology, Entomology and Malaria Control in Cambodia, namely Dr. Leang Rithea, Dr. Mey Bouth Denis, Dr. Chivv Lim, Va Soch, and Oung Chavvin, the staff from Phnom Dék and Pramoy Health Centre in Cambodia, the staff from the Ifakara Health Institute in Tanzania, namely Dr. Hassan Mshinda, Aggrey Malila and Fidelis Mbena, the staff of the Kibanoi Health Centre in Tanzania and all the drivers for their help in the sample collection. The computations were performed at the Vital-IT (http://www.vital-it.ch) Center for high-performance computing of the Swiss Institute of Bioinformatics and the authors are grateful for their support.
Support for this work was provided by grant No 320000-112479 from the Swiss National Science Foundation. The authors declare that they have no competing interests.
EMSH contributed to the design of the study, the acquisition of patient samples and laboratory data, the analysis and interpretation of data and to drafting of the manuscript. DS and AMK made substantial contributions to the acquisition of patient samples and the coordination of the study. FA was involved in the acquisition of patient samples and the assay development for the pharmacogenetic analysis. BZ, TM and LAD where involved in the assay development and in the quantitation of the drug concentration levels. MG and CC carried out the population pharmacokinetic analysis and made substantial contributions to interpretation of data and to drafting and revising the manuscript. TB contributed to interpretation of pharmacokinetic data and the revising of the manuscript. BG, HPB and PO conceived the study, participated in its design and coordination, helped to draft the manuscript and revised it critically for important intellectual content. All authors have given final approval of the version to be published.
Blaise Genton and Chantal Csajka contributed equally to this work.
Electronic supplementary material
Additional file 1: Number of samples per time point. The table provided summarizes the number of patient samples included in the population pharmacokinetic model of each anti-malarial drug (and its metabolite where applicable) for every sampling time point. (PDF 61 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Staehli Hodel, E.M., Guidi, M., Zanolari, B. et al. Population pharmacokinetics of mefloquine, piperaquine and artemether-lumefantrine in Cambodian and Tanzanian malaria patients. Malar J 12, 235 (2013). https://doi.org/10.1186/1475-2875-12-235
- Population pharmacokinetics
- Artemisinin-based combination therapy
- Nonlinear mixed-effects modelling