Computational toxicological investigation on the mechanism and pathways of xenobiotics metabolized by cytochrome P450: a case of BDE-47

Distinct Species-Specific and Toxigenic Metabolic Profiles for 6PPD and 6PPD Quinone by P450 Enzymes: Insights from In Vitro and In Silico Studies

The tire rubber antioxidant N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine (6PPD) and its quinone product (6PPDQ) are prevalent emerging contaminants, yet their biotransformation profiles remain poorly understood, hampering the assessment of environmental and health risks. This study investigated the phase-I metabolism of 6PPD and 6PPDQ across aquatic and mammalian species through in vitro liver microsome (LM) incubations and in silico simulations. A total of 40 metabolites from seven pathways were identified using the highly sensitive nano-electrospray ionization mass spectrometry. Notably, 6PPDQ was consistently detected as a 6PPD metabolite with an approximate 2% yield, highlighting biotransformation as a neglected indirect exposure pathway for 6PPDQ in organisms. 6PPDQ was calculated to form through a facile two-step phenyl hydroxylation of 6PPD, catalyzed by cytochrome P450 enzymes. Distinct species-specific metabolic kinetics were observed, with fish LM demonstrating retarded biotransformation rates for 6PPD and 6PPDQ compared to mammalian LM, suggesting the vulnerability of aquatic vertebrates to these contaminants. Intriguingly, two novel coupled metabolites were identified for 6PPD, which were predicted to exhibit elevated toxicity compared to 6PPDQ and result from C-N oxidative coupling by P450s. These unveiled metabolic profiles offer valuable insights for the risk assessment of 6PPD and 6PPDQ, which may inform future studies and regulatory actions.

Distribution, transformation, and toxic effects of the flame retardant tetrabromobisphenol S and its derivatives in the environment: A review

As widely used alternative brominated flame retardants, tetrabromobisphenol S (TBBPS) and its derivatives have attracted increasing amounts of attention in the field of environmental science. Previous studies have shown that TBBPS and its derivatives easily accumulate in environmental media and may cause risks to environmental safety and human health. Therefore, to explore the environmental behaviours of TBBPS and its derivatives, in this paper, we summarized relevant research on the distribution of these compounds in water, the atmosphere, soil and food/biota, as well as their transformation mechanisms (biological and nonbiological) and toxic effects. The summary results show that TBBPS and its derivatives have been detected in water, the atmosphere, soil, and food/biota globally, making them a ubiquitous pollutant. These compounds may be subject to adsorption, photolysis or biological degradation after being released into the environment, which in turn increases their ecological risk. TBBPS and its derivatives can cause a series of toxic effects, such as neurotoxicity, hepatotoxicity, cytotoxicity, thyrotoxicity, genotoxicity and phytotoxicity, to cells or living organisms in in vitro and in vivo exposure. Toxicological studies suggest that TBBPS as an alternative to TBBPA is not entirely environmentally friendly. Finally, we propose future directions for research on TBBPS and its derivatives, including the application of new technologies in studies on the migration, transformation, toxicology and human exposure risk assessment of TBBPS and its derivatives in the environment. This review provides useful information for obtaining a better understanding of the behaviour and potential toxic effects of TBBPS and its derivatives in the environment.

ROS-dependent catalytic mechanism of melatonin metabolism and its application in the measurement of reactive oxygen

Melatonin (Mel) is an endogenous active molecule whose metabolism progress significantly influences its bioactivity. However, the detailed metabolic pathway of Mel in the pathological state has not yet been fully illustrated. In this study, 16 metabolites of Mel in cancer cells and human liver microsomes were identified, of which seven novel metabolites were newly discovered. Among them, 2-hydroxymelatonin (2-O-Mel), as the major metabolite in cancer cells, was revealed for the first time, which was different from the metabolite found in the human liver. Furthermore, CYP1A1/1A2- and reactive oxygen species (ROS)-mediated 2-hydroxylation reactions of Mel were verified to be the two metabolic pathways in the liver and cancer cells, respectively. ROS-dependent formation of 2-O-Mel was the major pathway in cancer cells. Furthermore, the underlying catalytic mechanism of Mel to 2-O-Mel in the presence of ROS was fully elucidated using computational chemistry analysis. Therefore, the generation of 2-O-Mel from Mel could serve as another index for the endogenous reactive oxygen level. Finally, based on the ROS-dependent production of 2-O-Mel, Mel was successfully used for detecting the oxygen-carrying capacity of hemoglobin in human blood. Our investigation further enriched the metabolic pathway of Mel, especially for the ROS-dependent formation of 2-O-Mel that serves as a diagnostic and therapeutic target for the rational use of Mel in clinics.

2,2',4,4'-Tetrabromodiphenyl ether and cadmium co-exposure activates aryl hydrocarbon receptor pathway to induce ROS and GSDME-dependent pyroptosis in renal tubular epithelial cells

We have previously found that a mixture exposure of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) and cadmium (Cd) causes kidney damage; however, the mechanism was not fully understood. The aryl hydrocarbon receptor (AhR) is a ligand-receptor transcription factor that plays an important role in the adaptive response or metabolic detoxification of environmental toxins. Thus, this study aimed to examine the role of AhR in kidney toxicity. BDE-47 (50 μM) or Cd (5 μM) exposure reduced cell viability in renal tubular epithelial cells (HKC), with a larger effect observed in co-treatment. The cell morphology presented pyroptotic changes, including swollen cells, large bubbles, and plasma membrane pore formation. The gene expressions of AhR, heat shock protein 90 (Hsp90), AhR nuclear translocator (ARNT), and cytochrome P450 1B1 (CYP1B1) were increased, while CYP1A1 was decreased. Reactive oxygen species (ROS) were generated, which was reduced by the AhR antagonist CH223191. The apoptosis, necrosis, and intracellular lactated hydrogenase (LDH) release was elevated, and this was attenuated by N-acetylcysteine (NAC). Furthermore, the pyroptosis pathway was activated with increased protein levels of cleaved-caspase-3 and gasdermin E N-terminal (GSDME-NT), while caspase-8, caspase-3, and GSDME were decreased. These effects were alleviated by NAC and CH223191. Our data demonstrate a combined effect of BDE-47 and Cd on nephrotoxicity by activating AhR to induce ROS contributing to GSDME-dependent pyroptosis, and retardation of the AhR pathway could reduce this toxicity.

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (ε_bulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (ε_bulk^') from intrinsic KIEs on the rate-determining steps and the observed ε_bulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and S_N2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

Origin of metabolites diversity and selectivity of P450 catalyzed benzo[a]pyrene metabolic activation

Polycyclic Aromatic Hydrocarbon (PAHs) presents one of the most abundant class of environmental pollutants. Recent study shows a lab-synthesized PAHs derivative, helicenium, can selectively kill cancer cells rather than normal cells, calling for the in-depth understanding of the metabolic process. However, the origin of metabolites diversity and selectivity of P450 catalyzed PAHs metabolic activation is still unclear to a great extent. Here we systematically investigated P450 enzymes catalyzed activation mechanism of a representative PAHs, benzo[a]pyrene (BaP), and found the corresponding activation process mainly involves two elementary steps: electrophilic addition and epoxidation. Electrophilic addition step is evidenced to be rate determining step. Two representative binding modes of BaP with P450 were found, which enables the electrophilic addition of Heme (FeO) to almost all the carbons of BaP. This electrophilic addition was proposed to be accelerated by the P450 enzyme environment when compared with the gas phase and water solvent. To dig deeper on the origin of metabolites diversity, we built several linear regression models to explore the structural-energy relationships. The selectivity was eventually attributed to the integrated effects of structural (e.g. O-C distance and O-C-Fe angle) and electrostatic parameters (e.g. charge of C and O) from both BaP and P450.

Computational Insight into Biotransformation Profiles of Organophosphorus Flame Retardants to Their Diester Metabolites by Cytochrome P450

Biotransformation of organophosphorus flame retardants (OPFRs) mediated by cytochrome P450 enzymes (CYPs) has a potential correlation with their toxicological effects on humans. In this work, we employed five typical OPFRs including tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(1-chloro-2-propyl) phosphate (TCIPP), tri(2-chloroethyl) phosphate (TCEP), triethyl phosphate (TEP), and 2-ethylhexyl diphenyl phosphate (EHDPHP), and performed density functional theory (DFT) calculations to clarify the CYP-catalyzed biotransformation of five OPFRs to their diester metabolites. The DFT results show that the reaction mechanism consists of Cα-hydroxylation and O-dealkylation steps, and the biotransformation activities of five OPFRs may follow the order of TCEP ≈ TEP ≈ EHDPHP > TCIPP > TDCIPP. We further performed molecular dynamics (MD) simulations to unravel the binding interactions of five OPFRs in the CYP3A4 isoform. Binding mode analyses demonstrate that CYP3A4-mediated metabolism of TDCIPP, TCIPP, TCEP, and TEP can produce the diester metabolites, while EHDPHP metabolism may generate para-hydroxyEHDPHP as the primary metabolite. Moreover, the EHDPHP and TDCIPP have higher binding potential to CYP3A4 than TCIPP, TCEP, and TEP. This work reports the biotransformation profiles and binding features of five OPFRs in CYP, which can provide meaningful clues for the further studies of the metabolic fates of OPFRs and toxicological effects associated with the relevant metabolites.

Insights into the metabolic mechanism of PBDEs catalyzed by cytochrome P450 enzyme 3A4: A QM/MM study

Elucidating the metabolic mechanism and the derivatives of polybrominated diphenyl ethers (PBDEs) is significant to risk assessment. This study delineated the metabolic mechanism of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) catalyzed by P450 enzymes using a combination of molecular dynamic (MD), quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT). The calculation results reveal that the electrophilic addition is the main pathway for the biotransformation of BDE-47 catalyzed by P450 enzymes. 6-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (6-OH-BDE-47) is a more kinetically preferable product than 5-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (5-OH-BDE-47). Electrophilic addition reaction can lead to the formation of polybrominated dibenzo-p-dioxins and dibenzofurans (PBDD/Fs). The ecotoxicity assessment indicates that the final products of BDE-47 are still toxic to aquatic organisms, but the solubility increase of the hydroxylated products can accelerate their excretion from the body. We expect that the established metabolic mechanism and the derivatives will be used to predict the biotransformation of other PBDE congeners catalyzed by P450 enzymes in human livers.

Mechanistic insight into hydroxylation of 2,2',4,4'-tetrabromodiphenyl ether during biodegradation by typical aerobic bacteria: Experimental and computational studies

Polybrominated diphenyl ethers (PBDEs) are a group of persistent pollutants in the environment. Though aerobic biodegradation of PBDEs have been extensively studied, the involved hydroxylation mechanism decisive for whole biotransformation is not clear yet. During the effective biodegradation of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) by B. xenovorans LB400, the depletion of endogenous ∙OH by scavenger could bring about the significant decrease of biodegradation efficiency whereas ·O₂^- was nearly not influential. Given the importance of ∙OH in hydroxylation, the reaction mechanisms along major pathways of electrophilic addition and hydrogen abstraction were theoretically examined by density functional theory (DFT). For the less demand of activation energy, the relative preference of electrophilic addition was shown at aromatic C3-site. When the secondary reaction was considered after addition at C4-site, the barrierless association of ∙OH at C3-site and deprotonation by H₂O was validated as the energetically-favorable pathway that may cause dihydroxylation of BDE-47 into 3,4-dihydroxyl-BDE-17. The electrophilic addition followed by seconary barrierless trans-association of ∙OH and then dehydration seemed favorable for monohydroxylation as regards energetic barrier merely up to 194.01 kJ mol^-1, while the hydrogen abstraction by ∙OH from C5-site was more privileged actually. The theoretical insights would help well understand the hydroxylation mechanism of PBDEs by aerobes.

Ferrate(VI) Oxidation Mechanism of Substituted Anilines: A Density Functional Theory Investigation

Ferrate(VI) (Fe(VI)) is a promising oxidant coagulant and disinfectant for the degradation of organic micropollutants. However, it is hard to elucidate the detailed oxidation mechanism through the current experimental approaches. Substituted anilines (SANs) are important chemical compounds that are widely used in many industries. This paper presents the use of density functional theory (DFT) to understand the oxidation mechanism of SANs by Fe(VI) and the effect of substituents. The calculation results revealed that the primary oxidations of SANs follow the hydrogen atom transfer (HAT) mechanism. Interestingly, the hydroxyl oxygen of HFeO₄ ^- is more reactive than the carbonyl oxygen when reacting with SANs. The formation of the SAN radical is crucial, and all of the products are formed from it. Azobenzene is more favorable to generate the above products. In addition, the obtained results indicate that this kind of substituent has a much greater influence on the reaction rather than the position. Thus, the present study provides a valuable insight into the transformation pathways of SANs in the Fe(VI) oxidation process and the effects of the substituent on oxidation. These results will advance the understanding of Fe(VI) involved in wastewater treatment.

Investigating the molecular mechanism of hydroxylated bromdiphenyl ethers to inhibit the thyroid hormone sulfotransferase SULT1A1

Hydroxylated bromodiphenyl ethers (OH-BDEs) have raised great concern due to their potential endocrine disrupting effects on humans. In vitro experiments have indicated OH-BDEs can inhibit the activity of thyroid hormone (TH) sulfotransferases (SULTs); however, the molecular mechanism has not been investigated in depth. In this work, we employed 17 OH-BDEs with five or fewer Br atoms, and performed integrated computational simulations to unravel the possible inhibition mechanism of OH-BDEs on human SULT1A1. The molecular docking results demonstrate that OH-BDEs form hydrogen bonds with residues Lys106 and His108, and the neutral OH-BDEs show comparable binding energies with their anionic counterparts. The further hybrid quantum mechanical/molecular mechanical (QM/MM) calculations unravel a metabolic mechanism of OH-BDEs comprised by proton abstraction and sulfation steps. This mechanism is involved in the SULT1A1 inhibition by some OH-BDEs comprised of three or fewer Br atoms, while other OH-BDEs likely only form ternary complexes to competitively inhibit SULT1A1 activity. Moreover, the effect of the hydroxyl group of OH-BDEs on SULT1A1 inhibition potential follows the order of ortho-OH BDE > meta-OH BDE > para-OH BDE. These results provide an insight into the inhibition mechanism of OH-BDEs to SULT1A1 at the molecular level, which are beneficial in illuminating the molecular initiating events involved in the TH disruption of OH-BDEs.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

Reductive Debromination of Polybrominated Diphenyl Ethers: Dependence on Br Number of the Br-Rich Phenyl Ring

Reductive debromination has been widely studied for the degradation of polybrominated diphenyl ethers (PBDEs), although the reaction mechanisms are not so clear. In the present study, the photocatalytic degradation and debromination of ten PBDEs were carried out with CuO/TiO₂ nanocomposites as the photocatalyst under an anaerobic condition. The pseudo-first-order rate constants were obtained for the photocatalytic debromination of PBDEs, and their relative rate constants ( k_R) were evaluated against k_R= 1 for BDE209. Unlike the generally accepted summary that k_R is dependent on the total Br number ( N) of PBDEs, k_R is found to depend on the Br number on a phenyl ring with more Br atoms than the other one. In other words, a phenyl ring substituted by more Br is more reactive for the reductive debromination. The calculated LUMO energies ( E_LUMO) of all PBDEs are well correlated to the more reactive phenyl ring with more Br, compared with the N of two phenyl rings. The result was explained by LUMO localization on the Br-rich phenyl ring, suggesting that the reductive debromination occurs on the phenyl ring.

Roles of maize cytochrome P450 (CYP) enzymes in stereo-selective metabolism of hexabromocyclododecanes (HBCDs) as evidenced by in vitro degradation, biological response and in silico studies

In vitro biotransformation of HBCDs by maize cytochrome P450 (CYP) enzymes, responses of CYPs to HBCDs at protein and transcription levels, and in silico simulation of interactions between CYPs and HBCDs were investigated in order to elucidate the roles of CYPs in the metabolism of HBCDs in maize. The results showed that degradation reactions of HBCDs by maize microsomal CYPs followed the first-order kinetics and were stereo-selective, with the metabolic rates following the order (-)γ- > (+)γ- > (+)α- > (-)α-HBCD. The hydroxylated metabolites OH-HBCDs, OH-PBCDs and OH-TBCDs were detected. (+)/(-)-α-HBCDs significantly decreased maize CYP protein content and inhibited CYP enzyme activity, whereas (+)/(-)-γ-HBCDs had obvious effects on the induction of CYPs. HBCDs selectively mediated the gene expression of maize CYPs, including the isoforms of CYP71C3v2, CYP71C1, CYP81A1, CYP92A1 and CYP97A16. Molecular docking results suggested that HBCDs could bind with these five CYPs, with binding affinity following the order CYP71C3v2 < CYP81A1 < CYP97A16 < CYP92A1 < CYP71C1. The shortest distances between the Br-unsubstituted C atom of HBCD isomers and the iron atom of heme in CYPs were 4.18-11.7 Å, with only the distances for CYP71C3v2 being shorter than 6 Å (4.61-5.38 Å). These results suggested that CYP71C3v2 was an efficient catalyst for degradation of HBCDs. For (+)α- and (-)γ-HBCDs, their binding affinities to CYPs were lower and the distances to the iron atom of heme in CYPs were shorter than their corresponding antipodes, consistent with the in vitro experimental results showing that they had shorter half-lives and were more easily hydroxylated. This study provides solid evidence for the roles of maize CYPs in the metabolism of HBCDs, and gives insight into the molecular mechanisms of the enantio-selective metabolism of HBCDs by plant CYPs.

Carcinogenic Metabolic Activation Process of Naphthalene by the Cytochrome P450 Enzyme 1B1: A Computational Study

The metabolic activation and transformation of naphthalene by the cytochrome P450 enzyme (CYP 1B1) plays an important role in its potential carcinogenicity. The process has been explored by a quantum mechanics/molecular mechanics (QM/MM) computational method. Molecular dynamic simulations were performed to explore the interaction between naphthalene and CYP 1B1. Naphthalene involves α- and β-carbon, the electrophilic addition of which would result in different reaction pathways. Our computational results show that both additions on α- and β-carbon can generate naphthalene 1,2-oxide. The activation barrier for the addition on β-carbon is higher than that for the α-carbon by 2.6 kcal·mol^-1, which is possibly caused by the proximity between β-carbon and the iron-oxo group of Cpd I in the system. We also found that naphthalene 1,2-oxide is unstable and the O-C bond cleavage easily occurs via cellular hydronium ion, hydroxyl radical/anion; then it will convert to the potential ultimate carcinogen 1,2-naphthoquinone. The results demonstrate and inform a detailed process of generating naphthalene 1,2-oxide and new predictions for its conversion.

Metabolic Mechanism of Aryl Phosphorus Flame Retardants by Cytochromes P450: A Combined Experimental and Computational Study on Triphenyl Phosphate

Understanding metabolic mechanisms is critical and remains a difficult task in the risk assessment of emerging pollutants. Triphenyl phosphate (TPHP), a widely used aryl phosphorus flame retardant (aryl-PFR), has been frequently detected in the environment, and its major metabolite was considered as diphenyl phosphate (DPHP). However, knowledge of the mechanism for TPHP leading to DPHP and other metabolites is lacking. Our in vitro study shows that TPHP is metabolized into its diester metabolite DPHP and mono- and dihydroxylated metabolites by cytochromes P450 (CYP) in human liver microsomes, while CYP1A2 and CYP2E1 isoforms are mainly involved in such processes. Molecular docking gives the conformation for TPHP binding with the active species Compound I (an iron IV-oxo heme cation radical) in specific CYP isoforms, showing that the aromatic ring of TPHP is likely to undergo metabolism. Quantum chemical calculations have shown that the dominant reaction channel is the O-addition of Compound I onto the aromatic ring of TPHP, followed by a hydrogen-shuttle mechanism leading to ortho-hydroxy-TPHP as the main monohydroxylated metabolite; the subsequent H-abstraction-OH-rebound reaction acting on ortho-hydroxy-TPHP yields the meta- and ipso-position quinol intermediates, while the former of which can be metabolized into dihydroxy-TPHP by fast protonation, and the latter species needs to go through type-I ipso-substitution and fast protonation to be evolved into DPHP. We envision that the identified mechanisms may give inspiration for studying the metabolism of several other aryl-PFRs by CYP.

Quantum chemical simulations revealed the toxicokinetic mechanisms of organic phosphorus flame retardants catalyzed by P450 enzymes

The metabolic fate and toxicokinetics of organic phosphorus flame retardants catalyzed by cytochrome P450 enzymes (CYPs) are here investigated by in silico simulations, leveraging an active center model to mimic the CYPs, triphenyl phosphate (TPHP), tris(2-butoxyethyl) phosphate and tris(1,3-dichloro-2-propyl) phosphate as substrates. Our calculations elucidated key main pathways and predicted products, which were corroborated by current in vitro data. Results showed that alkyl OPFRs are eliminated faster than aryl and halogenated alkyl-substituted OPFRs. In addition, we discovered a proton shuttle pathway for aryl hydroxylation of TPHP and P = O bond-assisted H-transfer mechanisms (rather than nonenzymatic hydrolysis) that lead to O-dealkylation/dearylation of phosphotriesters.

Computational Insight into the Activation Mechanism of Carcinogenic N'-Nitrosonornicotine (NNN) Catalyzed by Cytochrome P450

Tobacco-specific N'-nitrosonornicotine (NNN), a genotoxic nitrosamine classified as Group 1 carcinogen, is also present in atmospheric particulate matter and has even been detected as a new disinfection byproduct in wastewaters. NNN generally requires metabolic activation by cytochrome P450 enzymes to exert its genotoxicity, but the respective biotransformation pathways have not been described in detail. In this work, we performed density functional theory (DFT) calculations to unravel possible NNN activation pathways including α-hydroxylation, β-hydroxylation, pyridine N-oxidation, and norcotinine formation. The results reveal an initial rate-determining H_α-atom abstraction step for α-hydroxylation, followed by an unexpected kinetic competition between denitrosation and OH rebound, leading to ( iso-)myosmine as a detoxified product and α-hydroxyNNNs as the precursor of carcinogenic diazohydroxides, respectively. Further detoxification routes are given by β-hydroxylation with relative high reaction barrier and N-oxidation with comparable barrier to the toxifying α-hydroxylation. Moreover, we show for the first time how norcotinine can be generated as a minor NNN metabolite that is formed from iso-myosmine through a unique porphyrin-assisted H atom 1,2-transfer mechanism. These results demonstrate that the carcinogenic potential of NNN is subject to a kinetic competition between activating and deactivating metabolic routes, and identify respective biomarkers to inform about the individual risk associated with NNN exposure.

Identification of Unknown Brominated Bisphenol S Congeners in Contaminated Soils as the Transformation Products of Tetrabromobisphenol S Derivatives

Compared with tetrabromobisphenol A (TBBPA) and its derivatives, the skeletally similar chemicals tetrabromobisphenol S (TBBPS) and its derivatives have been rarely studied, and very little is known about their structures, environmental occurrence, and behaviors. In this study, a total of 84 soil samples from a chemical industrial park have been collected and analyzed to investigate the occurrence of TBBPS and its derivatives and to identify novel TBBPS analogs. TBBPS, TBBPS bis(2,3-dibromopropyl ether) (TBBPS-BDBPE), and three byproducts, TBBPS mono(allyl ether) (TBBPS-MAE), TBBPS mono(2-bromoallyl ether) (TBBPS-MBAE), and TBBPS mono(2,3-dibromopropyl ether) (TBBPS-MDBPE), have been detected with contents ranging from below detection limits to 1934.6 ng/g dw and with detection frequencies of 21.4-97.6%. In addition, another 5 unknown TBBPS analogs, tribromobisphenol S (TriBBPS), 2,2',6'-TriBBPS-MAE (TriBBPS-MAE_3.2), 2,6,2'-TriBBPS-MAE (TriBBPS-MAE_3.4), 2',6'-DBBPS-MAE (DBBPS-MAE_2.0), and 2,6-DBBPS-MAE (DBBPS-MAE_2.6), have been identified in these soil samples by untargeted mass spectrometry screening. These unknown analogs have also been observed in laboratory transformation experiments of TBBPS-MDBPE conducted under reducing conditions. TriBBPS-MAE_3.4 and DBBPS-MAE_2.6 were more likely to be produced than TriBBPS-MAE_3.2 and DBBPS-MAE_2.0 due to the stereoselectivity of the transformation. TriBBPS-MAE_3.4 and DBBPS-MAE_2.0 were more stable, resulting in higher detection frequencies of these compounds in soil samples. Ether bond breakage and debromination contributed to the generation of these novel products. The results provide new information on the behaviors of TBBPS and its derivatives in the environment.

Plant Uptake and Metabolism of 2,4-Dibromophenol in Carrot: In Vitro Enzymatic Direct Conjugation

Plants can extensively uptake organic contaminants from soil and subsequently transform them into various products. Those compounds containing hydroxyl may undergo direct conjugation with endogenous biomolecules in plants, and potentially be preserved as conjugates, thus enabling overlooked risk via consumptions of food crops. In this study, we evaluated the uptake and metabolism of 2,4-dibromophenol (DBP) by both carrot cells and whole plant. DBP was completely removed from cell cultures with a half-life of 10.8 h. Four saccharide conjugates, three amino acid conjugates, and one phase I metabolite were identified via ultraperformance liquid chromatography quadrupole time-of-flight mass spectrometry analysis. The dibromophenol glucopyranoside (glucose conjugate) was quantitated by synthesized standard and accounted for 9.3% of the initial spiked DBP at the end of incubation. The activity of glycosyltransferase was positively related to the production of 2,4-dibromophenol glucopyranoside ( p = 0.02, R² = 0.86), implying the role of enzymatic catalysis involved in phase II metabolism.

Oxidation reactivity of 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE) by Compound I model of cytochrome P450s

Alternative brominated flame retardants (BFRs) have become prevalent as a consequence of restrictions on the use of polybrominated diphenyl ethers (PBDEs). For risk assessment of these alternatives, knowledge of their metabolism via cytochrome P450 enzymes is needed. We have previously proved that density functional theory (DFT) is able to predict the metabolism of PBDEs by revealing the molecular mechanisms. In the current study, the reactivity of 1,2-bis(2,4,6-tribromophenoxy)ethane and structurally similar chemicals with the Compound I model representing the active site of P450 enzymes was investigated. The DFT calculations delineated reaction pathways which lead to reasonable explanations for products that were detected by wet experiments, meanwhile intermediates which cannot be determined were also proposed. Results showed that alkyl hydrogen abstraction will lead to bis(2,4,6-tribromophenoxy)ethanol, which may undergo hydrolysis yielding 2,4,6-tribromophenol, a neurotoxic compound. In addition, a general pattern of oxidation reactivity regarding the 2,4,6-tribromophenyl moiety was observed among several model compounds. Our study has provided insights for convenient evaluation of the metabolism of other structurally similar BFRs.

Identification of Emerging Brominated Chemicals as the Transformation Products of Tetrabromobisphenol A (TBBPA) Derivatives in Soil

In contrast to the extensive investigation already conducted on tetrabromobisphenol A (TBBPA), the metabolism of TBBPA derivatives is still largely unknown. In this paper, we characterized unknown brominated compounds detected in 84 soil samples collected from sites around three brominated flame retardant production plants to determine possible transformation products of TBBPA derivatives. In addition to tribromobisphenol A (TriBBPA), dibromobisphenol A (DBBPA), and TBBPA, six novel transformation products, TriBBPA mono(allyl ether) (TriBBPA-MAE), DBBPA-MAE, hydroxyl TriBBPA-MAE, TBBPA mono(2-bromo-3-hydroxypropyl ether) (TBBPA-MBHPE), TBBPA mono(2,3-dihydroxypropyl ether) (TBBPA-MDHPE), and TBBPA mono(3-hydroxypropyl ether) (TBBPA-MHPE) were identified. The detection frequencies of these identified chemicals in soil samples ranged from 17% to 89%, indicating the widespread presence of the transformation products. To uncover the possible TBBPA derivative transformation pathways involved, super-reduced vitamin B12 (cyanocobalamin, (CCAs)) was used to treat TBBPA derivative and transformation products in this process were characterized. To our knowledge, this is the first study examining the transformation of TBBPA derivatives and the first to report several novel associated TBBPA and bisphenol A derivatives as transformation products. Our research suggests that ether bond breakage and debromination contribute to the transformation of TBBPA derivatives and the existence of the novel transformation products. These data provide new insights into the fate of TBBPA derivatives in environmental compartments.

Effects of Atmospheric Water on ·OH-initiated Oxidation of Organophosphate Flame Retardants: A DFT Investigation on TCPP

Tris(2-chloroisopropyl) phosphate (TCPP), a widely used organophosphate flame retardant, has been recognized as an important atmospheric pollutant. It is notable that TCPP has potential for long-range atmospheric transport. However, its atmospheric fate is unknown, restricting its environmental risk assessment. Herein we performed quantum chemical calculations to investigate the atmospheric transformation mechanisms and kinetics of TCPP initiated by ·OH in the presence of O₂/NO/NO₂, and the effects of ubiquitous water on these reactions. Results show the H-abstraction pathways are the most favorable for the titled reaction. The calculated gaseous rate constant and lifetime at 298 K are 1.7 × 10^-10 cm³molecule^-1 s^-1 and 1.7 h, respectively. However, when considering atmospheric water, the corresponding lifetime is about 0.5-20.2 days. This study reveals for the first time that water has a negative role in the ·OH-initiated degradation of TCPP by modifying the stabilities of prereactive complexes and transition states via forming hydrogen bonds, which unveils one underlying mechanism for the observed persistence of TCPP in the atmosphere. Water also influences secondary reaction pathways of selected TCPP radicals formed from the primary H-abstraction. These results demonstrate the importance of water in the evaluation of the atmospheric fate of newly synthesized chemicals and emerging pollutants.

Ferrate(vi) initiated oxidative degradation mechanisms clarified by DFT calculations: a case for sulfamethoxazole

Ferrate(vi) is an efficient and environmentally friendly oxidant for the degradation of organic micropollutants. However, the related mechanism for the degradation is ambiguous and can hardly be elucidated empirically due to the rapid oxidation process and unstable intermediates for experimental trapping. Herein we performed density function theory (DFT) calculations to unveil the mechanism of ferrate(vi)-mediated degradation, taking sulfamethoxazole as a model compound. The results show that nucleophilic attack (rather than electrophilic attack) of HFeO₄^- on the isoxazole moiety of sulfamethoxazole initiates the subsequent degradations, and ferrate(vi) rather than the water molecule provides O atoms for the oxidation of the nitroso group and isoxazole moiety. Electron delocalization from the Fe atom to the isoxazole moiety is crucial for the ring-opening of isoxazole, and organometallic intermediates suggested previously are not the necessary ones in the oxidation of sulfamethoxazole by HFeO₄^-. Thus, this study has theoretically clarified the ferrate(vi) oxidation mechanisms for a representative sulfonamide, which were also partially corroborated by the intermediates and products observed in the previous experimental studies for phosphite and tryptophan. This study provides an exemplification on the application of quantum chemical calculations to clarify the degradation pathways of organic micropollutants, which is important for the prediction of degradation products needed in their engineering design.

How PBDEs Are Transformed into Dihydroxylated and Dioxin Metabolites Catalyzed by the Active Center of Cytochrome P450s: A DFT Study

Predicting metabolism of chemicals and potential toxicities of relevant metabolites remains a vital and difficult task in risk assessment. Recent findings suggested that polybrominated diphenyl ethers (PBDEs) can be transformed into dihydroxylated and dioxin metabolites catalyzed by cytochrome P450 enzymes (CYPs), whereas the mechanisms pertinent to these transformations remain largely unknown. Here, by means of density functional theory (DFT) calculations, we probed the metabolic pathways of 2,2',4,4'-tetraBDE (BDE-47) using the active center model of CYPs (Compound I). Results show that BDE-47 is first oxidized to monohydroxylated products (HO-BDEs), wherein a keto-enol tautomerism is identified for rearrangement of the cyclohexenone intermediate. Dihydroxylation with HO-BDEs as precursors, has a unique phenolic H-abstraction and hydroxyl rebound pathway that is distinct from that for monohydroxylation, which accounts for the absence of epoxides in in vitro studies. Furthermore, we found only dihydroxylated PBDEs with heterophenyl -OH substituents ortho- and meta- to the ether bond serve as precursors for dioxins, which are evolved from aryl biradical coupling of diketone intermediates that are produced from dehydrogenation of the dihydroxylated PBDEs by Compound I. This study may enlighten the development of computational models that afford mechanism-based prediction of the xenobiotic biotransformation catalyzed by CYPs.

Modeling of Toxicity-Relevant Electrophilic Reactivity for Guanine with Epoxides: Estimating the Hard and Soft Acids and Bases (HSAB) Parameter as a Predictor

According to the electrophilic theory in toxicology, many chemical carcinogens in the environment and/or their active metabolites are electrophiles that exert their effects by forming covalent bonds with nucleophilic DNA centers. The theory of hard and soft acids and bases (HSAB), which states that a toxic electrophile reacts preferentially with a biological macromolecule that has a similar hardness or softness, clarifies the underlying chemistry involved in this critical event. Epoxides are hard electrophiles that are produced endogenously by the enzymatic oxidation of parent chemicals (e.g., alkenes and PAHs). Epoxide ring opening proceeds through a SN2-type mechanism with hard nucleophile DNA sites as the major facilitators of toxic effects. Thus, the quantitative prediction of chemical reactivity would enable a predictive assessment of the molecular potential to exert electrophile-mediated toxicity. In this study, we calculated the activation energies for reactions between epoxides and the guanine N7 site for a diverse set of epoxides, including aliphatic epoxides, substituted styrene oxides, and PAH epoxides, using a state-of-the-art density functional theory (DFT) method. It is worth noting that these activation energies for diverse epoxides can be further predicted by quantum chemically calculated nucleophilic indices from HSAB theory, which is a less computationally demanding method than the exacting procedure for locating the transition state. More importantly, the good qualitative/quantitative correlations between the chemical reactivity of epoxides and their bioactivity suggest that the developed model based on HSAB theory may aid in the predictive hazard evaluation of epoxides, enabling the early identification of mutagenicity/carcinogenicity-relevant SN2 reactivity.

Eco-toxicity and human estrogenic exposure risks from OH-initiated photochemical transformation of four phthalates in water: A computational study

Transformation products (TPs) of emerging organic contaminates (EOCs) in water are still rarely considered in environmental risk assessment, although some have been found to be concern. OH is believed as an important reactive species both in indirect phototransformation and advanced oxidation technology. Thus, eco-toxicity and human estrogenic exposure risks of four phthalates and TPs during the OH-initiated photochemical process were investigated using computational approach. Four phthalates can be degraded through OH-addition and H-transfer pathways. The OH-addition TPs were predominant for dimethyl phthalates, while H-transfer TPs were predominant for other three phthalates. Compared with phthalates, OH-addition TPs (o-OH-phthalates) were one level more toxic to aquatic organisms, and m-OH-phthalates exhibit higher estrogenic activity. Although H-transfer TPs were less harmful than OH-addition TPs, some of them still have aquatic toxicity and estrogenic activity. Therefore, more attentions should be paid to photochemical TPs and original EOCs, particularly those exhibiting high estrogenic activity to humans.

Modeling and predicting pKa values of mono-hydroxylated polychlorinated biphenyls (HO-PCBs) and polybrominated diphenyl ethers (HO-PBDEs) by local molecular descriptors

Hydroxylated polychlorinated biphenyls (HO-PCBs) and polybrominated diphenyl ethers (HO-PBDEs) are attracting considerable concerns because of their multiple endocrine-disrupting effects and wide existence in environment and organisms. The hydroxyl groups enable these chemicals to be ionizable, and dissociation constant, pKa, becomes an important parameter for investigating their environmental behavior and biological activities. In this study, a new pKa prediction model was developed using local molecular descriptors. The dataset contains 21 experimental pKa values of HO-PCBs and HO-PBDEs. The optimized geometries by ab initio HF/6-31G(∗∗) algorithm were used to calculate the site-specific molecular readiness to accept or donate electron charges. The developed model obtained good statistical performance, which significantly outperformed commercial software ACD and SPARC. Mechanism analysis indicates that pKa values increase with the charge-limited donor energy EQocc on hydroxyl oxygen atom and decrease with the energy-limited acceptor charge QEvac on hydroxyl hydrogen atom. The regression model was also applied to calculate pKa values for all 837 mono-hydroxylated PCBs and PBDEs in each class, aiming to provide basic data for the ecological risk assessment of these chemicals.

Detoxification of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) by cytochrome P450 enzymes: A theoretical investigation

Two types of detoxification routes, N-demethylation to form 4-phenyl-1,2,3,6-tetrahydropyridine (PTP) and aromatic hydroxylation to generate 4-(4'-hydroxyphenyl)-1-methyl-1,2,3,6-tetrahydropyridine (MPTP-OH), for 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) mediated by Compound I (Cpd I) of cytochrome P450 are investigated theoretically using hybrid density functional calculations. Quantum chemical results reveal that for the N-demethylation, the initial C-H bond activation is achieved via a hydrogen atom transfer (HAT) mechanism. This is followed by a subsequent O-rebound to yield the carbinolamine intermediate. Due to the nature of pericyclic reaction, the generated carbinolamine decomposes in a non-enzymatic aqueous environment with the assistance of water molecules, forming amine and hydrated formaldehyde. For the aromatic hydroxylation, an initial addition of Cpd I to the substrate occurs mainly through a side-on approach with a subsequent proton shuttle to form the phenol product. A comparison of the energy barriers for both routes indicates that the N-demethylation (7.5/5.7kcal/mol for the quartet/doublet state in solvent) is thermodynamically more favorable than the aromatic hydroxylation process (14.9/14.8kcal/mol for the quartet/doublet state in solvent). This trend is in good agreement with the experimental product distribution, viz., the N-demethylation product PTP is more than the aromatic hydroxylation product MPTP-OH. Taken together, these observations not only enrich our knowledge on the mechanistic details of the N-dealkylation and the aromatic hydroxylation by P450s, but also provide certain insights into the metabolism of other analogous toxins.

Transformation pathways of MeO-PBDEs catalyzed by active center of P450 enzymes: a DFT investigation employing 6-MeO-BDE-47 as a case

Recent in vivo and in vitro experiments indicated that methoxylated polybrominated diphenyl ethers (MeO-PBDEs) can be biotransformed into hydroxylated PBDEs (HO-PBDEs) that are more toxic than PBDEs and MeO-PBDEs. Nevertheless, the enzymatic transformation mechanism is not clear. We hypothesized that cytochrome P450 enzymes (CYPs) play a key role in the transformation and employed the density functional theory calculations to unveil the mechanism. The transformation of a model compound, 6-MeO-BDE-47, catalyzed by the active center of CYPs (Compound I), was computed. For the first time, our results show that the energy barriers for the addition of Compound I to the C atoms on the phenyl of 6-MeO-BDE-47 are much higher than that for hydroxylation of the methoxyl, indicating that O-demethylation is a dominating metabolic pathway. This is in line with experimental observations performed by others. The pathways for the transformation of 6-MeO-BDE-47 catalyzed by Compound I were clarified. A C-H bond of the methoxyl is activated by Compound I, followed by radical rebound to form carbinol intermediates, then the carbinols decompose to form 6-HO-BDE-47 with the assistance of water molecules. The computational method can be potentially employed to develop models that predict biotransformation of xenobiotics catalyzed by CYPs.

Evaluation of BDE-47 hydroxylation metabolic pathways based on a strong electron-withdrawing pentafluorobenzoyl derivatization gas chromatography/electron capture negative ionization quadrupole mass spectrometry

Understanding the metabolic pathways of polybrominated diphenyl ethers (PBDEs) is a key issue in the evaluation of their cytotoxicity after they enter the biota. In order to obtain more information concerning the metabolic pathways of PBDEs, we developed a strong electron-withdrawing pentafluorobenzoyl (PFBoyl) derivatization capillary gas chromatography/electron capture negative ionization quadrupole mass spectrometry (GC/ECNI-qMS). PFBoyl esterification greatly improves separation of the metabolites of PBDEs such as hydroxylated PBDEs (OH-PBDEs) and bromophenols (BPs) metabolites in rat liver microsomes (RLMs). On the other hand, the strong electron-withdrawing property of PFBoyl derivatized on OH-PBDEs and/or BPs makes cleavage of the ester bond on ECNI easier resulting in higher abundance of the structure-informative characteristic fragment ions at a high m/z region, which facilitate the identification of OH-PBDEs metabolites. Subsequent quantification can be performed by monitoring not only 79Br- (or 81Br-) but also their characteristic fragment ions, achieving more accurate isotope dilution quantification using GC/ECNI-qMS. These merits allow us to identify totally 12 metabolites of BDE-47, a typical example of PBDEs, in the RLMs in vitro incubation systems. In addition to the already known metabolites of BDE-47, one dihydroxylated 3,6-di-OH-BDE-47 and one dihydroxylated 3,5-di-OH-tetrabrominated dioxin were found. Moreover, the second hydroxylation took place on the same bromophenyl ring, where the first hydroxyl group was located, and was further confirmed via the identification of the dihydroxylated 2',6'-di-OH-BDE-28 of an asymmetric 2'-OH-BDE-28. This methodological development and its subsequent findings of the metabolic pathways of BDE-47 provided experimental evidence for understanding its dioxin-like behavior and endocrine disrupting risk.

A computational approach predicting CYP450 metabolism and estrogenic activity of an endocrine disrupting compound (PCB-30)

Endocrine disrupting chemicals influence growth and development through interactions with the hormone system, often through binding to hormone receptors such as the estrogen receptor. Computational methods can predict endocrine disrupting chemical activity of unmodified compounds, but approaches predicting activity following metabolism are lacking. The present study uses a well-known environmental contaminant, PCB-30 (2,4,6-trichlorobiphenyl), as a prototype endocrine disrupting chemical and integrates predictive (computational) and experimental methods to determine its metabolic transformation by cytochrome P450 3A4 (CYP3A4) and cytochrome P450 2D6 (CYP2D6) into estrogenic byproducts. Computational predictions suggest that hydroxylation of PCB-30 occurs at the 3- or 4-phenol positions and leads to metabolites that bind more strongly than the parent molecule to the human estrogen receptor alpha (hER-α). Gas chromatography-mass spectrometry experiments confirmed that the primary metabolite for CYP3A4 and CYP2D6 is 4-hydroxy-PCB-30, and the secondary metabolite is 3-hydroxy-PCB-30. Cell-based bioassays (bioluminescent yeast expressing hER-α) confirmed that hydroxylated metabolites are more estrogenic than PCB-30. These experimental results support the applied model's ability to predict the metabolic and estrogenic fate of PCB-30, which could be used to identify other endocrine disrupting chemicals involved in similar pathways.

Fast liquid chromatographic-tandem mass spectrometric method using mixed-mode phase chromatography and solid phase extraction for the determination of 12 mono-hydroxylated brominated diphenyl ethers in human serum

Hydroxylated polybrominated diphenyl ethers (OH-PBDEs) are formed from the oxidative metabolism of polybrominated diphenyl ethers (PBDEs) in humans, rats and mice, but their quantitation in human blood and other matrices with liquid chromatography-mass spectrometric techniques has been a challenge. In this study, a novel analytical method was developed and validated using only 250 μL of human serum for the quantitation of twelve OH-PBDEs, fully chromatographically separated in a 15 min analytical run. This method includes two novel approaches: an enzymatic hydrolysis procedure and a chromatographic separation using a mixed mode chromatography column. The enzymatic hydrolysis (EH) was found critical for 4'-OH-BDE17, which was not detectable without it. For the sample clean up, a solid phase extraction protocol was developed and validated for the extraction of the 12 congeners from human serum. In addition, for the first time baseline resolution of two components was achieved that correspond to a single peak previously identified as 6'-OH-BDE99. The method was validated for linearity, accuracy, precision, matrix effects, limit of quantification, limit of detection, sample stability and overall efficiency. Recoveries (absolute and relative) ranged from 66 to 130% with relative standard deviations <21% for all analytes. Limit of detection and quantitation ranged from 4 to 90 pg mL(-1) and 6-120 pg mL(-1), respectively, with no carry over effects. This method was applied in ten commercially available human serum samples from the general US population. The mean values of the congeners detected in all samples are 4'-OH-BDE17 (34.2 pg mL(-1)), 4-OH-BDE42 (33.9 pg mL(-1)), 5-OH-BDE47 (17.5 pg mL(-1)) and 4'-OH-BDE49 (12.4 pg mL(-1)).

Comparative antioxidant status in freshwater fish Carassius auratus exposed to eight imidazolium bromide ionic liquids: a combined experimental and theoretical study

Imidazolium bromide ionic liquids such as 1-alkyl-3-methylimidazolium bromides ([AMIm]Br) and 1-alkyl-2,3-dimethylimidazolium bromides ([AMMIm]Br) are common-use organic salts. However, data on comparative toxicological effects of these ILs are lacking for fish. In this study, a combined experimental and theoretical approach was applied to compare and analyze the effects of these ILs on biochemical biomarkers in liver of Carassius auratus treated with different concentrations (2 and 20mg/L) for 3 and 16d. Changes in the activities of superoxide dismutase, catalase, glutathione peroxidase, and in the levels of reduced glutathione and malondialdehyde were detected, indicating that these ILs exhibit potential biotoxicity. The integrated biomarker response (IBR) index suggested that 1-hexyl-3-methylimidazolium bromide ([HMIm]Br), 1-octyl-3-methylimidazolium bromide ([OMIm]Br), 1-hexyl-2,3-dimethylimidazolium bromide ([HMMIm]Br), and 1-octyl-2,3-dimethylimidazolium bromide ([OMMIm]Br) showed the highest biotoxicity under different concentrations or exposure time, while 1-ethyl-3-methylimidazolium bromide ([EMIm]Br) always showed the least stressful power towards the test organism. Quantum chemical calculations (electronic parameters, frontier molecular orbitals, and Wiberg bond order) were also conducted to interpret the experimental results. Notably, some descriptors were correlated with the toxicity order. In addition, theoretical calculations provided some valuable information on metabolic pathways of these ILs, which may help to get better understanding on their environmental behavior and fate. In general, the toxicological determination and analysis of these ILs were performed with a combined experimental and theoretical method, which may contribute to the future ecotoxicological studies.

Transcriptomic responses of European flounder (Platichthys flesus) liver to a brominated flame retardant mixture

Male European flounder (Platichthys flesus) were exposed to a technical mixture of brominated diphenyl ethers (PDBEs, DE-71, Pentamix) that had been purified to remove contaminating dioxins. Controls were exposed to carrier solvent alone. Fish were exposed to decadally increasing concentrations of Pentamix via both sediment and spiked food. The GENIPOL P. flesus cDNA microarray, differentially expressed gene profiling (DEG) and quantitative PCR were employed to detect hepatic transcriptional differences between exposed fish and controls. Gene transcriptional changes were more sensitive to Pentamix exposure than biomarkers measured previously. Pentamix exposure induced transcripts coding for enzymes of xenobiotic metabolism (CYP1A, aldo-keto reductases) and elicited endocrine disruption (vitellogenin and thyroid hormone receptor alpha), with effects on CYP1A and VTG occurring at the highest exposure. Ontology analysis clearly showed dose-responsive changes indicative of oxidative stress, induction of mitochondrial dysfunction, and apoptosis. We conclude that exposure to PBDEs in both sediment and food has a significant adverse effect on a broad range of crucial biochemical processes in the livers of this widely distributed estuarine fish species, the flounder.

Comparative antioxidant status in freshwater fish Carassius auratus exposed to six current-use brominated flame retardants: a combined experimental and theoretical study

Decabromodiphenyl ether (BDE-209) and several non-polybrominated diphenyl ether (PBDE) brominated flame retardants (BFRs), such as tetrabromobisphenol A (TBBPA), hexabromocyclododecane (HBCD), decabromodiphenyl ethane (DBDPE), hexabromobenzene (HBB) and pentabromotoluene (PBT), are persistent halogenated contaminants ubiquitously detected in aquatic systems. However, data on comparative toxicological effects of these BFRs are lacking for fish. In this study, a combined experimental and theoretical approach was used to compare and analyze the effects of these BFRs on biochemical biomarkers in liver of Carassius auratus injected intraperitoneally with different doses (10 and 100mg/kg) for 7, 14 and 30 days. Oxidative stress was evoked evidently for the prolonged exposure, represented by the significantly altered indices (superoxide dismutase, catalase, glutathione peroxidase, reduced glutathione, and malondialdehyde). The integrated biomarker response (IBR) index ranked biotoxicity as: PBT>HBB>HBCD>TBBPA>BDE-209>DBDPE. Quantum chemical calculations (electronic parameters, frontier molecular orbitals, and Wiberg bond order) were performed for theoretical analysis. Notably, some descriptors were correlated with the toxicity order, probably implying the existence of a potential structure-activity relationship when more BFRs were included. Besides, theoretical calculations also provided some valuable information regarding the molecular characteristics and metabolic pathways of these current-use BFRs, which may facilitate the understanding on their environmental behavior and fate. Overall, this study adopted a combined experimental and theoretical method for the toxicological determination and analysis of the BFRs, which may also be considered in future ecotoxicological studies.

OH-initiated oxidation mechanisms and kinetics of 2,4,4'-Tribrominated diphenyl ether

2,4,4'-Tribromodiphenyl ether (BDE-28) was selected as a typical congener of polybrominated diphenyl ethers (PBDEs) to examine its fate both in the atmosphere and in water solution. All the calculations were obtained at the ground state. The mechanism result shows that the oxidations between BDE-28 and OH radicals are highly feasible especially at the less-brominated phenyl ring. Hydroxylated dibrominated diphenyl ethers (OH-PBDEs) are formed through direct bromine-substitution reactions (P1∼P3) or secondary reactions of OH-adducts (P4∼P8). Polybrominated dibenzo-p-dioxins (PBDDs) resulting from o-OH-PBDEs are favored products compared with polybrominated dibenzofurans (PBDFs) generated by bromophenols and their radicals. The complete degradation of OH adducts in the presence of O2/NO, which generates unsaturated ketones and aldehydes, is less feasible compared with the H-abstraction pathways by O2. Aqueous solution reduces the feasibility between BDE-28 and the OH radical. The rate constant of BDE-28 and the OH radical is determined to be 1.79 × 10(-12) cm(3) molecule(-1) s(-1) with an atmospheric lifetime of 6.7 days.

Distinct photolytic mechanisms and products for different dissociation species of ciprofloxacin

As many antibiotics are ionizable and may have different dissociation forms in the aquatic environment, we hypothesized that the different dissociation species have disparate photolytic pathways, products, and kinetics, and adopted ciprofloxacin (CIP) as a case to test this hypothesis. Simulated sunlight experiments and matrix calculations were performed to differentiate the photolytic reactivity for each dissociation species (H4CIP(3+), H3CIP(2+), H2CIP(+), HCIP(0), and CIP(-)). The results prove that the five dissociation species do have dissimilar photolytic kinetics and products. H4CIP(3+) mainly undergoes stepwise cleavage of the piperazine ring, while H2CIP(+) mainly undergoes defluorination. For H3CIP(2+), HCIP(0), and CIP(-), the major photolytic pathway is oxidation. By density functional theory calculation, we clarified the defluorination mechanisms for the five dissociation species at the excited triplet states: All the five species can defluorinate by reaction with hydroxide ions (OH(-)) to form hydroxylated products, and H2CIP(+) can also undergo C-F bond cleavage to produce F(-) and a carbon-centered radical. This study is a first attempt to differentiate the photolytic products and mechanisms for different dissociation species of ionizable compounds. The results imply that for accurate ecological risk assessment of ionizable emerging pollutants, it is necessary to investigate the environmental photochemical behavior of all dissociation species.