Computational study on the metabolic activation mechanism of PeCDD by Cytochrome P450 1A1

A down-regulated cytochrome P450 in Neoseiulus barkeri Hughes (Acari: Phytoseiidae) can dechlorinate and hydroxylate chlorpyrifos without producing chlorpyrifos-oxon

Selection of chemical-resistant predatory mites is a good alternative to balance the contradiction between chemical control and biological control. Previously, a resistant strain of Neoseiulus barkeri for chlorpyrifos was obtained. In the current study, two up-regulated (NbCYP3A6, NbCYP3A16) and one down-regulated (NbCYP3A24) P450s were screened through differential expression analysis and other detoxification-related genes such as CCEs, GST, etc. were not found. 3D modelling and molecular docking indicated that the chlorine at position 5 on the pyridine ring of chlorpyrifos, as well as a methyl group, were closest to the heme iron of the enzymes (less than 5 Å). Three active recombinant P450 proteins were heterologously expressed and metabolized with chlorpyrifos in vitro. HPLC assay showed that only NbCYP3A24 could metabolize chlorpyrifos, with a metabolism rate of 21.60 %. Analysis of the m/z of metabolites by LC-MS/MS showed that chlorine at the 5C position of chlorpyrifos was stripped and hydroxylated, whereas chlorpyrifos-oxon, a common product of oxidation by P450, was not found. Knockdown of the NbCYP3A24 gene in the susceptiblestrain did reduce the susceptibility of N. barkeri to chlorpyrifos, suggesting that the biological activity of the metabolite may be similar to chlorpyrifos-oxon, thus enhancing the inhibitory effect on the target.

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.

Computational study on the endocrine-disrupting metabolic activation of Benzophenone-3 catalyzed by cytochrome P450 1A1: A QM/MM approach

Predicting the metabolic activation mechanism and potential hazardous metabolites of environmental endocrine-disruptors is a challenging and significant task in risk assessment. Here the metabolic activation mechanism of benzophenone-3 catalyzed by P450 1A1 was investigated by using Molecular Dynamics, Quantum Mechanics/Molecular Mechanics and Density Functional Theory approaches. Two elementary reactions involved in the metabolic activation of BP-3 with P450 1A1: electrophilic addition and hydrogen abstraction reactions were both discussed. Further conversion reactions of epoxidation products, ketone products and the formaldehyde formation reaction were investigated in the non-enzymatic environment based on previous experimental reports. Binding affinities analysis of benzophenone-3 and its metabolites to sex hormone binding globulin indirectly demonstrates that they all exhibit endocrine-disrupting property. Toxic analysis shows that the eco-toxicity and bioaccumulation values of the benzophenone-3 metabolites are much lower than those of benzophenone-3. However, the metabolites are found to have skin-sensitization effects. The present study provides a deep insight into the biotransformation process of benzophenone-3 catalyzed by P450 1A1 and alerts us to pay attention to the adverse effects of benzophenone-3 and its metabolites in human livers.

Toxicity assessment of dioxins and their transformation by-products from inferred degradation pathways

Due to the significant POPs characteristics, dioxins caused concern in public health and environmental protection. Evaluating the toxicity risk of dioxin degradation pathways is critical. OCDD, 1,2,3,4,6,7,8-HpCDD, and 1,2,3,4,6,7,8-HpCDF, which are highly abundant in the environment and have strong biodegradation capabilities, were selected as precursor molecules in this study. Firstly, their transformation pathways were deduced during the metabolism of biometabolism, microbial aerobic, microbial anaerobic, and photodegradation pathways, and density function theory (DFT) was used to calculate the Gibbs free energy to infer the possibility of the occurrence of the transformation pathway. Secondly, the carcinogenic potential of the precursors and their degradation products was evaluated using the TOPKAT modeling method. With the help of the positive indicator (0-1) normalization method and heat map analysis, a significant increase in the toxic effect of some of the transformation products was found, and it was inferred that it was related to the structure of the transformation products. Meanwhile, the strength of the endocrine disrupting effect of dioxin transformation products was quantitatively assessed using molecular docking and subjective assignment methods, and it was found that dioxin transformation products with a higher content of chlorine atoms and molecules similar to those of thyroid hormones exhibited a higher risk of endocrine disruption. Finally, the environmental health risks caused by each degradation pathway were comprehensively assessed with the help of the negative indicator (1-2) standardization method, which provides a theoretical basis for avoiding the toxicity risks caused by dioxin degradation transformation. In addition, the 3D-QSAR model was used to verify the necessity and rationality of this study. This paper provides theoretical support and reference significance for the toxicity assessment of dioxin degradation by-products from inferred degradation pathways.

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.

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

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

Endocrine-disrupting metabolic activation of 2-nitrofluorene catalyzed by human cytochrome P450 1A1: A QM/MM approach

Nitropolycyclic aromatic hydrocarbons (NPAHs) present one of the most important airborne pollutants. Recent studies have shown that one of the most abundant NPAHs, 2-Nitrofluorene (NF), was supposed to be converted to endocrine-disrupting metabolites by cytochrome P450 1A1 (CYP1A1) in human cells. However, the mechanism is still largely unexplored. Here the metabolic activation and transformation mechanism of NF catalyzed by CYP1A1 were systematically studied with the aid of Molecular Dynamics, Density Functional Theory and Quantum Mechanics/Molecular Mechanics techniques. We evidence that CYP1A1 can activate NF through two elementary processes: (i) electrophilic addition (12.4 kcal·mol-1) or hydrogen abstraction (38.2 kcal·mol-1) and (ii) epoxidation (5.9 and 8.7 kcal·mol-1) or NIH shift (12.5 and 14.9 kcal·mol-1) or proton shuttle (12.1 kcal·mol-1). Electrophilic addition was found to be the rate-determining step while epoxidation rather than NIH shift or proton shuttle is the more feasible pathway after electrophilic addition. Metabolites 6,7-epoxide-2-nitrofluorene and 7,8-epoxide-2-nitrofluorene were identified as the major epoxidation products. Epoxides are unstable and easy to react with hydrated hydrogen ions and hydroxyls to produce endocrine disrupter 7-hydroxy-2-nitrofluorene. Toxic analysis shows that some of the metabolites are more toxic to model aquatic organisms (e.g. Green algea) than NF. Binding affinity analysis to human sex hormone binding globulin reveals that NF metabolites all have endocrine-disrupting potential. This study provides a comprehensive understanding on the biotransformation process of NF and may aid future studies on various NPAHs activation catalyzed by human P450 enzyme.

Human cytochrome P450 3A-mediated two-step oxidation metabolism of dimethomorph: Implications in the mechanism-based enzyme inactivation

Dimethomorph (DMM), an effective and broad-spectrum fungicide applied in agriculture, is toxic to environments and living organisms due to the hazardous nature of its toxic residues. This study aims to investigate the human cytochrome P450 enzyme (CYP)-mediated oxidative metabolism of DMM by combining experimental and computational approaches. Dimethomorph was metabolized predominantly through a two-step oxidation process mediated by CYPs, and CYP3A was identified as the major contributor to DMM sequential oxidative metabolism. Meanwhile, DMM elicited the mechanism-based inactivation (MBI) of CYP3A in a suicide manner, and the iminium ion and epoxide reactive intermediates generated in DMM metabolism were identified as the culprits of MBI. Furthermore, three common pesticides, prochloraz (PCZ), difenoconazole (DFZ) and chlorothalonil (CTL), could significantly inhibit CYP3A-mediated DMM metabolism, and consequently trigger elevated exposure to DMM in vivo. Computational studies elucidated that the differentiation effects in charge distribution and the interaction pattern played crucial roles in DMM-induced MBI of CYP3A4 during sequential oxidative metabolism. Collectively, this study provided a global view of the two-step metabolic activation process of DMM mediated by CYP3A, which was beneficial for elucidating the environmental fate and toxicological mechanism of DMM in humans from a new perspective.

Mechanism of Sugar Ring Contraction and Closure Catalyzed by UDP-d-apiose/UDP-d-xylose Synthase (UAXS)

Uridine diphosphate (UDP)-apiose/UDP-xylose synthase (UAXS) is a member of the short-chain dehydrogenase/reductase superfamily (SDR), which catalyzes the ring contraction and closure of UDP-d-glucuronic acid (UDP-GlcA), affording UDP-apiose and UDP-xylose. UAXS is a special enzyme that integrates ring-opening, decarboxylation, rearrangement, and ring closure/contraction in a single active site. Recently, the ternary complex structure of UAXS was crystallized from Arabidopsis thaliana. In this work, to gain insights into the detailed formation mechanism of UDP-apiose and UDP-xylose, an enzyme-substrate reactant model has been constructed and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that the reaction starts from the C4-OH oxidation, which is accompanied by the conformational transformation of the sugar ring from chair type to boat type. The sugar ring-opening is prior to decarboxylation, and the deprotonation of the C2-OH group is the prerequisite for sugar ring-opening. Moreover, the keto-enol tautomerization of the decarboxylated intermediate is a necessary step for ring closure/contraction. Based on our calculation results, more UDP-apiose product was expected, which is in line with the experimental observation. Three titratable residues, Tyr185, Cys100, and Cys140, steer the reaction by proton transfer from or to UDP-GlcA, and Arg182, Glu141, and D337 constitute a proton conduit for sugar C2-OH deprotonation. Although Thr139 and Tyr105 are not directly involved in the enzymatic reaction, they are responsible for promoting the catalysis by forming hydrogen-bonding interactions with GlcA. Our calculations may provide useful information for understanding the catalysis of the SDR family.