Next Article in Journal
Effect of Pyridine on the Mesophase of Teraryl Liquid Crystals: A New Series of Nematic Liquid Crystals Named 2-(4-Alkoxybiphen-4′-yl)-5-methylpyridines
Next Article in Special Issue
Confocal Spectroscopy to Study Dimerization, Oligomerization and Aggregation of Proteins: A Practical Guide
Previous Article in Journal
Oncostatic-Cytoprotective Effect of Melatonin and Other Bioactive Molecules: A Common Target in Mitochondrial Respiration
Previous Article in Special Issue
Insights into Protein–Ligand Interactions: Mechanisms, Models, and Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Toxic Pyrrolizidine Alkaloids and Their Common Hepatotoxicity Mechanism

1
School of Life Sciences and Technology, Tongji University, Shanghai 200092, China
2
School of Biomedical Informatics, University of Texas Health Science Center, Houston, TX 77030, USA
3
The Ministry of Education (MOE) Key Laboratory for Standardization of Chinese Medicines, Institute of Chinese Materia Medica, Shanghai University of Traditional Chinese Medicine, Shanghai 201203, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(3), 318; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030318
Submission received: 30 December 2015 / Revised: 19 February 2016 / Accepted: 24 February 2016 / Published: 7 March 2016
(This article belongs to the Special Issue Proteins and Protein-Ligand Interactions)

Abstract

:
Pyrrolizidine Alkaloids (PAs) are currently one of the most important botanical hepatotoxic ingredients. Glutathion (GSH) metabolism is the most reported pathway involved in hepatotoxicity mechanism of PAs. We speculate that, for different PAs, there should be a common mechanism underlying their hepatotoxicity in GSH metabolism. Computational methods were adopted to test our hypothesis in consideration of the limitations of current experimental approaches. Firstly, the potential targets of 22 PAs (from three major PA types) in GSH metabolism were identified by reverse docking; Secondly, glutathione S-transferase A1 (GSTA1) and glutathione peroxidase 1 (GPX1) targets pattern was found to be a special characteristic of toxic PAs with stepwise multiple linear regressions; Furthermore, the molecular mechanism underlying the interactions within toxic PAs and these two targets was demonstrated with the ligand-protein interaction analysis; Finally, GSTA1 and GPX1 were proved to be significant nodes in GSH metabolism. Overall, toxic PAs could be identified by GSTA1 and GPX1 targets pattern, which suggests their common hepatotoxicity mechanism: the interfering of detoxication in GSH metabolism. In addition, all the strategies developed here could be extended to studies on toxicity mechanism of other toxins.

Graphical Abstract

1. Introduction

Pyrrolizidine Alkaloids (PAs) are well known as a type of the most hepatotoxic toxins in nature. Currently, more than 660 kinds of isolated PAs and their nitrogen oxide derivatives have been identified in over 6000 plants [1]. PAs mainly appear in four plant families (Boraginaceae, Asteraceae, Leguminosae and Orchidaceae [2,3,4]), and more than half of PAs produce toxicity [3]. Because of their abundance in nature and potent toxicity, including serious hepatotoxicity, carcinogenicity, pneumotoxicity, neurotoxicity and embrotoxicityc [5,6], PAs have received extensive attention.
Decades of efforts have been made to study the hepatotoxicity mechanism of PAs. In vitro experiment was regarded as an important approach in studying metabolic toxicity, and cell models including Primary cultured hepatocytes [7,8], L-02 liver cells [9,10], Human hepatoma cell HepG2 [11,12] and Human hepatoma cells HuH-7 [13,14] have been established. In the light of studies on the existing models, it is generally considered that PAs exhibit hepatotoxicity via their activated metabolite, pyrrolic ester, which was catalyzed by the hepatic cytochrome P450 (CYP450) enzyme system [15,16].
The PAs induced hepatotoxicity has long been considered to associate with many factors, such as reactive oxygen species (ROS), CYPs and GSH-metabolizing enzymes [17]. Most recently, accumulating evidences have demonstrated that GSH metabolism plays an important role in regulating the hepatotoxicity induced by PAs. In vivo experiments in our team reported that intracellular GSH effect on regulating the hepatotoxicity induced by Adon, Sene and Mone [18,19,20]. In addition, cellular GSH was showed with the ability to prevent hepatotoxicity induced by clivorine [18]. Some studies suggested that oxidative stress plays a crucial role in PAs induced hepatotoxicity [21]. Moreover, it has been showed that the activities of GSH-related antioxidant enzymes were changed by PAs. For instance, isoline induced the change of glutathione S-transferase (GST), glutathione peroxidase (GPX) and glycine N-methyltransferase (GNMT) [22]. As components of the intracellular GSH system, the activities of GPX, GST and glutathione reductase (GR) decreased significantly on clivorine-treated primary cultured rat hepatocytes [23]. Nine differentially expressed proteins were identified in injured liver of isoline-treated mice. Furthermore, two potential protein targets (GST-pi and ACADM) of isoline were identified as possible signaling molecules related to isoline-induced liver injury [24]. Besides, decreases of GR, GPX and GST activities were induced by monocrotaline [25,26].
However, due to the limitation of current experimental approaches, there are still many PAs of which hepatotoxicity mechanism remains unclear under GSH antioxidant system. If the relationship within the PAs and enzymes involved in the GSH metabolism can be determined systematically, it will lay a solid foundation for future studies concerning PAs hepatotoxicity, and promote the development of related drugs.
In this study, we aim to investigate the common hepatotoxicity mechanism of a large scale of PAs from the three major PAs types in GSH metabolism. Inverse docking and multiple stepwise linear regression were used in mass screening of the potential protein targets of 22 PAs from three major PA types in GSH metabolism. We expect to find the boundary which can distinguish toxic PAs and non-toxic compounds, and consequently uncover the relationship between PAs and liver enzymes in the GSH metabolism. Compared with nontoxic compounds, toxic PAs can be identified by a potential targets pattern, named as GSTA1 and GPX1. It is a possible molecular pattern which can accurately identify the toxic PAs. The molecular mechanisms of the interactions within these two targets and PAs were detected respectively, and the importance of them in GSH metabolism was further confirmed in this study. GSTA1 and GPX1 targets pattern indicates there is a common hepatotoxicity mechanism for toxic PAs, that these toxic PAs could break the balance of free radicals in human body and decrease the excretion of the exogenous toxins by simultaneously interfering two important detoxication pathways related to GSTA1 and GPX1, thus indirectly causing the toxicity on hepatocytes.

2. Results

2.1. Identifying the Potential Targets of PAs

In order to identify the potential targets of PAs, two datasets of the small molecules were used to detect the putative targets by reverse docking (Table S1). The positive dataset consists of 22 toxic PAs from three PAs types. For each PA, information on the prototype and its hydrolysis products are included. The negative dataset is composed of five non-toxic open loop PAs and 8 non-PA compounds. These eight non-PA compounds are flavonoids with the same weight to PAs and were extracted as by-products in the extraction of PAs in the laboratory. There are 51 proteases involved in the glutathione metabolism in KEGG [27], 34 of which have 3D structures (Table S2) [28]. Therefore, 35 small molecules were reversely docked to 34 proteases, and 14 proteases were hit. The results were illustrated in Figure 1.
The targets number distribution of the toxic PAs and non-toxic compounds are different (Figure 1b). The potential targets of non-toxic compounds are decentralized, and the median of their targets number is far lower than these of toxic PAs. In contrast, targets of toxic PAs are concentrated, and focused on GSTA1, GPX1, glutathione S-transferase pi 1(GSTP1), leucine amiopeptidase (LAP3) and isocitrate dehydrogenase [NADP] (IDH2).

2.2. Targets Enrichment of Toxic PAs

To find some commonalities of targets to distinguish toxic PAs from non-toxic compounds, specific targets of toxic PAs were enriched with four stepwise multiple linear regression models. There are four targets enriched in the best model (Table 1), according to the standardized coefficients, GSTA1 and GPX1 weight 50.2% and 69.7% of the toxic prediction precision respectively, while GR and LAP3 weight 68.3% and 33.4% of the non-toxic prediction precision respectively. These targets were combined into different patterns to distinguish the toxic and non-toxic compounds. GSTA1 and GPX1 pattern reached the highest precision with 94.29% (Table 2), indicating toxic PAs’ preference for binding to GSTA1 and GPX1.

2.3. Molecular Interactions within PAs and GSTA1, GPX1

To better understand the interaction of PAs and GSTA1 and GPX1 pattern, active necine bases dehydroretronecine (DHR) (R-6,7-dihydro-7-hydroxy-1-hydroxymethyl-5H-pyrrolizine) were used to dock GSTA1 and GPX1 respectively (Figure 2a1,a2). The molecular interactions were shown in Figure 2b1,b2. GSTA1-DHR interaction shows that the O on 1-hydroxymethyl of DHR accepts a hydrogen from the Arg45, while the H on the 1-hydroxymethyl donates a hydrogen to Gln53 (Figure 2b1), and OH7 bridges the hydrogen bond with Phe220 side chain via a water molecule. GPX1-DHR interaction shows that the OH on 1-hydroxymethyl accepts two hydrogens from the Arg180 and Thr143, and it also bridges hydrogen bonds with them via water molecules. The OH7 not only accepts hydrogen from Arg179, it also bridges a hydrogen with Arg179 by a water molecule. There also exists a π-π stacking interaction between pyrrole and Trp160 (Figure 2b2). These indicate the molecular mechanisms underlying how toxic PAs interact with GSTA1 and GPX1, and consequently interfere with the activities and functions of them.

2.4. Importance of the GSTA1 and GPX1 in GSH Metabolism

GSTA1 and GPX1 were analyzed in the protein-protein interaction (PPI) network to confirm their importance in GSH metabolism. A PPI network with 91 proteins and 194 interactions was produced by linking each pair of 44 proteases in GSH metabolism with the shortest path, and the importance of these proteins was measured by betweenness, degree and pagerank (Figure 3a). Under these three measurements, it’s interesting to find that GSTA1 and GPX1 were also the top ranking ones among the 44 proteins in GSH pathway (Figure 3b–d). This suggests that GSTA1 and GPX1 are the crucial nodes in the network, if they were interfered, the interaction networks would be heavily perturbed, which would lead to dysfunction of GSH metabolism.

3. Discussion

3.1. Toxicity Mechanism Analysis of Toxic PAs Based on GSTA1 and GPX1 Pattern

Our study identified the potential targets of PAs in GSH metabolism, and demonstrated the GSTA1 and GPX1 targets pattern can distinguish toxic PAs from non-toxic compounds with the accuracy over 94%. This pattern suggests toxic PAs probably bind to both these two proteins and influence their function. Furthermore, we have detected potential molecular interactions between secondary metabolites of toxic PAs and GSTA1, GPX1 respectively, which helps understand how PAs affect these two important enzymes.
In GSH metabolism (Figure 4), as an important intracellular anti-oxidant, cellular GSH directly takes part in the detoxication of free radicals, carcinogens and peroxides induced oxidative stress [29,30,31,32]. The detoxication effect of GSH could be partly attributed to its role as a cofactor with anti-oxidant enzymes (such as GPX and GST) in accelerating the excretion of toxins [33].
Intracellular GSH could serve as the electron donor in the reduction of peroxides catalyzed by GPX [34]. As a glutathione peroxidase, GPX1 protects cells from oxidative stress by catalyzing the reduction of hydrogen peroxide to water [35,36]. Indirectly, GPX1 is also related to the NADP+/NADPH redox couple, helping adjust the GSH/GSSG ratios to normal level. Generally, GPX1 is considered as a crucial antioxidant enzyme against oxidative stress, and has been found to be effective at removing intracellular peroxides under many physiological conditions [37,38].
Besides, under the catalysis of GST, GSH can conjugate the metabolites of exogenous toxins [39]. GSTs are a class of phase II drug metabolism enzymes, which act on the detoxication of electrophilic compounds capable of generating reactive oxygen species [40,41,42]. They catalyze the conjugation of GSH to these exogenous substrates via a sulfhydryl group-electrophilic center which makes the compounds more water-soluble [43,44]. In this way, they accelerate the excretion of the toxics [45] and decrease the level of oxidative stress, so as to achieve the goal of antioxidant defense and detoxication [46]. GSTA1 has been proved associated with the oxidative stress mediated by hydroxyl radicals (zOH) and the GSTA1 inducers with pro-oxidative potential [47]. It has also been found as a marker of hepatocyte injury in transplantation, toxicity and viral infections [48,49,50].
Taken together, toxic PAs could bind with GSTA1 and GPX1 in liver and interfere two important detoxication pathways in GSH metabolism by decreasing the activities of these two major cellular antioxidant enzymes. This would disrupt the balance of free radicals in the body and hinders the excretion of toxic compounds out of the cell, therefore indirectly resulting in the toxicity on hepatocytes.

3.2. Limitations of Reverse Docking

Reverse docking is an important method to calculate and predict the targets of small molecules in a large scale as well as in a non-laboratory environment. However, it yet has two main limitations. On the one hand, it may generate false positives from computing models and docking procedures. On the other hand, there are still no available 3D structures for some receptors, which is the foundation of docking study. Here, only receptors with 3D structures were selected for molecular docking, which may inevitably generate biased results. Besides, there is a need to further optimize computation models and statistical models to decrease false positives. Despite all this, our study could still lay a foundation to future PAs’ toxicity mechanism research, especially to the researches on the detoxication mechanisms regarding the GSTA1 and GPX1.

4. Materials and Methods

4.1. Dataset

The positive and negative datasets of PAs (Table S1) are retrievable from The Ministry of Education (MOE) Key Laboratory of Standardization of Chinese Medicines, Institute of Chinese Materia Medica, Shanghai University of Traditional Chinese Medicine, Shanghai, China.
The 3D structures of 34 proteins involved in the Glutathione metabolism in KEGG are acquired from protein data bank (PDB) (Table S2).

4.2. Identification of Potential Protein Targets of PAs

INVDOCK [51] was used as a reverse docking tool to find the potential targets of PAs. The targets of a small molecule were identified by computer-automatically simulation, by docking the molecule into all cavities of each protein in PDB database (http://www.rcsb.org/) [52,53,54]. Potential targets were selected based on the Affinity (ΔELP) between proteins and ligands according to a specific scoring function. A protein could be considered as a potential target of certain small molecule when ΔELP is simultaneously less than two thresholds [53,55]. The first threshold is the experience score (ΔE Threshold) obtained by statistical analysis of a large number of “protein-ligand” binding energy in PDB database. The other one is the energy score (ΔE Competitor) calculated by scoring function on natural ligand-protein crystal complex.
In this study, only when a protein was identified as the target of all the three states of a same PA, could it be considered as the potential target of this PA. As a result, 14 out of 34 proteins were hit by PAs.

4.3. Targets Enrichment of Toxic PAs

Stepwise multiple linear regression analysis was applied to seek the specific targets to distinguish toxic PAs (positive dataset) and non-toxic compounds (negative dataset). Four linear regression models were established based on the targets information, and the differences among them lie in the targets set of positive dataset, shown as follows:
Model 1: Targets set only contains the targets of the PAs prototype.
Model 2: Targets set only contains the targets of DHP.
Model 3: Targets set contains the targets of all small molecules, whether they are PAs prototype, DHP or DHR.
Model 4: Targets set contains the targets which were hit by all the three states of certain PA at the same time.
Key targets filtered out from each of the four models were further combined to predict whether a molecule is a toxic PA. The information of other three models and prediction precisions of different combinations of features are shown in Supplementary Table S3.

4.4. Molecular Interactions within Toxic PAs and GSTA1, GPX1

The binding sites and possible molecular interactions within toxic PAs and GSTA1, GPX1 were identified by Simulations module in the commercial software Molecular Operating Environment (MOE).

4.5. Importance of the Targets Pattern in GSH Metabolism

In order to investigate the connection of proteases related to the GSH metabolism on the human protein-protein interaction network links (PPI), in particular the importance of GSTA1 and GPx1 pattern, we constructed a subnetwork for GSH metabolism.
Firstly, we constructed a human PPI which composed of 55,363 edges and 10,523 nodes based on 6 PPI databases, including HPRD (Human Protein Reference Database) [56], MINT (the Molecular INTeraction database) [57], BioGRID (Biological General Repository for Interaction Datasets) [58], IntAct (IntAct Molecular Interaction Database) [59], DIP (Database of Interacting Proteins) [60] and MIPS (Mammalian Protein-Protein Interaction Database) [61].
Secondly, we mapped 44 proteases in GSH metabolism on the human PPI and made them interconnected by remaining some background nodes. Under this condition, a GSH metabolism network with 91 nodes and 194 edges was constructed (Figure 3a). Finally, we investigated the importance of the nodes with the measurements of degree, betweenness and pagerank.

5. Conclusions

Overall, our study suggests that toxic PAs and their metabolites are most likely to be identified by GSTA1 and GPX1 targets pattern. After binding with GSTA1 and GPX1, toxic PAs could further exert toxicity by breaking the oxidant-antioxidant balance and intercepting detoxication of the exogenous toxins. This pattern could provide some support in the future studies on PAs’ hepatotoxicity mechanism. Moreover, the computational workflow employed in this study may throw light upon the toxic mechanisms studies on other toxins.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/17/3/318/s1.

Acknowledgments

This work was supported in part by grants from the National Natural Science Foundation of China (31200986, 41530105, and 31171272) and the Fundamental Research Funds for the Central Universities (10247201546, 2000219083). The funders played no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author Contributions

Zhiwei Cao conceived and designed the project; Xinmiao Yan, Hong Kang, Jun Feng and Yiyan Yang collected data and carried out the analytical procedures; Xinmiao Yan, Hong Kang, Jun Feng, Yiyan Yang, Kailin Tang, Ruixin Zhu, Li Yang, Zhengtao Wang and Zhiwei Cao interpreted the results; Xinmiao Yan drafted the manuscript; Hong Kang, Jun Feng, Yiyan Yang, Kailin Tang, Ruixin Zhu, Li Yang, Zhengtao Wang and Zhiwei Cao revised the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chou, M.W.; Fu, P.P. Formation of DHP-derived DNA adducts in vivo from dietary supplements and chinese herbal plant extracts containing carcinogenic pyrrolizidine alkaloids. Toxicol. Ind. Health 2006, 22, 321–327. [Google Scholar] [PubMed]
  2. Roeder, E. Medicinal plants in Europe containing pyrrolizidine alkaloids. Die Pharm. 1995, 50, 83–98. [Google Scholar]
  3. Roeder, E. Medicinal plants in China containing pyrrolizidine alkaloids. Die Pharm. 2000, 55, 711–726. [Google Scholar]
  4. Stegelmeier, B.L.; Edgar, J.A.; Colegate, S.M.; Gardner, D.R.; Schoch, T.K.; Coulombe, R.A.; Molyneux, R.J. Pyrrolizidine alkaloid plants, metabolism and toxicity. J. Nat. Toxins 1999, 8, 95–116. [Google Scholar] [PubMed]
  5. Fu, P.P.; Chou, M.W.; Xia, Q.; Yang, Y.C.; Yan, J.; Doerge, D.R.; Chan, P.C. Genotoxic pyrrolizidine alkaloids and pyrrolizidine alkaloidn-oxides—Mechanisms leading to DNA adduct formation and tumorigenicity. J. Environ. Sci. Health Part C 2001, 19, 353–385. [Google Scholar] [CrossRef]
  6. Seaman, J.T. Pyrrolizidine alkaloid poisoning of sheep in New South Wales. Aust. Vet. J. 1987, 64, 164–167. [Google Scholar] [CrossRef] [PubMed]
  7. Green, C.E.; Segall, J.H.; Byard, J.L. Metabolism, cytotoxicity, and genotoxicity of the pyrrolizidine alkaloid senecionine in primary cultures of rat hepatocytes. Toxicol. Appl. Pharmacol. 1981, 60, 176–185. [Google Scholar] [CrossRef]
  8. Williams, G.M.; Mori, H.; Hirono, I.; Nagao, M. Genotoxicity of pyrrolizidine alkaloids in the hepatocyte primary culture/DNA-repair test. Mutat. Res. 1980, 79, 1–5. [Google Scholar] [CrossRef]
  9. Ji, L.L.; Zhang, M.; Sheng, Y.C.; Wang, Z.T. Pyrrolizidine alkaloid clivorine induces apoptosis in human normal liver L-02 cells and reduces the expression of p53 protein. Toxicol. In Vitro 2005, 19, 41–46. [Google Scholar] [CrossRef] [PubMed]
  10. Ji, L.L.; Zhao, X.G.; Chen, L.; Zhang, M.; Wang, Z.T. Pyrrolizidine alkaloid clivorine inhibits human normal liver L-02 cells growth and activates p38 mitogen-activated protein kinase in L-02 cells. Toxicon 2002, 40, 1685–1690. [Google Scholar] [CrossRef]
  11. Uhl, M.; Helma, C.; Knasmuller, S. Evaluation of the single cell gel electrophoresis assay with human hepatoma (Hep G2) cells. Mutat. Res. 2000, 468, 213–225. [Google Scholar] [CrossRef]
  12. Rizo, W.F.; Ferreira, L.E.; Colnaghi, V.; Martins, J.S.; Franchi, L.P.; Takahashi, C.S.; Beleboni, R.O.; Marins, M.; Pereira, P.S.; Fachin, A.L. Cytotoxicity and genotoxicity of coronaridine from Tabernaemontana catharinensis A.DC in a human laryngeal epithelial carcinoma cell line (Hep-2). Genet. Mol. Biol. 2013, 36, 105–110. [Google Scholar] [CrossRef] [PubMed]
  13. Zuckerman, M.; Steenkamp, V.; Stewart, M.J. Hepatic veno-occlusive disease as a result of a traditional remedy: Confirmation of toxic pyrrolizidine alkaloids as the cause, using an in vitro technique. J. Clin. Pathol. 2002, 55, 676–679. [Google Scholar] [CrossRef] [PubMed]
  14. Steenkamp, V.; Stewart, M.J.; van der Merwe, S.; Zuckerman, M.; Crowther, N.J. The effect of Senecio latifolius a plant used as a South African traditional medicine, on a human hepatoma cell line. J. Ethnopharmacol. 2001, 78, 51–58. [Google Scholar] [CrossRef]
  15. Gordon, G.J.; Coleman, W.B.; Grisham, J.W. Induction of cytochrome P450 enzymes in the livers of rats treated with the pyrrolizidine alkaloid retrorsine. Exp. Mol. Pathol. 2000, 69, 17–26. [Google Scholar] [CrossRef] [PubMed]
  16. Fu, P.P.; Xia, Q.; Lin, G.; Chou, M.W. Pyrrolizidine alkaloids—Genotoxicity, metabolism enzymes, metabolic activation, and mechanisms. Drug Metab. Rev. 2004, 36, 1–55. [Google Scholar] [CrossRef] [PubMed]
  17. Neuman, M.G.; Cohen, L.; Opris, M.; Nanau, R.M.; Hyunjin, J. Hepatotoxicity of pyrrolizidine alkaloids. J. Pharm. Pharm. Sci. 2015, 18, 825–843. [Google Scholar] [PubMed]
  18. Ji, L.; Chen, Y.; Wang, Z. Protection of S-adenosyl methionine against the toxicity of clivorine on hepatocytes. Environ. Toxicol. Pharmacol. 2008, 26, 331–335. [Google Scholar] [CrossRef] [PubMed]
  19. Chen, Y.; Ji, L.; Wang, H.; Wang, Z. Intracellular glutathione plays important roles in pyrrolizidine alkaloids-induced growth inhibition on hepatocytes. Environ. Toxicol. Pharmacol. 2009, 28, 357–362. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, Y.; Ji, L.; Xiong, A.; Yang, L.; Wang, Z. Involvement of intracellular glutathione in regulating isoline-induced cytotoxicity in human normal liver L-02 cells. Toxicol. Ind. Health 2013, 29, 567–575. [Google Scholar] [CrossRef] [PubMed]
  21. Liang, Q.; Sheng, Y.; Jiang, P.; Ji, L.; Xia, Y.; Min, Y.; Wang, Z. The gender-dependent difference of liver GSH antioxidant system in mice and its influence on isoline-induced liver injury. Toxicology 2011, 280, 61–69. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, T.Y.; Chen, Y.; Wang, Z.Y.; Ji, L.L.; Wang, Z.T. Pyrrolizidine alkaloid isoline-induced oxidative injury in various mouse tissues. Exp. Toxicol. Pathol. 2010, 62, 251–257. [Google Scholar] [CrossRef] [PubMed]
  23. Ji, L.; Liu, T.; Wang, Z. Pyrrolizidine alkaloid clivorine induced oxidative injury on primary cultured rat hepatocytes. Hum. Exp. Toxicol. 2010, 29, 303–309. [Google Scholar] [PubMed]
  24. Wang, Z.Y.; Kang, H.; Ji, L.L.; Yang, Y.Q.; Liu, T.Y.; Cao, Z.W.; Morahan, G.; Wang, Z.T. Proteomic characterization of the possible molecular targets of pyrrolizidine alkaloid isoline-induced hepatotoxicity. Environ. Toxicol. Pharmacol. 2012, 34, 608–617. [Google Scholar] [CrossRef] [PubMed]
  25. Amin, K.A.; Hashem, K.S.; Al-muzafar, H.M.; Taha, E.M. Oxidative hepatotoxicity effects of monocrotaline and its amelioration by lipoic acid, S-adenosyl methionine and vitamin E. J. Complement. Integr. Med. 2014, 11, 35–41. [Google Scholar] [CrossRef] [PubMed]
  26. Amin, K.A.; Hassan, M.S.; el Awad, S.T.; Hashem, K.S. The protective effects of cerium oxide nanoparticles against hepatic oxidative damage induced by monocrotaline. Int. J. Nanomed. 2011, 6, 143–149. [Google Scholar] [CrossRef] [PubMed]
  27. Kanehisa, M. The KEGG database. Silico Simul. Biol. Process. 2002, 247, 91–103. [Google Scholar]
  28. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The protein data bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed]
  29. Park, E.J.; Park, K. Induction of reactive oxygen species and apoptosis in BEAS-2B cells by mercuric chloride. Toxicol. In Vitro 2007, 21, 789–794. [Google Scholar] [CrossRef] [PubMed]
  30. Sakurai, T.; Ochiai, M.; Kojima, C.; Ohta, T.; Sakurai, M.H.; Takada, N.O.; Qu, W.; Waalkes, M.P.; Himeno, S.; Fujiwara, K. Preventive mechanism of cellular glutathione in monomethylarsonic acid-induced cytolethality. Toxicol. Appl. Pharmacol. 2005, 206, 54–65. [Google Scholar] [CrossRef] [PubMed]
  31. Yadav, U.C.; Ramana, K.V.; Awasthi, Y.C.; Srivastava, S.K. Glutathione level regulates HNE-induced genotoxicity in human erythroleukemia cells. Toxicol. Appl. Pharmacol. 2008, 227, 257–264. [Google Scholar] [CrossRef] [PubMed]
  32. Dukhande, V.V.; Malthankar-Phatak, G.H.; Hugus, J.J.; Daniels, C.K.; Lai, J.C. Manganese-induced neurotoxicity is differentially enhanced by glutathione depletion in astrocytoma and neuroblastoma cells. Neurochem. Res. 2006, 31, 1349–1357. [Google Scholar] [CrossRef] [PubMed]
  33. Hayes, J.D.; Pulford, D.J. The glutathione S-transferase supergene family: Regulation of GST and the contribution of the isoenzymes to cancer chemoprotection and drug resistance. Crit. Rev. Biochem. Mol. Biol. 1995, 30, 445–600. [Google Scholar] [CrossRef] [PubMed]
  34. Lash, L.H.; Putt, D.A.; Zalups, R.K. Influence of exogenous thiols on inorganic mercury-induced injury in renal proximal and distal tubular cells from normal and uninephrectomized rats. J. Pharmacol. Exp. Ther. 1999, 291, 492–502. [Google Scholar] [PubMed]
  35. Brigelius-Flohe, R. Tissue-specific functions of individual glutathione peroxidases. Free Radic. Biol. Med. 1999, 27, 951–965. [Google Scholar] [CrossRef]
  36. Higashi, Y.; Pandey, A.; Goodwin, B.; Delafontaine, P. Insulin-like growth factor-1 regulates glutathione peroxidase expression and activity in vascular endothelial cells: Implications for atheroprotective actions of insulin-like growth factor-1. Biochim. Biophys. Acta 2013, 1832, 391–399. [Google Scholar] [CrossRef] [PubMed]
  37. Antunes, F.; Han, D.; Cadenas, E. Relative contributions of heart mitochondria glutathione peroxidase and catalase to H2O2 detoxification in in vivo conditions. Free Radic. Biol. Med. 2002, 33, 1260–1267. [Google Scholar] [CrossRef]
  38. Cohen, G.; Hochstein, P. Glutathione Peroxidase: The Primary agent for the elimination of hydrogen peroxide in erythrocytes. Biochemistry 1963, 2, 1420–1428. [Google Scholar] [CrossRef] [PubMed]
  39. Miura, T.; Muraoka, S.; Fujimoto, Y. Lipid peroxidation induced by phenylbutazone radicals. Life Sci. 2002, 70, 2611–2621. [Google Scholar] [CrossRef]
  40. Sheehan, D.; Meade, G.; Foley, V.M.; Dowd, C.A. Structure, function and evolution of glutathione transferases: Implications for classification of non-mammalian members of an ancient enzyme superfamily. Biochem. J. 2001, 360, 1–16. [Google Scholar] [CrossRef] [PubMed]
  41. Leaver, M.J.; George, S.G. A piscine glutathione S-transferase which efficiently conjugates the end-products of lipid peroxidation. Mar. Environ. Res. 1998, 46, 71–74. [Google Scholar] [CrossRef]
  42. Litwack, G.; Ketterer, B.; Arias, I.M. Ligandin: A hepatic protein which binds steroids, bilirubin, carcinogens and a number of exogenous organic anions. Nature 1971, 234, 466–467. [Google Scholar] [CrossRef] [PubMed]
  43. Douglas, K.T. Mechanism of action of glutathione-dependent enzymes. Adv. Enzymol. Relat. Areas Mol. Biol. 1987, 59, 103–167. [Google Scholar] [PubMed]
  44. Oakley, A. Glutathione transferases: A structural perspective. Drug Metab. Rev. 2011, 43, 138–151. [Google Scholar] [CrossRef] [PubMed]
  45. Mannervik, B.; Alin, P.; Guthenberg, C.; Jensson, H.; Tahir, M.K.; Warholm, M.; Jornvall, H. Identification of three classes of cytosolic glutathione transferase common to several mammalian species: Correlation between structural data and enzymatic properties. Proc. Natl. Acad. Sci. USA 1985, 82, 7202–7206. [Google Scholar] [CrossRef] [PubMed]
  46. Sies, H. Oxidative stress: From basic research to clinical application. Am. J. Med. 1991, 91, 31S–38S. [Google Scholar] [CrossRef]
  47. Pinkus, R.; Weiner, L.M.; Daniel, V. Role of oxidants and antioxidants in the induction of AP-1, NF-κB, and glutathione S-transferase gene expression. J. Biol. Chem. 1996, 271, 13422–13429. [Google Scholar] [PubMed]
  48. Beckett, G.J.; Chapman, B.J.; Dyson, E.H.; Hayes, J.D. Plasma glutathione S-transferase measurements after paracetamol overdose: Evidence for early hepatocellular damage. Gut 1985, 26, 26–31. [Google Scholar] [CrossRef] [PubMed]
  49. Hughes, V.F.; Trull, A.K.; Gimson, A.; Friend, P.J.; Jamieson, N.; Duncan, A.; Wight, D.G.D.; Prevost, A.T.; Alexander, G.J.M. Randomized trial to evaluate the clinical benefits of serumα-glutathione S-transferase concentration monitoring after liver transplantation1. Transplantation 1997, 64, 1446–1452. [Google Scholar] [CrossRef] [PubMed]
  50. Loguercio, C.; Caporaso, N.; Tuccillo, C.; Morisco, F.; del Vecchio Blanco, G.; del Vecchio Blanco, C. Alpha-glutathione transferases in HCV-related chronic hepatitis: A new predictive index of response to interferon therapy? J. Hepatol. 1998, 28, 390–395. [Google Scholar] [CrossRef]
  51. Chen, Y.Z.; Ung, C.Y. Computer automated prediction of potential therapeutic and toxicity protein targets of bioactive compounds from Chinese medicinal plants. Am. J. Chin. Med. 2002, 30, 139–154. [Google Scholar] [CrossRef] [PubMed]
  52. Cao, Z.W.; Chen, Y.Z. Hydrogen-bond disruption probability in proteins by a modified self-consistent harmonic approach. Biopolymers 2001, 58, 319–328. [Google Scholar] [CrossRef]
  53. Chen, Y.Z.; Ung, C.Y. Prediction of potential toxicity and side effect protein targets of a small molecule by a ligand-protein inverse docking approach. J. Mol. Graph. Model. 2001, 20, 199–218. [Google Scholar] [CrossRef]
  54. Rose, P.W.; Bi, C.; Bluhm, W.F.; Christie, C.H.; Dimitropoulos, D.; Dutta, S.; Green, R.K.; Goodsell, D.S.; Prlic, A.; Quesada, M.; et al. The RCSB Protein Data Bank: New resources for research and education. Nucleic Acids Res. 2013, 41, D475–D482. [Google Scholar] [CrossRef] [PubMed]
  55. Chen, Y.Z.; Zhi, D.G. Ligand-protein inverse docking and its potential use in the computer search of protein targets of a small molecule. Proteins-Struct. Funct. Genet. 2001, 43, 217–226. [Google Scholar] [CrossRef]
  56. Prasad, T.S.K.; Kandasamy, K.; Pandey, A. Human protein reference database and human proteinpedia as discovery tools for systems biology. In Reverse Chemical Genetics; Koga, H., Ed.; Humana Press: New York, NY, USA, 2009; Volume 577, pp. 67–79. [Google Scholar]
  57. Licata, L.; Briganti, L.; Peluso, D.; Perfetto, L.; Iannuccelli, M.; Galeota, E.; Sacco, F.; Palma, A.; Nardozza, A.P.; Santonico, E.; et al. MINT, the molecular interaction database: 2012 update. Nucleic Acids Res. 2012, 40, D857–D861. [Google Scholar] [CrossRef] [PubMed]
  58. Chatr-Aryamontri, A.; Breitkreutz, B.J.; Heinicke, S.; Boucher, L.; Winter, A.; Stark, C.; Nixon, J.; Ramage, L.; Kolas, N.; O’Donnell, L.; et al. The BioGRID interaction database: 2013 update. Nucleic Acids Res. 2013, 41, D816–D823. [Google Scholar] [CrossRef] [PubMed]
  59. Orchard, S.; Ammari, M.; Aranda, B.; Breuza, L.; Briganti, L.; Broackes-Carter, F.; Campbell, N.H.; Chavali, G.; Chen, C.; del-Toro, N.; et al. The MIntAct project—IntAct as a common curation platform for 11 molecular interaction databases. Nucleic Acids Res. 2014, 42, D358–D363. [Google Scholar] [CrossRef] [PubMed]
  60. Salwinski, L.; Miller, C.S.; Smith, A.J.; Pettit, F.K.; Bowie, J.U.; Eisenberg, D. The Database of Interacting Proteins: 2004 update. Nucleic Acids Res. 2004, 32, D449–D451. [Google Scholar] [CrossRef] [PubMed]
  61. Pagel, P.; Kovac, S.; Oesterheld, M.; Brauner, B.; Dunger-Kaltenbach, I.; Frishman, G.; Montrone, C.; Mark, P.; Stumpflen, V.; Mewes, H.W.; et al. The MIPS mammalian protein-protein interaction database. Bioinformatics 2005, 21, 832–834. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Potential protein targets of Pyrrolizidine Alkaloids (PAs): (a) Targets sets of toxic PAs and non-toxic compounds; (b) The distribution of the targets number.
Figure 1. Potential protein targets of Pyrrolizidine Alkaloids (PAs): (a) Targets sets of toxic PAs and non-toxic compounds; (b) The distribution of the targets number.
Ijms 17 00318 g001
Figure 2. Molecular interactions between dehydroretronecine (DHR) and glutathione S-transferase A1 (GSTA1), glutathione peroxidase 1 (GPX1). (a1) Docking demonstration of DHR and GSTA1; (a2) Docking demonstration of DHR and GPX1. Red ribbon: helix; blue ribbon: loop; yellow ribbon: beta-strand; pink and green: surface; (b1) Interaction between DHR and GSTA1; (b2) Interaction between DHR and GPX1.
Figure 2. Molecular interactions between dehydroretronecine (DHR) and glutathione S-transferase A1 (GSTA1), glutathione peroxidase 1 (GPX1). (a1) Docking demonstration of DHR and GSTA1; (a2) Docking demonstration of DHR and GPX1. Red ribbon: helix; blue ribbon: loop; yellow ribbon: beta-strand; pink and green: surface; (b1) Interaction between DHR and GSTA1; (b2) Interaction between DHR and GPX1.
Ijms 17 00318 g002
Figure 3. Network analysis of GSTA1 and GPX1 targets pattern. (a) protein-protein interaction (PPI) network constructed by 44 enzymes in GSH metabolism; (bd) measurements of proteins in network, blue represents the proteases in GSH metabolism; red represents GSTA1 and GPX1.
Figure 3. Network analysis of GSTA1 and GPX1 targets pattern. (a) protein-protein interaction (PPI) network constructed by 44 enzymes in GSH metabolism; (bd) measurements of proteins in network, blue represents the proteases in GSH metabolism; red represents GSTA1 and GPX1.
Ijms 17 00318 g003aIjms 17 00318 g003b
Figure 4. Toxicity mechanism of PAs in glutathione metabolism: downward red arrows stand for the decrease of the activities, upward red arrows for the increase of the activities, black arrows for the direction of biological reactions and red double lines mean this pathway is blocked.
Figure 4. Toxicity mechanism of PAs in glutathione metabolism: downward red arrows stand for the decrease of the activities, upward red arrows for the increase of the activities, black arrows for the direction of biological reactions and red double lines mean this pathway is blocked.
Ijms 17 00318 g004
Table 1. Evaluation of the fourth model.
Table 1. Evaluation of the fourth model.
Coefficients a
ModelParameterUnstandardized CoefficientsStandardized CoefficientstSig.
BStd. ErrorBeta
1(Constant)8.47 × 10170.100 0.0001.000
GSTA10.8460.1170.7777.2070.000
2(Constant)0.1190.079 1.5160.139
GSTA10.7960.0890.7318.9290.000
GR−0.5970.115−0.423−5.1710.000
3(Constant)−0.0670.067 −1.0050.322
GSTA10.4710.0880.4335.3400.000
GR−0.7670.090−0.544−8.5260.000
GPX10.5520.1010.4885.4520.000
4(Constant)−0.0780.063 −1.2340.227
GSTA10.5460.0890.5026.1070.000
GR−0.9620.121−0.683−7.9640.000
GPX10.8580.1660.6975.1830.000
LAP3−0.3640.161−0.334−2.2610.031
a Dependent Variable: TOXIC.
Table 2. Prediction precision of different targets pattern.
Table 2. Prediction precision of different targets pattern.
Targets PatternToxic Prediction PrecisionNon-Toxic Predictive PrecisionTotal Prediction Precision
GSTA122/2210/1391.43%
GPX122/227/1382.86%
LAP30/224/1311.43%
GR22/225/1377.14%
LAP3, GR0/222/135.71%
GSTA1, GPX122/2211/1394.29%
GSTA1, GPX1, GR22/220/1362.86%
GSTA1, GPX1, LAP3, GR0/220/130.00%

Share and Cite

MDPI and ACS Style

Yan, X.; Kang, H.; Feng, J.; Yang, Y.; Tang, K.; Zhu, R.; Yang, L.; Wang, Z.; Cao, Z. Identification of Toxic Pyrrolizidine Alkaloids and Their Common Hepatotoxicity Mechanism. Int. J. Mol. Sci. 2016, 17, 318. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030318

AMA Style

Yan X, Kang H, Feng J, Yang Y, Tang K, Zhu R, Yang L, Wang Z, Cao Z. Identification of Toxic Pyrrolizidine Alkaloids and Their Common Hepatotoxicity Mechanism. International Journal of Molecular Sciences. 2016; 17(3):318. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030318

Chicago/Turabian Style

Yan, Xinmiao, Hong Kang, Jun Feng, Yiyan Yang, Kailin Tang, Ruixin Zhu, Li Yang, Zhengtao Wang, and Zhiwei Cao. 2016. "Identification of Toxic Pyrrolizidine Alkaloids and Their Common Hepatotoxicity Mechanism" International Journal of Molecular Sciences 17, no. 3: 318. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop