Bile duct ligation (BDL) is an experimental procedure that mimics obstructive cholestatic disease. One of the early consequences of BDL in rodents is the appearance of so-called bile infarcts that correspond to Charcot-Gombault necrosis in human cholestasis. The mechanisms causing bile infarcts and their pathophysiological relevance are unclear. Therefore, intravital two photon-based imaging of BDL mice was performed with fluorescent bile salts (BS) and non-BS organic anion analogues. Key findings were followed up by matrix-assisted laser desorption ionization imaging, clinical chemistry, immunostaining, and gene expression analyses. In the acute phase, 1-3 days after BDL, BS concentrations in bile increased and single-cell bile microinfarcts occurred in dispersed hepatocytes throughout the liver caused by the rupture of the apical hepatocyte membrane. This rupture occurred after loss of mitochondrial membrane potential, followed by entry of bile, cell death, and a "domino effect" of further death events of neighboring hepatocytes. Bile infarcts provided a trans-epithelial shunt between bile canaliculi and sinusoids by which bile constituents leaked into blood. In the chronic phase, ≥21 days after BDL, uptake of BS tracers at the sinusoidal hepatocyte membrane was reduced. This contributes to elevated concentrations of BS in blood and decreased concentrations in the biliary tract. Conclusion: Bile microinfarcts occur in the acute phase after BDL in a limited number of dispersed hepatocytes followed by larger infarcts involving neighboring hepatocytes, and they allow leakage of bile from the BS-overloaded biliary tract into blood, thereby protecting the liver from BS toxicity; in the chronic phase after BDL, reduced sinusoidal BS uptake is a dominant protective factor, and the kidney contributes to the elimination of BS until cholemic nephropathy sets in.

Atorvastatin δ-lactone, a major, pharmacologically inactive metabolite, has been associated with toxicity. In a previous study we showed that polymorphisms of UGT1A3 influence atorvastatin δ-lactone formation. Here we investigated the reverse reaction, atorvastatin δ-lactone hydrolysis, in a human liver bank. Screening of microarray data revealed paraoxonases PON1 and PON3 among 17 candidate esterases. Microsomal δ-lactone hydrolysis was significantly correlated to PON1 and PON3 protein (r(s) = 0.60; r(s) = 0.62, respectively; P < 0.0001). PON1 and PON3 were strongly correlated to each other (r(s) = 0.60) but PON1 was shown to be more extensively glycosylated than PON3. In addition a novel splice-variant of PON3 was identified. Genotyping of 40 polymorphisms within the PON-locus identified PON1 promoter polymorphisms (-108T > C, -832G > A, -1741G > A) and a tightly linked group of PON3 polymorphisms (-4984A > G, -4105G > A, -1091A > G, -746C > T, and F21F) to be associated with changes in atorvastatin δ-lactone hydrolysis and expression of PON1 but not PON3. However, carriers of the common PON1 polymorphisms L55M or Q192R showed no difference in δ-lactone hydrolysis or PON expression. Haplotype analysis revealed decreased δ-lactone hydrolysis in carriers of the most common haplotype *1 compared to carriers of haplotypes *2, *3, *4, and *7. Analysis of non-genetic factors showed association of hepatocellular and cholangiocellular carcinoma with decreased PON1 and PON3 expression, respectively. Increased C-reactive protein and γ-glutamyl transferase levels were associated with decreased protein expression of both enzymes, and increased bilirubin levels, cholestasis, and presurgical exposure to omeprazole or pantoprazole were related to decreased PON3 protein. In conclusion, PON-locus polymorphisms affect PON1 expression whereas non-genetic factors have an effect on PON1 and PON3 expression. This may influence response to therapy or adverse events in statin treatment.

Efflux transporters like MDR1 and MRP2 may modulate the pharmacokinetics of about 50 % of all drugs. It is currently unknown how much of the variation in the activities of important drug membrane transporters like MDR1 or MRP2 is determined by genetic or by environmental factors. In this study we assessed the heritability of the pharmacokinetics of talinolol as a putative probe drug for MDR1 and possibly other membrane transporters.

The Hedgehog (Hh) and Wnt/β-Catenin (Wnt) cascades are morphogen pathways whose pronounced influence on adult liver metabolism has been identified in recent years. How both pathways communicate and control liver metabolic functions are largely unknown. Detecting core components of Wnt and Hh signaling and mathematical modeling showed that both pathways in healthy liver act largely complementary to each other in the pericentral (Wnt) and the periportal zone (Hh) and communicate mainly by mutual repression. The Wnt/Hh module inversely controls the spatiotemporal operation of various liver metabolic pathways, as revealed by transcriptome, proteome, and metabolome analyses. Shifting the balance to Wnt (activation) or Hh (inhibition) causes pericentralization and periportalization of liver functions, respectively. Thus, homeostasis of the Wnt/Hh module is essential for maintaining proper liver metabolism and to avoid the development of certain metabolic diseases. With caution due to minor species-specific differences, these conclusions may hold for human liver as well.

Women are more prone to develop either hypothyroidism or cholesterol gallstones than men. However, a male predominance in cholesterol gallstones under hypothyroidism was reported. Recently, a novel pathogenic link between thyroid hormone (TH) deficiency and cholesterol gallstones has been described in male mice. Here, we investigate if TH deficiency impacts cholesterol gallstone formation in females by the same mechanism. Three-month-old C57BL/6J mice were randomly divided into a control, a TH deficient, a lithogenic, and a lithogenic + TH deficient group and diet-treated for two, four, and six weeks. Gallstone prevalence, liver function tests, bile composition, hepatic gene expression, and gallbladder aquaporin expression and localization were investigated. Cholesterol gallstones were observed in lithogenic + TH deficient but not lithogenic only female mice. Diminished hydrophilicity of primary bile acids due to decreased gene expression of hepatic detoxification phase II enzymes was observed. A sex-specific expression and localization of hepatobiliary aquaporins involved in transcellular water and glycerol permeability was observed under TH deficient and lithogenic conditions. TH deficiency promotes cholesterol gallstone formation in female C57BL/6J mice by the same mechanism as observed in males. However, cholesterol gallstone prevalence was lower in female than male C57BL/6J mice. Interestingly, the sex-specific expression and localization of hepatobiliary aquaporins could protect female C57BL/6J mice to cholestasis and could reduce biliary water transport in male C57BL/6J mice possibly contributing to the sex-dependent cholesterol gallstone prevalence under TH deficiency.

The extracellular signal-regulated kinase (ERK) pathway regulates cell growth, and is hyper-activated and associated with drug resistance in hepatocellular carcinoma (HCC). Metabolic pathways are profoundly dysregulated in HCC. Whether an altered metabolic state is linked to activated ERK pathway and drug response in HCC is unaddressed.

The Hedgehog signaling pathway regulates many processes during embryogenesis and the homeostasis of adult organs. Recent data suggest that central metabolic processes and signaling cascades in the liver are controlled by the Hedgehog pathway and that changes in hepatic Hedgehog activity also affect peripheral tissues, such as the reproductive organs in females. Here, we show that hepatocyte-specific deletion of the Hedgehog pathway is associated with the dramatic expansion of adipose tissue in mice, the overall phenotype of which does not correspond to the classical outcome of insulin resistance-associated diabetes type 2 obesity. Rather, we show that alterations in the Hedgehog signaling pathway in the liver lead to a metabolic phenotype that is resembling metabolically healthy obesity. Mechanistically, we identified an indirect influence on the hepatic secretion of the fibroblast growth factor 21, which is regulated by a series of signaling cascades that are directly transcriptionally linked to the activity of the Hedgehog transcription factor GLI1. The results of this study impressively show that the metabolic balance of the entire organism is maintained via the activity of morphogenic signaling pathways, such as the Hedgehog cascade. Obviously, several pathways are orchestrated to facilitate liver metabolic status to peripheral organs, such as adipose tissue.

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.

Introduction: Bulevirtide is a first-in-class antiviral drug to treat chronic hepatitis B/D. We investigated the drug-drug interaction potential and pharmacokinetics of high-dose subcutaneous bulevirtide (5 mg twice daily) with organic anion transporting polypeptide 1B1 (OATP1B1) and cytochrome P450 (CYP) 3A4. Methods: This was a single-center, open-label, fixed-sequence drug-drug interaction trial in 19 healthy volunteers. Before and at bulevirtide steady state, participants ingested a single 40 mg dose of pravastatin. A midazolam microdose was applied to quantify CYP3A4 activity. Results: At bulevirtide steady state, pravastatin area under the concentration-time curve (AUC0-∞) increased 1.32-fold (90% CI 1.08-1.61). The 5 mg bulevirtide twice-daily treatment resulted in a mean AUC0-12 of 1210 h*ng/ml (95% CI 1040-1408) and remained essentially unchanged under the influence of pravastatin. CYP3A4 activity did not change to a clinically relevant extent. As expected, total bile acids increased substantially (35-fold) compared to baseline during bulevirtide treatment. All study medication was well tolerated. Discussion: The study demonstrated that high-dose bulevirtide inhibited OATP1B-mediated hepatic uptake of the marker substrate pravastatin but the extent is considered clinically not relevant. Changes in CYP3A4 activity were also not clinically relevant. In conclusion, this study suggests that OATP1B substrate drugs as well as CYP3A4 substrates may safely be used without dose adjustment in patients treated with bulevirtide. However, in patients using high statin doses and where concomitant factors potentially further increase statin exposure, caution may be required when using bulevirtide.

In the progress of allergic and irritant contact dermatitis, chemicals that cause the generation of reactive oxygen species trigger a heat shock response in keratinocytes. In this study, an optical sensor cell line based on cultured human keratinocytes (HaCaT cells) expressing green fluorescent protein (GFP) under the control of the stress-inducible HSP70B' promoter were constructed. Exposure of HaCaT sensor cells to 25 µM cadmium, a model substance for oxidative stress induction, provoked a 1.7-fold increase in total glutathione and a ~300-fold induction of transcript level of the gene coding for heat shock protein HSP70B'. An extract of Arnica montana flowers resulted in a strong induction of the HSP70B' gene and a pronounced decrease of total glutathione in keratinocytes. The HSP70B' promoter-based sensor cells conveniently detected cadmium-induced stress using GFP fluorescence as read-out with a limit of detection of 6 µM cadmium. In addition the sensor cells responded to exposure of cells to A. montana extract with induction of GFP fluorescence. Thus, the HaCaT sensor cells provide a means for the automated detection of the compromised redox status of keratinocytes as an early indicator of the development of human skin disorders and could be applied for the prediction of skin irritation in more complex in vitro 3D human skin models and in the development of micro-total analysis systems (µTAS) that may be utilized in dermatology, toxicology, pharmacology and drug screenings.

PGRMC1 (progesterone receptor membrane component 1) is a highly conserved heme binding protein, which is overexpressed especially in hormone receptor-positive breast cancer and plays an important role in breast carcinogenesis. Nevertheless, little is known about the mechanisms by which PGRMC1 drives tumor progression. The aim of our study was to investigate the involvement of PGRMC1 in cholesterol metabolism to detect new mechanisms by which PGRMC1 can increase lipid metabolism and alter cancer-related signaling pathways leading to breast cancer progression.

We recently reported that miR-146a is differentially expressed in ALK+ and ALK- anaplastic large cell lymphoma (ALCL). In this study, the downstream targets of miR-146a in ALK+ ALCL were investigated by transcriptome analysis, identifying CD147 as potential target gene. Because CD147 is differentially expressed in ALK+ ALCL versus ALK- ALCL and normal T cells, this gene emerged as a strong candidate for the pathogenesis of this tumor. Here we demonstrate that CD147 is a direct target of miR-146 and contributes to the survival and proliferation of ALK+ ALCL cells in vitro and to the engraftment and tumor growth in vivo in an ALK+ ALCL-xenotransplant mouse model. CD147 knockdown in ALK+ ALCL cells resulted in loss of monocarboxylate transporter 1 (MCT1) expression, reduced glucose consumption and tumor growth retardation, as demonstrated by [18F]FDG-PET/MRI analysis. Investigation of metabolism in vitro and in vivo supported these findings, revealing reduced aerobic glycolysis and increased basal respiration in CD147 knockdown. In conclusion, our findings indicate that CD147 is of vital importance for ALK+ ALCL to maintain the high energy demand of rapid cell proliferation, promoting lactate export, and tumor growth. Furthermore, CD147 has the potential to serve as a novel therapeutic target in ALK+ ALCL, and warrants further investigation.

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants produced by incomplete combustion of organic matter. They induce their own metabolism by upregulating xenobiotic-metabolizing enzymes such as cytochrome P450 monooxygenase 1A1 (CYP1A1) by activating the aryl hydrocarbon receptor (AHR). However, previous studies showed that individual PAHs may also interact with the constitutive androstane receptor (CAR). Here, we studied ten PAHs, different in carcinogenicity classification, for their potential to activate AHR- and CAR-dependent luciferase reporter genes in human liver cells. The majority of investigated PAHs activated AHR, while non-carcinogenic PAHs tended to activate CAR. We further characterized gene expression, protein abundancies and activities of the AHR targets CYP1A1 and 1A2, and the CAR target CYP2B6 in human HepaRG hepatoma cells. Enzyme induction patterns strongly resembled the profiles obtained at the receptor level, with AHR-activating PAHs inducing CYP1A1/1A2 and CAR-activating PAHs inducing CYP2B6. In summary, this study provides evidence that beside well-known activation of AHR, some PAHs also activate CAR, followed by subsequent expression of respective target genes. Furthermore, we found that an increased PAH ring number is associated with AHR activation as well as the induction of DNA double-strand breaks, whereas smaller PAHs activated CAR but showed no DNA-damaging potential.

HepaRG cells are increasingly accepted as model for human drug metabolism and other hepatic functions. We used lentiviral transduction of undifferentiated HepaRG cells to deliver Cas9 and two alternative sgRNAs targeted at NADPH:cytochrome P450 oxidoreductase (POR), the obligate electron donor for microsomal cytochromes P450 (CYP). Cas9-expressing HepaRGVC (vector control) cells were phenotypically similar to wild type HepaRG cells and could be differentiated into hepatocyte-like cells by DMSO. Genetic POR-knockout resulted in phenotypic POR knockdown of up to 90% at mRNA, protein, and activity levels. LC-MS/MS measurement of seven CYP-activities showed differential effects of POR-knockdown with CYP2C8 being least and CYP2C9 being most affected. Further studies on cytochrome b5 (CYB5), an alternative NADH-dependent electron donor indicated particularly strong support of CYP2C8-dependent amodiaquine N-deethylation by CYB5 and this was confirmed by genetic CYB5 single- and POR/CYB5 double-knockout. POR-knockdown also affected CYP expression on mRNA and protein level, with CYP1A2 being induced severalfold, while CYP2C9 was strongly downregulated. In summary our results show that POR/NADPH- and CYB5/NADH-electron transport systems influence human drug metabolizing CYPs differentially and differently than mouse Cyps. Our Cas9-expressing HepaRGVC cells should be suitable to study the influence of diverse genes on drug metabolism and other hepatic functions.

Acetaminophen (APAP) overdose remains a frequent cause of acute liver failure, which is generally accompanied by increased levels of serum bile acids (BAs). However, the pathophysiological role of BAs remains elusive. Herein, we investigated the role of BAs in APAP-induced hepatotoxicity.

Okadaic acid (OA) is an algae-produced lipophilic marine biotoxin that accumulates in the fatty tissue of filter-feeding shellfish. Ingestion of contaminated shellfish leads to the diarrheic shellfish poisoning syndrome. Furthermore, several other effects of OA like genotoxicity, liver toxicity and tumor-promoting properties have been observed, probably linked to the phosphatase-inhibiting properties of the toxin. It has been shown that at high doses OA can disrupt the physical barrier of the intestinal epithelium. As the intestine and the liver do not only constitute a physical, but also a metabolic barrier against xenobiotic exposure, we here investigated the impact of OA on the expression of cytochrome P450 (CYP) enzymes and transporter proteins in human HepaRG cells liver cells in vitro at non-cytotoxic concentrations. The interplay of OA with known CYP inducers was also studied. Data show that the expression of various xenobiotic-metabolizing CYPs was downregulated after exposure to OA. Moreover, OA was able to counteract the activation of CYPs by their inducers. A number of transporters were also mainly downregulated. Overall, we demonstrate that OA has a significant effect on xenobiotic metabolism barrier in liver cells, highlighting the possibility for interactions of OA exposure with the metabolism of drugs and xenobiotics.

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 CCl4 and by a Western diet, respectively. Genes consistently up-regulated between the chronic CCl4 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/.

Current therapies for metastatic clear cell renal cell carcinoma (ccRCC) show limited efficacy. Drug efficacy, typically investigated in preclinical cell line models during drug development, is influenced by pharmacogenes involved in targeting and disposition of drugs. Here we show through genome-wide DNA methylation profiling, that methylation patterns are concordant between primary ccRCC and macro-metastases irrespective of metastatic sites (rs ≥ 0.92). However, 195,038 (41%) of all investigated CpG sites, including sites within pharmacogenes, were differentially methylated (adjusted P < 0.05) in five established RCC cell lines compared to primary tumors, resulting in altered transcriptional expression. Exemplarily, gene-specific analyses of DNA methylation, mRNA and protein expression demonstrate lack of expression of the clinically important drug transporter OCT2 (encoded by SLC22A2) in cell lines due to hypermethylation compared to tumors or metastases. Our findings provide evidence that RCC cell lines are of limited benefit for prediction of drug effects due to epigenetic alterations. Similar epigenetic landscape of ccRCC-metastases and tumors opens new avenue for future therapeutic strategies.

Organic cation transporters (OCTs) determine not only physiological processes but are also involved in the cellular uptake of anticancer agents. Based on microarray analyses in hepatocellular carcinoma (HCC), SLC22A1/OCT1 mRNA seems to be downregulated, but systematic protein expression data are currently missing. Moreover, the underlying molecular mechanisms responsible for altered SLC22A1 expression in HCC are not fully understood. Therefore, we investigated the role of DNA methylation in the transcriptional regulation of the family members SLC22A1/OCT1, SLC22A2/OCT2 and SLC22A3/OCT3 in HCC.

Metformin, an oral insulin-sensitizing drug, is actively transported into cells by organic cation transporters (OCT) 1, 2, and 3 (encoded by SLC22A1, SLC22A2, or SLC22A3), which are tissue specifically expressed at significant levels in various organs such as liver, muscle, and kidney. Because metformin does not undergo hepatic metabolism, drug-drug interaction by inhibition of OCT transporters may be important. So far, comprehensive data on the interaction of proton pump inhibitors (PPIs) with OCTs are missing although PPIs are frequently used in metformin-treated patients. Using in silico modeling and computational analyses, we derived pharmacophore models indicating that PPIs (i.e. omeprazole, pantoprazole, lansoprazole, rabeprazole, and tenatoprazole) are potent OCT inhibitors. We then established stably transfected cell lines expressing the human uptake transporters OCT1, OCT2, or OCT3 and tested whether these PPIs inhibit OCT-mediated metformin uptake in vitro. All tested PPIs significantly inhibited metformin uptake by OCT1, OCT2, and OCT3 in a concentration-dependent manner. Half-maximal inhibitory concentration values (IC(50)) were in the low micromolar range (3-36 µM) and thereby in the range of IC(50) values of other potent OCT drug inhibitors. Finally, we tested whether the PPIs are also transported by OCTs, but did not identify PPIs as OCT substrates. In conclusion, PPIs are potent inhibitors of the OCT-mediated metformin transport in vitro. Further studies are needed to elucidate the clinical relevance of this drug-drug interaction with potential consequences on metformin disposition and/or efficacy.