Four Decades of Cytochrome P450 2B Research: From Protein Adducts to Protein Structures and Beyond

This article features selected findings from the senior author and colleagues dating back to 1978 and covering approximately three-fourths of the 60 years since the discovery of cytochrome P450. Considering the vast number of P450 enzymes in this amazing superfamily and their importance for so many fields of science and medicine, including drug design and development, drug therapy, environmental health, and biotechnology, a comprehensive review of even a single topic is daunting. To make a meaningful contribution to the 50th anniversary of Drug Metabolism and Disposition, we trace the development of the research in a single P450 laboratory through the eyes of seven individuals with different backgrounds, perspectives, and subsequent career trajectories. All co-authors are united in their fascination for the structural basis of mammalian P450 substrate and inhibitor selectivity and using such information to improve drug design and therapy. An underlying theme is how technological advances enable scientific discoveries that were impossible and even inconceivable to prior generations. The work performed spans the continuum from: 1) purification of P450 enzymes from animal tissues to purification of expressed human P450 enzymes and their site-directed mutants from bacteria; 2) inhibition, metabolism, and spectral studies to isothermal titration calorimetry, deuterium exchange mass spectrometry, and NMR; 3) homology models based on bacterial P450 X-ray crystal structures to rabbit and human P450 structures in complex with a wide variety of ligands. Our hope is that humanizing the scientific endeavor will encourage new generations of scientists to make fundamental new discoveries in the P450 field. SIGNIFICANCE STATEMENT: The manuscript summarizes four decades of work from Dr. James Halpert's laboratory, whose investigations have shaped the cytochrome P450 field, and provides insightful perspectives of the co-authors. This work will also inspire future drug metabolism scientists to make critical new discoveries in the cytochrome P450 field.

The Evolution of Drug Metabolism and Disposition: A Perspective From the Editors

This article was solicited to commemorate the 50th anniversary of Drug Metabolism and Disposition (DMD) and features perspectives from five former editors spanning the years 1994 to 2020. During that time frame the journal underwent significant changes in manuscript submission and processing as well as multiple generational changes in the composition of the editorial board and associate editors. A constant, however, has been the commitment to be the premier journal for publications of articles in the areas of drug metabolism, absorption, distribution, excretion, and pharmacokinetics. Advances in some of those areas during the past 3 decades have been monumental. Two cases in point involve cytochromes P450 and drug transporters. In 1994 rigorous characterization of human cytochrome P450 enzymes was in its infancy, there were no proven selective inhibitors, and the idea of solving a human P450 X-ray crystal structure was just a fantasy. Likewise, little was known about individual drug transporters. Today, detailed knowledge of individual human P450 enzymes and drug transporters is integral in drug design and drug discovery and in avoiding drug interactions. In the face of these huge advances in knowledge, each editor has been charged with maintaining the caliber and significance of the journal and its financial solvency while serving the needs of individual authors. We present 5 individual perspectives on the challenges and rewards of serving as DMD editor and hope that, by humanizing the job, we will encourage others to assume positions of responsibility in publication of society journals. SIGNIFICANCE STATEMENT: The 5 most recent former editors of DMD describe their experiences and perspectives on the position in the context of constantly changing scientific emphases, technology, and publishing practices. The article offers subscribers, authors, and future editors and editorial board members valuable insights into the inner workings of the journal.

Mammalian cytochrome P450 biodiversity: Physiological importance, function, and protein and genomic structures of cytochromes P4502B in multiple species of woodrats with different dietary preferences

The vast diversity of cytochrome P450 enzymes in mammals has been proposed to result in large measure from plant-animal warfare, whereby evolution of chemical defenses such as phenolics and terpenoids in plants led to duplication and divergence of P450 genes in herbivores. Over evolutionary time, natural selection is predicted to have produced P450s with high affinity and enhanced metabolism of substrates that are ingested regularly by herbivores. Interestingly, however, almost all knowledge of the interactions of mammalian P450 enzymes with substrates stems from studies of the metabolism of drugs and model compounds rather than studies on wild mammalian herbivores and their respective PSMs. A question of particular interest centers on the role of individual P450 enzymes in the ability of certain herbivores to specialize on plants that are lethal to most other species, including those from the same genus as the specialists. We tackled this intricate problem using a tractable natural system (herbivorous woodrats, genus Neotoma) focusing on comparisons of the specialist N. stephensi, the facultative specialist N. lepida, and the generalist N. albigula, and employing a cross-disciplinary approach involving ecology, biochemistry, pharmacology, structural biology, and genomics. Based on multiple findings suggesting the importance of CYP2B enzymes for ingestion of juniper and a major constituent, α-pinene, we characterized the structure, function and activity of several CYP2B enzymes in woodrats with different dietary habits. Results to date suggest that differences in CYP2B gene copy number may contribute to differential tolerance of PSMs among woodrat species, although additional work is warranted to firmly link gene copy number to juniper tolerance.

Strategies in herbivory by mammals revisited: The role of liver metabolism in a juniper specialist (Neotoma stephensi) and a generalist (Neotoma albigula)

Although herbivory is widespread among mammals, few species have adopted a strategy of dietary specialization. Feeding on a single plant species often exposes herbivores to high doses of plant secondary metabolites (PSMs), which may exceed the animal's detoxification capacities. Theory predicts that specialists will have unique detoxification mechanisms to process high levels of dietary toxins. To evaluate this hypothesis, we compared liver microsomal metabolism of a juniper specialist, Neotoma stephensi (diet >85% juniper), to a generalist, N. albigula (diet ≤30% juniper). Specifically, we quantified the concentration of a key detoxification enzyme, cytochrome P450 2B (CYP2B) in liver microsomes, and the metabolism of α-pinene, the most abundant terpene in the juniper species consumed by the specialist woodrat. In both species, a 30% juniper diet increased the total CYP2B concentration (2-3×) in microsomes and microsomal α-pinene metabolism rates (4-fold). In N. stephensi, higher levels of dietary juniper (60% and 100%) further induced CYP2B and increased metabolism rates of α-pinene. Although no species-specific differences in metabolism rates were observed at 30% dietary juniper, total microsomal CYP2B concentration was 1.7× higher in N. stephensi than in N. albigula (p < .01), suggesting N. stephensi produces one or more variant of CYP2B that is less efficient at processing α-pinene. In N. stephensi, the rates of α-pinene metabolism increased with dietary juniper and were positively correlated with CYP2B concentration. The ability of N. stephensi to elevate CYP2B concentration and rate of α-pinene metabolism with increasing levels of juniper in the diet may facilitate juniper specialization in this species.

So many roads traveled: A career in science and administration

I have traveled many roads during my career. After spending my first 19 years in Los Angeles, I became somewhat of an academic nomad, studying and/or working in six universities in the United States and three in Sweden. In chronological order, I have a B.A. in Scandinavian languages and literature from UCLA, a Ph.D. in biochemistry from Uppsala University, and an M.S. in toxicology from the Karolinska Institute. I have been in schools of natural science, pharmacy, and medicine and have worked in multiple basic science departments and one clinical department. I have served as a research-track and tenured faculty member, department chair, associate dean, and dean. My research has spanned toxinology, biochemistry, toxicology, and pharmacology. Through all the moves, I have gained much and lost some. For the past 40 years, my interest has been cytochrome P450 structure-function and structure-activity relationships. My lab has focused on CYP2B enzymes using X-ray crystallography, site-directed mutagenesis, deuterium-exchange MS, isothermal titration calorimetry, and computational methods in conjunction with a variety of functional assays. This comprehensive approach has enabled detailed understanding of the structural basis of the remarkable substrate promiscuity of CYP2B enzymes. We also have investigated the mechanisms of CYP3A4 allostery using biophysical and advanced spectroscopic techniques, and discovered a pivotal role of P450-P450 interactions and of multiple-ligand binding. A major goal of this article is to provide lessons that may be useful to scientists in the early and middle stages of their careers and those more senior scientists contemplating an administrative move.

Cytochrome P450 2B diversity in a dietary specialist—the red tree vole (Arborimus longicaudus)

Although herbivores rely on liver enzymes to biotransform plant secondary metabolites ingested in plant-based diets, only a few enzymes from a handful of species have been characterized at the genomic level. In this study, we examined cytochrome P450 2B (CYP2B) sequence diversity and gene copy number in a conifer specialist, the red tree vole (Arborimus longicaudus). We fed captive individuals exclusively Douglas-fir (Pseudotsuga menziesii) foliage, cloned and sequenced their liver CYP2B cDNA, and estimated CYP2B gene copy number. We identified 21 unique CYP2B nucleotide sequences, and 20 unique CYP2B amino acid sequences. Gene copy number of CYP2B was estimated at 7.7 copies per haploid genome. We compared red tree vole CYP2B with CYP2B sequences of a generalist, the prairie vole (Microtus ochrogaster), found in GenBank. Our study revealed that the CYP2B enzymes of red tree voles possess unique sequences compared to CYP2B enzymes of other herbivorous species. The unique combination of amino acid residues at key substrate recognition sites of CYP2B enzymes may underlie the ability of the red tree vole to specialize on a highly toxic diet of Douglas-fir.

Use of Phenoxyaniline Analogues To Generate Biochemical Insights into the Interactio n of Polybrominated Diphenyl Ether with CYP2B Enzymes

Human hepatic cytochromes P450 (CYP) are integral to xenobiotic metabolism. CYP2B6 is a major catalyst of biotransformation of environmental toxicants, including polybrominated diphenyl ethers (PBDEs). CYP2B substrates tend to contain halogen atoms, but the biochemical basis for this selectivity and for species specific determinants of metabolism has not been identified. Spectral binding titrations and inhibition studies were performed to investigate interactions of rat CYP2B1, rabbit CYP2B4, and CYP2B6 with a series of phenoxyaniline (POA) congeners that are analogues of PBDEs. For most congeners, there was a <3-fold difference between the spectral binding constants (K_S) and IC₅₀ values. In contrast, large discrepancies between these values were observed for POA and 3-chloro-4-phenoxyaniline. CYP2B1 was the enzyme most sensitive to POA congeners, so the Val-363 residue from that enzyme was introduced into CYP2B4 or CYP2B6. This substitution partially altered the protein-ligand interaction profiles to make them more similar to that of CYP2B1. Addition of cytochrome P450 oxidoreductase (POR) to titrations of CYP2B6 with POA or 2'4'5'TCPOA decreased the affinity of both ligands for the enzyme. Addition of cytochrome b₅ to a recombinant enzyme system containing POR and CYP2B6 increased the POA IC₅₀ value and decreased the 2'4'5'TCPOA IC₅₀ value. Overall, the inconsistency between K_S and IC₅₀ values for POA versus 2'4'5'TCPOA is largely due to the effects of redox partner binding. These results provide insight into the biochemical basis of binding of diphenyl ethers to human CYP2B6 and changes in CYP2B6-mediated metabolism that are dependent on POA congener and redox partner identity.

Halogen-π Interactions in the Cytochrome P450 Active Site: Structural Insights into Human CYP2B6 Substrate Selectivity

Numerous cytochrome P450 (CYP) 2B6 substrates including drugs and environmental chemicals are halogenated. To assess the role of halogen-π bonds in substrate selectivity and orientation in the active site, structures of four CYP2B6 monoterpenoid complexes were solved by X-ray crystallography. Bornyl bromide exhibited dual orientations in the active site with the predominant orientation revealing a bromine-π bond with the Phe108 side chain. Bornane demonstrated two orientations with equal occupancy; in both, the C2 atom that bears the bromine in bornyl bromide was displaced by more than 2.5 Å compared with the latter complex. The bromine in myrtenyl bromide π-bonded with Phe297 in CYP2B6, whereas the two major orientations in the active site mutant I114V exhibited bromine-π interactions with two additional residues, Phe108 and Phe115. Analysis of existing structures suggests that halogen-π interactions may be unique to the CYP2B enzymes within CYP family 2 but are also important for CYP3A enzymes.

Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

We used high hydrostatic pressure as a tool for exploring the conformational landscape of human cytochrome P450 3A4 (CYP3A4) by electron paramagnetic resonance and fluorescence spectroscopy. Site-directed incorporation of a luminescence resonance energy transfer donor-acceptor pair allowed us to identify a pressure-dependent equilibrium between two states of the enzyme, where an increase in pressure increased the spatial separation between the two distantly located fluorophores. This transition is characterized by volume change (ΔV°) and P1/2 values of -36.8 ± 5.0 mL/mol and 1.45 ± 0.33 kbar, respectively, which corresponds to a Keq° of 0.13 ± 0.06, so that only 15% of the enzyme adopts the pressure-promoted conformation at ambient pressure. This pressure-promoted displacement of the equilibrium is eliminated by the addition of testosterone, an allosteric activator. Using site-directed spin labeling, we demonstrated that the pressure- and testosterone-sensitive transition is also revealed by pressure-induced changes in the electron paramagnetic resonance spectra of a nitroxide side chain placed at position 85 or 409 of the enzyme. Furthermore, we observed a pressure-induced displacement of the emission maxima of a solvatochromic fluorophore (7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin) placed at the same positions, which suggests a relocation to a more polar environment. Taken together, the results reveal an effector-dependent conformational equilibrium between open and closed states of CYP3A4 that involves a pronounced change at the interface between the region of α-helices A/A' and the meander loop of the enzyme, where residues 85 and 409 are located. Our study demonstrates the high potential of pressure-perturbation strategies for studying protein conformational landscapes.

The draft genome sequence and annotation of the desert woodrat Neotoma lepida

We present the de novo draft genome sequence for a vertebrate mammalian herbivore, the desert woodrat (Neotoma lepida). This species is of ecological and evolutionary interest with respect to ingestion, microbial detoxification and hepatic metabolism of toxic plant secondary compounds from the highly toxic creosote bush (Larrea tridentata) and the juniper shrub (Juniperus monosperma). The draft genome sequence and annotation have been deposited at GenBank under the accession LZPO01000000.

Effect of detergent binding on cytochrome P450 2B4 structure as analyzed by X-ray crystallography and deuterium-exchange mass spectrometry

Multiple crystal structures of CYP2B4 have demonstrated the binding of the detergent 5-cyclohexyl-1-pentyl-β-D-maltoside (CYMAL-5) in a peripheral pocket located adjacent to the active site. To explore the consequences of detergent binding, X-ray crystal structures of the peripheral pocket mutant CYP2B4 F202W were solved in the presence of hexaethylene glycol monooctyl ether (C8E6) and CYMAL-5. The structure in the presence of CYMAL-5 illustrated a closed conformation indistinguishable from the previously solved wild-type. In contrast, the F202W structure in the presence of C8E6 revealed a detergent molecule that coordinated the heme-iron and extended to the protein surface through the substrate access channel 2f. Despite the overall structural similarity of these detergent complexes, remarkable differences were observed in the A, A', and H helices, the F-G cassette, the C-D and β4 loop region. Hydrogen-deuterium exchange mass spectrometry (DXMS) was employed to probe these differences and to test the effect of detergents in solution. The presence of either detergent increased the H/D exchange rate across the plastic regions, and the results obtained by DXMS in solution were consistent in general with the relevant structural snapshots. The study provides insight into effect of detergent binding and the interpretation of associated conformational dynamics of CYP2B4.

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

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

Structure-Function Analysis of Mammalian CYP2B Enzymes Using 7-Substituted Coumarin Derivatives as Probes: Utility of Crystal Structures and Molecular Modeling in Understanding Xenobiotic Metabolism

Crystal structures of CYP2B35 and CYP2B37 from the desert woodrat were solved in complex with 4-(4-chlorophenyl)imidazole (4-CPI). The closed conformation of CYP2B35 contained two molecules of 4-CPI within the active site, whereas the CYP2B37 structure demonstrated an open conformation with three 4-CPI molecules, one within the active site and the other two in the substrate access channel. To probe structure-function relationships of CYP2B35, CYP2B37, and the related CYP2B36, we tested the O-dealkylation of three series of related substrates-namely, 7-alkoxycoumarins, 7-alkoxy-4-(trifluoromethyl)coumarins, and 7-alkoxy-4-methylcoumarins-with a C1-C7 side chain. CYP2B35 showed the highest catalytic efficiency (kcat/KM) with 7-heptoxycoumarin as a substrate, followed by 7-hexoxycoumarin. In contrast, CYP2B37 showed the highest catalytic efficiency with 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), followed by 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC). CYP2B35 had no dealkylation activity with 7-MFC or 7-EFC. Furthermore, the new CYP2B-4-CPI-bound structures were used as templates for docking the 7-substituted coumarin derivatives, which revealed orientations consistent with the functional studies. In addition, the observation of multiple -Cl and -NH-π interactions of 4-CPI with the aromatic side chains in the CYP2B35 and CYP2B37 structures provides insight into the influence of such functional groups on CYP2B ligand binding affinity and specificity. To conclude, structural, computational, and functional analysis revealed striking differences between the active sites of CYP2B35 and CYP2B37 that will aid in the elucidation of new structure-activity relationships.

Using the Specialization Framework to Determine Degree of Dietary Specialization in a Herbivorous Woodrat

To be considered a dietary specialist, mammalian herbivores must consume large quantities of a plant species considered "difficult" with respect to nutrient or toxin content, and possess specialized adaptations to deal with plant defensive compounds or low nutritional content. Populations of Neotoma lepida in the Great Basin consume Juniperus osteosperma, a plant heavily defended by terpenes, but a detailed dietary analysis of this population is lacking. Therefore, we investigated the extent of dietary specialization in this species in comparison with the better-studied specialist species, N. stephensi. Microhistological analysis of feces from N. lepida revealed that greater than 90% of their diet in nature was comprised of juniper. In laboratory tolerance trials, N. lepida tolerated a diet of 80% J. osteosperma, similar to that observed for N. stephensi. There was no difference in the abilities of N. lepida and N. stephensi to metabolize hexobarbital, a proxy compound for terpene metabolism. In preference tests of native and non-native juniper species, N. lepida did not exhibit a preference for its native or co-occurring juniper, J. osteosperma, over the non-native species, J. monosperma, whereas N. stephensi preferred its native or co-occurring juniper J. monosperma over non-native J. osteosperma. Behavioral and habitat differences between these woodrat species lead to the categorization of N. stephensi as an obligate juniper specialist with a small range that overlaps that of its preferred food, J. monosperma, and N. lepida as a facultative juniper specialist with a large range, and only a portion of its distribution containing populations that feed extensively on J. osteosperma.

Functional importance of a peripheral pocket in mammalian cytochrome P450 2B enzymes

The functional importance of a peripheral pocket found in previously published X-ray crystal structures of CYP2B4 and CYP2B6 was probed using a biophysical approach. Introduction of tryptophan within the pocket of CYP2B4 at F202 or I241 leads to marked impairment of 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC) or 7-benzyloxyresorufin O-dealkylation efficiency; a similar substitution at F195, near the surface access to the pocket, does not affect these activities. The analogous CYP2B6 F202W mutant is inactive in the 7-EFC O-dealkylation assay. The stoichiometry of 7-EFC deethylation suggested that the decreased activity of F202W and I241W in CYP2B4 and lack of activity of F202W in CYP2B6 coincided with a sharp increase in the flux of reducing equivalents through the oxidase shunt to produce excess water. The results indicate that the chemical identity of residues within this peripheral pocket, but not at the mouth of the pocket, is important in substrate turnover and redox coupling, likely through effects on active site topology.

Structural and biophysical characterization of human cytochromes P450 2B6 and 2A6 bound to volatile hydrocarbons: analysis and comparison

X-ray crystal structures of complexes of cytochromes CYP2B6 and CYP2A6 with the monoterpene sabinene revealed two distinct binding modes in the active sites. In CYP2B6, sabinene positioned itself with the putative oxidation site located closer to the heme iron. In contrast, sabinene was found in an alternate conformation in the more compact CYP2A6, where the larger hydrophobic side chains resulted in a significantly reduced active-site cavity. Furthermore, results from isothermal titration calorimetry indicated a much more substantial contribution of favorable enthalpy to sabinene binding to CYP2B6 as opposed to CYP2A6, consistent with the previous observations with (+)-α-pinene. Structural analysis of CYP2B6 complexes with sabinene and the structurally similar (3)-carene and comparison with previously solved structures revealed how the movement of the F206 side chain influences the volume of the binding pocket. In addition, retrospective analysis of prior structures revealed that ligands containing -Cl and -NH functional groups adopted a distinct orientation in the CYP2B active site compared with other ligands. This binding mode may reflect the formation of Cl-π or NH-π bonds with aromatic rings in the active site, which serve as important contributors to protein-ligand binding affinity and specificity. Overall, the findings from multiple techniques illustrate how drugs metabolizing CYP2B6 and CYP2A6 handle a common hydrocarbon found in the environment. The study also provides insight into the role of specific functional groups of the ligand that may influence the binding to CYP2B6.

Interactions among cytochromes P450 in microsomal membranes: oligomerization of cytochromes P450 3A4, 3A5, and 2E1 and its functional consequences

The body of evidence of physiologically relevant P450-P450 interactions in microsomal membranes continues to grow. Here we probe oligomerization of human CYP3A4, CYP3A5, and CYP2E1 in microsomal membranes. Using a technique based on luminescence resonance energy transfer, we demonstrate that all three proteins are subject to a concentration-dependent equilibrium between the monomeric and oligomeric states. We also observed the formation of mixed oligomers in CYP3A4/CYP3A5, CYP3A4/CYP2E1, and CYP3A5/CYP2E1 pairs and demonstrated that the association of either CYP3A4 or CYP3A5 with CYP2E1 causes activation of the latter enzyme. Earlier we hypothesized that the intersubunit interface in CYP3A4 oligomers is similar to that observed in the crystallographic dimers of some microsomal drug-metabolizing cytochromes P450 (Davydov, D. R., Davydova, N. Y., Sineva, E. V., Kufareva, I., and Halpert, J. R. (2013) Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone. Biochem. J. 453, 219-230). Here we report the results of intermolecular cross-linking of CYP3A4 oligomers with thiol-reactive bifunctional reagents as well as the luminescence resonance energy transfer measurements of interprobe distances in the oligomers of labeled CYP3A4 single-cysteine mutants. The results provide compelling support for the physiological relevance of the dimer-specific peripheral ligand-binding site observed in certain CYP3A4 structures. According to our interpretation, these results reveal an important general mechanism that regulates the activity and substrate specificity of the cytochrome P450 ensemble through interactions between multiple P450 species. As a result of P450-P450 cross-talk, the catalytic properties of the cytochrome P450 ensemble cannot be predicted by simple summation of the properties of the individual P450 species.

Concurrent cooperativity and substrate inhibition in the epoxidation of carbamazepine by cytochrome P450 3A4 active site mutants inspired by molecular dynamics simulations

Cytochrome P450 3A4 (CYP3A4) is the major human P450 responsible for the metabolism of carbamazepine (CBZ). To explore the mechanisms of interactions of CYP3A4 with this anticonvulsive drug, we carried out multiple molecular dynamics (MD) simulations, starting with the complex of CYP3A4 manually docked with CBZ. On the basis of these simulations, we engineered CYP3A4 mutants I369F, I369L, A370V, and A370L, in which the productive binding orientation was expected to be stabilized, thus leading to increased turnover of CBZ to the 10,11-epoxide product. In addition, we generated CYP3A4 mutant S119A as a control construct with putative destabilization of the productive binding pose. Evaluation of the kinetics profiles of CBZ epoxidation demonstrate that CYP3A4-containing bacterial membranes (bactosomes) as well as purified CYP3A4 (wild-type and mutants I369L/F) exhibit substrate inhibition in reconstituted systems. In contrast, mutants S119A and A370V/L exhibit S-shaped profiles that are indicative of homotropic cooperativity. MD simulations with two to four CBZ molecules provide evidence that the substrate-binding pocket of CYP3A4 can accommodate more than one molecule of CBZ. Analysis of the kinetics profiles of CBZ metabolism with a model that combines the formalism of the Hill equation with an allowance for substrate inhibition demonstrates that the mechanism of interactions of CBZ with CYP3A4 involves multiple substrate-binding events (most likely three). Despite the retention of the multisite binding mechanism in the mutants, functional manifestations reveal an exquisite sensitivity to even minor structural changes in the binding pocket that are introduced by conservative substitutions such as I369F, I369L, and A370V.

The role of cytochrome P450 2B6 and 2B4 substrate access channel residues predicted based on crystal structures of the amlodipine complexes

Recent X-ray crystal structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with amlodipine showed two bound ligand molecules, one in the active site and one in the substrate access channel. Based on the X-ray crystal structures, we investigated the interactions of P450 2B4 and 2B6 with amlodipine using absorbance spectroscopy, and determined the steady-state kinetics of 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin oxidation by some access channel mutants to evaluate the functional role of these residues in substrate turnover. The results of absorbance titrations are consistent with a simple mechanism with two parallel binding events that result in the formation of the enzyme complex with two molecules of amlodipine. Using this model we were able to resolve two separate ligand-binding events, which are characterized by two distinct KD values in each enzyme. The access channel mutants R73K in P450 2B6 and R73K, V216W, L219W, and F220W in P450 2B4 showed a significant decrease in kcat/KM with the both substrates. Overall, the results suggest that P450 2B4 and 2B6 form an enzyme complex with two molecules of amlodipine in solution, and R73, V216, L219 and F220 in P450 2B4 may play an important role in substrate metabolism.

A large-scale allosteric transition in cytochrome P450 3A4 revealed by luminescence resonance energy transfer (LRET)

Effector-induced allosteric transitions in cytochrome P450 3A4 (CYP3A4) were investigated by luminescence resonance energy transfer (LRET) between two SH-reactive probes attached to various pairs of distantly located cysteine residues, namely the double-cysteine mutants CYP3A4(C64/C468), CYP3A4(C377/C468) and CYP3A4(C64/C121). Successive equimolar labeling of these proteins with the phosphorescent probe erythrosine iodoacetamide (donor) and the near-infrared fluorophore DY-731 maleimide (acceptor) allowed us to establish donor/acceptor pairs sensitive to conformational motions. The interactions of all three double-labeled mutants with the allosteric activators α-naphthoflavone and testosterone resulted in an increase in the distance between the probes. A similar effect was elicited by cholesterol. These changes in distance vary from 1.3 to 8.5 Å, depending on the position of the donor/acceptor pair and the nature of the effector. In contrast, the changes in the interprobe distance caused by such substrates as bromocriptine or 1-pyrenebutanol were only marginal. Our results provide a decisive support to the paradigm of allosteric modulation of CYP3A4 and indicate that the conformational transition caused by allosteric effectors increases the spatial separation between the beta-domain of the enzyme (bearing residues Cys₆₄ and Cys₃₇₇) and the alpha-domain, where Cys₁₂₁ and Cys₄₆₈ are located.

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

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

An in vivo assay for elucidating the importance of cytochromes P450 for the ability of a wild mammalian herbivore (Neotoma lepida) to consume toxic plants

An in vivo assay using the cytochrome P450 (P450) suicide inhibitor 1-aminobenzotriazole (ABT) and 24-h food intake was developed to determine the relative importance of P450s in two populations of Neotoma lepida with respect to biotransformation of plant secondary compounds in the animals' natural diets. The efficacy of ABT as a P450 inhibitor was first validated using hypnotic-state assays with and without pretreatment with ABT. Pretreatment with 100 mg/kg ABT by gavage increased hexobarbital sleep times 3-4-fold in both populations, showing effective inhibition of P450s in woodrats. Next, the Great Basin population was fed a terpene-rich juniper diet, and the Mojave population was fed a phenolic-rich creosote diet, with rabbit chow serving as the control diet in each group. Treatment with ABT inhibited food intake in the Great Basin population fed the juniper diet to a greater extent (35%) than the Great Basin population fed the control diet (19%) or the Mojave population fed the creosote diet (16%). The food intake of the Mojave population fed the control diet was not significantly inhibited by ABT. The findings suggest that the biotransformation of terpenes in juniper relies more heavily on P450s than that of phenolics in creosote. This assay provides an inexpensive and noninvasive method to explore the relative importance of P450s in the biotransformation strategies of wild mammalian herbivores.

Structural and thermodynamic basis of (+)-α-pinene binding to human cytochrome P450 2B6

Despite recent advances in atomic-level understanding of drug and inhibitor interactions with human cytochromes P450, the decades-old questions of chemical and structural determinants of hydrocarbon binding are still unanswered. (+)-α-Pinene is a monoterpene hydrocarbon that is widely distributed in the environment and a potent P450 2B inhibitor. Therefore, a combined biophysical and structural analysis of human P450 2B6 interactions with (+)-α-pinene was undertaken to elucidate the basis of the very high affinity binding. Binding of (+)-α-pinene to the P450 active site was demonstrated by a Type I spectral shift. Thermodynamics of ligand binding were explored using isothermal titration calorimetry and compared to those of P450 2A6, which is much less flexible than 2B6 based on comparison of multiple X-ray crystal structures. Consistent with expectation, entropy is the major driving force for hydrocarbon binding to P450 2A6, as evidenced by the calorimetric results. However, formation of the 2B6-(+)-α-pinene complex has a significant enthalpic component. A 2.0 Å resolution crystal structure of this enzyme-ligand complex reveals that the highly plastic 2B6 utilizes previously unrecognized rearrangements of protein motifs. The results indicate that the specific components of enthalpic contribution to ligand binding are closely tied to the degree of enzyme flexibility.

A structural snapshot of CYP2B4 in complex with paroxetine provides insights into ligand binding and clusters of conformational states

An X-ray crystal structure of CYP2B4 in complex with the drug paroxetine [(3S,4R)-3-[(2H-1,3-benzodioxol-5-yloxy)methyl]-4-(4-fluorophenyl)piperidine] was solved at 2.14 Å resolution. The structure revealed a conformation intermediate to that of the recently solved complex with amlodipine and that of the more compact complex with 4-(4-chlorophenyl)imidazole in terms of the placement of the F-G cassette. Moreover, comparison of the new structure with 15 previously solved structures of CYP2B4 revealed some new insights into the determinants of active-site size and shape. The 2B4-paroxetine structure is nearly superimposable on a previously solved closed structure in a ligand-free state. Despite the overall conformational similarity among multiple closed structures, the active-site cavity volume of the paroxetine complex is enlarged. Further analysis of the accessible space and binding pocket near the heme reveals a new subchamber that resulted from the movement of secondary structural elements and rearrangements of active-site side chains. Overall, the results from the comparison of all 16 structures of CYP2B4 demonstrate a cluster of protein conformations that were observed in the presence or absence of various ligands.

CYP261 enzymes from deep sea bacteria: a clue to conformational heterogeneity in cytochromes P450

We have explored the adaptation of the cytochromes P450 (P450) of deep-sea bacteria to high hydrostatic pressures. Strict conservation of the protein fold and functional importance of protein-bound water make P450 a unique subject for the studies of high-pressure adaptation. Earlier, we expressed and purified a fatty-acid binding P450 from the deep-sea bacteria Photobacterium profundum SS9 (CYP261C1). Here, we report purification and initial characterization of its mesophilic ortholog from the shallow-water P. profundum 3TCK (CYP261C2), as well as another piezophilic enzyme, CYP261D1, from deep-sea Moritella sp. PE36. Comparison of the three enzymes revealed a striking peculiarity of the piezophilic enzymes. Both CYP261C1 and CYP261D1 possess an apparent pressure-induced conformational toggle actuated at the pressures commensurate with the physiological pressure of habitation of the host bacteria. Furthermore, in contrast to CYP261C2, the piezophilic CYP261 enzymes may be chromatographically separated into two fractions with different properties, and different thermodynamic parameters of spin equilibrium in particular. According to our concept, the changes in the energy landscape that evolved in pressure-tolerant enzymes must stabilize the less-hydrated, closed conformers, which may be transient in the catalytic mechanisms of nonpiezophilic enzymes. The studies of enzymes of piezophiles should help unravel the mechanisms that control water access during the catalytic cycle.

Potent mechanism-based inactivation of cytochrome P450 2B4 by 9-ethynylphenanthrene: implications for allosteric modulation of cytochrome P450 catalysis

The mechanism-based inactivation of cytochrome P450 2B4 (CYP2B4) by 9-ethynylphenanthrene (9EP) has been investigated. The partition ratio and k(inact) are 0.2 and 0.25 min(-1), respectively. Intriguingly, the inactivation exhibits sigmoidal kinetics with a Hill coefficient of 2.5 and an S(50) of 4.5 μM indicative of homotropic cooperativity. Enzyme inactivation led to an increase in mass of the apo-CYP2B4 by 218 Da as determined by electrospray ionization liquid chromatography and mass spectrometry, consistent with covalent protein modification. The modified CYP2B4 was purified to homogeneity and its structure determined by X-ray crystallography. The structure showed that 9EP is covalently attached to Oγ of Thr 302 via an ester bond, which is consistent with the increased mass of the protein. The presence of the bulky phenanthrenyl ring resulted in inward rotations of Phe 297 and Phe 206, leading to a compact active site. Thus, binding of another molecule of 9EP in the active site is prohibited. However, results from the quenching of 9EP fluorescence by unmodified or 9EP-modified CYP2B4 revealed at least two binding sites with distinct affinities, with the low-affinity site being the catalytic site and the high-affinity site on the protein periphery. Computer-aided docking and molecular dynamics simulations with one or two ligands bound revealed that the high-affinity site is situated at the entrance of a substrate access channel surrounded by the F' helix, β1-β2 loop, and β4 loop and functions as an allosteric site to enhance the efficiency of activation of the acetylenic group of 9EP and subsequent covalent modification of Thr 302.

X-ray crystal structure of the cytochrome P450 2B4 active site mutant F297A in complex with clopidogrel: insights into compensatory rearrangements of the binding pocket

Prior X-ray crystal structures of cytochrome P450 2B4 revealed the pivotal role of rearrangement of the side chains of residues F206 and F297 in the active site in accommodating various inhibitors or substrates. To explore the role of these residues, 2B4 F206A and F297A were created by site-directed mutagenesis and characterized functionally. The structure of F297A with clopidogrel demonstrated the reorientation of the ligand such that the methyl ester group is oriented toward the heme, whereas the thiophene moiety now extends to the additional void in the F297A mutant. Most interestingly, movement of the I helix and several amino acid side chains within the active site was observed in apparent response to the altered binding orientation. Results of flexible docking using the 2B4 wild type or the F297A-virtual mutant positioned either the thiophene or chlorophenyl group closer to heme. However, docking of clopidogrel using the real F297A mutant or a virtual mutant with the I-helix re-positioned oriented clopidogrel preferentially with either the methyl ester or the chlorophenyl group closest to heme. The study provides insight into how the altered active site adapts to accommodate and interact with the substrate in a distinct orientation while maintaining the overall closed protein conformation.

Stereoselective formation and metabolism of 4-hydroxy-retinoic Acid enantiomers by cytochrome p450 enzymes

All-trans-retinoic acid (atRA), the major active metabolite of vitamin A, plays a role in many biological processes, including maintenance of epithelia, immunity, and fertility and regulation of apoptosis and cell differentiation. atRA is metabolized mainly by CYP26A1, but other P450 enzymes such as CYP2C8 and CYP3As also contribute to atRA 4-hydroxylation. Although the primary metabolite of atRA, 4-OH-RA, possesses a chiral center, the stereochemical course of atRA 4-hydroxylation has not been studied previously. (4S)- and (4R)-OH-RA enantiomers were synthesized and separated by chiral column HPLC. CYP26A1 was found to form predominantly (4S)-OH-RA. This stereoselectivity was rationalized via docking of atRA in the active site of a CYP26A1 homology model. The docked structure showed a well defined niche for atRA within the active site and a specific orientation of the β-ionone ring above the plane of the heme consistent with stereoselective abstraction of the hydrogen atom from the pro-(S)-position. In contrast to CYP26A1, CYP3A4 formed the 4-OH-RA enantiomers in a 1:1 ratio and CYP3A5 preferentially formed (4R)-OH-RA. Interestingly, CYP3A7 and CYP2C8 preferentially formed (4S)-OH-RA from atRA. Both (4S)- and (4R)-OH-RA were substrates of CYP26A1 but (4S)-OH-RA was cleared 3-fold faster than (4R)-OH-RA. In addition, 4-oxo-RA was formed from (4R)-OH-RA but not from (4S)-OH-RA by CYP26A1. Overall, these findings show that (4S)-OH-RA is preferred over (4R)-OH-RA by the enzymes regulating atRA homeostasis. The stereoselectivity observed in CYP26A1 function will aid in better understanding of the active site features of the enzyme and the disposition of biologically active retinoids.

Conformational adaptation of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 revealed upon binding multiple amlodipine molecules

Structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with two molecules of the calcium channel blocker amlodipine have been determined by X-ray crystallography. The presence of two drug molecules suggests clear substrate access channels in each P450. According to a previously established nomenclature, amlodipine molecules were trapped in access pathway 2f in P450 2B6 and in pathway 2a or 2f in P450 2B4. These pathways overlap for part of the length and then diverge as they extend toward the protein surface. A previously described solvent channel was also found in each enzyme. The results indicate that key residues located on the surface and at the entrance of the substrate access channels in each of these P450s may play a crucial role in guiding substrate entry. In addition, the region of P450 2B6 and 2B4 involving helices B', F, F', and G' and part of helix G is substantially more open in the amlodipine complexes than in the corresponding 4-(4-chlorophenyl)imidazole complexes. The increased active site volume observed results from the major retraction of helices F, F', and B' and the β4 sheet region located close to the binding cavity to accommodate amlodipine. These structures demonstrate novel insight into distinct conformational states not observed with previous P450 2B structures and provide clear evidence of the substrate access channels in two drug-metabolizing P450s. In addition, the structures exhibit the versatility that can be exploited via in silico studies with other P450 2B6 ligands as large as raloxifene and itraconazole.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Peripheral ligand-binding site in cytochrome P450 3A4 located with fluorescence resonance energy transfer (FRET)

The mechanisms of ligand binding and allostery in the major human drug-metabolizing enzyme cytochrome P450 3A4 (CYP3A4) were explored with fluorescence resonance energy transfer (FRET) using a laser dye, fluorol-7GA (F7GA), as a model substrate. Incorporation into the enzyme of a thiol-reactive FRET probe, pyrene iodoacetamide, allowed us to monitor the binding by FRET from the pyrene donor to the F7GA acceptor. Cooperativity of the interactions detected by FRET indicates that the enzyme possesses at least two F7GA-binding sites that have different FRET efficiencies and are therefore widely separated. To probe spatial localization of these sites, we studied FRET in a series of mutants bearing pyrene iodoacetamide at different positions, and we measured the distances from each of the sites to the donor. Our results demonstrate the presence of a high affinity binding site at the enzyme periphery. Analysis of the set of measured distances complemented with molecular modeling and docking allowed us to pinpoint the most probable peripheral site. It is located in the vicinity of residues 217-220, similar to the position of the progesterone molecule bound at the distal surface of the CYP3A4 in a prior x-ray crystal structure. Peripheral binding of F7GA causes a substantial spin shift and serves as a prerequisite for the binding in the active site. This is the first indication of functionally important ligand binding outside of the active site in cytochromes P450. The findings strongly suggest that the mechanisms of CYP3A4 cooperativity involve a conformational transition triggered by an allosteric ligand.

Investigation by site-directed mutagenesis of the role of cytochrome P450 2B4 non-active-site residues in protein-ligand interactions based on crystal structures of the ligand-bound enzyme

Residues located outside the active site of cytochromes P450 2B have exhibited importance in ligand binding, structural stability and drug metabolism. However, contributions of non-active-site residues to the plasticity of these enzymes are not known. Thus, a systematic investigation was undertaken of unique residue-residue interactions found in crystal structures of P450 2B4 in complex with 4-(4-chlorophenyl)imidazole (4-CPI), a closed conformation, or in complex with bifonazole, an expanded conformation. Nineteen mutants distributed over 11 sites were constructed, expressed in Escherichia coli and purified. Most mutants showed significantly decreased expression, especially in the case of interactions found in the 4-CPI structure. Six mutants (H172A, H172F, H172Q, L437A, E474D and E474Q) were chosen for detailed functional analysis. Among these, the K(s) of H172F for bifonazole was ∼ 20 times higher than for wild-type 2B4, and the K(s) of L437A for 4-CPI was ∼ 50 times higher than for wild-type, leading to significantly altered inhibitor selectivity. Enzyme function was tested with the substrates 7-ethoxy-4-(trifluoromethyl)coumarin, 7-methoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin (7-BR). H172F was inactive with all three substrates, and L437A did not turn over 7-BR. Furthermore, H172A, H172Q, E474D and E474Q showed large changes in k(cat)/K(M) for each of the three substrates, in some cases up to 50-fold. Concurrent molecular dynamics simulations yielded distances between some of the residues in these putative interaction pairs that are not consistent with contact. The results indicate that small changes in the protein scaffold lead to large differences in solution behavior and enzyme function.

The structural basis for homotropic and heterotropic cooperativity of midazolam metabolism by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) metabolizes a significant portion of clinically relevant drugs and often exhibits complex steady-state kinetics that can involve homotropic and heterotropic cooperativity between bound ligands. In previous studies, the hydroxylation of the sedative midazolam (MDZ) exhibited homotropic cooperativity via a decrease in the ratio of 1'-OH-MDZ to 4-OH-MDZ at higher drug concentrations. In this study, MDZ exhibited heterotropic cooperativity with the antiepileptic drug carbamazepine (CBZ) with characteristic decreases in the 1'-OH-MDZ to 4-OH-MDZ ratios. To unravel the structural basis of MDZ cooperativity, we probed MDZ and CBZ bound to CYP3A4 using longitudinal T(1) nuclear magnetic resonance (NMR) relaxation and molecular docking with AutoDock 4.2. The distances calculated from longitudinal T(1) NMR relaxation were used during simulated annealing to constrain the molecules to the substrate-free X-ray crystal structure of CYP3A4. These simulations revealed that either two MDZ molecules or an MDZ molecule and a CBZ molecule assume a stacked configuration within the CYP3A4 active site. In either case, the proton at position 4 of the MDZ molecule was closer to the heme than the protons of the 1'-CH(3) group. In contrast, molecular docking of a single molecule of MDZ revealed that the molecule was preferentially oriented with the 1'-CH(3) position closer to the heme than position 4. This study provides the first detailed molecular analysis of heterotropic and homotropic cooperativity of a human cytochrome P450 from an NMR-based model. Cooperativity of ligand binding through direct interaction between stacked molecules may represent a common motif for homotropic and heterotropic cooperativity.

Structures of cytochrome P450 2B6 bound to 4-benzylpyridine and 4-(4-nitrobenzyl)pyridine: insight into inhibitor binding and rearrangement of active site side chains

The biochemical, biophysical, and structural analysis of the cytochrome P450 2B subfamily of enzymes has provided a wealth of information regarding conformational plasticity and substrate recognition. The recent X-ray crystal structure of the drug-metabolizing P450 2B6 in complex with 4-(4-chlorophenyl)imidazole (4-CPI) yielded the first atomic view of this human enzyme. However, knowledge of the structural basis of P450 2B6 specificity and inhibition has remained limited. In this study, structures of P450 2B6 were determined in complex with the potent inhibitors 4-benzylpyridine (4-BP) and 4-(4-nitrobenzyl)pyridine (4-NBP). Comparison of the present structures with the previous P450 2B6-4-CPI complex showed that reorientation of side chains of the active site residue Phe206 on the F-helix and Phe297 on the I-helix was necessary to accommodate the inhibitors. However, P450 2B6 does not require any major side chain rearrangement to bind 4-NBP compared with 4-BP, and the enzyme provides no hydrogen-bonding partners for the polar nitro group of 4-NBP within the hydrophobic active site. In addition, on the basis of these new structures, substitution of residue 172 with histidine as observed in the single nucleotide polymorphism Q172H and in P450 2B4 may contribute to a hydrogen bonding network connecting the E- and I-helices, thereby stabilizing active site residues on the I-helix. These results provide insight into the role of active site side chains upon inhibitor binding and indicate that the recognition of the benzylpyridines in the closed conformation structure of P450 2B6 is based solely on hydrophobicity, size, and shape.

Institutional Profile: University of California San Diego Pharmacogenomics Education Program (PharmGenEd™): bridging the gap between science and practice

Clinical application of evidence-based pharmacogenomics information has the potential to help healthcare professionals provide safe and effective medication management to patients. However, there is a gap between the advances of pharmacogenomics discovery and the health professionals' knowledge regarding pharmacogenomics testing and therapeutic uses. Furthermore, pharmacogenomics education materials for healthcare professionals have not been readily available or accessible. Pharmacogenomics Education Program (PharmGenEd™) is an evidence-based pharmacogenomics education program developed at the University of California San Diego Skaggs School of Pharmacy and Pharmaceutical Sciences and the School of Medicine (CA, USA), with funding support from the Centers for Disease Control and Prevention. Program components include continuing education modules, train-the-trainer materials and shared curriculum modules based on therapeutic topics, and virtual communities with online resources.

Structural analysis of mammalian cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene: insight into partial enzymatic activity

A combined structural and computational analysis of rabbit cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene (tBPA) has yielded insight into how the enzyme retains partial activity. Since conjugation to tBPA modifies a highly conserved active site residue, the residual activity of tBPA-labeled 2B4 observed in previous studies was puzzling. Here we describe the first crystal structures of a modified mammalian P450, which show an oxygenated metabolite of tBPA conjugated to Thr 302 of helix I. These results are consistent with previous studies that identified Thr 302 as the site of conjugation. In each structure, the core of 2B4 remains unchanged, but the arrangement of plastic regions differs. This results in one structure that is compact and closed. In this conformation, tBPA points toward helix B', making a 31° angle with the heme plane. This conformation is in agreement with previously performed in silico experiments. However, dimerization of 2B4 in the other structure, which is caused by movement of the B/C loop and helices F through G, alters the position of tBPA. In this case, tBPA lies almost parallel to the heme plane due to the presence of helix F' of the opposite monomer entering the active site to stabilize the dimer. However, docking experiments using this open form show that tBPA is able to rotate upward to give testosterone and 7-ethoxy-4-trifluoromethylcoumarin access to the heme, which could explain the previously observed partial activity.