类风湿关节炎（rheumatoid arthritis,RA）是一种难治性免疫性疾病，临床主要表现为侵蚀性、对称性多关节炎损害[1,2]。流行病学资料显示，RA的发病率女性和男性分别在70～74岁和75～79岁年龄组最高，男女比例为1∶4，我国患病率为0.32%～0.36%[3]。RA发生发展的生物过程复杂，确切的发病机制尚不清楚，认为主要与遗传因素、生活环境、细菌及病毒感染有关[4]。西医针对RA的治疗主要包括非甾体类抗炎药、生物制剂及糖皮质激素等。中医近年来在RA的治疗中发挥了重要作用，中医认为RA属“痹症”范畴，其主要病机是“外邪侵袭肢体，经络闭阻，不通则痛”，治疗以“驱邪通络”为主[5]，多使用具有祛风除湿功效的药物，这些药物对缓解患者疼痛、改善患者生活质量、延缓病情进展等有重要作用，但具体作用机制尚不明确，而网络药理学为分析药物与疾病关联性的有效工具，故本研究应用网络药理学方法探索药物与疾病之间的相互作用机制，有助于更好地了解药物发挥治疗作用的机理。

海桐皮Aralia chinensis L.又名鹊不踏、鸟不宿等，为豆科植物刺桐的干皮，为常用的中草药，性温，味微咸，入肝、心、肾3经，具有镇痛消炎、祛风行气、祛湿活血之效，是治疗风湿痹痛的关键药物之一，为临床常用药物[6]。

网络药理学[7,8]是一门新兴学科，其原理是基于系统生物学和生物网络平衡阐释疾病发生发展过程、从改善或恢复生物网络平衡的整体观角度认识药物与机体的相互作用并指导发现新药。本研究基于大数据平台及计算机技术，结合网络药理学原理，探讨海桐皮治疗RA的可能作用靶点及通路，为中药挖掘提供思路。

1、研究方法

1.1海桐皮药物有效成分及相关靶点收集

以“海桐皮”为关键词在中药系统药理学数据库与分析平台（traditional Chinese medicine systems pharmacology database and analysis platform,TCMSP)(http：//lsp.nwsuaf.edu.cn/tcmsp.php）中进行检索，按照口服生物利用度（oral bioavailability,OB）≥30%和药物类药性（drug-likeness,DL）≥0.18进行筛选，最终得到海桐皮的有效成分。导出海桐皮药物相关靶点，最终得到海桐皮相关靶点数据。

1.2 RA疾病数据收集

以“rheumatoid arthritis”为关键词分别在人类基因数据库（the human gene database,GeneCards）、Disgenet数据库进行查找，物种来源选择“Homo sapiens”，下载保存合并数据后去除重复基因，最终获得RA的相关疾病基因。

1.3海桐皮相关靶点基因注释

运用Perl脚本，以uniprot蛋白数据库为参照进行海桐皮药物相关靶点基因注释，得到海桐皮相关药物靶点基因。

1.4海桐皮-RA交集基因

运用R软件中的Venn Diagrams程序包取RA疾病基因和海桐皮药物相关靶点基因的交集基因，绘制韦恩图。

1.5药物-疾病作用网络构建

将药物-疾病交集基因运用perl脚本进行映射，获得药物-疾病之间相互作用的有效成分及相关基因靶点，随后将所得有效成分及相关基因靶点数据导入Cytoscape3.7.1软件进行可视化网络构建，展现药物-有效成分-基因靶点-疾病网络。

1.6海桐皮-RA蛋白互作网络构建

利用在线String网站（https：//string-db.org/）对药物-疾病交集基因进行药物-疾病之间蛋白相互作用网络构建，筛选标准：检测物种选择人（Homo sapiens），置信值（medium confidence）取0.04，去除网络中离散或无关联的蛋白，最终得出海桐皮-RA相互作用蛋白网络图。

1.7基因本体论功能富集分析（gene ontology,GO)和KEGG富集分析（kyoto encyclopedia of genes and genomes,KEGG）

利用R软件中的RSQLite、BiocManager、org.Hs.eg.db程序包处理所得交集基因，获得交集基因的ID，将获得的基因ID再运用R软件中的colorspace、stringi、DOSE、clusterProfiler、pathview程序包进行GO生物过程富集分析及KEGG通路富集分析。

2、结果

2.1海桐皮有效成分

最终获得海桐皮有效成分6个，见表1。

表1海桐皮有效成分

2.2海桐皮相关靶点基因注释结果

将海桐皮225个相关靶点以Uniprot数据库作为筛选库，对海桐皮相关靶点基因进行注释，最终得到46个药物相关靶点基因，主要包括PTGS1、NOS2、MAPK14、GSK3B、SLC6A3等。

2.3 RA疾病基因

在GeneCards数据库中找到RA基因253个，在Disgenet数据库找到RA基因430个，去除重复基因，最后得到RA基因450个，这些基因主要有IL-6、HLA-DRB1、TNF、PTPN22、IL-10、STAT4、MIF、IL-2RA、IRF5、IL-1B、HLA-B等。

2.4海桐皮药物靶点基因与RA疾病基因交集结果

对“2.2”项所得46个有关海桐皮药物靶点基因与“2.3”项所得450个RA疾病基因运用R×643.5.3软件进行计算分析，结果得到14个交集基因，分别为ADRA2B、ADRA1B、ADRB2、ADRA1D、NOS2、NR3C2、ADRA2A、SLC6A2、ADRB1、ADRA1A、ADRA2C、SLC6A4、DRD2、PPARG，绘制Venny图，结果见图1。

图1药物-疾病靶点基因Venny图

2.5药物-疾病相互作用网络

将获得的14个交集基因运用Perl脚本处理后导入Cytoscape3.7.1软件进行可视化网络构建，展现药物、疾病、基因、分子相互作用网络，见图2。

图2海桐皮主要活性成分-靶点-药物-疾病网络

2.6海桐皮-RA蛋白互作网络

最终得出的相应PPI网络有13个节点，32条边，平均每个节点的度为4.57，平均局部聚类系数为0.61,PPI enrichment P<1.0e-16，见图3。

图3交集基因PPI网络图

2.7海桐皮-RA蛋白互作网络核心蛋白

在海桐皮-RA蛋白相互作用网络中，下载其string＿interactions文件，并以tsv格式保存，运用R×643.5.3软件绘制海桐皮在该网络中核心蛋白条形图，结果显示ALB、IL-6、NOS3、VEGFA、FOS等在该网络中起关键作用，见图4。

图4 PPI网络核心蛋白条形图

2.8 GO与KEGG

将所得14个交集基因导入R×643.5.3软件进行GO生物过程富集分析和KEGG通路分析，结果显示海桐皮治疗RA的主要生物过程包括了肾上腺素受体活性、G蛋白偶联的胺受体活性、儿茶酚胺结合、肾上腺素受体结合等过程，KEGG信号通路富集分析显示主要涉及的信号通路包括了cGMP-PKG信号通路、神经活性配体-受体相互作用信号通路及cAMP等信号通路等。GO分析结果见图5,KEGG通路富集分析见图6—8及表2。

图5 GO富集分析条形图

图6 KEGG富集分析条形图

图7 cGMP-PKG信号通路图

图8神经活性配体-受体相互作用信号通路图

表2海桐皮治疗RA的KEGG通路富集分析结果

3、讨论

RA是一种慢性自身免疫性疾病，其特征性表现主要是进行性、对称性、多滑膜关节炎和关节外表现，最终导致关节畸形、功能障碍，严重影响患者生活质量[9]。对于RA的发病机制，目前研究认为其主要与免疫紊乱、细胞因子网络、遗传因素及环境因素等有关[10,11,12,13,14]。在分子机制方面，认为多种细胞因子如IL-B12、微小RNA(microRNAs,miRNA）、基因多态性及信号传导通路在RA的滑膜和骨质破坏过程中发挥作用[15,16,17,18,19,20,21]。近年来，中医药在治疗RA中发挥重要作用，中药对降低促炎性细胞因子产生及减轻患者疼痛与炎症反应有重要作用，对RA具有治疗作用[22,23,24]。海桐皮是临床常用于治疗风湿痹痛的有效药物，海桐皮中含有刺桐灵碱、氨基酸和有机酸等，由海桐皮组成的海桐皮汤化裁常用于治疗骨关节炎、骨折等，对骨关节修复具有重要作用[25,26,27]。

本研究通过TCMSP筛选得到海桐皮有效成分6个，分别为赤藓碱（Erythraline）、异巴巴胆（isobavachin）、豆甾醇（Stigmasterol）、赤藓碱（Erysotrine）、菜豆蛋白（Phaseolin）、菜豆素（Phaseollidin）。其中赤藓碱（Erythrina）属的第一种化学物质在1930年被报道，其中分离出类胡萝卜素，一种螺环状生物碱，显示出咖喱样活性[28]，赤藓碱生物碱通过调节烟碱样受体对中枢神经系统产生影响，其在中枢神经系统中具有重要生物活性，其对疼痛、焦虑、失眠有良好调控作用[29,30]，虽然目前未见其治疗RA的报道，但其治疗疼痛的特性可能缓解RA疼痛。豆甾醇又名豆固醇、甾醇属植物甾醇，具有较强的生理活性和表面活性，可抗炎、抗氧化、抗癌、降胆固醇等[31],CHEN等[32]研究发现，豆甾醇能抑制基质金属蛋白酶表达，抑制软骨降解。Erysotrine属赤藓碱类，赤藓碱的mulungu提取物能降低IL-4和IL-5水平，提高IL-13和干扰素γ(interferonγ,IFN-γ）水平，表明其对炎症治疗有一定潜在价值。异巴巴胆、菜豆蛋白、菜豆素目前尚未见相关报道，可为进一步研究海桐皮治疗RA提供新方向。

此外，本研究发现了14个海桐皮与RA的交集基因，通过STRING数据库得出海桐皮治疗RA的核心蛋白包括ALB、IL-6、NOS3、VEGFA、FOS、PPARG、NR3C1、CRP等20个核心基因，其中ALB基因编码是人类血液中最丰富的蛋白质。该蛋白在调节血浆胶体渗透压中起作用，并充当多种内源性分子（包括激素、脂肪酸和代谢产物）以及外源性药物的载体蛋白。该蛋白还表现出广泛底物特异性的酯酶样活性。IL-6为在炎症和B细胞成熟中起作用的细胞因子。另外，已证明编码蛋白质是患有自身免疫性疾病或感染者发烧的内源热原，该蛋白主要在急性和慢性炎症部位产生，并通过IL-6受体α诱导转录炎症反应。有研究表明，IL-6与RA有密切关系，参与了RA骨破环、炎症反应等[33,34]。一氧化氮合酶3(nitric oxide synthase,NOS3）是一种蛋白质编码基因，在多种过程中充当生物介质，包括神经传递以及抗微生物和抗肿瘤活性。NOS3基因参与了抑制血管内皮表面血小板和白细胞的黏附、聚集及血管平滑肌细胞增生[35]。同时，本研究发现，肾上腺素受体活性、G蛋白偶联的胺受体活性、儿茶酚胺结合等均参与了海桐皮治疗RA的过程，其中肾上腺受体作为联系神经系统与免疫系统的重要受体，可能参与了对RA的神经免疫调节，介导了RA的炎症反应过程及细胞凋亡过程[36,37,38],G蛋白偶联的胺受体（G proteincoupled receptors,GPCRs）在免疫应答调节中起重要作用，GPCRs家族包括β2肾上腺素受体、趋化因子受体，其中趋化因子在RA致炎机制中一方面通过迅速刺激依赖凝集素的T淋巴细胞黏附，趋化特定的细胞到炎症部位，导致免疫细胞穿过血管到达关节，另一方面通过其他生物功能调节因子，包括淋巴细胞活化、分化以及酶和细胞因子的释放、细胞毒作用及其他起作用[39,40,41]。儿茶酚胺被认为是神经细胞或内分泌细胞分泌的神经递质或激素，有研究表明，RA的炎症机制与免疫细胞中内源性儿茶酚胺的活性增强相关，并通过儿茶酚胺调控促炎性细胞或炎症抑制细胞增殖与功能实现[42]。

cGMP-PKG信号通路、神经活性配体-受体相互作用、钙信号通路、cAMP信号通路等同样参与了海桐皮治疗RA等生物过程，其中cGMP-PKG信号通路是指环鸟苷酸/蛋白激酶G信号通路[43]，其主要功能与平滑肌松弛，抑制非选择性离子通道的Ca2+电导率有关。神经活性配体受体相互作用信号通路是质膜上所有与细胞内外信号通路相关的受体和配体，与环境信息处理、信号分子与相互作用有关[44]，钙信号通路在细胞存活、增殖、代谢中发挥重要作用，其作用机理为Ca2+进入细胞内，激活各种Ca2+替代酶而产生作用[45,46]。研究发现，多种钙通道对于破骨细胞的分化和功能必不可少，破骨细胞分化和功能的阻断对抑制RA骨破坏有效[47]。cAMP信号通路又称环磷酸腺苷信号通路，典型的cAMP信号通路是由细胞外激动剂结合特异G蛋白偶联受体导致G蛋白三聚体分解为G蛋白α亚单位和G蛋白βγ亚单位，cAMP与其调控亚单位结合释放催化亚单位，催化亚单位进一步磷酸化位于胞质或胞核中的底物，包括酶、离子通道、染色体蛋白和转录因子等[48]。

综上所述，本研究基于网络药理学方法对海桐皮进行化学成分数据库构建、靶点预测、药物-靶点网络构建和蛋白相互作用网络分析GO生物过程和KEGG通路富集分析，初步阐述了海桐皮治疗RA的作用机制，但具体疗效尚需进一步验证。

参考文献:

[1]葛均波,徐永健.内科学[M].北京:人民卫生出版社,2015:808.

[2]李冀,李想,高彦宇.中医药治疗类风湿性关节炎研究进展[J].辽宁中医药大学学报,2019,21(12):5-8.

[3]袁林,吴金玉,唐婕,等.白虎加桂枝汤联合西药治疗类风湿关节炎风湿热痹证临床观察[J].辽宁中医药大学学报,2019,21(12):168-171.

[4]韩隆胤,王强,钱凯,等.基于网络药理学探讨南续益母片治疗类风湿关节炎的分子机制[J].中华中医药学刊,2020,38(4):149-152.

[5]吴勉华,王新月.中医内科学[M].北京:中国中医药出版社,2012:7-427.

[6]陈蔚文.中药学[M].北京:人民卫生出版社,2014:150.

[7]黄秀芳,庾国桢,童晶晶.基于网络药理学分析陈皮的药理作用机制[J].中成药,2019,41(12):3038-3045.

[8]宋红新,马旭冉,王敦方,等.基于网络药理学的黄芩汤治疗溃疡性结肠炎的潜在机制研究[J].药学学报,2020,55(2):247-255.

[9]姚红梅,何佳欢,姜锦林.类风湿关节炎的中医药治疗近况[J].中医临床研究,2019,11(31):145-148.

[10]韩海燕,沈敬华.类风湿关节炎发病机制的研究进展[J].世界最新医学信息文摘,2019,19(93):66-69.

[11]林基勇,庞宇舟,方刚,等.基于GoPubMed类风湿关节炎发病机制的文献计量学分析[J].现代预防医学,2019,46(3):563-567.

[12]叶霖,王友莲.类风湿关节炎的中西医发病机制综述[J].江西医药,2017,52(1):88-91.

[13]王绍文,余林健,万晓春,等.类风湿性关节炎免疫发病机制的研究进展[J].集成技术,2015,4(4):64-74.

[14]魏蕾,姜林娣.类风湿关节炎病因和发病机制研究进展[J].医学综述,2015,21(9):1548-1551.

[15]伍斌,鲁延富,姜凤良.类风湿关节炎发病机制的研究进展[J].医学综述,2014,20(23):4249-4251.

[16]祝静,邢艳.成纤维样滑膜细胞microRNA与类风湿关节炎[J].中国免疫学杂志,2015,31(2):277-281.

[21]田照,庞宇舟,袁德培,等.壮药龙钻通痹方影响类风湿关节炎大鼠成骨细胞转录因子Smad5的时序性分析[J].中国老年学杂志,2019,39(21):5307-5311.

[22]蔡萧君,颉彦鹏,王磊.基于三妙丸抗类风湿性关节炎作用机制的网络药理学研究[J].辽宁中医杂志,2020,47(5):148-152.

基金资助:国家自然科学基金(81973976);国家中医药管理局第六批全国老中医药专家学术经验继承项目(国中医药人教发[2017]29号);广西中医药大学广西一流学科建设开放课题(2018XK065,2018XK055);

文章来源:罗业浩,唐秀松,许栋涵等.基于网络药理学的海桐皮治疗类风湿关节炎机制研究[J].西部中医药,2023,36(11):1-7.