Action mechanism of seaweed-kunbu in treatment of breast cancer based on network pharmacology
-
摘要:目的
通过网络药理学方法分析海藻-昆布治疗乳腺癌的作用机制。
方法通过TCMSP数据库检索海藻-昆布的活性成分及其靶点。通过GeneCards和OMIM数据库获取乳腺癌的相关靶标, 并通过Cytoscape软件构建“药物-疾病靶点”网络图。利用String平台进行蛋白质-蛋白质相互作用(PPI)网络分析。采用R语言进行基因本体论(GO)功能注释和京都基因与基因组百科全书(KEGG)通路富集分析。采用分子对接技术验证关键化合物与靶点的结合能力,并通过体外实验进行功能验证。
结果本研究共获得海藻-昆布的有效化合物11个。海藻-昆布与乳腺癌高度相关的交集靶标为150个,筛选出丝氨酸/苏氨酸激酶1(AKT1)、p53肿瘤蛋白(TP53)、肿瘤坏死因子(TNF)、白细胞介素6(IL6)、血管内皮生长因子A(VEGFA)、JUN、胱天蛋白酶3(CASP3)、白细胞介素1B(IL1B)、MYC和人表皮生长因子受体(EGFR)共10个核心基因。GO分析表明,这些化合物可能参与信号转导、细胞凋亡、细胞增殖、炎症反应等生物过程。KEGG通路分析表明,主要富集信号通路为PI3K/Akt、MAPK和TNF等信号通路。分子对接和体外实验结果表明,海藻-昆布的主要化学成分岩藻甾醇对乳腺癌三阴型细胞BT-549具有抗增殖作用。
结论海藻-昆布可能通过多靶点、多通路来治疗乳腺癌,其主要化学成分岩藻甾醇可有效抑制BT-549细胞增殖活性。
Abstract:ObjectiveTo analyze the mechanism of action of seaweed-kunbu in the treatment of breast cancer through network pharmacology.
MethodsThe active ingredients and their targets of seaweed-kunbu were retrieved through the TCMSP database. Relevant targets for breast cancer were obtained through the GeneCards and OMIM databases, and "drug-disease target" network diagram was constructed using Cytoscape software. Protein-protein interaction (PPI) network analysis was performed using the String platform. Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using R language. Molecular docking technology was used to verify the binding ability of key compounds to targets, and functional verification was carried out through in vitro experiments.
ResultsA total of 11 effective compounds were obtained from seaweed-kunbu in this study, among which 10 were related to breast cancer. There were 150 highly relevant intersecting targets between seaweed-kunbu and breast cancer, and 10 core genes including Serine/threonine kinase 1 (AKT1), p53 tumor protein (TP53), tumor necrosis factor (TNF), interleukin6 (IL6), vascular endothelial growth factor A (VEGFA), JUN, caspase 3 (CASP3), interleukin1B (IL1B), MYC and human epidermal growth factor receptor (EGFR) were screened out. GO analysis indicated that these compounds may be involved in biological processes such as signal transduction, apoptosis, cell proliferation, and inflammatory response. KEGG pathway analysis showed that the main enriched signaling pathways were PI3K/Akt, MAPK and TNF. Molecular docking and in vitro experimental results demonstrated that fucosterol, the main chemical component of seaweed-kunbu, had an antiproliferative effect on triple-negative breast cancer BT-549 cells.
ConclusionSeaweed-kunbu may treat breast cancer through multiple targets and pathways, and its main chemical component fucosterol can effectively inhibit the proliferation activity of BT-549 cells.
-
Keywords:
- seaweed-kunbu /
- breast cancer /
- network pharmacology /
- molecular docking /
- signaling pathway /
- fucosterol
-
阿尔茨海默病(AD)是一种以β-淀粉样蛋白沉积和微管相关蛋白tau异常为特征的进行性神经退行性疾病,临床上表现为记忆和认知功能障碍、语言和行为能力缺失,患者后期表现出严重的失忆症和运动功能减弱,最终导致死亡[1]。尽管AD的特征性病理改变已明确,但靶向这些病理蛋白的药物疗效欠佳,因此进一步筛查有效靶点非常必要。本研究基于生物信息学分析筛选影响AD发生的核心基因和相关信号通路,现报告如下。
阿尔茨海默病(AD)是一种以β-淀粉样蛋白沉积和微管相关蛋白tau异常为特征的进行性神经退行性疾病,临床上表现为记忆和认知功能障碍、语言和行为能力缺失,患者后期表现出严重的失忆症和运动功能减弱,最终导致死亡[1]。尽管AD的特征性病理改变已明确,但靶向这些病理蛋白的药物疗效欠佳,因此进一步筛查有效靶点非常必要。本研究基于生物信息学分析筛选影响AD发生的核心基因和相关信号通路,现报告如下。
1. 数据与方法
1. 数据与方法
1.1 数据来源
以“Alzheimer′s disease”为关键词搜索基因表达综合(GEO)数据库,并下载GSE227221[2]和GSE162873[3]这2个RNA-seq原始数据集。其中, GSE227221数据集包含217个AD组织样本和226个对照神经组织样本, GSE162873数据集则包含8个AD组织样本。所有样本均从原始资料库下载,且病理资料完整。
1.1 数据来源
以“Alzheimer′s disease”为关键词搜索基因表达综合(GEO)数据库,并下载GSE227221[2]和GSE162873[3]这2个RNA-seq原始数据集。其中, GSE227221数据集包含217个AD组织样本和226个对照神经组织样本, GSE162873数据集则包含8个AD组织样本。所有样本均从原始资料库下载,且病理资料完整。
1.2 生物信息学分析方法
1.2 生物信息学分析方法
1.2.1 差异表达基因(DEGs)的筛选与分析
将GSE227221和GSE162873数据集导入GEO2R在线分析软件(该软件为交互式在线工具,能够比较2组或多组GEO序列数据,明确实验条件下哪些基因发生了显著变化)中,筛选出AD组织样本和正常脑组织样本之间的DEGs, 筛选标准为log2(FC)>1(FC为差异倍数)和P<0.05。
1.2.1 差异表达基因(DEGs)的筛选与分析
将GSE227221和GSE162873数据集导入GEO2R在线分析软件(该软件为交互式在线工具,能够比较2组或多组GEO序列数据,明确实验条件下哪些基因发生了显著变化)中,筛选出AD组织样本和正常脑组织样本之间的DEGs, 筛选标准为log2(FC)>1(FC为差异倍数)和P<0.05。
1.2.2 DEGs的疾病本体论(DO)、基因本体论(GO)和京都基因与基因组百科全书(KEGG)富集分析
使用R软件DOSE包对DEGs进行DO基因富集分析。利用在线数据分析工具DAVID对筛选后的DEGs进行GO功能分析,包括生物学过程、细胞成分和分子功能共3个方面的内容,这有助于理解这些DEGs的生物学功能及其在细胞中的定位。利用在线分析工具DAVID进行KEGG通路富集分析,明确这些DEGs在生物学过程中的角色及其参与的生物学通路。分析结果以柱状图和气泡图形式展示, P<0.05表示差异有统计学意义。
1.2.2 DEGs的疾病本体论(DO)、基因本体论(GO)和京都基因与基因组百科全书(KEGG)富集分析
使用R软件DOSE包对DEGs进行DO基因富集分析。利用在线数据分析工具DAVID对筛选后的DEGs进行GO功能分析,包括生物学过程、细胞成分和分子功能共3个方面的内容,这有助于理解这些DEGs的生物学功能及其在细胞中的定位。利用在线分析工具DAVID进行KEGG通路富集分析,明确这些DEGs在生物学过程中的角色及其参与的生物学通路。分析结果以柱状图和气泡图形式展示, P<0.05表示差异有统计学意义。
1.2.3 枢纽基因的筛选与分析
将筛选出的DEGs导入STRING在线数据库,进行蛋白质交互作用分析。STRING数据库提供了蛋白质之间已知和预测的相互作用信息,这些信息来源于共轭、实验室测试、共表达以及文献记录。本研究利用Cytoscape(版本3.7.2)软件构建蛋白质-蛋白质相互作用(PPI)网络。通过MCODE插件进行枢纽基因筛选,设置参数为关联程度阈值(Degree Cutoff)=2, 节点评分阈值(Node Score Cutoff)=0.2, 最小核心阈值(K-Core)=2, 最大深度(Max Depth)=100。这些参数有助于确定网络内高度相互连接的区域。
1.2.3 枢纽基因的筛选与分析
将筛选出的DEGs导入STRING在线数据库,进行蛋白质交互作用分析。STRING数据库提供了蛋白质之间已知和预测的相互作用信息,这些信息来源于共轭、实验室测试、共表达以及文献记录。本研究利用Cytoscape(版本3.7.2)软件构建蛋白质-蛋白质相互作用(PPI)网络。通过MCODE插件进行枢纽基因筛选,设置参数为关联程度阈值(Degree Cutoff)=2, 节点评分阈值(Node Score Cutoff)=0.2, 最小核心阈值(K-Core)=2, 最大深度(Max Depth)=100。这些参数有助于确定网络内高度相互连接的区域。
1.2.4 核心基因的筛选
根据MCODE评分,对枢纽基因进一步分析并筛选出核心基因。筛选条件为MCODE评分>10, 且评分由高到低排序。选取排名前6位的基因作为核心基因,推测其可能在AD的发生、进展及细胞凋亡过程中发挥关键作用。
1.2.4 核心基因的筛选
根据MCODE评分,对枢纽基因进一步分析并筛选出核心基因。筛选条件为MCODE评分>10, 且评分由高到低排序。选取排名前6位的基因作为核心基因,推测其可能在AD的发生、进展及细胞凋亡过程中发挥关键作用。
2. 结果
2. 结果
2.1 DEGs的筛选与分析
本研究从GSE227221和GSE162873数据集中共筛选出9 705个相关基因,其中GSE227221数据集筛选出7 293个相关基因, GSE162873数据集筛选出3 785个相关基因。将这2个数据集筛选出的基因取交集,共筛选出1 373个DEGs, 这些基因具有相同变化趋势,见图 1。
2.1 DEGs的筛选与分析
本研究从GSE227221和GSE162873数据集中共筛选出9 705个相关基因,其中GSE227221数据集筛选出7 293个相关基因, GSE162873数据集筛选出3 785个相关基因。将这2个数据集筛选出的基因取交集,共筛选出1 373个DEGs, 这些基因具有相同变化趋势,见图 1。
2.2 DO基因富集分析、GO基因富集分析和KEGG通路富集分析结果
DO基因富集分析结果显示,这些DEGs与系统性红斑狼疮、红斑狼疮和动脉硬化这3种疾病最为相关,见图 2。GO功能富集分析结果显示,这些DEGs主要富集于3条信号通路,即经典Wnt信号通路、磷脂酶C-活化G蛋白质-耦合受体通路和等离子体外侧膜通路,见图 3。KEGG信号通路富集分析结果显示,这些DEGs主要富集于以下通路,包括细胞因子-细胞因子受体相互作用、神经元活性配体-受体相互作用、癌症相关的转录失调,见图 4。
2.2 DO基因富集分析、GO基因富集分析和KEGG通路富集分析结果
DO基因富集分析结果显示,这些DEGs与系统性红斑狼疮、红斑狼疮和动脉硬化这3种疾病最为相关,见图 2。GO功能富集分析结果显示,这些DEGs主要富集于3条信号通路,即经典Wnt信号通路、磷脂酶C-活化G蛋白质-耦合受体通路和等离子体外侧膜通路,见图 3。KEGG信号通路富集分析结果显示,这些DEGs主要富集于以下通路,包括细胞因子-细胞因子受体相互作用、神经元活性配体-受体相互作用、癌症相关的转录失调,见图 4。
2.3 枢纽基因的筛选与分析
将DEGs导入STRING数据库,并利用Cytoscape软件构建PPI网络,该网络包含114个节点蛋白和76条边,见图 5。采用MCODE插件分析并挖掘网络中与AD诊断相关的基因,筛选出显著的相互作用模块,称为枢纽基因。本研究筛选出的枢纽基因网络包含61条边及20个节点蛋白,进一步筛选出的共同枢纽基因为 GIMAP1、GIMAP4、GIMAP5、GIMAP6、GIMAP7 及 GIMAP1-GIMAP5 。
2.3 枢纽基因的筛选与分析
将DEGs导入STRING数据库,并利用Cytoscape软件构建PPI网络,该网络包含114个节点蛋白和76条边,见图 5。采用MCODE插件分析并挖掘网络中与AD诊断相关的基因,筛选出显著的相互作用模块,称为枢纽基因。本研究筛选出的枢纽基因网络包含61条边及20个节点蛋白,进一步筛选出的共同枢纽基因为 GIMAP1、GIMAP4、GIMAP5、GIMAP6、GIMAP7 及 GIMAP1-GIMAP5 。
2.4 核心基因的筛选
应用Cytoscape软件中的MCODE插件,筛选条件设定为MCODE评分>10, 并按评分由高至低排序,选取前6个基因作为本研究的核心基因,即 GIMAP1、GIMAP4、GIMAP5、GIMAP6、GIMAP7 及 GIMAP1-GIMAP5 , 见图 6。
2.4 核心基因的筛选
应用Cytoscape软件中的MCODE插件,筛选条件设定为MCODE评分>10, 并按评分由高至低排序,选取前6个基因作为本研究的核心基因,即 GIMAP1、GIMAP4、GIMAP5、GIMAP6、GIMAP7 及 GIMAP1-GIMAP5 , 见图 6。
3. 讨论
痴呆是一种以记忆力、语言能力、执行能力和空间视觉能力进行性下降为特征的临床综合征,最终可导致个性和行为改变,丧失生活自理能力。目前, AD是痴呆的主要病因(占60%~80%)[4],这一现象部分归因于人口老龄化。据估计,全球范围内AD每年影响超过5 000万老年人,同时间接影响数千万AD患者配偶的生活,造成沉重的家庭负担和社会负担[5]。AD的发病机制尚未完全阐明。一小部分AD患者发病与3个基因的显性突变有关,分别为淀粉样前体蛋白(APP)、早老素1(PSEN1)和早老素2(PSEN2), 这些患者通常在65岁前发病[6]。大多数患者呈现散发趋势,发病年龄超过65岁。尽管AD没有明显的遗传性,但已有证据支持存在多种遗传风险因素,例如载脂蛋白E的E4等位基因存在于约16%的AD患者中。针对β-淀粉样蛋白开发的药物仅能缓解症状[7], 因此寻找新的治疗靶点至关重要。
本研究基于生物信息学分析对2个数据集的基因进行比对,共筛选出参与AD发生和进展的1 373个相同趋势的DEGs, 并筛选验证出核心基因GIMAP。人GIMAP基因全称GTP酶免疫相关蛋白,定位于第7号染色体,大小约为500 kb,包含7个功能性基因(GIMAP1、GIMAP2、GIMAP4、GIMAP5、GIMAP6、GIMAP7、GIMAP8)和1个假基因[8]。该基因此前被称为免疫相关核苷酸结合蛋白(IANs), 是一种小的GTP结合分子,具有共同的“AIG”GTP结合域,存在于高等植物和多细胞动物[9]。哺乳动物GIMAP基因(小鼠存在8个亚型,人类存在7个亚型[10])紧密聚集在一个单染色体位置上,这些基因优先表达于造血细胞和淋巴细胞,但近年来有研究[10]提示这些基因亦可表达于非血液细胞。本研究生物信息学分析也发现, GIMAP5 基因在AD神经组织中异常表达,提示AD与血液细胞或免疫细胞之间可能存在某种联系。
早在上世纪50年代,人们就观察到免疫功能异常与痴呆之间可能存在某种联系。1953年, GROSCH H[11]报道1名9岁女孩在接种咳嗽疫苗后,出现皮质下痴呆伴运动功能障碍。然而,随后的70年间,学者们并未发现痴呆或AD与免疫功能之间的确切关联。虽然免疫细胞对神经细胞有吞噬现象,可能对神经细胞造成免疫损伤[12-13],但并未发现这种损伤与AD病理改变如神经纤维缠结及Tau蛋白异常之间的直接联系。
GIMAP蛋白与GTPase在氨基末端序列相似,均包含鸟嘌呤核苷酸结合域[14-15]。这些蛋白大多参与了淋巴细胞的维持与发展。在小鼠模型中, GIMAP5的缺陷导致外周T细胞、B细胞和自然杀伤细胞(NK细胞)的数量减少[16-17]。GIMAP1对T细胞增殖的维持和B细胞功能的成熟至关重要[18-19]。GIMAP4可能促进T细胞的凋亡[20]。敲除GIMAP6使得Jurkat细胞系对凋亡诱导剂变得敏感[21]。鉴于GIMAP基因家族对维持免疫细胞功能的重要性[22], GIMAP基因的异常会导致免疫细胞功能障碍,从而减少免疫细胞对老化及异常神经元的清除和吞噬,增加痴呆或AD的发病风险。
本研究对筛选出的DEGs进行3种类型的富集分析,包括DO基因富集分析、GO功能富集分析和KEGG信号通路富集分析。DO基因富集分析结果显示,这些DEGs主要富集于系统性红斑狼疮、红斑狼疮和动脉硬化这3种疾病,而这3种疾病均与免疫功能异常有关[23-25]。系统性红斑狼疮和红斑狼疮均为自身免疫性疾病,与免疫功能异常关系密切,尽管动脉硬化与免疫功能异常无直接关系,但近年来的研究显示免疫功能异常参与动脉硬化的病理过程。GO功能富集分析结果显示, DEGs主要富集于经典Wnt信号通路、磷脂酶C-活化G蛋白质-耦合受体通路和等离子体外侧膜通路。值得注意的是,这3条信号通路均与GIMAP基因有关,存在一定的调控关系。KEGG信号通路富集分析结果显示,排名前3位的疾病通路分别为细胞因子-细胞因子受体相互作用、神经元活性配体-受体相互作用、癌症相关的转录失调。由此提示, AD与GIMAP基因及这些通路之间存在某种联系。
综上所述, GIMAP基因是AD发病的核心基因,其可能处于AD其他病理改变的上游。鉴于GIMAP基因与免疫功能维持的紧密关系,本研究推测AD是一种由免疫功能紊乱导致的疾病。
3. 讨论
痴呆是一种以记忆力、语言能力、执行能力和空间视觉能力进行性下降为特征的临床综合征,最终可导致个性和行为改变,丧失生活自理能力。目前, AD是痴呆的主要病因(占60%~80%)[4],这一现象部分归因于人口老龄化。据估计,全球范围内AD每年影响超过5 000万老年人,同时间接影响数千万AD患者配偶的生活,造成沉重的家庭负担和社会负担[5]。AD的发病机制尚未完全阐明。一小部分AD患者发病与3个基因的显性突变有关,分别为淀粉样前体蛋白(APP)、早老素1(PSEN1)和早老素2(PSEN2), 这些患者通常在65岁前发病[6]。大多数患者呈现散发趋势,发病年龄超过65岁。尽管AD没有明显的遗传性,但已有证据支持存在多种遗传风险因素,例如载脂蛋白E的E4等位基因存在于约16%的AD患者中。针对β-淀粉样蛋白开发的药物仅能缓解症状[7], 因此寻找新的治疗靶点至关重要。
本研究基于生物信息学分析对2个数据集的基因进行比对,共筛选出参与AD发生和进展的1 373个相同趋势的DEGs, 并筛选验证出核心基因GIMAP。人GIMAP基因全称GTP酶免疫相关蛋白,定位于第7号染色体,大小约为500 kb,包含7个功能性基因(GIMAP1、GIMAP2、GIMAP4、GIMAP5、GIMAP6、GIMAP7、GIMAP8)和1个假基因[8]。该基因此前被称为免疫相关核苷酸结合蛋白(IANs), 是一种小的GTP结合分子,具有共同的“AIG”GTP结合域,存在于高等植物和多细胞动物[9]。哺乳动物GIMAP基因(小鼠存在8个亚型,人类存在7个亚型[10])紧密聚集在一个单染色体位置上,这些基因优先表达于造血细胞和淋巴细胞,但近年来有研究[10]提示这些基因亦可表达于非血液细胞。本研究生物信息学分析也发现, GIMAP5 基因在AD神经组织中异常表达,提示AD与血液细胞或免疫细胞之间可能存在某种联系。
早在上世纪50年代,人们就观察到免疫功能异常与痴呆之间可能存在某种联系。1953年, GROSCH H[11]报道1名9岁女孩在接种咳嗽疫苗后,出现皮质下痴呆伴运动功能障碍。然而,随后的70年间,学者们并未发现痴呆或AD与免疫功能之间的确切关联。虽然免疫细胞对神经细胞有吞噬现象,可能对神经细胞造成免疫损伤[12-13],但并未发现这种损伤与AD病理改变如神经纤维缠结及Tau蛋白异常之间的直接联系。
GIMAP蛋白与GTPase在氨基末端序列相似,均包含鸟嘌呤核苷酸结合域[14-15]。这些蛋白大多参与了淋巴细胞的维持与发展。在小鼠模型中, GIMAP5的缺陷导致外周T细胞、B细胞和自然杀伤细胞(NK细胞)的数量减少[16-17]。GIMAP1对T细胞增殖的维持和B细胞功能的成熟至关重要[18-19]。GIMAP4可能促进T细胞的凋亡[20]。敲除GIMAP6使得Jurkat细胞系对凋亡诱导剂变得敏感[21]。鉴于GIMAP基因家族对维持免疫细胞功能的重要性[22], GIMAP基因的异常会导致免疫细胞功能障碍,从而减少免疫细胞对老化及异常神经元的清除和吞噬,增加痴呆或AD的发病风险。
本研究对筛选出的DEGs进行3种类型的富集分析,包括DO基因富集分析、GO功能富集分析和KEGG信号通路富集分析。DO基因富集分析结果显示,这些DEGs主要富集于系统性红斑狼疮、红斑狼疮和动脉硬化这3种疾病,而这3种疾病均与免疫功能异常有关[23-25]。系统性红斑狼疮和红斑狼疮均为自身免疫性疾病,与免疫功能异常关系密切,尽管动脉硬化与免疫功能异常无直接关系,但近年来的研究显示免疫功能异常参与动脉硬化的病理过程。GO功能富集分析结果显示, DEGs主要富集于经典Wnt信号通路、磷脂酶C-活化G蛋白质-耦合受体通路和等离子体外侧膜通路。值得注意的是,这3条信号通路均与GIMAP基因有关,存在一定的调控关系。KEGG信号通路富集分析结果显示,排名前3位的疾病通路分别为细胞因子-细胞因子受体相互作用、神经元活性配体-受体相互作用、癌症相关的转录失调。由此提示, AD与GIMAP基因及这些通路之间存在某种联系。
综上所述, GIMAP基因是AD发病的核心基因,其可能处于AD其他病理改变的上游。鉴于GIMAP基因与免疫功能维持的紧密关系,本研究推测AD是一种由免疫功能紊乱导致的疾病。
-
表 1 海藻-昆布的活性成分
分子编码 分子名称 分子质量 OB/% DL MOL010615 海藻甾醇 426.31 43.48 0.62 MOL010616 鹅掌菜酚 372.30 87.06 0.63 MOL010617 1553-41-9 302.50 45.66 0.21 MOL010625 24-二四亚甲基胆固醇 412.77 43.54 0.76 MOL001439 花生四烯酸 304.52 45.57 0.20 MOL000953 CLR 386.73 37.87 0.68 MOL009622 岩藻甾醇 412.77 43.78 0.76 MOL010578 N-[(1S)-1-(benzyl)-2-[[(1S)-1-(benzyl)-2-hydroxy-ethyl]amino]-2-keto-ethyl]benzamide 402.53 45.76 0.43 MOL010580 Diglycol dibenzoate 314.36 59.22 0.27 MOL005440 异岩藻甾醇 412.77 43.78 0.76 MOL000098 槲皮素 302.25 46.43 0.28 表 2 分子对接结合能
分子名称 基因 结合能 槲皮素 MYC -8.1 海藻甾醇 TP53 -7.2 海藻甾醇 CASP3 -6.9 槲皮素 EGFR -6.8 CLR TNF -6.7 岩藻甾醇 VEGFA -6.6 鹅掌菜酚 JUN -6.4 N-[(1S)-1-(benzyl)-2-[[(1S)-1-(benzyl)-2-hydroxy-ethyl]amino]-2-keto-ethyl]benzamide IL1B -6.3 异岩藻甾醇 AKT1 -6.1 海藻甾醇 IL6 -6.1 -
[1] SUNG H, FERLAY J, SIEGEL R L, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin, 2021, 71(3): 209-249. doi: 10.3322/caac.21660
[2] BARZAMAN K, KARAMI J, ZAREI Z, et al. Breast cancer: biology, biomarkers, and treatments[J]. Int Immunopharmacol, 2020, 84: 106535.
[3] 韦明婵, 林江, 易劲苍, 等. 含海洋中药的抗肿瘤方剂用药规律分析[J]. 北京中医药大学学报, 2018, 41(3): 253-258. doi: 10.3969/j.issn.1006-2157.2018.03.013 [4] HSIAO H H, WU T C, TSAI Y H, et al. Effect of oversulfation on the composition, structure, and in vitro anti-lung cancer activity of fucoidans extracted from Sargassum aquifolium[J]. Mar Drugs, 2021, 19(4): 215.
[5] TSOU M H, LEE C C, WU Z Y, et al. Bioactivity of crude fucoidan extracted from Sargassum ilicifolium (Turner) C. Agardh[J]. Sci Rep, 2022, 12(1): 15916.
[6] PRADHAN B, PATRA S, NAYAK R, et al. Multifunctional role of fucoidan, sulfated polysaccharides in human health and disease: a journey under the sea in pursuit of potent therapeutic agents[J]. Int J Biol Macromol, 2020, 164: 4263-4278.
[7] WANG X, WANG Z Y, ZHENG J H, et al. TCM network pharmacology: a new trend towards combining computational, experimental and clinical approaches[J]. Chin J Nat Med, 2021, 19(1): 1-11.
[8] NAMVAR F, MOHAMAD R, BAHARARA J, et al. Antioxidant, antiproliferative, and antiangiogenesis effects of polyphenol-rich seaweed (Sargassum muticum)[J]. Biomed Res Int, 2013, 2013: 604787.
[9] SHARMA A, KASHYAP D, SAK K, et al. Therapeutic charm of quercetin and its derivatives: a review of research and patents[J]. Pharm Pat Anal, 2018, 7(1): 15-32.
[10] YARLA N S, BISHAYEE A, SETHI G, et al. Targeting arachidonic acid pathway by natural products for cancer prevention and therapy[J]. Semin Cancer Biol, 2016, 40/41: 48-81.
[11] MONTEIRO P, LOMARTIRE S, COTAS J, et al. Call the eckols: present and future potential cancer therapies[J]. Mar Drugs, 2022, 20(6): 387.
[12] LI X, ZHOU D Y, LI F T, et al. Saringosterol acetate isolated from Sargassum fusiformis induces mitochondrial-mediated apoptosis in MCF-7 breast cancer cells[J]. Chem Biodivers, 2022, 19(3): e202100848.
[13] 杨玲, 郑超群, 吴福珍, 等. 岩藻甾醇对卵巢癌细胞抑制作用的体外研究[J]. 北方药学, 2020, 17(5): 145-146, 196. [14] MAO Z F, SHEN X L, DONG P, et al. Fucosterol exerts antiproliferative effects on human lung cancer cells by inducing apoptosis, cell cycle arrest and targeting of Raf/MEK/ERK signalling pathway[J]. Phytomedicine, 2019, 61: 152809.
[15] XU H, LIN F J, WANG Z L, et al. CXCR2 promotes breast cancer metastasis and chemoresistance via suppression of AKT1 and activation of COX2[J]. Cancer Lett, 2018, 412: 69-80.
[16] LI C W, XIA W Y, LIM S O, et al. AKT1 inhibits epithelial-to-mesenchymal transition in breast cancer through phosphorylation-dependent Twist1 degradation[J]. Cancer Res, 2016, 76(6): 1451-1462.
[17] JEOUNG D I, TANG B, SONENBERG M. Effects of tumor necrosis factor-alpha on antimitogenicity and cell cycle-related proteins in MCF-7 cells[J]. J Biol Chem, 1995, 270(31): 18367-18373.
[18] SLATTERY M L, HERRICK J S, TORRES-MEJIA G, et al. Genetic variants in interleukin genes are associated with breast cancer risk and survival in a genetically admixed population: the Breast Cancer Health Disparities Study[J]. Carcinogenesis, 2014, 35(8): 1750-1759.
[19] JOHNSON D E, O'KEEFE R A, GRANDIS J R. Targeting the IL-6/JAK/STAT3 signalling axis in cancer[J]. Nat Rev Clin Oncol, 2018, 15(4): 234-248.
[20] KAJAL K, BOSE S, PANDA A K, et al. Transcriptional regulation of VEGFA expression in T-regulatory cells from breast cancer patients[J]. Cancer Immunol Immunother, 2021, 70(7): 1877-1891.
[21] QIAO Y C, HE H, JONSSON P, et al. AP-1 is a key regulator of proinflammatory cytokine TNFα-mediated triple-negative breast cancer progression[J]. J Biol Chem, 2016, 291(35): 18309.
[22] PU X, STORR S J, ZHANG Y M, et al. Caspase-3 and caspase-8 expression in breast cancer: caspase-3 is associated with survival[J]. Apoptosis, 2017, 22(3): 357-368.
[23] LIU Y Q, ZHU C J, TANG L, et al. MYC dysfunction modulates stemness and tumorigenesis in breast cancer[J]. Int J Biol Sci, 2021, 17(1): 178-187.
[24] ZHAO L N, QIU T, JIANG D W, et al. SGCE promotes breast cancer stem cells by stabilizing EGFR[J]. Adv Sci, 2020, 7(14): 1903700.
[25] LI K, ZHANG T T, ZHAO C X, et al. Faciogenital Dysplasia 5 supports cancer stem cell traits in basal-like breast cancer by enhancing EGFR stability[J]. Sci Transl Med, 2021, 13(586): eabb2914.
[26] SONG X, LIU Z Y, YU Z Y. EGFR promotes the development of triple negative breast cancer through JAK/STAT3 signaling[J]. Cancer Manag Res, 2020, 12: 703-717.
[27] ZHAO Y Y, MA J, FAN Y L, et al. TGF-β transactivates EGFR and facilitates breast cancer migration and invasion through canonical Smad3 and ERK/Sp1 signaling pathways[J]. Mol Oncol, 2018, 12(3): 305-321.