Oxidized metabolites from cyclopenta-fused polycyclic aromatic hydrocarbons (CP-PAHs). A DFT model study of their carbocations formed by epoxide ring opening

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.

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

The metabolic activation and transformation of naphthalene by the cytochrome P450 enzyme (CYP 1B1) plays an important role in its potential carcinogenicity. The process has been explored by a quantum mechanics/molecular mechanics (QM/MM) computational method. Molecular dynamic simulations were performed to explore the interaction between naphthalene and CYP 1B1. Naphthalene involves α- and β-carbon, the electrophilic addition of which would result in different reaction pathways. Our computational results show that both additions on α- and β-carbon can generate naphthalene 1,2-oxide. The activation barrier for the addition on β-carbon is higher than that for the α-carbon by 2.6 kcal·mol^-1, which is possibly caused by the proximity between β-carbon and the iron-oxo group of Cpd I in the system. We also found that naphthalene 1,2-oxide is unstable and the O-C bond cleavage easily occurs via cellular hydronium ion, hydroxyl radical/anion; then it will convert to the potential ultimate carcinogen 1,2-naphthoquinone. The results demonstrate and inform a detailed process of generating naphthalene 1,2-oxide and new predictions for its conversion.

Computational toxicological investigation on the mechanism and pathways of xenobiotics metabolized by cytochrome P450: a case of BDE-47

Understanding the transformation mechanism and products of xenobiotics catalyzed by cytochrome P450 enzymes (CYPs) is vital to risk assessment. By density functional theory computation with the B3LYP functional, we simulated the reaction of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) catalyzed by the active species of CYPs (Compound I). The enzymatic and aqueous environments were simulated by the polarizable continuum model. The results reveal that the addition of Compound I to BDE-47 is the rate-determining step. The addition of Compound I to the ipso and nonsubstituted C atoms forms tetrahedral σ-adducts that further transform into epoxides. Hydroxylation of the epoxides leads to hydroxylated polybrominated diphenyl ethers and 2,4-dibromophenol. The addition to the Br-substituted C2 and C4 atoms has a higher barrier than addition to the nonsubstituted C atoms, forming phenoxide and cyclohexadienone which subsequently undergo debromination/hydroxylation. A novel mechanism was identified in which the approach of Compound I to C2 led to formation of a phenoxide and an expelled Br(-) ion. The predicted products were consistent with the metabolites identified by others. As a first attempt to simulate the enzymatic transformation of a polycyclic compound, this study may enlighten a computational method to predict the biotransformation of xenobiotics catalyzed by CYPs.