Next Article in Journal
Antimicrobial Activity and Barrier Properties against UV Radiation of Alkaline and Enzymatically Treated Linen Woven Fabrics Coated with Inorganic Hybrid Material
Previous Article in Journal
Synthesis of [18F]F-γ-T-3, a Redox-Silent γ-Tocotrienol (γ-T-3) Vitamin E Analogue for Image-Based In Vivo Studies of Vitamin E Biodistribution and Dynamics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 23
10.3390/molecules25235697
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Quantum Chemical Investigation of Polychlorinated Dibenzodioxins, Dibenzofurans and Biphenyls: Relative Stability and Planarity Analysis

by
Sopanant Datta
and
Taweetham Limpanuparb
*
Science Division, Mahidol University International College, Mahidol University, Salaya, Phutthamonthon, Nakhon Pathom 73170, Thailand
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(23), 5697; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25235697
Submission received: 16 November 2020 / Revised: 1 December 2020 / Accepted: 1 December 2020 / Published: 3 December 2020
(This article belongs to the Section Computational and Theoretical Chemistry)
Browse Figures
Versions Notes

Abstract

:
All the possible polychlorinated aromatic compounds in the classes of dibenzodioxins (PCDDs), dibenzofurans (PCDFs), and biphenyls (PCBs) were studied by the quantum chemical methods of HF/6-311++G(d,p), B3LYP/6-311++G(d,p), and MP2/cc-pVTZ. The calculated stabilities and structures of these compounds were compared with the available data on their abundance and toxicity. Prediction models for trends in energy and planarity among these congeners were proposed. The results discussed here can help contribute to the understanding of the role of dioxin-like compounds (DLCs) in the environment.

Graphical Abstract

1. Introduction

The Stockholm Convention on Persistent Organic Pollutants (POPs) [1], initially signed in 2001, includes 12 classes of compounds, three of which can be regarded as dioxin-like compounds (DLCs). (see Figure 1). A subset of DLCs was assigned toxic equivalency factors (TEFs) by the World Health Organization (WHO) [2] to indicate their potential health impacts on humans and mammals. These compounds may accumulate and bioconcentrate through the food chain and, as a result, food products in some jurisdictions such as the European Union are strictly monitored for the presence of these compounds [3]. Quantum chemical methods, semi-empirical calculations, and other related computational techniques have previously been applied to polychlorinated dibenzodioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and polychlorinated biphenyls (PCBs) [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. The goals of these studies were, primarily, to predict their chemical properties and understand their mechanisms of action or degradation in biotic systems. A summary can be found in Supplementary Table S1.
This communication describes the semi-automated codes for generating chemical structures, to be used as inputs for quantum chemical calculations in a standard package and to extract information from outputs for further analysis. Three different representative quantum chemical methods were selected for comparison: Hartree-Fock (HF), Density Functional Theory (DFT), and second-order Møller–Plesset perturbation theory (MP2). The lists of 76 PCDDs, 136 PCDFs, and 210 PCBs (inclusive of three non-chlorinated base structures) along with their systematic numberings [37,38,39] are shown in Supplementary Table S2–S4, respectively. The use of numberings herein will be in accordance with those in the lists. (This is not to be confused with PubChem ID in the Supplementary Materials and elution order number which is instrument-specific.) As shown in Supplementary Table S5, the congeners can be categorized into homologue groups by their empirical formulas or degrees of substitution.
To the best of our knowledge, in comparison to the existing literature this is the most comprehensive theoretical study in terms of the number of compounds in the study and the level of theory employed for the investigation. The computational results are presented in terms of the relative stability and planarity of compounds. These calculated results are compared to the experimental results in the literature in terms of their isomer distribution and toxicity. A better understanding of the nature of these compounds can lead us to a better management of and more sustainable solutions for POPs in the future.

2. Results

2.1. Relative Stability

The association between the substitution patterns of each compound and its relative stability, with respect to its isomers, was assessed. The substitution pattern parameters include the degree of chlorination (using the number of H atoms, or the number of Cl atoms subtracted from 10 for PCBs and 8 for PCDDs and PCDFs) and the intra-ring and cross-ring interactions between two substituents. Intra-ring interactions are quantified by the number of pairs of Cl and non-hydrogen substituents (Cl, O, and C) at o-, m- and p- positions on each phenyl ring. Cross-ring interactions are defined as the number of pairs of Cl substitution over an O-bridge (CR-O) and a C-C bond (CR-C). The possible parameters for intra-ring and cross-ring interactions are unique for each class of compounds, as shown in Table 1. The total number of o-, m- and p- pairs were also calculated from the tally of substituent pairs mentioned earlier. Overall, as listed in the table, there are a total of 11, 15, and 11 potential parameters for PCDDs, PCDFs, and PCBs, respectively. The following simple linear model was proposed for a relative stability prediction based on HF, B3LYP ∆G at 300 K, and MP2 electronic energy S = x0 + 102x1 + 104x2 + 106x3 + 108x4 + … where S is the energy score and x0, x1, x2, … are possible parameters. Coefficients of 102n are used as some variables can take up to two-digit values (0 to 12 for example). The association between energy scores and the calculated energy values was assessed using Spearman’s correlation coefficient (ρ). Inspired by the knapsack problem, parameters that give the highest ρ values were added to the model one by one, until the increase in the ρ value became less than 0.0001.
Table 2 shows the resulting parameters, along with the corresponding ρ values. The results from all methodologies gave ρ > 0.99, and were mostly in agreement with each other. The major predictor for all compounds, in addition to the degree of substitution, are the intra-ring interactions. The predictor with the highest priority is the number of substituent pairs at o- positions, regardless of substituent type, followed by substitutions at the m- and p- positions, respectively. This supports the expectations of steric hindrance from the substituents (bulkiness: Cl > O ≈ C > H). Further, a cross-ring interaction over a C-C bond is an important parameter, especially for PCDFs, whereas a cross-ring interaction over an O-bridge is of less priority for PCDDs. With these results, the stability ranks of these compounds, regardless of methodology, can be sufficiently explained by their substituent pair interactions. For example, the energy scores S of PCDDs-13 and 14, when fitted according to the B3LYP results, can be calculated as:
SPCDD-13 = x0 + 102x1 + 104x2 + 106x3 + 108x4 + 1010x5
= 1 + 102(0) + 104(1) + 106(2) + 108(3) + 1010(5)
= 50,302,010,001
SPCDD-14 = x1 + 102x2 + 104x3+ 106x4+ 108x5 + 1010x6
= 2 + 102(1) + 104(1) + 106(1) + 108(3) + 1010(5)
= 50,301,010,102, respectively.
PCDD-13 ranks higher than PCDD-14 and is predicted to have a higher energy and less stability.
The weights used in this rank modelling are to provide a non-overlapping priority of consideration for each predictor variable. Qualitative interpretation then becomes obvious. Multiple linear regression for the prediction of energy values, following a similar knapsack approach, was also explored and is shown in Supplementary Table S6. In this case, the priorities of the predictors overlap and the already high ρ values can be increased further. Our computational model is able to predict the relative stabilities of all of these classes of compounds with a stronger correlation than the linear models proposed earlier [16,17,18,19], while not increasing the numbers of parameters.

2.2. Isomer Distribution

The energetic data in the previous subsection were used to predict the distribution of isomers within each homologue group. For this analysis, energy values were used to evaluate the distribution of isomers using the Boltzmann distribution at 300, 600, and 900 K [40]. The results were compared to available data on the abundance of these compounds in nature [37,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56]. A brief review of these data in the literature is shown in Supplementary Table S7. Our PCDDs and PCDFs results agree with a limited number of reports on experimental isomer distribution from incineration sources, while no clear trends were observed for PCBs (see Supplementary Table S8. The median ρ values calculated for each homologue group are 0.5833, 0.6555, and −0.1044 for PCDDs, PCDFs, and PCBs, respectively). For all three classes, the distribution of compounds found in nature usually differ due to the different accumulation and decomposition mechanisms in biotic and abiotic sources. For PCBs, it was argued that the distribution of isomers from the source is kinetically controlled rather than thermodynamically controlled [26].

2.3. Planarity

Planarity is widely claimed to be an important factor for these compounds to bind with the aryl hydrocarbon receptor (AhR) [57]. Therefore, geometric data were extracted and selected dihedral angles were calculated to represent the coplanarity of two phenyl rings. The designated carbon atoms for dihedral angle calculation are shown in Figure 1 and the list below.
PCDDPCDFPCB
Ca,Cc,Cd,Ch
Cb,Cd,Cc,Cg
Ca,Ce,Cf,Ch
Cb,Cf,Ce,Cg		Ca,Cd,Cf,Cg
Cb,Cc,Ce,Ch		Ci,Ck,Cl,Cp
Cj,Cl,Ck,Co
Ci,Cm,Cn,Cp
Ci,Cn,Cm,Co		Ci,Cl,Cn,Co
Cj,Ck,Cm,Cp		Cq,Cs,Ct,Cu
Cq,Cs,Ct,Cv
Cr,Cs,Ct,Cu
Cr,Cs,Ct,Cv
The lists for PCDDs and PCDFs are separated into two groups representing two orthogonal axes of rotation. Absolute values of the acute angle representations of the angles were used for the mean calculation of each group.
Our data show that most PCDDs are planar, and therefore the data are insufficient for further prediction.
For PCDFs, with a limited number of compounds exhibiting non-planarity, there appears to be a moderate trend (R values ranging from 0.61 to 0.65) between the substitution pattern and coplanarity. The significant parameters (p < 0.05) involved are the presence or absence of a substituent pair at positions 1 and 9 and at positions 2 and 8. The reasons for this finding may be unclear due to a limited number of data points (n = 13) for analysis and moderate correlation. The prediction models were derived from multiple linear regression and are shown in Table S9.
For PCBs, a dihedral angle prediction, based on chlorine substitution position, was also assessed using multiple linear regression. In this model, the dihedral angle A of a compound is predicted by the equation:
A = C0 + C1(x2 + x2′ + x6 + x6′) + C2(x2 + x6)(x2′ + x6′) + C3(x3 + x3′ + x5 + x5′) + C4(x4 + x4′)
where; C1, C2, C3, and C4 are the coefficients to be determined; and xi is 1 if position i is occupied by a Cl atom, and 0 if otherwise.
Table 3 shows the coefficients of significant parameters (p < 0.05) and the corresponding correlation coefficients for the predictions. The results from all methodologies show that the most important parameter is Cl substitution at any one of the four positions ortho to the C-C bond connecting two rings, as mentioned in earlier findings [57]. This is then followed by the same kind of substitution on both rings (product of x2 + x6 and x2′ + x6′). Meta substitutions according to all methods and para substitutions, according to MP2 results, are less important parameters, as C3 and C4 are relatively small.

2.4. Toxicity

When comparing our geometric data to the toxic equivalency factor (TEF) values of these compounds [2], as shown in Supplementary Table S10, the following points are observed.
  • The twelve known dioxin-like PCBs with TEF values are among the ones with relatively small dihedral angles or are closer to being planar.
  • According to the B3LYP and MP2 results, there are three non-coplanar PCDDs (PCDDs-46, 65, and 71) from both methods, and two additional (PCDDs-42 and 45) from B3LYP only. All of these are not substituted at all four lateral positions (2,3,7,8), which are the substitution positions of PCDDs known to be toxic [2]. One of these, PCDD-46, is substituted at positions exactly opposite to that of 2,3,7,8-tetrachlorodibenzodioxin (TCDD), the most toxic dioxin-like compound [2].

3. Discussion, Conclusions and Future Work

Overall, this study provides high-quality quantum chemical calculations for the energetic and geometric properties of PCDDs, PCDFs, and PCBs. The comprehensive data provided for all of these three related classes of compounds can be used for further investigations. From our computational results, we have proposed prediction models for the important characteristics of these compounds. The relative stability prediction extends beyond earlier prediction models [16,17,18,19] with a high accuracy and with the allowance of qualitative interpretations. Prediction models have also been proposed for planarity, an important structural parameter for the toxic activities of these compounds. Furthermore, we have made associations between theoretical calculations (relative stability and planarity) and experimental data in the literature (abundance and toxicity).
Our work can be used for further analyses on the persistence and activities of the toxic compounds in nature. In terms of methodology, semi-automated codes in the Supplementary Materials can be extended further for the construction and analysis of other related compounds. All three methods, HF, B3LYP, and MP2, yielded similar trends for the classes of compounds in our study. Our results suggest that a computationally expensive treatment for correlation energy may not be needed for this purpose and a low-level method such as HF, in most cases, can sufficiently provide insightful findings on the molecular properties of the studied classes of compounds. The further investigation of different variations of DFT methods may also be useful. For example, range-separated hybrid functionals such as ωB97XD can be used for comparison with the results presented here.

4. Materials and Methods

We extended our exhaustive combinatorial investigation approach [58,59,60] to all members of the three classes of compounds. The Becke three-parameter hybrid functional combined with the Lee-Yang-Parr correlation functional (B3LYP), the most popular variation of DFT, was selected for this work. In total, three methods, HF/6-311++G(d,p), B3LYP/6-311++G(d,p), and MP2/cc-pVTZ [61], were used on Q-Chem 5.1 (developer version) [62] for geometry optimization. The HF- and B3LYP-optimised geometries were confirmed to be minima by verifying the absence of imaginary frequencies. All the raw output files and incomplete attempts made with MP2/6-311++G(d,p) are provided in the Supplementary Materials. The process of the structure enumeration and analysis of molecular properties, particularly relative energy and planarity, were completed by semi-automated scripts written in Mathematica for this project. The Mathematica notebook can be found in the Supplementary Materials. These Mathematica codes can be used for the construction and analysis of other compounds and for teaching purposes in classes relating to molecular modelling and cheminformatics. Our computational results were compared to the relative abundance of the compounds in nature and their estimated toxicity in humans and animals. The analyses presented are derived from the energetic and geometric data extracted from the computational results.
For energetic information, the electronic energy of each compound was obtained from the results of all methodologies and, additionally, the enthalpy, entropy, and Gibbs energy values were obtained from the frequency jobs of HF and B3LYP. Due to the free rotation about C-C bonds connecting the two phenyl rings in PCBs, 78 potential atropisomers [38] were explored and only the structures with the lower energy values were used for further calculations.

Supplementary Materials

The following are available online, Table S1: Summary of computational data and analyses on PCDDs, PCDFs, and PCBs in the literature; Table S2: List of all PCDDs by substitution positions and numbering, proposed by Ballschmiter et al., revised by Dorofeeva et al.; Table S3: List of all PCDFs by substitution positions and numbering, proposed by Ballschmiter et al.; Table S4: List of all PCBs by substitution positions and numbering, first proposed by Ballschmiter et al., later revised by IUPAC; Table S5: List of congener groups of PCDDs, PCDFs and PCBs by degree of substitution; Table S6: Values for constant C0 and coefficients C1, C2, C3,…, parameters x1, x2, x3,… and correlation coefficients for the energy prediction model E = C0 + C1x1 + C2x2 + C3x3 + …, where E is the predicted HF, B3LYP ∆G at 300 K in Hartree and MP2 electronic energy in Hartree; Table S7: Summary of experimental data on abundance of PCDDs, PCDFs and PCBs in the literature; Table S8: Rank correlation coefficients (ρ values) for association between predicted distribution (using Boltzmann distribution) of PCDDs and PCDFs based on relative stability and observed distribution in fly ash samples from municipal waste incinerators (MWI); Table S9: Values for constant C0, coefficients C1, C2, and correlation coefficients for the PCDFs coplanarity prediction model; Table S10: List of toxic equivalency factor (TEF) values of 7 PCDDs, 10 PCDFs and 12 PCBs adapted from Van den Berg et al.

Author Contributions

Conceptualization, S.D. and T.L.; methodology, S.D. and T.L.; software, T.L.; validation, S.D. and T.L.; formal analysis, S.D.; investigation, T.L.; resources, T.L.; data curation, S.D. and T.L.; writing—original draft preparation, S.D. and T.L.; writing—review and editing, S.D. and T.L.; visualization, S.D.; supervision, T.L.; project administration, T.L.; funding acquisition, T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Council of Thailand (NRCT) and Mahidol University.

Acknowledgments

We thank Wanutcha Lorpaiboon for constructive comments on the manuscript and are grateful to Wirote Chansri (UAE-IDEA Advance Analytical Company Limited) for an opportunity to learn about dioxin analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. United Nations. Treaties and International Agreements Registered or Filed and Recorded with the Secretariat of the United Nations; United Nations: New York, NY, USA, 2006. [Google Scholar]
  2. Van den Berg, M.; Birnbaum, L.S.; Denison, M.; De Vito, M.; Farland, W.; Feeley, M.; Fiedler, H.; Hakansson, H.; Hanberg, A.; Haws, L.; et al. The 2005 World Health Organization reevaluation of human and mammalian toxic equivalency factors for dioxins and dioxin-like compounds. Toxicol. Sci. 2006, 93, 223–241. [Google Scholar] [CrossRef] [Green Version]
  3. European Commission. Commission Regulation (EC) No 1881/2006 of 19 December 2006 setting maximum levels for certain contaminants in foodstuffs. Off. J. Eur. Union 2006, 364, 5–24. [Google Scholar]
  4. Arulmozhiraja, S.; Fujii, T.; Tokiwa, H. Electron affinity for the most toxic 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD):  A density functional theory study. J. Phys. Chem. A 2000, 104, 7068–7072. [Google Scholar] [CrossRef]
  5. Dorofeeva, O.V.; Moiseeva, N.F.; Yungman, V.S.; Novikov, V.P. Ideal gas thermodynamic properties of biphenyl. Thermochim. Acta 2001, 374, 7–11. [Google Scholar] [CrossRef]
  6. Iino, F.; Tsuchiya, K.; Imagawa, T.; Gullett, B.K. An isomer prediction model for PCNs, PCDD/Fs, and PCBs from municipal waste incinerators. Environ. Sci. Technol. 2001, 35, 3175–3181. [Google Scholar] [CrossRef]
  7. Arulmozhiraja, S.; Fujii, T.; Morita, M. Density functional theory studies on radical ions of selected polychlorinated biphenyls. J. Phys. Chem. A 2002, 106, 10590–10595. [Google Scholar] [CrossRef]
  8. Chana, A.; Concejero, M.A.; de Frutos, M.; González, M.J.; Herradón, B. Computational studies on biphenyl derivatives. Analysis of the conformational mobility, molecular electrostatic potential, and dipole moment of chlorinated biphenyl:  Searching for the rationalization of the selective toxicity of polychlorinated biphenyls (PCBs). Chem. Res. Toxicol. 2002, 15, 1514–1526. [Google Scholar] [CrossRef]
  9. León, L.A.; Notario, R.; Quijano, J.; Sánchez, C. Structures and enthalpies of formation in the gas phase of the most toxic polychlorinated dibenzo-p-dioxins. A DFT study. J. Phys. Chem. A 2002, 106, 6618–6627. [Google Scholar] [CrossRef]
  10. Dorofeeva, O.V.; Yungman, V.S. Enthalpies of formation of dibenzo-p-dioxin and polychlorinated dibenzo-p-dioxins calculated by density functional theory. J. Phys. Chem. A 2003, 107, 2848–2854. [Google Scholar] [CrossRef]
  11. Lee, J.-E.; Choi, W. DFT calculation on the electron affinity of polychlorinated dibenzo-p-dioxins. Bull. Korean Chem. Soc. 2003, 24, 792–796. [Google Scholar] [CrossRef] [Green Version]
  12. Lee, J.E.; Choi, W.; Mhin, B.J. DFT calculation on the thermodynamic properties of polychlorinated dibenzo-p-dioxins:  intramolecular Cl−Cl repulsion effects and their thermochemical implications. J. Phys. Chem. A 2003, 107, 2693–2699. [Google Scholar] [CrossRef]
  13. Arulmozhiraja, S.; Morita, M. Electron affinities and reductive dechlorination of toxic polychlorinated dibenzofurans: A density functional theory study. J. Phys. Chem. A 2004, 108, 3499–3508. [Google Scholar] [CrossRef]
  14. Dorofeeva, O.V.; Moiseeva, N.F.; Yungman, V.S. Thermodynamic properties of polychlorinated biphenyls in the gas phase. J. Phys. Chem. A 2004, 108, 8324–8332. [Google Scholar] [CrossRef]
  15. Ryu, J.-Y.; Mulholland, J.A.; Dunn, J.E.; Iino, F.; Gullett, B.K. Potential role of chlorination pathways in PCDD/F formation in a municipal waste incinerator. Environ. Sci. Technol. 2004, 38, 5112–5119. [Google Scholar] [CrossRef]
  16. Wang, Z.-Y.; Zhai, Z.-C.; Wang, L.-S.; Chen, J.-L.; Kikuchi, O.; Watanabe, T. Prediction of gas phase thermodynamic function of polychlorinated dibenzo-p-dioxins using DFT. J. Mol. Struc. Theochem 2004, 672, 97–104. [Google Scholar] [CrossRef]
  17. Wang, Z.-Y.; Zhai, Z.-C.; Wang, L.-S. Prediction of gas phase thermodynamic properties of polychlorinated dibenzo-furans by DFT. J. Mol. Struc. Theochem 2005, 725, 55–62. [Google Scholar] [CrossRef]
  18. Wang, Z.Y.; Han, X.Y.; Zhai, Z.C.; Wang, L.S. Study on the thermodynamic property and relative stability of a series of polychlorinated biphenyls by density functional theory. Acta Chim. Sinica 2005, 63, 964–972. [Google Scholar]
  19. Zhai, Z.-C.; Wang, Z.-Y.; Chen, X.-H.; Wang, L.-S. DFT calculation on 204 polychlorinated biphenyls: Their thermodynamic function and implication of Cl substitute position. J. Mol. Struc. Theochem 2005, 714, 123–131. [Google Scholar] [CrossRef]
  20. Dorofeeva, O.V.; Vishnevskiy, Y.V.; Moiseeva, N.F. Assessment of Gaussian-3X theory for chlorinated organic molecules. Enthalpies of formation of chlorobenzenes and predictions for polychlorinated aromatic compounds. Struct. Chem. 2006, 17, 383–392. [Google Scholar] [CrossRef]
  21. Thompson, D.; Ewan, B.C.R. A group additivity algorithm for polychlorinated dibenzofurans derived from selected DFT analyses. J. Phys. Chem. A 2007, 111, 5043–5047. [Google Scholar] [CrossRef]
  22. Wang, L.; Heard, D.E.; Pilling, M.J.; Seakins, P. A Gaussian-3X prediction on the enthalpies of formation of chlorinated phenols and dibenzo-p-dioxins. J. Phys. Chem. A 2008, 112, 1832–1840. [Google Scholar] [CrossRef] [PubMed]
  23. Jansson, S.; Antti, H.; Marklund, S.; Tysklind, M. Multivariate relationships between molecular descriptors and isomer distribution patterns of PCDD/Fs formed during MSW combustion. Environ. Sci. Technol. 2009, 43, 7032–7038. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, L.; He, Y.-L. The enthalpies of formation for polychlorinated dibenzofurans with use of G3XMP2 model chemistry and density functional theory. J. Phys. Chem. A 2009, 113, 238–245. [Google Scholar] [CrossRef] [PubMed]
  25. Ewan, B.C.R.; Thompson, D. Predicted specific heat behaviour of the polychlorinated dibenzofuran family from DFT analysis. Thermochim. Acta 2013, 558, 46–52. [Google Scholar] [CrossRef]
  26. Jansson, S.; Grabic, R. Multivariate relationships between molecular descriptors and isomer distribution patterns of PCBs formed during household waste incineration. Environ. Sci. Pollut. Res. 2014, 21, 3082–3090. [Google Scholar] [CrossRef]
  27. Zhou, Q.; Su, X.; Yong, Y.; Ju, W.; Fu, Z.; Li, X. Adsorption of 2, 3, 7, 8-tetrachlorodibenzao-p-dioxin (TCDD) on graphane decorated with Ni and Cu: A DFT study. Vacuum 2018, 149, 53–59. [Google Scholar] [CrossRef]
  28. Abronin, I.A.; Volkova, L.V. On the relationship between the energy characteristics of the isodesmic reactions of polychlorinated dioxins and their toxicity. Russ. Chem. Bull. 2019, 68, 867–869. [Google Scholar] [CrossRef]
  29. Bai, N.; Wang, W.; Zhao, Y.; Feng, W.; Li, P. Theoretical insights into the reaction mechanism between 2,3,7,8-tetrachlorodibenzofuran and hydrogen peroxide: A DFT study. ACS Omega 2019, 4, 358–367. [Google Scholar] [CrossRef]
  30. Hou, S.; Altarawneh, M.; Kennedy, E.M.; Mackie, J.C.; Weber, R.; Dlugogorski, B.Z. Formation of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/F) from oxidation of 4,4′-dichlorobiphenyl (4,4′-DCB). Proc. Combust. Inst. 2019, 37, 1075–1082. [Google Scholar] [CrossRef]
  31. Behjatmanesh-Ardakani, R.; Heydari, A. Molecular and dissociative adsorption of tetrachlorodibenzodioxin on M-doped graphenes (M = B, Al, N, P): Pure DFT and DFT + VdW calculations. J. Mol. Model. 2020, 26, 164. [Google Scholar] [CrossRef]
  32. Diao, J.; Li, Y.; Shi, S.; Sun, Y.; Sun, Y. QSAR models for predicting toxicity of polychlorinated dibenzo-p-dioxins and dibenzofurans using quantum chemical descriptors. Bull. Environ. Contam. Toxicol. 2010, 85, 109–115. [Google Scholar] [CrossRef]
  33. Larsson, M.; Kumar Mishra, B.; Tysklind, M.; Linusson, A.; Andersson, P.L. On the use of electronic descriptors for QSAR modelling of PCDDs, PCDFs and dioxin-like PCBs. SAR QSAR Environ. Res. 2013, 24, 461–479. [Google Scholar] [CrossRef]
  34. Hirano, M.; Hwang, J.-H.; Park, H.-J.; Bak, S.-M.; Iwata, H.; Kim, E.-Y. In silico analysis of the interaction of avian aryl hydrocarbon receptors and dioxins to decipher isoform-, ligand-, and species-specific activations. Environ. Sci. Technol. 2015, 49, 3795–3804. [Google Scholar] [CrossRef]
  35. Larsson, M.; Fraccalvieri, D.; Andersson, C.D.; Bonati, L.; Linusson, A.; Andersson, P.L. Identification of potential aryl hydrocarbon receptor ligands by virtual screening of industrial chemicals. Environ. Sci. Pollut. Res. 2018, 25, 2436–2449. [Google Scholar] [CrossRef] [Green Version]
  36. Zeinali, N.; Oluwoye, I.; Altarawneh, M.; Dlugogorski, B.Z. Destruction of dioxin and furan pollutants via electrophilic attack of singlet oxygen. Ecotoxicol. Environ. Saf. 2019, 184, 109605. [Google Scholar] [CrossRef]
  37. Ballschmiter, K.; Buchert, H.; Niemczyk, R.; Munder, A.; Swerev, M. Automobile exhausts versus municipal-waste incineration as sources of the polychloro-dibenzodioxins (PCDD) and -furans (PCDF) found in the environment. Chemosphere 1986, 15, 901–915. [Google Scholar] [CrossRef]
  38. Agudo, A.; Aronson, K.J.; Bonefeld-Jorgensen, E.C.; Cocco, P.; Cogliano, V.; Cravedi, J.-P.; Esch, H.; Fiedler, H.; Glauert, H.P.; Guo, Y.-L.L. Polychlorinated Biphenyls and Polybrominated Biphenyls; International Agency for Research on Cancer: Lyon, France, 2016. [Google Scholar]
  39. Ballschmiter, K.; Zell, M. Analysis of polychlorinated biphenyls (PCB) by glass capillary gas chromatography. Fresenius Z. Anal. Chem. 1980, 302, 20–31. [Google Scholar] [CrossRef]
  40. Dickson, L.C.; Lenoir, D.; Hutzinger, O. Quantitative comparison of de novo and precursor formation of polychlorinated dibenzo-p-dioxins under simulated municipal solid waste incinerator postcombustion conditions. Environ. Sci. Technol. 1992, 26, 1822–1828. [Google Scholar] [CrossRef]
  41. Tiernan, T.O.; Taylor, M.L.; Garrett, J.H.; VanNess, G.F.; Solch, J.G.; Deis, D.A.; Wagel, D.J. Chlorodibenzodioxins, chlorodibenzofurans and related compounds in the effluents from combustion processes. Chemosphere 1983, 12, 595–606. [Google Scholar] [CrossRef]
  42. Schwartz, T.R.; Campbell, R.D.; Stalling, D.L.; Little, R.L.; Petty, J.D.; Hogan, J.W.; Kaiser, E.M. Laboratory data base for isomer-specific determination of polychlorinated biphenyls. Anal. Chem. 1984, 56, 1303–1308. [Google Scholar] [CrossRef]
  43. Ballschmiter, K.; Swerev, M. Reaction pathways for the formation of polychlorodibenzodioxins (PCDD) and -furans (PCDF) in combustion processes I. Fresenius Z. Anal. Chem. 1987, 328, 125–127. [Google Scholar] [CrossRef]
  44. Yasuhara, A.; Ito, H.; Morita, M. Isomer-Specific Determination of polychlorinated dibenzo-p-dioxins and dibenzofurans in incinerator-related environmental samples. Environ. Sci. Technol. 1987, 21, 971–979. [Google Scholar] [CrossRef]
  45. Faroon, O.M.; Samuel Keith, L.; Smith-Simon, C.; De Rosa, C.T. Polychlorinated Biphenyls: Human Health Aspects; World Health Organization: Geneva, Switzerland, 2003. [Google Scholar]
  46. Hanari, N.; Horii, Y.; Okazawa, T.; Falandysz, J.; Bochentin, I.; Orlikowska, A.; Puzyn, T.; Wyrzykowska, B.; Yamashita, N. Dioxin-like compounds in pine needles around Tokyo Bay, Japan in 1999. J. Environ. Monit. 2004, 6, 305–312. [Google Scholar] [CrossRef]
  47. Kim, K.S.; Hirai, Y.; Kato, M.; Urano, K.; Masunaga, S. Detailed PCB congener patterns in incinerator flue gas and commercial PCB formulations (Kanechlor). Chemosphere 2004, 55, 539–553. [Google Scholar] [CrossRef]
  48. Haws, L.C.; Su, S.H.; Harris, M.; DeVito, M.J.; Walker, N.J.; Farland, W.H.; Finley, B.; Birnbaum, L.S. Development of a refined database of mammalian relative potency estimates for dioxin-like compounds. Toxicol. Sci. 2006, 89, 4–30. [Google Scholar] [CrossRef] [Green Version]
  49. Ishikawa, Y.; Noma, Y.; Yamamoto, T.; Mori, Y.; Sakai, S.-i. PCB decomposition and formation in thermal treatment plant equipment. Chemosphere 2007, 67, 1383–1393. [Google Scholar] [CrossRef]
  50. Srogi, K. Levels and congener distributions of PCDDs, PCDFs and dioxin-like PCBs in environmental and human samples: A review. Environ. Chem. Lett. 2008, 6, 1–28. [Google Scholar] [CrossRef] [Green Version]
  51. Sundqvist, K.L.; Tysklind, M.; Geladi, P.; Hopke, P.K.; Wiberg, K. PCDD/F source apportionment in the Baltic sea using positive matrix factorization. Environ. Sci. Technol. 2010, 44, 1690–1697. [Google Scholar] [CrossRef]
  52. Jansson, S.; Lundin, L.; Grabic, R. Characterisation and fingerprinting of PCBs in flue gas and ash from waste incineration and in technical mixtures. Chemosphere 2011, 85, 509–515. [Google Scholar] [CrossRef]
  53. Barahona-Urbina, C.; Nunez-Gonzalez, S.; Gomez-Jeria, J.S. Model-based quantum-chemical study of the uptake of some polychlorinated pollutant compounds by zucchini subspecies. J. Chil. Chem. Soc. 2012, 57, 1497–1503. [Google Scholar] [CrossRef] [Green Version]
  54. Hinwood, A.L.; Callan, A.C.; Heyworth, J.; Rogic, D.; de Araujo, J.; Crough, R.; Mamahit, G.; Piro, N.; Yates, A.; Stevenson, G.; et al. Polychlorinated biphenyl (PCB) and dioxin concentrations in residential dust of pregnant women. Environ. Sci. Process. Impacts 2014, 16, 2758–2763. [Google Scholar] [CrossRef] [Green Version]
  55. Thompson, D.; Ooi, T.C.; Anderson, D.R.; Fisher, R.; Ewan, B.C.R. The polychlorinated dibenzofuran fingerprint of iron ore sinter plant: Its persistence with suppressant and alternative fuel addition. Chemosphere 2016, 154, 138–147. [Google Scholar] [CrossRef]
  56. Huang, R.; Wang, P.; Zhang, J.; Chen, S.; Zhu, P.; Huo, W.; Jiang, Y.; Chen, Z.; Peng, J. The human body burden of polychlorinated dibenzo-p-dioxins/furans (PCDD/Fs) and dioxin-like polychlorinated biphenyls (DL-PCBs) in residents’ human milk from Guangdong Province, China. Toxicol. Res. 2019, 8, 552–559. [Google Scholar] [CrossRef]
  57. The Eightieth Meeting of the Joint FAO/WHO Expert Committee on Food Additives (JECFA). Safety Evaluation of Certain Food Additives and Contaminants, Supplement 1: Non-Dioxin-Like Polychlorinated Biphenyls; World Health Organization: Geneva, Switzerland, 2016. [Google Scholar]
  58. Datta, S.; Limpanuparb, T. Geometric and energetic data from quantum chemical calculations of halobenzenes and xylenes. Data Brief 2020, 30, 105386. [Google Scholar] [CrossRef]
  59. Limpanuparb, T.; Datta, S.; Chinsukserm, K.; Teeraniramitr, P. In silico geometric and energetic data of all possible simple rotamers made of non-metal elements. Data Brief 2020, 30, 105442. [Google Scholar] [CrossRef]
  60. Chinsukserm, K.; Lorpaiboon, W.; Teeraniramitr, P.; Limpanuparb, T. Geometric and energetic data from ab initio calculations of haloethene, haloimine, halomethylenephosphine, haloiminophosphine, halodiazene, halodiphosphene and halocyclopropane. Data Brief 2019, 27, 104738. [Google Scholar] [CrossRef]
  61. Riley, K.E.; Hobza, P. Assessment of the MP2 method, along with several basis sets, for the computation of interaction energies of biologically relevant hydrogen bonded and dispersion bound complexes. J. Phys. Chem. A 2007, 111, 8257–8263. [Google Scholar] [CrossRef]
  62. Shao, Y.; Gan, Z.; Epifanovsky, E.; Gilbert, A.T.; Wormit, M.; Kussmann, J.; Lange, A.W.; Behn, A.; Deng, J.; Feng, X. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol. Phys. 2015, 113, 184–215. [Google Scholar] [CrossRef] [Green Version]
Molecules 25 05697 g001 550
Figure 1. Substitution position numberings (in red) and selected carbon atoms for dihedral angle calculation (in blue) for PCDDs, PCDFs, and PCBs.
Figure 1. Substitution position numberings (in red) and selected carbon atoms for dihedral angle calculation (in blue) for PCDDs, PCDFs, and PCBs.
Molecules 25 05697 g001
Table
Table 1. Possible numbers of substituent pairs for each compound.
Table 1. Possible numbers of substituent pairs for each compound.
Substitution PositionsSubstituent PairPCDDPCDFPCB
Intra-Ring Interactions
o-Cl-Cl0–6 10–60–8
Cl-O0–40–2N/A
Cl-CN/A0–20–4
Total0–100–100–12
m-Cl-Cl0–40–40–8
Cl-O0–80–4N/A
Cl-CN/A0–40–4
Total0–120–120–12
p-Cl-Cl0–20–20–4
Cl-O0–40–2N/A
Cl-CN/A0–20–2
Total0–60–60–6
Cross-ring interactions
Over O-bridgeCl-Cl0–2 20–1 3N/A
Over C-C bondN/A0–1 30–4 4
1 Example: A fully substituted PCDD (PCDD-75) has six o-Cl-Cl pairs at positions (1,2), (2,3), (3,4), (6,7), (7,8), and (8,9). Possible pairs are substitutions at positions: 2 (1,9) and (4,6); 3 (4,6) over O-bridge and (1,9) over C-C bond; 4 (2,2′), (2,6′), (2′,6), and (6,6′) due to the free rotation about the C-C bond connecting two rings. N/A indicates a case whereby such a substituent pair is not possible.
Table
Table 2. Parameters (in order of x0, x1, x2, …) for relative stability prediction and the corresponding correlation coefficients (ρ).
Table 2. Parameters (in order of x0, x1, x2, …) for relative stability prediction and the corresponding correlation coefficients (ρ).
PCDDHFB3LYPMP2
ParametersCR-OCR-OCR-O
p-Cl-Clp-Cl-Clp-Cl-Cl
m-Cl-Clm-Cl-Clm-Cl-Cl
o-Cl-Clo-Cl-Clo-Cl-Cl
o-totalo-totalo-total
total Htotal Htotal H
ρ0.99980.99970.9996
PCDFHFB3LYPMP2
Parametersp-Cl-Om-Cl-Clm-Cl-Cl
m-Cl-Clp-Cl-Op-Cl-Cl
o-Cl-Oo-Cl-Op-Cl-O
o-Cl-Clo-Cl-Clo-Cl-O
CR-CCR-Co-Cl-Cl
total Ho-totalCR-C
total Htotal H
ρ0.99970.99620.9978
PCBHFB3LYPMP2
Parametersp-totalm-totalp-total
m-totalo-Cl-Clm-total
o-Cl-Clo-totalo-Cl-Cl
o-totaltotal Htotal H
total H
ρ0.99930.99900.9978
Table
Table 3. Constants, coefficients, and Pearson’s (R) and Spearman’s (ρ) correlation coefficients for the PCB planarity prediction model.
Table 3. Constants, coefficients, and Pearson’s (R) and Spearman’s (ρ) correlation coefficients for the PCB planarity prediction model.
CoefficientsRρ
C0C1C2C3C4
HF50.85820.325−10.7510.724N/A0.93540.8735
B3LYP37.92624.996−11.4190.529N/A0.95520.9599
MP233.30824.202−9.4990.668−0.5290.98000.9560
N/A indicates a case whereby the coefficient is not required and can be regarded as zero.
Sample Availability: Samples of the compounds are not available from the authors.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Datta, S.; Limpanuparb, T. Quantum Chemical Investigation of Polychlorinated Dibenzodioxins, Dibenzofurans and Biphenyls: Relative Stability and Planarity Analysis. Molecules 2020, 25, 5697. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25235697

AMA Style

Datta S, Limpanuparb T. Quantum Chemical Investigation of Polychlorinated Dibenzodioxins, Dibenzofurans and Biphenyls: Relative Stability and Planarity Analysis. Molecules. 2020; 25(23):5697. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25235697

Chicago/Turabian Style

Datta, Sopanant, and Taweetham Limpanuparb. 2020. "Quantum Chemical Investigation of Polychlorinated Dibenzodioxins, Dibenzofurans and Biphenyls: Relative Stability and Planarity Analysis" Molecules 25, no. 23: 5697. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25235697

Article Metrics

Back to TopTop