Department of Chemistry and Biochemistry, 238 Brown Lab, University of Delaware, Newark, Delaware 19716, USA.
J Phys Chem B. 2010 Sep 2;114(34):11076-92. doi: 10.1021/jp101597r.
Building upon the nonadditive electrostatic force field for alcohols based on the CHARMM charge equilibration (CHEQ) formalism, we introduce atom-pair specific solute-solvent Lennard-Jones (LJ) parameters for alcohol-water interaction force fields targeting improved agreement with experimental hydration free energies of a series of small molecule linear alcohols as well as ab initio water-alcohol geometries and energetics. We consider short-chain, linear alcohols from methanol to butanol as they are canonical small-molecule organic model compounds to represent the hydroxyl chemical functionality for parametrizing biomolecular force fields for proteins. We discuss molecular dynamics simulations of dilute aqueous solutions of methanol and ethanol in TIP4P-FQ water, with particular discussion of solution densities, structure defined in radial distribution functions, electrostatic properties (dipole moment distributions), hydrogen bonding patterns of water, as well as a Kirkwood-Buff (KB) integral analysis. Calculation of the latter provides an assessment of how well classical force fields parametrized to at least semiquantitatively match experimental hydration free energies capture the microscopic structures of dilute alcohol solutions; the latter translate into macroscopic thermodynamic properties through the application of KB analysis. We find that the CHEQ alcohol force fields of this work semiquantitatively match experimental KB integrals for methanol and ethanol mole fractions of 0.1 and 0.2. The force field combination qualitatively captures the concentration dependence of the alcohol-alcohol and water-water KB integrals, but due to inadequacies in the representation of the microscopic structures in such systems (which cannot be parametrized in any systematic fashion), a priori quantitative description of alcohol-water KB integrals remains elusive.
基于 CHARMM 电荷平衡 (CHEQ) 形式的醇非加和静电势场,我们引入了针对醇-水相互作用力场的原子对特定溶剂化 Lennard-Jones (LJ) 参数,以提高与一系列小分子线性醇的实验水合自由能、从头算水-醇几何形状和能量的一致性。我们考虑甲醇到丁醇的短链、线性醇,因为它们是典型的小分子有机模型化合物,用于为蛋白质生物分子力场参数化代表羟基化学官能团。我们讨论了甲醇和乙醇在 TIP4P-FQ 水中稀水溶液的分子动力学模拟,特别讨论了溶液密度、径向分布函数定义的结构、静电特性(偶极矩分布)、水的氢键模式,以及 Kirkwood-Buff (KB) 积分分析。计算后者可评估至少半定量匹配实验水合自由能的经典力场参数化如何捕捉稀醇溶液的微观结构;通过 KB 分析的应用,将后者转化为宏观热力学性质。我们发现,本工作的 CHEQ 醇力场对半定量匹配甲醇和乙醇的实验 KB 积分,摩尔分数为 0.1 和 0.2。力场组合定性地捕捉到了醇-醇和水-水 KB 积分的浓度依赖性,但由于这些系统中微观结构的表示存在不足(无法以任何系统的方式参数化),醇-水 KB 积分的先验定量描述仍然难以捉摸。