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

Multidisciplinary Insights into the Structure-Function Relationship of the CYP2B6 Active Site

Cytochrome P450 2B6 (CYP2B6) is a highly polymorphic human enzyme involved in the metabolism of many clinically relevant drugs, environmental toxins, and endogenous molecules with disparate structures. Over the last 20-plus years, in silico and in vitro studies of CYP2B6 using various ligands have provided foundational information regarding the substrate specificity and structure-function relationship of this enzyme. Approaches such as homology modeling, X-ray crystallography, molecular docking, and kinetic activity assays coupled with CYP2B6 mutagenesis have done much to characterize this originally neglected monooxygenase. However, a complete understanding of the structural details that make new chemical entities substrates of CYP2B6 is still lacking. Surprisingly little in vitro data has been obtained about the structure-function relationship of amino acids identified to be in the CYP2B6 active site. Since much attention has already been devoted to elucidating the function of CYP2B6 allelic variants, here we review the salient findings of in silico and in vitro studies of the CYP2B6 structure-function relationship with a deliberate focus on the active site. In addition to summarizing these complementary approaches to studying structure-function relationships, we note gaps/challenges in existing data such as the need for more CYP2B6 crystal structures, molecular docking results with various ligands, and data coupling CYP2B6 active site mutagenesis with kinetic parameter measurement under standard expression conditions. Harnessing in silico and in vitro techniques in tandem to understand the CYP2B6 structure-function relationship will likely offer further insights into CYP2B6-mediated metabolism. SIGNIFICANCE STATEMENT: The apparent importance of cytochrome P450 2B6 (CYP2B6) in the metabolism of various xenobiotics and endogenous molecules has grown since its discovery with many in silico and in vitro studies offering a partial description of its structure-function relationship. Determining the structure-function relationship of CYP2B6 is difficult but may be aided by thorough biochemical investigations of the CYP2B6 active site that provide a more complete pharmacological understanding of this important enzyme.

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.

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.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Structure of Cytochrome P450 2C9*2 in Complex with Losartan: Insights into the Effect of Genetic Polymorphism

The human CYP2C9 plays a crucial role in the metabolic clearance of a wide range of clinical therapeutics. The *2 allele is a prevalent genetic variation in CYP2C9 that is found in various populations. A marked reduction of catalytic activity toward many important drug substrates has been demonstrated by CYP2C9*2, which represents an amino acid variation at position 144 from arginine to cysteine. The crystal structure of CYP2C9*2 in complex with an antihypertensive drug losartan was solved using X-ray crystallography at 3.1-Å resolution. The Arg144Cys variation in the *2 complex disrupts the hydrogen-bonding interactions that were observed between the side chain of arginine and neighboring residues in the losartan complex of CYP2C9 and the wild-type (WT) ligand-free structure. The conformation of several secondary structural elements is affected, thereby altering the binding and orientation of drug and important amino acid side chains in the distal active site cavity. The new structure revealed distinct interactions of losartan in the compact active site of CYP2C9*2 and differed in occupancy at the other binding sites previously identified in the WT-losartan complex. Furthermore, the binding studies in solution using losartan illustrated lower activity of the CYP2C9*2 compared with the WT. Together, the findings yield valuable insights into the decreased hydroxylation activity of losartan in patients carrying CYP2C9*2 allele and provide a useful framework to investigate the effect of a single-nucleotide polymorphism that leads to altered metabolism of diverse drug substrates. SIGNIFICANCE STATEMENT: The *2 allele of the human drug-metabolizing enzyme CYP2C9 is found in different populations and results in significantly reduced activity toward various drug substrates. How the CYP2C9*2 variant induces altered drug metabolism is poorly understood given that the Arg144Cys variation is located far away from the active site. This work yield insight into the effect of distal variation using multitude of techniques that include X-ray crystallography, isothermal titration calorimetry, enzymatic characterization, and computational studies.

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.

Coumarins and P450s, Studies Reported to-Date

Cytochrome P450 enzymes (CYPs) are important phase I enzymes involved in the metabolism of endogenous and xenobiotic compounds mainly through mono-oxygenation reactions into more polar and easier to excrete species. In addition to their role in detoxification, they play important roles in the biosynthesis of endogenous compounds and the bioactivation of xenobiotics. Coumarins, phytochemicals abundant in food and commonly used in fragrances and cosmetics, have been shown to interact with P450 enzymes as substrates and/or inhibitors. In this review, these interactions and their significance in pharmacology and toxicology are discussed in detail.

Crystal Structure of CYP2B6 in Complex with an Efavirenz Analog

The over two dozen CYP2B structures of human, rabbit, and woodrat enzymes solved in the last decade have significantly enhanced our understanding of the structure-function relationships of drug metabolizing enzymes. More recently, an important role has emerged for halogen-π interactions in the CYP2B6 active site in substrate selectivity, explaining in part the preference for halogenated ligands as substrates. The mechanism by which such ligands interact with CYP2B enzymes involves conserved phenylalanine side chains, in particular F108, F115, or F297, in the active site, which form π bonds with halogens. To illustrate such halogen-π interactions using drugs that are major substrates of CYP2B6, we present here a crystal structure of CYP2B6 in complex with an analog of the widely used anti-HIV drug efavirenz, which contains a methyl group in place of the carbonyl oxygen. The chlorine of the efavirenz analog forms a π bond with the aromatic ring of F108, whereas the putative metabolism site on the distal end of the molecule is oriented towards the heme iron. The crystal structure showcases how CYP2B6 accommodates this important drug analog of considerable size in the active site by movement of various side chains without substantially increasing the active site volume. Furthermore, the CYP2B6-efavirenz analog complex provides a useful platform to investigate computationally as well as biophysically the effect of genetic polymorphisms on binding of the widely studied efavirenz.

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.

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.