Exploring the proton pump mechanism of cytochrome c oxidase in real time

Mitochondrial dysfunction is a hallmark of woody breast myopathy in broiler chickens

The woody breast (WB) myopathy poses significant economic and welfare concerns to the poultry industry, however, there is no effective strategy to mitigate this pathology due to its unknown etiology. After showing previously that hypoxia is a key factor in WB progression, we used here various techniques demonstrating dysregulated mitochondria (morphology, biogenesis, tethering, function, and bioenergetics) in WB-affected muscles and in hypoxic myoblasts compared to healthy tissues and normoxic cells, respectively. The increased levels of calcium (Ca2+) in both WB-affected tissues and hypoxic myoblasts suggested that mitochondrial Ca2+ overload is likely a leading cause for mitochondrial dysfunction that merits further in-depth investigation. These findings are the first, to the best of our knowledge, to provide fundamental insights into the underlying molecular mechanisms of WB and open new vistas for understanding the interplay between calcium, mitochondrial (dys)function, and avian muscle health for subsequent development of effective preventative/corrective strategies.

Long-range charge transfer mechanism of the III2IV2 mycobacterial supercomplex

Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III2IV2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria.

Effect of VBIT-4 on the functional activity of isolated mitochondria and cell viability

VBIT-4 is a new inhibitor of the oligomerization of VDAC proteins of the outer mitochondrial membrane preventing the development of oxidative stress, mitochondrial dysfunction, and cell death in various pathologies. However, as a VDAC inhibitor, VBIT-4 may itself cause mitochondrial dysfunction in healthy cells. The article examines the effect of VBIT-4 on the functional activity of rat liver mitochondria and cell cultures. We have demonstrated that high concentrations of VBIT-4 (15-30 μM) suppressed mitochondrial respiration in state 3 and 3UDNP driven by substrates of complex I and II. VBIT-4 induced depolarization of organelles fueled by substrates of complex I but not complex II of the respiratory chain. VBIT-4 has been found to inhibit the activity of complexes I, III, and IV of the respiratory chain. Molecular docking demonstrated that VBIT-4 interacts with the rotenone-binding site in complex I with similar affinity. 15-30 μM VBIT-4 caused an increase in H2O2 production in mitochondria, decreased the Ca2+ retention capacity, but increased the time of Ca2+-dependent mitochondrial swelling. We have found that the incubation of breast adenocarcinoma (MCF-7) with 30 μM VBIT-4 for 48 h led to the decrease of the mitochondrial membrane potential, an increase in ROS production and death of MCF-7 cells. The mechanism of action of VBIT-4 on mitochondria and cells is discussed.

Mitochondria at the Nanoscale: Physics Meets Biology-What Does It Mean for Medicine?

Mitochondria are commonly perceived as "cellular power plants". Intriguingly, power conversion is not their only function. In the first part of this paper, we review the role of mitochondria in the evolution of eukaryotic organisms and in the regulation of the human body, specifically focusing on cancer and autism in relation to mitochondrial dysfunction. In the second part, we overview our previous works, revealing the physical principles of operation for proton-pumping complexes in the inner mitochondrial membrane. Our proposed simple models reveal the physical mechanisms of energy exchange. They can be further expanded to answer open questions about mitochondrial functions and the medical treatment of diseases associated with mitochondrial disorders.

Coupled chemical reactions: Effects of electric field, diffusion, and boundary control

Chemical reactions involve the movement of charges, and this paper presents a mathematical model for describing chemical reactions in electrolytes. The model is developed using an energy variational method that aligns with classical thermodynamics principles. It encompasses both electrostatics and chemical reactions within consistently defined energetic and dissipative functionals. Furthermore, the energy variation method is extended to account for open systems that involve the input and output of charge and mass. Such open systems have the capability to convert one form of input energy into another form of output energy. In particular, a two-domain model is developed to study a reaction system with self-regulation and internal switching, which plays a vital role in the electron transport chain of mitochondria responsible for ATP generation-a crucial process for sustaining life. Simulations are conducted to explore the influence of electric potential on reaction rates and switching dynamics within the two-domain system. It shows that the electric potential inhibits the oxidation reaction while accelerating the reduction reaction.

The catalytic reaction of cytochrome c oxidase probed by in situ gas titrations and FTIR difference spectroscopy

Cytochrome c oxidase (CcO) is a transmembrane heme‑copper metalloenzyme that catalyzes the reduction of O2 to H2O at the reducing end of the respiratory electron transport chain. To understand this reaction, we followed the conversion of CcO from Rhodobacter sphaeroides between several active-ready and carbon monoxide-inhibited states via attenuated total reflection Fourier-transform infrared (ATR FTIR) difference spectroscopy. Utilizing a novel gas titration setup, we prepared the mixed-valence, CO-inhibited R2CO state as well as the fully-reduced R4 and R4CO states and induced the "active ready" oxidized state OH. These experiments are performed in the dark yielding FTIR difference spectra exclusively triggered by exposure to O2, the natural substrate of CcO. Our data demonstrate that the presence of CO at heme a3 does not impair the catalytic oxidation of CcO when the cycle starts from the fully-reduced states. Interestingly, when starting from the R2CO state, the release of the CO ligand upon purging with inert gas yield a product that is indistinguishable from photolysis-induced states. The observed changes at heme a3 in the catalytic binuclear center (BNC) result from the loss of CO and are unrelated to electronic excitation upon illumination. Based on our experiments, we re-evaluate the assignment of marker bands that appear in time-resolved photolysis and perfusion-induced experiments on CcO.

Crystallographic cyanide-probing for cytochrome c oxidase reveals structural bases suggesting that a putative proton transfer H-pathway pumps protons

Cytochrome c oxidase (CcO) reduces O2 in the O2-reduction site by sequential four-electron donations through the low-potential metal sites (CuA and Fea). Redox-coupled X-ray crystal structural changes have been identified at five distinct sites including Asp51, Arg438, Glu198, the hydroxyfarnesyl ethyl group of heme a, and Ser382, respectively. These sites interact with the putative proton-pumping H-pathway. However, the metal sites responsible for each structural change have not been identified, since these changes were detected as structural differences between the fully reduced and fully oxidized CcOs. Thus, the roles of these structural changes in the CcO function are yet to be revealed. X-ray crystal structures of cyanide-bound CcOs under various oxidation states showed that the O2-reduction site controlled only the Ser382-including site, while the low-potential metal sites induced the other changes. This finding indicates that these low-potential site-inducible structural changes are triggered by sequential electron-extraction from the low-potential sites by the O2-reduction site and that each structural change is insensitive to the oxidation and ligand-binding states of the O2-reduction site. Because the proton/electron coupling efficiency is constant (1:1), regardless of the reaction progress in the O2-reduction site, the structural changes induced by the low-potential sites are assignable to those critically involved in the proton pumping, suggesting that the H-pathway, facilitating these low-potential site-inducible structural changes, pumps protons. Furthermore, a cyanide-bound CcO structure suggests that a hypoxia-inducible activator, Higd1a, activates the O2-reduction site without influencing the electron transfer mechanism through the low-potential sites, kinetically confirming that the low-potential sites facilitate proton pump.

Examination of 'high-energy' metastable state of the oxidized (OH) bovine cytochrome c oxidase: Proton uptake and reaction with H2O2

Reoxidized cytochrome c oxidase appears to be in a 'high-energy' metastable state (OH) in which part of the energy released in the redox reactions is stored. The OH is supposed to relax to the resting 'as purified' oxidized state (O) in a time exceeding 200 ms. The catalytic heme a3-CuB center of these two forms should differ in a protonation and ligation state and the transition of OH-to-O is suggested to be associated with a proton transfer into this center. Employing a stopped-flow and UV-Vis absorption spectroscopy we investigated a proton uptake during the predicted relaxation of OH. It is shown, using a pH indicator phenol red, that from the time when the oxidation of the fully reduced CcO is completed (∼25 ms) up to ∼10 min, there is no uptake of a proton from the external medium (pH 7.8). Moreover, interactions of the assumed OH, generated 100 ms after oxidation of the fully reduced CcO, and the O with H2O2 (1 mM), result in the formation of two ferryl intermediates of the catalytic center, P and F, with very similar kinetics and the amounts of the formed ferryl states in both cases. These results implicate that the relaxation time of the catalytic center during the OH-to-O transition is either shorter than 100 ms or there is no difference in the structure of heme a3-CuB center of these two forms.

Investigation of the Mechanism of Membrane Potential Generation by Heme-Copper Respiratory Oxidases in a Real Time Mode

Heme-copper respiratory oxidases are highly efficient molecular machines. These membrane enzymes catalyze the final step of cellular respiration in eukaryotes and many prokaryotes: the transfer of electrons from cytochromes or quinols to molecular oxygen and oxygen reduction to water. The free energy released in this redox reaction is converted by heme-copper respiratory oxidases into the transmembrane gradient of the electrochemical potential of hydrogen ions H+). Heme-copper respiratory oxidases have a unique mechanism for generating H+, namely, a redox-coupled proton pump. A combination of direct electrometric method for measuring the kinetics of membrane potential generation with the methods of prestationary kinetics and site-directed mutagenesis in the studies of heme-copper oxidases allows to obtain a unique information on the translocation of protons inside the proteins in real time. The review summarizes the data of studies employing time-resolved electrometry to decipher the mechanisms of functioning of these important bioenergetic enzymes.

Protonation-State Dependence of Hydration and Interactions in the Two Proton-Conducting Channels of Cytochrome c Oxidase

Cytochrome c Oxidase (CcO), a membrane protein of the respiratory chain, pumps protons against an electrochemical gradient by using the energy of oxygen reduction to water. The ("chemical") protons required for this reaction and those pumped are taken up via two distinct channels, named D-channel and K-channel, in a step-wise and highly regulated fashion. In the reductive phase of the catalytic cycle, both channels transport protons so that the pumped proton passes the D-channel before the "chemical" proton has crossed the K-channel. By performing molecular dynamics simulations of CcO in the O→E redox state (after the arrival of the first reducing electron) with various combinations of protonation states of the D- and K-channels, we analysed the effect of protonation on the two channels. In agreement with previous work, the amount of water observed in the D-channel was significantly higher when the terminal residue E286 was not (yet) protonated than when the proton arrived at this end of the D-channel and E286 was neutral. Since a sufficient number of water molecules in the channel is necessary for proton transport, this can be understood as E286 facilitating its own protonation. K-channel hydration shows an even higher dependence on the location of the excess proton in the K-channel. Also in agreement with previous work, the K-channel exhibits a very low hydration level that likely hinders proton transfer when the excess proton is located in the lower part of the K-channel, that is, on the N-side of S365. Once the proton has passed S365 (towards the reaction site, the bi-nuclear centre (BNC)), the amount of water in the K-channel provides hydrogen-bond connectivity that renders proton transfer up to Y288 at the BNC feasible. No significant direct effect of the protonation state of one channel on the hydration level, hydrogen-bond connectivity, or interactions between protein residues in the other channel could be observed, rendering proton conductivity in the two channels independent of each other. Regulation of the order of proton uptake and proton passage in the two channels such that the "chemical" proton leaves its channel last must, therefore, be achieved by other means of communication, such as the location of the reducing electron.

Recent progress in experimental studies on the catalytic mechanism of cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces molecular oxygen (O₂) to water, coupled with a proton pump from the N-side to the P-side, by receiving four electrons sequentially from the P-side to the O₂-reduction site-including Fe_a3 and Cu_B-via the two low potential metal sites; Cu_A and Fe_a. The catalytic cycle includes six intermediates as follows, R (Fe_a3 ²⁺, Cu_B ¹⁺, Tyr244OH), A (Fe_a3 ²⁺-O₂, Cu_B ¹⁺, Tyr244OH), P_m (Fe_a3 ⁴⁺ = O^2-, Cu_B ²⁺-OH^-, Tyr244O•), F (Fe_a3 ⁴⁺ = O^2-, Cu_B ²⁺-OH^-, Tyr244OH), O (Fe_a3 ³⁺-OH^-, Cu_B ²⁺-OH^-, Tyr244OH), and E (Fe_a3 ³⁺-OH^-, Cu_B ¹⁺-H₂O, Tyr244OH). CcO has three proton conducting pathways, D, K, and H. The D and K pathways connect the N-side surface with the O₂-reduction site, while the H-pathway is located across the protein from the N-side to the P-side. The proton pump is driven by electrostatic interactions between the protons to be pumped and the net positive charges created during the O₂ reduction. Two different proton pump proposals, each including either the D-pathway or H-pathway as the proton pumping site, were proposed approximately 30 years ago and continue to be under serious debate. In our view, the progress in understanding the reaction mechanism of CcO has been critically rate-limited by the resolution of its X-ray crystallographic structure. The improvement of the resolutions of the oxidized/reduced bovine CcO up to 1.5/1.6 Å resolution in 2016 provided a breakthrough in the understanding of the reaction mechanism of CcO. In this review, experimental studies on the reaction mechanism of CcO before the appearance of the 1.5/1.6 Å resolution X-ray structures are summarized as a background description. Following the summary, we will review the recent (since 2016) experimental findings which have significantly improved our understanding of the reaction mechanism of CcO including: 1) redox coupled structural changes of bovine CcO; 2) X-ray structures of all six intermediates; 3) spectroscopic findings on the intermediate species including the Tyr244 radical in the P_m form, a peroxide-bound form between the A and Pm forms, and F_r, a one-electron reduced F-form; 4) time resolved X-ray structural changes during the photolysis of CO-bound fully reduced CcO using XFEL; 5) a simulation analysis for the Pm→Pr→F transition.

Interplay of Hydration and Protonation Dynamics in the K-Channel of Cytochrome c Oxidase

Cytochrome c oxidase is a membrane protein of the respiratory chain that consumes protons and molecular oxygen to produce water and uses the resulting energy to pump protons across the membrane. Our molecular dynamics simulations with an excess proton located at different positions in one of the proton-conducting channels, the K-channel, show a clear dependence of the number of water molecules inside the channel on the proton position. A higher hydration level facilitates the formation of hydrogen-bonded chains along which proton transfer can occur. However, a sufficiently high hydration level for such proton transport is observed only when the excess proton is located above S365, i.e., the lower third of the channel. From the channel entrance up to this point, proton transport is via water molecules as proton carriers. These hydronium ions move with their surrounding water molecules, up to K362, filling and widening the channel. The conformation of K362 depends on its own protonation state and on the hydration level, suggesting its role to be proton transport from a hydronium ion at the height of K362 to the upper part of the channel via a conformational change. The protonation-dependent conformational dynamics of E101 at the bottom of the channel renders proton transfer via E101 unlikely. Instead, its role is rather that of an amplifier of H96’s proton affinity, suggesting H96 as the initial proton acceptor.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O2 reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a3 D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Serial crystallography captures dynamic control of sequential electron and proton transfer events in a flavoenzyme

Flavin coenzymes are universally found in biological redox reactions. DNA photolyases, with their flavin chromophore (FAD), utilize blue light for DNA repair and photoreduction. The latter process involves two single-electron transfers to FAD with an intermittent protonation step to prime the enzyme active for DNA repair. Here we use time-resolved serial femtosecond X-ray crystallography to describe how light-driven electron transfers trigger subsequent nanosecond-to-microsecond entanglement between FAD and its Asn/Arg-Asp redox sensor triad. We found that this key feature within the photolyase-cryptochrome family regulates FAD re-hybridization and protonation. After first electron transfer, the FAD^•- isoalloxazine ring twists strongly when the arginine closes in to stabilize the negative charge. Subsequent breakage of the arginine-aspartate salt bridge allows proton transfer from arginine to FAD^•-. Our molecular videos demonstrate how the protein environment of redox cofactors organizes multiple electron/proton transfer events in an ordered fashion, which could be applicable to other redox systems such as photosynthesis.

Complex Interplay of Heme-Copper Oxidases with Nitrite and Nitric Oxide

Nitrite and nitric oxide (NO), two active and critical nitrogen oxides linking nitrate to dinitrogen gas in the broad nitrogen biogeochemical cycle, are capable of interacting with redox-sensitive proteins. The interactions of both with heme-copper oxidases (HCOs) serve as the foundation not only for the enzymatic interconversion of nitrogen oxides but also for the inhibitory activity. From extensive studies, we now know that NO interacts with HCOs in a rapid and reversible manner, either competing with oxygen or not. During interconversion, a partially reduced heme/copper center reduces the nitrite ion, producing NO with the heme serving as the reductant and the cupric ion providing a Lewis acid interaction with nitrite. The interaction may lead to the formation of either a relatively stable nitrosyl-derivative of the enzyme reduced or a more labile nitrite-derivative of the enzyme oxidized through two different pathways, resulting in enzyme inhibition. Although nitrite and NO show similar biochemical properties, a growing body of evidence suggests that they are largely treated as distinct molecules by bacterial cells. NO seemingly interacts with all hemoproteins indiscriminately, whereas nitrite shows high specificity to HCOs. Moreover, as biologically active molecules and signal molecules, nitrite and NO directly affect the activity of different enzymes and are perceived by completely different sensing systems, respectively, through which they are linked to different biological processes. Further attempts to reconcile this apparent contradiction could open up possible avenues for the application of these nitrogen oxides in a variety of fields, the pharmaceutical industry in particular.

Proton Pumping and Non-Pumping Terminal Respiratory Oxidases: Active Sites Intermediates of These Molecular Machines and Their Derivatives

Terminal respiratory oxidases are highly efficient molecular machines. These most important bioenergetic membrane enzymes transform the energy of chemical bonds released during the transfer of electrons along the respiratory chains of eukaryotes and prokaryotes from cytochromes or quinols to molecular oxygen into a transmembrane proton gradient. They participate in regulatory cascades and physiological anti-stress reactions in multicellular organisms. They also allow microorganisms to adapt to low-oxygen conditions, survive in chemically aggressive environments and acquire antibiotic resistance. To date, three-dimensional structures with atomic resolution of members of all major groups of terminal respiratory oxidases, heme-copper oxidases, and bd-type cytochromes, have been obtained. These groups of enzymes have different origins and a wide range of functional significance in cells. At the same time, all of them are united by a catalytic reaction of four-electron reduction in oxygen into water which proceeds without the formation and release of potentially dangerous ROS from active sites. The review analyzes recent structural and functional studies of oxygen reduction intermediates in the active sites of terminal respiratory oxidases, the features of catalytic cycles, and the properties of the active sites of these enzymes.

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.

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

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

In Silico Modeling of the Mitochondrial Pumping Complexes with Markov State Models

The mechanism of proton pumping by the mitochondrial electron transport chain complexes is still enigmatic after decades of research. Recently, there has been interest in in silico Markov state models to model the mitochondrial pumping complexes at the microscopic level, and this chapter describes the methods of constructing and simulating such models.

Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations

Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O₂ reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R₄CO, R₂CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R₄CO to R₂CO. We performed potential jump experiments to accumulate the R₂CO and R₄CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a₃. Our experiments also suggest propionate A of heme a₃ changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site.

Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling

Heme-copper oxidases (HCO), nitric oxide reductases (NOR), and sulfite reductases (SiR) catalyze the multi-electron and multi-proton reductions of O2, NO, and SO32-, respectively. Each of these reactions is important to drive cellular energy production through respiratory metabolism and HCO, NOR, and SiR evolved to contain heteronuclear active sites containing heme/copper, heme/nonheme iron, and heme-[4Fe-4S] centers, respectively. The complexity of the structures and reactions of these native enzymes, along with their large sizes and/or membrane associations, make it challenging to fully understand the crucial structural features responsible for the catalytic properties of these active sites. In this review, we summarize progress that has been made to better understand these heteronuclear metalloenzymes at the molecular level though study of the native enzymes along with insights gained from biomimetic models comprising either small molecules or proteins. Further understanding the reaction selectivity of these enzymes is discussed through comparisons of their similar heteronuclear active sites, and we offer outlook for further investigations.

A Resonance Raman Marker Band Characterizes the Slow and Fast Form of Cytochrome c Oxidase

Cytochrome c oxidase (CcO) in its as-isolated form is known to exist in a slow and fast form, which differ drastically in their ability to bind oxygen and other ligands. While preparation methods have been established that yield either the fast or the slow form of the protein, the underlying structural differences have not been identified yet. Here, we have performed surface enhanced resonance Raman (SERR) spectroscopy of CcO immobilized on electrodes in both forms. SERR spectra obtained in resonance with the heme a₃ metal-to-ligand charge transfer (MLCT) transition at 650 nm displayed a sharp vibrational band at 748 or 750 cm^-1 when the protein was in its slow or fast form, respectively. DFT calculations identified the band as a mode of the His-419 ligand that is sensitive to the oxygen ligand and the protonation state of Tyr-288 within the binuclear complex. Potential-dependent SERR spectroscopy showed a redox-induced change of this band around 525 mV versus Ag/AgCl exclusively for the fast form, which coincides with the redox potential of the Tyr-O/Tyr-O^- transition. Our data points to a peroxide ligand in the resting state of CcO for both forms. The observed frequencies and redox sensitivities of the Raman marker band suggest that a radical Tyr-288 is present in the fast form and a protonated Tyr-288 in the slow form.

Protonation Dynamics in the K-Channel of Cytochrome c Oxidase Estimated from Molecular Dynamics Simulations

Proton transfer reactions are one of the most fundamental processes in biochemistry. We present a simplistic approach for estimating proton transfer probabilities in a membrane protein, cytochrome c oxidase. We combine short molecular dynamics simulations at discrete protonation states with a Monte Carlo approach to exchange between those states. Requesting for a proton transfer the existence of a hydrogen-bonded connection between the two source and target residues of the exchange, restricts the acceptance of transfers to only those in which a proton-relay is possible. Together with an analysis of the hydrogen-bonded connectivity in one of the proton-conducting channels of cytochrome c oxidase, this approach gives insight into the protonation dynamics of the hydrogen-bonded networks. The connectivity and directionality of the networks are coupled to the conformation of an important protein residue in the channel, K362, rendering proton transfer in the entire channel feasible in only one of the two major conformations. Proton transport in the channel can thus be regulated by K362 not only through its possible role as a proton carrier itself, but also by allowing or preventing proton transport via water residues.

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

A quantum chemical approach for the mechanisms of redox-active metalloenzymes

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes. In our group we have developed an approach that has been applied to a large number of such systems. Hybrid density functional theory (hybrid DFT) has from the start of these investigations been the leading electronic structure tool. An understanding of how the method works in practice has significantly improved the accuracy and applicability. During the past ten years, it has been found that the results for redox enzymes mainly depend on the chosen fraction of exact exchange in the functional, and that a choice of 15% has worked best. The idea has therefore been to vary that fraction over a reasonable range and study the relative energy dependence. For modeling the enzymes, a cluster approach has been developed. In the present review the development of the method we used is described from its start in work on photosystem II, fifteen years ago. Examples from a few recent applications are described, where the metals have been iron, nickel, copper, cobalt or manganese. The results are in excellent agreement with available experiments, and a large number of new predictions have been made.

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes.

Probing the Proton-Loading Site of Cytochrome C Oxidase Using Time-Resolved Fourier Transform Infrared Spectroscopy

Crystal structure analyses at atomic resolution and FTIR spectroscopic studies of cytochrome c oxidase have yet not revealed protonation or deprotonation of key sites of proton transfer in a time-resolved mode. Here, a sensitive technique to detect protolytic transitions is employed. In this work, probing a proton-loading site of cytochrome c oxidase from Paracoccus denitrificans with time-resolved Fourier transform infrared spectroscopy is presented for the first time. For this purpose, variants with single-site mutations of N131V, D124N, and E278Q, the key residues in the D-channel, were studied. The reaction of mutated CcO enzymes with oxygen was monitored and analyzed. Seven infrared bands in the “fast” kinetic spectra were found based on the following three requirements: (1) they are present in the “fast” phases of N131V and D124N mutants, (2) they have reciprocal counterparts in the “slow” kinetic spectra in these mutants, and (3) they are absent in “fast” kinetic spectra of the E278Q mutant. Moreover, the double-difference spectra between the first two mutants and E278Q revealed more IR bands that may belong to the proton-loading site protolytic transitions. From these results, it is assumed that several polar residues and/or water molecule cluster(s) share a proton as a proton-loading site. This site can be propionate itself (holding only a fraction of H+), His403, and/or water cluster(s).

A Water Molecule Residing in the Fea33+···CuB2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study

Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.

X-ray structures of catalytic intermediates of cytochrome c oxidase provide insights into its O2 activation and unidirectional proton-pump mechanisms

Cytochrome c oxidase (CcO) reduces O2 to water, coupled with a proton-pumping process. The structure of the O2-reduction site of CcO contains two reducing equivalents, Fe a32+ and CuB1+, and suggests that a peroxide-bound state (Fe a33+-O--O--CuB2+) rather than an O2-bound state (Fe a32+-O2) is the initial catalytic intermediate. Unexpectedly, however, resonance Raman spectroscopy results have shown that the initial intermediate is Fe a32+-O2, whereas Fe a33+-O--O--CuB2+ is undetectable. Based on X-ray structures of static noncatalytic CcO forms and mutation analyses for bovine CcO, a proton-pumping mechanism has been proposed. It involves a proton-conducting pathway (the H-pathway) comprising a tandem hydrogen-bond network and a water channel located between the N- and P-side surfaces. However, a system for unidirectional proton-transport has not been experimentally identified. Here, an essentially identical X-ray structure for the two catalytic intermediates (P and F) of bovine CcO was determined at 1.8 Å resolution. A 1.70 Å Fe-O distance of the ferryl center could best be described as Fe a34+ = O2-, not as Fe a34+-OH- The distance suggests an ∼800-cm-1 Raman stretching band. We found an interstitial water molecule that could trigger a rapid proton-coupled electron transfer from tyrosine-OH to the slowly forming Fe a33+-O--O--CuB2+ state, preventing its detection, consistent with the unexpected Raman results. The H-pathway structures of both intermediates indicated that during proton-pumping from the hydrogen-bond network to the P-side, a transmembrane helix closes the water channel connecting the N-side with the hydrogen-bond network, facilitating unidirectional proton-pumping during the P-to-F transition.

Molecular mechanism of mitochondrial respiratory chain assembly and its relation to mitochondrial diseases

The mitochondrial respiratory chain (MRC) is comprised of ~92 nuclear and mitochondrial DNA-encoded protein subunits that are organized into five different multi-subunit respiratory complexes. These complexes produce 90% of the ATP required for cell sustenance. Specific sets of subunits are assembled in a modular or non-modular fashion to construct the MRC complexes. The complete assembly process is gradually chaperoned by a myriad of assembly factors that must coordinate with several other prosthetic groups to reach maturity, makingthe entire processextensively complicated. Further, the individual respiratory complexes can be integrated intovarious giant super-complexes whose functional roles have yet to be explored. Mutations in the MRC subunits and in the related assembly factors often give rise to defects in the proper assembly of the respiratory chain, which then manifests as a group of disorders called mitochondrial diseases, the most common inborn errors of metabolism. This review summarizes the current understanding of the biogenesis of individual MRC complexes and super-complexes, and explores how mutations in the different subunits and assembly factors contribute to mitochondrial disease pathology.

The structure of the oxidized state of cytochrome c oxidase - experiments and theory compared

Cytochrome c oxidase (CcO), the terminal enzyme in the respiratory chain, reduces molecular oxygen to water. Experimental data on the midpoint potentials of the heme iron/copper active site cofactors do not match the overall reaction energetics, and are also in conflict with the observed efficiency of energy conservation in CcO. Therefore it has been postulated that the ferric/cupric intermediate (the oxidized state) exists in two forms. One form, labelled O_H, is presumably involved during catalytic turnover, and should have a high Cu_B midpoint potential due to a metastable high energy structure. When no more electrons are supplied, the O_H state supposedly relaxes to the resting form, labelled O, with a lower energy and a lower midpoint potential. It has been suggested that there is a pure geometrical difference between the O_H and O states, obtained by moving a water molecule inside the active site. It is shown here that the difference between the two forms of the oxidized state must be of a more chemical nature. The reason is that all types of geometrically relaxed structures of the oxidized intermediate have similar energies, all with a high proton coupled reduction potential in accordance with the postulated O_H state. One hypothesized chemical modification of the O_H state is the transfer of an extra proton, possibly internal, into the active site. Such a protonated state has several properties that agree with experimental data on the relaxed oxidized state, including a decreased midpoint potential.

Consolidation of Nath's torsional mechanism of ATP synthesis and two-ion theory of energy coupling in oxidative phosphorylation and photophosphorylation

In a recent publication, Manoj raises criticisms against consensus views on the ATP synthase. The radical statements and assertions are shown to contradict a vast body of available knowledge that includes i) pioneering single-molecule biochemical and biophysical studies from the respected experimental groups of Kinosita, Yoshida, Noji, Börsch, Dunn, Gräber, Frasch, and Dimroth etc., ii) state-of-the-art X-ray and EM/cryo-EM structural information garnered over the decades by the expert groups of Leslie-Walker, Kühlbrandt, Mueller, Meier, Rubinstein, Sazanov, Duncan, and Pedersen on ATP synthase, iii) the pioneering energy-based computer simulations of Warshel, and iv) the novel theoretical and experimental works of Nath. Valid objections against Mitchell's chemiosmotic theory and Boyer's binding change mechanism put forth by Manoj have been addressed satisfactorily by Nath's torsional mechanism of ATP synthesis and two-ion theory of energy coupling and published 10 to 20 years ago, but these papers are not cited by him. This communication shows conclusively and in great detail that none of his objections apply to Nath's mechanism/theory. Nath's theory is further consolidated based on its previous predictive record, its consistency with biochemical evidence, its unified nature, its application to other related energy transductions and to disease, and finally its ability to guide the design of new experiments. Some constructive suggestions for high-resolution structural experiments that have the power to delve into the heart of the matter and throw unprecedented light on the nature of coupled ion translocation in the membrane-bound F_O portion of F₁F_O-ATP synthase are made.

Active Site Midpoint Potentials in Different Cytochrome c Oxidase Families: A Computational Comparison

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain, reducing molecular oxygen to water. The binuclear active site in C cO comprises a high-spin heme associated with a Cu_B complex and a redox active tyrosine. The electron transport in the respiratory chain is driven by increasing midpoint potentials of the involved cofactors, resulting in a release of free energy, which is stored by coupling the electron transfer to proton translocation across a membrane, building up an electrochemical gradient. In this context, the midpoint potentials of the active site cofactors in the C cOs are of special interest, since they determine the driving forces for the individual oxygen reduction steps and thereby affect the efficiency of the proton pumping. It has been difficult to obtain useful information on some of these midpoint potentials from experiments. However, since each of the reduction steps in the catalytic cycle of oxygen reduction to water corresponds to the formation of an O-H bond, they can be calculated with a reasonably high accuracy using quantum chemical methods. From the calculated O-H bond strengths, the proton-coupled midpoint potentials of the active site cofactors can be estimated. Using models representing the different families of C cO's (A, B, and C), the calculations give midpoint potentials that should be relevant during catalytic turnover. The calculations also suggest possible explanations for why some experimentally measured potentials deviate significantly from the calculated ones, i.e., for Cu_B in all oxidase families, and for heme b₃ in the C family.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

Control of transmembrane charge transfer in cytochrome c oxidase by the membrane potential

The respiratory chain in mitochondria is composed of membrane-bound proteins that couple electron transfer to proton translocation across the inner membrane. These charge-transfer reactions are regulated by the proton electrochemical gradient that is generated and maintained by the transmembrane charge transfer. Here, we investigate this feedback mechanism in cytochrome c oxidase in intact inner mitochondrial membranes upon generation of an electrochemical potential by hydrolysis of ATP. The data indicate that a reaction step that involves proton uptake to the catalytic site and presumably proton translocation is impaired by the potential, but electron transfer is not affected. These results define the order of electron and proton-transfer reactions and suggest that the proton pump is regulated by the transmembrane electrochemical gradient through control of internal proton transfer rather than by control of electron transfer.

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.

The electron distribution in the "activated" state of cytochrome c oxidase

Cytochrome c oxidase catalyzes reduction of O₂ to H₂O at a catalytic site that is composed of a copper ion and heme group. The reaction is linked to translocation of four protons across the membrane for each O₂ reduced to water. The free energy associated with electron transfer to the catalytic site is unequal for the four electron-transfer events. Most notably, the free energy associated with reduction of the catalytic site in the oxidized cytochrome c oxidase (state O) is not sufficient for proton pumping across the energized membrane. Yet, this electron transfer is mechanistically linked to proton pumping. To resolve this apparent discrepancy, a high-energy oxidized state (denoted O_H) was postulated and suggested to be populated only during catalytic turnover. The difference between states O and O_H was suggested to be manifested in an elevated midpoint potential of Cu_B in the latter. This proposal predicts that one-electron reduction of cytochrome c oxidase after its oxidation would yield re-reduction of essentially only Cu_B. Here, we investigated this process and found ~5% and ~6% reduction of heme a₃ and Cu_B, respectively, i.e. the apparent redox potentials for heme a₃ and Cu_B are lower than that of heme a.

Spectral components of detergent-solubilized bovine cytochrome oxidase

Cytochrome oxidase is the terminal oxidase of the mitochondrial electron transport chain and pumps 4 protons per oxygen reduced to water. Spectral shifts in the α-band of heme a have been observed in multiple studies and these shifts have the potential to shed light on the proton pumping intermediates. Previously we found that heme a had two spectral components in the α-band during redox titrations in living RAW 264.7 mouse macrophage cells, the classical 605 nm form and a blue-shifted 602 nm form. To confirm these spectral changes were not an artifact due to the complex milieu of the living cell, redox titrations were performed in the isolated detergent-solubilized bovine enzyme from both the Soret- and α-band using precise multiwavelength spectroscopy. This data verified the presence of the 602 nm form in the α-band, revealed a similar shift of heme a in the Soret-band and ruled out the reversal of calcium binding as the origin of the blue shift. The 602 nm form was found to be stabilized at high pH or by binding of azide, which is known to blue shift the α-band of heme a. Azide also stabilized the 602 nm form in the living cells. It is concluded there is a form of cytochrome oxidase in which heme a undergoes a blue shift to a 602 nm form and that redox titrations can be successfully performed in living cells where the oxidase operates in its authentic environment and in the presence of a proton motive force.

Dielectric Properties of a Protein Probed by Reversal of a Buried Ion Pair

Thirty years ago, Hwang and Warshel suggested that a microenvironment preorganized to stabilize an ion pair would be incapable of reorganizing to stabilize the reverse ion pair. The implications were that (1) proteins have a limited capacity to reorganize, even under the influence of strong interactions, such as those present when ionizable groups are buried in the hydrophobic interior of a protein, and (2) the inability of proteins to tolerate the reversal of buried ion pairs demonstrates the limitations inherent to continuum electrostatic models of proteins. Previously we showed that when buried individually in the interior of staphylococcal nuclease, Glu23 and Lys36 have p K_a values near pH 7, but when buried simultaneously, they establish a strong interaction of ∼5 kcal/mol and have p K_a values shifted toward more normal values. Here, using equilibrium thermodynamic measurements, crystal structures, and NMR spectroscopy experiments, we show that although the reversed, individual substitutions-Lys23 and Glu36-also have p K_a values near 7, when buried together, they neither establish a strong interaction nor promote reorganization of their microenvironment. These experiments both confirm Warshel's original hypothesis and expand it by showing that it applies to reorganization, as demonstrated by our artificial ion pairs, as well as to preorganization as is commonly argued for motifs that stabilize naturally occurring ion pairs in polar microenvironments. These data constitute a challenging benchmark useful to test the ability of structure-based algorithms to reproduce the compensation between self-energy, Coulomb and polar interactions in hydrophobic environments of proteins.

Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase (CcO) from the cellular inside to the binuclear redox centre as well as proton pumping through the membrane takes place through proton entrance via two distinct pathways, the D- and K-channel. Both channels show a dependence of their hydration level on the protonation states of their key residues, K362 for the K-channel, and E286 or D132 for the D-channel. In the oxidative half of CcO's catalytic cycle the D-channel is the proton-conducting path. For this channel, an interplay of protonation state of the D-channel residues with the water and hydrogen-bond dynamics has been observed in molecular dynamics simulations of the CcO protein, embedded in a lipid bi-layer, modelled in different protonation states. Protonation of residue E286 at the end of the D-channel results in a hydrogen-bonded network pointing from E286 to N139, that is against proton transport, and favouring N139 conformations which correspond to a closed asparagine gate (formed by residues N121 and N139). Consequently, the hydration level is lower than with unprotonated E286. In those models, the Asn gate is predominantly open, allowing water molecules to pass and thus increase the hydration level. The hydrogen-bonded network in these states exhibits longer life times of the Asn residues with water than other models and shows the D-channel to be traversable from the entrance, D132, to exit, E286. The D-channel can thus be regarded as auto-regulated with respect to proton transport, allowing proton passage only when required, that is the proton is located at the lower part of the D-channel (D132 to Asn gate) and not at the exit (E286).