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.
Method
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.
Results
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.
Conclusions
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 [1]. Pregnant women are at an increased risk of developing a malaria infection [2] 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 [6]. 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 [24]. 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 [19]. 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 [35]. 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.
Methods
Study design
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 [36]. 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.
Drug regimen
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.
Drug analysis
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 [37]. 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) [38]. 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 [42], and Xpose, version 4.04 [43], library for R was used for model diagnostics. Perl-speaks-NONMEM (PsN), version 3.4.2, [44] 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 [45]. 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 [46].
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].
(1)
where Pi 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].
(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 [50] 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) [51] was conducted based on results from a previously published population analysis [19]. 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.
Results
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.
Table 1
Admission demographic data of study population
Non-pregnant women
Pregnant women
P-value
Median (Range)
Median (Range)
Number of patients
12
12
-
Daily piperaquine (phosphate) dose (mg/kg)
18.1 (15.1-24.2)
17.6 (13.6-21.7)
0.884
Daily piperaquine (base) dose (mg/kg)
10.5 (8.71-13.9)
10.2 (7.83-12.5)
0.884
Age (years)
21.0 (16.0-43.0)
26.0 (18.0-33.0)
0.977
Body weight (kg)
53.0 (44.0-81.0)
59.0 (50.0-72.0)
0.544
Height (cm)
163 (150–174)
166 (150–174)
0.908
Estimated gestational age (weeks)
-
32.0 (15.3–40.1)
-
Parasitemia (parasites/μL)
13200 (936–68700)
12900 (624–118000)
0.488
Days of fever
2.5 (1–6)
3 (1–6)
0.095
Fever (°C)
38.3 (36.7-40.0)
38.2 (36.6-39.9)
0.795
Haemoglobin (g/dL)
8.70 (7.60-11.5)
9.65 (8.00-12.0)
0.099
Urea (mg/dL)
26.0 (24.0-28.0)
25.0 (24.0-28.0)
0.036
Serum glutamic pyruvic transaminase (IU/L)
5.00 (2.00-10.0)
4.50 (2.00-9.00)
0.445
Serum glutamic oxaloacetic transaminase (IU/L)
13.5 (2.00-18.0)
13.0 (3.00-21.0)
0.727
The p-values are calculated using the Mann–Whitney U-test.
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).
Table 2
Final parameter estimates describing the piperaquine population pharmacokinetics in women with uncomplicatedP.falciparummalaria
Population estimates [RSE %]
95% CI
BSV/BOV† [RSE %]
95% CI
CL/F (L/h)
44.6 [9.90]
37.3-53.8
22.5 [31.5]
14.6-29.0
VC/F (L)
1820 [11.5]
1450-2240
-
-
Q1/F (L/h)
47.7 [19.0]
32.4-69.2
-
-
VP1/F (L)
15900 [12.3]
12600-20400
-
-
Q2/F (L/h)
352 [11.1]
283-431
-
-
VP2/F (L)
7520 [17.1]
5520-10500
-
-
MTT (h)
1.70 [8.05]
1.45-2.00
60.7 [22.8] †
44.5-76.7
RUV
0.0973 [5.90]
0.0753-0.120
-
-
No. of trans comp
3 fix
-
-
-
F (%)
100 fix
-
34.7 [59.2]
9.52-54.9
64.8 [15.2] †
52.8-76.0
Coefficient of variation (%CV) for between-subject variability (BSV) and (†) between occasion variability (BOV) were calculated as the [exp(estimated variance)-1]1/2. Relative standard errors (RSE) and the 95% confidence intervals (CI) were based on 868 successful stratified bootstrap runs (out of 1000) and presented as 100×(standard deviation/mean value) and as 2.5 to 97.5 percentiles, respectively.
CL/F is the apparent elimination clearance. Vc/F is the apparent volume of distribution of the central compartment. Q1/F and Q2/F is the inter-compartment clearance between the central and first and second peripheral compartment, respectively. VP1/F and VP2/F is the apparent volume of distribution of first and second peripheral compartment, respectively. MTT is the mean transit time of the absorption model. RUV is the variance of the residual variability. No. of trans comp is the number of transit compartments used in the absorption phase. F represents the relative bioavailability.
Figure 1
The structural model for piperaquine pharmacokinetics. ktr is the 1st order absorption rate. Vc, Vp1 and Vp2 are the apparent volume of the central compartment, the first peripheral compartment and the second peripheral compartment, respectively. CL is the apparent elimination clearance and Q1 and Q2 is the apparent inter-compartmental clearances.
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.
Table 3
Secondary parameters of piperaquine pharmacokinetics in pregnant and non-pregnant women with uncomplicatedP.falciparummalaria
Secondary parameters
Total
Non-pregnant women
Pregnant women
p-value
CMAX (ng/mL)
158 [40.6-363]
102 [40.6-235]
185 [109–363]
0.021
TMAX (hours)
2.62 [0.887-4.64]
1.48 [0.887-4.18]
3.07 [1.65-4.64]
0.018
Half-life (days)
23.4 [19.1-33.3]
25.7 [20.9-33.3]
22.1 [19.1-25.8]
0.001
Day 7 concentration (ng/mL)
58.3 [16.6-146]
55.4 [16.6-146]
60.7 [40.1-103]
0.671
Day 28 concentration (ng/mL)
15.9 [4.85-38.6]
15.4 [4.85-38.6]
16.1 [9.68-26.8]
0.840
AUC0->90 (ng*h/mL)
40600 [12400–100000]
38000 [12400–100000]
42700 [27100–68700]
0.799
AUC48h->90 (ng*h/mL)
36400 [10600–90300]
35300 [10600–90300]
37700 [23500–63200]
0.887
Secondary parameters are predicted using the final model and values are presented as median [range]. The p-values are calculated with a Mann–Whitney U-test.
CMAX is the predicted maximum concentration after the first dose and TMAX is the time to CMAX. Half-life is the estimated terminal elimination half-life. Day 7 and 28 concentrations are the model predicted plasma concentrations of piperaquine at day 7 and 28, respectively. AUC0->90 is the predicted area under the concentration-time curve from time zero extrapolated to day 90. AUC48h->90 is the predicted area under the concentration-time curve from 48 hours extrapolated to day 90.
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.
Figure 2
Monte-Carlo Mapped Power (MCMP) curve for identifying pregnancy as a covariate. Triangles represents the power curve for indentifying pregnancy as covariate on apparent elimination clearance and circles is the power curve for identifying pregnancy as a covariate on the relative bioavailability. The dotted black line represents 80% power. The inserted numbers are the total number of subjects needed to identify pregnancy as a covariate given the used model and study sampling procedure.
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).
Figure 3
Box (25thto 75thpercentile) and whisker (1.5*interquartile range) plot of the full pregnancy-covariate model for piperaquine. Pregnancy was included as a categorical covariate and the solid black zero-line represents no covariate effect and the dotted black lines represent a covariate effect of ±20%. MTT is the mean transit time, Vc/F, Vp1/F and Vp2/F is the apparent volume of the central compartment, the first peripheral compartment and the second peripheral compartment respectively. CL/F is the apparent elimination clearance and Q1/F and Q2/F is the apparent inter-compartmental clearances.
Figure 4
Box (25thto 75thpercentile) and whisker (1.5*interquartile range) plot of the full covariate model with estimated gestational age implemented as a linear continuous covariate. Panel A displays the percent change per week of gestational age where the solid black zero-line represents no covariate effect. MTT is the mean transit time, Vc/F, Vp1/F and Vp2/F is the apparent volume of the central compartment, the first peripheral compartment ant the second peripheral compartment, respectively. CL/F is the apparent elimination clearance and Q1/F and Q2/F is the apparent inter-compartmental clearances. Panel B displays the elimination clearance estimates for non-pregnant women, women in the second trimester (week 20) and women in the third trimester (week 32).
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) [52]. 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).
Figure 5
Basic goodness-of-fit plots for the final piperaquine model. Observations plotted against population predicted concentrations (A) and against individual predicted concentrations (B). Conditional weighted residuals plotted against population predicted concentrations (C) and time (D). The solid line is the identity line and the broken line is the locally weighted least square regression line. The concentrations were transformed into their logarithms (base 10). Conditional weighted residuals were fixed between 4 and −4, excluding one data point.
Figure 6
The visual predictive check of the final piperaquine model.. The circles represent the observations, the solid line represents the 5th, 50th and 95th percentile of the observed data. The shaded areas represent the simulated 95% confidence interval around the 5th, 50th and 95th percentile. Concentrations were transformed to their logarithms (base 10). Predicted performance during the three first days is inserted in the top right corner.
Discussion
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 [19]. 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 [36], 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 [19].
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 [19]. 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.
Table 4
A literature comparison of the pharmacokinetic properties of piperaquine
Study population
Age (years)
Study size (Males/Females)
t1/2 (days)
CL/F (L/h/kg)
VD/F (L/kg)
Method
Reference
Piperaquine pharmacokinetics in pregnant women
Pregnant Sudanese women with uncomplicated P. falciparum malaria
18-33
12 (0/12)
22.1
0.678
384
Pop PK
This study
Pregnant Thai and Karen women with uncomplicated P. falciparum malaria
CL/F is the apparent elimination clearance, t1/2 is the elimination half-life, VD/F is the apparent volume of distribution, Pop PK represents a population pharmacokinetic analysis and NCA represents a non-compartmental analysis. Age is given as a range or (†) as mean ± standard deviation. Other parameters are given as median or as mean when indicated (*).
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.
Declarations
Acknowledgment
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.
Authors’ Affiliations
(1)
Unit for Pharmacokinetics and Drug Metabolism, Department of Pharmacology, University of Gothenburg
(2)
Faculty of Medicine, University of Khartoum
(3)
Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University
(4)
Centre for Tropical Medicine, Department of Clinical Medicine, Churchill Hospital
(5)
Shoklo Malaria Research Unit
References
WHOS: World malaria report: 2011. 2011, Switzerland: World Health OrganizationGoogle Scholar
Lindsay S, Ansell J, Selman C, Cox V, Hamilton K, Walraven G: Effect of pregnancy on exposure to malaria mosquitoes. Lancet. 1972, 2000: 355-Google Scholar
Schantz-Dunn J, Nour NM: Malaria and pregnancy: a global health perspective. Rev Obstet Gynecol. 2009, 2: 186-192.PubMed CentralPubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
Smereck J: Malaria in pregnancy: update on emergency management. J Emerg Med. 2010, 40: 393-396.View ArticlePubMedGoogle Scholar
Dafallah SE, El-Agib FH, Bushra GO: Maternal mortality in a teaching hospital in Sudan. Saudi Med J. 2003, 24: 369-372.PubMedGoogle Scholar
Nosten F, Brasseur P: Combination therapy for malaria: the way forward?. Drugs. 2002, 62: 1315-1329. 10.2165/00003495-200262090-00003.View ArticlePubMedGoogle Scholar
White NJ: Delaying antimalarial drug resistance with combination chemotherapy. Parasitologica. 1999, 41: 301-308.Google Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
Anderson GD: Pregnancy-induced changes in pharmacokinetics: a mechanistic-based approach. Clin Pharmacokinet. 2005, 44: 989-1008. 10.2165/00003088-200544100-00001.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.View ArticleGoogle Scholar
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.View ArticlePubMedGoogle Scholar
Beal SLLBS, Boeckmann AJ: NONMEM Users Guides. 2006, City, Maryland, USA: Ellicott City, Maryland, USAGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
Beal AJB SL, Sheiner LB: NONMEM users guides. University of Califonia. 1992, San Fransisco, CA: NONMEM project groupGoogle Scholar
Wilkins JJ: NONMEMory: a run management tool for NONMEM. Comput Methods Programs Biomed. 2005, 78: 259-267. 10.1016/j.cmpb.2005.02.003.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
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.View ArticlePubMedGoogle Scholar
Holford N: A degenerative predictive check [abstract]. 14th annual Meeting Population approach group Europe. 2005, 738: 14-Google Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
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.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.