Identifying multiscale translational safety biomarkers using a network-based systems approach.

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

The human hepatocyte TXG-MAPr: gene co-expression network modules to support mechanism-based risk assessment.

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

Contextualization of causal regulatory networks from toxicogenomics data applied to drug-induced liver injury.

In recent years, network-based methods have become an attractive analytical approach for toxicogenomics studies. They can capture not only the global changes of regulatory gene networks but also the relationships between their components. Among them, a causal reasoning approach depicts the mechanisms of regulation that connect upstream regulators in signaling networks to their downstream gene targets. In this work, we applied CARNIVAL, a causal network contextualisation tool, to infer upstream signaling networks deregulated in drug-induced liver injury (DILI) from gene expression microarray data from the TG-GATEs database. We focussed on six compounds that induce observable histopathologies linked to DILI from repeated dosing experiments in rats. We compared responses in vitro and in vivo to identify potential cross-platform concordances in rats as well as network preservations between rat and human. Our results showed similarities of enriched pathways and network motifs between compounds. These pathways and motifs induced the same pathology in rats but not in humans. In particular, the causal interactions "LCK activates SOCS3, which in turn inhibits TFDP1" was commonly identified as a regulatory path among the fibrosis-inducing compounds. This potential pathology-inducing regulation illustrates the value of our approach to generate hypotheses that can be further validated experimentally.

Integration of temporal single cell cellular stress response activity with logic-ODE modeling reveals activation of ATF4-CHOP axis as a critical predictor of drug-induced liver injury.

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

Converting networks to predictive logic models from perturbation signalling data with CellNOpt.

SUMMARY: The molecular changes induced by perturbations such as drugs and ligands are highly informative of the intracellular wiring. Our capacity to generate large datasets is increasing steadily. A useful way to extract mechanistic insight from the data is by integrating them with a prior knowledge network of signalling to obtain dynamic models. CellNOpt is a collection of Bioconductor R packages for building logic models from perturbation data and prior knowledge of signalling networks. We have recently developed new components and refined the existing ones to keep up with the computational demand of increasingly large datasets, including (i) an efficient integer linear programming, (ii) a probabilistic logic implementation for semi-quantitative datasets, (iii) the integration of a stochastic Boolean simulator, (iv) a tool to identify missing links, (v) systematic post-hoc analyses and (vi) an R-Shiny tool to run CellNOpt interactively. AVAILABILITY AND IMPLEMENTATION: R-package(s): https://github.com/saezlab/cellnopt. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Quantitative Systems Toxicology Modeling To Address Key Safety Questions in Drug Development: A Focus of the TransQST Consortium.

Omics data have been increasingly generated with limited demonstrated value in drug safety assessment. The TransQST consortium was launched to use omics and other data in mechanistic-based quantitative systems toxicology (QST) models to evaluate their potential use in species translation.

Dysregulated mesenchymal PDGFR-ß drives kidney fibrosis.

Kidney fibrosis is characterized by expansion and activation of platelet-derived growth factor receptor-ß (PDGFR-ß)-positive mesenchymal cells. To study the consequences of PDGFR-ß activation, we developed a model of primary renal fibrosis using transgenic mice with PDGFR-ß activation specifically in renal mesenchymal cells, driving their pathological proliferation and phenotypic switch toward myofibroblasts. This resulted in progressive mesangioproliferative glomerulonephritis, mesangial sclerosis, and interstitial fibrosis with progressive anemia due to loss of erythropoietin production by fibroblasts. Fibrosis induced secondary tubular epithelial injury at later stages, coinciding with microinflammation, and aggravated the progression of hypertensive and obstructive nephropathy. Inhibition of PDGFR activation reversed fibrosis more effectively in the tubulointerstitium compared to glomeruli. Gene expression signatures in mice with PDGFR-ß activation resembled those found in patients. In conclusion, PDGFR-ß activation alone is sufficient to induce progressive renal fibrosis and failure, mimicking key aspects of chronic kidney disease in humans. Our data provide direct proof that fibrosis per se can drive chronic organ damage and establish a model of primary fibrosis allowing specific studies targeting fibrosis progression and regression.

From expression footprints to causal pathways: contextualizing large signaling networks with CARNIVAL.

While gene expression profiling is commonly used to gain an overview of cellular processes, the identification of upstream processes that drive expression changes remains a challenge. To address this issue, we introduce CARNIVAL, a causal network contextualization tool which derives network architectures from gene expression footprints. CARNIVAL (CAusal Reasoning pipeline for Network identification using Integer VALue programming) integrates different sources of prior knowledge including signed and directed protein-protein interactions, transcription factor targets, and pathway signatures. The use of prior knowledge in CARNIVAL enables capturing a broad set of upstream cellular processes and regulators, leading to a higher accuracy when benchmarked against related tools. Implementation as an integer linear programming (ILP) problem guarantees efficient computation. As a case study, we applied CARNIVAL to contextualize signaling networks from gene expression data in IgA nephropathy (IgAN), a condition that can lead to chronic kidney disease. CARNIVAL identified specific signaling pathways and associated mediators dysregulated in IgAN including Wnt and TGF-ß, which we subsequently validated experimentally. These results demonstrated how CARNIVAL generates hypotheses on potential upstream alterations that propagate through signaling networks, providing insights into diseases.

Embracing the Dark Side: Computational Approaches to Unveil the Functionality of Genes Lacking Biological Annotation in Drug-Induced Liver Injury.

In toxicogenomics, functional annotation is an important step to gain additional insights into genes with aberrant expression that drive pathophysiological mechanisms. Nevertheless, there exists a gap on annotation of these genes which often hampers the interpretation of results and limits their applicability in translational medicine. In this study, we evaluated the coverage of functional annotations of differentially expressed genes (DEGs) induced by 10 selected compounds from the TG-GATEs database identified as high- or no-risk in causing drug-induced liver injury (most-DILI or no-DILI, respectively) using in vitro human data. Functional roles of DEGs not present in the most common biological annotation databases - termed "dark genes" - were unveiled via literature mining and via the identification of shared regulatory transcription factors or signaling pathways. Our results demonstrated that there were approximately 13% of dark genes induced by these compounds in vitro and we were able to obtain additional relevant information for up to 76% of those. Using interactome data from several sources, we have uncovered genes such as LRBA, and WDR26 as highly connected in the protein network that play roles in drug response. Genes such as MALAT1, H19, and MIR29C - whose links to hepatotoxicity have been confirmed - were identified as markers for the most-DILI group and appeared as top hits across all literature-based mining methods. Furthermore, we investigated the potential impact of dark genes on liver toxicity by identifying their rat orthologs in combination with their correlation to drug-induced liver pathologies observed in vivo following chemical exposure. We identified a set of important regulatory transcription factors of dark genes for all most-DILI compounds including E2F1 and JUND with supporting evidences in literature and we found Magee1 correlated with chemically induced bile duct hyperplasia and adverse responses at 29 days in rats in vivo. In conclusion, in this study we show the potential role of these poorly annotated genes in mechanisms underlying hepatotoxicity and offer a number of computational approaches that may help to minimize current gaps in gene annotation and highlight their values as potential biomarkers in toxicological studies.

FALCON: a toolbox for the fast contextualization of logical networks.

MOTIVATION: Mathematical modelling of regulatory networks allows for the discovery of knowledge at the system level. However, existing modelling tools are often computation-heavy and do not offer intuitive ways to explore the model, to test hypotheses or to interpret the results biologically. RESULTS: We have developed a computational approach to contextualize logical models of regulatory networks with biological measurements based on a probabilistic description of rule-based interactions between the different molecules. Here, we propose a Matlab toolbox, FALCON, to automatically and efficiently build and contextualize networks, which includes a pipeline for conducting parameter analysis, knockouts and easy and fast model investigation. The contextualized models could then provide qualitative and quantitative information about the network and suggest hypotheses about biological processes. AVAILABILITY AND IMPLEMENTATION: FALCON is freely available for non-commercial users on GitHub under the GPLv3 licence. The toolbox, installation instructions, full documentation and test datasets are available at https://github.com/sysbiolux/FALCON. FALCON runs under Matlab (MathWorks) and requires the Optimization Toolbox. CONTACT: thomas.sauter@uni.lu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A Probabilistic Boolean Network Approach for the Analysis of Cancer-Specific Signalling: A Case Study of Deregulated PDGF Signalling in GIST.

BACKGROUND: Signal transduction networks are increasingly studied with mathematical modelling approaches while each of them is suited for a particular problem. For the contextualisation and analysis of signalling networks with steady-state protein data, we identified probabilistic Boolean network (PBN) as a promising framework which could capture quantitative changes of molecular changes at steady-state with a minimal parameterisation. RESULTS AND CONCLUSION: In our case study, we successfully applied the PBN approach to model and analyse the deregulated Platelet-Derived Growth Factor (PDGF) signalling pathway in Gastrointestinal Stromal Tumour (GIST). We experimentally determined a rich and accurate dataset of steady-state profiles of selected downstream kinases of PDGF-receptor-alpha mutants in combination with inhibitor treatments. Applying the tool optPBN, we fitted a literature-derived candidate network model to the training dataset consisting of single perturbation conditions. Model analysis suggested several important crosstalk interactions. The validity of these predictions was further investigated experimentally pointing to relevant ongoing crosstalk from PI3K to MAPK signalling in tumour cells. The refined model was evaluated with a validation dataset comprising multiple perturbation conditions. The model thereby showed excellent performance allowing to quantitatively predict the combinatorial responses from the individual treatment results in this cancer setting. The established optPBN pipeline is also widely applicable to gain a better understanding of other signalling networks at steady-state in a context-specific fashion.

L-plastin Ser5 phosphorylation in breast cancer cells and in vitro is mediated by RSK downstream of the ERK/MAPK pathway.

Deregulated cell migration and invasion are hallmarks of metastatic cancer cells. Phosphorylation on residue Ser5 of the actin-bundling protein L-plastin activates L-plastin and has been reported to be crucial for invasion and metastasis. Here, we investigate signal transduction leading to L-plastin Ser5 phosphorylation using 4 human breast cancer cell lines. Whole-genome microarray analysis comparing cell lines with different invasive capacities and corresponding variations in L-plastin Ser5 phosphorylation level revealed that genes of the ERK/MAPK pathway are differentially expressed. It is noteworthy that in vitro kinase assays showed that ERK/MAPK pathway downstream ribosomal protein S6 kinases α-1 (RSK1) and α-3 (RSK2) are able to directly phosphorylate L-plastin on Ser5. Small interfering RNA- or short hairpin RNA-mediated knockdown and activation/inhibition studies followed by immunoblot analysis and computational modeling confirmed that ribosomal S6 kinase (RSK) is an essential activator of L-plastin. Migration and invasion assays showed that RSK knockdown led to a decrease of up to 30% of migration and invasion of MDA-MB-435S cells. Although the presence of L-plastin was not necessary for migration/invasion of these cells, immunofluorescence assays illustrated RSK-dependent recruitment of Ser5-phosphorylated L-plastin to migratory structures. Altogether, we provide evidence that the ERK/MAPK pathway is involved in L-plastin Ser5 phosphorylation in breast cancer cells with RSK1 and RSK2 kinases able to directly phosphorylate L-plastin residue Ser5.

optPBN: an optimisation toolbox for probabilistic Boolean networks.

BACKGROUND: There exist several computational tools which allow for the optimisation and inference of biological networks using a Boolean formalism. Nevertheless, the results from such tools yield only limited quantitative insights into the complexity of biological systems because of the inherited qualitative nature of Boolean networks. RESULTS: We introduce optPBN, a Matlab-based toolbox for the optimisation of probabilistic Boolean networks (PBN) which operates under the framework of the BN/PBN toolbox. optPBN offers an easy generation of probabilistic Boolean networks from rule-based Boolean model specification and it allows for flexible measurement data integration from multiple experiments. Subsequently, optPBN generates integrated optimisation problems which can be solved by various optimisers. In term of functionalities, optPBN allows for the construction of a probabilistic Boolean network from a given set of potential constitutive Boolean networks by optimising the selection probabilities for these networks so that the resulting PBN fits experimental data. Furthermore, the optPBN pipeline can also be operated on large-scale computational platforms to solve complex optimisation problems. Apart from exemplary case studies which we correctly inferred the original network, we also successfully applied optPBN to study a large-scale Boolean model of apoptosis where it allows identifying the inverse correlation between UVB irradiation, NFκB and Caspase 3 activations, and apoptosis in primary hepatocytes quantitatively. Also, the results from optPBN help elucidating the relevancy of crosstalk interactions in the apoptotic network. SUMMARY: The optPBN toolbox provides a simple yet comprehensive pipeline for integrated optimisation problem generation in the PBN formalism that can readily be solved by various optimisers on local or grid-based computational platforms. optPBN can be further applied to various biological studies such as the inference of gene regulatory networks or the identification of the interaction's relevancy in signal transduction networks.

Recent development and biomedical applications of probabilistic Boolean networks.

Probabilistic Boolean network (PBN) modelling is a semi-quantitative approach widely used for the study of the topology and dynamic aspects of biological systems. The combined use of rule-based representation and probability makes PBN appealing for large-scale modelling of biological networks where degrees of uncertainty need to be considered.A considerable expansion of our knowledge in the field of theoretical research on PBN can be observed over the past few years, with a focus on network inference, network intervention and control. With respect to areas of applications, PBN is mainly used for the study of gene regulatory networks though with an increasing emergence in signal transduction, metabolic, and also physiological networks. At the same time, a number of computational tools, facilitating the modelling and analysis of PBNs, are continuously developed.A concise yet comprehensive review of the state-of-the-art on PBN modelling is offered in this article, including a comparative discussion on PBN versus similar models with respect to concepts and biomedical applications. Due to their many advantages, we consider PBN to stand as a suitable modelling framework for the description and analysis of complex biological systems, ranging from molecular to physiological levels.