- 1College of Pharmaceutical Sciences, Zhejiang University of Technology, Hangzhou, China
- 2School of Life Sciences, Taizhou University, Taizhou, China
Background: Non-alcoholic fatty liver disease (NAFLD) is a burgeoning health problem but no drug has been approved for its treatment. Animal experiments and clinical trials have demonstrated the beneficial of saffron on NAFLD. However, the bioactive ingredients and therapeutic targets of saffron on NAFLD are unclear.
Purpose: This study aimed to identify the bioactive ingredients of saffron responsible for its effects on NAFLD and explore its therapy targets through network pharmacology combined with experimental tests.
Methods: Various network databases were searched to identify bioactive ingredients of saffron and identify NAFLD-related targets. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment were conducted to enrich functions and molecular pathways of common targets and the STRING database was used to establish a protein-protein interaction network (PPI). The effect of crocetin (CCT) on NAFLD was evaluated in a mouse model of NAFLD by measuring the biomarkers of lipid, liver and renal function, oxidative stress, and inflammation. Liver histopathology was performed to evaluate liver injury. Nuclear factor erythroid-related factor (Nrf2) and hemeoxygenase-1 (HO-1) were examined to elucidate underlying mechanism for the protective effect of saffron against NAFLD.
Results: A total of nine bioactive ingredients of saffron, including CCT, with 206 common targets showed therapeutic effects on NAFLD. Oxidative stress and diabetes related signaling pathways were identified as the critical signaling pathways mediating the therapeutic effects of the active bioactive ingredients on NAFLD. Treatment with CCT significantly reduced the activities of aspartate aminotransferase (AST), alanine transaminase (ALT), and the levels of total cholesterol (TC), triglyceride (TG), malondialdehyde (MDA), blood urea nitrogen (BUN), creatinine (CR), and uric acid (UA). CCT significantly increased the activities of superoxide dismutase (SOD), and catalase (CAT). Histological analysis showed that CCT suppressed high-fat diet (HFD) induced fat accumulation, steatohepatitis, and renal dysfunctions. Results of ELISA assay showed that CCT decreased the expression of tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), interleukin-1β (IL-1β), and increased the expression of HO-1 and Nrf2.
Conclusion: This study shows that CCT is a potential bioactive ingredient of saffron that treats NAFLD. Its mechanism of action involves suppressing of oxidative stress, mitigating inflammation, and upregulating Nrf2 and HO-1 expression.
Introduction
Non-alcoholic fatty liver disease (NAFLD) is a clinical syndrome characterized by lipid accumulation in the hepatocyte, steatosis, cellular injury, and hepatic infiltration of inflammatory cells (1, 2). The global prevalence of NAFLD is estimated to be ~25% (Younossi et al., 2018). In recent years, NAFLD pathogenesis has been explained using the “one-hit,” “two-hit,” and “multiple-hit” hypotheses (3). Traditionally, a high-fat diet (HFD) was considered an important factor in NAFLD development (4). However, growing evidence has revealed that other factors such as race, ethnicity, gender, age, and genetic pre-disposition are associated with NAFLD (5–7). Further elucidation of the molecular mechanism of NAFLD has revealed that oxidative stress, endoplasmic reticulum stress, perturbation of autophagy, mitochondrial dysfunction, hepatocellular apoptosis, and inflammatory responses cause liver damage in NAFLD (8). Among them, the role of oxidative stress in the pathogenesis of NAFLD has attracted significant attention (9). Previous studies have suggested multiple molecular pathways related to oxidative stress that is involved in NAFLD, including IRs-1/Akt, CYP2E1/JNK, SIRT1/SIRT3-FOXO3a, A(10)MPK/Nrf2, Nrf2/HO-1, and AMPK signaling pathways (10–13).
Herbal medicine has been used in NAFLD treatment, specifically on multiple mechanisms underlying lipid metabolism and inflammation (14). Saffron is one of the most expensive herbs in the world, obtained from dried stigmas of Crocus sativus L. (15). Islamic traditional medicine literature introduced saffron as an astringent, resolvent, and concoctive drug, with therapeutic effects on gastro-hepatoprotective, urogenital disorders, antidepressants, ocular disorders, etc. (16). The biological activity of saffron has been widely studied, specifically focusing on its immunomodulatory, anti-inflammatory, and antioxidant activities, and in the treatment of cardiovascular and depression-anxiety (17–20). Besides, recent research demonstrates a hepatoprotective effect of saffron in diabetic rats (21). Results from clinical trials have indicated that saffron supplementation has beneficial effects on the serum levels of inflammatory, oxidative stress, and adipokines biomarkers in patients with NAFLD (22). Existing research provides evidence on the effect of saffron in the treatment of NAFLD, however, further studies are needed to better elucidate the main bioactive ingredients and signaling pathways.
In the current study, network pharmacology was used to predict the main bioactive ingredients, potential therapeutic targets, and key pathways responsible for the effect of saffron in the treatment of NAFLD. Subsequently, the effects and possible signaling pathways of CCT on NAFLD treatment were evaluated in HFD fed mice. The schematic diagram of this study is shown in Figure 1.
Materials and Methods
Herb Materials, Reagents, and Animals
Saffron samples were collected from Shenghuo Agricultural Development Co., Ltd (Hangzhou, China). CCT reference substance was purchased from Chengdu Desite Biotechnology Co., Ltd (the content of CCT ≥ 98%, Chengdu, China). The total cholesterol (TC) assay kit, Triglyceride (TG) assay kit, Urea (BUN) assay kit, Creatinine (Cr) assay kit, Uric acid (UA) assay kit, Malondialdehyde (MDA) assay kit, Catalase (CAT) assay kit, Superoxide Dismutase (SOD) assay kit, Oil Red O dye, Haematoxylin, and Eosin dye were purchased from Jiancheng Biotechnology Institute Company (Nanjing, China). RIPA lysis buffer, Bicinchonininc acid (BCA) protein determination kit, and Nuclear and Cytoplasmic Protein Extraction kit were purchased from Beyotime Company (Shanghai, China). Nuclear factor erythroid-related factor (Nrf2) transcription factor assay kit, Hemeoxygenase-1 (HO-1) ELISA kit, Tumor necrosis factor-α (TNF-α) ELISA kit, Interleukin-1β (IL-1β) ELISA kit and Interleukin-6 (IL-6) ELISA kit were purchased from Mlbio Biotechnology Company Limited (Shanghai, China). The HFD (saccharose 20%, lard 15%, cholesterol 1.2%, sodium cholate 0.2%) was purchased from Trophic Animal Feed High-Tech Co., Ltd. (Nantong, China). C57BL/6J mice (male, 5 months old, bodyweight 27 ± 2 g; Shanghai Slac Laboratory Animal CO., Ltd., SCXK 2012-0002) were housed individually in ventilated cages (five per cage) and controlled environmental conditions in a clean grade facility at Zhejiang University of Technology (temperature 22–25°C; relative humidity 60%). The mice had free access to standard laboratory chow and tap water and were maintained under normal conditions of a 12/12 h light/dark cycle. All experiments were conducted at the Animal experimental research Center of Zhejiang University of Technology. The study was approved by the Animal Protection Research Ethics Committee of Zhejiang University of Technology (20190603064).
Screening for Bioactive Ingredients in Saffron
The bioactive ingredients of Saffron were obtained from the following databases: the Traditional Chinese Medicine Systems Pharmacology (TCMSP, http://tcmspw.com/tcmsp.php) (23), TCM database @taiwan (http://tcm.cmu.edu.tw) (24), and Traditional Chinese Medicine Integrated Database (TCMID, (http://www.megabionet.org/tcmid/) (25). The parameters oral bioavailability (OB) ≥ 30% and drug-likeness (DL) ≥ 0.18 were set as thresholds to identify the potential bioactive ingredients in TCMSP (26). To obtain comprehensive information of bioactive ingredients, recent literature and relevant international standards such as ISO 3632-1:2011 were reviewed (27–30).
Correlation Analysis Between Saffron and NAFLD-Related Targets
The corresponding targets of the bioactive ingredients of saffron were collected from the TCMSP database and SwissTargetPrediction database (http://www.swisstargetprediction.ch/). The GeneCards database (http://www.genecards.org/) and DisGeNET database (https://www.disgenet.org/) were searched using the terms “NAFLD” or “non-alcoholic fatty liver disease” (UMLS CUI: C0400966) to identify the potential therapeutic targets of NAFLD. The UniProt database (https://www.uniprot.org/) was used to standardize gene information, and further to remove duplicate genes and pseudogenes. The common targets of the bioactive ingredients were manually screened and the interaction network was visualized using Cytoscape 3.7.2. software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by using the R package “clusterProfile.” The visualization and integration of enrichment results were carried out by using Enrichlot and ggplot2 R package. The above steps were completed by R software 3.6.2 (x64) (31). Data were presented using bar charts and a bubble chart. P-value < 0.05 or Q-value < 0.05 were considered statistically significant. The construction of the protein-protein interaction (PPI) network for the common targets was generated through the String database (https://string-db.org/cgi/input.pl). Subsequently, the parameters of the PPI network, including the node degree (Degree), closeness centrality (CC), and betweenness centrality (BC) were calculated using the Network Analyzer tool function of the String database, and the results presented in a bar chart and 3D scatter plot.
Extraction and Purification of CCT
The extraction and purification of CCT from saffron were performed as reported in our previous research, but with slight modifications (32). Ultrasonic technology was applied in the extraction of CCT and 65% ethanol was used as the extracting solvent. The optimal extraction parameters were as follows: extraction temperature was 50°C, a 1:14 raw material to liquid ratio and ultrasonic power of 100 W. Extraction was performed 3 times for 1 h each. The D101 macroporous resin was used in extract purification and the purification parameters were as follows: eluents were water, 25% ethanol and 60% ethanol, elution volume of 8-bed volumes. The water and 25% ethanol eluent were discarded, while the 60% ethanol eluent was collected. NaOH solution (2 M) was used to adjust the pH of the solution to about 12, and the mixture was allowed to stand at room temperature for 12 h. H2SO4 solution (1 M) was used to adjust the pH to about 2, and the mixture precipitated overnight at 4°C. The sediment, which was a crude extract of CCT was collected by centrifugation. The crude extract was washed twice using methanol and 0.5% H2SO4 and washed twice with double-distilled water. Finally, the purified CCT was obtained by centrifugation and dried overnight in a vacuum oven at 60°C. The purity and structure of CCT were elucidated using HPLC (Figures 2A,B), 1H-NMR (Figure 2C) (600 MHz) and (-) ESI-MS/MS (Figure 2D) as shown in Figure 2.
Figure 2. The purity and structure analysis of CCT. (A) HPLC chromatogram of CCT reference substance. (B) HPLC chromatogram of CCT. (C) 1H-NMR of CCT. (D) (-) ESI-MS/MS of CCT.
Animal Treatment and Sample Collection
After 7 days of adaptive feeding, mice were randomly assigned to seven groups (n = 10 in each group) as follows: (1) Mice in the control group (control) were fed normal chow diets and given 0.5% sodium carboxyl methyl cellulose (CMC-Na) daily by oral gavage; (2) The high-fat diet group (HFD) mice were fed with HFD and given 0.5% CMC-Na daily by oral gavage; (3) The biphenyldimethylesterate group (DDB) mice were fed with HFD and given 200 mg/kg of biphenyldimethylesterate daily by oral gavage; (4) The Essentiale N group (Essential N) mice were fed with HFD and given 200 mg/kg of Essentiale N daily by oral gavage; (5) The Low-dose CCT group (CCT-L) mice were fed with HFD and given 10 mg/kg of CCT daily by oral gavage; (6) The Medium-dose CCT group (CCT-M) mice were fed with HFD and given 30 mg/kg of CCT daily by oral gavage; (7) The High-dose CCT group (CCT-H) mice were fed with HFD and given 50 mg/kg of CCT daily by oral gavage. DDB and Essential N, two drugs previously shown to be effective in the treatment of NAFLD and liver injury, were used as positive controls during the validation experiment (33–35).
The body weight was recorded and blood samples were collected from the mice orbit after an overnight fast at 2, 4, 6, and 10 weeks. All blood samples were centrifuged (Thermo, American) at 3,500 rpm for 15 min and the supernatant was collected and stored at −80°C for further analysis. Mice were isolated into single cages for 4 h for collection of feces 10 weeks after the last gavage dose. The mice were sacrificed to collect the liver tissues. The liver weight was measured and the liver index calculated using the following formula: liver index = liver weight/body weight ×100% (g/g). Each of the liver samples was divided in two; one portion was frozen and stored at −80°C for further analysis and the other portion was fixed in 10% neutralized formalin for subsequent paraffin embedding and histological analysis.
Biochemistry Analysis
Fasting serum AST, ALT, TC, TG, BUN, CR, and UA were analyzed using an automatic biochemical analyzer (model 7600; Hitachi Ltd., Japan). Feces and liver samples were placed in normal saline and homogenized using a tissue homogenizer (Scientz-48, Ningbo, China) for 60 s at 60 Hz. The supernatants were obtained by centrifugation for TC and TG analysis. MDA, CAT, and SOD in serum and homogenized liver tissue were analyzed using commercial assay kits according to the manufacturer's instructions.
Histological Analysis
Haematoxylin and eosin staining (H&E) and Oil Red O staining were used to evaluate liver histopathology. Briefly, liver specimens fixed in 10% neutral formalin were embedded in paraffin, sliced at 5 μm thickness, and stained with haematoxylin and eosin. The NAFLD activity score (NAS) ranged from 0 to 8 was used to evaluate histological liver damage. More specifically, NAS includes three histological scores as follows: steatosis (0–3), lobular inflammation (0–3), and ballooning degeneration (0–2) (36). Hepatic lipid accumulation was determined using Oil Red O staining. The images were captured using an inverted optical microscope (Olympus IX81, Japan, magnification, ×200). Isopropyl alcohol was used to elute Oil red O stain and quantification was performed at 535 nm to determine the degree of lipid deposition. Data were expressed as the fold change to control.
Inflammation Biochemical Assays
The levels of TNF-α, IL-1β, and IL-6 in the liver of each group were analyzed using commercial enzyme-linked immunosorbent assay (ELISA) kits according to the manufactures' instructions. The standard curve was used to calculate the concentration.
Enzyme-linked Immunosorbent Assay of Nrf2 and HO-1Nrf2 and HO-1 were determined by an enzyme-linked immunosorbent assay as described below. Part of the liver was first taken for total protein extraction by RIPA lysis buffer and the total protein concentration was determined by BCA protein determination kit. Nuclear extracts were isolated with a Nuclear and Cytoplasmic Protein Extraction kit and further used for the determination of Nrf2 and HO-1 by using commercial kits according to the manufactures' instructions. The standard curve was used to calculate the concentration.
Statistical Analysis
All data in this study were presented as means ± SEM and analyzed using SPSS 17.0 statistical software. One-way ANOVA was used to statistically analyze the data. P < 0.05 was considered to be statistically significant (37).
Results
Bioactive Ingredients of Saffron
A total of 70 bioactive ingredients were retrieved from TCMSP, TCM @taiwan, and TCMID databases. However, according to the screening criteria of OB (≥30%) and DL (≥0.18), 5 bioactive ingredients were identified, including quercetin, kaempferol, isorhamnetin, CCT, and n-heptanal. Besides, another four bioactive ingredients were identified from the literature and included in the final results for subsequent analysis and included, crocin I, crocin II, safranal, and picrocrocin (27–30). Detailed information on the bioactive ingredients is listed in Table 1. n-heptanal was not included in the follow-up analysis because of its unclear structure and chemical abstracts service (CAS) number.
NAFLD-Related and Bioactive Ingredients of Saffron Targets
A total of 167 corresponding targets of the bioactive ingredients of saffron were extracted from TCMSP. However, two similar bioactive ingredients were identified as having similar targets in the SwissTargetPrediction database, so the number of targets extracted in the SwissTargetPrediction database was up to 800. A total of 485 targets of the bioactive ingredients of saffron were identified after merging and removing the duplicates using two databases. A total of 1,058 NAFLD-related targets with a cut-off value larger than 20 were obtained from the DisGeNet database, while 1,017 targets were obtained from the Genecards database. Finally, a total of 1,654 NAFLD-related targets were identified after merging and removing the duplicate records from the two databases. After the manual screening, the common targets of the bioactive ingredients of saffron and NAFLD-related were limited to 206 as shown in Figure 3A. More detailed information of 206 common targets was listed in Supplementary Table 1. Thereafter, the bioactive ingredients targets-NAFLD-related targets network was constructed as shown in Figure 3B. The bioactive ingredients of saffron were found to interfere with multiple therapeutic targets related to NAFLD. The number of NAFLD-related targets associated with quercetin, kaempferol, safranal, crocin I, crocin II, and CCT was 106, 63, 49, 41, 39, and 38, respectively, while the total number of NAFLD-related targets associated with crocin I, crocin II and CCT was 118. These findings indicated that CCT and its glucosyl ester derivatives (crocin I and crocin II) play an important role in the treatment of NAFLD.
Figure 3. (A) Venn diagram of the targets of Saffron in NAFLD. (B) A network of bioactive ingredients and NAFLD-related targets. The red node represent saffron and the eight yellow nodes represent the bioactive ingredients in saffron; the green nodes represent common targets of the bioactive ingredients of saffron and NAFLD-related. The edges denote that nodes can interact with each other.
GO and KEGG Enrichment Analysis
The GO enrichment analysis of 206 common targets was shown in Figure 4A. The key biological processes implicated in the treatment of NAFLD, included reactive oxygen species, metabolic process, response to molecules of bacterial origin, cellular response to chemical stress, response to lipopolysaccharide, response to an antibiotic, response to a steroid hormone, etc. In terms of cellular components, they were involved in membrane raft, membrane microdomain, membrane region, vesicle lumen, secretory granule lumen, and cytoplasmic vesicle lumen, etc. Additionally, nuclear receptor activity, ligand-activated transcription factor activity, steroid hormone receptor activity, heme binding, phosphatase binding, and protein phosphatase binding were associated with the molecular functions in the treatment of NAFLD. Figure 4B was a simplified bubble chart showing KEGG enrichment, which showed that the underlying mechanism of NAFLD mainly involved diabetes-related pathways, apoptosis, oxidative stress, and inflammatory pathways. For clearer presentation, the specific position and function of common targets are colored red in the signaling pathway in Figure 4C and the specific targets in signaling pathways were listed in Supplementary Table 2.
Figure 4. GO and KEGG enrichment analysis. (A) GO enrichment analysis of the 206 common targets associated with NAFLD. The X -axis represents significant enrichment in the counts of these terms; The Y-axis showing the categories of biological process (BP), cellular component (CC), and molecular function (MF) in the GO of the target genes (P < 0.05). (B) KEGG enrichment analysis. The X-axis showing the counts of the target symbols in each pathway; the Y-axis showing the main pathways (P < 0.05). (C) The specific position and function of common targets in in signaling pathways. The red region represents the specific position of these common targets in the pathways.
PPT Network of Common Targets
The PPI network was constructed based on the NAFLD-related and bioactive ingredients of saffron targets. As shown in Figure 5A, the PPI network contained 184 nodes and 1,594 edges. The light blue edges and lavender edges represent the known interactions from curated databases and experimentally determined, respectively. The green edges, red edges, and dark blue edges represent predicted interactions with the gene neighborhood, gene fusions, and gene co-occurrence, respectively. The yellow edges, black edges, and light blue edges represent predicted interactions from textmining, co-expresson, and protein homology, respectively. Figure 5B showed that the top 10 core targets were AKT1, MAPK1, STAT3, PIK3CA, PIK3R1, RELA, TNF, APP, JUN, and HSP90AA1. Details of the PPI network are shown in Supplementary Table 3.
Figure 5. (A) Protein-protein interaction network. (B) The bar plot of the protein-protein interaction network. The X-axis represents the number of neighboring proteins of the target protein. The Y-axis represents the target protein.
Effects of CCT on Body Weight, Liver Weight, and Liver Index
Compared with the control group, the body weight gain of the HDF group significantly increased at week 10. CCT was found to cause a reduction in body weight gain induced by HFD. In HFD mice, the liver weight was significantly increased. Except for the Essentiale N group, the liver index in all treatment groups significantly decreased compared with the HDF group (Figures 6A,B). The liver of the HDF group was canary yellow and its surface had a distinct white dot shape representing hepatic steatosis. In contrast, the liver of the DDB group, CCT-M group, and CCT-H group was significantly improved as shown in Figure 6C.
Figure 6. (A) Body weight and liver weight. (B) Liver index. (C) Liver photographs. ##P < 0.01 vs. Control group; *P < 0.05, **P < 0.01 vs. HFD group.
Effects of CCT on Lipid, Liver, and Renal Function Biomarkers
As shown in Figure 7, the levels of serum and liver TC and TG were significantly increased in the HFD group, while were significantly decreased in all treatment groups compared with the HDF group (Figures 7A,B). The levels of fecal TC and TG were significantly increased in three dose groups of CCT compared with the HDF group (Figure 7C). These findings indicated that CCT can increase lipid excretion to alleviate lipid accumulation caused by NAFLD. The activities of ALT and AST, two markers of liver function, were also significantly increased in the HDF group compared with the control group. The effect of CCT in reducing ALT and AST activities was observed after treatment, and also in the DBB and Essential N groups (Figure 7G). For renal function, the serum levels of BUN, CR, and UA were significantly increased in the HDF group compared with the control group. However, a significant decrease was observed in three dose groups of CCT compared with the HDF group (Figures 7D–F).
Figure 7. Effects of CCT on lipid, liver and renal function biomarkers. (A) Effects of CCT on serum TC and TG. (B) Effects of CCT on liver TC and TG. (C) Effects of CCT on fecal TC and TG. (D) Effects of CCT on BUN. (E) Effects of CCT on CR. (F) Effects of CCT on UA. (G) Effects of CCT on serum AST and ALT. ##P < 0.01 vs. Control group; **P < 0.01 vs. HFD group.
Effects of CCT on Biomarkers of Oxidative Stress
As shown in Figure 8, the serum and liver tissue MDA levels were significantly increased in the HDF group compared with the control group, and there was a significant decrease in the three-dose groups of CCT and DDB group compared with the HDF group (Figure 8A). Similar trends were observed in the activities of SOD and CAT indicating that CCT can alleviate oxidative stress induced by NAFLD (Figures 8B,C).
Figure 8. Effects of CCT on oxidative stress biomarkers. (A) Effects of CCT on serum MDA and liver TG. (B) Effects of CCT on serum SOD and liver SOD. (C) Effects of CCT on serum CAT and liver CAT.##P < 0.01 vs. Control group; *P < 0.05, **P < 0.01 vs. HFD group.
Effects of CCT on Biomarkers of Inflammation
To determine whether CCT could inhibit inflammation in an HFD-induced NAFLD-mouse model, the levels of TNF-α, IL-1β, and IL-6 were measured in the liver tissue. As shown in Figure 9, TNF-α, IL-1β, and IL-6 levels were significantly increased in the HDF group compared with the control group, but significantly decreased in CCT groups, the DDB group and Essential N group.
Figure 9. Effects of CCT on inflammation biomarkers. (A) Effects of CCT on TNF-α. (B) Effects of CCT on IL-1β. (C) Effects of CCT on IL-6. ##P < 0.01 vs. Control group; **P < 0.01 vs. HFD group.
Effects of CCT on Liver Histopathology
In the present study, 10 photomicrographs of HE-staining were selected for use in calculating the NAS score at a magnification of ×200. The NAS score in the HFD group was significantly increased compared with the control group and decreased after CCT treatment. This indicated that CCT significantly attenuated inflammation, steatosis, and ballooning (Figures 10A,B). Besides, Oil Red O staining (Figures 10C,D) showed increased accumulation of lipid droplets in the liver of the HFD group, however, CCT treatment led to a significant decrease in the lipid droplets in the liver.
Figure 10. Effects of CCT on liver histopathology. (A) Representative photomicrographs of HE staining (×200). The blue arrow represents steatosis, the green arrow represents lobular inflammation, the black arrow represents ballooning degeneration. (B) The NAS scores. (C) Representative photomicrographs of Oil Red O staining (×200). (D) The OD of Oil Red O staining(fold change). ##P < 0.01 vs. Control group; *P < 0.05, **P < 0.01 vs. HFD group.
Effects of CCT on Nrf2 and HO-1
As shown in Figure 11, the expression of Nrf2 and HO-1 significantly decreased in the HDF group compared with the control group. After CCT treatment, the expression of Nrf2 and HO-1 significantly increased in all the groups, except for the DDB group compared with the HDF group.
Figure 11. Effects of CCT on Nrf2 and HO-1. (A) Effects of CCT on Nrf2. (B) Effects of CCT on HO-1. ##P < 0.01 vs. Control group; *P < 0.05, **P < 0.01 vs. HFD group.
Discussion
Network pharmacology is a newly developing interdisciplinary science based on biology, chemical informatics, and bioinformatics, and is also an effective tool used to explore bioactive ingredients from complex components (38). In recent years, network pharmacology has been widely used in traditional Chinese medicine to investigate the potential therapeutic targets and pharmacological mechanisms in different diseases, however, very few studies have used network pharmacology to study the pharmacological effects of saffron. Consistent with our results, a comparative study of the anti-thrombotic effects of saffron and carthami flos based on network pharmacology indicated that CCT was the key bioactive ingredient of saffron (39). Two other studies on systems pharmacology of saffron emphasized the importance of CCT (40, 41). CCT has various biological activities such as anti-tumor, neuroprotective, cardioprotective, hepatoprotective, antidepressant, and improves asthma, Alzheimer's, and diabetes mellitus, etc. (42). In studies of liver diseases, CCT showed comparative beneficial effects on CCL4-induced and dengue virus-induced liver damage via induction of antioxidant defense (43, 44). In this study, the therapeutic targets of crocin I, crocin II, and CCT were combined, because orally administered crocins are hydrolyzed to CCT before they are incorporated into the blood circulation, and these findings have been confirmed in several pharmacokinetic studies (45–47). As noted above, CCT was reasonable to reckon to be the critical bioactive ingredient of saffron and further used in the verification experiment.
Although several decades of NAFLD research have resulted in major scientific advancements in pathogenesis and therapeutic target, effective drugs are currently under clinical development (48). The main challenge is that the efficacy of a single target is limited because the pathogenesis and progression of NAFLD involve multiple pathogenic pathways, such as insulin resistance, lipotoxicity, oxidative stress, altered immune/cytokine/mitochondrial functioning, and apoptosis (48–50). In the current study, the possible signaling pathways of saffron intervention in NAFLD were analyzed according to the results obtained from GO, KEGG, and PPI analysis. The results of GO enrichment indicated that saffron may interfere with oxidative stress and nuclear receptor activity in the treatment of NAFLD. The results of KEGG enrichment indicated multiple signaling pathways, including AGE-RAGE signaling pathway in diabetic complications, HIF−1, TNF, adipocytokine, NF-kappa B, PPAR, and other diabetes-related signaling pathways involved in the treatment of NAFLD by using saffron. However, the PPI analysis gave very different results, incicating that AKT1, MAPK1, STAT3, PIK3CA, PIK3R1, RELA, TNF, APP, and JUN were the top 10 core targets. In view of the inconsistent results from KEGG enrichment and PPI, it is necessary to systematically analyze the targets and signaling pathways of saffron in the treatment of NAFLD. Taking into consideration that the results of oxidative stress and nuclear receptor activity from GO enrichment, the 206 common targets were reevaluated, of which NFE2L2 (Nrf2) and HMOX1 (HO-1) have attracted considerable attention. Nrf2 is a tran-scription factor that regulates the expression of several anti-oxidant genes. Previous studies show that Nrf2/HO-1-antioxidant response element (ARE)-antioxidant enzyme play a central mechanistic role in the regulation of the complex antioxidant system in the human body and are key nodes in multiple signaling pathways, including PI3K/AKT/MAPK and IL-6/JAK/STAT3 (51–54). Under normal physiological condition, Kelch-like-ECH-associated protein 1 (Keap1), an Nrf2 repressor, binds to Nrf2 to form a complex and resides in the cytoplasm. However, under stimulation with oxidative and electrophilic chemical signals, Nrf2 is released from Keap1 and transferred to the nucleus, where it binds to the ARE. Here, Keap1 plays an important role in cellular responses to chemical and oxidative stress by regulating the stability and nuclear translocation of Nrf2 protein (55, 56). As one of the key target genes of nuclear Nrf2, HO-1 is a downstream antioxidant enzyme and has a broad cytoprotective effect in various diseases (57). It can be inferred from the above results that CCT may play a role in the treatment of NAFLD by co-intervening the expression of Nrf2 and HO-1. The evidences in biomedical literature in support of this hypothesis have been identified. The antioxidant effect of scutellarin on NAFLD is dependent on PI3K/AKT activation with subsequent Nrf2 nuclear translocation, which increases the expression of HO-1 (58). Resolvin D1 mitigates non-alcoholic steatohepatitis by suppressing the TLR4-MyD88-mediated NF-kappa B and MAPK pathways and activates the Nrf2 pathway in mice (59). In view of the above considerations, Nrf2 and HO-1 were identified as key target genes of the possible signaling pathway and further validated in vivo in the current study.
In the validation experiment, HDF was used to induce NAFLD mice to further ascertain whether CCT can ameliorate NAFLD. The results showed that CCT regulates the levels (activities) of TC, TG ALT, AST, and liver index. In the diagnosis of NAFLD, pathological biopsy remains the “golden standard” (60). H&E stained liver sections revealed that CCT can ameliorate hepatocyte ballooning, steatosis, and inflammation to prevent NAFLD. Lipid droplets in the liver were stained using Oil red O, and CCT decreased the lipid droplets. Recent studies have found that NAFLD not only leads to abnormal liver function but also causes extrahepatic disease, such as renal damage (61, 62). Liver steatosis alters the gut barrier function and microbial composition to accumulate toxic metabolites resulting in renal damage (63, 64). In this study, the serum levels of BUN, CR, and UA were analyzed to determine whether CCT can ameliorate renal function caused by NAFLD. The results showed that feeding mice with HFD significantly increased the levels of BUN, CR, and UA, which is consistent with previous clinical trials (65). After CCT treatment, BUN, CR, and UA significantly decreased which was associated with efficient glomerular filtration function and improved NAFLD (66). MDA, a product of lipid oxidation, is used as a marker of oxidative stress. SOD and CAT are enzymatic antioxidant systems, which play important roles in protection against the deleterious effects of hydrogen peroxide and lipid peroxidation in diseases related to oxidative stress (67, 68). In this study, a significant increase in SOD and CAT activities, accompanied by a decrease in MDA level was observed after CCT treatment. These results indicate that CCT can reduce the levels of oxidative stress in the HFD-induced NAFLD model. There is overwhelming evidence that the release of adipocytokines like TNF-α, IL-1β, and IL-6 result in hepatocytes suffering a characteristic variation in their structure developing lobular inflammation and balloon degeneration associated with different degrees of scarring or fibrosis (69). In this study, CCT was found to significantly reduce the levels of TNF-α, IL-1β, and IL-6 which exerts a beneficial effect on NAFLD. After CCT treatment, Nrf2 and HO-1 expression was significantly increased indicating that CCT may activate multiple anti-oxidative signaling pathways to protect against liver lipid peroxidation and hepatitis.
Conclusion
In conclusion, this study explored, for the first time, the bioactive ingredients of saffron and their therapeutic targets in NAFLD using network pharmacology and animal experiments. CCT was identified as the bioactive compound that ameliorates HFD-induced NAFLD by decreasing level of antioxidants and proinflammatory factors. This study provides a scientific basis for further analysis of the clinical application of CCT in NAFLD.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Ethics Statement
All experiments were conducted at the Animal experimental research Center of Zhejiang University of Technology. The study was approved by the Animal Protection Research Ethics Committee of Zhejiang University of Technology (20190603064).
Author Contributions
PW and YT supervised the project, conceived, and designed the research methods. ZX, SL, QL, and YY performed the experiments. YJ and PF performed network pharmacology analysis. ZX wrote the manuscript. YC revised the manuscript. All authors reviewed and approved the final manuscript.
Funding
This work was financially supported by the Special Project of International Technology Cooperation of One Belt and One Road (Grant No. 2017C04009); the key projects of international scientific and technological innovation cooperation between governments (Grant No. 2017YFE0130100).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.681391/full#supplementary-material
Abbreviations
ALT, alanine transaminase; AST, aspartate aminotransferase; BUN, blood urea nitrogen; CAT, catalase; CCT, crocetin; CR, creatinine; DL, drug-likeness; DAVID, visualization and Integrated Discovery; GO, Gene ontology; H&E, haematoxylin and eosin staining; HFD, high-fat diet; HO-1, hemeoxygenase-1; IL-1β, interleukin-1β; IL-6, interleukin-6; KEGG, Kyoto Encyclopedia of Genes and Genomes; MDA, malondialdehyde; NAFLD, non-alcoholic fatty liver disease; NAS, NAFLD activity score; Nrf2, nuclear factor erythroid-related factor; OB, oral bioavailability; PPI, protein-protein interaction network; SOD, superoxide dismutase; TC, total cholesterol; TCMSP, traditional Chinese Medicine Systems Pharmacology; TCMID, traditional Chinese Medicine Integrated Database; TG, triglyceride; TNF-α, tumor necrosis factor-α; UA, uric acid.
References
1. Pappachan JM, Babu S, Krishnan B, Ravindran NC. Non-alcoholic fatty liver disease: a clinical update. J Clin Transl Hepatol. (2017) 5:384–93. doi: 10.14218/JCTH.201700013
2. Saeed A, Dullaart R, Schreuder T, Blokzijl H, Faber KN. Disturbed vitamin a metabolism in non-alcoholic fatty liver disease (NAFLD). Nutrients. (2017) 10:29. doi: 10.3390/nu10010029
3. Buzzetti E, Pinzani M, Tsochatzis EA. The multiple-hit pathogenesis of non-alcoholic fatty liver disease (NAFLD). Metabolism. (2016) 65:1038–48. doi: 10.1016/j.metabol.2015.12012
4. Lian CY, Zhai ZZ, Li ZF, Wang L. High fat diet-triggered non-alcoholic fatty liver disease: A review of proposed mechanisms. Chem Biol Interact. (2020) 330:109199. doi: 10.1016/j.cbi.2020109199
5. Ballestri S, Nascimbeni F, Baldelli E, Marrazzo A, Romagnoli D, Lonardo A. NAFLD as a sexual dimorphic disease: role of gender and reproductive status in the development and progression of nonalcoholic fatty liver disease and inherent cardiovascular risk. Adv Ther. (2017) 34:1291–326. doi: 10.1007/s12325-017-0556-1
6. Gong Z, Tas E, Yakar S, Muzumdar R. Hepatic lipid metabolism and non-alcoholic fatty liver disease in aging. Mol Cell Endocrinol. (2017) 455:115–30. doi: 10.1016/j.mce.2016.12022
7. Hotta K, Kitamoto T, Kitamoto A, Ogawa Y, Honda Y, Kessoku T, et al. Identification of the genomic region under epigenetic regulation during non-alcoholic fatty liver disease progression. Hepatol Res. (2018) 48:E320–34. doi: 10.1111/hepr12992
8. Berlanga A, Guiu-Jurado E, Porras JA, Auguet T. Molecular pathways in non-alcoholic fatty liver disease. Clin Exp Gastroenterol. (2014) 7:221–39. doi: 10.2147/CEGS62831
9. Dallio M, Sangineto M, Romeo M, Villani R, Romano AD, Loguercio C, et al. Immunity as cornerstone of non-alcoholic fatty liver disease: the contribution of oxidative stress in the disease progression. Int J Mol Sci. (2021) 22:436. doi: 10.3390/ijms22010436
10. Tong W, Ju L, Qiu M, Xie Q, Chen Y, Shen W, et al. Liraglutide ameliorates non-alcoholic fatty liver disease by enhancing mitochondrial architecture and promoting autophagy through the SIRT1/SIRT3-FOXO3a pathway. Hepatol Res. (2016) 46:933–43. doi: 10.1111/hepr12634
11. Qu LL, Yu B, Li Z, Jiang WX, Jiang JD, Kong WJ. gastrodin ameliorates oxidative stress and proinflammatory response in nonalcoholic fatty liver disease through the AMPK/Nrf2 pathway. Phytother Res. (2016) 30:402–11. doi: 10.1002/ptr5541
12. Jian T, Ding X, Wu Y, Ren B, Li W, Lv H, et al. Hepatoprotective effect of loquat leaf flavonoids in PM2.5-induced non-alcoholic fatty liver disease via regulation of IRs-1/Akt and CYP2E1/JNK pathways. Int J Mol Sci. (2018) 19:3005. doi: 10.3390/ijms19103005
13. Shen B, Zhao C, Wang Y, Peng Y, Cheng J, Li Z, et al. Aucubin inhibited lipid accumulation and oxidative stress via Nrf2/HO-1 and AMPK signalling pathways. J Cell Mol Med. (2019) 23:4063–75. doi: 10.1111/jcmm14293
14. Xu Y, Guo W, Zhang C, Chen F, Tan HY, Li S, et al. Herbal medicine in the treatment of non-alcoholic fatty liver diseases-efficacy, action mechanism, clinical application. Front Pharmacol. (2020) 11:601. doi: 10.3389/fphar.202000601
15. Mollazadeh H, Emami SA, Hosseinzadeh H. Razi's Al-Hawi and saffron (Crocus sativus): a review. Iran J Basic Med Sci. (2015) 18:1153–66. Available online at: https://www.researchgate.net/profile/Hamid-Mollazadeh/publication/290497291_Razi%27s_Al-Hawi_and_saffron_Crocus_sativus_A_review/links/5699fbfa08aea1476944482e/Razis-Al-Hawi-and-saffron-Crocus-sativus-A-review.pdf
16. Javadi B, Sahebkar A, Emami SA. A survey on saffron in major islamic traditional medicine books. Iran J Basic Med Sci. (2013) 16:1–11. Available online at: https://www.researchgate.net/profile/Behjat-Javadi/publication/236604674_A_Survey_on_Saffron_in_Major_Islamic_Traditional_Medicine_Books/links/568cc03708aeb488ea2fdff0/A-Survey-on-Saffron-in-Major-Islamic-Traditional-Medicine-Books.pdf
17. Boskabady MH, Farkhondeh T. Antiinflammatory, antioxidant, and immunomodulatory effects of Crocus sativus L. and its main constituents. Phytother Res. (2016) 30:1072–94. doi: 10.1002/ptr5622
18. Ghaffari S, Roshanravan N. Saffron; an updated review on biological properties with special focus on cardiovascular effects. Biomed Pharmacother. (2019) 109:21–7. doi: 10.1016/j.biopha.2018.10031
19. Zeinali M, Zirak MR, Rezaee SA, Karimi G, Hosseinzadeh H. Immunoregulatory and anti-inflammatory properties of Crocus sativus (Saffron) and its main active constituents: a review. Iran J Basic Med Sci. (2019) 22:334–44. https://www.researchgate.net/publication/333249258_Immunoregulatory_and_anti-inflammatory_properties_of_Crocus_sativus_Saffron_and_its_main_active_constituents_A_review
20. Milajerdi A, Jazayeri S, Shirzadi E, Hashemzadeh N, Azizgol A, Djazayery A, et al. The effects of alcoholic extract of saffron (Crocus satious L.) on mild to moderate comorbid depression-anxiety, sleep quality, and life satisfaction in type 2 diabetes mellitus: a double-blind, randomized and placebo-controlled clinical trial. Complement Ther Med. (2018) 41:196–202. doi: 10.1016/j.ctim.2018.09023
21. Konstantopoulos P, Doulamis IP, Tzani A, Korou ML, Agapitos E, Vlachos IS. Metabolic effects of Crocus sativus and protective action against non-alcoholic fatty liver disease in diabetic rats. Biomed Rep. (2017) 6:513–8. doi: 10.3892/br.2017884
22. Pour FK, Aryaeian N, Mokhtare M, Mirnasrollahi PR, Jannani L, Agah S, et al. The effect of saffron supplementation on some inflammatory and oxidative markers, leptin, adiponectin, and body composition in patients with nonalcoholic fatty liver disease: a double-blind randomized clinical trial. Phytother Res. (2020) 34:3367–78. doi: 10.1002/ptr6791
23. Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. (2014) 6:13. doi: 10.1186/1758-2946-6-13
24. Sanderson K. Databases aim to bridge the East-West divide of drug discovery. Nat Med. (2011) 17:1531. doi: 10.1038/nm1211-1531a
25. Xue R, Fang Z, Zhang M, Yi Z, Wen C, Shi T. TCMID: Traditional Chinese Medicine integrative database for herb molecular mechanism analysis. Nucleic Acids Res. (2013) 41:D1089–95. doi: 10.1093/nar/gks1100
26. An L, Lin Y, Li L, Kong M, Lou Y, Wu J, et al. Integrating network pharmacology and experimental validation to investigate the effects and mechanism of astragalus flavonoids against hepatic fibrosis. Front Pharmacol. (2020) 11:618262. doi: 10.3389/fphar.2020618262
27. Kabiri M, Rezadoost H, Ghassempour A. A comparative quality study of saffron constituents through HPLC and HPTLC methods followed by isolation of crocins and picrocrocin. LWT. (2017) 84:1–9. doi: 10.1016/j.lwt.2017.05033
28. Lambrianidou A, Koutsougianni F, Papapostolou I, Dimas K. Recent advances on the anticancer properties of saffron (Crocus sativus L.) and its major constituents. Molecules. (2020) 26:86. doi: 10.3390/molecules26010086
29. Cerda-Bernad D, Valero-Cases E, Pastor JJ, Frutos MJ. Saffron bioactives crocin, crocetin and safranal: effect on oxidative stress and mechanisms of action. Crit Rev Food Sci Nutr. (2020) 24:1–18. doi: 10.1080/10408398.20201864279
30. Bononi M, Milella P, Tateo F. Gas chromatography of safranal as preferable method for the commercial grading of saffron (Crocus sativus L.). Food Chem. (2015) 176:17–21. doi: 10.1016/j.foodchem.2014.12047
31. Chen S, Xu H, Hu F, Wang T. Identification of key players involved in CoCl2 hypoxia induced pulmonary artery hypertension in vitro. Front Genet. (2020) 11:232. doi: 10.3389/fgene.202000232
32. Lin S, Li Q, Jiang S, Xu Z, Jiang Y, Liu L, et al. Crocetin ameliorates chronic restraint stress-induced depression-like behaviors in mice by regulating MEK/ERK pathways and gut microbiota. J Ethnopharmacol. (2021) 268:113608. doi: 10.1016/j.jep.2020113608
33. Dajani AI, Popovic B. Essential phospholipids for nonalcoholic fatty liver disease associated with metabolic syndrome: a systematic review and network meta-analysis. World J Clin Cases. (2020) 8:5235–49. doi: 10.12998/wjcc.v8.i215235
34. Lee MM, Kim HG, Lee SB, Lee JS, Kim WY, Choi SH, et al. CGplus, a standardized herbal composition ameliorates non-alcoholic steatohepatitis in a tunicamycin-induced mouse model. Phytomedicine. (2018) 41:24–32. doi: 10.1016/j.phymed.2018.01020
35. Morsy MA, Ibrahim MA, Abd-Elghany MI. Dimethyl dimethoxy biphenyl dicarboxylate attenuates hepatic and metabolic alterations in high fructose-fed rats. Toxicol Ind Health. (2016) 32:59–67. doi: 10.1177/0748233713498445
36. Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, et al. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. (2005) 41:1313–21. doi: 10.1002/hep20701
37. Curtis MJ, Bond RA, Spina D, Ahluwalia A, Alexander SP, Giembycz MA, et al. Experimental design and analysis and their reporting: new guidance for publication in BJP. Br J Pharmacol. (2015) 172:3461–71. doi: 10.1111/bph12856
38. Lai X, Wang X, Hu Y, Su S, Li W, Li S. Editorial: network pharmacology and traditional medicine. Front Pharmacol. (2020) 11:1194. doi: 10.3389/fphar.202001194
39. Jiang J, Lin S, Li Q, Jiang S, Hu Y, Liu L, et al. Comparative studies of the anti-thrombotic effects of saffron and HongHua based on network pharmacology. Trop J Pharm Res. (2020) 7:1441–8. doi: 10.4314/tjpr.v19i716
40. Liu J, Liu J, Shen F, Qin Z, Jiang M, Zhu J, et al. Systems pharmacology analysis of synergy of TCM: an example using saffron formula. Sci Rep. (2018) 8:380. doi: 10.1038/s41598-017-18764-2
41. Liu J, Mu J, Zheng C, Chen X, Guo Z, Huang C, et al. Systems-pharmacology dissection of traditional chinese medicine compound saffron formula reveals multi-scale treatment strategy for cardiovascular diseases. Sci Rep. (2016) 6:19809. doi: 10.1038/srep19809
42. Hashemi M, Hosseinzadeh H. A comprehensive review on biological activities and toxicology of crocetin. Food Chem Toxicol. (2019) 130:44–60. doi: 10.1016/j.fct.2019.05017
43. Chen P, Chen Y, Wang Y, Cai S, Deng L, Liu J, et al. Comparative evaluation of hepatoprotective activities of geniposide, crocins and crocetin by CCl4-induced liver injury in mice. Biomol Ther (Seoul). (2016) 24:156–62. doi: 10.4062/biomolther.2015094
44. Sreekanth GP, Chuncharunee A, Yenchitsomanus PT, Limjindaporn T. Crocetin improves dengue virus-induced liver injury. Viruses. (2020) 12:825. doi: 10.3390/v12080825
45. Asai A, Nakano T, Takahashi M, Nagao A. Orally administered crocetin and crocins are absorbed into blood plasma as crocetin and its glucuronide conjugates in mice. J Agric Food Chem. (2005) 53:7302–6. doi: 10.1021/jf0509355
46. Xi L, Qian Z, Du P, Fu J. Pharmacokinetic properties of crocin (crocetin digentiobiose ester) following oral administration in rats. Phytomedicine. (2007) 14:633–6. doi: 10.1016/j.phymed.2006.11028
47. Christodoulou E, Grafakou ME, Skaltsa E, Kadoglou N, Kostomitsopoulos N, Valsami G. Preparation, chemical characterization and determination of crocetin's pharmacokinetics after oral and intravenous administration of saffron (Crocus sativus L.) aqueous extract to C57/BL6J mice. J Pharm Pharmacol. (2019) 71:753–64. doi: 10.1111/jphp13055
48. Friedman SL, Neuschwander-Tetri BA, Rinella M, Sanyal AJ. Mechanisms of NAFLD development and therapeutic strategies. Nat Med. (2018) 24:908–22. doi: 10.1038/s41591-018-0104-9
49. Younossi ZM, Loomba R, Rinella ME, Bugianesi E, Marchesini G, Neuschwander-Tetri BA, et al. Current and future therapeutic regimens for nonalcoholic fatty liver disease and nonalcoholic steatohepatitis. Hepatology. (2018) 68:361–71. doi: 10.1002/hep29724
50. Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M, et al. Global burden of NAFLD and NASH: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol. (2018) 15:11–20. doi: 10.1038/nrgastro.2017109
51. Singh S, Vrishni S, Singh BK, Rahman I, Kakkar P. Nrf2-ARE stress response mechanism: a control point in oxidative stress-mediated dysfunctions and chronic inflammatory diseases. Free Radic Res. (2010) 44:1267–88. doi: 10.3109/10715762.2010507670
52. Hannan MA, Dash R, Sohag A, Haque MN, Moon IS. Neuroprotection against oxidative stress: phytochemicals targeting trkb signaling and the Nrf2-ARE antioxidant system. Front Mol Neurosci. (2020) 13:116. doi: 10.3389/fnmol.202000116
53. Kumar H, Kim IS, More SV, Kim BW, Choi DK. Natural product-derived pharmacological modulators of Nrf2/ARE pathway for chronic diseases. Nat Prod Rep. (2014) 31:109–39. doi: 10.1039/C3NP70065H
54. Wang DJ, Cai YQ, Pan SZ, Zhang LZ, Chen YX, Chen FM, et al. Effect of total flavone of haw leaves on nuclear factor erythroid-2 related factor and other related factors in nonalcoholic steatohepatitis rats. Chin J Integr Med. (2018) 24:265–71. doi: 10.1007/s11655-016-2450-0
55. Chai J, Luo L, Hou F, Fan X, Yu J, Ma W, et al. Agmatine reduces lipopolysaccharide-mediated oxidant response via activating PI3K/Akt pathway and up-regulating Nrf2 and HO-1 expression in macrophages. PLoS ONE. (2016) 11:e0163634. doi: 10.1371/journal.pone0163634
56. He Y, Jiang J, He B, Shi Z. Chemical activators of the Nrf2 signaling pathway in nonalcoholic fatty liver disease. Nat Prod Commun. (2021) 16:1–9. doi: 10.1177/1934578X20987095
57. Nikam A, Ollivier A, Rivard M, Wilson JL, Mebarki K, Martens T, et al. Diverse Nrf2 activators coordinated to cobalt carbonyls induce heme oxygenase-1 and release carbon monoxide in vitro and in vivo. J Med Chem. (2016) 59:756–62. doi: 10.1021/acs.jmedchem5b01509
58. Fan H, Ma X, Lin P, Kang Q, Zhao Z, Wang L, et al. Scutellarin prevents nonalcoholic fatty liver disease (NAFLD) and hyperlipidemia via PI3K/AKT-dependent activation of nuclear factor (Erythroid-Derived 2)-like 2 (Nrf2) in rats. Med Sci Monit. (2017) 23:5599–612. doi: 10.12659/MSM907530
59. Li J, Deng X, Bai T, Wang S, Jiang Q, Xu K. Resolvin D1 mitigates non-alcoholic steatohepatitis by suppressing the TLR4-MyD88-mediated NF-kappaB and MAPK pathways and activating the Nrf2 pathway in mice. Int Immunopharmacol. (2020) 88:106961. doi: 10.1016/j.intimp.2020106961
60. Tapper EB, Lok A. Use of liver imaging and biopsy in clinical practice. N Engl J Med. (2017) 377:2296–7. doi: 10.1056/NEJMra1610570
61. Mantovani A, Zusi C, Dalbeni A, Grani G, Buzzetti E. Risk of kidney dysfunction in NAFLD. Curr Pharm Des. (2020) 26:1045–61. doi: 10.2174/1381612825666191026113119
62. Targher G, Byrne CD. Non-alcoholic fatty liver disease: an emerging driving force in chronic kidney disease. Nat Rev Nephrol. (2017) 13:297–310. doi: 10.1038/nrneph.201716
63. Musso G, Cassader M, Cohney S, Pinach S, Saba F, Gambino R. Emerging liver-kidney interactions in nonalcoholic fatty liver disease. Trends Mol Med. (2015) 21:645–62. doi: 10.1016/j.molmed.2015.08005
64. Paolella G, Mandato C, Pierri L, Poeta M, Stasi Di. M, Vajro P. Gut-liver axis and probiotics: their role in non-alcoholic fatty liver disease. World J Gastroenterol. (2014) 20:15518–31. doi: 10.3748/wjg.v20.i4215518
65. Dong B, Mao Y, Li Z, Yu F. The value of the atherogenic index of plasma in non-obese people with non-alcoholic fatty liver disease: a secondary analysis based on a cross-sectional study. Lipids Health Dis. (2020) 19:148. doi: 10.1186/s12944-020-01319-2
66. Shen Z, Munker S, Luo F, Ma H, Yu C, Li Y. Effect of non-alcoholic fatty liver disease on estimated glomerular filtration rate could be dependent on age. PLoS ONE. (2015) 10:e01306146. doi: 10.1371/journal.pone0130614
67. Krolow R, Arcego DM, Noschang C, Weis SN, Dalmaz C. Oxidative imbalance and anxiety disorders. Curr Neuropharmacol. (2014) 12:193–204. doi: 10.2174/1570159X11666131120223530
68. Zhu R, Wang Y, Zhang L, Guo Q. Oxidative stress and liver disease. Hepatol Res. (2012) 42:741–9. doi: 10.1111/j.1872-034X.2012.00996x
Keywords: saffron, crocetin, NAFLD, network pharmacology, Nrf2, HO-1
Citation: Xu Z, Lin S, Gong J, Feng P, Cao Y, Li Q, Jiang Y, You Y, Tong Y and Wang P (2021) Exploring the Protective Effects and Mechanism of Crocetin From Saffron Against NAFLD by Network Pharmacology and Experimental Validation. Front. Med. 8:681391. doi: 10.3389/fmed.2021.681391
Received: 25 March 2021; Accepted: 18 May 2021;
Published: 09 June 2021.
Edited by:
Águeda González Rodríguez, Princess University Hospital, SpainReviewed by:
Kunkai Su, Zhejiang University, ChinaAnabel Fernández-Iglesias, Institut de Recerca Biomèdica August Pi i Sunyer (IDIBAPS), Spain
Copyright © 2021 Xu, Lin, Gong, Feng, Cao, Li, Jiang, You, Tong and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ping Wang, wangping45@aliyun.com; Yingpeng Tong, fish166@126.com
†These authors have contributed equally to this work and share first authorship