Primary and secondary structural patterns in eukaryotic cytochrome P-450 families correspond to structures of the helix-rich domain of Pseudomonas putida cytochrome P-450cam. Indications for a similar overall topology

An extensive sequence analysis of the eukaryotic cytochrome P-450 (P-450) protein families was conducted with a view to identifying conserved regions that might be related to secondary structural features in the Pseudomonas putida camphor hydroxylase (P-450cam). All sequences available on-line were collected, classified and aligned within families. Distinctively different sequences were chosen from each of seven eukaryotic families, and an unbiased multi-alignment was constructed. Profile patterns of the most conserved regions were generated and screened against the sequence of P-450cam, the structure of which has been elucidated by X-ray crystallography. While some of these profiles did not map on the P-450cam sequence, the structurally most important helices were clearly identified and the correlations were found to be statistically significant. Our analysis suggests that the helix-rich domain with the cysteine pocket and the oxygen-binding site is conserved in all P-450 forms. Helices I and L from P-450cam can be easily identified in all eukaryotic P-450 forms. Other helices which seem to exist in all P-450 forms include helices C, D, G and J. K. In the helix-poor domain of P-450cam, only structures b3/b4 seem to have been conserved. The obvious sequence conservation throughout the helix-rich domain of the P-450cam protein might be expected for a molecular class whose overall topology is preserved. Additional support for the conservation of structure between eukaryotic cytochromes P-450 and P-450cam comes from secondary structure prediction of the eukaryotic sequences.     

Identification and location of alpha-helices in mammalian cytochromes P450

A model of the alpha-helical structure of mammalian cytochromes P450 is proposed. The location and sequence of alpha-helices in mammalian cytochromes P450 were predicted from their homology with those of cytochrome P450cam, and these sequences were generally confirmed as helical in nature by using a secondary structure prediction method. These analyses were applied to 26 sequences in 6 gene families of cytochrome P450. Mammalian cytochromes P450 consist of approximately 100 amino acid residues more than cytochrome P450cam. This difference was accounted for by three major areas of insertion: (1) at the N-terminus, (2) between helices C and D and between helices D and E, and (3) between helices J and K. Insertion 1 has been suggested by others as a membrane anchoring sequence, but the apparent insertions at 2 and 3 are novel observations; it is suggested that they may be involved in the binding of cytochrome P450 reductase. Only the mitochondrial cytochrome P450 family appeared to show a major variation from this pattern, as insertion 2 was absent, replaced by an insertion between helices G and H and between helices H and I. This may reflect the difference in electron donor proteins that bind to members of this cytochrome P450 family. Other than these differences the model of mammalian cytochromes P450 proposed maintains the general structure of cytochrome P450cam as determined by its alpha-helical composition.     

Secondary structure and membrane topology of cytochrome P450s

The secondary structure prediction of 19 microsomal cytochrome P450s from two different families was made on the basis of their amino acid sequences. It was shown that there is structural similarity between the heme-binding sites in these enzymes and those in the bacterial P450cam. An average predicted secondary structure of cytochrome P450 proteins with 70% accuracy contains about 46% alpha-helices, 12% beta-sheets, 9% beta-turns, and 33% random coils. In the region of residues 35-120 in microsomal P450s two adjacent beta alpha beta-units (the Rossmann domain), were recognized and may be available to interact with the NADPH-cytochrome P450 reductase. Using the procedure for identification of hydrophobic and membrane-associated alpha-helical segments, only one N-terminal transmembrane anchor was predicted. Also the heme-binding site may include the surface-bound helix. A model for vertebrate microsomal P450s having an amphipathic membrane protein located on the cytoplasmic side of the endoplasmic reticulum membrane, with their active center lying outside or on the bilayer border, is proposed.     

Role of protein and substrate dynamics in catalysis by Pseudomonas putida cytochrome P450cam

The role of protein structural flexibility and substrate dynamics in catalysis by cytochrome P450 enzymes is an area of current interest. We have addressed these in cytochrome P450(cam) (P450(cam)) and its Y96A mutant with camphor and its related compounds using fluorescence spectroscopy. Previously [Prasad et al. (2000) FEBS Lett. 477, 157-160], we provided experimental support to dynamic fluctuations in P450(cam), and substrate access into the active site region via the channel next to the flexible F-G helix-loop-helix segment. In the investigation described here, we show that the dynamic fluctuations in the enzyme are substrate dependent as reflected by tryptophan fluorescence quenching experiments. The orientation of tryptophan relative to heme (kappa(2)) for W42 obtained from time-resolved tryptophan fluorescence measurements show variation with type of substrate bound to P450(cam) suggesting regions distant from heme-binding site are affected by physicochemical and steric characteristics/protein-substrate interactions of P450(cam) active site. We monitored substrate dynamics in the active site region of P450(cam) by time-resolved substrate anisotropy measurements. The anisotropy decay of substrates bound to P450(cam) indicate that mobility of substrates is modulated by physicochemical and steric characteristics/protein-substrate interactions of local active site structure, and provides an understanding of factors controlling observed hydroxylated products for substrate bound P450(cam) complexes. The present study shows that P450(cam) local and peripheral structural flexibility and heterogeneity along with substrate mobility play an important role in regulating substrate binding orientation during catalysis and accommodating diverse range of substrates within P450(cam) heme pocket.     

Molecular modeling of the 3-D structure of cytochrome P-450scc.

Sequence-alignment studies of the bovine mitochondrial cholesterol side-chain cleavage enzyme cytochrome P-450scc with the bacterial cytochrome P-450cam (camphor hydroxylating enzyme) have been undertaken. Our novel alignment of the sequences revealed 69 identical residues and many highly conserved regions. The results of the sequence alignment studies were used to model the 3-D structure of P-450scc based on the available crystal structure of P-450cam. The major insertions in the sequence are found mainly on four external-loop regions of the molecule, while the core structure of P-450cam is retained with subtle internal modifications. The most hydrophobic of these four external loops is proposed as a candidate for membrane attachment.     

Crystal structure of P450cin in a complex with its substrate, 1,8-cineole, a close structural homologue to D-camphor, the substrate for P450cam

Cytochrome P450cin catalyzes the monooxygenation of 1,8-cineole, which is structurally very similar to d-camphor, the substrate for the most thoroughly investigated cytochrome P450, cytochrome P450cam. Both 1,8-cineole and d-camphor are C(10) monoterpenes containing a single oxygen atom with very similar molecular volumes. The cytochrome P450cin-substrate complex crystal structure has been solved to 1.7 A resolution and compared with that of cytochrome P450cam. Despite the similarity in substrates, the active site of cytochrome P450cin is substantially different from that of cytochrome P450cam in that the B' helix, essential for substrate binding in many cytochrome P450s including cytochrome P450cam, is replaced by an ordered loop that results in substantial changes in active site topography. In addition, cytochrome P450cin does not have the conserved threonine, Thr252 in cytochrome P450cam, which is generally considered as an integral part of the proton shuttle machinery required for oxygen activation. Instead, the analogous residue in cytochrome P450cin is Asn242, which provides the only direct protein H-bonding interaction with the substrate. Cytochrome P450cin uses a flavodoxin-like redox partner to reduce the heme iron rather than the more traditional ferredoxin-like Fe(2)S(2) redox partner used by cytochrome P450cam and many other bacterial P450s. It thus might be expected that the redox partner docking site of cytochrome P450cin would resemble that of cytochrome P450BM3, which also uses a flavodoxin-like redox partner. Nevertheless, the putative docking site topography more closely resembles cytochrome P450cam than cytochrome P450BM3.     

Effect of alcohols on binding of camphor to cytochrome P450cam: spectroscopic and stopped flow transient kinetic studies

Addition of alcohols to cytochrome P450cam (CYP101) was shown to release the substrate camphor from the heme pocket of the enzyme. The release of the substrate was found to be caused both due to increased solubility of the substrate in solution in presence of alcohol and due to change in the tertiary structure of the active site of the enzyme. The far-UV CD and near-UV CD spectra reveal that addition of alcohols to cytochrome P450cam cause a small change in the secondary structural elements but a significant change in the tertiary structural organization of this enzyme. The CD spectra at the heme region at various concentrations of alcohols indicate a substantial change in the tertiary structural organization around the heme moiety too. The equilibrium constant associated with the binding of camphor to Cyt P450cam is strongly dependent on the concentration of alcohols and the corresponding free energy associated with the binding is found to scale linearly with the concentration of alcohols. Kinetic experiments on binding of camphor to Cyt P450cam show that both k(on) and k(off) rate constants are strongly affected by addition of alcohols suggesting that alcohol expel camphor out of the heme cavity of Cyt P450cam by affecting tertiary structure of Cyt P450cam as well as by modifying the solubility properties of camphor in aqueous medium.     

Continuous B-epitope maps of cytochrome P450cam (CYP101) obtained by peptide scanning: correlation to spatial structure

Protein continuous B-epitopes can be revealed using short synthetic peptides that overlap a known protein sequence. Since the whole protein surface is considered to possess antigenic properties, a question that arises is whether a set of linear B-epitopes determined by peptide scanning correlates with a protein spatial structure. We have chosen cytochrome P450cam (CYP101) of Pseudomonas putida, with known 3D structure, as a template. Sera of two rabbits and antibody egg yolk preparations from three chickens were produced against the P450cam molecule. These polyclonals were analyzed separately in ELISA with 409 overlapping P450cam hexapeptides. The whole set of continuous antigenic sites of P450cam covered about 45% of the P450cam sequence. However, immunodominant sites (those revealed with more than 50% antibody preparations), the so-called "antigenic core," represent only 9% of the protein sequence. While the amount of water-accessible residues in the total antigenic map (42%) was close to that in the whole native P450cam molecule (39%), the amount of water-accessible residues in the antigenic core was significantly higher (64%). These results led to the conclusion that antigenic core epitopes can be associated to the molecular surface, whereas epitopes with low detection frequency may partly correspond to unfolded regions of the protein molecule.     

How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms

Cytochrome P450s form a ubiquitous protein family with functions including the synthesis and degradation of many physiologically important compounds and the degradation of xenobiotics. Cytochrome P450cam from Pseudomonas putida has provided a paradigm for the structural understanding of cytochrome P450s. However, the mechanism by which camphor, the natural substrate of cytochrome P450cam, accesses the buried active site is a long-standing puzzle. While there is recent crystallographic and simulation evidence for opening of a substrate-access channel in cytochrome P450BM-3, for cytochrome P450cam, no such conformational changes have been observed either in different crystal structures or by standard molecular dynamics simulations. Here, a novel simulation method, random expulsion molecular dynamics, is presented, in which substrate-exit channels from the buried active site are found by imposing an artificial randomly oriented force on the substrate, in addition to the standard molecular dynamics force field. The random expulsion molecular dynamics method was tested in simulations of the substrate-bound structure of cytochrome P450BM-3, and then applied to complexes of cytochrome P450cam with different substrates and with product. Three pathways were identified, one of which corresponds to a channel proposed earlier on the basis of crystallographic and site-directed mutagenesis data. Exit via the water-filled channel, which was previously suggested to be a product exit channel, was not observed. The pathways obtained by the random expulsion molecular dynamics method match well with thermal motion pathways obtained by an analysis of crystallographic B-factors. In contrast to large backbone motions (up to 4 A) observed in cytochrome P450BM-3 for the exit of palmitoleic acid, passage of camphor through cytochrome P450cam only requires small backbone motions (less than 2.4 A) in conjunction with side-chain rotations. Concomitantly, in almost all the exit trajectories, salt-links that have been proposed to act as ionic tethers between secondary structure elements of the protein, are perturbed.     

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide‐range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of protein classification are based on the canonical paradigm that attributes proteins’ function to their three‐dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in C_α atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.     

Cholesterol Ester Oxidation by Mycobacterial Cytochrome P450

Mycobacteria share a common cholesterol degradation pathway initiated by oxidation of the alkyl side chain by enzymes of cytochrome P450 (CYP) families 125 and 142. Structural and sequence comparisons of the two enzyme families revealed two insertions into the N-terminal region of the CYP125 family (residues 58–67 and 100–109 in the CYP125A1 sequence) that could potentially sterically block the oxidation of the longer cholesterol ester molecules. Catalytic assays revealed that only CYP142 enzymes are able to oxidize cholesteryl propionate, and although CYP125 enzymes could oxidize cholesteryl sulfate, they were much less efficient at doing so than the CYP142 enzymes. The crystal structure of CYP142A2 in complex with cholesteryl sulfate revealed a substrate tightly fit into a smaller active site than was previously observed for the complex of CYP125A1 with 4-cholesten-3-one. We propose that the larger CYP125 active site allows for multiple binding modes of cholesteryl sulfate, the majority of which trigger the P450 catalytic cycle, but in an uncoupled mode rather than one that oxidizes the sterol. In contrast, the more unhindered and compact CYP142 structure enables enzymes of this family to readily oxidize cholesteryl esters, thus providing an additional source of carbon for mycobacterial growth.     

Active site topologies of bacterial cytochromes P450101 (P450cam), P450108 (P450terp), and P450102 (P450BM-3). In situ rearrangement of their phenyl-iron complexes.

The reactions of cytochromes P450101 (P450cam), P450108 (P450terp), and P450102 (P450BM-3) with phenyldiazene result in the formation of phenyl-iron complexes with absorption maxima at 474-478 nm. Treatment of the cytochrome P450 complexes with K3Fe(CN)6 decreases the 474-478 nm absorbance and shifts the phenyl group from the iron to the porphyrin nitrogens. Acidification and extraction of the prosthetic group from each of the ferricyanide-treated enzymes yields a different mixture of the four possible N-phenylprotoporphyrin IX regioisomers. The ratios of the regioisomers with the phenyl ring on pyrrole rings B, A, C, and D (in order of elution from the high performance liquid chromatography column) are, respectively: cytochrome P450cam, 0:0:1:4; P450terp, 0:0:0:1; and P450BM-3, 2:10:2:1. The isomer ratio for recombinant cytochrome P450BM-3 without the cytochrome P450 reductase domain (2:9:2:1) shows that the reductase domain does not detectably perturb the active site topology of cytochrome P450BM-3. Potassium ions modulate the intensity of the spectrum of the phenyl-iron complex of cytochrome P450cam, but do not alter the N-phenyl isomer ratio. Computer graphics analysis of the crystal structure of the cytochrome P450cam phenyl-iron complex indicates that the active site of cytochrome P450cam is open above pyrrole ring D and, to a small extent, pyrrole ring C, in complete agreement with the observed N-phenylprotoporphyrin IX regioisomer pattern. The regioisomer ratios indicate that the active site of cytochrome P450terp is only open above pyrrole ring D, whereas that of cytochrome P450BM-3 is open to some extent above all the pyrrole rings but particularly above pyrrole ring A. The bacterial enzymes thus have topologies distinct from each other and from those of the mammalian enzymes so far investigated, which have active sites that are open to a comparable extent above pyrrole rings A and D.     

Protein engineering of cytochrome p450(cam) (CYP101) for the oxidation of polycyclic aromatic hydrocarbons

Mutations of the active site residues F87 and Y96 greatly enhanced the activity of cytochrome P450(cam) (CYP101) from Pseudomonas putida for the oxidation of the polycyclic aromatic hydrocarbons phenanthrene, fluoranthene, pyrene and benzo[a]pyrene. Wild-type P450(cam) had low (<0.01 min(-1)) activity with these substrates. Phenanthrene was oxidized to 1-, 2-, 3- and 4-phenanthrol, while fluoranthene gave mainly 3-fluoranthol. Pyrene was oxidized to 1-pyrenol and then to 1,6- and 1,8-pyrenequinone, with small amounts of 2-pyrenol also formed with the Y96A mutant. Benzo[a]pyrene gave 3-hydroxybenzo[a]pyrene as the major product. The NADH oxidation rate of the mutants with phenanthrene was as high as 374 min(-1), which was 31% of the camphor oxidation rate by wild-type P450(cam), and with fluoranthene the fastest rate was 144 min(-1). The oxidation of phenanthrene and fluoranthene were highly uncoupled, with highest couplings of 1.3 and 3.1%, respectively. The highest coupling efficiency for pyrene oxidation was a reasonable 23%, but the NADH turnover rate was slow. The product distributions varied significantly between mutants, suggesting that substrate binding orientations can be manipulated by protein engineering, and that genetic variants of P450(cam) may be useful for studying the oxidation of polycyclic aromatic hydrocarbons by P450 enzymes.     

High-resolution crystal structure of cytochrome P450cam

The crystal structure of Pseudomonas putida cytochrome P450cam with its substrate, camphor, bound has been refined to R = 0.19 at a normal resolution of 1.63 A. While the 1.63 A model confirms our initial analysis based on the 2.6 A model, the higher resolution structure has revealed important new details. These include a more precise assignment of sequence to secondary structure, the identification of three cis-proline residues, and a more detailed picture of substrate-protein interactions. In addition, 204 ordered solvent molecules have been found, one of which appears to be a cation. The cation stabilizes an unfavorable polypeptide conformation involved in forming part of the active site pocket, suggesting that the cation may be the metal ion binding site associated with the well-known ability of metal ions to enhance formation of the enzyme-substrate complex. Another unusual polypeptide conformation forms the proposed oxygen-binding pocket. A localized distortion and widening of the distal helix provides a pocket for molecular oxygen. An intricate system of side-chain to backbone hydrogen bonds aids in stabilizing the required local disruption in helical geometry. Sequence homologies strongly suggest a common oxygen-binding pocket in all P450 species. Further sequence comparisons between P450 species indicate common three-dimensional structures with changes focused in a region of the molecule postulated to be associated with the control of substrate specificity.     

Substrate Access to Cytochrome P450cam: A Comparison of a Thermal Motion Pathway Analysis with Molecular Dynamics Simulation Data

Cytochrome P450 enzymes are hemoprotein monooxygenases that catalyse the oxidation of a variety of compounds. The mechanism by which camphor, the natural substrate of Cytochrome P450cam (P450cam), accesses the active site is a long-standing puzzle, although putative access channels have been proposed. A thermal motion pathway (TMP) analysis was performed on the crystal structure of P450cam with camphor bound. Hereby, three distinct thermal motion pathway families (TMPFs) were found. Possible substrate access channels obtained by this analysis based on B-factors are compared with exit channels explored by molecular dynamics simulations (MDS) by imposing an artificial expulsion force on the substrate in addition to the standard MD force field. Two out of three TMPFs are supported by results obtained with the random expulsion MDS method. However, the pathway found by the TMP method to have the highest average B-factor could not be observed by MDS. The pathway proposed from crystallographic data, which is a small opening above the active site located near residues 185, 87 and 395 corresponds to the TMPF with the second highest average B-factor.     

Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences.

The substrate recognition regions in cytochrome P450 family 2 (CYP2) proteins were inferred by group-to-group alignment of CYP2 sequences and those of bacterial P450s, including Pseudomonas putida P450 101A (P450cam), whose substrate-binding residues have been definitely identified by x-ray crystallography of a substrate-bound form (Poulos T. L., Finzel, B. C., and Howard, A. J. (1987) J. Mol. Biol. 195, 687-700). The six putative substrate recognition sites, SRSs, thus identified are dispersively located along the primary structure and constitute about 16% of the total residues. All the reported point mutations and chimeric fragments that significantly affect the substrate specificities of the parental CYP2 enzymes fell within or overlapped some of the six SRSs. Analysis of nucleotide substitution patterns in closely related members in four subfamilies, CYP2A, 2B, 2C, and 2D, consistently indicated that the SRSs have accumulated more nonsynonymous (amino acid-changing) substitutions than the rest of the sequence. This observation supports the idea that diversification of duplicate genes of drug-metabolizing P450s occurs primarily in substrate recognition regions to cope with an increasing number of foreign compounds.     

Common active site architecture and binding strategy of four phenylpropanoid P450s from Arabidopsis thaliana as revealed by molecular modeling

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.     

Resonance Raman studies of cytochrome P450BM3 and its complexes with exogenous ligands.

Resonance Raman spectra are reported for both the heme domain and holoenzyme of cytochrome P450BM3 in the resting state and for the ferric NO, ferrous CO, and ferrous NO adducts in the absence and presence of the substrate, palmitate. Comparison of the spectrum of the palmitate-bound form of the heme domain with that of the holoenzyme indicates that the presence of the flavin reductase domain alters the structure of the heme domain in such a way that water accessibility to the distal pocket is greater for the holoenzyme, a result that is consistent with analogous studies of cytochrome P450cam. The data for the exogenous ligand adducts are compared to those previously reported for corresponding derivatives of cytochrome P450cam and document significant and important differences for the two proteins. Specifically, while the binding of substrate induces relatively dramatic changes in the nu(Fe-XY) modes of the ferrous CO, ferric NO, and ferrous NO derivatives of cytochrome P450cam, no significant changes are observed for the corresponding derivatives of cytochrome P450BM3 upon binding of palmitate. In fact, the spectral data for substrate-free cytochrome P450BM3 provide evidence for distortion of the Fe-XY fragment, even in the absence of substrate. This apparent distortion, which is nonexistent in the case of substrate-free cytochrome P450cam, is most reasonably attributed to interaction of the Fe-XY fragment with the F87 phenylalanine side chain. This residue is known to lie very close to the heme iron in the substrate-free derivative of cytochrome P450BM3 and has been suggested to prevent hydroxylation of the terminal, omega, position of long-chain fatty acids.     

A molecular model for the enzyme cytochrome P450(17 alpha), a major target for the chemotherapy of prostatic cancer

The enzyme cytochrome P450(17 alpha) catalyses two key steps in the biosynthesis of the androgens from pregnanes: the 17 alpha hydroxylation step and the subsequent 17-20 lyase reaction. Using a variety of techniques, including sequence alignment, secondary structure prediction, molecular mechanics and molecular dynamics, we have constructed a model for the three-dimensional structure of P450(17 alpha) based on that of P450cam, the only cytochrome P450 enzyme for which the crystal structure is known. The model suggests the possibility of two modes of binding of steroid substrates at the active site, perhaps reflecting the dual functionality of the enzyme.     

Structural evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.