Search

Voltage-sensitive neuronal Kv1 channels composed of four α subunits and four associated auxiliary β subunits control neuronal excitability and neurotransmission. Limited information exists on the combinations of α subunit isoforms (i.e. Kv1.1-1.6) or their positions in the oligomers, and how these affect sensitivity to blockers. It is known that TEA (tetraethylammonium) inhibits Kv1.1 channels largely due to binding a critical tyrosine (Tyr379) in the pore, whereas Val381 at the equivalent location in Kv1.2 makes it insensitive. With the eventual aim of developing blockers for therapeutic purposes, Kv1.1 and 1.2 α subunit genes were concatenated to form combinations representing those in central neurons, followed by surface expression in HEK (human embryonic kidney)-293 cells as single-chain functional proteins. Patch-clamp recordings demonstrated the influences of the ratios and positioning of these α subunits on the biophysical and pharmacological properties of oligomeric K+ channels. Raising the ratio of Kv1.1 to Kv1.2 in Kv1.2-1.2-1.1-1.2 led to the resultant channels being more sensitive to TEA and also affected their biophysical parameters. Moreover, mutagenesis of one or more residues in the first Kv1.2 to resemble those in Kv1.1 increased TEA sensitivity only when it is adjacent to a Kv1.1 subunit, whereas placing a non-interactive subunit between these two diminished susceptibility. The findings of the present study support the possibility of α subunits being precisely arranged in Kv1 channels, rather than being randomly assembled. This is important in designing drugs with abilities to inhibit particular oligomeric Kv1 subtypes, with the goal of elevating neuronal excitability and improving neurotransmission in certain diseases.

The Drosophila Shaker K(+) channel is the first cloned voltage-gated potassium channel and has, therefore, played an important role in structural and functional studies of those channels. While such a role is well justified for ion permeation, it is not clear whether this also extends to ligand binding. Despite the high degree of homology among Shaker and Kv1 channels, κ-conotoxin PVIIA (κ-PVIIA) binds to Shaker with high affinity but not to Kv1 channels. Here we address this issue by studying binding of κ-PVIIA to Shaker and Kv1 channels using molecular dynamics (MD) simulations. The structures of the channel-toxin complexes are constructed via docking and refinement with MD. The binding mode of each complex is characterized and compared to available mutagenesis data to validate the complex models. The potential of mean force for dissociation of the Shaker-κ-PVIIA complex is calculated from umbrella sampling MD simulations, and the corresponding binding free energy is determined, which provides further validation of the complex structure. Comparison of the Shaker and Kv1 complex models shows that a few mutations in the turret and extended regions are sufficient to abolish the observed sensitivity of Shaker to κ-PVIIA. This study demonstrates that Shaker is not always a good model for Kv1 channels for ligand binding. It also provides insights into the binding of the toxin to potassium channels that will be useful for improving affinity and selectivity properties of Kv1 channels.

The pharmacology and regulation of Transient Receptor Potential Ankyrin 1 (TRPA1) ion channel activity is intricate due to the physiological function as an integrator of multiple chemical, mechanical, and temperature stimuli as well as differences in species pharmacology. In this study, we describe and compare the current inhibition efficacy of human TRPA1 on three different TRPA1 antagonists. We used a homology model of TRPA1 based on Kv1.2 to select pore vestibule residues available for interaction with ligands entering the vestibule. Site-directed mutation constructs were expressed in Xenopus oocytes and their functionality and pharmacology assessed to support and improve our homology model. Based on the functional pharmacology results we propose an antagonist-binding site in the vestibule of the TRPA1 ion channel. We use the results to describe the proposed intravestibular ligand-binding site in TRPA1 in detail. Based on the single site substitutions, we designed a human TRPA1 receptor by substituting several residues in the vestibule and adjacent regions from the rat receptor to address and explain observed species pharmacology differences. In parallel, the lack of effect on HC-030031 inhibition by the vestibule substitutions suggests that this molecule interacts with TRPA1 via a binding site not situated in the vestibule., (Copyright © 2013 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Voltage-sensitive potassium ion channels are essential for life, but the molecular basis of their ion conduction is not well understood. In particular, the impact of ion concentration on ion conduction has not been fully studied. We performed several micro-second molecular dynamics simulations of the pore domain of the Kv1.2 potassium channel in KCl solution at four different ion concentrations, and scrutinized each of the conduction events, based on graphical representations of the simulation trajectories. As a result, we observed that the conduction mechanism switched with different ion concentrations: at high ion concentrations, potassium conduction occurred by Hodgkin and Keynes' knock-on mechanism, where the association of an incoming ion with the channel is tightly coupled with the dissociation of an outgoing ion, in a one-step manner. On the other hand, at low ion concentrations, ions mainly permeated by a two-step association/dissociation mechanism, in which the association and dissociation of ions were not coupled, and occurred in two distinct steps. We also found that this switch was triggered by the facilitated association of an ion from the intracellular side within the channel pore and by the delayed dissociation of the outermost ion, as the ion concentration increased.

Developing an understanding of the mechanism of voltage-gated ion channels in molecular terms requires knowledge of the structure of the active and resting conformations. Although the active-state conformation is known from x-ray structures, an atomic resolution structure of a voltage-dependent ion channel in the resting state is not currently available. This has motivated various efforts at using computational modeling methods and molecular dynamics (MD) simulations to provide the missing information. A comparison of recent computational results reveals an emerging consensus on voltage-dependent gating from computational modeling and MD simulations. This progress is highlighted in the broad context of preexisting work about voltage-gated channels.

The open state of voltage-gated potassium (Kv) channels is associated with an increased stability relative to the pre-open closed states and is reflected by a slowing of OFF gating currents after channel opening. The basis for this stabilization is usually assigned to intrinsic structural features of the open pore. We have studied the gating currents of Kv1.2 channels and found that the stabilization of the open state is instead conferred largely by the presence of cations occupying the inner cavity of the channel. Large impermeant intracellular cations such as N-methyl-d-glucamine (NMG(+)) and tetraethylammonium cause severe slowing of channel closure and gating currents, whereas the smaller cation, Cs(+), displays a more moderate effect on voltage sensor return. A nonconducting mutant also displays significant open state stabilization in the presence of intracellular K(+), suggesting that K(+) ions in the intracellular cavity also slow pore closure. A mutation in the S6 segment used previously to enlarge the inner cavity (Kv1.2-I402C) relieves the slowing of OFF gating currents in the presence of the large NMG(+) ion, suggesting that the interaction site for stabilizing ions resides within the inner cavity and creates an energetic barrier to pore closure. The physiological significance of ionic occupation of the inner cavity is underscored by the threefold slowing of ionic current deactivation in the wild-type channel compared with Kv1.2-I402C. The data suggest that internal ions, including physiological concentrations of K(+), allosterically regulate the deactivation kinetics of the Kv1.2 channel by impairing pore closure and limiting the return of voltage sensors. This may represent a primary mechanism by which Kv channel deactivation kinetics is linked to ion permeation and reveals a novel role for channel inner cavity residues to indirectly regulate voltage sensor dynamics.

Voltage-dependent potassium (Kv) channels provide the repolarizing power that shapes the action potential duration and helps control the firing frequency of neurons. The K(+) permeation through the channel pore is controlled by an intracellularly located bundle-crossing (BC) gate that communicates with the voltage-sensing domains (VSDs). During prolonged membrane depolarizations, most Kv channels display C-type inactivation that halts K(+) conduction through constriction of the K(+) selectivity filter. Besides triggering C-type inactivation, we show that in Shaker and Kv1.2 channels (expressed in Xenopus laevis oocytes), prolonged membrane depolarizations also slow down the kinetics of VSD deactivation and BC gate closure during the subsequent membrane repolarization. Measurements of deactivating gating currents (reporting VSD movement) and ionic currents (BC gate status) showed that the kinetics of both slowed down in two distinct phases with increasing duration of the depolarizing prepulse. The biphasic slowing in VSD deactivation and BC gate closure was strongly correlated in time and magnitude. Simultaneous recordings of ionic currents and fluorescence from a probe tracking VSD movement in Shaker directly demonstrated that both processes were synchronized. Whereas the first slowing originates from a stabilization imposed by BC gate opening, the subsequent slowing reflects the rearrangement of the VSD toward its relaxed state (relaxation). The VSD relaxation was observed in the Ciona intestinalis voltage-sensitive phosphatase and in its isolated VSD. Collectively, our results show that the VSD relaxation is not kinetically related to C-type inactivation and is an intrinsic property of the VSD. We propose VSD relaxation as a general mechanism for depolarization-induced slowing of BC gate closure that may enable Kv1.2 channels to modulate the firing frequency of neurons based on the depolarization history.

The human ether-a-go-go related gene 1 (hERG1) K ion channel is a key element for the rapid component of the delayed rectified potassium current in cardiac myocytes. Since there are no crystal structures for hERG channels, creation and validation of its reliable atomistic models have been key targets in molecular cardiology for the past decade. In this study, we developed and vigorously validated models for open, closed, and open-inactivated states of hERG1 using a multistep protocol. The conserved elements were derived using multiple-template homology modeling utilizing available structures for Kv1.2, Kv1.2/2.1 chimera, and KcsA channels. Then missing elements were modeled with the ROSETTA De Novo protein-designing suite and further refined with all-atom molecular dynamics simulations. The final ensemble of models was evaluated for consistency to the reported experimental data from biochemical, biophysical, and electrophysiological studies. The closed state models were cross-validated against available experimental data on toxin footprinting with protein-protein docking using hERG state-selective toxin BeKm-1. Poisson-Boltzmann calculations were performed to determine gating charge and compare it to electrophysiological measurements. The validated structures offered us a unique chance to assess molecular mechanisms of state-dependent drug binding in three different states of the channel.

Fullerene C(70) is known to partition into lipid membranes and change their physical properties. Together with gallic acid (GA), C(70) induces cell contraction and cell death. How C(70) and GA-induced perturbations of lipid membranes affect cellular function and membrane protein activity is not understood, though. Meanwhile, fullerene is also known to interfere with the activity of potassium channel proteins, but the mechanisms of protein inhibition are not known. Here we consider the possibility that membrane protein function would be inhibited by C(70) and/or GA through direct contact or through lipid-mediated interactions. To this end, we use microsecond time scale atomistic simulations to explore (a) modifications of membrane properties in the presence of C(70) and/or GA, and (b) the possible conformational changes in Kv1.2, a voltage-gated potassium channel, upon exposure to C(70), or GA, or both. C(70) is found to have an observable effect on structural and elastic properties of protein-free membranes, while the effects of GA on the membrane are less evident. Fullerene–GA interaction is strong and affects significantly the partitioning of C(70) in the membrane, stabilizing C(70) in the aqueous phase. When Kv1.2 is exposed to the solutes, only small conformational changes are observed on the microsecond time scale – comparable to the fluctuations observed in the absence of any solute. Blocking of the channel entrance is not observed, as fullerene binds mainly to hydrophobic residues, both in the water-exposed loops and in the transmembrane helices. The tilt angle of transmembrane helices in the voltage-sensing domain appears to be affected by direct contact with fullerene, but a generic effect due to the small increase in membrane thickness might also play a role. A small rotation of the S3 and S4 helices in the voltage-sensing domain is noticed when C(70) is embedded in the membrane. The interpretation of the observed conformational changes is not straightforward due to the associated time scales, which are difficult to sample with state-of-the-art computing resources. We cannot exclude that both membrane-mediated interactions and specific protein–solute interactions affect the conformation of the protein.

Voltage sensing confers conversion of a change in membrane potential to signaling activities underlying the physiological processes. For an ion channel, voltage sensitivity is usually experimentally measured by fitting electrophysiological data to Boltzmann distributions. In our study, a two-state model of the ion channel and equilibrium statistical mechanics principle were used to test the hypothesis of empirically calculating the overall voltage sensitivity of an ion channel on the basis of its closed and open conformations, and determine the contribution of individual residues to the voltage sensing. We examined the theoretical paradigm by performing experimental measurements with Kv1.2 channel and a series of mutants. The correlation between the calculated values and the experimental values is at respective level, R(2) = 0.73. Our report therefore provides in silico prediction of key conformations and has identified additional residues critical for voltage sensing., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The mechanism of ion channel voltage gating-how channels open and close in response to voltage changes-has been debated since Hodgkin and Huxley's seminal discovery that the crux of nerve conduction is ion flow across cellular membranes. Using all-atom molecular dynamics simulations, we show how a voltage-gated potassium channel (KV) switches between activated and deactivated states. On deactivation, pore hydrophobic collapse rapidly halts ion flow. Subsequent voltage-sensing domain (VSD) relaxation, including inward, 15-angstrom S4-helix motion, completes the transition. On activation, outward S4 motion tightens the VSD-pore linker, perturbing linker-S6-helix packing. Fluctuations allow water, then potassium ions, to reenter the pore; linker-S6 repacking stabilizes the open pore. We propose a mechanistic model for the sodium/potassium/calcium voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.

Understanding the detailed mechanism of the activation of voltage-gated ion channels has been a problem of great current interest. Reliable molecular simulations of voltage effects present a major challenge because meaningful converging microscopic simulations are not yet available and macroscopic treatments involve major uncertainties regarding the dielectric constant used and other key features. The current work has overcome some of the above challenges by using our recently developed coarse-grained (CG) model in simulating the activation of the Kv1.2 channel. The CG model has allowed us to explore problems that cannot be addressed at present by fully microscopic simulations, while providing insights on some features that are not usually considered in continuum models, including the distribution of the electrolytes between the membrane and the electrodes during the activation process and thus the nature of the gating current. Furthermore, the clear connection to microscopic descriptions combined with the power of CG modeling offers a powerful tool for exploring the energy balance between the protein conformational energy and the interaction with the external potential in voltage-activated channels. Our simulations have reproduced the observed experimental trend of the gating charge and, most significantly, the correct trend in the free energies, where the closed channel is more stable at negative potential and the open channel is more stable at positive potential. Moreover, we provide a unique view of the activation landscape and the time dependence of the activation process.

Eighteen new limonoids, chubularisins A-R (1-18), along with eleven known analogues, were isolated from the stem bark of Chukrasia tabularis. The structures of 1-18 were elucidated on the basis of spectroscopic data and chemical evidence. Compound 1 represented the first example of 8,9,12-orthoester of phragmalin limonoids. Interestingly, compounds 4, 8, and 22 exhibited potent and selective inhibition against the delayed rectifier (I(K)) K(+) current with IC(50) values of 0.61, 2.03, and 2.15 μM, respectively.

The voltage-gated proton channel (Hv1) is homologous to the voltage-sensing domain (VSD) of voltage-gated potassium (Kv) channels but lacks a separate pore domain. The Hv1 monomer has dual functions: it gates the proton current and also serves as the proton conduction pathway. To gain insight into the structure and dynamics of the yet unresolved proton permeation pathway, we performed all-atom molecular dynamics simulations of two different Hv1 homology models in a lipid bilayer in excess water. The structure of the Kv1.2-Kv2.1 paddle-chimera VSD was used as template to generate both models, but they differ in the sequence alignment of the S4 segment. In both models, we observe a water wire that extends through the membrane, whereas the corresponding region is dry in simulations of the Kv1.2-Kv2.1 paddle-chimera. We find that the kinetic stability of the water wire is dependent upon the identity and location of the residues lining the permeation pathway, in particular, the S4 arginines. A measurement of water transport kinetics indicates that the water wire is a relatively static feature of the permeation pathway. Taken together, our results suggest that proton conduction in Hv1 may occur via Grotthuss hopping along a robust water wire, with exchange of water molecules between inner and outer ends of the permeation pathway minimized by specific water-protein interactions. This article is part of a Special Issue entitled: Membrane protein structure and function., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The 34-residue polypeptide maurotoxin (MTx) isolated from scorpion venoms selectively inhibits the current of the voltage-gated potassium channel Kv1.2 by occluding the ion conduction pathway. Here using molecular dynamics simulation as a docking method, the binding modes of MTx to three closely related channels (Kv1.1, Kv1.2 and Kv1.3) are examined. We show that MTx forms more favorable electrostatic interactions with the outer vestibule of Kv1.2 compared to Kv1.1 and Kv1.3, consistent with the selectivity of MTx for Kv1.2 over Kv1.1 and Kv1.3 observed experimentally. One salt bridge in the bound complex of MTx-Kv1.2 forms and breaks in a simulation period of 20 ns, suggesting the dynamic nature of toxin-channel interactions. The toxin selectivity likely arises from the differences in the shape of the channel outer vestibule, giving rise to distinct orientations of MTx on block. Potential of mean force calculations show that MTx blocks Kv1.1, Kv1.2 and Kv1.3 with an IC(50) value of 6 µM, 0.6 nM and 18 µM, respectively.

Precise localization of axonal ion channels is crucial for proper electrical and chemical functions of axons. In myelinated axons, Kv1 (Shaker) voltage-gated potassium (Kv) channels are clustered in the juxtaparanodal regions flanking the node of Ranvier. The clustering can be disrupted by deletion of various proteins in mice, including contactin-associated protein-like 2 (Caspr2) and transient axonal glycoprotein-1 (TAG-1), a glycosylphosphatidylinositol-anchored cell adhesion molecule. However, the mechanism and function of Kv1 juxtaparanodal clustering remain unclear. Here, using a new myelin coculture of hippocampal neurons and oligodendrocytes, we report that tyrosine phosphorylation plays a critical role in TAG-1-mediated clustering of axonal Kv1.2 channels. In the coculture, myelin specifically ensheathed axons but not dendrites of hippocampal neurons and clustered endogenous axonal Kv1.2 into internodes. The trans-homophilic interaction of TAG-1 was sufficient to position Kv1.2 clusters on axonal membranes in a neuron/HEK293 coculture. Mutating a tyrosine residue (Tyr⁴⁵⁸) in the Kv1.2 C terminus or blocking tyrosine phosphorylation disrupted myelin- and TAG-1-mediated clustering of axonal Kv1.2. Furthermore, Kv1.2 voltage dependence and activation threshold were reduced by TAG-1 coexpression. This effect was eliminated by the Tyr⁴⁵⁸ mutation or by cholesterol depletion. Taken together, our studies suggest that myelin regulates both trafficking and activity of Kv1 channels along hippocampal axons through TAG-1.

The response of a membrane-bound Kv1.2 ion channel to an applied transmembrane potential has been studied using molecular dynamics simulations. Channel deactivation is shown to involve three intermediate states of the voltage sensor domain (VSD), and concomitant movement of helix S4 charges 10-15 Å along the bilayer normal; the latter being enabled by zipper-like sequential pairing of S4 basic residues with neighboring VSD acidic residues and membrane-lipid head groups. During the observed sequential transitions S4 basic residues pass through the recently discovered charge transfer center with its conserved phenylalanine residue, F(233). Analysis indicates that the local electric field within the VSD is focused near the F(233) residue and that it remains essentially unaltered during the entire process. Overall, the present computations provide an atomistic description of VSD response to hyperpolarization, add support to the sliding helix model, and capture essential features inferred from a variety of recent experiments.

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K(+) ions in various configurations. For both channels, we find that a pair of K(+) ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K(+) ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Background: Scorpion toxins are invaluable tools for ion channel research and are potential drugs for human channelopathies. However, it is still an open task to determine the molecular basis underlying the diverse interactions between toxin peptides and ion channels. The inhibitory peptide Maurotoxin (MTX) recognized the distantly related IK(Ca) and Kv1.2 channel with approximately the same potency and using the same functional residues, their differential binding mechanism remain elusive. In this study, we applied computational methods to explore the differential binding modes of MTX to Kv1.2 and IK(Ca) channels, which would help to understand the diversity of channel-toxin interactions and accelerate the toxin-based drug design., Results: A reasonably stable MTX-IK(Ca) complex was obtained by combining various computational methods and by in-depth comparison with the previous model of the MTX-Kv1.2 complex. Similarly, MTX adopted the β-sheet structure as the interacting surface for binding both channels, with Lys23 occluding the pore. In contrast, the other critical residues Lys27, Lys30, and Tyr32 of MTX adopted distinct interactions when associating with the IK(Ca) channel. In addition, the residues Gln229, Ala230, Ala233, and Thr234 on the IK(Ca) channel turret formed polar and non-polar interactions with MTX, whereas the turret of Kv1.2 was almost not involved in recognizing MTX. In all, the pairs of interacting residues on MTX and the IK(Ca) channel of the bound complex indicated that electrostatic and Van der Waal interactions contributed equally to the formation of a stable MTX-IK(Ca) complex, in contrast to the MTX-Kv1.2 binding that is dominantly mediated by electrostatic forces., Conclusions: Despite sharing similar pharmacological profiles toward both IK(Ca) and Kv1.2 channels, MTX adopted totally diverging modes in the two association processes. All the molecular information unveiled here could not only offer a better understanding about the structural differences between the IK(Ca) and Kv1.2 channels, but also provide novel structural clues that will help in the designing of more selective molecular probes to discriminate between these two channels.

K(V)1.2 is a potassium channel protein whose electrophysiological properties, including an oxygen sensing response, are modulated by auxiliary beta subunits. The beta subunits are homologous to oxidoreductases, supporting a hypothesis that the coupling of the two subunits helps connect the redox state of the cell to the electrical activity of the membrane. However the exact mechanism of the coupling has not been discovered to-date. We apply evolutionary correlation analysis to infer previously unknown components of the interaction network regulating the response to hypoxia of the K(V)1.2/β₂ complex. Briefly, evolutionary correlation analysis involves finding correlated amino acid substitutions in functionally equivalent proteins (for both subunits) across a range of species. The method thus depends on a reliable method of inferring functionally equivalent (orthologous) proteins in different species, which method we describe in a paper recently published in PLoS ONE. One key finding is the characterization of a network of motif interactions between the α and the β subunits. By significance testing, we show that the likelihood of the correlations shown in the K(V)1.2-β two interaction motifs arising by chance is less than 0.0003, which shows that the correlations are statistically highly significant. We therefore believe that the correlations are likely to be biologically relevant. Other major findings are correlations between specific motifs in the K(V)1.2-β₂ complex and motifs in other proteins such as RACK1 and Eif36sip, which in turn are connected to the hypoxia response. Our paper combines this correlation with the literature evidence on potassium channel inactivation and hypoxia response to identify specific motifs to serve as experimental targets for studies focused on this response. This work aims to add to our general understanding of two major issues in ion channel science: 1) how multi-protein complexes including ion channels function in coordinated fashion, and 2) how ion channels mediate the conversation between the intracellular and extracellular environments. We also aim to apply to ion channel science the principle that domain-domain interactions in proteins can be inferred from correlation of amino-acid substitutions in sets of functionally equivalent proteins.

Shaker-related Kv1 channels contain four channel-forming alpha subunits. Subfamily member Kv1.1 often occurs oligomerized with Kv1.2 alpha subunits in synaptic membranes, and so information was sought on the influence of their positions within tetramers on the channels' properties. Kv1.1 and 1.2 alpha genes were tandem linked in various arrangements, followed by expression as single-chain proteins in mammalian cells. As some concatenations reported previously seemed not to reliably position Kv1 subunits in their assemblies, the identity of expressed channels was methodically evaluated. Surface protein, isolated by biotinylation of intact transiently transfected HEK-293 cells, gave Kv1.1/1.2 reactivity on immunoblots with electrophoretic mobilities corresponding to full-length concatenated tetramers. There was no evidence of protein degradation, indicating that concatemers were delivered intact to the plasmalemma. Constructs with like genes adjacent (Kv1.1-1.1-1.2-1.2 or Kv1.2-1.2-1.1-1.1) yielded delayed-rectifying, voltage-dependent K(+) currents with activation parameters and inactivation kinetics slightly different from the diagonally positioned genes (Kv1.1-1.2-1.1-1.2 or 1.2-1.1-1.2-1.1). Pore-blocking petidergic toxins, alpha dendrotoxin, agitoxin-1, tityustoxin-Kalpha, and kaliotoxin, were unable to distinguish between the adjacent and diagonal concatamers. Unprecedentedly, external application of the pore-blocker tetraethylammonium (TEA) differentially inhibited the adjacent versus diagonal subunit arrangements, with diagonal constructs having enhanced susceptibility. Concatenation did not directly alter the sensitivities of homomeric Kv1.1 or 1.2 channels to TEA or the toxins. TEA inhibition of currents generated by channels made up from dimers (Kv1.1-1.2 and/or Kv1.2-1.1) was similar to the adjacently arranged constructs. These collective findings indicate that assembly of alpha subunits can be directed by this optimized concatenation, and that subunit arrangement in heteromeric Kv channels affects TEA affinity.

Concerted actions of various ion channels that are precisely targeted along axons are crucial for action potential initiation and propagation, and neurotransmitter release. However, the dynamics of channel protein transport in axons remain unknown. Here, using time-lapse imaging, we found fluorescently tagged Kv1.2 voltage-gated K(+) channels (YFP-Kv1.2) moved bi-directionally in discrete puncta along hippocampal axons. Expressing Kvbeta2, a Kv1 accessory subunit, markedly increased the velocity, the travel distance, and the percentage of moving time of these puncta in both anterograde and retrograde directions. Suppressing the Kvbeta2-associated protein, plus-end binding protein EB1 or kinesin II/KIF3A, by siRNA, significantly decreased the velocity of YFP-Kv1.2 moving puncta in both directions. Kvbeta2 mutants with disrupted either Kv1.2-Kvbeta2 binding or Kvbeta2-EB1 binding failed to increase the velocity of YFP-Kv1.2 puncta, confirming a central role of Kvbeta2. Furthermore, fluorescently tagged Kv1.2 and Kvbeta2 co-moved along axons. Surprisingly, when co-moving with Kv1.2 and Kvbeta2, EB1 appeared to travel markedly faster than its plus-end tracking. Finally, using fission yeast S. pombe expressing YFP-fusion proteins as reference standards to calibrate our microscope, we estimated the numbers of YFP-Kv1.2 tetramers in axonal puncta. Taken together, our results suggest that proper amounts of Kv1 channels and their associated proteins are required for efficient transport of Kv1 channel proteins along axons.

Voltage-dependent potassium channels (Kv) are homotetramers composed of four voltage sensors and one pore domain. Because of high-level structural flexibility, the first mammalian Kv structure, Kv1.2 at 2.9 A, has about 37% molecular mass of the transmembrane portion not resolved. In this study, by applying a novel normal-mode-based X-ray crystallographic refinement method to the original diffraction data and structural model, we established the structure of full-length Kv1.2 in its native form. This structure offers mechanistic insights into voltage sensing. Particularly, it shows a hydrophobic layer of about 10 A at the midpoint of the membrane bilayer, which is likely the molecular basis for the observed "focused electric field" of Kv1.2 between the internal and external solutions. This work also demonstrated the potential of the refinement method in bringing up large chunks of missing densities, thus beneficial to structural refinement of many difficult systems.

Voltage-gated potassium (Kv) channels, such as Kv1.2, are involved in the generation and propagation of action potentials. The Kv channel is a homotetramer, and each monomer is composed of a voltage-sensing domain (VSD) and a pore domain (PD). We analyzed the fluctuations of a model structure of Kv1.2 using elastic network models. The analysis suggested a network of coupled fluctuations of eight rigid structural units and seven hinges that may control the transition between the active and inactive states of the channel. For the most part, the network is composed of amino acids that are known to affect channel activity. The results suggested allosteric interactions and cooperativity between the subunits in the coupling between the motion of the VSD and the selectivity filter of the PD, in accordance with recent empirical data. There are no direct contacts between the VSDs of the four subunits, and the contacts between these and the PDs are loose, suggesting that the VSDs are capable of functioning independently. Indeed, they manifest many inherent fluctuations that are decoupled from the rest of the structure. In general, the analysis suggests that the two domains contribute to the channel function both individually and cooperatively., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Voltage sensors regulate the conformations of voltage-dependent ion channels and enzymes. Their nearly switchlike response as a function of membrane voltage comes from the movement of positively charged amino acids, arginine or lysine, across the membrane field. We used mutations with natural and unnatural amino acids, electrophysiological recordings, and x-ray crystallography to identify a charge transfer center in voltage sensors that facilitates this movement. This center consists of a rigid cyclic "cap" and two negatively charged amino acids to interact with a positive charge. Specific mutations induce a preference for lysine relative to arginine. By placing lysine at specific locations, the voltage sensor can be stabilized in different conformations, which enables a dissection of voltage sensor movements and their relation to ion channel opening.

We present the first atomic-resolution observations of permeation and gating in a K(+) channel, based on molecular dynamics simulations of the Kv1.2 pore domain. Analysis of hundreds of simulated permeation events revealed a detailed conduction mechanism, resembling the Hodgkin-Keynes "knock-on" model, in which translocation of two selectivity filter-bound ions is driven by a third ion; formation of this knock-on intermediate is rate determining. In addition, at reverse or zero voltages, we observed pore closure by a novel "hydrophobic gating" mechanism: A dewetting transition of the hydrophobic pore cavity-fastest when K(+) was not bound in selectivity filter sites nearest the cavity-caused the open, conducting pore to collapse into a closed, nonconducting conformation. Such pore closure corroborates the idea that voltage sensors can act to prevent pore collapse into the intrinsically more stable, closed conformation, and it further suggests that molecular-scale dewetting facilitates a specific biological function: K(+) channel gating. Existing experimental data support our hypothesis that hydrophobic gating may be a fundamental principle underlying the gating of voltage-sensitive K(+) channels. We suggest that hydrophobic gating explains, in part, why diverse ion channels conserve hydrophobic pore cavities, and we speculate that modulation of cavity hydration could enable structural determination of both open and closed channels.

The transmembrane part of the S6 inner helix of the Kv1.2 potassium channel is a pivotal part in sustaining channel activity. However, the role of its extreme C-terminal end, which is located on the cytoplasmic side of the channel, is largely unknown. Here, we investigated the role of the extreme C-terminal end of the S6 inner helix (the HRET region) by constructing a series of C-terminal-truncated mutations related to this region in the C-terminal-truncated Kv1.2 channel. Mutations on Thr421 or Glu420 significantly altered the activation of the truncated channel. Mutations on Arg419 demonstrated that neutral or basic, but not acidic amino acid, is essential at the position for the truncated channel activation, and no functional channel was observed when the channel was truncated from His418. Hence, our results indicate that the extreme C-terminal end of the S6 inner helix plays an important regulatory role in the activation of the C-terminal-truncated Kv1.2 channel.

The Kv1.4 potassium channel is reported to exhibit higher cell surface expression than the Kv1.1 potassium channel when expressed as a homomer in cell lines. Kv1.4 also shows highly efficient trans-Golgi glycosylation whereas Kv1.1 is not glycosylated. The surface expression and glycosylation of Kv1.2 is intermediate between those of Kv1.1 and Kv1.4. Amino acid determinants controlling the surface expression of Kv1 channels were localized to the highly conserved pore region and both positive and negative determinants of Kv1.1 and Kv1.4 trafficking have been reported. In this study, we analyzed the effect of substituting amino acids in the pore region of Kv1.2 with the corresponding amino acid present in Kv1.1 or Kv1.4 on glycosylation and trafficking of Kv1.2. Mutations in the outer pore region of Kv1.2 of Arg(354) to Pro (corresponding to Kv1.4) and to Ala (corresponding to Kv1.1) enhanced and reduced, respectively, cell surface expression of Kv1.2. Mutations in a different outer pore region of Val(381) to Lys (Kv1.4) and Tyr (Kv1.1) both reduced the cell surface expression. In contrast, mutation in the deep pore region of Ser(371) to Thr (Kv1.4) markedly enhanced cell surface expression. These results suggest that the cell surface expression of Kv1.2 is regulated by specific amino acids in the pore region in a similar manner to Kv1.1 and Kv1.4, and that the cell surface expression of Kv1.2, a channel intermediate between Kv1.1 and Kv1.4, can be attributed to these specific residues.

Recent structures of the potassium channel provide an essential beginning point for explaining how the pore is gated between open and closed conformations by changes in membrane voltage. Yet, the molecular details of this process and the connections to transmembrane gradients are not understood. To begin addressing how changes within a membrane environment lead to the channel's ability to sense shifts in membrane voltage and to gate, we performed double-bilayer simulations of the Kv1.2 channel. These double-bilayer simulations enable us to simulate realistic voltage drops from resting potential conditions to depolarized conditions by changes in the bath conditions on each side of the bilayer. Our results show how the voltage sensor domain movement responds to differences in transmembrane potential. The initial voltage sensor domain movement, S4 in particular, is modulated by the gating charge response to changes in voltage and is initially stabilized by the lipid headgroups. We show this response is directly coupled to the initial stages of pore domain motion. Results presented here provide a molecular model for how the pre-gating process occurs in sequential steps: Gating charge response, movement and stabilization of the S4 voltage sensor domain, and movement near the base of the S5 region to close the pore domain.

Cytosolic K(V)beta1 subunits co-assemble with transmembrane K(V)1 channel alpha-subunits and have complex effects on channel function. Fast inactivation, the most obvious effect conferred, is due to fast open channel block resulting from the binding of the N-terminus within the inner mouth of the pore. K(V)beta1 subunits also slow current deactivation, enhance slow inactivation and shift channel activation to more negative voltages, but the mechanisms underlying these actions are not known. Here we use voltage clamp fluorimetry at sites near the extracellular end of the S4 helix, the channel's primary voltage sensor, in combination with voltage clamp electrophysiology, to independently track the movement of the S4 helix along with ionic current, and thus identify the structural and mechanistic means by which the K(V)beta1.2 subunit confers its actions on the K(V)1.2 channel. We show that the negative shift in current activation is not due to direct actions of K(V)beta1.2 on the S4 segment. Instead, this shift results from an apparent saturation of channel activation at depolarized potentials as the extent of open channel block by the K(V)beta1.2 N-terminus progressively increases. The return of fluorescence to baseline is slowed along with current deactivation. According to our data, this is due to an inability of the activation gate to close while the K(V)beta1.2 N-terminus occupies the pore and strong coupling of the gate with the S4 segment. Together with data from previous studies, our findings provide a complete and coherent picture of the functional and structural interactions between K(V)beta1.2 and K(V)1.2.

N-type inactivation of the Kv1.2 voltage-gated potassium channel is a process in which the N-terminal of the protein (its first 20 amino acids) binds to the open-channel surface, extends and occludes its pore. This process has been experimentally studied in both intact and ShBDelta6-46 channels in which the inactivating peptides are supplied in the bath solution. In this work we provide a qualitative description of N-type inactivation by simulating the random walk of charged inactivating peptides in the electrostatic field that originates from the charges present in the channel and in the cellular membrane. Our results give a deeper insight into the previously reported influence of electrostatics on the rate of N-type inactivation of ShBDelta6-46. We also show how the enchaining of the peptides, i.e., considering the intact form of the channel, influences the N-type inactivation with different charges of those peptides.

Voltage-gated K(+) channels undergo a voltage-dependent conductance change that plays a key role in modulating cellular excitability. While the Open state is captured in crystal structures of Kv1.2 and a chimeric Kv1.2/Kv2.1 channel, the Close state and the mechanism of this transition are still a subject of debate. Here, we propose a model based on mutagenesis combined with measurements of both ionic and gating currents which is consistent with the idea that the Open state is the default state, the energy of the electric field being used to keep the channel closed. Our model incorporates an 'Activated state' where the bulk of sensor movement is completed without channel opening. The model accounts for the well characterized electrophysiology of the 'V2' and 'ILT' mutations in Shaker, where sensor movement and channel opening occur over distinct voltage ranges. Moreover, the model proposes relatively small protein rearrangements in going from the Activated to the Open state, consistent with the rapid transitions observed in single channel records of Shaker type channels at zero millivolts.

The recently determined crystal structure of a chimeric Kv1.2-Kv2.1 Kv channel at 2.4 A resolution motivated this molecular-dynamics simulation study of the chimeric channel and its mutants embedded in a DPPC membrane. For the channel protein, we used two types of C-terminus: E+ and Eo. E+ contains, and Eo lacks, the EGEE residue quartet located distal to the S6 helix. For both E+ and Eo, the following trend was observed: When S4 helices were restrained at the same position as in the x-ray structure (S4high), the S6 gate remained open for 12 ns. The results were similar when the S4 helices were pulled downward 7 A (S4low). However, S4middle (or S4low) facilitated the S6 gate-narrowing for the following mutated channels (shown in order of increasing effect): 1), E395W; 2), E395W-F401A-F402A; and 3), E395W-F401A-F402A-V478W. The amino acid numbering system is that used for the Shaker channel. Even though all four subunits were set at S4low, S6 gate-narrowing was often brought about by movements of only two opposing S6 helices toward the central axis of the pore, resulting in a twofold symmetry-like structure. A free-energy profile analysis over the ion conduction pathway shows that the two opposing S6 helices whose peptide backbones are approximately 10.4 A distant from each other lead to an energetic barrier of approximately 25 kJ/mol. S6 movement was coupled with translocation of the S4-S5 linker toward the central axis of the same subunit, and the coupling was mediated by salt bridges formed between the inner (intracellular side) end of S4 and that of S6. Simulations in which S4 of only one subunit was pulled down to S4low showed that a weak intersubunit coordination is present for S5 movement, whereas the coupling between the S4-S5 linker and S6 is largely an intrasubunit one. In general, whereas subunit-based behavior appears to be dominant and to permit heteromeric conformations of the pore domain, direct intersubunit coupling of S5 or S6 is weak. Therefore, the "concerted transition" of the pore domain that has been predicted based on electrophysiological analyses is likely to be mediated mainly by the dual effects of S4 and the S4-S5 linker; these segments of one subunit can interact with both S5 of the same subunit and that of the adjacent subunit.

Potassium channels are among the core functional elements of life because they underpin essential cellular functions including excitability, homeostasis, and secretion. We present here a series of multivalent calix[4]arene ligands that bind to the surface of voltage-dependent potassium channels (K(v)1.x) in a reversible manner. Molecular modeling correctly predicts the best candidates with a conical C(4) symmetry for optimal binding, and the effects on channel function are assessed electrophysiologically. Reversible inhibition was observed, without noticeable damage of the oocytes, for tetraacylguanidinium or tetraarginine members of the series with small lower rim O-substituents. Apparent binding constants were in the low micromolar range and had Hill coefficients of 1, consistent with a single site of binding. Suppression of current amplitude was accompanied by a positive shift in the voltage dependence of gating and slowing of both voltage sensor motion and channel opening. These effects are in keeping with expectations for docking in the central pore and interaction with the pore domain "turret."

Vascular smooth muscle Kv1 delayed rectifier K+ channels (KDR) containing Kv1.2 control membrane potential and thereby regulate contractility. Vasodilatory agonists acting via protein kinase A (PKA) enhance vascule smooth muscle Kv1 activity, but the molecular basis of this regulation is uncertain. We characterized the role of a C-terminal phosphorylation site, Ser-449, in Kv1.2 expressed in HEK 293 cells by biochemical and electrophysiological methods. We found that 1) in vitro phosphorylation of Kv1.2 occurred exclusively at serine residues, 2) one major phosphopeptide that co-migrated with 449pSASTISK was generated by proteolysis of in vitro phosphorylated Kv1.2, 3) the peptide 445KKSRSASTISK exhibited stoichiometric phosphorylation by PKA in vitro, 4) matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectroscopy (MS) and MS/MS confirmed in vitro Ser-449 phosphorylation by PKA, 5) in situ phosphorylation at Ser-449 was detected in HEK 293 cells by MALDI-TOF MS followed by MS/MS. MIDAS (multiple reaction monitoring-initiated detection and sequencing) analysis revealed additional phosphorylated residues, Ser-440 and Ser-441, 6) in vitro 32P incorporation was significantly reduced in Kv1.2-S449A, Kv1.2-S449D, and Kv1.2-S440A/S441A/S449A mutant channels, but Kv1.2-S440A/S441A was identical to wild-type Kv1.2 (Kv1.2-WT), and 7) bath applied 8-Br-cAMP or dialysis with PKA catalytic subunit (cPKA) increased Kv1.2-WT but not Kv1.2-S449A current amplitude. cPKA increased Kv1.2-WT current in inside-out patches. Rp-CPT-cAMPS reduced Kv1.2-WT current, blocked the increase due to 8-Br-cAMP, but had no effect on Kv1.2-S449A. cPKA increased current due to double mutant Kv1.2-S440A/S441A but had no effect on Kv1.2-S449D or Kv1.2-S440A/S441A/S449A. We conclude that Ser-449 in Kv1.2 is a site of PKA phosphorylation and a potential molecular mechanism for Kv1-containing KDR channel modulation by agonists via PKA activation.

A one-on-one mapping of protein functionality across different species is a critical component of comparative analysis. This paper presents a heuristic algorithm for discovering the Most Likely Functional Counterparts (MoLFunCs) of a protein, based on simple concepts from network theory. A key feature of our algorithm is utilization of the user's knowledge to assign high confidence to selected functional identification. We show use of the algorithm to retrieve functional equivalents for 7 membrane proteins, from an exploration of almost 40 genomes form multiple online resources. We verify the functional equivalency of our dataset through a series of tests that include sequence, structure and function comparisons. Comparison is made to the OMA methodology, which also identifies one-on-one mapping between proteins from different species. Based on that comparison, we believe that incorporation of user's knowledge as a key aspect of the technique adds value to purely statistical formal methods.

Voltage-gated potassium (Kv) channels control the electrical excitability of neurons and muscles. Despite this key role, how these channels open and close or gate is not fully understood. Gating is usually attributed to the bending and straightening of pore-lining helices at glycine and proline residues. In this work we focused on the role of proline in the Pro-Val-Pro (PVP) motif of the inner S6 helix in the Kv1.2 channel. We started by developing a simple hinged-rod model to fully explore the configurational space of bent helices and we related these configurations to the degree of pore opening. We then carried out fully atomistic simulations of the S6 helices and compared these simulations to the hinged-rod model. Both methods suggest that Kv1 channels are not tightly closed when the inner helices are straight, unlike what is seen in the non-PVP containing channels KcsA and KirBac. These results invite the possibility that the S6 helices may be kinked when Kv1 channels are closed. Our simulations indicate that the wild-type helix adopts multiple spatially distinct configurations, which is consistent with its role in adopting a closed state and an open state. The two most dominant configurational basins correspond to a 6 A movement of the helix tail accompanied by the PVP region undergoing a local alpha-helix to 3(10)-helix transition. We explored how single point mutations affect the propensity of the S6 helix to adopt particular configurations. Interestingly, mutating the first proline, P405 (P473 in Shaker), to alanine completely removed the bistable nature of the S6 helix possibly explaining why this mutation compromises the channel. Next, we considered four other mutations in the area known to affect channel gating and we saw similarly dramatic changes to the helix's dynamics and range of motion. Our results suggest a possible mechanism of helix pore closure and they suggest differences in the closed state of glycine-only channels, like KcsA, and PVP containing channels.

Ryanodine receptor type 1 (RyR1) produces spatially and temporally defined Ca2+ signals in several cell types. How signals received in the cytoplasmic domain are transmitted to the ion gate and how the channel gates are unknown. We used EGTA or neuroactive PCB 95 to stabilize the full closed or open states of RyR1. Single-channel measurements in the presence of FKBP12 indicate that PCB 95 inverts the thermodynamic stability of RyR1 and locks it in a long-lived open state whose unitary current is indistinguishable from the native open state. We analyzed two datasets of 15,625 and 18,527 frozen-hydrated RyR1-FKBP12 particles in the closed and open conformations, respectively, by cryo-electron microscopy. Their corresponding three-dimensional structures at 10.2 A resolution refine the structure surrounding the ion pathway previously identified in the closed conformation: two right-handed bundles emerging from the putative ion gate (the cytoplasmic "inner branches" and the transmembrane "inner helices"). Furthermore, six of the identifiable transmembrane segments of RyR1 have similar organization to those of the mammalian Kv1.2 potassium channel. Upon gating, the distal cytoplasmic domains move towards the transmembrane domain while the central cytoplasmic domains move away from it, and also away from the 4-fold axis. Along the ion pathway, precise relocation of the inner helices and inner branches results in an approximately 4 A diameter increase of the ion gate. Whereas the inner helices of the K+ channels and of the RyR1 channel cross-correlate best with their corresponding open/closed states, the cytoplasmic inner branches, which are not observed in the K+ channels, appear to have at least as important a role as the inner helices for RyR1 gating. We propose a theoretical model whereby the inner helices, the inner branches, and the h1 densities together create an efficient novel gating mechanism for channel opening by relaxing two right-handed bundle structures along a common 4-fold axis., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

Although tertiary folding of whole protein domains is prohibited by the cramped dimensions of the ribosomal tunnel, dynamic tertiary interactions may permit folding of small elementary units within the tunnel. To probe this possibility, we used a beta-hairpin and an alpha-helical hairpin from the cytosolic N terminus of a voltage-gated potassium channel and determined a probability of folding for each at defined locations inside and outside the tunnel. Minimalist tertiary structures can form near the exit port of the tunnel, a region that provides an entropic window for initial exploration of local peptide conformations. Tertiary subdomains of the nascent peptide fold sequentially, but not independently, during translation. These studies offer an approach for diagnosing the molecular basis for folding defects that lead to protein malfunction and provide insight into the role of the ribosome during early potassium channel biogenesis.

The cytoplasmic C-terminus plays regulatory roles in the gating of many ion channels. However, lack of structural information on the C-terminus prevents the elucidation of how the C-terminal domain interacts with the gating machinery to exert its effects on the channel gating. In this report, we investigated the regulatory role of the C-terminus with functional study and structural modeling of a succession of C-terminal truncations of the Kv1.2 and Kv1.2(427)-KcsA(112-160) chimeric channels. Functional study demonstrated a length-dependent shift of the activation curves for the C-terminal truncations of the Kv1.2 channel. Structural modeling indicated that the C-terminus of one subunit could dynamically interact with the S4-S5 linker of a neighboring subunit and the probability of interaction was dependent on the length of the C-terminal truncated Kv1.2 channels. In contrast, no length-dependent shift of the activation curve and probability of interaction between C-terminus and the neighboring S4-S5 linker were observed for the truncations of the Kv1.2-KcsA chimeric channel, suggesting that the native C-terminus of the Kv1.2 channel is essential for the interaction. Furthermore, surface plasmon resonance measurements indicated that there is direct interaction between the C-terminal domain and the S4-S5 linker of the Kv1.2 channel. These results imply that the dynamic interaction of the C-terminus with the S4-S5 linker from a neighboring subunit of the Kv1.2 channel provides a mechanism for its C-terminus to regulate the channel activation.

Early transition events of the voltage sensor (VS) of Kv1.2 potassium channel embedded in a lipid membrane are triggered using full atomistic molecular dynamics (MD) simulations. When subject to an applied hyperpolarized transmembrane (TM) voltage, the VS undergoes conformational changes and reaches a stable kinetic intermediate state, beta', within 20 ns. The gating charge ( approximately 2e) associated with this fast transition results mainly from salt-bridge rearrangements involving negative charges in S2 and S3 and all but the two top residues R(294) and R(297) of S4. Interactions of the latter with phosphomoieties of the lipid head groups appear to stabilize the kinetic state beta'.

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins-not to mention numerous applications in drug design. Here, we present a full 1 micros atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120 degrees rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation ( approximately 35 degrees ) of the extracellular end of all S4 segments is present also in a reference 0.5 micros simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 3(10) helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4-lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5-1 micros). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations.

Voltage-gated potassium Kv1 channels have three extracellular linkers, the S1-S2, the S3-S4, and the S5-P. The S1-S2 is the only linker that has an N-glycan and it is at a conserved position on this linker on Kv1.1-Kv1.5 and Kv1.7 channels. We hypothesize that an N-glycan is found at only this position due to its effect on folding, trafficking, and/or function of these channels. To investigate this hypothesis, N-glycosylation sites were engineered at different positions on the extracellular linkers of Kv1.2 to determine the effects of N-glycans on channel surface protein expression and function. Our data suggest that for Kv1 channels, (1) placing an N-glycan at non-native positions on the S1-S2 linker decreased cell surface protein expression but the N-glycan still affected function similarly as if it were at its native position, (2) placing a non-native N-glycan on the S3-S4 linker significantly altered function, and (3) placing a non-native N-glycan on the S5-P linker disrupted both trafficking and function. We suggest that Kv1 channels have an N-glycan at a conserved position on only the S1-S2 linker to overcome the constraints for proper folding, trafficking, and function that appear to occur if the N-glycan is moved from this position.

We performed 45 ns atomistic force field based molecular dynamics (MD) simulations to explore the structure and dynamics of the Kvl.2 voltage-dependent potassium ion channel in both open and closed state conformations. The Kv1.2 crystal structure (PDB 2A79) based open state model and homologically derived closed state model are embedded in a hydrated lipid membrane. We present a detailed analysis of the protein stability, the environment of the gating-charge-carrying residues, the salt bridge within the voltage sensor, S4 segment helix movement, helix interaction between S4-S5 linker and S6, interaction between subunits, and interaction between protein and lipid membrane. The four subunits of the channel lost the symmetry from the starting structure during the simulation, especially the voltage sensors. According to our comparative analysis of the open and closed state conformations, S4 in our simulation behaves more near the screw helix model except for the tilt action. The S4 segment azimuthal rotation angle difference between two conformations shifts about 60 degrees after being embedded in a membrane environment. The S4-S5 linker helix remains stable. The contact area between the S4 and S5 from the adjacent subunit increases after transition to the closed state from the open state. The current investigation provided valuable information to our understanding of the structure and dynamics of membrane-associated helices and the gating mechanism of the voltage-dependent potassium ion channel of the Kv channel.

The Shaker family voltage-dependent potassium channels (Kv1) are expressed in a wide variety of cells and are essential for cellular excitability. In humans, loss-of-function mutations of Kv1 channels lead to hyperexcitability and are directly linked to episodic ataxia and atrial fibrillation. All Kv1 channels assemble with beta subunits (Kv betas), and certain Kv betas, for example Kv beta 1, have an N-terminal segment that closes the channel by the N-type inactivation mechanism. In principle, dissociation of Kv beta 1, although never reported, should eliminate inactivation and thus potentiate Kv1 current. We found that cortisone increases rat Kv1 channel activity by binding to Kv beta 1. A crystal structure of the Kv beta-cortisone complex was solved to 1.82-A resolution and revealed novel cortisone binding sites. Further studies demonstrated that cortisone promotes dissociation of Kv beta. The new mode of channel modulation may be explored by native or synthetic ligands to fine-tune cellular excitability.

Voltage-dependent K(+) (Kv) channels play key roles in shaping electrical signaling in both excitable and nonexcitable cells. These channels open and close in response to the voltage changes across the cell membrane. Many studies have been carried out in order to understand the voltage-sensing mechanism. Our laboratory recently determined the atomic structures of a mammalian Kv channel Kv1.2 and a mutant of Kv1.2 named the 'paddle chimera' channel, in which the voltage sensor paddle was transferred from Kv2.1 to Kv1.2. These two structures provide atomic descriptions of voltage-dependent channels with unprecedented clarity. Until now, the functional integrity of these two channels biosynthesized in yeast cells has not been assessed. Here, we report the electrophysiological and pharmacological properties of Kv1.2 and the paddle chimera channels in planar lipid bilayers. We demonstrate that Pichia yeast produce 'normally functioning' mammalian Kv channels with qualitatively similar features to the Shaker K(+) channel in the absence of the N-terminal inactivation gate and that the paddle chimera mutant channel functions as well as Kv1.2. We find, however, that in several respects, the Kv1.2 channel exhibits functional properties that are distinct from Kv1.2 channels reported in the literature.

The transient receptor potential vanilloid subtype 1 (TRPV1) is a member of the TRP family gated by vanilloids, heat, and protons. Structurally, TRPV1 subunits have a modular architecture underlying different functionalities, namely stimuli recognition, channel gating, ion selectivity, subunit oligomerization, and regulation by intracellular signaling molecules. Considering modular organization and recent structural information in the ion channel field, we have modeled a full-length TRPV1 by assembly of its major modules: the cytosolic N-terminal, C-terminal, and membrane-spanning region. For N-terminal, we used the ankyrin repeat structure fused with the N-end segment. The membrane domain was modeled with the structure of the eukaryotic, voltage-gated Kv1.2 K+ channel. The C-terminus was cast using the coordinates of HCN channels. The extensive structure-function data available for TRPV1 was used to validate the models in terms of the location of molecular determinants of function in the structure. Additionally, the current information allowed the modeling of the vanilloid receptor in the closed and desensitized states. The closed state shows the N-terminal module highly exposed and accessible to adenosine triphosphate and the C-terminal accessible to phosphoinositides. In contrast, the desensitized state depicts the N-terminal and C-terminal modules close together, compatible with an interaction mediated by Ca2+ -calmodulin complex. These models identify potential previously unrecognized intra- and interdomain interactions that may play an important functional role. Although the molecular models should be taken with caution, they provide a helpful tool that yields testable hypothesis that further our understanding on ion channels work in terms of underlying protein structure.

In voltage-gated K(+) channels (Kv), membrane depolarization promotes a structural reorganization of each of the four voltage sensor domains surrounding the conducting pore, inducing its opening. Although the crystal structure of Kv1.2 provided the first atomic resolution view of a eukaryotic Kv channel, several components of the voltage sensors remain poorly resolved. In particular, the position and orientation of the charged arginine side chains in the S4 transmembrane segments remain controversial. Here we investigate the proximity of S4 and the pore domain in functional Kv1.2 channels in a native membrane environment using electrophysiological analysis of intersubunit histidine metallic bridges formed between the first arginine of S4 (R294) and residues A351 or D352 of the pore domain. We show that histidine pairs are able to bind Zn(2+) or Cd(2+) with high affinity, demonstrating their close physical proximity. The results of molecular dynamics simulations, consistent with electrophysiological data, indicate that the position of the S4 helix in the functional open-activated state could be shifted by approximately 7-8 A and rotated counterclockwise by 37 degrees along its main axis relative to its position observed in the Kv1.2 x-ray structure. A structural model is provided for this conformation. The results further highlight the dynamic and flexible nature of the voltage sensor.