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

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Specific inhibition of proton pumping by the T315V mutation in the K channel of cytochrome ba3 from Thermus thermophilus

Cytochrome ba₃ from Thermus thermophilus belongs to the B family of heme-copper oxidases and pumps protons across the membrane with an as yet unknown mechanism. The K channel of the A family heme-copper oxidases provides delivery of a substrate proton from the internal water phase to the binuclear heme-copper center (BNC) during the reductive phase of the catalytic cycle, while the D channel is responsible for transferring both substrate and pumped protons. By contrast, in the B family oxidases there is no D-channel and the structural equivalent of the K channel seems to be responsible for the transfer of both categories of protons. Here we have studied the effect of the T315V substitution in the K channel on the kinetics of membrane potential generation coupled to the oxidative half-reaction of the catalytic cycle of cytochrome ba₃. The results suggest that the mutated enzyme does not pump protons during the reaction of the fully reduced form with molecular oxygen in a single turnover. Specific inhibition of proton pumping in the T315V mutant appears to be a consequence of inability to provide rapid (τ ~ 100 μs) reprotonation of the internal transient proton donor(s) of the K channel. In contrast to the A family, the K channel of the B-type oxidases is necessary for the electrogenic transfer of both pumped and substrate protons during the oxidative half-reaction of the catalytic cycle.

Activation of O2 and NO in heme-copper oxidases - mechanistic insights from computational modelling

Heme-copper oxidases are transmembrane enzymes involved in aerobic and anaerobic respiration. The largest subgroup contains the cytochrome c oxidases (CcO), which reduce molecular oxygen to water. A significant part of the free energy released in this exergonic process is conserved as an electrochemical gradient across the membrane, via two processes, electrogenic chemistry and proton pumping. A deviant subgroup is the cytochrome c dependent NO reductases (cNOR), which reduce nitric oxide to nitrous oxide and water. This is also an exergonic reaction, but in this case none of the released free energy is conserved. Computational studies applying hybrid density functional theory to cluster models of the bimetallic active sites in the heme-copper oxidases are reviewed. To obtain a reliable description of the reaction mechanisms, energy profiles of the entire catalytic cycles, including the reduction steps have to be constructed. This requires a careful combination of computational results with certain experimental data. Computational studies have elucidated mechanistic details of the chemical parts of the reactions, involving cleavage and formation of covalent bonds, which have not been obtainable from pure experimental investigations. Important insights regarding the mechanisms of energy conservation have also been gained. The computational studies show that the reduction potentials of the active site cofactors in the CcOs are large enough to afford electrogenic chemistry and proton pumping, i.e. efficient energy conservation. These results solve a conflict between different types of experimental data. A mechanism for the proton pumping, involving a specific and crucial role for the active site tyrosine, conserved in all CcOs, is suggested. For the cNORs, the calculations show that the low reduction potentials of the active site cofactors are optimized for fast elimination of the toxic NO molecules. At the same time, the low reduction potentials lead to endergonic reduction steps with high barriers. To prevent even higher barriers, which would lead to a too slow reaction, when the electrochemical gradient across the membrane is present, the chemistry must occur in a non-electrogenic manner. This explains why there is no energy conservation in cNOR.

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

A Water Molecule Residing in the Fea33+···CuB2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study

Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.

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

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

DFT Fea3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase

Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.

Geometric and Electronic Structure Contributions to O-O Cleavage and the Resultant Intermediate Generated in Heme-Copper Oxidases

This study investigates the mechanism of O-O bond cleavage in heme-copper oxidase (HCO) enzymes, combining experimental and computational insights from enzyme intermediates and synthetic models. It is determined that HCOs undergo a proton-initiated O-O cleavage mechanism where a single water molecule in the active site enables proton transfer (PT) from the cross-linked tyrosine to a peroxo ligand bridging the heme Fe^III and Cu^II, and multiple H-bonding interactions lower the tyrosine p K_a. Due to sterics within the active site, the proton must either transfer initially to the O(Fe) (a high-energy intermediate), or from another residue over a ∼10 Å distance to reach the O(Cu) atom directly. While the distance between the H⁺ donor (Tyr) and acceptor (O(Cu)) results in a barrier to PT, this separation is critical for the low barrier to O-O cleavage as it enhances backbonding from Fe into the O₂^2- σ* orbital. Thus, PT from Tyr precedes O-O elongation and is rate-limiting, consistent with available kinetic data. The electron transfers from tyrosinate after the barrier via a superexchange pathway provided by the cross-link, generating intermediate P_M. P_M is evaluated using available experimental data. The geometric structure contains an Fe^IV═O that is H-bonded to the Cu^II-OH. The electronic structure is a singlet, where the Fe^IV and Cu^II are antiferromagnetically coupled through the H-bond between the oxo(Fe) and hydroxo(Cu) ligands, while the Cu^II and Tyr^• are ferromagnetically coupled due their delocalization into orthogonal magnetic orbitals on the cross-linked His residue. These findings provide critical insights into the mechanism of efficient O₂ reduction in HCOs, and the nature of the P_M intermediate that couples this reaction to proton pumping.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.

A Water Dimer Shift Activates a Proton Pumping Pathway in the PR → F Transition of ba3 Cytochrome c Oxidase

Broken-symmetry density functional calculations have been performed on the [Fe_a3⁴⁺,Cu_B²⁺] state of the dinuclear center (DNC) for the P_R → F part of the catalytic cycle of ba₃ cytochrome c oxidase (CcO) from Thermus thermophilus (Tt), using the OLYP-D3-BJ functional. The calculations show that the movement of the H₂O molecules in the DNC affects the pK_a values of the residue side chains of Tyr237 and His376⁺, which are crucial for proton transfer/pumping in ba₃ CcO from Tt. The calculated lowest energy structure of the DNC in the [Fe_a3⁴⁺,Cu_B²⁺] state (state F) is of the form Fe_a3⁴⁺═O^2-···Cu_B²⁺, in which the H₂O ligand that resulted from protonation of the OH^- ligand in the P_R state is dissociated from the Cu_B²⁺ site. The calculated Fe_a3⁴⁺═O^2- distance in F (1.68 Å) is 0.03 Å longer than that in P_R (1.65 Å), which can explain the different Fe_a3⁴⁺═O^2- stretching modes in P (804 cm^-1) and F (785 cm^-1) identified by resonance Raman experiments. In this F state, the Cu_B²⁺···O^2- (ferryl-oxygen) distance is only around 2.4 Å. Hence, the subsequent O_H state [Fe_a3³⁺-OH^--Cu_B²⁺] with a μ-hydroxo bridge can be easily formed, as shown by our calculations.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

ns-μs Time-Resolved Step-Scan FTIR of ba₃ Oxidoreductase from Thermus thermophilus: Protonic Connectivity of w941-w946-w927

Time-resolved step-scan FTIR spectroscopy has been employed to probe the dynamics of the ba₃ oxidoreductase from Thermus thermophilus in the ns-μs time range and in the pH/pD 6–9 range. The data revealed a pH/pD sensitivity of the D372 residue and of the ring-A propionate of heme a₃. Based on the observed transient changes a model in which the protonic connectivity of w941-w946-927 to the D372 and the ring-A propionate of heme a₃ is described.

Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.

Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. This article presents parameters for the cofactors of ba₃-type CcO that are compatible with the all-atom Amber ff12SB and ff14SB force fields. Specifically, parameters were developed for the Cu_A pair, heme b, and the dinuclear center that consists of heme a₃ and Cu_B bridged by a hydroperoxo group. The data includes geometries in XYZ coordinate format for cluster models that were employed to compute proton transfer energies and derive bond parameters and point charges for the force field using density functional theory. Also included are the final parameter files that can be employed with the Amber leap program to generate input files for molecular dynamics simulations with the Amber software package. Based on the high resolution (1.8 Å) X-ray crystal structure of the ba₃-type CcO from Thermus thermophilus (Protein Data Bank ID number PDB: 3S8F), we built a model that is embedded in a POPC lipid bilayer membrane and solvated with TIP3P water molecules and counterions. We provide PDB data files of the initial model and the equilibrated model that can be used for further studies.

A broken-symmetry density functional study of structures, energies, and protonation states along the catalytic O-O bond cleavage pathway in ba3 cytochrome c oxidase from Thermus thermophilus

Broken-symmetry density functional calculations have been performed on the [Fea3, CuB] dinuclear center (DNC) of ba3 cytochrome c oxidase from Thermus thermophilus in the states of [Fea3(3+)-(HO2)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], using both PW91-D3 and OLYP-D3 functionals. Tyr237 is a special tyrosine cross-linked to His233, a ligand of CuB. The calculations have shown that the DNC in these states strongly favors the protonation of His376, which is above propionate-A, but not of the carboxylate group of propionate-A. The energies of the structures obtained by constrained geometry optimizations along the O-O bond cleavage pathway between [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] have also been calculated. The transition of [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] → [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] shows a very small barrier, which is less than 3.0/2.0 kcal mol(-1) in PW91-D3/OLYP-D3 calculations. The protonation state of His376 does not affect this O-O cleavage barrier. The rate limiting step of the transition from state A (in which O2 binds to Fea3(2+)) to state PM ([Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], where the O-O bond is cleaved) in the catalytic cycle is, therefore, the proton transfer originating from Tyr237 to O-O to form the hydroperoxo [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] state. The importance of His376 in proton uptake and the function of propionate-A/neutral-Asp372 as a gate to prevent the proton from back-flowing to the DNC are also shown.

Synthetic heme/copper assemblies: toward an understanding of cytochrome c oxidase interactions with dioxygen and nitrogen oxides

Our long-time niche in synthetic biological inorganic chemistry has been to design ligands and generate coordination complexes of copper or iron ions or both, those reacting with dioxygen (O2) or nitrogen oxides (e.g., nitric oxide (NO(g)) and nitrite (NO2(-))) or both. As inspiration for this work, we turn to mitochondrial cytochrome c oxidase, which is responsible for dioxygen consumption and is also the predominant target for NO(g) and nitrite within mitochondria. In this Account, we highlight recent advances in studying synthetic heme/Cu complexes in two respects. First, there is the design, synthesis, and characterization of new O2 adducts whose further study will add insights into O2 reductive cleavage chemistry. Second, we describe how related heme/Cu constructs reduce nitrite ion to NO(g) or the reverse, oxidize NO(g) to nitrite. The reactions of nitrogen oxides occur as part of CcO's function, which is intimately tied to cellular O2 balance. We had first discovered that reduced heme/Cu compounds react with O2 giving μ-oxo heme-Fe(III)-O-Cu(II)(L) products; their properties are discussed. The O-atom is derived from dioxygen, and interrogations of these systems led to the construction and characterization of three distinctive classes of heme-peroxo complexes, two high-spin and one low-spin species. Recent investigations include a new approach to the synthesis of low-spin heme-peroxo-Cu complexes, employing a "naked" synthon, where the copper ligand denticity and geometric types can be varied. The result is a collection of such complexes; spectroscopic and structural features (by DFT calculations) are described. Some of these compounds are reactive toward reductants/protons effecting subsequent O-O cleavage. This points to how subtle improvements in ligand environment lead to a desired local structure and resulting optimized reactivity, as known to occur at enzyme active sites. The other sector of research is focused on heme/Cu assemblies mediating the redox interplay between nitrite and NO(g). In the nitrite reductase chemistry, the cupric center serves as a Lewis acid, while the heme is the redox active center providing the electron. The orientation of nitrite in approaching the ferrous heme center and N-atom binding are important. Also, detailed spectroscopic and kinetic studies of the NO(g) oxidase chemistry, in excellent agreement with theoretical calculations, reveal the intermediates and key mechanistic steps. Thus, we suggest that both chemical and biochemical heme/Cu-mediated nitrite reductase and NO(g) oxidase chemistry require N-atom binding to a ferrous heme along with cupric ion O-atom coordination, proceeding via a three-membered O-Fe-N chelate ring transition state. These important mechanistic features of heme/Cu systems interconverting NO(g) and nitrite are discussed for the first time.

Broken Symmetry DFT Calculations/Analysis for Oxidized and Reduced Dinuclear Center in Cytochrome c Oxidase: Relating Structures, Protonation States, Energies, and Mössbauer Properties in ba3 Thermus thermophilus

The Fea3(3+)···CuB(2+) dinuclear center (DNC) structure of the as-isolated oxidized ba3 cytochrome c oxidase (CcO) from Thermus thermophilus (Tt) is still not fully understood. When the proteins are initially crystallized in the oxidized state, they typically become radiolyticly reduced through X-ray irradiation. Several X-ray crystal structures of reduced ba3 CcO from Tt are available. However, depending on whether the crystals were prepared in a lipidic cubic phase environment or in detergent micelles, and whether the CcO's were chemically or radiolyticly reduced, the X-ray diffraction analysis of the crystals showed different Fea3(2+)···CuB(+) DNC structures. On the other hand, Mössbauer spectroscopic experiments on reduced and oxidized ba3 CcOs from Tt (Zimmermann et al., Proc. Natl. Acad. Sci. USA 1988, 85, 5779-5783) revealed multiple (57)Fea3(2+) and (57)Fea3(3+) components. Moreover, one of the (57)Fea3(3+) components observed at 4.2 K transformed from a proposed "low-spin" state to a different high-spin species when the temperature was increased above 190 K, whereas the other high-spin (57)Fea3(3+) component remained unchanged. In the current Article, in order to understand the heterogeneities of the DNC in both Mössbauer spectra and X-ray crystal structures, the spin crossover of one of the (57)Fea3(3+) components, and how the coordination and spin states of the Fea3(3+/2+) and Cu(2+/1+) sites relate to the heterogeneity of the DNC structures, we have applied density functional OLYP calculations to the DNC clusters established based on the different X-ray crystal structures of ba3 CcO from Tt. As a result, specific oxidized and reduced DNC structures related to the observed Mössbauer spectra and to spectral changes with temperature have been proposed. Our calculations also show that, in certain intermediate states, the His233 and His283 ligand side chains may dissociate from the CuB(+) site, and they may become potential proton loading sites during the catalytic cycle.

Proton transfer in the K-channel analog of B-type Cytochrome c oxidase from Thermus thermophilus

A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pKA value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.

An electrogenic nitric oxide reductase

Nitric oxide reductases (Nors) are members of the heme-copper oxidase superfamily that reduce nitric oxide (NO) to nitrous oxide (N₂O). In contrast to the proton-pumping cytochrome oxidases, Nors studied so far have neither been implicated in proton pumping nor have they been experimentally established as electrogenic. The copper-A-dependent Nor from Bacillus azotoformans uses cytochrome c₅₅₁ as electron donor but lacks menaquinol activity, in contrast to our earlier report (Suharti et al., 2001). Employing reduced phenazine ethosulfate (PESH) as electron donor, the main NO reduction pathway catalyzed by Cu(A)Nor reconstituted in liposomes involves transmembrane cycling of the PES radical. We show that Cu(A)Nor reconstituted in liposomes generates a proton electrochemical gradient across the membrane similar in magnitude to cytochrome aa₃, highlighting that bacilli using Cu(A)Nor can exploit NO reduction for increased cellular ATP production compared to organisms using cNor.

Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria

The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.

Mutation of a single residue in the ba3 oxidase specifically impairs protonation of the pump site

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound protein complex that couples electron transfer to O2 to proton translocation across the membrane. To elucidate the mechanism of the redox-driven proton pumping, we investigated the kinetics of electron and proton transfer in a structural variant of the ba3 oxidase where a putative "pump site" was modified by replacement of Asp372 by Ile. In this structural variant, proton pumping was uncoupled from internal electron transfer and O2 reduction. The results from our studies show that proton uptake to the pump site (time constant ∼65 μs in the wild-type cytochrome c oxidase) was impaired in the Asp372Ile variant. Furthermore, a reaction step that in the wild-type cytochrome c oxidase is linked to simultaneous proton uptake and release with a time constant of ∼1.2 ms was slowed to ∼8.4 ms, and in Asp372Ile was only associated with proton uptake to the catalytic site. These data identify reaction steps that are associated with protonation and deprotonation of the pump site, and point to the area around Asp372 as the location of this site in the ba3 cytochrome c oxidase.

A "naked" Fe(III)-(O₂²⁻)-Cu(II) species allows for structural and spectroscopic tuning of low-spin heme-peroxo-Cu complexes

Here we describe a new approach for the generation of heme-peroxo-Cu compounds, using a “naked” complex synthon, [(F₈)Fe^III-(O₂^2–)-Cu^II(MeTHF)₃]⁺ (MeTHF = 2-methyltetrahydrofuran; F₈ = tetrakis(2,6-difluorophenyl)porphyrinate). Addition of varying ligands (L) for Cu allows the generation and spectroscopic characterization of a family of high- and low-spin Fe^III-(O₂^2–)-Cu^II(L) complexes. These possess markedly varying Cu^II coordination geometries, leading to tunable Fe-O, O-O, and Cu-O bond strengths. DFT calculations accompanied by vibrational data correlations give detailed structural insights.

Detection of functional hydrogen-bonded water molecules with protonated/deprotonated key carboxyl side chains in the respiratory enzyme ba3-oxidoreductase

In this work we report FTIR data on the detection of functional hydrogen-bonded water molecules with protonated/deprotonated key carboxyl side chains in the respiratory enzyme ba₃-oxidoreductase.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.

Linking chemical electron-proton transfer to proton pumping in cytochrome c oxidase: broken-symmetry DFT exploration of intermediates along the catalytic reaction pathway of the iron-copper dinuclear complex

After a summary of the problem of coupling electron and proton transfer to proton pumping in cytochrome c oxidase, we present the results of our earlier and recent density functional theory calculations for the dinuclear Fe-a3-CuB reaction center in this enzyme. A specific catalytic reaction wheel diagram is constructed from the calculations, based on the structures and relative energies of the intermediate states of the reaction cycle. A larger family of tautomers/protonation states is generated compared to our earlier work, and a new lowest-energy pathway is proposed. The entire reaction cycle is calculated for the new smaller model (about 185-190 atoms), and two selected arcs of the wheel are chosen for calculations using a larger model (about 205 atoms). We compare the structural and redox energetics and protonation calculations with available experimental data. The reaction cycle map that we have built is positioned for further improvement and testing against experiment.

Density functional study for the bridged dinuclear center based on a high-resolution X-ray crystal structure of ba3 cytochrome c oxidase from Thermus thermophilus

Strong electron density for a peroxide type dioxygen species bridging the Fea3 and CuB dinuclear center (DNC) was observed in the high-resolution (1.8 Å) X-ray crystal structures (PDB entries 3S8G and 3S8F) of ba3 cytochrome c oxidase (CcO) from Thermus thermophilus. The crystals represent the as-isolated X-ray photoreduced CcO structures. The bridging peroxide was proposed to arise from the recombination of two radiation-produced HO(•) radicals formed either very near to or even in the space between the two metals of the DNC. It is unclear whether this peroxide species is in the O2(2-), O2(•)(-), HO2(-), or the H2O2 form and what is the detailed electronic structure and binding geometry including the DNC. In order to answer what form of this dioxygen species was observed in the DNC of the 1.8 Å X-ray CcO crystal structure (3S8G), we have applied broken-symmetry density functional theory (BS-DFT) geometric and energetic calculations (using OLYP potential) on large DNC cluster models with different Fea3-CuB oxidation and spin states and with O2(2-), O2(•)(-), HO2(-), or H2O2 in the bridging position. By comparing the DFT optimized geometries with the X-ray crystal structure (3S8G), we propose that the bridging peroxide is HO2(-). The X-ray crystal structure is likely to represent the superposition of the Fea3(2+)-(HO2(-))-CuB(+) DNC's in different states (Fe(2+) in low spin (LS), intermediate spin (IS), or high spin (HS)) with the majority species having the proton of the HO2(-) residing on the oxygen atom (O1) which is closer to the Fea3(2+) site in the Fea3(2+)-(HO-O)(-)-CuB(+) conformation. Our calculations show that the side chain of Tyr237 is likely trapped in the deprotonated Tyr237(-) anion form in the 3S8G X-ray crystal structure.

Single mutations that redirect internal proton transfer in the ba3 oxidase from Thermus thermophilus

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound proton pump. Results from earlier studies have shown that with the aa3-type oxidases proton uptake to the catalytic site and "pump site" occurs simultaneously. However, with ba3 oxidase the pump site is loaded before proton transfer to the catalytic site because the proton transfer to the latter is slower than that with the aa3 oxidases. In addition, the timing of formation and decay of catalytic intermediates is different in the two types of oxidases. In the present study, we have investigated two mutant ba3 CytcOs in which residues of the proton pathway leading to the catalytic site as well as the pump site were exchanged, Thr312Val and Tyr244Phe. Even though ba3 CytcO uses only a single proton pathway for transfer of the substrate and "pumped" protons, the amino-acid residue substitutions had distinctly different effects on the kinetics of proton transfer to the catalytic site and the pump site. The results indicate that the rates of these reactions can be modified independently by replacement of single residues within the proton pathway. Furthermore, the data suggest that the Thr312Val and Tyr244Phe mutations interfere with a structural rearrangement in the proton pathway that is rate limiting for proton transfer to the catalytic site.

Redundancy and modularity in membrane-associated dissimilatory nitrate reduction in Bacillus

The genomes of two phenotypically denitrifying type strains of the genus Bacillus were sequenced and the pathways for dissimilatory nitrate reduction were reconstructed. Results suggest that denitrification proceeds in the periplasmic space and in an analogous fashion as in Gram-negative organisms, yet with the participation of proteins that tend to be membrane-bound or membrane-associated. A considerable degree of functional redundancy was observed with marked differences between B. azotoformans LMG 9581(T) and B. bataviensis LMG 21833(T). In addition to the already characterized menaquinol/cyt c-dependent nitric oxide reductase (Suharti et al., 2001, 2004) of which the encoding genes could be identified now, evidence for another novel nitric oxide reductase (NOR) was found. Also, our analyses confirm earlier findings on branched electron transfer with both menaquinol and cytochrome c as reductants. Quite unexpectedly, both bacilli have the disposal of two parallel pathways for nitrite reduction enabling a life style as a denitrifier and as an ammonifying bacterium.

Spectroscopic and kinetic investigation of the fully reduced and mixed valence states of ba3-cytochrome c oxidase from Thermus thermophilus: a Fourier transform infrared (FTIR) and time-resolved step-scan FTIR study

The complete understanding of a molecular mechanism of action requires the thermodynamic and kinetic characterization of different states and intermediates. Cytochrome c oxidase reduces O(2) to H(2)O, a reaction coupled to proton translocation across the membrane. Therefore, it is necessary to undertake a thorough characterization of the reduced form of the enzyme and the determination of the electron transfer processes and pathways between the redox-active centers. In this study Fourier transform infrared (FTIR) and time-resolved step-scan FTIR spectroscopy have been applied to study the fully reduced and mixed valence states of cytochrome ba(3) from Thermus thermophilus. We used as probe carbon monoxide (CO) to characterize both thermodynamically and kinetically the cytochrome ba(3)-CO complex in the 5.25-10.10 pH/pD range and to study the reverse intramolecular electron transfer initiated by the photolysis of CO in the two-electron reduced form. The time-resolved step-scan FTIR data revealed no pH/pD dependence in both the decay of the transient Cu(B)(1+)-CO complex and rebinding to heme a(3) rates, suggesting that no structural change takes place in the vicinity of the binuclear center. Surprisingly, photodissociation of CO from the mixed valence form of the enzyme does not lead to reverse electron transfer from the reduced heme a(3) to the oxidized low-spin heme b, as observed in all the other aa(3) and bo(3) oxidases previously examined. The heme b-heme a(3) electron transfer is guaranteed, and therefore, there is no need for structural rearrangements and complex synchronized cooperativities. Comparison among the available structures of ba(3)- and aa(3)-cytochrome c oxidases identifies possible active pathways involved in the electron transfer processes and key structural elements that contribute to the different behavior observed in cytochrome ba(3).

Mobility of Xe atoms within the oxygen diffusion channel of cytochrome ba(3) oxidase

We use a form of "freeze-trap, kinetic crystallography" to explore the migration of Xe atoms away from the dinuclear heme a(3)/Cu(B) center in Thermus thermophilus cytochrome ba(3) oxidase. This enzyme is a member of the heme-copper oxidase superfamily and is thus crucial for dioxygen-dependent life. The mechanisms involved in the migration of oxygen, water, electrons, and protons into and/or out of the specialized channels of the heme-copper oxidases are generally not well understood. Pressurization of crystals with Xe gas previously revealed a O(2) diffusion channel in cytochrome ba(3) oxidase that is continuous, Y-shaped, 18-20 Å in length and comprised of hydrophobic residues, connecting the protein surface within the bilayer to the a(3)-Cu(B) center in the active site. To understand movement of gas molecules within the O(2) channel, we performed crystallographic analysis of 19 Xe laden crystals freeze-trapped in liquid nitrogen at selected times between 0 and 480 s while undergoing outgassing at room temperature. Variation in Xe crystallographic occupancy at five discrete sites as a function of time leads to a kinetic model revealing relative degrees of mobility of Xe atoms within the channel. Xe egress occurs primarily through the channel formed by the Xe1 → Xe5 → Xe3 → Xe4 sites, suggesting that ingress of O(2) is likely to occur by the reverse of this process. The channel itself appears not to undergo significant structural changes during Xe migration, thereby indicating a passive role in this important physiological function.

Bioenergetics at extreme temperature: Thermus thermophilus ba(3)- and caa(3)-type cytochrome c oxidases

Seven years into the completion of the genome sequencing projects of the thermophilic bacterium Thermus thermophilus strains HB8 and HB27, many questions remain on its bioenergetic mechanisms. A key fact that is occasionally overlooked is that oxygen has a very limited solubility in water at high temperatures. The HB8 strain is a facultative anaerobe whereas its relative HB27 is strictly aerobic. This has been attributed to the absence of nitrate respiration genes from the HB27 genome that are carried on a mobilizable but highly-unstable plasmid. In T. thermophilus, the nitrate respiration complements the primary aerobic respiration. It is widely known that many organisms encode multiple biochemically-redundant components of the respiratory complexes. In this minireview, the presence of the two cytochrome c oxidases (CcO) in T. thermophilus, the ba(3)- and caa(3)-types, is outlined along with functional considerations. We argue for the distinct evolutionary histories of these two CcO including their respective genetic and molecular organizations, with the caa(3)-oxidase subunits having been initially 'fused'. Coupled with sequence analysis, the ba(3)-oxidase crystal structure has provided evolutionary and functional information; for example, its subunit I is more closely related to archaeal sequences than bacterial and the substrate-enzyme interaction is hydrophobic as the elevated growth temperature weakens the electrostatic interactions common in mesophiles. Discussion on the role of cofactors in intra- and intermolecular electron transfer and proton pumping mechanism is also included.

Structural changes that occur upon photolysis of the Fe(II)(a3)-CO complex in the cytochrome ba(3)-oxidase of Thermus thermophilus: a combined X-ray crystallographic and infrared spectral study demonstrates CO binding to Cu(B)

The purpose of the work was to provide a crystallographic demonstration of the venerable idea that CO photolyzed from ferrous heme-a(3) moves to the nearby cuprous ion in the cytochrome c oxidases. Crystal structures of CO-bound cytochrome ba(3)-oxidase from Thermus thermophilus, determined at ~2.8-3.2Å resolution, reveal a Fe-C distance of ~2.0Å, a Cu-O distance of 2.4Å and a Fe-C-O angle of ~126°. Upon photodissociation at 100K, X-ray structures indicate loss of Fe(a3)-CO and appearance of Cu(B)-CO having a Cu-C distance of ~1.9Å and an O-Fe distance of ~2.3Å. Absolute FTIR spectra recorded from single crystals of reduced ba(3)-CO that had not been exposed to X-ray radiation, showed several peaks around 1975cm(-1); after photolysis at 100K, the absolute FTIR spectra also showed a significant peak at 2050cm(-1). Analysis of the 'light' minus 'dark' difference spectra showed four very sharp CO stretching bands at 1970cm(-1), 1977cm(-1), 1981cm(-1), and 1985cm(-1), previously assigned to the Fe(a3)-CO complex, and a significantly broader CO stretching band centered at ~2050cm(-1), previously assigned to the CO stretching frequency of Cu(B) bound CO. As expected for light propagating along the tetragonal axis of the P4(3)2(1)2 space group, the single crystal spectra exhibit negligible dichroism. Absolute FTIR spectrometry of a CO-laden ba(3) crystal, exposed to an amount of X-ray radiation required to obtain structural data sets before FTIR characterization, showed a significant signal due to photogenerated CO(2) at 2337cm(-1) and one from traces of CO at 2133cm(-1); while bands associated with CO bound to either Fe(a3) or to Cu(B) in "light" minus "dark" FTIR difference spectra shifted and broadened in response to X-ray exposure. In spite of considerable radiation damage to the crystals, both X-ray analysis at 2.8 and 3.2Å and FTIR spectra support the long-held position that photolysis of Fe(a3)-CO in cytochrome c oxidases leads to significant trapping of the CO on the Cu(B) atom; Fe(a3) and Cu(B) ligation, at the resolutions reported here, are otherwise unaltered.

Exploring the proton pump and exit pathway for pumped protons in cytochrome ba3 from Thermus thermophilus

The heme-copper oxygen reductases are redox-driven proton pumps. In the current work, the effects of mutations in a proposed exit pathway for pumped protons are examined in the ba(3)-type oxygen reductase from Thermus thermophilus, leading from the propionates of heme a(3) to the interface between subunits I and II. Recent studies have proposed important roles for His376 and Asp372, both of which are hydrogen-bonded to propionate-A of heme a(3), and for Glu126(II) (subunit II), which is hydrogen-bonded to His376. Based on the current results, His376, Glu126(II), and Asp372 are not essential for either oxidase activity or proton pumping. In addition, Tyr133, which is hydrogen-bonded to propionate-D of heme a(3), was also shown not to be essential for function. However, two mutations of the residues hydrogen-bonded to propionate-A, Asp372Ile and His376Asn, retain high electron transfer activity and normal spectral features but, in different preparations, either do not pump protons or exhibit substantially diminished proton pumping. It is concluded that either propionate-A of heme a(3) or possibly the cluster of groups centered about the conserved water molecule that hydrogen-bonds to both propionates-A and -D of heme a(3) is a good candidate to be the proton loading site.

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.

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.

Adaptation of aerobic respiration to low O2 environments

Aerobic respiration in bacteria, Archaea, and mitochondria is performed by oxygen reductase members of the heme-copper oxidoreductase superfamily. These enzymes are redox-driven proton pumps which conserve part of the free energy released from oxygen reduction to generate a proton motive force. The oxygen reductases can be divided into three main families based on evolutionary and structural analyses (A-, B- and C-families), with the B- and C-families evolving after the A-family. The A-family utilizes two proton input channels to transfer protons for pumping and chemistry, whereas the B- and C-families require only one. Generally, the B- and C-families also have higher apparent oxygen affinities than the A-family. Here we use whole cell proton pumping measurements to demonstrate differential proton pumping efficiencies between representatives of the A-, B-, and C-oxygen reductase families. The A-family has a coupling stoichiometry of 1 H(+)/e(-), whereas the B- and C-families have coupling stoichiometries of 0.5 H(+)/e(-). The differential proton pumping stoichiometries, along with differences in the structures of the proton-conducting channels, place critical constraints on models of the mechanism of proton pumping. Most significantly, it is proposed that the adaptation of aerobic respiration to low oxygen environments resulted in a concomitant reduction in energy conservation efficiency, with important physiological and ecological consequences.

High resolution structure of the ba3 cytochrome c oxidase from Thermus thermophilus in a lipidic environment

The fundamental chemistry underpinning aerobic life on Earth involves reduction of dioxygen to water with concomitant proton translocation. This process is catalyzed by members of the heme-copper oxidase (HCO) superfamily. Despite the availability of crystal structures for all types of HCO, the mode of action for this enzyme is not understood at the atomic level, namely how vectorial H⁺ and e^- transport are coupled. Toward addressing this problem, we report wild type and A120F mutant structures of the ba₃-type cytochrome c oxidase from Thermus thermophilus at 1.8 Å resolution. The enzyme has been crystallized from the lipidic cubic phase, which mimics the biological membrane environment. The structures reveal 20 ordered lipid molecules that occupy binding sites on the protein surface or mediate crystal packing interfaces. The interior of the protein encloses 53 water molecules, including 3 trapped in the designated K-path of proton transfer and 8 in a cluster seen also in A-type enzymes that likely functions in egress of product water and proton translocation. The hydrophobic O₂-uptake channel, connecting the active site to the lipid bilayer, contains a single water molecule nearest the Cu_B atom but otherwise exhibits no residual electron density. The active site contains strong electron density for a pair of bonded atoms bridging the heme Fe_a3 and Cu_B atoms that is best modeled as peroxide. The structure of ba₃-oxidase reveals new information about the positioning of the enzyme within the membrane and the nature of its interactions with lipid molecules. The atomic resolution details provide insight into the mechanisms of electron transfer, oxygen diffusion into the active site, reduction of oxygen to water, and pumping of protons across the membrane. The development of a robust system for production of ba₃-oxidase crystals diffracting to high resolution, together with an established expression system for generating mutants, opens the door for systematic structure-function studies.

Kinetic design of the respiratory oxidases

Energy conservation in all kingdoms of life involves electron transfer, through a number of membrane-bound proteins, associated with proton transfer across the membrane. In aerobic organisms, the last component of this electron-transfer chain is a respiratory heme-copper oxidase that catalyzes reduction of O(2) to H(2)O, linking this process to transmembrane proton pumping. So far, the molecular mechanism of proton pumping is not known for any system that is driven by electron transfer. Here, we show that this problem can be addressed and elucidated in a unique cytochrome c oxidase (cytochrome ba(3)) from a thermophilic bacterium, Thermus thermophilus. The results show that in this oxidase the electron- and proton-transfer reactions are orchestrated in time such that previously unresolved proton-transfer reactions could be directly observed. On the basis of these data we propose that loading of the proton pump occurs upon electron transfer, but before substrate proton transfer, to the catalytic site. Furthermore, the results suggest that the pump site alternates between a protonated and deprotonated state for every second electron transferred to the catalytic site, which would explain the noninteger pumping stoichiometry (0.5 H(+)/e(-)) of the ba(3) oxidase. Our studies of this variant of Nature's palette of mechanistic solutions to a basic problem offer a route toward understanding energy conservation in biological systems.

Stabilization of the peroxy intermediate in the oxygen splitting reaction of cytochrome cbb(3)

The proton-pumping cbb(3)-type cytochrome c oxidases catalyze cell respiration in many pathogenic bacteria. For reasons not yet understood, the apparent dioxygen (O(2)) affinity in these enzymes is very high relative to other members of the heme-copper oxidase (HCO) superfamily. Based on density functional theory (DFT) calculations on intermediates of the oxygen scission reaction in active-site models of cbb(3)- and aa(3)-type oxidases, we find that a transient peroxy intermediate (I(P), Fe[III]-OOH(-)) is ~6kcal/mol more stable in the former case, resulting in more efficient kinetic trapping of dioxygen and hence in a higher apparent oxygen affinity. The major molecular basis for this stabilization is a glutamate residue, polarizing the proximal histidine ligand of heme b(3) in the active site.

Modeling the MoFe nitrogenase system with broken symmetry density functional theory

Density functional theory (DFT) represents a unified framework for gaining molecular level insight into molybdenum-iron (MoFe) nitrogenase. However, accurately describing the electronic structure of the spin-polarized and spin-coupled iron-molybdenum cofactor (FeMo-co) where N(2) reduction occurs within MoFe nitrogenase is challenging. Therefore, the enhancement of DFT to include broken symmetry (BS-DFT) plus approximate spin projection has proven valuable because it provides a procedure to compute reliable geometries, energies, redox potentials, and quantities relevant to Mössbauer and ENDOR spectroscopies. After describing the theoretical tools necessary to obtain this information, we show by way of examples how BS-DFT is a very powerful partner to experiment. We expect that quantitative quantum chemical theory of this type will play an ever-increasing role in helping to decipher complex bioinorganic systems like those found in MoFe nitrogenase.