Identification of residues 99, 220, and 221 of human cytochrome P450 2C19 as key determinants of omeprazole activity

Deep mutational scanning of CYP2C19 in human cells reveals a substrate specificity-abundance tradeoff

The cytochrome P450s enzyme family metabolizes ∼80% of small molecule drugs. Variants in cytochrome P450s can substantially alter drug metabolism, leading to improper dosing and severe adverse drug reactions. Due to low sequence conservation, predicting variant effects across cytochrome P450s is challenging. Even closely related cytochrome P450s like CYP2C9 and CYP2C19, which share 92% amino acid sequence identity, display distinct phenotypic properties. Using variant abundance by massively parallel sequencing, we measured the steady-state protein abundance of 7,660 single amino acid variants in CYP2C19 expressed in cultured human cells. Our findings confirmed critical positions and structural features essential for cytochrome P450 function, and revealed how variants at conserved positions influence abundance. We jointly analyzed 4,670 variants whose abundance was measured in both CYP2C19 and CYP2C9, finding that the homologs have different variant abundances in substrate recognition sites within the hydrophobic core. We also measured the abundance of all single and some multiple wild type amino acid exchanges between CYP2C19 and CYP2C9. While most exchanges had no effect, substitutions in substrate recognition site 4 reduced abundance in CYP2C19. Double and triple mutants showed distinct interactions, highlighting a region that points to differing thermodynamic properties between the 2 homologs. These positions are known contributors to substrate specificity, suggesting an evolutionary tradeoff between stability and enzymatic function. Finally, we analyzed 368 previously unannotated human variants, finding that 43% had decreased abundance. By comparing variant effects between these homologs, we uncovered regions underlying their functional differences, advancing our understanding of this versatile family of enzymes.

Identification of cytochrome P450 2C18 and 2C76 in tree shrews: P450 2C18 effectively oxidizes typical human P450 2C9/2C19 chiral substrates warfarin and omeprazole with less stereoselectivity

Cytochromes P450 (P450s or CYPs), especially the CYP2C family, are important drug-metabolizing enzymes that play major roles in drug metabolism. Tree shrews, a non-rodent primate-like species, are used in various fields of biomedical research, notably hepatitis virus infection; however, its drug-metabolizing enzymes have not been fully investigated. In this study, tree shrew CYP2C18, CYP2C76a, CYP2C76b, and CYP2C76c cDNAs were identified and contained open reading frames of 489 or 490 amino acids with high sequence identities (70-78 %) to human CYP2Cs. Tree shrew CYP2C76a, CYP2C76b, and CYP2C76c showed higher sequence identities (79-80 %) to cynomolgus CYP2C76 and were not orthologous to any human CYP2C. Phylogenetic analysis revealed that tree shrew CYP2C18 and CYP2C76s were closely related to rat CYP2Cs and cynomolgus CYP2C76, respectively. Tree shrew CYP2C genes formed a gene cluster similar to human CYP2C genes. All four tree shrew CYP2C mRNAs showed predominant expressions in liver, among the tissue types examined; expression of CYP2C18 mRNA was also detected in small intestine. In liver, CYP2C18 mRNA was the most abundant among the tree shrew CYP2C mRNAs. In metabolic assays using human CYP2C substrates, all tree shrew CYP2Cs showed metabolic activities toward diclofenac, R,S-omeprazole, paclitaxel, and R,S-warfarin, with the activity of CYP2C18 exceeding that of the other CYP2Cs. Moreover, tree shrew CYP2C76 enzymes metabolized progesterone more efficiently than human, cynomolgus, or marmoset CYP2Cs. Therefore, these novel tree shrew CYP2Cs are expressed abundantly in liver, encode functional enzymes that metabolize human CYP2C substrates, and are likely responsible for drug clearances.

Novel Cytochrome P450 2C119 Enzymes in Cynomolgus and Rhesus Macaques Metabolize Progesterone, Diclofenac, and Omeprazole

Cynomolgus and rhesus macaques are used in drug metabolism studies due to their evolutionary and phylogenetic closeness to humans. Cytochromes P450 (P450s or CYPs), including the CYP2C family enzyme, are important endogenous and exogenous substrate-metabolizing enzymes and play major roles in drug metabolism. In cynomolgus and rhesus macaques, six CYP2Cs have been identified and characterized, namely, CYP2C8, CYP2C9, CYP2C18, CYP2C19, CYP2C76, and CYP2C93. In this study, CYP2C119, a new CYP2C, was identified and characterized in cynomolgus and rhesus macaques. Cynomolgus and rhesus CYP2C119 contained open reading frames of 489 amino acids with high sequence identities to human CYP2C8 and to cynomolgus and rhesus CYP2C8. Phylogenetic analysis showed that cynomolgus and rhesus CYP2C119 were closely related to cynomolgus and rhesus CYP2C8. In cynomolgus and rhesus genomes, CYP2C genes, including CYP2C119, form a cluster. Among the tissues analyzed, cynomolgus CYP2C119 mRNA was predominantly expressed in liver. Hepatic expressions of CYP2C119 mRNA in four cynomolgus and two rhesus macaques varied, with no expression in one rhesus macaque. Among the CYP2C mRNAs, CYP2C119 mRNA was expressed less abundantly than CYP2C8, CYP2C9, CYP2C19, and CYP2C76 mRNAs but more abundantly than CYP2C18 mRNA. Recombinant cynomolgus and rhesus CYP2C119 catalyzed progesterone 16α-, 17α-, and 21-hydroxylation and diclofenac and omeprazole oxidations, indicating that CYP2C119 is a functional enzyme. Therefore, the novel CYP2C119 gene, expressed in macaque liver, encodes a functional enzyme that metabolizes human CYP2C substrates and is likely responsible for drug clearances. SIGNIFICANCE STATEMENT: Cytochrome P450 2C119 was found in cynomolgus and rhesus macaques, in addition to the known P450 2C8, 2C9, 2C18, 2C19, 2C76, and 2C93. Cynomolgus and rhesus CYP2C119 contain open reading frames of 489 amino acids with high sequence identity to human CYP2C8. Cynomolgus CYP2C119 mRNA is predominantly expressed in the liver. Recombinant CYP2C119 catalyzed progesterone hydroxylation and diclofenac and omeprazole oxidations. Therefore, the novel CYP2C119 gene expressed in the macaque liver encodes a functional enzyme that metabolizes human CYP2C substrates.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Comprehensive in vitro and in silico assessments of metabolic capabilities of 24 genomic variants of CYP2C19 using two different substrates

Introduction: Most hepatically cleared drugs are metabolized by cytochromes P450 (CYPs), and Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines provide curated clinical references for CYPs to apply individual genome data for optimized drug therapy. However, incorporating novel pharmacogenetic variants into guidelines takes considerable time. Methods: We comprehensively assessed the drug metabolizing capabilities of CYP2C19 variants discovered through population sequencing of two substrates, S-mephenytoin and omeprazole. Results: Based on established functional assays, 75% (18/24) of the variants not yet described in Pharmacogene Variation (PharmVar) had significantly altered drug metabolizing capabilities. Of them, seven variants with inappreciable protein expression were evaluated as protein damaging by all three in silico prediction algorithms, Sorting intolerant from tolerant (SIFT), Polymorphism Phenotyping v2 (PolyPhen-2), and Combined annotation dependent depletion (CADD). The five variants with decreased metabolic capability (<50%) of wild type for either substrates were evaluated as protein damaging by all three in silico prediction algorithms, except CADD exact score of NM_000769.4:c.593T>C that was 19.68 (<20.0). In the crystal structure of the five polymorphic proteins, each altered residue of all those proteins was observed to affect the key structures of drug binding specificity. We also identified polymorphic proteins indicating different tendencies of metabolic capability between the two substrates (5/24). Discussion: Therefore, we propose a methodology that combines in silico prediction algorithms and functional assays on polymorphic CYPs with multiple substrates to evaluate the changes in the metabolism of all possible genomic variants in CYP genes. The approach would reinforce existing guidelines and provide information for prescribing appropriate medicines for individual patients.

A comprehensive analysis of six forms of cytochrome P450 2C (CYP2C) in pigs

Pigs are an important species used in drug metabolism studies; however, the cytochromes P450 (P450s or CYPs) have not been fully investigated in pigs.In this study, pig CYP2C32, CYP2C33, CYP2C34, CYP2C36, CYP2C42, and CYP2C49 cDNAs were isolated and found to contain open reading frames of 490 or 494 amino acids that shared 64-82% sequence identity with human CYP2C8/9/18/19.Pig CYP2C genes formed a gene cluster in a genomic region that corresponded to that of the human CYP2C cluster; an additional gene cluster was formed by pig CYP2C33a and CYP2C33b distant from the first cluster but located in the same chromosome.Among the tissues analysed, these pig CYP2C mRNAs were preferentially expressed in liver, small intestine, and/or kidney; pig CYP2C49, CYP2C32, CYP2C34, and CYP2C33 mRNAs were the most abundant CYP2C mRNAs in liver, jejunum, ileum, and kidney, respectively.Metabolic assays showed that pig CYP2C proteins (heterologously expressed in Escherichia coli) metabolised typical human CYP2C substrates diclofenac, warfarin, and/or omeprazole.The results suggest that these pig CYP2Cs are functional enzymes able to metabolise human CYP2C substrates in liver and small intestine, just as human CYP2Cs do.

Outcomes of a multi-ethnic Asian population on combined treatment with clopidogrel and omeprazole in 12,440 patients

Omeprazole is commonly co-prescribed with clopidogrel. Clopidogrel requires bio-activation by cytochrome P450 CYP2C19. Omeprazole may reduce clopidogrel's antithrombotic efficacy by inhibiting CYP2C19. Studies in Caucasians receiving omeprazole with clopidogrel showed no significant increase in death and myocardial infarction with this drug-drug interaction. There are limited large-scale studies in Asians, who may have a greater prevalence of CYP2C19 loss-of-function polymorphisms. A single centre retrospective cohort study was undertaken based on a review of medication records and prescription data. Patients prescribed clopidogrel from 2009 to 2012 were followed-up with until December 2012 (median:29 months). The primary outcome was all-cause mortality and secondary outcomes were myocardial infarction (MI), cerebrovascular accidents, and subsequent coronary interventions. Of 12,440 patients prescribed clopidogrel, 62%(n = 7714) were on omeprazole (63.8% Chinese, 13.9% Malay, 12.4% Indian, 10.0% others), and 38%(n = 4726) were not on omeprazole or other proton pump inhibitors (62.6% Chinese, 13.5% Malay, 10.7% Indian, 13.2% others). Mortality after co-prescription occurred in 14.3%(n = 1101) of patients, compared to 6.3%(n = 300) of patients prescribed clopidogrel only. Multivariate analysis using propensity score adjusted analysis showed no significant increase in all-cause mortality with co-prescription (adjusted hazards ratio [AHR] 1.13, [95%CI 0.95-1.35]). Patients on co-prescription had a higher risk of subsequent MI (16% vs 3.8%; AHR 2.03 [95%CI 1.70-2.44]), but not of cerebrovascular accidents (5.0% vs 2.0%; AHR 0.98 [95%CI 0.76-1.27]) or coronary interventions (1.7% vs 0.7%; AHR 1.28 [95%CI 0.83-1.96]). The risk of a subsequent MI was higher in the Malay (AHR 2.43 [95%CI 1.68-3.52]) and Chinese (AHR 2.06 [95%CI 1.63-2.60]) population as compared to the Indian (AHR 1.56 [95%CI 1.06-2.31]) population. In conclusion, the use of clopidogrel with omeprazole is associated with an increased risk of MI, but not mortality or stroke, in this multi-ethnic Asian population. These risks appear to vary among different ethnic groups.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Investigation on drug-binding in heme pocket of CYP2C19 with UV-visible and resonance Raman spectroscopies

Cytochrome P450 (CYP) is a class of heme-containing enzymes which mainly catalyze a monooxygenation reaction of various chemicals, and hence CYP plays a key role in the drug metabolism. Although CYP2C19 isoform is a minor hepatic CYP, it metabolizes clinically important drugs such as omeprazole and S‑mephenytoin. In this work, the interaction of purified CYP2C19 WT (CYP2C19) with seven drugs (phenytoin, S‑mephenytoin, omeprazole, lansoprazole, cimetidine, propranolol, and warfarin) was investigated using spectroscopic methods. The binding of each drug and the induced structural change in the heme distal environment were evaluated. Ferric form of CYP2C19 was revealed to contain a six-coordinate low-spin heme with a water molecule as a sixth ligand in a distal site, and the addition of each drug caused varied minor fraction of five-coordinate heme. It was suggested that the ligated water molecule was partly moved away from the heme distal environment and that the degree of water removal was dependent on the type of drugs. The effect on the coordination was varied with the studied drugs with wide variation in the dissociation constants from 2.6 μM for lansoprazole to 5400 μM for warfarin. Phenytoin and S‑mephenytoin showed that binding to CYP2C19 occurred in a stepwise manner and that the coordination of a water molecule was facilitated in the second binding step. In the ferrous CO-bound state, ν(FeCO) stretching mode was clearly observed at 471 cm^-1 in the absence of drugs. The Raman line was greatly up-shifted by omeprazole (487 cm^-1) and lansoprazole (477 cm^-1) but was minimally affected by propranolol, phenytoin, and S‑mephenytoin. These results indicate that slight chemical modification of a drug greatly affects the heme distal environments upon binding.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.

In Vitro Functional Characterisation of Cytochrome P450 (CYP) 2C19 Allelic Variants CYP2C19*23 and CYP2C19*24

Cytochrome P450 (CYP) 2C19 is essential for the metabolism of clinically used drugs including omeprazole, proguanil, and S-mephenytoin. This hepatic enzyme exhibits genetic polymorphism with inter-individual variability in catalytic activity. This study aimed to characterise the functional consequences of CYP2C19*23 (271 G>C, 991 A>G) and CYP2C19*24 (991 A>G, 1004 G>A) in vitro. Mutations in CYP2C19 cDNA were introduced by site-directed mutagenesis, and the CYP2C19 wild type (WT) as well as variants proteins were subsequently expressed using Escherichia coli cells. Catalytic activities of CYP2C19 WT and those of variants were determined by high performance liquid chromatography-based essay employing S-mephenytoin and omeprazole as probe substrates. Results showed that the level of S-mephenytoin 4'-hydroxylation activity of CYP2C19*23 (V _max 111.5 ± 16.0 pmol/min/mg, K _m 158.3 ± 88.0 μM) protein relative to CYP2C19 WT (V _max 101.6 + 12.4 pmol/min/mg, K _m 123.0 ± 19.2 μM) protein had no significant difference. In contrast, the K _m of CYP2C19*24 (270.1 ± 57.2 μM) increased significantly as compared to CYP2C19 WT (123.0 ± 19.2 μM) and V _max of CYP2C19*24 (23.6 ± 2.6 pmol/min/mg) protein was significantly lower than that of the WT protein (101.6 ± 12.4 pmol/min/mg). In vitro intrinsic clearance (CL_int = V _max/K _m) for CYP2C19*23 protein was 85.4 % of that of CYP2C19 WT protein. The corresponding CL_int value for CYP2C19*24 protein reduced to 11.0 % of that of WT protein. These findings suggested that catalytic activity of CYP2C19 was not affected by the corresponding amino acid substitutions in CYP2C19*23 protein; and the reverse was true for CYP2C19*24 protein. When omeprazole was employed as the substrate, K _m of CYP2C19*23 (1911 ± 244.73 μM) was at least 100 times higher than that of CYP2C19 WT (18.37 ± 1.64 μM) and V _max of CYP2C19*23 (3.87 ± 0.74 pmol/min/mg) dropped to 13.4 % of the CYP2C19 WT (28.84 ± 0.61 pmol/min/mg) level. Derived from V _max/K _m, the CL_int value of CYP2C19 WT was 785 folds of CYP2C19*23. K _m and V _max values could not be determined for CYP2C19*24 due to its low catalytic activity towards omeprazole 5'-hydroxylation. Therefore, both CYP2C19*23 and CYP2C19*24 showed marked reduced activities of metabolising omeprazole to 5-hydroxyomeprazole. Hence, carriers of CYP2C19*23 and CYP2C19*24 allele are potentially poor metabolisers of CYP2C19-mediated substrates.

Novel Marmoset Cytochrome P450 2C19 in Livers Efficiently Metabolizes Human P450 2C9 and 2C19 Substrates, S-Warfarin, Tolbutamide, Flurbiprofen, and Omeprazole

The common marmoset (Callithrix jacchus), a small New World monkey, has the potential for use in human drug development due to its evolutionary closeness to humans. Four novel cDNAs, encoding cytochrome P450 (P450) 2C18, 2C19, 2C58, and 2C76, were cloned from marmoset livers to characterize P450 2C molecular properties, including previously reported P450 2C8. The deduced amino acid sequence showed high sequence identities (>86%) with those of human P450 2Cs, except for marmoset P450 2C76, which has a low sequence identity (∼70%) with any human P450 2Cs. Phylogenetic analysis showed that marmoset P450 2Cs were more closely clustered with those of humans and macaques than other species investigated. Quantitative polymerase chain reaction analysis showed that all of the marmoset P450 2C mRNAs were predominantly expressed in liver as opposed to the other tissues tested. Marmoset P450 2C proteins were detected in liver by immunoblotting using antibodies against human P450 2Cs. Among marmoset P450 2Cs heterologously expressed in Escherichia coli, marmoset P450 2C19 efficiently catalyzed human P450 2C substrates, S-warfarin, diclofenac, tolbutamide, flurbiprofen, and omeprazole. Marmoset P450 2C19 had high Vmax and low Km values for S-warfarin 7-hydroxylation that were comparable to those in human liver microsomes, indicating warfarin stereoselectivity similar to findings in humans. Faster in vivo S-warfarin clearance than R-warfarin after intravenous administration of racemic warfarin (0.2 mg/kg) to marmosets was consistent with the in vitro kinetic parameters. These results indicated that marmoset P450 2C enzymes had functional characteristics similar to those of humans, and that P450 2C-dependent metabolic properties are likewise similar between marmosets and humans.

Correlating structure and function of drug-metabolizing enzymes: progress and ongoing challenges

This report summarizes a symposium sponsored by the American Society for Pharmacology and Experimental Therapeutics at Experimental Biology held April 20-24 in Boston, MA. Presentations discussed the status of cytochrome P450 (P450) knowledge, emphasizing advances and challenges in relating structure with function and in applying this information to drug design. First, at least one structure of most major human drug-metabolizing P450 enzymes is known. However, the flexibility of these active sites can limit the predictive value of one structure for other ligands. A second limitation is our coarse-grain understanding of P450 interactions with membranes, other P450 enzymes, NADPH-cytochrome P450 reductase, and cytochrome b5. Recent work has examined differential P450 interactions with reductase in mixed P450 systems and P450:P450 complexes in reconstituted systems and cells, suggesting another level of functional control. In addition, protein nuclear magnetic resonance is a new approach to probe these protein/protein interactions, identifying interacting b5 and P450 surfaces, showing that b5 and reductase binding are mutually exclusive, and demonstrating ligand modulation of CYP17A1/b5 interactions. One desired outcome is the application of such information to control drug metabolism and/or design selective P450 inhibitors. A final presentation highlighted development of a CYP3A4 inhibitor that slows clearance of human immunodeficiency virus drugs otherwise rapidly metabolized by CYP3A4. Although understanding P450 structure/function relationships is an ongoing challenge, translational advances will benefit from continued integration of existing and new biophysical approaches.

Key mutations alter the cytochrome P450 BM3 conformational landscape and remove inherent substrate bias

Cytochrome P450 monooxygenases (P450s) have enormous potential in the production of oxychemicals, due to their unparalleled regio- and stereoselectivity. The Bacillus megaterium P450 BM3 enzyme is a key model system, with several mutants (many distant from the active site) reported to alter substrate selectivity. It has the highest reported monooxygenase activity of the P450 enzymes, and this catalytic efficiency has inspired protein engineering to enable its exploitation for biotechnologically relevant oxidations with structurally diverse substrates. However, a structural rationale is lacking to explain how these mutations have such effects in the absence of direct change to the active site architecture. Here, we provide the first crystal structures of BM3 mutants in complex with a human drug substrate, the proton pump inhibitor omeprazole. Supported by solution data, these structures reveal how mutation alters the conformational landscape and decreases the free energy barrier for transition to the substrate-bound state. Our data point to the importance of such "gatekeeper" mutations in enabling major changes in substrate recognition. We further demonstrate that these mutants catalyze the same 5-hydroxylation reaction as performed by human CYP2C19, the major human omeprazole-metabolizing P450 enzyme.

Structural characterization of human cytochrome P450 2C19: active site differences between P450s 2C8, 2C9, and 2C19

To identify the structural features underlying the distinct substrate and inhibitor profiles of P450 2C19 relative to the closely related human enzymes, P450s 2C8 and 2C9, the atomic structure (Protein Data Bank code 4GQS) of cytochrome P450 2C19 complexed with the inhibitor (2-methyl-1-benzofuran-3-yl)-(4-hydroxy-3,5-dimethylphenyl)methanone (Protein Data Bank chemical component 0XV) was determined to 2.87 Å resolution by x-ray crystallography. The conformation of the peptide backbone of P450 2C19 is most similar to that of P450 2C8, but the substrate-binding cavity of P450 2C8 is much larger than that of P450 2C19 due to differences in the amino acid residues that form the substrate-binding cavities of the two enzymes. In contrast, the substrate-binding cavity of P450 2C19 is much more similar in size to that of the structure of the P450 2C9 flurbiprofen complex than to that of a modified P450 2C9 or that of P450 2C8. The cavities of the P450 2C19 0XV complex and the P450 2C9 flurbiprofen complex differ, however, because the helix B-C loops of the two enzymes are dissimilar. These conformational differences reflect the effects of adjacent structural elements that interact with the B-C loops and that differ between the two enzymes. The availability of a structure for 2C19 will facilitate computational approaches for predictions of substrate and inhibitor binding to this enzyme.

Phylogenomics of the benzoxazinoid biosynthetic pathway of Poaceae: gene duplications and origin of the Bx cluster

Background

The benzoxazinoids 2,4-dihydroxy-1,4-benzoxazin-3-one (DIBOA) and 2,4-dihydroxy-7- methoxy-1,4-benzoxazin-3-one (DIMBOA), are key defense compounds present in major agricultural crops such as maize and wheat. Their biosynthesis involves nine enzymes thought to form a linear pathway leading to the storage of DI(M)BOA as glucoside conjugates. Seven of the genes (Bx1-Bx6 and Bx8) form a cluster at the tip of the short arm of maize chromosome 4 that includes four P450 genes (Bx2-5) belonging to the same CYP71C subfamily. The origin of this cluster is unknown.

Results

We show that the pathway appeared following several duplications of the TSA gene (α-subunit of tryptophan synthase) and of a Bx2-like ancestral CYP71C gene and the recruitment of Bx8 before the radiation of Poaceae. The origins of Bx6 and Bx7 remain unclear. We demonstrate that the Bx2-like CYP71C ancestor was not committed to the benzoxazinoid pathway and that after duplications the Bx2-Bx5 genes were under positive selection on a few sites and underwent functional divergence, leading to the current specific biochemical properties of the enzymes. The absence of synteny between available Poaceae genomes involving the Bx gene regions is in contrast with the conserved synteny in the TSA gene region.

Conclusions

These results demonstrate that rearrangements following duplications of an IGL/TSA gene and of a CYP71C gene probably resulted in the clustering of the new copies (Bx1 and Bx2) at the tip of a chromosome in an ancestor of grasses. Clustering favored cosegregation and tip chromosomal location favored gene rearrangements that allowed the further recruitment of genes to the pathway. These events, a founding event and elongation events, may have been the key to the subsequent evolution of the benzoxazinoid biosynthetic cluster.

Ontogenesis of phase I hepatic drug metabolic enzymes in sheep

Cytochrome P450 (CYP) enzymes are important for the metabolism of many drugs. While there is information on their identity and ontogeny in humans and rodents, similar data in sheep are lacking. In the present study, cDNA sequences of several CYP enzymes (CYP2A6, CYP2C19, CYP2D6) were cloned by rapid amplification of cDNA ends. In adult, newborn and fetal sheep the mRNA and protein levels of these CYPs and the regulatory factor, hepatic nuclear factor 4α (HNF4α) were determined in liver samples using real-time PCR and western blotting. The effect of antenatal glucocorticoid on these enzymes was also studied by i.v. infusion of cortisol (0.45 mg h–1; 80 h) to another group of fetuses. The mRNA and protein levels of the CYPs and HNF4α were low or absent in the fetus, followed by increasing levels in the newborn and adult. Fetal cortisol administration significantly increased the mRNA and protein levels of CYP2D6. Moreover, the correlation observed between the CYP and HNF4α mRNA levels suggests a possible regulatory role for this transcription factor. The findings suggest that fetal and newborn lambs have a low ability to metabolise drugs that are substrates of these enzymes, and that this ability increases with advancing postnatal age, similar to the situation in humans.

Evaluation of the effects of 20 nonsynonymous single nucleotide polymorphisms of CYP2C19 on S-mephenytoin 4'-hydroxylation and omeprazole 5'-hydroxylation

CYP2C19 is a highly polymorphic enzyme that affects the metabolism of a wide range of therapeutic drugs. Almost all the identified alleles of CYP2C19 are derived from nonsynonymous single nucleotide polymorphisms (nsSNPs). The objective of this study was to functionally characterize 20 nsSNPs of CYP2C19, distributed throughout the entire coding region, most of which have not been thoroughly characterized. cDNAs of these variants were constructed and expressed in yeast cells. All variants had similar levels of apoprotein and holoprotein expression, except for CYP2C19.16 and D360N, which had significantly lower holoprotein levels than the wild-type (WT) CYP2C19 enzyme, and CYP2C19.5B, which showed only apoprotein. The activity of the CYP2C19 variants was investigated using two substrates, S-mephenytoin and omeprazole, and six different kinetic parameters were measured. CYP2C19.5B, CYP2C19.6, and CYP2C19.8 were found to be catalytically inactive. The entire dataset of the remaining 17 variants, together with the WT, was analyzed by multivariate analysis. This analysis indicated that CYP2C19.9, CYP2C19.10, CYP2C19.16, CYP2C19.18, CYP2C19.19, A161P, W212C, and D360N were substantially altered in catalytic properties in comparison with the WT, with each of these variants exhibiting either dramatically decreased catalytic activities or higher K(m) values. These results not only generally confirmed the function of previously reported variants but also identified additional reduced-function variants. These findings will greatly extend our understanding of CYP2C19 genetic polymorphisms in humans as well as facilitate the structure-function study of the CYP2C19 protein.

Functional characterization of two novel CYP2C19 variants (CYP2C19*18andCYP2C19*19) found in a Japanese population

Cytochrome P450 2C19 (CYP2C19) plays an important role in the metabolism of a wide range of therapeutic drugs and exhibits genetic polymorphism with interindividual differences in metabolic activity. We have previously described two CYP2C19 allelic variants, namely CYP2C19*18 and CYP2C19*19 with Arg329His/Ile331Val and Ser51Gly/Ile331Val substitutions, respectively. In order to investigate precisely the effect of amino acid substitutions on CYP2C19 function, CYP2C19 proteins of the wild-type (CYP2C19.1B having Ile331Val) and variants (CYP2C19.18 and CYP2C19.19) were heterologously expressed in yeast cells, and their S-mephenytoin 4'-hydroxylation activities were determined. The K(m) value of CYP2C19.19 for S-mephenytoin 4'-hydroxylation was significantly higher (3.0-fold) than that of CYP2C19.1B. Although no significant differences in V(max) values on the basis of microsomal and functional CYP protein levels were observed between CYP2C19.1B and CYP2C19.19, the V(max)/K(m) values of CYP2C19.19 were significantly reduced to 29-47% of CYP2C19.1B. By contrast, the K(m), V(max) or V(max)/K(m) values of CYP2C19.18 were similar to those of CYP2C19.1B. These results suggest that Ser51Gly substitution in CYP2C19.19 decreases the affinity toward S-mephenytoin of CYP2C19 enzyme, and imply that the genetic polymorphism of CYP2C19*19 also causes variations in the clinical response to drugs metabolized by CYP2C19.

Few alterations in clinical pathology and histopathology observed in a CYP2C18&19 humanized mice model

Background

This study was performed to characterize a gene-addition transgenic mouse containing a BAC (bacterial artificial chromosome) clone spanning the human CYP2C18&19 genes (tg-CYP2C18&19).

Methods

Hemizygous tg-CYP2C18&19, 11 week old mice were compared with wild-type littermates to obtain information regarding clinical status, clinical pathology and anatomical pathology. After one week of clinical observations, blood samples were collected, organs weighed, and tissues collected for histopathology.

Results

In males, the tissue weights were lower in tg-CYP2C18&19 than in wild-type mice for brain (p ≤ 0.05), adrenal glands (p ≤ 0.05) and brown fat deposits (p ≤ 0.001) while the heart weight was higher (p ≤ 0.001). In female tg-CYP2C18&19, the tissue weights were lower for brain (p ≤ 0.001) and spleen (p ≤ 0.001) compared to wild-type females. Male tg-CYP2C18&19 had increased blood glucose levels (p ≤ 0.01) while females had decreased blood triglyceride levels (p ≤ 0.01).

Conclusion

Despite the observed alterations, tg-CYP2C18&19 did not show any macroscopic or microscopic pathology at the examined age. Hence, these hemizygous transgenic mice were considered to be viable and healthy animals.

Lipophilicity Relationships in Inhibitors of CYP2C9 and CYP2C19 Enzymes

Quantitative structure-activity relationships (QSARs) within a series of cytochrome P450 2C9 (CYP2C9) and cytochrome P450 2C19 (CYP2C19) inhibitors are reported. In particular, it is noted that compound lipophilicity, in the form of log P values (where P is the octanol/water partition coefficient), is an important factor in explaining the variation in inhibitory potency within these series of compounds, many of which also act as substrates for the respective enzymes. In addition, there is a role for hydrogen bonding and pi-pi stacking interactions within the P450 active site which represent secondary factors in the binding processes of these compounds.

In silico description of differential enantioselectivity in methoxychlor O-demethylation by CYP2C enzymes

Methoxychlor undergoes metabolism by cytochrome P450 (CYP) enzymes forming a chiral mono-phenolic derivative (Mono-OH-M) as main metabolite. In the current study, members of the CYP2C family were examined for their chiral preference in Mono-OH-M formation. CYP2C9 and CYP2C19 possessed high enantioselectivity favoring the formation of S-Mono-OH-M; CYP2C3 showed no enantioselectivity, whereas CYP2C5 slightly favored the formation of R-Mono-OH-M. Molecular modeling calculations were utilized in order to explain the observed differences in chiral preference of CYP2C enzymes. Molecular docking calculations could describe neither the existence of chiral preference in metabolism, nor the enantiomer which is preferentially formed. Molecular dynamic calculations were also carried out and were found to be useful for accurate description of chiral preference in biotransformation of methoxychlor by CYP2C enzymes. An in silico model capable of predicting chiral preference in cytochrome P450 enzymes in general can be developed based on the analysis of the stability and rigidity parameters of interacting partners during molecular dynamic simulation.

Important amino acid residues that confer CYP2C19 selective activity to CYP2C9

Although CYP2C9 and CYP2C19 display 91% sequence identity at the amino acid level, the two enzymes have distinct substrate specificities for compounds such as diclofenac, progesterone and (S)-mephenytoin. Amino acid substitutions in CYP2C9 were made based on an alignment of CYP2C9, CYP2C19 and monkey CYP2C43 sequences. Mutants of CYP2C9 were expressed in Escherichia coli. Sixteen amino acids, which are common to both CYP2C19 and CYP2C43 but different between CYP2C9 and CYP2C19, were substituted in CYP2C9 (CYP2C9-16aa). Next, the mutated amino acids in CYP2C9-16aa were individually reverted to those of CYP2C9 to examine the effect of each substitution on the enzymatic activity for CYP2C marker substrates. In addition, the role of the F-G loop in CYP2C9 and CYP2C19 was examined for substrate specificity and enzymatic activity. Our results showed: (i) CYP2C9-16aa displays 11% (S)-mephenytoin 4'-hydroxylase and full omeprazole 5-hydroxylase activity compared with that of CYP2C19; (ii) residue 286 is important for conferring CYP2C9-like enzyme activity on CYP2C9-16aa and residue 442 in CYP2C19 may be involved in the interaction with NADPH-P450 reductase; (iii) substitution of the F-G loop in CYP2C9 to that of CYP2C19 enhances tolbutamide p-methyhydroxylase and diclofenac 4'-hydroxylase activities and confers partial (S)-mephenytoin 4'-hydroxylase and omeprazole 5-hydroxylase activities, which are attributed to CYP2C19.

Engineering human cytochrome P450 enzymes into catalytically self-sufficient chimeras using molecular Lego

The membrane-bound human cytochrome P450s have essential roles in the metabolism of endogenous compounds and drugs. Presented here are the results on the construction and characterization of three fusion proteins containing the N-terminally modified human cytochrome P450s CYP2C9, CY2C19 and CYP3A4 fused to the soluble NADPH-dependent oxidoreductase domain of CYP102A1 from Bacillus megaterium. The constructs, CYP2C9/BMR, CYP2C19/BMR and CYP3A4/BMR are well expressed in Escherichia coli as holo proteins. The chimeras can be purified in the absence of detergent and the purified enzymes are both active and correctly folded in the absence of detergent, as demonstrated by circular dichroism and functional studies. Additionally, in comparison with the parent P450 enzyme, these chimeras have greatly improved solubility properties. The chimeras are catalytically self-sufficient and present turnover rates similar to those reported for the native enzymes in reconstituted systems, unlike previously reported mammalian cytochrome P450 fusion proteins. Furthermore the specific activities of these chimeras are not dependent on the enzyme concentration present in the reaction buffer and they do not require the addition of accessory proteins, detergents or phospholipids to be fully active. The solubility, catalytic self-sufficiency and wild-type like activities of these chimeras would greatly simplify the studies of cytochrome P450 mediated drug metabolism in solution.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Functional Characterization of Novel Allelic Variants of CYP2C9 Recently Discovered in Southeast Asians

CYP2C9 was recently resequenced in 150 Asian subjects from Singapore. Several new coding variants were reported, and these variants are now named CYP2C9*14 (R125H), CYP2C9*15 (S162X), CYP2C9*16 (T299A), CYP2C9*17 (P382S), CYP2C9*18 (D397A), and CYP2C9*19 (Q454H). The CYP2C9*18 variant also contained an I359L change previously associated with the CYP2C9*3 allele. In this study, we assessed the functional consequences of the new coding changes. cDNAs containing each of the new coding changes were constructed by site-directed mutagenesis and expressed in a bacterial cDNA expression system, the allelic proteins were partially purified, and their ability to hydroxylate a prototype CYP2C9 substrate was assayed. Expression of cDNAs in Escherichia coli containing either the D397A change or the S162X (premature stop codon) could not be detected either spectrally or at the apoprotein level. CYP2C9.14 and CYP2C9.16 exhibited 80 to 90% lower catalytic activity toward tolbutamide at two substrate concentrations compared with wild-type CYP2C9.1. Kinetic analysis confirmed that CYP2C9.14 and CYP2C9.16 have a higher Km and a >90% lower intrinsic clearance of tolbutamide compared with wild-type CYP2C9.1. Both CYP2C9.17 and CYP2C9.19 proteins exhibited modest 30 to 40% decreases in catalytic activity toward tolbutamide. Thus, CYP2C9*15 and CYP2C9*18 may represent null alleles, whereas CYP2C9*14 and CYP2C9*16 allelic variants produce proteins that are clearly catalytically defective in vitro, indicating the existence of new defective putative alleles of CYP2C9 in Asians.

Analysis of human cytochrome P450 2C8 substrate specificity using a substrate pharmacophore and site-directed mutants

The structural determinants of substrate specificity of human liver cytochrome P450 2C8 (CYP2C8) were investigated using site-directed mutants chosen on the basis of a preliminary substrate pharmacophore and a three-dimensional (3D) model. Analysis of the structural features common to CYP2C8 substrates exhibiting a micromolar K(m) led to a substrate pharmacophore in which the site of oxidation by CYP2C8 is 12.9, 8.6, 4.4, and 3.9 A from features that could establish ionic or hydrogen bonds, and hydrophobic interactions with protein amino acid residues. Comparison of this pharmacophore with a 3D model of CYP2C8 constructed using the X-ray structure of CYP2C5 suggested potential CYP2C8 amino acid residues that could be involved in substrate recognition. Twenty CYP2C8 site-directed mutants were constructed and expressed in yeast to compare their catalytic activities using five CYP2C8 substrates that exhibit different structures and sizes [paclitaxel, fluvastatin, retinoic acid, a sulfaphenazole derivative (DMZ), and diclofenac]. Mutation of arginine 241 had marked effects on the hydroxylation of anionic substrates of CYP2C8 such as retinoic acid and fluvastatin. Serine 100 appears to be involved in hydrogen bonding interactions with a polar site of the CYP2C8 substrate pharmacophore, as shown by the 3-4-fold increase in the K(m) of paclitaxel and DMZ hydroxylation after the S100A mutation. Residues 114, 201, and 205 are predicted to be in close contact with substrates, and their mutations lead either to favorable hydrophobic interactions or to steric clashes with substrates. For instance, the S114F mutant was unable to catalyze the 6alpha-hydroxylation of paclitaxel. The S114F and F205A mutants were the best catalysts for retinoic acid and paclitaxel (or fluvastatin) hydroxylation, respectively, with k(cat)/K(m) values 5 and 2.1 (or 2.4) times higher, respectively, than those found for CYP2C8. Preliminary experiments of docking of the substrate into the experimentally determined X-ray structure of substrate-free CYP2C8, which became available quite recently [Schoch, G. A., et al. (2004) J. Biol. Chem. 279, 9497], were consistent with key roles for S100, S114, and F205 residues in substrate binding. The results suggest that the effects of mutation of arginine 241 on anionic substrate hydroxylation could be indirect and result from alterations of the packing of helix G with helix B'.

Charge and Substituent Effects on Affinity and Metabolism of Benzbromarone-Based CYP2C19 Inhibitors

Human cytochrome P450 (CYP) 2C19 is one of the most important CYP2C family members responsible for metabolizing commonly prescribed drugs. This research describes synthetic modifications to benzbromarone (Bzbr) to create the most potent CYP2C19 inhibitor ever reported. The most important features enabling analogues of Bzbr to bind to CYP2C19 with high affinity are low acidity (high pK(a) or nonionizability) and hydrophobic substituents adjacent to the phenol moiety. Though CYP2C19 was known to prefer neutral substrates, the extent was perhaps not realized until the anionic, parent compound Bzbr (K(i) = 3.7 microM) was compared to a less acidic dimethyl analogue (K(i) = 0.033 microM). However, differences in affinity for anionic and neutral Bzbr analogues did not appear to affect the regiospecificity of their metabolism by CYP's 2C19 and 2C9. In addition, some Bzbr analogues were metabolized both on the phenol and benzofuran rings. By using a substrate with a methyl ether in place of the Bzbr phenol, it was shown that some Bzbr analogues must be able to freely reposition after binding and oxidize the more energetically favorable position. Normally, O-demethylation of this methyl ether is favored over benzofuran hydroxylation based on ion current from LC/MS. Deuterium substitution of the methyl ether results in an inverse isotope effect on benzofuran hydroxylation (i.e. increased oxidation of this less favorable site). Likewise, Bzbr-based CoMFA models of CYP2C19 demonstrated no clear preference for any one ligand alignment. This suggests results from this modeling method must be interpreted carefully for each CYP isoform. In summary, Bzbr analogues have demonstrated they can be adapted to other CYP2C enzymes in order to probe isoform-specific properties.

Discovery of new potentially defective alleles of human CYP2C9

CYP2C9 is a clinically important enzyme, responsible for the metabolism of numerous clinically important therapeutic drugs. In the present study, we discovered 38 single nucleotide polymorphisms in CYP2C9 by resequencing of genomic DNA from 92 individuals from three different racial groups. Haplotype analysis predicted that there are at least 21 alleles of CYP2C9 in this group of individuals. Six new alleles were identified that contained coding changes: L19I (CYP2C9*7), R150H (CYP2C9*8), H251R (CYP2C9*9), E272G (CYP2C9*10), R335W(CYP2C9*11) and P489S (CYP2C9*12). When expressed in a bacterial cDNA expression system, several alleles exhibited altered catalytic activity. CYP2C9*11 appeared to be a putative poor metabolizer allele, exhibiting a three-fold increase in the Km and more than a two-fold decrease in the intrinsic clearance for tolbutamide. Examination of the crystal structure of human CYP2C9 reveals that R335 is located in the turn between the J and J' helices and forms a hydrogen-bonding ion pair with D341 from the J' helix. Abolishing this interaction in CYP2C9*11 individuals could destabilize the secondary structure and alter the substrate affinity. This new putative poor metabolizer (PM) allele was found in Africans. A second potentially PM allele CYP2C9*12 found in a racially unidentified sample also exhibited a modest decrease in the Vmax and the intrinsic clearance for tolbutamide in a recombinant system. Further clinical studies are needed to determine the effect of these new polymorphisms on the metabolism of CYP2C9 substrates.

The structure of human cytochrome P450 2C9 complexed with flurbiprofen at 2.0-A resolution

The structure of human P450 2C9 complexed with flurbiprofen was determined to 2.0 A by x-ray crystallography. In contrast to other structurally characterized P450 2C enzymes, 2C5, 2C8, and a 2C9 chimera, the native catalytic domain of P450 2C9 differs significantly in the conformation of the helix F to helix G region and exhibits an extra turn at the N terminus of helix A. In addition, a distinct conformation of the helix B to helix C region allows Arg-108 to hydrogen bond with Asp-293 and Asn-289 on helix I and to interact directly with the carboxylate of flurbiprofen. These interactions position the substrate for regioselective oxidation in a relatively large active site cavity and are likely to account for the high catalytic efficiency exhibited by P450 2C9 for the regioselective oxidation of several anionic non-steroidal anti-inflammatory drugs. The structure provides a basis for interpretation of a number of observations regarding the substrate selectivity of P450 2C9 and the observed effects of mutations on catalysis.

Characterization of the CYP2C8 Active Site by Homology Modeling

To compare the features of the active sites of CYP2C8, CYP2C9, and CYP2C19, homology modeling was performed based on the crystallographic coordinates of mammalian CYP2C5. It was found that CYP2C8 has a much larger pocket than the other forms due to the existence of an additional pocket. The approach to the additional pocket is comprised of Ile102, Ser114, Leu208, Val366, and Ile476, and the side chains of Ser114, Val366, and Ile476, which are smaller than the corresponding residues in the other CYPs, enable access to the pocket. The general features of the active site in the CYP2C8 model are similar to those of the previously constructed CYP3A4 model, which may account for the 2 CYPs sharing some of their substrates. The CYP2C8 model was validated by examining the bound orientation of paclitaxel and showing that it is consistent with the formation of the 6-beta hydroxylated derivative during metabolism. Docked paclitaxel was found to form a hydrogen bond with the side chain of Asn 99, which is a characteristic residue of CYP2C8 and is located in the additional pocket. Descriptors for CYP2C8 and CYP2C9 substrates were also examined with the molecular operating environment (MOE). The descriptor by which CYP2C8 and CYP2C9 substrates were classified most distinctly was found to be molar refractivity, which might be related to the longer shape and more polar nature of the active site of CYP2C8 in the CYP2C subfamily.

Crystal structure of human cytochrome P450 2C9 with bound warfarin

Cytochrome P450 proteins (CYP450s) are membrane-associated haem proteins that metabolize physiologically important compounds in many species of microorganisms, plants and animals. Mammalian CYP450s recognize and metabolize diverse xenobiotics such as drug molecules, environmental compounds and pollutants. Human CYP450 proteins CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4 are the major drug-metabolizing isoforms, and contribute to the oxidative metabolism of more than 90% of the drugs in current clinical use. Polymorphic variants have also been reported for some CYP450 isoforms, which has implications for the efficacy of drugs in individuals, and for the co-administration of drugs. The molecular basis of drug recognition by human CYP450s, however, has remained elusive. Here we describe the crystal structure of a human CYP450, CYP2C9, both unliganded and in complex with the anti-coagulant drug warfarin. The structure defines unanticipated interactions between CYP2C9 and warfarin, and reveals a new binding pocket. The binding mode of warfarin suggests that CYP2C9 may undergo an allosteric mechanism during its function. The newly discovered binding pocket also suggests that CYP2C9 may simultaneously accommodate multiple ligands during its biological function, and provides a possible molecular basis for understanding complex drug-drug interactions.

A rational approach to Re-engineer cytochrome P450 2B1 regioselectivity based on the crystal structure of cytochrome P450 2C5

The regioselectivity for progesterone hydroxylation by cytochrome P450 2B1 was re-engineered based on the x-ray crystal structure of cytochrome P450 2C5. 2B1 is a high K(m) progesterone 16alpha-hydroxylase, whereas 2C5 is a low K(m) progesterone 21-hydroxylase. Initially, nine individual 2B1 active-site residues were changed to the corresponding 2C5 residues, and the mutants were purified from an Escherichia coli expression system and assayed for progesterone hydroxylation. At 150 microm progesterone, I114A, F297G, and V363L showed 5-15% of the 21-hydroxylase activity of 2C5, whereas F206V showed high activity for an unknown product and a 13-fold decrease in K(m). Therefore, a quadruple mutant, I114A/F206V/F297G/V363L (Q), was constructed that showed 60% of 2C5 progesterone 21-hydroxylase activity and 57% regioselectivity. Based on their 2C5-like testosterone hydroxylation profiles, S294D and I477F alone and in combination were added to the quadruple mutant. All three mutants showed enhanced regioselectivity (70%) for progesterone 21-hydroxylation, whereas only Q/I477F had a higher k(cat). Finally, the remaining three single mutants, V103I, V367L, and G478V, were added to Q/I477F and Q/S294D/I477F, yielding seven additional multiple mutants. Among these, Q/V103I/S294D/I477F showed the highest k(cat) (3-fold higher than that of 2C5) and 80% regioselectivity for progesterone 21-hydroxylation. Docking of progesterone into a three-dimensional model of this mutant indicated that 21-hydroxylation is favored. In conclusion, a systematic approach to convert P450 regioselectivity across subfamilies suggests that active-site residues are mainly responsible for regioselectivity differences between 2B1 and 2C5 and validates the reliability of 2B1 models based on the crystal structure of 2C5.

Mutation of tyrosine 190 to alanine eliminates the inactivation of cytochrome P450 2B1 by peroxynitrite

We have previously reported that cytochrome P450 2B1 was inactivated by peroxynitrite and that the decrease in the catalytic activity correlated with an increase in the nitration of tyrosine. Digestion of the peroxynitrite-treated P450 2B1 with Lys C followed by amino acid sequencing of the major nitrotyrosine-containing peptide demonstrated that it spanned residues 160-225. This peptide contains two tyrosine residues at positions 190 and 203. In this study, we mutated Tyr 190 to Ala (Y190A) and Tyr 203 to Ala (Y203A) in wild-type recombinant P450 2B1 (WT) in order to identify the specific residue(s) that is nitrated and to determine whether nitrotyrosine formation is reponsible for the peroxynitrite-mediated inactivation of P450 2B1. All three P450s were expressed in Escherichia coli, purified to homogeneity, and characterized. The catalytic activities for four different substrates of P450 2B1 increased approximately 2-fold for the Y203A mutant, but decreased by about 60% for the Y190A mutant when compared to WT. The addition of peroxynitrite to the P450s resulted in concentration-dependent decreases in the catalytic activities of WT and Y203A, but no loss of the catalytic activities of Y190A. The extent of tyrosine nitration of Y190A by peroxynitrite decreased by approximately 75% as compared with WT or the Y203A protein. Following digestion of the peroxynitrite-modified proteins with Lys C, a major nitrotyrosine-containing peptide was detected from WT and Y203A, but not from Y190A. Collectively, these results indicate that Tyr 190 is the target residue for peroxynitrite-mediated nitration and that nitration of this tyrosine is a responsible for the inactivation of P450 2B1. Modeling studies suggest that Tyr 190 may play a structural role in maintaining the integrity of the protein for maximal activity through hydrogen bonding with Glu 149.

Substrate selectivity of human cytochrome P450 2C9: importance of residues 476, 365, and 114 in recognition of diclofenac and sulfaphenazole and in mechanism-based inactivation by tienilic acid

A series of six site-directed mutants of CYP 2C9 were constructed with the aim to better define the amino acid residues that play a critical role in substrate selectivity of CYP 2C9, particularly in three distinctive properties of this enzyme: (i) its selective mechanism-based inactivation by tienilic acid (TA), (ii) its high affinity and hydroxylation regioselectivity toward diclofenac, and (iii) its high affinity for the competitive inhibitor sulfaphenazole (SPA). The S365A mutant exhibited kinetic characteristics for the 5-hydroxylation of TA very similar to those of CYP 2C9; however, this mutant did not undergo any detectable mechanism-based inactivation by TA, which indicates that the OH group of Ser 365 could be the nucleophile forming a covalent bond with an electrophilic metabolite of TA in TA-dependent inactivation of CYP 2C9. The F114I mutant was inactive toward the hydroxylation of diclofenac; moreover, detailed analyses of its interaction with a series of SPA derivatives by difference visible spectroscopy showed that the high affinity of SPA to CYP 2C9 (K(s)=0.4 microM) was completely lost when the phenyl substituent of Phe 114 was replaced with the alkyl group of Ile (K(s)=190+/-20 microM), or when the phenyl substituent of SPA was replaced with a cyclohexyl group (K(s)=120+/-30 microM). However, this cyclohexyl derivative of SPA interacted well with the F114I mutant (K(s)=1.6+/-0.5 microM). At the opposite end, the F94L and F110I mutants showed properties very similar to those of CYP 2C9 toward TA and diclofenac. Finally, the F476I mutant exhibited at least three main differences compared to CYP 2C9: (i) big changes in the k(cat) and K(m) values for TA and diclofenac hydroxylation, (ii) a 37-fold increase of the K(i) value found for the inhibition of CYP 2C9 by SPA, and (iii) a great change in the regioselectivity of diclofenac hydroxylation, the 5-hydroxylation of this substrate by CYP 2C9 F476I exhibiting a k(cat) of 28min(-1). These data indicate that Phe 114 plays an important role in recognition of aromatic substrates of CYP 2C9, presumably via Pi-stacking interactions. They also provide the first experimental evidence showing that Phe 476 plays a crucial role in substrate recognition and hydroxylation by CYP 2C9.

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.

Identification and functional characterization of new potentially defective alleles of human CYP2C19

CYP2C19 is a clinically important enzyme responsible for the metabolism of a number of therapeutic drugs, such as mephenytoin, omeprazole, diazepam, proguanil, propranolol and certain antidepressants. Genetic polymorphisms in this enzyme result in poor metabolizers of these drugs. There are racial differences in the incidence of the poor metabolizer trait, which represents 13-23% of Asians but only 3-5% of Caucasians. In this study, single nucleotide polymorphisms (SNPs) in CYP2C19 were identified by direct sequencing of genomic DNA from 92 individuals from three different racial groups of varied ethnic background, including Caucasians, Asians and blacks. Several new alleles were identified containing the coding changes Arg114 His (CYP2C19*9), Pro227 Leu (CYP2C19*10), Arg150 His (CYP2C19*11), stop491 Cys (CYP2C19*12), Arg410 Cys (CYP2C19*13), Leu17 Pro (CYP2C19*14) and Ile19 Leu (CYP2C19*15). When expressed in a bacterial cDNA expression system, CYP2C19*9 exhibited a modest decrease in the V(max) for 4'-hydroxylation of -mephenytoin, and no alteration in its affinity for reductase. CYP2C19*10 exhibited a dramatically higher K(m) and lower V(max) for mephenytoin. CYP2C19*12was unstable and expressed poorly in a bacterial cDNA expression system. Clinical studies will be required to confirm whether this allele is defective in vivo. CYP2C19*9, CYP2C19*10 and CYP2C19*12 all occurred in African-Americans, or individuals of African descent, and represent new potentially defective alleles of CYP2C19 which are predicted to alter risk of these populations to clinically important drugs.

Substrate routes to the buried active site may vary among cytochromes P450: mutagenesis of the F-G region in P450 2B1

Until recently, all known structures of bacterial cytochromes P450 suggested that substrate access to the buried active site occurred via the F-G region, a surface loop distal to the heme cavity. However, the structure of P450 51 indicates a large opening from the protein surface along the I helix N-terminus, at right angles to the F-G channel. The single available microsomal P450 structure (2C5) does not obviously favor one potential access route over the other. To determine whether the F-G region forms part of the substrate access channel in the microsomal cytochrome P450 2B1, 11 residues between positions 208 and 230 were substituted with smaller and larger side chains in a highly expressed truncated form of the enzyme. Steady-state kinetic parameters were determined with the substrates testosterone, 7-ethoxy-4-trifluoromethylcoumarin (7-EFC), and 7-benzyoxyresorufin (7-BR). The largest changes, 2-6-fold increases in k(cat) with testosterone and 7-EFC, were observed for L209A, which also exhibits an altered testosterone metabolite profile and probably forms part of the active site roof. F219W demonstrated little or no activity with any of the three substrates examined, although the K(s) value for benzphetamine binding was unaltered. S221F showed little activity with 7-BR. No significant changes were observed in K(m)(testosterone) or S(50)(7-EFC) values for any of the mutants, in stark contrast to the 10-fold and 100-fold changes in K(m) observed for mutants in this region of other cytochromes P450. The minimal changes in 2B1 do not support access via the F-G region of 2B1 and suggest the alternate access route identified in P450 51.

Modelling human cytochromes P450 involved in drug metabolism from the CYP2C5 crystallographic template

A historical background to homology modelling of human P450s involved in drug metabolism is outlined, showing that the progress in crystallographic studies of bacterial forms of enzyme and, latterly, determination of a mammalian P450 crystal structure, has enabled the production of increasingly satisfactory models of human P450 enzymes. The methodology for the generation of P450 models by homology with crystallographic template structures is summarized, and recent results of CYP2C5-constructed models of P450s are described. These indicate that selective substrates are able to fit within the putative active sites of each enzyme, where key contacts with complementary amino acid residues are largely consistent with the results of site-directed mutagenesis experiments and metabolic studies. Consequently, the CYP2C5 crystal structure can be regarded at the current paradigm for homology modelling of the drug metabolizing P450s, especially those from the CYP2 family.

Amino acid residues affecting the activities of human cytochrome P450 2C9 and 2C19

The amino acid residues affecting the substrate specificity of human cytochrome P450 CYP2C9 and CYP2C19 for their metabolic activities were investigated using chimeras and mutant enzymes, which were constructed by replacing the corresponding residues. Although CYP2C19 showed nearly the same tolbutamide hydroxylase activity as CYP2C9, the activities for the CYP2C19(H99I) mutant and the chimeras that replaced residues 1-212 were much lower than those for CYP2C19. The activities of the CYP2C19(H99I) mutant and the chimeras that replaced residues 228-340 were lower than those for CYP2C19 toward S-mephenytoin, aminopyrine, and testosterone. These results suggest that residues in substrate recognition site (SRS) 3 and 4 are important for the substrate specificity, whereas His99 is important in the substrate binding of CYP2C19. For the 4'-hydroxylation of diclofenac, CYP2C9 had a lower K(m) and a higher V(max) than CYP2C19. Although the V(max) values for the CYP2C9(1-288)/CYP2C19(289-490) chimera and the CYP2C9(I99H, V292A, F295L, I331V) mutant were comparable to those for CYP2C9, its K(m) value was comparable to that for CYP2C19. The V(max) and K(m) values for the CYP2C19(1-288)/CYP2C9(289-490) chimera were comparable to those for CYP2C19, and the activity by CYP2C9(1-230)/CYP2C19(231-490) chimera was negligible. These results suggest that the residues 292, 295, and/or 331 of CYP2C9 are essential for the recognition of substrate in CYP2C9 and that the residues of 231-288 including SRS 3 are important for the metabolizing capacity of CYP2C9.

Development of a combined protein and pharmacophore model for cytochrome P450 2C9

A combined protein and pharmacophore model for cytochrome P450 2C9 (CYP2C9) has been derived using various computational chemistry techniques. A combination of pharmacophore modeling (using 31 metabolic pathways for 27 substrates), protein modeling (using the rabbit CYP2C5/3 crystal structure), and molecular orbital calculations was used to derive a model that incorporated steric, electronic, and chemical stability properties. The initial pharmacophore model (based on a subset of 17 metabolic pathways for 16 substrates) and the protein model used to construct the combined model were derived independently and showed a large degree of complementarity. The combined model is in agreement with experimental results concerning the substrates used to derive the model and with site-directed mutagenesis data available for CYP2C9. The model has been successfully used to predict the metabolism of substrates not used to construct the model, of which four examples are discussed in detail. The model has also been successful in explaining the differences in substrate specificity between CYP2C9 and CYP2C19.

Summary of information on human CYP enzymes: human P450 metabolism data

This chapter is an update of the data on substrates, reactions, inducers, and inhibitors of human CYP enzymes published previously by Rendic and DiCarlo (1), now covering selection of the literature through 2001 in the reference section. The data are presented in a tabular form (Table 1) to provide a framework for predicting and interpreting the new P450 metabolic data. The data are formatted in an Excel format as most suitable for off-line searching and management of the Web-database. The data are presented as stated by the author(s) and in the case when several references are cited the data are presented according to the latest published information. The searchable database is available either as an Excel file (for information contact the author), or as a Web-searchable database (Human P450 Metabolism Database, www.gentest.com) enabling the readers easy and quick approach to the latest updates on human CYP metabolic reactions.

Analysis of selective regions in the active sites of human cytochromes P450, 2C8, 2C9, 2C18, and 2C19 homology models using GRID/CPCA

This study demonstrates a selectivity analysis using the GRID/CPCA strategy on four human cytochrome P450 2C homology models (CYP2C8, 2C9, 2C18, and 2C19). Although the four enzymes share more than 80% amino acid sequence identity, the substrate specificity differs. To investigate the selectivity of the enzymes and the amino acids that determine the specificity of each CYP2C enzyme, a selectivity analysis was made using GRID/CPCA. In the GRID calculations 10 probes were used covering hydrophobic, steric, and hydrogen bond acceptor and donor interactions. The selectivity analysis showed that the most important determinants of selectivity among the CYP2C models are the geometrical features of the active sites and the hydrophobic interactions. The selectivity analysis singled out CYP2C8 as the most different of the four CYP2C enzymes with amino acids with distinct properties in positions 114, 205, and 476 (Ser, Phe, and Ile, respectively) compared to the other enzymes. An inverse pharmacophore model for CYP2C9 was constructed from the selective regions, and the model agreed with the docking of diclofenac where the properties of the ligand overlapped with the pharmacophoric points in the model.

Common and uncommon cytochrome P450 reactions related to metabolism and chemical toxicity

Cytochrome P450 (P450) enzymes catalyze a variety of reactions and convert chemicals to potentially reactive products as well as make compounds less toxic. Most of the P450 reactions are oxidations. The majority of these can be rationalized in the context of an FeO(3+) intermediate and odd electron abstraction/rebound mechanisms; however, other iron-oxygen complexes are possible and alternate chemistries can be considered. Another issue regarding P450-catalyzed reactions is the delineation of rate-limiting steps in the catalytic cycle and the contribution to reaction selectivity. In addition to the rather classical oxidations, P450s also catalyze less generally discussed reactions including reduction, desaturation, ester cleavage, ring expansion, ring formation, aldehyde scission, dehydration, ipso attack, one-electron oxidation, coupling reactions, rearrangement of fatty acid and prostaglandin hydroperoxides, and phospholipase activity. Most of these reactions are rationalized in the context of high-valent iron-oxygen intermediates and Fe(2+) reductions, but others are not and may involve acid-base catalysis. Some of these transformations are involved in the bioactivation and detoxication of xenobiotic chemicals.

Importance of amino acid residue 474 for substrate specificity of canine and human cytochrome p450 3A enzymes

Canine cytochromes P450 3A12 and 3A26 are identical in sequence at 481 of 503 amino-acid positions but exhibit different substrate specificities. A recent study utilizing chimeric enzymes and site-directed mutagenesis identified three residues (187, 368, and 369) that contribute to differences in steroid hydroxylation and also indicated the presence of additional determinants of specificity among the 44 carboxyl terminal residues. Therefore, three 3A26 multiple mutants (I187T-S368P-V369I-S467P, I187T-S368P-V369I-S474P, and I187T-S368P-V369I-R476K-I477L-T479A-R480Q) were constructed. Insertion of 3A12 residue Pro-474 into 3A26 I187T-S368P-V369I resulted in metabolite profiles with testosterone, androstenedione, and progesterone very similar to 3A12. Substitution of Pro-474 with Ser in P450 3A12 or human 3A4 significantly increased 2beta-hydroxylase activity with all three steroids. Residue 474 was also found to be an important contributor to diazepam metabolism by the canine and human enzymes. The results provide further evidence for the role of steric constraints exerted by the enzyme in P450 3A-mediated oxidations.