Mechanistic Systems Modeling to Improve Understanding and Prediction of Cardiotoxicity Caused by Targeted Cancer Therapeutics

Tyrosine kinase inhibitors (TKIs) are highly potent cancer therapeutics that have been linked with serious cardiotoxicity, including left ventricular dysfunction, heart failure, and QT prolongation. TKI-induced cardiotoxicity is thought to result from interference with tyrosine kinase activity in cardiomyocytes, where these signaling pathways help to control critical processes such as survival signaling, energy homeostasis, and excitation–contraction coupling. However, mechanistic understanding is limited at present due to the complexities of tyrosine kinase signaling, and the wide range of targets inhibited by TKIs. Here, we review the use of TKIs in cancer and the cardiotoxicities that have been reported, discuss potential mechanisms underlying cardiotoxicity, and describe recent progress in achieving a more systematic understanding of cardiotoxicity via the use of mechanistic models. In particular, we argue that future advances are likely to be enabled by studies that combine large-scale experimental measurements with Quantitative Systems Pharmacology (QSP) models describing biological mechanisms and dynamics. As such approaches have proven extremely valuable for understanding and predicting other drug toxicities, it is likely that QSP modeling can be successfully applied to cardiotoxicity induced by TKIs. We conclude by discussing a potential strategy for integrating genome-wide expression measurements with models, illustrate initial advances in applying this approach to cardiotoxicity, and describe challenges that must be overcome to truly develop a mechanistic and systematic understanding of cardiotoxicity caused by TKIs

Systematic Evaluation of the Descriptive and Predictive Performance of Paediatric Morphine Population Models

Purpose: A framework for the evaluation of paediatric population models is proposed and applied to two different paediatric population pharmacokinetic models for morphine. One covariate model was based on a systematic covariate analysis, the other on fixed allometric scaling principles. Methods: The six evaluation criteria in the framework were 1) number of parameters and condition number, 2) numerical diagnostics, 3) prediction-based diagnostics, 4) η-shrinkage, 5) simulation-based diagnostics, 6) diagnostics of individual and population parameter estimates versus covariates, including measurements of bias and precision of the population values compared to the observed individual values. The framework entails both an internal and external model evaluation procedure. Results: The application of the framework to the two models resulted in the detection of overparameterization and misleading diagnostics based on individual predictions caused by high shrinkage. The diagnostic of individual and population parameter estimates versus covariates proved to be highly informative in assessing obtained covariate relationships. Based on the framework, the systematic covariate model proved to be superior over the fixed allometric model in terms of predictive performance. Conclusions: The proposed framework is suitable for the evaluation of paediatric (covariate) models and should be applied to corroborate the descriptive and predictive properties of these models

Effects of antibiotic resistance, drug target attainment, bacterial pathogenicity and virulence, and antibiotic access and affordability on outcomes in neonatal sepsis: an international microbiology and drug evaluation prospective substudy (BARNARDS)

Background Sepsis is a major contributor to neonatal mortality, particularly in low-income and middle-income countries (LMICs). WHO advocates ampicillin–gentamicin as first-line therapy for the management of neonatal sepsis. In the BARNARDS observational cohort study of neonatal sepsis and antimicrobial resistance in LMICs, common sepsis pathogens were characterised via whole genome sequencing (WGS) and antimicrobial resistance profiles. In this substudy of BARNARDS, we aimed to assess the use and efficacy of empirical antibiotic therapies commonly used in LMICs for neonatal sepsis. Methods In BARNARDS, consenting mother–neonates aged 0–60 days dyads were enrolled on delivery or neonatal presentation with suspected sepsis at 12 BARNARDS clinical sites in Bangladesh, Ethiopia, India, Pakistan, Nigeria, Rwanda, and South Africa. Stillborn babies were excluded from the study. Blood samples were collected from neonates presenting with clinical signs of sepsis, and WGS and minimum inhibitory concentrations for antibiotic treatment were determined for bacterial isolates from culture-confirmed sepsis. Neonatal outcome data were collected following enrolment until 60 days of life. Antibiotic usage and neonatal outcome data were assessed. Survival analyses were adjusted to take into account potential clinical confounding variables related to the birth and pathogen. Additionally, resistance profiles, pharmacokinetic–pharmacodynamic probability of target attainment, and frequency of resistance (ie, resistance defined by in-vitro growth of isolates when challenged by antibiotics) were assessed. Questionnaires on health structures and antibiotic costs evaluated accessibility and affordability. Findings Between Nov 12, 2015, and Feb 1, 2018, 36 285 neonates were enrolled into the main BARNARDS study, of whom 9874 had clinically diagnosed sepsis and 5749 had available antibiotic data. The four most commonly prescribed antibiotic combinations given to 4451 neonates (77·42%) of 5749 were ampicillin–gentamicin, ceftazidime–amikacin, piperacillin–tazobactam–amikacin, and amoxicillin clavulanate–amikacin. This dataset assessed 476 prescriptions for 442 neonates treated with one of these antibiotic combinations with WGS data (all BARNARDS countries were represented in this subset except India). Multiple pathogens were isolated, totalling 457 isolates. Reported mortality was lower for neonates treated with ceftazidime–amikacin than for neonates treated with ampicillin–gentamicin (hazard ratio [adjusted for clinical variables considered potential confounders to outcomes] 0·32, 95% CI 0·14–0·72; p=0·0060). Of 390 Gram-negative isolates, 379 (97·2%) were resistant to ampicillin and 274 (70·3%) were resistant to gentamicin. Susceptibility of Gram-negative isolates to at least one antibiotic in a treatment combination was noted in 111 (28·5%) to ampicillin–gentamicin; 286 (73·3%) to amoxicillin clavulanate–amikacin; 301 (77·2%) to ceftazidime–amikacin; and 312 (80·0%) to piperacillin–tazobactam–amikacin. A probability of target attainment of 80% or more was noted in 26 neonates (33·7% [SD 0·59]) of 78 with ampicillin–gentamicin; 15 (68·0% [3·84]) of 27 with amoxicillin clavulanate–amikacin; 93 (92·7% [0·24]) of 109 with ceftazidime–amikacin; and 70 (85·3% [0·47]) of 76 with piperacillin–tazobactam–amikacin. However, antibiotic and country effects could not be distinguished. Frequency of resistance was recorded most frequently with fosfomycin (in 78 isolates [68·4%] of 114), followed by colistin (55 isolates [57·3%] of 96), and gentamicin (62 isolates [53·0%] of 117). Sites in six of the seven countries (excluding South Africa) stated that the cost of antibiotics would influence treatment of neonatal sepsis

Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies

Early prediction, quantification and translation of cardiovascular hemodynamic drug effects is essential in pre-clinical drug development. In this study, a novel hemodynamic cardiovascular systems (CVS) model was developed to support these goals. The model consisted of distinct system- and drug-specific parameter, and uses data for heart rate (HR), cardiac output (CO), and mean atrial pressure (MAP) to infer drug mode-of-action (MoA). To support further application of this model in drug development, we conducted a systematic analysis of the estimation performance of the CVS model to infer drug- and system-specific parameters. Specifically, we focused on the impact on model estimation performance when considering differences in available readouts and the impact of study design choices. To this end, a practical identifiability analysis was performed, evaluating model estimation performance for different combinations of hemodynamic endpoints, drug effect sizes, and study design characteristics. The practical identifiability analysis showed that MoA of drug effect could be identified for different drug effect magnitudes and both system- and drug-specific parameters can be estimated precisely with minimal bias. Study designs which exclude measurement of CO or use a reduced measurement duration still allow the identification and quantification of MoA with acceptable performance. In conclusion, the CVS model can be used to support the design and inference of MoA in pre-clinical CVS experiments, with a future potential for applying the uniquely identifiable systems parameters to support inter-species scaling

Anti-cancer treatment schedule optimization based on tumor dynamics modelling incorporating evolving resistance

Quantitative characterization of evolving tumor resistance under targeted treatment could help identify novel treatment schedules, which may improve the outcome of anti-cancer treatment. In this study, a mathematical model which considers various clonal populations and evolving treatment resistance was developed. With parameter values fitted to the data or informed by literature data, the model could capture previously reported tumor burden dynamics and mutant KRAS levels in circulating tumor DNA (ctDNA) of patients with metastatic colorectal cancer treated with panitumumab. Treatment schedules, including a continuous schedule, intermittent schedules incorporating treatment holidays, and adaptive schedules guided by ctDNA measurements were evaluated using simulations. Compared with the continuous regimen, the simulated intermittent regimen which consisted of 8-week treatment and 4-week suspension prolonged median progression-free survival (PFS) of the simulated population from 36 to 44 weeks. The median time period in which the tumor size stayed below the baseline level (T-TS<TS0) was prolonged from 52 to 60 weeks. Extending the treatment holiday resulted in inferior outcomes. The simulated adaptive regimens showed to further prolong median PFS to 56-64 weeks and T-TS<TS0 to 114-132 weeks under different treatment designs. A prospective clinical study is required to validate the results and to confirm the added value of the suggested schedules

Population Pharmacokinetics of Vancomycin in Premature Malaysian Neonates: Identification of Predictors for Dosing Determination▿

The present study determined the pharmacokinetic profile of vancomycin in premature Malaysian infants. A one-compartment infusion model with first-order elimination was fitted to serum vancomycin concentration data (n = 835 points) obtained retrospectively from the drug monitoring records of 116 premature newborn infants. Vancomycin concentrations were estimated by a fluorescence polarization immunoassay. Population and individual estimates of clearance and distribution volume and the factors which affected the variability observed for the values of these parameters were obtained using a population pharmacokinetic modeling approach. The predictive performance of the population model was evaluated by visual inspections of diagnostic plots and nonparametric bootstrapping with replacement. Dosing guidelines targeting a value of ≥400 for the area under the concentration-time curve over 24 h in the steady state divided by the MIC (AUC24/MIC ratio) were explored using Monte Carlo simulation. Body size (weight), postmenstrual age, and small-for-gestational-age status are important factors explaining the between-subject variability of vancomycin pharmacokinetic parameter values for premature neonates. The typical population parameter estimates of clearance and distribution volume for a 1-kg premature appropriate-for-gestational-age neonate with a postmenstrual age of 30 weeks were 0.0426 liters/h and 0.523 liters, respectively. There was a 20% reduction in clearance for small-for-gestational-age infants compared to the level for the appropriate-for-gestational-age control. Dosage regimens based on a priori target response values were formulated. In conclusion, the pharmacokinetic parameter values for vancomycin in premature Malaysian neonates were estimated. Improved dosage regimens based on a priori target response values were formulated by incorporating body size, postmenstrual age, and small-for-gestational-age status, using Monte Carlo simulations with the model-estimated pharmacokinetic parameter values

Severe skin toxicity in pediatric oncology patients treated with voriconazole and concomitant methotrexate

We report the occurrence of skin toxicities in pediatric oncology patients on concomitant treatment with voriconazole and methotrexate (MTX). Of 23 patients who received this combination, 11 patients suffered from cheilitis and/or photosensitivity. In contrast, only in 1 of 9 patients who received voriconazole without MTX was photosensitivity observed. A mechanism of action was not able to be identified. We describe two cases with severe skin toxicities. Caution is warranted when using voriconazole and concomitant MT

Quantitative modeling of tumor dynamics and development of drug resistance in non-small cell lung cancer patients treated with erlotinib

Insight into the development of treatment resistance can support the optimization of anticancer treatments. This study aims to characterize the tumor dynamics and development of drug resistance in patients with non-small cell lung cancer treated with erlotinib, and investigate the relationship between baseline circulating tumor DNA (ctDNA) data and tumor dynamics. Data obtained for the analysis included (1) intensively sampled erlotinib concentrations from 29 patients from two previous pharmacokinetic (PK) studies, and (2) tumor sizes, ctDNA measurements, and sparsely sampled erlotinib concentrations from 18 patients from the START-TKI study. A two-compartment population PK model was first developed which well-described the PK data. The PK model was subsequently applied to investigate the exposure-tumor dynamics relationship. To characterize the tumor dynamics, models accounting for intra-tumor heterogeneity and acquired resistance with or without primary resistance were investigated. Eventually, the model assumed acquired resistance only resulted in an adequate fit. Additionally, models with or without exposure-dependent treatment effect were explored, and no significant exposure-response relationship for erlotinib was identified within the observed exposure range. Subsequently, the correlation of baseline ctDNA data on EGFR and TP53 variants with tumor dynamics’ parameters was explored. The analysis indicated that higher baseline plasma EGFR mutation levels correlated with increased tumor growth rates, and the inclusion of ctDNA measurements improved model fit. This result suggests that quantitative ctDNA measurements at baseline have the potential to be a predictor of anticancer treatment response. The developed model can potentially be applied to design optimal treatment regimens that better overcome resistance.</p

Pooled Population Pharmacokinetic Analysis for Exploring Ciprofloxacin Pharmacokinetic Variability in Intensive Care Patients

Background and Objective: Previous pharmacokinetic (PK) studies of ciprofloxacin in intensive care (ICU) patients have shown large differences in estimated PK parameters, suggesting that further investigation is needed for this population. Hence, we performed a pooled population PK analysis of ciprofloxacin after intravenous administration using individual patient data from three studies. Additionally, we studied the PK differences between these studies through a post-hoc analysis. Methods: Individual patient data from three studies (study 1, 2, and 3) were pooled. The pooled data set consisted of 1094 ciprofloxacin concentration–time data points from 140 ICU patients. Nonlinear mixed-effects modeling was used to develop a population PK model. Covariates were selected following a stepwise covariate modeling procedure. To analyze PK differences between the three original studies, random samples were drawn from the posterior distribution of individual PK parameters. These samples were used for a simulation study comparing PK exposure and the percentage of target attainment between patients of these studies. Results: A two-compartment model with first-order elimination best described the data. Inter-individual variability was added to the clearance, central volume, and peripheral volume. Inter-occasion variability was added to clearance only. Body weight was added to all parameters allometrically. Estimated glomerular filtration rate on ciprofloxacin clearance was identified as the only covariate relationship resulting in a drop in inter-individual variability of clearance from 58.7 to 47.2%. In the post-hoc analysis, clearance showed the highest deviation between the three studies with a coefficient of variation of 14.3% for posterior mean and 24.1% for posterior inter-individual variability. The simulation study showed that following the same dose regimen of 400 mg three times daily, the area under the concentration–time curve of study 3 was the highest with a mean area under the concentration–time curve at 24 h of 58 mg·h/L compared with that of 47.7 mg·h/L for study 1 and 47.6 mg·h/L for study 2. Similar differences were also observed in the percentage of target attainment, defined as the ratio of area under the concentration–time curve at 24 h and the minimum inhibitory concentration. At the epidemiological cut-off minimum inhibitory concentration of Pseudomonas aeruginosa of 0.5 mg/L, percentage of target attainment was only 21%, 18%, and 38% for study 1, 2, and 3, respectively. Conclusions: We developed a population PK model of ciprofloxacin in ICU patients using pooled data of individual patients from three studies. A simple ciprofloxacin dose recommendation for the entire ICU population remains challenging owing to the PK differences within ICU patients, hence dose individualization may be needed for the optimization of ciprofloxacin treatment

Why we should sample sparsely and aim for a higher target: Lessons from model-based therapeutic drug monitoring of vancomycin in intensive care patients

Aims: To explore the optimal data sampling scheme and the pharmacokinetic (PK) target exposure on which dose computation is based in the model-based therapeutic drug monitoring (TDM) practice of vancomycin in intensive care (ICU) patients. Methods: We simulated concentration data for 1 day following four sampling schemes, Cmin, Cmax + Cmin, Cmax + Cmid-interval + Cmin, and rich sampling where a sample was drawn every hour within a dose interval. The datasets were used for Bayesian estimation to obtain PK parameters, which were used to compute the doses for the next day based on five PK target exposures: AUC24 = 400, 500, and 600 mg·h/L and Cmin = 15 and 20 mg/L. We then simulated data for the next day, adopting the computed doses, and repeated the above procedure for 7 days. Thereafter, we calculated the percentage error and the normalized root mean square error (NRMSE) of estimated against “true” PK parameters, and the percentage of optimal treatment (POT), defined as the percentage of patients who met 400 ≤ AUC24 ≤ 600 mg·h/L and Cmin ≤ 20 mg/L. Results: PK parameters were unbiasedly estimated in all investigated scenarios and the 6-day average NRMSE were 32.5%/38.5% (CL/V, where CL is clearance and V is volume of distribution) in the trough sampling scheme and 27.3%/26.5% (CL/V) in the rich sampling scheme. Regarding POT, the sampling scheme had marginal influence, while target exposure showed clear impacts that the maximum POT of 71.5% was reached when doses were computed based on AUC24 = 500 mg·h/L. Conclusions: For model-based TDM of vancomycin in ICU patients, sampling more frequently than taking only trough samples adds no value and dosing based on AUC24 = 500 mg·h/L lead to the best POT