Drug-induced liver injury (DILI) remains the main reason for drug development attritions largely due to poor mechanistic understanding. Toxicogenomic to interrogate the mechanism of DILI has been... Show moreDrug-induced liver injury (DILI) remains the main reason for drug development attritions largely due to poor mechanistic understanding. Toxicogenomic to interrogate the mechanism of DILI has been broadly performed. Gene co-regulation network-based transcriptome analysis is a bioinformatics approach that potentially contributes to improve mechanistic interpretation of toxicogenomic data. Here we performed an extensive concentration time course response-toxicogenomic study in the HepG2 cell line exposed to 20 DILI compounds, 7 reference compounds for stress response pathways, and 10 agonists for cytokines and growth factor receptors. We performed whole transcriptome targeted RNA sequencing to more than 500 conditions to and applied weighted gene co-regulated network analysis (WGCNA) to the transcriptomics data followed by identification of gene co-regulated networks (modules) that were strongly modulated upon the exposure of DILI compounds. Preservation analysis on the module responses of HepG2 and PHH demonstrated highly preserved adaptive stress response gene co-regulated networks. We correlated gene co-regulated networks with cell death onset and causal relationships of 67 critical target genes of these modules with onset of cell death was evaluated using RNA interference screening. We identified GTPBP2, HSPA1B, IRF1, SIRT1 and TSC22D3 as essential modulators of DILI compound-induced cell death. These genes were also induced by DILI compounds in PHH. Altogether, we demonstrate the application of large transcriptome datasets combined with network-based analysis and biological validation to uncover the candidate determinants of DILI. Show less
Callegaro, G.; Schimming, J.P.; Piñero González, J.; Kunnen, S.J.; Wijaya, L.S.; Trairatphisan, P.; ... ; Water, B. van de 2023
Animal testing is the current standard for drug and chemicals safety assessment, but hazards translation to human is uncertain. Human in vitro models can address the species translation but might... Show moreAnimal testing is the current standard for drug and chemicals safety assessment, but hazards translation to human is uncertain. Human in vitro models can address the species translation but might not replicate in vivo complexity. Herein, we propose a network-based method addressing these translational multiscale problems that derives in vivo liver injury biomarkers applicable to in vitro human early safety screening. We applied weighted correlation network analysis (WGCNA) to a large rat liver transcriptomic dataset to obtain co-regulated gene clusters (modules). We identified modules statistically associated with liver pathologies, including a module enriched for ATF4-regulated genes as associated with the occurrence of hepatocellular single-cell necrosis, and as preserved in human liver in vitro models. Within the module, we identified TRIB3 and MTHFD2 as a novel candidate stress biomarkers, and developed and used BAC-eGFPHepG2 reporters in a compound screening, identifying compounds showing ATF4-dependent stress response and potential early safety signals. Show less
Callegaro, G.; Kunnen, S.J.; Trairatphisan, P.; Grosdidier, S.; Niemeijer, M.; Hollander, W. den; ... ; Water, B. van de 2021
Mechanism-based risk assessment is urged to advance and fully permeate into current safety assessment practices, possibly at early phases of drug safety testing. Toxicogenomics is a promising... Show moreMechanism-based risk assessment is urged to advance and fully permeate into current safety assessment practices, possibly at early phases of drug safety testing. Toxicogenomics is a promising source of mechanisms-revealing data, but interpretative analysis tools specific for the testing systems (e.g. hepatocytes) are lacking. In this study, we present the TXG-MAPr webtool (available at https://txg-mapr.eu/WGCNA_PHH/TGGATEs_PHH/ ), an R-Shiny-based implementation of weighted gene co-expression network analysis (WGCNA) obtained from the Primary Human Hepatocytes (PHH) TG-GATEs dataset. The 398 gene co-expression networks (modules) were annotated with functional information (pathway enrichment, transcription factor) to reveal their mechanistic interpretation. Several well-known stress response pathways were captured in the modules, were perturbed by specific stressors and showed preservation in rat systems (rat primary hepatocytes and rat in vivo liver), with the exception of DNA damage and oxidative stress responses. A subset of 87 well-annotated and preserved modules was used to evaluate mechanisms of toxicity of endoplasmic reticulum (ER) stress and oxidative stress inducers, including cyclosporine A, tunicamycin and acetaminophen. In addition, module responses can be calculated from external datasets obtained with different hepatocyte cells and platforms, including targeted RNA-seq data, therefore, imputing biological responses from a limited gene set. As another application, donors' sensitivity towards tunicamycin was investigated with the TXG-MAPr, identifying higher basal level of intrinsic immune response in donors with pre-existing liver pathology. In conclusion, we demonstrated that gene co-expression analysis coupled to an interactive visualization environment, the TXG-MAPr, is a promising approach to achieve mechanistic relevant, cross-species and cross-platform evaluation of toxicogenomic data. Show less
Wijaya, L.S.; Trairatphisan, P.; Gabor, A.; Niemeijer, M.C.; Keet, J.; Alcalà Morera, A.; ... ; Water, B. van de 2021
Drug-induced liver injury (DILI) is the most prevalent adversity encountered in drug development and clinical settings leading to urgent needs to understand the underlying mechanisms. In this study... Show moreDrug-induced liver injury (DILI) is the most prevalent adversity encountered in drug development and clinical settings leading to urgent needs to understand the underlying mechanisms. In this study, we have systematically investigated the dynamics of the activation of cellular stress response pathways and cell death outcomes upon exposure of a panel of liver toxicants using live cell imaging of fluorescent reporter cell lines. We established a comprehensive temporal dynamic response profile of a large set of BAC-GFP HepG2 cell lines representing the following components of stress signaling: i) unfolded protein response (UPR) [ATF4, XBP1, BIP and CHOP]; ii) oxidative stress [NRF2, SRXN1, HMOX1]; iii) DNA damage [P53, P21, BTG2, MDM2]; and iv) NF-κB pathway [A20, ICAM1]. We quantified the single cell GFP expression as a surrogate for endogenous protein expression using live cell imaging over > 60 h upon exposure to 14 DILI compounds at multiple concentrations. Using logic-based ordinary differential equation (Logic-ODE), we modelled the dynamic profiles of the different stress responses and extracted specific descriptors potentially predicting the progressive outcomes. We identified the activation of ATF4-CHOP axis of the UPR as the key pathway showing the highest correlation with cell death upon DILI compound perturbation. Knocking down main components of the UPR provided partial protection from compound-induced cytotoxicity, indicating a complex interplay among UPR components as well as other stress pathways. Our results suggest that a systematic analysis of the temporal dynamics of ATF4-CHOP axis activation can support the identification of DILI risk for new candidate drugs. Show less
Copple, I.M.; Hollander, W. den; Callegaro, G.; Mutter, F.E.; Maggs, J.L.; Schofield, A.L.; ... ; Park, B.K. 2018
The transcription factor NRF2, governed by its repressor KEAP1, protects cells against oxidative stress. There is interest in modelling the NRF2 response to improve the prediction of clinical... Show moreThe transcription factor NRF2, governed by its repressor KEAP1, protects cells against oxidative stress. There is interest in modelling the NRF2 response to improve the prediction of clinical toxicities such as drug-induced liver injury (DILI). However, very little is known about the makeup of the NRF2 transcriptional network and its response to chemical perturbation in primary human hepatocytes (PHH), which are often used as a translational model for investigating DILI. Here, microarray analysis identified 108 transcripts (including several putative novel NRF2-regulated genes) that were both downregulated by siRNA targeting NRF2 and upregulated by siRNA targeting KEAP1 in PHH. Applying weighted gene co-expression network analysis (WGCNA) to transcriptomic data from the Open TG-GATES toxicogenomics repository (representing PHH exposed to 158 compounds) revealed four co-expressed gene sets or 'modules' enriched for these and other NRF2-associated genes. By classifying the 158 TG-GATES compounds based on published evidence, and employing the four modules as network perturbation metrics, we found that the activation of NRF2 is a very good indicator of the intrinsic biochemical reactivity of a compound (i.e. its propensity to cause direct chemical stress), with relatively high sensitivity, specificity, accuracy and positive/negative predictive values. We also found that NRF2 activation has lower sensitivity for the prediction of clinical DILI risk, although relatively high specificity and positive predictive values indicate that false positive detection rates are likely to be low in this setting. Underpinned by our comprehensive analysis, activation of the NRF2 network is one of several mechanism-based components that can be incorporated into holistic systems toxicology models to improve mechanistic understanding and preclinical prediction of DILI in man. Show less
Water, B. van de; Tijdens, R.B.; Verbrugge, A.; Huigsloot, M.; Dihal, A.A.; Stevens, J.L.; ... ; Mulder, G.J. 2000
Decreased phosphorylation of focal adhesion kinase and paxillin is associated with loss of focal adhesions and stress fibers and precedes the onset of apoptosis (van de Water, B., Nagelkerke, J. F... Show moreDecreased phosphorylation of focal adhesion kinase and paxillin is associated with loss of focal adhesions and stress fibers and precedes the onset of apoptosis (van de Water, B., Nagelkerke, J. F., and Stevens, J. L. (1999) J. Biol. Chem. 274, 13328-13337). The cortical actin cytoskeletal network is also lost during apoptosis, yet little is known about the temporal relationship between altered phosphorylation of proteins that are critical in the regulation of this network and their potential cleavage by caspases during apoptosis. Adducins are central in the cortical actin network organization. Cisplatin caused apoptosis of renal proximal tubular epithelial cells, which was associated with the cleavage of alpha-adducin into a 74-kDa fragment; this was blocked by a general caspase inhibitor benzyloxycarbonyl-Val-Ala-Asp-fluoromethyl ketone (z-VAD-fmk). Hemagglutinin-tagged human alpha-adducin was cleaved into a similar 74-kDa fragment by caspase-3 in vitro but not by caspase-6 or -7. Asp-Arg-Val-Asp(29)-Glu, Asp-Ile-Val-Asp(208)-Arg, and Asp-Asp-Ser-Asp(633)-Ala were identified as the principal caspase-3 cleavage sites; Asp-Asp-Ser-Asp(633)-Ala was key in the formation of the 74-kDa fragment. Cisplatin also caused an increased phosphorylation of alpha-adducin and gamma-adducin in the MARCKS domain that preceded alpha-adducin cleavage and was associated with loss of adducins from adherens junctions; this was not affected by z-VAD-fmk. In conclusion, the data support a model in which increased phosphorylation of alpha-adducin due to cisplatin leads to dissociation from the cytoskeleton, a situation rendered irreversible by caspase-3-mediated cleavage of alpha-adducin at Asp-Asp-Ser-Asp(633)-Ala. Show less
Water, B. van de; Nagelkerke, J.F.; Stevens, J.L. 1999
The relationship between focal adhesion protein (FAK) activity and loss of cell-matrix contact during apoptosis is not entirely clear nor has the role of FAK in chemically induced apoptosis been... Show moreThe relationship between focal adhesion protein (FAK) activity and loss of cell-matrix contact during apoptosis is not entirely clear nor has the role of FAK in chemically induced apoptosis been studied. We investigated the status of FAK phosphorylation and cleavage in renal epithelial cells during apoptosis caused by the nephrotoxicant dichlorovinylcysteine (DCVC). DCVC treatment caused a loss of cell-matrix contact which was preceded by a dissociation of FAK from the focal adhesions and tyrosine dephosphorylation of FAK. Paxillin was also dephosphorylated at tyrosine. DCVC treatment activated caspase-3 which was associated with cleavage of FAK. However, FAK cleavage occurred after cells had already lost focal adhesions indicating that cleavage of FAK by caspases is not responsible for loss of FAK from focal adhesions. Accordingly, although inhibition of caspase activity with zVAD-fmk blocked activation of caspase-3, FAK cleavage, and apoptosis, it neither affected dephosphorylation nor translocation of FAK or paxillin. However, zVAD-fmk completely blocked the cell detachment caused by DCVC treatment. Orthovanadate prevented DCVC-induced tyrosine dephosphorylation of both FAK and paxillin; however, it did not inhibit DCVC-induced apoptosis and actually potentiated focal adhesion disorganization and cell detachment. Thus, FAK dephosphorylation and loss of focal adhesions are not due to caspase activation; however, caspases are required for FAK proteolysis and cell detachment. Show less
Liu, H.; Miller, E.; Water, B. van de; Stevens, J.L. 1998