Different binding modes of two flaviolin substrate molecules in cytochrome P450 158A1 (CYP158A1) compared to CYP158A2

Cytochrome P450 158A2 (CYP158A2) has been shown to catalyze an unusual oxidative C-C coupling reaction to polymerize flaviolin and form highly conjugated pigments (three isomers of biflaviolin and one triflaviolin) in Streptomyces coelicolor A3(2) which protect the soil bacterium from deleterious effects of UV irradiation (Zhao B. et al. (2005) J. Biol. Chem. 280, 11599-11607). The present studies demonstrate that the subfamily partner CYP158A1, sharing 61% amino acid identity with CYP158A2, can also catalyze the same flaviolin dimerization reactions, but it generates just two of the three isomers of biflaviolin that CYP158A2 produces. Furthermore, the two CYP158A1 products have very different molar ratios compared with the corresponding CYP158A2 products, indicating that each enzyme maintains its own stereo- and regiospecificity. To find an explanation for these differences, three CYP158A1 structures have been solved by X-ray crystallography and have been compared with those for CYP158A2. The structures reveal surprising differences. Particularly, only one flaviolin molecule is present close to the heme iron in CYP158A1, and the second flaviolin molecule binds at the entrance of the putative substrate access channel on the protein distal surface 9 A away. Our work describes two members of the same P450 subfamily, which produce the same products by oxidative C-C coupling yet show very different structural orientations of substrate molecules in the active site.     

Role of active site water molecules and substrate hydroxyl groups in oxygen activation by cytochrome P450 158A2: a new mechanism of proton transfer

From the x-ray crystal structure of CYP158A2 (Zhao, B., Guengerich, F. P., Bellamine, A., Lamb, D. C., Izumikawa, M., Lei, L., Podust, L. M., Sundaramoorthy, M., Reddy, L. M., Kelly, S. L., Kalaitzis, J. A., Stec, D., Voehler, M., Falck, J. R., Moore, B. S., Shimada, T., and Waterman, M. R. (2005) J. Biol. Chem. 280, 11599-11607), one of 18 cytochrome P450 (CYP) genes in the actinomycete Streptomyces coelicolor, ordered active site water molecules (WAT505, WAT600, and WAT640), and hydroxyl groups of the substrate flaviolin were proposed to participate in proton transfer and oxygen cleavage in this monooxygenase. To probe their roles in catalysis, we have studied the crystal structures of a substrate analogue (2-hydroxy-1,4-naphthoquinone) complex with ferric CYP158A2 (2.15 A) and the flaviolin ferrous dioxygen-bound CYP158A2 complex (1.8 A). Catalytic activity toward 2-hydroxy-1,4-naphthoquinone was approximately 70-fold less than with flaviolin. In the ferrous dioxygen-bound flaviolin complex, the three water molecules in the ferric flaviolin complex still occupy the same positions and form hydrogen bonds to the distal dioxygen atom. These findings suggest that CYP158A2 utilizes substrate hydroxyl groups to stabilize active site water and further assist in the iron-linked dioxygen activation. A continuous hydrogen-bonded water network connecting the active site to the protein surface (bulk solvent) not present in the other two ferrous dioxygen-bound P450 structures (CYP101A1/P450cam and CYP107A1/P450eryF) is proposed to participate in the proton-delivery cascade, leading to dioxygen bond scission. This ferrous-dioxygen structure suggests two classes of P450s based on the pathway of proton transfer, one using the highly conserved threonine in the I-helix (CYP101A1) and the other requiring hydroxyl groups of the substrate molecules either directly transferring protons (CYP107A1) or stabilizing a water pathway for proton transfer (CYP158A2).     

The role of Ile87 of CYP158A2 in oxidative coupling reaction

Both CYP158A1 and CYP158A2 are able to catalyze an oxidative C-C coupling reaction producing biflaviolin or triflaviolin in Streptomyces coelicolor A3(2). The substrate-bound crystal structures of CYP158A2 and CYP158A1 reveal that the side chain of Ile87 in CYP158A2 points to the active site contacting the distal flaviolin molecule, however, the bulkier side chain of Lys90 in CYP158A1 (corresponding to Ile87 in CYP158A2) is toward the distal surface of the protein. These results suggest that these residues could be important in determining product regiospecificity. In order to explore the role of the two residues in catalysis, the reciprocal mutants, Ile87Lys and Lys90Ile, of CYP158A2 and CYP158A1, respectively, were generated and characterized. The mutant Ile87Lys enzyme forms two isomers of biflaviolin instead of three isomers of biflaviolin in wild-type CYP158A2. CYP158A1 containing the substitution of lysine with isoleucine has the same catalytic activity compared with the wild-type CYP158A1. The crystal structure of Ile87Lys showed that the BC loop in the mutant is in a very different orientation compared with the BC loop in both CYP158A1/A2 structures. These results shed light on the mechanism of the oxidative coupling reaction catalyzed by cytochrome P450.     

Ligand-assisted inhibition in cytochrome P450 158A2 from Streptomyces coelicolor A3(2)

Cytochrome P450 158A2 (CYP158A2) can polymerize flaviolin to red-brown pigments, which may afford physical protection to the organism, possibly against the deleterious effects of UV radiation. We have found that the small molecule malonic acid enables cocrystallization of this mixed function oxidase with the azole inhibitor 4-phenylimidazole. The presence of malonate molecules affects the behavior of the binding of 4-phenylimidazole to CYP158A2 and increases inhibition potency up to 2-fold compared to 4-phenylimidazole alone. We report here the crystal structure of the 4-phenylimidazole/malonate complex of CYP158A2 at 1.5 A. Two molecules of malonate used in crystallization are found above the single inhibitor molecule in the active site. Those two molecules are linked between the BC loop and beta 1-4/beta 6-1 strands via hydrogen bond interactions to stabilize the conformational changes of the BC loop and beta strands that take place upon inhibitor binding compared to the ligand-free structure we have reported previously. 4-Phenylimidazole can launch an extensive hydrogen-bonding network in the region of the F/G helices which may stabilize the conformational changes. Our findings clearly show that two molecules of malonate assist the inhibitor 4-phenylimidazole to assume a specific location producing more inhibition in the enzyme catalytic activity.     

Investigating conservation of the albaflavenone biosynthetic pathway and CYP170 bifunctionality in streptomycetes

Albaflavenone, a tricyclic sesquiterpene antibiotic, is biosynthesized in Streptomyces coelicolor A3(2) by enzymes encoded in a two-gene operon. Initially, sesquiterpene cyclase catalyzes the cyclization of farnesyl diphosphate to the terpenoid epi-isozizaene, which is oxidized to the final albaflavenone by cytochrome P450 (CYP)170A1. Additionally, this CYP is a bifunctional enzyme, being able to also generate farnesene isomers from farnesyl diphosphate, owing to a terpene synthase active site moonlighting on the CYP molecule. To explore the functionality of this operon in other streptomycetes, we have examined culture extracts by GC/MS and established the presence of albaflavenone in five Streptomyces species. Bioinformatics examination of the predicted CYP170 primary amino acid sequences revealed substitutions in the CYP terpene synthase active site. To examine whether the terpene synthase site was catalytically active in another CYP170, we characterized the least related CYP170 orthologue from Streptomyces albus (CYP170B1). Following expression and purification, CYP170B1 showed a normal reduced CO difference spectrum at 450 nm, in contrast to the unusual 440-nm peak observed for S. coelicolor A3(2) CYP170A1. CYP170B1 can catalyze the conversion of epi-isozizaene to albaflavenone, but was unable to catalyze the conversion of farnesyl diphosphate to farnesene. Molecular modeling with our crystal structure of CYP170A1 suggests that the absence of key amino acids for binding the essential terpene synthase cofactor Mg(2+) may be the explanation for the loss of CYP170B1 bifunctionality.     

Crystal Structure of Albaflavenone Monooxygenase Containing a Moonlighting Terpene Synthase Active Site

Albaflavenone synthase (CYP170A1) is a monooxygenase catalyzing the final two steps in the biosynthesis of this antibiotic in the soil bacterium, Streptomyces coelicolor A3(2). Interestingly, CYP170A1 shows no stereo selection forming equal amounts of two albaflavenol epimers, each of which is oxidized in turn to albaflavenone. To explore the structural basis of the reaction mechanism, we have studied the crystal structures of both ligand-free CYP170A1 (2.6 Å) and complex of endogenous substrate (epi-isozizaene) with CYP170A1 (3.3 Å). The structure of the complex suggests that the proximal epi-isozizaene molecules may bind to the heme iron in two orientations. In addition, much to our surprise, we have found that albaflavenone synthase also has a second, completely distinct catalytic activity corresponding to the synthesis of farnesene isomers from farnesyl diphosphate. Within the cytochrome P450 α-helical domain both the primary sequence and x-ray structure indicate the presence of a novel terpene synthase active site that is moonlighting on the P450 structure. This includes signature sequences for divalent cation binding and an α-helical barrel. This barrel is unusual because it consists of only four helices rather than six found in all other terpene synthases. Mutagenesis establishes that this barrel is essential for the terpene synthase activity of CYP170A1 but not for the monooxygenase activity. This is the first bifunctional P450 discovered to have another active site moonlighting on it and the first time a terpene synthase active site is found moonlighting on another protein.     

Structures of human cytochrome P-450 2E1. Insights into the binding of inhibitors and both small molecular weight and fatty acid substrates

Human microsomal cytochrome P-450 2E1 (CYP2E1) monooxygenates > 70 low molecular weight xenobiotic compounds, as well as much larger endogenous fatty acid signaling molecules such as arachidonic acid. In the process, CYP2E1 can generate toxic or carcinogenic compounds, as occurs with acetaminophen overdose, nitrosamines in cigarette smoke, and reactive oxygen species from uncoupled catalysis. Thus, the diverse roles that CYP2E1 has in normal physiology, toxicity, and drug metabolism are related to its ability to metabolize diverse classes of ligands, but the structural basis for this was previously unknown. Structures of human CYP2E1 have been solved to 2.2 angstroms for an indazole complex and 2.6 angstroms for a 4-methylpyrazole complex. Both inhibitors bind to the heme iron and hydrogen bond to Thr303 within the active site. Complementing its small molecular weight substrates, the hydrophobic CYP2E1 active site is the smallest yet observed for a human cytochrome P-450. The CYP2E1 active site also has two adjacent voids: one enclosed above the I helix and the other forming a channel to the protein surface. Minor repositioning of the Phe478 aromatic ring that separates the active site and access channel would allow the carboxylate of fatty acid substrates to interact with conserved 216QXXNN220 residues in the access channel while positioning the hydrocarbon terminus in the active site, consistent with experimentally observed omega-1 hydroxylation of saturated fatty acids. Thus, these structures provide insights into the ability of CYP2E1 to effectively bind and metabolize both small molecule substrates and fatty acids.     

Structural comparison of cytochromes P450 2A6, 2A13, and 2E1 with pilocarpine

Human xenobiotic-metabolizing cytochrome P450 (CYP) enzymes can each bind and monooxygenate a diverse set of substrates, including drugs, often producing a variety of metabolites. Additionally, a single ligand can interact with multiple CYP enzymes, but often the protein structural similarities and differences that mediate such overlapping selectivity are not well understood. Even though the CYP superfamily has a highly canonical global protein fold, there are large variations in the active site size, topology, and conformational flexibility. We have determined how a related set of three human CYP enzymes bind and interact with a common inhibitor, the muscarinic receptor agonist drug pilocarpine. Pilocarpine binds and inhibits the hepatic CYP2A6 and respiratory CYP2A13 enzymes much more efficiently than the hepatic CYP2E1 enzyme. To elucidate key residues involved in pilocarpine binding, crystal structures of CYP2A6 (2.4 ?), CYP2A13 (3.0 ?), CYP2E1 (2.35 ?), and the CYP2A6 mutant enzyme, CYP2A6?I208S/I300F/G301A/S369G (2.1 ?) have been determined with pilocarpine in the active site. In all four structures, pilocarpine coordinates to the heme iron, but comparisons reveal how individual residues lining the active sites of these three distinct human enzymes interact differently with the inhibitor pilocarpine.     

Nicotine and 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone Binding and Access Channel in Human Cytochrome P450 2A6 and 2A13 Enzymes

Cytochromes P450 (CYP) from the 2A subfamily are known for their roles in the metabolism of nicotine, the addictive agent in tobacco, and activation of the tobacco procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK). Although both the hepatic CYP2A6 and respiratory CYP2A13 enzymes metabolize these compounds, CYP2A13 does so with much higher catalytic efficiency, but the structural basis for this has been unclear. X-ray structures of nicotine complexes with CYP2A13 (2.5 Å) and CYP2A6 (2.3 Å) yield a structural rationale for the preferential binding of nicotine to CYP2A13. Additional structures of CYP2A13 with NNK reveal either a single NNK molecule in the active site with orientations corresponding to metabolites known to form DNA adducts and initiate lung cancer (2.35 Å) or with two molecules of NNK bound (2.1 Å): one in the active site and one in a more distal staging site. Finally, in contrast to prior CYP2A structures with enclosed active sites, CYP2A13 conformations were solved that adopt both open and intermediate conformations resulting from an ∼2.5 Å movement of the F to G helices. This channel occurs in the same region where the second, distal NNK molecule is bound, suggesting that the channel may be used for ligand entry and/or exit from the active site. Altogether these structures provide multiple new snapshots of CYP2A13 conformations that assist in understanding the binding and activation of an important human carcinogen, as well as critical comparisons in the binding of nicotine, one of the most widely used and highly addictive drugs in human use.     

Structure of the human lung cytochrome P450 2A13

The human lung cytochrome P450 2A13 (CYP2A13) activates the nicotine-derived procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) into DNA-altering compounds that cause lung cancer. Another cytochrome P450, CYP2A6, is also present in human lung, but at much lower levels. Although these two enzymes are 93.5% identical, CYP2A13 metabolizes NNK with much lower K(m) values than does CYP2A6. To investigate the structural differences between these two enzymes the structure of CYP2A13 was determined to 2.35A by x-ray crystallography and compared with structures of CYP2A6. As expected, the overall CYP2A13 and CYP2A6 structures are very similar with an average root mean square deviation of 0.5A for the Calpha atoms. Like CYP2A6, the CYP2A13 active site cavity is small and highly hydrophobic with a cluster of Phe residues composing the active site roof. Active site residue Asn(297) is positioned to hydrogen bond with an adventitious ligand, identified as indole. Amino acid differences between CYP2A6 and CYP2A13 at positions 117, 300, 301, and 208 relate to different orientations of the ligand plane in the two protein structures and may underlie the significant variations observed in binding and catalysis of many CYP2A ligands. In addition, docking studies suggest that residues 365 and 366 may also contribute to differences in NNK metabolism.     

Modelling of three-dimensional structures of cytochromes P450 11B1 and 11B2.

The final steps of the biosynthesis of glucocorticoids and mineralocorticoids in the adrenal cortex require the action of two different cytochromes P450--CYP11B1 and CYP11B2. Homology modelling of the three-dimensional structures of these cytochromes was performed based on crystallographic coordinates of two bacterial P450s, CYP102 (P450BM-3) and CYP108 (P450terp). Principal attention was given to the modelling of the active sites and a comparison of the active site structures of CYP11B1 and CYP11B2 was performed. It can be demonstrated that key residue contacts within the active site appear to depend on the orientation of the heme. The obtained 3D structures of CYP11B1 and CYP11B2 were used for investigation of structure-function relationships of these enzymes. Previously obtained results on naturally occurring mutants and on mutants obtained by site-directed mutagenesis are discussed.     

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.     

Structures of human microsomal cytochrome P450 2A6 complexed with coumarin and methoxsalen

Human microsomal cytochrome P450 2A6 (CYP2A6) contributes extensively to nicotine detoxication but also activates tobacco-specific procarcinogens to mutagenic products. The CYP2A6 structure shows a compact, hydrophobic active site with one hydrogen bond donor, Asn297, that orients coumarin for regioselective oxidation. The inhibitor methoxsalen effectively fills the active site cavity without substantially perturbing the structure. The structure should aid the design of inhibitors to reduce smoking and tobacco-related cancers.     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) 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.     

Differences in flexibility of active sites of cytochromes P450 probed by resonance Raman and UV-Vis absorption spectroscopy.

Spectroscopic methods reveal differences in flexibility and stability of P450 forms. Among microsomal P450s, the most flexible active site has been found in the CYP3A4 enzyme as it is compressible and the heme vinyl side chains may adopt two different conformations. On the other hand, active site of this enzyme denatures quite easily upon hydrostatic pressure. The most rigid active site able to withstand the effect of high pressure has CYP1A2. The bacterial CYP102 (BM3) flavocytochrome has also a rather stable, but flexible active site. The differences between CYP3A4 and CYP1A2 active sites apparently reflect their ability to bind various substrates: whereas the CYP3A4 binds a vast variety of structures, the CYP1A2 preferentially binds planar, aromatic structures and its substrate specificity is relatively narrow.     

Computational solvent mapping reveals the importance of local conformational changes for broad substrate specificity in mammalian cytochromes P450

Computational solvent mapping moves small organic molecules as probes around a protein surface, finds favorable binding positions, clusters the conformations, and ranks the clusters on the basis of their average free energy. Prior mapping studies of enzymes, crystallized in either substrate-free or substrate-bound form, have shown that the largest number of solvent probe clusters invariably overlaps in the active site. We have applied this method to five cytochromes P450. As expected, the mapping of two bacterial P450s, P450 cam (CYP101) and P450 BM-3 (CYP102), identified the substrate-binding sites in both ligand-bound and ligand-free P450 structures. However, the mapping finds the active site only in the ligand-bound structures of the three mammalian P450s, 2C5, 2C9, and 2B4. Thus, despite the large cavities seen in the unbound structures of these enzymes, the features required for binding small molecules are formed only in the process of substrate binding. The ability of adjusting their binding sites to substrates that differ in size, shape, and polarity is likely to be responsible for the broad substrate specificity of these mammalian P450s. Similar behavior was seen at "hot spots" of protein-protein interfaces that can also bind small molecules in grooves created by induced fit. In addition, the binding of S-warfarin to P450 2C9 creates a high-affinity site for a second ligand, which may help to explain the prevalence of drug-drug interactions involving this and other mammalian P450s.     

Human Cytochrome P450 2E1 Structures with Fatty Acid Analogs Reveal a Previously Unobserved Binding Mode

Human microsomal cytochrome P450 (CYP) 2E1 is widely known for its ability to oxidize >70 different, mostly compact, low molecular weight drugs and other xenobiotic compounds. In addition CYP2E1 oxidizes much larger C9–C20 fatty acids that can serve as endogenous signaling molecules. Previously structures of CYP2E1 with small molecules revealed a small, compact CYP2E1 active site, which would be insufficient to accommodate medium and long chain fatty acids without conformational changes in the protein. In the current work we have determined how CYP2E1 can accommodate a series of fatty acid analogs by cocrystallizing CYP2E1 with ω-imidazolyl-octanoic fatty acid, ω-imidazolyl-decanoic fatty acid, and ω-imidazolyl-dodecanoic fatty acid. In each structure direct coordination of the imidazole nitrogen to the heme iron mimics the position required for native fatty acid substrates to yield the ω-1 hydroxylated metabolites that predominate experimentally. In each case rotation of a single Phe²⁹⁸ side chain merges the active site with an adjacent void, significantly altering the active site size and topology to accommodate fatty acids. The binding of these fatty acid ligands is directly opposite the channel to the protein surface and the binding observed for fatty acids in the bacterial cytochrome P450 BM3 (CYP102A1) from Bacillus megaterium. Instead of the BM3-like binding mode in the CYP2E1 channel, these structures reveal interactions between the fatty acid carboxylates and several residues in the F, G, and B′ helices at successive distances from the active site.     

Analysis of Species-Selectivity of Human,Mouse and Rat Cytochrome P450 1A and 2B Subfamily Enzymes using Molecular Modeling,Docking and Dynamics Simulations

Cytochrome P450 (CYP) 1A and 2B subfamily enzymes are important drug metabolizing enzymes, and are highly conserved across species in terms of sequence homology. However, there are major to minor structural and macromolecular differences which provide for species-selectivity and substrate-selectivity. Therefore, species-selectivity of CYP1A and CYP2B subfamily proteins across human, mouse and rat was analyzed using molecular modeling, docking and dynamics simulations when the chiral molecules quinine and quinidine were used as ligands. The three-dimensional structures of 17 proteins belonging to CYP1A and CYP2B subfamilies of mouse and rat were predicted by adopting homology modeling using the available structures of human CYP1A and CYP2B proteins as templates. Molecular docking and dynamics simulations of quinine and quinidine with CYP1A subfamily proteins revealed the existence of species-selectivity across the three species. On the other hand, in the case of CYP2B subfamily proteins, no role for chirality of quinine and quinidine in forming complexes with CYP2B subfamily proteins of the three species was indicated. Our findings reveal the roles of active site amino acid residues of CYP1A and CYP2B subfamily proteins and provide insights into species-selectivity of these enzymes across human, mouse, and rat.