Molecular modeling of cytochrome P450 1A1: enzyme-substrate interactions and substrate binding affinities

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives. In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyresorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.     

Regiospecificity of Human Cytochrome P450 1A1-Mediated Oxidations: The Role of Steric Effects

Abstract

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,85- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine. Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.     

Regiospecificity of human cytochrome P450 1A1-mediated oxidations: the role of steric effects

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,8S- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine, Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.     

Functional analysis and crystallographic structure of clotrimazole bound OleP,a cytochrome P450 epoxidase from Streptomyces antibioticus involved in oleandomycin biosynthesis

BackgroundOleP is a cyt P450 from Streptomyces antibioticus carrying out epoxigenation of the antibiotic oleandomycin during its biosynthesis. The timing of its reaction has not been fully clarified, doubts remain regarding its substrate and catalytic mechanism.MethodsThe crystal structure of OleP in complex with clotrimazole, an inhibitor of P450s used in therapy, was solved and the complex formation dynamics was characterized by equilibrium and kinetic binding studies and compared to ketoconazole, another azole differing for the N1-substituent.ResultsClotrimazole coordinates the heme and occupies the active site. Most of the residues interacting with clotrimazole are conserved and involved in substrate binding in MycG, the P450 epoxigenase with the highest homology with OleP. Kinetic characterization of inhibitor binding revealed OleP to follow a simple bimolecular reaction, without detectable intermediates.ConclusionsClotrimazole-bound OleP adopts an open form, held by a π-π stacking chain that fastens helices F and G and the FG loop. Affinity is affected by the interactions of the N1 substituent within the active site, given the one order of magnitude difference of the off-rate constants between clotrimazole and ketoconazole. Based on structural similarities with MycG, we propose a binding mode for both oleandomycin intermediates, that are the candidate substrates of OleP.General significanceAmong P450 epoxigenases OleP is the only one that introduces an epoxide on a non-activated C–C bond. The data here presented are necessary to understand the rare chemistry carried out by OleP, to engineer it and to design more selective and potent P450-targeted drugs.     

The role of cytochrome P450 BM3 phenylalanine-87 and threonine-268 in binding organic hydroperoxides

BackgroundCytochrome P450 (P450) BM3, from Bacillus megaterium, catalyzes a wide range of chemical reactions and is routinely used as a model system to study mammalian P450 reactions and structure.MethodsThe metabolism of 2,6-di-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadienone (BHTOOH) and 2-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadien-1-one (BMPOOH) was examined with P450 BM3 and with the conserved T268 and F87 residues mutated to investigate their effects on organic hydroperoxide metabolism. To determine the effects of the mutations on the active site volume and architecture, the X-ray crystal structure of the F87A/T268A P450 BM3 heme domain (BMP) was determined and compared to previous structures. To investigate the interactions of the substrates with the F87 and T268 residues, BHTOOH and BMPOOH were docked into the BMP X-ray crystal structures.ResultsLower metabolism of BHTOOH and BMPOOH was observed in the WT P450 BM3 and the T268A P450 BM3 mutant than in the F87A and F87A/T268A P450 BM3 mutants. Large differences were found in the F–G loop regions and active site cavity volumes for the F87A mutated structures.ConclusionsAnalysis of the metabolism, X-ray crystal structures, and molecular docking simulations suggests that P450 BM3 activity toward BHTOOH and BMPOOH is mediated through substrate recognition by T268 and F87, and the active site cavity volume. Based on this information, a simplified representation is presented with the relative orientation of organic hydroperoxides in the P450 BM3 active site.General significanceThe metabolism results and structural analysis of this model P450 allowed us to rationalize the structural factors that influence organic hydroperoxide metabolism.     

Molecular Modeling of Cytochrome P450 2B1: Mode of Membrane Insertion and Substrate Specificity

A molecular model of a mammalian membrane-bound cytochrome P450, rat P450 2B1, was constructed in order to elucidate its mode of attachment to the endoplasmic reticulum and the structural basis of substrate specificity. The model was primarily derived from the structure of P450BM-3, which as a class II P450 is the most functionally similar P450 of known structure. However, model development was also guided by the conserved core regions of P450cam and P450terp. To optimally align the P450 2B1 and P450BM-3 sequences, multiple alignment was performed using sequences of five P450s in the II family, followed by minor adjustments on the basis of secondary structure predictions. The resulting P450 2B1 homology model structure was refined by molecular dynamics heating, equilibration, simulation, and energy minimization. The model suggests that the F–G loop serves as both a hydrophobic membrane anchor and entrance channel for hydrophobic substrates from the membrane to the P450 active site. To assess the mode of substrate binding, benzphetamine, testosterone, and benzo[a]pyrene were docked into the active site. The hydrophobic substrate-binding pocket is consistent with the preferences of this P450 toward hydrophobic substrates, while the presence of an acidic Glu-105 in this pocket is consistent with the preference of this P450 for the cationic substrate benzphetamine. This model is thus consistent with several known experimental properties of this P450, such as membrane attachment and substrate selectivity.     

Role of protein and substrate dynamics in catalysis by Pseudomonas putida cytochrome P450cam

The role of protein structural flexibility and substrate dynamics in catalysis by cytochrome P450 enzymes is an area of current interest. We have addressed these in cytochrome P450(cam) (P450(cam)) and its Y96A mutant with camphor and its related compounds using fluorescence spectroscopy. Previously [Prasad et al. (2000) FEBS Lett. 477, 157-160], we provided experimental support to dynamic fluctuations in P450(cam), and substrate access into the active site region via the channel next to the flexible F-G helix-loop-helix segment. In the investigation described here, we show that the dynamic fluctuations in the enzyme are substrate dependent as reflected by tryptophan fluorescence quenching experiments. The orientation of tryptophan relative to heme (kappa(2)) for W42 obtained from time-resolved tryptophan fluorescence measurements show variation with type of substrate bound to P450(cam) suggesting regions distant from heme-binding site are affected by physicochemical and steric characteristics/protein-substrate interactions of P450(cam) active site. We monitored substrate dynamics in the active site region of P450(cam) by time-resolved substrate anisotropy measurements. The anisotropy decay of substrates bound to P450(cam) indicate that mobility of substrates is modulated by physicochemical and steric characteristics/protein-substrate interactions of local active site structure, and provides an understanding of factors controlling observed hydroxylated products for substrate bound P450(cam) complexes. The present study shows that P450(cam) local and peripheral structural flexibility and heterogeneity along with substrate mobility play an important role in regulating substrate binding orientation during catalysis and accommodating diverse range of substrates within P450(cam) heme pocket.     

Exploring CYP1A1 as anticancer target: homology modeling and in silico inhibitor design

Microsomal cytochrome P450 family 1 enzymes has great importance in the bioactivation of mutagens. P450 catalyzed reactions involve a wide range of substrates, and this versatality is reflected in a structural diversity, evident in the active sites of available P450 structures. This structure offers a template to study further structure-function relationships of alternative substrates and other cytochrome P450 family 1 members. In this paper, we document a homology model of CYP P450 1A1 from Homo sapiens, developed on the basis of template crystal structure of human microsomal P450 1a2 in complex with inhibitor (PDB Id: 2HI4). Homology modeling is performed at the programs, both in the commercial and public realms. We tried to explore CYP1A1 as a potential target for anticancer chemotherapy. To gain an insight into the binding of ligands with enzyme, protein-ligand complex was developed by including information about the known ligand as spatial restraints and optimizing the mutual interactions between the ligand and the binding site. Active site characterization and the study for involvement of specific aminoacids in binding with ligand, facilitates structure based inhibitor design. This study should prove useful in the design and development of potential novel anticancer agents. Figure The refined structure of homology model of CYP1A1     

Structural Analysis of the Streptomyces avermitilis CYP107W1-Oligomycin A Complex and Role of the Tryptophan 178 Residue

CYP107W1 from Streptomyces avermitilis is a cytochrome P450 enzyme involved in the biosynthesis of macrolide oligomycin A. A previous study reported that CYP107W1 regioselectively hydroxylated C12 of oligomycin C to produce oligomycin A, and the crystal structure of ligand free CYP107W1 was determined. Here, we analyzed the structural properties of the CYP107W1-oligomycin A complex and characterized the functional role of the Trp178 residue in CYP107W1. The crystal structure of the CYP107W1 complex with oligomycin A was determined at a resolution of 2.6 Å. Oligomycin A is bound in the substrate access channel on the upper side of the prosthetic heme mainly by hydrophobic interactions. In particular, the Trp178 residue in the active site intercalates into the large macrolide ring, thereby guiding the substrate into the correct binding orientation for a productive P450 reaction. A Trp178 to Gly mutation resulted in the distortion of binding titration spectra with oligomycin A, whereas binding spectra with azoles were not affected. The Gly178 mutant’s catalytic turnover number for the 12-hydroxylation reaction of oligomycin C was highly reduced. These results indicate that Trp178, located in the open pocket of the active site, may be a critical residue for the productive binding conformation of large macrolide substrates.     

High-resolution crystal structure of cytochrome P450cam

The crystal structure of Pseudomonas putida cytochrome P450cam with its substrate, camphor, bound has been refined to R = 0.19 at a normal resolution of 1.63 A. While the 1.63 A model confirms our initial analysis based on the 2.6 A model, the higher resolution structure has revealed important new details. These include a more precise assignment of sequence to secondary structure, the identification of three cis-proline residues, and a more detailed picture of substrate-protein interactions. In addition, 204 ordered solvent molecules have been found, one of which appears to be a cation. The cation stabilizes an unfavorable polypeptide conformation involved in forming part of the active site pocket, suggesting that the cation may be the metal ion binding site associated with the well-known ability of metal ions to enhance formation of the enzyme-substrate complex. Another unusual polypeptide conformation forms the proposed oxygen-binding pocket. A localized distortion and widening of the distal helix provides a pocket for molecular oxygen. An intricate system of side-chain to backbone hydrogen bonds aids in stabilizing the required local disruption in helical geometry. Sequence homologies strongly suggest a common oxygen-binding pocket in all P450 species. Further sequence comparisons between P450 species indicate common three-dimensional structures with changes focused in a region of the molecule postulated to be associated with the control of substrate specificity.     

Binding and oxidation of alkyl 4-nitrophenyl ethers by rabbit cytochrome P450 1A2: evidence for two binding sites

Although most cytochrome P450 (P450) reactions demonstrate saturation kinetics that fit to the standard Michaelis-Menten equation, there are important exceptions where sigmoidal or nonhyperbolic behavior is observed and have been fit instead to kinetic models involving two binding sites. To assess these models, we demonstrate the consistency of a two binding site model to interpret both steady-state kinetics and binding events. Rates of 4-nitrophenol and formaldehyde production from the O-demethylation of 1-methoxy-4-nitrobenzene by P450 1A2 isolated from rabbit liver produced biphasic plots, when plotted against substrate concentration. Experiments confirmed the absence of the further oxidation of the products. Recombinant rabbit P450 1A2 yielded the same maximal velocity and more marked biphasicity. Overall, these steady-state data fit well to kinetic models involving two binding sites. Steady-state studies of substrates with bulkier O-ethyl or O-isopropoxy groups indicated decreased affinity for the second site. Based on binding studies, the affinity of P450 1A2 for these substrates increased 200-fold with the larger alkyl groups. To analyze the single binding site model, competition studies were conducted with 1,4-phenyldiisocyanide and the alkyl 4-nitrophenyl ethers. Although the observed dissociation constants and the competing titrant demonstrated a linear dependence, the affinity for the competing titrant depended on the presence of the other titrant, which violates the single binding site model. Alternatively, we applied a two binding site model to these data to obtain dissociation constants for the binary and ternary complexes. The agreement between the dissociation constants for the heterogeneous complexes supports the appropriateness of the two binding site model. This novel finding for P450 1A2 may be more common than originally perceived for P450s.     

Combining substrate dynamics, binding statistics, and energy barriers to rationalize regioselective hydroxylation of octane and lauric acid by CYP102A1 and mutants

Hydroxylations of octane and lauric acid by Cytochrome P450-BM3 (CYP102A1) wild-type and three active site mutants--F87A, L188Q/A74G, and F87V/L188Q/A74G--were rationalized using a combination of substrate orientation from docking, substrate binding statistics from molecular dynamics simulations, and barrier energies for hydrogen atom abstraction from quantum mechanical calculations. Wild-type BM3 typically hydroxylates medium- to long-chain fatty acids on subterminal (omega-1, omega-2, omega-3) but not the terminal (omega) positions. The known carboxylic anchoring site Y51/R47 for lauric acid, and hydrophobic interactions and steric exclusion, mainly by F87, for octane as well as lauric acid, play a role in the binding modes of the substrates. Electrostatic interactions between the protein and the substrate strongly modulate the substrate's regiodependent activation barriers. A combination of the binding statistics and the activation barriers of hydrogen-atom abstraction in the substrates is proposed to determine the product formation. Trends observed in experimental product formation for octane and lauric acid by wild-type BM3 and the three active site mutants were qualitatively explained. It is concluded that the combination of substrate binding statistics and hydrogen-atom abstraction barrier energies is a valuable tool to rationalize substrate binding and product formation and constitutes an important step toward prediction of product ratios.     

Inhibition of human cytochrome P450 1A2 by flavones: A molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone >3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.     

Structural characterization of the complex between alpha-naphthoflavone and human cytochrome P450 1B1

The atomic structure of human P450 1B1 was determined by x-ray crystallography to 2.7 Å resolution with α-naphthoflavone (ANF) bound in the active site cavity. Although the amino acid sequences of human P450s 1B1 and 1A2 have diverged significantly, both enzymes exhibit narrow active site cavities, which underlie similarities in their substrate profiles. Helix I residues adopt a relatively flat conformation in both enzymes, and a characteristic distortion of helix F places Phe²³¹ in 1B1 and Phe²²⁶ in 1A2 in similar positions for π-π stacking with ANF. ANF binds in a distinctly different orientation in P450 1B1 from that observed for 1A2. This reflects, in part, divergent conformations of the helix B′-C loop that are stabilized by different hydrogen-bonding interactions in the two enzymes. Additionally, differences between the two enzymes for other amino acids that line the edges of the cavity contribute to distinct orientations of ANF in the two active sites. Thus, the narrow cavity is conserved in both P450 subfamily 1A and P450 subfamily 1B with sequence divergence around the edges of the cavity that modify substrate and inhibitor binding. The conservation of these P450 1B1 active site amino acid residues across vertebrate species suggests that these structural features are conserved.     

Inhibition of human cytochrome P450 1A2 by flavones: A molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone >3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.     

Tight-binding inhibition by alpha-naphthoflavone of human cytochrome P450 1A2

Human cytochrome P450 (P450) enzymes exhibit remarkable diversity in their substrate specificities, participating in oxidation reactions of a wide range of xenobiotic drugs. Previously, we reported that alpha-naphthoflavone (ANF) is bound to the recombinant P450 1A2 tightly and stabilizes an overall enzyme conformation. The present study is designed to determine the type of P450 1A2 inhibition exerted by ANF, using two different substrates of P450 1A2, 7-ethoxycoumarin (EOC) and 7-ethoxyresorufin (EOR). ANF is generally known as a competitive inhibitor of the enzyme. However, in our tight-binding enzyme kinetics study, ANF acts as noncompetitive inhibitor in 7-ethoxycoumarin O-deethylation (ECOD) (K(i)=55.0 nM), but as competitive inhibitor in 7-ethoxyresorufin O-deethylation (EROD) (K(i)=1.4 nM). Based on homology modeling studies, ANF is positioned to bind to a hydrophobic cavity next to the active site where it may cause a direct effect on substrate binding. It is agreed with the predicted binding site of ANF in P450 3A4, in which ANF is rather known as a stimulating modulator. Our results suggest that ANF binds near the active site of P450 1A2 and exhibits differential inhibition mechanisms, possibly depending on the molecular structure of the substrate.     

Cooperativity in oxidation reactions catalyzed by cytochrome P450 1A2: highly cooperative pyrene hydroxylation and multiphasic kinetics of ligand binding

Rabbit liver cytochrome P450 (P450) 1A2 was found to catalyze the 5,6-epoxidation of alpha-naphthoflavone (alphaNF), 1-hydroxylation of pyrene, and the subsequent 6-, 8-, and other hydroxylations of 1-hydroxy (OH) pyrene. Plots of steady-state rates of product formation versus substrate concentration were hyperbolic for alphaNF epoxidation but highly cooperative (Hill n coefficients of 2-4) for pyrene and 1-OH pyrene hydroxylation. When any of the three substrates (alphaNF, pyrene, 1-OH pyrene) were mixed with ferric P450 1A2 using stopped-flow methods, the changes in the heme Soret spectra were relatively slow and multiphasic. Changes in the fluorescence of all of the substrates were much faster, consistent with rapid initial binding to P450 1A2 in a manner that does not change the heme spectrum. For binding of pyrene to ferrous P450 1A2, the course of the spectra revealed sequential changes in opposite directions, consistent with P450 1A2 being involved in a series of transitions to explain the kinetic multiphasicity as opposed to multiple, slowly interconverting populations of enzyme undergoing the same event at different rates. Models of rabbit P450 1A2 based on a published crystal structure of a human P450 1A2-alphaNF complex show active site space for only one alphaNF or for two pyrenes. The spectral changes observed for binding and hydroxylation of pyrene and 1-OH pyrene could be fit to a kinetic model in which hydroxylation occurs only when two substrates are bound. Elements of this mechanism may be relevant to other cases of P450 cooperativity.     

The structure of human microsomal cytochrome P450 3A4 determined by X-ray crystallography to 2.05-A resolution

The structure of P450 3A4 was determined by x-ray crystallography to 2.05-A resolution. P450 3A4 catalyzes the metabolic clearance of a large number of clinically used drugs, and a number of adverse drug-drug interactions reflect the inhibition or induction of the enzyme. P450 3A4 exhibits a relatively large substrate-binding cavity that is consistent with its capacity to oxidize bulky substrates such as cyclosporin, statins, taxanes, and macrolide antibiotics. Family 3A P450s also exhibit unusual kinetic characteristics that suggest simultaneous occupancy by smaller substrates. Although the active site volume is similar to that of P450 2C8 (PDB code: 1PQ2), the shape of the active site cavity differs considerably due to differences in the folding and packing of portions of the protein that form the cavity. Compared with P450 2C8, the active site cavity of 3A4 is much larger near the heme iron. The lower constraints on the motions of small substrates near the site of oxygen activation may diminish the efficiency of substrate oxidation, which may, in turn, be improved by space restrictions imposed by the presence of a second substrate molecule. The structure of P450 3A4 should facilitate a better understanding of the substrate selectivity of the enzyme.     

Structure of microsomal cytochrome P450 2B4 complexed with the antifungal drug bifonazole: insight into P450 conformational plasticity and membrane interaction

To better understand ligand-induced structural transitions in cytochrome P450 2B4, protein-ligand interactions were investigated using a bulky inhibitor. Bifonazole, a broad spectrum antifungal agent, inhibits monooxygenase activity and induces a type II binding spectrum in 2B4dH(H226Y), a modified enzyme previously crystallized in the presence of 4-(4-chlorophenyl)imidazole (CPI). Isothermal titration calorimetry and tryptophan fluorescence quenching indicate no significant burial of protein apolar surface nor altered accessibility of Trp-121 upon bifonazole binding, in contrast to recent results with CPI. A 2.3 A crystal structure of 2B4-bifonazole reveals a novel open conformation with ligand bound in the active site, which is significantly different from either the U-shaped cleft of ligand-free 2B4 or the small active site pocket of 2B4-CPI. The O-shaped active site cleft of 2B4-bifonazole is widely open in the middle but narrow at the top. A bifonazole molecule occupies the bottom of the active site cleft, where helix I is bent approximately 15 degrees to accommodate the bulky ligand. The structure also defines unanticipated interactions between helix C residues and bifonazole, suggesting an important role of helix C in azole recognition by mammalian P450s. Comparison of the ligand-free 2B4 structure, the 2B4-CPI structure, and the 2B4-bifonazole structure identifies structurally plastic regions that undergo correlated conformational changes in response to ligand binding. The most plastic regions are putative membrane-binding motifs involved in substrate access or substrate binding. The results allow us to model the membrane-associated state of P450 and provide insight into how lipophilic substrates access the buried active site.