最小距离分布函数（MDDF）解析水苏糖水合结构：致密首层、渗流氢键网络与动态水合核心

《Molecules》：Atomistic Insights into the Hydration Structure of Stachyose: A Minimum-Distance Distribution Function Approach

　　

水苏糖（stachyose）是豆科植物中丰富的棉子糖家族四糖（raffinose-family tetrasaccharide），具有显著的益生（prebiotic）和抗营养（anti-nutritional）活性，这些活性与其水相溶剂化行为密切相关。然而，水分子如何围绕这一非球形、柔性四糖组织，其原子尺度细节仍知之甚少，主要原因是传统径向分布函数（radial distribution function, RDF）对缺乏明确几何中心的分子并不适用。在此，研究人员采用分子动力学模拟（molecular dynamics simulations, MD；3个独立300 ns副本，GLYCAM06j/TIP3P）结合最小距离分布函数（minimum-distance distribution function, MDDF）框架，在多个长度和时间尺度上剖析水苏糖的水合结构。结果揭示了一个密集占据的第一水合层（3.5 ?），包含69.5 ± 0.3个水分子，并在1.86 ?处出现显著的MDDF峰，具有O–H···O氢键特征。直接氢键分析识别出每帧平均21.1 ± 0.1个溶质–水氢键，而驻留时间计算显示11.1 ± 0.1个水分子在5.5 ?壳层内停留超过200 ps，形成动态水合核心。值得注意的是，对第一壳层水进行DBSCAN聚类分析揭示了一个渗流氢键网络：最大瞬时簇涵盖约70个第一壳层水中的53–58个，但该网络在组成上具有流动性，单个水分子以约40 ps时间尺度交换，且任一时刻仅约16%的壳层水处于孤立状态。构象聚类进一步解析出四个主要家族（占全部帧的89.2%），末端残基（2CU和0LA）表现出比内部残基更高的溶剂暴露。令人意外的是，果糖末端（0LA）显示出与半乳糖末端（28.6%）基本相等的动态水合贡献（29.0%），尽管其亲水分数更小（65.9% vs. 半乳糖末端的70.0%），这凸显了仅基于几何的描述符的局限性。这些发现提供了水苏糖溶剂化的定量、残基分辨图景，并建立了一个适用于一般复杂碳水化合物的计算工作流程。

水苏糖（stachyose）是豆科植物中丰富的棉子糖家族四糖（raffinose-family tetrasaccharide）。棉子糖家族寡糖（raffinose-family oligosaccharides, RFOs）包括水苏糖、棉子糖和毛蕊花糖等，是重要益生元（prebiotic），可选择性促进有益肠道菌群生长；但人类小肠缺乏α-半乳糖苷酶，其在食品基质中的存在也与胀气和矿物质生物利用度降低有关。RFOs的生物学命运不仅由共价结构决定，也受水相溶剂化架构（aqueous solvation architecture）影响，因为水合层塑造构象集合、糖苷键对水解酶的暴露以及羟基对氢键伙伴的可及性。水合壳（hydration shell）也是寡糖与生物环境相互作用的主要介质，碳水化合物–蛋白识别常由水介导或调节，结合水可占结合界面相当比例并可通过桥接相互作用赋予特异性；寡糖呈现给受体的构象集合也受溶剂调控。糖类对不稳定生物分子的保护作用同样与水合有关，海藻糖（trehalose）的高水合数和扰动周围水氢键网络的能力是其生物保护 efficacy 的部分基础。因此，定量、原子尺度描述水如何围绕特定糖组织，是理解碳水化合物在食品科学、肠道生理和植物逆境生物学中功能的前提。实验上，碳水化合物水合已通过¹⁷O NMR弛豫、NMR与分子模拟结合、太赫兹吸收光谱以及体积或水活度测定水合数等方法研究；这些方法一致认为每个溶解糖分子动态扰动数十至数百个水分子，但空间分辨率有限，不能将水合特征归因于单个残基或官能团。MD模拟可填补此空白，但将原子尺度分析扩展到水苏糖这类较大柔性寡糖存在方法学困难：传统RDF需要选择几何中心（如质心），而该中心不对应化学意义上的“表面”；对非球形、柔性分子，中心参照RDF首峰可能反映分子平均形状而非真实溶剂–溶质接触距离，对线性或支化寡糖尤为严重，因为不同残基向溶剂突出的程度差异很大。最小距离分布函数（minimum-distance distribution function, MDDF）通过计算每帧任意溶质原子与任意溶剂原子之间的最小距离分布，规避该问题；所得分布本质形状无关，直接联系Kirkwood–Buff溶剂化热力学（Kirkwood–Buff solvation thermodynamics），并可分解为原子或残基水平贡献，提供自然的“溶剂壳”视角。研究人员因此结合900 ns显式溶剂MD模拟（3个独立300 ns副本）、MDDF分析、氢键谱、驻留时间测量和构象聚类，构建水苏糖水合的多层、原子尺度图景，并具体回答：第一水合层的厚度与数量；四个单糖残基溶剂暴露差异；水合水的寿命与连通模式；分子柔性如何调节SASA提供的静态几何图景。研究得出：水苏糖具有薄而致密的第一水合层，氢键直接主导；首层水形成拓扑渗流网络但组成快速交换，无长寿命有序簇；存在约11个长驻留水的动态水合核心；末端残基暴露高于内部残基，但静态SASA不能完全解释动态水合排序，果糖末端的呋喃环柔性带来超出亲水分数预期的动态水合贡献；总几何水合数对骨架波动较稳健。该工作提供了残基分辨的定量溶剂化图景，并建立了可推广至复杂碳水化合物的计算流程，对理解RFOs的溶解度、酶可及性和界面识别具有意义。

关键技术方法（≤250字）：采用3个独立300 ns、总计900 ns的显式溶剂MD（GLYCAM06j/TIP3P）；用ComplexMixtures.jl计算MDDF与Kirkwood–Buff积分并按残基/原子类型分解；用cpptraj按供体–受体距离3.0 ?、角度135°统计氢键；用MDAnalysis脚本在5.5 ?壳层计算驻留时间（10 ps保存，≥200 ps为长驻留，桥接≤2帧出壳）；用DBSCAN对3.5 ?首层水氧位置聚类（ε=3.0 ?，min_samples=2）；用FreeSASA/cpptraj计算SASA，并以平均链接RMSD进行构象聚类。无实验样本队列；体系为单个水苏糖分子溶于约2579个TIP3P水分子，未加抗衡离子。

研究结果如下。3.1 总体水合结构：MDDF与配位数（Overall Hydration Structure: MDDF and Coordination Numbers）。通过MDDF和watershell计数，发现水苏糖MDDF在1.86 ?出现尖锐首峰、幅度2.24，首层约比体相水密2.3倍，首峰位置与O–H···O氢键O···O间距一致；第一极小约3.51 ?，被取为第一水合层边界。3.5 ?内水数为69.5 ± 0.3，5.5 ?内为158.6 ± 0.7；首层约相当于每个羟基五个接触水或每个极性氧约三个水。传统质心参照RDF在相同距离范围仅弱而宽，无明确首峰，说明其对伸长四糖会抹平真实接触距离。MDDF在约5.5 ?外衰减至1且无明显结构，KB积分在8 ?外达平台，替换bulk范围12–16 ?或15–20 ?只改变排除体积渐近值而不改变短程特征；8.0 ?处KB积分为?412 ?³，原文解释为大体中性溶质的排除体积主导，物理有意义的水合数应取短程MDDF积分或直接几何计数，而非渐近KB值。3.2 方向性氢键与水桥（Directional Hydrogen Bonding and Water Bridges）。通过直接氢键统计，发现每帧平均21.1 ± 0.1个溶质–水氢键，其中水作供体11.8 ± 0.1、水苏糖羟基向水供体9.3 ± 0.1；水→溶质供给略多源于四个环氧和三个糖苷氧只能作受体。平均每个自由羟基约1.5个氢键，落在单糖和二糖报道范围内；绝对数较高反映四糖极性位点多（21个氧受体和14个羟基供体），GLYCAM06j/TIP3P基本以残基加和方式再现局部氢键统计。通过DBSCAN首层水聚类，发现按所有帧簇出现计，孤立水占51.4%、二聚桥19.0%、三聚9.3%、四水及以上20.3%；单帧内平均仅10.9 ± 0.1个首层水孤立，4.02 ± 0.04个二聚桥跨邻羟基，其余约45个水组成少数较大聚集体；最大瞬时簇为53–58个水，最高约占首层80%，表明首层拓扑渗流为单一氢键网络而非残基局部碎片。但驻留时间显示个体水约40 ps交换，网络组成流动；按水分子计任一时刻仅约16%首层水孤立。原文强调这种连通不等于结构有序，并与蛋白疏水腔中受限域的有序长寿命水簇形成对比；水苏糖小、均匀亲水、凸形、首层薄（约2 ?）且无疏水斑块，水–水缔合几何约束少，故渗流几乎始终存在，同时网络被水的大幅跳跃机制持续重组。3.3 动态水合核心：驻留时间（Dynamic Hydration Core: Residence Times）。通过5.5 ?壳层驻留时间计算，发现几何计数虽为158.6个水，但仅11.1 ± 0.1个水驻留超过200 ps，构成“动态水合核心”，约占几何壳7%、首层16%；驻留时间分布强烈偏斜，平均段约42 ps、中位30 ps，仅1.0%段超过200 ps，超过500 ps极罕见。该时间尺度处于生物分子界面水从暴露凸位点几皮秒到凹形多位点氢键数百皮秒的范围内；水苏糖无深腔，处于快交换端。原文指出TIP3P已知低估水黏度和驻留时间，故驻留数值应视为下界；多数5.5 ?壳层水为近体相交换的瞬态访客，约11个长驻留水构成更有序层，可能主导太赫兹和NMR弛豫等谱学可观测。3.4 残基水平溶剂化不对称（Residue-Level Solvation Asymmetry）。通过MDDF残基分解，发现末端半乳糖2CU_1和果糖0LA_4合计约占首层MDDF积分58%，虽仅占溶质原子50%；内部半乳糖6GA_2和葡萄糖6LA_3分别贡献19.2%和23.2%。这源于(1→6)连接寡糖的伸展“天线状”拓扑：内部残基受邻居空间遮蔽，末端残基自由伸入溶剂；按原子分解显示高暴露位点主要是羟基氧及其氢，环氧和C–H表面贡献小。通过SASA计算，发现四残基静态暴露意外均匀，每残基SASA仅162.2–172.3 ?²，末端半乳糖2CU_1略大（172.3 ?²）；所有残基均以亲水O/OH表面为主，疏水C/CH贡献均低于40%；亲水分数排序为2CU_1（70.0%）>6GA_2（67.2%）>0LA_4（65.9%）>6LA_3（62.5%）。但静态描述符不能预测动态水合排序（0LA_4≈2CU_1>6LA_3>6GA_2）：果糖末端0LA_4动态水合贡献29.0%与半乳糖末端28.6%基本相等，尽管亲水分数更低；6LA_3亲水分数最低但动态水合23.2%仍高于6GA_2的19.2%。原文将0LA_4异常高的动态水合归因于β-呋喃果糖环的构象可塑性：呋喃环可在低能景观上于信封和扭曲构象间自由假旋转，褶皱波动瞬时暴露平均结构中几何埋藏的羟基，从而增强动态溶剂可及性；并据此指出静态SASA等几何描述符不足以排序碳水化合物水合倾向，存在柔性环或连接时必须考虑动态采样。3.5 构象采样与水合稳健性（Conformational Sampling and Hydration Robustness）。通过900 ns联合集合的构象聚类，发现四个主家族共占89.2%帧：Cluster 0（49.4%）为伸展构象（重原子回转半径约5.1 ?，端到端距离约10.7 ?），两端残基伸入溶剂；Cluster 1（17.8%）和Cluster 2（11.8%）为弯曲态（Rg约4.5–4.8 ?，端到端约7 ?），源于中央α-1,6连接折叠；Cluster 3（10.2%）最紧凑（Rg约4.3 ?，两个骨架弯曲角>100°）且簇内定义最紧（簇内平均距离1.22 ?）；其余六小簇各<5%，为罕见瞬态。伸展构象占优与(1→6)连接额外ω扭转带来的更大构象自由度和水溶液中偏好的伸展几何一致。总几何水合数在各副本围绕约159个水波动，与全局构象转变无显著相关；瞬时壳层占据降至约130个水在各簇和副本中均出现，反映局部交换而非构象门控。原文讨论常规MD不能穷尽柔性四糖糖苷构象空间，但每个副本访问所有四个主家族，两个(1→6)连接的ω扭转每副本互变76–165次；与MM3 1→6连接立体能图比较，采样(φ, ψ, ω)区域对应所有三个ω旋转异构体的交错极小，但真空和晶体中偏好的trans-ψ区域在水溶液中未布居；因此采样足以支持所研究的局部水合可观测，但不能捕捉多连接慢重排。GLYCAM06j/TIP3P经蔗糖对照验证：晶体构象（108°, 305°）位于全局自由能最低内，参考GLYCAM06景观的四个特征basin在5°内重现；GLYCAM06对第二蔗糖basin M2的过稳定（ΔF(M2?M1)≈0.01 kcal mol^?1）不影响本工作结论，因为所报水合可观测为局部最小距离量，对次级糖苷basin相对布居不敏感。葡吡喃糖环全程保持⁴C₁（Cremer–Pople θ=9.8 ± 5.4°，99.8%帧低于30°），果呋喃糖假旋转相位围绕晶体值波动（圆均值?11.2° vs. 晶体?5°），与既往GLYCAM06蔗糖模拟一致。

讨论与结论总结：原文讨论强调，水苏糖首层水虽拓扑渗流为单一氢键网络，但快速组成周转和缺乏长寿命有序簇表明其表面均匀亲水、无疏水口袋，与蛋白限域水合位点不同；网络连通不能误作结构有序。TIP3P低估水黏度和驻留时间，故驻留数据为下界，但“小持久核心嵌入快交换群体”的定性结论稳健。常规MD不能穷尽柔性四糖构象空间，且力场对溶质–溶质相互作用和部分电荷处理存在差异；本研究所用GLYCAM06j/TIP3P经溶液相碳水化合物验证及蔗糖对照支持，且所报局部最小距离水合量对全局构象状态不敏感，因此两个 caveat 不影响核心结论：水苏糖水合主要由残基身份和末端暴露决定，而非全局骨架构象。研究结论译述如下：本研究首次利用MDDF框架给出具有生理意义四糖水苏糖的原子尺度水合图；水苏糖具有致密而薄的第一水合层（3.5 ?内69.5 ± 0.3个水），由直接O–H···O氢键主导（每帧21.1 ± 0.1个氢键），与单糖和二糖已建立氢键化学计量一致；首层水虽渗流为单一氢键网络，但快速组成周转（平均驻留约42 ps）和缺乏长寿命有序簇表明均匀亲水表面无疏水口袋，与蛋白限域水合位点相反；更宽几何壳内存在11.1 ± 0.1个长驻留（>200 ps）水的动态水合核心，这些水可能介导太赫兹和NMR弛豫实验观察到的谱学和热力学信号，其余壳层水为瞬态、类体相旁观者；末端残基暴露高于内部残基，但静态SASA不能捕捉全貌：果糖末端0LA动态水合贡献与半乳糖末端相当（29.0% vs. 28.6%），尽管亲水分数更低（65.9% vs. 70.0%），源于呋喃环柔性诱导的瞬时暴露；总几何水合数对适度骨架波动稳健，提示糖苷连接旋转异构体的环境扰动只会弱改变水苏糖溶解度和酶可及性；此处展示的计算工作流程——结合MDDF分解、方向性氢键分析、驻留时间测量和构象聚类——可方便转移至其他复杂碳水化合物、多糖和糖缀合物，适用于传统RDF方法不足的情形。展望上，副本交换等增强采样以及跨棉子糖家族（棉子糖、水苏糖、毛蕊花糖）和脱氧类似物的系统比较，可将局部水合图景连接至RFOs的全局构象热力学。