CHAP: A Versatile Tool for the Structural and Functional Annotation of Ion Channel Pores

The control of ion channel permeation requires the modulation of energetic barriers or “gates” within their pores. However, such barriers are often simply identified from the physical dimensions of the pore. Such approaches have worked well in the past, but there is now evidence that the unusual behavior of water within narrow hydrophobic pores can produce an energetic barrier to permeation without requiring steric occlusion of the pathway. Many different ion channels have now been shown to exploit “hydrophobic gating” to regulate ion flow, and it is clear that new tools are required for more accurate functional annotation of the increasing number of ion channel structures becoming available. We have previously shown how molecular dynamics simulations of water can be used as a proxy to predict hydrophobic gates, and we now present a new and highly versatile computational tool, the Channel Annotation Package (CHAP) that implements this methodology.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

Molecular dynamics study of the KcsA potassium channel

The structural, dynamical, and thermodynamic properties of a model potassium channel are studied using molecular dynamics simulations. We use the recently unveiled protein structure for the KcsA potassium channel from Streptomyces lividans. Total and free energy profiles of potassium and sodium ions reveal a considerable preference for the larger potassium ions. The selectivity of the channel arises from its ability to completely solvate the potassium ions, but not the smaller sodium ions. Self-diffusion of water within the narrow selectivity filter is found to be reduced by an order of magnitude from bulk levels, whereas the wider hydrophobic section of the pore maintains near-bulk self-diffusion. Simulations examining multiple ion configurations suggest a two-ion channel. Ion diffusion is found to be reduced to approximately (1)/(3) of bulk diffusion within the selectivity filter. The reduced ion mobility does not hinder the passage of ions, as permeation appears to be driven by Coulomb repulsion within this multiple ion channel.     

Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.     

Homology modeling and molecular dynamics simulations of transmembrane domain structure of human neuronal nicotinic acetylcholine receptor

A three-dimensional model of the transmembrane domain of a neuronal-type nicotinic acetylcholine receptor (nAChR), (alpha4)2(beta2)3, was constructed from a homology structure of the muscle-type nAChR recently determined by cryo-electron microscopy. The neuronal channel model was embedded in a fully hydrated DMPC lipid bilayer, and molecular-dynamics simulations were performed for 5 ns. A comparative analysis of the neuronal- versus muscle-type nAChR models revealed many conserved pore-lining residues, but an important difference was found near the periplasmic mouth of the pore. A flickering salt-bridge of alpha4-E266 with its adjacent beta2-K260 was observed in the neuronal-type channel during the course of the molecular-dynamics simulations. The narrowest region, with a pore radius of approximately 2 A formed by the salt-bridges, does not seem to be the restriction site for a continuous water passage. Instead, two hydrophobic rings, formed by alpha4-V259, alpha4-L263, and the homologous residues in the beta2-subunits, act as the gates for water flow, even though the region has a slightly larger pore radius. The model offers new insight into the water transport across the (alpha4)2(beta2)3 nAChR channel, and may lead to a better understanding of the structures, dynamics, and functions of this family of ion channels.     

Probing structure-function relationships and gating mechanisms in the CorA Mg2+ transport system

Recent crystal structures of the CorA Mg(2+) transport protein from Thermotoga maritima (TmCorA) revealed an unusually long ion pore putatively gated by hydrophobic residues near the intracellular end and by universally conserved asparagine residues at the periplasmic entrance. A conformational change observed in an isolated funnel domain structure also led to a proposal for the structural basis of gating. Because understanding the molecular mechanisms underlying ion channel and transporter gating remains an important challenge, we have undertaken a structure-guided engineering approach to probe structure-function relationships in TmCorA. The intracellular funnel domain is shown to constitute an allosteric regulatory module that can be engineered to promote an activated or closed state. A periplasmic gate centered about a proline-induced kink of the pore-lining helix is described where "helix-straightening" mutations produce a dramatic gain-of-function. Mutation to the narrowest constriction along the pore demonstrates that a hydrophobic gate is operational within this Mg(2+)-selective transport protein and likely forms an energetic barrier to ion flux. We also provide evidence that highly conserved acidic residues found in the short periplasmic loop are not essential for TmCorA function or Mg(2+) selectivity but may be required for proper protein folding and stability. This work extends our gating model for the CorA-Alr1-Mrs2 superfamily and reveals features that are characteristic of an ion channel. Aspects of these results that have broader implications for a range of channel and transporter families are highlighted.     

Aggregation and porin-like channel activity of a beta sheet peptide

Beta sheet peptides (e.g., amyloid beta) are known to form ion channels in lipid bilayers possibly through aggregation, though the channel structure is not clear. We have recently reported that a short beta sheet peptide, (xSxG)(6), forms porin-like voltage-gated channels in lipid bilayers [Thundimadathil et al. (2005) Biochem. Biophys. Res. Commun. 330, 585-590]. To account for the porin-like activity, oligomerization of the peptide into a beta barrel-like structure was proposed. In this work, peptide aggregation in aqueous and membrane environments and a detailed study of channel properties were performed to gain insight into the mechanism of channel formation. The complex nature of the channel was revealed by kinetic analysis and the occurrence of interconverting multiple conductance states. Ion channels were inhibited by Congo red, suggesting that the peptide aggregates are the active channel species. Peptide aggregation and fibril formation in water were confirmed by electron microscopy (EM) and Congo red binding studies. Furthermore, oligomeric structures in association with lipid bilayers were detected. Circular dichroism of peptide-incorporated liposomes and peptide-lipid binding studies using EM suggest a lipid-induced beta sheet aggregation. Gel electrophoresis of peptide-incorporated liposomes showed dimeric and multimeric structures. Taken together, this work indicates insertion of (xSxG)(6) as oligomers into the lipid bilayer, followed by rearrangement into a beta barrel-like pore structure. A large peptide pore comprising several individual beta sheets or smaller beta sheet aggregates is expected to have a complex behavior in membranes. A dyad repeat sequence and the presence of glycine, serine, and hydrophobic residues in a repeated pattern in this peptide may be providing a favorable condition for the formation of a beta barrel-like structure in lipid bilayers.     

Accelerated diffusion of Na(+) in a hydrophobic region revealed by molecular dynamics simulations of a synthetic ion channel

To get insight into the significance of the hydrophobic lining on the ion permeation, we performed molecular dynamics simulations on a Na(+) permeation through a de novo synthetic hydrophobic channel. Electrophysiological study has suggested that the channel is formed from a tail-to-tail associated dimer of a cyclic octa-peptide coupled with hydrophobic acyl chains. The acyl chains line the channel pore while the cyclic peptide forms the channel entrance [Z. Qi, M. Sokabe, K. Donowaki, H. Ishida, Biophys. J. 76 (1999) 631]. Molecular dynamics simulation of water in the channel indicated that the inferred structure is physically reasonable [Z. Qi, M. Sokabe, Biophys. Chem. 71 (1998) 35]. In the present study, the potential energy profile of the Na(+) and the energy contributions from each component of the system at different positions along the channel axis were calculated. An energy well instead of a peak is located at the central hydrophobic cavity of the channel, due to its ability of accommodating at least five water molecules to hydrate the ion. Interestingly, the ion diffuses much faster in the hydrophobic acyl chain region, particularly in the central hydrophobic cavity, than it does in the peptide ring region and even surprisingly faster than that in the bulk phase. These results provide a physical basis for an idea that the hydrophobic lining of the K(+) channel [D.A. Doyle, J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, R. MacKinnon, Science 280 (1998) 69] plays an active role to facilitate the ion permeation through the channel pore.     

Wide nanoscopic pore of maxi-anion channel suits its function as an ATP-conductive pathway

The newly proposed function of the maxi-anion channel as a conductive pathway for ATP release requires that its pore is sufficiently large to permit passage of a bulky ATP(4-) anion. We found a linear relationship between relative permeability of organic anions of different size and their relative ionic mobility (measured as the ratio of ionic conductance) with a slope close to 1, suggesting that organic anions tested with radii up to 0.49 nm (lactobionate) move inside the channel by free diffusion. In the second approach, we, for the first time, succeeded in pore sizing by the nonelectrolyte exclusion method in single-channel patch-clamp experiments. The cutoff radii of PEG molecules that could access the channel from intracellular (1.16 nm) and extracellular (1.42 nm) sides indicated an asymmetry of the two entrances to the channel pore. Measurements by symmetrical two-sided application of PEG molecules yielded an average functional pore radius of approximately 1.3 nm. These three estimates are considerably larger than the radius of ATP(4-) (0.57-0.65 nm) and MgATP(2-) (approximately 0.60 nm). We therefore conclude that the nanoscopic maxi-anion channel pore provides sufficient room to accommodate ATP and is well suited to its function as a conductive pathway for ATP release in cell-to-cell communication.     

Probing structural features of Alzheimer's amyloid-β pores in bilayers using site-specific amino acid substitutions

A current hypothesis for the pathology of Alzheimer's disease (AD) proposes that amyloid-β (Aβ) peptides induce uncontrolled, neurotoxic ion flux across cellular membranes. The mechanism of ion flux is not fully understood because no experiment-based Aβ channel structures at atomic resolution are currently available (only a few polymorphic states have been predicted by computational models). Structural models and experimental evidence lend support to the view that the Aβ channel is an assembly of loosely associated mobile β-sheet subunits. Here, using planar lipid bilayers and molecular dynamics (MD) simulations, we show that amino acid substitutions can be used to infer which residues are essential for channel structure. We created two Aβ(1-42) peptides with point mutations: F19P and F20C. The substitution of Phe19 with Pro inhibited channel conductance. MD simulation suggests a collapsed pore of F19P channels at the lower bilayer leaflet. The kinks at the Pro residues in the pore-lining β-strands induce blockage of the solvated pore by the N-termini of the chains. The cysteine mutant is capable of forming channels, and the conductance behavior of F20C channels is similar to that of the wild type. Overall, the mutational analysis of the channel activity performed in this work tests the proposition that the channels consist of a β-sheet rich organization, with the charged/polar central strand containing the mutation sites lining the pore, and the C-terminal strands facing the hydrophobic lipid tails. A detailed understanding of channel formation and its structure should aid studies of drug design aiming to control unregulated Aβ-dependent ion fluxes.     

In channelrhodopsin-2 Glu-90 is crucial for ion selectivity and is deprotonated during the photocycle

The light-activated microbial ion channel channelrhodopsin-2 (ChR2) is a powerful tool to study cellular processes with high spatiotemporal resolution in the emerging field of optogenetics. To customize the channel properties for optogenetic experiments, a detailed understanding of its molecular reaction mechanism is essential. Here, Glu-90, a key residue involved in the gating and selectivity mechanism of the ion channel is characterized in detail. The deprotonation of Glu-90 during the photocycle is elucidated by time-resolved FTIR spectroscopy, which seems to be part of the opening mechanism of the conductive pore. Furthermore, Glu-90 is crucial to ion selectivity as also revealed by mutation of this residue combined with voltage clamp experiments. By dynamic homology modeling, we further hypothesized that the conductive pore is flanked by Glu-90 and located between helices A, B, C, and G.     

Structural basis of ligand activation in a cyclic nucleotide regulated potassium channel

Here we describe the initial functional characterization of a cyclic nucleotide regulated ion channel from the bacterium Mesorhizobium loti and present two structures of its cyclic nucleotide binding domain, with and without cAMP. The domains are organized as dimers with the interface formed by the linker regions that connect the nucleotide binding pocket to the pore domain. Together, structural and functional data suggest the domains form two dimers on the cytoplasmic face of the channel. We propose a model for gating in which ligand binding alters the structural relationship within a dimer, directly affecting the position of the adjacent transmembrane helices.     

Roles of the Hydrophobic Gate and Exit Channel in Vigna radiata Pyrophosphatase Ion Translocation

Membrane-embedded pyrophosphatase (M-PPase) hydrolyzes pyrophosphate to drive ion (H⁺ and/or Na⁺) translocation. We determined crystal structures and functions of Vigna radiata M-PPase (VrH⁺-PPase), the VrH⁺-PPase–2P_i complex and mutants at hydrophobic gate (residue L555) and exit channel (residues T228 and E225). Ion pore diameters along the translocation pathway of three VrH⁺-PPases complexes (P_i-, 2P_i- and imidodiphosphate-bound states) present a unique wave-like profile, with different pore diameters at the hydrophobic gate and exit channel, indicating that the ligands induced pore size alterations. The 2P_i-bound state with the largest pore diameter might mimic the hydrophobic gate open. In mutant structures, ordered waters detected at the hydrophobic gate among VrH⁺-PPase imply the possibility of solvation, and numerous waters at the exit channel might signify an open channel. A salt-bridge, E225–R562 is at the way out of the exit channel of VrH⁺-PPase; E225A mutant makes the interaction eliminated and reveals a decreased pumping ability. E225–R562 might act as a latch to regulate proton release. A water wire from the ion gate (R-D-K-E) through the hydrophobic gate and into the exit channel may reflect the path of proton transfer.