Different types of biological proton transfer reactions studied by quantum chemical methods

Different types of proton transfer occurring in biological systems are described with examples mainly from ribonucleotide reductase (RNR) and cytochrome c oxidase (CcO). Focus is put on situations where electron and proton transfer are rather strongly coupled. In the long range radical transfer in RNR, it is shown that the presence of hydrogen atom transfer (HAT) is the most logical explanation for the experimental observations. In another example from RNR, it is shown that a transition state for concerted motion of both proton and electron can be found even if the donors are separated by a quite long distance. In CcO, the essential proton transfer for the OO bond cleavage, and the most recent modelings of proton translocation are described, indicating a few remaining major problems.     

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.     

Design of photoactive ruthenium complexes to study electron transfer and proton pumping in cytochrome oxidase

This review describes the development and application of photoactive ruthenium complexes to study electron transfer and proton pumping reactions in cytochrome c oxidase (CcO). CcO uses four electrons from Cc to reduce O(2) to two waters, and pumps four protons across the membrane. The electron transfer reactions in cytochrome oxidase are very rapid, and cannot be resolved by stopped-flow mixing techniques. Methods have been developed to covalently attach a photoactive tris(bipyridine)ruthenium group [Ru(II)] to Cc to form Ru-39-Cc. Photoexcitation of Ru(II) to the excited state Ru(II*), a strong reductant, leads to rapid electron transfer to the ferric heme group in Cc, followed by electron transfer to Cu(A) in CcO with a rate constant of 60,000s(-1). Ruthenium kinetics and mutagenesis studies have been used to define the domain for the interaction between Cc and CcO. New ruthenium dimers have also been developed to rapidly inject electrons into Cu(A) of CcO with yields as high as 60%, allowing measurement of the kinetics of electron transfer and proton release at each step in the oxygen reduction mechanism.     

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.     

A role for the protein in internal electron transfer to the catalytic center of cytochrome c oxidase

Internal electron transfer (ET) to heme a(3) during anaerobic reduction of oxidized bovine heart cytochrome c oxidase (CcO) was studied under conditions where heme a and Cu(A) were fully reduced by excess hexaamineruthenium. The data show that ET to heme a(3) is controlled by the state of ionization of a single protolytic residue with a pK(a) of 6.5 +/- 0.2. On the basis of the view that ET to the catalytic site is limited by coupled proton transfer, this pK(a) was attributed to Glu60 which is located at the entrance of the proton-conducting K channel on the matrix side of CcO. It is proposed that Glu60 controls proton entry into the channel. However, even with this channel open, there is the second factor that regulates ET, and this is ascribed to the rate of proton diffusion in the channel. In addition, it is concluded that proton transfer in the K channel is reversibly inhibited by the detergent Triton X-100. It is also found that the rate of ET to heme a(3) in the as-isolated resting enzyme and in CcO "activated" by reaction of fully reduced enzyme with O(2) is the same, implying that the catalytic sites of these two forms of oxidized enzyme are essentially identical.     

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.     

Transmembrane proton translocation by cytochrome c oxidase

Respiratory heme-copper oxidases are integral membrane proteins that catalyze the reduction of molecular oxygen to water using electrons donated by either quinol (quinol oxidases) or cytochrome c (cytochrome c oxidases, CcOs). Even though the X-ray crystal structures of several heme-copper oxidases and results from functional studies have provided significant insights into the mechanisms of O2 -reduction and, electron and proton transfer, the design of the proton-pumping machinery is not known. Here, we summarize the current knowledge on the identity of the structural elements involved in proton transfer in CcO. Furthermore, we discuss the order and timing of electron-transfer reactions in CcO during O2 reduction and how these reactions might be energetically coupled to proton pumping across the membrane.     

Assembly factors and ATP-dependent proteases in cytochrome c oxidase biogenesis

Eukaryotic cytochrome c oxidase (CcO), the terminal enzyme of the energy-transducing mitochondrial electron transport chain is a hetero-oligomeric, heme–copper oxidase complex composed of both mitochondrially and nuclear-encoded subunits. It is embedded in the inner mitochondrial membrane where it couples the transfer of electrons from reduced cytochrome c to molecular oxygen with vectorial proton translocation across the membrane. The biogenesis of CcO is a complicated sequential process that requires numerous specific accessory proteins, so-called assembly factors, which include translational activators, translocases, molecular chaperones, copper metallochaperones and heme a biosynthetic enzymes. Besides these CcO-specific protein factors, the correct biogenesis of CcO requires an even greater number of proteins with much broader substrate specificities. Indeed, growing evidence indicates that mitochondrial ATP-dependent proteases might play an important role in CcO biogenesis. Out of the four identified energy-dependent mitochondrial proteases, three were shown to be directly involved in proteolysis of CcO subunits. In addition to their well-established protein-quality control function these oligomeric proteolytic complexes with chaperone-like activities may function as molecular chaperones promoting productive folding and assembly of subunit proteins. In this review, we summarize the current knowledge of the functional involvement of eukaryotic CcO-specific assembly factors and highlight the possible significance for CcO biogenesis of mitochondrial ATP-dependent proteases.     

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.     

Alternative initial proton acceptors for the D pathway of Rhodobacter sphaeroides cytochrome c oxidase

To characterize protein structures that control proton uptake, we assayed forms of cytochrome c oxidase (CcO) containing a carboxyl or a thiol group in line with the initial, internal waters of the D pathway for proton transfer in the presence and absence of subunit III. Subunit III provides approximately half of the protein surrounding the entry region of the D pathway. The N139D/D132N mutant contains a carboxyl group 6 ? within the D pathway and lacks the normal, surface-exposed proton acceptor, Asp-132. With subunit III, the steady-state activity of this mutant is slow, but once subunit III is removed, its activity is the same as that of wild-type CcO lacking subunit III (～1800 H+/s). Thus, a carboxyl group～25% within the pathway enhances proton uptake even though the carboxyl has no direct contact with bulk solvent. Protons from solvent apparently move to internal Asp-139 through a short file of waters, normally blocked by subunit III. Cys-139 also supports rapid steady-state proton uptake, demonstrating that an anion other than a carboxyl can attract and transfer protons into the D pathway. When both Asp-132 and Asp/Cys-139 are present, the removal of subunit III increases CcO activity to rates greater than that of normal CcO because of simultaneous proton uptake by two initial acceptors. The results show how the environment of the initial proton acceptor for the D pathway in these CcO forms dictates the pH range of CcO activity, with implications for the function of Asp-132, the normal proton acceptor.     

The role of tryptophan 272 in the Paracoccus denitrificans cytochrome c oxidase

The mechanism of electron coupled proton transfer in cytochrome c oxidase (CcO) is still poorly understood. The P(M)-intermediate of the catalytic cycle is an oxoferryl state whose generation requires one additional electron, which cannot be provided by the two metal centres. The missing electron has been suggested to be donated to this binuclear site by a tyrosine residue. A tyrosine radical species has been detected in the P(M) and F* intermediates (formed by addition of H2O2) of the Paraccocus denitrificans CcO using electron paramagnetic resonance (EPR) spectroscopy. From the study of conserved variants its origin was determined to be Y167 which is surprising as this residue is not part of the active site. Upon inspection of the active site it becomes evident that W272 could be the actual donor of the missing electron, which can then be replenished from Y167 or from the Y280-H276 cross link in the natural cycle. To address the question, whether such a direct electron transfer pathway to the binuclear centre exists two tryptophan 272 variants in subunit I have been generated. These variants are characterised by their turnover rates as well as using EPR and optical spectroscopy. From these experiments it is concluded, that W272 is an important intermediate in the formation of the radical species appearing in P(M) and F* intermediates produced with hydrogen peroxide. The significance of this finding for the catalytic function of the enzyme is discussed.     

Structural studies on bovine heart cytochrome c oxidase

Among the X-ray structures of bovine heart cytochrome c oxidase (CcO), reported thus far, the highest resolution is 1.8?. CcO includes 13 different protein subunits, 7 species of phospholipids, 7 species of triglycerides, 4 redox-active metal sites (Cu(A), heme a (Fe(a)), Cu(B), heme a(3) (Fe(a3))) and 3 redox-inactive metal sites (Mg(2+), Zn(2+) and Na(+)). The effects of various O(2) analogs on the X-ray structure suggest that O(2) molecules are transiently trapped at the Cu(B) site before binding to Fe(a3)(2+) to provide O(2)(-). This provides three possible electron transfer pathways from Cu(B), Fe(a3) and Tyr244 via a water molecule. These pathways facilitate non-sequential 3 electron reduction of the bound O(2)(-) to break the OO bond without releasing active oxygen species. Bovine heart CcO has a proton conducting pathway that includes a hydrogen-bond network and a water-channel which, in tandem, connect the positive side phase with the negative side phase. The hydrogen-bond network forms two additional hydrogen-bonds with the formyl and propionate groups of heme a. Thus, upon oxidation of heme a, the positive charge created on Fe(a) is readily delocalized to the heme peripheral groups to drive proton-transport through the hydrogen-bond network. A peptide bond in the hydrogen-bond network and a redox-coupled conformational change in the water channel are expected to effectively block reverse proton transfer through the H-pathway. These functions of the pathway have been confirmed by site-directed mutagenesis of bovine CcO expressed in HeLa cells.     

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.     

Elementary steps of proton translocation in the catalytic cycle of cytochrome oxidase

Proton translocation in the catalytic cycle of cytochrome c oxidase (CcO) proceeds sequentially in a four-stroke manner. Every electron donated by cytochrome c drives the enzyme from one of four relatively stable intermediates to another, and each of these transitions is coupled to proton translocation across the membrane, and to uptake of another proton for production of water in the catalytic site. Using cytochrome c oxidase from Paracoccus denitrificans we have studied the kinetics of electron transfer and electric potential generation during several such transitions, two of which are reported here. The extent of electric potential generation during initial electron equilibration between CuA and heme a confirms that this reaction is not kinetically linked to vectorial proton transfer, whereas oxidation of heme a is kinetically coupled to the main proton translocation events during functioning of the proton pump. We find that the rates and amplitudes in multiphase heme a oxidation are different in the OH-->EH and PM-->F steps of the catalytic cycle, and that this is reflected in the kinetics of electric potential generation. We discuss this difference in terms of different driving forces and relate our results, and data from the literature, to proposed mechanisms of proton pumping in cytochrome c oxidase.     

Biogenesis of eukaryotic cytochrome c oxidase

Eukaryotic cytochrome c oxidase (CcO), the terminal component of the mitochondrial electron transport chain is a heterooligomeric complex that belongs to the superfamily of heme-copper containing terminal oxidases. The enzyme, composed of both mitochondrially and nuclear encoded subunits, is embedded in the inner mitochondrial membrane, where it catalyzes the transfer of electrons form reduced cytochrome c to dioxygen, coupling this reaction with vectorial proton pumping across the inner membrane. Due to the complexity of the enzyme, the biogenesis of CcO involves a multiplicity of steps, carried out by a number of highly specific gene products. These include mainly proteins that mediate the delivery and insertion of copper ions, synthesis and incorporation of heme moieties and membrane-insertion and topogenesis of constituent protein subunits. Isolated CcO deficiency represents one of the most frequently recognized causes of respiratory chain defects in humans, associated with severe, often fatal clinical phenotype. Here we review recent advancements in the understanding of this intricate process, with a focus on mammalian enzyme.     

Ribonucleotide reductase inhibition by p-alkoxyphenols studied by molecular docking and molecular dynamics simulations

Ribonucleotide reductase (RNR) is necessary for production of the precursor deoxyribonucleotides for DNA synthesis. Class Ia RNR functions via a stable free radical in one of the two components protein R2. The enzyme mechanism involves long range (proton coupled) electron transfer between protein R1 and the tyrosyl radical in protein R2. Earlier experimental studies showed that p-alkoxyphenols inhibit RNR. Here, molecular docking and molecular dynamics simulations involving protein R2 suggest an inhibition mechanism for p-alkoxyphenols . A low energy binding pocket is identified in protein R2. The preferred configuration provides a structural basis explaining their specific binding to the Escherichia coli and mouse R2 proteins. Trp48 (E. coli numbering), on the electron transfer pathway, is involved in the interactions with the inhibitors. The relative order of the binding energies calculated for the phenol derivatives to protein R2 is correlated with earlier experimental data on inhibition efficiency, in turn related to increasing size of the hydrophobic alkyl substituents. Using the configuration identified by molecular docking as a starting point for molecular dynamics simulations, we find that the p-allyloxyphenol interrupts the catalytic electron transfer pathway of the R2 protein by forming hydrogen bonds with Trp48 and Asp237, thus explaining the inhibitory activity of p-alkoxyphenols.     

Subunit III-depleted cytochrome c oxidase provides insight into the process of proton uptake by proteins

We review studies of subunit III-depleted cytochrome c oxidase (CcO III (-)) that elucidate the structural basis of steady-state proton uptake from solvent into an internal proton transfer pathway. The removal of subunit III from R. sphaeroides CcO makes proton uptake into the D pathway a rate-determining step, such that measurements of the pH dependence of steady-state O(2) consumption can be used to compare the rate and functional pK(a) of proton uptake by D pathways containing different initial proton acceptors. The removal of subunit III also promotes spontaneous suicide inactivation by CcO, greatly shortening its catalytic lifespan. Because the probability of suicide inactivation is controlled by the rate at which the D pathway delivers protons to the active site, measurements of catalytic lifespan provide a second method to compare the relative efficacy of proton uptake by engineered CcO III (-) forms. These simple experimental systems have been used to explore general questions of proton uptake by proteins, such as the functional value of an initial proton acceptor, whether an initial acceptor must be surface-exposed, which side chains will function as initial proton acceptors and whether multiple acceptors can speed proton uptake.     

2D-SEIRA spectroscopy to highlight conformational changes of the cytochrome c oxidase induced by direct electron transfer

Potentiometric titrations of the cytochrome c oxidase (CcO) immobilized in a biomimetic membrane system were followed by two-dimensional surface-enhanced IR absorption spectroscopy (2D SEIRAS) in the ATR-mode. Direct electron transfer was employed to vary the redox state of the enzyme. The CcO was shown to undergo a conformational transition from a non-activated to an activated state after it was allowed to turnover in the presence of oxygen. Differences between the non-activated and activated state were revealed by 2D SEIRA spectra recorded as a function of potential. The activated state was characterized by a higher number of correlated transitions as well as a higher number of amino acids associated with electron transfer.     

Proton interactions with hemes a and a3 in bovine heart cytochrome c oxidase

Three forms of cytochrome c oxidase, fully oxidized CcO (CcO-O), oxidized CcO complexed with cyanide (CcO.CN), and mixed valence CcO, in which both heme a(3) and Cu(B) are reduced and stabilized by carbon monoxide (MV.CO), were investigated by optical spectroscopy, MCD, and stopped-flow for the pH sensitivity of spectral features. In the pH range between pH 5.7 and 9.0, both heme a and heme a(3) in CcO-O interact with a single protolytic group. From the variation of the position of the Soret peak with changes in pH, a pK(a) of 6.6 +/- 0.2 was determined for this group. The pH sensitivity of heme a(3) is lost in the CcO.CN complex, and only heme a responds to pH changes. In MV.CO the spectra of both hemes are almost independent of pH between 5.7 and 11.0. The stoichiometry of proton uptake in the conversion of CcO-O both to MV.CO and to fully reduced CcO was determined between pH 5.8 and pH 8.2. Formation of MV.CO from CcO-O was accompanied by the uptake of approximately two protons, and this value was almost independent of pH. Full reduction of oxidized CcO was associated with the uptake of approximately 2 H(+) at basic pH, and this value increases with decreasing pH. On the basis of these proton uptake measurements, it is concluded that the pK(a) of the group is independent of the redox state of CcO. It is suggested that Glu60 of subunit II, located at the entrance of the proton conducting K-channel, is the protolytic residue that interacts with both hemes through a hydrogen-bonding network.