Active-site topologies of human CYP2D6 and its aspartate-301 --> glutamate, asparagine, and glycine mutants

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.

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.

Topological role of cytochrome P450 2D6 active site residues

Recent reports have identified Phe120, Asp301, Thr309, and Glu216 as important residues in cytochrome P450 2D6 (CYP2D6) substrate binding and catalysis. Complementary homology models have located these amino acids within the binding pocket of CYP2D6 and in the present study we have used aryldiazenes to test these models and gain further insight in the role these amino acids have in maintaining the integrity of the active site cavity. When Phe120 was replaced to alanine, there was a significant increase in probe migration to pyrrole nitrogens C and D, in agreement with homology models which have located the phenyl side-chain of Phe120 above these two pyrrole rings. No changes in topology were observed with the D301Q mutant, supporting claims that in this mutant the electrostatic interactions with the B/C-loop are largely maintained and the loop retains its native orientation. The T309V mutation resulted in significant topological alteration suggesting that, in addition to its potential role in dioxygen activation, Thr309 plays an important structural role within the active site crevice. Replacement of Ile106 with Glu, engineered to cause electrostatic repulsion with Glu216, had a profound topological effect in the higher region within the active site cavity and impaired the catalytic activity towards CYP2D6 probe substrates.

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.

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.

TOPOLOGICAL CHANGES IN THE CYP3A4 ACTIVE SITE PROBED WITH PHENYLDIAZENE: EFFECT OF INTERACTION WITH NADPH-CYTOCHROME P450 REDUCTASE AND CYTOCHROMEB5AND OF SITE-DIRECTED MUTAGENESIS

The active site topology of heterologously expressed CYP3A4 purified from an Escherichia coli expression system was examined using phenyldiazene. Incubation of CYP3A4 with phenyldiazene and subsequent oxidation yielded all four potential N-phenylprotoporphyrin IX regioisomers derived from attack on an available nitrogen atom in pyrrole rings B, A, C, or D (N(B):N(A):N(C):N(D) = 6:73:7: 13). Further study using 28 active site mutants showed that substitution of residues closer to the heme, Ala-305, Thr-309, or Ala-370, with a larger residue caused the most drastic changes in regioisomer formation, which reflected the location of each amino acid residue replaced in a CYP3A4 homology model. Previous studies have suggested a conformational change in CYP3A4 upon binding of NADPH-cytochrome P450 reductase (CPR) or cytochrome b(5) (b(5)). Therefore, regioisomer formation was also compared in the absence of redox partners and in the presence of CPR, b(5), or both. Formation of all four regioisomers in CYP3A4 wild type, particularly the minor ones, was reduced in the presence of b(5). CPR also greatly decreased the three minor isomers but increased the major isomer significantly. The presence of b(5) and CPR restored minor isomer formation and suppressed the enhancement of N(A) formation caused by CPR alone. Interestingly, the effects of the redox partners differed among representative active site mutants. In particular, the increase in N(C) upon substitution of Ala-370 with Phe was significantly reversed in the presence of redox partners, strongly suggesting that a conformational change occurs around pyrrole ring C due to protein-protein interactions between CYP3A4 and CPR or b(5).

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.

Influence of lipophilicity on the interactions of N-alkyl-4-phenyl-1,2,3,6-tetrahydropyridines and their positively charged N-alkyl-4-phenylpyridinium metabolites with cytochrome P450 2D6

The relationship between lipophilicity and CYP2D6 affinity of cyclic tertiary (N-alkyl-4-phenyl-1,2,3,6-tetrahydropyridines) and quaternary (N-alkyl-4-phenylpyridinium) amines was examined. The 1,2,3,6-tetrahydropyridine scaffold was chosen due to its common occurrence in the structures of CYP2D6 ligands such as the Parkinsonian neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) and the dehydrated haloperidol metabolite N-[4-(4-fluorophenyl)-4-oxobutyl]-4-(4-chlorophenyl)-1,2,3,6-tetrahydropyridine (HPTP). Likewise, the pyridinium framework is found in and 4-(4-chlorophenyl)-1-[4-(4-fluorophenyl)-4-oxobutyl]pyridinium and N-methyl-4-phenylpyridinium (MPP(+)), the positively charged metabolites of MPTP and haloperidol. The lack of CYP2D6 inhibition by MPTP and its pyridinium metabolite MPP(+) was due to their hydrophilic nature since higher N-alkyl homologs revealed substantial increases in inhibitory potency against recombinant CYP2D6-mediated bufuralol-1'-hydroxylation. The reasonable correlation between lipophilicity and CYP2D6 inhibition by pyridiniums and 1,2,3,6-tetrahydropyridines was only limited to straight chain N-alkyl analogs, since certain N-alkylaryl analogs of lower lipophilicity were better CYP2D6 inhibitors. CYP2D6 substrate properties of straight chain N-alkyltetrahydropyridines were also governed by lipophilicity, and N-heptyl-4-phenyl-1,2,3,6-tetrahydropyridine was the optimal substrate (K(mapp) = 0.63 microM). Metabolism studies indicated that the N-heptyl analog underwent monohydroxylation on the aromatic ring and on the N-heptyl group suggesting that 1,2,3,6-tetrahydropyridines can bind in more than one conformation in the CYP2D6 active site. Increased lipophilicity of haloperidol metabolites did not correlate with inhibitory potency since the more lipophilic HPTP metabolite was less potent as an inhibitor than reduced-haloperidol and reduced-HPTP. Furthermore, HPTP and reduced-HPTP, of comparable lipophilicity to the N-heptyltetrahydropyridine analog were inactive as CYP2D6 substrates. This observation suggests that steric constraints rather than lipophilicity are responsible for the lack of CYP2D6 substrate properties of cyclic tertiary amines tethered to bulky N-substituents. This phenomenon appears to be a common theme among several cyclic tertiary amine-containing anti-depressants and should be taken into consideration when designing central nervous system agents devoid of CYP2D6 substrate properties.

Role of conserved Asp293 of cytochrome P450 2C9 in substrate recognition and catalytic activity

Human cytochrome P450 2C9 (CYP2C9) is important in the metabolism of non-steroidal anti-inflammatory compounds such as diclofenac, the antidiabetic agent tolbutamide and other clinically important drugs, many of which are weakly acidic. Multiple sequence alignment of CYPs identified CYP2C9 Asp(293) as corresponding to Asp(301) of CYP2D6, which has been suggested to play a role in the binding of basic substrates to the latter enzyme. Replacement of Asp(293) with Ala (D293A) decreased activity by more than 90%, and led to an approx. 3- to 10-fold increase in K (m) values for the three test substrates tolbutamide, dextromethorphan and diclofenac. Conservative replacement of the carboxyl side chain in a Glu (D293E) mutant produced no significant changes in K (m) values and slight increases in k (cat) values. Changes in regiospecificity were observed for both the Ala and Glu substitutions; low levels of both dextromethorphan O- and N-demethylation were observed in the D293A mutant, whereas increased preference for O-demethylation was observed for the D293E mutant. Expression of constructs coding for Asn (D293N) and Gln (D293Q) substitutions failed to form a P450 correctly. Our analysis suggests a structural role for the carboxyl side chain of Asp(293) in CYP2C9 substrate binding and catalysis. The conservation of an Asp residue in other CYP families in a position equivalent to Asp(293) indicates a common mechanism for maintaining the active-site architecture.

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.

Expression, purification, and biochemical characterization of a human cytochrome P450 CYP2D6-NADPH cytochrome P450 reductase fusion protein

Cytochrome P450 CYP2D6 metabolizes a wide range of pharmaceutical compounds. A CYP2D6 fusion enzyme (CYP2D6F), containing an amino-terminal human CYP2D6 sequence and a carboxyterminal human NADPH-cytochrome P450 oxidoreductase (CPR) moiety, was constructed. High levels of expression were achieved in Escherichia coli (60-100 nmol/liter) and the enzyme was catalytically active with optimal activities achieved in the presence of the antioxidant, GSH. Turnover values for bufuralol 1'-hydroxylation, metoprolol alpha-hydroxylation, O-desmethylation, and dextromethorphan O-demethylation, using membranes expressing the fusion enzyme, were 5.6, 0.4, 0.72, and 6.19 min(-1), respectively. These values were similar to E. coli membranes which coexpressed human CYP2D6 and CPR (CYP2D6/R). The K(m) and k(cat) values for bufuralol metabolism were estimated to be 10.2 microM and 4.1 min(-1), respectively. The enzyme was purified using ion-exchange chromatography, affinity chromatography (2'-5' ADP-Sepharose), and gel filtration. Estimated turnover rates for bufuralol 1'-hydroxylation, metoprolol alpha-hydroxylation, O-desmethylation, and dextromethorphan O-demethylation were 1.2, 0.52, 0.79, and 0.76 min(-1), respectively. Bufuralol 1'-hydroxylase activity by purified CYP2D6F was enhanced by phospholipids and added CPR. The CYP2D6F enzyme was able to stimulate CYP3A4 testosterone 6beta-hydroxylase activity in a reconstitution system indicating that electron transfer may be largely intermolecular. The catalytically self-sufficient CYP2D6F enzyme will facilitate investigations of P450-CPR interactions and the development of new biocatalysts.

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.

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.

Polymorphic cytochromes P450 and drugs used in psychiatry

1. The cytochrome P450 monooxygenases, CYP2D6, CYP2C19, and CYP2C9, display polymorphism. CYP2D6 and CYP2C19 have been studied extensively, and despite their low abundance in the liver, they catalyze the metabolism of many drugs. 2. CYP2D6 has numerous allelic variants, whereas CYP2C19 has only two. Most variants are translated into inactive, truncated protein or fail to express protein. 3. CYP2C9 is expressed as the wild-type enzyme and has two variants, in each of which one amino acid residue has been replaced. 4. The nucleotide base sequences of the cDNAs of the three polymorphic genes and their variants have been determined, and the proteins derived from these genes have been characterized. 5. An absence of CYP2D6 and/or CYP2C19 in an individual produces a poor metabolizer (PM) of drugs that are substrates of these enzymes. 6. When two drugs that are substrates for a polymorphic CYP enzyme are administered concomitantly, each will compete for that enzyme and competitively inhibit the metabolism of the other substrate. This can result in toxicity. 7. Patients can be readily phenotyped or genotyped to determine their CYP2D6 or CYP2C19 enzymatic status. Poor metabolizers (PMs), extensive metabolizers (EMs), and ultrarapid metabolizers (URMs) can be identified. 8. Numerous substrates and inhibitors of CYP2D6, CYP2C19, and CYP2C9 are identified. 9. An individual's diet and age can influence CYP enzyme activity. 10. CYP2D6 polymorphism has been associated with the risk of onset of various illnesses, including cancer, schizophrenia, Parkinson's disease, Alzheimer's disease, and epilepsy.

Pyrimidinylimidazole inhibitors of CSBP/p38 kinase demonstrating decreased inhibition of hepatic cytochrome P450 enzymes

Pyrimidine analogs of the pyrimidinylimidazole class of CSBP/p38 kinase inhibitors were prepared in an effort to reduce the potent inhibition of hepatic cytochrome P450 observed for the pyridinyl compounds. The substitution of pyrimidin-4-yl, 2-methoxypyrimidin-4-yl, or 2-methylaminopyrimidin-4-yl for pyridin-4-yl effectively dissociates CSBP/p38 kinase from P450 inhibition for this series and furthermore achieves an increase in oral activity.

Structural determinants of cytochrome P450 substrate specificity, binding affinity and catalytic rate

The structural characteristics of cytochrome P450 substrates are summarised, showing that molecular descriptors can discriminate between chemicals of differing P450 isozyme specificity. Procedures for the estimation of P450 substrate binding interaction energies and rates of metabolism are described, providing specific examples in both individual compounds binding to P450s, including those of known crystal structure, and within series of structurally related chemicals. It is demonstrated that binding energy components are primarily hydrophobic/desolvation and electrostatic/hydrogen-bonded in nature, whereas electronic factors are of importance in determining variations in reaction rates. It is thus shown that the prediction of P450 substrate binding affinities and catalytic rates may be feasible, provided that sufficient structural information is available for the relevant enzyme-substrate complex.

Cytochrome P450cam substrate specificity: relationship between structure and catalytic oxidation of alkylbenzenes

The oxidation by cytochrome P450cam (CYP101) of ethylbenzene and a series of substrates derived from it by addition of one, two, three, or four carbon atoms has been examined. For each of the 18 substrates, the shift in spin state due to substrate binding, the extent of coupled turnover to give organic products and uncoupled turnover to give H2O2 and H2O, and the identities of the organic products have been determined. The same studies have been carried out with the T185L and T185F mutants of P450cam in which the active site volume is decreased. The results show that no detectable correlation exists between the observed spin state change and any other parameter studied. For substrates of equal size, coupled and uncoupled turnover vary widely but both are maximized when the phenyl ring bears one large alkyl substituent or a methyl ortho to the largest alkyl substituent. The presence of substituents in addition to these decreases activity. In the absence of other changes, coupled turnover is correlated with the size of the largest substituent, but no such correlation exists for uncoupled turnover. Decreasing the size of the active site cavity by a T185L mutation generally increases coupled turnover without altering the dependence on the alkyl group size. A T185F mutation causes too great an active site perturbation for structure-activity studies. Substrate oxidation occurs preferentially at 2 degrees or 3 degrees C-H bonds of the largest substituent or on the benzylic methyl ortho to it. Aromatic hydroxylation only competes with the oxidation of nonbenzylic 1 degree C-H bonds. The extent of coupled turnover is a function of substrate shape, substrate size, and cavity size, but still elusive parameters control the extent of uncoupled turnover.

Analysis of four residues within substrate recognition site 4 of human cytochrome P450 3A4: role in steroid hydroxylase activity and alpha-naphthoflavone stimulation

Sequence alignment of human cytochrome P450 3A4 with bacterial enzymes of known structure has provided a basis from which to predict residues involved in substrate oxidation. Substitutions were made at four residues (I301, F304, A305, and T309) predicted to be located within the highly conserved substrate recognition site 4. Site-directed mutants engineered to contain carboxy-terminal histidine tags were expressed in Escherichia coli and purified on a metal affinity column. The integrity of each protein was assessed by SDS-polyacrylamide gel electrophoresis and immunoblotting. Functional analysis was performed using progesterone and testosterone as substrates and alpha-naphthoflavone as an activator. In testosterone hydroxylase assays, all of the mutants displayed rates of total product formation similar to wild-type 3A4, with several mutants showing small differences in specific products formed. However, with progesterone as the substrate, mutants F304A, A305V, and T309A exhibited altered product ratios and/or changes in the rates of product formation. F304A and A305V also displayed altered flavonoid stimulation that resulted in product ratios dramatically different from wild-type 3A4. Therefore, the kinetics of progesterone hydroxylation of these mutants and the wild-type enzyme were further assessed, and the data were analyzed with the Hill equation. Results with wild-type 3A4 and F304A indicated that at high progesterone concentrations, hydroxylation rates and product ratios are independent of the presence of alpha-NF. This suggests that progesterone may be equivalent to alpha-NF as an activator. In contrast, A305V exhibited autoactivation by progesterone but inhibition by alpha-NF.

Ab initio molecular modeling in the study of drug metabolism

We discuss the use of ab initio quantum mechanical methods in drug metabolism studies. These methods require only the positions and atomic numbers of the atoms to be specified and offer greater transferability than conventional molecular modeling techniques. This fact, coupled with the accuracy of our approach, permits 'computational experiments' to be performed, allowing details of reaction mechanisms to be understood. We review the application of these methods to the cytochrome P450 superfamily of enzymes. There is much interest in understanding the mechanisms of these enzymes due to their participation in a wide range of metabolic processes including drug activation/deactivation. We find that our methods accurately reproduce the low- to high-spin transition of the haem Fe on binding of a substrate. Furthermore, we identify a new mechanism for the suppression of this spin transition, namely the shortening of the bond between the Fe atom and the coordinated S atom of the cysteine axial ligand. These results indicate that ab initio molecular modeling may be usefully applied in the study of drug metabolism and that further study of intermediate states in the P450 reaction cycle would be beneficial, particularly those which are not accessible using conventional experimental approaches.