论文部分内容阅读
随着化学工业的发展,越来越多新型化学品进入了人类的生活。根据美国化学文摘(Chemical Abstracts Service, CAS)的最新统计数据,目前注册的已发现和合成的化学物质总数已超过1亿种,而其中很多物质对环境和人类健康是有害的。面对数目如此庞大的化学物质,获得其毒性信息从而进行合理的监管和使用,是保护人类和环境免受危害的根本方法。然而根据统计,目前处于监管状态下的化学品仅有30多万种,绝大多数化学品缺少毒性信息,存在着巨大的数据缺口。而传统生物实验方法耗时长、花费大,无法满足大量化学品急需测试的需要,且与目前提倡的动物保护“3R(替代(Replacement)、减少(Reduction)和优化(Refinement))原则”相悖,因此能快速高效地获取化学品生物毒性数据的计算机模拟预测方法一定量构效关系(Quantitative structure-activity relationship, QSAR)方法已被越来越多的研究者所接受和使用。丰富的实验数据是进行QSAR研究的基础,只有在数据充足的前提下才能够借助计算机、应用统计学方法构建出合理的构效关系方程,然后结合生物化学及物理化学等相关知识,总结出化合物性质变化和分子结构特征变化之间的联系及规律,进而通过已知的分子结构特征对化合物的未知性质进行预测。因此,QSAR的研究对象应该是具有丰富实验数据的性质变量。目前在环境及生态毒理学领域常见的性质变量有急性毒性、慢性毒性、致癌性和皮肤敏感性等。本论文选择其中的急性毒性和致癌性进行了QSAR研究。QSAR研究成功的一个重要因素是获取合理准确的分子结构及相应的结构特征变量,为此需要对分子结构进行合理的几何优化。本论文中的分子结构优化都是在密度泛函理论(Density Functional Theory, DFT)方法的MPWB1K/6-31+G(d,p)水平下进行的,并且通过研究OH自由基引发的灭多威在大气中的降解机制对该水平的优化方法进行了实际效用的验证,在此基础上开展了有机物对啮齿动物致癌性和有机物对水生生物急性毒性的研究。本论文的主要研究内容及成果总结如下：(1)基于DFT方法研究OH自由基引发的灭多威在大气中降解机理灭多威是一种氨基甲酸酯类农药,主要用作广谱杀虫剂,但其对动物和人类也会产生毒害,因此研究其在大气中的降解机制具有重要的实际意义。本文研究了灭多威两种同分异构体在大气中O2、NO和H20存在的情况下由OH自由基引发的降解反应。在DFT方法的MPWB1K/6-31+G(d,p)水平下对反应物、中间体和产物进行了几何结构优化和频率计算,构建了反应的势能剖面图并讨论了所有可能的降解路径。采用小曲率隧道效应(small-curvature tunneling,SCT)校正的正则变分过渡态理论(Canonical variational theory, CVT)计算了每条反应路径在200-400K温度范围内基元反应的速率常数。研究发现,灭多威两种同分异构体的反应状况非常相似,既可以反生OH自由基加成反应,也可以发生H原子抽提反应。根据反应速率常数的分支比可以看出,对于Z式异构体,常温下抽提反应占较大比重,随着温度上升,加成反应比重上升。而对于E式异构体,常温下加成反应占较大比重,并随着温度上升比重不断增大。除此之外,还计算了两种同分异构体的大气寿命,目前还未见其他有关报道。(2)基于DFT方法研究有机物对啮齿动物致癌性的QSAR预测模型本论文基于CPDB肿瘤实验数据库,选择数据最多的肝脏肿瘤实验数据作为分析研究对象。通过对比雌雄两种性别啮齿生物的肿瘤实验数据,发现多种化合物存在明显差异,而且总体而言,雌性实验动物对致癌物表现出比雄性更强的耐受性。因此,综合考虑分子的结构特征,构建了有环化合物对雄性大鼠、无环化合物对雄性大鼠、有环化合物对雌性大鼠、无环化合物对雌性大鼠、有环化合物对雄性小鼠、无环化合物对雄性小鼠、有环化合物对雌性小鼠、无环化合物对雌性小鼠等八个定量预测模型。然后根据模型、结合已有知识对致癌机理进行了讨论。研究所构建的预测模型可以为研究性别对致癌作用的影响以及危险物质的管理提供重要的理论支持。(3)基于DFT方法研究有机物对水生生物急性毒性QSAR预测模型本章选择有机物对不同营养级典型水生生物(包括鱼类和水蚤)的急性毒性作为研究对象,基于DFT方法优化的分子结构分别构建了QSAR预测模型并进行了对比分析。在鱼类的急性毒性实验中,黑头呆鱼是最常用的实验对象之一,因此具有丰富的实验数据,为利用QSAR方法研究有机物的急性毒性预测模型提供了重要的数据支撑。本论文搜集整理了963种有机物对黑头呆鱼的96小时急性毒性实验数据,在DFT方法的MPWB1K/6-31+G(d,p)水平上对所有化合物的分子进行了几何结构优化,选择Dragon描述符和辛醇水分配系数logKow作为建模的备选变量,借助遗传算法和多元线性回归构建了目前数据集最大的黑头呆鱼急性毒性预测通用模型。除此之外,根据毒性比率及有机物毒性作用方式的划分,分别构建了麻醉性毒物和过量毒性毒物的急性毒性局域预测模型。大型蚤属于世界种,是毒性实验标准测试生物,具有丰富的实验数据。本章搜集整理了479种有机物对大型蚤的48小时急性毒性实验数据,同样在DFT方法的MPWB1K/6-31+G(d,p)水平上进行了化合物的几何结构优化,选择Dragon描述符及logKow作为建模的备选变量,借助遗传算法和多元线性回归构建了通用的急性毒性预测模型。除此之外，参考Verhaar区分框架,利用毒性比率将所研究化合物进行了毒性作用机制的分类并分别构建了QSAR模型,借助模型讨论了每种机制化合物的作用机理。最后,将有机物对黑头呆鱼和大型蚤呈现出的不同毒性作用机制进行了对比。研究发现,大型蚤对毒物的敏感性要高于黑头呆鱼；除了苯胺类物质,硫代磷酸酯类物质对黑头呆鱼和大型蚤的毒性也存在明显差异,其对黑头呆鱼仅表现为麻醉毒性,对大型蚤却呈现出过量毒性。