大鲵活性肽对D-半乳糖致小鼠机体氧化损伤的修复作用

关百婷，李 伟，赵 菲，佟长青

（大连海洋大学食品科学与工程学院，辽宁大连 116023）

大鲵（Andrias davidianus），属于两栖纲、有尾目、隐鳃鲵科，是我国特有的珍稀两栖动物，属于国家二类保护动物[1]。大鲵不仅有丰富的营养，而且还体虚、贫血、痢疾等疾病[3]。研究表明，大鲵含有丰富的蛋白质、不饱和脂肪酸、多糖、矿物质元素等成分[4]。这些营养物质构成了大鲵活性成分的物质基础。大鲵肌肉中蛋白质含量高，氨基酸比例合理，种类丰富。大鲵肌肉蛋白质经过酶解过程可以获得具有抗氧化、免疫调节等生物活性的大鲵活性肽[5−7]，适用于大鲵肽相关产品的开发。

衰老是多细胞生物体在发育成熟后受多因素作用的生理机能退行性过程，主要有自由基、免疫、端粒细胞与神经内分泌等学说解释机体的衰老原因，但主流的解释为自由基攻击等氧化反应导致机体衰老[8]。因机体内脏器受氧化物攻击，产生有毒有害物质无法及时清除，产生一系列生理生化应激反应，进而导致表观上、生理上的退行性衰老。研究表明，长期大量对动物注射D-半乳糖溶液，过多的D-半乳糖在醛糖还原酶的作用下生成半乳糖醇，半乳糖醇不能被机体代谢，堆积在细胞内的半乳糖醇可以导致细胞的氧化应激反应，进而导致细胞内活性氧（reactive oxygen species,ROS）增加[9]。过多的ROS引起细胞损害，导致机体多器官、多系统功能减退[10−11]。D-半乳糖可以使机体的很多组织出现类似衰老的退行性改变，出现与衰老相关的毛色变暗、行动迟缓等体征、生化改变、免疫功能下降、退行性神经功能紊乱，还出现胰腺损伤、脑损伤以及肾脏损伤[2−15]。因此，D-半乳糖致小鼠衰老模型被广泛用来评价各种生物活性成分的抗衰老作用。

目前，对大鲵肌肉的营养成分进行分析的研究报道较多。对大鲵肽的酶解制备工艺及其具有抗氧化活性的研究也有报道[16−17]。但鲜有对大鲵肽的抗衰老作用的深入研究。本实验采用D-半乳糖皮下注射建立小鼠亚急性衰老模型以及非定量标记蛋白质组学技术研究大鲵肽的抗衰老作用，以期从分子水平上阐明大鲵活性肽抗衰老机理，为大鲵的高值化开发利用提供科学依据。

1 材料与方法

1.1 材料与仪器

雄性SPF小鼠 体重18～22 g，大连医科大学实验动物中心，生产许可证SCXK（辽）2013-0003；大鲵活性肽 分子量1400～2000 Da，由辽宁省水产品加工及综合利用重点实验室制备；D-半乳糖 国药集团化学试剂有限公司；总抗氧化能力（Total antioxidant capacity，T-AOC）、超氧化物歧化酶（Superoxide dismutase，SOD）、丙二醛（Malondialdehyde，MDA）测定试剂盒 南京建成生物工程研究所。

5430R低温高速离心机、真空离心浓缩仪 Eppendorf公司；Q exactive质谱仪Thermo Scientific公司；Triple TOF 6600+质谱仪AB SCIEX公司；Agilent 1290 Infinity LC超高压液相色谱仪 Agilent公司。

1.2 实验方法

1.2.1 D-半乳糖致小鼠衰老模型的建立、给药和处理 将60只SPF级昆明小鼠（雄性）分为5组，分别为正常组、衰老模型组，大鲵活性肽低、中、高剂量组，对应剂量分别为：100、200、400 mg/（kg·d）。除正常组每天注射等体积生理盐水外，其他各组每天都单次皮下组织注射5%D-半乳糖0.5 mL/d（生理盐水配制）。每周记录1次小鼠体重。大鲵活性肽低、中、高剂量组每日按0.1 mL大鲵活性肽水溶液/20 g小鼠体重分别灌胃一次，正常组和衰老模型组以生理盐水代替大鲵活性肽水溶液进行灌胃，同时观察小鼠毛色与行动。连续灌胃6周后，称重，摘取眼球取血，立即离心（3000 r/min，10 min），取上层血清，−20℃冰箱冻存备用。

1.2.2 小鼠组织脏器指数测定 小鼠称重后处死，快速取出脑、肝脏、胸腺和脾4种脏器，用滤纸吸干血水，电子天平上称重。按照式（1）计算肝脏指数、脑指数、胸腺指数和脾指数。

1.2.3 生化指标测定 小鼠血清中总抗氧化（T-AOC）能力、超氧化物歧化酶（SOD）活性和丙二醛（MDA）含量的测定，按试剂盒说明书的操作步骤进行。

1.2.4 蛋白质组样品制备 取50μL大鲵活性肽中剂量组和D-半乳糖致小鼠衰老模型组小鼠血清，均采用去血清高丰度亲和色谱柱Agilent Multiple Affinity Removal LC Column-Mouse 3（4.6 mm×100 mm，1.66 mL），获得小鼠血清中低丰度蛋白组分。各样品血清低丰度蛋白组分经过超滤浓缩后加入等体积SDT裂解液，沸水浴裂解15 min，于14000×g条件下离心20 min。采用BCA法对上清液进行定量。大鲵活性肽中剂量组小鼠血清样品标记为ADB-T，D-半乳糖致小鼠衰老模型组小鼠血清样品标记为ADB-M。

1.2.5 SDS-PAGE电泳 取裂解后的各样品分别取20μg，分别加入等体积2×上样缓冲液，沸水浴5 min，取20μL进行12.5%SDS-PAGE电泳后，考马斯亮蓝染色后脱色液脱色，判断样品裂解情况。

1.2.6 样品酶解 采用Filter aided proteome preparation （FASP）方法对各样品进行酶解[18]。分别取样品ADB-T及ADB-M各30μL蛋白质溶液，分别加入0.003 mmol DTT，沸水浴5 min后冷却至室温。加入200μL尿素缓冲液（8 mol/L尿素，150 mmol/L Tris-HCl，p H8.0）混匀后，在10 k D超滤离心管中于14000×g离心15 min，弃滤液（重复该步骤一次）。加入100μL 100 mmol/L IAA缓冲液（以尿素缓冲液溶解），在600 r/min条件下振荡1 min，室温避光30 min后于相同离心条件下进行离心。加入100μL尿素缓冲液，相同离心条件下离心，该步骤重复两次。加入100μL 25 mmol/L NH4HCO3溶液，相同离心条件下离心两次。加入40μL 0.1 g/L胰蛋白酶缓冲液（以100 mmol/L NH4HCO3溶液配制），600 r/min振荡反应1 min，37℃静置16 h。换新收集管，相同离心条件下离心；再加入40 μL 25 mmol/L NH4HCO3，在相同离心条件下离心并收集滤液。通过C18Cartridge脱盐冻干后以40 μL 0.1%甲酸溶液复溶，测定280 nm吸光值[19]。

1.2.7 LC-MS/MS分析 根据FASP酶解后OD280定量结果，各取2μg酶解后的ADB-T及ADB-M进行LC-MS/MS分析。样品先经过RP-C18毛细管高效液相色谱柱（Thermo EASY column SC200 150 μm×100 mm）的纳升流速HPLC液相系统EASY-nLC1000，以甲酸乙腈水溶液进行梯度分离。预柱为RP-C18Thermo EASY column SC001 traps 150μm×20 mm。A液为0.1%甲酸乙腈水溶液（乙腈为2%），B液为0.1%甲酸乙腈水溶液（乙腈为84%）。梯度条件为：0～100 min，B液线性梯度0%～45%；100～108 min，B液线性梯度45%～100%；108～120 min，B液维持在100%[19]。分离后用Q-Exactive质谱仪进行质谱分析[15]。预分析时长60 min，正式分析时长120 min。

1.2.8 Max Quant的数据分析 质谱分析原始数据为RAW文件，采用Max Quant软件（版本号1.3.0.5）[20−21]进行查库鉴定，并根据Label free算法进行定量分析[22]。

1.2.9 Gene Ontology（GO）功能注释 利用Blast-2GO[23]对目标蛋白质集合进行GO注释。利用NCBI BLAST+ （ncbi-blast-2.2.28+-win32.exe）比对小鼠血清数据库，保留E值小于等于0.001的前10条。再依次通过提取、注释。

1.2.10 KEGG通路注释KEGG Automatic Annotation Server软件[24]获取目标蛋白质序列参与的通路信息。通过Fisher精确检验评价GO term或KEGG通路蛋白质富集度的显著性水平。

1.2.11 蛋白质聚类分析（Clustering） 对目标蛋白质定量信息归一化处理后，进行蛋白质聚类分析。

1.3 数据处理

数据处理采用SPSS 19.0软件进行单因素方差分析，结果表示为平均值±标准差，P<0.05为差异显著。

2 结果与分析

2.1 小鼠体质变化情况及脏器指数

各组小鼠在实验前，皮毛光滑，喜攀爬、啮咬树枝。衰老模型组小鼠连续6周每天皮下组织注射5%D-半乳糖0.5 mL后，小鼠毛色变暗，行动迟缓，攀爬行为减少等变化。该结果与黄杰等[13]建立D-半乳糖致衰老小鼠的反应与状态变化结果相似。衰老模型组小鼠体重仅增加56.18%，正常小鼠体重增60.30%，增加量低于正常组（图1）。灌胃低、中、高浓度大鲵活性肽水溶液各组小鼠平均体重高于模型组，趋近于正常组，但未表现出剂量依赖性，表明大鲵活性肽在低剂量情况下具有促进衰老小鼠恢复生长作用。

图1 大鲵活性肽对小鼠体质量的影响Fig.1 Effect of giant salamander bioactive peptideson body weight of mice

各组SPF级雄性昆明小鼠灌胃对应大鲵活性肽溶液6周后，小鼠脏器指数变化如表1。由表1可以看出，各灌胃大鲵活性肽组的脑指数、胸腺指数与正常组相比不显著（P>0.05）；中剂量组的脾指数显著高于正常组（P<0.05）；低、高剂量大鲵活性肽组的肝指数均极显著低于模型组（P<0.01）。灌胃大鲵肽后，小鼠的脾脏指数增加，表明大鲵肽在一定程度上可以增强小鼠机体免疫功能。

表1 大鲵活性肽对小鼠脏器指数的影响Table 1 Effect of giant salamander bioactive peptides on organ indexes for each group of mice

2.2 生化指标检测

对小鼠血清中MDA含量、SOD活力和T-AOC进行测定，结果见表2。由表2可见，与模型组相比，灌胃大鲵活性肽中、高剂量组小鼠血清中MDA含量极显著降低（P<0.01），并且SOD活力和T-AOC显著升高（P<0.05）。注射D-半乳糖致小鼠衰老后，小鼠血清中MDA含量极显著上升，SOD活力和T-AOC极显著降低（P<0.01），而在灌胃大鲵活性肽后，显著改善小鼠这些生化指标，但未发现剂量依赖性。因此，大鲵活性肽具有抗D-半乳糖致小鼠衰老作用。

表2 大鲵活性肽对小鼠生化指标的影响Table 2 Effect of giant salamander bioactive peptides on the levels of serum markers in SPFmice

2.3 血清样品SDS-PAGE和质谱检测

SDS-PAGE（图2A）结果表明，未经去除血清高丰度亲和色谱柱处理的中剂量组小鼠血清蛋白（ADB-T2-3原）中60 kDa附近蛋白含量过高而影响低丰度蛋白质分析，而去除高丰度血清蛋白后模型组（ADB-M-1、ADB-M-2、ADB-M-3）和的大鲵活性肽中剂量组（ADB-T2-1、ADB-T2-2、ADB-T2-3）电泳条带清晰，平行性好，可进行质谱分析。质谱分析Basepeak结果（图2B和图2C）表明，组间样本平行性好，酶解效果良好，可进行非标记定量质谱测定。

图2 小鼠血清样本平行性Fig.2 Samples parallelism of proteinsin serum in mice

2.4 蛋白质鉴定结果统计

模型组与中剂量组共鉴定到蛋白质425个，其中差异表达蛋白26个，通过大鲵活性肽治疗后，差异表达蛋白质表达上调13个（表3），下调13个（表4）。其中胞质四氢叶酸合成酶（存取编号：Q922D8）和肌糖原磷酸化酶显著增加，黄嘌呤脱氢酶（氧化酶）（存取编号：Q00519）显著减少。四氢叶酸作为生物体中的一碳单位转移酶的辅酶，在氨基酸代谢和核酸代谢中起重要作用。糖原磷酸化酶催化糖原转变为1-磷酸葡萄糖，补充体液中葡萄糖。黄嘌呤脱氢酶在核酸代谢过程中起重要作用，可催化黄嘌呤转变为尿酸。

表3 上调的差异表达蛋白质Table3 Up-regulated proteinswith significant differences

表4 下调的差异表达蛋白质Table 4 Down-regulated proteins with significant differences

以治疗组对模型组表达倍数与t检验P值作图，显示两组间蛋白质表达差异显著性的火山图如图3所示，图中横坐标大于1且纵坐标大于1.3的区域为蛋白质表达显著上调区域，横坐标小于−1且纵坐标大于1.3的区域为蛋白质表达显著下调区域。从图3中可以看出，13个表达上调的蛋白质中，最大表达倍数为4倍多（log24.49），而表达下调的13个蛋白质中有3个表达下调量极大，分别为0.067、0.021、0.004倍。

图3 蛋白质表达差异显著性的火山图Fig.3 Volcano plot of proteins with significant differences

对26个差异表达蛋白质进行聚类分析，结果如图4所示。这些差异表达蛋白质主要分为两类，即表达量上调和下调。在MCL-T处理组的上调表达的蛋白质中，胞质四氢叶酸合成酶在各组中的差异较大，在下调表达的蛋白质中，热休克蛋白HSP 90-β亚基、黄嘌呤脱氢酶和碳酸酐酶在各组中表达量差异较大。

图4 差异表达蛋白质聚类结果Fig.4 Cluster of proteins with significant differences

通过Blast2Go（https://www.blast2go.com/）软件对鉴定到的425种蛋白质进行基因本体（Gene ontology，GO）功能注释，2水平注释结果如图5所示。从图5中可以看出，小鼠低丰度蛋白参与了14个生理过程、具有6种分子功能、存在于7种细胞组分中。在这些生理过程中，代谢过程、细胞生成过程、多细胞组织过程、刺激应答、生理调节等过程都有10个以上差异表达蛋白；这些差异表达蛋白中，有15个具有催化活性、15个具有结合功能、6个具有分子功能调节活性；其中胞外区有21个、细胞中有13个、高分子复合物有14个、细胞器中有10个、胞外区部分有18个。

图5 差异表达蛋白基因本体2水平统计Fig.5 GOlevel 2 statistics

对参与的生理过程、具有的分子功能和所处的细胞组分对小鼠血清低丰度表达蛋白质进行分类，再通过Fisher精确检验方法对小鼠血清低丰度表达的26个差异表达蛋白进行GO功能富集分析，结果如图6和表5所示。结果表明，在小鼠血清低丰度表达蛋白质中，与ADB-M组相比，ADB-T组有5、3、2、2种显著表达的蛋白质分别参与了细胞趋化性、凋亡过程中半胱氨酸型肽链内切酶活性、转化生长因子β受体信号通路正向调节、嘌呤碱基代谢过程等4个生物过程，其富集度分别为0.227、0.333、0.667、0.667；而具有对C-N键（非肽键）作用的水解酶活性、丝氨酸水解酶活性、趋化活性等3个生物功能的蛋白质分别有10、9、2种表达显著；存在于微管和胞浆组分中的差异表达蛋白各有3个，其富集度分别为0.6和0.375。

表5 GO功能富集度Table 5 Richfactor of enriched GOterms

图6 差异表达蛋白GO功能富集结果Fig.6 Enriched GO termsof proteinswith significant differences

生物体内蛋白质在行使其功能时，受多方面因素调节，包括与其它蛋白质之间的协调作用。这些相互作用构成代谢通路。通过观察大鲵活性肽治疗前后蛋白质在KEGG通路中表达的变化，可以在蛋白质层面获得大鲵活性肽对治疗小鼠衰老的机理。如图7所示，通过Fisher精确检验方法对差异表达蛋白质进行KEGG通路富集分析，结果表明，内分泌调节的及其它因素调节的钙再吸收，肾素-血管紧张素系统，IL-17信号通路和坏死性凋亡等重要通路发生了显著变化。

图7 KEGG通路富集结果Fig.7 Enriched KEGG pathways

3 讨论

衰老是多细胞生命体生长发育成熟后在遗传、环境与生活方式等多种因素共同作用下引起生理机能逐步退行性下降的过程，小鼠衰老时表观学上为毛色灰暗，体重增长缓慢，精神萎靡等[19]。已有文献表明，机体衰老与体内自由基的量有关[9,13]，抗衰老活性物质一般具有清除自由基活性。具有抗衰老活性的生物活性物质可令代谢紊乱小鼠代谢趋向正常化，其体内相关抗氧化酶和还原性物质含量得到改善，例如抗氧化蛋白[25]、超氧化物歧化酶[26]等，也可以降低衰老小鼠MDA含量。

在本研究中，通过对小鼠连续6周注射D-半乳糖后，模型组与灌胃大鲵活性肽组相比，小鼠血清低丰度蛋白质表达差异的有26个，包括四氢叶酸合成酶、糖原磷酸化酶和黄嘌呤脱氢酶等蛋白质。四氢叶酸合成酶是叶酸代谢途径中的关键酶，在一碳单位转移过程起作用，主要参与核酸代谢、蛋白质与核酸的甲基化与修复反应[27]，为生物进行正常生长代谢、细胞损伤修复提供了必要的物质基础。本研究中衰老小鼠接受大鲵活性肽后，胞质四氢叶酸合成酶表达量上调，有利于D-半乳糖致衰老小鼠恢复正常。糖原磷酸化酶是糖原分解为1-磷酸葡萄糖的限速酶，其量的增加或降低影响生物的供能、高分子化合物合成代谢重要中间化合物来源[28]。本研究中的衰老小鼠经大鲵活性肽的调理后，糖原磷酸化酶表达量上升，有利于衰老小鼠恢复血糖浓度，进入正常生长状态。黄嘌呤脱氢酶参与嘌呤代谢，催化产生高氧化能力的过氧化物，使得机体易受氧化损伤[29−31]。灌胃大鲵活性肽后，衰老小鼠黄嘌呤脱氢酶表达量降低，在一定程度上可以减少过氧化物的含量，使得衰老小鼠脏器免受过氧化物的攻击。因此，对D-半乳糖致衰小鼠灌胃大鲵生物活性肽，可以改善小鼠体内TAOC、SOD和MDA，并使多个代谢关键酶的表达量增加或减少，改善小鼠衰老状态，因此本研究所用大鲵生物活性肽具有抗衰老活性。

4 结论

本研究通过D-半乳糖致小鼠衰老实验，发现大鲵活性肽可以使衰老小鼠体重和脾指数恢复增加，并引起SOD和T-AOC升高，降低肝指数和MDA含量，增强机体免疫，改善衰老小鼠紊乱的生化指标。衰老小鼠血清低丰度蛋白质中胞质四氢叶酸合成酶和肌糖原磷酸化酶表达量增加，黄嘌呤脱氢酶表达量减少，这些蛋白质表达量的变化，可以降低机体中自由基，延缓衰老进程。这些结果表明，分量为1400～2000 Da的大鲵活性肽具有抗小鼠衰老作用。