A physiologically-based pharmacokinetic precision dosing approach to manage dasatinib drug-drug interactions.

Dasatinib, a second-generation tyrosine kinase inhibitor, is approved for treating chronic myeloid and acute lymphoblastic leukemia. As a sensitive cytochrome P450 (CYP) 3A4 substrate and weak base with strong pH-sensitive solubility, dasatinib is susceptible to enzyme-mediated drug-drug interactions (DDIs) with CYP3A4 perpetrators and pH-dependent DDIs with acid-reducing agents. This work aimed to develop a whole-body physiologically-based pharmacokinetic (PBPK) model of dasatinib to describe and predict enzyme-mediated and pH-dependent DDIs, to evaluate the impact of strong and moderate CYP3A4 inhibitors and inducers on dasatinib exposure and to support optimized dasatinib dosing. Overall, 63 plasma profiles from perorally administered dasatinib in healthy volunteers and cancer patients were used for model development. The model accurately described and predicted plasma profiles with geometric mean fold errors (GMFEs) for area under the concentration-time curve from the first to the last timepoint of measurement (AUClast) and maximum plasma concentration (Cmax) of 1.27 and 1.29, respectively. Regarding the DDI studies used for model development, all (8/8) predicted AUClast and Cmax ratios were within twofold of observed ratios. Application of the PBPK model for dose adaptations within various DDIs revealed dasatinib dose reductions of 50%-80% for strong and 0%-70% for moderate CYP3A4 inhibitors and a 2.3-3.1-fold increase of the daily dasatinib dose for CYP3A4 inducers to match the exposure of dasatinib administered alone. The developed model can be further employed to personalize dasatinib therapy, thereby help coping with clinical challenges resulting from DDIs and patient-related factors, such as elevated gastric pH.

Physiologically based pharmacokinetic modeling of imatinib and N-desmethyl imatinib for drug-drug interaction predictions.

The first-generation tyrosine kinase inhibitor imatinib has revolutionized the development of targeted cancer therapy and remains among the frontline treatments, for example, against chronic myeloid leukemia. As a substrate of cytochrome P450 (CYP) 2C8, CYP3A4, and various transporters, imatinib is highly susceptible to drug-drug interactions (DDIs) when co-administered with corresponding perpetrator drugs. Additionally, imatinib and its main metabolite N-desmethyl imatinib (NDMI) act as inhibitors of CYP2C8, CYP2D6, and CYP3A4 affecting their own metabolism as well as the exposure of co-medications. This work presents the development of a parent-metabolite whole-body physiologically based pharmacokinetic (PBPK) model for imatinib and NDMI used for the investigation and prediction of different DDI scenarios centered around imatinib as both a victim and perpetrator drug. Model development was performed in PK-Sim® using a total of 60 plasma concentration-time profiles of imatinib and NDMI in healthy subjects and cancer patients. Metabolism of both compounds was integrated via CYP2C8 and CYP3A4, with imatinib additionally transported via P-glycoprotein. The subsequently developed DDI network demonstrated good predictive performance. DDIs involving imatinib and NDMI were simulated with perpetrator drugs rifampicin, ketoconazole, and gemfibrozil as well as victim drugs simvastatin and metoprolol. Overall, 12/12 predicted DDI area under the curve determined between first and last plasma concentration measurements (AUClast) ratios and 12/12 predicted DDI maximum plasma concentration (Cmax) ratios were within twofold of the respective observed ratios. Potential applications of the final model include model-informed drug development or the support of model-informed precision dosing.

Physiologically-based pharmacokinetic modeling of quinidine to establish a CYP3A4, P-gp, and CYP2D6 drug-drug-gene interaction network.

The antiarrhythmic agent quinidine is a potent inhibitor of cytochrome P450 (CYP) 2D6 and P-glycoprotein (P-gp) and is therefore recommended for use in clinical drug-drug interaction (DDI) studies. However, as quinidine is also a substrate of CYP3A4 and P-gp, it is susceptible to DDIs involving these proteins. Physiologically-based pharmacokinetic (PBPK) modeling can help to mechanistically assess the absorption, distribution, metabolism, and excretion processes of a drug and has proven its usefulness in predicting even complex interaction scenarios. The objectives of the presented work were to develop a PBPK model of quinidine and to integrate the model into a comprehensive drug-drug(-gene) interaction (DD(G)I) network with a diverse set of CYP3A4 and P-gp perpetrators as well as CYP2D6 and P-gp victims. The quinidine parent-metabolite model including 3-hydroxyquinidine was developed using pharmacokinetic profiles from clinical studies after intravenous and oral administration covering a broad dosing range (0.1-600 mg). The model covers efflux transport via P-gp and metabolic transformation to either 3-hydroxyquinidine or unspecified metabolites via CYP3A4. The 3-hydroxyquinidine model includes further metabolism by CYP3A4 as well as an unspecific hepatic clearance. Model performance was assessed graphically and quantitatively with greater than 90% of predicted pharmacokinetic parameters within two-fold of corresponding observed values. The model was successfully used to simulate various DD(G)I scenarios with greater than 90% of predicted DD(G)I pharmacokinetic parameter ratios within two-fold prediction success limits. The presented network will be provided to the research community and can be extended to include further perpetrators, victims, and targets, to support investigations of DD(G)Is.

Physiologically based pharmacokinetic modeling of tacrolimus for food-drug and CYP3A drug-drug-gene interaction predictions.

The immunosuppressant and narrow therapeutic index drug tacrolimus is metabolized mainly via cytochrome P450 (CYP) 3A4 and CYP3A5. For its pharmacokinetics (PK), high inter- and intra-individual variability can be observed. Underlying causes include the effect of food intake on tacrolimus absorption as well as genetic polymorphism in the CYP3A5 gene. Furthermore, tacrolimus is highly susceptible to drug-drug interactions, acting as a victim drug when coadministered with CYP3A perpetrators. This work describes the development of a whole-body physiologically based pharmacokinetic model for tacrolimus as well as its application for investigation and prediction of (i) the impact of food intake on tacrolimus PK (food-drug interactions [FDIs]) and (ii) drug-drug(-gene) interactions (DD[G]Is) involving the CYP3A perpetrator drugs voriconazole, itraconazole, and rifampicin. The model was built in PK-Sim® Version 10 using a total of 37 whole blood concentration-time profiles of tacrolimus (training and test) compiled from 911 healthy individuals covering the administration of tacrolimus as intravenous infusions as well as immediate-release and extended-release capsules. Metabolism was incorporated via CYP3A4 and CYP3A5, with varying activities implemented for different CYP3A5 genotypes and study populations. The good predictive model performance is demonstrated for the examined food effect studies with 6/6 predicted FDI area under the curve determined between first and last concentration measurements (AUClast ) and 6/6 predicted FDI maximum whole blood concentration (Cmax ) ratios within twofold of the respective observed ratios. In addition, 7/7 predicted DD(G)I AUClast and 6/7 predicted DD(G)I Cmax ratios were within twofold of their observed values. Potential applications of the final model include model-informed drug discovery and development or the support of model-informed precision dosing.

Prediction of Drug-Drug-Gene Interaction Scenarios of (E)-Clomiphene and Its Metabolites Using Physiologically Based Pharmacokinetic Modeling.

Clomiphene, a selective estrogen receptor modulator (SERM), has been used for the treatment of anovulation for more than 50 years. However, since (E)-clomiphene ((E)-Clom) and its metabolites are eliminated primarily via Cytochrome P450 (CYP) 2D6 and CYP3A4, exposure can be affected by CYP2D6 polymorphisms and concomitant use with CYP inhibitors. Thus, clomiphene therapy may be susceptible to drug-gene interactions (DGIs), drug-drug interactions (DDIs) and drug-drug-gene interactions (DDGIs). Physiologically based pharmacokinetic (PBPK) modeling is a tool to quantify such DGI and DD(G)I scenarios. This study aimed to develop a whole-body PBPK model of (E)-Clom including three important metabolites to describe and predict DGI and DD(G)I effects. Model performance was evaluated both graphically and by calculating quantitative measures. Here, 90% of predicted Cmax and 80% of AUClast values were within two-fold of the corresponding observed value for DGIs and DD(G)Is with clarithromycin and paroxetine. The model also revealed quantitative contributions of different CYP enzymes to the involved metabolic pathways of (E)-Clom and its metabolites. The developed PBPK model can be employed to assess the exposure of (E)-Clom and its active metabolites in as-yet unexplored DD(G)I scenarios in future studies.

Physiologically Based Pharmacokinetic Modeling to Describe the CYP2D6 Activity Score-Dependent Metabolism of Paroxetine, Atomoxetine and Risperidone.

The cytochrome P450 2D6 (CYP2D6) genotype is the single most important determinant of CYP2D6 activity as well as interindividual and interpopulation variability in CYP2D6 activity. Here, the CYP2D6 activity score provides an established tool to categorize the large number of CYP2D6 alleles by activity and facilitates the process of genotype-to-phenotype translation. Compared to the broad traditional phenotype categories, the CYP2D6 activity score additionally serves as a superior scale of CYP2D6 activity due to its finer graduation. Physiologically based pharmacokinetic (PBPK) models have been successfully used to describe and predict the activity score-dependent metabolism of CYP2D6 substrates. This study aimed to describe CYP2D6 drug-gene interactions (DGIs) of important CYP2D6 substrates paroxetine, atomoxetine and risperidone by developing a substrate-independent approach to model their activity score-dependent metabolism. The models were developed in PK-Sim®, using a total of 57 plasma concentration-time profiles, and showed good performance, especially in DGI scenarios where 10/12, 5/5 and 7/7 of DGI AUClast ratios and 9/12, 5/5 and 7/7 of DGI Cmax ratios were within the prediction success limits. Finally, the models were used to predict their compound's exposure for different CYP2D6 activity scores during steady state. Here, predicted DGI AUCss ratios were 3.4, 13.6 and 2.0 (poor metabolizers; activity score = 0) and 0.2, 0.5 and 0.95 (ultrarapid metabolizers; activity score = 3) for paroxetine, atomoxetine and risperidone active moiety (risperidone + 9-hydroxyrisperidone), respectively.

Physiologically-based pharmacokinetic modeling of dextromethorphan to investigate interindividual variability within CYP2D6 activity score groups.

This study provides a whole-body physiologically-based pharmacokinetic (PBPK) model of dextromethorphan and its metabolites dextrorphan and dextrorphan O-glucuronide for predicting the effects of cytochrome P450 2D6 (CYP2D6) drug-gene interactions (DGIs) on dextromethorphan pharmacokinetics (PK). Moreover, the effect of interindividual variability (IIV) within CYP2D6 activity score groups on the PK of dextromethorphan and its metabolites was investigated. A parent-metabolite-metabolite PBPK model of dextromethorphan, dextrorphan, and dextrorphan O-glucuronide was developed in PK-Sim and MoBi. Drug-dependent parameters were obtained from the literature or optimized. Plasma concentration-time profiles of all three analytes were gathered from published studies and used for model development and model evaluation. The model was evaluated comparing simulated plasma concentration-time profiles, area under the concentration-time curve from the time of the first measurement to the time of the last measurement (AUClast ) and maximum concentration (Cmax ) values to observed study data. The final PBPK model accurately describes 28 population plasma concentration-time profiles and plasma concentration-time profiles of 72 individuals from four cocktail studies. Moreover, the model predicts CYP2D6 DGI scenarios with six of seven DGI AUClast and seven of seven DGI Cmax ratios within the acceptance criteria. The high IIV in plasma concentrations was analyzed by characterizing the distribution of individually optimized CYP2D6 kcat values stratified by activity score group. Population simulations with sampling from the resulting distributions with calculated log-normal dispersion and mean parameters could explain a large extent of the observed IIV. The model is publicly available alongside comprehensive documentation of model building and model evaluation.

Novel models for the prediction of drug-gene interactions.

INTRODUCTION: Adverse drug reactions (ADRs) are among the leading causes of death, and frequently associated with drug-gene interactions (DGIs). In addition to pharmacogenomic programs for implementation of genetic preemptive testing into clinical practice, mathematical modeling can help to understand, quantify and predict the effects of DGIs in vivo. Moreover, modeling can contribute to optimize prospective clinical drug trial activities and to reduce DGI-related ADRs. AREAS COVERED: Approaches and challenges of mechanistical DGI implementation and model parameterization are discussed for population pharmacokinetic and physiologically based pharmacokinetic models. The broad spectrum of published DGI models and their applications is presented, focusing on the investigation of DGI effects on pharmacology and model-based dose adaptations. EXPERT OPINION: Mathematical modeling provides an opportunity to investigate complex DGI scenarios and can facilitate the development process of safe and efficient personalized dosing regimens. However, reliable DGI model input data from in vivo and in vitro measurements are crucial. For this, collaboration among pharmacometricians, laboratory scientists and clinicians is important to provide homogeneous datasets and unambiguous model parameters. For a broad adaptation of validated DGI models in clinical practice, interdisciplinary cooperation should be promoted and qualification toolchains must be established.

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.