解析pH对牛α-乳白蛋白分子构象与界面性质关系的影响：一种分子动力学方法

《Food Chemistry》：Deciphering pH effects on relationships between molecular conformation and interfacial properties of bovine α-lactalbumin: a molecular dynamics approach

　　

原子尺度分子动力学（Molecular Dynamics, MD）揭示pH和钙结合状态如何塑造α-乳白蛋白（α-lactalbumin, α-LA）的构象。Asp82质子化先于Holo-LA在pH 3.0时Ca²⁺释放。酸性pH增加局部去折叠和两亲性暴露。Apo-LA在pH 3.0时泡沫衰减比Holo-LA更慢。尽管界面张力值相似，但起泡性能存在差异。

**论文解读：pH调控牛α-乳白蛋白构象与界面性质的分子机制**

**1. 研究背景与现存问题**

乳清蛋白作为食品工业中重要的界面活性成分，广泛用于稳定泡沫和乳液体系。α-乳白蛋白（α-Lactalbumin, α-LA）是牛乳清中第二大丰度蛋白，分子量14.1 kDa，等电点（pHi）4.8，其天然功能是参与乳糖合成。既往研究表明，α-LA能够在油-水或空气-水界面形成薄膜，降低界面张力，因而具有乳化性和起泡性。然而，这些技术功能性高度依赖介质条件（如pH、温度、离子存在）。实验研究虽已观察到pH变化显著改变α-LA的起泡和乳化能力，并初步归因于二级结构变化（如酸性pH下α-螺旋含量增加、无规卷曲减少），但所用荧光光谱、圆二色光谱和量热法只能给出整体结构变化的信号，无法精确定位结构改变发生的氨基酸残基，也无法识别具体分子内相互作用的形成或断裂。因此，pH如何通过调控α-LA的分子构象来影响其界面性质，其分子机制仍不清楚。分子动力学（Molecular Dynamics, MD）模拟能够提供原子级分辨率，揭示蛋白质在界面组装中的构象行为。然而，以往计算研究多聚焦温度或配体结合，尚未有研究系统比较无钙（Apo）与含钙（Holo）α-LA在多个pH条件下的构象差异并耦合残基特异性质子化状态。

**2. 研究目的与意义**

本研究旨在通过MD模拟结合体外实验，系统考察牛α-LA在其等电点（pH 4.8）及上下各1.8个单位（pH 3.0和6.6）条件下、有无天然Ca²⁺结合时的构象行为，并关联其起泡性质。该研究首次在同一比较矩阵中评估Apo-和Holo-α-LA，整合钙位点稳定性、整体与局部构象重排、溶剂可及表面积（SASA）和局部离子分布，并与实验测定的起泡行为相关联，从而建立从原子细节到技术功能的解释框架。论文发表在《Food Chemistry》。

**3. 主要技术方法**

研究人员采用GROMACS 2018软件包，基于PDB晶体结构（Holo-LA：1F6S；Apo-LA：1F6R）构建六个模拟体系。利用PlayMolecule服务器中的PROPKA预测残基特异性pKa，在pH 3.0、4.8、6.6下分配各可电离残基的质子化状态。使用GROMOS 53A6力场，SPC水模型，0.15 M NaCl离子强度，25°C。每个体系经过多阶段能量最小化和15 ns NPT平衡后，进行200 ns生产模拟。轨迹分析包括均方根偏差（RMSD）、均方根涨落（RMSF）、二级结构演变（DSSP）、SASA、径向分布函数（RDF）和离子密度分布。体外实验采用滴重法测定界面张力，通过注射器往复发泡测定起泡容量（FC）和泡沫稳定性（FS），并使用光学显微镜观察泡沫微观结构。所有实验在25°C、蛋白质浓度70.5 μM（约0.1% w/v）下进行。

**4. 研究结果**

**4.1 全局构象稳定性与二级结构演变**

在pH 4.8和6.6下，Apo-和Holo-LA的RMSD均较低（最大<3.5 ?），结构稳定。pH 3.0时，两种形式均出现显著增大的构象波动；Apo-LA出现两次大的RMSD跃迁（约4.5 ?和5.5 ?），Holo-LA在80 ns后逐渐转变至约120 ns稳定。残基水平分析显示，pH 3.0下H1、H2螺旋和S1、S2折叠片是主要波动区域。DSSP分析表明，Apo-LA在pH 3.0下失去310螺旋h1b和h3c，但h2b发生结构增益转变为稳定α-螺旋；H3螺旋延伸纳入额外残基（Lys98-Trp104）。Holo-LA在pH 3.0下不可逆失去h1b，S1-S2 β-折叠在125-160 ns部分去稳定后重新折叠；钙解离后de novo形成α-螺旋（Ile101-Lys103）。

**4.2 钙位点稳定性与Ca²⁺解离**

Ca²⁺仅在pH 3.0时从其初始配位几何中解离。轨迹分析表明，约20 ns后，质子化的Asp82不再与Ca²⁺保持静电作用，Ca²⁺-氧距离从约2.4 ?增至>20 ?。Asp82质子化引起与Lys79的静电排斥，触发配位位点去稳定并最终释放钙离子。pH 4.8和6.6下Ca²⁺始终稳定结合。

**4.3 溶剂可及性、水合与离子分布**

RDF分析显示，水分子分布在所有pH下相似。离子分布呈现明显的pH依赖：pH 3.0时Cl^?在蛋白质表面富集（第一峰g(r)≈20），pH 6.6时Na⁺富集（g(r)≈22），pH 4.8时两种离子分布较均衡。SASA分析表明，pH 3.0下亲水和疏水残基的SASA均增大，指示两亲性暴露增加；pH 4.8下疏水SASA下降约5%，蛋白质更紧凑；pH 6.6下疏水SASA增大尤其是在Apo形式中。

**4.4 起泡性与界面张力**

起泡容量不受钙结合状态显著影响，但受pH显著影响：pH 4.8时最低，pH 3.0时高于pH 6.6。拟合模型为FC% = 0.073·pH2 ? 0.3094·pH + 0.39。界面张力在所有条件下无显著差异（p>0.05）。泡沫稳定性方面，pH 3.0时Apo-LA产生的泡沫比Holo-LA衰减更慢；pH 4.8和6.6下无显著差异。显微镜图像显示pH 4.8时气泡数量少且尺寸小，pH 3.0时Apo-LA气泡更多更大。

**5. 讨论与结论**

讨论部分指出，pH 3.0下α-LA构象灵活性增加和两亲性残基暴露是其起泡性能改善的分子基础。pH 4.8时蛋白质更致密，需要更高表面浓度覆盖界面，因此起泡能力下降。pH 6.6时疏水暴露增加有助于中等起泡能力。界面张力无显著差异表明该参数单独不能解释起泡行为差异，构象重排和暴露的疏水区域是更关键的因素。

研究结论：本研究结合MD模拟与体外实验，阐明了pH对牛α-LA结构动力学和起泡性质的影响。pH 3.0下α-LA发生显著构象变化，增加亲水和疏水残基的SASA，这与实测的更高泡沫容量和稳定性相关。Holo-LA在pH 3.0下失去钙离子，伴随310螺旋h1b不可逆丢失和β-折叠部分去稳定，但新α-螺旋的形成增加灵活性。pH 4.8下SASA变化最小，起泡性降低。pH 6.6下疏水残基暴露增加，尤其Apo形式，改善泡沫稳定性。Apo-LA在pH 3.0时展现出更优的泡沫稳定性，归因于酸性环境诱导的更大构象灵活性。该研究为食品工业中选择α-LA作为界面配料时的pH和钙负载条件提供了分子层面的指导，例如用于酸奶充气产品、酸性起泡饮料和乳-植物蛋白混合基质。