Exploration of Therapeutic Drugs for Gastric Cancer Using Drug Repositioning Strategy

在NCBI的GEO数据库中，下载与胃癌相关的3组基因表达谱数据集GSE54129、GSE65801和GSE26899。其中，GSE54129数据集的平台为GPL570，包含111例胃癌组织样本和21例正常对照组织样本；GSE65801数据集的平台为GPL14550，包含32例胃癌组织样本和32例配对的正常对照组织样本；GSE26899数据集的平台为GPL6947，包含96例胃癌组织样本和12例正常对照组织样本。

使用“Limma”R包鉴定在胃癌组织和正常对照组织样本之间显著差异表达的基因^[5]，鉴定标准为|log₂(fold change)|≥1且校正后的P<0.05。进一步，通过对差异表达基因(differentially expressed genes, DEGs)绘制韦恩图，筛选在上述3个数据集中均上调表达或下调表达的基因。

利用“clusterProfiler”R包^[6]，对在3个数据集中重叠的差异表达基因分别进行基因本体论(gene ontology, GO)功能富集分析及京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)通路富集分析，探究差异表达基因显著富集的条目(校正后的P值小于0.05)。

将在GSE54129、GSE65801和GSE26899这3个数据集中均显著上调的基因和均显著下调基因上传到CMap数据库中的“Query”模块进行分析^[7]。结果文件中，连通度分数(connectivity score)表示胃癌相关基因的表达特征与CMap中收录的小分子化合物之间的相关性。连通度分数越接近−100，表示药物的治疗潜力越大。本研究选择连通度分数小于−80作为筛选治疗胃癌候选小分子化合物的阈值。

利用STRING数据库对3个数据集中重叠差异表达基因之间的相互作用关系进行分析，置信分数(confidence score)代表蛋白质相互作用的强弱^[8]。本研究以置信分数大于0.400为阈值筛选可靠的蛋白质相互作用关系，并构建PPIN。进一步，使用Cytoscape软件(3.6.1版本)^[9]进行PPIN的可视化和拓扑性质分析，鉴定连接度(degree)值排名前10%的结点为网络中的枢纽结点。

为了探究候选药物的作用靶点，利用AutoDock软件(4.2.6版本)^[10]对PPIN中枢纽基因编码的蛋白质与候选药物进行分子对接。具体地，在PDB数据库( https://www.rcsb.org/)^[11]中分别检索和下载PPIN中枢纽基因编码蛋白质的三维结构文件，同时在PubChem数据库^[12]中搜索和下载候选药物的三维结构文件。进一步，在对蛋白质和药物配体结构进行预处理之后，利用AutoDock软件对候选药物和蛋白质进行分子对接，每一个蛋白质与相应药物要进行50次对接，选出结合能最低的对接结果。当结合能小于−5.0 kcal/mol时，表明配体和受体具有较好的结合活性。使用Pymol软件对亲和力最佳的对接结果进行可视化和分析^[13]。

为了探究潜在药物靶点的预后价值，对其与胃癌患者生存时间的相关性进行分析。具体地，使用在线工具Kaplan-Meier plotter^[14]，分析靶点基因高表达组和低表达组之间生存时间(包括总体生存期和无进展生存期)是否存在显著的差异。统计学上，使用对数秩检验(P<0.05)来评估两组间生存时间差异的显著性。

使用“Limma”R软件包分别对3个数据集GSE54129、GSE26899和GSE65801中肿瘤组织和正常对照组织样本间差异表达的基因进行鉴定。结果显示，以|log₂(fold change)|≥1和校正后的P<0.05为筛选标准，在GSE54129、GSE26899和GSE65801数据集中分别鉴定了2 533、521和2 403个差异表达基因(图1a～图1c)。值得注意的是，76个上调基因和127个下调基因在3个数据集中均显著差异表达(图1d～图1e)，表明这些基因在胃癌致病过程中可能发挥着重要的作用。

为了进一步研究上述3个胃癌数据集中重叠差异表达基因的功能，使用“clusterProfiler”R包分别进行了GO功能富集分析和KEGG代谢通路富集分析。如图2a所示，重叠的差异表达基因在细胞外基质组织(extracellular matrix organization)、组织稳态(tissue homeostasis)和激素代谢过程(hormone metabolic process)等生物过程条目显著富集，在含胶原蛋白的胞外基质(collagen-containing extracellular matrix)和内质网腔(endoplasmic reticulum lumen)等细胞组分条目显著富集，在细胞外基质结构成分(extracellular matrix structural constituent)、氧化还原酶活性(oxidoreductase activity)等分子功能条目显著富集。KEGG结果也显示重叠的差异表达基因在蛋白质的消化吸收(Protein digestion and absorption)、药物代谢−细胞色素P450通路(Drug metabolism-cytochrome P450)和细胞外基质−受体相互作用(ECM-receptor interaction)等与胃癌的发生和发展密切相关的代谢通路显著富集，如图2b所示。

为了探索治疗胃癌的候选小分子药物，本研究采用CMap在线分析工具对胃癌相关基因表达特征与小分子药物之间的联系进行了分析。具体地，通过将重叠的76个上调基因和127个下调基因上传到CMap数据库，以连接分数<−80分为标准，鉴定了10种治疗胃癌的候选药物，如表1所示。

其中，PD-0325901、菲达替尼(fedratinib，TG-101348)、NVP-AUY922、西地尼布和布雷非德菌素A(brefeldin A)为已报道的具有胃癌治疗价值的小分子药物，表明本实验的结果具有可靠性。另外，左炔诺孕酮、TWS-119、Tyrphostin-AG-1478、VU-0415374-1和L-690488为本研究新发现的治疗胃癌的候选小分子化合物。

蛋白质间的相互作用是细胞内蛋白质执行生命功能的重要方式，在包括癌症发生发展等多种生命活动中发挥着至关重要的作用。本研究利用STRING数据库，对在上述3个数据集中鉴定的203个差异表达基因进行相互作用分析，构建了由152个节点和480条边组成的PPIN，如图3a所示。鉴于PPIN中枢纽基因具有成为药物靶点的潜力，利用Cytoscape软件对构建的PPI网络进行可视化和拓扑性质分析，以度(degree)值排名前10%为筛选标准，鉴定了16个枢纽基因，如图3b所示。其中，TIMP1为基质金属蛋白酶组织抑制剂(tissue inhibitor of matrix metalloproteinase, TIMP)家族重要的成员之一，在抑制基质金属蛋白酶和细胞信号通路(如细胞死亡、增殖和血管生成)等方面具有重要作用^[15]；PDGFRB是血小板衍生生长因子受体(platelet-derived growth factor receptor, PDGFR)，在抑制血管生成和肿瘤细胞增殖过程中起关键作用^[16]。而且，TIMP1和PDGFR均可作为胃癌预后生物标志物^[16-17]，提示枢纽基因在胃癌发病过程中可能发挥着关键作用。

为了预测胃癌候选小分子药物的靶点，本研究利用Autodock软件对PPIN中枢纽基因编码蛋白质与基于CMap数据库鉴定的候选药物进行了分子对接分析。结果显示，左炔诺孕酮和西地尼布可分别与TIMP1和PDGFRB结合。具体地，西地尼布可与PDGFRB蛋白的口袋结合，结合位点分别为第63位和第112位谷氨酸(GLU)，结合能为−6.28 kcal/mol，如图4a所示，这与文献[18]报道的西地尼布能够结合PDGFR家族成员相符合；左炔诺孕酮与TIMP1第100位的丝氨酸(SER)、第102位的缬氨酸(VAL)、第112位的谷氨酰胺(GLN)和第162位的精氨酸(ARG)作用，结合能为−9.24 kcal/mol，如图4b所示。因此，本研究新鉴定的胃癌候选药物左炔诺孕酮可能通过结合TIMP1的口袋，达到治疗胃癌的效果。

首先，为了验证候选药物潜在靶点基因在肿瘤组织和正常对照组织中的表达模式，通过箱线图展示了TIMP1和PDGFRB在GSE54129、GSE26899和GSE65801 这3个数据集中的表达模式，并利用GEPIA2在线分析工具^[19]对药物潜在靶点TIMP1和PDGFRB在TCGA数据库中的表达水平进行了分析。如图5a~5b所示，靶点基因TIMP1和PDGFRB在胃癌肿瘤组织中的表达水平均显著上调(P<0.05)。然后，通过ROC曲线分析评估TIMP1和PDGFRB来区分胃癌患者和正常对照的敏感性和特异性，AUC值越高则代表区分能力越强。如图5c所示，在GSE54129、GSE26899和GSE65801数据集中，TIMP1的AUC值在0.960～0.994之间，基因PDGFRB的AUC值在0.786～0.974之间，表明TIMP1和PDGFRB对胃癌具有较好的诊断效果。进一步，利用Kaplan-Meier Plotter工具进行生存分析来研究TIMP1和PDGFRB分别与胃癌患者总体生存时间的相关性。结果显示，与TIMP1低表达组相比，TIMP1高表达组患者的总体生存率显著降低(P<0.05)；类似地，PDGFRB高表达组患者的总体生存率也显著低于PDGFRB低表达组患者的总体生存率(P<0.05)，如图5d所示，表明靶点基因TIMP1和PDGFRB的表达水平与胃癌患者的预后生存时间密切相关。

胃癌是世界范围内最常见的癌症类型之一，大多数患者到医院检查就诊时已发展到中晚期，手术治疗效果不佳，现有药物的治疗效果有限，因此，迫切需要研发具有更好疗效的靶向药物，提高胃癌的临床治疗效果。

本研究基于胃癌的基因表达谱数据，采用CMap的药物重定位方法，鉴定了10种治疗胃癌的候选小分子化合物。其中，PD-0325901、菲达替尼、NVP-AUY922、西地尼布和布雷非德菌素A为已报道的具有胃癌治疗潜力的小分子化合物。如PD-0325901作为MEK抑制剂能够抑制胃癌细胞的生长^[20]；西地尼布为KIT和VEGFR抑制剂，具有抗血管生成的作用，文献[21]对西地尼布联合顺铂和氟嘧啶治疗日本初治晚期胃癌患者进行了临床I期研究，发现20毫克/天西地尼布联合顺铂和氟嘧啶是可以耐受的，具有抗肿瘤的活性。此外，左炔诺孕酮、TWS-119、tyrphostin-AG-1478、VU-0415374-1和L-690488为本研究新发现的治疗胃癌的候选小分子化合物，值得深入研究。

近年来，越来越多的研究表明，PPIN是探究癌症致病机制和预测候选药物靶标的重要工具。通过鉴定PPIN中的枢纽节点并将其分别与上述筛选得到的小分子药物进行分子对接，发现小分子药物西地尼布可与PDGFRB结合，结合能为−6.28 kcal/mol，这与文献[18]报道的西地尼布能够结合PDGFR家族成员的结果是一致的。另外，分子对接结果也显示，左炔诺孕酮可与TIMP1的口袋结合，结合能为−9.24 kcal/mol。但是左炔诺孕酮对胃癌的治疗作用尚未报道，为本研究的新发现。

左炔诺孕酮是一种口服避孕药，主要通过抑制排卵、改变子宫内膜和宫颈黏液，达到避孕的效果。近年来，越来越多的研究表明，左炔诺孕酮或左炔诺孕酮宫内节育系统(levonorgestrel intrauterine system，LNG-IUS)还具有抗癌作用，尤其对于子宫内膜癌治疗有一定的疗效，但是它们对癌症的影响仍存在争议。如左炔诺孕酮宫内节育系统可通过降低异位子宫内膜中腺体和间质雌激素受体α(estrogen receptor α, ERα)、ERβ和孕激素受体(progesterone receptor, PR)的表达，改善轻度至中度子宫内膜异位症女性的症状^[22]。文献[23]发现左炔诺孕酮宫内节育系统能够影响宫腺肌症患者子宫内膜中类固醇受体共调节因子的表达水平。文献[24]的研究结果显示，左炔诺孕酮能够通过雄激素受体抑制人子宫内膜内皮细胞中组织纤溶酶原激活物的分泌。另外，文献[25]的研究结果显示，孕酮宫内节育系统与低于预期的子宫内膜癌、卵巢癌、胰腺癌和肺癌发病率相关，但与高于预期的乳腺癌发病率相关^[25]。需要注意的是，左炔诺孕酮作为一种激素类的药物，其对于男性患者和女性患者的治疗效果可能存在较大的差异。另外，长期使用左炔诺孕酮可能会产生抗癌和促癌的效应。一方面，左炔诺孕酮可在一定程度上降低子宫内膜癌等癌症的发生风险；另一方面，长期使用左炔诺孕酮可能会增加患乳腺癌的风险^[26]。因此，在使用左炔诺孕酮时，需要注意其副作用和风险，并根据个体情况选择合适的用药方案。

综上，本研究利用基于CMap数据库的药物重定位方法鉴定了10种治疗胃癌的候选药物，其中小分子化合物左炔诺孕酮可与PPIN中的枢纽节点TIMP1结合，且结合能低于已报道的药物西地尼布与其靶点PDGFRB的结合能，提示左炔诺孕酮可作为新型小分子药物来治疗胃癌，但具体的作用机制还需通过实验进一步验证。