CJI (Traditional Chinese Medicine)
基于网络药理学和细胞验证探讨甘草抗肺癌机制
张尚龙1,闫潇1,杨秀娟1,杨志军1,连小龙1,张楠1,叶礼巧1,杨必乾 1,
殷亭湄1,付晓艳1,马趣环1,邓毅1,2
1.甘肃中医药大学,甘肃 兰州 730000;2.甘肃省中药药理与毒理学重点实验室,甘肃 兰州 730000摘要:目的 基于网络药理学探究甘草对肺癌的作用机制,并进行细胞实验验证。方法 通过TCMSP筛选甘草药效活性成分,整理并校正其相应作用靶点,通过GeneCards、OMIM、DrugBank数据库收集肺癌相关基因靶点,筛选甘草治疗肺癌的关键靶点,利用STRING、DAVID数据库对关键靶点进行蛋白相互作用(PPI)网络分析及GO和KEGG通路富集分析,采用Cytoscape构建靶点-化合物网络。检测不同浓度甘草提取物对A549细胞存活率的影响,采用Western blot 检测A549细胞Bax、Bcl-2、HIF-1α 及GLUT1蛋白表达。结果 筛选甘草活性成分88个,作用靶点213个,肺癌靶点23 458个,匹配映射得到关键靶点209个。通过疾病-药物活性成分-作用靶点网络和PPI网络分析筛选出与甘草关联程度较大的靶点43个。GO富集分析得到生物过程1 202条、细胞组分46条、分子功能64条。KEGG通路富集分析得到甘草治疗肺癌的通路226条,主要包括癌症相关通路、HIF-1信号通路、动脉粥样硬化、小细胞肺癌、幽门螺杆菌感染的上皮细胞信号转导等。体外实验表明,甘草提取物对A549肺癌细胞增殖具有显著抑制作用,可上调Bax蛋白表达,下调Bcl-2、HIF-1α及GLUT1蛋白表达。结论 本研究通过网络药理学分析得到甘草治疗肺癌的活性成分、作用靶点及信号通路,可为甘草的抗癌活性及作用机制研究奠定基础。
关键词:甘草;网络药理学;肺癌;A549细胞
中图分类号:R273.42;R285 文献标识码:A 文章编号:1005-5304(2023)10-0036-05 DOI:10.19879/j.cnki.1005-5304.202301396 开放科学(资源服务)标识码(OSID): Mechanism of Glycyrrhizae Radix et Rhizoma against Lung Cancer
Based on Network Pharmacology and Cell Verification
ZHANG Shanglong1, YAN Xiao1, YANG Xiujuan1, YANG Zhijun1, LIAN Xiaolong1, ZHANG Nan1, YE Liqiao1, YANG Biqian1, YIN Tingmei1, FU Xiaoyan1, MA Quhuan1, DENG Yi1,2
1. Gansu University of Chinese Medicine, Lanzhou 730000, China; 2. Key Laboratory of Pharmacology and Toxicology of Traditional Chinese Medicine in Gansu Province, Lanzhou 730000, China
Abstract: Objective To explore the mechanism of Glycyrrhizae Radix et Rhizoma against lung cancer based on network pharmacology, and verify by cell experiment. Methods The active components of Glycyrrhizae Radix et Rhizoma were screened through TCMSP, and the corresponding action targets were corrected. The gene targets related to lung cancer were collected by GeneCards, OMIM and Drugbank databases, and the key targets of Glycyrrhizae Radix et Rhizoma in the treatment of lung cancer were screened. The key targets were transmitted to STRING database and DAVID database for PPI network analysis, GO enrichment analysis and KEGG pathway analysis, and constructed by Cytoscape. The effects of different concentrations of Glycyrrhizae Radix et Rhizoma extract on the survival rate of A549 cells were detected. The expressions of Bax, Bcl-2, HIF-1α and GLUT1 in A549 cells were detected by Western blot. Results Totally 88 active components of Glycyrrhizae Radix et Rhizoma, 213 action targets and 23 458 lung cancer targets were screened, and 209 key targets were obtained by matching mapping.基金项目:国家自然科学基金地区基金(81960723);甘肃省科技计划项目(21JR11RA145);兰州市科技计划项目(2022-3-21);甘肃省中医药管理局课题(GZPK-2022-37);甘肃中医药大学研究生创新基金(2021CX46)通讯作者:邓毅,E-mail:dengyi@gszy.edu.cn
Through the analysis of disease-drug active components-action target network and PPI network, 43 targets with high correlation with Glycyrrhizae Radix et Rhizoma were selected. GO enrichment analysis revealed 1 202 biological processes, 46 cell components, and 64 molecular functions. KEGG pathway analysis showed that there were 226 main pathways of Glycyrrhizae Radix et Rhizoma in the treatment of lung cancer, including pathways in cancer, HIF-1 signaling pathway, atherosclerosis, small cell lung cancer, Helicobacter pylori infection and so on. In vitro experiments showed that Glycyrrhizae Radix et Rhizoma extract could significantly inhibit the proliferation of A549 lung cancer cells, up-regulate the expression of Bax protein and down-regulate the expression of Bcl-2, HIF-1α and GLUT1 protein. Conclusion Based on network pharmacology, this study preliminarily discusses the active components, action targets and signaling pathways of Glycyrrhizae Radix et Rhizoma in the treatment of lung cancer, and lays the foundation for the study of anticancer activity and mechanism of Glycyrrhizae Radix et Rhizoma.
Keywords: Glycyrrhizae Radix et Rhizoma; network pharmacology; lung cancer; A549 cells
肺癌是全球发病率最高的恶性肿瘤之一,治疗方法主要有手术治疗、放射治疗、化学治疗、分子靶向治疗等。由于肺癌早期症状不明显,恶性程度较高,很多患者在确诊时已发生肿瘤转移和扩散,无法接受手术治疗,故化疗的作用不可替代。然而大剂量或重复低剂量化疗药物会引起较强的不良反应,且导致患者产生耐药性。甘草具有补脾益气、清热解毒、祛痰1.1.4 蛋白相互作用网络构建
将共有靶点导入STRING (https://string-db.org/)数据库,在“organism”中选择“Homo Sapiens”进行搜索,获得蛋白相互作用(PPI)网络,导出TSV格式文件,将其导入Cytoscape3.8.0(http://www.cytoscape. org/)软件筛选关键靶点,绘制PPI网络。
1.1.5 GO和KEGG通路富集分析
通过Metascape(http://metascape.org/)在线软件对获得的关键靶点进行GO和KEGG通路富集分析。GO富集分析包括生物过程(BP)、细胞组分(CC)、分子功能(MF)。将分析结果前10条进行可视化。1.1.6 共有靶点-化合物网络构建
将共有靶点及化合物通过Cytoscape3.8.0软件进行可视化,构建共有靶点-化合物网络图。
1.2 实验验证
1.2.1 药物与细胞
甘草采自甘肃省金塔县,经甘肃中医药大学中药鉴定教研室鉴定为豆科植物甘草Glycyrrhiza uralensis Fisch.的干燥根和根茎。甘草提取物制备方法:称取甘草粉末100 g,加50%乙醇(1∶10,W/V)浸提30 min,超声提取(功率100 W,频率60 kHz)30 min,真空抽滤后,减压回收溶剂至50 mL,水浴加热至干浸膏,快速研磨成细粉[8],4℃冰箱保存,实验时用DMEM高糖培养基配制成50 mg/mL母液。
人肺腺癌A549细胞,甘肃省肿瘤医院惠赠。
1.2.2 试剂与仪器
DMEM 高糖培养基( Cytiva 公司,批号AH29027412),胎牛血清(浙江天杭生物科技股份有限公司,批号21030707),噻唑蓝(MTT,北京索莱宝科技有限公司,批号917Q051),胰蛋白酶、青链霉素(武汉赛维尔生物科技有限公司,批号GP2211017、GA22020011094),二甲基亚砜(DMSO,北京索莱宝
β-actin、羊抗科技有限公司,批号1121E0331),兔抗兔二抗(Abbkine公司,批号分别为ATVMA0302、
HIF-1α及ATVE16011),兔抗Bax、兔抗Bcl-2、兔抗
兔抗 GLUT1 (Abmart 公司,批号分别为 T40051、T40056、TA1009、T55360),BCA蛋白浓度测定试剂盒(北京索莱宝科技有限公司,批号20211213)。
CO2恒温培养箱(型号HF151,力康发展有限公司),酶标仪(型号RT-6100,深圳雷杜生命科学股份有限公司),电泳仪(型号DYY-6D,北京六一生物科技有限公司),超净工作台(型号B2HFsafe,上海力申科学有限公司),低温冷冻离心机(型号H1650R,湖南湘仪实验室仪器开发有限公司),荧光倒置显微镜(型号MoticamPro2058,上海捷辰仪器有限公司)。1.2.3 细胞培养
将冻存的A549细胞迅速解冻,加入完全培养基, 1 000 r/min离心5 min,弃上清液,转移至培养瓶,加入5 mL配制好的完全培养液(10%胎牛血清、1%青链霉素混合液、89%高糖DEME培养基),轻轻吹打使细胞悬浮,转移悬浮液至细胞培养瓶,置37 ℃、5%CO2培养箱进行培养,待培养瓶中细胞融合至80%~90%时,用胰酶消化并转移至新的培养瓶中。
1.2.4 MTT法检测细胞存活率
将细胞消化悬浮后以5×104个/mL接种于96孔板,待细胞贴壁生长24 h后,分为调零组(不加细胞和药物)、对照组(不加药物)及给药组,每组3个复孔,给药组用不同浓度(1、2、3、4、5 mg/mL)甘草提取
μL
物处理24 h,加入20 MTT溶液(5 mg/mL),4 h
μL
后加入150 DMSO,混匀,于酶标仪波长570 nm处检测吸光度(OD值),计算细胞存活率。存活率(%)=
给药组-OD 对照组-OD
(OD 调零组)÷(OD 调零组)×100%。1.2.5 Western blot 检测 Bax、Bcl-2、HIF-1α、GLUT1蛋白表达
将A549细胞分为对照组和甘草低、中、高剂量组,以1×105个/mL接种于6孔板中,待其贴壁后,甘草低、中、高剂量组加入不同浓度(3、4、5 mg/mL)甘草提取物处理24 h,对照组加入不含药物的培养基,收集细胞,于冰上加入细胞裂解液,分别提取总蛋白,加入上样缓冲液,100 ℃金属浴变性。根据分子量不同进行SDS-PAGE,转PVDF膜,5%脱脂牛奶封闭2 h,加入Bax、Bcl-2、HIF-1α、GLUT1一抗(均为
1∶2 000), 4℃孵育过夜,1×TBST清洗3次,每次10 min,加二抗(1∶5 000),室温孵育2h,1×TBST清洗3次,每次10 min,曝光显影,用Image J软件计算灰度值。1.2.6 统计学方法
采用SPSS22.0统计软件进行分析。计量资料用x±s表示,组间比较采用独立样本t检验,多组间比较采用方差分析。P<0.05表示差异有统计学意义。
2 结果
2.1 甘草活性成分及靶点
通过TCMSP得到甘草化合物及靶点,删去无靶点的化合物,最终得到甘草活性成分88个,靶点213个。2.2 肺癌靶点
将GeneCards、OMIM、DrugBank数据库收集到的肺癌靶点进行合并,删除重复靶点,最终得到肺癌靶点23 458个。
2.3 甘草抗肺癌靶点
将检索到的213个甘草活性成分的靶点与23 458个肺癌靶点比较分析后取重复值,得到共同作用靶点209个,作为甘草治疗肺癌的潜在靶点。
2.4 核心靶点蛋白相互作用网络
将209个潜在作用靶点导入STRING数据库,获得的PPI网络导入 Cytoscape3.8.0软件进行分析,根据Degree值与中介中心性(betweenness centrality,BC)、接近中心性(closeness centrality)高于平均值筛选出排名靠前的核心靶点43个。Degree值排序居前列的靶点有 MAPK1、 MAPK3、 MMP2、 MMP9、 AKT1、VEGFA及HIF1A等。这些靶点主要参与细胞周期进程、细胞增殖、分化、侵袭、迁移、凋亡、细胞血管生成等过程。因此,蛋白质分子之间的相关性表明这些靶点可能是甘草对肺癌起到核心调控作用的靶点。见图1。
GO和KEGG通路富集分析结果对甘草抗肺癌209个潜在靶点进行GO富集分析,共得到1 312条结果。其中BP 1 202条,主要涉及凋亡信号通路、细胞对脂质的反应、蛋白磷酸化的正调控等;CC 46条,主要涉及转录调节复合物、过氧化物酶体、囊泡腔、细胞器外膜等;MF 64条,主要涉及细胞因子受体结合、转录因子的结合、一氧化氮合酶的调节活性、蛋白质丝氨酸/苏氨酸/酪氨酸激酶活性等。各前10条结果见图2。
对甘草抗肺癌209个潜在靶点进行KEGG通路富集分析,共得到226条通路,主要涉及癌症相关通路、动脉粥样硬化、HIF-1信号通路、非酒精性脂肪肝、小细胞肺癌、幽门螺杆菌感染的上皮细胞信号转导等,前10条通路见图3。Bcl-2家族是与细胞凋亡密切相关的一类蛋白,包括Bcl-2等凋亡抑制因子和Bax等促凋亡因子。甘草主要活性成分甘草查尔酮A[9]、甘草次酸[10]、异甘草素[11]等通过调节Bcl-2/Bax表达促进肿瘤细胞凋亡,从而发挥抗肿瘤作用。HIF-1信号通路是缺氧微环境的关键通路,HIF-1α
作为调节氧稳态的核转录因子,可调节70多个基因,如血管内皮生长因子、基质金属蛋白酶等,调控恶性肿瘤的细胞增殖、血管生成、侵袭和转移[12]。因此,推测甘草可能通过影响
HIF-1α肿瘤细胞的凋亡通路及 信号通路发挥抗肿瘤作用。2.7 实验验证结果
2.7.1 不同浓度甘草提取物对细胞存活率的影响
与对照组比较,1、2、3、4、5 mg/mL甘草提取物对A549细胞增殖均有抑制作用(P<0.01),表明甘草对A549细胞存活率的影响有剂量依赖性。见表1。图3 甘草抗肺癌潜在靶点KEGG通路富集分析2.6 共有靶点-化合物网络甘草与肺癌共有靶点-化合物网络见图4。该网络包含297个节点、1 495条边。2.7.2 甘草对A549细胞 Bax、Bcl-2、HIF-1α、GLUT1蛋白表达的影响
HIF-1α、与对照组比较,甘草各剂量组Bcl-2、GLUT1表达下调,Bax表达上调(P<0.01),表明甘草HIF-1α
对 信号通路有抑制作用。见表2、图5。
注:A.对照组;B.甘草低剂量组;C.甘草中剂量组;D.甘草高剂量组图5 各组A549细胞Bax、Bcl-2、HIF-1α和GLUT1蛋白免疫印迹3 讨论目前肺癌治疗以化疗为主,而化疗药物的不良反应限制其临床应用。甘草作为临床常用药物,可通过诱导细胞凋亡、干扰细胞周期、诱导自噬、抑制糖酵解、调节免疫、调控miRNA和多种信号通路等作用机制发挥抗肿瘤作用[2]。但甘草抗肺癌的作用机制尚未完全阐明。本研究通过网络药理学对甘草抗肺癌的机制进行研究,获得88个药物活性成分、209个潜在作用靶点,通过PPI网络进一步分析获得MAPK1、MAPK3、MMP2、MMP9及HIF1A等43个核心靶点。此外,通过对209个关键靶蛋白进行GO和KEGG通路富集分析,推测甘草活性成分靶点通过癌症相关通路、HIF-1信号通路调控肿瘤细胞增殖、侵袭、迁移及凋亡。
本研究筛选出甘草提取物3、4、5 mg/mL 3种给药浓度作用于人肺癌A549细胞,并选择富集程度较高的HIF-1α
信号通路进行验证。与对照组比较,不同浓度甘草提取物均可抑制A549细胞增殖且具有浓度依赖性,甘草各剂量组能上调促凋亡蛋白Bax表达,下调抗
Bcl-2、HIF-1α及
凋亡蛋白 GLUT1表达。
综上所述,本研究通过网络药理学及实验验证方法,对甘草抗肺癌机制进行研究,验证了甘草抗肺癌多成分、多靶点、多途径的作用特点,提示甘草抗肺癌的作用机制可能与抑制肺癌细胞增殖、促进细胞凋
HIF-1α
亡及调控 信号通路有关,相关机制仍需体内实验进一步验证。本研究可为甘草的临床应用及抗肿瘤新药的开发提供依据。今后课题组将利用代谢组学及蛋白质组学等多种手段结合网络药理学对甘草抗肺癌的作用机制进行深入研究。
参考文献:
[1] 韩维维,方东军,李陆军,等.甘草化学成分及生物活性研究进展[J].
化学工程师,2022,36(2):56-58,67.
[2] 连小龙,令颖,张尚龙,等.甘草有效成分抗肿瘤作用机制及联合抗肿
瘤研究进展[J].中国实验方剂学杂志,2022,28(11):238-245.
[3] 范敏慧,何晓其,吴亚萍.新甘草酚对HepG2人肝癌细胞增殖的影响及
机制研究[J].浙江中西医结合杂志,2020,30(12):972-976.
[4] He J, Chen K, Hu Z M, et al. UGT73F17, a new glycosyltransferase from Glycyrrhiza uralensis, catalyzes the regiospecific glycosylation of pentacyclic triterpenoids[J]. Chem Commun (Camb),2018,54(62):8594-8597.
[5] 孙宇,卿德刚,张娟,等.胀果甘草中查尔酮类物质对结肠癌的作用研
究[J].中国民族医药杂志,2020,26(1):29-31.
[6] 彭江华,马建永.甘草次酸对人乳腺癌MDA-MB-231细胞的影响研究[J].
中国临床药理学杂志,2019,35(20):2616-2619.
[7] 方诗琦,冷康,段金廒,等.甘草药渣中黄酮类成分及其抗氧化活性的
研究[J].中成药,2015,37(11):2443-2448.
[8] 李越峰,辛二旦,司昕蕾,等.甘草提取工艺研究[J].中兽医医药杂志,
2019,38(6):44-47.
[9] 袁永贵,张夏炎,朱晓俊,等.甘草查尔酮A对人乳腺癌MDA-MB-231细
胞凋亡的影响[J].中国实验方剂学杂志,2021,27(20):95-100.
[10] 刘景华,赵阳,李宾,等.甘草次酸对胃癌细胞SGC7901增殖、凋亡及PI3K-Akt-mTOR通路的影响[J]. 药物评价研究, 2020, 43(12): 24442449.
[11] 覃月穆,姜泽群,马艳霞,等.异甘草素对人肺癌细胞株H460增殖、凋亡及NF-κB信号通路的影响[J].天然产物研究与开发,2018,30(5): 847-855.
[12] Lee S H, Golinska M, Griffiths J R. HIF-1-independent mechanisms regulating metabolic adaptation in hypoxic cancer cells[J]. Cells,2021,10(9):2371.
(收稿日期:2023-01-31) (修回日期:2023-02-17;编辑:陈静)