基于生物信息学技术探讨补骨脂致肝损伤的作用机制

余 丹，迟骁玮 （. 南京中医药大学医学院·整合医学学院, 江苏 南京 003；. 国家卫生健康委医药卫生科技发展研究中心, 北京 00044）

补骨脂（Psoralea corylifoliaLinn）为豆科植物，以干燥成熟果实入药，性温，味辛、苦，归肾、脾经。有温肾助阳，纳气平喘，温脾止泻等功效，主治骨质疏松症、骨软化症、关节痛、哮喘、白癜风、银屑病、遗尿尿频、肾虚作喘等疾病，在中医临床治疗中被广泛应用于脾肾阳虚的治疗，但长期或者超量使用补骨脂及其制剂可导致肝损伤，临床表现主要为皮肤巩膜黄染、尿黄、乏力和食欲不振，同时伴有丙氨酸转移酶（ALT）、总胆红素（T-Bil）和直接胆红素（D-Bil）等生化指标不同程度升高，肝脏肿大或肝脏弥漫性病变。实验研究也表明，长期或大剂量服用补骨脂可引起大鼠不同程度的肝损伤。然而，补骨脂肝损伤成分和潜在机制至今尚未得到阐明。在本研究中，我们通过生物信息学工具评价补骨脂诱导肝损伤的潜在作用机制，为补骨脂肝损伤的临床和基础研究提供参考依据。

1 材料与方法

1.1 软件与数据库

TCMIP (http://lsp.nwu.edu.cn/, Version 2.3)；Pharm Mapper (http://59.78.96.61/pharmmapper/)；Uniprot 数据库中UniprotKB 检索功能（http://www.uniprot.org/）；String 数 据 库（https://string-db.org/,Version 4.5.0）；Disgenet（https://www.disgenet.org/）；DAVIA 基因富集分析网站 (http://david.ncifcrf.gov/summary.jsp)；Cytoscape 软 件(Version 2.7.1)；Chembiodraw Ultra 软件 (Version14.0)

1.2 补骨脂肝损伤成分筛选与靶点预测

通过TCMIP 数据库获取补骨脂所有化学成分，再结合文献挖掘与整理，筛选出补骨脂的主要成分，并对成分进行预测。

登陆Pharm Mapper 服务器，上传上述补骨脂成分，获得虚拟筛选结果。利用UniProt 数据库中UniProtKB 搜索功能，输入蛋白名称，限定物种为人，预测补骨脂成分的作用靶点。

在Disgenet 数据库中搜索药物性肝损伤相关的基因，去除重复基因和假阳性基因，与上述Pharm Mapper 服务器中返回的靶点相匹配，得到补骨脂成分诱导肝损伤的潜在作用靶点。

1.3 网络构建与分析

将补骨脂诱导肝损伤的作用靶点导入String数据库，利用Multiple proteins 工具，选择“Homo Sapiens”为生物体，获取蛋白相互作用。将补骨脂成分和作用靶点信息以源节点和节点直接导入Cytoscape3.7.1，构建药材-成分-靶标-疾病网络。借助Network Analyzer 工具对网络进行拓扑参数分析，并筛选出补骨脂诱导肝损伤的关键靶蛋白。

1.4 基因本体功能（GO）注释和京都基因与基因组百科全书（KEGG）富集分析

用DAVID 6.7 (http://david.ncifcrf.gov/summary.jsp) 对补骨脂药物性肝损伤PPI 网络中的节点蛋白进行GO 注释和KEGG 通路分析，并采用R 软件对结果进行图形化表示。

2 结果

2.1 补骨脂成分的相关信息

从TCMIP 数据库中共筛选出的32 个相关化合物，其中补骨脂素、补骨脂酚、补骨脂查尔酮、异补骨脂素和甲氧补骨脂素等22 个符合Lipinski5规则，而补骨脂素、异补骨脂素和补骨脂酚为补骨脂主要成分。因此，上述22 个化合物全部纳入进一步研究，化合物名称及详细信息见表1。

表1 补骨脂成分及相关信息

2.2 补骨脂成分肝损伤的靶点预测

将补骨脂22 个成分输入Pharm Mapper 中，得到所有潜在蛋白质靶点，经UniProt 转化为基因ID，并与Disgenet 数据库中8 类药物性肝损伤挖掘的疾病基因ID 根据匹配度由高到低进行排序，去除重复后，匹配共得到31 个补骨脂潜在的肝损伤作用靶点，见表2。

表2 补骨脂成分潜在肝损伤作用靶点

2.3 补骨脂成分-靶点-信号通路-肝损伤网络构建与分析

采用Cytoscape 软件的Merge 功能构建补骨脂成分-靶点-信号通路-肝损伤网络模型，如图1 所示（绿色椭圆代表靶点，黄色椭圆代表成分，粉红色菱形代表疾病，紫色椭圆表示补骨脂），节点代表活性或靶点，边代表成分、作用靶点和疾病间的相互关联，其中共有60 个节点，327 个边。节点的大小与度的大小呈正比关系，节点越大表示该节点的度越大。依据拓扑分析，13 个成分的度值大于10，分别为补骨脂素、补骨脂定、补骨脂酚、异新补骨脂查尔酮、异补骨脂黄酮、补骨脂异黄酮醛、异补骨脂定、双羟异补骨脂定、新补骨脂查尔酮、补骨脂二氢黄酮甲醚、补骨脂乙素、补骨脂查尔酮、新补骨脂异黄酮；具有较高度值的靶点蛋白有血清白蛋白（ALB）、谷胱甘肽S-转移酶P（GSTP1）、运甲状腺素蛋白（TTR）和过氧化物酶体增殖物激活受体γ（PPARG）。

图1 补骨脂成分-靶点-肝损伤网络

采用DAVID 在线富集平台进行生物学功能分析，获得103 个GO 注释（黄色表示生物学过程、绿色表示细胞组成、蓝色表示分子功能）。GO 分析表明，靶蛋白大多富集于RNA 聚合酶Ⅱ启动子启动转录、类固醇激素介导的信号通路、信号转录、DNA 模板转录正调控、RNA 聚合酶Ⅱ启动子转录的正调控、氧化还原反应、固有免疫应答、细胞对脂多糖（LPS）的反应、细胞氧化解毒等生物学过程，见图2。此外，KEGG 富集分析表明，5 条通路（P＜0.05）受补骨脂的影响，包括化学物致癌、AMPK信号通路、PPAR 信号通路、肝药酶P450（CYP）有害物代谢和谷胱甘肽代谢通路等，结果见图3。

图2 补骨脂肝损伤靶点的GO 富集图

图3 补骨脂肝损伤靶点的KEGG 富集图

3 讨论

在临床上，补骨脂被广泛应用于治疗骨质疏松、银屑病等多种疾病，其所导致的药物性肝损伤已引起人们关注，但对于其肝损伤成分及分子机制的研究尚不深入。为了探究补骨脂所致肝损伤的作用机制，本研究从补骨脂的成分出发，借助生物信息学技术构建了补骨脂药材-成分-靶点-信号通路-肝损伤网络，分析补骨脂中成分与靶点的作用关系，为其导致肝损伤作用机制研究提供了参考。

通过数据库检索，我们发现补骨脂主要含有补骨脂素、异补骨脂素、补骨脂酚、补骨脂定、新补骨脂异黄酮和补骨脂二氢黄酮甲醚等22 个成分[1-5]，其中补骨脂素、异补骨脂素和补骨脂酚等为主要成分[5-6]。补骨脂中个别成分已被发现具有导致肝损伤的作用，如研究发现补骨脂素与补骨脂酚可通过降低BSEP 蛋白的表达，升高NTCP 和CYP7A1 蛋白的表达，使得进入肝细胞的胆汁酸增多，进而肝细胞中胆汁酸的含量异常升高，最终导致胆汁淤积性肝损伤[7-9]；此外，补骨脂素通过上调细胞周期蛋白E1 和p27 蛋白水平并下调细胞周期蛋白D1 的水平，使肝细胞周期阻滞进而导致肝再生和代偿能力下降而诱发小鼠轻度肝损伤和L02 肝细胞损伤；而异补骨脂素通过抑制MRP2 和MRP3 导致HepG2的细胞肝损伤[10-11]。补骨脂酚在体外的生物利用度远低于补骨脂素和异补骨脂素，但在体内补骨脂酚的CYP450 代谢作用较强[12]，补骨脂酚在体内复杂环境中的吸收和代谢可能是导致其在体内外肝损伤差异的主要原因。

通过GO 分析整合预测，得到补骨脂导致肝损伤生物功能包括氧化还原反应、固有免疫应答、细胞对LPS 的反应、细胞氧化解毒等。目前已知，氧化还原反应与免疫应答在药物性肝损伤发病机制中起关键作用。有研究发现，可溶性T 细胞免疫球蛋白结构域和含黏蛋白结构域分子-3（sTIM-3）通过促进自噬和调节单核/巨噬细胞功能，减轻炎症反应和肝损伤[13]。氧化应激介导的脂质过氧化可导致胶原合成增加，氧化还原稳态失衡会干扰线粒体功能障碍进而导致药物引起的肝损伤。ROS 可影响过氧化物酶体增殖物激活的受体调节脂肪酸氧化以及丝裂原活化蛋白激酶（MAPK）和相关的应激敏感激酶活力而调节肝细胞凋亡[14]。LPS 参与了肝损伤中炎症或免疫破坏，可通过激活TLR-2和TLR-4 途径促进炎症活动，并且LPS 诱导的短暂缺血破坏了呼吸链复合物的活性，增强了ROS，从而进一步增强肠道高通透性，并能显著诱导肝损伤[15-16]。

KEGG 富集分析结果表明，补骨脂所致肝损伤可能涉及AMPK、PPAR、CYP 等多条信号通路。AMPK 是NF-κB 的上游蛋白，被认为是维持ADP/AMP/ATP 水平的关键细胞代谢传感器，AMPK 信号通路的激活可抑制非典型NF-κB 途径，从而保护肝脏免受损伤[17-18]。PPAR 信号通路在肝损伤和肝纤维化中起着关键作用，在活化的HSCs 中，PPAR-γ活性降低，从而导致HSCs 从维甲酸储存状态向产生ECM 的肌成纤维细胞转化，可促进CCl4诱导的肝损伤和纤维化，然而，抑制PPARG 可阻止Fsp27α mRNAs 的表达，进而减少乙醇诱导小鼠线粒体活性氧生成所导致的肝损伤[19-20]。肝脏CYP 酶在药物代谢和清除中发挥关键作用，例如CYP2E1 可增加活性氧的数量，进而导致氧化应激肝损伤[21-22]。基于上述结果和理论基础，有必要开展补骨脂AMPK、PPAR 等信号通路的作用研究，从而为其肝损伤机制提供理论依据。

综上所述，本研究通过生物信息学技术构建补骨脂药材-成分-靶点-信号通路-肝损伤网络对补骨脂成分所致肝损伤作用机制进行分析预测，研究结果充分表明了补骨脂多成分-多靶点-多通路的肝损伤作用特点，为进一步开展补骨脂致肝损伤作用机制的研究提供了新思路。