Quantitative relationship study for the
structure and biodegradability of substituted benzenes
Lu
Guanghua, Jiang Shiquan#, Zhao Yuanhui
(Department of Environmental Science, Northeast Normal University, ChangChun, 130024; #Design
Research Institute of Oil Chemical, Jilin Province, Changchun )
Received Sep. 29, 2000; Supported by the National Natural Science Foundation of China
(29877004).
Abstract The heat of formation(DHf), the
energy of the highest occupied molecular orbital(EHOMO), the molecular
weight(MW) of 52 compounds of substituted benzenes were calculated by the quantum
chemical method MOPAC6.0-AM1. The quantitative structure-biodegradability relationship
studies were performed with the maximum specific removal rates (QTOD). Through
multiple linear regression analysis, one equation was obtained as follows: QTOD=-0.614MW-9.704EHOMO-0.129DHf , n=52,
R2=0.860, SE=14.89, F=100.21, P=0.000. The QSBR equation was used to predict
biodegradability and biodegradation mechanism was discussed.
Keywords structure parameter, biodegradation correlation, QSBR model
1. INTRODUCTION
The number of synthetic organic chemicals
produced is large and increasing every year. The presence of many of these chemicals in
the ecosystem is a serious public health problem. Biodegradation is an important mechanism
for removing them from natural ecosystems. The high diversity of species and the metabolic
efficiency of microorganisms suggest that they play a main role in the ultimate
degradation of these chemicals [1]. Biodegradation can eliminate hazardous
chemicals by biotransforming them into innocuous forms, degrading by mineralization to
carbon dioxide and water.
Information regarding the extent and the rate of biodegradation of
organic chemicals is very important in evaluating relative persistence of the chemical in
the environment, which is important for regulating its manufacture and use. However,
gathering this information is labor intensive, time consuming, and expensive due to the
large number of these chemicals. Therefore, it is necessary to develop correlation and
predictive techniques in order to estimate their biodegradability [2]. At
present, the quantitative structure-biodegradability relationships (QSBRs) have been used
to predict the fate of organic chemicals and analyse biodegradation mechanism [3].
2. SOURCE OF DATA
The maximum specific removal rates of aromatic organic
compounds were obtained using an adapted mixed culture (activated sludge). This microbial
culture was cultivated semicontinuously on glucose and peptone and subsequently on the
compound under study at a mean biomass retention time (MBRT) of approximately 5 days. This
MBRT corresponds to biological treatment in the classical activated sludge wastewater
treatment process. During the biodegradation experiment the substance tested was the sole
source of organic carbon for the microbes of the culture. The initial concentration of the
compounds tested corresponded to the oxygen equivalent (TOD-theoretical oxygen demand) of
200 mg/l. The maximum specific removal rates are expressed as mg TOD of the compound
removed per gram of the initial solids of the microbial culture per hour (QTOD).
QTOD values of 52 compounds of substituted benzenes were cited from [4].
3.CALCULATION METHODS
MW, DHf
and EHOMO of the substituted benzenes were calculated by the quantum
chemical method MOPAC6.0-AM1 on energy-minimized structures. This method can automatically
optimize the bond length, the bond angle and the twist angle, and yield a lot of structure
information. Octanol-water partition coefficient (logP) was obtained from Biobyte
package. The parameter values of 52 compounds were listed in Table 1. All statistical
analyses were carried out using SPSS package.
4.RESULTS AND
DISCUSSION
Biodegradation of a chemical in the aquatic environment is predominantly through microbial
attack, and thus, we may presume, via enzymatic processes [3]. In
general, the factors determining the rate of biodegradation can be divided into two kinds.
(1) The uptake rates and transport rates (e.g., the uptake rates by microbial cells or
transport rates within the cell to the relevant enzymes). (2) The rates binding to the
active site of an enzyme, and/or by the rate at which they undergo enzymatic
transformation. In the absence of a specific uptake mechanism, organic compounds are
probably transported into bacterial cells by passive diffusion through the lipid membrane.
If the diffusion of chemicals in cell membrane and water belong to partition process, the
diffusion coefficient should be a direct proportion to logP. Therefore,
biodegradation rates should be related macroscopic hydrophobic parameters if diffusion and
uptake are rate-limiting step of biodegradation. The enzyme-catalyzed transformation of a
compound occurs by its binding to the site of the enzyme through the formation of hydrogen
or covalent bonds. The strength of this interaction is influenced by the electronic
structure of compound and the steric structure of compound coinciding with the active site
of enzyme. So, if bind to enzyme or transformation is rate-limiting step, biodegradation
rate of a compounds should be related to the factors influencing the binding or reacting
with enzyme (e.g. steric and/or electronic parameters).
Table 1.The experimental and
predicted biodegradability and structure parameters of 52 compounds of substituted
benzenes
Compound |
-EHOMO (eV) |
-DHf
(kJ/mol ) |
MW |
logP |
qTOD(mg/g·h) |
Exp.
|
Pre.* |
Acetanilid |
8.77 |
63.95 |
135.17 |
1.16 |
15 |
10.4 |
4-Aminoacetanilide |
8.42 |
66.63 |
150.18 |
0.08 |
11 |
-1.9 |
Aniline |
8.52 |
85.69 |
93.13 |
0.90 |
19 |
36.6 |
Benzoic acid |
10.08 |
284.07 |
122.12 |
1.87 |
86 |
60.0 |
Nitrobenzene |
10.56 |
105.71 |
123.11 |
1.85 |
14 |
40.5 |
2-Aminobenzoic acid |
8.78 |
295.74 |
137.14 |
1.21 |
27 |
39.2 |
2-Aminophenol |
8.36 |
97.35 |
109.13 |
0.62 |
21 |
26.7 |
2-Chloroaniline |
8.63 |
54.84 |
127.57 |
1.90 |
11 |
12.5
|
2-Chlorophenol |
9.19 |
111.48 |
128.56 |
2.15 |
14 |
24.6 |
2-Methylphenol |
9.00 |
123.06 |
108.14 |
1.95 |
54 |
36.8 |
2-Nitrobenzaldehyde |
10.77 |
3.89 |
151.12 |
1.74 |
14 |
12.2 |
2-Nitrobenzoic acid |
10.91 |
240.43 |
167.12 |
1.47 |
20 |
34.3 |
2-Nitrophenol |
9.95 |
66.00 |
139.11 |
1.79 |
14 |
19.7 |
2-Nitrotoluene |
10.17 |
79.09 |
137.14 |
2.30 |
33 |
24.7 |
Phthalic acid |
10.47 |
630.64 |
166.13 |
0.73 |
78 |
81.0 |
Pyrocatechol |
8.88 |
277.22 |
110.11 |
0.88 |
56 |
54.4 |
Salicylic acid |
9.47 |
469.66 |
138.12 |
2.26 |
95 |
67.6 |
2-Methylaniline |
8.44 |
-55.09 |
107.16 |
1.32 |
15 |
9.0 |
3-Aminobenzoic acid |
8.88 |
288.80 |
137.14 |
0.98 |
7 |
39.3 |
3-Aminophenol |
8.57 |
100.70 |
109.13 |
0.21 |
11 |
29.2 |
3-Chloroaniline |
8.73 |
-55.55 |
127.57 |
1.88 |
6 |
-0.8 |
3-Chlorophenol |
9.34 |
121.39 |
128.56 |
2.50 |
32 |
27.4 |
3-Methylphenol |
9.02 |
124.73 |
108.14 |
1.96 |
55 |
37.2 |
Isophthalic acid |
10.52 |
655.47 |
166.13 |
1.66 |
76 |
84.7 |
3-Nitrobenzaldehyde |
10.88 |
17.85 |
151.12 |
1.47 |
10 |
15.1 |
3-Nitrophenol |
9.95 |
66.00 |
139.11 |
2.00 |
18 |
19.7 |
3-Nitrotoluene |
10.20 |
-77.83 |
137.14 |
2.42 |
21 |
4.7 |
Resorcinol |
9.05 |
279.22 |
110.11 |
0.80 |
58 |
56.3 |
3-Methylaniline |
8.48 |
-54.42 |
107.16 |
1.40 |
30 |
9.5 |
4-Aminophenol |
8.27 |
95.89 |
109.13 |
0.04 |
17 |
25.6 |
4-Chloroaniline |
8.58 |
54.67 |
127.57 |
1.88 |
6 |
12.0 |
4-Chlorophenol |
9.13 |
122.47 |
128.56 |
2.39 |
40 |
25.5 |
4-Methylphenol |
8.88 |
124.52 |
108.14 |
1.94 |
55 |
35.9 |
Hydroquinone |
8.72 |
274.67 |
110.11 |
0.59 |
54 |
52.5 |
4-Hydroxybenzoic acid |
9.61 |
471.75 |
138.12 |
1.58 |
100 |
69.4 |
4-Nitrobenzaldehyde |
10.83 |
14.63 |
151.12 |
1.56 |
14 |
14.2 |
4-Nitrobenzoic acid |
10.90 |
259.16 |
167.12 |
1.89 |
20 |
36.6 |
4-Nitrophenol |
10.07 |
81.72 |
139.11 |
1.91 |
16 |
22.9 |
4-Nitrotoluene |
10.30 |
-72.11 |
137.14 |
2.37 |
33 |
6.4 |
4-Methylaniline |
8.35 |
-54.13 |
107.16 |
1.39 |
20 |
8.2 |
2-chloro-4-nitrophenol |
9.91 |
176.77 |
173.56 |
2.33 |
5 |
12.4 |
2,4-Diaminophenol |
9.23 |
138.52 |
124.00 |
-0.61 |
12 |
31.4 |
2,4-Dichlorophenol |
9.23 |
138.53 |
163.00 |
2.96 |
11 |
7.4 |
2,5-Dihydroxybenzoic acid |
9.38 |
565.06 |
154.15 |
1.62 |
80 |
70.6 |
2,3-Dimethylaniline |
9.59 |
8.04 |
121.18 |
1.81 |
13 |
19.7 |
2,5-Dimethylaniline |
9.67 |
14.61 |
121.18 |
1.86 |
4 |
21.3 |
3,4-Dimethylaniline |
9.80 |
11.42 |
121.18 |
1.86 |
3 |
22.3 |
2,3-Dimethylphenol |
9.61 |
173.32 |
122.17 |
2.42 |
40.6 |
35 |
2,4-Dimethylphenol |
8.77 |
152.15 |
122.17 |
2.47 |
29.7 |
28 |
2,6-Dimethylphenol |
8.89 |
150.19 |
122.17 |
2.47 |
9 |
30.7 |
3,5-Dimethylphenol |
8.97 |
156.21 |
122.17 |
2.47 |
11 |
32.2 |
2,4-Dinitrophenol |
10.81 |
36.91 |
184.11 |
1.79 |
6 |
-3.4 |
*Pre. is the predicted value calculated from
equation (2).
MW, DHf, EHOMO
and logP are selected as the structure parameters to establish the QSBRs. Through
multiple linear regression analyses, two models were obtained as Table 2.
Table 2. The results of
regression analyses (including constant or not)
Model 1 |
Coefficient |
SE |
t |
Sig. |
Model 2 |
Coefficient |
SE |
t |
Sig. |
constant |
8.549 |
25.049 |
0.341 |
0.734 |
constant |
0 |
|
|
|
MW |
-0.606 |
0.142 |
-0.462 |
0.000 |
MW |
-0.614 |
0.139 |
-4.431 |
0.000 |
EHOMO |
-8.688 |
3.550 |
-2.447 |
0.018 |
EHOMO |
-9.704 |
1.916 |
-5.064 |
0.000 |
DHf |
-0.129 |
0.013 |
-9.875 |
0.000 |
DHf |
-0.129 |
0.013 |
-9.960 |
0.000 |
Model 1 and 2 can be expressed as equation
(1) and (2) respectively as follows:
QTOD=8.549-0.606MW-8.688EHOMO-0.129DHf , n=52,
R2=0.677, SE=15.03, F=33.48, P=0.000
(1)
QTOD=-0.614MW-9.704EHOMO-0.129DHf , n=52,
R2=0.860, SE=14.89, F=100.21, P=0.000
(2)
There, R2 is the square of
correlation coefficient; SE is the standard error; F is the mean square radio; P is the
significant level; and n is the number of compounds.The results in Table 2 show that SE of
the constant item of model 1 is very high; through comparing equation (1) and (2), R2
and F value of (2) are much higher, and SE of (2) is lower. Therefore, (2) was used to
estimate the biodegradability (Table 1).
The biodegradability of 52 substituted benzenes can be divided into
three kinds according to their QTOD values in Table 1: non-readily
biodegradable (QTOD<30mg/g·h), biodegradable (30-60 mg/g·h), and
readily biodegradable (QTOD>60 mg/g·h). The correct prediction rate of
model 2 is up to 80%.
The obtained QSBR models show that the biodegradation of studied
compounds is related mainly to steric parameter MW, and electronic parameter DHf and EHOMO.
MW is the molecular weight, and it can reflect the size of a molecule. The smaller
the molecule, the smaller the steric resistance in microorganisms, the easier the molecule
penetrate into the cell through its cell membrane and arrive at the active site of an
enzyme [5]. Therefore, the smaller the MW, the easier the compound is
biodegraded. The heat of formation (DHf) of a compound is a measure of its stability, and
hence it would be expected to reflect the ability of the compound to biodegrade [6].
The lower the DHf value, the easier the chemical is degraded by
microbes. EHOMO is the energy of the highest occupied molecular orbital,
and it is related to ionization potential [7]. The higher the -EHOMO,
the higher the QTOD values, which shows that the enzyme-catalyzed reaction of
the compounds is related to electronic transfer.
Poor correlation between QTOD and logP (r is only
0.01) adding above analyses demonstrates that the rate-limiting step within the overall
biodegradation process is the rates of binding with the active site of an enzyme and
enzymatic reaction instead of the rates of uptake and transport.
REFERENCES
[1] Alexander M. Science, 1980, 211: 132-138.
[2] Sterier M P. Environ. Sci. Technol., 1980, 14-28-31.
[3] Dearden J C. SAR and QSAR in Environmental Research, 1996, 5: 17-26.
[4] Petter P, Sykora V. Biodegradation Prediction, Peijnenburg W J G M and Damborsky
J (ed.) Academic Press, 1996,17-26.
[5] Boething R S. Environ. Toxicol. Chem., 1986, 5: 797-806.
[6] Dearden J C, Cronin M T D. Biodegradation Prediction, Peijnenburg W J G M and
Damborsky J (ed.) Academic Press, 1996, 93-105.
[7] Liu C Q. Quantum Biology and Application, Beijing: Higher Education Press, 1990, 16.
|