Water-gated mechanism of proton translocation by cytochrome c oxidase

The O(2) reduction and proton pumping gate mechanism of bovine heart cytochrome c oxidase

X-ray structures of bovine heart cytochrome c oxidase with bound respiratory inhibitors (O(2) analogues) have been determined at 1.8-2.05Å resolution to investigate the function of the O(2) reduction site which includes two metal sites (Fe(a3)(2+) and Cu(B)(1+)). The X-ray structures of the CO- and NO-bound derivatives indicate that although there are three possible electron donors that can provide electrons to the bound O(2), located in the O(2) reduction site, the formation of the peroxide intermediate is effectively prevented to provide an O(2)-bound form as the initial intermediate. The structural change induced upon binding of CN(-) suggests a non-sequential 3-electron reduction of the bound O(2)(-) for the complete reduction without release of any reactive oxygen species. The X-ray structure of the derivative with CO bound to Cu(B)(1+) after photolysis from Fe(a3)(2+) demonstrates weak side-on binding. This suggests that Cu(B) controls the O(2) supply to Fe(a3)(2+) without electron transfer to provide sufficient time for collection of protons from the negative side of the mitochondrial membrane. The proton-pumping pathway of bovine heart cytochrome c oxidase includes a hydrogen-bond network and a water channel located in tandem between the positive and negative side of the mitochondrial membrane. Binding of a strong ligand to Fe(a3) induces a conformational change which significantly narrows the water channel and effectively blocks the back-leakage of protons from the hydrogen bond network. The proton pumping mechanism proposed by these X-ray structural analyses has been functionally confirmed by mutagenesis analyses of bovine heart cytochrome c oxidase. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Water in the orchestration of the cell machinery. Some misunderstandings: a short review

Nowadays, biologists can explore the cell at the nanometre level. They discover an unsuspected world, amazingly overcrowded, complex and heterogeneous, in which water, also, is complex and heterogeneous. In the cell, statistical phenomena, such as diffusion, long considered as the main transport for water soluble substances, must be henceforth considered as inoperative to orchestrate cell activity. Results at this level are not yet numerous enough to give an exact representation of the cell machinery; however, they are sufficient to cease reasoning in terms of statistics (diffusion, law of mass action, pH, etc.) and encourage cytologists and biochemists to prospect thoroughly the huge panoply of the biophysical properties of macromolecule-water associations at the nanometre level. Our main purpose, here, is to discuss some of the more common misinterpretations due to the ignorance of these properties, and expose briefly the bases for a better approach of the cell machinery. Giorgio Careri, who demonstrated the correlation between proton currents at the surface of lysozyme and activity of this enzyme was one of the pioneers of this approach.

The D-channel of cytochrome oxidase: an alternative view

The D-pathway in A-type cytochrome c oxidases conducts protons from a conserved aspartate on the negatively charged N-side of the membrane to a conserved glutamic acid at about the middle of the membrane dielectric. Extensive work in the past has indicated that all four protons pumped across the membrane on reduction of O(2) to water are transferred via the D-pathway, and that it is also responsible for transfer of two out of the four "chemical protons" from the N-side to the binuclear oxygen reduction site to form product water. The function of the D-pathway has been discussed in terms of an apparent pK(a) of the glutamic acid. After reacting fully reduced enzyme with O(2), the rate of formation of the F state of the binuclear heme-copper active site was found to be independent of pH up to pH~9, but to drop off at higher pH with an apparent pK(a) of 9.4, which was attributed to the glutamic acid. Here, we present an alternative view, according to which the pH-dependence is controlled by proton transfer into the aspartate residue at the N-side orifice of the D-pathway. We summarise experimental evidence that favours a proton pump mechanism in which the proton to be pumped is transferred from the glutamic acid to a proton-loading site prior to proton transfer for completion of oxygen reduction chemistry. The mechanism is discussed by which the proton-pumping activity is decoupled from electron transfer by structural alterations of the D-pathway. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Water molecule reorganization in cytochrome c oxidase revealed by FTIR spectroscopy

Although internal electron transfer and oxygen reduction chemistry in cytochrome c oxidase are fairly well understood, the associated groups and pathways that couple these processes to gated proton translocation across the membrane remain unclear. Several possible pathways have been identified from crystallographic structural models; these involve hydrophilic residues in combination with structured waters that might reorganize to form transient proton transfer pathways during the catalytic cycle. To date, however, comparisons of atomic structures of different oxidases in different redox or ligation states have not provided a consistent answer as to which pathways are operative or the details of their dynamic changes during catalysis. In order to provide an experimental means to address this issue, FTIR spectroscopy in the 3,560-3,800 cm(-1) range has been used to detect weakly H-bonded water molecules in bovine cytochrome c oxidase that might change during catalysis. Full redox spectra exhibited at least four signals at 3,674(+), 3,638(+), 3,620(-), and 3,607(+) cm(-1). A more complex set of signals was observed in spectra of photolysis of the ferrous-CO compound, a reaction that mimics the catalytic oxygen binding step, and their D(2)O and H(2)(18)O sensitivities confirmed that they arose from water molecule rearrangements. Fitting with Gaussian components indicated the involvement of up to eight waters in the photolysis transition. Similar signals were also observed in photolysis spectra of the ferrous-CO compound of bacterial CcO from Paracoccus denitrificans. Such water changes are discussed in relation to roles in hydrophilic channels and proton/electron coupling mechanism.

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.

Spectral components of the α-band of cytochrome oxidase

Oxidative redox titrations of the mitochondrial cytochromes were performed in near-anoxic RAW 264.7 cells by inhibiting complex I. Cytochrome oxidation changes were measured with multi-wavelength spectroscopy and the ambient redox potential was calculated from the oxidation state of endogenous cytochrome c. Two spectral components were separated in the α-band range of cytochrome oxidase and they were identified as the difference spectrum of heme a when it has a high (a(H)) or low (a(L)) midpoint potential (E(m)) by comparing their occupancy during redox titrations carried out when the membrane potential (ΔΨ) was dissipated with a protonophore to that predicted by the neoclassical model of redox cooperativity. The difference spectrum of a(L) has a maximum at 605nm whereas the spectrum of a(H) has a maximum at 602nm. The ΔΨ-dependent shift in the E(m) of a(H) was too great to be accounted for by electron transfer from cytochrome c to heme a against ΔΨ but was consistent with a model in which a(H) is formed after proton uptake against ΔΨ suggesting that the spectral changes are the result of protonation. A stochastic simulation was implemented to model oxidation states, proton uptake and E(m) changes during redox titrations. The redox anti-cooperativity between heme a and heme a(3), and proton binding, could be simulated with a model where the pump proton interacted with heme a and the substrate proton interacted with heme a(3) with anti-cooperativity between proton binding sites, but not with a single proton binding site coupled to both hemes.

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.

Simulating proton transport through a simplified model for trans-membrane proteins

Ab initio MD simulations on a polyglycine helix and water-wire expressed under periodic boundary conditions have created a channel that supports proton transfer up to distances of 10.5 A. The effect of varying the density of water molecules in the channel has been investigated. A range of cationic states are identified with widely varying lifetimes. The mechanism of proton transport in this model shares some features with the simulations reported for bulk water, with, e.g., the hydrogen bond distance shortening in the time period leading up to successful proton transfer. However, there are also some important differences such as the observation of a heightened number of proton rattling events. We also observe that the helix plays an important role in directing the behavior of the water wire: the most active proton transport regions of the water-wire are found in areas where the helix is most tightly coiled. Finally, we report on the effects of different DFT functionals to model a water-wire and on the importance of including dispersion corrections to stabilize the alpha-helical structure.

Partial steps of charge translocation in the nonpumping N139L mutant of Rhodobacter sphaeroides cytochrome c oxidase with a blocked D-channel

The N139L substitution in the D-channel of cytochrome oxidase from Rhodobacter sphaeroides results in an approximately 15-fold decrease in the turnover number and a loss of proton pumping. Time-resolved absorption and electrometric assays of the F --> O transition in the N139L mutant oxidase result in three major findings. (1) Oxidation of the reduced enzyme by O(2) shows approximately 200-fold inhibition of the F --> O step (k approximately 2 s(-1) at pH 8) which is not compatible with enzyme turnover ( approximately 30 s(-1)). Presumably, an abnormal intermediate F(deprotonated) is formed under these conditions, one proton-deficient relative to a normal F state. In contrast, the F --> O transition in N139L oxidase induced by single-electron photoreduction of intermediate F, generated by reaction of the oxidized enzyme with H(2)O(2), decelerates to an extent compatible with enzyme turnover. (2) In the N139L mutant, the protonic phase of Deltapsi generation coupled to the flash-induced F --> O transition greatly decreases in rate and magnitude and can be assigned to the movement of a proton from E286 to the binuclear site, required for reduction of heme a(3) from the Fe(4+) horizontal lineO(2-) state to the Fe(3+)-OH(-) state. Electrogenic reprotonation of E286 from the inner aqueous phase is missing from the F --> O step in the mutant. (3) In the N139L mutant, the KCN-insensitive rapid electrogenic phase may be composed of two components with lifetimes of approximately 10 and approximately 40 mus and a magnitude ratio of approximately 3:2. The 10 mus phase matches vectorial electron transfer from Cu(A) to heme a, whereas the 40 mus component is assigned to intraprotein proton displacement across approximately 20% of the membrane dielectric thickness. This proton displacement might be triggered by rotation of the charged K362 side chain coupled to heme a reduction. The two components of the rapid electrogenic phase have been resolved subsequently with other D-channel mutants as well as with cyanide-inhibited wild-type oxidase. The finding helps to reconcile the unusually high relative contribution of the microsecond electrogenic phase in the bacterial enzyme ( approximately 30%) with the net electrogenicity of the F --> O transition coupled to transmembrane transfer of two charges per electron.

Structure and organization of mitochondrial respiratory complexes: a new understanding of an old subject

The enzymatic complexes of the mitochondrial respiratory chain have been extensively investigated in their structural and functional properties. A clear distinction is possible today between three complexes in which the difference in redox potential allows proton translocation (complexes I, III, and IV) and those having the mere function to convey electrons to the respiratory chain. We also have a clearer understanding of the structure and function of most respiratory complexes, of their biogenesis and regulation, and of their capacity to generate reactive oxygen species. Past investigations led to the conclusion that the complexes are randomly dispersed and functionally connected by diffusion of smaller redox components, coenzyme Q and cytochrome c. More-recent investigations by native gel electrophoresis and single-particle image processing showed the existence of supramolecular associations. Flux-control analysis demonstrated that complexes I and III in mammals and I, III, and IV in plants kinetically behave as single units, suggesting the existence of substrate channeling. This review discusses conditions affecting the formation of supercomplexes that, besides kinetic advantage, have a role in the stability and assembly of the individual complexes and in preventing excess oxygen radical formation. Disruption of supercomplex organization may lead to functional derangements responsible for pathologic changes.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

Quantum-chemical calculations of a long proton wire. Application of a harmonic model to analysis of the structure of an ionic defect in a water chain with an excess proton

Quantum-chemical calculations of molecular complexes (NH(3))(3)Zn(2+)...(H(2)O)(n)...NH(3) (C(n), n = 11, 16, 21, and 30) simulating a proton wire donor-water chain-acceptor were carried out. Earlier found periodicity in the length of the O-H bonds in water chain is explained within the framework of a one-component harmonic model. In complexes C(n), the geometry and electronic structure of ionic defect in water chain with an excess proton were studied. Calculations carried out at ab initio (B3LYP/6-31+G**) and semiempirical (PM3) levels of theory predict different patterns of distribution of the O-H bonds lengths and positive charge on the H-bond hydrogen atoms in the region of ionic defect. The obtained data show how a length of water chain and position of a protonated water link in the chain influence the ionic defect structure. To describe the observed structures of ionic defect, the harmonic model was used and the role of parameters of the H-bonded chain was investigated. The performed analysis explains different mechanisms (concerted and stepwise) of proton transfer along the H-bonded chain derived from ab initio and semiempirical calculation schemes.

Communication between R481 and Cu(B) in cytochrome bo(3) ubiquinol oxidase from Escherichia coli

The R481 residue of cytochrome bo(3) ubiquinol oxidase from E. coli is highly conserved in the heme-copper oxidase superfamily. It has been postulated to serve as part of a proton loading site that regulates proton translocation across the protein matrix of the enzyme. Along these lines, proton pumping efficiency has been demonstrated to be abolished in many R481 mutants. However, R481Q in bo(3) from E. coli has been shown to be fully functional, implying that the positive charge of the arginine is not required for proton translocation [ Puustinen , A. and Wikstrom , M. ( 1999 ) Proc. Natl. Acad. Sci. U.S.A. 96 , 35 - 37 ]. In an effort to delineate the structural role of R481 in the bo(3) oxidase, we used resonance Raman spectroscopy to compare the nonfunctional R481L mutant and the functional R481Q mutant, to the wild type protein. Resonance Raman data of the oxidized and reduced forms of the R481L mutant indicate that the mutation introduces changes to the heme o(3) coordination state, reflecting a change in position and/or coordination of the Cu(B) located on the distal side of heme o(3), although it is approximately 10 A away from R481. In the reduced-CO adduct of R481L, the frequencies of the Fe-CO and C-O stretching modes indicate that, unlike the wild type protein, the Cu(B) is no longer close to the heme-bound CO. In contrast, resonance Raman data obtained from the various oxidation and ligation states of the R481Q mutant are similar to those of the wild type protein, except that the mutation causes an enhancement of the relative intensity of the beta conformer of the CO-adduct, indicating a shift in the equilibrium between the alpha and beta conformers. The current findings, together with crystallographic structural data of heme-copper oxidases, indicate that R481 plays a keystone role in stabilizing the functional structure of the Cu(B) site through a hydrogen bonding network involving ordered water molecules. The implications of these data on the proton translocation mechanism are considered.

Modulation of the active site conformation by site-directed mutagenesis in cytochrome c oxidase from Paracoccus denitrificans

The structural and functional properties of active site mutants of cytochrome c oxidase from Paracoccus denitrificans (PdCcO) were investigated with resonance Raman spectroscopy. Based on the Fe-CO stretching modes and low frequency heme modes, two conformers (alpha- and beta-forms) were identified that are in equilibrium in the enzyme. The alpha-conformer, which is the dominant species in the wild-type enzyme, has a shorter heme a(3) iron-Cu(B) distance and a more distorted heme, as compared to the beta-conformer, which has a more relaxed and open distal pocket. In general, the mutations caused a decrease in the population of the alpha-conformer, which is concomitant with a decreased in the catalytic activity, indicating that the alpha-conformer is the active form of the enzyme. The data suggest that the native structure of the enzyme is in a delicate balance of intramolecular interactions. We present a model in which the mutations destabilize the alpha-conformer, with respect to the beta-conformer, and raise the activation barrier for the inter-conversion between the two conformers. The accessibility of the two conformers in the conformational space of CcO plausibly plays a critical role in coupling the redox reaction to proton translocation during the catalytic cycle of the enzyme.

Properties of Arg481 mutants of the aa3-type cytochrome c oxidase from Rhodobacter sphaeroides suggest that neither R481 nor the nearby D-propionate of heme a3 is likely to be the proton loading site of the proton pump

Cytochrome c oxidase utilizes the energy from electron transfer and reduction of oxygen to water and pumps protons across the membrane, generating a proton motive force. A large body of biochemical work has shown that all the pumped protons enter the enzyme through the D-channel, which is apparent in X-ray structures as a chain of water molecules connecting D132 at the cytoplasmic surface of the enzyme to E286, near the enzyme active site. The exit pathway utilized by pumped protons beyond this point and leading to the bacterial periplasm is not known. Also not known is the proton loading site (or sites) which undergoes changes in pKa in response to the chemistry at the enzyme active site and drives the proton pump mechanism. In this paper, we examine the role of R481, a highly conserved arginine that forms an ion pair with the D-propionate of heme a3. The R481H, R481N, R481Q, and R481L mutants were examined. The R481H mutant oxidase is approximately 18% active and pumps protons with approximately 40% of the stoichiometry of the wild type. The R481N, R481Q, and R481L mutants each retain only approximately 5% of the steady-state activity, and this is shown to be due to inhibition of steps in the reaction of O(2) with the reduced enzyme. Neither the R481N mutant nor the R481Q mutant oxidases pump protons, but remarkably, the R481L mutant does pump protons with the same efficiency as the R481H mutant. Since the proton pump is clearly operating in the R481L mutant, these results rule out an essential role in the proton pump mechanism for R481 or its hydrogen bond partner, the D-propionate of heme a3.

Proton transport pathway in the ClC Cl-/H+ antiporter

A fundamental question concerning the ClC Cl-/H+ antiporters is the nature of their proton transport (PT) pathway. We addressed this issue by using a novel computational methodology capable of describing the explicit PT dynamics in the ClC-ec1 protein. The main result is that the Glu203 residue delivers a proton from the intracellular solution to the core of ClC-ec1 via a rotation of its side chain and subsequent acid dissociation. After reorientation of the Glu203 side chain, a transient water-mediated PT pathway between Glu203 and Glu148 is established that is able to receive and translocate the proton via Grotthuss shuttling after deprotonation of Glu203. A molecular-dynamics simulation of an explicit hydrated excess proton in this pathway suggests that a negatively charged Glu148 and the central Cl- ion act together to drive H+ to the extracellular side of the membrane. This finding is consistent with the experimental result that Cl- binding to the central site facilitates the proton movement. A calculation of the PT free-energy barrier for the ClC-ec1 E203V mutant also supports the proposal that a dissociable residue is required at this position for efficient delivery of H+ to the protein interior, in agreement with recent experimental results.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

Redox-dependent conformational changes in cytochrome C oxidase suggest a gating mechanism for proton uptake

A role for conformational change in the coupling mechanism of cytochrome c oxidase is the subject of controversy. Relatively small conformational changes have been reported in comparisons of reduced and oxidized crystal structures of bovine oxidase but none in bacterial oxidases. Comparing the X-ray crystal structures of the reduced (at 2.15 A resolution) and oxidized forms of cytochrome c oxidase from Rhodobacter sphaeroides, we observe a displacement of heme a(3) involving both the porphyrin ring and the hydroxyl farnesyl tail, accompanied by protein movements in nearby regions, including the mid part of helix VIII of subunit I which harbors key residues of the K proton uptake path, K362 and T359. The conformational changes in the reduced form are reversible upon reoxidation. They result in an opening of the top of the K pathway and more ordered waters being resolved in that region, suggesting an access path for protons into the active site. In all high-resolution structures of oxidized R. sphaeroides cytochrome c oxidase, a water molecule is observed in the hydrophobic region above the top of the D path, strategically positioned to facilitate the connection of residue E286 of subunit I to the active site or to the proton pumping exit path. In the reduced and reduced plus cyanide structures, this water molecule disappears, implying disruption of proton conduction from the D path under conditions when the K path is open, thus providing a mechanism for alternating access to the active site.

The chemistry of the CuB site in cytochrome c oxidase and the importance of its unique His-Tyr bond

The CuB metal center is at the core of the active site of the heme-copper oxidases, comprising a copper atom ligating three histidine residues one of which is covalently bonded to a tyrosine residue. Using quantum chemical methodology, we have studied the CuB site in several redox and ligand states proposed to be intermediates of the catalytic cycle. The importance of the His-Tyr crosslink was investigated by comparing energetics, charge, and spin distributions between systems with and without the crosslink. The His-Tyr bond was shown to decrease the proton affinity and increase the electron affinity of both Tyr-244 and the copper. A previously unnoticed internal electronic equilibrium between the copper atom and the tyrosine was observed, which seems to be coupled to the unique structure of the system. In certain states the copper and Tyr-244 compete for the unpaired electron, the localization of which is determined by the oxygenous ligand of the copper. This electronic equilibrium was found to be sensitive to the presence of a positive charge 10 A away from the center, simulating the effect of Lys-319 in the K-pathway of proton transfer. The combined results provide an explanation for why the heme-copper oxidases need two pathways of proton uptake, and why the K-pathway is active only in the second half of the reaction cycle.

Molecular dynamics simulation of water in cytochrome c oxidase reveals two water exit pathways and the mechanism of transport

We have examined the network of connected internal cavities in cytochrome c oxidase along which water produced at the catalytic center is removed from the enzyme. Using combination of structural analysis, molecular dynamics simulations, and free energy calculations we have identified two exit pathways that connect the Mg2+ ion cavity to the outside of the enzyme. Each pathway has a well-defined bottleneck, which determines the overall rate of water traffic along the exit pathway, and a specific cooperative mechanism of passing it. One of the pathways is going via Arg438/439 (in bovine numbering) toward the CuA center, approaching closely its His204B ligand and Lys171B residue; and the other is going toward Asp364 and Thr294. Comparison of the pathways among different aa3-type enzymes shows that they are well conserved. Possible connections of the finding to redox-coupled proton pumping mechanism are discussed. We propose specific mutations near the bottlenecks of the exit pathways that can test some of our hypotheses.

Electron transfer and energy transduction in the terminal part of the respiratory chain - lessons from bacterial model systems

This review focuses on the terminal part of the respiratory chain where, macroscopically speaking, electron transfer (ET) switches from the two-electron donor, ubiquinol, to the single-electron carrier, cytochrome c, to finally reduce the four-electron acceptor dioxygen. With 3-D structures of prominent representatives of such multi-subunit membrane complexes known for some time, this section of the ET chain still leaves a number of key questions unanswered. The two relevant enzymes, ubiquinol:cytochrome c oxidoreductase and cytochrome c oxidase, appear as rather diverse modules, differing largely in their design for substrate interaction, internal ET, and moreover, in their mechanisms of energy transduction. While the canonical mitochondrial complexes have been investigated for almost five decades, the corresponding bacterial enzymes have been established only recently as attractive model systems to address basic reactions in ET and energy transduction. Lacking the intricate coding background and mitochondrial assembly pathways, bacterial respiratory enzymes typically offer a much simpler subunit composition, while maintaining all fundamental functions established for their complex "relatives". Moreover, related issues ranging from primary steps in cofactor insertion to supramolecular architecture of ET complexes, can also be favourably addressed in prokaryotic systems to hone our views on prototypic structures and mechanisms common to all family members.

Mechanisms of proton transfer in proteins: localized charge transfer versus delocalized soliton transfer

Proton translocation coupled to redox chemistry is ubiquitous for membrane enzymes involved in energy generation in cells. In such enzymes, proton transport occurs in special proton conducting channels, which consist of a series of protonatable groups of the protein connected by chains of mobile water molecules. Here we discuss two possible mechanisms of proton transport along such structures: diffusion of a localized charge and delocalized soliton transitions, in which several protons are collectively shifted along a chain of hydrogen bonds.

A mitochondrial DNA mutation linked to colon cancer results in proton leaks in cytochrome c oxidase

An increasing number of cancer types have been found to be linked to specific mutations in the mitochondrial DNA, which result in specific structural changes of the respiratory enzyme complexes. In this study, we have investigated the effect of 2 such mutations identified in colon cancer patients, leading to the amino acid substitutions Ser458Pro and Gly125Asp in subunit I of cytochrome c oxidase (complex IV) [Greaves et al. (2006) Proc Natl Acad Sci USA 103:714-719]. We introduced these mutations in Rhodobacter sphaeroides, which carries an oxidase that serves as a model of the mitochondrial counterpart. The lack of expression of the former variant indicates that the amino acid substitution results in severely altered overall structure of the enzyme. The latter mutation (Gly171Asp in the bacterial oxidase) resulted in a structurally intact enzyme, but with reduced activity (approximately 30%), mainly due to slowed reduction of the redox site heme a. Furthermore, even though the Gly171Asp CytcO pumps protons, an intrinsic proton leak was identified, which would lead to a decreased overall energy-conversion efficiency of the respiratory chain, and would also perturb transport processes such as protein, ion, and metabolite trafficking. Furthermore, the specific leak may act to alter the balance between the electrical and chemical components of the proton electrochemical gradient.

Cytochrome c oxidase: exciting progress and remaining mysteries

Cytochrome c oxidase generates a proton motive force by two separate mechanisms. The first mechanism is similar to that postulated by Peter Mitchell, and is based on electrons and protons used to generate water coming from opposite sides of the membrane. The second mechanism was not initially anticipated, but is now firmly established as a proton pump. A brief review of the current state of our understanding of the proton pump of cytochrome oxidase is presented. We have come a long way since the initial observation of the pump by Mårten Wikström in 1977, but a number of essential questions remain to be answered.

Inhibition of bacterial oxidases by formamide and analogs

The enzymatic activity of Paracoccus denitrificans cytochrome c oxidase (COX) and Escherichia coli cytochrome b(o) ubiquinol oxidase (QOX) was determined in the presence of formamide, N,N-dimethyl formamide and N,N-dimethyl acetamide. Formamide was found to inhibit the enzyme activity of the oxidases most significantly, whereas the other two compounds inhibited the activity to a lesser extent. The effects of formamide and analogs on enzyme activity were very similar for COX and QOX, indicating that the mechanism of inhibition might be the same for both of these oxidases. The inhibition kinetics followed a non-competitive mechanism. Studies using proteoliposomes of COX and QOX containing the electron entry site of the enzyme directed towards the outside of the vesicles showed that the effect of inhibition by formamide was higher when the inhibitor was present on the outside of the proteoliposome compared to when it was present only in the aqueous core. This indicates that inhibition of enzyme activity by formamide possibly predominantly involves blocking of the water exit pathway in the oxidases.

Toward a chemical mechanism of proton pumping by the B-type cytochrome c oxidases: application of density functional theory to cytochrome ba3 of Thermus thermophilus

A mechanism for proton pumping by the B-type cytochrome c oxidases is presented in which one proton is pumped in conjunction with the weakly exergonic, two-electron reduction of Fe-bound O 2 to the Fe-Cu bridging peroxodianion and three protons are pumped in conjunction with the highly exergonic, two-electron reduction of Fe(III)- (-)O-O (-)-Cu(II) to form water and the active oxidized enzyme, Fe(III)- (-)OH,Cu(II). The scheme is based on the active-site structure of cytochrome ba 3 from Thermus thermophilus, which is considered to be both necessary and sufficient for coupled O 2 reduction and proton pumping when appropriate gates are in place (not included in the model). Fourteen detailed structures obtained from density functional theory (DFT) geometry optimization are presented that are reasonably thought to occur during the four-electron reduction of O 2. Each proton-pumping step takes place when a proton resides on the imidazole ring of I-His376 and the large active-site cluster has a net charge of +1 due to an uncompensated, positive charge formally associated with Cu B. Four types of DFT were applied to determine the energy of each intermediate, and standard thermochemical approaches were used to obtain the reaction free energies for each step in the catalytic cycle. This application of DFT generally conforms with previously suggested criteria for a valid model (Siegbahn, P. E. M.; Blomberg, M. A. R. Chem. Rev. 2000, 100, 421-437) and shows how the chemistry of O 2 reduction in the heme a 3 -Cu B dinuclear center can be harnessed to generate an electrochemical proton gradient across the lipid bilayer.

Uncovering channels in photosystem II by computer modelling: current progress, future prospects, and lessons from analogous systems

Even prior to the publication of the crystal structures for photosystem II (PSII), it had already been suggested that water, O(2) and H(+) channels exist in PSII to achieve directed transport of these molecules, and to avoid undesirable side reactions. Computational efforts to uncover these channels and investigate their properties are still at early stages, and have so far only been based on the static PSII structure. The rationale behind the proposals for such channels and the computer modelling studies thus far are reviewed here. The need to take the dynamic protein into account is then highlighted with reference to the specific issues and techniques applicable to the simulation of each of the three channels. In particular, lessons are drawn from simulation studies on other protein systems containing similar channels.

Charge parameterization of the metal centers in cytochrome c oxidase

Reliable atomic point charges are of key importance for a correct description of the electrostatic interactions when performing classical, force field based simulations. Here, we present a systematic procedure for point charge derivation, based on quantum mechanical methodology suited for the systems at hand. A notable difference to previous procedures is to include an outer region around the actual system of interest. At the cost of increasing the system sizes, here up to 265 atoms, including the surroundings achieves near-neutrality for the systems as well as structural stability, important factors for reliable charge distributions. In addition, the common problem of converting between C--H bonds and C--C bonds at the border vanishes. We apply the procedure to the four redox-active metal centers of cytochrome c oxidase: Cu(A), haem a, haem a(3), and Cu(B). Several relevant charge and ligand states are considered. Charges for two different force fields, CHARMM and AMBER, are presented.

Ion exclusion by sub-2-nm carbon nanotube pores

Biological pores regulate the cellular traffic of a large variety of solutes, often with high selectivity and fast flow rates. These pores share several common structural features: the inner surface of the pore is frequently lined with hydrophobic residues, and the selectivity filter regions often contain charged functional groups. Hydrophobic, narrow-diameter carbon nanotubes can provide a simplified model of membrane channels by reproducing these critical features in a simpler and more robust platform. Previous studies demonstrated that carbon nanotube pores can support a water flux comparable to natural aquaporin channels. Here, we investigate ion transport through these pores using a sub-2-nm, aligned carbon nanotube membrane nanofluidic platform. To mimic the charged groups at the selectivity region, we introduce negatively charged groups at the opening of the carbon nanotubes by plasma treatment. Pressure-driven filtration experiments, coupled with capillary electrophoresis analysis of the permeate and feed, are used to quantify ion exclusion in these membranes as a function of solution ionic strength, pH, and ion valence. We show that carbon nanotube membranes exhibit significant ion exclusion that can be as high as 98% under certain conditions. Our results strongly support a Donnan-type rejection mechanism, dominated by electrostatic interactions between fixed membrane charges and mobile ions, whereas steric and hydrodynamic effects appear to be less important.

Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidase

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a(3) (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that "bridges" of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

Looking for the minimum common denominator in haem-copper oxygen reductases: towards a unified catalytic mechanism

Haem-copper oxygen reductases are transmembrane protein complexes that reduce dioxygen to water and pump protons across the mitochondrial or periplasmatic membrane, contributing to the transmembrane difference of electrochemical potential. Seven years ago we proposed a classification of these enzymes into three different families (A, B and C), based on the amino acid residues of their proton channels and amino acid sequence comparison, later supported by the so far identified characteristics of the catalytic centre of members from each family. The three families have in common the same general structural fold of the catalytic subunit, which contains the same or analogous prosthetic groups, and proton channels. These observations raise the hypothesis that the mechanisms for dioxygen reduction, proton pumping and the coupling of the two processes may be the same for all these enzymes. Under this hypothesis, they should be performed and controlled by the same or equivalent elements/events, and the identification of retained elements in all families will reveal their importance and may prompt the definition of the enzyme operating mode. Thus, we believe that the search for a minimum common denominator has a crucial importance, and in this article we highlight what is already established for the haem-copper oxygen reductases and emphasize the main questions still unanswered in a comprehensive basis.

Carboxyl group functions in the heme-copper oxidases: information from mid-IR vibrational spectroscopy

Carboxyl groups of possible functional importance in bovine and bacterial cytochrome c oxidases (CcO) are reviewed and assessed. A critical analysis is presented of available mid-infrared vibrational data that pertain to these functional carboxyl groups. These data and their interpretations are discussed in relation to current models of the mechanism of proton and electron coupling in the protonmotive CcO superfamily.

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

Aerobic life is based on a molecular machinery that utilizes oxygen as a terminal electron sink. The membrane-bound cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water in mitochondria and many bacteria. The energy released in this reaction is conserved by pumping protons across the mitochondrial or bacterial membrane, creating an electrochemical proton gradient that drives production of ATP. A crucial question is how the protons pumped by CcO are prevented from flowing backwards during the process. Here, we show by molecular dynamics simulations that the conserved glutamic acid 242 near the active site of CcO undergoes a protonation state-dependent conformational change, which provides a valve in the pumping mechanism. The valve ensures that at any point in time, the proton pathway across the membrane is effectively discontinuous, thereby preventing thermodynamically favorable proton back-leakage while maintaining an overall high efficiency of proton translocation. Suppression of proton leakage is particularly important in mitochondria under physiological conditions, where production of ATP takes place in the presence of a high electrochemical proton gradient.

Bubbles, gating, and anesthetics in ion channels

We suggest that bubbles are the bistable hydrophobic gates responsible for the on-off transitions of single channel currents. In this view, many types of channels gate by the same physical mechanism-dewetting by capillary evaporation-but different types of channels use different sensors to modulate hydrophobic properties of the channel wall and thereby trigger and control bubbles and gating. Spontaneous emptying of channels has been seen in many simulations. Because of the physics involved, such phase transitions are inherently sensitive, unstable threshold phenomena that are difficult to simulate reproducibly and thus convincingly. We present a thermodynamic analysis of a bubble gate using morphometric density functional theory of classical (not quantum) mechanics. Thermodynamic analysis of phase transitions is generally more reproducible and less sensitive to details than simulations. Anesthetic actions of inert gases-and their interactions with hydrostatic pressure (e.g., nitrogen narcosis)-can be easily understood by actions on bubbles. A general theory of gas anesthesia may involve bubbles in channels. Only experiments can show whether, or when, or which channels actually use bubbles as hydrophobic gates: direct observation of bubbles in channels is needed. Existing experiments show thin gas layers on hydrophobic surfaces in water and suggest that bubbles nearly exist in bulk water.

Molecular mechanism of proton translocation by cytochrome c oxidase

Cytochrome c oxidase (CcO) is a terminal protein of the respiratory chain in eukaryotes and some bacteria. It catalyzes most of the biologic oxygen consumption on earth done by aerobic organisms. During the catalytic reaction, CcO reduces dioxygen to water and uses the energy released in this process to maintain the electrochemical proton gradient by functioning as a redox-linked proton pump. Even though the structures of several terminal oxidases are known, they are not sufficient in themselves to explain the molecular mechanism of proton pumping. Thus, additional extensive studies of CcO by varieties of biophysical and biochemical approaches are involved to shed light on the mechanism of proton translocation. In this review, we summarize the current level of knowledge about CcO, including the latest model developed to explain the CcO proton-pumping mechanism.