The protonation state of the cross-linked tyrosine during the catalytic cycle of cytochrome c oxidase

Examination of 'high-energy' metastable state of the oxidized (OH) bovine cytochrome c oxidase: Proton uptake and reaction with H2O2

Reoxidized cytochrome c oxidase appears to be in a 'high-energy' metastable state (OH) in which part of the energy released in the redox reactions is stored. The OH is supposed to relax to the resting 'as purified' oxidized state (O) in a time exceeding 200 ms. The catalytic heme a3-CuB center of these two forms should differ in a protonation and ligation state and the transition of OH-to-O is suggested to be associated with a proton transfer into this center. Employing a stopped-flow and UV-Vis absorption spectroscopy we investigated a proton uptake during the predicted relaxation of OH. It is shown, using a pH indicator phenol red, that from the time when the oxidation of the fully reduced CcO is completed (∼25 ms) up to ∼10 min, there is no uptake of a proton from the external medium (pH 7.8). Moreover, interactions of the assumed OH, generated 100 ms after oxidation of the fully reduced CcO, and the O with H2O2 (1 mM), result in the formation of two ferryl intermediates of the catalytic center, P and F, with very similar kinetics and the amounts of the formed ferryl states in both cases. These results implicate that the relaxation time of the catalytic center during the OH-to-O transition is either shorter than 100 ms or there is no difference in the structure of heme a3-CuB center of these two forms.

Structural insights into functional properties of the oxidized form of cytochrome c oxidase

Cytochrome c oxidase (CcO) is an essential enzyme in mitochondrial and bacterial respiration. It catalyzes the four-electron reduction of molecular oxygen to water and harnesses the chemical energy to translocate four protons across biological membranes. The turnover of the CcO reaction involves an oxidative phase, in which the reduced enzyme (R) is oxidized to the metastable OH state, and a reductive phase, in which OH is reduced back to the R state. During each phase, two protons are translocated across the membrane. However, if OH is allowed to relax to the resting oxidized state (O), a redox equivalent to OH, its subsequent reduction to R is incapable of driving proton translocation. Here, with resonance Raman spectroscopy and serial femtosecond X-ray crystallography (SFX), we show that the heme a3 iron and CuB in the active site of the O state, like those in the OH state, are coordinated by a hydroxide ion and a water molecule, respectively. However, Y244, critical for the oxygen reduction chemistry, is in the neutral protonated form, which distinguishes O from OH, where Y244 is in the deprotonated tyrosinate form. These structural characteristics of O provide insights into the proton translocation mechanism of CcO.

Structural basis for functional properties of cytochrome c oxidase

Cytochrome c oxidase (CcO) is an essential enzyme in mitochondrial and bacterial respiration. It catalyzes the four-electron reduction of molecular oxygen to water and harnesses the chemical energy to translocate four protons across biological membranes, thereby establishing the proton gradient required for ATP synthesis1. The full turnover of the CcO reaction involves an oxidative phase, in which the reduced enzyme (R) is oxidized by molecular oxygen to the metastable oxidized OH state, and a reductive phase, in which OH is reduced back to the R state. During each of the two phases, two protons are translocated across the membranes2. However, if OH is allowed to relax to the resting oxidized state (O), a redox equivalent to OH, its subsequent reduction to R is incapable of driving proton translocation2,3. How the O state structurally differs from OH remains an enigma in modern bioenergetics. Here, with resonance Raman spectroscopy and serial femtosecond X-ray crystallography (SFX)4, we show that the heme a3 iron and CuB in the active site of the O state, like those in the OH state5,6, are coordinated by a hydroxide ion and a water molecule, respectively. However, Y244, a residue covalently linked to one of the three CuB ligands and critical for the oxygen reduction chemistry, is in the neutral protonated form, which distinguishes O from OH, where Y244 is in the deprotonated tyrosinate form. These structural characteristics of O provide new insights into the proton translocation mechanism of CcO.

The three-spin intermediate at the O-O cleavage and proton-pumping junction in heme-Cu oxidases

[Figure: see text].

A Resonance Raman Marker Band Characterizes the Slow and Fast Form of Cytochrome c Oxidase

Cytochrome c oxidase (CcO) in its as-isolated form is known to exist in a slow and fast form, which differ drastically in their ability to bind oxygen and other ligands. While preparation methods have been established that yield either the fast or the slow form of the protein, the underlying structural differences have not been identified yet. Here, we have performed surface enhanced resonance Raman (SERR) spectroscopy of CcO immobilized on electrodes in both forms. SERR spectra obtained in resonance with the heme a₃ metal-to-ligand charge transfer (MLCT) transition at 650 nm displayed a sharp vibrational band at 748 or 750 cm^-1 when the protein was in its slow or fast form, respectively. DFT calculations identified the band as a mode of the His-419 ligand that is sensitive to the oxygen ligand and the protonation state of Tyr-288 within the binuclear complex. Potential-dependent SERR spectroscopy showed a redox-induced change of this band around 525 mV versus Ag/AgCl exclusively for the fast form, which coincides with the redox potential of the Tyr-O/Tyr-O^- transition. Our data points to a peroxide ligand in the resting state of CcO for both forms. The observed frequencies and redox sensitivities of the Raman marker band suggest that a radical Tyr-288 is present in the fast form and a protonated Tyr-288 in the slow form.

Protonation Dynamics in the K-Channel of Cytochrome c Oxidase Estimated from Molecular Dynamics Simulations

Proton transfer reactions are one of the most fundamental processes in biochemistry. We present a simplistic approach for estimating proton transfer probabilities in a membrane protein, cytochrome c oxidase. We combine short molecular dynamics simulations at discrete protonation states with a Monte Carlo approach to exchange between those states. Requesting for a proton transfer the existence of a hydrogen-bonded connection between the two source and target residues of the exchange, restricts the acceptance of transfers to only those in which a proton-relay is possible. Together with an analysis of the hydrogen-bonded connectivity in one of the proton-conducting channels of cytochrome c oxidase, this approach gives insight into the protonation dynamics of the hydrogen-bonded networks. The connectivity and directionality of the networks are coupled to the conformation of an important protein residue in the channel, K362, rendering proton transfer in the entire channel feasible in only one of the two major conformations. Proton transport in the channel can thus be regulated by K362 not only through its possible role as a proton carrier itself, but also by allowing or preventing proton transport via water residues.

Probing the Proton-Loading Site of Cytochrome C Oxidase Using Time-Resolved Fourier Transform Infrared Spectroscopy

Crystal structure analyses at atomic resolution and FTIR spectroscopic studies of cytochrome c oxidase have yet not revealed protonation or deprotonation of key sites of proton transfer in a time-resolved mode. Here, a sensitive technique to detect protolytic transitions is employed. In this work, probing a proton-loading site of cytochrome c oxidase from Paracoccus denitrificans with time-resolved Fourier transform infrared spectroscopy is presented for the first time. For this purpose, variants with single-site mutations of N131V, D124N, and E278Q, the key residues in the D-channel, were studied. The reaction of mutated CcO enzymes with oxygen was monitored and analyzed. Seven infrared bands in the “fast” kinetic spectra were found based on the following three requirements: (1) they are present in the “fast” phases of N131V and D124N mutants, (2) they have reciprocal counterparts in the “slow” kinetic spectra in these mutants, and (3) they are absent in “fast” kinetic spectra of the E278Q mutant. Moreover, the double-difference spectra between the first two mutants and E278Q revealed more IR bands that may belong to the proton-loading site protolytic transitions. From these results, it is assumed that several polar residues and/or water molecule cluster(s) share a proton as a proton-loading site. This site can be propionate itself (holding only a fraction of H+), His403, and/or water cluster(s).

Comparison Between O and OH Intermediates of Cytochrome c Oxidase Studied by FTIR Spectroscopy

Cytochrome c oxidase is terminal enzyme in the respiratory chain of mitochondria and many aerobic bacteria. It catalyzes reduction of oxygen to water. During its catalysis, CcO proceeds through several quite stable intermediates (R, A, PR/M, O/OH, E/EH). This work is concentrated on the elucidation of the differences between structures of oxidized intermediates O and O H in different CcO variants and at different pH values. Oxidized intermediates of wild type and mutated CcO from Paracoccus denitrificans were studied by means of static and time-resolved Fourier-transform infrared spectroscopy in acidic and alkaline conditions in the infrared region 1800-1000 cm-1. No reasonable differences were found between all variants in these conditions, and in this spectral region. This finding means that the binuclear center of oxygen reduction keeps a very similar structure and holds the same ligands in the studied conditions. The further investigation in search of differences should be performed in the 4000-2000 cm-1 IR region where water ligands absorb.

The structure of the oxidized state of cytochrome c oxidase - experiments and theory compared

Cytochrome c oxidase (CcO), the terminal enzyme in the respiratory chain, reduces molecular oxygen to water. Experimental data on the midpoint potentials of the heme iron/copper active site cofactors do not match the overall reaction energetics, and are also in conflict with the observed efficiency of energy conservation in CcO. Therefore it has been postulated that the ferric/cupric intermediate (the oxidized state) exists in two forms. One form, labelled O_H, is presumably involved during catalytic turnover, and should have a high Cu_B midpoint potential due to a metastable high energy structure. When no more electrons are supplied, the O_H state supposedly relaxes to the resting form, labelled O, with a lower energy and a lower midpoint potential. It has been suggested that there is a pure geometrical difference between the O_H and O states, obtained by moving a water molecule inside the active site. It is shown here that the difference between the two forms of the oxidized state must be of a more chemical nature. The reason is that all types of geometrically relaxed structures of the oxidized intermediate have similar energies, all with a high proton coupled reduction potential in accordance with the postulated O_H state. One hypothesized chemical modification of the O_H state is the transfer of an extra proton, possibly internal, into the active site. Such a protonated state has several properties that agree with experimental data on the relaxed oxidized state, including a decreased midpoint potential.

Geometric and Electronic Structure Contributions to O-O Cleavage and the Resultant Intermediate Generated in Heme-Copper Oxidases

This study investigates the mechanism of O-O bond cleavage in heme-copper oxidase (HCO) enzymes, combining experimental and computational insights from enzyme intermediates and synthetic models. It is determined that HCOs undergo a proton-initiated O-O cleavage mechanism where a single water molecule in the active site enables proton transfer (PT) from the cross-linked tyrosine to a peroxo ligand bridging the heme Fe^III and Cu^II, and multiple H-bonding interactions lower the tyrosine p K_a. Due to sterics within the active site, the proton must either transfer initially to the O(Fe) (a high-energy intermediate), or from another residue over a ∼10 Å distance to reach the O(Cu) atom directly. While the distance between the H⁺ donor (Tyr) and acceptor (O(Cu)) results in a barrier to PT, this separation is critical for the low barrier to O-O cleavage as it enhances backbonding from Fe into the O₂^2- σ* orbital. Thus, PT from Tyr precedes O-O elongation and is rate-limiting, consistent with available kinetic data. The electron transfers from tyrosinate after the barrier via a superexchange pathway provided by the cross-link, generating intermediate P_M. P_M is evaluated using available experimental data. The geometric structure contains an Fe^IV═O that is H-bonded to the Cu^II-OH. The electronic structure is a singlet, where the Fe^IV and Cu^II are antiferromagnetically coupled through the H-bond between the oxo(Fe) and hydroxo(Cu) ligands, while the Cu^II and Tyr^• are ferromagnetically coupled due their delocalization into orthogonal magnetic orbitals on the cross-linked His residue. These findings provide critical insights into the mechanism of efficient O₂ reduction in HCOs, and the nature of the P_M intermediate that couples this reaction to proton pumping.

Snapshot of an oxygen intermediate in the catalytic reaction of cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces dioxygen to water and harnesses the chemical energy to drive proton translocation across the inner mitochondrial membrane by an unresolved mechanism. By using time-resolved serial femtosecond crystallography, we identified a key oxygen intermediate of bovine CcO. It is assigned to the P_R-intermediate, which is characterized by specific redox states of the metal centers and a distinct protein conformation. The heme a ₃ iron atom is in a ferryl (Fe⁴⁺ = O^2-) configuration, and heme a and Cu_B are oxidized while Cu_A is reduced. A Helix-X segment is poised in an open conformational state; the heme a farnesyl sidechain is H-bonded to S382, and loop-I-II adopts a distinct structure. These data offer insights into the mechanism by which the oxygen chemistry is coupled to unidirectional proton translocation.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.

Dewetting transitions coupled to K-channel activation in cytochrome c oxidase

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane.

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane. CcO employs the so-called D- and K-channels for proton uptake, but the molecular mechanism for activation of the K-channel has remained elusive for decades. We show here by combining large-scale atomistic molecular simulations with graph-theoretical water network analysis, and hybrid quantum/classical (QM/MM) free energy calculations, that the K-channel is activated by formation of a reactive oxidized intermediate in the binuclear heme a₃/Cu_B active site. This state induces electrostatic, hydration, and conformational changes that lower the barrier for proton transfer along the K-channel by dewetting pathways that connect the D-channel with the active site. Our combined results reconcile previous experimental findings and indicate that water dynamics plays a decisive role in the proton pumping machinery in CcO.

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.

Mitochondrial cytochrome c oxidase: catalysis, coupling and controversies

Mitochondrial cytochrome c oxidase is a member of a diverse superfamily of haem-copper oxidases. Its mechanism of oxygen reduction is reviewed in terms of the cycle of catalytic intermediates and their likely chemical structures. This reaction cycle is coupled to the translocation of protons across the inner mitochondrial membrane in which it is located. The likely mechanism by which this occurs, derived in significant part from studies of bacterial homologues, is presented. These mechanisms of catalysis and coupling, together with current alternative proposals of underlying mechanisms, are critically reviewed.

Phenol-Induced O-O Bond Cleavage in a Low-Spin Heme-Peroxo-Copper Complex: Implications for O2 Reduction in Heme-Copper Oxidases

This study evaluates the reaction of a biomimetic heme-peroxo-copper complex, {[(DCHIm)(F8)FeIII]-(O22-)-[CuII(AN)]}+ (1), with a phenolic substrate, involving a net H-atom abstraction to cleave the bridging peroxo O-O bond that produces FeIV═O, CuII-OH, and phenoxyl radical moieties, analogous to the chemistry carried out in heme-copper oxidases (HCOs). A 3D potential energy surface generated for this reaction reveals two possible reaction pathways: one involves nearly complete proton transfer (PT) from the phenol to the peroxo ligand before the barrier; the other involves O-O homolysis, where the phenol remains H-bonding to the peroxo OCu in the transition state (TS) and transfers the H+ after the barrier. In both mechanisms, electron transfer (ET) from phenol occurs after the PT (and after the barrier); therefore, only the interaction with the H+ is involved in lowering the O-O cleavage barrier. The relative barriers depend on covalency (which governs ET from Fe), and therefore vary with DFT functional. However, as these mechanisms differ by the amount of PT at the TS, kinetic isotope experiments were conducted to determine which mechanism is active. It is found that the phenolic proton exhibits a secondary kinetic isotope effect, consistent with the calculations for the H-bonded O-O homolysis mechanism. The consequences of these findings are discussed in relation to O-O cleavage in HCOs, supporting a model in which a peroxo intermediate serves as the active H+ acceptor, and both the H+ and e- required for O-O cleavage derive from the cross-linked Tyr residue present at the active site.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286_I, Y288_I, and K362_I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.

Modeling the detailed kinetics of mitochondrial cytochrome c oxidase: Catalytic mechanism and nitric oxide inhibition

Cytochrome c oxidase (CcO) catalyzes the exothermic reduction of O₂ to H₂O by using electrons from cytochrome c, and hence plays a crucial role in ATP production. Although details on the enzyme structure and redox centers involved in O₂ reduction have been known, there still remains a considerable ambiguity on its mechanism of action, e.g., the number of sequential electrons donated to O₂ in each catalytic step, the sites of protonation and proton pumping, and nitric oxide (NO) inhibition mechanism. In this work, we developed a thermodynamically constrained mechanistic mathematical model for the catalytic action of CcO based on available kinetic data. The model considers a minimal number of redox centers on CcO and couples electron transfer and proton pumping driven by proton motive force (PMF), and accounts for the inhibitory effects of NO on the reaction kinetics. The model is able to fit well all the available kinetic data under diverse experimental conditions with a physiologically realistic unique parameter set. The model predictions show that: 1) the apparent K_m of O₂ varies considerably and increases from fully reduced to fully oxidized cytochrome c depending on pH and the energy state of mitochondria, and 2) the intermediate enzyme states depend on pH and cytochrome c redox fraction and play a central role in coupling mitochondrial respiration to PMF. The developed CcO model can easily be integrated into existing mitochondrial bioenergetics models to understand the role of the enzyme in controlling oxidative phosphorylation in normal and disease conditions.

The role of the K-channel and the active-site tyrosine in the catalytic mechanism of cytochrome c oxidase

The active site of cytochrome c oxidase (CcO) comprises an oxygen-binding heme, a nearby copper ion (CuB), and a tyrosine residue that is covalently linked to one of the histidine ligands of CuB. Two proton-conducting pathways are observed in CcO, namely the D- and the K-channels, which are used to transfer protons either to the active site of oxygen reduction (substrate protons) or for pumping. Proton transfer through the D-channel is very fast, and its role in efficient transfer of both substrate and pumped protons is well established. However, it has not been fully clear why a separate K-channel is required, apparently for the supply of substrate protons only. In this work, we have analysed the available experimental and computational data, based on which we provide new perspectives on the role of the K-channel. Our analysis suggests that proton transfer in the K-channel may be gated by the protonation state of the active-site tyrosine (Tyr244) and that the neutral radical form of this residue has a more general role in the CcO mechanism than thought previously. This article is part of a Special Issue entitled 'EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2-6, 2016', edited by Prof. Paolo Bernardi.

Protonation of the binuclear active site in cytochrome c oxidase decreases the reduction potential of CuB

One of the remaining mysteries regarding the respiratory enzyme cytochrome c oxidase is how proton pumping can occur in all reduction steps in spite of the low reduction potentials observed in equilibrium titration experiments for two of the active site cofactors, CuB(II) and Fea3(III). It has been speculated that, at least the copper cofactor can acquire two different states, one metastable activated state occurring during enzyme turnover, and one relaxed state with lower energy, reached only when the supply of electrons stops. The activated state should have a transiently increased CuB(II) reduction potential, allowing proton pumping. The relaxed state should have a lower reduction potential, as measured in the titration experiments. However, the structures of these two states are not known. Quantum mechanical calculations show that the proton coupled reduction potential for CuB is inherently high in the active site as it appears after reaction with oxygen, which explains the observed proton pumping. It is suggested here that, when the flow of electrons ceases, a relaxed resting state is formed by the uptake of one extra proton, on top of the charge compensating protons delivered in each reduction step. The extra proton in the active site decreases the proton coupled reduction potential for CuB by almost half a volt, leading to agreement with titration experiments. Furthermore, the structure for the resting state with an extra proton is found to have a hydroxo-bridge between CuB(II) and Fea3(III), yielding a magnetic coupling that can explain the experimentally observed EPR silence.

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.

Infrared and Raman spectroscopic investigation of the reaction mechanism of cytochrome c oxidase

Recent progress in studies on the proton-pumping and O₂reduction mechanisms of cytochrome c oxidase (CcO) elucidated by infrared (IR) and resonance Raman (rR) spectroscopy, is reviewed. CcO is the terminal enzyme of the respiratory chain and its O₂reduction reaction is coupled with H⁺ pumping activity across the inner mitochondrial membrane. The former is catalyzed by heme a3 and its mechanism has been determined using a rR technique, while the latter used the protein moiety and has been investigated with an IR technique. The number of H⁺ relative to e⁻ transferred in the reaction is 1:1, and their coupling is presumably performed by heme a and nearby residues. To perform this function, different parts of the protein need to cooperate with each other spontaneously and sequentially. It is the purpose of this article to describe the structural details on the coupling on the basis of the vibrational spectra of certain specified residues and chromophores involved in the reaction. Recent developments in time-resolved IR and Raman technology concomitant with protein manipulation methods have yielded profound insights into such structural changes. In particular, the new IR techniques that yielded the breakthrough are reviewed and assessed in detail. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

Computational study of the activated O(H) state in the catalytic mechanism of cytochrome c oxidase

Complex IV in the respiratory chain of mitochondria and bacteria catalyzes reduction of molecular oxygen to water, and conserves much of the liberated free energy as an electrochemical proton gradient, which is used for the synthesis of ATP. Photochemical electron injection experiments have shown that reduction of the ferric/cupric state of the enzyme's binuclear heme a3/CuB center is coupled to proton pumping across the membrane, but only if oxidation of the reduced enzyme by O2 immediately precedes electron injection. In contrast, reduction of the binuclear center in the "as-isolated" ferric/cupric enzyme is sluggish and without linkage to proton translocation. During turnover, the binuclear center apparently shuttles via a metastable but activated ferric/cupric state (O(H)), which may decay into a more stable catalytically incompetent form (O) in the absence of electron donors. The structural basis for the difference between these two states has remained elusive, and is addressed here using computational methodology. The results support the notion that CuB[II] is either three-coordinated in the O(H) state or shares an OH(-) ligand with heme a3 in a strained μ-hydroxo structure. Relaxation to state O is initiated by hydration of the binuclear site. The redox potential of CuB is expected, and found by density functional theory calculations, to be substantially higher in the O(H) state than in state O. Our calculations also suggest that the neutral radical form of the cross-linked tyrosine in the binuclear site may be more significant in the catalytic cycle than suspected so far.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Proton transfer in the quinol-dependent nitric oxide reductase from Geobacillus stearothermophilus during reduction of oxygen

Bacterial nitric oxide reductases (NOR) are integral membrane proteins that catalyse the reduction of nitric oxide to nitrous oxide, often as a step in the process of denitrification. Most functional data has been obtained with NORs that receive their electrons from a soluble cytochrome c in the periplasm and are hence termed cNOR. Very recently, the structure of a different type of NOR, the quinol-dependent (q)-NOR from the thermophilic bacterium Geobacillus stearothermophilus was solved to atomic resolution [Y. Matsumoto, T. Tosha, A.V. Pisliakov, T. Hino, H. Sugimoto, S. Nagano, Y. Sugita and Y. Shiro, Nat. Struct. Mol. Biol. 19 (2012) 238-246]. In this study, we have investigated the reaction between this qNOR and oxygen. Our results show that, like some cNORs, the G. stearothermophilus qNOR is capable of O(2) reduction with a turnover of ~3electronss(-1) at 40°C. Furthermore, using the so-called flow-flash technique, we show that the fully reduced (with three available electrons) qNOR reacts with oxygen in a reaction with a time constant of 1.8ms that oxidises the low-spin heme b. This reaction is coupled to proton uptake from solution and presumably forms a ferryl intermediate at the active site. The pH dependence of the reaction is markedly different from a corresponding reaction in cNOR from Paracoccus denitrificans, indicating that possibly the proton uptake mechanism and/or pathway differs between qNOR and cNOR. This study furthermore forms the basis for investigation of the proton transfer pathway in qNOR using both variants with putative proton transfer elements modified and measurements of the vectorial nature of the proton transfer. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).

Active site intermediates in the reduction of O(2) by cytochrome oxidase, and their derivatives

The mechanism of dioxygen activation and reduction in cell respiration, as catalysed by cytochrome c oxidase, has a long history. The work by Otto Warburg, David Keilin and Britton Chance defined the dioxygen-binding heme iron centre, viz. das Atmungsferment, or cytochrome a(3). Chance brought the field further in the mid-1970's by ingenious low-temperature studies that for the first time identified the primary enzyme-substrate (ES) Michaelis complex of cell respiration, the dioxygen adduct of heme a(3), which he termed Compound A. Further work using optical, resonance Raman, EPR, and other sophisticated spectroscopic techniques, some of which with microsecond time resolution, has brought us to the situation today, where major principles of how O(2) reduction occurs in respiration are well understood. Nonetheless, some questions have remained open, for example concerning the precise structures, catalytic roles, and spectroscopic properties of the breakdown products of Compound A that have been called P, F (for peroxy and ferryl), and O (oxidised). This nomenclature has been known to be inadequate for some time already, and an alternative will be suggested here. In addition, the multiple forms of P, F and O states have been confusing, a situation that we endeavour to help clarifying. The P and F states formed artificially by reacting cytochrome oxidase with hydrogen peroxide are especially scrutinised, and some novel interpretations will be given that may account for previously unexplained observations.

Proton-pumping mechanism of cytochrome c oxidase: a kinetic master-equation approach

Cytochrome c oxidase is an efficient energy transducer that reduces oxygen to water and converts the released chemical energy into an electrochemical membrane potential. As a true proton pump, cytochrome c oxidase translocates protons across the membrane against this potential. Based on a wealth of experiments and calculations, an increasingly detailed picture of the reaction intermediates in the redox cycle has emerged. However, the fundamental mechanism of proton pumping coupled to redox chemistry remains largely unresolved. Here we examine and extend a kinetic master-equation approach to gain insight into redox-coupled proton pumping in cytochrome c oxidase. Basic principles of the cytochrome c oxidase proton pump emerge from an analysis of the simplest kinetic models that retain essential elements of the experimentally determined structure, energetics, and kinetics, and that satisfy fundamental physical principles. The master-equation models allow us to address the question of how pumping can be achieved in a system in which all reaction steps are reversible. Whereas proton pumping does not require the direct modulation of microscopic reaction barriers, such kinetic gating greatly increases the pumping efficiency. Further efficiency gains can be achieved by partially decoupling the proton uptake pathway from the active-site region. Such a mechanism is consistent with the proposed Glu valve, in which the side chain of a key glutamic acid shuttles between the D channel and the active-site region. We also show that the models predict only small proton leaks even in the absence of turnover. The design principles identified here for cytochrome c oxidase provide a blueprint for novel biology-inspired fuel cells, and the master-equation formulation should prove useful also for other molecular machines. .

Solvent effects on the physicochemical properties of the cross-linked histidine-tyrosine ligand of cytochrome c oxidase

Density functional theory was used to explore the effects of aqueous solvation on the structure, vibrational frequencies, and the electronic absorption spectrum of 2-(4-methylimidazol-1-yl)-phenol (Me-ImPhOH), a chemical analogue of the cross-linked histidine-tyrosine Cu(B) ligand of cytochrome c oxidase. In addition, the phenolic-OH pK(a), the anodic redox potential for the biring radical/anion couple, and the phenolic-OH bond dissociation energy were calculated relative to phenol using a series of isodesmic reactions. In the gas phase, the imidazole moiety stabilizes the biring anion for all the models and greatly decreases the phenolic-OH pK(a) relative to phenol. Moreover, the conductor-like polarizable continuum model (C-PCM)-water-solvated reactions predict Delta pK(a) values that are five times smaller than the gas-phase reactions, in agreement with the proposed role of the cross-linked histidine-tyrosine as a proton donor in the enzyme. For the neutral biring radical solvation models, the imidazole moiety induces a high degree of asymmetry into the phenol ring when compared to unmodified phenoxyl radical. The biring radical pi-bonds of the imidazole ring are more localized when compared to unmodified 1-methylimidazole and Me-ImPhOH solvation models, suggesting reduced aromaticity for all biring radical solvation models. The C-PCM-water-solvated reactions predict relative biring radical reduction potentials that are an order of magnitude smaller than the gas-phase reactions. The biring O-H bond is weakened relative to phenol by less than 4 kcal/mol for all the reactions studied, suggesting that the imidazole moiety does not facilitate H-atom abstraction in the enzyme. Together, these results demonstrate the sensitive nature of the proton and electron donating ability of the histidine-tyrosine cross-linked ligand in cytochrome c oxidase and suggest that for quantitative predictions of reaction energies and thermodynamic properties, models of this ligand should take care to account for changes in environment and, more specifically, hydrogen bonding interactions.

Partial steps of charge translocation in the nonpumping N139L mutant of Rhodobacter sphaeroides cytochrome c oxidase with a blocked D-channel

The N139L substitution in the D-channel of cytochrome oxidase from Rhodobacter sphaeroides results in an approximately 15-fold decrease in the turnover number and a loss of proton pumping. Time-resolved absorption and electrometric assays of the F --> O transition in the N139L mutant oxidase result in three major findings. (1) Oxidation of the reduced enzyme by O(2) shows approximately 200-fold inhibition of the F --> O step (k approximately 2 s(-1) at pH 8) which is not compatible with enzyme turnover ( approximately 30 s(-1)). Presumably, an abnormal intermediate F(deprotonated) is formed under these conditions, one proton-deficient relative to a normal F state. In contrast, the F --> O transition in N139L oxidase induced by single-electron photoreduction of intermediate F, generated by reaction of the oxidized enzyme with H(2)O(2), decelerates to an extent compatible with enzyme turnover. (2) In the N139L mutant, the protonic phase of Deltapsi generation coupled to the flash-induced F --> O transition greatly decreases in rate and magnitude and can be assigned to the movement of a proton from E286 to the binuclear site, required for reduction of heme a(3) from the Fe(4+) horizontal lineO(2-) state to the Fe(3+)-OH(-) state. Electrogenic reprotonation of E286 from the inner aqueous phase is missing from the F --> O step in the mutant. (3) In the N139L mutant, the KCN-insensitive rapid electrogenic phase may be composed of two components with lifetimes of approximately 10 and approximately 40 mus and a magnitude ratio of approximately 3:2. The 10 mus phase matches vectorial electron transfer from Cu(A) to heme a, whereas the 40 mus component is assigned to intraprotein proton displacement across approximately 20% of the membrane dielectric thickness. This proton displacement might be triggered by rotation of the charged K362 side chain coupled to heme a reduction. The two components of the rapid electrogenic phase have been resolved subsequently with other D-channel mutants as well as with cyanide-inhibited wild-type oxidase. The finding helps to reconcile the unusually high relative contribution of the microsecond electrogenic phase in the bacterial enzyme ( approximately 30%) with the net electrogenicity of the F --> O transition coupled to transmembrane transfer of two charges per electron.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

Blocking the K-pathway still allows rapid one-electron reduction of the binuclear center during the anaerobic reduction of the aa3-type cytochrome c oxidase from Rhodobacter sphaeroides

The K-pathway is one of the two proton-input channels required for function of cytochrome c oxidase. In the Rhodobacter sphaeroides cytochrome c oxidase, the K-channel starts at Glu101 in subunit II, which is at the surface of the protein exposed to the cytoplasm, and runs to Tyr288 at the heme a3/CuB active site. Mutations of conserved, polar residues within the K-channel block or inhibit steady state oxidase activity. A large body of research has demonstrated that the K-channel is required to fully reduce the heme/Cu binuclear center, prior to the reaction with O2, presumably by providing protons to stabilize the reduced metals (ferrous heme a3 and cuprous CuB). However, there are conflicting reports which raise questions about whether blocking the K-channel blocks both electrons or only one electron from reaching the heme/Cu center. In the current work, the rate and extent of the anaerobic reduction of the heme/Cu center were monitored by optical and EPR spectroscopies, comparing the wild type and mutants that block the K-channel. The new data show that when the K-channel is blocked, one electron will still readily enter the binuclear center. The one-electron reduction of the resting oxidized ("O") heme/Cu center of the K362M mutant, results in a partially reduced binuclear center in which the electron is distributed about evenly between heme a3 and CuB in the R. sphaeroides oxidase. Complete reduction of the heme/Cu center requires the uptake of two protons which must be delivered through the K-channel.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

Electron transfer and energy transduction in the terminal part of the respiratory chain - lessons from bacterial model systems

This review focuses on the terminal part of the respiratory chain where, macroscopically speaking, electron transfer (ET) switches from the two-electron donor, ubiquinol, to the single-electron carrier, cytochrome c, to finally reduce the four-electron acceptor dioxygen. With 3-D structures of prominent representatives of such multi-subunit membrane complexes known for some time, this section of the ET chain still leaves a number of key questions unanswered. The two relevant enzymes, ubiquinol:cytochrome c oxidoreductase and cytochrome c oxidase, appear as rather diverse modules, differing largely in their design for substrate interaction, internal ET, and moreover, in their mechanisms of energy transduction. While the canonical mitochondrial complexes have been investigated for almost five decades, the corresponding bacterial enzymes have been established only recently as attractive model systems to address basic reactions in ET and energy transduction. Lacking the intricate coding background and mitochondrial assembly pathways, bacterial respiratory enzymes typically offer a much simpler subunit composition, while maintaining all fundamental functions established for their complex "relatives". Moreover, related issues ranging from primary steps in cofactor insertion to supramolecular architecture of ET complexes, can also be favourably addressed in prokaryotic systems to hone our views on prototypic structures and mechanisms common to all family members.