摘要

目的:从网络药理学和分子对接的角度探讨白花蛇舌草干预骨肉瘤的分子作用机制。方法:通过中药系统药理学数据库与分析平台(TCMSP)预测白花蛇舌草的活性成分及相关靶基因;通过GEO-NCBI等数据库获得骨肉瘤相关的基因。将骨肉瘤相关基因和白花蛇舌草药物成分靶基因进行交集,得到白花蛇舌草干预骨肉瘤的共同靶基因。采用Cytoscape 3.8.1软件构建白花蛇舌草有效成分-骨肉瘤靶基因调控网络;将共同靶基因导入STRING数据库构建蛋白质-蛋白质相互作用(PPI)网络;运用DAVID在线数据库及R语言对共同靶点进行GO富集分析和KEGG富集分析;用Cytoscape 3.8.1软件筛选出白花蛇舌草干预骨肉瘤的核心靶点,然后将核心靶点与对应的白花蛇舌草活性成分进行分子对接以验证活性成分与核心靶点的结合能力。结果:共筛选出白花蛇舌草5个有效成分(MOL001659多孔甾醇、MOL001670 2-甲氧基-3-甲基-9,10蒽醌、MOL000449豆甾醇、MOL000358β-谷甾醇、MOL000098槲皮素)及其对应的靶基因共169个,3 178个骨肉瘤相关基因;靶基因和骨肉瘤相关基因交集得到白花蛇舌草干预骨肉瘤的预测靶点98个;经过PPI拓扑分析,提示相互作用紧密的8个核心基因分别是雌激素受体1(ESR1)、JUN、RAC-α丝氨酸/苏氨酸-蛋白激酶(AKT1)、G1/S特异性细胞周期蛋白-D1(CCND1)、FOS、丝裂原活化蛋白激酶1(MAPK1)和细胞肿瘤抗原p53(TP53)与MYC;GO富集分析得到P<0.05的结果共计2 205条,其中生物进程(BP)结果 2 037条,细胞组分(CC)结果 39条,分子功能(MF)结果 129条;KEGG富集分析得到P<0.05的结果共有155条,包括IL-17、细胞凋亡、肿瘤坏死因子、p53、细胞衰老等通路。分子对接结果显示AKT1、CCND1、FOS、MAPK1、TP53、MYC与槲皮素,ESR1与2-甲氧基-3-甲基-9,10蒽醌,JUN与β-谷甾醇之间亲和力均小于-7 kcal/mol,亲和力较高,显示出较好的结合能力。结论:白花蛇舌草可通过多个不同成分作用于骨肉瘤多个靶点,进而调控多条通路干预骨肉瘤发生发展过程,这可为白花蛇舌草临床应用和开发提供参考。

关键词