羽扇豆种子蛋白源多肽对精氨酸酶的抑制动力学与分子模拟研究

《Journal of Functional Foods》：Kinetics and molecular docking studies of the inhibition of arginase by lupin seed protein-derived peptides

　　

精氨酸酶（arginase）功能失调与多种病理状态相关，包括心血管疾病、勃起功能障碍和癌症，因此成为极具吸引力的治疗靶点。在本研究中，研究人员评估了白羽扇豆（Lupinus albus）和蓝羽扇豆（Lupinus angustifolius）分离蛋白的酶解产物对精氨酸酶的抑制活性。在所测试的水解产物中，白羽扇豆风味蛋白酶（flavourzyme）水解产物表现出最高的抑制效力（在800 μg/mL下约为35%），其次为胰蛋白酶和胃蛋白酶水解产物，而碱性蛋白酶（alcalase）水解产物未显示活性。研究人员通过反相高效液相色谱（RP-HPLC）分级分离并进行串联质谱分析，从最具活性的级分中鉴定出十个肽序列。两个结构不同的肽，即短疏水六肽DMGAVA和长聚阳离子十六肽VVALKKVGVLKKPGPT，表现出最高的精氨酸酶抑制作用，被选中用于机制研究。Lineweaver-Burk分析揭示了混合型抑制动力学，且具有不同的特征：DMGAVA显著影响K_m和V_max（K_i约为266 μM），而VVALKKVGVLKKPGPT主要表现为竞争性行为，V_max变化极小（K_i约为589 μM）。

精氨酸酶（arginase）是含锰金属酶，催化L-精氨酸水解生成尿素和L-鸟氨酸，在氨解毒、多胺与脯氨酸合成以及调控一氧化氮（NO）生物利用度中发挥关键作用。精氨酸酶功能异常与心血管疾病（如高血压、动脉粥样硬化）、勃起功能障碍、肺动脉疾病及癌症等密切相关。尽管已开发出nor-NOHA、ABH、BEC、CB-1158等合成抑制剂，但这些化合物普遍存在口服生物利用度低、代谢降解快以及缺乏亚型选择性等药代动力学缺陷，且进入临床研究的品种极少，进展缓慢。因此，寻找来源天然、安全性高的精氨酸酶抑制剂成为研究热点，特别是食源性生物活性肽。羽扇豆种子富含优质蛋白，其酶解产物已被证实具有抑制血管紧张素转换酶（ACE）、二肽基肽酶-IV（DPP-IV）等活性，但关于其精氨酸酶抑制肽的报道尚属空白。基于此，研究人员系统评估了白羽扇豆与蓝羽扇豆蛋白水解物的精氨酸酶抑制活性，并利用色谱、质谱、动力学及计算模拟手段鉴定和解析了活性肽的作用机制。研究表明，羽扇豆蛋白水解物能够释放结构多样的精氨酸酶抑制肽，其中疏水短肽和聚阳离子长肽均可通过互补性机制发挥抑制作用，为开发精氨酸酶相关疾病的膳食干预策略提供了新思路和新分子模板。

为完成该研究，主要采用了以下关键技术方法：以来自The Lupin Platform Inc.的白羽扇豆（Lupinus albus ‘Dieta’）和蓝羽扇豆（Lupinus angustifolius ‘Boregine’）种子为原料，分别制备蛋白分离物；选用风味蛋白酶、胰蛋白酶、胃蛋白酶、胰凝乳蛋白酶、碱性蛋白酶及胃蛋白酶-胰酶顺序消化进行酶解；采用反相高效液相色谱（RP-HPLC）对水解物进行两步分级纯化；利用串联质谱（MS/MS）鉴定最具活性级分中的肽序列；合成候选肽后，通过基于尿素的比色法测定精氨酸酶抑制活性，并进行Lineweaver-Burk动力学分析；借助HPEPDOCK进行分子对接、用PRODIGY计算结合亲和力；采用GROMACS和CHARMM36力场进行100 ns分子动力学（MD）模拟，并用MMGBSA方法计算结合自由能。

3.1 羽扇豆蛋白水解物的精氨酸酶抑制活性
研究人员用不同蛋白酶水解两种羽扇豆分离蛋白，在400–800 μg/mL浓度范围内评估水解物的抑制活性。结果显示，所有活性水解物呈剂量依赖性抑制；白羽扇豆风味蛋白酶水解物（WFLA）活性最高（800 μg/mL时约35%），其次为白羽扇豆胰蛋白酶水解物（WTRP，约32%），胃蛋白酶水解物活性中等，而碱性蛋白酶及胃蛋白酶-胰酶顺序消化产物无可检测活性。白羽扇豆水解物活性普遍高于蓝羽扇豆，且水解度（DH）与抑制活性无直接相关性，表明肽的组成类型比水解程度更关键。该结果验证了蛋白酶选择对生成活性肽的决定性影响。

3.2 HPLC分级和活性肽鉴定
为富集活性肽，研究人员对活性较高的白羽扇豆胃蛋白酶水解物进行两轮RP-HPLC分级：第一轮得到F1–F7七个级分，其中亲水性最强的F1活性最高；第二轮将F1进一步分为F1–1至F1–7七个亚级分，其中F1–1、F1–5、F1–6和F1–7活性较强。对这些亚级分进行串联质谱鉴定，共得到10个序列，长度6–16个氨基酸，平均局部置信度（ALC）为55%–63%，涵盖疏水、碱性及富含脯氨酸等不同类型肽段，表明精氨酸酶抑制可能经由多种相互作用模式实现。

3.3 合成肽的精氨酸酶抑制活性
研究人员固相合成上述10个肽并测定活性，发现短疏水六肽DMGAVA和长聚阳离子十六肽VVALKKVGVLKKPGPT活性最强，二者结构截然不同：前者由小分子疏水残基和一个天冬氨酸组成，分子量约562 Da；后者含四个赖氨酸，净正电荷高，分子量约1635 Da。这种结构差异提示精氨酸酶可被不同理化特征的肽通过不同方式调节。与已报道的基于精氨酸酶I C端S型基序设计的合成五肽REGNH（IC₅₀约2.4 mM）相比，本文肽在更低浓度范围即表现出抑制作用，且由筛选而非理性设计得到。

3.4 精氨酸酶抑制动力学
通过不同L-精氨酸浓度下的Lineweaver-Burk双倒数作图分析动力学参数。DMGAVA使K_m由36.74 mM升至50.53 mM（100 μM肽时增加37.5%），同时V_max降低约15%，呈典型混合型抑制，K_i约266 μM，提示其同时结合游离酶和酶-底物复合物。VVALKKVGVLKKPGPT使K_m增至42.98 mM（增加17.0%），V_max几乎不变，主要表现为竞争性抑制，K_i约589 μM，说明其可能结合于活性位点入口或外周阴离子区域，空间阻碍底物进入而不影响催化效率。与已知小分子抑制剂（如CB-1158 K_i 4.1 nM）相比，这两个肽的K_i高出4–5个数量级，表明其更适合作为先导序列而非药理学竞争性抑制剂。

3.5 分子对接分析
分子对接显示，DMGAVA深入精氨酸酶活性位点，与S2底物结合环上的ASP128、ASN130及催化中心附近的ASP181、ASP183和参与底物定位的THR246形成广泛氢键，并与PRO20存在疏水接触，其直接接触催化残基解释了V_max的降低。VVALKKVGVLKKPGPT则主要结合于活性位点入口边缘，与GLU151、LEU152、TYR188、TYR197等外周残基形成氢键，与THR131、ILE189等形成疏水接触，仅与ASP181有共同锚定，该外周结合模式与竞争性抑制为主的动力学特征一致。

3.6 分子动力学（MD）模拟
100 ns MD模拟评估复合物稳定性。两种肽-酶复合物的回旋半径（Rg）均稳定，未发生显著整体构象变化。DMGAVA复合物在15–20 ns内快速达到平衡，RMSD稳定在0.15–0.20 nm，且活性位点区域残基的RMSF降低，说明其稳定了活性位点构象。VVALKKVGVLKKPGPT复合物则需30–40 ns才达平衡，RMSD在0.20–0.25 nm区间波动，呈现更广泛的残基动力学改变，但活性位点内部残基灵活性相对较高，与其外周结合、不深入催化口袋的模型一致。

3.7 结合自由能分析（MMGBSA）
MMGBSA计算表明，DMGAVA-精氨酸酶复合物的总结合自由能为?35.49 kcal/mol，VVALKKVGVLKKPGPT复合物为?39.64 kcal/mol，二者在热力学上均属有利结合，数值相似。但能量分解显示机制截然不同：DMGAVA以适度的范德华力（?17.96 kcal/mol）和静电作用（?83.26 kcal/mol）为主，去溶剂化代价较小；VVALKKVGVLKKPGPT则依靠极强静电作用（?558.25 kcal/mol）和较大范德华贡献（?41.84 kcal/mol），但伴随巨大的极性溶剂化惩罚（568.20 kcal/mol），反映了多个赖氨酸残基去溶剂化的能量成本。这种焓-熵补偿现象解释了结构差异巨大的肽为何具有相近的结合亲和力。

讨论部分指出，尽管VVALKKVGVLKKPGPT在对接得分上更负，但PRODIGY计算的ΔG和K_d与DMGAVA相当，且其K_i更高，说明结合亲和力与抑制效率并不等同，活性位点内结合可能更有利于高效抑制。研究人员还强调，本研究的肽为胃蛋白酶水解产物，天然具有胃稳定性，但肠道稳定性、经肠吸收及跨膜能力均未验证；尤其是16肽VVALKKVGVLKKPGPT，其吸收受限，而六肽DMGAVA可能性较大。目前结果仅证明体外活性，尚不能外推至体内功效。

结论：本研究首次证明羽扇豆种子蛋白可释放精氨酸酶抑制肽，且不同结构肽可通过互补机制发挥抑制作用。白羽扇豆较蓝羽扇豆更具产生精氨酸酶抑制肽的潜力。动力学与计算模拟结合的分析表明，肽的尺寸、电荷分布和柔韧性共同决定其对精氨酸酶的调节方式。这些发现扩展了羽扇豆蛋白的功能应用价值，但转化为功能性食品或营养品仍需进一步研究肽在胃肠消化中的稳定性、完整形态吸收以及细胞和体内的有效性。未来研究应聚焦于胃肠道消化稳定性、生物可及性和体内生理相关性，以支持向功能性食品体系的转化。