Crystal structure and refinement of cytochrome P450terp at 2.3 A resolution

Spatial relationship between L-arginine and heme binding sites of endothelial nitric-oxide synthase

Binding of L-arginine and imidazole to the endothelial nitric-oxide synthase (eNOS) was characterized by direct heme spectral perturbation. L-Arginine is competitive with imidazole for binding to eNOS. Both equilibrium binding and kinetic binding were measured at 4 and 23 degrees C for these two ligands. Kd (imidazole) is 60 microM and 110 microM, kon (imidazole) is 2.5 x 10(5) M-1 s-1 and 1. 2 x 10(6) M-1 s-1, koff (imidazole) is 11.8 s-1 and 116 s-1 at 4 and 23 degrees C, respectively. Corresponding values for L-arginine are calculated from the data of binding competition with imidazole and computer modeling. Kd (L-arginine) is 0.5 microM and 2.0 microM, kon (L-arginine) is 2 x 10(5) M-1 s-1 and 8 x 10(5) M-1 s-1, koff (L-arginine) is 0.08 s-1 and 1.6 s-1 at 4 and 23 degrees C, respectively. It is suggested that binding of both ligands occurs through the same access channel to the heme site based on their similarly slow association rate constants. A series of potential heme ligands and amino acid analogs of L-arginine were evaluated for their binding and their effect on the heme structure. All ligands besides cyanide tested for binding inhibition are competitive with either L-arginine or imidazole. The space for the distal heme ligand was estimated to be approximately 6.3 x 6.7 A by three groups of rigid planar ligands: imidazole, pyridine, and pyrimidine. Results of the thiazole and amino acid ligand series permitted the conclusion that the guanidine group of L-arginine is critical for its binding affinity and its specific orientation relative to the heme. Such a specific conformation is essential for the oxygenase mechanism of eNOS.

A three-dimensional protein model for human cytochrome P450 2D6 based on the crystal structures of P450 101, P450 102, and P450 108

Cytochromes P450 (P450s) constitute a superfamily of phase I enzymes capable of oxidizing and reducing various substrates. P450 2D6 is a polymorphic enzyme, which is absent in 5-9% of the Caucasian population as a result of a recessive inheritance of gene mutations. This deficiency leads to impaired metabolism of a variety of drugs. All drugs metabolized by P450 2D6 contain a basic nitrogen atom, and a flat hydrophobic region coplanar to the oxidation site which is either 5 or 7 A away from the basic nitrogen atom. The aim of this study was to build a three-dimensional structure for the protein and more specifically for the active site of P450 2D6 in order to determine the amino acid residues possibly responsible for binding and/ or catalytic activity. Furthermore, the structural features of the active site can be implemented into the existing small molecule substrate model, thus enhancing its predictive value with respect to possible metabolism by P450 2D6. As no crystal structures are yet available for membrane-bound P450s (such as P450 2D6), the crystal structures of bacterial (soluble) P450 101 (P450cam), P450 102 (P450BM3), and P450 108 (P450terp) have been used to build a three-dimensional model for P450 2D6 with molecular modeling techniques. Several important P450 2D6 substrates were consecutively docked into the active site of the protein model. The energy optimized positions of the substrates in the protein agreed well with the original relative positions of the substrates within the substrate model. This confirms the usefulness of small molecule models in the absence of structural protein data. Furthermore, the derived protein model indicates new leads for experimental validation and extension of the substrate model.

Bisallylic hydroxylation and epoxidation of polyunsaturated fatty acids by cytochrome P450

Polyunsaturated fatty acids can be oxygenated by cytochrome P450 to hydroxy and epoxy fatty acids. Two major classes of hydroxy fatty acids are formed by hydroxylation of the omega-side chain and by hydroxylation of bisallylic methylene carbons. Bisallylic cytochrome P450-hydroxylases transform linoleic acid to 11-hydroxylinoleic acid, arachidonic acid to 13-hydroxyeicosa-5Z,8Z,11Z,14Z-tetraenoic acid, 10-hydroxyeicosa-5Z,8Z,11Z,14Z-tetraenoic acid and 7-hydroxyeicosa-5Z,8Z,11Z,14Z-tetraenoic acid and eicosapentaenoic acid to 16-hydroxyeicosa-5Z,8Z,11Z,14Z,17Z-pent aenoic acid, 13-hydroxyeicosa-5Z,8Z,11Z,14Z,17Z-pent aenoic acid and 10-hydroxyeicosa-5Z,8Z,11Z,14Z,17Z-pent aenoic acid as major metabolites. The bisallylic hydroxy fatty acids are chemically unstable and decompose rapidly to cis-trans conjugated hydroxy fatty acids during acidic extractive isolation. Bisallylic hydroxylase activity appears to be augmented in microsomes induced by the synthetic glucocorticoid dexamethasone and by some other agents, but the P450 gene families of these hydroxylases have yet to be determined. The fatty acid epoxides, which are formed by cytochrome P450, are chemically stable, but are hydrolyzed to diols by soluble epoxide hydrolases. Epoxidation of polyunsaturated fatty acids is a prominent pathway of metabolism in the liver and the renal cortex and epoxy-genase activity appears to be under homeostatic control in the kidney. Many arachidonate epoxygenases have been identified belonging to the CYP2C gene subfamily. Epoxygenases have also been found in the central nervous system, endocrine organs, the heart and endothelial cells. Epoxides of arachidonic acid have been found to exert pharmacological effects on many cells.

The highly stereoselective oxidation of polyunsaturated fatty acids by cytochrome P450BM-3

Cytochrome P450BM-3 catalyzes NADPH-dependent metabolism of arachidonic acid to nearly enantiomerically pure 18(R)-hydroxyeicosatetraenoic acid and 14(S), 15(R)-epoxyeicosatrienoic acid (80 and 20% of total products, respectively). P450BM-3 oxidizes arachidonic acid with a rate of 3.2 +/- 0.4 micromol/min/nmol at 30 degrees C, the fastest ever reported for an NADPH-dependent, P450-catalyzed reaction. Fatty acid, oxygen, and NADPH are utilized in an approximately 1:1:1 molar ratio, demonstrating efficient coupling of electron transport to monooxygenation. Eicosapentaenoic and eicosatrienoic acids, two arachidonic acid analogs that differ in the properties of the C-15-C-18 carbons, are also actively metabolized by P450BM-3 (1.4 +/- 0.2 and 2.9 +/- 0.1 micromol/min/nmol at 30 degrees C, respectively). While the 17,18-olefinic bond of eicosapentaenoic acid is epoxidized with nearly absolute regio- and stereochemical selectivity to 17(S),18(R)-epoxyeicosatetraenoic acid (>/=99% of total products, 97% optical purity), P450BM-3 is only moderately regioselective during hydroxylation of the eicosatrienoic acid omega-1, omega-2, and omega-3 sp3 carbons, with 17-, 18-, and 19-hydroxyeicosatrienoic acid formed in a ratio of 2.4:2.2:1, respectively. Based on the above and on a model of arachidonic acid-bound P450BM-3, we propose: 1) the formation by P450BM-3 of a single oxidant species capable of olefinic bond epoxidation and sp3 carbon hydroxylation and 2) that product chemistry and, thus, catalytic outcome are critically dependent on active site spatial coordinates responsible for substrate binding and productive orientation between heme-bound active oxygen and acceptor carbon bond(s).

Putidaredoxin reductase-putidaredoxin-cytochrome p450cam triple fusion protein. Construction of a self-sufficient Escherichia coli catalytic system

Fusion proteins of cytochrome P450cam with putidaredoxin (Pd) and putidaredoxin reductase (PdR), the two proteins required to transfer electrons from NADH to P450cam, were constructed by fusing cDNAs encoding the three proteins in the expression vector pCWori+. Several fusion proteins, in which the order of the three protein domains and the linkers between them were varied, were expressed in Escherichia coli, purified, and characterized. The highest activity (kcat = 30 min-1) was obtained with a PdR-Pd-P450cam construct in which the peptides TDGTASS and PLEL were used, respectively, to link the PdR to the Pd and the Pd to the P450cam domains. Oxygen and NADH consumption is tightly coupled to substrate oxidation in the fusion proteins. The rate-limiting step in the catalytic turnover of these fusion proteins is electron transfer from Pd to P450cam. This is indicated by high rates of electron transfer from the PdR and Pd domains to exogenous electron acceptors, by an increase in the activity of the P450cam domain upon addition of exogenous Pd, and by the high activity of wild-type P450cam when incubated with a PdR-Pd fusion protein. E. coli cells expressing the PdR-Pd-P450cam fusion protein efficiently oxidize camphor to 5-exo-hydroxycamphor and 5-oxocamphor. E. coli cells expressing the triple fusion protein thus constitute the first heterologous self-sufficient catalytic system for the oxidation of camphor and other substrates by P450cam.

New applications of bacterial systems to problems in toxicology

Bacterial systems have long been of use in toxicology. In addition to providing general models of enzymes and paradigms for biochemistry and molecular biology, they have been adapted to practical genotoxicity assays. More recently, bacteria also have been used in the production of mammalian enzymes of relevance to toxicology. Escherichia coli has been used to express cytochrome P450, NADPH-cytochrome P450 reductase, flavin-containing monooxygenase, glutathione S-transferase, quinone reductase, sulfotransferase, N-acetyltransferase, UDP-glucuronosyl transferase, and epoxide hydrolase enzymes from humans and experimental animals. The expressed enzymes have been utilized in a variety of settings, including coupling with bacterial genotoxicity assays. Another approach has involved expression of mammalian enzymes directly in bacteria for use in genotoxicity systems. Particularly with Salmonella typhimurium. Applications include both the reversion mutagenesis assay and a system using a chimera with an SOS-response indicator and a reporter.

Identification of thromboxane A2 synthase active site residues by molecular modeling-guided site-directed mutagenesis

Human thromboxane A2 synthase (TXAS) exhibits spectral characteristics of cytochrome P450 but lacks monooxygenase activity. Its distinctive amino acid sequence makes TXAS the sole member of family 5 in the P450 superfamily. To better understand the structure-function relationship of this unusual P450, we have recently constructed a three-dimensional model for TXAS using P450BM-3 as the template (Ruan, K.-H., Milfeld, K., Kulmacz, R. J., and Wu, K. K. (1994) Protein Eng. 7, 1345-1551) and have identified a potential active site region. The catalytic roles of several putative active site residues were evaluated using selectively mutated recombinant TXAS expressed in COS-1 cells. Mutation of Ala-408 to Glu or Arg-413 to Gly led to a complete loss of enzyme activity despite expression of mutant protein levels equivalent to that of the wild-type TXAS. Mutation of Ala-408 to Gly or Leu retained the enzyme activity at levels of 30 or 40%, respectively. This suggests that Ala-408 provides a hydrophobic environment for substrate binding. Mutation of Arg-413 to Lys or Gln completely abolished the enzyme activity, indicating that this residue is essential to catalytic activity and supports its identification as an active site residue. Mutation of Arg-410 to Gly or Glu-433 to Ala resulted in >50% reduction in the enzyme activity without appreciably altering mutant protein expression, consistent with a more subtle effect of these residues on TXAS catalytic efficiency. Mutation of residues predicted to be involved in binding the heme prosthetic group, including the heme thiolate ligand Cys-480, Arg-478, Phe-127, and Asn-110, each markedly reduced the expressed protein level and abolished enzyme activity. This suggests that proper heme binding is important to synthesis or stability of recombinant TXAS. Mutation of Ile-346, which corresponds to P450cam-Thr-252, an essential amino acid involved in dioxygen bond scission, to Thr increased the enzymatic activity by 40%, suggesting that oxygen bond cleavage is not a rate-limiting step in thromboxane A2 biosynthesis. The present results from site-directed mutagenesis support the overall structure of the TXAS active site predicted by homology modeling and have allowed refinement of the position of bound substrate.

Clusters of charged residues in protein three-dimensional structures

Statistically significant charge clusters (basic, acidic, or of mixed charge) in tertiary protein structures are identified by new methods from a large representative collection of protein structures. About 10% of protein structures show at least one charge cluster, mostly of mixed type involving about equally anionic and cationic residues. Positive charge clusters are very rare. Negative (or histidine-acidic) charge clusters often coordinate calcium, or magnesium or zinc ions [e.g., thermolysin (PDB code: 3tln), mannose-binding protein (2msb), aminopeptidase (1amp)]. Mixed-charge clusters are prominent at interchain contacts where they stabilize quaternary protein formation [e.g., glutathione S-transferase (2gst), catalase (8act), and fructose-1,6-bisphosphate aldolase (1fba)]. They are also involved in protein-protein interaction and in substrate binding. For example, the mixed-charge cluster of aspartate carbamoyl-transferase (8atc) envelops the aspartate carbonyl substrate in a flexible manner (alternating tense and relaxed states) where charge associations can vary from weak to strong. Other proteins with charge clusters include the P450 cytochrome family (BM-3, Terp, Cam), several flavocytochromes, neuraminidase, hemagglutinin, the photosynthetic reaction center, and annexin. In each case in Table 2 we discuss the possible role of the charge clusters with respect to protein structure and function.

Investigation of the proton-assisted pathway to formation of the catalytically active, ferryl species of P450s by molecular dynamics studies of P450eryF

The recently determined crystal structure of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1) in its ferric heme substrate-bound form has been used to address one of the most fundamental unresolved aspects of the mechanism of oxidation common to this ubiquitous family of metabolizing heme proteins, the pathway from the twice reduced dioxygen species to the putative catalytically active ferryl oxygen species. Both of these species are too transient to have been characterized experimentally, and the transformation from one to the other has been only partially characterized. The observed requirement of two protons and the formation of water in this transformation suggests a proton-assisted dioxygen bond cleavage as a plausible pathway. However, this pathway is difficult to establish by experiment alone, and the source of the protons in the largely hydrophobic binding pocket of the P450s remains unclear. In this work we have performed molecular dynamics simulations of the twice reduced dioxygen substrate-bound form of this isozyme in order to (i) determine the plausibility of the proposed pathway to compound I formation, a proton-assisted cleavage of the dioxygen bond, and (ii) investigate the possible source of these protons. The analysis of the molecular dynamics trajectories of this species does indeed provide further evidence for this pathway and points to a source of protons. Specifically, two dynamically stable hydrogen bonds to the distal oxygen atom of the dioxygen ligand, one by the substrate and the other by a bound water, are found, consistent with the proposed proton-assisted cleavage of the bond and formation of water. In addition, an extensive dynamically stable hydrogen bond network is formed that connects the distal oxygen to Glu 360, a well-conserved residue in a channel accessible to solvent that could be the ultimate source of protons. The simulations were done for both a protonated and unprotonated Glu and led to a proposed mechanism of proton transfer by it to the distal oxygen atom. In order to validate the procedures used for the simulation of this transient twice-reduced species, we have used these same procedures to perform molecular dynamics simulations of two other forms of P450eryF, the ferric and ferryl substrate-bound species, and compared the results with experiment. The results for the ferric substrate-bound species were assessed by comparisons to the experimentally determined X-ray structure and fluctuations, and good agreement was found. The simulations performed for the ferryl substrate-bound species led to the correct prediction of the observed regio- and stereospecific hydroxylation of its natural substrate, 6-deoxyerythronolide B (6-DEB) at the 6S position. The results of these two additional studies lend credibility to the important mechanistic inferences from the simulations of the transient twice reduced dioxygen species: further evidence for a proton-assisted pathway from it to the catalytically active ferryl species and a possible source of the protons.

Substrate-assisted catalysis in cytochrome P450eryF

A highly conserved threonine in the active site of cytochromes P450 has been proposed to participate in O2 binding and cleavage. Cytochrome P450eryF is unusual in having alanine in place of this threonine and an ordered active site water molecule (Wat 519) which is hydrogen bonded to the substrate 5-hydroxyl group and is in position to operate as an acid catalyst required for cleaving dioxygen. To asses the role of this alanine residue and Wat 519 in catalysis, two mutant forms of P450eryF (Ala --> Ser,Ala --> Thr) and a substrate analogue lacking a 5-hydroxyl group were examined using kinetic, spectral and crystallographic techniques. In each case decreased catalytic activity was correlated with a loss or repositioning of Wat 519. These findings suggest that P450eryF utilizes the substrate to assist in the acid-catalysed dioxygen bond cleavage reaction.

Bacterial cytochromes P-450

The cytochromes P-450 (P-450s) constitute an extremely large family ('superfamily') of haemoproteins that catalyse the oxidation of a wide range of physiological and non-physiological compounds. A remarkable feature of the P-450s is the manipulation of the same basic structure and chemistry to achieve an enormous range of functions in organisms as diverse as bacteria and man. Indeed, the P-450s have been described as 'the most versatile biological catalyst known'. Much research is focussed on mammalian P-450s, with their roles in such processes as steroid transformations and the metabolism of carcinogens and other xenobiotics. However, our knowledge of the structure and function of the P-450s has been advanced by analysis of a limited number of its bacterial members, primarily P-450cam from Pseudomonas putida. Four P-450 structures have been solved to date, all of which are from bacterial sources. The aim of this review is to assess current knowledge of the many bacterial P-450s, with emphasis on their diverse biological roles and on the advances in our knowledge of this extremely important enzyme class, which have been made feasible through their study.

Influence of amino acid residue 374 of cytochrome P-450 2D6 (CYP2D6) on the regio- and enantio-selective metabolism of metoprolol

Cytochrome P-450 2D6 (CYP2D6) is an important human drug-metabolizing enzyme responsible for the oxidation of more than 30 widely used therapeutic agents. The enzymes encoded by the published genomic [Kimura, Umeno, Skoda, Meyer and Gonzalez (1989) Am. J. Hum. Genet. 45, 889-904] and cDNA [Gonzalez, Skoda, Kimura, Umeno, Zanger, Nebert, Gelboin, Hardwick and Meyer (1988) Nature 331, 442-446] sequences of CYP2D6, and presumed to represent wild-type sequences, differ at residue 374 and encode valine (CYP2D6-Val) and methionine (CYP2D6-Met) respectively. The influence of this amino acid difference on cytochrome P-450 expression, ligand binding, catalysis and stereoselective oxidation of metoprolol was investigated by the heterologous expression of the corresponding cDNAs in the yeast Saccharomyces cerevisiae. The level of expression of apo- and holo-protein was similar with each form of CYP2D6 cDNA, and the binding affinities of a series of ligands to CYP2D6-Val and CYP2D6-Met were identical. The enantioselective O-demethylation and alpha-hydroxylation of metoprolol were also similar with each form of CYP2D6, O-demethylation being R-(+)- enantioselective (CYP2D6-Val: R/S, 1.6; CYP2D6-Met: R/S, 1.4), whereas alpha-hydroxylation showed a preference for S-(-)-metoprolol (CYP2D6-Val: R/S, 0.7; CYP2D6-Met: R/S, 0.8). However, although the favoured regiomer overall was O-demethylmetoprolol (ODM), the regioselectivity for O-demethylation of each metoprolol enantiomer was significantly greater for CYP2D6-Val [R-(+)-: ODM/alpha-hydroxymetoprolol (alpha OH), 5.9; S-(-)-: ODM/alpha OH, 2.5) than that observed for CYP2D6-Met [R-(+)-: ODM/alpha OH, 2.2; S-(-)-: ODM/alpha OH, 1.4]. The stereoselective properties of CYP2D6-Val were consistent with those observed for CYP2D6 in human liver microsomes. The difference in the stereoselective properties of CYP2D6-Val and CYP2D6-Met were rationalized with respect to a homology model of the active site of CYP2D6 based on an alignment with the crystal structure of the bacterial cytochrome P-450BM-3' CYP102.

Function and membrane topology of wild-type and mutated cytochrome P-450c21

We have studied membrane topology of cytochrome P-450c21 (P450c21) using the approaches of mutagenesis and protease digestion. P450c21 is located at the cytoplasm with an N-terminal hydrophobic domain integrated into microsomal membranes. When this hydrophobic domain was replaced by a secretory signal peptide, P450c21 was translocated into the lumen and lost enzymic activity. No other topogenic sequence was detected in the bulk of the P450c21 peptide. A mutant protein with Pro-30 replaced by Leu (L30) corresponding to the mutation found in the diseased state was created. L30 protein lost 90% of enzymic activity, while a double mutant (L30R32) with an additional Leu-32 to Arg mutation had slightly higher residual enzymic activity. Apart from lower activity, L30 was also present in the cell at a lower level than wild-type P450c21. This lower level is probably due to increased degradation, as L30 is synthesized at a normal rate. Both L30 and L30R32 proteins, however, were integrated into membranes normally. Therefore the Pro-30 --> Leu mutation did not affect membrane integration, but affected the abundance and enzymic activity of P450c21.

Conserved Glu318 at the cytochrome P450 1A2 distal site is crucial in the nitric oxide complex stability

Nitric oxide synthase (NOS) has a thiolate-coordinated heme active site similar to that of cytochrome P450 (P450). Both NOS and P450 form stable nitric oxide (NO)-ferric heme complexes, whereas an NO-ferric heme complex of methemoglobin, that has an imidazole-coordinated heme active site, is easily reduced. The NO complex stability of the thiolate-coordinated hemoproteins, however, appeared irreconcilable with the strong electron-donating capability of the cysteine thiolate. In the present study, NO bindings to cytochrome P450 1A2 (P450 1A2) distal mutants were studied in the presence of various substrates. We found that a mutation at Glu-318 to Ala in the putative distal site of P450 1A2, suggested to be important in the O2 activation of P450 reactions, markedly facilitates the reduction of the NO-ferric complex. Addition of 1,2:3,4-dibenzanthracene or phenanthrene almost abolished the mutation effect on the NO complex. Based on these results, together with other spectral and kinetic data, it is suggested that the NO-ferric complex stability of P450, and perhaps of NOS, is largely ascribed to an ionic bridge between NO and the distal carboxyl group.

A model for human cytochrome P450 2D6 based on homology modeling and NMR studies of substrate binding

The cytochrome P450 responsible for the debrisoquine/sparteine polymorphism (P450 2D6) has been produced in large quantities by expression of a modified cDNA in baculovirus. A polyhistidine extension was incorporated at the C-terminus of the expressed protein, which, after purification of the protein on a nickel-agarose column, could be removed proteolytically by treatment with thrombin. Purified yields of P450 2D6 were 2.4 mg from 700 mL of cell culture. The protein had a greater than 90% heme content and was fully active, having no residual absorbance at 420 nm in the reduced CO complex. The quantities produced allowed direct study of the interaction of the substrate codeine with the enzyme by paramagnetic relaxation effects on the NMR spectrum of the substrate. Distances between the heme iron atom and substrate protons were calculated from these experiments, and the orientation of the substrate in the binding pocket was determined. This showed that codeine was bound with the methoxy group of the molecule closest to the heme iron (iron-methyl proton distance of 3.1 +/- 0.1 A), consistent with the observed O-demethylation to morphine. A model of the complex Of P450 2D6 with codeine was built from a multiple sequence and structure alignment of the known crystal structures for P450s, incorporating the experimental constraints derived from the NMR studies. This showed that the overall fold Of P450 2D6 is more similar to that of P450 BM3 than to either P450 cam or P450 terp. Codeine binds to P450 2D6 so that the methoxy group is directly above the A ring of the heme, while the basic nitrogen interacts with the carboxylate of aspartate 301.

Engineering a mineralocorticoid- to a glucocorticoid-synthesizing cytochrome P450

Site-directed mutagenesis of a domain (amino acids 299-338) aligning to the I-helix region of P450cam, P450BM3 and P450terp was used to investigate the different regioselectivities displayed in the hydroxylation reactions performed by human aldosterone synthase (P450aldo) and 11beta-hydroxylase (P45011beta). The two enzymes are 93% identical and are essential for the synthesis of mineralocorticoids and glucocorticoids in the human adrenal gland. Single replacement of P450aldo residues for P45011 beta-specific residues at positions 296, 301, 302, 320, and 335 only gave rise to slightly increased 11beta-hydroxylase activities. However, a L301P/A320V double substitution increased 11beta-hydroxylase activity to 60% as compared with that of P45011 beta. Additionally substituting Ala-320 for Val-320 of P45011 beta further enhanced this activity to 85%. The aldosterone synthase activities of the mutant P450aldo proteins were suppressed to a varying degree, with triple replacement mutant L301P/E302D/A320V retaining only 10% and double replacement mutant L301P/A320V retaining only 13% of the P450aldo wild type activity. These results demonstrate a switch in regio- and stereoselectivities of the engineered P450aldo enzyme due to manipulation of residues at three critical positions, and we attribute the determination of these features in P450aldo to the structure of a region analogous to the I-helix in P450cam.

Structural flexibility and functional versatility of cytochrome P450 and rapid evolution

P450 represents a large group of heme-thiolate enzymes that exhibit remarkably diverse activities for the metabolism of numerous endogenous and exogenous chemicals. Recent site-directed mutagenesis studies indicate that a single mutation at any of the key residues can be enough to alter the substrate and/or product specificities in the P450 activities. Molecular modeling predicts that these key residues are located within the substrate heme pocket. Structural elements involved in diversifying P450 activity appear to correspond to the B' helix, the F helix and the F/G interhelical loop in the bacterial P450s. Structures represented by these regions are extremely variable despite the fact that the core of the P450 substrate pocket is well conserved. A mutation within these regions may result in a significant geometrical alteration of the pocket and lead to diversify the P450 activity. Phylogenetical analysis shows a relatively high rate of nonsynonymous substitution within these substrate binding regions. The functional versatility of P450 can thus be largely accounted for in terms of pocket change brought about by rapid mutations.

Antigenic mapping of bacterial and animal cytochromes P-450

A peptide scanning (PEPSCAN) approach was used for antigenic mapping of two hepatic microsomal cytochromes P450 (rab1A2 and rab2B4) and the microbial cytochrome from Pseudomonas putida (P450 101 or P450cam). This approach includes simultaneous synthesis of pin-linked overlapping hexapeptides covering the whole sequences of three P450s and testing them by ELISA with corresponding polyclonal antisera. Microsomal cytochrome P450 maps were shown to vary depending on an antiserum used for testing the peptides, however, the most active linear B-epitopes were revealed with antisera from two animal species used. P450 linear B-epitopes were classified into individual and group-specific epitopes. While almost all P450 101 linear antigenic determinants are unique for this protein, rab1A2 and rab2B4 contain epitopes both individual for each protein, and subfamily- or even family-specific epitopes. These results point out the possibility of producing both monospecific and group-specific antipeptide antibodies against different P450s. The antigenic map of P450 101 was superimposed on the structural-functional map of this protein. Its linear B-epitopes were shown to coincide with boundaries of secondary structure elements, with surface-located, water accessible regions and with sites responsible for intermolecular interactions in the Pseudomonas putida monooxygenase system. Several known or predicted functionally active sites in microsomal cytochrome P450 rab1A2 and rab2B4 were also shown to coincide with linear B-epitopes. The peculiarities of epitope locations in the protein tertiary structure will allow to predict antigenic regions starting from protein structural information and vice versa, to structural protein models in accordance with antigenic mapping results. Antigenic regions which coincide with sites responsible for intermolecular interactions in monooxygenase systems may be synthesized as separate peptides and used as blockers of such interactions.

Conformational dynamics in cytochrome P450-substrate interactions

There now are four known cytochrome P450 crystal structures. Two of these, P450cam and P450eryF, are substrate-bound while P450terp and the heme domain of P450BM-3 are substrate-free. Here we describe a preliminary analysis of the P450BM-3 heme domain complexed with the 16-carbon fatty acid substrate, palmitoleic acid. A comparison of the substrate-free and -bound structures shows that a large conformational change in the substrate access channel accompanies substrate binding. This new information, together with the substrate-bound structures of P450cam and P450eryF, reveals which regions of P450 are the most important in controlling the dynamics of substrate binding and recognition.

NMR studies of recombinant cytochrome P450cam mutants

In the active center of cytochrome P450cam, Thr-252 is one of the conserved amino acid residues in the cytochrome P450 superfamily and plays a key role in the hydroxylation of camphor. T252A mutant, in which Thr-252 is replaced by alanine, consumed O2 at a rate comparable to that of the wild-type enzyme, whereas the amount of exo-5-hydroxycamphor formed was less than 10% of that formed by the wild-typed enzyme and H2O2 is the main product in the hydroxylation reaction. H2O2 was also yielded by the valine mutant and the consumption rate of O2 was much lower than that for the wild-type enzyme (Imai et al (1989) Proc Natl Acad Sci USA 86, 7823-7827). On the basis of the 1H- and 15N-NMR spectra, it was revealed that the anionic nature of the axial thiolate and the heme-environmental structures were substantially affected in the absence of d-camphor by the amino acid substitution at 252 Thr. In T252A mutant, however, the binding of camphor reduced these conformational alterations in the heme vicinity, probably due to the formation of interactions between camphor and enzyme. On the other hand, T252V mutant still exhibited large reduction of the anionic nature of the axial ligand in the presence of d-camphor and structural changes around heme were also enhanced, since the affinity of the valine mutant to d-camphor was low. These results imply that the hydrophobic and/or steric effects of the valine residue at 252 interfere with interactions around heme and camphor binding sites, which corresponds to the larger functional defects for T252V mutant.

Computer modeling of 3D structures of cytochrome P450s

The understanding of structure-function relationship of enzymes requires detailed information of their three-dimensional structure. Protein structure determination by X-ray and NMR methods, the two most frequently used experimental procedures, are often difficult and time-consuming. Thus computer modeling of protein structures has become an increasingly active and attractive option for obtaining predictive models of three-dimensional protein structures. Specifically, for the ubiquitous metabolizing heme proteins, the cytochrome P450s, the X-ray structures of four isozymes of bacterial origin, P450cam, P450terp, P450BM-3 and P450eryF have now been determined. However, attempts to obtain the structure of mammalian forms by experimental means have thus far not been successful. Thus, there have been numerous attempts to construct models of mammalian P450s using homology modeling methods in which the known structures have been used to various extents and in various strategies to build models of P450 isozymes. In this paper, we review these efforts and then describe a strategy for structure building and assessment of 3D models of P450s recently developed in our laboratory that corrects many of the weaknesses in the previous procedures. The results are 3D models that for the first time are stable to unconstrained molecular dynamics simulations. The use of this method is demonstrated by the construction and validation of a 3D model for rabbit liver microsomal P450 isozyme 2B4, responsible for the oxidative metabolism of diverse xenobiotics including widely used inhalation anesthetics. Using this 2B4 model, the substrate access channel, substrate binding site and plausible surface regions for binding with P450 redox partners were identified.

The crystal structure of chloroperoxidase: a heme peroxidase--cytochrome P450 functional hybrid

BACKGROUND

Chloroperoxidase (CPO) is a versatile heme-containing enzyme that exhibits peroxidase, catalase and cytochrome P450-like activities in addition to catalyzing halogenation reactions. The structure determination of CPO was undertaken to help elucidate those structural features that enable the enzyme to exhibit these multiple activities.

RESULTS

Despite functional similarities with other heme enzymes, CPO folds into a novel tertiary structure dominated by eight helical segments. The catalytic base, required to cleave the peroxide O-O bond, is glutamic acid rather than histidine as in other peroxidases. CPO contains a hydrophobic patch above the heme that could be the binding site for substrates that undergo P450-like reactions. The crystal structure also shows extensive glycosylation with both N- and O-linked glycosyl chains.

CONCLUSIONS

The proximal side of the heme in CPO resembles cytochrome P450 because a cysteine residue serves as an axial heme ligand, whereas the distal side of the heme is 'peroxidase-like' in that polar residues form the peroxide-binding site. Access to the heme pocket is restricted to the distal face such that small organic substrates can interact with the iron-linked oxygen atom which accounts for the P450-like reactions catalyzed by chloroperoxidase.

Evidence that aspartic acid 301 is a critical substrate-contact residue in the active site of cytochrome P450 2D6

Model building studies have intimated a role for aspartic acid 301 in the substrate binding of cytochrome P450 2D6 (CYP2D6). We have tested this hypothesis by generating a range of CYP2D6 mutants substituting a variety of amino acids at this site. The mutant proteins, which included substitution with a negatively charged glutamic acid residue or neutral asparagine, alanine, or glycine residues, were expressed in Saccharomyces cerevisiae. In addition, a mutant where aspartic acid 301 was deleted was also tested. All the mutants expressed approximately equivalent amounts of recombinant apoprotein and, apart from the alanine 301 and the aspartic acid 301 deletion mutants, gave carbon monoxide difference spectra of similar magnitude to the wild type. In the cases of the alanine and deletion mutants, the amount of holoprotein was significantly reduced or absent relative to the amount of apoprotein, indicating restricted heme incorporation. The glutamic acid mutant was shown to have similar catalytic properties to the wild type enzyme toward the substrates debrisoquine and metoprolol; however, some differences in regioselectivity and ligand binding were observed. The mutants containing neutral amino acids at position 301 exhibited marked reductions in catalytic activity. At low substrate concentrations little, if any, activity toward debrisoquine and metoprolol was measured. However, at a higher substrate concentration (2 mM) some activity was observed (about 10-20% of wild type levels). Consistent with the above findings, the debrisoquine-induced spin changes in the mutant proteins were markedly reduced. These data collectively demonstrate that aspartic acid 301 plays an important role in determining the substrate specificity and activity of CYP2D6 and provide experimental evidence supporting the role of this amino acid in forming an electrostatic interaction between the basic nitrogen atom in CYP2D6 substrates and the carboxylate group of aspartic acid 301.

Cytochrome P450

Since 1993, three new cytochrome P450 X-ray structures have been determined, giving a total of four known structures. Two of the new structures are in the substrate-free form and one is substrate-bound. These new structures, together with a wealth of mutagenesis studies on various P450s, have provided considerable information on what structural features control substrate specificity in P450s. In addition, some important insights into the catalytic mechanism have been made.

A molecular model for the interaction between vorozole and other non-steroidal inhibitors and human cytochrome P450 19 (P450 aromatase)

In a previous study (Vanden Bossche et al., Breast Cancer Res. Treat. 30 (1994) 43) the interaction between (+)-S-vorozole and the I-helix of cytochrome P450 19 (P450 aromatase) has been reported. In the present study we extended the "I-helix model" by incorporating the C-terminus of P450 aromatase. The crystal structures of P450 101 (P450 cam), 102 (P450 BM-3) and 108 (P450 terp) reveal that the C-terminus is structurally conserved and forms part of their respective substrate binding pocket. Furthermore, the present study is extended to the interaction between P450 aromatase and its natural substrate androstenedione and the non-steroidal inhibitors (-)-R-vorozole, (-)-S-fadrozole, R-liarozole and (-)-R-aminoglutethimide. It is found that (+)-S-vorozole, (-)-S-fadrozole and R-liarozole bind in a comparable way to P450 aromatase and interact with both the I-helix (Glu302 and Asp309) and C-terminus (Ser478 and His480). The weak activity of (-)-R-aminoglutethimide might be attributed to a lack of interaction with the C-terminus.

A three-dimensional model of aromatase cytochrome P450

P450 hemeproteins comprise a large gene superfamily that catalyzes monooxygenase reactions in the presence of a redox partner. Because the mammalian members are, without exception, membrane-bound proteins, they have resisted structure-function analysis by means of X-ray crystallographic methods. Among P450-catalyzed reactions, the aromatase reaction that catalyzes the conversion of C19 steroids to estrogens is one of the most complex and least understood. Thus, to better understand the reaction mechanism, we have constructed a three-dimensional model of P450arom not only to examine the active site and those residues potentially involved in catalysis, but to study other important structural features such as substrate recognition and redox-partner binding, which require examination of the entire molecule (excepting the putative membrane-spanning region). This model of P450arom was built based on a "core structure" identified from the structures of the soluble, bacterial P450s (P450cam, P450terp, and P450BM-P) rather than by molecular replacement, after which the less conserved elements and loops were added in a rational fashion. Minimization and dynamic simulations were used to optimize the model and the reasonableness of the structure was evaluated. From this model we have postulated a membrane-associated hydrophobic region of aliphatic and aromatic residues involved in substrate recognition, a redox-partner binding region that may be unique compared to other P450s, as well as residues involved in active site orientation of substrates and an inhibitor of P450arom, namely vorozole. We also have proposed a scheme for the reaction mechanism in which a "threonine switch" determines whether oxygen insertion into the substrate molecule involves an oxygen radical or a peroxide intermediate.

Three-dimensional models of human and other mammalian microsomal P450s constructed from an alignment with P450102 (P450bm3)

1. A novel modelling alignment for P450s, utilizing NADPH-P450 reductase for electron transfer, is proposed on the basis of analysis of their amino acid sequences. 2. Information used to facilitate the alignment process includes: the recent X-ray crystal structure of P450102 (P450bm3), site-directed mutagenesis experiments, chemical modification of specific residues, and antibody recognition studies. 3. The alignment has been used to construct a number of microsomal P450s of relevance to xenobiotic and endogenous metabolism.

Expression of a recombinant derivative of human aromatase P450 in insect cells utilizing the baculovirus vector system

Aromatase P450 (P450arom) is the enzyme responsible for estrogen biosynthesis. Studies of the relationship of the function of this enzyme to its structure have been hampered by lack of a suitable preparation. In the present report we describe the expression of a recombinant derivative of P450arom in insect cells by means of the baculovirus vector system. This protein, which lacks the first 41 amino acids from the N-terminus, and hence the membrane-spanning region, has spectral properties and activity similar to that of the wildtype protein. Moreover, the presence of a hexameric histidine tag at the C-terminus permits its facile purification by means of nickel-agarose affinity chromatography. This system permits the synthesis of quantities of a biologically active derivative of P450arom suitable for studies designed to explore the relationship of function to structure.

Structure of cytochrome P450eryF involved in erythromycin biosynthesis

Cytochrome P450eryF catalyzes the 6S-hydroxylation of 6-deoxyerythronolide B, the initial reaction in a multistep pathway to convert 6-deoxyerythronolide B into the antibiotic, erythromycin. The overall structure of P450eryF is similar to that of P450cam but differs in the exact positioning of several alpha-helices. The largest difference occurs in the B' helix and results in the enlargement of the substrate-binding pocket of P450eryF. The substrate is positioned with the macrolide ring perpendicular to the haem plane and contacts seven hydrophobic residues and three solvent molecules. The substrate participates in a network of hydrogen bonds that may provide a proton shuttle pathway in the oxygen cleavage reaction.

Identification by in vitro mutagenesis of the interaction of two segments of C2MstC1, a chimera of cytochromes P450 2C2 and P450 2C1

A hybrid cytochrome P450, C2MstC1, with 306 N-terminal amino acids derived from cytochrome P450 2C2 sequence and 184 C-terminal amino acids from cytochrome P450 2C1 acquires a novel progesterone 21-hydroxylase activity which is absent in the parent enzymes. Extension of the cytochrome P450 2C2 sequence to residue 382 reduced progesterone hydroxylase activity to 5% of that of C2MstC1, while further extension to residue 411 or 462 increased activity back to about 30 or 40%, respectively. In the chimera with cytochrome P450 2C2 sequence to residue 382, substitution of cytochrome P450 2C1 amino acids at positions 368, 369, and 374 increased progesterone hydroxylase activity to a level equivalent to that of C2MstC1. In the chimera with cytochrome P450 2C2 sequence extending to residue 411, substitutions of P450 2C1 amino acids at positions 386 and 388, in addition those at 368, 369, and 374, were required to obtain activities equivalent to that of C2MstC1, which suggests an interaction between these two regions. The lauric acid hydroxylase activities of all chimeras and mutant cytochromes P450 differed by 2-fold or less, demonstrating that the changes in progesterone hydroxylase activity reflected altered interactions with the substrate. Alignment of cytochrome P450 2C1 sequence with cytochromes P450cam, P450BM-3, and P450terp predicts that residues 368/369 and 386/388 are in adjacent antiparallel strands of the same beta-sheet, in agreement with the experimental data suggesting an interaction between these two regions.

Structure and function of cytochromes P450: a comparative analysis of three crystal structures

BACKGROUND

Cytochromes P450 catalyze the oxidation of a variety of hydrophobic substrates. Sequence identities between P450 families are generally low (10-30%), and consequently, the structure-function correlations among P450s are not clear. The crystal structures of P450terp and the hemoprotein domain of P450BM-3 were recently determined, and are compared here with the previously available structure of P450cam.

RESULTS

The topology of all three enzymes is quite similar. The heme-binding core structure is well conserved, except for local differences in the I helices. The greatest variation is observed in the substrate-binding regions. The structural superposition of the proteins permits an improved sequence alignment of other P450s. The charge distribution in the three structures is similarly asymmetric and defines a molecular dipole.

CONCLUSIONS

Based on this comparison we believe that all P450s will be found to possess the same tertiary structure. The ability to precisely predict other P450 substrate-contact residues is limited by the extreme structural heterogeneity in the substrate-recognition regions. The central I-helix structures of P450terp and P450BM-3 suggest a role for helix-associated solvent molecules as a source of catalytic protons, distinct from the mechanism for P450cam. We suggest that the P450 molecular dipole might aid in both redox-partner docking and proton recruitment for catalysis.

Arylhydrazines as probes of hemoprotein structure and function

The reactions of arylhydrazines (ArNHNH2) or aryldiazenes (ArN = NH) with simple iron porphyrins or with hemoproteins that have relatively open active sites, including hemoglobin, myoglobin, cytochrome P450, chloroperoxidase, catalase, prostaglandin synthase, and indoleamine-2,3-dioxygenase yield sigma-bonded aryl-iron complexes. Denaturation of the protein complexes under aerobic, acidic conditions shifts the aryl group to the porphyrin nitrogens and produces mixtures of the four possible N-arylprotoporphyrin IX regioisomers. The regioisomers are obtained in approximately equal amounts if the iron-to-nitrogen shift occurs outside of the protein but the ratio of isomers differs if the rearrangement is controlled by the protein. Only in the case of cytochrome P450 enzymes can the shift be induced to occur without denaturation of the protein. The isomer ratios obtained when the shift occurs in the intact active site provide direct experimental information on the active site topology and dynamics. Topological information has thus been obtained for cytochromes P450 1A1, 1A2, 2B1, 2B2, 2B4, 2B10, 2B11, 2E1, 11A1, 51, 101, 102, and 108. In contrast to hemoproteins with open active sites, conventional peroxidases react with arylhydrazines to give delta-meso-aryl adducts and covalent protein adducts. Reaction with the delta-meso edge but not the heme iron provides key evidence that restricting access of substrates to the ferryl oxygen helps direct the reaction towards peroxidase rather than peroxygenase catalysis.

Chimeras of the human cytochrome P450 1A family produced in yeast. Accumulation in microsomal membranes, enzyme kinetics and stability

An expression library of hybrid cDNAs was constructed in vivo by homeologous recombination in yeast between human P450 1A1 and P450 1A2 sequences. Two clones exhibiting highly enhanced monooxygenase activities in vivo were selected. Chimera S12 includes the 88 N-terminal residues of P450 1A1 fused to the complementary part of the P450 1A2 sequence. Chimera S71 derives from P450 1A1 by the substitution of the 36 C-terminal amino acid residues by the corresponding 38 residues of the 1A2 sequence. Biochemical analysis on microsomal fractions indicated that S12 and S71 have the same substrate specificities as 1A2 and 1A1, respectively. The observed increase in the in vivo monooxygenase activity is related to a ninefold increase in the microsomal S12 content as compared to the 1A2 content. In contrast, the expression level of S71 is slightly reduced but its turnover numbers are increased as compared to 1A1. The folding stability of chimeric P450 enzymes was evaluated by thermal and chaotropic agent denaturation. No difference was found between S12 and 1A2, but S71 appeared slightly less stable than 1A1. In vivo experiments indicated that S12 mRNA accumulation and stability are quite similar to the stability of parental 1A2 and, for both chimeras and parental enzymes, the protein half-lives are longer than the cell doubling time. The surprising accumulation of chimera S12 in the microsomal membrane is discussed in terms of the relationship of protein folding with transport to the endoplasmic reticulum membrane and the apparent expression levels of human P450 enzymes produced in yeast.

Structural variation in heme enzymes: a comparative analysis of peroxidase and P450 crystal structures

New peroxidase and cytochrome P450 structures have been solved recently, making more detailed comparisons of the two types of enzyme possible.