- Open Access
A population pharmacokinetic model of piperaquine in pregnant and non-pregnant women with uncomplicated Plasmodium falciparum malaria in Sudan
Malaria Journal volume 11, Article number: 398 (2012)
Pregnancy is associated with an increased risk of developing a malaria infection and a higher risk of developing severe malaria. The pharmacokinetic properties of many anti-malarials are also altered during pregnancy, often resulting in a decreased drug exposure. Piperaquine is a promising anti-malarial partner drug used in a fixed-dose combination with dihydroartemisinin. The aim of this study was to investigate the population pharmacokinetics of piperaquine in pregnant and non-pregnant Sudanese women with uncomplicated Plasmodium falciparum malaria.
Symptomatic patients received a standard dose regimen of the fixed dose oral piperaquine-dihydroartemisinin combination treatment. Densely sampled plasma aliquots were collected and analysed using a previously described LC-MS/MS method. Data from 12 pregnant and 12 non-pregnant women were analysed using nonlinear mixed-effects modelling. A Monte Carlo Mapped Power (MCMP) analysis was conducted based on a previously published study to evaluate the power of detecting covariates in this relatively small study.
A three-compartment disposition model with a transit-absorption model described the observed data well. Body weight was added as an allometric function on all clearance and volume parameters. A statistically significant decrease in estimated terminal piperaquine half-life in pregnant compared with non-pregnant women was found, but there were no differences in post-hoc estimates of total piperaquine exposure. The MCMP analysis indicated a minimum of 13 pregnant and 13 non-pregnant women were required to identify pregnancy as a covariate on relevant pharmacokinetic parameters (80% power and p=0.05). Pregnancy was, therefore, evaluated as a categorical and continuous covariate (i.e. estimate gestational age) in a full covariate approach. Using this approach pregnancy was not associated with any major change in piperaquine elimination clearance. However, a trend of increasing elimination clearance with increasing gestational age could be seen.
The population pharmacokinetic properties of piperaquine were well described by a three-compartment disposition model in pregnant and non-pregnant women with uncomplicated malaria. The modelling approach showed no major difference in piperaquine exposure between the two groups and data presented here do not warrant a dose adjustment in pregnancy in this vulnerable population.
Malaria is one of the most important infectious diseases with an estimated 216 million people infected worldwide in 2010 . Pregnant women are at an increased risk of developing a malaria infection  and are at a higher risk of progressing to severe malaria [3–5]. Plasmodium falciparum malaria is a major contributor to maternal mortality in Sudan; around 37% of all maternal deaths between 1985 and 1999 at the Medani Teaching Hospital in Medani City, Sudan, were attributed to malaria . Malaria also has severe effects on the foetus causing both foetal loss and low birth weight.
Artemisinin-based combination therapy (ACT) is recommended as first-line treatment for P. falciparum malaria in all endemic areas. The artemisinin derivatives have a very rapid parasiticidal effect, which substantially reduces the parasite biomass during the first days of treatment. These drugs have a short terminal elimination half-life and are, therefore, used in combination with longer acting anti-malarials, with the aim of preventing recrudescence by killing residual parasites. Combination therapies consisting of drugs with different mechanisms of action also reduce the risk of the development of drug resistance [7, 8].
The oral fixed-dose combination of dihydroartemisinin and piperaquine has shown excellent efficacy in the treatment of P. falciparum malaria [9–13]. Piperaquine is highly bound to plasma proteins (>99.9%), has a large apparent volume of distribution, (103–874 L/kg), a low apparent elimination clearance (0.6-1.3 L/h/kg) and, therefore, a long terminal elimination half-life (12–28 days) [14–20].
Pregnancy has considerable effects on the pharmacokinetic properties of many drugs. Renal elimination, expression of metabolising enzymes, volume of body water and the degree of plasma protein binding all change during pregnancy [21–23]. This could result in lower drug plasma concentrations . Previously published studies have reported a decrease in drug exposure during the later stages of pregnancy for artesunate, artemether, dihydroartemisinin, lumefantrine, sulphadoxine, atovaquone, proguanil and cycloguanil [19, 25–31]. Other anti-malarial drugs (e.g. quinine and amodiaquine/desethylamodiaquine) show no differences in pharmacokinetic properties in pregnant women compared to non-pregnant women [32–34].
Only one previous study has investigated the impact of pregnancy on the pharmacokinetic properties of piperaquine in patients with uncomplicated P. falciparum malaria . Pregnancy was found to affect the elimination clearance and the bioavailability of piperaquine, but with no change in total drug exposure. This was further supported by a non-compartmental analysis of the same study . No published information is available on the population pharmacokinetic properties of piperaquine in pregnant or non-pregnant women in an African country.
The aim of this study was to describe the population pharmacokinetic properties of piperaquine in pregnant and non-pregnant women with uncomplicated P. falciparum malaria in Sudan.
The study was conducted at the New Halfa Teaching Hospital in New Halfa, Sudan. Clinical details and non-compartmental analysis results are reported in full elsewhere . The participating women received a written and oral explanation of the study in their own language. If the woman could not read, the explanation was read to her. Ethics approval for the study was given by the College of Medical and Technical Studies, Khartoum, Sudan.
Symptomatic pregnant women with uncomplicated P. falciparum malaria in their second or third trimester attending the antenatal clinic in New Halfa were eligible to participate in the study. Non-pregnant women with uncomplicated P. falciparum malaria were also recruited as controls.
All patients received dihydroartemisinin-piperaquine tetra-phosphate (Duo Cotecxin, 40 mg/320 mg tablets, Beijing Holley-Cotec Pharmaceuticals, Co., Ltd.) once daily for three days. Drug administration was directly observed and taken with a glass of water under fasting conditions. The number of tablets was based on the patient’s body weight to achieve a daily dose of 20 mg piperaquine tetra-phosphate/kg.
Blood samples were obtained by venous puncture or a three-way tap. PCR, haematology and biochemistry samples were drawn before the first dose and on day 14. Blood samples (2 mL) for pharmacokinetic analysis were drawn pre-dose and at 1.5, 4, 8, 24, 25.5, 28, 32, 48, 49, 50, 52, 56, 60, 72 h after the first dose and on days 5, 7, 14, 21, 28, 35, 42, 49, 56, 63 and 90. The actual time of dosing and sampling were noted and used in the pharmacokinetic analysis. Blood samples were centrifuged at 2000×g for 10 minutes and plasma samples stored in liquid nitrogen until the samples were transferred to Khartoum there they were stored in −80°C.
The chemical analysis was performed using a previously published method with separation and quantification by liquid chromatography (LC) and tandem mass spectrometry (MS/MS) detection . The LC-system was an Agilent 1200 system consisting of a binary LC pump, a vacuum degasser, autosampler and a column compartment. The MS-system was an API 5000 triple/quadruple mass spectrometer with a Turbo V ionization source. The lower limit of quantification (LLOQ) was set to 1.5 ng/mL and the lower limit of detection (LLOD) was set to 0.38 ng/mL. This method reported an intra- and inter-day precision of below 10% for all quality control samples.
Pharmacokinetic and statistical analysis
The data were analysed using nonlinear mixed-effects modelling as implemented in NONMEM version VI (Icon Development Solutions, Ellicott City, Maryland, USA) . Piperaquine plasma concentrations were transformed into their natural logarithms to increase the stability of the numeric analysis. Models were fitted to the data using the first-order conditional estimation (FOCE) method with interaction [39–41]. Census, version 1.1 , and Xpose, version 4.04 , library for R was used for model diagnostics. Perl-speaks-NONMEM (PsN), version 3.4.2,  was used to automate the modelling process and for model diagnostics.
Model discrimination was based the on the objective function value (OFV) computed by NONMEM as −2×log likelihood . The OFV is approximately χ2 distributed and a decrease in OFV of 3.84 and 6.64 is considered a significant drop with p<0.05 and p<0.01, respectively, when adding one additional parameter (one degree of freedom between two nested models).
Structural models with one-, two-, three- and four- disposition compartments were fitted to the data. Several alternative absorption models were investigated; first-order absorption, first-order absorption with lag time, zero-order absorption, sequential zero- and first-order absorption, sequential zero- and first-order absorption with lag time, parallel first- and zero-order absorption, parallel first-order absorption and transit compartment absorption with a fixed number of 1–10 transit compartments. The full implementation of the transit compartment absorption model, which allows the number of transit compartments to vary between patients, was also evaluated .
It was assumed that drug elimination took place from the central compartment and the base model was parameterized as elimination clearance, central volume of distribution, inter-compartmental clearance(s) and peripheral distribution volume(s). Bioavailability was added to the model and the population value was fixed to 100%.
The distribution of the individual parameters was assumed to be log-normal and between-subject variability (BSV) was investigated on all parameters as an exponential random effect [Equation 1].
where P i is the individual estimate for a model parameter (e.g. individual drug clearance) in the ith individual. θ p is the population mean of parameter P and ηi,P is individual ith deviation from the population mean. BSV is estimated from a normal distribution with variance ω2 and zero mean. Between dose occasion variability (BOV) was evaluated on absorption parameters [Equation 2].
Where Pij is the individual parameter estimate for the ith patient on the jth dose occasion, κ is the deviation from the population mean after each dose occasion, taken from a normal distribution with variance Π2 and zero mean. An additive residual error model was assumed since data were transformed into their natural logarithms (i.e. essentially equivalent to an exponential error model on an arithmetic scale). Body weight was tried in the model as a simultaneous incorporation of an allometric function on all clearance (power of 0.75) and volume parameters (power of 1), considering the strong biological prior of this covariate relationship [47–49].
Basic goodness-of-fit plots and simulation-based diagnostics were used to evaluate the final model. Visual and numerical predictive checks  were performed using 2000 simulations at each concentration-time point with binning based on protocol times. The 5th, 50th and 95th percentile of observed data were plotted over the simulated 95% confidence interval of the same percentiles to evaluate the models predictive performance. Bootstrap diagnostics were performed using 1000 re-sampled datasets, stratified on pregnancy.
A Monte Carlo Mapped Power analysis (MCMP)  was conducted based on results from a previously published population analysis . The current study design in terms of sampling times and doses were used to create a modelled data set with 408 pregnant and 408 non-pregnant patients. This data set and the published population pharmacokinetic model was used to estimate the minimum number of individuals needed in each group (i.e. pregnant and non-pregnant) to identify the described covariate relationships (i.e. a 45.0% increase in elimination clearance or a 46.8% increase in bioavailability) at a given power (80%) and significance level (p=0.05). Full power curves were produced by plotting number of patients needed against the power to detect the assumed covariate relationship.
Pregnancy was investigated in the final pharmacokinetic model utilizing a full covariate approach where pregnancy was included as a categorical covariate on all pharmacokinetic parameters. Estimated gestational age (EGA) was evaluated as a continuous covariate on individual parameters using a linear and a power function, and the most appropriate covariate relationship (lowest OFV) was incorporated into a full covariate model for EGA. These two full covariate models were bootstrapped (n=200) to investigate the impact of pregnancy. A pregnancy related change in the parameter estimate of more than 20% was deemed to have clinical relevance.
Fourteen non-pregnant and twelve pregnant women were recruited into the study but two non-pregnant women withdrew their consent and were excluded from the analysis. Full demographics are given in Table 1. The treatment was well-tolerated, none of the patients vomited after treatment and no severe adverse events were reported during the study. One non-pregnant woman had a PCR-confirmed new infection on day 35. None of the women presented recrudescent malaria during the nine weeks of follow-up.
Five hundred sixty-four (564) post-dose plasma samples of piperaquine were used in the pharmacokinetic analysis. A total of 7 (1.2%) of these samples were measured to be below the limit of quantification and omitted. A three-compartmental disposition model resulted in a significantly better model fit compared with a two-compartment model (ΔOFV=−32.8). An additional peripheral compartment (four-compartment disposition model) did not further improve the fit (ΔOFV=0). A transit-compartment (n=3) absorption model was superior to all other absorption models (ΔOFV≥−44.5). In the final model the absorption rate constant (ka) and the transit-compartment rate constant (ktr) were set equal to increase the stability of the model. An additive residual error model was adequate in describing the residual random variability. BSV could be estimated reliably on elimination clearance, one inter-compartmental clearance parameter, and one peripheral volume parameter. The addition of BSV on the bioavailability resulted in a significant improvement in model fit (ΔOFV=−70.9). BOV had considerable impact on both mean transit absorption time and bioavailability (ΔOFV=−340). The final parameter estimates and a schematic picture of the final structural model is presented in Table 2 and Figure 1. Body weight was incorporated with an allometric function on all clearance and volume parameters (ΔOFV=−6.14).
Secondary parameters (i.e. total drug exposure, maximum concentration after first dose, time to maximum concentration, elimination half-life and day 7 concentrations) were analysed using the Mann–Whitney U-test to identify differences between pregnant and non-pregnant women (Table 3). There was a significant difference in terminal elimination half-life (p=0.0014), time to maximum concentration (p=0.0177) and maximum concentration (p=0.0205) [median (range) in pregnant vs. non-pregnant women: 22.1 (19.1-25.8) vs. 25.7 (20.9-33.3) days, 3.07 (1.65-4.64) vs. 1.48 (0.887-4.18) hours and 185 (109–363) vs. 102 (40.6-235) ng/mL, respectively]. However, no significant differences were found in day 7 concentrations (p=0.67), day 28 concentrations (p=0.84) or the total drug exposure (p=0.80) between the pregnant and non-pregnant women.
The MCMP analysis (5% significance level and 80% power) indicated that 8+8 and 13+13 women are needed to detect pregnancy as a covariate on elimination clearance and bioavailability, respectively (Figure 2). A formal stepwise covariate search was therefore not performed since it might result in a biased covariate selection in this small population sample.
Two full covariate models were constructed from the final model to investigate the clinical relevance of pregnancy and EGA separately. Pregnancy had a relatively large impact on mean transit absorption-time, volume of distributions and inter-compartment clearances but no significant effect on elimination clearance (Figure 3). The inclusion of EGA as a power function or a linear function produced similar results. EGA was therefore implemented as a linear function for the full covariate approach and resulted in similar results compared to pregnancy as a categorical covariate (Figure 4A). Bootstrap results for elimination clearance, stratified by trimester (i.e. non-pregnant, second trimester at 20 weeks, and third trimester at 32 weeks) were also investigated and resulted in a non-significant trend of increasing clearance with increasing EGA (Figure 4B).
The final model resulted in good model diagnostic performance and reliable parameter estimates (Figure 5 and Table 2). Calculated epsilon-shrinkage was low (13.0%) which indicates that model diagnostics can be assessed reliably. However, eta-shrinkage was relatively high for certain parameters (CL/F=18.9%, MTT=13.2-52.8%, F=12.4-46.6%) (Table 2) and empirical Bayes estimates should therefore be interpreted with caution (Table 3) . The final model had good predictive performance (Figure 6) with 4.8% (95% CI. 1.4%-11%) of observed data below and 2.1% (95% CI. 1.4%-10%) of observed data above the simulated 90% prediction interval (the three observations at day 90 were excluded because too few patients were followed up to this time for reliable simulations).
In this study, the pharmacokinetic properties of piperaquine have been investigated using nonlinear mixed-effects modelling in pregnant and non-pregnant Sudanese women treated with piperaquine-dihydroartemisinin for uncomplicated P. falciparum malaria. Few studies have been performed to evaluate the effect of pregnancy on the population pharmacokinetics of piperaquine, and this is the first study conducted in an African population. The treatment was well-tolerated and none of the participating women had recrudescent malaria infections.
Previous studies of piperaquine pharmacokinetics have presented both two- and three-compartment disposition models depending on the amount of data included in the analysis [14, 15, 17, 19, 20, 53]. A three-compartment disposition model described the piperaquine concentration-time data adequately in this study. This supports the general finding that a three-compartment disposition is more appropriate than a two-compartment disposition when modelling data from patients followed for a sufficient period of time.
The absorption phase was best described with a transit compartment model with three transit compartments including random effects on bioavailability (BSV and BOV) and mean transit absorption time (BOV). The transit-compartment model provides a more physiological representation of the absorption process compared to the absorption models used in previous studies (i.e. first-order absorption and parallel first-order absorption with lag time) [15, 17]. Recently published studies have also implemented the transit-compartment model but with two- and five-transit compartments which support the absorption model presented here [19, 20]. Small variations in the number of transit compartments are to be expected when modelling different studies due to population differences, sampling schedules and study size. The inclusion of BSV and BOV in the absorption model improved the description of the data in the absorption phase and accommodated the large between-subject and between-occasion variability in these data. The data in the absorption phase was not rich enough to estimate separate absorption rates for ka and ktr, and they were therefore set to be equal. This is a common limitation and the same approach has been used in previous studies [19, 20]. Incorporating BOV resulted in an increasing median bioavailability (0.77, 1.19 and 1.40 at dose 1, 2 and 3, respectively) and mean transit absorption time (1.55, 1.95 and 2.05 at dose 1, 2 and 3, respectively) during the treatment regimen. Similar patterns have been identified in previously studies on piperaquine . This might be an effect of disease recovery or differences in the food intake over the course of the dose regimen. This cannot be verified since parasite densities were not counted at each dose and food intake was not monitored in this study.
Pregnant women had a shorter terminal half-life compared to non-pregnant women, which is in agreement with the non-compartmental analysis , and higher maximum concentrations after the first dose. However, there were no differences in total piperaquine exposure, day 7 concentrations or day 28 concentrations, which supports previously published findings in an Asian pregnant population .
The main aim of this study was to investigate the pharmacokinetic differences between pregnant and non-pregnant women, but the sample size was not large enough to make a conventional covariate search. The MCMP analysis resulted in a minimum of 13 patients needed in each group in order to identify the previously described covariate relationships with 80% power and a significance level of 0.05 (Figure 2). However, this is under the assumption of perfect sampling since the MCMP analysis was based on simulated data using protocol sampling times. In the present study, some patients were not sampled for the complete follow-period and some samples were randomly missing which might increase the number of patients needed to identify the assumed covariate-relationships. Therefore the final model did not include any covariates except body weight, which has a strong biological prior [47, 48] and in addition gave a drop in OFV when included in the model. The full covariate approach suffered from identifiability issues when incorporating the pregnancy covariate simultaneously on clearance parameters, volume parameters and relative bioavailability. The full covariate approach was therefore used to investigate the net effect of a potential covariate on all parameters except relative bioavailability. This approach resulted in a model with reduced volume of distribution and inter-compartment clearance in pregnant women compared with non-pregnant women but no net-effect on apparent clearance. This is in agreement with previously published results where pregnancy affected both clearance and relative bioavailability but in different directions . These covariate relationships would also explain the difference in terminal elimination half-life and the lack of difference in total drug exposure.
The model presented in this study was built on data from few patients and a single individual can therefore have a considerable impact on the results. Piperaquine population pharmacokinetic parameter estimates from the final model are in agreement with previous reports (Table 4). However, the elimination clearance presented in this study for an African population is lower compared to previous studies in non-pregnant and pregnant patients. This might suggest an ethnicity related effect on elimination clearance but this needs to be confirmed in a larger population.
In conclusion, this study presents the population pharmacokinetic properties of piperaquine in pregnant and non-pregnant women with uncomplicated P. falciparum malaria in Sudan. The terminal half-life was shorter in pregnant compared to non-pregnant women, but the total drug exposure was comparable between the two groups. This supports previous findings that no dose adjustments are needed on account of altered piperaquine pharmacokinetics in pregnancy.
WHOS: World malaria report: 2011. 2011, Switzerland: World Health Organization
Lindsay S, Ansell J, Selman C, Cox V, Hamilton K, Walraven G: Effect of pregnancy on exposure to malaria mosquitoes. Lancet. 1972, 2000: 355-
Schantz-Dunn J, Nour NM: Malaria and pregnancy: a global health perspective. Rev Obstet Gynecol. 2009, 2: 186-192.
Abrams ET, Meshnick SR: Malaria during pregnancy in endemic areas: a lens for examining maternal-fetal conflict. Am J Hum Biol. 2009, 21: 643-650. 10.1002/ajhb.20919.
Smereck J: Malaria in pregnancy: update on emergency management. J Emerg Med. 2010, 40: 393-396.
Dafallah SE, El-Agib FH, Bushra GO: Maternal mortality in a teaching hospital in Sudan. Saudi Med J. 2003, 24: 369-372.
Nosten F, Brasseur P: Combination therapy for malaria: the way forward?. Drugs. 2002, 62: 1315-1329. 10.2165/00003495-200262090-00003.
White NJ: Delaying antimalarial drug resistance with combination chemotherapy. Parasitologica. 1999, 41: 301-308.
Ashley EA, Krudsood S, Phaiphun L, Srivilairit S, McGready R, Leowattana W, Hutagalung R, Wilairatana P, Brockman A, Looareesuwan S, Nosten F, White NJ: Randomized, controlled dose-optimization studies of dihydroartemisinin-piperaquine for the treatment of uncomplicated multidrug-resistant falciparum malaria in Thailand. J Infect Dis. 2004, 190: 1773-1782. 10.1086/425015.
Ashley EA, McGready R, Hutagalung R, Phaiphun L, Slight T, Proux S, Thwai KL, Barends M, Looareesuwan S, White NJ, Nosten F: A randomized, controlled study of a simple, once-daily regimen of dihydroartemisinin-piperaquine for the treatment of uncomplicated, multidrug-resistant falciparum malaria. Clin Infect Dis. 2005, 41: 425-432. 10.1086/432011.
Denis MB, Davis TM, Hewitt S, Incardona S, Nimol K, Fandeur T, Poravuth Y, Lim C, Socheat D: Efficacy and safety of dihydroartemisinin-piperaquine (Artekin) in Cambodian children and adults with uncomplicated falciparum malaria. Clin Infect Dis. 2002, 35: 1469-1476. 10.1086/344647.
Karunajeewa H, Lim C, Hung TY, Ilett KF, Denis MB, Socheat D, Davis TM: Safety evaluation of fixed combination piperaquine plus dihydroartemisinin (Artekin) in Cambodian children and adults with malaria. Br J Clin Pharmacol. 2004, 57: 93-99.
Tran TH, Dolecek C, Pham PM, Nguyen TD, Nguyen TT, Le HT, Dong TH, Tran TT, Stepniewska K, White NJ, Farrar J: Dihydroartemisinin-piperaquine against multidrug-resistant Plasmodium falciparum malaria in Vietnam: randomised clinical trial. Lancet. 2004, 363: 18-22. 10.1016/S0140-6736(03)15163-X.
Hung TY, Davis TM, Ilett KF, Karunajeewa H, Hewitt S, Denis MB, Lim C, Socheat D: Population pharmacokinetics of piperaquine in adults and children with uncomplicated falciparum or vivax malaria. Br J Clin Pharmacol. 2004, 57: 253-262.
Roshammar D, Hai TN, Friberg Hietala S, Van Huong N, Ashton M: Pharmacokinetics of piperaquine after repeated oral administration of the antimalarial combination CV8 in 12 healthy male subjects. Eur J Clin Pharmacol. 2006, 62: 335-341. 10.1007/s00228-005-0084-9.
Sim IK, Davis TM, 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.
Tarning J, Ashley EA, Lindegardh N, Stepniewska K, Phaiphun L, Day NP, 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.
Annerberg A, Lwin KM, Lindegardh N, Khrutsawadchai S, Ashley E, Day NP, Singhasivanon P, Tarning J, White NJ, Nosten F: A small amount of fat does not affect piperaquine exposure in patients with malaria. Antimicrob Agents Chemother. 2011, 55: 3971-3976. 10.1128/AAC.00279-11.
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.
Tarning J, Zongo I, Somé 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.
Dawes M, Chowienczyk PJ: Drugs in pregnancy. Pharmacokinetics in pregnancy. Best Pract Res Clin Obstet Gynaecol. 2001, 15: 819-826. 10.1053/beog.2001.0231.
Dean M, Stock B, Patterson RJ, Levy G: Serum protein binding of drugs during and after pregnancy in humans. Clin Pharmacol Ther. 1980, 28: 253-261. 10.1038/clpt.1980.158.
Anderson GD: Pregnancy-induced changes in pharmacokinetics: a mechanistic-based approach. Clin Pharmacokinet. 2005, 44: 989-1008. 10.2165/00003088-200544100-00001.
McGready R, Tan SO, Ashley EA, Pimanpanarak M, Viladpai-Nguen J, Phaiphun L, Wustefeld K, Barends M, Laochan N, Keereecharoen L, Lindegardh N, Singhasivanon P, White NJ, Nosten F: A randomised controlled trial of artemether-lumefantrine versus artesunate for uncomplicated plasmodium falciparum treatment in pregnancy. PLoS Med. 2008, 5: e253-10.1371/journal.pmed.0050253.
Green MD, van Eijk AM, van Ter Kuile FO, Ayisi JG, Parise ME, Kager PA, Nahlen BL, Steketee R, Nettey H: Pharmacokinetics of sulfadoxine-pyrimethamine in HIV-infected and uninfected pregnant women in Western Kenya. J Infect Dis. 2007, 196: 1403-1408. 10.1086/522632.
McGready R, Ashley EA, Moo E, Cho T, Barends M, Hutagalung R, Looareesuwan S, White NJ, Nosten F: A randomized comparison of artesunate-atovaquone-proguanil versus quinine in treatment for uncomplicated falciparum malaria during pregnancy. J Infect Dis. 2005, 192: 846-853. 10.1086/432551.
McGready R, Stepniewska K, Edstein MD, Cho T, Gilveray G, Looareesuwan S, White NJ, Nosten F: The pharmacokinetics of atovaquone and proguanil in pregnant women with acute falciparum malaria. Eur J Clin Pharmacol. 2003, 59: 545-552. 10.1007/s00228-003-0652-9.
McGready R, Stepniewska K, Lindegardh N, Ashley EA, La Y, Singhasivanon P, White NJ, Nosten F: The pharmacokinetics of artemether and lumefantrine in pregnant women with uncomplicated falciparum malaria. Eur J Clin Pharmacol. 2006, 62: 1021-1031. 10.1007/s00228-006-0199-7.
McGready R, Stepniewska K, Ward SA, Cho T, Gilveray G, Looareesuwan S, White NJ, Nosten F: Pharmacokinetics of dihydroartemisinin following oral artesunate treatment of pregnant women with acute uncomplicated falciparum malaria. Eur J Clin Pharmacol. 2006, 62: 367-371. 10.1007/s00228-006-0118-y.
Tarning J, McGready R, Lindegardh N, Ashley EA, Pimanpanarak M, Kamanikom B, Annerberg A, Day NP, Stepniewska K, Singhasivanon P: 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.
Tarning J, Kloprogge F, Piola P, Dhorda M, Muwanga S, Turyakira E, Nuengchamnong N, Nosten F, Day NP, White NJ: 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.
Abdelrahim II, Adam I, Elghazali G, Gustafsson LL, Elbashir MI, Mirghani RA: Pharmacokinetics of quinine and its metabolites in pregnant Sudanese women with uncomplicated Plasmodium falciparum malaria. J Clin Pharm Ther. 2007, 32: 15-19. 10.1111/j.1365-2710.2007.00788.x.
Rijken MJ, McGready R, Jullien V, Tarning J, Lindegardh N, Phyo AP, Win AK, Hsi P, Cammas M, Singhasivanon P, White NJ, Nosten F: Pharmacokinetics of amodiaquine and desethylamodiaquine in pregnant and postpartum women with Plasmodium vivax malaria. Antimicrob Agents Chemother. 2011, 55: 4338-4342. 10.1128/AAC.00154-11.
Tarning J, Chotsiri P, Jullien V, Rijken MJ, Bergstrand M, Cammas M, McGready R, Singhasivanon P, Day NP, White NJ, Nosten F, Lindegardh N: Population pharmacokinetic and pharmacodynamic modeling of amodiaquine and desethylamodiaquine in women with Plasmodium vivax malaria during and after pregnancy. Antimicrob Agents Chemother. 2012, 56: 5764-5773. 10.1128/AAC.01242-12.
Rijken MJ, McGready R, Phyo AP, Lindegardh N, Tarning J, Laochan N, Than HH, Mu O, Win AK, Singhasivanon P, White N, Nosten F: Pharmacokinetics of dihydroartemisinin and piperaquine in pregnant and nonpregnant women with uncomplicated falciparum malaria. Antimicrob Agents Chemother. 2011, 55: 5500-5506. 10.1128/AAC.05067-11.
Adam I, Tarning J, Lindegardh N, Mahgoub H, McGready R, Nosten F: Pharmacokinetics of piperaquine in pregnant women in Sudan with uncomplicated Plasmodium falciparum malaria. AmJTrop Med Hyg. 2012, 87: 35-40. 10.4269/ajtmh.2012.11-0410.
Lindegardh N, Annerberg A, White NJ, Day NP: Development and validation of a liquid chromatographic-tandem mass spectrometric method for determination of piperaquine in plasma stable isotope labeled internal standard does not always compensate for matrix effects. J Chromatogr B Analyt Technol Biomed Life Sci. 2008, 862: 227-236. 10.1016/j.jchromb.2007.12.011.
Beal SLLBS, Boeckmann AJ: NONMEM Users Guides. 2006, City, Maryland, USA: Ellicott City, Maryland, USA
Wahlby U, Bouw MR, Jonsson EN, Karlsson MO: Assessment of type I error rates for the statistical sub-model in NONMEM. J Pharmacokinet Pharmacodyn. 2002, 29: 251-269. 10.1023/A:1020254823597.
Wahlby U, Jonsson EN, Karlsson MO: Assessment of actual significance levels for covariate effects in NONMEM. J Pharmacokinet Pharmacodyn. 2001, 28: 231-252. 10.1023/A:1011527125570.
Beal AJB SL, Sheiner LB: NONMEM users guides. University of Califonia. 1992, San Fransisco, CA: NONMEM project group
Wilkins JJ: NONMEMory: a run management tool for NONMEM. Comput Methods Programs Biomed. 2005, 78: 259-267. 10.1016/j.cmpb.2005.02.003.
Jonsson EN, Karlsson MO: Xpose–an S-PLUS based population pharmacokinetic/pharmacodynamic model building aid for NONMEM. Comput Methods Programs Biomed. 1999, 58: 51-64.
Lindbom L, Pihlgren P, Jonsson EN: PsN-Toolkit–a collection of computer intensive statistical methods for non-linear mixed effect modeling using NONMEM. Comput Methods Programs Biomed. 2005, 79: 241-257. 10.1016/j.cmpb.2005.04.005.
Beal SL, Sheiner LB: Estimating population kinetics. Crit Rev Biomed Eng. 1982, 8: 195-222.
Savic RM, Jonker DM, Kerbusch T, Karlsson MO: Implementation of a transit compartment model for describing drug absorption in pharmacokinetic studies. J Pharmacokinet Pharmacodyn. 2007, 34: 711-726. 10.1007/s10928-007-9066-0.
Anderson BJ, Holford NH: Mechanism-based concepts of size and maturity in pharmacokinetics. Annu Rev Pharmacol Toxicol. 2008, 48: 303-332. 10.1146/annurev.pharmtox.48.113006.094708.
Anderson BJ, Holford NH: Mechanistic basis of using body size and maturation to predict clearance in humans. Drug Metab Pharmacokinet. 2009, 24: 25-36. 10.2133/dmpk.24.25.
McLeay SC, Morrish GA, Kirkpatrick CM, Green B: The relationship between drug clearance and body size: systematic review and meta-analysis of the literature published from 2000 to 2007. Clin Pharmacokinet. 2012, 51: 319-330. 10.2165/11598930-000000000-00000.
Holford N: A degenerative predictive check [abstract]. 14th annual Meeting Population approach group Europe. 2005, 738: 14-
Vong C, Bergstrand M, Nyberg J, Karlsson MO: Rapid sample size calculations for a defined likelihood ratio test-based power in mixed-effects models. AAPS J. 2012, 14: 176-186. 10.1208/s12248-012-9327-8.
Savic RM, Karlsson MO: Importance of shrinkage in empirical bayes estimates for diagnostics: problems and solutions. AAPS J. 2009, 11: 558-569. 10.1208/s12248-009-9133-0.
Karunajeewa HA, Ilett KF, Mueller I, Siba P, Law I, Page-Sharp M, Lin E, Lammey J, Batty KT, Davis TM: Pharmacokinetics and efficacy of piperaquine and chloroquine in Melanesian children with uncomplicated malaria. Antimicrob Agents Chemother. 2008, 52: 237-243. 10.1128/AAC.00555-07.
Salman S, Page-Sharp M, Batty KT, Kose K, Griffin S, Siba PM, Ilett KF, Mueller I, Davis TM: Pharmacokinetic comparison of two piperaquine-containing artemisinin combination therapies in Papua New Guinean children with uncomplicated malaria. Antimicrob Agents Chemother. 2012, 56: 3288-3297. 10.1128/AAC.06232-11.
We sincerely thank all women for their participation in completing this study. We also thank the diligent staff from New Halfa Teaching Hospital. The study treatment was kindly provided by Bejing Holley-Cotec Pharmaceuticals, Co., Ltd. (Bejing, China). Drug assays were supported by the Malaria in Pregnancy (MIP) consortium, which is funded through a grant from the Bill and Melinda Gates Foundation to the Liverpool School of Tropical Medicine. This investigation was part of the Wellcome Trust-Mahidol University-Oxford Tropical Medicine Research Programme, and the PKPDia collaboration, both supported by the Wellcome Trust of Great Britain.
The Welcome Trust is a UK-based medical research charity and is independent of all drug companies. It has no financial links with the manufacturers of either the diagnostics tests or the drugs used in this study. The authors declare no conflict of interest.
IA, FN, ND, and NW conceived the project. WH, NL and JT quantified the drug concentrations. RH and JT performed the pharmacokinetic analysis and wrote the first draft of the manuscript. All authors revised the manuscript critically for important intellectual content and approved the final manuscript.