Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity

Use of engineered cytochromes P450 for accelerating drug discovery and development

Numerous steps in drug development, including the generation of authentic metabolites and late-stage functionalization of candidates, necessitate the modification of often complex molecules, such as natural products. While it can be challenging to make the required regio- and stereoselective alterations to a molecule using purely chemical catalysis, enzymes can introduce changes to complex molecules with a high degree of stereo- and regioselectivity. Cytochrome P450 enzymes are biocatalysts of unequalled versatility, capable of regio- and stereoselective functionalization of unactivated CH bonds by monooxygenation. Collectively they catalyze over 60 different biotransformations on structurally and functionally diverse organic molecules, including natural products, drugs, steroids, organic acids and other lipophilic molecules. This catalytic versatility and substrate range makes them likely candidates for application as potential biocatalysts for industrial chemistry. However, several aspects of the P450 catalytic cycle and other characteristics have limited their implementation to date in industry, including: their lability at elevated temperature, in the presence of solvents, and over lengthy incubation times; the typically low efficiency with which they metabolize non-natural substrates; and their lack of specificity for a single metabolic pathway. Protein engineering by rational design or directed evolution provides a way to engineer P450s for industrial use. Here we review the progress made to date toward engineering the properties of P450s, especially eukaryotic forms, for industrial application, and including the recent expansion of their catalytic repertoire to include non-natural reactions.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Compensatory changes in CYP expression in three different toxicology mouse models: CAR-null, Cyp3a-null, and Cyp2b9/10/13-null mice

Targeted mutant models are common in mechanistic toxicology experiments investigating the absorption, metabolism, distribution, or elimination (ADME) of chemicals from individuals. Key models include those for xenosensing transcription factors and cytochrome P450s (CYP). Here we investigated changes in transcript levels, protein expression, and steroid hydroxylation of several xenobiotic detoxifying CYPs in constitutive androstane receptor (CAR)-null and two CYP-null mouse models that have subfamily members regulated by CAR; the Cyp3a-null and a newly described Cyp2b9/10/13-null mouse model. Compensatory changes in CYP expression that occur in these models may also occur in polymorphic humans, or may complicate interpretation of ADME studies performed using these models. The loss of CAR causes significant changes in several CYPs probably due to loss of CAR-mediated constitutive regulation of these CYPs. Expression and activity changes include significant repression of Cyp2a and Cyp2b members with corresponding drops in 6α- and 16β-testosterone hydroxylase activity. Further, the ratio of 6α-/15α-hydroxylase activity, a biomarker of sexual dimorphism in the liver, indicates masculinization of female CAR-null mice, suggesting a role for CAR in the regulation of sexually dimorphic liver CYP profiles. The loss of Cyp3a causes fewer changes than CAR. Nevertheless, there are compensatory changes including gender-specific increases in Cyp2a and Cyp2b. Cyp2a and Cyp2b were down-regulated in CAR-null mice, suggesting activation of CAR and potentially PXR following loss of the Cyp3a members. However, the loss of Cyp2b causes few changes in hepatic CYP transcript levels and almost no significant compensatory changes in protein expression or activity with the possible exception of 6α-hydroxylase activity. This lack of a compensatory response in the Cyp2b9/10/13-null mice is probably due to low CYP2B hepatic expression, especially in male mice. Overall, compensatory and regulatory CYP changes followed the order CAR-null > Cyp3a-null > Cyp2b-null mice.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Repression of multiple CYP2D genes in mouse primary hepatocytes with a single siRNA construct

The Cyp2d subfamily is the second most abun-dant subfamily of hepatic drug-metabolizing CYPs. In mice, there are nine Cyp2d members that are believed to have redundant catalytic activity. We are testing and optimizing the ability of one short interfering RNA (siRNA) construct to knockdown the expression of multiple mouse Cyp2ds in primary hepatocytes. Expression of Cyp2d10, Cyp2d11, Cyp2d22, and Cyp2d26 was observed in the primary male mouse hepatocytes. Cyp2d9, which is male-specific and growth hormone-dependent, was not expressed in male primary hepatocytes, potentially because of its dependence on pulsatile growth hormone release from the anterior pituitary. Several different siRNAs at different concentrations and with different reagents were used to knockdown Cyp2d expression. siRNA constructs designed to repress only one construct often mildly repressed several Cyp2d isoforms. A construct designed to knockdown every Cyp2d isoform provided the best results, especially when incubated with transfection reagents designed specifically for primary cell culture. Interestingly, a construct designed to knockdown all Cyp2d isoforms, except Cyp2d10, caused a 2.5× increase in Cyp2d10 expression, presumably because of a compensatory response. However, while RNA expression is repressed 24 h after siRNA treatment, associated changes in Cyp2d-mediated metabolism are tenuous. Overall, this study provides data on the expression of murine Cyp2ds in primary cell lines, valuable information on designing siRNAs for silencing multiple murine CYPs, and potential pros and cons of using siRNA as a tool for repressing Cyp2d and estimating Cyp2d's role in murine xenobiotic metabolism.

Double site saturation mutagenesis of the human cytochrome P450 2D6 results in regioselective steroid hydroxylation

The human cytochrome P450 2D6 (CYP2D6) is one of the major human drug metabolizing enzymes and acts preferably on substrates containing a basic nitrogen atom. Testosterone - just as other steroids - is an atypical substrate and only poorly metabolized by CYP2D6. The present study intended to investigate the influence of the two active site residues 216 and 483 on the capability of CYP2D6 to hydroxylate steroids such as for example testosterone. All 400 possible combinatorial mutations at these two positions have been generated and expressed individually in Pichia pastoris. Employing whole-cell biotransformations coupled with HPLC-MS analysis the testosterone hydroxylase activity and regioselectivity of every single CYP2D6 variant was determined. Covering the whole sequence space, CYP2D6 variants with improved activity and so far unknown regio-preference in testosterone hydroxylation were identified. Most intriguingly and in contrast to previous literature reports about mutein F483I, the mutation F483G led to preferred hydroxylation at the 2β-position, while the slow formation of 6β-hydroxytestosterone, the main product of wild-type CYP2D6, was further reduced. Two point mutations have already been sufficient to convert CYP2D6 into a steroid hydroxylase with the highest ever reported testosterone hydroxylation rate for this enzyme, which is of the same order of magnitude as for the conversion of the standard substrate bufuralol by wild-type CYP2D6. Furthermore, this study is also an example for efficient human CYP engineering in P. pastoris for biocatalytic applications and to study so far unknown pharmacokinetic effects of individual and combined mutations in these key enzymes of the human drug metabolism.

Active site substitution A82W improves the regioselectivity of steroid hydroxylation by cytochrome P450 BM3 mutants as rationalized by spin relaxation nuclear magnetic resonance studies

Cytochrome P450 BM3 from Bacillus megaterium is a monooxygenase with great potential for biotechnological applications. In this paper, we present engineered drug-metabolizing P450 BM3 mutants as a novel tool for regioselective hydroxylation of steroids at position 16β. In particular, we show that by replacing alanine at position 82 with a tryptophan in P450 BM3 mutants M01 and M11, the selectivity toward 16β-hydroxylation for both testosterone and norethisterone was strongly increased. The A82W mutation led to a ≤42-fold increase in V(max) for 16β-hydroxylation of these steroids. Moreover, this mutation improves the coupling efficiency of the enzyme, which might be explained by a more efficient exclusion of water from the active site. The substrate affinity for testosterone increased at least 9-fold in M11 with tryptophan at position 82. A change in the orientation of testosterone in the M11 A82W mutant as compared to the orientation in M11 was observed by T(1) paramagnetic relaxation nuclear magnetic resonance. Testosterone is oriented in M11 with both the A- and D-ring protons closest to the heme iron. Substituting alanine at position 82 with tryptophan results in increased A-ring proton-iron distances, consistent with the relative decrease in the level of A-ring hydroxylation at position 2β.

Modeling human cytochrome P450 2D6 metabolism and drug-drug interaction by a novel panel of knockout and humanized mouse lines

The highly polymorphic human cytochrome P450 2D6 enzyme is involved in the metabolism of up to 25% of all marketed drugs and accounts for significant individual differences in response to CYP2D6 substrates. Because of the differences in the multiplicity and substrate specificity of CYP2D family members among species, it is difficult to predict pathways of human CYP2D6-dependent drug metabolism on the basis of animal studies. To create animal models that reflect the human situation more closely and that allow an in vivo assessment of the consequences of differential CYP2D6 drug metabolism, we have developed a novel straightforward approach to delete the entire murine Cyp2d gene cluster and replace it with allelic variants of human CYP2D6. By using this approach, we have generated mouse lines expressing the two frequent human protein isoforms CYP2D6.1 and CYP2D6.2 and an as yet undescribed variant of this enzyme, as well as a Cyp2d cluster knockout mouse. We demonstrate that the various transgenic mouse lines cover a wide spectrum of different human CYP2D6 metabolizer phenotypes. The novel humanization strategy described here provides a robust approach for the expression of different CYP2D6 allelic variants in transgenic mice and thus can help to evaluate potential CYP2D6-dependent interindividual differences in drug response in the context of personalized medicine.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Engineering and assaying of cytochrome P450 biocatalysts

Cytochrome P450s constitute a highly fascinating superfamily of enzymes which catalyze a broad range of reactions. They are essential for drug metabolism and promise industrial applications in biotechnology and biosensing. The constant search for cytochrome P450 enzymes with enhanced catalytic performances has generated a large body of research. This review will concentrate on two key aspects related to the identification and improvement of cytochrome P450 biocatalysts, namely the engineering and assaying of these enzymes. To this end, recent advances in cytochrome P450 development are reported and commonly used screening methods are surveyed.

Functional expression and comparative characterization of nine murine cytochromes P450 by fluorescent inhibition screening

The increasing number of transgenic or gene knockout mouse models generated for use in drug metabolism studies has meant that a greater understanding of the function and substrate specificities of murine cytochromes P450 (P450s) has become essential, particularly with the recent advances in "humanized" mouse models. In this study, we have heterologously expressed nine murine P450s--Cyp1a1, Cyp1a2, Cyp1b1, Cyp2a4, Cyp2b20, Cyp2c29, Cyp2d22, Cyp2e1, and Cyp3a11--individually with human P450 oxidoreductase to generate functional monooxygenase systems in Escherichia coli. We have identified a suitable fluorogenic probe for each P450 and determined the apparent kinetic parameters. These probes have enabled the screening of a panel of 31 test compounds classified as "drugs," "natural compounds," "endogenous compounds," and "pesticides" by measurement of IC(50), thus allowing the comparison of binding affinities. Human P450s CYP2C9, CYP2D6, and CYP3A4 were also included in the study to enable direct comparisons to be made with the mouse enzymes. Although there were general similarities between human and mouse P450s, perhaps the most significant finding in this study was the observation that, despite 77% amino acid identity, Cyp2d22 and CYP2D6 were remarkably dissimilar in a range of enzymatic properties, with potentially serious implications for pharmacokinetic studies using CYP2D substrates. The data presented in this study provide a solid foundation with which to assess the degree of similarity (or difference) between mouse and human P450s involved in xenobiotic metabolism and can be used as a basis for further studies.

Analysis of CYP2D6 substrate interactions by computational methods

Cytochrome P450 CYP2D6 is involved in the oxidation of well over 150 drugs and, in general, those which contain a basic nitrogen atom in the molecule. To clarify how the residues of CYP2D6 are utilized for orientating a wide range of its specific substrates and distinguishing them from a variety of other organic compounds, docking studies by AutoDock and molecular dynamics (MD) simulations were conducted. Specific ligands were docked to both the homology model and crystal structures optimally to estimate the site of reaction on the ligand molecule and the binding energy for the complex, which were generally in good agreement with the experimental data. MD simulation for the CYP2D6-propranolol complex was then carried out to reveal the amino acid residues interacting with the substrate at the active site. Phe-120, Glu-216, Asp-301, and Phe-483 are identified as the substrate-binding residues in agreement with previously reported site-directed mutagenesis data and the crystal structure reported recently (PDB code: 2F9Q). As well as these residues, our theoretical prediction suggests that Phe-219 and Glu-222 are also important residues for mediating oxidation of substrates, especially propranolol.

Insights into drug metabolism by cytochromes P450 from modelling studies of CYP2D6-drug interactions

The cytochromes P450 (CYPs) comprise a vast superfamily of enzymes found in virtually all life forms. In mammals, xenobiotic metabolizing CYPs provide crucial protection from the effects of exposure to a wide variety of chemicals, including environmental toxins and therapeutic drugs. Ideally, the information on the possible metabolism by CYPs required during drug development would be obtained from crystal structures of all the CYPs of interest. For some years only crystal structures of distantly related bacterial CYPs were available and homology modelling techniques were used to bridge the gap and produce structural models of human CYPs, and thereby obtain useful functional information. A significant step forward in the reliability of these models came seven years ago with the first crystal structure of a mammalian CYP, rabbit CYP2C5, followed by the structures of six human enzymes, CYP1A2, CYP2A6, CYP2C8, CYP2C9, CYP2D6 and CYP3A4, and a second rabbit enzyme, CYP2B4. In this review we describe as a case study the evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism. This work has led directly to the successful design of CYP2D6 mutants with novel activity-including creating a testosterone hydroxylase, converting quinidine from inhibitor to substrate, creating a diclofenac hydroxylase and creating a dextromethorphan O-demethylase. Our modelling-derived hypothesis-driven integrated interdisciplinary studies have given key insight into the molecular determinants of CYP2D6 and other important drug metabolizing enzymes.

IN SILICO AND IN VITRO SCREENING FOR INHIBITION OF CYTOCHROME P450 CYP3A4 BY COMEDICATIONS COMMONLY USED BY PATIENTS WITH CANCER

Cytochrome P450 3A4 (CYP3A4) is the major enzyme responsible for phase I drug metabolism of many anticancer agents. It is also a major route for metabolism of many drugs used by patients to treat the symptoms caused by cancer and its treatment as well as their other illnesses, for example, cardiovascular disease. To assess the ability to inhibit CYP3A4 of drugs most commonly used by our patients during cancer therapy, we have made in silico predictions based on the crystal structures of CYP3A4. From this set of 33 common comedicated drugs, 10 were predicted to be inhibitors of CYP3A4, with the antidiarrheal drug loperamide predicted to be the most potent. There was significant correlation (r(2) = 0.75-0.66) between predicted affinity and our measured IC(50) values, and loperamide was confirmed as a potent inhibitor (IC(50) of 0.050 +/- 0.006 microM). Active site docking studies predicted an orientation of loperamide consistent with formation of the major (N-demethylated) metabolite, where it interacts with the phenylalanine cluster and Arg-212 and Glu-374; experimental evidence for the latter interaction comes from the approximately 12-fold increase in K(M) for loperamide observed for the Glu-374-Gln mutant. The commonly prescribed drugs loperamide, amitriptyline, diltiazem, domperidone, lansoprazole, omeprazole, and simvastatin were identified by our in silico and in vitro screens as relatively potent inhibitors of CYP3A4 that have the potential to interact with cytotoxic agents to cause adverse effects, highlighting the likelihood of drug-drug interactions affecting chemotherapy treatment.

Why Is Quinidine an Inhibitor of Cytochrome P450 2D6?

We have previously shown that Phe(120), Glu(216), and Asp(301) in the active site of cytochrome P450 2D6 (CYP2D6) play a key role in substrate recognition by this important drug-metabolizing enzyme (Paine, M. J., McLaughlin, L. A., Flanagan, J. U., Kemp, C. A., Sutcliffe, M. J., Roberts, G. C., and Wolf, C. R. (2003) J. Biol. Chem. 278, 4021-4027 and Flanagan, J. U., Maréchal, J.-D., Ward, R., Kemp, C. A., McLaughlin, L. A., Sutcliffe, M. J., Roberts, G. C., Paine, M. J., and Wolf, C. R. (2004) Biochem. J. 380, 353-360). We have now examined the effect of mutations of these residues on interactions of the enzyme with the prototypical CYP2D6 inhibitor, quinidine. Abolition of the negative charge at either or both residues 216 and 301 decreased quinidine inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation by at least 100-fold. The apparent dissociation constants (K(d)) for quinidine binding to the wild-type enzyme and the E216D and D301E mutants were 0.25-0.50 microm. The amide substitution of Glu(216) or Asp(301) resulted in 30-64-fold increases in the K(d) for quinidine. The double mutant E216Q/D301Q showed the largest decrease in quinidine affinity, with a K(d) of 65 microm. Alanine substitution of Phe(120), Phe(481),or Phe(483) had only a minor effect on the inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation and on binding. In contrast to the wild-type enzyme, a number of the mutants studied were found to be able to metabolize quinidine. E216F produced O-demethylated quinidine, and F120A and E216Q/D301Q produced both O-demethylated quinidine and 3-hydroxyquinidine metabolites. Homology modeling and molecular docking were used to predict the modes of quinidine binding to the wild-type and mutant enzymes; these were able to rationalize the experimental observations.

The role of phenylalanine 483 in cytochrome P450 2D6 is strongly substrate dependent

The polymorphic cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 30% of the drugs currently prescribed, and is thus clinically relevant. Typical CYP2D6 substrates generally contain a basic nitrogen atom and an aromatic moiety adjacent to the site of metabolism. Recently, we demonstrated the importance of active site residue F120 in substrate binding and catalysis in CYP2D6. On the basis of protein homology models, it is claimed that another active site phenylalanine, F483, may also play an important role in the interaction with the aromatic moiety of CYP2D6 substrates. Experimental data to support this hypothesis, however, is not yet available. In fact, in the only study performed, mutation of F483 to isoleucine or tryptophan did not affect the 1'-hydroxylation of bufuralol at all [Smith G, Modi S, Pillai I, Lian LY, Sutcliffe MJ, Pritchard MP, et al., Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity. Biochem J 1998;331:783-92]. In the present study, the role of F483 in ligand binding and metabolism by CYP2D6 was examined experimentally using site-directed mutagenesis. Replacement of F483 by alanine resulted in a 30-fold lower V(max) for bufuralol 1'-hydroxylation, while the K(m) was hardly affected. The V(max) for 3,4-methylenedioxy-methylamphetamine O-demethylenation on the other hand decreased only two-fold, whereas the effect on the K(m) was much larger. For dextromethorphan, in addition to dextrorphan (O-demethylation) and 3-methoxymorphinan (N-demethylation), two other metabolites were formed that could not be detected for the wild-type. The substrate 7-methoxy-4-(aminomethyl)-coumarin was not metabolised at all by CYP2D6[F483A], a phenomenon that was reported also for CYP2D6[F120A]. The presented data show that next to F120, residue F483 plays a very important role in the metabolism of typical CYP2D6 substrates. The influence of F483 on metabolism was found to be strongly substrate-dependent.

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.

Two Novel Single Nucleotide Polymorphisms (SNPs) of the CYP2D6 Gene in Japanese Individuals

We analyzed all the exons and exon-intron junctions of the CYP2D6 gene from 286 Japanese individuals. We detected two novel single nucleotide polymorphisms (SNPs) 2556C>T in exon 5 (Thr261Ile) and 3835A>C in exon 8 (Lys404Gln). Both these SNPs showed a frequency of 0.002.

Influence of phenylalanine 120 on cytochrome P450 2D6 catalytic selectivity and regiospecificity: crucial role in 7-methoxy-4-(aminomethyl)-coumarin metabolism

The polymorphic human debrisoquine hydroxylase, cytochrome P450 2D6 (CYP2D6), is one of the most important phase I drug metabolising enzymes. It is responsible for metabolising a large number of compounds that mostly share similarity in having a basic N-atom and an aromatic moiety. In homology modelling studies, it has been suggested that in fixation of this aromatic moiety, there may be an important role for phenylalanine 120 (Phe(120)). In this study, the role of Phe(120) in ligand binding and catalysis was experimentally examined by mutating it into an alanine. Strikingly, this substitution led to a completely abolished 7-methoxy-4-(aminomethyl)-coumarin (MAMC) O-demethylating activity of CYP2D6. On the other hand, bufuralol metabolism was hardly affected (K(m) of 1-hydroxylation mutant: 1.2 microM, wild-type: 2.9 microM, 4-hydroxylation mutant: 1.5 microM, and wild-type: 3.2 microM) and neither was affected dextromethorphan O-demethylation (K(m) mutant: 1.2 microM, wild-type: 2 microM, k(cat) mutant: 4.5 min(-1), and wild-type: 3.3 min(-1)). However, the Phe(120)Ala mutant also formed 3-hydroxymorphinan, the double demethylated form of dextromethorphan, which was not detected using wild-type CYP2D6. 3,4-Methylenedioxymethamphetamine (MDMA) was demethylenated by both mutant and wild-type CYP2D6 to 3,4-dihydroxymethamphetamine (3,4-OH-MA K(m) of mutant: 55 microM and wild-type: 2 microM). In addition, the mutant formed two additional metabolites; 3,4-methylenedioxyamphetamine (MDA) and N-hydroxy-3,4-methylenedioxymethamphetamine (N-OH-MDMA). Inhibition experiments of dextromethorphan O-demethylation showed a decreased affinity of the Phe(120)Ala mutant for quinidine (IC(50) mutant: 240 nM and wild-type, 40 nM), while IC(50)s for quinine were equal (1 microM). These data indicate the importance of Phe(120) in the selectivity and regiospecificity in substrate binding and catalysis by CYP2D6.

Phe120 contributes to the regiospecificity of cytochrome P450 2D6: mutation leads to the formation of a novel dextromethorphan metabolite

Although the residues that determine the preference of CYP2D6 (cytochrome P450 2D6) for compounds containing a basic nitrogen are well characterized, the contribution of other active site residues to substrate binding and orientation is less well understood. Our structural model of CYP2D6 identifies the aromatic residue Phe120 as a likely major feature of the active site. To examine the role of Phe120, mutants of CYP2D6 in which this residue has been substituted by alanine, leucine, tyrosine, serine, histidine, tryptophan or methionine residues have been prepared in bacterial membranes co-expressing human cytochrome NADPH cytochrome P450 oxidoreductase. The mutants have been characterized using the prototypical bufuralol 1' hydroxylase and dextromethorphan O- and N-demethylase activities of CYP2D6. Larger effects on K(m) values are observed for dextromethorphan O-demethylation than for bufuralol 1' hydroxylation, indicating that the Phe120 side chain is more important in dextromethorphan than in bufuralol binding. A role for this side chain in determining the regiospecificity of substrate oxidation was indicated by changes in the relative rates of O- and N-demethylation of dextromethorphan and, notably, by the formation of 7-hydroxy dextromethrophan, a novel dextromethorphan metabolite, in mutants in which it had been substituted. Computational studies of dextromethorphan binding to the active site of the Phe120-->Ala mutant were carried out to throw light on the way in which the removal of this side chain leads to different modes of ligand binding.

Comparative molecular field analysis and QSAR on substrates binding to cytochrome P450 2D6

In this study, we utilized comparative molecular field analysis (CoMFA) to gain a better understanding of the steric and electrostatic features of the cytochrome p450 2D6 (CYP2D6) active site. The training set consists of 24 substrates with reported K(M) values from liver microsomal CYP2D6 spanning an activity range of almost three log units. The low energy conformers were fit by root mean square (RMS) to minaprine at the site of metabolism and to the protonated nitrogen. In this manner, we constructed two CoMFA models, one model with a distance constraint and another without. The model with the distance parameter (non-cross-validated R(2)=0.99) was approximately equal to the CoMFA without a distance parameter (non-cross-validated R(2)=0.98). Validation of our CoMFA was accomplished by predicting the K(M) values of 15 diverse CYP2D6 substrates not in the original training set resulting in a predictive R(2)=0.62. Finally, we also pursued correlations of pK(a) and log P with CYP2D6 substrate K(M) in an effort to investigate other physicochemical properties.

Behavioral effects of intraventricular injections of low doses of ethanol, acetaldehyde, and acetate in rats: studies with low and high rate operant schedules

Although ethanol is typically classed as a sedative-hypnotic, low doses of ethanol have been shown to stimulate locomotor activity in mice. However, in rats the typical response to peripheral administration of ethanol is a dose-dependent suppression of motor activity and operant responding. The present study was undertaken to determine the effects of intraventricular (ICV) infusions of ethanol, acetaldehyde, and acetate on operant performance in rats. ICV injections of ethanol, acetaldehyde, or acetate were given to rats previously trained on either a differential-reinforcement-of-low-rates-of-responding (DRL) 30-s schedule, which generates low rates of responding, or a fixed ratio 5 (FR5) schedule, which generates relatively high rates. Ethanol, acetaldehyde, and acetate all produced a rate-increasing effect in rats on the DRL 30-s schedule at moderate doses (2.8 and 1.4 micromol, respectively). Acetate also produced a rate-decreasing effect on the DRL 30-s schedule at a larger dose (8.8 micromol). Performance on the FR5 schedule was unaltered by ethanol and acetaldehyde, even at doses as high as 17.6 micromol. However, acetate produced a rate-decreasing effect on the FR5 schedule at doses of 4.4, 5.6, and 8.8 micromol. Central administration of low doses of ethanol and its metabolites can increase operant responding on some schedules in rats. Acetate is the substance that is most potent for producing rate-suppressing effects. These results indicate that the major metabolites of ethanol are pharmacologically active when injected into the brain, and suggest that acetate may mediate some of the rate-suppressing effects of ethanol, such as sedation, ataxia or motor slowing.

A discordance between cytochrome P450 2D6 genotype and phenotype in patients undergoing methadone maintenance treatment

AIMS

To assess CYP2D6 activity and genotype in a group of patients undergoing methadone maintenance treatment (MMT).

METHODS

Blood samples from 34 MMT patients were genotyped by a polymerase chain reaction-based method, and results were compared with CYP2D6 phenotype (n = 28), as measured by the molar metabolic ratio (MR) of dextromethorphan (DEX)/dextrorphan (DOR) in plasma.

RESULTS

Whereas 9% of patients (3/34) were poor metabolizers (PM) by genotype, 57% (16/28) were PM by phenotype (P < 0.005). Eight patients, who were genotypically extensive metabolizers (EM), were assigned as PM by their phenotype. The number of CYP2D6*4 alleles and sex were significant determinants of CYP2D6 activity in MMT patients, whereas other covariates (methadone dose, age, weight) did not contribute to variation in CYP2D6 activity.

CONCLUSIONS

There was a discordance between genotype and in vivo CYP2D6 activity in MMT patients. This finding is consistent with inhibition of CYP2D6 activity by methadone and may have implications for the safety and efficacy of other CYP2D6 substrates taken by MMT patients.

Cytochrome P450 inhibition using recombinant proteins and mass spectrometry/multiple reaction monitoring technology in a cassette incubation

Detailed cytochrome P450 (P450) inhibition profiles are now required for the registration of novel molecular entities. This method uses combined substrates (phenacetin, diclofenac, S-mephenytoin, bufuralol, and midazolam) with combined recombinant P450 enzymes (CYP1A2, 2C9, 2C19, 2D6, and 3A4) in an attempt to limit interactions with other more minor P450s and associated reductases. Kinetic analysis of single substrate with single P450 (sP450) yielded apparent Km values of 25, 2, 20, 9, and 3 microM, for CYP1A2, 2C9, 2C19, 2D6, and 3A4, respectively. Combined substrates with combined P450s (cP450) yielded apparent Km values of 65, 4, 19, 7, and 2 microM. Selectivity of the substrates for each P450 isoform was checked. Phenacetin proved to be the least selective substrate. However, the ratio of the various P450s was modified in the final assay such that metabolism of phenacetin by other enzymes was approximately 20% of the metabolism by CYP1A2. IC50 determinations with alpha-naphthoflavone (0.04 microM), sulfaphenazole (0.26 microM), tranylcypromine (9 microM), quinidine (0.02 microM), and ketoconazole (0.01 microM) were similar for sP450 and cP450 enzymes. The assay was further evaluated with 11 literature compounds and 52 in-house new chemical entities, and the data compared with radiometric/fluorescent values. The overall protein level of the assay was reduced from the original starting point, as this led to some artificially high IC50 measurements when compared with existing lower protein assays (radiometric/fluorometric). This method offers high throughput P450 inhibition profiling with potential advantages over current radiometric or fluorometric methods.

Residues glutamate 216 and aspartate 301 are key determinants of substrate specificity and product regioselectivity in cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes a wide range of therapeutic drugs. CYP2D6 substrates typically contain a basic nitrogen atom, and the active-site residue Asp-301 has been implicated in substrate recognition through electrostatic interactions. Our recent computational models point to a predominantly structural role for Asp-301 in loop positioning (Kirton, S. B., Kemp, C. A., Tomkinson, N. P., St.-Gallay, S., and Sutcliffe, M. J. (2002) Proteins 49, 216-231) and suggest a second acidic residue, Glu-216, as a key determinant in the binding of basic substrates. We have evaluated the role of Glu-216 in substrate recognition, along with Asp-301, by site-directed mutagenesis. Reversal of the Glu-216 charge to Lys or substitution with neutral residues (Gln, Phe, or Leu) greatly decreased the affinity (K(m) values increased 10-100-fold) for the classical basic nitrogen-containing substrates bufuralol and dextromethorphan. Altered binding was also manifested in significant differences in regiospecificity with respect to dextromethorphan, producing enzymes with no preference for N-demethylation versus O-demethylation (E216K and E216F). Neutralization of Asp-301 to Gln and Asn had similarly profound effects on substrate binding and regioselectivity. Intriguingly, removal of the negative charge from either 216 or 301 produced enzymes (E216A, E216K, and D301Q) with elevated levels (50-75-fold) of catalytic activity toward diclofenac, a carboxylate-containing CYP2C9 substrate that lacks a basic nitrogen atom. Activity was increased still further (>1000-fold) upon neutralization of both residues (E216Q/D301Q). The kinetic parameters for diclofenac (K(m) 108 microm, k(cat) 5 min(-1)) along with nifedipine (K(m) 28 microm, k(cat) 2 min(-1)) and tolbutamide (K(m) 315 microm, k(cat) 1 min(-1)), which are not normally substrates for CYP2D6, were within an order of magnitude of those observed with CYP3A4 or CYP2C9. Neutralizing both Glu-216 and Asp-301 thus effectively alters substrate recognition illustrating the central role of the negative charges provided by both residues in defining the specificity of CYP2D6 toward substrates containing a basic nitrogen.

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.

A MOLECULAR MODEL OF CYP2D6 CONSTRUCTED BY HOMOLOGY WITH THE CYP2C5 CRYSTALLOGRAPHIC TEMPLATE: INVESTIGATION OF ENZYME-SUBSTRATE INTERACTIONS

The results of homology modelling of CYP2D6 based on the mammalian P450 crystal structure of rabbit CYP2C5 are reported. It is found that many CYP2D6-selective substrates are able to fit closely within the putative active site of the enzyme where there are favourable contacts with complementary amino acid residues, including aspartate-301 which has been probed via site-directed mutagenesis. The homology model of CYP2D6 is consistent with available experimental evidence from selective substrate metabolism and site-specific mutation data. Quantitative structure-activity relationships (QSARs) with substrate binding affinity based on KD values and inhibition data (Ki values) demonstrate the importance of hydrogen bonding, pi-pi stacking and relative molecular mass in describing variations in avidity towards the CYP2D6 enzyme, although the compound lipophilicity (log D(7.4)) appears to be the most important single descriptor for CYP2D6 inhibition. Calculation of substrate binding affinity based on contributions from active site interactions and lipophilic character gives satisfactory agreement with experimentally determined KD values.

Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes approximately one third of the drugs in current clinical use. To gain insight into its structure and function, we have produced four different sets of comparative models of 2D6: one based on the structures of P450s from four different microorganisms (P450 terp, P450 eryF, P450 cam, and P450 BM3), another on the only mammalian P450 (2C5) structure available, and the other two based on alternative amino acid sequence alignments of 2D6 with all five of these structures. Principal component analysis suggests that inclusion of the 2C5 crystal structure has a profound effect on the modeling process, altering the general topology of the active site, and that the models produced differ significantly from all of the templates. The four models of 2D6 were also used in conjunction with molecular docking to produce complexes with the substrates codeine and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP); this identified Glu 216 [in the F-helix; substrate recognition site (SRS) 2] as a key determinant in the binding of the basic moiety of the substrate. Our studies suggest that both Asp 301 and Glu 216 are required for metabolism of basic substrates. Furthermore, they suggest that Asp 301 (I-helix, SRS-4), a residue thought from mutagenesis studies to bind directly to the basic moiety of substrates, may play a key role in positioning the B'-C loop (SRS-1) and that the loss of activity on mutating Asp 301 may therefore be the result of an indirect effect (movement of the B'-C loop) on replacing this residue.

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

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

Comparative modelling of cytochromes P450

The superfamily of enzymes known as the cytochromes P450 (P450s) comprises a wide-ranging class of proteins with diverse functions. They are known, amongst other things, to be involved in the hormonal regulation of metabolism and reproduction, as well as having a major clinical significance through their association with diseases such as cancer, diabetes and hepatitis. Knowledge of the three-dimensional (3D) structure of a protein gives insight into its function. The 3D structures of P450s are therefore of considerable scientific interest. A number of high-resolution structures of P450s have been determined by X-ray crystallography and studies of these structures have provided valuable insights into the mechanism of these enzymes. Only one of these structures is mammalian and as yet there is no structural information on human P450s in the public domain. Until such a structure is solved it is necessary to employ alternative methods to gain structural insight into how human P450s perform their biological function. Here we report on the use of comparative modelling to predict the structure of human P450s based on knowledge of their amino acid sequences plus the 3D structures of other (not human) P450s. As an illustrative example of these techniques we have modelled the structure of P450 2C5 using five bacterial P450 structures as templates. We examine the importance of selecting suitable templates, obtaining a good amino acid sequence alignment, and evaluating the models generated. To improve the quality of the models an iterative cycle of sequence alignment, model building, and model evaluation is employed. The result is a model with excellent stereochemistry, good amino acid side chain environment properties, and a Calpha trace similar to the crystal structure.

Detection of cytochrome P450 mRNA transcripts in prostate samples by RT-PCR

BACKGROUND

The expression of cytochrome P450 (CYP)-dependent mono-oxygenases in the prostate is important, as it will determine the rate of activation of potential carcinogens as well as the metabolism of hormones with implications in diseases of the prostate. In addition, the levels of cytochromes P450 in prostatic tumours may well be determinants of the outcome of therapy involving P450 substrates such as anti-androgens.

METHODS

The gene expression of 12 different CYP genes was measured by reverse transcription-polymerase chain reaction (RT-PCR) in a total of 28 human prostatic tumour and nontumour samples.

RESULTS

Intriguingly, a large number of CYP mRNAs were detected in the prostate samples, including CYP1A2, -1B1, -2C19, -2D6, -3A4, -3A5, -3A7 and -4B1. CYP1B1 was consistently expressed and CYP3A5 and CYP4B1 were expressed in a majority of the samples tested.

CONCLUSIONS

These data demonstrate a wide range of CYP genes being expressed in the prostate. The relative importance of these enzymes in the pathogenesis and treatment of prostatic disease remains an important theme for further study.

Heterologous expression of cytochrome P450 2D6 mutants, electron transfer, and catalysis of bufuralol hydroxylation: the role of aspartate 301 in structural integrity

Cytochrome P450 (P450) 2D6 is a polymorphic human enzyme involved in the oxidation of >50 drugs, most of which contain a basic nitrogen. In confirmation of previous work by others, substitutions at Asp301 decreased rates of substrate oxidation by P450 2D6. An anionic residue (Asp, Glu) at this position was found to be important in proper protein folding and heme incorporation, and positively charged residues were particularly disruptive in bacterial and also in baculovirus expression systems. Truncation of 20 N-terminal amino acids had no significant effect on catalytic activity except to attenuate P450 2D6 interaction with membranes and NADPH-P450 reductase. The truncation of the N-terminus increased the level of bacterial expression of wild-type P450 2D6 (Asp301) but markedly reduced expression of all codon 301 mutants, including Glu301. Reduction of ferric P450 2D6 by NADPH-P450 reductase was enhanced in the presence of the prototypic substrate bufuralol. Bacterial flavodoxin, an NADPH-P450 reductase homolog, binds tightly to P450 2D6 but is inefficient in electron transfer to the heme. These results collectively indicate that the acidic residue at position 301 in P450 2D6 has a structural role in addition to any in substrate binding and that the N-terminus of P450 2D6 is relatively unimportant to catalytic activity beyond a role in facilitating binding to NADPH-P450 reductase.

Progesterone oxidation by cytochrome P450 2D isoforms in the brain

The existence of cytochrome P450 2D isoforms in the brain has been demonstrated, although their physiological functions remain to be elucidated. In this study we demonstrated that recombinant rat cytochrome P450 2D1 and 2D4 and human cytochrome P450 2D6 possess progesterone 6 beta- and 16 alpha- hydroxylation activities; 2 beta- and 21-hydroxylation activities; and 2 beta-, 6 beta-, 16 alpha- and 21-hydroxylation activities, respectively. Cytochrome P450 2D4 had the lowest K(m) value and the highest maximum velocity value toward these activities. Progesterone 2 beta- and 21-hydroxylation activities were also detected in rat brain microsomes, and these activities were completely inhibited by anticytochrome P450 2D antibodies. The presence of endogenous 2 beta- and 21-hydroxyprogesterones in rat brain tissues was also demonstrated. The mRNAs of cytochrome P450 2D4, CYP11A, and 3 beta-hydroxysteroid dehydrogenase were detected in the rat brain, suggesting that progesterone was generated from cholesterol by CYP11A and 3 beta-hydroxysteroid dehydrogenase and then underwent hydroxylation to hydroxyprogesterones by cytochrome P450 2D4 in rat brain. Collectively, our findings support the idea that cytochrome P450 2D may be involved in the regulation (metabolism and/or synthesis) of endogenous neuroactive steroids, such as progesterone and its derivatives, in brain tissues.

Influence of phenylalanine-481 substitutions on the catalytic activity of cytochrome P450 2D6

Homology models of the active site of cytochrome P450 2D6 (CYP2D6) have identified phenylalanine 481 (Phe(481)) as a putative ligand-binding residue, its aromatic side chain being potentially capable of participating in pi-pi interactions with the benzene ring of ligands. We have tested this hypothesis by replacing Phe(481) with tyrosine (Phe(481)-->Tyr), a conservative substitution, and with leucine (Phe(481)-->Leu) or glycine (Phe(481)-->Gly), two non-aromatic residues, and have compared the properties of the wild-type and mutant enzymes in microsomes prepared from yeast cells expressing the appropriate cDNA-derived protein. The Phe(481)-->Tyr substitution did not alter the kinetics [K(m) (microM) and V(max) (pmol/min per pmol) respectively] of oxidation of S-metoprolol (27; 4.60), debrisoquine (46; 2.46) or dextromethorphan (2; 8.43) relative to the respective wild-type values [S-metoprolol (26; 3.48), debrisoquine (51; 3.20) and dextromethorphan (2; 8.16)]. The binding capacities [K(s) (microM)] of a range of CYP2D6 ligands to the Phe(481)-->Tyr enzyme (S-metoprolol, 22.8; debrisoquine, 12.5; dextromethorphan, 2.3; quinidine, 0.13) were also similar to those for the wild-type enzyme (S-metoprolol, 10.9; debrisoquine, 8.9; dextromethorphan, 3.1; quinidine, 0.10). In contrast, the Phe(481)-->Leu and Phe(481)-->Gly substitutions increased significantly (3-16-fold) the K(m) values of oxidation of the three substrates [S-metoprolol (120-124 microM), debrisoquine (152-184 microM) and dextromethorphan (20-31 microM)]. Similarly, the K(s) values of the ligands to Phe(481)-->Leu and Phe(481)-->Gly mutants were also increased 3 to 10-fold (S-metoprolol, 33.2-41.9 microM; debrisoquine, 85-90 microM; dextromethorphan, 15.7-18.8 microM; quinidine 0.35-0.53 microM). However, contrary to a recent proposal that Phe(481) has the dominant role in the binding of substrates that undergo CYP2D6-mediated N-dealkylation routes of metabolism, the Phe(481)-->Gly substitution did not substantially decrease the capacity of the enzyme to N-deisopropylate metoprolol (wild-type, 1.12 pmol/min per pmol of P450; Phe(481)-->Gly, 0.71), whereas an Asp(301)-->Gly substitution decreased the N-dealkylation reaction by 95% of the wild-type rate. Overall, our results are consistent with the proposal that Phe(481) is a ligand-binding residue in the active site of CYP2D6 and that the residue interacts with ligands via a pi-pi interaction between its phenyl ring and the aromatic moiety of the ligand. However, the relative importance of Phe(481) in binding is ligand-dependent; furthermore, its importance is secondary to that of Asp(301). Finally, contrary to predictions of a recent homology model, Phe(481) does not seem to have a primary role in CYP2D6-mediated N-dealkylation.

Porcine microsomal vitamin D(3) 25-hydroxylase (CYP2D25). Catalytic properties, tissue distribution, and comparison with human CYP2D6

The metabolic activation of the prohormone vitamin D(3) requires a 25-hydroxylation that has been reported to be catalyzed by both mitochondrial CYP27A and a microsomal vitamin D(3) 25-hydroxylase in the liver. CYP27A has been extensively studied, but its role as a physiologically important vitamin D(3) 25-hydroxylase has been questioned. The present paper reports that the microsomal vitamin D(3) 25-hydroxylase, purified from pig liver, converted vitamin D(3) into 25-hydroxyvitamin D(3) in substrate concentrations which are within the physiological range (apparent K(m) = 0.1 microm). The enzyme 25-hydroxylated vitamin D(3), 1 alpha-hydroxyvitamin D(3) and vitamin D(2) and also converted tolterodine, a substrate for human CYP2D6, into its 5-hydroxymethyl metabolite. Tolterodine inhibited the microsomal 25-hydroxylation, whereas quinidine, an inhibitor of CYP2D6, did not markedly inhibit the reaction. The primary structure of the microsomal vitamin D(3) 25-hydroxylase, designated CYP2D25, shows 77% identity with that of human CYP2D6. Northern blot and reverse transcription-polymerase chain reaction experiments revealed that CYP2D25 mRNA is expressed in higher levels in liver than in kidney and in small amounts in adrenals, brain, heart, intestine, lung, muscle, spleen, and thymus. Experiments with human liver microsomes and recombinantly expressed CYP2D6 strongly indicate that the microsomal 25-hydroxylation of vitamin D(3) in human liver is catalyzed by an enzyme different from CYP2D6.

Characterization of Cyp2d22, a novel cytochrome P450 expressed in mouse mammary cells

Endogenous steroids and numerous environmental agents have potent effects on mammary development and carcinogenesis. Locally produced cytochrome P450 enzymes that modify such molecules are therefore likely to be important regulators of these processes. Here we describe the characterization of a novel mouse gene, termed Cyp2d22, that is highly expressed in the mammary tumor derived cell line RIII/Prl. Cyp2d22 is expressed at intermediate levels in the weakly tumorigenic cell line RIII/MG, whereas expression is low or absent in all normal mouse mammary epithelial cell lines tested and three C3H mammary tumor derived cell lines. Immunoblot analysis of mouse tissues with highly specific antisera indicates that 2D22 protein levels are most abundant in liver, while intermediate levels of expression are seen in adrenal, ovary, and mammary gland. Immunohistochemical staining of liver sections with these antisera demonstrates that 2D22 is most abundant in the first layer or two of parenchymal cells surrounding the central vein, with virtually no expression detected in periportal cells. Interestingly, sequence similarity and functional data suggest that Cyp2d22 may be the mouse ortholog of human CYP2D6. These observations support the hypothesis that 2D22 mediates a distinct, biologically significant activity in relation to other mouse 2D family members.

The active site of the thermophilic CYP119 from Sulfolobus solfataricus

CYP119 from Sulfolobus solfataricus, the first thermophilic cytochrome P450, is stable at up to 85 degrees C. UV-visible and resonance Raman show the enzyme is in the low spin state and only modestly shifts to the high spin state at higher temperatures. Styrene only causes a small spin state shift, but T(1) NMR studies confirm that styrene is bound in the active site. CYP119 catalyzes the H(2)O(2)-dependent epoxidation of styrene, cis-beta-methylstyrene, and cis-stilbene with retention of stereochemistry. This catalytic activity is stable to preincubation at 80 degrees C for 90 min. Site-specific mutagenesis shows that Thr-213 is catalytically important and Thr-214 helps to control the iron spin state. Topological analysis by reaction with aryldiazenes shows that Thr-213 lies above pyrrole rings A and B and is close to the iron atom, whereas Thr-214 is some distance away. CYP119 is very slowly reduced by putidaredoxin and putidaredoxin reductase, but these proteins support catalytic turnover of the Thr-214 mutants. Protein melting curves indicate that the thermal stability of CYP119 does not depend on the iron spin state or the active site architecture defined by the threonine residues. Independence of thermal stability from active site structural factors should facilitate the engineering of novel thermostable catalysts.

Theoretical investigation of substrate specificity for cytochromes P450 IA2, P450 IID6 and P450 IIIA4

Three-dimensional models of the cytochromes P450 IA2, P450 IID6 and P450 IIIA4 were built by means of comparative modeling using the X-ray crystallographic structures of P450 CAM, P450 BM-3, P450 TERP and P450 ERYF as templates. The three cytochromes were analyzed both in their intrinsic structural features and in their interaction properties with fifty specific and non-specific substrates. Substrate/enzyme complexes were obtained by means of both automated rigid and flexible body docking. The comparative analysis of the three cytochromes and the selected substrates, in their free and bound forms, allowed for the building of semi-quantitative models of substrate specificity based on both molecular and intermolecular interaction descriptors. The results of this study provide new insights into the molecular determinants of substrate specificity for the three different eukaryotic P450 isozymes and constitute a useful tool for predicting the specificity of new compounds.