Search

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

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

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

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

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

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.

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

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

Infliximab is approved for treatment of various chronic inflammatory diseases including inflammatory bowel disease (IBD). However, high variability in infliximab trough levels has been associated with diverse response rates. Model-informed precision dosing (MIPD) with population pharmacokinetic models could help to individualize infliximab dosing regimens and improve therapy. The aim of this study was to evaluate the predictive performance of published infliximab population pharmacokinetic models for IBD patients with an external data set. The data set consisted of 105 IBD patients with 336 infliximab concentrations. Literature review identified 12 published models eligible for external evaluation. Model performance was evaluated with goodness-of-fit plots, prediction- and variability-corrected visual predictive checks (pvcVPCs) and quantitative measures. For anti-drug antibody (ADA)-negative patients, model accuracy decreased for predictions > 6 months, while bias did not increase. In general, predictions for patients developing ADA were less accurate for all models investigated. Two models with the highest classification accuracy identified necessary dose escalations (for trough concentrations < 5 µg/mL) in 88% of cases. In summary, population pharmacokinetic modeling can be used to individualize infliximab dosing and thereby help to prevent infliximab trough concentrations dropping below the target trough concentration. However, predictions of infliximab concentrations for patients developing ADA remain challenging.

A better understanding of their interaction with cell-based tissue is a fundamental prerequisite towards the safe production and application of engineered nanomaterials. Quantitative experimental data on the correlation between physicochemical characteristics and the interaction and transport of engineered nanomaterials across biological barriers, in particular, is still scarce, thus hampering the development of effective predictive non-testing strategies. Against this background, the presented study investigated the translocation of gold and silver nanoparticles across the gastrointestinal barrier along with related biological effects using an in vitro 3D-triple co-culture cell model. Standardized in vitro assays and quantitative polymerase chain reaction showed no significant influence of the applied nanoparticles on both cell viability and generation of reactive oxygen species. Transmission electron microscopy indicated an intact cell barrier during the translocation study. Single particle ICP-MS revealed a time-dependent increase of translocated nanoparticles independent of their size, shape, surface charge, and stability in cell culture medium. This quantitative data provided the experimental basis for the successful mathematical description of the nanoparticle transport kinetics using a non-linear mixed effects modeling approach. The results of this study may serve as a basis for the development of predictive tools for improved risk assessment of engineered nanomaterials in the future.

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.

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.

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.

Fentanyl is widely used for analgesia, sedation, and anesthesia both in adult and pediatric populations. Yet, only few pharmacokinetic studies of fentanyl in pediatrics exist as conducting clinical trials in this population is especially challenging. Physiologically-based pharmacokinetic (PBPK) modeling is a mechanistic approach to explore drug pharmacokinetics and allows extrapolation from adult to pediatric populations based on age-related physiological differences. The aim of this study was to develop a PBPK model of fentanyl and norfentanyl for both adult and pediatric populations. The adult PBPK model was established in PK-Sim® using data from 16 clinical studies and was scaled to several pediatric subpopulations. ~93% of the predicted AUClast values in adults and ~88% in pediatrics were within 2-fold of the corresponding value observed. The adult PBPK model predicted a fraction of fentanyl dose metabolized to norfentanyl of ~33% and a fraction excreted in urine of ~7%. In addition, the pediatric PBPK model was used to simulate differences in peak plasma concentrations after bolus injections and short infusions. The novel PBPK models could be helpful to further investigate fentanyl pharmacokinetics in both adult and pediatric populations.

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.

Buprenorphine plays a crucial role in the therapeutic management of pain in adults, adolescents and pediatric subpopulations. However, only few pharmacokinetic studies of buprenorphine in children, particularly neonates, are available as conducting clinical trials in this population is especially challenging. Physiologically-based pharmacokinetic (PBPK) modeling allows the prediction of drug exposure in pediatrics based on age-related physiological differences. The aim of this study was to predict the pharmacokinetics of buprenorphine in pediatrics with PBPK modeling. Moreover, the drug-drug interaction (DDI) potential of buprenorphine with CYP3A4 and P-glycoprotein perpetrator drugs should be elucidated. A PBPK model of buprenorphine and norbuprenorphine in adults has been developed and scaled to children and preterm neonates, accounting for age-related changes. One-hundred-percent of the predicted AUClast values in adults (geometric mean fold error (GMFE): 1.22), 90% of individual AUClast predictions in children (GMFE: 1.54) and 75% in preterm neonates (GMFE: 1.57) met the 2-fold acceptance criterion. Moreover, the adult model was used to simulate DDI scenarios with clarithromycin, itraconazole and rifampicin. We demonstrate the applicability of scaling adult PBPK models to pediatrics for the prediction of individual plasma profiles. The novel PBPK models could be helpful to further investigate buprenorphine pharmacokinetics in various populations, particularly pediatric subgroups.