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.     

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.     

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.     

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.     

Electron transfer pathways in cytochrome c oxidase

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

Oxygen and proton pathways in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump, which couples the reduction of oxygen to water to the translocation of protons across the membrane. The recently solved x-ray structures of cytochrome c oxidase permit molecular dynamics simulations of the underlying transport processes. To eventually establish the proton pump mechanism, we investigate the transport of the substrates, oxygen and protons, through the enzyme. Molecular dynamics simulations of oxygen diffusion through the protein reveal a well-defined pathway to the oxygen-binding site starting at a hydrophobic cavity near the membrane-exposed surface of subunit I, close to the interface to subunit III. A large number of water sites are predicted within the protein, which could play an essential role for the transfer of protons in cytochrome c oxidase. The water molecules form two channels along which protons can enter from the cytoplasmic (matrix) side of the protein and reach the binuclear center. A possible pumping mechanism is proposed that involves a shuttling motion of a glutamic acid side chain, which could then transfer a proton to a propionate group of heme α₃. Proteins 30:100–107, 1998. © 1998 Wiley-Liss, Inc.     

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 a3/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 H2O in soluble bovine heart cytochrome c oxidase; (ii) proton release/uptake associated with anaerobic oxidation/reduction of heme a/Cu(A) and heme a3/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 a3/Cu(B) with acid/base clusters, C1 and C2 respectively, and protonmotive steps of the reduction of O2 to water are involved in proton pumping.     

Molecular structure and function of bacterial nitric oxide reductase

The crystal structure of the membrane-integrated nitric oxide reductase cNOR from Pseudomonas aeruginosa was determined. The smaller NorC subunit of cNOR is comprised of 1 trans-membrane helix and a hydrophilic domain, where the heme c is located, while the larger NorB subunit consists of 12 trans-membrane helices, which contain heme b and the catalytically active binuclear center (heme b(3) and non-heme Fe(B)). The roles of the 5 well-conserved glutamates in NOR are discussed, based on the recently solved structure. Glu211 and Glu280 appear to play an important role in the catalytic reduction of NO at the binuclear center by functioning as a terminal proton donor, while Glu215 probably contributes to the electro-negative environment of the catalytic center. Glu135, a ligand for Ca(2+) sandwiched between two heme propionates from heme b and b(3), and the nearby Glu138 appears to function as a structural factor in maintaining a protein conformation that is suitable for electron-coupled proton transfer from the periplasmic region to the active site. On the basis of these observations, the possible molecular mechanism for the reduction of NO by cNOR is discussed. This article is part of a Special Issue entitled: Respiratory Oxidases.     

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N(2)O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666-70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb(3) oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.     

Structural Assignment of Spectra by Characterization of Conformational Substates in Bound MbCO

Residue motions of the distal heme pocket and bound CO ligand of carbonmonoxy Myoglobin are studied using a combination of molecular dynamics simulations and quantum chemical methods. Using mixed quantum mechanics/molecular mechanics calculations together with sampling from molecular dynamics simulations (QM/MM(MD)), the experimentally observed spectroscopic A₀ and A₁ substates of the bound CO ligand are assigned to the open and closed conformation of His⁶⁴ and the His_ɛ⁶⁴ tautomer, respectively. Several previously proposed origins of the A₃ substate, including rotamers of the doubly protonated His⁶⁴H⁺ side chain, His⁶⁴H⁺ inside the distal pocket, and cooperative motions with Arg⁴⁵, are investigated with QM/MM(MD). However, the signatures of the calculated infrared spectra do not agree with the experimentally observed ones. For additional insight on this, extensive molecular dynamics simulations are used together with improved electrostatics for the bound ligand. A CO fluctuating charge model is developed to describe the ab initio dipole and quadrupole moments of the bound ligand. CO absorption spectra are then obtained directly from the dynamics simulations. Finally, the electrostatics of the heme pocket is examined in detail in an attempt to determine the structural origins of the observed spectroscopic A-states from MD simulations. However, contrary to related simulations for unbound CO in myoglobin, the shifts and splittings for carbonmonoxy Myoglobin are generally small and difficult to relate to structural change. This suggests that coupling of the CO motion to other degrees of freedom, such as the Fe-CO stretching and bending, is important to correctly describe the dynamics of bound CO in myoglobin.     

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.     

Thermodynamic redox behavior of the heme centers of cbb3 heme-copper oxygen reductase from Bradyrhizobium japonicum

A comprehensive study of the thermodynamic redox behavior of the hemes from the cbb3 oxygen reductase from Bradyrhizobium japonicum was performed. This enzyme is a member of the C-type heme-copper oxygen reductase superfamily and has three subunits with six redox centers: four low-spin hemes and a high-spin heme and one copper ion, composing the site where oxygen is reduced. In this analysis, the visible spectra and redox properties of the five heme centers were deconvoluted. Their redox profiles and the pH dependence of the midpoint reduction potentials (redox-Bohr effect) were investigated. The reference reduction potentials (defined for a state where all centers are reduced) and homotropic interaction potentials were determined in the framework of a model of pairwise interacting redox centers. At pH 7.7, the reference reduction potentials for the three hemes c are 390, 300, and 220 mV, with low interaction potentials between them, weaker than -15 mV. For hemes b and b3, reference reduction potentials of 375 and 290 mV, respectively, were obtained; these two redox centers show an interaction potential weaker than -60 mV. The midpoint reduction potentials of all five hemes are pH-dependent. The study of these thermodynamic parameters is important in understanding the coupling mechanism of the redox and chemical processes during oxygen reduction. The analysis of the thermodynamic redox behavior of the cbb3 oxygen reductase contributes to the investigation of the mechanism of electron transfer and proton translocation by heme-copper oxygen reductases in general and indicates a thermodynamic coupling for the electron and proton transfer mechanisms.     

Water dynamics simulation as a tool for probing proton transfer pathways in a heptahelical membrane protein

The proton transfer pathway in a heptahelical membrane protein, the light-driven proton pump bacteriorhodopsin (BR), is probed by a combined approach of structural analysis of recent X-ray models and molecular dynamics (MD) simulations that provide the diffusion pathways of internal and external water molecules. Analyzing the hydrogen-bond contact frequencies of the water molecules with protein groups, the complete proton pathway through the protein is probed. Beside the well-known proton binding sites in the protein interior-the protonated Schiff base, Asp85 and Asp96, and the H(5)O(2) (+) complex stabilized by Glu204 and Glu194-the proton release and uptake pathways to the protein surfaces are described in great detail. Further residues were identified, by mutation of which the proposed pathways can be verified. In addition the diffusion pathway of water 502 from Lys216 to Asp96 is shown to cover the positions of the intruding waters 503 and 504 in the N-intermediate. The transiently established water chain in the N-state provides a proton pathway from Asp96 to the Schiff base in the M- to N-transition in a Grotthus-like mechanism, as concluded earlier from time-resolved Fourier transform infrared experiments [le Coutre et al., Proc Nat Acad Sci USA 1995;92:4962-4966].     

Proton pumping by cytochrome c oxidase – A 40 year anniversary

Cytochrome c oxidase is a remarkable energy transducer that seems to work almost purely by Coulombic principles without the need for significant protein conformational changes. In recent years it has become possible to follow key partial reactions of the catalytic cycle in real time, both with respect to electron and proton movements. These experiments have largely set the stage for the proton pump mechanism. The structures of the catalytic binuclear heme?copper site that is common to the huge family of heme?copper oxidases, are today well understood throughout the catalytic cycle of oxygen reduction to water based on both spectroscopic studies and quantum chemical calculations. Here, we briefly review this progress, and add some recent details into how the proton pump mechanism is protected from failure by leakage.     

Theoretical investigations on Azotobacter vinelandii ferredoxin I: effects of electron transfer on protein dynamics

Structural, energetic, and dynamical studies of Azotobacter vinelandii ferredoxin I are presented for native and mutant forms. The protein contains two iron-sulfur clusters, one of which ([3Fe-4S]) is believed to play a central role in the electron-coupled proton transfer. Different charge sets for the [3Fe-4S] cluster in its reduced and oxidized state are calculated with broken symmetry ab initio density functional theory methods and used in molecular dynamics (MD) simulations. The validity of the ab initio calculations is assessed by comparing partially optimized structures of the [3Fe-4S] clusters with x-ray structures. Possible proton transfer pathways between the protein and the iron-sulfur cluster are examined by both MD simulations and ab initio calculations. The MD simulations identify three main-chain hydrogen atoms--HN(13), HN(14), and HN(16)--that are within H-bonding distance of the [3Fe-4S] cluster throughout the MD simulations. They could thus play a role in the proton transfer from the protein to the iron-sulfur cluster. By contrast, the HD2(15) atom of the Asp-15 is seldom close enough to the [3Fe-4S] cluster to transfer a proton. Poisson-Boltzmann calculations indicate that there is a low, but nonzero probability, that Asp-15 is protonated at pH 7; this is a requirement for it to serve as a proton donor. Ab initio calculations with a fragment model for the protein find similar behavior for the transfer of a proton from the OH of the protonated side chain and the main-chain NH of Asp-15. The existence of a stable salt bridge between Asp-15 and Lys-84 in the D15E mutant, versus its absence in the wild-type, has been suggested as the cause of the difference in the rate of proton transfer. Extensive MD simulations were done to test this idea; the results do not support the proposal. The present findings, together with the available data, serve as the basis for an alternative proposal for the mechanism of the coupled electron-proton transfer reaction in ferredoxin I.     

How inter-subunit contacts in the membrane domain of complex I affect proton transfer energetics

The respiratory complex I is a redox-driven proton pump that employs the free energy released from quinone reduction to pump protons across its complete ca. 200?Å wide membrane domain. Despite recently resolved structures and molecular simulations, the exact mechanism for the proton transport process remains unclear. Here we combine large-scale molecular simulations with quantum chemical density functional theory (DFT) models to study how contacts between neighboring antiporter-like subunits in the membrane domain of complex I affect the proton transfer energetics. Our combined results suggest that opening of conserved Lys/Glu ion pairs within each antiporter-like subunit modulates the barrier for the lateral proton transfer reactions. Our work provides a mechanistic suggestion for key coupling effects in the long-range force propagation process of complex I.     

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 CuB 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 delta-nitrogen of His291 ligand of CuB 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 pKa of His291 for an extended model in which the whole Fe(a3)-CuB 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 pKa 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.     

Evolutionary migration of a post-translationally modified active-site residue in the proton-pumping heme-copper oxygen reductases

In the respiratory chains of aerobic organisms, oxygen reductase members of the heme-copper superfamily couple the reduction of O2 to proton pumping, generating an electrochemical gradient. There are three distinct families of heme-copper oxygen reductases: A, B, and C types. The A- and B-type oxygen reductases have an active-site tyrosine that forms a unique cross-linked histidine-tyrosine cofactor. In the C-type oxygen reductases (also called cbb3 oxidases), an analogous active-site tyrosine has recently been predicted by molecular modeling to be located within a different transmembrane helix in comparison to the A- and B-type oxygen reductases. In this work, Fourier-transform mass spectrometry is used to show that the predicted tyrosine forms a histidine-tyrosine cross-linked cofactor in the active site of the C-type oxygen reductases. This is the first known example of the evolutionary migration of a post-translationally modified active-site residue. It also verifies the presence of a unique cofactor in all three families of proton-pumping respiratory oxidases, demonstrating that these enzymes likely share a common reaction mechanism and that the histidine-tyrosine cofactor may be a required component for proton pumping.     

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.     

Quantum chemistry as a tool in bioenergetics

Recent developments of quantum chemical methods have made it possible to tackle crucial questions in bioenergetics. The most important systems, cytochrome c oxidase in cellular respiration and photosystem II (PSII) in photosynthesis will here be used as examples to illustrate the power of the quantum chemical tools. One main contribution from quantum chemistry is to put mechanistic suggestions onto an energy scale. Accordingly, free energy profiles can be constructed both for reduction of molecular oxygen in cytochrome c oxidase and water oxidation in PSII, including O-O bond cleavage and formation, and also proton pumping in cytochrome c oxidase. For the construction of the energy diagrams, the computational results sometimes have to be combined with experimental information, such as reduction potentials and rate constants for individual steps in the reactions.