Carcinogenic Metabolic Activation Process of Naphthalene by the Cytochrome P450 Enzyme 1B1: A Computational Study

Differential eco-toxicological responses toward Eisenia fetida exposed to soil contaminated with naphthalene and typical metabolites

Naphthalene (NAP) was frequently detected in polycyclic aromatic hydrocarbons (PAHs)-contaminated soil, and its residues may pose an eco-toxicological threat to soil organisms. The toxic effects of NAP were closely tied to phenolic and quinone metabolites in biological metabolism. However, the present knowledge concerning the eco-toxicological impacts of NAP metabolites at the animal level is scanty. Here, we assessed the differences in the eco-toxicological responses of Eisenia fetida (E. fetida) in NAP, 1-naphthol (1-NAO) or 1,4-naphthoquinone (1,4-NQ) contaminated soils. NAP, 1-NAO, and 1,4-NQ exposure triggered the onset of oxidative stress as evidenced by the destruction of the antioxidant enzyme system. The lipid peroxidation and DNA oxidative damage levels induced by 1-NAO and 1,4-NQ were higher than those of NAP. The elevation of DNA damage varied considerably depending on differences in oxidative stress and the direct mode of action of NAP or its metabolites with DNA. All three toxicants induced different degrees of physiological damage to the body wall, but only 1, 4-NQ caused the shedding of intestinal epithelial cells. The integrated biomarker response for different exposure times illustrated that the comprehensive toxicity at the animal level was 1,4-NQ > 1-NAO > NAP, and the time-dependent trends of oxidative stress responses induced by the three toxicants were similar. At the initial stage, the antioxidant system of E. fetida responded positively to the provocation, but the ability of E. fetida to resist stimulation decreased with the prolongation of time resulting in provocation oxidative damage. This study would provide new insights into the toxicological effects and biohazard of PAHs on soil animals.

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.

Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators

We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.

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.

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.

Mutagenicity risk prediction of PAH and derivative mixtures by in silico simulations oriented from CYP compound I-mediated metabolic activation

PAHs and their derivatives are the main sources of mutagenicity and carcinogenicity in airborne particular matter and cause serious public health and environmental problems. Risk assessment is challenging due to the mixed nature and deficiency of toxicity data of most PAHs and their derivatives. Cytochrome P450 enzymes (CYPs) play important roles in PAH-induced carcinogenicity via metabolic activation, and CYP conformations with compound I structures strongly influence metabolic sites and metabolite species. In this study, complexes of BaP with CYP1A1, CYP1B1 or CYP2C19 compound I were successfully simulated by QM/MM methods and verified by metabolic clearance, and the mutagenicity of chemicals was then predicted by the BaP-7,8-epoxide-related metabolic conformation fitness (MCF) approach, which was validated by Ames tests, showing satisfying accuracy (R² = 0.46-0.66). Furthermore, a prediction model of the mutagenicity risk of PAH and derivative mixtures was established based on the relative potential factor (RPF) approach and the RPF calculated from the mathematical relationship between the minimum MCF (MCF_min) and RPF, which was successfully validated by the mutagenesis of PAH and derivative mixture mimic-simulating PM2.5 samples collected in eastern China. This study provides fast reliable tools for assessing risk of the complex components of environmental PAHs and their derivatives.

Insights into the metabolic mechanism of PBDEs catalyzed by cytochrome P450 enzyme 3A4: A QM/MM study

Elucidating the metabolic mechanism and the derivatives of polybrominated diphenyl ethers (PBDEs) is significant to risk assessment. This study delineated the metabolic mechanism of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) catalyzed by P450 enzymes using a combination of molecular dynamic (MD), quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT). The calculation results reveal that the electrophilic addition is the main pathway for the biotransformation of BDE-47 catalyzed by P450 enzymes. 6-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (6-OH-BDE-47) is a more kinetically preferable product than 5-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (5-OH-BDE-47). Electrophilic addition reaction can lead to the formation of polybrominated dibenzo-p-dioxins and dibenzofurans (PBDD/Fs). The ecotoxicity assessment indicates that the final products of BDE-47 are still toxic to aquatic organisms, but the solubility increase of the hydroxylated products can accelerate their excretion from the body. We expect that the established metabolic mechanism and the derivatives will be used to predict the biotransformation of other PBDE congeners catalyzed by P450 enzymes in human livers.

Computational study on the detoxifying mechanism of DDT metabolized by cytochrome P450 enzymes

Predicting the detoxifying mechanism and potential toxic derivatives of xenobiotic substances is significant for risk assessment. The present study delineated the detoxifying mechanism of 1-chloro-4-[2,2,2-trichloro-1-(4-chlorophenyl)ethyl]benzene (DDT) metabolized by human P450 enzymes using a combination of molecular dynamic (MD), quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT). This study highlights that DDT can be metabolized by P450 enzymes through the hydrogen abstraction and electrophilic addition mechanism, and the main derivatives are epoxides (2,3-oxide-DDT and 3,4-oxide-DDT), DDE and dicofol. The epoxides are unstable and the C-O bond cleavage easily occurs by the reaction with hydronium ion or hydroxyl radicals, yielding endocrine disruptor hydroxylated DDT. The eco-toxicity evaluation indicates that the derivatives of DDT are less toxic than DDT, and the solubility increase of the derivatives can accelerate their excretion from the body. The study can provide an understanding of the biotransformation of DDT by the P450 enzymes in human livers.

Multiscale Simulations on the Catalytic Plasticity of CYP76AH1

The catalytic promiscuity and fidelity of cytochrome P450 enzymes are widespread in the skeletal modification of terpenoid natural products and have attracted much attention. CYP76AH1 is involved in key modification reactions in the biosynthetic pathway of tanshinone, a well-known medicinal norditerpenoid. In this work, classical molecular dynamic simulations, metadynamics, and DFT calculations were performed to investigate the protein conformational dynamics, ligand binding poses, and catalytic reaction mechanism in wide-type and mutant CYP76AH1. Our results not only reveal a plausible enzymatic mechanism for mutant CYP76AH1 leading to various products but also provide valuable guidance for rational protein engineering of the CYP76 family.

Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences

Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.

Cu(I)-Catalyzed Ligand-Free Tandem One-Pot or Sequential Annulation via Knoevenagel Intermediates: An Entry into Multifunctional Naphthalenes, Phenanthrenes, Quinolines, and Benzo[b]carbazoles

A simple but efficient one-pot or sequential copper-catalyzed protocol using 2-bromoaldehydes and active methylene group containing substrates that affords multifunctional naphthalenes, phenanthrenes, quinolines, and benzo[b]carbazoles via Knoevenagel condensation, C-arylation, and decarboxylation, followed by aromatization, is developed. The reaction utilizes the potential of Knoevenagel intermediates and does not require any ancillary ligand. The phenanthrene products thus obtained show moderate fluorescence activity. Structural elaboration of the products to obtain dihydrobenzoquinazolines is also highlighted.