Modeling ion channel regulation

Remarkable recent successes in structure determinations of voltage-gated channels, ligand-gated channels, mechanosensitive channels and proton channels have advanced our understanding of the molecular basis of ion channel gating substantially. Models have helped to clarify aspects of this process and are now being designed as sophisticated biomimetics for various technological applications.     

Theoretical models of the ion channel structure of amyloid beta-protein.

Theoretical methods are used to develop models for the ion channel structure of the membrane-bound amyloid beta-protein. This follows recent observations that the beta-protein forms cation-selective channels in lipid bilayers in vitro. Amyloid beta-protein is the main component of the extracellular plaques in the brain that are characteristic of Alzheimer's disease. Based on the amino acid sequence and the unique environment of the membrane, the secondary structure of the 40-residue beta-protein is predicted to form a beta-hairpin followed by a helix-turn-helix motif. The channel structures were-designed as aggregates of peptide subunits in identical conformations. Three types of models were developed that are distinguished by whether the pore is formed by the beta-hairpins, the middle helices, or by the more hydrophobic C-terminal helices. The latter two types can be converted back and forth by a simple conformational change, which would explain the variable conduction states observed for a single channel. It is also demonstrated how lipid headgroups could be incorporated into the pore lining, and thus affect the ion selectivity. The atomic-scale detail of the models make them useful for designing experiments to determine the real structure of the channel, and thus further the understanding of peptide channels in general. In addition, if beta-protein-induced channel activity is found to be the cause of cell death in Alzheimer's disease, then the models may be helpful in designing counteracting drugs.     

Modeling negative ion defect migration through the gramicidin A channel

The results of potential of mean force (PMF) calculations for the distinct stages of proton conduction through the gramicidin A channel, including proton migration, reorientation of the water file and negative ion defect migration, are presented. The negative ion defect migration mechanism was hypothesized in experimental studies but was not considered previously in molecular dynamics simulations. The model system consisted of the peptide chains constructed on the base of the structure PDBID:1JNO, the inner file of nine water molecules and external clusters of water molecules placed at both ends of the channel. Potential energy functions were computed with the CHARMM/PM6/TIP3P parameters. The results obtained for proton migration and water file reorientation are basically consistent with those reported previously by Pómès and Roux (Biophys J 82:2304, 2002) within the similar approach. For the newly considered mechanism of negative ion defect migration from the channel center to the end of the water file we obtain the energy 3.8 kcal mol⁻¹ which is not considerably different from the activation energy of water reorientation, 5.4 kcal mol⁻¹. Therefore this mechanism may principally compete for the rate-limiting step in proton conduction in gramicidin.     

Crumpled structure of the custom hydrophobic lytic peptide cecropin B3.

The solution structure of a custom lytic peptide, cecropin B3 (CB3), having two identical hydrophobic segments on both the N- and C-termini, was investigated by two-dimensional NMR spectroscopy. The need to determine the structure of this peptide is rooted in its specific ability to lyse lipid layers that have a high content of anionic lipid. The lytic activities of CB3 on cell membranes including cancer cells and bacteria is found to be less than cecropin B1. The results show that CB3 has four discrete segments forming alpha helical structures. The crumpled structure of CB3 provides evidence for the lysis of the lipid layer being via a pathway that differs from pore formation. The results in this study provide strong clues towards a rational design for a potent antimicrobial and antitumor peptide.     

Modeling of agonist binding to the ligand-gated ion channel superfamily of receptors

A generalized model is presented of agonist binding to ligand-gated ion channels (LGICs). Broad similarity in the structure of agonists suggests that the binding sites of LGICs may have evolved from a protobinding site. Aligned sequence data identified as a candidate for such a site a highly conserved 15 residue stretch of primary structure in the N-terminal extracellular region of all known LGIC subunits. We modeled this subregion, termed the cys-loop, as a rigid, amphiphilic beta-hairpin and propose that it may form a major determinant of a conserved structural binding cleft. In the model of the binding complex (1) an invariant aspartate residue at position 11 of the cys-loop is the anionic site interacting with the positively charged amine group of agonists, (2) a local dipole within the pi-electron system of agonists is favorably oriented in the electrostatic field of the invariant aspartate, (3) the epsilon ring-proton of a conserved aromatic residue at the turn of the cys-loop interacts orthogonally with the agonist pi-electron density at its electronegative center, and (4) selective recognition is partly a result of the type of amino acid residue at position 6 of the cys-loop. Additionally, formation of a hydrogen bond between the electronegative atom of the pi-electron system of agonist and a complementary group in the receptor may be important in the high-affinity binding of agonists.     

Open channel block by gadolinium ion of the stretch-inactivated ion channel in mdx myotubes.

Currents flowing through single stretch-inactivated ion channels were recorded from cell-attached patches on myotubes from mdx mice. Adding micromolar concentrations of gadolinium to patch electrodes containing normal saline produced rapid transitions in the single-channel current between the fully open and closed states. The kinetics of the current fluctuations followed the predictions of a simple model of open channel block in which the transitions in the current arise from the entry and exit of Gd from the channel pore: histograms of the open and closed times were well fit with single exponentials, the blocking rate depended linearly on the concentration of gadolinium in the patch electrode, and the unblocking rate was independent of the concentration of gadolinium. Hyperpolarizing the patch increased the rate of unblocking (approximately e-fold per 85 mV), suggesting the charged blocking particle can exit the channel into the cell under the influence of the applied membrane field. The rate of blocking was rapid and was independent of the patch potential, consistent with the rate of ion entry into the pore being determined by its rate of diffusion in solution. When channel open probability was reduced by applying suction to the electrode, the blocking kinetics were independent of the extent of inactivation, suggesting that mechanosensitive gating does not modify the structure of the channel pore.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Modeling of ion permeation in calcium and sodium channel selectivity filters

Structure-function studies have shown that it is possible to convert a sodium channel to a calcium-selective channel by a single amino acid substitution in the selectivity filter locus. Ion permeation through the "model selectivity filter" was modeled with a reduced set of functional groups representative of the constituent amino acid side chains. Force-field minimizations were conducted to obtain the energy profile of the cations as they get desolvated and bind to the "model selectivity filter." The calculations suggest that the ion selectivity in the calcium channel is due to preferential binding, whereas in the sodium channel it is due to exclusion. Energetics of displacement of a bound cation from the calcium "model selectivity filter" by another cation suggest that "multi-ion mechanism" reduces the activation barrier for ion permeation. Thus, the simple model captures qualitatively most of the conduction characteristics of sodium and calcium channels. However, the computed barriers for permeation are fairly large, suggesting that ion interaction with additional residues along the transport path may be essential to effect desolvation.     

Insight toward epithelial Na+ channel mechanism revealed by the acid-sensing ion channel 1 structure

The epithelial Na(+) channel/degenerin (ENaC/DEG) protein family includes a diverse group of ion channels, including nonvoltage-gated Na(+) channels of epithelia and neurons, and the acid-sensing ion channel 1 (ASIC1). In mammalian epithelia, ENaC helps regulate Na(+) and associated water transport, making it a critical determinant of systemic blood pressure and pulmonary mucosal fluidity. In the nervous system, ENaC/DEG proteins are related to sensory transduction. While the importance and physiological function of these ion channels are established, less is known about their structure. One hallmark of the ENaC/DEG channel family is that each channel subunit has only two transmembrane domains connected by an exceedingly large extracellular loop. This subunit structure was recently confirmed when Jasti and colleagues determined the crystal structure of chicken ASIC1, a neuronal acid-sensing ENaC/DEG channel. By mapping ENaC to the structural coordinates of cASIC1, as we do here, we hope to provide insight toward ENaC structure. ENaC, like ASIC1, appears to be a trimeric channel containing 1alpha, 1beta, and 1gamma subunit. Heterotrimeric ENaC and monomeric ENaC subunits within the trimer possibly contain many of the major secondary, tertiary, and quaternary features identified in cASIC1 with a few subtle but critical differences. These differences are expected to have profound effects on channel behavior. In particular, they may contribute to ENaC insensitivity to acid and to its constitutive activity in the absence of time- and ligand-dependent inactivation. Experiments resulting from this comparison of cASIC1 and ENaC may help clarify unresolved issues related to ENaC architecture, and may help identify secondary structures and residues critical to ENaC function.     

Effects of membrane lipids on ion channel structure and function

Biologic membranes are not simply inert physical barriers, but complex and dynamic environments that affect membrane protein structure and function. Residing within these environments, ion channels control the flux of ions across the membrane through conformational changes that allow transient ion flux through a central pore. These conformational changes may be modulated by changes in transmembrane electrochemical potential, the binding of small ligands or other proteins, or changes in the local lipid environment. Ion channels play fundamental roles in cellular function and, in higher eukaryotes, are the primary means of intercellular signaling, especially between excitable cells such as neurons. The focus of this review is to examine how the composition of the bilayer affects ion channel structure and function. This is an important consideration because the bilayer composition varies greatly in different cell types and in different organellar membranes. Even within a membrane, the lipid composition differs between the inner and outer leaflets, and the composition within a given leaflet is both heterogeneous and highly dynamic. Differential packing of lipids (and proteins) leads to the formation of microdomains, and lateral diffusion of these microdomains or "lipid rafts" serve as mobile platforms for the clustering and organization of bilayer constituents including ion channels. The structure and function of these channels are sensitive to specific chemical interactions with neighboring components of the membrane and also to the biophysical properties of their membrane microenvironment (e.g., fluidity, lateral pressure profile, and bilayer thickness). As specific examples, we have focused on the K+ ion channels and the ligand-gated nicotinicoid receptors, two classes of ion channels that have been well-characterized structurally and functionally. The responsiveness of these ion channels to changes in the lipid environment illustrate how ion channels, and more generally, any membrane protein, may be regulated via cellular control of membrane composition.     

The NMR solution structure of the ion channel peptaibol chrysospermin C bound to dodecylphosphocholine micelles.

Chrysospermin C is a 19-residue peptaibol capable of forming transmembrane ion channels in phospholipid bilayers. The conformation of chrysospermin C bound to dodecylphosphocholine micelles has been solved using heteronuclear NMR spectroscopy. Selective 15N-labeling and 13C-labeling of specific alpha-aminoisobutyric acid residues was used to obtain complete stereospecific assignments for all eight alpha-aminoisobutyric acid residues. Structures were calculated using 339 distance constraints and 40 angle constraints obtained from NMR data. The NMR structures superimpose with mean global rmsd values to the mean structure of 0. 27 A (backbone heavy atoms) and 0.42 A (all heavy atoms). Chrysospermin C bound to decylphosphocholine micelles displays two well-defined helices at the N-terminus (residues Phe1-Aib9) and C-terminus (Aib13-Trp-ol19). A slight bend preceding Pro14, i.e. encompassing residues 10-12, results in an angle of approximately 38 degrees between the mean axes of the two helical regions. The bend structure observed for chrysospermin C is compatible with the sequences of all 18 long peptaibols and may represent a common 'active' conformation. The structure of chrysospermin C shows clear hydrophobic and hydrophilic surfaces which would be appropriate for the formation of oligomeric ion channels.     

Emerging models of glutamate receptor ion channel structure and function

Excitatory synaptic transmission in the brain is mediated by ligand-gated ion channels (iGluRs) activated by glutamate. Distinct from other neurotransmitter receptors, the extracellular domains of iGluRs are loosely packed assemblies with two clearly distinct layers, each of which has both local and global 2-fold axes of symmetry. By contrast, the iGluR transmembrane segments have 4-fold symmetry and share a conserved pore loop architecture found in tetrameric voltage-gated ion channels. The striking layered architecture of iGluRs revealed by the 3.6?? resolution structure of an AMPA receptor homotetramer likely arose from gene fusion events that occurred early in evolution. Although this modular design has greatly facilitated biophysical and structural studies on individual iGluR domains, and suggested conserved mechanisms for iGluR gating, recent work is beginning to reveal unanticipated diversity in the structure, allosteric regulation, and assembly of iGluR subtypes.     

Determination of the structure of a membrane-incorporated ion channel. Solid-state nuclear magnetic resonance studies of gramicidin A.

Solid-state nuclear magnetic resonance (NMR) measurements on 13C-labeled analogues of the ion channel-forming peptide, gramicidin A, have been used to directly determine the structure of this peptide in lipid membranes. Seven gramicidin analogues, each labeled in a single carbonyl group of gly2, L-ala3, D-leu4, L-val7, D-leu10, D-leu12, or D-leu14 were synthesized by the solid-phase method. These gramicidin analogues were incorporated into aligned multilayers of dimyristoylphosphatidylcholine, or diether lipid bearing 14- or 16-carbon chains, at a 1:15 peptide:lipid mole ratio. Proton-enhanced, 13C, solid-state spectra were obtained at several temperatures and over a range of sample orientations with respect to the spectrometer magnetic field to permit accurate measurement of the chemical shift anisotropies. The observed anisotropies indicate that all of the labeled carbonyl bonds are oriented almost parallel to the molecular long axis and perpendicular to the lipid bilayer plane. These orientations are consistent with gramicidin forming a beta 6.3 single-strand helix that is oriented parallel to the methylene chains of the lipid molecules. Comparison of the linewidths from labeled residues that are in the innermost turn of the helix (gly2, ala3, and D-leu4), in the center of the molecule (val7), and in the turn nearest the lipid bilayer surface (D-leu10, D-leu12, and D-leu14) suggests that although the peptide behaves largely as a rigid barrel, segments of the peptide close to the membrane surface possess greater motional freedom.(ABSTRACT TRUNCATED AT 250 WORDS)     

Reptation theory of ion channel gating.

Reptation theory is a highly successful approach for describing polymer dynamics in entangled systems. In turn, this molecular process is the basis of viscoelasticity. We apply a modified version of reptation dynamics to develop an actual physical model of ion channel gating. We show that at times longer than microseconds these dynamics predict an alpha-helix-screw motion for the amphipathic protein segment that partially lines the channel pore. Such motion has been implicated in several molecular mechanics studies of both voltage-gated and transmitter-gated channels. The experimental probability density function (pdf) for this process follows t-3/2 which has been observed in several experimental systems. Reptation theory predicts that channel gating will occur on the millisecond time scale and this is consistent with experimental results from single-channel recording. We examine the consequences of reptation over random barriers and we show that, to first order, the pdf remains unchanged. In the case of a charged helix undergoing reptation in the presence of a transmembrane potential we show that the tail of the pdf will be exponential. We provide a list of practical experimental predictions to test the validity of this physical theory.     

Amyloid beta ion channel: 3D structure and relevance to amyloid channel paradigm

Alzheimer's disease (AD) is a protein misfolding disease. Early hypothesis of AD pathology posits that 39-43 AA long misfolded amyloid beta (Abeta) peptide forms a fibrillar structure and induces pathophysiological response by destabilizing cellular ionic homeostasis. Loss of cell ionic homeostasis is believed to be either indirectly due to amyloid beta-induced oxidative stress or directly by its interaction with the cell membrane and/or activating pathways for ion exchange. Significantly though, no Abeta specific cell membrane receptors are known and oxidative stress mediated pathology is only partial and indirect. Most importantly, recent studies strongly indicate that amyloid fibrils may not by themselves cause AD pathology. Subsequently, a competing hypothesis has been proposed wherein amyloid derived diffusible ligands (ADDLs) that are large Abeta oligomers (approximately >60 kDa), mediate AD pathology. No structural details, however, of these large globular units exist nor is there any known suitable mechanism by which they would induce AD pathology. Experimental data indicate that they alter cell viability by non-specifically changing the plasma membrane stability and increasing the overall ionic leakiness. The relevance of this non-specific mechanism for AD-specific pathology seems limited. Here, we provide a viable new paradigm: AD pathology mediated by amyloid ion channels made of small Abeta oligomers (trimers to octamers). This review is focused to 3D structural analysis of the Abeta channel. The presence of amyloid channels is consistent with electrophysiological and cell biology studies summarized in companion reviews in this special issue. They show ion channel-like activity and channel-mediated cell toxicity. Amyloid ion channels with defined gating and pharmacological agents would provide a tangible target for designing therapeutics for AD pathology.     

Amyloid beta ion channel: 3D structure and relevance to amyloid channel paradigm

Alzheimer's disease (AD) is a protein misfolding disease. Early hypothesis of AD pathology posits that 39-43 AA long misfolded amyloid beta (Aβ) peptide forms a fibrillar structure and induces pathophysiological response by destabilizing cellular ionic homeostasis. Loss of cell ionic homeostasis is believed to be either indirectly due to amyloid beta-induced oxidative stress or directly by its interaction with the cell membrane and/or activating pathways for ion exchange. Significantly though, no Aβ specific cell membrane receptors are known and oxidative stress mediated pathology is only partial and indirect. Most importantly, recent studies strongly indicate that amyloid fibrils may not by themselves cause AD pathology. Subsequently, a competing hypothesis has been proposed wherein amyloid derived diffusible ligands (ADDLs) that are large Aβ oligomers (∼ > 60 kDa), mediate AD pathology. No structural details, however, of these large globular units exist nor is there any known suitable mechanism by which they would induce AD pathology. Experimental data indicate that they alter cell viability by non-specifically changing the plasma membrane stability and increasing the overall ionic leakiness. The relevance of this non-specific mechanism for AD-specific pathology seems limited. Here, we provide a viable new paradigm: AD pathology mediated by amyloid ion channels made of small Aβ oligomers (trimers to octamers). This review is focused to 3D structural analysis of the Aβ channel. The presence of amyloid channels is consistent with electrophysiological and cell biology studies summarized in companion reviews in this special issue. They show ion channel-like activity and channel-mediated cell toxicity. Amyloid ion channels with defined gating and pharmacological agents would provide a tangible target for designing therapeutics for AD pathology.     

Molecular dynamics simulation of the antiamoebin ion channel: linking structure and conductance

Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K⁺ and Cl⁻ with their first solvation shells intact. The free energy barrier encountered by K⁺ is only 2.2 kcal/mol whereas Cl⁻ encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.