分子价连接性指数[14]是被公认为应用最广、最优秀的一类拓扑指数,是由Randic于1975年提出,并由Kier和Hall扩展的第一类拓扑指数,从而形成了分子连接指数系列。分子连接性指数都是用顶点度计算,用分子中所有边两端顶点的顶点度乘积平方根倒数的和,用公式表示为:
mXtV = Σ(δiv δjvδkv…)-0.5
上式中的m表示阶数, t表示拓扑指数子图的类型。子图类型t的取值有四种,对应四种子图类型,它们分别为:
①路径项:路径项以t=p表示,即mXp。路径项对应“直链”的子图。
②簇项:路径项为“直链”子图,当“直链”子图中出现一条或两条分支时,就用簇项来表示。簇项以t=c表示,即mXc。簇项对应的子图中至少有一个顶点与三条边或四条边相连,且所有顶点不与两条边相连。
③路径/簇项:上述的路径项描述纯粹的直链,簇项描述纯粹的分支,但在有机化合物中,直链和分支结构经常是混杂在一起的,为了描述这种结构,分子连接性指数还包含了路径/簇项。用t=pc表示,即 mXpc。它们的子图可以看作是簇项结构中增加了一条或若干条边。
④链项:针对环状化合物的特点,分子连接性指数中定义了链项来专门描述环状化合物。链项以t=ch表示,即mXch。p、c、pc及ch见图1, 分别与a、b、c、d对应。
分子结构的基本子图类型
1.2.2 原子型电拓扑态指数
根据分子中的隐氢图,由于每个原子处于不同的拓扑状态,其成键电子状态也不一样,据此计算得到的所有非氢原子的一组不变量即为原子型电拓扑态指数(Ek) [15]。Ek包含两个方面:一方面是由每个非氢原子类型k自身形成的原子特征值,以“Hk”表示;另一方面是被其它非氢原子作用所产生的特征值增量,用“ΔHk”表示。定义原子类型k的电拓扑态指数Ek为:
Ek=∑(Hk+ΔHk)j 文献综述
上式中的j是原子类型k的数目。
为了研究多氟有机化合物分子结构与其蒸汽压之间的构效关系,由Chem
DrawUltra建立35种多氟有机化合物的分子结构,以mol文件的形式保存,再用Matlab调出其结构,最后,采用修正的程序[15]计算得到得到原子型电拓扑态指数(Ek)、Kier分子价连接性指数(mXtV)。
1.3 试验方法
将用上述方法计算得到的12个分子连接性指数、46个原子电性拓扑指数,总共56个参数设置成自变量,因变量是对应的蒸汽压( logVP),应用最佳变量子集回归法筛选变量,选择其中最好的变量组合,得到QSPR预测模型。
2结果讨论
2.1 多氟有机化合物QSPR模型的构建
将上述计算得到的35种多氟有机化合物的两类拓扑指数,总共56个参数和对应的蒸汽压一起输入SPSS 13.0计算机软件,置信区间在95%内进行统计回归分析,应用最佳变量子集回归法筛选变量,选择其中最好的变量组合。回归所得到的各个QSRR模型参数见表1。回归方程的相关系数表示为R,判定系数表示为R 2,校正判定系数表示为R 2adj,估计标准误差表示为S,Fischer检验值表示为F。