Sonoda Milton T, Skaf Munir S
Institute of Chemistry, State University of Campinas-UNICAMP, Campinas, SP, Brazil.
J Phys Chem B. 2007 Oct 18;111(41):11948-56. doi: 10.1021/jp0749120. Epub 2007 Sep 22.
We use molecular dynamics simulations to investigate structure and dynamics of fructose aqueous solutions in the 1-5 M concentration range at ambient conditions. We analyze hydration structures, H-bond statistics, and size distribution of H-bonded carbohydrate clusters as functions of concentration. We find that the local tetrahedral order of water is reasonably well-preserved and that the solute tends to appear as scattered "isolated" molecules at low concentrations and as H-bonded clusters for less diluted solutions. The sugar cluster size distribution exhibits a sharp transition to a percolated cluster between 3.5 and 3.8 M. The percolated cluster forms an intertwined network of H-bonded saccharides that imprisons water. For the dynamics, we find good agreement between simulation and available experimental results for the self-diffusion coefficients. Water librational dynamics is little affected by sugar concentration, whereas reorientational relaxation is described by a concentration-independent bulk-like component attributed to noninterfacial water molecules and a slower component (strongly concentration dependent) that arises from interfacial solvent molecules and, hence, depends on the dynamics of the cluster structure itself. Analysis of H-bonding survival probability functions indicates that the formation of carbohydrate clusters upon increasing concentration enhances the H-bond relaxation time and slows down the entire system dynamics. We find that multiexponential or stretched-exponential fits alone cannot describe the H-bond survival probabilities for the entire postlibrational time span of our data (0.1-100 ps), as opposed to a combined stretched-plus-biexponential function, which provides excellent fits. Our results suggest that water dynamics in concentrated fructose solutions resembles in many ways that of protein hydration water.
我们使用分子动力学模拟来研究环境条件下浓度范围为1 - 5 M的果糖水溶液的结构和动力学。我们分析了水合结构、氢键统计以及氢键连接的碳水化合物簇的尺寸分布与浓度的函数关系。我们发现水的局部四面体有序性得到了较好的保留,并且溶质在低浓度时倾向于以分散的“孤立”分子形式出现,而在稀释度较低的溶液中则以氢键连接的簇形式出现。糖簇尺寸分布在3.5至3.8 M之间呈现出向渗流簇的急剧转变。渗流簇形成了一个由氢键连接的糖类交织网络,将水困住。对于动力学,我们发现模拟结果与自扩散系数的现有实验结果吻合良好。水的摆动动力学受糖浓度的影响较小,而重取向弛豫由归因于非界面水分子的与浓度无关的本体样成分以及由界面溶剂分子产生的较慢成分(强烈依赖浓度)来描述,因此取决于簇结构本身的动力学。氢键生存概率函数分析表明,随着浓度增加碳水化合物簇的形成会增加氢键弛豫时间并减缓整个系统的动力学。我们发现,与提供出色拟合的拉伸加双指数函数不同,单独的多指数或拉伸指数拟合无法描述我们数据整个摆动后时间跨度(0.1 - 100 ps)内的氢键生存概率。我们的结果表明浓果糖溶液中的水动力学在许多方面类似于蛋白质水合水的动力学。