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

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.

Nitro-PAHs: Occurrences, ecological consequences, and remediation strategies for environmental restoration

Nitrated polycyclic aromatic hydrocarbons (nitro-PAHs) are persistent pollutants that have been introduced into the environment as a result of human activities. They are produced when PAHs undergo oxidation and are highly resistant to degradation, resulting in prolonged exposure and significant health risks for wildlife and humans. Nitro-PAHs' potential to induce cancer and mutations has raised concerns about their harmful effects. Furthermore, their ability to accumulate in the food chain seriously threatens the ecosystem and human health. Moreover, nitro-PAHs can disrupt the normal functioning of the endocrine system, leading to reproductive and developmental problems in humans and other organisms. Reducing nitro-PAHs in the environment through source management, physical removal, and chemical treatment is essential to mitigate the associated environmental and human health risks. Recent studies have focused on improving nitro-PAHs' phytoremediation by incorporating microorganisms and biostimulants. Microbes can break down nitro-PAHs into less harmful substances, while biostimulants can enhance plant growth and metabolic activity. By combining these elements, the effectiveness of phytoremediation for nitro-PAHs can be increased. This study aimed to investigate the impact of introducing microbial and biostimulant agents on the phytoremediation process for nitro-PAHs and identify potential solutions for addressing the environmental risks associated with these pollutants.

Pyrene and its derivatives increase lung adverse effects by activating aryl hydrocarbon receptor transcription

Derivatives of polycyclic aromatic hydrocarbons (PAHs) pose significant threat to environment and human health due to their widespread and potential hazards. However, adverse effects and action mechanisms of PAH derivatives on human health have not been attempted yet. Herein, we chose pyrene and its derivatives (1-hydroxypyrene, 1-nitropyrene, and 1-methylpyrene) to investigate adverse effect mechanism to human lungs using in vitro and in vivo methods. Results showed that pyrene derivatives have higher lung health risks than original pyrene. They can activate AhR, subsequently affecting expression of downstream target genes CYP1A1 and CYP1B1. The binding energies of pyrene and its derivatives ranged from -16.07 to -27.25 kcal/mol by molecular dynamics simulations, implying that pyrene and its derivatives acted as agonists of AhR and increased adverse effects on lungs. Specifically, 1-nitropyrene exhibited stabler binding conformation and stronger AhR expression. In addition, sensitivity of pyrene and its derivatives to AhR activation was attributed to type and number of key amino acids in AhR, that is, pyrene (Leu293), 1-nitropyrene (Cys333, Met348, and Val381), 1-hydroxypyrene (Leu293 and Phe287), and 1-methylpyrene (Met348). In summary, we provide a universal approach for understanding action mechanisms of PAH derivatives on human health, and their adverse effects should be taken seriously.

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.

New insight into biodegradation mechanism of phenylurea herbicides by cytochrome P450 enzymes: Successive N-demethylation mechanism

Phenylurea herbicides (PUHs) present one of the most important herbicides, which have cause serious effects on ecological environment and humans. Nowadays enzyme strategy shows great advantages in degradation of PUHs. Here density functional theory (DFT), quantitative structure - activity relationship (QSAR) and quantum mechanics/molecular mechanics (QM/MM) approaches are used to investigate the degradation mechanism of PUHs catalyzed by P450 enzymes. Two successive N-demethylation pathways are identified and two hydrogen abstraction (H-abstraction) reaction pathways are identified as the rate-determining step through high-throughput DFT calculations. The Boltzmann-weighted average energy barrier of the second H-abstraction pathway (19.95 kcal/mol) is higher than that of the first H-abstraction pathway (16.80 kcal/mol). Two QSAR models are established to predict the energy barriers of the two H-abstraction pathways based on the quantum chemical descriptors and mordred molecular descriptors. The determination coefficient (R2) values of QSAR models are > 0.9, which reveal that the established QSAR models have great predictive capability. QM/MM calculations indicate that human P450 enzymes are more efficient in degradation of PUHs than crop and weed P450 enzymes. Correlations between energy barriers and key structural/charge parameters are revealed and key parameters that have influence on degradation efficiency of PUHs are identified. This study provides lateral insights into the biodegradation strategy and removal method of PUHs and valuable information for designing or engineering of highly efficient degradation enzymes and genetically modified crops.

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.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

Transformation and toxicity dynamics of polycyclic aromatic hydrocarbons in a novel biological-constructed wetland-microalgal wastewater treatment process

In this study, a novel wastewater treatment process combining sequencing batch reactor, constructed wetland and microalgal membrane photobioreactor (BCM process) was proposed, and its performance on removal, transformation and toxicity reduction of polycyclic aromatic hydrocarbons (PAHs) was intensively explored. Satisfactory PAHs removal (90.58%-97.50%) was achieved and molecular weight had significant impact on the removal pathways of different PAHs. Adsorption dominated the removal of high molecular weight PAHs, while the contribution ratio of microbial degradation increased with the decrease of molecular weight of PAHs. More importantly, it was reported for the first time that substituted PAHs (SPAHs) produced by microbial degradation of PAHs would lead to increased toxicity during the BCM process. High PAHs (75.37%-88.52%) and SPAHs removal (99.56%-100.00%) were achieved in the microalgae unit due to its abundant cytochrome P450 enzyme, which decreased the bacterial toxicity by 90.93% and genotoxicity by 93.08%, indicating that microalgae played significance important role in ensuring water security. In addition, the high quantitative relationship (R² = 0.98) between PAHs, SPAHs and toxicity exhibited by regression model analysis proved that more attention should be paid to the ecotoxicity of derivatives of refractory organic matters in wastewater treatment plants.