Exploring pathways and barriers for coupled ET/PT in cytochrome c oxidase: a general framework for examining energetics and mechanistic alternatives

Gaining a detailed understanding of the energetics of the proton pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Although there are several current mechanistic proposals, most of these ideas have not been subjected to consistent structure-function considerations. In particular most works have not related the activation barriers for different mechanistic proposals to the protein structure. The present work describes a general approach for exploring the energetics of different feasible models of the action of CcO, using the observed protein structure, established simulation methods and a modified Marcus' formulation. We start by reviewing our methods for evaluation of the energy diagrams for different proton translocation paths and then present a systematic analysis of various constraints that should be imposed on any energy diagram for the pumping process. After the general analysis we turn to the actual computational study, where we construct energy diagrams for forward and backward paths, using the estimated calculated reduction potentials and pK(a) values of all the relevant sites (including internal water molecules). We then explore the relationship between the calculated energy diagrams and key experimental constraints. This comparison allows us to identify some barriers that are not fully consistent with the overall requirement for an efficient pumping. In particular we identify back leakage channels, which are hard to block without stopping the forward channels. This helps to identify open problems that will require further experimental and theoretical studies. We also consider reasonable adjustments of the calculated barriers that may lead to a working pump. Although the present analysis does not establish a unique and workable model for the mechanism of CcO, it presents what is probably the most consistent current analysis of the barriers for different feasible pathways. Perhaps more importantly, the framework developed here should provide a general way for examining any proposal for the action of CcO as well as for the analysis of further experimental findings about the action of this fascinating system.     

Exploration of the cytochrome c oxidase pathway puzzle and examination of the origin of elusive mutational effects

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is a problem of great current interest. Despite promising mechanistic proposals, so far, a physically consistent model that would reproduce all the relevant barriers needed to create a working pump has not been presented. In addition, there are major problems in elucidating the origin of key mutational effects and in understanding the nature of the apparent pK(a) values associated with the pH dependencies of specific proton transfer (PT) reactions in CcO. This work takes a key step in resolving the above problems, by considering mutations, such as the Asn139Asp replacement, that blocks proton pumping without affecting PT to the catalytic site. We first introduce a formulation that makes it possible to relate the apparent pK(a) of Glu286 to different conformational states of this residue. We then use the new formulation along with the calculated pK(a) values of Glu286 at these different conformations to reproduce the experimentally observed apparent pK(a) of the residue. Next, we take the X-ray structures of the native and Asn139Asp mutant of the Paracoccus denitrificans CcO (N131D in this system) and reproduce for the first time the change in the primary PT pathways (and other key features) based on simulations that start with the observed structural changes. We also consider the competition between proton transport to the catalytic site and the pump site, as a function of the bulk pH, as well as the H/D isotope effect, and use this information to explore the relative height of the two barriers. The paper emphasizes the crucial role of energy-based considerations that include the PT process, and the delicate control of PT in CcO.     

Energy diagrams and mechanism for proton pumping in cytochrome c oxidase

The powerful technique of energy diagrams has been used to analyze the mechanism for proton pumping in cytochrome c oxidase. Energy levels and barriers are derived starting out from recent kinetic experiments for the O to E transition, and are then refined using general criteria and a few additional experimental facts. Both allowed and non-allowed pathways are obtained in this way. A useful requirement is that the forward and backward rate should approach each other for the full membrane gradient. A key finding is that an electron on heme a (or the binuclear center) must have a significant lowering effect on the barrier for proton uptake, in order to prevent backflow from the pump-site to the N-side. While there is no structural gating in the present mechanism, there is thus an electronic gating provided by the electron on heme a. A quantitative analysis of the energy levels in the diagrams, leads to Prop-A of heme a₃ as the most likely position for the pump-site, and the Glu278 region as the place for the transition state for proton uptake. Variations of key redox potentials and pK_a values during the pumping process are derived for comparison to experiments.     

Understanding the cytochrome c oxidase proton pump: thermodynamics of redox linkage.

The cytochrome c oxidase complex (CcO) catalyzes the four-electron reduction of dioxygen to water by using electrons from ferrocytochrome c. Redox free energy released in this highly exergonic process is utilized to drive the translocation of protons across a transmembrane electrochemical gradient. Although numerous chemical models of proton pumping have been developed, few attempts have been made to explain the stepwise transfer of energy in the context of proposed protein conformational changes. A model is described that seeks to clarify the thermodynamics of the proton pumping function of CcO and that illustrates the importance of electron and proton gating to prevent the occurrence of the more exergonic electron leak and proton slip reactions. The redox energy of the CcO-membrane system is formulated in terms of a multidimensional energy surface projected into two dimensions, a nuclear coordinate associated with electron transfer and a nuclear coordinate associated with elements of the proton pump. This model provides an understanding of how a transmembrane electrochemical gradient affects the efficiency of the proton pumping process. Specifically, electron leak and proton slip reactions become kinetically viable as a result of the greater energy barriers that develop for the desired reactions in the presence of a transmembrane potential.     

Network analysis of a proposed exit pathway for protons to the P-side of cytochrome c oxidase

Cytochrome c Oxidase (CcO) reduces O₂, the terminal electron acceptor, to water in the aerobic, respiratory electron transport chain. The energy released by O₂ reductions is stored by removing eight protons from the high pH, N-side, of the membrane with four used for chemistry in the active site and four pumped to the low pH, P-side. The proton transfers must occur along controllable proton pathways that prevent energy dissipating movement towards the N-side. The CcO N-side has well established D- and K-channels to deliver protons to the protein interior. The P-side has a buried core of hydrogen-bonded protonatable residues designated the Proton Loading Site cluster (PLS cluster) and many protonatable residues on the P-side surface, providing no obvious unique exit. Hydrogen bond pathways were identified in Molecular Dynamics (MD) trajectories of Rb. sphaeroides CcO prepared in the P_R state with the heme a₃ propionate and Glu286 in different protonation states. Grand Canonical Monte Carlo sampling of water locations, polar proton positions and residue protonation states in trajectory snapshots identify a limited number of water mediated, proton paths from PLS cluster to the surface via a (P-exit) cluster of residues. Key P-exit residues include His93, Ser168, Thr100 and Asn96. The hydrogen bonds between PLS cluster and P-exit clusters are mediated by a water wire in a cavity centered near Thr100, whose hydration can be interrupted by a hydrophobic pair, Leu255B (near Cu_A) and Ile99. Connections between the D channel and PLS via Glu286 are controlled by a second, variably hydrated cavity.

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's Cu_B center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model—the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from δ-nitrogen of His291 ligand of Cu_B center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pK_a of His291 for an extended model in which the whole Fe_a3-Cu_B center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pK_a units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.     

Expanding the view of proton pumping in cytochrome c oxidase through computer simulation

In cytochrome c oxidase (CcO), a redox-driven proton pump, protons are transported by the Grotthuss shuttling via hydrogen-bonded water molecules and protonatable residues. Proton transport through the D-pathway is a complicated process that is highly sensitive to alterations in the amino acids or the solvation structure in the channel, both of which can inhibit proton pumping and enzymatic activity. Simulations of proton transport in the hydrophobic cavity showed a clear redox state dependence. To study the mechanism of proton pumping in CcO, multi-state empirical valence bond (MS-EVB) simulations have been conducted, focusing on the proton transport through the D-pathway and the hydrophobic cavity next to the binuclear center. The hydration structures, transport pathways, effects of residues, and free energy surfaces of proton transport were revealed in these MS-EVB simulations. The mechanistic insight gained from them is herein reviewed and placed in context for future studies.     

From static structure to living protein: computational analysis of cytochrome c oxidase main-chain flexibility

Crystallographic structure and deuterium accessibility comparisons of CcO in different redox states have suggested conformational changes of mechanistic significance. To predict the intrinsic flexibility and low energy motions in CcO, this work has analyzed available high-resolution crystallographic structures with ProFlex and elNémo computational methods. The results identify flexible regions and potential conformational changes in CcO that correlate well with published structural and biochemical data and provide mechanistic insights. CcO is predicted to undergo rotational motions on the interior and exterior of the membrane, driven by transmembrane helical tilting and bending, coupled with rocking of the β-sheet domain. Consequently, the proton K-pathway becomes sufficiently flexible for internal water molecules to alternately occupy upper and lower parts of the pathway, associated with conserved Thr-359 and Lys-362 residues. The D-pathway helices are found to be relatively rigid, with a highly flexible entrance region involving the subunit I C-terminus, potentially regulating the uptake of protons. Constriction and dilation of hydrophobic channels in RsCcO suggest regulation of the oxygen supply to the binuclear center. This analysis points to coupled conformational changes in CcO and their potential to influence both proton and oxygen access.     

Quantum chemistry as a tool in bioenergetics

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

Energy diagrams and mechanism for proton pumping in cytochrome c oxidase

The powerful technique of energy diagrams has been used to analyze the mechanism for proton pumping in cytochrome c oxidase. Energy levels and barriers are derived starting out from recent kinetic experiments for the O to E transition, and are then refined using general criteria and a few additional experimental facts. Both allowed and non-allowed pathways are obtained in this way. A useful requirement is that the forward and backward rate should approach each other for the full membrane gradient. A key finding is that an electron on heme a (or the binuclear center) must have a significant lowering effect on the barrier for proton uptake, in order to prevent backflow from the pump-site to the N-side. While there is no structural gating in the present mechanism, there is thus an electronic gating provided by the electron on heme a. A quantitative analysis of the energy levels in the diagrams, leads to Prop-A of heme a(3) as the most likely position for the pump-site, and the Glu278 region as the place for the transition state for proton uptake. Variations of key redox potentials and pK(a) values during the pumping process are derived for comparison to experiments.     

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's CuB center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model-the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from delta-nitrogen of His291 ligand of CuB center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pKa of His291 for an extended model in which the whole Fe(a3)-CuB center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pKa units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.     

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.     

Two conformational states of Glu242 and pKas in bovine cytochrome c oxidase.

Cytochrome c oxidase (CcO) is the terminal enzyme in the respiratory electron transport chain of aerobic organisms. It catalyses the reduction of atmospheric oxygen to water, and couples this reaction to proton pumping across the membrane; this process generates the electrochemical gradient that subsequently drives the synthesis of ATP. The molecular details of the mechanism by which electron transfer is coupled to proton pumping in CcO is poorly understood. Recent calculations from our group indicate that His291, a ligand of the Cu(B) center of the enzyme, may play the role of the pumping element. In this paper we describe calculations in which a DFT/continuum electrostatic method is used to explore the coupling of the conformational changes of Glu242 residue, the main proton donor of both chemical and pump protons, to its pKa, and the pKa of His291, a putative proton loading site of our pumping model. The computations are done for several redox states of metal centers, different protonation states of Glu242 and His291, and two well-defined conformations of the Glu242 side chain. Thus, in addition to equilibrium redox/protonation states of the catalytic cycle, we also examine the transient and intermediate states. Different dielectric models are employed to investigate the robustness of the results, and their viability in the light of the proposed proton pumping mechanism of CcO. The main results are in agreement with the experimental measurements and support the proposed pumping mechanism. Additionally, the present calculations indicate a possibility of gating through conformational changes of Glu242; namely, in the pumping step, we find that Glu242 needs to be reprotonated before His291 can eject a proton to the P-site of membrane. As a result, the reprotonation of Glu can control proton release from the proton loading site.     

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.     

Critical roles of the CuB site in efficient proton pumping as revealed by crystal structures of mammalian cytochrome c oxidase catalytic intermediates

Mammalian cytochrome c oxidase (CcO) reduces O₂ to water in a bimetallic site including Fe_a3 and Cu_B giving intermediate molecules, termed A-, P-, F-, O-, E-, and R-forms. From the P-form on, each reaction step is driven by single-electron donations from cytochrome c coupled with the pumping of a single proton through the H-pathway, a proton-conducting pathway composed of a hydrogen-bond network and a water channel. The proton-gradient formed is utilized for ATP production by F-ATPase. For elucidation of the proton pumping mechanism, crystal structural determination of these intermediate forms is necessary. Here we report X-ray crystallographic analysis at ∼1.8 Å resolution of fully reduced CcO crystals treated with O₂ for three different time periods. Our disentanglement of intermediate forms from crystals that were composed of multiple forms determined that these three crystallographic data sets contained ∼45% of the O-form structure, ∼45% of the E-form structure, and ∼20% of an oxymyoglobin-type structure consistent with the A-form, respectively. The O- and E-forms exhibit an unusually long Cu_B²⁺-OH⁻ distance and Cu_B¹⁺-H₂O structure keeping Fe_a3³⁺-OH⁻ state, respectively, suggesting that the O- and E-forms have high electron affinities that cause the O→E and E→R transitions to be essentially irreversible and thus enable tightly coupled proton pumping. The water channel of the H-pathway is closed in the O- and E-forms and partially open in the R-form. These structures, together with those of the recently reported P- and F-forms, indicate that closure of the H-pathway water channel avoids back-leaking of protons for facilitating the effective proton pumping.     

The protonation state of a heme propionate controls electron transfer in cytochrome c oxidase

In cytochrome c oxidase (CcO), exergonic electron transfer reactions from cytochrome c to oxygen drive proton pumping across the membrane. Elucidation of the proton pumping mechanism requires identification of the molecular components involved in the proton transfer reactions and investigation of the coupling between internal electron and proton transfer reactions in CcO. While the proton-input trajectory in CcO is relatively well characterized, the components of the output pathway have not been identified in detail. In this study, we have investigated the pH dependence of electron transfer reactions that are linked to proton translocation in a structural variant of CcO in which Arg481, which interacts with the heme D-ring propionates in a proposed proton output pathway, was replaced with Lys (RK481 CcO). The results show that in RK481 CcO the midpoint potentials of hemes a and a(3) were lowered by approximately 40 and approximately 15 mV, respectively, which stabilizes the reduced state of Cu(A) during reaction of the reduced CcO with O(2). In addition, while the pH dependence of the F --> O rate in wild-type CcO is determined by the protonation state of two protonatable groups with pK(a) values of 6.3 and 9.4, only the high-pK(a) group influences this rate in RK481 CcO. The results indicate that the protonation state of the Arg481 heme a(3) D-ring propionate cluster having a pK(a) of approximately 6.3 modulates the rate of internal electron transfer and may act as an acceptor of pumped protons.     

Cytochrome c oxidase loses catalytic activity and structural integrity during the aging process in Drosophila melanogaster

The hypothesis, that structural deterioration of cytochrome c oxidase (CcO) is a causal factor in the age-related decline in mitochondrial respiratory activity and an increase in H₂O₂ generation, was tested in Drosophila melanogaster. CcO activity and the levels of seven different nuclear DNA-encoded CcO subunits were determined at three different stages of adult life, namely, young-, middle-, and old-age. CcO activity declined progressively with age by 33%. Western blot analysis, using antibodies specific to Drosophila CcO subunits IV, Va, Vb, VIb, VIc, VIIc, and VIII, indicated that the abundance these polypeptides decreased, ranging from 11% to 40%, during aging. These and previous results suggest that CcO is a specific intra-mitochondrial site of age-related deterioration, which may have a broad impact on mitochondrial physiology.     

The proton-pumping site of cytochrome c oxidase: a model of its structure and mechanism

Cytochrome c oxidase is an electron-transfer driven proton pump. In this paper, we propose a complete chemical mechanism for the enzyme's proton-pumping site. The mechanism achieves pumping with chemical reaction steps localized at a redox center within the enzyme; no indirect coupling through protein conformational changes is required. The proposed mechanism is based on a novel redox-linked transition metal ligand substitution reaction. The use of this reaction leads in a straightforward manner to explicit mechanisms for achieving all of the processes previously determined (Blair, D.F., Gelles, J. and Chan, S.I. (1986) Biophys. J. 50, 713-733) to be needed to accomplish redox-linked proton pumping. These processes include: (1) modulation of the energetics of protonation/deprotonation reactions and modulation of the energetics of redox reactions by the structural state of the pumping site; (2) control of the rates of the pump's redox reactions with its electron-transfer partners during the turnover cycle (gating of electrons); and (3) regulation of the rates of the protonation/deprotonation reactions between the pumping site and the aqueous phases on the two sides of the membrane during the reaction cycle (gating of protons). The model is the first proposed for the cytochrome oxidase proton pump which is mechanistically complete and sufficiently specific that a realistic assessment can be made of how well the model pump would function as a redox-linked free-energy transducer. This assessment is accomplished via analyses of the thermodynamic properties and steady-state kinetics expected of the model. These analyses demonstrate that the model would function as an efficient pump and that its behavior would be very similar to that observed of cytochrome oxidase both in the mitochondrion and in purified preparations. The analysis presented here leads to the following important general conclusions regarding the mechanistic features of the oxidase proton pump. (1) A workable proton-pump mechanism does not require large protein conformational changes. (2) A redox-linked proton pump need not display a pH-dependent midpoint potential, as has frequently been assumed. (3) Mechanisms for redox-linked proton pumps that involve transition metal ligand exchange reactions are quite attractive because such reactions readily lend themselves to the linked gating processes necessary for proton pumping.     

Structural and energetic determinants of primary proton transfer in bacteriorhodopsin.

In the light-driven bacteriorhodopsin proton pump, the first proton transfer step is from the retinal Schiff base to a nearby carboxylate group. The mechanism of this transfer step is highly controversial, in particular whether a direct proton jump is allowed. Here, we review the structural and energetic determinants of the direct proton transfer path computed by using a combined quantum mechanical/molecular mechanical approach. Both protein flexibility and electrostatic interactions play an important role in shaping the proton transfer energy profile. Detailed analysis of the energetics of putative transitions in the first half of the photocycle focuses on two elements that determine the likelihood that a given configuration of the active site is populated during the proton-pumping cycle. First, the rate-limiting barrier for proton transfer must be consistent with the kinetics of the photocycle. Second, the active-site configuration must be compatible with a productive overall pumping cycle.     

Molecular dynamics and structural models of the cyanobacterial NDH-1 complex

NDH-1 is a gigantic redox-driven proton pump linked with respiration and cyclic electron flow in cyanobacterial cells. Based on experimentally resolved X-ray and cryo-EM structures of the respiratory complex I, we derive here molecular models of two isoforms of the cyanobacterial NDH-1 complex involved in redox-driven proton pumping (NDH-1L) and CO₂-fixation (NDH-1MS). Our models show distinct structural and dynamic similarities to the core architecture of the bacterial and mammalian respiratory complex I. We identify putative plastoquinone-binding sites that are coupled by an electrostatic wire to the proton pumping elements in the membrane domain of the enzyme. Molecular simulations suggest that the NDH-1L isoform undergoes large-scale hydration changes that support proton-pumping within antiporter-like subunits, whereas the terminal subunit of the NDH-1MS isoform lacks such structural motifs. Our work provides a putative molecular blueprint for the complex I-analogue in the photosynthetic energy transduction machinery and demonstrates that general mechanistic features of the long-range proton-pumping machinery are evolutionary conserved in the complex I-superfamily.