应用分子表面静电势预测羟基甲氧基多溴联苯醚正辛醇-空气分配系数.docVIP

应用分子表面静电势预测羟基/甲氧基多溴联苯醚的正辛醇-空气分配系数 Prediction of octanol-air partition coefficient of OH/MeO-PBDE congeners using molecular surface electrostatic potentials 徐文斌1 王 维2 许惠英1,＊ (1浙江树人大学生物与环境工程学院, 杭州 310015; 2镇海环境监测站, 宁波 315200; E-mail: xuhy65@163.com) 1 WANG Wei2 XU Hui-Ying1,* (1College of Biology Environment Engineering, Zhejiang Shuren University, Hangzhou 310015; 2Zhenhai Environmental Monitoring Station, ningbo 315200; E-mail: xuhy65@163.com) 引言 羟基多溴联苯醚（OH- PBDEs）和甲氧基多溴联苯醚（MeO-PBDEs）（PBDEs）OH/MeO-PBDE化合物材料与方法 首先，用HyperChem7.0软件构建并用其内置的AM1半经验量子化学方法，优化29个羟基甲氧基多溴联苯醚分子的初始几何参数。其后，用Gaussian98对各个分子在HF/6-31G*水平上进一步优化，并在优化结构基础上提取结构参数和进行静电势的计算。对29个羟基甲氧基多溴联苯醚分子的表面积（S）和体积（Vmc）进行了计算。在所计算的参数中与本研究有关的是：Vmc为分子体积；Vs, max ：分子表面上最正的静电势；Vs, min ：分子表面上最负的静电势。最后用统计软件SPSS.10.0中的多元线性回归分析模块建立起羟基甲氧基多溴联苯醚分子结构与正辛醇-空气分配系数的定量关系。 3．结果与讨论羟基甲氧基多溴联苯醚分子的正辛醇-空气分配系数取自文献采用逐步回归方法建立了羟基甲氧基多溴联苯醚分子结构与正辛醇-空气分配系数的定量关系logKOA = 3.334×10-2×Vmc + 0.556×Vs, max + 1.514×Vs, min + 3.990 N = 29 R = 0.952 SD = 0.3617 F = 79.962 式中，N为样本数，R为复相关系数，SD为标准偏差，F为相关关系式的F分布值。 我们发现，影响羟基甲氧基多溴联苯醚logKOA的最主要因素为分子体积(Vmc)（R 为 0.884）。通常而言，分子越大，分子与分子间的非反应性色散作用越强，更易分配在正辛醇有机相中，因而OH/MeO-PBDE化合物logKOA值大。此外，中还引入了2个分子静电势参数：Vs, max 和 Vs, min 。前者表示分子表面上最正的静电势分子表面上最负的静电势。从文献定义和以往研究的情况看，Vs, min 体现了一个分子静电氢键贡献大小，反映了分子接受质子形成氢键的能力（即氢键碱度），其值越小（绝对值越大），接受质子形成氢键的能力越强。Vs, max 反映了分子提供质子形成氢键的能力（即氢键酸度），其值越大，提供质子形成氢键的能力越强。，Vs, min和Vs, max与正辛醇-空气分配系数呈正相关，表明OH/MeO-PBDE化合物越容易分配在中分子提供质子形成氢键的能力越强，OH/MeO-PBDE化合物越容易分配在中。为了检验模型的稳健性，本研究对模型中各变量的自相关程度进行了分析，得到模型各变量的变异膨胀系数（VIF） (Vmc)、1.120 (Vs, max) 和 1.742 (Vs, min)。变异膨胀因子（VIF）的取值大于等于1。当 =1.0，表示各变量间无自相关；当 =1.0~5.0，表示相关方程可以接受；当10，表示回归方程不稳定，需要重新检验。VIF值均，因此从统计学角度表明，所建模型是完全可以接受的。 图1 的预测值与实验值的关系 Fig.1 Relationships between predicted and experimental为了进一步验证模型的预测能力，用“留一法”做交叉验证分析，Rcv）达到0.908。图1清晰地显示了logKOA的实验值和预测值之间偏差很小，两者之间很好的一致性表明所建立的模型具有很强的预测能力。4. 结论 从上述对羟基甲氧基多溴联苯醚分子结构与正辛醇-空气分配系数的QSPR研究看，分子表面静电势参数结合分子可以好地用于定量描述羟基甲氧基多溴联苯醚的分子结构与性质之间的相关关系。结果表明，正辛醇-空气分配系数logK