Prediction of the Active Ingredients and Mechanism of ASH against Liver Cancer Based on Network Pharmacology and Molecular Docking ()
1. Introduction
Artemisiae scopariae Herba (ASH) also called Artemisia capillaris (AC), belongs to the family of Asteraceae and the genus Artemisia [1] . As a therapeutic traditional medicine, it showed the anti-inflammatory effects in chronic hepatitis B virus infection and liver cirrhosis [2] . Moreover, the major constituents of ASH such as capillin and scoparone exhibit anti-cancer effects in liver, prostate, and lung cancers [2] . The extract of ASH (AC68) not only induced apoptosis but also inhibited cell growth, migration, and invasion of liver cancer cells by blocking the PI3K/AKT pathway [3] . ASH has certain effects on the treatment of hepatitis, also inhibits cell growth, invasion, and metastasis in liver cancer, together with induction of apoptosis, as well as hepatoprotective property [3] [4] .
Liver cancer (LC) is the second most common cause of cancer-related death worldwide, it ranks fifth in terms of global cases and second in terms of deaths for males, hepatocellular carcinoma (HCC) is the most common type of LC worldwide [5] [6] . Lack of suitable biomarkers for early detection and limited treatment strategies are the major causes of high mortality [7] . A variety of risk factors have been associated with the development of LC, including hepatitis viruses, cirrhosis obesity and fatty liver disease [8] . The tumor microenvironment (TME) plays an important role in tumor progression and metastasis which contributes to tumor cell proliferation, survival, migration, and invasion [9] . More and more studies have revealed that TME has critical roles in the progression of LC [10] .
There are several treatment options for LC including chemotherapy, surgery, radiation and immunotherapy [11] . Unfortunately, each of the treatment options suffers some drawbacks. Chinese medicines (CMs) have potential to both prevent LC occurrence and retard LC progression. The actions of CMs on LC may include tumor growth inhibition, antimetastatic activities, anti-inflammation, anti-LC stem cells, reversal on multi-drug resistance and induction/reduction of oxidative stress [12] . As a traditional Chinese medicine, ASH shows hepatoprotective property, and there is no effective therapy available to treat LC at present, which is expected to have a certain effect in the treatment of LC.
2. Methods
2.1. The Putative Targets of ASH
“Artemisiae scopariae Herba” was used as a keyword in the traditional Chinese medicine system pharmacology database and analysis platform (TCMSP, http://tcmspw.com/tcmsp.php) to get the ingredients of ASH and the parameters for the selection of active ingredients were set as follows: oral bioavailability (OB) ≥ 30% and drug-likeness (DL) ≥ 0.18 [13] . In addition, potential targets of active ingredients were obtained from the TCMSP database, then the target proteins were imported into the UniProt database (https://www.uniprot.org/) to obtain the gene names.
2.2. Related Targets of LC and Prediction of Potential Targets of ASH against LC
Liver cancer related genes were retrieved from Genecards (https://www.genecards.org/), OMIM (https://omim.org/) and TTD (http://db.idrblab.net/ttd/) database. The search results from each database were combined and duplicates were removed.
Online software Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html) was used to obtain the common targets between ASH and LC. The “drug-ingredient-disease-target” network of ASH anti-LC was constructed by using Cytoscape software (v3.7.2).
2.3. Functional Enrichment Analysis
The common targets were used for GO and KEGG pathway enrichment with the Cluster Profiler package in R software (v3.6.3), and the “ggplot2” package was used to visualize the GO and KEGG enrichment results [14] .
2.4. Protein-Protein Interaction (PPI) Network Construction and Hub Genes Screening
The PPI network was retrieved from STRING Version 11.0 (https://string-db.org/) by selecting Homo sapiens as the organism, and a confidence score > 0.9 (highest confidence) was set as significant [15] . PPI network was then visualized by Cytoscape software (v3.7.2), CytoNCA plugin of Cytoscape was used to screen hub genes based on three criteria: degree centrality (DC), betweenness centrality (BC) and closeness centrality (CC) [16] .
2.5. Molecular Docking Simulation
The molecular docking was performed to further investigate interactions between hub targets and their correspondent active ingredients. The structures of hub target proteins were obtained from the PDB database (https://www.rcsb.org/) and processed with PyMOL software. The 2D structure of the active ingredients was downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov) and energy minimised using MM2 in Chem3D. Afterwards, both the ligand and the receptor were converted to the PDBQT format using AutoDock Tools (ADT 1.5.6), and the docking is performed by running Vina. If the binding energy is less than −5 kJ∙mol−1, it indicated that the target has certain binding activity with the ingredient, the lower the binding energy value, the stronger the binding to the target protein [17] .
3. Results
3.1. Active Compounds and Targets of ASH
A total of 13 active ingredients of ASH were screened from TCMSP database, the correspondent target proteins of each active ingredient were also obtained (Table 1). After removing duplicated targets, we obtained 169 standard gene names (based on UniProt annotation).
3.2. Gene Targets of LC and Potential Targets of ASH against LC
A total of 1004 targets for LC were collected from the GeneCards database by setting relevance score > 20, we identify 493 and 13 LC-related targets from the OMIM and TTD databases. After merging and deleting the duplicate genes, we collected 1410 LC-related target genes.
Venn diagram was generated through the online tool Venny 2.1 and acquired 103 genes about ASH against LC (Figure 1). To further discover the mechanism of ASH against LC, the PPI network of targets for ASH anti-LC was shown as Figure 2.
Table 1. Active ingredients in Artemisiae scopariae Herba (ASH).
Figure 1. Venn diagram of the targets in LC and ASH. Blue indicates LC targets, yellow indicates ASH targets, the intersection in yellow shaded indicates the common targets between LC and ASH.
Figure 2. PPI network of targets for ASH against LC (yellow diamonds represent 13 active ingredients in ASH; blue nodes represent intersection targets between ASH and LC targets).
3.3. Enrichment Analysis
1775 GO terms were obtained with p.adjust < 0.01, there are 1665 terms of biological process (BP), 24 terms of cell composition (CC) and 86 terms of molecular function (MF), accounting for 93.8%, 1.4% and 4.8% respectively. Then the top 5 BP, CC, MF terms are visualized (Table 2, Figure 3).
In total, 131 KEGG pathways were significantly enriched (p.adjust < 0.01), and the top 10 were visualized (Table 3, Figure 4), involved in the prostate cancer, hepatitis B, bladder cancer, kaposi sarcoma-associated herpesvirus infection, hepatitis C, AGE-RAGE signaling pathway in diabetic complications, small cell lung cancer, fluid shear stress and atherosclerosis, IL-17 signaling pathway, human cytomegalovirus infection.
Table 2. Gene Ontology (GO) enrichment.
Figure 3. GO enrichment analysis (adjusted p value < 0.01). The y-axis represents GO categories identified in the GO analysis, and the x-axis represents the significance [−log10 (p.adjust)]. Blue, red, and green boxes represent GO BP, CC, and MF, respectively.
Figure 4. KEGG pathway enrichment analysis (adjusted p value < 0.01). The y-axis shows the pathway term, and the x-axis shows the Gene Ratio of each KEGG pathway. The color indicates the significance [−log10 (p.adjust)], and the size of the circle represents the number of genes enriched in the corresponding pathways.
Table 3. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment.
3.4. PPI Network and Hub Genes
The PPI network was constructed by STRING at a confidence value of 0.9 and it included 100 nodes and 434 edges (Figure 5, Table 4). Nodes which had high degree were identified as the hub nodes in the PPI network, hub genes with degree ≥ 20 were selected (AKT1, TP53, JUN, MAPK1, TNF, RELA, IL6).
Figure 5. Protein-protein interaction (PPI) network of common genes. The nodes represent genes, and the size and color represent the degree value (the higher the degree, the larger the node and the brighter of the color). The thickness and color of the edges display the combined-score between two genes, the higher the score, the wider the edge and the darker of the color.
Table 4. Topological information of 7 hub targets.
3.5. Molecular Docking
Molecular docking results showed that hub protein targets and active ingredients showed good binding interactions (Table 5). Among them, the compound quercetin and MAPK1 displayed the lowest binding energy (−8.9), which suggests that quercetin demonstrated the best docking score against MAPK1 (Figure 6).
Table 5. The binding energy of the hub targets bound to the active ingredients.
Figure 6. MAPK1-quercetin docking. Quercetin forms 4 hydrophobic interactions with VAL37, ILE82, LEU154, forms 5 H-bonds with LYS52, GLN103, MET106, LYS112, ASP165.
4. Discussion
In this study, a total of 13 active ingredients were obtained from Artemisiae scopariae Herba, including flavonoids, chromones and phytosterols. Flavonoids are well known for their physiological anti-inflammatory and antitumor activities, flavonoids may modulate almost all key processes involved in carcinogenesis including apoptosis, proliferation, angiogenesis and metastatic progression [18] [19] . In recent years, various flavonoids have been recognized as having potential protective activity against artificially induced-liver damage [20] . Chromones have been reported to possess antimicrobial, antiviral, and antitumoral activities and the ability to inhibit several enzymes, it also affects the function and activity of liver-metabolizing enzymes [21] . Phytosterols possess hepatoprotective effect, and the anti-cancer effect of phytosterols are achieved by inhibition of cell cycle progression, promotion of cellular apoptosis, inhibition of cell invasion, migration and adhesion, as well as stimulation of the immune function [22] .
The PPI network showed that the targets of ASH against LC do not work alone but instead is a complex interconnected network, according to degree value we finally screened out 7 hub genes (AKT1, TP53, JUN, MAPK1, TNF, RELA, IL6). AKT1 belongs to the family of serine/threonine protein kinases (AKT1, AKT2, and AKT3) known as AKT kinases, AKT is closely associated with cell survival, proliferation, apoptosis, migration and angiogenesis in hepatocellular carcinoma (HCC) [23] [24] . AKT1 participates in the initiation, progression and metastasis of malignant tumors, silencing AKT1 significantly stimulated apoptosis and suppressed the cell cycle, whereas increasing AKT1 expression promoted HCC cells proliferation [24] .
TP53 is the most widely studied tumor suppressor gene, playing an important role in inhibiting tumor development, the function of it is to inhibit cell proliferation in response to DNA damage. By regulating target genes, TP53 induces a variety of cellular responses, including growth arrest, senescence, and apoptosis [25] . JUN encodes c-Jun protein which has essential influence in cell proliferation, survival, and death [26] [27] . Hepatitis C virus infection stimulates c-Jun signaling via protein kinase R to promote proliferation of HCC [27] .
MAPK1 plays a votal role in cancer progression, especially in cancer metastasis, and in HCC development, simultaneous activation of the MAPK1 pathways has been shown to enhance cell-cycle progression [28] . Tumor necrosis factor (TNF) is a mediator of the acute phase response in the liver and can initiate proliferation and cause cell death in hepatocytes, participates in many forms of hepatic pathology, including ischemia/reperfusion injury, alcoholic and viral hepatitis, and injury by hepatotoxins [29] [30] . Study has shown that TNF-α expression in HCC is significantly higher than that in normal hepatic tissue, positively related with the proliferation and invasion ability of HCC cells [31] .
RELA, a member of the NF-κB family, work as a potential factor in the onset and progression of cancers through regulating the expression of genes linked to cell proliferation, migration, invasion, etc. [32] [33] . Clinically, RelA expression has been associated with a lower degree of apoptosis and cirrhosis in the livers of patients with hepatitis C, and liver RelA mRNA levels were inversely associated with severe liver damage and mortality [34] . IL-6 is a cytokine produced by various cells serve a key function in the proliferation, apoptosis, recurrence and metastasis of liver cancer cells [35] . High IL6 level linked with patients’ mortality in cirrhotic patients caused by hepatitis B virus and hepatitis C virus infection. Similarly, high serum level of IL6 was associated with liver-related mortality in chronic HCV patients [36] .
5. Conclusion
In this study, we adopted network pharmacology and molecular docking technology to explore the mechanism of ASH anti-LC. The results indicated that ASH may interact with hub genes such as AKT1, TP53, JUN, etc., regulates hepatitis B, hepatitis C and other signaling pathways, which exerts anticancer effects. Although we lacked experimental validation, it also provides theoretical basis for the treatment of LC in the future.
Acknowledgements
The study was supported by the Basic Research Foundation of School of Modern College of Humanities and Sciences of Shanxi Normal University (Grant No. 2020JCYJ19).