Hierarchical cluster analysis and nonlinear mixed-effects modelling for candidate biomarker detection in preclinical models of cancer.

BACKGROUND: Preclinical models of cancer can be of translational benefit when assessing how different biomarkers are regulated in response to particular treatments. Detection of molecular biomarkers in preclinical models of cancer is difficult due inter-animal variability in responses, combined with limited accessibility of longitudinal data. METHODS: Nonlinear mixed-effects modelling (NLME) was used to analyse tumour growth data based on expected tumour growth rates observed 7 days after initial doses (DD7) of Radiotherapy (RT) and Combination of RT with DNA Damage Response Inhibitors (DDRi). Cox regression was performed to confirm an association between DD7 and survival. Hierarchical Cluster Analysis (HCA) was then used to identify candidate biomarkers impacting responses to RT and RT/DDRi and these were validated using NLME. RESULTS: Cox regression confirmed significant associations between DD7 and survival. HCA of RT treated samples, combined with NLME confirmed significant associations between DD7 and Cluster specific CD8+ Ki67 MFI, as well as DD7 and cluster specific Natural Killer cell density in RT treated mice. CONCLUSION: Application of NLME, as well as HCA of candidate biomarkers may provide additional avenues to assess the effect of RT in MC38 syngeneic tumour models. Additional studies would need to be conducted to confirm association between DD7 and biomarkers in RT/DDRi treated mice.

Covariates in population pharmacokinetic studies of critically ill adults receiving ß-lactam antimicrobials: a systematic review and narrative synthesis.

Introduction: Population pharmacokinetic studies of ß-lactam antimicrobials in critically ill patients derive models that inform their dosing. In non-linear mixed-effects modelling, covariates are often used to improve model fit and explain variability. We aimed to investigate which covariates are most commonly assessed and which are found to be significant, along with global patterns of publication. Methods: We conducted a systematic review, searching MEDLINE, Embase, CENTRAL and Web of Science on 01 March 2023, including studies of critically ill adults receiving ß-lactam antimicrobials who underwent blood sampling for population pharmacokinetic studies. We extracted and categorized all reported covariates and assessed reporting quality using the ClinPK checklist. Results: Our search identified 151 studies with 6018 participants. Most studies reported observational cohorts (120 studies, 80%), with the majority conducted in high-income settings (136 studies, 90%). Of the 1083 identified covariate instances, 237 were unique; the most common categories were patient characteristics (n = 404), biomarkers (n = 206) and physiological parameters (n = 163). Only seven distinct commonly reported covariates (CLCR, weight, glomerular filtration rate, diuresis, need for renal replacement, serum albumin and C-reactive protein) were significant more than 20% of the time. Conclusions: Covariates are most commonly chosen based on biological plausibility, with patient characteristics and biomarkers the most frequently investigated. We developed an openly accessible database of reported covariates to aid investigators with covariate selection when designing population pharmacokinetic studies. Novel covariates, such as sepsis subphenotypes, have not been explored yet, leaving a research gap for future work.

Mixed effects modeling of radiotherapy in combination with immune checkpoint blockade or inhibitors of the DNA damage response pathway.

Dosage optimization to maximize efficacy and minimize toxicity is a potential issue when administering radiotherapy (RT) in combination with immune checkpoint blockade (ICB) or inhibitors of the DNA Damage Response Pathway (DDRi) in the clinic. Preclinical models and mathematical modeling can help identify ideal dosage schedules to observe beneficial effects of a tri-therapy. The aim of this study is to describe a mathematical model to capture the impact of RT in combination with inhibitors of the DNA Damage Response Pathway or blockade of the immune checkpoint protein - programmed death ligand 1 (PD-L1). This model describes how RT mediated activation of antigen presenting cells can induce an increase in cytolytic T cells capable of targeting tumor cells, and how combination drugs can potentiate the immune response by inhibiting the rate of T cell exhaustion. The model was fitted using preclinical data, where MC38 tumors were treated in vivo with RT alone or in combination with anti-PD-L1 as well as with either olaparib or the ataxia telangiectasia mutated (ATM) inhibitor-AZD0156. The model successfully described the observed data and goodness-of-fit, using visual predictive checks also confirmed a successful internal model validation for each treatment modality. The results demonstrated that the anti-PD-L1 effect in combination with RT was maximal in vivo and any additional benefit of DDRi at the given dosage and schedule used was undetectable. Model fit results indicated AZD0156 to be a more potent DDRi than olaparib. Simulations of alternative doses indicated that reducing efficacy of anti-PD-L1 by 68% would potentially provide evidence for a benefit of ATM inhibition in combination with ICB and increase the relative efficacy of tri-therapy.

Exploring the Boundaries for In Vitro-In Vivo Extrapolation: Use of Isolated Rat Hepatocytes in Co-culture and Impact of Albumin Binding Properties in the Prediction of Clearance of Various Drug Types.

Prediction of hepatic clearance of drugs (via uptake or metabolism) from in vitro systems continues to be problematic, particularly when plasma protein binding is high. The following work explores simultaneous assessment of both clearance processes, focusing on a commercial hepatocyte-fibroblast co-culture system (HµREL) over a 24-hour period using six probe drugs (ranging in metabolic and transporter clearance and low-to-high plasma protein binding). A rat hepatocyte co-culture assay was established using drug depletion (measuring both medium and total concentrations) and cell uptake kinetic analysis, both in the presence and absence of plasma protein (1% bovine serum albumin). Secretion of endogenous albumin was monitored as a marker of viability, and this reached 0.004% in incubations (at a rate similar to in vivo synthesis). Binding to stromal cells was substantial and required appropriate correction factors. Drug concentration-time courses were analyzed both by conventional methods and a mechanistic cell model prior to in vivo extrapolation. Clearance assayed by drug depletion in conventional suspended rat hepatocytes provided a benchmark to evaluate co-culture value. Addition of albumin appeared to improve predictions for some compounds (where fraction unbound in the medium is less than 0.1); however, for high-binding drugs, albumin significantly limited quantification and thus predictions. Overall, these results highlight ongoing challenges concerning in vitro hepatocyte system complexity and limitations of practical expediency. Considering this, more reliable measurement of hepatically cleared compounds seems possible through judicious use of available hepatocyte systems, including co-culture systems, as described herein; this would include those compounds with low metabolic turnover but high active uptake clearance. SIGNIFICANCE STATEMENT: Co-culture systems offer a more advanced tool than standard hepatocytes, with the ability to be cultured for longer periods of time, yet their potential as an in vitro tool has not been extensively assessed. We evaluate the strengths and limitations of the HµREL system using six drugs representing various metabolic and transporter-mediated clearance pathways with various degrees of albumin binding. Studies in the presence/absence of albumin allow in vitro-in vivo extrapolation and a framework to maximize their utility.

Radiation in Combination with Immune Checkpoint Blockade and DNA Damage Response Inhibitors in Mice: Dosage Optimization in MC38 Syngeneic Tumors via Modelling and Simulation.

Clinical trials assessing the impact of radiotherapy (RT) in combination with DNA damage response pathway inhibitors (DDRis) and/or immune checkpoint blockade are currently ongoing. However, current methods for optimizing dosage and schedule are limited. A mathematical model was developed to capture the impacts of RT in combination with DDRi and/or anti-PD-L1 [immune checkpoint inhibitor (ICI)] on tumor immune interactions in the MC38 syngeneic tumor model. The model was fitted to datasets that assessed the impact of RT in combination with the DNA protein kinase inhibitor (DNAPKi) AZD7648. The model was further fitted to datasets from studies that were used to assess both RT/ICI combinations as well as RT/ICI combinations followed by concurrent administration of the poly ADP ribose polymerase inhibitor (PARPi) olaparib. Nonlinear mixed-effects modeling was performed followed by internal validation with visual predictive checks (VPC). Simulations of alternative dosage regimens and scheduling were performed to identify optimal candidate dosage regimens of RT/DNAPKi and RT/PARPi/ICI. Model fits and VPCs confirmed a successful internal validation for both datasets and demonstrated very small differences in the median, lower, and upper percentile values of tumor diameters between RT/ICI and RT/PARPi/ICI, which indicated that the triple combination of RT/PARPi/ICI at the given dosage and schedule does not provide additional benefit compared with ICI in combination with RT. Simulation of alternative dosage regimens indicated that lowering the dosage of ICI to between 2 and 4 mg/kg could induce similar benefits to the full dosage regimen, which could be of translational benefit. SIGNIFICANCE STATEMENT: This work provides a mixed-effects model framework to quantify the effects of combination radiotherapy/DNA damage response pathway inhibitors/immune checkpoint inhibitors in preclinical tumor models and identify optimal dosage regimens, which could be of translational benefit.

Estimation of FMO3 Ontogeny by Mechanistic Population Pharmacokinetic Modelling of Risdiplam and Its Impact on Drug-Drug Interactions in Children.

BACKGROUND AND OBJECTIVE: Spinal muscular atrophy (SMA) is a progressive neuromuscular disease caused by insufficient levels of survival motor neuron (SMN) protein. Risdiplam (EvrysdiTM) increases SMN protein and is approved for the treatment of SMA. Risdiplam has high oral bioavailability and is primarily eliminated through hepatic metabolism by flavin-containing monooxygenase3 (FMO3) and cytochrome P450 (CYP) 3A, by 75% and 20%, respectively. While the FMO3 ontogeny is critical input data for the prediction of risdiplam pharmacokinetics (PK) in children, it was mostly studied in vitro, and robust in vivo FMO3 ontogeny is currently lacking. We derived in vivo FMO3 ontogeny by mechanistic population PK modelling of risdiplam and investigated its impact on drug-drug interactions in children. METHODS: Population and physiologically based PK (PPK and PBPK) modelling conducted during the development of risdiplam were integrated into a mechanistic PPK (Mech-PPK) model to estimate in vivo FMO3 ontogeny. A total of 10,205 risdiplam plasma concentration-time data from 525 subjects aged 2 months-61 years were included. Six different structural models were examined to describe the in vivo FMO3 ontogeny. Impact of the newly estimated FMO3 ontogeny on predictions of drug-drug interaction (DDI) in children was investigated by simulations for dual CYP3A-FMO3 substrates including risdiplam and theoretical substrates covering a range of metabolic fractions (fm) of CYP3A and FMO3 (fmCYP3A:fmFMO3 = 10%:90%, 50%:50%, 90%:10%). RESULTS: All six models consistently predicted higher FMO3 expression/activity in children, reaching a maximum at the age of 2 years with an approximately threefold difference compared with adults. Different trajectories of FMO3 ontogeny in infants < 4 months of age were predicted by the six models, likely due to limited observations for this age range. Use of this  in vivo FMO3 ontogeny function improved prediction of risdiplam PK in children compared to in vitro FMO3 ontogeny functions. The simulations of theoretical dual CYP3A-FMO3 substrates predicted comparable or decreased CYP3A-victim DDI propensity in children compared to adults across the range of fm values. Refinement of FMO3 ontogeny in the risdiplam model had no impact on the previously predicted low CYP3A-victim or -perpetrator DDI risk of risdiplam in children. CONCLUSION: Mech-PPK modelling successfully estimated in vivo FMO3 ontogeny from risdiplam data collected from 525 subjects aged 2 months-61 years. To our knowledge, this is the first investigation of in vivo FMO3 ontogeny by population approach using comprehensive data covering a wide age range. Derivation of a robust in vivo FMO3 ontogeny function has significant implications on the prospective prediction of PK and DDI in children for other FMO3 substrates in the future, as illustrated in the current study for FMO3 and/or dual CYP3A-FMO3 substrates. CLINICAL TRIAL REGISTRY NUMBERS: NCT02633709, NCT03032172, NCT02908685, NCT02913482, NCT03988907.

Use of In Vivo Imaging and Physiologically-Based Kinetic Modelling to Predict Hepatic Transporter Mediated Drug-Drug Interactions in Rats.

Gadoxetate, a magnetic resonance imaging (MRI) contrast agent, is a substrate of organic-anion-transporting polypeptide 1B1 and multidrug resistance-associated protein 2. Six drugs, with varying degrees of transporter inhibition, were used to assess gadoxetate dynamic contrast enhanced MRI biomarkers for transporter inhibition in rats. Prospective prediction of changes in gadoxetate systemic and liver AUC (AUCR), resulting from transporter modulation, were performed by physiologically-based pharmacokinetic (PBPK) modelling. A tracer-kinetic model was used to estimate rate constants for hepatic uptake (khe), and biliary excretion (kbh). The observed median fold-decreases in gadoxetate liver AUC were 3.8- and 1.5-fold for ciclosporin and rifampicin, respectively. Ketoconazole unexpectedly decreased systemic and liver gadoxetate AUCs; the remaining drugs investigated (asunaprevir, bosentan, and pioglitazone) caused marginal changes. Ciclosporin decreased gadoxetate khe and kbh by 3.78 and 0.09 mL/min/mL, while decreases for rifampicin were 7.20 and 0.07 mL/min/mL, respectively. The relative decrease in khe (e.g., 96% for ciclosporin) was similar to PBPK-predicted inhibition of uptake (97-98%). PBPK modelling correctly predicted changes in gadoxetate systemic AUCR, whereas underprediction of decreases in liver AUCs was evident. The current study illustrates the modelling framework and integration of liver imaging data, PBPK, and tracer-kinetic models for prospective quantification of hepatic transporter-mediated DDI in humans.

Development of a sensitive biochemical assay for the detection of tofacitinib adherence.

Rheumatoid arthritis (RA) is a chronic autoimmune inflammatory disease. Tofacitinib is a Janus Kinase inhibitor licensed for the treatment of RA that, unlike biologic anti-rheumatic drugs, is administered orally, but studies of long-term treatment adherence rates are lacking. The measurement of adherence, however, is challenging and there is currently no gold standard test for adherence. Here, we developed a novel HPLC MS/MS assay for the quantification of tofacitinib. The assay demonstrated a LLOQ for tofacitinib of 0.1 ng ml-1, within run accuracy was 81-85% at LLOQ and 91-107% at all other levels. To investigate the ability of the assay to detect adherence, tofacitinib was measured in a random selection of serum samples (n = 10) of tofacitinib treated RA patients who self-reported adherent behaviour. The assay measured tofacitinib in all samples above the LLOQ demonstrating the potential of the assay to sensitively measure biochemical adherence in real-world patient samples. This method for detection of adherence has the potential to be a more objective measure that could be used in the future in the clinic but will require further studies to explore factors that may influence measurement of drug levels, such as clinical characteristics of patients.

Investigation of simplified physiologically-based pharmacokinetic models in rat and human.

Whole-body physiologically-based pharmacokinetic (PBPK) models have many applications in drug research and development. It is often necessary to inform these models with animal or clinical data, updating model parameters, and making the model more predictive for future applications. This provides an opportunity and a challenge given the large number of parameters of such models. The aim of this work was to propose new mechanistic model structures with reduced complexity allowing for parameter optimization. These models were evaluated for their ability to estimate realistic values for unbound tissue to plasma partition coefficients (Kpu) and simulate observed pharmacokinetic (PK) data. Two approaches are presented: using either established kinetic lumping methods based on tissue time constants (drug-dependent) or a novel clustering analysis to identify tissues sharing common Kpu values or Kpu scalars based on similarities of tissue composition (drug-independent). PBPK models derived from these approaches were assessed using PK data of diazepam in rats and humans. Although the clustering analysis produced minor differences in tissue grouping depending on the method used, two larger groups of tissues emerged. One including the kidneys, liver, spleen, gut, heart, and lungs, and another including bone, brain, muscle, and pancreas whereas adipose and skin were generally considered distinct. Overall, a subdivision into four tissue groups appeared most physiologically relevant in terms of tissue composition. Several models were found to have similar abilities to describe diazepam i.v. data as empirical models. Comparability of estimated Kpus to experimental Kpu values for diazepam was one criterion for selecting the appropriate PK model structure.

A "middle-out approach" for the prediction of human drug disposition from preclinical data using simplified physiologically based pharmacokinetic (PBPK) models.

Simplified physiologically based pharmacokinetic (PBPK) models using estimated tissue-to-unbound plasma partition coefficients (Kpus) were previously investigated by fitting them to in vivo pharmacokinetic (PK) data. After optimization with preclinical data, the performance of these models for extrapolation of distribution kinetics to human were evaluated to determine the best approach for the prediction of human drug disposition and volume of distribution (Vss) using PBPK modeling. Three lipophilic bases were tested (diazepam, midazolam, and basmisanil) for which intravenous PK data were available in rat, monkey, and human. The models with Kpu scalars using k-means clustering were generally the best for fitting data in the preclinical species and gave plausible Kpu values. Extrapolations of plasma concentrations for diazepam and midazolam using these models and parameters obtained were consistent with the observed clinical data. For diazepam and midazolam, the human predictions of Vss after optimization in rats and monkeys were better compared with the Vss estimated from the traditional PBPK modeling approach (varying from 1.1 to 3.1 vs. 3.7-fold error). For basmisanil, the sparse preclinical data available could have affected the model performance for fitting and the subsequent extrapolation to human. Overall, this work provides a rational strategy to predict human drug distribution using preclinical PK data within the PBPK modeling strategy.

Coproporphyrin I as an Endogenous Biomarker to Detect Reduced OATP1B Activity and Shift in Elimination Route in Chronic Kidney Disease.

Coproporphyrin I (CPI) is an endogenous biomarker of organic anion transporting polypeptide 1B transporter (OATP1B). CPI plasma baseline was reported to increase with severity of chronic kidney disease (CKD). Further, ratio of CPI area under the plasma concentration-time curve (AUCR) in the presence/absence of OATP1B inhibitor rifampin was higher in patients with CKD compared with healthy participants, in contrast to pitavastatin (a clinical OATP1B probe). This study investigated mechanism(s) contributing to altered CPI baseline in patients with CKD by extending a previously developed physiologically-based pharmacokinetic (PBPK) model to this patient population. CKD-related covariates were evaluated in a stepwise manner on CPI fraction unbound in plasma (fu,p ), OATP1B-mediated hepatic uptake clearance (CLactive ), renal clearance (CLR ), and endogenous synthesis (ksyn ). The CPI model successfully recovered increased baseline and rifampin-mediated AUCR in patients with CKD by accounting for the following disease-related changes: 13% increase in fu,p , 29% and 39% decrease in CLactive in mild and moderate to severe CKD, respectively, decrease in CLR proportional to decline in glomerular filtration rate, and 27% decrease in ksyn in severe CKD. Almost complete decline in CPI renal elimination in severe CKD increased its fraction transported by OATP1B, rationalizing differences in the CPI-rifampin interaction observed between healthy participants and patients with CKD. In conclusion, mechanistic modeling performed here supports CKD-related decrease in OATP1B function to inform prospective PBPK modeling of OATP1B-mediated drug-drug interaction in these patients. Monitoring of CPI allows detection of CKD-drug interaction risk for OATP1B drugs with combined hepatic and renal elimination which may be underestimated by extrapolating the interaction risk based on pitavastatin data in healthy participants.

In Vitro to In Vivo Extrapolation Linked to Physiologically Based Pharmacokinetic Models for Assessing the Brain Drug Disposition.

Drug development for the central nervous system (CNS) is a complex endeavour with low success rates, as the structural complexity of the brain and specifically the blood-brain barrier (BBB) poses tremendous challenges. Several in vitro brain systems have been evaluated, but the ultimate use of these data in terms of translation to human brain concentration profiles remains to be fully developed. Thus, linking up in vitro-to-in vivo extrapolation (IVIVE) strategies to physiologically based pharmacokinetic (PBPK) models of brain is a useful effort that allows better prediction of drug concentrations in CNS components. Such models may overcome some known aspects of inter-species differences in CNS drug disposition. Required physiological (i.e. systems) parameters in the model are derived from quantitative values in each organ. However, due to the inability to directly measure brain concentrations in humans, compound-specific (drug) parameters are often obtained from in silico or in vitro studies. Such data are translated through IVIVE which could be also applied to preclinical in vivo observations. In such exercises, the limitations of the assays and inter-species differences should be adequately understood in order to verify these predictions with the observed concentration data. This report summarizes the state of IVIVE-PBPK-linked models and discusses shortcomings and areas of further research for better prediction of CNS drug disposition. Graphical abstract.

Evaluation of area under the concentration curve adjusted by the terminal-phase as a metric to reduce the impact of variability in bioequivalence testing.

AIM: To quantify the utility of a terminal-phase adjusted area under the concentration curve method in increasing the probability of a correct and conclusive outcome of a bioequivalence (BE) trial for highly variable drugs when clearance (CL) varies more than the volume of distribution (V). METHODS: Data from a large population of subjects were generated with variability in CL and V, and used to simulate a two-period, two-sequence crossover BE trial. The 90% confidence interval for formulation comparison was determined following BE assessment using the area under the concentration curve (AUC) ratio test, and the proposed terminal-phase adjusted AUC ratio test. An outcome of bioequivalent, nonbioequivalent or inconclusive was then assigned according to predefined BE limits. RESULTS: When CL varied more than V, the proposed approach enhanced the probability of correctly assigning bioequivalent or nonbioequivalent and reduced the risk of an inconclusive trial. For a hypothetical drug with between-subject variability of 35% for CL and 10% for V, when the true test-reference ratio of bioavailability was 1.15, a crossover study of n = 14 subjects analysed by the proposed method would have 80% or 20% probability of claiming bioequivalent or nonbioequivalent, compared to 22%, 46% or 32% probability of claiming bioequivalent, nonbioequivalent or inconclusive using the standard AUC ratio test. CONCLUSIONS: The terminal-phase adjusted AUC ratio test represents a simple and readily applicable approach to enhance the BE assessment of drug products when CL varies more than V.

Model-Based Drug-Drug Interaction Extrapolation Strategy From Adults to Children: Risdiplam in Pediatric Patients With Spinal Muscular Atrophy.

Risdiplam (Evrysdi) improves motor neuron function in patients with spinal muscular atrophy (SMA) and has been approved for the treatment of patients ≥2 months old. Risdiplam exhibits time-dependent inhibition of cytochrome P450 (CYP) 3A in vitro. While many pediatric patients receive risdiplam, a drug-drug interaction (DDI) study in pediatric patients with SMA was not feasible. Therefore, a novel physiologically-based pharmacokinetic (PBPK) model-based strategy was proposed to extrapolate DDI risk from healthy adults to children with SMA in an iterative manner. A clinical DDI study was performed in healthy adults at relevant risdiplam exposures observed in children. Risdiplam caused an 1.11-fold increase in the ratio of midazolam area under the curve with and without risdiplam (AUCR)), suggesting an 18-fold lower in vivo CYP3A inactivation constant compared with the in vitro value. A pediatric PBPK model for risdiplam was validated with independent data and combined with a validated midazolam pediatric PBPK model to extrapolate DDI from adults to pediatric patients with SMA. The impact of selected intestinal and hepatic CYP3A ontogenies on the DDI susceptibility in children relative to adults was investigated. The PBPK analysis suggests that primary CYP3A inhibition by risdiplam occurs in the intestine rather than the liver. The PBPK-predicted risdiplam CYP3A inhibition risk in pediatric patients with SMA aged 2 months-18 years was negligible (midazolam AUCR of 1.09-1.18) and included in the US prescribing information of risdiplam. Comprehensive evaluation of the sensitivity of predicted CYP3A DDI on selected intestinal and hepatic CYP3A ontogeny functions, together with PBPK model-based strategy proposed here, aim to guide and facilitate DDI extrapolations in pediatric populations.

Reduced physiologically-based pharmacokinetic model of dabigatran etexilate-dabigatran and its application for prediction of intestinal P-gp-mediated drug-drug interactions.

BACKGROUND: Dabigatran etexilate (DABE) has been suggested as a clinical probe for intestinal P-glycoprotein (P-gp)-mediated drug-drug interaction (DDI) studies and, as an alternative to digoxin. Clinical DDI data with various P-gp inhibitors demonstrated a dose-dependent inhibition of P-gp with DABE. The aims of this study were to develop a joint DABE (prodrug)-dabigatran reduced physiologically-based-pharmacokinetic (PBPK) model and to evaluate its ability to predict differences in P-gp DDI magnitude between a microdose and a therapeutic dose of DABE. METHODS: A joint DABE-dabigatran PBPK model was developed with a mechanistic intestinal model accounting for the regional P-gp distribution in the gastrointestinal tract. Model input parameters were estimated using DABE and dabigatran pharmacokinetic (PK) clinical data obtained after administration of DABE alone or with a strong P-gp inhibitor, itraconazole, and over a wide range of DABE doses (from 375 µg to 400 mg). Subsequently, the model was used to predict extent of DDI with additional P-gp inhibitors and with different DABE doses. RESULTS: The reduced DABE-dabigatran PBPK model successfully described plasma concentrations of both prodrug and metabolite following administration of DABE at different dose levels and when co-administered with itraconazole. The model was able to capture the dose dependency in P-gp mediated DDI. Predicted magnitude of itraconazole P-gp DDI was higher at the microdose (predicted vs. observed median fold-increase in AUC+inh/AUCcontrol (min-max) = 5.88 (4.29-7.93) vs. 6.92 (4.96-9.66) ) compared to the therapeutic dose (predicted median fold-increase in AUC+inh/AUCcontrol = 3.48 (2.37-4.84) ). In addition, the reduced DABE-dabigatran PBPK model predicted successfully the extent of DDI with verapamil and clarithromycin as P-gp inhibitors. Model-based simulations of dose staggering predicted the maximum inhibition of P-gp when DABE microdose was concomitantly administered with itraconazole solution; simulations also highlighted dosing intervals required to minimise the DDI risk depending on the DABE dose administered (microdose vs. therapeutic). CONCLUSIONS: This study provides a modelling framework for the evaluation of P-gp inhibitory potential of new molecular entities using DABE as a clinical probe. Simulations of dose staggering and regional differences in the extent of intestinal P-gp inhibition for DABE microdose and therapeutic dose provide model-based guidance for design of prospective clinical P-gp DDI studies.

Physiologically Based Pharmacokinetic Modeling of Transporter-Mediated Hepatic Disposition of Imaging Biomarker Gadoxetate in Rats.

Physiologically based pharmacokinetic (PBPK) models are increasingly used in drug development to simulate changes in both systemic and tissue exposures that arise as a result of changes in enzyme and/or transporter activity. Verification of these model-based simulations of tissue exposure is challenging in the case of transporter-mediated drug-drug interactions (tDDI), in particular as these may lead to differential effects on substrate exposure in plasma and tissues/organs of interest. Gadoxetate, a promising magnetic resonance imaging (MRI) contrast agent, is a substrate of organic-anion-transporting polypeptide 1B1 (OATP1B1) and multidrug resistance-associated protein 2 (MRP2). In this study, we developed a gadoxetate PBPK model and explored the use of liver-imaging data to achieve and refine in vitro-in vivo extrapolation (IVIVE) of gadoxetate hepatic transporter kinetic data. In addition, PBPK modeling was used to investigate gadoxetate hepatic tDDI with rifampicin i.v. 10 mg/kg. In vivo dynamic contrast-enhanced (DCE) MRI data of gadoxetate in rat blood, spleen, and liver were used in this analysis. Gadoxetate in vitro uptake kinetic data were generated in plated rat hepatocytes. Mean (%CV) in vitro hepatocyte uptake unbound Michaelis-Menten constant (Km,u) of gadoxetate was 106 µM (17%) (n = 4 rats), and active saturable uptake accounted for 94% of total uptake into hepatocytes. PBPK-IVIVE of these data (bottom-up approach) captured reasonably systemic exposure, but underestimated the in vivo gadoxetate DCE-MRI profiles and elimination from the liver. Therefore, in vivo rat DCE-MRI liver data were subsequently used to refine gadoxetate transporter kinetic parameters in the PBPK model (top-down approach). Active uptake into the hepatocytes refined by the liver-imaging data was one order of magnitude higher than the one predicted by the IVIVE approach. Finally, the PBPK model was fitted to the gadoxetate DCE-MRI data (blood, spleen, and liver) obtained with and without coadministered rifampicin. Rifampicin was estimated to inhibit active uptake transport of gadoxetate into the liver by 96%. The current analysis highlighted the importance of gadoxetate liver data for PBPK model refinement, which was not feasible when using the blood data alone, as is common in PBPK modeling applications. The results of our study demonstrate the utility of organ-imaging data in evaluating and refining PBPK transporter IVIVE to support the subsequent model use for quantitative evaluation of hepatic tDDI.

Investigating the impact of target lesion selection on drug effect evaluation and tumour growth rate determination using tumour growth inhibition models: Example of malignant pleural mesothelioma patients treated with cisplatin alone or in combination with pemetrexed.

In the last update of the RECIST criteria in 2009, it was proposed that the number of target lesions to be followed over time for response-to-treatment assessment be reduced from 10 to 5 lesions maximum, with up to 2 per organ. We explored the impact of reducing the number of target lesion on the assessment of drug effect in a randomised phase III clinical trial using a tumour growth inhibition (TGI) model. Tumour size measurements from 441 (out of 456) patients were used to build two datasets for which observations were the sum of longest diameters of all measurable lesions (ALL dataset) or following the RECIST 1.1 recommendations (R1.1 dataset). TGI models incorporating a categorical covariate for treatment group or a pharmacokinetic metric (i.e. dose; simulated area under the curve) were used to describe the longitudinal tumour size kinetics. Drug exposure was not superior to treatment group at describing drug effect. ALL and R1.1 individual estimates of drug effect appeared to be strongly correlated (r2=0.88). Including pharmacokinetic metrics in TGI models should be conducted carefully when no pharmacokinetic samples are available. Reducing the number of target lesion did not seem to compromise the determination of drug effect using TGI models.

Population pharmacokinetic modeling and simulation to support qualification of pyridoxic acid as endogenous biomarker of OAT1/3 renal transporters.

Renal clearance of many drugs is mediated by renal organic anion transporters OAT1/3 and inhibition of these transporters may lead to drug-drug interactions (DDIs). Pyridoxic acid (PDA) and homovanillic acid (HVA) were indicated as potential biomarkers of OAT1/3. The objective of this study was to develop a population pharmacokinetic model for PDA and HVA to support biomarker qualification. Simultaneous fitting of biomarker plasma and urine data in the presence and absence of potent OAT1/3 inhibitor (probenecid, 500 mg every 6 h) was performed. The impact of study design (multiple vs. single dose of OAT1/3 inhibitor) and ability to detect interactions in the presence of weak/moderate OAT1/3 inhibitors was investigated, together with corresponding power calculations. The population models developed successfully described biomarker baseline and PDA/HVA OAT1/3-mediated interaction data. No prominent effect of circadian rhythm on PDA and HVA individual baseline levels was evident. Renal elimination contributed greater than 80% to total clearance of both endogenous biomarkers investigated. Estimated probenecid unbound in vivo OAT inhibitory constant was up to 6.4-fold lower than in vitro values obtained with PDA as a probe. The PDA model was successfully verified against independent literature reported datasets. No significant difference in power of DDI detection was found between multiple and single dose study design when using the same total daily dose of 2000 mg probenecid. Model-based simulations and power calculations confirmed sensitivity and robustness of plasma PDA data to identify weak, moderate, and strong OAT1/3 inhibitors in an adequately powered clinical study to support optimal design of prospective clinical OAT1/3 interaction studies.

Model-Informed Precision Dosing: Background, Requirements, Validation, Implementation, and Forward Trajectory of Individualizing Drug Therapy.

Model-informed precision dosing (MIPD) has become synonymous with modern approaches for individualizing drug therapy, in which the characteristics of each patient are considered as opposed to applying a one-size-fits-all alternative. This review provides a brief account of the current knowledge, practices, and opinions on MIPD while defining an achievable vision for MIPD in clinical care based on available evidence. We begin with a historical perspective on variability in dose requirements and then discuss technical aspects of MIPD, including the need for clinical decision support tools, practical validation, and implementation of MIPD in health care. We also discuss novel ways to characterize patient variability beyond the common perceptions of genetic control. Finally, we address current debates on MIPD from the perspectives of the new drug development, health economics, and drug regulations.

Does "Birth" as an Event Impact Maturation Trajectory of Renal Clearance via Glomerular Filtration? Reexamining Data in Preterm and Full-Term Neonates by Avoiding the Creatinine Bias.

Glomerular filtration rate (GFR) is an important measure of renal function. Various models for its maturation have recently been compared; however, these have used markers, which are subject to different renal elimination processes. Inulin clearance data (a purer probe of GFR) collected from the literature were used to determine age-related changes in GFR aspects of renal drug excretion in pediatrics. An ontogeny model was derived using a best-fit model with various combinations of covariates such as postnatal age, gestational age at birth, and body weight. The model was applied to the prediction of systemic clearance of amikacin, gentamicin, vancomycin, and gadobutrol. During neonatal life, GFR increased as a function of both gestational age at birth and postnatal age, hence implying an impact of birth and a discrepancy in GFR for neonates with the same postmenstrual age depending on gestational age at birth (ie, neonates who were outside the womb longer had higher GFR, on average). The difference in GFR between pre-term and full-term neonates with the same postmenstrual age was negligible from beyond 1.25 years. Considering both postnatal age and gestational age at birth in GFR ontogeny models is important because postmenstrual age alone ignores the impact of birth. Most GFR models use covariates of body size in addition to age. Therefore, prediction from these models will also depend on the change in anthropometric characteristics with age. The latter may not be similar in various ethnic groups, and this makes the head-to-head comparison of models very challenging.