Crystal Structure of Human Cytochrome P450 2D6

Functional analysis of phenylalanine residues in the active site of cytochrome P450 2C9

The two published crystal structures of cytochrome P450 2C9, complexed with ( S)-warfarin or flurbiprofen, implicate a cluster of three active site phenylalanine residues (F100, F114, F476) in ligand binding. However, these three residues appear to interact differently with these two ligands based on the static crystal structures. To elucidate the importance of CYP2C9's active site phenylalanines on substrate binding, orientation, and catalytic turnover, a series of leucine and tryptophan mutants were constructed and their interactions with ( S)-warfarin and ( S)-flurbiprofen examined. The F100-->L mutation had minor effects on substrate binding and metabolism of each substrate. In contrast, the F114L and F476L mutants exhibited substantially reduced ( S)-warfarin metabolism and altered hydroxy metabolite profiles but only modestly decreased nonsteroidal antiinflammatory drug (NSAID) turnover while maintaining product regioselectivity. The F114-->W and F476-->W mutations also had opposing effects on ( S)-warfarin versus NSAID turnover. Notably, the F476W mutant increased the efficiency of ( S)-warfarin metabolism 5-fold, yet decreased the efficiency of ( S)-flurbiprofen turnover 20-fold. (1)H NMR T 1 relaxation studies suggested a slightly closer positioning of ( S)-warfarin to the heme in the F476W mutant relative to the wild-type enzyme, and stoichiometry studies indicated enhanced coupling of reducing equivalents to product formation for ( S)-warfarin, again in contrast to effects observed with ( S)-flurbiprofen. These data demonstrate that F114 and F476, but not F100, influence ( S)-warfarin's catalytic orientation. Differential interactions of F476 mutants with the two substrates suggest that their catalytically productive binding modes are not superimposable.

Human cytochromes P450 in the metabolism of drugs: new molecular models of enzyme-substrate interactions

The overall predictive ability of molecular modelling, as applied to the cytochrome P450 (CYP) system, is analysed in the light of current developments in a variety of techniques, including X-ray crystallography, molecular biology, enzyme kinetics, molecular mechanics and dynamics, in relation to its relevance to drug metabolism in humans. This review demonstrates that it is possible to generate realistic models for the major human CYPs, which metabolise xenobiotics that compare favourably with crystal structures, and thus may be used to derive substrate binding energies that agree closely with experimental K(m) values obtained from enzyme kinetics.

Quantitation of human cytochrome P450 2D6 protein with immunoblot and mass spectrometry analysis

Accurate quantification of cytochrome P450 (P450) protein contents is essential for reliable assessment of drug safety, including the prediction of in vivo clearance from in vitro metabolism data, which may be hampered by the use of uncharacterized standards and existence of unknown allelic isozymes. Therefore, this study aimed to delineate the variability in absolute quantification of polymorphic CYP2D6 drug-metabolizing enzyme and compare immunoblot and nano liquid chromatography coupled to mass spectrometry (nano-LC/MS) methods in identification and relative quantification of CYP2D6.1 and CYP2D6.2 allelic isozymes. Holoprotein content of in-house purified CYP2D6 isozymes was determined according to carbon monoxide difference spectrum, and total protein was quantified with bicinchoninic acid protein assay. Holoprotein/total CYP2D6 protein ratio was markedly higher for purified CYP2D6.1 (71.0%) than that calculated for CYP2D6.1 Supersomes (35.5%), resulting in distinct linear calibration range (0.05-0.50 versus 0.025-0.25 pmol) that was determined by densitometric analysis of immunoblot bands. Likewise, purified CYP2D6.2 and CYP2D6.10 and the CYP2D6.10 Supersomes all showed different holoprotein/total CYP2D6 protein ratios and distinct immunoblot linear calibration ranges. In contrast to immunoblot, nano-LC/MS readily distinguished CYP2D6.2 (R296C and S486T) from CYP2D6.1 by isoform-specific proteolytic peptides that contain the altered amino acid residues. In addition, relative quantitation of the two allelic isozymes was successfully achieved with label-free protein quantification, consistent with the nominated ratio. Because immunoblot and nano-LC/MS analyses measure total P450 protein (holoprotein and apoprotein) in a sample, complete understanding of holoprotein and apoprotein contents in P450 standards is desired toward reliable quantification. Our data also suggest that nano-LC/MS not only facilitates P450 quantitation but also provides genotypic information.

Structures of human cytochrome P-450 2E1. Insights into the binding of inhibitors and both small molecular weight and fatty acid substrates

Human microsomal cytochrome P-450 2E1 (CYP2E1) monooxygenates > 70 low molecular weight xenobiotic compounds, as well as much larger endogenous fatty acid signaling molecules such as arachidonic acid. In the process, CYP2E1 can generate toxic or carcinogenic compounds, as occurs with acetaminophen overdose, nitrosamines in cigarette smoke, and reactive oxygen species from uncoupled catalysis. Thus, the diverse roles that CYP2E1 has in normal physiology, toxicity, and drug metabolism are related to its ability to metabolize diverse classes of ligands, but the structural basis for this was previously unknown. Structures of human CYP2E1 have been solved to 2.2 angstroms for an indazole complex and 2.6 angstroms for a 4-methylpyrazole complex. Both inhibitors bind to the heme iron and hydrogen bond to Thr303 within the active site. Complementing its small molecular weight substrates, the hydrophobic CYP2E1 active site is the smallest yet observed for a human cytochrome P-450. The CYP2E1 active site also has two adjacent voids: one enclosed above the I helix and the other forming a channel to the protein surface. Minor repositioning of the Phe478 aromatic ring that separates the active site and access channel would allow the carboxylate of fatty acid substrates to interact with conserved 216QXXNN220 residues in the access channel while positioning the hydrocarbon terminus in the active site, consistent with experimentally observed omega-1 hydroxylation of saturated fatty acids. Thus, these structures provide insights into the ability of CYP2E1 to effectively bind and metabolize both small molecule substrates and fatty acids.

AmineDB: Large scale docking of amines with CYP2D6 and scoring for druglike properties—towards defining the scope of the chemical defense against foreign amines in humans

Organic amines are prevalent in nature and in drugs, especially the psychotherapeutic agents, and a major defense against potentially toxic amines is metabolism by CYP2D6. In order to understand better the constraints on the broad specificity of CYP2D6, 4207 amines were docked into the binding site of this enzyme. Docking poses were found predominantly with the positively charged amino groups closest to Asp301, with aromatic rings close to Phe120 and sometimes extending as far as Phe483. Organic amines that bind best to CYP2D6 tend to have larger molecular weights and logP values. Organic amines that score highly as being druglike, based on a Bayesian model constructed using a 5223-drug training set, are least likely to bind to CYP2D6. This correlation suggests that the set of known drugs, which have been largely designed or selected to avoid high affinity CYP binding, partially encodes the binding site preferences (or rather anti-preferences) of CYP2D6. Finally, in order to benchmark our docking and druglike scoring procedures, an analysis of psychotherapeutic agents is presented. All of these data, including the 4207 AM1-optimized ligand structures in proper ionization states, docking poses and scores, Druglike Scores and Lipinski properties, can be viewed from an online database, the AmineDB.

Modes of heme binding and substrate access for cytochrome P450 CYP74A revealed by crystal structures of allene oxide synthase

Cytochrome P450s exist ubiquitously in all organisms and are involved in many biological processes. Allene oxide synthase (AOS) is a P450 enzyme that plays a key role in the biosynthesis of oxylipin jasmonates, which are involved in signal and defense reactions in higher plants. The crystal structures of guayule (Parthenium argentatum) AOS (CYP74A2) and its complex with the substrate analog 13(S)-hydroxyoctadeca-9Z,11E-dienoic acid have been determined. The structures exhibit a classic P450 fold but possess a heme-binding mode with an unusually long heme binding loop and a unique I-helix. The structures also reveal two channels through which substrate and product may access and leave the active site. The entrances are defined by a loop between beta3-2 and beta3-3. Asn-276 in the substrate binding site may interact with the substrate's hydroperoxy group and play an important role in catalysis, and Lys-282 at the entrance may control substrate access and binding. These studies provide both structural insights into AOS and related P450s and a structural basis to understand the distinct reaction mechanism.

CYP53A15 of Cochliobolus lunatus, a target for natural antifungal compounds

A novel cytochrome P450, CYP53A15, was identified in the pathogenic filamentous ascomycete Cochliobolus lunatus. The protein, classified into the CYP53 family, was capable of para hydroxylation of benzoate. Benzoate is a key intermediate in the metabolism of aromatic compounds in fungi and yet basically toxic to the organism. To guide functional analyses, protein structure was predicted by homology modeling. Since many naturally occurring antifungal phenolic compounds are structurally similar to CYP53A15 substrates, we tested their putative binding into the active site of CYP53A15. Some of these compounds inhibited CYP53A15. Increased antifungal activity was observed when tested in the presence of benzoate. Some results suggest that CYP53A15 O-demethylation activity is important in detoxification of other antifungal substances. With the design of potent inhibitors, CYP53 enzymes could serve as alternative antifungal drug targets.

Application of molecular modeling for prediction of substrate specificity in cytochrome P450 1A2 mutants

Molecular dynamics (MD) simulations of 7-ethoxy- and 7-methoxyresorufin bound in the active site of cytochrome P450 (P450) 1A2 wild-type and various mutants were used to predict changes in substrate specificity of the mutants. A total of 26 multiple mutants representing all possible combinations of five key amino acid residues, which are different between P450 1A1 and 1A2, were examined. The resorufin substrates were docked in the active site of each enzyme in the productive binding orientation, and MD simulations were performed on the enzyme-substrate complex. Ensembles collected from MD trajectories were then scored on the basis of geometric parameters relating substrate position with respect to the activated oxoheme cofactor. The results showed a high correlation between the previous experimental data on P450 1A2 wild-type and single mutants with respect to the ratio between 7-ethoxyresorufin-O-deethylase (EROD) and 7-methoxyresorufin-O-demethylase (MROD) activities and the equivalent in silico "E/M scores" (the ratio of hits obtained with 7-ethoxyresorufin to those obtained with 7-methoxyresorufin). Moreover, this correlation served to establish linear regression models used to evaluate E/M scores of multiple P450 1A2 mutants. Seven mutants, all of them incorporating the L382V substitution, were predicted to shift specificity to that of P450 1A1. The predictions were then verified experimentally. The appropriate P450 1A2 multiple mutants were constructed by site-directed mutagenesis, expressed in Escherichia coli, and assayed for EROD and MROD activities. Of six mutants, five demonstrated an increased EROD/MROD ratio, confirming modeling predictions.

The endocannabinoid anandamide is a substrate for the human polymorphic cytochrome P450 2D6

Members of the cytochrome P450 (P450) family of drug-metabolizing enzymes are present in the human brain, and they may have important roles in the oxidation of endogenous substrates. The polymorphic CYP2D6 is one of the major brain P450 isoforms and has been implicated in neurodegeneration, psychosis, schizophrenia, and personality traits. The objective of this study was to determine whether the endocannabinoid arachidonoylethanolamide (anandamide) is a substrate for CYP2D6. Anandamide is the endogenous ligand to the cannabinoid receptor CB1, which is also activated by the main psychoactive component in marijuana. Signaling via the CB1 receptor alters sensory and motor function, cognition, and emotion. Recombinant CYP2D6 converted anandamide to 20-hydroxyeicosatetraenoic acid ethanolamide and 5,6-, 8,9-, 11,12-, and 14,15-epoxyeicosatrienoic acid ethanolamides (EET-EAs) with low micromolar K(m) values. CYP2D6 further metabolized the epoxides of anandamide to form novel dioxygenated derivatives. Human brain microsomal and mitochondrial preparations metabolized anandamide to form hydroxylated and epoxygenated products, respectively. An inhibitory antibody against CYP2D6 significantly decreased the mitochondrial formation of the EET-EAs. To our knowledge, anandamide and its epoxides are the first eicosanoid-like molecules to be identified as CYP2D6 substrates. Our study suggests that anandamide may be a physiological substrate for brain mitochondrial CYP2D6, implicating this polymorphic enzyme as a potential component of the endocannabinoid system in the brain. This study also offers support to the hypothesis that neuropsychiatric phenotype differences among individuals with genetic variations in CYP2D6 could be ascribable to interactions of this enzyme with endogenous substrates.

The development of drug metabolism research as expressed in the publications of ASPET: Part 3, 1984-2008

The dramatic changes in drug metabolism research in the last 25 years are well documented in the publications of the American Society for Pharmacology and Experimental Therapeutics (ASPET). New analytical tools combined with modern molecular biological techniques have provided unprecedented access to the workings of the cell. A field that concentrated on only a handful of primary enzymes now has a list of hundreds in its purview. Genetic variation, environmental impact, and molecular diversity have all come under study in attempts to follow the fate of drugs and chemicals. Examples from ASPET journals will be used to illustrate the dramatic advancements in the field.

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.

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

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

Crystallization and preliminary X-ray analysis of allene oxide synthase, cytochrome P450 CYP74A2, from Parthenium argentatum

Oxylipins are oxygenated derivatives of fatty acids and pivotal signaling molecules in plants and animals. Allene oxide synthase (AOS) is a key cytochrome P450 CYP74 enzyme involved in the biosynthesis of plant oxylipin jasmonates to convert 13(S)-hydroperoxide to allene oxide. Guayule (Parthenium argentatum) AOS, CYP74A2, was expressed in Escherichia coli. Protein was purified using affinity chromatography and size exclusion chromatography, and then crystallized. Two different crystal forms were obtained from 0.2 M (NH(4))H(2)PO(4), 50% MPD, 0.1 M Tris, pH 8.5 at 277 K using the hanging-drop vapor-diffusion method. Preliminary X-ray analysis was carried out, and the crystals were found to belong to the tetragonal space group I422 with cell parameters a = b = 126.5, c = 163.9 A, and the monoclinic space group C2 with cell parameters a = 336.5, b = 184.2, c = 159.0 A, beta = 118.6 degrees . Diffraction data were collected to 2.4 A resolution from a tetragonal form of crystal using a home X-ray source.

Efficient macromolecular crystallization using microfluidics and randomized design of screening reagents

Microfluidic technologies enable a relatively new approach to macromolecular crystallization, but offer several significant advantages over more traditional techniques. Microfluidic devices provide significant savings in the amount of material required to complete a set of experiments, although recent innovations with vapor diffusion and microbatch methods have also greatly reduced their material requirements. When compared with these other methods, microfluidic approaches still consume 5-100x less material. In addition, comparisons in one set of experiments suggest that microfluidic free-interface diffusion may also offer substantially higher success rates than sitting drop vapor diffusion. Microfluidic methods also provide opportunities for experimental strategies involving testing multiple samples in parallel. When combined with randomized design of screening reagents, microfluidic devices provide a highly efficient method for sampling crystallization space. Commercial microfluidic crystallization chips have been in circulation for a number of years now and stable protocols for their use, tips and tricks, and data on their success and failure are now available.

Molecular characterization of CYP2B6 substrates

CYP2B6 has not been as fully characterized at the molecular level as other members of the human cytochrome P450 family. As more widely used in vitro probes for characterizing the involvement of this enzyme in the metabolism of xenobiotics have become available, the number of molecules identified as CYP2B6 substrates has increased. In this study we have analyzed the available kinetic data generated by multiple laboratories with human recombinant expressed CYP2B6 and along with calculated molecular properties derived from the ChemSpider database, we have determined the molecular features that appear to be important for CYP2B6 substrates. In addition we have applied 2D and 3D QSAR methods to generate predictive pharmacophore and 2D models. For 28 molecules with K(m) data, the molecular weight (mean +/- SD) is 253.78+/-74.03, ACD/logP is 2.68+/-1.51, LogD(pH 5.5) is 1.51+/-1.43, LogD(pH 7.4) is 2.02+/-1.25, hydrogen bond donor (HBD) count is 0.57 +/-0.57, hydrogen bond acceptor (HBA) count is 2.57+/-1.37, rotatable bonds is 3.50+/-2.71 and total polar surface area (TPSA) is 27.63+/-19.42. A second set of 15 molecules without K(m) data possessed similar mean molecular property values. These properties are comparable to those of a set of 21 molecules used in a previous pharmacophore modeling study (Ekins et al., J Pharmacol Exp Ther 288 (1), 21-29, 1999). Only the LogD and HBD values were statistically significantly different between these different datasets. We have shown that CYP2B6 substrates are generally small hydrophobic molecules that are frequently central nervous system active, which may be important for drug discovery research.

Analysis of the cryptophycin P450 epoxidase reveals substrate tolerance and cooperativity

Cryptophycins are potent anticancer agents isolated from Nostoc sp. ATCC 53789 and Nostoc sp. GSV 224. The most potent natural cryptophycin analogues retain a beta-epoxide at the C2'-C3' position of the molecule. A P450 epoxidase encoded by c rpE recently identified from the cryptophycin gene cluster was shown to install this key functional group into cryptophycin-4 (Cr-4) to produce cryptophycin-2 (Cr-2) in a regio- and stereospecific manner. Here we report a detailed characterization of the CrpE epoxidase using an engineered maltose binding protein (MBP)-CrpE fusion. The substrate tolerance of the CrpE polypeptide was investigated with a series of structurally related cryptophycin analogues generated by chemoenzymatic synthesis. The enzyme specifically installed a beta-epoxide between C2' and C3' of cyclic cryptophycin analogues. The kcat/Km values of the enzyme were determined to provide further insights into the P450 epoxidase catalytic efficiency affected by substrate structural variation. Finally, binding analysis revealed cooperativity of MBP-CrpE toward natural and unnatural desepoxy cryptophycin substrates.

Generation of a set of simple, interpretable ADMET rules of thumb

A set of simple, consistent structure-property guides have been determined from an analysis of a number of key ADMET assays run within GSK: solubility, permeability, bioavailability, volume of distribution, plasma protein binding, CNS penetration, brain tissue binding, P-gp efflux, hERG inhibition, and cytochrome P450 1A2/2C9/2C19/2D6/3A4 inhibition. The rules have been formulated using molecular properties that chemists intuitively know how to alter in a molecule, namely, molecular weight, logP, and ionization state. The rules supplement the more predictive black-box models available to us by clearly illustrating the key underlying trends, which are in line with reports in the literature. It is clear from the analyses reported herein that almost all ADMET parameters deteriorate with either increasing molecular weight, logP, or both, with ionization state playing either a beneficial or detrimental affect depending on the parameter in question. This study re-emphasizes the need to focus on a lower molecular weight and logP area of physicochemical property space to obtain improved ADMET parameters.

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.

Thermodynamics of ligand binding to P450 2B4 and P450eryF studied by isothermal titration calorimetry

Structural plasticity and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. As more mammalian P450 crystal structures have emerged, computational modeling has become a major tool to predict drug metabolism and interactions. There is a need for real solution thermodynamic data to support modeling and crystallographic observations. Using isothermal titration calorimetry (ITC) we successfully evaluated the conformational flexibility of P450 2B4 in binding imidazole inhibitors of different size and chemistry and dissected the stoichiometry and energetics of ligand binding allostery in P450eryF. Thermodynamic signatures obtained by ITC nicely correlated with structural and modeling results. Thus, ITC is a powerful tool to study structure-function relationships in P450s.

Crystal structure of cytochrome P450 MoxA from Nonomuraea recticatena (CYP105)

Cytochrome P450 MoxA (P450moxA) from a rare actinomycete Nonomuraea recticatena belongs to the CYP105 family and exhibits remarkably broad substrate specificity. Here, we demonstrate that P450moxA acts on several luciferin derivatives, which were originally identified as substrates of the human microsomal P450s. We also describe the crystal structure of P450moxA in substrate-free form. Structural comparison with various bacterial and human microsomal P450s reveals that the P450moxA structure is most closely related to that of the fungal nitric oxide reductase P450nor (CYP55A1). Final refined model of P450moxA comprises almost all the residues, including the "BC-loop" and "FG-loop" regions pivotal for substrate recognition, and the current structure thus defines a well-ordered substrate-binding pocket. Clear electron density map reveals that the MES molecule is bound to the substrate-binding site, and the sixth coordination position of the heme iron is not occupied by a water molecule, probably due to the presence of MES molecule in the vicinity of the heme. The unexpected binding of the MES molecule might reflect the ability of P450moxA to accommodate a broad range of structurally diverse compounds.

Unusual regioselectivity and active site topology of human cytochrome P450 2J2

The oxidation of six derivatives of terfenadone by recombinant human CYP2J2 (CYP = cytochrome P450) was studied by high-performance liquid chromatography coupled to mass spectrometry (MS) using tandem MS techniques and by 1H NMR spectroscopy. CYP2J2 exhibited a surprising regioselectivity in favor of the hydroxylation of the substrate terminal chain at the weakly reactive homobenzylic position. In contrast, hydroxylation of the same substrates by CYP3A4 mainly occurred on the most chemically reactive sites of the substrates (N-oxidation and benzylic hydroxylation). A 3D homology model of CYP2J2 was constructed using recently published structures of CYP2A6, CYP2B4, CYP2C8, CYP2C9, and CYP2D6 as templates. In contrast with other CYP2 structures, it revealed an active site cavity with a severely restricted access of substrates to the heme through a narrow hydrophobic channel. Dynamic docking of terfenadone derivatives in the CYP2J2 active site allowed one to interpret the unexpected regioselectivity of the hydroxylation of these substrates by CYP2J2, which is mainly based on this restricted access to the iron. The structural features that have been found to be important for recognition of substrates or inhibitors by CYP2J2 were also interpreted on the basis of CYP2J2-substrate interactions in this model.

Stereoselective hexobarbital 3'-hydroxylation by CYP2C19 expressed in yeast cells and the roles of amino acid residues at positions 300 and 476

We examined the enzymatic function of recombinant CYP2C19 in enantiomeric hexobarbital (HB) 3'-hydroxylation, and searched the roles of amino acid residues, such as Phe-100, Phe-114, Asp-293, Glu-300, and Phe-476 of CYP2C19 in the stereoselective HB 3'-hydroxylation, using a yeast cell expression system and site-directed mutagenesis method. CYP2C19 wild-type exerted substrate enantioselectivity of (R)-HB>>(S)-HB and metabolite diastereoselectivity of 3'(R)<3'(S) in 3'-hydroxylation of HB enantiomers. The substitution of Asp-293 by alanine failed to yield an observable peak at 450 nm in its reduced carbon monoxide-difference spectrum. CYP2C19-E300A and CYP2C19-E300V with alanine and valine, respectively, in place of Glu-300 exerted total HB 3'-hydroxylation activities of 45 and 108%, respectively, that of the wild-type. Interestingly, these two mutants showed substrate enantioselectivity of (R)-HB<(S)-HB, which is opposite to that of the wild-type, while metabolite diasteroselectivity remained unchanged. The replacement of Phe-476 by alanine increased total HB 3'-hydroxylation activity to approximately 3-fold that of the wild-type. Particularly, 3'(S)-OH-(S)-HB-forming activity elevated to 7-fold that of the wild-type, resulting in the reversal of the substrate enantioselectivity. In contrast, the substitution of phenylalanine at positions 100 and 114 by alanine did not produce a remarkable change in the total activity or the substrate enantioselectivity. These results indicate that Glu-300 and Phe-476 are important in stereoselective oxidation of HB enantiomers by CYP2C19.

Molecular mechanism of phase I and phase II drug-metabolizing enzymes: implications for detoxification

Enzymes that catalyze the biotransformation of drugs and xenobiotics are generally referred to as drug-metabolizing enzymes (DMEs). DMEs can be classified into two main groups: oxidative or conjugative. The NADPH-cytochrome P450 reductase (P450R)/cytochrome P450 (P450) electron transfer systems are oxidative enzymes that mediate phase I reactions, whereas the UDP-glucuronosyltransferases (UGTs) are conjugative enzymes that mediate phase II enzymes. Both enzyme systems are localized to the endoplasmic reticulum (ER) where a number of drugs are sequentially metabolized. DMEs, including P450s and UGTs, generally have a highly plastic active site that can accommodate a wide variety of substrates. The P450 and UGT genes constitute a supergene family, in which UGT proteins are encoded by distinct genes and a complex gene. Both the P450 and UGT genes have evolved to diversify their functions. This chapter reviews advances in understanding the structure and function of the P450R/P450 and UGT enzyme systems. In particular, the coordinate biotransformation of xenobiotics by phase I and II enzymes in the ER membrane is examined.

State-of-the-art tools for computational site of metabolism predictions: comparative analysis, mechanistical insights, and future applications

In drug design, it is crucial to have reliable information on how a chemical entity behaves in the presence of metabolizing enzymes. This requires substantial experimental efforts. Consequently, being able to predict the likely site/s of metabolism in any compound, synthesized or virtual, would be highly beneficial and time efficient. In this work, six different methodologies for predictions of the site of metabolism (SOM) have been compared and validated using structurally diverse data sets of drug-like molecules with well-established metabolic pattern in CYP3A4, CYP2C9, or both. Three of the methods predict the SOM based on the ligand's chemical structure, two additional methods use structural information of the enzymes, and the sixth method combines structure and ligand similarity and reactivity. The SOM is correctly predicted in 50 to 90% of the cases, depending on method and enzyme, which is an encouraging rate. We also discuss the underlying mechanisms of cytochrome P450 metabolism in the light of the results from this comparison.

Inhibition of recombinant cytochrome P450 isoforms 2D6 and 2C9 by diverse drug-like molecules

The affinities of a diverse set of 500 drug-like molecules to cytochrome P450 isoforms 2C9 and 2D6 were measured using recombinant expressed enzyme. The dose-response curve of each compound was fitted with a series of equations representing typical or various types of atypical kinetics. Atypical kinetics was identified where the Akaike Information Criterion, plus other criteria, suggested the kinetics was more complex than expected for a Michaelis-Menten model. Approximately 20% of the compounds were excluded due to poor solubility, and approximately 15% were excluded due to fluorescence interference. Of the remaining compounds, roughly half were observed to bind with an affinity of 200 microM or lower for each of the two isoforms. Atypical kinetics was observed in 18% of the compounds that bind to cytochrome 2C9, but less than 2% for 2D6. The resulting collection of competitive inhibitors and inactive compounds were analyzed for trends in binding affinity. For CYP2D6, a clear relationship between polar surface area and charge was observed, with the most potent inhibitors having a formal positive charge and a low percent polar surface area. For CYP2C9, no clear trend between activity and physicochemical properties could be seen for the group as a whole; however, certain classes of compounds have altered frequencies of activity and atypical kinetics.

In silico pharmacology for drug discovery: applications to targets and beyond

Computational (in silico) methods have been developed and widely applied to pharmacology hypothesis development and testing. These in silico methods include databases, quantitative structure-activity relationships, similarity searching, pharmacophores, homology models and other molecular modeling, machine learning, data mining, network analysis tools and data analysis tools that use a computer. Such methods have seen frequent use in the discovery and optimization of novel molecules with affinity to a target, the clarification of absorption, distribution, metabolism, excretion and toxicity properties as well as physicochemical characterization. The first part of this review discussed the methods that have been used for virtual ligand and target-based screening and profiling to predict biological activity. The aim of this second part of the review is to illustrate some of the varied applications of in silico methods for pharmacology in terms of the targets addressed. We will also discuss some of the advantages and disadvantages of in silico methods with respect to in vitro and in vivo methods for pharmacology research. Our conclusion is that the in silico pharmacology paradigm is ongoing and presents a rich array of opportunities that will assist in expediating the discovery of new targets, and ultimately lead to compounds with predicted biological activity for these novel targets.

Extending the capabilities of nature's most versatile catalysts: directed evolution of mammalian xenobiotic-metabolizing P450s

Cytochrome P450 enzymes are amongst the most versatile enzymatic catalysts known. The ability to introduce a single atom of oxygen into an organic substrate has led to the diversification and exploitation of these enzymes throughout nature. Nowhere is this versatility more apparent than in the mammalian liver, where P450 monooxygenases catalyze the metabolic clearance of innumerate drugs and other environmental chemicals. In addition to the aromatic and aliphatic hydroxylations, N- and O-dealkylations, and heteroatom oxidations that are common in drug metabolism, many more unusual reactions catalyzed by P450s have been discovered, including reductions, group transfers and other biotransformations not typically associated with monooxygenases. A research area that shows great potential for development over the next few decades is the directed evolution of P450s as biocatalysts. Mammalian xenobiotic-metabolizing P450s are especially well suited to such protein engineering due to their ability to interact with relatively wide ranges of substrates with marked differences in structure and physicochemical properties. Typical characteristics, such as the low turnover rates and poor coupling seen during the metabolism of xenobiotics, as well as the enzyme specificity towards particular substrates and reactions, can be improved by directed evolution. This mini-review will cover the fundamental enabling technologies required to successfully engineer P450s, examine the work done to date on the directed evolution of mammalian forms, and provide a perspective on what will be required for the successful implementation of engineered enzymes.

Structural insight into the altered substrate specificity of human cytochrome P450 2A6 mutants

Human P450 2A6 displays a small active site that is well adapted for the oxidation of small planar substrates. Mutagenesis of CYP2A6 resulted in an increased catalytic efficiency for indole biotransformation to pigments and conferred a capacity to oxidize substituted indoles (Wu, Z.-L., Podust, L.M., Guengerich, F.P. J. Biol. Chem. 49 (2005) 41090-41100.). Here, we describe the structural basis that underlies the altered metabolic profile of three mutant enzymes, P450 2A6 N297Q, L240C/N297Q and N297Q/I300V. The Asn297 substitution abolishes a potential hydrogen bonding interaction with substrates in the active site, and replaces a structural water molecule between the helix B'-C region and helix I while maintaining structural hydrogen bonding interactions. The structures of the P450 2A6 N297Q/L240C and N297Q/I300V mutants provide clues as to how the protein can adapt to fit the larger substituted indoles in the active site, and enable a comparison with other P450 family 2 enzymes for which the residue at the equivalent position was seen to function in isozyme specificity, structural integrity and protein flexibility.