Structure-function analysis of human cytochrome P450 3A4 using 7-alkoxycoumarins as active-site probes

Systematic Optimization of HPO-CPR to Boost (+)-Nootkatone Synthesis in Engineered Saccharomyces cerevisiae

As an important and expensive natural sesquiterpene compound in grapefruit, the interest in (+)-nootkatone is stimulated by its strong grapefruit-like odor and physiological activities, which induce efforts for its microbial production. However, the low catalytic efficiency of the cytochrome P450-P450 reductase (HPO-CPR) system is the main challenge. We developed a high-throughput screening (HTS) method using the principle of the color reaction between carbonyl compounds and 2,4-dinitrophenylhydrazine (DNPH), which could rapidly screen the activity of candidate HPO mutants. After optimizing the pairing of HPO and CPR and through semirational design, the optimal mutant HPO_M18 with catalytic performance 2.54 times that of the initial was obtained. An encouraging (+)-nootkatone titer of 2.39 g/L was achieved through two-stage fed-batch fermentation after metabolic engineering and endoplasmic reticulum engineering, representing the highest titer reported to date. Our findings lay the foundation for the development of an economically viable bioprocess for (+)-nootkatone.

CYP3A4 mutation causes vitamin D-dependent rickets type 3

Genetic forms of vitamin D-dependent rickets (VDDRs) are due to mutations impairing activation of vitamin D or decreasing vitamin D receptor responsiveness. Here we describe two unrelated patients with early-onset rickets, reduced serum levels of the vitamin D metabolites 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D, and deficient responsiveness to parent and activated forms of vitamin D. Neither patient had a mutation in any genes known to cause VDDR; however, using whole exome sequencing analysis, we identified a recurrent de novo missense mutation, c.902T>C (p.I301T), in CYP3A4 in both subjects that alters the conformation of substrate recognition site 4 (SRS-4). In vitro, the mutant CYP3A4 oxidized 1,25-dihydroxyvitamin D with 10-fold greater activity than WT CYP3A4 and 2-fold greater activity than CYP24A1, the principal inactivator of vitamin D metabolites. As CYP3A4 mutations have not previously been linked to rickets, these findings provide insight into vitamin D metabolism and demonstrate that accelerated inactivation of vitamin D metabolites represents a mechanism for vitamin D deficiency.

Biochemical Mechanisms for Geographical Adaptations to Novel Toxin Exposures in Butterflyfish

Some species of butterflyfish have had preyed upon corals for millions of years, yet the mechanism of butterflyfish specialized coral feeding strategy remains poorly understood. Certain butterflyfish have the ability to feed on allelochemically rich soft corals, e.g. Sinularia maxima. Cytochrome P450 (CYP) is the predominant enzyme system responsible for the detoxification of dietary allelochemicals. CYP2-like and CYP3A-like content have been associated with butterflyfish that preferentially consumes allelochemically rich soft corals. To investigate the role of butterflyfish CYP2 and CYP3A enzymes in dietary preference, we conducted oral feeding experiments using homogenates of S. maxima and a toxin isolated from the coral in four species of butterflyfish with different feeding strategies. After oral exposure to the S. maxima toxin 5-episinulaptolide (5ESL), which is not normally encountered in the Hawaiian butterflyfish diet, an endemic specialist, Chaetodon multicinctus experienced 100% mortality compared to a generalist, Chaetodon auriga, which had significantly more (3–6 fold higher) CYP3A-like basal content and catalytic activity. The specialist, Chaetodon unimaculatus, which preferentially feed on S. maxima in Guam, but not in Hawaii, had 100% survival, a significant induction of 8–12 fold CYP3A-like content, and an increased ability (2-fold) to metabolize 5ESL over other species. Computer modeling data of CYP3A4 with 5ESL were consistent with microsomal transformation of 5ESL to a C15-16 epoxide from livers of C. unimaculatus. Epoxide formation correlated with CYP3A-like content, catalytic activity, induction, and NADPH-dependent metabolism of 5ESL. These results suggest a potentially important role for the CYP3A family in butterflyfish-coral diet selection through allelochemical detoxification.

Coumarin Derivatives as Substrate Probes of Mammalian Cytochromes P450 2B4 and 2B6: Assessing the Importance of 7-Alkoxy Chain Length, Halogen Substitution, and Non-Active Site Mutations

Using a combined structural and biochemical approach, the functional importance of a recently described peripheral pocket bounded by the E-, F-, G-, and I-helices in CYP2B4 and 2B6 was probed. Three series of 4-substituted-7-alkoxycoumarin derivatives with -H, -CH3, or -CF3 at the 4 position of the coumarin core were used initially to monitor functional differences between CYP2B4 and 2B6. 7-Ethoxy-4-(trifluoromethyl)coumarin (7-EFC) displayed the highest catalytic efficiency among these substrates. Mutants were made to alter side-chain polarity (V/E194Q) or bulk (F/Y244W) to alter access to the peripheral pocket. Modest increases in catalytic efficiency of 7-EFC O-deethylation by the mutants were magnified considerably by chlorination or bromination of the substrate ethoxy chain. A structure of CYP2B6 Y244W in complex with (+)-α-pinene was solved at 2.2 Å and showed no CYMAL-5 in the peripheral pocket. A ligand free structure of CYP2B4 F244W was solved at 3.0 Å with CYMAL-5 in the peripheral pocket. In both instances, comparison of the respective wild-type and mutant CYP2B enzymes revealed that CYMAL-5 occupancy of the peripheral pocket had little effect on the topology of active site residue side-chains, despite the fact that the peripheral pocket and active site are located on opposite sides of the I-helix. Analysis of available CYP2B structures suggest that the effect of the amino acid substitutions within the peripheral pocket derive from altered interactions between the F and G helices.

Structure-function analysis of porcine cytochrome P450 3A29 in the hydroxylation of T-2 toxin as revealed by docking and mutagenesis studies

T-2 toxin, one of the type A trichothecenes, presents a potential hazard to human and animal health. Our previous work demonstrated that porcine cytochrome P450 3A29 (CYP3A29) played an important role in the hydroxylation of T-2 toxin. To identify amino acids involved in this metabolic process, T-2 toxin was docked into a homology model of CYP3A29 based on a crystal structure of CYP3A4 using AutoDock 4.0. Nine residues of CYP3A29, Arg105, Arg106, Phe108, Ser119, Lys212, Phe213, Phe215, Arg372 and Glu374, which were found within 5 Å around T-2 toxin were subjected to site-directed mutagenesis. In the oxidation of nifedipine, the CL_int value of R106A was increased by nearly two-folds compared with the wild-type CYP3A29, while the substrate affinities and CL_int values of S119A and K212A were significantly reduced. In the hydroxylation of T-2 toxin, the generation of 3′-OH-T-2 by R105A, S119A and K212A was significantly less than that by the wild-type, whereas R106A slightly increased the generation of 3′-OH-T-2. These results were further confirmed by isothermal titration calorimetry analysis, suggesting that these four residues are important in the hydroxylation of T-2 toxin and Arg105 may be a specific recognition site for the toxin. Our study suggests a possible structure-function relationship of CYP3A29 in the hydroxylation of T-2 toxin, providing with new insights into the mechanism of CYP3A enzymes in the biotransformation of T-2 toxin.

Characterization of the active site properties of CYP4F12

Cytochrome P450 4F12 is a drug-metabolizing enzyme that is primarily expressed in the liver, kidney, colon, small intestine, and heart. The properties of CYP4F12 that may impart an increased catalytic selectivity (decreased promiscuity) were explored through in vitro metabolite elucidation, kinetic isotope effect experiments, and computational modeling of the CYP4F12 active site. By using astemizole as a probe substrate for CYP4F12 and CYP3A4, it was observed that although CYP4F12 favored astemizole O-demethylation as the primary route of metabolism, CYP3A4 was capable of metabolizing astemizole at multiple sites on the molecule. Deuteration of astemizole at the site of O-demethylation resulted in an isotope effect of 7.1 as well as an 8.3-fold decrease in the rate of clearance for astemizole by CYP4F12. Conversely, although an isotope effect of 3.8 was observed for the formation of the O-desmethyl metabolite when deuterated astemizole was metabolized by CYP3A4, there was no decrease in the clearance of astemizole. Development of a homology model of CYP4F12 based on the crystal structure of cytochrome P450 BM3 predicted an active site volume for CYP4F12 that was approximately 76% of the active site volume of CYP3A4. As predicted, multiple favorable binding orientations were available for astemizole docked into the active site of CYP3A4, but only a single binding orientation with the site of O-demethylation oriented toward the heme was identified for CYP4F12. Overall, it appears that although CYP4F12 may be capable of binding similar ligands to other cytochrome P450 enzymes such as CYP3A4, the ability to achieve catalytically favorable orientations may be inherently more difficult because of the increased steric constraints of the CYP4F12 active site.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

Cytochrome P450 3A5 plays a prominent role in the oxidative metabolism of the anti-human immunodeficiency virus drug maraviroc

Maraviroc is an anti-human immunodeficiency virus drug that acts by blocking viral entry into target cells. With use of ultra-performance liquid chromatography-mass spectrometry several monooxygenated, dioxygenated, and glucuronidated metabolites of maraviroc were identified both in vitro and in vivo. Characterization of the enzymes involved in the production of these metabolites determined that cytochrome P450 3A5 was the principal enzyme responsible for the formation of an abundant metabolite of maraviroc that resulted from oxygenation of the dichlorocyclohexane ring. For the formation of this metabolite, the V(max) values for CYP3A4 and CYP3A5 were 0.04 and 0.93 pmol · min⁻¹ · pmol P450⁻¹, and the K(m) values were 11.1 and 48.9 μM, respectively. Furthermore, human liver microsomes isolated from donors homozygous for the loss-of-function CYP3A5*3 allele exhibited a 79% decrease in formation of this metabolite compared with those homozygous for the wild-type CYP3A5*1 allele. To probe which divergent residues between CYP3A4 and CYP3A5 might play a role in the differential activities of these enzymes toward maraviroc, mutations were introduced into both enzymes and metabolism of maraviroc was measured. A CYP3A5 L57F mutant exhibited a 61% decrease in the formation of this metabolite, whereas formation by a CYP3A4 F57L mutant was increased by 337% compared with that of the wild type. Taken together, these data provide novel insights into the biotransformation of maraviroc as well as the potential role of CYP3A4 and CYP3A5 divergent residues in the enzymatic activities of these two highly homologous enzymes.

Activation of CYP3A-mediated testosterone 6β-hydroxylation by tanshinone IIA and midazolam 1-hydroxylation by cryptotanshinone in human liver microsomes

This study evaluated the in vitro activation of CYP3A-mediated midazolam 1-hydroxylation and testosterone 6β-hydroxylation by tanshinone I, tanshinone IIA, and cryptotanshinone. The abilities of tanshinones to activate CYP3A-mediated midazolam 1-hydroxylation and testosterone 6β-hydroxylation in human liver microsomes (HLMs) were tested. Substrate- and effector-dependent activation of CYP3A by tanshinones were both observed. Cryptotanshinone was shown to activate CYP3A-mediated midazolam 1-hydroxylation in a concentration-dependent manner. In contrast, tanshinone IIA and tanshinone I did not activate this hydroxylation reaction. In addition, tanshinone IIA activated CYP3A-mediated testosterone 6β-hydroxylation, whereas cryptotanshinone and tanshinone I did not. The results from our study enhance the understanding of CYP3A activation by tanshinone IIA and cryptotanshinone in HLMs. Additionally, these data allow for an accurate prediction of the magnitude and likelihood of Danshen-drug interactions.

Homology modelling of CYP3A4 from the CYP2C5 crystallographic template: analysis of typical CYP3A4 substrate interactions

1. The results of homology modelling of cytochrome P4503A4 (CYP3A4), which is a human enzyme of major importance for the Phase 1 metabolism of drug substrates, from the CYP2C5 crystal structure is reported. 2. The overall homology between the two protein sequences was generally good (46%) with 24% of amino acid residues being identical and a 22% similarity between matched pairs in the CYP3A4 and CYP2C5 aligned sequences, thus indicating that CYP2C5 represents a viable template for modelling CYP3A4 by homology. 3. The CYP3A4 model appears to show consistency with the reported findings from the extensive site-directed mutagenesis studies already published. 4. Typical CYP3A4 substrates, such as midazolam, testosterone, nifedipine and verapamil, are shown to fit the putative active site of the enzyme structure in a manner consistent with their known positions of metabolism.

CYP3 phylogenomics: evidence for positive selection of CYP3A4 and CYP3A7

OBJECTIVE

CYP3A metabolizes 50% of currently prescribed drugs and is frequently involved in clinically relevant drug interactions. The understanding of roles and regulations of the individual CYP3A genes in pharmacology and physiology is incomplete.

METHODS

Using genomic sequences from 16 species we investigated the evolution of CYP3 genomic loci over a period of 450 million years.

RESULTS

CYP3A genes in amniota evolved from two ancestral CYP3A genes. Upon the emergence of eutherian mammals, one of them was lost, whereas, the other acquired a novel genomic environment owing to translocation. In primates, CYP3A underwent rapid evolutionary changes involving multiple gene duplications, deletions, pseudogenizations, and gene conversions. The expansion of CYP3A in catarrhines (Old World monkeys, great apes, and humans) differed substantially from New World primates (e.g. common marmoset) and strepsirrhines (e.g. galago). We detected two recent episodes of particularly strong positive selection acting on primate CYP3A protein-coding sequence: (i) on CYP3A7 early in hominoid evolution, which was accompanied by a restriction of its hepatic expression to fetal period and (ii) on human CYP3A4 following the split of the chimpanzee and human lineages. In agreement with these findings, three out of four positively selected amino acids investigated in previous biochemical studies of CYP3A affect the activity and regioselectivity.

CONCLUSIONS

CYP3A7 and CYP3A4 may have acquired catalytic functions especially important for the evolution of hominoids and humans, respectively.

Estrogenic modulation of CYP3A38, CYP3A40, and CYP19 in mature male medaka (Oryzias latipes)

We examined cytochrome P450 production and activity and circulating hormone concentrations in male medaka exposed to 17beta-estradiol (E2) or 17alpha-ethinylestradiol (EE2). Intraperitoneal injection of E2 at 1, 10, or 100 microg/g-fish completely suppressed CYP3A38 protein production and suppressed CYP3A40 protein levels by 89%, 52%, or 47%, respectively. CYP3A38 and CYP3A40 mRNA expression was unaltered, and CYP3A enzymatic activity initially increased and then decreased with increasing E2 dose. Males co-cultured with females were exposed to a markedly high concentration (43 ng/L) of E2 secreted by females. CYP3A protein levels in co-cultured males were suppressed. Serum testosterone (TE) and 11keto-testosterone levels in co-cultured males were downregulated to 40% of pre-exposure levels. Serum E2 levels increased in co-cultured males or males exposed to EE2. Testicular CYP19, which converts TE to E2, increased by 9.5 times in males exposed to 50 ng/L EE2 and by 21.5 times in those exposed to 100 ng/L EE2. Male medaka exposed to EE2 showed increased serum Vtg levels. Estrogenic exposure induced Vtg production, suppressed CYP3A protein production, downregulated TE metabolism, and enhanced CYP19 activity. Serum E2 endogenously induced by CYP19 could contribute to Vtg induction in male medaka.

Exposure to p,p'-DDE or dieldrin during the reproductive season alters hepatic CYP expression in largemouth bass (Micropterus salmoides)

Largemouth bass (LMB) in Central Florida living on sites with high levels of organochlorine pesticides (OCPs) have exhibited poor reproductive success and altered steroid profiles. The mechanism underlying these changes is unknown, however changes in the rate of steroid metabolism could alter steroid homeostasis. Members of the CYP2 and CYP3A families play a significant role in the metabolism of many xenobiotics and endogenous compounds, including sex steroids. Therefore, the goal of this study was to identify members of the CYP2 and CYP3A families in LMB and characterize the effects of OCP exposure on their expression. Full-length clones of two CYP3A isoforms were obtained from LMB liver, CYP3A68 and 3A69, which exhibited significant sequence divergence. Full-length clones for CYP2N14 and CYP2P11 were also obtained from LMB liver. Steady-state mRNA levels of each of these CYPs increased in both sexes between early reproductive phase (December) and peak reproductive phase (March). Expression of CYP3A68 and CYP2P11 was sexually dimorphic during peak reproductive phase with 2-fold higher expression in females and males, respectively. Foodborne exposure to 46 ppm p,p'-DDE or 0.8 ppm dieldrin for 30 days did not have a significant effect on expression of CYPs. However, 4 months exposure to p,p'-DDE induced CYP3A68 and 3A69 expression in both sexes, while dieldrin produced weak induction of CYP3A68 and suppressed CYP3A69 expression in females, but had no effect on males. Neither p,p'-DDE nor dieldrin significantly altered the expression of CYP2P11 or CYP2N14. This work demonstrates that there are significant changes in CYP expression that occur during LMB reproduction which can be modified by exposure to OCPs.

Functional characterization of medaka CYP3A38 and CYP3A40: kinetics and catalysis by expression in a recombinant baculovirus system

Phylogenic analysis of the teleost genomic lineages has demonstrated the precedent for multiple genome duplications. Among many of the genes duplicated, cytochrome P450 genes have undergone independent diversification, which can be traced to a single ancestral gene. In teleosts, cytochrome P450s, from all major families, have been identified. Among these, the CYP3A family has been cloned in several teleost species and demonstrated to contain multiple paralogs differing in gene expression patterns and tissue distribution. Herein we characterized the catalytic and kinetic activities of two medaka CYP3A paralogs (CYP3A38 and CYP3A40) with benzyloxyresorufin (BFC), a fluorescent 3A-selective substrate, and testosterone, a known metabolic substrate for CYP3A enzymes. Recombinant CYP3A was produced using the baculovirus expression vector system in Spodoptera frugiperda (Sf9) and Trichoplusia ni (Tn5) insect cells and accounted for up to 24% of total cellular protein. Following addition of a heme-albumin conjugate to log phase cells, spectral P450 content reached a maximum of 560 and 2350 pmol/mg microsomal protein for CYP3A38 and CYP3A40, respectively. Incubations containing recombinant CYP3A, human NADPH-cytochrome P-450 oxidoreductase reductase, human cytochrome b5, and a NADPH generation system catalyzed the dealkylation of BFC and hydroxylation of testosterone with a high degree of stereoselectivity. However, efficiencies and specificities were significantly different between the two isoforms. Km and Vmax activities based on BFC-catalysis were 0.116 and 0.363 muM, and 7.95 and 7.77 nmol/min/nmol P450 for CYP3A38 and CYP3A40, respectively. CYP3A38 preferentially catalyzed testosterone hydroxylation at the 6beta-, 2beta- and 16beta-positions with minor hydroxylation at other positions within the steroid nucleus. Testosterone catalysis with CYP3A40 was limited predominantly to the 6beta- and 2beta-positions. Putative identification of CYP3A substrate recognition sites (SRS) 1-6 indicates that 12 of the 49 amino acid differences between CYP3A38 and CYP3A40 OFRs occur in SRS regions previously known to be associated with steroid hydroxylation. We suggest that differences in kinetics and catalytic activities are a result of amino acid substitutions in SRS regions 1, 3 and 5 within the CYP3A38 and CYP3A40 protein sequence.

Genetically engineered bacterial cells co-expressing human cytochrome P450 with NADPH-cytochrome P450 reductase: prediction of metabolism and toxicity of drugs in humans

Genetically engineered bacterial cells expressing human cytochrome P450 (CYP) have been developed as new tools to predict the metabolism and toxicity of drugs in humans. There are various host cells for the heterologous expression of a form of CYP. Among them, bacterial cells such as Escherichia coli (E. coli) have advantages with regard to ease of use and high yield of protein. CYP protein could be first expressed by the modification of the N-terminal amino acid sequence in E. coli cells in 1991. Since then, many forms of human CYP have been successfully expressed in E. coli cells. Since the E. coli cells do not possess endogeneous electron transport systems to support the full catalytic activity of CYP, E. coli strains co-expressing both human CYP and NADPH-cytochrome P450 reductase (OR) have been established. Each form of CYP expressed in the E. coli cells efficiently catalyzed the oxidation of a representative substrate at an efficient rate, indicating that the OR was sufficiently expressed to support the catalytic activity of CYP. According to the studies performed so far, the modification of the N-terminal amino acid sequence of CYP did not seem to affect the catalytic properties of CYP. The human CYP expressed in the E. coli cells were applicable for studies to determine a metabolic pathway(s) of drugs and to estimate kinetic parameters of drug metabolism by human CYP. Drug-drug interactions caused by inhibition of the metabolism of drugs by human CYP could also be examined by in vitro inhibition studies with CYP expressed in the E. coli cells. Recently, human CYP was co-expressed with the OR in Salmonella typhimurium (S. typhimurium) cells used for mutation assay (Ames test) by applying the technology for the expression of human CYP and the OR in E. coli cells, to evaluate whether chemicals including drugs are metabolically activated by human CYP and show mutagenicity. These strains of bacteria are considered as useful tools to study the metabolism and the toxicity of drugs in humans.

Kinetic and molecular analysis of 5-epiaristolochene 1,3-dihydroxylase, a cytochrome P450 enzyme catalyzing successive hydroxylations of sesquiterpenes

The final step of capsidiol biosynthesis is catalyzed by 5-epiaristolochene dihydroxylase (EAH), a cytochrome P450 enzyme that catalyzes the regio- and stereospecific insertion of two hydroxyl moieties into the bicyclic sesquiterpene 5-epiaristolochene (EA). Detailed kinetic studies using EA and the two possible monohydroxylated intermediates demonstrated the release of 1beta-hydroxy-EA ((OH)EA) at high EA concentrations and a 10-fold catalytic preference for 1beta(OH)EA versus 3alpha(OH)EA, indicative of a preferred reaction order of hydroxylation at C-1, followed by that at C-3. Sequence alignments and homology modeling identified active-site residues tested for their contribution to substrate specificity and overall enzymatic activity. Mutants EAH-S368C and EAH-S368V exhibited wild-type catalytic efficiencies for 1beta(OH)EA biosynthesis, but were devoid of the successive hydroxylation activity for capsidiol biosynthesis. In contrast to EAH-S368C, EAH-S368V catalyzed the relative equal biosynthesis of 1beta(OH)EA, 2beta(OH)EA, and 3beta(OH)EA from EA with wild-type efficiency. Moreover, EAH-S368V converted approximately 1.5% of these monohydroxylated products to their respective ketone forms. Alanine and threonine mutations at position 368 were significantly compromised in their conversion rates of EA to capsidiol and correlated with 3.6- and 5.7-fold increases in their Km values for the 1beta(OH)EA intermediate, respectively. A role for Ile486 in the successive hydroxylations of EA was also suggested by the EAH-I468A mutant, which produced significant amounts 1beta(OH)EA, but negligible amounts of capsidiol from EA. The altered product profile of the EAH-I486A mutant correlated with a 3.6-fold higher Km for EA and a 4.4-fold slower turnover rate (kcat) for 1beta(OH)EA. These kinetic and mutational studies were correlated with substrate docking predictions to suggest how Ser368 and Ile486 might contribute to active-site topology, substrate binding, and substrate presentation to the oxo-Fe-heme reaction center.

CYP3A4 is a vitamin D-24- and 25-hydroxylase: analysis of structure function by site-directed mutagenesis

Studies were performed to identify the microsomal enzyme that 24-hydroxylates vitamin D, whether 25-hydroxylation occurs, and structure function of the enzyme. Sixteen hepatic recombinant microsomal cytochrome P450 enzymes expressed in baculovirus-infected insect cells were screened for 24-hydroxylase activity. CYP3A4, a vitamin D-25-hydroxylase, and CYP1A1 had the highest 24-hydroxylase activity with 1 alpha-hydroxyvitamin D(2) (1 alpha OHD(2)) as substrate. The ratio of rates of 24-hydroxylation of 1 alpha-hydroxyvitamin D(3) (1 alpha OHD(3)), 1 alpha OHD(2), and vitamin D(2) by CYP3A4 was 3.6/2.8/1.0. Structures of 24-hydroxyvitamin D(2), 1,24(S)-dihydroxyvitamin D(2), and 1,24-dihydroxyvitamin D(3) were confirmed by HPLC and gas chromatography retention time and mass spectroscopy. In characterized human liver microsomes, 24-hydroxylation of 1 alpha OHD(2) by CYP3A4 correlated significantly with 6 beta-hydroxylation of testosterone, a marker of CYP3A4 activity. 24-Hydroxylase activity in recombinant CYP3A4 and pooled human liver microsomes showed dose-dependent inhibition by ketoconazole, troleandomycin, alpha-naphthoflavone, and isoniazid, known inhibitors of CYP3A4. Rates of 24- and 25-hydroxylation of 1 alpha OHD(2) and 1 alpha OHD(3) were determined in recombinant wild-type CYP3A4 and site-directed mutants and naturally occurring variants expressed in Escherichia coli. Substitution of residues showed the most prominent alterations of function at residues 119, 120, 301, 305, and 479. Thus, CYP3A4 is both a 24- and 25-hydroxylase for vitamin D(2), 1 alpha OHD(2), and 1 alpha OHD(3).

Characterization of the CYP3A4 active site by homology modeling

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

Analysis of medaka cytochrome P450 3A homotropic and heterotropic cooperativity

We have previously demonstrated that medaka CYP3A is associated with metabolism of several endobiotics including steroids and bile acids. In this study, we demonstrate that medaka CYP3A catalysis exhibits unusual kinetic behaviors consistent with allosteric interaction which cannot be described by hyperbolic kinetic models. Using 7-benzyloxy-4-(trifluoromethyl)-coumarin (BFC) and nonylphenol as CYP3A substrates, we describe both homotropic and heterotropic cooperative interactions. Given the role of CYP3A in maintaining the homeostatic balance for numerous endobiotics, enzymatic activation/inhibition by endocrine disruptors (EDCs) represents a putative (non-genomic) mechanism for endocrine disruption.

Phylogenetic analysis of the cytochrome P450 3 (CYP3) gene family

Cytochrome P450 genes (CYP) constitute a superfamily with members known from the Bacteria, Archaea, and Eukarya. The CYP3 gene family includes the CYP3A and CYP3B subfamilies. Members of the CYP3A subfamily represent the dominant CYP forms expressed in the digestive and respiratory tracts of vertebrates. The CYP3A enzymes metabolize a wide variety of chemically diverse lipophilic organic compounds. To understand vertebrate CYP3 diversity better, we determined the killifish (Fundulus heteroclitus) CYP3A30 and CYP3A56 and the ball python (Python regius) CYP3A42 sequences. We performed phylogenetic analyses of 45 vertebrate CYP3 amino acid sequences using a Bayesian approach. Our analyses indicate that teleost, diapsid, and mammalian CYP3A genes have undergone independent diversification and that the ancestral vertebrate genome contained a single CYP3A gene. Most CYP3A diversity is the product of recent gene duplication events. There is strong support for placement of the guinea pig CYP3A genes within the rodent CYP3A diversification. The rat, mouse, and hamster CYP3A genes are mixed among several rodent CYP3A subclades, indicative of a complex history involving speciation and gene duplication.

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.

Differential oxidation of mifepristone by cytochromes P450 3A4 and 3A5: selective inactivation of P450 3A4

The principal enzyme involved in the oxidation of mifepristone is cytochrome P450 3A4 (CYP3A4), which undergoes mechanism-based inactivation by the drug. However, no information is available on the interaction with CYP3A5, the second most abundant CYP3A enzyme in adult human liver. Oxidation of mifepristone by recombinant CYP3A4 produced mono- and didemethylated products and one C-hydroxylated metabolite, as reported previously. However, CYP3A5 produced only the demethylated metabolites. The apparent V(max) and K(M) values for formation of the monodemethylated product by CYP3A4 and CYP3A5 were 46 and 30 nmol/min/nmol P450, and 36 and 16 microM, respectively. Unlike CYP3A4, CYP3A5 was not inactivated by mifepristone. The basis of this differential susceptibility was explored using site-directed mutants in which a CYP3A4 residue was converted to its 3A5 counterpart. Surprisingly, none of these replacements caused a significant decrease in CYP3A4 inactivation by mifepristone. Examination of selected CYP3A4 mutants at 20 other positions indicated that the relative formation rate of the C-hydroxylated product could not account for the differential susceptibility of CYP3A4 and 3A5. Together these results indicate that mifepristone fails to orient itself in the CYP3A5 active site in such a way that its propylenic group is accessible for oxidation, thus rendering CYP3A5 unable to produce the C-hydroxylated product or putative ketene that leads to enzyme inactivation. Identification of mifepristone as a selective mechanism-based inactivation of CYP3A4 may be very useful in distinguishing between the two major CYP3A enzymes collectively responsible for the oxidative metabolism of over half of the drugs currently in use.

Structural determinants of aromatase cytochrome p450 inhibition in substrate recognition site-1

The porcine gonadal form of aromatase cytochrome P450 (P450arom) exhibits higher sensitivity to inhibition by the imidazole, etomidate, than the placental isozyme. The residue(s) responsible for this functional difference was mapped using chimeragenesis and point mutation analysis of the placental isozyme, and the kinetic analysis was conducted on native and mutant enzymes after overexpression in insect cells. The etomidate sensitivity of the placental isozyme was markedly increased by substitution of the predicted substrate recognition site-1 (SRS-1) and essentially reproduced that of the gonadal isozyme by substitution of SRS-1 and the predicted B helix. A single isoleucine (I) to methionine (M) substitution at position 133 of the placental isozyme (I(133)M) was proven to be the critical residue within SRS-1. Residue 133 is located in the B'-C loop and has been shown to be equally important in other steroid-metabolizing P450s. Single point mutations (including residues 110, 114, 120, 128, 137, and combinations thereof among others) and mutation of the entire B and C helixes were without marked effect on etomidate inhibitory sensitivity. The same mutation (I(133)M) introduced into human P450arom also markedly increased etomidate sensitivity. Mutation of Ile(133) to either alanine (I(133)A) or tyrosine (I(133)Y) decreased apparent enzyme activity, but the I(133)A mutant was sensitive to etomidate inhibition, suggesting that it is Ile(133) that decreases etomidate binding rather than Met(133) increasing enzyme sensitivity. Androstenedione turnover and affinity were similar for the I(133)M mutant and the native placental isozyme. These data suggest that Ile(133) is a contact residue in SRS-1 of P450arom, emphasize the functional conservation that exists in SRS-1 of a number of steroid-hydroxylating P450 enzymes, and suggest that substrate and inhibitor binding are dependent on different contact points to varying degrees.

Modifying rates of reductive elimination of leaving groups from indolequinone prodrugs: a key factor in controlling hypoxia-selective drug release

3-(4-Methylcoumarin-7-yloxy)methylindole-4,7-diones were synthesised as model prodrugs in order to investigate the correlation between rates of reductive elimination from the (indolyl-3-yl)methyl position with reductive metabolism by hypoxic tumor cells and NADPH: cytochrome P450. Rates of elimination of the chromophore/fluorophore (7-hydroxy-4-methylcoumarin) following one-electron reduction of indolequinones to their semiquinone radicals (Q*-) was measured by pulse radiolysis utilising spectrophotometric and fluorometric detection. Incorporation of a thienyl or methyl substituent at the (indol-3-yl)CHR-position (where R=thienyl or methyl adjacent to the phenolic ether linking bond) significantly shortened the half-life of reductive elimination from 87 to 6 and 2 ms, respectively. Elimination from the methyl substituted analogue can thus compete effectively with the reaction of the semiquinone radical with oxygen at levels typically present in tumours (half-life approximately 1.8 ms at 0.5% O2). Chemical kinetic predictions were confirmed by metabolism in breast tumour MCF-7 cells between 0-2.1% O2. Rates of reductive release of the fluorophore from the non-fluorescent parent indolequinones (R=H, Me, thienyl) were similar under anoxia ( approximately 1.7 nmol coumarinmin(-1)mg protein(-1)) reflecting the similarity in one-electron reduction potential. Whereas coumarin release from the indolequinone (R=H) was completely inhibited above 0.5% O2, the enhanced rate of reductive elimination when R=thienyl or Me increased the metabolic rate of release to approximately 0.35 and 0.7 nmol coumarinmin(-1)mg protein(-1), respectively at 0.5% O2; complete inhibition occurring by 2.1% O2. Similar 'oxygen profiles' of release were observed with NADPH: cytochrome P450 reductase. In conclusion, it is possible to modify rates of reductive elimination from indolequinones to control the release of drugs over a range of tumour hypoxia.

CYP3A4 active site volume modification by mutagenesis of leucine 211

The leucine 211 --> phenylalanine (L211F) and leucine 211 --> tyrosine (L211Y) mutant forms of cytochrome P450 3A4 have been generated by site-directed mutagenesis and expressed functionally in Escherichia coli. Substrate binding affinities (S50 values) for testosterone and 7-benzyloxy-4-trifluoromethylcoumarin (BFC) were similar for the mutants and wild-type CYP3A4 (49 and 21 microM for L211F, 35 and 20 microM for L211Y, and 33 and 20 microM for the wild type, respectively). For erythromycin, however, the K(m) values determined for the L211F and L211Y mutants were 2.4- and 10.5-fold higher than for the wild type. Furthermore, IC50 values for the inhibition of testosterone 6 beta-hydroxylation by erythromycin and troleandomycin for L211F were 2.4- and 3.7-fold higher, and those for L211Y were 3.4- and 9.2-fold higher than those measured for the wild type. Conversely, small inhibitors, such as diazepam, exhibited no significant difference in IC50 values between the wild type and the L211F and L211Y mutants. It is proposed that large substrates bound in the catalytic center of CYP3A4 with molecular volumes greater than approximately 600 A(3) were less well accommodated in the altered active sites, resulting in lower association energies and increased IC50 values.

Midazolam oxidation by cytochrome P450 3A4 and active-site mutants: an evaluation of multiple binding sites and of the metabolic pathway that leads to enzyme inactivation

Midazolam (MDZ) oxidation by recombinant CYP3A4 purified from Escherichia coli and 30 mutants generated at 15 different substrate recognition site positions has been studied to determine the role of individual residues in regioselectivity and to investigate the possible existence of multiple binding sites. Initial results showed that oxidation of MDZ by CYP3A4 causes time- and concentration-dependent enzyme inactivation with K(I) and k(inact) values of 5.8 microM and 0.15 min(-1), respectively. The different time courses of MDZ hydroxylation by mutants that predominantly formed 1'-OH MDZ as opposed to 4-OH MDZ provided strong evidence that the 1'-OH MDZ pathway leads to CYP3A4 inactivation. Correlational analysis of 1'-OH formation versus 4-OH formation by the mutants supports the inference that the two metabolites result from the binding of MDZ at two separate sites. Thus, substitution of residues Phe-108, Ile-120, Ile-301, Phe-304, and Thr-309 with a larger amino acid caused an increase in the ratio of 1'-OH/4-OH MDZ formation, whereas substitution of residues Ser-119, Ile-120, Leu-210, Phe-304, Ala-305, Tyr-307, and Thr-309 with a smaller amino acid decreased this ratio. Kinetic analyses of nine key mutants revealed that the alteration in regioselectivity is caused by a change in kinetic parameters (V(max) and K(M)) for the formation of both metabolites in most cases. The study revealed the role of various active-site residues in the regioselectivity of MDZ oxidation, identified the metabolic pathway that leads to enzyme inactivation, and provided an indication that the two proposed MDZ binding sites in CYP3A4 may be partially overlapping.

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.

A comparative approach to structure-function studies of mammalian aromatases

To date, structure--function studies of aromatase cytochrome P450 (P450arom) have been advanced by point mutation analyses utilizing almost exclusively the human enzyme, in conjunction with computer-generated models of the three-dimensional form of the enzyme based on prokaryotic cytochromes P450. Recent studies have identified duplicated isozymes of porcine P450arom, the gonadal and placental forms of which appear to differ substantially in substrate utilization and inhibitor sensitivity. We present a comparative approach to define regions of P450arom responsible for specific functional characteristics using complimentary DNAs encoding the porcine isozymes. Constructs encoding the native and chimeric porcine and human P450arom enzymes were transiently expressed and activity was assessed using the tritiated water assay. Sensitivity to inhibition by the imidazole etomidate was investigated, and P450arom expression was assessed by immunoblot analysis. All constructs yielded active P450arom, suggesting that exchanging entire structural elements does not preclude catalytic function. The activity of the gonadal isozyme was shown to be inhibited by etomidate at concentrations 185 and 300-fold lower than those required to induce a similar inhibition of the placental and human enzymes, respectively. In contrast, there was only a two-fold difference in the sensitivity of the gonadal and placental isozymes to inhibition by CGS16949A. Analysis of chimeric constructs indicated that the sensitivity to etomidate was associated with residues in the B, B' and C helices of the gonadal P450arom encompassing only one of six putative substrate recognition sites. Additionally, sensitivity to etomidate was not correlated with enzyme activity among the chimeric enzymes. Therefore, it appears that residues of the porcine gonadal P450arom that are responsible for etomidate binding may be distinct from those involved in substrate recognition and metabolism. These data support the notion that a comparative approach employing the use of chimeric enzymes provides a useful tool in directing point mutational analysis to determine residues important in inhibitor and perhaps substrate recognition of P450 enzymes such as P450arom. These studies are currently in progress.

The role of cytochrome 2B1 substrate recognition site residues 115, 294, 297, 298, and 362 in the oxidation of steroids and 7-alkoxycoumarins

At least two substitutions were made at each of five amino acid residues in rat cytochrome P450 2B1 that align to residues of known importance in other P450s. The mutants were histidine tagged for purification from Escherichia coli, and the proteins were assessed for testosterone and 7-alkoxycoumarin oxidation. Alteration of each of the sites studied, Phe-115, Ser-294, Phe-297, Ala-298, and Leu-362, was found to affect overall enzyme activity or the metabolite profile. In particular, most of the mutants, excluding F297A, A298G, and L362F, exhibited significantly altered ratios of 16alpha-hydroxytestosterone:16beta-hydroxytestosterone, with the most dramatic alteration being displayed by A298V. Four 7-butoxycoumarin metabolites were produced by CYP2B1, of which two, 7-hydroxycoumarin and 7-(3-hydroxybutoxy)coumarin, were formed at nearly equal rates. Several mutants, F115A, F297A, F297I, and A298V, exhibited an increased predominance of one of the metabolites. The results from this study illustrate the conservation of functionally important residues across P450 subfamilies and families.

Identification and functional characterization of eight CYP3A4 protein variants

The genetic component of the inter-individual variability in CYP3A4 activity has been estimated to be between 60% and 90%, but the underlying genetic factors remain largely unknown. A study of 213 Middle and Western European DNA samples resulted in the identification of 18 new CYP3A4 variants, including eight protein variants. A total of 7.5% of the population studied was found to be heterozygous for one of these variants. In a bacterial heterologous expression system, two mutants, R130Q and P416L, did not result in detectable P450 holoprotein. One mutant, T363M, expressed at significantly lower levels than wild-type CYP3A4. G56D, V170I, D174H and M445T were not significantly different when compared with wild-type CYP3A4 in expression or steroid hydroxylase activity. L373F displayed a significantly altered testosterone metabolite profile and a four-fold increase in the Km value for 1'-OH midazolam formation. The results suggest a limited contribution of CYP3A4 protein variants to the inter-individual variability of CYP3A4 activity in Caucasians. Some variants may, however, play a role in the atypical response to drugs or altered sensitivity to carcinogens.

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.

Cloning and tissue distribution of a novel human cytochrome p450 of the CYP3A subfamily, CYP3A43

On the basis of the detection of an expressed sequence tag ('EST') similar to the human cytochrome P450 3A4 cDNA, we have identified a novel member of the human cytochrome P450 3A subfamily. The coding region is 1512-bp long and shares 84, 83, and 82% sequence identity on the cDNA level with CYP3A4, 3A5, and 3A7, respectively, with a corresponding amino acid identity of 76, 76, and 71%. Quantitative real time based mRNA analysis revealed CYP3A43 expression levels at about 0.1% of CYP3A4 and 2% of CYP3A5 in the liver, with significant expression in 70% of the livers examined. Gene specific PCR of cDNA from extrahepatic tissues showed, with the exception of the testis, only low levels of CYP3A43 expression. The CYP3A43 cDNA was heterologously expressed in yeast, COS-1 cells, mouse hepatic H2.35 cells and in human embryonic kidney (HEK) 293 cells, but in contrast to CYP3A4 which was formed in all cell types, no detectable CYP3A43 protein was produced. This indicates a nonfunctional protein or specific conditions required for proper folding. It is concluded that CYP3A43 mRNA is expressed mainly in liver and testis and that the protein would not contribute significantly to human drug metabolism.