Molecular simulation reveals structural determinants of the hanatoxin binding in Kv2.1 channels

The carboxyl terminus of the S3 segment (S3C) in voltage-gated potassium channels was suggested to be the binding site of gating modifier toxins like hanatoxin. It has also been proposed to have a helical secondary structural arrangement. The currently available structures in high resolution for such channel molecules are restricted to regions illustrating the pore function. Therefore no further direct experimental data to elucidate the detailed mechanism for such toxin binding can be derived. In order to examine the putative three-dimensional structure of S3C and to analyze the residues required for hanatoxin binding, molecular simulation and docking were performed, based on the solution structure of hanatoxin and the structural information from mutational scanning data for the S3C fragment in Kv2.1. Our results indicate that hydrophobic and electrostatic interactions are both utilized to stabilize the toxin binding. Precise docking residues and the appropriate orientation for binding regarding amphipathic environments are also described. Compared with the functional data proposed by previous studies, the helical structural arrangement for the C-terminus of the S3 segment in voltage-gated potassium channels can therefore be further emphasized and analyzed. The possible location/orientation for toxin binding with respect to membrane distribution around the S3C segment is also discussed in this paper.     

Helical structure of the COOH terminus of S3 and its contribution to the gating modifier toxin receptor in voltage-gated ion channels

The voltage-sensing domains in voltage-gated K(+) channels each contain four transmembrane (TM) segments, termed S1 to S4. Previous scanning mutagenesis studies suggest that S1 and S2 are amphipathic membrane spanning alpha-helices that interface directly with the lipid membrane. In contrast, the secondary structure of and/or the environments surrounding S3 and S4 are more complex. For S3, although the NH(2)-terminal part displays significant helical character in both tryptophan- and alanine-scanning mutagenesis studies, the structure of the COOH-terminal portion of this TM is less clear. The COOH terminus of S3 is particularly interesting because this is where gating modifier toxins like Hanatoxin interact with different voltage-gated ion channels. To further examine the secondary structure of the COOH terminus of S3, we lysine-scanned this region in the drk1 K(+) channel and examined the mutation-induced changes in channel gating and Hanatoxin binding affinity, looking for periodicity characteristic of an alpha-helix. Both the mutation-induced perturbation in the toxin-channel interaction and in gating support the presence of an alpha-helix of at least 10 residues in length in the COOH terminus of S3. Together with previous scanning mutagenesis studies, these results suggest that, in voltage-gated K(+) channels, the entire S3 segment is helical, but that it can be divided into two parts. The NH(2)-terminal part of S3 interfaces with both lipid and protein, whereas the COOH-terminal part interfaces with water (where Hanatoxin binds) and possibly protein. A conserved proline residue is located near the boundary between the two parts of S3, arguing for the presence of a kink in this region. Several lines of evidence suggest that these structural features of S3 probably exist in all voltage-gated ion channels.     

A possible molecular mechanism of hanatoxin binding-modified gating in voltage-gated K+-channels

While S4 is known as the voltage sensor in voltage-gated potassium channels, the carboxyl terminus of S3 (S3C) is of particular interest concerning the site for gating modifier toxins like hanatoxin. The thus derived helical secondary structural arrangement for S3C, as well as its surrounding environment, has since been intensively and vigorously debated. Our previous structural analysis based on molecular simulation has provided sufficient information to describe reasonable docking conformation and further experimental designs (Lou et al., 2002. J. Mol. Recognit. 15: 175-179). However, if one only relies on such information, more advanced structure-functional interpretations for the roles S3C may play in the modification of gating behavior upon toxin binding will remain unknown. In order to have better understanding of the molecular details regarding this issue, we have performed the docking simulation with the S3C sequence from the hanatoxin-insensitive K+-channel, shaker, and analyzed the conformational changes resulting from such docking. Compared with other functional data from previous studies with respect to the proximity of the S3-S4 linker region, we suggested a significant movement of drk1 S3C, but not shaker S3C, in the direction presumably towards S4, which was comprehended as a possible factor interfering with S4 translocation during drk1 gating in the presence of toxin. In combination with the discussions for structural roles of the length of the S3-S4 linker, a possible molecular mechanism to illustrate the hanatoxin binding-modified gating is proposed.     

Localization and molecular determinants of the Hanatoxin receptors on the voltage-sensing domains of a K(+) channel

Hanatoxin inhibits voltage-gated K(+) channels by modifying the energetics of activation. We studied the molecular determinants and physical location of the Hanatoxin receptors on the drk1 voltage-gated K(+) channel. First, we made multiple substitutions at three previously identified positions in the COOH terminus of S3 to examine whether these residues interact intimately with the toxin. We also examined a region encompassing S1-S3 using alanine-scanning mutagenesis to identify additional determinants of the toxin receptors. Finally, guided by the structure of the KcsA K(+) channel, we explored whether the toxin interacts with the peripheral extracellular surface of the pore domain in the drk1 K(+) channel. Our results argue for an intimate interaction between the toxin and the COOH terminus of S3 and suggest that the Hanatoxin receptors are confined within the voltage-sensing domains of the channel, at least 20-25 A away from the central pore axis.     

Modeling of benzocaine analog interactions with the D4S6 segment of NaV4.1 voltage-gated sodium channels

Local anesthetics (LAs) are compounds that inhibit the propagation of action potentials in excitable tissues by blocking voltage-gated Na+ channels. Mutagenesis studies have demonstrated that several amino acid residues are important sites of LA interaction with the channel, but these studies provide little information regarding the molecular forces that govern drug-binding interactions, including the binding orientation of drugs. We used computational methods to construct a simple model of benzocaine analog binding with the D4S6 segment of rat skeletal muscle (NaV4.1) sodium channels. The model revealed that four hydrophobic residues form a binding cavity for neutral LAs, and docking studies indicated that increasing hydrophobicity among the benzocaine analogs allowed a better fit within the binding cavity. The similarities between our simple model and published experimental data suggested that modeling of LA interactions with sodium channels, along with experimental approaches, could further enhance our understanding of LA interactions with sodium 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.     

The helical character of the S6 segment of TRPV1 channels

The era of the molecular structure of ion channels has revealed that their transmembrane segments are alpha helices, as was suspected from hydropathy analysis and experimental data. TRP channels are recent additions to the known families of ion channels and little structural data is available. In a recent work, we explored the conformational changes occurring at the putative S6 segment of TRPV1 channels and observed a periodicity of chemical modification of residues suggestive of an alpha helical structure. Further analysis of the periodicity of the disposition of hydrophobic residues in the S6 segment, suggests that the general architecture of the TRPV1 S6 segment, is very similar to that of voltage-dependent channels of known structure—an aqueous cavity lined by an amphipathic alpha helix, with most of the hydrophobic residues pointing into it.     

Solution structure of omega-grammotoxin SIA,a gating modifier of P/Q and N-type Ca(2+) channel

omega-Grammotoxin SIA (GrTx) is a 36 amino acid residue protein toxin from spider venom that inhibits P/Q and N-type voltage-gated Ca(2+) channels by modifying voltage-dependent gating. We determined the three-dimensional structure of GrTx using NMR spectroscopy. The toxin adopts an "inhibitor cystine knot" motif composed of two beta-strands (Leu19-Cys21 and Cys30-Trp32) and a beta-bulge (Trp6, Gly7-Cys30) with a +2x, -1 topology, which are connected by four chain reversals. Although GrTx was originally identified as an inhibitor of voltage-gated Ca(2+) channel, it also binds to K(+) channels with lower affinity. A similar cross-reaction was observed for Hanatoxin1 (HaTx), which binds to the voltage-sensing domains of K(+) and Ca(2+) channels with different affinities. A detailed comparison of the GrTx and HaTx structures identifies a conserved face containing a large hydrophobic patch surrounded by positively charged residues. The slight differences in the surface shape, which result from the orientation of the surface aromatic residues and/or the distribution of the charged residues, may explain the differences in the binding affinity of these gating modifiers with different voltage-gated ion channels.     

Homology model and targeted mutagenesis identify critical residues for arachidonic acid inhibition of Kv4 channels

Polyunsaturated fatty acids such as arachidonic acid (AA) exhibit inhibitory modulation of Kv4 potassium channels. Molecular docking approaches using a Kv4.2 homology model predicted a membrane-embedded binding pocket for AA comprised of the S4-S5 linker on one subunit and several hydrophobic residues within S3, S5 and S6 from an adjacent subunit. The pocket is conserved among Kv4 channels. We tested the hypothesis that modulatory effects of AA on Kv4.2/KChIP channels require access to this site. Targeted mutation of a polar residue (K318) and a nonpolar residue (G314) within the S4-S5 linker as well as a nonpolar residue in S3 (V261) significantly impaired the effects of AA on K⁺ currents in Xenopus oocytes. These residues may be important in stabilizing (K318) or regulating access to (V261, G314) the negatively charged carboxylate moiety on the fatty acid. Structural specificity was supported by the lack of disruption of AA effects observed with mutations at residues located near, but not within the predicted binding pocket. Furthermore, we found that the crystal structure of the related Kv1.2/2.1 chimera lacks the structural features present in the proposed AA docking site of Kv4.2 and the Kv1.2/2.1 K⁺ currents were unaffected by AA. We simulated the mutagenic substitutions in our Kv4.2 model to demonstrate how specific mutations may disrupt the putative AA binding pocket. We conclude that AA inhibits Kv4 channel currents and facilitates current decay by binding within a hydrophobic pocket in the channel in which K318 within the S4-S5 linker is a critical residue for AA interaction.     

Molecular basis for Kv1.5 channel block: conservation of drug binding sites among voltage-gated K+ channels

Kv1.5 channels conduct the ultrarapid delayed rectifier current (IKur) that contributes to action potential repolarization of human atrial myocytes. Block of these channels has been proposed as a treatment for atrial arrhythmias. Here we report a novel and potent inhibitor of Kv1.5 potassium channels, N-benzyl-N-pyridin-3-yl-methyl-2-(toluene-4-sulfonylamino)-benzamide hydrochloride (S0100176), which exhibits features consistent with preferential block of the open state. The IC50 of S0100176 for Kv1.5 expressed in Xenopus oocytes was 0.7 microm. Ala-scanning mutagenesis within the pore helix and the S6 segment, regions that form the walls of the central cavity, was combined with voltage clamp analysis to identify point mutations that altered drug affinity. This approach identified Thr-479, Thr-480, Val-505, Ile-508, and Val-512 as the most important residues for block by S0100176. Mutations of these key residues to Ala or other amino acids caused marked changes in the IC50 of S0100176 (p<0.01). For example, the IC50 of S0100176 increased 362-fold for T480A, 26-fold for V505A, 150-fold for I508A, and 99-fold for V512A. We used modeling to dock S0100176 into the inner cavity of a Kv1.5 pore homology model that was generated based on the crystal structure of KcsA. The docking predicted that the five residues identified by the Ala scan were positioned less than 4.5 A from the compound. Based on the homology models, the positions of the five amino acids identified to interact with S0100176 face toward the central cavity and overlap with putative binding sites for other blockers and voltage-gated potassium channels.     

Cysteines control the N- and C-linker-dependent gating of KCNH1 potassium channels

KCNH1 (EAG1) is a member of the Kv family of voltage-gated potassium channels. However, KCNH1 channels also show some amino-acid sequence similarity to cyclic-nucleotide-regulated channels: they harbor an N-terminal PAS domain, a C-terminal cyclic nucleotide binding homology domain (cNBHD), and N- and C-terminal binding sites for calmodulin. Another notable feature is the channels' high sensitivity toward oxidative modification. Using human KCNH1 expressed in Xenopus oocytes and HEK 293 cells we investigated how oxidative modification alters channel function. Intracellular application of H₂O₂ or cysteine-specific modifiers potently inhibited KCNH1 channels in two phases. Our systematic cysteine mutagenesis study showed that the rapid and dominant phase was attributed to a right-shift in the voltage dependence of activation, caused by chemical modification of residues C145 and C214. The slow component depended on the C-terminal residues C532 and C562. The cysteine pairs are situated at structural elements linking the transmembrane S1 segment with the PAS domain (N-linker) and the transmembrane channel gate S6 with the cNBH domain (C-linker), respectively. The functional state of KCNH1 channels is determined by the oxidative status of these linkers that provide an additional dimension of channel regulation.     

Evaluation of lipid‐binding properties of the N‐terminal helical segments in human apolipoprotein A‐I using fragment peptides

Although the N‐terminal region in human apolipoprotein (apo) A‐I is thought to stabilize the lipid‐free structure of the protein, its role in lipid binding is unknown. Using synthetic fragment peptides, we examined the lipid‐binding properties of the first 43 residues (1–43) of apoA‐I in comparison with residues 44–65 and 220–241, which have strong lipid affinity in the molecule. Circular dichroism measurements demonstrated that peptides corresponding to each segment have potential propensity to form α‐helical structure in trifluoroethanol. Spectroscopic and thermodynamic measurements revealed that apoA‐I (1–43) peptide has the strong ability to bind to lipid vesicles and to form α‐helical structure comparable to apoA‐I (220–241) peptide. Substitution of Tyr‐18 located at the center of the most hydrophobic region in residues 1–43 with a helix‐breaking proline resulted in the impaired lipid binding, indicating that the α‐helical structure in this region is required to trigger the lipid binding. In contrast, apoA‐I (44–65) peptide exhibited a lower propensity to form α‐helical structure upon binding to lipid, and apoA‐I (44–65/S55P) peptide exhibited diminished, but not completely impaired, lipid binding, suggesting that the central region of residues 44–65 is not pivotally involved in the formation of the α‐helical structure and lipid binding. These results indicate that the most N‐terminal region of apoA‐I molecule, residues 1–43, contributes to the lipid interaction of apoA‐I through the hydrophobic helical residues. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular mapping of general anesthetic sites in a voltage-gated ion channel

Several voltage-gated ion channels are modulated by clinically relevant doses of general anesthetics. However, the structural basis of this modulation is not well understood. Previous work suggested that n-alcohols and inhaled anesthetics stabilize the closed state of the Shaw2 voltage-gated (Kv) channel (K-Shaw2) by directly interacting with a discrete channel site. We hypothesize that the inhibition of K-Shaw2 channels by general anesthetics is governed by interactions between binding and effector sites involving components of the channel's activation gate. To investigate this hypothesis, we applied Ala/Val scanning mutagenesis to the S4-S5 linker and the post-PVP S6 segment, and conducted electrophysiological analysis to evaluate the energetic impact of the mutations on the inhibition of the K-Shaw2 channel by 1-butanol and halothane. These analyses identified residues that determine an apparent binding cooperativity and residue pairs that act in concert to modulate gating upon anesthetic binding. In some instances, due to their critical location, key residues also influence channel gating. Complementing these results, molecular dynamics simulations and in silico docking experiments helped us visualize possible anesthetic sites and interactions. We conclude that the inhibition of K-Shaw2 by general anesthetics results from allosteric interactions between distinct but contiguous binding and effector sites involving inter- and intrasubunit interfaces.     

The investigation of interactions of kappa-Hefutoxin1 with the voltage-gated potassium channels: a computational simulation

Kappa-Hefutoxin1 is a K(+) channel-blocking toxin from the scorpion Heterometrus fluvipes. It is a 22-residue protein that adapts a novel fold of two parallel helices linked by two disulfide bridges without beta-sheets. Recognition of interactions of kappa-Hefutoxin1 with the voltage-gated potassium channels, Kv1.1, Kv1.2, and Kv1.3, was studied by 3D-Dock software package. All structures of kappa-Hefutoxin1 were considered during the simulations, which indicated that even small changes in the structure of kappa-Hefutoxin1 considerably affected both the recognition and the binding between kappa-Hefutoxin1 and the Kv1 channels. kappa-Hefutoxin1 is located around the extracellular part of the Kv1 channels, making contacts with its helices. Lys 19, Tyr 5, Arg 6, Trp 9, or Arg 10 in the toxin and residues Asp 402, His 404, Thr 407,Gly 401, and Asp 386 in each subunit of the Kv potassium channel are the key residues for the toxin-channel recognition. Moreover, the simulation result demonstrates that the hydrophobic interactions are important in interaction of negatively charged toxins with potassium channels. The results of our docking/molecular dynamics simulations indicate that our 3D model structure of the kappa-Hefutoxin1-complex is both reasonable and acceptable and could be helpful for smarter drug design and the blocking agents of Kv1 channels.     

Interaction between extracellular Hanatoxin and the resting conformation of the voltage-sensor paddle in Kv channels

In voltage-activated potassium (Kv) channels, basic residues in S4 enable the voltage-sensing domain to move in response to membrane depolarization and thereby trigger the activation gate to open. In the X-ray structure of the KvAP channel, the S4 helix is located near the intracellular boundary of the membrane where it forms a "voltage-sensor paddle" motif with the S3b helix. It has been proposed that the paddle is lipid-exposed and that it translocates through the membrane as it activates. We studied the interaction of externally applied Hanatoxin with the voltage-sensor paddle in Kv channels and show that the toxin binds tightly even at negative voltages where the paddle is resting and the channel is closed. Moreover, measurements of gating charge movement suggest that Hanatoxin interacts with and stabilizes the resting paddle. These findings point to an extracellular location for the resting conformation of the voltage-sensor paddle and constrain its transmembrane movements during activation.     

Structural and functional consequences of the presence of a fourth disulfide bridge in the scorpion short toxins: solution structure of the potassium channel inhibitor HsTX1

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide-bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide-bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.     

Structure–function studies on the voltage-gated sodium channel

Recent research on structure–function relationships aspects of voltage-gated sodium channels (VGSCs) are reviewed. Data issued from the literature are summarized and compared, including results from our own studies. The latter deal with the effects of drug binding, deglycosylation and the role of hydrophobic residues in the voltage sensors. Methods mainly consist of circular dichroism (CD) to asses the channel's secondary structure and conductance measurements after reconstitution into planar lipid bilayers. Molecular modelling was also used to tentatively explain experimental data. Since 30% of the channel's mass are glycoconjugates, the effects of removing them were first investigated. Then, the effects of the neurotoxin Batrachotoxin and the anticonvulsant Lamotrigine were studied. Both drugs induced a significant increase in the channel's helical content and a molecular model shows that lamotrigine interacts with residues previously identified as forming the binding sites in the pore. Finally, the role of hydrophobic residues with long sidechains in the voltage sensors (S4s) was investigated. Recent research on related studies on VGSCs are discussed.     

Cooperative, non-specific binding of a zinc finger peptide to DNA.

The DNA binding and structural properties of Xfin-31 (Lee, M.S., Gippert, G.P., Soman, K.V., Case, D.A. and Wright, P.E., 1989, Science 245, 635-637), a twenty five amino acid zinc finger peptide, in the reduced, oxidized and zinc complex forms, as well as the fourteen residue helical segment of the zinc finger (residues 12-25) have been compared using affinity coelectrophoresis (ACE) and circular dichroism (CD) spectroscopy. The zinc complex and oxidized peptides bind cooperatively to DNA although the cooperativity factor, omega, is more than 15-fold greater for the zinc complex. The reduced peptide in the absence of zinc and the helical segment do not bind cooperatively (omega = 1). Hence, the binding constant for singly contiguous sites (K omega) ranges over 100-fold for the various peptides even though the intrinsic binding constants (K) are similar. An increase in binding order and affinity for the other forms of Xfin-31 is correlated with an increasing similarity of the CD spectrum to that of the Xfin-31 zinc complex. The surprising DNA binding activity of the oxidized peptide may result from hydrophobic interactions between the amino-terminal loop formed by the Cys3-Cys6 disulfide bond and conserved hydrophobic residues in the carboxyl-terminal segment. Xfin-31 may be a particularly useful model for studying several poorly understood aspects of cooperative, non-specific DNA binding since it is small, has a stable, well-defined structure, and structures of zinc fingers bound to DNA have been determined.     

Potential role of the membrane in hERG channel functioning and drug-induced long QT syndrome

The human ether-à-go-go related gene (hERG) potassium channels are located in the myocardium cell membrane where they ensure normal cardiac activity. The binding of drugs to this channel, a side effect known as drug-induced (acquired) long QT syndrome (ALQTS), can lead to arrhythmia or sudden cardiac death. The hERG channel is a unique member of the family of voltage-gated K⁺ channels because of the long extracellular loop connecting its transmembrane S5 helix to the pore helix in the pore domain. Considering the proximal position of the S5-P linker to the membrane surface, we have investigated the interaction of its central segment I⁵⁸³-Y⁵⁹⁷ with bicelles. Liquid and solid-state NMR experiments as well as circular dichroism results show a strong affinity of the I⁵⁸³-Y⁵⁹⁷ segment for the membrane where it would sit on the surface with no defined secondary structure. A structural dependence of this segment on model membrane composition was observed. A helical conformation is favoured in detergent micelles and in the presence of negative charges. Our results suggest that the interaction of the S5-P linker with the membrane could participate in the stabilization of transient channel conformations, but helix formation would be triggered by interactions with other hERG domains. Because potential drug binding sites on the S5-P linker have been identified, we have explored the role of this segment in ALQTS. Four LQTS-liable drugs were studied which showed more affinity for the membrane than this hERG segment. Our results, therefore, identify two possible roles for the membrane in channel functioning and ALQTS.     

Helical structure and packing orientation of the S2 segment in the Shaker K+ channel

Six transmembrane segments, S1-S6, cluster around the central pore-forming region in voltage-gated K+ channels. To investigate the structural characteristics of the S2 segment in the Shaker K+ channel, we replaced each residue in S2 singly with tryptophan (or with alanine for the native tryptophan). All but one of the 23 Trp mutants expressed voltage-dependent K+ currents in Xenopus oocytes. The effects of the mutations were classified as being of low or high impact on channel gating properties. The periodicity evident in the effects of these mutations supports an alpha-helical structure for the S2 segment. The high- and low-impact residues cluster onto opposite faces of a helical wheel projection of the S2 segment. The low-impact face is also tolerant of single mutations to asparagine. All results are consistent with the idea that the low-impact face projects toward membrane lipids and that changes in S2 packing occur upon channel opening. We conclude that the S2 segment is a transmembrane alpha helix and that the high-impact face packs against other transmembrane segments in the functional channel.