Research on the Potential Mechanism of Gentiopicroside Against Gastric Cancer Based on Network Pharmacology.

BACKGROUND: Gastric cancer was still one of the commonly diagnosed cancer types and the third-most common cause of cancer-related death in the world. Gentiopicroside, which is extracted from the Gentianella acuta, is commonly used in both traditional treatment and modern clinical care; therefore, its anticancer effects have been attracted more attention. However, the systematic analysis of action mechanism of Gentiopicroside on gastric cancer (GC) has not yet been carried out. AIM: A network pharmacology-based strategy combined with molecular docking studies and in vitro validation was employed to investigate potential targets and molecular mechanism of Gentiopicroside against GC.
MATERIALS AND METHODS: Potential targets of Gentiopicroside, as well as related genes of GC, were acquired from public databases. Potential targets, and signaling pathways were determined through bioinformatic analysis, including protein-protein interaction (PPI), the Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Subsequently, molecular docking and cell experiments were performed to further verify the above findings.
RESULTS: Our findings revealed that the anticancer activity of Gentiopicroside potentially involves 53 putative identified target genes. In addition, GO, KEGG, and network analyses revealed that these targets were associated with cell proliferation, metabolic process, and other physiological processes. Furthermore, we have proved that critical compound affected the expression of CCND1, CCNE1, p-AKT and p-P38 at protein levels. These findings provide an overview of the anticancer action of Gentiopicroside from a network perspective; meanwhile, it might also set an example for future studies of other materials used in traditional Chinese medicine (TCM).
CONCLUSION: This study comprehensively illuminated the potential targets and molecular mechanism of Gentiopicroside against GC. It also provided a promising approach to uncover the scientific basis and therapeutic mechanism of TCM treating for disease.

Introduction

According to the estimates from the International Agency for Research on Cancer (IARC) GLOBOCAN project, gastric cancer was still the fifth-most commonly diagnosed cancer type in 2018 and the third-most common cause of cancer-related death.1 And it is a multifactorial disease concerning Helicobacter pylori infection, lifestyle, genetic, or other risk factors.2 Which has posed a formidable potential threat to public health owing to its high morbidity and mortality.3 As we all know, the therapies for GC have been well developed in recent decades, the treatment strategies for GC mainly include surgery, chemotherapy, radiotherapy, targeted therapies, and immunotherapy; nevertheless, its mortality remains relatively high as a result of frequent recurrence and metastasis.4 Traditional Chinese medicine mainly from natural products is one of the most abundant resources of active compounds for drug discovery, playing an important role in the history of disease treatments, such as artemisinin used to cure malaria,5 arsenic trioxide in the treatment of acute promyelocytic leukaemia,6,7 and berberine in the treatment of type 2 diabetes mellitus,8,9 et al. Besides, TCM has been widely used in anticancer therapeutic prescription for thousands of years in China, basing on multiple mechanisms of anticancer effects, including apoptosis induction, proliferation inhibition, metastasis suppression, multidrug resistance reversal and immune function regulation.10–12 The advantage of using TCM is related to its ability on affecting multiple signaling pathways and molecular targets cancer-related, while causing few side effects.13 Gentiopicroside isolated from Gentianella acuta, which has been shown to exert a wide spectrum of pharmacological activities, such as ameliorates diabetic nephropathy,14 reducing Pain, Pruritus, and Corticosteroid Addictive Dermatitis,15 as well as antidepressant-like effects,16 and so on. Close attention has been paid to Gentiopicroside due to its potential role in the treatment of cancer.17–19 Accordingly, the molecular mechanisms which Gentiopicroside induces and the resulting changes in cellular phenotypes are widely studied. Meanwhile, employment of computational methodologies for the identification of drug targets and the underlying mechanisms is becoming mainstream for further study. In particular, computational target identification and following the molecular mechanisms can accelerate the drug discovery and drug design process.20–23

Therefore, we elucidated the pharmacological actions of Gentiopicroside on gastric cancer systematically using computational methodologies. Then, we used the molecular docking and a series of in vitro assays to validate its effects on GC. A schematic diagram of the analysis procedures for Gentiopicroside target gene prediction is shown in Figure 1.

Figure 1

Workflow for the identification potential mechanism of Gentiopicroside in treating gastric cancer.

Materials and Methods

Computational Targets Fishing

Structural information of Gentiopicroside (Figure 2A, CAS:20831–76-9) was obtained from the NCBI PubChem (). At the same time, the protein targets of this compound were retrieved both from the TCMSP () and Pharmmapper (). Data on the Gastric cancer-related targets were obtained from two databases: (i) We obtained 137 Gastric cancer-related targets from the Online Mendelian Inheritance in Man24 (), an information platform mainly designed for human genes and genetic diseases. (ii) We obtained 18,223 Gastric cancer-related targets from the GeneCards: The Human Gene Database25 (). Common targets for Gentiopicroside and GC were obtained using Venny 2.1 (). Common targets for Gentiopicroside and GC were obtained using Venny 2.1 ().

Figure 2

The venn diagram of the targets both in gastric cancer targets and Gentiopicroside targets.

Protein–Protein Interaction (PPI) Network Construction

The common targets related to Gentiopicroside treatment of GC were imported into the STRING database () to obtain their interaction relationship, “Homo sapiens” was chosen and a scoring value > 0.9 was selected as the high confidence basis for protein interactions. Then, the common targets PPI network files were imported into Cytoscape 3.6.0 to merge and got the intersection of PPI network and core genes.

Gene Function and Pathway Enrichment Analysis

Gene ontology (Go) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by WebGestalt ().

Network Construction

In order to understand the complex relationships among compound, targets, and diseases, we used Cytoscape26 () to construct and analyze the networks. It provides a powerful set of data integration, analysis, and visualization functions to analyze the complicated networks.

Molecular Docking Analysis

Autodock vina 1.1.2 was used for network pharmacology-based prediction and analysis and illustrated the mechanism of ligand acting on a complex molecular network by applying high-precision docking simulation and molecular pathway map. We used Autodock vina to evaluate the binding potential between selected target and thecompound Gentiopicroside.

Cell Culture

The human gastric cancer cell lines HGC-27, MGC-803, BGC-823, AGS, and GES-1 were purchased from Chinese Academy of Sciences Shanghai Cell Bank and were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum (Invitrogen, Carlsbad, CA, USA).

Cell Proliferation Assay

The cell viability of GC cells was measured by MTT assay, cells were plated into 96-well plates at a density of 2000 cells/well. After an overnight culture, cells were treated in DMEM with 0, 25, 50, 100, 200, 400µM of Gentiopicroside for 48hrs. Then, the cells were incubated with 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT) for 4 h at 37°C. The absorbance at 490 nm was measured.

Flow Cytometry

Cell cycle analysis was performed 48 h after transfection with propidium iodide staining as described previously. Cell apoptosis was also analyzed 48 h by Annexin V and propidium iodide staining as described previously. All the samples were assayed in triplicate.

Western Blot Analysis

Protein were separated by 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE), and transferred to 0.22 μm nitrocellulose membranes (Sigma) and incubated with specific primary antihuman antibodies. They were then incubated with specific antibodies according to the manufacturer’s protocol. The GAPDH antibody was used as the control.

Statistical Analysis

Comparisons between groups were analyzed using the Student’s t-test and one-way analysis of variance (ANOVA) by GraphPad Prism 7.0 (GraphPad Software Inc., California, USA). P < 0.05 was considered to indicate a statistically significant difference.

Results

Targets Fishing

There are a total of 148 drug-related targets were retrieved. Afterwards, from OMIM and GeneCards database, there were 18,324 gastric caner-related target genes. Then, we inputted drug targets and disease targets into Venny 2.1, after eliminating the redundancy we found 135 common targets between Gentiopicroside and GC, the common targets are listed in Figure 2B.

Target Protein Interaction Network Analysis

The common targets were imported into the STRING database to obtain their interaction relationship, and a scoring value > 0.9 was selected as the high confidence basis for protein interactions. The interaction network has 135 nodes and 64 edges (Figure 3A). TSV data were imported into Cytoscape, and the results showed a total of 61 nodes and 64 edges. Then, 10 targets were further screened as the center genes, including CAT, GOT1, LDHB, OCRL, CBL, CCND1, RAC3, GPI, B2M, TALDO1 (Figure 3B).

Figure 3

Construction of the protein–protein interaction network of Gentiopicroside and gastric cancer’ associated targets. (A) Clusters of screened PPI networks. (B) The top 10 potential effective targets in gastric cancer.

GO and Pathway Analysis

In order to obtain a more detailed understanding of the potential mechanisms of Gentiopicroside acting on gastric cancer, GO function and Pathway analysis were performed using WebGestalt. As shown in Figure 4A, the top seven functions were used as metabolic process, biological regulation, response to stimulus, cellular component, cell proliferation, reproduction, growth. These functional terms are highly relevant to anti-tumor activities. As for pathway analysis, the 53 targets participate in 10 KEGG pathways with significant P-value < 0.01 including Metabolic pathways, PI3K-Akt signaling pathway and so on, which are shown in Figure 4B.

Figure 4

GO analysis and KEGG analysis for the major targets of Gentiopicroside. (A) GO map of putative target genes. (B) KEGG pathway analysis of putative target genes.

Network Analysis

Based on target and pathway analyses, the networks of the active compound and associated targets showed the relationships between compound and targets. As shown in Figure 5, this compound, disease, targets and pathways interaction network has 65 nodes and 168 edges. The red, blue, pink and green nodes represent the disease, Pathways, targets and compound; respectively, the edges represent the interactions among them.

Figure 5

Drug-GC-targets-Pathways PPI network (The red, blue, pink and green nodes represent the disease, Pathways, targets and compound; respectively, the edges represent the interactions among them).

Docking Results Analysis

We selected CCND1 as the target gene for molecular docking. This gene was mapped to the compound for molecular docking. Affinity was the score for the molecular docking, and when the score was lower, the binding affinity was stronger. An affinity <−7 indicated a stronger binding activity.27 The vina score of Gentiopicroside for treating gastric cancer was −5.4 kcal mol−1. In Figure 6, we have shown the combination of the compound and gene in a 3D graph. The results indicated that Gentiopicroside had a good binding activity to CCND1.

Figure 6

Molecular docking. (A) Stereogram of the interaction between CCND1 and Gentiopicroside. (B) Plane chart of the interaction between CCND1 and Gentiopicroside.

The Effects of  Gentiopicroside  on the Regulation of GC

We assessed the effects of Gentiopicroside on GC cell viability. Four human GC cell lines and GES-1 were exposed to different concentrations of Gentiopicroside for 48hrs. As shown in Figure 7A, Gentiopicroside exerted obvious growth inhibition of GC cells with IC50 values of about 62–205µM for 48hrs, of which HGC-27 was the most sensitive to the cytotoxic effect of Gentiopicroside (IC50 value: 62.03±7.63µM), whereas GES-1 was the least sensitive. The results indicate that at a certain dosage range, Gentiopicroside might selectively kill human GC cells rather than normal human gastric epithelial cells. And in the HGC-27 cell, the vitality of different concentration of Gentiopicroside was significantly inhibited compared with control groups (Figure 7B). Thus, we selected HGC-27 cell in the next assays.

Figure 7

(A) Gentiopicroside displayed a significant growth inhibition in GC cells with IC50 values. (B) Gentiopicroside suppressed HGC-27 cells proliferation significantly under the concentration of 0, 25, 50,100, 200, 400 μM with 24 h, 48 h and 72h incubation, respectively. (*P<0.05, **P<0.01).

Gentiopicroside Induces Apoptosis and Cell Cycle Arrest

Flow cytometry analysis was performed to evaluate the effects of Gentiopicroside on apoptosis and cell cycle distribution of HGC-27 cells. The results indicated that the proportion of these cells in apoptosis was significantly increased after treating by Gentiopicroside (Figure 8A). Cell cycle distribution showed that the proportion of HGC-27 cells was increased in the G0/G1 phase, but decreased in the S phase in the Gentiopicroside-treated group (Figure 8B).

Figure 8

(A) Gentiopicroside suppressed HGC-27 cell apoptosis significantly under the concentration of 100 μM with 48 h incubation, respectively. (B) Gentiopicroside suppressed HGC-27 cell cycle significantly under the concentration of 100 μM with 48 h incubation, respectively. (*P<0.05, **P<0.01).

Effects of Gentiopicroside on Potential Targets in vitro

In this study, we selected CCND1, CCNE1, AKT, P38 as the potential targets of Gentiopicroside combined with the basis of network pharmacology. Western blotting analysis was used to evaluate these protein levels in vitro. As shown in Figure 9, p-P38 protein levels were significantly increased in the Gentiopicroside treated group compared with the control group. However, the protein expression of p-AKT, CCND1, CCNE1 was significantly attenuated after Gentiopicroside treatment (Figure 9A). What is more, after using LY294002 to inhibit the PI3K/AKT pathway, the level of p-AKT inhibited further, and it also has a significant effect on p-P38 (Figure 9B).

Figure 9

Effects of Gentiopicroside and LY294002 on the expression of PI3K/AKT pathway in HGC-27 cells. (A) WB confirmed that Gentiopicroside reduced CCND1, CCNE1 and p-AKT protein expression. While p-P38 significantly increased after stimulation with Gentiopicroside for 48 h. (B) HGC-27 cells were treated with Gentiopicroside and LY294002 as indicated, AKT, P38, p-P38 and p-AKT protein levels were determined using Western blotting.

Discussion

Traditional Chinese medicine (TCM) has more than 2000 years of history and has gained widespread clinical applications. However, the explicit role of TCM in preventing and treating GC remains unclear due to lack of sound scientific evidence. Mostly it is very hard to clarify the effective mechanism of TCM because of its multiple ingredients and targets. In recent years, virtual pharmacology research was regarded as the fastest and most effective screening method in the early study of drug effectiveness.21,23,28–32 Gentianella acuta is an old drug commonly used for treating hepatitis and carcinoma. Therefore, in this study, Gentiopicroside was employed to study the gastric cancer enriching mechanism by network pharmacology approach. We integrated information from publicly available databases to predict interactions between Gentiopicroside and its potential targets in GC, as well as the signaling pathways and networks involved. In this study, 53 targets were screened and participated in 10 pathways based on P-values of each enriched pathway and their relationship to GC, KEGG analysis of the core PPI network indicated that Gentiopicroside may result in multiple effects for gastric cancer treatment by regulating multiple pathways. Furthermore, CCND1 with higher score, which was closely related to the occurrence and development of GC, was selected for molecular docking simulation by the LeDock. And the possible mechanism of Gentiopicroside inhibited tumor development through regulating PI3K-Akt signaling pathway. Accumulating evidence indicates that the PI3K/AKT signaling pathway, which is one of the most frequently activated signal-transduction pathways in cancer, is involved in cell cycle, cell proliferation, growth, migration, angiogenesis, and apoptosis. These results are consistent with several studies in vitro suggesting that Gentiopicroside exhibited anticancer effects on ovarian cancer cell lines through induction of apoptosis,17 while inhibited Cell Growth and Migration on Cervical Cancer via the Reciprocal MAPK/Akt Signaling Pathways.33 And to the best of our knowledge, this is the first study to report that Gentiopicroside exhibited anticancer effect on gastric cancer.

This study reports multicomponent therapeutics by Gentiopicroside to provide multitarget, multipathway regulation of GC activity. Although all of the findings in this paper gave evidence to prove the major active compound in Gentianella acuta to gastric cancer, the limitation of our study is that the molecular mechanism of this study was too simple, only a series of simple assays in vitro were carried. In future, more experiments are needed to provide support for our findings, and future studies should focus on providing experimental evidence and expanding the role of Gentianella acuta on gastric cancer.