水消毒副产物卤代苯醌的结构与毒性分析

doi:10.13481/j.1671-587X.20230312

摘要/Abstract

摘要：

目的分析卤代苯醌（HBQs）对人肝癌HepG2细胞的半数抑制浓度（IC₅₀）与结构参数间的定量构效关系，探讨可能影响HBQs细胞毒性的结构参数。方法体外培养HepG2细胞，采用不同浓度HBQs溶液作用细胞，同时设空白对照组［10%胎牛血清（FBS）+DMEM培养基］和阴性对照组［10% FBS+DMEM培养基+1.33% 甲醇（MeOH）］。MTS法和中性红摄取试验（NRU）法确定HBQs的作用浓度，采用MTS法和NRU法检测各组细胞存活率并计算HBQs对人肝癌HepG2细胞的IC₅₀。采用量子化学计算模块MOPAC计算每种HBQs的结构参数。采用单因素线性回归分析方法，根据HBQs的IC₅₀和HBQs结构参数构建定量结构-毒性关系（QSTR）模型。结果 MTS法检测2，6-二氯-1，4-苯醌（2，6-DCBQ）的作用浓度为75、100、125、150、175、200、225和250 μmol·L^－1，2，6-二氯-3-甲基-1，4-苯醌（DCMBQ）的作用浓度为175、200、225、250、275、300和325 μmol·L^－1，2，3，6-三氯-1，4-苯醌（TCBQ）的作用浓度为200、225、250、275、300和325 μmol·L^－1，2，5-二溴-1，4-苯醌（2，5-DBBQ）的作用浓度为75、100、125、150、175和200 μmol·L^－1，2，6-二溴-1，4-苯醌（2，6-DBBQ）的作用浓度为200、225、250、275、300和325 μmol·L^－1；NRU法检测，2，6-DCBQ的作用浓度为25、50、75、100和125 μmol·L^－1，DCMBQ的作用浓度为125、150、175、200、225和250 μmol·L^－1，TCBQ的作用浓度为175、200、225、250、275和300 μmol·L^－1，2，5-DBBQ的作用浓度为75、100、125、150和175 μmol·L^－1，2，6-DBBQ的作用浓度为125、150、175、200、225和250 μmol·L^－1。各组HepG2 细胞存活率随HBQs浓度的升高明显降低，呈现明显剂量-反应关系。MTS法检测HBQs对HepG2细胞毒性作用为2，6-DCBQ>2，5-DBBQ>DCMBQ>TCBQ>2，6-DBBQ；NRU法检测HBQs对HepG2细胞毒性作用为2，6-DCBQ>2，5-DBBQ>DCMBQ>2，6-DBBQ>TCBQ。单因素多元线性回归分析，MTS法和NRU法检测HBQs 的IC₅₀与HBQs结构参数间相关性分析比较差异均无统计学意义（P>0.05）。结论卤代原子的取代位置、种类和数量可影响HBQs的毒性作用，且呈现卤素取代个数越多其毒性越小的趋势。

关键词: 卤代苯醌, 细胞毒性, 结构参数, 定量结构-毒性关系

Abstract:

Objective To analyze the quantitative constitutive relationship between the median-inhibitory concentration （IC₅₀） of halobenzoquinones （HBQs） to the human hepatocellular carcinoma HepG2 cells and structural parameters， and to explore the structural parameters that may affect the cytotoxicity of HBQs. Methods The HepG2 cells were cultured in vitro and cultured with different contentrations of HBQs solutions，at the same time， blank control group ［10% fetal bovine serum（FBS）+DMEM culture medium］ and negative control group ［（10% FBS+DMEM culture medium +1.33% methanol （MeOH）］ were set up. The action concentrations of HBQs were comfirmed with MTS and neural red uptake（NRU） methods，respectively.The survival rates of cells in various groups were detected by MTS and NRU methods， and the IC₅₀ of HBQs was calculated. The quantum chemical calculation module MOPAC was used to calculate the structural parameters of every kind of HBQs. Single-factor linear regression analysis was used to construct the quantitative structure-toxicity relationship （QSTR） models based on the IC₅₀ and structural parameters of the HBQs. Results The MTS method determination results showed that the action concentrations of 2，6-dichloro-1，4-benzoquinone（2，6-DCBQ） were 75， 100， 125， 150， 175， 200， 225， and 250 μmol·L^－1，the action concentrations of 2，6-dichloro-3-methyl-1，4-benzoquinone（DCMBQ） were 175， 200， 225， 250， 275， 300，and 325 μmol·L^－1， the action concentrations of 2，3，6-trichlorocyclohexa-1，4- benzoquinone（TCBQ） were 200， 225， 250， 275， 300， and 325 μmol·L^－1， the action concentrations of 2，5-dibromo-1，4-benzoquinone（2，5-DBBQ） were 75， 100， 125， 150， 175， and 200 μmol·L^－1，and the action concentrations of 2，6-dibromo-1，4-benzoquinone（2，6-DBBQ） were 200， 225， 250， 275， 300， 325 μmol·L^－1.The NRU method determination results showed that the action concentrations of 2，6-DCBQ were 25， 50， 75， 100， and 125 μmol·L^－1， the action concentrations of DCMBQ were 125， 150， 175， 200， 225，and 250 μmol·L^－1， the action concentrations of TCBQ were 175， 200， 225， 250， 275，and 300 μmol·L^－1，the action concentrations of 2，5-DBBQ were 75， 100， 125， 150，and 175 μmol·L^－1，the action concentrations of 2，6-DBBQ were 125， 150，175， 200， 225， and 250 μmol·L^－1.The survival rates of HepG2 cells in various groups were significantly decreased with the increase of HBQs concentration， showing a dose-response relationship. The cytotoxic effects of HBQs on the HepG2 cells detected by MTS method were 2，6-DCBQ>2，5-DBBQ>DCMBQ>TCBQ>2，6-DBBQ；and the cytotoxic effects of HBQs on HepG2 cells detected by NRU method were 2，6-DCBQ>2，5-DBBQ>DCMBQ>2，6-DBBQ>TCBQ.The single-factor multiple linear regression analysis results showed that the correlation analysis between IC₅₀ and structural parameters of the HBQs detected by MTS and NRU methods had no statistically significant differences（P>0.05）. Conclusion The substitution position， type and number of halogen atoms can affect the toxic effects of HBQs， and there is a trend that the more halogen substitution， the less toxic it is.

Key words: Halobenzoquinone, Cytotoxicity, Structural parameters, Quantitative structure-toxicity relationship

中图分类号:

R114

杜海英,姜永莉,韦英惯,李金华,卢日峰. 水消毒副产物卤代苯醌的结构与毒性分析[J]. 吉林大学学报(医学版), 2023, 49(3): 640-646.

Haiying DU,Yongli JIANG,Yingguan WEI,Jinhua LI,Rifeng LU. Analysis on structure and toxicity of water disinfection by-products halobenzoquinones[J]. Journal of Jilin University(Medicine Edition), 2023, 49(3): 640-646.

图/表 7

图1

表1

图2

表2

表3

表4

表5

参考文献 22

1	HUANG R F， WANG W， QIAN Y C， et al. Ultra pressure liquid chromatography-negative electrospray ionization mass spectrometry determination of twelve halobenzoquinones at ng/L levels in drinking water［J］. Anal Chem， 2013， 85（9）： 4520-4529.
2	WANG W， QIAN Y C， BOYD J M， et al. Halobenzoquinones in swimming pool waters and their formation from personal care products［J］. Environ Sci Technol， 2013， 47（7）： 3275-3282.
3	QIN F， ZHAO Y Y， ZHAO Y L， et al. A toxic disinfection by-product， 2， 6-dichloro-1， 4-benzoquinone， identified in drinking water［J］. Angew Chem Int Ed Engl， 2010， 49（4）： 790-792.
4	LI J H， WANG W， MOE B， et al. Chemical and toxicological characterization of halobenzoquinones， an emerging class of disinfection byproducts［J］. Chem Res Toxicol， 2015， 28（3）： 306-318.
5	BULL R J， RECKHOW D A， LI X F， et al. Potential carcinogenic hazards of non-regulated disinfection by-products： haloquinones， halo-cyclopentene and cyclohexene derivatives， N-halamines， halonitriles， and heterocyclic amines［J］. Toxicology，2011， 286（1/2/3）： 1-19.
6	LI J H， MOE B， VEMULA S， et al. Emerging disinfection byproducts， halobenzoquinones： effects of isomeric structure and halogen substitution on cytotoxicity， formation of reactive oxygen species， and genotoxicity［J］. Environ Sci Technol， 2016， 50（13）： 6744-6752.
7	ZUO Y T， HU Y， LU W W， et al. Toxicity of 2， 6-dichloro-1， 4-benzoquinone and five regulated drinking water disinfection by-products for the Caenorhabditis elegans nematode［J］. J Hazard Mater， 2017， 321： 456-463.
8	LECURIEUX F， MARZIN D， ERB F. Study of the genotoxic activity of five chlorinated propanones using the SOS chromotest， the Ames-fluctuation test and the newt micronucleus test［J］. Mutat Res， 1994， 341（1）： 1-15.
9	PROCHÁZKA E， ESCHER B I， PLEWA M J， et al. In vitro cytotoxicity and adaptive stress responses to selected haloacetic acid and halobenzoquinone water disinfection byproducts［J］. Chem Res Toxicol， 2015， 28（10）： 2059-2068.
10	PLEWA M J， WAGNER E D， RICHARDSON S D， et al. Chemical and biological characterization of newly discovered iodoacid drinking water disinfection byproducts［J］. Environ Sci Technol， 2004， 38（18）： 4713-4722.
11	ATTENE-RAMOS M S， WAGNER E D， PLEWA M J.Comparative human cell toxicogenomic analysis of monohaloacetic acid drinking water disinfection byproducts［J］. Environ Sci Technol， 2010， 44（19）： 7206-7212.
12	KARELSON M， LOBANOV V S， KATRITZKY A R. Quantum-chemical descriptors in QSAR/QSPR studies［J］. Chem Rev， 1996， 96（3）： 1027-1044.
13	MÉNDEZ-HERNÁNDEZ D D， TARAKESHWAR P， GUST D， et al. Simple and accurate correlation of experimental redox potentials and DFT-calculated HOMO/LUMO energies of polycyclic aromatic hydrocarbons［J］. J Mol Model， 2013，19（7）：2845-2848.
14	YEOM Y， HAN J， ZHANG X，et al.A review on the degradation efficiency，DBP formation，and toxicity variation in the UV/chlorine treatment of micropollutants［J］. Chem Eng J， 2021， 424： 130053.
15	LIU J Q， ZHANG X R. Comparative toxicity of new halophenolic DBPs in chlorinated saline wastewater effluents against a marine alga： halophenolic DBPs are generally more toxic than haloaliphatic ones［J］. Water Res， 2014， 65： 64-72.
16	ALAMIERY A. Corrosion inhibition effect of 2-N-phenylamino-5-（3-phenyl-3-oxo-1-propyl）-1， 3， 4-oxadiazole on mild steel in 1 M hydrochloric acid medium： insight from gravimetric and DFT investigations［J］. Mater Sci Energy Technol， 2021， 4： 398-406.
17	HUNG S， MOHAN A， RECKHOW D A， et al. Assessment of the in vitro toxicity of the disinfection byproduct 2， 6-dichloro-1， 4-benzoquinone and its transformed derivatives［J］. Chemosphere， 2019， 234： 902-908.
18	YANG M T， ZHANG X R. Comparative developmental toxicity of new aromatic halogenated DBPs in a chlorinated saline sewage effluent to the marine polychaete Platynereis dumerilii［J］. Environ Sci Technol， 2013， 47（19）： 10868-10876.
19	RYBICKA-JASIŃSKA K， VUIIEV V I. Molecular electrets-Why do dipoles matter for charge transfer and excited-state dynamics？［J］. J Photochem Photobiol A Chem， 2020， 401： 112779.
20	CHAN K T， JENSEN N， O’BRIEN P J. Structure-activity relationships for thiol reactivity and rat or human hepatocyte toxicity induced by substituted p-benzoquinone compounds［J］. J Appl Toxicol， 2008， 28（5）： 608-620.
21	SIRAKI A G， CHAN T S， O'BRIEN P J. Application of quantitative structure-toxicity relationships for the comparison of the cytotoxicity of 14 p-benzoquinone congeners in primary cultured rat hepatocytes versus PC12 cells［J］. Toxicol Sci， 2004， 81（1）： 148-159.
22	LIU Z Y， DANG K， GAO J H， et al. Toxicity prediction of 1， 2， 4-triazoles compounds by QSTR and inter species QSTTR models［J］. Ecotoxicol Environ Saf， 2022， 242： 113839.

HBQs	IC₅₀［c_B/（μmol·L^－1）]	95% Confidence interval (95%CI)
2,6-DCBQ	132.40±0.02	123.10－142.40
DCMBQ	219.70±0.01	212.00－227.60
TCBQ	229.30±0.01	221.80－237.00
2,6-DBBQ	274.40±0.01	263.10－286.20
2,5-DBBQ	148.10±0.01	141.60－154.80

HBQs	IC₅₀［c_B/（μmol·L^－1）]	95%CI
2,6-DCBQ	72.40±0.02	66.10－79.40
DCMBQ	170.30±0.02	153.20－189.20
TCBQ	195.20±0.01	186.50－204.20
2,6-DBBQ	170.70±0.02	155.00－188.00
2,5-DBBQ	94.90±0.01	88.60－101.70

HBQs	Molecular formula	logP	pKa	R	S	V	μ	E_HOMO	E_LUMO
2,6-DCBQ	C₆H₂Cl₂O₂	0.65	6.30	37.0	151.60	113.1	0.881	－10.977	－2.476
2,6-DBBQ	C₆H₂Br₂O₂	0.98	6.00	42.7	159.29	110.7	0.848	－10.464	－2.421
2,5-DBBQ	C₆H₂Br₂O₂	0.98	5.90	42.7	59.34	110.7	0.270	－10.462	－2.424
DCMBQ	C₇H₄Cl₂O₂	1.77	5.80	44.0	90.95	136.1	1.449	－9.989	－2.631
TCBQ	C₆HCl₃O₂	1.02	5.70	41.6	166.30	124.6	0.971	－10.598	－2.634

Parameter	MTS method		NRU method
Parameter	P value	R²	P value	R²
logP	0.441	0.208	0.303	0.339
pKa	0.311	0.330	0.109	0.629
R	0.245	0.410	0.188	0.490
V	0.575	0.116	0.333	0.306
S	0.465	0.189	0.292	0.351
μ	0.435	0.212	0.377	0.263
E_HOMO	0.357	0.282	0.274	0.373
E_LUMO	0.612	0.096	0.321	0.319

Parameter	MTS method		NRU method
Parameter	P value	R²	P value	R²
logP	0.441	0.208	0.303	0.339
pKa	0.311	0.330	0.109	0.629
R	0.245	0.410	0.188	0.490
V	0.575	0.116	0.333	0.306
S	0.465	0.189	0.292	0.351
μ	0.435	0.212	0.377	0.263
E_HOMO	0.357	0.282	0.274	0.373
E_LUMO	0.612	0.096	0.321	0.319