《Foods》：Preparation, Binding Behavior and Molecular Simulation of Binary Complexes of Phloridzin with Whey Protein Isolate
编辑推荐：
乳清分离蛋白（WPI）可通过非共价相互作用与类黄酮组装成超分子复合物，但其潜在的结合机制尚不完全清楚。在本研究中，研究人员采用实验与计算相结合的方法系统地研究了WPI–根皮苷（PHL）复合物的形成机制。高效液相色谱（HPLC）定量分析了PHL的结合含量为1.3
乳清分离蛋白（WPI）可通过非共价相互作用与类黄酮组装成超分子复合物，但其潜在的结合机制尚不完全清楚。在本研究中，研究人员采用实验与计算相结合的方法系统地研究了WPI–根皮苷（PHL）复合物的形成机制。高效液相色谱（HPLC）定量分析了PHL的结合含量为1.3%（w/w）。等温滴定量热法（ITC）表明该过程是熵驱动的，且主要由疏水和静电相互作用主导。互补的圆二色光谱（CD）和分子动力学模拟（MD）揭示了复合作用诱导蛋白质二级结构发生适度的构象调整。综合而言，这种多尺度分析为WPI–PHL复合物的动态形成提供了机制性见解，并为蛋白质–类黄酮识别提供了理论启示。
论文解读
研究背景与意义：根皮苷（Phloridzin, PHL）是一种主要存在于蔷薇科和杜鹃花科植物中的多酚化合物，在苹果属植物中含量丰富。人类长期通过饮食摄入PHL，但多酚在胃肠道中的稳定性和靶向递送能力是其生物活性的关键。多酚在酸性条件下较稳定，而在小肠近中性pH环境中易降解，并可能与消化酶结合或被氧化降解。乳清分离蛋白（Whey Protein Isolate, WPI）具有氨基、羧基和羟基等基团，可通过非共价键与多酚结合，从而改变蛋白质结构并影响类黄酮的可及性和生物利用度。已有研究表明蛋白质可作为多酚的保护剂和载体，防止其在胃肠道消化过程中被酶解和氧化降解。然而，目前对类黄酮与WPI二元复合物的研究主要集中在制备和表征上，而对原子水平上分子间相互作用机制的系统研究仍有限，尤其是对于具有独特二氢查尔酮结构的PHL与WPI的分子相互作用机制，结合分子动力学（Molecular Dynamics, MD）模拟和可视化技术进行系统阐明的研究相对较少。因此，本研究采用实验与计算相结合的综合策略，阐明WPI与PHL之间的微观结合行为，旨在为蛋白质-类黄酮识别提供理论基础，并拓展PHL在功能食品工业中的应用。
主要关键技术与方法：研究人员采用高效液相色谱（HPLC）定量分析WPI–PHL复合物中PHL的结合含量；利用等温滴定量热法（ITC）测定热力学参数，包括结合常数、焓变和熵变；通过圆二色光谱（CD）分析蛋白质二级结构变化；采用分子对接（Autodock 4.2）和分子动力学模拟（AMBER99SB-ILDN力场，100 ns）计算结合自由能（MM-PBSA）和残基波动；并运用独立梯度模型（IGM）可视化非共价相互作用类型及空间分布。主要样本包括商业来源的WPI（蛋白质含量93.77%）、α-乳白蛋白（α-La, 纯度≥85%）和β-乳球蛋白（β-Lg, 纯度≥90%），PHL纯度98%。ITC实验在所有样品上以25°C进行，MD模拟采用1:1结合模型。
研究结果：
3.1 高效液相色谱分析：通过外标法建立PHL浓度与峰面积的标准曲线（y = 17,398x + 25.078），测定复合物样品中PHL质量百分比为1.3 ± 0.12%。该数值反映了经充分透析去除未结合PHL后与WPI稳定结合的PHL量，表明透析步骤有效去除未结合PHL。
3.2 圆二色光谱分析：CD光谱显示WPI和WPI–PHL复合物在210–220 nm处有负宽峰，表明富含β-折叠。计算二级结构含量发现，结合PHL后α-螺旋和β-折叠比例略有下降（如反平行β-折叠从27.2%降至25.7%），β-转角和无规卷曲略微增加（无规卷曲从33.3%增至34.56%），表明蛋白质发生适度的构象调整，有利于PHL进入疏水口袋，但整体二级结构骨架基本保持不变。
3.3 等温滴定量热法分析：分别研究α-La和β-Lg与PHL的结合。ITC结果显示结合过程为放热过程（ΔH < 0），熵变为正（β-Lg: +58.6 J·mol?1·K?1；α-La: +54.6 J·mol?1·K?1），表明范德华力、氢键以及疏水相互作用共同驱动结合。通过ΔG = –RT ln K
a计算得到β-Lg/PHL和α-La/PHL的ΔG分别为–22.6 kJ·mol?1和–23.4 kJ·mol?1，均为负值，表明结合自发进行。结合化学计量数N分别为1.7±0.2和1.5±0.4，提示每个蛋白分子平均结合多于一个PHL分子，但由于数据质量不支持复杂模型拟合，后续MD采用1:1结合模型聚焦于高亲和力位点。
3.4 分子动力学模拟结果：对分子对接获得的初始构象进行100 ns MD模拟，系统在约60 ns后达到动态平衡。均方根波动（RMSF）分析显示，结合PHL后两种蛋白质的大多数氨基酸残基更稳定；α-La中Glu49、Tyr103、Trp104和Lys108（位于疏水结合口袋）波动较大，β-Lg中Ile84、Phe105和Val118（位于水运输口袋）波动较大，表明这些残基积极参与与PHL的分子相互作用。溶剂可及表面积（SASA）分析表明，PHL结合后两种蛋白的SASA均增加，β-Lg的增加幅度（3.95 nm2）大于α-La（1.04 nm2），说明β-Lg在结合过程中发生更大程度的结构松弛。结合自由能计算（MM-PBSA）显示，PHL与α-La的结合自由能（-485.0 ± 0.8 kJ/mol）比与β-Lg（-77.5 ± 0.5 kJ/mol）更负，表明PHL对α-La具有更强的亲和力，与ITC趋势一致。
3.5 IGM分子间相互作用分析：IGM散点图显示在sign(λ?)ρ=0附近有密集的δg
inter散点，表明复合物中存在氢键和范德华力。三维等值面图显示，PHL分子的酚羟基、糖苷羟基和芳香环基团是关键的原子基团，与附近氨基酸残基形成相互作用，其中疏水氨基酸原子和PHL酚基团贡献较大。IGM方法清晰可视化了蛋白-PHL复合系统中参与相互作用的关键原子位置和强度。
讨论与结论：讨论部分指出，研究结果证实了PHL主要通过非共价相互作用与WPI形成稳定超分子复合物，这与氢键和范德华力主导蛋白-多酚结合的普遍认识一致。糖苷基团（C-5位）增强了与WPI的相互作用，与近期谱学研究发现糖基化类黄酮具有更高结合常数相符。MD模拟中SASA增加支持了“蛋白质诱导契合”模型，CD数据中有序结构减少暗示了蛋白构象调整以容纳PHL。ITC的负ΔG、放热ΔH和有利ΔS表明结合过程兼具焓和熵驱动特征，α-La更强的结合亲和力可能归因于其较高比例的正电荷残基。研究局限性包括：所有实验在可控实验室条件（pH 6.8, 25°C）下进行，未模拟胃肠道复杂环境；商业WPI为多种蛋白混合物，本研究仅单独研究α-La和β-Lg，未探索协同或竞争作用；MD采用1:1简化模型可能忽略协同结合效应；缺乏NMR或SAXS等进一步光谱验证。未来方向包括：模拟消化和生物利用度研究、AI驱动的蛋白质工程、多组分系统优化、冷冻电镜或固态NMR高级结构表征、安全性及感官评价。
研究结论部分翻译：结果表明，在根皮苷（PHL）与WPI的结合过程中，蛋白质的有序结构略有减少，伴随着无规卷曲含量适度增加（从33.28%增至34.56%），这可能是由于相互作用过程中蛋白质结构灵活性增强，有利于小分子进入疏水口袋并形成复合物。此外，PHL与WPI中的α-La和β-Lg的结合均为能自发进行的放热过程，PHL的羟基与氨基酸残基形成氢键，芳香环贡献疏水相互作用。由ITC和分子模拟获得的结合自由能和结合常数表明，PHL对α-La表现出更大（更负）的结合自由能和更高的结合常数，意味着对α-La具有更强的结合亲和力。这些发现表明WPI–PHL复合物具有作为功能性食品成分的潜力，其中WPI可作为有效载体提高PHL的稳定性和生物利用度。未来研究应评估该复合物在模拟胃肠道条件下的保护效果，并评估其在真实食品基质中的性能，同时通过定点突变确定次级结合位点的精确位置。
