Structures of cytochrome P450 2B4 complexed with the antiplatelet drugs ticlopidine and clopidogrel

Natural Products Inhibition of Cytochrome P450 2B6 Activity and Methadone Metabolism

Methadone is cleared predominately by hepatic cytochrome P450 (CYP) 2B6-catalyzed metabolism to inactive metabolites. CYP2B6 also catalyzes the metabolism of several other drugs. Methadone and CYP2B6 are susceptible to pharmacokinetic drug-drug interactions. Use of natural products such as herbals and other botanicals is substantial and growing, and concomitant use of prescription medicines and non-prescription herbals is common and may result in interactions, often precipitated by CYP inhibition. Little is known about herbal product effects on CYP2B6 activity, and CYP2B6-catalyzed methadone metabolism. We screened a family of natural product compounds used in traditional medicines, herbal teas, and synthetic analogs of compounds found in plants, including kavalactones, flavokavains, chalcones and gambogic acid, for inhibition of expressed CYP2B6 activity and specifically inhibition of CYP2B6-mediated methadone metabolism. An initial screen evaluated inhibition of CYP2B6-catalyzed 7-ethoxy-4-(trifluoromethyl) coumarin O-deethylation. Hits were further evaluated for inhibition of racemic methadone metabolism, including mechanism of inhibition and kinetic constants. In order of decreasing potency, the most effective inhibitors of methadone metabolism were dihydromethysticin (competitive, K i 0.074 µM), gambogic acid (noncompetitive, K i 6 µM), and 2,2'-dihydroxychalcone (noncompetitive, K i 16 µM). Molecular modeling of CYP2B6-methadone and inhibitor binding showed substrate and inhibitor binding position and orientation and their interactions with CYP2B6 residues. These results show that CYP2B6 and CYP2B6-catalyzed methadone metabolism are inhibited by certain natural products, at concentrations which may be clinically relevant. SIGNIFICANCE STATEMENT: This investigation identified several natural product constituents which inhibit in vitro human recombinant CYP2B6 and CYP2B6-catalyzed N-demethylation of the opioid methadone. The most potent inhibitors (K i) were dihydromethysticin (0.074 µM), gambogic acid (6 µM) and 2,2'-dihydroxychalcone (16 µM). Molecular modeling of ligand interactions with CYP2B6 found that dihydromethysticin and 2,2'-dihydroxychalcone bound at the active site, while gambogic acid interacted with an allosteric site on the CYP2B6 surface. Natural product constituents may inhibit CYP2B6 and methadone metabolism at clinically relevant concentrations.

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.

Inhibition of CYP2C8 by Acyl Glucuronides of Gemfibrozil and Clopidogrel: Pharmacological Significance, Progress and Challenges

The lipid-regulating drug gemfibrozil is a useful medication for reducing high cholesterol and triglycerides in the blood. In addition to oxidation, it undergoes extensive glucuronidation to produce gemfibrozil acyl glucuronide, which is a known mechanism-based inactivator of cytochrome P450 (CYP) 2C8. Such selective and time-dependent inhibition results in clinically important drug–drug interactions (DDI) with the drugs metabolized by CYP2C8. Similarly, the acyl glucuronide of clopidogrel, a widely used antiplatelet agent, is a potent time-dependent inhibitor of CYP2C8 that demonstrated significant DDI with the substrates of CYP2C8. Current progress in atomic-level understanding mostly involves studying how different drugs bind and undergo oxidation in the active site of CYPs. It is not clear how an acyl glucuronide metabolite of the drug gemfibrozil or clopidogrel interacts in the active site of CYP2C8 and selectively inhibit the enzyme. This mini-review summarizes the current knowledge on some of the important clinical DDI caused by gemfibrozil and clopidogrel due to the inhibition of CYP2C8 by acyl glucuronide metabolites of these drugs. Importantly, it examines recent developments and potential applications of structural biology tools to elucidate the binding and orientation of gemfibrozil acyl glucuronide and clopidogrel acyl glucuronide in the active site near heme that contributes to the inhibition and inactivation of CYP2C8.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

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.

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.

Capturing the photo-signaling state of a photoreceptor in a steady-state fashion by binding a transition metal complex

Binding a small molecule to proteins causes conformational changes, but often to a limited extent. Here, we demonstrate that the interaction of a CO-releasing molecule (CORM3) with a photoreceptor photoactive yellow protein (PYP) drives large structural changes in the latter. The interaction of CORM3 and a mutant of PYP, Met100Ala, not only trigger the isomerization of its chromophore, p-coumaric acid, from its anionic trans configuration to a protonated cis configuration, but also increases the content of β-sheet at the cost of α-helix and random coil in the secondary structure of the protein. The CORM3 derived Met100Ala is found to highly resemble the signaling state, which is one of the key photo-intermediates of this photoactive protein, in both protein local conformation and chromophore configuration. The organometallic reagents hold promise as protein engineering tools. This work highlights a novel approach to structurally accessing short lived intermediates of proteins in a steady-state fashion.

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.

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.

Bifunctional CYP81AA proteins catalyse identical hydroxylations but alternative regioselective phenol couplings in plant xanthone biosynthesis

Xanthones are natural products present in plants and microorganisms. In plants, their biosynthesis starts with regioselective cyclization of 2,3′,4,6-tetrahydroxybenzophenone to either 1,3,5- or 1,3,7-trihydroxyxanthones, catalysed by cytochrome P450 (CYP) enzymes. Here we isolate and express CYP81AA-coding sequences from Hypericum calycinum and H. perforatum in yeast. Microsomes catalyse two consecutive reactions, that is, 3′-hydroxylation of 2,4,6-trihydroxybenzophenone and C–O phenol coupling of the resulting 2,3′,4,6-tetrahydroxybenzophenone. Relative to the inserted 3′-hydroxyl, the orthologues Hc/HpCYP81AA1 cyclize via the para position to form 1,3,7-trihydroxyxanthone, whereas the paralogue HpCYP81AA2 directs cyclization to the ortho position, yielding the isomeric 1,3,5-trihydroxyxanthone. Homology modelling and reciprocal mutagenesis reveal the impact of S375, L378 and A483 on controlling the regioselectivity of HpCYP81AA2, which is converted into HpCYP81AA1 by sextuple mutation. However, the reciprocal mutations in HpCYP81AA1 barely affect its regiospecificity. Product docking rationalizes the alternative C–O phenol coupling reactions. Our results help understand the machinery of bifunctional CYPs.

Xanthones are pharmacologically and biosynthetically intriguing compounds. Here, the authors identify two cytochrome P450 enzymes, which hydroxylate and cyclize the benzophenone precursor to either 1,3,7- or 1,3,5-trihydroxyxanthones, and pinpoint residues that determine the alternative regioselectivities.

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.

A Comparative Study of Molecular Structure, pKa, Lipophilicity, Solubility, Absorption and Polar Surface Area of Some Antiplatelet Drugs

Theoretical chemistry methods have been used to study the molecular properties of antiplatelet agents (ticlopidine, clopidogrel, prasugrel, elinogrel, ticagrelor and cangrelor) and several thiol-containing active metabolites. The geometries and energies of most stable conformers of these drugs have been computed at the Becke3LYP/6-311++G(d,p) level of density functional theory. Computed dissociation constants show that the active metabolites of prodrugs (ticlopidine, clopidogrel and prasugrel) and drugs elinogrel and cangrelor are completely ionized at pH 7.4. Both ticagrelor and its active metabolite are present at pH = 7.4 in neutral undissociated form. The thienopyridine prodrugs ticlopidine, clopidogrel and prasugrel are lipophilic and insoluble in water. Their lipophilicity is very high (about 2.5–3.5 logP values). The polar surface area, with regard to the structurally-heterogeneous character of these antiplatelet drugs, is from very large interval of values of 3–255 Ų. Thienopyridine prodrugs, like ticlopidine, clopidogrel and prasugrel, with the lowest polar surface area (PSA) values, exhibit the largest absorption. A high value of polar surface area (PSA) of cangrelor (255 Ų) results in substantial worsening of the absorption in comparison with thienopyridine drugs.

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.

Insights into the role of substrates on the interaction between cytochrome b5 and cytochrome P450 2B4 by NMR

Mammalian cytochrome b5 (cyt b5) is a membrane-bound protein capable of donating an electron to cytochrome P450 (P450) in the P450 catalytic cycle. The interaction between cyt b5 and P450 has been reported to be affected by the substrates of P450; however, the mechanism of substrate modulation on the cyt b5-P450 complex formation is still unknown. In this study, the complexes between full-length rabbit cyt b5 and full-length substrate-free/substrate-bound cytochrome P450 2B4 (CYP2B4) are investigated using NMR techniques. Our findings reveal that the population of complexes is ionic strength dependent, implying the importance of electrostatic interactions in the complex formation process. The observation that the cyt b5-substrate-bound CYP2B4 complex shows a weaker dependence on ionic strength than the cyt b5-substrate-free CYP2B4 complex suggests the presence of a larger fraction of steoreospecific complexes when CYP2B4 is substrate-bound. These results suggest that a CYP2B4 substrate likely promotes specific interactions between cyt b5 and CYP2B4. Residues D65, V66, T70, D71 and A72 are found to be involved in specific interactions between the two proteins due to their weak response to ionic strength change. These findings provide insights into the mechanism underlying substrate modulation on the cyt b5-P450 complexation process.

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.

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

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.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.

Trendspotting in the Protein Data Bank

The Protein Data Bank (PDB) was established in 1971 as a repository for the three dimensional structures of biological macromolecules. Since then, more than 85000 biological macromolecule structures have been determined and made available in the PDB archive. Through analysis of the corpus of data, it is possible to identify trends that can be used to inform us abou the future of structural biology and to plan the best ways to improve the management of the ever-growing amount of PDB data.

Steroid-based facial amphiphiles for stabilization and crystallization of membrane proteins

Amphiphile selection is a critical step for structural studies of membrane proteins (MPs). We have developed a family of steroid-based facial amphiphiles (FAs) that are structurally distinct from conventional detergents and previously developed FAs. The unique FAs stabilize MPs and form relatively small protein-detergent complexes (PDCs), a property considered favorable for MP crystallization. We attempted to crystallize several MPs belonging to different protein families, including the human gap junction channel protein connexin 26, the ATP binding cassette transporter MsbA, the seven-transmembrane G protein-coupled receptor-like bacteriorhodopsin, and cytochrome P450s (peripheral MPs). Using FAs alone or mixed with other detergents or lipids, we obtained 3D crystals of the above proteins suitable for X-ray crystallographic analysis. The fact that FAs enhance MP crystallizability compared with traditional detergents can be attributed to several properties, including increased protein stability, formation of small PDCs, decreased PDC surface flexibility, and potential to mediate crystal lattice contacts.

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.

Jahn-Teller dynamics in a series of high-symmetry Co(II) chelates determine paramagnetic relaxation enhancements

NMR paramagnetic relaxation enhancements (PREs) of a series of structurally characterized, trigonal bis-trispyrazolylborate (Tp) chelates of high-spin Co(II), spanning 100-850 MHz in field, are reported. Prior knowledge of the metal-nucleus distances allows numerical extraction of position-dependent electron spin relaxation rates (τ(c)(-1)) from direct measurement of the individual PREs of the four symmetry distinct protons in Co(Tp)(2), using available closed-form expressions. The data for this electronically complex system where spin-orbit coupling defines the ground state electronic structure are analyzed in terms of the Solomon-Bloembergen-Morgan (SBM) relations, as well as available zero-field splitting limit theories. A simple angular correction is shown to be sufficient to reconcile the individual τ(c)(T) data for the four classes of protons. The data identify a previously unrecognized dynamic Jahn-Teller effect in these historically important complexes, with a barrier of ~230 cm(-1), pointing to a level of dynamics in trispyrazolylborate chemistry that has not been described before, and further show that it is the Jahn-Teller that is responsible for the PREs in fluid solution. A field-dependent component is also identified for the two protons nearest g(//), which is suggested to arise due to Zeeman mixing of excited state character into the ground level.

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.

Preferred binding orientations of phenacetin in CYP1A1 and CYP1A2 are associated with isoform-selective metabolism

Human cytochromes P450 1A1 and 1A2 play important roles in drug metabolism and chemical carcinogenesis. Although these two enzymes share high sequence identity, they display different substrate specificities and inhibitor susceptibilities. In the present studies, we investigated the structural basis for these differences with phenacetin as a probe using a number of complementary approaches, such as enzyme kinetics, stoichiometric assays, NMR, and molecular modeling. Kinetic and stoichiometric analyses revealed that substrate specificity (k(cat)/K(m)) of CYP1A2 was approximately 18-fold greater than that of CYP1A1, as expected. Moreover, despite higher H₂O₂ production, the coupling efficiency of reducing equivalents to acetaminophen formation in CYP1A2 was tighter than that in CYP1A1. CYP1A1, in contrast to CYP1A2, displayed much higher uncoupling, producing more water. The subsequent NMR longitudinal (T₁) relaxation studies with the substrate phenacetin and its product acetaminophen showed that both compounds displayed similar binding orientations within the active site of CYP1A1 and CYP1A2. However, the distance between the OCH₂ protons of the ethoxy group (site of phenacetin O-deethylation) and the heme iron was 1.5 Å shorter in CYP1A2 than in CYP1A1. The NMR findings are thus consistent with our kinetic and stoichiometric results, providing a likely molecular basis for more efficient metabolism of phenacetin by CYP1A2.

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.

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.

New amphiphiles for membrane protein structural biology

A challenging requirement for X-ray crystallography or NMR structure determination of membrane proteins (MPs), in contrast to soluble proteins, is the necessary use of amphiphiles to mimic the hydrophobic environment of membranes. A number of new detergents, lipids and non-detergent-like amphiphiles have been developed that stabilize MPs, and these have contributed to increased success in MP structural determinations in recent years. Despite some successes, currently available reagents are inadequate and there remains a pressing need for new amphiphiles. Literature examples and some new developments are selected here as a framework for discussing desirable properties of new amphiphiles for MP structural biology.

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.

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.

Structure and function of cytochromes P450 2B: from mechanism-based inactivators to X-ray crystal structures and back

This article reviews work from the author dating back to 1978 and focuses on the structural basis of cytochrome P450 (P450) function using available contemporary techniques. Early studies used mechanism-based inactivators that bound to the protein moiety of hepatic P450s to try to localize the active site. Subsequent studies used cDNA cloning, heterologous expression, site-directed mutagenesis, and homology modeling based on multiple bacterial P450 X-ray crystal structures to predict the active sites of CYP2B enzymes with considerable accuracy. Breakthroughs in engineering and expression of mammalian P450s enabled us to determine our first X-ray crystal structure of ligand-free rabbit CYP2B4. To date, we have solved 11 CYP2B4 and three human CYP2B6 structures, which represent four significantly different conformations. The plasticity of CYP2B4 has been confirmed by deuterium exchange mass spectrometry and is substantiated by molecular dynamics simulations. In addition to major movement of secondary structure elements, more subtle reorientation of active site side chains, especially Phe206, Phe297, and Glu301, contributes to the ability of CYP2B enzymes to bind various ligands. Isothermal titration calorimetry has proven to be a useful tool for studying the thermodynamics of ligand binding to CYP2B4 and CYP2B6, and NMR has enabled study of ligand binding orientation in solution as an adjunct to X-ray crystallography. A major challenge remains to harness the power of the various approaches to facilitate prediction of CYP2B specificity and inhibition.

NMR-derived models of amidopyrine and its metabolites in complexes with rabbit cytochrome P450 2B4 reveal a structural mechanism of sequential N-dealkylation

To understand the molecular basis of sequential N-dealkylation by cytochrome P450 2B enzymes, we studied the binding of amidopyrine (AP) as well as the metabolites of this reaction, desmethylamidopyrine (DMAP) and aminoantipyrine (AAP), using the X-ray crystal structure of rabbit P450 2B4 and two nuclear magnetic resonance (NMR) techniques: saturation transfer difference (STD) spectroscopy and longitudinal (T(1)) relaxation NMR. Results of STD NMR of AP and its metabolites bound to P450 2B4 were similar, suggesting that they occupy similar niches within the enzyme's active site. The model-dependent relaxation rates (R(M)) determined from T(1) relaxation NMR of AP and DMAP suggest that the N-linked methyl is closest to the heme. To determine the orientation(s) of AP and its metabolites within the P450 2B4 active site, we used distances calculated from the relaxation rates to constrain the metabolites to the X-ray crystal structure of P450 2B4. Simulated annealing of the complex revealed that the metabolites do indeed occupy similar hydrophobic pockets within the active site, while the N-linked methyls are free to rotate between two binding modes. From these bound structures, a model of N-demethylation in which the N-linked methyl functional groups rotate between catalytic and noncatalytic positions was developed. This study is the first to provide a structural model of a drug and its metabolites complexed to a cytochrome P450 based on NMR and to provide a structural mechanism for how a drug can undergo sequential oxidations without unbinding. The rotation of the amide functional group might represent a common structural mechanism for N-dealkylation reactions for other drugs such as the local anesthetic lidocaine.

Comparison of in vitro metabolism of ticlopidine by human cytochrome P450 2B6 and rabbit cytochrome P450 2B4

A recent X-ray crystal structure of a rabbit cytochrome P450 2B4 (CYP2B4)-ticlopidine complex indicated that the compound could be modeled with either the thiophene or chlorophenyl group oriented toward the heme prosthetic group. Subsequent NMR relaxation and molecular docking studies suggested that orientation with the chlorophenyl ring closer to the heme was the preferred one. To evaluate the predictive value of these findings, the oxidation of ticlopidine by reconstituted CYP2B4 was studied and compared with CYP2B6, in which the thiophene portion of the molecule likely orients toward the heme. In vitro incubation of ticlopidine with both enzymes yielded the same set of metabolites: 7-hydroxyticlopidine (M1), 2-oxoticlopidine (M2), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M3), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M4), ticlopidine N-oxide (M5), and ticlopidine S-oxide dimer, a dimerization product of ticlopidine S-oxide (M6). The rates of metabolite formation deviated markedly from linearity with time, consistent with the known inactivation of CYP2B6 by ticlopidine. Fitting to a first-order equation yielded similar rate constants (k(obs)) for both enzymes. However, the amplitude (R(max)) of M1 and M6 formation was 4 to 5 times higher for CYP2B6 than CYP2B4, indicating a greater residence time of ticlopidine with its thiophene ring closer to heme in CYP2B6. In contrast, CYP2B4 formed M4 and M5 in more abundance than CYP2B6, indicating an alternate orientation. Overall, the results suggest that the preferential orientation of ticlopidine in the active site of CYP2B4 predicted by X-ray crystallography and NMR studies is unproductive and that ticlopidine likely reorients within CYP2B4 to a more productive mode.