Phenotyping of CYP450 in human liver microsomes using the cocktail approach

The cocktail approach is an advantageous strategy used to monitor the activities of several cytochromes P450 (CYPs) in a single test to increase the throughput of in vitro phenotyping studies. In this study, a cocktail mixture was developed with eight CYP-specific probe substrates to simultaneously evaluate the activity of the most important CYPs, namely, CYP1A2, CYP2A6, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and the CYP3A subfamily. After cocktail incubation in the presence of human liver microsomes (HLMs), the eight selected substrates and their specific metabolites were analyzed by ultra-high-pressure liquid chromatography and electrospray ionization quadrupole time-of-flight mass spectrometry. Qualitative and quantitative data were simultaneously acquired to produce an overview of the extended phase I biotransformation routes for each probe substrate in the HLMs and to generate phenotypic profiles of various HLMs. A comparison of the cocktail strategy with an individual substrate assay for each CYP produced similar results. Moreover, the cocktail was tested on HLMs with different allelic variants and/or in the presence of selective inhibitors. The results were in agreement with the genetic polymorphisms of the CYPs and the expected effect of the alterations. All of these experiments confirmed the reliability of this cocktail assay for phenotyping of the microsomal CYPs.     

Optimization and validation of a label-free MRM method for the quantification of cytochrome P450 isoforms in biological samples

Cytochromes P450 (CYPs) play critical roles in oxidative metabolism of many endogenous and exogenous compounds. Protein expression levels of CYPs in liver provide relevant information for a better understanding of the importance of CYPs in pharmacology and toxicology. This work aimed at establishing a simple method to quantify six CYPs (CYP3A4, CYP3A5, CYP1A2, CYP2D6, CYP2C9, and CYP2J2) in various biological samples without isotopic labeling. The biological matrix was spiked with the standard peptides prior to the digestion step to realize a label-free quantification by mass spectrometry. The method was validated and applied to quantify these six isoforms in both human liver microsomes and mitochondria, but also in recombinant expression systems such as baculosomes and the HepG2 cell line. The results showed intra-assay and interassay accuracy and precision within 16 % and 5 %, respectively, at the low quality control level, and demonstrated the advantages of the method in terms of reproducibility and cost.

Flexibility of human cytochromes P450: molecular dynamics reveals differences between CYPs 3A4, 2C9, and 2A6, which correlate with their substrate preferences

Molecular dynamics (MD) simulations at normal and high temperature were used to study the flexibility and malleability of three microsomal cytochromes P450 (CYPs): CYP3A4, CYP2C9, and CYP2A6. Comparison of B-factors (describing the atomic fluctuations) between X-ray and MD data shows that the X-ray B-factors are significantly lower in the regions where the crystal contacts occur than for other regions. Consequently, the conclusions about CYP flexibility based solely on the X-ray data might be misleading. Comparison of flexibility patterns of the three CYPs enabled common features and variations in flexibility and malleability of the studied CYPs to be identified. The previously described pattern of flexibility in topological elements of microsomal CYPs (a rigid heme binding core, a malleable distal side and intermediately flexible proximal side) was confirmed. These topological features provide an important combination of high stereo- and regio-specificity (mediated by the relative rigidity in the neighborhood of the heme), together with high substrate promiscuity due to the more flexible active site and the malleability of the distal side. The data acquired here show that the malleability of the three studied CYPs correlates with their substrate specificity: CYP2A6 has a narrow substrate range and is the most rigid, CYP3A4 is the most promiscuous CYP known and is the most malleable, and CYP2C9 is intermediate in terms of both its substrate specificity and malleability. Thus, the malleability of CYPs is probably a major determinant of their substrate specificity.     

Prediction of cytochrome P450 3A4, 2D6, and 2C9 inhibitors and substrates by using support vector machines

Statistical learning methods have been used in developing filters for predicting inhibitors of two P450 isoenzymes, CYP3A4 and CYP2D6. This work explores the use of different statistical learning methods for predicting inhibitors of these enzymes and an additional P450 enzyme, CYP2C9, and the substrates of the three P450 isoenzymes. Two consensus support vector machine (CSVM) methods, "positive majority" (PM-CSVM) and "positive probability" (PP-CSVM), were used in this work. These methods were first tested for the prediction of inhibitors of CYP3A4 and CYP2D6 by using a significantly higher number of inhibitors and noninhibitors than that used in earlier studies. They were then applied to the prediction of inhibitors of CYP2C9 and substrates of the three enzymes. Both methods predict inhibitors of CYP3A4 and CYP2D6 at a similar level of accuracy as those of earlier studies. For classification of inhibitors of CYP2C9, the best CSVM method gives an accuracy of 88.9% for inhibitors and 96.3% for noninhibitors. The accuracies for classification of substrates and nonsubstrates of CYP3A4, CYP2D6, and CYP2C9 are 98.2 and 90.9%, 96.6 and 94.4%, and 85.7 and 98.8%, respectively. Both CSVM methods are potentially useful as filters for predicting inhibitors and substrates of P450 isoenzymes. These methods generally give better accuracies than single SVM classification systems, and the performance of the PP-CSVM method is slightly better than that of the PM-CSVM method.     

Coexpression of CPR from Various Origins Enhances Biotransformation Activity of Human CYPs in S. pombe

Cytochrome P450 enzymes (CYPs or P450s) are the most important enzymes involved in the phase I metabolism of drugs (and other xenobiotics) in humans, and the corresponding drug metabolites are needed as reference substances for their structural confirmation and for pharmacological or toxicological characterization. We have previously shown that biotechnological synthesis of such metabolites is feasible by whole-cell biotransformation with human CYPs recombinantly expressed in the fission yeast Schizosaccharomyces pombe. It was the aim of this study to compare the activity of seven human microsomal CYPs (CYP2C9, CYP2D6, CYP3A4, CYP3A5, CYP3A7, CYP17, and CYP21) upon coexpression with NADPH-cytochrome P450 oxidoreductases (CPRs) from various origins, namely, human CPR (hCPR) and its homologues from fission yeast (ccr1) and the bishop’s weed Ammi majus (AmCPR), respectively. For this purpose, 28 recombinant strains were needed, with five of them having been constructed previously and 23 strains being newly constructed. Bioconversion experiments showed that coexpression of a CPR does not only influence the reaction rate but, in some cases, also exerts an influence on the metabolite pattern. For CYP3A enzymes, coexpression of hCPR yielded the best results, while for another two, hCPR was equally helpful as ccr1 (both CYP17 and CYP21) or AmCPR (CYP17 only), respectively. Interestingly, CYP2D6 displayed its highest activity when coexpressed with ccr1 and CYP2C9 with AmCPR. These results corroborate the view of CPR as a well-suited bio-brick in synthetic biology for the construction of artificial enzyme complexes.     

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B′ and F–G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.     

Direct inhibition of Re Du Ning Injection and its active compounds on human liver cytochrome P450 enzymes by a cocktail method

The aim of this study was to investigate the direct inhibitory effects of Re Du Ning Injection (RDN) and its active compounds on the major cytochrome P450 enzyme (CYP) isoforms (CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) of human liver microsomes by ‘a cocktail method’. The activity of each CYP isform was represented as the formation rate of the specific metabolite from relevant substrate. Then a sensitive and specific ultra‐performance liquid chromatography–tandem mass spectrometry (UPLC–MS/MS) method was developed and validated to simultaneously analyze the seven metabolites. RDN (0.035–2.26 mg/mL) showed a strong inhibitiory effect on CYP2C8, followed by CYP2C9, CYP2B6, CYP2C19, CYP1A2 and CYP3A4. The IC₅₀ value for each enzyme was 0.19, 0.66, 0.72, 1.27, 1.66 and 2.13 mg/mL, respectively. RDN competitively inhibited the activities of CYP1A2 (K _i = 1.22 mg/mL), CYP2B6 (K _i = 0.65 mg/mL) and CYP3A4 (K _i = 0.88 mg/mL); it also exhibited mixed inhibition of CYP2C8, CYP2C9 and CYP2C19 with a K _i value of 0.26, 0.64 and 0.82 mg/mL, respectively. However, the activity of CYP2D6 was not significantly inhibited even by 2.26 mg/mL RDN. Moreover, the data of nine active compounds on the CYPs showed that cryptochlorogenin acid, sochlorogenic acid B and sochlorogenic acid C were the major contributors to the inhibitory effect of RDN on CYP2C8, while the inhibitory effect of RDN on CYP2C9 might be caused by sochlorogenic acid A and sochlorogenic acid C. Moreover, neochlorogenic acid might be the major contributor to the inhibitory effect on CYP2B6. All of the findings suggested that drug–drug interactions may occur and great caution should be taken when RDN is combined with drugs metabolized by these CYPs.     

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.     

Development and validation of a liquid-chromatography high-resolution tandem mass spectrometry approach for quantification of nine cytochrome P450 (CYP) model substrate metabolites in an in vitro CYP inhibition cocktail

Knowledge about the cytochrome P450 (CYP) inhibition potential of new drug candidates is important for drug development because of its risk of interactions. For novel psychoactive substances (NPS), corresponding data are not available. For developing a general drug inhibition cocktail assay, a liquid-chromatography high-resolution tandem mass spectrometry multi-analyte approach was developed and validated for quantifying low concentrations of O-diethyl phenacetin for CYP 1A2, 7-hydroxy coumarin for CYP 2A6, 4-hydroxy bupropion for CYP 2B6, N-diethyl amodiaquine for CYP 2C8, 4-hydroxy diclofenac for CYP 2C9, 5-hydroxy omeprazole for CYP 2C19, O-dimethyl dextromethorphan for CYP 2D6, 6-hydroxy chlorzoxazone for CYP 2E1, and 6-beta-hydroxy testosterone for CYP 3A in the incubation mixture in the presence of substrates and inhibitors. The tested matrix effects ranged from 63 to 141 % and the recoveries from 95 to 110 %. Time-saving one-point calibration allowed sufficient quantification, although some of the validation results for 7-hydroxy coumarin, 4-hydroxy bupropion, 4-hydroxy diclofenac, and 6-beta-hydroxy testosterone were outside the acceptance criteria (AC) but without influence of the IC₅₀ calculation. Validation showed also that the approach was sensitive and selective using mass spectral multiplexing. In conclusion, the presented assay was suitable for the quantification of the model substrate metabolites and could be used for the development of a CYP inhibition assay for testing most CYPs and a wide range of drugs of abuse.     

Capillary electrophoresis to assess drug metabolism induced in vitro using single CYP450 enzymes (Supersomes): application to the chiral metabolism of mephenytoin and methadone

Capillary electrophoresis (CE) with multiwavelength absorbance detection is demonstrated to be an effective tool for the assessment of in vitro drug metabolism studies using microsomes containing single human cytochrome P450 enzymes (CYPs) expressed in baculovirus-infected insect cells (Supersomes). Mephenytoin (MEPH), dextromethorphan, diclofenac, caffeine, and methadone (MET) were successfully applied as test substrates for CYP2C19, CYP2D6*1, CYP2C9*1, CYP1A2, and CYP3A4, respectively. For each system, the CE-based assay could be shown to permit the simultaneous analysis of the parent drug and its targeted metabolite. Using a chiral micellar electrokinetic capillary chromatography assay, the aromatic hydroxylation of MEPH catalyzed by CYP2C19 could thereby be confirmed to be highly stereoselective, an aspect that is in agreement with data obtained via urinary analysis after intake of racemic MEPH by extensive metabolizer phenotypes. The MET to 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) conversion was investigated with a chiral zone electrophoresis assay. Incubation of racemic and nonracemic MET with CYP3A4 revealed no stereoselectivity for the transformation to EDDP, whereas no EDDP formation was observed with CYP1A2. CYP2C9 and CYP2C19 provided enhanced formation of R-EDDP and CYP2D6 incubation resulted in the preferential conversion to S-EDDP. Investigations using racemic MET and human liver microsomes revealed a modest stereoselectivity with an R/S EDDP ratio < 1 which is similar to the in vivo findings in urine.     

Identification of cytochrome P450 isoforms involved in the metabolism of artocarpin and assessment of its drug–drug interaction

Artocarpin isolated from an agricultural plant Artocarpus communis has shows anti‐inflammation and anticancer activities. In this study, we utilized recombinant human UDP‐glucuronosyltransferasesupersomes (UGTs) and human liver microsomes to explore its inhibitory effect on UGTs and cytochrome p450 enzymes (CYPs). Chemical inhibition studies and screening assays with recombinant human CYPs were used to identify if CYP isoform is involved in artocarpin metabolism. Artocarpin showed strong inhibition against UGT1A3, UGT1A6, UGT1A7, UGT1A8, UGT1A9, UGT1A10, UGT2B7, CYP2C8 and CYP3A4. In particular, artocarpin exhibited competitive inhibition against CYP3A4 and noncompetitive inhibition against UGT1A3 and UGT1A7. The half inhibition concentration values for CYP3A4, UGT1A3 and UGT1A7 were 4.67, 3.82 and 4.82 μm , and the inhibition kinetic parameters for them were 0.78, 2.67 and 3.14 μm , respectively. After artocarpin was incubated in human liver microsomes and determined by HPLC, we observed its main metabolites (M1 and M2). In addition, we proved that CYP2D6 played the key role in the biotransformation of artocarpin in human liver microsomes. The result of molecular docking further confirmed that artocarpin interacted with CYP2D6, CYP2C8 and CYP3A4 through hydrogen bonds. This study provided preliminary results for further research on artocarpin or artocarpin‐containing herbs.     

Insights on cytochrome p450 enzymes and inhibitors obtained through QSAR studies

The cytochrome P450 (CYP) superfamily of heme enzymes play an important role in the metabolism of a large number of endogenous and exogenous compounds, including most of the drugs currently on the market. Inhibitors of CYP enzymes have important roles in the treatment of several disease conditions such as numerous cancers and fungal infections in addition to their critical role in drug-drug interactions. Structure activity relationships (SAR), and three-dimensional quantitative structure activity relationships (3D-QSAR) represent important tools in understanding the interactions of the inhibitors with the active sites of the CYP enzymes. A comprehensive account of the QSAR studies on the major human CYPs 1A1, 1A2, 1B1, 2A6, 2B6, 2C9, 2C19, 2D6, 2E1, 3A4 and a few other CYPs are detailed in this review which will provide us with an insight into the individual/common characteristics of the active sites of these enzymes and the enzyme-inhibitor interactions.     

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.     

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.     

QSAR models for P450 (2D6) substrate activity

Human Cytochrome P450 (CYP) is a large group of enzymes that possess an essential function in metabolising different exogenous and endogenous compounds. Humans have more than 50 different genes encoding CYP enzymes, among these a gene encoding for the CYP isoenzyme 2D6, a CYP able to metabolise drugs and other chemicals. A training set of 747 chemicals primarily based on in vivo human data for the CYP isoenzyme 2D6 was collected from the literature. QSAR models focusing on substrate/non-substrate activity were constructed by the use of MultiCASE, Leadscope and MDL quantitative structure–activity relationship (QSAR) modelling systems. They cross validated (leave-groups-out) with concordances of 71%, 81% and 82%, respectively. Discrete organic European Inventory of Existing Commercial Chemical Substances (EINECS) chemicals were screened to predict an approximate percentage of CYP 2D6 substrates. These chemicals are potentially present in the environment. The biological importance of the CYP 2D6 and the use of the software mentioned above were discussed.