Water-gated mechanism of proton translocation by cytochrome c oxidase

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.

Superprotonic Conductivity in a Metalloporphyrin-Based SMOF (Supramolecular Metal-Organic Framework)

Metal-organic frameworks and supramolecular metal-organic frameworks (SMOFs) exhibit great potential for a broad range of applications taking advantage of the high surface area and pore sizes and tunable chemistry. In particular, metalloporphyrin-based MOFs and SMOFs are becoming of great importance in many fields due to the bioessential functions of these macrocycles that are being mimicked. On the other hand, during the last years, proton-conducting materials have aroused much interest, and those presenting high conductivity values are potential candidates to play a key role in some solid-state electrochemical devices such as batteries and fuel cells. In this way, using metalloporphyrins as building units we have obtained a new crystalline material with formula [H(bipy)]2[(MnTPPS)(H2O)2]·2bipy·14H2O, where bipy is 4,4'-bipyidine and TPPS4- is the meso-tetra(4-sulfonatephenyl) porphyrin. The crystal structure shows a zig-zag water chain along the [100] direction located between the sulfonate groups of the porphyrin. Taking into account those structural features, the compound was tested for proton conduction by complex electrochemical impedance spectroscopy (EIS). The as-obtained conductivity is 1 × 10-2 S·cm-1 at 40 °C and 98% relative humidity, which is a remarkably high value.

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.

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.

Transient water wires mediate selective proton transport in designed channel proteins

Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.

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.

Effects and Influence of External Electric Fields on the Equilibrium Properties of Tautomeric Molecules

In this review, we have attempted to briefly summarize the influence of an external electric field on an assembly of tautomeric molecules and to what experimentally observable effects this interaction can lead to. We have focused more extensively on the influence of an oriented external electric field (OEEF) on excited-state intramolecular proton transfer (ESIPT) from the studies available to date. The possibilities provided by OEEF for regulating several processes and studying physicochemical processes in tautomers have turned this direction into an attractive area of research due to its numerous applications.

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 O₂ 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 a₃ 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.

Mechanism of Side Chain-Controlled Proton Conductivity in Bioinspired Peptidic Nanostructures

Bioinspired peptide assemblies are promising candidates for use as proton-conducting materials in electrochemical devices and other advanced technologies. Progress toward applications requires establishing foundational structure-function relationships for transport in these materials. This experimental-theoretical study sheds light on how the molecular structure and proton conduction are linked in three synthetic cyclic peptide nanotube assemblies that comprise the three canonical basic amino acids (lysine, arginine, and histidine). Experiments find an order of magnitude higher proton conductivity for lysine-containing peptide assemblies compared to histidine and arginine containing assemblies. The simulations indicate that, upon peptide assembly, the basic amino acid side chains are close enough to enable direct proton transfer. The proton transfer kinetics is determined in the simulations to be governed by the structure and flexibility of the side chains. Together, experiments and theory indicate that the proton mobility is the main determinant of proton conductivity, critical for the performance of peptide-based devices.

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.

Architecture of bacterial respiratory chains

Bacteria power their energy metabolism using membrane-bound respiratory enzymes that capture chemical energy and transduce it by pumping protons or Na⁺ ions across their cell membranes. Recent breakthroughs in molecular bioenergetics have elucidated the architecture and function of many bacterial respiratory enzymes, although key mechanistic principles remain debated. In this Review, we present an overview of the structure, function and bioenergetic principles of modular bacterial respiratory chains and discuss their differences from the eukaryotic counterparts. We also discuss bacterial supercomplexes, which provide central energy transduction systems in several bacteria, including important pathogens, and which could open up possible avenues for treatment of disease.

The Redox-Active Tyrosine Is Essential for Proton Pumping in Cytochrome c Oxidase

Cellular respiration involves electron transport via a number of enzyme complexes to the terminal Cytochrome c oxidase (CcO), in which molecular oxygen is reduced to water. The free energy released in the reduction process is used to establish a transmembrane electrochemical gradient, via two processes, both corresponding to charge transport across the membrane in which the enzymes are embedded. First, the reduction chemistry occurring in the active site of CcO is electrogenic, which means that the electrons and protons are delivered from opposite sides of the membrane. Second, the exergonic chemistry is coupled to translocation of protons across the entire membrane, referred to as proton pumping. In the largest subfamily of the CcO enzymes, the A-family, one proton is pumped for every electron needed for the chemistry, making the energy conservation particularly efficient. In the present study, hybrid density functional calculations are performed on a model of the A-family CcOs. The calculations show that the redox-active tyrosine, conserved in all types of CcOs, plays an essential role for the energy conservation. Based on the calculations a reaction mechanism is suggested involving a tyrosyl radical (possibly mixed with tyrosinate character) in all reduction steps. The result is that the free energy released in each reduction step is large enough to allow proton pumping in all reduction steps without prohibitively high barriers when the gradient is present. Furthermore, the unprotonated tyrosine provides a mechanism for coupling the uptake of two protons per electron in every reduction step, i.e. for a secure proton pumping.

The modified Q-cycle: A look back at its development and forward to a functional model

On looking back at a lifetime of research, it is interesting to see, in the light of current progress, how things came to be, and to speculate on how things might be. I am delighted in the context of the Mitchell prize to have that excuse to present this necessarily personal view of developments in areas of my interests. I have focused on the Q-cycle and a few examples showing wider ramifications, since that had been the main interest of the lab in the 20 years since structures became available, - a watershed event in determining our molecular perspective. I have reviewed the evidence for our model for the mechanism of the first electron transfer of the bifurcated reaction at the Q_o-site, which I think is compelling. In reviewing progress in understanding the second electron transfer, I have revisited some controversies to justify important conclusions which appear, from the literature, not to have been taken seriously. I hope this does not come over as nitpicking. The conclusions are important to the final section in which I develop an internally consistent mechanism for turnovers of the complex leading to a state similar to that observed in recent rapid-mix/freeze-quench experiments, reported three years ago. The final model is necessarily speculative but is open to test.

Spatially dependent H-bond dynamics at interfaces of water/biomimetic self-assembled lattice materials

Understanding hydrogen-bond interactions in self-assembled lattice materials is crucial for preparing such materials, but the role of hydrogen bonds (H bonds) remains unclear. To gain insight into H-bond interactions at the materials' intrinsic spatial scale, we investigated ultrafast H-bond dynamics between water and biomimetic self-assembled lattice materials (composed of sodium dodecyl sulfate and β-cyclodextrin) in a spatially resolved manner. To accomplish this, we developed an infrared pump, vibrational sum-frequency generation (VSFG) probe hyperspectral microscope. With this hyperspectral imaging method, we were able to observe that the primary and secondary OH groups of β-cyclodextrin exhibit markedly different dynamics, suggesting distinct H-bond environments, despite being separated by only a few angstroms. We also observed another ultrafast dynamic reflecting a weakening and restoring of H bonds between bound water and the secondary OH of β-cyclodextrin, which exhibited spatial uniformity within self-assembled domains, but heterogeneity between domains. The restoration dynamics further suggest heterogeneous hydration among the self-assembly domains. The ultrafast nature and meso- and microscopic ordering of H-bond dynamics could contribute to the flexibility and crystallinity of the material--two critically important factors for crystalline lattice self-assemblies--shedding light on engineering intermolecular interactions for self-assembled lattice materials.

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.

Modulation of the electron-proton coupling at cytochrome a by the ligation of the oxidized catalytic center in bovine cytochrome c oxidase

Cytochrome a was suggested as the key redox center in the proton pumping process of bovine cytochrome c oxidase (CcO). Recent studies showed that both the structure of heme a and its immediate vicinity are sensitive to the ligation and the redox state of the distant catalytic center composed of iron of cytochrome a₃ (Fe_a3) and copper (Cu_B). Here, the influence of the ligation at the oxidized Fe_a3³⁺-Cu_B²⁺ center on the electron-proton coupling at heme a was examined in the wide pH range (6.5-11). The strength of the coupling was evaluated by the determination of pH dependence of the midpoint potential of heme a (E_m(a)) for the cyanide (the low-spin Fe_a3³⁺) and the formate-ligated CcO (the high-spin Fe_a3³⁺). The measurements were performed under experimental conditions when other three redox centers of CcO are oxidized. Two slightly differing linear pH dependencies of E_m(a) were found for the CN- and the formate-ligated CcO with slopes of -13 mV/pH unit and -23 mV/pH unit, respectively. These linear dependencies indicate only a weak and unspecific electron-proton coupling at cytochrome a in both forms of CcO. The lack of the strong electron-proton coupling at the physiological pH values is also substantiated by the UV-Vis absorption and electron-paramagnetic resonance spectroscopy investigations of the cyanide-ligated oxidized CcO. It is shown that the ligand exchange at Fe_a³⁺ between His-Fe_a³⁺-His and His-Fe_a³⁺-OH^- occurs only at pH above 9.5 with the estimated pK >11.0.

On the beneficent thickness of water

In the 1930s, Lars Onsager published his famous 'reciprocal relations' describing free energy conversion processes. Importantly, these relations were derived on the assumption that the fluxes of the processes involved in the conversion were proportional to the forces (free energy gradients) driving them. For chemical reactions, however, this condition holds only for systems operating close to equilibrium-indeed very close; nominally requiring driving forces to be smaller than k B T. Fairly soon thereafter, however, it was quite inexplicably observed that in at least some biological conversions both the reciprocal relations and linear flux-force dependency appeared to be obeyed no matter how far from equilibrium the system was being driven. No successful explanation of how this 'paradoxical' behaviour could occur has emerged and it has remained a mystery. We here argue, however, that this anomalous behaviour is simply a gift of water, of its viscosity in particular; a gift, moreover, without which life almost certainly could not have emerged. And a gift whose appreciation we primarily owe to recent work by Prof. R. Dean Astumian who, as providence has kindly seen to it, was led to the relevant insights by the later work of Onsager himself.

Proton transfer in uncoupled variants of cytochrome c oxidase

Cytochrome c oxidase is a membrane-bound redox-driven proton pump that harbors two proton-transfer pathways, D and K, which are used at different stages of the reaction cycle. Here, we address the question if a D pathway with a modified energy landscape for proton transfer could take over the role of the K pathway when the latter is blocked by a mutation. Our data indicate that structural alterations near the entrance of the D pathway modulate energy barriers that influence proton transfer to the proton-loading site. The data also suggest that during reduction of the catalytic site, its protonation has to occur via the K pathway and that this proton transfer to the catalytic site cannot take place through the D pathway.

Long-range proton-coupled electron transfer in biological energy conversion: towards mechanistic understanding of respiratory complex I

Biological energy conversion is driven by efficient enzymes that capture, store and transfer protons and electrons across large distances. Recent advances in structural biology have provided atomic-scale blueprints of these types of remarkable molecular machinery, which together with biochemical, biophysical and computational experiments allow us to derive detailed energy transduction mechanisms for the first time. Here, I present one of the most intricate and least understood types of biological energy conversion machinery, the respiratory complex I, and how its redox-driven proton-pump catalyses charge transfer across approximately 300 Å distances. After discussing the functional elements of complex I, a putative mechanistic model for its action-at-a-distance effect is presented, and functional parallels are drawn to other redox- and light-driven ion pumps.

Electric field modulated redox-driven protonation and hydration energetics in energy converting enzymes

Biological energy conversion is catalysed by proton-coupled electron transfer (PCET) reactions that form the chemical basis of respiratory and photosynthetic enzymes. Despite recent advances in structural, biophysical, and computational experiments, the mechanistic principles of these reactions still remain elusive. Based on common functional features observed in redox enzymes, we study here generic mechanistic models for water-mediated long-range PCET reactions. We show how a redox reaction within a buried protein environment creates an electric field that induces hydration changes between the proton acceptor and donor groups, and in turn, lowers the reaction barrier and increases the thermodynamic driving forces for the water-mediated PCET process. We predict linear free energy relationships, and discuss the proposed mechanism in context of PCET in cytochrome c oxidase.

Metalloproteins in the Biology of Heterocysts

Cyanobacteria are photoautotrophic microorganisms present in almost all ecologically niches on Earth. They exist as single-cell or filamentous forms and the latter often contain specialized cells for N₂ fixation known as heterocysts. Heterocysts arise from photosynthetic active vegetative cells by multiple morphological and physiological rearrangements including the absence of O₂ evolution and CO₂ fixation. The key function of this cell type is carried out by the metalloprotein complex known as nitrogenase. Additionally, many other important processes in heterocysts also depend on metalloproteins. This leads to a high metal demand exceeding the one of other bacteria in content and concentration during heterocyst development and in mature heterocysts. This review provides an overview on the current knowledge of the transition metals and metalloproteins required by heterocysts in heterocyst-forming cyanobacteria. It discusses the molecular, physiological, and physicochemical properties of metalloproteins involved in N₂ fixation, H₂ metabolism, electron transport chains, oxidative stress management, storage, energy metabolism, and metabolic networks in the diazotrophic filament. This provides a detailed and comprehensive picture on the heterocyst demands for Fe, Cu, Mo, Ni, Mn, V, and Zn as cofactors for metalloproteins and highlights the importance of such metalloproteins for the biology of cyanobacterial heterocysts.

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.

Cytochrome c oxidase structures suggest a four-state stochastic pump mechanism

A simple stochastic model for a cytochrome c oxidase proton pump.

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.

Alkane guest packing drives switching between multimeric deep-cavity cavitand assembly states

Alkane guest transfer into aqueous dimeric, tetrameric, hexameric, and octameric assemblies of the deep-cavity cavitand TEMOA is examined using molecular simulations. The experimental transitions between aggregation states strongly correlate with calculated alkane transfer free energy minima, demonstrating the guiding role of guest packing on stabilizing multimeric complexes. The predictive simulation approach described affords a salient rationale as to why octameric assemblies have yet to be experimentally observed.

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).

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

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.

Proton Wire Dynamics in the Green Fluorescent Protein

Inside proteins, protons move on proton wires (PWs). Starting from the highest resolution X-ray structure available, we conduct a 306 ns molecular dynamics simulation of the (A-state) wild-type (wt) green fluorescent protein (GFP) to study how its PWs change with time. We find that the PW from the chromophore via Ser205 to Glu222, observed in all X-ray structures, undergoes rapid water molecule insertion between Ser205 and Glu222. Sometimes, an alternate Ser205-bypassing PW exists. Side chain rotations of Thr203 and Ser205 play an important role in shaping the PW network in the chromophore region. Thr203, with its bulkier side chain, exhibits slower transitions between its three rotameric states. Ser205 experiences more frequent rotations, slowing down when the Thr203 methyl group is close by. The combined states of both residues affect the PW probabilities. A random walk search for PWs from the chromophore reveals several exit points to the bulk, one being a direct water wire (WW) from the chromophore to the bulk. A longer WW connects the "bottom" of the GFP barrel with a "water pool" (WP1) situated below Glu222. These two WWs were not observed in X-ray structures of wt-GFP, but their analogues have been reported in related fluorescent proteins. Surprisingly, the high-resolution X-ray structure utilized herein shows that Glu222 is protonated at low temperatures. At higher temperatures, we suggest ion pairing between anionic Glu222 and a proton hosted in WP1. Upon photoexcitation, these two recombine, while a second proton dissociates from the chromophore and either exits the protein using the short WW or migrates along the GFP-barrel axis on the long WW. This mechanism reconciles the conflicting experimental and theoretical data on proton motion within GFP.

Atomistic determinants of co-enzyme Q reduction at the Qi-site of the cytochrome bc1 complex

The cytochrome (cyt) bc₁ complex is an integral component of the respiratory electron transfer chain sustaining the energy needs of organisms ranging from humans to bacteria. Due to its ubiquitous role in the energy metabolism, both the oxidation and reduction of the enzyme’s substrate co-enzyme Q has been studied vigorously. Here, this vast amount of data is reassessed after probing the substrate reduction steps at the Q_i-site of the cyt bc₁ complex of Rhodobacter capsulatus using atomistic molecular dynamics simulations. The simulations suggest that the Lys251 side chain could rotate into the Q_i-site to facilitate binding of half-protonated semiquinone – a reaction intermediate that is potentially formed during substrate reduction. At this bent pose, the Lys251 forms a salt bridge with the Asp252, thus making direct proton transfer possible. In the neutral state, the lysine side chain stays close to the conserved binding location of cardiolipin (CL). This back-and-forth motion between the CL and Asp252 indicates that Lys251 functions as a proton shuttle controlled by pH-dependent negative feedback. The CL/K/D switching, which represents a refinement to the previously described CL/K pathway, fine-tunes the proton transfer process. Lastly, the simulation data was used to formulate a mechanism for reducing the substrate at the Q_i-site.

Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.

A Protonated Water Cluster as a Transient Proton-Loading Site in Cytochrome c Oxidase

Cytochrome c oxidase (CcO) is a redox-driven proton pump that powers aerobic respiratory chains. We show here by multi-scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton-loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps.

Analyzing the electrogenicity of cytochrome c oxidase

Measurements of voltage changes in response to charge separation within membrane proteins can offer fundamental information on spectroscopically "invisible" steps. For example, results from studies of voltage changes associated with electron and proton transfer in cytochrome c oxidase could, in principle, be used to discriminate between different theoretical models describing the molecular mechanism of proton pumping. Earlier analyses of data from these measurements have been based on macroscopic considerations that may not allow for exploring the actual molecular mechanisms. Here, we have used a coarse-grained model describing the relation between observed voltage changes and specific charge-transfer reactions, which includes an explicit description of the membrane, the electrolytes, and the electrodes. The results from these calculations offer mechanistic insights at the molecular level. Our main conclusion is that previously assumed mechanistic evidence that was based on electrogenic measurements is not unique. However, the ability of our calculations to obtain reliable voltage changes means that we have a tool that can be used to describe a wide range of electrogenic charge transfers in channels and transporters, by combining voltage measurements with other experiments and simulations to analyze new mechanistic proposals.

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.

pK(a) of a Proton-Conducting Water Chain in Photosystem II

Recent high-resolution crystal structures of the water-oxidizing enzyme photosystem II (PSII) show that O4 of the catalytic Mn4CaO5 cluster forms an H-bond with a water molecule W539, which belongs to a chain of water molecules (O4-water chain). Oxidation of Mn4CaO5 to S1 resulted in elongation of the O-H bonds and decrease in pKa(O-H/O(-)) in the [O4-H···OW539-H···OW538-H···OW393] region along the O4-water chain. In S1, removal of all water molecules from the O4-water chain, except W539, resulted in a significant pKa upshift at O4; this suggests that the proton-conducting water chain serves as a conducting media for protons and significantly decreases the donor pKa, leading to a downhill proton transfer. The absence of a corresponding proton-conducting channel is disadvantageous for release of protons from the proton-releasing site, as in the case of O5 that has no H-bond partner.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Computationally Efficient Multiscale Reactive Molecular Dynamics to Describe Amino Acid Deprotonation in Proteins

An important challenge in the simulation of biomolecular systems is a quantitative description of the protonation and deprotonation process of amino acid residues. Despite the seeming simplicity of adding or removing a positively charged hydrogen nucleus, simulating the actual protonation/deprotonation process is inherently difficult. It requires both the explicit treatment of the excess proton, including its charge defect delocalization and Grotthuss shuttling through inhomogeneous moieties (water and amino residues), and extensive sampling of coupled condensed phase motions. In a recent paper (J. Chem. Theory Comput. 2014, 10, 2729-2737), a multiscale approach was developed to map high-level quantum mechanics/molecular mechanics (QM/MM) data into a multiscale reactive molecular dynamics (MS-RMD) model in order to describe amino acid deprotonation in bulk water. In this article, we extend the fitting approach (called FitRMD) to create MS-RMD models for ionizable amino acids within proteins. The resulting models are shown to faithfully reproduce the free energy profiles of the reference QM/MM Hamiltonian for PT inside an example protein, the ClC-ec1 H(+)/Cl(-) antiporter. Moreover, we show that the resulting MS-RMD models are computationally efficient enough to then characterize more complex 2-dimensional free energy surfaces due to slow degrees of freedom such as water hydration of internal protein cavities that can be inherently coupled to the excess proton charge translocation. The FitRMD method is thus shown to be an effective way to map ab initio level accuracy into a much more computationally efficient reactive MD method in order to explicitly simulate and quantitatively describe amino acid protonation/deprotonation in proteins.

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.

Mechanism of Oxygen Reduction in Cytochrome c Oxidase and the Role of the Active Site Tyrosine

Cytochrome c oxidase, the terminal enzyme in the respiratory chain, reduces molecular oxygen to water and stores the released energy through electrogenic chemistry and proton pumping across the membrane. Apart from the heme-copper binuclear center, there is a conserved tyrosine residue in the active site (BNC). The tyrosine delivers both an electron and a proton during the O-O bond cleavage step, forming a tyrosyl radical. The catalytic cycle then occurs in four reduction steps, each taking up one proton for the chemistry (water formation) and one proton to be pumped. It is here suggested that in three of the reduction steps the chemical proton enters the center of the BNC, leaving the tyrosine unprotonated with radical character. The reproprotonation of the tyrosine occurs first in the final reduction step before binding the next oxygen molecule. It is also suggested that this reduction mechanism and the presence of the tyrosine are essential for the proton pumping. Density functional theory calculations on large cluster models of the active site show that only the intermediates with the proton in the center of the BNC and with an unprotonated tyrosyl radical have a high electron affinity of similar size as the electron donor, which is essential for the ability to take up two protons per electron and thus for the proton pumping. This type of reduction mechanism is also the only one that gives a free energy profile in accordance with experimental observations for the amount of proton pumping in the working enzyme.