NP analysis
Bioactive Components and Associated Targets in YHS
TCMSP was used to search and collect 77 components in total. A total of 48 components were identified based on YSH components (OB ≥ 30%, DL ≥ 0.18) (Table 1). There were 183 targets estimated in YHS via the SymMap database.
Table 1
Effective components of Corydalis yanhusuo.
Molecule ID | Molecule name | Molecular weight | OB (%) | DL |
MOL001454 | berberine | 336.39 | 36.86 | 0.78 |
MOL001458 | coptisine | 320.34 | 30.67 | 0.86 |
MOL001460 | Cryptopin | 369.45 | 78.74 | 0.72 |
MOL001461 | Dihydrochelerythrine | 349.41 | 32.73 | 0.81 |
MOL001463 | Dihydrosanguinarine | 333.36 | 59.31 | 0.86 |
MOL001474 | sanguinarine | 332.35 | 37.81 | 0.86 |
MOL000217 | (S)-Scoulerine | 327.41 | 32.28 | 0.54 |
MOL002670 | Cavidine | 353.45 | 35.64 | 0.81 |
MOL002903 | (R)-Canadine | 339.42 | 55.37 | 0.77 |
MOL000359 | sitosterol | 414.79 | 36.91 | 0.75 |
MOL004071 | Hyndarin | 355.47 | 73.94 | 0.64 |
MOL004190 | (-)-alpha-N-methylcanadine | 354.46 | 45.06 | 0.8 |
MOL004191 | Capaurine | 371.47 | 62.91 | 0.69 |
MOL004193 | Clarkeanidine | 327.41 | 86.65 | 0.54 |
MOL004195 | CORYDALINE | 369.5 | 65.84 | 0.68 |
MOL004196 | Corydalmine | 340.45 | 52.5 | 0.59 |
MOL004197 | Corydine | 341.44 | 37.16 | 0.55 |
MOL004198 | 18797–79–0 | 367.43 | 46.06 | 0.85 |
MOL004199 | Corynoloxine | 365.41 | 38.12 | 0.6 |
MOL004200 | methyl-[2-(3,4,6,7-tetramethoxy-1-phenanthryl)ethyl]amine | 355.47 | 61.15 | 0.44 |
MOL004202 | dehydrocavidine | 351.43 | 38.99 | 0.81 |
MOL004203 | Dehydrocorybulbine | 352.44 | 46.97 | 0.63 |
MOL004204 | dehydrocorydaline | 366.47 | 41.98 | 0.68 |
MOL004205 | Dehydrocorydalmine | 338.41 | 43.9 | 0.59 |
MOL004208 | demethylcorydalmatine | 327.41 | 38.99 | 0.54 |
MOL004209 | 13-methyldehydrocorydalmine | 352.44 | 35.94 | 0.63 |
MOL004210 | (1S,8'R)-6,7-dimethoxy-2-methylspiro [3,4-dihydroisoquinoline-1,7'-6,8-dihydrocyclopenta [g] [1, 3]benzodioxole]-8'-ol | 369.45 | 43.95 | 0.72 |
MOL004214 | isocorybulbine | 368.51 | 40.18 | 0.66 |
MOL004215 | leonticine | 327.46 | 45.79 | 0.26 |
MOL004216 | 13-methylpalmatrubine | 352.44 | 40.97 | 0.63 |
MOL004220 | N-methyllaurotetanine | 341.44 | 41.62 | 0.56 |
MOL004221 | norglaucing | 341.44 | 30.35 | 0.56 |
MOL004224 | pontevedrine | 381.41 | 30.28 | 0.71 |
MOL004225 | pseudocoptisine | 320.34 | 38.97 | 0.86 |
MOL004226 | 24240–05–9 | 353.4 | 53.75 | 0.83 |
MOL004228 | saulatine | 396.47 | 42.74 | 0.79 |
MOL004230 | stylopine | 323.37 | 48.25 | 0.85 |
MOL004231 | Tetrahydrocorysamine | 337.4 | 34.17 | 0.86 |
MOL004232 | tetrahydroprotopapaverine | 329.43 | 57.28 | 0.33 |
MOL004233 | ST057701 | 341.44 | 31.87 | 0.56 |
MOL004234 | 2,3,9,10-tetramethoxy-13-methyl-5,6-dihydroisoquinolino [2,1-b]isoquinolin-8-one | 381.46 | 76.77 | 0.73 |
MOL000449 | Stigmasterol | 412.77 | 43.83 | 0.76 |
MOL000785 | palmatine | 352.44 | 64.6 | 0.65 |
MOL000787 | Fumarine | 353.4 | 59.26 | 0.83 |
MOL000790 | Isocorypalmine | 341.44 | 35.77 | 0.59 |
MOL000791 | bicuculline | 367.38 | 69.67 | 0.88 |
MOL000793 | C09367 | 325.39 | 47.54 | 0.69 |
MOL000098 | quercetin | 302.25 | 46.43 | 0.28 |
HCC Targets Prediction as well as Shared Targets
Based on the GeneCards database, 17613 HCC-related targets were identified. We showed 88 genes to be potential therapeutic targets using scores ≥ 55 as the screening condition. The YHS targets were compared with HCC targets, and the Venn diagram shows that there were 12 shared targets (Fig. 2A). The obtained shared targets were YHS’s potential anti-HCC targets.
PPI Network Construction and Drug-Components-Targets
These 12 shared targets were entered into the String database to create the target PPI prediction network (Fig. 2B). With the elimination of isolated target proteins for obtaining protein interaction data, the confidence level was > 0.4. The drug-ingredient-target network was created by Cytoscape to obtain relationships of effective components with shared targets and HCC (Fig. 3). The drug-ingredient-target network had 231 nodes and 702 edges. Quercetin (MOL000098, degree = 150), Isocorypalmine (MOL000790, degree = 21), Hyndarin (MOL004071 degree = 20), (S)-Scoulerine (MOL000217, degree = 19), leonticine (MOL004215, degree = 18), and (R)-Canadine (MOL002903, degree = 17) were identified as an important active compound. Important anti-HCC therapeutic targets identified were PTGS2 (degree = 47), PTGS1 (degree = 43), SCN5A (degree = 42), KCNH2 (degree = 39), CALM3 (degree = 36), RXRA (degree = 31), ADRB2 (degree = 29), OPRD1 (degree = 28), and OPRM1 (degree = 28). Based on the above findings, the targets were closely linked to others in the network, which could have significant implications for HCC.
GO and KEGG Analysis
We used the DAVID database to analyze GO and KEGG enrichment on 12 shared targets. Altogether 118 biological processes (BP), 10 cellular components (CC), and 19 molecular functions (MF) terms were enriched, and they met the EASE scores ≤ 0.05 and Count ≥ 2. Figure 4A-4C shows 10 significant GO terms. Typically, the above genes were enriched in BP terms such as cellular response to organic cyclic compounds, positive regulation of gene expression and transcription using DNA as a template, protein complexes, mitochondria, cytosol, transcription factor binding, and protein binding. Therefore, YHS can be used to treat HCC by regulating cancer cell proliferation through various gene biological activities. KEGG enrichment on the target was done to identify the underlying pathways related to the treatment of HCC by YHS. Figure 4D shows the top ten pathways enriched using KEGG. The enriched pathways were closely related to YSH, like, for example, TNF and MAPK pathways.
Target-Pathway Network and Analysis
To characterize the YSH-related mechanism in the treatment of HCC, we built a target-pathway network with 11 proteins (NFE2L2 removed) and the associated pathways (Fig. 5). It covered 21 nodes and 70 edges (11 and 10 for proteins and pathways separately). The Hepatitis B pathway, cancer pathways, TNF pathway, and MAPK pathway were the most significant, with the highest degree values (Table 2). AKT1, MAPK1, CASP3, and TNF had high levels and had important effects on cell growth, making them a critical therapeutic target in YHS for treating HCC.
Table 2
Topological structure analysis of the targets-pathways network of Corydalis yanhusuo targeting HCC.
Node | Degree | Betweenness | Closeness |
Hepatitis B | 9 | 0.071848 | 0.606061 |
AKT1 | 10 | 0.098886 | 0.666667 |
MAPK1 | 10 | 0.098886 | 0.666667 |
CASP3 | 8 | 0.055394 | 0.588235 |
IL6 | 5 | 0.015682 | 0.47619 |
TNF | 7 | 0.039278 | 0.526316 |
CASP8 | 6 | 0.023456 | 0.5 |
TP53 | 6 | 0.029093 | 0.526316 |
MYC | 6 | 0.029093 | 0.526316 |
TGFB1 | 9 | 0.077105 | 0.625 |
Proteoglycans in cancer | 8 | 0.131336 | 0.571429 |
IGF2 | 1 | 0 | 0.37037 |
Pathways in cancer | 9 | 0.116906 | 0.606061 |
CDKN2A | 2 | 0.001548 | 0.4 |
Colorectal cancer | 6 | 0.01908 | 0.512821 |
Chronic myeloid leukemia | 6 | 0.04642 | 0.512821 |
Tuberculosis | 7 | 0.033935 | 0.540541 |
Chagas disease (American trypanosomiasis) | 6 | 0.02507 | 0.512821 |
TNF signaling pathway | 6 | 0.024763 | 0.512821 |
Toxoplasmosis | 6 | 0.020358 | 0.512821 |
MAPK signaling pathway | 7 | 0.031336 | 0.540541 |
Molecular Docking
The network pharmacology analysis revealed that quercetin and hyndarin interact with TNF-α and P38-MAPK at different levels (Fig. 6). The lower vina scores indicated a stable and potent compound-protein interaction. Vina scores for quercetin and hyndarin were less than − 0.6, indicating that quercetin and hyndarin had a stronger and more stable binding affinity for TNF-α and P38-MAPK. Consequently, quercetin and hyndarin were the suitable interacting partners for TNF-α and P38-MAPK protein.
Experimental Validation
Using network pharmacology analysis, quercetin and hyndarin were identified as important components of Corydalis yanhusuo in the treatment of HCC. Based on NP analysis, quercetin, and hyndarin were obtained. These two molecules were then validated for their role in resisting cancer cell proliferation.
For validation, HepG2 cells were treated for 48 h with quercetin and hyndarin. Following treatment, the IC50 values were found to be 178.8 µg/mL, and 413.4 µM, respectively. Quercetin and hyndarin also showed dose-dependent growth inhibition of HepG2 cells (Fig. 7A). Additionally, with the increase in doses of quercetin and hyndarin, the HepG2 clone number declined (p < 0.001) (Fig. 7B-C). Similarly, the HepG2 cell apoptosis degree was found to increase with an increase in quercetin and hyndarin concentration. Upon statistical comparison, the treatment group had increased apoptosis levels compared with the control group (p < 0.001) (Fig. 8). In vitro experiments yielded similar NP analysis results.
After 48-h of quercetin and hyndarin treatment of HepG2 cells, MAPK protein expression was down-regulated (p < 0.01), but MAPK and TNF protein levels within the HepG2 cells of the quercetin group were not statistically significant (Fig. 9), indicating that hyndarin has inhibitory effects on MAPK signaling pathways, and quercetin may have inhibitory effects on other pathways.