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 (C_max) 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.

The anticonvulsant carbamazepine is frequently used in the long-term therapy of epilepsy and is a known substrate and inducer of cytochrome P450 (CYP) 3A4 and CYP2B6. Carbamazepine induces the metabolism of various drugs (including its own); on the other hand, its metabolism can be affected by various CYP inhibitors and inducers. The aim of this work was to develop a physiologically based pharmacokinetic (PBPK) parent−metabolite model of carbamazepine and its metabolite carbamazepine-10,11-epoxide, including carbamazepine autoinduction, to be applied for drug–drug interaction (DDI) prediction. The model was developed in PK-Sim, using a total of 92 plasma concentration−time profiles (dosing range 50–800 mg), as well as fractions excreted unchanged in urine measurements. The carbamazepine model applies metabolism by CYP3A4 and CYP2C8 to produce carbamazepine-10,11-epoxide, metabolism by CYP2B6 and UDP-glucuronosyltransferase (UGT) 2B7 and glomerular filtration. The carbamazepine-10,11-epoxide model applies metabolism by epoxide hydroxylase 1 (EPHX1) and glomerular filtration. Good DDI performance was demonstrated by the prediction of carbamazepine DDIs with alprazolam, bupropion, erythromycin, efavirenz and simvastatin, where 14/15 DDI AUC_last ratios and 11/15 DDI C_max ratios were within the prediction success limits proposed by Guest et al. The thoroughly evaluated model will be freely available in the Open Systems Pharmacology model repository.

Background Metformin is a widely prescribed antidiabetic BCS Class III drug (low permeability) that depends on active transport for its absorption and disposition. It is recommended by the US Food and Drug Administration as a clinical substrate of organic cation transporter 2/multidrug and toxin extrusion protein for drug–drug interaction studies. Cimetidine is a potent organic cation transporter 2/multidrug and toxin extrusion protein inhibitor. Objective The objective of this study was to provide mechanistic whole-body physiologically based pharmacokinetic models of metformin and cimetidine, built and evaluated to describe the metformin-SLC22A2 808G>T drug–gene interaction, the cimetidine-metformin drug–drug interaction, and the impact of renal impairment on metformin exposure. Methods Physiologically based pharmacokinetic models were developed in PK-Sim^® (version 8.0). Thirty-nine clinical studies (dosing range 0.001–2550 mg), providing metformin plasma and urine data, positron emission tomography measurements of tissue concentrations, studies in organic cation transporter 2 polymorphic volunteers, drug–drug interaction studies with cimetidine, and data from patients in different stages of chronic kidney disease, were used to develop the metformin model. Twenty-seven clinical studies (dosing range 100–800 mg), reporting cimetidine plasma and urine concentrations, were used for the cimetidine model development. Results The established physiologically based pharmacokinetic models adequately describe the available clinical data, including the investigated drug–gene interaction, drug–drug interaction, and drug–drug–gene interaction studies, as well as the metformin exposure during renal impairment. All modeled drug–drug interaction area under the curve and maximum concentration ratios are within 1.5-fold of the observed ratios. The clinical data of renally impaired patients shows the expected increase in metformin exposure with declining kidney function, but also indicates counter-regulatory mechanisms in severe renal disease; these mechanisms were implemented into the model based on findings in preclinical species. Conclusions Whole-body physiologically based pharmacokinetic models of metformin and cimetidine were built and qualified for the prediction of metformin pharmacokinetics during drug–gene interaction, drug–drug interaction, and different stages of renal disease. The model files will be freely available in the Open Systems Pharmacology model repository. Current guidelines for metformin treatment of renally impaired patients should be reviewed to avoid overdosing in CKD3 and to allow metformin therapy of CKD4 patients.

Purpose: To build a physiologically based pharmacokinetic (PBPK) model of the clinical OATP1B1/OATP1B3/BCRP victim drug rosuvastatin for the investigation and prediction of its transporter-mediated drug-drug interactions (DDIs). Methods: The Rosuvastatin model was developed using the open-source PBPK software PK-Sim®, following a middle-out approach. 42 clinical studies (dosing range 0.002–80.0 mg), providing rosuvastatin plasma, urine and feces data, positron emission tomography (PET) measurements of tissue concentrations and 7 different rosuvastatin DDI studies with rifampicin, gemfibrozil and probenecid as the perpetrator drugs, were included to build and qualify the model. Results: The carefully developed and thoroughly evaluated model adequately describes the analyzed clinical data, including blood, liver, feces and urine measurements. The processes implemented to describe the rosuvastatin pharmacokinetics and DDIs are active uptake by OATP2B1, OATP1B1/OATP1B3 and OAT3, active efflux by BCRP and Pgp, metabolism by CYP2C9 and passive glomerular filtration. The available clinical rifampicin, gemfibrozil and probenecid DDI studies were modeled using in vitro inhibition constants without adjustments. The good prediction of DDIs was demonstrated by simulated rosuvastatin plasma profiles, DDI AUC_last ratios (AUC_last during DDI/AUC_last without co-administration) and DDI C_max ratios (C_max during DDI/C_max without co-administration), with all simulated DDI ratios within 1.6-fold of the observed values. Conclusions: A whole-body PBPK model of rosuvastatin was built and qualified for the prediction of rosuvastatin pharmacokinetics and transporter-mediated DDIs. The model is freely available in the Open Systems Pharmacology model repository, to support future investigations of rosuvastatin pharmacokinetics, rosuvastatin therapy and DDI studies during model-informed drug discovery and development (MID3).

The beta-blocker metoprolol (the sixth most commonly prescribed drug in the USA in 2017) is subject to considerable drug–gene interaction (DGI) effects caused by genetic variations of the CYP2D6 gene. CYP2D6 poor metabolizers (5.7% of US population) show approximately five-fold higher metoprolol exposure compared to CYP2D6 normal metabolizers. This study aimed to develop a whole-body physiologically based pharmacokinetic (PBPK) model to predict CYP2D6 DGIs with metoprolol. The metoprolol (R)- and (S)-enantiomers as well as the active metabolite α-hydroxymetoprolol were implemented as model compounds, employing data of 48 different clinical studies (dosing range 5–200 mg). To mechanistically describe the effect of CYP2D6 polymorphisms, two separate metabolic CYP2D6 pathways (α-hydroxylation and O-demethylation) were incorporated for both metoprolol enantiomers. The good model performance is demonstrated in predicted plasma concentration–time profiles compared to observed data, goodness-of-fit plots, and low geometric mean fold errors of the predicted AUC_last (1.27) and C_max values (1.23) over all studies. For DGI predictions, 18 out of 18 DGI AUC_last ratios and 18 out of 18 DGI C_max ratios were within two-fold of the observed ratios. The newly developed and carefully validated model was applied to calculate dose recommendations for CYP2D6 polymorphic patients and will be freely available in the Open Systems Pharmacology repository.