基于网络药理学与分子对接探讨铁皮石斛抗炎性衰老的作用

石苒羲 ，杨 涛 ，程 庆 ，赵梁辰 ，郭嘉晖 ，王利梅

（1）昆明医科大学第二临床医学院;2）第一临床医学院;3）实验动物学部，云南 昆明 650500）

衰老是生物体不可逆转的自然进程，随年龄增长各系统功能出现退行性改变，其中便有免疫功能的退化[1]，使免疫细胞分化、增殖能力下降，出现表型改变和功能缺陷，促炎因子产生增加，血清中炎症介质水平升高[2-4]。2000 年，Franceschi[5]第一次提出“炎性衰老”的概念，随着年龄的增长体内各组织发生相关低度炎性反应，被称为炎性衰老[2，4]。在炎性衰老过程中有两个重要的炎性因子IL-6 和IL-10[6-7]。IL-6 作为介导急性期炎症反应的重要细胞因子，是炎症介质网络的关键成员[7-8]，通过促炎作用对机体造成影响，引起衰老带来的如老年性痴呆、代谢紊乱、动脉粥样硬化等年龄相关的健康问题[9]，它的水平对多种疾病的诊断和严重程度判断具有预测价值[8，10]，能作为炎症感染的严重程度、预后判断及老年人炎性衰老的预测指标[6，11-12]。IL-10 是公认的免疫抑制因子，具有很强的抗炎作用，可以抑制单核巨噬细胞的增殖和活性，使炎症因子的合成与释放得到有效抑制[13]，可以阻断各种炎症通路[14]，在慢性炎症的发展过程中能起到一定的延缓作用，对减轻机体的炎症损伤具有重要意义。老年小鼠对肺炎链球菌的先天免疫应答与炎症加剧中有数据表明，炎症加剧与IL-10 介导的免疫调节减弱有关，具体表现在IL-10 浓度的降低[15]。

铁皮石斛是兰科石斛属草本植物，有抗衰老、抗肿瘤、降血糖、抗氧化、增强免疫等药理作用，能够有效防止亚健康，有延缓衰老的效果[16-18]。铁皮石斛是一种有效的免疫调节剂，对免疫功能有增强的作用，能降低TNF-α、IL-6 的产生，抑制巨噬细胞炎症[19-21]。衰老时体内细胞因子改变使机体处于慢性低度炎性状态，调节细胞因子、减少炎症反应可达到抗炎性衰老的目的。石斛具有优越的药理作用，具有无毒、无致突变、无遗传毒性的“三无”作用。石斛具有提高免疫、抗衰老的功效，但作用机制尚未被阐明。因此，本研究基于网络药理学方法，利用分子对接技术验证关键成分与靶点结合情况，预测石斛抗衰老的潜在机制，通过实验进一步验证，为石斛抗衰老的后续研究及临床应用提供参考。

1 材料与方法

1.1 石斛抗炎性衰老网络药理学研究

1.1.1 石斛主要成分检索及靶点预测利用ETCM 检索石斛主要成分，通过Pubchem 查询化合物SMILES 号，用SMILES 号在Swiss TargetPrediction数据库进行成分靶点的预测，获得石斛主要成分作用靶点。将成分靶点导入到Cytoscape 中，构建主要成分靶点网络。

1.1.2 石斛与衰老大鼠作用靶点筛选从Rat Genome Database（https://rgd.mcw.edu/）网站下载SD 衰老大鼠基因。将石斛主要成分作用靶点和衰老大鼠基因进行交集，利用Venny 网站输入二者信息，得到石斛主要成分作用于衰老大鼠的靶点。

1.1.3 蛋白-蛋白交联互作网络构建在STRING数据库，输入交集基因，构建蛋白-蛋白相互作用（protein-protein interaction，PPI）网络，将结果导入Cytoscape 软件，根据Degree 值作图。

1.1.4 靶点GO 和KEGG 富集通路分析利用DAVID 数据库进行GO 富集分析和KEGG 富集通路分析，得到富集信息。选择P≤ 0.01 且位靠前列的20 条显著通路，用富集的名字、GeneRatio、FDR、count 值通过微生信平台构建富集气泡图。

1.2 石斛活性成分-核心靶点分子对接

利用Cytoscape 软件中Centiscape2.2 对关键靶点进行筛选，主要参考 Betweenness、Closeness 和Degree 值选取靶点。利用Cytoscape软件中cytoNCA 对活性成分进行筛选，degree ＞100，并参考相关文献。从ZINC 数据库下载主要药物活性成分的2D 结构，保存为MOL2 格式，从UniProt 数据库找到相应蛋白并在AlphaFold Protein Structure Database 网站下载靶点蛋白的3D 结构，保存为PDB 格式，用pymol 对大分子进行前处理。利用Autodock 软件去水、加氢，对受体和配体进行分子对接，根据对接数值评价其活性，并利用pymol 进行可视化展示。

1.3 石斛抗炎性衰老的动物体内实验

1.3.1 实验动物雄性SD 大鼠30 只，20 月龄，体重（300±20）g，购自昆明医科大学实验动物学部，使用许可证号SYXK（滇）K2015-0002。经昆明医科大学实验动物伦理审查委员会审核批准，伦理审查批准号：KMMU20220834。

1.3.2 实验试剂及配制铁皮石斛粉末（含30%石斛多糖），陕西西安百川生物科技有限公司购买，将蒸馏水与石斛粉末按1∶1 配制为溶液，4 ℃冰箱保存。百令胶囊0.5 g/粒，生产厂家杭州中美华东制药有限公司，1∶20 配制为溶液，4 ℃冰箱保存。IL-6、IL-10 酶联试剂盒，购自上海酶联生物科技有限公司。

1.3.3 实验动物分组和给药自然衰老SD 雄性大鼠30 只，随机分为5 组。分别为空白对照组、石斛组（100 mg/kg、200 mg/kg、400 mg/kg）、百令胶囊阳性对照组（40 mg/kg）。空白对照组每天灌胃蒸馏水，每组每天定时灌胃1 次，连续灌胃10 周。

1.3.4 取材及检测方法心脏取血，血样离心后取上清液待测。ELISA 测定时每孔加入50 µL 血清和100 µL HRP 酶标记，37 ℃温育60 min，稀释洗涤液洗涤5 次并拍干。加入显色剂A 和显色剂B，37 ℃避光显色15 min。最后加入终止液50 µL，以空白组调零，450 nm 波长下依次测量各孔的吸光度。

1.4 统计学处理

网络药理数据通过 Genecard、KEGG PATHWAY Database 等数据库进一步分析。IL-6与IL-10 表达量数据用SPSS 软件分析，结果用单因素方差检验，最小显著差异法（LSD）进行两两比较，P＜0.05 为差异具有统计学意义。

2 结果

2.1 石斛抗炎性衰老网络药理学研究结果

2.1.1 石斛的主要活性成分、作用靶点、蛋白互作网络检索到石斛主要成分28 项，得到492 个作用靶点，见图1A。绿色菱形代表活性成分，紫色椭圆形代表靶点。通过Venny 软件，输入成分靶点492 个，衰老大鼠基因413 个，绘制交集韦恩图，见图1B，获得石斛对衰老大鼠的作用靶点共64 个。在STRING 数据库输入交集基因靶点绘制PPI 网络图，见图1C，节点代表作用靶点，边代表两两之间互相作用关系，节点颜色越蓝大小越大说明越重要。显示石斛抗炎性衰老的关键核心靶点有AKT1、STAT3、NOS3、GSK3B。

图1 网络药理学分析结果Fig.1 Network pharmacological analysis results

2.1.2 石斛主要成分对衰老大鼠GO 富集分析和KEGG 富集通路分析利用DAVID 数据库对交集靶点进行生物富集分析得到GO 功能富集条目见图2A～2C。结果显示药物作用于细胞质、质膜、核、线粒体等，与ATP 结合、蛋白质丝氨酸/苏氨酸/酪氨酸激酶活性、蛋白激酶活性、跨膜受体蛋白酪氨酸激酶活性、血管内皮生长因子结合、蛋白激酶C 活性等分子功能有关，参与蛋白质磷酸化、凋亡过程、有机环状化合物反应、对脂多糖反应、EPK1 和EPK2 级联反应正向调节、肽基丝氨酸磷酸化、蛋白质自磷酸化、细胞迁移正向调节、肽基苏氨酸磷酸化等生物过程，其中石斛参与生物过程中衰老最为显著。

图2 GO 分析和KEGG 分析结果Fig.2 GO analysis and KEGG analysis result

KEGG 富集到较为显著20 条信号通路见图2D，主要有：PI3K-AKT、Calcium、EGFR 酪氨酸激酶抑制剂耐药性、HIF-1、Apelin、VEGF等信号通路，通路上多个基因与炎症的发生发展关系密切。PI3K-AKT 是炎症反应中的关键通路。VEGF、HIF-1、Apelin 3 个信号通路中多个基因参与NO 的合成，间接影响炎性因子合成。在PI3K-AKT、RAP1 信号通路中，通过PI3K→AKT→NF-κB调节细胞因子合成促进炎症。在AGE-RAGE signaling pathway in diabetic complications 中，可能通过PKC 和AKT 调节NF-κB（相关通路PI3K-AKT 通路、钙通路），间接调节IL-6、TNF-α、IL-1 合成。

2.2 主要活性成分和核心靶点的分子对接

对靶点和成分进行筛选，见表1、表2，Degree（节点的度值）表示与此节点有联系的节点的数量。对靶点AKT1、STAT3、NOS3、GSK3B和成分齿叶黄皮素（Dentatin）、吴茱萸次碱（Rutaecarpine）、紫花前胡（Nodakenetin）、吴茱萸碱（Evodaiamine）分别进行分子对接，见图3、图4、表3，紫色分子为靶点，绿色部分为靶点与活性成分结合部位，蓝色分子为活性成分。以结合能作为判断成分与靶点结合能力大小的标准，结合能小于0 表示可以自由结合，越小更易发生结合。结合后出现氢键表示二者结合稳定，提示较好的结合关系。结果显示各成分与各靶点间均能形成数量不等的氢键，AKT1 与各活性成分、GSK3B与除紫花前胡外的成分、STAT3 和NOS3 与齿叶黄皮素、吴茱萸碱结合能均小于-5 KJ/mol，表示结合情况良好。

图3 分子对接可视化局部图Fig.3 Partial diagram of molecular docking

图4 分子对接可视化整体图Fig.4 Overall diagram of molecular docking

2.3 石斛对IL-6 与IL-10 蛋白表达的影响

ELISA 测得IL-6、IL-10 浓度，使用SPSS 软件分析，见图5。

图5 不同给药对IL-6、IL-10 的影响Fig.5 The effect of different susstances on the expression of IL-6、IL-10

中剂量组血清中IL-6 含量较空白组（P＜0.01）、百令组（P＜0.01）、低剂量组（P＜0.05）相比均降低，差异具有统计学意义，提示中剂量石斛可以更好的降低IL-6 含量。低剂量组、高剂量组、百令组血清与空白组相比、高剂量组与低剂量组相比，IL-6 的含量也有所下降，但差异无统计学意义（P＞ 0.05）。

低剂量组较空白组血清中IL-10 含量升高，差异具统计学意义（P＜0.05），较其它组含量升高但不具有统计学意义，提示低剂量石斛能更好的增加IL-10 含量。中剂量组、高剂量组、百令组与空白组相比、中剂量组与高剂量组与百令组相比，IL-10 含量增加但差异无统计学意义（P＞ 0.05）。

3 讨论

本研究通过网络药理学筛选石斛主要活性成分，预测其作用靶点和抗衰老作用机制，对石斛抗衰老的机制进行了初步探讨，并通过分子对接和动物实验进行验证，阐明了石斛抗炎性衰老的物质基础和作用机制，可在实践中为石斛的应用提供新的启示。

本研究检索到石斛活性成分28 项，获得潜在靶点492 个、交集靶点64 个。通过PPI 网络得到核心靶点有AKT1、STAT3、NOS3、GSK3B，这4种蛋白都对炎症通路具有调节作用。NF-κB可调节IL-6 的生成[22]，AKT1（蛋白激酶Bα）是其上游基因[23]，是PI3K 通路的中心环节，抑制AKT会使IL-6 表达下降[24]，而活化的STAT3（信号转导和转录激活因子3）可激活NF-κB信号通路，加剧炎症反应[25]；GSK3B 是糖原合成酶激酶3β，其活性受AKT 的调节[26]，上升可以促进IL-10合成增加；NOS3（一氧化碳合酶3）可诱导NO 合成来促进IL-6 的生成[27]。对石斛成分进行筛选得到齿叶黄皮素、紫花前胡、吴茱萸次碱、吴茱萸碱四种主要活性成分。紫花前胡可通过抑制NFκB信号通路而对炎症有抑制作用[28]。吴茱萸碱是生物碱类的一种，可抑制巨噬细胞、单核细胞、小胶质细胞抗炎因子的产生，并减少中性粒细胞中活性氧的产生，具有抗炎活性[29]。吴茱萸次碱可改善炎症小鼠模型的炎症反应，调节免疫平衡[30]。

在KEGG 富集到的通路中，PI3K-AKT 是炎症反应中的关键通路，激活可以促进炎症。PI3KAKT、RAP1、AEG-RAGE、EGFR 等信号通路中抑制AKT 可对NF-κB调控使IL-6 表达下降[31]。在FOCAL ADHESION、PI3K-AKT、LIPID AND ATHEROSLEROSIS 等信号通路中，提升GSK3B活性可增加IL-10 转录。同时，下调AKT 和NFκB的活性也可促进IL-10 的转录，使体内炎症状态得以抑制[26]。在 Apelin、LIPID AND ATHEROSLEROSIS、EGFR 等信号通路，抑制STAT3 活性可减少IL-6 生成[25]。VEGF、Apelin、HIF-1、Calcium 等信号通路 中NOS1、NOS2、NOS3 可促进NO 合成，增强促炎介质如IL6 和IL8 的合成，导致炎症严重程度增加[32]，对NO进行调节，可在一定程度上减缓炎症反应的进程，改变体内炎性状态。

分子对接结果发现，AKT1 与各活性成分、GSK3B 与除紫花前胡外的其它成分、STAT3 和NOS3 与齿叶黄皮素、吴茱萸碱对接结果良好，表明石斛主要通过以上四种成分发挥抗炎性作用，主要作用靶点是AKT1。吴茱萸碱和吴茱萸次碱为抑炎物质[29-30，33]与GSK3B 结合可促进IL-10产生，吴茱萸碱还可与STAT3 作用抑制IL-6 生成。紫花前胡作用于AKT1 减少IL-6 的产生，也可间接调节GSK3B 活性进而调节IL-10 表达[28]。结合动物实验，石斛组能降低血清中IL-6 的含量、增加血清中IL-10 的含量，石斛可以调节细胞因子的含量来改善机体炎性状态，进而对炎性衰老起延缓作用，与网络药理学预测结果一致。但IL-6 含量与IL-10 含量变化没有直接关系，说明不同剂量石斛对IL-6、IL-10 的调节作用不同。

综上，石斛可通过多通路、多基因对细胞因子进行调节进而影响炎症反应而具有抗炎性衰老的功能。主要通过齿叶黄皮素、吴茱萸次碱、紫花前胡、吴茱萸碱等活性成分调控AKT1、STAT3、NOS3、GSK3B 等关键基因以及 PI3K-Akt、Calcium、EGFR 酪氨酸激酶抑制剂耐药性、HIF-1、Apelin、VEGF 信号通路调节IL-6、IL-10 表达量来达到抗炎、抗衰老的作用。但动物实验样本数据较小，结果是否存在普遍性还有待探究，同时机体免疫调节系统是一个极其复杂的网络系统，不能较好地完全反映免疫功能的变化。实验针对铁皮石斛的免疫调节研究仍处于比较浅显的阶段，仅证明了石斛下调促炎因子IL-6 的含量，上调抑炎因子IL-10 含量。但更深入的机制仍有待于继续深入研究。