A semiphysiological pharmacokinetic model for artemisinin in healthy subjects incorporating autoinduction of metabolism and saturable first-pass hepatic extraction

A Semi-Mechanistic Population Pharmacokinetic Model of Noscapine in Healthy Subjects Considering Hepatic First-Pass Extraction and CYP2C9 Genotypes

INTRODUCTION

Noscapine is a commonly used cough suppressant, with ongoing research on its anti-inflammatory and anti-tumor properties. The drug has a pronounced pharmacokinetic variability.

OBJECTIVE

This evaluation aims to describe the pharmacokinetics of noscapine using a semi-mechanistic population pharmacokinetic model and to identify covariates that could explain inter-individual pharmacokinetic variability.

METHODS

Forty-eight healthy volunteers (30 men and 18 women, mean age 33 years) were enrolled in a randomized, two-period, two-stage, crossover bioequivalence study of noscapine in two different liquid formulations. Noscapine plasma concentrations following oral administration of noscapine 50 mg were evaluated by a non-compartmental analysis and by a population pharmacokinetic model separately.

RESULTS

Compared to the reference formulation, the test formulation exhibited ratios (with 94.12% confidence intervals) of 0.784 (0.662-0.929) and 0.827 (0.762-0.925) for peak plasma concentrations and area under the plasma concentration-time curve, respectively. Significant differences in p values (< 0.01) were both observed when comparing peak plasma concentrations and area under the plasma concentration-time curve between CYP2C9 genotype-predicted phenotypes. A three-compartmental model with zero-order absorption and first-order elimination process best described the plasma data. The introduction of a liver compartment was able to describe the profound first-pass effect of noscapine. Total body weight and the CYP2C9 genotype-predicted phenotype were both identified as significant covariates on apparent clearance, which was estimated as 958 ± 548 L/h for extensive metabolizers (CYP2C9*1/*1 and *1/*9), 531 ± 304 L/h for intermediate metabolizers with an activity score of 1.5 (CYP2C9*1/*2), and 343 ± 197 L/h for poor metabolizers and intermediate metabolizers with an activity score of 1.0 (CYP2C9*1/*3, *2/*3, and*3/*3).

CONCLUSION

The current work is expected to facilitate the future pharmacokinetic/pharmacodynamic development of noscapine. This study was registered prior to starting at "Deutsches Register Klinischer Studien" under registration no. DRKS00017760.

CYP2D6 Phenotype Influences Pharmacokinetic Parameters of Venlafaxine: Results from a Population Pharmacokinetic Model in Older Adults with Depression

In this study, we aimed to improve upon a published population pharmacokinetic (PK) model for venlafaxine (VEN) in the treatment of depression in older adults, then investigate whether CYP2D6 metabolizer status affected model-estimated PK parameters of VEN and its active metabolite O-desmethylvenlafaxine. The model included 325 participants from a clinical trial in which older adults with depression were treated with open-label VEN (maximum 300 mg/day) for 12 weeks and plasma levels of VEN and O-desmethylvenlafaxine were assessed at weeks 4 and 12. We fitted a nonlinear mixed-effect PK model using NONMEM to estimate PK parameters for VEN and O-desmethylvenlafaxine adjusted for CYP2D6 metabolizer status and age. At both lower doses (up to 150 mg/day) and higher doses (up to 300 mg/day), CYP2D6 metabolizers impacted PK model-estimated VEN clearance, VEN exposure, and active moiety (VEN + O-desmethylvenlafaxine) exposure. Specifically, compared with CYP2D6 normal metabolizers, (i) CYP2D6 ultra-rapid metabolizers had higher VEN clearance; (ii) CYP2D6 intermediate metabolizers had lower VEN clearance; (iii) CYP2D6 poor metabolizers had lower VEN clearance, higher VEN exposure, and higher active moiety exposure. Overall, our study showed that including a pharmacogenetic factor in a population PK model could increase model fit, and this improved model demonstrated how CYP2D6 metabolizer status affected VEN-related PK parameters, highlighting the importance of genetic factors in personalized medicine.

Population pharmacokinetic rationale for intravenous contezolid acefosamil followed by oral contezolid dosage regimens

Contezolid is a novel oxazolidinone antibiotic with a promising safety profile. Oral contezolid and its intravenous (IV) prodrug contezolid acefosamil (CZA) are in development for treatment of diabetic foot and acute bacterial skin and skin structure infections (ABSSSI). The prodrug CZA is converted to active contezolid via intermediate MRX-1352. This study aimed to provide the pharmacokinetic rationale for safe, effective, and flexible dosage regimens with initial IV CZA followed by oral contezolid. We simultaneously modeled plasma concentrations from 110 healthy volunteers and 74 phase 2 patients with ABSSSI via population pharmacokinetics (using the importance sampling estimation algorithm), and optimized dosage regimens by Monte Carlo simulations. This included data on MRX-1352, contezolid, and its metabolite MRX-1320 from 66 healthy volunteers receiving intravenous CZA (150-2400 mg) for up to 28 days, and 74 patients receiving oral contezolid [800 mg every 12 h (q12h)] for 10 days. The apparent total clearance for 800 mg oral contezolid with food was 16.0 L/h (23.4% coefficient of variation) in healthy volunteers and 17.7 L/h (53.8%) in patients. CZA was rapidly converted to MRX-1352, which subsequently transformed to contezolid. The proposed dosage regimen used an IV CZA 2000 mg loading dose with 1000 mg IV CZA q12h as maintenance dose(s), followed by 800 mg oral contezolid q12h (with food). During each 24-h period, Monte Carlo simulations predicted this regimen to achieve consistent areas under the curve of 91.9 mg·h/L (range: 76.3-106 mg·h/L) under all scenarios. Thus, this regimen was predicted to reliably achieve efficacious contezolid exposures independent of timing of switch from IV CZA to oral contezolid.IMPORTANCEThis study provides the population pharmacokinetic rationale for the dosage regimen of the intravenous (IV) prodrug contezolid acefosamil (CZA) followed by oral contezolid. We developed the first integrated population model for the pharmacokinetics of the MRX-1352 intermediate prodrug, active contezolid, and its main metabolite MRX-1320 based on data from three clinical studies in healthy volunteers and phase 2 patients. The proposed regimen was predicted to reliably achieve efficacious contezolid exposures independent of timing of switch from IV CZA to oral contezolid.

Dose assessment of melatonin in sepsis (DAMSEL2) study: Pharmacokinetics of two doses of oral melatonin in patients with sepsis

Sepsis is defined as a dysregulated host response to infection, and high-dose melatonin has been proposed as a treatment due to its antioxidant and anti-inflammatory properties. However, there are no data describing the pharmacokinetics of high-dose oral melatonin in critically ill patients. We undertook an open-label trial to determine the tolerance of melatonin administration in these patients and pharmacokinetic analysis, to inform a planned randomised controlled trial. Two cohorts of critically ill patients with sepsis due to community-acquired pneumonia received either 20 or 50 mg oral melatonin liquid as a single dose. Blood samples and clinical measures were analysed over the next 24 h. Melatonin was well tolerated and there were no adverse events. Pharmacokinetic modelling showed that a semiphysiological model, which incorporates saturable first-pass hepatic extraction, was a good fit for our data. Maximum levels of melatonin were extremely high in patients receiving the 50 mg dose and levels of the major metabolite were much lower than expected and not different from those seen after 20 mg, suggesting saturation at the higher dose. We conclude that 20 mg seems a suitable dose of liquid melatonin in patients with sepsis.

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.

Effect of autoinduction and food on the pharmacokinetics of furmonertinib and its active metabolite characterized by a population pharmacokinetic model

Furmonertinib (AST2818) is a novel third-generation irreversible EGFR TKI and recently has been approved in China for the treatment of non-small cell lung cancer (NSCLC) with EGFR-sensitizing and T790M resistance mutations. In the current study, we developed a semi-mechanistic population pharmacokinetic model to characterize the nonstationary pharmacokinetics (PK) of the furmonertinib and its active metabolite AST5902 simultaneously. The PK data of furmonertinib and AST5902 were obtained from 38 NSCLC patients and 16 healthy volunteers receiving 20-240 mg furmonertinib in three clinical trials. A nonlinear mixed-effects modeling approach was used to describe the PK data. The absorption process of furmonertinib was described by a transit compartment model. The disposition of both furmonertinib and AST5902 was described by a two-compartment model. An indirect response model accounted for the autoinduction of furmonertinib metabolism mediated by CYP3A4. The model-based simulation suggested that furmonertinib clearance was increased in one cycle of treatment (orally once daily for 21 days) compared to baseline, ranging from 1.1 to 1.8 fold corresponding to the dose range of 20-240 mg. The concentration of furmonertinib was decreased over time whereas that of AST5902 was increased. Interestingly, the concentration of the total active compounds (furmonertinib and AST5902) appeared to be stable. The food intake, serum alkaline phosphatase and body weight were identified as statistically significant covariates. The mechanism of food effect on PK was investigated, where the food intake might increase the bioavailability of furmonertinib via increasing the splanchnic blood flow. Overall, a population PK model was successfully developed to characterize the nonstationary PK of furmonertinib and AST5902 simultaneously. The concentrations of total active compounds were less affected by the autoinduction of furmonertinib metabolism.

Population Pharmacokinetic Model and Alternative Dosing Regimens for Dolutegravir Coadministered with Rifampicin

Dolutegravir-based regimens are recommended as first-line therapy for HIV in low- and middle-income countries where tuberculosis is the most common opportunistic infection. Concurrent HIV/tuberculosis treatment is challenging because of drug-drug interactions. Our analysis aimed to characterize dolutegravir's population pharmacokinetics when coadministered with rifampicin and assess alternative dolutegravir dosing regimens. We developed a population pharmacokinetic model of dolutegravir in NONMEM with data from two healthy-volunteer studies (RADIO and ClinicalTrials.gov identifier NCT01231542) and validated it with data from the INSPIRING study, which consisted of participants living with HIV. The model was developed with 817 dolutegravir plasma concentrations from 41 participants. A 2-compartment model with first-order elimination and lagged absorption best described dolutegravir's pharmacokinetics. For a typical 70-kg individual, we estimated a clearance, absorption rate constant, central volume, and peripheral volume of 1.03 L/h, 1.61 h^-1, 12.7 L, and 3.85 L, respectively. Rifampicin coadministration increased dolutegravir clearance by 144% (95% confidence interval [CI], 126 to 161%). Simulations showed that when 50 or 100 mg once-daily dolutegravir is coadministered with rifampicin in 70-kg individuals, 71.7% and 91.5% attain trough concentrations above 0.064 mg/L, the protein-adjusted 90% inhibitory concentration (PA-IC₉₀), respectively. The model developed from healthy-volunteer data describes patient data reasonably well but underpredicts trough concentrations. Although 50 mg of dolutegravir given twice daily achieves target concentrations in more than 99% of individuals cotreated with rifampicin, 100 mg of dolutegravir, once daily, in the same population is predicted to achieve satisfactory pharmacokinetic target attainment. The efficacy of this regimen should be investigated since it presents an opportunity for treatment simplification.

Pharmacokinetics and Dose Optimization Strategies of Para-Aminosalicylic Acid in Children with Rifampicin-Resistant Tuberculosis

Treatment options for children with Rifampicin-resistant tuberculosis (RR-TB) remain limited, and para-aminosalicylic acid (PAS) is still a relevant component of treatment regimens. Prevention of resistance to companion drugs by PAS is dose related, and at higher concentrations, PAS may exhibit significant bactericidal activity in addition to its bacteriostatic properties. The optimal dosing of PAS in children is uncertain, specifically for delayed-release granule preparations, which are the most used. A population pharmacokinetic model was developed describing PAS pharmacokinetics in children receiving routine RR-TB treatment. Model-based simulations evaluated current World Health Organization (WHO) weight-band doses against the adult pharmacokinetic target of 50 to 100 mg/liter for peak concentrations. Of 27 children included, the median (range) age and weight were 3.87 (0.58 to 13.7) years and 13.3 (7.15 to 30.5) kg, respectively; 4 (14.8%) were HIV positive. PAS followed one-compartment kinetics with first-order elimination and transit compartment absorption. The typical clearance in a 13-kg child was 9.79 liters/h. Increased PAS clearance was observed in both pharmacokinetic profiles from the only patient receiving efavirenz. No effect of renal function, sex, ethnicity, nutritional status, HIV status, antiretrovirals (lamivudine, abacavir, and lopinavir-ritonavir), or RR-TB drugs was detected. In simulations, target concentrations were achieved only using the higher WHO dose range of 300 mg/kg once daily. A transit compartment adequately describes absorption for the slow-release PAS formulation. Children should be dosed at the higher range of current WHO-recommended PAS doses and in a once-daily dose to optimize treatment.

Population Pharmacokinetics of Antimalarial Naphthoquine in Combination with Artemisinin in Tanzanian Children and Adults: Dose Optimization

The combination antimalarial therapy of artemisinin-naphthoquine (ART-NQ) was developed as a single-dose therapy, aiming to improve adherence relative to the multiday schedules of other artemisinin combination therapies. The pharmacokinetics of ART-NQ has not been well characterized, especially in children. A pharmacokinetic study was conducted in adults and children over 5 years of age (6 to 10, 11 to 17, and ≥18 years of age) with uncomplicated malaria in Tanzania. The median weights for the three age groups were 20, 37.5, and 55 kg, respectively. Twenty-nine patients received single doses of 20 mg/kg of body weight for artemisinin and 8 mg/kg for naphthoquine, and plasma drug concentrations were assessed at 13 time points over 42 days from treatment. We used nonlinear mixed-effects modeling to interpret the data, and allometric scaling was employed to adjust for the effect of body size. The pharmacokinetics of artemisinin was best described by one-compartment model and that of naphthoquine by a two-compartment disposition model. Clearance values for a typical patient (55-kg body weight and 44.3-kg fat-free mass) were estimated as 66.7 L/h (95% confidence interval [CI], 57.3 to 78.5 L/h) for artemisinin and 44.2 L/h (95% CI, 37.9 to 50.6 L/h) for naphthoquine. Nevertheless, we show via simulation that patients weighing ≥70 kg achieve on average a 30% lower day 7 concentration compared to a 48-kg reference patient at the doses tested, suggesting dose increases may be warranted to ensure adequate exposure. (This study has been registered at ClinicalTrials.gov under identifier NCT01930331.).

Posaconazole bioavailability of the solid oral tablet is reduced during severe intestinal mucositis

OBJECTIVES

To describe the absolute oral bioavailability of the solid oral formulation of posaconazole and the impact of severe intestinal mucositis in haematology patients. Further, to describe posaconazole protein binding in haematology patients.

METHODS

A pharmacokinetic study was performed where patients receiving induction chemotherapy or a haematopoietic cell transplantation were randomized to receive 7 days of intravenous posaconazole therapy followed by 9 days of oral therapy or vice versa. Patients received posaconazole licensed dose until day 12, thereafter a reduced once daily dose of 200 mg was given. At days 7, 12, and 16 blood samples were obtained for pharmacokinetic curves, trough samples were collected on all other days. Total and unbound posaconazole pharmacokinetics were analysed by population pharmacokinetic modelling. The presence of severe intestinal mucositis was assessed by plasma citrulline levels and analysed as a binary covariate using 10 μmol/L as cut-off. Monte Carlo simulations were performed to simulate posaconazole exposure at steady state.

RESULTS

Twenty-three patients were included for analysis with 581 total posaconazole concentrations and 91 paired unbound concentrations. Absolute bioavailability in the final model was estimated at 51.4% (%relative standard error (RSE) 56.5) and 67.6% (%RSE 75.0) in patients with and without severe intestinal mucositis, respectively. Posaconazole unbound fraction was estimated at 2.7% (%RSE 3.9).

CONCLUSIONS

Posaconazole bioavailability is reduced in haematological patients with severe intestinal mucositis, requiring an increase in oral posaconazole dose to 400 mg twice daily on day 1 followed by 400 mg once daily or a switch to intravenous therapy.

Artemisinin Cocrystals for Bioavailability Enhancement. Part 2: In Vivo Bioavailability and Physiologically Based Pharmacokinetic Modeling

We report the evaluation and prediction of the pharmacokinetic (PK) performance of artemisinin (ART) cocrystal formulations, that is, 1:1 artemisinin/orcinol (ART-ORC) and 2:1 artemisinin/resorcinol (ART²-RES), using in vivo murine animal and physiologically based pharmacokinetic (PBPK) models. The efficacy of the ART cocrystal formulations along with the parent drug ART was tested in mice infected with Plasmodium berghei. When given at the same dose, the ART cocrystal formulation showed a significant reduction in parasitaemia at day 4 after infection compared to ART alone. PK parameters including C_max (maximum plasma concentration), T_max (time to C_max), and AUC (area under the curve) were obtained by determining drug concentrations in the plasma using liquid chromatography-high-resolution mass spectrometry (LC-HRMS), showing enhanced ART levels after dosage with the cocrystal formulations. The dose-response tests revealed that a significantly lower dose of the ART cocrystals in the formulation was required to achieve a similar therapeutic effect as ART alone. A PBPK model was developed using a PBPK mouse simulator to accurately predict the in vivo behavior of the cocrystal formulations by combining in vitro dissolution profiles with the properties of the parent drug ART. The study illustrated that information from classical in vitro and in vivo experimental investigations of the parent drug of ART formulations can be coupled with PBPK modeling to predict the PK parameters of an ART cocrystal formulation in an efficient manner. Therefore, the proposed modeling strategy could be used to establish in vitro and in vivo correlations for different cocrystals intended to improve dissolution properties and to support clinical candidate selection, contributing to the assessment of cocrystal developability and formulation development.

Population pharmacokinetic modeling of molibresib and its active metabolites in patients with solid tumors: A semimechanistic autoinduction model

Molibresib (GSK525762) is an investigational, orally bioavailable, small‐molecule bromodomain and extraterminal (BET) protein inhibitor for the treatment of advanced solid tumors. Molibresib was initially evaluated in a first‐time‐in‐human (FTIH) study BET115521 consisting of two parts: Part 1 of the study (dose escalation) was conducted in 94 patients with nuclear protein in testis midline carcinoma and other solid tumors, and Part 2 (expansion cohort) was conducted in 99 patients with different solid tumor types. Molibresib is metabolized by cytochrome P450 3A4 enzymes to produce two major active metabolites that are equipotent to the parent molecule. The metabolites are measured together after full conversion of one to the other and reported as an active metabolite composite (GSK3529246). The molibresib pharmacokinetic (PK) profile has been characterized by a decrease in exposure over time, with the decrease more pronounced at higher doses, and accompanied by a slight increase of the metabolite concentrations. Autoinduction of molibresib metabolism was suspected and confirmed in vitro. Here we report the development of a semimechanistic liver‐compartment population PK model using PK data from the FTIH study, which adequately describes the autoinduction of molibresib clearance and the PK of both molibresib and GSK3529246. Covariate analysis indicated body weight had a significant effect on the volume of distribution of molibresib and GSK3529246, and higher levels of aspartate aminotransferase resulted in the lower clearance of GSK3529246. This model was used to simulate individual patient exposures based on covariate information for use in future alternative dosing strategies and exposure–response analyses.

Mycobacterium tuberculosis precursor rRNA as a measure of treatment-shortening activity of drugs and regimens

There is urgent need for new drug regimens that more rapidly cure tuberculosis (TB). Existing TB drugs and regimens vary in treatment-shortening activity, but the molecular basis of these differences is unclear, and no existing assay directly quantifies the ability of a drug or regimen to shorten treatment. Here, we show that drugs historically classified as sterilizing and non-sterilizing have distinct impacts on a fundamental aspect of Mycobacterium tuberculosis physiology: ribosomal RNA (rRNA) synthesis. In culture, in mice, and in human studies, measurement of precursor rRNA reveals that sterilizing drugs and highly effective drug regimens profoundly suppress M. tuberculosis rRNA synthesis, whereas non-sterilizing drugs and weaker regimens do not. The rRNA synthesis ratio provides a readout of drug effect that is orthogonal to traditional measures of bacterial burden. We propose that this metric of drug activity may accelerate the development of shorter TB regimens.

Pharmacokinetics and Drug-Drug Interactions of Isoniazid and Efavirenz in Pregnant Women Living With HIV in High TB Incidence Settings: Importance of Genotyping

The World Health Organization guidelines recommend that individuals living with HIV receive ≥ 6 months of isoniazid preventive therapy, including pregnant women. Yet, plasma isoniazid exposure during pregnancy, in the antiretroviral therapy era, has not been well-described. We investigated pregnancy-induced and pharmacogenetic-associated pharmacokinetic changes and drug-drug interactions between isoniazid and efavirenz in pregnant women. Eight hundred forty-seven women received isoniazid for 28 weeks, either during pregnancy or at 12 weeks postpartum, and 786 women received efavirenz. After adjusting for NAT2 and CYP2B6 genotype and weight, pregnancy increased isoniazid and efavirenz clearance by 26% and 15%, respectively. Isoniazid decreased efavirenz clearance by 7% in CYP2B6 normal metabolizers and 13% in slow and intermediate metabolizers. Overall, both isoniazid and efavirenz exposures were reduced during pregnancy, but the main determinants of drug concentration were NAT2 and CYP2B6 genotypes, which resulted in a five-fold difference for both drugs between rapid and slow metabolizers.

QSAR-Based Virtual Screening of Natural Products Database for Identification of Potent Antimalarial Hits

With about 400,000 annual deaths worldwide, malaria remains a public health burden in tropical and subtropical areas, especially in low-income countries. Selection of drug-resistant Plasmodium strains has driven the need to explore novel antimalarial compounds with diverse modes of action. In this context, biodiversity has been widely exploited as a resourceful channel of biologically active compounds, as exemplified by antimalarial drugs such as quinine and artemisinin, derived from natural products. Thus, combining a natural product library and quantitative structure-activity relationship (QSAR)-based virtual screening, we have prioritized genuine and derivative natural compounds with potential antimalarial activity prior to in vitro testing. Experimental validation against cultured chloroquine-sensitive and multi-drug-resistant P. falciparum strains confirmed the potent and selective activity of two sesquiterpene lactones (LDT-597 and LDT-598) identified in silico. Quantitative structure-property relationship (QSPR) models predicted absorption, distribution, metabolism, and excretion (ADME) and physiologically based pharmacokinetic (PBPK) parameters for the most promising compound, showing that it presents good physiologically based pharmacokinetic properties both in rats and humans. Altogether, the in vitro parasite growth inhibition results obtained from in silico screened compounds encourage the use of virtual screening campaigns for identification of promising natural compound-based antimalarial molecules.

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.

Predicting the Disposition of the Antimalarial Drug Artesunate and Its Active Metabolite Dihydroartemisinin Using Physiologically Based Pharmacokinetic Modeling

Artemisinin-based combination therapies (ACTs) have proven to be effective in helping to combat the global malaria epidemic. To optimally apply these drugs, information about their tissue-specific disposition is required, and one approach to predict these pharmacokinetic characteristics is physiologically based pharmacokinetic (PBPK) modeling. In this study, a whole-body PBPK model was developed to simulate the time-dependent tissue concentrations of artesunate (AS) and its active metabolite, dihydroartemisinin (DHA). The model was developed for both rats and humans and incorporated drug metabolism of the parent compound and major metabolite. Model calibration was conducted using data from the literature in a Bayesian framework, and model verification was assessed using separate sets of data. Results showed good agreement between model predictions and the validation data, demonstrating the capability of the model in predicting the blood, plasma, and tissue pharmacokinetics of AS and DHA. It is expected that such a tool will be useful in characterizing the disposition of these chemicals and ultimately improve dosing regimens by enabling a quantitative assessment of the tissue-specific drug levels critical in the evaluation of efficacy and toxicity.

Auto-Induction of Intestinal First-Pass Effect Related Time-Dependent Pharmacokinetics of Artemisinin Rather than Dihydroartemisinin

Artemisinin (ART) drugs showed declining plasma concentrations after repeated oral dosing, known as time-dependent pharmacokinetics (PK). ART and dihydroartemisinin (DHA) were adopted as representatives to evaluate the roles of first-pass effects and systemic metabolism in time-dependent PK by comparison of oral versus intravenous administration and 1 dose versus 5 consecutive doses PK in rats and dogs, respectively. The hepatic extraction ratio (ER_h) and the intestinal elimination changes were further investigated in rats to distinguish the roles of hepatic first-pass effect or intestinal first-pass effect. The induction capacities of ARTs to cytochrome P450 (CYP450) in rats and human cells were evaluated as well. For ART, only the oral groups showed time-dependent PK. A fairly high ER_h that obtained for ART was not sensitive to multiple oral doses. An increased elimination and CYP450 expression have also been found in the intestine. For DHA, though a significant CYP450 induction was observed, neither time-dependent PK nor changes in the first-pass effects was found. In conclusion, time-dependent PK of ART was mainly caused by the increased intestinal first-pass effect rather than hepatic first-pass effect or systemic metabolism. DHA was not involved in auto-induction elimination, thus showing no time-dependent PK.

A semiphysiological population pharmacokinetic model of agomelatine and its metabolites in Chinese healthy volunteers

AIMS

Agomelatine is an antidepressant for major depressive disorders. It undergoes extensive first-pass hepatic metabolism and displays irregular absorption profiles and large interindividual variability (IIV) and interoccasion variability of pharmacokinetics. The objective of this study was to characterize the complex pharmacokinetics of agomelatine and its metabolites in healthy subjects.

METHODS

Plasma concentration-time data of agomelatine and its metabolites were collected from a 4-period, cross-over bioequivalence study, in which 44 healthy subjects received 25 mg agomelatine tablets orally. Nonlinear mixed effects modelling was used to characterize the pharmacokinetics and variability of agomelatine and its metabolites. Deterministic simulations were carried out to investigate the influence of pathological changes due to liver disease on agomelatine pharmacokinetics.

RESULTS

A semiphysiological pharmacokinetic model with parallel first-order absorption and a well-stirred liver compartment adequately described the data. The estimated IIV and interoccasion variability of the intrinsic clearance of agomelatine were 130.8% and 28.5%, respectively. The IIV of the intrinsic clearance turned out to be the main cause of the variability of area under the curve-based agomelatine exposure. Simulations demonstrated that a reduction in intrinsic clearance or liver blood flow, and an increase in free drug fraction had a rather modest influence on agomelatine exposures (range: -50 to 200%). Portosystemic shunting, however, substantially elevated agomelatine exposure by 12.6-109.1-fold.

CONCLUSIONS

A semiphysiological pharmacokinetic model incorporating first-pass hepatic extraction was developed for agomelatine and its main metabolites. The portosystemic shunting associated with liver disease might lead to significant alterations of agomelatine pharmacokinetics, and lead to substantially increased exposure.

A Microbial Transformation Model for Simulating Mammal Metabolism of Artemisinin

Artemisinin (ART) is a highly effective antimalarial agent isolated from the traditional Chinese herb Qinghao. Metabolism of ART and its derivatives in the body is one of the most pressing issues for pharmaceutical scientists. Herein, an efficient in vitro microorganism model for simulation of metabolism of ART in vivo was developed employing Cunninghamella elegans. Metabolites in the microbial transformation system and plasma of mice pre-administrated ART orally were analyzed by ultra-performance liquid chromatography (UPLC)-electrospray ionization (ESI)-quadrupole time-of-flight (Q-TOF)-mass spectrometry (MS^E) combined with UNIFI software. Thirty-two metabolites were identified in vitro and 23 were identified in vivo. After comparison, 16 products were found to be common to both models including monohydroxylated ART, dihydroxylated ART, deoxyartemisinin, hydroxylated deoxyartemisinin, hydroxylated dihydroartemisinin (DHA), and hydroxylated deoxy-DHA. These results revealed that C. elegans CICC 40250 functioned as an appropriate model to mimic ART metabolism in vivo. Moreover, an overall description of metabolites of ART from C. elegans CICC 40250 has been provided. Notably, DHA was detected and identified as a metabolite of ART in mouse plasma for the first time.

Effect of efavirenz-based antiretroviral therapy and high-dose rifampicin on the pharmacokinetics of isoniazid and acetyl-isoniazid

Objectives

To describe the pharmacokinetics of isoniazid and acetyl-isoniazid in TB/HIV-coinfected patients, and assess the effects of efavirenz co-administration and a 50% increase in the dose of rifampicin on the pharmacokinetic parameters of isoniazid and acetyl-isoniazid.

Methods

TB/HIV-coinfected patients participating in the three-treatment-arm RAFA randomized controlled trial conducted in West Africa were recruited into the pharmacokinetics sub-study. Five serial blood samples were collected on a single visit between 4 and 8 weeks after initiation of antituberculosis treatment. Concentration-time data for isoniazid and acetyl-isoniazid were analysed using non-linear mixed-effects models.

Results

Isoniazid concentrations from 150 patients were available for analysis, and 79 of these (53%) also had concentrations of acetyl-isoniazid. Isoniazid pharmacokinetics was best described with a two-compartment disposition model with lagged first-order absorption and elimination using a semi-mechanistic model describing hepatic extraction. The model identified two elimination pathways, separating formation of acetyl-isoniazid from other routes of metabolism. The predicted AUC0-24 is reduced by 29% in patients who are fast acetylators of isoniazid and receiving efavirenz-based ART (6.73 versus 4.68 mg·h/L). In slow acetylators, efavirenz-based ART had no effect on isoniazid exposure (AUC0-24 = 17.5 mg·h/L).

Conclusions

Efavirenz-based ART affects the acetylation metabolic pathway amongst rapid acetylators, resulting in reduced exposure to isoniazid. Pharmacokinetics of isoniazid and acetyl-isoniazid were not influenced by the 50% increase in rifampicin dose.

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 Artemether, Dihydroartemisinin, and Lumefantrine in Rwandese Pregnant Women Treated for Uncomplicated Plasmodium falciparum Malaria

The artemisinin-based combination therapy artemether-lumefantrine is commonly used in pregnant malaria patients. However, the effect of pregnancy-related changes on exposure is unclear, and pregnancy has been associated with decreased efficacy in previous studies. This study aimed to characterize the population pharmacokinetics of artemether, its active metabolite dihydroartemisinin, and lumefantrine in 22 Rwandese pregnant women in their second (n = 11) or third (n = 11) trimester with uncomplicated Plasmodium falciparum malaria. These patients were enrolled from Rwamagana district hospital and received the standard fixed oral dose combination of 80 mg of artemether and 480 mg of lumefantrine twice daily for 3 days. Venous plasma concentrations were quantified for all three analytes using liquid chromatography coupled with tandem mass spectroscopy, and data were analyzed using nonlinear mixed-effects modeling. Lumefantrine pharmacokinetics was described by a flexible but highly variable absorption, with a mean absorption time of 4.04 h, followed by a biphasic disposition model. The median area under the concentration-time curve from 0 h to infinity (AUC_0-∞) for lumefantrine was 641 h · mg/liter. Model-based simulations indicated that 11.7% of the study population did not attain the target day 7 plasma concentration (280 ng/ml), a threshold associated with increased risk of recrudescence. The pharmacokinetics of artemether was time dependent, and the autoinduction of its clearance was described using an enzyme turnover model. The turnover half-life was predicted to be 30.4 h. The typical oral clearance, which started at 467 liters/h, increased 1.43-fold at the end of treatment. Simulations suggested that lumefantrine pharmacokinetic target attainment appeared to be reassuring in Rwandese pregnant women, particularly compared to target attainment in Southeast Asia. Larger cohorts will be required to confirm this finding.

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.

Effect of rifampicin and efavirenz on moxifloxacin concentrations when co-administered in patients with drug-susceptible TB

Objectives

We compared the pharmacokinetics of moxifloxacin during rifampicin co-treatment or when dosed alone in African patients with drug-susceptible recurrent TB.

Methods

Patients in the intervention arm of the Improving Retreatment Success (IMPRESS) randomized controlled TB trial received 400 mg of moxifloxacin, with rifampicin, isoniazid and pyrazinamide in the treatment regimen. Moxifloxacin concentrations were measured in plasma during rifampicin-based TB treatment and again 4 weeks after treatment completion, when given alone as a single dose. Moxifloxacin concentration-time data were analysed using non-linear mixed-effects models.

Results

We included 58 patients; 42 (72.4%) were HIV co-infected and 40 (95%) of these were on efavirenz-based ART. Moxifloxacin pharmacokinetics was best described using a two-compartment disposition model with first-order lagged absorption and elimination using a semi-mechanistic model describing hepatic extraction. Oral clearance (CL/F) of moxifloxacin during rifampicin-based TB treatment was 24.3 L/h for a typical patient (fat-free mass of 47 kg), resulting in an AUC of 16.5 mg·h/L. This exposure was 7.8% lower than the AUC following the single dose of moxifloxacin given alone after TB treatment completion. In HIV-co-infected patients taking efavirenz-based ART, CL/F of moxifloxacin was increased by 42.4%, resulting in a further 30% reduction in moxifloxacin AUC.

Conclusions

Moxifloxacin clearance was high and plasma concentrations low in our patients overall. Moxifloxacin AUC was further decreased by co-administration of efavirenz-based ART and, to a lesser extent, rifampicin. The clinical relevance of the low moxifloxacin concentrations for TB treatment outcomes and the need for moxifloxacin dose adjustment in the presence of rifampicin and efavirenz co-treatment need further investigation.

An automated sampling importance resampling procedure for estimating parameter uncertainty

Quantifying the uncertainty around endpoints used for decision-making in drug development is essential. In nonlinear mixed-effects models (NLMEM) analysis, this uncertainty is derived from the uncertainty around model parameters. Different methods to assess parameter uncertainty exist, but scrutiny towards their adequacy is low. In a previous publication, sampling importance resampling (SIR) was proposed as a fast and assumption-light method for the estimation of parameter uncertainty. A non-iterative implementation of SIR proved adequate for a set of simple NLMEM, but the choice of SIR settings remained an issue. This issue was alleviated in the present work through the development of an automated, iterative SIR procedure. The new procedure was tested on 25 real data examples covering a wide range of pharmacokinetic and pharmacodynamic NLMEM featuring continuous and categorical endpoints, with up to 39 estimated parameters and varying data richness. SIR led to appropriate results after 3 iterations on average. SIR was also compared with the covariance matrix, bootstrap and stochastic simulations and estimations (SSE). SIR was about 10 times faster than the bootstrap. SIR led to relative standard errors similar to the covariance matrix and SSE. SIR parameter 95% confidence intervals also displayed similar asymmetry to SSE. In conclusion, the automated SIR procedure was successfully applied over a large variety of cases, and its user-friendly implementation in the PsN program enables an efficient estimation of parameter uncertainty in NLMEM.

A Semi-Physiological Population Model to Quantify the Effect of Hematocrit on Everolimus Pharmacokinetics and Pharmacodynamics in Cancer Patients

Introduction and Objective

Everolimus (a drug from the class of mammalian target of rapamycin [mTOR] inhibitors) is associated with frequent toxicity-related dose reductions. Everolimus accumulates in erythrocytes, but the extent to which hematocrit affects everolimus plasma pharmacokinetics and pharmacodynamics is unknown. We aimed to investigate the everolimus pharmacokinetics/pharmacodynamics and the influence of hematocrit in cancer patients.

Methods

A semi-physiological pharmacokinetic model for everolimus was developed from pharmacokinetic data from 73 patients by non-linear mixed-effects modeling. Using a simulation study with a known pharmacodynamic model describing S6K1 (a downstream mTOR effector) inhibition, we investigated the impact of hematocrit.

Results

The apparent volume of distribution of the central and peripheral compartment were estimated to be 207 L with a relative standard error (RSE) of 5.0 % and 485 L (RSE 4.2 %), respectively, with an inter-compartmental clearance of 72.1 L/h (RSE 3.2 %). The apparent intrinsic clearance was 198 L/h (RSE 4.3 %). A decrease in hematocrit from 45 % to 20 % resulted in a predicted reduction in whole-blood exposure of ~50 %, but everolimus plasma pharmacokinetics and pharmacodynamics were not affected. The predicted S6K1 inhibition was at a plateau level in the approved dose of 10 mg once daily.

Conclusions

A population pharmacokinetic model was developed for everolimus in cancer patients. Hematocrit influenced whole-blood pharmacokinetics, but not plasma pharmacokinetics or pharmacodynamics. Everolimus whole-blood concentrations should always be corrected for hematocrit. Since predicted mTOR inhibition was at a plateau level in the approved dose, dose reductions may have only a limited impact on mTOR inhibition.

Population Pharmacokinetic and Pharmacodynamic Modeling of Artemisinin Resistance in Southeast Asia

Orally administered artemisinin-based combination therapy is the first-line treatment against uncomplicated P. falciparum malaria worldwide. However, the increasing prevalence of artemisinin resistance is threatening efforts to treat and eliminate malaria in Southeast Asia. This study aimed to characterize the exposure-response relationship of artesunate in patients with artemisinin sensitive and resistant malaria infections. Patients were recruited in Pailin, Cambodia (n = 39), and Wang Pha, Thailand (n = 40), and received either 2 mg/kg/day of artesunate mono-therapy for 7 consecutive days or 4 mg/kg/day of artesunate monotherapy for 3 consecutive days followed by mefloquine 15 and 10 mg/kg for 2 consecutive days. Plasma concentrations of artesunate and its active metabolite, dihydroartemisinin, and microscopy-based parasite densities were measured and evaluated using nonlinear mixed-effects modeling. All treatments were well tolerated with minor and transient adverse reactions. Patients in Cambodia had substantially slower parasite clearance compared to patients in Thailand. The pharmacokinetic properties of artesunate and dihydroartemisinin were well described by transit-compartment absorption followed by one-compartment disposition models. Parasite density was a significant covariate, and higher parasite densities were associated with increased absorption. Dihydroartemisinin-dependent parasite killing was described by a delayed sigmoidal Emax model, and a mixture function was implemented to differentiate between sensitive and resistant infections. This predicted that 84% and 16% of infections in Cambodia and Thailand, respectively, were artemisinin resistant. The final model was used to develop a simple diagnostic nomogram to identify patients with artemisinin-resistant infections. The nomogram showed > 80% specificity and sensitivity, and outperformed the current practice of day 3 positivity testing.

Age, Weight, and CYP2D6 Genotype Are Major Determinants of Primaquine Pharmacokinetics in African Children

Low-dose primaquine is recommended to prevent Plasmodium falciparum malaria transmission in areas threatened by artemisinin resistance and areas aiming for malaria elimination. Community treatment campaigns with artemisinin-based combination therapy in combination with the gametocytocidal primaquine dose target all age groups, but no studies thus far have assessed the pharmacokinetics of this gametocytocidal drug in African children. We recruited 40 children participating in a primaquine efficacy trial in Burkina Faso to study primaquine pharmacokinetics. These children received artemether-lumefantrine and either a 0.25- or a 0.40-mg/kg primaquine dose. Seven blood samples were collected from each participant for primaquine and carboxy-primaquine plasma levels determinations: one sample was collected before primaquine administration and six after primaquine administration according to partially overlapping sampling schedules. Physiological population pharmacokinetic modeling was used to assess the impact of weight, age, and CYP2D6 genotype on primaquine and carboxy-primaquine pharmacokinetics. Despite linear weight normalized dosing, the areas under the plasma concentration-time curves and the peak concentrations for both primaquine and carboxy-primaquine increased with age and body weight. Children who were CYP2D6 poor metabolizers had higher levels of the parent compound, indicating a lower primaquine CYP2D6-mediated metabolism. Our data indicate that primaquine and carboxy-primaquine pharmacokinetics are influenced by age, weight, and CYP2D6 genotype and suggest that dosing strategies may have to be reconsidered to maximize the transmission-blocking properties of primaquine. (This study has been registered at ClinicalTrials.gov under registration no. NCT01935882.)

A Mechanism-Based Pharmacokinetic Enzyme Turnover Model for Dichloroacetic Acid Autoinhibition in Rats

Dichloroacetic acid (DCA), a halogenated organic acid, is a pyruvate dehydrogenase kinase inhibitor that has been used to treat congenital or acquired lactic acidosis and is currently in early-phase clinical trials for cancer treatment. DCA was found to inhibit its own metabolism by irreversibly inactivating glutathione transferase zeta 1 (GSTZ1-1), resulting in nonlinear kinetics and abnormally high accumulation ratio after repeated dosing. In this analysis, a semi-mechanistic pharmacokinetic enzyme turnover model was developed for the first time to capture DCA autoinhibition, gastrointestinal region-dependent absorption, and time-dependent change in bioavailability in rats. The maximum rate constant for DCA-induced GSTZ1-1 inactivation is estimated to be 0.96/h, which is 110 times that of the rate constant for GSTZ1-1 natural degradation (0.00875/h). The model-predicted DCA concentration that corresponds to 50% of maximum enzyme inhibition (EC₅₀) is 4.32 mg/L. The constructed pharmacokinetic enzyme turnover model, when applied to human data, could be used to predict the accumulation of DCA after repeated oral dosing, guide selection of dosing regimens in clinical studies, and facilitate clinical development of DCA.

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.

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.

Integrated semi-physiological pharmacokinetic model for both sunitinib and its active metabolite SU12662

AIMS

Previously published pharmacokinetic (PK) models for sunitinib and its active metabolite SU12662 were based on a limited dataset or lacked important elements such as correlations between sunitinib and its metabolite. The current study aimed to develop an improved PK model that circumvented these limitations and to prove the utility of the PK model in treatment optimization in clinical practice.

METHODS

One thousand two hundred and five plasma samples from 70 cancer patients were collected from three PK studies with sunitinib and SU12662. A semi-physiological PK model for sunitinib and SU12662 was developed incorporating pre-systemic metabolism using non-linear mixed effects modelling (nonmem). Allometric scaling based on body weight was applied. The final model was used for simulation of the PK of different treatment regimens.

RESULTS

Sunitinib and SU12662 PK were best described by a one and two compartment model, respectively. Introduction of pre-systemic formation of SU12662 strongly improved model fit, compared with solely systemic metabolism. The clearance of sunitinib and SU12662 was estimated at 35.7 (relative standard error (RSE) 5.7%) l h(-1) and 17.1 (RSE 7.4%) l h(-1), respectively for 70 kg patients. Correlation coefficients were estimated between inter-individual variability of both clearances, both volumes of distribution and between clearance and volume of distribution of SU12662 as 0.53, 0.48 and 0.45, respectively. Simulation of the PK model predicted correctly the ratio of patients who did not reach proposed PK targets for efficacy.

CONCLUSIONS

A semi-physiological PK model for sunitinib and SU12662 in cancer patients was presented including pre-systemic metabolism. The model was superior to previous PK models in many aspects.

Model-Based Evaluation of Higher Doses of Rifampin Using a Semimechanistic Model Incorporating Autoinduction and Saturation of Hepatic Extraction

Rifampin is a key sterilizing drug in the treatment of tuberculosis (TB). It induces its own metabolism, but neither the onset nor the extent of autoinduction has been adequately described. Currently, the World Health Organization recommends a rifampin dose of 8 to 12 mg/kg of body weight, which is believed to be suboptimal, and higher doses may potentially improve treatment outcomes. However, a nonlinear increase in exposure may be observed because of saturation of hepatic extraction and hence this should be taken into consideration when a dose increase is implemented. Intensive pharmacokinetic (PK) data from 61 HIV-TB-coinfected patients in South Africa were collected at four visits, on days 1, 8, 15, and 29, after initiation of treatment. Data were analyzed by population nonlinear mixed-effects modeling. Rifampin PKs were best described by using a transit compartment absorption and a well-stirred liver model with saturation of hepatic extraction, including a first-pass effect. Autoinduction was characterized by using an exponential-maturation model: hepatic clearance almost doubled from the baseline to steady state, with a half-life of around 4.5 days. The model predicts that increases in the dose of rifampin result in more-than-linear drug exposure increases as measured by the 24-h area under the concentration-time curve. Simulations with doses of up to 35 mg/kg produced results closely in line with those of clinical trials.

Development of a physiology-directed population pharmacokinetic and pharmacodynamic model for characterizing the impact of genetic and demographic factors on clopidogrel response in healthy adults

Clopidogrel (Plavix®), is a widely used antiplatelet agent, which shows high inter-individual variability in treatment response in patients following the standard dosing regimen. In this study, a physiology-directed population pharmacokinetic/pharmacodynamic (PK/PD) model was developed based on clopidogrel and clopidogrel active metabolite (clop-AM) data from the PAPI and the PGXB2B studies using a step-wise approach in NONMEM (version 7.2). The developed model characterized the in vivo disposition of clopidogrel, its bioactivation into clop-AM in the liver and subsequent platelet aggregation inhibition in the systemic circulation reasonably well. It further allowed the identification of covariates that significantly impact clopidogrel's dose-concentration-response relationship. In particular, CYP2C19 intermediate and poor metabolizers converted 26.2% and 39.5% less clopidogrel to clop-AM, respectively, compared to extensive metabolizers. In addition, CES1 G143E mutation carriers have a reduced CES1 activity (82.9%) compared to wild-type subjects, which results in a significant increase in clop-AM formation. An increase in BMI was found to significantly decrease clopidogrel's bioactivation, whereas increased age was associated with increased platelet reactivity. Our PK/PD model analysis suggests that, in order to optimize clopidogrel dosing on a patient-by-patient basis, all of these factors have to be considered simultaneously, e.g. by using quantitative clinical pharmacology tools.

Vatalanib population pharmacokinetics in patients with myelodysplastic syndrome: CALGB 10105 (Alliance)

AIMS

Vatalanib is an oral anti-angiogenesis agent that inhibits vascular endothelial growth factor receptor tyrosine kinases, which in patients showed auto induction of metabolism and variability in pharmacokinetic (PK) disposition. The objective was to characterize the population PK and time-dependent change in vatalanib clearance and assess exposure-toxicity relationship in patients with myelodysplastic syndrome (MDS).

METHODS

This was an open-label phase II study of vatalanib in MDS patients receiving 750-1250 mg once daily in 28-day cycles. Serial blood samples were obtained and plasma vatalanib concentrations measured by HPLC. Population PK analysis was performed using nonmem 7.2 with FO estimation since FOCE failed. The final model was evaluated using goodness-of-fit plots, bootstrap analysis, and visual predictive check.

RESULTS

Pharmacokinetic data were complete for 137 patients (86 M, 51 F), of median age 70 years (range 20-91). A one-compartment model with lagged first-order absorption and time-dependent change in oral clearance was fitted to the vatalanib plasma concentration versus time data. The population means for pre-induction and post-induction oral clearance were 24.1 l h(-1) (range: 9.6-45.5) and 54.9 l h(-1) (range: 39.8-75.6), respectively. The apparent oral clearance increased 2.3-fold, (range: 1.7-4.1-fold) from first dose to steady state. Our data did not identify a significant relationship of the predefined covariates with vatalanib pharmacokinetics, although power to detect such a relationship was limited.

CONCLUSIONS

Vatalanib pharmacokinetics were highly variable and the extent of auto induction was not determined to correlate with any of the pre-defined covariates.

Opposite malaria and pregnancy effect on oral bioavailability of artesunate - a population pharmacokinetic evaluation

Aim

The aim was to compare the pharmacokinetic properties of artesunate and dihydroartemisinin in the same women: i) pregnant with acute uncomplicated malaria on day 1 and 2, ii) pregnant with convalescent malaria on day 7 and iii) in a healthy state 3 months post-partum on day 1, 2 and 7.

Methods

Non-linear mixed-effects modelling was used to compare plasma concentration–time profiles of artesunate and dihydroartemisinin over 7 days of treatment following oral and intravenous artesunate administration to pregnant women with uncomplicated Plasmodium falciparum malaria during their second or third trimesters of pregnancy. The same women were restudied 3 months after delivery when fully recovered. Non-compartmental results of the same study have been published previously.

Results

Twenty pregnant patients on the Thailand-Myanmar border were studied and 15 volunteered to be restudied 3 months post-partum. Malaria and pregnancy had no effect on the pharmacokinetic properties of artesunate or dihydroartemisinin after intravenous artesunate administration. However, malaria and pregnancy had opposite effects on the absorption of orally administered artesunate. Malaria increased the absolute oral bioavailability of artesunate by 87%, presumably by inhibiting first pass effect, whereas pregnancy decreased oral bioavailability by 23%.

Conclusions

The population pharmacokinetic analysis demonstrated opposite effects of malaria and pregnancy on the bioavailability of orally administered artesunate. Lower drug exposures during the second and third trimesters of pregnancy may contribute to lower cure rates and thus the development of drug resistance. Dose optimization studies are required for artesunate containing artemisinin-based combination therapies (ACTs) in later pregnancy.

Pharmacokinetics of efavirenz and treatment of HIV-1 among pregnant women with and without tuberculosis coinfection

BACKGROUND

Pregnancy and tuberculosis treatment or prophylaxis can affect efavirenz pharmacokinetics, maternal human immunodeficiency virus type 1 (HIV-1) treatment outcomes, and mother-to-child transmission (MTCT) risk.

METHODS

We evaluated a prospective cohort of pregnant, HIV-infected women with and without tuberculosis in Soweto, South Africa. Pharmacokinetic sampling was performed at gestation week 37 and during the postpartum period. Efavirenz trough concentrations (Cmin) were predicted using population pharmacokinetic models. HIV-viral load was measured at delivery for mothers and at 6 weeks of age for infants.

RESULTS

Ninety-seven women participated; 44 had tuberculosis. Median efavirenz Cmin during pregnancy was 1.35 µg/mL (interquartile range [IQR], 0.90-2.07 µg/mL; 27% had an efavirenz Cmin of < 1 µg/mL), compared with a median postpartum value of 2.00 µg/mL (IQR, 1.40-3.59 µg/mL; 13% had an efavirenz Cmin of < 1 µg/mL). A total of 72% of pregnant women with extensive CYP2B6 genotypes had an efavirenz Cmin of <1 µg/mL. Rifampin did not reduce the efavirenz Cmin. Isoniazid (for prophylaxis or treatment), though, reduced the rate of efavirenz clearance. At delivery, median durations of ART were 13 weeks (IQR, 9-18 weeks) and 21 weeks (IQR, 13-64 weeks) for women with and those without tuberculosis, respectively; 55% and 83%, respectively, had a viral load of <20 copies/mL (P = .021). There was 1 case of MTCT.

CONCLUSIONS

Pregnancy increased the risk of low efavirenz concentrations, but MTCT was rare. A detectable HIV-viral load at delivery was more common among pregnant women with tuberculosis, in whom ART was generally initiated later.

Population pharmacokinetics of artesunate and dihydroartemisinin during long-term oral administration of artesunate to patients with metastatic breast cancer

PURPOSE

The purpose of this study were firstly to characterize the population pharmacokinetics of artesunate (ARS) and its active metabolite dihydroartemisinin (DHA) in patients with metastatic breast cancer during long-term (>3 weeks) daily oral ARS administration and secondly to study the relationship between salivary and plasma concentrations of DHA.

METHODS

Drug concentration-time data from 23 patients, receiving oral ARS (100, 150, or 200 mg OD), was analyzed using nonlinear mixed effects modeling. A combined drug-metabolite population pharmacokinetic model was developed to describe the plasma pharmacokinetics of ARS and DHA in plasma. Saliva drug concentrations were incorporated as being directly proportional to plasma concentrations.

RESULTS

A first-order absorption model for ARS linked to a combined two-compartment disposition model for ARS and one-compartment disposition model for DHA provided the best fit to the data. No covariates were identified that could explain between-subject variability. A time-dependent increase in apparent elimination clearance of DHA was observed. Salivary DHA concentrations were proportionally correlated with total DHA plasma concentrations, with an estimated slope factor of 0.116.

CONCLUSIONS

Population pharmacokinetics of ARS and DHA in patients with breast cancer was well described by a combined drug-metabolite model without any covariates and with an increase in apparent elimination clearance of DHA over time. The estimated DHA saliva/plasma ratio was in good agreement with the reported DHA unbound fraction in human plasma. Saliva ARS concentrations correlated poorly with plasma concentrations. This suggests the use of saliva sampling for therapeutic drug monitoring of DHA. However, further studies are warranted to investigate the robustness of this approach.

Reprint of "Pharmacokinetic modelling of the anti-malarial drug artesunate and its active metabolite dihydroartemisinin"

A four compartment mechanistic mathematical model is developed for the pharmacokinetics of the commonly used anti-malarial drug artesunate and its principle metabolite dihydroartemisinin following oral administration of artesunate. The model is structurally unidentifiable unless additional constraints are imposed. Combinations of mechanistically derived constraints are considered to assess their effects on structural identifiability and on model fits. Certain combinations of the constraints give rise to locally or globally identifiable model structures. Initial validation of the model under various combinations of the constraints leading to identifiable model structures was performed against a dataset of artesunate and dihydroartemisinin concentration-time profiles of 19 malaria patients. When all the discussed constraints were imposed on the model, the resulting globally identifiable model structure was found to fit reasonably well to those patients with normal drug absorption profiles. However, there is wide variability in the fitted parameters and further investigation is warranted.

Population pharmacokinetic modelling to assess clinical drug-drug interaction between AZD7325 and midazolam

WHAT IS KNOWN AND OBJECTIVE

AZD7325 is a selective gamma-amino-butyric acid (GABAA )α2, 3 receptor modulator. The aims of this analysis were to develop population pharmacokinetic (PPK) models of AZD7325 and midazolam and to assess the induction effect of AZD7325 on CYP3A4 with midazolam as a substrate.

METHODS

Drug-drug interaction data of AZD7325 and midazolam from 24 healthy subjects were available for model development. PPK models were developed in a sequential manner using NONMEM. Both AZD7325 and midazolam pharmacokinetics were described by two-compartment models, and a transit compartment absorption model and a first-order absorption model were applied for the absorption of AZD7325 and midazolam, respectively. The induction of CYP3A by AZD7325 was described by a transit enzyme model, where the elimination of midazolam was proportionally linked to the enzyme amount. Simulations were performed to predict dosing regimens to account for the induction of CYP3A4.

RESULTS AND DISCUSSION

The population estimates for AZD7325 clearance, intercompartmental clearance, central and peripheral volume were 36, 29·2 L/h, 169 and 392 L, respectively, with interindividual variability (IIV) of 35% and 24% for clearance and central volume. Midazolam clearance, intercompartmental clearance, central and peripheral volume were estimated to be 62·7, 34·7 L/h, 133 and 146 L, respectively, with 43% IIV for clearance. The estimated mean transit time for induction of the CYP3A4 enzyme was 197 h, with 57% IIV.

WHAT IS NEW AND CONCLUSION

The PPK models developed adequately described the clinical observation of AZD7325-mediated CYP3A4 enzyme induction with midazolam as a probe. The model could provide basis for the rational dosing of AZD7325 in clinical practice.

A semiphysiological population pharmacokinetic model for dynamic inhibition of liver and gut wall cytochrome P450 3A by voriconazole

BACKGROUND

Accurate predictions of cytochrome P450 (CYP) 3A-mediated drug-drug interactions (DDIs) account for dynamic changes of CYP3A activity at both major expression sites (liver and gut wall) by considering the full pharmacokinetic profile of the perpetrator and the substrate. Physiological-based in vitro-in vivo extrapolation models have become of increasing interest. However, due to discrepancies between the predicted and observed magnitude of DDIs, the role of models fully based on in vivo data is still essential.

OBJECTIVE

The primary objective of this study was to develop a coupled dynamic model for the interaction of the CYP3A inhibitor voriconazole and the prototypical CYP3A substrate midazolam.

METHODS

Raw concentration data were obtained from a DDI study. Ten subjects were given either no pretreatment (control) or voriconazole twice daily orally. Midazolam was given either intravenously or orally after the last voriconazole dose and during control phases. Data analysis was performed by the population pharmacokinetic approach using non-linear mixed effects modelling (NONMEM 7.2.0). Model evaluation was performed using visual predictive checks and bootstrap analysis.

RESULTS

A semiphysiological model was able to describe the pharmacokinetics of midazolam, its major metabolite and voriconazole simultaneously. By considering the temporal disposition of all three substances in the liver and gut wall, a time-varying CYP3A inhibition process was implemented. Only the incorporation of hypothetical enzyme site compartments resulted in an adequate fit, suggesting a sustained inhibitory effect through accumulation. Novel key features of this analysis are the identification of (1) an apparent sustained inhibitory effect by voriconazole due to a proposed quasi accumulation at the enzyme site, (2) a significantly reduced inhibitory potency of intravenous voriconazole for oral substrates, (3) voriconazole as a likely uridine diphosphate glucuronosyltransferase (UGT) 2B inhibitor and (4) considerable sources of interindividual variability.

CONCLUSION

The proposed semiphysiological modelling approach generated a mechanistic description of the complex DDI occurring at major CYP3A expression sites and thus may serve as a powerful tool to maximise information acquired from clinical DDI studies. The model has been shown to draw precise and accurate predictions. Therefore, simulations based on this kind of models may be used for various clinical scenarios to improve pharmacotherapy.

Clinical population pharmacokinetics and toxicodynamics of linezolid

Thrombocytopenia is a common side effect of linezolid, an oxazolidinone antibiotic often used to treat multidrug-resistant Gram-positive bacterial infections. Various risk factors have been suggested, including linezolid dose and duration of therapy, baseline platelet counts, and renal dysfunction; still, the mechanisms behind this potentially treatment-limiting toxicity are largely unknown. A clinical study was conducted to investigate the relationship between linezolid pharmacokinetics and toxicodynamics and inform strategies to prevent and manage linezolid-associated toxicity. Forty-one patients received 42 separate treatment courses of linezolid (600 mg every 12 h). A new mechanism-based, population pharmacokinetic/toxicodynamic model was developed to describe the time course of plasma linezolid concentrations and platelets. A linezolid concentration of 8.06 mg/liter (101% between-patient variability) inhibited the synthesis of platelet precursor cells by 50%. Simulations predicted treatment durations of 5 and 7 days to carry a substantially lower risk than 10- to 28-day therapy for platelet nadirs of <100 ×10(9)/liter. The risk for toxicity did not differ noticeably between 14 and 28 days of therapy and was significantly higher for patients with lower baseline platelet counts. Due to the increased risk of toxicity after longer durations of linezolid therapy and large between-patient variability, close monitoring of patients for development of toxicity is important. Dose individualization based on plasma linezolid concentration profiles and platelet counts should be considered to minimize linezolid-associated thrombocytopenia. Overall, oxazolidinone therapy over 5 to 7 days even at relatively high doses was predicted to be as safe as 10-day therapy of 600 mg linezolid every 12 h.

Pharmacokinetic modelling of the anti-malarial drug artesunate and its active metabolite dihydroartemisinin

A four compartment mechanistic mathematical model is developed for the pharmacokinetics of the commonly used anti-malarial drug artesunate and its principle metabolite dihydroartemisinin following oral administration of artesunate. The model is structurally unidentifiable unless additional constraints are imposed. Combinations of mechanistically derived constraints are considered to assess their effects on structural identifiability and on model fits. Certain combinations of the constraints give rise to locally or globally identifiable model structures. Initial validation of the model under various combinations of the constraints leading to identifiable model structures was performed against a dataset of artesunate and dihydroartemisinin concentration-time profiles of 19 malaria patients. When all the discussed constraints were imposed on the model, the resulting globally identifiable model structure was found to fit reasonably well to those patients with normal drug absorption profiles. However, there is wide variability in the fitted parameters and further investigation is warranted.

Population pharmacokinetics of fusidic acid: rationale for front-loaded dosing regimens due to autoinhibition of clearance

The objectives of this analysis were to develop a population pharmacokinetic (PK) model to describe the absorption and disposition of fusidic acid after single and multiple doses and to determine the effect of food on the rate and extent of bioavailability. Plasma PK data from three phase 1 studies (n = 75; n = 14 with and without food) in which healthy subjects received sodium fusidate (500 to 2,200 mg) as single or multiple oral doses every 8 h (q8h) or q12h for up to 7 days were modeled using S-ADAPT (MCPEM algorithm). Accumulation of fusidic acid after multiple doses was more than that predicted based on single-dose data. The PK of fusidic acid was best described using a time-dependent mixed-order absorption process, two disposition compartments, and a turnover process to describe the autoinhibition of clearance. The mean total clearance (% coefficient of variation) was 1.28 liters/h (33%) and the maximum extent of autoinhibition was 71.0%, with a 50% inhibitory concentration (IC(50)) of 46.3 mg/liter (36%). Food decreased the extent of bioavailability by 18%. As a result of the autoinhibition of clearance, steady state can be achieved earlier with dosing regimens that contain higher doses (after 8 days for 750 mg q12h and 1 day for 1,500 mg q12h on day 1 followed by 600 mg q12h versus 3 weeks for 500 mg q12h). Given that large initial doses autoinhibit the clearance of fusidic acid, this characteristic provides a basis for the administration of front-loaded dosing regimens of sodium fusidate which would allow for effective concentrations to be achieved early in therapy.

Antihistomonal effects of artemisinin and Artemisia annua extracts in vitro could not be confirmed by in vivo experiments in turkeys and chickens

Five different Artemisia annua-derived materials (i.e. dry leaves, pure artemisinin, and hexane, dichloromethane or methanol extracts of leaves) were screened for their in vitro activities against six clonal cultures of Histomonas meleagridis. Except for the methanol extract, all tested materials displayed in vitro activity against all tested protozoal clones. Neither the dry plant material, extracts nor artemisinin showed any antibacterial activity against the xenic bacteria accompanying the six H. meleagridis clones at concentration levels identical to the antihistomonal setting. The dichloromethane extract of dry leaves (Ext-DCM) (minimal lethal concentration=1.0 mg/ml) and artemisinin (half-maximal inhibitory concentration=1.295 mg/ml) had the most promising antihistomonal properties and were therefore subsequently tested in a standardized experimental infection model in both turkeys and chickens infected with clonal H. meleagridis. There were no differences between treatment groups, where all infected turkeys showed severe clinical histomonosis and demonstrated severe typhlohepatitis typical for histomonosis. Consistent with the infection model used, the infected chickens did not show any adverse clinical signs but contracted severe lesions in their caeca 7 and 10 days post infection (d.p.i.), liver lesions were absent to mild after 7 d.p.i. and progressed to severe lesions at 10 d.p.i.; thus no differences between treatment groups were observed. In conclusion, neither artemisinin nor Ext-DCM was able to prevent experimental histomonosis in turkeys and chickens at the given concentrations, which is contrary to the antihistomonal effect noticed in vitro even though the same clonal culture was used. The results of this study therefore clearly demonstrate the importance of defined in vivo experimentation in order to assess and verify in vitro results.