14 September 2020: Database Analysis
Identification of Target Genes Related to Sulfasalazine in Triple-Negative Breast Cancer Through Network PharmacologyHaochen Yu1ABCEF, Ke Hu2ABCF, Tao Zhang1BDF, Haoyu Ren1ACEFG*
Med Sci Monit 2020; 26:e926550
BACKGROUND: The anti-inflammatory drug sulfasalazine (SAS) has been confirmed to inhibit the growth of triple-negative breast cancer (TNBC), but the mechanism is not clear. The aim of this study was to use network pharmacology to find relevant pathways of SAS in TNBC patients.
MATERIAL AND METHODS: Through screening of the GeneCards, CTD, and ParmMapper databases, potential genes related to SAS and TNBC were identified. In addition, gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed using the R programming language. Protein–protein interaction networks were constructed using Cytoscape. The Kaplan-Meier plotter screened genes related to TNBC prognosis. TNBC patient gene expression profiles and clinical data were downloaded from The Cancer Genome Atlas database. A heatmap was generated using the R programming language that presents the expression of potential target genes in patients with TNBC.
RESULTS: Eighty potential target genes were identified through multiple databases. The bioinformatical analyses predicted the interrelationships, potential pathways, and molecular functions of the genes from multiple aspects, which are associated with physiological processes such as the inflammatory response, metabolism of reactive oxygen species (ROS), and regulation of proteins in the matrix metalloproteinase (MMP) family. Survival analysis showed that 12 genes were correlated with TNBC prognosis. Heatmapping showed that genes such as those encoding members of the MMP family were differentially expressed in TNBC tissues and normal tissues.
CONCLUSIONS: Our analysis revealed that the main reasons for the inhibitory effect of SAS on TNBC cells may be inhibition of the inflammatory response and MMP family members and activation of ROS.
Keywords: Gene Regulatory Networks, Matrix Metalloproteinases, Reactive Oxygen Species, sulfasalazine, Triple Negative Breast Neoplasms, Anti-Inflammatory Agents, Non-Steroidal, Computational Biology, Pharmacology, Clinical, Protein Interaction Maps
In 2020, among the top 10 major cancer types represented by new cancer cases and the number of estimated deaths according to gender, breast cancer (BC) will rank first as the cause of new cases of tumors in females, and as a cause of death, it will rank second . BC is the most common type of malignant tumor in females and includes numerous subtypes with high heterogeneity . Approximately one-fourth of breast tumors are “triple-negative” (ER− [estrone receptor]/PR− [progesterone receptor]/HER2− [HER2 receptor tyrosine kinase]) breast cancer (TNBC) and usually have a basal molecular phenotype . Due to chemotherapy resistance and the lack of TNBC-targeted therapeutics, patient prognosis is grim . Meanwhile, an increasing number of drugs not used to treated tumors have been shown to have critical antitumor effects. For example, aspirin can inhibit metastasis of colon cancer cells  and metformin can control the oxidative phosphorylation rate and growth of thyroid cancer in vivo and in vitro . There is a critical need to explore more anti-resistance mechanisms and the anticancer effects of some non-cancer drugs in the field of TNBC. In our previous research, we found that the anti-inflammatory drug SAS could inhibit growth of BC cells, and TNBC cells show particularly sensitivity to SAS. The main mechanism underlying the effect of SAS is production of a large number of reactive oxygen species (ROS), inducing ferroptosis . However, the study also found that in addition to the ferroptosis pathway, there may be other mechanisms that can induce the death of TNBC cells, such as autophagy and apoptosis. To investigate this further, we decided to reexamine the specific reasons why SAS inhibits proliferation of TNBC cells by using a network pharmacological method.
In 2007, Hopkins at the University of Dundee in the United Kingdom proposed the term “network pharmacology”  and then described network pharmacology as a “next-generation drug development model” . Network pharmacology is based on the “disease gene-target drug” interaction network. It systematically and comprehensively observes the influence of drugs on the disease network, thereby revealing the underlying mechanisms of drugs that act on humans. Because of the development of network pharmacology, we have a better understanding of the pharmacological effects of many natural medicines and some traditional Chinese medicines [10,11]. Through network pharmacology, some traditional Chinese medicines also have been discovered for treatment of BC . However, there have been few reports about use of network pharmacology for TNBC and the role of the SAS drug network.
Our previous studies confirmed that SAS leads to a more sensitive growth-inhibiting phenotype in TNBC cells . Moreover, the possibility of multipurpose mechanism pathways was also discovered. Therefore, screening of multiple databases through the new method of network pharmacology is helpful to find new mechanisms underlying the phenotype.
In our research, we used different databases to mine some genes related to SAS and TNBC. A drug network with these genes was constructed and analyzed. New potential targets and pathways involved in the effects of SAS on TNBC were discovered. Some specific genes play important roles in the phenotype. Our research mainly demonstrated the effects of SAS on three major types of genes. The first type inhibits expression of inflammation-related genes to reduce the effect of inflammation on TNBC. Inhibition of the matrix metalloproteinase (MMP) family and increased accumulation of ROS are the main causes of inhibition of growth and proliferation of TNBC. Such findings help us to further study the mechanism of the effects of SAS on TNBC and focus on developing more precise treatment options. This is also a model for exploring the anticancer effects of non-tumor drugs. A schematic diagram showing the analysis procedures used for SAS target gene prediction is shown in Figure 1.
Material and Methods
DATA ACQUISITION AND DATABASE APPLICATION:
The molecular structure of SAS was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/). The GeneCards (https://www.genecards.org/) knowledgebase automatically integrates gene-centric data from approximately 150 web sources, including genomic, transcriptomic, proteomic, genetic, clinical, and functional information. Information on a total of 3 885 TNBC-related genes was retrieved from the GeneCards database using “triple-negative breast cancer” as key words. According to a GeneCards Inferred Functionality Score (GIFtS) greater than 20, we selected the top 19.9% of genes, which comprised a total of 773 genes. The GIFtS algorithm uses the wealth of information in GeneCards annotations to produce scores aimed at predicting the degree of functionality of a gene. The PharmMapper (http://www.lilab-ecust.cn/pharmmapper/) server is a freely accessible web server designed to identify potential target candidates for the given small molecules (drugs, natural products, or other newly discovered compounds with unidentified binding targets) using a pharmacophore mapping approach . The PharmMapper database was used for molecular structural analysis to screen 69 genes related to SAS. The Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) is a robust, publicly available database with toxicogenomic information. It provides manually curated key information about chemical gene/protein interactions and chemical-disease and gene-disease relationships derived from peer-reviewed scientific literature . The database was used to rescreen data for 153 genes related to SAS. Data for patients with TNBC was obtained from The Cancer Genome Atlas (TCGA) database (https://www.cancer.gov/). After analysis, 80 potential targets related to both TNBC and SAS were identified.
GENE FUNCTION AND PATHWAY ENRICHMENT ANALYSIS:
Functional enrichment analyses of 80 potential target genes were performed by the R programming language packages “clusterProfiler,” “org.Hs.eg.db,” “enrichplot,” and “ggplot2.” The gene ontology (GO) categories consisted of biological processes (BPs) and molecular functions (MFs). KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis was performed to determine the pathways significantly associated with 80 potential target genes. GO analysis is a commonly used approach for annotating genes and gene products according to molecular functions, biological pathways, and cellular components . KEGG is a useful resource for systematic analysis of gene functions and related high-level genome functional information.
PPI (PROTEIN–PROTEIN INTERACTION) NETWORK CONSTRUCTION AND MODULE ANALYSIS:
The STRING database version 11.0 (http://string-db.org/) aims to collect and integrate interactions between proteins, including direct (physical) interactions and indirect (functional) interactions , and a combined score >0.95 was set as the cut-off criterion. Then, we used Cytoscape version 3.6.0  to visualize the biomolecular interaction networks of the potential genes. Molecular Complex Detection (MCODE), a plug-in used to produce the best results for calculating correlation levels, was subsequently used to identify the clusters in the network.
ANALYSIS OF POTENTIAL TARGET GENE EXPRESSION IN TNBC:
To analyze the expression of 80 potential target genes in TNBC, we downloaded data regarding expression of genes in 220 TNBC samples from the TCGA database. The “heatmap” package in the R programming language was used for analysis, and the expression heatmap was generated.
SURVIVAL ANALYSIS OF POTENTIAL TARGET GENES:
The Kaplan-Meier plotter (https://kmplot.com/analysis/) is capable of assessing the effect of 54k genes on survival in 21 cancer types. The system includes gene chip and RNA-seq data, and the sources for the databases include Gene Expression Omnibus (GEO), The European Genome-Phenome Archive (EGA), and TCGA . In this system, we set the status of ER, PR and HER2 to negative and identified a total of 255 patients with TNBC. Relapse-free survival (RFS) in patients with TNBC was assessed by Kaplan–Meier plots and the log rank P value; log rank P<0.05 was set as the cut-off criterion.
All statistical tests were performed with R version 3.5. P<0.05 was considered to indicate a statistically significant difference.
SCREENING AND IDENTIFICATION OF TNBC- AND SAS-RELATED GENES:
The molecular structure of SAS was obtained from the PubChem database (Figure 2A). Information a total of 773 TNBC-related genes was retrieved from the GeneCards database. PharmMapper was applied to identify 69 SAS-related genes through pharmacophores. Information on genes related to SAS was sorted from the CTD database a second time, and 153 potential target genes were identified. We first summarized the data for SAS-related target genes in the PharmMapper database and the CTD database, and these data were further intersected with the TNBC-related target genes selected by GeneCards (Figure 2B). Finally, 80 target genes related to SAS and TNBC were identified (Table 1).
GENE FUNCTION AND PATHWAY ENRICHMENT ANALYSES:
As shown in Figure 3A, the top 10 biological processes (BP) included the ROS metabolic process, response to steroid hormone, regulation of ROS metabolic process, cellular response to drugs, response to lipopolysaccharide, regulation of the inflammatory response, cellular response to oxidative stress, response to molecules of bacterial origin, response to oxidative stress, and response to ROS. These functional terms are highly relevant to regulation and metabolism of ROS. This is consistent with our previous findings that TNBC death may be closely related to metabolism of ROS. Moreover, Figure 3A also shows the top 10 cellular components (CC), including the membrane raft, membrane microdomain, membrane region, vesicle lumen, secretory granule lumen, extracellular matrix, cytoplasmic vesicle lumen, protein kinase complex, cyclin-dependent protein kinase holoenzyme complex, and collagen-containing extracellular matrix. This suggests that many cytokines and production of proteins may be involved. In addition, in Figure 3A, we can also see that the molecular functions (MFs) of the 80 genes with relatively high correlations included cytokine activity, receptor ligand activity, cytokine receptor binding, growth factor activity, protease binding, and growth factor receptor binding. Based on the MF results, these 80 genes are closely associated with activation and regulation of various receptors. The relevant signal pathways were discovered by using KEGG analysis. The two pathways that were most significant were the PI3K-Akt signaling pathway and the interleukin (IL)-17 signaling pathway (Figure 3B). The phosphatidylinositol-3-kinase -Akt pathway promotes cell metabolism, proliferation, and survival in response to extracellular signals. The most significant role of IL-17 is its involvement in induction and mediation of pro-inflammatory responses. Both of these important pathways are involved in multiple physiological processes in TNBC cells. The results of the above analysis lead us to believe that SAS inhibits growth of TNBC cells for a variety of possible reasons. The large accumulation of ROS, activation of inflammation, and regulation of cellular metabolism could be among the reasons. Of course, KEGG analysis gives many other pathways as well, which are calculable results, but for the SAS and TNBC projects that we focused on, we felt that the pathways that could cause cell death and inflammation production made the most sense.
PPI NETWORK ANALYSIS:
With interaction and pathway exploration mining among 80 targets based on the STRING database, an entire interaction network was constructed using Cytoscape and a total of 61 nodes and 152 edges were obtained from 80 genes. The yellow nodes correspond to the most relevant target genes (Figure 4A). We used MCODE in Cytoscape to calculate the optimal grouping of the entire network, producing two optimal subgroups (Figure 4B, 4C) that were significantly related to theMMP and IL families. Based on the STRING statistics, the top 15 nodes with the most edges are presented in Figure 4D. From the results of PPI Network analysis, there are more potential pathways between these 15 target genes and other molecules. The MMP and IL families in the two subgroups also showed significant differences in expression in heat map analysis. Among them, CXCL8 and TNF were statistically significant in survival analysis, which proved that they might be associated with prognosis. However, CXCL8 and TNF did not differ significantly in expression of TNBC by heatmap analysis.
IDENTIFICATION OF 80 PROGNOSIS-ASSOCIATED GENES:
Survival analysis was performed with the median as the cut-off criterion. As a result, 12 genes (ABCC2, BCL2, CASP9, CDH1, CDK2, CXCL8, EGFR, GSK3B, IFENG, IGFBP3, TNF, and TYMP) were found to be associated with RFS in patients with TNBC. As shown in Figure 5, patients with high expression of ABCC2, CDH1, CDK2, CXCL8, GSK3B, and IGFBP3 had shorter RFS. In contrast, low expression of the remaining six genes predicted poor survival. Among these prognosis-related genes, BCL2, EGFR,and IFENG were shown to be associated with ROS metabolism in GO analysis and TNF and EGFR were shown to be associated with inflammatory pathways. The KEGG shows that BCL2 is AKT pathway-relevant and regulates cell proliferation, cycle, and various metabolic processes. Expression of BCL2 and CASP9 also differed significantly between TNBC and normal tissues, but there was little difference in expression of other prognosis-related genes screened.
DIFFERENTIAL EXPRESSION OF 80 GENES IN TNBC AND NORMAL TISSUES:
Differential expression analysis of genes in 220 TNBC samples was performed using the “heatmap” software package in the R programming language. As shown in Figure 6, we found that several genes were significantly differentially expressed. MMP family genes were significantly overexpressed in TNBC tissues. CCNE1 and ERIC5 were also highly overexpressed in TNBC. However, levels of expression of peroxisome proliferator activated receptor gamma (PPARG), BCL2, TLR4, and FABP4 were significantly lower in TNBC than in normal tissues, suggesting that they may be involved in suppressing cancer development. From the results of the expression analysis, we can once again focus on BCL2, who also showed some differences in the survival analysis. Similarly, various gene expressions in the MMP family showed significant differences. Together with results of previous PPI analyses, we have more reason to believe that the MMP family has an important role in cell death caused by SAS.
Considering that TNBC lacks ER expression and effective therapeutic targets, chemotherapy is indispensable in clinical practice. The chemotherapy regimens used for TNBC mainly include paclitaxel and anthracyclines . However, chemotherapy leads to many inevitable problems, such as resistance and bone marrow suppression. Therefore, it is worth investigating the new use of old drugs. Our previous research found that SAS has a significant inhibitory effect on TNBC cells, mainly by activating the ferroptosis pathway to produce a large amount of ROS to cause cell death . In addition, it has also been shown that there may be many other potential pathways involved in generation of ROS. In this study, a rigorous network pharmacological analysis of 80 genes related to SAS and TNBC was conducted, which provides an unprecedented opportunity to determine some of the relevant pathways and indicates an appropriate direction for future research.
An increasing number of studies have shown that inflammation plays an important role in tumorigenesis, and TNBC is no exception . Interestingly, we found that some of the 80 genes identified in this study, such as TNF and NF-κB are significantly related to inflammation, as are many IL family genes. In general, SAS is broadly used to treat chronic inflammation in the gut, joints, and retina. In addition to the obvious inhibitory effect of SAS on inflammation, some studies have also shown a significant inhibitory effect on the NF-κB signaling pathway . The NF-κB transcription factor family members isare well-documented critical components involved in cell survival and proliferation, immunity, and inflammation, thereby acting as important mediators of cancer emergence, progression, and metastasis . In a cancer setting, tumors are often characterized by elevated levels of cytokines that are produced and secreted by the classical pathway, including two proinflammatory cytokines, TNF and IL-1B, which commonly trigger tumor activation. Indeed, emerging evidence suggests that TNF fuels progression of breast cancer by promoting proliferation, transformation, angiogenesis, invasion, and metastasis  as well as by being vital for viability of TNBC cell lines . Inhibition of TNF by SAS was proven many years ago , so expression or activity of some molecules downstream of TNF, such as RELA, also will be affected . On the other hand, SAS has a significant inhibitory effect on both IL-1B [28,29] and IL-6 [30,31]. Therefore, we can generate a scientific hypothesis stating that SAS inhibits the expression or activity of inflammation-related factors and thus inhibits the occurrence and development of TNBC. However, not enough basic research has been conducted to confirm this hypothesis.
MMPs are a group of Zn-dependent neutral peptide endonucleases characterized by transmembrane transport and secretion , the main role of which are to degrade the extracellular matrix . Recently, mechanisms by which MMPs promote tumorigenesis and development have become relatively well known, such as degradation of the extracellular matrix, activation of cell growth factors, and regulation of angiogenesis. Interactions between various mechanisms jointly stimulate the occurrence and development of tumors. MMPs also exert many necessary functions in breast cancer and have negative prognostic value for patients with breast cancer . This is consistent with the results we obtained in the heatmap analysis. MMPs are mainly produced by interstitial fibroblasts, and tumor progression and poor prognosis are related to expression of MMP1, MMP7, and MMP12 in the interstitium . MMPs have been studied in breast cancer subtypes but rarely have been studied in TNBC. In TNBC, patients with high expression of CD147 and MMP9 have worse prognosis than patients with low expression of those 2 genes . Many relevant studies have been performed on MMPs and SAS, and we found that SAS can inhibit expression and activity of MMPs. SAS inhibits expression and activity of MMP1, MMP2, MMP3, and MMP9, which are activated by IL-1B . Lakey RL et al. found that mRNA expression of MMP13 increased under IL-1A cotreatment with oncostatin M (OSM). However, SAS can inhibit this phenomenon . Therefore, there are reasons to speculate that SAS can affect degradation of the extracellular matrix, activation of cell growth factors, and regulation of angiogenesis by inhibiting expression and activity of MMPs in TNBC.
A large number of ROS-related biological processes were found to be associated with the 80 genes by GO analysis. ROS produced by peroxidation can affect various cell functions. Therefore, we focused on the impact of ROS on TNBC. ROS are known to play a controversial role during cancer initiation and progression. Studies of ROS have not come to consistent conclusions regarding TNBC. A previous study found increased ROS levels in TNBC cell lines and a dependency on ROS for survival, because antioxidant treatment induced cell death in TNBC cells but not in an ER -positive cell line . In contrast, there is also literature showing that the decrease in TNBC growth is due to the large accumulation of ROS . Our study found that SAS promotes expression of TFRC by inhibiting the expression of SLC7A11, thereby inducing ferroptosis in MDA-MB-231 TNBC cells and producing a large amount of ROS. It was also found in our previous study that reducing ROS accumulation in SAS-treated TNBC cells with ferroptosis inhibitors can reverse the decrease in cell proliferation .
ROS also plays an important role in TNBC radiotherapy. In combination with radiation, niclosamide treatment resulted in a significant increase in ROS generation and induction of apoptosis in parental and radioresistant TNBC cells in vitro and TNBC xenograft tumors in vivo . Moreover, ROS have been reported to be associated with chemotherapy resistance in TNBC. Lu Lu et al. demonstrated that MYC and MCL1 cooperated in maintaining chemotherapy-resistant cancer stem cells (CSCs) in TNBC. MYC and MCL1 increased mitochondrial oxidative phosphorylation (mtOXPHOS) and generation of ROS, which are involved in maintenance of CSCs . Previous results have shown that SAS increases accumulation of ROS, which is the main cause of inhibition of growth and proliferation of TNBC. However, some specific mechanisms are still unclear. From our GO analysis, approximately 25% of genes were shown to be associated with regulation of ROS. Here, we selected some genes from among 80 potential target genes for analysis. From the results, ANXA1, NOS2 and NQO1 were the genes that appeared most frequently in ROS-related biological processes. Giovanna Leoni et al. described a novel signaling pathway that is activated by an endogenous formyl peptide receptor ligand and ANXA1, which mediate activation of ROS by an epithelial nicotinamide adenine dinucleotide phosphate oxidase . Bhaskar S et al. found that copper sulfate results in increased lipoprotein oxidation which results in increased expression of NOS2 mRNA . Studies on NQO1 and ROS have shown a continued increasing trend. In 2018, Luo S et al. found in a study of glioblastoma multiforme that knocking out NQO1 increased the number of ROS and decreased cell proliferation. In contrast, overexpression of NQO1 attenuates ROS and increases cell proliferation . Treatment with the PPARG agonist pioglitazone can trigger changes in metabolic patterns, change the cellular energy supply pattern from glucose oxidation to fatty acid oxidation, and inhibit pyruvate oxidation to reduce glutathione levels. These metabolic changes cause a significant increase in ROS levels in cancer cells . Similar to pioglitazone, SAS can act as an agonist of PPARG, increasing PPARG expression . Therefore, this also may be the explanation for the large amount of ROS produced in SAS-treated tumor cells. Both nuclear factor, erythroid 2 such as 2 (NFELE2) and NAD(P)H quinone dehydrogenase 1 (NQO1) are important factors involved in the oxidative peroxidation reaction, and SAS can generate a large amount of ROS by increasing their expression and activity . Consequently, SAS can damage cells by peroxidation by NFELE2 and NQO1. This could be another reason for accumulation of ROS in SAS-treated cells. In short, SAS may stimulate accumulation of ROS in cells through different pathways, thereby inhibiting the growth of TNBC cells.
However, we must admit that our research still has some limitations. First, our research lacks the support of basic experiments. In addition, there are not enough data supporting clinical use of SAS in patients with TNBC. Furthermore, in the survival analysis, we did not have sufficient data to explore the effect of the 80 genes on the rate of overall survival (OS) in patients with TNBC. We chose RFS as the indicator. RFS mainly reflects recurrence and progression of the disease, but the best indicator for the survival time of patients is still OS.
In summary, we first used different databases to identify 80 genes linked with a noncancer drug (SAS) and the BC type with the worst prognosis (TNBC). Finally, we found that SAS may reduce development of TNBC by inhibiting expression of inflammation-related genes and MMP-related genes. Simultaneously stimulating expression of genes related to the oxidative stress pathway and increasing the accumulation of intracellular ROS may also be one of the potential mechanisms by which SAS suppresses cancer (Figure 7). Therefore, our research can set the direction for follow-up research on SAS and provide a theoretical basis for clinical treatment of SAS in TNBC.
FiguresFigure 1. Flowchart of screening of 80 genes. Figure 2. (A) Chemical structure of sulfasalazine that was downloaded from the PubChem database. (B) Identification of potential target genes associated with both SAS and TNBC. Figure 3. (A, B) Functional and pathway enrichment analysis of genes related to SAS and TNBC. Figure 4. (A) Protein–protein interaction (PPI) network of 80 genes associated with SAS and TNBC. The 8 yellow nodes indicate the 8 genes that are connected to the highest number of other genes. (B, C) From the entire PPI network, 2 significant subgroups were identified by MCODE in Cytoscape. (D) The top 15 nodes with the most connections among 80 nodes. Figure 5. The 12 genes identified by the Kaplan-Meier survival analysis were significantly associated with relapse-free survival (RFS) among the 80 genes. Figure 6. Heatmap of 80 differentially expressed genes. Red color represents upregulated genes, and green represents downregulated genes in TNBC patients and normal patients. Figure 7. Schematic diagram of potential pathways for SAS to inhibit growth and proliferation of TNBC cells.
1. Siegel RL, Miller KD, Jemal A, Cancer statistics: Cancer J Clin, 2020; 70; 7-30
2. Polyak K, Heterogeneity in breast cancer: J Clin Invest, 2011; 121; 3786-88
3. Weigelt B, Reis-Filho JS, Histological and molecular types of breast cancer: Is there a unifying taxonomy?: Nat Rev Clin Oncol, 2009; 6; 718-30
4. Kim C, Gao R, Sei E, Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing: Cell, 2018; 173; 879-893.e13
5. Kagbo-Kue S, Ajose T, Bakinde N, Aspirin inhibited the metastasis of colon cancer cells by inhibiting the expression of toll-like receptor 4: Cell Biosci, 2018; 8; 36
6. Thakur S, Daley B, Gaskins K: Clin Cancer Res, 2018; 24; 4030-43
7. Yu H, Yang C, Jian J, Sulfasalazine induced ferroptosis in breast cancer cells is reduced by the inhibitory effect of estrogen receptor on the transferrin receptor: Oncol Rep, 2019; 42; 826-38
8. Hopkins AL, Network pharmacology: Nat Biotechnol, 2007; 25; 1110-11
9. Hopkins AL, Network pharmacology: The next paradigm in drug discovery: Nat Chem Biol, 2008; 4; 682-90
10. Zhang Z, Yi P, Yang J, Integrated network pharmacology analysis and serum metabolomics to reveal the cognitive improvement effect of Bushen Tiansui formula on Alzheimer’s disease: J Ethnopharmacol, 2020; 249; 112371
11. Luo TT, Lu Y, Yan SK, Network pharmacology in research of chinese medicine formula: methodology, application and prospective: Chin J Integr Med, 2020; 26; 72-80
12. Zeng L, Yang K, Exploring the pharmacological mechanism of Yanghe Decoction on HER2-positive breast cancer by a network pharmacology approach: J Ethnopharmacol, 2017; 199; 68-85
13. Wang X, Shen Y, Wang S, PharmMapper 2017 update: A web server for potential drug target identification with a comprehensive target pharmacophore database: Nucleic Acids Res, 2017; 45; W356-60
14. Davis AP, Grondin CJ, Johnson RJ, The comparative toxicogenomics database: Update 2017: Nucleic Acids Res, 2017; 45; D972-78
15. The Gene Ontology Consortium, Expansion of the Gene Ontology knowledgebase and resources: Nucleic Acids Res, 2017; 45; D331-38
16. Kanehisa M, Furumichi M, Tanabe M, KEGG: New perspectives on genomes, pathways, diseases and drugs: Nucleic Acids Res, 2017; 45; D353-61
17. Szklarczyk D, Morris JH, Cook H, The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible: Nucleic Acids Res, 2017; 45; D362-68
18. Shannon P, Markiel A, Ozier O, Cytoscape: A software environment for integrated models of biomolecular interaction networks: Genome Res, 2003; 13; 2498-504
19. Gyorffy B, Lanczky A, Eklund AC, An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients: Breast Cancer Res Treat, 2010; 123; 725-31
20. Sikov WM, Berry DA, Perou CM, Impact of the addition of carboplatin and/or bevacizumab to neoadjuvant once-per-week paclitaxel followed by dose-dense doxorubicin and cyclophosphamide on pathologic complete response rates in stage II to III triple-negative breast cancer: CALGB 40603 (Alliance): J Clin Oncol, 2015; 33; 13-21
21. Matsumoto H, Koo SL, Dent R, Role of inflammatory infiltrates in triple negative breast cancer: J Clin Pathol, 2015; 68; 506-10
22. Wahl C, Liptay S, Adler G, Schmid RM, Sulfasalazine: A potent and specific inhibitor of nuclear factor kappa B: J Clin Invest, 1998; 101; 1163-74
23. Shostak K, Chariot A, NF-kappaB, stem cells and breast cancer: The links get stronger: Breast Cancer Res, 2011; 13; 214
24. Barbie TU, Alexe G, Aref AR, Targeting an IKBKE cytokine network impairs triple-negative breast cancer growth: J Clin Invest, 2014; 124; 5411-23
25. Yu M, Zhou X, Niu L, Targeting transmembrane TNF-alpha suppresses breast cancer growth: Cancer Res, 2013; 73; 4061-74
26. Pileczki V, Braicu C, Gherman CD, Berindan-Neagoe I, TNF-alpha gene knockout in triple negative breast cancer cell line induces apoptosis: Int J Mol Sci, 2012; 14; 411-20
27. Bissonnette EY, Enciso JA, Befus AD, Inhibitory effects of sulfasalazine and its metabolites on histamine release and TNF-alpha production by mast cells: J Immunol, 1995; 156; 218-23
28. Soliman NA, Keshk WA, Rizk FH, Ibrahim MA, The possible ameliorative effect of simvastatin versus sulfasalazine on acetic acid induced ulcerative colitis in adult rats: Chem Biol Interact, 2019; 298; 57-65
29. Soubh AA, Abdallah , El-Abhar HS, Geraniol ameliorates TNBS-induced colitis: Involvement of Wnt/beta-catenin, p38MAPK, NFkappaB, and PPARgamma signaling pathways: Life Sci, 2015; 136; 142-50
30. Niu X, Zhang H, Li W, Protective effect of cavidine on acetic acid-induced murine colitis via regulating antioxidant, cytokine profile and NF-kappaB signal transduction pathways: Chem Biol Interact, 2015; 239; 34-45
31. Keerthy HK, Mohan CD, Sivaraman Siveen K: J Biol Chem, 2014; 289; 31879-90
32. Martin MD, Matrisian LM, The other side of MMPs: Protective roles in tumor progression: Cancer Metastasis Rev, 2007; 26; 717-24
33. Kessenbrock K, Plaks V, Werb Z, Matrix metalloproteinases: Regulators of the tumor microenvironment: Cell, 2010; 141; 52-67
34. Cid S, Eiro N, Fernandez B, Prognostic influence of tumor stroma on breast cancer subtypes: Clin Breast Cancer, 2018; 18; e123-33
35. Finak G, Bertos N, Pepin F, Stromal gene expression predicts clinical outcome in breast cancer: Nat Med, 2008; 14; 518-27
36. Zhao S, Ma W, Zhang M, High expression of CD147 and MMP-9 is correlated with poor prognosis of triple-negative breast cancer (TNBC) patients: Med Oncol, 2013; 30; 335
37. Lu Y, Fukuda K, Li Q, Role of nuclear factor-kappaB in interleukin-1-induced collagen degradation by corneal fibroblasts: Exp Eye Res; 83; 560-68 20016
38. Lakey RL, Cawston TE, Sulfasalazine blocks the release of proteoglycan and collagen from cytokine stimulated cartilage and down-regulates metalloproteinases: Rheumatology (Oxford), 2009; 48; 1208-12
39. Sarmiento-Salinas FL, Delgado-Magallon A, Montes-Alvarado JB, Breast cancer subtypes present a differential production of reactive oxygen species (ROS) and susceptibility to antioxidant treatment: Front Oncol, 2019; 9; 480
40. Gong X, Smith JR, Swanson HM, Rubin LP, Carotenoid lutein selectively inhibits breast cancer cell growth and potentiates the effect of chemotherapeutic agents through ROS-mediated mechanisms: Molecules, 2018; 23; 905
41. Lu L, Dong J, Wang L, Activation of STAT3 and Bcl-2 and reduction of reactive oxygen species (ROS) promote radioresistance in breast cancer and overcome of radioresistance with niclosamide: Oncogene, 2018; 37; 5292-304
42. Warne JP, John CD, Christian HC, Gene deletion reveals roles for annexin A1 in the regulation of lipolysis and IL-6 release in epididymal adipose tissue: Am J Physiol Endocrinol Metab, 2006; 291; E1264-73
43. Luo S, Lei K, Xiang D, Ye K, NQO1 is regulated by PTEN in Glioblastoma, mediating cell proliferation and oxidative stress: Oxid Med Cell Longev, 2018 9146528
44. Lee KM, Giltnane JM, Balko JM, MYC and MCL1 cooperatively promote chemotherapy-resistant breast cancer stem cells via regulation of mitochondrial oxidative phosphorylation: Cell Metab, 2017; 26; 633-47.e7
45. Srivastava N, Kollipara RK, Singh DK, Inhibition of cancer cell proliferation by PPARgamma is mediated by a metabolic switch that increases reactive oxygen species levels: Cell Metab, 2014; 20; 650-61
46. Martinez-Augustin O, Lopez-Posadas R, Gonzalez R, Genomic analysis of sulfasalazine effect in experimental colitis is consistent primarily with the modulation of NF-kappaB but not PPAR-gamma signaling: Pharmacogenet Genomics, 2009; 19; 363-72
16 September 2022 : Review articleEffects of Physiotherapy on Rehabilitation and Quality of Life in Patients Hospitalized for COVID-19: A Rev...
Med Sci Monit In Press; DOI: 10.12659/MSM.938141
03 October 2022 : Clinical ResearchEffect of Vitamin D Concentration on Course of COVID-19
Med Sci Monit 2022; 28:e937741
01 October 2022 : EditorialEditorial: Rebound COVID-19 and Cessation of Antiviral Treatment for SARS-CoV-2 with Paxlovid and Molnupiravir
Med Sci Monit 2022; 28:e938532
30 Sep 2022 : Clinical ResearchCorrelation Between Trough Level of Abiraterone and Prostate-Specific Antigen (PSA) Response in Metastatic ...
Med Sci Monit In Press; DOI: 10.12659/MSM.938091
30 Sep 2022 : Review articleA Review of the Biological Mechanisms of Dexmedetomidine for Postoperative Neurocognitive Disorders
Med Sci Monit In Press; DOI: 10.12659/MSM.937862
29 Sep 2022 : Animal ResearchEffect of Antilogous Platelet-Rich Plasma on the Revascularization of Rabbit Prefabricated Flap
Med Sci Monit In Press; DOI: 10.12659/MSM.937718
Most Viewed Current Articles
30 Dec 2021 : Clinical ResearchRetrospective Study of Outcomes and Hospitalization Rates of Patients in Italy with a Confirmed Diagnosis o...
Med Sci Monit 2021; 27:e935379
13 Nov 2021 : Clinical ResearchAcceptance of COVID-19 Vaccination and Its Associated Factors Among Cancer Patients Attending the Oncology ...
Med Sci Monit 2021; 27:e932788
08 Mar 2022 : Review articleA Review of the Potential Roles of Antioxidant and Anti-Inflammatory Pharmacological Approaches for the Man...
Med Sci Monit 2022; 28:e936292