Reply to “Comment on Acidity of a Cu-bound Histidine in the Binuclear Center of Cytochrome c Oxidase”

The mechanism for proton pumping in cytochrome c oxidase from an electrostatic and quantum chemical perspective

The mechanism for proton pumping in cytochrome c oxidase in the respiratory chain, has for decades been one of the main unsolved problems in biochemistry. However, even though several different suggested mechanisms exist, many of the steps in these mechanisms are quite similar and constitute a general consensus framework for discussing proton pumping. When these steps are analyzed, at least three critical gating situations are found, and these points are where the suggested mechanisms in general differ. The requirements for gating are reviewed and analyzed in detail, and a mechanism is suggested, where solutions for all the gating situations are formulated. This mechanism is based on an electrostatic analysis of a kinetic experiment fior the O to E transition. The key component of the mechanism is a positively charged transition state. An electron on heme a opens the gate for proton transfer from the N-side to a pump loading site (PLS). When the negative charge of the electron is compensated by a chemical proton, the positive transition state prevents backflow from the PLS to the N-side at the most critical stage of the pumping process. The mechanism has now been tested by large model DFT calculations, and these calculations give strong support for the suggested mechanism.

B3LYP study on reduction mechanisms from O2 to H2O at the catalytic sites of fully reduced and mixed-valence bovine cytochrome c oxidases

Reduction mechanisms of oxygen molecule to water molecules in the fully reduced (FR) and mixed-valence (MV) bovine cytochrome c oxidases (CcO) have been systematically examined based on the B3LYP calculations. The catalytic cycle using four electrons and four protons has been also shown consistently. The MV CcO catalyses reduction to produce one water molecule, while the FR CcO catalyses to produce two water molecules. One water molecule is added into vacant space between His240 and His290 in the catalytic site. This water molecule constructs the network of hydrogen bonds of Tyr244, farnesyl ethyl, and Thr316 that is a terminal residue of the K-pathway. It plays crucial roles for the proton transfer to the dioxygen to produce the water molecules in both MV and FR CcOs. Tyr244 functions as a relay of the proton transfer from the K-pathway to the added water molecule, not as donors of a proton and an electron to the dioxygen. The reduction mechanisms of MV and FR CcOs are strictly distinguished. In the FR CcO, the Cu atom at the Cu(B) site maintains the reduced state Cu(I) during the process of formation of first water molecule and plays an electron storage. At the final stage of formation of first water molecule, the Cu(I) atom releases an electron to Fe-O. During the process of formation of second water molecule, the Cu atom maintains the oxidized state Cu(II). In contrast with experimental proposals, the K-pathway functions for formation of first water molecule, while the D-pathway functions for second water molecule. The intermediates, P(M), P(R), F, and O, obtained in this work are compared with those proposed experimentally.

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

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

Toward a chemical mechanism of proton pumping by the B-type cytochrome c oxidases: application of density functional theory to cytochrome ba3 of Thermus thermophilus

A mechanism for proton pumping by the B-type cytochrome c oxidases is presented in which one proton is pumped in conjunction with the weakly exergonic, two-electron reduction of Fe-bound O 2 to the Fe-Cu bridging peroxodianion and three protons are pumped in conjunction with the highly exergonic, two-electron reduction of Fe(III)- (-)O-O (-)-Cu(II) to form water and the active oxidized enzyme, Fe(III)- (-)OH,Cu(II). The scheme is based on the active-site structure of cytochrome ba 3 from Thermus thermophilus, which is considered to be both necessary and sufficient for coupled O 2 reduction and proton pumping when appropriate gates are in place (not included in the model). Fourteen detailed structures obtained from density functional theory (DFT) geometry optimization are presented that are reasonably thought to occur during the four-electron reduction of O 2. Each proton-pumping step takes place when a proton resides on the imidazole ring of I-His376 and the large active-site cluster has a net charge of +1 due to an uncompensated, positive charge formally associated with Cu B. Four types of DFT were applied to determine the energy of each intermediate, and standard thermochemical approaches were used to obtain the reaction free energies for each step in the catalytic cycle. This application of DFT generally conforms with previously suggested criteria for a valid model (Siegbahn, P. E. M.; Blomberg, M. A. R. Chem. Rev. 2000, 100, 421-437) and shows how the chemistry of O 2 reduction in the heme a 3 -Cu B dinuclear center can be harnessed to generate an electrochemical proton gradient across the lipid bilayer.