The role of cytochrome 2B1 substrate recognition site residues 115, 294, 297, 298, and 362 in the oxidation of steroids and 7-alkoxycoumarins

Functional Characterization of Allelic Variations of Human Cytochrome P450 2C8 (V181I, I244V, I331T, and L361F)

The human cytochrome P450 2C8 is responsible for the metabolism of various clinical drugs as well as endogenous fatty acids. Allelic variations can significantly influence the metabolic outcomes. In this study, we characterize the functional effects of four nonsynonymous single nucleotide polymorphisms *15, *16, *17, and *18 alleles recently identified in cytochrome P450 2C8. The recombinant allelic variant enzymes V181I, I244V, I331T, and L361F were successfully expressed in Escherichia coli and purified. The steady-state kinetic analysis of paclitaxel 6-hydroxylation revealed a significant reduction in the catalytic activities of the V181I, I244V, and L361F variants. The calculated catalytic efficiency (k_cat/K_m) of these variants was 5-26% of that of the wild-type enzyme. The reduced activities were due to both decreased k_cat values and increased K_m values of the variants. The epoxidation of arachidonic acid by the variants was analyzed. The L361F variant only exhibited 4-6% of the wild-type catalytic efficiency in ω-9- and ω-6-epoxidation reactions to produce 11,12-epoxyeicosatrienoic acid (EET) and 14,15-EET, respectively. These reductions were mainly due to a decrease in the k_cat value of the L361F variant. The binding titration analysis of paclitaxel and arachidonic acid showed that all variants had similar affinities to those of the wild-type (10-14 μM for paclitaxel and 20-49 μM for arachidonic acid). The constructed paclitaxel docking model of the variant enzyme suggests that the L361F substitution leads to the incorrect orientation of paclitaxel in the active site, with the 6'C of paclitaxel displaced from the productive catalytic location. This study suggests that individuals carrying the newly identified P450 2C8 allelic variations are likely to have an altered metabolism of clinical medicines and production of fatty acid-derived signal molecules.

Newly identified tree shrew cytochrome P450 2B6 (CYP2B6) and pig CYP2B6b are functional drug-metabolising enzymes

Tree shrews have high phylogenetic affinity to humans and are used in various fields of biomedical research, especially hepatitis virus infection; however, cytochromes P450 (P450s or CYPs) have not been investigated in this species.In this study, tree shrew CYP2B6 and pig CYP2B6b were newly identified and had amino acid sequences highly identical (80% and 78%, respectively) to human CYP2B6, containing sequence motifs characteristic of P450s.Phylogenetic analysis revealed that novel tree shrew CYP2B6 was more closely related to known human CYP2B6 than dog, pig, or rat CYP2Bs are.Among the tissue types analysed, tree shrew CYP2B6 mRNA was preferentially expressed in liver and lung, whereas pig CYP2B6b mRNA was preferentially expressed in jejunum and lung.Tree shrew CYP2B6 and pig CYP2B6b proteins heterologously expressed in Escherichia coli metabolised human CYP2B6 substrates efavirenz, ethoxycoumarin, propofol, and testosterone, suggesting that these novel CYP2Bs are functional drug-metabolizing enzymes in liver and/or lung.

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 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.

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.

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.

Molecular dynamics analysis reveals structural insights into mechanism of nicotine N-demethylation catalyzed by tobacco cytochrome P450 mono-oxygenase

CYP82E4, a cytochrome P450 monooxygenase, has nicotine N-demethylase (NND) activity, which mediates the bioconversion of nicotine into nornicotine in senescing tobacco leaves. Nornicotine is a precursor of the carcinogen, tobacco-specific nitrosamine. CYP82E3 is an ortholog of CYP82E4 with 95% sequence identity, but it lacks NND activity. A recent site-directed mutagenesis study revealed that a single amino acid substitution, i.e., cysteine to tryptophan at the 330 position in the middle of protein, restores the NND activity of CYP82E3 entirely. However, the same amino acid change caused the loss of the NND activity of CYP82E4. To determine the mechanism of the functional turnover of the two molecules, four 3D structures, i.e., the two molecules and their corresponding cys–trp mutants were modeled. The resulting structures exhibited that the mutation site is far from the active site, which suggests that no direct interaction occurs between the two sites. Simulation studies in different biological scenarios revealed that the mutation introduces a conformation drift with the largest change at the F-G loop. The dynamics trajectories analysis using principal component analysis and covariance analysis suggests that the single amino acid change causes the opening and closing of the transfer channels of the substrates, products, and water by altering the motion of the F-G and B-C loops. The motion of helix I is also correlated with the motion of both the F-G loop and the B-C loop and; the single amino acid mutation resulted in the curvature of helix I. These results suggest that the single amino acid mutation outside the active site region may have indirectly mediated the flexibility of the F-G and B-C loops through helix I, causing a functional turnover of the P450 monooxygenase.

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.

Effect of conformational dynamics on substrate recognition and specificity as probed by the introduction of a de novo disulfide bond into cytochrome P450 2B1

The conformational dynamics of cytochrome P450 2B1 (CYP2B1) were investigated through the introduction of a disulfide bond to link the I- and K-helices by generation of a double Cys variant, Y309C/S360C. The consequences of the disulfide bonding were examined both experimentally and in silico by molecular dynamics simulations. Under high hydrostatic pressures, the partial inactivation volume for the Y309C/S360C variant was determined to be -21 cm3mol(-1), which is more than twice as much as those of the wild type (WT) and single Cys variants (Y309C, S360C). This result indicates that the engineered disulfide bond has substantially reduced the protein plasticity of the Y309C/S360C variant. Under steady-state turnover conditions, the S360C variant catalyzed the N-demethylation of benzphetamine and O-deethylation of 7-ethoxy-trifluoromethylcoumarin as the WT did, whereas the Y309C variant retained only 39% of the N-demethylation activity and 66% of the O-deethylation activity compared with the WT. Interestingly, the Y309C/S360C variant restored the N-demethylation activity to the same level as that of the WT but decreased the O-deethylation activity to only 19% of the WT. Furthermore, the Y309C/S360C variant showed increased substrate specificity for testosterone over androstenedione. Molecular dynamics simulations revealed that the engineered disulfide bond altered substrate access channels. Taken together, these results suggest that protein dynamics play an important role in regulating substrate entry and recognition.

Mechanism-based inactivation of CYP2B1 and its F-helix mutant by two tert-butyl acetylenic compounds: covalent modification of prosthetic heme versus apoprotein

The mechanism-based inactivation of cytochrome CYP2B1 [wild type (WT)] and its Thr205 to Ala mutant (T205A) by tert-butylphenylacetylene (BPA) and tert-butyl 1-methyl-2-propynyl ether (BMP) in the reconstituted system was investigated. The inactivation of WT by BPA exhibited a k(inact)/K(I) value of 1343 min(-1)mM(-1) and a partition ratio of 1. The inactivation of WT by BMP exhibited a k(inact)/K(I) value of 33 min(-1)mM(-1) and a partition ratio of 10. Liquid chromatography/tandem mass spectrometry analysis (LC/MS/MS) of the WT revealed 1) inactivation by BPA resulted in the formation of a protein adduct with a mass increase equivalent to the mass of BPA plus one oxygen atom, and 2) inactivation by BMP resulted in the formation of multiple heme adducts that all exhibited a mass increase equivalent to BMP plus one oxygen atom. LC/MS/MS analysis indicated the formation of glutathione (GSH) conjugates by the reaction of GSH with the ethynyl moiety of BMP or BPA with the oxygen being added to the internal or terminal carbon. For the inactivation of T205A by BPA and BMP, the k(inact)/K(I) values were suppressed by 100- and 4-fold, respectively, and the partition ratios were increased 9- and 3.5-fold, respectively. Only one major heme adduct was detected following the inactivation of the T205A by BMP. These results show that the Thr205 in the F-helix plays an important role in the efficiency of the mechanism-based inactivation of CYP2B1 by BPA and BMP. Homology modeling and substrate docking studies were presented to facilitate the interpretation of the experimental results.

Single mutations change CYP2F3 from a dehydrogenase of 3-methylindole to an oxygenase

Pulmonary cytochrome P450 2F3 (CYP2F3) catalyzes the dehydrogenation of the pneumotoxin 3-methylindole (3MI) to an electrophilic intermediate, 3-methyleneindolenine, which is responsible for the toxicity of the parent compound. Members of the CYP2F subfamily are the only enzymes known to exclusively dehydrogenate 3MI, without detectable formation of oxygenation products. Thus, CYP2F3 is an attractive model to study dehydrogenation mechanisms. The purpose of this study was to identify specific residues that could facilitate 3MI dehydrogenation. Both single and double mutations were constructed to study the molecular mechanisms that direct dehydrogenation. Double mutations in substrate recognition sites (SRS) 1 produced an inactive enzyme, while double mutants in SRS 4 did not alter 3MI metabolism. However, double mutations in SRS 5 and SRS 6 successfully introduced oxygenase activity to CYP2F3. Single mutations in SRS 5, SRS 6 and near SRS 2 also introduced 3MI oxygenase activity. Mutants S474H and D361T oxygenated 3MI but also increased dehydrogenation rates, while G214L, E215Q and S475I catalyzed 3MI oxygenation exclusively. A homology model of CYP2F3 was precisely consistent with specific dehydrogenation of 3MI via initial hydrogen atom abstraction from the methyl group. In addition, intramolecular kinetic deuterium isotope studies demonstrated an isotope effect ( K H/ K D) of 6.8. This relatively high intramolecular deuterium isotope effect confirmed the initial hydrogen abstraction step; a mutant (D361T) that retained the dehydrogenation reaction exhibited the same deuterium isotope effect. The results showed that a single alteration, such as a serine to isoleucine change at residue 475, dramatically switched catalytic preference from dehydrogenation to oxygenation.

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

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

Structure-function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure-function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.

Metabolism of coumarin by human P450s: a molecular modelling study

The oxidative metabolism of coumarin via several human cytochrome P450 (CYP) enzymes from families CYP1, CYP2 and CYP3 is rationalized in terms of molecular modelling studies carried out on the key interactions with various amino acid residues in the relevant active sites. The findings from modelling by homology with the CYP2C5 crystallographic template are in agreement with the known metabolism of coumarin in human P450s from the CYP1, CYP2 and CYP3 families, which has been published recently, and with independently reported information from site-directed mutagenesis studies.

Roles of the threonine 407, aspartic acid 417, and threonine 419 residues in P450 2B1 in metabolism

We have previously observed that the quadruple (S407T-N417D-A419T-K473M) and triple (S407T-N17D-A419T) mutants of the chimeric construct of P450 2B1/2B2 do not undergo mechanism-based inactivation by 17alpha-ethynylestradiol (17EE) and tert-butyl 1-methyl-2-propynyl ether (tBMP). The ability of these mutants to metabolize 17EE, benzphetamine, and testosterone has been investigated. The profile for 17EE metabolism by both mutants was characteristic of both wild-types. The two mutants metabolized testosterone to form androstenedione with no formation of the hydroxy products as was seen with both the wild-types. Benzphetamine metabolism by the mutants showed that both mutants exhibited an increased tendency to catalyze demethylation rather than debenzylation. In the presence of the alternate oxidants cumene hydroperoxide and tert-butyl hydroperoxide, the wild-type 2B1 was not inactivated by 17EE. Metabolism of 17EE by 2B1 supported by these alternate oxidants revealed differences in the metabolites that may be related to the inability of 2B1 to be inactivated under these conditions.

POSSIBLE PATHWAY(S) OF TESTOSTERONE EGRESS FROM THE ACTIVE SITE OF CYTOCHROME P450 2B1: A STEERED MOLECULAR DYNAMICS SIMULATION

To probe the possible substrate exit channel(s) in cytochrome P450 (P450) 2B1 and to clarify the role of residues previously identified by site-directed mutagenesis, a homology model was constructed based on the X-ray crystal structure of a P450 2B4-inhibitor complex. Testosterone was docked into the active site of P450 2B1 and was then pulled out through three putative channels using steered molecular dynamics simulations. The results indicated that of the three channels, the "solvent channel," lined by helices E, F, and I and the beta3 hairpin, required the largest rupture force and backbone motion, which rendered it unlikely as an exit route. The relatively small rupture forces and backbone motions for the other two channels suggested them as possible candidates for testosterone passage. The opening of channel 1, located between helices G and I and the B'-C loop, is characterized by rotation of the aromatic ring of Phe297 together with a bending of the B'-C loop. The opening of channel 2, penetrating through the B'-C loop/B' helix, is achieved by an expansion of this region and a small displacement of the backbone. Interestingly, during the egress of testosterone along channel 1, Phe297 and Phe108 appear to act as two clamps to stabilize testosterone binding and prevent it from leaving the active site. Phe115 acts as a gatekeeper for channel 2. These results are in agreement with previous site-directed mutagenesis experiments.

Functional role of residues in the helix B' region of cytochrome P450 2B1

Comparison of several recently determined X-ray crystal structures of mammalian cytochrome P450 family 2 enzymes suggests considerable movement of helix B' when ligands bind. To investigate the functional role of helix B' in P450 2B1, residues 100-109 were substituted with alanine and phenylalanine. Kinetic properties were examined with the typical 2B substrates 7-benzyloxyresorufin, 7-ethoxy-4-trifluoromethylcoumarin, benzphetamine, and testosterone. Several mutants showed 2- to 3-fold changes in k(cat) values and significant differences in catalytic efficiencies among the substrates examined, consistent with structural information suggesting that the helix B' region can adopt multiple conformations with different contact residues depending on the substrate. Homology modeling of P450 2B1 was performed based on an inhibitor-bound P450 2B4 structure, and the docking analyses were consistent with experimental results. The findings suggest that residues in the helix B' region affect regio- and stereoselective oxidation in P450 family 2 enzymes as well as substrate entry.

Genetically engineered bacterial cells co-expressing human cytochrome P450 with NADPH-cytochrome P450 reductase: prediction of metabolism and toxicity of drugs in humans

Genetically engineered bacterial cells expressing human cytochrome P450 (CYP) have been developed as new tools to predict the metabolism and toxicity of drugs in humans. There are various host cells for the heterologous expression of a form of CYP. Among them, bacterial cells such as Escherichia coli (E. coli) have advantages with regard to ease of use and high yield of protein. CYP protein could be first expressed by the modification of the N-terminal amino acid sequence in E. coli cells in 1991. Since then, many forms of human CYP have been successfully expressed in E. coli cells. Since the E. coli cells do not possess endogeneous electron transport systems to support the full catalytic activity of CYP, E. coli strains co-expressing both human CYP and NADPH-cytochrome P450 reductase (OR) have been established. Each form of CYP expressed in the E. coli cells efficiently catalyzed the oxidation of a representative substrate at an efficient rate, indicating that the OR was sufficiently expressed to support the catalytic activity of CYP. According to the studies performed so far, the modification of the N-terminal amino acid sequence of CYP did not seem to affect the catalytic properties of CYP. The human CYP expressed in the E. coli cells were applicable for studies to determine a metabolic pathway(s) of drugs and to estimate kinetic parameters of drug metabolism by human CYP. Drug-drug interactions caused by inhibition of the metabolism of drugs by human CYP could also be examined by in vitro inhibition studies with CYP expressed in the E. coli cells. Recently, human CYP was co-expressed with the OR in Salmonella typhimurium (S. typhimurium) cells used for mutation assay (Ames test) by applying the technology for the expression of human CYP and the OR in E. coli cells, to evaluate whether chemicals including drugs are metabolically activated by human CYP and show mutagenicity. These strains of bacteria are considered as useful tools to study the metabolism and the toxicity of drugs in humans.

Structure of mammalian cytochrome P450 2B4 complexed with 4-(4-chlorophenyl)imidazole at 1.9-A resolution: insight into the range of P450 conformations and the coordination of redox partner binding

A 1.9-A molecular structure of the microsomal cytochrome P450 2B4 with the specific inhibitor 4-(4-chlorophenyl)imidazole (CPI) in the active site was determined by x-ray crystallography. In contrast to the previous experimentally determined 2B4 structure, this complex adopted a closed conformation similar to that observed for the mammalian 2C enzymes. The differences between the open and closed structures of 2B4 were primarily limited to the lid domain of helices F through G, helices B' and C, the N terminus of helix I, and the beta(4) region. These large-scale conformational changes were generally due to the relocation of conserved structural elements toward each other with remarkably little remodeling at the secondary structure level. For example, the F' and G' helices were maintained with a sharp turn between them but are placed to form the exterior ceiling of the active site in the CPI complex. CPI was closely surrounded by residues from substrate recognition sites 1, 4, 5, and 6 to form a small, isolated hydrophobic cavity. The switch from open to closed conformation dramatically relocated helix C to a more proximal position. As a result, heme binding interactions were altered, and the putative NADPH-cytochrome P450 reductase binding site was reformed. This suggests a structural mechanism whereby ligand-induced conformational changes may coordinate catalytic activity. Comparison of the 2B4/CPI complex with the open 2B4 structure yields insights into the dynamics involved in substrate access, tight inhibitor binding, and coordination of substrate and redox partner binding.

Mutagenesis and molecular dynamics suggest structural and functional roles for residues in the N-terminal portion of the cytochrome P450 2B1 I helix

To investigate their potential roles in ligand access, binding, and subsequent metabolism, residues in the N-terminal portion of the cytochrome P450 2B1 I helix were mutated to alanine and phenylalanine. Of the 18 mutants from E286 to S294 only 7 yielded holoprotein in an Escherichia coli expression system. Substitutions at positions 289, 290, 292, and 294 caused >/= 2-fold changes in kcat and/or Km for two or more of the 2B1 substrates examined, testosterone, 7-ethoxy-4-trifluoromethylcoumarin, 7-benzyloxyresorufin, and benzphetamine. I290 substitutions had the largest effects on steady-state parameters for three substrates and increased benzphetamine affinity. Steered molecular dynamics simulations of testosterone egress along the I helix identified hydrophobic interactions with I290, L293, and S294 and water bridges to E286 and S294. Sensitivity of holoprotein formation to substitution and effects on substrate binding and metabolism suggest structural and functional roles for residues in the N-terminus of the cytochrome P450 2B1 I helix.

Key substrate recognition residues in the active site of a plant cytochrome P450, CYP73A1. Homology guided site-directed mutagenesis

CYP73 enzymes are highly conserved cytochromes P450 in plant species that catalyse the regiospecific 4-hydroxylation of cinnamic acid to form precursors of lignin and many other phenolic compounds. A CYP73A1 homology model based on P450 experimentally solved structures was used to identify active site residues likely to govern substrate binding and regio-specific catalysis. The functional significance of these residues was assessed using site-directed mutagenesis. Active site modelling predicted that N302 and I371 form a hydrogen bond and hydrophobic contacts with the anionic site or aromatic ring of the substrate. Modification of these residues led to a drastic decrease in substrate binding and metabolism without major perturbation of protein structure. Changes to residue K484, which is located too far in the active site model to form a direct contact with cinnamic acid in the oxidized enzyme, did not influence initial substrate binding. However, the K484M substitution led to a 50% loss in catalytic activity. K484 may affect positioning of the substrate in the reduced enzyme during the catalytic cycle, or product release. Catalytic analysis of the mutants with structural analogues of cinnamic acid, in particular indole-2-carboxylic acid that can be hydroxylated with different regioselectivities, supports the involvement of N302, I371 and K484 in substrate docking and orientation.

Threonine-205 in the F helix of p450 2B1 contributes to androgen 16 beta-hydroxylation activity and mechanism-based inactivation

Four mutants of Thr-205 in cytochrome p450 2B1 were constructed and expressed in Escherichia coli. The Ser-, Ala-, and Val-mutants displayed stable reduced CO difference spectra and were able to metabolize 7-ethoxy-4-(trifluoromethyl)coumarin, testosterone, androstenedione, and benzphetamine. The Arg-mutant displayed an unstable reduced CO difference spectrum at 450 nm, was concomitantly converted to a denatured form with a peak at 422 nm, and showed no catalytic activity with any of the four substrates tested. The Ser-mutant displayed activity and metabolite profiles for testosterone and androstenedione similar to those of the wild-type p450 2B1 (WT). Substitution of Thr-205 with Ala or Val markedly suppressed the 16 beta-hydroxylation activity but exhibited little effect on the 16 alpha-hydroxylation activity for testosterone and androstenedione. Because 16 beta-hydroxylation activity of androgens is a specific p450 2B subfamily marker and residue 205 is located in the F helix, which forms the ceiling of the active site, we postulate that the gamma-hydroxyl side chain of Thr may play an important role in directing the 16 beta-face of testosterone and androstenedione toward the active site. Surprisingly, the Val-mutant retained full activity for benzphetamine demethylation. When mechanism-based inactivators for p450 2B1 were used to evaluate the susceptibility to inactivation, the Val-mutant was resistant to inactivation by 17 alpha-ethynylestradiol and less sensitive to inactivation by 2-ethynylnaphthalene compared with the WT enzyme. Our results demonstrate the importance of Thr-205 in determining substrate specificity and product formation as well as in influencing the susceptibility of p450 2B1 to mechanism-based inactivators.

Molecular modelling of CYP2B6 based on homology with the CYP2C5 crystal structure: analysis of enzyme-substrate interactions

The results of homology modelling of CYP2B6 based on the CYP2C5 crystal structure is described in terms of substrates and inhibitors binding within the putative active site. In general these results are in agreement with currently available evidence from substrate metabolism, mode of inhibitor action and site-directed mutagenesis experiments within the CYP2B subfamily of enzymes. Consequently, the model based on the CYP2C5 template represents an advance on those models produced from bacterial P450s, such as CYP101 and CYP102. Quantitative Structure-Activity Relationships (QSARs) for substrates binding to CYP2B6 indicate a key role for hydrogen bonding, and lipophilic character, as determined by the log P parameter (where P is the octanol/water partition coefficient), is also of importance for explaining the variation in experimental binding affinity for CYP2B6 substrates. It is possible to estimate the binding energies for typical CYP2B6 substrates based on their properties and interactions with the enzyme, which show good concordance with experimental data in the form of apparent Km values.

A rational approach to Re-engineer cytochrome P450 2B1 regioselectivity based on the crystal structure of cytochrome P450 2C5

The regioselectivity for progesterone hydroxylation by cytochrome P450 2B1 was re-engineered based on the x-ray crystal structure of cytochrome P450 2C5. 2B1 is a high K(m) progesterone 16alpha-hydroxylase, whereas 2C5 is a low K(m) progesterone 21-hydroxylase. Initially, nine individual 2B1 active-site residues were changed to the corresponding 2C5 residues, and the mutants were purified from an Escherichia coli expression system and assayed for progesterone hydroxylation. At 150 microm progesterone, I114A, F297G, and V363L showed 5-15% of the 21-hydroxylase activity of 2C5, whereas F206V showed high activity for an unknown product and a 13-fold decrease in K(m). Therefore, a quadruple mutant, I114A/F206V/F297G/V363L (Q), was constructed that showed 60% of 2C5 progesterone 21-hydroxylase activity and 57% regioselectivity. Based on their 2C5-like testosterone hydroxylation profiles, S294D and I477F alone and in combination were added to the quadruple mutant. All three mutants showed enhanced regioselectivity (70%) for progesterone 21-hydroxylation, whereas only Q/I477F had a higher k(cat). Finally, the remaining three single mutants, V103I, V367L, and G478V, were added to Q/I477F and Q/S294D/I477F, yielding seven additional multiple mutants. Among these, Q/V103I/S294D/I477F showed the highest k(cat) (3-fold higher than that of 2C5) and 80% regioselectivity for progesterone 21-hydroxylation. Docking of progesterone into a three-dimensional model of this mutant indicated that 21-hydroxylation is favored. In conclusion, a systematic approach to convert P450 regioselectivity across subfamilies suggests that active-site residues are mainly responsible for regioselectivity differences between 2B1 and 2C5 and validates the reliability of 2B1 models based on the crystal structure of 2C5.

Analysis of differential substrate selectivities of CYP2B6 and CYP2E1 by site-directed mutagenesis and molecular modeling

Human CYP2B6 and CYP2E1 were used to investigate the extent to which differential substrate selectivities between cytochrome P450 subfamilies reflect differences in active-site residues as opposed to distinct arrangement of the backbone of the enzymes. Reciprocal CYP2B6 and CYP2E1 mutants at active-site positions 103, 209, 294, 363, 367, and 477 (numbering according to CYP2B6) were characterized using the CYP2B6-selective substrate 7-ethoxy-4-trifluoromethylcoumarin, the CYP2E1-selective substrate p-nitrophenol, and the common substrates 7-ethoxycoumarin, 7-butoxycoumarin, and arachidonic acid. This report is the first to study the active site of CYP2E1 by systematic site-directed mutagenesis. One of the most intriguing findings was that substitution of CYP2E1 Phe-477 with valine from CYP2B6 resulted in significant 7-ethoxy-4-trifluoromethylcoumarin deethylation. Use of three-dimensional models of CYP2B6 and CYP2E1 based on the crystal structure of CYP2C5 suggested that deethylation of 7-ethoxy-4-trifluoromethylcoumarin by CYP2E1 is impeded by van der Waals overlaps with the side chain of Phe-477. Interestingly, none of the CYP2B6 mutants acquired enhanced ability to hydroxylate p-nitrophenol. Substitution of residue 363 in CYP2E1 and CYP2B6 resulted in significant alterations of the metabolite profile for the side chain hydroxylation of 7-butoxycoumarin. Probing of CYP2E1 mutants with arachidonic acid indicated that residues Leu-209 and Phe-477 are critical for substrate orientation in the active site. Overall, the study revealed that differences in the side chains of active-site residues are partially responsible for differential substrate selectivities across cytochrome P450 subfamilies. However, the relative importance of active-site residues appears to be dependent on the structural similarity of the compound to other substrates of the enzyme.

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.

Combined three-dimensional quantitative structure-activity relationship analysis of cytochrome P450 2B6 substrates and protein homology modeling

Understanding the basis of the substrate specificity of cytochrome P450 2B6 (CYP2B6) is important for determining the role of this enzyme in drug metabolism and for predicting new substrates. Pharmacophores were generated for 16 structurally diverse CYP2B6 substrates with Catalyst after overlapping the reaction sites. Two pharmacophores were determined for the CYP2B6 binding site. Both include two hydrophobes and one hydrogen bond acceptor. The three-dimensional structure of CYP2B6 was then modeled based on the crystal structure of CYP2C5. Benzyloxyresorufin and 7-ethoxy-4-trifluoromethylcoumarin, the two lowest K(m) substrates in the training set, were then docked in the active site of CYP2B6. The pharmacophores were combined with the CYP2B6 model by comparing the docking results and the mapping of the two substrates with the pharmacophores. The results indicated that the active site of CYP2B6 complements the pharmacophores. The pharmacophores and the CYP2B6 model were used in conjunction to predict the K(m) values of substrates in a test set of five compounds and yielded satisfactory predictions for benzphetamine, cinnarizine, bupropion, and verapamil but not lidocaine. The CYP2B6 model, the pharmacophores, and the combination of the model with these pharmacophores provide insight into the interactions of CYP2B6 with substrates. The pharmacophores may be used as queries to search a database to predict new substrates for CYP2B6 when the reaction site is known (N- or O-dealkylation). For C-hydroxylation, the CYP2B6 model is helpful in evaluating the possible reaction sites in order for the pharmacophores to predict corresponding K(m) values.