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.

Novel mathematical approach to accurately quantify 3D endothelial cell morphology and vessel geometry based on fluorescently marked endothelial cell contours: Application to the dorsal aorta of wild-type and Endoglin-deficient zebrafish embryos.

Endothelial cells, which line the lumen of blood vessels, locally sense and respond to blood flow. In response to altered blood flow dynamics during early embryonic development, these cells undergo shape changes that directly affect vessel geometry: In the dorsal aorta of zebrafish embryos, elongation of endothelial cells in the direction of flow between 48 and 72 hours post fertilization (hpf) reduces the vessel's diameter. This remodeling process requires Endoglin; excessive endothelial cell growth in the protein's absence results in vessel diameter increases. To understand how these changes in vessel geometry emerge from morphological changes of individual endothelial cells, we developed a novel mathematical approach that allows 3D reconstruction and quantification of both dorsal aorta geometry and endothelial cell surface morphology. Based on fluorescently marked endothelial cell contours, we inferred cross-sections of the dorsal aorta that accounted for dorsal flattening of the vessel. By projection of endothelial cell contours onto the estimated cross-sections and subsequent triangulation, we finally reconstructed 3D surfaces of the individual cells. By simultaneously reconstructing vessel cross-sections and cell surfaces, we found in an exploratory analysis that morphology varied between endothelial cells located in different sectors of the dorsal aorta in both wild-type and Endoglin-deficient zebrafish embryos: In wild-types, ventral endothelial cells were smaller and more elongated in flow direction than dorsal endothelial cells at both 48 hpf and 72 hpf. Although dorsal and ventral endothelial cells in Endoglin-deficient embryos had similar sizes at 48 hpf, dorsal endothelial cells were much larger at 72 hpf. In Endoglin-deficient embryos, elongation in flow direction increased between 48 hpf and 72 hpf in ventral endothelial cells but hardly changed in dorsal endothelial cells. Hereby, we provide evidence that dorsal endothelial cells contribute most to the disparate changes in dorsal aorta diameter in wild-type and Endoglin-deficient embryos between 48 hpf and 72 hpf.

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.

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 AUC 24 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.

A quantitative modeling framework to understand the physiology of the hypothalamic-pituitary-adrenal axis and interaction with cortisol replacement therapy.

Congenital adrenal hyperplasia (CAH) is characterized by impaired adrenal cortisol production. Hydrocortisone (synthetic cortisol) is the drug-of-choice for cortisol replacement therapy, aiming to mimic physiological cortisol circadian rhythm. The hypothalamic-pituitary-adrenal (HPA) axis controls cortisol production through the pituitary adrenocorticotropic hormone (ACTH) and feedback mechanisms. The aim of this study was to quantify key mechanisms involved in the HPA axis activity regulation and their interaction with hydrocortisone therapy. Data from 30 healthy volunteers was leveraged: Endogenous ACTH and cortisol concentrations without any intervention as well as cortisol concentrations measured after dexamethasone suppression and single dose administration of (i) 0.5-10 mg hydrocortisone as granules, (ii) 20 mg hydrocortisone as granules and intravenous bolus. A stepwise model development workflow was used: A newly developed model for endogenous ACTH and cortisol was merged with a refined hydrocortisone pharmacokinetic model. The joint model was used to simulate ACTH and cortisol trajectories in CAH patients with varying degrees of enzyme deficiency, with or without hydrocortisone administration, and healthy individuals. Time-dependent ACTH-driven endogenous cortisol production and cortisol-mediated feedback inhibition of ACTH secretion processes were quantified and implemented in the model. Comparison of simulated ACTH and cortisol trajectories between CAH patients and healthy individuals showed the importance of administering hydrocortisone before morning ACTH secretion peak time to suppress ACTH overproduction observed in untreated CAH patients. The developed framework allowed to gain insights on the physiological mechanisms of the HPA axis regulation, its perturbations in CAH and interaction with hydrocortisone administration, paving the way towards cortisol replacement therapy optimization.

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.

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.

Combination of Pharmacokinetic and Pathogen Susceptibility Information To Optimize Meropenem Treatment of Gram-Negative Infections in Critically Ill Patients.

Meropenem is one of the most frequently used antibiotics to treat life-threatening infections in critically ill patients. This study aimed to develop a meropenem dosing algorithm for the treatment of Gram-negative infections based on intensive care unit (ICU)-specific resistance data. Antimicrobial susceptibility testing of Gram-negative bacteria obtained from critically ill patients was carried out from 2016 to 2020 at a tertiary care hospital. Based on the observed MIC distribution, stochastic simulations (n = 1,000) of an evaluated pharmacokinetic meropenem model, and a defined pharmacokinetic/pharmacodynamic target (100%T>4×MIC while minimum concentrations were <44.5 mg/L), dosing recommendations for patients with varying renal function were derived. Pathogen-specific MIC distributions were used to calculate the cumulative fraction of response (CFR), and the overall MIC distribution was used to calculate the local pathogen-independent mean fraction of response (LPIFR) for the investigated dosing regimens. A CFR/LPIFR of >90% was considered adequate. The observed MIC distribution significantly differed from the EUCAST database. Based on the 6,520 MIC values included, a three-level dosing algorithm was developed. If the pathogen causing the infection is unknown (level 1), known (level 2), known to be neither Pseudomonas aeruginosa nor Acinetobacter baumannii, or classified as susceptible (level 3), a continuous infusion of 1.5 g daily reached sufficient target attainment independent of renal function. In all other cases, dosing needs to be adjusted based on renal function. ICU-specific susceptibility data should be assessed regularly and integrated into dosing decisions. The presented workflow may serve as a blueprint for other antimicrobial settings.

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 (CLCRCG ) using different body size descriptors were compared to the fosfomycin clearance (CLFOF ) from 30 surgical patients (body mass index = 20.1-52.0 kg m-2 ), receiving 8000 mg as intravenous infusion. RESULTS: The concordance between CLFOF and creatinine clearance predictions was highest for CLCRCG employing either ideal body weight or adjusted body weight (if body mass >1.3 ideal body weight; CLCRCG_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 CLFOF (P30 = 72.3-76.7%) was only marginally lower than the reported P30 in the original CKD-EPI and MDRD publications (P30 = 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 CLCRCG_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).

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.

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 AUCNO/AUCVRC 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.

Target Site Pharmacokinetics of Meropenem: Measurement in Human Explanted Lung Tissue by Bronchoalveolar Lavage, Microdialysis, and Homogenized Lung Tissue.

Pneumonia is one of the most common infections in intensive care patients, and it is often treated with beta-lactam antibiotics. Even if therapeutic drug monitoring in blood is available, it is unclear whether sufficient concentrations are reached at the target site: the lung. The present study was initiated to fill this knowledge gap. Various compartments from 10 patients' explanted lungs were subjected to laboratory analysis. Meropenem was quantified in serum, bronchoalveolar lavage (BAL) fluid, microdialysate, and homogenized lung tissue with isotope dilution liquid chromatography tandem mass spectrometry (ID-LC-MS/MS). BAL fluid represents diluted epithelial lining fluid (ELF), and microdialysate represents interstitial fluid (IF). Differences between target site and blood concentrations were investigated. The median meropenem concentration in blood, ELF, IF, and tissue were 26.8, 18.0, 12.1, and 9.1 mg/liter, respectively. A total of 37.5% of the target site ELF and IF meropenem concentrations were below the clinical EUCAST breakpoint of 8 mg/liter. The median ELF/serum quotient was 61.8% (interquartile range [IQR], 24.8% to 87.6%), the median IF/serum quotient was 35.4% (IQR, 23.8% to 54.3%), and the median tissue/serum quotient was 34.2% (IQR, 28.3% to 38.2%). We observed a substantial interindividual variability between the blood and the compartments (ELF and IF), whereas the intraindividual variability was relatively low. Target site measurement in different lung compartments was feasible and successfully applied in a clinical setting. A relevant amount of 37.5% of the target site concentrations were below the clinical EUCAST breakpoint, indicating subtherapeutic dosing in high-risk patients receiving perioperative antibiotic prophylaxis in lung transplantation. (This study has been registered at ClinicalTrials.gov under identifier NCT03970265.).

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.

Quantitative relationship between infliximab exposure and inhibition of C-reactive protein synthesis to support inflammatory bowel disease management.

AIM: Quantitative and kinetic insights into the drug exposure-disease response relationship might enhance our knowledge on loss of response and support more effective monitoring of inflammatory activity by biomarkers in patients with inflammatory bowel disease (IBD) treated with infliximab (IFX). This study aimed to derive recommendations for dose adjustment and treatment optimisation based on mechanistic characterisation of the relationship between IFX serum concentration and C-reactive protein (CRP) concentration. METHODS: Data from an investigator-initiated trial included 121 patients with IBD during IFX maintenance treatment. Serum concentrations of IFX, antidrug antibodies (ADA), CRP, and disease-related covariates were determined at the mid-term and end of a dosing interval. Data were analysed using a pharmacometric nonlinear mixed-effects modelling approach. An IFX exposure-CRP model was generated and applied to evaluate dosing regimens to achieve CRP remission. RESULTS: The generated quantitative model showed that IFX has the potential to inhibit up to 72% (9% relative standard error [RSE]) of CRP synthesis in a patient. IFX concentration leading to 90% of the maximum CRP synthesis inhibition was 18.4 µg/mL (43% RSE). Presence of ADA was the most influential factor on IFX exposure. With standard dosing strategy, ≥55% of ADA+ patients experienced CRP nonremission. Shortening the dosing interval and co-therapy with immunomodulators were found to be the most beneficial strategies to maintain CRP remission. CONCLUSIONS: With the generated model we could for the first time establish a robust relationship between IFX exposure and CRP synthesis inhibition, which could be utilised for treatment optimisation in IBD patients.

Role of TDM-based dose adjustments for taxane anticancer drugs.

The classical taxanes (paclitaxel, docetaxel), the newer taxane cabazitaxel and the nanoparticle-bound nab-paclitaxel are among the most widely used anticancer drugs. Still, the optimal use and the value of pharmacological personalization of the taxanes is still controversial. We give an overview on the pharmacological properties of the taxanes, including metabolism, pharmacokinetics-pharmacodynamic relations and aspects in the clinical use of taxanes. The latter includes the ongoing debate on the most effective and safe regimen, the recommended initial dose, and pharmacological dosing individualization. The taxanes are among the most widely used anticancer drugs in patients with solid malignancies. Despite their longtime use in clinical routine, the optimal dosing strategy (weekly versus 3-weekly) or optimal average dose (cabazitaxel, nab-paclitaxel) has not been fully resolved, as it may differ according to tumour entity and line of treatment. The value of pharmacological individualization of the taxanes (TDM, TCI) has been partly explored for 3-weekly paclitaxel and docetaxel, but remains mostly unexplored for cabazitaxel and nab-paclitaxel at present.

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 AUCISF 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%CVPK, integral-CA = 26.4%-72.6%CVPK) from recovery-related variability only in integral-CA (41.0%-50.3%CVrecovery). This also led to increased variability of AUCplasma 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.

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.

CYP2D6 phenotype explains reported yohimbine concentrations in four severe acute intoxications.

The indole alkaloid yohimbine is an alpha-2 receptor antagonist used for its sympathomimetic effects. Several cases of yohimbine intoxication have been reported and the most recent one involved four individuals taking a yohimbine-containing drug powder. All individuals developed severe intoxication symptoms and were admitted to the hospital. Even though all individuals were assumed to have taken the same dose of the drug powder, toxicology analyses revealed yohimbine blood concentrations of 249-5631 ng/mL, amounting to a 22-fold difference. The reason for this high variability remained to be elucidated. We used recently reported knowledge on the metabolism of yohimbine together with state-of-the art nonlinear mixed-effects modelling and simulation and show that a patient's cytochrome P450 2D6 (CYP2D6) phenotype can explain the large differences observed in the measured concentration after intake of the same yohimbine dose. Our findings can be used both for the identification of safe doses in therapeutic use of yohimbine and for an explanation of individual cases of overdosing.

Time-to-Event Analysis of Paclitaxel-Associated Peripheral Neuropathy in Advanced Non-Small-Cell Lung Cancer Highlighting Key Influential Treatment/Patient Factors.

Paclitaxel-associated peripheral neuropathy (PN), a major dose-limiting toxicity, significantly impacts patients' quality of life/treatment outcome. Evaluation of risk factors often ignores time of PN onset, precluding the impact of time-dependent factors, e.g., drug exposure, needed to comprehensively characterize PN. We employed parametric time-to-event (TTE) analysis to describe the time course of risk of first occurrence of clinically relevant PN grades ≥2 (PN2+, n = 105, common terminology criteria v4.0) and associated patient/treatment characteristics, leveraging data from 365 patients (1454 cycles) receiving paclitaxel every 3 weeks (plus carboplatin AUC = 6 or cisplatin 80 mg/m2) for ≤6 cycles. Paclitaxel was intravenously administered (3 hours) as standard 200-mg/m2 doses (n = 182) or as pharmacokinetic-guided dosing (n = 183). A cycle-varying hazard TTE model linking surge in hazard of PN2+ to paclitaxel administration [PN2+ proportions (i.e., cases per 1000 patients), 1st day, cycle 1: 4.87 of 1000; cycle 6: 7.36 of 1000] and linear decline across cycle (last day, cycle 1: 1.64 of 1000; cycle 6: 2.48 of 1000) adequately characterized the time-varying hazard of PN2+. From joint covariate evaluation, PN2+ proportions (1st day, cycle 1) increased by 1.00 per 1000 with 5-µmol·h/l higher paclitaxel exposure per cycle (AUC between the start and end of a cycle, most relevant covariate), 0.429 per 1000 with 5-year higher age, 1.31 per 1000 (smokers vs. nonsmokers), and decreased by 0.670 per 1000 (females vs. males). Compared to 200 mg/m2 dosing every 3 weeks, model-predicted cumulative risk of PN2+ was significantly higher (42%) with 80 mg/m2 weekly dosing but reduced by 11% with 175 mg/m2 dosing every 3 weeks. The established TTE modeling framework enables quantification and comparison of patient's cumulative risks of PN2+ for different clinically relevant paclitaxel dosing schedules, sparing patients PN2+ to improve paclitaxel therapy. SIGNIFICANCE STATEMENT: Characterization of risk factors of paclitaxel-associated peripheral neuropathy (PN) typically involves time-independent comparison of PN odds in patient subpopulations, concealing the impact of time-dependent factors, e.g., changing paclitaxel exposure, required to comprehensively characterize PN. We developed a parametric time-to-event model describing the time course in risk of clinically relevant paclitaxel-associated PN, identifying the highest risk in older male smokers with higher paclitaxel area under the plasma concentration-time curve between the start and end of a cycle. The developed framework enabled quantification of patient's risk of PN for clinically relevant paclitaxel dosing schedules, facilitating future dosing decisions.