Pentabromoethylbenzene Exposure Induces Transcriptome Aberration and Thyroid Dysfunction: In Vitro, in Silico, and in Vivo Investigations

Selective cascade activation of polycyclic aromatic hydrocarbons in human cells: Role of enzyme's intrinsic electric field

Polycyclic aromatic hydrocarbons (PAHs) are major environmental organic pollutants. Some metabolites of PAHs show greater toxicity to humans while the others do not. It is highly important to decipher PAHs' regioselective activation mechanism and identify the major metabolites to accurately evaluate their public health risk. Here, we have performed a thorough computational study of benzo[a]anthracene (BA) metabolized by P450 1A1 by employing molecular docking, molecular dynamics simulations, quantum chemical calculation, and quantum mechanics/molecular mechanics calculations. Our findings show that highly-reactive species such as 3,4-epoxide, 8,9-epoxide, 3,4-diol-1,2-epoxide, and 8,9-diol-10,11-epoxide were major metabolites, which can efficiently react with guanine and damage DNA with extremely low energy barrier, therefore, supports the regioselective metabolism of BA. The origin of this selective activation is mainly contributed to both the oxygen‑carbon distance and previously overlooked enzyme's intrinsic electric field. Consequently, based on the resolved activation selectivity of BA. We built a high-throughput strategy to efficiently predict the metabolites of other PAHs. The accuracy of the strategy is validated by studying 16 PAHs on the priority control list. Hopefully this will aid the accurate evaluation of public health risks associated with PAH emissions.

Natural and anthropogenic organic brominated compounds in the southwestern Atlantic ocean: Bioaccumulation in coastal and oceanic dolphin species

Marine pollution is considered a current driver of change in the oceans and despite the urgency to develop more studies, there is limited information in the southern hemisphere. This study aimed to analyze the levels and profiles of natural (MeO-PBDEs) and anthropogenic (BFRs: PBDEs, HBB, PBEB) organic brominated compounds in adipose tissue of two species of dolphins with different distribution and trophic requirements from the Southwestern Atlantic Ocean; the short-beaked common dolphin (Delphinus delphis) and the Fraser's dolphin (Lagenodelphis hosei). In addition, we aim to investigate maternal transfer and biological pattern relationship (sex, age, sexual maturity) in short-beaked common dolphin bioaccumulation. The levels of both groups of contaminants were in the same order of magnitude as those reported for other marine mammals on both a regional and global scale. BFRs profiles were dominated by BDE 28 and BDE 47 in short-beaked common dolphin and Fraser's dolphin, respectively, whereas 2-MeO-BDE 68 was the most abundant natural compound in both species. Evidence of maternal transfer, temporary increase in BDE 154 levels and no influence of sex, age, or sexual maturity on brominated compound concentration was observed in short-beaked common dolphin. This study fills a gap in the knowledge of the Southwestern Atlantic Ocean providing new information on emerging organic pollutants bioavailability for dolphins and, therefore, for the different trophic webs. In addition, it serves as a baseline for further contamination assessments.

Computational Risk Assessment of Persistence, Bioaccumulation, and Toxicity of Novel Flame-Retardant Chemicals

Novel brominated flame retardants (NBFRs) have emerged as chemicals of environmental concern, as they have been widely used as an alternative to polybrominated diphenyl ethers (PBDEs). Considering the similar structural features of NBFRs and PBDEs necessitates a comprehensive investigation to understand the physicochemical relationships of these compounds and their ability to alter biological functions. In this study, we investigated the persistent nature of NBFRs in terms of thyroid-disrupting potential by understanding the structure-stability aspects using density functional theory (DFT)-based reactivity parameters and interactions via molecular docking and molecular dynamics (MD) simulations. The results indicate that the DFT-based stability descriptor (chemical hardness) is associated with the persistent nature of NBFRs. The computed molecular interaction profile revealed prominent interactions between thyroid receptor-β (TR-β) and NBFRs. Stable trajectory and interactions with TR-β were obtained with ATE, p-TBX, PBT, PBEB, and TBBPA-DBPE during 100 ns of MD simulation. The results of these studies have suggested that the presence of a higher number of halogenated atoms increases the stability vis-à-vis the persistence and endocrine disruption potential of NBFRs.

Bisphenol analogues induce thyroid dysfunction via the disruption of the thyroid hormone synthesis pathway

Bisphenol analogues are widely used in industrial and daily-use consumer products having imperfect thyroid hormones (THs) structures. Widespread exposure interferes with thyroid-related health outcomes in human. The mechanisms of disruption on TH synthesis and subsequent thyroid dysfunction by different bisphenol analogues remain unclear. Here, we evaluated bisphenol-induced thyroid endocrine disruption in C57BL/6 mice at doses of 0.002, 0.02, 2, and 20 mg/kg body weight/day (BW/d) for five consecutive weeks. Administration of 20 mg/kg BW/d bisphenol S (BPS) and 2 mg/kg BW/d tetrabromobisphenol S (TBBPS) significantly increased serum thyrotropin (TSH) levels to 1.21-fold and 1.20-fold of control group, respectively, indicating that bisphenols induced thyroid dysfunction in mice. Height of the thyroid follicle epithelium significantly increased to 1.27-, 1.24-, 1.26-, and 1.36-fold compared to control group with BPA, BPS, TBBPA, and TBBPS at 20 mg/kg BW/d, respectively, indicating impairment of the thyroid gland structure, and TBBPS showed potent effect. Exposure to bisphenol analogues of 0.02 mg/kg BW/d downregulated the protein expression levels of thyrotropin receptor, the sodium/iodide symporter, thyroperoxidase. The TH-dependent effects were further determined using the T-Screen assay at 10-11 M to 10-5 M concentrations. Bisphenol analogues significantly decreased TH-dependent GH3 cell proliferation, indicating the antagonistic activity of bisphenol analogues. The gene responsible for THs synthesis of thyrotropin releasing hormone receptor and TSH were upregulated, but downregulation of thyroid receptor β was observed. Our results suggest that bisphenol analogues distinctly induce thyroid dysfunction via TH synthesis, implying adverse effect of bisphenol analogues on TH homeostasis and subsequent physiological processes.

Alternative and legacy flame retardants in marine mammals from three northern ocean regions

Flame retardants are globally distributed contaminants that have been linked to negative health effects in humans and wildlife. As top predators, marine mammals bioaccumulate flame retardants and other contaminants in their tissues which is one of many human-imposed factors threatening population health. While some flame retardants, such as the polybrominated diphenyl ethers (PBDE), have been banned because of known toxicity and environmental persistence, limited data exist on the presence and distribution of current-use alternative flame retardants in marine mammals from many industrialized and remote regions of the world. Therefore, this study measured 44 legacy and alternative flame retardants in nine marine mammal species from three ocean regions: the Northwest Atlantic, the Arctic, and the Baltic allowing for regional, species, age, body condition, temporal, and tissue comparisons to help understand global patterns. PBDE concentrations were 100-1000 times higher than the alternative brominated flame retardants (altBFRs) and Dechloranes. 2,2',4,5,5'-pentabromobiphenyl (BB-101) and hexabromobenzene (HBBZ) were the predominant altBFRs, while Dechlorane-602 was the predominant Dechlorane. This manuscript also reports only the second detection of hexachlorocyclopentadienyl-dibromocyclooctane (HCDBCO) in marine mammals. The NW Atlantic had the highest PBDE concentrations followed by the Baltic and Arctic which reflects greater historical use of PBDEs in North America compared to Europe and greater industrialization of North America and Baltic countries compared to the Arctic. Regional patterns for other compounds were more complicated, and there were significant interactions among species, regions, body condition and age class. Lipid-normalized PBDE concentrations in harbor seal liver and blubber were similar, but HBBZ and many Dechloranes had higher concentrations in liver, indicating factors other than lipid dynamics affect the distribution of these compounds. The health implications of contamination by this mixture of compounds are of concern and require further research.

The mechanism underlying pentabromoethylbenzene-induced adipogenesis and the obesogenic outcome in both cell and mouse model

Convergent evidence links traditional brominated flame retardants (BFRs) exposure to weight gain, while the obesogenic potency of new BFRs (NBFRs) remain largely unknown. Aiding by luciferase-reporter gene assay, the present study revealed only pentabromoethylbenzene (PBEB), an alternative for penta-BDEs, binds with retinoid X receptor α (RXRα) but not peroxisomeproliferator receptor γ (PPARγ) among the seven testing NBFRs. An apparent induction of adipogenesis in 3T3-L1 cells was observed at nanomolar of PBEB, much lower than penta-BFRs. Mechanistic research uncovered PBEB initiated the adipogenesis by demethylated CpG sites in the PPARγ promoter region. Specifically, activation RXRα by PBEB strengthened the activity of RXRα/PPARγ heterodimer, tightened the interaction between the heterodimer and PPAR response elements, and further enhanced adipogenesis. RNA sequencing combined with k-means clustering analysis exposed adenosine 5'-monophosphate (AMP)-activated protein kinase and phosphoinositide-3-kinase (PI3K)/protein kinase B (AKT) signaling as two predominant pathways that enriched in PBEB-induced lipogenesis. The obesogenic outcome was further corroborated in offspring mice when the maternal mice exposed to environmental relevant doses of PBEB. We found the male offspring exhibited adipocyte hypertrophy and increased weight gain in the epididymal white adipose tissue (eWAT). Consistent with in vitro findings, the reduction in protein phosphorylation of both AMPK and PI3K/AKT were observed within eWAT. Thus, we posited PBEB disrupts the pathways controlling adipogenesis and adipose tissue maintenance, supporting its potential as an environmental obesogen.

TRIAC disrupts cerebral thyroid hormone action via negative feedback and heterogenous distribution among organs

As 3,3',5-triiodothyroacetic acid (TRIAC), a metabolite of thyroid hormones (THs), was previously detected in sewage effluent, we aimed to investigate exogenous TRIAC's potential for endocrine disruption. We administered either TRIAC or 3,3',5-triiodo-L-thyronine (LT3) to euthyroid mice and 6-propyl-2-thiouracil-induced hypothyroid mice. In hypothyroid mice, TRIAC administration suppressed the hypothalamus-pituitary-thyroid (HPT) axis and upregulated TH-responsive genes in the pituitary gland, the liver, and the heart. We observed that, unlike LT3, TRIAC administration did not upregulate cerebral TH-responsive genes. Measurement of TRIAC contents suggested that TRIAC was not efficiently trafficked into the cerebrum. By analyzing euthyroid mice, we found that cerebral TRIAC content did not increase despite TRIAC administration at higher concentrations, whereas serum levels and cerebral contents of THs were substantially decreased. Disruption by TRIAC is due to the additive effects of circulating endogenous THs being depleted via a negative feedback loop involving the HPT axis and heterogeneous distribution of TRIAC among different organs.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Machine Learning Model for Screening Thyroid Stimulating Hormone Receptor Agonists Based on Updated Datasets and Improved Applicability Domain Metrics

Machine learning (ML) models for screening endocrine-disrupting chemicals (EDCs), such as thyroid stimulating hormone receptor (TSHR) agonists, are essential for sound management of chemicals. Previous models for screening TSHR agonists were built on imbalanced datasets and lacked applicability domain (AD) characterization essential for regulatory application. Herein, an updated TSHR agonist dataset was built, for which the ratio of active to inactive compounds greatly increased to 1:2.6, and chemical spaces of structure-activity landscapes (SALs) were enhanced. Resulting models based on 7 molecular representations and 4 ML algorithms were proven to outperform previous ones. Weighted similarity density (ρ_s) and weighted inconsistency of activities (I_A) were proposed to characterize the SALs, and a state-of-the-art AD characterization methodology AD_SAL{ρ_s, I_A} was established. An optimal classifier developed with PubChem fingerprints and the random forest algorithm, coupled with AD_SAL{ρ_s ≥ 0.15, I_A ≤ 0.65}, exhibited good performance on the validation set with the area under the receiver operating characteristic curve being 0.984 and balanced accuracy being 0.941 and identified 90 TSHR agonist classes that could not be found previously. The classifier together with the AD_SAL{ρ_s, I_A} may serve as efficient tools for screening EDCs, and the AD characterization methodology may be applied to other ML models.

Investigating the adverse outcome pathways (AOP) of neurotoxicity induced by DBDPE with a combination of in vitro and in silico approaches

Current studies have shown an association between DBDPE and neurotoxicity. In this study, the adverse outcome pathway (AOP) and mechanistic analysis of DBDPE-induced neurotoxicity were explored by a combination of in vitro and in silico approaches in SK-N-SH cells. DBDPE-induced oxidative stress caused DNA strand breaks, resulting in the activation of poly (ADP-ribose) (PAR) polymerase-1 (PARP-1). Activation of PARP1 could cause toxic damage in various organ systems, especially in the nervous system. DBDPE-induced apoptosis via the caspase-dependent intrinsic mitochondrial pathway and the PARP1-dependent pathway. Activation of PARP1 by DBDPE was deemed the initiating event, thereby affecting the key downstream biochemical events (e.g., ROS production, DNA damage, membrane potential changes, and ATP reduction), which induced apoptosis. Furthermore, excessive activation of PARP1 was accompanied by the translocation of the apoptosis-inducing factor (AIF), which was associated with PARP1-dependent cell death. The inhibition of PARP1 by PJ34 reduced DBDPE-induced apoptosis and maintained cellular ATP levels. PJ34 also prevented the translocation of AIF from the mitochondria to the nucleus. These findings improve the understanding of the mechanism of DBDPE-induced neurotoxic effects and provide a theoretical basis for the ecological risk of DBDPE.

Potential Disrupting Effects of Wastewater-Derived Disinfection Byproducts on Chinese Rare Minnow (Gobiocypris rarus) Transthyretin: An In Vitro and In Silico Study

The available information about whether wastewater-derived disinfection byproducts (DBPs) could elicit potential endocrine-related detrimental effects on aquatic organisms was scarce. Herein, the potential disrupting effects and underlying binding mechanism of 14 wastewater-derived aliphatic and aromatic DBPs and 12 other substances on Chinese rare minnow (Gobiocypris rarus) transthyretin (CrmTTR) were tested and revealed by in vitro and in silico methods. The amino acid sequences of CrmTTR were determined, and the recombinant CrmTTR with a molecular mass of 66.3 kDa was expressed and purified. In vitro assay results indicated that eight selected aromatic DBPs exhibited detectable CrmTTR disrupting ability. Meanwhile, six aliphatic DBPs were not CrmTTR binders. Molecular modeling results implied that hydrophobic hydrogen bonds and/or ionic pair interactions were non-negligible. Four binary classification models with high classification performance were constructed. A significant positive linear relationship was observed for the binding affinity data from CrmTTR and human TTR (n = 18, r = 0.922, p < 0.0001). However, the binding affinity for 13 out of 18 tested compounds with CrmTTR was higher than that with human TTR. All the results highlighted that some wastewater-derived DBPs may be potential disruptors on the aquatic organism endocrine system, and interspecies variation should not be neglected in future determination of the potential endocrine disrupting effects of wastewater-derived DBPs.

Perfluorooctanoic acid and perfluorooctanesulfonic acid induce immunotoxicity through the NF-κB pathway in black-spotted frog (Rana nigromaculata)

Perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) are widely detected in the environment and wild animals, thus posing a threat to wildlife and public health; however, knowledge about their immunotoxicity and the underlying mechanism remains limited. In the present study, male black-spotted frogs (Rana nigromaculata) were exposed to environmentally relevant concentrations (0, 1, and 10 μg/L) of PFOA or PFOS for 21 days; subsequently, biochemical analysis, molecular docking, and gene expression determination were conducted. The results indicated that exposure to 10 μg/L PFOA decreased the serum levels of immunoglobulin A. PFOS exposure significantly increased the hepatic levels of interleukin-1β, interleukin-6, tumor necrosis factor-α, interferon-γ, and nitric oxide; but PFOA significantly increased the levels of only tumor necrosis factor-α. Furthermore, PFOA and PFOS exposure significantly decreased the activity of inducible nitric oxide synthase and total nitric oxide synthase. IBRv2 analysis indicated that PFOA and PFOS had a similar effect on these immune indicators, but PFOS was more toxic than PFOA. Molecular docking revealed that PFOA and PFOS can bind to nuclear factor-κB (NF-κB) by forming stable hydrogen bonds. PFOA and PFOS exposure upregulated the gene expression of NF-κB and its downstream genes. Significant correlations between the expression of genes involved in the NF-κB pathway and immune-related indicators suggests that PFOA- and PFOS-induced immunotoxicity was associated with the activation of NF-κB. Our findings provide novel insights into the potential role of NF-κB in immunotoxicity induced by PFOA and PFOS in frogs.

Promotion of Bladder Cancer Cell Metastasis by 2-Mercaptobenzothiazole via Its Activation of Aryl Hydrocarbon Receptor Transcription: Molecular Dynamics Simulations, Cell-Based Assays, and Machine Learning-Driven Prediction

2-Mercaptobenzothiazole (MBT) is an industrial chemical widely used for rubber products, corrosion inhibitors, and polymer materials with multiple environmental and exposure pathways. A growing body of evidence suggests its potential bladder cancer (BC) risk as a public health concern; however, the molecular mechanism remains poorly understood. Herein, we demonstrate the activation of the aryl hydrocarbon receptor (AhR) by MBT and reveal key events in carcinogenesis associated with BC. MBT alters conformational changes of AhR ligand binding domain (LBD) as revealed by 500 ns molecular dynamics simulations and activates AhR transcription with upregulation of AhR-target genes CYP1A1 and CYP1B1 to approximately 1.5-fold. MBT upregulates the expression of MMP1, the cancer cell metastasis biomarker, to 3.2-fold and promotes BC cell invasion through an AhR-mediated manner. MBT is further revealed to induce differentially expressed genes (DEGs) most enriched in cancer pathways by transcriptome profiling. The exposure of MBT at environmentally relevant concentrations induces BC risk via AhR signaling disruption, transcriptome aberration, and malignant cell metastasis. A machine learning-based model with an AUC value of 0.881 is constructed to successfully predict 31 MBT analogues. Overall, we provide molecular insight into the BC risk of MBT and develop an effective tool for rapid screening of AhR agonists.

Azole-Induced Color Vision Deficiency Associated with Thyroid Hormone Signaling: An Integrated In Vivo, In Vitro, and In Silico Study

Azoles that are used in pesticides, pharmaceuticals, and personal care products can have toxic effects on fish. However, there is no information regarding azole-induced visual disorder associated with thyroid disruption. We evaluated changes in retinal morphology, optokinetic response, transcript abundance of the genes involved in color perception and hypothalamic-pituitary-thyroid (HPT) axis, and thyroid hormone (TH) levels in zebrafish larvae exposed to common azoles, such as climbazole (CBZ, 0.1 and 10 μg/L) and triadimefon (TDF, 50 and 500 μg/L), at environmentally relevant and predicted worst-case environmental concentrations. Subsequently, the effect of azoles on TH-dependent GH3 cell proliferation and thyroid receptor (TR)-regulated transcriptional activity, as well as the in silico binding affinity between azoles and TR isoforms, was investigated. Azole exposure decreased cell densities of the ganglion cell layer, inner nuclear layer, and photoreceptor layer. Zebrafish larvae exposed to environmentally relevant concentrations of CBZ and TDF showed a decrease in optokinetic response to green-white and red-white stripes but not blue-white stripes, consistent with disturbance in the corresponding opsin gene expression. Azole exposure also reduced triiodothyronine levels and concomitantly increased HPT-related gene expression. Molecular docking analysis combined with in vitro TR-mediated transactivation and dual-luciferase reporter assays demonstrated that CBZ and TDF exhibited TR antagonism. These results are comparable to those obtained from a known TR antagonist, namely, TR antagonist 1, as a positive control. Therefore, damage to specific color perception by azoles appears to result from lowered TH signaling, indicating the potential threat of environmental TH disruptors to the visual function of fish.

Concentration-dependent toxicokinetics of novel PFOS alternatives and their chronic combined toxicity in adult zebrafish

The increasing use of perfluorooctanesulfonate (PFOS) alternatives has led to their release into the aquatic environment. This study sought to determine the effects of exposure concentration on the toxicokinetics of PFOS and its alternatives, including perfluorobutanesulfonic acid (PFBS), perfluorohexanesulfonic acid (PFHxS), chlorinated polyfluorinated ether sulfonate (F-53B) and sodium p-perfluorous nonenoxybenzenesulfonate (OBS) in adult zebrafish by exposure to mixtures of the five per- and polyfluoroalkyl substances (PFAS) at 1, 10, and 100 ng/mL for 28-day, followed by a 14-day depuration. PFAS predominantly accumulated in the blood and liver, and the bioconcentration factor (BCF) decreased in the order of F-53B > PFOS > OBS ≫ PFHxS > PFBS in whole-fish homogenates. The uptake rate constants and BCF of the short-chain PFAS (≤C6) positively correlated with increasing exposure concentration, while the long-chain PFAS (≥C8) exhibited a pattern of first increasing and then decreasing. A consistent increase in the elimination rate constants of short- and long-chain PFAS was observed with increasing exposure concentration. All PFAS form tight conformations with ZSA and ZL-FABP via hydrogen bonding as revealed by molecular docking analysis. Furthermore, chronic combined exposure to PFAS induced the occurrence of vacuolation and oxidative stress in the zebrafish liver. Our findings uniquely inform the concentration-dependent bioconcentration potential and health risks to aquatic organisms of these PFOS alternatives in the environment.

Thyroid endocrine disruption and hepatotoxicity induced by bisphenol AF: Integrated zebrafish embryotoxicity test and deep learning

Bisphenol AF (BPAF) is an emerging contaminant prevalent in the environment as one of main substitutes of bisphenol A (BPA). It was found that BPAF exhibited estrogenic effects in zebrafish larvae in our previous study, while little is known about its effects on the thyroid and liver. A 7 d zebrafish embryotoxicity test was conducted to study the potential thyroid disruption and hepatotoxicity of BPAF. BPAF decreased levels of thyroid hormones and deiodinases but increased expressions of transthyretin at 12.5 and 125 μg/L after 7 d exposure, indicating that both the metabolism and transport of thyroid hormones were perturbed. The thyroid hormone receptor (TR) levels decreased significantly upon exposure to ≥12.5 μg/L BPAF, implying that BPAF acts as a TR antagonist, which coincided well with the prediction from the Direct Message Passing Neural Network. The liver impairment (mainly cell necrosis of hepatocytes) and apoptosis were triggered by 125 μg/L and ≥12.5 μg/L BPAF respectively, accompanied by the increased activities of caspase 3 and caspase 9. Thus BPAF might not be a safe alternative to BPA given the thyroid and liver toxicity. DMPNN appears useful to screen for thyroid disrupting activity from molecular structures.

Nonmetal doped carbon nitride nanosheet as photocatalyst for degradation of 4, 5-dichloroguaiacol

Metal-free photocatalysts for high efficient photocatalytic degradation of pollutants have attracted growing concern in recent years. Herein, relying on density functional theory (DFT) calculations, boron and phosphorus doped C₂N layers were explored for the potential of utilization as photocatalysts for 4, 5-dichloroguaiacol (4, 5-DCG) removal. Our computations revealed that the adsorption energy of 4, 5-DCG on B@N-doped C₂N layers were 26.56 kcal mol^-1, and the ΔG^≠ of initial reactions of 4, 5-DCG with OH were also reduced onto the B@N-doped C₂N substrates. The band gap of B@N-doped C₂N was 2.27 eV. The obtained results showed that the doping of boron atom into C₂N layer narrows bandgap, and retains well catalytic performance and adsorption properties. Hence, B@N-doped C₂N layer is a promising photocatalyst for organic pollutants removal. Possible degradation pathways of 4, 5-DCG and aquatic toxicity assessment during degradation were also carried out. Products with higher toxicity would be formed and the transformation products were still toxic to three nutrient levels of aquatic organisms (green algae, fish, and daphnia).

Computational biotransformation of polyethylene terephthalate by depolymerase: A QM/MM approach

Despite increasing environmental concerns on ever-lasting Polyethylene Terephthalate (PET), its global production is continuously growing. Effective strategies that can completely remove PET from environment are urgently desired. Here biotransformation processes of PET by one of the most effective enzymes, leaf-branch compost cutinase (LCC), were systematically explored with Molecular Dynamics and Quantum Mechanics/Molecular Mechanics approaches. We found that four concerted steps are required to complete the whole catalytic cycle. The last concerted step, deacylation, was determined as the rate-determining step with Boltzmann-weighted average barrier of 13.6 kcal/mol and arithmetic average of 16.1 ± 2.9 kcal/mol. Interestingly, unprecedented fluctuations of hydrogen bond length during LCC catalyzed transformation process toward PET were found. This fluctuation was also observed in enzyme IsPETase, indicating that it may widely exist in other catalytic triad (Ser-His-Asp) containing enzymes as well. In addition, possible features (bond, angle, dihedral angle and charge) that influence the catalytic reaction were identified and correlations between activation energies and key features were established. Our results present new insights into catalytic mechanism of hydrolases and shed light on the efficient recycling of the ever-lasting PET.

The interactive effects of ethinylestradiol and progesterone on transcriptional expression of genes along the hypothalamus-pituitary-thyroid axis in embryonic zebrafish (Danio rerio)

Progestins and estrogens are widespread in various aquatic environments and their potential endocrine disruption effects to aquatic organisms have drawn growing concern. However, their combined effects in aquatic organisms remain elusive. The aim of the present study was to assess the effects of the binary mixtures of gestodene (GES) and 17α-ethinylestradiol (EE2) on the hypothalamic-pituitary-thyroid (HPT) axis of zebrafish (Danio rerio) using the eleuthero-embryos. Embryos were exposed to GES and EE2 alone or in combination at concentrations ranging from 41 to 5329 ng L^-1 (nominal ones from 50 to 5000 ng L^-1) for 48 h, 96 h and 144 h post fertilization (hpf). The results showed that the transcripts of the genes along the HPT axis displayed pronounced alterations. There was no clear pattern in the change of the transcripts of these genes over time and with concentrations. However, in general, the transcripts of the genes were inversely affected by EE2 (increase 0.5 to 4.2-folds) and GES (inhibition 0.4 to 4.9-folds), and their mixtures showed interactive effects in embryonic zebrafish. In addition, physiological data (mortality, malformation, body length and heart rate etc.) denoted higher toxicity of the two chemicals in combination than alone based on the developmental toxicity and neurotoxicity (locomotor behavior). These results indicated that the interactive effects of these two chemicals might be different between at the transcriptional level and at the whole organismal level. In summary, GES and EE2 affect the HPT axis (related genes expression and thyroid hormones (THs) levels) and exhibit developmental toxicity and neurotoxicity.

PFOA regulate adenosine receptors and downstream concentration-response cAMP-PKA pathway revealed by integrated omics and molecular dynamics analyses

As an important pollutant, perfluorooctane acid (PFOA) has been widely concerned and reported by thousands of times, while less is known about the concentration-response pathway of PFOA. The aim of the present work was to reveal the concentration-response mechanism of PFOA in human cells. Omics results showed that calcium-related pathways play key roles in PFOA injury mechanisms. The results of GO and KEGG analyses showed that the cAMP signaling pathway was presented as the top one in all of the regulatory patterns and concentrations groups of PFOA. In the cAMP signaling pathway, the adenosine A₁ receptor (ADORA1) recognized the low concentration of PFOA and induced pathway "G_i-cAMP-PKA" to decrease the concentration of cAMP. This indicated that the low concentration of PFOA may promote breast hyperplasia and inhibit lactation. While adenosine A_2A receptor (ADORA2A) recognized the high concentration of PFOA and induced pathway "G_S-AC-cAMP-RKA" to increase the concentration of cAMP, induce cell damage and may lead to the deterioration of breast cancer. The results of molecular dynamics simulation showed that PFOA could bind to ADORA1 and ADORA2A, thus cause subsequent signal transduction. Furthermore, considering the strong binding ability of PFOA with ADORA1, PFOA tends to bind to ADORA1 at a low concentration. On the other side, PFOA at high concentration will continue to bind to another receptor protein, ADORA2A, and activate subsequent signaling pathways. Combined analyses of transcriptomic and proteomic revealed that different concentrations of PFOA regulate cellular calcium-related pathways. The cAMP pathway showed a concentration-response effect of PFOA. After treatment with different concentrations of PFOA, ADORA1 and ADORA2A were activated respectively, showing opposite cellular effects, leading to kinds of breast lesions. In the nervous system, PFOA might induce a variety of nervous system diseases. The present work was an exploration on the toxicological mechanism of PFOA, providing important information on the health impacts of PFOA in humans.

Carcinogenic Risk of 2,6-Di-tert-Butylphenol and Its Quinone Metabolite 2,6-DTBQ Through Their Interruption of RARβ: In Vivo, In Vitro, and In Silico Investigations

Thousands of contaminants are used worldwide and eventually released into the environment, presenting a challenge of health risk assessment. The identification of key toxic pathways and characterization of interactions with target biomacromolecules are essential for health risk assessments. The adverse outcome pathway (AOP) incorporates toxic mechanisms into health risk assessment by emphasizing the relationship among molecular initiating events (MIEs), key events (KEs), and adverse outcome (AO). Herein, we attempted the use of AOP to decipher the toxic effects of 2,6-di-tert-butylphenol (2,6-DTBP) and its para-quinone metabolite 2,6-di-tert-butyl-1,4-benzoquinone (2,6-DTBQ) based on integrated transcriptomics, molecular modeling, and cell-based assays. Through transcriptomics and quantitative real-time PCR validation, we identified retinoic acid receptor β (RARβ) as the key target biomacromolecule. The epigenetic analysis and molecular modeling revealed RARβ interference as one MIE, including DNA methylation and conformational changes. In vitro assays extended subsequent KEs, including altered protein expression of p-Erk1/2 and COX-2, and promoted cancer cell H4IIE proliferation and metastasis. These toxic effects altogether led to carcinogenic risk as the AO of 2,6-DTBP and 2,6-DTBQ, in line with chemical carcinogenesis identified from transcriptome profiling. Overall, our simplified AOP network of 2,6-DTBP and 2,6-DTBQ facilitates relevant health risk assessment.

Toxicity pathways of lipid metabolic disorders induced by typical replacement flame retardants via data-driven analysis, in silico and in vitro approaches

Endocrine-disrupting chemicals can interfere with hormone action via various pathways, thereby increasing the risk of adverse health outcomes. Organophosphorus ester (OPEs) retardants, a group of new emerging endocrine disruption chemicals, have been referred to as metabolism disruptors and reported to induce chronic health problems. However, the toxicity pathways were mainly focused on nuclear receptor signaling pathways. Significantly, the membrane receptor pathway (such as G protein-coupled estrogen receptor 1 (GPER) signaling pathway) had been gradually realized as the important role in respond more effective to lipid metabolism disorder than traditional nuclear receptors, whereas the detailed mechanism was unclear yet. Therefore, this study innovatively integrated the bibliometric analysis, in silico and in vitro approach to develop toxicity pathways for the mechanism interpretation. Bibliometric analysis found that the typical OPEs - triphenyl phosphate was a major concern of lipid metabolism abnormality. Results verified that TPP could damage the structures of cell membranes and exert an agonistic effect of GPER as the molecular initiating event. Then, the activated GPER could trigger the PI3K-Akt/NCOR1 and mTOR/S6K2/PPARα transduction pathways as key event 1 (KE1) and affect the process of lipid metabolism and synthesis (CPT1A, CPT2, SREBF2 and SCD) as KE2. As a result, these alterations led to lipid accumulation as adverse effect at cellular-levels. Furthermore, the potential outcomes (such as immunity damage, weight change and steatohepatitis) at high biological levels were expanded. These findings improved knowledge to deeply understand toxicity pathways of phosphorus flame retardants and then provided a theoretical basis for risk assessments.

Benzophenone-1 induced aberrant proliferation and metastasis of ovarian cancer cells via activated ERα and Wnt/β-catenin signaling pathways

Benzophenone-1 (BP-1) belongs to personal care product-related contaminants of emerging concern and has been recently reported to induce xenoestrogenic effects. However, the underlying mechanisms leading to the activation of target receptors and subsequent various adverse outcomes remain unclear, which is beneficial to safety and health risk assessment of benzophenone-type ultraviolet filters with their widespread occurrence. Herein, we investigated disrupting effects of BP-1 at environmentally relevant concentrations (10^-9-10^-6 M) on estrogen receptor (ER) α-associated signaling pathways. Molecular dynamics simulations together with yeast-based assays revealed the steady binding of BP-1 to ERα ligand binding domain (LBD) and hence the observed agonistic activity. BP-1 triggered interaction between ERα and β-catenin in human SKOV3 ovarian cancer cells and caused translocation of β-catenin from the cytoplasm to the nucleus, leading to aberrant activation of Wnt/β-catenin. BP-1 consequently induced dissemination of SKOV3 via regulating epithelial-mesenchymal transitions (EMT) biomarkers including minimally downregulating ZO-1 gene to 78.0 ± 10.1% and maximally upregulating MMP9 gene to 144.1 ± 29.7% and promoted 1.03-1.83 fold proliferation, migration and invasion of SKOV3. We provide the first evidence that the BP-1 activated ERα triggers crosstalk between ERα and Wnt/β-catenin pathway, leading to the abnormal stimulation and progression of SKOV3 cancer cells.

Characterizing the potentially neuronal acetylcholinesterase reactivity toward chiral pyraclofos: Enantioselective insights from spectroscopy, in silico docking, molecular dynamics simulation and per-residue energy decomposition studies

Chiral organophosphorus agents are distributed ubiquitously in the environment, but the neuroactivity of these asymmetric chemicals to humans remains uncertain. This scenario was to explore the stereoselective neurobiological response of human acetylcholinesterase (AChE) to chiral pyraclofos at the enantiomeric scale, and then decipher the microscopic basis of enantioselective neurotoxicity of pyraclofos enantiomers. The results indicated that (R)-/(S)-pyraclofos can form the bioconjugates with AChE with a stoichiometric ratio of 1:1, but the neuronal affinity of (R)-pyraclofos (K = 6.31 × 10⁴ M^-1) with AChE was larger than that of (S)-pyraclofos (K = 1.86 × 10⁴ M^-1), and significant enantioselectivity was existed in the biochemical reaction. The modes of neurobiological action revealed that pyraclofos enantiomers were situated at the substrate binding domain, and the strength of the overall noncovalent bonds between (S)-pyraclofos and the residues was weaker than that of (R)-pyraclofos, resulting in the high inhibitory effect of (R)-pyraclofos toward the activity of AChE. Dynamic enantioselective biointeractions illustrated that the intervention of inherent conformational flexibility in the AChE-(R)-pyraclofos was greater than that of the AChE-(S)-pyraclofos, which arises from the big spatial displacement and the conformational flip of the binding domain composed of the residues Thr-64～Asn-89, Gly-122～Asp-134, and Thr-436～Tyr-449. Energy decomposition exhibited that the Gibbs free energies of the AChE-(R)-/(S)-pyraclofos were ΔG° = －37.4/－30.2 kJ mol^-1, respectively, and the disparity comes from the electrostatic energy during the stereoselective neurochemical reactions. Quantitative conformational analysis further confirmed the atomic-scale computational chemistry conclusions, and the perturbation of (S)-pyraclofos on the AChE's ordered conformation was lower than that of (R)-pyraclofos, which is germane to the interaction energies of the crucial residues, e.g. Tyr-124, Tyr-337, Asp-74, Trp-86, and Tyr-119. Evidently, this attempt will contribute mechanistic information to uncovering the neurobiological effects of chiral organophosphates on the body.

Human transthyretin binding affinity of halogenated thiophenols and halogenated phenols: An in vitro and in silico study

Serious harmful effects have been reported for thiophenols, which are widely used industrial materials. To date, little information is available on whether such chemicals can elicit endocrine-related detrimental effects. Herein the potential binding affinity and underlying mechanism of action between human transthyretin (hTTR) and seven halogenated-thiophenols were examined experimentally and computationally. Experimental results indicated that the halogenated-thiophenols, except for pentafluorothiophenol, were powerful hTTR binders. The differentiated hTTR binding affinity of halogenated-thiophenols and halogenated-phenols were observed. The hTTR binding affinity of mono- and di-halo-thiophenols was higher than that of corresponding phenols; while the opposite relationship was observed for tri- and penta-halo-thiophenols and phenols. Our results also confirmed that the binding interactions were influenced by the degree of ligand dissociation. Molecular modeling results implied that the dominant noncovalent interactions in the molecular recognition processes between hTTR and halogenated-thiophenols were ionic pair, hydrogen bonds and hydrophobic interactions. Finally, a model with acceptable predictive ability was developed, which can be used to computationally predict the potential hTTR binding affinity of other halogenated-thiophenols and phenols. Taken together, our results highlighted that more research is needed to determine their potential endocrine-related harmful effects and appropriate management actions should be taken to promote their sustainable use.

Transcriptome aberration in mice uterus associated with steroid hormone response and inflammation induced by dioxybenzone and its metabolites

Benzophenone-type UV filters have been implicated in multiple adverse reproductive outcomes, yet the underlying processes and molecular targets on the female reproductive tract remain largely unknown. Herein, we investigated the effect of dioxybenzone, one of the widely used congeners, and its demethylated (M1) and hydroxylated (M2) metabolites on transcriptome profiles of ICR mice uterus and identified potential cellular targets in human endometrial stromal cells (HESCs) separated from normal endometrium tissues. Dioxybenzone, M1 and M2 (20 mg/kg bw/d) significantly induced transcriptome aberration with the induction of 683, 802, and 878 differentially expressed genes mainly involved in cancer, reproductive system disease and inflammatory disease. Compared to dioxybenzone, M1 and M2 exhibited a transcriptome profile more similar to estradiol in mice uterus, and subsequently promoted thicker endometrial columnar epithelial layer through upregulation of estrogen receptor target genes-Sprr2s. Dioxybenzone, M1 and M2 (0.1 or 1 μM) also exhibited estrogenic disrupting effect via increasing the mRNA expressions and production of the growth factors responsible for epithelial proliferation, including Fgfs and Igf-1 in HESCs. Additionally, the mRNA expressions of several inflammatory cytokines especially IL-1β in mice uterus and HESCs was significantly upregulated by dioxybenzone and its metabolites. Overall, we revealed that dioxybenzone and its metabolites triggered transcriptome perturbation dually associated with abnormal steroid hormone response and inflammation, both as key determinants to reproductive health risks.

Severe contamination and time trends of legacy and novel halogenated flame retardants in multiple environmental media from Lake Shihwa, Korea: Effectiveness of regulatory action

Novel halogenated flame retardants (HFRs) were introduced to industrial markets as alternatives to legacy brominated FRs (BFRs), such as polybrominated diphenyl ethers (PBDEs). In the present study, PBDEs and their brominated and chlorinated alternatives, novel BFRs (NBFRs) and dechlorane plus (DP), were measured in multiple environmental matrices in a highly industrialized lake in Korea. Legacy and novel HFRs were detected in multiple samples, indicating ubiquitous contamination. Concentrations of HFRs in water and sediment observed in creeks running through machine, textiles, and automobile industrial complexes were significantly higher than those observed in inside and outside of the lake. Higher bioaccumulation levels of HFRs were observed in inshore compared with offshore waters. Results suggest that multi-matrix distribution of legacy and novel HFRs was dependent on the geographical proximity to industrial sources. Compared with previous studies, the highest levels of PBDEs and NBFRs were recorded in water samples on a global scale, implying on-going emissions from industrial activities. Decabromodiphenyl ethane (DBDPE) was a dominant compound in water samples, whereas the concentrations of PBDEs, NBFRs, and DP in sediment were similar to each other. This suggests a shift in consumption from legacy to novel HFRs, preferentially in water environments. A significant declining trend in PBDEs was observed in water and sediment collected between 2008 and 2015, indicating the effectiveness of regulatory actions. Based on their environmental occurrence and bioaccumulation potential, pentabromoethylbenzene and bis(2-ethylhexyl) tetrabromophthalate may pose emerging concerns regarding contamination of aquatic environments.

New insights into the mechanism of hepatocyte apoptosis induced by typical organophosphate ester: An integrated in vitro and in silico approach

Apoptosis is one of the typical features of liver diseases, therefore molecular targets of hepatic apoptosis and regulatory mechanisms need to be further investigated. The caspases play important functions in the execution of apoptosis and many studies have focused on classical caspase-dependent cell death pathways. However, other types of cell death pathways (such as mitochondrial poly (ADP-ribose) polymerase-1 (PARP1) pathway) are suggested to be also as important as the caspase-mediated pathways in reflection of early toxic effects in hepatocytes, which requires additional research. In this work, an approach integrated in silico and in vitro was used to investigate the underlying toxicological mechanisms of hepatocyte apoptosis through the PARP1 dependent cell death pathway induced by triphenyl phosphate (TPP). Docking view showed that TPP could interact with helix αJ to affect the activation of PARP1 as a molecular initial event. In vitro assays suggested some biochemical events downstream of PARP1 activation, such as mitochondrial injury, apoptosis inducing factor (AIF) release, reactive oxygen species (ROS) production, and DNA damage. Moreover, the apoptosis was alleviated when cells were pretreated with PJ34 hydrochloride (PARP1 inhibitor), suggesting the mitochondrial PARP1 dependent pathway played a pivotal role in L02 cells apoptosis. This study indicated that PARP1 was an important molecular target in this process. And it also helped to understand the mechanism of hepatocytes apoptosis, early hepatic toxicity, and even liver diseases.

Obesogenic and developmental effects of TBT on the gene expression of juvenile Japanese medaka (Oryzias latipes)

The widely used antifoulant tributyltin chloride (TBT) is highly toxic to aquatic organisms. In the present study, four-week-old Japanese medaka (Oryzias latipes) juveniles were orally exposed to TBT at 1 and 10 ng/g bw/d for 1, 2, and 4 weeks, respectively. Half of the tested medaka juveniles showed bone morphology alterations in both 1 and 10 ng/g bw/d TBT 4-week exposure groups. Nile Red (NR) staining showed that the juveniles exposed to 1 ng/g bw/d TBT for 2 and 4 weeks had significantly enlarged adipocyte areas. The mRNA-Seq analysis indicated that 1 ng/g bw/d TBT exposure for 2 weeks affected bone morphology through developmental processes. The GO and KEGG analyses suggested that the adipogenic effect of TBT observed in this study may be induced by metabolic processes, oxidative phosphorylation, and fatty acid degradation and metabolism pathways. Therefore, both morphological observation and mRNA-Seq analysis showed obesogenic effects and developmental toxicity of TBT to juvenile Japanese medaka.

2,6-Di-tert-butylphenol and its quinone metabolite trigger aberrant transcriptional responses in C57BL/6 mice liver

2,6-Di-tert-butylphenol (2,6-DTBP) is used as an antioxidant with wide commercial applications and its residues have been detected in various environmental matrices. 2,6-DTBP may enter human body via ingestion, inhalation or other exposure pathways. However, its susceptibility to biotransformation and potential of the metabolic products to trigger aberrant transcriptional responses remain unclear. Here, we investigated in vitro and in vivo biotransformation of 2,6-DTBP and characterized the RNA-Seq based transcriptional profiling of C57BL/6 mice liver after the exposure to 2,6-DTBP and its metabolites. 2,6-DTBP was metabolized into hydroxylated (2,6-DTBH) and para-quinone (2,6-DTBQ) products with residues detected in serum and liver of C57BL/6 mice. 2,6-DTBP and 2,6-DTBQ induced the aberrant transcription in C57BL/6 mice liver featured with 373-2861 differentially expressed genes (DEGs). They also up-regulated 1.09-2.92 fold mRNA expression of carcinogenesis-related genes such as Ccnd1, TGFβ1 and FOS in C57BL/6 mice liver. Our study indicated potential carcinogenic risk of 2,6-DTBP and its metabolites, beneficial to further evaluation of health risk of TBPs-related contaminants.

Effect of short-chain chlorinated paraffins (SCCPs) on lipid membranes: Combination of molecular dynamics and membrane damage experiments

In recent years, more attention has been paid to the biological effects of short-chain chlorinated paraffin (SCCP). Studies have shown that SCCPs exposure could cause metabolic damage and lipid metabolic damage. In the present work, based on E. coli membrane damage experiments and molecular dynamics (MD) simulation, the effects of SCCPs on the membrane structure and membrane properties were studied to explore the possible toxic damage effects of SCCPs on cell membrane. Experiments results showed that SCCPs had a significant inhibitory effect on E. coli. The E. coli cell membrane of the bacteria was broken and the macromolecules of the cell flowed out when exposed to SCCPs. SCCPs would lead to the decrease and depolarization of cell membrane potential, and then affect the integrity and permeability of cell membrane. The further molecular dynamic simulation revealed that SCCP molecules can easily enter the lipid DPPC membranes from the aqueous phase and tended to aggregate inside bilayer stably. The bound of SCCPs could lead to significant variations in DPPC bilayer with a less dense, more disorder and rougher layer, which thus made the damage of cell membrane. In a word, although the overall toxicity of SCCPs to cell was relatively weak, the damage to the cell membrane may be one of the mechanisms of its toxicity. MAIN FINDING OF THE WORK: The exposure of SCCPs could cause structural change of cell membrane in E. coli, which verified the damage to the cell membrane may be one of the mechanisms of its toxicity.

Enhanced Disrupting Effect of Benzophenone-1 Chlorination Byproducts to the Androgen Receptor: Cell-Based Assays and Gaussian Accelerated Molecular Dynamics Simulations

Benzophenone-1 (BP-1), one of the commonly used ultraviolet filters, has caused increasing public concern due to frequently detected residues in environmental and recreational waters. Its susceptibility to residual chlorine and the potential to subsequently trigger endocrine disruption remain unknown. We herein investigated the chlorination of BP-1 in swimming pool water and evaluated the endocrine disruption toward the human androgen receptor (AR). The structures of monochlorinated (P1) and dichlorinated (P2) products were separated and characterized by mass spectrometry and ¹H-¹H NMR correlation spectroscopy. P1 and P2 exhibited significantly higher antiandrogenic activity in yeast two-hybrid assays (EC₅₀, 6.13 μM and 9.30 μM) than did BP-1 (12.89 μM). Our 350 ns Gaussian accelerated molecular dynamics simulations showed the protein dynamics in a long-time scale equilibrium, and further energy calculations revealed that although increased hydrophobic interactions are primarily responsible for enhanced binding affinities between chlorinated products and the AR ligand binding domain, the second chloride in P2 still hinders the complex motion because of the solvation penalty. The mixture of BP-1-P1-P2 elicited additive antiandrogenic activity, well fitted by the concentration addition model. P1 and P2 at 1 μM consequently downregulated the mRNA expression of AR-regulated genes, NKX3.1 and KLK3, by 1.7-9.1-fold in androgen-activated LNCaP cells. Because chlorination of BP-1 occurs naturally by residual chlorine in aquatic environments, our results regarding enhanced antiandrogenic activity and disturbed AR signaling provided evidence linking the use of personal care products with potential health risks.