The coupling of electron transfer and proton translocation: electrostatic calculations on Paracoccus denitrificans cytochrome c oxidase

Investigation of the Mechanism of Membrane Potential Generation by Heme-Copper Respiratory Oxidases in a Real Time Mode

Heme-copper respiratory oxidases are highly efficient molecular machines. These membrane enzymes catalyze the final step of cellular respiration in eukaryotes and many prokaryotes: the transfer of electrons from cytochromes or quinols to molecular oxygen and oxygen reduction to water. The free energy released in this redox reaction is converted by heme-copper respiratory oxidases into the transmembrane gradient of the electrochemical potential of hydrogen ions H+). Heme-copper respiratory oxidases have a unique mechanism for generating H+, namely, a redox-coupled proton pump. A combination of direct electrometric method for measuring the kinetics of membrane potential generation with the methods of prestationary kinetics and site-directed mutagenesis in the studies of heme-copper oxidases allows to obtain a unique information on the translocation of protons inside the proteins in real time. The review summarizes the data of studies employing time-resolved electrometry to decipher the mechanisms of functioning of these important bioenergetic enzymes.

Effect of Ca2+ on the redox potential of heme a in cytochrome c oxidase

Subunit I of cytochrome c oxidase (CcO) from mitochondria and many bacteria contains a cation binding site (CBS) located at the outer positively charged aqueous phase not far from heme a. Binding of Ca²⁺ with the CBS in bovine CcO inhibits activity of the enzyme 2-3 -fold [Vygodina, T., Kirichenko, A. & Konstantinov A.A. (2013) Direct Regulation of Cytochrome c Oxidase by Calcium Ions, PLoS One.8 e74436]. Here we show that binding of Ca²⁺ at CBS of bovine CcO shifts E_m of heme a to the positive by 15-20 mV. Na⁺ ions that bind to the same site and compete with Ca²⁺ do not affect E_m of heme a and also prevent and reverse the effect of Ca²⁺. No effect of Ca²⁺ or EGTA is observed on E_m of heme a with the wild type bacterial oxidases from R.sphaeroides or P.denitrificans that contain tightly-bound calcium at the site. In the D477A mutant CcO from P. denitrificans that binds Ca²⁺ reversibly like the mitochondrial CcO, calcium shifts redox titration curve of heme a to the positive by ∼35-50 mV that is in good agreement with the results of electrostatic calculations; however, as shown earlier, it does not inhibit CcO activity of the mutant enzyme. Therefore the data do not support the proposal that the inhibitory effect of Ca²⁺ on CcO activity may be explained by the Ca²⁺-induced shift of E_m of heme a. Rather, Ca²⁺ retards electron transfer by inhibition of charge dislocation in the exit part of the proton channel H in mammalian CcO, that is absent in the bacterial oxidases.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286_I, Y288_I, and K362_I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.

Proton Wire Dynamics in the Green Fluorescent Protein

Inside proteins, protons move on proton wires (PWs). Starting from the highest resolution X-ray structure available, we conduct a 306 ns molecular dynamics simulation of the (A-state) wild-type (wt) green fluorescent protein (GFP) to study how its PWs change with time. We find that the PW from the chromophore via Ser205 to Glu222, observed in all X-ray structures, undergoes rapid water molecule insertion between Ser205 and Glu222. Sometimes, an alternate Ser205-bypassing PW exists. Side chain rotations of Thr203 and Ser205 play an important role in shaping the PW network in the chromophore region. Thr203, with its bulkier side chain, exhibits slower transitions between its three rotameric states. Ser205 experiences more frequent rotations, slowing down when the Thr203 methyl group is close by. The combined states of both residues affect the PW probabilities. A random walk search for PWs from the chromophore reveals several exit points to the bulk, one being a direct water wire (WW) from the chromophore to the bulk. A longer WW connects the "bottom" of the GFP barrel with a "water pool" (WP1) situated below Glu222. These two WWs were not observed in X-ray structures of wt-GFP, but their analogues have been reported in related fluorescent proteins. Surprisingly, the high-resolution X-ray structure utilized herein shows that Glu222 is protonated at low temperatures. At higher temperatures, we suggest ion pairing between anionic Glu222 and a proton hosted in WP1. Upon photoexcitation, these two recombine, while a second proton dissociates from the chromophore and either exits the protein using the short WW or migrates along the GFP-barrel axis on the long WW. This mechanism reconciles the conflicting experimental and theoretical data on proton motion within GFP.

Amazing structure of respirasome: unveiling the secrets of cell respiration

Respirasome, a huge molecular machine that carries out cellular respiration, has gained growing attention since its discovery, because respiration is the most indispensable biological process in almost all living creatures. The concept of respirasome has renewed our understanding of the respiratory chain organization, and most recently, the structure of respirasome solved by Yang's group from Tsinghua University (Gu et al. Nature 237(7622):639-643, 2016) firstly presented the detailed interactions within this huge molecular machine, and provided important information for drug design and screening. However, the study of cellular respiration went through a long history. Here, we briefly showed the detoured history of respiratory chain investigation, and then described the amazing structure of respirasome.

Molecular simulation and modeling of complex I

Molecular modeling and molecular dynamics simulations play an important role in the functional characterization of complex I. With its large size and complicated function, linking quinone reduction to proton pumping across a membrane, complex I poses unique modeling challenges. Nonetheless, simulations have already helped in the identification of possible proton transfer pathways. Simulations have also shed light on the coupling between electron and proton transfer, thus pointing the way in the search for the mechanistic principles underlying the proton pump. In addition to reviewing what has already been achieved in complex I modeling, we aim here to identify pressing issues and to provide guidance for future research to harness the power of modeling in the functional characterization of complex I. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.

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.

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.

Calculating pH-dependent free energy of proteins by using Monte Carlo protonation probabilities of ionizable residues

Protein folding, stability, and function are usually influenced by pH. And free energy plays a fundamental role in analysis of such pH-dependent properties. Electrostatics-based theoretical framework using dielectric solvent continuum model and solving Poisson-Boltzmann equation numerically has been shown to be very successful in understanding the pH-dependent properties. However, in this approach the exact computation of pH-dependent free energy becomes impractical for proteins possessing more than several tens of ionizable sites (e.g. > 30), because exact evaluation of the partition function requires a summation over a vast number of possible protonation microstates. Here we present a method which computes the free energy using the average energy and the protonation probabilities of ionizable sites obtained by the well-established Monte Carlo sampling procedure. The key feature is to calculate the entropy by using the protonation probabilities. We used this method to examine a well-studied protein (lysozyme) and produced results which agree very well with the exact calculations. Applications to the optimum pH of maximal stability of proteins and protein-DNA interactions have also resulted in good agreement with experimental data. These examples recommend our method for application to the elucidation of the pH-dependent properties of proteins.

GMCT : a Monte Carlo simulation package for macromolecular receptors

Generalized Monte Carlo titration (GMCT) is a versatile suite of computer programs for the efficient simulation of complex macromolecular receptor systems as for example proteins. The computational model of the system is based on a microstate description of the receptor and an average description of its surroundings in terms of chemical potentials. The receptor can be modeled in great detail including conformational flexibility and many binding sites with multiple different forms that can bind different ligand types. Membrane embedded systems can be modeled including electrochemical potential gradients. Overall properties of the receptor as well as properties of individual sites can be studied with a variety of different Monte Carlo (MC) simulation methods. Metropolis MC, Wang-Landau MC and efficient free energy calculation methods are included. GMCT is distributed as free open source software at www.bisb.uni-bayreuth.de under the terms of the GNU Affero General Public License.

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.

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(3). In MD simulations we observe formation of stable water chains that connect the highly conserved Glu323 residue on the proximal side of heme b(3) 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(3). Protonation of Glu323 provides a mechanism to tune the redox potential of heme b(3) with possible implications for proton pumping.

Proton-pumping mechanism of cytochrome c oxidase: a kinetic master-equation approach

Cytochrome c oxidase is an efficient energy transducer that reduces oxygen to water and converts the released chemical energy into an electrochemical membrane potential. As a true proton pump, cytochrome c oxidase translocates protons across the membrane against this potential. Based on a wealth of experiments and calculations, an increasingly detailed picture of the reaction intermediates in the redox cycle has emerged. However, the fundamental mechanism of proton pumping coupled to redox chemistry remains largely unresolved. Here we examine and extend a kinetic master-equation approach to gain insight into redox-coupled proton pumping in cytochrome c oxidase. Basic principles of the cytochrome c oxidase proton pump emerge from an analysis of the simplest kinetic models that retain essential elements of the experimentally determined structure, energetics, and kinetics, and that satisfy fundamental physical principles. The master-equation models allow us to address the question of how pumping can be achieved in a system in which all reaction steps are reversible. Whereas proton pumping does not require the direct modulation of microscopic reaction barriers, such kinetic gating greatly increases the pumping efficiency. Further efficiency gains can be achieved by partially decoupling the proton uptake pathway from the active-site region. Such a mechanism is consistent with the proposed Glu valve, in which the side chain of a key glutamic acid shuttles between the D channel and the active-site region. We also show that the models predict only small proton leaks even in the absence of turnover. The design principles identified here for cytochrome c oxidase provide a blueprint for novel biology-inspired fuel cells, and the master-equation formulation should prove useful also for other molecular machines. .

Probing protonation/deprotonation of tyrosine residues in cytochrome ba3 oxidase from Thermus thermophilus by time-resolved step-scan Fourier transform infrared spectroscopy

Elucidating the properties of the heme Fe-Cu(B) binuclear center and the dynamics of the protein response in cytochrome c oxidase is crucial to understanding not only the dioxygen activation and bond cleavage by the enzyme but also the events related to the release of the produced water molecules. The time-resolved step-scan FTIR difference spectra show the ν(7a)(CO) of the protonated form of Tyr residues at 1247 cm(-1) and that of the deprotonated form at 1301 cm(-1). By monitoring the intensity changes of the 1247 and 1301 cm(-1) modes as a function of pH, we measured a pK(a) of 7.8 for the observed tyrosine. The FTIR spectral changes associated with the tyrosine do not belong to Tyr-237 but are attributed to the highly conserved in heme-copper oxidases Tyr-136 and/or Tyr-133 residue (Koutsoupakis, K., Stavrakis, S., Pinakoulaki, E., Soulimane, T., and Varotsis, C. (2002) J. Biol. Chem. 277, 32860-32866). The oxygenation of CO by the mixed-valence form of the enzyme revealed the formation of the ∼607 nm P (Fe(IV)=O) species in the pH 6-9 range and the return to the oxidized form without the formation of the 580 nm F form. The data indicate that Tyr-237 is not involved in the proton transfer pathway in the oxygenation of CO by the mixed-valence form of the enzyme. The implication of these results with respect to the role of Tyr-136 and Tyr-133 in proton transfer/gating along with heme a(3) ring D propionate-H(2)O-ring A propionate-Asp-372 site to the exit/output proton channel (H(2)O pool) is discussed.

Redox Bohr effects and the role of heme a in the proton pump of bovine heart cytochrome c oxidase

Structural and functional observations are reviewed which provide evidence for a central role of redox Bohr effect linked to the low-spin heme a in the proton pump of bovine heart cytochrome c oxidase. Data on the membrane sidedness of Bohr protons linked to anaerobic oxido-reduction of the individual metal centers in the liposome reconstituted oxidase are analysed. Redox Bohr protons coupled to anaerobic oxido-reduction of heme a (and Cu(A)) and Cu(B) exhibit membrane vectoriality, i.e. protons are taken up from the inner space upon reduction of these centers and released in the outer space upon their oxidation. Redox Bohr protons coupled to anaerobic oxido-reduction of heme a(3) do not, on the contrary, exhibit vectorial nature: protons are exchanged only with the outer space. A model of the proton pump of the oxidase, in which redox Bohr protons linked to the low-spin heme a play a central role, is described. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Interheme electron tunneling in cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme of the respiratory chain that catalyzes respiratory reduction of dioxygen (O(2)) to water in all eukaryotes and many aerobic bacteria. CcO, and its homologs among the heme-copper oxidases, has an active site composed of an oxygen-binding heme and a copper center in the vicinity, plus another heme group that donates electrons to this site. In most oxidoreduction enzymes, electron transfer (eT) takes place by quantum-mechanical electron tunneling. Here we show by independent molecular dynamics and quantum-chemical methods that the heme-heme eT in CcO differs from the majority of cases in having an exceptionally low reorganization energy. We show that the rate of interheme eT in CcO may nevertheless be predicted by the Moser-Dutton equation if reinterpreted as the average of the eT rates between all individual atoms of the donor and acceptor weighed by the respective packing densities between them. We argue that this modification may be necessary at short donor/acceptor distances comparable to the donor/acceptor radii.

Mechanisms underlying dioxygen reduction in laccases. Structural and modelling studies focusing on proton transfer

BACKGROUND

Laccases are enzymes that couple the oxidation of substrates with the reduction of dioxygen to water. They are the simplest members of the multi-copper oxidases and contain at least two types of copper centres; a mononuclear T1 and a trinuclear that includes two T3 and one T2 copper ions. Substrate oxidation takes place at the mononuclear centre whereas reduction of oxygen to water occurs at the trinuclear centre.

RESULTS

In this study, the CotA laccase from Bacillus subtilis was used as a model to understand the mechanisms taking place at the molecular level, with a focus in the trinuclear centre. The structures of the holo-protein and of the oxidised form of the apo-protein, which has previously been reconstituted in vitro with Cu(I), have been determined. The former has a dioxygen moiety between the T3 coppers, while the latter has a monoatomic oxygen, here interpreted as a hydroxyl ion. The UV/visible spectra of these two forms have been analysed in the crystals and compared with the data obtained in solution. Theoretical calculations on these and other structures of CotA were used to identify groups that may be responsible for channelling the protons that are needed for reduction of dioxygen to water.

CONCLUSIONS

These results present evidence that Glu 498 is the only proton-active group in the vicinity of the trinuclear centre. This strongly suggests that this residue may be responsible for channelling the protons needed for the reduction. These results are compared with other data available for these enzymes, highlighting similarities and differences within laccases and multicopper oxidases.

Functional interactions between membrane-bound transporters and membranes

One key role of many cellular membranes is to hold a transmembrane electrochemical ion gradient that stores free energy, which is used, for example, to generate ATP or to drive transmembrane transport processes. In mitochondria and many bacteria, the gradient is maintained by proton-transport proteins that are part of the respiratory (electron-transport) chain. Even though our understanding of the structure and function of these proteins has increased significantly, very little is known about the specific role of functional protein-membrane and membrane-mediated protein-protein interactions. Here, we have investigated the effect of membrane incorporation on proton-transfer reactions within the membrane-bound proton pump cytochrome c oxidase. The results show that the membrane acts to accelerate proton transfer into the enzyme's catalytic site and indicate that the intramolecular proton pathway is wired via specific amino acid residues to the two-dimensional space defined by the membrane surface. We conclude that the membrane not only acts as a passive barrier insulating the interior of the cell from the exterior solution, but also as a component of the energy-conversion machinery.

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.

MCCE2: improving protein pKa calculations with extensive side chain rotamer sampling

Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pK(a)s. Explicit changes in side chain conformations throughout a titration create a position dependent, heterogeneous dielectric response giving a more accurate picture of coupled ionization and position changes. The MCCE2 methods for choosing a group of input heavy atom and proton positions are described. The pK(a)s calculated with different isosteric conformers, heavy atom rotamers and proton positions, with different degrees of optimization are tested against a curated group of 305 experimental pK(a)s in 33 proteins. QUICK calculations, with rotation around Asn and Gln termini, sampling His tautomers and torsion minimum hydroxyls yield an RMSD of 1.34 with 84% of the errors being <1.5 pH units. FULL calculations adding heavy atom rotamers and side chain optimization yield an RMSD of 0.90 with 90% of the errors <1.5 pH unit. Good results are also found for pK(a)s in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pK(a)s by creating more neutral than ionized conformers. Methods for correcting these errors are introduced. Calculations are compared with multiple X-ray and NMR derived structures in 36 soluble proteins. Calculations with X-ray structures give significantly better pK(a)s. Results with the default protein dielectric constant of 4 are as good as those using a value of 8. The MCCE2 program can be downloaded from http://www.sci.ccny.cuny.edu/~mcce.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: new insights into the active site and the proton transfer pathways

The structure of the two-subunit cytochrome c oxidase from Paracoccus denitrificans has been refined using X-ray cryodata to 2.25 A resolution in order to gain further insights into its mechanism of action. The refined structural model shows a number of new features including many additional solvent and detergent molecules. The electron density bridging the heme a(3) iron and Cu(B) of the active site is fitted best by a peroxo-group or a chloride ion. Two waters or OH(-) groups do not fit, one water (or OH(-)) does not provide sufficient electron density. The analysis of crystals of cytochrome c oxidase isolated in the presence of bromide instead of chloride appears to exclude chloride as the bridging ligand. In the D-pathway a hydrogen bonded chain of six water molecules connects Asn131 and Glu278, but the access for protons to this water chain is blocked by Asn113, Asn131 and Asn199. The K-pathway contains two firmly bound water molecules, an additional water chain seems to form its entrance. Above the hemes a cluster of 13 water molecules is observed which potentially form multiple exit pathways for pumped protons. The hydrogen bond pattern excludes that the Cu(B) ligand His326 is present in the imidazolate form.

Intramolecular electron transfer processes in Cu(B)-deficient cytochrome bo studied by pulse radiolysis

The Escherichia coli cytochrome bo is a heme-copper terminal ubiquinol oxidase, and functions as a redox-driven proton pump. We applied pulse radiolysis technique for studying the one-electron reduction processes in the Cu(B)-deficient mutant, His333Ala. We found that the Cu(B) deficiency suppressed the heme b-to-heme o electron transfer two order of the magnitude (4.0 x 10(2) s(-1)), as found for the wild-type enzyme in the presence of 1 mM KCN (3.0 x 10(2) s(-1)). Potentiometric analysis of the His333Ala mutant revealed the 40 mV decrease in the E(m) value for low-spin heme b and the 160 mV increase in the E(m) value of high-spin heme o. Our results indicate that Cu(B) not only serves as one-electron donor to the bound dioxygen upon the O-O bond cleavage, but also facilitates dioxygen reduction at the heme-copper binuclear centre by modulating the E(m) value of heme o through magnetic interactions. And the absence of a putative OH(-) bound to Cu(B) seems not to affect the uptake of the first chemical proton via K-channel in the His333Ala mutant.

Charge parameterization of the metal centers in cytochrome c oxidase

Reliable atomic point charges are of key importance for a correct description of the electrostatic interactions when performing classical, force field based simulations. Here, we present a systematic procedure for point charge derivation, based on quantum mechanical methodology suited for the systems at hand. A notable difference to previous procedures is to include an outer region around the actual system of interest. At the cost of increasing the system sizes, here up to 265 atoms, including the surroundings achieves near-neutrality for the systems as well as structural stability, important factors for reliable charge distributions. In addition, the common problem of converting between C--H bonds and C--C bonds at the border vanishes. We apply the procedure to the four redox-active metal centers of cytochrome c oxidase: Cu(A), haem a, haem a(3), and Cu(B). Several relevant charge and ligand states are considered. Charges for two different force fields, CHARMM and AMBER, are presented.

Looking for the minimum common denominator in haem-copper oxygen reductases: towards a unified catalytic mechanism

Haem-copper oxygen reductases are transmembrane protein complexes that reduce dioxygen to water and pump protons across the mitochondrial or periplasmatic membrane, contributing to the transmembrane difference of electrochemical potential. Seven years ago we proposed a classification of these enzymes into three different families (A, B and C), based on the amino acid residues of their proton channels and amino acid sequence comparison, later supported by the so far identified characteristics of the catalytic centre of members from each family. The three families have in common the same general structural fold of the catalytic subunit, which contains the same or analogous prosthetic groups, and proton channels. These observations raise the hypothesis that the mechanisms for dioxygen reduction, proton pumping and the coupling of the two processes may be the same for all these enzymes. Under this hypothesis, they should be performed and controlled by the same or equivalent elements/events, and the identification of retained elements in all families will reveal their importance and may prompt the definition of the enzyme operating mode. Thus, we believe that the search for a minimum common denominator has a crucial importance, and in this article we highlight what is already established for the haem-copper oxygen reductases and emphasize the main questions still unanswered in a comprehensive basis.

Carboxyl group functions in the heme-copper oxidases: information from mid-IR vibrational spectroscopy

Carboxyl groups of possible functional importance in bovine and bacterial cytochrome c oxidases (CcO) are reviewed and assessed. A critical analysis is presented of available mid-infrared vibrational data that pertain to these functional carboxyl groups. These data and their interpretations are discussed in relation to current models of the mechanism of proton and electron coupling in the protonmotive CcO superfamily.

Glutamic acid 242 is a valve in the proton pump of cytochrome c oxidase

Aerobic life is based on a molecular machinery that utilizes oxygen as a terminal electron sink. The membrane-bound cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water in mitochondria and many bacteria. The energy released in this reaction is conserved by pumping protons across the mitochondrial or bacterial membrane, creating an electrochemical proton gradient that drives production of ATP. A crucial question is how the protons pumped by CcO are prevented from flowing backwards during the process. Here, we show by molecular dynamics simulations that the conserved glutamic acid 242 near the active site of CcO undergoes a protonation state-dependent conformational change, which provides a valve in the pumping mechanism. The valve ensures that at any point in time, the proton pathway across the membrane is effectively discontinuous, thereby preventing thermodynamically favorable proton back-leakage while maintaining an overall high efficiency of proton translocation. Suppression of proton leakage is particularly important in mitochondria under physiological conditions, where production of ATP takes place in the presence of a high electrochemical proton gradient.

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

Redox properties of Thermus thermophilus ba3: different electron-proton coupling in oxygen reductases?

A comprehensive study of the thermodynamic redox behavior of the hemes of the ba3 enzyme from Thermus thermophilus, a B-type heme-copper oxygen reductase, is presented. This enzyme, in contrast to those having a single type of heme, allows the B- and A-type hemes to be monitored separately by visible spectroscopy and the reduction potential of each heme to be determined unequivocally. The relative order of the midpoint reduction potentials of each center changed in the pH range from 6 to 8.4, and both hemes present a significant redox-Bohr effect. For instance, at pH 7, the midpoint reduction potentials of the hemes B and A3 are 213 mV and 285 mV, respectively, whereas at pH 8.4, the order is reversed: 246 mV for heme B and 199 mV for heme A3. The existence of redox anticooperativity was established by introducing a redox interaction parameter in a model of pairwise interacting redox centers.

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.

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 a3, including Asp399 and the CuB 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.

Storage of an excess proton in the hydrogen-bonded network of the d-pathway of cytochrome C oxidase: identification of a protonated water cluster

The mechanism of proton transport in the D-pathway of cytochrome c oxidase (CcO) is further elucidated through examining a protonated water/hydroxyl cluster inside the channel. The second generation multi-state empirical valence bond (MS-EVB2) model was employed in a molecular dynamics study based on a high-resolution X-ray structure to simulate the interaction of the excess proton with the channel environment. Our results indicate that a hydrogen-bonded network consisting of about 5 water molecules surrounded by three side chains and two backbone groups (S197, S200, S201, F108) is involved in storage and translocation of an excess proton to the extracellular side of CcO.

Energy transduction: proton transfer through the respiratory complexes

A series of metalloprotein complexes embedded in a mitochondrial or bacterial membrane utilize electron transfer reactions to pump protons across the membrane and create an electrochemical potential (DeltamuH+). Current understanding of the principles of electron-driven proton transfer is discussed, mainly with respect to the wealth of knowledge available from studies of cytochrome c oxidase. Structural, experimental, and theoretical evidence supports the model of long-distance proton transfer via hydrogen-bonded water chains in proteins as well as the basic concept that proton uptake and release in a redox-driven pump are driven by charge changes at the membrane-embedded centers. Key elements in the pumping mechanism may include bound water, carboxylates, and the heme propionates, arginines, and associated water above the hemes. There is evidence for an important role of subunit III and proton backflow, but the number and nature of gating mechanisms remain elusive, as does the mechanism of physiological control of efficiency.

Calculated proton uptake on anaerobic reduction of cytochrome C oxidase: is the reaction electroneutral?

Cytochrome c oxidase is a transmembrane proton pump that builds an electrochemical gradient using chemical energy from the reduction of O(2). Ionization states of all residues were calculated with Multi-Conformation Continuum Electrostatics (MCCE) in seven anaerobic oxidase redox states ranging from fully oxidized to fully reduced. One long-standing problem is how proton uptake is coupled to the reduction of the active site binuclear center (BNC). The BNC has two cofactors: heme a(3) and Cu(B). If the protein needs to maintain electroneutrality, then 2 protons will be bound when the BNC is reduced by 2 electrons in the reductive half of the reaction cycle. The effective pK(a)s of ionizable residues around the BNC are evaluated in Rhodobacter sphaeroides cytochrome c oxidase. At pH 7, only a hydroxide coordinated to Cu(B) shifts its pK(a) from below 7 to above 7 and so picks up a proton when heme a(3) and Cu(B) are reduced. Glu I-286, Tyr I-288, His I-334, and a second hydroxide on heme a(3) all have pK(a)s above 7 in all redox states, although they have only 1.6-3.5 DeltapK units energy cost for deprotonation. Thus, at equilibrium, they are protonated and cannot serve as proton acceptors. The propionic acids near the BNC are deprotonated with pK(a)s well below 7. They are well stabilized in their anionic state and do not bind a proton upon BNC reduction. This suggests that electroneutrality in the BNC is not maintained during the anaerobic reduction. Proton uptake on reduction of Cu(A), heme a, heme a(3), and Cu(B) shows approximately 2.5 protons bound per 4 electrons, in agreement with prior experiments. One proton is bound by a hydroxyl group in the BNC and the rest to groups far from the BNC. The electrochemical midpoint potential (E(m)) of heme a is calculated in the fully oxidized protein and with 1 or 2 electrons in the BNC. The E(m) of heme a shifts down when the BNC is reduced, which agrees with prior experiments. If the BNC reduction is electroneutral, then the heme a E(m) is independent of the BNC redox state.

Electrostatic environment of hemes in proteins: pK(a)s of hydroxyl ligands

The pK(a)s of ferric aquo-heme and aquo-heme electrochemical midpoints (E(m)s) at pH 7 in sperm whale myoglobin, Aplysia myoblogin, hemoglobin I, heme oxygenase 1, horseradish peroxidase and cytochrome c oxidase were calculated with Multi-Conformation Continuum Electrostatics (MCCE). The pK(a)s span 3.3 pH units from 7.6 in heme oxygenase 1 to 10.9 in peroxidase, and the E(m)s range from -250 mV in peroxidase to 125 mV in Aplysia myoglobin. Proteins with higher in situ ferric aquo-heme pK(a)s tend to have lower E(m)s. Both changes arise from the protein stabilizing a positively charged heme. However, compared with values in solution, the protein shifts the aquo-heme E(m)s more than the pK(a)s. Thus, the protein has a larger effective dielectric constant for the protonation reaction, showing that electron and proton transfers are coupled to different conformational changes that are captured in the MCCE analysis. The calculations reveal a breakdown in the classical continuum electrostatic analysis of pairwise interactions. Comparisons with DFT calculations show that Coulomb's law overestimates the large unfavorable interactions between the ferric water-heme and positively charged groups facing the heme plane by as much as 60%. If interactions with Cu(B) in cytochrome c oxidase and Arg 38 in horseradish peroxidase are not corrected, the pK(a) calculations are in error by as much as 6 pH units. With DFT corrected interactions calculated pK(a)s and E(m)s differ from measured values by less than 1 pH unit or 35 mV, respectively. The in situ aquo-heme pK(a) is important for the function of cytochrome c oxidase since it helps to control the stoichiometry of proton uptake coupled to electron transfer [Song, Michonova-Alexova, and Gunner (2006) Biochemistry 45, 7959-7975].