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.

Personalized Chronomodulated 5-Fluorouracil Treatment: A Physiologically-Based Pharmacokinetic Precision Dosing Approach for Optimizing Cancer Therapy.

The discovery of circadian clock genes greatly amplified the study of diurnal variations impacting cancer therapy, transforming it into a rapidly growing field of research. Especially, use of chronomodulated treatment with 5-fluorouracil (5-FU) has gained significance. Studies indicate high interindividual variability (IIV) in diurnal variations in dihydropyrimidine dehydrogenase (DPD) activity - a key enzyme for 5-FU metabolism. However, the influence of individual DPD chronotypes on chronomodulated therapy remains unclear and warrants further investigation. To optimize precision dosing of chronomodulated 5-FU, this study aims to: (i) build physiologically-based pharmacokinetic (PBPK) models for 5-FU, uracil, and their metabolites, (ii) assess the impact of diurnal variation on DPD activity, (iii) estimate individual DPD chronotypes, and (iv) personalize chronomodulated 5-FU infusion rates based on a patient's DPD chronotype. Whole-body PBPK models were developed with PK-Sim(R) and MoBi(R). Sinusoidal functions were used to incorporate variations in enzyme activity and chronomodulated infusion rates as well as to estimate individual DPD chronotypes from DPYD mRNA expression or DPD enzymatic activity. Four whole-body PBPK models for 5-FU, uracil, and their metabolites were established utilizing data from 41 5-FU and 10 publicly available uracil studies. IIV in DPD chronotypes was assessed and personalized chronomodulated administrations were developed to achieve (i) comparable 5-FU peak plasma concentrations, (ii) comparable 5-FU exposure, and (iii) constant 5-FU plasma levels via "noise cancellation" chronomodulated infusion. The developed PBPK models capture the extent of diurnal variations in DPD activity and can help investigate individualized chronomodulated 5-FU therapy through testing alternative personalized dosing strategies.

A dynamic time-to-event model for prediction of acute graft-versus-host disease in patients after allogeneic hematopoietic stem cell transplantation.

BACKGROUND: Acute graft-versus-host disease (aGvHD) is a major cause of death for patients following allogeneic hematopoietic stem cell transplantation (HSCT). Effective management of moderate to severe aGvHD remains challenging despite recent advances in HSCT, emphasizing the importance of prophylaxis and risk factor identification. METHODS: In this study, we analyzed data from 1479 adults who underwent HSCT between 2005 and 2017 to investigate the effects of aGvHD prophylaxis and time-dependent risk factors on the development of grades II-IV aGvHD within 100 days post-HSCT. RESULTS: Using a dynamic longitudinal time-to-event model, we observed a non-monotonic baseline hazard overtime with a low hazard during the first few days and a maximum hazard at day 17, described by Bateman function with a mean transit time of approximately 11 days. Multivariable analysis revealed significant time-dependent effects of white blood cell counts and cyclosporine A exposure as well as static effects of female donors for male recipients, patients with matched related donors, conditioning regimen consisting of fludarabine plus total body irradiation, and patient age in recipients of grafts from related donors on the risk to develop grades II-IV aGvHD. Additionally, we found that higher cumulative hazard on day 7 after allo-HSCT are associated with an increased incidence of grades II-IV aGvHD within 100 days indicating that an individual assessment of the cumulative hazard on day 7 could potentially serve as valuable predictor for later grades II-IV aGvHD development. Using the final model, stochastic simulations were performed to explore covariate effects on the cumulative incidence over time and to estimate risk ratios. CONCLUSION: Overall, the presented model showed good descriptive and predictive performance and provides valuable insights into the interplay of multiple static and time-dependent risk factors for the prediction of aGvHD.

Motivations for Adolescent COVID-19 Vaccination: A Comparative Study of Adolescent and Caregiver Perspectives in Germany.

Given the crucial role of vaccination in halting the COVID-19 pandemic, it is imperative to understand the factors that motivate adolescents to get vaccinated. We surveyed adolescents and their accompanying guardians scheduled to receive a COVID-19 vaccination (Comirnaty) in an urban region in Germany in mid-2021 regarding their motivation for getting vaccinated and collected data on their sociodemographic characteristics, medical history, vaccination status, and any history of COVID-19 infection in the family. We also queried information strategies related to the SARS-CoV-2 pandemic. Motivations for getting vaccinated were similar among adolescents and their parents. The primary reasons for vaccination were protection against SARS-CoV-2-related illness and gaining access to leisure facilities. This was not influenced by gender, health status, migration background, or the presence of chronic or acute diseases. The percentage of parents who had received SARS-CoV-2 immunization and the proportion of parents with a high level of education were higher among study participants than in the general population. Adolescents were especially willing to be vaccinated if they came from a better educational environment and had a high vaccination rate in the family. Emphasizing the importance of vaccination among all segments of the population and removing barriers to vaccines may lead to an ameliorated acceptance of COVID-19 vaccines.

Physiologically Based Pharmacokinetic Modeling of Bergamottin and 6,7-Dihydroxybergamottin to Describe CYP3A4 Mediated Grapefruit-Drug Interactions.

Grapefruit is a moderate to strong inactivator of CYP3A4, which metabolizes up to 50% of marketed drugs. The inhibitory effect is mainly attributed to furanocoumarins present in the fruit, irreversibly inhibiting preferably intestinal CYP3A4 as suicide inhibitors. Effects on CYP3A4 victim drugs can still be measured up to 24 hours after grapefruit juice (GFJ) consumption. The current study aimed to establish a physiologically-based pharmacokinetic (PBPK) grapefruit-drug interaction model by modeling the relevant CYP3A4 inhibiting ingredients of the fruit to simulate and predict the effect of GFJ consumption on plasma concentration-time profiles of various CYP3A4 victim drugs. The grapefruit model was developed in PK-Sim and coupled with previously developed PBPK models of CYP3A4 substrates that were publicly available and already evaluated for CYP3A4-mediated drug-drug interactions. Overall, 43 clinical studies were used for model development. Models of bergamottin (BGT) and 6,7-dihydroxybergamottin (DHB) as relevant active ingredients in GFJ were established. Both models include: (i) CYP3A4 inactivation informed by in vitro parameters, (ii) a CYP3A4 mediated clearance estimated during model development, as well as (iii) passive glomerular filtration. The final model successfully describes interactions of GFJ ingredients with 10 different CYP3A4 victim drugs, simulating the effect of the CYP3A4 inactivation on the victims' pharmacokinetics as well as their main metabolites. Furthermore, the model sufficiently captures the time-dependent effect of CYP3A4 inactivation as well as the effect of grapefruit ingestion on intestinal and hepatic CYP3A4 concentrations.

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.

Significant impact of time-of-day variation on metformin pharmacokinetics.

AIMS/HYPOTHESIS: The objective was to investigate if metformin pharmacokinetics is modulated by time-of-day in humans using empirical and mechanistic pharmacokinetic modelling techniques on a large clinical dataset. This study also aimed to generate and test hypotheses on the underlying mechanisms, including evidence for chronotype-dependent interindividual differences in metformin plasma and efficacy-related tissue concentrations. METHODS: A large clinical dataset consisting of individual metformin plasma and urine measurements was analysed using a newly developed empirical pharmacokinetic model. Causes of daily variation of metformin pharmacokinetics and interindividual variability were further investigated by a literature-informed mechanistic modelling analysis. RESULTS: A significant effect of time-of-day on metformin pharmacokinetics was found. Daily rhythms of gastrointestinal, hepatic and renal processes are described in the literature, possibly affecting drug pharmacokinetics. Observed metformin plasma levels were best described by a combination of a rhythm in GFR, renal plasma flow (RPF) and organic cation transporter (OCT) 2 activity. Furthermore, the large interindividual differences in measured metformin concentrations were best explained by individual chronotypes affecting metformin clearance, with impact on plasma and tissue concentrations that may have implications for metformin efficacy. CONCLUSIONS/INTERPRETATION: Metformin's pharmacology significantly depends on time-of-day in humans, determined with the help of empirical and mechanistic pharmacokinetic modelling, and rhythmic GFR, RPF and OCT2 were found to govern intraday variation. Interindividual variation was found to be partly dependent on individual chronotype, suggesting diurnal preference as an interesting, but so-far underappreciated, topic with regard to future personalised chronomodulated therapy in people with type 2 diabetes.

A Physiologically Based Pharmacokinetic Model of Ketoconazole and Its Metabolites as Drug-Drug Interaction Perpetrators.

The antifungal ketoconazole, which is mainly used for dermal infections and treatment of Cushing's syndrome, is prone to drug-food interactions (DFIs) and is well known for its strong drug-drug interaction (DDI) potential. Some of ketoconazole's potent inhibitory activity can be attributed to its metabolites that predominantly accumulate in the liver. This work aimed to develop a whole-body physiologically based pharmacokinetic (PBPK) model of ketoconazole and its metabolites for fasted and fed states and to investigate the impact of ketoconazole's metabolites on its DDI potential. The parent-metabolites model was developed with PK-Sim® and MoBi® using 53 plasma concentration-time profiles. With 7 out of 7 (7/7) DFI AUClast and DFI Cmax ratios within two-fold of observed ratios, the developed model demonstrated good predictive performance under fasted and fed conditions. DDI scenarios that included either the parent alone or with its metabolites were simulated and evaluated for the victim drugs alfentanil, alprazolam, midazolam, triazolam, and digoxin. DDI scenarios that included all metabolites as reversible inhibitors of CYP3A4 and P-gp performed best: 26/27 of DDI AUClast and 21/21 DDI Cmax ratios were within two-fold of observed ratios, while DDI models that simulated only ketoconazole as the perpetrator underperformed: 12/27 DDI AUClast and 18/21 DDI Cmax ratios were within the success limits.

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.

Pharmacometric Modeling of the Impact of Azelastine Nasal Spray on SARS-CoV-2 Viral Load and Related Symptoms in COVID-19 Patients.

The histamine-1 receptor antagonist azelastine was recently found to impact SARS-CoV-2 viral kinetics in a Phase 2 clinical trial (CARVIN). Thus, we investigated the relationship between intranasal azelastine administrations and viral load, as well as symptom severity in COVID-19 patients and analyzed the impact of covariates using non-linear mixed-effects modeling. For this, we developed a pharmacokinetic (PK) model for the oral and intranasal administration of azelastine. A one-compartment model with parallel absorption after intranasal administration described the PK best, covering both the intranasal and the gastro-intestinal absorption pathways. For virus kinetic and symptoms modeling, viral load and symptom records were gathered from the CARVIN study that included data of 82 COVID-19 patients receiving placebo or intranasal azelastine. The effect of azelastine on viral load was described by a dose-effect model targeting the virus elimination rate. An extension of the model revealed a relationship between COVID-19 symptoms severity and the number of infected cells. The analysis revealed that the intranasal administration of azelastine led to a faster decline in viral load and symptoms severity compared to placebo. Moreover, older patients showed a slower decline in viral load compared to younger patients and male patients experienced higher peak viral loads than females.

Model-Based Analysis of SARS-CoV-2 Infections, Hospitalization and Outcome in Germany, the Federal States and Districts.

The coronavirus disease 2019 (COVID-19) pandemic challenged many national health care systems, with hospitals reaching capacity limits of intensive care units (ICU). Thus, the estimation of acute local burden of ICUs is critical for appropriate management of health care resources. In this work, we applied non-linear mixed effects modeling to develop an epidemiological SARS-CoV-2 infection model for Germany, with its 16 federal states and 400 districts, that describes infections as well as COVID-19 inpatients, ICU patients with and without mechanical ventilation, recoveries, and fatalities during the first two waves of the pandemic until April 2021. Based on model analyses, covariates influencing the relation between infections and outcomes were explored. Non-pharmaceutical interventions imposed by governments were found to have a major impact on the spreading of SARS-CoV-2. Patient age and sex, the spread of variant B.1.1.7, and the testing strategy (number of tests performed weekly, rate of positive tests) affected the severity and outcome of recorded cases and could reduce the observed unexplained variability between the states. Modeling could reasonably link the discrepancies between fine-grained model simulations of the 400 German districts and the reported number of available ICU beds to coarse-grained COVID-19 patient distribution patterns within German regions.

Human papillomavirus vaccination of girls in the German model region Saarland: Insurance data-based analysis and identification of starting points for improving vaccination rates.

In Germany, the incidence of cervical cancer, a disease caused by human papillomaviruses (HPV), is higher than in neighboring European countries. HPV vaccination has been recommended for girls since 2007. However, it continues to be significantly less well received than other childhood vaccines, so its potential for cancer prevention is not fully realized. To find new starting points for improving vaccination rates, we analyzed pseudonymized routine billing data from statutory health insurers in the PRÄZIS study (prevention of cervical carcinoma and its precursors in women in Saarland) in the federal state Saarland serving as a model region. We show that lowering the HPV vaccination age to 9 years led to more completed HPV vaccinations already in 2015. Since then, HPV vaccination rates and the proportion of 9- to 11-year-old girls among HPV-vaccinated females have steadily increased. However, HPV vaccination rates among 15-year-old girls in Saarland remained well below 50% in 2019. Pediatricians vaccinated the most girls overall, with a particularly high proportion at the recommended vaccination age of 9-14 years, while gynecologists provided more HPV catch-up vaccinations among 15-17-year-old girls, and general practitioners compensated for HPV vaccination in Saarland communities with fewer pediatricians or gynecologists. We also provide evidence for a significant association between attendance at the children´s medical check-ups "U11" or "J1" and HPV vaccination. In particular, participation in HPV vaccination is high on the day of U11. However, obstacles are that U11 is currently not financed by all statutory health insurers and there is a lack of invitation procedures for both U11 and J1, resulting in significantly lower participation rates than for the earlier U8 or U9 screenings, which are conducted exclusively with invitations and reminders. Based on our data, we propose to restructure U11 and J1 screening in Germany, with mandatory funding for U11 and organized invitations for HPV vaccination at U11 or J1 for both boys and girls.

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.

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

Renal Transporter-Mediated Drug-Biomarker Interactions of the Endogenous Substrates Creatinine and N1 -Methylnicotinamide: A PBPK Modeling Approach.

Endogenous biomarkers for transporter-mediated drug-drug interaction (DDI) predictions represent a promising approach to facilitate and improve conventional DDI investigations in clinical studies. This approach requires high sensitivity and specificity of biomarkers for the targets of interest (e.g., transport proteins), as well as rigorous characterization of their kinetics, which can be accomplished utilizing physiologically-based pharmacokinetic (PBPK) modeling. Therefore, the objective of this study was to develop PBPK models of the endogenous organic cation transporter (OCT)2 and multidrug and toxin extrusion protein (MATE)1 substrates creatinine and N1 -methylnicotinamide (NMN). Additionally, this study aimed to predict kinetic changes of the biomarkers during administration of the OCT2 and MATE1 perpetrator drugs trimethoprim, pyrimethamine, and cimetidine. Whole-body PBPK models of creatinine and NMN were developed utilizing studies investigating creatinine or NMN exogenous administration and endogenous synthesis. The newly developed models accurately describe and predict observed plasma concentration-time profiles and urinary excretion of both biomarkers. Subsequently, models were coupled to the previously built and evaluated perpetrator models of trimethoprim, pyrimethamine, and cimetidine for interaction predictions. Increased creatinine plasma concentrations and decreased urinary excretion during the drug-biomarker interactions with trimethoprim, pyrimethamine, and cimetidine were well-described. An additional inhibition of NMN synthesis by trimethoprim and pyrimethamine was hypothesized, improving NMN plasma and urine interaction predictions. To summarize, whole-body PBPK models of creatinine and NMN were built and evaluated to better assess creatinine and NMN kinetics while uncovering knowledge gaps for future research. The models can support investigations of renal transporter-mediated DDIs during drug development.

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.

In Vitro-In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier.

Static in vitro permeation experiments are commonly used to gain insights into the permeation properties of drug substances but exhibit limitations due to missing physiologic cell stimuli. Thus, fluidic systems integrating stimuli, such as physicochemical fluxes, have been developed. However, as fluidic in vitro studies display higher complexity compared to static systems, analysis of experimental readouts is challenging. Here, the integration of in silico tools holds the potential to evaluate fluidic experiments and to investigate specific simulation scenarios. This study aimed to develop in silico models that describe and predict the permeation and disposition of two model substances in a static and fluidic in vitro system. For this, in vitro permeation studies with a 16HBE cellular barrier under both static and fluidic conditions were performed over 72 h. In silico models were implemented and employed to describe and predict concentration-time profiles of caffeine and diclofenac in various experimental setups. For both substances, in silico modeling identified reduced apparent permeabilities in the fluidic compared to the static cellular setting. The developed in vitro-in silico modeling framework can be expanded further, integrating additional cell tissues in the fluidic system, and can be employed in future studies to model pharmacokinetic and pharmacodynamic drug behavior.