Xpose--an S-PLUS based population pharmacokinetic/pharmacodynamic model building aid for NONMEM

Semi-physiological Pharmacokinetic Model of Clozapine and Norclozapine in Healthy, Non-smoking Volunteers: The Impact of Race and Genetics

BACKGROUND AND OBJECTIVES: Clozapine is the medication of choice for treatment-resistant schizophrenia. However, it has a complex metabolism and unexplained interindividual variability. The current work aims to develop a pharmacokinetic model of clozapine and norclozapine in non-smokers and assess the impact of demographic and genetic predictors.

METHODS

Healthy volunteers were recruited in a population pharmacokinetic study. Blood samples were collected at 30 min and 1, 2, 3, 5 and 8 h following a single flat dose of clozapine (12.5 mg). The clozapine and norclozapine concentrations were measured via high-performance liquid chromatography-ultraviolet method. A semi-physiological pharmacokinetic model of clozapine and norclozapine was developed using nonlinear mixed-effects modeling. Clinical and genetic predictors were evaluated, including CYP1A2 (rs762551) and ABCB1 (rs2032582), using restriction fragment length polymorphism.

RESULTS

A total of 270 samples were collected from 33 participants. The data were best described using a two-compartment model for clozapine and a two-compartment model for norclozapine with first-order absorption and elimination and pre-systemic metabolism. The estimated (relative standard error) clearance of clozapine and norclozapine were 27 L h-1 (31.5 %) and 19.6 L h-1 (30%), respectively. Clozapine clearance was lower in sub-Saharan Africans (n = 4) and higher in Caucasians (n = 9) than Asians (n = 20). Participants with CYP1A2 (rs762551) (n = 18) and ABCB1 (rs2032582) (n = 12) mutant alleles had lower clozapine clearance in the univariate analysis.

CONCLUSIONS

This is the first study to develop a semi-physiological pharmacokinetic model of clozapine and norclozapine accounting for the pre-systemic metabolism. Asians required lower doses of clozapine as compared with Caucasians, while clozapine pharmacokinetics in sub-Saharan Africans should be further investigated in larger trials.

Tenofovir, emtricitabine, lamivudine and dolutegravir concentrations in plasma and urine following drug intake cessation in a randomized controlled directly observed pharmacokinetic trial to aid point-of-care testing

BACKGROUND

Poor adherence to ART and pre-exposure prophylaxis (PrEP) can impact patient and public health. Point-of-care testing (POCT) may aid monitoring and adherence interventions.

OBJECTIVES

We report the pharmacokinetics of tenofovir [dosed as tenofovir disoproxil (TDF) and tenofovir alafenamide (TAF)], emtricitabine (FTC), lamivudine (3TC) and dolutegravir (DTG) in plasma and urine following drug cessation to evaluate adherence targets in urine for POCT.

METHODS

Subjects were randomized (1:1) to receive DTG/FTC/TAF or DTG/3TC/TDF for 15 days. Plasma and spot urine were collected on Day 15 (0-336 h post final dose). Drug concentrations were quantified using LC-MS, and non-linear mixed-effects models applied to determine drug disposition between matrices and relationship with relevant plasma [dolutegravir protein-adjusted 90% inhibitory concentration (PA-IC90 = 64 ng/mL) and minimum effective concentration (MEC = 324 ng/mL)] and urinary thresholds [tenofovir disoproxil fumarate 1500 ng/mL].

RESULTS

Of 30 individuals enrolled, 29 were included (72% female at birth, 90% Caucasian). Median (range) predicted time to plasma dolutegravir PA-IC90 and MEC were 83.5 (41.0-152) and 49.0 h (23.7-78.9), corresponding to geometric mean (90%) urine concentrations of 5.42 (4.37-6.46) and 27.4 ng/mL (22.1-32.7). Tenofovir in urine reached 1500 ng/mL by 101 h (58.6-205) with an equivalent plasma concentration of 6.20 ng/mL (4.21-8.18).

CONCLUSIONS

These data support use of a urinary tenofovir threshold of <1500 ng/mL (tenofovir disoproxil fumarate-based regimens) as a marker of three or more missed doses for a POCT platform. However, due to low dolutegravir concentrations in urine, POCT would be limited to a readout of recent dolutegravir intake (one missed dose).

Integration of individual preclinical and clinical anti-infective PKPD data to predict clinical study outcomes

The AIDA randomized clinical trial found no significant difference in clinical failure or survival between colistin monotherapy and colistin-meropenem combination therapy in carbapenem-resistant Gram-negative infections. The aim of this reverse translational study was to integrate all individual preclinical and clinical pharmacokinetic-pharmacodynamic (PKPD) data from the AIDA trial in a pharmacometric framework to explore whether individualized predictions of bacterial burden were associated with the trial outcomes. The compiled dataset included for each of the 207 patients was (i) information on the infecting Acinetobacter baumannii isolate (minimum inhibitory concentration, checkerboard assay data, and fitness in a murine model), (ii) colistin plasma concentrations and colistin and meropenem dosing history, and (iii) disease scores and demographics. The individual information was integrated into PKPD models, and the predicted change in bacterial count at 24 h for each patient, as well as patient characteristics, was correlated with clinical outcomes using logistic regression. The in vivo fitness was the most important factor for change in bacterial count. A model-predicted growth at 24 h of ≥2-log10 (164/207) correlated positively with clinical failure (adjusted odds ratio, aOR = 2.01). The aOR for one unit increase of other significant predictors were 1.24 for SOFA score, 1.19 for Charlson comorbidity index, and 1.01 for age. This study exemplifies how preclinical and clinical anti-infective PKPD data can be integrated through pharmacodynamic modeling and identify patient- and pathogen-specific factors related to clinical outcomes - an approach that may improve understanding of study outcomes.

Saliva monitoring of prednisolone in children with first onset steroid-sensitive nephrotic syndrome: Is it possible?

AIMS

Prednisolone is the cornerstone of treatment for idiopathic nephrotic syndrome in children, but is associated with marked side-effects. Therapeutic drug monitoring using saliva would be a patient-friendly option to monitor prednisolone treatment. To assess the feasibility of saliva monitoring, we described the pharmacokinetics (PK) of unbound prednisolone in plasma and saliva of children with first onset steroid-sensitive nephrotic syndrome (SSNS).

METHODS

Children (age 2-16 years) with SSNS participating in a randomized, placebo-controlled trial with levamisole were treated with an 18-week tapering schedule of prednisolone. Five serial samples were collected at 4 (saliva) and 8 weeks (saliva and plasma) after first onset. A nonlinear mixed-effects model was used to estimate the PK parameters of unbound prednisolone and the saliva-to-plasma ratio. Monte Carlo simulations were performed to assess the predictive performance of saliva monitoring.

RESULTS

From 39 children, 109 plasma and 275 saliva samples were available. Estimates (relative squared error) of unbound plasma clearance and volume of distribution were 93 (5%) L h-1 70 kg-1 and 158 (7%) L 70 kg-1, respectively. Typical saliva-to-plasma ratio was 1.30 (8%). Monte Carlo simulations demonstrated that on basis of 4 saliva samples and a single plasma sample unbound plasma area-under-the-concentration-time curve can be predicted within 20% imprecision in 79% of the patients compared to 87% based on 4 plasma samples.

CONCLUSION

Saliva proved to be a reliable and patient-friendly option to determine prednisolone plasma exposure in children with SSNS. This opens opportunities for further PK and pharmacodynamics studies of prednisolone in a variety of paediatric conditions.

Population Pharmacokinetics of Cabotegravir Following Oral Administration and Long-Acting Intramuscular Injection in Real-World People with HIV

Long-acting cabotegravir has been studied mainly in the stringent framework of clinical trials, which does not necessarily reflect the situation of people with HIV (PWH) in routine clinical settings. The present population pharmacokinetic analysis aims to build real-world reference percentile curves of cabotegravir concentrations, accounting for patient-related factors that may affect cabotegravir exposure. The second objective is to simulate whether dosing interval adjustments of cabotegravir could be considered in specific subpopulations. Overall, 238 PWH contributed to 1,038 cabotegravir levels (186 during the initial oral administration phase and 852 after intramuscular injection). Cabotegravir pharmacokinetics was best described using a one-compartment model with distinct first order-absorption for oral and intramuscular administrations, and identical volume and clearance. Our model showed almost 40% faster absorption and 30% higher clearance than previously reported, resulting in a time to steady-state of 8 months and an elimination half-life of 4.6 weeks for long-acting cabotegravir. Sex and body mass index significantly influenced absorption, and bodyweight affected clearance. Model-based simulations showed that cabotegravir trough concentrations in females were 25% lower 4 weeks after the intramuscular loading dose, but 42% higher during the late maintenance phase. Finally, simulations indicated that in females, despite significantly higher cabotegravir concentrations, longer intervals between injections may not consistently ensure levels above the 4-fold protein-adjusted 90% inhibitory target concentration.

Population pharmacokinetic modelling of primaquine exposures in lactating women and breastfed infants

Current guidelines advise against primaquine treatment for breastfeeding mothers to avoid the potential for haemolysis in infants with G6PD deficiency. To predict the haemolytic risk, the amount of drug received from the breast milk and the resulting infant drug exposure need to be characterised. Here, we develop a pharmacokinetic model to describe the drug concentrations in breastfeeding women using venous, capillary, and breast milk data. A mother-to-infant model is developed to mimic the infant feeding pattern and used to predict their drug exposures. Primaquine and carboxyprimaquine exposures in infants are <1% of the exposure in mothers. Therefore, even in infants with the most severe G6PD deficiency variants, it is highly unlikely that standard doses of primaquine (0.25-1 mg base/kg once daily given to the mother for 1-14 days) would cause significant haemolysis. After the neonatal period, primaquine should not be restricted for breastfeeding women (Clinical Trials Registration: NCT01780753).

Pharmacokinetic-pharmacodynamic analysis of drug liking blockade by buprenorphine subcutaneous depot (CAM2038) in participants with opioid use disorder

Buprenorphine is used to treat opioid use disorder (OUD). Weekly and monthly subcutaneous long-acting buprenorphine injections (CAM2038) provide more stable buprenorphine plasma levels and reduce the treatment burden, misuse, and diversion associated with sublingual transmucosal buprenorphine formulations. To characterize the pharmacokinetic/pharmacodynamic (PK/PD) relationship, a maximum inhibition (Imax) model was developed relating CAM2038 buprenorphine plasma concentration to drug liking maximum effect (Emax) visual analog scale (VAS; bipolar) score after intramuscular hydromorphone administration. Data included time-matched observations of buprenorphine plasma concentration and drug liking Emax VAS score after hydromorphone 18 mg administration in 47 non-treatment-seeking adults with moderate to severe OUD in a phase 2 study. Analysis used non-‍linear mixed-effects modeling (NONMEM®). The final Imax model adequately described the PK/PD relationship between buprenorphine plasma concentration and drug liking Emax VAS score. Simulations showed drug liking was effectively blocked at low buprenorphine plasma concentrations (0.4 ng/mL) where the upper 95% confidence interval of the drug liking Emax VAS score was below the pre-defined 11-point complete blockade threshold. The buprenorphine plasma concentration required to achieve 90% of the maximal effect (IC90) of drug liking was 0.675 ng/mL. Interindividual variability in responses to buprenorphine was observed; some participants experienced fluctuating responses, and a few did not achieve drug liking blockade even with higher buprenorphine plasma concentrations. This affirms the need to individualize treatment and titrate doses for optimal treatment outcomes. PK/PD models were also developed for desire to use VAS and Clinical Opiate Withdrawal Scale (COWS) scores, with results aligned to those for drug liking.

Confirmation of the utility of the CA-125 elimination rate (KELIM) as an indicator of the chemosensitivity in advanced-stage ovarian cancer in a "real-life setting"

OBJECTIVE

The modeled CA-125 ELIMination rate constant K (KELIM) has been validated as a marker of response to chemotherapy in >12,000 patients with advanced epithelial ovarian carcinoma (EOC) treated in first-line setting enrolled in >12 clinical trials. Patient KELIM is calculable online https://www.biomarker-kinetics.org/presentation. The objective was to investigate the prognostic value of KELIM in a large real-life national cancer registry with non-selected patients.

METHODS

We investigated 4,025 EOC patients from the Netherlands Cancer Registry treated with neoadjuvant chemotherapy (NACT) ± followed by interval debulking surgery (IDS). Patient KELIM values were calculated in patients with ≥ 3 CA-125 measurements during NACT. KELIM was standardized with a pre-specified cut-off and scored as unfavorable/favorable (<1.0/≥1.0). KELIM's prognostic value regarding radiological response, completeness of IDS, progression-free survival (PFS), and overall survival (OS) was assessed using univariate/multivariate analyses.

RESULTS

The data from 1,582 patients treated with heterogeneous chemotherapy regimens and sequences were assessable. KELIM was prognostic for radiological response and the likelihood of complete IDS after NACT (odds ratio=2.59; 95% confidence interval [CI]=2.04-3.29). Moreover, KELIM was independently associated with PFS (hazard ratio [HR]=0.76; 95% CI=0.66-0.87), and OS (HR=0.79; 95% CI=0.69-0.91). Combining KELIM with the completeness of the IDS resulted in 3 prognostic groups (satisfactory, intermediate, and poor) with significant OS differences, namely a good, intermediate, and poor survival respectively.

CONCLUSION

The value of KELIM, as a pragmatic indicator of response to chemotherapy, was maintained in a large real-life population-based cohort, highlighting its applicability in routine conditions.

A Population Pharmacokinetic Analysis of L-Glutamine Exposure in Patients with Sickle Cell Disease: Evaluation of Dose and Food Effects

BACKGROUND AND OBJECTIVE

L-Glutamine is a treatment for children and adults with sickle cell disease. A comprehensive evaluation of the pharmacokinetics of L-glutamine in sickle cell disease has not been conducted. We aimed to assess the effects of long-term dosing, multiple dose levels, and food intake on L-glutamine exposure in patients with sickle cell disease compared to normal participants.

METHODS

We conducted an open-label dose-ascending trial of L-glutamine in pediatric and adult participants with sickle cell disease (N = 8) and adult healthy volunteers (N = 4), providing a total of 400 plasma L-glutamine concentrations. Each participant received three ascending oral doses (0.1 and 0.3 g/kg twice daily and 0.6 g/kg once daily) over 3 weeks. Plasma L-glutamine concentrations were quantified using ion exchange chromatography. Both a non-compartmental pharmacokinetic analysis and a population pharmacokinetic analysis were performed.

RESULTS

L-glutamine had rapid absorption and elimination, and there was no significant change in the baseline (pre-dose) L-glutamine concentration throughout the study, indicating no drug accumulation. Pharmacokinetics was best described by a one-compartment model with first-order kinetics. The dose-normalized peak concentration decreased with dose escalation, indicating the capacity-limited non-linear pharmacokinetics of oral L-glutamine. A covariate analysis showed that baseline L-glutamine concentrations correlated negatively with glutamine clearance, whereas dose positively correlated with volume of distribution. Food intake did not significantly affect glutamine clearance, indicating that L-glutamine can be taken with or without food.

CONCLUSIONS

We report the first pharmacokinetic study of multiple-dose, long-term oral L-glutamine therapy and the first population pharmacokinetic analysis of L-glutamine for sickle cell disease. These findings may permit optimized dosing of L-glutamine for patients with sickle cell disease to maximize treatment benefits.

CLINICAL TRIAL REGISTRATION

This trial is registered at ClinicalTrials.gov (NCT04684381).

Model-Informed Precision Dosing to Reduce Vincristine-Induced Peripheral Neuropathy in Pediatric Patients: A Pharmacokinetic and Pharmacodynamic Modeling and Simulation Analysis

BACKGROUND

Vincristine-induced peripheral neuropathy (VIPN) is a common adverse effect of vincristine, a drug often used in pediatric oncology. Previous studies demonstrated large inter- and intrapatient variability in vincristine pharmacokinetics (PK). Model-informed precision dosing (MIPD) can be applied to calculate patient exposure and individualize dosing using therapeutic drug monitoring (TDM) measurements. This study set out to investigate the PK/pharmacodynamic (PKPD) relationship of VIPN and determine the utility of MIPD to support clinical decisions regarding dose selection and individualization.

METHODS

Data from 35 pediatric patients were utilized to quantify the relationship between vincristine dose, exposure and the development of VIPN. Measurements of vincristine exposure and VIPN (Common Terminology Criteria for Adverse Events [CTCAE]) were available at baseline and for each subsequent dosing occasions (1-5). A PK and PKPD analysis was performed to assess the inter- and intraindividual variability in vincristine exposure and VIPN over time. In silico trials were performed to portray the utility of vincristine MIPD in pediatric subpopulations with a certain age, weight and cytochrome P450 (CYP) 3A5 genotype distribution.

RESULTS

A two-compartmental model with linear PK provided a good description of the vincristine exposure data. Clearance and distribution parameters were related to bodyweight through allometric scaling. A proportional odds model with Markovian elements described the incidence of Grades 0, 1 and ≥ 2 VIPN overdosing occasions. Vincristine area under the curve (AUC) was the most significant exposure metric related to the development of VIPN, where an AUC of 50 ng⋅h/mL was estimated to be related to an average VIPN probability of 40% over five dosing occasions. The incidence of Grade ≥ 2 VIPN reduced from 62.1 to 53.9% for MIPD-based dosing compared with body surface area (BSA)-based dosing in patients. Dose decreases occurred in 81.4% of patients with MIPD (vs. 86.4% for standard dosing) and dose increments were performed in 33.4% of patients (no dose increments allowed for standard dosing).

CONCLUSIONS

The PK and PKPD analysis supports the use of MIPD to guide clinical dose decisions and reduce the incidence of VIPN. The current work can be used to support decisions with respect to dose selection and dose individualization in children receiving vincristine.

Tumor growth inhibition modeling in patients with second line biliary tract cancer and first line non-small cell lung cancer based on bintrafusp alfa trials

This analysis aimed to quantify tumor dynamics in patients receiving either bintrafusp alfa (BA) or pembrolizumab, by population pharmacokinetic (PK)-pharmacodynamic modeling, and investigate clinical and molecular covariates describing the variability in tumor dynamics by pharmacometric and machine-learning (ML) approaches. Data originated from two clinical trials in patients with biliary tract cancer (BTC; NCT03833661) receiving BA and non-small cell lung cancer (NSCLC; NCT03631706) receiving BA or pembrolizumab. Individual drug exposure was estimated from previously developed population PK models. Population tumor dynamics models were developed for each drug-indication combination, and covariate evaluations performed using nonlinear mixed-effects modeling (NLME) and ML (elastic net and random forest models) approaches. The three tumor dynamics' model structures all included linear tumor growth components and exponential tumor shrinkage. The final BTC model included the effect of drug exposure (area under the curve) and several covariates (demographics, disease-related, and genetic mutations). Drug exposure was not significant in either of the NSCLC models, which included two, disease-related, covariates in the BA arm, and none in the pembrolizumab arm. The covariates identified by univariable NLME and ML highly overlapped in BTC but showed less agreement in NSCLC analyses. Hyperprogression could be identified by higher tumor growth and lower tumor kill rates and could not be related to BA exposure. Tumor size over time was quantitatively characterized in two tumor types and under two treatments. Factors potentially related to tumor dynamics were assessed using NLME and ML approaches; however, their net impact on tumor size was considered as not clinically relevant.

Population Pharmacokinetics of Digoxin in Nonagenarian Patients: Optimization of the Dosing Regimen

OBJECTIVE

The aim of this study was to develop a population pharmacokinetic model of digoxin in patients over 90 years old and to propose an equation for adjusting digoxin dose in this population.

METHODS

We included 326 nonagenarian patients admitted to Severo Ochoa University Hospital (Spain) who received digoxin and were under therapeutic drug monitoring. All data were retrospectively collected, and population modeling was performed with non-linear mixed-effect modeling software (NONMEM®). One- and two-compartment models were tested to calculate digoxin clearance (Cl), volume of distribution (Vd), absorption rate constant (Ka), and bioavailability (bioavailable fraction, F). The covariates were evaluated by stepwise covariate model building, and the final model was internally validated by bootstrap analysis with 1000 resamples. External validation was performed with another population of 95 patients with the same characteristics as the modeling group.

RESULTS

The population was 26% males, with a mean age of 93.2 years (90-103 years), mean creatinine 1.11 mg/dL (0.42-3.81 mg/dL), and mean total body weight 61.2 kg (40-100 kg). The pharmacokinetics of digoxin were best described by a one-compartment model (ADVAN2 TRANS2), with first-order conditional estimation with interaction. The covariates with influence on our model were creatinine clearance based on the Cockcroft-Gault equation (CG), serum potassium (K), co-administration of loop diuretics, and sex: Cl/F = 4.55 · (CG/36.4)0.468 · 0.83LD · 1.21SEX; Vd/F = 355 · (K/4.3)-0.849; Ka = 1.22 h-1 [where LD indicates loop diuretics (1 for administered, 0 for otherwise) and SEX indicates patient sex (1 for male, 0 for female)]. Based on our results, we proposed an equation to adjust the digoxin dosing regimen in nonagenarian patients: dose (mg) = 0.144 · (CG/36.4)0.468 · 0.83LD · 1.21SEX.

CONCLUSIONS

The greatest influence on digoxin clearance came from renal function calculated by the Cockcroft-Gault equation. Vd was decreased by K. The model developed showed a precise predictive performance to be applied for therapeutic drug monitoring.

Combining pharmacometric models with predictive and prognostic biomarkers for precision therapy in Crohn's disease: A case study of brazikumab

Pharmacometric models were used to investigate the utility of biomarkers in predicting the efficacy (Crohn's Disease Activity Index [CDAI]) of brazikumab and provide a data-driven framework for precision therapy for Crohn's disease (CD). In a phase IIa trial in patients with moderate to severe CD, treatment with brazikumab, an anti-interleukin 23 monoclonal antibody, was associated with clinical improvement. Brazikumab treatment effect was determined to be dependent on the baseline IL-22 (BIL22) or baseline C-reactive protein (BCRP; predictive biomarkers), and placebo effect was found to be correlated with the baseline CDAI (a prognostic biomarker). A maximal total inhibition on CDAI input function of 50.6% and 42.4% was predicted for patients with extremely high BIL22 or BCRP, compared to a maximal total inhibition of 20.9% and 17.8% for patients with extremely low BIL22 or BCRP, respectively, which were mainly due to the placebo effect. We demonstrated that model-derived baseline biomarker levels that achieve 50% of maximum unbound systemic concentration of 22.8 pg/mL and 8.03 mg/L for BIL22 and BCRP as the cutoffs to select subpopulations can effectively identify high-response subgroup patients with improved separation of responders when compared to using the median values as the cutoff. This work exemplifies the utility of pharmacometrics to quantify biomarker-driven responses in biologic therapies and distinguish between predictive and prognostic biomarkers, complementing clinical efforts of identifying subpopulations with higher likelihood of response to brazikumab.

Evaluation of the effect of erenumab on migraine-specific questionnaire in patients with chronic and episodic migraine

Erenumab is a fully human anti-canonical calcitonin gene-related peptide receptor monoclonal antibody approved for migraine prevention. The Migraine-Specific Quality-of-Life Questionnaire (MSQ) is a 14-item patient-reported outcome instrument that measures the impact of migraine on health-related quality of life. Erenumab data from four phase II/III clinical trials were used to develop an item response theory (IRT) model within a nonlinear mixed effects framework, (i) evaluate the MSQ item information with respect to patient disability, (ii) characterize the longitudinal progression of the MSQ, and (iii) quantify the effect of erenumab on the MSQ in patients with migraine. The majority (80%) of information was found to be contained in 9 out of 14 items, extending the current knowledge on the reliability of the MSQ as a psychometric tool. Simulations across three MSQ domains show significant improvement from baseline, exceeding minimally important differences. Overall, the IRT model platform developed herein allows for systematic quantification of the effect of erenumab on the MSQ in patients with migraine.

Pharmacokinetic analysis of acute and dietary exposure to piperonyl butoxide in the mouse

Piperonyl butoxide (PBO) is a popular insecticide synergist present in thousands of commercial, agricultural, and household products. PBO inhibits cytochrome P450 activity, impairing the ability of insects to detoxify insecticides. PBO was recently discovered to also inhibit Sonic hedgehog signaling, a pathway required for embryonic development, and rodent studies have demonstrated the potential for in utero PBO exposure to cause structural malformations of the brain, face, and limbs, or more subtle neurodevelopmental abnormalities. The current understanding of the pharmacokinetics of PBO in mice is limited, particularly with respect to dosing paradigms associated with developmental toxicity. To establish a pharmacokinetic (PK) model for oral exposure, PBO was administered to female C57BL/6J mice acutely by oral gavage (22-1800 mg/kg) or via diet (0.09 % PBO in chow). Serum and adipose samples were collected, and PBO concentrations were determined by HPLC-MS/MS. The serum concentrations of PBO were best fit by a linear one-compartment model. PBO concentrations in visceral adipose tissue greatly exceeded those in serum. PBO concentrations in both serum and adipose tissue decreased quickly after cessation of dietary exposure. The elimination half-life of PBO in the mouse after gavage dosing was 6.5 h (90 % CI 4.7-9.5 h), and systemic oral clearance was 83.3 ± 20.5 mL/h. The bioavailability of PBO in chow was 41 % that of PBO delivered in olive oil by gavage. Establishment of this PK model provides a foundation for relating PBO concentrations that cause developmental toxicity in the rodent models to Sonic hedgehog signaling pathway inhibition.

Population Pharmacokinetic Analysis Supports Initiation Treatment and Bridging from Sublingual Buprenorphine to Subcutaneous Administration of a Buprenorphine Depot (CAM2038) in the Treatment of Opioid Use Disorder

BACKGROUND AND OBJECTIVE

In treating opioid use disorder (OUD), subcutaneous (SC) extended-release buprenorphine (BPN) depots, e.g., CAM2038, have been shown to provide smaller and less frequent fluctuations in BPN plasma concentrations and pharmacodynamic responses, improve outcomes, reduce treatment burden, and lower risks of misuse and diversion compared to daily sublingual (SL) BPN. This analysis characterized the pharmacokinetics (PK) of BPN following intravenous and SL administration, and administration of SC CAM2038 weekly and monthly.

METHODS

Pharmacokinetic data from two Phase 1 and two Phase 2 trials in healthy participants and participants with OUD, respectively, were used to develop a population PK model using non-linear mixed effects modelling. The analysis included data from 252 participants and 10,658 BPN observations.

RESULTS

The disposition of BPN was best described by a three-compartment model with first-order elimination, and absorption of SL BPN and SC CAM2038 weekly and monthly by dual parallel absorption pathways. Model diagnostics indicated good predictive performance of BPN concentrations. Buprenorphine plasma concentration-time profiles were simulated for treatment initiation, switching from SL BPN to CAM2038 weekly and monthly, and tapering after interrupting treatment with CAM2038. Simulations predicted CAM2038 weekly and monthly doses that provided BPN plasma maximum concentration (Cmax) and trough concentration (Ctrough) values at steady state within those observed following SL BPN administration.

CONCLUSIONS

This population PK model supports the use of CAM2038 doses as individualized treatment for OUD across different treatment stages, including initiation, switching from SL BPN according to established dose conversion schedules, and tapering.

TRIAL REGISTRATIONS

ISRCTN41550730 (05/19/2014), ISRCTN24987553 (07/29/2014), NCT02611752 (11/23/2015), NCT02710526 (03/16/2016).

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.

Evaluation of postprandial total triglycerides within the TIGG model for characterizing postprandial response of glucose, insulin, and GLP-1

The TIGG model is the first model to integrate glucose and insulin regulation, incretin effect, and triglyceride (TG) response in the lipoprotein subclasses of chylomicrons and VLDL-V6. This model described the response following a high-fat meal in individuals who are lean, obese, and very obese and provided insights into the possible regulation of glucose homeostasis in the extended period following a meal. Often, total TGs are analyzed within clinical studies, instead of lipoprotein subclasses. We extended the existing TIGG model to capture the observed total TGs and determined if this model could be used to predict the postprandial TG response of chylomicron and VLDL-V6 when only total TGs are available. To assess if the lipoprotein distinction was important for the model, a second model (tTIGG) was developed using only the postprandial response in total TGs, instead of postprandial TG response in chylomicrons and VLDL-V6. The two models were compared on their predictability to characterize the postprandial response of glucose, insulin, and active GLP-1. Both models were able to characterize the postprandial TG response in individuals who are lean, obese, or very obese following a high-fat meal. The extended TIGG model resulted in a better model fit of the glucose data compared to the tTIGG model, indicating that chylomicron and VLDL-V6 provided additional information compared to total TGs. Furthermore, the expanded TIGG model was able to predict the postprandial TG response of chylomicrons and VLDL-V6 using the total TGs and could therefore be used in studies where only total TGs were collected.

Population pharmacokinetic analysis of lopinavir in HIV negative individuals exposed to SARS-CoV-2: a COPEP (COronavirus Post-Exposure Prophylaxis) sub-study

BACKGROUND

Lopinavir/ritonavir (LPV/r) is a drug traditionally used for the treatment of HIV that has been repurposed as a potential post-exposure prophylaxis agent against COVID-19 in the COronavirus Post-Exposure Prophylaxis (COPEP) study. The present analysis aims to evaluate LPV levels in individuals exposed to SARS-CoV-2 versus people living with HIV (PLWH) by developing a population pharmacokinetic (popPK) model, while characterizing external and patient-related factors that might affect LPV exposure along with dose-response association.

METHODS

We built a popPK model on 105 LPV concentrations measured in 105 HIV-negative COPEP individuals exposed to SARS-CoV-2, complemented with 170 LPV concentrations from 119 PLWH followed in our routine therapeutic drug-monitoring programme. Published LPV popPK models developed in PLWH and in COVID-19 patients were retrieved and validated in our study population by mean prediction error (MPE) and root mean square error (RMSE). The association between LPV model-predicted residual concentrations (Cmin) and the appearance of the COVID-19 infection in the COPEP participants was investigated.

RESULTS

A one-compartment model with linear absorption and elimination best described LPV concentrations in both our analysis and in the majority of the identified studies. Globally, similar PK parameters were found in all PK models, and provided close MPEs (from -19.4% to 8.0%, with a RMSE of 3.4% to 49.5%). No statistically significant association between Cmin and the occurrence of a COVID-19 infection could be detected.

CONCLUSION

Our analysis indicated that LPV circulating concentrations were similar between COPEP participants and PLWH, and that published popPK models described our data in a comparable way.

Population pharmacokinetic modelling of fibrinogen in patients with congenital or acquired-chronic or acute-hypofibrinogenaemia

AIMS

Fibrinogen is the key substrate for coagulation. Fibrinogen pharmacokinetics (PK) after single doses of fibrinogen concentrate (FC), using modelling approaches, has only been evaluated in congenital afibrinogenaemic patients. The aims of this study are to characterize the fibrinogen PK in patients with acquired-chronic (cirrhosis) or acute-hypofibrinogenaemia (critical haemorrhage), showing endogenous production. Influencing factors of differences on the fibrinogen PK between subpopulations will be identified.

METHODS

A total of 428 time-concentration values from 132 patients were recorded. Eighty-two out of 428 values were from 41 cirrhotic patients administered with placebo and 90 out of 428 were from 45 cirrhotic patients that were given FC, 161 out of 428 values were from 14 afibrinogenaemic patients and 95 out of 428 values were from 32 severe acute trauma haemorrhagic patients. A turnover model that accounted for endogenous production and exogenous dose was fitted using NONMEM74. The production rate (Ksyn), distribution volume (V), plasma clearance (CL) and concentration yielding to 50% of maximal fibrinogen production (EC50) were estimated.

RESULTS

Fibrinogen disposition was described by a one-compartment model with CL and V values of 0.0456 L·h-1 and 4.34 L·70 kg-1 , respectively. Body weight was statistically significant in V. Three different Ksyn values were identified that increased from 0.00439 g·h-1 (afibrinogenaemia), to 0.0768 g·h-1 (cirrhotics) and 0.1160 g·h-1 (acute severe trauma). EC50 was of 0.460 g·L-1 .

CONCLUSIONS

This model will be key as a support tool for dose calculation to achieve specified target fibrinogen concentrations, in each of the studied populations.

Prediction of relative change in free nerve growth factor following subcutaneous administration of tanezumab, a novel monoclonal antibody to nerve growth factor

Tanezumab is a monoclonal antibody against nerve growth factor (NGF). We investigated tanezumab pharmacokinetic (PK)-NGF relationships and predicted the extent of systemic free NGF suppression with target-mediated drug disposition (TMDD) modeling using data from three pivotal phase III interventional studies (NCT02697773, NCT02709486, and NCT02528188) in patients with osteoarthritis. Patients received tanezumab 2.5 mg or 5 mg every 8 weeks (q8w) subcutaneously. A TMDD model using a previously established population PK model was used to describe plasma tanezumab and serum total NGF concentration data, and simulations were performed to predict "unobserved" free NGF versus time profiles and dose-response relationships for free NGF. A total of 2992 patients had available data for plasma tanezumab or serum total NGF concentrations and were included in the analysis; 706 of these had data for both tanezumab and total NGF concentrations. The model generally performed well to predict observed total NGF concentrations up to ~24 weeks after each dose. Simulations suggested free NGF concentration would be suppressed by ~75% (median) near the peak of tanezumab concentration and by less than 5% (median) around the trough tanezumab concentration with a tanezumab 2.5 mg q8w regimen. Free NGF concentration was predicted to return to baseline level at ~8 weeks (95% prediction interval: 5-16 weeks) after the last tanezumab dose. This model adequately described plasma tanezumab and serum total NGF concentrations following s.c. administration of tanezumab 2.5 or 5 mg q8w, allowed prediction of relative change in systemic free NGF following s.c. administration of tanezumab.

Exposure-response analysis of alemtuzumab in pediatric allogeneic HSCT for nonmalignant diseases: the ARTIC study

Alemtuzumab (anti-CD52 antibody) is frequently prescribed to children with nonmalignant diseases undergoing allogeneic hematopoietic stem cell transplantation (HSCT) to prevent graft failure (GF) and acute graft-versus-host disease (aGVHD). The aim of this multicenter study was the characterization of alemtuzumab population pharmacokinetics to perform a novel model-based exposure-response analysis in 53 children with nonmalignant immunological or hematological disease and a median age of 4.4 years (interquartile range [IQR], 0.8-8.7). The median cumulative alemtuzumab dose was 0.6 mg/kg (IQR, 0.6-1) administered over 2 to 7 days. A 2-compartment population pharmacokinetics model with parallel linear and nonlinear elimination including allometrically scaled bodyweight (median, 17.50 kg; IQR, 8.76-33.00) and lymphocyte count at baseline (mean, 2.24 × 109/L; standard deviation ± 1.87) as significant pharmacokinetic predictors was developed using nonlinear mixed effects modeling. Based on the model-estimated median concentration at day of HSCT (0.77 μg/mL; IQR, 0.33-1.82), patients were grouped into a low- (≤0.77 μg/mL) or high- (>0.77 μg/mL) exposure groups. High alemtuzumab exposure at day of HSCT correlated with delayed CD4+ and CD8+ T-cell reconstitution (P value < .0001) and increased risk of GF (P value = .043). In contrast, alemtuzumab exposure did not significantly influence the incidence of aGVHD grade ≥2, mortality, chimerism at 1 year, viral reactivations, and autoimmunity at a median follow-up of 3.3 years (IQR, 2.5-8.0). In conclusion, this novel population pharmacokinetics model is suitable for individualized intravenous precision dosing to predict alemtuzumab exposure in pediatric allogeneic HSCT for nonmalignant diseases, aiming at the achievement of early T-cell reconstitution and prevention of GF in future prospective studies.

Model-Informed Precision Dosing Guidance of Ethosuximide Developed from a Randomized Controlled Clinical Trial of Childhood Absence Epilepsy

Ethosuximide was identified as the optimal option for new-onset childhood absence epilepsy (CAE) in a randomized, two-phase dose escalation comparative effectiveness trial of ethosuximide, lamotrigine, and valproic acid. However, 47% of ethosuximide initial monotherapy participants experienced short-term treatment failure. This study aimed to characterize the initial monotherapy ethosuximide exposure-response relationship and to propose model-informed precision dosing guidance. Dose titration occurred over a 16-20-week period until patients experienced seizure freedom or intolerable side effects. Subjects with initial monotherapy failure were randomized to one of the other two medications and dose escalation was repeated. A population pharmacokinetic model was created using plasma concentration data (n = 1,320), collected at 4-week intervals from 211 unique participants during both the initial and second monotherapy phases. A logistic regression analysis was performed on the initial monotherapy cohort (n = 103) with complete exposure-response data. Eighty-four participants achieved seizure freedom with a wide range of ethosuximide area under the curves (AUC) ranging from 420 to 2,420 μg·h/mL. AUC exposure estimates for achieving a 50% and 75% probability of seizure freedom were 1,027 and 1,489 μg·h/mL, respectively, whereas the corresponding cumulative frequency of intolerable adverse events was 11% and 16%. Monte Carlo Simulation indicated a daily dose of 40 and 55 mg/kg to achieve 50% and 75% probability of seizure freedom in the overall population, respectively. We identified the need for adjusted mg/kg dosing in different body weight cohorts. This ethosuximide proposed model-informed precision dosing guidance to achieve seizure freedom carries promise to optimize initial monotherapy success for patients with CAE.

Model Based Estimation of Posaconazole Tablet and Suspension Bioavailability in Hospitalized Children Using Real-World Therapeutic Drug Monitoring Data in Patients Receiving Intravenous and Oral Dosing

Invasive fungal infections are a major cause of morbidity and mortality for immunocompromised patients. Posaconazole is approved for treatment and prophylaxis of invasive fungal infection in adult patients, with intravenous, oral suspension, and gastroresistant/delayed-released tablet formulations available. In Europe, until very recently, posaconazole was used off-label in children, although a new delayed-release suspension approved for pediatric use is expected to become available soon. A population pharmacokinetic model was developed which uses posaconazole therapeutic drug monitoring data following intravenous and oral dosing in hospitalized children, thus enabling estimation of pediatric suspension and tablet oral bioavailability. In total, 297 therapeutic drug monitoring plasma levels from 104 children were included in this analysis. The final model was a one-compartment model with first-order absorption and nonlinear elimination. Allometric scaling on clearance and volume of distribution was included a priori. Tablet bioavailability was estimated to be 66%. Suspension bioavailability was estimated to decrease with increasing doses, ranging from 3.8% to 32.2% in this study population. Additionally, concomitant use of proton pump-inhibitors was detected as a significant covariate, reducing suspension bioavailability by 41.0%. This is the first population pharmacokinetic study to model posaconazole data from hospitalized children following intravenous, tablet, and suspension dosing simultaneously. The incorporation of saturable posaconazole clearance into the model has been key to the credible joint estimation of tablet and suspension bioavailability. To aid rational posaconazole dosing in children, this model was used alongside published pharmacodynamic targets to predict the probability of target attainment using typical pediatric dosing regimen.

Population Pharmacokinetic Modeling of Free Plasma and Free Brain Concentrations of Ceftaroline in Healthy and Methicillin-Resistant Staphylococcus aureus-Infected Wistar Rats

A population pharmacokinetic model was developed to describe alterations in ceftaroline brain disposition caused by meningitis in healthy and methicillin-resistant Staphylococcus aureus (MRSA)-infected rats. Blood and brain microdialysate samples were obtained after a single bolus dose of ceftaroline fosamil (20 mg/kg) administered intravenously. Plasma data were modeled as one compartment, and brain data were added to the model as a second compartment, with bidirectional drug transport between plasma and brain (Qin and Qout). The cardiac output (CO) of the animals showed a significant correlation with the relative recovery (RR) of plasma microdialysis probes, with animals with greater CO presenting lower RR values. The Qin was approximately 60% higher in infected animals, leading to greater brain exposure to ceftaroline. Ceftaroline brain penetration was influenced by MRSA infection, increasing from 17% (Qin/Qout) in healthy animals to 27% in infected animals. Simulations of a 2-h intravenous infusion of 50 mg/kg every 8 h achieved >90% probability of target attainment (PTA) in plasma and brain for the modal MRSA MIC (0.25 mg/L), suggesting that the drug should be considered an option for treating central nervous system infections.

Plan E, Weber B. Bounded integer model-based analysis of psoriasis area and severity index in patients with moderate-to-severe plaque psoriasis receiving BI 730357

BI 730357 is investigated as an oral treatment of plaque psoriasis. We analyzed the impact of three dosage regimens on the Psoriasis Area and Severity Index (PASI) response with modeling based on phase I and II data from 109 healthy subjects and 274 patients with moderate-to-severe plaque psoriasis. The pharmacokinetics (PK) was characterized by a two-compartment model with dual absorption paths and a first-order elimination. Higher baseline C-reactive protein was associated with lower clearance and patients generally had lower clearance compared with healthy subjects. A bounded integer PK/pharmacodynamic model characterized the effect on the observed PASI. The maximum drug effect was largest for patients with no prior biologic use, smaller for patients with prior use of non-interleukin-17 inhibitors, and smallest for patients with prior interleukin-17 inhibitor use. The models allowed robust simulation of large patient populations, predicting a plateau in PASI outcomes for BI 730357 exposure above 2000 nmol/L.

Comparison of Three Renal Function Formulas for Ganciclovir/Valganciclovir Dose Individualization in CMV-Infected Solid Organ Transplantation Patients Using a Population Approach

BACKGROUND AND OBJECTIVE

The gold standard treatment of established cytomegalovirus infection or prevention in solid organ transplantation is the intravenous administration of ganciclovir (GCV) or oral administration of valganciclovir (VGCV), both adjusted to the renal function. In both instances, there is a high interindividual pharmacokinetic variability, mainly owing to the wide range of variation of both the renal function and body weight. Therefore, accurate estimation of the renal function is crucial for GCV/VGCV dose optimization. This study aimed to compare three different formulas for estimating the renal function in solid organ transplantation patients with cytomegalovirus infection, for individualizing antiviral therapy with GCV/VGCV, using a population approach.

METHODS

A population pharmacokinetic analysis was performed using NONMEM 7.4. A total of 650 plasma concentrations obtained after intravenous GCV and oral VGCV administrations were analyzed, from intensive and sparse sampling designs. Three different population pharmacokinetic models were built with the renal function given by Cockcroft-Gault, Modification of Diet in Renal Disease, or Chronic Kidney Disease EPIdemiology Collaboration (CKD-EPI) formulas. Pharmacokinetic parameters were allometrically scaled to body weight.

RESULTS

The CKD-EPI formula was identified as the best predictor of between-patient variability in GCV clearance. Internal and external validation techniques showed that the CKD-EPI model had better stability and performed better compared with the others.

CONCLUSIONS

The model based on the more accurate estimation of the renal function with the CKD-EPI formula and body weight as a size metric most used in the clinical practice can refine initial dose recommendations and contribute to GCV and VGCV dose individualization when required in the prevention or treatment of cytomegalovirus infection in solid organ transplantation patients.

Population plasma and urine pharmacokinetics and the probability of target attainment of fosfomycin in healthy male volunteers

PURPOSE

A population pharmacokinetic model of fosfomycin was developed in healthy volunteers after intravenous administration, and different dosing regimens were evaluated in terms of the probability of target attainment for Escherichia coli using both plasma and urinary pharmacokinetic/pharmacodynamic targets.

METHODS

Eight healthy men received fosfomycin as both intermittent 8 g q8h and continuous infusion 1 g/h with a loading dose of 8 g in a crossover study design. Dense sampling was conducted during both regimens. Population pharmacokinetic modelling was performed using NONMEM. Monte Carlo simulations were conducted to evaluate the Probability of Target Attainment (PTA) of different dosing regimens using bactericidal (AUC_24h/MIC of 83 and 75%T_>MIC) and bacteriostatic (AUC_24h/MIC of 25) plasma targets and bacteriostatic (AUC_24h/MIC of 3994) urine target.

RESULTS

A total of 176 plasma and 86 urine samples were available for PK analysis. A two-compartment model with a urine compartment best described the data. Glomerular filtration rate (GFR) showed a significant correlation with renal clearance and was implemented in the final model. Simulation results show that the dose of 4 g q8h reached 100% of PTA using bactericidal and bacteriostatic targets for MIC up to 16 mg/L.

CONCLUSION

For the clinical breakpoint of 32 mg/L, the standard dosing regimen (4 g q8h) might not be sufficient to reach the bactericidal target. Higher dosing of 8 g q8h as an intermittent infusion or 0.75 g/h as a continuous infusion might be required. Continuous infusion resulted in better attainment of the %T_>MIC target than intermittent infusion.

Benefit From Fractionated Dose-Dense Chemotherapy in Patients With Poor Prognostic Ovarian Cancer: ICON-8 Trial

PURPOSE

An international meta-analysis identified a group of patients with advanced epithelial ovarian cancer (EOC) with a very poor survival because of two unfavorable features: (1) a poor chemosensitivity defined by an unfavorable modeled CA-125 ELIMination rate constant K (KELIM) score <1.0 with the online calculator CA-125-Biomarker Kinetics, and (2) an incomplete debulking surgery. We assumed that patients belonging to this poor prognostic group would benefit from a fractionated densified chemotherapy regimen.

METHODS

The data set of ICON-8 phase III trial (ClinicalTrials.gov identifier: NCT01654146), where patients with EOC were treated with the standard three-weekly, or the weekly dose-dense, carboplatin-paclitaxel regimens and debulking primary surgery (immediate primary surgery [IPS] or delayed primary [or interval] surgery [DPS]), was investigated. The association between treatment arm efficacy, standardized KELIM (scored as favorable ≥1.0, or unfavorable <1.0), and surgery completeness was assessed by univariate/multivariate analyses in IPS and DPS cohorts.

RESULTS

Of 1,566 enrolled patients, KELIM was calculated with the online model in 1,334 with ≥3 CA-125 available values (85%). As previously reported, both KELIM and surgery completeness were complementary prognostic covariates, and could be combined into three prognostic groups with large OS differences: (1) good if favorable KELIM and complete surgery; (2) intermediate if either unfavorable KELIM or incomplete surgery; and (3) poor if unfavorable KELIM and incomplete surgery. Weekly dose-dense chemotherapy was associated with PFS/OS improvement in the poor prognostic group in both the IPS cohort (PFS: hazard ratio [HR], 0.50; 95% CI, 0.31 to 0.79; OS: HR, 0.58; 95% CI, 0.35 to 0.95) and the DPS cohort (PFS: HR, 0.53; 95% CI, 0.37 to 0.76; OS: HR, 0.57; 95% CI, 0.39 to 0.82).

CONCLUSION

Fractionated dose-dense chemotherapy might be beneficial for patients belonging to the poor prognostic group characterized by lower tumor chemosensitivity assessed with the online calculator CA-125-Biomarker Kinetics and incomplete debulking surgery. Further investigation in the future SALVOVAR trial is warranted.

Population Pharmacokinetic Study of Benzylpenicillin in Critically Unwell Adults

Pharmacokinetics are highly variable in critical illness, and suboptimal antibiotic exposure is associated with treatment failure. Benzylpenicillin is a commonly used beta-lactam antibiotic, and pharmacokinetic data of its use in critically ill adults are lacking. We performed a pharmacokinetic study of critically unwell patients receiving benzylpenicillin, using data from the ABDose study. Population pharmacokinetic modelling was undertaken using NONMEM version 7.5, and simulations using the final model were undertaken to optimize the pharmacokinetic profile. We included 77 samples from 12 participants. A two-compartment structural model provided the best fit, with allometric weight scaling for all parameters and a creatinine covariate effect on clearance. Simulations (n = 10,000) demonstrated that 25% of simulated patients receiving 2.4 g 4-hourly failed to achieve a conservative target of 50% of the dosing interval with free drug above the clinical breakpoint MIC (2 mg/L). Simulations demonstrated that target attainment was improved with continuous or extended dosing. To our knowledge, this study represents the first full population PK analysis of benzylpenicillin in critically ill adults.

Joint longitudinal model-based meta-analysis of FEV1 and exacerbation rate in randomized COPD trials

Model-based meta-analysis (MBMA) is an approach that integrates relevant summary level data from heterogeneously designed randomized controlled trials (RCTs). This study not only evaluated the predictability of a published MBMA for forced expiratory volume in one second (FEV₁) and its link to annual exacerbation rate in patients with chronic obstructive pulmonary disease (COPD) but also included data from new RCTs. A comparative effectiveness analysis across all drugs was also performed. Aggregated level data were collected from RCTs published between July 2013 and November 2020 (n = 132 references comprising 156 studies) and combined with data used in the legacy MBMA (published RCTs up to July 2013 - n = 142). The augmented data (n = 298) were used to evaluate the predictive performance of the published MBMA using goodness-of-fit plots for assessment. Furthermore, the model was extended including drugs that were not available before July 2013, estimating a new set of parameters. The legacy MBMA model predicted the post-2013 FEV₁ data well, and new estimated parameters were similar to those of drugs in the same class. However, the exacerbation model overpredicted the post-2013 mean annual exacerbation rate data. Inclusion of year when the study started on the pre-treatment placebo rate improved the model predictive performance perhaps explaining potential improvements in the disease management over time. The addition of new data to the legacy COPD MBMA enabled a more robust model with increased predictability performance for both endpoints FEV₁ and mean annual exacerbation rate.

Ustekinumab Trough Concentrations Are Associated with Biochemical Outcomes in Patients with Crohn's Disease

OBJECTIVE

It is unknown whether ustekinumab (UST) levels can predict clinical outcomes in Crohn's disease (CD) patients. We assessed the exposure-response relationship of UST trough concentrations with biochemical outcomes at week 24 in a prospective, real-world setting.

METHODS

We performed a prospective study in patients with CD starting UST in four academic centres in the Netherlands. All patients received a weight-adjusted intravenous (IV) UST induction dose, followed by one subcutaneous (SC) dose of 90 mg UST at 8 weeks. Maintenance therapy consisted of 90 mg subcutaneous UST every 8 or 12 weeks. Individual UST concentration time course during treatment were estimated using a population pharmacokinetic (PK) model. Quartile analysis and logistic regression were performed to analyse if UST concentrations at week 8 were associated with biochemical remission rates at week 24 (C-reactive protein (CRP) ≤ 5 mg/L and / or faecal calprotectin (FC) ≤ 250 mg/kg).

RESULTS

In total, 124 patients with CD were included. Patients achieving biochemical remission at week 12 and 24 had significantly higher UST levels at week 8 compared to patients without biochemical remission (6.6 µg/mL versus 3.9 µg/mL, P < 0.01 and 6.3 µg/mL versus 3.9 µg/mL, P < 0.01, respectively). In quartile analysis, patients with UST levels in the highest quartile (≥ 6.3 µg/mL at week 8) had higher biochemical remission rates at week 12 and week 24. There was no association between UST levels at and corticosteroid-free clinical remission rates.

CONCLUSION

In this real-world cohort of patients with CD, UST levels in the highest quartile (≥ 6.3 µg/mL) at week 8 were associated with higher biochemical remission rates at week 24.

Predictive Factors of Piperacillin Exposure and the Impact on Target Attainment after Continuous Infusion Administration to Critically Ill Patients

Critically ill patients undergo significant pathophysiological changes that affect antibiotic pharmacokinetics. Piperacillin/tazobactam administered by continuous infusion (CI) improves pharmacokinetic/pharmacodynamic (PK/PD) target attainment. This study aimed to characterize piperacillin PK after CI administration of piperacillin/tazobactam in critically ill adult patients with preserved renal function and to determine the empirical optimal dosing regimen. A total of 218 piperacillin concentrations from 106 patients were simultaneously analyzed through the population PK approach. A two-compartment linear model best described the data. Creatinine clearance (CL_CR) estimated by CKD-EPI was the covariate, the most predictive factor of piperacillin clearance (CL) interindividual variability. The mean (relative standard error) parameter estimates for the final model were: CL: 12.0 L/h (6.03%); central and peripheral compartment distribution volumes: 20.7 L (8.94%) and 62.4 L (50.80%), respectively; intercompartmental clearance: 4.8 L/h (26.4%). For the PK/PD target of 100% fT_>1×MIC, 12 g of piperacillin provide a probability of target attainment > 90% for MIC < 16 mg/L, regardless of CL_CR, but higher doses are needed for MIC = 16 mg/L when CL_CR > 100 mL/min. For 100% fT_>4×MIC, the highest dose (24 g/24 h) was not sufficient to ensure adequate exposure, except for MICs of 1 and 4 mg/L. Our model can be used as a support tool for initial dose guidance and during therapeutic drug monitoring.

Mathematical modeling of the early modeled CA-125 longitudinal kinetics (KELIM-PARP) as a pragmatic indicator of rucaparib efficacy in patients with recurrent ovarian carcinoma in ARIEL2 & STUDY 10

BACKGROUND

PARP inhibitors (PARPi) have revolutionized the management of advanced ovarian carcinoma, and were investigated as forefront treatment in recurrent disease. The objective was to explore if mathematical modeling of the early longitudinal CA-125 kinetics could be used as a pragmatic indicator of the subsequent rucaparib efficacy, like it is for platinum-based chemotherapy.

METHODS

The datasets of ARIEL2 and Study 10 involving recurrent HGOC patients treated with rucaparib were retrospectively investigated. The same strategy as those successfully developed for platinum chemotherapy, based on CA-125 ELIMination rate constant K (KELIM™), was implemented. Individual values of rucaparib-adjusted KELIM (KELIM-PARP) were estimated based on the longitudinal CA-125 kinetics during the first 100 treatment days, and then scored as favorable (KELIM-PARP ≥1.0) or unfavorable (KELIM-PARP <1.0). The prognostic value of KELIM-PARP regarding treatment efficacy (radiological response, and progression-free survival (PFS)) was assessed using univariable/multivariable analyses, with respect to platinum-sensitivity and homologous recombination deficiency (HRD) status.

FINDINGS

The data from 476 patients were assessed. The CA-125 longitudinal kinetics during the first 100-treatment days could be accurately assessed using the KELIM-PARP model. In patients with platinum-sensitive diseases, BRCA mutational status KELIM-PARP score and were associated with subsequent complete/partial radiological responses (KELIM-PARP: odds-ratio = 2.81, 95% CI 1.86-4.52), and PFS (KELIM-PARP: hazard-ratio = 0.67, 95% CI 0.50-0.91). The patients with BRCA-wild type cancer and favorable KELIM-PARP experienced long PFS with rucaparib regardless of HRD. In platinum-resistant disease patients, KELIM-PARP was associated with subsequent radiological response (odds-ratio = 2.80, 95% CI 1.82-4.72).

INTERPRETATION

This proof-of-concept study confirms the early CA-125 longitudinal kinetics during rucaparib in recurrent HGOC patients are assessable by mathematical modeling, to generate individual a KELIM-PARP score associated with the subsequent efficacy. This pragmatic strategy might be useful for selecting the patients for PARPi-based combination regimens, when identifying efficacy biomarker is challenging. Further assessment of this hypothesis is warranted.

FUNDING

The present study was supported by Clovis Oncology with a grant to academic research association.

Single dose tocilizumab for COVID-19 associated cytokine storm syndrome: Less is more

AIMS

We aim to evaluate the clinical pharmacokinetics of a single dose interleukin-6 (IL-6) antibody tocilizumab (TCZ) in methylprednisolone (MP)-treated COVID-19 patients with cytokine storm syndrome (CSS).

METHODS

MP pre-treated patients with COVID-19-associated CSS, defined as at least two elevations of C-reactive protein (CRP) >100 mg/L, ferritin >900 μg/L or D-dimers >1500 μg/L, received intravenous TCZ (8 mg/kg, max. 800 mg) upon clinical deterioration. A nonlinear-mixed effects model was developed based on TCZ serum concentrations and dosing information. Population pharmacokinetic parameters were estimated and concentration-time profiles were plotted against individual predicted values. Fixed dose simulations were subsequently performed based on the final model.

RESULTS

In total 40 patients (mean [SD] age: 62 [12] years, 20% female, body weight: 87 [17] kg) with COVID-19 induced CSS were evaluated on pharmacokinetics and laboratory parameters. A biphasic elimination of TCZ serum concentration was described by a homogeneous population pharmacokinetic model. Serum TCZ concentrations above the 1 μg/L target saturation threshold were covered for 16 days in all evaluated patients treated with a single dose of 8 mg/kg. In a simulation with TCZ 400 mg fixed dose, this condition of full IL-6 receptor occupancy at minimum serum concentration was also met.

CONCLUSIONS

A single dose (8 mg/kg, max. 800 mg) is sufficient to cover a period of 16 days of IL-6-mediated hyperinflammation in COVID-19-induced CSS in MP-treated patients. Based on body weight PK simulations, a fixed-dose tocilizumab of 400 mg should be considered to prevent overtreatment, future drug shortage and unnecessary drug expenditure.

Pharmacometric modeling of drug adverse effects: an application of mixture models in schizophrenia spectrum disorder patients treated with clozapine

Clozapine has superior efficacy to other antipsychotics yet is underutilized due to its adverse effects, such as neutropenia, weight gain, and tachycardia. The current investigation aimed to introduce a pharmacometric approach to simultaneously model drug adverse effects, with examples from schizophrenia spectrum patients receiving clozapine. The adverse drug effects were represented as a function of time by incorporating a mixture model to describe individual susceptibility to the adverse effects. Applications of the proposed method were presented by analyzing retrospective data from patients' medical records in Psychiatric Clinic, Penang General Hospital. Tachycardia, weight gain, and absolute neutrophils count (ANC) decrease were best described by an offset, a piecewise linear, and a transient surge function, respectively. 42.9% of the patients had all the adverse effects, including weight gain (0.01 kg/m² increase every week over a baseline of 24.7 kg/m² until stabilizing at 279 weeks), ANC decrease (20% decrease from 4540 cells/µL week 12-20.8), and tachycardia (14% constant increase over a baseline of 87.9 bpm for a clozapine maintenance dose of 450 mg daily). 32.5% of the patients had only tachycardia, while the remaining 24.6% had none of the adverse effects. A new pharmacometric approach was proposed to describe adverse drug effects with examples of clozapine-induced weight gain, ANC drop, and tachycardia. The current approach described the longitudinal time changes of continuous data while assessing patient susceptibility. Furthermore, the model revealed the possible co-existence of ANC drop and weight gain; thus, neutrophil monitoring might predict future changes in body weight.

An innovative ethosuximide granule formulation designed for pediatric use: Comparative pharmacokinetics, safety, tolerability, and palatability profile versus reference syrup

Ethosuximide, the first-line therapy for childhood absence epilepsy, is currently formulated as a syrup (Zarontin®, Pfizer) with a bitter taste and high sugar content, poorly adapted to children, and a ketogenic diet. The collaborative European FP7 project KIEKIDS aimed at developing an innovative sugar-free, tasteless formulation convenient for pediatric use. This dual Phase-I study evaluated two granule formulations based on lipid multiparticulate (LMP) technology. Two panels of 6 healthy adult volunteers underwent a randomized, placebo-controlled, partly blinded, 3-way cross-over trial, comparing ethosuximide granules A or B with placebo granules and syrup at single 10 mg/kg doses. Corresponding plasma pharmacokinetic profiles of ethosuximide were compared, along with palatability, safety, and tolerability. The LMP granule A proved suboptimal due to bitterness and adherence to beaker walls, while the optimized granule B revealed excellent palatability, similar to placebo granules, and low adherence to glass. The relative bioavailability of granules A versus syrup, based on dose-normalized Cmax and AUC0-∞ was 93.7% [90% CI: 76.3-115.1] and 96.1% [91.0-101.5], respectively. For granules B it was 87.6% [81.6-94.0] and 92.5% [88.5-96.6], respectively, with slightly delayed tmax of 0.75 h [0.5-4.05] compared to syrup 0.5 h [0.3-0.8]. Tolerability visual analog scales revealed a trend for statistically non-significant improvement versus syrup at peak (30 min) for transient dizziness (both granules), fatigue (granules A), and anxiety (granules B). The innovative ethosuximide granule formulation B achieves a suitable profile for pediatric use, being sugar-free, tasteless, bioequivalent, and well-tolerated while enabling precise adjustment to body weight.

Pharmacokinetic analysis of linezolid for multidrug resistant tuberculosis at a tertiary care centre in Mumbai, India

Linezolid is an oxazolidinone used to treat multidrug-resistant tuberculosis (MDR-TB), including in the recently-endorsed shorter 6-month treatment regimens. Due to its narrow therapeutic index, linezolid is often either dose-adjusted or discontinued due to intolerance or toxicity during treatment, and the optimal balance between linezolid efficacy and toxicity remains unclear. India carries a significant burden of MDR-TB cases in the world, but limited information on the pharmacokinetics of linezolid and minimum inhibitory concentration (MIC) distribution is available from Indian MDR-TB patients. We enrolled participants from a tertiary care centre in Mumbai, India, treated for MDR-TB and receiving linezolid daily doses of 600 or 300 mg. Pharmacokinetic visits were scheduled between 1 and 15 months after treatment initiation to undergo intensive or sparse blood sampling. Linezolid concentration versus time data were analysed using non-linear mixed-effects modelling, with simulations to evaluate doses for different scenarios. We enrolled 183 participants (121 females), with a median age of 26 years (interquartile range [IQR] 21-35), weight 55.0 kg (IQR 45.6-65.8), and fat-free mass 38.7 kg (IQR 32.7-46.0). Linezolid pharmacokinetics was best described by a one-compartment model with first-order elimination allometrically scaled by fat-free mass and transit compartment absorption. The typical clearance value was 3.81 L/h. Simulations predicted that treatment with 300 mg daily achieves a high probability of target attainment (PTA) when linezolid MIC was ≤0.25 mg/L (61.5% of participant samples tested), while 600 mg daily would be required if MIC were 0.5 mg/L (29% of samples). While linezolid 300 mg daily is predicted to achieve effective targets for the majority of adults with MDR-TB, it failed to achieve the therapeutic target for 21% participants. A dose of 600 mg had a PTA >90% for all susceptible samples, but with a higher likelihood of exceeding toxicity thresholds (31% vs 9.6%). These data suggest potential benefit to individualized dosing taking host and microbial characteristics into account to improve the likelihood of treatment efficacy while minimizing risk of toxicity from linezolid for the treatment of MDR-TB. Further prospective evaluation in different clinical settings is urgently needed to inform safety and efficacy of these lower doses.

Influence of UGT1A1 and SLC22A6 polymorphisms on the population pharmacokinetics and pharmacodynamics of raltegravir in HIV-infected adults: a NEAT001/ANRS143 sub-study

Using concentration-time data from the NEAT001/ARNS143 study (single sample at week 4 and 24), we determined raltegravir pharmacokinetic parameters using nonlinear mixed effects modelling (NONMEM v.7.3; 602 samples from 349 patients) and investigated the influence of demographics and SNPs (SLC22A6 and UGT1A1) on raltegravir pharmacokinetics and pharmacodynamics. Demographics and SNPs did not influence raltegravir pharmacokinetics and no significant pharmacokinetic/pharmacodynamic relationships were observed. At week 96, UGT1A1*28/*28 was associated with lower virological failure (p = 0.012), even after adjusting for baseline CD4 count (p = 0.048), but not when adjusted for baseline HIV-1 viral load (p = 0.082) or both (p = 0.089). This is the first study to our knowledge to assess the influence of SNPs on raltegravir pharmacodynamics. The lack of a pharmacokinetic/pharmacodynamic relationship is potentially an artefact of raltegravir's characteristic high inter and intra-patient variability and also suggesting single time point sampling schedules are inadequate to thoroughly assess the influence of SNPs on raltegravir pharmacokinetics.

Rituximab pharmacokinetic and pharmacokinetic-pharmacodynamic evaluation based on a study in diffuse large B-cell lymphoma: Influence of tumor size on pharmacokinetic and assessment of pharmacokinetic similarity

Dr. Reddy's Laboratories rituximab (DRL_RI; Dr. Reddy's Laboratories SA, Basel, Switzerland) is under development as a rituximab biosimilar. Study RI-01-002 (Clinical Trials Registry - India/2012/11/003129), comparing DRL_RI to the reference medicinal product (RMP) MabThera® (Roche, Grenzach-Wyhlen, Germany), demonstrated pharmacokinetic (PK) equivalence and showed comparable pharmacodynamic, efficacy, safety, and immunogenicity profiles. We used data from the same study to perform population PK and PK-pharmacodynamic analyses: first exploring possible factors influencing the PK similarity assessment between products and then performing simulations to investigate the impact of tumor size on rituximab PK. Nonlinear mixed-effects models for PK, tumor size, tumor size-PK, and tumor response were developed independently. The final PK model included drug product as a dose-scaling parameter and predicted a 6.75% higher dose reaching the system in RMP-treated patients. However, when tumor size was included in the tumor size-PK model, the drug product effect was no longer observed. The model rather indicated that patients with larger tumor size have higher clearance. Further simulations confirmed that higher baseline tumor size is associated to slightly lower rituximab exposure. Tumor response, described by a continuous-time Markov model, did not differ between drug products. Both had higher effects during the first 20 weeks of treatment. Also, the model described a subpopulation of nonresponders to treatment (42%) with faster transitions to a worse state. The different rituximab exposure initially detected between drug products (6.75%) was shown using PK/PK-pharmacodynamic analysis to be attributed to a tumor size imbalance between treatment groups. PK/PK-pharmacodynamic analyses may contribute to PK similarity assessments.

An adaptive phase II/III safety and efficacy randomized controlled trial of single day or three-day fixed-dose albendazole-ivermectin co-formulation versus albendazole for the treatment of Trichuris trichiura and other STH infections. ALIVE trial protocol

Background: Soil-transmitted helminths (STH) are targeted for control through mass drug-administration campaigns to prevent morbidity affecting at-risk groups in endemic regions. Although broadly successful, the use of albendazole and mebendazole achieved variable progress, with deficiencies against Trichuris trichiura and a predictable low efficacy against Strongyloides stercoralis. Novel drug combinations offer a potential solution, providing they can be delivered safely and maintain efficacy against all STH species. Here we present the protocol of a clinical trial to evaluate a fixed-dose combination (FDC) tablet containing albendazole and ivermectin that will be compared against albendazole against STH . Methods: An adaptive phase II/III randomized controlled trial will be undertaken in STH endemic sites in Ethiopia, Kenya and Mozambique to evaluate an oral FDC of 400 mg albendazole and either 9- or 18 mg ivermectin. FDC will be administered as a single dose or single doses over three-consecutive days and assessed against a single dose of 400 mg albendazole. In the phase II trial, 126 T. trichiura-infected children weighting 15 to 45 kg will be treated in a dose-escalation manner to determine safety objectives. In the phase III trial, 1097 participants aged 5 to 18 years old infected with T. trichiura, hookworm and S. stercoralis will be recruited to determine safety and efficacy. The trial will be open-label with blinded outcome assessors. Cure rate measured 21-days after-treatment in duplicate Kato-Katz is the primary efficacy outcome. Secondary objectives include efficacy evaluation by quantitative polymerase chain reaction (PCR) as an outcome measurement, description of pharmacokinetic parameters, palatability and acceptability evaluations, and monitoring of anthelmintic resistance. Conclusions: This trial with registrational goals seeks to evaluate an innovative fixed-dose combination of albendazole and ivermectin co-formulated tablets, with the goal of providing an anthelmintic regimen with improved efficacy and spectrum of coverage against STH. ClinicalTrials.gov registration: NCT05124691 (18/11/2021).

Identification of Patients With Ovarian Cancer Experiencing the Highest Benefit From Bevacizumab in the First-Line Setting on the Basis of Their Tumor-Intrinsic Chemosensitivity (KELIM): The GOG-0218 Validation Study

PURPOSE

In patients with high-grade ovarian cancer, predictors of bevacizumab efficacy in first-line setting are needed. In the ICON-7 trial, a poor tumor intrinsic chemosensitivity (defined by unfavorable modeled cancer antigen-125 [CA-125] ELIMination rate constant K [KELIM] score) was a predictive biomarker. Only the patients with high-risk disease (suboptimally resected stage III, or stage IV) exhibiting unfavorable KELIM score < 1.0 had overall survival (OS) benefit from bevacizumab (median: 29.7 v 20.6 months; hazard ratio [HR], 0.78). An external validation study in the GOG-0218 trial was performed.

METHODS

In GOG-0218, 1,873 patients were treated with carboplatin-paclitaxel ± concurrent-maintenance bevacizumab/placebo. Patient KELIM values were calculated with CA-125 kinetics during the first 100 chemotherapy days by the Lyon University team. The association between KELIM score (favorable ≥ 1.0, or unfavorable < 1.0) and bevacizumab benefit for progression-free survival (PFS)/OS was independently assessed by NGR-GOG using univariate/multivariate analyses.

RESULTS

KELIM was assessable in 1,662 patients with ≥ 3 CA-125 available values. An unfavorable KELIM score was associated with bevacizumab benefit compared with placebo (PFS: HR, 0.70; 95% CI, 0.59 to 0.82; OS: HR, 0.87; 95% CI, 0.73 to 1.03), whereas a favorable KELIM was not (PFS: HR, 0.96; 95% CI, 0.79 to 1.17; OS: HR, 1.11; 95% CI, 0.89 to 1.39). The highest benefit was observed in patients with a high-risk disease exhibiting unfavorable KELIM, for PFS (median: 9.1 v 5.6 months; HR, 0.64; 95% CI, 0.53 to 0.78), and for OS (median: 35.1 v 29.1 months; HR, 0.79; 95% CI, 0.65 to 0.97).

CONCLUSION

This GOG-0218 trial investigation validates ICON-7 findings about the association between poor tumor chemosensitivity and benefit from concurrent-maintenance bevacizumab, suggesting that bevacizumab may mainly be effective in patients with poorly chemosensitive disease. Bevacizumab may be prioritized in patients with a high-risk and poorly chemosensitive disease to improve their PFS/OS (patient KELIM score calculator available on the Biomarker Kinetics website).

Pharmacokinetic Model Based on Stochastic Simulation and Estimation for Therapeutic Drug Monitoring of Tacrolimus in Korean Adult Transplant Recipients

BACKGROUND

Tacrolimus shows high variability in inter- and intraindividual pharmacokinetics (PK); therefore, it is important to develop an appropriate model for accurate therapeutic drug monitoring (TDM) procedures. This study aimed to develop a pharmacokinetic model for tacrolimus that can be used for TDM procedures in Korean adult transplant recipients by integrating published models with acquired real-world TDM data and evaluating clinically meaningful covariates.

METHODS

Clinical data of 1829 trough blood samples from 269 subjects were merged with simulated data sets from published models and analyzed using a nonlinear mixed-effect model. The stochastic simulation and estimation (SSE) method was used to obtain the final parameter estimates.

RESULTS

The final estimated values for apparent clearance, the volume of distribution, and absorption rate were 21.2 L/h, 510 L, and 3.1/h, respectively. The number of postoperative days, age, body weight, and type of transplant organs were the major clinical factors affecting tacrolimus PK.

CONCLUSIONS

A tacrolimus PK model that can incorporate published PK models and newly collected data from the Korean population was developed using the SSE method. Despite the limitations in model development owing to the nature of TDM data, the SSE method was useful in retrieving complete information from the TDM data by integrating published PK models while maintaining the variability of the model.

Population Pharmacokinetic Modeling and Probability of Target Attainment of Ceftaroline in Brain and Soft Tissues

Ceftaroline, approved to treat skin infections and pneumonia due to methicillin-resistant Staphylococcus aureus (MRSA), has been considered for the treatment of central nervous system (CNS) infections. A population pharmacokinetic (popPK) model was developed to describe ceftaroline soft tissue and cerebrospinal fluid (CSF) distributions and investigate the probability of target attainment (PTA) of the percentage of the dosing interval that the unbound drug concentration exceeded the MIC (%fT_>MIC) to treat MRSA infections. Healthy subjects' plasma and microdialysate concentrations from muscle and subcutaneous tissue following 600 mg every 12 h (q12h) and q8h and neurosurgical patients' plasma and CSF concentrations following single 600-mg dosing were used. Plasma concentrations were described by a two-compartment model, and tissue concentrations were incorporated as three independent compartments linked to the central compartment by bidirectional transport (clearance in [CL_in] and CL_out). Apparent volumes were fixed to physiological interstitial values. Healthy status and body weight were identified as covariates for the volume of the central compartment, and creatinine clearance was identified for clearance. The CSF glucose concentration (GLUC) was inversely correlated with CL_in,CSF. Simulations showed a PTA of >90% in plasma and soft tissues for both regimens assuming an MIC of 1 mg/L and a %fT_>MIC of 28.8%. Using the same target, patients with inflamed meninges (0.5 < GLUC ≤ 2 mmol/L) would reach PTAs of 99.8% and 97.2% for 600 mg q8h and q12h, respectively. For brain infection with mild inflammation (2 < GLUC ≤ 3.5 mmol/L), the PTAs would be reduced to 34.3% and 9.1%, respectively. Ceftaroline's penetration enhanced by meningeal inflammation suggests that the drug could be a candidate to treat MRSA CNS infections.

Population Pharmacokinetic Modeling and Simulation of Rifapentine Supports Concomitant Antiretroviral Therapy with Efavirenz and Non-Weight Based Dosing

The Brief Rifapentine-Isoniazid Efficacy for TB Prevention/A5279 trial demonstrated a 1-month daily regimen of rifapentine and isoniazid was noninferior to 9 months of isoniazid alone for preventing TB in persons living with HIV (PLWH). Our objective was to evaluate rifapentine pharmacokinetics in trial participants receiving antiretroviral therapy (ART) and perform simulations to compare weight-based rifapentine dosing with a standard, fixed dose. Nonlinear mixed effect modeling was used to estimate rifapentine and 25-desacetyl rifapentine population pharmacokinetic characteristics. The pharmacokinetic model was validated using a nonparametric bootstrap and visual predictive checks. Monte Carlo simulations were performed to compare weight-based and fixed dose regimens. Rifapentine and 25-desacetyl rifapentine concentrations (347 of each; 185 participants) were each described with a one-compartment model with one-way conversion between rifapentine and 25-desacetyl rifapentine. The absorption rate was nearly doubled in fed versus fasting states. Rifapentine clearance was increased 31% in those receiving efavirenz (EFV)-based versus nevirapine-based ART. Metabolite clearance was allometrically scaled with fat-free mass. Simulations showed lower rifapentine exposures with weight-based compared with fixed dosing. With 10 mg/kg weight-based regimens, 26% and 62% of simulated exposures in <35 kg and 35-45 kg weight classes were above target (AUC0 to 24 h of 257 mg*hr/L); 85% of simulated exposures across all weight classes with fixed dosing were above target. These data support fixed dosing with rifapentine 600 mg daily for TB prevention regardless of weight for PLWH 13 years or older receiving the 4-week regimen and no need for dose adjustment when given with EFV-based ART. Clinical Trials Registration. NCT01404312.

Population Pharmacokinetics and Dose Optimization Based on Renal Function of Rivaroxaban in Thai Patients with Non-Valvular Atrial Fibrillation

Low-dose rivaroxaban has been used in Asian patients with direct oral anticoagulants (DOACs) eligible for atrial fibrillation (AF). However, there are few pharmacokinetic (PK) data in Thai patients to support precise dosing. This study aimed to develop a population PK model and determine the optimal rivaroxaban doses in Thai patients. A total of 240 Anti-Xa levels of rivaroxaban from 60 Thai patients were analyzed. A population PK model was established using the nonlinear mixed-effect modeling approach. Monte Carlo simulations were used to predict drug exposures at a steady state for various dosages. Proportions of patients having rivaroxaban exposure within typical exposure ranges were determined. A one-compartment model with first-order absorption best described the data. Creatinine clearance (CrCl) and body weight significantly affected CL/F and V/F, respectively. Regardless of body weight, a higher proportion of patients with CrCl < 50 mL/min receiving the 10-mg once-daily dose had rivaroxaban exposures within the typical exposure ranges. In contrast, a higher proportion of patients with CrCl ≥ 50 mL/min receiving the 15-mg once-daily dose had rivaroxaban exposures within the typical exposure ranges. The study’s findings suggested that low-dose rivaroxaban would be better suited for Thai patients and suggested adjusting the medication’s dose in accordance with renal function.

Improved Confidence in a Confirmatory Stage by Application of Item-Based Pharmacometrics Model: Illustration with a Phase III Active Comparator-Controlled Trial in COPD Patients

PURPOSE

The current study aimed to illustrate how a non-linear mixed effect (NLME) model-based analysis may improve confidence in a Phase III trial through more precise estimates of the drug effect.

METHODS

The FULFIL clinical trial was a Phase III study that compared 24 weeks of once daily inhaled triple therapy with twice daily inhaled dual therapy in patients with chronic obstructive pulmonary disease (COPD). Patient reported outcome data, obtained by using The Evaluating Respiratory Symptoms in COPD (E-RS:COPD) questionnaire, from the FULFIL study were analyzed using an NLME item-based response theory model (IRT). The change from baseline (CFB) in E-RS:COPD total score over 4-week intervals for each treatment arm was obtained using the IRT and compared with published results obtained with a mixed model repeated measures (MMRM) analysis.

RESULTS

The IRT included a graded response model characterizing item parameters and a Weibull function combined with an offset function to describe the COPD symptoms-time course in patients receiving either triple therapy (n = 907) or dual therapy (n = 894). The IRT improved precision of the estimated drug effect compared to MMRM, resulting in a sample size of at least 3.64 times larger for the MMRM analysis to achieve the IRT precision in the CFB estimate.

CONCLUSION

This study shows the advantage of IRT over MMRM with a direct comparison of the same primary endpoint for the two analyses using the same observed clinical trial data, resulting in an increased confidence in Phase III.

Model-Informed Drug Development of New Cefoperazone Sodium and Sulbactam Sodium Combination (3:1): Pharmacokinetic/Pharmacodynamic Analysis and Antibacterial Efficacy Against Enterobacteriaceae

Objective: Cefoperazone/sulbactam is a commonly used antibiotic combination against the extended-spectrum beta-lactamases (ESBLs)-producing bacteria. The objective of this study was to evaluate the efficacy of a new cefoperazone/sulbactam combination (3:1) for Enterobacteriaceae infection via model-informed drug development (MIDD) approaches.

Methods: Sulperazon [cefoperazone/sulbactam (2:1)] was used as a control. Pharmacokinetic (PK) data was collected from a clinical phase I trial. Minimum inhibitory concentrations (MICs) were determined using two-fold broth microdilution method. The percent time that the free drug concentration exceeded the minimum inhibitory concentration (%fT_>MIC) was used as the pharmacokinetic/pharmacodynamic indicator correlated with efficacy. Models were developed to characterize the PK profile of cefoperazone and sulbactam. Monte Carlo simulations were employed to determine the investigational regimens of cefoperazone/sulbactam (3:1) for the treatment of infections caused by Enterobacteriaceae based on the probability of target attainment (PTA) against the tested bacteria.

Results: Two 2-compartment models were developed to describe the PK profiles of cefoperazone and sulbactam. Simulation results following the single-dose showed that the regimens of cefoperazone/sulbactam combinations in the ratios of 3:1 and 2:1 achieved similar PTA against the tested bacteria. Simulation results from the multiple-dose showed that the dosing regimen of cefoperazone/sulbactam (4 g, TID, 3 g:1 g) showed slightly better antibacterial effect than cefoperazone/sulbactam (6 g, BID, 4 g:2 g) against the Escherichia coli (ESBL⁻) and Klebsiella pneumoniae (ESBL⁻). For the other tested bacteria, the above regimens achieved a similar PTA.

Conclusions: Cefoperazone/sulbactam (3:1) showed similar bactericidal activity to sulperazon [cefoperazone/sulbactam (2:1)] against the tested bacteria. For the ESBL-producing and cefoperazone-resistant E. coli and K. pneumoniae, Cefoperazone/sulbactam (3:1) did not exhibit advantage as anticipated. Our study indicated that further clinical trials should be carried out cautiously to avoid the potential risks of not achieving the expected target.

Daptomycin Population Pharmacokinetics in Patients Affected by Severe Gram-Positive Infections: An Update

Daptomycin pharmacokinetics may not depend on renal function only and it significantly differs between healthy volunteers and severely ill patients. Herein, we propose a population pharmacokinetics model based on 424 plasma daptomycin concentrations collected from 156 patients affected by severe Gram-positive infections during a routine therapeutic drug monitoring protocol. Model building and validation were performed using NONMEM 7.2 (ICON plc), Xpose4 and Perl-speaks-to-NONMEM. The final pop-PK model was a one-compartment first-order elimination model, with a 2.7% IIV for drug clearance (Cl), influence of creatinine clearance on drug clearance and of sex on distribution volume. After model validation, we simulated 10,000 patients with the Monte-Carlo method to predict the efficacy and tolerability of different daptomycin daily dosages. For the most common 6 mg/kg daily dose, the simulated probability of overcoming the toxic minimum concentration (24.3 mg/L) was 14.8% and the efficacy (expressed as a cumulative fraction of response) against methicillin-resistant S. aureus, S. pneumoniae and E. faecium was 95.77%, 99.99% and 68%, respectively. According to the model-informed precision dosing paradigm, pharmacokinetic models such as ours could help clinicians to perform patient-tailored antimicrobial dosing and maximize the odds of therapy success without neglecting toxicity risks.

Development of a population pharmacokinetic/pharmacodynamic model for various oral paclitaxel formulations co-administered with ritonavir and thrombospondin-1 based on data from early phase clinical studies

Purpose

Orally administered paclitaxel offers increased patient convenience while providing a method to prolong exposure without long continuous, or repeated, intravenous infusions. The oral bioavailability of paclitaxel is improved through co-administration with ritonavir and application of a suitable pharmaceutical formulation, which addresses the dissolution-limited absorption of paclitaxel. We aimed to characterize the pharmacokinetics of different paclitaxel formulations, co-administered with ritonavir, and to investigate a pharmacodynamic relationship between low-dose metronomic (LDM) treatment with oral paclitaxel and the anti-angiogenic marker thrombospondin-1 (TSP-1).

Methods

Fifty-eight patients treated with different oral paclitaxel formulations were included for pharmacokinetic analysis. Pharmacodynamic data was available for 36 patients. All population pharmacokinetic/pharmacodynamic modelling was performed using non-linear mixed-effects modelling.

Results

A pharmacokinetic model consisting of gut, liver, central, and peripheral compartments was developed for paclitaxel. The gastrointestinal absorption rate was modelled with a Weibull function. Relative gut bioavailabilities of the tablet and capsule formulations, as fractions of the gut bioavailability of the drinking solution, were estimated to be 0.97 (95%CI: 0.67–1.33) and 0.46 (95%CI: 0.34–0.61), respectively. The pharmacokinetic/pharmacodynamic relationship between paclitaxel and TSP-1 was modelled using a turnover model with paclitaxel plasma concentrations driving an increase in TSP-1 formation rate following an E_max relationship with an EC₅₀ of 284 ng/mL (95%CI: 122–724).

Conclusion

The developed pharmacokinetic model adequately described the paclitaxel plasma concentrations for the different oral formulations co-administered with ritonavir. This model, and the established pharmacokinetic/pharmacodynamic relationship with TSP-1, may facilitate future development of oral paclitaxel.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-022-04445-z.