This review encompasses the most important advances in liver functions and hepatotoxicity and analyzes which mechanisms can be studied in vitro. In a complex architecture of nested, zonated lobules, the liver consists of approximately 80 % hepatocytes and 20 % non-parenchymal cells, the latter being involved in a secondary phase that may dramatically aggravate the initial damage. Hepatotoxicity, as well as hepatic metabolism, is controlled by a set of nuclear receptors (including PXR, CAR, HNF-4α, FXR, LXR, SHP, VDR and PPAR) and signaling pathways. When isolating liver cells, some pathways are activated, e.g., the RAS/MEK/ERK pathway, whereas others are silenced (e.g. HNF-4α), resulting in up- and downregulation of hundreds of genes. An understanding of these changes is crucial for a correct interpretation of in vitro data. The possibilities and limitations of the most useful liver in vitro systems are summarized, including three-dimensional culture techniques, co-cultures with non-parenchymal cells, hepatospheres, precision cut liver slices and the isolated perfused liver. Also discussed is how closely hepatoma, stem cell and iPS cell-derived hepatocyte-like-cells resemble real hepatocytes. Finally, a summary is given of the state of the art of liver in vitro and mathematical modeling systems that are currently used in the pharmaceutical industry with an emphasis on drug metabolism, prediction of clearance, drug interaction, transporter studies and hepatotoxicity. One key message is that despite our enthusiasm for in vitro systems, we must never lose sight of the in vivo situation. Although hepatocytes have been isolated for decades, the hunt for relevant alternative systems has only just begun.

The heat shock response (HSR) is essential to survive acute proteotoxic stress and has been studied extensively in unicellular organisms and tissue culture cells, but to a lesser extent in intact metazoan animals. To identify the regulatory pathways that control the HSR in Caenorhabditis elegans, we performed a genome-wide RNAi screen and identified 59 genes corresponding to 7 positive activators required for the HSR and 52 negative regulators whose knockdown leads to constitutive activation of the HSR. These modifiers function in specific steps of gene expression, protein synthesis, protein folding, trafficking, and protein clearance, and comprise the metazoan heat shock regulatory network (HSN). Whereas the positive regulators function in all tissues of C. elegans, nearly all of the negative regulators exhibited tissue-selective effects. Knockdown of the subunits of the proteasome strongly induces HS reporter expression only in the intestine and spermatheca but not in muscle cells, while knockdown of subunits of the TRiC/CCT chaperonin induces HS reporter expression only in muscle cells. Yet, both the proteasome and TRiC/CCT chaperonin are ubiquitously expressed and are required for clearance and folding in all tissues. We propose that the HSN identifies a key subset of the proteostasis machinery that regulates the HSR according to the unique functional requirements of each tissue. The heat shock response (HSR) is an essential stress response that functions to maintain protein folding homeostasis, or proteostasis, and whose critical role in human diseases is recently becoming apparent. Previously, most of our understanding of the HSR has come from cultured cells and unicellular organisms. Here we present the identification of the heat shock regulatory network (HSN) in Caenorhabditis elegans, an intact, multicellular organism, using genome-wide RNAi screening. We identify 59 positive and negative regulators of the HSR, all of which have a previously established role in proteostasis, linking the function of the HSR to its regulation. Some HSN genes were previously established in other systems, many were indirectly linked to HSR, and others are novel. Unexpectedly, almost all negative regulators of the HSR act in distinct, tissue-selective patterns, despite their broad expression and universal cellular requirements. Therefore, our data indicate that the HSN consists of a specific subset of the proteostasis machinery that functions to link the proteostasis network to HSR regulation in a tissue-selective manner.

Nuclear receptor activation in liver leads to coordinated alteration of the expression of multiple gene products with attendant phenotypic changes of hepatocytes. Peroxisome proliferators including endogenous fatty acids, environmental chemicals, and drugs induce a multi-enzyme metabolic response that affects lipid and fatty acid processing. We studied the signaling network for the peroxisome proliferator-associated receptor alpha (PPARα) in primary human hepatocytes using the selective PPARα ligand, GW7647. We measured gene expression over multiple concentrations and times and conducted ChIP-seq studies at 2 and 24h to assess genomic binding of PPARα. Over all treatments there were 192 genes differentially expressed. Of these only 51% showed evidence of PPARα binding-either directly at PPARα response elements or via alternative mechanisms. Almost half of regulated genes had no PPARα binding. We then developed two novel bioinformatics methods to visualize the dose-dependent activation of both the transcription factor circuitry for PPARα and the downstream metabolic network in relation to functional annotation categories. Available databases identified several key transcription factors involved with the non-genomic targets after GW7647 treatment, including SP1, STAT1, ETS1, ERα, and HNF4α. The linkage from PPARα binding through gene expression likely requires intermediate protein kinases to activate these transcription factors. We found enrichment of functional annotation categories for organic acid metabolism and cell lipid metabolism among the differentially expressed genes. Lipid transport processes showed enrichment at the highest concentration of GW7647 (10 μM). While our strategy for mapping transcriptional networks is evolving, these approaches are necessary in moving from toxicogenomic methods that derive signatures of activity to methods that establish pathway structure, showing the coordination of the activated nuclear receptor with other signaling pathways.

We provide an overview of computational systems biology approaches as applied to the study of chemical- and drug-induced toxicity. The concept of “toxicity pathways” is described in the context of the 2007 US National Academies of Science report, “Toxicity testing in the 21st Century: A Vision and A Strategy.” Pathway mapping and modeling based on network biology concepts are a key component of the vision laid out in this report for a more biologically based analysis of dose-response behavior and the safety of chemicals and drugs. We focus on toxicity of the liver (hepatotoxicity) – a complex phenotypic response with contributions from a number of different cell types and biological processes. We describe three case studies of complementary multi-scale computational modeling approaches to understand perturbation of toxicity pathways in the human liver as a result of exposure to environmental contaminants and specific drugs. One approach involves development of a spatial, multicellular “virtual tissue” model of the liver lobule that combines molecular circuits in individual hepatocytes with cell–cell interactions and blood-mediated transport of toxicants through hepatic sinusoids, to enable quantitative, mechanistic prediction of hepatic dose-response for activation of the aryl hydrocarbon receptor toxicity pathway. Simultaneously, methods are being developing to extract quantitative maps of intracellular signaling and transcriptional regulatory networks perturbed by environmental contaminants, using a combination of gene expression and genome-wide protein-DNA interaction data. A predictive physiological model (DILIsym™) to understand drug-induced liver injury (DILI), the most common adverse event leading to termination of clinical development programs and regulatory actions on drugs, is also described. The model initially focuses on reactive metabolite-induced DILI in response to administration of acetaminophen, and spans multiple biological scales.

A major goal of systems biology is to understand how organism-level behavior arises from a myriad of molecular interactions. Often this involves complex sets of rules describing interactions among a large number of components. As an alternative, we have developed a simple, macro-level model to describe how chronic temperature stress affects reproduction in C. elegans. Our approach uses fundamental engineering principles, together with a limited set of experimentally derived facts, and provides quantitatively accurate predictions of performance under a range of physiologically relevant conditions. We generated detailed time-resolved experimental data to evaluate the ability of our model to describe the dynamics of C. elegans reproduction. We find considerable heterogeneity in responses of individual animals to heat stress, which can be understood as modulation of a few processes and may represent a strategy for coping with the ever-changing environment. Our experimental results and model provide quantitative insight into the breakdown of a robust biological system under stress and suggest, surprisingly, that the behavior of complex biological systems may be determined by a small number of key components. Dynamic response to changing conditions in the environment is an essential property of all biological systems. Whereas extensive research over the last several decades has elucidated numerous molecular responses to environmental stress, there is much less known how these translate into organismal-level responses. Two types of modeling approaches are often used to bridge this gap. Fine-grained models seek to explain phenomena as resulting from interactions of large numbers of individual components. This approach demands a highly detailed knowledge of the underlying molecular mechanisms and has an inherent difficulty in crossing spatial scales and organizational hierarchies. As an alternative, here we present a macro-level model of reproduction in C. elegans that uses fundamental engineering principles, together with a limited set of experimentally derived facts, to provide quantitatively accurate predictions of performance under a range of physiologically relevant conditions. One important finding is that individuals within a population display considerable heterogeneity in their response to heat stress. This could be a reflection of different strategies for coping with the ever-changing environment. Our study further demonstrates that dynamic behaviors of systems may be determined by a small number of key components that lead to the emergence of organismal phenomena.

In the past 10 years, the public, private, and non-profit sectors have found agreement that hazard identification and risk assessment should capitalize on the explosion of knowledge in the biological sciences, moving away from in life animal testing toward more human-relevant in vitro and in silico methods, collectively referred to as new approach meth­odologies (NAMs). The goals for implementation of NAMs are to efficiently identify possible chemical hazards and to gather dose-response data to inform more human-relevant safety assessment. While work proceeds to develop NAMs, there has been less emphasis on creating decision criteria or showing how risk context should guide selection and use of NAMs. Here, we outline application scenarios for NAMs in different risk contexts and place different NAMs and conven­tional testing approaches into four broad levels. Level 1 relies solely on computational screening; Level 2 consists of high throughput in vitro screening with human cells intended to provide broad coverage of possible responses; Level 3 focuses on fit-for-purpose assays selected based on presumptive modes of action (MOA) and designed to provide more quan­titative estimates of relevant dose responses; Level 4 has a variety of more complex multi-dimensional or multi-cellular assays and might include targeted in vivo studies to further define MOA. Each level also includes decision-appropriate exposure assessment tools. Our aims here are to (1) foster discussion about context-dependent applications of NAMs in relation to risk assessment needs and (2) describe a functional roadmap to identify where NAMs are expected to be ade­quate for chemical safety decision-making.