Physiologically Based Pharmacokinetic Models of Probenecid and Furosemide to Predict Transporter Mediated Drug-Drug Interactions.

PURPOSE: To provide whole-body physiologically based pharmacokinetic (PBPK) models of the potent clinical organic anion transporter (OAT) inhibitor probenecid and the clinical OAT victim drug furosemide for their application in transporter-based drug-drug interaction (DDI) modeling. METHODS: PBPK models of probenecid and furosemide were developed in PK-Sim®. Drug-dependent parameters and plasma concentration-time profiles following intravenous and oral probenecid and furosemide administration were gathered from literature and used for model development. For model evaluation, plasma concentration-time profiles, areas under the plasma concentration-time curve (AUC) and peak plasma concentrations (Cmax) were predicted and compared to observed data. In addition, the models were applied to predict the outcome of clinical DDI studies. RESULTS: The developed models accurately describe the reported plasma concentrations of 27 clinical probenecid studies and of 42 studies using furosemide. Furthermore, application of these models to predict the probenecid-furosemide and probenecid-rifampicin DDIs demonstrates their good performance, with 6/7 of the predicted DDI AUC ratios and 4/5 of the predicted DDI Cmax ratios within 1.25-fold of the observed values, and all predicted DDI AUC and Cmax ratios within 2.0-fold. CONCLUSIONS: Whole-body PBPK models of probenecid and furosemide were built and evaluated, providing useful tools to support the investigation of transporter mediated DDIs.

Data Digitizing: Accurate and Precise Data Extraction for Quantitative Systems Pharmacology and Physiologically-Based Pharmacokinetic Modeling.

In quantitative systems pharmacology (QSP) and physiologically-based pharmacokinetic (PBPK) modeling, data digitizing is a valuable tool to extract numerical information from published data presented as graphs. To quantify their relevance, a literature search revealed a remarkable mean increase of 16% per year in publications citing digitizing software together with QSP or PBPK. Accuracy, precision, confounder influence, and variability were investigated using scaled median symmetric accuracy (ζ), thus finding excellent accuracy (mean ζ = 0.99%). Although significant, no relevant confounders were found (mean ζ ± SD circles = 0.69% ± 0.68% vs. triangles = 1.3% ± 0.62%). Analysis of 181 literature peak plasma concentration values revealed a considerable discrepancy between reported and post hoc digitized data with 85% having ζ > 5%. Our findings suggest that data digitizing is precise and important. However, because the greatest pitfall comes from pre-existing errors, we recommend always making published data available as raw values.

Comprehensive Parent-Metabolite PBPK/PD Modeling Insights into Nicotine Replacement Therapy Strategies.

BACKGROUND: Nicotine, the pharmacologically active substance in both tobacco and many electronic cigarette (e-cigarette) liquids, is responsible for the addiction that sustains cigarette smoking. With 8 million deaths worldwide annually, smoking remains one of the major causes of disability and premature death. However, nicotine also plays an important role in smoking cessation strategies. OBJECTIVES: The aim of this study was to develop a comprehensive, whole-body, physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model of nicotine and its major metabolite cotinine, covering various routes of nicotine administration, and to simulate nicotine brain tissue concentrations after the use of combustible cigarettes, e-cigarettes, nicotine gums, and nicotine patches. METHODS: A parent-metabolite, PBPK/PD model of nicotine for a non-smoking and a smoking population was developed using 91 plasma and brain tissue concentration-time profiles and 11 heart rate profiles. Among others, cytochrome P450 (CYP) 2A6 and 2B6 enzymes were implemented, including kinetics for CYP2A6 poor metabolizers. RESULTS: The model is able to precisely describe and predict both nicotine plasma and brain tissue concentrations, cotinine plasma concentrations, and heart rate profiles. 100% of the predicted area under the concentration-time curve (AUC) and maximum concentration (Cmax) values meet the twofold acceptance criterion with overall geometric mean fold errors of 1.12 and 1.15, respectively. The administration of combustible cigarettes, e-cigarettes, nicotine patches, and nicotine gums was successfully implemented in the model and used to identify differences in steady-state nicotine brain tissue concentration patterns. CONCLUSIONS: Our PBPK/PD model may be helpful in further investigations of nicotine dependence and smoking cessation strategies. As the model represents the first nicotine PBPK/PD model predicting nicotine concentration and heart rate profiles after the use of e-cigarettes, it could also contribute to a better understanding of the recent increase in youth e-cigarette use.

Physiologically-Based Pharmacokinetic Models for CYP1A2 Drug-Drug Interaction Prediction: A Modeling Network of Fluvoxamine, Theophylline, Caffeine, Rifampicin, and Midazolam.

This study provides whole-body physiologically-based pharmacokinetic models of the strong index cytochrome P450 (CYP)1A2 inhibitor and moderate CYP3A4 inhibitor fluvoxamine and of the sensitive CYP1A2 substrate theophylline. Both models were built and thoroughly evaluated for their application in drug-drug interaction (DDI) prediction in a network of perpetrator and victim drugs, combining them with previously developed models of caffeine (sensitive index CYP1A2 substrate), rifampicin (moderate CYP1A2 inducer), and midazolam (sensitive index CYP3A4 substrate). Simulation of all reported clinical DDI studies for combinations of these five drugs shows that the presented models reliably predict the observed drug concentrations, resulting in seven of eight of the predicted DDI area under the plasma curve (AUC) ratios (AUC during DDI/AUC control) and seven of seven of the predicted DDI peak plasma concentration (Cmax ) ratios (Cmax during DDI/Cmax control) within twofold of the observed values. Therefore, the models are considered qualified for DDI prediction. All models are comprehensively documented and publicly available, as tools to support the drug development and clinical research community.

PBPK Models for CYP3A4 and P-gp DDI Prediction: A Modeling Network of Rifampicin, Itraconazole, Clarithromycin, Midazolam, Alfentanil, and Digoxin.

According to current US Food and Drug Administration (FDA) and European Medicines Agency (EMA) guidance documents, physiologically based pharmacokinetic (PBPK) modeling is a powerful tool to explore and quantitatively predict drug-drug interactions (DDIs) and may offer an alternative to dedicated clinical trials. This study provides whole-body PBPK models of rifampicin, itraconazole, clarithromycin, midazolam, alfentanil, and digoxin within the Open Systems Pharmacology (OSP) Suite. All models were built independently, coupled using reported interaction parameters, and mutually evaluated to verify their predictive performance by simulating published clinical DDI studies. In total, 112 studies were used for model development and 57 studies for DDI prediction. 93% of the predicted area under the plasma concentration-time curve (AUC) ratios and 94% of the peak plasma concentration (Cmax ) ratios are within twofold of the observed values. This study lays a cornerstone for the qualification of the OSP platform with regard to reliable PBPK predictions of enzyme-mediated and transporter-mediated DDIs during model-informed drug development. All presented models are provided open-source and transparently documented.

A physiologically based pharmacokinetic (PBPK) parent-metabolite model of the chemotherapeutic zoptarelin doxorubicin-integration of in vitro results, Phase I and Phase II data and model application for drug-drug interaction potential analysis.

PURPOSE: Zoptarelin doxorubicin is a fusion molecule of the chemotherapeutic doxorubicin and a luteinizing hormone-releasing hormone receptor (LHRHR) agonist, designed for drug targeting to LHRHR positive tumors. The aim of this study was to establish a physiologically based pharmacokinetic (PBPK) parent-metabolite model of zoptarelin doxorubicin and to apply it for drug-drug interaction (DDI) potential analysis. METHODS: The PBPK model was built in a two-step procedure. First, a model for doxorubicin was developed, using clinical data of a doxorubicin study arm. Second, a parent-metabolite model for zoptarelin doxorubicin was built, using clinical data of three different zoptarelin doxorubicin studies with a dosing range of 10-267 mg/m2, integrating the established doxorubicin model. DDI parameters determined in vitro were implemented to predict the impact of zoptarelin doxorubicin on possible victim drugs. RESULTS: In vitro, zoptarelin doxorubicin inhibits the drug transporters organic anion-transporting polypeptide 1B3 (OATP1B3) and organic cation transporter 2 (OCT2). The model was applied to evaluate the in vivo inhibition of these transporters in a generic manner, predicting worst-case scenario decreases of 0.5% for OATP1B3 and of 2.5% for OCT2 transport rates. Specific DDI simulations using PBPK models of simvastatin (OATP1B3 substrate) and metformin (OCT2 substrate) predict no significant changes of the plasma concentrations of these two victim drugs during co-administration. CONCLUSIONS: The first whole-body PBPK model of zoptarelin doxorubicin and its active metabolite doxorubicin has been successfully established. Zoptarelin doxorubicin shows no potential for DDIs via OATP1B3 and OCT2.

A physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) model of the histone deacetylase (HDAC) inhibitor vorinostat for pediatric and adult patients and its application for dose specification.

PURPOSE: This study aimed at recommending pediatric dosages of the histone deacetylase (HDAC) inhibitor vorinostat and potentially more effective adult dosing regimens than the approved standard dosing regimen of 400 mg/day, using a comprehensive physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling approach. METHODS: A PBPK/PD model for vorinostat was developed for predictions in adults and children. It includes the maturation of relevant metabolizing enzymes. The PBPK model was expanded by (1) effect compartments to describe vorinostat concentration-time profiles in peripheral blood mononuclear cells (PBMCs), (2) an indirect response model to predict the HDAC inhibition, and (3) a thrombocyte model to predict the dose-limiting thrombocytopenia. Parameterization of drug and system-specific processes was based on published and unpublished in silico, in vivo, and in vitro data. The PBPK modeling software used was PK-Sim and MoBi. RESULTS: The PBPK/PD model suggests dosages of 80 and 230 mg/m2 for children of 0-1 and 1-17 years of age, respectively. In comparison with the approved standard treatment, in silico trials reveal 11 dosing regimens (9 oral, and 2 intravenous infusion rates) increasing the HDAC inhibition by an average of 31%, prolonging the HDAC inhibition by 181%, while only decreasing the circulating thrombocytes to a tolerable 53%. The most promising dosing regimen prolongs the HDAC inhibition by 509%. CONCLUSIONS: Thoroughly developed PBPK models enable dosage recommendations in pediatric patients and integrated PBPK/PD models, considering PD biomarkers (e.g., HDAC activity and platelet count), are well suited to guide future efficacy trials by identifying dosing regimens potentially superior to standard dosing regimens.

Clarithromycin, Midazolam, and Digoxin: Application of PBPK Modeling to Gain New Insights into Drug-Drug Interactions and Co-medication Regimens.

Clarithromycin is a substrate and mechanism-based inhibitor of cytochrome P450 (CYP) 3A4 as well as a substrate and competitive inhibitor of P-glycoprotein (P-gp) and organic anion-transporting polypeptides (OATP) 1B1 and 1B3. Administered concomitantly, clarithromycin causes drug-drug interactions (DDI) with the victim drugs midazolam (CYP3A4 substrate) and digoxin (P-gp substrate). The objective of the presented study was to build a physiologically based pharmacokinetic (PBPK) DDI model for clarithromycin, midazolam, and digoxin and to exemplify dosing adjustments under clarithromycin co-treatment. The PBPK model development included an extensive literature search for representative PK studies and for compound characteristics of clarithromycin, midazolam, and digoxin. Published concentration-time profiles were used for model development (training dataset), and published and unpublished individual profiles were used for model evaluation (evaluation dataset). The developed single-compound PBPK models were linked for DDI predictions. The full clarithromycin DDI model successfully predicted the metabolic (midazolam) and transporter (digoxin) DDI, the acceptance criterion (0.5 ≤ AUCratio,predicted/AUCratio,observed ≤ 2) was met by all predictions. During co-treatment with 250 or 500 mg clarithromycin (bid), the midazolam and digoxin doses should be reduced by 74 to 88% and by 21 to 22%, respectively, to ensure constant midazolam and digoxin exposures (AUC). With these models, we provide highly mechanistic tools to help researchers understand and characterize the DDI potential of new molecular entities and inform the design of DDI studies with potential CYP3A4 and P-gp substrates.