Abstract
Tanshinone IIA (Tan-IIA) is the main bioactive component of Chinese herbal medicine salvia miltiorrhiza (Danshen). Sodium sulfonate of Tan-IIA is widely used in the treatment of cardiovascular and cerebrovascular diseases. Tan-IIA also has inhibitory effects on tumor cells such as gastric cancer, but its therapeutic effect and mechanism on human neuroblastoma have not been evaluated, so its pharmacological mechanism is systematically evaluated by the combined method of network pharmacology and molecular docking. PharmMapper and SwissTargetPrediction predicted 331 potential Tan-IIA-related targets, and 1,152 potential neuroblastoma-related targets were obtained from GeneCards, DisGeNET, DrugBank, OMIM and Therapeutic Target databases (TTD), 107 common targets for Tan-IIA and neuroblastoma. Through gene ontology (GO) functional annotation, Kyoto Encyclopedia of Genes and Genomesa (KEGG) pathway enrichment, protein-protein interaction (PPI) network and cytoHubba plug-in, 10 related signal pathways (Pathways in cancer, PI3K-Akt signaling pathway, Prostate cancer, etc.) and 10 hub genes were identified. The results of molecular docking showed that Tan-IIA could interact with 10 targets: GRB2, SRC, EGFR, PTPN1, ESR1, IGF1, MAPK1, PIK3R1, AKT1 and IGF1R. This study analyzed the related pathways and targets of Tan-IIA in the treatment of human neuroblastoma, as well as the potential anticancer and anti-tumor targets and related signaling pathways of Tan-IIA, which provides a reference for us to find and explore effective drugs for the treatment of human neuroblastoma.
1 Introduction
Neuroblastoma (NB) is the most common extracranial solid tumor in childhood and the most common cancer in infancy. It originates from the neural crest progenitor cells of the developing sympathetic nervous system, and includes ganglioneuroblastoma and ganglioneuroma, accounting for 15% of cancer mortality in children (; ). Neuroblastoma usually occurs in the paravertebral sympathetic ganglia or pelvic ganglia of the neck, chest and abdomen, as well as in the chromaffin cells of the adrenal medulla. Its etiological types are different in biological, genetic, clinical and morphological characteristics (). As a common solid tumor in children, NB usually affects preschool children, but studies have shown that in some cases it may also occur in adolescents and adults. Among children with NB, the survival rate has improved slightly over the past few years, but the mortality rate is still close to 50 per cent. Currently, the main treatment methods of neuroblastoma include surgery, radiotherapy, chemotherapy, bone marrow transplantation and so on, but the prognosis is poor (; ), The risk of severe disability and chronic disease is still very high after treatment, and chemotherapy-related toxicity is often observed (). In this case, there is an urgent need to develop new therapeutic methods with high efficiency and low toxicity. Therefore, it is of great significance to deeply study the pathogenesis of NB and find the key targets of its treatment.
Tanshinone IIA (Tan-IIA) is a natural diterpene quinone isolated from the root of salvia miltiorrhiza and is one of the most pharmacologically active components in salvia miltiorrhiza (). In recent years, Tan-IIA has been used in animal and in vitro studies to treat cardiovascular diseases, postmenopausal syndrome, angina pectoris, myocardial infarction, hypertension, hyperlipidemia, acute ischemic stroke, chronic renal diabetes, Alzheimer’s disease and cancer (; ; ; ), Tan-IIA sulfonate is also widely used in the treatment of cardiovascular and cerebrovascular diseases (). In addition, previous studies have shown that Tan-IIA can effectively inhibit the proliferation of many kinds of tumor cells, and some studies have proved that Tan-IIA has a neuroprotective effect (). Currently, studies on Tan-IIA in the treatment of NB have shown that it can inhibit the proliferation of tumor cells, but there is no study on the specific molecular mechanism of its effect. Therefore, based on the network pharmacology research method, the molecular targets and pathways of Tan-IIA in the treatment of NB were predicted and verified by molecular docking. The flow chart of the research route is shown in Figure 1.
FIGURE 1
2 Materials and methods
2.1 Materials
TCMSP (https://tcmsp-e.com/tcmsp.php) (), PubChem (https://pubchem.ncbi.nlm.nih.gov/) (), PharmMapper (http://www.lilab-ecust.cn/pharmmapper/) (; ; Wang et al., 2017), SwissTargetPrediction (http://www.swisstargetprediction.ch/) (; ), GeneCards (https://www.genecards.org/) (), DisGeNET (https://www.disgenet.org/) (), DrugBank (https://go.drugbank.com/) (Wishart et al., 2018), OMIM (https://www.omim.org/) (), TherapeuticTargetDatabase (http://db.idrblab.net/ttd/) (Zhou et al., 2022), STRING (https://cn.string-db.org/) (), Metascape (https://metascape.org/gp/index.html) (Zhou et al., 2019), jVeen (http://jvenn.toulouse.inra.fr/app/index.html) (), Cytoscape (https://cytoscape.org/) (), AutoDock 1.5.6 (https://autodock.scripps.edu/) (Zhang and Sanner, 2019), RCSB Protein Data Bank (http://www.rcsb.org/) (), PyMOL (https://pymol.org/2/) ().
2.2 Collection of pharmacokinetic information of Tan-IIA by TCMSP
TCMSP is a systematic pharmacology database providing information of traditional Chinese medicine about compounds, ADME-related (absorption, distribution, metabolism and excretion) characteristics, targets and diseases. We can find Chinese herbal medicines with potential biological effects through the database. We used ‘Tanshinone IIA’ as the key word and searched for its pharmacological and molecular properties in TCMSP.
2.3 Prediction of related targets of Tan-IIA
PubChem is a public database of small molecular bioactivity data, from which the related files of Tan-IIA are downloaded and uploaded to SwissTargetPrediction and PharmMapper databases to predict their molecular targets. By comparing the molecular similarity with the active compounds of Tan-IIA, the targets related to Tan-IIA were screened. The common potential drug targets predicted by the two databases were selected for further verification.
2.4 Prediction of related targets in human neuroblastoma
Acquisition of related targets of human neuroblastoma through GeneCards、DisGeNET、DrugBank、OMIM、TherapeuticTargetDatabase database. These databases contain rich and cutting-edge disease-related goals. All databases were searched with ‘Neuroblastoma’ as the keyword, and common genes related to human neuroblastoma were selected as candidate targets.
2.5 PPI network construction
UniProt protein database was used to standardize the disease target and the potential target of drug composition as GeneSymbol, and then the two were mapped to obtain the potential target of Tan-IIA to interfere with human neuroblastoma. The common targets between drugs and diseases were screened by jveen, and the co-expression targets were obtained.
STRING aims to collect and integrate information about all functional interactions between expressed proteins. The Tan-IIA-related goals and neuroblastoma-related goals obtained from the above steps were analyzed and compared. Then the overlapping genes of the target were extracted and introduced into STRING 11.5, and “Homo sapiens” was selected from the optional species, and the minimum interaction score 0.7 was high confidence score. Then the drug-disease target PPI network was obtained.
2.6 Enrichment analysis of gene function and pathway
The list of coexpressed target genes of Tan-IIA in the treatment of human neuroblastoma was inputted into the Metascape database, and the enrichment analysis of GO and KEGG was carried out.
2.7 Constructing Tan-IIA-targetpathway network
Through the information obtained in “2.5” and “2.6” steps, the Cytoscape3.8.2 software was used to establish the relationship between the active compound tanshinone IIA, candidate targets and the related KEGG pathway (drug-target-pathway network), to intuitively understand the interaction of each node. In addition, the first 10 hub genes were predicted using the cell Hubba plug-in in Cytoscape 3.8.2 and the Tan-IIA-target-pathway network was formed.
2.8 Molecular docking
Molecular docking can simulate the interaction between small molecular ligands and biological macromolecular receptors, and predict the binding mode and affinity between them, to select the leading drugs that can bind to the target and provide guidance for the reasonable optimization of drug molecular structure. It is of great significance for drug research and development. The top 10 gene targets of Tan-IIA- NB were selected and analyzed. The crystal structure of the corresponding gene target protein was downloaded from RCSB Protein Data Bank (PDB) database and preserved in pdb format. The 3D structure of Tan-IIA was obtained by PubChem database and converted into mol2 format file by Chem3D. Removal of water molecules and heteroatoms from receptors with AutoDock Tools 1.5.6 and PyMol 2.5.2, and hydrogen atoms and calculated charges were added, which were stored in pdbqt format (). Then, the conformation of ligand-receptor binding was predicted by AutoDockTools 1.5.6. Binding energy is one of the results of molecular docking, which can be used to screen the potential of ligand-receptor binding, in which the docking result ≤ -5 kcal/mol can be considered to have binding activity between ligand and nuclear receptors (Yang et al., 2021). The conformation with the best binding energy was selected as the final conformation and saved as a pdb format file for imaging with PyMOL 2.5.2.
2.9 Experimental verification
2.9.1 Cell culture and treatment
Human neuroblastoma SH-SY5Y cell line was cultured in DMEM complete medium containing 10% fetal bovine serum and 1% penicillin/streptomycin in 5% CO2 and 95% air humidification incubator. At 37°C, SH-SY5Y cells were adherent cells, and the medium was changed every day until the cells reached 60%–80% density. Digest with 0.25% trypsin solution, and then add DMEM to stop digestion, during passage, the cells were resuscitated to the appropriate density and transferred to the culture flask to continue culture.
2.9.2 Cell proliferation assay
Human neuroblastoma cell SHSY-5Y viability was measured by MTT assay. One day before administration, SHSY-5Y cells in logarithmic growth phase were prepared into a single cell suspension with DMEM complete medium. 3,000 cells per well were inoculated into 96-well plate and cultured at 37°C in 5% CO2 cell incubator. After being treated with Tan ⅡA for 24 h, the supernatant was discarded after being incubated with 5 mg·mL-1 MTT solution 10 μL at 37°C for 4 h. DMSO 150 μL was added to each well. 10 min was shaken at low speed on the shaker, and the absorbance was determined by wavelength 570 nm with multi-function enzyme labeling instrument.
2.9.3 Wound healing migration assay
SHSY-5Y was allowed to grow into full confluence in 6-well plates. A line of SHSY-5Y cells was then scraped away in each well using a pipette tip. Subsequently, cells were washed twice with PBS to remove detached cells. Fresh DMEM medium containing different concentrations of Tan ⅡA (0, 20 and 50 μmolL-1) was added to the scratched monolayers. Images were taken using an inverted microscope (Eclipse TS100, Nikon, Japan) after 12 h of incubation. The migrated cells were observed from three randomly chosen fields and quantified by manual counting. Inhibition percentage was expressed as percentage of the untreated cells (100%).
2.9.4 Detection of lipid peroxide content in tumor cells
BODIPY-C11 581/591 is a fat-soluble fluorescent probe that can indicate the presence of lipid peroxides in living cells. SH-SY5Y cells were inoculated in 24-well plates and incubated at the concentration of 5 μM BODIPY-C11 581/591 in 5% CO2 incubator at 37°C for 30 min. Wash it gently with PBS twice at room temperature, and then use BODIPY-C11 probe to detect the concentration of iron ion and the level of lipid peroxidation under fluorescence microscope or confocal microscope. According to the quantitative analysis of the measured fluorescence signal, the excitation wavelengths of BODIPY-C11 are 581 nm (reduced, red) and 500 nm (oxidized, green).
3 Results and discussion
3.1 Pharmacokinetic information of Tan-IIA
The pharmacological and molecular properties of Tan-IIA are obtained using TCMSP (Table 1), including molecular weight (MW), low lipid-water partition coefficient (AlogP), hydrogen bond donor (Hdon), hydrogen bond receptor (Hacc), oral bioavailability (OB), Caco-2 permeability (Caco-2), blood-brain barrier (BBB) and drug-like properties (DL). Generally speaking, according to the parameter information and screening criteria of TCMSP database, the DL of effective drugs should be greater than 0.18. Clearly, Tan-IIA shows good drug similarity, which is worthy of further study.
TABLE 1
| Name | MW | AlogP | Hdon | Hacc | OB (%) | Caco-2 (nm/s) | BBB | DL |
|---|---|---|---|---|---|---|---|---|
| Tanshinone IIA | 294.37 | 4.66 | 0 | 3 | 49.89 | 1.05 | 0.70 | 0.40 |
Pharmacological and molecular properties of Tanshinone IIA.
3.2 Overlapping targets of Tan-IIA and NB
SwissTargetPrediction and PharmaMapper predicted 331 potential Tan-IIA related targets for further verification and obtained 1152 NB related targets from GeneCards, DisGeNET, DrugBank, OMIM and TTD databases. 107 overlapping genes of the target were identified as coexpression targets (Figure 2).
FIGURE 2
3.3 PPI network analysis
The PPI network analyzed by STRING 11.5 is shown in Figure 3. The proteins in the PPI network are represented by circular nodes with protein 3D structure inside, and the lines between the nodes represent the interaction between proteins. The more lines, the stronger the relationship between the two. A total of 107 genes were input, of which 101 nodes interacted with each other, and the other six nodes (NPY2R, PREP, MMP7, NPY5R, MCHR1 and SLC6A2) were not associated with other nodes.
FIGURE 3
3.4 GO and KEGG pathways enrichment analysis
The functions and pathways of potential genes were analyzed by Metascape database to enrich information. The results of GO enrichment analysis are shown in Figure 4. With p < 0.001 as the screening condition, the first 30 GO items were screened, including 10 biological processes, 10 cellular components and 10 molecular functions. In addition, these items mainly focus on protein phosphorylation, response to hormones, membrane raft, kinase activity and so on. For KEGG pathway analysis, Table 2; Figure 5 list the first 10 related signaling pathways p < 0.001), which are mainly concentrated in proteoglycans in cancer, pathways in cancer, PI3K-Akt signaling pathways, endocrine resistance and so on.
FIGURE 4
TABLE 2
| Pathway | Genes | p-Value |
|---|---|---|
| Pathways in cancer | ABL1, AKT1, AKT2, XIAP, AR, BCL2L1, BMP2, CASP3, CCNA2, CCNE1, CDK2, CDK6, NQO1, EGFR, ESR1, FGFR1, GRB2, GSK3B, GSTM1, GSTP1, HDAC1, HDAC2, HMOX1, HSP90AA1, IGF1, IGF1R, IL2, JAK2, KIT, MDM2, MET, MMP2, MMP9, NOS2, PIK3R1, PPARG, MAPK1, MAPK8, MAPK10, MAP2K1, PTGS2, RARA, RARB, RXRA, TERT | <0.001 |
| PI3K-Akt signaling pathway | AKT1, AKT2, BCL2L1, CCNE1, CDK2, CDK6, EGFR, ERBB4, FGFR1, GRB2, GSK3B, NR4A1, HSP90AA1, IGF1, IGF1R, IL2, INSR, JAK2, KDR, KIT, MDM2, MET, NOS3, PIK3CG, PIK3R1, MAPK1, MAP2K1, RXRA | <0.001 |
| Prostate cancer | AKT1, AKT2, AR, CCNE1, CDK2, EGFR, FGFR1, GRB2, GSK3B, GSTP1, HSP90AA1, IGF1, IGF1R, MDM2, MMP9, PIK3R1, MAPK1, MAP2K1 | <0.001 |
| Proteoglycans in cancer | AKT1, AKT2, CASP3, MAPK14, EGFR, ERBB4, ESR1, FGFR1, GRB2, IGF1, IGF1R, KDR, MDM2, MET, MMP2, MMP9, PIK3R1, MAPK1, MAP2K1, PTPN11, SRC | <0.001 |
| Endocrine resistance | AKT1, AKT2, MAPK14, EGFR, ESR1, GRB2, IGF1, IGF1R, MDM2, MMP2, MMP9, PIK3R1, MAPK1, MAPK8, MAPK10, MAP2K1, SRC | <0.001 |
| Chemical carcinogenesis - reactive oxygen species | ABL1, AKT1, AKT2, CBR1, MAPK14, NQO1, EGFR, GRB2, GSTM1, HMOX1, MET, PIK3R1, MAPK1, MAPK8, MAPK10, MAP2K1, PTPN11, SOD2, SRC, AKR1C3 | <0.001 |
| EGFR tyrosine kinase inhibitor resistance | AKT1, AKT2, BCL2L1, EGFR, GRB2, GSK3B, IGF1, IGF1R, JAK2, KDR, MET, PIK3R1, MAPK1, MAP2K1, SRC | <0.001 |
| Progesterone-mediated oocyte maturation | AKT1, AKT2, CCNA2, CDK2, MAPK14, HSP90AA1, IGF1, IGF1R, PDE3B, PIK3R1, PLK1, MAPK1, MAPK8, MAPK10, MAP2K1, AURKA | <0.001 |
| Ras signaling pathway | ABL1, AKT1, AKT2, BCL2L1, EGFR, FGFR1, GRB2, IGF1, IGF1R, INSR, KDR, KIT, MET, PIK3R1, PLA2G2A, MAPK1, MAPK8, MAPK10, MAP2K1, PTPN11 | <0.001 |
| FoxO signaling pathway | AKT1, AKT2, CDK2, MAPK14, EGFR, GRB2, IGF1, IGF1R, INSR, MDM2, PIK3R1, PLK1, MAPK1, MAPK8, MAPK10, MAP2K1, SOD2 | <0.001 |
KEGG enrichment analysis and related genes.
FIGURE 5
3.5 Construction of Tan-IIA-target-NB network
For KEGG pathway analysis, Table 2; Figure 5 list the first 10 related signaling pathways p < 0.001), after deduplicating related targets, the Tan-IIA target pathway network contains a total of 72 nodes (Figure 6). Nodes of different colors and shapes represent different active compounds, targets and pathways, including 1 compound (blue diamond), 61 targets (purple circle) and 10 KEGG pathways (green square). The number of connections between nodes represents the importance of nodes in the network. At the same time, the first 10 hub genes are calculated using the cytoHubba plug-in in Cytoscape3.8.2 and visualized in Figure 7. As shown in Figure 7, the node color represents the node score, and the darker and redder nodes mean the higher score.
FIGURE 6
FIGURE 7
3.6 Molecular docking
Tan-IIA was docked with the first 10 Hub gene targets GRB2, SRC, EGFR, PTPN1, ESR1, IGF1, MAPK1, PIK3R1, AKT1 and IGF1R using AutoDockTools 1.5.6, and 10 groups of ligand-receptors were obtained. Results the binding energy between Tan-IIA and target protein were between-6.47∼-9.2 kcal/mol, the binding energy of each group was less than -5 kcal/mol, and the combination of less than -7 kcal/mol were found in 7 groups, accounting for 70%.
According to the docking results, the binding energies of 10 target proteins to Tan-IIA were all less than -5 kcal/mol (Table 3) and were shown in docking mode 3D map (Figure 8). There are multiple hydrogen bonds between the core target protein and Tan-IIA (yellow dotted lines represent hydrogen bonds). As shown in Figure 8, Tan-IIA binds to the amino acid residues ARG-8 and PRO-6 sites of GRB2; LYS-198 sites of SRC, ASP-238 sites of EGFR; LYS-197 sites of PTPN1; THR-347 sites of ESR1; ASP-59 and CYS-62 sites of IGF1; THR-49, VAL-219 and CYS-221 sites of IGF1R; LYS-54 and GLU-71 sites of MAPK1; PRO-107 sites of PIK3R1; GLN-47, ARG-41 and GLU-40 sites of AKT1; CYS-221, VAL-219 THR-49 site binding.
TABLE 3
| Gene | Binding energy | Gene | Binding energy |
|---|---|---|---|
| ESR1 | −9.2 | IGF1 | −7.36 |
| EGFR | −8.84 | SRC | −7.09 |
| IGF1R | −8.09 | PIK3R1 | −6.96 |
| MAPK1 | −7.79 | AKT1 | −6.93 |
| PTPN1 | −7.39 | GRB2 | −6.47 |
Docking information of Tanshinone IIA with 10 target proteins.
FIGURE 8
3.7 Tan-IIA inhibits the proliferation of SHSY-5Y
The effect of Tan ⅡA on the activity of human neuroblastoma SHSY-5Y cells was detected by MMT assays. Compared with control group, the survival rate of SHSY-5Y cells treated with different doses of Tan ⅡA (0, 20, 50 μmol L-1) for 24 h decreased with the increase of drug concentration, indicating that Tan ⅡA significantly inhibited the proliferation of human neuroblastoma SHSY-5Y cells. The 50% inhibitory concentration (IC50) of Tan ⅡA on SHSY-5Y cells determined by GraphPad Prism 9.0 was 34.98 μmol L-1 (Figure 9). Under the condition of Tan ⅡA intervention for 24 h, the concentration gradient (0, 20, 50 μmol L-1) was selected as the concentration of the follow-up experiment.
FIGURE 9
3.8 Tan-ⅡA inhibits SHSY-5Y migration
We carried out wound healing assays to investigate the effects of Tan ⅡA on cell migration to evaluate the ability of SHSY-5Y in the presence of various concentrations of Tan ⅡA (0, 20 and 50 μmol L-1). The results showed that Tan ⅡA could significantly inhibit the migration of SHSY-5Y in a dose-dependent manner (Figure 10).
FIGURE 10
3.9 Tan-IIA reduces lipid peroxidation of SHSY-5Y
BODIPY-C11 581/591 is commonly used to detect intracellular lipid peroxidation ROS levels. EGCG can reduce the incidence of cancer, neurodegenerative diseases induced by oxidative stress, diabetes and cytokine-induced inflammation in vivo (Yu et al., 2023). In addition, EGCG can significantly reduce the level of intracellular lipid peroxidation by reducing the activity of LDH and the level of ROS, increasing the content of NO and SOD and reducing lipid peroxidation, so it was chosen as the positive control (). SHSY-5Y cells were treated with different concentrations of Tan-IIA (20, 50 μmol L-1) and EGCG (25 μmol L-1). The results in Figure 11 showed that the level of lipid peroxidation in the model group was significantly higher than that in the blank control group (p < 0.01). Compared with the model group, the level of lipid peroxidation in each treatment group was significantly lower than that in the model group (p < 0.05), the effect of 50 μM Tan ⅡA was similar to that of positive control drug EGCG (Table 4).
FIGURE 11
TABLE 4
| Control | Model | Tan ⅡA 20 μM | Tan ⅡA 50 μM | EGCG | |
|---|---|---|---|---|---|
| Lipid peroxidation level | 1.00 | 6.13 ± 0.21 | 4.26 ± 0.35 | 3.01 ± 0.41 | 2.97 ± 0.36 |
Effect of Tan ⅡA on Lipid Peroxidation in SHSY-5Y cells.
4 Discussion
Salvia miltiorrhiza is a kind of traditional Chinese medicine for promoting blood circulation and removing blood stasis. Recent studies have found that Salvia miltiorrhiza plays an important role in improving nervous system function and brain metabolism.
In this study, the network pharmacology and molecular docking methods were used to screen the mechanism of Tan-IIA in the treatment of human neuroblastoma. First of all, by combining multiple drug databases and disease databases, we initially obtained 107 overlapping targets (Figure 2) for the construction of PPI network, and we found that there were protein-protein interactions in 101 nodes (Figure 3). KEGG analysis showed that the pathway of Tan-IIA in the treatment of human neuroblastoma included Pathways in cancer、PI3K-Akt signaling pathway、Prostate cancer、Proteoglycans in cancer、Endocrine resistance、Chemical carcinogenesis - reactive oxygen species、EGFR tyrosine kinase inhibitor resistance、Progesterone-mediated oocyte maturation、Ras signaling pathway、FoxO signaling pathway and so on. It can be seen that these pathways are mainly related to anti-cancer, anti-apoptosis, antioxidant stress and other physiological and pathological processes (), but also closely related to the pathological process of NB. Tan-IIA is associated with many disease targets in the anti-cancer and anti-tumor target network. Through PPI protein network analysis and Cytoscape software visualization, we can see that the key targets of Tan-IIA on human neuroblastoma are GRB2, SRC, EGFR, PTPN1, ESR1, IGF1, MAPK1, PIK3R1, AKT1, IGF1R and so on. GRB2 plays an important role in physiological processes such as cell proliferation and differentiation (; ). At the same time, it is highly expressed in many kinds of tumor tissues, such as breast cancer, leukemia, pancreatic cancer and so on (; ; Yin et al., 2013), and can be located in the nucleus of tumor cells. SRC is a proto-oncogene, when a tumor occurs, SRC will show obvious distribution and shift in the cell, from cytoplasmic distribution to nuclear distribution. At the same time, it plays an important regulatory role in colon, breast, prostate, pancreatic and other tumor tissues (; ). EGFR, the epidermal growth factor receptor, is a member of HER family. EGFR is a receptor tyrosine kinase that participates in the regulation of cell growth, differentiation, survival and metastasis (), Studies have shown that many tumorigenesis, such as lung cancer, lung adenocarcinoma and cholangiocarcinoma, are closely related to overexpression or mutation of EGFR (; ). Similarly, the abnormal expression of PTPN1 in gastric cancer, ovarian cancer, colon cancer and other tumors can promote cell proliferation, colony formation and migration, while reduce the apoptosis of cancer cell lines (; Yu et al., 2019), and the overexpression of PTPN1 is associated with poor prognosis (). ESR1 is a ligand-dependent transcription factor, which can mediate non-genotypic membrane signal pathways to enhance cell proliferation, which is not conducive to the survival of cancer patients. ESR1 is related to the occurrence of gynecological tumors such as breast cancer, endometrial cancer, uterine leiomyoma and so on (Yin et al., 2014; ). IGF1 is an insulin-like growth factor, which can promote scar formation and wound healing, regulate cell proliferation and metabolism, promote tumor cell growth, inhibit tumor cell apoptosis, and promote tumor occurrence and development (). MAPK1 plays an important role in regulating the migration and growth of nerve cells, activating or overexpressing MAPK1 can promote the migration and invasion of tumor cells, increase cell viability and participate in epithelial-mesenchymal transformation (; ). PIK3R1 is a potential protective target, which can negatively regulate the growth and invasion of tumor cells, and is closely related to the occurrence and development of tumor. The decrease of PIK3R1 is significantly correlated with tumor stage, metastasis and prognosis (; ). AKT1 plays an important role in the development, growth and metabolism of mammals. It has been found that the abnormal expression of AKT1 is related to schizophrenia, bipolar disorder, Parkinson’s disease and other diseases. Phosphorylated AKT1 is closely related to the growth, invasion, metastasis, apoptosis and autophagy of tumor cells (; ; ). IGF1R is a tyrosine kinase that is overexpressed in various types of tumor cells and plays an important role in the proliferation and survival of tumor cells. At the same time, the expression of IGF1R in the heart plays an important physiological role in regulating cardiac contractile function and cell proliferation and differentiation (; ). GO enrichment analysis showed that the anticancer effect of Tan-IIA was related to plasma membrane, kinase activity and protein phosphorylation. The plasma membrane is a very thin membrane that surrounds the cell surface, which is mainly composed of membrane lipids and membrane proteins. The basic function of plasma membrane is to maintain the relative stability of intracellular microenvironment. These pathways play a crucial role in cell survival, proliferation, differentiation, and tumorigenesis.
The results of molecular docking of 10 hub genes screened by cytoHubba plug-in showed that Tan-IIA could effectively interact with 10 targets: GRB2, SRC, EGFR, PTPN1, ESR1, IGF1, MAPK1, PIK3R1, AKT1 and IGF1R. The results showed that they all had binding activity, and Tan-IIA had the highest binding or docking score to ESR1 (−9.2 kcal/mol). By analyzing the interaction mode between proteins and ligands, it can be concluded that Tan-IIA can bind these selected target proteins well and obtain lower binding energy mainly by forming hydrogen bonds. Molecular docking models can provide evidence for how compounds act on anti-cancer targets. The results of molecular docking preliminarily confirmed that Tan-IIA can play a therapeutic role in cancer through key targets.
Through the cell experiment, it was observed that Tan-IIA could effectively inhibit the proliferation of SHSY-5Y cells. In order to evaluate the effect of Tan-IIA on the survival of SHSY-5Y cells, its 24 h IC50 was determined to be 34.98 μmol L-1 by MTT experiment, which indicated that Tan-IIA could be a potential drug for the treatment of human neuroblastoma. In addition, scratch assay showed that Tan-IIA inhibited the migration ability of SHSY-5Y cells in a dose-dependent manner. Iron death represents a mode of programmed cell death that has been discovered in recent years, which is different from apoptosis. It is characterized by lipid peroxidation and iron accumulatio (). Lipid peroxidation plays a key role in the initiation and progression of iron death. This process involves oxidative degradation of lipids driven by ROS, which eventually leads to cell damage and rupture of cell membrane. Lipid peroxidation and iron death promote each other, forming a complex interaction, in which ROS promotes iron accumulation, and the same iron accumulation enhances oxidative stress (). The level of intracellular lipid peroxidation was evaluated by BODIPY581/591C11 lipid peroxidation probe (). Iron death plays different roles in different stages of tumor development and inhibits tumor cell proliferation and differentiation in intracellular homeostasis regulation. Therefore, regulating iron death in human neuroblastoma has become a promising treatment strategy. Tan-IIA plays the role of anti-human neuroblastoma through multi-component, multi-target and multi-pathway, so it is an anti-tumor drug with high research value. This study provides a new insight into the possibility of iron death of human neuroblastoma cells induced by Tan-IIA, and provides a new mechanism for potential intervention in the treatment of human neuroblastoma, which needs further molecular function and clinical research.
5 Conclusion
In this study, the mechanism of Tan-IIA in the treatment of human neuroblastoma was studied systematically and comprehensively for the first time by means of network pharmacology and molecular docking. Drug target association prediction preliminarily analyzed the potential anticancer, anti-tumor targets and related signal pathways of Tan-IIA. 10 central genes and their related biological functions, as well as 10 main KEGG pathways were studied by network pharmacology. Ten potential targets for effective binding to Tan-IIA were obtained by molecular docking, which provides a basic theoretical basis for its anti-tumor research, basic experimental verification a basic basis for future experimental verification, and provides a better reference for the development and integration of anti-tumor research of traditional Chinese medicine and modern medicine. This study provides an idea for the research and development of new drugs for the treatment of human neuroblastoma and provides a theoretical basis for the clinical transformation of human neuroblastoma against human neuroblastoma.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
NT: Writing–original draft, Visualization, Validation, Data curation, Conceptualization. YW: Writing–original draft, Funding acquisition, Formal Analysis, Data curation. JM: Writing–original draft, Conceptualization. YZ: Writing–review and editing, Supervision, Software. YC: Writing–review and editing, Supervision, Software. WS: Writing–original draft, Visualization, Data curation. JZ: Writing–original draft, Visualization, Data curation. HS: Writing–review and editing, Funding acquisition. BL: Writing–review and editing, Validation, Funding acquisition, Formal Analysis.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by National Science Fund for Young Scholars of China (Grant No. 82204594).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbraityteA.AukrustP.KouL.AnandI. S.YoungJ.McmurrayJ. J. V.et al (2019). T cell and monocyte/macrophage activation markers associate with adverse outcome, but give limited prognostic value in anemic patients with heart failure: results from RED-HF. Clin. Res. Cardiol.108 (2), 133–141. 10.1007/s00392-018-1331-2
2
AmbergerJ. S.HamoshA. (2017). Searching online mendelian inheritance in man (OMIM): a knowledgebase of human genes and genetic phenotypes. Curr. Protoc. Bioinforma.58, 1. 10.1002/cpbi.27
3
BardouP.MarietteJ.EscudiéF.DjemielC.KloppC. (2014). jvenn: an interactive Venn diagram viewer. BMC Bioinforma.15 (1), 293. 10.1186/1471-2105-15-293
4
BodzinA. S.WeiZ.HurttR.GuT.DoriaC. (2012). Gefitinib resistance in HCC mahlavu cells: upregulation of CD133 expression, activation of IGF-1R signaling pathway, and enhancement of IGF-1R nuclear translocation. J. Cell Physiol.227 (7), 2947–2952. 10.1002/jcp.23041
5
BrizzolaraA.GarbatiP.VellaS.CalderoniM.QuattroneA.ToniniG. P.et al (2020). Co-administration of fendiline hydrochloride enhances chemotherapeutic efficacy of cisplatin in neuroblastoma treatment. Molecules25 (22), 5234. 10.3390/molecules25225234
6
BurleyS. K.BhikadiyaC.BiC.BittrichS.ChenL.CrichlowG. V.et al (2021). RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res.49 (D1), D437–D451. 10.1093/nar/gkaa1038
7
CheungN. K.DyerM. A. (2013). Neuroblastoma: developmental biology, cancer genomics and immunotherapy. Nat. Rev. Cancer13 (6), 397–411. 10.1038/nrc3526
8
DainaA.MichielinO.ZoeteV. (2019). SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res.47 (W1), W357–W364. 10.1093/nar/gkz382
9
DalyR. J.BinderM. D.SutherlandR. L. (1994). Overexpression of the Grb2 gene in human breast cancer cell lines. Oncogene9 (9), 2723–2727.
10
DhahriA. A.AhmadR.RaoA.BhattiD.AhmadS. H.GhufranS.et al (2021). Use of prophylactic steroids to prevent hypocalcemia and voice dysfunction in patients undergoing thyroidectomy: a randomized clinical trial. JAMA Otolaryngol. Head. Neck Surg.147 (10), 866–870. 10.1001/jamaoto.2021.2190
11
Di FulvioM.HenkelsK. M.Gomez-CambroneroJ. (2007). Short-hairpin RNA-mediated stable silencing of Grb2 impairs cell growth and DNA synthesis. Biochem. Biophys. Res. Commun.357 (3), 737–742. 10.1016/j.bbrc.2007.04.013
12
FellerS. M.TuchschererG.VossJ. (2003). High affinity molecules disrupting GRB2 protein complexes as a therapeutic strategy for chronic myelogenous leukaemia. Leuk. Lymphoma44 (3), 411–427. 10.1080/1042819021000037930
13
FengB.XuZ.ZhouF.YuH.SunQ.WangD.et al (2015). Near infrared light-actuated gold nanorods with cisplatin-polypeptide wrapping for targeted therapy of triple negative breast cancer. Nanoscale7 (36), 14854–14864. 10.1039/c5nr03693c
14
GaillardS. L.AndreanoK. J.GayL. M.SteinerM.JorgensenM. S.DavidsonB. A.et al (2019). Constitutively active ESR1 mutations in gynecologic malignancies and clinical response to estrogen-receptor directed therapies. Gynecol. Oncol.154 (1), 199–206. 10.1016/j.ygyno.2019.04.010
15
GaoS.LiuZ.LiH.LittleP. J.LiuP.XuS. (2012). Cardiovascular actions and therapeutic potential of tanshinone IIA. Atherosclerosis220 (1), 3–10. 10.1016/j.atherosclerosis.2011.06.041
16
GfellerD.MichielinO.ZoeteV. (2013). Shaping the interaction landscape of bioactive molecules. Bioinformatics29 (23), 3073–3079. 10.1093/bioinformatics/btt540
17
HowellK. R.LawA. J. (2020). Neurodevelopmental concepts of schizophrenia in the genome-wide association era: AKT/mTOR signaling as a pathological mediator of genetic and environmental programming during development. Schizophr. Res.217, 95–104. 10.1016/j.schres.2019.08.036
18
JavedS.BhattacharyyaS.BaggaR.SrinivasanR. (2020). Insulin growth factor-1 pathway in cervical carcinoma cancer stem cells. Mol. Cell Biochem.473 (1-2), 51–62. 10.1007/s11010-020-03807-6
19
KajarabilleN.Latunde-DadaG. O. (2019). Programmed cell-death by ferroptosis: antioxidants as mitigators. Int. J. Mol. Sci.20, 4968. 10.3390/ijms20194968
20
KasaharaK.NakayamaY.KiharaA.MatsudaD.IkedaK.KugaT.et al (2007). Rapid trafficking of c-Src, a non-palmitoylated Src-family kinase, between the plasma membrane and late endosomes/lysosomes. Exp. Cell Res.313 (12), 2651–2666. 10.1016/j.yexcr.2007.05.001
21
KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2021). PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res.49 (D1), D1388–D1395. 10.1093/nar/gkaa971
22
LeslieN. R.FotiM. (2011). Non-genomic loss of PTEN function in cancer: not in my genes. Trends Pharmacol. Sci.32 (3), 131–140. 10.1016/j.tips.2010.12.005
23
LiuR.LiaoX. Y.PanM. X.TangJ. C.ChenS. F.ZhangY.et al (2019). Glycine exhibits neuroprotective effects in ischemic stroke in rats through the inhibition of M1 microglial polarization via the NF-κB p65/hif-1α signaling pathway. J. Immunol.202 (6), 1704–1714. 10.4049/jimmunol.1801166
24
LiuX.OuyangS.YuB.LiuY.HuangK.GongJ.et al (2010). PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach. Nucleic Acids Res.38, W609–W614. 10.1093/nar/gkq300
25
LouandreC.EzzoukhryZ.GodinC.BarbareJ. C.MaziereJ. C.ChauffertB.et al (2013). Iron-dependent cell death of hepatocellular carcinoma cells exposed to sorafenib. Int. J. Cancer133, 1732–1742. 10.1002/ijc.28159
26
LuB. L.LiJ.ZhouJ.LiW. W.WuH. F. (2016). Tanshinone IIA decreases the levels of inflammation induced by Aβ1-42 in brain tissues of Alzheimer's disease model rats. Neuroreport27 (12), 883–893. 10.1097/WNR.0000000000000618
27
LukschR.CastellaniM. R.ColliniP.De BernardiB.ConteM.GambiniC.et al (2016). Neuroblastoma (Peripheral neuroblastic tumours). Crit. Rev. Oncol. Hematol.107, 163–181. 10.1016/j.critrevonc.2016.10.001
28
MarisJ. M. (2010). Recent advances in neuroblastoma. N. Engl. J. Med.362 (23), 2202–2211. 10.1056/NEJMra0804577
29
MarisJ. M.HogartyM. D.BagatellR.CohnS. L. (2007). Neuroblastoma. Lancet369 (9579), 2106–2120. 10.1016/S0140-6736(07)60983-0
30
MatthayK. K.MarisJ. M.SchleiermacherG.NakagawaraA.MackallC. L.DillerL.et al (2016). Neuroblastoma. Nat. Rev. Dis. Prim.2, 16078. 10.1038/nrdp.2016.78
31
MillischerV.MathesonG. J.MartinssonL.Römer EkI.SchallingM.LavebrattC.et al (2020). AKT1 and genetic vulnerability to bipolar disorder. Psychiatry Res.284, 112677. 10.1016/j.psychres.2019.112677
32
NoiM.MukaishoK. I.MurakamiS.KoshinumaS.MachidaY.YamoriM.et al (2020). Expressions of ezrin, ERK, STAT3, and AKT in tongue cancer and association with tumor characteristics and patient survival. Clin. Exp. Dent. Res.6 (4), 420–427. 10.1002/cre2.293
33
OrsolicN.SirovinaD.GajskiG.Garaj-VrhovacV.Jazvinscak JembrekM.KosalecI. (2013). Assessment of DNA damage and lipid peroxidation in diabetic mice: effects of propolis and epigallocatechin gallate (EGCG). Mutat. Res.757, 36–44. 10.1016/j.mrgentox.2013.04.022
34
PapanastasiouA. D.SirinianC.PlakoulaE.ZolotaV.ZarkadisI. K.KalofonosH. P. (2017). RANK and EGFR in invasive breast carcinoma. Cancer Genet.216-217, 61–66. 10.1016/j.cancergen.2017.07.004
35
Pineiro-HermidaS.Alfaro-ArnedoE.GregoryJ. A.TorrensR.Ruíz-MartínezC.AdnerM.et al (2017). Characterization of the acute inflammatory profile and resolution of airway inflammation after Igf1r-gene targeting in a murine model of HDM-induced asthma. PLoS One12 (12), e0190159. 10.1371/journal.pone.0190159
36
PineroJ.SaüchJ.SanzF.FurlongL. I. (2021). The DisGeNET cytoscape app: exploring and visualizing disease genomics data. Comput. Struct. Biotechnol. J.19, 2960–2967. 10.1016/j.csbj.2021.05.015
37
RadtkeD.LacherS. M.SzumilasN.SandrockL.AckermannJ.NitschkeL.et al (2016). Grb2 is important for T cell development, Th cell differentiation, and induction of experimental autoimmune encephalomyelitis. J. Immunol.196 (7), 2995–3005. 10.4049/jimmunol.1501764
38
RuJ.LiP.WangJ.ZhouW.LiB.HuangC.et al (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform6, 13. 10.1186/1758-2946-6-13
39
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res.13 (11), 2498–2504. 10.1101/gr.1239303
40
ShiW.LiC.WartmannT.KahlertC.DuR.PerrakisA.et al (2022). Sensory ion channel candidates inform on the clinical course of pancreatic cancer and present potential targets for repurposing of FDA-approved agents. J. Pers. Med.12 (3), 478. 10.3390/jpm12030478
41
SigismundS.AvanzatoD.LanzettiL. (2018). Emerging functions of the EGFR in cancer. Mol. Oncol.12 (1), 3–20. 10.1002/1878-0261.12155
42
StelzerG.RosenN.PlaschkesI.ZimmermanS.TwikM.FishilevichS.et al (2016). The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinforma.54, 1 30 1–1 30 33. 10.1002/cpbi.5
43
SuL. J.ZhangJ. H.GomezH.MuruganR.HongX.XuD.et al (2019). Reactive oxygen species-induced lipid peroxidation in apoptosis, autophagy, and ferroptosis. Oxid. Med. Cell Longev.2019, 5080843. 10.1155/2019/5080843
44
SummyJ. M.GallickG. E. (2003). Src family kinases in tumor progression and metastasis. Cancer Metastasis Rev.22 (4), 337–358. 10.1023/a:1023772912750
45
SwainS. S.SinghS. R.SahooA.PandaP. K.HussainT.PatiS. (2022). Integrated bioinformatics-cheminformatics approach toward locating pseudo-potential antiviral marine alkaloids against SARS-CoV-2-Mpro. Proteins90 (9), 1617–1633. 10.1002/prot.26341
46
SzklarczykD.GableA. L.NastouK. C.LyonD.KirschR.PyysaloS.et al (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res.49 (D1), D605–D612. 10.1093/nar/gkaa1074
47
TaniguchiC. M.WinnayJ.KondoT.BronsonR. T.GuimaraesA. R.AlemánJ. O.et al (2010). The phosphoinositide 3-kinase regulatory subunit p85alpha can exert tumor suppressor properties through negative regulation of growth factor signaling. Cancer Res.70 (13), 5305–5315. 10.1158/0008-5472.CAN-09-3399
48
ThamrongwaranggoonU.DetaryaM.SeubwaiW.SaengboonmeeC.HinoS.KogaT.et al (2022). Lactic acidosis promotes aggressive features of cholangiocarcinoma cells via upregulating ALDH1A3 expression through EGFR axis. Life Sci.302, 120648. 10.1016/j.lfs.2022.120648
49
WangK.JiW.YuY.LiZ.NiuX.XiaW.et al (2020). Correction: FGFR1-ERK1/2-SOX2 axis promotes cell proliferation, epithelial-mesenchymal transition, and metastasis in FGFR1-amplified lung cancer. Oncogene39 (42), 6619–6620. 10.1038/s41388-020-01441-6
50
WangX.PanC.GongJ.LiuX.LiH. (2016). Enhancing the enrichment of pharmacophore-based target prediction for the polypharmacological profiles of drugs. J. Chem. Inf. Model56 (6), 1175–1183. 10.1021/acs.jcim.5b00690
51
WangX.ShenY.WangS.LiS.ZhangW.LiuX.et al (2017). PharmMapper 2017 update: a web server for potential drug target identification with a comprehensive target pharmacophore database. Nucleic Acids Res.45 (W1), W356–W360. 10.1093/nar/gkx374
52
WishartD. S.FeunangY. D.GuoA. C.LoE. J.MarcuA.GrantJ. R.et al (2018). DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res.46 (D1), D1074–D1082. 10.1093/nar/gkx1037
53
YangY.HeY.WeiX.WanH.DingZ.YangJ.et al (2021). Network pharmacology and molecular docking-based mechanism study to reveal the protective effect of salvianolic acid C in a rat model of ischemic stroke. Front. Pharmacol.12, 799448. 10.3389/fphar.2021.799448
54
YinJ.CaiZ.ZhangL.ZhangJ.HeX.DuX.et al (2013). A recombined fusion protein PTD-Grb2-SH2 inhibits the proliferation of breast cancer cells in vitro. Int. J. Oncol.42 (3), 1061–1069. 10.3892/ijo.2013.1768
55
YinL.ZhangX. T.BianX. W.GuoY. M.WangZ. Y. (2014). Disruption of the ER-α36-EGFR/HER2 positive regulatory loops restores tamoxifen sensitivity in tamoxifen resistance breast cancer cells. PLoS One9 (9), e107369. 10.1371/journal.pone.0107369
56
YuM.LiuZ.LiuY.ZhouX.SunF.LiuY.et al (2019). PTP1B markedly promotes breast cancer progression and is regulated by miR-193a-3p. FEBS J.286 (6), 1136–1153. 10.1111/febs.14724
57
YuQ.ZhangN.GanX.ChenL.WangR.LiangR.et al (2023). EGCG attenuated acute myocardial infarction by inhibiting ferroptosis via miR-450b-5p/ACSL4 axis. Phytomedicine119, 154999. 10.1016/j.phymed.2023.154999
58
ZhangY.SannerM. F. (2019). AutoDock CrankPep: combining folding and docking to predict protein-peptide complexes. Bioinformatics35 (24), 5121–5127. 10.1093/bioinformatics/btz459
59
ZhouY.ZhangY.LianX.LiF.WangC.ZhuF.et al (2022). Therapeutic target database update 2022: facilitating drug discovery with enriched comparative data of targeted agents. Nucleic Acids Res.50 (D1), D1398–D1407. 10.1093/nar/gkab953
60
ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun.10 (1), 1523. 10.1038/s41467-019-09234-6
Summary
Keywords
neuroblastoma, tanshinone IIA, network pharmacology, molecular docking, gene ontology
Citation
Tang N, Wang Y, Miao J, Zhao Y, Cao Y, Sun W, Zhang J, Sui H and Li B (2024) Potential pharmacological mechanisms of tanshinone IIA in the treatment of human neuroblastoma based on network pharmacological and molecular docking Technology. Front. Pharmacol. 15:1363415. doi: 10.3389/fphar.2024.1363415
Received
10 January 2024
Accepted
29 February 2024
Published
11 March 2024
Volume
15 - 2024
Edited by
Yurong Lai, Gilead), United States
Reviewed by
Kaiyun Guan, University of North Carolina at Chapel Hill, United States
Jiachen Liu, Washington University in St. Louis, United States
Updates
Copyright
© 2024 Tang, Wang, Miao, Zhao, Cao, Sun, Zhang, Sui and Li.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bing Li, 790651120@qq.com; Hua Sui, dianasuihua@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.