Jianxiong Wu
Xietao Ye
Songhong Yang
Huan Yu*
Lingyun Zhong- School of Pharmacy, Jiangxi University of Traditional Chinese Medicine, Nanchang, China
Liver diseases are mostly triggered by oxidative stress and inflammation, leading to extracellular matrix overproduction and prone to develop into liver fibrosis, cirrhosis and hepatocellular carcinoma. Liver injury (LI) refers to various pathogenic factors leading to the destruction of stem cells that then affect the liver’s normal function, causing a series of symptoms and abnormal liver function indicators. Citri Reticulatae Pericarpium (CRP) is one of the most commonly used traditional Chinese medicines; it contains flavonoids including hesperidin, nobiletin, and tangeretin. CRP has antibacterial, antioxidant, and antitumor effects that reduce cholesterol, prevent atherosclerosis and decrease LI. Here we analyzed the components of CRP and their targets of action in LI treatment and assessed the relationships between them using a systems pharmacology approach. Twenty-five active ingredients against LI were selected based on ultra-performance liquid chromatography-quadrupole/time-of-flight mass spectrometry results and databases. The drug targets and disease-related targets were predicted. The 117 common targets were used to construct a protein-protein interaction network. We identified 1719 gene ontology items in LI treatment, including 1,525 biological processes, 55 cellular components, and 139 molecular functions. These correlated with 49 Kyoto Encyclopedia of Genes and Genomes pathways. These findings suggest that CRP may counteract LI by affecting apoptotic, inflammatory, and energy metabolism modules. In vitro experiments suggested that the mechanism may involve hesperidin and naringenin acting on CASP3, BAX, and BCL2 to affect the apoptosis pathway, attenuating liver fibrosis. Naringenin significantly inhibited AKT1 phosphorylation, which in turn mediated activation of the phosphoinositide 3-kinase-Akt signaling pathways against LI. This study provides a reference for systematically exploring the mechanism of CRP’s anti-LI action and is also expands of the application of systems pharmacology in the study of traditional Chinese medicine.
Introduction
The pathogenetic events leading to liver injury (LI) include oxidative damage, apoptosis, disordered lipid metabolism, and mitochondrial dysfunction. LI refers to one or more results of organ damage, inflammation, bleeding, infection, or other symptoms. Alcohol use is a major predisposing factor for LI (Yang et al., 2020). The annual proportion of alcohol-induced LI has been increasing (Tsutsumi et al., 1994; Vatsalya et al., 2016). Excessive alcohol consumption often causes fatty acid imbalances that trigger liver inflammation.
Traditional Chinese medicine (TCM) is used to treat various diseases. Numerous studies have demonstrated that flavonoids and alkaloids components in TCM prevent or treat LI (Liu, 2016; Zhang et al., 2018). Citri Reticulatae Pericarpium (CRP, orange peel) is a dried and mature pericarp Citrus × aurantium L (Citrus × reticulata Blanco; Rutaceae) (Prc, 2015). CRP extract has excellent biological activity, and the citrus-rich extract is rich in polymethoxy flavonoids that regulate metabolism. CRP extract has excellent biological activity, and citrus-rich extracts rich in polymethoxy flavonoids reduce malnutrition (Zeng et al., 2020). CRP is used in TCM to treat nausea, vomiting, and anemia (Yu et al., 2018). Volatile oils and flavonoids in CRP are considered major components; these include narirutin, nobiletin, and hesperidin (Cui et al., 2019). These components act alone or in combination to combat inflammatory responses and lipid peroxidation, in turn attenuating oxidative stress-induced hepatocyte damage (Lin et al., 2019).
According to the TCM theory, LI is caused by damp heat; viruses invade the liver and lodge in the blood, resulting in liver loss and catharsis regulation. The spleen loses its health and transport capabilities, causing disordered qi and blood. CRP is bitter and pungent in taste, warm in nature; it regulates qi and dampness, eliminating phlegm and cough, stopping vomiting, eliminating food from the stomach, invigorating the spleen, soothing the liver, and promoting bile discharge (Prc, 2015). A typical TCM formula for treating liver disease is “Chaihu Shugan Powder,” one of CRP’s main components (Nie et al., 2020). TCM treatment reflects three points: overall, dynamic, and syndrome differentiation. Medications regulate bodily functions and balance them by eliminating pathogenic factors and strengthening vital energy.
Systems pharmacology is based on the theory of systems biology, integrating multi-disciplinary technologies such as multi-directional pharmacology, bioinformatics, and computer science to construct multi-level “disease-gene-target-drug” networks. This study analyzed and elucidated the links between CRP and LI using a modern systems pharmacology approach and in vitro experiments (Figure 1).
FIGURE 1. A comprehensive strategy diagram for the study of the mechanism of CRP action on LI.
Materials and Methods
Instruments and Herb Materials
While collecting and analyzing compounds, we used several instruments, including an ultra-pure liquid chromatograph (UPLC; Nexera X2 LC-30A, Shimadzu Corp. Kyoto, Japan)-hybrid triple-quadruple time-of-flight mass spectrometer (Triple TOF™5600+, AB Sciex, Toronto, Ontario, Canada) with an electrospray ionization source (ESI) and KQ-500 E ultrasonic cleaner (power 500 W, frequency 40 kHz Kunshan Ultrasonic Instrument Co., Ltd. Kunshan, China). The UPLC Titank C18 column (2.1 × 100 mm × 1.8 µm) was purchased from FLM Scientific Instruments Co., Ltd. (Guangzhou, China). Methanol, acetonitrile, and formic acid were purchased from ACS (Washington D.C. United States). Experimental water was ultrapure, obtained using a Milli-QB system (Bedford, MA, United States).
We obtained CRP medicinal materials (batch no. A162002–2) from the Jiangxi Xingan CRP planting base. They were formally identified by Professor Gong Qianfeng of Jiangxi University of TCM as dried and mature pericarps of Rutaceae plant oranges. The remaining samples were stored in the laboratory of the Jiangxi University of TCM. We used 22 standards (purity ≥98%) as references and in subsequent validation studies. Hesperidin (batch no. CHB18023), hesperetin (batch no. CHB180524), nobiletin (batch no. CHB180529), 3′-demethylnobiletin (batch no. CHB180131), 5-o-demethylnobiletin (batch no. CHB180321), narirutin (batch no. CHB:1809917), naringenin (batch no. CHB180914), sinensetin (batch no. CHB180126), isosinensetin (batch no. CHB180127), tangeretin (batch no. CHB190125), 5-hydroxy-4′,6,7,8-tetramethoxyflavon (batch no. CHB181119), 3,3′,4′,5,6,7,8-heptamethoxyflavone (batch no. CHB180120), 5,6,7,4′-tetramethoxyflavone (batch no. CHB190106), 5,7,4′-trimethoxyflavone (batch no. CHB181127), diosimin (batch no. CHB180125), diosmetin (batch no. CHB180312), poncirin (batch no. CHB180625), limonin (batch no. CHB180124), nomilin (batch no. CHB180315), naringenin chalcone (batch no. CHB180913), auraptene (batch no. CHB180529), and naringenin-7-O-β-D-glucoside (batch no. CHB181105) were all purchased from Chengdu Chroma-Biotechnology Co., Ltd. (Sichuan, China).
Screening of Active Compounds
Compounds Collected From UPLC-Q-TOF-MS/MS Analysis
Identification technology of compounds and metabolites in CRP has developed rapidly, and mass spectrometry identification is a common technical method (Zheng et al., 2013; Zeng et al., 2017; Li et al., 2019a). We added an appropriate amount of methanol to prepare standard solutions of 50–100 μg/ml and accurately weighed CRP powder (1.0 g), placed it in a conical flask containing 80 ml methanol, mixed it well, treated it with an ultrasound machine for 30 min, passed the solution through filter paper, then increased the volume to 100 ml, and finally filtered the solution through a 0.22 μm microporous membrane to obtain the CRP extraction solution. Ultra-Performance Liquid Chromatography-Quadrupole-Time-of-Flight Tandem Mass Spectrometry (UPLC-Q-TOF-MS/MS) conditions were as follows: chromatographic separation was performed at a flow rate of 0.25 ml/min at 40°C, and elution was performed with a linear gradient program according to various periods. The mobile phase system included solvent A (100% acetonitrile, v/v) and solvent B (0.01% formic acid in water, v/v). The process was solvent A (5–12%) for 5 min, (12–21%) for 10 min, (21–35%) for 10 min, (35–45%) for 11 min, (45–48%) for 4 min, (48–60%) for 1 min, (60–85%) for 9 min, (85–95%) for 1 min, (95%) for 2 min, (95–5%) for 1 min, and isocratic elution at 5% for 2.9 min. The ion source gases 1 and 2 were both set to 50 psi, curtain gas was set to 40 psi, ion spray voltage floating was set to 5500 V in the positive mode and 4500 V in the negative mode, ion source temperature was 500°C, collision energy was 40 V, collision energy spread was 15 V, declustering potential was 100 V, and nitrogen was used as a nebulizer and auxiliary gas. Samples were analyzed in both positive and negative ionization modes with a scanning mass-to-charge (m/z) range from 50 to 1,000. Data were collected in information-dependent acquisition mode and were analyzed using the Analyst TF 1.6 data acquisition workstation and PeakView®1.2 software. We compared the fragment ion information obtained for various compounds with that of the 22 standards and identified several CRP compounds. Also, we studied and analyzed the literature documenting fragment ion information of CRP components in detail carefully, matched it with our information of multiple compounds, and identified additional compounds.
Compounds Collected From Databases
To collect drug ingredients more thoroughly, we used The Encyclopedia of Traditional Chinese Medicine (http://www.tcmip.cn/ETCM), the Traditional Chinese Medicine Systems Pharmacology database and Analysis Platform (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php), and the TCM database Taiwan (http://tcm.cmu.edu.tw). All three online databases were searched for compounds, using the Chinese name “ChenPi” and the Latin name Citrus reticulata Blanco.
Compound Identification Analysis
We integrated the collected compounds and established a CRP chemical composition. First, the concept of oral bioavailability (OB) was introduced. OB refers to the percentage of unmodified drug entering the circulatory system after oral administration, representing drug utilization efficiency. Higher OB correlates with a higher possibility of clinical usefulness (Lv W. et al., 2020). Drug-likeness (DL) refers to the potential of a compound to become a drug (Huang et al., 2017). Molecules with OB ≥ 30% or DL ≥ 0.18 are considered to have better pharmacological effects (Zhang M. et al., 2020). We referred to this principle to select active ingredients for the subsequent step of the analysis. Considering the special biological activities of some compounds, we included them in the final candidate component database.
Identification of the Related Targets and Gene Symbols of CRP Compounds
We used information in the TCMSP (http://lsp.nwu.edu.cn/tcmsp.php) and Swiss Target prediction (STP, https://www.swisstargetprediction.ch/) to extract protein targets for each CRP component (Mu et al., 2020). After removing redundant information, the targets that could interact with each CRP component were retained.
Acquisition of Liver Injury Gene Targets
The gene targets for LI was obtained from the GeneCards database (https://www.genecards.org/, version 4.9.0) and the Online Mendelian Inheritance in Man database (OMIM, http://www.omim.org/, updated on February 28, 2019) (Mu et al., 2020; Xu et al., 2020). We performed keyword searches using “liver injury” and carefully screened the results. Finally, we obtained genes associated with LI-related diseases.
Drug-Compounds-Genes-Disease Network Construction
To obtain overlapping targets, we first cross-contrasted CRP-related targets with LI-related targets. We then built a complex information network (D-C-G-D) based on the interactions among drugs (CRP), compounds, gene symbols, and disease (LI) using Cytoscape software (version 3.8.0) and drew a schematic diagram.
Protein-Protein Interaction Network Construction
PPI data were obtained from the STRING database (https://string-db.org/, Version 11.0, updated on January 19, 2019). The interaction relationships provided by STRING are based on the confidence score and can be used to filter and assess functional genomics data (Tao et al., 2020). The organism species was set to Homo sapiens (Human), and the merged gene symbols were analyzed. We also used the BioGPS database (https://biogps.org) to identify proteins’ higher expression in some major organs (Gu et al., 2020). The heat map was generated using Heml software (version 1.0).
Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using R software (version 4.0.2). We used the R package “BiocGenerics,” a classic and practical tool specially designed for biologists to conduct an in-depth analysis of these targets (Yu et al., 2012). All enriched entries were determined using the Bioconductor database (http://bioconductor.org/).
Computational Validation of Ingredients-Targets Interactions
Molecular docking is a drug design method that considers receptors’ characteristics and the mode of interaction between receptors and drug molecules (Yang et al., 2019). We docked a variety of active components and key targets (BAX, BCL2, and caspase-3 [CASP3]) for a total of six component-target interactions. We obtained data from the RCSB Protein Data Bank (PDB, www.rcsb.org); the X-ray crystal structures of the key targets BAX, BCL2, and CASP3 were obtained and confirmed concerning the relevant literature; their PDB IDs were 1F16, 1G5M, and 5IBC, respectively. After processing targets and compounds using PyMOL software (version 1.3), docking work was performed using AUTODOCK VINA software (Version 1.1.2, Scripps Research, San Diego, CA, United States). The required input files for the AutoDock program were prepared using AutoDock tools.
Experimental Validation
Cell Culture
The LX-2 cell line was purchased from Beina Chuanglian Biotechnology Research Institute (HeNan, China). Cells were cultured in high-glucose Dulbecco’s minimum essential medium (Solarbio, Beijing, China) containing 10% fetal bovine serum (Tianhang Biotechnology, Zhejiang, China) in a 37°C, 5% CO2 incubator. Dimethyl sulfoxide was used to dissolve treatment compounds.
CCK-8 Assay for Cell Viability
Logarithmic growth phase cells were plated in 96-well plates (5 × 103 cells per well) and cultured overnight in a 37°C, 5% CO2 incubator. The following day, various hesperidin and naringenin (0, 25, 50, 75, 100, and 125 μM) were added to fresh medium. After 12 h, the model and administration groups were incubated with 5 ng/ml transforming growth factor (TGF)-β1 for 24 h, and then 10 μL CCK-8 assay solution (Solarbio) was added to each well. A microplate reader (Nanjing Detie Experimental Equipment Co., Ltd. Nanjing, China) was used to measure the optical density (OD) value around 450 nm. Cell survival was calculated as:
Real-Time Quantitative Polymerase Chain Reaction (qRT-PCR)
Total RNA was extracted using TRIzol® Reagent (Thermo Scientific, Waltham, MA, USA) and was reverse-transcribed with oligo-DT using HiScript™ Reverse transcriptase (Thermo Scientific) according to manufacturer instructions. The primers used were synthesized by Tsingke (Beijing, China). The sequences (forward and reverse, respectively), were as follows:
5′-CATGGGCTGGACATTGGACT-3′ and 5′-AAAGTAGGAGAGGAGGCCGT-3′ for BAX; 5′-TGAGTGCTCGCAGCTCATAC-3′ and 5′-TTCCCTGAGGTTTGCTGCAT-3′ for CASP3; 5′-CTTTGAGTTCGGTGGGGTCA-3′ and 5′-GAAATCAAACAGAGGCCGCA-3′ for BCL2; and 5′-GAAAGCCTGCCGGTGACTAA-3′ and 5′-TTCCCGTTCTCAGCCTTGAC-3′ for the internal control glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The GenBank accession numbers of BAX, CASP3, BCL2, and GAPDH used were NM_001291428.2, NM_001354777.2, NM_000633.3, and NM_001256799.3, respectively. Reverse transcription reaction conditions were: 30°C 10 min, 42°C 60 min, 99°C 5 min, and 4°C 5 min. The synthesized cDNA was stored at −20°C for later use. Gene-expression data were normalized to that of the endogenous control GAPDH. The 2-△△ Ct method was used as the basis for relative gene expression.
Western Blot
LX-2 cells were collected, lyzed in radioimmunoprecipitation assay buffer, and then centrifuged at 13,000 rpm for 10 min at 4°C. Bicinchoninic acid protein analysis kits (Solarbio) were used to calculate the samples’ protein concentrations based on standard curves. Protein samples (30 μg/sample) were added to the loading buffer, heated at 97°C for 6 min, and centrifuged at room temperature. After loading the samples, they were separated on 4% sodium dodecyl sulfate-polyacrylamide gels with voltage 150 V. The proteins were transferred to polyvinylidene membranes (Millipore, MA, United States) using a electroporation device. Next, the membranes were blocked with 5% skimmed milk powder in phosphate-buffered saline for 30 min. After incubation with a single antibody, a ChemiScope Mini 3,300 chemiluminescence imaging system (Shanghai, China) was used to image the blots. The primary antibodies were from Dingguo Changsheng Biotechnology Co., Ltd. (Beijing, China). Band intensity was analyzed using Quantity One software (version 4.6.6).
Statistical Analyses
The significance of results was determined based on one-way analysis of variance using Prism 8.0.1 software (Graph Pad Inc. San Diego, CA, United States). All experiments were performed in triplicate, and data are presented as mean ± SD. Differences were considered significant at p < 0.05.
Results
Screening of Active Compounds Collected From UPLC-Q-TOF-MS/MS Analysis
We showed a typical total ion chromatogram of alcohol-soluble components extracted from CRP (Figures 2A,B) and the mass spectrometric cleavage law of hesperidin and nobiletin (Figures 2C,D). By sampling the primary mass spectrum, we determined each component’s relative molecular weight and then obtained fragmentation information based on the secondary mass spectrum. Based on the retention time and mass spectrometric information of the standard material and the compounds reported in the literature, we identified a total of 38 compounds. Twenty-two of these (8, 9, 11, 12, 14, 16, 17, 18, 19, 21, 22, 24, 25, 26, 27, 28, 29, 31, 32, 34, 36, and 37) were identified through careful comparison of information between the samples and standard pure substances. The details are shown in Supplementary Table S1. Through the comparison of the multi-level fragment ion information recorded in the relevant literature and the analysis of the fragmentation law, the other 16 components were identified. Table 1 presents the details of the 38 compounds.
FIGURE 2. (A, B) Total ion chromatograms of CRP obtained by UPLC/Q-TOF-MS/MS in positive and negative ion mode. (C, D) Fragmentation Law of hesperidin and nobiletin by Mass Spectrometry.
TABLE 1. CRP compounds identified by UPLC/Q-TOF-MS/MS.
Compounds Collected From Databases
We searched online databases commonly used for TCM ingredients, The Encyclopedia of Traditional Chinese Medicine (http://www.tcmip.cn/ETCM), Traditional Chinese Medicine Systems Pharmacology database and Analysis Platform (http://lsp.nwu.edu.cn/tcmsp.php), and the TCM database Taiwan (http://tcm.cmu.edu.tw) to produce a more comprehensive compound list (Yu et al., 2020; Zhu and Hou, 2020). We collected a total of 65 compounds from these databases as candidates for CRP studies (Supplementary Table S2).
Compound Identification Analysis
We collated the results of previous work, combined the components obtained using UPLC/Q-TOF-MS/MS analysis and all CRP components identified in online databases. After removing repeated items, we established a complete CRP internal compound composition library. To screen out the active ingredients of CRP, the classical absorption, distribution, metabolism, and excretion (ADME) parameters of OB and DL were used to screen them. We found five ingredients that met these conditions; however, this result was far below our expectations. CRP contains a large number of flavonoids, including hesperidin and naringin (Dong et al., 2010). Flavonoid glycosides are more easily absorbed after hydrolysis into flavonoid aglycones (Hostetler et al., 2012), and hesperidin is one of the typical constituents (Nectoux et al., 2019). The OB of hesperidin was 13.33%, which was <30%; however, it is an active compound that is significantly effective against LI (Adefegha et al., 2017; Elhelaly et al., 2019). Narirutin has an OB of 8.15% and a DL of 0.75. However, it alleviates LI by regulating protein phosphorylation in the mitogen-activated protein kinase (MAPK) pathway, which may be related to lipid metabolism (Jianping et al., 2019; Gao et al., 2020). Tangeretin has an OB of 21.38% and a DL of 0.43, but it can ameliorate oxygen-glucose deprivation-induced LI through the JNK signaling pathway (Wu et al., 2019). Rosmarinic acid has an OB of 1.38% and a DL of 0.35; it can activate the AMP-activated protein kinase (AMPK) pathway and inhibit inflammation through nuclear factor changes erythroid 2-related factor and nuclear factor-κB signaling and ameliorates obesity and LI (Gyhye et al., 2019). Another mechanism by which compounds counteract LI is antioxidation, and the candidate p-Cymen-8-ol has both antioxidant effects and antibacterial activity (Pires et al., 2006; Santana et al., 2009). Although some components have relatively low kinetic values, they have excellent biological activity; therefore, we still consider them active CRP components. For these reasons, we propose that as long as the candidate components in CRP have good biological activity and intersect with LI’s target, they can be considered active CRP compounds for this indication. We selected 25 components as the final active compounds in the candidate compound library (Table 2).
TABLE 2. |Final selected compounds per the details of the active compounds in CRP.
Identification of Related Targets and Gene Symbols of Compounds in CRP
After we collected the relevant proteins from the TCMSP and STP databases, they were converted in the UniProt database, and redundant terms were removed to obtain 25 components in CRP and the 126 known associated target symbols (Supplementary Table S3).
Acquisition of Known Therapeutic Gene Targets for Liver Injury
A subset of the LI-related targets was obtained from the GeneCards database. We obtained 7372 known LI target symbols. Also, 92 known therapeutic target data for LI were obtained from OMIM. After eliminating redundancies, a total of 7442 known therapeutic targets for LI were collected (Supplementary Table S4).
Drug-Compounds-Genes-Disease Network Construction
We contrasted the obtained drug targets with disease-related genes to obtain their common cross genes. We constructed Venn diagrams (Figure 3A) with 117 overlaps of 7442 disease gene symbols and 126 drug-gene symbols. To intuitively demonstrate how CRP counteracts LI, we constructed a D-C-G-D network using Cytoscape software (Figure 3B). D-C-G-D details are in Supplementary Table S5.
FIGURE 3. (A) Venn diagram of related targets of CRP and LI. (B) D-C-G-D network. Green and red nodes indicate CRP and LI, respectively. Twenty-five orange nodes represent active ingredients in CRP; 117 blue nodes represent overlapping gene symbols between diseases and drugs, with edges indicating that nodes can interact; red edges indicate the action of LI with genes, green edges indicate the interaction of CRP with active ingredients, and gray edges indicate the interaction of active ingredients with genes. (C) The PPI network was obtained from the STRING database platform. (D) The PPI network is arranged according to degree value. (E) The expression of the core target in vital organs. Red and blue indicate higher and lower expression, respectively. (F) Bar graph of GO function enrichment of overlapping targets. (G) Bar graph of KEGG enrichment of overlapping targets. The Y-axis represents GO terms or KEGG pathways. The X-axis indicates the number of genes enriched in this pathway. The redder the color, the smaller the p. adjust value; it also indicates the reliability and importance. The bluer the color, the greater the p. adjust value.
Protein-Protein Interaction Network Construction
An action relationship map between genes was obtained after 117 overlapping genes were entered into the STRING online database (Figure 3C). These genes may be critical targets for LI treatment using CRP. The gene names are listed in Supplementary Table S6. After removing two free genes, we calculated the topological indices DC, BC, and CC for 115 genes. Degree (DC) indicates how much a gene is linked to other genes, suggesting that the gene may be necessary for the process of CRP counteracting LI. Detailed information about genes and topological indices is provided in Supplementary Table S7. We calculated the mean value of the 115 overlapping gene degrees to be 29.42. All genes with degree values >29.42 were selected, resulting in 25 core targets for interaction networks with other genes. After processing these associated targets using Cytoscape software, the inner loop of Figure 3D shows 25 genes with higher DC values. The deeper red color and large nodes indicate stronger correlations, and these may be critical genes for the CRP treatment of LI. To investigate LI and other major organs’ connection, we considered these core targets’ expression in various organs. According to the presented results of the heat map (Figure 3E), the 25 core targets’ expression was generally similar in the main organs. For example, AKT1 and MAPK3 were highly expressed in the liver, heart, lung, kidney, small intestine, and colon, while expression levels of PPARG and IL4 were low in these organs. This finding suggests that, in addition to the liver, LI is also closely related to other organs, consistent with the fact that TCM considers the human body as a whole. The specific expression of these core targets in various organs is shown in Supplementary Table S8.
Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment
To better summarize their specific functions, we performed GO enrichment analysis of 117 consensus genes and classified and elaborated them according to three modules: biological process (BP), cellular component (CC), and molecular function (MF), including 1525 GO terms for BP terms, 55 GO terms for CC terms, and 139 GO terms for MF terms, for a total of 1719 enriched GO terms. Based on the p. adjust value (p < 0.01), we selected the first ten terms from small to large for the enrichment results of BP, CC, and MF (Figure 3F). The BP analysis results revealed essential ways in which genes participate in bodily processes. These processes involve responses to nutrient levels (GO: 0031667), aging (GO: 0007568), responses to lipopolysaccharide (GO: 0032496), responses to a molecule of bacterial origin (GO: 0002237), and responses to metal ion (GO: 0010038). CC analysis then identified the sites where these intersection genes function intracellularly and extracellularly, primarily involving membrane raft (GO: 0045121), membrane microdomain (GO: 0098857), membrane region (GO: 0098589), caveola (GO: 0005901), and plasma membrane raft (GO: 0044853). Similarly, from MF analysis, we understand the forms of protein and target binding, primarily involving G protein-coupled amine receptor activity (GO: 0008227), adrenergic receptor activity (GO: 0004935), catecholamine binding (GO: 1901338), heme-binding (GO: 0020037), and tetrapyrrole binding (GO: 0046906). Overall, GO enrichment analysis results suggest that most of the effects of CRP on the liver occur in cell proliferation, regulation of enzyme activity, and energy metabolism. The detailed results are in Supplementary Table S9.
To further identify CRP’s potential pathways against LI, we also performed KEGG pathway enrichment analysis on 177 consensus genes, and the screen identified 49 pathways. For multiple pathways obtained by enrichment according to p. adjust (p < 0.01), we selected the top 20 paths for presentation (Figure 3G). The details of all KEGG pathway enrichment analyses are provided in Supplementary Table S10. Among these pathways, the interleukin (IL)-17 signaling pathway (hsa04657), endocrine resistance (hsa01522), the T cell receptor signaling pathway (hsa04660), Th1 and Th2 cell differentiation (hsa04658), platinum drug resistance (hsa01524), and others are critical pathways by which CRP acts on LI. Some of these are inflammatory pathways, some are apoptotic differentiation pathways, and others are receptor resistance pathways. Activation of inflammatory pathways such as IL-17 signaling (hsa04657) and tumor necrosis factor signaling (hsa04668) can significantly affect hepatocyte function (Yin and Feng, 2017; Doğanyiğit et al., 2020). Normal hepatocytes accelerate apoptosis after liver damage (Zhang T. et al., 2020). When CRP is administered to humans, the apoptosis pathway (hsa04210) may be affected by proteins such as BAX, BCL2, and CASP3 to maintain normal function. Liver function is also affected by the endocrine system, as endocrine resistance causes toxin-associated steatohepatitis (Clair, et al,2018).
Computational Validation of Ingredients-Targets Interactions
According to the GO and KEGG results, we found that CRP anti-LI is a biological process with multiple targets and pathways. Hesperidin is highly abundant in CRP and is a major active component (Xia et al., 2006). Nobiletin has suitable ADME parameters, suggesting that it may exhibit good drug performance during LI treatment. Based on PPI analysis, we know that some genes linked to naringenin have high metrics. We believe that components such as hesperidin, nobiletin, and naringenin simultaneously affect LI through common targets. Therefore, these compounds were subjected to molecular docking with apoptosis regulators BAX, CASP3, and BCL2 to test our hypothesis. The lower binding energy between molecules represents more potent force binding energy between them, and <0 kcal/mol favors the binding reaction (Shi et al., 2020).
According to the docking results, the binding energies of hesperidin to BAX and CASP3 were −8.2 and −8.6 kcal/mol, respectively; the binding energies of nobiletin to BAX and BCL2 were −5.9 and −7.7 kcal/mol, respectively; and the binding energies of naringenin to BCL2 and CASP3 were −7.7 and −7.2 kcal/mol, respectively. These results suggest that these compounds and their corresponding proteins have good binding. Figure 4 shows the combinations. Multiple linkages were formed between hesperidin and three amino acid residues in BAX, ARG-37, LEU-45, and ASP-48, including a C-O double bond and strong binding of the amino group (Figure 4A). Hesperidin interacts with two amino acid residues ASN-208 and ARG-207, in CASP3 (Figure 4B). These results suggest that a compound simultaneously forms multiple bonds with different amino acids in a protein, allowing the drug to affect the protein. The molecular docking results for nobiletin with BAX and BCL2 are shown in Figures 4C,D, respectively. The C-O double bond on nobiletin binds to the -OH group on SER-55 and the other O atom is connected to the H atom on ASP-519 (Figure 4C). The two O atoms on nobiletin and the amino acid residue TYR-9 on BCL2 and the H atom on ASP-196 form two hydrogen bonds (Figure 4D). Figures 4E,F show the docking results of naringenin with BCL2 and CASP3, respectively. The naringenin molecule binds to three residues (ASN-11, TYR-9, and HIS-186) on BCL2 to form multiple hydrogen bonds (Figure 4E). Naringenin also forms strong hydrophobic binding with residues ARG-207, SER-205, and SER-251 in CASP3 (Figure 4F). Our analysis showed that hesperidin, nobiletin, and naringenin effectively acts on protein targets closely related to apoptotic effects.
FIGURE 4. Selected compounds’ interactions with the targets. (A) Hesperidin with BAX, (B) hesperidin with CASP3, (C) nobiletin with BAX, (D) nobiletin with BCL2, (E) naringenin with BCL2, (F) naringenin with CASP3. The molecule is represented in a ball-stick model with atoms C, O, and N in green, red, and blue, respectively. Dashed lines indicate hydrogen bonds, and the numbers above represent distances in angstroms (Å).
Experimental Validation
Cell Viability
CCK-8 assays were performed to determine the effects of different doses of hesperidin and naringenin on LX-2 cell viability (Figures 5A,D). When hesperidin or naringenin was added at concentrations >75 μmol/L, cell proliferation was not significant, and viability remained high (100–130%). Based on this, three concentrations (75, 100, and 125 μmol/L) were selected for subsequent experiments.
FIGURE 5. Effect of hesperidin or naringenin on LX-2 cells. The effects of hesperidin (A) or naringenin (D) on LX-2 cell viability using the CCK-8 assay. mRNA expression of BAX (B), CASP3 (C, E), and BCL2 (F) were determined by qRT-PCR. Protein levels of CASP3 (G), BCL2 (H), AKT1 (I), pAKT1 (J) were determined by western blot, with results of each group (K). #p< 0.05, ##p< 0.01 versus blank control group; *p< 0.05, **p< 0.01 vs. TGF-β1-treated group.
Target Validation
TGF-β1-induced activation of LX-2 cells predisposes to liver fibrosis, resulting from chronic LI developing into cirrhosis that ultimately leads to liver failure (Bestion et al., 2020; Xiang et al., 2020). We chose hesperidin and naringenin to evaluate our prediction results further to validate changes in several relevant proteins in LX-2 cells with or without TGF-β1 (5 ng/ml) treatment. We performed qRT-PCR to determine the effect of hesperidin on BAX and CASP3 mRNA levels and that of naringenin on CASP3 and BCL2 mRNA levels. Western blot was used to measure naringenin’s effect on CASP3, BCL2, ATK1, and pAKT1 protein levels. We found that mRNA levels of BAX and CASP3 were significantly lower after TGF-β1 treatment (p < 0.01, Figure 5B,C,E), BCL2 mRNA was significantly elevated (p < 0.01, Figure 5F), and mRNA levels of CASP3 and BCL2 were consistent with their protein expression results (Figures 5G,H). Compared with the group that only received TGF-β1 stimulation, The levels of multiple index proteins were adjusted to different extents after treatment with different concentrations of hesperidin or naringenin, which made the results obtained after compound treatment more similar to the normal group. Western blotting revealed that ATK1 was only weakly affected by TGF-β1 and naringenin; protein expression changes in each group were relatively stable (Figure 5I). Protein levels of pAKT1 increased rapidly after TGF-β1 treatment but decreased significantly after naringenin treatment (Figure 5J). This suggests that naringenin may have a more pronounced inhibitory effect on AKT1 phosphorylation. These data collectively suggest that CRP reduces LI by mediating cell apoptosis and protein phosphorylation and regulating BAX expression, CASP3, BCL2, and ATK1, all of which are consistent with our molecular docking results. The results support our predicted results based on systems pharmacology.
Discussion
The liver is both a metabolic organ and an important mediator of immune function. LI may cause cancer, and the use of anti-cancer drugs may lead to more severe liver damage (Björnsson et al., 2020). TCM originates from nature and has few side effects, but its application and development have some limitations due to the complexity of the compositions and mechanisms of action. CRP is a pure, natural TCM that can also be used as a food, and a variety of compounds in CRP have been shown to have anti-proliferative effects on hepatoma cells (Chu et al., 2017).
In this study, we identified the active components and related targets of CRP. We linked them to LI-related targets to construct a D-C-G-D network that depicts links among 117 consensus genes and diseases. The PPI analysis showed complex interactions between 117 overlapping genes, and topological value analysis further indicated 25 core targets. We also examined these targets’ expression relationships in the liver, heart, lung, kidney, small intestine, and colon. The results showed that their expression levels were similar in these organs that belong to the TCM category of “five viscera and six entrails.” The theory of TCM posits that the viscera is an organic whole, the center of human life activities, and the organs work together to maintain the living body.
The GO analysis results revealed biological connections between drugs and diseases, and the BP enrichment clarified the response to the stimulation of external substances, such as response to acid chemical (GO:0001101). Hepatocytes are more sensitive to exogenous substance stimulation, and liver lesions develop when individuals regularly work in environments containing chemicals (Zhao et al., 2020). The CC analysis revealed the cellular environment in which the intersection protein or its product occurred, including membrane raft (GO: 0045121). Alcohol alters plasma membranes in vitro; preventing membrane raft oxidation is one approach to inhibit LI caused by excessive alcohol consumption (Chen et al., 2018). MF results suggest how intersection proteins or their products bind or function in vivo after the compound performs its action. Among them, catecholamine binding (GO: 1901338) is related to the expression of genes, and catecholamine receptors’ regulation prevents acute stress-induced LI (Zhu et al., 2014). We also found that a large number of GO terms were related to mitochondria, including regulation of mitochondrial membrane potential (hsa04210), apoptotic mitochondrial changes (GO: 0008637), and mitochondrial outer membrane permeabilization (GO: 0097345). This finding suggests that one of the critical mechanisms by which CRP ameliorates LI is to alter cellular energy metabolism (Xu et al., 2017; Dilberger et al., 2019).
The KEGG pathway enrichment results revealed that inflammatory and apoptotic pathways might be the primary CRP action pathways. CRP mediates the IL-17 signaling pathway by regulating targets such as PTGS2, MAPK3, and CASP3. PTGS2, also known as cyclooxygenase 2, extensively participates in hepatic inflammation induced by xenobiotics (Zhang et al., 2019). Through the D-C-G-D network, we know that hesperidin, nobiletin, and tangeretin are all involved in regulating PTGS2. The inhibition of PTGS2 mRNA expression is one of the manifestations of LI’s acute inflammation resolution (Li et al., 2020). It also suggests that a protein can be affected by several compounds simultaneously. LX-2 cells are prone to proliferate faster, and this can lead to liver fibrosis (Bestion et al., 2020). Our in vitro results suggest that CRP can induce apoptosis in LX-2 cells by regulating targets such as BAX, BCL2, and CASP3, consistent with our finding that the apoptosis pathway (hsa04210) is affected by CRP. The liver has a unique regenerative capacity (Nejak-Bowen and Monga, 2011). Some cancer pathways are associated with uncontrolled cell growth, and their dysregulation can lead to liver cancer. The phosphoinositide 3-kinase (PI3K)-Akt signaling pathway (hsa04151) can be activated by a variety of cell stimulants and toxins and is potentially correlated with liver cancer (Shi et al., 2020). Our prediction results suggest that AKT1 has a high degree in this pathway, and the phosphorylated form is a biomarker for cancer and tumor biology (Balasuriya et al., 2020). Our in vitro validation experiments demonstrated that naringenin significantly inhibited AKT1 phosphorylation, which in turn mediated PI3K-Akt signaling pathways to counteract LI. This interesting finding will be the focus of our future study. Nobiletin in CRP inhibits hepatoma metastasis resulting from PI3K-Akt signaling and may become a new compound for liver cancer treatment (Shi et al., 2013). The western blot results indicated that naringenin significantly reduced pAKT1 levels, suggesting that it prevents LI exacerbation. Increased cytochrome P450 expression improves hepatocyte function (Lewis et al., 2020); xenobiotics’ metabolism by cytochrome P450 (HSA00980) pathway was analyzed. We found that it resulted from naringenin, tangeretin, narirutin, and other compounds by affecting GSTP1, CYP1A1, and CBR1 targets.
We propose the following mechanism by which CRP counteracts LI. Hesperidin and naringenin in CRP affect the apoptosis pathway by acting on CASP3, BAX, and BCL2, reducing the likelihood of liver fibrosis. Naringenin significantly inhibits AKT1 phosphorylation, in turn mediating PI3K-Akt signaling to counteract LI. Systems pharmacology is a useful strategy to predict therapeutic mechanisms and identify new research directions. We acknowledge that the predictions and actual situations may be biased; however, they complement the in vitro findings, consistent with the molecular docking results that test the system pharmacology-based screening strategy’s reliability from a different perspective.
Conclusion
Previous drug development models mainly followed the concept of “one drug, one gene, and one disease.” In this study, modern technologies such as UPLC-Q-TOF-MS/MS and molecular docking were used to clarify the specific modes of the synergistic effects of multiple CRP components against LI, and provided convincing evidence. System pharmacology of TCM is an emerging multi-field interdisciplinary, but its development still has certain limitations. In addition to being an important reference for CRP in the treatment of liver diseases, our research also has many shortcomings. The specific interactions between all the drug components, proteins and multiple signal pathways involved in it need to be further studied. This approach of Systems pharmacology is a powerful new way to elucidate the mystery of TCM.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
JW and QG conceived and designed the study. SY conceived and designed the in vitro experimental validation. JW, XY, and SY contributed equally to this work and shared first authorship. All authors have read and approved the final manuscript.
Funding
This study was supported by the National Traditional Chinese Medicine Standardization Project (ZYBZH-Y-JX-27).
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.
Acknowledgments
The authors are thankful to Jiangxi University of Traditional Chinese Medicine for assistance in conducting this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.618846/full#supplementary-material.
References
Adefegha, S. A., Rosa Leal, D. B., Olabiyi, A. A., Oboh, G., and Castilhos, L. G. (2017). Hesperidin attenuates inflammation and oxidative damage in pleural exudates and liver of rat model of pleurisy. Redox Rep. 22 (6), 563–571. doi:10.1080/13510002.2017.1344013
Balasuriya, N., Davey, N. E., Johnson, J. L., Liu, H., Biggar, K. K., Cantley, L. C., et al. (2020). Phosphorylation-dependent substrate selectivity of protein kinase B (AKT1). J. Biol. Chem. 295 (24), 8120–8134. doi:10.1074/jbc.ra119.012425
Bestion, E., Jilkova, Z., Mège, J., Novello, M., Kurma, K., Ahmad Pour, S. T., et al. (2020). GNS561 acts as a potent anti-fibrotic and pro-fibrolytic agent in liver fibrosis through TGF-β1 inhibition[J]. Ther. Adv. chronic Dis. 11, 2040622320942042. doi:10.1177/2040622320942042
Björnsson, H. K., Sverrisdottir, A., and Bjornsson, E. (2020). 1088 drug-induced liver injury due to chemotherapeutic agents and other cancer medications. A. Population Based Study[j]. Gastroenterol. 158 (6).
Chen, Y., He, X., Yuan, X., Hong, J., Bhat, O., Li, G., et al. (2018). NLRP3 inflammasome formation and activation in nonalcoholic steatohepatitis: therapeutic target for antimetabolic syndrome remedy FTZ. Oxid Med. Cel Longev 2018, 2901871. doi:10.1155/2018/2901871
Chu, C-C., Chen, S-Y., Chyau, C-C., and Duh, P-D. (2017). Antiproliferative effect of sweet orange peel and its bioactive compounds against human hepatoma cells, in vitro and in vivo[J]. J. Funct. Foods 33. doi:10.1016/j.jff.2017.03.051
Clair, H. B., Pinkston, C. M., Rai, S. N., Pavuk, M., Dutton, N. D., et al. (2018). Liver disease in a residential cohort with elevated polychlorinated biphenyl exposures[J]. Toxicol. Sci. : official J. Soc. Toxicol. 164 (1), 39–49. doi:10.1093/toxsci/kfy076
Cui, Y., Zhao, J., Zhou, J., Tan, G., Zhao, Q., Zhang, Y., et al. (2019). Development of a sensitive monoclonal antibody-based indirect competitive enzyme-linked immunosorbent assay for analysing nobiletin in citrus and herb samples. Food Chem. 293, 144–150. doi:10.1016/j.foodchem.2019.04.101
Dilberger, B., Baumanns, S., Schmitt, F., Schmiedl, T., Hardt, M., Wenzel, U., et al. (2019). Mitochondrial oxidative stress impairs energy metabolism and reduces stress resistance and longevity of[J]. Oxidative medicine and cellular longevity of C. elegans. Oxid Med. Cel Longev 2019, 6840540. doi:10.1155/2019/6840540
Dong, P., Qiu, P., Zhu, Y., Li, S., Ho, C.-T., McClements, D. J., et al. (2010). Simultaneous determination of four 5-hydroxy polymethoxyflavones by reversed-phase high performance liquid chromatography with electrochemical detection. J. Chromatogr. A 1217 (5), 642–647. doi:10.1016/j.chroma.2009.11.097
Doğanyiğit, Z., Okan, A., Kaymak, E., Pandır, D., and Silici, S. (2020).Investigation of protective effects of apilarnil against lipopolysaccharide induced liver injury in rats via TLR 4/HMGB-1/NF-κB pathway. Biomedicine & Pharmacotherapy, 125, 109967. doi:10.1016/j.biopha.2020.109967
Elhelaly, A. E., Albasher, G., Alfarraj, S., Almeer, R., Bahbah, E. I., Fouda, M. M. A., et al. (2019). Protective effects of hesperidin and diosmin against acrylamide-induced liver, kidney, and brain oxidative damage in rats. Environ. Sci. Pollut. Res. 26 (34), 35151–35162. doi:10.1007/s11356-019-06660-3
Feumba Dibanda, R., Panyoo Akdowa, E., Rani P., A., Metsatedem Tongwa, Q., and Mbofung F., C. M. (2020). Effect of microwave blanching on antioxidant activity, phenolic compounds and browning behaviour of some fruit peelings. Food Chem. 302, 125308. doi:10.1016/j.foodchem.2019.125308
Gao, Y., Tian, Y., Zhang, X., Zhang, X., Duan, Z., Ren, F., et al. (2020). Magnesium isoglycyrrhizinate ameliorates concanavalin A-induced liver injury via the p38 and JNK MAPK pathway. Immunopharmacology and immunotoxicology 42 (5), 445–455. doi:10.1080/08923973.2020.1808984
Gu, S., Xue, Y., Gao, Y., Shen, S., Zhang, Y., Chen, K., et al. (2020). Mechanisms of indigo naturalis on treating ulcerative colitis explored by GEO gene chips combined with network pharmacology and molecular docking[J]. Scientific Rep. 10 (1), 15204. doi:10.1038/s41598-020-71030-w
Gyhye, Y., Myungsuk, K., Ahmad, R., Yang-Ju, S., Hong, C. R., Kim, S. M., et al. (2019). Lemon balm extract and its major chemical compound, rosmarinic acid, alleviate damages of liver in an animal model of nonalcoholic steatohepatitis (NASH) (P06-093-19)[J]. Curr. Dev. Nutr. 3 (Suppl. 1).
Hostetler, G., Riedl, K., Cardenas, H., Diosa-Toro, M., Arango, D., Schwartz, S., et al. (2012). Flavone deglycosylation increases their anti-inflammatory activity and absorption. Mol. Nutr. Food Res. 56 (4), 558–569. doi:10.1002/mnfr.201100596
Huang, J., Cheung, F., Tan, H.-Y., Hong, M., Wang, N., Yang, J., et al. (2017). Identification of the active compounds and significant pathways of yinchenhao decoction based on network pharmacology. Mol. Med. Rep. 16 (4), 4583–4592. doi:10.3892/mmr.2017.7149
Jianping, J., Li, Y., Zheng, S., Lixia, W., Letian, S., and Thomas, E. (2019). Hepatoprotective and anti-inflammatory effects of total flavonoids of Qu Zhi Ke (peel of Citrus changshan-huyou) on non-alcoholic fatty liver disease in rats via modulation of NF-κB and MAPKs[J]. Phytomedicine : Int. J. phytotherapy phytopharmacology 64.
Johnson, A., Hjelmeland, A., Heymann, H., and Ebeler, S. (2019). GC-Recomposition-Olfactometry (GC-R) and multivariate study of three terpenoid compounds in the aroma profile of Angostura bitters[J]. Scientific Rep. 9 (1), 7633. doi:10.1038/s41598-019-44064-y
Kim, J. H., Cha, J. Y., Shin, T. S., and Chun, S. S. (2018). Volatile flavor components of blended tea with fermented tea and herbs. pnf 23 (3), 245–253. doi:10.3746/pnf.2018.23.3.245
Lewis, B., Paula, S., Paine, S. W., Mellor, H. R., Abernethy, L., Littlewood, P., et al. (2020). Hydrostatic pressure regulates CYP1A2 expression in human hepatocytes via a mechanosensitive aryl hydrocarbon receptor-dependent pathway[J]. Am. J. Physiol. Cel Physiol. 318 (5).
Li, S.-Z., Guan, X.-M., Gao, Z., Lan, H.-C., Yin, Q., Chu, C., et al. (2019a). A simple method to discriminate Guangchenpi and Chenpi by high-performance thin-layer chromatography and high-performance liquid chromatography based on analysis of dimethyl anthranilate. J. Chromatogr. B 1126-1127, 121736. doi:10.1016/j.jchromb.2019.121736
Li, S.-Z., Zeng, S.-L., Wu, Y., Zheng, G.-D., Chu, C., Yin, Q., et al. (2019b). Cultivar differentiation of Citri Reticulatae Pericarpium by a combination of hierarchical three-step filtering metabolomics analysis, DNA barcoding and electronic nose. Analytica Chim. Acta 1056, 62–69. doi:10.1016/j.aca.2019.01.004
Li, Y., Chen, Y., Huang, X., Huang, D., Gan, H., Yao, N., et al. (2020). Tanshinol A ameliorates triton‐1339W‐induced hyperlipidemia and liver injury in C57bl/6J mice by regulating mRNA expression of lipemic‐oxidative injury genes. Lipids 55 (2), 127–140. doi:10.1002/lipd.12217
Lin, Z.-H., Chan, Y.-F., Pan, M.-H., Tung, Y.-C., and Su, Z.-Y. (2019). Aged citrus peel (chenpi) prevents acetaminophen-induced hepatotoxicity by epigenetically regulating Nrf2 pathway. Am. J. Chin. Med. 47 (8), 1833–1851. doi:10.1142/s0192415x19500939
Liu, P. (2016). Inhibition of pathological angiogenesis of Chinese medicine against liver fibrosis. Chin. J. Integr. Med. 22 (8), 569–572. doi:10.1007/s11655-016-2468-3
Luo, Y., Zeng, W., Huang, K.-E., Li, D.-X., Chen, W., Yu, X.-Q., et al. (2019). Discrimination of Citrus reticulata Blanco and Citrus reticulata 'Chachi' as well as the Citrus reticulata 'Chachi' within different storage years using ultra high performance liquid chromatography quadrupole/time-of-flight mass spectrometry based metabolomics approach. J. Pharm. Biomed. Anal. 171, 218–231. doi:10.1016/j.jpba.2019.03.056
Lv, W., Lin, T., Ren, Z., Jiang, Y., Zhang, J., Bi, F., et al. (2020a). Rapid discrimination of Citrus reticulata 'Chachi' by headspace-gas chromatography-ion mobility spectrometry fingerprints combined with principal component analysis. Food Res. Int. 131, 108985. doi:10.1016/j.foodres.2020.108985
Lv, X., Xu, Z., Xu, G., Li, H., Wang, C., Chen, J., et al. (2020b). Investigation of the active components and mechanisms of Schisandra chinensis in the treatment of asthma based on a network pharmacology approach and experimental validation. Food Funct. 11 (4), 3032–3042. doi:10.1039/d0fo00087f
Mu, C., Sheng, Y., Wang, Q., Amin, A., Li, X., and Xie, Y. (2020). Potential compound from herbal food of rhizoma polygonati for treatment of COVID-19 analyzed by network pharmacology and molecular docking technology[J]. J. Funct. Foods 77, 104149. doi:10.1016/j.jff.2020.104149
Nectoux, A. M., Abe, C., Huang, S.-W., Ohno, N., Tabata, J., Miyata, Y., et al. (2019). Absorption and metabolic behavior of hesperidin (rutinosylated hesperetin) after single oral administration to sprague-dawley rats. J. Agric. Food Chem. 67 (35), 9812–9819. doi:10.1021/acs.jafc.9b03594
Nejak-Bowen, K. N., and Monga, S. P. S. (2011). Beta-catenin signaling, liver regeneration and hepatocellular cancer: sorting the good from the bad. Semin. Cancer Biol. 21 (1), 44–58. doi:10.1016/j.semcancer.2010.12.010
Nie, H., Deng, Y., Zheng, C., Pan, M., Xie, J., Zhang, Y., et al. (2020). A network pharmacology‐based approach to explore the effects of Chaihu Shugan powder on a non‐alcoholic fatty liver rat model through nuclear receptors. J. Cel Mol Med 24 (9), 5168–5184. doi:10.1111/jcmm.15166
Pires, A. M. L., Albuquerque, M. R. J. R., Nunes, E. P., Melo, V. M. M., Silveira, E. R., and Pessoa, O. D. L., (2006). Chemical composition and antibacterial activity of the essential oils of blainvillea rhomboidea (asteraceae)[J]. Nat. Product. Commun. 1 (5). doi:10.1177/1934578x0600100510
Prc, P. C. O. (2015). Pharmacopoeia of the people's Republic of China[M]. 1. Beijing: China Medical Science Press.
Qiao, Y., Xie, B., Zhang, Y., Zhang, Y., Fan, G., Yao, X., et al. (2008). Characterization of aroma active compounds in fruit juice and peel oil of jinchen sweet orange fruit (citrus sinensis (L.) osbeck) by GC-MS and GC-O. Molecules 13 (6), 1333–1344. doi:10.3390/molecules13061333
Run-Xia, G., Yu, Z., Xing-He, T., Feng, W., Wen, C., and Xiao-an, H. (2013). Change in essential oil composition of ponkan peels under different storage conditions[J]. Beijing: Food ence.
Santana, A. I., Vila, R., Espinosa, A., Olmedo, D., Gupta, M. P., and Cañiguera, S. (2009). Composition and biological activity of essential oils from Protium confusum[J]. Nat. Prod. Commun. 4 (10), 1401–1406. doi:10.1177/1934578x0900401017
Shi, M.-D., Liao, Y.-C., Shih, Y.-W., and Tsai, L.-Y. (2013). Nobiletin attenuates metastasis via both ERK and PI3K/Akt pathways in HGF-treated liver cancer HepG2 cells. Phytomedicine 20, 743–752. doi:10.1016/j.phymed.2013.02.004
Shi, P., Xie, Y., Xie, R., Lin, Z., Yao, H., and Wu, S. (2020). Acanthopanax senticosusAn integrated pharmacokinetic study of an extract preparation by combination of virtual screening, systems pharmacology, and multi-component pharmacokinetics in rats[J]. Front. Pharmacol. 11, 1295. doi:10.3389/fphar.2020.01295
Shi, Q., Xue, C., Yuan, X., and Yu, Z. (2020). Gene signatures and prognostic values of m1A-related regulatory genes in hepatocellular carcinoma[J]. Scientific Rep. 10 (1), 15083. doi:10.1038/s41598-020-72178-1
Tao, Q., Du, J., Li, X., Zeng, J., Tan, B., Xu, J., et al. (2020). Network pharmacology and molecular docking analysis on molecular targets and mechanisms of Huashi Baidu formula in the treatment of COVID-19. Drug Dev. Ind. Pharm. 46 (8), 1345–1353. doi:10.1080/03639045.2020.1788070
Tsutsumi, M., Wang, J.-S., and Takada, A. (1994). Microheterogeneity of serum glycoproteins in alcoholics: is desialo-transferrin the marker of chronic alcohol drinking or alcoholic liver injury?. Alcohol. Clin. Exp. Res. 18 (2), 392–397. doi:10.1111/j.1530-0277.1994.tb00031.x
Vatsalya, V., Song, M., Schwandt, M. L., Cave, M. C., Barve, S. S., George, D. T., et al. (2016). Effects of sex, drinking history, and omega-3 and omega-6 fatty acids dysregulation on the onset of liver injury in very heavy drinking alcohol-dependent patients. Alcohol. Clin. Exp. Res. 40 (10), 2085–2093. doi:10.1111/acer.13197
Wang, Y., Yi, L., Liang, Y., Li, H., Yuan, D., Gao, H., et al. (2008). Comparative analysis of essential oil components in Pericarpium Citri Reticulatae Viride and Pericarpium Citri Reticulatae by GC-MS combined with chemometric resolution method. J. Pharm. Biomed. Anal. 46 (1), 66–74. doi:10.1016/j.jpba.2007.08.030
Wu, C., Zhao, J., Chen, Y., Li, T., Zhu, R., Zhu, B., et al. (2019). Tangeretin protects human brain microvascular endothelial cells against oxygen‐glucose deprivation‐induced injury. J. Cel Biochem 120 (4), 4883–4891. doi:10.1002/jcb.27762
Xia, J., Kotani, A., Hakamata, H., and Kusu, F. (2006). Determination of hesperidin in Pericarpium Citri Reticulatae by semi-micro HPLC with electrochemical detection. J. Pharm. Biomed. Anal. 41 (4), 1401–1405. doi:10.1016/j.jpba.2006.02.030
Xiang, D., Zou, J., Zhu, X., Chen, X., Luo, J., Kong, L., et al. (2020). Physalin D attenuates hepatic stellate cell activation and liver fibrosis by blocking TGF-β/Smad and YAP signaling. Phytomedicine 78, 153294. doi:10.1016/j.phymed.2020.153294
Xiao, H., Yang, C. S., Li, S., Jin, H., Ho, C.-T., and Patel, T. (2009). Monodemethylated polymethoxyflavones from sweet orange (Citrus sinensis) peel Inhibit growth of human lung cancer cells by apoptosis. Mol. Nutr. Food Res. 53 (3), 398–406. doi:10.1002/mnfr.200800057
Xu, F., Liu, Y., Zhao, H., Yu, K., Song, M., Zhu, Y., et al. (2017). Aluminum chloride caused liver dysfunction and mitochondrial energy metabolism disorder in rat. J. Inorg. Biochem. 174, 55–62. doi:10.1016/j.jinorgbio.2017.04.016
Xu, H.-H., Li, S.-M., Xu, R., Fang, L., Xu, H., and Tong, P.-J. (2020). Predication of the underlying mechanism of Bushenhuoxue formula acting on knee osteoarthritis via network pharmacology-based analyses combined with experimental validation. J. ethnopharmacology 263, 113217. doi:10.1016/j.jep.2020.113217
Yang, Q., Han, B., Xue, J., Lv, Y., Li, S., Liu, Y., et al. (2020). Hexavalent chromium induces mitochondrial dynamics disorder in rat liver by inhibiting AMPK/PGC-1α signaling pathway. Environ. Pollut. 265, 114855. doi:10.1016/j.envpol.2020.114855
Yang, S., Zhang, J., Yan, Y., Yang, M., Li, C., Li, J., et al. (2019). Network pharmacology-based strategy to investigate the pharmacologic mechanisms of atractylodes macrocephala koidz. For the treatment of chronic gastritis. Front. Pharmacol. 10, 1629. doi:10.3389/fphar.2019.01629
Yi, L., Dong, N., Liu, S., Yi, Z., and Zhang, Y. (2015). Chemical features of Pericarpium Citri Reticulatae and Pericarpium Citri Reticulatae viride revealed by GC-MS metabolomics analysis. Food Chem. 186, 192–199. doi:10.1016/j.foodchem.2014.07.067
Yin, S., and Feng, D. (2017). IL-22 in liver inflammation, injury and repair[J]. Liver Research. 1 (3), 181-185. doi:10.1016/j.livres.2017.09.003
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for comparing biological themes among gene clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yu, J.-W., Yuan, H.-W., Bao, L.-D., and Si, L.-G. (2020). Interaction between piperine and genes associated with sciatica and its mechanism based on molecular docking technology and network pharmacology. Mol. Divers. 25 (1), 233-248. doi:10.1007/s11030-020-10055-9
Yu, X., Sun, S., Guo, Y., Liu, Y., Yang, D., Li, G., et al. (2018). Citri Reticulatae Pericarpium (Chenpi): botany, ethnopharmacology, phytochemistry, and pharmacology of a frequently used traditional Chinese medicine. J. ethnopharmacology 220, 265–282. doi:10.1016/j.jep.2018.03.031
Zeng, S., Li, S., Xiao, P., Cai, Y., Chu, C., Chen, B., et al. (2020). Citrus polymethoxyflavones attenuate metabolic syndrome by regulating gut microbiome and amino acid metabolism[J]. Sci. Adv. 6 (1), eaax6208. doi:10.1126/sciadv.aax6208
Zeng, S.-L., Duan, L., Chen, B.-Z., Li, P., and Liu, E.-H. (2017). Chemicalome and metabolome profiling of polymethoxylated flavonoids in Citri Reticulatae Pericarpium based on an integrated strategy combining background subtraction and modified mass defect filter in a Microsoft Excel Platform. J. Chromatogr. A 1508, 106–120. doi:10.1016/j.chroma.2017.06.015
Zhang, H. Y., Wang, H. L., Zhong, G. Y., and Zhu, J. X. (2018). Molecular mechanism and research progress on pharmacology of traditional Chinese medicine in liver injury. Pharm. Biol. 56 (1), 594–611. doi:10.1080/13880209.2018.1517185
Zhang, L-Y., Zhan, D-L., Chen, Y-Y., Wang, W-H., He, C-Y., Lin, Y., et al. (2019). Aflatoxin B1 enhances pyroptosis of hepatocytes and activation of Kupffer cells to promote liver inflammatory injury via dephosphorylation of cyclooxygenase-2: an in vitro, ex vivo and in vivo study[J]. Arch. Toxicol. 93 (11). doi:10.1007/s00204-019-02572-w
Zhang, M., Yuan, Y., Zhou, W., Qin, Y., Xu, K., Men, J., et al. (2020a). Network pharmacology analysis of Chaihu Lizhong Tang treating non-alcoholic fatty liver disease. Comput. Biol. Chem. 86, 107248. doi:10.1016/j.compbiolchem.2020.107248
Zhang, T., Li, H., Wang, K., Xu, B., Chen, Z-N., and Bian, H. (2020b). Deficiency of CD147 attenuated non-alcoholic steatohepatitis progression in an NLRP3-dependent manner[J]. Front. Cel. Dev. Biol. 8, 784. doi:10.3389/fcell.2020.00784
Zhao, X., Chen, T., Wang, D., Du, Y., Wang, Y., Zhu, W., et al. (2020). Polybrominated diphenyl ethers and decabromodiphenyl ethane in paired hair/serum and nail/serum from corresponding chemical manufacturing workers and their correlations to thyroid hormones, liver and kidney injury markers. Sci. Total Environ. 729, 139049. doi:10.1016/j.scitotenv.2020.139049
Zheng, G., Chao, Y., Luo, M., Xie, B., Zhang, D., Hu, P., et al. (2020a). Construction and chemical profile on "activity fingerprint" of Citri Reticulatae Pericarpium from different cultivars based on HPLC-UV, LC/MS-IT-TOF, and principal component analysisEvid. Based Complement. Alternat Med., 2020, 4736152. doi:10.1155/2020/4736152
Zheng, G., Liu, M., Chao, Y., Yang, Y., Zhang, D., Tao, Y., et al. (2020b). Identification of lipophilic components in Citri Reticulatae Pericarpium cultivars by supercritical CO fluid extraction with ultra-high-performance liquid chromatography-Q Exactive Orbitrap tandem mass spectrometry[J]. J. Separat. Sci. 43 (17), 3421–3440. doi:10.1002/jssc.202000490
Zheng, G.-D., Zhou, P., Yang, H., Li, Y.-s., Li, P., and Liu, E.-H. (2013). Rapid resolution liquid chromatography-electrospray ionisation tandem mass spectrometry method for identification of chemical constituents in Citri Reticulatae Pericarpium. Food Chem. 136 (2), 604–611. doi:10.1016/j.foodchem.2012.08.040
Zheng, Y. y., Zeng, X., Peng, W., Wu, Z., and Su, W. w. (2019). Characterisation and classification of Citri Reticulatae Pericarpium varieties based on UHPLC‐Q‐TOF‐MS/MS combined with multivariate statistical analyses. Phytochem. Anal. 30 (3), 278–291. doi:10.1002/pca.2812
Zhu, N., and Hou, J. (2020). Exploring the mechanism of action Xianlingubao Prescription in the treatment of osteoporosis by network pharmacology. Comput. Biol. Chem. 85, 107240. doi:10.1016/j.compbiolchem.2020.107240
Keywords: systems pharmacology, experiment verification, liver injury, mechanism of action, citri reticulatae pericarpium
Citation: Wu J, Ye X, Yang S, Yu H, Zhong L and Gong Q (2021) Systems Pharmacology Study of the Anti-Liver Injury Mechanism of Citri Reticulatae Pericarpium. Front. Pharmacol. 12:618846. doi: 10.3389/fphar.2021.618846
Received: 18 October 2020; Accepted: 03 March 2021;
Published: 12 April 2021.
Edited by:
Xuezhou Zhou, Beijing Jiaotong University, ChinaReviewed by:
Min Wu, China Academy of Chinese Medical Sciences, ChinaE-Hu Liu, China Pharmaceutical University, China
Copyright © 2021 Wu, Ye, Yang, Yu, Zhong and Gong. 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: Huan Yu, yuhuanhebei@163.com; Lingyun Zhong, ly1638163@163.com; Qianfeng Gong, gongqf2002@163.com
†These authors have contributed equally to this work