Open Systems Pharmacology Research Articles

PBPK-based translation from preclinical species to humans for the full-size IgG therapeutic efalizumab.

Animal models play a vital role in pharmaceutical research and development by supporting the planning and design of later clinical studies. To improve confidence and reliability of first in human dose estimates it is essential to assess the comparability of animal studies with the human situation. In the context of large molecules, it is particularly important to evaluate the cross-species-translatability of parameters related to neonatal fragment crystallizable receptor (FcRn) binding and target mediated drug disposition (TMDD), as they greatly influence distribution and disposition of proteins in the body of an organism. Plasma pharmacokinetic data of the therapeutic protein efalizumab were obtained from literature. Physiologically based pharmacokinetic (PBPK) models were built for three different species (rabbit, non-human primate (NHP), human). Target binding was included in the NHP and human models. The assumption of similar target turnover and target-binding in NHP and human was explored, to gain insights into how these parameters might be translated between species. Efalizumab PBPK models were successfully developed for three species and concentration-time-profiles could be described appropriately across different intravenously administered doses. The final NHP and human models feature a common set of parameters for target turnover and drug-target-complex internalization, as well as comparable target-binding parameters. Our analyses show that different parameter values for FcRn affinity are crucial to accurately describe the concentration-time profiles. Based on the available data in rabbits, NHP and humans, parameters for FcRn affinity cannot be translated between species, but parameters related to target mediated drug disposition can be translated from NHP to human. The inclusion of additional pharmacokinetic (PK) data including different efalizumab doses would further support and confirm our findings on identifying TMDD and, thus, binding kinetics of efalizumab in NHPs. Furthermore, we suggest that information on target expression and internalization rates could make it possible to develop comprehensive human PBPK models with minimal animal testing. In this project, we compared the pharmacokinetics of a therapeutic protein in rabbit, NHP and human using an open PBPK modeling platform (Open Systems Pharmacology Suite, http://www.open-systems-pharmacology.org). Our findings could support similar translatory studies for first in human dose predictions in the future.

In-Depth Analysis of the Selection of PBPK Modeling Tools: Bibliometric and Social Network Analysis of the Open Systems Pharmacology Community.

Since the Open Source Initiative laid the foundation for the open source software environment in 1998, the popularity of free and open source software has been steadily increasing. Model-informed drug discovery and development (MID3), a key component of pharmaceutical research and development, heavily makes use of computational models which can be developed using various software including the Open Systems Pharmacology (OSP) software (PK-Sim/MoBi), a free and open source software tool for physiologically based pharmacokinetic (PBPK) modeling. In this study, we aimed to investigate the impact, application areas, and reach of the OSP software as well as the relationships and collaboration patterns between organizations having published OSP-related articles between 2017 and 2023. Therefore, we conducted a bibliometric analysis of OSP-related publications and a social network analysis of the organizations with which authors of OSP-related publications were affiliated. On several levels, we found evidence for a significant growth in the size of the OSP community as well as its visibility in the MID3 community since OSP's establishment in 2017. Specifically, the annual publication rate of PubMed-indexed PBPK-related articles using the OSP software outpaced that of PBPK-related articles using any software. Our bibliometric analysis and network analysis demonstrated that the expansion of the OSP community was predominantly driven by new authors and organizations without prior connections to the community involving the generation of research clusters de novo and an overall diversification of the network. These findings suggest an ongoing evolution of the OSP community toward a more segmented, diverse, and inclusive network.

P05-23 Learning by doing: modelling and simulation for risk assessment in the open-source Open Systems Pharmacology Suite (PK-Sim® and MoBi®)

P05-23 Learning by doing: modelling and simulation for risk assessment in the open-source Open Systems Pharmacology Suite (PK-Sim® and MoBi®)

Physiologically Based Pharmacokinetic Modeling of Nilotinib for Drug-Drug Interactions, Pediatric Patients, and Pregnancy and Lactation.

Nilotinib is a second-generation BCR-ABL tyrosine kinase inhibitor for the treatment of Philadelphia chromosome-positive chronic myeloid leukemia in both adult and pediatric patients. The pharmacokinetics (PK) of nilotinib in specific populations such as pregnant and lactating people remain poorly understood. Therefore, the objectives of the current study were to develop a physiologically based pharmacokinetic (PBPK) model to predict nilotinib PK in virtual drug-drug interaction (DDI) studies, as well as in pediatric, pregnant, and lactating populations. The nilotinib PBPK model was built in PK-Sim, which is part of the free and open-source software Open Systems Pharmacology. The observed clinical data for the validation of the nilotinib models were obtained from the literature. The model reasonably predicted nilotinib concentrations in the adult population; the DDIs between nilotinib and rifampin or ketoconazole in the adult population; and the PK in the pediatric, pregnant, and lactating populations, although in the latter 2 populations plasma concentrations were slightly underestimated. The ratio of predicted versus observed PK parameters for the adult model ranged from 0.71 to 1.11 for area under the concentration-time curve and 0.55 to 0.95 for maximum concentration. For the DDI, the predicted area under the concentration-time curve ratio and maximum concentration ratio fell within the Guest criterion. The current study demonstrated the utility of using PBPK modeling to understand the mechanistic basis of PK differences between adults and specific populations, such as pediatrics, and pregnant and lactating individuals, indicating that this technology can potentially inform or optimize dosing conditions in specific populations.

Development of a Generic Fetal Physiologically Based Pharmacokinetic Model and Prediction of Human Maternal and Fetal Organ Concentrations of Cefuroxime.

Physiologically based pharmacokinetic (PBPK) models for pregnant women have recently been successfully used to predict maternal and umbilical cord pharmacokinetics (PK). Because there is very limited opportunity for conducting clinical and PK investigations for fetal drug exposure, PBPK models may provide further insights. The objectives of this study were to extend a whole-body pregnancy PBPK model by multiple compartments representing fetal organs, and to predict the PK of cefuroxime in the maternal and fetal plasma, the amniotic fluid, and several fetal organs. To this end, a previously developed pregnancy PBPK model for cefuroxime was updated using the open-source software Open Systems Pharmacology (PK-Sim®/MoBi®). Multiple compartments were implemented to represent fetal organs including brain, heart, liver, lungs, kidneys, the gastrointestinal tract (GI), muscles, and fat tissue, as well as another compartment lumping organs and tissues not explicitly represented. This novel PBPK model successfully predicted cefuroxime concentrations in maternal blood, umbilical cord, amniotic fluid, and several fetal organs including heart, liver, and lungs. Further model validation with additional clinical PK data is needed to build confidence in the model. Being developed with an open-source software, the presented generic model can be freely re-used and tailored to address specific questions at hand, e.g., to assist the design of clinical studies in the context of drug research or to predict fetal organ concentrations of chemicals in the context of fetal health risk assessment.

Translation of Monoclonal Antibodies Pharmacokinetics from Animal to Human Using Physiologically Based Modeling in Open Systems Pharmacology (OSP) Suite: A Retrospective Analysis of Bevacizumab.

Interspecies translation of monoclonal antibodies (mAbs) pharmacokinetics (PK) in presence of target-mediated drug disposition (TMDD) is particularly challenging. Incorporation of TMDD in physiologically based PK (PBPK) modeling is recent and needs to be consolidated and generalized to provide better prediction of TMDD regarding inter-species translation during preclinical and clinical development steps of mAbs. The objective of this study was to develop a generic PBPK translational approach for mAbs using the open-source software (PK-Sim® and Mobi®). The translation of bevacizumab based on data in non-human primates (NHP), healthy volunteers (HV), and cancer patients was used as a case example for model demonstration purpose. A PBPK model for bevacizumab concentration-time data was developed using data from literature and the Open Systems Pharmacology (OSP) Suite version 10. PK-sim® was used to build the linear part of bevacizumab PK (mainly FcRn-mediated), whereas MoBi® was used to develop the target-mediated part. The model was first developed for NHP and used for a priori PK prediction in HV. Then, the refined model obtained in HV was used for a priori prediction in cancer patients. A priori predictions were within 2-fold prediction error (predicted/observed) for both area under the concentration-time curve (AUC) and maximum concentration (Cmax) and all the predicted concentrations were within 2-fold average fold error (AFE) and average absolute fold error (AAFE). Sensitivity analysis showed that FcRn-mediated distribution and elimination processes must be accounted for at all mAb concentration levels, whereas the lower the mAb concentration, the more significant the target-mediated elimination. This project is the first step to generalize the full PBPK translational approach in Model-Informed Drug Development (MIDD) of mAbs using OSP Suite.

Generic Workflow to Predict Medicine Concentrations in Human Milk Using Physiologically-Based Pharmacokinetic (PBPK) Modelling-A Contribution from the ConcePTION Project.

Women commonly take medication during lactation. Currently, there is little information about the exposure-related safety of maternal medicines for breastfed infants. The aim was to explore the performance of a generic physiologically-based pharmacokinetic (PBPK) model to predict concentrations in human milk for ten physiochemically diverse medicines. First, PBPK models were developed for "non-lactating" adult individuals in PK-Sim/MoBi v9.1 (Open Systems Pharmacology). The PBPK models predicted the area-under-the-curve (AUC) and maximum concentrations (Cmax) in plasma within a two-fold error. Next, the PBPK models were extended to include lactation physiology. Plasma and human milk concentrations were simulated for a three-months postpartum population, and the corresponding AUC-based milk-to-plasma (M/P) ratios and relative infant doses were calculated. The lactation PBPK models resulted in reasonable predictions for eight medicines, while an overprediction of human milk concentrations and M/P ratios (>2-fold) was observed for two medicines. From a safety perspective, none of the models resulted in underpredictions of observed human milk concentrations. The present effort resulted in a generic workflow to predict medicine concentrations in human milk. This generic PBPK model represents an important step towards an evidence-based safety assessment of maternal medication during lactation, applicable in an early drug development stage.

Gene expression databases for physiologically based pharmacokinetic modeling of humans and animal species.

In drug research, developing a sound understanding of the key mechanistic drivers of pharmacokinetics (PK) for new molecular entities is essential for human PK and dose predictions. Here, characterizing the absorption, distribution, metabolism, and excretion (ADME) processes is crucial for a mechanistic understanding of the drug-target and drug-body interactions. Sufficient knowledge on ADME processes enables reliable interspecies and human PK estimations beyond allometric scaling. The physiologically based PK (PBPK) modeling framework allows the explicit consideration of organ-specific ADME processes. The sum of all passive and active ADME processes results in the observed plasma PK. Gene expression information can be used as surrogate for protein abundance and activity within PBPK models. The absolute and relative expression of ADME genes can differ between species and strains. This is affecting both, the PK and pharmacodynamics and is therefore posing a challenge for the extrapolation from preclinical findings to humans. We developed an automated workflow that generates whole-body gene expression databases for humans and other species relevant in drug development, animal health, nutritional sciences, and toxicology. Solely, bulk RNA-seq data curated and provided by the Swiss Institute of Bioinformatics from healthy, normal, and untreated primary tissue samples were considered as an unbiased reference of normal gene expression. The databases are interoperable with the Open Systems Pharmacology Suite (PK-Sim and MoBi) and enable seamless access to a central source of curated cross-species gene expression data. This will increase data transparency, increase reliability and reproducibility of PBPK model simulations, and accelerate mechanistic PBPK model development in the future.

A generic avian physiologically-based kinetic (PBK) model and its application in three bird species

Physiologically-based kinetic (PBK) models are effective tools for designing toxicological studies and conducting extrapolations to inform hazard characterization in risk assessment by filling data gaps and defining safe levels of chemicals. In the present work, a generic avian PBK model for male and female birds was developed using PK-Sim and MoBi from the Open Systems Pharmacology Suite (OSPS). The PBK model includes an ovulation model (egg development) to predict concentrations of chemicals in eggs from dietary exposure. The model was parametrized for chicken (Gallus gallus), bobwhite quail (Colinus virginianus) and mallard duck (Anas platyrhynchos) and was tested with nine chemicals for which in vivo studies were available. Time-concentration profiles of chemicals reaching tissues and egg compartment were simulated and compared to in vivo data. The overall accuracy of the PBK model predictions across the analyzed chemicals was good. Model simulations were found to be in the range of 22–79% within a 3-fold and 41–89% were within 10- fold deviation of the in vivo observed data. However, for some compounds scarcity of in-vivo data and inconsistencies between published studies allowed only a limited goodness of fit evaluation. The generic avian PBK model was developed following a “best practice” workflow describing how to build a PBK model for novel species. The credibility and reproducibility of the avian PBK models were scored by evaluation according to the available guidance documents from WHO (2010), and OECD (2021), to increase applicability, confidence and acceptance of these in silico models in chemical risk assessment.

A Physiologically Based Pharmacokinetic and Pharmacodynamic Model of the CYP3A4 Substrate Felodipine for Drug-Drug Interaction Modeling.

The antihypertensive felodipine is a calcium channel blocker of the dihydropyridine type, and its pharmacodynamic effect directly correlates with its plasma concentration. As a sensitive substrate of cytochrome P450 (CYP) 3A4 with high first-pass metabolism, felodipine shows low oral bioavailability and is susceptible to drug–drug interactions (DDIs) with CYP3A4 perpetrators. This study aimed to develop a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) parent–metabolite model of felodipine and its metabolite dehydrofelodipine for DDI predictions. The model was developed in PK-Sim® and MoBi® using 49 clinical studies (94 plasma concentration–time profiles in total) that investigated different doses (1–40 mg) of the intravenous and oral administration of felodipine. The final model describes the metabolism of felodipine to dehydrofelodipine by CYP3A4, sufficiently capturing the first-pass metabolism and the subsequent metabolism of dehydrofelodipine by CYP3A4. Diastolic blood pressure and heart rate PD models were included, using an Emax function to describe the felodipine concentration–effect relationship. The model was tested in DDI predictions with itraconazole, erythromycin, carbamazepine, and phenytoin as CYP3A4 perpetrators, with all predicted DDI AUClast and Cmax ratios within two-fold of the observed values. The model will be freely available in the Open Systems Pharmacology model repository and can be applied in DDI predictions as a CYP3A4 victim drug.

Physiologically Based Pharmacokinetic (PBPK) Modelling to Predict Breast Milk Exposure of Valproic Acid

Breastfeeding has many positive effects on the health and well‐being of mother and baby. During lactation, a pharmacological therapy may be required. Currently there is a knowledge gap regarding the safety of medicines during lactation, which may lead to unnecessary cessation of breastfeeding or incur unknown risks to the child when continuing breastfeeding while on medication. The Innovative Medicines Initiative (IMI) ConcePTION aims to address this and other issues related to medication use during pregnancy and lactation.The objective of this study performed in the IMI ConcePTION framework was to explore the utility of Physiologically‐Based Pharmacokinetic (PBPK) modelling in the generation of information regarding the safety of medicines (e.g. valproic acid, VPA) during breastfeeding.To achieve this goal, PK‐Sim® (v 9.1) (Open Systems Pharmacology) was used to first develop a reference PBPK model for VPA in adults using 7 training datasets. The median fold errors of the Area Under the Curve (AUC) and maximum plasma concentration (Cmax) for the 17 validation datasets were 0.89 (range: 0.74‐0.134) and 0.90 (range: 0.60‐1.36), indicating a good predictive performance. In a next step, this PBPK model was extended with a breast compartment (including a subcompartment representing the breast milk) according to the work from Dallmann et al. (2018). Transfer of medicines between the milk and the intracellular compartment was parametrized with a clearance for secretion to milk (Clsec), and a clearance for reuptake (Clre). These clearance values were determined based on physicochemical properties using the empirical model from Koshimichi et al. (2011). Clsec and Clre were 235.05 and 1452.78 mL/h. The PBPK model for lactating women was used to predict exposure of VPA in plasma and human breast milk. The simulated milk‐to‐plasma (M/P) ratio was 0.016, similar to reported M/P ratio’s (range: 0.007‐0.25). Finally, the daily infant dose was calculated based on the simulated concentration in the milk and used as input for an infant PBPK model. The relative infant dose was 1.4 %. Average steady‐state infant plasma concentration was estimated to be 0.83 mg/L (relative infant exposure: 1.2 %), compared to a wide range of observations starting from undetectable concentrations up to 41 mg/L.PBPK modelling was successfully used to predict the transfer of VPA into the human breast milk, and to estimate subsequent infant exposure via breastfeeding. Based on ongoing application to 11 other medicines, a generic workflow for PBPK‐based breast milk transfer of medicines has been developed. Current efforts to determine the in vitro permeability of VPA across the blood‐milk barrier may obviate the current need for (sparse) milk concentrations. This study illustrated the potential of PBPK modelling as a tool to fill the knowledge gap regarding the exposure to medicines in lactating women and their infants.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Clopidogrel and Its Four Relevant Metabolites for CYP2B6, CYP2C8, CYP2C19, and CYP3A4 Drug-Drug-Gene Interaction Predictions.

The antiplatelet agent clopidogrel is listed by the FDA as a strong clinical index inhibitor of cytochrome P450 (CYP) 2C8 and weak clinical inhibitor of CYP2B6. Moreover, clopidogrel is a substrate of—among others—CYP2C19 and CYP3A4. This work presents the development of a whole-body physiologically based pharmacokinetic (PBPK) model of clopidogrel including the relevant metabolites, clopidogrel carboxylic acid, clopidogrel acyl glucuronide, 2-oxo-clopidogrel, and the active thiol metabolite, with subsequent application for drug–gene interaction (DGI) and drug–drug interaction (DDI) predictions. Model building was performed in PK-Sim® using 66 plasma concentration-time profiles of clopidogrel and its metabolites. The comprehensive parent-metabolite model covers biotransformation via carboxylesterase (CES) 1, CES2, CYP2C19, CYP3A4, and uridine 5′-diphospho-glucuronosyltransferase 2B7. Moreover, CYP2C19 was incorporated for normal, intermediate, and poor metabolizer phenotypes. Good predictive performance of the model was demonstrated for the DGI involving CYP2C19, with 17/19 predicted DGI AUClast and 19/19 predicted DGI Cmax ratios within 2-fold of their observed values. Furthermore, DDIs involving bupropion, omeprazole, montelukast, pioglitazone, repaglinide, and rifampicin showed 13/13 predicted DDI AUClast and 13/13 predicted DDI Cmax ratios within 2-fold of their observed ratios. After publication, the model will be made publicly accessible in the Open Systems Pharmacology repository.

Physiologically Based Pharmacokinetic Modeling of Rosuvastatin to Predict Transporter-Mediated Drug-Drug Interactions

PurposeTo 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).MethodsThe 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.ResultsThe 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 AUClast ratios (AUClast during DDI/AUClast without co-administration) and DDI Cmax ratios (Cmax during DDI/Cmax without co-administration), with all simulated DDI ratios within 1.6-fold of the observed values.ConclusionsA 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).

A generic framework for the physiologically‐based pharmacokinetic platform qualification of PK‐Sim and its application to predicting cytochrome P450 3A4–mediated drug–drug interactions

The success of applications of physiologically‐based pharmacokinetic (PBPK) modeling in drug development and drug labeling has triggered regulatory agencies to demand rigorous demonstration of the predictive capability of the specific PBPK platform for a particular intended application purpose. The effort needed to comply with such qualification requirements exceeds the costs for any individual PBPK application. Because changes or updates of a PBPK platform would require (re‐)qualification, a reliable and efficient generic qualification framework is needed. We describe the development and implementation of an agile and sustainable technical framework for automatic PBPK platform (re‐)qualification of PK‐Sim® embedded in the open source and open science GitHub landscape of Open Systems Pharmacology. The qualification approach enables the efficient assessment of all aspects relevant to the qualification of a particular purpose and provides transparency and traceability for all stakeholders. As a showcase example for the power and versatility of the qualification framework, we present the qualification of PK‐Sim® for the intended purpose of predicting cytochrome P450 3A4 (CYP3A4)–mediated drug–drug interactions (DDIs). Several perpetrator PBPK models featuring various degrees of CYP3A4 modulation and different types of mechanisms (competitive inhibition, mechanism‐based inactivation, and induction) were coupled with a set of PBPK models of sensitive CYP3A4 victim drugs. Simulations were compared to a comprehensive data set of 135 observations from published clinical DDI studies. The platform's overall predictive performance showed reasonable accuracy and precision (geometric mean fold error of 1.4 for both area under the plasma concentration‐time curve ratios and peak plasma concentration ratios with/without perpetrator) and suggests that PK‐Sim® can be applied to quantitatively assess CYP3A4‐mediated DDI in clinically untested scenarios.

Physiologically Based Pharmacokinetic Modeling of Bupropion and Its Metabolites in a CYP2B6 Drug-Drug-Gene Interaction Network.

The noradrenaline and dopamine reuptake inhibitor bupropion is metabolized by CYP2B6 and recommended by the FDA as the only sensitive substrate for clinical CYP2B6 drug–drug interaction (DDI) studies. The aim of this study was to build a whole-body physiologically based pharmacokinetic (PBPK) model of bupropion including its DDI-relevant metabolites, and to qualify the model using clinical drug–gene interaction (DGI) and DDI data. The model was built in PK-Sim® applying clinical data of 67 studies. It incorporates CYP2B6-mediated hydroxylation of bupropion, metabolism via CYP2C19 and 11β-HSD, as well as binding to pharmacological targets. The impact of CYP2B6 polymorphisms is described for normal, poor, intermediate, and rapid metabolizers, with various allele combinations of the genetic variants CYP2B6*1, *4, *5 and *6. DDI model performance was evaluated by prediction of clinical studies with rifampicin (CYP2B6 and CYP2C19 inducer), fluvoxamine (CYP2C19 inhibitor) and voriconazole (CYP2B6 and CYP2C19 inhibitor). Model performance quantification showed 20/20 DGI ratios of hydroxybupropion to bupropion AUC ratios (DGI AUCHBup/Bup ratios), 12/13 DDI AUCHBup/Bup ratios, and 7/7 DDGI AUCHBup/Bup ratios within 2-fold of observed values. The developed model is freely available in the Open Systems Pharmacology model repository.

Pharmacokinetics of the CYP3A4 and CYP2B6 Inducer Carbamazepine and Its Drug-Drug Interaction Potential: A Physiologically Based Pharmacokinetic Modeling Approach.

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 AUClast ratios and 11/15 DDI Cmax 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.

Physiologically Based Pharmacokinetic Modeling of Metoprolol Enantiomers and α-Hydroxymetoprolol to Describe CYP2D6 Drug-Gene Interactions.

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 AUClast (1.27) and Cmax values (1.23) over all studies. For DGI predictions, 18 out of 18 DGI AUClast ratios and 18 out of 18 DGI Cmax 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.

A Physiologically-Based Pharmacokinetic Model of Trimethoprim for MATE1, OCT1, OCT2, and CYP2C8 Drug-Drug-Gene Interaction Predictions.

Trimethoprim is a frequently-prescribed antibiotic and therefore likely to be co-administered with other medications, but it is also a potent inhibitor of multidrug and toxin extrusion protein (MATE) and a weak inhibitor of cytochrome P450 (CYP) 2C8. The aim of this work was to develop a physiologically-based pharmacokinetic (PBPK) model of trimethoprim to investigate and predict its drug–drug interactions (DDIs). The model was developed in PK-Sim®, using a large number of clinical studies (66 plasma concentration–time profiles with 36 corresponding fractions excreted in urine) to describe the trimethoprim pharmacokinetics over the entire published dosing range (40 to 960 mg). The key features of the model include intestinal efflux via P-glycoprotein (P-gp), metabolism by CYP3A4, an unspecific hepatic clearance process, and a renal clearance consisting of glomerular filtration and tubular secretion. The DDI performance of this new model was demonstrated by prediction of DDIs and drug–drug–gene interactions (DDGIs) of trimethoprim with metformin, repaglinide, pioglitazone, and rifampicin, with all predicted DDI and DDGI AUClast and Cmax ratios within 1.5-fold of the clinically-observed values. The model will be freely available in the Open Systems Pharmacology model repository, to support DDI studies during drug development.

Physiologically-based pharmacokinetic (PBPK) exploration of extrinsic factors influencing vericiguat pharmacokinetics

Abstract Introduction Vericiguat is a once daily, novel oral stimulator of soluble guanylate cyclase (sGC) that showed clinical benefit in the Phase III VICTORIA study in heart failure patients with reduced ejection fraction (HFrEF, NCT02861534). Nonclinical and clinical studies demonstrated that the primary route of elimination of vericiguat was glucuronidation to an inactive metabolite M-1 (N-glucuronide). This glucuronidation was catalyzed by uridine 5'-diphospho-glucuronosyltransferases (UGT)1A9 as well as UGT1A1, thus vericiguat may have a potential for victim drug-drug interaction (DDI) when co-administered with potent UGT inhibitors. Purpose In a clinical DDI study with mefenamic acid as an UGT1A9 inhibitor no clinically relevant increase in vericiguat exposure in healthy subjects was observed (EudraCT 2014–000764–17). This analysis aims to prospectively investigate as extrinsic factors the DDI potential with atazanavir as a selective UGT1A1 inhibitor via full dynamic physiologically-based pharmacokinetic (PBPK) modelling. Methods A PBPK model for vericiguat and M-1 in healthy adults was built with PK-Sim (PBPK platform as part of the Open Systems Pharmacology Suite) by integrating physicochemical, in vitro metabolism and transporter data as well as PK data from clinical pharmacology studies in order to assess the victim DDI potential of vericiguat when co-administered with UGT inhibitors. First, PBPK models for mefenamic acid and atazanavir were separately developed and verified using published literature data. The PBPK model for vericiguat was then verified with regard to its fraction of metabolism by UGTs by comparing simulated and observed data of the clinical mefenamic acid DDI study. Finally, the UGT1A1 DDI potential of vericiguat was prospectively predicted by simulating an in silico study between the UGT1A1 inhibitor atazanavir and vericiguat. Results In line with the results of the clinical DDI study with mefenamic acid, an increase in total vericiguat exposure by 14% (area under the concentration time curve ratio (AUCR) of 1.14 (geoCV 5.3%; 90% population interval: 1.06 to 1.25) and peak exposure increase by 6% (CmaxR of 1.06; geoCV 5.9%; 90% population interval: 1.01 to 1.20) was simulated using the PBPK model. A prospective prediction of a virtual DDI trial between the UGT1A1 inhibitor atazanavir yielded an AUCR of 1.12 (geoCV 2.9%; 90% population interval: 1.07 to 1.17) and a CmaxR of 1.04 (geoCV 1.1%; 90% population interval: 1.03 to 1.06). The proposed population intervals for AUCR and CmaxR for both DDI studies lie within the default no-effect boundary of 0.80 to 1.25 according to the to January 2020 FDA DDI guideline. Conclusion(s) Results of UGT1A9-DDI simulations were consistent with those of the clinical study-The prospective UGT1A1-DDI simulation results suggest a low potential for vericiguat to be subject to DDI when co-administered with UGT1A1 inhibitors. Funding Acknowledgement Type of funding source: Private company. Main funding source(s): Funding for this research was provided by Bayer and Merck Sharp & Dohme Corp., a subsidiary of Merck & Co., Inc., Kenilworth, NJ, USA

A Mechanistic, Enantioselective, Physiologically Based Pharmacokinetic Model of Verapamil and Norverapamil, Built and Evaluated for Drug-Drug Interaction Studies.

The calcium channel blocker and antiarrhythmic agent verapamil is recommended by the FDA for drug–drug interaction (DDI) studies as a moderate clinical CYP3A4 index inhibitor and as a clinical Pgp inhibitor. The purpose of the presented work was to develop a mechanistic whole-body physiologically based pharmacokinetic (PBPK) model to investigate and predict DDIs with verapamil. The model was established in PK-Sim®, using 45 clinical studies (dosing range 0.1–250 mg), including literature as well as unpublished Boehringer Ingelheim data. The verapamil R- and S-enantiomers and their main metabolites R- and S-norverapamil are represented in the model. The processes implemented to describe the pharmacokinetics of verapamil and norverapamil include enantioselective plasma protein binding, enantioselective metabolism by CYP3A4, non-stereospecific Pgp transport, and passive glomerular filtration. To describe the auto-inhibitory and DDI potential, mechanism-based inactivation of CYP3A4 and non-competitive inhibition of Pgp by the verapamil and norverapamil enantiomers were incorporated based on in vitro literature. The resulting DDI performance was demonstrated by prediction of DDIs with midazolam, digoxin, rifampicin, and cimetidine, with 21/22 predicted DDI AUC ratios or Ctrough ratios within 1.5-fold of the observed values. The thoroughly built and qualified model will be freely available in the Open Systems Pharmacology model repository to support model-informed drug discovery and development.