Homology modeling of the three membrane proteins of the dhurrin metabolon: catalytic sites, membrane surface association and protein-protein interactions

Formation of metabolons (macromolecular enzyme complexes) facilitates the channelling of substrates in biosynthetic pathways. Metabolon formation is a dynamic process in which transient structures mediated by weak protein-protein interactions are formed. In Sorghum, the cyanogenic glucoside dhurrin is derived from l-tyrosine in a pathway involving the two cytochromes P450 (CYPs) CYP79A1 and CYP71E1, a glucosyltransferase (UGT85B1), and the redox partner NADPH-dependent cytochrome P450 reductase (CPR). Experimental evidence suggests that the enzymes of this pathway form a metabolon. Homology modeling of the three membrane bound proteins was carried out using the Sybyl software and available relevant crystal structures. Residues involved in tight positioning of the substrates and intermediates in the active sites of CYP79A1 and CYP71E1 were identified. In both CYPs, hydrophobic surface domains close to the N-terminal trans-membrane anchor and between the F′ and G helices were identified as involved in membrane anchoring. The proximal surface of both CYPs showed positively charged patches complementary to a negatively charged bulge on CPR carrying the FMN domain. A patch of surface exposed, positively charged amino acid residues positioned on the opposite face of the membrane anchor was identified in CYP71E1 and might be involved in binding UGT85B1 via a hypervariable negatively charged loop in this protein.     

The signal-anchor sequence of CYP2C1 inserts into the membrane as a hairpin structure

Microsomal cytochrome P450s (CYPs) are anchored to the endoplasmic reticulum membrane by the N-terminal signal-anchor sequence which is predicted to insert into the membrane as a type 1 transmembrane helix with a luminally located N-terminus. We have mapped amino acids of the CYP2C1 signal-anchor, fused to Cys-free glutathione S-transferase, within the membrane by Cys-specific labeling with membrane-impermeant maleimide polyethylene glycol. At the C-terminal end of the signal-anchor, Trp-20 was mapped to the membrane–cytosol interface and Leu-19 was within the membrane. Unexpectedly, at the N-terminal end, Glu-2 and Pro-3 were mapped to the cytoplasmic side of the membrane rather than the luminal side as expected of a type 1 transmembrane helix. Similar results were observed for the N-terminal amino acids of the signal-anchor sequences of CYP3A4 and CYP2E1. These observations indicate that contrary to the current model of the signal-anchor of CYPs as a type 1 transmembrane helix, CYP2C1, CYP2E1, and CYP3A4 are monotopic membrane proteins with N-terminal signal-anchors that have a hairpin or wedge orientation in the membrane.     

Antibody binding modulates the dynamics of the membrane-bound prion protein

Misfolding of the cellular prion protein (PrP^C) is associated with lethal neurodegeneration. PrP^C consists of a flexible tail (residues 23–123) and a globular domain (residues 124–231) whose C-terminal end is anchored to the cell membrane. The neurotoxic antibody POM1 and the innocuous antibody POM6 recognize the globular domain. Experimental evidence indicates that POM1 binding to PrP^C emulates the influence on PrP^C of the misfolded prion protein (PrP^Sc) while the binding of POM6 has the opposite biological response. Little is known about the potential interactions between flexible tail, globular domain, and the membrane. Here, we used atomistic simulations to investigate how these interactions are modulated by the binding of the Fab fragments of POM1 and POM6 to PrP^C and by interstitial sequence truncations to the flexible tail. The simulations show that the binding of the antibodies restricts the range of orientations of the globular domain with respect to the membrane and decreases the distance between tail and membrane. Five of the six sequence truncations influence only marginally this distance and the contact patterns between tail and globular domain. The only exception is a truncation coupled to a charge inversion mutation of four N-terminal residues, which increases the distance of the flexible tail from the membrane. The interactions of the flexible tail and globular domain are modulated differently by the two antibodies.     

Inhibition of cytochrome P450 activities by oleanolic acid and ursolic acid in human liver microsomes

Oleanolic acid (OA) and ursolic acid (UA), triterpene acids having numerous pharmacological activities including anti-inflammatory, anti-cancer, and hepato-protective effects, were tested for their ability to modulate the activities of several cytochrome P450 (CYP) enzymes using human liver microsomes. OA competitively inhibited CYP1A2-catalyzed phenacetin O-deethylation and CYP3A4-catalyzed midazolam 1-hydroxylation, the major human drug metabolizing CYPs, with IC50 (Ki) values of 143.5 (74.2) microM and 78.9 (41.0) microM, respectively. UA competitively inhibited CYP2C19-catalyzed S-mephenytoin 4'-hydroxylation with an IC50 (Ki) value of 119.7 (80.3) microM. However, other CYPs tested showed no or weak inhibition by both OA and UA. The present study demonstrates that OA and UA have inhibitory effects on CYP isoforms using human liver microsomes. It is thus likely that consumption of herbal medicines containing OA or UA, or administration of OA or UA, can cause drug interactions in humans when used concomitantly with drugs that are metabolized primarily by CYP isoforms. In addition, it appears that the inhibitory effect of OA on CYP1A2 is, in part, related to its anti-inflammatory and anticancer activities.     

Polymorphic Cytochrome P450 Enzymes (CYPs) and Their Role in Personalized Therapy

The cytochrome P450 (CYP) enzymes are major players in drug metabolism. More than 2,000 mutations have been described, and certain single nucleotide polymorphisms (SNPs) have been shown to have a large impact on CYP activity. Therefore, CYPs play an important role in inter-individual drug response and their genetic variability should be factored into personalized medicine. To identify the most relevant polymorphisms in human CYPs, a text mining approach was used. We investigated their frequencies in different ethnic groups, the number of drugs that are metabolized by each CYP, the impact of CYP SNPs, as well as CYP expression patterns in different tissues. The most important polymorphic CYPs were found to be 1A2, 2D6, 2C9 and 2C19. Thirty-four common allele variants in Caucasians led to altered enzyme activity. To compare the relevant Caucasian SNPs with those of other ethnicities a search in 1,000 individual genomes was undertaken. We found 199 non-synonymous SNPs with frequencies over one percent in the 1,000 genomes, many of them not described so far. With knowledge of frequent mutations and their impact on CYP activities, it may be possible to predict patient response to certain drugs, as well as adverse side effects. With improved availability of genotyping, our data may provide a resource for an understanding of the effects of specific SNPs in CYPs, enabling the selection of a more personalized treatment regimen.     

Dynamics of CYP51: implications for function and inhibitor design

Sterol 14α‐demethylase (cytochrome P450 family 51 (CYP51)) is an essential enzyme occurring in all biological kingdoms. In eukaryotes, it is located in the membrane of the endoplasmic reticulum. Selective inhibitors of trypanosomal CYP51s that do not affect the human CYP51 have been discovered in vitro and found to cure acute and chronic mouse Chagas disease without severe side effects in vivo. Crystal structures indicate that CYP51 may be more rigid than most CYPs, and it has been proposed that this property may facilitate antiparasitic drug design. Therefore, to investigate the dynamics of trypanosomal CYP51, we built a model of membrane‐bound Trypanosoma brucei CYP51 and then performed molecular dynamics simulations of T. brucei CYP51 in membrane‐bound and soluble forms. We compared the dynamics of T. brucei CYP51 with those of human CYP51, CYP2C9, and CYP2E1. In the simulations, the CYP51s display low mobility in the buried active site although overall mobility is similar in all the CYPs studied. The simulations suggest that in CYP51, pathway 2f serves as the major ligand access tunnel, and both pathways 2f (leading to membrane) and S (leading to solvent) can serve as ligand egress tunnels. Compared with the other CYPs, the residues at the entrance of the ligand access tunnels in CYP51 have higher mobility that may be necessary to facilitate the passage of its large sterol ligands. The water (W) tunnel is accessible to solvent during most of the simulations of CYP51, but its width is affected by the conformations of the heme's two propionate groups. These differ from those observed in the other CYPs studied because of differences in their hydrogen‐bonding network. Our simulations give insights into the dynamics of CYP51 that complement the available experimental data and have implications for drug design against CYP51 enzymes. Copyright © 2015 John Wiley & Sons, Ltd.     

Cytochrome allelic variants and clopidogrel metabolism in cardiovascular diseases therapy

Clopidogrel and aspirin are among the most prescribed dual antiplatelet therapies to treat the acute coronary syndrome and heart attacks. However, their potential clinical impacts are a subject of intense debates. The therapeutic efficiency of clopidogrel is controlled by the actions of hepatic cytochrome P450 (CYPs) enzymes and impacted by individual genetic variations. Inter-individual polymorphisms in CYPs enzymes affect the metabolism of clopidogrel into its active metabolites and, therefore, modify its turnover and clinical outcome. So far, clinical trials fail to confirm higher or lower adverse cardiovascular effects in patients treated with combinations of clopidogrel and proton pump inhibitors, compared with clopidogrel alone. Such inconclusive findings may be due to genetic variations in the cytochromes CYP2C19 and CYP3A4/5. To investigate potential interactions/effects of these cytochromes and their allele variants on the treatment of acute coronary syndrome with clopidogrel alone or in combination with proton pump inhibitors, we analyze recent literature and discuss the potential impact of the cytochrome allelic variants on cardiovascular events and stent thrombosis treated with clopidogrel. The diversity of CYP2C19 polymorphisms and prevalence span within various ethnic groups, subpopulations and demographic areas are also debated.     

An unusual TOM20/TOM22 bypass mechanism for the mitochondrial targeting of cytochrome P450 proteins containing N-terminal chimeric signals

Previously we showed that xenobiotic-inducible cytochrome P450 (CYP) proteins are bimodally targeted to the endoplasmic reticulum and mitochondria. In the present study, we investigated the mechanism of delivery of chimeric signal-containing CYP proteins to the peripheral and channel-forming mitochondrial outer membrane translocases (TOMs). CYP+33/1A1 and CYP2B1 did not require peripheral TOM70, TOM20, or TOM22 for translocation through the channel-forming TOM40 protein. In contrast, CYP+5/1A1 and CYP2E1 were able to bypass TOM20 and TOM22 but required TOM70. CYP27, which contains a canonical cleavable mitochondrial signal, required all of the peripheral TOMs for its mitochondrial translocation. We investigated the underlying mechanisms of bypass of peripheral TOMs by CYPs with chimeric signals. The results suggested that interaction of CYPs with Hsp70, a cytosolic chaperone involved in the mitochondrial import, alone was sufficient for the recognition of chimeric signals by peripheral TOMs. However, sequential interaction of chimeric signal-containing CYPs with Hsp70 and Hsp90 resulted in the bypass of peripheral TOMs, whereas CYP27 interacted only with Hsp70 and was not able to bypass peripheral TOMs. Our results also show that delivery of chimeric signal-containing client proteins by Hsp90 required the cytosol-exposed N-terminal 143 amino acids of TOM40. TOM40 devoid of this domain was unable to bind CYP proteins. These results suggest that, compared with the unimodal mitochondria-targeting signals, the chimeric mitochondria-targeting signals are highly evolved and dynamic in nature.     

Targeting cytochrome P450 enzymes: a new approach in anti-cancer drug development

Cytochrome P450s (CYPs) represent a large class of heme-containing enzymes that catalyze the metabolism of multitudes of substrates both endogenous and exogenous. Until recently, however, CYPs have been largely overlooked in cancer drug development, acknowledged only for their role in phase I metabolism of chemotherapeutics. The first successful strategy targeting CYP enzymes in cancer therapy was the development of potent inhibitors of CYP19 (aromatase) for the treatment of breast cancer. Aromatase inhibitors ushered in a new era in hormone ablation therapy for estrogen dependent cancers, and have paved the way for similar strategies (i.e., inhibition of CYP17) that combat androgen dependent prostate cancer. Identification of CYPs involved in the inactivation of anti-cancer metabolites of vitamin D(3) and vitamin A has triggered development of agents that target these enzymes as well. The discovery of the over-expression of exogenous metabolizing CYPs, such as CYP1B1, in cancer cells has roused interest in the development of inhibitors for chemoprevention and of prodrugs designed to be activated by CYPs only in cancer cells. Finally, the expression of CYPs within tumors has been utilized in the development of bioreductive molecules that are activated by CYPs only under hypoxic conditions. This review offers the first comprehensive analysis of strategies in drug development that either inhibit or exploit CYP enzymes for the treatment of cancer.     

Antitumor drug ellipticine inhibits the activities of rat hepatic cytochromes P450

Ellipticine is a potent antineoplastic agent, whose mode of action is considered to be based mainly on DNA intercalation and/or inhibition of topoisomerase II. Recently, we found that ellipticine also forms the cytochrome P450 (CYP)-mediated covalent DNA adducts. Here, we study the effect of ellipticine on CYP enzymes in rat hepatic microsomes, studying its binding to the enzymes and its potential to inhibit the CYP activities measured with their selective substrates. Although ellipticine was reported to be a selective and strong inhibitor of CYP1A1/2, we found that its inhibitory potential is non-specific. Ellipticine is the most potent inhibitor for CYP3A-dependent 6beta-hydroxylation of progesterone, followed by CYP1A1/2-dependent ethoxyresorufin O-deethylation and CYP2B-mediated pentoxyresorufin O-depentylation. Lower inhibition was detected for 1'-hydroxylation of bufurarol, 21-hydroxylation of progesterone and 6-hydroxylation of chlorzoxazone catalyzed by CYP2D, CYP2C and CYP2E1, respectively. Ellipticine binds to several CYPs of rat hepatic microsomes. The binding titration of ellipticine typically give reverse type I spectrum with CYPs in rat hepatic microsomes. The results indicate that inhibition of CYPs by ellipticine cannot be explained only by its differential potency to bind to individual CYPs.     

Searching the Cytochrome P450 Enzymes for the Metabolism of Meranzin Hydrate: A Prospective Antidepressant Originating from Chaihu-Shugan-San

Meranzin hydrate (MH), an absorbed bioactive compound from the Traditional Chinese Medicine (TCM) Chaihu-Shugan-San (CSS), was first isolated in our laboratory and was found to possess anti-depression activity. However, the role of cytochrome P450s (CYPs) in the metabolism of MH was unclear. In this study, we screened the CYPs for the metabolism of MH in vitro by human liver microsomes (HLMs) or human recombinant CYPs. MH inhibited the enzyme activities of CYP1A2 and CYP2C19 in a concentration-dependent manner in the HLMs. The K_m and V_max values of MH were 10.3±1.3 µM and 99.1±3.3 nmol/mg protein/min, respectively, for the HLMs; 8.0±1.6 µM and 112.4±5.7 nmol/nmol P450/min, respectively, for CYP1A2; and 25.9±6.6 µM and 134.3±12.4 nmol/nmol P450/min, respectively, for CYP2C19. Other human CYP isoforms including CYP2A6, CYP2C9, CYP2D6, CYP2E1 and CYP3A4 showed minimal or no effect on MH metabolism. The results suggested that MH was simultaneously a substrate and an inhibitor of CYP1A2 and CYP2C9, and MH had the potential to perpetrate drug-drug interactions with other CYP1A2 and CYP2C19 substrates.     

Structure and dynamics of the membrane-bound cytochrome P450 2C9

The microsomal, membrane-bound, human cytochrome P450 (CYP) 2C9 is a liver-specific monooxygenase essential for drug metabolism. CYPs require electron transfer from the membrane-bound CYP reductase (CPR) for catalysis. The structural details and functional relevance of the CYP-membrane interaction are not understood. From multiple coarse grained molecular simulations started with arbitrary configurations of protein-membrane complexes, we found two predominant orientations of CYP2C9 in the membrane, both consistent with experiments and conserved in atomic-resolution simulations. The dynamics of membrane-bound and soluble CYP2C9 revealed correlations between opening and closing of different tunnels from the enzyme's buried active site. The membrane facilitated the opening of a tunnel leading into it by stabilizing the open state of an internal aromatic gate. Other tunnels opened selectively in the simulations of product-bound CYP2C9. We propose that the membrane promotes binding of liposoluble substrates by stabilizing protein conformations with an open access tunnel and provide evidence for selective substrate access and product release routes in mammalian CYPs. The models derived here are suitable for extension to incorporate other CYPs for oligomerization studies or the CYP reductase for studies of the electron transfer mechanism, whereas the modeling procedure is generally applicable to study proteins anchored in the bilayer by a single transmembrane helix.     

Differential inhibition of CYP17A1 and CYP21A2 activities by the P450 oxidoreductase mutant A287P

P450 oxidoreductase (POR) has a pivotal role in facilitating electron transfer from nicotinamide adenine dinucleotide phosphate to microsomal cytochrome P450 (CYP) enzymes, including the steroidogenic enzymes CYP17A1 and CYP21A2. Mutations in POR have been shown recently to cause congenital adrenal hyperplasia with apparent combined CYP17A1 and CYP21A2 deficiency that comprises a variable clinical phenotype, including glucocorticoid deficiency, ambiguous genitalia, and craniofacial malformations. To dissect structure-function relationships potentially explaining this phenotypic diversity, we investigated whether specific POR mutations have differential effects on CYP17A1 and CYP21A2. We compared the impact of missense mutations encoding for single amino acid changes in three distinct regions of the POR molecule: 1), Y181D and H628P close to the central electron transfer area, 2) S244C located within the hinge close to the flavin adenine dinucleotide and flavin mononucleotide domains of POR, and 3) A287P that is clearly distant from the two other regions. Functional analysis using a yeast microsomal assay with coexpression of human CYP17A1 or CYP21A2 with wild-type or mutant human POR revealed equivalent decreases in CYP17A1 and CYP21A2 activities by Y181D, H628P, and S244C. In contrast, A287P had a differential inhibitory effect, with decreased catalytic efficiency (Vmax/Km) for CYP17A1, whereas CYP21A2 retained near normal activity. In vivo analysis of urinary steroid excretion by gas chromatography/mass spectrometry in 11 patients with POR mutations showed that A287P homozygous patients had the highest corticosterone/cortisol metabolite ratios, further indicative of preferential inhibition of CYP17A1. These findings provide novel mechanistic insights into the redox regulation of human steroidogenesis. Differential interaction of POR with electron-accepting CYP enzymes may explain the phenotypic variability in POR deficiency, with additional implications for hepatic drug metabolism by POR-dependant CYP enzymes.     

Identifying Cytochrome P450 Functional Networks and Their Allosteric Regulatory Elements

Cytochrome P450 (CYP) enzymes play key roles in drug metabolism and adverse drug-drug interactions. Despite tremendous efforts in the past decades, essential questions regarding the function and activity of CYPs remain unanswered. Here, we used a combination of sequence-based co-evolutionary analysis and structure-based anisotropic thermal diffusion (ATD) molecular dynamics simulations to detect allosteric networks of amino acid residues and characterize their biological and molecular functions. We investigated four CYP subfamilies (CYP1A, CYP2D, CYP2C, and CYP3A) that are involved in 90% of all metabolic drug transformations and identified four amino acid interaction networks associated with specific CYP functionalities, i.e., membrane binding, heme binding, catalytic activity, and dimerization. Interestingly, we did not detect any co-evolved substrate-binding network, suggesting that substrate recognition is specific for each subfamily. Analysis of the membrane binding networks revealed that different CYP proteins adopt different membrane-bound orientations, consistent with the differing substrate preference for each isoform. The catalytic networks were associated with conservation of catalytic function among CYP isoforms, whereas the dimerization network was specific to different CYP isoforms. We further applied low-temperature ATD simulations to verify proposed allosteric sites associated with the heme-binding network and their role in regulating metabolic fate. Our approach allowed for a broad characterization of CYP properties, such as membrane interactions, catalytic mechanisms, dimerization, and linking these to groups of residues that can serve as allosteric regulators. The presented combined co-evolutionary analysis and ATD simulation approach is also generally applicable to other biological systems where allostery plays a role.     

Carcinogenic pollutants o-nitroanisole and o-anisidine are substrates and inducers of cytochromes P450

2-Methoxyaniline (o-anisidine) and 2-methoxynitrobenzene (o-nitroanisole) are important pollutants and potent carcinogens for rodents. o-Anisidine is oxidized by microsomes of rats and rabbits to N-(2-methoxyphenyl)hydroxylamine that is also formed as the reduction metabolite of o-nitroanisole. o-Anisidine is a promiscuity substrate of rat and rabbit cytochrome P450 (CYP) enzymes, because CYPs of 1A, 2B, 2E and 3A subfamilies oxidize o-anisidine. Using purified CYP enzymes, reconstituted with NADPH: CYP reductase, rabbit CYP2E1 was the most efficient enzyme oxidizing o-anisidine, but the ability of CYP1A1, 1A2, 2B2, 2B4 and 3A6 to participate in o-anisidine oxidation was also proved. Utilizing Western blotting and consecutive immunoquantification employing chicken polyclonal anti bodies raised against various CYPs, the effect of o-anisidine and o-nitroanisole on the expression of the CYP enzymes was investigated. The expression of CYP1A1/2 was found to be strongly induced in rats treated with either compounds. In addition, 7-ethoxyresorufin O-deethylation, a marker activity for both CYP1A1 and 1A2, was significantly increased in rats treated with either carcinogen. The data demonstrate the participation of different rat and rabbit CYP enzymes in o-anisidine oxidation and indicate that both experimental animal species might serve as suitable models to mimic the o-anisidine oxidation in human. Furthermore, by induction of rat hepatic and renal CYP1A1/2, both o-nitroanisole and o-anisidine influence their carcinogenic effects, modifying their detoxification and/or activation pathways.     

Comparison of the 1.85 A structure of CYP154A1 from Streptomyces coelicolor A3(2) with the closely related CYP154C1 and CYPs from antibiotic biosynthetic pathways

The genus Streptomyces produces two-thirds of microbially derived antibiotics. Polyketides form the largest and most diverse group of these natural products. Antibiotic diversity of polyketides is generated during their biosynthesis by several means, including postpolyketide modification performed by oxidoreductases, a broad group of enzymes including cytochrome P450 monooxygenases (CYPs). CYPs catalyze site-specific oxidation of macrolide antibiotic precursors significantly affecting antibiotic activity. Efficient manipulation of Streptomyces CYPs in generating new antibiotics will require identification and/or engineering of monooxygenases with activities toward a diverse array of chemical substrates. To begin to link structure to function of CYPs involved in secondary metabolic pathways of industrially important species, we determined the X-ray structure of Streptomyces coelicolor A3(2) CYP154A1 at 1.85 A and analyzed it in the context of the closely related CYP154C1 and more distant CYPs from polyketide synthase (EryF) and nonribosomal peptide synthetase (OxyB) biosynthetic pathways. In contrast to CYP154C1, CYP154A1 reveals an active site inaccessible from the molecular surface, and an absence of catalytic activities observed for CYP154C1. Systematic variations in the amino acid patterns and length of the surface HI loop correlate with degree of rotation of the F and G helices relative to the active site in CYP154A1-related CYPs, presumably regulating the degree of active site accessibility and its dimensions. Heme in CYP154A1 is in a 180 degrees flipped orientation compared with most other structurally determined CYPs.     

Identification of universal selectivity‐determining positions in cytochrome P450 monooxygenases by systematic sequence‐based literature mining

Cytochrome P450 monooxygenases (CYPs) are a large, highly diverse protein family with a common fold. The sequences, structures, and functions of CYPs have been extensively studied resulting in more than 53,000 scientific articles. A sequence‐based literature mining algorithm was designed to systematically analyze this wealth of information on SNPs, designed mutations, structural interactions, or functional roles of individual residues. Structurally corresponding positions in different CYPs were compared and universal selectivity‐determining positions were identified. Based on the Cytochrome P450 Engineering Database ( www.CYPED.BioCatNet.de ) and a standard numbering scheme for all CYPs, 4000 residues in 168 CYPs mentioned in 2400 articles could be assigned to 440 structurally corresponding standard positions of the CYP fold, covering 96% of all standard positions. Seventeen individual standard positions were mentioned in the context of more than 32 different CYPs. The majority of these most frequently mentioned positions are located on the six substrate recognition sites and are involved in control of selectivity, such as the well‐studied position 87 in CYP102A1 (P450_BM‐3) which was mentioned in the articles on 63 different CYPs. The recurrent citation of the 17 frequently mentioned positions for different CYPs suggests their universal functional relevance. Proteins 2015; 83:1593–1603. © 2015 Wiley Periodicals, Inc.     

SDS‐PAGE‐free protocol for comprehensive identification of cytochrome P450 enzymes and uridine diphosphoglucuronosyl transferases in human liver microsomes

Comprehensive identification of cytochrome P450 enzymes (CYPs) and uridine diphosphoglucuronosyl transferases (UGTs) in human liver microsomes (HLMs) was performed with an SDS‐PAGE‐free protocol. HLMs were solubilized with 5% v/v ionic liquid, 1‐butyl‐3‐methyl imidazolium tetrafluoroborate, followed by tryptic digestion, and 2D‐SCX‐RPLC‐ESI‐MS/MS (LTQ XL) analysis in triplicate. In total, 27 CYPs and 12 UGTs were confidently identified with average sequence coverage as 30.99 and 25.07%, average peptide number as 14 and 13, and average unique peptide number as 7 and 4, respectively. The highly similar isoforms of CYP3A, CYP2C, and CYP4F subfamilies could be unambiguously differentiated from each other, despite the fact that the sequence similarity of CYP2C9 and CYP2C19 is 91%. In addition, protein spectral count was used to approximately evaluate the relative abundance of identified CYPs and UGTs, and the results agreed with previous immunochemistry reports.