Alagona Giuliano, Ghio Caterina
CNR-IPCF, Institute for Physical Chemistry Processes, Molecular Modeling Lab, Via Moruzzi 1, I-56124 Pisa, Italy.
J Phys Chem A. 2006 Jan 19;110(2):647-59. doi: 10.1021/jp053612k.
All possible combinations of stable dihedral values have been considered in vacuo at the B3LYP/6-31G level for 3,9-dihydroxy-4,8-diprenylpterocarpan (erybraedin C), whose hydroxy out-out conformation had been examined earlier together with the conformational preferences of 3,9-dimethoxy-4-prenylpterocarpan (bitucarpin A) at the same level (Phys. Chem. Chem. Phys. 2004, 6, 2849). The structure with O5 trans with respect to H6a (O(t)) is about 2 kcal/mol less stable in vacuo than that with one of the H6 trans to it (H(t)); in aqueous solution its energy gap is nearly conserved. The in-in arrangement of the hydroxyl groups of erybraedin turns out to be preferred in vacuo (even considering zero point and thermal effects), where pseudo H-bonds are formed between hydroxy hydrogens and pi electron distributions of prenyl groups. The continuum solvent effect (water) at the IEF-PCM/B3LYP/6-31G level on the relative stability of the various rotamers is very limited both on bitucarpin and erybraedin. Considering the dihydrated derivatives, significant differences in the solvation energy are found between the distinct hydration sites, increasing in the order: methoxy O, ring O, hydroxy O, and hydroxy H. In hydroxy-water interactions, in fact, water prefers to behave as an H-bond acceptor unless nearby bulky groups prevent its approach. Interestingly enough, a bridging water molecule between the hydroxy H of erybraedin and the prenyl group can be found. The inclusion of BSSE corrections in hydroxy-water interactions decidedly favors out-out hydrated arrangements, followed by out-in and in-out ones. Bulk solvent effects with IEF-PCM about the dihydrated systems almost invert the stability order found in vacuo. When a four-water cluster is considered using QM methods, waters gather in H-bonded pairs around the solute OH groups. MD simulations, carried out on a pterocarpan solute (J. Phys. Chem. B 2005, 109, 16918), supply water adducts consistent with a liquid state that have also been embedded in the continuum solvent.
在真空条件下,已在B3LYP/6 - 31G水平上考虑了3,9 - 二羟基 - 4,8 - 二异戊烯基紫檀素(刺桐定C)稳定二面角值的所有可能组合,其羟基向外构象已与3,9 - 二甲氧基 - 4 - 异戊烯基紫檀素(比图卡品A)在相同水平下的构象偏好一起进行了早期研究(《物理化学化学物理》,2004年,6卷,2849页)。相对于H6a,O5呈反式(O(t))的结构在真空中比其中一个H6与其呈反式(H(t))的结构稳定性低约2千卡/摩尔;在水溶液中,其能量差距几乎保持不变。刺桐定羟基的向内排列在真空中(即使考虑零点和热效应)被证明是更有利的,此时羟基氢与异戊烯基的π电子分布之间形成了假氢键。在IEF - PCM/B3LYP/6 - 31G水平上,连续介质溶剂效应(水)对比图卡品和刺桐定各种旋转异构体相对稳定性的影响非常有限。考虑二水合衍生物时,不同水合位点的溶剂化能存在显著差异,按以下顺序增加:甲氧基O、环O、羟基O和羟基H。实际上,在羟基 - 水相互作用中,除非附近有庞大基团阻止其靠近,水更倾向于作为氢键受体。有趣的是,可以发现刺桐定的羟基H与异戊烯基之间存在一个桥连水分子。在羟基 - 水相互作用中包含基组重叠误差校正明显有利于向外 - 向外水合排列,其次是向外 - 向内和向内 - 向外排列。使用IEF - PCM对二水合体系的体相溶剂效应几乎颠倒了在真空中发现的稳定性顺序。当使用量子力学方法考虑四水簇时,水在溶质OH基团周围以氢键对的形式聚集。对一种紫檀素溶质进行的分子动力学模拟(《物理化学杂志B》,2005年,109卷,16918页)提供了与液态一致的水加合物,这些加合物也已嵌入连续介质溶剂中。