Spatio-temporal visualization of the distribution of acetaminophen as well as its metabolites and adducts in mouse livers by MALDI MSI

MALDI Imaging Mass Spectrometry Visualizes the Distribution of Antidepressant Duloxetine and Its Major Metabolites in Mouse Brain, Liver, Kidney, and Spleen Tissues

Imaging mass spectrometry (IMS) is a powerful tool for mapping the spatial distribution of unlabeled drugs and metabolites that may find application in assessing drug delivery, explaining drug efficacy, and identifying potential toxicity. This study focuses on determining the spatial distribution of the antidepressant duloxetine, which is widely prescribed despite common adverse effects (liver injury, constant headaches) whose mechanisms are not fully understood. We used high-resolution IMS with matrix-assisted laser desorption/ionization to examine the distribution of duloxetine and its major metabolites in four mouse organs where it may contribute to efficacy or toxicity: brain, liver, kidney, and spleen. In none of these tissues is duloxetine or its metabolites homogeneously distributed, which has implications for both efficacy and toxicity. We found duloxetine to be similarly distributed in spleen red pulp and white pulp but differentially distributed in different anatomic regions of the liver, kidney, and brain, with dose-dependent patterns. Comparison with hematoxylin and eosin staining of tissue sections reveals that the ion images of endogenous lipids help delineate anatomic regions in the brain and kidney, while heme ion images assist in differentiating regions within the spleen. These endogenous metabolites may serve as a valuable resource for examining the spatial distribution of other drugs in tissues when staining images are not available. These findings may facilitate future mechanistic studies of the therapeutic and adverse effects of duloxetine. In the current work, we did not perform absolute quantification of duloxetine, which will be reported in due course. SIGNIFICANCE STATEMENT: The study utilized imaging mass spectrometry to examine the spatial distribution of duloxetine and its primary metabolites in mouse brain, liver, kidney, and spleen. These results may pave the way for future investigations into the mechanisms behind duloxetine's therapeutic and adverse effects. Furthermore, the mass spectrometry images of specific endogenous metabolites such as heme could be valuable in analyzing the spatial distribution of other drugs within tissues in scenarios where histological staining images are unavailable.

Acetaminophen overdose causes a breach of the blood-bile barrier in mice but not in rats

Acetaminophen (APAP) is known to cause a breach of the blood-bile barrier in mice that, via a mechanism called futile bile acid (BA) cycling, increases BA concentrations in hepatocytes above cytotoxic thresholds. Here, we compared this mechanism in mice and rats, because both species differ massively in their susceptibility to APAP and compared the results to available human data. Dose and time-dependent APAP experiments were performed in male C57BL6/N mice and Wistar rats. The time course of BA concentrations in liver tissue and in blood was analyzed by MALDI-MSI and LC-MS/MS. APAP and its derivatives were measured in the blood by LC-MS. APAP-induced liver damage was analyzed by histopathology, immunohistochemistry, and by clinical chemistry. In mice, a transient increase of BA in blood and in peri-central hepatocytes preceded hepatocyte death. The BA increase coincided with oxidative stress in liver tissue and a compromised morphology of bile canaliculi and immunohistochemically visualized tight junction proteins. Rats showed a reduced metabolic activation of APAP compared to mice. However, even at very high doses that caused cell death of hepatocytes, no increase of BA concentrations was observed neither in liver tissue nor in the blood. Correspondingly, no oxidative stress was detectable, and the morphology of bile canaliculi and tight junction proteins remained unaltered. In conclusion, different mechanisms cause cell death in rats and mice, whereby oxidative stress and a breach of the blood-bile barrier are seen only in mice. Since transient cholestasis also occurs in human patients with APAP overdose, mice are a clinically relevant species to study APAP hepatotoxicity but not rats.

A liver digital twin for in silico testing of cellular and inter-cellular mechanisms in regeneration after drug-induced damage

This communication presents a mathematical mechanism-based model of the regenerating liver after drug-induced pericentral lobule damage resolving tissue microarchitecture. The consequence of alternative hypotheses about the interplay of different cell types on regeneration was simulated. Regeneration dynamics has been quantified by the size of the damage-induced dead cell area, the hepatocyte density and the spatial-temporal profile of the different cell types. We use deviations of observed trajectories from the simulated system to identify branching points, at which the systems behavior cannot be explained by the underlying set of hypotheses anymore. Our procedure reflects a successful strategy for generating a fully digital liver twin that, among others, permits to test perturbations from the molecular up to the tissue scale. The model simulations are complementing current knowledge on liver regeneration by identifying gaps in mechanistic relationships and guiding the system toward the most informative (lacking) parameters that can be experimentally addressed.

Biomarker discovery in acetaminophen hepatotoxicity: leveraging single-cell transcriptomics and mechanistic insight

INTRODUCTION

Acetaminophen (APAP) overdose is the leading cause of drug-induced liver injury and can cause a rapid progression to acute liver failure (ALF). Therefore, the identification of prognostic biomarkers to determine which patients will require a liver transplant is critical for APAP-induced ALF.

AREAS COVERED

We begin by relating the mechanistic investigations in mouse models of APAP hepatotoxicity to the human APAP overdose pathophysiology. We draw insights from the established sequence of molecular events in mice to understand the progression of events in the APAP overdose patient. Through this mechanistic understanding, several new biomarkers, such as CXCL14, have recently been evaluated. We also explore how single-cell RNA sequencing, spatial transcriptomics, and other omics approaches have been leveraged for identifying novel biomarkers and how these approaches will continue to push the field of biomarker discovery forward.

EXPERT OPINION

Recent investigations have elucidated several new biomarkers or combination of markers such as CXCL14, a regenerative miRNA signature, a cell death miRNA signature, hepcidin, LDH, CPS1, and FABP1. While these biomarkers are promising, they all require further validation. Larger cohort studies analyzing these new biomarkers in the same patient samples, while adding these candidate biomarkers to prognostic models will further support their clinical utility.

Bioactivation and reactivity research advances - 2022 year in review‡

With the 50th year mark since the launch of Drug Metabolism and Disposition journal, the field of drug metabolism and bioactivation has advanced exponentially in the past decades (Guengerich 2023).This has, in a major part, been due to the continued advances across the whole spectrum of applied technologies in hardware, software, machine learning (ML), and artificial intelligence (AI). LC-MS platforms continue to evolve to support key applications in the field, and automation is also improving the accuracy, precision, and throughput of these supporting assays. In addition, sample generation and processing is being aided by increased diversity and quality of reagents and bio-matrices so that what is being analyzed is more relevant and translatable. The application of in silico platforms (applied software, ML, and AI) is also making great strides, and in tandem with the more traditional approaches mentioned previously, is significantly advancing our understanding of bioactivation pathways and how these play a role in toxicity. All of this continues to allow the area of bioactivation to evolve in parallel with associated fields to help bring novel or improved medicines to patients with urgent or unmet needs.Shuai Wang and Cyrus Khojasteh, on behalf of the authors.

Basic concepts of mixture toxicity and relevance for risk evaluation and regulation

Exposure to multiple substances is a challenge for risk evaluation. Currently, there is an ongoing debate if generic "mixture assessment/allocation factors" (MAF) should be introduced to increase public health protection. Here, we explore concepts of mixture toxicity and the potential influence of mixture regulation concepts for human health protection. Based on this analysis, we provide recommendations for research and risk assessment. One of the concepts of mixture toxicity is additivity. Substances may act additively by affecting the same molecular mechanism within a common target cell, for example, dioxin-like substances. In a second concept, an "enhancer substance" may act by increasing the target site concentration and aggravating the adverse effect of a "driver substance". For both concepts, adequate risk management of individual substances can reliably prevent adverse effects to humans. Furthermore, we discuss the hypothesis that the large number of substances to which humans are exposed at very low and individually safe doses may interact to cause adverse effects. This commentary identifies knowledge gaps, such as the lack of a comprehensive overview of substances regulated under different silos, including food, environmentally and occupationally relevant substances, the absence of reliable human exposure data and the missing accessibility of ratios of current human exposure to threshold values, which are considered safe for individual substances. Moreover, a comprehensive overview of the molecular mechanisms and most susceptible target cells is required. We conclude that, currently, there is no scientific evidence supporting the need for a generic MAF. Rather, we recommend taking more specific measures, which focus on compounds with relatively small ratios between human exposure and doses, at which adverse effects can be expected.

Detection of Distinct Distributions of Acetaminophen and Acetaminophen-Cysteine in Kidneys up to 10 μm Resolution and Identification of a Novel Acetaminophen Metabolite Using an AP-MALDI Imaging Mass Microscope

Drug distribution studies in tissue are crucial for understanding the pharmacokinetics and potential toxicity of drugs. Recently, matrix-assisted laser desorption ionization-mass spectrometry imaging (MALDI-MSI) has gained attention for drug distribution studies due to its high sensitivity, label-free nature, and ability to distinguish between parent drugs, their metabolites, and endogenous molecules. Despite these advantages, achieving high spatial resolution in drug imaging is challenging. Importantly, many drugs and metabolites are rarely detectable by conventional vacuum MALDI-MSI because of their poor ionization efficiency. It has been reported that acetaminophen (APAP) and one of its major metabolites, APAP-Cysteine (APAP-CYS), cannot be detected by vacuum MALDI-MSI without derivatization. In this context, we showed the distribution of both APAP and APAP-CYS in kidneys at high spatial resolution (25 and 10 μm) by employing an atmospheric pressure-MALDI imaging mass microscope without derivatization. APAP was highly accumulated in the renal pelvis 1 h after drug administration, while APAP-CYS exhibited characteristic distributions in the outer medulla and renal pelvis at both 30 min and 1 h after administration. Interestingly, cluster-like distributions of APAP and APAP-CYS were observed in the renal pelvis at 10 μm spatial resolution. Additionally, a novel APAP metabolite, tentatively coined as APAP-butyl sulfate (APAP-BS), was identified in the kidney, brain, and liver by combining MSI and tandem MSI. For the first time, our study revealed differential distributions of APAP, APAP-CYS (in kidneys), and APAP-BS (in kidney, brain, and liver) and is believed to enhance the understanding of the pharmacokinetics and potential nephrotoxicity of this drug.

In vitro to in vivo acetaminophen hepatotoxicity extrapolation using classical schemes, pharmacodynamic models and a multiscale spatial-temporal liver twin

In vitro to in vivo extrapolation represents a critical challenge in toxicology. In this paper we explore extrapolation strategies for acetaminophen (APAP) based on mechanistic models, comparing classical (CL) homogeneous compartment pharmacodynamic (PD) models and a spatial-temporal (ST), multiscale digital twin model resolving liver microarchitecture at cellular resolution. The models integrate consensus detoxification reactions in each individual hepatocyte. We study the consequences of the two model types on the extrapolation and show in which cases these models perform better than the classical extrapolation strategy that is based either on the maximal drug concentration (Cmax) or the area under the pharmacokinetic curve (AUC) of the drug blood concentration. We find that an CL-model based on a well-mixed blood compartment is sufficient to correctly predict the in vivo toxicity from in vitro data. However, the ST-model that integrates more experimental information requires a change of at least one parameter to obtain the same prediction, indicating that spatial compartmentalization may indeed be an important factor.

Human Wharton's Jelly-derived mesenchymal stem cells prevent acetaminophen-induced liver injury in a mouse model unlike human dermal fibroblasts

The persistence of hepatotoxicity induced by N-acetyl-para-aminophenol (Acetaminophen or Paracetamol, abbreviated as APAP) as the most common cause of acute liver failure in the United States, despite the availability of N-acetylcysteine, illustrates the clinical relevance of additional therapeutic approaches. While human mesenchymal stem cells (MSCs) have shown protection in mouse models of liver injury, the MSCs used are generally not cleared for human use and it is unclear whether these effects are due to xenotransplantation. Here we evaluated GMP manufactured clinical grade human Wharton's Jelly mesenchymal stem cells (WJMSCs), which are currently being investigated in human clinical trials, in a mouse model of APAP hepatotoxicity in comparison to human dermal fibroblasts (HDFs) to address these issues. C57BL6J mice were treated with a moderate APAP overdose (300 mg/kg) and WJMSCs were administered 90 min later. Liver injury was evaluated at 6 and 24 h after APAP. WJMSCs treatment reduced APAP-induced liver injury at both time points unlike HDFs, which showed no protection. APAP-induced JNK activation as well as AIF and Smac release from mitochondria were prevented by WJMSCs treatment without influencing APAP bioactivation. Mechanistically, WJMSCs treatment upregulated expression of Gclc and Gclm to enhance recovery of liver GSH levels to attenuate mitochondrial dysfunction and accelerated recovery of pericentral hepatocytes to re-establish liver zonation and promote liver homeostasis. Notably, preventing GSH resynthesis with buthionine sulfoximine prevented the protective effects of WJMSCs. These data indicate that these GMP-manufactured WJMCs could be a clinically relevant therapeutic approach in the management of APAP hepatotoxicity in humans.

Desorption Electrospray Ionization Mass Spectrometry Imaging Allows Spatial Localization of Changes in Acetaminophen Metabolism in the Liver after Intervention with 4-Methylpyrazole

Acetaminophen (APAP) overdose is the most common cause of acute liver failure in the US, and hepatotoxicity is initiated by a reactive metabolite which induces characteristic centrilobular necrosis. The only clinically available antidote is N-acetylcysteine, which has limited efficacy, and we have identified 4-methylpyrazole (4MP, Fomepizole) as a strong alternate therapeutic option, protecting against generation and downstream effects of the cytotoxic reactive metabolite in the clinically relevant C57BL/6J mouse model and in humans. However, despite the regionally restricted necrosis after APAP, our earlier studies on APAP metabolites in biofluids or whole tissue homogenate lack the spatial information needed to understand region-specific consequences of reactive metabolite formation after APAP overdose. Thus, to gain insight into the regional variation in APAP metabolism and study the influence of 4MP, we established a desorption electrospray ionization mass spectrometry imaging (DESI-MSI) platform for generation of ion images for APAP and its metabolites under ambient air, without chemical labeling or a prior coating of tissue which reduces chemical interference and perturbation of small molecule tissue localization. The spatial intensity and distribution of both oxidative and nonoxidative APAP metabolites were determined from mouse liver sections after a range of APAP overdoses. Importantly, exclusive differential signal intensities in metabolite abundance were noted in the tissue microenvironment, and 4MP treatment substantially influenced this topographical distribution.

Inhibition of cytochrome P450 enhances the nephro- and hepatotoxicity of ochratoxin A

The mycotoxin ochratoxin A (OTA) is a contaminant in food that causes nephrotoxicity and to a minor degree hepatotoxicity. Recently, we observed that OTA induces liver damage preferentially to the cytochrome P450 (CYP)-expressing pericentral lobular zone, similar to hepatotoxic substances known to be metabolically toxified by CYP, such as acetaminophen or carbon tetrachloride. To investigate whether CYP influences OTA toxicity, we used a single dose of OTA (7.5 mg/kg; intravenous) with and without pre-treatment with the pan CYP-inhibitor 1-aminobenzotriazole (ABT) 2 h before OTA administration. Blood, urine, as well as liver and kidney tissue samples were collected 24 h after OTA administration for biochemical and histopathological analyses. Inhibition of CYPs by ABT strongly increased the nephro- and hepatotoxicity of OTA. The urinary kidney damage biomarkers kidney injury molecule-1 (KIM-1) and neutrophil gelatinase-associated lipocalin (NGAL) were increased > 126-fold and > 20-fold, respectively, in mice treated with ABT and OTA compared to those receiving OTA alone. The blood biomarkers of liver damage, alanine transaminase (ALT) and aspartate transaminase (AST) both increased > 21- and 30-fold, respectively, when OTA was administered to ABT pre-treated mice compared to the effect of OTA alone. Histological analysis of the liver revealed a pericentral lobular damage induced by OTA despite CYP-inhibition by ABT. Administration of ABT alone caused no hepato- or nephrotoxicity. Overall, the results presented are compatible with a scenario where CYPs mediate the detoxification of OTA, yet the mechanisms responsible for the pericental liver damage pattern still remain to be elucidated.

Advances in mass spectrometry imaging for toxicological analysis and safety evaluation of pharmaceuticals

Safety issues caused by pharmaceuticals have frequently occurred worldwide, posing a tremendous threat to human health. As an essential part of drug development, the toxicological analysis and safety evaluation is of great significance. In addition, the risk of pharmaceuticals accumulation in the environment and the monitoring of the toxicity from natural medicines have also received ongoing concerns. Due to a lack of spatial distribution information provided by common analytical methods, analyses that provide spatial dimensions could serve as complementary safety evaluation methods for better prediction and evaluation of drug toxicity. With advances in technical solutions and software algorithms, mass spectrometry imaging (MSI) has received increasing attention as a popular analytical tool that enables the simultaneous implementation of qualitative, quantitative, and localization without complex sample pretreatment and labeling steps. In recent years, MSI has become more attractive, powerful, and sensitive and has been applied in several scientific fields that can meet the safety assessment requirements. This review aims to cover a detailed summary of the various MSI technologies utilized in the biomedical and pharmaceutical area, including technical principles, advantages, current status, and future trends. Representative applications and developments in the safety-related issues of different pharmaceuticals and natural medicines are also described to provide a reference for pharmaceutical research, improve rational clinical medicine use, and ensure public safety.

Hypoalbuminemia affects the spatio-temporal tissue distribution of ochratoxin A in liver and kidneys: consequences for organ toxicity

Hypoalbuminemia (HA) is frequently observed in systemic inflammatory diseases and in liver disease. However, the influence of HA on the pharmacokinetics and toxicity of compounds with high plasma albumin binding remained insufficiently studied. The ‘lack-of-delivery-concept’ postulates that HA leads to less carrier mediated uptake of albumin bound substances into hepatocytes and to less glomerular filtration; in contrast, the ‘concept-of-higher-free-fraction’ argues that increased concentrations of non-albumin bound compounds facilitate hepatocellular uptake and enhance glomerular filtration. To address this question, we performed intravital imaging on livers and kidneys of anesthetized mice to quantify the spatio-temporal tissue distribution of the mycotoxin ochratoxin A (OTA) based on its auto-fluorescence in albumin knockout and wild-type mice. HA strongly enhanced the uptake of OTA from the sinusoidal blood into hepatocytes, followed by faster secretion into bile canaliculi. These toxicokinetic changes were associated with increased hepatotoxicity in heterozygous albumin knockout mice for which serum albumin was reduced to a similar extent as in patients with severe hypoalbuminemia. HA also led to a shorter half-life of OTA in renal capillaries, increased glomerular filtration, and to enhanced uptake of OTA into tubular epithelial cells. In conclusion, the results favor the ‘concept-of-higher-free-fraction’ in HA; accordingly, HA causes an increased tissue uptake of compounds with high albumin binding and increased organ toxicity. It should be studied if this concept can be generalized to all compounds with high plasma albumin binding that are substrates of hepatocyte and renal tubular epithelial cell carriers.

Design and Measurement of Drug Tissue Concentration Asymmetry and Tissue Exposure-Effect (Tissue PK-PD) Evaluation

The "free drug hypothesis" assumes that, in the absence of transporters, the steady state free plasma concentrations equal to that at the site of action that elicit pharmacologic effects. While it is important to utilize the free drug hypothesis, exceptions exist that the free plasma exposures, either at C_max, C_trough, and C_average, or at other time points, cannot represent the corresponding free tissue concentrations. This "drug concentration asymmetry" in both total and free form can influence drug disposition and pharmacological effects. In this review, we first discuss options to assess total and free drug concentrations in tissues. Then various drug design strategies to achieve concentration asymmetry are presented. Last, the utilities of tissue concentrations in understanding exposure-effect relationships and translational projections to humans are discussed for several therapeutic areas and modalities. A thorough understanding in plasma and tissue exposures correlation with pharmacologic effects can provide insightful guidance to aid drug discovery.

Spatial analysis of drug absorption, distribution, metabolism, and toxicology using mass spectrometry imaging

Mass spectrometry imaging (MSI) is emerging as a powerful analytical tool for detection, quantification, and simultaneous spatial molecular imaging of endogenous and exogenous molecules via in situ mass spectrometry analysis of thin tissue sections without the requirement of chemical labeling. The MSI generates chemically specific and spatially resolved ion distribution information for administered drugs and metabolites, which allows numerous applications for studies involving various stages of drug absorption, distribution, metabolism, excretion, and toxicity (ADMET). MSI-based pharmacokinetic imaging analysis provides a histological context and cellular environment regarding dynamic drug distribution and metabolism processes, and facilitates the understanding of the spatial pharmacokinetics and pharmacodynamic properties of drugs. Herein, we discuss the MSI's current technological developments that offer qualitative, quantitative, and spatial location information of small molecule drugs, antibody, and oligonucleotides macromolecule drugs, and their metabolites in preclinical and clinical tissue specimens. We highlight the macro and micro drug-distribution in the whole-body, brain, lung, liver, kidney, stomach, intestine tissue sections, organoids, and the latest applications of MSI in pharmaceutical ADMET studies.

From Seeing to Simulating: A Survey of Imaging Techniques and Spatially-Resolved Data for Developing Multiscale Computational Models of Liver Regeneration

Liver regeneration, which leads to the re-establishment of organ mass, follows a specifically organized set of biological processes acting on various time and length scales. Computational models of liver regeneration largely focused on incorporating molecular and signaling detail have been developed by multiple research groups in the recent years. These modeling efforts have supported a synthesis of disparate experimental results at the molecular scale. Incorporation of tissue and organ scale data using noninvasive imaging methods can extend these computational models towards a comprehensive accounting of multiscale dynamics of liver regeneration. For instance, microscopy-based imaging methods provide detailed histological information at the tissue and cellular scales. Noninvasive imaging methods such as ultrasound, computed tomography and magnetic resonance imaging provide morphological and physiological features including volumetric measures over time. In this review, we discuss multiple imaging modalities capable of informing computational models of liver regeneration at the organ-, tissue- and cellular level. Additionally, we discuss available software and algorithms, which aid in the analysis and integration of imaging data into computational models. Such models can be generated or tuned for an individual patient with liver disease. Progress towards integrated multiscale models of liver regeneration can aid in prognostic tool development for treating liver disease.

Metabolite alteration analysis of acetaminophen-induced liver injury using a mass microscope

Acetaminophen (APAP)-induced liver injury (APAP-ILI), which occurs during APAP overdose, has been extensively studied. The production of N-acetyl-p-benzoquinone imine (NAPQI), the reactive metabolite of APAP, primarily contributes to liver injury. However, the mechanism underlying APAP-ILI has not been fully characterized. For further clarification, it is important to consider drug localization and endogenous substances in the injured liver. Herein, we show the localization of NAPQI metabolites and the injury site-specific changes in endogenous substances in the rat liver following APAP overdose using a mass microscope. Our results of on-tissue derivatization matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) showed that the glutathione metabolite of APAP, a detoxified metabolite of NAPQI, localized in the damaged central vein region in the rat liver following APAP administration. Moreover, in the conventional MALDI-MSI, the intensities of some phospholipids, phosphocreatine, and ceramides decreased or increased in the damaged regions compared with those in non-damaged regions. Phosphocreatine was localized in the damaged cells, whereas its related mitochondrial creatine kinase was localized in the non-damaged cells. These results are expected to contribute to further elucidation of the mechanisms underlying APAP-ILI. Our findings illustrate the localization of NAPQI-related metabolites and endogenous molecules associated with APAP-ILI, which may be related to apoptosis or metabolic adaptation ultimately protecting the cells. As MALDI-MSI can analyze and differentiate regions with tissue damage, it is a valuable tool for analyzing the mechanism underlying drug-induced liver injury to identify novel biomarkers.

Automated Detection of Portal Fields and Central Veins in Whole-Slide Images of Liver Tissue

Many physiological processes and pathological phenomena in the liver tissue are spatially heterogeneous. At a local scale, biomarkers can be quantified along the axis of the blood flow, from portal fields (PFs) to central veins (CVs), i.e., in zonated form. This requires detecting PFs and CVs. However, manually annotating these structures in multiple whole-slide images is a tedious task. We describe and evaluate a fully automated method, based on a convolutional neural network, for simultaneously detecting PFs and CVs in a single stained section. Trained on scans of hematoxylin and eosin-stained liver tissue, the detector performed well with an F1 score of 0.81 compared to annotation by a human expert. It does, however, not generalize well to previously unseen scans of steatotic liver tissue with an F1 score of 0.59. Automated PF and CV detection eliminates the bottleneck of manual annotation for subsequent automated analyses, as illustrated by two proof-of-concept applications: We computed lobulus sizes based on the detected PF and CV positions, where results agreed with published lobulus sizes. Moreover, we demonstrate the feasibility of zonated quantification of biomarkers detected in different stainings based on lobuli and zones obtained from the detected PF and CV positions. A negative control (hematoxylin and eosin) showed the expected homogeneity, a positive control (glutamine synthetase) was quantified to be strictly pericentral, and a plausible zonation for a heterogeneous F4/80 staining was obtained. Automated detection of PFs and CVs is one building block for automatically quantifying physiologically relevant heterogeneity of liver tissue biomarkers. Perspectively, a more robust and automated assessment of zonation from whole-slide images will be valuable for parameterizing spatially resolved models of liver metabolism and to provide diagnostic information.

Applications of spatially resolved omics in the field of endocrine tumors

Endocrine tumors derive from endocrine cells with high heterogeneity in function, structure and embryology, and are characteristic of a marked diversity and tissue heterogeneity. There are still challenges in analyzing the molecular alternations within the heterogeneous microenvironment for endocrine tumors. Recently, several proteomic, lipidomic and metabolomic platforms have been applied to the analysis of endocrine tumors to explore the cellular and molecular mechanisms of tumor genesis, progression and metastasis. In this review, we provide a comprehensive overview of spatially resolved proteomics, lipidomics and metabolomics guided by mass spectrometry imaging and spatially resolved microproteomics directed by microextraction and tandem mass spectrometry. In this regard, we will discuss different mass spectrometry imaging techniques, including secondary ion mass spectrometry, matrix-assisted laser desorption/ionization and desorption electrospray ionization. Additionally, we will highlight microextraction approaches such as laser capture microdissection and liquid microjunction extraction. With these methods, proteins can be extracted precisely from specific regions of the endocrine tumor. Finally, we compare applications of proteomic, lipidomic and metabolomic platforms in the field of endocrine tumors and outline their potentials in elucidating cellular and molecular processes involved in endocrine tumors.

Innovation in drug toxicology: Application of mass spectrometry imaging technology

Mass spectrometry imaging (MSI) is a powerful molecular imaging technology that can obtain qualitative, quantitative, and location information by simultaneously detecting and mapping endogenous or exogenous molecules in biological tissue slices without specific chemical labeling or complex sample pretreatment. This article reviews the progress made in MSI and its application in drug toxicology research, including the tissue distribution of toxic drugs and their metabolites, the target organs (liver, kidney, lung, eye, and central nervous system) of toxic drugs, the discovery of toxicity-associated biomarkers, and explanations of the mechanisms of drug toxicity when MSI is combined with the cutting-edge omics methodologies. The unique advantages and broad prospects of this technology have been fully demonstrated to further promote its wider use in the field of pharmaceutical toxicology.

Study of the Distribution of Acetaminophen and Its Metabolites in Rats, from the Whole-Body to Isolated Organ Levels, by Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging after On-Tissue Chemical Derivatization

During drug development, detailed investigations of the pharmacokinetic profile of the drug are required to characterize its absorption, distribution, metabolism, and excretion properties. Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) is an established technique for studies of the distribution of drugs and their metabolites. It has advantages over autoradiography, which is conventionally used for distribution studies: it does not require the radiolabeling of drugs and can distinguish between the drug and its metabolites directly in the tissue. However, its lack of sensitivity in certain cases remains challenging. Novel procedures, such as on-tissue chemical derivatization (OTCD), could be developed to increase sensitivity. We used OTCD to enhance the sensitivity of MALDI-MSI for one of the most widely used drugs, acetaminophen, and to study its distribution in tissues. Without derivatization, this drug and some of its metabolites are undetectable by MALDI-MSI in the tissues of treated rats. We used 2-fluoro-1-methylpyridinium p-toluene sulfonate as a derivatization reagent, to increase the ionization yield of acetaminophen and some of its metabolites. The OTCD protocol made it possible to study the distribution of acetaminophen and its metabolites in whole-body sections at a spatial resolution of 400 μm and in complex anatomical structures, such as the testis and epididymis, at a spatial resolution <50 μm. The OTCD is also shown to be compatible with the quantification of acetaminophen by MALDI-MSI in whole-body tissues. This protocol could be applied to other molecules bearing phenol groups and presenting a low ionization efficiency.

Spatio-Temporal Multiscale Analysis of Western Diet-Fed Mice Reveals a Translationally Relevant Sequence of Events during NAFLD Progression

Mouse models of non-alcoholic fatty liver disease (NAFLD) are required to define therapeutic targets, but detailed time-resolved studies to establish a sequence of events are lacking. Here, we fed male C57Bl/6N mice a Western or standard diet over 48 weeks. Multiscale time-resolved characterization was performed using RNA-seq, histopathology, immunohistochemistry, intravital imaging, and blood chemistry; the results were compared to human disease. Acetaminophen toxicity and ammonia metabolism were additionally analyzed as functional readouts. We identified a sequence of eight key events: formation of lipid droplets; inflammatory foci; lipogranulomas; zonal reorganization; cell death and replacement proliferation; ductular reaction; fibrogenesis; and hepatocellular cancer. Functional changes included resistance to acetaminophen and altered nitrogen metabolism. The transcriptomic landscape was characterized by two large clusters of monotonously increasing or decreasing genes, and a smaller number of 'rest-and-jump genes' that initially remained unaltered but became differentially expressed only at week 12 or later. Approximately 30% of the genes altered in human NAFLD are also altered in the present mouse model and an increasing overlap with genes altered in human HCC occurred at weeks 30-48. In conclusion, the observed sequence of events recapitulates many features of human disease and offers a basis for the identification of therapeutic targets.

Transcriptomic Cross-Species Analysis of Chronic Liver Disease Reveals Consistent Regulation Between Humans and Mice

Mouse models are frequently used to study chronic liver diseases (CLDs). To assess their translational relevance, we quantified the similarity of commonly used mouse models to human CLDs based on transcriptome data. Gene‐expression data from 372 patients were compared with data from acute and chronic mouse models consisting of 227 mice, and additionally to nine published gene sets of chronic mouse models. Genes consistently altered in humans and mice were mapped to liver cell types based on single‐cell RNA‐sequencing data and validated by immunostaining. Considering the top differentially expressed genes, the similarity between humans and mice varied among the mouse models and depended on the period of damage induction. The highest recall (0.4) and precision (0.33) were observed for the model with 12‐months damage induction by CCl₄ and by a Western diet, respectively. Genes consistently up‐regulated between the chronic CCl₄ model and human CLDs were enriched in inflammatory and developmental processes, and mostly mapped to cholangiocytes, macrophages, and endothelial and mesenchymal cells. Down‐regulated genes were enriched in metabolic processes and mapped to hepatocytes. Immunostaining confirmed the regulation of selected genes and their cell type specificity. Genes that were up‐regulated in both acute and chronic models showed higher recall and precision with respect to human CLDs than exclusively acute or chronic genes. Conclusion: Similarly regulated genes in human and mouse CLDs were identified. Despite major interspecies differences, mouse models detected 40% of the genes significantly altered in human CLD. The translational relevance of individual genes can be assessed at https://saezlab.shinyapps.io/liverdiseaseatlas/.

Spatial Reconstruction of the Early Hepatic Transcriptomic Landscape After an Acetaminophen Overdose Using Single-Cell RNA-Sequencing

An acetaminophen (APAP) overdose is the most common cause of acute liver failure in the United States. A hallmark characteristic of APAP hepatotoxicity is centrilobular necrosis. General, innate mechanisms such as lower amounts of GSH and higher cytochrome P450 2e1 expression in pericentral (PC) hepatocytes are known to contribute to the differences in susceptibility to cell injury between periportal (PP) hepatocytes and PC hepatocytes. Although a sequence of molecular events involving formation of the reactive metabolite N-acetyl-p-benzoquinone imine, GSH depletion, oxidative stress, and c-Jun N-terminal kinase activation define the early cell stress trajectory following APAP exposure, their activation in PC versus PP hepatocytes is not well characterized. By using single-cell RNA-sequencing, we provide the first reconstruction of the early transcriptomic APAP liver lobule after validation of our methodology using human liver single-cell RNA-sequencing data. Two hours after APAP treatment, we find that PP hepatocytes progress along the APAP stress axis to oxidative stress, before resolving injury due to innate and adaptive mechanisms. However, PC hepatocytes continue along this stress axis as indicated by activation of mitogen-activated protein kinase genes, which is absent in PP hepatocytes. We also identify a population of glutamine synthetase enriched PC hepatocytes in close proximity to the central vein, where a stepwise induction of a stress program culminated in cell death. Collectively, these findings elucidate a molecular sequence of events distinguishing the differential response to APAP exposure between PP and PC hepatocytes and identify a subset of uniquely susceptible PC hepatocytes.

Mitochondrial Membrane Potential Drives Early Change in Mitochondrial Morphology After Acetaminophen Exposure

Mitochondrial morphology plays a critical role in regulating mitochondrial and cellular function. It is well established that oxidative stress and mitochondrial injury are central to acetaminophen (APAP) hepatotoxicity. However, the role of mitochondrial dynamics, namely the remodeling of mitochondrial morphology through fusion and fission, has largely gone unexplored. To investigate this, we used primary mouse hepatocytes treated with APAP which allowed for real-time visualization of mitochondrial morphology using mitotracker green. We found that alterations in mitochondrial morphology were dose dependent, with a biphasic response in mitochondrial shape at higher APAP doses. Importantly, these two distinct mitochondrial morphologies corresponded with differences in mitochondrial respiratory function and polarization. The early change in mitochondrial morphology can be reversible and appears to be an adaptive response caused by alterations in membrane potential, which ultimately help preserve mitochondrial function. The later delayed change in mitochondrial morphology is irreversible and is driven by loss of mitochondrial membrane potential, decreased canonical fusion proteins, and alterations in mitochondrial lipid composition. Collectively, these later changes tilt the scales toward mitochondrial fission resulting in fragmented mitochondria with reduced functionality. This work provides evidence of adaptive early changes in mitochondrial morphology, which results in functional consequences that are dictated by the severity of APAP overdose.

Subcellular spatio-temporal intravital kinetics of aflatoxin B1 and ochratoxin A in liver and kidney

Local accumulation of xenobiotics in human and animal tissues may cause adverse effects. Large differences in their concentrations may exist between individual cell types, often due to the expression of specific uptake and export carriers. Here we established a two-photon microscopy-based technique for spatio-temporal detection of the distribution of mycotoxins in intact kidneys and livers of anesthetized mice with subcellular resolution. The mycotoxins ochratoxin A (OTA, 10 mg/kg b.w.) and aflatoxin B1 (AFB1, 1.5 mg/kg b.w.), which both show blue auto-fluorescence, were analyzed after intravenous bolus injections. Within seconds after administration, OTA was filtered by glomeruli, and enriched in distal tubular epithelial cells (dTEC). A striking feature of AFB1 toxicokinetics was its very rapid uptake from sinusoidal blood into hepatocytes (t1/2 ~ 4 min) and excretion into bile canaliculi. Interestingly, AFB1 was enriched in the nuclei of hepatocytes with zonal differences in clearance. In the cytoplasm of pericentral hepatocytes, the half-life (t1/2~ 63 min) was much longer compared to periportal hepatocytes of the same lobules (t1/2 ~ 9 min). In addition, nuclear AFB1 from periportal hepatocytes cleared faster compared to the pericentral region. These local differences in AFB1 clearance may be due to the pericentral expression of cytochrome P450 enzymes that activate AFB1 to protein- and DNA-binding metabolites. In conclusion, the present study shows that large spatio-temporal concentration differences exist within the same tissues and its analysis may provide valuable additional information to conventional toxicokinetic studies.

Multi-Purpose Ontology-Based Visualization of Spatio-Temporal Data: A Case Study on Silk Heritage

Due to the increasing use of data analytics, information visualization is getting more and more important. However, as data get more complex, so does visualization, often leading to ad hoc and cumbersome solutions. A recent alternative is the use of the so-called knowledge-assisted visualization tools. In this paper, we present STMaps (Spatio-Temporal Maps), a multipurpose knowledge-assisted ontology-based visualization tool of spatio-temporal data. STMaps has been (originally) designed to show, by means of an interactive map, the content of the SILKNOW project, a European research project on silk heritage. It is entirely based on ontology support, as it gets the source data from an ontology and uses also another ontology to define how data should be visualized. STMaps provides some unique features. First, it is a multi-platform application. It can work embedded in an HTML page and can also work as a standalone application over several computer architectures. Second, it can be used for multiple purposes by just changing its configuration files and/or the ontologies on which it works. As STMaps relies on visualizing spatio-temporal data provided by an ontology, the tool could be used to visualize the results of any domain (in other cultural and non-cultural contexts), provided that its datasets contain spatio-temporal information. The visualization mechanisms can also be changed by changing the visualization ontology. Third, it provides different solutions to show spatio-temporal data, and also deals with uncertain and missing information. STMaps has been tested to browse silk-related objects, discovering some interesting relationships between different objects, showing the versatility and power of the different visualization tools proposed in this paper. To the best of our knowledge, this is also the first ontology-based visualization tool applied to silk-related heritage.

Intestinal Dysbiosis Amplifies Acetaminophen-Induced Acute Liver Injury

BACKGROUND & AIMS

Acute liver failure (ALF) represents an unmet medical need in Western countries. Although the link between intestinal dysbiosis and chronic liver disease is well-established, there is little evidence for a functional role of gut-liver interaction during ALF. Here we hypothesized that intestinal dysbiosis may affect ALF.

METHODS

To test this hypothesis, we assessed the association of proton pump inhibitor (PPI) or long-term antibiotics (ABx) intake, which have both been linked to intestinal dysbiosis, and occurrence of ALF in the 500,000 participants of the UK BioBank population-based cohort. For functional studies, male Nlrp6-/- mice were used as a dysbiotic mouse model and injected with a sublethal dose of acetaminophen (APAP) or lipopolysaccharide (LPS) to induce ALF.

RESULTS

Multivariate Cox regression analyses revealed a significantly increased risk (odds ratio, 2.3-3) for developing ALF in UK BioBank participants with PPI or ABx. Similarly, dysbiotic Nlrp6-/- mice displayed exacerbated APAP- and LPS-induced liver injury, which was linked to significantly reduced gut and liver tissue microbiota diversity and correlated with increased intestinal permeability at baseline. Fecal microbiota transfer (FMT) from Nlrp6-/- mice into wild-type (WT) mice augmented liver injury on APAP treatment in recipient WT mice, resembling the inflammatory phenotype of Nlrp6-/- mice. Specifically, FMT skewed monocyte polarization in WT mice toward a Ly6Chi inflammatory phenotype, suggesting a critical function of these cells as sensors of gut-derived signals orchestrating the inflammatory response.

CONCLUSIONS

Our data show an important yet unknown function of intestinal microbiota during ALF. Intestinal dysbiosis was transferrable to healthy WT mice via FMT and aggravated liver injury. Our study highlights intestinal microbiota as a targetable risk factor for ALF.

Research mechanisms of and pharmaceutical treatments for ferroptosis in liver diseases

Regulated cell death (RCD) is a universal process in living organisms that is essential for tissue homeostasis or to the restoration of biological equilibrium following stress. Ferroptosis is a specific nonapoptotic cell death that is dependent on iron and is very different from other forms of RCD. Ferroptosis can affect the development of liver diseases such as drug-induced liver injury (DILI), liver fibrosis, and hepatocellular carcinoma (HCC) by regulating the level of intracellular iron, the production of intracellular reactive oxygen species, and lipid peroxides. In this review, we summarize the current knowledge of ferroptosis, in terms of discovery, history, characteristics, mechanism, and the factors regulating liver diseases. We discuss how these factors and signaling pathways change in the context of liver disease. Furthermore, we focus on delineating the roles of effective therapeutic drugs or compounds in liver disease. In summary, we discuss the role of ferroptosis in liver disease, providing a strategy and new ideas for preventing liver disease, finding new therapeutic targets, and reducing liver damage.

Aryl Hydrocarbon Receptor Activity in Hepatocytes Sensitizes to Hyperacute Acetaminophen-Induced Hepatotoxicity in Mice

BACKGROUND & AIMS

Acetaminophen (APAP)-induced liver injury is one of the most common causes of acute liver failure, however, a clear definition of sensitizing risk factors is lacking. Here, we investigated the role of the ligand-activated transcription factor aryl hydrocarbon receptor (Ahr) in APAP-induced liver injury. We hypothesized that Ahr, which integrates environmental, dietary, microbial and metabolic signals into complex cellular transcriptional programs, might act as a rheostat for APAP-toxicity.

METHODS

Wildtype or conditional Ahr knockout mice lacking Ahr in hepatocytes (Alb^Δ/ΔAhr) or myeloid cells (LysM^Δ/ΔAhr) were treated with the specific Ahr ligand 2-(1'H-indole-3'-carbonyl)-thiazole-4-carboxylic acid methyl ester (ITE) together with APAP.

RESULTS

Ahr activation by ITE, which by itself was non-toxic, exacerbated APAP-induced hepatotoxicity compared to vehicle-treated controls, causing 80% vs. 0% mortality after administration of a normally sublethal APAP overdose. Of note, Ahr activation induced hepatocyte death even at APAP doses within the therapeutic range. Aggravated liver injury was associated with significant neutrophil infiltration; however, lack of Ahr in myeloid cells did not protect LysM^Δ/ΔAhr mice from exacerbated APAP hepatotoxicity. In contrast, Alb^Δ/ΔAhr mice were largely protected from ITE-induced aggravated liver damage, indicating that Ahr activation in hepatocytes, but not in myeloid cells, was instrumental for disease exacerbation. Mechanistically, Ahr activation fueled hepatic accumulation of toxic APAP metabolites by up-regulating expression of the APAP-metabolizing enzyme Cyp1a2, a direct Ahr downstream target.

CONCLUSIONS

Ahr activation in hepatocytes potentiates APAP-induced hepatotoxicity. Thus, individual exposition to environmental Ahr ligands might explain individual sensitivity to hyperacute liver failure.

Controlled Functional Zonation of Hepatocytes In Vitro by Engineering of Wnt Signaling

Key liver functions, including protein synthesis, carbohydrate metabolism, and detoxification, are performed by specific populations of hepatocytes that are defined by their relative positions within the liver lobules. On a molecular level, the functional heterogeneity with periportal and pericentral phenotypes, so-called metabolic liver zonation, is mainly established by a gradient of canonical Wnt signaling activity. Since the relevant physiological cues are missing in in vitro liver models, they fail to reflect the functional heterogeneity and thus lack many liver functions. We synthetically re-engineered Wnt signaling in murine and human hepatocytes using a doxycycline-dependent cassette for externally controlled digital expression of stabilized β-catenin. Thereby, we achieved adjustable mosaic-like activation of Wnt signaling in in vitro-cultured hepatocytes that was resistant to negative-feedback loops. This allowed the establishment of long-term-stable periportal-like and pericentral-like phenotypes that mimic the heterogeneity observed in vivo. The in vitro-zonated hepatocytes show differential expression of drug-metabolizing enzymes and associated differential toxicity and higher levels of autophagy. Furthermore, recombinant adeno-associated virus and hepatitis C virus preferentially transduce the pericentral-like zonation phenotype, suggesting a bias of these viruses that has been unappreciated to date. These tightly controlled in vivo-like systems will be important for studies evaluating aspects of liver zonation and for the assessment of drug toxicity for mouse and man.

A mitochondrial journey through acetaminophen hepatotoxicity

Acetaminophen (APAP) overdose is the leading cause of acute liver failure in the United States and APAP-induced hepatotoxicity is initiated by formation of a reactive metabolite which depletes hepatic glutathione and forms protein adducts. Studies over the years have established the critical role of c-Jun N terminal kinase (JNK) and its mitochondrial translocation, as well as mitochondrial oxidant stress and subsequent induction of the mitochondrial permeability transition in APAP pathophysiology. However, it is now evident that mitochondrial responses to APAP overdose are more nuanced than appreciated earlier, with multiple levels of control, for example, to dose of APAP. In addition, mitochondrial dynamics, as well as the organelle's importance in recovery and regeneration after APAP-induced liver injury is also being recognized, which are exciting new areas with significant therapeutic potential. Thus, this review examines the temporal course of hepatocyte mitochondrial responses to an APAP overdose with an emphasis on mechanistic response to various trigger checkpoints such as NAPQI-mitochondrial protein adduct formation and activated JNK translocation. Mitochondrial dynamics, the organelle's role in recovery after APAP and emerging areas of research which promise to provide further insight into modulation of APAP pathophysiology by these fascinating organelles will also be discussed.

Cortex dictamni-induced liver injury in mice: The role of P450-mediated metabolic activation of furanoids

Many furan containing compounds have been reported to be toxic resulted from the metabolic activation of the furan ring to reactive metabolite (RM). Cortex Dictamni (CD), a widely used herbal medicine, has been reported to cause severe even fatal hepatotoxicity. The injurious components and mechanism of CD-induced liver injury remain unclear. Our preliminary study showed that dictamnine, one major furanoid in CD, caused mouse liver injury via its reactive epoxide metabolite. Besides dictamnine, the major components of CD are series of bioactivation-alerting furanoids. Thus, we hypothesize that series of furanoids in CD may undergo metabolic activation and play a key role in CD-induced liver injury. Here, a single oral dose of 60 g/kg ethanol extract of CD (ECD) caused severe hepatocellular necrosis in mice at 24 h post-dose. ECD-induced liver injury showed a dose- and time-dependent manner. The hepatotoxic effects could be completely abolished by P450 nonselective inhibitor 1-aminobenzotriazole (ABT) and strongly modulated by other P450 modulators. The furanoids-concentrated fraction of ECD was responsible for the hepatotoxicity. At least ten furanoids with high abundance in ECD, such as obakunone, dictamnine, fraxinellone, limonin, were found to be metabolized to reactive epoxide or cis-enedione. The RM levels were consistent with the liver injury degree. Multiple furanoids, rather than single one, cooperatively contributed to the hepatotoxicity. ECD-induced liver injury could be reproduced by a mixture of pure furanoids. In summary, this study provides toxic component profiles of CD and demonstrates that P450-mediated bioactivation of multiple furanoids is responsible for CD-induced liver injury.

Holland C, Zaza A, Murad W, Hassan R, A. Ahmed Y, Abbas T, A. Abdelrahim E, Schneider KM, Matz-Soja M, Reinders J, Gebhardt R, Berres ML, Hatting M, Drasdo D, Saez-Rodriguez J, Trautwein C, G. Hengstler J. Influence of Liver Fibrosis on Lobular Zonation

Little is known about how liver fibrosis influences lobular zonation. To address this question, we used three mouse models of liver fibrosis, repeated CCl4 administration for 2, 6 and 12 months to induce pericentral damage, as well as bile duct ligation (21 days) and mdr2-/- mice to study periportal fibrosis. Analyses were performed by RNA-sequencing, immunostaining of zonated proteins and image analysis. RNA-sequencing demonstrated a significant enrichment of pericentral genes among genes downregulated by CCl4; vice versa, periportal genes were enriched among the upregulated genes. Immunostaining showed an almost complete loss of pericentral proteins, such as cytochrome P450 enzymes and glutamine synthetase, while periportal proteins, such as arginase 1 and CPS1 became expressed also in pericentral hepatocytes. This pattern of fibrosis-associated 'periportalization' was consistently observed in all three mouse models and led to complete resistance to hepatotoxic doses of acetaminophen (200 mg/kg). Characterization of the expression response identified the inflammatory pathways TGFβ, NFκB, TNFα, and transcription factors NFKb1, Stat1, Hif1a, Trp53, and Atf1 among those activated, while estrogen-associated pathways, Hnf4a and Hnf1a, were decreased. In conclusion, liver fibrosis leads to strong alterations of lobular zonation, where the pericentral region adopts periportal features. Beside adverse consequences, periportalization supports adaptation to repeated doses of hepatotoxic compounds.

Propagation of Pericentral Necrosis During Acetaminophen-Induced Liver Injury: Evidence for Early Interhepatocyte Communication and Information Exchange

Acetaminophen (APAP)-induced liver injury is clinically significant, and APAP overdose in mice often serves as a model for drug-induced liver injury in humans. By specifying that APAP metabolism, reactive metabolite formation, glutathione depletion, and mitigation of mitochondrial damage within individual hepatocytes are functions of intralobular location, an earlier virtual model mechanism provided the first concrete multiattribute explanation for how and why early necrosis occurs close to the central vein (CV). However, two characteristic features could not be simulated consistently: necrosis occurring first adjacent to the CV, and subsequent necrosis occurring primarily adjacent to hepatocytes that have already initiated necrosis. We sought parsimonious model mechanism enhancements that would manage spatiotemporal heterogeneity sufficiently to enable meeting two new target attributes and conducted virtual experiments to explore different ideas for model mechanism improvement at intrahepatocyte and multihepatocyte levels. For the latter, evidence supports intercellular communication via exosomes, gap junctions, and connexin hemichannels playing essential roles in the toxic effects of chemicals, including facilitating or counteracting cell death processes. Logic requiring hepatocytes to obtain current information about whether downstream and lateral neighbors have triggered necrosis enabled virtual hepatocytes to achieve both new target attributes. A virtual hepatocyte that is glutathione-depleted uses that information to determine if it will initiate necrosis. When a less-stressed hepatocyte is flanked by at least two neighbors that have triggered necrosis, it too will initiate necrosis. We hypothesize that the resulting intercellular communication-enabled model mechanism is analogous to the actual explanation for APAP-induced hepatotoxicity at comparable levels of granularity.