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

Long-range charge transfer mechanism of the III2IV2 mycobacterial supercomplex

Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III2IV2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria.

Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H2 gas powered by a [NiFe] active site and transduces the free energy into proton pumping and Na+/H+ exchange across the membrane. Despite recent structural advances, the mechanistic principles of H2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H2 and how the catalysis couples to conformational dynamics in the membrane domain of Mbh. By combining large-scale quantum chemical density functional theory (DFT) and correlated ab initio wave function methods with atomistic molecular dynamics simulations, we show that the proton transfer reactions required for the catalysis are gated by electric field effects that direct the protons by water-mediated reactions from Glu21L toward the [NiFe] site, or alternatively along the nearby His75L pathway that also becomes energetically feasible in certain reaction steps. These local proton-coupled electron transfer (PCET) reactions induce conformational changes around the active site that provide a key coupling element via conserved loop structures to the ion transport activity. We find that H2 forms in a heterolytic proton reduction step, with spin crossovers tuning the energetics along key reaction steps. On a general level, our work showcases the role of electric fields in enzyme catalysis and how these effects are employed by the [NiFe] active site of Mbh to drive PCET reactions and ion transport.

Radical Transport Facilitated by a Proton Transfer Network at the Subunit Interface of Ribonucleotide Reductase

Ribonucleotide reductases (RNRs) play an essential role in the conversion of nucleotides to deoxynucleotides in all organisms. The Escherichia coli class Ia RNR requires two homodimeric subunits, α and β. The active form is an asymmetric αα'ββ' complex. The α subunit houses the site for nucleotide reduction initiated by a thiyl radical (C439•), and the β subunit houses the diferric-tyrosyl radical (Y122•) that is essential for C439• formation. The reactions require a highly regulated and reversible long-range proton-coupled electron transfer pathway involving Y122•[β] ↔ W48?[β] ↔ Y356[β] ↔ Y731[α] ↔ Y730[α] ↔ C439[α]. In a recent cryo-EM structure, Y356[β] was revealed for the first time and it, along with Y731[α], spans the asymmetric α/β interface. An E52[β] residue, which is essential for Y356 oxidation, allows access to the interface and resides at the head of a polar region comprising R331[α], E326[α], and E326[α'] residues. Mutagenesis studies with canonical and unnatural amino acid substitutions now suggest that these ionizable residues are important in enzyme activity. To gain further insights into the roles of these residues, Y356• was photochemically generated using a photosensitizer covalently attached adjacent to Y356[β]. Mutagenesis studies, transient absorption spectroscopy, and photochemical assays monitoring deoxynucleotide formation collectively indicate that the E52[β], R331[α], E326[α], and E326[α'] network plays the essential role of shuttling protons associated with Y356 oxidation from the interface to bulk solvent.

Dissecting the Ultrafast Stepwise Bidirectional Proton Relay in a Blue-Light Photoreceptor

Proton relays through H-bond networks are essential in realizing the functionality of protein machines such as in photosynthesis and photoreceptors. It has been challenging to dissect the rates and energetics of individual proton-transfer steps during the proton relay. Here, we have designed a proton rocking blue light using a flavin (BLUF) domain with the flavin mononucleotide (FMN)-glutamic acid (E)-tryptophan (W) triad and have resolved the four individual proton-transfer steps kinetically using ultrafast spectroscopy. We have found that after the photo-induced charge separation forming FMN^·-/E-COOH/WH^·+, the proton first rapidly jumps from the bridging E-COOH to FMN^- (τ_fPT2 = 3.8 ps; KIE = 1.0), followed by a second proton transfer from WH^·+ to E-COO^- (τ_fPT1 = 336 ps; KIE = 2.6) which immediately rocks back to W^· (τ_rPT1 = 85 ps; KIE = 6.7), followed by a proton return from FMNH^· to E-COO^- (τ_rPT2 = 34 ps; KIE = 3.3) with the final charge recombination between FMN^·- and WH^·+ to close the reaction cycle. Our results revisited the Grotthuss mechanism on the ultrafast timescale using the BLUF domain as a paradigm protein.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Amino acid deprotonation rates from classical force fields

Acid ionization constants (pKa's) of titratable amino acid side chains have received a large amount of experimental and theoretical attention. In many situations, however, the rates of protonation and deprotonation, kon & koff, may also be important, for example in understanding the mechanism of action of proton channels or membrane proteins that couple proton transport to other processes. Protonation and deprotonation involve making and breaking of covalent bonds, which cannot be studied by classical force fields. However, environment effects on the rates should be captured by such methods. Here we present an approach for estimating deprotonation rates based on Warshel's extension of Marcus's theory of electron transfer, with input from molecular simulations. The missing bond dissociation energy is represented by a constant term determined by fitting the pKa value in solution. The statistics of the energy gap between protonated and deprotonated states is used to compute free energy curves of the two states and thus free energy barriers, from which the rate can be deduced. The method is applied to Glu, Asp, and His in bulk solution and select membrane proteins: the M2 proton channel, bacteriorhodopsin, and cytochrome c oxidase.

Structural basis for safe and efficient energy conversion in a respiratory supercomplex

Proton-translocating respiratory complexes assemble into supercomplexes that are proposed to increase the efficiency of energy conversion and limit the production of harmful reactive oxygen species during aerobic cellular respiration. Cytochrome bc complexes and cytochrome aa₃ oxidases are major drivers of the proton motive force that fuels ATP generation via respiration, but how wasteful electron- and proton transfer is controlled to enhance safety and efficiency in the context of supercomplexes is not known. Here, we address this question with the 2.8 Å resolution cryo-EM structure of the cytochrome bcc-aa₃ (III₂-IV₂) supercomplex from the actinobacterium Corynebacterium glutamicum. Menaquinone, substrate mimics, lycopene, an unexpected Q_c site, dioxygen, proton transfer routes, and conformational states of key protonable residues are resolved. Our results show how safe and efficient energy conversion is achieved in a respiratory supercomplex through controlled electron and proton transfer. The structure may guide the rational design of drugs against actinobacteria that cause diphtheria and tuberculosis.

Aerobic energy metabolism is driven by proton-pumping respiratory supercomplexes. The study reports the structural basis for energy conversion in such supercomplex. It may aid metabolic engineering and drug design against diphtheria and tuberculosis.

Functional Dynamics of an Ancient Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) is a redox-driven Na⁺/H⁺ transporter that employs the energy from hydrogen gas (H₂) production to catalyze proton pumping and Na⁺/H⁺ exchange across cytoplasmic membranes of archaea. Despite a recently resolved structure of this ancient energy-transducing enzyme [Yu et al. Cell2018, 173, 1636–1649], the molecular principles of its redox-driven ion-transport mechanism remain puzzling and of major interest for understanding bioenergetic principles of early cells. Here we use atomistic molecular dynamics (MD) simulations in combination with data clustering methods and quantum chemical calculations to probe principles underlying proton reduction as well as proton and sodium transport in Mbh from the hyperthermophilic archaeon Pyrococcus furiosus. We identify putative Na⁺ binding sites and proton pathways leading across the membrane and to the NiFe-active center as well as conformational changes that regulate ion uptake. We suggest that Na⁺ binding and protonation changes at a putative ion-binding site couple to proton transfer across the antiporter-like MbhH subunit by modulating the conformational state of a conserved ion pair at the subunit interface. Our findings illustrate conserved coupling principles within the complex I superfamily and provide functional insight into archaeal energy transduction mechanisms.

Beating Heart of Cytochrome c Oxidase: The Shared Proton of Heme a3 Propionates

Cytochrome c oxidase (CcO) pumps protons from the N-side to the P-side and consumes electrons from the P-side of the mitochondrial membrane driven by energy gained from reduction of dioxygen to water. ATP synthesis uses the resulting proton gradient and electrostatic potential difference. Since the distance a proton travels through CcO is too large for a one-step transfer process, proton-loading sites (PLS) that can carry protons transiently are necessary. One specific pump-active PLS couples to the redox reaction, thus energizing the proton to move across the membrane against electric potential and proton gradient. The PLS should also prevent proton backflow. Therefore, the propionates of the two redox-active hemes in CcO were suggested as PLS candidates although, according to CcO crystal structures, none of the four propionates can be protonated on account of strong H-bonds. Here, we show that modeling the local structure around heme a₃ propionates enhances significantly their capability of carrying a proton jointly. This was not possible for the propionates of heme a. The modeled structures are stable in molecular dynamics simulations (MDS) and are energetically similar to the crystal structure. Precise electrostatic energy computations of MDS data are used to estimate the pK_A values of all titratable residues in CcO. For the modeled structures, the heme a₃ propionates have pK_A values high enough to host a proton transiently but not too high to fix the proton permanently. The change in pK_A throughout the redox reaction is sufficient to push the proton to the P-side of the membrane and to provide the protons with the necessary amount of energy for ATP synthesis.

Electron Transfer Coupled to Conformational Dynamics in Cell Respiration

Cellular respiration is a fundamental process required for energy production in many organisms. The terminal electron transfer complex in mitochondrial and many bacterial respiratory chains is cytochrome c oxidase (CcO). This converts the energy released in the cytochrome c/oxygen redox reaction into a transmembrane proton electrochemical gradient that is used subsequently to power ATP synthesis. Despite detailed knowledge of electron and proton transfer paths, a central question remains as to whether the coupling between electron and proton transfer in mammalian mitochondrial forms of CcO is mechanistically equivalent to its bacterial counterparts. Here, we focus on the conserved span between H376 and G384 of transmembrane helix (TMH) X of subunit I. This conformationally-dynamic section has been suggested to link the redox activity with the putative H pathway of proton transfer in mammalian CcO. The two helix X mutants, Val380Met (V380M) and Gly384Asp (G384D), generated in the genetically-tractable yeast CcO, resulted in a respiratory-deficient phenotype caused by the inhibition of intra-protein electron transfer and CcO turnover. Molecular aspects of these variants were studied by long timescale atomistic molecular dynamics simulations performed on wild-type and mutant bovine and yeast CcOs. We identified redox- and mutation-state dependent conformational changes in this span of TMH X of bovine and yeast CcOs which strongly suggests that this dynamic module plays a key role in optimizing intra-protein electron transfers.

Architecture of bacterial respiratory chains

Bacteria power their energy metabolism using membrane-bound respiratory enzymes that capture chemical energy and transduce it by pumping protons or Na⁺ ions across their cell membranes. Recent breakthroughs in molecular bioenergetics have elucidated the architecture and function of many bacterial respiratory enzymes, although key mechanistic principles remain debated. In this Review, we present an overview of the structure, function and bioenergetic principles of modular bacterial respiratory chains and discuss their differences from the eukaryotic counterparts. We also discuss bacterial supercomplexes, which provide central energy transduction systems in several bacteria, including important pathogens, and which could open up possible avenues for treatment of disease.

Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations

Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O₂ reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R₄CO, R₂CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R₄CO to R₂CO. We performed potential jump experiments to accumulate the R₂CO and R₄CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a₃. Our experiments also suggest propionate A of heme a₃ changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site.

The Redox-Active Tyrosine Is Essential for Proton Pumping in Cytochrome c Oxidase

Cellular respiration involves electron transport via a number of enzyme complexes to the terminal Cytochrome c oxidase (CcO), in which molecular oxygen is reduced to water. The free energy released in the reduction process is used to establish a transmembrane electrochemical gradient, via two processes, both corresponding to charge transport across the membrane in which the enzymes are embedded. First, the reduction chemistry occurring in the active site of CcO is electrogenic, which means that the electrons and protons are delivered from opposite sides of the membrane. Second, the exergonic chemistry is coupled to translocation of protons across the entire membrane, referred to as proton pumping. In the largest subfamily of the CcO enzymes, the A-family, one proton is pumped for every electron needed for the chemistry, making the energy conservation particularly efficient. In the present study, hybrid density functional calculations are performed on a model of the A-family CcOs. The calculations show that the redox-active tyrosine, conserved in all types of CcOs, plays an essential role for the energy conservation. Based on the calculations a reaction mechanism is suggested involving a tyrosyl radical (possibly mixed with tyrosinate character) in all reduction steps. The result is that the free energy released in each reduction step is large enough to allow proton pumping in all reduction steps without prohibitively high barriers when the gradient is present. Furthermore, the unprotonated tyrosine provides a mechanism for coupling the uptake of two protons per electron in every reduction step, i.e. for a secure proton pumping.

Gated Proton Release during Radical Transfer at the Subunit Interface of Ribonucleotide Reductase

The class Ia ribonucleotide reductase of Escherichia coli requires strict regulation of long-range radical transfer between two subunits, α and β, through a series of redox-active amino acids (Y₁₂₂•[β] ↔ W₄₈?[β] ↔ Y₃₅₆[β] ↔ Y₇₃₁[α] ↔ Y₇₃₀[α] ↔ C₄₃₉[α]). Nowhere is this more precarious than at the subunit interface. Here, we show that the oxidation of Y₃₅₆ is regulated by proton release involving a specific residue, E₅₂[β], which is part of a water channel at the subunit interface for rapid proton transfer to the bulk solvent. An E₅₂Q variant is incapable of Y₃₅₆ oxidation via the native radical transfer pathway or non-native photochemical oxidation, following photosensitization by covalent attachment of a photo-oxidant at position 355[β]. Substitution of Y₃₅₆ for various F_nY analogues in an E₅₂Q-photoβ₂, where the side chain remains deprotonated, recovered photochemical enzymatic turnover. Transient absorption and emission data support the conclusion that Y₃₅₆ oxidation requires E₅₂ for proton management, suggesting its essential role in gating radical transport across the protein-protein interface.

Proton transfer in uncoupled variants of cytochrome c oxidase

Cytochrome c oxidase is a membrane-bound redox-driven proton pump that harbors two proton-transfer pathways, D and K, which are used at different stages of the reaction cycle. Here, we address the question if a D pathway with a modified energy landscape for proton transfer could take over the role of the K pathway when the latter is blocked by a mutation. Our data indicate that structural alterations near the entrance of the D pathway modulate energy barriers that influence proton transfer to the proton-loading site. The data also suggest that during reduction of the catalytic site, its protonation has to occur via the K pathway and that this proton transfer to the catalytic site cannot take place through the D pathway.

Explaining leak states in the proton pump of heme-copper oxidases observed in single-molecule experiments

Heme-copper oxidases couple the exergonic oxygen reduction with the endergonic proton translocation. Redox-linked structural changes have been localized in deeply buried regions of the protein, near the low-potential heme. How these movements can modulate distant gating events along the intramolecular proton path, where the entry (exit) of pumped proton occurs, is a major concern for the proton pump models. Generally, these models associate, more or less directly, all translocation events with redox transitions. Although they can account for many phenomenological aspects of the pump, evidences from single-molecules experiments about leak states of the pump represent a formidable challenge. Disconnecting the redox-linked pK_a shifts of the proton loading site from the external barriers, we obtain a simple stochastic mechanism which behaves similarly to the real enzyme, able to reverse the flow of the proton transfer.

Allosteric Cooperativity in Proton Energy Conversion in A1-Type Cytochrome c Oxidase

Cytochrome c oxidase (CcO), the Cu_A, heme a, heme a₃, Cu_B enzyme of respiratory chain, converts the free energy released by aerobic cytochrome c oxidation into a membrane electrochemical proton gradient (ΔμH⁺). ΔμH⁺ derives from the membrane anisotropic arrangement of dioxygen reduction to two water molecules and transmembrane proton pumping from a negative (N) space to a positive (P) space separated by the membrane. Spectroscopic, potentiometric, and X-ray crystallographic analyses characterize allosteric cooperativity of dioxygen binding and reduction with protonmotive conformational states of CcO. These studies show that allosteric cooperativity stabilizes the favorable conformational state for conversion of redox energy into a transmembrane ΔμH⁺.

Snapshot of an oxygen intermediate in the catalytic reaction of cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces dioxygen to water and harnesses the chemical energy to drive proton translocation across the inner mitochondrial membrane by an unresolved mechanism. By using time-resolved serial femtosecond crystallography, we identified a key oxygen intermediate of bovine CcO. It is assigned to the P_R-intermediate, which is characterized by specific redox states of the metal centers and a distinct protein conformation. The heme a ₃ iron atom is in a ferryl (Fe⁴⁺ = O^2-) configuration, and heme a and Cu_B are oxidized while Cu_A is reduced. A Helix-X segment is poised in an open conformational state; the heme a farnesyl sidechain is H-bonded to S382, and loop-I-II adopts a distinct structure. These data offer insights into the mechanism by which the oxygen chemistry is coupled to unidirectional proton translocation.

Mechanism of proton transfer through the KC proton pathway in the Vibrio cholerae cbb3 terminal oxidase

The heme‑copper oxidases (HCuOs) are terminal components of the respiratory chain, catalyzing oxygen reduction coupled to the generation of a proton motive force. The C-family HCuOs, found in many pathogenic bacteria under low oxygen tension, utilize a single proton uptake pathway to deliver protons both for O2 reduction and for proton pumping. This pathway, called the KC-pathway, starts at Glu-49P in the accessory subunit CcoP, and connects into the catalytic subunit CcoN via the polar residues Tyr-(Y)-227, Asn (N)-293, Ser (S)-244, Tyr (Y)-321 and internal water molecules, and continues to the active site. However, although the residues are known to be functionally important, little is known about the mechanism and dynamics of proton transfer in the KC-pathway. Here, we studied variants of Y227, N293 and Y321. Our results show that in the N293L variant, proton-coupled electron transfer is slowed during single-turnover oxygen reduction, and moreover it shows a pH dependence that is not observed in wildtype. This suggests that there is a shift in the pKa of an internal proton donor into an experimentally accessible range, from >10 in wildtype to ~8.8 in N293L. Furthermore, we show that there are distinct roles for the conserved Y321 and Y227. In Y321F, proton uptake from bulk solution is greatly impaired, whereas Y227F shows wildtype-like rates and retains ~50% turnover activity. These tyrosines have evolutionary counterparts in the K-pathway of B-family HCuOs, but they do not have the same roles, indicating diversity in the proton transfer dynamics in the HCuO superfamily.

Redox-Driven Proton Pumps of the Respiratory Chain

In aerobic cells, the proton gradient that drives ATP synthesis is created by three different proton pumps-membrane enzymes of the respiratory electron transport chain known as complex I, III, and IV. Despite the striking dissimilarity of structures and apparent differences in molecular mechanisms of proton pumping, all three enzymes have much in common and employ the same universal physical principles of converting redox energy to proton pumping. In this study, we describe a simple mathematical model that illustrates the general principles of redox-driven proton pumps and discuss their implementation in complex I, III, and IV of the respiratory chain.

Dewetting transitions coupled to K-channel activation in cytochrome c oxidase

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane.

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane. CcO employs the so-called D- and K-channels for proton uptake, but the molecular mechanism for activation of the K-channel has remained elusive for decades. We show here by combining large-scale atomistic molecular simulations with graph-theoretical water network analysis, and hybrid quantum/classical (QM/MM) free energy calculations, that the K-channel is activated by formation of a reactive oxidized intermediate in the binuclear heme a₃/Cu_B active site. This state induces electrostatic, hydration, and conformational changes that lower the barrier for proton transfer along the K-channel by dewetting pathways that connect the D-channel with the active site. Our combined results reconcile previous experimental findings and indicate that water dynamics plays a decisive role in the proton pumping machinery in CcO.

Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase (CcO) from the cellular inside to the binuclear redox centre as well as proton pumping through the membrane takes place through proton entrance via two distinct pathways, the D- and K-channel. Both channels show a dependence of their hydration level on the protonation states of their key residues, K362 for the K-channel, and E286 or D132 for the D-channel. In the oxidative half of CcO's catalytic cycle the D-channel is the proton-conducting path. For this channel, an interplay of protonation state of the D-channel residues with the water and hydrogen-bond dynamics has been observed in molecular dynamics simulations of the CcO protein, embedded in a lipid bi-layer, modelled in different protonation states. Protonation of residue E286 at the end of the D-channel results in a hydrogen-bonded network pointing from E286 to N139, that is against proton transport, and favouring N139 conformations which correspond to a closed asparagine gate (formed by residues N121 and N139). Consequently, the hydration level is lower than with unprotonated E286. In those models, the Asn gate is predominantly open, allowing water molecules to pass and thus increase the hydration level. The hydrogen-bonded network in these states exhibits longer life times of the Asn residues with water than other models and shows the D-channel to be traversable from the entrance, D132, to exit, E286. The D-channel can thus be regarded as auto-regulated with respect to proton transport, allowing proton passage only when required, that is the proton is located at the lower part of the D-channel (D132 to Asn gate) and not at the exit (E286).

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pK_a is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.

The reductive phase of Rhodobacter sphaeroides cytochrome c oxidase disentangled by CO ligation

Cytochrome c oxidase (CcO) is a membrane protein of the respiratory chain that catalytically reduces molecular oxygen (O2) to water while translocating protons across the membrane. The enzyme hosts two copper and two heme iron moieties (heme a/heme a3). The atomic details of the sequential steps that go along with this redox-driven proton translocation are a matter of debate. Particularly for the reductive phase of CcO that precedes oxygen binding experimental data are scarce. Here, we use CcO under anaerobic conditions where carbon monoxide (CO) is bound to heme a3 which in tandem with CuB forms the binuclear center (BNC). Fourier-transform infrared (FTIR) absorption spectroscopy is combined with electro-chemistry to probe different redox and protonation states populated by variation of the external electrostatic potential. With this approach, the redox behavior of heme a and the BNC could be separated and the corresponding redox potentials were determined. We also infer the protonation of one of the propionate side chains of heme a3 to correlate with the oxidation of heme a. Experimental changes in the local electric field surrounding CO bound to heme a3 are determined by their vibrational Stark effect and agree well with electrostatic computations. The comparison of experimental and computational results indicates that changes of the heme a3/CuB redox state are coupled to proton transfer towards heme a3. The latter supports the role of the heme a3 propionate D as proton loading site.

Cavity hydration dynamics in cytochrome c oxidase and functional implications

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O₂ reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O₂ reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively [Formula: see text]4 [Formula: see text]s) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a₃ and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of [Formula: see text]2-[Formula: see text]s trajectories suggests that hydration-level change occurs on the timescale of 100-200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.

Protonation-State-Dependent Communication in Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase from the cellular inside to the binuclear redox center (BNC) can occur through two distinct pathways, the D- and K-channels. For the protein to function as both a redox enzyme and a proton pump, proton transfer into the protein toward the BNC or toward a proton loading site (and ultimately through the membrane) must be highly regulated. The P_R → F transition is the first step in a catalytic cycle that requires proton transfer from the bulk at the N-side to the BNC. Molecular dynamics simulations of the P_R → F intermediate of this transition, with 16 different combinations of protonation states of key residues in the D- and K-channel, show the impact of the K-channel on the D-channel to be protonation-state dependent. Strength as well as means of communication, correlations in positions, or communication along the hydrogen-bonded network depends on the protonation state of the K-channel residue K362. The conformational and hydrogen-bond dynamics of the D-channel residue N139 is regulated by an interplay of protonation in the D-channel and K362. N139 thus assumes a gating function by which proton passage through the D-channel toward E286 is likely facilitated for states with protonated K362 and unprotonated E286. In contrast, proton passage through the D-channel is hindered by N139's preference for a closed conformation in situations with protonated E286.

Role of water and protein dynamics in proton pumping by respiratory complex I

Membrane bound respiratory complex I is the key enzyme in the respiratory chains of bacteria and mitochondria, and couples the reduction of quinone to the pumping of protons across the membrane. Recently solved crystal or electron microscopy structures of bacterial and mitochondrial complexes have provided significant insights into the electron and proton transfer pathways. However, due to large spatial separation between the electron and proton transfer routes, the molecular mechanism of coupling remains unclear. Here, based on atomistic molecular dynamics simulations performed on the entire structure of complex I from Thermus thermophilus, we studied the hydration of the quinone-binding site and the membrane-bound subunits. The data from simulations show rapid diffusion of water molecules in the protein interior, and formation of hydrated regions in the three antiporter-type subunits. An unexpected water-protein based connectivity between the middle of the Q-tunnel and the fourth proton channel is also observed. The protonation-state dependent dynamics of key acidic residues in the Nqo8 subunit suggest that the latter may be linked to redox-coupled proton pumping in complex I. We propose that in complex I the proton and electron transfer paths are not entirely separate, instead the nature of coupling may in part be ‘direct’.

Understanding the essential proton-pumping kinetic gates and decoupling mutations in cytochrome c oxidase

Cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water and uses the released free energy to pump protons against the transmembrane proton gradient. To better understand the proton-pumping mechanism of the wild-type (WT) CcO, much attention has been given to the mutation of amino acid residues along the proton translocating D-channel that impair, and sometimes decouple, proton pumping from the chemical catalysis. Although their influence has been clearly demonstrated experimentally, the underlying molecular mechanisms of these mutants remain unknown. In this work, we report multiscale reactive molecular dynamics simulations that characterize the free-energy profiles of explicit proton transport through several important D-channel mutants. Our results elucidate the mechanisms by which proton pumping is impaired, thus revealing key kinetic gating features in CcO. In the N139T and N139C mutants, proton back leakage through the D-channel is kinetically favored over proton pumping due to the loss of a kinetic gate in the N139 region. In the N139L mutant, the bulky L139 side chain inhibits timely reprotonation of E286 through the D-channel, which impairs both proton pumping and the chemical reaction. In the S200V/S201V double mutant, the proton affinity of E286 is increased, which slows down both proton pumping and the chemical catalysis. This work thus not only provides insight into the decoupling mechanisms of CcO mutants, but also explains how kinetic gating in the D-channel is imperative to achieving high proton-pumping efficiency in the WT CcO.

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

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

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.

A Protonated Water Cluster as a Transient Proton-Loading Site in Cytochrome c Oxidase

Cytochrome c oxidase (CcO) is a redox-driven proton pump that powers aerobic respiratory chains. We show here by multi-scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton-loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase

Monte Carlo (MC) simulations of conformational changes and protonation of Glu-242, a key residue that shuttles protons in cytochrome c oxidase (CcO), are reported. Previous studies suggest that this residue may play a role of the valve of the enzyme proton pump. Here we examine how sensitive the results of simulations are to the computational method used. We applied both molecular mechanic (MM) and hybrid quantum mechanic:molecular mechanic (QM:MM) methods and find that the results are qualitatively different. The results indicate that the mechanism for proton gating in CcO is still an open issue.

Molecular simulation and modeling of complex I

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

Single Enzyme Experiments Reveal a Long-Lifetime Proton Leak State in a Heme-Copper Oxidase

Heme-copper oxidases (HCOs) are key enzymes in prokaryotes and eukaryotes for energy production during aerobic respiration. They catalyze the reduction of the terminal electron acceptor, oxygen, and utilize the Gibbs free energy to transport protons across a membrane to generate a proton (ΔpH) and electrochemical gradient termed proton motive force (PMF), which provides the driving force for the adenosine triphosphate (ATP) synthesis. Excessive PMF is known to limit the turnover of HCOs, but the molecular mechanism of this regulatory feedback remains relatively unexplored. Here we present a single-enzyme study that reveals that cytochrome bo₃ from Escherichia coli, an HCO closely homologous to Complex IV in human mitochondria, can enter a rare, long-lifetime leak state during which proton flow is reversed. The probability of entering the leak state is increased at higher ΔpH. By rapidly dissipating the PMF, we propose that this leak state may enable cytochrome bo₃, and possibly other HCOs, to maintain a suitable ΔpH under extreme redox conditions.

Redox-induced activation of the proton pump in the respiratory complex I

Complex I functions as a redox-linked proton pump in the respiratory chains of mitochondria and bacteria, driven by the reduction of quinone (Q) by NADH. Remarkably, the distance between the Q reduction site and the most distant proton channels extends nearly 200 Å. To elucidate the molecular origin of this long-range coupling, we apply a combination of large-scale molecular simulations and a site-directed mutagenesis experiment of a key residue. In hybrid quantum mechanics/molecular mechanics simulations, we observe that reduction of Q is coupled to its local protonation by the His-38/Asp-139 ion pair and Tyr-87 of subunit Nqo4. Atomistic classical molecular dynamics simulations further suggest that formation of quinol (QH2) triggers rapid dissociation of the anionic Asp-139 toward the membrane domain that couples to conformational changes in a network of conserved charged residues. Site-directed mutagenesis data confirm the importance of Asp-139; upon mutation to asparagine the Q reductase activity is inhibited by 75%. The current results, together with earlier biochemical data, suggest that the proton pumping in complex I is activated by a unique combination of electrostatic and conformational transitions.

Structural Changes and Proton Transfer in Cytochrome c Oxidase

In cytochrome c oxidase electron transfer from cytochrome c to O2 is linked to transmembrane proton pumping, which contributes to maintaining a proton electrochemical gradient across the membrane. The mechanism by which cytochrome c oxidase couples the exergonic electron transfer to the endergonic proton translocation is not known, but it presumably involves local structural changes that control the alternating proton access to the two sides of the membrane. Such redox-induced structural changes have been observed in X-ray crystallographic studies at residues 423-425 (in the R. sphaeroides oxidase), located near heme a. The aim of the present study is to investigate the functional effects of these structural changes on reaction steps associated with proton pumping. Residue Ser425 was modified using site-directed mutagenesis and time-resolved spectroscopy was used to investigate coupled electron-proton transfer upon reaction of the oxidase with O2. The data indicate that the structural change at position 425 propagates to the D proton pathway, which suggests a link between redox changes at heme a and modulation of intramolecular proton-transfer rates.

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.

Respiratory conservation of energy with dioxygen: cytochrome C oxidase

Cytochrome c oxidase (CcO) is the terminal oxidase of cell respiration which reduces molecular oxygen (O₂) to H2O coupled with the proton pump. For elucidation of the mechanism of CcO, the three-dimensional location and chemical reactivity of each atom composing the functional sites have been extensively studied by various techniques, such as crystallography, vibrational and time-resolved electronic spectroscopy, since the X-ray structures (2.8 Å resolution) of bovine and bacterial CcO have been published in 1995.X-ray structures of bovine CcO in different oxidation and ligand binding states showed that the O₂reduction site, which is composed of Fe (heme a 3) and Cu (CuB), drives a non-sequential four-electron transfer for reduction of O₂to water without releasing any reactive oxygen species. These data provide the crucial structural basis to solve a long-standing problem, the mechanism of the O₂reduction.Time-resolved resonance Raman and charge translocation analyses revealed the mechanism for coupling between O₂reduction and the proton pump: O₂is received by the O₂reduction site where both metals are in the reduced state (R-intermediate), giving the O₂-bound form (A-intermediate). This is spontaneously converted to the P-intermediate, with the bound O₂fully reduced to 2 O²⁻. Hereafter the P-intermediate receives four electron equivalents from the second Fe site (heme a), one at a time, to form the three intermediates, F, O, and E to regenerate the R-intermediate. Each electron transfer step from heme a to the O₂reduction site is coupled with the proton pump.X-ray structural and mutational analyses of bovine CcO show three possible proton transfer pathways which can transfer pump protons (H) and chemical (water-forming) protons (K and D). The structure of the H-pathway of bovine CcO indicates that the driving force of the proton pump is the electrostatic repulsion between the protons on the H-pathway and positive charges of heme a, created upon oxidation to donate electrons to the O₂reduction site. On the other hand, mutational and time-resolved electrometric findings for the bacterial CcO strongly suggest that the D-pathway transfers both pump and chemical protons. However, the structure for the proton-gating system in the D-pathway has not been experimentally identified. The structural and functional diversities in CcO from various species suggest a basic proton pumping mechanism in which heme a pumps protons while heme a 3 reduces O₂as proposed in 1978.

Lysine 362 in cytochrome c oxidase regulates opening of the K-channel via changes in pKA and conformation

The metabolism of aerobic life uses the conversion of molecular oxygen to water as an energy source. This reaction is catalyzed by cytochrome e oxidase (CeO) consuming four electrons and four protons, which move along specific routes. While all four electrons are transferred via the same cofactors to the binuclear reaction center (BNC), the protons take two different routes in the A-type CeO, i.e., two of the four chemical protons consumed in the reaction arrive via the D-channel in the oxidative first half starting after oxygen binding. The other two chemical protons enter via the K-channel in the reductive second half of the reaction cycle. To date, the mechanism behind these separate proton transport pathways has not been understood. In this study, we propose a model that can explain the reaction-step specific opening and closing of the K-channel by conformational and pKA changes of its central lysine 362. Molecular dynamics simulations reveal an upward movement of Lys362 towards the BNC, which had already been supposed by several experimental studies. Redox state-dependent pKA calculations provide evidence that Lys362 may protonate transiently, thereby opening the K-channel only in the reductive second half of the reaction cycle. From our results, we develop a model that assigns a key role to Lys362 in the proton gating between the two proton input channels of the A-type CeO.

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.