Search

The tumor size ratio (TSR), time-to-tumor growth (TTG), and tumor growth rate (kG) are frequently suggested as model-based predictors of overall survival (OS) for different types of tumors. When the tumor metrics are applied in forecasting of the outcome for individual patients at an early stage, the tumor data might be sparse resulting in imprecise prediction. This simulation study aimed to investigate how the tumor follow-up data and estimation approaches influence the accuracy in the tumor size metrics and the predicted hazard of death for individual patients. Longitudinal tumor size and OS data were simulated using tumor growth inhibition and Weibull distribution models, respectively. Based on the model and increasing measurement durations, the accuracy (defined as 80–125% of the simulated “true” value) in individual metrics and hazard was computed. TSR week 6 (TSRw6) accuracy was adequate for 91% of the patients when tumor size was measured up to 12 weeks. For TTG and kG metrics, the highest accuracy observed was lower (43 and 77%, respectively) and occurred later (42 and 60 weeks, respectively). The simultaneous (joint) and sequential estimation approaches resulted in similar accuracies, however, in general, the sequential approach where individual tumor size parameters are fixed, demonstrated inferior estimation properties. The TSRw6 and the model-predicted tumor time course (absolute or relative change) had better forecasting properties than TTG or kG. The population pharmacokinetic (PK) parameters and data approach performed similarly or better than the simultaneous approach and had a better accuracy in estimating individuals' hazard of death than the individual PK parameters method. Title in thesis list of papers of Sreenath Madathil Krishnan: Bayesian forecasting of tumor size metrics and survival

Both initiation and suppression of inflammation are hallmarks of the immune response. If not balanced, the inflammation may cause extensive tissue damage, which is associated with common diseases, e.g., asthma and atherosclerosis. Anti-inflammatory drugs come with side effects that may be aggravated by high and fluctuating drug concentrations. To remedy this, an anti-inflammatory drug should have an appropriate pharmacokinetic half-life or better still, a sustained anti-inflammatory drug response. However, we still lack a quantitative mechanistic understanding of such sustained effects. Here, we study the anti-inflammatory response to a common glucocorticoid drug, dexamethasone. We find a sustained response 22 hours after drug removal. With hypothesis testing using mathematical modeling, we unravel the underlying mechanism-a slow release of dexamethasone from the receptor-drug complex. The developed model is in agreement with time-resolved training and testing data and is used to simulate hypothetical treatment schemes. This work opens up for a more knowledge-driven drug development to find sustained anti-inflammatory responses and fewer side effects. Funding Agencies|Swedish Research CouncilSwedish Research Council [Dnr 2019-03767]; Heart and Lung FoundationSwedish Heart-Lung Foundation; CENIIT (Center for Industrial Information Technology); Swedish Research CouncilSwedish Research Council; Horizon 2020; KAW and Sci Life Lab Covid-19 platform

Intradermal delivery of AZD8601, an mRNA designed to produce vascular endothelial growth factor A (VEGF-A), has previously been shown to accelerate cutaneous wound healing in a murine diabetic model. Here, we develop population pharmacokinetic and pharmacodynamic models aiming to quantify the effect of AZD8601 injections on the dynamics of wound healing. A dataset of 584 open wound area measurements from 131 mice was integrated from 3 independent studies encompassing different doses, dosing timepoints, and number of doses. Evaluation of several candidate models showed that wound healing acceleration is not likely driven directly by time-dependent VEGF-A concentration. Instead, we found that administration of AZD8601 induced a sustained acceleration of wound healing depending on the accumulated dose, with a dose producing 50% of the maximal effect of 92 mu g. Simulations with this model showed that a single dose of 200 mu g AZD8601 can reduce the time to reach 50% wound healing by up to 5 days. Funding Agencies|Integrated Cardiometabolic Center, Karolinska Institute, Stockholm, Sweden; Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden

Bile acids released postprandially can modify the rate and extent of lipophilic compounds’ absorption. This study aimed to predict the enterohepatic circulation (EHC) of total bile acids (TBAs) in response to caloric intake from their spillover in plasma. A model for TBA EHC was combined with a previously developed gastric emptying (GE) model. Longitudinal gallbladder volumes and TBA plasma concentration data from 30 subjects studied after ingestion of four different test drinks were supplemented with literature data. Postprandial gallbladder refilling periods were implemented to improve model predictions. The TBA hepatic extraction was reduced with the high‐fat drink. Basal and nutrient‐induced gallbladder emptying rates were altered by type 2 diabetes (T2D). The model was predictive of the central trend and the variability of gallbladder volume and TBA plasma concentration for all test drinks. Integration of this model within physiological pharmacokinetic modeling frameworks could improve the predictions for lipophilic compounds’ absorption considerably.

In antidiabetic drug development, phase I studies usually involve short-term glucose provocations. Multiple designs are available for these provocations (e.g., meal tolerance tests (MTTs) and graded glucose infusions (GGIs)). With a highly nonlinear, complex system as the glucose homeostasis, the various provocations will contribute with different information offering a rich choice. Here, we investigate the most appropriate study design in phase I for several hypothetical mechanisms of action of a study drug. Five drug effects in diabetes therapeutic areas were investigated using six study designs. Power to detect drug effect was assessed using the likelihood ratio test, whereas precision and accuracy of the quantification of drug effect was assessed using stochastic simulation and estimations. An overall summary was developed to aid designing the studies of antihyperglycemic drug development using model-based analysis. This guidance is to be used when the integrated glucose insulin model is used, involving the investigated drug mechanisms of action. Correction in: CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY, Volume: 8, Issue: 7, Pages: 520-520, DOI: 10.1002/psp4.12391

The relationships between exposure, biomarkers (vascular endothelial growth factor (VEGF), soluble VEGF receptors (sVEGFR)-1, -2, -3, and soluble stem cell factor receptor (sKIT)), tumor sum of longest diameters (SLD), diastolic blood pressure (dBP), and overall survival (OS) were investigated in a modeling framework. The dataset included 64 metastatic renal cell carcinoma patients (mRCC) treated with oral axitinib. Biomarker timecourses were described by indirect response (IDR) models where axitinib inhibits sVEGFR-1, -2, and -3 production, and VEGF degradation. No effect was identified on sKIT. A tumor model using sVEGFR-3 dynamics as driver predicted SLD data well. An IDR model, with axitinib exposure stimulating the response, characterized dBP increase. In a time-to-event model the SLD timecourse predicted OS better than exposure, biomarker- or dBP-related metrics. This type of framework can be used to relate pharmacokinetics, efficacy, and safety to long-term clinical outcome in mRCC patients treated with VEGFR inhibitors. (ClinicalTrial.gov identifier NCT00569946.).

Count data, or number of events per time interval, are discrete data arising from repeated time to event observations. Their mean count, or piecewise constant event rate, can be evaluated by discrete probability distributions from the Poisson model family. Clinical trial data characterization often involves population count analysis. This tutorial presents the basics and diagnostics of count modeling and simulation in the context of pharmacometrics. Consideration is given to overdispersion, underdispersion, autocorrelation, and inhomogeneity.