Jessica H. Hartman

CYP2E1 metabolism of styrene involves allostery.

We are the first to report allosterism during styrene oxidation by recombinant CYP2E1 and human liver microsomes. At low styrene concentrations, oxidation is inefficient because of weak binding to CYP2E1 (Ks = 830 μM). A second styrene molecule then binds CYP2E1 with higher affinity (Kss = 110 μM) and significantly improves oxidation to achieve a kcat of 6.3 nmol · min−1 · nmol CYP2E1−1. The transition between these metabolic cycles coincides with reported styrene concentrations in blood from exposed workers; thus, this CYP2E1 mechanism may be relevant in vivo. Scaled modeling of the in vitro-positive allosteric mechanism for styrene metabolism to its in vivo clearance led to significant deviations from the traditional model based on Michaelis-Menten kinetics. Low styrene levels were notably much less toxic than generally assumed. We interrogated the allosteric mechanism using the CYP2E1-specific inhibitor and drug 4-methylpyrazole, which we have shown binds two CYP2E1 sites. From the current studies, styrene was a positive allosteric effector on 4-methylpyrazole binding, based on a 10-fold increase in 4-methylpyrazole binding affinity from Ki 0.51 to Ksi 0.043 μM. The inhibitor was a negative allosteric effector on styrene oxidation, because kcat decreased 6-fold to 0.98 nmol · min−1 · nmol CYP2E1−1. Consequently, mixtures of styrene and other molecules can induce allosteric effects on binding and metabolism by CYP2E1 and thus mitigate the efficiency of their metabolism and corresponding effects on human health. Taken together, our elucidation of mechanisms for these allosteric reactions provides a powerful tool for further investigating the complexities of CYP2E1 metabolism of drugs and pollutants.

Structure of pyrazole derivatives impact their affinity, stoichiometry, and cooperative interactions for CYP2E1 complexes.

CYP2E1 plays a critical role in detoxification and carcinogenic activation of drugs, pollutants, and dietary compounds; however, these metabolic processes can involve poorly characterized cooperative interactions that compromise the ability to understand and predict CYP2E1 metabolism. Herein, we employed an array of ten azoles with an emphasis on pyrazoles to establish the selectivity of catalytic and cooperative CYP2E1 sites through binding and catalytic studies. Spectral binding studies for monocyclic azoles suggested two binding events, while bicyclic azoles suggested one. Pyrazole had moderate affinity toward the CYP2E1 catalytic site that improved when a methyl group was introduced at either position 3 or 4. The presence of methyl groups simultaneously at positions 3 and 5 blocked binding, and a phenyl group at position 3 did not improve binding affinity. In contrast, pyrazole fusion to a benzene or cyclohexane ring greatly increased affinity. The consequences of these binding events on CYP2E1 catalysis were studied through inhibition studies with 4-nitrophenol, a substrate known to bind both sites. Most pyrazoles shared a common mixed cooperative inhibition mechanism in which pyrazole binding rescued CYP2E1 from substrate inhibition. Overall, inhibitor affinities toward the CYP2E1 catalytic site were similar to those reported in binding studies, and the same trend was observed for binding at the cooperative site. Taken together, these studies identified key structural determinants in the affinity and stoichiometry of azole interactions with CYP2E1 and consequences on catalysis that further advance an understanding of the relationship between structure and function for this enzyme.

CYP2E1 hydroxylation of aniline involves negative cooperativity.

CYP2E1 plays a role in the metabolic activation and elimination of aniline, yet there are conflicting reports on its mechanism of action, and hence relevance, in aniline metabolism. Based on our work with similar compounds, we hypothesized that aniline binds two CYP2E1 sites during metabolism resulting in cooperative reaction kinetics and tested this hypothesis through rigorous in vitro studies. The kinetic profile for recombinant CYP2E1 demonstrated significant negative cooperativity based on a fit of data to the Hill equation (n=0.56). Mechanistically, the data were best explained through a two-binding site cooperative model in which aniline binds with high affinity (K(s)=30 μM) followed by a second weaker binding event (K(ss)=1100 uM) resulting in a threefold increase in the oxidation rate. Binding sites for aniline were confirmed by inhibition studies with 4-methylpyrazole. Inhibitor phenotyping experiments with human liver microsomes validated the central role for CYP2E1 in aniline hydroxylation and indicated minor roles for CYP2A6 and CYP2C9. Importantly, inhibition of minor metabolic pathways resulted in a kinetic profile for microsomal CYP2E1 that replicated the preferred mechanism and parameters observed with the recombinant enzyme. Scaled modeling of in vitro CYP2E1 metabolism of aniline to in vivo clearance, especially at low aniline levels, led to significant deviations from the traditional model based on non-cooperative, Michaelis-Menten kinetics. These findings provide a critical mechanistic perspective on the potential importance of CYP2E1 in the metabolic activation and elimination of aniline as well as the first experimental evidence of a negatively cooperative metabolic reaction catalyzed by CYP2E1.

Structural basis for cooperative binding of azoles to CYP2E1 as interpreted through guided molecular dynamics simulations.

CYP2E1 metabolizes a wide array of small, hydrophobic molecules, resulting in their detoxification or activation into carcinogens through Michaelis-Menten as well as cooperative mechanisms. Nevertheless, the molecular determinants for CYP2E1 specificity and metabolic efficiency toward these compounds are still unknown. Herein, we employed computational docking studies coupled to molecular dynamics simulations to provide a critical perspective for understanding a structural basis for cooperativity observed for an array of azoles from our previous binding and catalytic studies (Hartman et al., 2014). The resulting 28 CYP2E1 complexes in this study revealed a common passageway for azoles that included a hydrophobic steric barrier causing a pause in movement toward the active site. The entrance to the active site acted like a second sieve to restrict access to the inner chamber. Collectively, these interactions impacted the final orientation of azoles reaching the active site and hence could explain differences in their biochemical properties observed in our previous studies, such as the consequences of methylation at position 5 of the azole ring. The association of a second azole demonstrated significant differences in interactions stabilizing the bound complex than observed for the first binding event. Intermolecular interactions occurred between the two azoles as well as CYP2E1 residue side chains and backbone and involved both hydrophobic contacts and hydrogen bonds. The relative importance of these interactions depended on the structure of the respective azoles indicating the absence of specific defining criteria for binding unlike the well-characterized dominant role of hydrophobicity in active site binding. Consequently, the structure activity relationships described here and elsewhere are necessary to more accurately identify factors impacting the observation and significance of cooperativity in CYP2E1 binding and catalysis toward drugs, dietary compounds, and pollutants.

Subcellular localization of rat CYP2E1 impacts metabolic efficiency toward common substrates

Cytochrome P450 2E1 (CYP2E1) detoxifies or bioactivates many low molecular-weight compounds. Most knowledge about CYP2E1 activity relies on studies of the enzyme localized to endoplasmic reticulum (erCYP2E1); however, CYP2E1 undergoes transport to mitochondria (mtCYP2E1) and becomes metabolically active. We report the first comparison of in vitro steady-state kinetic profiles for erCYP2E1 and mtCYP2E1 oxidation of probe substrate 4-nitrophenol and pollutants styrene and aniline using subcellular fractions from rat liver. For all substrates, metabolic efficiency changed with substrate concentration for erCYP2E1 reflected in non-hyperbolic kinetic profiles but not for mtCYP2E1. Hyperbolic kinetic profiles for the mitochondrial enzyme were consistent with Michaelis-Menten mechanism in which metabolic efficiency was constant. By contrast, erCYP2E1 metabolism of 4-nitrophenol led to a loss of enzyme efficiency at high substrate concentrations when substrate inhibited the reaction. Similarly, aniline metabolism by erCYP2E1 demonstrated negative cooperativity as metabolic efficiency decreased with increasing substrate concentration. The opposite was observed for erCYP2E1 oxidation of styrene; the sigmoidal kinetic profile indicated increased efficiency at higher substrate concentrations. These mechanisms and CYP2E1 levels in mitochondria and endoplasmic reticulum were used to estimate the impact of CYP2E1 subcellular localization on metabolic flux of pollutants. Those models showed that erCYP2E1 mainly carries out aniline metabolism at all aniline concentrations. Conversely, mtCYP2E1 dominates styrene oxidation at low styrene concentrations and erCYP2E1 at higher concentrations. Taken together, subcellular localization of CYP2E1 results in distinctly different enzyme activities that could impact overall metabolic clearance and/or activation of substrates and thus impact the interpretation and prediction of toxicological outcomes.

Inhibitory potency of 4-carbon alkanes and alkenes toward CYP2E1 activity.

CYP2E1 has been implicated in the bioactivation of many small molecules into reactive metabolites which form adducts with proteins and DNA, and thus a better understanding of the molecular determinants of its selectivity are critical for accurate toxicological predictions. In this study, we determined the potency of inhibition of human CYP2E1 for various 4-carbon alkanes, alkenes and alcohols. In addition, known CYP2E1 substrates and inhibitors including 4-methylpyrazole, aniline, and dimethylnitrosamine were included to determine their relative potencies. Of the 1,3-butadiene-derived metabolites studied, 3,4-epoxy-1-butene was the strongest inhibitor with an IC50 of 110 μM compared to 1700 μM and 6600 μM for 1,2-butenediol and 1,2:3,4-diepoxybutane, respectively. Compared to known inhibitors, inhibitory potency of 3,4-epoxy-1-butene is between 4-methylpyrazole (IC50 = 1.8 μM) and dimethylnitrosamine (IC50 = 230 μM). All three butadiene metabolites inhibit CYP2E1 activity through a simple competitive mechanism. Among the 4-carbon compounds studied, the presence and location of polar groups seems to influence inhibitory potency. To further examine this notion, the investigation was extended to include structurally and chemically similar analogues, including propylene oxide and various butane alcohols. Those results demonstrated preferential recognition of CYP2E1 toward the type and location of polar and hydrophobic structural elements. Taken together, CYP2E1 metabolism may be modified in vivo by exposure to 4-carbon compounds, such as drugs, and nutritional constituents, a finding that highlights the complexity of exposure to mixtures.

Cooperative effects for CYP2E1 differ between styrene and its metabolites.

Abstract 1. Cooperative interactions are frequently observed in the metabolism of drugs and pollutants by cytochrome P450s; nevertheless, the molecular determinants for cooperativity remain elusive. Previously, we demonstrated that steady-state styrene metabolism by CYP2E1 exhibits positive cooperativity. 2. We hypothesized that styrene metabolites have lower affinity than styrene toward CYP2E1 and limited ability to induce cooperative effects during metabolism. To test the hypothesis, we determined the potency and mechanism of inhibition for styrene and its metabolites toward oxidation of 4-nitrophenol using CYP2E1 Supersomes® and human liver microsomes. 3. Styrene inhibited the reaction through a mixed cooperative mechanism with high affinity for the catalytic site (67 µM) and lower affinity for the cooperative site (1100 µM), while increasing substrate turnover at high concentrations. Styrene oxide and 4-vinylphenol possessed similar affinity for CYP2E1. Styrene oxide behaved cooperatively like styrene, but 4-vinylphenol decreased turnover at high concentrations. Styrene glycol was a very poor competitive inhibitor. Among all compounds, there was a positive correlation with binding and hydrophobicity. 4. Taken together, these findings for CYP2E1 further validate contributions of cooperative mechanisms to metabolic processes, demonstrate the role of molecular structure on those mechanisms and underscore the potential for heterotropic cooperative effects between different compounds.

1,3-Butadiene-induced mitochondrial dysfunction is correlated with mitochondrial CYP2E1 activity in Collaborative Cross mice

Cytochrome P450 2E1 (CYP2E1) metabolizes low molecular weight hydrophobic compounds, including 1,3-butadiene, which is converted by CYP2E1 to electrophilic epoxide metabolites that covalently modify cellular proteins and DNA. Previous CYP2E1 studies have mainly focused on the enzyme localized in the endoplasmic reticulum (erCYP2E1); however, active CYP2E1 has also been found in mitochondria (mtCYP2E1) and the distribution of CYP2E1 between organelles can influence an individuals response to exposure. Relatively few studies have focused on the contribution of mtCYP2E1 to activation of chemical toxicants. We hypothesized that CYP2E1 bioactivation of 1,3-butadiene within mitochondria adversely affects mitochondrial respiratory complexes I-IV. A population of Collaborative Cross mice was exposed to air (control) or 200ppm 1,3-butadiene. Subcellular fractions (mitochondria, DNA, and microsomes) were collected from frozen livers and CYP2E1 activity was measured in microsomes and mitochondria. Individual activities of mitochondrial respiratory complexes I-IV were measured using in vitro assays and purified mitochondrial fractions. In air- and 1,3-butadiene-exposed mouse samples, mtDNA copy numbers were assessed by RT-PCR, and mtDNA integrity was assessed through a PCR-based assay. No significant changes in mtDNA copy number or integrity were observed; however, there was a decrease in overall activity of mitochondrial respiratory complexes I, II, and IV after 1,3-butadiene exposure. Additionally, higher mtCYP2E1 (but not erCYP2E1) activity was correlated with decreased mitochondrial respiratory complex activity (in complexes I-IV) in the 1,3-butadiene-exposed (not control) animals. Together, these results represent the first in vivo link between mitochondrial CYP2E1 activity and mitochondrial toxicity.

Cooperativity in CYP2E1 metabolism of acetaminophen and styrene mixtures.

Risk assessment for exposure to mixtures of drugs and pollutants relies heavily on in vitro characterization of their bioactivation and/or metabolism individually and extrapolation to mixtures assuming no interaction. Herein, we demonstrated that in vitro CYP2E1 metabolic activation of acetaminophen and styrene mixtures could not be explained through the Michaelis-Menten mechanism or any models relying on that premise. As a baseline for mixture studies with styrene, steady-state analysis of acetaminophen oxidation revealed a biphasic kinetic profile that was best described by negative cooperativity (Hill coefficient=0.72). The best-fit mechanism for this relationship involved two binding sites with differing affinities (Ks=830μM and Kss=32mM). Introduction of styrene inhibited that reaction less than predicted by simple competition and thus provided evidence for a cooperative mechanism within the mixture. Likewise, acetaminophen acted through a mixed-type inhibition mechanism to impact styrene epoxidation. In this case, acetaminophen competed with styrene for CYP2E1 (Ki=830μM and Ksi=180μM for catalytic and effector sites, respectively) and resulted in cooperative impacts on binding and catalysis. Based on modeling of in vivo clearance, cooperative interactions between acetaminophen and styrene resulted in profoundly increased styrene activation at low styrene exposure levels and therapeutic acetaminophen levels. Current Michaelis-Menten based toxicological models for mixtures such as styrene and acetaminophen would fail to detect this concentration-dependent relationship. Hence, future studies must assess the role of alternate CYP2E1 mechanisms in bioactivation of compounds to improve the accuracy of interpretations and predictions of toxicity.

Predicting CYP2C19 catalytic parameters for enantioselective oxidations using artificial neural networks and a chirality code.

Cytochromes P450 (CYP for isoforms) play a central role in biological processes especially metabolism of chiral molecules; thus, development of computational methods to predict parameters for chiral reactions is important for advancing this field. In this study, we identified the most optimal artificial neural networks using conformation-independent chirality codes to predict CYP2C19 catalytic parameters for enantioselective reactions. Optimization of the neural networks required identifying the most suitable representation of structure among a diverse array of training substrates, normalizing distribution of the corresponding catalytic parameters (k(cat), K(m), and k(cat)/K(m)), and determining the best topology for networks to make predictions. Among different structural descriptors, the use of partial atomic charges according to the CHelpG scheme and inclusion of hydrogens yielded the most optimal artificial neural networks. Their training also required resolution of poorly distributed output catalytic parameters using a Box-Cox transformation. End point leave-one-out cross correlations of the best neural networks revealed that predictions for individual catalytic parameters (k(cat) and K(m)) were more consistent with experimental values than those for catalytic efficiency (k(cat)/K(m)). Lastly, neural networks predicted correctly enantioselectivity and comparable catalytic parameters measured in this study for previously uncharacterized CYP2C19 substrates, R- and S-propranolol. Taken together, these seminal computational studies for CYP2C19 are the first to predict all catalytic parameters for enantioselective reactions using artificial neural networks and thus provide a foundation for expanding the prediction of cytochrome P450 reactions to chiral drugs, pollutants, and other biologically active compounds.