目的基于生物信息学分析和网络药理学探究莫诺拉韦在合并乳腺癌和新型冠状病毒感染(BRCA/COVID-19)患者中的治疗靶点和作用机制。
方法从美国癌症基因组图谱计划(TCGA)数据库检索了乳腺浸润癌(BRCA)患者的RNA测序数据和临床信息,并从公共数据库获取相关的新型冠状病毒感染(COVID-19)靶点。利用单因素、多因素Cox回归分析构建BRCA/COVID-19患者的预后风险模型。通过多个指标的ROC曲线评估模型的可靠性和敏感度。根据中位风险评分将患者分为高、低风险组,分析其生存状态、分布和临床相关性。莫诺拉韦的靶点来自药物数据库,筛选出治疗BRCA/COVID-19的潜在靶基因集。生物信息学分析包括蛋白质-蛋白质相互作用网络分析(PPI)和基因本体(GO)/京都基因与基因组百科全书(KEGG)富集分析,探究中心基因和信号通路。分子对接技术分析莫诺拉韦与关键靶点基因之间的结合亲和力。
结果生物信息学方法确定BRCA/COVID-19患者筛选的潜在生物标志物。网络药理学方法揭示5个核心靶点(CDK1、EGFR、NCPAG、TOP2A和IL6)。富集分析表明,除了其抗病毒效果外,莫诺拉韦还影响细胞因子介导的信号通路、白细胞介素(IL)-17信号通路、Janus激酶信号转导和转录激活因子(JAK-STAT)信号通路、免疫调节和炎症调控。分子对接显示了莫诺拉韦与核心靶点以及主要COVID-19蛋白酶之间的良好结合亲和力。
结论本研究揭示了BRCA/COVID-19患者的预后标志物和临床特征,确定了莫诺拉韦治疗的潜在药理靶点和作用机制,并为这一领域的基础研究提供宝贵的见解。
ObjectiveTo explore the therapeutic targets and mechanisms of Molnupiravir in patients with combined breast cancer and COVID-19 (BRCA/COVID-19) using bioinformatics analysis and network pharmacology.
MethodsRNA sequencing data and clinical information of BRCA patients were retrieved from The Cancer Genome Atlas Program (TCGA) database, and relevant COVID-19 targets were obtained from public databases. A prognostic risk model for BRCA/COVID-19 patients was constructed using single-factor and multi-factor Cox re-gression analysis. The reliability and sensitivity of the model were evaluated through ROC curve analysis of multi-ple indicators. Patients were stratified into high- and low-risk groups based on median risk scores, and their sur-vival status, distribution, and clinical relevance were analyzed. Potential target genes for treating BRCA/COVID-19 were identified from drug databases. Bioinformatics analysis included protein-protein interaction (PPI) network analysis and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to ex-plore central genes and signaling pathways. Molecular docking analysis was used to analyze the binding affinity between Molnupiravir and key target genes.
ResultsBioinformatics methods identified potential biomarkers screened for BRCA/COVID-19 patients. Network pharmacology methods revealed 5 core targets (CDK1, EGFR, NCAPG, TOP2A, and IL6). Enrichment analysis indicated that besides its antiviral effects, Molnupiravir also af-fects cytokine-mediated signaling pathways, the IL-17 signaling pathway, JAK-STAT signaling pathway, immune regulation, and inflammation control. Molecular docking showed good binding affinity between Molnupiravir and core targets as well as the main COVID-19 protease.
ConclusionThis study revealed prognostic markers and clinical characteristics of BRCA/COVID-19 patients, identified potential pharmacological targets and mechanisms of Molnupiravir treatment, and provided valuable insights for basic research in this field.
叶若雷,苏燕萍,华建伟,等. 乳腺癌/新型冠状病毒感染患者中莫诺拉韦的潜在治疗靶点和机制:基于生物信息学和网络药理学分析[J]. 中国药物与临床,2025,25(06):387-397.
DOI:10.11655/zgywylc2025.06.009
注:BRCA,乳腺浸润癌;COVID-19,新型冠状病毒感染;a.BRCA交集基因的韦恩图;b.BRCA/COVID-19交集基因的韦恩图;c.差异表达基因的火山图。
基因 | K-M P值 | HR值(95% CI) | Cox P值 |
---|---|---|---|
SDC1 | 0.002 19 | 1.15(1.03,1.30) | 0.000 77 |
ATP6AP1 | 0.023 509 | 1.25(1.11,1.45) | 0.028 892 |
MRPL13 | 0.005 502 | 1.28(1.09,1.57) | 9.52×10 -5 |
ANO6 | 0.015 546 | 1.35(1.13,1.38) | 0.005 14 |
SLC35A2 | 0.001 036 | 1.15(1.04,1.19) | 0.000 579 |
LIMCH1 | 0.008 73 | 1.15(1.11,1.20) | 0.000 195 |
DCAF13 | 0.025 678 | 1.05(1.00,1.12) | 0.002 106 |
COMTD1 | 0.024 968 | 1.31(1.04,1.37) | 0.001 326 |
C2CD2 | 0.005 922 | 1.14(1.05,1.23) | 0.044 919 |
PHLDB2 | 0.039 036 | 1.26(1.12,1.40) | 0.025 387 |
PAX7 | 0.011 96 | 1.18(1.10,1.21) | 0.000 691 |
EPB41L4B | 0.005 672 | 1.29(1.12,1.44) | 4.35×10 -5 |
CENPI | 0.033 095 | 1.03(1.00,1.06) | 0.030 022 |
IL27 | 0.024 297 | 1.24(1.02,1.50) | 0.029 |
基因 | 回归系数 | HR值(95% CI) | Z值 | P值 |
---|---|---|---|---|
ANO6 | 0.289 59 | 1.34(1.02,1.76) | 2.069 | 0.038 |
ATP6AP1 | 0.242 26 | 1.27(0.97,1.67) | 1.755 | 0.079 |
COMTD1 | 0.257 62 | 1.29(1.06,1.58) | 2.49 | 0.012 |
EPB41L4B | 0.239 38 | 1.27(1.04,1.55) | 2.369 | 0.017 |
MRPL13 | 0.337 02 | 1.40(1.12,1.75) | 2.95 | 0.003 |
PAX7 | 0.132 58 | 1.14(1.01,1.29) | 2.182 | 0.029 |
SDC1 | 0.120 72 | 1.13(0.97,1.32) | 1.532 | 0.125 |
注:BRCA,乳腺浸润癌;COVID-19,新型冠状病毒感染;a.单因素Cox分析森林图;b.高、低风险组患者OS的K-M曲线;c.BRCA患者高低风险分布图;d.BRCA患者风险评分生存图;e.1年、3年、5年内风险评分和临床特征的ROC曲线。
注:a.单因素Cox分析;b.多因素Cox分析。
注:GO,基因本体富集分析;KEGG,京都基因与基因组百科全书富集分析。a.交叉基因的韦恩图;b.GO富集分析;c.KEGG富集分析。
基因 | 度值 | 介数中心性 | 邻域组件中心性 | 紧密中心性 | 核心网络 |
---|---|---|---|---|---|
EGFR | 49 | 34 552 | 44 | 189.383 3 | / |
IL6 | 46 | 14 926 | 44 | 176.416 7 | / |
CDK1 | 42 | 18 901 | 40 | 163.516 7 | MCODE Cluster1 |
TOP2A | 37 | 5 711 | 37 | 151.9 | MCODE Cluster1 |
NCAPG | 35 | 761 | 35 | 139.442 9 | MCODE Cluster1 |
注:PPI,蛋白质-蛋白质相互作用网络分析。
注:PPI,蛋白质-蛋白质相互作用网络分析;BRCA,乳腺浸润癌;COVID-19,新型冠状病毒感染;a.PPI网络图;b~d.基于MCODE插件进行聚类的PPI网络;e.基于cytoHubba和MCODE插件综合分析结果筛选出的5个中心基因。
注:BRCA,乳腺浸润癌;COVID-19,新型冠状病毒感染。
注:BRCA,乳腺浸润癌;COVID-19,新型冠状病毒感染。a.莫诺拉韦与COVID-19宿主蛋白5r84对接;b~f.莫诺拉韦与BRCA/COVID-19核心靶点对接。
靶点 | PDB编号 | 中心坐标( x, y, z) | 亲和力/(kCal/mol) |
---|---|---|---|
COVID-19 | 5R84 | 9.8,-0.3,20.8 | -6.1 |
EGFR | 1M17 | 21.7,0.3,52.1 | -6.8 |
IL-6 | 4CNI | 68.8,-50,-5.9 | -6.4 |
CDK1 | 6GU2 | 328.6,213.9,192.3 | -6.7 |
TOP2A | 4R1F | 49.1,-2.6,82.2 | -8.3 |
NCAPG | 6IGX | 356.7,37.4,433.4 | -5.6 |
叶若雷、徐艳艳:实验研究设计;苏燕萍:参与文章撰写;华建伟:采集处理数据;骆松梅:对文章的知识性内容作批评性审阅

你好,我可以帮助您更好的了解本文,请向我提问您关注的问题。