Molecular orbital studies on brominated diphenyl ethers. Part I--conformational properties

Polybrominated diphenyl ethers (PBDEs) are widely used as additive flame retardants and quantities in the environment are on the rise. Because they are structurally related to polychlorinated biphenyls and also to thyroid hormones, there is serious concern that PBDEs may pose a danger to human health. Knowledge of their conformational properties is key to assessing their environmental fate and risk. The conformational properties of PBDEs were investigated by quantum chemical methods including semiempirical self-consistent field molecular orbital (SCF-MO), ab initio SCF-MO and density functional theory (DFT). Conformational analyses of model congeners 2,2',4,6'-tetrabromodiphenyl ether and 2,3,4,4',5,6-hexabromodiphenyl ether, based on energy maps calculated by semiempirical AM1 method, may indicate that all PBDE congeners except those with the tetra-ortho-bromination are conformationally flexible (or soft) due to low energy barriers for interconversion of stable conformers. The results of the conformational analyses are in conformity with recently published X-ray crystallographic data. For comparison with the results of the semiempirical method, higher level ab initio and DFT models were applied as well. The optimized geometries all lie well inside low energy regions on the maps and thus also ascertain the semiempirical calculations. According to computed geometric parameters and net atomic charges, the model B3LYP/3-21G* seemed to give better results than B3LYP/6-31G* and HF/6-31G*.     

Crystal structure and density functional theory studies of toxic quinone metabolites of polychlorinated biphenyls

Lower chlorinated polychlorinated biphenyls (PCBs) are readily metabolized via hydroxylated metabolites to reactive PCB quinones. Although these PCB metabolites elicit biochemical changes by mechanisms involving cellular target molecules, such as the aryl hydrocarbon receptor, and toxicity by interacting with enzymes like topoisomerases, only few PCB quinones have been synthesized and their conformational properties investigated. Similar to the parent compounds, knowledge of the three-dimensional structure of PCB quinones may therefore be important to assess their fate and risk. To address this gap in our knowledge, the gas phase molecular structure of a series of PCB quinones was predicted using HF/3-21G, B3LYP/6-31G?? and UB3LYP/6-311G?? calculations and compared to the respective solid state structure. All three methods overestimated the Cl-C bond length, but otherwise provided a reasonable approximation of the solid state bond angles and bond lengths. Overall, the UB3LYP/6-311G?? level of theory yielded the best approximation of the molecular structure of PCB quinones in the solid state. Chlorine addition at the ortho position of both rings was found to increase the dihedral angle of the resulting quinone compound, which may have important implications for their interaction with cellular targets and, thus, their toxicity.     

A new potential toxaphene congener: synthesis, GC/EI-MS study, crystal structure, NMR analysis, and ab initio calculations of 3-endo,5-endo-dichloro-7,7-bis-chloromethyl-4-dichloromethyl-tricyclo[2.2.1.0(2,6)]heptane

A new potential toxaphene congener 3-endo,5-endo-dichloro-7,7-bis-chloromethyl-4-dichloromethyl-tricyclo[2.2.1.0(2,6)]heptane 2 has been isolated from reaction mixture obtained by the chlorination of 2-exo, 10,10-trichlorobornane 1. The X-ray structural analysis of 2 revealed an unusual tricyclic structure, where the two chlorine atoms occupying endo-positions are in close spatial proximity with each other and near to the neighbouring CHCl2 group. Further, it revealed that the symmetry of the molecule is distorted. The 1H and 13C NMR spectra of 2 have been assigned by means of 1H, 1H double-quantum filtered correlation spectroscopy (DQF COSY), PFG 1H, 13C HMQC (pulsed field gradient heteronuclear multiple-quantum coherence), 1H, 13C heteronuclear multiple bond correlation (HMBC) experiments, and computer aided 1H NMR spectral analysis. The asymmetry of 2 is also discernible on the 1H NMR parameters. In addition, gas chromatographic (GC) properties and electron impact (EI) mass spectrum of 2 has been studied. Ab initio Hartree-Fock (HF) method with the basis set 6-31G(d) has been used for the optimization of the equilibrium geometry and calculation of total energy for 2. The optimized geometry is in good agreement with the crystal structure. According to the rotation energy profile calculated at the HF/6-31G(d) level, rotation of the chloromethyl and dichloromethyl groups are highly unlikely at the room temperature.     

Theoretical studies of the conformations and 19F NMR spectra of linear and a branched perfluorooctanesulfonamide (PFOSAmide)

Technical perfluorooctanesulfonate (PFOS) and its derivatives, such as perfluorooctanesulfonamide (PFOSA), are not clean compounds but, instead, complex mixtures of linear and branched isomers, and other compounds including sulfonate homologues. Questions have been raised as to whether the linear and the branched isomers behave differently in the environment. However, little is known about the physical properties or the finer details of the structures of the individual branched isomers. This study sought an effective computational method to model the preferred conformations of PFOS derivatives, and the energy differences between them and to determine if these results can be used to explain the temperature dependence of their NMR spectra. Good predictions of the 19F chemical shifts were obtained for some PFOSA-type molecules with a computational approach [B3LYP-GIAO/6-31++G(d,p)//B3LYP/6-31G(d,p)] that is relatively inexpensive. Large 5JFF couplings found in one of the branched isomers could be rationalized on the basis of the relevant F-F distances in the optimized structure. At low temperatures, the splitting observed in the NMR spectrum at C-1 for these sulfonamides can be explained by the existence of the two conformers predicted by the computations.     

Quantitative structure-activity relationships of nitroaromatics toxicity to the algae (Scenedesmus obliguus)

The DFT-B3LYP method, with the basis set 6-311G( * *), was employed to calculate the molecular geometries and electronic structures of 25 nitroaromatics. The acute toxicity (-lgEC(50)) of these compounds to the algae (Scenedesmus obliguus) along with hydrophobicity described by logK(OW), and two quantum chemical parameters-energy of the lowest unoccupied molecular orbital, E(LUMO), and the charge of the nitro group, [ForQ(NO2), were used to establish the quantitative structure-activity relationships (QSARs). For 18 mononitro derivatives, the hydrophobicity parameter logK(OW) could interpret the toxic mechanism successfully. Dinitro aromatic compounds were susceptible to be reduced to aniline for their electrophilic nature. Their toxicity was controlled mainly by electronic factors instead of hydrophobicity. The electronic parameters, E(LUMO) and Q(NO2), were used to yield the following model: -lg EC(50) = 3.746 - 25.053 E(LUMO) + 6.481 Q(NO2) (n=22, R=0.926, SE=0.206, F=56.854, P<0.001). The predicted toxic values using the above equation are in good agreement with the experimental values.     

NMR and molecular modeling in environmental chemistry: prediction of 13C chemical shifts in selected C10-chloroterpenes employing DFT/GIAO theory

Accurate predictions of 13C NMR chemical shifts (standard error approximately 1.7 ppm) are achieved for a subset of chlorinated bornanes by empirical scaling of shifts from GIAO calculations with geometries obtained from HF/6-31G* calculations. The optimized molecular geometries were compared with X-ray structures for three of the toxaphene components most frequently detected in environmental samples (Parlar nos. 26, 50 and 62), and the concordance between the experimental and calculated values was found to be satisfactory. Taken overall, the results indicate that theoretical methods hold great promise for rationalizing 13C NMR chemical shifts in organohalogen compounds. However, it appeared that the DFT/GIAO shifts need to be empirically scaled to achieve good numerical agreement with experimental shifts in chlorinated bornanes. Obviously, there is a need to develop new computational methods to describe the large deshielding effects of chlorine atoms properly.     

A noval method of predicting soil organic carbon content normalized adsorption coefficients (Koc) of polybrominated biphenyls (PBBs)

By a soil column liquid chromatograph (SCLC) method, the soil organic carbon content normalized adsorption coefficients (K(oc)) of six polybrominated biphenyls (PBB15, PBB26, PBB31, PBB49, PBB103 and PBB153) are determined. Based on the similarity between the molecular structures of PBBs and PCBs, a simple linear predictive model has been developed with the correlation coefficient R=0.9812 and standard error SE=0.19. The logK(oc) values of any PBB congeners can be predicted by using the logK(oc) values of the corresponding PCBs according to this model. Using the published data for logK(oc) values of PCB congeners, logK(oc) values of all 209 PBB congeners have been for the first time predicted. Compared with the data obtained from the experiment, the results of prediction are very accurate.