Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel.

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.     

Molecular basis of proton blockage in aquaporins

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.     

Free energy profiles for H+ conduction along hydrogen-bonded chains of water molecules.

The molecular mechanism for proton conduction along hydrogen-bonded chains, or "proton wires," is studied with free energy simulations. The complete transport of a charge along a proton wire requires two complementary processes: 1) translocation of an excess proton (propagation of an ionic defect), and 2) reorientation of the hydrogen-bonded chain (propagation of a bonding defect). The potential of mean force profile for these two steps is computed in model systems comprising a single-file chain of nine dissociable and polarizable water molecules represented by the PM6 model of Stillinger and co-workers. Results of molecular dynamics simulations with umbrella sampling indicate that the unprotonated chain is preferably polarized, and that the inversion of its total dipole moment involves an activation free energy of 8 kcal/mol. In contrast, the rapid translocation of an excess H+ across a chain extending between two spherical solvent droplets is an activationless process. These results suggest that the propagation of a bonding defect constitutes a limiting step for the passage of several protons along single-file chains of water molecules, whereas the ionic translocation may be fast enough to occur within the lifetime of transient hydrogen-bonded water chains in biological membranes.     

Regeneration and inhibition of proton pumping activity of bacteriorhodopsin blue membrane by cationic amine anesthetics

Bacteriorhodopsin (bR) is the prototype of an integral membrane protein with seven membrane-spanning alpha-helices and serves as a model of the G-protein-coupled drug receptors. This study is aimed at reaching a greater understanding of the role of amine local anesthetic cations on the proton transport in the bR protein, and furthermore, the functional role of "the cation" in the proton pumping mechanism. The effect of the amine anesthetic cations on the proton pump in the bR blue membrane was compared with those by divalent (Ca2+, Mg2+ and Mn2+) and monovalent metal cations (Li+, Na+, K+ and Cs+), which are essential for the correct functioning of the proton pumping of the bR protein. The results suggest that the interacting site of the divalent cation to the bR membrane may differ from that of the monovalent metal cation. The electric current profile of the bR blue membrane in the presence of the amine anesthetic cations was biphasic, involving the generation and inhibition of the proton pumping activity in a concentration-dependent manner. The extent of the regeneration of the proton pump by the additives increased in the order of monovalent metal cation相似文献   

The formation and dynamics of proton wires in channel environments

By virtue of an accurate interaction model, the equilibrium and dynamical properties of an excess proton in aqueous systems are studied, in which the water and excess proton are confined to hydrophobic cylindrical channels. Solvation structures of the excess proton and its mobility along the channel are considered as a function of the channel radius. It is found that when the aqueous proton systems are sufficiently constricted there is a substantial increase in the diffusion of the excess proton charge accompanied by a decrease in the diffusion of water molecules along the channel. Such systems present clear evidence for the possible existence of "proton wires."     

Enhancement of water flow across a carbon nanotube

It is important to investigate how to enhance the flow rate of single-file water molecules across nanochannels. To our knowledge, all the existing methods are based on pressure gradients, external point charges or uniform/graded electric fields. Accordingly, these methods are all based on exogenous tools, and thus bring challenges for both energy-saving and miniaturisation. In contrast, here we manage to reveal an endogenously determined mechanism of flow enhancement. On the basis of molecular dynamics simulations, we investigate water permeation across a single-walled carbon nanotube (SWCNT A) in the presence of another SWCNT (SWCNT B). We find that the flow rate of the single-file water molecules across SWCNT A is enhanced for the case of unblocked SWCNT B compared with the rate for the case of blocked SWCNT B, and that this flow rate is determined by the separation between the two SWCNTs and the diameter of SWCNT B.     

Water alignment, dipolar interactions, and multiple proton occupancy during water-wire proton transport

A discrete multistate kinetic model for water-wire proton transport is constructed and analyzed using Monte Carlo simulations. In the model, each water molecule can be in one of three states: oxygen lone-pairs pointing leftward, pointing rightward, or protonated (H(3)O(+)). Specific rules for transitions among these states are defined as protons hop across successive water oxygens. Our model also includes water-channel interactions that preferentially align the water dipoles, nearest-neighbor dipolar coupling interactions, and Coulombic repulsion. Extensive Monte Carlo simulations were performed and the observed qualitative physical behaviors discussed. We find the parameters that allow the model to exhibit superlinear and sublinear current-voltage relationships, and show why alignment fields, whether generated by interactions with the pore interior or by membrane potentials, always decrease the proton current. The simulations also reveal a "lubrication" mechanism that suppresses water dipole interactions when the channel is multiply occupied by protons. This effect can account for an observed sublinear-to-superlinear transition in the current-voltage relationship.     

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.     

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(3)/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(2) reduction.     

Proteolipid of adenosinetriphosphatase from yeast mitochondria forms proton-selective channels in planar lipid bilayers

Proteolipid isolated from yeast mitochondrial adenosinetriphosphatase by butanol extraction is reincorporated into lipid vesicles from which planar membranes are formed. The proteolipid permits electric conductance through the membrane. This conductance occurs through membrane channels which are highly selective for protons. Proton channels in the membrane are directly observed at high proton concentrations in the aqueous phases. Channels open and close independently from each other; their open-state conductances and lifetimes are monodisperse but influenced by the applied voltage (12 pS and 3 s, respectively, at pH 2.2 and 100 mV). Proton channels do not occur in single proteolipid molecules; the conducting structure consists of at least two polypeptide chains since channels form in a (reversible) bimolecular reaction of nonconducting forms of proteolipid. The number of proton channels at a constant proteolipid concentration changes in sharp transitions and by orders of magnitudes upon critical changes of membrane composition and pH. These transitions are caused by transitions of proteolipid organization in the membrane from a dispersed state (equilibrium between channel-forming "dimers" and a large pool of "monomers") to a state of almost complete aggregation of proteolipid which stabilizes large proton-conducting structures (probably associates of channel-forming dimers). This self-association of isolated proteolipid into structures containing proton-selective channels suggests that the six proteolipids in the adenosinetriphosphatase complex exist as a self-associating entity containing most likely three proton channels.     

Membrane topology of Alzheimer's disease-related presenilin 1. Evidence for the existence of a molecular species with a seven membrane-spanning and one membrane-embedded structure.

A significant member of early-onset familial type of Alzheimer's disease cases has been shown to be caused by dominant mutations in either of the two genes encoding presenilin 1 (PS1) and presenilin 2 (PS2). These two proteins are highly homologous to each other and have been reported to be mainly localized to the membranes of intracellular compartments such as the endoplasmic reticulum. Information about the membrane topological structures of these proteins is indispensable for understanding their physiological and pathological roles. Although several models have been proposed previously, their precise membrane topologies remain unknown. In this study, we examined this issue in detail by expressing a series of C-terminally deleted PS1 mutants fused to the hydrophilic portion of Escherichia coli leader peptidase in vitro using a reticulocyte lysate in the presence of microsomal membranes. Our results predict that PS1 exists mainly in a seven membrane-spanning structure with its C-terminal end exposed to the luminal space. This was also confirmed by expressing these fusion proteins in cultured cells. We further showed that a ninth hydrophobic segment is tightly bound to the membrane without spanning it. Based on the above observations, we propose a novel "seven membrane-spanning and one membrane-embedded" topological model for presenilins.     

Redox-driven proton pumping by heme-copper oxidases

One of the key problems of molecular bioenergetics is the understanding of the function of redox-driven proton pumps on a molecular level. One such class of proton pumps are the heme-copper oxidases. These enzymes are integral membrane proteins in which proton translocation across the membrane is driven by electron transfer from a low-potential donor, such as, e.g. cytochrome c, to a high-potential acceptor, O(2). Proton pumping is associated with distinct exergonic reaction steps that involve gradual reduction of oxygen to water. During the process of O(2) reduction, unprotonated high pK(a) proton acceptors are created at the catalytic site. Initially, these proton acceptors become protonated as a result of intramolecular proton transfer from a residue(s) located in the membrane-spanning part of the enzyme, but removed from the catalytic site. This residue is then reprotonated from the bulk solution. In cytochrome c oxidase from Rhodobacter sphaeroides, the proton is initially transferred from a glutamate, E(I-286), which has an apparent pK(a) of 9.4. According to a recently published structure of the enzyme, the deprotonation of E(I-286) is likely to result in minor structural changes that propagate to protonatable groups on the proton output (positive) side of the protein. We propose that in this way, the free energy available from the O(2) reduction is conserved during the proton transfer. On the basis of the observation of these structural changes, a possible proton-pumping model is presented in this paper. Initially, the structural changes associated with deprotonation of E(I-286) result in the transfer of a proton to an acceptor for pumped protons from the input (negative) side of the membrane. After reprotonation of E(I-286) this acceptor releases a proton to the output side of the membrane.     

Oxygen-oxygen distances in protein-bound crystallographic water suggest the presence of protonated clusters

BackgroundThe availability of high-resolution X-ray structures has shown that proteins contain numerous water molecules, but their role is still not fully understood. Protonated and deprotonated water species are often involved in biochemical reactions. However protons are exceedingly difficult to detect directly because they are electron-poor species.MethodsThe oxygen‑oxygen distance of the crystallographic water molecules was analyzed in a large high-resolution data set. Moreover, a detailed analysis was carried out on the protein-bound water in the available structures of carbonic anhydrase II and cytochrome c oxidase, chosen as protein models in which protonated and deprotonated water species play a significant role.ResultsThe analysis shows an excess of water-water distances below the expected value for hydrogen bond. In the cavities and on the surface of the considered model proteins, clusters of water molecules are found, whose structure suggests the presence of chemical species deriving from self-ionization of water.ConclusionsThe presence of a small maximum below the hydrogen bond threshold in the oxygen‑oxygen distance distribution of crystallographic water molecules, along with the location of many of these water clusters, suggest the presence of Zundel-like structures in, or near, the proteins. Particularly significant is the presence of such structures in protein regions which have been identified as proton antennae or channels.General significanceThis work shows the possibilities, still unexplored, offered by this type of analysis in detecting in structures obtained by X-ray diffraction the presence of aqueous protons or hydroxide ions, which are chemical species as important as elusive.     

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.     

Determinants of Water Permeability through Nanoscopic Hydrophilic Channels

Naturally occurring pores show a variety of polarities and sizes that are presumably directly linked to their biological function. Many biological channels are selective toward permeants similar or smaller in size than water molecules, and therefore their pores operate in the regime of single-file water pores. Intrinsic factors affecting water permeability through such pores include the channel-membrane match, the structural stability of the channel, the channel geometry and channel-water affinity. We present an extensive molecular dynamics study on the role of the channel geometry and polarity on the water osmotic and diffusive permeability coefficients. We show that the polarity of the naturally occurring peptidic channels is close to optimal for water permeation, and that the water mobility for a wide range of channel polarities is essentially length independent. By systematically varying the geometry and polarity of model hydrophilic pores, based on the fold of gramicidin A, the water density, occupancy, and permeability are studied. Our focus is on the characterization of the transition between different permeation regimes in terms of the structure of water in the pores, the average pore occupancy and the dynamics of the permeating water molecules. We show that a general relationship between osmotic and diffusive water permeability coefficients in the single-file regime accounts for the time averaged pore occupancy, and that the dynamics of the permeating water molecules through narrow non single file channels effectively behaves like independent single-file columns.     

Topography of diphtheria Toxin's T domain in the open channel state

When diphtheria toxin encounters a low pH environment, the channel-forming T domain undergoes a poorly understood conformational change that allows for both its own membrane insertion and the translocation of the toxin's catalytic domain across the membrane. From the crystallographic structure of the water-soluble form of diphtheria toxin, a "double dagger" model was proposed in which two transmembrane helical hairpins, TH5-7 and TH8-9, anchor the T domain in the membrane. In this paper, we report the topography of the T domain in the open channel state. This topography was derived from experiments in which either a hexahistidine (H6) tag or biotin moiety was attached at residues that were mutated to cysteines. From the sign of the voltage gating induced by the H6 tag and the accessibility of the biotinylated residues to streptavidin added to the cis or trans side of the membrane, we determined which segments of the T domain are on the cis or trans side of the membrane and, consequently, which segments span the membrane. We find that there are three membrane-spanning segments. Two of them are in the channel-forming piece of the T domain, near its carboxy terminal end, and correspond to one of the proposed "daggers," TH8-9. The other membrane-spanning segment roughly corresponds to only TH5 of the TH5-7 dagger, with the rest of that region lying on or near the cis surface. We also find that, in association with channel formation, the amino terminal third of the T domain, a hydrophilic stretch of approximately 70 residues, is translocated across the membrane to the trans side.     

Converting fructose to 5-hydroxymethylfurfural: a quantum mechanics/molecular mechanics study of the mechanism and energetics

We studied the energetics of the closed-ring mechanism of the acid-catalysed dehydration of d-fructose to 5-hydroxymethylfurfural (HMF) by carrying out canonical ensemble free-energy calculations using bias-sampling, hybrid Quantum Mechanics/Molecular Mechanics Molecular Dynamics simulations with explicit water solvent at 363 K. The quantum mechanical calculations are performed at the PM3 theory level. We find that the reaction proceeds via intramolecular proton and hydride transfers. Solvent dynamics effects are analysed, and we show that the activation energy for the hydride transfers is due to re-organization of the polar solvent environment. We also find that in some instances intramolecular proton transfer is facilitated by mediating water, whereas in others the presence of quantum mechanical water has no effect. From a micro-kinetic point of view, we find that the rate-determining step of the reaction involves a hydride transfer prior to the third dehydration step, requiring an activation free energy of 31.8 kcal/mol, and the respective rate is found in good agreement with reported experimental values in zeolites. Thermodynamically, the reaction is exothermic by .