Abstract

Objective. To investigate the pharmacological mechanism of the iridoid glycosides from Fructus Gardeniae in Jiangxi province by network pharmacology. To provide a valuable research strategy for the rational use and in-depth research and development of Fructus Gardeniae from Jiangxi. Method. Previous research results of our group show that the contents of iridoid glycosides in Fructus Gardeniae from Jiangxi province have a significant difference compared with other regions (). Based on our previous experimental results, this study selected six characteristic high-content bioactive iridoid glycosides components of Fructus Gardeniae from Jiangxi province as candidate components. TCMSP database was used to obtain the process parameters of absorption, distribution, metabolism, and excretion (ADME) of candidate components. PubChem and SWISS online database were used to predict the related targets. Cytoscape software was used to the construct compound-target-disease (C-T-D) network of the Fructus Gardeniae iridoid glycosides ingredients. Furthermore, the GO biological process analysis and the pathway enrichment analysis were carried out using the CTD online analysis platform; then, an illustrated network that contains the main “chemicals-targets-pathway (C-T-P)” was constructed to analyze main biological pathways for obtaining the deep mechanism of Fructus Gardeniae in Jiangxi. Results. 6 iridoid glycosides, namely geniposide, gardenoside, geniposidic acid, genipin 1-gentiobioside, gardoside, and shanzhiside, from Fructus Gardeniae in Jiangxi province were obtained as candidate components through previous work and network pharmacology screening. 36 corresponding targets were acted, such as BCL2, MAPT, F2, BCL2L1, PRKCD, PRKCB, HIF1A, and PRKCA. These targets could joint in pathways, such as signaling by GPCR, neuroactive ligand-receptor interaction, inflammatory mediator regulation of TRP channels, and ion channel transport. Interestingly, these pathways were highly associated with liver diseases, neurological diseases, hypertension, neoplasms, hyperalgesia, and inflammation. Remarkably, we boldly speculate that the Fructus Gardeniae from Jiangxi province can play a pharmacological role in hepatic encephalopathy through regulating multiple signaling pathways in an integrated manner. Conclusion. The method based on system pharmacology could help to find the key targets of characteristic high-content chemical constituents of herb from different producing areas, the signaling pathway and disease network of TCM, and provide useful information and data support for giving a further study on traditional Chinese medicine resources in different regions of China.

1. Introduction

Gardenia jasminoides Ellis (Rubiaceae) is an evergreen shrub. The dried ripe fruits of Gardenia jasminoides have been recorded as Fructus Gardeniae (“zhi zi” in Chinese) in the Chinese Pharmacopoeia (2000–2015 edition) [1], and the fruits are used in traditional Chinese medicine because of their diuretic, antiphlogistic, anti-inflammatory, hemostatic, antipyretic, cholagogue, and laxative effects [24]. Fructus Gardeniae has been reported to protect the liver and gallbladder cells from infringement, treat cardiovascular and cerebrovascular diseases, and promote ligament cell proliferation and collagen synthesis effectively [5, 6]. Researchers found that the chemical composition of Fructus Gardeniae includes iridoid glycosides [5, 7], crocins [8], phenylpropanoids [9], glycoprotein [10], and polysaccharides. Among those, iridoid glycosides are the major active constituents in Fructus Gardeniae [1113].

It is well known that traditional Chinese medicine (TCM) are normally collected from multiple geographical regions [14]. Accumulating studies have demonstrated that the place of origin could significantly affect the quality of TCM since climate and environment influence biosynthesis and accumulation of secondary metabolites in an organism [15]. The geographical origins of Fructus Gardeniae are also multiple. The contents of chemical constituents in Fructus Gardeniae from different habitats are different, and the different chemical constituents can regulate the subjects through multiple targets and pathways, thus producing different pharmacological effects. In our previous studies, Fructus Gardeniae cultivated in 10 different provinces in China has different high-content bioactive components and therefore, their different pharmacological activities. Preexperimental results demonstrated that the contents of geniposide, gardenoside, geniposidic acid, genipin 1-gentiobioside, gardoside, and shanzhiside in Fructus Gardeniae from Jiangxi province were significantly higher than those in other regions () [16]; therefore, 6 iridoid glycosides were selected as the candidate components for the further pharmacological mechanism prediction. In this study, the key targets and signaling pathways of the characteristic high-content bioactive components group of Fructus Gardeniae from Jiangxi province were searched by network pharmacology, and the related diseases were predicted, and the chemicals-target-disease network was constructed. This study provides a basis for the in-depth study of the pharmacological mechanism of Fructus Gardeniae from Jiangxi province, China. The schematic illustration of this study is shown in Figure 1.

2. Materials and Methods

2.1. Candidate Components and ADME Parameters

Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform (http://lsp.nwu.edu.cn/index.php, TCMSP) is a database of systems pharmacology for drug discovery from herbal medicines [17]. A growing number of systems pharmacology-based studies have recently been conducted in efforts to decipher the active mechanisms of TCM. Chen et al. predicted the potential mechanism of Yinchensini decoction based on systemic pharmacology by using TCMSP database. We followed the methods of Chen et al. [18]. From the TCMSP database, geniposide, gardenoside, geniposidic acid, genipin 1-gentiobioside, gardoside, and shanzhiside in Fructus Gardeniae from Jiangxi province are retrieved. And the process parameters of absorption, distribution, metabolism, and excretion (ADME) of candidate components were obtained by the TCMSP database. Each candidate’s drug ability was analyzed according to its oral bioavailability (OB) and drug-likeness (DL) indices recommended by TCMSP. OB refers to the speed and degree of absorbing drugs into the circulatory system, which is a reliable indicator to evaluate the intrinsic quality of drugs objectively. The higher the OB of the compound is, the more likely the compound has the potential of being developed into drugs for clinical application [19]. DL refers to the comparative measure of the functional or physical properties of compounds with those of the majority of known drugs [20], which is the sum of the pharmacokinetic properties and safety. It can be used to analyze the results of the drug activity, predict in vivo pharmacokinetics, and optimize compounds. Bioactive molecules with high DL were often considered to exhibit relatively better pharmacologically.

2.2. Targets Screening

The targets prediction for the characteristic high-content bioactive components group of Fructus Gardeniae from Jiangxi province were performed using Swiss Target Prediction (http://www.swisstargetprediction.ch/) [21] by inputting the canonical SMILES into SMILES string (s) [22, 23], with the organism selected as “Homo sapiens”. The compound targets having no relationship with the compounds were deleted. Swiss is a free public resource to explore interactions between chemicals and gene targets. The potential targets are ranked by probability from high to low. And gene information including gene ID, name, and organism was identified using the UniProt database (http://www.uniprot.org/).

2.3. Chemicals-Target-Disease Network Construction

In this study, a “chemicals-target-disease (C-T-D)” network related to Fructus Gardeniae from Jiangxi province was constructed. In order to elucidate the relationship among candidate compounds, potential targets and related diseases, a “compound-target (C-T)” network, and a “disease-target (D-T)” network were generated for studying. (1) C-T network: compound-target interactions were visualized by the C-T network, in which all the active ingredients were connected to their corresponding targets. (2) T-D network: the target name was inputted to the CTD (http://ctdbase.org/) to predict the disease associated with the genes, diseases with genome frequency ≥1% were chosen and others were eliminated. Then, put the data into excel, plotting the network graph of the C-T network and D-T network with Cytoscape 3.7.0 software [24]. (3) C-T-D network: The merge function of Cytoscape software is used to generate the C-T network and the T-D network to construct a C-T-D network, that we used nodes to represent the compounds, targets, and diseases, and the lines between two nodes represented the interaction. We also used the degree to determine the size of each node and the betweenness to determine the thickness of lines and nodes and ensured that the graph would be obvious [25].

2.4. Gene Ontology and Pathway Analysis

The Comparative Toxicogenomics Database (http://ctdbase.org/, CTD) was adopted to conduct the Gene Ontology (GO) function analysis and the signaling pathway enrichment analysis. During this procedure, the significance level was set to 0.01, and organism was selected as Homo sapiens. The GO defines concepts related to gene function and the interrelationships among the functions of different genes [26]. It describes the functions of the candidate components in Fructus Gardeniae from Jiangxi province in terms of the molecular function, the cellular component involved, and the biological process affected [27]. In this study, pathway enrichment was additionally performed, for studying the biological effects and multidimensional pharmacological mechanism of Fructus Gardeniae from Jiangxi province at the pathway level. Then, an illustrated network that contains main “chemicals-targets-pathway (C-T-P)” of Fructus Gardeniae from Jiangxi province was constructed. The illustrated network that contains the main chemicals-targets-signaling pathway (C-T-P) of Fructus Gardeniae from Jiangxi province was established to understand their interaction. The results of the GO function analysis and the signaling pathway enrichment analysis were visualized via OmicShare platform2 (OmicShare, 2018).

3. Results

3.1. Target Identification of Candidate Components

The ADME parameters of candidate components obtained from TCMSP are shown in Table 1. Potential targets were predicted by SWISS servers. And based on the aforementioned target fishing approach, a total of 36 targets were predicted to interact with the 6 compounds identified in Fructus Gardeniae from Jiangxi province. The results of target name, gene name, UniProt ID, and matching probability are also shown in Table 1.

3.2. C-T-D Network Construction and Analysis

As shown in Figure 2, a chemicals-target-disease (C-T-D) network model graph was constructed with 74 nodes (6 candidate components, 36 potential targets, and 32 related diseases) and 416 edges. In the network, the yellow nodes (36) represent the potential targets related to the candidate compounds; the purple capsule-shaped nodes (6) represent the candidate compounds in Fructus Gardeniae from Jiangxi province; and the blue head-shaped nodes (32) represent the related diseases to the targets, while the edges represent the interactions between them. The node size varies according to degree, and edge thickness varies according to betweenness. Results of the network topology analysis are as follows: network density (0.168), network heterogeneity (0.527), and shortest paths (5402, 100%). The average degree of nodes is 12.27027, and there are 34 nodes larger than the average degree. The average betweenness centrality of nodes is 0.01763, and there are 26 nodes larger than the average betweenness centrality. The key core nodes (target or disease) are screened based on the topological properties of the degree of network nodes. The results showed that the eight targets of BCL2, MAPT, F2, BCL2L1, PRKCD, PRKCB, HIF1A, and PRKCA were in the front of degree ranking (degree > 20), and they were the pivotal nodes in the network, indicating that they might be the core targets of the pharmacological mechanism of candidate components; among the disease nodes, the nervous system diseases, cardiovascular diseases, mental disorders, neoplasms, brain diseases, endocrine system diseases, nutritional and metabolic diseases, digestive system diseases, chemically induced disorders, immune system diseases, metabolic diseases, and liver diseases rank ahead, suggesting that the candidate components may play a pharmacological role in these diseases.

3.3. Analysis of GO Enrichment

As highlighted in Figure 3, the gene ontology enrichment analysis consisted of three parts, BP (biological process), CC (cellular component), and MF (molecular function). Different categories of biological process, cellular component, and molecular function were represented by a purple, blue, and yellow bar, respectively. The height of the bar represented the number of genes observed in the category.

The enrichment results showed that there were 445 enrichment processes related to the biological processes, which cover response to the oxygen-containing compound, regulation of biological quality, response to organonitrogen compound, response to nitrogen compound, chemical synaptic transmission, and anterograde transsynaptic signaling; 73 enrichment results in the related items of cell composition, involving neuron projection, synapse part, neuron part, plasma membrane region, and postsynapse; 58 enrichment results are related to molecular function, which includes extracellularly the glycine-gated ion channel activity, the extracellularly glycine-gated chloride channel activity, the neurotransmitter receptor activity, protein phosphatase 2A binding, and chloride transmembrane transporter activity. Each value of enrichment results was calculated ( values < 0.01 were considered to be significantly enriched), ranking values according to the order from small to large. The top 20 enrichment results are displayed, and details are shown in Tables 24.

3.4. Analysis of Pathway

The pathway analysis result showed that 33 of the 36 (91.6%) potential targets were enriched and involved 181 signaling pathways, and 84 of these pathways were significantly correlated with the target genes (). The top 20 pathway with lower values and more genes enrichment are listed in Table 5, including signaling by GPCR, neuroactive ligand-receptor interaction, inflammatory mediator regulation of TRP channels, ion channel transport, dopaminergic synapse, amphetamine addiction, disinhibition of SNARE formation, EGFR tyrosine kinase inhibitor resistance, morphine addiction, and the HIF-1 signaling pathway. For the 20 listed pathways, 8 of them belonged to signal transduction, 3 belonged to human diseases, 3 belonged to transport of small molecules, 2 belonged to organismal systems, 2 belonged to hemostasis, 1 belonged to signaling molecules and interaction, while, the last 1 belonged to cellular processes. These pathways mainly involve liver diseases (liver cirrhosis and chemical- and drug-induced liver injury), neurological diseases (amphetamine-related disorders, cocaine-related disorders, autistic disorder, and schizophrenia), hypertension, neoplasms, hyperalgesia, inflammation, cardiomegaly, asthma, and diabetes mellitus. The illustrated network that contains the main chemicals-targets-signaling pathway (C-T-P) of Fructus Gardeniae from Jiangxi province was established to understand their interaction (Figure 4). The senior bubble map visually showed these significantly enriched pathways (Figure 5). The size and color of the nodes in the bubble graph were decided by the number of associated genes and the values. The size of the nodes indicated how many target genes are associated, and the colors from purple to yellow reflected the values from high to low.

4. Discussion

Bioactive components in TCM are the material basis to ensure its quality. From the biological point of view, the formation of the original plant varieties of TCM can be seen as the product of the interaction between genotype and habitat environment. The efficacy of TCM is closely related to its geographical origin. Newly Revised Materia Medica (Xin xiu ben cao), a materia medica commissioned by the government of the Tang Dynasty, stated: “if medicinal material is not produced from its native environment, the effect will be different.” That is to say, different geographical regions are suitable for growing different medicinal materials, and different medicinal materials also have different adaptability to the ecological environment of different producing areas. China has a vast territory, spanning tropical, subtropical, temperate, subfrigid, and plateau climate zones, and its terrain is complex. The conditions of sunshine, temperature, precipitation, and soil vary greatly in different regions, thus forming a complex and diverse ecosystem. There are also differences in the properties, chemical composition, and pharmacodynamics of traditional Chinese medicines bred in different ecosystems. The research on the correlation between the origin of TCM and the content of its chemical constituents has always been a worldwide scientific research hotspot in the field of traditional Chinese medicine resources. However, there are few experimental studies on pharmacodynamics of TCM from different geographical origins. The reason is that the chemical composition system of TCM is very complex. The contents of chemical constituents in medicinal materials from different habitats are different, and different chemical constituents can regulate the subjects through multiple targets and pathways, thus producing different pharmacological effects. Based on the above reasons, it is very difficult to carry out pharmacological experimental studies on traditional Chinese medicines from different origins because a large number of screening experiments need to be carried out from different geographical origins and different pharmacological effects. In the previous study, we collated the content data of Fructus Gardeniae from different geographical origins in scientific research literatures, and the content data were comprehensively evaluated by the single factor analysis of variance, principal component analysis, and cluster analysis to obtain the characteristic high-content bioactive components group in Fructus Gardeniae from different geographical origins. Based on network pharmacology, this study explored the key targets, signaling pathways, and disease networks of characteristic high-content chemical constituent groups of TCM from different geographical origins. The application of the network pharmacology method in the field of traditional Chinese medicine resources can solve the current research problems and provide powerful support for the research on pharmacological mechanism of traditional Chinese medicine from different geographical origins in China.

Preexperimental results demonstrated that among the 10 different provinces (Jiangxi, Sichuan, and Zhejiang, etc.) in China, the contents of six iridoid glycosides (geniposide, gardenoside, geniposidic acid, genipin 1-gentiobioside, gardoside, and shanzhiside) in Fructus Gardeniae from Jiangxi province were significantly higher than those in other regions (). Then, these six iridoid glycosides constitute the characteristic high-content chemical constituent group of Fructus Gardeniae from Jiangxi province. Iridoid glycosides have been demonstrated to be the major bioactive ingredients in Fructus Gardeniae [2830], among which, geniposide is the quality control standard recorded in the Chinese Pharmacopoeia (2000–2015 edition) [31], and gardenoside as a hydroxylation product of geniposide is also an important constituent in Fructus Gardeniae [5, 6]. This is similar to those recorded in ancient Chinese materia medica books. Jiangxi province is traditionally considered to be the “daodi” production region for Fructus Gardeniae (zhi zi), that is, Fructus Gardeniae from Jiangxi province is considered as “daodi medicinal materials” and is called “Jiang zhi zi” in Chinese. The term “daodi medicinal material” refers to a concept that has been widely recognized in the field of TCM for centuries. It is defined as “medicinal material that is produced and assembled in specific geographical regions with designated natural conditions and ecological environment [32]. These factors lead to quality and clinical effects surpass those of same botanical origin produced from other geographical regions, and thus, it is widely recognized and has long enjoyed a good reputation [33].” Although it is also known as genuine medicinal material, geo-authentic medicinal material, authentic and superior medicinal herbal, authentic medicinal, and geoherb [3436]. Based on our previous experimental results, this study selected the six characteristic high-content bioactive components of Fructus Gardeniae from Jiangxi province as candidate components, using Swiss Target Prediction to obtain 36 potential targets related to candidate components, and imported 18 targets into CTD database to predict target-related diseases. 32 diseases with a genome frequency ≥1% were selected to construct the C-T-D network. The nodes in the network ranked according to degree. The results showed that the eight targets of BCL2, MAPT, F2, BCL2L1, PRKCD, PRKCB, HIF1A, and PRKCA were in the front of degree ranking (degree >20), and among the disease nodes, the nervous system diseases, cardiovascular diseases, mental disorders, neoplasms, brain diseases, endocrine system diseases, nutritional and metabolic diseases, digestive system diseases, chemically induced disorders, immune system diseases, metabolic diseases, and liver diseases rank ahead, suggesting that characteristic high-content bioactive components group (geniposide, gardenoside, geniposidic acid, genipin 1-gentiobioside, gardoside, and shanzhiside) of Fructus Gardeniae from Jiangxi province may play a pharmacological role in these diseases through the potential targets such as BCL2, MAPT, F2, BCL2L1, and PRKCD. That is to say, in the field of traditional Chinese medicine resources, if we want to carry out pharmacodynamic experimental research on Fructus Gardeniae from Jiangxi province, we should start with antinervous diseases, anticardiovascular diseases, antimental disorders, antineoplasms, antibrain diseases, and antiendocrine system diseases; if we want to carry out the related pharmacological mechanism experimental research, we can start with BCL2, MAPT, F2, BCL2L1, PRKCD, PRKCB, HIF1A, PRKCA, and other 25 potential targets.

In this study, the GO enrichment analysis and the pathway enrichment analysis were carried out in order to explore the multidimensional pharmacological mechanism of the characteristic high-content bioactive components group of Fructus Gardeniae from Zhejiang province. The results showed that the pathways with higher number of related annotated genes count and lower the value mainly included signaling by GPCR, neuroactive ligand-receptor interaction, inflammatory mediator regulation of TRP channels, ion channel transport, dopaminergic synapse, amphetamine addiction, disinhibition of SNARE formation, EGFR tyrosine kinase inhibitor resistance, morphine addiction, and HIF-1 signaling pathway. Among these pathways, the GPCR pathway, platelet activation signaling-and-aggregation, and ion channel transport are mainly associated with liver cirrhosis and hypertension; neuroactive ligand-receptor interaction, dopamine receptors, dopaminergic synapse, amphetamine addiction, and morphine addiction are mainly associated with neurological diseases such as autism, schizophrenia, and cocaine-related diseases; EGFR tyrosine kinase inhibitor resistance, gap junction, and the HIF-1 signaling pathway are mainly associated with neoplasms and liver cirrhosis; inflammatory mediator regulation of TRP channels is mainly associated with hyperalgesia, inflammation, and liver cirrhosis; ligand-gated ion channel transport is mainly associated with autistic disorder and brain injuries. From the above analysis, we find an interesting phenomenon that among these diseases related to high-enrichment pathways, neurological diseases (including amphetamine-related disorders, cocaine-related disorders, autistic disorder, and schizophrenia.) and liver diseases (including liver cirrhosis, carcinoma hepatocellular, and chemical- and drug-induced liver injury.) seem to present a particularly high frequency. These results make us naturally associate with a complex disease, which is “hepatic encephalopathy”. Hepatic encephalopathy (HE) is a known neurologic complication of advanced cirrhosis [37]. Between 30% and 50% of hospitalization cases of cirrhosis are related to HE [38]. HE is a syndrome of spectrum of neuropsychiatric abnormalities caused by portosystemic venous shunting, with intrinsic liver disease. Related neurocomorbidities include depression, bipolar disorder, schizophrenia, drug abuse, and suicidal tendencies [39]. On the basis of the results obtained in this study, we boldly speculate that the Fructus Gardeniae from Jiangxi province can play a pharmacological role in hepatic encephalopathy through regulating multiple signaling pathways in an integrated manner.

This study reveals the related targets, key biological pathways, and main disease types of the characteristic high-content chemical constituent group of Fructus Gardeniae from Jiangxi province from the perspective of network pharmacology. In conclusion, our study applied the network pharmacology method to the field of traditional Chinese medicine resources innovatively, and it laid a foundation for the study of multidimensional pharmacological mechanism of Fructus Gardeniae from different geographical regions in China.

Data Availability

The data used to support the findings of our research are included within the article and the supplementary information files.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors thank the National Science and Technology Basic Project of the Ministry of Science and Technology of China (no. 2015FY111500-190) and the Independent Subject of China Academy of Chinese Medical Sciences (no. ZZ130311).

Supplementary Materials

Supplementary 1. Table S1: the analysis data of chemicals-target (C-T) network. Supplementary 2. Table S2: the analysis data of target-disease (T-D) network. Supplementary 3. Table S3: the analysis data of targets-signaling pathway (T-P) network. (Supplementary Materials)