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.

Multi-component signal transduction pathways and gene regulatory circuits underpin integrated cellular responses to perturbations. A recurring set of network motifs serve as the basic building blocks of these molecular signalling networks. This review focuses on ultrasensitive response motifs (URMs) that amplify small percentage changes in the input signal into larger percentage changes in the output response. URMs generally possess a sigmoid input–output relationship that is steeper than the Michaelis–Menten type of response and is often approximated by the Hill function. Six types of URMs can be commonly found in intracellular molecular networks and each has a distinct kinetic mechanism for signal amplification. These URMs are: (i) positive cooperative binding, (ii) homo-multimerization, (iii) multistep signalling, (iv) molecular titration, (v) zero-order covalent modification cycle and (vi) positive feedback. Multiple URMs can be combined to generate highly switch-like responses. Serving as basic signal amplifiers, these URMs are essential for molecular circuits to produce complex nonlinear dynamics, including multistability, robust adaptation and oscillation. These dynamic properties are in turn responsible for higher-level cellular behaviours, such as cell fate determination, homeostasis and biological rhythm.

The approaches to quantitatively assessing the health risks of chemical exposure have not changed appreciably in the past 50 to 80 years, the focus remaining on high-dose studies that measure adverse outcomes in homogeneous animal populations. This expensive, low-throughput approach relies on conservative extrapolations to relate animal studies to much lower-dose human exposures and is of questionable relevance to predicting risks to humans at their typical low exposures. It makes little use of a mechanistic understanding of the mode of action by which chemicals perturb biological processes in human cells and tissues. An alternative vision, proposed by the U.S. National Research Council (NRC) report Toxicity Testing in the 21^st Century: A Vision and a Strategy, called for moving away from traditional high-dose animal studies to an approach based on perturbation of cellular responses using well-designed in vitro assays. Central to this vision are (a) “toxicity pathways” (the innate cellular pathways that may be perturbed by chemicals) and (b) the determination of chemical concentration ranges where those perturbations are likely to be excessive, thereby leading to adverse health effects if present for a prolonged duration in an intact organism. In this paper we briefly review the original NRC report and responses to that report over the past 3 years, and discuss how the change in testing might be achieved in the U.S. and in the European Union (EU). EU initiatives in developing alternatives to animal testing of cosmetic ingredients have run very much in parallel with the NRC report. Moving from current practice to the NRC vision would require using prototype toxicity pathways to develop case studies showing the new vision in action. In this vein, we also discuss how the proposed strategy for toxicity testing might be applied to the toxicity pathways associated with DNA damage and repair.

The image of the "epigenetic landscape", with a series of branching valleys and ridges depicting stable cellular states and the barriers between those states, has been a popular visual metaphor for cell lineage specification - especially in light of the recent discovery that terminally differentiated adult cells can be reprogrammed into pluripotent stem cells or into alternative cell lineages. However the question of whether the epigenetic landscape can be mapped out quantitatively to provide a predictive model of cellular differentiation remains largely unanswered.

Pluripotent stem cell-derived organoids can recapitulate significant features of organ development in vitro. We hypothesized that creating human heart organoids by mimicking aspects of in utero gestation (e.g., addition of metabolic and hormonal factors) would lead to higher physiological and anatomical relevance. We find that heart organoids produced using this self-organization-driven developmental induction strategy are remarkably similar transcriptionally and morphologically to age-matched human embryonic hearts. We also show that they recapitulate several aspects of cardiac development, including large atrial and ventricular chambers, proepicardial organ formation, and retinoic acid-mediated anterior-posterior patterning, mimicking the developmental processes found in the post-heart tube stage primitive heart. Moreover, we provide proof-of-concept demonstration of the value of this system for disease modeling by exploring the effects of ondansetron, a drug administered to pregnant women and associated with congenital heart defects. These findings constitute a significant technical advance for synthetic heart development and provide a powerful tool for cardiac disease modeling.

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.