Makowski Mariusz, Liwo Adam, Scheraga Harold A
Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA.
J Phys Chem B. 2007 Mar 22;111(11):2910-6. doi: 10.1021/jp065916s. Epub 2007 Feb 27.
A physics-based model is proposed to derive approximate analytical expressions for the cavity component of the free energy of hydrophobic association of spherical and spheroidal solutes in water. The model is based on the difference between the number and context of the water molecules in the hydration sphere of a hydrophobic dimer and of two isolated hydrophobic solutes. It is assumed that the water molecules touching the convex part of the molecular surface of the dimer and those in the hydration spheres of the monomers contribute equally to the free energy of solvation, and those touching the saddle part of the molecular surface of the dimer result in a more pronounced increase in free energy because of their more restricted mobility (entropy loss) and fewer favorable electrostatic interactions with other water molecules. The density of water in the hydration sphere around a single solute particle is approximated by the derivative of a Gaussian centered on the solute molecule with respect to its standard deviation. On the basis of this approximation, the number of water molecules in different parts of the hydration sphere of the dimer is expressed in terms of the first and the second mixed derivatives of the two Gaussians centered on the first and second solute molecules, respectively, with respect to the standard deviations of these Gaussians, and plausible analytical expressions for the cavity component of the hydrophobic-association energy of spherical and spheroidal solutes are introduced. As opposed to earlier hydration-shell models, our expressions reproduce the desolvation maxima in the potentials of mean force of pairs of nonpolar solutes in water, and their advantage over the models based on molecular-surface area is that they have continuous gradients in the coordinates of solute centers.
提出了一种基于物理的模型,以推导水中球形和椭球形溶质疏水缔合自由能的空腔分量的近似解析表达式。该模型基于疏水二聚体和两个孤立疏水溶质水合球中水分子数量和环境的差异。假设接触二聚体分子表面凸部的水分子以及单体水合球中的水分子对溶剂化自由能的贡献相同,而接触二聚体分子表面鞍部的水分子由于其迁移率更受限(熵损失)以及与其他水分子的有利静电相互作用较少,导致自由能有更显著的增加。单个溶质粒子周围水合球中的水密度通过以溶质分子为中心的高斯函数相对于其标准差的导数来近似。基于这种近似,二聚体水合球不同部分的水分子数量分别用两个以第一和第二溶质分子为中心的高斯函数相对于这些高斯函数标准差的一阶和二阶混合导数来表示,并引入了球形和椭球形溶质疏水缔合能空腔分量的合理解析表达式。与早期的水合壳模型不同,我们的表达式重现了水中非极性溶质对平均力势中的去溶剂化最大值,并且它们相对于基于分子表面积的模型的优势在于它们在溶质中心坐标处具有连续梯度。