Redox-coupled proton pumping in cytochrome c oxidase: further insights from computer simulation

The membrane-bound enzyme cytochrome c oxidase, the terminal member in the respiratory chain, converts oxygen into water and generates an electrochemical gradient by coupling the electron transfer to proton pumping across the membrane. Here we have investigated the dynamics of an excess proton and the surrounding protein environment near the active sites. The multi-state empirical valence bond (MS-EVB) molecular dynamics method was used to simulate the explicit dynamics of proton transfer through the critically important hydrophobic channel between Glu242 (bovine notation) and the D-propionate of heme a3 (PRDa3) for the first time. The results from these molecular dynamics simulations indicate that the PRDa3 can indeed re-orientate and dissociate from Arg438, despite the high stability of such an ion pair, and has the ability to accept protons via bound water molecules. Any large conformational change of the adjacent heme a D-propionate group is, however, sterically blocked directly by the protein. Free energy calculations of the PRDa3 side chain isomerization and the proton translocation between Glu242 and the PRDa3 site have also been performed. The results exhibit a redox state-dependent dynamical behavior and indicate that reduction of the low-spin heme a may initiate internal transfer of the pumped proton from Glu242 to the PRDa3 site.     

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.     

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.     

Dynamic water networks in cytochrome cbb3 oxidase

Heme-copper oxidases (HCOs) are terminal electron acceptors in aerobic respiration. They catalyze the reduction of molecular oxygen to water with concurrent pumping of protons across the mitochondrial and bacterial membranes. Protons required for oxygen reduction chemistry and pumping are transferred through proton uptake channels. Recently, the crystal structure of the first C-type member of the HCO superfamily was resolved [Buschmann et al. Science 329 (2010) 327–330], but crystallographic water molecules could not be identified. Here we have used molecular dynamics (MD) simulations, continuum electrostatic approaches, and quantum chemical cluster calculations to identify proton transfer pathways in cytochrome cbb₃. In MD simulations we observe formation of stable water chains that connect the highly conserved Glu323 residue on the proximal side of heme b₃ both with the N- and the P-sides of the membrane. We propose that such pathways could be utilized for redox-coupled proton pumping in the C-type oxidases. Electrostatics and quantum chemical calculations suggest an increased proton affinity of Glu323 upon reduction of high-spin heme b₃. Protonation of Glu323 provides a mechanism to tune the redox potential of heme b₃ with possible implications for proton pumping.     

Effective Pumping Proton Collection Facilitated by a Copper Site (CuB) of Bovine Heart Cytochrome c Oxidase,Revealed by a Newly Developed Time-resolved Infrared System

X-ray structural and mutational analyses have shown that bovine heart cytochrome c oxidase (CcO) pumps protons electrostatically through a hydrogen bond network using net positive charges created upon oxidation of a heme iron (located near the hydrogen bond network) for O₂ reduction. Pumping protons are transferred by mobile water molecules from the negative side of the mitochondrial inner membrane through a water channel into the hydrogen bond network. For blockage of spontaneous proton back-leak, the water channel is closed upon O₂ binding to the second heme (heme a₃) after complete collection of the pumping protons in the hydrogen bond network. For elucidation of the structural bases for the mechanism of the proton collection and timely closure of the water channel, conformational dynamics after photolysis of CO (an O₂ analog)-bound CcO was examined using a newly developed time-resolved infrared system feasible for accurate detection of a single C=O stretch band of α-helices of CcO in H₂O medium. The present results indicate that migration of CO from heme a₃ to Cu_B in the O₂ reduction site induces an intermediate state in which a bulge conformation at Ser-382 in a transmembrane helix is eliminated to open the water channel. The structural changes suggest that, using a conformational relay system, including Cu_B, O₂, heme a₃, and two helix turns extending to Ser-382, Cu_B induces the conformational changes of the water channel that stimulate the proton collection, and senses complete proton loading into the hydrogen bond network to trigger the timely channel closure by O₂ transfer from Cu_B to heme a₃.     

Hydrogen-Bonded Networks Along and Bifurcation of the E-Pathway in Quinol:Fumarate Reductase

The E-pathway of transmembrane proton transfer has been demonstrated previously to be essential for catalysis by the diheme-containing quinol:fumarate reductase (QFR) of Wolinella succinogenes. Two constituents of this pathway, Glu-C180 and heme b_D ring C (b_D-C-) propionate, have been validated experimentally. Here, we identify further constituents of the E-pathway by analysis of molecular dynamics simulations. The redox state of heme groups has a crucial effect on the connectivity patterns of mobile internal water molecules that can transiently support proton transfer from the b_D-C-propionate to Glu-C180. The short H-bonding paths formed in the reduced states can lead to high proton conduction rates and thus provide a plausible explanation for the required opening of the E-pathway in reduced QFR. We found evidence that the b_D-C-propionate group is the previously postulated branching point connecting proton transfer to the E-pathway from the quinol-oxidation site via interactions with the heme b_D ligand His-C44. An essential functional role of His-C44 is supported experimentally by site-directed mutagenesis resulting in its replacement with Glu. Although the H44E variant enzyme retains both heme groups, it is unable to catalyze quinol oxidation. All results obtained are relevant to the QFR enzymes from the human pathogens Campylobacter jejuni and Helicobacter pylori.     

Redox-coupled proton transfer in the active site of cytochrome cbb3

Cytochrome cbb₃ is a distinct member of the superfamily of respiratory heme-copper oxidases, and is responsible for driving the respiratory chain in many pathogenic bacteria. Like the canonical heme-copper oxidases, cytochrome cbb₃ reduces oxygen to water and couples the released energy to pump protons across the bacterial membrane. Homology modeling and recent electron paramagnetic resonance (EPR) studies on wild type and a mutant cbb₃ enzyme [V. Rauhamäki et al. J. Biol. Chem. 284 (2009) 11301-11308] have led us to perform high-level quantum chemical calculations on the active site. These calculations bring molecular insight into the unique hydrogen bonding between the proximal histidine ligand of heme b₃ and a conserved glutamate, and indicate that the catalytic mechanism involves redox-coupled proton transfer between these residues. The calculated spin densities give insight in the difference in EPR spectra for the wild type and a recently studied E383Q-mutant cbb₃-enzyme. Furthermore, we show that the redox-coupled proton movement in the proximal cavity of cbb₃-enzymes contributes to the low redox potential of heme b₃, and suggest its potential implications for the high apparent oxygen affinity of these enzymes.     

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(3)/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(2) 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.     

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.     

Two conformational states of Glu242 and pKas in bovine cytochrome c oxidase.

Cytochrome c oxidase (CcO) is the terminal enzyme in the respiratory electron transport chain of aerobic organisms. It catalyses the reduction of atmospheric oxygen to water, and couples this reaction to proton pumping across the membrane; this process generates the electrochemical gradient that subsequently drives the synthesis of ATP. The molecular details of the mechanism by which electron transfer is coupled to proton pumping in CcO is poorly understood. Recent calculations from our group indicate that His291, a ligand of the Cu(B) center of the enzyme, may play the role of the pumping element. In this paper we describe calculations in which a DFT/continuum electrostatic method is used to explore the coupling of the conformational changes of Glu242 residue, the main proton donor of both chemical and pump protons, to its pKa, and the pKa of His291, a putative proton loading site of our pumping model. The computations are done for several redox states of metal centers, different protonation states of Glu242 and His291, and two well-defined conformations of the Glu242 side chain. Thus, in addition to equilibrium redox/protonation states of the catalytic cycle, we also examine the transient and intermediate states. Different dielectric models are employed to investigate the robustness of the results, and their viability in the light of the proposed proton pumping mechanism of CcO. The main results are in agreement with the experimental measurements and support the proposed pumping mechanism. Additionally, the present calculations indicate a possibility of gating through conformational changes of Glu242; namely, in the pumping step, we find that Glu242 needs to be reprotonated before His291 can eject a proton to the P-site of membrane. As a result, the reprotonation of Glu can control proton release from the proton loading site.     

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 study of the proposed proton transport pathways in [FeFe]-hydrogenase

Possible proton transport pathways in Clostridium pasteurianum (CpI) [FeFe]-hydrogenase were investigated with molecular dynamics simulations. This study was undertaken to evaluate the functional pathway and provide insight into the hydrogen bonding features defining an active proton transport pathway. Three pathways were evaluated, two of which consist of water wires and one of predominantly amino acid residues. Our simulations suggest that protons are not transported through water wires. Instead, the five-residue motif (Glu282, Ser319, Glu279, H₂O, Cys299) was found to be the likely pathway, consistent with previously made experimental observations. The pathway was found to have a persistent hydrogen bonded core (residues Cys299 to Ser319), with less persistent hydrogen bonds at the ends of the pathway for both H₂ release and H₂ uptake. Single site mutations of the four residues have been shown experimentally to deactivate the enzyme. The theoretical evaluation of these mutations demonstrates redistribution of the hydrogen bonds in the pathway, resulting in enzyme deactivation. Finally, coupling between the protein dynamics near the proton transport pathway and the redox partner binding regions was also found as a function of H₂ uptake and H₂ release states, which may be indicative of a correlation between proton and electron movement within the enzyme.     

The Protein Effect in the Structure of Two Ferryl-Oxo Intermediates at the Same Oxidation Level in the Heme Copper Binuclear Center of Cytochrome c Oxidase

Identification of the intermediates and determination of their structures in the reduction of dioxygen to water by cytochrome c oxidase (CcO) are particularly important to understanding both O₂ activation and proton pumping by the enzyme. In this work, we report the products of the rapid reaction of O₂ with the mixed valence form (Cu_A²⁺, heme a³⁺, heme a₃²⁺-Cu_B¹⁺) of the enzyme. The resonance Raman results show the formation of two ferryl-oxo species with characteristic Fe(IV)=O stretching modes at 790 and 804 cm⁻¹ at the peroxy oxidation level (P_M). Density functional theory calculations show that the protein environment of the proximal H-bonded His-411 determines the strength of the distal Fe(IV)=O bond. In contrast to previous proposals, the P_M intermediate is also formed in the reaction of Y167F with O₂. These results suggest that in the fully reduced enzyme, the proton pumping ν_Fe(IV)=O = 804 cm⁻¹ to ν_Fe(IV)=O = 790 cm⁻¹ transition (P→F, where P is peroxy and F is ferryl) is triggered not only by electron transfer from heme a to heme a₃ but also by the formation of the H-bonded form of the His-411-Fe(IV)=O conformer in the proximal site of heme a₃. The implications of these results with respect to the role of an O=Fe(IV)-His-411-H-bonded form to the ring A propionate of heme a₃-Asp-399-H₂O site and, thus, to the exit/output proton channel (H₂O) pool during the proton pumping P→F transition are discussed. We propose that the environment proximal to the heme a₃ controls the spectroscopic properties of the ferryl intermediates in cytochrome oxidases.     

Concerted involvement of cooperative proton–electron linkage and water production in the proton pump of cytochrome c oxidase

In cytochrome c oxidase, oxido-reductions of heme a/Cu_A and heme a₃/Cu_B are cooperatively linked to proton transfer at acid/base groups in the enzyme. H⁺/e⁻ cooperative linkage at Fe_a3/Cu_B is envisaged to be involved in proton pump mechanisms confined to the binuclear center. Models have also been proposed which involve a role in proton pumping of cooperative H⁺/e⁻ linkage at heme a (and Cu_A). Observations will be presented on: (i) proton consumption in the reduction of molecular oxygen to H₂O in soluble bovine heart cytochrome c oxidase; (ii) proton release/uptake associated with anaerobic oxidation/reduction of heme a/Cu_A and heme a₃/Cu_B in the soluble oxidase; (iii) H⁺ release in the external phase (i.e. H⁺ pumping) associated with the oxidative (R → O transition), reductive (O → R transition) and a full catalytic cycle (R → O → R transition) of membrane-reconstituted cytochrome c oxidase. A model is presented in which cooperative H⁺/e⁻ linkage at heme a/Cu_A and heme a₃/Cu_B with acid/base clusters, C₁ and C₂ respectively, and protonmotive steps of the reduction of O₂ to water are involved in proton pumping.     

Variable proton-pumping stoichiometry in structural variants of cytochrome c oxidase

Cytochrome c oxidase is a multisubunit membrane-bound enzyme, which catalyzes oxidation of four molecules of cytochrome c²⁺ and reduction of molecular oxygen to water. The electrons are taken from one side of the membrane while the protons are taken from the other side. This topographical arrangement results in a charge separation that is equivalent to moving one positive charge across the membrane for each electron transferred to O₂. In this reaction part of the free energy available from O₂ reduction is conserved in the form of an electrochemical proton gradient. In addition, part of the free energy is used to pump on average one proton across the membrane per electron transferred to O₂. Our understanding of the molecular design of the machinery that couples O₂ reduction to proton pumping in oxidases has greatly benefited from studies of so called “uncoupled” structural variants of the oxidases. In these uncoupled oxidases the catalytic O₂-reduction reaction may display the same rates as in the wild-type CytcO, yet the electron/proton transfer to O₂ is not linked to proton pumping. One striking feature of all uncoupled variants studied to date is that the (apparent) pK_a of a Glu residue, located deeply within a proton pathway, is either increased or decreased (from 9.4 in the wild-type oxidase). The altered pK_a presumably reflects changes in the local structural environment of the residue and because the Glu residue is found near the catalytic site as well as near a putative exit pathway for pumped protons these changes are presumably important for controlling the rates and trajectories of the proton transfer. In this paper we summarize data obtained from studies of uncoupled structural oxidase variants and present a hypothesis that in quantitative terms offers a link between structural changes, modulation of the apparent pK_a and uncoupling of proton pumping from O₂ reduction.     

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.     

Resolving the Negative Potential Side (n-side) Water-accessible Proton Pathway of F-type ATP Synthase by Molecular Dynamics Simulations

The rotation of F₁F_o-ATP synthase is powered by the proton motive force across the energy-transducing membrane. The protein complex functions like a turbine; the proton flow drives the rotation of the c-ring of the transmembrane F_o domain, which is coupled to the ATP-producing F₁ domain. The hairpin-structured c-protomers transport the protons by reversible protonation/deprotonation of a conserved Asp/Glu at the outer transmembrane helix (TMH). An open question is the proton transfer pathway through the membrane at atomic resolution. The protons are thought to be transferred via two half-channels to and from the conserved cAsp/Glu in the middle of the membrane. By molecular dynamics simulations of c-ring structures in a lipid bilayer, we mapped a water channel as one of the half-channels. We also analyzed the suppressor mutant cP24D/E61G in which the functional carboxylate is shifted to the inner TMH of the c-protomers. Current models concentrating on the “locked” and “open” conformations of the conserved carboxylate side chain are unable to explain the molecular function of this mutant. Our molecular dynamics simulations revealed an extended water channel with additional water molecules bridging the distance of the outer to the inner TMH. We suggest that the geometry of the water channel is an important feature for the molecular function of the membrane part of F₁F_o-ATP synthase. The inclination of the proton pathway isolates the two half-channels and may contribute to a favorable clockwise rotation in ATP synthesis mode.