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).     

Structural and functional insights into CYP2C8.3: A genetic polymorph of cytochrome P450 2C8

The cytochrome P450 (CYP) superfamily plays a key role in the oxidative metabolism of a wide range of exogenous chemicals. CYP2C8 is the principal enzyme responsible for the metabolism of the anti-cancer drug paclitaxel in the human liver, and carries out the oxidative metabolism of at least 5% of clinical drugs. Polymorphisms in CYP2C8 have been closely implicated in individualized medication. CYP2C8.3, a common polymorph of CYP2C8 with dual amino acid substitutions R139K and K399R, is found primarily in Caucasians. In this study, CYP2C8.3 and its wild type (WT) CYP2C8 were expressed in E. coli, and their purified proteins were characterized by UV-visible spectroscopy, mass spectrometry, and circular dichroism. Their thermal stability, substrate binding ability, and metabolic activity against paclitaxel were investigated. The electron transfer kinetics during paclitaxel metabolism by WT CYP2C8 or CYP2C8.3 was studied by stopped-flow kinetics. The results revealed that mutations in CYP2C8.3 did not greatly influence the heme active site or protein thermal stability and paclitaxel binding ability, but the metabolic activity against paclitaxel was significantly depressed to just 11% of that of WT CYP2C8. Electron transfer from CYP reductase to CYP2C8.3 was found to be significantly slower than that to WT CYP2C8 during catalysis, and this might be the main reason for the depressed metabolic activity. Since the polymorph CYP2C8.3 is defective in catalyzing substrates of CYP2C8 in vitro, it might be expected to have important clinical and pathophysiological consequences in homozygous individuals, and this study provides valuable information in this aspect.     

一种细胞色素P4502C8基因多态-CYP2C8.3的结构和功能

细胞色素P450超级家族在代谢众多的外源性化学物质方面发挥重要的作用.细胞色素P4502C8是人体肝脏中主要负责代谢抗癌药物紫杉醇的酶,它至少负责代谢5%的临床药物.细胞色素P450 2C8的基因多态性与用药个体化有着密切的关系.CYP2C8.3是常见的P450 2C8的基因多态之一,其发生了双点突变,分别是R139K...     

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.     

Characterization of the CYP3A4 active site by homology modeling

Human microsomal cytochrome P450s participate in drug metabolism and detoxification. Among them, CYP3A4 is the most important isoform for drug-drug interactions. To gain a better understanding of the active site, a homology model of CYP3A4 was constructed based on the crystallographic coordinates of mammalian CYP2C5. The putative active site is much larger than that of CYP2C5 and is divided into three parts (i.e. a proximal and two distal sites from the heme). Most residues reported to be important for ligand-binding are located in the active site of the model. Moreover, some inhibitors (paclitaxel etc.) docked into the model have complementary shapes to the pocket. Pharmacophore docking of 14 substrates was also performed using Ph4Dock of MOE. Calculated interaction energies showed a moderate correlation with the logarithm of apparent K(m) values. These results suggest that this model is reliable enough to be used in the design of compounds for removing undesirable CYP3A4 inhibition.     

Effective virtual screening protocol for CYP2C9 ligands using a screening site constructed from flurbiprofen and S-warfarin pockets

An effective virtual screening protocol was developed against an extended active site of CYP2C9, which was derived from X-ray structures complexed with flubiprofen and S-warfarin. Virtual screening has been effectively supported by our structure-based pharmacophore model. Importance of hot residues identified by mutation data and structural analysis was first estimated in an enrichment study. Key role of Arg108 and Phe114 in ligand binding was also underlined. Our screening protocol successfully identified 76% of known CYP2C9 ligands in the top 1% of the ranked database resulting 76-fold enrichment relative to random situation. Relevance of the protocol was further confirmed in selectivity studies, when 89% of CYP2C9 ligands were retrieved from a mixture of CYP2C9 and CYP2C8 ligands, while only 22% of CYP2C8 ligands were found applying the structure-based pharmacophore constraints. Moderate discrimination of CYP2C9 ligands from CYP2C18 and CYP2C19 ligands could also be achieved extending the application domain of our virtual screening protocol for the entire CYP2C family. Our findings further demonstrate the existence of an active site comprising of at least two binding pockets and strengthens the need of involvement of protein flexibility in virtual screening.     

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.     

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.     

Ligand- and protein-based modeling studies of the inhibitors of human cytochrome P450 2D6 and a virtual screening for potential inhibitors from the Chinese herbal medicine, Scutellaria baicalensis (Huangqin,Baikal Skullcap)

We have previously examined the binding patterns of various substrates to human cytochrome P450 2D6 (CYP2D6) using a series of molecular modeling methods. In this study, we further explored the binding modes of various types of inhibitors to CYP2D6 using a combination of ligand- and protein-based modeling approaches. Firstly, we developed and validated a pharmacophore model for CYP2D6 inhibitors, which consisted of two hydrophobic features and one hydrogen bond acceptor feature. Secondly, we constructed and validated a quantitative structure-activity relationship (QSAR) model for CYP2D6 inhibitors which gave a poor to moderate prediction accuracy. Thirdly, a panel of CYP2D6 inhibitors were subject to molecular docking into the active site of wild-type and mutated CYP2D6 enzyme. We demonstrated that 8 residues in the active site (Leu213, Glu216, Ser217, Gln244, Asp301, Ser304, Ala305, and Phe483) played an important role in the binding to the inhibitors via hydrogen bond formation and/or π-π stacking interaction. Apparent changes in the binding modes of the inhibitors have been observed with Phe120Ile, Glu216Asp, Asp301Glu mutations in CYP2D6. Finally, we screened for potential binders/inhibitors from the Chinese herbal medicine Scutellaria baicalensis (Huangqin, Baikal Skullcap) using the established pharmacophore model for CYP2D6 inhibitors and molecular docking approach. Overall, 18 out of 40 compounds from S. baicalensis were mapped to the pharmacophore model of CYP2D6 inhibitors and most herbal compounds from S. baicalensis could be docked into the active site of CYP2D6. Our study has provided insights into the molecular mechanisms of interaction of synthetic and herbal compounds with human CYP2D6 and further benchmarking studies are needed to validate our modeling and virtual screening results.     

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.     

Pyrene.pyrene complexes at the active site of cytochrome P450 3A4: evidence for a multiple substrate binding site

Cytochrome P450 monooxygenases (CYPs) metabolize nearly all drugs and toxins. Recently, it has become clear that CYPs exhibit both homotropic and heterotropic allosteric kinetics for many substrates. However, the mechanism of cooperative kinetics has not been established for any specific human CYP/substrate combination. Suggested mechanisms include binding of multiple substrates within distinct, static, subsites of a single large active site or binding of multiple substrates within a single fluid active site. CYP3A4 hydroxylates pyrene with positive cooperativity. Therefore, experiments were designed to exploit the fluorescence properties of pyrene, which diagnostically distinguish between pyrene.pyrene complexes versus spatially separated pyrene substrates. Pyrene complexes (excimers) yield an emission spectrum clearly distinct from pyrene monomers. In lipid-free aqueous/glycerol solutions of CYP3A4, addition of pyrene affords a concentration-dependent low-spin to high-spin conversion of the CYP3A4 heme prosthetic group, indicating occupancy of the active site by pyrene. Under the same conditions, in the presence of CYP3A4 but not other heme proteins, the excimer/monomer ratio (E/M) of pyrene was decreased in emission spectra, compared to pyrene alone. However, excitation spectra indicate a CYP3A4-dependent increase in the wavelength shift for the excimer excitation spectrum versus the monomer excitation spectrum, as well as changes in the excimer excitation peak shape and vibronic structure. These changes are reversed by the CYP3A4 substrate testosterone. Together, the results demonstrate that pyrene.pyrene ground-state complexes occupy the CYP3A4 active site, and they provide the first spectroscopic evidence for substrate complexes within a single fluid active site. Functional implications include the possibility that turnover rate, regioselectivity, and stereoselectivity of the reaction are determined by the substrate.substrate complex rather than individual substrates.     

Amino acids of ricin and its polypeptides

Ricin and its corresponding polypeptides (A & B chain) were purified from castor seed. The molecular weight of ricin subunits were 29,000 and 28,000 daltons. The amino acids in ricin determined were Asp45 The22 Ser40 Glu53 Cys4 Gly96 His5 Ile21 Leu33 Lys20 Met4 Phe13 Pro37 Tyr11 Ala45 Val23 Arg20 indicating that ricin contains approximately 516 amino acid residues. The amino acids of the two subunits of ricin A and B chains were Asp23 The12 Ser21 Glu29 Cys2 Gly48 His3 Ile12, Leu17 Lys10 Met2 Phe6 Pro17 Tyr7 Ala35 Val13 Arg13 while in B chain the amino acids were Asp22 The10 Ser19 Glu25 Cys2 Gly47 His1 Ile10, Leu15 Lys11 Met1 Phe7 Pro6 Tyr5 Ala32Val11 Arg10. The total helical content of ricin came around 53.6% which is a new observation.     

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.     

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.     

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.     

Structure-Activity relationship of MDSA and its derivatives against Staphylococcus aureus Ser/Thr phosphatase Stp1

Ser/thr phosphatase Stp1 is an important virulence factor for Staphylococcus aureus (S. aureus) and plays a key role in its infectivity, suggesting that it could serve as a potential target for treatment of S. aureus infection. Previous studies found that the activity of Stp1 was inhibited by MDSA and its derivatives. In this paper, we used molecular docking, molecular modeling, molecular dynamics simulations, binding free energy decomposition calculations, and hydrogen bond analyses to explore the structure-activity relationship. Energy decomposition indicated that MDSA, hydroxymethyl MDSA, carboxymethyl MDSA and methyl MDSA can bind to the catalytic pocket of Stp1. Furthermore, Met39, Ile163, Ile164, Val167, Gly195 and Asp233 were key residues in the Stp1-inhibitor complexes. Due to the lack of a double salicylate structure, salicylic acid cannot bind to the active site of Stpl, leading to loss of inhibitory activity. Based on these results, the structure-activity relationship at the atomic level was determined, which can promote the development of new and more effective anti-drug resistance inhibitors.     

Fourier transform infrared spectroscopy as a tool to study structural properties of cytochromes P450 (CYPs)

Cytochrome P450 proteins (CYPs) are a big class of heme proteins which are involved in various metabolic processes of living organisms. CYPs are the terminal catalytically active components of monooxygenase systems where the substrate binds and is hydroxylated. In order to be functionally competent, the protein structures of CYPs possess specific properties that must be explored in order to understand structure–function relationships and mechanistic aspects. Fourier transform infrared spectroscopy (FTIR) is one tool that is used to study these structural properties. The application of FTIR spectroscopy to the secondary structures of CYP proteins, protein unfolding, protein–protein interactions and the structure and dynamics of the CYP heme pocket is reviewed. A comparison with other thiolate heme proteins (nitric oxide synthase and chloroperoxidase) is also included. Figure The protein secondary structure, protein unfolding, redox-partner protein–protein interaction, structural changes induced by the reduction of the heme iron, and the structure and dynamics of the active site of cytochromes P450 (CYP) can be studied using Fourier transform infrared spectroscopy (FTIR). FTIR spectroscopy is a good approach for gaining a deeper insight into structure–function relationships in CYPs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.