基于LC-MS技术的二氢丹参酮Ⅰ抗肝纤维化肝脏代谢组学研究

陶朝阳，朱臻宇，邢心睿，曹 奇，王 辉 （. 西藏军区总医院药物制剂研究中心，西藏 拉萨 85007；. 海军军医大学药学院，上海 00433）

肝纤维化（hepatic fibrosis，HF）是以肝内细胞外基质（extracellular matrix protein，ECM）过度沉积和纤维瘢痕形成为特征的慢性肝病，在世界范围内具有较高的发病率和病死率，持续发展导致肝硬化、肝癌的发生[1]。肝纤维化发病机制复杂，病毒性肝炎、自身免疫性肝炎、非酒精性脂肪性肝病、胆汁淤积性肝病等均可能导致慢性肝炎症，并最终导致肝纤维化[2]。然而，除肝移植外，目前尚无有效治疗肝纤维化的方法，针对肝纤维化早期识别和治疗对于预防相关负面后果至关重要。鉴于肝脏是人体最大的器官和主要代谢枢纽，探讨肝纤维化的代谢特征有望发现新的标志物和治疗靶点。

二氢丹参酮Ⅰ（dihydrotanshinone I，DHI）是中药丹参中的亲脂性成分，被认为是一种潜在的治疗肝纤维化的药物，其肝脏保护作用、抗癌作用、抗流感活性、抗炎作用等生物学功能已多次报道[3-5]。在本课题组前期扶正化瘀方抗肝纤维化机制研究中发现DHI是扶正化瘀胶囊的重要药效物质基础，也在细胞活性实验中证明其能显著抑制细胞活性发挥抗肝纤维化作用[6]。然而，DHI对肝纤维化治疗作用的体内药效及机制尚不明确。本研究利用肝脏代谢组学方法研究肝纤维化密切相关的生物标志物，探索肝纤维化相关病理过程，同时采用DHI进行干预，研究其对肝纤维化的治疗作用及作用机制，为肝纤维化早期诊断、有效治疗提供科学依据。

1 材料

1.1 仪器

METTLER AE240 型电子天平（瑞士梅特勒公司）；FRESCO17台式冷冻离心机 (Thermo Fisher，美国)；DZG-6020真空干燥箱 (上海益恒实验仪器公司)；Agilent 1290 Infinity 液相色谱仪、Agilent 6538 Q-TOF/MS 质谱仪、Micro17高速离心机（Thermo Fisher Scientific，美国）；HSS T3柱（2.1 mm×100 mm，2.5μm）(Waters，美国)。

1.2 试药

DHI（纯度 98%，上海一飞生物科技有限公司）；硫代乙酰胺(thioacetamide，TAA，东京化成工业株式会社)；甲醇、乙腈（均为色谱纯，德国Merck公司），甲酸（色谱纯，ROE scientific INC，美国）；水为实验室制备的超纯水，其他试剂均为分析纯。

1.3 实验动物

SD大鼠，雄性，质量(200～250 g)，共28只，购自中国科学院上海实验动物中心，合格证号：SCXK(沪)2012-0002。饲养于海军军医大学实验动物中心，遵循动物实验的标准操作规范。饲养条件：温度（22±2）℃，相对湿度40%～60%，12 h昼夜交替循环的条件下笼养。

2 方法

2.1 动物模型的构建和治疗

将28只雄性SD大鼠随机分为4组：正常组、肝纤维化模型组、DHI低剂量组和DHI高剂量组，每组7只。24 h适应性饲养后，除正常组外，模型组及不同用药剂量干预组每周3次腹腔注射200 mg/kg TAA，持续给TAA造模8周；同时，自第5周起，按照给药剂量持续给药4周：正常组、模型组，每日给予生理盐水10 ml/kg；DHI低剂量组，每日给予DHI15 mg/kg；DHI高剂量组，每日给予30 mg/kg。

2.2 样品收集

最后一次给药24 h后，采用脊椎脱臼法处死大鼠。快速切除肝脏后，用0 ℃生理盐水冲洗并立即放于液氮中快速冷冻，储存于-80 ℃冰箱直至分析。

2.3 样品制备

将冻存的肝脏组织置于室温自然解冻，取约100 mg肝脏样本于匀浆管中，加入800μl甲醇，在60 Hz下充分匀浆至没有纤维颗粒，4 ℃离心15 min（14 500×g）后取上清液置于1.5 ml的离心管中，氮气吹干。在上述氮气吹干的样品残渣中加入300 μl含内标甲醇涡旋30 s（内标为L-2-氯苯丙氨酸，浓度为5 μg/ml），每个样品取10μl混匀作为质量控制样品。

2.4 LC-MS分析

色谱条件：Agilent 1 290 Infinity UHPLC，色谱柱：Waters XSelect HSS T3 column色谱柱，柱温：30 ℃；进样量：3 μl；流动相A：含0.1 %甲酸的水，流动相B：含0.1 %甲酸的乙腈。流速：0.4 ml/min；梯度洗脱条件：0～2 min，2 %B，2～17 min，2 %～98 % B，17～19 min，98 %B。质谱条件：Agilent 6530 Accurate-Mass QTOF/MS，离子源：电喷雾（ESI）离子源，正、负离子检测模式；干燥气温度：350 ℃，干燥气体流量：11 L/min；碎裂电压：120 V；毛细管电压：4 000 V（ESI+）/3 500 V（ESI-）；质谱扫描范围：50～1 500m/z。

2.5 数据分析

2.5.1 数据处理

将采集的质谱数据转换为mzData格式文件，然后通过R软件XCMS程序将质谱数据转化为含有保留时间、质核比、峰强度的数据矩阵。保留频数超过80%的质核比数据，并对所有峰面积以内标峰面积和肝组织质量进行归一化处理。

2.5.2 差异代谢物筛选

将上述获得的二维矩阵列表导入SIMCA 14.1软件中进行多元统计分析，采用正交偏最小二乘判别分析研究各组间差异，并得到变量VIP值。正常组和模型组组间比较采用独立样本t检验，采用VIP＞1，且选择差异具有显著性（P＜0.05）的变量作为潜在的差异代谢物，然后检索数据库（HMDB、METLIN、KEGG数据库），筛选出潜在的差异代谢物。

2.5.3 统计学分析

使用SPSS 21.0统计软件进行数据统计分析，组间数据采用单因素方差分析（ANOVA），两组样本分析采用独立样本t检验，分析数据差异的统计学意义，以P＜0.05为差异有统计学意义。

3 结果

3.1 基于LC-MS的肝组织代谢轮廓分析

本研究采用UPLC-Q-TOF/MS在正、负两种模式下对肝组织进行代谢组学分析，两种模式下典型的总离子流图见图1。代谢组学数据的稳定性、重现性对研究的可靠性非常重要，为考察实验的系统稳定性，本研究从每个肝组织样品中取10 μl混匀后作为质量控制样品。在样品序列一开始连续进样10针QC样品，并按每7个样品再进样一针，共进样14次。正、负离子模式下，QC样品均显示良好的聚集状况，证明该分析系统稳定可靠。

图1 肝脏正负离子模式下典型的LC-MS总离子流图

3.2 多变量统计分析

经过数据处理后，LC-MS数据集在正离子模式下得到368个离子，在负离子模式下得到249个离子。根据OPLS-DA方法对数据进行分析，OPLSDA是一种多因变量对多自变量的回归建模方法，最大特点是可以去除自变量和分类变量无关的分类变异，根据得分图可以揭示数据离散程度，发现异常值[7]。在得分图上具有相似的代谢物组成的样本，处在比较相似的位置，样本间距越远代表代谢物差异越大，样本间生理状态相差越大。将正常组、模型组、DHI低剂量组和DHI高剂量组的LCMS数据导入到SIMCA进行分析，其得分如图2A、2B所示，正常组与模型组样本各自聚为一类，并且完全分离，说明肝纤维化大鼠模型的肝组织代谢轮廓发生显著改变，代谢物的种类或水平发生了明显变化。DHI给药组（DHI低剂量组、DHI高剂量组）能够与模型组明显区分，并向正常组靠近，表明DHI给药组可恢复部分代谢物至正常水平。

图2 基于LC-MS正负离子模式下的OPLS-DA分析得分图

3.3 差异代谢物的筛选

OPLS-DA可用于寻找导致聚类间显著差异的变量，筛选正常组和肝纤维化大鼠模型组间潜在的差异代谢物。基于OPLS-DA模式下VIP＞1筛选两组间的差异代谢物，进行t检验后，筛出模型组和正常组间差异具有显著性的变量（P＜0.05）。基于VIP和t检验以及HMDB和KEGG等数据库比对，筛选出38个与TAA诱导肝纤维化大鼠模型相关的较为重要的差异代谢物作为潜在生物标志物（表1）。通过KEGG、HMDB等数据库查询，我们发现这些代谢物主要涉及谷胱甘肽代谢、褪黑素代谢、氨基酸代谢、脂质代谢、三羧酸循环等途径。

表1 TAA诱导肝纤维化相关的差异代谢物及其代谢通路

3.4 DHI的干预治疗

以肝纤维化相关的差异代谢物的相对含量作为检测指标可以评价DHI对肝纤维化的治疗作用，比较发现33种代谢物的含量发生明显逆转（图3），同时，有14种代谢物的含量回调与DHI的剂量成正相关性，即高剂量组比低剂量组回调更多（图3A、3B）。

图3 33种潜在标志物在正常组、肝纤维化模型组、DHI给药组中的相对含量变化

4 讨论

4.1 研究意义

肝纤维化是一种发病机制复杂，可导致肝硬化、肝衰竭甚至肝癌等严重肝脏疾病的慢性流行肝病，然而目前并没有药物可以治疗肝纤维化，中药及其活性成分由于其多靶点、多通路等特点被认为是肝纤维化治疗的潜在药物。代谢组学可以对生物体内小分子代谢产物进行动态分析，通过分析阐述代谢物与生理病理变化间的联系，对肝纤维化的机制进一步阐述，也能通过肝纤维化相关差异代谢物的含量相对变化分析药物对肝纤维化模型的药效作用。本研究用肝脏代谢组学方法分析了与TAA诱导的肝纤维化模型密切相关的38种代谢物，主要涉及谷胱甘肽代谢、褪黑素代谢、氨基酸代谢、脂质代谢、三羧酸循环等途径，DHI能够通过调节部分代谢通路发挥预防和治疗肝纤维化作用。

4.2 谷胱甘肽代谢

肝纤维化的产生伴随着肝细胞死亡和肝星状细胞(HSCs)的活化，HSCs约占正常人肝脏中非实质细胞的1/3和总驻留细胞的15%，HSCs从静态到激活会使得ECM过剩表达而导致纤维化，而氧化应激在HSCs激活和ECM形成中发挥重要的促进作用[1,8]。谷胱甘肽（GSH）是人体内含量最丰富的抗氧化剂，也是体内氧化防御体系的主要组成成分之一[9]。GSH可以在谷胱甘肽过氧化物酶和谷胱甘肽-S-转移酶以及谷胱甘肽还原酶的作用下与其氧化态相互转换，清除部分有机过氧化物，调节体内的氧化还原稳态，缓解氧化应激对组织造成的损伤[10-11]。有研究表明，在对乙酰氨基酚、四氯化碳、重金属砷等诱导下，大鼠肝组织中GSH含量下降[12]。本研究中，TAA诱导的肝纤维化大鼠模型肝组织中GSH显著下降，这与谷胱甘肽在其他肝损模型中下调的趋势一致。经过DHI干预后GSH回调，且高剂量DHI比低剂量组回调比例更大，这说明DHI可能调整体内谷胱甘肽代谢通路，通过提高体内谷胱甘肽含量恢复肝组织抗氧化功能，缓解氧化应激对肝脏的进一步损伤，从而起到抗肝纤维化作用。

4.3 褪黑素代谢

肝纤维化过程中，抑制HSC的激活和增殖是预防和治疗肝纤维的重要途径。血小板衍生生长因子(PDGF)可以激活JAK2/STAT3信号通路，致使HSC增殖，并抑制HSC凋亡[13]。已有研究证明，褪黑素通过抑制JAK2/STAT3相关信号通路或抗氧化机制来抑制HSC的激活和增殖从而发挥肝脏保护作用[14-15]。在本研究中，环6-羟基褪黑素及N-乙酰血清素硫酸盐均为褪黑素代谢产物，在模型组中含量显著增加，说明褪黑素被大量代谢，肝脏保护作用被抑制。经过DHI的干预后两种褪黑素代谢产物均回调，这提示DHI可能通过调整褪黑素代谢，回调褪黑素及其代谢产物的机体内含量发挥抗肝纤维化活性。

4.4 氨基酸代谢

肝脏是机体物质代谢的中枢器官，在氨基酸的新陈代谢和蛋白质的合成与分解中发挥重要作用。天冬氨酸是一种酸性氨基酸。天冬氨酸在哺乳动物中一般被认为是一种营养上非必需的氨基酸。然而，越来越多的文献表明，天冬氨酸在包括肝脏生理学在内的许多生物和生理过程中起着重要作用，如合成精氨酸以维持巨噬细胞应对免疫挑战，同时，也有证据表明天冬氨酸可以减轻肝损伤、增强肝脏功能[16]。本研究发现TAA诱导大鼠肝纤维化后，与正常组比较，模型组中丙氨酸、天冬氨酸和谷氨酸代谢通路被干扰，L-天冬氨酸含量显著变化，经过DHI的干预后回调，这提示DHI可能调整丙氨酸、天冬氨酸和谷氨酸代谢通路，改变体内天冬氨酸含量发挥肝保护作用。

4.5 脂质代谢

脂质不仅是细胞膜的组成成分，而且参与信号转导。磷脂酰胆碱、鞘脂、溶血性磷脂酰胆碱是血清脂蛋白和细胞膜的重要组成成分。溶血磷脂酰乙醇胺是磷脂酶A1水解磷脂酰乙醇胺(PE)失去一分子脂肪酸生成的产物，作为一种溶血卵磷脂，在肝脏中与卵磷脂一起在线粒体间进行转移[17]。有研究证明，PE通过N-甲基转移酶生成PC的通路约占PC产生量的30%，磷脂酰胆碱为肝脏重要营养来源，可以拮抗肝脏受病毒、药物、酒精及其他有毒物质的侵害，防止肝纤维化[9]。肝脏是脂质生成和脂肪酸氧化等脂类代谢的主要场所，肝纤维化时可以引起脂质合成、转运及分解代谢的普遍紊乱，而脂质代谢异常又会加重肝损害，引起肝脏的脂毒性[18]。在本研究中，4种LysoPE在模型组中均发生下降，另外一些磷脂代谢相关物质，如GlcCer(d18:1/12:0)、PE(16:0/P-16:0)、PE(18:0/15:0)、CL(i-13:0/i-22:0/i-12:0/i-13:0)和脂质代谢相关物质如二酰甘油（DG）、甘油磷酸甘油、牛磺酸熊去氧胆酸等在正常组与模型组大鼠中也存在显著差异，说明TAA诱导的大鼠肝纤维化模型脂质代谢的紊乱。经过DHI干预治疗后，代谢物水平不同程度得到回调，表明DHI可通过干预多个脂质代谢通路发挥其抗肝纤维化活性。

4.6 能量代谢

三羧酸循环是糖、脂肪、氨基酸三大营养素的最终代谢通路，在能量代谢、提供生物合成的前体中起重要作用。琥珀酸既是三羧酸循环的的重要中间产物，也是一种重要的细胞外信号分子，在信号传递、炎症反应、肝纤维化的发生发展中发挥着重要作用[19]。在本研究中，肝纤维化模型组大鼠琥珀酸的代谢变化趋势均与正常大鼠相反，而给予DHI后趋近正常，说明DHI能够改善肝纤维化大鼠的能量代谢。

5 结论

本研究建立了LC-MS代谢组学分析的方法，通过OPLS-DA筛选出与肝纤维化密切相关的38种差异代谢产物，主要涉及谷胱甘肽代谢、褪黑素代谢、氨基酸代谢、脂质代谢、三羧酸循环等途径。基于38种差异代谢物在正常组、模型组及DHI给药组之间相对含量的差异，对DHI预防治疗大鼠肝纤维化模型的药效进行了评价，结果表明DHI能够回调34种潜在差异代谢物的相对含量，调节部分失衡代谢通路，发挥抗肝纤维化作用。本研究为肝纤维化的新标志物和治疗靶点的发现以及DHI作为抗肝纤维化潜在药物的进一步开发应用提供实验依据。