Molecular Docking of the Scorpion Toxin Tc1 to the Structural Model of the Voltage-gated Potassium Channel Kv1.1 from Human Homo sapiens

Abstract

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K⁺ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

A set of homology models of pore loop domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6

Homology models of the pore loop domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were generated based on the crystallographic structure of KcsA. The results of amino acid sequence alignment indicate that these Kv channels are composed of two structurally and functionally independent domains: the N-terminal 'voltage sensor' domain and the C-terminal 'pore loop' domain. The homology models reveal that the pore loop domains of these Kv channels exhibit similar folds to those of KcsA. The structural features and specific packing of aromatic residues around the selectivity filter of these Kv channels are nearly identical to those of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. The distribution of polar and nonpolar side chains on the surfaces of the KcsA and Kv channels reveals that they exhibit a segregation of side chains common to most integral membrane proteins. As the hydrogen bond between Glu71 and Asp80 in KcsA plays an important role in stabilizing the channel, the substituted Val residue in the Kv family corresponding to Glu71 of KcsA stabilizes the channel by making hydrophobic contact with Tyr residue from the signature sequence of the selectivity filter. The homology models of these Kv channels provide particularly attractive subjects for further structure-based studies.     

KcsA crystal structure as framework for a molecular model of the Na(+) channel pore

The crystal structure of the pore-forming part of the KcsA bacterial K(+)-selective channel suggests a possible motif for related voltage-gated channels. We examined the hypothesis that the spacial orientation of the KcsA M1 and M2 alpha-helices also predicts the backbone location of S5 and S6 helices of the voltage-gated Na(+) channel. That channel's P region structure is expected to be different because selectivity is determined by side-chain interactions rather than by main-chain carbonyls, and its outer vestibule accommodates relatively large toxin molecules, tetrodotoxin (TTX) and saxitoxin (STX), which interact with selectivity ring residues. The Na(+) channel P loop was well-modeled by the alpha-helix-turn-beta-strand motif, which preserves the relationships for toxin interaction with the Na(+) channel found experimentally. This outer vestibule was docked into the extracellular part of the inverted teepee structure formed by the S5 and S6 helices that were spacially located by coordinates of the KcsA M1 and M2 helix main chains [Doyle et al. (1998) Science 280, 69-74], but populated with side chains of the respective S5 and S6 structures. van der Waals contacts were optimized with minimal adjustment of the S5, S6, and P loop structures, forming a densely packed pore structure. Nonregular external S5-P and P-S6 segments were not modeled here, except the P-S6 segment of domain II. The resulting selectivity region structure is consistent with Na(+) channel permeation properties, offering suggestions for the molecular processes involved in selectivity. The ability to construct a Na(+) channel pore model consistent with most of the available biophysical and mutational information suggests that the KcsA structural framework may be conserved in voltage-gated channels.     

Structural conservation of the pores of calcium-activated and voltage-gated potassium channels determined by a sea anemone toxin.

The structurally defined sea anemone peptide toxins ShK and BgK potently block the intermediate conductance, Ca(2+)-activated potassium channel IKCa1, a well recognized therapeutic target present in erythrocytes, human T-lymphocytes, and the colon. The well characterized voltage-gated Kv1.3 channel in human T-lymphocytes is also blocked by both peptides, although ShK has a approximately 1,000-fold greater affinity for Kv1.3 than IKCa1. To gain insight into the architecture of the toxin receptor in IKCa1, we used alanine-scanning in combination with mutant cycle analyses to map the ShK-IKCa1 interface, and compared it with the ShK-Kv1.3 interaction surface. ShK uses the same five core residues, all clustered around the critical Lys(22), to interact with IKCa1 and Kv1.3, although it relies on a larger number of contacts to stabilize its weaker interactions with IKCa1 than with Kv1.3. The toxin binds to IKCa1 in a region corresponding to the external vestibule of Kv1.3, and the turret and outer pore of the structurally defined bacterial potassium channel, KcsA. Based on the NMR structure of ShK, we deduce the toxin receptor in IKCa1 to have x-y dimensions of approximately 22 A, a diameter of approximately 31 A, and a depth of approximately 8 A; we estimate that the ion selectivity lies approximately 13 A below the outer lip of the toxin receptor. These dimensions are in good agreement with those of the KcsA channel determined from its crystal structure, and the inferred structure of Kv1.3 based on mapping with scorpion toxins. Thus, these distantly related channels exhibit architectural similarities in the outer pore region. This information could facilitate development of specific and potent modulators of the therapeutically important IKCa1 channel.     

Influence of pore residues on permeation properties in the Kv2.1 potassium channel. Evidence for a selective functional interaction of K+ with the outer vestibule

The Kv2.1 potassium channel contains a lysine in the outer vestibule (position 356) that markedly reduces open channel sensitivity to changes in external [K(+)]. To investigate the mechanism underlying this effect, we examined the influence of this outer vestibule lysine on three measures of K(+) and Na(+) permeation. Permeability ratio measurements, measurements of the lowest [K(+)] required for interaction with the selectivity filter, and measurements of macroscopic K(+) and Na(+) conductance, were all consistent with the same conclusion: that the outer vestibule lysine in Kv2.1 interferes with the ability of K(+) to enter or exit the extracellular side of the selectivity filter. In contrast to its influence on K(+) permeation properties, Lys 356 appeared to be without effect on Na(+) permeation. This suggests that Lys 356 limited K(+) flux by interfering with a selective K(+) binding site. Combined with permeation studies, results from additional mutagenesis near the external entrance to the selectivity filter indicated that this site was located external to, and independent from, the selectivity filter. Protonation of a naturally occurring histidine in the same outer vestibule location in the Kv1.5 potassium channel produced similar effects on K(+) permeation properties. Together, these results indicate that a selective, functional K(+) binding site (e.g., local energy minimum) exists in the outer vestibule of voltage-gated K(+) channels. We suggest that this site is the location of K(+) hydration/dehydration postulated to exist based on the structural studies of KcsA. Finally, neutralization of position 356 enhanced outward K(+) current magnitude, but did not influence the ability of internal K(+) to enter the pore. These data indicate that in Kv2.1, exit of K(+) from the selectivity filter, rather than entry of internal K(+) into the channel, limits outward current magnitude. We discuss the implications of these findings in relation to the structural basis of channel conductance in different K(+) channels.     

Modeling of the outer vestibule and selectivity filter of the L-type Ca2+ channel

Using the KcsA bacterial K+ channel crystal structure [Doyle, D. A., et al. (1998) Science 280, 69-74] and the model of the outer vestibule of the Na+ channel [Lipkind, G. M., and Fozzard, H. A. (2000) Biochemistry 39, 8161-8170] as structural templates, we propose a structural model of the outer vestibule and selectivity filter of the pore of the Ca2+ channel (alpha1C or Ca(v)1.2). The Ca2+ channel P loops were modeled by alpha-helix-turn-beta-strand motifs, with the glutamate residues of the EEEE motif located in the turns. P loops were docked in the extracellular part of the inverted teepee structure formed by S5 and S6 alpha-helices with backbone coordinates from the M1 and M2 helices of the KcsA crystal structure. This construction results in a conical outer vestibule that tapers to the selectivity filter at the bottom. The modeled selectivity ring forms a wide open pore ( approximately 6 A) in the absence of Ca2+. When Ca2+ is present ( approximately 1 microM), all four glutamate side chains move to the center and form a cage around the dehydrated Ca2+ ion, blocking the pore. In the millimolar concentration range, Ca2+ also interacts with two low-affinity sites located externally and internally, which were modeled by the same carboxylate groups of the selectivity filter. Calculation of the resulting electrostatic potentials show that the single Ca2+ ion is located in an electrostatic trap. Only when three Ca2+ ions are bound simultaneously in the high- and low-affinity sites of the selectivity filter is Ca2+ able to overcome electrostatic attraction, permitting Ca2+ flux.     

External tetraethylammonium as a molecular caliper for sensing the shape of the outer vestibule of potassium channels

External tetraethylammonium (TEA+) blocked currents through Kv1.1 channels in a voltage-independent manner between 0 and 100 mV. Lowering extracellular pH (pHo) increased the Kd for TEA+ block. A histidine at position 355 in the Kv1.1 channel protein (homologous to Shaker 425) was responsible for this pH-dependent reduction of TEA+ sensitivity, since the TEA+ effect became independent of pHo after chemical modification of the Kv1.1 channel at H355 and in the H355G and H355K mutant Kv1.1 channels. The Kd values for TEA+ block of the two mutant channels (0.34 +/- 0.06 mM, n = 7 and 0.84 +/- 0. 09 mM, n = 13, respectively) were as expected for a vestibule containing either no or a total of four positive charges at position 355. In addition, the pH-dependent TEA+ effect in the wt Kv1.1 channel was sensitive to the ionic strength of the solution. All our observations are consistent with the idea that lowering pHo increased protonation of H355. This increase in positive charge at H355 will repel TEA+ electrostatically, resulting in a reduction of the effective [TEA+]o at the receptor site. From this reduction we can estimate the distance between TEA+ and each of the four histidines at position 355 to be approximately 10 A, assuming fourfold symmetry of the channel and assuming that TEA+ binds in the central axis of the pore. This determination of the dimensions of the outer vestibule of Kv1.1 channels confirms and extends earlier reports on K+ channels using crystal structure data as well as peptide toxin/channel interactions and points out a striking similarity between vestibules of Kv1.1 and KcsA channels.     

Structure of the BgK-Kv1.1 complex based on distance restraints identified by double mutant cycles. Molecular basis for convergent evolution of Kv1 channel blockers

A structural model of BgK, a sea anemone toxin, complexed with the S5-S6 region of Kv1.1, a voltage-gated potassium channel, was determined by flexible docking under distance restraints identified by a double mutant cycles approach. This structure provides the molecular basis for identifying the major determinants of the BgK-Kv1.1 channel interactions involving the BgK dyad residues Lys(25) and Tyr(26). These interactions are (i) electrostatic interactions between the extremity of Lys(25) side chain and carbonyl oxygen atoms of residues from the channel selectivity filter that may be strengthened by solvent exclusion provided by (ii) hydrophobic interactions involving BgK residues Tyr(26) and Phe(6) and Kv1.1 residue Tyr(379) whose side chain protrudes in the channel vestibule. In other Kv1 channel-BgK complexes, these interactions are likely to be conserved, implicating both conserved and variable residues from the channels. The data suggest that the conservation in sea anemone and scorpion potassium channel blockers of a functional dyad composed of a lysine, and a hydrophobic residue reflects their use of convergent binding solutions based on a crucial interplay between these important conserved interactions.     

Engineering a potent and specific blocker of voltage-gated potassium channel Kv1.3, a target for autoimmune diseases

A polypeptide toxin extracted from scorpion venom, OSK1, is modified such that its potency is drastically enhanced in blocking one class of voltage-gated potassium channels, Kv1.3, which is a pharmacological target for immunosuppressive therapy. The bound complex of Kv1.3 and OSK1 reveals that one lysine residue of the toxin is in the proximity of another lysine residue on the external vestibule of the channel, just outside of the selectivity filter. This unfavorable electrostatic interaction is eliminated by interchanging the positions of two amino acids in the toxin. The potentials of mean force of the wild-type and mutant OSK1 bound to Kv1.1-Kv1.3 channels are constructed using molecular dynamics, and the half-maximal inhibitory concentration (IC(50)) of each toxin-channel complex is computed. We show that the IC(50) values predicted for three toxins and three channels match closely with experiment. Kv1.3 is half-blocked by 0.2 pM mutant OSK1; it is >10000-fold more specific for this channel than for Kv1.1 and Kv1.2.     

Structural differences of bacterial and mammalian K+ channels

Using a peptide toxin, kaliotoxin (KTX), we gained new insight into the topology of the pore region of a voltage-gated potassium channel, mKv1.1. In order to find new interactions between mKv1.1 and KTX, we investigated the pH dependence of KTX block which was stronger at pH(o) 6.2 compared with pH(o) 7.4. Using site-directed mutagenesis on the channel and the toxin, we found that protonation of His(34) in KTX caused the pH(o) dependence of KTX block. Glu(350) and Glu(353) in mKv1.1, which interact with His(34) in KTX, were calculated to be 4 and 7 A away from His(34)/KTX, respectively. Docking of KTX into a homology model of mKv1.1 based on the KcsA crystal structure using this and other known interactions as constraints showed structural differences between mKv1.1 and KcsA within the turret (amino acids 348-357). To satisfy our data, we would have to modify the KcsA crystal structure for the mKv1.1 channel orienting Glu(350) 7 A and Glu(353) 4 A more toward the center of the pore compared with KcsA. This would place Glu(350) 15 A and Glu(353) 11 A away from the center of the pore instead of the distances for the equivalent KcsA residues with 22 A for Gly(53) and 15 A for Gly(56), respectively. Bacterial and mammalian potassium channels may have structural differences regarding the turret of the outer pore vestibule. This topological difference between both channel types may have substantial influence on structure-guided development of new drugs for mammalian potassium channels by rational drug design.     

Looking over toxin-K(+) channel interactions. Clues from the structural and functional characterization of α-KTx toxin Tc32, a Kv1.3 channel blocker

α-KTx toxin Tc32, from the Amazonian scorpion Tityus cambridgei, lacks the dyad motif, including Lys27, characteristic of the family and generally associated with channel blockage. The toxin has been cloned and expressed for the first time. Electrophysiological experiments, by showing that the recombinant form blocks Kv1.3 channels of olfactory bulb periglomerular cells like the natural Tc32 toxin, when tested on the Kv1.3 channel of human T lymphocytes, confirmed it is in an active fold. The nuclear magnetic resonance-derived structure revealed it exhibits an α/β scaffold typical of the members of the α-KTx family. TdK2 and TdK3, all belonging to the same α-KTx 18 subfamily, share significant sequence identity with Tc32 but diverse selectivity and affinity for Kv1.3 and Kv1.1 channels. To gain insight into the structural features that may justify those differences, we used the recombinant Tc32 nuclear magnetic resonance-derived structure to model the other two toxins, for which no experimental structure is available. Their interaction with Kv1.3 and Kv1.1 has been investigated by means of docking simulations. The results suggest that differences in the electrostatic features of the toxins and channels, in their contact surfaces, and in their total dipole moment orientations govern the affinity and selectivity of toxins. In addition, we found that, regardless of whether the dyad motif is present, it is always a Lys side chain that physically blocks the channels, irrespective of its position in the toxin sequence.     

Sequence-function analysis of the K+-selective family of ion channels using a comprehensive alignment and the KcsA channel structure

Sequence-function analysis of K(+)-selective channels was carried out in the context of the 3.2 A crystal structure of a K(+) channel (KcsA) from Streptomyces lividans (Doyle et al., 1998). The first step was the construction of an alignment of a comprehensive set of K(+)-selective channel sequences forming the putative permeation path. This pathway consists of two transmembrane segments plus an extracellular linker. Included in the alignment are channels from the eight major classes of K(+)-selective channels from a wide variety of species, displaying varied rectification, gating, and activation properties. Segments of the alignment were assigned to structural motifs based on the KcsA structure. The alignment's accuracy was verified by two observations on these motifs: 1), the most variability is shown in the turret region, which functionally is strongly implicated in susceptibility to toxin binding; and 2), the selectivity filter and pore helix are the most highly conserved regions. This alignment combined with the KcsA structure was used to assess whether clusters of contiguous residues linked by hydrophobic or electrostatic interactions in KcsA are conserved in the K(+)-selective channel family. Analysis of sequence conservation patterns in the alignment suggests that a cluster of conserved residues is critical for determining the degree of K(+) selectivity. The alignment also supports the near-universality of the "glycine hinge" mechanism at the center of the inner helix for opening K channels. This mechanism has been suggested by the recent crystallization of a K channel in the open state. Further, the alignment reveals a second highly conserved glycine near the extracellular end of the inner helix, which may be important in minimizing deformation of the extracellular vestibule as the channel opens. These and other sequence-function relationships found in this analysis suggest that much of the permeation path architecture in KcsA is present in most K(+)-selective channels. Because of this finding, the alignment provides a robust starting point for homology modeling of the permeation paths of other K(+)-selective channel classes and elucidation of sequence-function relationships therein. To assay these applications, a homology model of the Shaker A channel permeation path was constructed using the alignment and KcsA as the template, and its structure evaluated in light of established structural criteria.     

Pharmacology and Surface Electrostatics of the K Channel Outer Pore Vestibule

In spite of a generally well-conserved outer vestibule and pore structure, there is considerable diversity in the pharmacology of K channels. We have investigated the role of specific outer vestibule charged residues in the pharmacology of K channels using tetraethylammonium (TEA) and a trivalent TEA analog, gallamine. Similar to Shaker K channels, gallamine block of Kv3.1 channels was more sensitive to solution ionic strength than was TEA block, a result consistent with a contribution from an electrostatic potential near the blocking site. In contrast, TEA block of another type of K channel (Kv2.1) was insensitive to solution ionic strength and these channels were resistant to block by gallamine. Neutralizing either of two lysine residues in the outer vestibule of these Kv2.1 channels conferred ionic strength sensitivity to TEA block. Kv2.1 channels with both lysines neutralized were sensitive to block by gallamine, and the ionic strength dependence of this block was greater than that for TEA. These results demonstrate that Kv3.1 (like Shaker) channels contain negatively charged residues in the outer vestibule of the pore that influence quaternary ammonium pharmacology. The presence of specific lysine residues in wild-type Kv2.1 channels produces an outer vestibule with little or no net charge, with important consequences for quaternary ammonium block. Neutralizing these key lysines results in a negatively charged vestibule with pharmacological properties approaching those of other types of K channels.     

Mutating a critical lysine in ShK toxin alters its binding configuration in the pore-vestibule region of the voltage-gated potassium channel,Kv1.3

The voltage-gated potassium channel in T lymphocytes, Kv1.3, an important target for immunosuppressants, is blocked by picomolar concentrations of the polypeptide ShK toxin and its analogue ShK-Dap22. ShK-Dap22 shows increased selectivity for Kv1.3, and our goal was to determine the molecular basis for this selectivity by probing the interactions of ShK and ShK-Dap22 with the pore and vestibule of Kv1.3. The free energies of interactions between toxin and channel residues were measured using mutant cycle analyses. These data, interpreted as approximate distance restraints, guided molecular dynamics simulations in which the toxins were docked with a model of Kv1.3 based on the crystal structure of the bacterial K(+)-channel KcsA. Despite the similar tertiary structures of the two ligands, the mutant cycle data imply that they make different contacts with Kv1.3, and they can be docked with the channel in configurations that are consistent with the mutant cycle data for each toxin but quite distinct from one another. ShK binds to Kv1.3 with Lys22 occupying the negatively charged pore of the channel, whereas the equivalent residue in ShK-Dap22 interacts with residues further out in the vestibule, producing a significant change in toxin orientation. The increased selectivity of ShK-Dap22 is achieved by strong interactions of Dap22 with His404 and Asp386 on Kv1.3, with only weak interactions between the channel pore and the toxin. Potent and specific blockade of Kv1.3 apparently occurs without insertion of a positively charged residue into the channel pore. Moreover, the finding that a single residue substitution alters the binding configuration emphasizes the need to obtain consistent data from multiple mutant cycle experiments in attempts to define protein interaction surfaces using these data.     

Structural similarities between glutamate receptor channels and K(+) channels examined by scanning mutagenesis

The pores of glutamate receptors and K(+) channels share sequence homology, suggesting a conserved secondary structure. Scanning mutagenesis with substitution of alanine and tryptophan in GluR6 channels was performed based on the structure of KcsA. Our assay used disruption of voltage-dependent polyamine block to test for changes in the packing of pore-forming regions. Alanine scanning from D567 to R603 revealed reduced rectification resulting from channel block in two regions. A periodic pattern from F575 to M589 aligned with the pore helix in KcsA, whereas a cluster of sensitive positions around Q590, a site regulated by RNA editing, mapped to the selectivity filter in KcsA. Tryptophan scanning from D567 to R603 revealed similar patterns, but with a complete disruption of spermine block for 7 out of the 37 positions and a pM dissociation constant for Q590W. Molecular modeling with KcsA coordinates showed that GluR6 pore helix mutants disrupting polyamine block pack against M1 and M2, and are not exposed in the ion channel pore. In the selectivity filter, tryptophan creates an aromatic cage consistent with the pM dissociation constant for Q590W. A scan with glutamate substitution was used to map the cytoplasmic entrance to the pore based on charge neutralization experiments, which established that E594 was uniquely required for high affinity polyamine block. In E594Q mutants, introduction of glutamate at positions S593-L600 restored polyamine block at positions corresponding to surface-exposed residues in KcsA. Our results reinforce proposals that the pore region of glutamate receptors contains a helix and pore loop analogous to that found in K(+) channels. At the cytoplasmic entrance of the channel, a negatively charged amino acid, located in an extended loop with solvent-exposed side chains, is required for high affinity polyamine block and probably attracts cations via a through space electrostatic mechanism.     

Control of outer vestibule dynamics and current magnitude in the Kv2.1 potassium channel

In Kv2.1 potassium channels, changes in external [K+] modulate current magnitude as a result of a K+-dependent interconversion between two outer vestibule conformations. Previous evidence indicated that outer vestibule conformation (and thus current magnitude) is regulated by the occupancy of a selectivity filter binding site by K+. In this paper, we used the change in current magnitude as an assay to study how the interconversion between outer vestibule conformations is controlled. With 100 mM internal K+, rapid elevation of external [K+] from 0 to 10 mM while channels were activated produced no change in current magnitude (outer vestibule conformation did not change). When channels were subsequently closed and reopened in the presence of elevated [K+], current magnitude was increased (outer vestibule conformation had changed). When channels were activated in the presence of low internal [K+], or when K+ flow into conducting channels was transiently interrupted by an internal channel blocker, increasing external [K+] during activation did increase current magnitude (channel conformation did change). These data indicate that, when channels are in the activated state under physiological conditions, the outer vestibule conformation remains fixed despite changes in external [K+]. In contrast, when channel occupancy is lowered, (by channel closing, an internal blocker or low internal [K+]), the outer vestibule can interconvert between the two conformations. We discuss evidence that the ability of the outer vestibule conformation to change is regulated by the occupancy of a nonselectivity filter site by K+. Independent of the outer vestibule-based potentiation mechanism, Kv2.1 was remarkably insensitive to K+-dependent processes that influence current magnitude (current magnitude changed by <7% at membrane potentials between -20 and 30 mV). Replacement of two outer vestibule lysines in Kv2.1 by smaller neutral amino acids made current magnitude dramatically more sensitive to the reduction in K+ driving force (current magnitude changed by as much as 40%). When combined, these outer vestibule properties (fixed conformation during activation and the presence of lysines) all but prevent variation in Kv2.1 current magnitude when [K+] changes during activation. Moreover, the insensitivity of Kv2.1 current magnitude to changes in K+ driving force promotes a more uniform modulation of current over a wide range of membrane potentials by the K+-dependent regulation of outer vestibule conformation.     

A model of sodium channels

We have explored the permeation and blockage of ions in sodium channels, relating the channel structure to function using electrostatic profiles and Brownian dynamics simulations. The model used resembles the KcsA potassium channel with an added external vestibule and a shorter selectivity filter. The electrostatic energy landscape seen by permeating ions is determined by solving Poisson's equation. The two charged amino acid rings of Glu-Glu-Asp-Asp (EEDD) and Asp-Glu-Lys-Ala (DEKA) around the selectivity filter region are seen to play a crucial role in making the channel sodium selective, and strongly binding calcium ions such that they block the channel. Our model closely reproduces a range of experimental data including the current-voltage curves, current-concentration curves and blockage of monovalent ions by divalent ions.     

An electrostatic interaction between TEA and an introduced pore aromatic drives spring-in-the-door inactivation in Shaker potassium channels

Slow inactivation of Kv1 channels involves conformational changes near the selectivity filter. We examine such changes in Shaker channels lacking fast inactivation by considering the consequences of mutating two residues, T449 just external to the selectivity filter and V438 in the pore helix near the bottom of the selectivity filter. Single mutant T449F channels with the native V438 inactivate very slowly, and the canonical foot-in-the-door effect of extracellular tetraethylammonium (TEA) is not only absent, but the time course of slow inactivation is accelerated by TEA. The V438A mutation dramatically speeds inactivation in T449F channels, and TEA slows inactivation exactly as predicted by the foot-in-the-door model. We propose that TEA has this effect on V438A/T449F channels because the V438A mutation produces allosteric consequences within the selectivity filter and may reorient the aromatic ring at position 449. We investigated the possibility that the blocker promotes the collapse of the outer vestibule (spring-in-the-door) in single mutant T449F channels by an electrostatic attraction between a cationic TEA and the quadrupole moments of the four aromatic rings. To test this idea, we used in vivo nonsense suppression to serially fluorinate the introduced aromatic ring at the 449 position, a manipulation that withdraws electrons from the aromatic face with little effect on the shape, net charge, or hydrophobicity of the aromatic ring. Progressive fluorination causes monotonically enhanced rates of inactivation. In further agreement with our working hypothesis, increasing fluorination of the aromatic gradually transforms the TEA effect from spring-in-the-door to foot-in-the-door. We further substantiate our electrostatic hypothesis by quantum mechanical calculations.     

A model of sodium channels

We have explored the permeation and blockage of ions in sodium channels, relating the channel structure to function using electrostatic profiles and Brownian dynamics simulations. The model used resembles the KcsA potassium channel with an added external vestibule and a shorter selectivity filter. The electrostatic energy landscape seen by permeating ions is determined by solving Poisson's equation. The two charged amino acid rings of Glu-Glu-Asp-Asp (EEDD) and Asp-Glu-Lys-Ala (DEKA) around the selectivity filter region are seen to play a crucial role in making the channel sodium selective, and strongly binding calcium ions such that they block the channel. Our model closely reproduces a range of experimental data including the current-voltage curves, current-concentration curves and blockage of monovalent ions by divalent ions.     

Tight steric closure at the intracellular activation gate of a voltage-gated K(+) channel.

In voltage-gated K(+) channels (Kv), an intracellular gate regulates access from the cytoplasm to the pore by organic channel blockers and by chemical modifiers. But is ion flow itself controlled instead by constriction of the narrow selectivity filter near the extracellular surface? We find that the intracellular gate of Kv channels is capable of regulating access even by the small cations Cd(2+) and Ag(+). It can also exclude small neutral or negatively charged molecules, indicating that the gate operates by steric exclusion rather than electrostatically. Just intracellular to the gated region, channel closure does not restrict access even to very large reagents. Either these Kv channels have a broader inner entrance than seen in the KcsA crystal, even in the closed state, or the region is highly flexible (but nevertheless remains very securely closed nearby).