Center for Drug Design and Development, The University of Toledo, Toledo, Ohio 43606-3390, USA.
J Phys Chem A. 2012 Jul 26;116(29):7726-41. doi: 10.1021/jp304164g. Epub 2012 Jul 17.
Theoretical calculations up to the ab initio IEF-PCM/CCSD(T)/CBS//IEF-PCM/B3LYP/6-311++G** and IEF-PCM/B97D/aug-cc-pvtz levels have been performed for 2X-ethanol and 2X-phenol systems with X = F, NH(2), NO(2) in chloroform and aqueous solution. The calculated relative free energies by means of the IEF-PCM continuum dielectric method do not differ very much at the DFT and ab initio levels. Application of explicit solvent models and the FEP/MC method for determining relative solvation free energies causes, however, large deviations in the predicted equilibrium compositions, although the predominant conformation for the solute is generally in agreement with that from the corresponding IEF-PCM calculations. Existence of an intramolecular hydrogen bond (HB structure) for species with the HO-C-C-X moiety is preferred compared with a conformation when the hydrogen bond is disrupted (NoHB) for the considered F- and NO(2)-substituted molecules both in chloroform and aqueous solution. For 2NH(2)-ethanol, the HB structure is predominant in chloroform, whereas the 93:7 ratio for the OCCN trans/gauche species was obtained in aqueous solution. 2NH(2)-phenol exhibits a subtle equilibrium of the HB and NoHB conformations in both solvents. Potential of mean force calculations predict about a 10% solute association for the trans 2NH(2)-ethanol solute even in the fairly dilute 0.22 molar solution, whereas direct MC simulations do not support the maintenance of a doubly hydrogen-bonded dimer. Aqueous solution characteristics, as coordination numbers and numbers of strongly bound water molecules to the solute at T = 298 K and p = 1 atm, correspond reasonably to the derived molecular structures.
已在氯仿和水溶液中对 2X-乙醇和 2X-苯酚体系(X = F、NH(2)、NO(2))进行了从头算IEF-PCM/CCSD(T)/CBS//IEF-PCM/B3LYP/6-311++G**和 IEF-PCM/B97D/aug-cc-pvtz 水平的理论计算。通过 IEF-PCM 连续介电方法计算的相对自由能在 DFT 和从头算水平上差异不大。然而,应用显式溶剂模型和 FEP/MC 方法确定相对溶剂化自由能会导致预测平衡组成出现较大偏差,尽管溶质的主要构象通常与相应的 IEF-PCM 计算一致。与氢键被破坏(无氢键,NoHB)的构象相比,HO-C-C-X 部分的物种存在分子内氢键(HB 结构),这在氯仿和水溶液中对于考虑的 F-和 NO(2)-取代分子都是有利的。对于 2NH(2)-乙醇,在氯仿中 HB 结构占主导地位,而在水溶液中则得到 OCCN 反式/ gauche 物种的 93:7 比例。2NH(2)-苯酚在两种溶剂中均表现出 HB 和 NoHB 构象的微妙平衡。平均力势计算预测,即使在相当稀的 0.22 摩尔溶液中,trans 2NH(2)-乙醇溶质的溶质缔合率约为 10%,而直接 MC 模拟不支持保持双氢键二聚体。水溶液特性,如在 T = 298 K 和 p = 1 atm 时的配位数和与溶质结合的强结合水分子数,与推导的分子结构相当吻合。