基于肽组学引导的牛初乳肽三维结构-机制研究：多靶点抗糖尿病与抗氧化潜力的计算机模拟及体外验证

《Food Chemistry: X》：Peptidomics-guided 3D structural-mechanistic investigation of bovine colostrum peptides reveals multi-targeted in silico and in vitro antidiabetic and antioxidant potential

　　

摘要专业翻译 糖尿病（Diabetes mellitus）是全球主要的健康负担，其进展与氧化应激及慢性并发症相关，这对当前的药物疗法持续构成挑战。现有抗糖尿病药物往往仅靶向单一代谢通路，且可能伴随不良反应，因此亟需安全、多功能的替代方案。在此背景下，食源性生物活性肽因其低毒性、良好的生物利用度以及调节葡萄糖代谢中多个生物学靶点的能力而受到广泛关注。 牛初乳（Bovine colostrum）富含多种蛋白质组分，包括乳清蛋白、酪蛋白及乳脂球膜蛋白，经水解后可产生具有多功能特性的生物活性肽。尽管牛初乳在营养、抗菌和抗氧化方面的特性已获公认，但其肽库在抗糖尿病应用方面的潜力尚未被充分发掘。 抑制α-淀粉酶（α-amylase, AM）、α-葡萄糖苷酶（α-glucosidase, AG）和人醛糖还原酶（human aldose reductase, HAR）在碳水化合物代谢调控中发挥关键作用。AM和AG催化复杂碳水化合物分解为葡萄糖，HAR则通过多元醇通路将过量葡萄糖分流为山梨醇，该通路与糖尿病并发症密切相关。这些酶的协同活性会加速活性氧（reactive oxygen species, ROS）生成和氧化应激，进而导致与糖尿病相关的细胞损伤。这种相互依赖性表明，需要同时实现酶抑制和自由基清除，而非单独处理任一机制。 既往关于初乳或乳源肽降血糖效应的研究多聚焦于对AM和AG的单一抑制，而非多酶靶向策略。此外，能够同时靶向AM、AG和HAR的肽的实验证据仍然稀缺。因此，鉴定具有多酶抑制潜力的天然肽，可能为同时控制餐后高血糖并靶向与长期糖尿病并发症相关的通路提供综合策略。 本研究采用整合计算与实验工作流程，考察了先前从印度牛初乳水解物中鉴定的肽的抗糖尿病潜力。研究人员基于预测的生物活性和理化性质对候选肽进行初步筛选，随后结合分子建模与计算方法进行评估。分子对接用于鉴定肽与选定靶标催化残基之间的有利结合相互作用，随后进行分子动力学（molecular dynamics, MD）模拟以评估肽-靶标相互作用的持续性和所得复合物的稳定性。此外，MM/GBSA分析作为结合能的近似计算估计，用于支持肽-靶标相互作用的比较评估。先导肽进一步经计算机模拟评估其药代动力学性质、稳定性和胃肠道易感性，并进行三维建模，预测亲水/疏水平衡、空间相容性、分子柔性和静电势。最后，对合成的肽通过体外酶抑制和酶动力学实验以及抗氧化实验进行验证，从而将计算预测与实验证据相结合，阐明其潜在生物活性及肽-酶相互作用背后的分子特征。

论文解读文章

一、研究背景与问题提出

糖尿病是全球范围内重大的公共卫生负担，其病程进展与氧化应激和多种慢性并发症密切相关，对现有药物治疗方案构成持续挑战。目前临床使用的抗糖尿病药物大多仅靶向单一代谢通路，且可能伴随不良副作用，这使得开发安全、多功能的替代疗法成为迫切需求。在这一背景下，食源性生物活性肽因其低毒性、良好的生物利用度以及能够同时调节葡萄糖代谢中多个生物学靶点的独特优势，引起了研究人员的广泛关注。

牛初乳是多种蛋白质组分的丰富来源，包括乳清蛋白、酪蛋白以及乳脂球膜蛋白，这些蛋白在水解后能够产生具有多功能特性的生物活性肽。尽管牛初乳在营养、抗菌和抗氧化方面的价值已被充分认识，但其肽库在抗糖尿病应用方面的潜力尚未得到系统开发。

从酶学机制来看，α-淀粉酶（AM）和α-葡萄糖苷酶（AG）催化复杂碳水化合物分解为葡萄糖，是餐后血糖升高的关键调控节点；人醛糖还原酶（HAR）则通过多元醇通路将过量葡萄糖转化为山梨醇，该通路与糖尿病微血管并发症的发生发展密切相关。三种酶的协同作用会加速活性氧（ROS）生成和氧化应激，造成与糖尿病相关的细胞损伤。这种功能上的相互依赖性表明，理想的抗糖尿病策略需要同时实现多酶抑制和自由基清除，而非仅针对单一机制。

然而，既往关于初乳或乳源肽降血糖效应的研究大多集中于对AM或AG的单一抑制，缺乏多酶靶向的系统性策略。同时，能够同时靶向AM、AG和HAR的天然肽的实验证据仍然十分有限。因此，鉴定具有多酶抑制潜力的天然肽，有望为同时控制餐后高血糖并干预长期糖尿病并发症相关通路提供综合性的治疗策略。

二、研究内容与主要技术方法

本研究采用整合计算与实验的双轨工作流程，系统考察了先前从印度牛初乳水解物中鉴定的肽段的抗糖尿病潜力。研究人员首先基于预测的生物活性和理化性质对候选肽进行虚拟筛选，随后综合运用多种关键技术手段开展深入研究。主要技术方法包括：利用分子对接技术鉴定肽与靶酶催化残基之间的结合模式；采用分子动力学（MD）模拟评估肽-靶标复合物的动态稳定性和相互作用持续性；运用MM/GBSA分析对结合自由能进行近似计算估计；借助PEP-FOLD 4.0服务器进行肽段三维结构建模；通过ADMETLab 3.0、DeepSTABp和ExPASy PeptideCutter等工具预测药代动力学性质、稳定性和胃肠道耐受性；最后通过固相肽合成获得候选肽，并利用体外酶抑制、酶动力学及自由基清除实验进行实验验证。

三、研究结果

3.1 肽组学文库构建与先导肽筛选

研究人员从先前鉴定的牛初乳肽数据集出发，初始数据集包括乳清蛋白组（WP）167条肽、初乳脂球蛋白膜蛋白组（CFGMP）64条肽和酪蛋白组（CP）82条肽。经过去除各组分内重复序列，并排除分子量大于3 kDa、长度超出5-20个氨基酸范围、具有致敏性、毒性、高度疏水或易于聚集的肽段后，WP保留69条、CFGMP保留25条、CP保留34条。进一步经理化性质分析和PeptideRanker活性预测（评分>0.5）筛选后，最终WP保留19条、CFGMP保留13条、CP保留14条肽进入分子对接评估。通过系统的虚拟筛选流程，最终确定三条先导肽：P1（HMWPGDIK）、P2（GPAGPQGPR）和P3（DKVLVFFK）用于后续深入研究。

3.2 分子对接与分子动力学模拟分析

研究人员利用Glide软件的超高精度（XP）模式，将三条先导肽分别对接到AM（PDB ID: 1DHK）、AG（PDB ID: 3A4A）和HAR（PDB ID: 1IEI）的催化活性位点。分子对接结果显示，三条肽均能与靶酶形成稳定的结合构象。随后进行的100 ns分子动力学模拟进一步证实，肽-酶复合物在整个模拟过程中保持结构稳定，RMSD值维持在合理范围内。MM/GBSA结合自由能计算显示，三条肽与三种靶酶的结合自由能均为负值，表明结合过程在能量上是有利的，其中部分肽-酶组合的结合能甚至优于阳性对照药物阿卡波糖（Acarbose）和泽那司他（Zenarestat）。

3.3 肽段三维建模与理化性质分析

通过PEP-FOLD 4.0服务器对三条先导肽进行三维结构预测，并结合CABS-flex 3.0服务器分析肽段残基的柔性和接触图谱。三维结构分析揭示了氨基酸组成对肽段亲水/疏水平衡、空间相容性、分子柔性和静电势的影响，这些分子特征共同决定了肽段与靶酶之间的相互作用模式。研究发现，肽段中特定氨基酸残基的空间排布和理化特性对其与酶催化口袋的匹配程度具有决定性作用。

3.4 药代动力学与稳定性预测

ADMETLab 3.0预测结果显示，三条先导肽均展现出良好的药代动力学特征。DeepSTABp稳定性预测和ExPASy PeptideCutter胃肠道酶切位点分析表明，这些肽段具有合理的稳定性特征和可预测的胃肠道代谢行为，为其口服给药潜力提供了初步支持。

3.5 体外多酶抑制活性与抑制动力学

体外酶抑制实验显示，三条肽对AM、AG和HAR均表现出剂量依赖性的抑制效应。其中，P3（DKVLVFFK）展现出最为突出的多酶抑制潜力，其对三种靶酶的半数抑制浓度（IC₅₀）值均处于较低水平。酶动力学分析通过Lineweaver-Burk双倒数作图揭示了抑制机制，结果表明这些肽段主要通过竞争性或混合型抑制模式发挥作用，并计算获得了相应的抑制常数（K_i）值。

3.6 抗氧化活性评估

DPPH、超氧阴离子和ABTS自由基清除实验结果显示，三条先导肽均表现出浓度依赖性的自由基清除能力，其中部分肽段的半数有效浓度（EC₅₀）与阳性对照丁基羟基茴香醚（BHA）相当或更优。这一结果证实了这些肽段兼具抗糖尿病和抗氧化的双重功能特性。

四、讨论与结论

本研究的讨论部分强调，这是首次系统报道牛初乳来源肽段能够同时靶向AM、AG和HAR三种关键糖尿病相关酶的研究。研究通过整合计算机模拟预测与体外实验验证，建立了从肽组学数据到功能验证的完整研究范式。三维结构-机制分析揭示了氨基酸组成与肽段生物活性之间的内在关联，为合理设计具有多靶点活性的食源性肽类药物提供了理论依据。

研究结论指出，从牛初乳水解物中鉴定的三条生物活性肽（P1、P2和P3）展现出同时抑制AM、AG和HAR的潜力，并兼具显著的抗氧化活性，具有双重功能特性。这些发现为开发基于牛初乳肽的新型多功能抗糖尿病候选药物奠定了实验基础，同时也为利用食物来源生物活性肽防治糖尿病及其并发症提供了新的研究思路和科学依据。