The intrinsic flexibility of the Kv voltage sensor and its implications for channel gating

Exploring Flexibility and Folding Patterns Throughout Time in Voltage Sensors

The voltage-sensing domain (VSD) is a module capable of responding to changes in the membrane potential through conformational changes and facilitating electromechanical coupling to open a pore gate, activate proton permeation pathways, or promote enzymatic activity in some membrane-anchored phosphatases. To carry out these functions, this module acts cooperatively through conformational changes. The VSD is formed by four transmembrane segments (S1-S4) but the S4 segment is critical since it carries positively charged residues, mainly Arg or Lys, which require an aqueous environment for its proper function. The discovery of this module in voltage-gated ion channels (VGICs), proton channels (Hv1), and voltage sensor-containing phosphatases (VSPs) has expanded our understanding of the principle of modularity in the voltage-sensing mechanism of these proteins. Here, by sequence comparison and the evaluation of the relationship between sequence composition, intrinsic flexibility, and structural analysis in 14 selected representatives of these three major protein groups, we report five interesting differences in the folding patterns of the VSD both in prokaryotes and eukaryotes. Our main findings indicate that this module is highly conserved throughout the evolutionary scale, however: (1) segments S1 to S3 in eukaryotes are significantly more hydrophobic than those present in prokaryotes; (2) the S4 segment has retained its hydrophilic character; (3) in eukaryotes the extramembranous linkers are significantly larger and more flexible in comparison with those present in prokaryotes; (4) the sensors present in the kHv1 proton channel and the ciVSP phosphatase, both of eukaryotic origin, exhibit relationships of flexibility and folding patterns very close to the typical ones found in prokaryotic voltage sensors; and (5) archaeal channels KvAP and MVP have flexibility profiles which are clearly contrasting in the S3-S4 region, which could explain their divergent activation mechanisms. Finally, to elucidate the obscure origins of this module, we show further evidence for a possible connection between voltage sensors and TolQ proteins.

Sweetening K-channels: what sugar taught us about permeation and gating

Because they enable for the modification of both viscosity and osmolarity, sugars have been used as a biophysical probe of voltage-gated K-channels for a while. Viscosity variations made it possible to measure the pore sizes in large and small conductance K-channels using techniques similar to those used in the 1980s to study the gramicidin A channel. These analyses led to the finding that the size of the internal mouth appears to be the primary cause of the conductance differences between Shaker-like channels and large conductance BK-channels. As an osmotic agent, adding sugar unilaterally causes streaming potentials that indicate H₂O/K⁺ cotransport across the BK-channel pore. Osmotic experiments on Shaker K-channels suggest that the pore gate operation and the slow inactivation displace comparable amounts of water. Functionally isolated voltage sensors allow estimation of individual osmotic work for each voltage sensing charge during voltage-activation, reporting dramatic internal and external remodeling of the Voltage Sensing Domain´s solvent exposed surfaces. Remarkably, each charge of the VSD appears to take a unique trajectory. Thus, manipulation of viscosity and osmolarity, together with 3D structures, brings in solid grounds to harmonize function and structure in membrane proteins such as K-channels and, in a wider scope, other structurally dynamic proteins.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Voltage vs. Ligand I: Structural basis of the intrinsic flexibility of S3 segment and its significance in ion channel activation

We systematically predict the internal flexibility of the S3 segment, one of the most mobile elements in the voltage-sensor domain. By analyzing the primary amino acid sequences of V-sensor containing proteins, including Hv1, TPC channels and the voltage-sensing phosphatases, we established correlations between the local flexibility and modes of activation for different members of the VGIC superfamily. Taking advantage of the structural information available, we also assessed structural aspects to understand the role played by the flexibility of S3 during the gating of the pore. We found that S3 flexibility is mainly determined by two specific regions: (1) a short NxxD motif in the N-half portion of the helix (S3a), and (2) a short sequence at the beginning of the so-called paddle motif where the segment has a kink that, in some cases, divide S3 into two distinct helices (S3a and S3b). A good correlation between the flexibility of S3 and the reported sensitivity to temperature and mechanical stretch was found. Thus, if the channel exhibits high sensitivity to heat or membrane stretch, local S3 flexibility is low. On the other hand, high flexibility of S3 is preferentially associated to channels showing poor heat and mechanical sensitivities. In contrast, we did not find any apparent correlation between S3 flexibility and voltage or ligand dependence. Overall, our results provide valuable insights into the dynamics of channel-gating and its modulation.

Impact of geometry changes in the channel pore by the gating movements on the channel's conductance

Kv 1.2 are voltage-dependent potassium channels of great biological importance. Despite the existence of many reports considering structure - function relations of the Kv 1.2 channel's quantitative domains, some details of the voltage gating remain ambiguous, or even unknown. One of the examples is the range of the S4-S6 domains motions involved in channel activation and gating. Another important question is to what extent the channel geometry influences the observable channel conductance at different voltages, and what mechanism stands behind. Does the narrowing of the pore reduce the conductance by ohmic resistance growth? The answer is surprisingly negative. But it can be explained in an alternative way by considering the fluctuations. To address these problems, we formulate geometric models that mimic the generic features of voltage sensor movement and trigger the movement of the other domains involved in gating. We carry out a complete simulation of S4-S6 domains translations and tilts. The obtained pore profiles allow to estimate the (ohmic) conductance dependency on the voltage. From a family of analysed models, we choose the one most concurring with the experimental data. The results allow to suggest the most probable scenario of S4-S6 domains movement during channel activation by membrane depolarization.

Conservation analysis of residues in the S4-S5 linker and the terminal part of the S5-P-S6 pore modulus in Kv and HCN channels: flexible determinants for the electromechanical coupling

Protein mobility is important to achieve protein function. Intrinsic flexibility associated with motion underlies this important issue and the analysis of side chain flexibility gives insights to understand it. In this work, the S5-P-S6 pore modulus (PM) of members of Kv and HCN channels was examined by a combination of sequence alignment, residue composition analysis, and intrinsic side chain flexibility. The PM sequences were organized as a database that was used to reveal and correlate the functional diversity of each analyzed family. Specifically, we focused our attention on the crucial role of the S4-S5 linker and its well-described interaction with the S6 T during the electromechanical coupling. Our analysis suggests the presence of a Gly-hinge in the middle of the S4-S5 linkers. This apparent Gly-hinge links a flexible N-terminal segment with a rigid C-terminal one, although in Kv7 channels, the latter segment is even more flexible. Instead, HCN channels exhibit a putative Thr-hinge and is rich in aromatic residues, in consequence, their linker is more rigid. Concerning S6, we confirm the presence of the two flexible kinks previously described and we provide the complete segmental flexibility profiles for the different families. Our results are discussed in terms of the relation between residue composition, conservation, and local conformational flexibility. This provides important insights to understand and differentiate the characteristic gating properties of these channels as well as their implications in cell physiology.

Detergent-mediated protein aggregation

Because detergents are commonly used to solvate membrane proteins for structural evaluation, much attention has been devoted to assessing the conformational bias imparted by detergent micelles in comparison to the native environment of the lipid bilayer. Here, we conduct six 500-ns simulations of a system with >600,000 atoms to investigate the spontaneous self assembly of dodecylphosphocholine detergent around multiple molecules of the integral membrane protein PagP. This detergent formed equatorial micelles in which acyl chains surround the protein's hydrophobic belt, confirming existing models of the detergent solvation of membrane proteins. In addition, unexpectedly, the extracellular and periplasmic apical surfaces of PagP interacted with the headgroups of detergents in other micelles 85 and 60% of the time, respectively, forming complexes that were stable for hundreds of nanoseconds. In some cases, an apical surface of one molecule of PagP interacted with an equatorial micelle surrounding another molecule of PagP. In other cases, the apical surfaces of two molecules of PagP simultaneously bound a neat detergent micelle. In these ways, detergents mediated the non-specific aggregation of folded PagP. These simulation results are consistent with dynamic light scattering experiments, which show that, at detergent concentrations ≥600 mM, PagP induces the formation of large scattering species that are likely to contain many copies of the PagP protein. Together, these simulation and experimental results point to a potentially generic mechanism of detergent-mediated protein aggregation.

Down-state model of the voltage-sensing domain of a potassium channel

Voltage-sensing domains (VSDs) of voltage-gated potassium (Kv) channels undergo a series of conformational changes upon membrane depolarization, from a down state when the channel is at rest to an up state, all of which lead to the opening of the channel pore. The crystal structures reported to date reveal the pore in an open state and the VSDs in an up state. To gain insights into the structure of the down state, we used a set of experiment-based restraints to generate a model of the down state of the KvAP VSD using molecular-dynamics simulations of the VSD in a lipid bilayer in excess water. The equilibrated VSD configuration is consistent with the biotin-avidin accessibility and internal salt-bridge data used to generate it, and with additional biotin-avidin accessibility data. In the model, both the S3b and S4 segments are displaced approximately 10 A toward the intracellular side with respect to the up-state configuration, but they do not move as a rigid body. Arginine side chains that carry the majority of the gating charge also make large excursions between the up and down states. In both states, arginines interact with water and participate in salt bridges with acidic residues and lipid phosphate groups. An important feature that emerges from the down-state model is that the N-terminal half of the S4 segment adopts a 3(10)-helical conformation, which appears to be necessary to satisfy a complex salt-bridge network.

NMR structural and dynamical investigation of the isolated voltage-sensing domain of the potassium channel KvAP: implications for voltage gating

The structure and dynamics of the isolated voltage-sensing domain (VSD) of the archaeal potassium channel KvAP was studied by high-resolution NMR. The almost complete backbone resonance assignment and partial side-chain assignment of the (2)H,(13)C,(15)N-labeled VSD were obtained for the protein domain solubilized in DPC/LDAO (2:1) mixed micelles. Secondary and tertiary structures of the VSD were characterized using secondary chemical shifts and NOE contacts. These data indicate that the spatial structure of the VSD solubilized in micelles corresponds to the structure of the domain in an open state of the channel. NOE contacts and secondary chemical shifts of amide protons indicate the presence of tightly bound water molecule as well as hydrogen bond formation involving an interhelical salt bridge (Asp62-R133) that stabilizes the overall structure of the domain. The backbone dynamics of the VSD was studied using (15)N relaxation measurements. The loop regions S1-S2 and S2-S3 were found mobile, while the S3-S4 loop (voltage-sensor paddle) was found stable at the ps-ns time scale. The moieties of S1, S2, S3, and S4 helices sharing interhelical contacts (at the level of the Asp62-R133 salt bridge) were observed in conformational exchange on the micros-ms time scale. Similar exchange-induced broadening of characteristic resonances was observed for the VSD solubilized in the membrane of lipid-protein nanodiscs composed of DMPC, DMPG, and POPC/DOPG lipids. Apparently, the observed interhelical motions represent an inherent property of the VSD of the KvAP channel and can play an important role in the voltage gating.

Voltage-dependent conformational changes of KVAP S4 segment in bacterial membrane environment

The nature and magnitude of voltage sensor conformational changes during ion channel activation are controversial. We have analyzed the topology of the K(V)AP voltage sensor domain in the absence and presence of a hyperpolarized voltage using native, right-side out membrane vesicles from E. coli. This approach does not disrupt the normal membrane environment of the channel protein and does not involve detergent solubilization. We found that voltage-dependent conformational changes are focused in the N-terminal half of the K(V)AP S4 segment, in excellent agreement with results obtained with Shaker. Homologous residues in the K(V)AP and Shaker S4 segments are transferred from the extracellular to the intracellular compartment upon hyperpolarization. Taken together with X-ray structures indicating that the K(V)AP S4 segment is outwardly displaced at 0 mV compared to S4 in a mammalian Shaker channel, our results are consistent with the idea that S4 moves further during voltage-dependent activation in K(V)AP than in Shaker.

Solution structure of the HsapBK K+ channel voltage-sensor paddle sequence

Voltage-gated potassium channels open and close in response to changes in the membrane potential. In this study, we have determined the NMR solution structure of the putative S3b-S4 voltage-sensor paddle fragment, the part that moves to mediate voltage gating, of the HsapBK potassium channel in dodecylphosphocholine (DPC) micelles. This paper presents the first structure of the S3b-S4 fragment from a BK channel. Diffusion coefficients as determined from PFG NMR experiments showed that a well-defined complex between the peptide and DPC molecules was formed. The structure reveals a helix-turn-helix motif, which is in agreement with crystal structures of other voltage-gated potassium channels, thus indicating that it is feasible to study the isolated fragment. The paddle motifs generally contain several basic residues, implicated in the gating. The critical Arg residues in this structure all reside on the surface, which is in agreement with crystal structures of K(v) channels. Similarities in the structure of the S3b-S4 fragment in BK and K(v) channels as well as important differences are seen, which may be important for explaining the details in paddle movement within a bilayer.

What makes a gate? The ins and outs of Kv-like K+ channels in plants

Gating of K(+) and other ion channels is 'hard-wired' within the channel protein. So it remains a puzzle how closely related channels in plants can show an unusually diverse range of biophysical properties. Gating of these channels lies at the heart of K(+) mineral nutrition, signalling, abiotic and biotic stress responses in plants. Thus, our knowledge of the molecular mechanics underpinning K(+) channel gating will be important for rational engineering of related traits in agricultural crops. Several key studies have added significantly to our understanding of channel gating in plants and have challenged current thinking about analogous processes found in animal K(+) channels. Such studies highlight how much of K(+) channel gating remains to be explored in plants.

Potassium channel opening: a subtle two-step

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.

Conformational changes and slow dynamics through microsecond polarized atomistic molecular simulation of an integral Kv1.2 ion channel

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 µs 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° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs 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₁₀ 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 µs). 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.

Proteins that transport ions across the cellular membrane are essential for cellular life. The proteins conducting positively charged potassium ions are key players in heart beat and nerve impulse generation because they are regulating the electrical excitability of the cell (together with proteins transporting other ions). These particular ion channels open and close in response to voltage changes across cellular membranes, but the details of this process are still not fully understood. It is, however, known that the main protein element responsible is a helical section containing several charges. Through new computer simulation methods, we have been able to run unprecedentedly long atomic simulations of an entire potassium channel embedded within a patch of membrane to help to shed new light on this gating process. Upon changing the voltage across the membrane, we observe a change in structure of this helical protein segment that appears to be an early sign of transition from the open to the closed state of the channel. This has also been previously proposed to be critical for the gating process. Understanding these structural changes on an atomic level is essential for both advancing basic science and enabling drug design targeting of voltage-regulated ion channels.

Double bilayers and transmembrane gradients: a molecular dynamics study of a highly charged peptide

The position and extent of movement of a charged peptide within a membrane bilayer provides much controversy. In our study, we have examined the nature of the highly charged helix-turn-helix motif (S3b and S4) to address how a highly charged peptide is stabilized within a bilayer in the presence of various transmembrane electrical potentials. Our double-bilayer simulation results show how the variation of the salt concentrations between the inner and outer bath establishes a transmembrane potential. Our results also show that important features of the peptide affected by changes in electrical potential are the center of mass depth, the swivel/kink degrees of conformation, and the hydrogen-bonding patterns. As the voltage gradient across the bilayer increased, the center of mass of the peptide shifted in a direction toward the outer bath. The peptide also has a higher percent helical content and the swivel/kink conformation is more rigid for nonpolarized systems where no voltage drop occurred between salt baths. Our results also provide some suggestions for how this domain may be affected by environmental changes as part of the voltage sensor in a K-channel.

Molecular dynamics simulation of Kv channel voltage sensor helix in a lipid membrane with applied electric field

In this article, we present the results of the molecular dynamics simulations of amphiphilic helix peptides of 13 amino-acid residues, placed at the lipid-water interface of dipalmitoylphosphatidylcholine bilayers. The peptides are identical with, or are derivatives of, the N-terminal segment of the S4 helix of voltage-dependent K channel KvAP, containing four voltage-sensing arginine residues (R1-R4). Upon changing the direction of the externally applied electric field, the tilt angle of the wild-type peptide changes relative to the lipid-water interface, with the N-terminus heading up with an outward electric field. These movements were not observed using an octane membrane in place of the dipalmitoylphosphatidylcholine membrane, and were markedly suppressed by 1), substituting Phe located one residue before the first arginine (R1) with a hydrophilic residue (Ser, Thr); or 2), changing the periodicity rule of Rs from at-every-third to at-every-fourth position; or 3), replacing R1 with a lysine residue (K). These and other findings suggest that the voltage-dependent movement requires deep positioning of Rs when the resting (inward) electric field is present. Later, we performed simulations of the voltage sensor domain (S1-S4) of Kv1.2 channel. In simulations with a strong electric field (0.1 V/nm or above) and positional restraints on the S1 and S2 helices, S4 movement was observed consisting of displacement along the S4 helix axis and a screwlike axial rotation. Gating-charge-carrying Rs were observed to make serial interactions with E183 in S1 and E226 in S2, in the outer water crevice. A 30-ns-backward simulation started from the open-state model gave rise to a structure similar to the recent resting-state model, with S4 moving vertically approximately 6.7 A. The energy landscape around the movement of S4 appears very ragged due to salt bridges formed between gating-charge-carrying residues and negatively charged residues of S1, S2, and S3 helices. Overall, features of S3 and S4 movements are consistent with the recent helical-screw model. Both forward and backward simulations show the presence of at least two stable intermediate structures in which R2 and R3 form salt bridges with E183 or E226, respectively. These structures are the candidates for the states postulated in previous gating kinetic models, such as the Zagotta-Hoshi-Aldrich model, to account for more than one transition step per subunit for activation.

The structure of the lipid-embedded potassium channel voltage sensor determined by double-electron-electron resonance spectroscopy

A four-pulse electron paramagnetic resonance experiment was used to measure long-range inter-subunit distances in reconstituted KvAP, a voltage-dependent potassium (Kv) channel. The measurements have allowed us to reach the following five conclusions about the native structure of the voltage sensor of KvAP. First, the S1 helix of the voltage sensor engages in a helix packing interaction with the pore domain. Second, the crystallographically observed antiparallel helix-turn-helix motif of the voltage-sensing paddle is retained in the membrane-embedded voltage sensor. Third, the paddle is oriented in such a way as to expose one face to the pore domain and the opposite face to the membrane. Fourth, the paddle and the pore domain appear to be separated by a gap that is sufficiently wide for lipids to penetrate between the two domains. Fifth, the critical voltage-sensing arginine residues on the paddle appear to be lipid exposed. These results demonstrate the importance of the membrane for the native structure of Kv channels, suggest that lipids are an integral part of their native structure, and place the voltage-sensing machinery into a complex lipid environment near the pore domain.

Dynamics of the Kv1.2 voltage-gated K+ channel in a membrane environment

All-atom molecular dynamics simulations are used to better understand the dynamic environment experienced by the Kv1.2 channel in a lipid membrane. The structure of the channel is stable during the trajectories. The pore domain keeps a well-defined conformation, whereas the voltage-sensing domains undergo important lateral fluctuations, consistent with their modular nature. A channel-like region at the center of the S1-S4 helical bundle fills rapidly with water, reminiscent of the concept of high-dielectric aqueous crevices. The first two arginines along S4 (R294 and R297) adopt an interfacial position where they interact favorably with water and the lipid headgroups. The following two arginines (R300 and R303) interact predominantly with water and E226 in S2. Despite the absence of a structurally permanent gating pore formed by protein residues and surrounding the S4 helix, as traditionally pictured, the charged residues are located in a favorable environment and are not extensively exposed to the membrane nonpolar region. Continuum electrostatic computations indicate that the transmembrane potential sensed by the charged residues in the voltage sensor varies abruptly over the outer half of the membrane in the arginine-rich region of S4; thus, the voltage gradient or membrane electric field is "focused". Interactions of basic residues with the lipid headgroups at the intracellular membrane-solution interface reduce the membrane thickness near the channel, resulting in an increased transmembrane field.

Accessibility of four arginine residues on the S4 segment of the Bacillus halodurans sodium channel

The voltage-gated Na(+) channel of Bacillus halodurans (NaChBac) is composed of six transmembrane segments (S1-S6), with a pore-forming region composed of segments S5 and S6 and a voltage-sensing domain composed of segments S1-S4. The S4 segment forms the core of the voltage sensor. We explored the accessibility of four arginine residues on the S4 segment of NaChBac, which are positioned at every third position from each other. These arginine residues on the S4 segment were replaced with cysteines using site-directed mutagenesis. Na(+) currents were recorded using the whole-cell configuration of the patch-clamp technique. We tested the effect of the sulfhydryl reagents applied from inside and outside the cellular space in the open and closed conformations. Structural models of the voltage sensor of NaChBac were constructed based on the recently crystallized KvAP and Kv1.2 K(+) channels to visualize arginine residue accessibility. Our results suggest that arginine accessibility did not change significantly between the open and closed conformations, supporting the idea of a small movement of the S4 segment during gating. Molecular modeling of the closed conformation also supported a small movement of S4, which is mainly characterized by a rotation and a tilt along the periphery of the pore. Interestingly, the second arginine residue of the S4 segment (R114) was accessible to sulfhydryl reagents from both sides of the membrane in the closed conformation and, based on our model, seemed to be at the junction of the intracellular and extracellular water crevices.

How does a voltage sensor interact with a lipid bilayer? Simulations of a potassium channel domain

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 × 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is “focused” across the center of the bilayer.

Bilayer deformation by the Kv channel voltage sensor domain revealed by self-assembly simulations

Coarse-grained molecular dynamics simulations are used to explore the interaction with a phospholipid bilayer of the voltage sensor (VS) domain and the S4 helix from the archaebacterial voltage-gated potassium (Kv) channel KvAP. Multiple 2-mus self-assembly simulations reveal that the isolated S4 helix may adopt either interfacial or transmembrane (TM) locations with approximately equal probability. In the TM state, the insertion of the voltage-sensing region of S4 is facilitated via local bilayer deformation that, combined with side chain "snorkeling," enables its Arg side chains to interact with lipid headgroups and water. Multiple 0.2-mus self-assembly simulations of the VS domain are also performed, along with simulations of MscL and KcsA, to permit comparison with more "canonical" integral membrane protein structures. All three stably adopt a TM orientation within a bilayer. For MscL and KcsA, there is no significant bilayer deformation. In contrast, for the VS, there is considerable local deformation, which is again primarily due to the lipid-exposed S4. It is shown that for both the VS and isolated S4 helix, the positively charged side chains of S4 are accommodated within the membrane through a combination of stabilizing interactions with lipid glycerol and headgroup regions, water, and anionic side chains. Our results support the possibility that bilayer deformation around key gating charge residues in Kv channels may result in "focusing" of the electrostatic field, and indicate that, when considering competing models of voltage-sensing, it is essential to consider the dynamics and structure of not only the protein but also of the local lipid environment.

The interaction of spider gating modifier peptides with voltage-gated potassium channels

Gating modifier peptides bind to ion channels and alter the gating process of these molecules. One of the most extensively studied peptides, Hanatoxin (HaTx), isolated from a Chilean tarantula, has been used to characterize the blocking properties of the voltage-gated potassium channel Kv2.1. These studies have provided some insight into the gating mechanism in Kv channels. In this review we will discuss the interaction of HaTx and related spider peptides with Kv channels illustrating the properties of the binding surface of these peptides, their membrane partitioning characteristics, and will provide a working hypothesis for how the peptides inhibit gating of Kv channels. Advanced simulation results support the concept of mutual conformational changes upon peptide binding to the S3b region of the channel which will restrict movement of S4 and compromise coupling of the gating machinery to opening of the pore.

A voltage-sensor water pore

Voltage-sensor (VS) domains cause the pore of voltage-gated ion channels to open and close in response to changes in transmembrane potential. Recent experimental studies suggest that VS domains are independent structural units. This independence is revealed dramatically by a voltage-dependent proton-selective channel (Hv), which has a sequence homologous to the VS domains of voltage-gated potassium channels (Kv). Here we show by means of molecular dynamics simulations that the isolated open-state VS domain of the KvAP channel in a lipid membrane has a configuration consistent with a water channel, which we propose as a common feature underlying the conductance of protons, and perhaps other cations, through VS domains.

Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation

Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a "ring" of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7-0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.