Titration behavior of residues at the entrance of the D-pathway of cytochrome c oxidase from paracoccus denitrificans investigated by continuum electrostatic calculations

Role of aspartate 132 at the orifice of a proton pathway in cytochrome c oxidase

Proton transfer across biological membranes underpins central processes in biological systems, such as energy conservation and transport of ions and molecules. In the membrane proteins involved in these processes, proton transfer takes place through specific pathways connecting the two sides of the membrane via control elements within the protein. It is commonly believed that acidic residues are required near the orifice of such proton pathways to facilitate proton uptake. In cytochrome c oxidase, one such pathway starts near a conserved Asp-132 residue. Results from earlier studies have shown that replacement of Asp-132 by, e.g., Asn, slows proton uptake by a factor of ∼5,000. Here, we show that proton uptake at full speed (∼10(4) s(-1)) can be restored in the Asp-132-Asn oxidase upon introduction of a second structural modification further inside the pathway (Asn-139-Thr) without compensating for the loss of the negative charge. This proton-uptake rate was insensitive to Zn(2+) addition, which in the wild-type cytochrome c oxidase slows the reaction, indicating that Asp-132 is required for Zn(2+) binding. Furthermore, in the absence of Asp-132 and with Thr at position 139, at high pH (>9), proton uptake was significantly accelerated. Thus, the data indicate that Asp-132 is not strictly required for maintaining rapid proton uptake. Furthermore, despite the rapid proton uptake in the Asn-139-Thr/Asp-132-Asn mutant cytochrome c oxidase, proton pumping was impaired, which indicates that the segment around these residues is functionally linked to pumping.

Proton uptake and pKa changes in the uncoupled Asn139Cys variant of cytochrome c oxidase

Cytochrome c oxidase (CytcO) is a membrane-bound enzyme that links electron transfer from cytochrome c to O(2) to proton pumping across the membrane. Protons are transferred through specific pathways that connect the protein surface with the catalytic site as well as the proton input with the proton output sides. Results from earlier studies have shown that one site within the so-called D proton pathway, Asn139, located ~10 Å from the protein surface, is particularly sensitive to mutations that uncouple the O(2) reduction reaction from the proton pumping activity. For example, none of the Asn139Asp (charged) or Asn139Thr (neutral) mutant CytcOs pump protons, although the proton-uptake rates are unaffected. Here, we have investigated the Asn139Cys and Asn139Cys/Asp132Asn mutant CytcOs. In contrast to other structural variants investigated to date, the Cys side chain may be either neutral or negatively charged in the experimentally accessible pH range. The data show that the Asn139Cys and Asn139Asp mutations result in the same changes of the kinetic and thermodynamic parameters associated with the proton transfer. The similarity is not due to introduction of charge at position 139, but rather introduction of a protonatable group that modulates the proton connectivity around this position. These results illuminate the mechanism by which CytcO couples electron transfer to proton pumping.

Electron transfer pathways in cytochrome c oxidase

Mixed quantum mechanical/molecular mechanics calculations were used to explore the electron pathway of the terminal electron transfer enzyme, cytochrome c oxidase. This enzyme catalyzes the reduction of molecular oxygen to water in a multiple step process. Density functional calculations on the three redox centers allowed for the characterization of the electron transfer mechanism, following the sequence Cu(A)→heme a→heme a(3). This process is largely affected by the presence of positive charges, confirming the possibility of a proton coupled electron transfer. An extensive mapping of all residues involved in the electron transfer, between the Cu(A) center (donor) and the O(2) reduction site heme a(3)-Cu(B) (receptor), was obtained by selectively activating/deactivating different quantum regions. The method employed, called QM/MM e-pathway, allowed the identification of key residues along the possible electron transfer paths, consistent with experimental data. In particular, the role of arginines 481 and 482 appears crucial in the Cu(A)→heme a and in the heme a→heme a(3) electron transfer processes. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Decoupling mutations in the D-channel of the aa(3)-type cytochrome c oxidase from Rhodobacter sphaeroides suggest that a continuous hydrogen-bonded chain of waters is essential for proton pumping

The aa(3)-type cytochrome c oxidase from Rhodobacter sphaeroides utilizes two proton-input channels to provide all the protons for chemistry (water formation) and proton pumping. The D-channel is responsible for the uptake of all pumped protons, four protons per O(2). Several substitutions of either N139 or N207, near the entrance of the D-channel, were previously reported to decouple the proton pump from oxidase activity. In this work, the characteristics of additional mutations in this region of the protein (N139, N207, N121, and S142) are determined to elucidate the mechanism of decoupling. With the exception of the substitution of a large, hydrophobic residue (N139L), all the mutations of N139 resulted in an enzyme with high oxidase activity but with a severely diminished proton pumping stoichiometry. Whereas N207D was previously shown to be decoupled, N207A and N207T exhibit nearly wild-type behavior. The new data display a pattern. Small, nonionizable substitutions of N139 or N121 result in decoupling of the proton pump but maintain high turnover rates. These residues are directly hydrogen bonded to two water molecules (Water6574 and Water6584) that are part of the single-file chain of water molecules within the D-channel leading to E286 at the top of the channel. The data suggest that the integrity of this water chain within the D-channel is critical for rapid proton transfer. The mechanism of decoupling is most likely due to the slowing of the rate of proton delivery below a threshold that is required for protonation of the putative proton loading site. Protons delivered outside this time window are delivered to the active site where they are consumed in the formation of water. The rate of proton delivery required to protonate the pump site must be significantly faster than the rate of delivery of protons to the catalytic site. For this reason, mutations can result in decoupling of the proton pump without slowing the catalytic turnover by the enzyme.

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.

A pathogenic mutation in cytochrome c oxidase results in impaired proton pumping while retaining O(2)-reduction activity

In this work we have investigated the effect of a pathogenic mitochondrial DNA mutation found in human colon cells, at a functional-molecular level. The mutation results in the amino-acid substitution Tyr19His in subunit I of the human CytcO and it is associated with respiratory deficiency. It was introduced into Rhodobacter sphaeroides, which carries a cytochrome c oxidase (cytochrome aa(3)) that serves as a model of the mitochondrial counterpart. The residue is situated in the middle of a pathway that is used to transfer substrate protons as well as protons that are pumped across the membrane. The Tyr33His (equivalent residue in the bacterial CytcO) structural variant of the enzyme was purified and its function was investigated. The results show that in the structurally altered CytcO the activity decreased due to slowed proton transfer; proton transfer from an internal proton donor, the highly-conserved Glu286, to the catalytic site was slowed by a factor of approximately 5, while reprotonation of the Glu from solution was slowed by a factor of approximately 40. In addition, in the structural variant proton pumping was completely impaired. These results are explained in terms of introduction of a barrier for proton transfer through the D pathway and changes in the coordination of water molecules surrounding the Glu286 residue. The study offers an explanation, at the molecular level, to the link between a specific amino-acid substitution and a pathogenic phenotype identified in human colon cells.

Association of cytochrome c with membrane-bound cytochrome c oxidase proceeds parallel to the membrane rather than in bulk solution

Electron transfer between the water-soluble cytochrome c and the integral membrane protein cytochrome c oxidase (COX) is the terminal reaction in the respiratory chain. The first step in this reaction is the diffusional association of cytochrome c toward COX, and it is still not completely clear whether cytochrome c diffuses in the bulk solution while encountering COX, or whether it prefers to diffuse laterally on the membrane surface. This is a rather crucial question, since in the latter case the association would be strongly dependent on the lipid composition and the presence of additional membrane proteins. We applied Brownian dynamics simulations to investigate the effect of an atomistically modeled dipalmitoyl phosphatidylcholine membrane on the association behavior of cytochrome c toward COX from Paracoccus denitrificans. We studied the negatively charged, physiological electron-transfer partner of COX, cytochrome c(552), and the positively charged horse-heart cytochrome c. As expected, both cytochrome c species prefer diffusion in bulk solution while associating toward COX embedded in a membrane, where the partial charges of the lipids were switched off, and the corresponding optimal association pathways largely overlap with the association toward fully solvated COX. Remarkably, after switching on the lipid partial charges, both cytochrome c species were strongly attracted by the inhomogeneous charge distribution caused by the zwitterionic lipid headgroups. This effect is particularly enhanced for horse-heart cytochrome c and is stronger at lower ionic strength. We therefore conclude that in the presence of a polar or even a charged membrane, cytochrome c diffuses laterally rather than in three dimensions.

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.

High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: new insights into the active site and the proton transfer pathways

The structure of the two-subunit cytochrome c oxidase from Paracoccus denitrificans has been refined using X-ray cryodata to 2.25 A resolution in order to gain further insights into its mechanism of action. The refined structural model shows a number of new features including many additional solvent and detergent molecules. The electron density bridging the heme a(3) iron and Cu(B) of the active site is fitted best by a peroxo-group or a chloride ion. Two waters or OH(-) groups do not fit, one water (or OH(-)) does not provide sufficient electron density. The analysis of crystals of cytochrome c oxidase isolated in the presence of bromide instead of chloride appears to exclude chloride as the bridging ligand. In the D-pathway a hydrogen bonded chain of six water molecules connects Asn131 and Glu278, but the access for protons to this water chain is blocked by Asn113, Asn131 and Asn199. The K-pathway contains two firmly bound water molecules, an additional water chain seems to form its entrance. Above the hemes a cluster of 13 water molecules is observed which potentially form multiple exit pathways for pumped protons. The hydrogen bond pattern excludes that the Cu(B) ligand His326 is present in the imidazolate form.

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.

Ultrafast dynamics and anionic active states of the flavin cofactor in cryptochrome and photolyase

We report here our systematic studies of the dynamics of four redox states of the flavin cofactor in both photolyases and insect type 1 cryptochromes. With femtosecond resolution, we observed ultrafast photoreduction of oxidized state flavin adenine dinucleotide (FAD) in subpicosecond and of neutral radical semiquinone (FADH(*)) in tens of picoseconds through intraprotein electron transfer mainly with a neighboring conserved tryptophan triad. Such ultrafast dynamics make these forms of flavin unlikely to be the functional states of the photolyase/cryptochrome family. In contrast, we find that upon excitation the anionic semiquinone (FAD(*-)) and hydroquinone (FADH(-)) have longer lifetimes that are compatible with high-efficiency intermolecular electron transfer reactions. In photolyases, the excited active state (FADH(-)*) has a long (nanosecond) lifetime optimal for DNA-repair function. In insect type 1 cryptochromes known to be blue-light photoreceptors the excited active form (FAD(*-)*) has complex deactivation dynamics on the time scale from a few to hundreds of picoseconds, which is believed to occur through conical intersection(s) with a flexible bending motion to modulate the functional channel. These unique properties of anionic flavins suggest a universal mechanism of electron transfer for the initial functional steps of the photolyase/cryptochrome blue-light photoreceptor family.

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.

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.

Impaired proton pumping in cytochrome c oxidase upon structural alteration of the D pathway

Cytochrome c oxidase is a membrane-bound enzyme, which catalyses the one-electron oxidation of four molecules of cytochrome c and the four-electron reduction of O(2) to water. Electron transfer through the enzyme is coupled to proton pumping across the membrane. Protons that are pumped as well as those that are used for O(2) reduction are transferred though a specific intraprotein (D) pathway. Results from earlier studies have shown that replacement of residue Asn139 by an Asp, at the beginning of the D pathway, results in blocking proton pumping without slowing uptake of substrate protons used for O(2) reduction. Furthermore, introduction of the acidic residue results in an increase of the apparent pK(a) of E286, an internal proton donor to the catalytic site, from 9.4 to ~11. In this study we have investigated intramolecular electron and proton transfer in a mutant cytochrome c oxidase in which a neutral residue, Thr, was introduced at the 139 site. The mutation results in uncoupling of proton pumping from O(2) reduction, but a decrease in the apparent pK(a) of E286 from 9.4 to 7.6. The data provide insights into the mechanism by which cytochrome c oxidase pumps protons and the structural elements involved in this process.

Electrostatic control of proton pumping in cytochrome c oxidase

As part of the mitochondrial respiratory chain, cytochrome c oxidase utilizes the energy produced by the reduction of O2 to water to fuel vectorial proton transport. The mechanism coupling proton pumping to redox chemistry is unknown. Recent advances have provided evidence that each of the four observable transitions in the complex catalytic cycle consists of a similar sequence of events. However, the physico-chemical basis underlying this recurring sequence has not been identified. We identify this recurring pattern based on a comprehensive model of the catalytic cycle derived from the analysis of oxygen chemistry and available experimental evidence. The catalytic cycle involves the periodic repetition of a sequence of three states differing in the spatial distribution of charge in the active site: [0|1], [1|0], and [1|1], where the total charge of heme a and the binuclear center appears on the left and on the right, respectively. This sequence recurs four times per turnover despite differences in the redox chemistry. This model leads to a simple, robust, and reproducible sequence of electron and proton transfer steps and rationalizes the pumping mechanism in terms of electrostatic coupling of proton translocation to redox chemistry. Continuum electrostatic calculations support the proposed mechanism and suggest an electrostatic origin for the decoupled and inactive phenotypes of ionic mutants in the principal proton-uptake pathway.

The proton pumping pathway of bovine heart cytochrome c oxidase

X-ray structures of bovine heart cytochrome c oxidase have suggested that the enzyme, which reduces O(2) in a process coupled with a proton pumping process, contains a proton pumping pathway (H-pathway) composed of a hydrogen bond network and a water channel located in tandem across the enzyme. The hydrogen bond network includes the peptide bond between Tyr-440 and Ser-441, which could facilitate unidirectional proton transfer. Replacement of a possible proton-ejecting aspartate (Asp-51) at one end of the H-pathway with asparagine, using a stable bovine gene expression system, abolishes the proton pumping activity without influencing the O(2) reduction function. Blockage of either the water channel by a double mutation (Val386Leu and Met390Trp) or proton transfer through the peptide by a Ser441Pro mutation was found to abolish the proton pumping activity without impairment of the O(2) reduction activity. These results significantly strengthen the proposal that H-pathway is involved in proton pumping.

Kinetic models of redox-coupled proton pumping

Cytochrome c oxidase, the terminal enzyme of the respiratory chain, pumps protons across the inner mitochondrial membrane against an opposing electrochemical gradient by reducing oxygen to water. To explore the fundamental mechanisms of such redox-coupled proton pumps, we develop kinetic models at the single-molecule level consistent with basic physical principles. We demonstrate that pumping against potentials >150 mV can be achieved purely through electrostatic couplings, given an asymmetric arrangement of charge centers; however, nonlinear gates are essential for highly efficient real enzymes. The fundamental requirements for proton pumping identified here highlight a possible evolutionary origin of cytochrome c oxidase pumping. The general design principles are relevant also for other molecular machines and suggest future applications in biology-inspired fuel cells.

Replacing Asn207 by aspartate at the neck of the D channel in the aa3-type cytochrome c oxidase from Rhodobacter sphaeroides results in decoupling the proton pump

Cytochrome oxidase catalyzes the reduction of O2 to water and conserves the considerable free energy available from this reaction in the form of a proton motive force. For each electron, one proton is electrogenically pumped across the membrane. Of particular interest is the mechanism by which the proton pump operates. Previous studies of the oxidase from Rhodobacter sphaeroides have shown that all of the pumped protons enter the enzyme through the D channel and that a point mutant, N139D, in the D channel completely eliminates proton pumping without reducing oxidase activity. N139 is one of three asparagines near the entrance of the D channel, where there is a narrowing or neck, through which a single file of water molecules pass. In the current work, it is shown that replacement of a second asparagine in this region by an asparate, N207D, also decouples the proton pump without altering the oxidase activity of the enzyme. Previous studies demonstrated that the N139D mutant results in an increase in the apparent pKa of E286, a functionally critical residue that is located 20 A away from N139 at the opposite end of the D channel. In the current work, it is shown that the N207 mutation also increases the apparent pKa of E286. This finding reinforces the proposal that the elimination of proton pumping is the result of an increase of the apparent proton affinity of E286, which, in turn, prevents the timely proton transfer to a proton accepter group within the exit channel of the proton pump.

Mutations which decouple the proton pump of the cytochrome c oxidase from Rhodobacter sphaeroides perturb the environment of glutamate 286

Mutants that decouple the proton pump of cytochrome c oxidase from Rhodobacter sphaeroides are postulated to do so by increasing the pK(a) of glutamate 286, which is 20 Angstrom away. The possibility that a conformational change near E286 is induced by the decoupling mutations (N139D and N207D) was investigated by FTIR difference spectroscopy. In both decoupled mutants, the reduced-minus-oxidized FTIR difference spectra show a shift of 2 cm(-1) to lower frequency of the band resulting from the absorbance of E286 in the oxidized enzyme. The decoupling mutants may influence E286 by altering the chain of water molecules which runs from the site of the mutations to E286.

Free energy profiles for H+ conduction in the D-pathway of Cytochrome c Oxidase: a study of the wild type and N98D mutant enzymes

The molecular mechanism for proton conduction in the D-pathway of Cytochrome c Oxidase (CcO) is investigated through the free energy profile, i.e., potential of mean force (PMF) calculations of both the native enzyme and the N98D mutant. The multistate empirical valence bond (MS-EVB) model was applied to simulate the interaction of an excess proton with the channel environment. In the study of the wild type enzyme, the PMF reveals the previously proposed proton trap inside the channel; it also shows a high free energy barrier against the passage of proton at the entry of the channel, where two conserved asparagines (ASN80/98) may be essential for the gating of proton uptake. We also present data from an investigation of the N98D mutant, which has been previously shown to completely eliminate proton pumping but significantly enhance the oxidase activity in Rhodobacter sphaeroides. These results suggest that mutating Asn98 to negatively charged aspartate will create an unfavorable energy barrier sufficiently high to prevent the overall proton uptake through the D-pathway, whereas with a protonated aspartic acid the proton conduction was found to be accelerated. Plausible explanations for the origin of the uncoupling of proton pumping from the oxidase activity will be discussed.

Monte Carlo simulations of proton pumps: on the working principles of the biological valve that controls proton pumping in cytochrome c oxidase

Gaining a detailed understanding of the proton-pumping process in cytochrome c oxidase (COX) is one of the challenges of modern biophysics. Recent mutation experiments have highlighted this challenge by showing that a single mutation (the N139D mutation) blocks the overall pumping while continuing to channel protons to the binuclear center without inhibiting the oxidase activity. Rationalizing this result has been a major problem because the mutation is quite far from E286, which is believed to serve as the branching point for the proton transport in the pumping process. In the absence of a reasonable explanation for this important observation, we have developed a Monte Carlo simulation method that can convert mutation and structural information to pathways for proton translocation and simulate the pumping process in COX on a millisecond and even subsecond time scale. This tool allows us to reproduce and propose a possible explanation to the effect of the N139D mutation and to offer a consistent model for the origin of the "valve effect" in COX, which is crucial for maintaining uphill proton pumping. Furthermore, obtaining the first structure-based simulation of proton pumping in COX, or in any other protein, indicates that our approach should provide a powerful tool for verification of mechanistic hypotheses about the action of proton transport proteins.