Understanding the role of the essential Asp251 in cytochrome p450cam using site-directed mutagenesis, crystallography, and kinetic solvent isotope effect

Molecular oxygen activation and proton transfer mechanisms in lanosterol 14alpha-demethylase catalysis

The CYP51 lanosterol 14alpha-demethylases are evolutionarily ancient enzymes ubiquitously distributed throughout the biological domains. The experimental X-ray crystal structure of Mycobacterium tuberculosis (Mtb) CYP51 is the first of an enzyme capable of catalyzing inert C-C bond cleavage. Amino acid sequence comparisons of CYP51 family members with other members of the CYP superfamily reveal the almost universally conserved "acid-alcohol" pair, putatively involved in proton transport and O(2) activation, is replaced with a His-Thr dyad. In this study, extended molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics calculations (QM/MM) are applied to characterize reactive oxygen intermediates and to unravel mechanisms of O(2) activation vis-a-vis proton transport for this important enzyme. MD confirms stable His259deltaH(+)-Thr260OH-O(2) (Mtb numbering) hydrogen bonding early in the simulations, suggesting these amino acids could function similarly to the Asp251-Thr252 pair in CYP101. QM/MM calculations support this dyad competently catalyzes the peroxo to Compound 0 (Cmpd 0) reaction, albeit an endothermic homolytic O-O scission mechanism affording Compound I (Cmpd I) was identified. Disruption of the His259H(+)-Thr260OH hydrogen bond in MD simulation divulges a second previously unidentified hydrogen-bond network, including three water molecules linking Glu173 in the CYP51 F-helix to the distal O(2) atom. Expansion of the QM region to contain these atoms unveils an unprecedented triradicaloid electronic structure of the peroxo intermediate characterized by spin polarization to the Glu173 side chain, attributable to the protein electrostatic environment. This amino acid, in concert with an active-site water network, catalyzes a facile protonation of the peroxo intermediate and offers a series of redundant heterolytic and homolytic mechanisms, affording exothermic formation of the ultimate oxidant Cmpd I. In summary, this study highlights the importance of the protein electrostatic environment to tune the electronic structure of CYP catalytic intermediates in addition to Cmpd I and illustrates the diversity of proton transport pathways available to these enzymes to drive catalysis.

How is the reactivity of cytochrome P450cam affected by Thr252X mutation? A QM/MM study for X = serine, valine, alanine, glycine

Proton transfer reactions play a vital role in the catalytic cycle of cytochrome P450cam and are responsible for the formation of the iron-oxo species called Compound I (Cpd I) that is supposed to be the active oxidant. Depending on the course of the proton transfer, protonation of the last observable intermediate (ferric hydroperoxo complex, Cpd 0) can lead to either the formation of Cpd I (coupling reaction) or the ferric resting state (uncoupling reaction). The ratio of these two processes is drastically affected by mutation of the Thr252 residue. In this work, we study the effect of Thr252X (X = serine, valine, alanine, glycine) mutations on the formation of Cpd I by means of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations and classical simulations. In the wild-type enzyme, the coupling reaction is favored since its rate-limiting barrier is 13 kcal/mol lower than that for uncoupling. This difference is reduced to 7 kcal/mol in the serine mutant. In the case of valine, alanine, and glycine mutants, an additional water molecule enters the active site and lowers the activation energy of the uncoupling reaction significantly. With the additional water molecule, coupling and uncoupling have similar barriers in the valine mutant, and the uncoupling reaction becomes favored in the alanine and glycine mutants. These findings agree very well with experimental results and thus confirm the assumption that uncontrolled proton delivery by solvent water networks is responsible for the uncoupling reaction. The present study provides a detailed mechanistic understanding of the role of the Thr252 residue.

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].

Oxygen activation by cytochrome P450 monooxygenase

Unlike photosystem II (PSII) that catalyzes formation of the O-O bond, the cytochromes P450 (P450), members of a superfamily of hemoproteins, catalyze the scission of the O-O bond of dioxygen molecules and insert a single oxygen atom into unactivated hydrocarbons through a hydrogen abstraction-oxygen rebound mechanism. Hydroxylation of the unactivated hydrocarbons at physiological temperatures is vital for many cellar processes such as the biosynthesis of many endogenous compounds and the detoxification of xenobiotics in humans and plants. Even though it carries out the opposite of the water splitting reaction, P450 may share similarities to PSII in proton delivery networks, oxygen and water access channels, and consecutive electron transfer processes. In this article, we review recent advances in understanding the molecular mechanisms by which P450 activates dioxygen.

EPR and ENDOR studies of cryoreduced compounds II of peroxidases and myoglobin. Proton-coupled electron transfer and protonation status of ferryl hemes

The nature of the [Fe(IV)-O] center in hemoprotein Compounds II has recently received considerable attention, as several experimental and theoretical investigations have suggested that this group is not necessarily the traditionally assumed ferryl ion, [Fe(IV)=O]2+, but can be the protonated ferryl, [Fe(IV)-OH]3+. We show here that cryoreduction of the EPR-silent Compound II by gamma-irradiation at 77 K produces Fe(III) species retaining the structure of the precursor [Fe(IV)=O]2+ or [Fe(IV)-OH]3+, and that the properties of the cryogenerated species provide a report on structural features and the protonation state of the parent Compound II when studied by EPR and 1H and 14N ENDOR spectroscopies. To give the broadest view of the properties of Compounds II we have carried out such measurements on cryoreduced Compounds II of HRP, Mb, DHP and CPO and on CCP Compound ES. EPR and ENDOR spectra of cryoreduced HRP II, CPO II and CCP ES are characteristic of low-spin hydroxy-Fe(III) heme species. In contrast, cryoreduced "globins", Mb II, Hb II, and DHP II, show EPR spectra having lower rhombicity. In addition the cryogenerated ferric "globin" species display strongly coupled exchangeable (1)H ENDOR signals, with A max approximately 20 MHz and a iso approximately 14 MHz, both substantially greater than for hydroxide/water ligand protons. Upon annealing at T > 180 K the cryoreduced globin compounds II relax to the low-spin hydroxy-ferric form with a solvent kinetic isotope effect, KIE > 6. The results presented here together with published resonance Raman and Mossbauer data suggest that the high-valent iron center of globin and HRP compounds II, as well as of CCP ES, is [Fe(IV)=O]2+, and that its cryoreduction produces [Fe(III)-O]+. Instead, as proposed by Green and co-workers, CPO II contains [Fe(IV)-OH]3+ which forms [Fe(III)-OH]2+ upon radiolysis. The [Fe(III)-O]+ generated by cryoreduction of HRP II and CCP ES protonate at 77 K, presumably because the heme is linked to a distal-pocket hydrogen bonding/proton-delivery network through an H-bond to the "oxide" ligand. The data also indicate that Mb and HRP compounds II exist as two major conformational substates.

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.

Combined QM/MM calculations of active-site vibrations in binding process of P450cam to putidaredoxin

Combined QM/MM calculations of the active-site of cytochrome P450cam have been performed before and after the binding of P450cam to putidaredoxin. The calculations were carried out for both a 5-coordinated and a 6-coordinated active-site of cytochrome P450cam, with either a water molecule or a carbon monoxide molecule as a 6th distal ligand. An experimentally observed increase in the Fe-S stretching frequency that occurs after cytochrome P450cam binds to putidaredoxin, has been reproduced in our study. Experimentally observed changes in the Fe-C and C-O vibration frequencies that occur after binding of both proteins, have also been reproduced in our study. The computed increase of the Fe-S and Fe-C stretching frequencies is correlated with a corresponding decrease of the Fe-S and Fe-C interatomic distances. According to our calculations, for the active-site with carbon monoxide in the triplet electronic state, the binding process increases the spin densities on the iron and sulfur atoms, which changes the Fe-C and C-O stretching frequencies in opposite directions, in agreement with experimental data.

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.

Proton assisted oxygen-oxygen bond splitting in cytochrome p450

Proton assisted O-O bond splitting of cytochromes' P450 hydroperoxo Compound 0 has been investigated by density functional theory, showing a barrier for the slightly endothermic formation of the iron-oxo Compound I. The barrier and the endothermicity increase with decreasing acidity of the distal proton source. Protonation of the proximal iron heme ligand favors the O-O bond scission and provides an important regulatory component in the catalytic cycle. The Compound 0 --> I conversion is slightly exothermic for the peroxidase and catalase models. Implications of the energetic relationship between the two reactive intermediates are discussed in terms of possible oxidative pathways.

Proton-shuffle mechanism of O-O activation for formation of a high-valent oxo-iron species of bleomycin

Bleomycins (BLMs) can utilize H2O2 to cleave DNA in the presence of ferric ions. DFT calculations were used to study the mechanism of O-O bond cleavage in the low-spin FeIII-hydroperoxo complex of BLM. The following alternative hypotheses were investigated using realistic structural models: (a) heterolytic cleavage of the O-O bond, generating a Compound I (Cpd I) like intermediate, formally BLM-FeV=O; (b) homolytic O-O cleavage, leading to a BLM-FeIV=O species and an OH* radical; and (c) a direct O-O cleavage/H-abstraction mechanism by ABLM. The calculations showed that (a) is a facile and viable mechanism; it involves acid-base proton reshuffle mediated by the side-chain linkers of BLM, causing thereby heterolytic cleavage of the O-O bond and generation of Cpd I. Formation of Cpd I is found to involve a barrier of 13.3 kcal/mol, which is lower than the barriers in the alternative mechanisms (b and c) that possess respective barriers of 31 and 17 kcal/mol. The so-formed Cpd I species with a radical on the side-chain linker, methylvalerate (V), adjacent to the BLM-FeIV=O complex, resembles the formation of the active species of cytochrome c peroxidase in the Poulos-Kraut proton-shuffle mechanism in heme peroxidases (Poulos, T. L.; Kraut, J. J. Biol. Chem. 1980, 255, 8199-8205). Experimental data are discussed and shown to be in accord with this proposal. It suggests that the high-valence Cpd I species of BLM participates in the DNA cleavage. This is an alternative mechanistic hypothesis to the exclusive reactivity scenario based on ABLM (FeIII-OOH).

Proton-coupled electron transfer: the mechanistic underpinning for radical transport and catalysis in biology

Charge transport and catalysis in enzymes often rely on amino acid radicals as intermediates. The generation and transport of these radicals are synonymous with proton-coupled electron transfer (PCET), which intrinsically is a quantum mechanical effect as both the electron and proton tunnel. The caveat to PCET is that proton transfer (PT) is fundamentally limited to short distances relative to electron transfer (ET). This predicament is resolved in biology by the evolution of enzymes to control PT and ET coordinates on highly different length scales. In doing so, the enzyme imparts exquisite thermodynamic and kinetic controls over radical transport and radical-based catalysis at cofactor active sites. This discussion will present model systems containing orthogonal ET and PT pathways, thereby allowing the proton and electron tunnelling events to be disentangled. Against this mechanistic backdrop, PCET catalysis of oxygen-oxygen bond activation by mono-oxygenases is captured at biomimetic porphyrin redox platforms. The discussion concludes with the case study of radical-based quantum catalysis in a natural biological enzyme, class I Escherichia coli ribonucleotide reductase. Studies are presented that show the enzyme utilizes both collinear and orthogonal PCET to transport charge from an assembled diiron-tyrosyl radical cofactor to the active site over 35A away via an amino acid radical-hopping pathway spanning two protein subunits.

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.

Sterol 14alpha-demethylase cytochrome P450 (CYP51), a P450 in all biological kingdoms

The CYP51 family is an intriguing subject for fundamental P450 structure/function studies and is also an important clinical drug target. This review updates information on the variety of the CYP51 family members, including their physiological roles, natural substrates and substrate preferences, and catalytic properties in vitro. We present experimental support for the notion that specific conserved regions in the P450 sequences represent a CYP51 signature. Two possible roles of CYP51 in P450 evolution are discussed and the major approaches for CYP51 inhibition are summarized.

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.

Gauging the Relative Oxidative Powers of Compound I, Ferric-Hydroperoxide, and the Ferric-Hydrogen Peroxide Species of Cytochrome P450 Toward C−H Hydroxylation of a Radical Clock Substrate

Density functional calculations were performed in response to the controversies regarding the identity of the oxidant species in cytochrome P450. The calculations were used to gauge the relative C-H hydroxylation reactivity of three potential oxidant species of the enzyme, the high-valent oxo-iron species Compound I (Cpd I), the ferric hydroperoxide Compound 0 (Cpd 0), and the ferric-hydrogen peroxide complex Fe(H(2)O(2)). The results for the hydroxylation of a radical probe substrate, 1, show the following trends: (a) Cpd I is the most reactive species; in its presence the other two reagents will be silent. (b) In the absence of Cpd I, substrate oxidation by Cpd 0 and Fe(H(2)O(2)) will take place via a stepwise mechanism that involves initial O-O homolysis followed by H-abstraction from 1. (c) Cpd 0 will undergo mostly porphyrin hydroxylation and only approximately 15% of substrate oxidation producing mostly the rearranged alcohol, 3 (Scheme 2). (d) Fe(H(2)O(2)) will generate mostly free hydrogen peroxide (uncoupling). A small fraction will perform substrate oxidation and lead mostly to 3. Reactivity probes for these reagents are kinetic isotope effect (KIE) and the product ratio of unrearranged to rearranged alcohols, [2/3]. Thus, for substrate oxidation by Cpd 0 or Fe(H(2)O(2)) KIE will be small, approximately 2, while Cpd I will have large KIE values. Typically both Cpd 0 and Fe(H(2)O(2)) will lead to a [2/3] ratio < 1, while Cpd I will lead to ratios > 1. In addition, the product isotope effect (KIE(2)/KIE(3) not equal 1) is expected from the reactivity of Cpd I.

Biodiversity of cytochrome P450 redox systems

P450s (cytochrome P450 mono-oxygenases) are a superfamily of haem-containing mono-oxygenase enzymes that participate in a wide range of biochemical pathways in different organisms from all of the domains of life. To facilitate their activity, P450s require sequential delivery of two electrons passed from one or more redox partner enzymes. Although the P450 enzymes themselves show remarkable similarity in overall structure, it is increasingly apparent that there is enormous diversity in the redox partner systems that drive the P450 enzymes. This paper examines some of the recent advances in our understanding of the biodiversity of the P450 redox apparatus, with a particular emphasis on the redox systems in the pathogen Mycobacterium tuberculosis.

Structural biology of heme monooxygenases

Over the past few years the number of crystal structures available for heme monooxygenases has substantially increased. Those most closely related to one another are cytochrome P450, nitric oxide synthase, and heme oxygenase. The present mini-review provides a summary of some recently published work on how crystallography and solution studies have provided new insights on function and especially the oxygen activation process. It now appears that in all three monooxygenases highly ordered solvent in the active site serves as direct proton donors to the iron-linked dioxygen; a requirement for splitting the O-O bond. This is in sharp contrast to the related peroxidase family of enzymes where strategically positioned amino acid side chains serve the function of shuttling protons. The P450cam-oxy-complex as well as various mutants in a complex with either oxygen or carbon monoxide have enabled a fairly detailed picture to be developed on the role of specific amino acids and conformational changes in both electron transfer and oxygen activation.

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

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

New Features in the Catalytic Cycle of Cytochrome P450 during the Formation of Compound I from Compound 0

Density functional theory (DFT) is applied to the dark section of the catalytic cycle of the enzyme cytochrome P450, namely, the formation of the active species, Compound I (Cpd I), from the ferric-hydroperoxide species (Cpd 0) by a protonation-assisted mechanism. The chosen 96-atom model includes the key functionalities deduced from experiment: Asp(251), Thr(252), Glu(366), and the water channels that relay the protons. The DFT model calculations show that (a) Cpd I is not formed spontaneously from Cpd 0 by direct protonation, nor is the process very exothermic. The process is virtually thermoneutral and involves a significant barrier such that formation of Cpd I is not facile on this route. (b) Along the protonation pathway, there exists an intermediate, a protonated Cpd 0, which is a potent oxidant since it is a ferric complex of water oxide. Preliminary quantum mechanical/molecular mechanical calculations confirm that Cpd 0 and Cpd I are of similar energy for the chosen model and that protonated Cpd 0 may exist as an unstable intermediate. The paper also addresses the essential role of Thr(252) as a hydrogen-bond acceptor (in accord with mutation studies of the OH group to OMe).

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.

Water-Assisted Oxo Mechanism for Heme Metabolism

A mechanism of heme metabolism by heme oxygenase (HO) is discussed from B3LYP density functional theory calculations. The concerted OH group attack to the alpha-carbon by the iron-hydroperoxo species is investigated using a model with full protoporphyrin IX to confirm our previous conclusion that this species does not have sufficient oxidizing power for heme oxidation (J. Am. Chem. Soc. 2004, 126, 3672). Calculated activation energies and structures of the intermediates and transition state for this process remain unchanged from those for a small model with porphine in the previous study, which shows that the inclusion of the side chain of the porphyrin ring is not essential in describing the OH group transfer. The activation barrier for a direct oxo attack to the alpha-carbon by an iron-oxo model is calculated to be 49.8 kcal/mol, the barrier height of which looks very high for the enzymatic reaction under physiological conditions. This large activation energy is due to a highly bent porphyrin structure in the transition state. However, a bridging water molecule plays an important role in reducing the porphyrin distortion in the transition state, resulting in a remarkable decrease of the activation barrier to 13.9 kcal/mol. A whole-enzyme model with about 4000 atoms is constructed to elucidate functions of the protein environment in this enzymatic reaction using QM/MM calculations. The key water molecule is fixed in the protein environment to ensure the low-barrier and regioselective heme oxidation. A water-assisted oxo mechanism of heme oxidation by heme oxygenase is proposed from these calculational results.

Oxoiron(IV) porphyrin π-cation radical complexes with a chameleon behavior in cytochrome P450 model reactions

There is an intriguing, current controversy on the involvement of multiple oxidizing species in oxygen transfer reactions by cytochromes P450 and iron porphyrin complexes. The primary evidence for the "multiple oxidants" theory was that products and/or product distributions obtained in the catalytic oxygenations were different depending on reaction conditions such as catalysts, oxidants, and solvents. In the present work, we carried out detailed mechanistic studies on competitive olefin epoxidation, alkane hydroxylation, and C=C epoxidation versus allylic C-H hydroxylation in olefin oxygenation with in situ generated oxoiron(IV) porphyrin pi-cation radicals (1) under various reaction conditions. We found that the products and product distributions were markedly different depending on the reaction conditions. For example, 1 bearing different axial ligands showed different product selectivities in competitive epoxidations of cis-olefins and trans-olefins and of styrene and para-substituted styrenes. The hydroxylation of ethylbenzene by 1 afforded different products, such as 1-phenylethanol and ethylbenzoquinone, depending on the axial ligands of 1 and substrates. Moreover, the regioselectivity of C=C epoxidation versus C-H hydroxylation in the oxygenation of cyclohexene by 1 changed dramatically depending on the reaction temperatures, the electronic nature of the iron porphyrins, and substrates. These results demonstrate that 1 can exhibit diverse reactivity patterns under different reaction conditions, leading us to propose that the different products and/or product distributions observed in the catalytic oxygenation reactions by iron porphyrin models might not arise from the involvement of multiple oxidizing species but from 1 under different circumstances. This study provides strong evidence that 1 can behave like a "chameleon oxidant" that changes its reactivity and selectivity under the influence of environmental changes.

The Final Catalytic Step of Cytochrome P450 Aromatase: A Density Functional Theory Study

B3LYP density functional theory calculations are used to unravel the mysterious third step of aromatase catalysis. The feasibility of mechanisms in which the reduced ferrous dioxygen intermediate mediates androgen aromatization is explored and determined to be unlikely. However, proton-assisted homolysis of the peroxo hemiacetal intermediate to produce P450 compound I and the C19 gem-diol likely proceeds with a low energetic barrier. Mechanisms for the aromatization and deformylation sequence which are initiated by 1beta-hydrogen atom abstraction by P450 compound I are considered. 1beta-Hydrogen atom abstraction from substrates in the presence of the 2,3-enol encounters strikingly low barriers (5.3-7.8 kcal/mol), whereas barriers for this same process rise to 17.0-27.1 kcal/mol in the keto tautomer. Transition states for 1beta-hydrogen atom abstraction from enolized substrates in the presence of the 19-gem-diol decayed directly to the experimentally observed products. If the C19 aldehyde remains unhydrated, aromatization occurs with concomitant decarbonylation and therefore does not support dehydration of the C19 aldehyde prior to the final catalytic step. On the doublet surface, the transition state connects to a potentially labile 1(10) dehydrogenated product, which may undergo rapid aromatization, as well as formic acid. Ab initio molecular dynamics confirmed that the 1beta-hydrogen atom abstraction and deformylation or decarbonylation occur in a nonsynchronous, coordinated manner. These calculations support a dehydrogenase behavior of aromatase in the final catalytic step, which can be summarized by 1beta-hydrogen atom abstraction followed by gem-diol deprotonation.

Models and mechanisms of O-O bond activation by cytochrome P450. A critical assessment of the potential role of multiple active intermediates in oxidative catalysis

Cytochrome P450 enzymes promote a number of oxidative biotransformations including the hydroxylation of unactivated hydrocarbons. Whereas the long-standing consensus view of the P450 mechanism implicates a high-valent iron-oxene species as the predominant oxidant in the radicalar hydrogen abstraction/oxygen rebound pathway, more recent studies on isotope partitioning, product rearrangements with 'radical clocks', and the impact of threonine mutagenesis in P450s on hydroxylation rates support the notion of the nucleophilic and/or electrophilic (hydro)peroxo-iron intermediate(s) to be operative in P450 catalysis in addition to the electrophilic oxenoid-iron entity; this may contribute to the remarkable versatility of P450s in substrate modification. Precedent to this mechanistic concept is given by studies with natural and synthetic P450 biomimics. While the concept of an alternative electrophilic oxidant necessitates C-H hydroxylation to be brought about by a cationic insertion process, recent calculations employing density functional theory favour a 'two-state reactivity' scenario, implicating the usual ferryl-dependent oxygen rebound pathway to proceed via two spin states (doublet and quartet); state crossing is thought to be associated with either an insertion or a radicalar mechanism. Hence, challenge to future strategies should be to fold the disparate and sometimes contradictory data into a harmonized overall picture.

Structural aspects of ligand binding to and electron transfer in bacterial and fungal P450s

Cytochrome P450 enzymes are heme-containing monooxygenases that are named after an absorption band at 450 nm when complexed with carbon monoxide. They catalyze a wide variety of reactions and are unique in their ability to hydroxylate nonactivated hydrocarbons. P450 enzymes are involved in numerous biological processes, which include the biosynthesis of lipids, steroids, antibiotics, and the degradation of xenobiotics. In line with the variety of reactions catalyzed, the size of their substrates varies significantly. Some P450s have open active sites (e.g., BM3), and some have shielded active sites that open only transiently (e.g., P450cam), whereas others bind the substrate only when attached to carrier proteins (e.g., Oxy proteins). Structural aspects of both organic and gaseous ligand binding and electron transfer are described.

Formation of transient oxygen complexes of cytochrome p450 BM3 and nitric oxide synthase under high pressure

The kinetics of formation and transformation of oxygen complexes of two heme-thiolate proteins (the F393H mutant of cytochrome P450 BM3 and the oxygenase domain of endothelial nitric oxide synthase, eNOS) were studied under high pressure. For BM3, oxygen-binding characteristics (rate and activation volume) matched those measured for CO-binding. In contrast, pressure revealed a different CO- and oxygen-binding mechanism for eNOS, suggesting that it is hazardous to take CO-binding as a model for oxygen-binding. With eNOS, a ferric NO complex is formed as an intermediate in the second reaction cycle. Here we report the pressure stability of this compound. Furthermore, in the presence of 4-amino-tetrahydrobiopterin (ABH(4)), an analog to the natural second electron donor tetrahydrobiopterin (BH(4)), biphasic pressure profiles of the oxygen-binding rates were observed, both in the first and the second reaction cycles, indicative of the formation of an additional reaction intermediate. This was confirmed by experiments where ABH(4) was replaced by ABH(2), a cofactor which cannot deliver an electron. Altogether, high pressure appears to be a useful tool to characterize elementary steps in the reaction cycle of heme-thiolate proteins.

Kinetic isotope effects on the rate-limiting step of heme oxygenase catalysis indicate concerted proton transfer/heme hydroxylation

Heme oxygenase (HO) catalyzes the O2 and NADPH/cytochrome P450 reductase-dependent conversion of heme to biliverdin, free iron ion, and CO through a process in which the heme participates as both dioxygen-activating prosthetic group and substrate. We earlier confirmed that the first step of HO catalysis is a monooxygenation in which the addition of one electron and two protons to the HO oxy-ferroheme produces ferric-alpha-meso-hydroxyheme (h). Cryoreduction/EPR and ENDOR measurements further showed that hydroperoxo-ferri-HO converts directly to h in a single kinetic step without formation of a Compound I. We here report details of that rate-limiting step. One-electron 77 K cryoreduction of human oxy-HO and annealing at 200 K generates a structurally relaxed hydroperoxo-ferri-HO species, denoted R. We here report the cryoreduction/annealing experiments that directly measure solvent and secondary kinetic isotope effects (KIEs) of the rate-limiting R --> h conversion, using enzyme prepared with meso-deuterated heme and in H2O/D2O buffers to measure the solvent KIE (solv-KIE), and the secondary KIE (sec-KIE) associated with the conversion. This approach is unique in that KIEs measured by monitoring the rate-limiting step are not susceptible to masking by KIEs of other processes, and these results represent the first direct measurement of the KIEs of product formation by a kinetically competent reaction intermediate in any dioxygen-activating heme enzyme.The observation of both solv-KIE(298) = 1.8 and sec-KIE(298) = 0.8 (inverse) indicates that the rate-limiting step for formation of h by HO is a concerted process: proton transfer to the hydroperoxo-ferri-heme through the distal-pocket H-bond network, likely from a carboxyl group acting as a general acid catalyst, occurring in synchrony with bond formation between the terminal hydroperoxo-oxygen atom and the alpha-meso carbon to form a tetrahedral hydroxylated-heme intermediate. Subsequent rearrangement and loss of H2O then generates h.

Compressibility and uncoupling of cytochrome P450cam: high pressure FTIR and activity studies

The effect of the hydrostatic pressure on the CO ligand stretch vibration in cytochrome P450cam-CO bound with various substrates is studied by FTIR. The vibration frequency is linearily shifted to lower values with increasing pressure. The slope of the shift gives the isothermal compressibility of the heme pocket and is found to be related to the high-spin state content in an opposite direction to that previously observed from the pressure-induced shift of the Soret band. This opposite behaviour is explained by the dual effect of heme pocket water molecules both on the CO ligand and on electrostatic potentials produced by the protein at the distal side. The latter effect disturbs ligand-distal side contacts which are needed for a specific proton transfer in oxygen activation when dioxygen is the ligand. Their loss results in uncoupled H(2)O(2) formation.

A theoretical study on the mechanism of camphor hydroxylation by compound I of cytochrome p450

Mechanistic and energetic aspects for the conversion of camphor to 5-exo-hydroxycamphor by the compound I iron-oxo species of cytochrome P450 are discussed from B3LYP DFT calculations. This reaction occurs in a two-step manner along the lines that the oxygen rebound mechanism suggests. The activation energy for the first transition state of the H atom abstraction at the C5 atom of camphor is computed to be more than 20 kcal/mol. This H atom abstraction is the rate-determining step in this hydroxylation reaction, leading to a reaction intermediate that involves a carbon radical species and the iron-hydroxo species. The second transition state of the rebound step that connects the reaction intermediate and the product alcohol complex lies a few kcal/mol below that for the H atom abstraction on the doublet and quartet potential energy surfaces. This energetic feature allows the virtually barrierless recombination in both spin states, being consistent with experimentally observed high stereoselectivity and brief lifetimes of the reaction intermediate. The overall energetic profile of the catalytic mechanism of camphor hydroxylation particularly with respect to why the high activation energy for the H atom abstraction is accessible under physiological conditions is also considered and calculated. According to a proton source model involving Thr252, Asp251, and two solvent water molecules (Biochemistry 1998, 37, 9211), the energetics for the conversion of the iron-peroxo species to compound I is studied. A significant energy over 50 kcal/mol is released in the course of this dioxygen activation process. The energy released in this chemical process is an important driving force in alkane hydroxylation by cytochrome P450. This energy is used for the access to the high activation energy for the H atom abstraction.

The 1.92-A structure of Streptomyces coelicolor A3(2) CYP154C1. A new monooxygenase that functionalizes macrolide ring systems

Evolutionary links between cytochrome P450 monooxygenases, a superfamily of extraordinarily divergent heme-thiolate proteins catalyzing a wide array of NADPH/NADH- and O(2)-dependent reactions, are becoming better understood because of availability of an increasing number of fully sequenced genomes. Among other reactions, P450s catalyze the site-specific oxidation of the precursors to macrolide antibiotics in the genus Streptomyces introducing regiochemical diversity into the macrolide ring system, thereby significantly increasing antibiotic activity. Developing effective uses for Streptomyces enzymes in biosynthetic processes and bioremediation requires identification and engineering of additional monooxygenases with activities toward a diverse array of small molecules. To elucidate the molecular basis for substrate specificity of oxidative enzymes toward macrolide antibiotics, the x-ray structure of CYP154C1 from Streptomyces coelicolor A3(2) was determined (Protein Data Bank code ). Relocation of certain common P450 secondary structure elements, along with a novel structural feature involving an additional beta-strand transforming the five-stranded beta-sheet into a six-stranded variant, creates an open cleft-shaped substrate-binding site between the two P450 domains. High sequence similarity to macrolide monooxygenases from other microbial species translates into catalytic activity of CYP154C1 toward both 12- and 14-membered ring macrolactones in vitro.

Protein side-chain motion and hydration in proton-transfer pathways. Results for cytochrome p450cam

Proton-transfer reactions form an integral part of bioenergetics and enzymatic catalysis. The identification of proton-conducting pathways inside a protein is a key to understanding the mechanisms of biomolecular proton transfer. Proton pathways are modeled here as hydrogen bonded networks of proton-conducting groups, including proton-exchanging groups of amino acid side chains and bound water molecules. We focus on the identification of potential proton-conducting pathways inside a protein of known structure. However, consideration of the static structure alone is often not sufficient to detect suitable proton-transfer paths, leading, for example, from the protein surface to the active site buried inside the protein. We include dynamic fluctuations of amino acid side chains and water molecules into our analysis. To illustrate the method, proton transfer into the active site of cytochrome P450cam is studied. The cooperative rotation of amino acids and motion of water molecules are found to connect the protein surface to the molecular oxygen. Our observations emphasize the intrinsic dynamical nature of proton pathways where critical connections in the network may be transiently provided by mobile groups.

The role of copper and protons in heme-copper oxidases: kinetic study of an engineered heme-copper center in myoglobin

To probe the role of copper and protons in heme-copper oxidase (HCO), we have performed kinetic studies on an engineered heme-copper center in sperm whale myoglobin (Leu-29 --> HisPhe-43 --> His, called Cu(B)Mb) that closely mimics the heme-copper center in HCO. In the absence of metal ions, the engineered Cu(B) center in Cu(B)Mb decreases the O(2) binding affinity of the heme. However, addition of Ag(I), a redox-inactive mimic of Cu(I), increases the O(2)-binding affinity. More importantly, copper ion in the Cu(B) center is essential for O(2) reduction, as no O(2) reduction can be observed in copper-free, Zn(II), or Ag(I) derivatives of Cu(B)Mb. Instead of producing a ferryl-heme as in HCO, the Cu(B)Mb generates verdoheme because the engineered Cu(B)Mb may lack a hydrogen bonding network that delivers protons to promote the heterolytic OO cleavage necessary for the formation of ferryl-heme. Reaction of oxidized Cu(B)Mb with H(2)O(2), a species equivalent in oxidation state to 2e(-), reduced O(2) but, possessing the extra protons, resulted in ferryl-heme formation, as in HCO. The results showed that the Cu(B) center plays a critical role in O(2) binding and reduction, and that proton delivery during the O(2) reduction is important to avoid heme degradation and to promote the HCO reaction.

Proton-coupled O-O activation on a redox platform bearing a hydrogen-bonding scaffold

Porphyrin architectures bearing a hydrogen-bonding scaffold have been synthesized. The H-bond pendant allows proton-coupled electron transfer (PCET) to be utilized as a vehicle for effecting catalytic O-O bond activation chemistry. Suzuki cross-coupling reactions provide a modular synthetic strategy for the attachment of porphyrins to a rigid xanthene or dibenzofuran pillar bearing the H-bond pendant. The resulting HPX (hanging porphyrin xanthene) and HPD (hanging porphyrin dibenzofuran) systems permit both the orientation and acid-base properties of the hanging H-bonding group to be controlled. Comparative reactivity studies for the catalase-like disproportionation of hydrogen peroxide and the epoxidation of olefins by the HPX and HPD platforms with acid and ester hanging groups reveal that the introduction of a proton-transfer network, properly oriented to a redox-active platform, can orchestrate catalytic O-O bond activation. For the catalase and epoxidation reaction types, a marked reactivity enhancement is observed for the xanthene-bridged platform appended with a pendant carboxylic acid group, establishing that this approach can yield superior catalysts to analogues that do not control both proton and electron inventories.

Molecular recognition in (+)-alpha-pinene oxidation by cytochrome P450cam

Oxygenated derivatives of the monoterpene (+)-alpha-pinene are found in plant essential oils and used as fragrances and flavorings. (+)-alpha-Pinene is structurally related to (+)-camphor, the natural substrate of the heme monooxygenase cytochrome P450(cam) from Pseudomonas putida. The aim of the present work was to apply the current understanding of P450 substrate binding and catalysis to engineer P450(cam) for the selective oxidation of (+)-alpha-pinene. Consideration of the structures of (+)-camphor and (+)-alpha-pinene lead to active-site mutants containing combinations of the Y96F, F87A, F87L, F87W, and V247L mutations. All mutants showed greatly enhanced binding and rate of oxidation of (+)-alpha-pinene. Some mutants had tighter (+)-alpha-pinene binding than camphor binding by the wild-type. The most active was the Y96F/V247L mutant, with a (+)-alpha-pinene oxidation rate of 270 nmol (nmol of P450(cam))(-)(1) min(-)(1), which was 70% of the rate of camphor oxidation by wild-type P450(cam). Camphor is oxidized by wild-type P450(cam) exclusively to 5-exo-hydroxycamphor. If the gem dimethyl groups of (+)-alpha-pinene occupied similar positions to those found for camphor in the wild-type structure, (+)-cis-verbenol would be the dominant product. All P450(cam) enzymes studied gave (+)-cis-verbenol as the major product but with much reduced selectivity compared to camphor oxidation by the wild-type. (+)-Verbenone, (+)-myrtenol, and the (+)-alpha-pinene epoxides were among the minor products. The crystal structure of the Y96F/F87W/V247L mutant, the most selective of the P450(cam) mutants initially examined, was determined to provide further insight into P450(cam) substrate binding and catalysis. (+)-alpha-Pinene was bound in two orientations which were related by rotation of the molecule. One orientation was similar to that of camphor in the wild-type enzyme while the other was significantly different. Analysis of the enzyme/substrate contacts suggested rationalizations of the product distribution. In particular competition rather than cooperativity between the F87W and V247L mutations and substrate movement during catalysis were proposed to be major factors. The crystal structure lead to the introduction of the L244A mutation to increase the selectivity of pinene oxidation by further biasing the binding orientation toward that of camphor in the wild-type structure. The F87W/Y96F/L244A mutant gave 86% (+)-cis-verbenol and 5% (+)-verbenone. The Y96F/L244A/V247L mutant gave 55% (+)-cis-verbenol but interestingly also 32% (+)-verbenone, suggesting that it may be possible to engineer a P450(cam) mutant that could oxidize (+)-alpha-pinene directly to (+)-verbenone. Verbenol, verbenone, and myrtenol are naturally occurring plant fragrance and flavorings. The preparation of these compounds by selective enzymatic oxidation of (+)-alpha-pinene, which is readily available in large quantities, could have applications in synthesis. The results also show that the protein engineering of P450(cam) for high selectivity of substrate oxidation is more difficult than achieving high substrate turnover rates because of the subtle and dynamic nature of enzyme-substrate interactions.

Crystal structure of OxyB, a cytochrome P450 implicated in an oxidative phenol coupling reaction during vancomycin biosynthesis

Gene-inactivation studies point to the involvement of OxyB in catalyzing the first oxidative phenol coupling reaction during glycopeptide antibiotic biosynthesis. The oxyB gene has been cloned and sequenced from the vancomycin producer Amycolatopsis orientalis, and the hemoprotein has been produced in Escherichia coli, crystallized, and its structure determined to 1.7-A resolution. OxyB gave UV-visible spectra characteristic of a P450-like hemoprotein in the low spin ferric state. After reduction to the ferrous state by dithionite or by spinach ferredoxin and ferredoxin reductase, the CO-ligated form gave a 450-nm peak in a UV-difference spectrum. Addition of putative heptapeptide substrates to resting OxyB produced type I changes to the UV spectrum, but no turnover was observed in the presence of ferredoxin and ferredoxin reductase, showing that either the peptides or the reduction system, or both, are insufficient to support a full catalytic cycle. OxyB exhibits the typical P450-fold, with helix L containing the signature sequence FGHGXHXCLG and Cys(347) being the proximal axial thiolate ligand of the heme iron. The structural similarity of OxyB is highest to P450nor, P450terp, CYP119, and P450eryF. In OxyB, the F and G helices are rotated out of the active site compared with P450nor, resulting in a much more open active site, consistent with the larger size of the presumed heptapeptide substrate.

Formation and decay of hydroperoxo-ferric heme complex in horseradish peroxidase studied by cryoradiolysis

Using radiolytic reduction of the oxy-ferrous horseradish peroxidase (HRP) at 77 K, we observed the formation and decay of the putative intermediate, the hydroperoxo-ferric heme complex, often called "Compound 0." This intermediate is common for several different enzyme systems as the precursor of the Compound I (ferryl-oxo pi-cation radical) intermediate. EPR and UV-visible absorption spectra show that protonation of the primary intermediate of radiolytic reduction, the peroxo-ferric complex, to form the hydroperoxo-ferric complex is completed only after annealing at temperatures 150-180 K. After further annealing at 195-205 K, this complex directly transforms to ferric HRP without any observable intervening species. The lack of Compound I formation is explained by inability of the enzyme to deliver the second proton to the distal oxygen atom of hydroperoxide ligand, shown to be necessary for dioxygen bond heterolysis on the "oxidase pathway," which is non-physiological for HRP. Alternatively, the physiological substrate H2O2 brings both protons to the active site of HRP, and Compound I is subsequently formed via rearrangement of the proton from the proximal to the distal oxygen atom of the bound peroxide.

Mechanistic enzymology of oxygen activation by the cytochromes P450

The P450 cytochromes represent a universal class of heme-monooxygenases. The detailed mechanistic understanding of their oxidative prowess is a critical theme in the studies of metabolism of a wide range of organic compounds including xenobiotics. Integral to the O2 bond cleavage mechanism by P450 is the enzyme's concerted use of protein and solvent-mediated proton transfer events to transform reduced dioxygen to a species capable of oxidative chemistry. To this end, a wide range of kinetic, structural, and mutagenesis data has been accrued. A critical role of conserved acid-alcohol residues in the P450 distal pocket, as well as stabilized waters, enables the enzyme to catalyze effective monooxygenation chemistry. In this review, we discuss the detailed mechanism of P450 dioxygen scission utilizing the CYP101 hydroxylation of camphor as a model system. The application of low-temperature radiolytic techniques has enabled a structural and spectroscopic analysis of the nature of critical intermediate states in the reaction.

Soluble methane monooxygenase: activation of dioxygen and methane

The mechanisms by which soluble methane monooxygenase uses dioxygen to convert methane selectively to methanol have come into sharp focus. Diverse techniques have clarified subtle details about each step in the reaction, from binding and activating dioxygen, to hydroxylation of alkanes and other substrates, to the electron transfer events required to complete the catalytic cycle.

Oxidation and electronic state dependence of proton transfer in the enzymatic cycle of cytochrome P450eryF

Hydrogen bond networks, consisting of hydrogen bonded waters anchored by polar/acidic amino acid sidechains, are often present in the vicinity of the oxygen binding clefts of P450s. Density functional and quantum dynamics calculations of a O(2) binding cleft network model of cytochrome P450eryF(CYP107A1) indicate that such structural motifs facilitate ultrafast proton transfer from network waters to the dioxygen of the reduced oxyferrous species via a multiple proton translocation mechanism with barriers of 7-10 kcal/mol on its doublet ground state, and that the energies of the proton transfer reactant and constrained proton transfer products have an electronic and oxidation state dependence [J. Am. Chem. Soc. 124 (2002) 1430]. In the present study, the origin of the oxidation state dependence is shown to have its roots in differential proton affinities while the electronic state dependence of the reduced oxyferrous heme has its origins in subtle differences in network topologies near the transition state of the initial proton transfer event. Relaxed potential surface scans and unconstrained proton transfer product optimizations indicate that the proton transfer product in both the singlet oxyferrous heme and the reduced oxyferrous heme species in a quartet state are not viable stable (bound) states relative to the reactant form. While the proton affinity of H(3)O(+) is sufficient for it to protonate both the oxyferrous and the reduced oxyferrous heme species, hydrogen bond network stabilized water is only capable of protonating the reduced oxyferrous form. This interpretation is substantiated by study of the NO bound reduced ferrous heme of P450nor, which is isoelectronic with the oxyferrous heme and has a similar proton affinity. Density functional calculations on a more extensive O(2) binding cleft model support the multiple proton translocation mechanism of transfer but indicates that the significant negative charge density on the bound dioxygen of the reduced oxyferrous heme species, in its doublet ground state, polarizes the associated hydrogen bond network sufficiently so as to result in short, strong, low-barrier hydrogen bonds. The computed O-H-O bond distances are less than 2.55 A and have a near degeneracy of the proton transfer reactant and initial (sudden) proton transfer products. These low-barrier hydrogen bond features, in addition to the finding of a (zero point uncorrected) barrier of 1.3 kcal/mol, indicate that proton transfer from water to the distal oxygen should be rapid, facile and may not require large curvature tunneling as originally suggested by use of a smaller model. An initial assessment of protonation of the reduced oxyferrous heme distal oxygen by a model of 6-deoxyerythronolide B (6-DEB) indicates it to be low barrier (3.8 kcal/mol) and exothermic (-2.9 kcal/mol). The combined results indicate the plausibility of simultaneous diprotonation of the distal oxygen of the reduced oxyferrous heme, leading to O-O bond scission, using the combined water network and 6-DEB substrate protonation agents.

Specific and non-specific effects of potassium cations on substrate-protein interactions in cytochromes P450cam and P450lin

Substrate binding to cytochrome P450cam is generally considered to be a two-step process. The first step corresponds to the entrance of the substrate, camphor, into the heme pocket. The second step corresponds to a spin transition (low spin-->high spin) of the iron in the protein-substrate complex. This spin transition is related to the mobility of the substrate inside the active site [Biochim Biophys Acta 1338 (1997) 77]. Potassium cations (K(+)) have a specific effect on the spin equilibrium. This is generally attributed to the K(+) ion-induced conformational change of tyrosine 96, the hydroxyl group of which is hydrogen bonded to the keto group of camphor and results in optimum substrate orientation and reduced mobility of this substrate in the active site. In the present paper, we show that K(+) not only affects the substrate-Tyr 96 couple, but acts more globally since K(+) effects are also observed in the Tyr96Phe mutant as well as in complexes with camphor-analogues. Large compounds, that fit well in the heme pocket and bind with higher affinity than camphor, display high spin contents that are less dependent on the presence of K(+). In contrast, K(+) has a significant effect on the high spin content of substrate-cytochrome P450cam complexes with looser interactions. We conclude that large compounds with higher affinities than camphor have more van der Waals contacts with the active site residues. Their mobilities are then reduced and less dependent on the presence of K(+). In this study, we also explored, for comparison, the K(+) effect on the spin transition state of another member of the P450 superfamily, cytochrome P450lin. This effect is not as strong as those observed for cytochrome P450cam. Even though the spin equilibrium does not change dramatically in the presence of K(+) or Na(+), the value of the dissociation constant (K(d)) for linalool binding is significantly affected by ionic strength. Analysis of the thermodynamic parameters for the linalool binding strongly suggests that, similarly to our previous finding for cytochrome P450cam, electrostatic gates participate in the control of substrate access.

Thermophilic cytochrome P450 (CYP119) from Sulfolobus solfataricus: high resolution structure and functional properties

Crystal structures of a thermostable cytochrome P450 (CYP119) and a site-directed mutant, (Phe24Leu), from the acidothermophilic archaea Sulfolobus solfataricus were determined at 1.5-2.0 A resolution. We identify important crystallographic waters in the ferric heme pocket, observe protein conformational changes upon inhibitor binding, and detect a unique distribution of surface charge not found in other P450s. An analysis of factors contributing to thermostability of CYP119 of these high resolution structures shows an apparent increase in clustering of aromatic residues and optimum stacking. The contribution of aromatic stacking was investigated further with the mutant crystal structure and differential scanning calorimetry.

The 'push' effect of the thiolate ligand in cytochrome P450: a theoretical gauging

The 'push' effect of the thiolate ligand in cytochrome P450 is investigated using density functional calculations. Theory supports Dawson's postulate that the 'push' effect is crucial for the heterolytic O-O bond cleavage of ferric-peroxide, as well as for controlling the Fe(III)/Fe(II) redox process and gating the catalytic cycle. Two energetic factors that contribute to the 'push' effect are revealed. The dominant one is the field factor (DeltaE(field)=54-103 kcal/mol) that accounts for the classical electrostatic repulsion with the negative charge of thiolate. The smaller factor is a quantum mechanical effect (DeltaE(QM)(sigma)=39 kcal/mol, DeltaE(QM)(pi)=4 kcal/mol), which is associated with the sigma- and pi-donor capabilities of thiolate. The effects of ligand replacement, changes in hydrogen bonding and dielectric screening are discussed in term of these quantities. In an environment with a dielectric constant of 5.7, the total 'push' effect is reduced to 29-33 kcal/mol. Manifestations of the 'push' effect on other properties of thiolate enzymes are discussed.

Substrates modulate the rate-determining step for CO binding in cytochrome P450cam (CYP101). A high-pressure stopped-flow study

The high-pressure stopped-flow technique is applied to study the CO binding in cytochrome P450cam (P450cam) bound with homologous substrates (1R-camphor, camphane, norcamphor and norbornane) and in the substrate-free protein. The activation volume DeltaV # of the CO on-rate is positive for P450cam bound with substrates that do not contain methyl groups. The kon rate constant for these substrate complexes is in the order of 3 x 10(6) M(-1) x s(-1). In contrast, P450cam complexed with substrates carrying methyl groups show a negative activation volume and a low kon rate constant of approximately 3 x 10(4) M(-1) x s(-1). By relating kon and DeltaV # with values for the compressibility and the influx rate of water for the heme pocket of the substrate complexes it is concluded that the positive activation volume is indicative for a loosely bound substrate that guarantees a high solvent accessibility for the heme pocket and a very compressible active site. In addition, subconformers have been found for the substrate-free and camphane-bound protein which show different CO binding kinetics.