Predictive models for cytochrome p450 isozymes based on quantitative high throughput screening data

The human cytochrome P450 (CYP450) isozymes are the most important enzymes in the body to metabolize many endogenous and exogenous substances including environmental toxins and therapeutic drugs. Any unnecessary interactions between a small molecule and CYP450 isozymes may raise a potential to disarm the integrity of the protection. Accurately predicting the potential interactions between a small molecule and CYP450 isozymes is highly desirable for assessing the metabolic stability and toxicity of the molecule. The National Institutes of Health Chemical Genomics Center (NCGC) has screened a collection of over 17,000 compounds against the five major isozymes of CYP450 (1A2, 2C9, 2C19, 2D6, and 3A4) in a quantitative high throughput screening (qHTS) format. In this study, we developed support vector classification (SVC) models for these five isozymes using a set of customized generic atom types. The CYP450 data sets were randomly split into equal-sized training and test sets. The optimized SVC models exhibited high predictive power against the test sets for all five CYP450 isozymes with accuracies of 0.93, 0.89, 0.89, 0.85, and 0.87 for 1A2, 2C9, 2C19, 2D6, and 3A4, respectively, as measured by the area under the receiver operating characteristic (ROC) curves. The important atom types and features extracted from the five models are consistent with the structural preferences for different CYP450 substrates reported in the literature. We also identified novel features with significant discerning power to separate CYP450 actives from inactives. These models can be useful in prioritizing compounds in a drug discovery pipeline or recognizing the toxic potential of environmental chemicals.     

In silico modeling of p450 substrates, inhibitors, activators, and inducers

Cytochrome P450 enzymes are the predominant mediators of phase I metabolism of exogenous small molecules. As a result of their extensive role in metabolism of xenobiotics, drug compounds, and endogenous compounds, as well as their wide tissue distribution, significant drug discovery resources are spent to avoid interacting with this class of enzymes. Here we review historical and recent in silico modeling of 7 cytochrome P450 enzymes of particular interest, specifically CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP3A4. For each we provide a brief biological background including known inhibitors, substrates, and inducers, as well as details of computational modeling efforts and advances in structural biology. We also provide similar details for 3 nuclear receptors known to regulate gene expression of these enzyme families.     

Comparison of metabolic soft spot predictions of CYP3A4, CYP2C9 and CYP2D6 substrates using MetaSite and StarDrop

Metabolite identification study plays an important role in determining the sites of metabolic liability of new chemical entities (NCEs) in drug discovery for lead optimization. Here we compare the two predictive software, MetaSite and StarDrop, available for this purpose. They work very differently but are used to predict the site of oxidation by major human cytochrome P450 (CYP) isoforms. Neither software can predict non-CYP catalyzed metabolism nor the rates of metabolism. For the purpose of comparing the two software packages, we tested known probe substrate for these enzymes, which included 12 substrates of CYP3A4 and 18 substrates of CYP2C9 and CYP2D6 were analyzed by each software and the results were compared. It is possible that these known substrates were part of the training set but we are not aware of it. To assess the performance of each software we assigned a point system for each correct prediction. The total points assigned for each CYP isoform experimentally were compared as a percentage of the total points assigned theoretically for the first choice prediction for all substrates for each isoform. Our results show that MetaSite and StarDrop are similar in predicting the correct site of metabolism by CYP3A4 (78% vs 83%, respectively). StarDrop appears to do slightly better in predicting the correct site of metabolism by CYP2C9 and CYP2D6 metabolism (89% and 93%, respectively) compared to MetaSite (63% and 70%, respectively). The sites of metabolism (SOM) from 34 in-house NCEs incubated in human liver microsomes or human hepatocytes were also evaluated using two prediction software packages and the results showed comparable SOM predictions. What makes this comparison challenging is that the contribution of each isoform to the intrinsic clearance (Clint) is not known. Overall the software were comparable except for MetaSite performing better for CYP2D6 and that MetaSite has a liver model that is absent in StarDrop that predicted with 82% accuracy.     

Liquid chromatography/tandem mass spectrometric determination of inhibition of human cytochrome P450 isozymes by resveratrol and resveratrol-3-sulfate

trans-Resveratrol, a phenolic phytoalexin occurring in grapes, wine, peanuts, and cranberries, has been reported to both have anticarcinogenic, antioxidative, phytoestrogenic, and cardioprotective activities, and to be a weak inhibitor of cytochrome P450 (CYP)3A4, which might have significance for drug-drug interactions. Since trans-resveratrol is rapidly converted in vivo to primarily trans-resveratrol-3-sulfate, a rapid, selective, and sensitive method using liquid chromatography/tandem mass spectrometry (LC/MS/MS) was developed to investigate human cytochrome P450 inhibition by trans-resveratrol-3-sulfate. Effects of trans-resveratrol and trans-resveratrol-3-sulfate on the metabolism of selective cytochrome P450 substrates (CYP1A2/ethoxyresorufin, CYP2C9/diclofenac, CYP2C19/(S)-mephenytoin, CYP2D6/bufuralol, CYP3A4/testosterone) were monitored using cDNA-expressed human recombinant isozymes. For method validation, LC/MS/MS was used to measure the inhibition of various cytochrome P450 isozymes by different concentrations (0-50 microM) of known selective inhibitors. IC(50) values of 3.2, 1.4, 8.9, 0.2, and 0.3 microM were obtained for the standard isozyme inhibitors CYP1A2/furafylline, CYP2C9/sulfaphenazole, CYP2C19/tranylcypromine, CYP2D6/quinidine, and CYP3A4/ketoconazole, respectively, which were in good agreement with literature values. trans-Resveratrol showed IC(50) values of 11.6 microM for CYP2C19 and 1.1 microM for CYP3A4, but the IC(50) values exceeded 50 microM for all the other CYP isozymes, which indicated no inhibition. No enzyme inhibition was observed for trans-resveratrol-3-sulfate. Our results indicate that trans-resveratrol is a marginal inhibitor of CYP3A4 and a weak inhibitor of CYP2C19, but its major metabolite trans-resveratrol-3-sulfate is not an inhibitor of any of the cytochrome P450 isozymes investigated.     

Cytochrome P450 reaction phenotyping and inhibition and induction studies of pinostrobin in human liver microsomes and hepatocytes

Pinostrobin (PI, 5‐hydroxy‐7‐methoxyflavanone) is a natural flavonoid known for its rich pharmacological activities. The objective of this study was to identify the human liver cytochrome P450 (CYP450) isoenzymes involved in the metabolism of PI. A single hydoxylated metabolite was obtained from PI after an incubation with pooled human liver microsomes (HLMs). The relative contributions of different CYP450s were evaluated using CYP450‐selective inhibitors in HLMs and recombinant human CYP450 enzymes, and the results revealed the major involvement of CYP1A2, CYP2C9 and CYP2E1 in PI metabolism. We also evaluated the ability of PI to inhibit and induce human cytochrome P450 enzymes in vitro . High‐performance liquid chromatography and liquid chromatography–tandem mass spectrometry analytical techniques were used to estimate the enzymatic activities of seven drug‐metabolizing CYP450 isozymes in vitro . In HLMs, PI did not inhibit CYP 1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6 or CYP3A4 (IC₅₀ > 100 μm ). In the induction studies, PI had minimal effects on CYP1A2, CYP2B6and CYP3A4 activity. Based on these results, PI would not be expected to cause clinically significant CYP450 inhibition or induction.     

Preparation, characterization, and substrate metabolism of gold-immobilized cytochrome P450 2C9

The cytochrome P450 enzymes represent an important class of heme-containing enzymes. There is considerable interest in immobilizing these enzymes on a surface so that interactions between a single enzyme and other species can be studied with respect to electron transfer, homodimer or heterodimer interactions, or for construction of biological-based chips for standardizing cytochrome P450 metabolism or for high-throughput screening of pharmaceutical agents. Previous studies have generally immobilized P450 enzymes in a matrix or on a surface. Here, we have attached CYP2C9 to gold substrates such that the resulting construct maintains the ability to bind and metabolize substrates in the presence of NADPH and cytochrome P450 reductase. The activity of these chips is directly dependent upon the linkers used to attach CYP2C9 and to the presence of key molecules in the active site during enzyme attachment. A novel method to detect substrate-enzyme binding, namely, superconducting quantum interference device (SQUID) magnetometry, was used to monitor the binding of substrates. Most significantly, conditions that allow measurable CYP2C9 metabolism to occur have been developed.     

2D SMARTCyp reactivity-based site of metabolism prediction for major drug-metabolizing cytochrome P450 enzymes

Cytochrome P450 (CYP) 3A4, 2D6, 2C9, 2C19, and 1A2 are the most important drug-metabolizing enzymes in the human liver. Knowledge of which parts of a drug molecule are subject to metabolic reactions catalyzed by these enzymes is crucial for rational drug design to mitigate ADME/toxicity issues. SMARTCyp, a recently developed 2D ligand structure-based method, is able to predict site-specific metabolic reactivity of CYP3A4 and CYP2D6 substrates with an accuracy that rivals the best and more computationally demanding 3D structure-based methods. In this article, the SMARTCyp approach was extended to predict the metabolic hotspots for CYP2C9, CYP2C19, and CYP1A2 substrates. This was accomplished by taking into account the impact of a key substrate-receptor recognition feature of each enzyme as a correction term to the SMARTCyp reactivity. The corrected reactivity was then used to rank order the likely sites of CYP-mediated metabolic reactions. For 60 CYP1A2 substrates, the observed major sites of CYP1A2 catalyzed metabolic reactions were among the top-ranked 1, 2, and 3 positions in 67%, 80%, and 83% of the cases, respectively. The results were similar to those obtained by MetaSite and the reactivity + docking approach. For 70 CYP2C9 substrates, the observed sites of CYP2C9 metabolism were among the top-ranked 1, 2, and 3 positions in 66%, 86%, and 87% of the cases, respectively. These results were better than the corresponding results of StarDrop version 5.0, which were 61%, 73%, and 77%, respectively. For 36 compounds metabolized by CYP2C19, the observed sites of metabolism were found to be among the top-ranked 1, 2, and 3 sites in 78%, 89%, and 94% of the cases, respectively. The computational procedure was implemented as an extension to the program SMARTCyp 2.0. With the extension, the program can now predict the site of metabolism for all five major drug-metabolizing enzymes with an accuracy similar to or better than that achieved by the best 3D structure-based methods. Both the Java source code and the binary executable of the program are freely available to interested users.     

Automated information extraction and structure-activity relationship analysis of cytochrome P450 substrates

Information on CYP-chemical interactions was comprehensively explored by a text-mining technique, to confirm our previous structure-activity relationship model for CYP substrates (Yamashita et al. J. Chem. Inf. Model. 2008, 48, 364-369). The text-mining technique is based on natural language processing and can extract chemical names and their interaction patterns according to sentence context. After chemicals were automatically extracted and classified into CYP substrates, inhibitors, and inducers, 709 substrates were retrieved from the PubChem database and categorized as 216, 145, 136, 217, 156, and 379 substrates for CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, respectively. Although the previous classification model was developed using data from only 161 compounds, the model classified the substrates found by text-mining analysis with reasonable accuracy. This confirmed the validity of both the multi-objective classification model for CYP substrates and the text-mining procedure.     

Convenient Gram-Scale Metabolite Synthesis by Engineered Fission Yeast Strains Expressing Functional Human P450 Systems

The growing need for the characterization of cytochrome P450 (P450) metabolites often necessitates their synthesis up to Gram-scale. This task may in principle be achieved by using various techniques including chemical synthesis, the use of laboratory animals, in vitro P450 systems or microbial biotransformation. However, these approaches are in many instances unfavorable due to low yields, laborious purification, costs of cofactors, or the formation of non-physiologic metabolites. The fission yeast Schizosaccharomyces pombe has previously been shown by others and us to be very well suited for the heterologous expression of human P450s. In this study, we demonstrate whole-cell biotransformation reactions carried out with fission yeast strains that coexpress human cytochrome P450 reductase (CPR) and one of the following P450 isoforms: CYP2B6, CYP2C9, CYP2C19, CYP2D6, or CYP3A4, respectively. These strains could successfully convert their respective standard substrates but showed different responses with respect to incubation pH, the presence of glucose, and temperature, respectively. In addition, the preparative of synthesis of 2.8?g of 4'-hydroxydiclofenac was achieved by whole-cell biotransformation of diclofenac using a CPR-CYP2C9 coexpressing fission yeast strain.     

流动注射-串联质谱法快速筛选植物提取物对药物代谢酶的抑制潜性

基于流动注射-串联质谱(FI-MS/MS)技术,综合运用微孔板低温孵育、低温样品相分离、六合一进样和多探针同时定量分析,建立了体外快速筛选复杂中草药组分抑制细胞色素P450同工酶活性的新方法。FI-MS/MS新技术不依赖色谱分离,也无需复杂的仪器设备和样品前处理,分析通量大幅提高。采用基质匹配标样和内标法定量,方法的准确度、精密度、灵敏度和线性范围均满足孵育体系中代谢物分析的需要。优化条件下测得各酶动力学参数和对照抑制剂的半抑制率(IC50)均与参考值一致;仅用两抑制剂浓度点即可获得满意的IC50近似值,两点IC50的测试重现性良好。结合两点IC50程序,FI-MS/MS方法成功应用于中药粗提液对人体主要P450酶抑制效应的快速筛选与分类,为多组分体系复杂药代相互作用的研究和评估提供了一种有效手段。     

CYPstrate: A Set of Machine Learning Models for the Accurate Classification of Cytochrome P450 Enzyme Substrates and Non-Substrates

The interaction of small organic molecules such as drugs, agrochemicals, and cosmetics with cytochrome P450 enzymes (CYPs) can lead to substantial changes in the bioavailability of active substances and hence consequences with respect to pharmacological efficacy and toxicity. Therefore, efficient means of predicting the interactions of small organic molecules with CYPs are of high importance to a host of different industries. In this work, we present a new set of machine learning models for the classification of xenobiotics into substrates and non-substrates of nine human CYP isozymes: CYPs 1A2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D6, 2E1, and 3A4. The models are trained on an extended, high-quality collection of known substrates and non-substrates and have been subjected to thorough validation. Our results show that the models yield competitive performance and are favorable for the detection of CYP substrates. In particular, a new consensus model reached high performance, with Matthews correlation coefficients (MCCs) between 0.45 (CYP2C8) and 0.85 (CYP3A4), although at the cost of coverage. The best models presented in this work are accessible free of charge via the “CYPstrate” module of the New E-Resource for Drug Discovery (NERDD).     

Simultaneous determination of bupropion,metroprolol, midazolam,phenacetin, omeprazole and tolbutamide in rat plasma by UPLC‐MS/MS and its application to cytochrome P450 activity study in rats

A specific ultra‐performance liquid chromatography tandem mass spectrometry method is described for the simultaneous determination of bupropion, metroprolol, midazolam, phenacetin, omeprazole and tolbutamide in rat plasma with diazepam as internal standard, which are the six probe drugs of the six cytochrome P450 isoforms CYP2B6, CYP2D6, CYP3A4, CYP1A2, CYP2C19 and CYP2C9. Plasma samples were protein precipitated with acetonitrile. The chromatographic separation was achieved using a UPLC® BEH C₁₈ column (2.1 × 100 mm, 1.7 µm). The mobile phase consisted of acetonitrile and water (containing 0.1% formic acid) with gradient elution. The triple quadrupole mass spectrometric detection was operated by multiple reaction monitoring in positive electrospray ionization. The precisions were <13%, and the accuracy ranged from 93.3 to 110.4%. The extraction efficiency was >90.5%, and the matrix effects ranged from 84.3 to 114.2%. The calibration curves in plasma were linear in the range of 2–2000 ng/mL, with correlation coefficient (r²) >0.995. The method was successfully applied to pharmacokinetic studies of the six probe drugs of the six CYP450 isoforms and used to evaluate the effects of erlotinib on the activities of CYP2B6, CYP2D6, CYP3A4, CYP1A2, CYP2C19 and CYP2C9 in rats. Erlotinib may inhibit the activity of CYP2B6 and CYP3A4, and may induce CYP2C9 of rats. Copyright © 2015 John Wiley & Sons, Ltd.     

Site-Directed Mutagenesis of Cytochrome P450 2D6 and 2C19 Enzymes: Expression and Spectral Characterization of Naturally Occurring Allelic Variants

Genetic polymorphism of the cytochrome P450 (CYP) genes particularly affects CYP2D6 and CYP2C19 to a functionally relevant extent, and it is therefore crucial to elucidate the enzyme kinetic and molecular basis for altered catalytic activity of these allelic variants. This study explored the expression and function of the reported alleles CYP2D6*2, CYP2D6*10, CYP2D6*17, CYP2C19*23, CYP2C19*24, and CYP2C19*25 with respect to gene polymorphisms. Site-directed mutagenesis (SDM) was carried out to generate these six alleles. After DNA sequencing, the CYP2D6 and CYP2C19 wild types alongside with their alleles were each independently co-expressed with NADPH-CYP oxidoreductase (OxR) in Escherichia coli. The expressed proteins were analyzed using Western blotting, reduced carbon monoxide (CO) difference spectral scanning, and cytochrome c reductase assay. Results from Western blot revealed the presence of all CYP wild-type and allelic proteins in E. coli membrane fractions. The reduced CO difference spectra scanning presented the distinct peak of absorbance at 450 nm, and the cytochrome c reductase assay has confirmed that spectrally active OxR was expressed in each protein preparation. As a conclusion, the results obtained from this study have proven the CYP variants to be immunoreactive and spectrally active and are suitable for use to examine biotransformation and interaction mechanism of the enzymes.     

High-throughput screening of inhibitory potential of nine cytochrome P450 enzymes in vitro using liquid chromatography/tandem mass spectrometry

The early detection of potential drug-drug interactions is an important issue of drug discovery that has led to the development of high-throughput screening (HTS) methods for potential drug interactions. We developed a HTS method for potential interactions of inhibitory drugs for nine human P450 enzymes using cocktail incubation and tandem mass spectrometry in vitro. This new method involves incubation of two cocktail doses and single cassette analysis. The two cocktail doses in vitro were developed to minimize solvent effects and mutual drug interactions among substrates: cocktail A was composed of phenacetin for CYP1A2, coumarin for CYP2A6, paclitaxel for CYP2C8, S-mephenytoin for CYP2C19, dextromethorphan for CYP2D6, and midazolam for CYP3A4; and cocktail B was composed of three substrates including bupropion for CYP2B6, tolbutamide for CYP2C9, and chlorzoxazone for CYP2E1. In the incubation study of these cocktails, the reaction mixtures were pooled and simultaneously analyzed using liquid chromatography/tandem mass spectrometry employing a fast gradient. The method was validated by comparing the inhibition data obtained from the incubation of each individual probe substrate alone with data from the new method. The IC50 value of each inhibitor in the cocktail agreed well with that of the individual probe drug as well as with values previously reported in the literature. As a HTS method for potential interactions of the inhibition of these nine P450 enzymes, this new method will be useful in the drug discovery process and for the mechanistic understanding of drug interactions.     

In vitro interaction cocktail assay for nine major cytochrome P450 enzymes with 13 probe reactions and a single LC/MSMS run: analytical validation and testing with monoclonal anti-CYP antibodies

A sensitive and rugged LC/MSMS method was developed for a comprehensive in vitro metabolic interaction screening assay with N-in-1 approach reported earlier. A cocktail consisting of ten cytochrome P450 (CYP)-selective probe substrates with known kinetic, metabolic and interaction properties in vivo was incubated in a pool of human liver microsomes, and metabolites of melatonin (CYP1A2), coumarin (CYP2A6), bupropion (CYP2B6), amodiaquine (CYP2C8) tolbutamide (CYP2C9), omeprazole (CYP2C19 and CYP3A4), dextromethorphan (CYP2D6), chlorzoxazone (CYP2E1), midazolam (CYP3A4) and testosterone (CYP3A4) were simultaneously analysed with a single LC/MSMS run. Altogether, 13 metabolites and internal standard phenacetin were analysed in multiple reaction mode. Polarity switching mode was utilized to acquire negative ion mode electrospray data for hydroxychlorzoxazone and positive ionization data for the rest of the analytes. Fast gradient elution was applied, giving total injection cycle of 8 min. The method was modified for two different LC/MSMS systems, and was validated for linear range, detection limit, accuracy and precision for each metabolite. In addition, cocktail inhibition system was further tested using monoclonal anti-CYP antibodies as inhibitors for each probe reaction.