A chemically explicit model for the mechanism of proton pumping in heme-copper oxidases

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

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.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

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.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

Modulation of the electron-proton coupling at cytochrome a by the ligation of the oxidized catalytic center in bovine cytochrome c oxidase

Cytochrome a was suggested as the key redox center in the proton pumping process of bovine cytochrome c oxidase (CcO). Recent studies showed that both the structure of heme a and its immediate vicinity are sensitive to the ligation and the redox state of the distant catalytic center composed of iron of cytochrome a₃ (Fe_a3) and copper (Cu_B). Here, the influence of the ligation at the oxidized Fe_a3³⁺-Cu_B²⁺ center on the electron-proton coupling at heme a was examined in the wide pH range (6.5-11). The strength of the coupling was evaluated by the determination of pH dependence of the midpoint potential of heme a (E_m(a)) for the cyanide (the low-spin Fe_a3³⁺) and the formate-ligated CcO (the high-spin Fe_a3³⁺). The measurements were performed under experimental conditions when other three redox centers of CcO are oxidized. Two slightly differing linear pH dependencies of E_m(a) were found for the CN- and the formate-ligated CcO with slopes of -13 mV/pH unit and -23 mV/pH unit, respectively. These linear dependencies indicate only a weak and unspecific electron-proton coupling at cytochrome a in both forms of CcO. The lack of the strong electron-proton coupling at the physiological pH values is also substantiated by the UV-Vis absorption and electron-paramagnetic resonance spectroscopy investigations of the cyanide-ligated oxidized CcO. It is shown that the ligand exchange at Fe_a³⁺ between His-Fe_a³⁺-His and His-Fe_a³⁺-OH^- occurs only at pH above 9.5 with the estimated pK >11.0.

Explaining leak states in the proton pump of heme-copper oxidases observed in single-molecule experiments

Heme-copper oxidases couple the exergonic oxygen reduction with the endergonic proton translocation. Redox-linked structural changes have been localized in deeply buried regions of the protein, near the low-potential heme. How these movements can modulate distant gating events along the intramolecular proton path, where the entry (exit) of pumped proton occurs, is a major concern for the proton pump models. Generally, these models associate, more or less directly, all translocation events with redox transitions. Although they can account for many phenomenological aspects of the pump, evidences from single-molecules experiments about leak states of the pump represent a formidable challenge. Disconnecting the redox-linked pK_a shifts of the proton loading site from the external barriers, we obtain a simple stochastic mechanism which behaves similarly to the real enzyme, able to reverse the flow of the proton transfer.

Kinetics and Efficiency of Energy-Transducing Enzymes

Complexes I to IV, with the exception of Complex II, are redox-driven proton pumps that convert redox energy of oxygen reduction to proton gradient across the mitochondrial or bacterial membrane; in turn, the created electrochemical gradient drives the adenosine triphosphate synthesis in the cells by utilizing complex V of the chain. Here we address a general question of the efficiency of such enzymes, considering them as molecular machines that couple endergonic and exergonic reactions and converting one form of free energy into another. One well-known example of the efficiency is given by Carnot's theorem for heat engines. Here we extend the concept to respiratory enzymes and specifically focus on the proton pumping by Complex I of the respiratory chain, nicotinamide adenine dinucleotide dehydrogenase. To discuss the efficiency issues, we develop a model of enzyme kinetics, which generalizes the Michaelis-Menten model. Our model includes several substrates and products and, in general, can be considered as Generalized Michaelis-Menten Kinetic model. The model might be useful for describing complex enzyme kinetics, regardless of the efficiency issues that are addressed in this paper.

Redox-Driven Proton Pumps of the Respiratory Chain

In aerobic cells, the proton gradient that drives ATP synthesis is created by three different proton pumps-membrane enzymes of the respiratory electron transport chain known as complex I, III, and IV. Despite the striking dissimilarity of structures and apparent differences in molecular mechanisms of proton pumping, all three enzymes have much in common and employ the same universal physical principles of converting redox energy to proton pumping. In this study, we describe a simple mathematical model that illustrates the general principles of redox-driven proton pumps and discuss their implementation in complex I, III, and IV of the respiratory chain.

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.

Response of Heme Symmetry to the Redox State of Bovine Cytochrome c Oxidase

Second-derivative absorption spectroscopy was employed to monitor the response of effective symmetry of cytochromes a and a₃ to the redox and ligation states of bovine cytochrome c oxidase (CcO). The Soret band π → π* electronic transitions were used to display the changes in symmetry of these chromophores induced by the reduction of CcO inhibited by the exogenous ligands and during catalytic turnover. The second derivative of the difference absorption spectra revealed only a single Soret band for the oxidized cytochromes a and a₃ and cyanide-ligated oxidized cytochrome a₃. In contrast, two absorption bands were resolved in ferrous cytochrome a and ferrous cytochrome a₃ ligated with cyanide. A transition from one-band spectrum to two-band spectrum indicates the lowering of symmetry of these hemes due to the alteration of their immediate surroundings. It is suggested that the changes in polarity occurring in the vicinity of these cofactors are main reason for the split of the Soret band of both ferrous cytochrome a and cyanide-bound ferrous cytochrome a₃.

Protonation-State-Dependent Communication in Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase from the cellular inside to the binuclear redox center (BNC) can occur through two distinct pathways, the D- and K-channels. For the protein to function as both a redox enzyme and a proton pump, proton transfer into the protein toward the BNC or toward a proton loading site (and ultimately through the membrane) must be highly regulated. The P_R → F transition is the first step in a catalytic cycle that requires proton transfer from the bulk at the N-side to the BNC. Molecular dynamics simulations of the P_R → F intermediate of this transition, with 16 different combinations of protonation states of key residues in the D- and K-channel, show the impact of the K-channel on the D-channel to be protonation-state dependent. Strength as well as means of communication, correlations in positions, or communication along the hydrogen-bonded network depends on the protonation state of the K-channel residue K362. The conformational and hydrogen-bond dynamics of the D-channel residue N139 is regulated by an interplay of protonation in the D-channel and K362. N139 thus assumes a gating function by which proton passage through the D-channel toward E286 is likely facilitated for states with protonated K362 and unprotonated E286. In contrast, proton passage through the D-channel is hindered by N139's preference for a closed conformation in situations with protonated E286.

Role of conformational change and K-path ligands in controlling cytochrome c oxidase activity

Given the central role of cytochrome c oxidase (CcO) in health and disease, it is an increasingly important question as to how the activity and efficiency of this key enzyme are regulated to respond to a variety of metabolic states. The present paper summarizes evidence for two modes of regulation of activity: first, by redox-induced conformational changes involving the K-proton uptake path; and secondly, by ligand binding to a conserved site immediately adjacent to the entrance of the K-path that leads to the active site. Both these phenomena highlight the importance of the K-path in control of CcO. The redox-induced structural changes are seen in both the two-subunit and a new four-subunit crystal structure of bacterial CcO and suggest a gating mechanism to control access of protons to the active site. A conserved ligand-binding site, first discovered as a bile salt/steroid site in bacterial and mammalian oxidases, is observed to bind an array of ligands, including nucleotides, detergents, and other amphipathic molecules. Highly variable effects on activity, seen for these ligands and mutations at the K-path entrance, can be explained by differing abilities to inhibit or stimulate K-path proton uptake by preventing or allowing water organization. A new mutant form in which the K-path is blocked by substituting the conserved carboxyl with a tryptophan clarifies the singularity of the K-path entrance site. Further study in eukaryotic systems will determine the physiological significance and pharmacological potential of ligand binding and conformational change in CcO.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Calcium ions inhibit reduction of heme a in bovine cytochrome c oxidase

The effect of Ca(2+) on the rate of heme a reduction by dithionite and hexaammineruthenium (RuAm) was studied in the cyanide-complexed bovine cytochrome oxidase (CcO). The rate of heme a reduction is proportional to RuAm concentration below 300 μM with kv of 0.53×10(6) M(-1) s(-1). Ca(2+) inhibits the rate of heme a reduction by dithionite by ∼25%. As the reaction speeds up with increased concentrations of RuAm, the inhibition by Ca(2+) disappears. The inhibition of heme a reduction may contribute to recently described partial inhibition of CcO by Ca(2+) in the enzymatic assays. The inhibitory effect of Ca(2+) on heme a reduction indicates that ET through heme a may be coupled to proton movement in the exit part of the proton channel H.

Lysine 362 in cytochrome c oxidase regulates opening of the K-channel via changes in pKA and conformation

The metabolism of aerobic life uses the conversion of molecular oxygen to water as an energy source. This reaction is catalyzed by cytochrome e oxidase (CeO) consuming four electrons and four protons, which move along specific routes. While all four electrons are transferred via the same cofactors to the binuclear reaction center (BNC), the protons take two different routes in the A-type CeO, i.e., two of the four chemical protons consumed in the reaction arrive via the D-channel in the oxidative first half starting after oxygen binding. The other two chemical protons enter via the K-channel in the reductive second half of the reaction cycle. To date, the mechanism behind these separate proton transport pathways has not been understood. In this study, we propose a model that can explain the reaction-step specific opening and closing of the K-channel by conformational and pKA changes of its central lysine 362. Molecular dynamics simulations reveal an upward movement of Lys362 towards the BNC, which had already been supposed by several experimental studies. Redox state-dependent pKA calculations provide evidence that Lys362 may protonate transiently, thereby opening the K-channel only in the reductive second half of the reaction cycle. From our results, we develop a model that assigns a key role to Lys362 in the proton gating between the two proton input channels of the A-type CeO.

Mutation of a single residue in the ba3 oxidase specifically impairs protonation of the pump site

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound protein complex that couples electron transfer to O2 to proton translocation across the membrane. To elucidate the mechanism of the redox-driven proton pumping, we investigated the kinetics of electron and proton transfer in a structural variant of the ba3 oxidase where a putative "pump site" was modified by replacement of Asp372 by Ile. In this structural variant, proton pumping was uncoupled from internal electron transfer and O2 reduction. The results from our studies show that proton uptake to the pump site (time constant ∼65 μs in the wild-type cytochrome c oxidase) was impaired in the Asp372Ile variant. Furthermore, a reaction step that in the wild-type cytochrome c oxidase is linked to simultaneous proton uptake and release with a time constant of ∼1.2 ms was slowed to ∼8.4 ms, and in Asp372Ile was only associated with proton uptake to the catalytic site. These data identify reaction steps that are associated with protonation and deprotonation of the pump site, and point to the area around Asp372 as the location of this site in the ba3 cytochrome c oxidase.

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.

Characterizing the proton loading site in cytochrome c oxidase

Cytochrome c oxidase (CcO) uses the energy released by reduction of O2 to H2O to drive eight charges from the high pH to low pH side of the membrane, increasing the electrochemical gradient. Four electrons and protons are used for chemistry, while four more protons are pumped. Proton pumping requires that residues on a pathway change proton affinity through the reaction cycle to load and then release protons. The protonation states of all residues in CcO are determined in MultiConformational Continuum Electrostatics simulations with the protonation and redox states of heme a, a3, Cu(B), Y288, and E286 used to define the catalytic cycle. One proton is found to be loaded and released from residues identified as the proton loading site (PLS) on the P-side of the protein in each of the four CcO redox states. Thus, the same proton pumping mechanism can be used each time CcO is reduced. Calculations with structures of Rhodobacter sphaeroides, Paracoccus denitrificans, and bovine CcO derived by crystallography and molecular dynamics show the PLS functions similarly in different CcO species. The PLS is a cluster rather than a single residue, as different structures show 1-4 residues load and release protons. However, the proton affinity of the heme a3 propionic acids primarily determines the number of protons loaded into the PLS; if their proton affinity is too low, less than one proton is loaded.

Unusual post-translational protein modifications: the benefits of sophistication

This review summarizes the “seemingly bizarre”, yet naturally occurring, covalent non-disulphide cross-links in enzymatic and scaffolding proteins and their functions.

Computational prediction and in vitro analysis of potential physiological ligands of the bile acid binding site in cytochrome c oxidase

A conserved bile acid site has been crystallographically defined in the membrane domain of mammalian and Rhodobacter sphaeroides cytochrome c oxidase (RsCcO). Diverse amphipathic ligands were shown previously to bind to this site and affect the electron transfer equilibrium between heme a and a3 cofactors by blocking the K proton uptake path. Current studies identify physiologically relevant ligands for the bile acid site using a novel three-pronged computational approach: ROCS comparison of ligand shape and electrostatics, SimSite3D comparison of ligand binding site features, and SLIDE screening of potential ligands by docking. Identified candidate ligands include steroids, nicotinamides, flavins, nucleotides, retinoic acid, and thyroid hormones, which are predicted to make key protein contacts with the residues involved in bile acid binding. In vitro oxygen consumption and ligand competition assays on RsCcO wildtype and its Glu101Ala mutant support regulatory activity and specificity of some of these ligands. An ATP analog and GDP inhibit RsCcO under low substrate conditions, while fusidic acid, cholesteryl hemisuccinate, retinoic acid, and T3 thyroid hormone are more potent inhibitors under both high and low substrate conditions. The sigmoidal kinetics of RsCcO inhibition in the presence of certain nucleotides is reminiscent of previously reported ATP inhibition of mammalian CcO, suggesting regulation involving the conserved core subunits of both mammalian and bacterial oxidases. Ligand binding to the bile acid site is noncompetitive with respect to cytochrome c and appears to arrest CcO in a semioxidized state with some resemblance to the "resting" state of the enzyme.

Functions of the hydrophilic channels in protonmotive cytochrome c oxidase

The structures and functions of hydrophilic channels in electron-transferring membrane proteins are discussed. A distinction is made between proton channels that can conduct protons and dielectric channels that are non-conducting but can dielectrically polarize in response to the introduction of charge changes in buried functional centres. Functions of the K, D and H channels found in A1-type cytochrome c oxidases are reviewed in relation to these ideas. Possible control of function by dielectric channels and their evolutionary relation to proton channels is explored.

Current advances in research of cytochrome c oxidase

The function of cytochrome c oxidase as a biomolecular nanomachine that transforms energy of redox reaction into protonmotive force across a biological membrane has been subject of intense research, debate, and controversy. The structure of the enzyme has been solved for several organisms; however details of its molecular mechanism of proton pumping still remain elusive. Particularly, the identity of the proton pumping site, the key element of the mechanism, is still open to dispute. The pumping mechanism has been for a long time one of the key unsolved issues of bioenergetics and biochemistry, but with the accelerating progress in this field many important details and principles have emerged. Current advances in cytochrome oxidase research are reviewed here, along with a brief discussion of the most complete proton pumping mechanism proposed to date, and a molecular basis for control of its efficiency.

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.

Redox-controlled proton gating in bovine cytochrome c oxidase

Cytochrome c oxidase is the terminal enzyme in the electron transfer chain of essentially all organisms that utilize oxygen to generate energy. It reduces oxygen to water and harnesses the energy to pump protons across the mitochondrial membrane in eukaryotes and the plasma membrane in prokaryotes. The mechanism by which proton pumping is coupled to the oxygen reduction reaction remains unresolved, owing to the difficulty of visualizing proton movement within the massive membrane-associated protein matrix. Here, with a novel hydrogen/deuterium exchange resonance Raman spectroscopy method, we have identified two critical elements of the proton pump: a proton loading site near the propionate groups of heme a, which is capable of transiently storing protons uploaded from the negative-side of the membrane prior to their release into the positive side of the membrane and a conformational gate that controls proton translocation in response to the change in the redox state of heme a. These findings form the basis for a postulated molecular model describing a detailed mechanism by which unidirectional proton translocation is coupled to electron transfer from heme a to heme a 3, associated with the oxygen chemistry occurring in the heme a 3 site, during enzymatic turnover.

A conserved amphipathic ligand binding region influences k-path-dependent activity of cytochrome C oxidase

A conserved, crystallographically defined bile acid binding site was originally identified in the membrane domain of mammalian and bacterial cytochrome c oxidase (CcO). Current studies show other amphipathic molecules including detergents, fatty acids, steroids, and porphyrins bind to this site and affect the already 50% inhibited activity of the E101A mutant of Rhodobacter sphaeroides CcO as well as altering the activity of wild-type and bovine enzymes. Dodecyl maltoside, Triton X100, C12E8, lysophophatidylcholine, and CHOBIMALT detergents further inhibit RsCcO E101A, with lesser inhibition observed in wild-type. The detergent inhibition is overcome in the presence of micromolar concentrations of steroids and porphyrin analogues including deoxycholate, cholesteryl hemisuccinate, bilirubin, and protoporphyrin IX. In addition to alleviating detergent inhibition, amphipathic carboxylates including arachidonic, docosahexanoic, and phytanic acids stimulate the activity of E101A to wild-type levels by providing the missing carboxyl group. Computational modeling of dodecyl maltoside, bilirubin, and protoporphyrin IX into the conserved steroid site shows energetically favorable binding modes for these ligands and suggests that a groove at the interface of subunit I and II, including the entrance to the K-path and helix VIII of subunit I, mediates the observed competitive ligand interactions involving two overlapping sites. Spectral analysis indicates that ligand binding to this region affects CcO activity by altering the K-path-dependent electron transfer equilibrium between heme a and heme a(3). The high affinity and specificity of a number of compounds for this region, and its conservation and impact on CcO activity, support its physiological significance.

From static structure to living protein: computational analysis of cytochrome c oxidase main-chain flexibility

Crystallographic structure and deuterium accessibility comparisons of CcO in different redox states have suggested conformational changes of mechanistic significance. To predict the intrinsic flexibility and low energy motions in CcO, this work has analyzed available high-resolution crystallographic structures with ProFlex and elNémo computational methods. The results identify flexible regions and potential conformational changes in CcO that correlate well with published structural and biochemical data and provide mechanistic insights. CcO is predicted to undergo rotational motions on the interior and exterior of the membrane, driven by transmembrane helical tilting and bending, coupled with rocking of the β-sheet domain. Consequently, the proton K-pathway becomes sufficiently flexible for internal water molecules to alternately occupy upper and lower parts of the pathway, associated with conserved Thr-359 and Lys-362 residues. The D-pathway helices are found to be relatively rigid, with a highly flexible entrance region involving the subunit I C-terminus, potentially regulating the uptake of protons. Constriction and dilation of hydrophobic channels in RsCcO suggest regulation of the oxygen supply to the binuclear center. This analysis points to coupled conformational changes in CcO and their potential to influence both proton and oxygen access.

Product-controlled steady-state kinetics between cytochrome aa(3) from Rhodobacter sphaeroides and equine ferrocytochrome c analyzed by a novel spectrophotometric approach

Cytochrome c oxidase (CcO) catalyzes the reduction of molecular oxygen to water using ferrocytochrome c (cyt c(2+)) as the electron donor. In this study, the oxidation of horse cyt c(2+) by CcO from Rhodobacter sphaeroides, was monitored using stopped-flow spectrophotometry. A novel analytic procedure was applied in which the spectra were deconvoluted into the reduced and oxidized forms of cyt c by a least-squares fitting method, yielding the reaction rates at various concentrations of cyt c(2+) and cyt c(3+). This allowed an analysis of the effects of cyt c(3+) on the steady-state kinetics between CcO and cyt c(2+). The results show that cyt c(3+) exhibits product inhibition by two mechanisms: competition with cyt c(2+) at the catalytic site and, in addition, an interaction at a second site which further modulates the reaction of cyt c(2+) at the catalytic site. These results are generally consistent with previous reports, indicating the reliability of the new procedure. We also find that a 6×His-tag at the C-terminus of the subunit II of CcO affects the binding of cyt c at both sites. The approach presented here should be generally useful in spectrophotometric studies of complex enzyme kinetics. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).

Exploring the proton pump and exit pathway for pumped protons in cytochrome ba3 from Thermus thermophilus

The heme-copper oxygen reductases are redox-driven proton pumps. In the current work, the effects of mutations in a proposed exit pathway for pumped protons are examined in the ba(3)-type oxygen reductase from Thermus thermophilus, leading from the propionates of heme a(3) to the interface between subunits I and II. Recent studies have proposed important roles for His376 and Asp372, both of which are hydrogen-bonded to propionate-A of heme a(3), and for Glu126(II) (subunit II), which is hydrogen-bonded to His376. Based on the current results, His376, Glu126(II), and Asp372 are not essential for either oxidase activity or proton pumping. In addition, Tyr133, which is hydrogen-bonded to propionate-D of heme a(3), was also shown not to be essential for function. However, two mutations of the residues hydrogen-bonded to propionate-A, Asp372Ile and His376Asn, retain high electron transfer activity and normal spectral features but, in different preparations, either do not pump protons or exhibit substantially diminished proton pumping. It is concluded that either propionate-A of heme a(3) or possibly the cluster of groups centered about the conserved water molecule that hydrogen-bonds to both propionates-A and -D of heme a(3) is a good candidate to be the proton loading site.

Gating and regulation of the cytochrome c oxidase proton pump

As a consumer of 95% of the oxygen we breathe, cytochrome c oxidase plays a major role in the energy balance of the cell. Regulation of its oxygen reduction and proton pumping activity is therefore critical to physiological function in health and disease. The location and structure of pathways for protons that are required to support cytochrome c oxidase activity are still under debate, with respect to their requirements for key residues and fixed waters, and how they are gated to prevent (or allow) proton backflow. Recent high resolution structures of bacterial and mammalian forms reveal conserved lipid and steroid binding sites as well as redox-linked conformational changes that provide new insights into potential regulatory ligands and gating modes. Mechanistic interpretation of these findings and their significance for understanding energy regulation is discussed.

Coupled electron and proton transfer reactions during the O→E transition in bovine cytochrome c oxidase

A combined DFT/electrostatic approach is employed to study the coupling of proton and electron transfer reactions in cytochrome c oxidase (CcO) and its proton pumping mechanism. The coupling of the chemical proton to the internal electron transfer within the binuclear center is examined for the O→E transition. The novel features of the His291 pumping model are proposed, which involve timely well-synchronized sequence of the proton-coupled electron transfer reactions. The obtained pK(a)s and E(m)s of the key ionizable and redox-active groups at the different stages of the O→E transition are consistent with available experimental data. The PT step from E242 to H291 is examined in detail for various redox states of the hemes and various conformations of E242 side-chain. Redox potential calculations of the successive steps in the reaction cycle during the O→E transition are able to explain a cascade of equilibria between the different intermediate states and electron redistribution between the metal centers during the course of the catalytic activity. All four electrometric phases are discussed in the light of the obtained results, providing a robust support for the His291 model of proton pumping in CcO.

Adaptation of aerobic respiration to low O2 environments

Aerobic respiration in bacteria, Archaea, and mitochondria is performed by oxygen reductase members of the heme-copper oxidoreductase superfamily. These enzymes are redox-driven proton pumps which conserve part of the free energy released from oxygen reduction to generate a proton motive force. The oxygen reductases can be divided into three main families based on evolutionary and structural analyses (A-, B- and C-families), with the B- and C-families evolving after the A-family. The A-family utilizes two proton input channels to transfer protons for pumping and chemistry, whereas the B- and C-families require only one. Generally, the B- and C-families also have higher apparent oxygen affinities than the A-family. Here we use whole cell proton pumping measurements to demonstrate differential proton pumping efficiencies between representatives of the A-, B-, and C-oxygen reductase families. The A-family has a coupling stoichiometry of 1 H(+)/e(-), whereas the B- and C-families have coupling stoichiometries of 0.5 H(+)/e(-). The differential proton pumping stoichiometries, along with differences in the structures of the proton-conducting channels, place critical constraints on models of the mechanism of proton pumping. Most significantly, it is proposed that the adaptation of aerobic respiration to low oxygen environments resulted in a concomitant reduction in energy conservation efficiency, with important physiological and ecological consequences.

Alternative initial proton acceptors for the D pathway of Rhodobacter sphaeroides cytochrome c oxidase

To characterize protein structures that control proton uptake, we assayed forms of cytochrome c oxidase (CcO) containing a carboxyl or a thiol group in line with the initial, internal waters of the D pathway for proton transfer in the presence and absence of subunit III. Subunit III provides approximately half of the protein surrounding the entry region of the D pathway. The N139D/D132N mutant contains a carboxyl group 6 Å within the D pathway and lacks the normal, surface-exposed proton acceptor, Asp-132. With subunit III, the steady-state activity of this mutant is slow, but once subunit III is removed, its activity is the same as that of wild-type CcO lacking subunit III (∼1800 H+/s). Thus, a carboxyl group∼25% within the pathway enhances proton uptake even though the carboxyl has no direct contact with bulk solvent. Protons from solvent apparently move to internal Asp-139 through a short file of waters, normally blocked by subunit III. Cys-139 also supports rapid steady-state proton uptake, demonstrating that an anion other than a carboxyl can attract and transfer protons into the D pathway. When both Asp-132 and Asp/Cys-139 are present, the removal of subunit III increases CcO activity to rates greater than that of normal CcO because of simultaneous proton uptake by two initial acceptors. The results show how the environment of the initial proton acceptor for the D pathway in these CcO forms dictates the pH range of CcO activity, with implications for the function of Asp-132, the normal proton acceptor.

Electron transfer pathways in cytochrome c oxidase

Mixed quantum mechanical/molecular mechanics calculations were used to explore the electron pathway of the terminal electron transfer enzyme, cytochrome c oxidase. This enzyme catalyzes the reduction of molecular oxygen to water in a multiple step process. Density functional calculations on the three redox centers allowed for the characterization of the electron transfer mechanism, following the sequence Cu(A)→heme a→heme a(3). This process is largely affected by the presence of positive charges, confirming the possibility of a proton coupled electron transfer. An extensive mapping of all residues involved in the electron transfer, between the Cu(A) center (donor) and the O(2) reduction site heme a(3)-Cu(B) (receptor), was obtained by selectively activating/deactivating different quantum regions. The method employed, called QM/MM e-pathway, allowed the identification of key residues along the possible electron transfer paths, consistent with experimental data. In particular, the role of arginines 481 and 482 appears crucial in the Cu(A)→heme a and in the heme a→heme a(3) electron transfer processes. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Crystallographic and online spectral evidence for role of conformational change and conserved water in cytochrome oxidase proton pump

Crystal structures in both oxidized and reduced forms are reported for two bacterial cytochrome c oxidase mutants that define the D and K proton paths, showing conformational change in response to reduction and the loss of strategic waters that can account for inhibition of proton transfer. In the oxidized state both mutants of the Rhodobacter sphaeroides enzyme, D132A and K362M, show overall structures similar to wild type, indicating no long-range effects of mutation. In the reduced state, the mutants show an altered conformation similar to that seen in reduced wild type, confirming this reproducible, reversible response to reduction. In the strongly inhibited D132A mutant, positions of residues and waters in the D pathway are unaffected except in the entry region close to the mutation, where a chloride ion replaces the missing carboxyl and a 2-Å shift in N207 results in loss of its associated water. In K362M, the methionine occupies the same position as the original lysine, but K362- and T359-associated waters in the wild-type structure are missing, likely accounting for the severe inhibition. Spectra of oxidized frozen crystals taken during X-ray radiation show metal center reduction, but indicate development of a strained configuration that only relaxes to a native form upon annealing. Resistance of the frozen crystal to structural change clarifies why the oxidized conformation is observable and supports the conclusion that the reduced conformation has functional significance. A mechanism is described that explains the conformational change and the incomplete response of the D-path mutant.

Similarity of cytochrome c oxidases in different organisms

Most of biological oxygen reduction is catalyzed by the heme-copper oxygen reductases. These enzymes are redox-driven proton pumps that take part in generating the proton gradient in both prokaryotes and mitochondria that drives synthesis of ATP. The enzymes have been divided into three evolutionarily-related groups: the A-, B-, and C-families. Recent comparative studies suggest that all oxygen reductases perform the same chemistry for oxygen reduction and comprise the same essential elements of the proton pumping mechanism, such as the proton loading and kinetic gating sites, which, however, appear to be different in different families. All species of the A-family, however, demonstrate remarkable similarity of the central processing unit of the enzyme, as revealed by their recent crystal structures. Here we demonstrate that cytochrome c oxidases (CcO) of such diverse organisms as a mammal (bovine heart mitochondrial CcO), photosynthetic bacteria (Rhodobacter sphaeroides CcO), and soil bacteria (Paracoccus denitrificans CcO) are not only structurally similar, but almost identical in microscopic electrostatics and thermodynamics properties of their key amino-acids. By using pK(a) calculations of some of the key residues of the catalytic site, D- and K- proton input, and putative proton output channels of these three different enzymes, we demonstrate that the microscopic properties of key residues are almost identical, which strongly suggests the same mechanism in these species. The quantitative precision with which the microscopic physical properties of these enzymes have remained constant despite different evolutionary routes undertaken is striking.