Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans.

Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration-time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration-time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.

Computational identification of natural senotherapeutic compounds that mimic dasatinib based on gene expression data.

The major risk factor for chronic disease is chronological age, and age-related chronic diseases account for the majority of deaths worldwide. Targeting senescent cells that accumulate in disease-related tissues presents a strategy to reduce disease burden and to increase healthspan. The senolytic combination of the tyrosine-kinase inhibitor dasatinib and the flavonol quercetin is frequently used in clinical trials aiming to eliminate senescent cells. Here, our goal was to computationally identify natural senotherapeutic repurposing candidates that may substitute dasatinib based on their similarity in gene expression effects. The natural senolytic piperlongumine (a compound found in long pepper), and the natural senomorphics parthenolide, phloretin and curcumin (found in various edible plants) were identified as potential substitutes of dasatinib. The gene expression changes underlying the repositioning highlight apoptosis-related genes and pathways. The four compounds, and in particular the top-runner piperlongumine, may be combined with quercetin to obtain natural formulas emulating the dasatinib + quercetin formula.

Increased sinusoidal export of drug glucuronides is a compensative mechanism in liver cirrhosis of mice.

Rationale: Liver cirrhosis is known to affect drug pharmacokinetics, but the functional assessment of the underlying pathophysiological alterations in drug metabolism is difficult. Methods: Cirrhosis in mice was induced by repeated treatment with carbon tetrachloride for 12 months. A cocktail of six drugs was administered, and parent compounds as well as phase I and II metabolites were quantified in blood, bile, and urine in a time-dependent manner. Pharmacokinetics were modeled in relation to the altered expression of metabolizing enzymes. In discrepancy with computational predictions, a strong increase of glucuronides in blood was observed in cirrhotic mice compared to vehicle controls. Results: The deviation between experimental findings and computational simulations observed by analyzing different hypotheses could be explained by increased sinusoidal export and corresponded to increased expression of export carriers (Abcc3 and Abcc4). Formation of phase I metabolites and clearance of the parent compounds were surprisingly robust in cirrhosis, although the phase I enzymes critical for the metabolism of the administered drugs in healthy mice, Cyp1a2 and Cyp2c29, were downregulated in cirrhotic livers. RNA-sequencing revealed the upregulation of numerous other phase I metabolizing enzymes which may compensate for the lost CYP isoenzymes. Comparison of genome-wide data of cirrhotic mouse and human liver tissue revealed similar features of expression changes, including increased sinusoidal export and reduced uptake carriers. Conclusion: Liver cirrhosis leads to increased blood concentrations of glucuronides because of increased export from hepatocytes into the sinusoidal blood. Although individual metabolic pathways are massively altered in cirrhosis, the overall clearance of the parent compounds was relatively robust due to compensatory mechanisms.

A physiologically based model of bile acid metabolism in mice.

Bile acid (BA) metabolism is a complex system that includes a wide variety of primary and secondary, as well as conjugated and unconjugated BAs that undergo continuous enterohepatic circulation (EHC). Alterations in both composition and dynamics of BAs have been associated with various diseases. However, a mechanistic understanding of the relationship between altered BA metabolism and related diseases is lacking. Computational modeling may support functional analyses of the physiological processes involved in the EHC of BAs along the gut-liver axis. In this study, we developed a physiologically based model of murine BA metabolism describing synthesis, hepatic and microbial transformations, systemic distribution, excretion, and EHC of BAs at the whole-body level. For model development, BA metabolism of specific pathogen-free (SPF) mice was characterized in vivo by measuring BA levels and composition in various organs, expression of transporters along the gut, and cecal microbiota composition. We found significantly different BA levels between male and female mice that could only be explained by adjusted expression of the hepatic enzymes and transporters in the model. Of note, this finding was in agreement with experimental observations. The model for SPF mice could also describe equivalent experimental data in germ-free mice by specifically switching off microbial activity in the intestine. The here presented model can therefore facilitate and guide functional analyses of BA metabolism in mice, e.g., the effect of pathophysiological alterations on BA metabolism and translation of results from mouse studies to a clinically relevant context through cross-species extrapolation.

Quantitative systems pharmacology of the eye: Tools and data for ocular QSP.

Good eyesight belongs to the most-valued attributes of health, and diseases of the eye are a significant healthcare burden. Case numbers are expected to further increase in the next decades due to an aging society. The development of drugs in ophthalmology, however, is difficult due to limited accessibility of the eye, in terms of drug administration and in terms of sampling of tissues for drug pharmacokinetics (PKs) and pharmacodynamics (PDs). Ocular quantitative systems pharmacology models provide the opportunity to describe the distribution of drugs in the eye as well as the resulting drug-response in specific segments of the eye. In particular, ocular physiologically-based PK (PBPK) models are necessary to describe drug concentration levels in different regions of the eye. Further, ocular effect models using molecular data from specific cellular systems are needed to develop dose-response correlations. We here describe the current status of PK/PBPK as well as PD models for the eyes and discuss cellular systems, data repositories, as well as animal models in ophthalmology. The application of the various concepts is highlighted for the development of new treatments for postoperative fibrosis after glaucoma surgery.

Blood Endotoxin Levels as Biomarker of Nonalcoholic Fatty Liver Disease: A Systematic Review and Meta-analysis.

BACKGROUND & AIMS: Growing evidence supports a role of gut-derived metabolites in nonalcoholic fatty liver disease (NAFLD), but the relation of endotoxin levels with gut permeability and NAFLD stage remains unclear. This systematic review with meta-analysis aims to provide further insights. METHODS: PubMed, Embase, and Cochrane Library were searched for studies published until January 2022 assessing blood endotoxins in patients with NAFLD. Meta-analyses and univariate/multivariate meta-regression, as well as correlation analyses, were performed for endotoxin values and potential relationships to disease stage, age, sex, parameters of systemic inflammation, and metabolic syndrome, as well as liver function and histology. RESULTS: Forty-three studies were included, of which 34 were used for meta-analyses. Blood endotoxin levels were higher in patients with simple steatosis vs liver-healthy controls (standardized mean difference, 0.86; 95% confidence interval, 0.62-1.11) as well as in patients with nonalcoholic steatohepatitis vs patients with nonalcoholic fatty liver/non-nonalcoholic steatohepatitis (standardized mean difference, 0.81; 95% confidence interval, 0.27-1.35; P = .0078). Consistently, higher endotoxin levels were observed in patients with more advanced histopathological gradings of liver steatosis and fibrosis. An increase of blood endotoxin levels was partially attributed to a body mass index rise in patients with NAFLD compared with controls. Nevertheless, significant increases of blood endotoxin levels in NAFLD retained after compensation for differences in body mass index, metabolic condition, or liver enzymes. Increases in blood endotoxin levels were associated with increases in C-reactive protein concentrations, and in most cases, paralleled a rise in markers for intestinal permeability. CONCLUSION: Our results support blood endotoxin levels as relevant diagnostic biomarker for NAFLD, both for disease detection as well as staging during disease progression, and might serve as surrogate marker of enhanced intestinal permeability in NAFLD. Registration number in Prospero: CRD42022311166.

Computational Models for Clinical Applications in Personalized Medicine-Guidelines and Recommendations for Data Integration and Model Validation.

The future development of personalized medicine depends on a vast exchange of data from different sources, as well as harmonized integrative analysis of large-scale clinical health and sample data. Computational-modelling approaches play a key role in the analysis of the underlying molecular processes and pathways that characterize human biology, but they also lead to a more profound understanding of the mechanisms and factors that drive diseases; hence, they allow personalized treatment strategies that are guided by central clinical questions. However, despite the growing popularity of computational-modelling approaches in different stakeholder communities, there are still many hurdles to overcome for their clinical routine implementation in the future. Especially the integration of heterogeneous data from multiple sources and types are challenging tasks that require clear guidelines that also have to comply with high ethical and legal standards. Here, we discuss the most relevant computational models for personalized medicine in detail that can be considered as best-practice guidelines for application in clinical care. We define specific challenges and provide applicable guidelines and recommendations for study design, data acquisition, and operation as well as for model validation and clinical translation and other research areas.

Early prediction of decompensation (EPOD) score: Non-invasive determination of cirrhosis decompensation risk.

BACKGROUND & AIMS: Decompensation is a hallmark of disease progression in cirrhotic patients. Early detection of a phase transition from compensated cirrhosis to decompensation would enable targeted therapeutic interventions potentially extending life expectancy. This study aims to (a) identify the predictors of decompensation in a large, multicentric cohort of patients with compensated cirrhosis, (b) to build a reliable prognostic score for decompensation and (c) to evaluate the score in independent cohorts. METHODS: Decompensation was identified in electronic health records data from 6049 cirrhosis patients in the IBM Explorys database training cohort by diagnostic codes for variceal bleeding, encephalopathy, ascites, hepato-renal syndrome and/or jaundice. We identified predictors of clinical decompensation and developed a prognostic score using Cox regression analysis. The score was evaluated using the IBM Explorys database validation cohort (N = 17662), the Penn Medicine BioBank (N = 1326) and the UK Biobank (N = 317). RESULTS: The new Early Prediction of Decompensation (EPOD) score uses platelet count, albumin, and bilirubin concentration. It predicts decompensation during a 3-year follow-up in three validation cohorts with AUROCs of 0.69, 0.69 and 0.77, respectively, and outperforms the well-known MELD and Child-Pugh score in predicting decompensation. Furthermore, the EPOD score predicted the 3-year probability of decompensation. CONCLUSIONS: The EPOD score provides a prediction tool for the risk of decompensation in patients with cirrhosis that outperforms well-known cirrhosis scores. Since EPOD is based on three blood parameters, only, it provides maximal clinical feasibility at minimal costs.

Algorithmic surveillance of ICU patients with acute respiratory distress syndrome (ASIC): protocol for a multicentre stepped-wedge cluster randomised quality improvement strategy.

INTRODUCTION: The acute respiratory distress syndrome (ARDS) is a highly relevant entity in critical care with mortality rates of 40%. Despite extensive scientific efforts, outcome-relevant therapeutic measures are still insufficiently practised at the bedside. Thus, there is a clear need to adhere to early diagnosis and sufficient therapy in ARDS, assuring lower mortality and multiple organ failure. METHODS AND ANALYSIS: In this quality improvement strategy (QIS), a decision support system as a mobile application (ASIC app), which uses available clinical real-time data, is implemented to support physicians in timely diagnosis and improvement of adherence to established guidelines in the treatment of ARDS. ASIC is conducted on 31 intensive care units (ICUs) at 8 German university hospitals. It is designed as a multicentre stepped-wedge cluster randomised QIS. ICUs are combined into 12 clusters which are randomised in 12 steps. After preparation (18 months) and a control phase of 8 months for all clusters, the first cluster enters a roll-in phase (3 months) that is followed by the actual QIS phase. The remaining clusters follow in month wise steps. The coprimary key performance indicators (KPIs) consist of the ARDS diagnostic rate and guideline adherence regarding lung-protective ventilation. Secondary KPIs include the prevalence of organ dysfunction within 28 days after diagnosis or ICU discharge, the treatment duration on ICU and the hospital mortality. Furthermore, the user acceptance and usability of new technologies in medicine are examined. To show improvements in healthcare of patients with ARDS, differences in primary and secondary KPIs between control phase and QIS will be tested. ETHICS AND DISSEMINATION: Ethical approval was obtained from the independent Ethics Committee (EC) at the RWTH Aachen Faculty of Medicine (local EC reference number: EK 102/19) and the respective data protection officer in March 2019. The results of the ASIC QIS will be presented at conferences and published in peer-reviewed journals. TRIAL REGISTRATION NUMBER: DRKS00014330.

A Model-Based Workflow to Benchmark the Clinical Cholestasis Risk of Drugs.

We present a generic workflow combining physiology-based computational modeling and in vitro data to assess the clinical cholestatic risk of different drugs systematically. Changes in expression levels of genes involved in the enterohepatic circulation of bile acids were obtained from an in vitro assay mimicking 14 days of repeated drug administration for 10 marketed drugs. These changes in gene expression over time were contextualized in a physiology-based bile acid model of glycochenodeoxycholic acid. The simulated drug-induced response in bile acid concentrations was then scaled with the applied drug doses to calculate the cholestatic potential for each compound. A ranking of the cholestatic potential correlated very well with the clinical cholestasis risk obtained from medical literature. The proposed workflow allows benchmarking the cholestatic risk of novel drug candidates. We expect the application of our workflow to significantly contribute to the stratification of the cholestatic potential of new drugs and to support animal-free testing in future drug development.

Data-driven personalization of a physiologically based pharmacokinetic model for caffeine: A systematic assessment.

Physiologically based pharmacokinetic (PBPK) models have been proposed as a tool for more accurate individual pharmacokinetic (PK) predictions and model-informed precision dosing, but their application in clinical practice is still rare. This study systematically assesses the benefit of using individual patient information to improve PK predictions. A PBPK model of caffeine was stepwise personalized by using individual data on (1) demography, (2) physiology, and (3) cytochrome P450 (CYP) 1A2 phenotype of 48 healthy volunteers participating in a single-dose clinical study. Model performance was benchmarked against a caffeine base model simulated with parameters of an average individual. In the first step, virtual twins were generated based on the study subjects' demography (height, weight, age, sex), which implicated the rescaling of average organ volumes and blood flows. The accuracy of PK simulations improved compared with the base model. The percentage of predictions within 0.8-fold to 1.25-fold of the observed values increased from 45.8% (base model) to 57.8% (Step 1). However, setting physiological parameters (liver blood flow determined by magnetic resonance imaging, glomerular filtration rate, hematocrit) to measured values in the second step did not further improve the simulation result (59.1% in the 1.25-fold range). In the third step, virtual twins matching individual demography, physiology, and CYP1A2 activity considerably improved the simulation results. The percentage of data within the 1.25-fold range was 66.15%. This case study shows that individual PK profiles can be predicted more accurately by considering individual attributes and that personalized PBPK models could be a valuable tool for model-informed precision dosing approaches in the future.

Subcellular spatio-temporal intravital kinetics of aflatoxin B1 and ochratoxin A in liver and kidney.

Local accumulation of xenobiotics in human and animal tissues may cause adverse effects. Large differences in their concentrations may exist between individual cell types, often due to the expression of specific uptake and export carriers. Here we established a two-photon microscopy-based technique for spatio-temporal detection of the distribution of mycotoxins in intact kidneys and livers of anesthetized mice with subcellular resolution. The mycotoxins ochratoxin A (OTA, 10 mg/kg b.w.) and aflatoxin B1 (AFB1, 1.5 mg/kg b.w.), which both show blue auto-fluorescence, were analyzed after intravenous bolus injections. Within seconds after administration, OTA was filtered by glomeruli, and enriched in distal tubular epithelial cells (dTEC). A striking feature of AFB1 toxicokinetics was its very rapid uptake from sinusoidal blood into hepatocytes (t1/2 ~ 4 min) and excretion into bile canaliculi. Interestingly, AFB1 was enriched in the nuclei of hepatocytes with zonal differences in clearance. In the cytoplasm of pericentral hepatocytes, the half-life (t1/2~ 63 min) was much longer compared to periportal hepatocytes of the same lobules (t1/2 ~ 9 min). In addition, nuclear AFB1 from periportal hepatocytes cleared faster compared to the pericentral region. These local differences in AFB1 clearance may be due to the pericentral expression of cytochrome P450 enzymes that activate AFB1 to protein- and DNA-binding metabolites. In conclusion, the present study shows that large spatio-temporal concentration differences exist within the same tissues and its analysis may provide valuable additional information to conventional toxicokinetic studies.

Network integration and modelling of dynamic drug responses at multi-omics levels.

Uncovering cellular responses from heterogeneous genomic data is crucial for molecular medicine in particular for drug safety. This can be realized by integrating the molecular activities in networks of interacting proteins. As proof-of-concept we challenge network modeling with time-resolved proteome, transcriptome and methylome measurements in iPSC-derived human 3D cardiac microtissues to elucidate adverse mechanisms of anthracycline cardiotoxicity measured with four different drugs (doxorubicin, epirubicin, idarubicin and daunorubicin). Dynamic molecular analysis at in vivo drug exposure levels reveal a network of 175 disease-associated proteins and identify common modules of anthracycline cardiotoxicity in vitro, related to mitochondrial and sarcomere function as well as remodeling of extracellular matrix. These in vitro-identified modules are transferable and are evaluated with biopsies of cardiomyopathy patients. This to our knowledge most comprehensive study on anthracycline cardiotoxicity demonstrates a reproducible workflow for molecular medicine and serves as a template for detecting adverse drug responses from complex omics data.

A workflow to build PBTK models for novel species.

Physiology-based pharmacokinetic and toxicokinetic (PBPK/TK) models allow us to simulate the concentration of xenobiotica in the plasma and different tissues of an organism. PBPK/TK models are therefore routinely used in many fields of life sciences to simulate the physiological concentration of exogenous compounds in plasma and tissues. The application of PBTK models in ecotoxicology, however, is currently hampered by the limited availability of models for focal species. Here, we present a best practice workflow that describes how to build PBTK models for novel species. To this end, we extrapolated eight previously established rabbit models for several drugs to six additional mammalian species (human, beagle, rat, monkey, mouse, and minipig). We used established PBTK models for these species to account for the species-specific physiology. The parameter sensitivity in the resulting 56 PBTK models was systematically assessed to rank the relevance of the parameters on overall model performance. Interestingly, more than 80% of the 609 considered model parameters showed a negligible sensitivity throughout all models. Only approximately 5% of all parameters had a high sensitivity in at least one of the PBTK models. This approach allowed us to rank the relevance of the various parameters on overall model performance. We used this information to formulate a best practice guideline for the efficient development of PBTK models for novel animal species. We believe that the workflow proposed in this study will significantly support the development of PBTK models for new animal species in the future.

A Physiology-Based Model of Human Bile Acid Metabolism for Predicting Bile Acid Tissue Levels After Drug Administration in Healthy Subjects and BRIC Type 2 Patients.

Drug-induced liver injury (DILI) is a matter of concern in the course of drug development and patient safety, often leading to discontinuation of drug-development programs or early withdrawal of drugs from market. Hepatocellular toxicity or impairment of bile acid (BA) metabolism, known as cholestasis, are the two clinical forms of DILI. Whole-body physiology-based modelling allows a mechanistic investigation of the physiological processes leading to cholestasis in man. Objectives of the present study were: (1) the development of a physiology-based model of the human BA metabolism, (2) population-based model validation and characterisation, and (3) the prediction and quantification of altered BA levels in special genotype subgroups and after drug administration. The developed physiology-based bile acid (PBBA) model describes the systemic BA circulation in humans and includes mechanistically relevant active and passive processes such as the hepatic synthesis, gallbladder emptying, transition through the gastrointestinal tract, reabsorption into the liver, distribution within the whole body, and excretion via urine and faeces. The kinetics of active processes were determined for the exemplary BA glycochenodeoxycholic acid (GCDCA) based on blood plasma concentration-time profiles. The robustness of our PBBA model was verified with population simulations of healthy individuals. In addition to plasma levels, the possibility to estimate BA concentrations in relevant tissues like the intracellular space of the liver enhance the mechanistic understanding of cholestasis. We analysed BA levels in various tissues of Benign Recurrent Intrahepatic Cholestasis type 2 (BRIC2) patients and our simulations suggest a higher susceptibility of BRIC2 patients toward cholestatic DILI due to BA accumulation in the liver. The effect of drugs on systemic BA levels were simulated for cyclosporine A (CsA). Our results confirmed the higher risk of DILI after CsA administration in healthy and BRIC2 patients. The presented PBBA model enhances our mechanistic understanding underlying cholestasis and drug-induced alterations of BA levels in blood and organs. The developed PBBA model might be applied in the future to anticipate potential risk of cholestasis in patients.

Prediction of human drug-induced liver injury (DILI) in relation to oral doses and blood concentrations.

Drug-induced liver injury (DILI) cannot be accurately predicted by animal models. In addition, currently available in vitro methods do not allow for the estimation of hepatotoxic doses or the determination of an acceptable daily intake (ADI). To overcome this limitation, an in vitro/in silico method was established that predicts the risk of human DILI in relation to oral doses and blood concentrations. This method can be used to estimate DILI risk if the maximal blood concentration (Cmax) of the test compound is known. Moreover, an ADI can be estimated even for compounds without information on blood concentrations. To systematically optimize the in vitro system, two novel test performance metrics were introduced, the toxicity separation index (TSI) which quantifies how well a test differentiates between hepatotoxic and non-hepatotoxic compounds, and the toxicity estimation index (TEI) which measures how well hepatotoxic blood concentrations in vivo can be estimated. In vitro test performance was optimized for a training set of 28 compounds, based on TSI and TEI, demonstrating that (1) concentrations where cytotoxicity first becomes evident in vitro (EC10) yielded better metrics than higher toxicity thresholds (EC50); (2) compound incubation for 48 h was better than 24 h, with no further improvement of TSI after 7 days incubation; (3) metrics were moderately improved by adding gene expression to the test battery; (4) evaluation of pharmacokinetic parameters demonstrated that total blood compound concentrations and the 95%-population-based percentile of Cmax were best suited to estimate human toxicity. With a support vector machine-based classifier, using EC10 and Cmax as variables, the cross-validated sensitivity, specificity and accuracy for hepatotoxicity prediction were 100, 88 and 93%, respectively. Concentrations in the culture medium allowed extrapolation to blood concentrations in vivo that are associated with a specific probability of hepatotoxicity and the corresponding oral doses were obtained by reverse modeling. Application of this in vitro/in silico method to the rat hepatotoxicant pulegone resulted in an ADI that was similar to values previously established based on animal experiments. In conclusion, the proposed method links oral doses and blood concentrations of test compounds to the probability of hepatotoxicity.

Bile Microinfarcts in Cholestasis Are Initiated by Rupture of the Apical Hepatocyte Membrane and Cause Shunting of Bile to Sinusoidal Blood.

Bile duct ligation (BDL) is an experimental procedure that mimics obstructive cholestatic disease. One of the early consequences of BDL in rodents is the appearance of so-called bile infarcts that correspond to Charcot-Gombault necrosis in human cholestasis. The mechanisms causing bile infarcts and their pathophysiological relevance are unclear. Therefore, intravital two photon-based imaging of BDL mice was performed with fluorescent bile salts (BS) and non-BS organic anion analogues. Key findings were followed up by matrix-assisted laser desorption ionization imaging, clinical chemistry, immunostaining, and gene expression analyses. In the acute phase, 1-3 days after BDL, BS concentrations in bile increased and single-cell bile microinfarcts occurred in dispersed hepatocytes throughout the liver caused by the rupture of the apical hepatocyte membrane. This rupture occurred after loss of mitochondrial membrane potential, followed by entry of bile, cell death, and a "domino effect" of further death events of neighboring hepatocytes. Bile infarcts provided a trans-epithelial shunt between bile canaliculi and sinusoids by which bile constituents leaked into blood. In the chronic phase, ≥21 days after BDL, uptake of BS tracers at the sinusoidal hepatocyte membrane was reduced. This contributes to elevated concentrations of BS in blood and decreased concentrations in the biliary tract. Conclusion: Bile microinfarcts occur in the acute phase after BDL in a limited number of dispersed hepatocytes followed by larger infarcts involving neighboring hepatocytes, and they allow leakage of bile from the BS-overloaded biliary tract into blood, thereby protecting the liver from BS toxicity; in the chronic phase after BDL, reduced sinusoidal BS uptake is a dominant protective factor, and the kidney contributes to the elimination of BS until cholemic nephropathy sets in.

Using quantitative systems pharmacology to evaluate the drug efficacy of COX-2 and 5-LOX inhibitors in therapeutic situations.

A quantitative analysis of dose-response relationships is essential in preclinical and clinical drug development in order to optimize drug efficacy and safety, respectively. However, there is a lack of quantitative understanding about the dynamics of pharmacological drug-target interactions in biological systems. In this study, a quantitative systems pharmacology (QSP) approach is applied to quantify the drug efficacy of cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX) inhibitors by coupling physiologically based pharmacokinetic models, at the whole-body level, with affected biological networks, at the cellular scale. Both COX-2 and 5-LOX are key enzymes in the production of inflammatory mediators and are known targets in the design of anti-inflammatory drugs. Drug efficacy is here evaluated for single and appropriate co-treatment of diclofenac, celecoxib, zileuton, and licofelone by quantitatively studying the reduction of prostaglandins and leukotrienes. The impact of rifampicin pre-treatment on prostaglandin formation is also investigated by considering pharmacokinetic drug interactions with diclofenac and celecoxib, finally suggesting optimized dose levels to compensate for the reduced drug action. Furthermore, a strong correlation was found between pain relief observed in patients as well as celecoxib- and diclofenac-induced decrease in prostaglandins after 6 h. The findings presented reveal insights about drug-induced modulation of cellular networks in a whole-body context, thereby describing complex pharmacokinetic/pharmacodynamic behavior of COX-2 and 5-LOX inhibitors in therapeutic situations. The results demonstrate the clinical benefit of using QSP to predict drug efficacy and, hence, encourage its use in future drug discovery and development programs.

Spatio-temporal visualization of the distribution of acetaminophen as well as its metabolites and adducts in mouse livers by MALDI MSI.

Acetaminophen (APAP) is one of the most intensively studied compounds that causes hepatotoxicity in the pericentral region of the liver lobules. However, spatio-temporal information on the distribution of APAP, its metabolites and GSH adducts in the liver tissue is not yet available. Here, we addressed the question, whether APAP-GSH adducts and GSH depletion show a zonated pattern and whether the distribution of APAP and its glucuronide as well as sulfate conjugates in liver lobules are zonated. For this purpose, a matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) technique was established, where the MSI images were superimposed onto CYP2E1 immunostained tissue. A time-dependent analysis (5, 15, 30, 60, 120, 240, 480 min) after intraperitoneal administration of 300 mg/kg APAP and a dose-dependent analysis (56 up to 500 mg APAP/kg) at 30 min were performed. The results demonstrate that the MALDI MSI technique allows the assignment of compounds and their metabolites to specific lobular zones. APAP-GSH adducts and GSH depletion occurred predominantly in the CYP2E1-positive zone of the liver, although GSH also decreased in the periportal region. In contrast, the parent compound, APAP sulfate and APAP glucuronide did not show a zonated pattern and tissue concentrations showed a similar time course as the corresponding analyses were performed with blood from the portal and liver veins. In conclusion, the present study is in agreement with the concept that pericentral CYPs form NAPQI that in the same cell binds to and depletes GSH but a lower level of GSH adducts is also observed in the periportal region. The results also provide further evidence of the recently published concept of 'aggravated loss of clearance capacity' according to which also liver tissue that survives intoxication may transiently show decreased metabolic capacity.

Bringing in vitro analysis closer to in vivo: Studying doxorubicin toxicity and associated mechanisms in 3D human microtissues with PBPK-based dose modelling.

Doxorubicin (DOX) is a chemotherapeutic agent of which the medical use is limited due to cardiotoxicity. While acute cardiotoxicity is reversible, chronic cardiotoxicity is persistent or progressive, dose-dependent and irreversible. While DOX mechanisms of action are not fully understood yet, 3 toxicity processes are known to occur in vivo: cardiomyocyte dysfunction, mitochondrial dysfunction and cell death. We present an in vitro experimental design aimed at detecting DOX-induced cardiotoxicity by obtaining a global view of the induced molecular mechanisms through RNA-sequencing. To better reflect the in vivo situation, human 3D cardiac microtissues were exposed to physiologically-based pharmacokinetic (PBPK) relevant doses of DOX for 2 weeks. We analysed a therapeutic and a toxic dosing profile. Transcriptomics analysis revealed significant gene expression changes in pathways related to "striated muscle contraction" and "respiratory electron transport", thus suggesting mitochondrial dysfunction as an underlying mechanism for cardiotoxicity. Furthermore, expression changes in mitochondrial processes differed significantly between the doses. Therapeutic dose reflects processes resembling the phenotype of delayed chronic cardiotoxicity, while toxic doses resembled acute cardiotoxicity. Overall, these results demonstrate the capability of our innovative in vitro approach to detect the three known mechanisms of DOX leading to toxicity, thus suggesting its potential relevance for reflecting the patient situation. Our study also demonstrated the importance of applying physiologically relevant doses during toxicological research, since mechanisms of acute and chronic toxicity differ.