Primary intermediate in the reaction of mixed-valence cytochrome c oxidase with oxygen

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

A Water Molecule Residing in the Fea33+···CuB2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study

Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.

DFT Fea3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase

Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.

A Water Dimer Shift Activates a Proton Pumping Pathway in the PR → F Transition of ba3 Cytochrome c Oxidase

Broken-symmetry density functional calculations have been performed on the [Fe_a3⁴⁺,Cu_B²⁺] state of the dinuclear center (DNC) for the P_R → F part of the catalytic cycle of ba₃ cytochrome c oxidase (CcO) from Thermus thermophilus (Tt), using the OLYP-D3-BJ functional. The calculations show that the movement of the H₂O molecules in the DNC affects the pK_a values of the residue side chains of Tyr237 and His376⁺, which are crucial for proton transfer/pumping in ba₃ CcO from Tt. The calculated lowest energy structure of the DNC in the [Fe_a3⁴⁺,Cu_B²⁺] state (state F) is of the form Fe_a3⁴⁺═O^2-···Cu_B²⁺, in which the H₂O ligand that resulted from protonation of the OH^- ligand in the P_R state is dissociated from the Cu_B²⁺ site. The calculated Fe_a3⁴⁺═O^2- distance in F (1.68 Å) is 0.03 Å longer than that in P_R (1.65 Å), which can explain the different Fe_a3⁴⁺═O^2- stretching modes in P (804 cm^-1) and F (785 cm^-1) identified by resonance Raman experiments. In this F state, the Cu_B²⁺···O^2- (ferryl-oxygen) distance is only around 2.4 Å. Hence, the subsequent O_H state [Fe_a3³⁺-OH^--Cu_B²⁺] with a μ-hydroxo bridge can be easily formed, as shown by our calculations.

A broken-symmetry density functional study of structures, energies, and protonation states along the catalytic O-O bond cleavage pathway in ba3 cytochrome c oxidase from Thermus thermophilus

Broken-symmetry density functional calculations have been performed on the [Fea3, CuB] dinuclear center (DNC) of ba3 cytochrome c oxidase from Thermus thermophilus in the states of [Fea3(3+)-(HO2)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], using both PW91-D3 and OLYP-D3 functionals. Tyr237 is a special tyrosine cross-linked to His233, a ligand of CuB. The calculations have shown that the DNC in these states strongly favors the protonation of His376, which is above propionate-A, but not of the carboxylate group of propionate-A. The energies of the structures obtained by constrained geometry optimizations along the O-O bond cleavage pathway between [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] have also been calculated. The transition of [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] → [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] shows a very small barrier, which is less than 3.0/2.0 kcal mol(-1) in PW91-D3/OLYP-D3 calculations. The protonation state of His376 does not affect this O-O cleavage barrier. The rate limiting step of the transition from state A (in which O2 binds to Fea3(2+)) to state PM ([Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], where the O-O bond is cleaved) in the catalytic cycle is, therefore, the proton transfer originating from Tyr237 to O-O to form the hydroperoxo [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] state. The importance of His376 in proton uptake and the function of propionate-A/neutral-Asp372 as a gate to prevent the proton from back-flowing to the DNC are also shown.

Proton translocation in cytochrome c oxidase: insights from proton exchange kinetics and vibrational spectroscopy

Cytochrome c oxidase is the terminal enzyme in the electron transfer chain. It reduces oxygen to water and harnesses the released energy to translocate protons across the inner mitochondrial membrane. The mechanism by which the oxygen chemistry is coupled to proton translocation is not yet resolved owing to the difficulty of monitoring dynamic proton transfer events. Here we summarize several postulated mechanisms for proton translocation, which have been supported by a variety of vibrational spectroscopic studies. We recently proposed a proton translocation model involving proton accessibility to the regions near the propionate groups of the heme a and heme a3 redox centers of the enzyme based by hydrogen/deuterium (H/D) exchange Raman scattering studies (Egawa et al., PLoS ONE 2013). To advance our understanding of this model and to refine the proton accessibility to the hemes, the H/D exchange dependence of the heme propionate group vibrational modes on temperature and pH was measured. The H/D exchange detected at the propionate groups of heme a3 takes place within a few seconds under all conditions. In contrast, that detected at the heme a propionates occurs in the oxidized but not the reduced enzyme and the H/D exchange is pH-dependent with a pKa of ~8.0 (faster at high pH). Analysis of the thermodynamic parameters revealed that, as the pH is varied, entropy/enthalpy compensation held the free energy of activation in a narrow range. The redox dependence of the possible proton pathways to the heme groups is discussed. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

Electronic structure of a low-spin heme/Cu peroxide complex: spin-state and spin-topology contributions to reactivity

This study details the electronic structure of the heme–peroxo–copper adduct {[(F8)Fe(DCHIm)]-O2-[Cu(AN)]}+ (LS(AN)) in which O2(2–) bridges the metals in a μ-1,2 or “end-on” configuration. LS(AN) is generated by addition of coordinating base to the parent complex {[(F8)Fe]-O2-[Cu(AN)]}+ (HS(AN)) in which the O2(2–) bridges the metals in an μ-η2:η2 or “side-on” mode. In addition to the structural change of the O2(2–) bridging geometry, coordination of the base changes the spin state of the heme fragment (from S = 5/2 in HS(AN) to S = 1/2 in LS(AN)) that results in an antiferromagnetically coupled diamagnetic ground state in LS(AN). The strong ligand field of the porphyrin modulates the high-spin to low-spin effect on Fe–peroxo bonding relative to nonheme complexes, which is important in the O–O bond cleavage process. On the basis of DFT calculations, the ground state of LS(AN) is dependent on the Fe–O–O–Cu dihedral angle, wherein acute angles (<~150°) yield an antiferromagnetically coupled electronic structure while more obtuse angles yield a ferromagnetic ground state. LS(AN) is diamagnetic and thus has an antiferromagnetically coupled ground state with a calculated Fe–O–O–Cu dihedral angle of 137°. The nature of the bonding in LS(AN) and the frontier molecular orbitals which lead to this magneto-structural correlation provide insight into possible spin topology contributions to O–O bond cleavage by cytochrome c oxidase.

Critical structural role of R481 in cytochrome c oxidase from Rhodobacter sphaeroides

The R481 residue in cytochrome c oxidase from Rhodobacter sphaeroides forms hydrogen bonds with the propionate groups of both heme a and heme a(3). It has been postulated that R481 is the proton loading site in the proton exit pathway essential for proton translocation. A recent functional study showed that the mutations of R481 to His, Leu and Gln cause the reduction of the activity to approximately 5-18% of the native level, and the absence of proton pumping in R481Q but retention of approximately 40% efficiency in R481H and R481L (H.J. Lee, L. Ojemyr, A. Vakkasoglu, P. Brzezinski and R. B. Gennis, manuscript submitted). To decipher the molecular mechanism underlying the perturbed functionalities, we have used resonance Raman spectroscopy to examine the structural properties of the three mutants. The data show that the frequencies of the formyl CO stretching modes of both the heme a and a(3) in the mutants are characteristic of formyl groups exposed to an aqueous environment, indicating that the mutations disrupt the native H-bonding interaction between the formyl group of heme a and R52, as well as the hydrophobic environment surrounding the formyl group of heme a(3). In addition to the change in the environments of heme a and a(3), the Raman data show that the mutations induce a partial conversion of the heme a(3) from a high-spin to a low-spin state, suggesting that the mutations are associated with the rearrangement of the Cu(B)-heme a(3) binuclear center. The Raman results reported here demonstrate that R481 plays a critical role in supporting efficient proton pumping, by holding the heme groups in a proper environment.

Cytochrome c oxidase: chemistry of a molecular machine

The plethora of proposed chemical models attempting to explain the proton pumping reactions catalyzed by the CcO complex, especially the number of recent models, makes it clear that the problem is far from solved. Although we have not discussed all of the models proposed to date, we have described some of the more detailed models in order to illustrate the theoretical concepts introduced at the beginning of this section on proton pumping as well as to illustrate the rich possibilities available for effecting proton pumping. It is clear that proton pumping is effected by conformational changes induced by oxidation/reduction of the various redox centers in the CcO complex. It is for this reason that the CcO complex is called a redox-linked proton pump. The conformational changes of the proton pump cycle are usually envisioned to be some sort of ligand-exchange reaction arising from unstable geometries upon oxidation/reduction of the various redox centers. However, simple geometrical rearrangements, as in the Babcock and Mitchell models are also possible. In any model, however, hydrogen bonds must be broken and reformed due to conformational changes that result from oxidation/reduction of the linkage site during enzyme turnover. Perhaps the most important point emphasized in this discussion, however, is the fact that proton pumping is a directed process and it is electron and proton gating mechanisms that drive the proton pump cycle in the forward direction. Since many of the models discussed above lack effective electron and/or proton gating, it is clear that the major difficulty in developing a viable chemical model is not formulating a cyclic set of protein conformational changes effecting proton pumping (redox linkage) but rather constructing the model with a set of physical constraints so that the proposed cycle proceeds efficiently as postulated. In our discussion of these models, we have not been too concerned about which electron of the catalytic cycle was entering the site of linkage, but merely whether an ET to the binuclear center played a role. However, redox linkage only occurs if ET to the activated binuclear center is coupled to the proton pump. Since all of the models of proton pumping presented here, with the exception of the Rousseau expanded model and the Wikström model, have a maximum stoichiometry of 1 H+/e-, they inadequately explain the 2 H+/e- ratio for the third and fourth electrons of the dioxygen reduction cycle (see Section V.B). One way of interpreting this shortfall of protons is that the remaining protons are pumped by an as yet undefined indirectly coupled mechanism. In this scenario, the site of linkage could be coupled to the pumping of one proton in a direct fashion and one proton in an indirect fashion for a given electron. For a long time, it was assumed that at least some elements of such an indirect mechanism reside in subunit III. While recent evidence argues against the involvement of subunit III in the proton pump, subunit III may still participate in a regulatory and/or structural capacity (Section II.E). Attention has now focused on subunits I and II in the search for residues intimately involved in the proton pump mechanism and/or as part of a proton channel. In particular, the role of some of the highly conserved residues of helix VIII of subunit I are currently being studied by site directed mutagenesis. In our opinion, any model that invokes heme alpha 3 or CuB as the site of linkage must propose a very effective means by which the presumedly fast uncoupling ET to the dioxygen intermediates is prevented. It is difficult to imagine that ET over the short distance from heme alpha 3 or CuB to the dioxygen intermediate requires more than 1 ns. In addition, we expect the conformational changes of the proton pump to require much more than 1 ns (see Section V.B).

Structures and proton-pumping strategies of mitochondrial respiratory enzymes

Enzymes of the mitochondrial respiratory chain serve as proton pumps, using the energy made available from electron transfer reactions to transport protons across the inner mitochondrial membrane and create an electrochemical gradient used for the production of ATP. The ATP synthase enzyme is reversible and can also serve as a proton pump by coupling ATP hydrolysis to proton translocation. Each of the respiratory enzymes uses a different strategy for performing proton pumping. In this work, the strategies are described and the structural bases for the action of these proteins are discussed in light of recent crystal structures of several respiratory enzymes. The mechanisms and efficiency of proton translocation are also analyzed in terms of the thermodynamics of the substrate transformations catalyzed by these enzymes.

Time dependence of the catalytic intermediates in cytochrome c oxidase

Cytochrome c oxidase, the terminal enzyme in the electron transfer chain, catalyzes the reduction of oxygen to water in a multiple step process by utilizing four electrons from cytochrome c. To study the reaction mechanism, the resonance Raman spectra of the intermediate states were measured during single turnover of the enzyme after catalytic initiation by photolysis of CO from the fully reduced CO-bound enzyme. By measuring the change in intensity of lines associated with heme a, the electron transfer steps were determined and found to be biphasic with apparent rate constants of approximately 40 x 10(3) s(-1) and approximately 1 x 10(3) s(-1). The time dependence for the oxidation of heme a and for the measured formation and decay of the oxy, the ferryl ("F"), and the hydroxy intermediates could be simulated by a simple reaction scheme. In this scheme, the presence of the "peroxy" ("P") intermediate does not build up a sufficient population to be detected because its decay rate is too fast in buffered H(2)O at neutral pH. A comparison of the change in the spin equilibrium with the formation of the hydroxy intermediate demonstrates that this intermediate is high spin. We also confirm the presence of an oxygen isotope-sensitive line at 355 cm(-1), detectable in the spectrum from 130 to 980 micros, coincident with the presence of the F intermediate.

Crystal structure of bovine heart cytochrome c oxidase at 2.8 A resolution

Thirteen different polypeptide subunits, each in one copy, five phosphatidyl ethanolamines and three phosphatidyl glycerols, two hemes A, three Cu ions, one Mg ion, and one Zn ion are detectable in the crystal structure of bovine heart cytochrome c oxidase in the fully oxidized form at 2.8 A resolution. A propionate of hems a, a peptide unit (-CO-NH-), and an imidazole bound to CuA are hydrogen-bonded sequentially, giving a facile electron transfer path from CUA to heme a. The O2 binding and reduction site, heme a3, is 4.7 A apart from CuB. Two possible proton transfer paths from the matrix side to the cytosolic side are located in subunit I, including hydrogen bonds and internal cavities likely to contain randomly oriented water molecules. Neither path includes the O2 reduction site. The O2 reduction site has a proton transfer path from the matrix side possibly for protons for producing water. The coordination geometry of CuB and the location of Tyr244 in subunit I at the end of the scalar proton path suggests a hydroperoxo species as the two electron reduced intermediate in the O2 reduction process.

Microsecond generation of oxygen-bound cytochrome c oxidase by rapid solution mixing

Current understanding of the oxygen reduction and proton translocation processes in cytochrome c oxidase is largely derived from the data obtained by a nonphysiological method for initiating the catalytic reaction: photolyzing carbon monoxide (CO) from the CO-inhibited enzyme in the presence of oxygen (O2). However, considerable evidence suggests that the use of CO introduces artifacts into the reaction mechanism. We have therefore developed a rapid solution mixer with a mixing time of 20 microseconds to study the catalytic reaction by directly mixing the enzyme with O2 without using CO. Unexpectedly, the resonance Raman scattering detected for the first 120 microseconds after the mixing show that the CO influences neither the structure of the primary oxy-intermediate, its rate of decay, nor the rate of oxidation of cytochrome a. This implies that CO has an effect on the later stages of the catalytic process, which may involve the proton translocation steps, and calls for the re-examination of the catalytic process by using the direct mixing method. In addition, these results demonstrate the feasibility of using the rapid mixing device for the study of biological reactions in the microsecond time domain.

The cytochrome oxidase superfamily of redox-driven proton pumps

Most respiratory oxidases of eukaryotic and prokaryotic organisms are members of a superfamily of enzymes that couple the redox energy available from the reduction of molecular oxygen to the mechanism of pumping protons across the membrane. The recent applications of site-directed mutagenesis and of a variety of spectroscopic techniques have allowed major advances in our understanding of the structure and function of these proteins.

Protein-heme interactions in hemoglobin from the mollusc Scapharca inaequivalvis: evidence from resonance Raman scattering

Resonance Raman spectra of the Scapharca inaequivalvis homodimeric hemoglobin (HbI) have been measured for the ligand-bound and ligand-free ferrous forms of the protein. In the deoxy derivative, the iron-histidine (Fe-His) stretching mode, proposed as a marker of the oxygen affinity and a conduit linking the hemes to the subunit interface, gives rise to a Raman peak centered at 203 cm-1, an unusually low frequency compared to that reported for other hemoglobins and myoglobins. In the CO-bound derivative, three isotope-sensitive lines at 517, 583, and 1945 cm-1 have been assigned to the Fe-CO stretching, Fe-C-O bending, and C-O stretching modes, respectively. From the frequencies of these modes and from their relative intensities, the Fe-C-O geometry appears to be tilted from axial coordination and shows a bending angle which has been estimated to be about 171 +/- 5 degrees. For the oxygen derivative, only one isotope-sensitive peak has been detected at 570 cm-1, in line with the values found for myoglobin and other hemoglobins. Resonance Raman spectra of HbI modified with p-(chloromercuri)benzoate (PMB) at Cys92 have been measured in parallel with those of the native protein. Despite the large increase in oxygen affinity produced by the PMB modification, the frequency of the Fe-His stretching mode is unshifted in the deoxy derivative.(ABSTRACT TRUNCATED AT 250 WORDS)

Probing heart cytochrome c oxidase structure and function by infrared spectroscopy

IR spectra directly probe specific vibrators in bovine heart cytochrome c oxidase, yielding quantitative as well as qualitative information on structures and reactions at these vibrators. C-O IR spectra reveal that CO binds to Fe2+ a3 as two conformers each in isolated immobile environments sensitive to Fea and/or CuA oxidation state but remarkably insensitive to pH, medium, anesthetics, and other factors that affect activity. C-N IR spectra reveal that the one CN- that binds to fully and partially oxidized enzyme can be in three different structures. These structures vary in relative amounts with redox level, thereby reflecting dynamic electron exchange among Fea, CuA, and CuB with associated changes in protein conformation of likely significance in O2 reduction and H(+)-pumping. Azide IR spectra also reflect redox-dependent long-range effects. The amide I IR bands, due to C-O vibrators of peptide linkages and composed of multiple bands derived from different secondary structures, reveal high levels of alpha-helix (approximately 60%) and subtle changes with redox level and exposure to anesthetics. N2O IR spectra reveal that these anesthetic molecules at clinically relevant levels occupy three sites of different polarity within the enzyme as the enzyme is reversibly, but only partially, inhibited.

Discrete steps in dioxygen activation--the cytochrome oxidase/O2 reaction

The kinetic constraints that are imposed on cytochrome oxidase in its dual function as the terminal oxidant in the respiratory process and as a redox-linked proton pump provide a unique opportunity to investigate the molecular details of biological O2 activation. By using flow/flash techniques, it is possible to visualize individual steps in the O2-binding and reduction process, and results from a number of spectroscopic investigations on the oxidation of reduced cytochrome oxidase by O2 are now available. In this article, we use these results to synthesize a reaction mechanism for O2 activation in the enzyme and to simulate time-concentration profiles for a number of intermediates that have been observed experimentally. Kinetic manifestation of the consequences of coupling exergonic electron transfer to endergonic proton translocation emerge from this analysis. Energetic efficiency in this process apparently requires that potentially toxic intermediate oxidation states of dioxygen accumulate to substantial concentration during the reduction reaction.

Proton translocation in cytochrome c oxidase: redox linkage through proximal ligand exchange on cytochrome a3

An analysis of resonance Raman scattering data from CO-bound cytochrome c oxidase and from the photodissociated enzyme indicates that histidine may not be coordinated to the iron atom of cytochrome a3 in the CO-bound form of the enzyme. Instead, the data suggest that either a water molecule or a different amino acid residue occupies the proximal ligand position. From these data, it is postulated that ligand exchange on cytochrome a3 can occur under physiological conditions. Studies of mutant hemoglobins have demonstrated that tyrosinate binds preferentially to histidine in the ferric forms of the proteins. In cytochrome c oxidase tyrosine residues are located near the histidine residues recently implicated in coordination to cytochrome a3 (Shapleigh et al., 1992; Hosler et al., this volume). Expanding on these concepts, we propose a model for proton translocation at the O2-binding site based on proximal ligand exchange between tyrosine and histidine on cytochrome a3. The pumping steps take place at the level of the peroxy intermediate and at the level of the ferryl intermediate in the catalytic cycle and are thereby consistent with the recent results of Wilkstrom (1989) who found that proton pumping occurs only at these two steps. It is shown that the model may be readily extended to account for the pumping of two protons at each of the steps.

Oxygen activation and the conservation of energy in cell respiration

Many of the membrane-associated oxidases that catalyse respiratory reduction of O2 to water simultaneously couple this exergonic reaction to the translocation of protons across the inner mitochondrial membrane, or the cell membrane in prokaryotes, a process by which metabolic energy is conserved for subsequent synthesis of ATP. The molecular mechanism of O2 reduction and its linkage to H+ translocation are now emerging. The bimetallic haem iron-copper reaction centre in this family of enzymes is the critical structure for catalysis of both these processes.

Ferryl and hydroxy intermediates in the reaction of oxygen with reduced cytochrome c oxidase

Cytochrome c oxidase catalyses the 4-electron reduction of dioxygen to water and translocates protons vectorially across the inner mitochondrial membrane. Proposed reaction pathways for the catalytic cycle of the O2 reduction are difficult to verify without knowing the structures of the intermediates, but we now have such information for the catalytic intermediates in the first steps of the reaction of O2 with cytochrome c oxidase from resonance Raman spectroscopy, a technique that enables iron-ligand stretching modes to be identified. Here we report on two more key intermediates: a ferryl-oxo (Fe4 = O2-) and a ferric-hydroxy (Fe3+--OH-) intermediate at the level of 3- and 4-electron reduction, respectively. We identified these intermediates by their characteristic iron-oxygen stretching frequencies (786 cm-1 for Fe4+ = O2-, and 450 cm-1 for Fe3+ -- OH-) and oxygen and deuterium isotope shifts. The oxo atom in the ferryl intermediate is hydrogen-bonded and the iron-oxygen bond in the hydroxy intermediate is anomalously weak. With the identification of the primary, ferryl and hydroxy intermediates, the predominant structures at almost all stages of O2 reduction are now known and the catalytic pathway can be described with more certainty.

Cytochrome c oxidase: decay of the primary oxygen intermediate involves direct electron transfer from cytochrome a

The decay of the primary intermediate generated in the reaction of oxygen with cytochrome c oxidase is nearly one order of magnitude faster in the fully reduced form of the enzyme than it is in the mixed valence form. To account for this observation, we propose a model describing the early molecular events in the reaction. In this model the decay of the primary Fe-O2 intermediate in the fully reduced enzyme is a consequence of direct electron transfer from cytochrome a. To test the model we measured the time dependence of the oxidation of cytochrome a by monitoring the resonance Raman scattering intensity of its vibrational modes. A rapid oxidation of cytochrome a was detected that quantitatively agrees with the model. These results indicate that the mechanism of oxygen reduction and proposed frameworks for proton translocation must be reexamined.