Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross‐substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily‐specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

Spectroelectrochemistry of cytochrome P450cam

The spectroelectrochemistry of camphor-bound cytochrome P450cam (P450cam) using gold electrodes is described. The electrodes were modified with either 4,4(')-dithiodipyridin or sodium dithionite. Electrolysis of P450cam was carried out when the enzyme was in solution, while at the same time UV-visible absorption spectra were recorded. Reversible oxidation and reduction could be observed with both 4,4(')-dithiodipyridin and dithionite modified electrodes. A formal potential (E(0')) of -373mV vs Ag/AgCl 1M KCl was determined. The spectra of P450cam complexed with either carbon monoxide or metyrapone, both being inhibitors of P450 catalysis, clearly indicated that the protein retained its native state in the electrochemical cell during electrolysis.     

Ethylbenzene hydroxylation by cytochrome P450cam.

The metabolism of ethylbenzene by cytochrome P450cam was analyzed by experimental and theoretical methods. The present experiments indicate that ethylbenzene is hydroxylated almost exclusively at the secondary ethyl carbon with about a 2:1 ratio of R:S product. Several molecular dynamics trajectories were performed with different starting conformations of ethylbenzene in the active site of P450cam. The stereochemistry of hydroxylation predicted from the molecular dynamics simulations was found to be in good agreement with the observed products.     

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.     

Enantiomeric discrimination of Ru-substrates by cytochrome P450cam

Molecules with photosensitizers attached to substrates (Wilker et al., Angew. Chem. Int. Ed. 38 (1999) 90-92) or cofactors (Hamachi et al., J. Am. Chem. Soc. 121 (1999) 5500-5506) can rapidly deliver redox equivalents to buried active sites. The structure of cytochrome P450cam (P450) co-crystallized with a prototypal sensitizer-substrate, [Ru-C9-Ad]Cl2, has been determined (Dmochowski et al., Proc. Natl. Acad. Sci. USA 96 (1999) 12987-12990); and, in separate UV-vis absorption and time-resolved luminescence experiments, the binding of the lambda and delta enantiomers of Ru-C9-Ad to P450 has been measured. The results, KD(delta/lambda) approximately 2, indicate that the bipyridyl ligands of the lambda isomer interact more favorably with hydrophobic residues at the entrance to the substrate channel. We conclude that enantiospecific interactions may be exploited in the design of enzyme-metallosubstrate conjugates.     

Superoxide anion production by the autoxidation of cytochrome P450cam

Chemiluminescence occurs on autoxidation of oxygenated ferrous cytochrome P450_cam and is abolished by reagents that scavenge free radicals, by superoxide dismutase and singlet oxygen quenchers. We attribute the chemiluminescence to the decay of an excited singlet oxygen which arises from a superoxide anion radical precursor.     

Identification of selectivity-determining residues in cytochrome P450 monooxygenases: a systematic analysis of the substrate recognition site 5

The large and diverse family of cytochrome P450 monooxygenases (CYPs) was systematically analyzed to identify selectivity- and specificity-determining residues in the substrate recognition site 5, which is located in close vicinity to the heme center. A positively charged heme-interacting residue was identified in the structures of 29 monooxygenases and in 97.7% of the 6379 CYP sequences investigated here. This heme-interacting residue restricts the conformation of the substrate recognition site 5 and is preferentially located at position 10 or 11 after the conserved ExxR motif (in 94.4% of the sequences), in 3.3% of the sequences at position 9 or 12. As a result, a classification by the position of the heme-interacting residue allows to predict residues that are closest to the heme center and restrict its accessibility. In 98.4% of all CYP sequences a preferentially hydrophobic residue is located at position 5 after the ExxR motif that is predicted to point close to the heme center. Replacing this residue by hydrophobic residues of different size has been shown to change substrate specificity and regioselectivity for CYPs of different superfamilies. Twenty-seven percent of all CYPs are predicted to contain a second selectivity-determining residue at position 9 after the ExxR motif that can be identified by the pattern EXXR-X(7)-{P}-x-P-[HKR].     

Spring-loading the active site of cytochrome P450cam

A hydrogen bond network has been identified that adjusts protein-substrate contacts in cytochrome P450(cam) (CYP101A1). Replacing the native substrate camphor with adamantanone or norcamphor causes perturbations in NMR-detected NH correlations assigned to the network, which includes portions of a β sheet and an adjacent helix that is remote from the active site. A mutation in this helix reduces enzyme efficiency and perturbs the extent of substrate-induced spin state changes at the haem iron that accompany substrate binding. In turn, the magnitude of the spin state changes induced by alternate substrate binding parallel the NMR-detected perturbations observed near the haem in the enzyme active site.     

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.     

A functional proline switch in cytochrome P450cam

The two-protein complex between putidaredoxin (Pdx) and cytochrome P450(cam) (CYP101) is the catalytically competent species for camphor hydroxylation by CYP101. We detected a conformational change in CYP101 upon binding of Pdx that reorients bound camphor appropriately for hydroxylation. Experimental evidence shows that binding of Pdx converts a single X-proline amide bond in CYP101 from trans or distorted trans to cis. Mutation of proline 89 to isoleucine yields a mixture of both bound camphor orientations, that seen in Pdx-free and that seen in Pdx-bound CYP101. A mutation in CYP101 that destabilizes the cis conformer of the Ile 88-Pro 89 amide bond results in weaker binding of Pdx. This work provides direct experimental evidence for involvement of X-proline isomerization in enzyme function.     

Electron transfer reactions in Zn-substituted cytochrome P450cam

We have investigated photoinduced electron transfer (ET) reactions between zinc-substituted cytochrome P450cam (ZnP450) and several inorganic reagents by using the laser flash photolysis method, to reveal roles of the electrostatic interactions in the regulation of the ET reactions. The laser pulse irradiation to ZnP450 yielded a strong reductant, the triplet excited state of ZnP450, (3)ZnP450, which was able to transfer one electron to anionic redox partners, OsCl(6)(2-) and Fe(CN)(6)(3-), with formation of the porphyrin pi-cation radical, ZnP450(+). In contrast, the ET reactions from (3)ZnP450 to cationic redox partners, such as Ru(NH(3))(6)(3+) and Co(phen)(3)(3+), were not observed even in the presence of 100-fold excess of the oxidant. One of the possible interpretations for the preferential ET to the anionic redox partner is that the cationic patch on the P450cam surface, a putative interaction site for the anionic reagents, is located near the heme (less than 10 A from the heme edge), while the anionic surface is far from the heme moiety (more than 16 A from the heme edge), which would yield 8000-fold faster ET rates through the cationic patch. The ET rate through the anionic patch to the cationic partner would be substantially slower than that of the phosphorescence process in (3)ZnP450, resulting in no ET reactions to the cationic reagents. These results demonstrate that the asymmetrical charge distribution on the protein surface is critical for the ET reaction in P450cam.     

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.     

Linking of cytochrome P450cam and putidaredoxin by a co-ordination bridge

Cytochrome P450_cam (CYP101) catalyzes the oxidation of D(+)-camphor at the 5 position. The enzyme couples the reduction of dioxygen to the oxidation of the substrate. To transfer electrons from the reductant (NADH) to the cytochrome, two additional proteins are required. These are putidaredoxin (PdX) and putidaredoxin reductase (PdR). We have chemically linked a form of PdX with a histidine tag at the C-terminus to the P450. To accomplish this, we have modified cysteine 334 on P450 with a bipyridinyl group, and co-ordinated the C-terminal histidine tag of PdX by the addition of Ni²⁺ or Ru³⁺. The Ru³⁺ complex was the most stable. The non-linked system gave mostly 5-ketocamphor, a product of two consecutive hydroxylations, and H₂O₂, a product of 2-electron uncoupling. The Ni²⁺ complex gave both 5-exo-hydroxycamphor and 5-ketocamphor, but it also uncoupled. The Ru³⁺ complex gave a single product (5-exo-hydroxycamphor) and did not uncouple at the optimal PdR concentration. Our results are consistent with other studies of this system, in that strong binding of PdX to P450 is crucial for good coupling and for release of 5-exo-hydroxycamphor.     

Model of the coordination site for cobalt-substituted cytochrome P450cam

Cobalt-substituted cytochrome P450_cam was recently reconstituted by Wagner et al (J. Biol. Chem. 256, 6266, (1981)). A model of its coordination site was constructed to determine the mode of axial coordination of the native enzyme. Complexes were prepared from cobalt porphyrins (cobalt-protoporphyrin IX (CoPPIX), cobalt-meso-tetraphenylporphyrin, cobalt-γ-laurylpyridyl triphenylporphyrin, and cobalt-octaethylporphyrin), thioglycolate ester, and tetramethylammonium hydroxide in organic solvents. Complexes prepared in an organic solvent such as CHC1₃ under air at room temperature exhibited a stable Soret hyperporphyrin spectrum characterized by split Soret bands, especially like that of the thiol-Co-P450_cam complex Comparison of the spectra of the hyperporphyrin spectral complexes titrated with various types of alcohol and imidazole, with the spectrum of Co-P450_cam in the oxidized state support the idea that an axial thiolate at the fifth position and a hydroxyl group of alcohol at the sixth position of the heme form the coordination site of Co-P450_cam The CoPPIX-thiolate-ethanol complex retaining S⁻-Co(III)-OH coordination is thought to be a possible model of Co-P450_cam in the oxidized state.     

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.     

Differential contribution of active site residues in substrate recognition sites 1 and 5 to cytochrome P450 2C8 substrate selectivity and regioselectivity

Selected active site residues in substrate recognition sites (SRS) 1 and 5 of cytochrome P450 2C8 (CYP2C8) were mutated to the corresponding amino acids present in CYP2C9 to investigate the contribution of these positions to the unique substrate selectivity and regioselectivity of CYP2C8. The effects of mutations, singly and in combination, were assessed from changes in the kinetics of paclitaxel 6alpha-hydroxylation, a CYP2C8-specific pathway, and the tolylmethyl and ring hydroxylations of torsemide, a mixed CYP2C9/CYP2C8 substrate. Within SRS1, the single mutation S114F abolished paclitaxel 6alpha-hydroxylation, while the I113V substitution resulted in modest parallel reductions in K(m) and V(max). Mutations in SRS5 (viz., V362L, G365S, and V366L) reduced paclitaxel intrinsic clearance (V(max)/K(m)) by 88-100%. Torsemide is preferentially metabolized by CYP2C9, and it was anticipated that the mutations in CYP2C8 might increase activity. However, methyl and ring hydroxylation intrinsic clearances were either unchanged or decreased by the mutations, although hydroxylation regioselectivity was often altered relative to wild-type CYP2C8. The mutations significantly increased (28-968%) K(m) values for both torsemide methyl and ring hydroxylation but had variable effects on V(max). The effects of the combined mutations in SRS1, SRS5, and SRS1 plus SRS5 were generally consistent with the changes produced by the separate mutations. Mutation of CYP2C8 at position 359 (S359I), a site of genetic polymorphism in CYP2C9, resulted in relatively minor changes in paclitaxel- and torsemide-hydroxylase activities. The results are consistent with multiple substrate binding orientations within the CYP2C8 active site and a differential contribution of active site residues to paclitaxel and torsemide binding and turnover.     

Three clusters of conformational states in p450cam reveal a multistep pathway for closing of the substrate access channel

Conformational changes in the substrate access channel have been observed for several forms of cytochrome P450, but the extent of conformational plasticity exhibited by a given isozyme has not been completely characterized. Here we present crystal structures of P450cam bound to a library of 12 active site probes containing a substrate analogue tethered to a variable linker. The structures provide a unique view of the range of protein conformations accessible during substrate binding. Principal component analysis of a total of 30 structures reveals three discrete clusters of conformations: closed (P450cam-C), intermediate (P450cam-I), and fully open (P450cam-O). Relative to P450cam-C, the P450cam-I state results predominantly from a retraction of helix F, while both helices F and G move in concert to reach the fully open P450cam-O state. Both P450cam-C and P450cam-I are well-defined states, while P450cam-O shows evidence of a somewhat broader distribution of conformations and includes the open form recently seen in the absence of substrate. The observed clustering of protein conformations over a wide range of ligand variants suggests a multistep closure of the enzyme around the substrate that begins by conformational selection from an ensemble of open conformations and proceeds through a well-defined intermediate, P450cam-I, before full closure to the P450cam-C state in the presence of small substrates. This multistep pathway may have significant implications for a full understanding of substrate specificity, kinetics, and coupling of substrate binding to P450 function.     

Uncovering the role of hydrophobic residues in cytochrome P450-cytochrome P450 reductase interactions

Cytochrome P450 (CYP or P450)-mediated drug metabolism requires the interaction of P450s with their redox partner, cytochrome P450 reductase (CPR). In this work, we have investigated the role of P450 hydrophobic residues in complex formation with CPR and uncovered novel roles for the surface-exposed residues V267 and L270 of CYP2B4 in mediating CYP2B4--CPR interactions. Using a combination of fluorescence labeling and stopped-flow spectroscopy, we have investigated the basis for these interactions. Specifically, in order to study P450--CPR interactions, a single reactive cysteine was introduced in to a genetically engineered variant of CYP2B4 (C79SC152S) at each of seven strategically selected surface-exposed positions. Each of these cysteine residues was modified by reaction with fluorescein-5-maleimide (FM), and the CYP2B4-FM variants were then used to determine the K(d) of the complex by monitoring fluorescence enhancement in the presence of CPR. Furthermore, the intrinsic K(m) values of the CYP2B4 variants for CPR were measured, and stopped-flow spectroscopy was used to determine the intrinsic kinetics and the extent of reduction of the ferric P450 mutants to the ferrous P450--CO adduct by CPR. A comparison of the results from these three approaches reveals that the sites on P450 exhibiting the greatest changes in fluorescence intensity upon binding CPR are associated with the greatest increases in the K(m) values of the P450 variants for CPR and with the greatest decreases in the rates and extents of reduced P450--CO formation.     

The role of cytochrome P450 lysine residues in the interaction between cytochrome P450IA1 and NADPH-cytochrome P450 reductase.

Cytochrome P450IA1 (purified from hepatic microsomes of beta-naphthoflavone-treated rats) has been covalently modified with the lysine-modifying reagent acetic anhydride. Different levels of lysine residue modification in cytochrome P450IA1 can be achieved by varying the concentration of acetic anhydride. Modification of lysine residues in P450IA1 greatly inhibits the interaction of P450IA1 with NADPH-cytochrome P450 reductase. Modification of 1.0 and 3.3 mol lysine residues per mole P450IA1 resulted in 30 and 95% decreases, respectively, in 7-ethoxycoumarin hydroxylation by a reconstituted P450IA1/reductase complex. However, modification of 3.3 mol lysine residues per mole P450IA1 decreased only cumene hydroperoxide-supported P450-dependent 7-ethoxycoumarin hydroxylation by 30%. Spectral and fluorescence studies showed no indication of global conformational change of P450IA1 even with up to 8.8 mol lysine residues modified per mole P450IA1. These data suggest that at least three lysine residues in P450IA1 may be involved in the interaction with reductase. Identification of lysine residues in P450IA1 possibly involved in this interaction was carried out by [14C]acetic anhydride modification, trypsin digestion, HPLC separation, and amino acid sequencing. The lysine residue candidates identified in this manner were K97, K271, K279, and K407.