Leveraging QSP Models for MIPD: A Case Study for Warfarin/INR

Warfarin dosing remains challenging due to substantial inter-individual variability, which can lead to unsafe or ineffective therapy with standard dosing. Model-informed precision dosing (MIPD) can help individualize warfarin dosing, requiring the selection of a suitable model. For models developed from clinical data, the dependence on the study design and population raises questions about generalizability. Quantitative system pharmacology (QSP) models promise better extrapolation abilities; however, their complexity and lack of validation on clinical data raise questions about applicability in MIPD. We have previously derived a mechanistic warfarin/international normalized ratio (INR) model from a blood coagulation QSP model. In this article, we evaluated the predictive performance of the warfarin/INR model in the context of MIPD using an external dataset with INR data from patients starting warfarin treatment. We assessed the accuracy and precision of model predictions, benchmarked against an empirically based reference model. Additionally, we evaluated covariate contributions and assessed the predictive performance separately in the more challenging outpatient data. The warfarin/INR model performed comparably to the reference model across various measures despite not being calibrated with warfarin initiation data. Including CYP2C9 and/or VKORC1 genotypes as covariates improved the prediction quality of the warfarin/INR model, even after assimilating 4 days of INR data. The outpatient INR exhibited higher unexplained variability, and predictions slightly exceeded observed values, suggesting that model adjustments might be necessary when transitioning from an inpatient to an outpatient setting. Overall, this research underscores the potential of QSP-derived models for MIPD, offering a complementary approach to empirical model development.

TLD-1, a novel liposomal doxorubicin, in patients with advanced solid tumors: Dose escalation and expansion part of a multicenter open-label phase I trial (SAKK 65/16)

BACKGROUND

TLD-1 is a novel liposomal doxorubicin that compared favorably to conventional doxorubicin liposomal formulations in preclinical models. This phase I first-in-human study aimed to define the maximum tolerated dose (MTD), recommended phase 2 dose (RP2D), safety and preliminary activity of TLD-1 in patients with advanced solid tumors.

PATIENTS AND METHODS

We recruited patients with advanced solid tumors who failed standard therapy and received up to 3 prior lines of palliative systemic chemotherapy. TLD-1 was administered intravenously every 3 weeks up to a maximum of 9 cycles (6 for patients with prior anthracyclines) from a starting dose of 10 mg/m2, according to an accelerated titration design followed by a modified continual reassessment method.

RESULTS

30 patients were enrolled between November 2018 and May 2021. No dose-limiting toxicities (DLT) were observed. Maximum administered dose of TLD-1 was 45 mg/m2, RP2D was defined at 40 mg/m2. Most frequent treatment-related adverse events (TRAE) of any grade included palmar-plantar erythrodysesthesia (PPE) (50% of patients), oral mucositis (50%), fatigue (30%) and skin rash (26.7%). Most common G3 TRAE included PPE in 4 patients (13.3%) and oral mucositis in 2 (6.7%). Overall objective response rate was 10% in the whole population and 23.1% among 13 patients with breast cancer; median time-to-treatment failure was 2.7 months. TLD-1 exhibit linear pharmacokinetics, with a median terminal half-life of 95 h.

CONCLUSIONS

The new liposomal doxorubicin formulation TLD-1 showed a favourable safety profile and antitumor activity, particularly in breast cancer. RP2D was defined at 40 mg/m2 administered every 3 weeks. (NCT03387917).

Bayesian framework for multi-source data integration-Application to human extrapolation from preclinical studies

In preclinical investigations, for example, in in vitro, in vivo, and in silico studies, the pharmacokinetic, pharmacodynamic, and toxicological characteristics of a drug are evaluated before advancing to first-in-man trial. Usually, each study is analyzed independently and the human dose range does not leverage the knowledge gained from all studies. Taking into account all preclinical data through inferential procedures can be particularly interesting in obtaining a more precise and reliable starting dose and dose range. Our objective is to propose a Bayesian framework for multi-source data integration, customizable, and tailored to the specific research question. We focused on preclinical results extrapolated to humans, which allowed us to predict the quantities of interest (e.g. maximum tolerated dose, etc.) in humans. We build an approach, divided into four steps, based on a sequential parameter estimation for each study, extrapolation to human, commensurability checking between posterior distributions and final information merging to increase the precision of estimation. The new framework is evaluated via an extensive simulation study, based on a real-life example in oncology. Our approach allows us to better use all the information compared to a standard framework, reducing uncertainty in the predictions and potentially leading to a more efficient dose selection.

PBTK model-based analysis of CYP3A4 induction and the toxicokinetics of the pyrrolizidine alkaloid retrorsine in man

Cytochrome P450 (CYP)3A4 induction by drugs and pesticides plays a critical role in the enhancement of pyrrolizidine alkaloid (PA) toxicity as it leads to increased formation of hepatotoxic dehydro-PA metabolites. Addressing the need for a quantitative analysis of this interaction, we developed a physiologically-based toxicokinetic (PBTK) model. Specifically, the model describes the impact of the well-characterized CYP3A4 inducer rifampicin on the kinetics of retrorsine, which is a prototypic PA and contaminant in herbal teas. Based on consumption data, the kinetics after daily intake of retrorsine were simulated with concomitant rifampicin treatment. Strongest impact on retrorsine kinetics (plasma AUC24 and C max reduced to 67% and 74% compared to the rifampicin-free reference) was predicted directly after withdrawal of rifampicin. At this time point, the competitive inhibitory effect of rifampicin stopped, while CYP3A4 induction was still near its maximum. Due to the impacted metabolism kinetics, the cumulative formation of intestinal retrorsine CYP3A4 metabolites increased to 254% (from 10 to 25 nmol), while the cumulative formation of hepatic CYP3A4 metabolites was not affected (57 nmol). Return to baseline PA toxicokinetics was predicted 14 days after stop of a 14-day rifampicin treatment. In conclusion, the PBTK model showed to be a promising tool to assess the dynamic interplay of enzyme induction and toxification pathways.

Nonlinear Mixed-Effects Model of Z-Endoxifen Concentrations in Tamoxifen-Treated Patients from the CEPAM Cohort

Tamoxifen is widely used in patients with hormone receptor-positive breast cancer. The polymorphic enzyme CYP2D6 is primarily responsible for metabolic activation of tamoxifen, resulting in substantial interindividual variability of plasma concentrations of its most important metabolite, Z-endoxifen. The Z-endoxifen concentration thresholds below which tamoxifen treatment is less efficacious have been proposed but not validated, and prospective trials of individualized tamoxifen treatment to achieve Z-endoxifen concentration thresholds are considered infeasible. Therefore, we aim to validate the association between Z-endoxifen concentration and tamoxifen treatment outcomes, and identify a Z-endoxifen concentration threshold of tamoxifen efficacy, using pharmacometric modeling and simulation. As a first step, the CYP2D6 Endoxifen Percentage Activity Model (CEPAM) cohort was created by pooling data from 28 clinical studies (> 7,000 patients) with measured endoxifen plasma concentrations. After cleaning, data from 6,083 patients were used to develop a nonlinear mixed-effect (NLME) model for tamoxifen and Z-endoxifen pharmacokinetics that includes a conversion factor to allow inclusion of studies that measured total endoxifen but not Z-endoxifen. The final parent-metabolite NLME model confirmed the primary role of CYP2D6, and contributions from body weight, CYP2C9 phenotype, and co-medication with CYP2D6 inhibitors, on Z-endoxifen pharmacokinetics. Future work will use the model to simulate Z-endoxifen concentrations in patients receiving single agent tamoxifen treatment within large prospective clinical trials with long-term survival to identify the Z-endoxifen concentration threshold below which tamoxifen is less efficacious. Identification of this concentration threshold would allow personalized tamoxifen treatment to improve outcomes in patients with hormone receptor-positive breast cancer.

Therapeutic drug monitoring of vedolizumab therapy in inflammatory bowel disease

BACKGROUND

Therapeutic drug monitoring is effective for optimizing anti-tumor necrosis factor therapies in inflammatory bowel disease, but for vedolizumab, a gut-selective leucocyte migration inhibitor, data are scarce.

METHODS

Observational cohort study including 116 bio-experienced inflammatory bowel disease patients treated with vedolizumab for active luminal disease. Biobanked trough blood samples (n = 676) covering 96% of patients were analyzed using a drug-binding immunofluorometric assay. Steroid-free treatment outcomes were classified by clinical disease activity indices and objective findings, primarily endoscopy.

RESULTS

Patients with clinical remission to vedolizumab induction therapy (37%) had significantly higher trough levels than those without at weeks 6 (mean 34.1 vs 28.0 μg/mL, P = 0.03) and 10 (34.8 vs 27.5 μg/mL, P = 0.01). Optimal thresholds for discrimination were 32.4 μg/mL (AUCROC 0.66, P = 0.04) and 23.5 (AUCROC 0.67, P = 0.01), respectively. This positive association persisted during maintenance phase with 11.9 μg/mL (AUCROC 0.69, P < 0.01) associated with clinical remission (37%) and 15.3 (AUCROC 0.74, P < 0.001) for objective remission (46%). Stratification by temporal evolution of treatment effects revealed higher induction and maintenance vedolizumab levels in persistent and slow responders as compared to secondary or persistent failures. Pharmacokinetics was influenced by rare formation of anti-vedolizumab antibodies (2%), and to a lesser extent gender and albumin during induction, but not disease severity, concomitant steroids, or thiopurine metabolites. Switching to subcutaneous administrations resulted in 2.3-fold increase in steady-state trough levels.

CONCLUSION

Our study supports maintaining adequate drug exposure being essential for sustained positive outcomes of vedolizumab and emphasizes individualized, therapeutic drug monitoring-based treatment regimens. Controlled trials and pharmacokinetic modeling are, however, needed.

Index analysis: An approach to understand signal transduction with application to the EGFR signalling pathway

In systems biology and pharmacology, large-scale kinetic models are used to study the dynamic response of a system to a specific input or stimulus. While in many applications, a deeper understanding of the input-response behaviour is highly desirable, it is often hindered by the large number of molecular species and the complexity of the interactions. An approach that identifies key molecular species for a given input-response relationship and characterises dynamic properties of states is therefore highly desirable. We introduce the concept of index analysis; it is based on different time- and state-dependent quantities (indices) to identify important dynamic characteristics of molecular species. All indices are defined for a specific pair of input and response variables as well as for a specific magnitude of the input. In application to a large-scale kinetic model of the EGFR signalling cascade, we identified different phases of signal transduction, the peculiar role of Phosphatase3 during signal activation and Ras recycling during signal onset. In addition, we discuss the challenges and pitfalls of interpreting the relevance of molecular species based on knock-out simulation studies, and provide an alternative view on conflicting results on the importance of parallel EGFR downstream pathways. Beyond the applications in model interpretation, index analysis is envisioned to be a valuable tool in model reduction.

Impact of Key Components of Intensified Ceftaroline Dosing on Pharmacokinetic/Pharmacodynamic Target Attainment

BACKGROUND AND OBJECTIVE

Ceftaroline fosamil is a β-lactam antibiotic approved as a 600 mg twice daily dose (≤1 h infusion, 'standard dosing') or a 600 mg thrice daily dose (2 h infusion) to treat complicated skin and soft tissue infections caused by Staphylococcus aureus (minimum inhibitory concentration [MIC] 2-4 mg/L). We sought to systematically evaluate the relative impact of the three key components of the intensified dosing regimen (i.e. shortened dosing interval, prolonged infusion duration and increased total daily dose [TDD]) on the pharmacokinetic/pharmacodynamic (PK/PD) target attainment given different grades of bacterial susceptibility.

METHODS

A population PK model was developed using data from 12 healthy volunteers (EudraCT-2012-005134-11) receiving standard or intensified dosing. PK/PD target attainment (ƒT>MIC = 35% and 100%) after 24 h was compared following systematically varied combinations of the (1) dosing interval (every 12 h [q12h]→ every 8 h [q8h]); (2) infusion duration (1 h→2 h); and (3) individual and total daily dose (400→900 mg, i.e. TDD 1200→1800 mg), as well as for varying susceptibility of S. aureus (MIC 0.032-8 mg/L).

RESULTS

A two-compartment model with linear elimination adequately described ceftaroline concentrations (n = 274). The relevance of the dosing components dosing interval/infusion duration/TDD for ƒT>MIC systematically changed with pathogen susceptibility. For susceptible pathogens with MIC ≤1 mg/L, shortened dosing intervals appeared as the main driver of the improved target attainment associated with the intensified dosing regimen, followed by increased TDD and infusion duration. For less susceptible pathogens, the advantage of q8h dosing and 2 h infusions declined, and increased TDD improved ƒT>MIC the most.

CONCLUSION

The analysis calls to mind consideration of dose increases when prolonging the infusion duration in the case of low bacterial susceptibility.

C-Reactive Protein as an Early Predictor of Efficacy in Advanced Non-Small-Cell Lung Cancer Patients: A Tumor Dynamics-Biomarker Modeling Framework

In oncology, longitudinal biomarkers reflecting the patient's status and disease evolution can offer reliable predictions of the patient's response to treatment and prognosis. By leveraging clinical data in patients with advanced non-small-cell lung cancer receiving first-line chemotherapy, we aimed to develop a framework combining anticancer drug exposure, tumor dynamics (RECIST criteria), and C-reactive protein (CRP) concentrations, using nonlinear mixed-effects models, to evaluate and quantify by means of parametric time-to-event models the significance of early longitudinal predictors of progression-free survival (PFS) and overall survival (OS). Tumor dynamics was characterized by a tumor size (TS) model accounting for anticancer drug exposure and development of drug resistance. CRP concentrations over time were characterized by a turnover model. An x-fold change in TS from baseline linearly affected CRP production. CRP concentration at treatment cycle 3 (day 42) and the difference between CRP concentration at treatment cycles 3 and 2 were the strongest predictors of PFS and OS. Measuring longitudinal CRP allows for the monitoring of inflammatory levels and, along with its reduction across treatment cycles, presents a promising prognostic marker. This framework could be applied to other treatment modalities such as immunotherapies or targeted therapies allowing the timely identification of patients at risk of early progression and/or short survival to spare them unnecessary toxicities and provide alternative treatment decisions.

A prognostic baseline blood biomarker and tumor growth kinetics integrated model in paclitaxel/platinum treated advanced non-small cell lung cancer patients

Paclitaxel/platinum chemotherapy, the backbone of standard first-line treatment of advanced non-small cell lung cancer (NSCLC), exhibits high interpatient variability in treatment response and high toxicity burden. Baseline blood biomarker concentrations and tumor size (sum of diameters) at week 8 relative to baseline (RS8) are widely investigated prognostic factors. However, joint analysis of data on demographic/clinical characteristics, blood biomarker levels, and chemotherapy exposure-driven early tumor response for improved prediction of overall survival (OS) is clinically not established. We developed a Weibull time-to-event model to predict OS, leveraging data from 365 patients receiving paclitaxel/platinum combination chemotherapy once every three weeks for ≤six cycles. A developed tumor growth inhibition model, combining linear tumor growth and first-order paclitaxel area under the concentration-time curve-induced tumor decay, was used to derive individual RS8. The median model-derived RS8 in all patients was a 20.0% tumor size reduction (range from -78% to +15%). Whereas baseline carcinoembryonic antigen, cytokeratin fragments, and thyroid stimulating hormone levels were not significantly associated with OS in a subset of 221 patients, and lactate dehydrogenase, interleukin-6 and neutrophil-to-lymphocyte ratio levels were significant only in univariate analyses (p value < 0.05); C-reactive protein (CRP) in combination with RS8 most significantly affected OS (p value < 0.01). Compared to the median population OS of 11.3 months, OS was 128% longer at the 5th percentile levels of both covariates and 60% shorter at their 95th percentiles levels. The combined paclitaxel exposure-driven RS8 and baseline blood CRP concentrations enables early individual prognostic predictions for different paclitaxel dosing regimens, forming the basis for treatment decision and optimizing paclitaxel/platinum-based advanced NSCLC chemotherapy.

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.

PBTK modeling of the pyrrolizidine alkaloid retrorsine to predict liver toxicity in mouse and rat

Retrorsine is a hepatotoxic pyrrolizidine alkaloid (PA) found in herbal supplements and medicines, food and livestock feed. Dose-response studies enabling the derivation of a point of departure including a benchmark dose for risk assessment of retrorsine in humans and animals are not available. Addressing this need, a physiologically based toxicokinetic (PBTK) model of retrorsine was developed for mouse and rat. Comprehensive characterization of retrorsine toxicokinetics revealed: both the fraction absorbed from the intestine (78%) and the fraction unbound in plasma (60%) are high, hepatic membrane permeation is dominated by active uptake and not by passive diffusion, liver metabolic clearance is 4-fold higher in rat compared to mouse and renal excretion contributes to 20% of the total clearance. The PBTK model was calibrated with kinetic data from available mouse and rat studies using maximum likelihood estimation. PBTK model evaluation showed convincing goodness-of-fit for hepatic retrorsine and retrorsine-derived DNA adducts. Furthermore, the developed model allowed to translate in vitro liver toxicity data of retrorsine to in vivo dose-response data. Resulting benchmark dose confidence intervals (mg/kg bodyweight) are 24.1-88.5 in mice and 79.9-104 in rats for acute liver toxicity after oral retrorsine intake. As the PBTK model was built to enable extrapolation to different species and other PA congeners, this integrative framework constitutes a flexible tool to address gaps in the risk assessment of PA.

Model-Informed Target Morning 17α-Hydroxyprogesterone Concentrations in Dried Blood Spots for Pediatric Congenital Adrenal Hyperplasia Patients

Monitoring cortisol replacement therapy in congenital adrenal hyperplasia (CAH) patients is vital to avoid serious adverse events such as adrenal crises due to cortisol underexposure or metabolic consequences due to cortisol overexposure. The less invasive dried blood spot (DBS) sampling is an advantageous alternative to traditional plasma sampling, especially in pediatric patients. However, target concentrations for important disease biomarkers such as 17α-hydroxyprogesterone (17-OHP) are unknown using DBS. Therefore, a modeling and simulation framework, including a pharmacokinetic/pharmacodynamic model linking plasma cortisol concentrations to DBS 17-OHP concentrations, was used to derive a target morning DBS 17-OHP concentration range of 2-8 nmol/L in pediatric CAH patients. Since either capillary or venous DBS sampling is becoming more common in the clinics, the clinical applicability of this work was shown by demonstrating the comparability of capillary and venous cortisol and 17-OHP concentrations collected by DBS sampling, using a Bland-Altman and Passing-Bablok analysis. The derived target morning DBS 17-OHP concentration range is a first step towards providing improved therapy monitoring using DBS sampling and adjusting hydrocortisone (synthetic cortisol) dosing in children with CAH. In the future, this framework can be used to assess further research questions, e.g., target replacement ranges for the entire day.

Deriving mechanism-based pharmacodynamic models by reducing quantitative systems pharmacology models: An application to warfarin

Quantitative systems pharmacology (QSP) models integrate comprehensive qualitative and quantitative knowledge about pharmacologically relevant processes. We previously proposed a first approach to leverage the knowledge in QSP models to derive simpler, mechanism-based pharmacodynamic (PD) models. Their complexity, however, is typically still too large to be used in the population analysis of clinical data. Here, we extend the approach beyond state reduction to also include the simplification of reaction rates, elimination of reactions, and analytic solutions. We additionally ensure that the reduced model maintains a prespecified approximation quality not only for a reference individual but also for a diverse virtual population. We illustrate the extended approach for the warfarin effect on blood coagulation. Using the model-reduction approach, we derive a novel small-scale warfarin/international normalized ratio model and demonstrate its suitability for biomarker identification. Due to the systematic nature of the approach in comparison with empirical model building, the proposed model-reduction algorithm provides an improved rationale to build PD models also from QSP models in other applications.

Development of a project for interprofessional collaboration between medical and pharmacy students to improve medication safety in polypharmacy (PILLE)

AIM

Interprofessional collaboration is particularly relevant to patient safety in outpatient care with polypharmacy. The educational project "PILLE" is meant to give medical and pharmacy students an understanding of the roles and competencies needed for cooperation in the provision of healthcare and to enable interprofessional learning.

METHOD

The curriculum is aimed at pharmacy and medical students and was developed in six steps according to the Kern cycle. It is comprised of an interprofessional seminar, a joint practical training in a simulated pharmacy, and a tandem job shadowing at a primary care practice. The project was implemented in three stages due to the pandemic: The interprofessional online seminar based on the ICAP model and the digital inverted classroom was held in the 2020 winter semester; the interprofessional practical training was added in the 2021 summer semester; and the interprofessional tandem job shadowing at a primary care practice in the 2021 winter semester. Attitudes toward interprofessional learning, among other things, was measured in the evaluation using the SPICE-2D questionnaire (Student Perceptions of Physician-Pharmacist Interprofessional Clinical Education).

RESULTS

In the first three semesters, a total of 105 students (46 pharmacy, 59 medicine) participated in the project, of which 78 participated in the evaluation (74% response rate). The students stated, in particular, that they had learned about the competencies and roles of the other profession and desired additional and more specific preparatory materials for the course sessions. The SPICE-2D questionnaire showed high values for both groups of students already in the pre-survey and these increased further as a result of the project.

CONCLUSION

Joint case-based learning could be implemented under the conditions imposed by the pandemic. Online teaching is a low-threshold means to enable interprofessional exchange.

Cell-level systems biology model to study inflammatory bowel diseases and their treatment options

To help understand the complex and therapeutically challenging inflammatory bowel diseases (IBDs), we developed a systems biology model of the intestinal immune system that is able to describe main aspects of IBD and different treatment modalities thereof. The model, including key cell types and processes of the mucosal immune response, compiles a large amount of isolated experimental findings from literature into a larger context and allows for simulations of different inflammation scenarios based on the underlying data and assumptions. In the context of a large and diverse virtual IBD population, we characterized the patients based on their phenotype (in contrast to healthy individuals, they developed persistent inflammation after a trigger event) rather than on a priori assumptions on parameter differences to a healthy individual. This allowed to reproduce the enormous diversity of predispositions known to lead to IBD. Analyzing different treatment effects, the model provides insight into characteristics of individual drug therapy. We illustrate for anti-TNF-α therapy, how the model can be used (i) to decide for alternative treatments with best prospects in the case of nonresponse, and (ii) to identify promising combination therapies with other available treatment options.

Insights in the maturational processes influencing hydrocortisone pharmacokinetics in congenital adrenal hyperplasia patients using a middle-out approach

Introduction: Hydrocortisone is the standard of care in cortisol replacement therapy for congenital adrenal hyperplasia patients. Challenges in mimicking cortisol circadian rhythm and dosing individualization can be overcome by the support of mathematical modelling. Previously, a non-linear mixed-effects (NLME) model was developed based on clinical hydrocortisone pharmacokinetic (PK) pediatric and adult data. Additionally, a physiologically-based pharmacokinetic (PBPK) model was developed for adults and a pediatric model was obtained using maturation functions for relevant processes. In this work, a middle-out approach was applied. The aim was to investigate whether PBPK-derived maturation functions could provide a better description of hydrocortisone PK inter-individual variability when implemented in the NLME framework, with the goal of providing better individual predictions towards precision dosing at the patient level. Methods: Hydrocortisone PK data from 24 adrenal insufficiency pediatric patients and 30 adult healthy volunteers were used for NLME model development, while the PBPK model and maturation functions of clearance and cortisol binding globulin (CBG) were developed based on previous studies published in the literature. Results: Clearance (CL) estimates from both approaches were similar for children older than 1 year (CL/F increasing from around 150 L/h to 500 L/h), while CBG concentrations differed across the whole age range (CBG_NLME stable around 0.5 μM vs. steady increase from 0.35 to 0.8 μM for CBG _PBPK). PBPK-derived maturation functions were subsequently included in the NLME model. After inclusion of the maturation functions, none, a part of, or all parameters were re-estimated. However, the inclusion of CL and/or CBG maturation functions in the NLME model did not result in improved model performance for the CL maturation function (ΔOFV > -15.36) and the re-estimation of parameters using the CBG maturation function most often led to unstable models or individual CL prediction bias. Discussion: Three explanations for the observed discrepancies could be postulated, i) non-considered maturation of processes such as absorption or first-pass effect, ii) lack of patients between 1 and 12 months, iii) lack of correction of PBPK CL maturation functions derived from urinary concentration ratio data for the renal function relative to adults. These should be investigated in the future to determine how NLME and PBPK methods can work towards deriving insights into pediatric hydrocortisone PK.

Microdialysis of Drug and Drug Metabolite: a Comprehensive In Vitro Analysis for Voriconazole and Voriconazole N-oxide

Purpose

Voriconazole is a therapeutically challenging antifungal drug associated with high interindividual pharmacokinetic variability. As a prerequisite to performing clinical trials using the minimally-invasive sampling technique microdialysis, a comprehensive in vitro microdialysis characterization of voriconazole (VRC) and its potentially toxic N-oxide metabolite (NO) was performed.

Methods

The feasibility of simultaneous microdialysis of VRC and NO was explored in vitro by investigating the relative recovery (RR) of both compounds in the absence and presence of the other. The dependency of RR on compound combination, concentration, microdialysis catheter and study day was evaluated and quantified by linear mixed-effects modeling.

Results

Median RR of VRC and NO during individual microdialysis were high (87.6% and 91.1%). During simultaneous microdialysis of VRC and NO, median RR did not change (87.9% and 91.1%). The linear mixed-effects model confirmed the absence of significant differences between RR of VRC and NO during individual and simultaneous microdialysis as well as between the two compounds (p > 0.05). No concentration dependency of RR was found (p = 0.284). The study day was the main source of variability (46.3%) while the microdialysis catheter only had a minor effect (4.33%). VRC retrodialysis proved feasible as catheter calibration for both compounds.

Conclusion

These in vitro microdialysis results encourage the application of microdialysis in clinical trials to assess target-site concentrations of VRC and NO. This can support the generation of a coherent understanding of VRC pharmacokinetics and its sources of variability. Ultimately, a better understanding of human VRC pharmacokinetics might contribute to the development of personalized dosing strategies.

Mutations in ALK signaling pathways conferring resistance to ALK inhibitor treatment lead to collateral vulnerabilities in neuroblastoma cells

BACKGROUND

Development of resistance to targeted therapies has tempered initial optimism that precision oncology would improve poor outcomes for cancer patients. Resistance mechanisms, however, can also confer new resistance-specific vulnerabilities, termed collateral sensitivities. Here we investigated anaplastic lymphoma kinase (ALK) inhibitor resistance in neuroblastoma, a childhood cancer frequently affected by activating ALK alterations.

METHODS

Genome-wide forward genetic CRISPR-Cas9 based screens were performed to identify genes associated with ALK inhibitor resistance in neuroblastoma cell lines. Furthermore, the neuroblastoma cell line NBLW-R was rendered resistant by continuous exposure to ALK inhibitors. Genes identified to be associated with ALK inhibitor resistance were further investigated by generating suitable cell line models. In addition, tumor and liquid biopsy samples of four patients with ALK-mutated neuroblastomas before ALK inhibitor treatment and during tumor progression under treatment were genomically profiled.

RESULTS

Both genome-wide CRISPR-Cas9-based screens and preclinical spontaneous ALKi resistance models identified NF1 loss and activating NRASQ61K mutations to confer resistance to chemically diverse ALKi. Moreover, human neuroblastomas recurrently developed de novo loss of NF1 and activating RAS mutations after ALKi treatment, leading to therapy resistance. Pathway-specific perturbations confirmed that NF1 loss and activating RAS mutations lead to RAS-MAPK signaling even in the presence of ALKi. Intriguingly, NF1 loss rendered neuroblastoma cells hypersensitive to MEK inhibition.

CONCLUSIONS

Our results provide a clinically relevant mechanistic model of ALKi resistance in neuroblastoma and highlight new clinically actionable collateral sensitivities in resistant cells.

Evaluation of a Meropenem and Piperacillin Monitoring Program in Intensive Care Unit Patients Calls for the Regular Assessment of Empirical Targets and Easy-to-Use Dosing Decision Tools

The drug concentrations targeted in meropenem and piperacillin/tazobactam therapy also depend on the susceptibility of the pathogen. Yet, the pathogen is often unknown, and antibiotic therapy is guided by empirical targets. To reliably achieve the targeted concentrations, dosing needs to be adjusted for renal function. We aimed to evaluate a meropenem and piperacillin/tazobactam monitoring program in intensive care unit (ICU) patients by assessing (i) the adequacy of locally selected empirical targets, (ii) if dosing is adequately adjusted for renal function and individual target, and (iii) if dosing is adjusted in target attainment (TA) failure. In a prospective, observational clinical trial of drug concentrations, relevant patient characteristics and microbiological data (pathogen, minimum inhibitory concentration (MIC)) for patients receiving meropenem or piperacillin/tazobactam treatment were collected. If the MIC value was available, a target range of 1–5 × MIC was selected for minimum drug concentrations of both drugs. If the MIC value was not available, 8–40 mg/L and 16–80 mg/L were selected as empirical target ranges for meropenem and piperacillin, respectively. A total of 356 meropenem and 216 piperacillin samples were collected from 108 and 96 ICU patients, respectively. The vast majority of observed MIC values was lower than the empirical target (meropenem: 90.0%, piperacillin: 93.9%), suggesting empirical target value reductions. TA was found to be low (meropenem: 35.7%, piperacillin 50.5%) with the lowest TA for severely impaired renal function (meropenem: 13.9%, piperacillin: 29.2%), and observed drug concentrations did not significantly differ between patients with different targets, indicating dosing was not adequately adjusted for renal function or target. Dosing adjustments were rare for both drugs (meropenem: 6.13%, piperacillin: 4.78%) and for meropenem irrespective of TA, revealing that concentration monitoring alone was insufficient to guide dosing adjustment. Empirical targets should regularly be assessed and adjusted based on local susceptibility data. To improve TA, scientific knowledge should be translated into easy-to-use dosing strategies guiding antibiotic dosing.

TLD-1, a novel liposomal doxorubicin, in patients (pts) with advanced solid tumors: Dose escalation and expansion part of a multicenter open-label phase I trial (SAKK 65/16)

3027

Background: TLD-1 is a novel liposomal doxorubicin that compared favorably to conventional liposomal formulations of doxorubicin in preclinical in vivo mouse breast cancer models. This phase I first-in-human trial is aiming to determine the recommended phase II dose (RP2D), toxicity profile, pharmacokinetics and preliminary activity. Methods: Patients with a maximum of 3 prior lines of systemic chemotherapy and preferably anthracycline-sensitive disease were eligible. TLD-1 was administered on day 1 iv over 60-90 minutes (depending on individual dose) q 21 days, for up to 6 or 9 cycles (according to prior anthracycline-exposure) with premedication of 8mg dexamethasone. Dose escalation with dose levels (DL) 1-7 of 10, 16, 23, 30, 35, 40 and 45mg/m2 started with an accelerated titration design, treating one pt at each DL up to DL6 (40mg/m2) followed by a modified continual reassessment method at DL7 due to observed toxicity. Results: 30 pts (F:M = 24:6) have been treated, one each at DLs 1-6, 15 pts at DL7 and an additional 9 pts at DL6. Most frequent tumor types included breast (n = 13), ovarian (n = 6), cervical cancer (n = 2) and cholangiocarcinoma (n = 2). Median age was 67.5 years (range:38-83), 13 pts were exposed to prior anthracyclines. The median number of cycles was 4 (range:1-9). No dose-limiting toxicities (DLT) occurred during cycle 1. At DLs 1 to 5, no treatment-related G3 AEs (TRAE) were observed. At DL6, there was one case of mucositis G3, one of palmar-plantar-erythrodysesthesia (PPE) G3 and one of anemia and neutropenia G3 each. One patient with pre-existing valvular cardiopathy developed symptoms of heart-failure G3 after 8 cycles. Echocardiography showed severe mitral regurgitation with normal LV-EF. In addition one case of urinary-tract infection G3 was seen. Dose-modifications or -delays due to AEs occurred in 7/50 cycles. At DL7, one case of mucositis G3, 3 events of PPE G3 and one case of fatigue G3 were reported. In addition, one case of infection with shingles occurred. Dose-modifications or -delays due to AEs occurred in 12/61 cycles. Shingles and heart failure were reported as SAEs. All toxicities listed above were categorized as TRAE. 29/30 pts were evaluable for response. Three breast cancer pts had a partial response, 2 at DL7 and 1 at DL6, 14 pts had stable disease. Conclusions: No DLT was observed up to DL7. RP2D was defined at 40mg/m2 due to cumulative PPE G3 at DL7. The trial is ongoing with a comparative PK-part evaluating the two iv liposomal formulations of doxorubicin TLD-1 and Caelyx. Clinical trial information: NCT03387917.

Does the cytokine adsorber CytoSorb® reduce vancomycin exposure in critically ill patients with sepsis or septic shock? a prospective observational study

BACKGROUND

Hemadsorption of cytokines is used in critically ill patients with sepsis or septic shock. Concerns have been raised that the cytokine adsorber CytoSorb^® unintentionally adsorbs vancomycin. This study aimed to quantify vancomycin elimination by CytoSorb^®.

METHODS

Critically ill patients with sepsis or septic shock receiving continuous renal replacement therapy and CytoSorb^® treatment during a prospective observational study were included in the analysis. Vancomycin pharmacokinetics was characterized using population pharmacokinetic modeling. Adsorption of vancomycin by the CytoSorb^® was investigated as linear or saturable process. The final model was used to derive dosing recommendations based on stochastic simulations.

RESULTS

20 CytoSorb^® treatments in 7 patients (160 serum samples/24 during CytoSorb^®-treatment, all continuous infusion) were included in the study. A classical one-compartment model, including effluent flow rate of the continuous hemodialysis as linear covariate on clearance, best described the measured concentrations (without CytoSorb^®). Significant adsorption with a linear decrease during CytoSorb^® treatment was identified (p < 0.0001) and revealed a maximum increase in vancomycin clearance of 291% (initially after CytoSorb^® installation) and a maximum adsorption capacity of 572 mg. For a representative patient of our cohort a reduction of the area under the curve (AUC) by 93 mg/L*24 h during CytoSorb^® treatment was observed. The additional administration of 500 mg vancomycin over 2 h during CytoSorb^® attenuated the effect and revealed a negligible reduction of the AUC by 4 mg/L*24 h.

CONCLUSION

We recommend the infusion of 500 mg vancomycin over 2 h during CytoSorb^® treatment to avoid subtherapeutic concentrations. Trial registration NCT03985605. Registered 14 June 2019, https://clinicaltrials.gov/ct2/show/NCT03985605.

Exploring Dried Blood Spot Cortisol Concentrations as an Alternative for Monitoring Pediatric Adrenal Insufficiency Patients: A Model-Based Analysis

Congenital adrenal hyperplasia (CAH) is the most common form of adrenal insufficiency in childhood; it requires cortisol replacement therapy with hydrocortisone (HC, synthetic cortisol) from birth and therapy monitoring for successful treatment. In children, the less invasive dried blood spot (DBS) sampling with whole blood including red blood cells (RBCs) provides an advantageous alternative to plasma sampling. Potential differences in binding/association processes between plasma and DBS however need to be considered to correctly interpret DBS measurements for therapy monitoring. While capillary DBS samples would be used in clinical practice, venous cortisol DBS samples from children with adrenal insufficiency were analyzed due to data availability and to directly compare and thus understand potential differences between venous DBS and plasma. A previously published HC plasma pharmacokinetic (PK) model was extended by leveraging these DBS concentrations. In addition to previously characterized binding of cortisol to albumin (linear process) and corticosteroid-binding globulin (CBG; saturable process), DBS data enabled the characterization of a linear cortisol association with RBCs, and thereby providing a quantitative link between DBS and plasma cortisol concentrations. The ratio between the observed cortisol plasma and DBS concentrations varies highly from 2 to 8. Deterministic simulations of the different cortisol binding/association fractions demonstrated that with higher blood cortisol concentrations, saturation of cortisol binding to CBG was observed, leading to an increase in all other cortisol binding fractions. In conclusion, a mathematical PK model was developed which links DBS measurements to plasma exposure and thus allows for quantitative interpretation of measurements of DBS samples.

Plasma concentration guided dosing of drugs used for the treatment of childhood leukaemias: protocol for a systematic review

INTRODUCTION

Childhood leukaemia is the most common type of cancer in children and represents among 25% of the diagnoses in children <15 years old. Childhood survival rates have significantly improved within the last 40 years due to a rapid advancement in therapeutic interventions. However, in high-risk groups, survival rates remain poor. Pharmacokinetic (PK) data of cancer medications in children are limited and thus current dosing regimens are based on studies with small sample sizes. In adults, large variability in PK is observed and dose individualisation (plasma concentration guided dosing) has been associated with improved clinical outcomes; whether this is true for children is still unknown. This provides an opportunity to explore this strategy in children to potentially reduce toxicities and ensure optimal dosing. This paper will provide a protocol to systematically review studies that have used dose individualisation of drugs used in the treatment of childhood leukaemias.

METHODS AND ANALYSIS

Systematic review methodology will be applied to identify, select and extract data from published plasma guided dosing studies conducted in a paediatric leukaemia cohort. Databases (eg, Ovid Embase, Ovid MEDLINE, Ovid Cochrane) and clinical trial registries (CENTRAL, ClinicalTrials.gov and ISRCTN) will be used to perform the systematic literature search (up until February 2021). Only full empirical studies will be included, with primary clinical outcomes (progression-free survival, toxicities, minimal residual disease status, complete cytogenetic response, partial cytogenetic response and major molecular response) being used to decide whether the study will be included. The quality of included studies will be undertaken, with a subgroup analysis where appropriate.

ETHICS AND DISSEMINATION

This systematic review will not require ethics approval as there will not be collection of primary data. Findings of this review will be made available through publications in peer-reviewed journals and conference presentations. Gaps will be identified in current literature to inform future-related research.

PROSPERO REGISTRATION NUMBER

CRD42021225045.

Microdialysis of Voriconazole and its N-Oxide Metabolite: Amalgamating Knowledge of Distribution and Metabolism Processes in Humans

PURPOSE

Voriconazole is an essential antifungal drug whose complex pharmacokinetics with high interindividual variability impedes effective and safe therapy. By application of the minimally-invasive sampling technique microdialysis, interstitial space fluid (ISF) concentrations of VRC and its potentially toxic N-oxide metabolite (NO) were assessed to evaluate target-site exposure for further elucidating VRC pharmacokinetics.

METHODS

Plasma and ISF samples of a clinical trial with an approved VRC dosing regimen were analyzed for VRC and NO concentrations. Concentration-time profiles, exposure assessed as area-under-the-curve (AUC) and metabolic ratios of four healthy adults in plasma and ISF were evaluated regarding the impact of multiple dosing and CYP2C19 genotype.

RESULTS

VRC and NO revealed distribution into ISF with AUC values being ≤2.82- and 17.7-fold lower compared to plasma, respectively. Intraindividual variability of metabolic ratios was largest after the first VRC dose administration while interindividual variability increased with multiple dosing. The CYP2C19 genotype influenced interindividual differences with a maximum 6- and 24-fold larger AUC_NO/AUC_VRC ratio between the intermediate and rapid metabolizer in plasma and ISF, respectively. VRC metabolism was saturated/auto-inhibited indicated by substantially decreasing metabolic concentration ratios with increasing VRC concentrations and after multiple dosing.

CONCLUSION

The feasibility of the simultaneous microdialysis of VRC and NO in vivo was demonstrated and provided new quantitative insights by leveraging distribution and metabolism processes of VRC in humans. The exploratory analysis suggested substantial dissimilarities of VRC and NO pharmacokinetics in plasma and ISF. Ultimately, a thorough understanding of target-site pharmacokinetics might contribute to the optimization of personalized VRC dosing regimens.

Quantification of the Time Course of CYP3A Inhibition, Activation, and Induction Using a Population Pharmacokinetic Model of Microdosed Midazolam Continuous Infusion

BACKGROUND

Cytochrome P450 (CYP) 3A contributes to the metabolism of many approved drugs. CYP3A perpetrator drugs can profoundly alter the exposure of CYP3A substrates. However, effects of such drug-drug interactions are usually reported as maximum effects rather than studied as time-dependent processes. Identification of the time course of CYP3A modulation can provide insight into when significant changes to CYP3A activity occurs, help better design drug-drug interaction studies, and manage drug-drug interactions in clinical practice.

OBJECTIVE

We aimed to quantify the time course and extent of the in vivo modulation of different CYP3A perpetrator drugs on hepatic CYP3A activity and distinguish different modulatory mechanisms by their time of onset, using pharmacologically inactive intravenous microgram doses of the CYP3A-specific substrate midazolam, as a marker of CYP3A activity.

METHODS

Twenty-four healthy individuals received an intravenous midazolam bolus followed by a continuous infusion for 10 or 36 h. Individuals were randomized into four arms: within each arm, two individuals served as a placebo control and, 2 h after start of the midazolam infusion, four individuals received the CYP3A perpetrator drug: voriconazole (inhibitor, orally or intravenously), rifampicin (inducer, orally), or efavirenz (activator, orally). After midazolam bolus administration, blood samples were taken every hour (rifampicin arm) or every 15 min (remaining study arms) until the end of midazolam infusion. A total of 1858 concentrations were equally divided between midazolam and its metabolite, 1'-hydroxymidazolam. A nonlinear mixed-effects population pharmacokinetic model of both compounds was developed using NONMEM^®. CYP3A activity modulation was quantified over time, as the relative change of midazolam clearance encountered by the perpetrator drug, compared to the corresponding clearance value in the placebo arm.

RESULTS

Time course of CYP3A modulation and magnitude of maximum effect were identified for each perpetrator drug. While efavirenz CYP3A activation was relatively fast and short, reaching a maximum after approximately 2-3 h, the induction effect of rifampicin could only be observed after 22 h, with a maximum after approximately 28-30 h followed by a steep drop to almost baseline within 1-2 h. In contrast, the inhibitory impact of both oral and intravenous voriconazole was prolonged with a steady inhibition of CYP3A activity followed by a gradual increase in the inhibitory effect until the end of sampling at 8 h. Relative maximum clearance changes were +59.1%, +46.7%, -70.6%, and -61.1% for efavirenz, rifampicin, oral voriconazole, and intravenous voriconazole, respectively.

CONCLUSIONS

We could distinguish between different mechanisms of CYP3A modulation by the time of onset. Identification of the time at which clearance significantly changes, per perpetrator drug, can guide the design of an optimal sampling schedule for future drug-drug interaction studies. The impact of a short-term combination of different perpetrator drugs on the paradigm CYP3A substrate midazolam was characterized and can define combination intervals in which no relevant interaction is to be expected.

CLINICAL TRIAL REGISTRATION

The trial was registered at the European Union Drug Regulating Authorities for Clinical Trials (EudraCT-No. 2013-004869-14).

A versatile high-performance LC-MS/MS assay for the quantification of voriconazole and its N-oxide metabolite in small sample volumes of multiple human matrices for biomedical applications

Voriconazole (VRC) pharmacokinetics, in particular its complex metabolism, is still not fully understood which challenges its optimal therapeutic use. To increase knowledge on the pharmacokinetics of this antifungal drug, it is essential to broaden the perspective and expand in vitro and clinical in vivo investigations in particular to aspects such as unbound plasma, target-site and metabolite concentrations. Innovative sampling approaches such as microdialysis, a minimally-invasive technique for the analysis of compound concentrations in target-site human tissue fluids, are associated with bioanalytical challenges, i.e. small sample volumes and low concentrations. Thus, a bioanalytical LC-MS/MS assay for the simultaneous quantification of VRC and its main N-oxide (NO) metabolite in human plasma, ultrafiltrate and microdialysate was developed and validated according to the European Medicines Agency guideline. Quantification was rapid, simple and feasible for clinically relevant concentrations from 5 to 5000 ng/mL in plasma and ultrafiltrate as well as from 4 to 4000 ng/mL in microdialysate. Due to the high sensitivity of the assay, only 20 µL of plasma or ultrafiltrate and 5 µL of microdialysate were required. For VRC and NO in all matrices, between-run accuracy was high with a maximum mean deviation of 7.0% from the nominal value and between-run precision was demonstrated by ≤ 11.8% coefficient of variation. Both compounds proved stable under various conditions. The assay suitability was demonstrated by the application to a clinical study quantifying simultaneously VRC and NO concentrations in plasma, ultrafiltrate and microdialysate. Additionally, the assay was successfully adapted for pharmacokinetic analyses in human tissue-derived in vitro experiments. Overall, by reducing the required sample volume, the bioanalytical method allows for an increased number of plasma samples in vulnerable populations, e.g. infants, and enables the generation of concentration-time profiles with a higher temporal resolution in microdialysis studies. Consequently, the developed assay is apt to elucidate the complex pharmacokinetics of VRC in clinical settings as prerequisite for therapy optimisation.

Comparative Plasma and Interstitial Tissue Fluid Pharmacokinetics of Meropenem Demonstrate the Need for Increasing Dose and Infusion Duration in Obese and Non-obese Patients

BACKGROUND AND OBJECTIVES

A quantitative evaluation of the PK of meropenem, a broad-spectrum β-lactam antibiotic, in plasma and interstitial space fluid (ISF) of subcutaneous adipose tissue of obese patients is lacking as of date. The objective of this study was the characterisation of meropenem population pharmacokinetics in plasma and ISF in obese and non-obese patients for identification of adequate dosing regimens via Monte-Carlo simulations.

METHODS

We obtained plasma and microdialysate concentrations after administration of meropenem 1000 mg to 15 obese and 15 non-obese surgery patients from a prospective clinical trial. After characterizing plasma- and microdialysis-derived ISF pharmacokinetics via population pharmacokinetic analysis, we simulated thrice-daily (TID) meropenem short-term (0.5 h), prolonged (3.0 h), and continuous infusions. Adequacy of therapy was assessed by the probability of pharmacokinetic/pharmacodynamic (PK/PD) target attainment (PTA) analysis based on time unbound concentrations exceeded minimum inhibitory concentrations (MIC) on treatment day 1 (%fT _{> MIC}) and the sum of PTA weighted by relative frequency of MIC values for infections by pathogens commonly treated with meropenem. To avoid interstitial tissue fluid concentrations below MIC for the entire dosing interval during continuous infusions, a more conservative PK/PD index was selected (%fT _{> 4 × MIC}).

RESULTS

Adjusted body weight (ABW) and calculated creatinine clearance (CLCR_{CG_ABW}) of all patients (body mass index [BMI] = 20.5-81.5 kg/m²) explained a considerable proportion of the between-patient pharmacokinetic variability (15.1-31.0% relative reduction). The ISF:plasma ratio of %fT _{> MIC} was relatively similar for MIC ≤ 2 mg/L but decreased for MIC = 8 mg/L over ABW = 60-120 kg (0.50-0.20). Steady-state concentrations were 2.68 times (95% confidence interval [CI] = 2.11-3.37) higher in plasma than in ISF, supporting PK/PD targets related to four times the MIC during continuous infusions to avoid suspected ISF concentrations constantly below the MIC. A 3000 mg/24 h continuous infusion was sufficient at MIC = 2 mg/L for patients with CLCR_{CG_ABW} ≤ 100 mL/min and ABW < 90 kg, whereas 2000 mg TID prolonged infusions were adequate for those with CLCR_{CG_ABW} ≤ 100 mL/min and ABW > 90 kg. For MIC = 2 mg/L and %fT_{> MIC} = 95, PTA was adequate in patients over the entire investigated range of body mass and renal function using a 6000 mg continuous infusion. A prolonged infusion of meropenem 2000 mg TID was sufficient for MIC ≤ 8 mg/L and all investigated ABW and CLCR_{CG_ABW} when employing the PK/PD target %fT _{> MIC} = 40. Short-term infusions of 1000 mg TID were sufficient for CLCR_{CG_ABW} ≤ 130 mL/min and distributions of MIC values for Escherichia coli, Citrobacter freundii, and Klebsiella pneumoniae but not for Pseudomonas aeruginosa.

CONCLUSIONS

This analysis indicated a need for higher doses (≥ 2000 mg) and prolonged infusions (≥ 3 h) for obese and non-obese patients at MIC ≥ 2 mg/L. Higher PTA was achieved with prolonged infusions in obese patients and with continuous infusions in non-obese patients.

TRIAL REGISTRATION

EudraCT: 2012-004383-22.

Development of a Model-Informed Dosing Tool to Optimise Initial Antibiotic Dosing-A Translational Example for Intensive Care Units

The prevalence and mortality rates of severe infections are high in intensive care units (ICUs). At the same time, the high pharmacokinetic variability observed in ICU patients increases the risk of inadequate antibiotic drug exposure. Therefore, dosing tailored to specific patient characteristics has a high potential to improve outcomes in this vulnerable patient population. This study aimed to develop a tabular dosing decision tool for initial therapy of meropenem integrating hospital-specific, thus far unexploited pathogen susceptibility information. An appropriate meropenem pharmacokinetic model was selected from the literature and evaluated using clinical data. Probability of target attainment (PTA) analysis was conducted for clinically interesting dosing regimens. To inform dosing prior to pathogen identification, the local pathogen-independent mean fraction of response (LPIFR) was calculated based on the observed minimum inhibitory concentrations distribution in the hospital. A simple, tabular, model-informed dosing decision tool was developed for initial meropenem therapy. Dosing recommendations achieving PTA > 90% or LPIFR > 90% for patients with different creatinine clearances were integrated. Based on the experiences during the development process, a generalised workflow for the development of tabular dosing decision tools was derived. The proposed workflow can support the development of model-informed dosing tools for initial therapy of various drugs and hospital-specific conditions.

EGb 761® Does Not Affect Blood Coagulation and Bleeding Time in Patients with Probable Alzheimer's Dementia-Secondary Analysis of a Randomized, Double-Blind Placebo-Controlled Trial

Following reports of bleeding upon Ginkgo intake, we assessed whether Ginkgo extract EGb 761^® affects coagulation or platelet function or increases the risk of bleeding. In a double-blind, placebo-controlled trial, prothrombin time, activated partial thromboplastin time, international normalized ratio and bleeding time were measured in patients with Alzheimer’s dementia at baseline, weeks 6 and 26. A total of 513 patients were randomized to 120 mg (n = 169) or 240 mg EGb 761^® (n = 170) or placebo (n = 174). No relevant changes were found for coagulation parameters and bleeding time. Numbers of bleeding-related adverse events were similar in all groups. Concomitant intake of acetylsalicylic acid was documented for 68 patients in the placebo group and 105 in the EGb 761^® groups. Within these groups, the means at baseline and week 26 differed by less than 1 unit for prothrombin time and bleeding time and less than 0.1 unit for international normalized ratio. Data on warfarin treatment in nine patients each taking placebo or EGb 761^® did not indicate enhancement of warfarin effects by EGb 761^®. No evidence was found that EGb 761^® affects hemostasis or increases the bleeding risk. No pharmacodynamic interactions with warfarin or acetylsalicylic acid were found.

A continued learning approach for model-informed precision dosing: updating models in clinical practice

Model-informed precision dosing (MIPD) is a quantitative dosing framework that combines prior knowledge on the drug-disease-patient system with patient data from therapeutic drug/biomarker monitoring (TDM) to support individualized dosing in ongoing treatment. Structural models and prior parameter distributions used in MIPD approaches typically build on prior clinical trials that involve only a limited number of patients selected according to some exclusion/inclusion criteria. Compared to the prior clinical trial population, the patient population in clinical practice can be expected to include also altered behavior and/or increased interindividual variability, the extent of which, however, is typically unknown. Here, we address the question of how to adapt and refine models on the level of the model parameters to better reflect this real-world diversity. We propose an approach for continued learning across patients during MIPD using a sequential hierarchical Bayesian framework. The approach builds on two stages to separate the update of the individual patient parameters from updating the population parameters. Consequently, it enables continued learning across hospitals or study centers, since only summary patient data (on the level of model parameters) need to be shared, but no individual TDM data. We illustrate this continued learning approach with neutrophil-guided dosing of paclitaxel. The present study constitutes an important step towards building confidence in MIPD and eventually establishing MIPD increasingly in everyday therapeutic use.

Evaluating prediction methods for glomerular filtration to optimise drug doses in obese and nonobese patients

AIMS

The most suitable method for predicting the glomerular filtration rate (GFR) in obesity is currently debated. Therefore, multiple GFR/creatinine clearance prediction methods were applied to (morbidly) obese and nonobese patients ranging from moderate renal impairment to glomerular hyperfiltration and their predictions were rated based on observed fosfomycin pharmacokinetics, as this model drug is exclusively eliminated via glomerular filtration.

METHODS

The GFR/creatinine clearance predictions via the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI), Modification of Diet in Renal Disease (MDRD; indexed and de-indexed by body surface area) and creatinine clearance via the Cockcroft-Gault formula (CLCR_CG ) using different body size descriptors were compared to the fosfomycin clearance (CL_FOF ) from 30 surgical patients (body mass index = 20.1-52.0 kg m^-2 ), receiving 8000 mg as intravenous infusion.

RESULTS

The concordance between CL_FOF and creatinine clearance predictions was highest for CLCR_CG employing either ideal body weight or adjusted body weight (if body mass >1.3 ideal body weight; CLCR_{CG_ABW-Schwartz} , concordance-correlation coefficient [95% confidence interval] = 0.474 [0.156; 0.703], CCC) and GFR predictions via the de-indexed MDRD equation (concordance-correlation coefficient = 0.452 [0.137; 0.685]). The proportion of predicted GFR values within ±30% of the observed CL_FOF (P₃₀  = 72.3-76.7%) was only marginally lower than the reported P₃₀ in the original CKD-EPI and MDRD publications (P₃₀  = 84.1-90.0%).

CONCLUSION

This analysis represents a successful proof-of-concept for evaluating GFR/creatinine clearance prediction methods: Across all body mass index classes CLCR_{CG_ABW-Schwartz} or the de-indexed MDRD were most suitable for predicting creatinine clearance/GFR also in (morbidly) obese, CKD stage <3B individuals in therapeutic use. Their application is proposed in optimising doses for vital therapies in obese patients requiring monitoring of renal function (e.g. methotrexate dosing).

Similar Piperacillin/Tazobactam Target Attainment in Obese versus Nonobese Patients despite Differences in Interstitial Tissue Fluid Pharmacokinetics

Precision dosing of piperacillin/tazobactam in obese patients is compromised by sparse information on target-site exposure. We aimed to evaluate the appropriateness of current and alternative piperacillin/tazobactam dosages in obese and nonobese patients. Based on a prospective, controlled clinical trial in 30 surgery patients (15 obese/15 nonobese; 0.5-h infusion of 4 g/0.5 g piperacillin/tazobactam), piperacillin pharmacokinetics were characterized in plasma and at target-site (interstitial fluid of subcutaneous adipose tissue) via population analysis. Thereafter, multiple 3-4-times daily piperacillin/tazobactam short-term/prolonged (recommended by EUCAST) and continuous infusions were evaluated by simulation. Adequacy of therapy was assessed by probability of pharmacokinetic/pharmacodynamic target-attainment (PTA ≥ 90%) based on time unbound piperacillin concentrations exceed the minimum inhibitory concentration (MIC) during 24 h (%fT>MIC). Lower piperacillin target-site maximum concentrations in obese versus nonobese patients were explained by the impact of lean (approximately two thirds) and fat body mass (approximately one third) on volume of distribution. Simulated steady-state concentrations were 1.43-times, 95%CI = (1.27; 1.61), higher in plasma versus target-site, supporting targets of %fT>2×MIC instead of %fT>4×MIC during continuous infusion to avoid target-site concentrations constantly below MIC. In all obesity and renally impairment/hyperfiltration stages, at MIC = 16 mg/L, adequate PTA required prolonged (thrice-daily 4 g/0.5 g over 3.0 h at %fT>MIC = 50) or continuous infusions (24 g/3 g over 24 h following loading dose at %fT>MIC = 98) of piperacillin/tazobactam.

Cost-effectiveness of oral anticancer drugs and associated individualised dosing approaches in patients with cancer: protocol for a systematic review

INTRODUCTION

Oral anticancer drugs (OADs) have rapidly expanded with more than 70 OADs targeting several molecular targets. Many of the OADs exert an exposure-response relationship but still, a 'one-size fits-all' dose is used, ignoring interindividual variability. Several of these OADs share similar mechanisms of actions and thus target the same cancer and has resulted in a substantial research focus on comparing the health benefit of each. However, significantly less is known about the cost-benefit associated with OADs. This paper will provide a protocol to systematically review studies that have evaluated the cost-effectiveness of OADs and their associated individualised dosing interventions.

METHODS AND ANALYSIS

Systematic review methodology will be applied to identify, select and extract data from published economic evaluation (costs and outcomes/benefits) studies of OADs and their associated individualised dosing interventions. Bibliographic databases (eg, Ovid EMBASE, Ovid MEDLINE) will be used to perform the systematic literature search (between 1 January 2000 and October 2020). Only full economic evaluations will be included, but no restrictions on study outcomes will be applied. The quality of included primary studies will be assessed using the Consolidated Health Economic Evaluation Reporting Standards checklist for reporting economic evaluations. Studies with low-quality evidence will be excluded. A narrative synthesis of the results from the included studies will be undertaken, with a subgroup analysis where appropriate.

ETHICS AND DISSEMINATION

This systematic review will not require ethics approval as there will not be any collection of primary data. Findings of this review will be disseminated through publications in peer-reviewed journals, presentations at workshops or conferences and sharing through a media release. Findings from this review will provide evidence to direct and inform policy-makers where cost-neutral strategies may be effective or where dose individualising strategies may be economically beneficial. Additionally, gaps will be identified in the current literature to inform future-related research.

PROSPERO REGISTRATION NUMBER

CRD42020218170.

ELECTRONIC SUPPLEMENTAL MATERIAL

The online version of this article contains supplemental material, which is available to authorised users.

Quantification of persister formation of Escherichia coli leveraging electronic cell counting and semi-mechanistic pharmacokinetic/pharmacodynamic modelling

BACKGROUND

Persister formation of Escherichia coli under fluoroquinolone exposure causes treatment failure and promotes emergence of resistant strains. Semi-mechanistic pharmacokinetic/pharmacodynamic modelling of data obtained from in vitro infection model experiments comprehensively characterizes exposure-effect relationships, providing mechanistic insights.

OBJECTIVES

To quantify persister formation of E. coli under levofloxacin exposure and to explain the observed growth-kill behaviour, leveraging electronic cell counting and pharmacokinetic/pharmacodynamic modelling.

METHODS

Three fluoroquinolone-resistant clinical E. coli isolates were exposed to levofloxacin in static and dynamic in vitro infection model experiments. Complementary to plate counting, bacterial concentrations over time were quantified by electronic cell counting and amalgamated in a semi-mechanistic pharmacokinetic/pharmacodynamic model (1281 bacterial and 394 levofloxacin observations).

RESULTS

Bacterial regrowth was observed under exposure to clinically relevant dosing regimens in the dynamic in vitro infection model. Electronic cell counting facilitated identification of three bacterial subpopulations: persisters, viable cells and dead cells. Two strain-specific manifestations of the levofloxacin effect were identified: a killing effect, characterized as a sigmoidal Emax model, and an additive increase in persister formation under levofloxacin exposure. Significantly different EC50 values quantitatively discerned levofloxacin potency for two isolates displaying the same MIC value: 8 mg/L [EC50 = 17.2 (95% CI = 12.6-23.8) mg/L and 8.46 (95% CI = 6.86-10.3) mg/L, respectively]. Persister formation was most pronounced for the isolate with the lowest MIC value (2 mg/L).

CONCLUSIONS

The developed pharmacokinetic/pharmacodynamic model adequately characterized growth-kill behaviour of three E. coli isolates and unveiled strain-specific levofloxacin potencies and persister formation. The mimicked dosing regimens did not eradicate the resistant isolates and enhanced persister formation to a strain-specific extent.

Perioperative administration of cefazolin and metronidazole in obese and non-obese patients: a pharmacokinetic study in plasma and interstitial fluid

OBJECTIVES

To assess plasma and tissue pharmacokinetics of cefazolin and metronidazole in obese patients undergoing bariatric surgery and non-obese patients undergoing intra-abdominal surgery.

PATIENTS AND METHODS

Fifteen obese and 15 non-obese patients received an IV short infusion of 2 g cefazolin and 0.5 g metronidazole for perioperative prophylaxis. Plasma and microdialysate from subcutaneous tissue were sampled until 8 h after dosing. Drug concentrations were determined by HPLC-UV. Pharmacokinetic parameters were calculated non-compartmentally.

RESULTS

In obese patients (BMI 39.5-69.3 kg/m2) compared with non-obese patients (BMI 18.7-29.8 kg/m2), mean Cmax of total cefazolin in plasma was lower (115 versus 174 mg/L) and Vss was higher (19.4 versus 14.2 L). The mean differences in t½ (2.7 versus 2.4 h), CL (5.14 versus 4.63 L/h) and AUC∞ (402 versus 450 mg·h/L) were not significant. The influence of obesity on the pharmacokinetics of metronidazole was similar (Cmax 8.99 versus 14.7 mg/L, Vss 73.9 versus 51.8 L, t½ 11.9 versus 9.1 h, CL 4.62 versus 4.13 L/h, AUC∞ 116 versus 127 mg·h/L). Regarding interstitial fluid (ISF), mean concentrations of cefazolin remained >4 mg/L until 6 h in both groups, and those of metronidazole up to 8 h in the non-obese group. In obese patients, the mean ISF concentrations of metronidazole were between 3 and 3.5 mg/L throughout the measuring interval.

CONCLUSIONS

During the time of surgery, cefazolin concentrations in plasma and ISF of subcutaneous tissue were lower in obese patients, but not clinically relevant. Regarding metronidazole, the respective differences were higher, and may influence dosing of metronidazole for perioperative prophylaxis in obese patients.

Rationale of a lower dexamethasone dose in prenatal congenital adrenal hyperplasia therapy based on pharmacokinetic modelling

CONTEXT

Prenatal dexamethasone therapy is used in female foetuses with congenital adrenal hyperplasia to suppress androgen excess and prevent virilisation of the external genitalia. The traditional dexamethasone dose of 20 µg/kg/day has been used since decades without examination in clinical trials and is thus still considered experimental.

OBJECTIVE

As the traditional dexamethasone dose potentially causes adverse effects in treated mothers and foetuses, we aimed to provide a rationale of a reduced dexamethasone dose in prenatal congenital adrenal hyperplasia therapy based on a pharmacokinetics-based modelling and simulation framework.

METHODS

Based on a published dexamethasone dataset, a nonlinear mixed-effects model was developed describing maternal dexamethasone pharmacokinetics. In stochastic simulations (n = 1000), a typical pregnant population (n = 124) was split into two dosing arms receiving either the traditional 20 µg/kg/day dexamethasone dose or reduced doses between 5 and 10 µg/kg/day. Target maternal dexamethasone concentrations, identified from the literature, served as a threshold to be exceeded by 90% of mothers at a steady state to ensure foetal hypothalamic-pituitary-adrenal axis suppression.

RESULTS

A two-compartment dexamethasone pharmacokinetic model was developed and subsequently evaluated to be fit for purpose. The simulations, including a sensitivity analysis regarding the assumed foetal:maternal dexamethasone concentration ratio, resulted in 7.5 µg/kg/day to be the minimum effective dose and thus our suggested dose.

CONCLUSIONS

We conclude that the traditional dexamethasone dose is three-fold higher than needed, possibly causing harm in treated foetuses and mothers. The clinical relevance and appropriateness of our recommended dose should be tested in a prospective clinical trial.

Microdialysis sampling to monitor target-site vancomycin concentrations in septic infants: a feasible way to close the knowledge gap

This work is dedicated to the memory of Hartmut Derendorf (1953-2020), a pioneer of modern pharmacokinetics and valued mentor of this project.

OBJECTIVES

Septic infants/neonates need effective antibiotic exposure, but dosing recommendations are challenging as the pharmacokinetics in this age are highly variable. For vancomycin, which is used as a standard treatment, comprehensive pharmacokinetic knowledge especially at the infection site is lacking. Hence, an exploratory clinical study was conducted to assess the feasibility and safety of microdialysis sampling for vancomycin monitoring at the target site.

METHODS

Nine infants/neonates with therapeutic indications for vancomycin treatment were administered 15 mg/kg as 1-hour infusions every 8-24 hours. Microdialysis catheters were implanted in the subcutaneous interstitial space fluid of the lateral thigh. Samples were collected every 30 minutes over 24 hours, followed by retrodialysis for catheter calibration. Prior in vitro investigations have evaluated impact factors on relative recovery and retrodialysis.

RESULTS

In vitro investigations showed the applicability of microdialysis for vancomycin monitoring. Microdialysis sampling was well tolerated in all infants/neonates (23-255 days) without major bleeding or other adverse events. Pharmacokinetic profiles were obtained and showed plausible vancomycin concentration-time courses.

CONCLUSIONS

Microdialysis as a minimally invasive technique for continuous longer-term sampling is feasible and safe in infants/neonates. Interstitial space fluid profiles were plausible and showed substantial interpatient variation. Hence, a larger microdialysis trial is warranted to further characterise the pharmacokinetics and variability of vancomycin at the target site and ultimately improve vancomycin dosing in these vulnerable patients.

The Biosynthetic Monophosphoryl Lipid A Enhances the Therapeutic Outcome of Antibiotic Therapy in Pneumococcal Pneumonia

Alternative treatment strategies against bacterial infections are required to decrease the use of antibiotics. This study tested the hypothesis that stimulation of the innate immune receptor Toll-like receptor 4 can be combined with antibiotics to improve the treatment of invasive pneumonia. The efficacy of the biosynthetic monophosphoryl lipid A (MPLA), a clinically approved Toll-like receptor 4 activator, was tested in a mouse model of Streptococcus pneumoniae respiratory infection. Interestingly, administration of amoxicillin or MPLA decreased 400- to 11 000-fold the bacterial load in the lungs and spleen but did not enhance survival compared to mock treatment. The single administration of a combination of MPLA and amoxicillin further reduced 10- to 18-fold the bacterial colonization and invasion and significantly improved protection against lethal disease. The combined administration of MPLA and amoxicillin in a context of infection was associated with transient increase of the serum concentrations of amoxicillin and pro-inflammatory cytokines and chemokines as well as the expression of immune genes in lung tissue. Remarkably, the systemic and lung immune activation extended beyond amoxicillin elimination, suggesting a two-step and cooperative anti-infective effect, i.e., rapid antibiotic-mediated alteration of bacteria and a long-lasting impact through mucosal and systemic immunity. Our proof-of-concept study demonstrated for the first time that boosting Toll-like receptor 4 signaling can synergize with antibiotics in order to increase the efficacy of therapy of bacterial pneumonia, thereby in fine reducing the dose or regimen of antibiotics.

Early Survival Prediction Framework in CD19-Specific CAR-T Cell Immunotherapy Using a Quantitative Systems Pharmacology Model

Chimeric antigen receptor (CAR)-T cell therapy has revolutionized treatment of relapsed/refractory non-Hodgkin lymphoma (NHL). However, since 36-60% of patients relapse, early response prediction is crucial. We present a novel population quantitative systems pharmacology model, integrating literature knowledge on physiology, immunology, and adoptive cell therapy together with 133 CAR-T cell phenotype, 1943 cytokine, and 48 metabolic tumor measurements. The model well described post-infusion concentrations of four CAR-T cell phenotypes and CD19+ metabolic tumor volume over 3 months after CAR-T cell infusion. Leveraging the model, we identified a low expansion subpopulation with significantly lower CAR-T cell expansion capacities amongst 19 NHL patients. Together with two patient-/therapy-related factors (autologous stem cell transplantation, CD4+/CD8+ T cells), the low expansion subpopulation explained 2/3 of the interindividual variability in the CAR-T cell expansion capacities. Moreover, the low expansion subpopulation had poor prognosis as only 1/4 of the low expansion subpopulation compared to 2/3 of the reference population were still alive after 24 months. We translated the expansion capacities into a clinical composite score (CCS) of 'Maximum naïve CAR-T cell concentrations/Baseline tumor burden' ratio and propose a CCSTN-value > 0.00136 (cells·µL-1·mL-1 as predictor for survival. Once validated in a larger cohort, the model will foster refining survival prediction and solutions to enhance NHL CAR-T cell therapy response.

Novel Pharmacokinetic/Pharmacodynamic Parameters Quantify the Exposure-Effect Relationship of Levofloxacin against Fluoroquinolone-Resistant Escherichia coli

Minimal inhibitory concentration-based pharmacokinetic/pharmacodynamic (PK/PD) indices are commonly applied to antibiotic dosing optimisation, but their informative value is limited, as they do not account for bacterial growth dynamics over time. We aimed to comprehensively characterise the exposure–effect relationship of levofloxacin against Escherichia coli and quantify strain-specific characteristics applying novel PK/PD parameters. In vitro infection model experiments were leveraged to explore the exposure–effect relationship of three clinical Escherichia coli isolates, harbouring different genomic fluoroquinolone resistance mechanisms, under constant levofloxacin concentrations or human concentration–time profiles (≤76 h). As an exposure metric, the ‘cumulative area under the levofloxacin–concentration time curve’ was determined. The antibiotic effect was assessed as the ‘cumulative area between the growth control and the bacterial-killing and -regrowth curve’. PK/PD modelling was applied to characterise the exposure–effect relationship and derive novel PK/PD parameters. A sigmoidal E_max model with an inhibition term best characterised the exposure–effect relationship and allowed for discrimination between two isolates sharing the same MIC value. Strain- and exposure-pattern-dependent differences were captured by the PK/PD parameters and elucidated the contribution of phenotypic adaptation to bacterial regrowth. The novel exposure and effect metrics and derived PK/PD parameters allowed for comprehensive characterisation of the isolates and could be applied to overcome the limitations of the MIC in clinical antibiotic dosing decisions, drug research and preclinical development.

Computational Treatment Simulations to Assess the Need for Personalized Tamoxifen Dosing in Breast Cancer Patients of Different Biogeographical Groups

Simple Summary

Tamoxifen is a drug often used to treat the most common type of breast cancer. Its metabolite endoxifen is formed by the liver enzyme CYP2D6, whose activity is variable and depends on a patient’s genetic profile. The frequency of CYP2D6 variants with different functional enzymatic activity varies largely between populations. To ensure sufficient efficacy of tamoxifen, a certain target concentration of endoxifen is needed, and 20% of White breast cancer patients have been shown not to reach this target concentration. However, little is known about the risk of not attaining the endoxifen target amongst other ethnic populations. This study investigated the risk for suboptimal endoxifen concentration in nine different biogeographical populations based on their distinct CYP2D6 genetic profile. The variability between the populations was high (up to three-fold), and East Asian breast cancer patients were identified as the population with the highest need for personalized tamoxifen dosing.

Abstract

Tamoxifen is used worldwide to treat estrogen receptor-positive breast cancer. It is extensively metabolized, and minimum steady-state concentrations of its metabolite endoxifen (C_{SS,min ENDX}) >5.97 ng/mL have been associated with favorable outcome. Endoxifen formation is mediated by the enzyme CYP2D6, and impaired CYP2D6 function has been associated with lower C_SS,min ENDX. In the Women’s Healthy Eating and Living (WHEL) study proposing the target concentration, 20% of patients showed subtarget C_{SS,min ENDX} at tamoxifen standard dosing. CYP2D6 allele frequencies vary largely between populations, and as 87% of the patients in the WHEL study were White, little is known about the risk for subtarget C_{SS,min ENDX} in other populations. Applying pharmacokinetic simulations, this study investigated the risk for subtarget C_{SS,min ENDX} at tamoxifen standard dosing and the need for dose individualization in nine different biogeographical groups with distinct CYP2D6 allele frequencies. The high variability in CYP2D6 allele frequencies amongst the biogeographical groups resulted in an up to three-fold difference in the percentages of patients with subtarget C_{SS,min ENDX}. Based on their CYP2D6 allele frequencies, East Asian breast cancer patients were identified as the population for which personalized, model-informed precision dosing would be most beneficial (28% of patients with subtarget C_{SS,min ENDX}).

The Use of Translational Modelling and Simulation to Develop Immunomodulatory Therapy as an Adjunct to Antibiotic Treatment in the Context of Pneumonia

The treatment of respiratory tract infections is threatened by the emergence of bacterial resistance. Immunomodulatory drugs, which enhance airway innate immune defenses, may improve therapeutic outcome. In this concept paper, we aim to highlight the utility of pharmacometrics and Bayesian inference in the development of immunomodulatory therapeutic agents as an adjunct to antibiotics in the context of pneumonia. For this, two case studies of translational modelling and simulation frameworks are introduced for these types of drugs up to clinical use. First, we evaluate the pharmacokinetic/pharmacodynamic relationship of an experimental combination of amoxicillin and a TLR4 agonist, monophosphoryl lipid A, by developing a pharmacometric model accounting for interaction and potential translation to humans. Capitalizing on this knowledge and associating clinical trial extrapolation and statistical modelling approaches, we then investigate the TLR5 agonist flagellin. The resulting workflow combines expert and prior knowledge on the compound with the in vitro and in vivo data generated during exploratory studies in order to construct high-dimensional models considering the pharmacokinetics and pharmacodynamics of the compound. This workflow can be used to refine preclinical experiments, estimate the best doses for human studies, and create an adaptive knowledge-based design for the next phases of clinical development.

Therapeutic drug monitoring of oral targeted antineoplastic drugs

PURPOSE

This review provides an overview of the current challenges in oral targeted antineoplastic drug (OAD) dosing and outlines the unexploited value of therapeutic drug monitoring (TDM). Factors influencing the pharmacokinetic exposure in OAD therapy are depicted together with an overview of different TDM approaches. Finally, current evidence for TDM for all approved OADs is reviewed.

METHODS

A comprehensive literature search (covering literature published until April 2020), including primary and secondary scientific literature on pharmacokinetics and dose individualisation strategies for OADs, together with US FDA Clinical Pharmacology and Biopharmaceutics Reviews and the Committee for Medicinal Products for Human Use European Public Assessment Reports was conducted.

RESULTS

OADs are highly potent drugs, which have substantially changed treatment options for cancer patients. Nevertheless, high pharmacokinetic variability and low treatment adherence are risk factors for treatment failure. TDM is a powerful tool to individualise drug dosing, ensure drug concentrations within the therapeutic window and increase treatment success rates. After reviewing the literature for 71 approved OADs, we show that exposure-response and/or exposure-toxicity relationships have been established for the majority. Moreover, TDM has been proven to be feasible for individualised dosing of abiraterone, everolimus, imatinib, pazopanib, sunitinib and tamoxifen in prospective studies. There is a lack of experience in how to best implement TDM as part of clinical routine in OAD cancer therapy.

CONCLUSION

Sub-therapeutic concentrations and severe adverse events are current challenges in OAD treatment, which can both be addressed by the application of TDM-guided dosing, ensuring concentrations within the therapeutic window.

A Model-Based Pharmacokinetic/Pharmacodynamic Analysis of the Combination of Amoxicillin and Monophosphoryl Lipid A Against S. pneumoniae in Mice

Combining amoxicillin with the immunostimulatory toll-like receptor 4 agonist monophosphoryl lipid A (MPLA) represents an innovative approach for enhancing antibacterial treatment success. Exploiting pharmacokinetic and pharmacodynamic data from an infection model of Streptococcus pneumoniae infected mice, we aimed to evaluate the preclinical exposure-response relationship of amoxicillin with MPLA coadministration and establish a link to survival. Antibiotic serum concentrations, bacterial numbers in lung and spleen and survival data of mice being untreated or treated with amoxicillin (four dose levels), MPLA, or their combination were analyzed by nonlinear mixed-effects modelling and time-to-event analysis using NONMEM® to characterize these treatment regimens. On top of a pharmacokinetic interaction, regarding the pharmacodynamic effects the combined treatment was superior to both monotherapies: The amoxicillin efficacy at highest dose was increased by a bacterial reduction of 1.74 log10 CFU/lung after 36 h and survival was increased 1.35-fold to 90.3% after 14 days both compared to amoxicillin alone. The developed pharmacometric pharmacokinetic/pharmacodynamic disease-treatment-survival models provided quantitative insights into a novel treatment option against pneumonia revealing a pharmacokinetic interaction and enhanced activity of amoxicillin and the immune system stimulator MPLA in combination. Further development of this drug combination flanked with pharmacometrics towards the clinical setting seems promising.

Constitutive Cell Proliferation Regulating Inhibitor of Protein Phosphatase 2A (CIP2A) Mediates Drug Resistance to Erlotinib in an EGFR Activating Mutated NSCLC Cell Line

Exploring mechanisms of drug resistance to targeted small molecule drugs is critical for an extended clinical benefit in the treatment of non-small cell lung cancer (NSCLC) patients carrying activating epidermal growth factor receptor (EGFR) mutations. Here, we identified constitutive cell proliferation regulating inhibitor of protein phosphatase 2A (CIP2A) in the HCC4006rErlo0.5 NSCLC cell line adapted to erlotinib as a model of acquired drug resistance. Constitutive CIP2A resulted in a constitutive activation of Akt signaling. The proteasome inhibitor bortezomib was able to reduce CIP2A levels, which resulted in an activation of protein phosphatase 2A and deactivation of Akt. Combination experiments with erlotinib and bortezomib revealed a lack of interaction between the two drugs. However, the effect size of bortezomib was higher in HCC4006rErlo0.5, compared to the erlotinib-sensitive HCC4006 cells, as indicated by an increase in Emax (0.911 (95%CI 0.867-0.954) vs. 0.585 (95%CI 0.568-0.622), respectively) and decrease in EC50 (52.4 µM (95%CI 46.1-58.8 µM) vs. 73.0 µM (95%CI 60.4-111 µM), respectively) in the concentration-effect model, an earlier onset of cell death induction, and a reduced colony surviving fraction (0.38 ± 0.18 vs. 0.95 ± 0.25, respectively, n = 3, p < 0.05). Therefore, modulation of CIP2A with bortezomib could be an interesting approach to overcome drug resistance to erlotinib treatment in NSCLC.

Which Analysis Approach Is Adequate to Leverage Clinical Microdialysis Data? A Quantitative Comparison to Investigate Exposure and Reponse Exemplified by Levofloxacin

Purpose

Systematic comparison of analysis methods of clinical microdialysis data for impact on target-site drug exposure and response.

Methods

39 individuals received a 500 mg levofloxacin short-term infusion followed by 24-h dense sampling in plasma and microdialysate collection in interstitial space fluid (ISF). ISF concentrations were leveraged using non-compartmental (NCA) and compartmental analysis (CA) via (ii) relative recovery correction at midpoint of the collection interval (midpoint-NCA, midpoint-CA) and (ii) dialysate-based integrals of time (integral-CA). Exposure and adequacy of community-acquired pneumonia (CAP) therapy via pharmacokinetic/pharmacodynamic target-attainment (PTA) analysis were compared between approaches.

Results

Individual AUC_ISF estimates strongly varied for midpoint-NCA and midpoint-CA (≥52.3%CV) versus integral-CA (≤32.9%CV) owing to separation of variability in PK parameters (midpoint-CA = 46.5%–143%CV_PK, integral-CA = 26.4%–72.6%CV_PK) from recovery-related variability only in integral-CA (41.0%–50.3%CV_recovery). This also led to increased variability of AUC_plasma for midpoint-CA (56.0%CV) versus midpoint-NCA and integral-CA (≤33.0%CV), and inaccuracy of predictive model performance of midpoint-CA in plasma (visual predictive check). PTA analysis translated into 33% of evaluated patient cases being at risk of incorrectly rejecting recommended dosing regimens at CAP-related epidemiological cut-off values.

Conclusions

Integral-CA proved most appropriate to characterise clinical pharmacokinetics- and microdialysis-related variability. Employing this knowledge will improve the understanding of drug target-site PK for therapeutic decision-making.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-021-02994-1.