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.     

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₃, 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.     

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.     

Water-gated mechanism of proton translocation by cytochrome c oxidase

Cytochrome c oxidase is essential for aerobic life as a membrane-bound energy transducer. O(2) reduction at the haem a(3)-Cu(B) centre consumes electrons transferred via haem a from cytochrome c outside the membrane. Protons are taken up from the inside, both to form water and to be pumped across the membrane (M.K.F. Wikstr?m, Nature 266 (1977) 271; M. Wikstr?m, K. Krab, M. Saraste, Cytochrome Oxidase, A Synthesis, Academic Press, London, 1981 ). The resulting electrochemical proton gradient drives ATP synthesis (P. Mitchell, Chemiosmotic Coupling in Oxidative and Photosynthetic Phosphorylation, Glynn Research, Bodmin, UK, 1966 ). Here we present a molecular mechanism for proton pumping coupled to oxygen reduction that is based on the unique properties of water in hydrophobic cavities. An array of water molecules conducts protons from a conserved glutamic acid, either to the Delta-propionate of haem a(3) (pumping), or to haem a(3)-Cu(B) (water formation). Switching between these pathways is controlled by the redox-state-dependent electric field between haem a and haem a(3)-Cu(B), which determines the water-dipole orientation, and therefore the proton transfer direction. Proton transfer via the propionate provides a gate to O(2) reduction. This pumping mechanism explains the unique arrangement of the metal cofactors in the structure. It is consistent with the large body of biochemical data, and is shown to be plausible by molecular dynamics simulations.     

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.     

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₃.     

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.     

13th International Congress of Biochemistry

Cytochrome caa₃ (cytochrome oxidase) from the thermophilic bacterium PS3 can exhibit full catalytic activity in the presence of ascorbate and TMPD or other electron donors and in the absence of added soluble c-type cytochromes. It appears to possess only a low-affinity and not a high-affinity site for the soluble cytochromes. Proteoliposomal cytochrome caa₃ develops an effective membrane potential in the presence of ascorbate and TMPD or PMS, in the absence of added soluble cytochrome c. Reduction of the a₃ centre is blocked in the presence of cyanide. During reductive titrations of the cyanide-inhibited enzyme, electrons initially equilibrate among three centres, the c haem, the a haem and one of the associated Cu atoms. During steady-state turnover, electrons probably enter the complex via the bound c haem; the a haem and perhaps an associated Cu_A atom are reduced next. It is concluded that, despite its size and hydrophobic association with the aa₃ complex, the haem c-containing subunit can behave in an analogous way to that of mammalian cytochrome c, bound at the high-affinity site of the eucaryotic enzyme.     

Cytochrome oxidase: pathways for electron tunneling and proton transfer

 Electrons from cytochrome c, the substrate of cytochrome oxidase, a redox-linked proton pump, are accepted by Cu_A in subunit II. From there they are transferred to the proton pumping machinery in subunit I, cytochrome a and cytochrome a ₃–Cu_B. The reduction of the latter site, which is the dioxygen reducing unit, is coupled to proton uptake. Dioxygen reduction involves a peroxide and a ferryl ion intermediate, and it is the transition between these and back to the resting oxidized enzyme that are coupled to proton pumping. The X-ray structures suggest electron–transfer pathways that can account for the observed rates provided that the reorganization energies are small. They also reveal two proton-transfer pathways, and mutagenesis experiments have shown that one is used for proton uptake during the initial reduction of cytochrome a ₃–Cu_B, whereas the other mediates transfer of the pumped protons. Received: 23 March 1998 / Accepted: 11 May 1998     

Critical roles of the CuB site in efficient proton pumping as revealed by crystal structures of mammalian cytochrome c oxidase catalytic intermediates

Mammalian cytochrome c oxidase (CcO) reduces O₂ to water in a bimetallic site including Fe_a3 and Cu_B giving intermediate molecules, termed A-, P-, F-, O-, E-, and R-forms. From the P-form on, each reaction step is driven by single-electron donations from cytochrome c coupled with the pumping of a single proton through the H-pathway, a proton-conducting pathway composed of a hydrogen-bond network and a water channel. The proton-gradient formed is utilized for ATP production by F-ATPase. For elucidation of the proton pumping mechanism, crystal structural determination of these intermediate forms is necessary. Here we report X-ray crystallographic analysis at ∼1.8 Å resolution of fully reduced CcO crystals treated with O₂ for three different time periods. Our disentanglement of intermediate forms from crystals that were composed of multiple forms determined that these three crystallographic data sets contained ∼45% of the O-form structure, ∼45% of the E-form structure, and ∼20% of an oxymyoglobin-type structure consistent with the A-form, respectively. The O- and E-forms exhibit an unusually long Cu_B²⁺-OH⁻ distance and Cu_B¹⁺-H₂O structure keeping Fe_a3³⁺-OH⁻ state, respectively, suggesting that the O- and E-forms have high electron affinities that cause the O→E and E→R transitions to be essentially irreversible and thus enable tightly coupled proton pumping. The water channel of the H-pathway is closed in the O- and E-forms and partially open in the R-form. These structures, together with those of the recently reported P- and F-forms, indicate that closure of the H-pathway water channel avoids back-leaking of protons for facilitating the effective proton pumping.     

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.     

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.     

The H channel is not a proton transfer path in yeast cytochrome c oxidase

Cytochrome c oxidases (CcOs) in the respiratory chains of mitochondria and bacteria are primary consumers of molecular oxygen, converting it to water with the concomitant pumping of protons across the membrane to establish a proton electrochemical gradient. Despite a relatively well understood proton pumping mechanism of bacterial CcOs, the role of the H channel in mitochondrial forms of CcO remains debated. Here, we used site-directed mutagenesis to modify a central residue of the lower span of the H channel, Q413, in the genetically tractable yeast Saccharomyces cerevisiae. Exchange of Q413 to several different amino acids showed no effect on rates and efficiencies of respiratory cell growth, and redox potential measurements indicated minimal electrostatic interaction between the 413 locus and the nearest redox active component heme a. These findings clearly exclude a primary role of this section of the H channel in proton pumping in yeast CcO. In agreement with the experimental data, atomistic molecular dynamics simulations and continuum electrostatic calculations on wildtype and mutant yeast CcOs highlight potential bottlenecks in proton transfer through this route. Our data highlight the preference for neutral residues in the 413 locus, precluding sufficient hydration for formation of a proton conducting wire.     

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 a₃ (PRDa₃) for the first time. The results from these molecular dynamics simulations indicate that the PRDa₃ 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 PRDa₃ side chain isomerization and the proton translocation between Glu242 and the PRDa₃ 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 PRDa₃ site.     

Mechanisms for enzymatic reduction of nitric oxide to nitrous oxide - A comparison between nitric oxide reductase and cytochrome c oxidase

Cytochrome c oxidases (CcO) reduce O₂ to H₂O in the respiratory chain of mitochondria and many aerobic bacteria. In addition, some species of CcO can also reduce NO to N₂O and water while others cannot. Here, the mechanism for NO-reduction in CcO is investigated using quantum mechanical calculations. Comparison is made to the corresponding reaction in a “true” cytochrome c-dependent NO reductase (cNOR). The calculations show that in cNOR, where the reduction potentials are low, the toxic NO molecules are rapidly reduced, while the higher reduction potentials in CcO lead to a slower or even impossible reaction, consistent with experimental observations. In both enzymes the reaction is initiated by addition of two NO molecules to the reduced active site, forming a hyponitrite intermediate. In cNOR, N₂O can then be formed using only the active-site electrons. In contrast, in CcO, one proton-coupled reduction step most likely has to occur before N₂O can be formed, and furthermore, proton transfer is most likely rate-limiting. This can explain why different CcO species with the same heme a₃-Cu active site differ with respect to NO reduction efficiency, since they have a varying number and/or properties of proton channels. Finally, the calculations also indicate that a conserved active site valine plays a role in reducing the rate of NO reduction in CcO.     

Studies of cytochrome synthesis in rat liver

The incorporation of radioactive amino acids and of δ-amino[2,3-³H₂]laevulinate into rat liver cytochromes b₅ and c and cytochrome oxidase has been examined with and without protein-synthesis inhibitors. Cycloheximide promptly inhibits labelling of both haem and protein for cytochrome c in parallel fashion. Although incorporation of ¹⁴C-labelled amino acid into microsomal cytochrome b₅ is also rapidly inhibited, cycloheximide incompletely inhibits haem labelling of cytochrome b₅ and cytochrome a+a₃, and inhibition occurs only after repeated antibiotic injections. The possibility of apo-protein pools, or of haem exchange, with a rapidly renewed `free' haem pool, is considered. Consistent with this model is the observation of non-enzymic haem exchange in vitro between cytochrome b₅ and methaemoglobin. Chloramphenicol, injected intravenously over 5h, results in a 20–40% decrease in incorporation of δ-amino[2,3-³H₂]laevulinate into haem a+a₃ and haem of cytochromes b₅ and c. With the dosage schedule of chloramphenicol studied, amino acid labelling of total liver protein and of cytochrome c was not inhibited. Similarly, ferrochelatase activity was not decreased.     

Cytochrome a-type terminal oxidase derived from Thiobacillus novellus. Molecular and enzymatic properties

Cytochrome a-type terminal oxidase was purified from Thiobacillus novellus to an electrophoretically homogeneous state and some of its properties were studied.The enzyme shows absorption peaks at 428 and 602 nm in the oxidized form, and at 442 and 602 nm in the reduced form. The CO compound of the reduced enzyme shows peaks at 431 and 599 nm. The enzyme has 1 mol of haem a and 1 g-atom of copper per 55 600 g and is composed of two kinds of subunit, of 32 000 and 23 000 daltons, respectively.The enzyme reacts rapidly with tuna, bonito and yeast cytochromes c as well as with T. novellus cytochrome c, while it reacts slowly with horse and cow cytochromes c. The reduction product of oxygen catalysed by the enzyme is water.     

The sequence of electron carriers in the reaction of cytochromec oxidase with oxygen

Kinetic studies of the electron transfer processes performed by cytochrome oxidase have assigned rates of electron transfer between the metal centers involved in the oxidation of ferrocytochromec by molecular oxygen. Transient-state studies of the reaction with oxygen have led to the proposal of a sequence of carriers from cytochromec, to Cu_A, to cytochromea, and then to the binuclear (i.e., cytochromea ₃-Cu_B) center. Electron exchange rates between these centers agree with relative center-to-center distances as follows; cytochromec to Cu_A 5–7 Å, cytochromec to cytochromea 20–25 Å, Cu_A to cytochromea 14–16 Å and cytochromea to cytochrome a₃-Cu_B 8–10 Å. It is proposed that the step from cytochromea to the binuclear center is the key control point in the reaction and that this step is one of the major points of energy transduction in the reaction cycle.