Theoretical and computational analysis of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme

We have developed theory and the computational scheme for the analysis of the kinetics of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme. The theory allows one to connect the charge motions inside the enzyme to the membrane potential observed in the experiments by using data from the "dielectric topography" map of the enzyme that we have created. The developed theory is applied for the analysis of the potentiometric data recently reported by the Wikstr?m group [I. Belevich, D.A. Bloch, N. Belevich, M. Wikstr?m and M.I. Verkhovsky, Exploring the proton pump mechanism of cytochrome c oxidase in real time, Proc. Natl. Acad. Sci. U. S. A. 104 (2007) 2685-2690] on the O to E transition in Paracoccus denitrificans oxidase. Our analysis suggests, that the electron transfer to the binuclear center is coupled to a proton transfer (proton loading) to a group just "above" the binuclear center of the enzyme, from which the pumped proton is subsequently expelled by the chemical proton arriving to the binuclear center. The identity of the pump site could not be determined with certainty, but could be localized to the group of residues His326 (His291 in bovine), propionates of heme a(3), Arg 473/474, and Trp164. The analysis also suggests that the dielectric distance from the P-side to Fe a is 0.4 or larger. The difficulties and pitfalls of quantitative interpretation of potentiometric data are discussed.     

Towards the mechanism of proton pumping by the haem-copper oxidases

The haem-copper oxidases comprise a large family of enzymes that is widespread among aerobic organisms. These remarkable membrane-bound proteins catalyse the respiratory reduction of dioxygen to water, and conserve free energy from this reaction by operating as proton pumps. The mechanism of redox-dependent proton translocation has been elusive despite the availability of high resolution crystal structures from several oxidases. Here, we discuss some recent as well as some older results that may shed light on this mechanism. We conclude that proton-pumping is initiated by vectorial proton transfer from a conserved glutamic acid (Glu242 in the bovine enzyme) to a proton acceptor above the haem groups, and that this primary event is mechanistically coupled to electron transfer from haem a to the binuclear haem a₃/Cu_B centre. Subsequently, Glu242 is reprotonated from the negatively charged side of the membrane. Next this proton is transferred to the binuclear site to complete the chemistry, Glu242 is reprotonated once more, and the “prepumped” proton is ejected on the opposite side of the membrane. The different kinetics of electron-coupled proton transfer in different steps of the catalytic cycle may be related to differences in the driving force due to different E_m values of the electron acceptor in the binuclear site.     

Mechanism and energetics of proton translocation by the respiratory heme-copper oxidases

Recent time-resolved optical and electrometric experiments have provided a sequence of events for the proton-translocating mechanism of cytochrome c oxidase. These data also set limits for the mechanistic, kinetic, and thermodynamic parameters of the proton pump, which are analysed here in some detail. The analysis yields limit values for the pK of the “pump site”, its modulation during the proton-pumping process, and suggests its identity in the structure. Special emphasis is made on side-reactions that may short-circuit the pump, and the means by which these may be avoided. We will also discuss the most prominent proton pumping mechanisms proposed to date in relation to these data.     

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's Cu_B center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model—the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from δ-nitrogen of His291 ligand of Cu_B center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pK_a of His291 for an extended model in which the whole Fe_a3-Cu_B center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pK_a units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.     

Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a₃/Cu_B site, and for proton pumping. When present, the side chain of Glu-242 is orientated “down” towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity “above” Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an “up” position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the “up” position is preferred energetically when the cavity contains four water molecules, but the “down” position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the “resting” state of the enzyme, and that the “up/down” equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O₂ reduction.     

Spectroscopic study on the communication between a heme a3 propionate, Asp399 and the binuclear center of cytochrome c oxidase from Paracoccus denitrificans

The proton pumping mechanism of cytochrome c oxidase on a molecular level is highly disputed. Recently theoretical calculations and real time electron transfer measurements indicated the involvement of residues in the vicinity of the ring A propionate of heme a₃, including Asp399 and the Cu_B ligands His 325, 326. In this study we probed the interaction of Asp399 with the binuclear center and characterize the protonation state of its side chain. Redox induced FTIR difference spectra of mutations at the site in direct comparison to wild type, indicate that below pH 5 Asp 399 displays signals typical for the deprotonation of the acidic residue with reduction of the enzyme. Interestingly at a pH higher than 5, no contributions from Asp 399 are evident. In order to probe the interaction of the site with the binuclear center we followed the rebinding of CO by infrared spectroscopy for mutations on residue Asp399 to Glu, Asn and Leu. Previously different CO conformers have been identified for bacterial cytochrome c oxidases, and its pH dependent behaviour discussed to be relevant for catalysis. Interestingly we observe the lack of this pH dependency and a strong influence on the observable conformers for all mutants studied here, clearly suggesting a communication of the site with the heme-copper center and the nearby histidine residues.     

Cytochrome oxidase as a proton pump

The general structure of cytochrome oxidase is reviewed and evidence that the enzyme acts as a redox-linked proton pump outlined. The overall H⁺/e ^– stoichiometry of the pump is discussed and results [Wikström (1989),Nature 338, 293] which suggest that only the final two electrons which reduce the peroxide adduct to water are coupled to protein translocated are considered in terms of the restrictions they place on pump mechanisms. Direct and indirect mechanisms for proton translocation are discussed in the context of evidence for redox-linked conformational changes in the enzyme, the role of subunit III, and the nature of the Cu_A site.     

Internal charge transfer in cytochrome c oxidase at a limited proton supply: Proton pumping ceases at high pH

Background

In the membrane-bound enzyme cytochrome c oxidase, electron transfer from cytochrome c to O₂ is linked to proton uptake from solution to form H₂O, resulting in a charge separation across the membrane. In addition, the reaction drives pumping of protons across the membrane.

Methods

In this study we have measured voltage changes as a function of pH during reaction of the four-electron reduced cytochrome c oxidase from Rhodobacter sphaeroides with O₂. These electrogenic events were measured across membranes containing purified enzyme reconstituted into lipid vesicles.

Results

The results show that the pH dependence of voltage changes (primarily associated with proton transfer) during O₂ reduction does not match that of the previously studied absorbance changes (primarily associated with electron transfer). Furthermore, the voltage changes decrease with increasing pH.

Conclusions

The data indicate that cytochrome c oxidase does not pump protons at high pH (10.5) (or protons are taken from the “wrong” side of the membrane) and that at this pH the net proton-uptake stoichiometry is ∼ 1/2 of that at pH 8. Furthermore, the results provide a basis for interpretation of results from studies of mutant forms of the enzyme.

General significance

These results provide new insights into the function of cytochrome c oxidase.     

The inside pH determines rates of electron and proton transfer in vesicle-reconstituted cytochrome c oxidase

Cytochrome c oxidase is the terminal enzyme in the respiratory chains of mitochondria and many bacteria where it translocates protons across a membrane thereby maintaining an electrochemical proton gradient. Results from earlier studies on detergent-solubilized cytochrome c oxidase have shown that individual reaction steps associated with proton pumping display pH-dependent kinetics. Here, we investigated the effect of pH on the kinetics of these reaction steps with membrane-reconstituted cytochrome c oxidase such that the pH was adjusted to different values on the inside and outside of the membrane. The results show that the pH on the inside of the membrane fully determines the kinetics of internal electron transfers that are linked to proton pumping. Thus, even though proton release is rate limiting for these reaction steps (Salomonsson et al., Proc. Natl. Acad. Sci. USA, 2005, 102, 17624), the transition kinetics is insensitive to the outside pH (in the range 6-9.5).     

A pathogenic mutation in cytochrome c oxidase results in impaired proton pumping while retaining O2-reduction activity

In this work we have investigated the effect of a pathogenic mitochondrial DNA mutation found in human colon cells, at a functional-molecular level. The mutation results in the amino-acid substitution Tyr19His in subunit I of the human CytcO and it is associated with respiratory deficiency. It was introduced into Rhodobacter sphaeroides, which carries a cytochrome c oxidase (cytochrome aa₃) that serves as a model of the mitochondrial counterpart. The residue is situated in the middle of a pathway that is used to transfer substrate protons as well as protons that are pumped across the membrane. The Tyr33His (equivalent residue in the bacterial CytcO) structural variant of the enzyme was purified and its function was investigated. The results show that in the structurally altered CytcO the activity decreased due to slowed proton transfer; proton transfer from an internal proton donor, the highly-conserved Glu286, to the catalytic site was slowed by a factor of ∼ 5, while reprotonation of the Glu from solution was slowed by a factor of ∼ 40. In addition, in the structural variant proton pumping was completely impaired. These results are explained in terms of introduction of a barrier for proton transfer through the D pathway and changes in the coordination of water molecules surrounding the Glu286 residue. The study offers an explanation, at the molecular level, to the link between a specific amino-acid substitution and a pathogenic phenotype identified in human colon cells.     

Molecular dynamics simulation of water in cytochrome c oxidase reveals two water exit pathways and the mechanism of transport

We have examined the network of connected internal cavities in cytochrome c oxidase along which water produced at the catalytic center is removed from the enzyme. Using combination of structural analysis, molecular dynamics simulations, and free energy calculations we have identified two exit pathways that connect the Mg²⁺ ion cavity to the outside of the enzyme. Each pathway has a well-defined bottleneck, which determines the overall rate of water traffic along the exit pathway, and a specific cooperative mechanism of passing it. One of the pathways is going via Arg438/439 (in bovine numbering) toward the Cu_A center, approaching closely its His204_B ligand and Lys171_B residue; and the other is going toward Asp364 and Thr294. Comparison of the pathways among different aa₃-type enzymes shows that they are well conserved. Possible connections of the finding to redox-coupled proton pumping mechanism are discussed. We propose specific mutations near the bottlenecks of the exit pathways that can test some of our hypotheses.     

Electrostatic control of proton pumping in cytochrome c oxidase

As part of the mitochondrial respiratory chain, cytochrome c oxidase utilizes the energy produced by the reduction of O₂ to water to fuel vectorial proton transport. The mechanism coupling proton pumping to redox chemistry is unknown. Recent advances have provided evidence that each of the four observable transitions in the complex catalytic cycle consists of a similar sequence of events. However, the physico-chemical basis underlying this recurring sequence has not been identified. We identify this recurring pattern based on a comprehensive model of the catalytic cycle derived from the analysis of oxygen chemistry and available experimental evidence. The catalytic cycle involves the periodic repetition of a sequence of three states differing in the spatial distribution of charge in the active site: [0|1], [1|0], and [1|1], where the total charge of heme a and the binuclear center appears on the left and on the right, respectively. This sequence recurs four times per turnover despite differences in the redox chemistry. This model leads to a simple, robust, and reproducible sequence of electron and proton transfer steps and rationalizes the pumping mechanism in terms of electrostatic coupling of proton translocation to redox chemistry. Continuum electrostatic calculations support the proposed mechanism and suggest an electrostatic origin for the decoupled and inactive phenotypes of ionic mutants in the principal proton-uptake pathway.     

Proton pumping mechanism and catalytic cycle of cytochrome c oxidase: Coulomb pump model with kinetic gating

Using electrostatic calculations, we have examined the dependence of the protonation state of cytochrome c oxidase from bovine heart on its redox state. Based on these calculations, we propose a possible scheme of redox-linked proton pumping. The scheme involves His291 - one of the ligands of the Cu(B) redox center - which plays the role of the proton loading site (PLS) of the pump. The mechanism of pumping is based on ET reaction between two hemes of the enzyme, which is coupled to a transfer of two protons. Upon ET, the first proton (fast reaction) is transferred to the PLS (His291), while subsequent transfer of the second "chemical" proton to the binuclear center (slow reaction) is accompanied by the ejection of the first (pumped) proton. Within the proposed model, we discuss the catalytic cycle of the enzyme.     

Molecular dynamics in cytochrome c oxidase Mössbauer spectra deconvolution

In this work low temperature molecular dynamics simulations of cytochrome c oxidase are used to predict an experimentally observable, namely Mössbauer spectra width. Predicted lineshapes are used to model Lorentzian doublets, with which published cytochrome c oxidase Mössbauer spectra were simulated. Molecular dynamics imposed constraints to spectral lineshapes permit to obtain useful information, like the presence of multiple chemical species in the binuclear center of cytochrome c oxidase. Moreover, a benchmark of quality for molecular dynamic simulations can be obtained. Despite the overwhelming importance of dynamics in electron–proton transfer systems, limited work has been devoted to unravel how much realistic are molecular dynamics simulations results. In this work, molecular dynamics based predictions are found to be in good agreement with published experimental spectra, showing that we can confidently rely on actual simulations. Molecular dynamics based deconvolution of Mössbauer spectra will lead to a renewed interest for application of this approach in bioenergetics.     

Inhibition of membrane-bound cytochrome c oxidase by zinc ions: High-affinity Zn-binding site at the P-side of the membrane

In the presence of the uncoupler, external zinc ions inhibit rapidly turnover of cytochrome c oxidase reconstituted in phospholipid vesicles or bound to the membrane of intact mitochondria. The effect is promoted by electron leaks into the oxidase during preincubation with Zn²⁺. Inhibition of liposome-bound bovine cytochrome oxidase by external Zn²⁺ titrates with a K_i of 1 ± 0.3 μM. Presumably, the Zn²⁺-binding group at the positively charged side is not reactive in the oxidized enzyme, but becomes accessible to the cation in some partially reduced state(s) of the oxidase; reduction of Cu_B is tentatively proposed to be responsible for the effect.     

Mechanism and energetics by which glutamic acid 242 prevents leaks in cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme of aerobic respiration. The energy released from the reduction of molecular oxygen to water is used to pump protons across the mitochondrial or bacterial membrane. The pump function introduces a mechanistic requirement of a valve that prevents protons from flowing backwards during the process. It was recently found that Glu-242, a key amino acid in transferring protons to be pumped across the membrane and to the site of oxygen reduction, fulfils the function of such a valve by preventing simultaneous contact to the pump site and to the proton-conducting D-channel (Kaila V.R.I. et al. Proc. Natl. Acad. Sci. USA 105, 2008). Here we have incorporated the valve model into the framework of the reaction mechanism. The function of the Glu valve is studied by exploring how the redox state of the surrounding metal centers, dielectric effects, and membrane potential, affects the energetics and leaks of this valve. Parallels are drawn between the dynamics of Glu-242 and the long-standing obscure difference between the metastable O_H and stable O states of the binuclear center. Our model provides a suggestion for why reduction of the former state is coupled to proton translocation while reduction of the latter is not.     

The proton pump of heme-copper oxidases.

Proton pumping heme-copper oxidases represent the terminal, energy-transfer enzymes of respiratory chains in prokaryotes and eukaryotes. The Cu_B-heme a₃ (or heme o) binuclear center, associated with the largest subunit I of cytochrome c and quinol oxidases, is directly involved in the coupling between dioxygen reduction and proton pumping. The role of the other subunits is less clear. The following aspects will be covered in this paper:i) the efficiency of coupling in the mitochondrial aa₃ cytochrome c oxidase. In particular, the effect of respiratory rate and protonmotive force on the H⁺/e^? stoichiometry and the role of subunit IV; ii) mutational analysis of the aa₃ quinol oxidase of Bacillus subtilis addressed to the role of subunit III, subunit IV and specific residues in subunit I; iii) possible models of the protonmotive catalytic cycle at the binuclear center. The observations available suggest that H⁺/e^?coupling is based on the combination of protonmotive redox catalysis at the binuclear center and co-operative proton transfer in the protein.     

Time-resolved OH → EH transition of the aberrant ba3 oxidase from Thermus thermophilus

The kinetics of single-electron injection into the oxidized nonrelaxed state (O_H → E_H transition) of the aberrant ba₃ cytochrome oxidase from Thermus thermophilus, noted for its lowered efficiency of proton pumping, was investigated by time-resolved optical spectroscopy. Two main phases of intraprotein electron transfer were resolved. The first component (τ ∼ 17 μs) reflects oxidation of Cu_A and reduction of the heme groups (low-spin heme b and high-spin heme a₃ in a ratio close to 50:50). The subsequent component (τ ∼ 420 μs) includes reoxidation of both hemes by Cu_B. This is in significant contrast to the O_H → E_H transition of the aa₃-type cytochrome oxidase from Paracoccus denitrificans, where the fastest phase is exclusively due to transient reduction of the low-spin heme a, without electron equilibration with the binuclear center. On the other hand, the one-electron reduction of the relaxed O state in ba₃ oxidase was similar to that in aa₃ oxidase and only included rapid electron transfer from Cu_A to the low-spin heme b. This indicates a functional difference between the relaxed O and the pulsed O_H forms also in the ba₃ oxidase from T. thermophilus.     

Knockdown of Human COX17 Affects Assembly and Supramolecular Organization of Cytochrome c Oxidase

Assembly of cytochrome c oxidase, the terminal enzyme of the mitochondrial respiratory chain, requires a concerted activity of a number of chaperones and factors for the insertion of subunits, accessory proteins, cofactors and prosthetic groups. It is now well accepted that the multienzyme complexes of the respiratory chain are organized in vivo as supramolecular functional structures, so-called supercomplexes. Here, we investigate the role of COX17 in the biogenesis of the respiratory chain in HeLa cells. In accordance with its predicted function as a copper chaperone and its role in formation of the binuclear copper centre of cytochrome c oxidase, COX17 siRNA knockdown affects activity and assembly of cytochrome c oxidase. While the abundance of cytochrome c oxidase dimers seems to be unaffected, blue native gel electrophoresis reveals the disappearance of COX-containing supercomplexes as an early response. We observe the accumulation of a novel ∼ 150 kDa complex that contains Cox1, but not Cox2. This observation may indicate that the absence of Cox17 interferes with copper delivery to Cox2, but not to Cox1. We suggest that supercomplex formation is not simply due to assembly of completely assembled complexes. An interdependent assembly scenario for the formation of supercomplexes that rather requires the coordinated synthesis and association of individual complexes, is proposed.