BACKGROUND\nOBJECTIVES\nMETHODS\nRESULTS\nDISCUSSION\nUnderstanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as... Show moreBACKGROUND\nOBJECTIVES\nMETHODS\nRESULTS\nDISCUSSION\nUnderstanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses.\nWe aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors.\n signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli.\n signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3.\nOverall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro-based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891. Show less
Braak, B. ter; Niemeijer, M.; Boon, R.; Parmentier, C.; Baze , A.; Richert. L.; ... ; Water, B. van de 2021
Various adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte... Show moreVarious adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte-like cells (HLCs) provide a promising tool to study cellular stress response pathways, but in this context there is limited insight on how HLCs compare to other in vitro liver models. Here, we systematically compared the transcriptomic profiles upon chemical activation in HLCs, hiPSC, primary human hepatocytes (PHH) and HepG2 liver cancer cells. We used targeted RNA-sequencing to map concentration transcriptional response using benchmark concentration modeling for the various stress responses in the different test systems. We found that HLCs are very sensitive towards oxidative stress and inflammation conditions as corresponding genes were activated at over 3 fold lower concentrations of the corresponding pathway inducing compounds as compared to PHH. PHH were the most sensitive model when studying UPR related effects. Due to the non-proliferative nature of PHH and HLCs, these do not pose a good/sensitive model to pick up DNA damage responses, while hiPSC and HepG2 were more sensitive in these conditions. We envision that this study contributes to a better understanding on how HLCs can contribute to the assessment of cell physiological stress response activation to predict hepatotoxic events. Show less
The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the... Show moreThe TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the "measured only" and "whole transcriptome" (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the "measured only" genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original "measured only" dataset. Furthermore, the inclusion of the extrapolated genes raised "tunicamycin" from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations. Show less
Idiosyncratic drug-induced liver injury (IDILI) is a severe disease that cannot be detected during drug development. It has been shown that hepatotoxicity of some compounds associated with IDILI... Show moreIdiosyncratic drug-induced liver injury (IDILI) is a severe disease that cannot be detected during drug development. It has been shown that hepatotoxicity of some compounds associated with IDILI becomes apparent when these are combined in vivo and in vitro with LPS or TNF. Among these compounds trovafloxacin (TVX) induced apoptosis in the liver and increased pro-inflammatory cytokines in mice exposed to LPS/TNF. The hepatocyte survival and the cytokine release after TNF/LPS stimulation relies on a pulsatile activation of NF-κB. We set out to evaluate the dynamic activation of NF-κB in response to TVX + TNF or LPS models, both in mouse and human cells. Remarkably, TVX prolonged the first translocation of NF-κB induced by TNF both in vivo and in vitro. The prolonged p65 translocation caused by TVX was associated with an increased phosphorylation of IKK and MAPKs and accumulation of inhibitors of NF-κB such as IκBα and A20 in HepG2. Coherently, TVX suppressed further TNF-induced NF-κB translocations in HepG2 leading to decreased transcription of ICAM-1 and inhibitors of apoptosis. TVX prolonged LPS-induced NF-κB translocation in RAW264.7 macrophages increasing the secretion of TNF. In summary, this study presents new, relevant insights into the mechanism of TVX-induced liver injury underlining the resemblance between mouse and human models. In this study we convincingly show that regularly used toxicity models provide a coherent view of relevant pathways for IDILI. We propose that assessment of the kinetics of activation of NF-κB and MAPKs is an appropriate tool for the identification of hepatotoxic compounds during drug development. Show less
Oppelt, A.; Kaschek, D.; Huppelschoten, S.; Sison-Young, R.; Zhang, F.; Buck-Wiese, M.; ... ; Klingmüller, U. 2018
Drug-induced liver injury (DILI) has become a major problem for patients and for clinicians, academics and the pharmaceutical industry. To date, existing hepatotoxicity test systems are only poorly... Show moreDrug-induced liver injury (DILI) has become a major problem for patients and for clinicians, academics and the pharmaceutical industry. To date, existing hepatotoxicity test systems are only poorly predictive and the underlying mechanisms are still unclear. One of the factors known to amplify hepatotoxicity is the tumor necrosis factor alpha (TNFα), especially due to its synergy with commonly used drugs such as diclofenac. However, the exact mechanism of how diclofenac in combination with TNFα induces liver injury remains elusive. Here, we combined time-resolved immunoblotting and live-cell imaging data of HepG2 cells and primary human hepatocytes (PHH) with dynamic pathway modeling using ordinary differential equations (ODEs) to describe the complex structure of TNFα-induced NFκB signal transduction and integrated the perturbations of the pathway caused by diclofenac. The resulting mathematical model was used to systematically identify parameters affected by diclofenac. These analyses showed that more than one regulatory module of TNFα-induced NFκB signal transduction is affected by diclofenac, suggesting that hepatotoxicity is the integrated consequence of multiple changes in hepatocytes and that multiple factors define toxicity thresholds. Applying our mathematical modeling approach to other DILI-causing compounds representing different putative DILI mechanism classes enabled us to quantify their impact on pathway activation, highlighting the potential of the dynamic pathway model as a quantitative tool for the analysis of DILI compounds. Show less
Wink, S.; Hiemstra, S.W.; Huppelschoten, S.; Klip, J.E.; Water, B. van de 2018
Drug-induced liver injury remains a concern during drug treatment and development. There is an urgent need for improved mechanistic understanding and prediction of DILI liabilities using in vitro... Show moreDrug-induced liver injury remains a concern during drug treatment and development. There is an urgent need for improved mechanistic understanding and prediction of DILI liabilities using in vitro approaches. We have established and characterized a panel of liver cell models containing mechanism-based fluorescent protein toxicity pathway reporters to quantitatively assess the dynamics of cellular stress response pathway activation at the single cell level using automated live cell imaging. We have systematically evaluated the application of four key adaptive stress pathway reporters for the prediction of DILI liability: SRXN1-GFP (oxidative stress), CHOP-GFP (ER stress/UPR response), p21 (p53-mediated DNA damage-related response) and ICAM1 (NF-κB-mediated inflammatory signaling). 118 FDA-labeled drugs in five human exposure relevant concentrations were evaluated for reporter activation using live cell confocal imaging. Quantitative data analysis revealed activation of single or multiple reporters by most drugs in a concentration and time dependent manner. Hierarchical clustering of time course dynamics and refined single cell analysis allowed the allusion of key events in DILI liability. Concentration response modeling was performed to calculate benchmark concentrations (BMCs). Extracted temporal dynamic parameters and BMCs were used to assess the predictive power of sub-lethal adaptive stress pathway activation. Although cellular adaptive responses were activated by non-DILI and severe-DILI compounds alike, dynamic behavior and lower BMCs of pathway activation were sufficiently distinct between these compound classes. The high-level detailed temporal- and concentration-dependent evaluation of the dynamics of adaptive stress pathway activation adds to the overall understanding and prediction of drug-induced liver liabilities. Show less
In previous studies at our laboratory it was demonstrated that drug exposure of HepG2 cells can lead to an altered TNFα-induced NF-κB oscillatory phenotype, concurrent with a synergistically... Show moreIn previous studies at our laboratory it was demonstrated that drug exposure of HepG2 cells can lead to an altered TNFα-induced NF-κB oscillatory phenotype, concurrent with a synergistically increased sensitivity for TNFα-induced apoptosis. We have also shown that synergistic drug/ TNFα-induced cell death is dependent on concurrent cellular stress responses. To monitor these stress responses, we have developed a fluorescent protein stress response reporter platform, feasible for high throughput approaches. In this thesis, we focus on an in vitro HepG2 cell model, in which the addition of TNFα reflects inflammatory stress. We studied inflammatory signaling, both for use in improved in vitro testing approaches and to gain mechanistic understanding of TNFα-induced signaling in hepatocytes. Show less
Sluis, T.C. van der; Duikeren, S. van; Huppelschoten, S.; Jordanova, E.S.; Nejad, E.B.; Sloots, A.; ... ; Melief, C.J.M. 2015