Theoretical Studies of the Grotthuss Mechanism in Biological Proton Wires

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Slowing Down of the Molecular Reorientation of Water in Concentrated Alkaline Solutions

It is generally accepted that the hydroxide ion (OH^-) is a strong hydrogen bond acceptor and that its anomalously high diffusion constant in water results from a Grotthuss-like structural diffusion mechanism. However, the spatial extent over which OH^- ions influence the dynamics of the hydrogen-bond network of water remained largely unclear. Here, we measure the ultrafast dynamics of OH groups of HDO molecules interacting with the deuterated hydroxide ion OD^-. For solutions with OD^- concentrations up to 4 M, we find that HDO molecules that are not directly interacting with the ions have a reorientation time constant of ∼2.7 ps, similar to that of pure liquid water. When the concentration of OD^- ions is increased, the reorientation time constant increases, indicating a strong slowing down of the structural dynamics of the solution.

Hydration interactions beyond the first solvation shell in aqueous phenolate solution

We investigate the orientational dynamics of water molecules solvating phenolate ions using ultrafast vibrational spectroscopy and density functional theory-based molecular dynamics simulations.

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.

ELECTRON TRANSFER REACTIONS COUPLED TO PROTON TRANSLOCATION: CYTOCHROME OXIDASE, PROTON PUMPS, AND BIOLOGICAL ENERGY TRANSDUCTION

Cytochrome oxidase (COX) is the terminal component of electron transport chain of the respiratory system in mitochondria, and one of the key enzymes responsible for energy generation in cells. COX functions as a proton pump that utilizes free energy of oxygen reduction for translocation of protons across the mitochondrion membrane. The proton gradient created in the process is later utilized to drive synthesis of ATP. Although the structure of COX has been recently resolved, the molecular mechanism of proton pumping remains unknown. In this paper, general principles and possible molecular mechanisms of energy transformations in this enzyme will be discussed. The main question is how exactly chemical energy of oxygen reduction and water formation is transformed into a proton gradient; or, how exactly electron transfer reactions are utilized to translocate protons across the mitochondrion membrane against the electrochemical gradient. A key to the solution of this problem is in understanding correlated transport of electrons and protons. Here, theoretical models are discussed for coupled electron and proton transfer reactions in which an electron is tunneling over long distance between two redox cofactors, and a coupled proton is moving along a proton conducting channel in a classical, diffusion-like random walk fashion. Such reactions are typical for COX and other enzymes involved in biological energy transformations.

Voltage-gated proton channels and other proton transfer pathways

Proton channels exist in a wide variety of membrane proteins where they transport protons rapidly and efficiently. Usually the proton pathway is formed mainly by water molecules present in the protein, but its function is regulated by titratable groups on critical amino acid residues in the pathway. All proton channels conduct protons by a hydrogen-bonded chain mechanism in which the proton hops from one water or titratable group to the next. Voltage-gated proton channels represent a specific subset of proton channels that have voltage- and time-dependent gating like other ion channels. However, they differ from most ion channels in their extraordinarily high selectivity, tiny conductance, strong temperature and deuterium isotope effects on conductance and gating kinetics, and insensitivity to block by steric occlusion. Gating of H(+) channels is regulated tightly by pH and voltage, ensuring that they open only when the electrochemical gradient is outward. Thus they function to extrude acid from cells. H(+) channels are expressed in many cells. During the respiratory burst in phagocytes, H(+) current compensates for electron extrusion by NADPH oxidase. Most evidence indicates that the H(+) channel is not part of the NADPH oxidase complex, but rather is a distinct and as yet unidentified molecule.

Computer simulation of water in cytochrome c oxidase

Statistical mechanics and molecular dynamics simulations have been carried out to study the distribution and dynamics of internal water molecules in bovine heart cytochrome c oxidase (CcO). CcO is found to be capable of holding plenty of water, which in subunit I alone amounts to about 165 molecules. The dynamic characterization of these water molecules is carried out. The nascent water molecules produced in the redox reaction at the heme a(3)-CuB binuclear site form an intriguing chain structure. The chain begins at the position of Glu242 at the end of the D channel, and has a fork structure, one branch of which leads to the binuclear center, and the other to the propionate d of heme a(3). The branch that leads to the binuclear center has dynamic access both to the site where the formation of water occurs, and to delta-nitrogen of His291. From the binuclear center, the chain continues to run into the K channel. The stability of this hydrogen bond network is examined dynamically. The catalytic site is located at the hydrophobic region, and the nascent water molecules are produced at the top of the energy hill. The energy gradient is utilized as the mechanism of water removal from the protein. The water exit channels are explored using high-temperature dynamics simulations. Two putative channels for water exit from the catalytic site have been identified. One is leading directly toward Mg(2+) site. However, this channel is only open when His291 is dissociated from CuB. If His291 is bound to CuB, the only channel for water exit is the one that originates at E242 and leads toward the middle of the membrane. This is the same channel that is presumably used for oxygen supply.

The role of Trp side chains in tuning single proton conduction through gramicidin channels

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.

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.

Crystal structure of F65A/Y131C-methylimidazole carbonic anhydrase V reveals architectural features of an engineered proton shuttle

The crystal structure of F65A/Y131C murine alpha-carbonic anhydrase V (CAV), covalently modified at cysteine residues with 4-chloromethylimidazole, is reported at 1.88 A resolution. This modification introduces a methylimidazole (MI) group at residue C131 in the active site with important consequences. F65A/Y131C-MI CAV exhibits an up to 3-fold enhancement of catalytic activity over that of wild-type CAV [Earnhardt, J. N., Wright, S. K., Qian, M., Tu, C., Laipis, P. J., Viola, R. E., and Silverman, D. N. (1999) Arch. Biochem. Biophys. 361, 264-270]. In this modified CAV variant, C131-MI acts as a proton shuttle, facilitating the deprotonation of a zinc-bound water molecule to regenerate the nucleophilic zinc-bound hydroxide ion. A network of three hydrogen-bonded water molecules, across which proton transfer likely proceeds, bridges the zinc-bound water molecule and the C131-MI imidazole group. The structure of F65A/Y131C-MI CAV is compared to structures of Y64H/F65A murine CAV, wild-type human alpha-carbonic anhydrase II, and the gamma-carbonic anhydrase from Methanosarcina thermophilain an effort to outline common features of catalytic proton shuttles.

Thermodynamic view of activation energies of proton transfer in various gramicidin A channels

The temperature dependencies (range: 5-45 degrees C) of single-channel proton conductances (g(H)) in native gramicidin A (gA) and in two diastereoisomers (SS and RR) of the dioxolane-linked gA channels were measured in glycerylmonooleate/decane (GMO) and diphytanoylphosphatidylcholine/decane (DiPhPC) bilayers. Linear Arrhenius plots (ln (g(H)) versus K(-1)) were obtained for the native gA and RR channels in both types of bilayers, and for the SS channel in GMO bilayers only. The Arrhenius plot for proton transfer in the SS channel in DiPhPC bilayers had a break in linearity around 20 degrees C. This break seems to occur only when protons are the permeating cations in the SS channel. The activation energies (E(a)) for proton transfer in various gA channels (approximately 15 kJ/mol) are consistent with the rate-limiting step being in the channel and/or at the membrane-channel/solution interface, and not in bulk solution. E(a) values for proton transfer in gA channels are considerably smaller than for the permeation of nonproton currents in gA as well as in various other ion channels. The E(a) values for proton transfer in native gA channels are nearly the same in both GMO and DiPhPC bilayers. In contrast, for the dioxolane linked gA dimers, E(a) values were strongly modulated by the lipid environment. The Gibbs activation free energies (Delta G(#)(o)) for protons in various gA channels are within the range of 27-29 kJ/mol in GMO bilayers and of 20-22 kJ/mol in DiPhPC bilayers. The largest difference between Delta G(#)(o) for proton currents occurs between native gA (or SS channels) and the RR channel. In general, the activation entropy (Delta S) is mostly responsible for the differences between g(H) values in various gA channels, and also in distinct bilayers. However, significant differences between the activation enthalpies (Delta H(#)(o)) for proton transfer in the SS and RR channels occur in distinct membranes.

Modulation of proton transfer in the water wire of dioxolane-linked gramicidin channels by lipid membranes

Proton conductance (g(H)) in single SS stereoisomers of dioxolane-linked gramicidin A (gA) channels were measured in different phospholipid bilayers at different HCl concentrations. In particular, measurements were obtained in bilayers made of 1,2-diphytanoyl 3-phosphocholine (DiPhPC) or its ethylated derivative 1,2-diphytanoyl 3-ethyl-phosphocholine (et-DiPhPC,). The difference between these phospholipids is that in et-DiPhPC one of the phosphate oxygens is covalently linked to an ethyl group and cannot be protonated. In relatively dilute acid solutions, g(H) in DiPhPC is significantly higher than in et-DiPhPC. At high acid concentrations, g(H) is the same in both diphytanoyl bilayers. Such differences in g(H) can be accounted for by surface charge effects at the membrane/solution interfaces. In the linear portion of the log g(H)-log [H] relationship, g(H) values in diphytanoyl bilayers were significantly larger (approximately 10-fold) than in neutral glyceryl monooleate (GMO) membranes. The slopes of the linear log-log relationships between g(H) and [H] in diphytanoyl and GMO bilayers are essentially the same (approximately 0.76). This slope is significantly lower than the slope of the log-log plot of proton conductivity versus proton concentration in aqueous solutions (approximately 1.00). Because the chemical composition of the membrane-channel/solution interface is strikingly different in GMO and diphytanoyl bilayers, the reduced slope in g(H)-[HCl] relationships may be a characteristic of proton transfer in the water wire inside the SS channel. Values of g(H) in diphytanoyl bilayers were also significantly larger than in membranes made of the more common biological phospholipids 1-palmitoyl 2-oleoyl phosphocholine (POPC) or 1-palmitoyl 2-oleoyl phosphoethanolamine (POPE). These differences, however, cannot be accounted for by different surface charge effects or by different internal dipole potentials. On the other hand, maximum g(H) measured in the SS channel does not depend on the composition of the bilayer and is determined essentially by the reduced mobility of protons in concentrated acid solutions. Finally, no experimental evidence was found in support of a lateral proton movement at the phospholipid/solution interface contributing to g(H) in single SS channels. Protein-lipid interactions are likely to modulate g(H) in the SS channel.

A combined molecular dynamics and diffusion model of single proton conduction through gramicidin

We develop a model for proton conduction through gramicidin based on the molecular dynamics simulations of Pomès and Roux (Biophys. J. 72:A246, 1997). The transport of a single proton through the gramicidin pore is described by a potential of mean force and diffusion coefficient obtained from the molecular dynamics. In addition, the model incorporates the dynamics of a defect in the hydrogen bonding structure of pore waters without an excess proton. Proton entrance and exit were not simulated by the molecular dynamics. The single proton conduction model includes a simple representation of these processes that involves three free parameters. A reasonable value can be chosen for one of these, and the other two can be optimized to yield a good fit to the proton conductance data of, Ann. N.Y. Acad. Sci. 339:8-20) for pH > or = 1.7. A sensitivity analysis shows the significance of this fit.

Proton mobilities in water and in different stereoisomers of covalently linked gramicidin A channels

Proton conductivities in bulk solution (lambda(H)) and single-channel proton conductances (g(H)) in two different stereoisomers of the dioxolane-linked gramicidin A channel (the SS and RR dimers) were measured in a wide range of bulk proton concentrations ([H], 0.1-8000 mM). Proton mobilities (micro(H)) in water as well as in the SS and RR dimers were calculated from the conductivity data. In the concentration range of 0.1-2000 mM, a straight line with a slope of 0.75 describes the log (g(H))-log ([H]) relationship in the SS dimer. At [H] > 2000 mM, saturation is followed by a decline in g(H). The g(H)-[H] relationship in the SS dimer is qualitatively similar to the [H] dependence of lambda(H). However, the slope of the straight line in the log(lambda(H))-log([H]) plot is 0.96, indicating that the rate-limiting step for proton conduction through the SS dimer is not the diffusion of protons in bulk solution. The significant difference between the slopes of those linear relationships accounts for the faster decline of micro(H) as a function of [H] in the SS dimer in relation to bulk solution. In the high range of [H], saturation and decline of g(H) in the SS dimer can be accounted for by the significant decrease of micro(H) in bulk solution. At any given [H], g(H) in the RR dimer is significantly smaller than in the SS. Moreover, the g(H)-[H] relationship in the RR stereoisomer is qualitatively different from that in the SS. Between 1 and 50 mM [H], g(H) can be fitted with an adsorption isotherm, suggesting the presence of a proton-binding site inside the pore (pK(a) approximately 2), which limits proton exit from the channel. At 100 mM < [H] < 3000 mM, g(H) increases linearly with [H]. The distinctive shape of the g(H)-[H] relationship in the RR dimer suggests that the channel can be occupied simultaneously by more than one proton. At higher [H], the saturation and decline of g(H) in the RR dimer reflect the properties of micro(H) in bulk solution. In the entire range of [H], protons seem to cross the SS and RR channels via a Grotthuss-like mechanism. The rate-limiting step for proton transfer in the SS dimer is probably the membrane-channel/bulk solution interface. It is also proposed that the smaller g(H) in the RR dimer is the consequence of a different organization and dynamics of the H-bonded network of water molecules inside the pore of the channel, resulting in a slower proton transfer and multiple pore occupancy by protons.