取代苯类化合物对绿藻的结构毒性定量关系研究Word格式.docx
- 文档编号:8165824
- 上传时间:2023-05-10
- 格式:DOCX
- 页数:28
- 大小:104.64KB
取代苯类化合物对绿藻的结构毒性定量关系研究Word格式.docx
《取代苯类化合物对绿藻的结构毒性定量关系研究Word格式.docx》由会员分享,可在线阅读,更多相关《取代苯类化合物对绿藻的结构毒性定量关系研究Word格式.docx(28页珍藏版)》请在冰点文库上搜索。
张爱茜等[15]用量子化学参数法进行了氯代芳香族化合物对羊角月牙藻(绿藻的一种)的毒性和QSAR分析的研究。
2取代苯类对绿藻的毒性及结构参数的数据来源及意义
2.1材料与实验方法
实验藻种是由中国科学院武汉水生所提供的斜生栅列藻(ScendesmusOblligquus),室内培养。
四十取代苯类化合物均为分析纯。
实验按照OECD藻类阻碍生长实验标准进行,实验数据取自文献[16]。
实验器皿选用250ml锥形瓶,每个化合物设六组(包括空白),作一平行。
向每个锥形瓶中加入30ml两倍浓度的培养液,按对数间距0.2的浓度梯度加入测试化合物的标准溶液,摇匀后,加入1ml约1*104个/ml处于对数生长初期的藻种总体积用蒸馏水补充至60ml。
实验样品置于光照培养箱中静止培养,温度与酸碱度如下:
温度为20±
1。
;
PH:
7.2±
0.2;
照明:
明暗对照:
12h:
12h,平均光照强度为4000lx。
将受试化合物浓度的对数与绿藻生长阻碍率作一元线性回归,求出半抑制浓度EC50(mol/l)值,列于表2.1-1。
2.2参数计算
使用量子化学方法计算得到所研究化合物的分子最低空轨道能(Elumo)和分子最高占有轨道能(Ehomo)、范德华面积(Svdw)及分子量(Mw)。
使用软件SRC-lgKOWforWindows得到分子辛醇/水分配系数的对数lgP(实验值)所研究化合物的结构参数见表2.1-1。
化合物的结构参数数据由陆老师提供。
量子化学各参数意义如下:
Elumo:
分子最低空轨道能,它描述的是分子与亲核试剂之间的相互作用的能力,因此直接与分子对电子的亲合力有关;
Ehomo:
分子最高占有轨道能,都是负值,绝对值越大则分子该轨道上的电子越稳定;
它与分子电离势相关,可以作为分子给出电子能力的量度,Ehomo越高,分子给电子能力越强;
Elumo和Ehomo都是反映化合物发生由轨道控制的反应的倾向[17]。
Svdw:
分子的范德华面积,它可以定量反映分子的空间形状和大小,与分子量Mw一样,可以作为立体效应参数。
lgP:
辛醇/水分配系数,它是宏观疏水性参数,能反映化合物脂溶性的大小,lgP值越高,疏水性越强,化合物越容易富集到生物的体内。
表2.1-1取代苯类化合物对绿藻的毒性及结构参数
序
号
化合物
Ehomo
(ev)
Elumo
Svdw
(A2)
Mw
lgP
lgl/EC50
(mol/L)
1
间硝基氯苯
-10.02
-1.15
82.14
157.56
2.49
3.95
2
对硝基氯苯
-10.07
-1.18
82.88
2.35
4.01
3
3,4—二氯硝基苯
-10.17
-1.62
91.43
192.00
3.12
4.52
4
2,5—二氯硝基苯
-10.13
-1.60
90.93
3.09
4.31
5
邻氯苯胺
-9.17
-0.19
74.14
127.57
1.90
2.89
6
间氯苯胺
-9.16
-0.21
75.34
1.88
2.79
7
2,4—二氯苯胺
-9.28
-0.45
85.17
162.02
2.78
3.74
8
2,5—二氯苯胺
-9.86
-0.62
85.25
2.75
3.82
9
2,4,6—三溴苯胺
-9.87
-0.71
106.91
329.82
4.03
4.37
10
3—溴苯胺
-9.14
79.43
172.03
2.10
2.80
11
4—氨基苯乙醚
-8.92
-0.09
83.14
137.18
1.24
3.10
12
间溴硝基苯
-9.98
-1.42
86.33
202.01
2.64
4.32
13
对溴硝基苯
-10.01
-1.50
86.45
2.55
3.88
14
邻硝基苯胺
-9.54
-1.22
75.71
138.13
1.85
3.33
15
对硝基苯胺
-9.31
-1.36
76.38
1.37
3.48
16
2,4—硝基苯胺
-9.78
-1.32
89.59
183.12
1.84
4.68
17
邻甲苯胺
-9.01
0.04
71.30
107.15
1.32
2.34
18
苯胺
-0.13
64.37
93.13
0.90
2.56
19
2,4—二氯酚
-9.38
83.75
163.00
3.06
3.62
20
2,4,6—三氯酚
-9.51
-0.66
94.02
197.45
3.69
3.81
21
五氯酚
-10.12
-1.27
110.13
266.34
5.04
4.63
22
间硝基苯酚
-1.38
74.93
139.11
2.00
3.75
23
对硝基苯酚
-9.55
-1.33
76.61
1.91
3.57
24
苯酚
-9.06
0.13
62.66
94.11
1.46
2.46
25
硝基苯
-9.88
-1.06
71.92
123.11
3.26
续表2.1-1取代苯类化合物对绿藻的毒性及结构参数
27
2,6—二硝基甲苯
-10.33
-1.19
90.34
182.14
4.06
28
对硝基甲苯
-9.89
78.22
137.14
2.37
29
邻二硝基苯
-10.45
-1.63
83.07
168.11
1.69
30
间二硝基苯
-10.34
-1.31
85.26
1.49
4.85
31
2—硝基苯甲醚
-9.47
-0.96
83.67
153.14
1.73
3.44
32
4—硝基苯甲醚
83.68
2.03
3.65
33
对二硝基苯
-10.18
-1.75
84.31
4.96
34
邻硝基氯苯
-10.04
-1.11
82.05
2.26
3.94
35
3—硝基苯甲醚
-9.81
-1.09
84.48
2.16
3.71
36
-9.36
76.56
1.39
3.40
37
邻硝基苯酚
-9.40
-1.29
73.83
1.79
3.51
38
2,4—二硝基苯酚
-9.96
-1.52
87.38
184.11
1.67
4.16
39
2,3—二氯苯胺
-9.80
-0.46
83.90
2.82
3.98
40
对甲苯胺
-8.94
-0.02
71.37
3.19
3统计分析过程及结果
为了建立最后的线性回归模型,需要对化合物进行分类这样能在每一组化合物内建立线性回归模型以便使其可能具有比四十个化合物样本一起建立的线性回归模型具有更好的拟合程度。
但是,在此次研究中分子结构参数五个,为了使分类的依据更加科学有效需要对变量进行因子分析以提出主因子作为分类的依据。
在此之前首先对化合物的五个结构参数与半抑制浓度进行相关分析求出相关系数矩阵,以观察各个变量的相关性,为进一步的分析做准备。
所有的分析(包括最后线性模型)均用SPSS11.0完成。
3.1相关分析
五个化合物结构参数(Ehomo、Elumo、Svdw、Mw与lgP)与半抑制浓度之间的相关系数矩阵见表3.1-1。
从表3.1-1中我们可以看出Svdw与Mw的相关程度最高,这是比较容易解释的:
取代苯类化合物有同一个母体即苯环,分子量大往往既意味着分子比较大自然作为反映分子空间形状的范德华面积也应该较大。
表3.1-1还反映出一个问题:
Ehomo与Elumo的相关程度较高,这也容易理解:
这两个参数都与分子的轨道能有关它们之间应该有某种相关性。
值得注意的是分子最高占有轨道能和分子最空轨道能以及范德华面积与半抑制浓度的相关程度很好,这可以理解为取代苯类化合物对生物产生毒性的反应大都是由轨道控制的反应;
此外,为凡是化学反应都是通过电子的得失与成键作用完成的,因此表征分子对电子控制能力的分子最高占有轨道能与分子最低空轨道能这两个参数就与表示
表3.1-1取代苯类化合物结构与毒性的相关系数矩阵
变量
1.000
0.770
-0.600
-0.545
-0.364
-0.862
-0.390
-0.362
-0.143
-0.748
0.932
0.786
0.719
0.778
0.656
0.413
生物毒性的半抑制浓度有密切的关系。
而范德华面积则表征着分子在空间的大小也对决定化合物分子在生物体内的输移起着重要的作用,因此,Svdw与反应化合物生物毒性的半抑制浓度应该有较大的关系。
最后,表2还反应出辛醇/水分配系数与Svdw和Mw有相当高的相关程度。
3.2因子分析
因子分析是从研究相关矩阵内部的依赖关系出发,把一些具有错综复杂关的变量归结为少数几个综合因子的方法[18]。
因子分析的基本思想是根据相关性大小把变量分组,使得同组内的变量之间相关性较高,但不同组内的变量相关性较低。
每组变量代表一个基本结构,这个基本结构称为公共因子。
对于所研究的问题就可以用最少个数的不可测的所谓公共因子的线性函数与特殊因子之和来描述原来观测的每一个分量。
因子分析还可以用来对变量或样品进行分类处理。
可根据因子得分值,在因子轴所构成的空间中把变量或样品点画出,以达到分类的目的。
在本次研究中将化合物的分子最低空轨道能、分子最高占有轨道能、范德华面积、分子量以及分子的辛醇/水分配系数的对数作为原始变量输入。
经过计算后得到总方差分析结果见表3.2-1,而因子得分系数矩阵见表3.2-2;
由于原始的得分系数矩阵意义不明显,故对得分系数矩阵进行旋转,旋转后的得分系数矩阵见表3.2-3。
从总方差分析的结果表3.2-1中可以得到公因子和特征值,五个特征值由大到小依次排列,第一特征值与最大公因子相对应,其余依次类推。
根据公因子方差累积百分比大于百分之八十的原则,从五个公因子中提取前两个做主因子是完全可以的,且可以概括绝大部分信息。
表3.2-1总方差分析
因子
序号
初始特征值
旋转方差贡献总和
总方差
方差贡献率(%)
累积方差贡献率(%)
1
3.322
66.446
2.640
52.795
2
1.193
23.856
90.302
1.875
37.506
3
0.223
4.463
94.768
0.198
3.952
98.720
5
0.064
1.280
100.000
表3.2-2初始因子得分系数矩阵
变量
Ehomo
-0.783
0.522
-0.608
0.737
0.938
0.231
0.918
0.271
Lgp
0.785
0.500
选取主因子后,计算得到初始因子载荷(表3.2-2)可以看出初始因子载荷各个因子上都比较均匀,这样即难以得到主因子的物理意义对以后的分析是十分不利的;
因此,将载荷进行旋转,初始因子载荷经过方差最大正交旋转后,每个因子的载荷会向两极分化,即每个公因子只在少数几个因子上有高载荷(见表3.2-3)。
表3.2-3旋转后因子得分矩阵
因子
-0.350
0.873
-0.084
0.953
0.904
-0.342
0.910
-0.297
0.930
-3.24E-02
从表3.2-3中可以看出:
第一主因子在范德华面积、分子量、辛醇/水分配系数上的载荷很高;
而第二主因子在分子最低空轨道能、分子最高占有轨道能这两项上的载荷很高。
因此,可以认为第一主因子主要代表化合物的空间大小和形状(lgP与Svdw以及Mw关系很大——这一点可以从三者的相关系数很大看出)。
第二主因子可以认为是主要与分子的电性化参数有关可以反映化合物分子的得失电子的能力。
各个化合物的因子得分数据见表3.2-4,各化合物因子得分的绝对值越高表明该因子代表的化学性质在该化合物上表现的越强烈,因此因子分析有助于对化合物进行分类。
根据各化合物的因子得分作出因子的得分散点图以初步观察化合物的分类情况,因子得分散点图见图3.2-1。
从图3.2-1中可以看到多数的化合物的因子得分的绝对值都在1以内,只有四至五个点的绝对值超过1,而其中有三个点离群的很厉害,因此,可以预期在聚类分析中这三种化合物会被聚类出大组而独立。
表3.2-4化合物因子得分系数表
因子1
因子2
0.00867
-0.49653
-0.05498
-0.62177
0.77899
-0.99918
0.75696
-0.94013
-0.39537
1.34344
-0.36815
1.32709
0.65094
1.18311
0.50692
0.34062
3.33622
0.67841
0.27707
1.45988
-0.26874
1.56855
0.51827
-0.68179
0.44272
-0.82735
-0.69276
-0.21852
-0.94425
-0.24202
0.11847
-0.56621
-0.87263
1.61296
-1.52507
1.30612
0.74289
1.14492
1.64530
0.95336
3.22519
-0.00699
-0.67248
-0.17862
-0.66103
-0.33575
-1.19622
1.68485
-0.93931
-0.42976
26
2,4—二硝基甲苯
0.14302
-1.17882
0.25528
-0.97002
-0.43512
-0.65262
-0.51626
-1.70788
-0.41307
-1.28169
-0.23829
0.12739
-0.29631
-0.77391
对二
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 取代 化合物 绿藻 结构 毒性 定量 关系 研究