Integration of data from multiple sources for simultaneous modelling analysis: experience from nevirapine population pharmacokinetics

Evaluating pediatric tuberculosis dosing guidelines: A model-based individual data pooled analysis

BACKGROUND

The current World Health Organization (WHO) pediatric tuberculosis dosing guidelines lead to suboptimal drug exposures. Identifying factors altering the exposure of these drugs in children is essential for dose optimization. Pediatric pharmacokinetic studies are usually small, leading to high variability and uncertainty in pharmacokinetic results between studies. We pooled data from large pharmacokinetic studies to identify key covariates influencing drug exposure to optimize tuberculosis dosing in children.

METHODS AND FINDINGS

We used nonlinear mixed-effects modeling to characterize the pharmacokinetics of rifampicin, isoniazid, and pyrazinamide, and investigated the association of human immunodeficiency virus (HIV), antiretroviral therapy (ART), drug formulation, age, and body size with their pharmacokinetics. Data from 387 children from South Africa, Zambia, Malawi, and India were available for analysis; 47% were female and 39% living with HIV (95% on ART). Median (range) age was 2.2 (0.2 to 15.0) years and weight 10.9 (3.2 to 59.3) kg. Body size (allometry) was used to scale clearance and volume of distribution of all 3 drugs. Age affected the bioavailability of rifampicin and isoniazid; at birth, children had 48.9% (95% confidence interval (CI) [36.0%, 61.8%]; p < 0.001) and 64.5% (95% CI [52.1%, 78.9%]; p < 0.001) of adult rifampicin and isoniazid bioavailability, respectively, and reached full adult bioavailability after 2 years of age for both drugs. Age also affected the clearance of all drugs (maturation), children reached 50% adult drug clearing capacity at around 3 months after birth and neared full maturation around 3 years of age. While HIV per se did not affect the pharmacokinetics of first-line tuberculosis drugs, rifampicin clearance was 22% lower (95% CI [13%, 28%]; p < 0.001) and pyrazinamide clearance was 49% higher (95% CI [39%, 57%]; p < 0.001) in children on lopinavir/ritonavir; isoniazid bioavailability was reduced by 39% (95% CI [32%, 45%]; p < 0.001) when simultaneously coadministered with lopinavir/ritonavir and was 37% lower (95% CI [22%, 52%]; p < 0.001) in children on efavirenz. Simulations of 2010 WHO-recommended pediatric tuberculosis doses revealed that, compared to adult values, rifampicin exposures are lower in most children, except those younger than 3 months, who experience relatively higher exposure for all drugs, due to immature clearance. Increasing the rifampicin doses in children older than 3 months by 75 mg for children weighing <25 kg and 150 mg for children weighing >25 kg could improve rifampicin exposures. Our analysis was limited by the differences in availability of covariates among the pooled studies.

CONCLUSIONS

Children older than 3 months have lower rifampicin exposures than adults and increasing their dose by 75 or 150 mg could improve therapy. Altered exposures in children with HIV is most likely caused by concomitant ART and not HIV per se. The importance of the drug-drug interactions with lopinavir/ritonavir and efavirenz should be evaluated further and considered in future dosing guidance.

TRIAL REGISTRATION

ClinicalTrials.gov registration numbers; NCT02348177, NCT01637558, ISRCTN63579542.

Towards Model-Informed Precision Dosing of Voriconazole: Challenging Published Voriconazole Nonlinear Mixed-Effects Models with Real-World Clinical Data

BACKGROUND AND OBJECTIVES

Model-informed precision dosing (MIPD) frequently uses nonlinear mixed-effects (NLME) models to predict and optimize therapy outcomes based on patient characteristics and therapeutic drug monitoring data. MIPD is indicated for compounds with narrow therapeutic range and complex pharmacokinetics (PK), such as voriconazole, a broad-spectrum antifungal drug for prevention and treatment of invasive fungal infections. To provide guidance and recommendations for evidence-based application of MIPD for voriconazole, this work aimed to (i) externally evaluate and compare the predictive performance of a published so-called 'hybrid' model for MIPD (an aggregate model comprising features and prior information from six previously published NLME models) versus two 'standard' NLME models of voriconazole, and (ii) investigate strategies and illustrate the clinical impact of Bayesian forecasting for voriconazole.

METHODS

A workflow for external evaluation and application of MIPD for voriconazole was implemented. Published voriconazole NLME models were externally evaluated using a comprehensive in-house clinical database comprising nine voriconazole studies and prediction-/simulation-based diagnostics. The NLME models were applied using different Bayesian forecasting strategies to assess the influence of prior observations on model predictivity.

RESULTS

The overall best predictive performance was obtained using the aggregate model. However, all NLME models showed only modest predictive performance, suggesting that (i) important PK processes were not sufficiently implemented in the structural submodels, (ii) sources of interindividual variability were not entirely captured, and (iii) interoccasion variability was not adequately accounted for. Predictive performance substantially improved by including the most recent voriconazole observations in MIPD.

CONCLUSION

Our results highlight the potential clinical impact of MIPD for voriconazole and indicate the need for a comprehensive (pre-)clinical database as basis for model development and careful external model evaluation for compounds with complex PK before their successful use in MIPD.

Efficient and relevant stepwise covariate model building for pharmacometrics

Covariate modeling is an important opportunity for pharmacometrics to influence decision making in drug development. The stepwise covariate model (SCM) building procedure is the most common method for covariate model development. Despite its advantages, the traditional SCM method is known to have long runtimes and the suboptimal ability to select relevant covariates, especially in more complex phase III settings. In this work, two alternative approaches are presented: SCM+, which introduces the “adaptive scope reduction” and changes to general estimation settings, and “stage‐wise filtering,” which groups covariates into categories based on their importance (mechanistic, structural, and exploratory). The three methods (SCM, SCM+, and SCM+ with stage‐wise filtering) are applied to data from a simulated phase III population pharmacokinetic study and are compared in terms of efficiency and relevance. The two SCM+ methods were considerably more efficient than the traditional SCM: the number of function evaluations was reduced by 70% for SCM+ and by 76% for SCM+ with stage‐wise filtering compared to SCM; the corresponding number of executed models was reduced by 44% for SCM+ and 70% for SCM+ with stage‐wise filtering. In addition, among the three methods, SCM+ with stage‐wise filtering selected the highest number of relevant covariates. Given the improved efficiency and ability to select relevant covariates shown in this work, the use of SCM+ and stage‐wise filtering can greatly increase the efficiency of covariate modeling in drug development, which will ultimately facilitate more timely support for decision making.

Pharmacokinetics of standard versus high-dose isoniazid for treatment of multidrug-resistant tuberculosis

BACKGROUND

The WHO-endorsed shorter-course regimen for MDR-TB includes high-dose isoniazid. The pharmacokinetics of high-dose isoniazid within MDR-TB regimens has not been well described.

OBJECTIVES

To characterize isoniazid pharmacokinetics at 5-15 mg/kg as monotherapy or as part of the MDR-TB treatment regimen.

METHODS

We used non-linear mixed-effects modelling to evaluate the combined data from INHindsight, a 7 day early bactericidal activity study with isoniazid monotherapy, and PODRtb, an observational study of patients on MDR-TB treatment including terizidone, pyrazinamide, moxifloxacin, kanamycin, ethionamide and/or isoniazid.

RESULTS

A total of 58 and 103 participants from the INHindsight and PODRtb studies, respectively, were included in the analysis. A two-compartment model with hepatic elimination best described the data. N-acetyltransferase 2 (NAT2) genotype caused multi-modal clearance, and saturable first-pass was observed beyond 10 mg/kg dosing. Saturable isoniazid kinetics predicted an increased exposure of approximately 50% beyond linearity at 20 mg/kg dosing. Participants treated with the MDR-TB regimen had a 65.6% lower AUC compared with participants on monotherapy. Ethionamide co-administration was associated with a 29% increase in isoniazid AUC.

CONCLUSIONS

Markedly lower isoniazid exposures were observed in participants on combination MDR-TB treatment compared with monotherapy. Isoniazid displays saturable kinetics at doses >10 mg/kg. The safety implications of these phenomena remain unclear.

Population pharmacokinetics of ethambutol in African children: a pooled analysis

OBJECTIVES

Ethambutol protects against the development of resistance to co-administered drugs in the intensive phase of first-line anti-TB treatment in children. It is especially relevant in settings with a high prevalence of HIV or isoniazid resistance. We describe the population pharmacokinetics of ethambutol in children with TB to guide dosing in this population.

METHODS

We pooled data from 188 intensively sampled children from the DATiC, DNDi and SHINE studies, who received 15-25 mg/kg ethambutol daily according to WHO guidelines. The median (range) age and weight of the cohort were 1.9 (0.3-12.6) years and 9.6 (3.9-34.5) kg, respectively. Children with HIV (HIV+; n = 103) received ART (lopinavir/ritonavir in 92%).

RESULTS

Ethambutol pharmacokinetics were best described by a two-compartment model with first-order elimination and absorption transit compartments. Clearance was estimated to reach 50% of its mature value by 2 months after birth and 99% by 3 years. Typical steady-state apparent clearance in a 10 kg child was 15.9 L/h. In HIV+ children on lopinavir/ritonavir, bioavailability was reduced by 32% [median (IQR) steady-state Cmax = 0.882 (0.669-1.28) versus 1.66 (1.21-2.15) mg/L). In young children, bioavailability correlated with age. At birth, bioavailability was 73.1% of that in children 3.16 years or older.

CONCLUSIONS

To obtain exposure within the 2-6 mg/L recommended range for Cmax, the current doses must be doubled (or tripled with HIV+ children on lopinavir/ritonavir) for paediatric patients. This raises concerns regarding the potential for ocular toxicity, which would require evaluation.

Evaluation of the estimation and classification performance of NONMEM when applying mixture model for drug clearance

Maximum likelihood estimation of parameters involving mixture model is known to have significant and specific patterns of errors. Population pharmacokinetic (PopPK) modeling using NONMEM is no exception. A few relevant studies on estimation and classification performance were done, but a comprehensive study was not yet available. The current study aims to evaluate performance and likelihood ratio test (LRT)‐based true covariate detection rate when fitting a bimodal mixture of drug clearance (CL) in NONMEM. A large number of PopPK datasets with various settings were simulated and then estimated. The estimates were compared to the simulated values and summarized. The separation between the CL distributions of the two subpopulations is systematically overestimated. The major factor associated with the performance is the change in the minimum objective function value after removing the mixture component (dOFV). Other significant factors include estimated disparity index (DI), estimated mixing proportion, and number of subjects in the dataset. Small dOFV and large estimated DI are associated with the worst performance. Omitting a true mixture resulted in reduced true covariate detection rates. It is recommended that on top of routinely generated standard errors and model diagnostics, dOFV, and other factors when necessary, should be taken into account for the evaluation of performance when fitting mixture model using NONMEM. In addition, when fitting mixture model for CL is intended, the mixture component should be introduced prior to LRT‐based covariate model development for CL.

Abacavir pharmacokinetics in African children living with HIV: A pooled analysis describing the effects of age, malnutrition and common concomitant medications

Aims

Abacavir is part of WHO‐recommended regimens to treat HIV in children under 15 years of age. In a pooled analysis across four studies, we describe abacavir population pharmacokinetics to investigate the influence of age, concomitant medications, malnutrition and formulation.

Methods

A total of 230 HIV‐infected African children were included, with median (range) age of 2.1 (0.1–12.8) years and weight of 9.8 (2.5–30.0) kg. The population pharmacokinetics of abacavir was described using nonlinear mixed‐effects modelling.

Results

Abacavir pharmacokinetics was best described by a two‐compartment model with first‐order elimination, and absorption described by transit compartments. Clearance was predicted around 54% of its mature value at birth and 90% at 10 months. The estimated typical clearance at steady state was 10.7 L/h in a child weighing 9.8 kg co‐treated with lopinavir/ritonavir, and was 12% higher in children receiving efavirenz. During coadministration of rifampicin‐based antituberculosis treatment and super‐boosted lopinavir in a 1:1 ratio with ritonavir, abacavir exposure decreased by 29.4%. Malnourished children living with HIV had higher abacavir exposure initially, but this effect waned with nutritional rehabilitation. An additional 18.4% reduction in clearance after the first abacavir dose was described, suggesting induction of clearance with time on lopinavir/ritonavir‐based therapy. Finally, absorption of the fixed dose combination tablet was 24% slower than the abacavir liquid formulation.

Conclusion

In this pooled analysis we found that children on lopinavir/ritonavir or efavirenz had similar abacavir exposures, while concomitant TB treatment and super‐boosted lopinavir gave significantly reduced abacavir concentrations.

A Model-Informed Method for the Purpose of Precision Dosing of Isoniazid in Pulmonary Tuberculosis

Background and Objective

This study aimed to develop and evaluate a population pharmacokinetic model and limited sampling strategy for isoniazid to be used in model-based therapeutic drug monitoring.

Methods

A population pharmacokinetic model was developed based on isoniazid and acetyl-isoniazid pharmacokinetic data from seven studies with in total 466 patients from three continents. Three limited sampling strategies were tested based on the available sampling times in the dataset and practical considerations. The tested limited sampling strategies sampled at 2, 4, and 6 h, 2 and 4 h, and 2 h after dosing. The model-predicted area under the concentration–time curve from 0 to 24 h (AUC₂₄) and the peak concentration from the limited sampling strategies were compared to predictions using the full pharmacokinetic curve. Bias and precision were assessed using the mean error (ME) and the root mean square error (RMSE), both expressed as a percentage of the mean model-predicted AUC₂₄ or peak concentration on the full pharmacokinetic curve.

Results

Performance of the developed model was acceptable and the uncertainty in parameter estimations was generally low (the highest relative standard error was 39% coefficient of variation). The limited sampling strategy with sampling at 2 and 4 h was determined as most suitable with an ME of 1.1% and RMSE of 23.4% for AUC₂₄ prediction, and ME of 2.7% and RMSE of 23.8% for peak concentration prediction. For the performance of this strategy, it is important that data on both isoniazid and acetyl-isoniazid are used. If only data on isoniazid are available, a limited sampling strategy using 2, 4, and 6 h can be employed with an ME of 1.7% and RMSE of 20.9% for AUC₂₄ prediction, and ME of 1.2% and RMSE of 23.8% for peak concentration prediction.

Conclusions

A model-based therapeutic drug monitoring strategy for personalized dosing of isoniazid using sampling at 2 and 4 h after dosing was successfully developed. Prospective evaluation of this strategy will show how it performs in a clinical therapeutic drug monitoring setting.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-020-00971-2.

A population pharmacokinetics analysis assessing the exposure of raltegravir once-daily 1200 mg in pregnant women living with HIV

Once-daily two 600 mg tablets (1200 mg q.d.) raltegravir offers an easier treatment option compared to the twice-daily regimen of one 400 mg tablet. No pharmacokinetic, efficacy, or safety data of the 1200 mg q.d. regimen have been reported in pregnant women to date as it is challenging to collect these clinical data. This study aimed to develop a population pharmacokinetic (PopPK) model to predict the pharmacokinetic profile of raltegravir 1200 mg q.d. in pregnant women and to discuss the expected pharmacodynamic properties of raltegravir 1200 mg q.d. during pregnancy based on previously reported concentration-effect relationships. Data from 11 pharmacokinetic studies were pooled (n = 221). A two-compartment model with first-order elimination and absorption through three sequential transit compartments best described the data. We assessed that the bio-availability of the 600 mg tablets was 21% higher as the 400 mg tablets, and the bio-availability in pregnant women was 49% lower. Monte-Carlo simulations were performed to predict the pharmacokinetic profile of 1200 mg q.d. in pregnant and nonpregnant women. The primary criteria for efficacy were that the lower bound of the 90% confidence interval (CI) of the concentration before next dose administration (C_trough ) geometric mean ratio (GMR) of simulated pregnant/nonpregnant women had to be greater than 0.75. The simulated raltegravir C_trough GMR (90% CI) was 0.51 (0.41-0.63), hence not meeting the primary target for efficacy. Clinical data from two pregnant women using 1200 mg q.d. raltegravir showed a similar C_trough ratio pregnant/nonpregnant. Our pharmacokinetic results support the current recommendation of not using the raltegravir 1200 mg q.d. regimen during pregnancy until more data on the exposure-response relationship becomes available.

Population pharmacokinetic meta-analysis of individual data to design the first randomized efficacy trial of vancomycin in neonates and young infants

OBJECTIVES

In the absence of consensus, the present meta-analysis was performed to determine an optimal dosing regimen of vancomycin for neonates.

METHODS

A 'meta-model' with 4894 concentrations from 1631 neonates was built using NONMEM, and Monte Carlo simulations were performed to design an optimal intermittent infusion, aiming to reach a target AUC0-24 of 400 mg·h/L at steady-state in at least 80% of neonates.

RESULTS

A two-compartment model best fitted the data. Current weight, postmenstrual age (PMA) and serum creatinine were the significant covariates for CL. After model validation, simulations showed that a loading dose (25 mg/kg) and a maintenance dose (15 mg/kg q12h if <35 weeks PMA and 15 mg/kg q8h if ≥35 weeks PMA) achieved the AUC0-24 target earlier than a standard 'Blue Book' dosage regimen in >89% of the treated patients.

CONCLUSIONS

The results of a population meta-analysis of vancomycin data have been used to develop a new dosing regimen for neonatal use and to assist in the design of the model-based, multinational European trial, NeoVanc.

Development and Application of a Mechanistic Population Modeling Approach to Describe Abemaciclib Pharmacokinetics

Abemaciclib is an oral anticancer drug that inhibits cyclin dependent kinases 4 and 6 and is metabolized by cytochrome P450 3A in the intestines and liver to active metabolites. The objectives were (1) to develop a mechanistic model to characterize the pharmacokinetics (PK) of the active moieties and investigate the effect of patient factors and (2) apply the model to dat from two phase III breast cancer trials of abemaciclib in combination with endocrine therapy. To develop the model, data from seven phase I studies and two phase II studies including 421 patients with cancer and 65 healthy individuals were pooled for nonlinear mixed effects modeling. The PK was similar between patients and healthy subjects, and the effects of diarrhea, formulation, race, and patient covariates on exposure were negligible. Application of the model confirmed its predictive performance and that abemaciclib PK did not change when coadministered with endocrine therapy.

Evaluating darunavir/ritonavir dosing regimens for HIV-positive pregnant women using semi-mechanistic pharmacokinetic modelling

BACKGROUND

Darunavir 800 mg once (q24h) or 600 mg twice (q12h) daily combined with low-dose ritonavir is used to treat HIV-positive pregnant women. Decreased total darunavir exposure (17%-50%) has been reported during pregnancy, but limited data on unbound exposure are available.

OBJECTIVES

To evaluate total and unbound darunavir exposures following standard darunavir/ritonavir dosing and to explore the value of potential optimized darunavir/ritonavir dosing regimens for HIV-positive pregnant women.

PATIENTS AND METHODS

A population pharmacokinetic analysis was conducted based on data from 85 women. The final model was used to simulate total and unbound darunavir AUC0-τ and Ctrough during the third trimester of pregnancy, as well as to assess the probability of therapeutic exposure.

RESULTS

Simulations predicted that total darunavir exposure (AUC0-τ) was 24% and 23% lower in pregnancy for standard q24h and q12h dosing, respectively. Unbound darunavir AUC0-τ was 5% and 8% lower compared with post-partum for standard q24h and q12h dosing, respectively. The probability of therapeutic exposure (unbound) during pregnancy was higher for standard q12h dosing (99%) than for q24h dosing (94%).

CONCLUSIONS

The standard q12h regimen resulted in maximal and higher rates of therapeutic exposure compared with standard q24h dosing. Darunavir/ritonavir 600/100 mg q12h should therefore be the preferred regimen during pregnancy unless (adherence) issues dictate q24h dosing. The value of alternative dosing regimens seems limited.

An Individual Participant Data Population Pharmacokinetic Meta-analysis of Drug-Drug Interactions between Lumefantrine and Commonly Used Antiretroviral Treatment

Treating malaria in HIV-coinfected individuals should consider potential drug-drug interactions. Artemether-lumefantrine is the most widely recommended treatment for uncomplicated malaria globally. Lumefantrine is metabolized by CYP3A4, an enzyme that commonly used antiretrovirals often induce or inhibit. A population pharmacokinetic meta-analysis was conducted using individual participant data from 10 studies with 6,100 lumefantrine concentrations from 793 nonpregnant adult participants (41% HIV-malaria-coinfected, 36% malaria-infected, 20% HIV-infected, and 3% healthy volunteers).

Treating malaria in HIV-coinfected individuals should consider potential drug-drug interactions. Artemether-lumefantrine is the most widely recommended treatment for uncomplicated malaria globally. Lumefantrine is metabolized by CYP3A4, an enzyme that commonly used antiretrovirals often induce or inhibit. A population pharmacokinetic meta-analysis was conducted using individual participant data from 10 studies with 6,100 lumefantrine concentrations from 793 nonpregnant adult participants (41% HIV-malaria-coinfected, 36% malaria-infected, 20% HIV-infected, and 3% healthy volunteers). Lumefantrine exposure increased 3.4-fold with coadministration of lopinavir-ritonavir-based antiretroviral therapy (ART), while it decreased by 47% with efavirenz-based ART and by 59% in the patients with rifampin-based antituberculosis treatment. Nevirapine- or dolutegravir-based ART and malaria or HIV infection were not associated with significant effects. Monte Carlo simulations showed that those on concomitant efavirenz or rifampin have 49% and 80% probability of day 7 concentrations <200 ng/ml, respectively, a threshold associated with an increased risk of treatment failure. The risk of achieving subtherapeutic concentrations increases with larger body weight. An extended 5-day and 6-day artemether-lumefantrine regimen is predicted to overcome these drug-drug interactions with efavirenz and rifampin, respectively.

A Mechanism-Based Population Pharmacokinetic Analysis Assessing the Feasibility of Efavirenz Dose Reduction to 400 mg in Pregnant Women

BACKGROUND

Reducing the dose of efavirenz can improve safety, reduce costs, and increase access for patients with HIV infection. According to the World Health Organization, a similar dosing strategy for all patient populations is desirable for universal roll-out; however, it remains unknown whether the 400 mg daily dose is adequate during pregnancy.

METHODS

We developed a mechanistic population pharmacokinetic model using pooled data from women included in seven studies (1968 samples, 774 collected during pregnancy). Total and free efavirenz exposure (AUC24 and C12) were predicted for 400 (reduced) and 600 mg (standard) doses in both pregnant and non-pregnant women.

RESULTS

Using a 400 mg dose, the median efavirenz total AUC24 and C12 during the third trimester of pregnancy were 91 and 87% of values among non-pregnant women, respectively. Furthermore, the median free efavirenz C12 and AUC24 were predicted to increase during pregnancy by 11 and 15%, respectively.

CONCLUSIONS

It was predicted that reduced-dose efavirenz provides adequate exposure during pregnancy. These findings warrant prospective confirmation.

Population Pharmacokinetics of the Antimalarial Amodiaquine: a Pooled Analysis To Optimize Dosing

Amodiaquine plus artesunate is the recommended antimalarial treatment in many countries where malaria is endemic. However, pediatric doses are largely based on a linear extrapolation from adult doses.

Amodiaquine plus artesunate is the recommended antimalarial treatment in many countries where malaria is endemic. However, pediatric doses are largely based on a linear extrapolation from adult doses. We pooled data from previously published studies on the pharmacokinetics of amodiaquine, to optimize the dose across all age groups. Adults and children with uncomplicated malaria received daily weight-based doses of amodiaquine or artesunate-amodiaquine over 3 days. Plasma concentration-time profiles for both the parent drug and the metabolite were characterized using nonlinear mixed-effects modeling. Amodiaquine pharmacokinetics were adequately described by a two-compartment disposition model, with first-order elimination leading to the formation of desethylamodiaquine, which was best described by a three-compartment disposition model. Body size and age were the main covariates affecting amodiaquine clearance. After adjusting for the effect of weight, clearance rates for amodiaquine and desethylamodiaquine reached 50% of adult maturation at 2.8 months (95% confidence interval [CI], 1.5 to 3.7 months) and 3.9 months (95% CI, 2.6 to 5.3 months) after birth, assuming that the baby was born at term. Bioavailability was 22.4% (95% CI, 15.6 to 31.9%) lower at the start of treatment than during convalescence, which suggests a malaria disease effect. Neither the drug formulation nor the hemoglobin concentration had an effect on any pharmacokinetic parameters. Results from simulations showed that current manufacturer dosing recommendations resulted in low desethylamodiaquine exposure in patients weighing 8 kg, 15 to 17 kg, 33 to 35 kg, and >62 kg compared to that in a typical 50-kg patient. We propose possible optimized dosing regimens to achieve similar drug exposures among all age groups, which require further validation.

Population Pharmacokinetic Properties of Sulfadoxine and Pyrimethamine: a Pooled Analysis To Inform Optimal Dosing in African Children with Uncomplicated Malaria

Sulfadoxine-pyrimethamine with amodiaquine is recommended by the World Health Organization as seasonal malaria chemoprevention for children aged 3 to 59 months in the sub-Sahel regions of Africa. Suboptimal dosing in children may lead to treatment failure and increased resistance. Pooled individual patient data from four previously published trials on the pharmacokinetics of sulfadoxine and pyrimethamine in 415 pediatric and 386 adult patients were analyzed using nonlinear mixed-effects modeling to evaluate the current dosing regimen and, if needed, to propose an optimized dosing regimen for children under 5 years of age. The population pharmacokinetics of sulfadoxine and pyrimethamine were both best described by a one-compartment disposition model with first-order absorption and elimination. Body weight, age, and nutritional status (measured as the weight-for-age Z-score) were found to be significant covariates. Allometric scaling with total body weight and the maturation of clearance in children by postgestational age improved the model fit. Underweight-for-age children were found to have 15.3% and 26.7% lower bioavailabilities of sulfadoxine and pyrimethamine, respectively, for each Z-score unit below -2. Under current dosing recommendations, simulation predicted that the median day 7 concentration was below the 25th percentile for a typical adult patient (50 kg) for sulfadoxine for patients in the weight bands of 8 to 9, 19 to 24, 46 to 49, and 74 to 79 kg and for pyrimethamine for patients in the weight bands of 8 to 9, 14 to 24, and 42 to 49 kg. An evidence-based dosing regimen was constructed that would achieve sulfadoxine and pyrimethamine exposures in young children and underweight-for-age young children that were similar to those currently seen in a typical adult.

Effect of diurnal variation, CYP2B6 genotype and age on the pharmacokinetics of nevirapine in African children

OBJECTIVES

To characterize the effects of CYP2B6 polymorphisms, diurnal variation and demographic factors on nevirapine pharmacokinetics in African children.

METHODS

Non-linear mixed-effects modelling conducted in NONMEM 7.3 described nevirapine plasma concentration-time data from 414 children aged 0.3-15 years.

RESULTS

Nevirapine pharmacokinetics was best described using a one-compartment disposition model with elimination through a well-stirred liver model accounting for a first-pass effect and transit-compartment absorption. Intrinsic clearance was affected by diurnal variation (characterized using a cosine function with peak amplitude 29% at 12 noon) and CYP2B6 metabolizer status [extensive metabolizer (EM) 516GG|983TT, reference; intermediate metabolizer (IM) 516GT|983TT or 516GG|983TC, 17% lower; slow metabolizer (SM) 516TT|983TT or 516GT|983TC, 50% lower; ultra-slow metabolizer (USM) 516GG|983CC, 68% lower]. Age was found to affect pre-hepatic bioavailability: 31.7% lower at birth and increasing exponentially. Median (90% CI) evening Cmin values in the different metabolizer groups were 5.01 (3.01-7.47), 6.55 (3.65-13.32), 11.59 (5.44-22.71) and 12.32 (12.32-27.25) mg/L, respectively. Evening Cmin values were <3 mg/L in 43% of EM weighing <6 kg and 26% of IM weighing <6 kg, while 73% of SM and 88% of USM in all weight-bands had evening Cmin values >8 mg/L. Cmin was not markedly affected by administration time, but was altered by unequal splitting of the daily dose.

CONCLUSIONS

Diurnal variation does not greatly affect nevirapine exposure. However, when daily doses cannot be split equally, the larger dose should be given in the morning. To achieve homogeneous exposures, nevirapine doses for SM and USM should be reduced by 50%, and children weighing <6 kg with EM or IM metabolizer status should receive the same dose as children weighing 6-10 kg.

Optimization of the strength of the efavirenz/lamivudine/abacavir fixed-dose combination for paediatric patients

BACKGROUND

Child-friendly, low-cost, solid, oral fixed-dose combinations (FDCs) of efavirenz with lamivudine and abacavir are urgently needed to improve clinical management and drug adherence for children.

METHODS

Data were pooled from several clinical trials and therapeutic drug monitoring datasets from different countries. The number of children/observations was 505/3667 for efavirenz. Population pharmacokinetic analyses were performed using a non-linear mixed-effects approach. For abacavir and lamivudine, data from 187 and 920 subjects were available (population pharmacokinetic models previously published). Efavirenz/lamivudine/abacavir FDC strength options assessed were (I) 150/75/150, (II) 120/60/120 and (III) 200/100/200 mg. Monte Carlo simulations of the different FDC strengths were performed to determine the optimal dose within each of the WHO weight bands based on drug efficacy/safety targets.

RESULTS

The probability of being within the target efavirenz concentration range 12 h post-dose (1-4 mg/L) varied between 56% and 60%, regardless of FDC option. Option I provided a best possible balance between efavirenz treatment failure and toxicity risks. For abacavir and lamivudine, simulations showed that for option I >75% of subjects were above the efficacy target.

CONCLUSIONS

According to simulations, a paediatric efavirenz/lamivudine/abacavir fixed-dose formulation of 150 mg efavirenz, 75 mg lamivudine and 150 mg abacavir provided the most effective and safe concentrations across WHO weight bands, with the flexibility of dosage required across the paediatric population.

The impact of genetic polymorphisms on the pharmacokinetics of efavirenz in African children

Aims

Using a model‐based approach, the efavirenz steady‐state pharmacokinetics in African children is characterized, quantifying demographic and genotypic effects on the drug's disposition. Simulations are also conducted allowing prediction of optimized doses of efavirenz in this population.

Methods

We modelled the steady‐state population pharmacokinetics of efavirenz in Ugandan and Zambian children using nonlinear mixed‐effects modelling. Individual mid‐dose efavirenz concentrations were derived and simulations explored genotype‐based dose optimization strategies.

Results

A two‐compartment model with absorption through transit compartments well described 2086 concentration‐time points in 169 children. The combined effect of single nucleotide polymorphisms (SNPs) 516G>T and 983T>C explained 44.5% and 14.7% of the variability in efavirenz clearance and bioavailability, respectively. The detected frequencies of composite CYP2B6 genotype were 0.33 for 516GG|983TT, 0.35 for 516GT|983TT, 0.06 for 516GG|983TC, 0.18 for 516TT|983TT, 0.07 516GT|983TC and 0.01 for 516GG|983CC. The corresponding estimated clearance rates were 6.94, 4.90, 3.93, 1.92, 1.36, and 0.74 l h⁻¹ for a 15.4 kg child and median (95% CI) observed mid‐dose concentrations 1.55 (0.51–2.94), 2.20 (0.97–4.40), 2.03 (1.19–4.53), 7.55 (2.40–14.74), 7.79 (3.66–24.59) and 18.22 (11.84–22.76) mg l⁻¹, respectively. Simulations showed that wild‐type individuals had exposures at the bottom of therapeutic range, while slower metabolizers were overexposed.

Conclusions

Dosage guidelines for African children should take into consideration the combined effect of SNPs CYP2B6 516G>T and 983T>C.

Improved prediction of tacrolimus concentrations early after kidney transplantation using theory-based pharmacokinetic modelling

Aims

The aim was to develop a theory-based population pharmacokinetic model of tacrolimus in adult kidney transplant recipients and to externally evaluate this model and two previous empirical models.

Methods

Data were obtained from 242 patients with 3100 tacrolimus whole blood concentrations. External evaluation was performed by examining model predictive performance using Bayesian forecasting.

Results

Pharmacokinetic disposition parameters were estimated based on tacrolimus plasma concentrations, predicted from whole blood concentrations, haematocrit and literature values for tacrolimus binding to red blood cells. Disposition parameters were allometrically scaled to fat free mass. Tacrolimus whole blood clearance/bioavailability standardized to haematocrit of 45% and fat free mass of 60 kg was estimated to be 16.1 l h⁻¹ [95% CI 12.6, 18.0 l h⁻¹]. Tacrolimus clearance was 30% higher (95% CI 13, 46%) and bioavailability 18% lower (95% CI 2, 29%) in CYP3A5 expressers compared with non-expressers. An E_max model described decreasing tacrolimus bioavailability with increasing prednisolone dose. The theory-based model was superior to the empirical models during external evaluation displaying a median prediction error of −1.2% (95% CI −3.0, 0.1%). Based on simulation, Bayesian forecasting led to 65% (95% CI 62, 68%) of patients achieving a tacrolimus average steady-state concentration within a suggested acceptable range.

Conclusion

A theory-based population pharmacokinetic model was superior to two empirical models for prediction of tacrolimus concentrations and seemed suitable for Bayesian prediction of tacrolimus doses early after kidney transplantation.

Population pharmacokinetics of nevirapine in Thai HIV-infected patients

BACKGROUND

Nevirapine-based antiretroviral therapy is widely used as a first-line treatment for HIV-infected patients in resource-limited settings. Nevirapine plasma concentration has been shown to be associated with virological response and treatment failure. Therefore, identifying sources of the variability of nevirapine pharmacokinetics is important for dose optimization. The purpose of the current study was to determine the population mean pharmacokinetic parameters and identify factors that influence pharmacokinetic parameters of nevirapine in Thai HIV-infected patients.

METHODS

The model was developed by a non-linear mixed-effects modelling approach using NONMEM. Model validation was performed using bootstrap analysis and external validation. Additionally, nevirapine plasma concentrations of 200 mg twice daily (NVPBID) and 400 mg once daily (NVPOD) were simulated using the final model to investigate the impact of the covariates and different dosage regimens on nevirapine steady state concentrations.

RESULTS

The apparent clearance (CL/F) of nevirapine estimated from this population was 2.51 l/h which is lower than the values previously reported in other populations. The concomitant use of rifampicin increased CL/F by 20%. Simulated nevirapine plasma concentrations from NVPBID were superior to the NVPOD regimen.

CONCLUSIONS

This population-based pharmacokinetic model can be used for optimizing nevirapine dosage regimens for individual patients to improve efficacy and safety of nevirapine therapy in this population.