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.     

The mechanism of proton exclusion in the aquaporin-1 water channel

Aquaporins are efficient, yet strictly selective water channels. Remarkably, proton permeation is fully blocked, in contrast to most other water-filled pores which are known to conduct protons well. Blocking of protons by aquaporins is essential to maintain the electrochemical gradient across cellular and subcellular membranes. We studied the mechanism of proton exclusion in aquaporin-1 by multiple non-equilibrium molecular dynamics simulations that also allow proton transfer reactions. From the simulations, an effective free energy profile for the proton motion along the channel was determined with a maximum-likelihood approach. The results indicate that the main barrier is not, as had previously been speculated, caused by the interruption of the hydrogen-bonded water chain, but rather by an electrostatic field centered around the fingerprint Asn-Pro-Ala (NPA) motif. Hydrogen bond interruption only forms a secondary barrier located at the ar/R constriction region. The calculated main barrier height of 25-30 kJ mol(-1) matches the barrier height for the passage of protons across pure lipid bilayers and, therefore, suffices to prevent major leakage of protons through aquaporins. Conventional molecular dynamics simulations additionally showed that negatively charged hydroxide ions are prevented from being trapped within the NPA region by two adjacent electrostatic barriers of opposite polarity.     

Cation-selective mutations in the M2 domain of the inhibitory glycine receptor channel reveal determinants of ion-charge selectivity

Ligand-gated ion channel receptors mediate neuronal inhibition or excitation depending on their ion charge selectivity. An investigation into the determinants of ion charge selectivity of the anion-selective alpha1 homomeric glycine receptor (alpha1 glycine receptor [GlyR]) was undertaken using point mutations to residues lining the extra- and intracellular ends of the ion channel. Five mutant GlyRs were studied. A single substitution at the intracellular mouth of the channel (A-1'E GlyR) was sufficient to convert the channels to select cations over anions with P(Cl)/P(Na) = 0.34. This result delimits the selectivity filter and provides evidence that electrostatic interactions between permeating ions and pore residues are a critical factor in ion charge selectivity. The P-2'Delta mutant GlyR retained its anion selectivity (P(Cl)/P(Na) = 3.81), but it was much reduced compared with the wild-type (WT) GlyR (P(Cl)/P(Na) = 27.9). When the A-1'E and the P-2'Delta mutations were combined (selectivity double mutant [SDM] GlyR), the relative cation permeability was enhanced (P(Cl)/P(Na) = 0.13). The SDM GlyR was also Ca(2+) permeable (P(Ca)/P(Na) = 0.29). Neutralizing the extracellular mouth of the SDM GlyR ion channel (SDM+R19'A GlyR) produced a more Ca(2+)-permeable channel (P(Ca)/P(Na) = 0.73), without drastically altering monovalent charge selectivity (P(Cl)/P(Na) = 0.23). The SDM+R19'E GlyR, which introduces a negatively charged ring at the extracellular mouth of the channel, further enhanced Ca(2+) permeability (P(Ca)/P(Na) = 0.92), with little effect on monovalent selectivity (P(Cl)/P(Na) = 0.19). Estimates of the minimum pore diameter of the A-1'E, SDM, SDM+R19'A, and SDM+R19'E GlyRs revealed that these pores are larger than the alpha1 GlyR, with the SDM-based GlyRs being comparable in diameter to the cation-selective nicotinic acetylcholine receptors. This result provides evidence that the diameter of the ion channel is also an important factor in ion charge selectivity.     

Peregrination of the selectivity filter delineates the pore of the human voltage-gated proton channel hHV1

Extraordinary selectivity is crucial to all proton-conducting molecules, including the human voltage-gated proton channel (hH_V1), because the proton concentration is >10⁶ times lower than that of other cations. Here we use “selectivity filter scanning” to elucidate the molecular requirements for proton-specific conduction in hH_V1. Asp¹¹², in the middle of the S1 transmembrane helix, is an essential part of the selectivity filter in wild-type (WT) channels. After neutralizing Asp¹¹² by mutating it to Ala (D112A), we introduced Asp at each position along S1 from 108 to 118, searching for “second site suppressor” activity. Surprisingly, most mutants lacked even the anion conduction exhibited by D112A. Proton-specific conduction was restored only with Asp or Glu at position 116. The D112V/V116D channel strikingly resembled WT in selectivity, kinetics, and ΔpH-dependent gating. The S4 segment of this mutant has similar accessibility to WT in open channels, because R211H/D112V/V116D was inhibited by internally applied Zn²⁺. Asp at position 109 allowed anion permeation in combination with D112A but did not rescue function in the nonconducting D112V mutant, indicating that selectivity is established externally to the constriction at F150. The three positions that permitted conduction all line the pore in our homology model, clearly delineating the conduction pathway. Evidently, a carboxyl group must face the pore directly to enable conduction. Molecular dynamics simulations indicate reorganization of hydrogen bond networks in the external vestibule in D112V/V116D. At both positions where it produces proton selectivity, Asp frequently engages in salt linkage with one or more Arg residues from S4. Surprisingly, mean hydration profiles were similar in proton-selective, anion-permeable, and nonconducting constructs. That the selectivity filter functions in a new location helps to define local environmental features required to produce proton-selective conduction.     

Single amino acids in the carboxyl terminal domain of aquaporin-1 contribute to cGMP-dependent ion channel activation

Background  

Aquaporin-1 (AQP1) functions as an osmotic water channel and a gated cation channel. Activation of the AQP1 ion conductance by intracellular cGMP was hypothesized to involve the carboxyl (C-) terminus, based on amino acid sequence alignments with cyclic-nucleotide-gated channels and cGMP-selective phosphodiesterases.     

Origins of peptide selectivity and phosphoinositide binding revealed by structures of disabled-1 PTB domain complexes

Formation of the mammalian six-layered neocortex depends on a signaling pathway that involves Reelin, the very low-density lipoprotein receptor, the apolipoprotein E receptor-2 (ApoER2), and the adaptor protein Disabled-1 (Dab1). The 1.5 A crystal structure of a complex between the Dab1 phosphotyrosine binding (PTB) domain and a 14-residue peptide from the ApoER2 tail explains the unusual preference of Dab1 for unphosphorylated tyrosine within the NPxY motif of the peptide. Crystals of the complex soaked with the phosphoinositide PI-4,5P(2) (PI) show that PI binds to conserved basic residues on the PTB domain opposite the peptide binding groove. This finding explains how the Dab1 PTB domain can simultaneously bind PI and the ApoER2 tail. Recruitment of the Dab1 PTB domain to PI-rich regions of the plasma membrane may facilitate association with the Reelin receptor cytoplasmic tails to transduce a critical positional cue to migrating neurons.     

Anomalous proton selectivity in a large channel: colicin A

Some of the bactericidal proteins known as colicins exert their toxic action by forming a large, nonselective channel in the inner membrane of target bacteria. The structure of this channel is unknown. It conducts large ions but has a much smaller conductance than would be expected for a channel of its deduced size. Here we report that the colicin channel, particularly the colicin A channel, is selective for protons over other cations (and anions) by many orders of magnitude. This was deduced from measurements of reversal potentials in pH gradients across planar lipid bilayers containing these channels. For example, in symmetric 0.1 M KCl with a pH 5/pH 8 gradient across the membrane, the reversal potential of colicin A is -21 mV, rather than 0. Such a result would be unremarkable for a narrow channel but is beyond explanation by current understanding of permeation for a channel of its diameter. For this reason, we re-examined the issue of the diameter of the channel lumen and confirmed that the lumen is indeed "too large" ( approximately 10 A) to select for protons by the amount that we measure. We are thus compelled to propose that an unorthodox mechanism is at work in this protein.     

Charge delocalization in proton channels, II: the synthetic LS2 channel and proton selectivity

In this study, the minimalist synthetic LS2 channel is used as a prototype to examine the selectivity of protons over other cations. The free-energy profiles along the transport pathway of LS2 are calculated for three cation species: a realistic delocalized proton (including Grotthuss shuttling)--H(+), a classical (nonshuttling) hydronium--H(3)O(+), and a potassium cation--K(+). The overall barrier for K(+) is approximately twice as large as that for H(+), explaining the >100 times larger maximal ion conductance for the latter, in qualitative agreement with the experimental result. The profile for the classical hydronium is quantitatively intermediate between those of H(+) and K(+) and qualitatively more similar to that of H(+), for which the locations of the peaks are well correlated with the troughs of the pore radius profile. There is a strong correlation between the free-energy profiles and the very different characteristic hydration structures of the three cation species. This work suggests that the passage of various cations through ion channels cannot always be explained by simple electrostatic desolvation considerations.     

The transmembrane domain and the proton channel in proton-pumping transhydrogenases

Proton-pumping nicotinamide nucleotide transhydrogenases are composed of three main domains, the NAD(H)-binding and NADP(H)-binding hydrophilic domains I (dI) and III (dIII), respectively, and the hydrophobic domain II (dII) containing the assumed proton channel. dII in the Escherichia coli enzyme has recently been characterised with regard to topology and a packing model of the helix bundle in dII is proposed. Extensive mutagenesis of conserved charged residues of this domain showed that important residues are betaHis91 and betaAsn222. The pH dependence of betaH91D, as well as betaH91C (unpublished), when compared to that of wild type shows that reduction of 3-acetylpyridine-NAD(+) by NADPH, i.e., the reverse reaction, is optimal at a pH essentially coinciding with the pK(a) of the residue in the beta91 position. It is therefore concluded that the wild-type transhydrogenase is regulated by the degree of protonation of betaHis91. The mechanisms of the interactions between dI+dIII and dII are suggested to involve pronounced conformational changes in a 'hinge' region around betaR265.     

Computational studies of proton transport through the M2 channel

The M2 ion channel is an essential component of the influenza A virus. This low-pH gated channel has a high selectivity for protons. Evidence from various experimental data has indicated that the essential structure responsible for the channel is a parallel homo-tetrameric alpha-helix bundle having a left-handed twist with each helix tilted with respect to the membrane normal. A backbone structure has been determined by solid state nuclear magnetic resonance (NMR). Though detailed structures for the side chains are not available yet, evidence has indicated that His37 and Trp41 in the alpha-helix are implicated in the local molecular structure responsible for the selectivity and channel gate. More definitive conformations for the two residues were recently suggested based on the known backbone structure and recently obtained NMR data. While two competitive proton-conductance mechanisms have been proposed, the actual proton-conductance mechanism remains an unsolved problem. Computer simulations of an excess proton in the channel and computational studies of the His37/Trp41 conformations have provided insights into these structural and mechanism issues.     

Backbone structure of the amantadine-blocked trans-membrane domain M2 proton channel from Influenza A virus

Amantadine is known to block the M2 proton channel of the Influenza A virus. Here, we present a structure of the M2 trans-membrane domain blocked with amantadine, built using orientational constraints obtained from solid-state NMR polarization-inversion-spin-exchange-at-the-magic-angle experiments. The data indicates a kink in the monomer between two helical fragments having 20 degrees and 31 degrees tilt angles with respect to the membrane normal. This monomer structure is then used to construct a plausible model of the tetrameric amantadine-blocked M2 trans-membrane channel. The influence of amantadine binding through comparative cross polarization magic-angle spinning spectra was also observed. In addition, spectra are shown of the amantadine-resistant mutant, S31N, in the presence and absence of amantadine.     

Interaction of divalent metal ions with the carboxyl-terminal domain of human voltage-gated proton channel Hv1

The voltage-gated proton channel Hv1 functions as a dimer, in which the intracellular C-terminal domain of the protein is responsible for the dimeric architecture and regulates proton permeability. Although it is well known that divalent metal ions have effect on the proton channel activity, the interaction of divalent metal ions with the channel in detail is not well elucidated. Herein, we investigated the interaction of divalent metal ions with the C-terminal domain of human Hv1 by CD spectra and fluorescence spectroscopy. The divalent metal ions binding induced an obvious conformational change at pH 7 and a pH-sensitive reduction of thermostability in the C-terminal domain. The interactions were further estimated by fluorescence spectroscopy experiments. There are at least two binding sites for divalent metal ions binding to the C-terminal domain of Hv1, either of which is close to His²⁴⁴ or His²⁶⁶ residue. The binding of Zn²⁺ to the two sites both enhanced the fluorescence of the protein at pH 7, whereas the binding of other divalent metal ions to the two sites all resulted fluorescence quenching. The orders of the strength of divalent metal ions binding to the two sites from strong to weak are both Co²⁺, Ca²⁺, Ni²⁺, Mg²⁺, and Mn²⁺. The strength of Ca²⁺, Co²⁺, Mg²⁺, Mn²⁺ and Ni²⁺ binding to the site close to His²⁴⁴ is stronger than that of these divalent metal ions binding to the site close to His²⁶⁶.     

Ion selectivity of the Kat1 K+ channel pore

Kat1 is a highly selective inward-rectifying K⁺ channel that opens for extended periods under conditions of extreme hyperpolarization. Over 200 point mutants in the pore region of the Kat1 K⁺ channel were generated and examined in the yeast Saccharomyces cerevisiae and Xenopus oocytes to assess the effect of the mutations on ion selectivity. Substitutions at the tyrosine of the signature sequence G-Y-G resulted in the most significant alterations in ion selectivity, consistent with its role in the selectivity filter. However, greater than 80% of the mutations throughout the greater pore region also conferred a defect in selectivity demonstrating that the entire pore of Kat1 contributes to the ion selectivity of this channel. Surprisingly, we identified a novel class of mutant channel that conferred enhanced selectivity of K⁺ over Na⁺. Mutants of this class frequently displayed sensitivity to the competing ion Cs⁺. This finding has led us to speculate that the Kat1 channel pore has evolved to balance not only K⁺/Na⁺ selectivity, but selectivity over Cs⁺, and possibly a wide spectrum of potential competing ions.     

Kinetics of proton transport into influenza virions by the viral M2 channel

M2 protein of influenza A viruses is a tetrameric transmembrane proton channel, which has essential functions both early and late in the virus infectious cycle. Previous studies of proton transport by M2 have been limited to measurements outside the context of the virus particle. We have developed an in vitro fluorescence-based assay to monitor internal acidification of individual virions triggered to undergo membrane fusion. We show that rimantadine, an inhibitor of M2 proton conductance, blocks the acidification-dependent dissipation of fluorescence from a pH-sensitive virus-content probe. Fusion-pore formation usually follows internal acidification but does not require it. The rate of internal virion acidification increases with external proton concentration and saturates with a pK(m) of ～4.7. The rate of proton transport through a single, fully protonated M2 channel is approximately 100 to 400 protons per second. The saturating proton-concentration dependence and the low rate of internal virion acidification derived from authentic virions support a transporter model for the mechanism of proton transfer.     

Hypertonic induction of aquaporin-1 water channel independent of transcellular osmotic gradient

Aquaporin-1 (AQP1) water channel plays a critical role for water reabsorption in the urinary concentrating mechanism. AQP1 expression in renal cells is upregulated by hypertonicity, but not urea, suggesting the requirement of an osmotic gradient. To investigate whether AQP1 expression is regulated by apical and/or basolateral hypertonicity, murine renal medullary mIMCD-K2 cells grown on permeable support were exposed to hypertonic medium. When the medium on the apical or basolateral membrane side was switched to hypertonic, the transcellular osmotic gradient was dissipated within 8h. Basolateral hypertonicity increased AQP1 expression more than apical hypertonicity. Comparable apical and basolateral hypertonicity without a transcellular hypertonic gradient, however, increased AQP1 expression. Cell surface biotinylation experiments revealed that hypertonicity promoted AQP1 trafficking to both plasma cell membranes. These results indicate that AQP1 expression is predominantly mediated by basolateral hypertonicity but a transcellular osmotic gradient is not necessary for its induction.     

External Ba2+ block of the two-pore domain potassium channel TREK-1 defines conformational transition in its selectivity filter

TREK-1 is a member of the two-pore domain potassium channel family that is known as a leak channel and plays a key role in many physiological and pathological processes. The conformational transition of the selectivity filter is considered as an effective strategy for potassium channels to control the course of potassium efflux. It is well known that TREK-1 is regulated by a large volume of extracellular and intracellular signals. However, until now, little was known about the selectivity filter gating mechanism of the channel. In this research, it was found that Ba(2+) blocked the TREK-1 channel in a concentration- and time-dependent manner. A mutagenesis analysis showed that overlapped binding of Ba(2+) at the assumed K(+) binding site 4 (S4) within the selectivity filter was responsible for the inhibitory effects on TREK-1. Then, Ba(2+) was used as a probe to explore the conformational transition in the selectivity filter of the channel. It was confirmed that collapsed conformations were induced by extracellular K(+)-free and acidification at the selectivity filters, leading to nonconductive to permeable ions. Further detailed characterization demonstrated that the two conformations presented different properties. Additionally, the N-terminal truncated isoform (ΔN41), a product derived from alternative translation initiation, was identified as a constitutively nonconductive variant. Together, these results illustrate the important role of selectivity filter gating in the regulation of TREK-1 by the extracellular K(+) and proton.     

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.     

Proton transport behavior through the influenza A M2 channel: insights from molecular simulation

The structural properties of the influenza A virus M2 transmembrane channel in dimyristoylphosphatidylcholine bilayer for each of the four protonation states of the proton-gating His-37 tetrad and their effects on proton transport for this low-pH activated, highly proton-selective channel are studied by classical molecular dynamics with the multistate empirical valence-bond (MS-EVB) methodology. The excess proton permeation free energy profile and maximum ion conductance calculated from the MS-EVB simulation data combined with the Poisson-Nernst-Planck theory indicates that the triply protonated His-37 state is the most likely open state via a significant side-chain conformational change of the His-37 tetrad. This proposed open state of M2 has a calculated proton permeation free energy barrier of 7 kcal/mol and a maximum conductance of 53 pS compared to the experimental value of 6 pS. By contrast, the maximum conductance for Na(+) is calculated to be four orders of magnitude lower, in reasonable agreement with the experimentally observed proton selectivity. The pH value to activate the channel opening is estimated to be 5.5 from dielectric continuum theory, which is also consistent with experimental results. This study further reveals that the Ala-29 residue region is the primary binding site for the antiflu drug amantadine (AMT), probably because that domain is relatively spacious and hydrophobic. The presence of AMT is calculated to reduce the proton conductance by 99.8% due to a significant dehydration penalty of the excess proton in the vicinity of the channel-bound AMT.