利用互补的AIMD与经典MD模拟解析果糖互变异构体的复杂氢键和水合模式

《ACS Physical Chemistry Au》：Unraveling the Complex Hydrogen-Bonding and Solvation Patterns of Fructose Tautomers via Complementary AIMD and Classical MD Simulations

　　

理解果糖互变异构体在分子层面的水合作用仍是一项挑战，这源于其快速的相互转化和结构复杂性。研究人员采用了一种双尺度方法，将长时间经典分子动力学（MD）和从头算分子动力学（AIMD）相结合，以统计学和电子层面的准确性阐明吡喃糖和呋喃糖互变异构体氢键网络的差异。模拟结果突出了由羟基供体和受体位点数量有限所驱动的分子内与分子间氢键之间的内在竞争。关键在于，这些内部键的形成受柔性呋喃糖中环折叠（ring puckering）以及两类互变异构体中羟甲基空间取向的驱动。这种结构上的相互影响间接作用于整个分子间溶剂化壳层，该壳层同时包含亲水区和疏水区。通过将这些相互竞争的分子内和分子间模式与Bader电荷及水偶极矩相关联，研究人员为碳水化合物溶剂化建立了定量的电子结构基础。这些见解提供了一个全面的原子级框架，加深了研究人员对生物相关过程（如代谢酶中的识别和对接效率）的理解。

果糖是生物体内重要的单糖，其在水溶液中的行为直接影响生物系统的识别与代谢过程。水分子介导果糖各结构形式之间的连续平衡，使糖能够适应不同的化学环境。果糖在水溶液中存在五种互变异构体，其中β-d-吡喃果糖占68.23%，β-d-呋喃果糖占22.35%，此外还有少量α-异构体和开链形式。吡喃环采用稳定的椅式构象，而呋喃环则经历持续的假旋转，具有较高的构象灵活性。果糖的五个羟基和环氧提供了丰富的氢键位点，但位点数量有限，导致分子内与分子间氢键存在竞争。尽管已有多种实验和计算方法研究果糖水溶液，但关于局部氢键环境的关键问题仍未解决，特别是羟甲基取向对氢键网络的影响、亲水与疏水区域的相互作用机制等。为此，研究人员采用经典分子动力学（MD）与从头算分子动力学（AIMD）相结合的双尺度计算策略，旨在从统计和电子结构两个层面全面解析果糖互变异构体的水合行为。

研究人员构建了四种果糖异构体（α/β-d-呋喃果糖和α/β-d-吡喃果糖）的水溶液模型，每个体系包含一个果糖分子和99个水分子。首先进行经典MD模拟，采用GROMACS软件包、CHARMM36力场和SPC/E水模型，在NVT系综下于298 K运行600 ns，获得充分的构象采样。随后从经典MD轨迹中提取代表性构象作为AIMD模拟的起点，使用CP2K软件包在DFT/BLYP-D3水平上计算力，采用GTH赝势和TZV2P基组，平面波截断能300 Ry，模拟至少150 ps。通过Bader电荷分析和Wannier偶极矩计算表征电子性质，并利用空间密度函数分析水分子三维分布。

**构象分析**：通过经典MD模拟计算二面角分布和环折叠参数。结果表明，呋喃果糖的环折叠呈现多样性，α-和β-异构体分别偏好不同的假旋转状态，且羟甲基取向与环折叠强烈耦合；吡喃果糖则分别采用⁵C₂和²C₅椅式构象，刚性较强。羟甲基的空间取向直接决定了分子内氢键的可能位点。

**分子内氢键性质**：基于经典MD模拟，采用几何判据（O–H···O距离小于2.5 ?，角度小于60°）分析分子内氢键概率。呋喃糖中，两个羟甲基（OH1和OH6）以及OH6与环氧（O5）的相互作用是主要分子内氢键来源；吡喃糖中，分子内氢键较少且概率较低，β-d-吡喃果糖在严格判据下甚至无法检测到分子内氢键。环折叠和羟甲基取向共同决定了分子内氢键的模式。

**分子间氢键性质**：通过经典MD和AIMD模拟计算亲水壳层配位数和各氧位点的氢键数。结果显示，β-d-吡喃果糖与水的氢键数显著多于α-异构体，且所有OH位点几乎均匀参与供体和受体作用；呋喃糖中β-异构体也表现出更好的受体能力。环氧原子在所有体系中最弱的受体，这与Bader电荷分析一致。经典MD与AIMD结果定性相符，定量差异源于电子极化效应和采样时间不同。

**果糖的电子性质**：基于AIMD模拟的Bader电荷分析显示，环氧原子携带最小的负电荷，故其氢键接受能力最弱。C2位电荷最小，归因于异头效应。α-异构体中OH2位负电荷略大于β-异构体，这与异头中心的立体电子构型相关。

**水的电子性质**：通过Wannier偶极矩分析水分子。四配位水分子的偶极矩大于三配位水分子。当果糖作为氢键供体时，与之直接成键的水分子偶极矩通常大于作为受体时，与糖羟基酸性强于水的实验事实一致。氢键类型（供体/受体数量）影响偶极矩，而第一配位层水与体相水之间无显著差异。

**水合层性质**：通过空间密度函数分析发现，亲水壳层的形状受羟甲基取向和环折叠影响，β-d-吡喃果糖的亲水壳层呈不对称分布，与OH4的赤道取向相关。疏水区域的水分子通过范德华力和CH···O弱相互作用聚集在碳氢骨架周围，但仍与周围水形成氢键网络。

讨论部分指出，经典MD与AIMD互补：经典MD提供纳秒尺度的统计采样，AIMD能准确描述电子极化和动态电荷转移。β-吡喃糖具有最高的氢键效率，这可能解释其更高的溶解度和宏观丰度。局部水合网络对生物识别和对接过程具有关键作用。结论部分总结：该项工作全面表征了果糖的复杂氢键模式，证明双尺度方法能在不牺牲统计显著性的前提下获得准确的电子结构描述；分子内与分子间氢键网络存在紧密耦合，特定位点参与分子内键后即无法与周围水作用；通过Bader电荷和水偶极矩关联，为碳水化合物溶剂化提供了定量的电子结构基础。这些位点特异性的水合模式有助于理解环折叠和羟甲基取向如何影响碳水化合物在代谢酶中的识别与对接效率。该研究发表于《ACS Physical Chemistry Au》。