Functional interaction between amino-acid residues 242 and 290 in cytochromes P-450 2B1 and 2B11

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

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

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

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

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

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

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

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

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

Decreased susceptibility of the cytochrome P450 2B6 variant K262R to inhibition by several clinically important drugs

Cytochrome P450 (P450) 2B6 metabolizes a number of clinically relevant drugs and is one of the most highly polymorphic human P450 enzymes, with the Lys(262)-->Arg substitution being especially common in several genetic variants. Therefore, K262R (2B6*4) was created in the CYP2B6dH background (N-terminal-modified and C-terminal His-tagged) and expressed in Escherichia coli. The recombinant CYP2B6dH and K262R were purified and studied to investigate the effect of the Lys(262)-->Arg substitution with six of the most potent drug inhibitors of CYP2B6, namely, clopidogrel, clotrimazole, itraconazole, raloxifene, sertraline, and ticlopidine. K262R showed a >3-fold increase in the K(i) values with clopidogrel, itraconazole, and raloxifene and approximately 6-fold increase in K(i) with sertraline compared with CYP2B6dH. Likewise, K262R showed 2-, 4-, and >20-fold higher K(s) values than CYP2B6dH with clopidogrel, sertraline, and itraconazole, respectively. In contrast, when tested with several known type II inhibitors of CYP2B enzymes, K262R showed a 10-fold lower IC(50) with 4-(phenyl)pyridine and approximately 2-fold lower IC(50) with 4-(4-nitrobenzyl)pyridine or 1-(4-phenyl)benzylimidazole than CYP2B6dH. Subsequent analysis predicted possible in vivo drug-drug interactions between the CYP2B6 substrate efavirenz and drug inhibitors clopidogrel, clotrimazole, itraconazole, sertraline, and ticlopidine. Furthermore, Q172H/K262R (2B6*6), which is the most common genetic variant of CYP2B6 harboring K262R, was created in CYP2B6dH, expressed, purified, and characterized for inhibition. Q172H/K262R showed a >6-fold increase in K(i) with sertraline and clopidogrel compared with CYP2B6dH. The results suggest that individuals, especially homozygotes, with the 2B6*4 or 2B6*6 allele might be less susceptible to drug interactions resulting from P450 inhibition.

Identification and analysis of conserved sequence motifs in cytochrome P450 family 2. Functional and structural role of a motif 187RFDYKD192 in CYP2B enzymes

Using a multiple alignment of 175 cytochrome P450 (CYP) family 2 sequences, 20 conserved sequence motifs (CSMs) were identified with the program PCPMer. Functional importance of the CSM in CYP2B enzymes was assessed from available data on site-directed mutants and genetic variants. These analyses suggested an important role of the CSM 8, which corresponds to(187)RFDYKD(192) in CYP2B4. Further analysis showed that residues 187, 188, 190, and 192 have a very high rank order of conservation compared with 189 and 191. Therefore, eight mutants (R187A, R187K, F188A, D189A, Y190A, K191A, D192A, and a negative control K186A) were made in an N-terminal truncated and modified form of CYP2B4 with an internal mutation, which is termed 2B4dH/H226Y. Function was examined with the substrates 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC), 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), 7-benzyloxy-4-(trifluoromethyl)coumarin (7-BFC), and testosterone and with the inhibitors 4-(4-chlorophenyl)imidazole (4-CPI) and bifonazole (BIF). Compared with the template and K186A, the mutants R187A, R187K, F188A, Y190A, and D192A showed > or =2-fold altered substrate specificity, k(cat), K(m), and/or k(cat)/K(m) for 7-MFC and 7-EFC and 3- to 6-fold decreases in differential inhibition (IC(50,BIF)/IC(50,4-CPI)). Subsequently, these mutants displayed 5-12 degrees C decreases in thermal stability (T(m)) and 2-8 degrees C decreases in catalytic tolerance to temperature (T(50)) compared with the template and K186A. Furthermore, when R187A and D192A were introduced in CYP2B1dH, the P450 expression and thermal stability were decreased. In addition, R187A showed increased activity with 7-EFC and decreased IC(50,BIF)/IC(50,4-CPI) compared with 2B1dH. Analysis of long range residue-residue interactions in the CYP2B4 crystal structures indicated strong hydrogen bonds involving Glu(149)-Asn(177)-Arg(187)-Tyr(190) and Asp(192)-Val(194), which were significantly-reduced/abolished by the Arg(187)-->Ala and Asp(192)-->Alasubstitutions, respectively.

Re-engineering cytochrome P450 2B11dH for enhanced metabolism of several substrates including the anti-cancer prodrugs cyclophosphamide and ifosfamide

Based on recent directed evolution of P450 2B1, six P450 2B11 mutants at three positions were created in an N-terminal modified construct termed P450 2B11dH and characterized for enzyme catalysis using five substrates. Mutant I209A demonstrated a 3.2-fold enhanced k(cat)/K(m) for 7-ethoxy-4-trifluoromethylcourmarin O-deethylation, largely due to a dramatic decrease in K(m) (0.72 microM vs. 18 microM). I209A also demonstrated enhanced selectivity for testosterone 16beta-hydroxylation over 16alpha-hydroxylation. In contrast, V183L showed a 4-fold increased k(cat) for 7-benzyloxyresorufin debenzylation and a 4.7-fold increased k(cat)/K(m) for testosterone 16alpha-hydroxylation. V183L also displayed a 1.7-fold higher k(cat)/K(m) than P450 2B11dH with the anti-cancer prodrugs cyclophosphamide and ifosfamide, resulting from a approximately 4-fold decrease in K(m). Introduction of the V183L mutation into full-length P450 2B11 did not enhance the k(cat)/K(m). Overall, the re-engineered P450 2B11dH enzymes exhibited enhanced catalytic efficiency with several substrates including the anti-cancer prodrugs.

Directed Evolution of Mammalian Cytochrome P450 2B1

Cytochrome P450 2B1 has been subjected to directed evolution to investigate the role of amino acid residues outside of the active site and to engineer novel, more active P450 catalysts. A high throughput screening system was developed to measure H(2)O(2)-supported oxidation of the marker fluorogenic substrate 7-ethoxy-4-trifluoromethylcoumarin (7-EFC). Random mutagenesis by error-prone polymerase chain reaction and activity screening were optimized using the L209A mutant of P450 2B1 in an N-terminally modified construct with a C-terminal His tag (P450 2B1dH). Two rounds of mutagenesis and screening and one subcloning step yielded the P450 2B1 quadruple mutant V183L/F202L/L209A/S334P, which demonstrated a 6-fold higher k(cat) than L209A. Further random or site-directed mutagenesis did not improve the activity. When assayed in an NADPH-supported reconstituted system, V183L/L209A demonstrated lower 7-EFC oxidation than L209A. Therefore, F202L/L209A/S334P was generated, which showed a 2.5-fold higher k(cat)/K(m) for NADPH-dependent 7-EFC oxidation than L209A. F202L/L209A/S334P also showed enhanced catalytic efficiency with 7-benzyloxyresorufin, benzphetamine, and testosterone, and a 10-fold increase in stereoselectivity for testosterone 16alpha-versus 16beta-hydroxylation compared with 2B1dH. Enhanced catalytic efficiency of F202L/L209A/S334P was also retained in the full-length P450 2B1 background with 7-EFC and testosterone as substrates. Finally, the individual mutants were tested for metabolism of the anti-cancer prodrugs cyclophosphamide and ifosfamide. Several of the mutants showed increased metabolism via the therapeutically beneficial 4-hydroxylation pathway, with L209A/S334P showing 2.8-fold enhancement of k(cat)/K(m) with cyclophosphamide and V183L/L209A showing 3.5-fold enhancement with ifosfamide. Directed evolution can thus be used to enhance P450 2B1 catalytic efficiency across a panel of substrates and to identify functionally important residues distant from the active site.

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.

Activation of the anticancer prodrugs cyclophosphamide and ifosfamide: identification of cytochrome P450 2B enzymes and site-specific mutants with improved enzyme kinetics

Cyclophosphamide (CPA) and ifosfamide (IFA) are oxazaphosphorine anticancer prodrugs metabolized by two alternative cytochrome P450 (P450) pathways, drug activation by 4-hydroxylation and drug inactivation by N-dechloroethylation, which generates the neurotoxic and nephrotoxic byproduct chloroacetaldehyde. CPA and IFA metabolism catalyzed by P450s 2B1, 2B4, 2B5, and seven site-specific 2B1 mutants was studied in a reconstituted Escherichia coli expression system to identify residues that contribute to the unique activities and substrate specificities of these enzymes. The catalytic efficiency of CPA 4-hydroxylation by rat P450 2B1 was 10- to 35-fold higher than that of rabbit P450 2B4 or 2B5. With IFA, approximately 50% of metabolism proceeded via N-dechloroethylation for 2B1 and 2B4, whereas CPA N-dechloroethylation corresponded to only approximately 3% of total metabolism (2B1) or was absent (2B4, 2B5). Improved catalytic efficiency of CPA and IFA 4-hydroxylation was obtained upon substitution of 2B1 Ile-114 by Val, and replacement of Val-363 by Leu or Ile selectively suppressed CPA N-dechloroethylation >or=90%. P450 2B1-V367A, containing the Ala replacement found in 2B5, exhibited only approximately 10% of wild-type 2B1 activity for both substrates. Canine P450 2B11, which has Val-114, Leu-363, and Val-367, was therefore predicted to be a regioselective CPA 4-hydroxylase with high catalytic efficiency. Indeed, P450 2B11 was 7- to 8-fold more active as a CPA and IFA 4-hydroxylase than 2B1, exhibited a highly desirable low K(m) (80-160 microM), and catalyzed no CPA N-dechloroethylation. These findings provide insight into the role of specific P450 2B residues in oxazaphosphorine metabolism and pave the way for gene therapeutic applications using P450 enzymes with improved catalytic activity toward these anticancer prodrug substrates.

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.

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.

Structure-function relationships of rat liver CYP3A9 to its human liver orthologs: site-directed active site mutagenesis to a progesterone dihydroxylase

CYP3A9 is an estrogen-inducible ortholog of human liver CYP3A4 with 76.5% sequence identity to CYP3A4. Unlike CYP3A4, it is a very poor testosterone 6beta- and 2beta-hydroxylase, but a relatively better catalyst of progesterone monohydroxylation largely at 6beta, 16alpha, and 21 positions with negligible 6beta, 21-dihydroxylation. We reasoned that such differences in substrate catalyses must be due to differences in the active site architecture of each CYP3A enzyme. Indeed, alignment of CYP3A4 substrate recognition sites (SRSs) with the corresponding regions of CYP3A9 sequence revealed that of the 22 fully divergent residues, 4 reside in SRS regions [P107N (SRS-1), M371G (SRS-5), and L479K and G480Q (SRS-6)]. Accordingly, we substituted these and other divergent CYP3A9 SRS residues with the corresponding residues of CYP3A4 and/or CYP3A5. Our findings of the influence of these site-directed mutations of the CYP3A9 active site on its catalysis of testosterone and three other established but structurally different CYP3A substrates (progesterone, imipramine, and carbamazepine) are described. These findings revealed that some mutations (N107P, N107S, V207T, G371M, and Q480G) not only improved the ability of CYP3A9 to hydroxylate testosterone at the 6beta and 2beta positions, but also converted it into a robust progesterone 6beta, 21-dihydroxylase. The latter in the case of CYP3A9N107P was accompanied by a shift from sigmoidal to hyperbolic enzyme-substrate kinetics. In contrast, the catalytic potential of CYP3A9 mutants K206N, K206S, M240V, and K479L/Q480G was either relatively unchanged or negligible to nonexistent. Together these findings attest to the unique substrate-active site fit of each CYP3A enzyme.

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.

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

At least two substitutions were made at each of five amino acid residues in rat cytochrome P450 2B1 that align to residues of known importance in other P450s. The mutants were histidine tagged for purification from Escherichia coli, and the proteins were assessed for testosterone and 7-alkoxycoumarin oxidation. Alteration of each of the sites studied, Phe-115, Ser-294, Phe-297, Ala-298, and Leu-362, was found to affect overall enzyme activity or the metabolite profile. In particular, most of the mutants, excluding F297A, A298G, and L362F, exhibited significantly altered ratios of 16alpha-hydroxytestosterone:16beta-hydroxytestosterone, with the most dramatic alteration being displayed by A298V. Four 7-butoxycoumarin metabolites were produced by CYP2B1, of which two, 7-hydroxycoumarin and 7-(3-hydroxybutoxy)coumarin, were formed at nearly equal rates. Several mutants, F115A, F297A, F297I, and A298V, exhibited an increased predominance of one of the metabolites. The results from this study illustrate the conservation of functionally important residues across P450 subfamilies and families.

Folding Requirements Are Different between Sterol 14α-Demethylase (CYP51) from Mycobacterium tuberculosis and Human or Fungal Orthologs

Upon sequence alignment of CYP51 sterol 14alpha-demethylase from animals, plants, fungi, and bacteria, arginine corresponding to Arg-448 of CYP51 in Mycobacterium tuberculosis (MT) is conserved near the C terminus of all family members. In MTCYP51 Arg-448 forms a salt bridge with Asp-287, connecting beta-strand 3-2 with helix J. Deletion of the three C-terminal residues of MTCYP51 has little effect on expression of P450 in Escherichia coli. However, truncation of the fourth amino acid (Arg-448) completely abolishes P450 expression. We have investigated whether Arg-448 has other structural or functional roles in addition to folding and whether its conservation reflects conservation of a common folding pathway in the CYP51 family. Characterization of wild type protein and three mutants, R448K, R448I, and R448A, including examination of catalytic activity, secondary and tertiary structure analysis by circular dichroism and tryptophan fluorescence, and studies of both equilibrium and temporal MTCYP51 unfolding behavior, shows that Arg-448 does not play any role in P450 function or maintenance of the native structure. C-terminal truncation of Candida albicans and human CYP51 orthologs reveals that, despite conservation in sequence, the requirement for arginine at the homologous C-terminal position in folding in E. coli is not conserved. Thus, despite similar spatial folds, functionally related but evolutionarily distinct P450s can follow different folding pathways.

2,2',3,3',6,6'-hexachlorobiphenyl hydroxylation by active site mutants of cytochrome P450 2B1 and 2B11

The structural basis of species differences in cytochrome P450 2B-mediated hydroxylation of 2,2',3,3',6,6'-hexachlorobiphenyl (236HCB) was evaluated by using 14 site-directed mutants of cytochrome P450 2B1 and three point mutants of 2B11 expressed in Escherichia coli. To facilitate metabolite identification, seven possible products, including three hydroxylated and four dihydroxylated hexachlorobiphenyls, were synthesized by direct functionalization of precursors and Ullmann and crossed Ullmann reactions. HPLC and GC/MS analysis and comparison with authentic standards revealed that 2B1, 2B11, and all their mutants produced 4, 5-dihydroxy-236HCB and 5-hydroxy-236HCB, while 2B11 L363V and 2B1 I114V mutants also catalyzed hydroxylation at the 4-position. The amount of products formed by 2B1 mutants I114V, F206L, L209A, T302S, V363A, V363L, V367A, I477A, I477L, G478S, I480A, and I480L was smaller than that of the wild type. I477V exhibited unaltered 236HCB metabolism, and I480V produced twice as much dihydroxy product as the wild type. For 2B11, substitution of Val-114 or Asp-290 with Ile decreased the product yields. Replacement of Leu-363 with Val dramatically altered the profile of 236HCB metabolites. In addition to an increase in the overall level of hydroxylation, the mutant mainly catalyzed hydroxylation at the 4-position. Incubation of P450 2B1 with 5-hydroxy-236HCB produced 4,5-dihydroxy-236HCB, which indicates that 4,5-dihydroxy-236HCB may be formed by a direct hydroxylation of 5-hydroxy-236HCB. The findings from this study demonstrate the importance of residues 114, 206, 209, 302, 363, 367, 477, 478, and 480 in 2B1 and 114, 290, and 363 in 2B11 for 236HCB metabolism.

Molecular modelling of CYP2B6, the human CYP2B isoform, by homology with the substrate-bound CYP102 crystal structure: evaluation of CYP2B6 substrate characteristics, the cytochrome b5 binding site and comparisons with CYP2B1 and CYP2B4

1. Molecular modelling studies of CYP2B isoforms from rat (CYP2B1), rabbit (CYP2B4) and man (CYP2B6) are reported, with particular emphasis on substrate interactions with the human CYP2B isoform, CYP2B6. 2. The findings represent an advance on our previous study that focused primarily on the rat CYP2B isoform, CYP2B1, and involved homology modelling with substrate-free CYP102. 3. The current work utilizes the recently published substrate-bound CYP102 crystal structure as a template for construction of the CYP2B subfamily isoforms and shows, in particular, that known CYP2B6 substrate specificity and regioselectivity can be rationalized by putative active site interactions.

Identification of amino acid substitutions that confer a high affinity for sulfaphenazole binding and a high catalytic efficiency for warfarin metabolism to P450 2C19

Human cytochrome P450s 2C9 and 2C19 metabolize many important drugs including tolbutamide, phenytoin, and (S)-warfarin. Although they differ at only 43 of 490 amino acids, sulfaphenazole (SFZ) is a potent and selective inhibitor of P450 2C9 with an IC50 and a spectrally determined binding constant, KS, of <1 microM. P450 2C19 is not affected by SFZ at concentrations up to 100 microM. A panel of CYP2C9/2C19 chimeric proteins was constructed in order to identify the sequence differences that underlie this difference in SFZ binding. Replacement of amino acids 227-338 in 2C19 with the corresponding region of 2C9 resulted in high-affinity SFZ binding (KS approximately 4 microM) that was not seen when a shorter fragment of 2C9 was substituted (227-282). However, replacement of amino acids 283-338 resulted in extremely low holoenzyme expression levels in Escherichia coli, indicating protein instability. A single mutation, E241K, which homology modeling indicated would restore a favorable charge pair interaction between K241 in helix G and E288 in helix I, led to successful expression of this chimera that exhibited a KS < 10 microM for SFZ. Systematic replacement of the remaining differing amino acids revealed that two amino acid substitutions in 2C19 (N286S, I289N) confer high-affinity SFZ binding (KS < 5 microM). When combined with a third substitution, E241K, the resulting 2C19 triple mutant exhibited a high cataltyic efficiency for warfarin metabolism with the relaxed stereo- and regiospecificity of 2C19 and a lower KM for (S)-warfarin metabolism (<10 microM) typical of 2C9.

Identification of residues 286 and 289 as critical for conferring substrate specificity of human CYP2C9 for diclofenac and ibuprofen

Specificity of human CYP2C9 for two substrates, diclofenac and ibuprofen, was studied using chimeras and site-directed mutants of CYP2C9 and the highly related CYP2C19 expressed in Escherichia coli. Data were correlated with the presence of putative substrate recognition sites (SRS). A CYP2C19 chimera containing residues 228-340 (SRS 3 and 4) of 2C9 conferred both diclofenac hydroxylation and 2- and 3-hydroxylation of ibuprofen. The regiospecificity of this construct for metabolism of ibuprofen differed from that of CYP2C9 by favoring 2-hydroxylation over 3-hydroxylation. A CYP2C9 construct containing residues 228-340 of CYP2C19 lacked both diclofenac and ibuprofen hydroxylase activities. When residues 228-282 (containing SRS 3) of CYP2C9 were replaced by those of CYP2C19, the chimera retained appreciable activity for diclofenac and ibuprofen, and tolbutamide activity was inhibited by a specific CYP2C9 inhibitor, sulfaphenazole. This suggested that SRS 3 is not important in conferring specificity. CYP2C9 and CYP2C19 differ in five residues within the region 283-340 (within SRS 4). Mutations to analyze SRS 4 were made on a CYP2C19 chimera containing residues 228-282 of CYP2C9. A single I289N mutation conferred a dramatic increase in diclofenac hydroxylation and a small increase in ibuprofen 2-hydroxylation. A second mutation (N286S and I289N) increased diclofenac hydroxylation and conferred a dramatic increase in ibuprofen 2-hydroxylation. A V288E mutation did not increase activity toward either substrate and decreased activity toward the two substrates in combination with the I289N or the N286S, I289N mutants. Therefore residues 286 and 289 of CYP2C9 are important in conferring specificity for diclofenac and ibuprofen.

Analysis of human cytochrome P450 3A4 cooperativity: construction and characterization of a site-directed mutant that displays hyperbolic steroid hydroxylation kinetics

Cytochrome P450 3A4 is generally considered to be the most important human drug-metabolizing enzyme and is known to catalyze the oxidation of a number of substrates in a cooperative manner. An allosteric mechanism is usually invoked to explain the cooperativity. Based on a structure-activity study from another laboratory using various effector-substrate combinations and on our own studies using site-directed mutagenesis and computer modeling of P450 3A4, the most likely location of effector binding is in the active site along with the substrate. Our study was designed to test this hypothesis by replacing residues Leu-211 and Asp-214 with the larger Phe and Glu, respectively. These residues were predicted to constitute a portion of the effector binding site, and the substitutions were designed to mimic the action of the effector by reducing the size of the active site. The L211F/D214E double mutant displayed an increased rate of testosterone and progesterone 6beta-hydroxylation at low substrate concentrations and a decreased level of heterotropic stimulation elicited by alpha-naphthoflavone. Kinetic analyses of the double mutant revealed the absence of homotropic cooperativity with either steroid substrate. At low substrate concentrations the steroid 6beta-hydroxylase activity of the wild-type enzyme was stimulated by a second steroid, whereas L211F/D214E displayed simple substrate inhibition. To analyze L211F/D214E at a more mechanistic level, spectral binding studies were carried out. Testosterone binding by the wild-type enzyme displayed homotropic cooperativity, whereas substrate binding by L211F/D214E displayed hyperbolic behavior.