# Mechanism of Jujube (Ziziphus jujuba Mill.) Fruit in the Appetite Regulation Based on Network Pharmacology and Molecular Docking Method

Authors: Yu Zhu; Ju Huang; Tao Shen; Rensong Yue
Journal: Contrast Media & Molecular Imaging
 (2022)

---

## Abstract

Objective. To investigate the mechanism of jujube (Ziziphus jujuba Mill.) in appetite regulation based on network pharmacology. Methods. The active components and action targets of jujube were retrieved through the TCMSP and TCMID databases. GeneCards, DisGeNet, Therapeutic Target Database, and OMIM were used to screen the related targets for appetite, appetite suppression, and appetite regulation, and the intersection target of the two was selected. A protein-protein interaction (PPI) network was constructed. Important protein nodes and subnets were predicted based on the cytoHubba plug-in, and the hub gene was screened. Additionally, GO and KEGG pathway analyses were performed to obtain potential biological processes and signaling pathways of key targets. And the active ingredient-target-action pathway diagram was constructed. Results. A total of 16 active components were screened from jujube, including 131 action targets related to appetite and appetite regulation. Three key targets (MAOA, MMP2, and HSPB1) were screened out by MCODE analysis. KEGG enrichment analysis was mainly enriched in neuroactive ligand-receptor interaction, serotonin-containing synapse, gap junction, cAMP signaling pathway, and dopaminergic synapse. Molecular docking results showed that the components coclaurine, (−)-catenin, (+)-stepholidine, berberine, cianidanol, coclaurine, and moupinamide in jujube had strong binding activity to the main targets (ESR1, ADRA2C, and MMP2). Conclusion. Based on network pharmacology, the appetite modulating effects of jujube on multiple components, targets, and channels were explored, and the main active components of jujube were predicted to act on multiple signaling pathways to regulate appetite. The molecular docking results showed that the components in jujube had strong binding activity to the main targets, which provided new ideas and methods to further investigate the mechanisms of appetite regulation by jujube.


---

## Body

## 1. Introduction

Hyperphagia is one of the most common and intractable symptoms in diabetic patients and is an important culprit in disrupting glucose homeostasis [1]. Studies have shown that recurrent blood glucose fluctuations are an independent risk factor for increased diabetic complications and cardiovascular mortality, and hyperphagia is a key driver of blood glucose drift and disease progression. The existing diabetes treatment drugs GLP-1 agonists can play a role in suppressing appetite in diabetic patients by inhibiting the appetite center and delaying gastric emptying [2, 3]. However, these drugs are expensive and can cause varying degrees of vomiting, headache, nasopharyngitis, and significant weight loss.Fatty and sweet foods can trap the spleen and stomach, which will lead to spleen deficiency for a long time and eventually become diabetic, manifesting as easy hunger. Supplementing spleen deficiency and responding to the sweet nature of the spleen is one of the main tools in the clinical treatment of hyperphagia in diabetes in Chinese medicine [4, 5]. Spleen deficiency and overflow of Qi are the keys to the pathogenesis of hyperphagia in type 2 diabetes (T2DM), and the main treatment rule is to “conform to the preference of the sick, use sweet herbs to treat the spleen” [4]. Clinical practice has shown that large doses of jujube (Ziziphus jujuba Mill.) decoction can produce the symptoms of gastric fullness in patients and then achieve good results in appetite control without producing more obvious blood glucose fluctuations, but can also assist other glucose-lowering treatment options to maintain blood glucose homeostasis [4, 5]. Some basic studies suggest that jujube components have antioxidant [6], improving insulin resistance [7], inducing cancer cell apoptosis [8], inhibiting α-glucosidase activity [9], and sedative-hypnotic effects [10]. However, there is a relative lack of research on the pharmacological mechanisms of action of jujube on appetite regulation.Network pharmacology is a research method that uses high-throughput screening in databases, network visualization, and data analysis techniques to reveal the complex biological network relationships among drugs, targets, and diseases and to analyze and predict the pharmacological mechanisms of drugs [11]. Molecular docking is a powerful tool for predicting the affinity and binding mode of proteins and ligands. For a given protein and ligand (protein, DNA/RNA, or small molecule), their binding mode and free energy of binding can be predicted to investigate their functions and mechanisms of action [12]. Virtual screening based on molecular docking methods has become one of the required processes for drug development against specific target proteins [12]. By constructing a database based on the main components of jujube and establishing a Chinese medicine target dataset through target prediction, then constructing a compound-target network, constructing a protein-protein interaction (PPI) network [13], performing gene ontology (GO) function enrichment analysis [14], and performing pathway-based enrichment analysis, we can explain the effects of jujube on appetite regulation at the molecular level by constructing a multidimensional network of Chinese medicine chemical composition-target of action-disease target-PPI network. By constructing a multidimensional network of Chinese herbal chemical constituents-targets-disease targets-PPI network, we can explain the effects of jujube on appetite regulation at the molecular level.


## 2. Materials and Methods

### 2.1. Composition of Jujube

The Chinese medicine name jujube was used to obtain drug composition information using the TCMSP (https://tcmspw.com/tcmsp.php) [15] database, TCMID (https://www.megabionet.org/tcmid/) [16] database, and other Chinese medicine composition databases, including the number of ingredients, molecular name, and molecular mass. The structures were imported into PubChem (https://pubchem.ncbi.nlm.nih.gov/) [17] for searching and normalization, supplementing their PubChem CID, and downloading the SDF structures.


### 2.2. Screening of Active Ingredients in Jujube

ADMET [18] is the absorption, distribution, metabolism, excretion, and toxicity of a drug. ADME is the study of drug metabolism kinetics and is commonly used in contemporary drug design and screening. In this study, we used the ADMET Descriptors module of Discovery Studio 2017R2 to predict the ADMET parameters of herbal ingredients, based on which human intestinal absorption (ADMET_Absorption_Level) and ADMET-Aqueous Solubility was used for the screening of the Chinese herbal ingredients. Compounds with ADMET_Absorption_Level 0, 1, 2 and ADMET_Solubility_Level 1, 2, 3, 4 were selected for inclusion in the study.


### 2.3. Prediction Screening of Target Proteins of Active Ingredients of Jujube

The predictive screening of the active ingredient target proteins of jujube was performed using the following two steps: the smile structures of the screened active ingredients were entered into DrugBank (https://go.drugbank.com/) [19], Therapeutic Target Database (https://db.idrblab.net/ttd/) [20], and Swiss Target Prediction platform (https://www.swisstargetprediction.ch/) [21] to predict the relevant targets of the active ingredients of the herbs and select “Homo sapiens” for the screening. The targets of “Homo sapiens” species were selected for screening, and a database of active ingredient targets of jujube ingredients was constructed.The target prediction was carried out according to the method of Fu et al. [22], and all targets in the target database were scored by deep learning and a Bayesian network algorithm. The network topology parameters were calculated according to the scores, and the targets of the active ingredients of jujube were screened for subsequent research.


### 2.4. Disease Target Screening

The search term “Appetite Depressants/Appetite/Appetite Regulation” was set and the GeneCards database (https://www.GeneCards.org/) [23], the DisGeNet database (https://www.disgenet.org/) [24], the Therapeutic Target Database, and the OMIM database (https://www.omim.org) [25] were used to obtain the disease-related targets. The GeneCards database was applied to screen genes with a score greater than 2.3, and the DisGeNet database was queried for genes derived from the CTD (https://ctdbase.org/)-human [26] database. The OMIM and Therapeutic Target Database were applied to collect disease-associated genes. The data obtained from the four databases were merged to take the intersection, and the duplicate or invalid genes were removed to build the disease target database.


### 2.5. PPI Network Construction and Screening of Hub Genes

STRING database [27] is a database for searching known proteins and predicting protein-protein interactions, which contains 2031 species containing 9.6 million proteins and 138 million protein-protein interactions. We used the STRING database to construct a PPI network for the intersection of compounds and disease targets. The intersection was taken for the targets of the active ingredient of the compound and the disease targets. The intersection targets were uploaded to the STRING database (https://string-db.org/). The information of the constructed PPI network was imported into Cytoscape 3.8.2 software [28], based on the cytoHubba [29] plug-in topological algorithm to predict the important protein nodes and subnetworks in the network. And this study used DEGREE (Degree Correlation), MNC (Maximum Neighborhood Component), MCC (Maximal Clique Centrality), EPC (Edge Percolated Component), CLONESS (Closeness Centrality), and visualized hub gene. Cluster analysis was performed using the plug-in MCODE [30] to find gene clusters by using the correlation between proteins in the network, derive subnetworks, extract the differential genes contained in each gene cluster, and analyze the subnetworks. The main biological processes of the targets in the subnetworks are analyzed by finding gene clusters using correlations between proteins in the network, deriving subnetworks, and extracting the differential genes contained in each gene cluster.


### 2.6. Functional Enrichment and Disease Enrichment Analysis

GO functional annotation analysis is a common approach to perform large-scale functional enrichment studies of genes, including biological process (BP), molecular function (MF), and cellular component (CC). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway is a widely used database for storing information about genomes, biological pathways, diseases, and drugs. The screened hub genes were imported into the DAVID 6.8 database (https://david.ncifcrf.gov/) [31], and the species selection “Homo sapiens” was used for GO analysis and KEGG pathway analysis (P<0.05) to analyze the key targets involved in the relevant biological processes. Signaling pathways were analyzed and visualized using Cytoscape software.


### 2.7. Construction of “Active Ingredient-Potential Target-Action Pathway” Network

Cytoscape 3.8.2 software [28] was used to construct the “active ingredient-potential target-action pathway” network. The network consists of three parts, namely, active ingredients, target proteins, and pathways, to analyze and explore the multicomponent-multitarget-multipathway mechanism of Chinese medicine for the treatment of diseases.


### 2.8. Molecular Docking of Key Targets and Components

Molecular docking was performed using CDOCK in the Receptor-Ligand International module of the software Discovery Studio 2017 R2 [32] to precisely dock and analyze the key targets and the main components of Chinese medicine. The 3D structures of the small molecule compounds of the main active ingredients of traditional Chinese medicine were downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/) according to their PubChem_ID numbers and imported into Discovery Studio 2017 R2. The high-resolution crystal structures of the targets were downloaded from the PDB [33] (https://www.rcsb.org/pdb/home/home.do) protein database, and the active sites of the proteins were centered on the active amino acid sites of the original ligand action labeled in the crystal structure itself, and the corresponding “active pockets” were constructed, so that the system searches for “active pockets” near the active site and finally locates the “active pocket” information to the target “active pocket.”The parameters of the CDOCKER algorithm [34] module were set as follows: Pose Cluster Radius was set to 0.5, Random Conformations was set to 10, Orientations to Refine was set to 10, and the rest of the default parameters were kept unchanged. The process is shown in Figure 1.Figure 1
Flowchart.


## 2.1. Composition of Jujube

The Chinese medicine name jujube was used to obtain drug composition information using the TCMSP (https://tcmspw.com/tcmsp.php) [15] database, TCMID (https://www.megabionet.org/tcmid/) [16] database, and other Chinese medicine composition databases, including the number of ingredients, molecular name, and molecular mass. The structures were imported into PubChem (https://pubchem.ncbi.nlm.nih.gov/) [17] for searching and normalization, supplementing their PubChem CID, and downloading the SDF structures.


## 2.2. Screening of Active Ingredients in Jujube

ADMET [18] is the absorption, distribution, metabolism, excretion, and toxicity of a drug. ADME is the study of drug metabolism kinetics and is commonly used in contemporary drug design and screening. In this study, we used the ADMET Descriptors module of Discovery Studio 2017R2 to predict the ADMET parameters of herbal ingredients, based on which human intestinal absorption (ADMET_Absorption_Level) and ADMET-Aqueous Solubility was used for the screening of the Chinese herbal ingredients. Compounds with ADMET_Absorption_Level 0, 1, 2 and ADMET_Solubility_Level 1, 2, 3, 4 were selected for inclusion in the study.


## 2.3. Prediction Screening of Target Proteins of Active Ingredients of Jujube

The predictive screening of the active ingredient target proteins of jujube was performed using the following two steps: the smile structures of the screened active ingredients were entered into DrugBank (https://go.drugbank.com/) [19], Therapeutic Target Database (https://db.idrblab.net/ttd/) [20], and Swiss Target Prediction platform (https://www.swisstargetprediction.ch/) [21] to predict the relevant targets of the active ingredients of the herbs and select “Homo sapiens” for the screening. The targets of “Homo sapiens” species were selected for screening, and a database of active ingredient targets of jujube ingredients was constructed.The target prediction was carried out according to the method of Fu et al. [22], and all targets in the target database were scored by deep learning and a Bayesian network algorithm. The network topology parameters were calculated according to the scores, and the targets of the active ingredients of jujube were screened for subsequent research.


## 2.4. Disease Target Screening

The search term “Appetite Depressants/Appetite/Appetite Regulation” was set and the GeneCards database (https://www.GeneCards.org/) [23], the DisGeNet database (https://www.disgenet.org/) [24], the Therapeutic Target Database, and the OMIM database (https://www.omim.org) [25] were used to obtain the disease-related targets. The GeneCards database was applied to screen genes with a score greater than 2.3, and the DisGeNet database was queried for genes derived from the CTD (https://ctdbase.org/)-human [26] database. The OMIM and Therapeutic Target Database were applied to collect disease-associated genes. The data obtained from the four databases were merged to take the intersection, and the duplicate or invalid genes were removed to build the disease target database.


## 2.5. PPI Network Construction and Screening of Hub Genes

STRING database [27] is a database for searching known proteins and predicting protein-protein interactions, which contains 2031 species containing 9.6 million proteins and 138 million protein-protein interactions. We used the STRING database to construct a PPI network for the intersection of compounds and disease targets. The intersection was taken for the targets of the active ingredient of the compound and the disease targets. The intersection targets were uploaded to the STRING database (https://string-db.org/). The information of the constructed PPI network was imported into Cytoscape 3.8.2 software [28], based on the cytoHubba [29] plug-in topological algorithm to predict the important protein nodes and subnetworks in the network. And this study used DEGREE (Degree Correlation), MNC (Maximum Neighborhood Component), MCC (Maximal Clique Centrality), EPC (Edge Percolated Component), CLONESS (Closeness Centrality), and visualized hub gene. Cluster analysis was performed using the plug-in MCODE [30] to find gene clusters by using the correlation between proteins in the network, derive subnetworks, extract the differential genes contained in each gene cluster, and analyze the subnetworks. The main biological processes of the targets in the subnetworks are analyzed by finding gene clusters using correlations between proteins in the network, deriving subnetworks, and extracting the differential genes contained in each gene cluster.


## 2.6. Functional Enrichment and Disease Enrichment Analysis

GO functional annotation analysis is a common approach to perform large-scale functional enrichment studies of genes, including biological process (BP), molecular function (MF), and cellular component (CC). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway is a widely used database for storing information about genomes, biological pathways, diseases, and drugs. The screened hub genes were imported into the DAVID 6.8 database (https://david.ncifcrf.gov/) [31], and the species selection “Homo sapiens” was used for GO analysis and KEGG pathway analysis (P<0.05) to analyze the key targets involved in the relevant biological processes. Signaling pathways were analyzed and visualized using Cytoscape software.


## 2.7. Construction of “Active Ingredient-Potential Target-Action Pathway” Network

Cytoscape 3.8.2 software [28] was used to construct the “active ingredient-potential target-action pathway” network. The network consists of three parts, namely, active ingredients, target proteins, and pathways, to analyze and explore the multicomponent-multitarget-multipathway mechanism of Chinese medicine for the treatment of diseases.


## 2.8. Molecular Docking of Key Targets and Components

Molecular docking was performed using CDOCK in the Receptor-Ligand International module of the software Discovery Studio 2017 R2 [32] to precisely dock and analyze the key targets and the main components of Chinese medicine. The 3D structures of the small molecule compounds of the main active ingredients of traditional Chinese medicine were downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/) according to their PubChem_ID numbers and imported into Discovery Studio 2017 R2. The high-resolution crystal structures of the targets were downloaded from the PDB [33] (https://www.rcsb.org/pdb/home/home.do) protein database, and the active sites of the proteins were centered on the active amino acid sites of the original ligand action labeled in the crystal structure itself, and the corresponding “active pockets” were constructed, so that the system searches for “active pockets” near the active site and finally locates the “active pocket” information to the target “active pocket.”The parameters of the CDOCKER algorithm [34] module were set as follows: Pose Cluster Radius was set to 0.5, Random Conformations was set to 10, Orientations to Refine was set to 10, and the rest of the default parameters were kept unchanged. The process is shown in Figure 1.Figure 1
Flowchart.


## 3. Results

### 3.1. Screening of Active Compounds in Jujube

A total of 21 chemical components were obtained from the database collected for the jujube (DZ), with the main structural types being flavonoids, phenylpropanoids, alkaloids, terpenoids, etc. Chinese medicine contains a large number of chemical components, and the DS software was used to predict the ADMET parameters of the chemical components contained in the compound based on their chemical structures, which helps to find the possible active components quickly. Finally, a total of 16 active ingredients were screened, and the corresponding information about the screened active ingredients is detailed in Table1.Table 1
Chemical composition and ADME parameters of jujube.
No.CompoundPubChem CIDADMET_Absorption_LevelADMET_Solublity_LevelHerbDZ1(−)-Catechin7316003JujubeDZ2(+)-Stepholidine1244299902JujubeDZ3Berberine235302JujubeDZ4Betulinic acid6497121JujubeDZ5Ceanothic acid16135221JujubeDZ6Cianidanol906403JujubeDZ7Coclaurine16048703JujubeDZ8Coumestrol528170702JujubeDZ9Malkangunin9047315503JujubeDZ10Mauritine D644302612JujubeDZ11Moupinamide528053703JujubeDZ12Nuciferine1014602JujubeDZ13Protopine497002JujubeDZ14Quercetin528034313JujubeDZ15Spiradine A44175603JujubeDZ16Stepharine9845503Jujube


### 3.2. Screening of Appetite-Related Targets in Jujube

The targets withp value >0.9 were screened from the prediction results of the target database as active ingredient targets, and a total of 194 active ingredient targets were obtained from the Chinese herbal compound. A larger number of targets were screened in the GeneCards database according to the search term “Appetite,” and 3330 targets were selected based on the criterion of score >2.3. 224 targets were obtained from OMIM database, and 275 targets were obtained from DisGeNEt. A total of 3552 targets were obtained after combining and deweighting, with 275 targets in OMIM database, 275 targets in DisGeNEt database, and 1 target in TTD database. The 194 potential targets of Chinese herbal ingredients were intersected with 3552 targets of disease targets, and a Wayne diagram was drawn (Figure 2(a)), and 131 potential targets of Chinese herbal compounds were initially obtained, and a compound name-Chinese herbal medicine-drug target interaction network was constructed (Figure 2(b)).Figure 2
Interaction network diagram of jujube for appetite treatment. (a) Venn diagram of jujube action targets and disease targets; (b) jujube-component-target interactions network. Green hexagon is Chinese medicine, red circle is component, and blue diamond is target.
(a)(b)


### 3.3. Construction of the PPI Network and Screening of Key Targets

The screened potential targets were input into the STRING database to obtain the target protein interaction information and imported into Cytoscape to construct the PPI network (Figure3(a)), which had 130 nodes (target proteins) and 1004 edges (protein interactions). It indicated that among the predicted disease-related targets, the more targets could have effective interactions with that target. Using the five parameters of MNC, DEGREE, MCC, CLONESS, and EPC for screening (Figure 3(b)), the algorithm’s computational analysis of the network structure and weighted linkage between nodes could screen out important key genes. The intersection of the top 30 results of each algorithm was taken to obtain 11 key targets (Figure 3(c), Tables 2 and 3).Figure 3
Protein interaction network diagram. (a) All target protein interaction network; (b) top 30 target Venn diagram of MNC, DEGREE, MCC, CLONESS, and EPC; (c) key target protein interaction network diagram.
(a)(b)(c)Table 2
cytoHubba five-algorithm calculation results ranking.
RankMCCCLONESSMNCEPCDEGREE1MAOAMAPK3MAPK3MAPK3MAPK32DRD2EGFREGFRSRCEGFR3HTR1ASRCSRCEGFRSRC4DRD4HSP90AA1HSP90AA1APPHSP90AA15DRD3ESR1TNFHSP90AA1TNF6SLC6A3TNFESR1ESR1ESR17SLC18A2APPAPPPIK3CAAPP8ADRA2CPIK3CAPIK3CATNFPIK3CA9ADRA2ABRCA1BRCA1RELABRCA110HTR3ARELAGRIN2BBRCA1GRIN2B11HTR1BGSK3BSLC6A3GSK3BSLC6A312CHRNA4GRIN2BHTR3AGRIN2BHTR3A13ACHEMAPTRELANTRK2RELA14PIK3CAPRKCACHRNA4NFKB1CHRNA415HTR2APPARGACHECHRNA4ACHE16ADRA1BNTRK2CDK1DRD2GSK3B17ADRA1ANFKB1DRD2SLC6A3CDK118ADRA1DCDK1MAOAPRKCADRD219HTR2CDRD2NTRK2ACHEMAOA20HTR2BACHEGSK3BMAPTPPARG21APPHIF1APPARGIGF1RNTRK222OXTRCHRNA4HIF1ACDK1HIF1A23F2CDK5MAPTHIF1ASLC18A224NTSR1GRIN1MMP2HTR3AGRIN125SRCIGF1RTOP2ASLC18A2MAPT26EGFRABCB1NFKB1PPARGMMP227MAPK3PTPN1SLC18A2ADRA1BTOP2A28HSP90AA1MMP2GRIN1MAOANFKB129ESR1TOP2APRKCAMCL1PRKCA30RELAF2DRD4GRIN1DRD4Table 3
Analysis of topological parameters of key targets.
NameClosenessBetweennessDegreeACHE0.4867920.01486926APP0.5657890.06205839CHRNA40.4813430.01189227DRD20.4904940.01091126EGFR0.6113740.07756852ESR10.5707960.05457942HSP90AA10.5863640.07245646MAPK30.6417910.14001158PIK3CA0.560870.05258838RELA0.5330580.01169727SRC0.6056340.07343649


### 3.4. Subnetwork Analysis

MCODE subnetwork analysis can discover more closely connected groups or genes in the network. It is calculated by weighting the points with the highest weight and set as SEED, from SEED, recursively move outward to find nodes that can join the subnetwork. Subnetwork 1 is centered on MAOA, and the important targets connected to it, such as SLC6A3, DRD3, DRD4, SLC18A2, and HTR1A, are all dopaminergic synapse-related targets, indicating that subnetwork 1 is closely related to dopaminergic synapses (Figure4(a)). The core of subnetwork 2 is MMP2, and the important targets connected with it, such as TNF, are all targets related to inflammatory processes, indicating that subnetwork 2 is closely related to the regulation of tryptophan channels by inflammatory mediators (Figure 4(b)). The core of subnetwork 3 is HSPB1, and the important targets connected with it, such as CHRNA4, BRCA1, and HTR3A, are all targets related to cAMP signaling pathway, indicating that subnetwork 3 is closely related to cAMP signaling channels (Figure 4(c)).Figure 4
Subnetwork diagram. (a) Subnetwork 1; (b) subnetwork 2; (c) subnetwork 3. The larger the node, the darker the color means the higher the degree value of the target point.
(a)(b)(c)


### 3.5. GO Enrichment Analysis

In order to explore the functional distribution of key targets, 131 key targets were entered into the DAVID 6.8 database for GO enrichment analysis. The results showed 343 biological processes, 103 molecular functions, and 58 cellular components. Combined with the literature, the key targets in the biological process were filtered byP<0.5 and the number of enriched targets was high, and the key targets were concentrated in response to stimulus, signaling, cell proliferation, positive regulation of biological process, negative regulation of biological process, etc. (Figure 5(a) and Table 4). Among the molecular functions, molecular transducer activity, catalytic activity, transporter activity, and transcriptional regulation activity are mainly involved (Figure 5(b) and Table 4). Among the cellular components, the membrane is the most involved target, followed by synapse, organelle part, and cell junctions (Figure 5(c) and Table 4).Figure 5
GO enrichment analysis of the key targets. (a) Biological process (BP) analysis result; (b) molecular functions (MF) analysis result; (c) cellular components (CC) analysis result.
(a)(b)(c)Table 4
GO analysis table.
ClassGOTermCountP valueMolecular functionGO:0060089Molecular transducer activity542.19E − 24GO:0003824Catalytic activity744.44E − 09GO:0005215Transporter activity219.23E − 05GO:0005488Binding1283.18E − 03GO:0104005Hijacked molecular function31.70E − 02GO:0016209Antioxidant activity32.42E − 02GO:0098772Molecular function regulator204.35E − 02GO:0140110Transcription regulator activity172.23E − 01GO:0005198Structural molecule activity47.77E − 01Cellular componentGO:0045202Synapse495.58E − 23GO:0044456Synapse part343.96E − 15GO:0030054Cell junction302.44E − 08GO:0044425Membrane part814.76E − 08GO:0032991Protein-containing complex682.16E − 07GO:0016020Membrane966.75E − 07GO:0031974Membrane-enclosed lumen661.02E − 06GO:0044422Organelle part962.78E − 05GO:0005623Cell1301.52E − 04GO:0044464Cell part1301.52E − 04Biological processGO:0050896Response to stimulus1253.86E − 29GO:0023052Signaling1081.98E − 26GO:0032501Multicellular organismal process1104.81E − 22GO:0048518Positive regulation of biological process995.94E − 21GO:0007610Behavior332.62E − 20GO:0008283Cell proliferation531.19E − 18GO:0048511Rhythmic process235.07E − 17GO:0048519Negative regulation of biological process885.95E − 17GO:0032502Developmental process914.46E − 15GO:0065007Biological regulation1267.53E − 15


### 3.6. KEGG Pathway Analysis

The KEGG pathway enrichment analysis of potential targets by the DAVID 6.8 data platform (P<0.05) is shown in Figures 6(a)–6(d). The top 10 pathways were Neuroactive ligand-receptor interaction, Serotonergic synapse , Gap junction, cAMP signaling pathway, Dopaminergic synapse, Calcium signaling pathway, Hypoxia-inducible factor pathway (HIF-1 signaling pathway), Prolactin signaling pathway, Thyroid hormone signaling pathway, and Inflammatory mediator regulation of TRP channels (Table 5). It is suggested that jujube components may exert appetite modulating effects through the above pathways.Figure 6
KEGG pathway enrichment results. (a) KEGG pathway network diagram; (b) KEGG enrichment pathway annotated classification results; (c) KEGG pathway enrichment circle diagram; (d) KEGG enrichment result bar graph.
(a)(b)(c)(d)Table 5
KEGG pathway enrichment information.
PathwayCountP valueNeuroactive ligand-receptor interaction301.28E − 15Serotonergic synapse184.07E − 14Gap junction153.49E − 12cAMP signaling pathway193.98E − 10Dopaminergic synapse141.04E − 08Calcium signaling pathway186.88E − 08HIF-1 signaling pathway109.55E − 06Prolactin signaling pathway81.60E − 05Thyroid hormone signaling pathway101.68E − 05Inflammatory mediator regulation of TRP channels91.73E − 05


### 3.7. “Active Ingredient-Key Target-Pathway” Network Construction for Jujube

The active ingredients, potential targets, and selected signaling pathways of jujube were imported into Cytoscape 3.8.2 software to construct the “jujube-active-ingredient-target-action pathway” diagram (Figure7(a)). The MCC algorithm of cytoHubba was used to further calculate the closest association of each component with key targets (Figure 7(b)). The results concluded that 5-hydroxytryptamine (5-HT)-containing synapses was the key pathway of action.Figure 7
Dates-component-target-pathway diagram. (a) Dazao-ingredient-target-pathway diagram; (b) the key target-ingredient diagram. The green hexagon is the name of a single herbal medicine, the red circle is the ingredient, the blue diamond is the target, and the purple arrow is the pathway name.
(a)(b)


### 3.8. Molecular Docking Simulation of Target Interaction with Related Compounds

The docking of key target genes and their related compounds revealed that the docking binding energies of the targets and their related compounds were all negative (Table6), indicating that the related compounds bind well to the targets. Among them, coclaurine had the lowest docking binding energy of −8.38 with ESR1, (−)-Catechin and ESR1, (+)-Stepholidine and ADRA2C, Berberine and ADRA2C, Cianidanol and ESR1, Coclaurine and ESR1, Moupinamide and MMP2, Quercetin and SCR, and Quercetin and TNF all had binding energies <−5 kcal/mol, and the binding patterns are shown in Figures 8(a)–8(h).Table 6
Molecular docking results of key targets and their related compounds.
CompoundTargetCombined energy (kcal/mol)(−)-CatechinESR1−7.01ADRA2C−5.47(+)-StepholidineADRA2C−7.10BerberineADRA2C−7.80CianidanolESR1−7.14ADRA2C−5.40CoclaurineESR1−7.02ADRA2C−4.82CoumestrolESR1−8.38MoupinamideMMP2−7.87EGFR−4.03NuciferineADRA2C−6.79QuercetinTNF−7.61SRC−6.95MMP2−6.72MAPK3−5.24HSP90AA1−6.25ESR1−7.36EGFR−4.83APP−4.57Figure 8
Molecular docking simulations of targets and compounds. (a) 3D and 2D plots of molecular docking of (−)-catechin-ESR1; (b) 3D and 2D plots of molecular docking of (+)-stepholidine-ADRA2C; (c) molecular docking 3D and 2D diagrams of berberine-ADRA2C; (d) molecular docking 3D and 2D diagrams of cianidanol-ESR1; (e) molecular docking 3D and 2D plots of coclaurine-ESR1; (f) molecular docking 3D and 2D plots of moupinamide-MMP2 target; (g) molecular docking 3D and 2D diagrams of quercetin-SCR; (h) molecular docking 3D and 2D diagrams of quercetin-TNF.
(a)(b)(c)(d)(e)(f)(g)(h)


## 3.1. Screening of Active Compounds in Jujube

A total of 21 chemical components were obtained from the database collected for the jujube (DZ), with the main structural types being flavonoids, phenylpropanoids, alkaloids, terpenoids, etc. Chinese medicine contains a large number of chemical components, and the DS software was used to predict the ADMET parameters of the chemical components contained in the compound based on their chemical structures, which helps to find the possible active components quickly. Finally, a total of 16 active ingredients were screened, and the corresponding information about the screened active ingredients is detailed in Table1.Table 1
Chemical composition and ADME parameters of jujube.
No.CompoundPubChem CIDADMET_Absorption_LevelADMET_Solublity_LevelHerbDZ1(−)-Catechin7316003JujubeDZ2(+)-Stepholidine1244299902JujubeDZ3Berberine235302JujubeDZ4Betulinic acid6497121JujubeDZ5Ceanothic acid16135221JujubeDZ6Cianidanol906403JujubeDZ7Coclaurine16048703JujubeDZ8Coumestrol528170702JujubeDZ9Malkangunin9047315503JujubeDZ10Mauritine D644302612JujubeDZ11Moupinamide528053703JujubeDZ12Nuciferine1014602JujubeDZ13Protopine497002JujubeDZ14Quercetin528034313JujubeDZ15Spiradine A44175603JujubeDZ16Stepharine9845503Jujube


## 3.2. Screening of Appetite-Related Targets in Jujube

The targets withp value >0.9 were screened from the prediction results of the target database as active ingredient targets, and a total of 194 active ingredient targets were obtained from the Chinese herbal compound. A larger number of targets were screened in the GeneCards database according to the search term “Appetite,” and 3330 targets were selected based on the criterion of score >2.3. 224 targets were obtained from OMIM database, and 275 targets were obtained from DisGeNEt. A total of 3552 targets were obtained after combining and deweighting, with 275 targets in OMIM database, 275 targets in DisGeNEt database, and 1 target in TTD database. The 194 potential targets of Chinese herbal ingredients were intersected with 3552 targets of disease targets, and a Wayne diagram was drawn (Figure 2(a)), and 131 potential targets of Chinese herbal compounds were initially obtained, and a compound name-Chinese herbal medicine-drug target interaction network was constructed (Figure 2(b)).Figure 2
Interaction network diagram of jujube for appetite treatment. (a) Venn diagram of jujube action targets and disease targets; (b) jujube-component-target interactions network. Green hexagon is Chinese medicine, red circle is component, and blue diamond is target.
(a)(b)


## 3.3. Construction of the PPI Network and Screening of Key Targets

The screened potential targets were input into the STRING database to obtain the target protein interaction information and imported into Cytoscape to construct the PPI network (Figure3(a)), which had 130 nodes (target proteins) and 1004 edges (protein interactions). It indicated that among the predicted disease-related targets, the more targets could have effective interactions with that target. Using the five parameters of MNC, DEGREE, MCC, CLONESS, and EPC for screening (Figure 3(b)), the algorithm’s computational analysis of the network structure and weighted linkage between nodes could screen out important key genes. The intersection of the top 30 results of each algorithm was taken to obtain 11 key targets (Figure 3(c), Tables 2 and 3).Figure 3
Protein interaction network diagram. (a) All target protein interaction network; (b) top 30 target Venn diagram of MNC, DEGREE, MCC, CLONESS, and EPC; (c) key target protein interaction network diagram.
(a)(b)(c)Table 2
cytoHubba five-algorithm calculation results ranking.
RankMCCCLONESSMNCEPCDEGREE1MAOAMAPK3MAPK3MAPK3MAPK32DRD2EGFREGFRSRCEGFR3HTR1ASRCSRCEGFRSRC4DRD4HSP90AA1HSP90AA1APPHSP90AA15DRD3ESR1TNFHSP90AA1TNF6SLC6A3TNFESR1ESR1ESR17SLC18A2APPAPPPIK3CAAPP8ADRA2CPIK3CAPIK3CATNFPIK3CA9ADRA2ABRCA1BRCA1RELABRCA110HTR3ARELAGRIN2BBRCA1GRIN2B11HTR1BGSK3BSLC6A3GSK3BSLC6A312CHRNA4GRIN2BHTR3AGRIN2BHTR3A13ACHEMAPTRELANTRK2RELA14PIK3CAPRKCACHRNA4NFKB1CHRNA415HTR2APPARGACHECHRNA4ACHE16ADRA1BNTRK2CDK1DRD2GSK3B17ADRA1ANFKB1DRD2SLC6A3CDK118ADRA1DCDK1MAOAPRKCADRD219HTR2CDRD2NTRK2ACHEMAOA20HTR2BACHEGSK3BMAPTPPARG21APPHIF1APPARGIGF1RNTRK222OXTRCHRNA4HIF1ACDK1HIF1A23F2CDK5MAPTHIF1ASLC18A224NTSR1GRIN1MMP2HTR3AGRIN125SRCIGF1RTOP2ASLC18A2MAPT26EGFRABCB1NFKB1PPARGMMP227MAPK3PTPN1SLC18A2ADRA1BTOP2A28HSP90AA1MMP2GRIN1MAOANFKB129ESR1TOP2APRKCAMCL1PRKCA30RELAF2DRD4GRIN1DRD4Table 3
Analysis of topological parameters of key targets.
NameClosenessBetweennessDegreeACHE0.4867920.01486926APP0.5657890.06205839CHRNA40.4813430.01189227DRD20.4904940.01091126EGFR0.6113740.07756852ESR10.5707960.05457942HSP90AA10.5863640.07245646MAPK30.6417910.14001158PIK3CA0.560870.05258838RELA0.5330580.01169727SRC0.6056340.07343649


## 3.4. Subnetwork Analysis

MCODE subnetwork analysis can discover more closely connected groups or genes in the network. It is calculated by weighting the points with the highest weight and set as SEED, from SEED, recursively move outward to find nodes that can join the subnetwork. Subnetwork 1 is centered on MAOA, and the important targets connected to it, such as SLC6A3, DRD3, DRD4, SLC18A2, and HTR1A, are all dopaminergic synapse-related targets, indicating that subnetwork 1 is closely related to dopaminergic synapses (Figure4(a)). The core of subnetwork 2 is MMP2, and the important targets connected with it, such as TNF, are all targets related to inflammatory processes, indicating that subnetwork 2 is closely related to the regulation of tryptophan channels by inflammatory mediators (Figure 4(b)). The core of subnetwork 3 is HSPB1, and the important targets connected with it, such as CHRNA4, BRCA1, and HTR3A, are all targets related to cAMP signaling pathway, indicating that subnetwork 3 is closely related to cAMP signaling channels (Figure 4(c)).Figure 4
Subnetwork diagram. (a) Subnetwork 1; (b) subnetwork 2; (c) subnetwork 3. The larger the node, the darker the color means the higher the degree value of the target point.
(a)(b)(c)


## 3.5. GO Enrichment Analysis

In order to explore the functional distribution of key targets, 131 key targets were entered into the DAVID 6.8 database for GO enrichment analysis. The results showed 343 biological processes, 103 molecular functions, and 58 cellular components. Combined with the literature, the key targets in the biological process were filtered byP<0.5 and the number of enriched targets was high, and the key targets were concentrated in response to stimulus, signaling, cell proliferation, positive regulation of biological process, negative regulation of biological process, etc. (Figure 5(a) and Table 4). Among the molecular functions, molecular transducer activity, catalytic activity, transporter activity, and transcriptional regulation activity are mainly involved (Figure 5(b) and Table 4). Among the cellular components, the membrane is the most involved target, followed by synapse, organelle part, and cell junctions (Figure 5(c) and Table 4).Figure 5
GO enrichment analysis of the key targets. (a) Biological process (BP) analysis result; (b) molecular functions (MF) analysis result; (c) cellular components (CC) analysis result.
(a)(b)(c)Table 4
GO analysis table.
ClassGOTermCountP valueMolecular functionGO:0060089Molecular transducer activity542.19E − 24GO:0003824Catalytic activity744.44E − 09GO:0005215Transporter activity219.23E − 05GO:0005488Binding1283.18E − 03GO:0104005Hijacked molecular function31.70E − 02GO:0016209Antioxidant activity32.42E − 02GO:0098772Molecular function regulator204.35E − 02GO:0140110Transcription regulator activity172.23E − 01GO:0005198Structural molecule activity47.77E − 01Cellular componentGO:0045202Synapse495.58E − 23GO:0044456Synapse part343.96E − 15GO:0030054Cell junction302.44E − 08GO:0044425Membrane part814.76E − 08GO:0032991Protein-containing complex682.16E − 07GO:0016020Membrane966.75E − 07GO:0031974Membrane-enclosed lumen661.02E − 06GO:0044422Organelle part962.78E − 05GO:0005623Cell1301.52E − 04GO:0044464Cell part1301.52E − 04Biological processGO:0050896Response to stimulus1253.86E − 29GO:0023052Signaling1081.98E − 26GO:0032501Multicellular organismal process1104.81E − 22GO:0048518Positive regulation of biological process995.94E − 21GO:0007610Behavior332.62E − 20GO:0008283Cell proliferation531.19E − 18GO:0048511Rhythmic process235.07E − 17GO:0048519Negative regulation of biological process885.95E − 17GO:0032502Developmental process914.46E − 15GO:0065007Biological regulation1267.53E − 15


## 3.6. KEGG Pathway Analysis

The KEGG pathway enrichment analysis of potential targets by the DAVID 6.8 data platform (P<0.05) is shown in Figures 6(a)–6(d). The top 10 pathways were Neuroactive ligand-receptor interaction, Serotonergic synapse , Gap junction, cAMP signaling pathway, Dopaminergic synapse, Calcium signaling pathway, Hypoxia-inducible factor pathway (HIF-1 signaling pathway), Prolactin signaling pathway, Thyroid hormone signaling pathway, and Inflammatory mediator regulation of TRP channels (Table 5). It is suggested that jujube components may exert appetite modulating effects through the above pathways.Figure 6
KEGG pathway enrichment results. (a) KEGG pathway network diagram; (b) KEGG enrichment pathway annotated classification results; (c) KEGG pathway enrichment circle diagram; (d) KEGG enrichment result bar graph.
(a)(b)(c)(d)Table 5
KEGG pathway enrichment information.
PathwayCountP valueNeuroactive ligand-receptor interaction301.28E − 15Serotonergic synapse184.07E − 14Gap junction153.49E − 12cAMP signaling pathway193.98E − 10Dopaminergic synapse141.04E − 08Calcium signaling pathway186.88E − 08HIF-1 signaling pathway109.55E − 06Prolactin signaling pathway81.60E − 05Thyroid hormone signaling pathway101.68E − 05Inflammatory mediator regulation of TRP channels91.73E − 05


## 3.7. “Active Ingredient-Key Target-Pathway” Network Construction for Jujube

The active ingredients, potential targets, and selected signaling pathways of jujube were imported into Cytoscape 3.8.2 software to construct the “jujube-active-ingredient-target-action pathway” diagram (Figure7(a)). The MCC algorithm of cytoHubba was used to further calculate the closest association of each component with key targets (Figure 7(b)). The results concluded that 5-hydroxytryptamine (5-HT)-containing synapses was the key pathway of action.Figure 7
Dates-component-target-pathway diagram. (a) Dazao-ingredient-target-pathway diagram; (b) the key target-ingredient diagram. The green hexagon is the name of a single herbal medicine, the red circle is the ingredient, the blue diamond is the target, and the purple arrow is the pathway name.
(a)(b)


## 3.8. Molecular Docking Simulation of Target Interaction with Related Compounds

The docking of key target genes and their related compounds revealed that the docking binding energies of the targets and their related compounds were all negative (Table6), indicating that the related compounds bind well to the targets. Among them, coclaurine had the lowest docking binding energy of −8.38 with ESR1, (−)-Catechin and ESR1, (+)-Stepholidine and ADRA2C, Berberine and ADRA2C, Cianidanol and ESR1, Coclaurine and ESR1, Moupinamide and MMP2, Quercetin and SCR, and Quercetin and TNF all had binding energies <−5 kcal/mol, and the binding patterns are shown in Figures 8(a)–8(h).Table 6
Molecular docking results of key targets and their related compounds.
CompoundTargetCombined energy (kcal/mol)(−)-CatechinESR1−7.01ADRA2C−5.47(+)-StepholidineADRA2C−7.10BerberineADRA2C−7.80CianidanolESR1−7.14ADRA2C−5.40CoclaurineESR1−7.02ADRA2C−4.82CoumestrolESR1−8.38MoupinamideMMP2−7.87EGFR−4.03NuciferineADRA2C−6.79QuercetinTNF−7.61SRC−6.95MMP2−6.72MAPK3−5.24HSP90AA1−6.25ESR1−7.36EGFR−4.83APP−4.57Figure 8
Molecular docking simulations of targets and compounds. (a) 3D and 2D plots of molecular docking of (−)-catechin-ESR1; (b) 3D and 2D plots of molecular docking of (+)-stepholidine-ADRA2C; (c) molecular docking 3D and 2D diagrams of berberine-ADRA2C; (d) molecular docking 3D and 2D diagrams of cianidanol-ESR1; (e) molecular docking 3D and 2D plots of coclaurine-ESR1; (f) molecular docking 3D and 2D plots of moupinamide-MMP2 target; (g) molecular docking 3D and 2D diagrams of quercetin-SCR; (h) molecular docking 3D and 2D diagrams of quercetin-TNF.
(a)(b)(c)(d)(e)(f)(g)(h)


## 4. Discussion

Chinese medicine is difficult to elucidate molecular mechanisms due to the complexity of their chemical composition. Network pharmacology has been increasingly applied to the study of Chinese pharmaceutical preparations in recent years [11]. It transforms drug research from a “single target, single drug” model to a “network target, multicomponent therapy” model [35]. Studies have shown a favorable effect of jujube on glycosylated hemoglobin and some antioxidant effects in patients with T2DM [36]. However, for the time being, no studies have been conducted to analyze the network pharmacology of the active ingredients of jujube. Therefore, in this study, based on network pharmacology and molecular docking research methods, we constructed a multidimensional network through target prediction and protein interaction networks to elucidate the principle of action of jujube in treating diabetic hunger and regulating appetite from molecular prediction level.In the present study, a total of 16 components including coclaurine, (−)-catenin, (+)-stepholidine, berberine, cianidanol, coclaurine, and moupinamide were identified as potential active ingredients of jujube. These active ingredients include 131 targets of action related to appetite and appetite regulation. Through PPI network analysis of jujube in appetite regulation, we identified 11 key targets: MAPK3, EGFR, SRC, HSP90AA1, and so on. For further screening by MCODE analysis, three key targets (MAOA, MMP2, and HSPB1) were identified. The molecular docking results showed that the main components of jujube had strong binding activity to the main targets (ESR1, ADRA2C, and MMP2).MAOA can metabolize monoamine neurotransmitters [37]. Studies have shown that MAOA regulates food intake and energy expenditure [38]. Gardner et al. [39] similarly showed that MAOA is involved in regulating appetite and food intake related to obesity genes. HSP27 regulates actin dynamics and thus cell motility [40]. The inhibition of feeding by fibroblast growth factor (FGF)-1 is accompanied by the induction of HSP27 in periventricular astrocytes [41]. MMP2 is an endopeptidase that reduces the basement membrane around adipocytes, thus promoting the development of adipocyte hypertrophy [42]. Studies have shown high levels of MMP2 gene expression in patients with T2DM [43]. Previous studies suggested that the significant anorexigenic effect of estradiol in male rats would be related to ESR1 present in the lateral hypothalamic region [44]. The results of molecular docking suggested that coclaurine had the lowest docking binding energy with ESR1 at -8.38. And the stronger binding activity of coumestrol to ESR1 in jujube may be related to the estrogen-like effect of coumestrol. The accuracy of network prediction was reflected from the side.Further, GO enrichment and KEGG pathway enrichment analysis showed that key genes act on signaling pathways such as neuroactive ligand-receptor interaction, serotonergic (5-hydroxytryptamine; 5-HT) synapse, cAMP signaling pathway, dopaminergic synapse, calcium signaling pathway, and hypoxia-inducible factor pathway. Among them, the key pathway of 5-HT synapses was closely related to appetite control and the treatment of related diseases. The key role of 5-HT in appetite control was formally proposed almost 30 years ago [45]. Studies have shown that the biogenic amine neurotransmitter 5-HT is negatively correlated with food intake and that a decrease in food intake is associated with 5-HT [46]. Agonists of 5-HT improve obesity and glycemic control in the population [47].However, it should be noted that network pharmacology is based on existing databases and results for network modeling, and there is a certain false-positive rate of predicted results due to the differences in raw experimental data under different experimental conditions.This paper presents a predictive analysis of the appetite regulation mechanism of jujube based on the theoretical level, and we hope that the results of this analysis can provide new ideas for the next in-depth research, and we also expect better research basis at the level of new drug development. Our research team will continue to focus on the progress of pharmacological research on jujube and will use the results of this paper as a reference to conduct relevant animal and cellular experiments to further investigate the effects and regulatory mechanisms of jujube on the appetite center of the hypothalamus.



---

*Source: 5070086-2022-04-10.xml*