Structure-activity relationship in the passage of different pyrrolizidine alkaloids through the gastrointestinal barrier: ABCB1 excretes heliotrine and echimidine

Risk assessment of short-term intake of pyrrolizidine alkaloids in food: derivation of an acute reference dose

Pyrrolizidine alkaloids (PA) are phytochemicals that are known to act as human hepatotoxins and are also considered to be genotoxic carcinogens. Several plant-derived foods are frequently contaminated with PA, like teas and herbal infusions, spices and herbs or certain food supplements. With respect to the chronic toxicity of PA, the carcinogenic potential of PA is generally regarded as the critical toxicological effect. The risk assessment of the short-term toxicity of PA, however, is internationally less consistent. The characteristic pathological syndrome of acute PA toxicity is hepatic veno-occlusive disease. High PA exposure levels may lead to liver failure and even death as documented by several case reports. In the present report, we suggest a risk assessment approach for the derivation of an acute reference dose (ARfD) for PA of 1 µg/kg body weight per day based on a sub-acute animal toxicity study in rats after oral PA administration. The derived ARfD value is further supported by several case reports describing acute human poisoning following accidental PA intake. The here derived ARfD value may be used for PA risk assessment in cases where the short-term toxicity of PA is of interest in addition to the assessment of the long-term risks.

Structure-Dependent Toxicokinetics of Selected Pyrrolizidine Alkaloids In Vitro

Phytochemicals like pyrrolizidine alkaloids (PAs) can affect the health of humans and animals. PAs can occur for example in tea, honey or herbs. Some PAs are known to be cytotoxic, genotoxic, and carcinogenic. Upon intake of high amounts, hepatotoxic and pneumotoxic effects were observed in humans. This study aims to elucidate different toxicokinetic parameters like the uptake of PAs and their metabolism with in vitro models. We examined the transport rates of differently structured PAs (monoester, open-chained diester, cyclic diester) over a model of the intestinal barrier. After passing the intestinal barrier, PAs reach the liver, where they are metabolized into partially instable electrophilic metabolites interacting with nucleophilic centers. We investigated this process by the usage of human liver, intestinal, and lung microsomal preparations for incubation with different PAs. These results are completed with the detection of apoptosis as indicator for bioactivation of the PAs. Our results show a structure-dependent passage of PAs over the intestinal barrier. PAs are structure-dependently metabolized by liver microsomes and, to a smaller extent, by lung microsomes. The detection of apoptosis of A549 cells treated with lasiocarpine and monocrotaline following bioactivation by human liver or lung microsomes underlines this result. Conclusively, our results help to shape the picture of PA toxicokinetics which could further improve the knowledge of molecular processes leading to observed effects of PAs in vivo.

Organic Cation Transporter I and Na+ /taurocholate Co-Transporting Polypeptide are Involved in Retrorsine- and Senecionine-Induced Hepatotoxicity in HepaRG cells

SCOPE

1,2-unsaturated pyrrolizidine alkaloids (PAs) are secondary plant metabolites that are found in many plant species throughout the world. They are of concern for risk assessment as consumption of contaminated foodstuff can cause severe liver damage. Of late, transporter-mediated uptake and transport has advanced as a vital determinant of PA toxicity. In this study, the authors investigate a transporter-mediated uptake of PAs and its implications in PA toxicity.

METHODS AND RESULTS

We show that transporter expression levels are significantly affected by treatment with the PAs senecionine (Sc) and retrorsine (Re) in the human hepatoma cell line HepaRG. Furthermore, the specific contribution to PA uptake of the two transporters Na⁺ /taurocholate co-transporting polypeptide (SLC10A1) and organic cation transporter I (SLC22A1), both belonging to the heterogeneous solute carrier super family, is investigated by means of a siRNA-mediated knockdown approach. Knockdown of both uptake transporters result in reduced uptake of Re and Sc in a time-dependent manner and attenuated PA-mediated cytotoxic effects in HepaRG cells.

CONCLUSION

Our results confirm previous findings of active transport mechanisms of PAs into hepatocytes and highlight the importance of toxicokinetic studies for the risk assessment of PAs.

Active Transport of Hepatotoxic Pyrrolizidine Alkaloids in HepaRG Cells

1,2-unsaturated pyrrolizidine alkaloids (PAs) are secondary plant metabolites occurring as food contaminants that can cause severe liver damage upon metabolic activation in hepatocytes. However, it is yet unknown how these contaminants enter the cells. The role of hepatic transporters is only at the beginning of being recognized as a key determinant of PA toxicity. Therefore, this study concentrated on assessing the general mode of action of PA transport in the human hepatoma cell line HepaRG using seven structurally different PAs. Furthermore, several hepatic uptake and efflux transporters were targeted with pharmacological inhibitors to identify their role in the uptake of the PAs retrorsine and senecionine and in the disposition of their N-oxides (PANO). For this purpose, PA and PANO content was measured in the supernatant using LC-MS/MS. Also, PA-mediated cytotoxicity was analyzed after transport inhibition. It was found that PAs are taken up into HepaRG cells in a predominantly active and structure-dependent manner. This pattern correlates with other experimental endpoints such as cytotoxicity. Pharmacological inhibition of the influx transporters Na⁺/taurocholate co-transporting polypeptide (SLC10A1) and organic cation transporter 1 (SLC22A1) led to a reduced uptake of retrorsine and senecionine into HepaRG cells, emphasizing the relevance of these transporters for PA toxicokinetics.

Metabolic Pattern of Hepatotoxic Pyrrolizidine Alkaloids in Liver Cells

Contamination with 1,2-unsaturated pyrrolizidine alkaloids (PAs) is a serious problem for certain phytomedicines, foods, and animal feeds. Several of these PAs are genotoxic and carcinogenic, primarily in the liver, upon cytochrome P450 (CYP)-catalyzed activation into reactive (pyrrolic and pyrrole-like) metabolites. Here we investigated the metabolism of selected PAs (echimidine, europine, lasiocarpine, lycopsamine, retrorsine, and senecionine) in rat hepatocytes in primary culture and in human CYP3A4-transfected HepG2 cells. The open-chained diesters echimidine and lasiocarpine and the cyclic diester senecionine were extensively metabolized in rat hepatocytes into a broad spectrum of products released into the medium. A large portion of unidentified, possibly irreversibly bound, products remained in the cells while detectable amounts of reactive and other metabolites were found in the incubation media. In HepG2-CYP3A4 cells, lasiocarpine was more extensively metabolized than echimidine and senecionine which also gave rise to the release of pyrrolic metabolites. In human cells, no pyrrolic metabolites were detected in retrorsine or lycopsamine incubations, while no such metabolites were detected from europine in both cell types. Other types of metabolic changes comprised modifications such as side chain demethylation or oxygenation reactions like the formation of N-oxides. The latter, considered as a detoxification step, was a major pathway with cyclic diesters, was less distinctive for echimidine and lycopsamine and almost negligible for lasiocarpine and europine. Our data are in agreement with previously published cyto- and genotoxicity findings and suggests that the metabolic pattern may contribute substantially to the specific toxic potency of a certain congener. In addition, marked differences were found for certain congeners between rat hepatocytes and transfected human HepG2 cells, whereby a high level of bioactivation was found for lasiocarpine, whereas a very low level of bioactivation was observed for monoesters, in particular in human cells.

Toxic Prediction of Pyrrolizidine Alkaloids and Structure-Dependent Induction of Apoptosis in HepaRG Cells

Pyrrolizidine alkaloids (PAs) are common phytotoxins and could cause liver genotoxicity/carcinogenicity following metabolic activation. However, the toxicity of different structures remains unclear due to the wide variety of PAs. In this study, the absorption, distribution, metabolism, excretion, and toxicity (ADMET) of 40 PAs were analyzed, and their toxicity was predicted by Komputer Assisted Technology (TOPKAT) using Discovery Studio software. The in silico results showed that all PAs except retronecine had good intestinal absorption, and all PAs were predicted to have different toxicity ranges. To verify the predictive results, 4 PAs were selected to investigate cell injury and possible mechanisms of the differentiation in HepaRG cells, including retronecine type of twelve-membered cyclic diester (retrorsine), eleven-membered cyclic diester (monocrotaline), noncyclic diester (retronecine), and platynecine type (platyphylline). After 24 h exposure, retronecine-type PAs exhibited concentration-dependent cytotoxicity. The high-content screening assay showed that cell oxidative stress, mitochondrial damage, endoplasmic reticulum stress, and the concentration of calcium ions increased, and neutral lipid metabolism was changed notably in HepaRG cells. Induced apoptosis by PAs was indicated by cell cycle arrest in the G2/M phase, disrupting the mitochondrial membrane potential. Overall, our study revealed structure-dependent cytotoxicity and apoptosis after PA exposure, suggesting that the prediction results of in silico have certain reference values for compound toxicity. A 1,2-membered cyclic diester seems to be a more potent apoptosis inducer than other PAs.

In vitro biotransformation of pyrrolizidine alkaloids in different species: part II-identification and quantitative assessment of the metabolite profile of six structurally different pyrrolizidine alkaloids

Pyrrolizidine alkaloids (PA) exert their toxic effects only after bioactivation. Although their toxicity has already been studied and metabolic pathways including important metabolites were described, the quantification of the latter revealed a large unknown portion of the metabolized PA. In this study, the qualitative and quantitative metabolite profiles of structurally different PAs in rat and human liver microsomes were investigated. Between five metabolites for europine and up to 48 metabolites for lasiocarpine were detected. Proposals for the chemical structure of each metabolite were derived based on fragmentation patterns using high-resolution mass spectrometry. The metabolite profiles of the diester PAs showed a relatively good agreement between both species. The metabolic reactions were summarized into three groups: dehydrogenation, oxygenation, and shortening of necic acid(s). While dehydrogenation of the necine base is considered as bioactivation, both other routes are considered as detoxification steps. The most abundant changes found for open chained diesters were dealkylations, while the major metabolic pathway for cyclic diesters was oxygenation especially at the nitrogen atom. In addition, all diester PAs formed several dehydrogenation products, via the insertion of a second double bond in the necine base, including the formation of glutathione conjugates. In rat liver microsomes, all investigated PAs formed dehydropyrrolizidine metabolites with the highest amount formed by lasiocarpine, whereas in human liver microsomes, these metabolites could only be detected for diesters. Our findings demonstrate that an extensive analysis of PA metabolism can provide the basis for a better understanding of PA toxicity and support future risk assessment.

Hepatotoxic pyrrolizidine alkaloids induce DNA damage response in rat liver in a 28-day feeding study

Pyrrolizidine alkaloids (PA) are secondary plant metabolites that occur as food and feed contaminants. Acute and subacute PA poisoning can lead to severe liver damage in humans and animals, comprising liver pain, hepatomegaly and the development of ascites due to occlusion of the hepatic sinusoids (veno-occlusive disease). Chronic exposure to low levels of PA can induce liver cirrhosis and liver cancer. However, it is not well understood which transcriptional changes are induced by PA and whether all hepatotoxic PA, regardless of their structure, induce similar responses. Therefore, a 28-day subacute rat feeding study was performed with six structurally different PA heliotrine, echimidine, lasiocarpine, senecionine, senkirkine, and platyphylline, administered at not acutely toxic doses from 0.1 to 3.3 mg/kg body weight. This dose range is relevant for humans, since consumption of contaminated tea may result in doses of ~ 8 µg/kg in adults and cases of PA ingestion by contaminated food was reported for infants with doses up to 3 mg/kg body weight. ALT and AST were not increased in all treatment groups. Whole-genome microarray analyses revealed pronounced effects on gene expression in the high-dose treatment groups resulting in a set of 36 commonly regulated genes. However, platyphylline, the only 1,2-saturated and, therefore, presumably non-hepatotoxic PA, did not induce significant expression changes. Biological functions identified to be affected by high-dose treatments (3.3 mg/kg body weight) comprise cell-cycle regulation associated with DNA damage response. These functions were found to be affected by all analyzed 1,2-unsaturated PA.

In conclusion, 1,2-unsaturated hepatotoxic PA induced cell cycle regulation processes associated with DNA damage response. Similar effects were observed for all hepatotoxic PA. Effects were observed in a dose range inducing no histopathological alterations and no increase in liver enzymes. Therefore, transcriptomics studies identified changes in expression of genes known to be involved in response to genotoxic compounds at PA doses relevant to humans under worst case exposure scenarios.

Structure-dependent hepato-cytotoxic potencies of selected pyrrolizidine alkaloids in primary rat hepatocyte culture

Contamination of food, feed and herbal medicines with plants containing pyrrolizidine alkaloids (PA) leads to measurable amounts of PA in many products. Since a number of PA are hepatotoxic in humans and animals and hepato-carcinogenic in animal experiments, the assessment of the relative toxic potencies of widely occurring PA contaminants warrants detailed investigation. Here, we studied the hepato-cytotoxic potencies of a number of relevant PA congeners in rat hepatocytes in primary culture. It was found that cyclic and open di-esters were much more toxic than mono-esters. Furthermore, the hepatocellular levels of cytochrome P450-catalyzed 7-benzoxyresorufin O-dealkylase (BROD) activity decreasing over time in culture, played an important role for activation of PA into cytotoxic metabolites. With a highly toxic PA (lasiocarpine), inhibition of BROD activity with ketoconazole markedly reduced toxicity while this was not obvious with the less toxic congener lycopsamine. Depletion of cellular glutathione with buthionine sulfoximine had no significant influence on the effects of highly toxic PA whereas it slightly increased toxicity of less potent congeners. Overall, our data partially confirm previously published structure-dependent interim Relative Potency (iREP) factors although for echimidine and monocrotaline in particular, substantial deviations were found, possibly due to specific toxicokinetic properties of these congeners.

Development of a sensitive analytical method for determining 44 pyrrolizidine alkaloids in teas and herbal teas via LC-ESI-MS/MS

Pyrrolizidine alkaloids (PA) and PA-N-oxides (PANO) are a large group of secondary plant metabolites comprising more than 660 compounds. Exhibiting geno- and hepatotoxic properties, they are responsible for multiple cases of food and feed poisoning over the last 100 years. For food and feed safety reasons, relevant PA/PANO should be monitored extensively in the main sources of PA/PANO intake. In this study, a sensitive analytical method was developed for detecting a broad range of 44 commercially available PA/PANO compounds, and in-house validation procedures were performed for several (herbal) teas. Various extraction solvents and procedures, as well as solid phase extraction materials for sample clean-up and analyte concentration, were tested to establish the methods' efficiency and effectiveness. Chromatographic conditions were optimised to obtain the best possible separation of isomers for the 44 PA/PANO analytes. The final method was proven very sensitive and accurate, with detection limits ranging from 0.1 to 7.0 μg/kg and precisions between 0.7 and 16.1%. For 40 of the analytes, the recovery rates ranged from 60.7 to 128.8%. The applicability and trueness of the method were examined by analysing tea samples from a local supermarket and comparing them to a reference material. At least one PA/PANO analyte was detected in 17 of the 18 samples under investigation, and the sum contents of the samples ranged from 0.1 to 47.9 μg/kg. Knowledge of the PA/PANO composition in a sample can be used to indicate the botanical origin of the impurity and, thus, the geographical region of cultivation.

Prediction of in vivo genotoxicity of lasiocarpine and riddelliine in rat liver using a combined in vitro-physiologically based kinetic modelling-facilitated reverse dosimetry approach

Pyrrolizidine alkaloids (PAs) are naturally occurring genotoxic compounds, and PA-containing plants can pose a risk to humans through contaminated food sources and herbal products. Upon metabolic activation, PAs can form DNA adducts, DNA and protein cross links, chromosomal aberrations, micronuclei, and DNA double-strand breaks. These genotoxic effects may induce gene mutations and play a role in the carcinogenesis of PAs. This study aims to predict in vivo genotoxicity for two well-studied PAs, lasiocarpine and riddelliine, in rat using in vitro genotoxicity data and physiologically based kinetic (PBK) modelling-based reverse dosimetry. The phosphorylation of histone protein H2AX was used as a quantitative surrogate endpoint for in vitro genotoxicity of lasiocarpine and riddelliine in primary rat hepatocytes and human HepaRG cells. The in vitro concentration-response curves obtained from primary rat hepatocytes were subsequently converted to in vivo dose-response curves from which points of departure (PoDs) were derived that were compared to available in vivo genotoxicity data. The results showed that the predicted PoDs for lasiocarpine and riddelliine were comparable to in vivo genotoxicity data. It is concluded that this quantitative in vitro-in silico approach provides a method to predict in vivo genotoxicity for the large number of PAs for which in vivo genotoxicity data are lacking by integrating in vitro genotoxicity assays with PBK modelling-facilitated reverse dosimetry.

Determination of genotoxic potencies of pyrrolizidine alkaloids in HepaRG cells using the γH2AX assay

Pyrrolizidine alkaloids (PAs) are secondary metabolites from plants that have been found in substantial amounts in herbal supplements, infusions and teas. Several PAs cause cancer in animal bioassays, mediated via a genotoxic mode of action, but for the majority of the PAs, carcinogenicity data are lacking. It is assumed in the risk assessment that all PAs have the same potency as riddelliine, which is considered to be one of the most potent carcinogenic PAs in rats. This may overestimate the risks, since many PAs are expected to have lower potencies. In this study we determined the concentration-dependent genotoxicity of 37 PAs representing different chemical classes using the γH2AX in cell western assay in HepaRG human liver cells. Based on these in vitro data, PAs were grouped into different potency classes. The group with the highest potency consists particularly of open diester PAs and cyclic diester PAs (including riddelliine). The group of the least potent or non-active PAs includes the monoester PAs, non-esterified necine bases, PA N-oxides, and the unsaturated PA trachelanthamine. This study reveals differences in in vitro genotoxic potencies of PAs, supporting that the assumption that all PAs have a similar potency as riddelliine is rather conservative.

Intrinsic relative potency of a series of pyrrolizidine alkaloids characterized by rate and extent of metabolism

1,2-Unsaturated pyrrolizidine alkaloids (PAs) are sometimes present in foods or herbal supplements/medicines as impurities and pose potential concerns for liver genotoxicity/carcinogenicity. PAs display a strong structure toxicity relationship, however, current regulatory approaches to risk assessment take the precautionary approach of assuming all PAs display the same potency as the most toxic congeners lasiocarpine (LAS) and riddelliine (RID). Here we explore the relative potencies of a series of structurally diverse PAs by measuring DNA adduct formation in vitro in a rat sandwich culture hepatocyte (SCH) cell system. The adducts generated are consistent with those identified in vivo as biomarkers of PA exposure and potential liver-tumor formation. DNA reactive PAs require metabolic activation to form intermediates that bind DNA, therefore, adduct formation is a direct reflection of reactive metabolite formation. Since the area under the concentration versus time curve (AUC) for the depletion of parent PA from the extracellular media is a measure of PA exposure, the ratio of adducts/AUC provides a measure of hepatocyte exposure to DNA-binding metabolites corresponding to an intrinsic potency for DNA adduct formation. Intrinsic potencies relative to potencies for LAS compare well with existing relative potency data further affirming that PA toxicity varies considerably with chemical structure.

Relative potency of fifteen pyrrolizidine alkaloids to induce DNA damage as measured by micronucleus induction in HepaRG human liver cells

Plant-based 1,2-unsaturated Pyrrolizidine Alkaloids (PAs) can be found as contaminants in foods like teas, herbs and honey. PAs are responsible for liver genotoxicity/carcinogenicity following metabolic activation, making them a relevant concern for safety assessment. Current regulatory risk assessments take a precautionary approach and assume all PAs are as potent as the known most potent representatives: lasiocarpine and riddelliine. Our study investigated whether genotoxicity potency differed as a consequence of structural differences, assessing micronuclei in vitro in HepaRG cells which express metabolising enzymes at levels similar to primary human hepatocytes. Benchmark Dose (BMD) analysis was used to calculate the critical effect dose for 15 PAs representing 6 structural classes. When BMD confidence intervals were used to rank PAs, lasiocarpine was the most potent PA and plotted distinctly from all other PAs examined. PA-N-oxides were least potent, notably less potent than their corresponding parent PA's. The observed genotoxic potency compared favorably with existing in vitro data when metabolic competency was considered. Although further consideration of biokinetics will be needed to develop a robust understanding of relative potencies for a realistic risk assessment of PA mixtures, these data facilitate understanding of their genotoxic potencies and affirm that not all PAs are created equal.

Editor's Highlight: Identification of Any Structure-Specific Hepatotoxic Potential of Different Pyrrolizidine Alkaloids Using Random Forests and Artificial Neural Networks

Pyrrolizidine alkaloids (PAs) are characteristic metabolites of some plant families and form a powerful defense mechanism against herbivores. More than 600 different PAs are known. PAs are ester alkaloids composed of a necine base and a necic acid, which can be used to divide PAs in different structural subcategories. The main target organs for PA metabolism and toxicity are liver and lungs. Additionally, PAs are potentially genotoxic, carcinogenic and exhibit developmental toxicity. Only for very few PAs, in vitro and in vivo investigations have characterized their toxic potential. However, these investigations suggest that structural differences have an influence on the toxicity of single PAs. To investigate this structural relationship for a large number of PAs, a quantitative structural-activity relationship (QSAR) analysis for hepatotoxicity of over 600 different PAs was performed, using Random Forest- and artificial Neural Networks-algorithms. These models were trained with a recently established dataset specific for acute hepatotoxicity in humans. Using this dataset, a set of molecular predictors was identified to predict the hepatotoxic potential of each compound in validated QSAR models. Based on these models, the hepatotoxic potential of the 602 PAs was predicted and the following hepatotoxic rank order in 3 main categories defined (1) for necine base: otonecine > retronecine > platynecine; (2) for necine base modification: dehydropyrrolizidine ≫ tertiary PA = N-oxide; and (3) for necic acid: macrocyclic diester ≥ open-ring diester > monoester. A further analysis with combined structural features revealed that necic acid has a higher influence on the acute hepatotoxicity than the necine base.

Structure-dependent induction of apoptosis by hepatotoxic pyrrolizidine alkaloids in the human hepatoma cell line HepaRG: Single versus repeated exposure

Pyrrolizidine alkaloids (PA) are secondary plant compounds. PA intoxication in humans causes severe acute and chronic hepatotoxicity. However, the molecular mechanisms of PA hepatotoxicity in humans are not well understood yet. Therefore, we investigated cell death parameters in human HepaRG cells following either single (24 h) or repeated dose treatment (14 d) with structurally different PA of the retronecine (echimidine, senecionine), heliotridine (heliotrine), and otonecine type (senkirkine). After 24 h of exposure only retronecine-type PA were cytotoxic in HepaRG cells and induced apoptosis indicated by a loss of membrane asymmetry, disruption of the mitochondrial membrane potential, and increased pro-caspase and PARP cleavage. In contrast, after 14 d all four PA exerted the aforementioned effects. Furthermore, the apoptotic events caspase 3, 8 and 9 activation as well as nuclear condensation and DNA fragmentation were only detected for the retronecine-type PA after single exposure (6 h). Overall, our studies revealed a time- and structure-dependent apoptosis after PA exposure, suggesting that retronecine-type PA seem to be more potent apoptosis inducers than heliotridine- or otonecine-type PA. Furthermore, our results suggest that PA-induced apoptosis in HepaRG cells occur most probably by involving both, the extrinsic death receptor pathway as well as the intrinsic mitochondrial pathway.

1,3-Oxazin-6-one Derivatives and Bohemamine-Type Pyrrolizidine Alkaloids from a Marine-Derived Streptomyces spinoverrucosus

Two new 1,3-oxazin-6-one derivatives (1 and 2) and six new bohemamine-type pyrrolizidine alkaloids (3-8) were isolated from the marine-derived Streptomyces spinoverrucosus strain SNB-048. Their structures including the absolute configurations were fully elucidated on the basis of spectroscopic analysis, ECD spectra, quantum chemical calculations, and chemical methods. Compounds 1 and 2 possess a γ-lactam moiety and a 1,3-oxazin-6-one system.

Disturbance of gene expression in primary human hepatocytes by hepatotoxic pyrrolizidine alkaloids: A whole genome transcriptome analysis

1,2-unsaturated pyrrolizidine alkaloids (PA) are plant metabolites predominantly occurring in the plant families Asteraceae and Boraginaceae. Acute and chronic PA poisoning causes severe hepatotoxicity. So far, the molecular mechanisms of PA toxicity are not well understood. To analyze its mode of action, primary human hepatocytes were exposed to a non-cytotoxic dose of 100 μM of four structurally different PA: echimidine, heliotrine, senecionine, senkirkine. Changes in mRNA expression were analyzed by a whole genome microarray. Employing cut-off values with a |fold change| of 2 and a q-value of 0.01, data analysis revealed numerous changes in gene expression. In total, 4556, 1806, 3406 and 8623 genes were regulated by echimidine, heliotrine, senecione and senkirkine, respectively. 1304 genes were identified as commonly regulated. PA affected pathways related to cell cycle regulation, cell death and cancer development. The transcription factors TP53, MYC, NFκB and NUPR1 were predicted to be activated upon PA treatment. Furthermore, gene expression data showed a considerable interference with lipid metabolism and bile acid flow. The associated transcription factors FXR, LXR, SREBF1/2, and PPARα/γ/δ were predicted to be inhibited. In conclusion, though structurally different, all four PA significantly regulated a great number of genes in common. This proposes similar molecular mechanisms, although the extent seems to differ between the analyzed PA as reflected by the potential hepatotoxicity and individual PA structure.