Proton-transport mechanisms in cytochrome c oxidase revealed by studies of kinetic isotope effects

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

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

The proton pumping bo oxidase from Vitreoscilla

The cytochrome bo3 quinol oxidase from Vitreoscilla (vbo3) catalyses oxidation of ubiquinol and reduction of O2 to H2O. Data from earlier studies suggested that the free energy released in this reaction is used to pump sodium ions instead of protons across a membrane. Here, we have studied the functional properties of heterologously expressed vbo3 with a variety of methods. (i) Following oxygen consumption with a Clark-type electrode, we did not observe a measurable effect of Na+ on the oxidase activity of purified vbo3 solubilized in detergent or reconstituted in liposomes. (ii) Using fluorescent dyes, we find that vbo3 does not pump Na+ ions, but H+ across the membrane, and that H+-pumping is not influenced by the presence of Na+. (iii) Using an oxygen pulse method, it was found that 2 H+/e- are ejected from proteoliposomes, in agreement with the values found for the H+-pumping bo3 oxidase of Escherichia coli (ecbo3). This coincides with the interpretation that 1 H+/e- is pumped across the membrane and 1 H+/e- is released during quinol oxidation. (iv) When the electron transfer kinetics of vbo3 upon reaction with oxygen were followed in single turnover experiments, a similar sequence of reaction steps was observed as reported for the E. coli enzyme and none of these reactions was notably affected by the presence of Na+. Overall the data show that vbo3 is a proton pumping terminal oxidase, behaving similarly to the Escherichia coli bo3 quinol oxidase.

The regulation of oxidative phosphorylation pathway on Vibrio alginolyticus adhesion under adversities

Vibrio alginolyticus is one of the most important pathogens in mariculture and leading to heavy losses. After treatment with Cu²⁺, Pb²⁺, and low pH, the expression of oxidative phosphorylation pathway genes, including coxA, coxB, coxC, ccoN, ccoO, and ccoQ, was found commonly downregulated by RNA‐seq as well as quantitative real‐time PCR. RNAi significantly reduced the expression of coxA, coxB, coxC, ccoN, ccoO, and ccoQ in V. alginolyticus. Compared with the wild‐type strain, the adhesion abilities of RNAi strains of V. alginolyticus were significantly impaired, as well as their cytochrome C oxidase activity. ccoQ appeared to be more important in the regulation of bacterial adhesion in these target genes, while ccoO was relatively weak in the regulation of the adhesion. Meanwhile, the changes of temperature, salinity, pH, and starvation affected coxA, coxB, coxC, ccoN, ccoO, and ccoQ expression remarkably. These findings indicated that: the oxidative phosphorylation pathway is a critical regulator of adhesion in V. alginolyticus; coxA, coxB, coxC, ccoN, ccoO, and ccoQ regulate the bacterial adhesion in response to environmental changes such as temperature, salinity, pH, and starvation.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.

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

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

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

Origin of the Non-Arrhenius Behavior of the Rates of Enzymatic Reactions

The origin of the non-Arrhenius behavior of the rate constant for hydride transfer enzymatic reactions has been a puzzling problem since its initial observation. This effect has been used originally to support the idea that enzymes work by dynamical effects and more recently to suggest an entropy funnel model. Our analysis, however, has advanced the idea that the reason for the non-Arrhenius trend reflects the temperature dependence of the rearrangements of the protein polar groups in response to the change in the charge distribution of the reacting system during the transition from the ground state (GS) to the transition state (TS). Here we examine the validity of our early proposal by simulating the catalytic reaction of alcohol dehydrogenase (ADH) and determine the microscopic origin of the entropic and enthalpic contributions to the activation barrier. The corresponding analysis establishes the origin of the non-Arrhenius behaviors and quantifies our original suggestion that the classical effect is due to the entropic contributions of the environment. We also find that the quantum effects reflect in part the temperature dependence of the donor-acceptor distance.

Redox induced protonation of heme propionates in cytochrome c oxidase: Insights from surface enhanced resonance Raman spectroscopy and QM/MM calculations

Understanding the coupling between heme reduction and proton translocation in cytochrome c oxidase (CcO) is still an open problem. The propionic acids of heme a₃ have been proposed to act as a proton loading site (PLS) in the proton pumping pathway, yet this proposal could not be verified by experimental data so far. We have set up an experiment where the redox states of the two hemes in CcO can be controlled via external electrical potential. Surface enhanced resonance Raman (SERR) spectroscopy was applied to simultaneously monitor the redox state of the hemes and the protonation state of the heme propionates. Simulated spectra based on QM/MM calculations were used to assign the resonant enhanced CH₂ twisting modes of the propionates to the protonation state of the individual heme a and heme a₃ propionates respectively. The comparison between calculated and measured H₂OD₂O difference spectra allowed a sound band assignment. In the fully reduced enzyme at least three of the four heme propionates were found to be protonated whereas in the presence of a reduced heme a and an oxidized heme a₃ only protonation of one heme a₃ propionates was observed. Our data supports the postulated scenario where the heme a₃ propionates are involved in the proton pathway.

Multiscale simulations reveal key features of the proton-pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces oxygen to water and uses the released free energy to pump protons across the membrane. We have used multiscale reactive molecular dynamics simulations to explicitly characterize (with free-energy profiles and calculated rates) the internal proton transport events that enable proton pumping during first steps of oxidation of the fully reduced enzyme. Our results show that proton transport from amino acid residue E286 to both the pump loading site (PLS) and to the binuclear center (BNC) are thermodynamically driven by electron transfer from heme a to the BNC, but that the former (i.e., pumping) is kinetically favored whereas the latter (i.e., transfer of the chemical proton) is rate-limiting. The calculated rates agree with experimental measurements. The backflow of the pumped proton from the PLS to E286 and from E286 to the inside of the membrane is prevented by large free-energy barriers for the backflow reactions. Proton transport from E286 to the PLS through the hydrophobic cavity and from D132 to E286 through the D-channel are found to be strongly coupled to dynamical hydration changes in the corresponding pathways and, importantly, vice versa.

Bio-functionalized hybrid nanocomposite membranes for direct methanol fuel cells

Featured methanol-blocking characteristics of nanocomposite membrane.

Structural Changes and Proton Transfer in Cytochrome c Oxidase

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

Kinetic and vibrational isotope effects of proton transfer reactions in channelrhodopsin-2

Channelrhodopsins (ChRs) are light-gated cation channels. After blue-light excitation, the protein undergoes a photocycle with different intermediates. Here, we have recorded transient absorbance changes of ChR2 from Chlamydomonas reinhardtii in the visible and infrared regions with nanosecond time resolution, the latter being accomplished using tunable quantum cascade lasers. Because proton transfer reactions play a key role in channel gating, we determined vibrational as well as kinetic isotope effects (VIEs and KIEs) of carboxylic groups of various key aspartic and glutamic acid residues by monitoring their C=O stretching vibrations in H2O and in D2O. D156 exhibits a substantial KIE (>2) in its deprotonation and reprotonation, which substantiates its role as the internal proton donor to the retinal Schiff base. The unusual VIE of D156, upshifted from 1736 cm(-1) to 1738 cm(-1) in D2O, was scrutinized by studying the D156E variant. The C=O stretch of E156 shifted down by 8 cm(-1) in D2O, providing evidence for the accessibility of the carboxylic group. The C=O stretching band of E90 exhibits a VIE of 9 cm(-1) and a KIE of ∼2 for the de- and the reprotonation reactions during the lifetime of the late desensitized state. The KIE of 1 determined in the time range from 20 ns to 5 ms is incompatible with early deprotonation of E90.

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

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

Mutation of a single residue in the ba3 oxidase specifically impairs protonation of the pump site

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound protein complex that couples electron transfer to O2 to proton translocation across the membrane. To elucidate the mechanism of the redox-driven proton pumping, we investigated the kinetics of electron and proton transfer in a structural variant of the ba3 oxidase where a putative "pump site" was modified by replacement of Asp372 by Ile. In this structural variant, proton pumping was uncoupled from internal electron transfer and O2 reduction. The results from our studies show that proton uptake to the pump site (time constant ∼65 μs in the wild-type cytochrome c oxidase) was impaired in the Asp372Ile variant. Furthermore, a reaction step that in the wild-type cytochrome c oxidase is linked to simultaneous proton uptake and release with a time constant of ∼1.2 ms was slowed to ∼8.4 ms, and in Asp372Ile was only associated with proton uptake to the catalytic site. These data identify reaction steps that are associated with protonation and deprotonation of the pump site, and point to the area around Asp372 as the location of this site in the ba3 cytochrome c oxidase.

Site-specific protonation kinetics of acidic side chains in proteins determined by pH-dependent carboxyl (13)C NMR relaxation

Proton-transfer dynamics plays a critical role in many biochemical processes, such as proton pumping across membranes and enzyme catalysis. The large majority of enzymes utilize acid-base catalysis and proton-transfer mechanisms, where the rates of proton transfer can be rate limiting for the overall reaction. However, measurement of proton-exchange kinetics for individual side-chain carboxyl groups in proteins has been achieved in only a handful of cases, which typically have involved comparative analysis of mutant proteins in the context of reaction network modeling. Here we describe an approach to determine site-specific protonation and deprotonation rate constants (kon and koff, respectively) of carboxyl side chains, based on (13)C NMR relaxation measurements as a function of pH. We validated the method using an extensively studied model system, the B1 domain of protein G, for which we measured rate constants koff in the range (0.1-3) × 10(6) s(-1) and kon in the range (0.6-300) × 10(9) M(-1) s(-1), which correspond to acid-base equilibrium dissociation constants (Ka) in excellent agreement with previous results determined by chemical shift titrations. Our results further reveal a linear free-energy relationship between log kon and pKa, which provides information on the free-energy landscape of the protonation reaction, showing that the variability among residues in these parameters arises primarily from the extent of charge stabilization of the deprotonated state by the protein environment. We find that side-chain carboxyls with extreme values of koff or kon are involved in hydrogen bonding, thus providing a mechanistic explanation for the observed stabilization of the protonated or deprotonated state.

Redox changes accompanying inorganic carbon limitation in Synechocystis sp. PCC 6803

Inorganic carbon (Ci) is the major sink for photosynthetic reductant in organisms capable of oxygenic photosynthesis. In the absence of abundant Ci, the cyanobacterium Synechocystis sp. strain PCC6803 expresses a high affinity Ci acquisition system, the CO2-concentrating mechanisms (CCM), controlled by the transcriptional regulator CcmR and the metabolites NADP+ and α-ketoglutarate, which act as co-repressors of CcmR by modulating its DNA binding. The CCM thus responds to internal cellular redox changes during the transition from Ci-replete to Ci-limited conditions. However, the actual changes in the metabolic state of the NADPH/NADP+ system that occur during the transition to Ci-limited conditions remain ill-defined. Analysis of changes in the redox state of cells experiencing Ci limitation reveals systematic changes associated with physiological adjustments and a trend towards the quinone and NADP pools becoming highly reduced. A rapid and persistent increase in F0 was observed in cells reaching the Ci-limited state, as was the induction of photoprotective fluorescence quenching. Systematic changes in the fluorescence induction transients were also observed. As with Chl fluorescence, a transient reduction of the NADPH pool ('M' peak), is assigned to State 2→State 1 transition associated with increased electron flow to NADP+. This was followed by a characteristic decline, which was abolished by Ci limitation or inhibition of the Calvin-Benson-Bassham (CBB) cycle and is thus assigned to the activation of the CBB cycle. The results are consistent with the proposed regulation of the CCM and provide new information on the nature of the Chl and NADPH fluorescence induction curves.

Characterizing the proton loading site in cytochrome c oxidase

Cytochrome c oxidase (CcO) uses the energy released by reduction of O2 to H2O to drive eight charges from the high pH to low pH side of the membrane, increasing the electrochemical gradient. Four electrons and protons are used for chemistry, while four more protons are pumped. Proton pumping requires that residues on a pathway change proton affinity through the reaction cycle to load and then release protons. The protonation states of all residues in CcO are determined in MultiConformational Continuum Electrostatics simulations with the protonation and redox states of heme a, a3, Cu(B), Y288, and E286 used to define the catalytic cycle. One proton is found to be loaded and released from residues identified as the proton loading site (PLS) on the P-side of the protein in each of the four CcO redox states. Thus, the same proton pumping mechanism can be used each time CcO is reduced. Calculations with structures of Rhodobacter sphaeroides, Paracoccus denitrificans, and bovine CcO derived by crystallography and molecular dynamics show the PLS functions similarly in different CcO species. The PLS is a cluster rather than a single residue, as different structures show 1-4 residues load and release protons. However, the proton affinity of the heme a3 propionic acids primarily determines the number of protons loaded into the PLS; if their proton affinity is too low, less than one proton is loaded.

Current advances in research of cytochrome c oxidase

The function of cytochrome c oxidase as a biomolecular nanomachine that transforms energy of redox reaction into protonmotive force across a biological membrane has been subject of intense research, debate, and controversy. The structure of the enzyme has been solved for several organisms; however details of its molecular mechanism of proton pumping still remain elusive. Particularly, the identity of the proton pumping site, the key element of the mechanism, is still open to dispute. The pumping mechanism has been for a long time one of the key unsolved issues of bioenergetics and biochemistry, but with the accelerating progress in this field many important details and principles have emerged. Current advances in cytochrome oxidase research are reviewed here, along with a brief discussion of the most complete proton pumping mechanism proposed to date, and a molecular basis for control of its efficiency.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

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.

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

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

The membrane modulates internal proton transfer in cytochrome c oxidase

The functionality of membrane proteins is often modulated by the surrounding membrane. Here, we investigated the effect of membrane reconstitution of purified cytochrome c oxidase (CytcO) on the kinetics and thermodynamics of internal electron and proton-transfer reactions during O(2) reduction. Reconstitution of the detergent-solubilized enzyme in small unilamellar soybean phosphatidylcholine vesicles resulted in a lowering of the pK(a) in the pH dependence profile of the proton-uptake rate. This pK(a) change resulted in decreased proton-uptake rates in the pH range of ~6.5-9.5, which is explained in terms of lowering of the pK(a) of an internal proton donor within CytcO. At pH 7.5, the rate decreased to the same extent when vesicles were prepared from the pure zwitterionic lipid 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or the anionic lipid 1,2-dioleoyl-sn-glycero-3-phospho(1-rac-glycerol) (DOPG). In addition, a small change in the internal Cu(A)-heme a electron equilibrium constant was observed. This effect was lipid-dependent and explained in terms of a lower electrostatic potential within the membrane-spanning part of the protein with the anionic DOPG lipids than with the zwitterionic DOPC lipids. In conclusion, the data show that the membrane significantly modulates internal charge-transfer reactions and thereby the function of the membrane-bound enzyme.

Insights into the mechanism of proton transport in cytochrome c oxidase

Cytochrome c oxidase (CcO), known as complex IV of the electron transport chain, plays several important roles in aerobic cellular respiration. Electrons transferred from cytochrome c to CcO's catalytic site reduce molecular oxygen and produce a water molecule. These electron transfers also drive active proton pumping from the matrix (N-side) to intermembrane region (P-side) in mitochondria; the resultant proton gradient activates ATP synthase to produce ATP from ADP. Although the existence of the coupling between the electron transfer and the proton transport (PT) is established experimentally, its mechanism is not yet fully understood at the molecular level. In this work, it is shown why the reduction of heme a is essential for proton pumping. This is demonstrated via novel reactive molecular dynamics (MD) simulations that can describe the Grotthuss shuttling associated with the PT as well as the dynamic delocalization of the excess proton electronic charge defect. Moreover, the "valve" role of the Glu242 residue (bovine CcO notation) and the gate role of d-propionate of heme a(3) (PRDa3) in the explicit PT are explicitly demonstrated for the first time. These results provide conclusive evidence for the CcO proton transporting mechanism inferred from experiments, while deepening the molecular level understanding of the CcO proton switch.