Kinetic gating of the proton pump in cytochrome c oxidase

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Thermodynamically consistent determination of free energies and rates in kinetic cycle models

Kinetic and thermodynamic models of biological systems are commonly used to connect microscopic features to system function in a bottom-up multiscale approach. The parameters of such models-free energy differences for equilibrium properties and in general rates for equilibrium and out-of-equilibrium observables-have to be measured by different experiments or calculated from multiple computer simulations. All such parameters necessarily come with uncertainties so that when they are naively combined in a full model of the process of interest, they will generally violate fundamental statistical mechanical equalities, namely detailed balance and an equality of forward/backward rate products in cycles due to T. Hill. If left uncorrected, such models can produce arbitrary outputs that are physically inconsistent. Here we develop a maximum likelihood approach (named multibind ) based on the so-called potential graph to combine kinetic or thermodynamic measurements to yield state resolved models that are thermodynamically consistent while being most consistent with the provided data and their uncertainties. We demonstrate the approach with two theoretical models, a generic two-proton binding site and a simplified model of a sodium/proton antiporter. We also describe an algorithm to use the multibind approach to solve the inverse problem of determining microscopic quantities from macroscopic measurements and as an example we predict the microscopic p K a s and protonation states of a small organic molecule from 1D NMR data. The multibind approach is applicable to any thermodynamic or kinetic model that describes a system as transitions between well-defined states with associated free energy differences or rates between these states. A Python package multibind , which implements the approach described here, is made publicly available under the MIT Open Source license.

WHY IT MATTERS

The increase in computational efficiency and rapid advances in methodology for quantitative free energy and rate calculations has allowed for the construction of increasingly complex thermodynamic or kinetic "bottom-up" models of chemical and biological processes. These multi-scale models serve as a framework for analyzing aspects of cellular function in terms of microscopic, molecular properties and provide an opportunity to connect molecular mechanisms to cellular function. The underlying model parameters-free energy differences or rates-are constrained by thermodynamic identities over cycles of states but these identities are not necessarily obeyed during model construction, thus potentially leading to inconsistent models. We address these inconsistencies through the use of a maximum likelihood approach for free energies and rates to adjust the model parameters in such a way that they are maximally consistent with the input parameters and exactly fulfill the thermodynamic cycle constraints. This approach enables formulation of thermodynamically consistent multi-scale models from simulated or experimental measurements.

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., they use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. We show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near-reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

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.

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.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

The redox-coupled proton-channel opening in cytochrome c oxidase

The interplay of cytochrome c oxidase's cofactor electrostatics, long-range conformational changes, H-bond rearrangement, and water dynamics enables transient proton-channel activation.

Kinetics and Efficiency of Energy-Transducing Enzymes

Complexes I to IV, with the exception of Complex II, are redox-driven proton pumps that convert redox energy of oxygen reduction to proton gradient across the mitochondrial or bacterial membrane; in turn, the created electrochemical gradient drives the adenosine triphosphate synthesis in the cells by utilizing complex V of the chain. Here we address a general question of the efficiency of such enzymes, considering them as molecular machines that couple endergonic and exergonic reactions and converting one form of free energy into another. One well-known example of the efficiency is given by Carnot's theorem for heat engines. Here we extend the concept to respiratory enzymes and specifically focus on the proton pumping by Complex I of the respiratory chain, nicotinamide adenine dinucleotide dehydrogenase. To discuss the efficiency issues, we develop a model of enzyme kinetics, which generalizes the Michaelis-Menten model. Our model includes several substrates and products and, in general, can be considered as Generalized Michaelis-Menten Kinetic model. The model might be useful for describing complex enzyme kinetics, regardless of the efficiency issues that are addressed in this paper.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K_v1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of −70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Redox-Driven Proton Pumps of the Respiratory Chain

In aerobic cells, the proton gradient that drives ATP synthesis is created by three different proton pumps-membrane enzymes of the respiratory electron transport chain known as complex I, III, and IV. Despite the striking dissimilarity of structures and apparent differences in molecular mechanisms of proton pumping, all three enzymes have much in common and employ the same universal physical principles of converting redox energy to proton pumping. In this study, we describe a simple mathematical model that illustrates the general principles of redox-driven proton pumps and discuss their implementation in complex I, III, and IV of the respiratory chain.

Tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor

We present a tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor. Our derivation is based on the continuous diffusion equation for overdamped Brownian motion on a time-independent tilted periodic potential in two dimensions. The free-energy potential is nonseparable to allow coupling between the chemical and mechanical degrees of freedom. We formally discretize the chemical coordinate by expanding in Wannier states that are localized along the chemical coordinate and parametrized along the mechanical coordinate. A discrete-continuous equation is derived that is valid for anisotropic systems with weak mechanochemical coupling and deep potential wells along the chemical coordinate. The discrete-continuous description is consistent with established theoretical models of molecular motors with discrete chemical states but is constrained by the underlying continuous two-dimensional potential. In particular, we derive analytic expressions for the effective potential along the mechanical coordinate and for the rate of thermal hopping between chemical states. We determine the thermodynamic efficiency of energy conversion and find that, for a molecular motor with one chemical state per cycle, the derived discrete-continuous equation can accurately describe mechanochemical coupling but cannot describe energy conversion.

Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pK_a is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.

Mitochondrial cytochrome c oxidase: catalysis, coupling and controversies

Mitochondrial cytochrome c oxidase is a member of a diverse superfamily of haem-copper oxidases. Its mechanism of oxygen reduction is reviewed in terms of the cycle of catalytic intermediates and their likely chemical structures. This reaction cycle is coupled to the translocation of protons across the inner mitochondrial membrane in which it is located. The likely mechanism by which this occurs, derived in significant part from studies of bacterial homologues, is presented. These mechanisms of catalysis and coupling, together with current alternative proposals of underlying mechanisms, are critically reviewed.

Submolecular Gates Self-Assemble for Hot-Electron Transfer in Proteins

Redox reactions play key roles in fundamental biological processes. The related spatial organization of donors and acceptors is assumed to undergo evolutionary optimization facilitating charge mobilization within the relevant biological context. Experimental information from submolecular functional sites is needed to understand the organization strategies and driving forces involved in the self-development of structure-function relationships. Here we exploit chemically resolved electrical measurements (CREM) to probe the atom-specific electrostatic potentials (ESPs) in artificial arrays of bacteriochlorophyll (BChl) derivatives that provide model systems for photoexcited (hot) electron donation and withdrawal. On the basis of computations we show that native BChl's in the photosynthetic reaction center (RC) self-assemble at their ground-state as aligned gates for functional charge transfer. The combined computational and experimental results further reveal how site-specific polarizability perpendicular to the molecular plane enhances the hot-electron transport. Maximal transport efficiency is predicted for a specific, ∼5 Å, distance above the center of the metalized BChl, which is in remarkably close agreement with the distance and mutual orientation of corresponding native cofactors. These findings provide new metrics and guidelines for analysis of biological redox centers and for designing charge mobilizing machines such as artificial photosynthesis.

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.

Intrinsic irreversibility limits the efficiency of multidimensional molecular motors

We consider the efficiency limits of Brownian motors able to extract work from the temperature difference between reservoirs or from external thermodynamic forces. These systems can operate in a variety of modes, including as isothermal engines, heat engines, refrigerators, and heat pumps. We derive analytical results showing that certain classes of multidimensional Brownian motor, including the Smoluchowski-Feynman ratchet, are unable to attain perfect efficiency (Carnot efficiency for heat engines). This demonstrates the presence of intrinsic irreversibilities in their operating mechanism. We present numerical simulations showing that in some cases the loss process that limits efficiency is associated with vortices in the probability current.

Molecular simulation and modeling of complex I

Molecular modeling and molecular dynamics simulations play an important role in the functional characterization of complex I. With its large size and complicated function, linking quinone reduction to proton pumping across a membrane, complex I poses unique modeling challenges. Nonetheless, simulations have already helped in the identification of possible proton transfer pathways. Simulations have also shed light on the coupling between electron and proton transfer, thus pointing the way in the search for the mechanistic principles underlying the proton pump. In addition to reviewing what has already been achieved in complex I modeling, we aim here to identify pressing issues and to provide guidance for future research to harness the power of modeling in the functional characterization of complex I. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.

RhNRG-1β Protects the Myocardium against Irradiation-Induced Damage via the ErbB2-ERK-SIRT1 Signaling Pathway

Radiation-induced heart disease (RIHD), which is a serious side effect of the radiotherapy applied for various tumors due to the inevitable irradiation of the heart, cannot be treated effectively using current clinical therapies. Here, we demonstrated that rhNRG-1β, an epidermal growth factor (EGF)-like protein, protects myocardium tissue against irradiation-induced damage and preserves cardiac function. rhNRG-1β effectively ameliorated irradiation-induced myocardial nuclear damage in both cultured adult rat-derived cardiomyocytes and rat myocardium tissue via NRG/ErbB2 signaling. By activating ErbB2, rhNRG-1β maintained mitochondrial integrity, ATP production, respiratory chain function and the Krebs cycle status in irradiated cardiomyocytes. Moreover, the protection of irradiated cardiomyocytes and myocardium tissue by rhNRG-1β was at least partly mediated by the activation of the ErbB2-ERK-SIRT1 signaling pathway. Long-term observations further showed that rhNRG-1β administered in the peri-irradiation period exerts continuous protective effects on cardiac pump function, the myocardial energy metabolism, cardiomyocyte volume and interstitial fibrosis in the rats receiving radiation via NRG/ErbB2 signaling. Our findings indicate that rhNRG-1β can protect the myocardium against irradiation-induced damage and preserve cardiac function via the ErbB2-ERK-SIRT1 signaling pathway.

Proton translocation in cytochrome c oxidase: insights from proton exchange kinetics and vibrational spectroscopy

Cytochrome c oxidase is the terminal enzyme in the electron transfer chain. It reduces oxygen to water and harnesses the released energy to translocate protons across the inner mitochondrial membrane. The mechanism by which the oxygen chemistry is coupled to proton translocation is not yet resolved owing to the difficulty of monitoring dynamic proton transfer events. Here we summarize several postulated mechanisms for proton translocation, which have been supported by a variety of vibrational spectroscopic studies. We recently proposed a proton translocation model involving proton accessibility to the regions near the propionate groups of the heme a and heme a3 redox centers of the enzyme based by hydrogen/deuterium (H/D) exchange Raman scattering studies (Egawa et al., PLoS ONE 2013). To advance our understanding of this model and to refine the proton accessibility to the hemes, the H/D exchange dependence of the heme propionate group vibrational modes on temperature and pH was measured. The H/D exchange detected at the propionate groups of heme a3 takes place within a few seconds under all conditions. In contrast, that detected at the heme a propionates occurs in the oxidized but not the reduced enzyme and the H/D exchange is pH-dependent with a pKa of ~8.0 (faster at high pH). Analysis of the thermodynamic parameters revealed that, as the pH is varied, entropy/enthalpy compensation held the free energy of activation in a narrow range. The redox dependence of the possible proton pathways to the heme groups is discussed. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

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.

The causes of reduced proton-pumping efficiency in type B and C respiratory heme-copper oxidases, and in some mutated variants of type A

The heme-copper oxidases may be divided into three categories, A, B, and C, which include cytochrome c and quinol-oxidising enzymes. All three types are known to be proton pumps and are found in prokaryotes, whereas eukaryotes only contain A-type cytochrome c oxidase in their inner mitochondrial membrane. However, the bacterial B- and C-type enzymes have often been reported to pump protons with an H(+)/e(-) ratio of only one half of the unit stoichiometry in the A-type enzyme. We will show here that these observations are likely to be the result of difficulties with the measuring technique together with a higher sensitivity of the B- and C-type enzymes to the protonmotive force that opposes pumping. We find that under optimal conditions the H(+)/e(-) ratio is close to unity in all the three heme-copper oxidase subfamilies. A higher tendency for proton leak in the B- and C-type enzymes may result from less efficient gating of a proton pump mechanism that we suggest evolved before the so-called D-channel of proton transfer. There is also a discrepancy between results using whole bacterial cells vs. phospholipid vesicles inlaid with oxidase with respect to the observed proton pumping after modification of the D-channel residue asparagine-139 (Rhodobacter sphaeroides numbering) to aspartate in A-type cytochrome c oxidase. This discrepancy might also be explained by a higher sensitivity of proton pumping to protonmotive force in the mutated variant. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

The role of protein dynamics and thermal fluctuations in regulating cytochrome c/cytochrome c oxidase electron transfer

In this overview we present recent combined electrochemical, spectroelectrochemical, spectroscopic and computational studies from our group on the electron transfer reactions of cytochrome c and of the primary electron acceptor of cytochrome c oxidase, the CuA site, in biomimetic complexes. Based on these results, we discuss how protein dynamics and thermal fluctuations may impact on protein ET reactions, comment on the possible physiological relevance of these results, and finally propose a regulatory mechanism that may operate in the Cyt/CcO electron transfer reaction in vivo. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Mutations in the D-channel of cytochrome c oxidase causes leakage of the proton pump

It has experimentally been found that certain mutations close to the entry point of the proton transfer channel in cytochrome c oxidase stop proton translocation but not the oxygen reduction chemistry. This effect is termed uncoupling. Since the mutations are 20Å away from the catalytic center, this is very surprising. A new explanation for this phenomenon is suggested here, involving a local effect at the entry point of the proton channel, rather than the long range effects suggested earlier.

Energy transfer in a molecular motor in the Kramers regime

We present a theoretical treatment of energy transfer in a molecular motor described in terms of overdamped Brownian motion on a multidimensional tilted periodic potential. The tilt represents a thermodynamic force driving the system out of equilibrium and, for nonseparable potentials, energy transfer occurs between degrees of freedom. For deep potential wells, the continuous theory transforms to a discrete master equation that is tractable analytically. We use this master equation to derive formal expressions for the hopping rates, drift and diffusion, and the efficiency and rate of energy transfer in terms of the thermodynamic force. These results span both strong and weak coupling between degrees of freedom, describe the near and far from equilibrium regimes, and are consistent with generalized detailed balance and the Onsager relations. We thereby derive a number of diverse results for molecular motors within a single theoretical framework.

Redox-state sensing by hydrogen bonds in the CuA center of cytochrome c oxidase

Cytochrome c oxidases (CcO) couple electron transfer to active proton translocation through a gated mechanism that minimizes energy losses by preventing protons from flowing backwards or leaking. Such a complex mechanism requires that information about the redox and protonation states of the different centers be transmitted between different parts of the oxidase. Here we report a network of residues located around the electron entry point of CcO, the CuA site in subunit II, that experience collective pH equilibria around neutral pH. This network starts at the occluded side of the CuA site and extends to the interface between subunits I and II of the CcO, where the proton exit is located and through which electrons flow into subunit I. One of the residues in this network is directly involved in a hydrogen bond to one of the CuA ligands, whose strength is highly sensitive to the redox state of the metal center. We propose that this interaction mediates the transmission of redox changes from ET centers to other functional regions of the oxidase, and possibly also in other similar machineries, as part of their gating and regulatory mechanisms.

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.

Optical spectroscopy in turbid media using an integrating sphere: mitochondrial chromophore analysis during metabolic transitions

Recent evidence suggests that the activity of mitochondrial oxidative phosphorylation complexes (MOPCs) is modulated at multiple sites. Here, a method of optically monitoring electron distribution within and between MOPCs is described using a center-mounted sample in an integrating sphere (to minimize scattering effects) with a rapid-scanning spectrometer. The redox-sensitive MOPC absorbances (∼465-630 nm) were modeled using linear least squares analysis with individual chromophore spectra. Classical mitochondrial activity transitions (e.g., ADP-induced increase in oxygen consumption) were used to characterize this approach. Most notable in these studies was the observation that intermediates of the catalytic cycle of cytochrome oxidase are dynamically modulated with metabolic state. The MOPC redox state, along with measurements of oxygen consumption and mitochondrial membrane potential, was used to evaluate the conductances of different sections of the electron transport chain. This analysis then was applied to mitochondria isolated from rabbit hearts subjected to ischemia/reperfusion (I/R). Surprisingly, I/R resulted in an inhibition of all measured MOPC conductances, suggesting a coordinated down-regulation of mitochondrial activity with this well-established cardiac perturbation.

Effect of calcium on the oxidative phosphorylation cascade in skeletal muscle mitochondria

Calcium is believed to regulate mitochondrial oxidative phosphorylation, thereby contributing to the maintenance of cellular energy homeostasis. Skeletal muscle, with an energy conversion dynamic range of up to 100-fold, is an extreme case for evaluating the cellular balance of ATP production and consumption. This study examined the role of Ca(2+) in the entire oxidative phosphorylation reaction network in isolated skeletal muscle mitochondria and attempted to extrapolate these results back to the muscle, in vivo. Kinetic analysis was conducted to evaluate the dose-response effect of Ca(2+) on the maximal velocity of oxidative phosphorylation (V(maxO)) and the ADP affinity. Force-flow analysis evaluated the interplay between energetic driving forces and flux to determine the conductance, or effective activity, of individual steps within oxidative phosphorylation. Measured driving forces [extramitochondrial phosphorylation potential (ΔG(ATP)), membrane potential, and redox states of NADH and cytochromes b(H), b(L), c(1), c, and a,a(3)] were compared with flux (oxygen consumption) at 37 °C; 840 nM Ca(2+) generated an ~2-fold increase in V(maxO) with no change in ADP affinity (~43 μM). Force-flow analysis revealed that Ca(2+) activation of V(maxO) was distributed throughout the oxidative phosphorylation reaction sequence. Specifically, Ca(2+) increased the conductance of Complex IV (2.3-fold), Complexes I and III (2.2-fold), ATP production/transport (2.4-fold), and fuel transport/dehydrogenases (1.7-fold). These data support the notion that Ca(2+) activates the entire muscle oxidative phosphorylation cascade, while extrapolation of these data to the exercising muscle predicts a significant role of Ca(2+) in maintaining cellular energy homeostasis.

Hydrogen-bonded networks along and bifurcation of the E-pathway in quinol:fumarate reductase

The E-pathway of transmembrane proton transfer has been demonstrated previously to be essential for catalysis by the diheme-containing quinol:fumarate reductase (QFR) of Wolinella succinogenes. Two constituents of this pathway, Glu-C180 and heme b(D) ring C (b(D)-C-) propionate, have been validated experimentally. Here, we identify further constituents of the E-pathway by analysis of molecular dynamics simulations. The redox state of heme groups has a crucial effect on the connectivity patterns of mobile internal water molecules that can transiently support proton transfer from the b(D)-C-propionate to Glu-C180. The short H-bonding paths formed in the reduced states can lead to high proton conduction rates and thus provide a plausible explanation for the required opening of the E-pathway in reduced QFR. We found evidence that the b(D)-C-propionate group is the previously postulated branching point connecting proton transfer to the E-pathway from the quinol-oxidation site via interactions with the heme b(D) ligand His-C44. An essential functional role of His-C44 is supported experimentally by site-directed mutagenesis resulting in its replacement with Glu. Although the H44E variant enzyme retains both heme groups, it is unable to catalyze quinol oxidation. All results obtained are relevant to the QFR enzymes from the human pathogens Campylobacter jejuni and Helicobacter pylori.

Energy-converting respiratory Complex I: on the way to the molecular mechanism of the proton pump

In respiring organisms the major energy transduction flux employs the transmembrane electrochemical proton gradient as a physical link between exergonic redox reactions and endergonic ADP phosphorylation. Establishing the gradient involves electrogenic, transmembrane H(+) translocation by the membrane-embedded respiratory complexes. Among others, Complex I (NADH:ubiquinone oxidoreductase) is the most structurally complex and functionally enigmatic respiratory enzyme; its molecular mechanism is as yet unknown. Here we highlight recent progress and discuss the catalytic events during Complex I turnover in relation to their role in energy conversion and to the enzyme structure.

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

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.

Kinetics of the terminal electron transfer step in cytochrome c oxidase

Cytochrome c oxidase (cco) catalyzes the oxygen reduction reaction in most aerobically respiring organisms. Decades of research have uncovered many aspects relating to structure and function of this enzyme. However, the origin of the unusually fast terminal electron transfer step from heme a to heme a(3) in cco has been the subject of intense discussions over recent years. Yet, no satisfactory consensus has been achieved. Carrying out large-scale molecular dynamics simulation of the protein embedded in a solvated membrane, we obtain a reorganization free energy λ = 0.57 eV. Evaluation of the quantized single-mode rate equation using the experimental rate and the computed reorganization free energy gives a value of 1.5 meV for the average electronic coupling (H(ab)) between heme a and heme a(3). Thus, according to our calculations, the nanosecond electron transfer (ET) is due to a small but significant activation barrier (ΔG(‡) = 0.12 eV) in combination with effective electronic coupling between the two cofactors. The activation free energy is caused predominantly by collective reorganization of protein residues. We show that our results are consistent with the weak temperature dependence observed in experiment if one allows for very minor variations in the donor-acceptor distance as the temperature changes.

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.

Glu-286 rotation and water wire reorientation are unlikely the gating elements for proton pumping in cytochrome C oxidase

One of the key unresolved issues regarding proton pumping in cytochrome c oxidase (CcO) is the identity of the gating element that prevents the backflow of protons. In this study, we analyze two popular proposals for this element: isomerization of the key branching residue (Glu-286) and (re)orientation of water molecules in the hydrophobic cavity. Using a multifaceted set of computational analyses that involve CcO embedded in either an implicit or explicit treatment of lipid membrane, we show that neither Glu-286 nor active-site water likely constitutes the gating element. Detailed energetic and structural analyses of the simulation results indicate that the gating-relevant properties of these structural motifs observed in previous work are likely a result of the simplified computational models employed in those studies.

Proton-pumping mechanism of cytochrome c oxidase: a kinetic master-equation approach

Cytochrome c oxidase is an efficient energy transducer that reduces oxygen to water and converts the released chemical energy into an electrochemical membrane potential. As a true proton pump, cytochrome c oxidase translocates protons across the membrane against this potential. Based on a wealth of experiments and calculations, an increasingly detailed picture of the reaction intermediates in the redox cycle has emerged. However, the fundamental mechanism of proton pumping coupled to redox chemistry remains largely unresolved. Here we examine and extend a kinetic master-equation approach to gain insight into redox-coupled proton pumping in cytochrome c oxidase. Basic principles of the cytochrome c oxidase proton pump emerge from an analysis of the simplest kinetic models that retain essential elements of the experimentally determined structure, energetics, and kinetics, and that satisfy fundamental physical principles. The master-equation models allow us to address the question of how pumping can be achieved in a system in which all reaction steps are reversible. Whereas proton pumping does not require the direct modulation of microscopic reaction barriers, such kinetic gating greatly increases the pumping efficiency. Further efficiency gains can be achieved by partially decoupling the proton uptake pathway from the active-site region. Such a mechanism is consistent with the proposed Glu valve, in which the side chain of a key glutamic acid shuttles between the D channel and the active-site region. We also show that the models predict only small proton leaks even in the absence of turnover. The design principles identified here for cytochrome c oxidase provide a blueprint for novel biology-inspired fuel cells, and the master-equation formulation should prove useful also for other molecular machines. .

Reversibility and efficiency in electrocatalytic energy conversion and lessons from enzymes

Enzymes are long established as extremely efficient catalysts. Here, we show that enzymes can also be extremely efficient electrocatalysts (catalysts of redox reactions at electrodes). Despite being large and electronically insulating through most of their volume, some enzymes, when attached to an electrode, catalyze electrochemical reactions that are otherwise extremely sluggish (even with the best synthetic catalysts) and require a large overpotential to achieve a useful rate. These enzymes produce high electrocatalytic currents, displayed in single bidirectional voltammetric waves that switch direction (between oxidation and reduction) sharply at the equilibrium potential for the substrate redox couple. Notoriously irreversible processes such as CO(2) reduction are thereby rendered electrochemically reversible--a consequence of molecular evolution responding to stringent biological drivers for thermodynamic efficiency. Enzymes thus set high standards for the catalysts of future energy technologies.