Et tu, Grotthuss! and other unfinished stories

This review article is divided into three sections. In Section 1, a short biographical note on Freiherr von Grotthuss is followed by a detailed summary of the main findings and ideas present in his 1806 paper. Attempts to place Grotthuss contribution in the context of the science done at his time were also made. In Section 2, the modern version of the Grotthuss mechanism is reviewed. The classical Grotthuss model has been recently questioned and new mechanisms and ideas regarding proton transfer are briefly discussed. The last section discusses the significance of a classical Grotthuss mechanism for proton transfer in water chains inside protein cavities. This has been an interesting new twist in the ongoing history of the Grotthuss mechanism. A summary and discussion of what was learned from probably the simplest currently available experimental models of proton transfer in water wires in semi-synthetic ion channels are critically presented. This review ends discussing some of the questions that need to be addressed in the near future.     

Proton transport in carbonic anhydrase: Insights from molecular simulation

This article reviews the insights gained from molecular simulations of human carbonic anhydrase II (HCA II) utilizing non-reactive and reactive force fields. The simulations with a reactive force field explore protein transfer and transport via Grotthuss shuttling, while the non-reactive simulations probe the larger conformational dynamics that underpin the various contributions to the rate-limiting proton transfer event. Specific attention is given to the orientational stability of the His64 group and the characteristics of the active site water cluster, in an effort to determine both of their impact on the maximal catalytic rate. The explicit proton transfer and transport events are described by the multistate empirical valence bond (MS-EVB) method, as are alternative pathways for the excess proton charge defect to enter/leave the active site. The simulation results are interpreted in light of experimental results on the wild-type enzyme and various site-specific mutations of HCA II in order to better elucidate the key factors that contribute to its exceptional efficiency.     

Lewis-inspired representation of dissociable water in clusters and Grotthuss chains

Proton transfer to and from water is critical to the function of water in many settings. However, it has been challenging to model. Here, we present proof-of-principle for an efficient yet robust model based on Lewis-inspired submolecular particles with interactions that deviate from Coulombic at short distances to take quantum effects into account. This “LEWIS” model provides excellent correspondence with experimental structures for water molecules and water clusters in their neutral, protonated and deprotonated forms; reasonable values for the proton affinities of water and hydroxide; a good value for the strength of the hydrogen bond in the water dimer; the correct order of magnitude for the stretch and bend force constants of water; and the expected time course for Grotthuss transport in water chains.     

The mechanism of the proton transfer: an outline

A brief summary of the principal notions of the quantum-mechanical theory of the charge transfer reactions has been presented. In the framework of this theory, the mechanism of the proton transfer consists in the classical medium reorganization that equalizes the proton energy levels in the initial and final states, and a consequent proton transfer via a quantum-mechanical underbarrier transition. On the basis of this mechanism, factors influencing the proton transfer probability, and hence kinetic isotope effect, have been discussed; among them are the optimum tunneling distance, the involvement of the excited vibrational states, etc. Semi-classical and quantum-mechanical treatments of the Swain-Schaad relations have been compared. Some applications to enzymatic proton-transfer reactions have been described.     

A modeling and simulation perspective on the mechanism and function of respiratory complex I

Respiratory complex I is a giant redox-driven proton pump, and central to energy production in mitochondria and bacteria. It catalyses the reduction of quinone to quinol, and converts the free energy released into the endergonic proton translocation across the membrane. The proton pumping sets up the proton electrochemical gradient, which propels the synthesis of ATP. Despite the availability of extensive biochemical, biophysical and structural data on complex I, the mechanism of coupling between the electron and proton transfer reactions remain uncertain. In this work, we discuss current state-of-the-art in the field with particular emphasis on the molecular mechanism of respiratory complex I, as deduced from computational modeling and simulation approaches, but in strong alliance with the experimental data. This leads to novel synthesis of mechanistic ideas on a highly complex enzyme of the electron transport chain that has been associated with a number of mitochondrial and neurodegenerative disorders.     

Charge delocalization in proton channels, I: the aquaporin channels and proton blockage

The explicit contribution to the free energy barrier and proton conductance from the delocalized nature of the excess proton is examined in aquaporin channels using an accurate all-atom molecular dynamics computer simulation model. In particular, the channel permeation free energy profiles are calculated and compared for both a delocalized (fully Grotthuss shuttling) proton and a classical (nonshuttling) hydronium ion along two aquaporin channels, Aqp1 and GlpF. To elucidate the effects of the bipolar field thought to arise from two alpha-helical macrodipoles on proton blockage, free energy profiles were also calculated for computational mutants of the two channels where the bipolar field was eliminated by artificially discharging the backbone atoms. Comparison of the free energy profiles between the proton and hydronium cases indicates that the magnitude of the free energy barrier and position of the barrier peak for the fully delocalized and shuttling proton are somewhat different from the case of the (localized) classical hydronium. The proton conductance through the two aquaporin channels is also estimated using Poisson-Nernst-Planck theory for both the Grotthuss shuttling excess proton and the classical hydronium cation.     

Protonation of excited state pyrene-1-carboxylate by phosphate and organic acids in aqueous solution studied by fluorescence spectroscopy

Pyrene-1-carboxylic acid has a pK of 4.0 in the ground state and 8.1 in the singlet electronic excited state. In the pH range of physiological interest (pH approximately 5-8), the ground state compound is largely ionized as pyrene-1-carboxylate, but protonation of the excited state molecule occurs when a proton donor reacts with the carboxylate during the excited state lifetime of the fluorophore. Both forms of the pyrene derivatives are fluorescent, and in this work the protonation reaction was measured by monitoring steady-state and time-resolved fluorescence. The rate of protonation of pyrene-COO(-) by acetic, chloroacetic, lactic, and cacodylic acids is a function of DeltapK, as predicted by Marcus theory. The rate of proton transfer from these acids saturates at high concentration, as expected for the existence of an encounter complex. Trihydrogen-phosphate is a much better proton donor than dihydrogen- and monohydrogen-phosphate, as can be seen by the pH dependence. The proton-donating ability of phosphate does not saturate at high concentrations, but increases with increasing phosphate concentration. We suggest that enhanced rate of proton transfer at high phosphate concentrations may be due to the dual proton donating and accepting nature of phosphate, in analogy to the Grotthuss mechanism for proton transfer in water. It is suggested that in molecular structures containing multiple phosphates, such as membrane surfaces and DNA, proton transfer rates will be enhanced by this mechanism.     

Water channel structures analysed by electron crystallography

Background

The mechanisms underlying water transport through aquaporin (AQP) have been debated for two decades. The water permeation phenomenon of AQP seems inexplicable because the Grotthuss mechanism does not allow for simultaneous fast water permeability and inhibition of proton transfer through the hydrogen bonds of water molecules.

Scope of review

The AQP1 structure determined by electron crystallography provided the first insights into the proton exclusion mechanism despite fast water permeation. Although several studies have provided clues about the mechanism based on the AQP structure, each proposed mechanism remains incomplete. The present review is focused on AQP function and structure solved by electron crystallography in an attempt to fill the gaps between the findings in the absence and presence of lipids.

Major conclusions

Many AQP structures can be superimposed regardless of the determination method. The AQP fold is preserved even under conditions lacking lipids, but the water arrangement in the channel pore differs. The differences might be explained by dipole moments formed by the two short helices in the lipid bilayer. In addition, structure analyses of double-layered two-dimensional crystals of AQP suggest an array formation and cell adhesive function.

General significance

Electron crystallography findings not only have contributed to resolve some of the water permeation mechanisms, but have also elucidated the multiple functions of AQPs in the membrane. The roles of AQPs in the brain remain obscure, but their multiple activities might be important in the regulation of brain and other biological functions. This article is part of a Special Issue entitled Aquaporins.     

Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a₃/Cu_B site, and for proton pumping. When present, the side chain of Glu-242 is orientated “down” towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity “above” Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an “up” position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the “up” position is preferred energetically when the cavity contains four water molecules, but the “down” position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the “resting” state of the enzyme, and that the “up/down” equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O₂ reduction.     

Active proton leak in mitochondria: A new way to regulate substrate oxidation

The main function of mitochondria is energy transduction, from substrate oxidation to the free energy of ATP synthesis, through oxidative phosphorylation. For physiological reasons, the degree of coupling between these two processes must be modulated in order to adapt redox potential and ATP turnover to cellular needs. Such a modulation leads to energy wastage. To this day, two energy wastage mechanisms have been described: the membrane passive proton conductance (proton leak) and the decrease in the coupling efficiency between electrons transfer and proton extrusion at the proton pumps level (redox or proton slipping). In this paper, we describe a new energy wastage mechanism of interest. We show that in isolated yeast mitochondria, the membrane proton conductance is strictly dependent on the external dehydrogenases activity. An increase in their activity leads to an increase in the membrane proton conductance. This proton permeability is independent of the respiratory chain and ATP synthase proton pumps. We propose to name this new mechanism “active proton leak.” Such a mechanism could allow a wide modulation of substrate oxidation in response to cellular redox constraints.     

Enhancement of proton conductance by mutations of the selectivity filter of aquaporin-1

Prevention of cation permeation in wild-type aquaporin-1 (AQP1) is believed to be associated with the Asn-Pro-Ala (NPA) region and the aromatic/arginine selectivity filter (SF) domain. Previous work has suggested that the NPA region helps to impede proton permeation due to the protein backbone collective macrodipoles that create an environment favoring a directionally discontinuous channel hydrogen-bonded water chain and a large electrostatic barrier. The SF domain contributes to the proton permeation barrier by a spatial restriction mechanism and direct electrostatic interactions. To further explore these various effects, the free-energy barriers and the maximum cation conductance for the permeation of various cations through the AQP1-R195V and AQP1-R195S mutants are predicted computationally. The cations studied included the hydrated excess proton that utilizes the Grotthuss shuttling mechanism, a model “classical” charge localized hydronium cation that exhibits no Grotthuss shuttling, and a sodium cation. The hydrated excess proton was simulated using a specialized multi-state molecular dynamics method including a proper physical treatment of the proton shuttling and charge defect delocalization. Both AQP1 mutants exhibit a surprising cooperative effect leading to a reduction in the free-energy barrier for proton permeation around the NPA region due to altered water configurations in the SF region, with AQP1-R195S having a higher conductance than AQP1-R195V. The theoretical predictions are experimentally confirmed in wild-type AQP1 and the mutants expressed in Xenopus oocytes. The combined results suggest that the SF domain is a specialized structure that has evolved to impede proton permeation in aquaporins.     

Molecular mechanism of H+ conduction in the single-file water chain of the gramicidin channel

The conduction of protons in the hydrogen-bonded chain of water molecules (or "proton wire") embedded in the lumen of gramicidin A is studied with molecular dynamics free energy simulations. The process may be described as a "hop-and-turn" or Grotthuss mechanism involving the chemical exchange (hop) of hydrogen nuclei between hydrogen-bonded water molecules arranged in single file in the lumen of the pore, and the subsequent reorganization (turn) of the hydrogen-bonded network. Accordingly, the conduction cycle is modeled by two complementary steps corresponding respectively to the translocation 1) of an ionic defect (H+) and 2) of a bonding defect along the hydrogen-bonded chain of water molecules in the pore interior. The molecular mechanism and the potential of mean force are analyzed for each of these two translocation steps. It is found that the mobility of protons in gramicidin A is essentially determined by the fine structure and the dynamic fluctuations of the hydrogen-bonded network. The translocation of H+ is mediated by spontaneous (thermal) fluctuations in the relative positions of oxygen atoms in the wire. In this diffusive mechanism, a shallow free-energy well slightly favors the presence of the excess proton near the middle of the channel. In the absence of H+, the water chain adopts either one of two polarized configurations, each of which corresponds to an oriented donor-acceptor hydrogen-bond pattern along the channel axis. Interconversion between these two conformations is an activated process that occurs through the sequential and directional reorientation of water molecules of the wire. The effect of hydrogen-bonding interactions between channel and water on proton translocation is analyzed from a comparison to the results obtained previously in a study of model nonpolar channels, in which such interactions were missing. Hydrogen-bond donation from water to the backbone carbonyl oxygen atoms lining the pore interior has a dual effect: it provides a coordination of water molecules well suited both to proton hydration and to high proton mobility, and it facilitates the slower reorientation or turn step of the Grotthuss mechanism by stabilizing intermediate configurations of the hydrogen-bonded network in which water molecules are in the process of flipping between their two preferred, polarized states. This mechanism offers a detailed molecular model for the rapid transport of protons in channels, in energy-transducing membrane proteins, and in enzymes.     

Energy transport mechanism in the form of proton soliton in a one-dimensional hydrogen-bonded polypeptide chain

The dynamics of protons in a one-dimensional hydrogen-bonded (HB) polypeptide chain (PC) is investigated theoretically. A new Hamiltonian is formulated with the inclusion of higher-order molecular interactions between peptide groups (PGs). The wave function of the excitation state of a single particle is replaced by a new wave function of a two-quanta quasi-coherent state. The dynamics is governed by a higher-order nonlinear Schrödinger equation and the energy transport is performed by the proton soliton. A nonlinear multiple-scale perturbation analysis has been performed and the evolution of soliton parameters such as velocity and amplitude is explored numerically. The proton soliton is thermally stable and very robust against these perturbations. The energy transport by the proton soliton is more appropriate to understand the mechanism of energy transfer in biological processes such as muscle contraction, DNA replication, and neuro-electric pulse transfer on biomembranes.     

Bacteriorhodopsin: a high-resolution structural view of vectorial proton transport

Recent 3-D structures of several intermediates in the photocycle of bacteriorhodopsin (bR) provide a detailed structural picture of this molecular proton pump in action. In this review, we describe the sequence of conformational changes of bR following the photoisomerization of its all-trans retinal chromophore, which is covalently bound via a protonated Schiff base to Lys216 in helix G, to a 13-cis configuration. The initial changes are localized near the protein's active site and a key water molecule is disordered. This water molecule serves as a keystone for the ground state of bR since, within the framework of the complex counter ion, it is important both for stabilizing the structure of the extracellular half of the protein, and for maintaining the high pK_a of the Schiff base (the primary proton donor) and the low pK_a of Asp85 (the primary proton acceptor). Subsequent structural rearrangements propagate out from the active site towards the extracellular half of the protein, with a local flex of helix C exaggerating an early movement of Asp85 towards the Schiff base, thereby facilitating proton transfer between these two groups. Other coupled rearrangements indicate the mechanism of proton release to the extracellular medium. On the cytoplasmic half of the protein, a local unwinding of helix G near the backbone of Lys216 provides sites for water molecules to order and define a pathway for the reprotonation of the Schiff base from Asp96 later in the photocycle. A steric clash of the photoisomerized retinal with Trp182 in helix F drives an outward tilt of the cytoplasmic half of this helix, opening the proton transport channel and enabling a proton to be taken up from the cytoplasm. Although bR is the first integral membrane protein to have its catalytic mechanism structurally characterized in detail, several key results were anticipated in advance of the structural model and the general framework for vectorial proton transport has, by and large, been preserved.     

The quinone chemistry of bc complexes

The quinone chemistry that gives rise to the rather unusual strict bifurcation of electron transfer at the Q(o) site of the cytochrome bc complexes remains controversial. In this article, I review recent ideas and propose a "logic-gated" binding mechanism that combines classical quinone electrochemistry with specific hydrogen bonding requirements and results in a reversible reaction that minimizes unwanted side-reactions that could otherwise undermine the efficiency of the Q-cycle proton/electron coupling mechanism.     

The dynamics and energetics of water permeation and proton exclusion in aquaporins

Aquaporins and aquaglyceroporins are passive membrane channels that, in many species, facilitate highly efficient yet strictly selective permeation of water and small solutes across lipid bilayers. Their ability to block proton flux is particularly remarkable, because other aqueous pores and water efficiently conduct protons, via the so-called Grotthuss mechanism. How efficient water permeation is achieved and how it is reconciled with the seemingly contradictory task of strict proton exclusion have been long-standing puzzles. Because neither the dynamics of the water molecules nor the mobility of protons inside the aquaporin channel could be experimentally accessed so far, several groups addressed this challenge using a variety of atomistic computer simulation methods.     

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.     

Net proton uptake is preceded by multiple proton transfer steps upon electron injection into cytochrome c oxidase

Cytochrome c oxidase (COX), the last enzyme of the respiratory chain of aerobic organisms, catalyzes the reduction of molecular oxygen to water. It is a redox-linked proton pump, whose mechanism of proton pumping has been controversially discussed, and the coupling of proton and electron transfer is still not understood. Here, we investigated the kinetics of proton transfer reactions following the injection of a single electron into the fully oxidized enzyme and its transfer to the hemes using time-resolved absorption spectroscopy and pH indicator dyes. By comparison of proton uptake and release kinetics observed for solubilized COX and COX-containing liposomes, we conclude that the 1-μs electron injection into Cu(A), close to the positive membrane side (P-side) of the enzyme, already results in proton uptake from both the P-side and the N (negative)-side (1.5 H(+)/COX and 1 H(+)/COX, respectively). The subsequent 10-μs transfer of the electron to heme a is accompanied by the release of 1 proton from the P-side to the aqueous bulk phase, leaving ～0.5 H(+)/COX at this side to electrostatically compensate the charge of the electron. With ～200 μs, all but 0.4 H(+) at the N-side are released to the bulk phase, and the remaining proton is transferred toward the hemes to a so-called "pump site." Thus, this proton may already be taken up by the enzyme as early as during the first electron transfer to Cu(A). These results support the idea of a proton-collecting antenna, switched on by electron injection.     

Understanding the mechanism of proton movement linked to oxygen reduction in cytochrome c oxidase: lessons from other proteins

Cytochrome c oxidase is a large intrinsic membrane protein designed to use the energy of electron transfer and oxygen reduction to pump protons across a membrane. The molecular mechanism of the energy conversion process is not understood. Other proteins with simpler, better resolved structures have been more completely defined and offer insight into possible mechanisms of proton transfer in cytochrome c oxidase. Important concepts that are illustrated by these model systems include the ideas of conformational change both close to and at a distance from the triggering event, and the formation of a transitory water-linked proton pathway during a catalytic cycle. Evidence for the applicability of these concepts to cytochrome c oxidase is discussed.     

Proton transfers in the bacteriorhodopsin photocycle

The steps in the mechanism of proton transport in bacteriorhodopsin include examples for most kinds of proton transfer reactions that might occur in a transmembrane pump: proton transfer via a bridging water molecule, coupled protonation/deprotonation of two buried groups separated by a considerable distance, long-range proton migration over a hydrogen-bonded aqueous chain, and capture as well as release of protons at the membrane-water interface. The conceptual and technical advantages of this system have allowed close examination of many of these model reactions, some at an atomic level.