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

Bacterial N-acyl-homoserine lactone enhances the degradation of sulfamethoxazole by microalgae and the associated metabolic regulatory mechanisms

Recent studies have revealed that N-acyl-homoserine lactones (AHLs), common quorum-sensing (QS) signal molecules in Gram-negative bacteria, can also influence microalgal cells. However, their role in regulating the metabolism of pollutants, such as antibiotics, within microalgae remains poorly understood. This study investigated the effects of N-hexanoyl-L-homoserine lactone (C6-HSL) on the degradation of sulfamethoxazole (SMX) in aquaculture wastewater by Chlorella vulgaris. The addition of 0.5 μM C6-HSL resulted in the highest biomass accumulation and the maximum SMX removal efficiency (95.6 %). At optimal concentrations, C6-HSL effectively modulated key secondary messenger signaling pathways including reactive oxygen species (ROS), nitric oxide (NO), and calcium ions (Ca2+) in microalgal cells. Additionally, it upregulated the activity of detoxification enzymes such as glutathione S-transferase (GST) and cytochrome P450 (CYP450), thereby altering SMX degradation pathways and significantly enhancing its removal. Transcriptomic analysis further demonstrated that exogenous C6-HSL upregulated critical genes associated with ROS, Ca2+, and NO signaling, along with genes encoding antioxidant enzymes and those involved in SMX metabolism. These findings indicated that C6-HSL, as a bacterial QS signal, could enhance microalgal tolerance and antibiotic degradation, offering a novel strategy to improve microalgae-based antibiotic removal in wastewater treatment.

Health Risks of Low-Dose Dietary Exposure to Triphenyl Phosphate and Diphenyl Phosphate in Mice: Insights from the Gut-Liver Axis

Aryl phosphate esters have been detected throughout the natural environment and in human blood samples, making it important to determine the health risks associated with exposure to triphenyl phosphate (TPHP) and its metabolite diphenyl phosphate (DPHP). Here, C57BL/6J male mice were exposed to TPHP or DPHP for 12 weeks at estimated daily intake doses of 0.1 and 7 μg/kg bw/day. TPHP intake affected the levels of short-chain fatty acids and bile acids in the gut, enhancing the production of 29 medium- and long-chain fatty acids in the liver by 3.72-fold and significantly increasing hepatic lipid and cholesterol levels. Metabolomic and molecular analysis confirmed that elevated liver cholesterol levels persisted after an 8 week recovery period. Gut microbiota-dependent cholesterol alterations were the toxic end points observed in TPHP-fed mice, as supported by the results of fecal microbiota transplantation. In DPHP-fed mice, serotonergic and glutamatergic synapses were simultaneously altered in the liver and intestine, corresponding to the reduction of five brain neurotransmitters (15.4-60.8%). Decreased liver carbohydrate levels and insulin resistance were observed in the DPHP-fed mice. These results suggest that TPHP and DPHP affect metabolism via different toxic modes, mediated through the gut-liver axis, providing novel insights into the mechanisms of organophosphate-ester-mediated metabolic disruption.

Potential emission risks of organic compounds from cement kilns

The extensive application of cement kiln industry results in substantial stack gas emissions, posing a potential risk of discharging organic pollutants. Cement industry is not considered as a primary contributor to persistent organic pollutants like polychlorinated dibenzo-p-dioxins and furans (PCDD/Fs), owing to its extremely low emission factor. However, knowledge on the previously unrecognized chemicals that may possess higher emission factors from cement industry is lacking. In this study, non-target screening of pollutants from cement industries were conducted using gas chromatography-quadrupole time-of-flight mass spectrometry. Priority pollutants from cement industry were identified. Alkylated polycyclic aromatic hydrocarbons and chlorobenzene were noteworthy pollutants from cement industry. By comparing the toxicity data and semi-quantitative of the pollutants, 9 chemical compounds, including 1-Methylphenanthrene, were identified as priority pollutants. The emission factors of these selected pollutants from cement industries were about 3-6 orders higher than PCDD/Fs, suggesting a potentially high risk of emission. The nationwide emission estimation indicates a significant release of approximately 300 kg of 1-Methylphenanthrene emanating from the cement industry, necessitating heightened attention. The results are helpful for comprehensive risk assessment of organic pollutants from cement industry.

Theoretical Investigation of Cytochrome P450 Enzyme-Mediated Biotransformation Mechanism of BHPF: Unveiling the Metabolic Safety Aspects of an Alternative to BPA

Fluorene-9-bisphenol (BHPF), emerging as an alternative to bisphenol A (BPA), is extensively utilized in industry and consumer goods. BHPF exhibits antiestrogenic effects and potential reproductive toxicity. Similar to BPA, BHPF can closely access the active site of the cytochrome P450 (CYP450) enzyme, interact with the Fe=O moiety, and potentially initiate metabolic reactions. Using density functional theory (DFT) calculations, we explored the mechanisms underlying BHPF activation using a CYP450 compound I (Cpd I) model, identifying several plausible products. Compared with the higher energy barriers associated with phenyl ring addition reactions, the formation of a phenoxyl-type radical through phenolic hydrogen atom abstraction, followed by OH rebound or radical coupling, represents an energetically favorable pathway. The OH rebound process yields three primary products: 9-(3,4-dihydroxyphenyl)-9-(4-hydroxyphenyl)fluorene (PRD1), semiquinone radical anion (PRD2), and 9-(4-hydroxyphenyl)fluorene carbocation (PRD3), constituting the major outcomes of the BHPF metabolic reaction. Importantly, a lipophilic ether metabolite, BHPF-O-BHPF (PRD4), formed through the coupling of phenoxyl radicals, reflects a widespread metabolic pathway observed in phenolic molecules. Despite constituting a minor proportion, the toxicity of this product necessitates increased attention. These findings contribute significantly to an enhanced understanding of the potential hazards associated with BHPF and other unknown chemical entities.

The investigation of the toxicity of organophosphorus flame retardants (OPFRs) by using in silico toxicity prediction platform ProTox- 3.0

From the past to the present, many chemicals have been used for the purpose of flame retardant. Due to PBDEs' (Polybrominated diphenyl ether) lipophilic and accumulative properties, some of them are banned from the market. As an alternative to these chemicals, OPFRs (organophosphorus flame retardants) have started to be used as flame retardants. In this article, acute toxicity profiles, mutagenicity, carcinogenicity, blood-brain barrier permeability, ecotoxicity and nutritional toxicity as also AHR, ER affinity and MMP, aromatase affinity, CYP2C9, CYP3A4 interaction of the of 16 different compounds of the OPFRs were investigated using a computational toxicology method; ProTox- 3.0. According to our results, eight compounds were found to be active in terms of carcinogenic effect, whereas two compounds were found to be active for mutagenicity. On the other hand, all compounds were found to be active in terms of blood-barrier permeability. Fourteen compounds and four compounds are found to have ecotoxic and nutritional toxic potency, respectively. Eight compounds were determined as active to AhR, and four chemicals were found to be active in Estrogen Receptor alpha. Eight chemicals were found to be active in terms of mitochondrial membrane potency. Lastly, three chemicals were found to be active in aromatase enzymes. In terms of CYP interaction potencies, eight compounds were found to be active in both CYP2C9 and CYP3A4. This research provided novel insights into the potential toxic effects of OPFRs. However, further studies are needed to evaluate their toxicity. Moreover, these findings lay the groundwork for in vitro and in vivo toxicity research.

Stereoselective In Vitro Metabolism, Hepatotoxicity, and Cytotoxic Effects of Four Enantiomers of the Fungicide Propiconazole

Propiconazole (PRO) is a chiral triazole fungicide that has been widely used for several years. However, its metabolic characteristics and hepatotoxicity in the chiral level environment remain unclear. In this study, the stereoselective behavior of PRO was investigated by using liver microsome incubation, cell viability assay, inhalation exposure, and molecular docking. Our results demonstrated that the isomers trans (-)-2R,4R-PRO and cis (+)-2R,4S-PRO exhibited slower metabolic rates in rat liver microsomes. The cytochrome P450 family 1 subfamily A polypeptide 2 enzyme was found to play a key role in the metabolism of PRO, contributing to its stereoselective behavior. Histopathological and cell viability results showed that exposure to rac-PRO could induce severe hepatotoxicity in mice. This effect might be related to the accumulation of cis (+)-2R,4S-PRO in the liver, which has a slow metabolism and is highly toxic. Our findings indicate that avoiding the application of cis (+)-2R,4S-PRO in agriculture can significantly reduce adverse effects on nontarget organisms.

Evidence that cadmium aggravate the toxicity of triphenyl phosphate in aquatic sediments to Corbicula fluminea

The ubiquitous co-existence of triphenyl phosphate (TPhP) and heavy metals in sediments raises significant biotoxicity concerns. However, uncertainty still exists regarding their combined toxicity to benthic organisms. Therefore, this research was conducted to elucidate the influences of cadmium (Cd) on TPhP toxicity to Corbicula fluminea (C. fluminea) in sediments. As a result, Cd promoted the accumulation of TPhP in C. fluminea and enhanced TPhP toxicity, manifested by damaged cell membranes and pronounced histological alterations. Molecular docking revealed that TPhP-Cd complexes exhibit greater binding affinity to cytochrome P4501A1 (CYP1A1) compared to TPhP alone. With the activity of CYP1A1 increasing, the biotransformation of TPhP was promoted in low-TPhP+Cd treatments (T5C0/T5C5/T5C35). Additionally, metabolites related to antioxidant defence and repair processes were reinforced to alleviate the toxicity of TPhP and Cd. However, excessive oxidative stress impaired the CYP1A1 activity in high-TPhP+Cd treatments (T35C0/T35C5/T35C35). Furthermore, metabolic pathway analysis revealed significant perturbations in the citrate cycle, alanine, aspartate and glutamate metabolism, purine metabolism, and pyrimidine metabolism. These disruptions weakened the repair capacity and aggravated apoptosis in digestive glands, potentially contributing to the synergistic toxicity of TPhP and Cd. The results highlight the ecological risks posed by TPhP in combination with heavy metals to benthic organisms.

Evidence for the formation of 6PPD-quinone from antioxidant 6PPD by cytochrome P450

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine (6PPD) as a rubber antioxidant has attracted global concern, since its ozone-oxidation product 6PPD-quinone (6PPDQ) was found to be the primary toxicant responsible for urban runoff mortality syndrome in coho salmon. However, the biotransformation fate and associated toxicological mechanism of 6PPD have not received much study yet. In this work, the in vitro assays showed 6PPD can be transformed into 6PPDQ by cytochromes P450 (CYP450) in human liver microsomes (HLMs) with 0.98 % production rate, and the adducts of 6PPDQ with calf thymus DNA and the N-N coupling product between 6PPD and 6PPDQ were further identified after 6PPD incubation in HLMs. Further evidence for the 6PPDQ formation can be obtained from the in vivo assays that the 6PPDQ-DNA adducts and 6PPD-N-N-6PPDQ dimer were detected in mice by oral gavage with 6PPD, and the latter dimer species was detected as well in 6PPD exposure to zebrafish larvae. Especially, the bioaccumulation property and high reactivity of 6PPDQ result in the continuous formation of the significant DNA adducts and 6PPD-N-N-6PPDQ dimer even in case of low production rate of biotransformation of 6PPD to 6PPDQ, which may provide potentially effective biomarkers for such process. DFT computations revealed the formation mechanism of 6PPDQ is the (N)H-abstraction of 6PPD by CYP450, followed by amino radical rebound at the nearby ortho-carbon, yielding a quinol intermediate due to spin delocalization, that might readily undergo further oxidation by CYP450 into 6PPDQ.

Species-specific metabolism of triphenyl phosphate and its mono-hydroxylated product by human and rat CYP2E1 and the carp ortholog

Organophosphorus flame retardants (PFRs) are a class of flame retardants and environmental pollutants with various biological effects. Recentstudies have evidenced activation of some PFRs by human CYP enzymes (including CYP2E1) for genotoxic effects. However, the activity of CYPs in fish species toward PFR metabolism remains unclear. This study was aimed on comparing the metabolism of triphenyl phosphate (TPHP) and 4-OH-TPHP in human, rat, and common carp, and the involvement of human CYP2E1 and its orthologs in the metabolism, by using fomepizole (4-MP, CYP2E1 inhibitor) as a modulator, in silico molecular docking and dynamics analyses. The rate of TPHP metabolism was apparently faster with human and rat, microsomes than with fish microsomes, the major metabolites were phosphodiester and hydroxylated phosphate, with 30-80 % of TPHP forming unidentified metabolites in the system of each species. 4-OH-TPHP was readily metabolized by both human and rat microsomes, whereas it was hardly metabolized in carp assays. Meanwhile, with 4-MP the transformation of TPHP to 4-OH-TPHP was enhanced in the human/rat systems while suppressed in the carp system. Moreover, the formation of unidentified metabolites in human and rat systems was mostly inhibited by 4-MP. Through molecular dynamics analysis TPHP and its primary metabolites showed high affinity for human and rat CYP2E1, as well as the carp ortholog (CYP2G1-like enzyme), however, the 4-OH-TPHP bond to the latter was too far from the heme to permit a biochemical reaction. This study suggests that the metabolism/activation of TPHP might be favored in mammals rather than carp, a fish species.

Biotransformation, Bioaccumulation, and Bioelimination of Triphenyl Phosphate and Its Dominant Metabolite Diphenyl Phosphate In Vivo

Aryl phosphorus flame retardants (aryl-PFRs), such as triphenyl phosphate (TPHP) and diphenyl phosphate (DPHP), are widely used worldwide. Understanding the fates of aryl-PFRs in vivo is crucial to assessing their toxicity and the risks they pose. Seven TPHP metabolites, including Phase I hydrolysis and hydroxylation and Phase II glucuronidation products, were identified in C57BL/6J male mice following subacute dietary exposure to aryl-PFRs (70 μg/kg body weight (bw)/day) for 7 days. TPHP was almost completely metabolized by mice (∼97%), with DPHP the major metabolite formed (34%-58%). In addition, mice were exposed to aryl-PFRs (7 μg/kg bw/day) for 12 weeks. Both TPHP and DPHP occurred at higher concentrations in the digestive tract (intestine and stomach), liver and heart. The total concentration of DPHP in all organs was 3.55-fold greater than that of TPHP. Recovery analysis showed that the rate of TPHP elimination from mouse organs reached 38%, while only 3%-5% of DPHP was removed, suggesting that the rates of degradation and elimination of DPHP were slower than TPHP and its bioaccumulation potential was higher. These results highlight the critical role of DPHP in the biotransformation, bioaccumulation, and bioelimination of TPHP, providing valuable insights into the fate of aryl-PFRs in vivo.

Organophosphate esters (OPEs) pollution characteristics, bioaccumulation and human consumption implication in wild marine organisms from the Yellow River Estuary, China

As the substitutes of polybrominated diphenyl ethers, organophosphate esters (OPEs) with high concentrations have accumulated in the estuaries, bays, and harbors. However, limited information is available about the OPEs in the estuary organism categories, especially under the multiple industrial pressure. This study investigated the occurrence, bioaccumulation and human consumption implication in wild marine organisms from the Yellow River Estuary, where located many petroleum and chemical manufacturing industries. This study found that concentrations of Σ13OPEs ranged from 547 ng/L to 1164 ng/L in seawater (median: 802 ng/L), from 384 to 1366 ng/g dw in the sediment (median: 601 ng/g dw), and from 419 to 959 ng/g dw (median: 560 ng/g dw) in the marine organisms. The congener compositions in the organisms were dominated by alkyl-OPEs (80.7 %), followed by halogenated-OPEs (18.8 %) and aryl-OPEs (0.5 %). Based on the principal component analysis, petrochemical pollution, and industrial wastewater discharge were distinguished as the main plausible sources of OPEs to the YRE ecosystem. Most OPEs had potential or strong bioaccumulation capacity on the organisms, with a positive correlation between log BAF (Bioaccumulation Factor) and log Kow of OPEs. The highest estimated daily intake value of OPEs was tri-n-propyl phosphate, exceeding 300 ng/kg·bw/day via consuming fish. The highest hazard quotients from OPEs ranged from 0.001 to 0.1, indicating a low risk to human health by consuming marine organisms in the YRE. As the consumption of OPEs increases year by year, the risks of OPEs still cannot be ignored.

Organophosphate esters in milk across thirteen countries from 2020 to 2023: Concentrations, sources, temporal trends and ToxPi priority to humans

Recent increases in organophosphate ester (OPE) application have led to their widespread presence, yet little is known about their temporal trends in food. This study collected milk samples from 13 countries across three continents during 2020-2023, finding detectable OPEs in all samples (range: 2.25-19.7; median: 7.06 ng/g ww). Although no statistical temporal differences were found for the total OPEs during the 4-year sampling campaign, it was interesting to observe significant variations in the decreasing trend for Cl-OPEs and concentration variations for aryl-OPEs and alkyl-OPEs (p < 0.05), indicating changing OPE use patterns. Packaged milk exhibited significant higher OPE levels than those found in directly collected raw unpackaged milk, and milk with longer shelf-life showed higher OPE levels, revealing packaging material as a contamination source. No significant geographical differences were observed in milk across countries (p > 0.05), but Shandong Province, a major OPE production site in China, showed relatively higher OPE concentrations. The Monte Carlo simulation of estimated daily intakes indicated no exposure risk from OPEs through milk consumption. The molecular docking method was used to assess human hormone binding affinity with OPEs, amongst which aryl-OPEs had the highest binding energies. The Toxicological-Priority-Index method which integrated chemical property, detection frequency, risk quotients, hazardous quotients and endocrine-disrupting effects was employed to prioritize OPEs. Aryl-OPEs showed the highest scores, deserving attention in the future.

Specific human CYP enzymes-dependent mutagenicity of tris(2-butoxyethyl) phosphate (an organophosphorus flame retardant) in human and hamster cell lines

Tris(2-butoxyethyl) phosphate (TBOEP) is an organophosphorus flame retardant ubiquitously present in the environment and even the human body. TBOEP is toxic in multiple tissues, which forms dealkylated and hydroxylated metabolites under incubation with human hepatic microsomes; however, the impact of TBOEP metabolism on its toxicity, particularly mutagenicity (typically requiring metabolic activation), is left unidentified. In this study, the mutagenicity of TBOEP in human hepatoma cell lines (HepG2 and C3A) and the role of specific CYPs were studied. Through molecular docking, TBOEP bound to human CYP1A1, 1B1, 2B6 and 3A4 with energies and conformations favorable for catalyzing reactions, while the conformations of its binding with human CYP1A2 and 2E1 appeared unfavorable. In C3A cells (endogenous CYPs being substantial), TBOEP exposing for 72 h (2-cell cycle) at low micromolar levels induced micronucleus, which was abolished by 1-aminobenzotriazole (inhibitor of CYPs); in HepG2 cells (CYPs being insufficient) TBOEP did not induce micronucleus, whose effect was however potentiated by pretreating the cells with PCB126 (CYP1A1 inducer) or rifampicin (CYP3A4 inducer). TBOEP induced micronucleus in Chinese hamster V79-derived cell lines genetically engineered for stably expressing human CYP1A1 and 3A4, but not in cells expressing the other CYPs. In C3A cells, TBOEP selectively induced centromere protein B-free micronucleus (visualized by immunofluorescence) and PIG-A gene mutations, and elevated γ-H2AX rather than p-H3 (by Western blot) which indicated specific double-strand DNA breaks. Therefore, this study suggests that TBOEP may induce DNA/chromosome breaks and gene mutations in human cells, which requires metabolic activation by CYPs, primarily CYP1A1 and 3A4.

Metabolomics Simultaneously Derives Benchmark Dose Estimates and Discovers Metabolic Biotransformations in a Rat Bioassay

Benchmark dose (BMD) modeling estimates the dose of a chemical that causes a perturbation from baseline. Transcriptional BMDs have been shown to be relatively consistent with apical end point BMDs, opening the door to using molecular BMDs to derive human health-based guidance values for chemical exposure. Metabolomics measures the responses of small-molecule endogenous metabolites to chemical exposure, complementing transcriptomics by characterizing downstream molecular phenotypes that are more closely associated with apical end points. The aim of this study was to apply BMD modeling to in vivo metabolomics data, to compare metabolic BMDs to both transcriptional and apical end point BMDs. This builds upon our previous application of transcriptomics and BMD modeling to a 5-day rat study of triphenyl phosphate (TPhP), applying metabolomics to the same archived tissues. Specifically, liver from rats exposed to five doses of TPhP was investigated using liquid chromatography-mass spectrometry and 1H nuclear magnetic resonance spectroscopy-based metabolomics. Following the application of BMDExpress2 software, 2903 endogenous metabolic features yielded viable dose-response models, confirming a perturbation to the liver metabolome. Metabolic BMD estimates were similarly sensitive to transcriptional BMDs, and more sensitive than both clinical chemistry and apical end point BMDs. Pathway analysis of the multiomics data sets revealed a major effect of TPhP exposure on cholesterol (and downstream) pathways, consistent with clinical chemistry measurements. Additionally, the transcriptomics data indicated that TPhP activated xenobiotic metabolism pathways, which was confirmed by using the underexploited capability of metabolomics to detect xenobiotic-related compounds. Eleven biotransformation products of TPhP were discovered, and their levels were highly correlated with multiple xenobiotic metabolism genes. This work provides a case study showing how metabolomics and transcriptomics can estimate mechanistically anchored points-of-departure. Furthermore, the study demonstrates how metabolomics can also discover biotransformation products, which could be of value within a regulatory setting, for example, as an enhancement of OECD Test Guideline 417 (toxicokinetics).

Mechanistic and predictive studies on the oxidation of furans by cytochrome P450: A DFT study

Furan-containing compounds distribute widely in food, herbal medicines, industrial synthetic products, and environmental media. These compounds can undergo oxidative metabolism catalyzed by cytochrome P450 enzymes (CYP450) within organisms, which may produce reactive products, possibly reacting with biomolecules to induce toxic effects. In this work, we performed DFT calculations to investigate the CYP450-mediated metabolic mechanism of furan-ring oxidation using 2-methylfuran as a model substrate, meanwhile, we studied the regioselective competition of another hydroxylation reaction involving methyl group of 2-methylfuran. As a result, we found the toxicological-relevant cis-enedione product can be produced from O-addition directly via a concerted manner without formation of an epoxide intermediate as traditionally believed. Moreover, our calculations demonstrate the kinetic and thermodynamic feasibility of both furan-ring oxidation and methyl hydroxylation pathways, although the former pathway is a bit more favorable. We then constructed a linear model to predict the rate-limiting activation energies (ΔE*) of O-addition with 11 diverse furan substates based on their adiabatic ionization potentials (AIPs) and condensation Fukui functions (CFFs). The results show a good predictive ability (R2=0.94, Q2CV=0.87). Therefore, AIP and CFF with clear physichem meanings relevant to the mechanism, emerge as pivotal molecular descriptors to enable the fast prediction of furan-ring oxidation reactivities for quick insight into the toxicological risk of furans, using just ground-state calculations.

Microalgae enhanced co-metabolism of sulfamethoxazole using aquacultural feedstuff components: Co-metabolic pathways and enzymatic mechanisms

The application of antibiotics in freshwater aquaculture leads to increased contamination of aquatic environments. However, limited information is available on the co-metabolic biodegradation of antibiotics by microalgae in aquaculture. Feedstuffs provide multiple organic substrates for microalgae-mediated co-metabolism. Herein, we investigated the co-metabolism of sulfamethoxazole (SMX) by Chlorella pyrenoidosa when adding main components of feedstuff (glucose and lysine). Results showed that lysine had an approximately 1.5-fold stronger enhancement on microalgae-mediated co-metabolism of SMX than glucose, with the highest removal rate (68.77% ± 0.50%) observed in the 9-mM-Lys co-metabolic system. Furthermore, we incorporated reactive sites predicted by density functional theory calculations, 14 co-metabolites identified by mass spectrometry, and the roles of 18 significantly activated enzymes to reveal the catalytic reaction mechanisms underlying the microalgae-mediated co-metabolism of SMX. In lysine- and glucose-treated groups, five similar co-metabolic pathways were proposed, including bond breaking on the nucleophilic sulfur atom, ring cleavage and hydroxylation at multiple free radical reaction sites, together with acylation and glutamyl conjugation on electrophilic nitrogen atoms. Cytochrome P450, serine hydrolase, and peroxidase play crucial roles in catalyzing hydroxylation, bond breaking, and ring cleavage of SMX. These findings provide theoretical support for better utilization of microalgae-driven co-metabolism to reduce sulfonamide antibiotic residues in aquaculture.

Uptake, translocation, and metabolism of organophosphate esters (OPEs) in plants and health perspective for human: A review

Plant uptake, accumulation, and transformation of organophosphate esters (OPEs) play vital roles in their geochemical cycles and exposure risks. Here we reviewed the recent research advances in OPEs in plants. The mean OPE concentrations based on dry/wet/lipid weight varied in 4.80-3,620/0.287-26.8/12,000-315,000 ng g-1 in field plants, and generally showed positive correlations with those in plant habitats. OPEs with short-chain substituents and high hydrophilicity, particularly the commonly used chlorinated OPEs, showed dominance in most plant samples, whereas some tree barks, fruits, seeds, and roots demonstrated dominance of hydrophobic OPEs. Both hydrophilic and hydrophobic OPEs can enter plants via root and foliar uptake, and the former pathway is mainly passively mediated by various membrane proteins. After entry, different OPEs undergo diverse subcellular distributions and acropetal/basipetal/intergenerational translocations, depending on their physicochemical properties. Hydrophilic OPEs mainly exist in cell sap and show strong transferability, hydrophobic OPEs demonstrate dominant distributions in cell wall and limited migrations owing to the interception of Casparian strips and cell wall. Additionally, plant species, transpiration capacity, growth stages, commensal microorganisms, and habitats also affect OPE uptake and transfer in plants. OPE metabolites derived from various Phase I transformations and Phase II conjugations are increasingly identified in plants, and hydrolysis and hydroxylation are the most common metabolic processes. The metabolisms and products of OPEs are closely associated with their structures and degradation resistance and plant species. In contrast, plant-derived food consumption contributes considerably to the total dietary intakes of OPEs by human, particularly the cereals, and merits specifical attention. Based on the current research limitations, we proposed the research perspectives regarding OPEs in plants, with the emphases on their behavior and fate in field plants, interactions with plant-related microorganisms, multiple uptake pathways and mechanisms, and comprehensive screening analysis and risk evaluation.

Mechanistic insight into biotransformation of novel triazine-based flame retardant 1,3,5-tris(2,3-dibromopropyl)-1,3,5-triazinane-2,4,6-trione by human cytochrome P450s

The escalating focus on the environmental occurrence and toxicology of emerging pollutants underscores the imperative need for a profound exploration of their metabolic transformations mediated by human CYP450 enzymes. Such investigations have the potential to unravel the intricate metabolite profiles, substantially altering the toxicological outcomes. In this study, we integrated the computational simulations with in vitro metabolism experiments to investigate the metabolic activity and mechanism of an emerging pollutant, 1,3,5-tris(2,3-dibromopropyl)-1,3,5-triazinane-2,4,6-trione (TDBP-TAZTO), catalyzed by human CYP450s. The results highlight the important contributions of CYP2E1, 3A4 and 2C9 to the biotransformation of TDBP-TAZTO, leading to the identification of four distinct metabolites. The effective binding conformations governing biotransformation reactions of TDBP-TAZTO within active CYP450s are unveiled. Structural instability of primary hydroxyTDBP-TAZTO products suggests three potential outcomes: (1) generation of an alcohol metabolite through successive debromination and reduction reactions, (2) formation of a dihydroxylated metabolite through secondary hydroxylation by CYP450, and (3) production of an N-dealkylated metabolite via decomposition and isomerization reactions in the aqueous environment. The formation of a desaturated debrominated metabolite may arise from H-abstraction and barrier-free Br release during the primary oxidation, potentially competing with the generation of hydroxyTDBP-TAZTO. These findings provide detailed mechanistic insight into TDBP-TAZTO biotransformation by CYP450s, which can enrich our understanding of the metabolic fate and associated health risk of this chemical.

In vitro metabolism of the emerging contaminant 6PPD-quinone in human and rat liver microsomes: Kinetics, pathways, and mechanism

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone (6PPD-Q) is an ozonation product of the rubber antioxidant N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine (6PPD). 6PPD-Q has recently been detected in various environmental media, which may enter the human body via inhalation and skin contact pathways. However, the human metabolism of 6PPD-Q has remained unknown. This study investigated the in vitro Cytochrome P450-mediated metabolism of 6PPD-Q in human and rat liver microsomes (HLMs and RLMs). 6PPD-Q was significantly metabolized at lower concentrations but slowed at high concentrations. The intrinsic clearance (CLint) of 6PPD-Q was 21.10 and 18.58 μL min-1 mg-1 protein of HLMs and RLMs, respectively, suggesting low metabolic ability compared with other reported pollutants. Seven metabolites and one intermediate were identified, and metabolites were predicted immunotoxic or mutagenic toxicity. Mono- and di-oxygenation reactions were the main phase I in vitro metabolic pathways. Enzyme inhibition experiments and molecular docking techniques were further used to reveal the metabolic mechanism. CYP1A2, 3A4, and 2C19, especially CYP1A2, play critical roles in 6PPD-Q metabolism in HLMs, whereas 6PPD-Q is extensively metabolized in RLMs. Our study is the first to demonstrate the in vitro metabolic profile of 6PPD-Q in HLMs and RLMs. The results will significantly contribute to future human health management targeting the emerging pollutant 6PPD-Q.

Unravelling bioaccumulation, depletion and metabolism of organophosphate triesters in laying hens: Insight of in vivo biotransformation assisted by diester metabolites

Organophosphate triesters (tri-OPEs) threaten human health through dietary exposure, but little is known about their feed-to-food transfer and in vivo behavior in farm animals. Herein 135 laying hens were fed with contaminated feed (control group, low-level group and high-level group) to elucidate the bioaccumulation, distribution, and metabolism of the six most commonly reported tri-OPEs. The storage (breast muscle), metabolism and mobilization (liver and blood) and non-invasive (feather) tissues were collected. The exposure-increase (D1∼14) and depuration-decrease (D15∼42) trends indicated that feed exposure caused tri-OPE accumulation in animal tissues. Tissue-specific and moiety-specific behavior was observed for tri-OPEs. The highest transfer factor (TF) and transfer rate (TR) were observed in liver (TF: 14.8%∼82.3%; TR: 4.40%∼24.5%), followed by feather, breast muscle, and blood. Tris(2-chloroisopropyl) phosphate (TCIPP) had the longest half-life in feather (72.2 days), while triphenyl phosphate (TPhP) showed the shortest half-life in liver (0.41 days). Tri-OPEs' major metabolites (organophosphate diesters, di-OPEs) were simultaneously studied, which exhibited dose-dependent and time-dependent variations following administration. In breast muscle, the inclusion of di-OPEs resulted in TF increases of 735%, 1108%, 798%, and 286% than considering TCIPP, tributyl phosphate, tris(2-butoxyethyl) phosphate and tris(2-ethylhexyl) phosphate alone. Feather was more of a proxy of birds' long-term exposure to tri-OPEs, while short-term exposure was better reflected by di-OPEs. Both experimental and in silico modeling methods validated aryl-functional group facilitated the initial accumulation and metabolism of TPhP in the avian liver compared to other moiety-substituted tri-OPEs.

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.

Lipid metabolic disturbance induced by triphenyl phosphate and hydroxy metabolite in HepG2 cells

Triphenyl phosphate (TPHP) has been widely used as flame retardants and been detected with increasing frequency in environment. TPHP can transform into mono-hydroxylated phosphate (OH-TPHP) and diester diphenyl phosphate (DPHP) through biotransformation. So far, information on the cytotoxicity and molecular regulatory mechanisms of TPHP metabolites are still limit. This study investigated the adverse effects of TPHP, OH-TPHP, and DPHP in HepG2 cells in terms of cell proliferation, lactate dehydrogenase release, reactive oxygen species generation, and mitochondrial membrane potential. The transcriptomic changes were measured using RNA sequencing, and bioinformatics characteristics including biological functions, signal pathways and protein-protein interaction were analyzed to explore the potential molecular mechanisms. Results displayed that the order of cytotoxicity was OH-TPHP> TPHP> DPHP. The prioritized biological functions changes induced by TPHP and OH-TPHP were correlated with lipid metabolism. Significant lipid accumulation was observed as confirmed by increased total cholesterol and triglycerides contents, and enhanced oil red O staining. Enrichment of PPARα/γ and down-stream genes suggested the participation of PPARs signal pathway in lipid metabolism disorder. In addition, TPHP and OH-TPHP induced endoplasmic reticulum stress (ERS), which was further confirmed by the ERS inhibitor experiment. In general, TPHP and OH-TPHP had obvious cytotoxic effects in HepG2 cells. PPARs signal pathway and endoplasmic reticulum stress may be involved in the lipid metabolism disorder induced by TPHP and OH-TPHP.

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.

New insight into molecular mechanism of P450-Catalyzed metabolism of emerging contaminants and its consequence for human health: A case study of preservative methylparaben

Hydroxylated metabolites in the living body are considered as a potential biomarker of exposure to emerging contaminations (ECs) and breast cancer, but their formation mechanism has not received enough attention. Besides, the adverse impacts of metabolites during the metabolic transformation of ECs largely remain unknown. In this study, we employed a density functional calculation combing with in-vitro incubation of human liver microsomes to explore the bio-transformation of preservative methylparaben (MPB) in human bodies. Our results showed that hydroxylated metabolites of MPB (OH-MPB) were observed experimentally, while a formation mechanism was revealed at the molecular level. That is, hydroxylated metabolite was exclusively formed via the hydrogen abstraction from the phenolic hydroxyl group of MPB followed by the OH-rebound pathway, rather than the direct hydroxylation on the benzene ring. The increasing of hydroxyl groups on ECs could improve the metabolisms. This was confirmed in the metabolism of ECs without hydroxyl group and with multiple-hydroxyl groups, respectively. Furthermore, toxicity assessments show that compared to parent MPB, the hydroxylated metabolites have increased negative impacts on the gastrointestinal system and liver. A semiquinone product exhibits potential damage in the cardiovascular system and epoxides are toxic to the blood and gastrointestinal system. The findings deepen our insight into the biotransformation of parabens in human health, especially by providing health warnings about the potential impacts caused by semiquinone and epoxides.

An Insight into the Metabolism of 2,5-Disubstituted Monotetrazole Bearing Bisphenol Structures: Emerging Bisphenol A Structural Congeners

The non-estrogenic 2,5-disubstituted tetrazole core-bearing bisphenol structures (TbB) are being researched as emerging structural congeners of Bisphenol A, an established industrial endocrine disruptor. However, there is no understanding of TbB's adverse effects elicited via metabolic activation. Therefore, the current study aimed to investigate the metabolism of TbB ligands, with in silico results serving as a guide for in vitro studies. The Cytochrome P450 enzymes (CYP) inhibitory assay of TbB ligands on the seven human liver CYP isoforms (i.e., 1A2, 2A6, 2D6, 2C9, 2C8, 2C19, and 3A4) using human liver microsomes (HLM) revealed TbB ligand 223-3 to have a 50% inhibitory effect on all the CYP isoforms at a 10 μM concentration, except 1A2. The TbB ligand 223-10 inhibited 2B6 and 2C8, whereas the TbB ligand 223-2 inhibited only 2C9. The first-order inactivity rate constant (Kobs) studies indicated TbB ligands 223-3, 223-10 to be time-dependent (TD) inhibitors, whereas the TbB 223-2 ligand did not show such a significant effect. The 223-3 exhibited a TD inhibition for 2C9, 2C19, and 1A2 with Kobs values of 0.0748, 0.0306, and 0.0333 min-1, respectively. On the other hand, the TbB ligand 223-10 inhibited 2C9 in a TD inhibition manner with Kobs value 0.0748 min-1. However, the TbB ligand 223-2 showed no significant TD inhibition effect on the CYPs. The 223-2 ligand biotransformation pathway by in vitro studies in cryopreserved human hepatocytes suggested the clearance via glucuronidation with the predominant detection of only 223-2 derived mono glucuronide as a potential inactive metabolite. The present study demonstrated that the 223-2 ligand did not elicit any significant adverse effect via metabolic activation, thus paving the way for its in vivo drug-drug interactions (DDI) studies.

Triphenyl phosphate proved more potent than its metabolite diphenyl phosphate in inducing hepatic insulin resistance through endoplasmic reticulum stress

Triphenyl phosphate (TPHP) is a widely used flame retardant and plasticizer and has been detected extensively in environmental media, wildlife and human bodies. Several epidemiological and animal studies have revealed that TPHP exposure is positively associated with glucose homeostasis disruption and diabetes. However, the effects of TPHP on hepatic glucose homeostasis and the underlying mechanisms remain unclear. The present work aimed to investigate the cytotoxicity and glucose metabolism disruption of TPHP and its metabolite diphenyl phosphate (DPHP) within hepatocytes. The cell viability assay undertaken on human normal liver (L02) cells showed that TPHP exhibited more potent hepatotoxicity than DPHP. RNA sequencing (RNA-seq) data showed that TPHP and DPHP presented different modes of toxic action. Insulin resistance is one of the predominant toxicities for TPHP, but not for DPHP. The insulin-stimulated glucose uptake and glycogen synthesis were impaired by TPHP, while DPHP exhibited no significant impairment on these factors. TPHP exposure induced endoplasmic reticulum (ER) stress, and the ER stress antagonist 4-PBA restored the impairment of insulin-stimulated glucose uptake and glycogen synthesis induced by TPHP. TPHP could also induce liver ER stress and insulin resistance in mice. Taken together, the results suggested that TPHP induces more potent insulin resistance through ER stress than its metabolite DPHP.

First Insight into the Formation of In Vivo Transformation Products of 2-Ethylhexyl diphenyl phosphate in Zebrafish and Prediction of Their Potential Toxicities

As a frequently detected organophosphorus flame retardant in the environment, 2-ethylhexyl diphenyl phosphate (EHDPHP) is vulnerable to biotransformation, while the transformation mechanisms and potential toxicities of its transformation products remain unclear. In the present study, in vivo transformation products of EHDPHP in exposed zebrafish for 21d were analyzed by suspect screening and identified by mass spectrometry. Fifteen metabolites were identified, including 10 phase I and 5 phase II products with monohydroxylated products being primary, among which 5-OH-EHDPHP was the most predominant. Two sulfation products and one terminal desaturation metabolite of EHDPHP were reported for the first time. A density functional calculation coupled with molecular docking disclosed that the specific conformation of EHDPHP docked in the protein pockets favored the primary formation of 5-OH-EHDPHP, which was fortified to be a more suitable biomarker of EHDPHP exposure. The in vitro tests suggested that EHDPHP transformation took place not only in liver but also in intestine, where gut microbes played an important role. Due to lack of standards, in silico toxicity prediction combined with molecular docking indicated that several metabolites potentially cause higher toxicities than EHDPHP. The results provide deep insight into the potential health risks due to specific in vivo transformation of EHDPHP.

Triphenyl phosphate induces clastogenic effects potently in mammalian cells, human CYP1A2 and 2E1 being major activating enzymes

As a new-type flame retardant and toxic substance, triphenyl phosphate (TPP) is a ubiquitous pollutant present even in human blood. TPP is transformed by human CYP enzymes to oxidized/dealkylated metabolites. The impact of TPP metabolism on its toxicity, however, remains unclear. In this study, the genotoxicity of TPP in several mammalian cell lines and its relevance to CYP/sulfortransferase (SULT) activities were investigated. The results indicated that TPP induced micronucleus formation at ≥1 μM concentrations in a human hepatoma (C3A, endogenous CYPs being substantial) cell line, which was abolished by 1-aminobenzotriazole (CYPs inhibitor). In cell line HepG2 (parental to C3A with lower CYP expression) TPP was inactive up to 10 μM, while pretreatment with ethanol (CYP2E1 inducer), PCB 126 (CYP1A inducer), or rifampicin (CYP3A inducer) led to micronucleus formation by TPP. In V79-Mz and V79-derived cells expressing human CYP1A1 TPP was inactive (up to 32 μM), and in cells expressing human CYP1B1, 2B6 and 3A4 it induced micronucleus weakly (positive only at 32 μM). However, TPP induced micronucleus potently in V79-derived cells expressing human CYP1A2, while this effect was drastically reduced by human SULT1A1 co-expression; likewise, TPP was inactive in cells expressing both human CYP2E1 and SULT1A1, but became positive with pentachlorophenol (inhibitor of SULT1) co-exposure. Moreover, in C3A cells TPP selectively induced centromere-free micronucleus (immunofluorescent assay), and TPP increased γ-H2AX (by Western blot, indicating double-strand DNA breaks). In conclusion, this study suggests that TPP is potently clastogenic, human CYP1A2 and 2E1 being major activating enzymes while SULT1A1 involved in detoxification.

Accumulation and translocation of traditional and novel organophosphate esters and phthalic acid esters in plants during the whole life cycle

Organophosphate esters (OPEs) and phthalic acid esters (PAEs) are widespread contaminants in the environment. The variations of these chemicals in plants throughout their life cycle is little known. In this study, OPEs, OPE metabolites, and PAEs in peanut and corn grown under field conditions, soil, and air were measured to understand the uptake and translocation, distributions in the plant compartments, and metabolism in the plants. The soil concentrations showed an enrichment effect of OPEs onto the rhizosphere soil but a depletion effect of PAEs on rhizosphere soils. The PAE concentrations between peanut (with a mean of 1295 ng/g dw) and corn (3339 ng/g dw) were significantly different, but the OPE concentrations were not significantly different (with means of 15.6 and 19.2 ng/g dw, respectively). OPE metabolites were also detected in the plants, with lower concentrations and detection rates. Similarities and differences in the temporal variations of the concentrations of traditional OPEs, novel OPEs, and PAEs in plants during their growth were observed. The variations were dependent on both plant species and particular tissues. The leaf compartment is the most important reservoir of OPEs and PAEs (but not OPE metabolites) for both species, highlighting the importance of an aerial uptake pathway. The chemicals have a low potential to be translocated into peanut and corn kernels, reducing their risks via food consumption. Less hydrophobic compounds have higher root concentration factors in this study. These observations differ from those of previous hydroponic experiments.

Biotransformation of sulfamethoxazole by microalgae: Removal efficiency, pathways, and mechanisms

Recently, the biotransformation of sulfamethoxazole (SMX) by microalgae has attracted increasing interest. In particular, cytochrome P450 (CYP450) has been suggested to be the main enzymatic contributor to this biodegradation. However, the molecular evidence of CYP450 enzymes being involved in SMX biodegradation remains relatively unclear, hindering its applicability. Herein, the biodegradation of SMX by Chlorella sorokiniana (C. sorokiniana) was investigated, and comprehensively elucidated the reaction mechanism underlying CYP450-mediated SMX metabolism. C. sorokiniana was able to efficiently remove over 80% of SMX mainly through biodegradation, in which CYP450 enzymes responded substantially to metabolize SMX in cells. Additionally, screening of transformation products (TPs) revealed that N⁴-hydroxylation-SMX (TP270) was the main TP in the SMX biodegradation pathway of microalgae. Molecular dynamics (MD) simulation suggested that the aniline of SMX was the most prone to undergo metabolism, while density functional theory (DFT) indicated that SMX was metabolized by CYP450 enzymes through H-abstraction-OH-rebound reaction. Collectively, this work reveals key details of the hydroxylamine group of SMX, elucidates the SMX biodegradation pathway involving CYP450 in microalgae in detail, and accelerates the development of using microalgae-mediated CYP450 to eliminate antibiotics.

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.

Using Physical Organic Chemistry Knowledge to Predict Unusual Metabolites of Synthetic Phenolic Antioxidants by Cytochrome P450

Biotransformation, especially by human CYP450 enzymes, plays a crucial role in regulating the toxicity of organic compounds in organisms, but is poorly understood for most emerging pollutants, as their numerous "unusual" biotransformation reactions cannot retrieve examples from the textbooks. Therefore, in order to predict the unknown metabolites with altering toxicological profiles, there is a realistic need to develop efficient methods to reveal the "unusual" metabolic mechanism of emerging pollutants. Combining experimental work with computational predictions has been widely accepted as an effective approach in studying complex metabolic reactions; however, the full quantum chemical computations may not be easily accessible for most environmentalists. Alternatively, this work practiced using the concepts from physical organic chemistry for studying the interrelationships between structure and reactivity of organic molecules, to reveal the "unusual" metabolic mechanism of synthetic phenolic antioxidants catalyzed by CYP450, for which the simple pencil-and-paper and property-computation methods based on physical organic chemistry were performed. The phenol-coupling product of butylated hydroxyanisole (BHA) (based on spin aromatic delocalization) and ipso-addition quinol metabolite of butylated hydroxytoluene (BHT) (based on hyperconjugative effect) were predicted as two "unusual" metabolites, which were further confirmed by our in vitro analysis. We hope this easily handled approach will promote environmentalists to attach importance to physical organic chemistry, with an eye to being able to use the knowledge gained to efficiently predict the fates of substantial unknown synthesized organic compounds in the future.

Twenty-four hours of Thiamethoxam: In vivo and molecular dynamics simulation study on the toxicokinetic and underlying mechanisms in quails (Coturnix japonica)

Neonicotinoids is the most widely used insecticide, its contamination has led to sustained bird population declines. However, the toxicokinetic and underlying mechanisms of neonicotinoid toxicity in birds are largely unknown. Thiamethoxam (TMX), as a representative neonicotinoid insecticide, is now widely detected in most environmental medium and animal bodies. In this study, 5 mg/kg body weight TMX (potential environmental intake level) were orally administrated to male Japanese quails (Coturnix japonica). We found a rapid absorption, distribution, metabolism and elimination of TMX in quails in a period of 24 h, with the main metabolite, clothianidin (CLO), being extensively distributed and rapidly eliminated from tissues as well. The maximum plasm concentration of CLO was consistent with wild birds. Metabolomics analysis and followed determination of liver enzymes mRNA expression indicated the rapid metabolism was mediated mainly by CYPs and GSTs that involved riboflavin metabolism and glutathione metabolism pathways upon TMX exposure. Molecular dynamic simulation showed the strongest binding interaction in quail CYP2H1-TMX and CYP3A12-CLO complexes among a set of CYPs-substrate. The present study elucidated toxicokinetic and underlying metabolic mechanisms of TMX in quails at environmentally-relevant concentration, the findings would facilitate the understanding of potential risks of TMX and its metabolites to birds.

Enantioselective Metabolic Mechanism and Metabolism Pathway of Pydiflumetofen in Rat Liver Microsomes: In Vitro and In Silico Study

Pydiflumetofen (PYD) has been used worldwide. However, the enantioselective fate of PYD within mammals is not clear. Thus, the enantioselective metabolism and its potential mechanisms of PYD were explored via in vitro and in silico. Consistent results were observed between metabolism and enzyme kinetics experiments, with S-PYD metabolizing faster than R-PYD in rat liver microsomes. Moreover, CYP3A1 and carboxylesterase 1 were found to be major enzymes participating in the metabolism of PYD. Based on the computational results, S-PYD bound with CYP3A1 and carboxylesterase 1 more tightly with lower binding free energy than R-PYD, explaining the mechanism of enantioselective metabolism. Nine phase I metabolites of PYD were identified, and metabolic pathways of PYD were speculated. This study is the first to clarify the metabolism of PYD in mammals, and further research to evaluate the toxicological implications of these metabolites will help in assessing the risk of PYD.

Computational Investigation of the Bisphenolic Drug Metabolism by Cytochrome P450: What Factors Favor Intramolecular Phenol Coupling

Intramolecular phenol coupling reactions of alkaloids can lead to active metabolites catalyzed by the mammalian cytochrome P450 enzyme (P450); however, the mechanistic knowledge of such an "unusual" process is lacking. This work performs density functional theory computations to reveal the P450-mediated metabolic pathway leading from R-reticuline to the morphine precursor salutaridine by exploring possible intramolecular phenol coupling mechanisms involving diradical coupling, radical addition, and electron transfer. The computed results show that the outer-sphere electron transfer with a high barrier (>20.0 kcal/mol) is unlikely to happen. However, for inter-sphere intramolecular phenol coupling, it reveals that intramolecular phenol coupling of R-reticuline proceeds via the diradical mechanism consecutively by compound I and protonated compound II of P450 rather than the radical addition mechanism. The existence of a much higher radical rebound barrier than that of H-abstraction in the quartet high-spin state can endow the R-reticuline phenoxy radical with a sufficient lifetime to enable intramolecular phenol coupling, while the H-abstraction/radical rebound mode with a negligible rebound barrier leading to phenol hydroxylation can only happen in the doublet low-spin state. Therefore, the ratio [coupling]/[hydroxylation] can be approximately reflected by the relative yield of the high-spin and low-spin H-abstraction by P450, which thus can provide a theoretical ratio of 16:1 for R-reticuline, which is in accordance with previous experimental results. Especially, the high rebound barrier of the phenoxy radical derived from the weak electron-donating ability of the phenoxy radical is revealed as an intrinsic nature. Therefore, the revealed intramolecular phenol coupling mechanism can be potentially extended to several other bisphenolic drugs to infer groups of unexpected metabolites in organisms.

Computational Study of the Peroxygenase Mechanism Catalyzed by Hemoglobin Dehaloperoxidase Involved in the Degradation of Chlorophenols

The biochemical evidence showed that hemoglobin dehaloperoxidase (DHP B) from Amphitrite Ornata is a multifunctional hemoprotein that catalyzes both dehalogenation and hydroxylation of halophenols via the peroxidase and peroxygenase mechanism, respectively, which sets the basis for the degradation of halophenols. In the peroxygenase mechanism, the reaction was previously suggested to be triggered either by the hydrogen atom abstraction by the Fe═O center or by the proton abstraction by His55. To illuminate the peroxygenase mechanism of DHP B at the atomistic level, on the basis of the high-resolution crystal structure, computational models were constructed, and a series of quantum mechanical/molecular mechanical calculations have been performed. According to the calculation results, the pathway (Path a) initiated by the H-abstraction by the Fe═O center is feasible. In another pathway (Path b), His55 can abstract the proton from the hydroxyl group of the substrate (4-Cl-o-cresol) to initiate the reaction; however, its feasibility depends on the prior electron transfer from the substrate to the porphyrin group. The rate-limiting step of Path a is the OH-rebound, which corresponds to an energy barrier of 14.7 kcal/mol at the quartet state. His55 acts as an acid-base catalyst and directly involves in the catalysis. Our mutant study indicates that His55 can be replaced by other titratable residues. These findings may provide useful information for further understanding of the catalytic reaction of DHP B and for the design of enzymes in the degradation of pollutants, in particular, halophenols.

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.

Regioselective hydroxylation of carbendazim by mammalian cytochrome P450: A combined experimental and computational study

Carbendazim (CBZ), a broad-spectrum pesticide frequently detected in fruits and vegetables, could trigger potential toxic risks to mammals. To facilitate the assessment of health risks, this study aimed to characterize the cytochrome P450 (CYPs)-mediated metabolism profiles of CBZ by a combined experimental and computational study. Our results demonstrated that CYPs-mediated region-selective hydroxylation was a major metabolism pathway for CBZ in liver microsomes from various species including rat, mouse, minipig, dog, rabbit, guinea pig, monkey, cow and human, and the metabolite was biosynthesized and well-characterized as 6-OH-CBZ. CYP1A displayed a predominant role in the region-selective hydroxylation of CBZ that could attenuate its toxicity through converting it into a less toxic metabolite. Meanwhile, five other common pesticides including chlorpyrifos-methyl, prochloraz, chlorfenapyr, chlorpyrifos, and chlorothalonil could significantly inhibit the region-selective hydroxylation of CBZ, and consequently remarkably increased CBZ exposure in vivo. Furthermore, computational study clarified the important contribution of the key amino acid residues Ser122, and Asp313 in CYP1A1, as well as Asp320 in CYP1A2 to the hydroxylation of CBZ through hydrogen bonds. These results would provide some useful information for the metabolic profiles of CBZ by mammalian CYPs, and shed new insights into CYP1A-mediated metabolic detoxification of CBZ and its health risk assessment.

Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate

Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcohol_β-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and βC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the βC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcohol_β-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehyde_β-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.

Human internal exposure to organophosphate esters: A short review of urinary monitoring on the basis of biological metabolism research

As alternatives to traditional brominated flame retardants, organophosphate flame retardants (OPFRs), especially for organophosphate esters (OPEs) -- the most widely used and investigated OPFRs, have raised people's concern on their environmental and health-related risks over the years. Considering their extensive environmental occurrence and potential adverse effects, precise estimation on the human body burden of OPEs will be conducive to the restrictions on the usage of these compounds scientifically. Biomonitoring research can provide precise information on human exposure to OPEs as it reveals the degree of external exposure from all exposure routes. Knowledge on biotransformation and metabolism of OPEs in the biosystems is of great significance for our understanding of the internal exposure to these compounds. In this study, the biological metabolic processes of nine OPEs prevalent in the environment, involving tris(2-chloroethyl) phosphate (TCEP), tris(1-chloro-2-propyl) phosphate (TCIPP), tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tripropyl phosphate (TPrP), tri-n-butyl phosphate (TnBP), tris(2-butoxyethyl) phosphate (TBOEP), triphenyl phosphate (TPhP), 2-ethylhexyl diphenyl phosphate (EHDPP), and tricresyl phosphate (TCrP), are comprehensively reviewed. Specifically, the metabolic pathway, kinetics and mechanism of OPEs are depicted in detail. Under this context, the advances and limitations on biomonitoring of OPE metabolites in human urine are summarized. The requirements of specificity, quantitative stability, high detection frequency/concentration are needed for OPE metabolites to be considered and validated as biomarkers. Thus far, deeper elucidations on the metabolic processes and identification of biomarkers of OPEs are urgently required, given that some OPEs have no suitable biomarkers in human biomonitoring. For better assessment of the body burden of OPEs in humans, reliable and effective methodologies for urine sampling and estimation on internal exposure to OPEs need to be further developed in the future.

Stereostructure-activity mechanism of cyproconazole by cytochrome P450 in rat liver microsomes: A combined experimental and computational study

Cyproconazole (CPZ), representing the chiral triazole fungicides, is widely used in the pharmaceutical and agricultural fields. To clarify its potential adverse effects on the generalized CYP-mediated processes within mammalian, a comparative experimental and computational approach was employed to investigate the CYP-mediated metabolism processes of CPZ stereoisomers in rat liver microsomes (RLMs). The depletion rate of CPZ stereoisomers in vitro incubation system with RLMs followed the order RR-> SS-> SR-> RS-CPZ. The results of kinetic assays were in line with the depletion rate results. Further inhibition assay confirmed the stereoselective metabolism of CPZ stereoisomers by different CYP isoforms. Molecular dynamics (MD) simulation revealed the stereoselective metabolism mechanism. Several hydrogen bonds and π-stacking restrict the position of CPZ isomers in the active cavity of CYPs so that the 4'-nitrogen on the triazole ring can bind closely to the heme of CYP, which results in the metabolism of CPZ isomers. By combining the computational and experimental approaches, the structure-activity relationship of CPZ and CYP was elucidated, and this method can be further applied to predict the degree of uncertainty in the process of xenobiotic biotransformation of triazole fungicides and serve as a basis for risk assessment.

Long-Range Transport, Trophic Transfer, and Ecological Risks of Organophosphate Esters in Remote Areas

Organophosphate esters (OPEs) have been a focus in the field of environmental science due to their large volume production, wide range of applications, ubiquitous occurrence, potential bioaccumulation, and worrisome ecological and health risks. Varied physicochemical properties among OPE analogues represent an outstanding scientific challenge in studying the environmental fate of OPEs in recent years. There is an increasing number of studies focusing on the long-range transport, trophic transfer, and ecological risks of OPEs. Therefore, it is necessary to conclude the OPE pollution status on a global scale, especially in the remote areas with vulnerable and fragile ecosystems. The present review links together the source, fate, and environmental behavior of OPEs in remote areas, integrates the occurrence and profile data, summarizes their bioaccumulation, trophic transfer, and ecological risks, and finally points out the predominant pollution burden of OPEs among organic pollutants in remote areas. Given the relatively high contamination level and bioaccumulation/biomagnification behavior of OPEs, in combination with the sensitivity of endemic species in remote areas, more attention should be paid to the potential ecological risks of OPEs.

Molecular neural crest cell markers enable discrimination of organophosphates in the murine cardiac embryonic stem cell test

•

Organophosphates induced distinctive effects on neural crest cells within the ESTc.

•

Neural crest gene transcripts were of added value to the original ESTc read-out.

•

Mechanistic information adds value to the applicability of the ESTc.

The cardiac embryonic stem cell test (ESTc) originally used the differentiation of beating cardiomyocytes for embryotoxicity screenings of compounds. However, the ESTc consists of a heterogeneous cell population, including neural crest (NC) cells, which are important contributors to heart development in vivo. Molecular markers for NC cells were investigated to explore if this approach improved discrimination between structurally related chemicals, using the three organophosphates (OP): chlorpyrifos (CPF), malathion (MLT), and triphenyl phosphate (TPP). To decrease the test duration and to improve the objective quantification of the assay read-out, gene transcript biomarkers were measured on study day 4 instead of the traditional cardiomyocyte beating assessment at day 10. Gene expression profiling and immunocytochemistry were performed using markers for pluripotency, proliferation and cardiomyocyte and NC differentiation. Cell proliferation was also assessed by measurements of embryoid body (EB) size and total protein quantification (day 7). Exposure to the OPs resulted in similar patterns of inhibition of beating cardiomyocyte differentiation and of myosin protein expression on day 10. However, these three chemically related compounds induced distinctive effects on NC cell differentiation, indicated by changes in expression levels of the NC precursor (Msx2), NC marker (Ap2α), and epithelial to mesenchymal transition (EMT; Snai2) gene transcripts. This study shows that investigating NC markers can provide added value for ESTc outcome profiling and may enhance the applicability of this assay for the screening of structurally related test chemicals.

Phosphorus Deficiency Promoted Hydrolysis of Organophosphate Esters in Plants: Mechanisms and Transformation Pathways

The biotransformation of organophosphate esters (OPEs) in white lupin (Lupinus albus) and wheat (Triticum aestivum L.) was investigated in hydroponic experiments with different phosphorus (P)-containing conditions. The hydrolysis rates of OPEs followed the order of triphenyl phosphate (TPHP) > tri-n-butyl phosphate (TnBP) > tris(1,3-dichloro-2-propyl) phosphate (TDCPP). Hydrolysis of OPEs was accelerated at P-deficient conditions, and faster hydrolysis took place in white lupin than in wheat. Coincidingly, the production of acid phosphatase (ACP) in both plants was promoted, and much higher intracellular and extracellular ACPs were observed in white lupin under P-deficient conditions. In vitro experiments revealed that ACP was a key enzyme to hydrolyze OPEs. The hydrolysis rates of OPEs were significantly correlated with the Hirshfeld charges, calculated by density functional theory, of the oxygen atom in the single P-O bond. Using ultra-high-performance liquid chromatography coupled with Orbitrap Fusion mass spectrometer, 30 metabolites were successfully identified. Some of these metabolites, such as sulfate-conjugated products, hydration of cysteine-conjugated products of TPHP, and reductively dechlorinated metabolites of TDCPP, were observed for the first time in plants. It is noteworthy that OPEs may transform into many hydroxylated metabolites, and special attention should be paid to their potential environmental effects.

Flame Retardants-Mediated Interferon Signaling in the Pathogenesis of Nonalcoholic Fatty Liver Disease

Nonalcoholic fatty liver disease (NAFLD) is a growing concern worldwide, affecting 25% of the global population. NAFLD is a multifactorial disease with a broad spectrum of pathology includes steatosis, which gradually progresses to a more severe condition such as nonalcoholic steatohepatitis (NASH), fibrosis, cirrhosis, and eventually leads to hepatic cancer. Several risk factors, including exposure to environmental toxicants, are involved in the development and progression of NAFLD. Environmental factors may promote the development and progression of NAFLD by various biological alterations, including mitochondrial dysfunction, reactive oxygen species production, nuclear receptors dysregulation, and interference in inflammatory and immune-mediated signaling. Moreover, environmental contaminants can influence immune responses by impairing the immune system’s components and, ultimately, disease susceptibility. Flame retardants (FRs) are anthropogenic chemicals or mixtures that are being used to inhibit or delay the spread of fire. FRs have been employed in several household and outdoor products; therefore, human exposure is unavoidable. In this review, we summarized the potential mechanisms of FRs-associated immune and inflammatory signaling and their possible contribution to the development and progression of NAFLD, with an emphasis on FRs-mediated interferon signaling. Knowledge gaps are identified, and emerging pharmacotherapeutic molecules targeting the immune and inflammatory signaling for NAFLD are also discussed.

Increased adverse effects during metabolic transformation of short-chain chlorinated paraffins by cytochrome P450: A theoretical insight into 1-chlorodecane

Short-chain chlorinated paraffins (SCCPs), frequently detected in human tissues or organs, can result in threat to human health by disturbing normal metabolism. However, their metabolism mechanisms and fates are largely unclear. Therefore, to better understand the impacts of SCCPs and their metabolites on the human health, the metabolic mechanism and kinetics of SCCPs by cytochrome P450 enzymes (CYPs) were explored using density functional theory employed 1-chlorodecane as a model SCCPs. The results show that 1-chlorodecane could be readily metabolized by CYPs, and the rate constant reaches up 42.3 s^-1 in human body. Dechlorination of 1-chlorodecane is unlikely to occur and hydroxylation is dominated via H-abstraction pathways, especially from the intermediate C atom of 1-chlorodecane. The toxicity assessments suggest that the two metabolites, 10-chloro-decan-5-ol and 1-chlorodecanol could exhibit higher bioaccumulation, carcinogenicity and more serious damage on cardiovascular system after the metabolism of 1-chlorodecane. To our knowledge, this is the first study from the viewpoint of theoretical analysis to explore the metabolism of typical SCCPs in human body. It may provide deep insight into the metabolic transformation mechanism of SCCPs and cause the concerns about the adverse effects of their metabolites in human body.

Transplacental Behaviors of Organophosphate Tri- and Diesters Based on Paired Human Maternal and Cord Whole Blood: Efficiencies and Impact Factors

Organophosphate tri- and diesters (tri-OPEs and di-OPEs) were quantified in 63 paired maternal and cord whole blood samples collected in Hubei, China, in which tri-o-cresyl phosphate (ToCP) was predominant. The transplacental transfer efficiencies (expressed as cord blood to maternal blood (C:M) concentration ratios) of aryl-tri-OPEs, such as ToCP (1.61) and triphenyl phosphate (TPHP) (1.06), were higher than those of alkyl-tri-OPEs (0.66-0.76). For the target tri-OPEs and some traditional organic compounds, the C:M ratios first increased with log K_ow in the range of 1.63-5.23 and then decreased, showing a parabolic relationship. However, ToCP, with a log K_ow of 6.34, deviated from this relationship and displayed the highest C:M ratio (1.61). Molecular docking indicated a very strong binding affinity between ToCP and transthyretin, suggesting that ToCP might be actively transported by transthyretin in the placenta. The di-OPE levels in the blood samples were significantly lower than the corresponding tri-OPE levels, and those in the cord blood were influenced not only by their transplacental behaviors but also by their low excretion rates and the metabolic characteristics of their parent compounds in the fetus. This study provides useful information for accurately assessing the health risks posed by tri-OPEs to pregnant women and fetuses.

Enzymes, Reacting with Organophosphorus Compounds as Detoxifiers: Diversity and Functions

Organophosphorus compounds (OPCs) are able to interact with various biological targets in living organisms, including enzymes. The binding of OPCs to enzymes does not always lead to negative consequences for the body itself, since there are a lot of natural biocatalysts that can catalyze the chemical transformations of the OPCs via hydrolysis or oxidation/reduction and thereby provide their detoxification. Some of these enzymes, their structural differences and identity, mechanisms, and specificity of catalytic action are discussed in this work, including results of computational modeling. Phylogenetic analysis of these diverse enzymes was specially realized for this review to emphasize a great area for future development(s) and applications.