Random mutagenesis of human cytochrome p450 2A6 and screening with indole oxidation products.

Cytochrome P450 (P450) 2A6 mutants from randomized libraries generated in the substrate recognition sequence (SRS) regions were screened in Escherichia coli on the basis of indole metabolism. SRS 3 and 4 libraries yielded colonies that produced indigo at least as well as wild-type (WT) P450 2A6, and some colonies were consistently more blue upon replating. One mutant, F209T, showed indole 3-hydroxylation WT. The double mutant L240C/N297Q consistently produced very blue colonies. Five mutants yielded mixtures of pigments from indole different than WT, as judged by visible spectra and HPLC of products. When bacteria expressing the mutants were grown in the presence of each of 26 substituted indoles, a variety of patterns of formation of different dyes was seen with several of the mutants. This approach has potential value in understanding P450 2A6 function and generating new dyestuffs and other products.     

Analysis of coumarin 7-hydroxylation activity of cytochrome P450 2A6 using random mutagenesis

Cytochrome P450 (P450) 2A6 is an important human enzyme involved in the metabolism of many xenobiotic chemicals including coumarin, indole, nicotine, and carcinogenic nitrosamines. A combination of random mutagenesis and high-throughput screening was used in the analysis of P450 2A6, utilizing a fluorescent coumarin 7-hydroxylation assay. The steady-state kinetic parameters (k(cat) and Km) for coumarin 7-hydroxylation by wild-type P450 2A6 and 35 selected mutants were measured and indicated that mutants throughout the coding region can have effects on activity. Five mutants showing decreased catalytic efficiency (k(cat)/Km) were further analyzed for substrate selectivity and binding affinities and showed reduced catalytic activities for 7-methoxycoumarin O-demethylation, tert-butyl methyl ether O-demethylation, and indole 3-hydroxylation. All mutants except one (K476E) showed decreased coumarin binding affinities (and also higher Km values), indicating that this is a major basis for the decreased enzymatic activities. A recent x-ray crystal structure of P450 2A6 bound to coumarin (Yano, J. K., Hsu, M. H., Griffin, K. J., Stout, C. D., and Johnson, E. F. (2005) Nat. Struct. Mol. Biol. 12, 822-823) indicates that the recovered A481T and N297S mutations appear to be close to coumarin, suggesting direct perturbation of substrate interaction. The decreased enzymatic activity of the K476E mutant was associated with decreases both in NADPH oxidation and the reduction rate of the ferric P450 2A6-coumarin complex. The attenuation is caused in part to lower binding affinity for NADPH-P450 reductase, but the K476E mutant did not achieve the wild-type coumarin 7-hydroxylation activity even at high reductase concentrations.     

Cytochrome P450 3A4-catalyzed testosterone 6beta-hydroxylation stereochemistry, kinetic deuterium isotope effects, and rate-limiting steps

Testosterone 6beta-hydroxylation is a prototypic reaction of cytochrome P450 (P450) 3A4, the major human P450. Biomimetic reactions produced a variety of testosterone oxidation products with 6beta-hydroxylation being only a minor reaction, indicating that P450 3A4 has considerable control over the course of steroid hydroxylation because 6beta-hydroxylation is not dominant in a thermodynamically controlled oxidation of the substrate. Several isotopically labeled testosterone substrates were prepared and used to probe the catalytic mechanism of P450 3A4: (i) 2,2,4,6,6-(2)H(5); (ii) 6,6-(2)H(2); (iii) 6alpha-(2)H; (iv) 6beta-(2)H; and (v) 6beta-(3)H testosterone. Only the 6beta-hydrogen was removed by P450 3A4 and not the 6alpha, indicating that P450 3A4 abstracts hydrogen and rebounds oxygen only at the beta face. Analysis of the rates of hydroxylation of 6beta-(1)H-, 6beta-(2)H-, and 6beta-(3)H-labeled testosterone and application of the Northrop method yielded an apparent intrinsic kinetic deuterium isotope effect ((D)k) of 15. The deuterium isotope effects on k(cat) and k(cat)/K(m) in non-competitive reactions were only 2-3. Some "switching" to other hydroxylations occurred because of 6beta-(2)H substitution. The high (D)k value is consistent with an initial hydrogen atom abstraction reaction. Attenuation of the high (D)k in the non-competitive experiments implies that C-H bond breaking is not a dominant rate-limiting step. Considerable attenuation of a high (D)k value was also seen with a slower P450 3A4 reaction, the O-dealkylation of 7-benzyloxyquinoline. Thus P450 3A4 is an enzyme with regioselective flexibility but also considerable regioselectivity and stereoselectivity in product formation, not necessarily dominated by the ease of C-H bond breaking.     

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.     

Formation of indigo by recombinant mammalian cytochrome P450

The development of bicistronic systems for coexpression of recombinant human cytochrome P450 enzymes (P450s) with their redox partner, NADPH-cytochrome P450 reductase (NPR), has enabled P450 activity to be reconstituted within bacterial cells. During expression of recombinant P450 2E1 and some other forms, we observed the formation of a blue pigment in bacterial cultures. The pigment was extracted from cultures and shown to comigrate with standard indigo on TLC. UV-visible spectroscopy and mass spectrometric analysis provided further support for identification of the pigment as indigo. Indigo is known to form following the spontaneous oxidation of 3-hydroxyindole. Accordingly, we speculated that indole, formed as a breakdown product of tryptophan in bacteria, was hydroxylated by the P450 system, leading to indigo formation. Bacterial membranes containing recombinant P450 2E1 and human NPR were incubated in vitro with indole and shown to catalyze formation of a blue pigment in a time- and cofactor-dependent manner. These studies suggest potential applications of mammalian P450 enzymes in industrial indigo production or in the development of novel colorimetric assays based on indole hydroxylation.     

Investigation of the role of cytochrome P450 2B4 active site residues in substrate metabolism based on crystal structures of the ligand-bound enzyme

Based on the X-ray crystal structures of 4-(4-chlorophenyl)imidazole (4-CPI)- and bifonazole (BIF)-bound P450 2B4, eight active site mutants at six positions were created in an N-terminal modified construct termed 2B4dH and characterized for enzyme inhibition and catalysis. I363A showed a >4-fold decrease in differential inhibition by BIF and 4-CPI (IC(50,BIF)/IC(50,4-CPI)). F296A, T302A, I363A, V367A, and V477A showed a 2-fold decreased k(cat) for 7-ethoxy-4-trifluoromethylcoumarin O-deethylation, whereas V367A and V477F showed an altered K(m). T302A, V367L, and V477A showed >4-fold decrease in total testosterone hydroxylation, whereas I363A, V367A, and V477F showed altered stereo- and regioselectivity. Interestingly, I363A showed a 150-fold enhanced k(cat)/K(m) with testosterone, and yielded a new metabolite. Furthermore, testosterone docking into three-dimensional models of selected mutants based on the 4-CPI-bound structure suggested a re-positioning of residues 363 and 477 to yield products. In conclusion, our results suggest that the 4-CPI-bound 2B4dH/H226Y crystal structure is an appropriate model for predicting enzyme catalysis.     

A P450 BM-3 mutant hydroxylates alkanes, cycloalkanes, arenes and heteroarenes

P450 monooxygenases from microorganisms, similar to those of eukaryotic mitochondria, display a rather narrow substrate specificity. For native P450 BM-3, no other substrates than fatty acids or an indolyl-fatty acid derivative have been reported (Li, Q.S., Schwaneberg, U., Fischer, P., Schmid, R.D., 2000. Directed evolution of the fatty-acid hydroxylase P450BM-3 into an indole-hydroxylating catalyst. Chem. Eur. J. 6 (9), 1531-1536). Engineering the substrate specificity of Bacillus megaterium cytochrome P-450 BM3: hydroxylation of alkyl trimethylammonium compounds. Biochem. J. 327, 537-544). We thus were quite surprised to observe, in the course of our investigations on the rational evolution of this enzyme towards mutants, capable of hydroxylating shorter-chain fatty acids, that a triple mutant P450 BM-3 (Phe87Val, Leu188-Gln, Ala74Gly, BM-3 mutant) could efficiently hydroxylate indole, leading to the formation of indigo and indirubin (Li, Q.S., Schwaneberg, U., Fischer, P., Schmid, R.D., 2000. Directed evolution of the fatty-acid hydroxylase P450BM-3 into an indole-hydroxylating catalyst. Chem. Eur. J. 6 (9), 1531-1536). Indole is not oxidized by the wild-type enzyme; it lacks the carboxylate group by which the proper fatty acid substrates are supposed to be bound at the active site of the native enzyme, via hydrogen bonds to the charged amino acid residues Arg47 and Tyr51. Our attempts to predict the putative binding mode of indole to P450 BM-3 or the triple mutant by molecular dynamics simulations did not provide any useful clue. Encouraged by the unexpected activity of the triple mutant towards indole, we investigated in a preliminary, but systematic manner several alkanes, alicyclic, aromatic, and heterocyclic compounds, all of which are unaffected by the native enzyme, for their potential as substrates. We here report that this triple mutant indeed is capable to hydroxylate a respectable range of other substrates, all of which bear little or no resemblance to the fatty acid substrates of the native enzyme.     

Hydroxylation of indole by laboratory-evolved 2-hydroxybiphenyl 3-monooxygenase

Directed enzyme evolution of 2-hydroxybiphenyl 3-monooxygenase (HbpA; EC ) from Pseudomonas azelaica HBP1 resulted in an enzyme variant (HbpA(ind)) that hydroxylates indole and indole derivatives such as hydroxyindoles and 5-bromoindole. The wild-type protein does not catalyze these reactions. HbpA(ind) contains amino acid substitutions D222V and V368A. The activity for indole hydroxylation was increased 18-fold in this variant. Concomitantly, the K(d) value for indole decreased from 1.5 mm to 78 microm. Investigation of the major reaction products of HbpA(ind) with indole revealed hydroxylation at the carbons of the pyrrole ring of the substrate. Subsequent enzyme-independent condensation and oxidation of the reaction products led to the formation of indigo and indirubin. The activity of the HbpA(ind) mutant monooxygenase for the natural substrate 2-hydroxybiphenyl was six times lower than that of the wild-type enzyme. In HbpA(ind), there was significantly increased uncoupling of NADH oxidation from 2-hydroxybiphenyl hydroxylation, which could be attributed to the substitution D222V. The position of Asp(222) in HbpA, the chemical properties of this residue, and the effects of its substitution indicate that Asp(222) is involved in substrate activation in HbpA.     

Protein engineering of toluene ortho-monooxygenase of Burkholderia cepacia G4 for regiospecific hydroxylation of indole to form various indigoid compounds

Previous work showed that random mutagenesis produced a mutant of toluene ortho-monooxygenase (TOM) of Burkholderia cepacia G4 containing the V106A substitution in the hydroxylase -subunit (TomA3) that changed the color of the cell suspension from wild-type brown to green in rich medium. Here, DNA shuffling was used to isolate a random TOM mutant that turned blue due to mutation TomA3 A113V. To better understand the TOM reaction mechanism, we studied the specificity of indole hydroxylation using a spectrum of colored TOM mutants expressed in Escherichia coli TG1 and formed as a result of saturation mutagenesis at TomA3 positions A113 and V106. Colonies expressing these altered enzymes ranged in color from blue through green and purple to orange; and the enzyme products were identified using thin-layer chromatography, high performance liquid chromatography, and liquid chromatography–mass spectroscopy. Derived from the single TOM template, enzymes were identified that produced primarily isoindigo (wild-type TOM), indigo (A113V), indirubin (A113I), and isatin (A113H and V106A/A113G). The discovery that wild-type TOM formed isoindigo via C-2 hydroxylation of the indole pyrrole ring makes this the first oxygenase shown to form this compound. Variant TOM A113G was unable to form indigo, indirubin, or isoindigo (did not hydroxylate the indole pyrrole ring), but produced 4-hydroxyindole and unknown yellow compounds from C-4 hydroxylation of the indole benzene ring. Mutations at V106 in addition to A113G restored C-3 indole oxidation, so along with C-2 indole oxidation, isatin, indigo, and indirubin were formed. Other TomA3 V106/A113 mutants with hydrophobic, polar, or charged amino acids in place of the Val and/or Ala residues hydroxylated indole at the C-3 and C-2 positions, forming isatin, indigo, and indirubin in a variety of distributions. Hence, for the first time, a single enzyme was genetically modified to produce a wide range of colors from indole.     

Rate-determining steps in phenacetin oxidations by human cytochrome P450 1A2 and selected mutants

Mutants with altered activities were obtained from random libraries of human cytochrome P450 (P450) 1A2 with the putative substrate recognition sequences (SRS) mutated [Parikh, A., Josephy, P. D., and Guengerich, F. P. (1999) Biochemistry 38, 5283-5289]. Six mutants from SRS 2 (E225I, E225N, F226I, and F226Y) and 4 (D320A and V322A) regions were expressed as oligohistidine-tagged proteins, purified to homogeneity, and used to analyze kinetics of individual steps in the catalytic cycle, to determine which reaction steps have been altered. When the wild-type, E225I, E225N, F226I, F226Y, D320A, and V322A proteins were reconstituted with NADPH-P450 reductase, rates of 7-ethoxyresorufin O-deethylation and phenacetin O-deethylation were in accord with those expected from membrane preparations. Within each assay, the values of k(cat)/K(m) varied by 2-3 orders of magnitude, and in the case of E225I and E225N, these parameters were 7-8-fold higher than for the wild-type enzyme. The coupling efficiency obtained from the rates of product formation and NADPH oxidation was low (<20%) in all enzymes. No correlation was found between activities and several individual steps in the catalytic cycle examined, including substrate binding, reduction kinetics, NADPH oxidation, and H(2)O(2) formation. Quench reactions did not show a burst for either phenacetin O-deethylation or formation of the acetol, a minor product, indicating that rate-determining steps occur prior to product formation. Inter- and intramolecular kinetic deuterium isotope effects for phenacetin O-deethylation were 2-3. In the case of phenacetin acetyl hydroxylation (acetol formation), large isotope effects [(D)k(cat) or (D)(k(cat)/K(m)) > 10] were observed, providing evidence for rate-limiting C-H bond cleavage. We suggest that the very high isotope effect for acetol formation reflects rate-limiting hydrogen atom abstraction; the lower isotope effect for O-deethylation may be a consequence of a 1-electron transfer pathway resulting from the low oxidation potential of the substrate phenacetin. These pre-steady-state, steady-state, and kinetic hydrogen isotope effect studies indicate that the rate-limiting steps are relatively unchanged over an 800-fold range of catalytic activity. We hypothesize that these SRS mutations alter steps leading to the formation of the activated Michaelis complex following the introduction of the first electron.     

Oxidations of p-alkoxyacylanilides catalyzed by human cytochrome P450 1A2: structure-activity relationships and simulation of rate constants of individual steps in catalysis

Human cytochrome P450 (P450) 1A2 is involved in the oxidation of many important drugs and carcinogens. The prototype substrate phenacetin is oxidized to an acetol as well as the O-dealkylation product [Yun, C.-H., Miller, G. P., and Guengerich, F. P. (2000) Biochemistry 39, 11319-11329]. In an effort to improve rates of catalysis of P450 1A2 enzymes, we considered a set of p-alkoxyacylanilide analogues of phenacetin and found that variations in the O-alkyl and N-acyl substituents altered the rates of the two oxidation reactions and the ratio of acetol/phenol products. Moving one methylene group of phenacetin from the O-alkyl group to the N-acyl moiety increased rates of both oxidations approximately 5-fold and improved the coupling efficiency (oxidation products formed/NADPH consumed) from 6% to 38%. Noncompetitive kinetic deuterium isotope effects of 2-3 were measured for all O-dealkylation reactions examined with wild-type P450 1A2 and the E225I mutant, which has 6-fold higher activity. A trend of decreasing kinetic deuterium isotope effect for E225I > wild-type > mutant D320A was observed for O-demethylation of p-methoxyacetanilide, which follows the trend for k(cat). The set of O-dealkylation and acetol formation results for wild-type P450 1A2 and the E225I mutant with several of the protiated and deuterated substrates were fit to a model developed for the basic catalytic cycle and a set of microscopic rate constants in which the only variable was the rate of product formation (substrate oxygenation, including hydrogen abstraction). In this model, k(cat) is considerably less than any of the microscopic rate constants and is affected by several individual rate constants, including the rate of formation of the oxygenating species, the rate of substrate oxidation by the oxygenating species, and the rates of generation of reduced oxygen species (H(2)O(2), H(2)O). This analysis of the effects of the individual rate constants provides a framework for consideration of other P450 reactions and rate-limiting steps.     

Threonine 201 in the diiron enzyme toluene 4-monooxygenase is not required for catalysis

The diiron enzyme toluene 4-monooxygenase from Pseudomonas mendocina KR1 catalyzes the NADH- and O(2)-dependent hydroxylation of toluene. A combination of sequence alignments and spectroscopic studies indicate that T4MO has an active site structure closely related to the crystallographically characterized methane monooxygenase hydroxylase. In the methane monooxygenase hydroxylase, active site residue T213 has been proposed to participate in O(2) activation by analogy to certain proposals made for cytochrome P450. In this work, mutagenesis of the comparable residue in the toluene 4-monooxygenase hydroxylase, T201, has been used to investigate the role of an active site hydroxyl group in catalysis. Five isoforms (T201S, T201A, T201G, T201F, and T201K) that retain catalytic activity based on an in vivo indigo formation assay were identified, and detailed characterizations of the purified T201S, T201A, and T201G variants are reported. These isoforms have k(cat) values of 1.2, 1.0, and 0.6 s(-)(1), respectively, and k(cat)/K(M) values that vary by only approximately 4-fold relative to that of the native isoform. Moreover, these isoforms exhibit 80-90% coupling efficiency, which also compares favorably to the >94% coupling efficiency determined for the native isoform. For the T201S, T201A, and T201G isoforms, the regiospecificity of toluene hydroxylation was nearly identical to that of the natural isoform, with p-cresol representing 90-95% of the total product distribution. In contrast, the T201F isoform caused a substantial shift in the product distribution, and gave o- and p-cresol in a 1:1 ratio. In addition, the amount of benzyl alcohol was increased approximately 10-fold with the T201F isoform. For reaction with p-xylene, previous studies have shown that the native isoform reacted to give 4-methybenzyl alcohol and 2, 5-dimethylphenol in a 4:1 ratio [Pikus, J. D., Studts, J. M., McClay, K., Steffan, R. J., and Fox, B. G. (1997) Biochemistry 36, 9283-9289]. For comparison, the T201S, T201A, and T201F isoforms gave a slightly relaxed 3:1 ratio of these products, while the T201G isoform gave a dramatically relaxed 1:1 ratio. On the basis of these studies, we conclude that the hydroxyl group of T201 is not essential to maintaining the turnover rate or the coupling of the toluene 4-monooxygenase complex. However, changing the volume occupied by the side chain at the position of T201 can lead to alterations in the regiospecificity of the hydroxylation, presumably by producing different orientations for substrate binding during catalysis.     

Expansion of substrate specificity of cytochrome P450 2A6 by random and site-directed mutagenesis

The natural product indole is a substrate for cytochrome P450 2A6. Mutagenesis of P450 2A6 was done to expand its capability in the oxidization of bulky substituted indole compounds, which are not substrates for the wild-type enzyme or the double mutant L240C/N297Q, as determined in our previous work (Wu, Z.-L., Aryal, P., Lozach, O., Meijer, L., and Guengerich, F. P. (2005) Chem. Biodivers. 2, 51-65). Error-prone PCR and site-directed mutagenesis led to the identification of two critical amino acid residue changes (N297Q and I300V) that achieve the purpose. The new mutant (N297Q/I300V) was able to oxidize both 4- and 5-benzyloxy(OBzl)indoles to form colored products. Both changes were required for oxidation of these bulky substrates. The colored product derived from 5-OBzl-indole was mainly 5,5'-di-OBzl-indirubin, whereas the dominant blue dye isolated upon incubations with 4-OBzl-indole was neither an indigo nor an indirubin. Two-dimensional NMR experiments led to assignment of the structure as 4-OBzl-2-(4'-OBzl-1',7'-dihydro-7'-oxo-6'H-indol-6'-ylidene)indolin-3-one, in which a pyrrole ring and a benzene ring are connected with a double bond instead of the pyrrole-pyrrole connection of other indigoids. Monomeric oxidation products were also isolated and characterized; three phenols (4-OBzl-1H-indol-5-ol, 4-OBzl-1H-indol-6-ol, and 4-OBzl-1H-indol-7-ol) and one quinone (4-OBzl-1H-indole-6,7-dione, the postulated immediate precursor of the final blue dye) were identified. The results are interpreted in the context of a crystal structure of a P450 2A6-coumarin complex. The I300V change opens an additional pocket to accommodate the OBzl bulk. The N2297Q change is postulated to generate a hydrogen bond between Gln and the substrate oxygen. Thus, the substrate specificity of P450 2A6 was expanded, and new products were obtained in this study.     

Selection of human cytochrome P450 1A2 mutants with enhanced catalytic activity for heterocyclic amine N-hydroxylation

Cytochrome P450 (P450) 1A2 is the major enzyme involved in the metabolism of 2-amino-3,5-dimethylimidazo[4,5-f]quinoline (MeIQ) and other heterocyclic arylamines and their bioactivation to mutagens. Random mutant libraries of human P450 1A2, in which mutations were made throughout the entire open reading frame, were screened with Escherichia coli DJ3109pNM12, a strain designed to bioactivate MeIQ and detect mutagenicity of the products. Mutant clones with enhanced activity were confirmed using quantitative measurement of MeIQ N-hydroxylation. Three consecutive rounds of random mutagenesis and screening were performed and yielded a highly improved P450 1A2 mutant, SF513 (E225N/Q258H/G437D), with >10-fold increased MeIQ activation based on the E. coli genotoxicity assay and 12-fold enhanced catalytic efficiency (k(cat)/K(m)) in steady-state N-hydroxylation assays done with isolated membrane fractions. SF513 displayed selectively enhanced activity for MeIQ compared to other heterocyclic arylamines. The enhanced catalytic activity was not attributed to changes in any of several individual steps examined, including substrate binding, total NADPH oxidation, or H(2)O(2) formation. Homology modeling based on an X-ray structure of rabbit P450 2C5 suggested that the E225N and Q258H mutations are located in the F-helix and G-helix, respectively, and that the G437D mutation is in the "meander" region, apparently rather distant from the substrate. In summary, the approach generated a mutant enzyme with selectively elevated activity for a single substrate, even to the extent of a difference of a single methyl group, and several mutations had interacting roles in the development of the selected mutant protein.     

Kinetics of cytochrome P450 2E1-catalyzed oxidation of ethanol to acetic acid via acetaldehyde.

The P450 2E1-catalyzed oxidation of ethanol to acetaldehyde is characterized by a kinetic deuterium isotope effect that increases K(m) with no effect on k(cat), and rate-limiting product release has been proposed to account for the lack of an isotope effect on k(cat) (Bell, L. C., and Guengerich, F. P. (1997) J. Biol. Chem. 272, 29643-29651). Acetaldehyde is also a substrate for P450 2E1 oxidation to acetic acid, and k(cat)/K(m) for this reaction is at least 1 order of magnitude greater than that for ethanol oxidation to acetaldehyde. Acetic acid accounts for 90% of the products generated from ethanol in a 10-min reaction, and the contribution of this second oxidation has been overlooked in many previous studies. The noncompetitive intermolecular kinetic hydrogen isotope effects on acetaldehyde oxidation to acetic acid ((H)(k(cat)/K(m))/(D)(k(cat)/K(m)) = 4.5, and (D)k(cat) = 1.5) are comparable with the isotope effects typically observed for ethanol oxidation to acetaldehyde, and k(cat) is similar for both reactions, suggesting a possible common catalytic mechanism. Rapid quench kinetic experiments indicate that acetic acid is formed rapidly from added acetaldehyde (approximately 450 min(-1)) with burst kinetics. Pulse-chase experiments reveal that, at a subsaturating concentration of ethanol, approximately 90% of the acetaldehyde intermediate is directly converted to acetic acid without dissociation from the enzyme active site. Competition experiments suggest that P450 2E1 binds acetic acid and acetaldehyde with relatively high K(d) values, which preclude simple tight binding as an explanation for rate-limiting product release. The existence of a rate-determining step between product formation and release is postulated. Also proposed is a conformational change in P450 2E1 occurring during the course of oxidation and the discrimination of P450 2E1 between acetaldehyde and its hydrated form, the gem-diol. This multistep P450 reaction is characterized by kinetic control of individual reaction steps and by loose binding of all ligands.     

Directed evolution of mammalian cytochrome P450 2B1: mutations outside of the active site enhance the metabolism of several substrates, including the anticancer prodrugs cyclophosphamide and ifosfamide

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.     

Oxidation of methoxyphenethylamines by cytochrome P450 2D6. Analysis of rate-limiting steps

Cytochrome P450 (P450) 2D6 is involved in the oxidation of a large fraction ( approximately 30%) of drugs used by humans and also catalyzes the O-demethylation of the model substrates 3- and 4-methoxyphenethylamine followed by subsequent ring hydroxylation to dopamine. Burst kinetics were not observed; rate-limiting step(s) must occur prior to product formation. Rates of reduction of ferric P450 2D6 were stimulated by 3- or 4-methoxyphenethylamine or the inhibitor quinidine; reduction is not the most rate-limiting step. The non-competitive intramolecular deuterium isotope effect, an estimate of the intrinsic isotope effect, for 4-methoxyphenethylamine O-demethylation was 9.6. Intermolecular non-competitive deuterium isotope effects of 3.1-3.8 were measured for k(cat) and k(cat)/K(m) for both O-demethylation reactions, implicating at least partially rate-limiting C-H bond breaking. Simulation of steady-state kinetic data yielded a catalytic mechanism dominated by the rates of (i) Fe(2+)O(2)(-) protonation (plus O-O bond scission) and (ii) C-H bond breaking, consistent with the appearance of the spectral intermediates in the steady state, attributed to iron-oxygen complexes. However, all the rates of individual steps (or rates of combined steps) are considerably higher than k(cat), and the contributions of several steps must be considered in understanding rates of the P450 2D6 reactions.     

Reconstructing the substrate for uracil DNA glycosylase: tracking the transmission of binding energy in catalysis.

The DNA repair enzyme uracil DNA glycosylase (UDG) is a powerful N-glycohydrolase that cleaves the glycosidic bond of deoxyuridine in DNA. We have investigated the role of substrate binding energy in catalysis by systematically dismantling the optimal substrate Ap(+1)UpA(-1)pA(-2) by replacing the nucleotides at the +1, -1, or -2 position with a tetrahydrofuran abasic site nucleotide (D), a 3-hydroxypropyl phosphodiester spacer (S), a phosphate monoester (p), or a hydroxyl group (h). Contrary to previous reports, the minimal substrate for UDG is 2'-deoxyuridine (hUh). UDG has a significant catalytic efficiency (CE) for hUh of 4 x 10(7) M(-1) [CE = (k(cat)/K(m))(1/k(non)), where k(non) is the rate of the spontaneous hydrolysis reaction of hUh at 25 degrees C]. Addition of +1 and -1 phosphate monoanions to form pUp increases k(cat)/K(m) by 45-fold compared to that of hUh. The k(cat)/K(m) for pUp, but not pU or Up, is found to decrease by 20-fold over the pH range of 6-9 with a pK(a) of 7.1, which is identical to the pK(a) values for deprotonation of the +1 and -1 phosphate groups determined by the pH dependence of the (31)P NMR chemical shifts. This pH dependence indicates that binding of the pUp tetraanion is disfavored, possibly due to unfavorable desolvation or electrostatic properties of the highly charged +1 and -1 phosphate groups. Addition of flexible hydroxypropyl groups to the +1 and -1 positions to make SpUpS increases k(cat)/K(m) by more than 10(5)-fold compared to that of hUh, which is a 20-fold greater effect than observed with rigid D substituents in these positions (i.e., DpUpD). The -2 phosphoester or nucleotide is found to increase the reactivity of trimer substrates with rigid furanose rings or nucleotides in the +1 and -1 positions by 1300-270000-fold (i.e., DpUpD --> DpUpDpA or ApUpA --> ApUpApA). In contrast, the -2 nucleotide provides only an 8-fold rate enhancement when appended to the substrate containing the more flexible +1 and -1 S substituents (SpUpS --> SpUpSpA). These context-dependent effects of a -2 nucleotide are interpreted in terms of a mechanism in which the binding energy of this "handle" is used drive the rigid +1 and -1 A or D substituents into their binding pockets, resulting in a net catalytic benefit of -4.3 to -7.5 kcal/mol. Taken together, these results systematically track how UDG uses distant site binding interactions to produce an overall four billion-fold increase in CE compared to that of the minimal substrate hUh.     

Monooxygenation of an aromatic ring by F43W/H64D/V68I myoglobin mutant and hydrogen peroxide. Myoglobin mutants as a model for P450 hydroxylation chemistry

Myoglobin (Mb) is used as a model system for other heme proteins and the reactions they catalyze. The latest novel function to be proposed for myoglobin is a P450 type hydroxylation activity of aromatic carbons (Watanabe, Y., and Ueno, T. (2003) Bull. Chem. Soc. Jpn. 76, 1309-1322). Because Mb does not contain a specific substrate binding site for aromatic compounds near the heme, an engineered tryptophan in the heme pocket was used to model P450 hydroxylation of aromatic compounds. The monooxygenation product was not previously isolated because of rapid subsequent oxidation steps (Hara, I., Ueno, T., Ozaki, S., Itoh, S., Lee, K., Ueyama, N., and Watanabe, Y. (2001) J. Biol. Chem. 276, 36067-36070). In this work, a Mb variant (F43W/H64D/V68I) is used to characterize the monooxygenated intermediate. A modified (+16 Da) species forms upon the addition of 1 eq of H2O2. This product was digested with chymotrypsin, and the modified peptide fragments were isolated and characterized as 6-hydroxytryptophan using matrix-assisted laser desorption ionization time-of-flight tandem mass spectroscopy and 1H NMR. This engineered Mb variant represents the first enzyme to preferentially hydroxylate the indole side chain of Trp at the C6 position. Finally, heme extraction was used to demonstrate that both the formation of the 6-hydroxytryptophan intermediate (+16 Da) and subsequent oxidation to form the +30 Da final product are catalyzed by the heme cofactor, most probably via the compound I intermediate. These results provide insight into the mechanism of hydroxylation of aromatic carbons by heme proteins, demonstrating that non-thiolate-ligated heme enzymes can perform this function. This establishes Mb compound I as a model for P450 type aromatic hydroxylation chemistry.     

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.