Integrative model-based comparison of target site-specific antimicrobial effects: A case study with ceftaroline and lefamulin

OBJECTIVE

Predictions of antimicrobial effects typically rely on plasma-based pharmacokinetic-pharmacodynamic (PK-PD) targets, ignoring target-site concentrations and potential differences in tissue penetration between antibiotics. In this study, we applied PK-PD modelling to compare target site-specific effects of antibiotics by integrating clinical microdialysis data, in vitro time-kill curves, and antimicrobial susceptibility distributions. As a case study, we compared the effect of lefamulin and ceftaroline against methicillin-resistant Staphylococcus aureus (MRSA) at soft-tissue concentrations.

METHODS

A population PK model describing lefamulin concentrations in plasma, subcutaneous adipose and muscle tissue was developed. For ceftaroline, a similar previously reported PK model was adopted. In vitro time-kill experiments were performed with six MRSA isolates and a PD model was developed to describe bacterial growth and antimicrobial effects. The clinical PK and in vitro PD models were linked to compare antimicrobial effects of ceftaroline and lefamulin at the different target sites.

RESULTS

Considering minimum inhibitory concentration (MIC) distributions and standard dosages, ceftaroline showed superior anti-MRSA effects compared to lefamulin both at plasma and soft-tissue concentrations. Looking at the individual antibiotics, lefamulin effects were highest at soft-tissue concentrations, while ceftaroline effects were highest at plasma concentrations, emphasising the importance of considering target-site PK-PD in antibiotic treatment optimisation.

CONCLUSION

Given standard dosing regimens, ceftaroline appeared more effective than lefamulin against MRSA at soft-tissue concentrations. The PK-PD model-based approach applied in this study could be used to compare or explore the potential of antibiotics for specific indications or in populations with unique target-site PK.

Virtual Patient Simulation Using Copula Modeling

Virtual patient simulation is increasingly performed to support model-based optimization of clinical trial designs or individualized dosing strategies. Quantitative pharmacological models typically incorporate individual-level patient characteristics, or covariates, which enable the generation of virtual patient cohorts. The individual-level patient characteristics, or covariates, used as input for such simulations should accurately reflect the values seen in real patient populations. Current methods often make unrealistic assumptions about the correlation between patient's covariates or require direct access to actual data sets with individual-level patient data, which may often be limited by data sharing limitations. We propose and evaluate the use of copulas to address current shortcomings in simulation of patient-associated covariates for virtual patient simulations for model-based dose and trial optimization in clinical pharmacology. Copulas are multivariate distribution functions that can capture joint distributions, including the correlation, of covariate sets. We compare the performance of copulas to alternative simulation strategies, and we demonstrate their utility in several case studies. Our work demonstrates that copulas can reproduce realistic patient characteristics, both in terms of individual covariates and the dependence structure between different covariates, outperforming alternative methods, in particular when aiming to reproduce high-dimensional covariate sets. In conclusion, copulas represent a versatile and generalizable approach for virtual patient simulation which preserve relationships between covariates, and offer an open science strategy to facilitate re-use of patient data sets.

Predictions of Bedaquiline Central Nervous System Exposure in Patients with Tuberculosis Meningitis Using Physiologically based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The use of bedaquiline as a treatment option for drug-resistant tuberculosis meningitis (TBM) is of interest to address the increased prevalence of resistance to first-line antibiotics. To this end, we describe a whole-body physiologically based pharmacokinetic (PBPK) model for bedaquiline to predict central nervous system (CNS) exposure.

METHODS

A whole-body PBPK model was developed for bedaquiline and its metabolite, M2. The model included compartments for brain and cerebrospinal fluid (CSF). Model predictions were evaluated by comparison to plasma PK time profiles following different dosing regimens and sparse CSF concentrations data from patients. Simulations were then conducted to compare CNS and lung exposures to plasma exposure at clinically relevant dosing schedules.

RESULTS

The model appropriately described the observed plasma and CSF bedaquiline and M2 concentrations from patients with pulmonary tuberculosis (TB). The model predicted a high impact of tissue binding on target site drug concentrations in CNS. Predicted unbound exposures within brain interstitial exposures were comparable with unbound vascular plasma and unbound lung exposures. However, unbound brain intracellular exposures were predicted to be 7% of unbound vascular plasma and unbound lung intracellular exposures.

CONCLUSIONS

The whole-body PBPK model for bedaquiline and M2 predicted unbound concentrations in brain to be significantly lower than the unbound concentrations in the lung at clinically relevant doses. Our findings suggest that bedaquiline may result in relatively inferior efficacy against drug-resistant TBM when compared with efficacy against drug-resistant pulmonary TB.

Longitudinal metabolite profiling of Streptococcus pneumoniae-associated community-acquired pneumonia

INTRODUCTION

Longitudinal biomarkers in patients with community-acquired pneumonia (CAP) may help in monitoring of disease progression and treatment response. The metabolic host response could be a potential source of such biomarkers since it closely associates with the current health status of the patient.

OBJECTIVES

In this study we performed longitudinal metabolite profiling in patients with CAP for a comprehensive range of metabolites to identify potential host response biomarkers.

METHODS

Previously collected serum samples from CAP patients with confirmed Streptococcus pneumoniae infection (n = 25) were used. Samples were collected at multiple time points, up to 30 days after admission. A wide range of metabolites was measured, including amines, acylcarnitines, organic acids, and lipids. The associations between metabolites and C-reactive protein (CRP), procalcitonin, CURB disease severity score at admission, and total length of stay were evaluated.

RESULTS

Distinct longitudinal profiles of metabolite profiles were identified, including cholesteryl esters, diacyl-phosphatidylethanolamine, diacylglycerols, lysophosphatidylcholines, sphingomyelin, and triglycerides. Positive correlations were found between CRP and phosphatidylcholine (34:1) (cor = 0.63) and negative correlations were found for CRP and nine lysophosphocholines (cor = - 0.57 to - 0.74). The CURB disease severity score was negatively associated with six metabolites, including acylcarnitines (tau = - 0.64 to - 0.58). Negative correlations were found between the length of stay and six triglycerides (TGs), especially TGs (60:3) and (58:2) (cor = - 0.63 and - 0.61).

CONCLUSION

The identified metabolites may provide insight into biological mechanisms underlying disease severity and may be of interest for exploration as potential treatment response monitoring biomarker.

Model-based dose optimization framework for bedaquiline, pretomanid and linezolid for the treatment of drug-resistant tuberculosis

AIMS

Bedaquiline, pretomanid and linezolid (BPaL) combination treatment against Mycobacterium tuberculosis is promising, yet safety and adherence concerns exist that motivate exploration of alternative dosing regimens. We developed a mechanistic modelling framework to compare the efficacy of the current and alternative BPaL treatment strategies.

METHODS

Pharmacodynamic models for each drug in the BPaL combination treatment were developed using in vitro time-kill data. These models were combined with pharmacokinetic models, incorporating body weight, lesion volume, site-of-action distribution, bacterial susceptibility and pharmacodynamic interactions to assemble the framework. The model was qualified by comparing the simulations against the observed clinical data. Simulations were performed evaluating bedaquiline and linezolid approved (bedaquiline 400 mg once daily [QD] for 14 days followed by 200 mg three times a week, linezolid 1200 mg QD) and alternative dosing regimens (bedaquiline 200 mg QD, linezolid 600 mg QD).

RESULTS

The framework adequately described the observed antibacterial activity data in patients following monotherapy for each drug and approved BPaL dosing. The simulations suggested a minor difference in median time to colony forming unit (CFU)-clearance state with the bedaquiline alternative compared to the approved dosing and the linezolid alternative compared to the approved dosing. Median time to non-replicating-clearance state was predicted to be 15 days from the CFU-clearance state.

CONCLUSIONS

The model-based simulations suggested that comparable efficacy can be achieved using alternative bedaquiline and linezolid dosing, which may improve safety and adherence in drug-resistant tuberculosis patients. The framework can be utilized to evaluate treatment optimization approaches, including dosing regimen and duration of treatment predictions to eradicate both replicating- and non-replicating bacteria from lung and lesions.

Semi-mechanistic modeling of resistance development to β-lactam and β-lactamase-inhibitor combinations

The use of β-lactam (BL) and β-lactamase inhibitor (BLI) combinations, such as piperacillin-tazobactam (PIP-TAZ) is an effective strategy to combat infections by extended-spectrum β-lactamase-producing bacteria. However, in Gram-negative bacteria, resistance (both mutational and adaptive) to BL-BLI combination can still develop through multiple mechanisms. These mechanisms may include increased β-lactamase activity, reduced drug influx, and increased drug efflux. Understanding the relative contribution of these mechanisms during resistance development helps identify the most impactful mechanism to target in designing a treatment to counter BL-BLI resistance. This study used semi-mechanistic mathematical modeling in combination with antibiotic sensitivity assays to assess the potential impact of different resistance mechanisms during the development of PIP-TAZ resistance in a Klebsiella pneumoniae isolate expressing CTX-M-15 and SHV-1 β-lactamases. The mathematical models were used to evaluate the potential impact of several cellular changes as a sole mediator of PIP-TAZ resistance. Our semi-mechanistic model identified 2 out of the 13 inspected mechanisms as key resistance mechanisms that may independently support the observed magnitude of PIP-TAZ resistance, namely porin loss and efflux pump up-regulation. Simulation using the resulting models also suggested the possible adjustment of PIP-TAZ dose outside its commonly used 8:1 dosing ratio. The current study demonstrated how theory-based mechanistic models informed by experimental data can be used to support hypothesis generation regarding potential resistance mechanisms, which may guide subsequent experimental studies.

Unraveling the Effects of Acute Inflammation on Pharmacokinetics: A Model-Based Analysis Focusing on Renal Glomerular Filtration Rate and Cytochrome P450 3A4-Mediated Metabolism

BACKGROUND AND OBJECTIVES: Acute inflammation caused by infections or sepsis can impact pharmacokinetics. We used a model-based analysis to evaluate the effect of acute inflammation as represented by interleukin-6 (IL-6) levels on drug clearance, focusing on renal glomerular filtration rate (GFR) and cytochrome P450 3A4 (CYP3A4)-mediated metabolism.

METHODS

A physiologically based model incorporating renal and hepatic drug clearance was implemented. Functions correlating IL-6 levels with GFR and in vitro CYP3A4 activity were derived and incorporated into the modeling framework. We then simulated treatment scenarios for hypothetical drugs by varying the IL-6 levels, the contribution of renal and hepatic drug clearance, and protein binding. The relative change in observed area under the concentration-time curve (AUC) was computed for these scenarios.

RESULTS

Inflammation showed opposite effects on drug exposure for drugs eliminated via the liver and kidney, with the effect of inflammation being inversely proportional to the extraction ratio (ER). For renally cleared drugs, the relative decrease in AUC was close to 30% during severe inflammation. For CYP3A4 substrates, the relative increase in AUC could exceed 50% for low-ER drugs. Finally, the impact of inflammation-induced changes in drug clearance is smaller for drugs with a larger unbound fraction.

CONCLUSION

This analysis demonstrates differences in the impact of inflammation on drug clearance for different drug types. The effects of inflammation status on pharmacokinetics may explain the inter-individual variability in pharmacokinetics in critically ill patients. The proposed model-based analysis may be used to further evaluate the effect of inflammation, i.e., by incorporating the effect of inflammation on other drug-metabolizing enzymes or physiological processes.

Predictions of Bedaquiline and Pretomanid Target Attainment in Lung Lesions of Tuberculosis Patients using Translational Minimal Physiologically Based Pharmacokinetic Modeling

BACKGROUND

Site-of-action concentrations for bedaquiline and pretomanid from tuberculosis patients are unavailable. The objective of this work was to predict bedaquiline and pretomanid site-of-action exposures using a translational minimal physiologically based pharmacokinetic (mPBPK) approach to understand the probability of target attainment (PTA).

METHODS

A general translational mPBPK framework for the prediction of lung and lung lesion exposure was developed and validated using pyrazinamide site-of-action data from mice and humans. We then implemented the framework for bedaquiline and pretomanid. Simulations were conducted to predict site-of-action exposures following standard bedaquiline and pretomanid, and bedaquiline once-daily dosing. Probabilities of average concentrations within lesions and lungs greater than the minimum bactericidal concentration for non-replicating (MBCNR) and replicating (MBCR) bacteria were calculated. Effects of patient-specific differences on target attainment were evaluated.

RESULTS

The translational modeling approach was successful in predicting pyrazinamide lung concentrations from mice to patients. We predicted that 94% and 53% of patients would attain bedaquiline average daily PK exposure within lesions (Cavg-lesion) > MBCNR during the extensive phase of bedaquiline standard (2 weeks) and once-daily (8 weeks) dosing, respectively. Less than 5% of patients were predicted to achieve Cavg-lesion > MBCNR during the continuation phase of bedaquiline or pretomanid treatment, and more than 80% of patients were predicted to achieve Cavg-lung >MBCR for all simulated dosing regimens of bedaquiline and pretomanid.

CONCLUSIONS

The translational mPBPK model predicted that the standard bedaquiline continuation phase and standard pretomanid dosing may not achieve optimal exposures to eradicate non-replicating bacteria in most patients.

Associations Between Symptoms, Donor Characteristics and IgG Antibody Response in 2082 COVID-19 Convalescent Plasma Donors

Many studies already reported on the association between patient characteristics on the severity of COVID-19 disease outcome, but the relation with SARS-CoV-2 antibody levels is less clear. To investigate this in more detail, we performed a retrospective observational study in which we used the IgG antibody response from 11,118 longitudinal antibody measurements of 2,082 unique COVID convalescent plasma donors. COVID-19 symptoms and donor characteristics were obtained by a questionnaire. Antibody responses were modelled using a linear mixed-effects model. Our study confirms that the SARS-CoV-2 antibody response is associated with patient characteristics like body mass index and age. Antibody decay was faster in male than in female donors (average half-life of 62 versus 72 days). Most interestingly, we also found that three symptoms (headache, anosmia, nasal cold) were associated with lower peak IgG, while six other symptoms (dry cough, fatigue, diarrhoea, fever, dyspnoea, muscle weakness) were associated with higher IgG concentrations.

Physiologically Based Modelling Framework for Prediction of Pulmonary Pharmacokinetics of Antimicrobial Target Site Concentrations

BACKGROUND AND OBJECTIVES

Prediction of antimicrobial target-site pharmacokinetics is of relevance to optimize treatment with antimicrobial agents. A physiologically based pharmacokinetic (PBPK) model framework was developed for prediction of pulmonary pharmacokinetics, including key pulmonary infection sites (i.e. the alveolar macrophages and the epithelial lining fluid).

METHODS

The modelling framework incorporated three lung PBPK models: a general passive permeability-limited model, a drug-specific permeability-limited model and a quantitative structure-property relationship (QSPR)-informed perfusion-limited model. We applied the modelling framework to three fluoroquinolone antibiotics. Incorporation of experimental drug-specific permeability data was found essential for accurate prediction.

RESULTS

In the absence of drug-specific transport data, our QSPR-based model has generic applicability. Furthermore, we evaluated the impact of drug properties and pathophysiologically related changes on pulmonary pharmacokinetics. Pulmonary pharmacokinetics were highly affected by physiological changes, causing a shift in the main route of diffusion (i.e. paracellular or transcellular). Finally, we show that lysosomal trapping can cause an overestimation of cytosolic concentrations for basic compounds when measuring drug concentrations in cell homogenate.

CONCLUSION

The developed lung PBPK model framework constitutes a promising tool for characterization of pulmonary exposure of systemically administrated antimicrobials.

Targeted temperature management after cardiac arrest is associated with reduced metabolism of pantoprazole - A probe drug of CYP2C19 metabolism

OBJECTIVE

Targeted temperature management (TTM) is part of standard post-resuscitation care. TTM may downregulate cytochrome enzyme activity and thus impact drug metabolism. This study compared the pharmacokinetics (PK) of pantoprazole, a probe drug of CYP2C19-dependent metabolism, at different stages of TTM following cardiac arrest.

METHODS

This prospective controlled study was performed at the Medical University of Vienna and enrolled 16 patients following cardiac arrest. The patients completed up to three study periods (each lasting 24 h) in which plasma concentrations of pantoprazole were quantified: (P1) hypothermia (33 °C) after admission, (P2) normothermia after rewarming (36 °C, intensive care), and (P3) normothermia during recovery (normal ward, control group). PK was analysed using non-compartmental analysis and nonlinear mixed-effects modelling.

RESULTS

16 patients completed periods P1 and P2; ten completed P3. The median half-life of pantoprazole was 2.4 h (quartiles: 1.8-4.8 h) in P1, 2.8 h (2.1-6.8 h, p = 0.046 vs. P1, p = 0.005 vs. P3) in P2 and 1.2 h (0.9 - 2.3 h, p = 0.007 vs. P1) in P3. A two-compartment model described the PK data best. Typical values for clearance were estimated separately for each study period, indicating 40% and 29% reductions during P1 and P2, respectively, compared to P3. The central volume of distribution was estimated separately for P2, indicating a 64% increase compared to P1 and P3.

CONCLUSION

CYP2C19-dependent drug metabolism is downregulated during TTM following cardiac arrest. These results may influence drug choice and dosing of similarly metabolized drugs and may be helpful for designing studies in similar clinical situations.

Host-directed therapies for tuberculosis: quantitative systems pharmacology approaches

Host-directed therapies (HDTs) that modulate host-pathogen interactions offer an innovative strategy to combat Mycobacterium tuberculosis (Mtb) infections. When combined with tuberculosis (TB) antibiotics, HDTs could contribute to improving treatment outcomes, reducing treatment duration, and preventing resistance development. Translation of the interplay of host-pathogen interactions leveraged by HDTs towards therapeutic outcomes in patients is challenging. Quantitative understanding of the multifaceted nature of the host-pathogen interactions is vital to rationally design HDT strategies. Here, we (i) provide an overview of key Mtb host-pathogen interactions as basis for HDT strategies; and (ii) discuss the components and utility of quantitative systems pharmacology (QSP) models to inform HDT strategies. QSP models can be used to identify and optimize treatment targets, to facilitate preclinical to human translation, and to design combination treatment strategies.

Identification of antibiotic collateral sensitivity and resistance interactions in population surveillance data

BACKGROUND

Collateral effects of antibiotic resistance occur when resistance to one antibiotic agent leads to increased resistance or increased sensitivity to a second agent, known respectively as collateral resistance (CR) and collateral sensitivity (CS). Collateral effects are relevant to limit impact of antibiotic resistance in design of antibiotic treatments. However, methods to detect antibiotic collateral effects in clinical population surveillance data of antibiotic resistance are lacking.

OBJECTIVES

To develop a methodology to quantify collateral effect directionality and effect size from large-scale antimicrobial resistance population surveillance data.

METHODS

We propose a methodology to quantify and test collateral effects in clinical surveillance data based on a conditional t-test. Our methodology was evaluated using MIC data for 419 Escherichia coli strains, containing MIC data for 20 antibiotics, which were obtained from the Pathosystems Resource Integration Center (PATRIC) database.

RESULTS

We demonstrate that the proposed approach identifies several antibiotic combinations that show symmetrical or non-symmetrical CR and CS. For several of these combinations, collateral effects were previously confirmed in experimental studies. We furthermore provide insight into the power of our method for multiple collateral effect sizes and MIC distributions.

CONCLUSIONS

Our proposed approach is of relevance as a tool for analysis of large-scale population surveillance studies to provide broad systematic identification of collateral effects related to antibiotic resistance, and is made available to the community as an R package. This method can help mapping CS and CR, which could guide combination therapy and prescribing in the future.

Design principles of collateral sensitivity-based dosing strategies

Collateral sensitivity (CS)-based antibiotic treatments, where increased resistance to one antibiotic leads to increased sensitivity to a second antibiotic, may have the potential to limit the emergence of antimicrobial resistance. However, it remains unclear how to best design CS-based treatment schedules. To address this problem, we use mathematical modelling to study the effects of pathogen- and drug-specific characteristics for different treatment designs on bacterial population dynamics and resistance evolution. We confirm that simultaneous and one-day cycling treatments could supress resistance in the presence of CS. We show that the efficacy of CS-based cycling therapies depends critically on the order of drug administration. Finally, we find that reciprocal CS is not essential to suppress resistance, a result that significantly broadens treatment options given the ubiquity of one-way CS in pathogens. Overall, our analyses identify key design principles of CS-based treatment strategies and provide guidance to develop treatment schedules to suppress resistance.

Identification of high-dimensional omics-derived predictors for tumor growth dynamics using machine learning and pharmacometric modeling

Pharmacometric modeling can capture tumor growth inhibition (TGI) dynamics and variability. These approaches do not usually consider covariates in high-dimensional settings, whereas high-dimensional molecular profiling technologies ("omics") are being increasingly considered for prediction of anticancer drug treatment response. Machine learning (ML) approaches have been applied to identify high-dimensional omics predictors for treatment outcome. Here, we aimed to combine TGI modeling and ML approaches for two distinct aims: omics-based prediction of tumor growth profiles and identification of pathways associated with treatment response and resistance. We propose a two-step approach combining ML using least absolute shrinkage and selection operator (LASSO) regression with pharmacometric modeling. We demonstrate our workflow using a previously published dataset consisting of 4706 tumor growth profiles of patient-derived xenograft (PDX) models treated with a variety of mono- and combination regimens. Pharmacometric TGI models were fit to the tumor growth profiles. The obtained empirical Bayes estimates-derived TGI parameter values were regressed using the LASSO on high-dimensional genomic copy number variation data, which contained over 20,000 variables. The predictive model was able to decrease median prediction error by 4% as compared with a model without any genomic information. A total of 74 pathways were identified as related to treatment response or resistance development by LASSO, of which part was verified by literature. In conclusion, we demonstrate how the combined use of ML and pharmacometric modeling can be used to gain pharmacological understanding in genomic factors driving variation in treatment response.

Clinical Pharmacokinetics and Pharmacodynamics of Immune Checkpoint Inhibitors

Immune checkpoint inhibitors (ICIs) have demonstrated significant clinical impact in improving overall survival of several malignancies associated with poor outcomes; however, only 20–40% of patients will show long-lasting survival. Further clarification of factors related to treatment response can support improvements in clinical outcome and guide the development of novel immune checkpoint therapies. In this article, we have provided an overview of the pharmacokinetic (PK) aspects related to current ICIs, which include target-mediated drug disposition and time-varying drug clearance. In response to the variation in treatment exposure of ICIs and the significant healthcare costs associated with these agents, arguments for both dose individualization and generalization are provided. We address important issues related to the efficacy and safety, the pharmacodynamics (PD), of ICIs, including exposure–response relationships related to clinical outcome. The unique PK and PD aspects of ICIs give rise to issues of confounding and suboptimal surrogate endpoints that complicate interpretation of exposure–response analysis. Biomarkers to identify patients benefiting from treatment with ICIs have been brought forward. However, validated biomarkers to monitor treatment response are currently lacking.

Host-response biomarkers for the diagnosis of bacterial respiratory tract infections

Appropriate antibiotic treatment for respiratory tract infections (RTIs) necessitates rapid and accurate diagnosis of microbial etiology, which remains challenging despite recent innovations. Several host response-based biomarkers due to infection have been suggested to allow discrimination of bacterial and non-bacterial microbial RTI etiology. This review provides an overview of clinical studies that investigated the diagnostic performance of host-response proteomic biomarkers to identify RTI microbial etiology. Procalcitonin and C-reactive protein have been studied most extensively; whereof procalcitonin has demonstrated the strongest diagnostic performance compared to other biomarkers. Proadrenomedullin, soluble triggering receptor expressed on myeloid cells-1, neopterin and pentraxin-3 need more studies to confirm their diagnostic value. For syndecan-4 and lipocalin-2 currently insufficient evidence exists. Common limitations in several of the studies were the relatively small scale setting, heterogeneous patient population and the absence of statistical power calculation.

A Population Pharmacokinetic Model of Oral Docetaxel Coadministered With Ritonavir to Support Early Clinical Development

Oral administration of docetaxel is an attractive alternative for conventional intravenous (IV) administration. The low bioavailability of docetaxel, however, hinders the application of oral docetaxel in the clinic. The aim of the current study was to develop a population pharmacokinetic (PK) model for docetaxel and ritonavir based on the phase 1 studies and to support drug development of this combination treatment. PK data were collected from 191 patients who received IV docetaxel and different oral docetaxel formulations (drinking solution, ModraDoc001 capsule, and ModraDoc006 tablet) coadministered with ritonavir. A PK model was first developed for ritonavir. Subsequently, a semiphysiological PK model was developed for docetaxel, which incorporated the inhibition of docetaxel metabolism by ritonavir. The uninhibited intrinsic clearance of docetaxel was estimated based on data on IV docetaxel as 1980 L/h (relative standard error, 11%). Ritonavir coadministration extensively inhibited the hepatic metabolism of docetaxel to 9.3%, which resulted in up to 12-fold higher docetaxel plasma concentrations compared to oral docetaxel coadministered without ritonavir. In conclusion, a semiphysiological PK model for docetaxel and ritonavir was successfully developed. Coadministration of ritonavir resulted in increased plasma concentrations of docetaxel after administration of the oral formulations of ModraDoc. Furthermore, the oral ModraDoc formulations showed lower variability in plasma concentrations between and within patients compared to the drinking solution. Comparable exposure could be reached with the oral ModraDoc formulations compared to IV administration.

Disruption of podocyte cytoskeletal biomechanics by dasatinib leads to nephrotoxicity

Nephrotoxicity is a critical adverse event that leads to discontinuation of kinase inhibitor (KI) treatment. Here we show, through meta-analyses of FDA Adverse Event Reporting System, that dasatinib is associated with high risk for glomerular toxicity that is uncoupled from hypertension, suggesting a direct link between dasatinib and podocytes. We further investigate the cellular effects of dasatinib and other comparable KIs with varying risks of nephrotoxicity. Dasatinib treated podocytes show significant changes in focal adhesions, actin cytoskeleton, and morphology that are not observed with other KIs. We use phosphoproteomics and kinome profiling to identify the molecular mechanisms of dasatinib-induced injury to the actin cytoskeleton, and atomic force microscopy to quantify impairment to cellular biomechanics. Furthermore, chronic administration of dasatinib in mice causes reversible glomerular dysfunction, loss of stress fibers, and foot process effacement. We conclude that dasatinib induces nephrotoxicity through altered podocyte actin cytoskeleton, leading to injurious cellular biomechanics.

Systems Pharmacology: Defining the Interactions of Drug Combinations

The majority of diseases are associated with alterations in multiple molecular pathways and complex interactions at the cellular and organ levels. Single-target monotherapies therefore have intrinsic limitations with respect to their maximum therapeutic benefits. The potential of combination drug therapies has received interest for the treatment of many diseases and is well established in some areas, such as oncology. Combination drug treatments may allow us to identify synergistic drug effects, reduce adverse drug reactions, and address variability in disease characteristics between patients. Identification of combination therapies remains challenging. We discuss current state-of-the-art systems pharmacology approaches to enable rational identification of combination therapies. These approaches, which include characterization of mechanisms of disease and drug action at a systems level, can enable understanding of drug interactions at the molecular, cellular, physiological, and organismal levels. Such multiscale understanding can enable precision medicine by promoting the rational development of combination therapy at the level of individual patients for many diseases.

Crowdsourced Asparagus Urinary Odor Population Kinetics

The consumption of asparagus is associated with the production of malodorous urine with considerable interindividual variability (IIV). To characterize the urinary odor kinetics after consumption of asparagus spears, we conducted a study with consenting attendees from two American Society for Clinical Pharmacology and Therapeutics (ASCPT) meetings. Subjects were randomized to eat a specific number of asparagus spears, and then asked to report their urinary odor perception. Eighty‐seven subjects were included in the final analysis. A mixed effect proportional odds model was developed that adequately characterized the dose‐response relationship. We estimated the half‐life of the asparagus effect on malodorous urine to be 4.7 hours (relative standard error (RSE) = 13.2%), and identified a dose‐response slope term with good precision (24.3%). Age was found as the predictor of IIV in slope estimates. This study design and tools can be used as a demonstration “crowdsourcing” project for studying population kinetics in organizational and educational settings.

Population pharmacokinetics of oxcarbazepine and its metabolite 10-hydroxycarbazepine in healthy subjects

Oxcarbazepine is indicated for the treatment of partial or generalised tonic-clonic seizures. Most of the absorbed oxcarbazepine is converted into its active metabolite, 10-hydroxycarbazepine (MHD), which can exist as R-(-)- and S-(+)-MHD enantiomers. Here we describe the influence of the P-glycoprotein (P-gp) inhibitor verapamil, on the disposition of oxcarbazepine and MHD enantiomers, both of which are P-gp substrates. Healthy subjects (n=12) were randomised to oxcarbazepine or oxcarbazepine combined with verapamil at doses of 300mg b.i.d. and 80mg t.i.d., respectively. Blood samples (n=185) were collected over a period of 12h post oxcarbazepine dose. An integrated PK model was developed using nonlinear mixed effects modelling using a meta-analytical approach. The pharmacokinetics of oxcarbazepine was described by a two-compartment model with absorption transit compartments and first-order elimination. The concentration-time profiles of both MHD enantiomers were characterised by a one-compartment distribution model. Clearance estimates (95% CI) were 84.9L/h (69.5-100.3) for oxcarbazepine and 2.0L/h (1.9-2.1) for both MHD enantiomers. The volume of distribution was much larger for oxcarbazepine (131L (97-165)) as compared to R-(-)- and S-(+)-MHD (23.6L (14.4-32.8) vs. 31.7L (22.5-40.9), respectively). Co-administration of verapamil resulted in a modest increase of the apparent bioavailability of oxcarbazepine by 12% (10-28), but did not affect parent or metabolite clearances. Despite the evidence of comparable systemic levels of OXC and MHD following administration of verapamil, differences in brain exposure to both moieties cannot be excluded after P-glycoprotein inhibition.

Evidence-based drug treatment for special patient populations through model-based approaches

The majority of marketed drugs remain understudied in some patient populations such as pregnant women, paediatrics, the obese, the critically-ill, and the elderly. As a consequence, currently used dosing regimens may not assure optimal efficacy or minimal toxicity in these patients. Given the vulnerability of some subpopulations and the challenges and costs of performing clinical studies in these populations, cutting-edge approaches are needed to effectively develop evidence-based and individualized drug dosing regimens. Five key issues are presented that are essential to support and expedite the development of drug dosing regimens in these populations using model-based approaches: 1) model development combined with proper validation procedures to extract as much valid information from available study data as possible, with limited burden to patients and costs; 2) integration of existing data and the use of prior pharmacological and physiological knowledge in study design and data analysis, to further develop knowledge and avoid unnecessary or unrealistic (large) studies in vulnerable populations; 3) clinical proof-of-principle in a prospective evaluation of a developed drug dosing regimen, to confirm that a newly proposed regimen indeed results in the desired outcomes in terms of drug concentrations, efficacy, and/or safety; 4) pharmacodynamics studies in addition to pharmacokinetics studies for drugs for which a difference in disease progression and/or in exposure-response relation is anticipated compared to the reference population; 5) additional efforts to implement developed dosing regimens in clinical practice once drug pharmacokinetics and pharmacodynamics have been characterized in special patient populations. The latter remains an important bottleneck, but this is essential to truly realize evidence-based and individualized drug dosing for special patient populations. As all tools required for this purpose are available, we have the moral and societal obligation to make safe and effective pharmacotherapy available for these patients too.

Evaluating renal function and age as predictors of amikacin clearance in neonates: model-based analysis and optimal dosing strategies

AIMS

We aimed to compare the performance of renal function and age as predictors of inter-individual variability (IIV) in clearance of amikacin in neonates through parallel development of population pharmacokinetic (PK) models and their associated impact on optimal dosing regimens.

METHODS

Amikacin concentrations were retrospectively collected for 149 neonates receiving amikacin (post-natal age (PNA) between 4-89 days). Two population PK models were developed in parallel, considering at least as predictors current body weight (WT), in combination with either creatinine clearance (CLcr ) or age descriptors. Using stochastic simulations for both renal function or age-based dosing, we identified optimal dosing strategies that were based on attainment of optimal peak- (PCC) and trough target concentration coverage (TCC) windows associated with efficacy and toxicity.

RESULTS

The CLcr and age-based population PK models both included current body weight (WT) on CL, central distribution volume and intercompartmental clearance, in combination with either CLcr or PNA as predictors for IIV of clearance (CL). The WT-CLcr model explained 6.9% more IIV in CL compared with the WT-PNA model. Both models successfully described an external dataset (n = 53) of amikacin PK. The simulation analysis of optimal dose regimens suggested similar performance of either CLcr or PNA based dosing.

CONCLUSION

CLcr predicted more IIV in CL, but did not translate into clinically relevant improvements of target concentrations. Our optimized dose regimens can be considered for further evaluation to optimize initial treatment with amikacin.

Pronounced between-subject and circadian variability in thymidylate synthase and dihydropyrimidine dehydrogenase enzyme activity in human volunteers

AIMS

The enzymatic activity of dihydropyrimidine dehydrogenase (DPD) and thymidylate synthase (TS) are important for the tolerability and efficacy of the fluoropyrimidine drugs. In the present study, we explored between-subject variability (BSV) and circadian rhythmicity in DPD and TS activity in human volunteers.

METHODS

The BSVs in DPD activity (n = 20) in peripheral blood mononuclear cells (PBMCs) and in plasma, measured by means of the dihydrouracil (DHU) and uracil (U) plasma levels and DHU : U ratio (n = 40), and TS activity in PBMCs (n = 19), were examined. Samples were collected every 4 h throughout 1 day for assessment of circadian rhythmicity in DPD and TS activity in PBMCs (n = 12) and DHU : U plasma ratios (n = 23). In addition, the effects of genetic polymorphisms and gene expression on DPD and TS activity were explored.

RESULTS

Population mean (± standard deviation) DPD activity in PBMCs and DHU : U plasma ratio were 9.2 (±2.1) nmol mg(-1) h(-1) and 10.6 (±2.4), respectively. Individual TS activity in PBMCs ranged from 0.024 nmol mg(-1) h(-1) to 0.596 nmol mg(-1) h(-1) . Circadian rhythmicity was demonstrated for all phenotype markers. Between 00:30 h and 02:00 h, DPD activity in PBMCs peaked, while the DHU : U plasma ratio and TS activity in PBMCs showed trough activity. Peak-to-trough ratios for DPD and TS activity in PBMCs were 1.69 and 1.62, respectively. For the DHU : U plasma ratio, the peak-to-trough ratio was 1.43.

CONCLUSIONS

BSV and circadian variability in DPD and TS activity were demonstrated. Circadian rhythmicity in DPD might be tissue dependent. The results suggested an influence of circadian rhythms on phenotype-guided fluoropyrimidine dosing and supported implications for chronotherapy with high-dose fluoropyrimidine administration during the night.

Pooled population pharmacokinetic model of imipenem in plasma and the lung epithelial lining fluid

AIMS

Several clinical trials have confirmed the therapeutic benefit of imipenem for treatment of lung infections. There is however no knowledge of the penetration of imipenem into the lung epithelial lining fluid (ELF), the site of action relevant for lung infections. Furthermore, although the plasma pharmacokinetics (PK) of imipenem has been widely studied, most studies have been based on selected patient groups. The aim of this analysis was to characterize imipenem plasma PK across populations and to quantify imipenem ELF penetration.

METHODS

A population model for imipenem plasma PK was developed using data obtained from healthy volunteers, elderly subjects and subjects with renal impairment, in order to identify predictors for inter-individual variability (IIV) of imipenem PK. Subsequently, a clinical study which measured plasma and ELF concentrations of imipenem was included in order to quantify lung penetration.

RESULTS

A two compartmental model best described the plasma PK of imipenem. Creatinine clearance and body weight were included as subject characteristics predictive for IIV on clearance. Typical estimates for clearance, central and peripheral volume, and inter-compartmental clearance were 11.5 l h(-1) , 9.37 l, 6.41 l, 13.7 l h(-1) , respectively (relative standard error (RSE) <8%). The distribution of imipenem into ELF was described using a time-independent penetration coefficient of 0.44 (RSE 14%).

CONCLUSION

The identified lung penetration coefficient confirms the clinical relevance of imipenem for treatment of lung infections, while the population PK model provided insights into predictors of IIV for imipenem PK and may be of relevance to support dose optimization in various subject groups.

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

AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Towards integrative systems pharmacology models in oncology drug development

Quantitative systems pharmacology (QSP) modeling represents an emerging area of value to further streamline knowledge integration and to better inform decision making processes in drug development. QSP models reside at the interface between systems biology models and pharmacological models, yet their concrete implementation still needs to be established further. This review outlines key modeling techniques in both of these areas and to subsequently discuss challenges and opportunities for further integration, in oncology drug development.

Design of a drug-drug interaction study of vincristine with azole antifungals in pediatric cancer patients using clinical trial simulation

BACKGROUND

The aim of the current work was to perform a clinical trial simulation (CTS) analysis to optimize a drug-drug interaction (DDI) study of vincristine in children who also received azole antifungals, taking into account challenges of conducting clinical trials in this population, and, to provide a motivating example of the application of CTS in the design of pediatric oncology clinical trials.

PROCEDURE

A pharmacokinetic (PK) model for vincristine in children was used to simulate concentration-time profiles. A continuous model for body surface area versus age was defined based on pediatric growth curves. Informative sampling time windows were derived using D-optimal design. The CTS framework was used to different magnitudes of clearance inhibition (10%, 25%, or 40%), sample size (30-500), the impact of missing samples or sampling occasions, and the age distribution, on the power to detect a significant inhibition effect, and in addition, the relative estimation error (REE) of the interaction effect.

RESULTS

A minimum group specific sample size of 38 patients with a total sample size of 150 patients was required to detect a clearance inhibition effect of 40% with 80% power, while in the case of a lower effect of clearance inhibition, a substantially larger sample size was required. However, for the majority of re-estimated drug effects, the inhibition effect could be estimated precisely (REE < 25%) in even smaller sample sizes and with lower effect sizes.

CONCLUSION

This work demonstrated the utility of CTS for the evaluation of PK clinical trial designs in the pediatric oncology population.

Population pharmacokinetic-pharmacodynamic analysis for eribulin mesilate-associated neutropenia

AIMS

Eribulin mesilate is an inhibitor of microtubule dynamics that is approved for the treatment of late-stage metastatic breast cancer. Neutropenia is one of the major dose-limiting adverse effects of eribulin. The objective of this analysis was to develop a population pharmacokinetic-pharmacodynamic model for eribulin-associated neutropenia.

METHODS

A combined data set of 12 phase I, II and III studies for eribulin mesilate was analysed. The population pharmacokinetics of eribulin was described using a previously developed model. The relationship between eribulin pharmacokinetic and neutropenia was described using a semi-physiological lifespan model for haematological toxicity. Patient characteristics predictive of increased sensitivity to develop neutropenia were evaluated using a simulation framework.

RESULTS

Absolute neutrophil counts were available from 1579 patients. In the final covariate model, the baseline neutrophil count (ANC0) was estimated to be 4.03 × 10(9) neutrophils l(-1) [relative standard error (RSE) 1.2%], with interindividual variability (IIV, 37.3 coefficient of variation % [CV%]). The mean transition time was estimated to be 109 h (RSE 1.8%, IIV 13.9CV%), the feedback constant (γ) was estimated to be 0.216 (RSE 1.4%, IIV 12.2CV%), and the linear drug effect coefficient (SLOPE) was estimated to be 0.0451 μg l(-1) (RSE 3.2%, IIV 54CV%). Albumin, aspartate transaminase and receival of granulocyte colony-stimulating factor (G-CSF) were identified as significant covariates on SLOPE, and albumin, bilirubin, G-CSF, alkaline phosphatase and lactate dehydrogenase were identified as significant covariates on mean transition time.

CONCLUSIONS

The developed model can be applied to investigate optimal treatment strategies quantitatively across different patient groups with respect to neutropenia. Albumin was identified as the most clinically important covariate predictive of interindividual variability in the neutropenia time course.

The status of pharmacometrics in pregnancy: highlights from the 3(rd) American conference on pharmacometrics

Physiological changes during pregnancy may alter drug pharmacokinetics. Therefore, mechanistic understanding of these changes and, ultimately, clinical studies in pregnant women are necessary to determine if and how dosing regimens should be adjusted. Because of the typically limited number of patients who can be recruited in this patient group, efficient design and analysis of these studies is of special relevance. This paper is a summary of a conference session organized at the American Conference of Pharmacometrics in April 2011, around the topic of applying pharmacometric methodology to this important problem. The discussion included both design and analysis of clinical studies during pregnancy and in silico predictions. An overview of different pharmacometric methods relevant to this subject was given. The impact of pharmacometrics was illustrated using a range of case examples of studies around pregnancy.

Model-based evaluation and optimization of cardiac monitoring protocols for adjuvant treatment of breast cancer with trastuzumab

PURPOSE

Trastuzumab treatment is associated with occurrence of cardiac toxicity, for which monitoring of the left ventricular ejection fraction (LVEF) is indicated. The performance of the currently used monitoring protocol as defined in the summary of product characteristics (SPC) is however unknown. The objective of this analysis was to develop a model-based framework for evaluation and optimization of cardiac monitoring strategies.

METHODS

The model-based framework comprised a previously developed exposure-response model for trastuzumab induced changes in LVEF, and a protocol-execution model that allowed incorporation of treatment interventions as described by a monitoring protocol. Metrics for evaluation of toxicity, dose intensity and monitoring burden were defined to allow evaluation and optimization of cardiac monitoring protocols.

RESULTS

The success of a protocol-defined dose reduction was improved from 40% for the SPC-based protocol, to 79% for a scoring-based protocol, thereby decreasing the observed severity of cardiotoxicity. Including adaptation based on risk-profile allowed reduction of the mean number of LVEF measurements by 19%.

CONCLUSIONS

This model-based evaluation approach enabled evaluation and optimization of cardiac monitoring protocols that would be difficult to evaluate in a clinical setting. This approach can potentially be applied for other drugs that use repeated evaluation of continuous biomarkers for toxicity.

Leveraging physiological data from literature into a pharmacokinetic model to support informative clinical study design in pregnant women

PURPOSE

Physiological changes during pregnancy can effect pharmacokinetic (PK) parameters, which may lead to altered dose requirements. We aimed to leverage literature-based physiological changes during pregnancy into a PK model and compare its performance to a published reference model in pregnant women and to use the literature-based model to determine informative PK sampling times for a clinical study that aims to quantify the PK of enoxaparin throughout pregnancy.

METHODS

Changes in total body water (BW) and creatinine clearance (CRCL) during pregnancy were described using regression models. BW and CRCL were linked to a PK model of enoxaparin in non-pregnant women. Performance of the literature-based PK model was compared to a previously published empirical reference model. D-optimal sampling times were determined and evaluated for literature-based and reference models.

RESULTS

The literature-based model adequately predicted anti-Xa plasma concentrations when compared to reference model predictions. An informative sampling design was successfully developed, with parameters expected with good precision (RSE < 36.4%).

CONCLUSION

A literature-based model describing enoxaparin PK during pregnancy was developed and evaluated. The modelling framework could be used to support development of informative designs in pregnancy when prior models are unavailable.

Pharmacokinetic-pharmacodynamic relationships of central nervous system effects of scopolamine in healthy subjects

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT

• The cholinergic system is important for different central nervous system functions, including memory, learning and attention. Scopolamine, a centrally active muscarinic antagonist, has been used to model dementia and to demonstrate the pharmacological effects of cholinergic drugs, but for most effects the concentration-effect relationships are unknown.

WHAT THIS STUDY ADDS

• We determined the pharmacokinetic-pharmacodynamic relationships of scopolamine using a multidimensional central nervous system test battery in a large group of healthy volunteers. The results suggested there are various functional cholinergic systems with different pharmacological characteristics, which can be used to study the effects of drugs that directly or indirectly modify cholinergic systems. The design of such studies should take the different concentration-effect relationships into account. AIM(S) Although scopolamine is a frequently used memory impairment model, the relationships between exposure and corresponding central nervous system (CNS) effects are mostly unknown. The aim of our study was to characterize these using pharmacokinetic-pharmacodynamic (PK-PD) modelling.

METHODS

In two double-blind, placebo-controlled, four-way crossover studies, 0.5-mg scopolamine was administered i.v. to 90 healthy male subjects. PK and PD/safety measures were monitored pre-dose and up to 8.5 h after administration. PK-PD relationships were modelled using non-linear mixed-effect modelling.

RESULTS

Most PD responses following scopolamine administration in 85 subjects differed significantly from placebo. As PD measures lagged behind the plasma PK profile, PK-PD relationships were modelled using an effect compartment and arbitrarily categorized according to their equilibration half-lives (t(1/2) k(eo) ; hysteresis measure). t(1/2) k(eo) for heart rate was 17 min, saccadic eye movements and adaptive tracking 1-1.5 h, body sway, smooth pursuit, visual analogue scales alertness and psychedelic 2.5-3.5 h, pupil size, finger tapping and visual analogue scales feeling high more than 8 h.

CONCLUSIONS

Scopolamine affected different CNS functions in a concentration-dependent manner, which based on their distinct PK-PD characteristics seemed to reflect multiple distinct functional pathways of the cholinergic system. All PD effects showed considerable albeit variable delays compared with plasma concentrations. The t(1/2) k(eo) of the central effects was longer than of the peripheral effects on heart rate, which at least partly reflects the long CNS retention of scopolamine, but possibly also the triggering of independent secondary mechanisms. PK-PD analysis can optimize scopolamine administration regimens for future research and give insight into the physiology and pharmacology of human cholinergic systems.