http://www.chemistrymag.org/cji/2001/031003pe.htm

  Jan. 5, 2001  Vol.3 No.1 P.3 Copyright cij17logo.gif (917 bytes)


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.129
DHf , n=52, R2=0.677, SE=15.03, F=33.48, P=0.000              (1)
QTOD=-0.614MW-9.704EHOMO-0.129
DHf , 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.


[ Back ] [ Up ] [ Next ] Mirror Site in USA  Europe  China  CSTNet ChinaNet