Drug toxicity is often caused by electrophilic reactive metabolites that covalently bind to proteins. Consequently, the quantitative strength of a molecule’s reactivity with glutathione (GSH) is a frequently used indicator of its toxicity. Through cysteine, GSH (and proteins) scavenges reactive molecules to form conjugates in the body. GSH conjugates to specific atoms in reactive molecules: their sites of reactivity. The value of knowing a molecule’s sites of reactivity is unexplored in the literature. This study tests the value of site of reactivity data that identifies the atoms within 1213 reactive molecules that conjugate to GSH and builds models to predict molecular reactivity with glutathione. An algorithm originally written to model sites of cytochrome P450 metabolism (called XenoSite) finds clear patterns in molecular structure that identify sites of reactivity within reactive molecules with 90.8% accuracy and separate reactive and unreactive molecules with 80.6% accuracy. Furthermore, the model output strongly correlates with quantitative GSH reactivity data in chemically diverse, external data sets. Site of reactivity data is nearly unstudied in the literature prior to our efforts, yet it contains a strong signal for reactivity that can be utilized to more accurately predict molecule reactivity and, eventually, toxicity.

•Basis for azole cooperativity binding and catalytic inhibition of CYP2E1 is unknown.•Guided molecular dynamics simulations yielded 28 CYP2E1–azole complexes.•Active site contacts dominated by common hydrophobic and steric interactions.•Cooperative site interactions are more diverse and dependent on azole structure.•CYP2E1 structure activity relationships revealed challenge to predict cooperativity.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.

Drug toxicity is frequently caused by electrophilic reactive metabolites that covalently bind to proteins. Epoxides comprise a large class of three-membered cyclic ethers. These molecules are electrophilic and typically highly reactive due to ring tension and polarized carbon–oxygen bonds. Epoxides are metabolites often formed by cytochromes P450 acting on aromatic or double bonds. The specific location on a molecule that undergoes epoxidation is its site of epoxidation (SOE). Identifying a molecule’s SOE can aid in interpreting adverse events related to reactive metabolites and direct modification to prevent epoxidation for safer drugs. This study utilized a database of 702 epoxidation reactions to build a model that accurately predicted sites of epoxidation. The foundation for this model was an algorithm originally designed to model sites of cytochromes P450 metabolism (called XenoSite) that was recently applied to model the intrinsic reactivity of diverse molecules with glutathione. This modeling algorithm systematically and quantitatively summarizes the knowledge from hundreds of epoxidation reactions with a deep convolution network. This network makes predictions at both an atom and molecule level. The final epoxidation model constructed with this approach identified SOEs with 94.9% area under the curve (AUC) performance and separated epoxidized and non-epoxidized molecules with 79.3% AUC. Moreover, within epoxidized molecules, the model separated aromatic or double bond SOEs from all other aromatic or double bonds with AUCs of 92.5% and 95.1%, respectively. Finally, the model separated SOEs from sites of sp2 hydroxylation with 83.2% AUC. Our model is the first of its kind and may be useful for the development of safer drugs. The epoxidation model is available at http://swami.wustl.edu/xenosite.

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 (kcat, Km, and kcat/Km), 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 (kcat and Km) were more consistent with experimental values than those for catalytic efficiency (kcat/Km). 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.

Coumadin (R/S-warfarin) is a commonly prescribed anticoagulant for over ∼20 million Americans. Although highly efficacious, positive clinical outcomes during warfarin therapy depend on maintaining a narrow therapeutic range for the drug. This goal is challenging due to large inter-individual variability in patient response, which has been attributed to diversity in drug metabolism. Warfarin is given as a racemic mixture and evidence suggest differences of R and S-warfarin in their therapeutic activities and metabolism. Previous investigation of warfarin metabolism has been hampered by the inability to quantify the individual enantiomers. To overcome this limitation a multi-mode LC–MS/MS method is reported. This strategy combines phenyl based reverse phase chromatography with chiral phase chromatography prior to quantitation by liquid chromatography tandem mass spectrometry. This approach was made possible through advances in UPLC technology producing narrow peaks suitable for transferring to a second column. The reported method separated individual R and S enantiomers of hydroxywarfarin and warfarin. All four possible isomers of 10-hydroxywarfarin were resolved to reveal unprecedented insights into the stereo-specific metabolism of warfarin. Characterization of the method demonstrated that it is robust and sensitive with inter-day coefficients of error between <7% and a detection limit of 2 nM in sample or 10 fmol on column for each analyte. Individual metabolites may be suitable surrogate biomarkers or predictive markers that predict warfarin dose, adverse interactions, or other important clinical outcomes during anticoagulant therapy. Consequently, the metabolite profiles obtained through this dual phase UPLC–MS/MS method are expected to increase our understanding of the role warfarin metabolism plays in patient response to therapy and yield new strategies to improve patient outcomes.

Coumadin (R/S-warfarin) anticoagulant therapy poses a risk to over 50 million Americans, in part due to interpersonal variation in drug metabolism. Consequently, it is important to understand how metabolic capacity is influenced among patients. Cytochrome P450s (P450 or CYP for a specific isoform) catalyze the first major step in warfarin metabolism to generate five hydroxywarfarins for each drug enantiomer. These primary metabolites are thought to reach at least 5-fold higher levels in plasma than warfarin. We hypothesized that hydroxywarfarins inhibit the hydroxylation of warfarin by CYP2C9, thereby limiting enzymatic capacity toward S-warfarin. To test this hypothesis, we investigated the ability of all five racemic hydroxywarfarins to block CYP2C9 activity toward S-warfarin using recombinant enzyme and human liver microsomes. We initially screened for the inhibition of CYP2C9 by hydroxywarfarins using a P450-Glo assay to determine IC50 values for each hydroxywarfarin. Compared to the substrate, CYP2C9 bound its hydroxywarfarin products with less affinity but retained high affinity for 10- and 4′-hydroxywarfarins, products from CYP3A4 reactions. S-Warfarin steady-state inhibition studies with recombinant CYP2C9 and pooled human liver microsomes confirmed that hydroxywarfarin products from CYP reactions possess the capacity to competitively inhibit CYP2C9 with biologically relevant inhibition constants. Inhibition of CYP2C9 by 7-hydroxywarfarin may be significant given its abundance in human plasma, despite its weak affinity for the enzyme. 10-Hydroxywarfarin, which has been reported as the second most abundant plasma metabolite, was the most potent inhibitor of CYP2C9, displaying approximately 3-fold higher affinity than S-warfarin. These results indicate that hydroxywarfarin metabolites produced by CYP2C9 and other CYPs may limit metabolic capacity toward S-warfarin through competitive inhibition. Subsequent processing of hydroxywarfarins to secondary metabolites, such as hydroxywarfarin glucuronides, could suppress product feedback inhibition, and therefore could play an important role in the modulation of metabolic pathways governing warfarin inactivation and elimination.

As a promiscuous redox partner, the biological role of cytochrome P450 reductase (CPR) depends significantly on protein–protein interactions. We tested a hypothesized CPR docking site by mutating D113, E115, and E116 to alanine and assaying activity toward various electron acceptors as a function of ionic strength. Steady-state cytochrome c studies demonstrated the mutations improved catalytic efficiency and decreased the impact of ionic strength on catalytic parameters when compared to wild type. Based on activity toward 7-ethoxy-4-trifluoro-methylcoumarin, CYP2B1 and CPR favored formation of an active CYP2B1•CPR complex and inactive (CYP2B1)2•CPR complex until higher ionic strength whereby only the binary complex was observed. The mutations increased dissociation constants only for the binary complex and suppressed the ionic strength effect. Studies with a non-binding substrate, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) suggest changes in activity toward cytochrome c and CYP2B1 reflect alterations in the route of electron transfer caused by the mutations. Electrostatic modeling of catalytic and binding parameters confirmed the importance of D113 and especially the double mutant E115 and E116 as mediators in forming charge–charge interactions between CPR and complex partners.

•Metabolically active CYP2E1 localizes to the endoplasmic reticulum and mitochondria.•CYP2E1 localization impacts substrate affinity, stoichiometry, and turnover.•Toxicological relevance of CYP2E1 subcellular localization may depend on substrate.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.Download high-res image (269KB)Download full-size image

•Synthesis of R- and S-7-hydroxywarfarin through green chemistry.•Two novel 7-hydroxywarfarin glucuronides generated by microsomal hepatic UGTs.•Prediction of glucuronide metabolite structures based on spectral properties.•Glucuronidation involved Michaelis–Menten and substrate inhibition kinetics.•Multiple UGTs in HLMs contribute to R- and S-7-hydroxywarfarin metabolism.The widely used anticoagulant Coumadin (R/S-warfarin) undergoes oxidation by cytochromes P450 into hydroxywarfarins that subsequently become conjugated for excretion in urine. Hydroxywarfarins may modulate warfarin metabolism transcriptionally or through direct inhibition of cytochromes P450 and thus, UGT action toward hydroxywarfarin elimination may impact levels of the parent drugs and patient responses. Nevertheless, relatively little is known about conjugation by UDP-glucuronosyltransferases in warfarin metabolism. Herein, we identified probable conjugation sites, kinetic mechanisms and hepatic UGT isoforms involved in microsomal glucuronidation of R- and S-7-hydroxywarfarin. Both compounds underwent glucuronidation at C4 and C7 hydroxyl groups based on elution properties and spectral characteristics. Their formation demonstrated regio- and enantioselectivity by UGTs and resulted in either Michaelis–Menten or substrate inhibition kinetics. Glucuronidation at the C7 hydroxyl group occurred more readily than at the C4 group, and the reaction was overall more efficient for R-7-hydroxywarfarin due to higher affinity and rates of turnover. The use of these mechanisms and parameters to model in vivo clearance demonstrated that contributions of substrate inhibition would lead to underestimation of metabolic clearance than that predicted by Michaelis–Menten kinetics. Lastly, these processes were driven by multiple UGTs indicating redundancy in glucuronidation pathways and ultimately metabolic clearance of R- and S-7-hydroxywarfarin.

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 = 32 mM). 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.Download full-size image

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 individual’s 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.