Inhibition of human cytochrome P450 1A2 by flavones: A molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone >3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.     

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.     

Molecular Modeling of Cytochrome P450 1A1: Enzyme-Substrate Interactions and Substrate Binding Affinities

Abstract

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 methoxyre- sorufin, 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.     

Tight-binding inhibition by alpha-naphthoflavone of human cytochrome P450 1A2

Human cytochrome P450 (P450) enzymes exhibit remarkable diversity in their substrate specificities, participating in oxidation reactions of a wide range of xenobiotic drugs. Previously, we reported that alpha-naphthoflavone (ANF) is bound to the recombinant P450 1A2 tightly and stabilizes an overall enzyme conformation. The present study is designed to determine the type of P450 1A2 inhibition exerted by ANF, using two different substrates of P450 1A2, 7-ethoxycoumarin (EOC) and 7-ethoxyresorufin (EOR). ANF is generally known as a competitive inhibitor of the enzyme. However, in our tight-binding enzyme kinetics study, ANF acts as noncompetitive inhibitor in 7-ethoxycoumarin O-deethylation (ECOD) (K(i)=55.0 nM), but as competitive inhibitor in 7-ethoxyresorufin O-deethylation (EROD) (K(i)=1.4 nM). Based on homology modeling studies, ANF is positioned to bind to a hydrophobic cavity next to the active site where it may cause a direct effect on substrate binding. It is agreed with the predicted binding site of ANF in P450 3A4, in which ANF is rather known as a stimulating modulator. Our results suggest that ANF binds near the active site of P450 1A2 and exhibits differential inhibition mechanisms, possibly depending on the molecular structure of the substrate.     

Flavones and flavone synthases

Within the secondary metabolite class of flavonoids which consist of more than 9000 known structures, flavones define one of the largest subgroups. Their natural distribution is demonstrated for almost all plant tissues. Various flavone aglyca and their O- or C-glycosides have been described in the literature. The diverse functions of flavones in plants as well as their various roles in the interaction with other organisms offer many potential applications, not only in plant breeding but also in ecology, agriculture and human nutrition and pharmacology. In this context, the antioxidative activity of flavones, their use in cancer prevention and treatment as well as the prevention of coronary heart disease should be emphasized. The therapeutic potential of flavones makes these compounds valuable targets for drug design, including recombinant DNA approaches. The biosynthesis of flavones in plants was found to be catalyzed by two completely different flavone synthase proteins (FNS), a unique feature within the flavonoids. The first, FNS I, a soluble dioxygenase, was only described for members of the Apiaceae family so far. The second, FNS II, a membrane bound cytochrome P450 enzyme, has been found in all other flavone accumulating tissues. This phenomenon is particularly of interest from the evolutionary point of view concerning the flavone biosynthesis and functions in plants. Recently, FNS I and FNS II genes have been cloned from a number of plant species. This now enables detailed biochemical and molecular characterizations and also the development of direct metabolic engineering strategies for modifications of flavone synthesis in plants to improve their nutritional and/or biopharmaceutical value.     

Oxidation of human cytochrome P450 1A2 substrates by Bacillus megaterium cytochrome P450 BM3

Cytochrome P450 enzymes (P450s or CYPs) are good candidates for biocatalysis in the production of fine chemicals, including pharmaceuticals. Despite the potential use of mammalian P450s in various fields of biotechnology, these enzymes are not suitable as biocatalysts due to their low stability, low catalytic activity, and limited availability. Recently, wild-type and mutant forms of bacterial P450 BM3 (CYP102A1) from Bacillus megaterium have been found to metabolize various. It has therefore been suggested that CYP102A1 may be used to generate the metabolites of drugs and drug candidates. In this report, we show that the oxidation reactions of typical human CYP1A2 substrates (phenacetin, ethoxyresorufin, and methoxyresorufin) are catalyzed by both wild-type and mutant forms of CYP102A1. In the case of phenacetin, CYP102A1 enzymes show only O-deethylation product, even though two major products are produced as a result of O-deethylation and 3-hydroxylation reactions by human CYP1A2. Formation of the metabolites was confirmed by HPLC analysis and LC–MS to compare the metabolites with the actual biological metabolites produced by human CYP1A2. The results demonstrate that CYP102A1 mutants can be used for cost-effective and scalable production of human CYP1A2 drug metabolites. Our computational findings suggest that a conformational change in the cavity size of the active sites of the mutants is dependent on activity change. The modeling results further suggest that the activity change results from the movement of several specific residues in the active sites of the mutants.     

Aromatase inhibition by flavonoids

Several synthetic flavones were found to inhibit the aromatization of androstenedione to estrone catalyzed by human placental microsomes. Twenty-one compounds were tested and the IC50 of the most active were: flavone, 10 microM; 7-hydroxyflavone, 0.5 microM; 7,4'-dihydroxyflavone, 2.0 microM; flavanone, 8.0 microM; and 4'-hydroxyflavanone, 10 microM. Most of the others had IC50 values ranging from 80 to greater than 200 microM. These findings show that 4'-hydroxylation results in either no change or very little change in IC50 for flavanone, isoflavone and isoflavanone as well as other ring A hydroxylated flavones. Derivatives of flavone with a hydroxyl substituent at position 5, 6 and 7 were also screened. 7-Hydroxyflavone (11) was the most effective competitive inhibitor (IC50 = 0.5 microM) with an apparent Ki value of 0.25 microM. Compound 11 also induced a change in the absorption spectrum of the aromatase cytochrome P-450 which is indicative of substrate displacement. The relative binding affinities of the flavonoid analogs were determined and only ring A adn ring B dihydroxylated analogs were found to bind to the estrogen receptor.     

Structural characterization of the complex between alpha-naphthoflavone and human cytochrome P450 1B1

The atomic structure of human P450 1B1 was determined by x-ray crystallography to 2.7 Å resolution with α-naphthoflavone (ANF) bound in the active site cavity. Although the amino acid sequences of human P450s 1B1 and 1A2 have diverged significantly, both enzymes exhibit narrow active site cavities, which underlie similarities in their substrate profiles. Helix I residues adopt a relatively flat conformation in both enzymes, and a characteristic distortion of helix F places Phe²³¹ in 1B1 and Phe²²⁶ in 1A2 in similar positions for π-π stacking with ANF. ANF binds in a distinctly different orientation in P450 1B1 from that observed for 1A2. This reflects, in part, divergent conformations of the helix B′-C loop that are stabilized by different hydrogen-bonding interactions in the two enzymes. Additionally, differences between the two enzymes for other amino acids that line the edges of the cavity contribute to distinct orientations of ANF in the two active sites. Thus, the narrow cavity is conserved in both P450 subfamily 1A and P450 subfamily 1B with sequence divergence around the edges of the cavity that modify substrate and inhibitor binding. The conservation of these P450 1B1 active site amino acid residues across vertebrate species suggests that these structural features are conserved.     

Differential mechanisms for the inhibition of human cytochrome P450 1A2 by apigenin and genistein

The inhibitory effects of flavonoids on the human cytochrome P450 1A2 (CYP1A2) were examined. Among flavonoids tested, galangin, kaempferol, chrysin, and apigenin were potent inhibitors. Although apigenin belonging to flavones and genistein belonging to isoflavones are similar in the chemical structures, the inhibitory potencies for CYP1A2 were distinguished markedly between these two flavonoids. In computer‐docking simulation, apigenin interacted with the same mode of cocrystallized α‐naphthoflavone in the active site of CYP1A2, and then the B ring of apigenin was placed close to the heme iron of the enzyme with a single orientation. In contrast, the docked genistein conformation showed two different binding modes, and the A ring of genistein was oriented to the heme iron of CYP1A2. Furthermore, the binding free energy of apigenin was lower than that of genistein. These results demonstrate a possible mechanism that causes the differential inhibitory potencies of apigenin and genistein for CYP1A2. © 2010 Wiley Periodicals, Inc. J Biochem Mol Toxicol 24:230–234, 2010; View this article online at wileyonlinelibrary.com . DOI 10.1002/jbt.20328     

Cytochrome P450s in flavonoid metabolism

In this review, cytochrome P450s characterized at the molecular level catalyzing aromatic hydroxylations, aliphatic hydroxylations and skeleton formation in the flavonoid metabolism are surveyed. They are involved in the biosynthesis of anthocyanin pigments and condensed tannin (CYP75, flavonoid 3′,5′-hydroxylase and 3′-hydroxylase), flavones [CYP93B, (2S)-flavanone 2-hydroxylase and flavone synthase II], and leguminous isoflavonoid phytoalexins [CYP71D9, flavonoid 6-hydroxylase; CYP81E, isoflavone 2′-hydroxylase and 3′-hydroxylase; CYP93A, 3,9-dihydroxypterocarpan 6a-hydroxylase; CYP93C, 2-hydroxyisoflavanone synthase (IFS)]. Other P450s of the flavonoid metabolism include methylenedioxy bridge forming enzyme, cyclases producing glyceollins, flavonol 6-hydroxylase and 8-dimethylallylnaringenin 2′-hydroxylase. Mechanistic studies on the unusual aryl migration by CYP93C, regulation of IFS expression in plant organs and its biotechnological applications are introduced, and flavonoid metabolisms by non-plant P450s are also briefly discussed.     

Substrate-modulated Cytochrome P450 17A1 and Cytochrome b5 Interactions Revealed by NMR

The membrane heme protein cytochrome b₅ (b₅) can enhance, inhibit, or have no effect on cytochrome P450 (P450) catalysis, depending on the specific P450, substrate, and reaction conditions, but the structural basis remains unclear. Here the interactions between the soluble domain of microsomal b₅ and the catalytic domain of the bifunctional steroidogenic cytochrome P450 17A1 (CYP17A1) were investigated. CYP17A1 performs both steroid hydroxylation, which is unaffected by b₅, and an androgen-forming lyase reaction that is facilitated 10-fold by b₅. NMR chemical shift mapping of b₅ titrations with CYP17A1 indicates that the interaction occurs in an intermediate exchange regime and identifies charged surface residues involved in the protein/protein interface. The role of these residues is confirmed by disruption of the complex upon mutagenesis of either the anionic b₅ residues (Glu-48 or Glu-49) or the corresponding cationic CYP17A1 residues (Arg-347, Arg-358, or Arg-449). Cytochrome b₅ binding to CYP17A1 is also mutually exclusive with binding of NADPH-cytochrome P450 reductase. To probe the differential effects of b₅ on the two CYP17A1-mediated reactions and, thus, communication between the superficial b₅ binding site and the buried CYP17A1 active site, CYP17A1/b₅ complex formation was characterized with either hydroxylase or lyase substrates bound to CYP17A1. Significantly, the CYP17A1/b₅ interaction is stronger when the hydroxylase substrate pregnenolone is present in the CYP17A1 active site than when the lyase substrate 17α-hydroxypregnenolone is in the active site. These findings form the basis for a clearer understanding of this important interaction by directly measuring the reversible binding of the two proteins, providing evidence of communication between the CYP17A1 active site and the superficial proximal b₅ binding site.     

Isoform-selective inhibition of chrysin towards human cytochrome P450 1A2. Kinetics analysis,molecular docking,and molecular dynamics simulations

Our kinetics studies demonstrated that the nature product chrysin exhibited a high inhibitory affinity of 54 nM towards human cytochrome P450 1A2 and was comparable to α-naphthoflavone (49 nM), whereas it represented a moderate affinity of 5225 nM against human cytochrome P450 2C9. However, it remains unclear how this inhibitor selectively binds 1A2. To better understand the isoform selectivity of chrysin, molecular docking and molecular dynamics simulations were performed. Chrysin formed a strong H-bond with Asp313 of 1A2. The stacking interactions with Phe226 also contributed to its tight binding to 1A2. The larger and much more open active site architectures of 2C9 may explain the weaker inhibitory affinity of chrysin towards 2C9. The predicted binding free energies suggest that chrysin preferred 1A2 (ΔG_{bind, pred} = ?23.11 kcal/mol) to 2C9 (?20.41 kcal/mol). Additionally, the present work revealed that 7-hydroxy-flavone bound to 1A2 in a similar pattern as chrysin and represented a slightly less negative predicted binding free energy, which was further validated by our kinetics analysis (IC₅₀ = 240 nM). Results of the study can provide insight for designing novel isoform-selective 1A2 inhibitors.     

Resveratrol is a selective human cytochrome P450 1A1 inhibitor.

Resveratrol (trans-3,4',5-trihydroxystilbene) is a phytoalexin compound found in juice and wine produced from dark-skinned grape cultivars and reported to have anti-inflammatory and anticarcinogenic activities. To investigate the mechanism of anticarcinogenic activities of resveratrol, the effects on cytochrome P450 (P450) were determined in human liver microsomes and Escherichia coli membranes coexpressing human P450 1A1 or P450 1A2 with human NADPH-P450 reductase (bicistronic expression system). Resveratrol slightly inhibited ethoxyresorufin O-deethylation (EROD) activity in human liver microsomes with an IC(50) of 1.1 mM. Interestingly, resveratrol exhibited potent inhibition of human P450 1A1 in a dose-dependent manner with IC(50) of 23 microM for EROD and IC(50) of 11 microM for methoxyresorufin O-demethylation (MROD). However, the inhibition of human P450 1A2 by resveratrol was not so strong (IC(50) 1.2 mM for EROD and 580 microM for MROD). Resveratrol showed over 50-fold selectivity for P450 1A1 over P450 1A2. The activities of human NADPH-P450 reductase were not significantly changed by resveratrol. In a human P450 1A1/reductase bicistronic expression system, resveratrol inhibited human P450 1A1 activity in a mixed-type inhibition (competitive-noncompetitive) with a K(i) values of 9 and 89 microM. These results suggest that resveratrol is a selective human P450 1A1 inhibitor, and may be considered for use as a strong cancer chemopreventive agent in humans.     

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.     

Homology model of human retinoic acid metabolising enzyme cytochrome P450 26A1 (CYP26A1): active site architecture and ligand binding

Homology models of cytochrome P450 RA1 (CYP26A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using MOE software the lowest energy CYP26A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP26A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP26A1 model. The natural substrate, all-trans-retinoic acid (atRA), and inhibitor R 15866, were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure. Multiple hydrophobic interactions, including the side chains of TRP112, PHE299, PHE222, PHE84, PHE374 and PRO371, are important for binding of atRA and R115866. Additional hydrogen bonding interactions were noted as follows: atRA-- C==O of the atRA carboxylate group and ARG86; R115866--benzothiazole nitrogen and the backbone NH of SER115.     

Exploring CYP1A1 as anticancer target: homology modeling and in silico inhibitor design

Microsomal cytochrome P450 family 1 enzymes has great importance in the bioactivation of mutagens. P450 catalyzed reactions involve a wide range of substrates, and this versatality is reflected in a structural diversity, evident in the active sites of available P450 structures. This structure offers a template to study further structure-function relationships of alternative substrates and other cytochrome P450 family 1 members. In this paper, we document a homology model of CYP P450 1A1 from Homo sapiens, developed on the basis of template crystal structure of human microsomal P450 1a2 in complex with inhibitor (PDB Id: 2HI4). Homology modeling is performed at the programs, both in the commercial and public realms. We tried to explore CYP1A1 as a potential target for anticancer chemotherapy. To gain an insight into the binding of ligands with enzyme, protein-ligand complex was developed by including information about the known ligand as spatial restraints and optimizing the mutual interactions between the ligand and the binding site. Active site characterization and the study for involvement of specific aminoacids in binding with ligand, facilitates structure based inhibitor design. This study should prove useful in the design and development of potential novel anticancer agents. Figure The refined structure of homology model of CYP1A1     

Targeting of the highly conserved threonine 302 residue of cytochromes P450 2B family during mechanism-based inactivation by aryl acetylenes

Cytochromes P450 (CYPs or P450s) contain a highly conserved threonine residue in the active site, which is referred to as Thr302 in the amino acid sequence of CYP2B4. Extensive biochemical and crystallographic studies have established that this Thr302 plays a critical role in activating molecular oxygen to generate Compound I, a putative iron(IV)-oxo porphyrin cation radical, that carries out the preliminary oxygenation of CYP substrates. Because of its proximity to the center of the P450 active site, this Thr302 is susceptible to mechanism-based inactivation under certain conditions. In this article, we review recent studies on the mechanism-based inactivation of three mammalian P450s in the 2B family, CYP2B1 (rat), 2B4 (rabbit) and 2B6 (human) by tert-butylphenylacetylene (tBPA). These studies showed that tBPA is a potent mechanism-based inactivator of CYP2B1, 2B4 and 2B6 with high k_inact/K_I ratios (0.23–2.3 min⁻¹ μM⁻¹) and low partition ratios (0–5). Furthermore, mechanistic studies revealed that tBPA inactivates these three CYP2B enzymes through the formation of a single ester adduct with the Thr302 in the active site. These inhibitory properties of tBPA allowed the preparation of a modified CYP2B4 where the Thr302 was covalently and stoichiometrically labeled by a reactive intermediate of tBPA in quantities large enough to permit spectroscopic and crystallographic studies of the consequences of covalent modification of Thr302. Molecular modeling studies revealed a unique binding mode of tBPA in the active site that may shed light on the potency of this inhibition. The results from these studies may serve as a basis for designing more specific and potent inhibitors for P450s by targeting this highly conserved threonine residue which is present in the active sites of most mammalian P450s.     

A survey of active site access channels in cytochromes P450

In cytochrome P450s, the active site is situated deep inside the protein next to the heme cofactor, and is often completely isolated from the surrounding solvent. To identify routes by which substrates may enter into and products exit from the active site, random expulsion molecular dynamics simulations were performed for three cytochrome P450s: CYP101, CYP102A1 and CYP107A1 [J. Mol. Biol. 303 (2000) 797; Proc. Natl. Acad. Sci. USA 99 (2002) 5361]. Amongst the different pathways identified, one pathway was found to be common to all three cytochrome P450s although the mechanism of ligand passage along it was different in each case and apparently adapted to the substrate specificity of the enzyme. Recently, a number of new crystal structures of cytochrome P450s have been solved. Here, we analyse the open channels leading to the active site that these structures reveal. We find that in addition to showing the common pathway, they provide experimental evidence for the existence of three additional channels that were identified by simulation. We also discuss how the location of xenon binding sites in CYP101 suggests a role for one of the pathways identified by molecular dynamics simulations as a route for gaseous species, such as oxygen, to access the active site.     

A proteolytically sensitive region common to several rat liver cytochromes P450: effect of cleavage on substrate binding.

Limited proteolysis of rat liver microsomes was used to probe the topography and structure of cytochrome P450 bound to the endoplasmic reticulum. Three cytochromes P450 from two families were examined. Monoclonal antibodies to cytochrome P450 forms 1A1, 2B1, and 2E1 were used to immunopurify these proteolyzed cytochromes P450 from microsomes from rats treated with 3-methylcholanthrene, phenobarbital, and acetone, respectively. Electrophoretic and immunoblot analysis of tryptic fragments revealed a highly sensitive cleavage site in all three cytochromes P450. N-Terminal sequencing was performed on the fragments after transfer onto poly(vinylidene difluoride) membranes and showed that this preferential cleavage site is at amino acid position 298 of P450 1A1, position 277 of P450 2B1, and position 278 of P450 2E1. Multiple sequence alignment revealed that these positions are at the amino terminal of a highly conserved region of these cytochromes P450. The important functional role implied by primary sequence conservation along with the proteolytic sensitivity at its amino terminal suggests that this region is a protein domain. Comparison with the known structure of the bacterial cytochrome P450cam predicts that this proteolytically sensitive site is within an interhelical turn region connected to the distal helix that partially encompasses the heme-containing active site. Substrate binding to the cleaved cytochromes P450 was examined in order to determine whether the newly added conformational freedom near the cleavage site functionally altered these cytochromes P450. Cleavage of P450 2B1 abolished benzphetamine binding, which indicates that the cleavage site contains an important structural determinant for binding this substrate. However, cleavage did not affect benzo[a]pyrene binding to P450 1A1.     

Binding characteristics of aromatase inhibitors and phytoestrogens to human aromatase

We have evaluated the binding characteristics of three steroidal inhibitors [4-hydroxyandrostene-dione (4-OHA), 7-(4′-amino)phenylthio-1,4-androstadiene-3,17-dione (7-APTADD), and bridge (2,19-methyleneoxy) androstene-3,17-dione (MDL 101,003)], four nonsteroidal inhibitors [aminoglutethimide (AG), CGS 20267, ICI D1033, and vorozole (R83842)], and two flavone phytoestrogens (chrysin, and 7,8-dihydroxyflavone) to aromatase through a combination of computer modeling and inhibitory profile studies on the wild-type and six aromatase mutants (I133Y, P308F, D309A, T310S, I395F, and I474Y). We have generated two aromatase models based on the x-ray structures of cytochrome P450-cam and cytochrome P450bm3, respectively. A major difference between the cytochrome P450cam-based and cytochrome P450bm3-based models is in the predicted lengths of helices F and G. In the cytochrome P450cam-based model, helices F and G lie antiparallel and extend across the active-site face of the molecule from one edge to the center, so that the carboxyl-terminal residues of helix F and the N-terminal residues of helix G make a major contribution to the structure of the active site. In the cytochrome P450bm3-based model, both helices are longer and so extend almost all the way across the active-site face of the molecule. Considering the size of the androgen substrate, we evaluated our results mainly based on the cytochrome P450cam model. The mutations involved in this study are thought to be at or near the proposed active site pocket. The inhibitory profile analysis has produced very interesting results and provided a molecular basis as to how seven aromatase inhibitors with different structures bind to the active site of aromatase. Furthermore, the investigation reveals that phytoestrogens bind to the active site of aromatase in a different orientation from that in the estrogen receptor.