Four Decades of Cytochrome P450 2B Research: From Protein Adducts to Protein Structures and Beyond.

This article features selected findings from the senior author and colleagues dating back to 1978 and covering approximately three-fourths of the 60 years since the discovery of cytochrome P450. Considering the vast number of P450 enzymes in this amazing superfamily and their importance for so many fields of science and medicine, including drug design and development, drug therapy, environmental health, and biotechnology, a comprehensive review of even a single topic is daunting. To make a meaningful contribution to the 50th anniversary of Drug Metabolism and Disposition, we trace the development of the research in a single P450 laboratory through the eyes of seven individuals with different backgrounds, perspectives, and subsequent career trajectories. All co-authors are united in their fascination for the structural basis of mammalian P450 substrate and inhibitor selectivity and using such information to improve drug design and therapy. An underlying theme is how technological advances enable scientific discoveries that were impossible and even inconceivable to prior generations. The work performed spans the continuum from: 1) purification of P450 enzymes from animal tissues to purification of expressed human P450 enzymes and their site-directed mutants from bacteria; 2) inhibition, metabolism, and spectral studies to isothermal titration calorimetry, deuterium exchange mass spectrometry, and NMR; 3) homology models based on bacterial P450 X-ray crystal structures to rabbit and human P450 structures in complex with a wide variety of ligands. Our hope is that humanizing the scientific endeavor will encourage new generations of scientists to make fundamental new discoveries in the P450 field. SIGNIFICANCE STATEMENT: The manuscript summarizes four decades of work from Dr. James Halpert's laboratory, whose investigations have shaped the cytochrome P450 field, and provides insightful perspectives of the co-authors. This work will also inspire future drug metabolism scientists to make critical new discoveries in the cytochrome P450 field.

Increased Phenacetin Oxidation upon the L382V Substitution in Cytochrome P450 1A2 is Associated with Altered Substrate Binding Orientation.

Leucine382 of cytochrome P450 1A2 (CYP1A2) plays an important role in binding and O-dealkylation of phenacetin, with the L382V mutation increasing substrate oxidation (Huang and Szklarz, 2010, Drug Metab. Dispos. 38:1039⁻1045). This was attributed to altered substrate binding orientation, but no direct experimental evidence had been available. Therefore, in the current studies, we employed nuclear magnetic resonance (NMR) longitudinal (T1) relaxation measurements to investigate phenacetin binding orientations within the active site of CYP1A2 wild type (WT) and mutants. Paramagnetic relaxation time (T1P) for each proton of phenacetin was calculated from the T1 value obtained from the enzymes in ferric and ferrous-CO state in the presence of phenacetin, and used to model the orientation of phenacetin in the active site. All aromatic protons of phenacetin were nearly equidistant from the heme iron (6.34⁻8.03 Å). In contrast, the distance between the proton of the ⁻OCH2⁻ group, which is abstracted during phenacetin oxidation, and the heme iron, was much shorter in the L382V (5.93 Å) and L382V/N312L (5.96 Å) mutants compared to the N312L mutant (7.84 Å) and the wild type enzyme (6.55 Å), consistent with modeling results. These studies provide direct evidence for the molecular mechanism underlying increased oxidation of phenacetin upon the L382V mutation.

Rational Re-Engineering of the O-Dealkylation of 7-Alkoxycoumarin Derivatives by Cytochromes P450 2B from the Desert Woodrat Neotoma lepida.

On the basis of recent functional and structural characterization of cytochromes P450 2B from the desert woodrat (Neotoma lepida), the 7-alkoxycoumarin and 7-alkoxy-4-(trifluoromethyl)coumarin O-dealkylation profiles of CYP2B35 and CYP2B37 were re-engineered. Point mutants interchanging residues at seven positions in the enzyme active sites were created and purified from an Escherichia coli expression system. In screens for O-dealkylation activity, wild-type CYP2B35 metabolized long-chain 7-alkoxycoumarins but not 7-alkoxy-4-(trifluoromethyl)coumarins or short-chain 7-alkoxycoumarins. Wild-type CYP2B37 metabolized short-chain substrates from both series of compounds. CYP2B35 A367V showed maximal activity with 7-butoxycoumarin as opposed to 7-heptoxycoumarin in the parental enzyme, and CYP2B35 A363I/A367V produced an activity profile like that generated by CYP2B37. CYP2B35 A363I/A367V/I477F showed 7-ethoxycoumarin and 7-ethoxy-4-(trifluoromethyl)coumarin O-dealkylation rates similar to those of CYP2B37 and higher than those of the double mutant. A CYP2B35 septuple mutant retained a CYP2B37-like activity profile. In contrast, the CYP2B37 septuple mutant produced very low rates of O-dealkylation of all substrates. As mutating residue 108 in either enzyme was detrimental, this change was removed from both septuple mutants. Remarkably, the CYP2B35 sextuple mutant produced an activity profile that was a hybrid of that of CYP2B35 and CYP2B37, whereas the CYP2B37 sextuple mutant had almost no O-dealkylation activity. Docking of 7-substituted coumarin derivatives into a model of the CYP2B35 sextuple mutant based on a previous crystal structure of the 4-(4-chlorophenyl)imidazole wild-type complex revealed how the mutant can exhibit activities of both CYP2B35 and CYP2B37.

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

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

Human cytochrome P450 1A1 structure and utility in understanding drug and xenobiotic metabolism.

Cytochrome P450 (CYP) 1A1 is an extrahepatic monooxygenase involved in the metabolism of endogenous substrates and drugs, as well as the activation of certain toxins and environmental pollutants. CYP1A1 is particularly well known for its ability to biotransform polycyclic aromatic hydrocarbons, such as benzo[a]pyrene in tobacco smoke, into carcinogens. CYP1A1 possesses functional similarities and differences with human CYP1A2 and CYP1B1 enzymes, but the structural basis for this has been unclear. We determined a 2.6 Å structure of human CYP1A1 with the inhibitor α-naphthoflavone. α-Naphthoflavone binds within an enclosed active site, with the planar benzochromen-4-one core packed flat against the I helix that composes one wall of the active site, and the 2-phenyl substituent oriented toward the catalytic heme iron. Comparisons with previously determined structures of the related cytochrome P450 1A2 and 1B1 enzymes reveal distinct features among the active sites that may underlie the functional variability of these enzymes. Finally, docking studies probed the ability of CYP1A structures to assist in understanding their known in vitro interactions with several typical substrates and inhibitors.

The aryl hydrocarbon receptor interacts with nuclear factor erythroid 2-related factor 2 to mediate induction of NAD(P)H:quinoneoxidoreductase 1 by 2,3,7,8-tetrachlorodibenzo-p-dioxin.

NAD(P)H:quinoneoxidoreductase 1 (NQO1) belongs to a group of the aryl hydrocarbon receptor (AhR) battery of drug-metabolizing enzymes that are characteristically induced by both AhR agonists and nuclear factor erythroid 2-related factor 2 (Nrf2) activators. We have previously reported that induction of Nqo1 by the AhR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in hepa1c1c7 cells involves Nrf2 (Ma et al., Biochem J 377, 205-213, 2004). Here we analyzed the molecular mechanism of induction. Induction required AhR and its DNA-binding partner Arnt because induction was not observed in AhR or Arnt-defective cells, but induction was restored upon reconstitution of the variant cells with functional AhR or Arnt. Induction also required Nrf2, as induction by benzo[a]pyrene was lost in the liver of Nrf2 knockout mice similarly to induction by butyl hydroxyanisol, demonstrating a cross-interaction between the AhR and Nrf2 pathways for induction in vivo. TCDD increased the protein level and induced the nuclear accumulation of Nrf2 with a delayed kinetics compared with activation of AhR. Chromatin immunoprecipitation revealed that TCDD recruited both AhR and Nrf2 to the Nqo1 promoter enhancer region containing a DRE and an ARE in time-dependent manners. Co-immunoprecipitation experiments revealed that, in addition to AhR-Arnt binding, TCDD induced an interaction between AhR and Nrf2 as well as Keap1. The findings reveal that TCDD induces multi protein complexes to mediate cross-interaction between the AhR and Nrf2 pathways, uncovering a novel mechanistic aspect of gene regulation by environmental chemicals through AhR and Nrf2.

Preferred binding orientations of phenacetin in CYP1A1 and CYP1A2 are associated with isoform-selective metabolism.

Human cytochromes P450 1A1 and 1A2 play important roles in drug metabolism and chemical carcinogenesis. Although these two enzymes share high sequence identity, they display different substrate specificities and inhibitor susceptibilities. In the present studies, we investigated the structural basis for these differences with phenacetin as a probe using a number of complementary approaches, such as enzyme kinetics, stoichiometric assays, NMR, and molecular modeling. Kinetic and stoichiometric analyses revealed that substrate specificity (k(cat)/K(m)) of CYP1A2 was approximately 18-fold greater than that of CYP1A1, as expected. Moreover, despite higher H2O2 production, the coupling efficiency of reducing equivalents to acetaminophen formation in CYP1A2 was tighter than that in CYP1A1. CYP1A1, in contrast to CYP1A2, displayed much higher uncoupling, producing more water. The subsequent NMR longitudinal (T1) relaxation studies with the substrate phenacetin and its product acetaminophen showed that both compounds displayed similar binding orientations within the active site of CYP1A1 and CYP1A2. However, the distance between the OCH2 protons of the ethoxy group (site of phenacetin O-deethylation) and the heme iron was 1.5 Å shorter in CYP1A2 than in CYP1A1. The NMR findings are thus consistent with our kinetic and stoichiometric results, providing a likely molecular basis for more efficient metabolism of phenacetin by CYP1A2.

Differential oxidation of two thiophene-containing regioisomers to reactive metabolites by cytochrome P450 2C9.

The uricosuric diuretic agent tienilic acid (TA) is a thiophene-containing compound that is metabolized by P450 2C9 to 5-OH-TA. A reactive metabolite of TA also forms a covalent adduct to P450 2C9 that inactivates the enzyme and initiates immune-mediated hepatic injury in humans, purportedly through a thiophene-S-oxide intermediate. The 3-thenoyl regioisomer of TA, tienilic acid isomer (TAI), is chemically very similar and is reported to be oxidized by P450 2C9 to a thiophene-S-oxide, yet it is not a mechanism-based inactivator (MBI) of P450 2C9 and is reported to be an intrinsic hepatotoxin in rats. The goal of the work presented in this article was to identify the reactive metabolites of TA and TAI by the characterization of products derived from P450 2C9-mediated oxidation. In addition, in silico approaches were used to better understand both the mechanisms of oxidation of TA and TAI and/or the structural rearrangements of oxidized thiophene compounds. Incubation of TA with P450 2C9 and NADPH yielded the well-characterized 5-OH-TA metabolite as the major product. However, contrary to previous reports, it was found that TAI was oxidized to two different types of reactive intermediates that ultimately lead to two types of products, a pair of hydroxythiophene/thiolactone tautomers and an S-oxide dimer. Both TA and TAI incorporated ¹8O from ¹8O2 into their respective hydroxythiophene/thiolactone metabolites indicating that these products are derived from an arene oxide pathway. Intrinsic reaction coordinate calculations of the rearrangement reactions of the model compound 2-acetylthiophene-S-oxide showed that a 1,5-oxygen migration mechanism is energetically unfavorable and does not yield the 5-OH product but instead yields a six-membered oxathiine ring. Therefore, arene oxide formation and subsequent NIH-shift rearrangement remains the favored mechanism for formation of 5-OH-TA. This also implicates the arene oxide as the initiating factor in TA induced liver injury via covalent modification of P450 2C9. Finally, in silico modeling of P450 2C9 active site ligand interactions with TA using the catalytically active iron-oxo species revealed significant differences in the orientations of TA and TAI in the active site, which correlated well with experimental results showing that TA was oxidized only to a ring carbon hydroxylated product, whereas TAI formed both ring carbon hydroxylated products and an S-oxide.

Significant increase in phenacetin oxidation on L382V substitution in human cytochrome P450 1A2.

Human CYP1A2 is an important drug-metabolizing enzyme, similar in sequence to CYP1A1 but with distinct substrate specificity. We have previously shown that residue 382 affected CYP1A1 and CYP1A2 specificities with alkoxyresorufins. To determine whether this residue is also important for the metabolism of other substrates, we have investigated phenacetin oxidation by single (T124S, T223N, V227G, N312L, and L382V) and multiple (L382V/T223N, L382V/N312L, L382V/T223N/N312L, and L382V/T124S/N312L) mutants of CYP1A2. The enzymes were expressed in Escherichia coli and purified. All the CYP1A2 mutants that contained the L382V substitution displayed much higher activities than the wild-type enzyme, with k(cat) values 3-fold higher, in contrast to other mutants, for which k(cat) decreased. Likewise, a significant increase in specificity, expressed as the k(cat)/K(m) ratio, was observed for the mutants containing the L382V substitution. The efficiency of coupling of reducing equivalents to acetaminophen formation was decreased for all the single mutants except L382V, for which the coupling increased. This effect was also observed with multiple CYP1A2 mutants containing the L382V substitution. Low activities of the four other single mutants were likely caused by dramatically increased uncoupling to water. In contrast, the increase in activity of the L382V-containing mutants resulted from decreased water formation. This finding is consistent with molecular dynamics results, which showed decreased phenacetin mobility leading to increased product formation. The results of these studies confirm the importance of residue 382 in CYP1A2-catalyzed oxidations and show that a single residue substitution can dramatically affect enzymatic activity.

Application of molecular modeling for prediction of substrate specificity in cytochrome P450 1A2 mutants.

Molecular dynamics (MD) simulations of 7-ethoxy- and 7-methoxyresorufin bound in the active site of cytochrome P450 (P450) 1A2 wild-type and various mutants were used to predict changes in substrate specificity of the mutants. A total of 26 multiple mutants representing all possible combinations of five key amino acid residues, which are different between P450 1A1 and 1A2, were examined. The resorufin substrates were docked in the active site of each enzyme in the productive binding orientation, and MD simulations were performed on the enzyme-substrate complex. Ensembles collected from MD trajectories were then scored on the basis of geometric parameters relating substrate position with respect to the activated oxoheme cofactor. The results showed a high correlation between the previous experimental data on P450 1A2 wild-type and single mutants with respect to the ratio between 7-ethoxyresorufin-O-deethylase (EROD) and 7-methoxyresorufin-O-demethylase (MROD) activities and the equivalent in silico "E/M scores" (the ratio of hits obtained with 7-ethoxyresorufin to those obtained with 7-methoxyresorufin). Moreover, this correlation served to establish linear regression models used to evaluate E/M scores of multiple P450 1A2 mutants. Seven mutants, all of them incorporating the L382V substitution, were predicted to shift specificity to that of P450 1A1. The predictions were then verified experimentally. The appropriate P450 1A2 multiple mutants were constructed by site-directed mutagenesis, expressed in Escherichia coli, and assayed for EROD and MROD activities. Of six mutants, five demonstrated an increased EROD/MROD ratio, confirming modeling predictions.

An Effective Approach to Teaching Pharmacogenomics in the First Year of Pharmacy Curriculum.

Objective. To develop an effective method in teaching pharmacogenomics as a part of a new course, Biopharmaceutics and Pharmacogenomics. Methods. Teaching effectiveness was measured by quizzes, retrospective pre- and post-surveys, team activities, and journal reflections. Four team activities were included in the course: genomic disease, patient case, genetic counselor and a debate about personalized medicine. Outcomes and course impact were evaluated at the end of the course. The evaluation methods included the assessment of knowledge, students' perceptions regarding the utility of team activities, the impact of the course on students' confidence to discuss pharmacogenomics with health care providers or patients, and long-term knowledge retention, measured in the following P2 semester. Results. Seventy-six students were enrolled in the course. Multiple assessments during the course demonstrated that students' knowledge of pharmacogenomics improved. The team activities had a positive impact on student learning, and the course improved their confidence level to discuss pharmacogenomics with another health care provider or a patient. While 86% of the students considered themselves "unconfident," "somewhat unconfident" or "neither confident nor unconfident" at the beginning of the course, 91% reported being "confident" or "somewhat confident" by the end of the course. This increase in confidence was statistically significant. Furthermore, students showed knowledge retention six months after taking the course. Conclusion. Implementation of a new course in pharmacogenomics was effective and well received by the students. It also prepared students for system-based therapeutics courses later in the curriculum.

Arachidonic and eicosapentaenoic acid metabolism by human CYP1A1: highly stereoselective formation of 17(R),18(S)-epoxyeicosatetraenoic acid.

Human cytochrome P450 1A1 (CYP1A1) and human NADPH-cytochrome P450 reductase were expressed and purified from Spodoptera frugiperda insect cells. A reconstituted enzymatically active system metabolized polyunsaturated fatty acids such as arachidonic (AA) and eicosapentaenoic acid (EPA). CYP1A1 was an AA hydroxylase which oxidizes this substrate at a rate of 650+/-10 pmol/min/nmol CYP1A1, with over 90% of metabolites accounted for by hydroxylation products and with 19-OH-AA as major product. Epoxyeicosatrienoic acid (EET), mainly 14,15-EET, accounted for about 7% of total metabolites. Unlike rat CYP1A1, the human enzyme exhibited no 20-OH-AA as product. In contrast, with EPA as substrate CYP1A1 was mainly an epoxygenase, oxidizing with over 68% of total metabolites EPA to 17(R),18(S)-epoxyeicosatetraenoic acid (17(R),18(S)-EETeTr). 19-OH-EPA accounted for about 31% of total metabolites. Significantly, the 17,18-olefinic bond of EPA was epoxidized to 17(R),18(S)-EETeTr with nearly absolute regio- and stereoselectivity. Molecular modeling analyses provided rationale for high efficiency of AA hydroxylation at C(19) and its gradual decrease down to C(14), as well as for the limited EPA 17(S),18(R) epoxidation due to unfavorable enzyme-substrate interactions. The absence of omega-hydroxylation for both substrates is not due to steric factors, but probably a consequence of different reactivities of omega and (omega-1) carbons for hydrogen abstraction. It is suggested that the capacity of human CYP1A1 to metabolize AA and EPA and its inducibility by polycyclic aromatic hydrocarbons may affect the production of physiologically active metabolites, in particular, in the cardiovascular system and other extrahepatic tissues including lung.

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

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

Application of 3-Dimensional Homology Modeling of Cytochrome P450 2B1 for Interpretation of Site-Directed Mutagenesis Results.

Abstract Three-dimensional structures of cytochrome P450 2B1 were modeled based on the crystallographic structure of P450(cam). The effect of the alignment, loop choice, and minimization with or without water was assessed. Although final models were similar in overall structure, the identity of active site residues depended upon the alignment. An example is Phe-206, which may or may not form part of the active site. The choice of the loop conformation had a lesser effect, while including water in the final minimization step was essential for preserving the shape and size of the active site. The best model (model 2) was in good agreement with the data from site-directed mutagenesis studies, and correctly predicted the effect of substitutions at 9 out of 10 amino acid positions. Thus, residues important for P450 2B1 activity, such as Ile- 114, Phe-206, Ile-290, Thr-302, Val-363, and Gly-478, constitute part of the active site and are able to interact with the substrate androstenedione through hydrophobic interactions. On the other hand, Ser-303, Ser-360 and Lys-473 are far from the active site and/or cannot interact with the substrate, in agreement with experimental data. The model indicates other residues likely to be important for enzyme function, such as Tyr- 111, Leu-209, Ile-477, and Ile- 480, which can be tested experimentally. The substrate may assume numerous binding orientations consistent with observed patterns of hydroxylation at C(5) and C(6). The replacement in the model of certain amino acid residues to mimic residue substitutions from site-directed mutagenesis studies and docking of the substrate into the modified active site allowed a plausible explanation for alterations in regio- and stereospecificities of some mutants of P450 2B1, such as Gly-478 → Ala or Val-363 Ala.

Characterization of substrate binding to cytochrome P450 1A1 using molecular modeling and kinetic analyses: case of residue 382.

Key residue Val-382 in P450 1A1 has been predicted to interact with the alkoxy chain of resorufin derivatives. Therefore, we undertook a detailed analysis of substrate mobility in the active site of the P450 1A1 homology model and assessed the effect of mutations at position 382. Dynamic trajectories of 7-methoxy-, 7-ethoxy-, and 7-pentoxyresorufin indicated that 7-ethoxyresorufin would be oxidized most efficiently by the wild-type enzyme. The Val-382-->Ala mutation would increase the O-dealkylation of 7-pentoxyresorufin but decrease the oxidation of other substrates. In the case of the V382L mutant, the large bulk of Leu would block alkoxyresorufins from productive binding orientations leading to lowered activities. Binding free energy calculations for three substrates with 1A1 WT and two mutants indicated that binding constants would be similar for all enzyme-substrate combinations. Modeling predictions were tested experimentally. The plasmid containing the cDNA for human P450 1A1 modified for bacterial expression was altered to include a C-terminal PCR-generated six histidine domain to facilitate enzyme purification. The V382A and V382L mutants were constructed by site-directed mutagenesis and Escherichia coli-expressed enzymes purified using Ni-NTA affinity chromatography. The activity of the WT 1A1 was highest toward 7-ethoxyresorufin and lowest toward 7-pentoxyresorufin. Both mutants displayed a decrease in V(max) with 7-methoxy- and 7-ethoxyresorufin, whereas for the V382A mutant, V(max) with 7-pentoxyresorufin was increased. No significant changes in K(m) were observed relative to the wild-type enzyme. The experimental results are thus in good agreement with modeling predictions.

The effect of reciprocal active site mutations in human cytochromes P450 1A1 and 1A2 on alkoxyresorufin metabolism.

Five reciprocal active site mutants of P450 1A1 and 1A2 and an additional mutant, Val/Leu-382 --> Ala, were constructed, expressed in Escherichia coli, and purified by Ni-NTA affinity chromatography. In nearly every case, the residue replacement led to loss of 7-methoxy- and 7-ethoxyresorufin O-dealkylase activity compared to the wild-type enzymes, except for the P450 1A1 S122T mutation which increased both activities. Mutations at position 382 in both P450 1A1 and 1A2 shifted substrate specificity from one enzyme to another, confirming the importance of this residue. Changes in activity of P450 1A enzymes upon amino acid replacement were, in general, consistent with molecular dynamics analyses of substrate motion in the active site of homology models.