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.     

Recombinant production of the human aquaporins in the yeast Pichia pastoris (Invited Review)

Abstract

Aquaporins are water facilitating proteins embedded in the cellular membranes. Such channels have been identified in almost every living organism – including humans. These proteins are vital molecules and their malfunction can lead to several severe disorders and diseases. Hence, an increased understanding of their structure, function and regulation is of the utmost importance for developing current and future drugs. Heading towards this goal, the first problem to overcome is to acquire the proteins in sufficient amounts to enable functional and structural characterization. Using a suitable host organism, large amounts of target molecules can possibly be produced, but for membrane proteins limitations are frequently encountered. In the work described here, we have produced the 13 human aquaporins (hAQPs) in one of the most successful hosts for recombinant overproduction of eukaryotic proteins; the yeast Pichia pastoris, in order to explore the underlying bottleneck to a successful membrane protein production experiment. Here we present exceptional yield of hAQP1, whereas some other hAQPs were below the threshold needed for scaled up production. In the overproduction process, we have established methods for efficient production screening as well as for accurate determination of the initial production yield. Furthermore, we have optimized the yield of low producing targets, enabling studies of proteins previously out of reach, exemplified with hAQP4 as well as the homologue PfAQP. Taken together, our results. present insight into factors directing high production of eukaryotic membrane proteins together with suggestions on ways to optimize the recombinant production in the yeast P. pastoris.     

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.     

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.     

Homology modeling of representative subfamilies of Arabidopsis major intrinsic proteins. Classification based on the aromatic/arginine selectivity filter

Major intrinsic proteins (MIPs) are a family of membrane channels that facilitate the bidirectional transport of water and small uncharged solutes such as glycerol. The 35 full-length members of the MIP family in Arabidopsis are segregated into four structurally homologous subfamilies: plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), nodulin 26-like intrinsic membrane proteins (NIPs), and small basic intrinsic proteins (SIPs). Computational methods were used to construct structural models of the putative pore regions of various plant MIPs based on homology modeling with the atomic resolution crystal structures of mammalian aquaporin 1 and the bacterial glycerol permease GlpF. Based on comparisons of the narrow selectivity filter regions (the aromatic/Arg [ar/R] filter), the members of the four phylogenetic subfamilies of Arabidopsis MIPs can be classified into eight groups. PIPs possess a uniform ar/R signature characteristic of high water transport aquaporins, whereas TIPs are highly diverse with three separate conserved ar/R regions. NIPs possess two separate conserved ar/R regions, one that is similar to the archetype, soybean (Glycine max) nodulin 26, and another that is characteristic of Arabidopsis NIP6;1. The SIP subfamily possesses two ar/R subgroups, characteristic of either SIP1 or SIP2. Both SIP ar/R residues are divergent from all other MIPs in plants and other kingdoms. Overall, these findings suggest that higher plant MIPs have a common fold but show distinct differences in proposed pore apertures, potential to form hydrogen bonds with transported molecules, and amphiphilicity that likely results in divergent transport selectivities.     

Structural basis for proton conduction and inhibition by the influenza M2 protein

The influenza M2 protein forms an acid‐activated and drug‐sensitive proton channel in the virus envelope that is important for the virus lifecycle. The functional properties and high‐resolution structures of this proton channel have been extensively studied to understand the mechanisms of proton conduction and drug inhibition. We review biochemical and electrophysiological studies of M2 and discuss how high‐resolution structures have transformed our understanding of this proton channel. Comparison of structures obtained in different membrane‐mimetic solvents and under different pH using X‐ray crystallography, solution NMR, and solid‐state NMR spectroscopy revealed how the M2 structure depends on the environment and showed that the pharmacologically relevant drug‐binding site lies in the transmembrane (TM) pore. Competing models of proton conduction have been evaluated using biochemical experiments, high‐resolution structural methods, and computational modeling. These results are converging to a model in which a histidine residue in the TM domain mediates proton relay with water, aided by microsecond conformational dynamics of the imidazole ring. These mechanistic insights are guiding the design of new inhibitors that target drug‐resistant M2 variants and may be relevant for other proton channels.     

Molecular modeling of the bacterial outer membrane receptor energizer, ExbBD/TonB, based on homology with the flagellar motor, MotAB

The MotA/MotB proteins serve as the motor that drives bacterial flagellar rotation in response to the proton motive force (pmf). They have been shown to comprise a transmembrane proton pathway. The ExbB/ExbD/TonB protein complex serves to energize transport of iron siderophores and vitamin B₁₂ across the outer membrane of the Gram-negative bacterial cell using the pmf. These two protein complexes have the same topology and are homologous. Based on molecular data for the MotA/MotB proteins, we propose simple three-dimensional channel structures for both MotA/MotB and ExbB/ExbD/TonB using modeling methods. Features of the derived channels are discussed, and two possible proton transfer pathways for the ExbBD/TonB system are proposed. These analyses provide a guide for molecular studies aimed at elucidating the mechanism by which chemiosmotic energy can be transferred either between two adjacent membranes to energize outer membrane transport or to the bacterial flagellum to generate torque.     

Molecular modeling of the bacterial outer membrane receptor energizer,ExbBD/TonB,based on homology with the flagellar motor,MotAB

The MotA/MotB proteins serve as the motor that drives bacterial flagellar rotation in response to the proton motive force (pmf). They have been shown to comprise a transmembrane proton pathway. The ExbB/ExbD/TonB protein complex serves to energize transport of iron siderophores and vitamin B12 across the outer membrane of the Gram-negative bacterial cell using the pmf. These two protein complexes have the same topology and are homologous. Based on molecular data for the MotA/MotB proteins, we propose simple three-dimensional channel structures for both MotA/MotB and ExbB/ExbD/TonB using modeling methods. Features of the derived channels are discussed, and two possible proton transfer pathways for the ExbBD/TonB system are proposed. These analyses provide a guide for molecular studies aimed at elucidating the mechanism by which chemiosmotic energy can be transferred either between two adjacent membranes to energize outer membrane transport or to the bacterial flagellum to generate torque.     

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.     

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.     

Phenotypes of aquaporin mutants in genetically altered mice

The ability to move water across lipid membranes is crucial for nutrient intake, energy generation, waste excretion, and a myriad of other functions associated with life. Aquaporins, a family of integral membrane proteins, are now recognized as the channels responsible for transporting hydrophilic molecules, including water, across relatively impervious, hydrophobic cell membranes. A tremendous amount of work has been published on characterizing these proteins, which have been found in all bacteria, yeast, plants, and animals examined to date. In addition, an increasing number of mouse models with genetically altered aquaporin expression are being reported. This article will briefly review the basic biochemistry of aquaporins and then evaluate the use (and misuse) of mice in the quest for understanding the comparative pathophysiology of aquaporins in humans.     

Ion channels as antivirus targets

Ion channels are membrane proteins that are found in a number of viruses and which are of crucial physiological importance in the viral life cycle.They have one common feature in that their action mode involves a change of electrochemical or proton gradient across the bilayer lipid membrane which modulates viral or cellular activity.We will discuss a group of viral channel proteins that belong to the viroproin family,and which participate in a number of viral functions including promoting the release of viral particles from cells.Blocking these channel-forming proteins may be"lethal",which can be a suitable and potential therapeutic strategy.In this review we discuss seven ion channels of viruses which can lead serious infections in human beings: M2 of influenza A,NB and BM2 of influenza B,CM2 of influenza C,Vpu of HIV-1,p7 of HCV and 2B of picornaviruses.     

Comparative genomics and site-directed mutagenesis support the existence of only one input channel for protons in the C-family (cbb3 oxidase) of heme-copper oxygen reductases

Oxygen reductase members of the heme-copper superfamily are terminal respiratory oxidases in mitochondria and many aerobic bacteria and archaea, coupling the reduction of molecular oxygen to water to the translocation of protons across the plasma membrane. The protons required for catalysis and pumping in the oxygen reductases are derived from the cytoplasmic side of the membrane, transferred via proton-conducting channels comprised of hydrogen bond chains containing internal water molecules along with polar amino acid side chains. Recent analyses identified eight oxygen reductase families in the superfamily: the A-, B-, C-, D-, E-, F-, G-, and H-families of oxygen reductases. Two proton input channels, the K-channel and the D-channel, are well established in the A-family of oxygen reductases (exemplified by the mitochondrial cytochrome c oxidases and by the respiratory oxidases from Rhodobacter sphaeroides and Paracoccus denitrificans). Each of these channels can be identified by the pattern of conserved polar amino acid residues within the protein. The C-family (cbb3 oxidases) is the second most abundant oxygen reductase family after the A-family, making up more than 20% of the sequences of the heme-copper superfamily. In this work, sequence analyses and structural modeling have been used to identify likely proton channels in the C-family. The pattern of conserved polar residues supports the presence of only one proton input channel, which is spatially analogous to the K-channel in the A-family. There is no pattern of conserved residues that could form a D-channel analogue or an alternative proton channel. The functional importance of the residues proposed to be part of the K-channel was tested by site-directed mutagenesis using the cbb3 oxidases from R. sphaeroides and Vibrio cholerae. Several of the residues proposed to be part of the putative K-channel had significantly reduced catalytic activity upon mutation: T219V, Y227F/Y228F, N293D, and Y321F. The data strongly suggest that in the C-family only one channel functions for the delivery of both catalytic and pumped protons. In addition, it is also proposed that a pair of acidic residues, which are totally conserved among the C-family, may be part of a proton-conducting exit channel for pumped protons. The residues homologous to these acidic amino acids are highly conserved in the cNOR family of nitric oxide reductases and have previously been implicated as part of a proton-conducting channel delivering protons from the periplasmic side of the membrane to the enzyme active site in the cNOR family. It is possible that the C-family contains a homologous proton-conducting channel that delivers pumped protons in the opposite direction, from the active site to the periplasm.     

Review. Proton-coupled gating in chloride channels

The physiologically indispensable chloride channel (CLC) family is split into two classes of membrane proteins: chloride channels and chloride/proton antiporters. In this article we focus on the relationship between these two groups and specifically review the role of protons in chloride-channel gating. Moreover, we discuss the evidence for proton transport through the chloride channels and explore the possible pathways that the protons could take through the chloride channels. We present results of a mutagenesis study, suggesting the feasibility of one of the pathways, which is closely related to the proton pathway proposed previously for the chloride/proton antiporters. We conclude that the two groups of CLC proteins, although in principle very different, employ similar mechanisms and pathways for ion transport.     

Mitochondrial uncoupling proteins: new insights from functional and proteomic studies

Mitochondria are the major sites of ATP synthesis through oxidative phosphorylation, a process that is weakened by proton leak. Uncoupling proteins are mitochondrial membrane proteins specialized in inducible proton conductance. They dissipate the proton electrochemical gradient established by the respiratory chain at the expense of reducing substrates. Several physiological roles have been suggested for uncoupling proteins, including roles in the control of the cellular energy balance and in preventive action against oxidative stress. This review focuses on new leads emerging from comparative proteomics about the involvement of uncoupling protein in the mitochondrial physiology. A brief overview on uncoupling proteins and on proteomics applied to mitochondria is also presented herein.     

Secondary Water Pore Formation for Proton Transport in a ClC Exchanger Revealed by an Atomistic Molecular-Dynamics Simulation

Several prokaryotic ClC proteins have been demonstrated to function as exchangers that transport both chloride ions and protons simultaneously in opposite directions. However, the path of the proton through the ClC exchanger, and how the protein brings about the coupled movement of both ions are still unknown. In this work, we use an atomistic molecular dynamics (MD) simulation to demonstrate that a previously unknown secondary water pore is formed inside an Escherichia coli ClC exchanger. The secondary water pore is bifurcated from the chloride ion pathway at E148. From the systematic simulations, we determined that the glutamate residue exposed to the intracellular solution, E203, plays an important role as a trigger for the formation of the secondary water pore, and that the highly conserved tyrosine residue Y445 functions as a barrier that separates the proton from the chloride ion pathways. Based on our simulation results, we conclude that protons in the ClC exchanger are conducted via a water network through the secondary water pore, and we propose a new mechanism for the coupled transport of chloride ions and protons. It has been reported that several members of ClC proteins are not just channels that simply transport chloride ions across lipid bilayers; rather, they are exchangers that transport both the chloride ion and proton in opposite directions. However, the ion transit pathways and the mechanism of the coupled movement of these two ions have not yet been unveiled. In this article, we report a new finding (to our knowledge) of a water pore inside a prokaryotic ClC protein as revealed by computer simulation. This water pore is bifurcated from the putative chloride ion, and water molecules inside the new pore connect two glutamate residues that are known to be key residues for proton transport. On the basis of our simulation results, we conclude that the water wire that is formed inside the newly found pore acts as a proton pathway, which enables us to resolve many problems that could not be addressed by previous experimental studies.     

Proton channel models: Filling the gap between experimental data and the structural rationale

Voltage-gated proton channels are integral membrane proteins with the capacity to permeate elementary particles in a voltage and pH dependent manner. These proteins have been found in several species and are involved in various physiological processes. Although their primary topology is known, lack of details regarding their structures in the open conformation has limited analyses toward a deeper understanding of the molecular determinants of their function and regulation. Consequently, the function-structure relationships have been inferred based on homology models. In the present work, we review the existing proton channel models, their assumptions, predictions and the experimental facts that support them. Modeling proton channels is not a trivial task due to the lack of a close homolog template. Hence, there are important differences between published models. This work attempts to critically review existing proton channel models toward the aim of contributing to a better understanding of the structural features of these proteins.     

Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study

The kinetic Monte Carlo reaction-path-following technique is applied to determine the lowest-energy water pathway and the coordinating amino acids in bAQP1 and GlpF channels, both treated as rigid. In bAQP1, water molecules pass through the pore between the asparagine-proline-alanine (NPA) and selectivity filter (SF) sites one at a time. The water chain is interrupted at the SF where one water forms three stable hydrogen bonds with protein atoms. In this SF, water's conformation depends on the protonation locus of H182. In GlpF, two water molecules bond simultaneously to the NPA asparagines and pass through the SF in zigzag fashion. No water single-file forms in rigid GlpF. To accommodate a single file of waters requires narrowing the GlpF pore. Our results reveal that in both proteins a proposed bipolar water arrangement is thermally disrupted in the NPA region, especially in the cytoplasmic part of the pore. The equilibrium hydrogen-bonded chain is occasionally interrupted in the hydrophobic zones adjacent to the NPA motifs. The permeation of alkali cations through bAQP1 and GlpF is barred due to a large free-energy barrier in the NPA region as well as a large energy barrier blocking entry from the cytoplasm. Permeation of halides is prevented due to two large energy barriers in the cytoplasmic and periplasmic pores as well as a large free-energy barrier barring entry from the periplasm. Our results, based on modeling charge permeation, support an electrostatic rather than orientational basis for proton exclusion. Binding within the aquaporin pore cannot compensate sufficiently for dehydration of the protonic charge; there is also an electrostatic barrier in the NPA region blocking proton transport. The highly ordered single file of waters, which is drastically interrupted at the SF of bAQP1, may also contribute to proton block.     

Membrane hydration and structure on a subnanometer scale as seen by high resolution solid state nuclear magnetic resonance: POPC and POPC/C12EO4 model membranes.

The position on a subnanometer scale and the dynamics of structurally important water in model membranes was determined using a combination of proton magic-angle spinning NMR (MAS) with two-dimensional NOESY NMR techniques. Here, we report studies on phosphocholine lipid bilayers that were then modified by the addition of a nonionic surfactant that is shown to dehydrate the lipid. These studies are supplemented by 13C magic-angle spinning NMR investigations to get information on the dynamics of segmental motions of the membrane molecules. It can be shown that the hydrophilic chain of the surfactant is positioned at least partially within the hydrophobic core of the lipid bilayer. With the above NMR approach, we are able to establish molecular contacts between water and the lipid headgroup as well as with certain groups of the hydrocarbon chains and the glycerol backbone. This is possible because high resolution proton and 13C-NMR spectra of multilamellar bilayer membranes are obtained using MAS. A phase-sensitive NOESY must also be applied to distinguish positive and negative cross-peaks in the two-dimensional plot. These studies have high potential to investigate membrane proteins hydration and structural organization in a natural lipid bilayer surrounding.     

Protein structural change at the cytoplasmic surface as the cause of cooperativity in the bacteriorhodopsin photocycle.

The effects of excitation light intensity on the kinetics of the bacteriorhodopsin photocycle were investigated. The earlier reported intensity-dependent changes at 410 and 570 nm are explained by parallel increases in two of the rate constants, for proton transfers to D96 from the Schiff base and from the cytoplasmic surface, without changes in the others, as the photoexcited fraction is increased. Thus, it appears that the pKa of D96 is raised by a cooperative effect within the purple membrane. This interpretation of the wild-type kinetics was confirmed by results with several mutant proteins, where the rates are well separated in time and a model-dependent analysis is unnecessary. Based on earlier results that demonstrated a structural change of the protein after deprotonation of the Schiff base that increases the area of the cytoplasmic surface, and the effects of high hydrostatic pressure and lowered water activity on the photocycle steps in question, we suggest that the pKa of D96 is raised by a lateral pressure that develops when other bacteriorhodopsin molecules are photoexcited within the two-dimensional lattice of the purple membrane. Expulsion of no more than a few water molecules bound near D96 by this pressure would account for the calculated increase of 0.6 units in the pKa.