Binding and Metabolism of Brominated Flame Retardant β-1,2-Dibromo-4-(1,2-dibromoethyl)cyclohexane in Human Microsomal P450 Enzymes: Insights from Computational Studies

Computational simulations uncover enantioselective metabolism of chiral triazole fungicides by human CYP450 enzymes: A case study of tebuconazole

Tebuconazole (TEB), a prominent chiral triazole fungicide, has been extensively utilized for plant pathogen control globally. Despite experimental evidence of TEB metabolism in mammals, the enantioselectivity in the biotransformation of R- and S-TEB enantiomers by specific CYP450s remains elusive. In this work, integrated in silico simulations were employed to unveil the binding interactions and enantioselective metabolic fate of TEB enantiomers within human CYP1A2, 2B6, 2E1, and 3A4. Molecular dynamics (MD) simulations clearly delineated the binding specificity of R- and S-TEB to the four CYP450s, crucially determining their differences in metabolic activity and enantioselectivity. The primary driving force for robust ligand binding was identified as van der Waals interactions with CYP450s, particularly involving the hydrophobic residues. Mechanistic insights derived from quantum mechanics/molecular mechanics (QM/MM) calculations established C2-methyl hydroxylation as the predominant route of R-/S-TEB metabolism, while C6-hydroxylation and triazol epoxidation were deemed kinetically infeasible pathways. Specifically, the resulting hydroxy-R-TEB metabolite primarily originates from R-TEB biotransformation by 1A2, 2E1 and 3A4, whereas hydroxy-S-TEB is preferentially produced by 2B6. These findings significantly contribute to our comprehension of the binding specificity and enantioselective metabolic fate of chiral TEB by CYP450s, potentially informing further research on human health risk assessment associated with TEB exposure.

Molecular Interactions Governing the Rat Aryl Hydrocarbon Receptor Activities of Polycyclic Aromatic Compounds and Predictive Model Development

Polycyclic aromatic compounds (PACs) exhibit rat aryl hydrocarbon receptor (rAhR) activities, leading to diverse biological or toxic effects. In this study, the key amino residues and molecular interactions that govern the rAhR activity of PACs were investigated using in silico strategies. The homology model of rAhR was first docked with 90 PACs to yield complexes, and the results of the molecular dynamics simulations of 16 typical complexes showed that the binding energies of the complexes range from -7.37 to -26.39 kcal/mol. The major contribution to the molecular interaction comes from van der Waals forces, and Pro295 and Arg316 become the key residues involved in most complexes. Two QSAR models were further developed to predict the rAhR activity of PACs (in terms of log IEQ for PACs without halogen substitutions and log%-TCDD-max for halogenated PACs). Both models have good predictive ability, robustness, and extrapolation ability. Molecular polarizability, electronegativity, size, and nucleophilicity are identified as the important factors affecting the rAhR activity of PACs. The developed models could be employed to predict the rAhR activity of other reactive PACs. This work provides insight into the mechanisms and interactions of the rAhR activity of PACs and assists in the assessment of their fate and risk in organisms.

Mechanism of 7H-Dibenzo[c,g]carbazole metabolism in cytochrome P450 1A1: Insights from computational studies

7H-Dibenzo[c,g]carbazole (DBC) is a prevalent environmental contaminant that induces tumorigenesis in several experimental animals. Recently, it has been utilized to develop high-performance solar cells and organic phosphorescent materials. It is imperative to strengthen investigations of DBC metabolism to understand its potential risks to human health. In this study, human CYP1A1 was employed as the metabolic enzyme to investigate the metabolic mechanism of DBC by molecular docking, molecular dynamics (MD) simulation, and quantum mechanical (QM) calculation. The results indicate that DBC binds to CYP1A1 in two modes (mode 1 and mode 2) mainly through nonpolar solvation energies (ΔGnonpolar). The formation of the two binding modes is attributed to the anchoring effect of the hydrogen bond formed by DBC with Asp320 (mode 1) or Ser116 (mode 2). Mode 1 is a "reactive" conformation, while mode 2 is not considered a "reactive" conformation. C5 is identified as the dominant site, and the pyrrole nitrogen cannot participate in the metabolism. DBC is metabolized mainly by a distinct electrophilic addition-rearrangement mechanism, with an energy barrier of 21.74 kcal/mol. The results provide meaningful insights into the biometabolic process of DBC and contribute to understanding its environmental effects and health risks.

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.

Pharmacokinetic characterisation of a valproate Autism Spectrum Disorder rat model in a context of co-exposure to α-Hexabromocyclododecane

Assessing the role of α-hexabromocyclododecane α-HBCDD as a factor of susceptibility for Autism Spectrum disorders by using valproic acid-exposed rat model (VPA) required characterizing VPA pharmacokinetic in the context of α-HBCDD-co-exposure in non-pregnant and pregnant rats. The animals were exposed to α-HBCDD by gavage (100 ng/kg/day) for 12 days. This was followed by a single intraperitoneal dose of VPA (500 mg/kg) or a daily oral dose of VPA (500 mg/kg) for 3 days. Exposure to α-HBCDD did not affect the pharmacokinetics of VPA in pregnant or non-pregnant rats. Surprisingly, VPA administration altered the pharmacokinetics of α-HBCDD. VPA also triggered higher foetal toxicity and lethality with the PO than IP route. α-HBCDD did not aggravate the embryotoxicity observed with VPA, regardless of the route of exposure. Based on this evidence, a single administration of 500 mg/kg IP is the most suitable VPA model to investigate α-HBCDD co-exposure.

Enantioselective metabolism of novel chiral insecticide Paichongding by human cytochrome P450 3A4: A computational insight

As a novel chiral neonicotinoid insecticide, Paichongding (IPP) has been widely applied in agriculture due to its excellent insecticidal activity. However, the enantioselective metabolism of IPP stereoisomers (5R7R-IPP, 5S7S-IPP, 5R7S-IPP, and 5S7R-IPP) mediated by enzymes in non-target organisms, especially the cytochrome P450s (CYPs), remains unknown. To address this knowledge gap, we developed an integrated computational framework to elucidate the binding interactions and enantioselective metabolism of IPP stereoisomers in human CYP3A4. The results reveal that 5R7R-IPP shows much stronger binding affinity to CYP3A4 than 5S7S-IPP, while enantiomers 5R7S-IPP and 5S7R-IPP have no essential difference in their binding potential, owing to their specific interactions with key CYP3A4 residues. Although enantiomers 5R7R-IPP and 5S7S-IPP feature distinct binding modes resulting from the chiral differences, their transformation activities are slightly different, with C5 and C13 being the primary metabolic sites, respectively. In contrast, CYP3A4 preferably metabolizes 5R7S-IPP over 5S7R-IPP. The metabolism of epimers 5R7R-IPP and 5R7S-IPP share C5-hydroxylation routes due to the conserved 5R-conformaitons, but differ with the transformation routes at C11/C13 and C3 sites. The 7R-chirality of 5S7R-IPP significantly reduces the metabolic potency compared to 5S7S-IPP. CYP3A4-catalyzed hydroxylation and desaturation of IPP stereoisomers generate various chiral metabolites, with C5- and C13-hydroxyIPPs further transforming into depropylated products. Furthermore, the toxicity assessment reveals that IPP, along with the majority of its hydroxylated, desaturated, and depropylated metabolites, can potentially induce adverse effects on human health, specifically hepatotoxicity, respiratory toxicity, and carcinogenicity. This study provides valuable insights into the enantioselective fate of chiral IPP metabolism by CYP3A4, and the identified metabolites can serve as potential biomarkers for monitoring IPP exposure and associated health risk in human body.

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.

In Silico simulation of Cytochrome P450-Mediated metabolism of aromatic amines: A case study of N-Hydroxylation

Aromatic amines, the widely used raw materials in industry, cause long-term exposure to human bodies. They can be metabolized by cytochrome P450 enzymes to form active electrophilic compounds, which will potentially react with nucleophilic DNA to exert carcinogenesis. The short lifetime and versatility of the oxidant (a high-valent iron (IV)-oxo species, compound I) of P450 enzymes prompts us to use theoretical methods to investigate the metabolism of aromatic amines. In this work, the density functional theory (DFT) has been employed to simulate the hydroxylation metabolism through H-abstraction and to calculate the activation energy of this reaction for 28 aromatic amines. The results indicate that the steric effects, inductive effects and conjugative effects greatly contribute to the metabolism activity of the chemicals. The further correlation reveals that the dissociation energy of -NH₂ (BDE_N-H) can successfully predict the time-consuming calculated activation energy (R² for aromatic and heteroaromatic amines are 0.93 and 0.86, respectively), so BDE_N-H can be taken as a key parameter to characterize the relative stability of aromatic amines in P450 enzymes and further to quickly assess their potential toxicity. The validation results prove such relationship has good statistical performance (q_cv² for aromatic and heteroaromatic amines are 0.95 and 0.90, respectively) and can be used to other aromatic amines in the application domain, greatly reducing computational cost and providing useful support for experimental research.

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.

Bioaccumulation and biotransformation of tetrabromoethylcyclohexane (TBECH) in maize (Zea mays L.): Stereoselective driving roles of plant biomacromolecules

The bioaccumulation and biotransformation of tetrabromoethylcyclohexane (TBECH) in maize were investigated. Furthermore, the roles of plant biomacromolecules such as lipid transfer proteins (LTPs), CYP and GST enzymes in driving the biological processes of TBECH stereoisomers were explored. The uptake and translocation of TBECH in maize were diastereo- and enantio-selective. Isomerization from α- to δ-TBECH and β- to γ-TBECH, and metabolites of debromination, hydroxylation and TBECH-GSH adducts were identified in maize roots. The gene expressions of LTPs, CYPs and GSTs were extensively changed in maize after exposure to technical TBECH. CYP and GST enzyme activities as well as GST31 and CYP71C3v2 gene expressions were selectively induced or inhibited by TBECH diastereomers over time. TBECH was able to dock into the active sites and bind with specific residues of the typical biomacromolecules ZmLTP1.6, GST31 and CYP71C3v2, indicating their roles in the bioaccumulation and metabolization of TBECH. Binding modes and affinities to biomacromolecules were significantly different between α- and β-TBECH, which contributed to their stereo-selectivity. This study provided a deep understanding of the biological fate of TBECH, and revealed the driving molecular mechanisms of the selectivity of TBECH stereoisomers in plants.

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 insight into biotransformation of halophenols by cytochrome P450: Mechanism and reactivity for epoxidation

Halophenols (XPs) have aroused great interests due to their high toxicity and low biodegradability. Previous experimental studies have shown that XPs can be catalytically transformed into epoxides and haloquinones by cytochrome P450 enzymes (CYPs). However, these metabolites have never been detected directly. Moreover, the effects of the reaction site and the type and number of halogen substituents on the biotransformation reactivity of halophenols still remain unknown. In this work, we performed density functional theory (DFT) calculations to simulate the CYP-mediated biotransformation of 36 XPs with mono-, di-, and tri-halogen (F, Cl, and Br) substitutions to unravel the mechanism and relevant kinetics of XPs epoxidation. The whole epoxidation process consists of initial rate-determining O-addition and subsequent ring-closure steps. The simulation results show that the epoxidation in low-spin (LS) state is kinetically preferred over that in high-spin (HS) state, and the formation of epoxide metabolite is strongly exothermic. For all XPs, the epoxidation reactivity follows the order of ortho/para O-addition > meta O-addition. Moreover, the O-addition with higher energy barriers roughly corresponds to chlorophenols and fluorophenols with more halogen atoms. Compared with dichlorophenols, the additional ortho-Cl substitution on trichlorophenols can slightly increase the energy barriers of meta O-addition. By contrast, the additional inclusion of an ortho-Cl to monochlorophenols enhances the meta O-addition reactivity of dichlorophenols. Overall, the present work clarifies the biotransformation routes of XPs to produce epoxides, and identifies the key factors affecting the epoxidation reactivity, which are beneficial in understanding comprehensively the metabolic fate and toxicity of XPs.

Sublethal effects of DBE-DBCH diastereomers on physiology, behavior, and gene expression of Daphnia magna

1,2-dibromo-4-(1,2-dibromoethyl)-cyclohexane (DBE-DBCH) is a brominated flame retardant used in commercial and industrial applications. The use of DBE-DBCH containing products has resulted in an increased release into the environment. However, limited information is available on the long-term effects of DBE-DBCH and its effects in aquatic invertebrates. Thus, the present study was aimed at determining how DBE-DBCH diastereomers (αβ and γδ) affects aquatic invertebrates using Daphnia magna as a model organism. Survival, reproduction, feeding, swimming behavior and toxicogenomic responses to environmental relevant concentrations of DBE-DBCH were analyzed. Chronic exposure to DBE-DBCH resulted in decreased lifespan, and reduced fecundity. Expression of genes involved in reproductive processes, vtg1 and jhe, were also inhibited. DBE-DBCH also induced hypoxia by inhibiting the transcription of genes involved in heme biosynthesis and oxygen transport. Furthermore, DBE-DBCH also inhibited feeding resulting in emptiness of the alimentary canal. Increased expression of the stress response biomarkers was observed following DBE-DBCH exposure. In addition, DBE-DBCH diastereomers also altered the swimming behavior of Daphnia magna. The present study demonstrates that DBE-DBCH cause multiple deleterious effects on Daphnia magna, including effects on reproduction and hormonal systems. These endocrine disrupting effects are in agreement with effects observed on vertebrates. Furthermore, as is the case in vertebrates, DBE-DBCH γδ exerted stronger effects than DBE-DBCH αβ on Daphnia magna. This indicate that DBE-DBCH γδ has properties making it more toxic to all so far studied animals than DBE-DBCH αβ.