A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment

Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 μs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

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.

Stereo-preference of camphor for H-bonding with phenol, methanol and chloroform: A combined matrix isolation IR spectroscopic and quantum chemical investigation

Camphor is known to be held in the substrate pocket of cytochrome P450cam enzyme via H-bond with a tyrosine residue of the enzyme in a unique orientation. This structural exclusivity results in regio- and stereo-specific hydroxylation of camphor by the enzyme. We have carried out a combined IR spectroscopic and quantum chemical investigation to shed light on the factors influencing the conformational exclusivity of 1R-(+)-camphor in the substrate pocket of Cytochrome P450cam, and to determine whether the selectivity is an inherent property of the substrate itself, or is imposed by the enzyme. For this purpose, complexes of camphor have been studied with three H-bond donors namely phenol, methanol and chloroform. Each of the three donors was found to form stable complexes with two distinct conformers; the one mimicking the conformation in enzyme substrate pocket was found to be more stable of the two, for all three donors. Experimentally, both conformers of the H-bonded complexes were identified separately for phenol and methanol in an argon matrix at 8 K, but not for chloroform due to very small energy barrier for interconversion of the two conformers. In room temperature solution phase spectra of camphor with all three donors, the differences in spectral attributes between the two isomeric H-bonded complexes were lost due to thermal motions.

Effect of redox partner binding on CYP101D1 conformational dynamics

We have compared the thermodynamics of substrate and redox partner binding of P450cam to its close homologue, CYP101D1, using isothermal titration calorimetry (ITC). CYP101D1 binds camphor about 10-fold more weakly than P450cam which is consistent with the inability of camphor to cause a complete low- to high-spin shift in CYP101D1. Even so molecular dynamics simulations show that camphor is very stable in the CYP101D1 active site similar to P450cam. ITC data on the binding of the CYP101D1 ferredoxin redox partner (abbreviated Arx) shows that the substrate-bound closed state of CYP101D1 binds Arx more tightly than the substrate-free open form. This is just the opposite to P450cam where Pdx (ferredoxin redox partner of P450cam) favors binding to the P450cam open state. In addition, CYP101D1-Arx binding has a large negative ΔS while the P450cam-Pdx has a much smaller ΔS indicating that interactions at the docking interface are different. The most obvious difference is that PDXD38 which forms an important ion pair with P450camR112 at the center of the interface is ArxL39 in Arx. This suggests that Arx may adopt a different orientation than Pdx in order to optimize nonpolar interactions with ArxL39.

Chemoenzymatic Synthesis of a Novel Borneol-Based Polyester

Terpenes are a class of natural compounds that have recently moved into the focus as a bio-based resource for chemical production, owing to their abundance, their mostly cyclic structures, and the presence of olefin or single hydroxy groups. To apply this raw material in new industrial fields, a second hydroxy group is inserted into borneol by cytochrome P450cam (CYP101) enzymes in a whole-cell catalytic biotransformation with Pseudomonas putida KT2440. Next, a semi-continuous batch system was developed to produce 5-exo-hydroxyborneol with a final concentration of 0.54 g L^-1 . The bifunctional terpene was then used for the synthesis of a bio-based polyester by a solvent-free polycondensation reaction. The resulting polymer showed a glass transition temperature of around 70 °C and a molecular weight in the range of 2000-4000 g mol^-1 (M_w ). These results show that whole-cell catalytic biotransformation of terpenes could lead to bio-based, higher-functionalized monomers, which might be basic raw materials for different fields of application, such as biopolymers.

Detection of substrate-dependent conformational changes in the P450 fold by nuclear magnetic resonance

Cytochrome P450 monooxygenases typically catalyze the insertion of one atom of oxygen from O2 into unactivated carbon-hydrogen and carbon-carbon bonds, with concomitant reduction of the other oxygen atom to H2O by NAD(P)H. Comparison of the average structures of the camphor hydroxylase cytochrome P450(cam) (CYP101) obtained from residual dipolar coupling (RDC)-restrained molecular dynamics (MD) in the presence and absence of substrate camphor shows structural displacements resulting from the essential collapse of the active site upon substrate removal. This collapse has conformational consequences that extend across the protein structure, none of which were observed in analogous crystallographic structures. Mutations were made to test the involvement of the observed conformational changes in substrate binding and recognition. All of the mutations performed based upon the NMR-detected perturbations, even those remote from the active site, resulted in modified substrate selectivity, enzyme efficiency and/or haem iron spin state. The results demonstrate that solution NMR can provide insights into enzyme structure-function relationships that are difficult to obtain by other methods.

Synergistic effects of mutations in cytochrome P450cam designed to mimic CYP101D1

A close orthologue to cytochrome P450cam (CYP101A1) that catalyzes the same hydroxylation of camphor to 5-exo-hydroxycamphor is CYP101D1. There are potentially important differences in and around the active site that could contribute to subtle functional differences. Adjacent to the heme iron ligand, Cys357, is Leu358 in P450cam, whereas this residue is Ala in CYP101D1. Leu358 plays a role in binding of the P450cam redox partner, putidaredoxin (Pdx). On the opposite side of the heme, about 15-20 Å away, Asp251 in P450cam plays a critical role in a proton relay network required for O2 activation but forms strong ion pairs with Arg186 and Lys178. In CYP101D1 Gly replaces Lys178. Thus, the local electrostatic environment and ion pairing are substantially different in CYP101D1. These sites have been systematically mutated in P450cam to the corresponding residues in CYP101D1 and the mutants analyzed by crystallography, kinetics, and UV-vis spectroscopy. Individually, the mutants have little effect on activity or structure, but in combination there is a major drop in enzyme activity. This loss in activity is due to the mutants being locked in the low-spin state, which prevents electron transfer from the P450cam redox partner, Pdx. These studies illustrate the strong synergistic effects on well-separated parts of the structure in controlling the equilibrium between the open (low-spin) and closed (high-spin) conformational states.

Effect of length of molecular recognition moiety on enzymatic activity switching

We site-specifically conjugated biotin-PEG derivatives with spacer arms of different lengths to mutant P450cam (3mD) and evaluated the activity of and structural changes in the conjugates as a first step toward clarifying the mechanism whereby the activity of the 3mD conjugate is inhibited. 3mD was prepared by site-specific mutation to inhibit its enzymatic activity artificially, after which the derivative compounds were conjugated to the enzyme. 3mD has one cysteine on its surface with a reactive thiol group that can react with compounds near the active site, where a conformational change will be induced after conjugation. The activity of 3mD was retained in the biotin-PEG₂-3mD conjugate, but was dramatically reduced in the biotin-PEG₁₁-3mD conjugate. To investigate the effect of poly(ethylene glycol) (PEG) length on the enzymatic activity after conjugation, PEGs of different lengths, exceeding that in biotin-PEG₁₁, and whose termini were not biotin, were conjugated to 3mD. The activity of 3mD decreased in all these conjugates. This indicates that the activity of 3mD in these conjugates decreased after its conjugation with PEG molecules that exceeded a certain length. The biotin-PEG₂-3mD, which retains enzymatic activity after conjugation, showed avidin responsiveness; the enzymatic activity decreased after avidin binding.

Structural and dynamic implications of an effector-induced backbone amide cis-trans isomerization in cytochrome P450cam

Experimental evidence has been provided for a functionally relevant cis-trans isomerization of the Ile88-Pro89 peptide bond in cytochrome P450(cam) (CYP101). The isomerization is proposed to be a key element of the structural reorganization leading to the catalytically competent form of CYP101 upon binding of the effector protein putidaredoxin (Pdx). A detailed comparison of the results of molecular dynamics simulations on the cis and trans conformations of substrate- and carbonmonoxy-bound ferrous CYP101 with sequence-specific Pdx-induced structural perturbations identified by nuclear magnetic resonance is presented, providing insight into the structural and dynamic consequences of the isomerization. The mechanical coupling between the Pdx binding site on the proximal face of CYP101 and the site of isomerization is described.

13C-Methyl isocyanide as an NMR probe for cytochrome P450 active sites

The cytochromes P450 (CYPs) play a central role in many biologically important oxidation reactions, including the metabolism of drugs and other xenobiotic compounds. Because they are often assayed as both drug targets and anti-targets, any tools that provide: (a) confirmation of active site binding and (b) structural data, would be of great utility, especially if data could be obtained in reasonably high throughput. To this end, we have developed an analog of the promiscuous heme ligand, cyanide, with a (13)CH(3)-reporter attached. This (13)C-methyl isocyanide ligand binds to bacterial (P450cam) and membrane-bound mammalian (CYP2B4) CYPs. It can be used in a rapid 1D experiment to identify binders, and provides a qualitative measure of structural changes in the active site.

Substrate binding induces structural changes in cytochrome P450cam

The binding of (+)-camphor to cytochrome P450cam (P450cam) expels a cluster of waters at the active site, raising the redox potential of the haem to an extent that allows reduction by the electron-transfer system. This binding was reported to involve no significant structural changes in the protein. Here, two ferric P450cam structures partially complexed with (+)-camphor were determined by X-ray crystallography at 1.30-1.35 A resolution, revealing the structures of the substrate-free and substrate-bound forms. (+)-Camphor binding induces rotation of Thr101 to form a hydrogen bond that acts as a hydrogen donor to a peripheral haem propionate. This bonding contributes to the redox-potential change.

Characterizing metabolic inhibition using electrochemical enzyme/DNA biosensors

Studies of metabolic enzyme inhibition are necessary in drug development and toxicity investigations as potential tools to limit or prevent appearance of deleterious metabolites formed, for example, by cytochrome (cyt) P450 enzymes. In this paper, we evaluate the use of enzyme/DNA toxicity biosensors as tools to investigate enzyme inhibition. We have examined DNA damage due to cyt P450cam metabolism of styrene using DNA/enzyme films on pyrolytic graphite (PG) electrodes monitored via Ru(bpy)(3)(2+)-mediated DNA oxidation. Styrene metabolism initiated by hydrogen peroxide was evaluated with and without the inhibitors, imidazole, imidazole-4-acetic acid, and sulconazole (in micromolar range) to monitor DNA damage inhibition. The initial rates of DNA damage decreased with increased inhibitor concentrations. Linear and nonlinear fits of Michaelis-Menten inhibition models were used to determine apparent inhibition constants (K(I)*) for the inhibitors. Elucidation of the best fitting inhibition model was achieved by comparing correlation coefficients and the sum of the square of the errors (SSE) from each inhibition model. Results confirmed the utility of the enzyme/DNA biosensor for metabolic inhibition studies. A simple competitive inhibition model best approximated the data for imidazole, imidazole-4-acetic acid and sulconazole with K(I)* of 268.2, 142.3, and 204.2 microM, respectively.

Electrochemiluminescent/voltammetric toxicity screening sensor using enzyme-generated DNA damage

Simultaneous optical and voltammetric detection of bioactivated genotoxicity is reported for the first time employing ultrathin films of DNA, model metabolic enzymes, and electrochemiluminescence (ECL) generating metallopolymer [Ru(bpy)2PVP10]2+ on pyrolytic graphite (PG) electrodes. Cytochrome P450cam and myoglobin were used as model monoxygenase enzymes to mimic in vivo processes. Sensor film growth and component amounts were monitored using a quartz crystal microbalance (QCM). Subsequent to the enzyme reaction, DNA damage in the sensor films was measured simultaneously using a simple apparatus combining a standard voltammetry cell coupled with an optical fiber and photomultiplier tube. The model enzyme reaction converted styrene to styrene oxide, which reacts with DNA nucleobases. ECL and SWV signals increased with enzyme reaction time on the scale of several min, and provided relative enzyme turnover rates for DNA damage suitable for toxicity screening applications. Within 1 min, the sensor detects approximately 3 damaged bases per 10,000 DNA bases using this simultaneous detection.

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.

Intramolecular electron transfer in a cytochrome P450cam system with a site-specific branched structure

Cytochrome P450 (P450) is an attractive oxygenase due to the diverse catalytic reactions and the broad substrate specificity. Class I P450s require an excess concentration (more than 10 times) of iron-sulfur proteins, which transfer electrons to P450s, to attain the maximum catalytic activity and this requirement is a critical bottleneck for practical applications. Here, we show a site-specific branched fusion protein of P450 with its electron transfer proteins using enzymatic cross-linking with transglutaminase. A branched fusion protein of P450 from Pseudomonas putida (P450cam), which was composed of one molecule each of P450cam, putidaredoxin (Pdx) and Pdx reductase, showed higher catalytic activity (306 min(-1)) and coupling efficiency (99%) than the equimolar reconstitution system due to the intramolecular electron transfer. The unique site-specific branched structure simply increased local concentration of proteins without denaturation of each protein. Therefore, enzymatic post-translational protein manipulation can be a powerful alternative to conventional strategies for the creation of multicomponent enzyme systems with novel proteinaceous architecture.

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.

FTIR studies of the redox partner interaction in cytochrome P450: The Pdx–P450cam couple

Recently we have developed a new approach to study protein-protein interactions using Fourier transform infrared spectroscopy in combination with titration experiments and principal component analysis (FTIR-TPCA). In the present paper we review the FTIR-TPCA results obtained for the interaction between cytochrome P450 and the redox partner protein in two P450 systems, the Pseudomonas putida P450cam (CYP101) with putidaredoxin (P450cam-Pdx), and the Bacillus megaterium P450BM-3 (CYP102) heme domain with the FMN domain (P450BMP-FMND). Both P450 systems reveal similarities in the structural changes that occur upon redox partner complex formation. These involve an increase in beta-sheets and alpha-helix content, a decrease in the population of random coil/3(10)-helix structure, a redistribution of turn structures within the interacting proteins and changes in the protonation states or hydrogen-bonding of amino acid carboxylic side chains. We discuss in detail the P450cam-Pdx interaction in comparison with literature data and conclusions drawn from experiments obtained by other spectroscopic techniques. The results are also interpreted in the context of a 3D structural model of the Pdx-P450cam complex.

Specific effects of potassium ion binding on wild-type and L358P cytochrome P450cam

The camphor monoxygenase cytochrome P450cam (CYP101) requires potassium ion (K+) to drive formation of the characteristic high-spin state of the heme Fe+3 upon substrate binding. Amide 1H, 15N correlations in perdeuterated [U-15N] CYP101 were monitored as a function of K+ concentration by 2D-TROSY-HSQC in both camphor-bound oxidized (CYP-S) and camphor- and CO-bound reduced CYP101 (CYP-S-CO). In both forms, K+-induced spectral perturbations are detected in the vicinity of the K+ binding site proposed from crystallographic structures, but are larger and more widespread structurally in CYP-S than in CYP-S-CO. In CYP-S-CO, K+-induced perturbations occur primarily near the proposed K+ binding site in the B-B' loop and B' helix, which are also perturbed by binding of effector, putidaredoxin (Pdx). The spectral effects of K+ binding in CYP-S-CO oppose those observed upon Pdxr titration. However, Pdxr titration of CYP-S-CO in the absence of K+ results in multiple conformations. The spin-state equilibrium in the L358P mutant of CYP101 is more sensitive to K+ concentration than WT CYP101, consistent with a hypothesis that L358P preferentially populates conformations enforced by Pdx binding in WT CYP101. Thallium(I), a K+ mimic, minimizes the effects of Pdx titration on the NMR spectrum of CYP-S-CO, but is competent to replace K+ in driving the formation of high-spin CYP-S. These observations suggest that the role of K+ is to stabilize conformers of CYP-S that drive the spin-state change prior to the first electron transfer, and that K+ stabilizes the CYP-S-CO conformer that interacts with Pdx. However, upon binding of Pdx, further conformational changes occur that disfavor K+ binding.

Mechanism of O2Activation by Cytochrome P450cam Studied by Isotope Effects and Transient State Kinetics†

The early steps in dioxygen activation by the monooxygenase cytochrome P450cam (CYP101) include binding of O2 to ferrous P450cam to yield the ferric-superoxo form (oxyP450cam) followed by an irreversible, long-range electron transfer from putidaredoxin to reduce the oxyP450cam. The steady state kinetic parameter kcat/Km(O2) has been studied by a variety of probes that indicate a small D2O solvent isotope effect (1.21 +/- 0.08), a very small solvent viscosogen effect, and a 16O/18O isotope effect of 1.0147 +/- 0.0007. This latter value, which can be compared with the 16O/18O equilibrium isotope effect of 1.0048 +/- 0.0003 measured for oxyP450cam formation, is attributed to a primarily rate-limiting outer-sphere electron transfer from the heme iron center as O2 that has prebound to protein approaches the active site cofactor. The electron transfer from putidaredoxin to oxyP450cam was investigated by rapid mixing at 25 degrees C to complement previous lower-temperature measurements. A rate of 390 +/- 23 s-1 (and a near-unity solvent isotope effect) supports the view that the long-range electron transfer from reduced putidaredoxin to oxyP450cam is rapid relative to dissociation of O2 from the enzyme. P450cam represents the first enzymatic reaction of O2 in which both equilibrium and kinetic 16O/18O isotope effects have been measured.

On the identity and reactivity patterns of the “second oxidant” of the T252A mutant of cytochrome P450cam in the oxidation of 5-methylenenylcamphor

Density functional calculations show that in the absence of Compound I, the primary oxidant species of P450, the precursor species, Compound 0 (FeOOH), can effect double bond activation of 5-methylenylcamphor (1). The mechanism is initiated by homolytic cleavage of the O-O bond and formation of an OH radical bound to the Compound II species by hydrogen bonding interactions. Subsequently, the so-formed OH radical can either activate the double bond of 1 or attack the meso position of the heme en route to heme degradation. The calculations show that double bond activation is preferred over attack on the heme. Past the double bond activation, the intermediate can either lead to epoxidation or to a glycol formation. The glycol formation is predicted to be preferred, although in the P450(cam) pocket the competition may be closer. Therefore, in the absence of Compound I, Compound 0 will be capable of epoxidizing double bonds. Previous studies [E. Derat, D. Kumar, H. Hirao, S. Shaik, J. Am. Chem. Soc. 128 (2006) 473-484] showed that in the case of a substrate that can undergo only C-H activation, the bound OH prefers heme degradation over hydrogen abstraction. Since the epoxidation barrier for Compound I is much smaller than that of Compound 0 (12.8 vs. 18.9kcal/mol), when Compound I is present in the cycle, Compound 0 will be silent. As such, our mechanism explains lucidly why T252A P450(cam) can epoxidize olefins like 5-methylenylcamphor but is ineffective in camphor hydroxylation [S. Jin, T.M. Makris, T. A. Bryson, S.G. Sligar, J.H. Dawson, J. Am. Chem. Soc. 125 (2003) 3406-3407]. Our calculations show that the glycol formation is a marker reaction of Compound 0 with 5-methylenylcamphor. If this product can be found in T252A P450(cam) or in similar mutants of other P450 isozymes, this will constitute a more definitive proof for the action of Cpd 0 in P450 enzymes.

Rapid kinetics investigations of peracid oxidation of ferric cytochrome P450cam: nature and possible function of compound ES

Previously, we reported spectroscopic properties of cytochrome P450cam compound I, (ferryl iron plus a porphyrin pi-cation radical (Fe(IV)=O/Por(+))), as well as compound ES (Fe(IV)=O/Tyr()) in reactions of substrate-free ferric enzyme with m-chloroperbenzoic acid [T. Spolitak, J.H. Dawson, D.P. Ballou, J. Biol. Chem. 280 (2005) 20300-9]. Compound ES arises by intramolecular electron transfer from nearby tyrosines to the porphyrin pi-cation radical of Compound I, and has been characterized by rapid-freeze-quench-Mössbauer/EPR spectroscopy; the tyrosyl radical was assigned to Tyr96 for wild type or to Tyr75 for the Tyr96Phe variant [V. Schünemann, F. Lendzian, C. Jung, J. Contzen, A.L. Barra, S.G. Sligar, A.X. Trautwein, J. Biol. Chem. 279 (2004) 10919-10930]. Here we report rapid-scanning stopped-flow studies of the reactions of peracids with substrate-free ferric Y75F, Y96F, and Y96F/Y75F P450cam variants, showing how these active site changes influence electron transfer from nearby tyrosines and affect formation of intermediates. Curiously, rates of generation of Compounds I and ES for both single mutants were not very different from wild type. Contrasting with the earlier EPR results, the Y96F/Y75F variant was also shown to form an ES-like species, but more slowly. When substrate is not present, or is improperly bound, compound I rapidly converts to compound ES, which can be reduced to form H(2)O and ferric P450, thus avoiding the modification of nearby protein groups or release of reactive oxygen species.

Differential behavior of the sub-sites of cytochrome 450 active site in binding of substrates, and products (implications for coupling/uncoupling)

The cytochrome P450 catalyzes hydroxylation of many substrates in the presence of O(2) and specific electron transport system. The ternary complex S-Fe(+)O(2) with substrate and O(2) bound to their respective sites on the reduced enzyme is an important intermediate in the formation of the hydroxylating species. Then the active site may be considered as having two sub-sites geared for entirely different types of functionally relevant interactions. The two sites are the substrate binding site, the specific protein residues (Site I), and the L(6) position of the iron (Site II) to which O(2) binds upon reduction. In the ferric enzyme, when substrate binds to Site I, the low spin six-coordinated P450 is converted to the readily reducible high spin five coordinated state. Certain amines and OH compounds, such as products of P450-catalyzed reactions, can bind to Site II resulting in six coordinated inhibited complexes. Then the substrate and product interactions with the two sub-sites can regulate the functional state of the enzyme during catalysis. Product interactions have received very little attention. CYP101 is the only P450 in which X-ray and spectroscopic data on all three structures, the substrate-free, camphor-bound and the 5-exo-OHcamphor-bound are available. The substrate-free CYP101 is low spin and six-coordinated with a water molecule ligated at the L(6) position of the iron. The substrate camphor binds to Site I, and releases the L(6) water despite its inability to bind to this site, indicating that Site I binding can inhibit Site II ligation. The product 5-exo-OHcamphor in addition to binding to Site I, binds to Site II through its -OH group forming Fe-O bond, resulting in the low spin six-coordinated complex. New temperature-jump relaxation kinetic data indicating that Site II ligation inhibits Site I binding are presented. It appears that the Site I and Site II function as interacting sub-sites. The inhibitory allosteric interactions between the two sub-sites are also reflected in the data on binding of the substrate camphor (S) in the presence of the product 5-exo-OH camphor (P) to CYP101 (E). The data are in accordance with the two-site model involving the ternary complex ESP. The affinity of the substrate to the product-bound enzyme as well as the affinity of the product to the substrate-bound enzyme decreased with increase in product concentration, which is consistent with mixed inhibition indicative of inhibitory allosteric interactions between the two sub-sites. Implications of these observations for coupling/uncoupling mechanisms are discussed in the light of the published findings consistent with the two-site behavior of the P450 active site. In addition, kinetic data indicating that the transient high spin intermediate may have to be taken into account for understanding how some P450s have been able to express appreciable hydroxylation activities in the absence of substrate-induced low to high spin transition, observable by the traditional static spectroscopy, are presented.

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.

QM/MM Study of Mechanisms for Compound I Formation in the Catalytic Cycle of Cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.

Cytochrome P450 active site plasticity: attenuation of imidazole binding in cytochrome P450cam by an L244A mutation

We have identified a P450(cam) mutation, L244A, that mitigates the affinity for imidazole and substituted imidazoles while maintaining a high affinity for the natural substrate camphor. The P450(cam) L244A crystal structure solved in the absence of any ligand reveals that the I-helix is displaced inwards by over 1 A in response to the cavity created by the change from leucine to alanine. Furthermore, the crystal structures of imidazole-bound P450(cam) and the 1-methylimidazole-bound P450(cam) L244A mutant reveal that the ligands have distinct binding modes in the two proteins. Whereas in wild-type P450(cam) the imidazole coordinates to the iron in an orientation roughly perpendicular to the plane of the heme, in the L244A mutant the rearranged I helix, and specifically residue Val247, forces the imidazole into an orientation almost parallel to the heme that impairs its ability to coordinate to the heme iron. As a result, the imidazole is much more weakly bound to the mutant than it is to the wild-type enzyme. Despite the constriction of the active site by the mutation, previous work with the L244A mutant has shown that it oxidizes larger substrates than the wild-type enzyme. This paradoxical situation, in which a mutation that nominally increases the active site cavity appears to decrease it, suggests that the mutation actually increases the active site maleability, allowing it to better expand to oxidize larger substrates.

Increasing the Catalytic Performance of a Whole Cell Biocatalyst Harboring a Cytochrome P450cam System by Stabilization of an Electron Transfer Component

Catalytic activity of a recombinant Escherichia coli whole cell biocatalyst harboring a cytochrome P450cam monooxygenase system from Pseudomonas putida coupled with enzymatic co-factor regeneration was investigated. About 0.7 micromol camphor was hydroxylated per mg dry cells at 4 degrees C in 50 mM Tris/HCl buffer (pH 7.4) when utilizing a stable putidaredoxin (Pdx) mutant, C73S/C85S-Pdx (Cys73Ser, Cys85Ser double mutant), instead of wild-type Pdx, which was about two-fold improvement in the substrate conversion. Ten-micromole camphor was completely hydroxylated at 20 degrees C in 6 h by 15 mg dry cell weight of whole cell biocatalyst including C73S/C85S-Pdx. Thus, modulation of protein-protein interaction in multicomponent enzymatic catalysis in whole cells is important.

The effect of heme environment on the hydrogen abstraction reaction of camphor in P450cam catalysis: a QM/MM study

The discrepancies between the published QM/MM studies (Schöneboom, J. C.; Cohen, S.; Lin, H.; Shaik, S.; Thiel, W. J. Am. Chem. Soc. 2004, 126, 4017; Guallar, V.; Friesner, R. A. J. Am. Chem. Soc. 2004, 126, 8501) on H-abstraction of camphor in P450cam have largely been resolved. The crystallographic water molecule 903 situated near the oxo atom of Compound I acts as a catalyst for H-abstraction, lowering the barrier by about 4 kcal/mol. Spin density at the A-propionate side chain of heme can occur in the case of incomplete screening but has no major effect on the computed barrier.

Comparison of the complexes formed by cytochrome P450cam with cytochrome b5 and putidaredoxin, two effectors of camphor hydroxylase activity

Structural perturbations in cytochrome P450cam (CYP101) induced by the soluble fragment of cytochrome b5, a nonphysiological effector of CYP101, were investigated by NMR spectroscopy and compared with the perturbations induced by the physiological reductant and effector putidaredoxin (Pdx). Chemical shifts of perdeuterated [U-15N]CYP101 backbone amide (NH) resonances were monitored as a function of cytochrome b5 concentration by 1H-15N TROSY-HSQC experiments. The association of cytochrome b5 with the reduced CYP101-camphor-carbon monoxide complex (CYP-S-CO) perturbs many of the same resonances that Pdx does, including regions of the CYP101 molecule implicated in substrate access and orientation. The perturbations are smaller in magnitude than those observed with Pdx(r) due to a lower binding affinity (a Kd of 13 +/- 3 mM, for the reduced cytochrome b5-CYP-S-CO complex compared to a Kd of 26 +/- 12 microM for the Pdx-CYP-S-CO complex). The results are in accord with our previous suggestion that the observed perturbations are related to effector activity and support the proposal that the primary role of the effector is to populate the active conformation of CYP101 to prevent uncoupling [Pochapsky, S. S., et al. (2003) Biochemistry 42, 5649-5656]. A titratable perturbation is observed at the 1H resonance of the 8-CH3 group of CYP101-bound camphor upon addition of cytochrome b5, a phenomenon also associated with the formation of the CYP101 x Pdx complex, albeit with larger perturbations [Wei, J. Y., et al. (2005) J. Am. Chem. Soc. 127, 6974-6976]. The effector activity of the particular rat cytochrome b5 construct used for NMR studies was confirmed by monitoring the enzymatic turnover that yielded 5-exo-hydroxycamphor using gas chromatography and mass spectrometry. Finally, the common features of the perturbations observed in the NMR spectra of the two complexes are discussed, and their relevance to effector activity is considered.

Radical intermediates in the catalytic oxidation of hydrocarbons by bacterial and human cytochrome P450 enzymes

Cytochromes P450cam and P450BM3 oxidize alpha- and beta-thujone into multiple products, including 7-hydroxy-alpha-(or beta-)thujone, 7,8-dehydro-alpha-(or beta-)thujone, 4-hydroxy-alpha-(or beta-)thujone, 2-hydroxy-alpha-(or beta-)thujone, 5-hydroxy-5-isopropyl-2-methyl-2-cyclohexen-1-one, 4,10-dehydrothujone, and carvacrol. Quantitative analysis of the 4-hydroxylated isomers and the ring-opened product indicates that the hydroxylation proceeds via a radical mechanism with a radical recombination rate ranging from 0.7 +/- 0.3 x 10(10) s(-1) to 12.5 +/- 3 x 10(10) s(-1) for the trapping of the carbon radical by the iron-bound hydroxyl radical equivalent. 7-[2H]-alpha-Thujone has been synthesized and used to amplify C-4 hydroxylation in situations where uninformative C-7 hydroxylation is the dominant reaction. The involvement of a carbon radical intermediate is confirmed by the observation of inversion of stereochemistry of the methyl-substituted C-4 carbon during the hydroxylation. With an L244A mutation that slightly increases the P450(cam) active-site volume, this inversion is observed in up to 40% of the C-4 hydroxylated products. The oxidation of alpha-thujone by human CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4 occurs with up to 80% C-4 methyl inversion, in agreement with a dominant radical hydroxylation mechanism. Three minor desaturation products are produced, with at least one of them via a cationic pathway. The cation involved is proposed to form by electron abstraction from a radical intermediate. The absence of a solvent deuterium isotope effect on product distribution in the P450cam reaction precludes a significant role for the P450 ferric hydroperoxide intermediate in substrate hydroxylation. The results indicate that carbon hydroxylation is catalyzed exclusively by a P450 ferryl species via radical intermediates whose detailed properties are substrate- and enzyme-dependent.

A recombinant Escherichia coli whole cell biocatalyst harboring a cytochrome P450cam monooxygenase system coupled with enzymatic cofactor regeneration

A cytochrome P450cam monooxygenase (P450cam) system from the soil bacterium Pseudomonas putida requires electron transfer among three different proteins and a cofactor, nicotinamide adenine dinucleotide (NADH), for oxygenation of its natural substrate, camphor. Herein, we report a facile way to significantly enhance the catalytic efficiency of the P450cam system by the coupling of its native electron transfer system with enzymatic NADH regeneration catalyzed by glycerol dehydrogenase (GLD) in Escherichia coli whole cell biocatalysts. Recombinant E. coli harboring the P450cam system, but lacking GLD, exhibited little activity for camphor hydroxylation. In contrast, coexpression of GLD with the proteinaceous electron transfer components of P450cam resulted in about tenfold improvement in the substrate conversion, implying that the whole cell biocatalyst utilized molecular oxygen, endogenous NADH, and glycerol in the cell for catalysis. The addition of glycerol to the reaction media further promoted camphor hydroxylation, suggesting that exogenous glycerol is also available for GLD in the host cell and actively participates in the catalytic cycle. These results clearly show the utility of GLD towards functional reconstruction of the native P450cam system. The present approach may also be useful for E. coli whole cell biocatalysts with the other NADH-dependent oxygenases and oxidoreductases.

Effect of oxygen on the dehalogenation of 1,2-Dibromo-3-chloropropane by cytochrome P450cam (CYP101)

Cytochromes P450 are known to exhibit diverse catalytic functions against a large number of hydrocarbon substrates. The determinants of specific activity(ies) that operate on specific substrates have not been widely explored. Earlier, we showed that dehalogenation of 1,2-dibromo-3-chloropropane (DBCP) by P450cam (CYP101) monooxygenase exhibits oxygen- and substrate-dependent product distributions and reaction rates (1). Bromochloroacetone was the major conversion product when incubation media were saturated with oxygen, whereas allyl chloride was the sole product accounting for virtually all of the DBCP converted in the absence of oxygen. In an effort to develop a quantitative understanding of the effect of oxygen on product distribution and reaction rate, we have identified first generation products and measured reaction rates at four oxygen levels ranging from 0.01% to 100% saturation. In addition to bromochloroacetone and allyl chloride, a number of bromochloropropene isomers were identified in the presence of oxygen and are thought to be formed by an elimination mechanism. These products accounted for greater than 97 mol % of the reacted DBCP, which was run to high conversion (60-100 mol % DBCP converted). These measurements suggest that P450cam acts on the DBCP substrate through hydroxylation to produce 1-bromo-3-chloroacetone, through reduction to produce allyl chloride, and through elimination to produce bromochloropropene, with oxygen concentration determining the extent of each activity. A global data fitting kinetic model that describes the time-varying concentrations of substrate and products was developed to quantify the controlling level of oxygen on these multiple activities. The parameters of the model were compared with independent measurements and data from the literature.

Freeze-quenched iron-oxo intermediates in cytochromes P450

Since the discovery of cytochromes P450 and their assignment to heme proteins a reactive iron-oxo intermediate as the hydroxylating species has been discussed. It is believed that the electronic structure of this intermediate corresponds to an iron(IV)-porphyrin-pi-cation radical system (Compound I). To trap this intermediate the reaction of P450 with oxidants (shunt pathway) has been used. The common approaches are stopped-flow experiments with UV-visible spectroscopic detection or rapid-mixing/freeze-quench studies with EPR and Mössbauer spectroscopic characterization of the trapped intermediate. Surprisingly, the two approaches seem to give conflicting results. While the stopped-flow data indicate the formation of a porphyrin-pi-cation radical, no such species is seen by EPR spectroscopy, although the Mössbauer data indicate iron(IV) for P450cam (CYP101) and P450BMP (CYP102). Instead, radicals on tyrosine and tryptophan residues are observed. These findings are reviewed and discussed with respect to intramolecular electron transfer from aromatic amino acids to a presumably transiently formed porphyrin-pi-cation radical.

Direct Observation of a Novel Perturbed Oxyferrous Catalytic Intermediate during Reduced Putidaredoxin-initiated Turnover of Cytochrome P-450-CAM

The single turnover of (1R)(+)-camphor-bound oxyferrous cytochrome P450-CAM with one equivalent of dithionite-reduced putidaredoxin (Pdx) was monitored for the appearance of transient intermediates at 3 degrees C by double mixing rapid scanning stopped-flow spectroscopy. With excess camphor, three successive species were observed after generating oxyferrous P450-CAM and reacting versus reduced Pdx: a perturbed oxyferrous derivative, a species that was a mixture of high and low spin Fe(III), and high spin ferric camphor-bound enzyme. The rates of the first two steps, approximately 140 and approximately 85 s(-1), were assigned to formation of the perturbed oxyferrous intermediate and to electron transfer from reduced Pdx, respectively. In the presence of stoichiometric substrate, three phases with similar rates were seen even though the final state is low spin ferric P450-CAM. This is consistent with substrate being hydroxylated during the reaction. The single turnover reaction initiated by adding dioxygen to a preformed reduced P450-CAM.Pdx complex with excess camphor also led to phases with similar rates. It is proposed that formation of the perturbed oxyferrous intermediate reflects alteration of H-bonding to the proximal Cys, increasing the reduction potential of the oxyferrous state and triggering electron transfer from reduced Pdx. This species may be a direct spectral signature of the effector role of Pdx on P450-CAM reactivity (i.e. during catalysis). The substrate-free oxyferrous enzyme also reacted readily with reduced Pdx, showing that the inability of substrate-free P450-CAM to accept electrons from reduced Pdx and function as an NADH oxidase is completely due to the incapacity of reduced Pdx to deliver the first but not the second electron.

Structure and redox properties of the haem centre in the C357M mutant of cytochrome P450cam

The effects of site-specific mutation of the axial cysteine (C357M) to a methionine residue in cytochrome P450cam on the enzyme's coordination geometry and redox potential have been investigated. The absorption spectra of the haem centre in the C357M mutant of the enzyme showed close similarity to those of cytochrome c both in the oxidised and reduced forms. A well-defined absorption peak at 695 nm, similar to that seen in the case of cytochrome c and characteristic of methionine ligation to the ferric haem, was observed. The results indicated that the haem of C357M cytochrome P450cam is possibly axially coordinated to a methionine and a histidine, analogously to cytochrome c. The circular dichroism spectra in the visible and the far-UV regions suggested that the tertiary structure of the haem cavity in the C357M mutant cytochrome P450cam was distinctly different from that in the wild-type enzyme or in cytochrome c, although the secondary structure of the mutant remained identical to that of the wild-type cytochrome P450cam. Comparison of the natures of the CD spectra in the 400 nm and 695 nm regions of the C357M mutant of cytochrome P450cam with those of horse cytochrome c suggested (R) chirality at the sulfur atom of the iron-bound methionine residue in the mutant. The redox potential of the haem centre, estimated by redox titration of the C357M mutant, was found to be +260 mV, which is much higher than that in the wild-type enzyme and similar to the redox potential of cytochrome c. This supported the concept that axial ligation of the haem plays the major role in tuning the redox potential of the haem centre in haem proteins.

An artificial electron donor supported catalytic cycle of Pseudomonas putida cytochrome P450cam

Putidaredoxin (PdX), the physiological effector of cytochrome P450cam (P450cam), serves to gate electron transfer into oxy-P450cam during the catalytic cycle of the enzyme. Redox-linked structural changes in PdX are necessary for the effective P450cam turnover reaction. PdX is believed to be difficult to be replaced by an artificial electron donor in the reaction pathway of P450cam. We demonstrate that the catalytic cycle of wild-type P450cam can be supported in the presence of an artificial reductant, potassium ferrocyanide. Upon rapid mixing of ferrocyanide ion with P450cam, we observed an intermediate with spectral features characteristic of compound I. The rate constant for the formation of compound I in the presence of ferrocyanide supported reaction cycle was found to be comparable to the ones observed for H2O2 supported compound I formation in wild-type P450cam, but was much lower than those observed for classical peroxidases. The results presented in this paper form the first kinetic analysis of this intermediate for an artificial electron-driven P450cam catalytic pathway in solution.

Crystallographic Study on the Dioxygen Complex of Wild-type and Mutant Cytochrome P450cam

Two key amino acids, Thr252 and Asp251, are known to be important for dioxygen activation by cytochrome P450cam. We have solved crystal structures of a critical intermediate, the ferrous dioxygen complex (Fe(II)-O2), of the wild-type P450cam and its mutants, D251N and T252A. The wild-type dioxygen complex structure is very much the same as reported previously (Schlichting, I., Berendzen, J., Chu, K., Stock, A. M., Maves, S. A., Benson, D. E., Sweet, R. M., Ringe, D., Petsko, G. A., and Sligar, S. G. (2000) Science 287, 1615-1622) with the exception of higher occupancy and a more ordered structure of the iron-linked dioxygen and two "catalytic" water molecules that form part of a proton relay system to the iron-linked dioxygen. Due to of the altered conformation of the I helix groove these two waters are missing in the D251N dioxygen complex which explains its lower catalytic activity and slower proton transfer to the dioxygen ligand. Similarly, the T252A mutation was expected to disrupt the active site solvent structure leading to hydrogen peroxide formation rather than substrate hydroxylation. Unexpectedly, however, the two "catalytic" waters are retained in the T252A mutant. Based on these findings, we propose that the Thr(252) accepts a hydrogen bond from the hydroperoxy (Fe(III)-OOH) intermediate that promotes the second protonation on the distal oxygen atom, leading to O-O bond cleavage and compound I formation.

Detection of a high-barrier conformational change in the active site of cytochrome P450cam upon binding of putidaredoxin

The orientation of the substrate camphor in the active site of reduced CO-bound cytochrome P450cam (CYP101) as a function of reduced putidaredoxin (Pdxr) addition has been examined by NMR using perdeuterated CYP101 and perdeuterated Pdx as well as isotopically labeled d-camphor. This permits the 1H resonances of CYP101-bound camphor to be observed without interference from the signals of CYP101 or Pdx and confirms assignments of the methyl signals of camphor in the bound form. The Cys4Fe2S2 ferredoxin Pdx is the physiological redox partner and effector of CYP101. The addition of Pdx to the reduced CYP101-camphor-CO complex results in a conformational selection that is slow on the chemical shift time scale with spectral effects observed primarily at the 8-CH3 group of the camphor. The camphor signals are ring current shifted by the heme, and for the 9- and 10-CH3 resonances, these shifts are reasonably well predicted by ring current calculations from the crystal structure of CO-bound CYP101. However, in the absence of Pdx, the 8-CH3 resonance of CYP101-bound camphor is observed at considerably higher field than predicted. Dynamic simulations using ring current shift restraints generated a structure with low chemical shift violations in which the hydrogen bond between the camphor carbonyl oxygen and the OH of Tyr96 is lost, and an expansion of the active site takes place that permits reorientation of the camphor within the active site.

The P450cam G248E mutant covalently binds its prosthetic heme group

Previous studies on mammalian peroxidases and cytochrome P450 family 4 enzymes have shown that a carboxylic group positioned close to a methyl group of the prosthetic heme is required for the formation of a covalent link between a protein carboxylic acid side chain and the heme. To determine whether there are additional requirements for covalent bond formation in the P450 enzymes, a glutamic acid or an aspartic acid has been introduced into P450(cam) close to the heme 5-methyl group. Spectroscopic and kinetic studies of the resulting G248E and G248D mutants suggest that the carboxylate group coordinates with the heme iron atom, as reported for a comparable P450(BM3) mutant [Girvan, H. M., Marshall, K. R., Lawson, R. J., Leys, D., Joyce, M. G., Clarkson, J., Smith, W. E., Cheesman, M. R., and Munro, A. W. (2004) J. Biol. Chem. 279, 23274-23286]. The two P450(cam) mutants have low catalytic activity, but in contrast to the P450(BM3) mutant, incubation of the G248E (but not G248D) mutant with camphor, putidaredoxin, putidaredoxin reductase, and NADH results in partial covalent binding of the heme to the protein. No covalent attachment is observed in the absence of camphor or any of the other reaction components. Pronase digestion of the G248E P450(cam) mutant after covalent attachment of the heme releases 5-hydroxyheme, establishing that the heme is covalently attached through its 5-methyl group as predicted by in silico modeling. The results establish that a properly positioned carboxyl group is the sole requirement for autocatalytic formation of a heme-protein link in P450 enzymes, but also show that efficient covalent binding requires placement of the carboxyl close to the methyl but in a manner that prevents strong coordination to the iron atom.

The influence of substrate on the spectral properties of oxyferrous wild-type and T252A cytochrome P450-CAM

To probe whether the nature of the substrate can directly influence the spectral properties of oxyferrous cytochrome P450-CAM, the complex has been investigated in the absence and in the presence of the natural substrate (1R)-camphor (camphor) and of several camphor analogs. The oxyferrous complex of T252A P450-CAM, a mutant lacking the hydroxyl group that forms a hydrogen bond to the heme iron-coordinated dioxygen, has also been studied to gauge the influence of this hydrogen bond. UV-visible absorption and magnetic circular dichroism (MCD) spectra of these oxyferrous adducts prepared and stabilized at -40 degrees C in 60% (v/v) ethylene glycol are generally similar, exhibiting absorption bands at approximately 355, approximately 420, approximately 554, and approximately 585 nm (shoulder) and a characteristic MCD trough at approximately 585 nm. The MCD spectrum of camphor-bound oxyferrous P450-CAM is similar to that of the substrate-free oxyferrous enzyme, but the spectrum of the oxyferrous enzyme differs detectably in the presence of substrate analogs. The spectra of the oxyferrous T252A mutant and wild-type enzyme are overall similar except for Soret band position blue shifts by 2-6 nm for the mutant. 5-Methylenylcamphor (epoxidation substrate) appears to have an anomalous binding mode for the mutant compared with that for the wild-type enzyme. The present results indicate that the structures of the camphor analogs can sensitively influence the physical (spectroscopic) properties of the P450 dioxygen complex and could also affect its reactivity. The ability of substrate to modulate the reactivity of P450 intermediates could be a relevant factor in explaining the remarkable diversity of reactions catalyzed by the enzyme.

Reaction of Ferric Cytochrome P450cam with Peracids

Reactions of substrate-free ferric cytochrome P450cam with peracids to generate Fe=O intermediates have previously been investigated with contradictory results. Using stopped-flow spectrophotometry, the reaction with m-chloroperoxybenzoic acid demonstrated an Fe(IV)=O + porphyrin pi-cation radical (Cpd I) (Egawa, T., Shimada, H., and Ishimura, Y. (1994) Biochem. Biophys. Res. Commun. 201, 1464-1469). By contrast, with peracetic acid, Fe(IV)=O plus a tyrosyl radical were observed by freeze-quench Mossbauer and EPR spectroscopy (Schunemann, V., Jung, C., Trautwein, A. X., Mandon, D., and Weiss, R. (2000) FEBS Lett. 479, 149-154). Our detailed kinetic studies have resolved these contradictory results. At pH >7, a significant fraction of Cpd I is formed transiently, whereas at low pH only a species with a Soret band at 406 nm, presumably Fe(IV)=O + tyrosyl radical, is observed. Evidence for formation of an acylperoxo complex en route to Cpd I was obtained. Because of rapid heme destruction, steps subsequent to formation of the highly oxidized forms could not be fully characterized. Heme destruction was avoided by including peroxidase substrates (e.g. guaiacol), which were oxidized to characteristic peroxidase products as the Fe(III)-P450 was regenerated. Addition of ascorbate to either of the high valent species also reforms the Fe(III) state with only a small loss of heme absorbance. These results indicate that typical peroxidase chemistry occurs with P450cam and offer an explanation for the contrasting results reported earlier. The delineation of improved conditions (pH, temperature, choice of peracid) for generating highly oxidized species with P450cam should be valuable for their further characterization.

Substrate modulation of the properties and reactivity of the oxy-ferrous and hydroperoxo-ferric intermediates of cytochrome P450cam as shown by cryoreduction-EPR/ENDOR spectroscopy

EPR/ENDOR studies have been carried out on oxyferrous cytochrome P450cam one-electron cryoreduced by gamma-irradiation at 77 K in the absence of substrate and in the presence of a variety of substrates including its native hydroxylation substrate, camphor (a), and the alternate substrates, 5-methylenyl-camphor (b), 5,5-difluorocamphor (c), norcamphor (d), and adamantanone (e); the equivalent experiments have been performed on the T252A mutant complexed with a and b. The present study shows that the properties and reactivity of the oxyheme and of both the primary and the annealed intermediates are modulated by a bound substrate. This includes alterations in the properties of the heme center itself (g tensor; (14)N, (1)H, hyperfine couplings). It also includes dramatic changes in reactivity: the presence of any substrate increases the lifetime of hydroperoxoferri-P450cam (2) no less than ca. 20-fold. Among the substrates, b stands out as having an exceptionally strong influence on the properties and reactivity of the P450cam intermediates, especially in the T252A mutant. The intermediate, 2(T252A)-b, does not lose H(2)O(2), as occurs with 2(T252A)-a, but decays with formation of the epoxide of b. Thus, these observations show that substrate can modulate the properties of both the monoxygenase active-oxygen intermediates and the proton-delivery network that encompasses them.

Proton transfer at helium temperatures during dioxygen activation by heme monooxygenases

In the first measurement of enzymatic proton transfer at liquid helium temperatures, we examine protonation of the peroxo-ferriheme state of heme oxygenase (HO) produced by in situ radiolytic cryoreduction of oxy-HO in H2O and D2O solvents at ca. 4 K and above, and compare these findings with analogous measurements for oxy-P450cam and for oxy-Mb. Proton transfer in HO occurs at helium temperatures in both solvents; it occurs in P450cam at approximately 50 K and higher; in Mb it does not occur until T > 170 K. For Mb, this transfer at 180 K is biphasic, and the majority phase shows a solvent kinetic isotope effect of 3.8. We discuss these results in the context of the picture of environmentally coupled tunneling, which links proton transfer to two classes of protein motions: environmental reorganization (lambda in Marcus-like equations) and protein fluctuations ("active dynamics"; gating) which modulate the distance of proton transfer.

Biotransformation of the sesquiterpene (+)-valencene by cytochrome P450cam and P450BM-3

The sesquiterpenoids are a large class of naturally occurring compounds with biological functions and desirable properties. Oxidation of the sesquiterpene (+)-valencene by wild type and mutants of P450cam from Pseudomonas putida, and of P450BM-3 from Bacillus megaterium, have been investigated as a potential route to (+)-nootkatone, a fine fragrance. Wild type P450cam did not oxidise (+)-valencene but the mutants showed activities up to 9.8 nmol (nmol P450)(-1) min(-1), with (+)-trans-nootkatol and (+)-nootkatone constituting >85% of the products. Wild type P450BM-3 and mutants had higher activities (up to 43 min(-1)) than P450cam but were much less selective. Of the many products, cis- and trans-(+)-nootkatol, (+)-nootkatone, cis-(+)-valencene-1,10-epoxide, trans-(+)-nootkaton-9-ol, and (+)-nootkatone-13S,14-epoxide were isolated from whole-cell reactions and characterised. The selectivity patterns suggest that (+)-valencene has one binding orientation in P450cam but multiple orientations in P450BM-3.

Cryoreduction EPR and 13C, 19F ENDOR study of substrate-bound substates and solvent kinetic isotope effects in the catalytic cycle of cytochrome P450cam and its T252A mutant

We recently used cryoreduction EPR/ENDOR techniques to show that a substrate can modulate the properties of both the monooxygenase active-oxygen intermediates and of the proton-delivery network which encompasses them. In the present report we use Q-band pulsed 19F ENDOR (Mims 3-pulse sequence) to examine the substrate binding geometries of camphor, through use of the 5,5'--difluorocamphor, and 13C ENDOR to examine the binding of 5-methylenyl camphor labeled with 13C at C11. These probes are examined in multiple states of the catalytic cycle of P450cam and its T252A mutant. As part of this investigation we further report a new cryoreduction reaction, the reduction of a ferroheme to the EPR-visible Fe(I) state, and use it to probe the substrate binding to the EPR-silent ferroheme state. Finally we report the solvent kinetic isotope effect on the decay of the camphor complex of the hydroperoxo-ferric intermediate, the first such measurement on an individual step within the P450cam reaction cycle. Following reduction of oxyferrous-P450cam, this step is the rate-limiting step in camphor hydroxylation, and its solv-KIE of 1.8 at 190 K establishes that it involves activation of the hydroperoxo moiety by transfer of the 'second' proton of catalysis. We suggest that the finding that the heme pocket can exist in multiple substates, including multiple substrate binding locations, even in P450cam, along with the established possibility that the hydroperoxo-ferriheme intermediate can react with substrate, may explain the formation of multiple products by P450s.

Identification of broad specificity P450CAM variants by primary screening against indole as substrate

High-throughput screening of cytochrome P450CAM libraries, for their ability to oxidise indole to indigo and indirubin, has resulted in the identification of variants with activity towards the structurally unrelated substrate diphenylmethane.

A predictive pattern of computed barriers for C-h hydroxylation by compound I of cytochrome p450

The communication presents DFT calculations of 10 different C-H hydroxylation barriers by the active species of the enzyme cytochrome P450. The work demonstrates the existence of an excellent barrier-bond energy correlation. The so-obtained equation of the straight line is demonstrated to be useful for predicting barriers of related C-H activation processes, as well as for assessing barrier heights within the protein environment. This facility is demonstrated be estimating the barrier of camphor hydroxylation by P450cam.

Structural stability and dynamics of hydrogenated and perdeuterated cytochrome P450cam (CYP101)

Perdeuterated and hydrogenated cytochrome P450cam (P450cam), from Pseudomonas putida, has been characterized concerning thermal stability and structural dynamics. For the first time, Fourier transform infrared (FTIR) spectroscopy was used to characterize a perdeuterated protein. The secondary structure compositions were determined from the fitted amide I' spectral region, giving band populations at 10 degrees C for the perdeuterated protein of 22% between 1605 and 1624 cm(-1) (beta-sheets), 47% between 1633 and 1650 cm(-1) (alpha-helix (29%) plus unordered/3(10)-helix (18%)), and 28% between 1657 and 1677 cm(-1) (turns) and for the hydrogenated protein of 22% between 1610 and 1635 cm(-1) (beta-sheets), 52% between 1640 and 1658 cm(-1) (alpha-helix (41%) plus unordered/3(10)-helix (11%)), and 24% between 1665 and 1680 cm(-1) (turns). Thermal unfolding experiments revealed that perdeuterated P450cam was less stable than the hydrogenated protein. The midpoint transition temperatures were 60.8 and 64.4 degrees C for the perdeuterated and hydrogenated P450cam, respectively. Step-scan time-resolved FTIR was applied to the P450cam-CO complex to study the ligand-rebinding process after flash photolysis. Rebinding of the ligand occurred with the same kinetics and rate constants k(on), 8.9 x 10(4) and 8.3 x 10(4) M(-1) s(-1) for the perdeuterated and hydrogenated P450cam, respectively.Perdeuterated P450cam was expressed for a neutron crystallographic study to determine the specific hydration states and hydrogen-bonding networks at the active site. The analyses presented here show that perdeuterated P450cam is structurally similar to its hydrogenated counterpart, despite its reduced thermal stability, suggesting that information obtained from the neutron structure will be representative of the normal hydrogenated P450cam.

Substrate binding favors enhanced NO binding to P450cam

Ferric cytochrome P450cam from Pseudomonas putida (P450cam) in buffer solution at physiological pH 7.4 reversibly binds NO to yield the nitrosyl complex P450cam(NO). The presence of 1R-camphor affects the dynamics of NO binding to P450cam and enhances the association and dissociation rate constants significantly. In the case of the substrate-free form of P450cam, subconformers are evident and the NO binding kinetics are much slower than in the presence of the substrate. The association and dissociation processes were investigated by both laser flash photolysis and stopped-flow techniques at ambient and high pressure. Large and positive values of S and V observed for NO binding to and release from the substrate-free P450cam complex are consistent with the operation of a limiting dissociative ligand substitution mechanism, where the lability of coordinated water dominates the reactivity of the iron(III)-heme center with NO. In contrast, NO binding to P450cam in the presence of camphor displays negative activation entropy and activation volume values that support a mechanism dominated by a bond formation process. Volume profiles for the binding of NO appear to be a valuable approach to explain the differences observed for P450cam in the absence and presence of the substrate and enable the clarification of the underlying reaction mechanisms at a molecular level. Changes in spin state of the iron center during the binding/release of NO contribute significantly to the observed volume effects. The results are discussed in terms of relevance for the biological function of cytochrome P450 and in context to other investigations of the related reactions between NO and imidazole- and thiolate-ligated iron(III) hemoproteins.

Radical rebound mechanism in cytochrome P-450-catalyzed hydroxylation of the multifaceted radical clocks alpha- and beta-thujone

Alpha-thujone (1alpha) and beta-thujone (1beta) were used to investigate the mechanism of hydrocarbon hydroxylation by cytochromes P-450(cam) (CYP101) and P-450(BM3) (CYP102). The thujones are hydroxylated by these enzymes at various positions, but oxidation at C-4 gives rise to both rearranged and unrearranged hydroxylation products. Rearranged products result from the formation of a radical intermediate that can undergo either inversion of stereochemistry or ring opening of the adjacent cyclopropane ring. Both of these rearrangements, as well as a C-4 desaturation reaction, are observed. The ring opening clock gives oxygen rebound rates that range from 0.2 x 10(10) to 2.8 x 10(10) s(-1) for the different substrate and enzyme combinations. The C-4 inversion reaction provides independent confirmation of a radical intermediate. The phenol product expected if a C-4 cationic rather than radical intermediate is formed is not detected. The results are consistent with a two-state process and provide support for a radical rebound but not a hydroperoxide insertion mechanism for cytochrome P-450 hydroxylation.