目的 通过生物信息学和网络药理学方法探讨红景天苷治疗三阴性乳腺癌（TNBC）的作用机制，阐明其产生治疗作用的主要靶点和信号通路。 方法 通过基因表达综合数据库（GEO）获取数据集 GSE45827，利用R软件包GSEABase进行基因集富集分析（GSEA），采用limma R软件包寻找相邻正常组织和TNBC组织之间的差异表达基因（DEGs），对DEGs进行基因本体论（GO）功能富集分析和京都基因与基因组百科全书（KEGG）信号通路富集分析，将DEGs与药物靶点结合，导入基因/蛋白相互作用检索搜查工具String数据库，形成蛋白-蛋白相互作用 （PPI）网络。使用MCODE插件对 PPI网络进行功能模块筛选，对SCORE值排名前2位的关键模块基因再次进行GO功能富集分析和KEGG信号通路富集分析。将2次KEGG富集分析所得通路与转录组数据GSEA富集分析结果取交集，获得红景天苷治疗TNBC的作用通路。使用CytoHubba插件计算出关键模块中最大团中心性（MCC）评分前 10 位的关键节点基因， 即为核心基因。 应用 AutoDock Vina 1.1.2 和PyMOL 2.3.0软件完成分子对接。 结果 KEGG与GSEA富集分析的结果取交集得到13条共同通路，涉及细胞周期、细胞衰老和p53信号通路等。GO功能富集分析结果中所涉及的有丝分裂、核分裂和姐妹染色单体分离等生物学过程与细胞周期有密切关联，与KEGG富集分析结果一致。SCORE值排名第1位的关键模块中包含5个红景天苷药物作用靶点，分别为重组人细胞周期蛋白A2（CCNA2）、细胞周期检查点激酶1（CHEK1）、驱动蛋白家族成员11（KIF11）、DNA拓扑异构酶2（TOP2A）和胸腺嘧啶酸合酶（TYMS），将上述蛋白与红景天苷进行分子对接，结果均表现出很强的结合能力（结合能<－7.0 kcal·mol^－1）。 结论 红景天苷的紧密结合靶标位于TNBC的DEGs关键功能模块中，可以与CCNA2蛋白结合产生直接的调控作用，与KIF11、TOPA2、CHEK1和TYMS蛋白结合可针对TNBC的关键节点基因产生间接的调控作用，红景天苷有可能成为TNBC的临床治疗药物。

Objective To discuss the mechanism of salidroside in the treatment of triple negative breast cancer （TNBC） by using the bioinformatics and network pharmacology methods， and to clarify the main targets and signaling pathways involved in the therapeutic effect. Methods The dataset GSE45827 was obtained from the Gene Expression Omnibus （GEO） database； the gene set enrichment analysis （GSEA） was performed by using the R software package GSEABase；the differentially expressed genes （DEGs） between the adjacent normal tissue and TNBC tissue were identified by limma R software package；the Gene Ontology （GO） functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes （KEGG） signaling pathway enrichment analysis were performed on the DEGs， and the DEGs were integrated with the drug targets to import into gene/protein interaction retrieval tool String database， and the protein-protein interaction （PPI） networks were constructed；the functional module screening of the PPI network was conducted by MCODE plugin， and the top 2 modules ranked by SCORE value were further subjected to GO functional enrichment analysis and KEGG signaling pathway analysis. The pathways obtained from the two rounds of KEGG enrichment analysis were intersected with the results of GSEA enrichment analysis to identify the pathways involved in the therapeutic effect of salidroside on TNBC. The top 10 key node genes in the highest scoring module determined by the maximum clique centrality （MCC） score caculated by CytoHubba plugi were considered as the core genes； the molecular docking was performed by AutoDock Vina1.1.2 and PyMOL2.3.0 Software. Results The intersection of KEGG and GSEA enrichment analysis results showed 13 singaling pathways， including the cell cycle， cellular senescence， and p53 signaling pathways，and so on. The biological processes involved in the GO functional analysis， such as mitosis， nuclear division， and sister chromatid separation， were closely related to the cell cycle and consistented with the results of the KEGG signaling pathway enrichment analysis. The top ranked module based on the SCORE value contained 5 drug target genes of Rhodiola glycoside，such as cyclin A2 （CCNA2）， checkpoint kinase 1 （CHEK1）， kinesin family member 11 （KIF11）， DNA topoisomerase 2-alpha （TOP2A）， and thymidylate synthase （TYMS）. The molecular docking results demonstrated strong binding affinities between the above proteins and Rhodiola glycoside （binding energy<－7.0 kcal·mol^－1）. Conclusion The tightly binding target of salidroside is located in the key functional modules of DEGs of TNBC， which can directly regulate by binding with CCNA2 and protein， and indirectly regulate the key differentially genes of TNBC by binding with KIF11， TOPA2， CHEK1 and TYMS proteins. Therefore， salidroside may be a potential clinical therapeutic drug for TNBC.