Identification of amino acid substitutions that confer a high affinity for sulfaphenazole binding and a high catalytic efficiency for warfarin metabolism to P450 2C19

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.

The Impact of Inorganic Systems and Photoactive Metal Compounds on Cytochrome P450 Enzymes and Metabolism: From Induction to Inhibition

While cytochrome P450 (CYP; P450) enzymes are commonly associated with the metabolism of organic xenobiotics and drugs or the biosynthesis of organic signaling molecules, they are also impacted by a variety of inorganic species. Metallic nanoparticles, clusters, ions, and complexes can alter CYP expression, modify enzyme interactions with reductase partners, and serve as direct inhibitors. This commonly overlooked topic is reviewed here, with an emphasis on understanding the structural and physiochemical basis for these interactions. Intriguingly, while both organometallic and coordination compounds can act as potent CYP inhibitors, there is little evidence for the metabolism of inorganic compounds by CYPs, suggesting a potential alternative approach to evading issues associated with rapid modification and elimination of medically useful compounds.

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.

Novel Cytochrome P450 2C94 Functionally Metabolizes Diclofenac and Omeprazole in Dogs

Cytochromes P450 (P450s or CYPs) are important drug-metabolizing enzymes. Because dogs are frequently used in drug metabolism studies, knowledge of dog CYP2C enzymes is essential because in humans these enzymes are abundant and play major roles in liver and intestine. The present study identified and characterized novel dog CYP2C94 along with previously identified dog CYP2C21 and CYP2C41. Dog CYP2C21, CYP2C41, and CYP2C94 cDNAs, respectively, contained open reading frames of 490, 489, and 496 amino acids and shared high-sequence identities (70%, 75%, and 58%) with human CYP2Cs. Dog CYP2C94 mRNA was preferentially expressed in liver, just as dog CYP2C21 and CYP2C41 mRNAs were. In dog liver, CYP2C21 mRNA was the most abundant, followed by CYP2C94 and CYP2C41 mRNAs. Moreover, the hepatic expressions of all three dog CYP2C mRNAs varied in four individual dogs, two of which did not express CYP2C41 mRNA. The three dog CYP2C genes had similar gene structures, and CYP2C94, although located on the same chromosome, was in a genomic region far from the gene cluster containing CYP2C21 and CYP2C41 Metabolic assays with recombinant proteins showed that dog CYP2C94, along with CYP2C21 and CYP2C41, efficiently catalyzed oxidations of diclofenac, warfarin, and/or omeprazole, indicating that dog CYP2C94 is a functional enzyme. Novel dog CYP2C94 is expressed abundantly in liver and encodes a functional enzyme that metabolizes human CYP2C substrates; it is, therefore, likely responsible for drug clearances in dogs. SIGNIFICANCE STATEMENT: Novel dog cytochrome P450 2C94 (CYP2C94) was identified and characterized along with dog CYP2C21 and CYP2C41. Dog CYP2C94, isolated from liver, had 58% sequence identity and a close phylogenetic relationship with its human homologs and was expressed in liver at the mRNA level. Dog CYP2C94 (and CYP2C21 and CYP2C41) catalyzed oxidations of diclofenac and omeprazole, human CYP2C9 and CYP2C19 substrates, respectively, but CYP2C41 also hydroxylated warfarin. CYP2C94 is therefore a functional drug-metabolizing enzyme likely responsible for drug clearances in dogs.

Construction of a fused grid-based CYP2C19-Template system and the application

A ligand-accessible space in the CYP2C19 active site was reconstituted as a fused grid-based Template with the use of structural data of the ligands. An evaluation system of CYP2C19-mediated metabolism has been developed on Template with the introduction of the idea of Trigger-residue initiated ligand-movement and fastening. Reciprocal comparison of the data of simulation on Template with experimental results suggested a unified way of the interaction of CYP2C19 and its ligands through the simultaneous plural-contact with Rear-wall of Template. CYP2C19 was expected to have a room for ligands between vertically standing parallel walls termed Facial-wall and Rear-wall, which were separated by a distance corresponding to 1.5-Ring (grid) diameter size. The ligand sittings were stabilized through contacts with Facial-wall and the left-side borders of Template including specific Position 29 or Left-end after Trigger-residue initiated ligand-movement. Trigger-residue movement is suggested to force ligands to stay firmly in the active site and then to initiate CYP2C19 reactions. Simulation experiments for over 450 reactions of CYP2C19 ligands supported the system established.

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.

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.

The X-Ray Crystal Structure of the Human Mono-Oxygenase Cytochrome P450 3A5-Ritonavir Complex Reveals Active Site Differences between P450s 3A4 and 3A5

The contributions of cytochrome P450 3A5 to the metabolic clearance of marketed drugs is unclear, but its probable role is to augment the metabolism of several drugs that are largely cleared by P450 3A4. Selective metabolism by 3A4 is often a concern in drug development owing to potential drug-drug interactions and the variability of 3A4 and 3A5 expression. The contribution of P450 3A5 to these clearance pathways varies between individuals owing to genetic differences and similarities and differences in the metabolic properties of 3A5 compared with 3A4. To better understand the structural differences between P450s 3A4 and 3A5, the structure of 3A5 complexed with ritonavir was determined by X-ray crystallography to a limiting resolution of 2.91 Å. The secondary and tertiary structures of 3A5 and 3A4 are similar, but the architectures of their active sites differ. The 3A5 active site is taller and narrower than that of 3A4. As a result, ritonavir adopts a distinctly different conformation to fit into the cavity of 3A5 than seen for 3A4. These structural changes reflect amino acid differences that alter the conformation of the helix F through helix G region in the upper portion of the cavity and ionic interactions between residues in the beta-sheet domain that reduce the width of the cavity. The structural differences exhibited by 3A4 and 3A5 suggest that the overlap of catalytic activities may reflect molecular flexibility that determines how alternative conformers fit into the different active site architectures of the two enzymes.

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.

Electrochemically driven drug metabolism via a CYP1A2-UGT1A10 bienzyme confined in a graphene nano-cage

A graphene nano-cage with regulatable space for the assembly of a cytochrome P450 1A2-UDP-glucuronosyltransferase 1A10 bienzyme complex has been constructed via a click reaction, and successfully used to study drug sequential metabolism using an electrochemically-driven method.

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.

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.

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.

Re-engineering of CYP2C9 to probe acid-base substrate selectivity

A common feature of many CYP2C9 ligands is their weak acidity. As revealed by crystallography, the structural basis for this behavior involves a charge-pairing interaction between an anionic moiety on the substrate and an active site R108 residue. In the present study we attempted to re-engineer CYP2C9 to better accept basic ligands by charge reversal at this key residue. We expressed and purified the R108E and R108E/D293N mutants and compared their ability with that of native CYP2C9 to interact with (S)-warfarin, diclofenac, pyrene, propranolol, and ibuprofen amine. As expected, the R108E mutant maintained all the native enzyme's pyrene 1-hydroxylation activity, but catalytic activity toward diclofenac and (S)-warfarin was abrogated. In contrast, the double mutant displayed much less selectivity in its behavior toward these control ligands. Neither of the mutants displayed significant enhancement of propranolol metabolism, and all three preparations exhibited a type II (inhibitor) rather than type I (substrate) spectrum with ibuprofen amine, although binding became progressively weaker with the single and double mutants. Collectively, these data underscore the importance of the amino acid at position 108 in the acid substrate selectivity of CYP2C9, highlight the accommodating nature of the CYP2C9 active site, and provide a cautionary note regarding facile re-engineering of these complex cytochrome P450 active sites.

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.

High-affinity binding of [3H]cimetidine to a heme-containing protein in rat brain

[(3)H]Cimetidine (3HCIM) specifically binds to an unidentified site in the rat brain. Because recently described ligands for this site have pharmacological activity, 3HCIM binding was characterized. 3HCIM binding was saturable, heat-labile, and distinct from the histamine H(2) receptor. To test the hypothesis that 3HCIM binds to a cytochrome P450 (P450), the effects of nonselective and isoform-selective P450 inhibitors were studied. The heme inhibitor KCN and the nonselective P450 inhibitor metyrapone both produced complete, concentration-dependent inhibition of 3HCIM binding (K(i) = 1.3 mM and 11.9 muM, respectively). Binding was largely unaffected by inhibitors of CYP1A2, 2B6, 2C8, 2C9, 2D6, 2E1, and 19A1 but was eliminated by inhibitors of CYP2C19 (tranylcypromine) and CYP3A4 (ketoconazole). Synthesis and testing of CC11 [4(5)-(benzylthiomethyl)-1H-imidazole] and CC12 [4(5)-((4-iodobenzyl)-thiomethyl)-1H-imidazole] confirmed both drugs to be high-affinity inhibitors of 3HCIM binding. On recombinant human P450s, CC12 was a potent inhibitor of CYP2B6 (IC(50) = 11.7 nM), CYP2C19 (51.4 nM), and CYP19A1 (140.7 nM) and had a range of activities (100-494 nM) on nine other isoforms. Although the 3HCIM binding site pharmacologically resembles some P450s, eight recombinant human P450s and three recombinant rat P450s did not exhibit 3HCIM binding. Inhibition by KCN and metyrapone suggests that 3HCIM binds to a heme-containing brain protein (possibly a P450). However, results with selective P450 inhibitors, recombinant P450 isoforms, and a P450 antibody did not identify a 3HCIM-binding P450 isoform. Finally, CC12 is a new, potent inhibitor of CYP2B6 and CYP2C19 that may be a valuable tool for P450 research.

Clinical and toxicological relevance of CYP2C9: drug-drug interactions and pharmacogenetics

CYP2C9 is a major cytochrome P450 enzyme that is involved in the metabolic clearance of a wide variety of therapeutic agents, including nonsteroidal antiinflammatories, oral anticoagulants, and oral hypoglycemics. Disruption of CYP2C9 activity by metabolic inhibition or pharmacogenetic variability underlies many of the adverse drug reactions that are associated with the enzyme. CYP2C9 is also the first human P450 to be crystallized, and the structural basis for its substrate and inhibitor selectivity is becoming increasingly clear. New, ultrapotent inhibitors of CYP2C9 have been synthesised that aid in the development of quantitative structure-activity relationship (QSAR) models to facilitate drug redesign, and extensive resequencing of the gene and studies of its regulation will undoubtedly help us understand interindividual variability in drug response and toxicity controlled by this enzyme.

Progress in cytochrome P450 active site modeling

Models capable of predicting the possible involvement of cytochromes P450 in the metabolism of drugs or drug candidates are important tools in drug discovery and development. Ideally, functional information would be obtained from crystal structures of all the cytochromes P450 of interest. Initially, only crystal structures of distantly related bacterial cytochromes P450 were available-comparative modeling techniques were used to bridge the gap and produce structural models of human cytochromes P450, and thereby obtain some useful functional information. A significant step forward in the reliability of these models came four years ago with the first crystal structure of a mammalian cytochrome P450, rabbit CYP2C5, followed by the structures of two human enzymes, CYP2C8 and CYP2C9, and a second rabbit enzyme, CYP2B4. The evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism, is presented as a case study.

The Structure of Human Microsomal Cytochrome P450 3A4 Determined by X-ray Crystallography to 2.05-Å Resolution

The structure of P450 3A4 was determined by x-ray crystallography to 2.05-A resolution. P450 3A4 catalyzes the metabolic clearance of a large number of clinically used drugs, and a number of adverse drug-drug interactions reflect the inhibition or induction of the enzyme. P450 3A4 exhibits a relatively large substrate-binding cavity that is consistent with its capacity to oxidize bulky substrates such as cyclosporin, statins, taxanes, and macrolide antibiotics. Family 3A P450s also exhibit unusual kinetic characteristics that suggest simultaneous occupancy by smaller substrates. Although the active site volume is similar to that of P450 2C8 (PDB code: 1PQ2), the shape of the active site cavity differs considerably due to differences in the folding and packing of portions of the protein that form the cavity. Compared with P450 2C8, the active site cavity of 3A4 is much larger near the heme iron. The lower constraints on the motions of small substrates near the site of oxygen activation may diminish the efficiency of substrate oxidation, which may, in turn, be improved by space restrictions imposed by the presence of a second substrate molecule. The structure of P450 3A4 should facilitate a better understanding of the substrate selectivity of the enzyme.

Active-site characteristics of CYP2C19 and CYP2C9 probed with hydantoin and barbiturate inhibitors

Three series of N-3 alkyl substituted phenytoin, nirvanol, and barbiturate derivatives were synthesized and their inhibitor potencies were tested against recombinant CYP2C19 and CYP2C9 to probe the interaction of these ligands with the active sites of these enzymes. All compounds were found to be competitive inhibitors of both enzymes, although the degree of inhibitory potency was generally much greater towards CYP2C19. Inhibitor stereochemistry did not markedly influence K(i) towards CYP2C9, and log P adequately predicted inhibitor potency for this enzyme. In contrast, stereochemistry was an important factor in determining inhibitor potency towards CYP2C19. (S)-(+)-N-3-Benzylnirvanol and (R)-(-)-N-3-benzylphenobarbital emerged as the most potent and selective CYP2C19 inhibitors, with K(i) values of < 250nM--at least two orders of magnitude greater inhibitor potency than towards CYP2C9. Both inhibitors were metabolized preferentially at their C-5 phenyl substituents, indicating that CYP2C19 prefers to orient the N-3 substituents away from the active oxygen species. These features were incorporated into expanded CoMFA models for CYP2C9, and a new, validated CoMFA model for CYP2C19.

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.

Mutagenesis and molecular dynamics suggest structural and functional roles for residues in the N-terminal portion of the cytochrome P450 2B1 I helix

To investigate their potential roles in ligand access, binding, and subsequent metabolism, residues in the N-terminal portion of the cytochrome P450 2B1 I helix were mutated to alanine and phenylalanine. Of the 18 mutants from E286 to S294 only 7 yielded holoprotein in an Escherichia coli expression system. Substitutions at positions 289, 290, 292, and 294 caused >/= 2-fold changes in kcat and/or Km for two or more of the 2B1 substrates examined, testosterone, 7-ethoxy-4-trifluoromethylcoumarin, 7-benzyloxyresorufin, and benzphetamine. I290 substitutions had the largest effects on steady-state parameters for three substrates and increased benzphetamine affinity. Steered molecular dynamics simulations of testosterone egress along the I helix identified hydrophobic interactions with I290, L293, and S294 and water bridges to E286 and S294. Sensitivity of holoprotein formation to substitution and effects on substrate binding and metabolism suggest structural and functional roles for residues in the N-terminus of the cytochrome P450 2B1 I helix.

Differential roles of Arg97, Asp293, and Arg108 in enzyme stability and substrate specificity of CYP2C9

CYP2C9 metabolizes a wide range of drugs, many of which are negatively charged at physiological pH. Therefore, it has been thought that complementarily charged amino acid(s) are critically involved in substrate binding. Previous studies have implicated arginine residues at positions 97, 105, and 108 and aspartate at position 293 in the normal catalytic function of the enzyme. To elucidate the role of these amino acids in the substrate specificity of CYP2C9, a series of mutants were constructed and analyzed for functional activity, thermal stability, and ligand binding. Charge-modifying mutations at positions 97, 105, and 293 decreased catalytic activity toward diclofenac, (S)-warfarin, and pyrene in a substrate-independent manner with Arg105 the least, and Arg97 the most, sensitive amino acids in this regard. Decreases in functional activity paralleled thermal instability of the mutants, suggesting that loss of function reflects more generalized structural changes rather than the absence of a specifically charged amino acid at these three positions. The R108H mutant was inactive toward all three substrates because of unexpected nitrogen ligation to the heme. Conversely, the R108F mutant exhibited substrate-dependent catalytic behavior, with almost complete loss of activity toward (S)-warfarin and diclofenac, but preservation of pyrene metabolism. In addition, the R108F mutation abrogated the Type I difference spectra induced by flurbiprofen and benzbromarone, obligate anions at physiological pH. These data identify critical roles for Arg97 and Asp293 in the structural stability of the enzyme and demonstrate a selective role for Arg108 in the binding and metabolism of negatively charged substrates of CYP2C9.

Structure-function relationships of rat liver CYP3A9 to its human liver orthologs: site-directed active site mutagenesis to a progesterone dihydroxylase

CYP3A9 is an estrogen-inducible ortholog of human liver CYP3A4 with 76.5% sequence identity to CYP3A4. Unlike CYP3A4, it is a very poor testosterone 6beta- and 2beta-hydroxylase, but a relatively better catalyst of progesterone monohydroxylation largely at 6beta, 16alpha, and 21 positions with negligible 6beta, 21-dihydroxylation. We reasoned that such differences in substrate catalyses must be due to differences in the active site architecture of each CYP3A enzyme. Indeed, alignment of CYP3A4 substrate recognition sites (SRSs) with the corresponding regions of CYP3A9 sequence revealed that of the 22 fully divergent residues, 4 reside in SRS regions [P107N (SRS-1), M371G (SRS-5), and L479K and G480Q (SRS-6)]. Accordingly, we substituted these and other divergent CYP3A9 SRS residues with the corresponding residues of CYP3A4 and/or CYP3A5. Our findings of the influence of these site-directed mutations of the CYP3A9 active site on its catalysis of testosterone and three other established but structurally different CYP3A substrates (progesterone, imipramine, and carbamazepine) are described. These findings revealed that some mutations (N107P, N107S, V207T, G371M, and Q480G) not only improved the ability of CYP3A9 to hydroxylate testosterone at the 6beta and 2beta positions, but also converted it into a robust progesterone 6beta, 21-dihydroxylase. The latter in the case of CYP3A9N107P was accompanied by a shift from sigmoidal to hyperbolic enzyme-substrate kinetics. In contrast, the catalytic potential of CYP3A9 mutants K206N, K206S, M240V, and K479L/Q480G was either relatively unchanged or negligible to nonexistent. Together these findings attest to the unique substrate-active site fit of each CYP3A enzyme.

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.

Quantitative liquid chromatography/mass spectrometry/mass spectrometry warfarin assay for in vitro cytochrome P450 studies

A sensitive assay using high-performance liquid chromatography tandem mass spectrometry (MS/MS) has been established for the quantitative analysis of cytochrome P450 form-specific activities using warfarin as a probe substrate. Four metabolites, 6-, 7-, 8-, and 10-hydroxywarfarin, were chromatographically resolved within 10 min using gradient mobile phases. The mass spectrometry was operated under negative ionization mode. The MS/MS product ion spectra of warfarin and the metabolites were generated using collision-activated dissociation and interpreted. The abundant product ions of the metabolites were selected for quantification applying multiple reaction monitoring. Quantification was based on a quadratic or power curve of the peak area ratio of the metabolite over the internal standard against the respective concentration of the metabolite. This assay has been validated from 2 to 1000 nM for 10-hydroxywarfarin and from 2 to 5000 nM for 6-, 7-, and 8-hydroxywarfarin and successfully applied to evaluate cytochrome P450-mediated drug-drug interactions in vitro using human hepatocytes and liver microsomal preparations.

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.

Competitive CYP2C9 inhibitors: enzyme inhibition studies, protein homology modeling, and three-dimensional quantitative structure-activity relationship analysis

This study describes the generation of a three-dimensional quantitative structure activity relationship (3D-QSAR) model for 29 structurally diverse, competitive CYP2C9 inhibitors defined experimentally from an initial data set of 73 compounds. In parallel, a homology model for CYP2C9 using the rabbit CYP2C5 coordinates was built. For molecules with a known interaction mode with CYP2C9, this homology model, in combination with the docking program GOLD, was used to select conformers to use in the 3D-QSAR analysis. The remaining molecules were docked, and the GRID interaction energies for all conformers proposed by GOLD were calculated. This was followed by a principal component analysis (PCA) of the GRID energies for all conformers of all compounds. Based on the similarity in the PCA plot to the inhibitors with a known interaction mode, the conformer to be used in the 3D-QSAR analysis was selected. The compounds were randomly divided into two groups, the training data set (n = 21) to build the model and the external validation set (n = 8). The PLS (partial least-squares) analysis of the interaction energies against the K(i) values generated a model with r(2) = 0.947 and a cross-validation of q(2) = 0.730. The model was able to predict the entire external data set within 0.5 log units of the experimental K(i) values. The amino acids in the active site showed complementary features to the grid interaction energies in the 3D-QSAR model and were also in agreement with mutagenesis studies.

Protein engineering of cytochromes P-450

The cytochromes P-450 are an immensely important superfamily of heme-containing enzymes. They catalyze the monooxygenation of an enormous range of substrates. In bacteria, cytochromes P-450 are known to catalyze the hydroxylation of environmentally significant substrates such as camphor, phenolic compounds and many herbicides. In eukaryotes, these enzymes perform key roles in the synthesis and interconversion of steroids, while in mammals hepatic cytochromes P-450 are vital for the detoxification of many drugs. As such, the cytochromes P-450 are of considerable interest in medicine and biotechnology and are obvious targets for protein engineering. The purpose of this article is to illustrate the ways in which protein engineering has been used to investigate and modify the properties of cytochromes P-450. Illustrative examples include: the manipulation of substrate selectivity and regiospecificity, the alteration of membrane binding properties, and probing the route of electron transfer.

Arginines 97 and 108 in CYP2C9 are important determinants of the catalytic function

Human cytochrome P450 2C9 (CYP2C9) is one of the major drug metabolising enzymes which exhibits a broad substrate specificity. The B-C loop is located in the active-site but has been difficult to model, owing to its diverse and flexible structure. To elucidate the function of the B-C loop we used homology modelling based on the Cyp102 structure in combination with functional studies of mutants using diclofenac as a model substrate for CYP2C9. The study shows the importance of the conserved arginine in position 97 and the arginine in position 108 for the catalytic function. The R97A mutant had a 13-fold higher K(m) value while the V(max) was in the same order as the wild type. The R108 mutant had a 100-fold lower activity with diclofenac compared to the wild-type enzyme. The other six mutants (S95A, F100A, L102A, E104A, R105A, and N107A) had kinetic parameters similar to the CYP2C9 wild-type. Our homology model based on the CYP102 structure as template indicates that R97, L102, and R105 are directed into the active site, whereas R108 is not. The change in catalytic function when arginine 97 was replaced with alanine and the orientation of this amino acid in our homology model indicates its importance for substrate interaction.

Mechanism-based inactivation of cytochrome P450 2B1 by 7-ethynylcoumarin: verification of apo-P450 adduction by electrospray ion trap mass spectrometry

7-Ethynylcoumarin was synthesized as a potential mechanism-based inhibitor, and it was found to be an effective inactivator of 7-ethoxy-4-(trifluoromethyl)coumarin (7EFC) O-deethylation catalyzed by purified, reconstituted P450 2B1. In contrast, 7-ethynylcoumarin demonstrated minimal inactivation of P450 2A6-mediated 7-hydroxycoumarin formation. The inactivation of P450 2B1 demonstrated pseudo-first-order kinetics and was NADPH- and inhibitor-dependent. The maximal rate constant for the inactivation of 2B1 was 0.39 min(-)(1) at 30 degrees C, and thus, the time required to inactivate 50% of the P450 2B1 that was present (t(1/2)) was 1.8 min. The estimated concentration which led to half-maximal inactivation (K(I)) was 25 microM. No protection from inactivation was seen in the presence of nucleophiles (glutathione and sodium cyanide), an iron chelator (deferroxamine), or superoxide dismutase and catalase. Addition of the substrate (7EFC) protected P450 2B1 from inactivation, in a concentration-dependent manner. The partition ratio for P450 2B1 was 25; i.e., the number of metabolic events was 25-fold higher than the number of inactivating events. Incubations of 7-ethynylcoumarin with P450 2B1 for 10 min resulted in an 80% loss in enzymatic activity, while 90% of the ability to form a reduced-CO complex remained. This activity loss was not recovered following dialysis, indicative of irreversible inactivation. Covalent attachment of the entire inhibitor and oxygen to apo-P450 2B1, in a 1:1 ratio, was shown via electrospray ion trap mass spectrometry. This method also verified the absence of modification to the heme or the cytochrome P450 reductase. Taken together, the characterization of the inhibition seen with P450 2B1 and 7-ethynylcoumarin was consistent with all of the criteria required to distinguish a mechanism-based inactivator. In addition, electrospray ion trap mass spectrometry has the potential to be applied to protein adducts above and beyond those associated with the mechanism-based inactivation of cytochrome P450s.

Homology modeling and substrate binding study of human CYP2C9 enzyme

The CYP2C subfamily of human liver P450 isozymes is of major importance in drug metabolism. The most abundant 2C isozyme, CYP2C9, regioselectively hydroxylates a wide variety of substrates. A major obstacle to understanding this specificity in human CYP2C9 is the absence of a 3D structure. A 3D model of CYP2C9 was built, assessed, and used to characterize explicit enzyme-substrate complexes using methods previously developed in our laboratory. The 3D model was assessed by determining its stability to unconstrained molecular dynamics and by comparison of specific properties with those of known protein structures. The CYP2C9 model was then used to characterize explicit enzyme complexes with three structurally and chemically diverse substrates: (S)-naproxen, phenytoin, and progesterone. Each substrate was found to bind to the enzyme with a favorable interaction energy and to remain in the binding site during unconstrained molecular dynamics. Moreover, the mode of binding of each substrate led to calculated preferred hydroxylation sites consistent with experiment. Binding-site residues identified for the models included Arg 105 and Arg97 as key cationic residues, as well as Asn 202, Asp 293, Pro 101, Leu 102, Gly 296, and Phe 476. Site-specific mutations are proposed for further integrated computational and experimental study.