Expansion of Substrate Specificity of Cytochrome P450 2A6 by ...

6 downloads 122 Views 698KB Size Report
Jul 27, 2005 - 4-OBzl-indole was neither an indigo nor an indirubin. Two-dimen- sional NMR experiments led to assignment of the structure as. 4-OBzl-2-(4 ...
THE JOURNAL OF BIOLOGICAL CHEMISTRY VOL. 280, NO. 49, pp. 41090 –41100, December 9, 2005 © 2005 by The American Society for Biochemistry and Molecular Biology, Inc. Printed in the U.S.A.

Expansion of Substrate Specificity of Cytochrome P450 2A6 by Random and Site-directed Mutagenesis*□ S

Received for publication, July 27, 2005, and in revised form, September 8, 2005 Published, JBC Papers in Press, October 7, 2005, DOI 10.1074/jbc.M508182200

Zhong-Liu Wu, Larissa M. Podust, and F. Peter Guengerich1 From the Department of Biochemistry and the Center in Molecular Toxicology, Vanderbilt University School of Medicine, Nashville, Tennessee 37232-0146 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.

Microsomal P4502 enzymes (also termed “heme thiolate P450” (1)) are well known for their remarkable capabilities in the catalysis of

* This work was supported by United States Public Health Service Grants R01 CA90426 and P30 ES00267, the Vanderbilt University Medical Center Discovery Grant Program, and a Vanderbilt-Meharry Center for AIDS research award (to L. M. P.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. □ S The on-line version of this article (available at http://www.jbc.org) contains supplemental Figs. 1S–9S, Table 1S, and “Analytical Data for Monomeric Metabolites Derived from 4-OBn Indole and 5-OBn Indole.” 1 To whom correspondence should be addressed: Dept. of Biochemistry and Center in Molecular Toxicology, Vanderbilt University School of Medicine, 638 Robinson Research Bldg., 23rd and Pierce Aves., Nashville, TN 37232-0146. Tel.: 615-322-2261; Fax: 615-322-3141; E-mail: [email protected]. 2 The abbreviations used are: P450, cytochrome P450; OBzl, OBn, and OCH2Ph, benzyloxy; HPLC, high performance liquid chromatography; APCI, atmospheric pressure chemical ionization; NOESY, nuclear Overhauser effect correlation spectroscopy; HMBC, heteronuclear multiple bond correlation; CIGAR, constant time inversedetected gradient accordion rescaled long-range; WT, wild-type; MS, mass spectrometry.

41090 JOURNAL OF BIOLOGICAL CHEMISTRY

diverse oxygenation reactions (2– 6). These enzymes have been studied mostly as the principal catalysts involved in sterol synthesis, drug metabolism, and xenobiotic disposition. The use of P450 enzymes in the area of biocatalysis and fine chemical production is still largely unexploited (7). Random mutagenesis is one of the main approaches in terms of developing P450 enzymes with new functions and can also be used to enhance the knowledge of structure-function relationships of these enzymes and thus provide more information for drug development. Random mutagenesis and molecular breeding approaches can be very useful when there is no crystal structure available to perform rational design, which is still the case for most P450 enzymes today. Furthermore, even rational designs based on crystal structures do have limitations because the enzymes possess both rigidity and flexibility. Predictions are even more complex in the case of oxidoreductase-catalyzed reactions when several factors are involved in the catalytic cycle. In addition to the active site (where substrate is bound) that is obvious in the crystal structure, other residues in the protein may also play important roles. With the availability of various high throughput screening methods and robotic systems, random mutagenesis has become a widely used tool in protein engineering both to improve catalytic efficiency and to investigate structure-function relationships (8). For P450 enzymes, random mutagenesis and molecular breeding are emerging areas. Bacterial P450 102A1 has been used as a starting point to develop catalysts that can hydroxylate alkanes in a regio- and enantioselective manner (9) or that have special features (10). This laboratory has developed P450 1A2 using several strategies to improve the catalytic efficiency toward several substrates (11–13). Molecular breeding of P450 2A6 was done in this laboratory (14), and a double mutant (L240C/N297Q) was selected from libraries constructed by random mutagenesis using randomized primers, each covering four positions of one of the substrate recognition site regions (2), and a staggered extension process method (15). The screening was based on the metabolism of indole, a substrate of P450 2A6 (16, 17) that can be hydroxylated to a product (indoxyl) that dimerizes to the blue compound indigo. The activity of the selected mutant was enhanced toward indole as well as some substituted indoles. This mutant has been used successfully to produce various substituted indigoids with enhanced activities as potent kinase inhibitors (18, 19). In the course of previous work on the biotransformation of 45 indole compounds catalyzed by the L240C/N297Q mutant, we found that this mutant showed limitations in using some substrates, among which were those bearing bulky groups (19). In the x-ray crystal structure of the P450 2A6-coumarin complex recently determined by Johnson and coworkers (20), the active site of P450 2A6 contains large aromatic residues that reduce the volume of the substrate-binding site. The active site is ⬃6-fold smaller than that of P450 2C8 (Protein Data Bank code 1PQ2) (21), consistent with the size of typical substrates such as coumarin (22), nicotine (23, 24), and indole. We are interested in expanding the substrate specificity of P450 2A6

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005

Expansion of Substrate Specificity of P450 2A6 using mutagenesis methods to augment the active pocket so that it can accommodate indoles with bulky substitution groups to increase the diversity of libraries of products generated for use as chemicals for screening of biological activity (19). In this work, error-prone PCR was used to construct randomized libraries, and a colorimetric colony-based method was used in the screening of mutants. Two critical amino acid residue changes that allow oxidization of both 4- and 5-benzyloxy(OBzl) indoles to form colored products were identified. The results are interpreted in the context of the x-ray crystal structure (20).

EXPERIMENTAL PROCEDURES Chemicals—Indole, 4-chloroindole, and 4- and 5-OBzl-indoles were purchased from Aldrich or Lancaster Synthesis (Windham, NH). All other reagents and solvents were obtained from general commercial suppliers and used without further purification. Spectroscopy—UV-visible spectra were recorded either with a Hewlett-Packard 1040A diode array detector (on-line with HPLC) or in Me2SO or H2O using a Cary 14/OLIS or Aminco DW2a/OLIS spectrophotometer for quantitative studies. Mass spectra were obtained in the Vanderbilt facility either on-line with HPLC or via direct infusion using a Finnigan TSQ 7000 triple quadrupole mass spectrometer (Thermo Electron Corp., San Jose, CA) equipped with a standard API-1 atmospheric pressure chemical ionization (APCI) source in the positive or negative ion mode. N2 was used as a sheath gas (50 p.s.i.). The vaporizer temperature was set to 500 °C, and the corona current was maintained at 5 ␮A. The capillary was set at 220 °C and 25 V (or ⫺25 V in the negative mode). The tube lens voltage was set to 80 V (or to ⫺96 V in the negative mode). The collision-induced dissociation was set to ⫺30 V in tandem mass spectrometry experiments. Data acquisition and spectral analysis were conducted with Finnigan ICIS software on a Digital Equipment Alpha workstation. All NMR experiments were carried out in the Vanderbilt facility at 298 K with Bruker AvanceTM DRX 500- and 400-MHz instruments. Samples were dissolved in d6-Me2SO or CD3CN, and semiquantitative estimates of the amounts of indigoids were made by adding known amounts of CH3CO2H and referencing to this internal –CH3 signal (␦ 1.91). Detailed structural information for four was obtained from four separate two-dimensional NMR experiments. Through-space 1H-1H proximity was observed using nuclear Overhauser effect correlation spectroscopy (NOESY) with a mixing time of 500 ms. Single C–H bond connectivity was obtained using heteronuclear multiple quantum correlation methods. Long-range C–H bond connectivity (two to three bonds) was obtained using heteronuclear multiple bond correlation (HMBC) with a 55-ms delay and 1H-13C CIGAR-HMBC with a delay for evolution of long-range couplings being modulated between 6 and 12 Hz and a low pass J-filter to suppress one-bond correlations. The NOESY and COSY spectra of compounds 6, 8, and 10 were recorded (mixing time of 600 ms for NOESY) to determine the hydroxylation positions. All NMR data were processed using XWIN-NMRTM software on an Octane workstation (Silicon Graphics, Mountain View, CA). Expression Systems—The system involved expression of human P450 2A6 in Escherichia coli (25). The P450 2A6 double mutant L240C/N297Q used to construct the library had been selected from libraries generated in selected regions of the P450 2A6 sequence on the basis of their ability to produce colored indigoids (14). All of the wild-type (WT) and mutant P450 2A6 enzymes and human NADPH-cytochrome P450 reductase were expressed using a bicistronic pCW⬘ vector (pCW2A6bc) (26) and a tryptophanase-negative strain (trnA⫺) of E. coli as the host bacteria (19). The expression medium contained 48 g of Terrific Broth/liter (BD Biosciences), 2.0 g of Bacto-peptone/liter (BD Biosciences), 0.4% (v/v)

DECEMBER 9, 2005 • VOLUME 280 • NUMBER 49

glycerol, 1.0 mM 5-aminolevulinic acid, 100 ␮g of ampicillin/liter, 1.0 mM thiamin, 1.0 mM isopropyl ␤-D-thiogalactopyranoside, a 0.025% (v/v) mixture of trace elements (27), and 1% (v/v) starter cultures that had been grown with each newly transformed E. coli cell in LB medium (BD Biosciences) containing 100 ␮g of ampicillin/liter for 16 h at 32 °C with gyratory shaking at 200 rpm. The solid medium used in library screening contained 40 g of LB agar/liter (BD Biosciences), 1.0 mM 5-aminolevulinic acid, 100 ␮g of ampicillin/liter, 1.0 mM thiamin, 1.0 mM isopropyl ␤-D-thiogalactopyranoside, a 0.025% (v/v) mixture of trace elements (27), and 0.5 mM 5-OBzl-indole and was cast in large bioassay Q tray plates (25 ⫻ 25 cm; Genetix, Hampshire, UK). Construction of P450 2A6 Libraries—The open reading frame region of the P450 2A6 gene (1.5 kb) was mutated using a low fidelity PCR method and amplified in a 50-␮l PCR mixture containing 250 ng of pCW2A6bc-L240C/N297Q or pCW2A6bc-WT plasmid, 120 ng of forward primer 5⬘-TAGGAGGTCATATGGCTGCT-3⬘ and reverse primer 5⬘-ATTTCTAGACCGGAAGGCTT-3⬘, 2.5 units of Mutazyme威 DNA polymerase (Stratagene, La Jolla, CA), 0.2 mM each dNTP, and 10⫻ Mutazyme威 reaction buffer (Stratagene). The 1.5-kb amplified PCR library fragment was purified by agarose gel electrophoresis and cloned into the pCW2A6bc vector using the NdeI and XbaI restriction sites. The ligation mixture was transformed into E. coli DH10B ultracompetent cells (Invitrogen) by electroporation. The library plasmid DNA was purified using a QIAprep威 miniprep kit (Qiagen Inc.). Phenotypic Selection Based on Color Formation—The bicistronic library DNA was transformed into the E. coli trnA⫺ strain by heat shock at 41 °C, plated on large bioassay Q tray plates, and allowed to grow at 30 °C for 2 days. Individual colonies with blue color were picked and grown in 3 ml of LB medium containing ampicillin (100 ␮g/ml) as starter cultures for 15 h at 32 °C (see Fig. 1). Plasmid DNA was prepared from the cultures using a QIAprep威 miniprep kit. Verification of the new mutants was performed by inoculating the starter cultures into 1 ml of Terrific Broth expression medium (1%, v/v) fortified with 1 mM 5-OBzl-indole on 24-well plates. The plates were then incubated at 30 °C for 30 h with gyratory shaking at 250 rpm. Clones able to transform 5-OBzl-indole to form blue cultures were selected. Nucleotide Sequence Analysis—Sequencing of the plasmid DNA in the clones that produced blue cultures (following verification) was performed in the Vanderbilt facility using an Applied Biosystems Model 3700 fluorescence sequencing unit with a Taq dye terminator kit (Applied Biosystems, Foster City, CA). The mutated sequences were verified by comparing the sequences of the sense and antisense strands (entire open reading frame). Construction of Site-directed Mutants—Site-directed mutagenesis was conducted using the QuikChange威 mutagenesis kit (Stratagene) following the supplier’s protocol. The primers were designed so that the mutation position was close to the middle with ⬃10 –15 bases of template complementary sequence on both sides and with a melting temperature of ⱖ75 °C (supplemental data, Table 1S). For practical purpose, a pBluescript SK(⫹) vector of the selected P450 2A6 mutant DNA (Bn1; ⬃4.5 kb) was first constructed and subjected to site-directed mutagenesis (see Fig. 1). The 1.5-kb mutated fragment was subcloned into the pCW2A6bc vector (⬃8.5 kb) using the NdeI and XbaI restriction sites. The ligation mixture was transformed into E. coli, and positive clones were confirmed following purification of plasmids. The constructed site-directed mutant plasmids were verified by sequencing both strands of DNA (see above). Whole Cell Assay for Oxidation of Indoles—Cells were incubated at 28 °C for 14 h with gyratory shaking at 200 rpm. Cell pellets containing

JOURNAL OF BIOLOGICAL CHEMISTRY

41091

Expansion of Substrate Specificity of P450 2A6 expressed mutant P450 2A6 and NADPH-cytochrome P450 reductase were harvested from 500 ml of Terrific Broth expression medium of the E. coli trnA⫺ strain by centrifugation at 5000 ⫻ g for 20 min, washed with 200 ml of M9 minimal medium (28), resuspended in 100 ml of M9 minimal medium, and used directly (after measuring the P450 concentration). The reaction mixtures consisted of 200 pmol of P450 (whole cells), 100 mM potassium phosphate buffer (pH 7.4), and varying concentrations of substrate in a total volume of 1 ml. The reactions were carried out at 37 °C in 24-well plates with gyratory shaking at 280 rpm and quenched by the addition of 0.1 ml of 0.5 M aqueous NaOH after 30 min (for indole and 4-chloroindole) or after 100 min (for 4- and 5-OBzlindoles). The formation of colored products (indigo and indirubin) was measured by recording absorbance spectra in the visible region, and the estimated concentrations were calculated based on the extinction coefficient of the major product (subsequently determined). Estimates of kcat and Km were made using nonlinear regression with GraphPad Prism software. Preparation of Membranes and Purification of Recombinant WT P450 2A6 and Mutants—Bacterial inner membranes containing mutant P450 2A6 and NADPH-cytochrome P450 reductase were isolated from 500 ml of Terrific Broth expression medium of the E. coli trnA⫺ strain. The pellets were resuspended in a final volume of 10 ml of 50 mM Tris acetate (pH 7.4) containing 250 mM sucrose and 0.25 mM EDTA and homogenized with a Dounce homogenizer. The homogenized membrane fractions were stored on ice until analysis. P450 enzymes were purified from cholate extracts of E. coli membranes as described previously for P450 1A2 mutants (12). Briefly, membranes were solubilized at 4 °C in 20 mM potassium phosphate buffer (pH 7.4) containing 20% (v/v) glycerol, 1.2 mM magnesium acetate, 10 mM ␤-mercaptoethanol, 0.36% (w/v) Tergitol NP-10, and 0.36% (w/v) sodium cholate. The solubilized proteins (with a C-terminal His5 tag) were purified using a nickel-nitrilotriacetic acid column (Qiagen Inc.) eluted with 400 mM imidazole. The recovered proteins were dialyzed to remove imidazole and stored at ⫺20 °C in 50 mM potassium phosphate buffer (pH 7.7) containing 1 mM EDTA and 20% (v/v) glycerol. The proteins were ⬎95% homogeneous as judged by SDS-PAGE; the apparent specific contents of WT P450 2A6 and the L240C/N297Q, Bn1, Bn1/1, and Bn1/4 mutants were 10.5, 10.3, 8.9, 11.2, and 9.6 nmol of P450/mg of protein, respectively, based on protein estimations done with a commercial BCA method (Pierce) used according to the manufacturer’s instructions. Estimation of Substrate Binding Affinity—The binding affinity of P450 2A6 (WT and mutants) for indole or 5-OBzl-indole was estimated by titration (29). Substrates were dissolved in Me2SO and added to 1.0 ml of 100 mM potassium phosphate buffer (pH 7.4) containing 2– 4 nmol of purified enzyme, and spectra were recorded following each addition (350 –500 nm). The absorbance changes at 390 and 420 nm were used to estimate a spectral Kd (Ks) using GraphPad Prism software, with fitting to a one-site binding hyperbola equation using nonlinear regression. Coumarin 7-Hydroxylation Assay—The reaction mixtures contained 20 pmol of P450 (cell membranes), 100 mM potassium phosphate buffer (pH 7.4), and varying concentrations of coumarin (from a 100-fold stock in H2O) in a total volume of 0.5 ml. The reactions were carried out at 37 °C and quenched with 0.1 ml of 2 M aqueous HCl after 5 min. The product was extracted as described (30) and assayed fluorometrically in a plate reader. Large Scale Biotransformation and Product Separation—The biotransformations were performed in 2.8-liter Fernbach flasks with 500 ml of Terrific Broth expression medium. The cultures were incubated for

41092 JOURNAL OF BIOLOGICAL CHEMISTRY

FIGURE 1. Scheme for preparation and analysis of P450 2A6 mutants. The course of experiments is briefly described. WT P450 2A6 was subjected to random mutagenesis and a staggered extension process (StEP) (3), with selection using enhancement of the ability to convert indole to blue indigo. The result was the double mutant L240C/N297Q (DM) (14). This mutant served as the starting material in this study. The cDNA was subjected to error-prone PCR, and the resulting proteins were screened for their ability to convert 5-OBzl(OBn)-indole to colored products. The result was a mutant (Bn1) that was found to have a total of five mutations relative to WT P450 2A6 (I140M, L240C, N297Q, I300V, and I366V). Site-directed mutagenesis was used to produce all combinations of these substitutions and to establish which are necessary to allow P450 2A6 to oxidize the bulky 4- and 5-OBzl-indoles. The result was Bn1/125, with only two substitutions (N297Q and I300V). SRS, substrate recognition site.

5 h at 30 °C with gyratory shaking at 250 rpm. At that time, the substrate (4- or 5-OBzl-indole) in Me2SO solution (1 M stocks) was added to the medium (0.5 ml), and shaking was continued for 22 h at 30 °C. The pelleted cells and supernatants were separated by centrifugation at 5000 ⫻ g for 30 min. The supernatant was extracted with ethyl acetate (3 ⫻ 200 ml) and dried with Na2SO4. The pellet (from 500 ml of bacterial culture) was stirred with 500 ml of (CH3)2CO and filtered through Celite威 (18, 19). The ethyl acetate and (CH3)2CO extracts were combined and concentrated in vacuo, and the resulting solids were fractionated by semipreparative HPLC. HPLC was done on a Beckman Ultra[sphere column (10 ⫻ 250 mm, 5 ␮m) with a CH3CN/H2O mobile phase at a flow rate of 3 ml/min with 40% (v/v) CH3CN (in H2O) for 10 min, a 40 – 60% CH3CN linear gradient over 20 min, a 60 – 80% CH3CN linear gradient over 10 min, and a 80 –100% CH3CN linear gradient over 5 min. To separate compounds 5 and 6, HPLC conditions of 30% (v/v) CH3CN (in H2O) for 10 min and a 30 –35% CH3CN linear gradient over 20 min were used with the same column. Detection was at 254 or 600 nm (single wavelength) or with the on-line HPLC/diode array detector system. Peaks were collected manually, and CH3CN was removed with a rotary evaporator. The remaining aqueous phase was lyophilized to dryness. In analytical systems, a Beckman Ultrasphere HPLC column (4.6 ⫻ 250 mm, 5 ␮m) was used with a CH3CN/H2O or CH3OH/H2O mobile phase at a flow rate of 1 ml/min. Substrate Docking and Modeling of the P450 2A6 N297Q/I300V Mutant—The coordinates of the P450 2A6 N297Q/I300V mutant were built by introducing substitutions into the WT P450 2A6 coordinates

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005

DECEMBER 9, 2005 • VOLUME 280 • NUMBER 49

b

a

0.062 ⫾ 0.016 0.033 ⫾ 0.009 0.037 ⫾ 0.007 0.074 ⫾ 0.059 0.092 ⫾ 0.027 0.053 ⫾ 0.006 0.19 ⫾ 0.03 0.052 ⫾ 0.011 0.090 ⫾ 0.019 0.12 ⫾ 0.02 0.092 ⫾ 0.032 0.013 ⫾ 0.002 0.046 ⫾ 0.010

149 ⫾ 35 159 ⫾ 38 75 ⫾ 14 1769 ⫾ 1055 96 ⫾ 26 164 ⫾ 17 62 ⫾ 10 129 ⫾ 26 71 ⫾ 14 56 ⫾ 10 196 ⫾ 59 145 ⫾ 25 280 ⫾ 59

min⫺1 ␮M⫺1

0.011 ⫾ 0.001 0.12 ⫾ 0.04 0.079 ⫾ 0.013

␮M

min⫺1

kcat/Km

157 ⫾ 11 222 ⫾ 66 152 ⫾ 23

Indole Km

kcat

None 1.7 ⫾ 0.1 L240C/N297Q 27 ⫾ 4 V140M/L240C/ 12 ⫾ 1 N297Q/I300V/ I366V L240C/N297Q/ 9.3 ⫾ 1.0 I300V/I366V V140M/N297Q/ 5.2 ⫾ 0.6 I300V/I366V V140M/L240C/ 2.8 ⫾ 0.2 I300V/I366V V140M/L240C/ 131 ⫾ 69 N297Q/I366V V140M/L240C/ 8.8 ⫾ 0.9 N297Q/I300V N297Q/I300V/ 8.7 ⫾ 0.4 I366V L240C/N297Q/ 11.9 ⫾ 0.6 I300V L240C/N297Q/ 6.7 ⫾ 0.6 I366V V140M/L240C/ 6.4 ⫾ 0.4 N297Q V140M/N297Q/ 6.5 ⫾ 0.3 I300V N297Q/I300V 18 ⫾ 3 I300V 1.9 ⫾ 0.1 N297Q 13 ⫾ 2

Mutations

14 ⫾ 3 1.9 ⫾ 0.2 8.3 ⫾ 0.9

8.9 ⫾ 1.5

8.0 ⫾ 0.9

14 ⫾ 3

43 ⫾ 21

54 ⫾ 58

3.4 ⫾ 0.2

51 ⫾ 19

6.3 ⫾ 0.5

3.6 ⫾ 0.3

7.0 ⫾ 0.5

1.2 ⫾ 0.1 15 ⫾ 2 14 ⫾ 1

min⫺1

kcat

652 ⫾ 233 174 ⫾ 37 260 ⫾ 47

383 ⫾ 97

126 ⫾ 33

304 ⫾ 98

1620 ⫾ 893

3039 ⫾ 3544

74 ⫾ 15

1212 ⫾ 552

276 ⫾ 34

87 ⫾ 19

83 ⫾ 16

55 ⫾ 5 215 ⫾ 51 186 ⫾ 32

␮M

4-Cl-indole Km

0.021 ⫾ 0.011 0.011 ⫾ 0.003 0.032 ⫾ 0.007

0.023 ⫾ 0.007

0.063 ⫾ 0.018

0.046 ⫾ 0.018

0.027 ⫾ 0.020

0.018 ⫾ 0.028

0.046 ⫾ 0.010

0.042 ⫾ 0.025

0.023 ⫾ 0.003

0.041 ⫾ 0.009

0.084 ⫾ 0.017

0.022 ⫾ 0.002 0.070 ⫾ 0.019 0.075 ⫾ 0.014

min⫺1 ␮M⫺1

kcat/Km

1.7 ⫾ 0.2 — —

1.6 ⫾ 0.2





3.0 ⫾ 0.5

1.7 ⫾ 0.2

0.60 ⫾ 0.04





0.50 ⫾ 0.02

0.38 ⫾ 0.02

— — 0.28 ⫾ 0.02

a

min⫺1

kcat

201 ⫾ 45 — —

288 ⫾ 50





418 ⫾ 98

285 ⫾ 63

54 ⫾ 1





55 ⫾ 7

43 ⫾ 10

— — 51 ⫾ 13

␮M

0.0085 ⫾ 0.0021 — —

0.0056 ⫾ 0.0012





0.0072 ⫾ 0.0021

0.0060 ⫾ 0.0015

0.011 ⫾ 0.001





0.0091 ⫾ 0.0012

0.0088 ⫾ 0.0021

— — 0.0055 ⫾ 0.0015

min⫺1 ␮M⫺1

5-OBzl-indole Km kcat/Km

—, no activity was detected. The numbers 1–5 indicate the positions that were recovered from the Bn1 mutant corresponding to V140M, L240C, N297Q, I300V, and I366V in sequence.

Bn1/125 Bn1/1253 Bn1/1254

Bn1/52

Bn1/45

Bn1/41

Bn1/15

Bn1/12

Bn1/5

Bn1/4

Bn1/3

Bn1/2

Bn1/1b

WT DM Bn1

P450 2A6

Steady-state kinetic parameters for WT P450 2A6 and mutants in the biotransformation of indoles DM, double mutant.

TABLE ONE

3.8 ⫾ 0.6 — —

5.3 ⫾ 2.3





4.8 ⫾ 0.6

1.2 ⫾ 0.1

1.5 ⫾ 0.2





0.91 ⫾ 0.09

0.48 ⫾ 0.03

— — 0.45 ⫾ 0.01

min⫺1

kcat

449 ⫾ 98 — —

864 ⫾ 472





486 ⫾ 86

122 ⫾ 27

139 ⫾ 36





108 ⫾ 25

42 ⫾ 12

— — 47 ⫾ 6

␮M

0.0085 ⫾ 0.0023 — —

0.0061 ⫾ 0.0043





0.0099 ⫾ 0.0021

0.0098 ⫾ 0.0023

0.011 ⫾ 0.003





0.0084 ⫾ 0.0021

0.011 ⫾ 0.003

— — 0.0096 ⫾ 0.0012

min⫺1 ␮M⫺1

4-OBzl-indole Km kcat/Km

Expansion of Substrate Specificity of P450 2A6

JOURNAL OF BIOLOGICAL CHEMISTRY

41093

Expansion of Substrate Specificity of P450 2A6

FIGURE 2. Steady-state kinetics of oxidation of indole compounds by P450 2A6 mutants. A, indole as substrate; B, 4-chloroindole as substrate; C, 5-OBzl(OBn)-indole as substrate; D, 4-OBzl-indole as substrate. Each point presented is the mean ⫾ S.D. (range) of duplicate assays. DM, double mutant L240C/N297Q.

(20) in silico using interactive graphics techniques in the program O (31), followed by manual docking of each of the two substrates into the binding site cavity. A rotamer of Gln capable of donating a hydrogen bond to the substrate bound in the active site was selected. Subsequent energy minimization (conjugate gradient method) was performed using a CNS software routine with no experimental energy term (32). Energy minimization resulted in molecular models of P450 2A6 with bound substrates having a root mean square deviation of 0.55 Å for C-␣ atoms compared with the original structure. Subsequently, ligand coordinates were removed, and protein coordinates were used as a rigid receptor for docking of these same substrates using the FlexX Suite module of Tripos SYBYL software (33), including the pharmacophore constraint of Gln297 in the enzyme active site as a hydrogen bond donor.

RESULTS Library Screening—Two randomized libraries were constructed using error-prone PCR following the manufacturer’s instructions (Stratagene) to adjust the mutation rate to zero to three mutations/1000 bp using either WT P450 2A6 or the double mutant L240C/N297Q as a template (Fig. 1). The screening method was a convenient colorimetric colony-based method that allows rapid screening of thousands of colonies by simply plating them on LB expression plates fortified with the target substrate (5-OBzl-indole in this case) and picking up the blue colonies directly from the plate. LB expression medium was used instead of Terrific Broth to decrease the background color. A tryptophanase-negative strain (trnA⫺) of E. coli (19) was used to avoid the contribution of indole to the products, i.e. to force the P450 2A6 mutants to use only OBzl-substituted indoles. Because of the temperature sensitivity of this strain (␭ prophage can be switched on at 42 °C) (34), a slightly lower temperature than 42 °C was used for the heat shock to achieve a sufficient number of colonies to work with. A total number of ⬃3 ⫻ 103 colonies were screened from each of the two libraries. No blue colonies were found upon screening the WT P450 2A6-derived library. The double mutant-derived library yielded three colonies, designated Bn1 (with a clear blue color), Bn2, and Bn3 (with a brown color), among which only Bn1 yielded a distinctive blue color when incubated in 1 ml of Terrific Broth expression medium fortified with 5-OBzlindole (Fig. 1). This mutant was determined to be P450 2A6 I140M/

41094 JOURNAL OF BIOLOGICAL CHEMISTRY

L240C/N297Q/I300V/I366V by sequence analysis, with three more mutation positions added to the template DNA. Critical Residue Identification and Kinetic Parameters for WT P450 2A6 and Mutants—To determine which mutation positions were critical in gaining the ability to oxidize 5-OBzl-indole, 13 site-directed mutants were constructed from the Bn1 mutant, with one to four mutation sites recovered (Fig. 1 and TABLE ONE, Bn1/1 to Bn1/1254). The hydroxylation of indole (the native substrate), 4-chloroindole (bearing a small substituent), and 4- and 5-OBzl-indoles (bearing bulky substituents) was measured based on the formation of colored dimers of the substrates (16). The low solubility of the 4- and 5-OBzl-indole substrates in H2O made vigorous shaking a necessity to increase diffusion; therefore, whole cells containing WT and mutant P450 2A6 were used instead of membranes or purified enzymes, which were more fragile under these conditions. The amounts of colored products were calculated based on their extinction coefficients, which were obtained either from the literature (indigo) (35) or from the isolation, identification, and quantification of the products as described under “Characterization of Colored Products Generated from 4-Chloroindole and 4- and 5-OBzlindoles” (see below). For the small substrates indole and 4-chloroindole, all of these mutants (as well as WT P450 2A6) had catalytic activity; the original double mutant (L240C/N297Q) was one of the most efficient in terms of catalytic efficiencies. No dramatic enhancement of activity was found for indole or 4-chloroindole oxidation catalyzed by the new mutants, which was not surprising because the screening method used in this work focused specifically on 5-OBzl-indole biotransformation. The mutations at the 3-position (N297Q) and 4-position (I300V) (Fig. 1) are critical in gaining activity for hydroxylation of bulky substituted substrates (TABLE ONE). Mutants containing both of these mutations (Bn1/1, Bn1/2, Bn1/5, etc.) all showed activity in colored product formation from bulky substituted substrates, whereas mutants lacking either of these (Bn1/3, Bn1/4, Bn1/41, etc.) showed no activity. Thus, a new double mutant, P450 2A6 N297Q/I300V (Bn1/125), was found to be the simplest mutant able to oxidize both 4- and 5-OBzlindoles to form colored products. Only small differences were observed in catalytic efficiencies among these gain-of-activity mutants, although the kcat values do differ, with the Bn1 mutant having the lowest kcat.

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005

Expansion of Substrate Specificity of P450 2A6

FIGURE 3. Biotransformation results of 4- and 5-OBzl-indole oxidations catalyzed by WT P450 2A6 and mutants. A, biotransformation of 4-OBzl(OBn)-indole. HPLC conditions were as follows: 40% (v/v) CH3CN (in H2O) for 10 min, followed by a 40 –90% (v/v) CH3CN linear gradient over 50 min. B, biotransformation of 5-OBzl-indole. HPLC conditions were as follows: 75– 80% (v/v) CH3OH (in H2O with 1% (v/v) CF3CO2H) over 10 min, followed by a 80 –100% CH3OH linear gradient over 20 min and then 100% CH3OH for 10 min. The panels are as follows (from top to bottom): HPLC profiles of substrate and cell extracts from pCW⬘-hor (vector with only human NADPH-cytochrome P450 reductase), from WT P450 2A6, from the P450 2A6 double mutant L240C/N297Q (2A6 DM), and from Bn1/125; enlarged HPLC profile of cell extracts from Bn1/125; and photographs showing the colored cell cultures. mAU, milli-absorbance units.

Steady-state kinetic plots of the biotransformation of indole compounds catalyzed by Bn1/125, Bn1, Bn1/1253, Bn1/1254, the double mutant L240C/N297Q, and WT P450 2A6 are shown in Fig. 2. Convenient as the colorimetric method is, the caveat exists that those mutants that showed no color derived from 4- or 5-OBzl-indole might catalyze the oxidation of the substrates to products without absorbance in the visible range. To address this possibility, HPLC profiles (monitored at 254 nm) were analyzed with crude extracts from the biotransformation of 4- and 5-OBzl-indoles catalyzed by WT P450 2A6, the original double mutant (L240C/N297Q), N297Q/I300V (Bn1/125), and the pCW⬘-hor vector (with only human NADPH-cytochrome P450 reductase) as negative controls (Fig. 3). The expression levels were normal 11 h after induction. In the HPLC profiles, the cell extracts from WT P450 2A6 and P450 2A6 L240C/N297Q showed the same patterns as negative controls, and no products were detected.

DECEMBER 9, 2005 • VOLUME 280 • NUMBER 49

Substrate Binding Affinities of WT P450 2A6 and Mutants—The binding affinities of WT P450 2A6 and some mutants for indole and 5-OBzl-indole were estimated (TABLE TWO). In general, there was no correlation between the catalytic activities and the apparent Kd values. WT P450 2A6 had the lowest catalytic efficiency toward indole among the enzymes tested; however, this enzyme showed the tightest binding of indole (apparent Kd of 13 ␮M). For 5-OBzl-indole, which was the substrate for the Bn1 and Bn1/1 enzymes but not for WT P450 2A6, the Kd values are extremely high for Bn1 (apparent Kd of 710 ␮M) and Bn1/1 (apparent Kd of 414 ␮M), but very low for the WT enzyme (apparent Kd of 8 ␮M). These discrepancies suggest that some of the apparent binding may be nonproductive. Coumarin 7-Hydroxylation Activity—Analysis of the rates of coumarin 7-hydroxylation with membranes prepared from WT P450 2A6 and the double mutants L240C/N297Q and N297Q/I300V (Bn1/125) indi-

JOURNAL OF BIOLOGICAL CHEMISTRY

41095

Expansion of Substrate Specificity of P450 2A6 TABLE TWO

Binding affinities of WT P450 2A6 and mutants for indoles ND, not determined. Kd a

P450 2A6

Indole

5-OBzl-indole ␮M

13 ⫾ 1 44 ⫾ 1 121 ⫾ 6 258 ⫾ 22 143 ⫾ 12

WT L240C/N297Q Bn1 Bn1/1 Bn1/4 a

8⫾1 ND 710 ⫾ 70 414 ⫾ 123 ND

Fitted with the general equation Y ⫽ Bmax ⫻ X/(Kd ⫹ X) (hyperbolic one-site binding) in GraphPad Prism. FIGURE 4. Colored dimers isolated from cell cultures incubated with indole compounds. OBn, OBzl.

TABLE THREE

Kinetic parameters for WT P450 2A6 and mutants in coumarin 7-hydroxylation P450 2A6

WT L240C/N297Q Bn1/125

kcat

Coumarin 7-hydroxylation Km kcat/Km

min⫺1

␮M

min⫺1 ␮M⫺1

7.2 ⫾ 0.2 45 ⫾ 1 34 ⫾ 1

4.4 ⫾ 0.5 19 ⫾ 0.7 18 ⫾ 1.4

1.6 ⫾ 0.2 2.4 ⫾ 0.1 1.9 ⫾ 0.2

cated that the kcat for both mutants was increased ⬃5-fold over that for the WT enzyme, whereas the catalytic efficiency indicated by kcat/Km was not enhanced due to an increase in Km (TABLE THREE). Characterization of Colored Products Generated from 4-Chloroindole and 4- and 5-OBzl-indoles—The colored products from cell cultures fortified with indole, 4-chloroindole, and 4- and 5-OBzl-indoles were analyzed by HPLC at a wavelength of 600 nm. The main product derived from indole was indigo (compound 1), as established by comparison with standard material. 4-Chloroindole yielded 4,4⬘-dichloroindigo (compound 2) as a single peak (Fig. 4). 5-OBzl-indole yielded mainly 5,5⬘-di-OBzl-indirubin (compound 3) (Fig. 4), an expected dimer based on our previous work (18, 19). The main product of 4-OBzl-indole was originally expected to be an indigo because the substituent group at the 4-position makes formation of an indirubin unfavorable due to its interference with the carbonyl group at the 3⬘-position (19, 36). However, the dominant blue dye isolated upon incubations with 4-OBzl-indole was neither indigo nor indirubin, although it had the same molecular mass as that of the indigoids and a UV-visible spectrum similar to that of indigo. 1H NMR clearly showed an unsymmetrical structure with one group of two doublets and one triplet (hydrogens from a 4-substituted benzene ring) plus another group of a singlet and two triplets with small coupling constants (hydrogens at the 2- and 3-positions of the pyrrole ring), suggesting that a pyrrole ring and a benzene ring are connected with a double bond instead of a pyrrole-pyrrole connection. 1H-1H COSY, 1H-1H NOESY, 1H-13C heteronuclear multiple quantum correlation, 1H-13C HMBC, and 1H-13C CIGAR-HMBC NMR experiments (supplemental data, Figs. 1S–5S) were used to assign the structure as 4-OBzl-2-(4⬘-OBzl-1⬘,7⬘-dihydro-7⬘-oxo-6⬘H-indol-6⬘-ylidene)indolin-3-one (compound 4) (Fig. 4). The formation of this compound was as reproducible as that of the other indigoids. The yield of this dye from cultures (13 ␮mol/liter of culture) is at the same level as that of 5,5⬘-di-OBzl-indirubin (9 ␮mol/ liter of culture) by the same mutant, P450 2A6 N297Q/I300V (Bn1/ 125), and similar to that reported in our previous work on 6-nitroindole and methyl indole-5-carboxylate converted by P450 2A6 L240C/

41096 JOURNAL OF BIOLOGICAL CHEMISTRY

N297Q (19). Other gain-of-activity mutants such as Bn1/52, Bn1/51, etc., also produced this dimer as the dominant colored product. The spectroscopic data of the four colored dimers follow. Indigo (compound 1): blue solid; ␭max (Me2SO) 605 nm, ␭max (H2O) 670 nm, ⑀670 ⫽ 3.9 mM⫺1 cm⫺1 (35); 1H NMR (400 MHz, d6-Me2SO) ␦ 6.95 (t, 2H, J ⫽ 7.6 Hz, H-5), 7.33 (d, 2H, J ⫽ 7.6 Hz, H-7), 7.51 (t, 2H, J ⫽ 7.6 Hz, H-6), 7.61 (d, 2H, J ⫽ 7.6 Hz, H-4), and 10.49 (s, 2H, NH); APCI mass spectrometry (MS), m/z 261 ([M ⫺ H]⫺). 4,4⬘-Dichloroindigo (compound 2) (37): blue solid; ␭max (Me2SO) 620 nm, ␭max (H2O) 640 nm, ⑀640 ⫽ 13.2 mM⫺1 cm⫺1; 1H NMR (400 MHz, d6-Me2SO) ␦ 6.93 (d, 2H, J ⫽ 7.2 Hz, H-7), 7.27 (d, 2H, J ⫽ 7.6 Hz, H-5), 7.47 (dd, 2H, J1 ⫽ J2 ⫽ 8.0 Hz, H-6), and 11.94 (br s, 2H, NH); APCI-MS, m/z 330 ([M]⫺). 5,5⬘-DiOBzl-indirubin (compound 3): purple solid; ␭max (Me2SO) 588 nm, ␭max (H2O) 588 nm, ⑀588 ⫽ 11.2 mM⫺1 cm⫺1; 1H NMR (400 MHz, d6-Me2SO) ␦ 5.10 (s, 2H, OCH2Ph), 5.13 (s, 2H, OCH2Ph), 6.79 (d, 1H, J ⫽ 8.4 Hz, H-7⬘), 6.92 (dd, 1H, J1 ⫽ 8.4 Hz, J2 ⫽ 2.8 Hz, H-6⬘), 7.24 (d, 1H, J ⫽ 2.4 Hz, H-4), 7.28 (dd, 1H, J1 ⫽ 8.4 Hz, J2 ⫽ 2.8 Hz, H-6), 7.33 (t, 2H, J ⫽ 7.2 Hz, OCH2Ph), 7.35 (d, 1H, J ⫽ 8.4 Hz, H-7), 7.39 (t, 2H, J ⫽ 7.2 Hz, OCH2Ph), 7.40 (t, 2H, J ⫽ 7.2 Hz, OCH2Ph), 7.45 (d, 2H, J ⫽ 7.2 Hz, OCH2Ph), 7.45 (d, 2H, J ⫽ 7.2 Hz, OCH2Ph), 8.55 (d, 1H, J ⫽ 2.4 Hz, H-4⬘), 10.68 (br s, 1H, 1⬘-NH), and 10.86 (br s, 1H, 1-NH); APCI-MS, m/z 475 ([M ⫹ H]⫹). 4-OBzl-2-(4⬘-OBzl-1⬘,7⬘-dihydro-7⬘-oxo-6⬘H-indol-6⬘-ylidene)indolin-3-one (compound 4): blue solid; ␭max (Me2SO) 622 nm, ␭max (H2O) 670 nm, ⑀670 ⫽ 11.4 mM⫺1 cm⫺1; 1H NMR (500 MHz, d6-Me2SO) ␦ 5.14 (s, 2H, OCH2Ph), 5.26 (s, 2H, OCH2Ph), 6.34 (t, 1H, J ⫽ 2.5 Hz, H-3⬘), 6.69 (d, 1H, J ⫽ 8.0 Hz, H-5), 6.96 (d, 1H, J ⫽ 8.0 Hz, H-7), 7.26 (t, 1H, J ⫽ 2.5 Hz, H-2⬘), 7.45 (dd, 1H, J1 ⫽ J2 ⫽ 8.0 Hz, H-6), 7.35 (t, 2H, J ⫽ 7.3 Hz, OCH2Ph), 7.41–7.44 (m, 4H, OCH2Ph), 7.51 (d, 2H, J ⫽ 7.5 Hz, OCH2Ph), 7.54 (d, 2H, J ⫽ 7.5 Hz, OCH2Ph), 7.64 (s, 1H, H-5⬘), 11.99 (br s, 1H, H-N(1)), and 12.21 (br s, 1H, H-N(1⬘)); 13C NMR (125 MHz, d6-Me2SO) ␦ 69.5 (OCH2Ph), 70.0 (OCH2Ph), 93.6 (C-5⬘), 103.5 (C-3⬘), 106.5 (C-7), 106.7 (C-5), 108.6 (C-3a), 113.9 (C-6⬘), 126.9 (C-3a⬘), 127.9 (C-7a⬘), 128.0 (C-2⬘), 128.5 (OCH2Ph), 128.7 (OCH2Ph), 129.0 (OCH2Ph), 136.9 (C-2), 137.1 (OCH2Ph), 138.9 (C-6), 151.4 (C-4⬘), 152.8 (C-7a), 158.1 (C-4), 176.7 (C-7⬘), and 187.0 (C-3); APCI-MS, m/z 475 ([M ⫹ H]⫹). Identification of Intermediates Generated in Vivo by Recombinant P450 2A6 N297Q/I300V—As shown in Fig. 3, many other products were found in the cell extracts in addition to the colored dimers that have absorbance at 254 nm. Considering the possible interference of other compounds from the E. coli system, liquid chromatography/tandem mass spectrometry (APCI positive mode) was done with the mixtures to help identify peaks that have a fragment at m/z 91 (benzyl group) as evidence of the origin from the substrate 4- or 5-OBzl-indole.

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005

Expansion of Substrate Specificity of P450 2A6

FIGURE 5. Oxidized products isolated from cell cultures incubated with 4- or 5-OBzl-indole. OBn, OBzl.

FIGURE 6. Postulated major reactions in pigment formation derived from 4-OBzl-indole catalyzed by P450 2A6 mutants. OBn, OBzl.

Peaks with qualifying mass spectra and reasonable abundance were collected by preparative HPLC and analyzed in detail. Compounds 5 and 6 were originally eluted together (tR ⫽ 15.4 min), but separated under modified elution conditions (tR ⫽ 24.9 and 26.7 min, respectively) (see “Large Scale Biotransformation and Product Separation” under “Experimental Procedures”). Analysis of the metabolites derived from 4-OBzlindole showed that the benzene ring was hydroxylated at all three positions (positions 5–7). Three phenols, 4-OBzl-1H-indol-5-ol (compound 8), 4-OBzl-1H-indol-6-ol (compound 9), and 4-OBzl-1Hindol-7-ol (compound 10), and one quinone, 4-OBzl-1H-indole-6,7dione (compound 6; which could be the immediate precursor of the final blue dye), were isolated, and their structures were established by MS, 1H NMR, 1H-1H COSY, and 1H-1H NOESY methods (supplemental data, Figs. 6S–9S). Two other products derived from hydroxylation on the pyrrole ring, 4-OBzl-indolin-2-one (compound 7) and 4-OBzlindoline-2,3-dione (compound 5), were also isolated, as expected (Fig. 5). The characterization data for these eight compounds are included in supplemental “Analytical Data for Monomeric Metabolites Derived from 4-OBn Indole and 5-OBn Indole” (where “OBn” is benzyloxy). Although there were several dimers detected by liquid chromatography/tandem mass spectrometry and 1H NMR, the amounts were insufficient for full characterization. 4-OBzl-indolin-2-one (compound 7) and 5-OBzl-indolin-2-one (compound 11) appear to be the primary oxidation products formed from each corresponding substituted indole, consistent with our previous study of indole oxidation (17).

DECEMBER 9, 2005 • VOLUME 280 • NUMBER 49

Postulated Mechanism for 4-OBzl-indole Biotransformation—On the basis of the above results and previous work (17), we propose that the mechanism involved in the formation of compound 4 follows the course shown in Fig. 6. The substrate 4-OBzl-indole is first hydroxylated at one of the five free positions to form compounds 7–10 and an indoxyl (compound 13a), which undergoes further oxidation to form isatin (compound 5) via a tautomer (compound 13b) and an intermediate compound (compound 13c). Accordingly, the primary compounds 9 and 10 can be oxidized to give compound 6, possibly via radical intermediates (compounds 9a and 10a, respectively). Alternatively, oxidation of compounds 9 and 10 to catechols can be proposed, followed by 2-electron oxidation to the o-quinone (compound 6). The final blue dimer (compound 4) then could be the result of combinational reaction of compounds 13c and 6 or compounds 10a and 5. Substrate Docking and Modeling of the P450 2A6 N297Q/I300V Mutant—Analysis of the active site in the crystal structure of the P450 2A6-coumarin complex indicates that coumarin occupies a cavity of ⬃359 Å3. The single mutation I300V, despite elimination of only one methylene group in Val compared with Leu, results in expansion of the active site to an adjacent hydrophobic cavity built by Val110, Ile208, Leu241, and Leu244 and adds an additional 83 Å3 to the total cavity volume, increasing the latter up to 440 Å3 (Fig. 7A). The functionally conservative substitution N297Q does not significantly affect the cavity volume (353 Å3). However, the longer side chain of glutamine presumably enables a new hydrogen bond with the oxygen atom in the substrate

JOURNAL OF BIOLOGICAL CHEMISTRY

41097

Expansion of Substrate Specificity of P450 2A6

FIGURE 7. P450 2A6 active site and putative binding of 4- and 5-OBzl-indoles. A, active-site cavities in coumarin-bound WT P450 2A6 (blue; 359 Å3) and in the I300V mutant (red; 440 Å3). Coumarin is shown is green in the space-filling mode. B, binding of 5-OBzl-indole in the (minimized) active site of the N297Q/I300V mutant. 5-OBzl-indole is shown in green in the space-filling mode. C, schematic representation of putative 4- and 5-OBzl(OBn)-indole orientations in the P450 2A6 binding site. The figure was generated using Tripos SYBYL software.

(20), which may stabilize a certain orientation of the latter in the active site. This proposed electronic-propelled interaction would be as critical as the cavity enlargement effect, and their cooperative contributions result in the gain-of-activity mutant. With both the I300V and N297Q substitutions, the estimated size of the cavity is 460 Å3, about one-third larger than in the native enzyme. Both critical mutation positions Asn297 and Ile300 are on the surface of the I helix facing the substrate-binding site and are separated by one-half of an ␣-helical turn. Orientation of the cavity axis at ⬃60° to the heme plane (Fig. 7A) suggests a likely substrate-binding model in the active site and explains the reaction products obtained upon conversion of both 4- and 5-OBzlindoles. Alignment of the long axis of 5-OBzl-indole with the cavity axis and putative flipping of the compound along this axis would always present the pyrrole ring to the oxygen scission site (Figs. 7C and 8A). In contrast, flipping of 4-OBzl-indole along its long axis would present both pyrrole and benzene rings to the oxygen scission site (Fig. 8, B and C). The representative binding modes for both substrates resulting from automated docking are shown (Fig. 8). In about one-half (15 of 30) of the highest scoring solutions 5-OBzl-indole was bound in the conformation represented by that shown in Fig. 8A, with the pyrrole ring presented to the oxygen scission site. As expected, automated docking of 4-OBzlindole returned a larger variety of docked complexes, with the highest scoring conformers represented in Fig. 8B, where the benzene ring is exposed to the oxygen scission site. An alternative family of conformers featuring the pyrrole ring of 4-OBzl-indole as a target for oxygenation is shown in Fig. 8C.

DISCUSSION We have reported the identification of P450 2A6 enzymes with expanded function using a process of sequential random and site-directed mutagenesis. New chemical reactions were also found for some products formed by the new mutants. These mutations provide for more diversity of product formation. The available x-ray structure of the P450 2A6-coumarin complex (20) has been used to explain the effects of the mutations.

41098 JOURNAL OF BIOLOGICAL CHEMISTRY

Although exhaustive systematic screening of the WT P450 2A6 and P450 2A6 L240C/N297Q (the original double mutant) libraries to cover 95% possible changes (38) was not performed, this strategy used here was successful. Searching for gain-of-activity mutants, for which a low frequency is expected, using the colony-based method described here should be considered more of a selection than a screening process. All blue colonies definitely indicate a gain-of-activity mutant. Any incorrect identification resulted from the uncertainty in reading the color as “blue” due to the disturbance of background color and the brown color of the colonies, which made the second confirmation necessary. From the results of site-directed mutagenesis, we know that one of the critical changes is the N297Q mutation, an existing mutation in the original double mutant enzyme. This mutant (L240C/N297Q) had been selected previously in a molecular breeding experiment (14) that was based on the use of randomized oligonucleotides (39) targeted to the substrate recognition site regions expected to comprise the active site (2). The nucleotide change is from AAC to CAA, a complex double transversion that would not likely be achieved by error-prone PCR. This could be the main reason why no gain-of-activity mutants were found in the WT P450 2A6 library. Analysis of the N297Q and I300V substitutions in silico using the available crystal structure revealed that, although the I300V mutation adds ⬃83 Å3 to the wild-type binding site, resulting in 440 Å3 of total space, it is not enough to accommodate either of the OBzl-substituted indoles used in this study (which would require a volume of ⬃600 Å3). Presumably substrate-induced fit adjustment to the enzyme active site is required to complement the favorable steric and electronic effects of both mutations as deduced from the substrate-docking model. Alternatively, the active site of the N297Q/I300V mutant may have changed enough from that of the WT enzyme to expand the space, although we prefer the former hypothesis regarding induced fit. The compound 4 colored dimer appears to be novel in that no such structure has been reported. In our previous work, ⬃20 indigoids were isolated from P450 2A6-catalyzed reactions (18, 19). All of them have parent rings of indigo or indirubin, in which both positions of the pyr-

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005

Expansion of Substrate Specificity of P450 2A6

FIGURE 8. Stereo views of 4- and 5-OBzl-indoles docked into the active site of the P450 2A6 N297Q/I300V mutant. Substrates are shown in yellow in the ball-and-stick mode. A, binding of 5-OBzl-indole. Coumarin (from the original P450 2A6 crystal structure) is shown in green. B and C, binding of 4-OBzl-indole in two different orientations. The figure was generated using the program SETOR.

role ring (positions 2 and 3) are oxidized and one of the positions is connected to another monomer, leaving the benzene ring intact, consistent with the reactivity of the indole. Indole is a ␲-excessive heteroaromatic, which means that the ␲-electron densities on carbon atoms of the pyrrole ring are greater than those of the benzene, and the pyrrole ring should be more readily oxidized. The most important resonance structures of indole maintain an unbroken benzene ring and bear the positive charge on nitrogen and the negative charge on the C-3 atom. The 3-position of the indole ring is highly reactive toward electrophilic reagents (40). Although it has been reported that indole can be hydroxylated at the 5- and 6-positions in P450-catalyzed oxidation systems, the amounts of these products are limited, and they apparently do not undergo further oxidation to form dimers (17, 41). The intermediate 4-OBzl-1H-indole-6,7-dione (compound 6) does have analogs in nature, e.g. the redox cofactors in methylamine dehydrogenase (tryptophan tryptophylquinone) (42) and amine dehydrogenase (cysteine tryp-

DECEMBER 9, 2005 • VOLUME 280 • NUMBER 49

tophylquinone) (43), which are also 6,7-diones derived from L-tryptophan and formed during post-translational modification. How the two quinone oxygens are incorporated into tryptophan tryptophylquinone is still unclear, although it is believed that other enzymes, probably monooxygenases, are required for the biogenesis instead of a self-processing mechanism (44). Tryptophan tryptophylquinone could be considered as an indole dimer as well in that the pyrrole ring (C-2) of Trp108 is connected to the benzene ring (C-4) of Trp57 with a single bond. Another report on the oxidation of indole derivatives includes structures involving formation of dimers with the benzene-benzene ring connected; in this case, however, the 3-position of pyrrole was protected and the benzene ring was dioxidized in advance (45). The isolation of six intermediates upon incubation with 4-OBzl-indole, especially those derived from hydroxylation on the benzene ring, not only supports the formation of compound 4, but also clearly indicates that the substrate is accommodated inside the active site with

JOURNAL OF BIOLOGICAL CHEMISTRY

41099

Expansion of Substrate Specificity of P450 2A6 variable orientations. On the other hand, the gain-of-activity mutants retained their normal activities (producing indigo or indirubin), like WT P450 2A6, with the smaller indole compounds tested (indole and 4-chloroindole), which are presumed to be more flexible in the active site. The lack of oxidation of the benzene ring of indole by the mutant might be explained partially by the electron-donating property of the OBzl group. The 4-OBzl substituent may activate the para-position to be hydroxylated to compound 8, which contributes to the whole process (Fig. 6). In addition, the rigidity of substrate docking of a bulky substrate compared with a substrate of smaller size may also contribute to the increased chance of hydroxylation on the benzene ring, as long as the substrate can be accommodated properly inside the active site with the benzene ring presented to the active oxygen species. In conclusion, we have isolated new mutants of P450 2A6 with expanded substrate specificity toward bulky substituted indoles from a randomized library constructed by error-prone PCR and identified the critical residues that contribute to the effect by site-directed mutagenesis. The products (dimers and monomeric oxidation products) were characterized, and we proposed a rational mechanism for the formation of an unusual dimer derived from 4-OBzl-indole. We interpreted the results in the context of a crystal structure of a P450 2A6-coumarin complex. The same approach, together with the analysis of the crystal structure, could be adapted to generate more diversity of products from the P450 2A6-catalyzed metabolism and could be used as a tool in the development of P450 enzymes with new functions. Acknowledgments—We thank E. F. Johnson and J. K. Yano for sharing the details of the P450 2A6 crystal structure in advance of publication; E. Dawson and the Vanderbilt Center for Structural Biology for consultation, training, and access to facilities for the computational modeling and analysis; M. Voehler and D. Stec for assistance in the NMR analysis; and K. Trisler for assistance in preparation of the manuscript.

REFERENCES 1. Palmer, G., and Reedijk, J. (1992) J. Biol. Chem. 267, 665– 677 2. Gotoh, O. (1992) J. Biol. Chem. 267, 83–90 3. Zhao, H., Giver, L., Shao, Z., Affholter, J. A., and Arnold, F. H. (1998) Nat. Biotechnol. 16, 258 –261 4. Ortiz de Montellano, P. R. (ed) (2005) Cytochrome P450: Structure, Mechanism, and Biochemistry, 3rd Ed., Kluwer Academic/Plenum Publishers, New York 5. Guengerich, F. P. (2001) Chem. Res. Toxicol. 14, 611– 650 6. Guengerich, F. P. (2005) in Cytochrome P450: Structure, Mechanism, and Biochemistry (Ortiz de Montellano, P. R., ed) 3rd Ed., pp. 377–531, Kluwer Academic/Plenum Publishers, New York 7. Guengerich, F. P. (2002) Nat. Rev. Drug Discov. 1, 359 –366 8. Arnold, F. H., and Georgiou, G. (2003) Methods Mol. Biol. 230 9. Peters, M. W., Meinhold, P., Glieder, A., and Arnold, F. H. (2003) J. Am. Chem. Soc. 125, 13442–13450 10. Seng Wong, T., Arnold, F. H., and Schwaneberg, U. (2004) Biotechnol. Bioeng. 85, 351–358 11. Parikh, A., Josephy, P. D., and Guengerich, F. P. (1999) Biochemistry 38, 5283–5289 12. Yun, C.-H., Miller, G. P., and Guengerich, F. P. (2000) Biochemistry 39, 11319 –11329

41100 JOURNAL OF BIOLOGICAL CHEMISTRY

13. Kim, D., and Guengerich, F. P. (2004) Biochemistry 43, 981–988 14. Nakamura, K., Martin, M. V., and Guengerich, F. P. (2001) Arch. Biochem. Biophys. 395, 25–31 15. Stemmer, W. P. (1994) Nature 370, 389 –391 16. Gillam, E. M. J., Aguinaldo, A. M., Notley, L. M., Kim, D., Mundkowski, R. G., Volkov, A., Arnold, F. H., Soucek, P., DeVoss, J. J., and Guengerich, F. P. (1999) Biochem. Biophys. Res. Commun. 265, 469 – 472 17. Gillam, E. M. J., Notley, L. M., Cai, H., DeVoss, J. J., and Guengerich, F. P. (2000) Biochemistry 39, 13817–13824 18. Guengerich, F. P., Sorrells, J. L., Schmitt, S., Krauser, J. A., Aryal, P., and Meijer, L. (2004) J. Med. Chem. 43, 3236 –3241 19. Wu, Z.-L, Aryal, P., Lozach, O., Meijer, L., and Guengerich, F. P. (2005) Chem. Biodivers. 2, 51– 65 20. Yano, J. K., Hsu, M. H., Griffin, K. J., Stout, C. D., and Johnson, E. F. (2005) Nat. Struct. Mol. Biol. 12, 822– 823 21. Schoch, G. A., Yano, J. K., Wester, M. R., Griffin, K. J., Stout, C. D., and Johnson, E. F. (2004) J. Biol. Chem. 279, 9497–9503 22. Yun, C.-H., Shimada, T., and Guengerich, F. P. (1991) Mol. Pharmacol. 40, 679 – 685 23. Yamazaki, H., Inoue, K., Hashimoto, M., and Shimada, T. (1999) Arch. Toxicol. 73, 65–70 24. Hecht, S. S., Hochalter, J. B., Villalta, P. W., and Murphy, S. E. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 12493–12497 25. Soucek, P. (1999) Arch. Biochem. Biophys. 370, 190 –200 26. Parikh, A., Gillam, E. M. J., and Guengerich, F. P. (1997) Nat. Biotechnol. 15, 784 –788 27. Sandhu, P., Baba, T., and Guengerich, F. P. (1993) Arch. Biochem. Biophys. 306, 443– 450 28. Sambrook, J., and Russell, D. W. (2001) Molecular Cloning: A Laboratory Manual, 3rd Ed., Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 29. Schenkman, J. B., Remmer, H., and Estabrook, R. W. (1967) Mol. Pharmacol. 3, 113–123 30. Guengerich, F. P. (2001) in Principles and Methods of Toxicology (Hayes, A. W., ed) 4th Ed., pp. 1625–1687, Taylor & Francis, Philadelphia 31. Jones, T. A., Zou, J. Y., Cowan, S. W., and Kjeldgaard, M. (1991) Acta Crystallogr. Sect. A 47, 110 –119 32. Brunger, A. T., Adams, P. D., Clore, G. M., DeLano, W. L., Gros, P., Grosse-Kunstleve, R. W., Jiang, J. S., Kuszewski, J., Nilges, M., Pannu, N. S., Read, R. J., Rice, L. M., Simonson, T., and Warren, G. L. (1998) Acta Crystallogr. Sect. D Biol. Crystallogr. 54, 905–921 33. Rarey, M., Kramer, B., Lengauer, T., and Klebe, G. (1996) J. Mol. Biol. 261, 470 – 489 34. Yu, D., Ellis, H. M., Lee, E. C., Jenkins, N. A., Copeland, N. G., and Court, D. L. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 5978 –5983 35. Li, Q. S., Schwaneberg, U., Fischer, P., and Schmid, R. D. (2000) Chem. Eur. J. 6, 1531–1536 36. Polychronopoulos, P., Magiatis, P., Skaltsounis, A. L., Myrianthopoulos, V., Mikros, E., Tarricone, A., Musacchio, A., Roe, S. M., Pearl, L., Leost, M., Greengard, P., and Meijer, L. (2004) J. Med. Chem. 47, 935–946 37. Weinstein, J., and Wyman, G. M. (1956) J. Am. Chem. Soc. 78, 2387–2390 38. Cirino, P. C., Mayer, K. M., and Umeno, D. (2003) Methods Mol. Biol. 231, 3–9 39. Parikh, A., and Guengerich, F. P. (1998) BioTechniques 24, 428 – 431 40. Remers, W. A., Houlihan, W. J., and Brown, R. K. (ed) (1971) The Chemistry of Heterocyclic Compounds, Vol. 25, Wiley-Interscience, New York 41. King, L. J., Parke, D. V., and Williams, R. T. (1966) Biochem. J. 98, 266 –277 42. McIntire, W. S. (1992) Essays Biochem. 27, 119 –134 43. Datta, S., Mori, Y., Takagi, K., Kawaguchi, K., Chen, Z. W., Okajima, T., Kuroda, S., Ikeda, T., Kano, K., Tanizawa, K., and Mathews, F. S. (2001) Proc. Natl. Acad. Sci. U. S. A. 98, 14268 –14273 44. Pearson, A. R., De La Mora-Rey, T., Graichen, M. E., Wang, Y., Jones, L. H., Marimanikkupam, S., Agger, S. A., Grimsrud, P. A., Davidson, V. L., and Wilmot, C. M. (2004) Biochemistry 43, 5494 –5502 45. Tabatabaie, T., and Dryhurst, G. (1992) J. Med. Chem. 35, 2261–2274

VOLUME 280 • NUMBER 49 • DECEMBER 9, 2005