硝基芳烃衍生物结构—毒性定量关系和致毒机理的理论研究

【摘要】： 硝基芳烃是一类重要的污染物，应用广泛且在环境中大量存在，对生态环境造成了严重污染，威胁到人们的生产、生活和身体健康，因此对硝基芳烃的毒性进行系统研究具有重要意义。定量结构-活性关系(Quantitative Structure-Activity Relationship，QSAR)是一种非常重要的毒性研究方法，已广泛用于硝基芳烃类化合物的毒性研究。目前，在建立硝基芳烃毒性QSAR模型过程中，普遍存在如何获得分子描述符的问题，先前的研究主要通过实验或采用半经验分子轨道计算获得，前者费时费力且误差较大，后者属于近似研究，结果欠准确。随着计算机和量子化学理论的飞速发展，考虑了电子相关的第一性原理方法密度泛函理论(Density Functional Theory，DFT)得以广泛应用，其计算结果较为精确，又比较节省机时，但迄今为止，尚未见有将DFT方法应用于硝基芳烃毒性研究的文献报道。此外，文献研究表明，2，4，6-三硝基甲苯(TNT)及其代谢产物可与血红蛋白(Hemoglobin，Hb)发生共价结合，且Hb加合物水平与动物的染毒剂量相关。其它硝基芳烃也有类似报道。但是至今亦未见有与此相关的动态理论计算工作报道。对硝基芳烃衍生物的致毒机理有待深入研究。 在此背景下，本论文围绕硝基芳烃衍生物的毒性进行了系统的理论计算研究，全部工作主要分为三章。 1．以硝基芳烃对黑呆头鱼的急性毒性研究为例，运用不同理论水平的量子化学计算方法，包括半经验AM1和PM3方法，从头算HF／6-31G~*和DFT-B3LYP／6-311G~(**)方法等计算28种硝基芳烃化合物的分子几何和电子结构，获得它们的分子结构描述符，结合硝基芳烃对黑呆头鱼的急性毒性实验值，建立了相应的QSAR模型。结果表明，基于DFT-B3LYP方法计算所得QSAR模型稳定性和显著性明显优于AM1和PM3方法；相比于HF方法，B3LYP模型能对硝基芳烃的致毒机理作出更为合理的解释。故就本例而言，B3LYP方法更适合于硝基芳烃急性毒性的QSAR研究。 2．运用DFT-B3LYP方法研究硝基芳烃的分子结构与其毒性之间的定量关系(QSAR)，寻找影响其毒性大小的关键结构因素，对同系物的已有毒性进行阐明、对未知毒性进行预测，并为硝基芳烃的毒性降解处理提供有意义的建议。 采用DFT-B3LYP方法，在6-311G~(**)基组水平下，对单(多)硝基苯、单(多)硝基甲苯、硝基苯胺、卤代硝基苯等多系列的硝基芳烃化合物进行分子几何全优化计算，获得相关量子化学参数作为分子描述符，结合硝基芳烃对梨形四膜虫、圆腹雅罗鱼和斜生栅列藻等不同物种的急性毒性实验数据，分类建立了毒性QSAR模型，讨论了硝基、氨基和卤素等取代基对化合物毒性的影响。结果表明，硝基是硝基芳烃衍生物的主要致毒基团，硝基数越多，化合物的毒性越大。单硝基取代芳烃和多硝基取代芳烃具有不同的致毒机理，前者为极性麻醉型毒性化合物，其毒性大小由疏水性和电子反应性共同决定，毒性作用主要为对细胞膜的穿透刺激作用；后者为反应性毒性化合物，其毒性大小主要由电子反应性决定，分子可能发生单电子还原或亲核取代等反应而致毒。苯环上的其它取代基视其电子特性不同对硝基芳烃的毒性影响也不同，给电子基团如氨基会削弱硝基芳烃的毒性；吸电子基团如卤素则会增强化合物的毒性。故采用合适的包覆剂将硝基包覆或用微生物对其进行降解应该是治理硝基芳烃类化合物污染的有效途径。文中给出的稳健QSAR模型可用于相应化合物对不同物种急性毒性的预测。 3．在前两章研究的基础上，以乙硫醇(C_2H_5SH)作为含还原性巯基的生物蛋白大分子(R-SH)的简单模型物，对TNT及其代谢产物、硝基苯和二硝基苯以及卤代硝基芳烃等不同类型的硝基芳烃的还原代谢中间体亚硝基芳烃与乙硫醇的反应，进行了量子化学动态计算，寻找反应过渡态，求得反应活化能，从理论上探讨了硝基芳烃的致毒机理。 以DFT-B3LYP／6-31G~*方法全优化计算了TNT、2-氨基-4，6-二硝基甲苯和4-氨基-2，6-二硝基甲苯在生物体内的代谢还原中间体亚硝基芳烃及其与乙硫醇反应各驻点的结构，求得质子化亚硝基芳烃中亚硝基氮与巯基硫间距离分别为2.005(?)、1.935(?)和1.955(?)时的反应过渡态(TS1，TS2和TS3)以及相应的活化能15.443 kJ／mol、36.286kJ／mol和45.994 kJ／mol。表明硝基数越多，反应越易进行；氨基具致钝作用；当亚硝基处于甲基邻位时，因受到位阻作用，故反应活化能增高。 以类似方法和步骤计算研究了硝基苯、1，3-二硝基苯和1，4-二硝基苯的亚硝基化合物与乙硫醇的反应机理，获得类似结果和结论。求得氮-硫原子间距分别等于1.966(?)、2.045(?)和2.054(?)的三个反应过渡态和相应的活化能35.707kJ／mol、17.783kJ／mol和12.702kJ／mol。表明反应速度受硝基数目和位置的影响：硝基数越多反应活化能越小；对二硝基苯快于间二硝基苯。三个反应活化能的大小顺序与硝基芳烃毒性大小排序恰好相反，预示硝基芳烃与生物蛋白大分子中还原性巯基间的反应可能是其致毒的关键步骤。这与第二部分QSAR研究结果相吻合。 最后，以类似方法和步骤对卤代硝基芳烃的致毒机理进行了计算研究。求得4-氟硝基苯、4-氯硝基苯和4-溴硝基苯的亚硝基代谢产物与乙硫醇的反应过渡态和相应的活化能。过渡态的氮-硫原子间距分别为1.940(?)、1.947(?)和1.941(?)，活化能依次为43.397kJ／mol、41.529kJ／mol和44.332kJ／mol。不同卤素在同一位置取代对反应速度几无影响，表明卤代硝基芳烃至少应存在两种致毒反应：一是其被还原代谢前，与卤素相连的苯环碳原子先受到亲核试剂进攻，使卤素原子被取代；二是其亚硝基代谢中间体与含巯基的蛋白质大分子发生了共价加合，形成了加合物。 总之，本文用第一性原理DFT-B3LYP方法，对硝基芳烃的分子结构与毒性的定量关系进行了系统研究，找到影响硝基芳烃毒性的关键因素，给出了较好的毒性预测模型(QSAR)，并提出了硝基芳烃污染物毒性降解的初步建议。首次对硝基芳烃的致毒机理进行了动态理论计算，建议了硝基芳烃在生物体内与蛋白质作用的简单模型和致毒的具体过程，找到关键基元反应的过渡态和活化能，探讨了取代基效应和位阻效应。如上关于硝基芳烃衍生物结构-毒性定量关系和致毒机理的理论研究，为理论化学、材料化学和环境化学的交叉研究和发展，提供了新的示例和结果，提供了丰富的原始数据和相应规律。