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

The differential impacts of equivalent gating-charge mutations in voltage-gated sodium channels

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling, and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function, loss-of-function effects, or both. However, the mechanism behind this functional diversity of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating property impacts in electrophysiology measurements. We conducted molecular dynamics (MD) simulations on wild-type channels and six mutants to elucidate the structural basis of their differential impacts. Our 120-µs MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis revealed how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study advances our understanding of structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

A Deep Learning Approach to Uncover Voltage-Gated Ion Channels' Intermediate States

Owing to recent advancements in cryo-electron microscopy, voltage-gated ion channels have gained a greater comprehension of their structural characteristics. However, a significant enigma remains unsolved for a large majority of these channels: their gating mechanism. This mechanism, which encompasses the conformational changes between open and closed states, is pivotal to their proper functioning. Beyond the binary states of open and closed, an ensemble of intermediate states defines the transition path in-between. Due to the lack of experimental data, one might resort to molecular dynamics simulations as an alternative to decipher these states and the transitions between them. However, the high-energy barriers and the colossal time scales involved hinder access to the latter. We present here an application of deep learning as a reliable pipeline for a comprehensive exploration of voltage-gated ion channel conformational rearrangements during gating. We showcase the pipeline performance specifically on the Kv1.2 voltage sensor domain and confront the results with existing data. We demonstrate how our physics-based deep learning approach contributes to the theoretical understanding of these channels and how it might provide further insights into the exploration of channelopathies.

ELUCIDATING THE DIFFERENTIAL IMPACTS OF EQUIVALENT GATING-CHARGE MUTATIONS IN VOLTAGE-GATED SODIUM CHANNELS

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function (GoF), loss-of-function (LoF) effects, or both. However, the mechanism behind this functional divergence of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino-acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating-property impacts in electrophysiology measurements. We conducted 120 μs molecular dynamics (MD) simulations on wild-type and six mutants to elucidate the structural basis of their differential impacts. Our μs-scale MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis elucidated how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study represents a significant leap forward in understanding structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Anionic lipids unlock the gates of select ion channels in the pacemaker family

Lipids play important roles in regulating membrane protein function, but the molecular mechanisms used are elusive. Here we investigated how anionic lipids modulate SthK, a bacterial pacemaker channel homolog, and HCN2, whose activity contributes to pacemaking in the heart and brain. Using SthK allowed the reconstitution of purified channels in controlled lipid compositions for functional and structural assays that are not available for the eukaryotic channels. We identified anionic lipids bound tightly to SthK and their exact binding locations and determined that they potentiate channel activity. Cryo-EM structures in the most potentiating lipids revealed an open state and identified a nonannular lipid bound with its headgroup near an intersubunit salt bridge that clamps the intracellular channel gate shut. Breaking this conserved salt bridge abolished lipid modulation in SthK and eukaryotic HCN2 channels, indicating that anionic membrane lipids facilitate channel opening by destabilizing these interactions. Our findings underline the importance of state-dependent protein-lipid interactions.

Interaction fingerprint of transmembrane segments in voltage sensor domains

Voltage sensor domain (VSD) in channel and non-channel membrane proteins shares a common function in the detection of changes in the transmembrane electric potential. The VSD is made of four helical transmembrane segments (S1-S4) that form a structurally conserved scaffold through inter-transmembrane residue-residue interactions. Details about these interactions are yet to be fully understood in the context of the unique structural and physical characteristics of the voltage sensor unit. In this study, molecular dynamics simulations were carried out to investigate transmembrane helix-helix interactions via residue-based nonbonding energies using the activated and resting state conformations of VSD from Hv1, CiVSP, KvAP and NavAb. Inter-transmembrane interaction energies within the VSD were determined. Analysis of electrostatic and van der Waals components revealed the strengths and weaknesses of the interactions between each pair of transmembrane segments. In all cases the S4 helix had the highest electrostatic contribution to favor the key role as the voltage sensitive segment. Electrostatic interactions for the S1-S2 pair as well as the S1-S3 pair were relatively weak. Van der Waal interaction energies between adjacent segments were on average greater than that between diagonally opposite segments. Salt bridge interactions between S4-arginines and the negatively charged residues in other segments appear to contribute more to stabilizing the energy than the van der Waals interactions between nonpolar residues. The overall behavior of residue-residue contacts is similar among the transmembrane domains, reflecting the common inter- transmembrane interaction pattern in the VSD. In addition, analysis of the residue positions suggested that subtle differences in the orientation of the salt-bridges can be attributed to the difference in the inter-transmembrane interaction strengths inside the VSDs.

Regulation of a pentameric ligand-gated ion channel by a semiconserved cationic lipid-binding site

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of noncanonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid-binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABAA receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABAA receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Tracking Membrane Protein Dynamics in Real Time

Abstract

Membrane proteins govern critical cellular processes and are central to human health and associated disease. Understanding of membrane protein function is obscured by the vast ranges of structural dynamics—both in the spatial and time regime—displayed in the protein and surrounding membrane. The membrane lipids have emerged as allosteric modulators of membrane protein function, which further adds to the complexity. In this review, we discuss several examples of membrane dependency. A particular focus is on how molecular dynamics (MD) simulation have aided to map membrane protein dynamics and how enhanced sampling methods can enable observing the otherwise inaccessible biological time scale. Also, time-resolved X-ray scattering in solution is highlighted as a powerful tool to track membrane protein dynamics, in particular when combined with MD simulation to identify transient intermediate states. Finally, we discuss future directions of how to further develop this promising approach to determine structural dynamics of both the protein and the surrounding lipids.

Graphic Abstract

Tracking the motion of the KV1.2 voltage sensor reveals the molecular perturbations caused by a de novo mutation in a case of epilepsy

KEY POINTS

KV1.2 channels, encoded by the KCNA2 gene, regulate neuronal excitability by conducting K+ upon depolarization. A new KCNA2 missense variant was discovered in a patient with epilepsy, causing amino acid substitution F302L at helix S4, in the KV1.2 voltage-sensing domain. Immunocytochemistry and flow cytometry showed that F302L does not impair KCNA2 subunit surface trafficking. Molecular dynamics simulations indicated that F302L alters the exposure of S4 residues to membrane lipids. Voltage clamp fluorometry revealed that the voltage-sensing domain of KV1.2-F302L channels is more sensitive to depolarization. Accordingly, KV1.2-F302L channels opened faster and at more negative potentials; however, they also exhibited enhanced inactivation: that is, F302L causes both gain- and loss-of-function effects. Coexpression of KCNA2-WT and -F302L did not fully rescue these effects. The proband's symptoms are more characteristic of patients with loss of KCNA2 function. Enhanced KV1.2 inactivation could lead to increased synaptic release in excitatory neurons, steering neuronal circuits towards epilepsy.

ABSTRACT

An exome-based diagnostic panel in an infant with epilepsy revealed a previously unreported de novo missense variant in KCNA2, which encodes voltage-gated K+ channel KV1.2. This variant causes substitution F302L, in helix S4 of the KV1.2 voltage-sensing domain (VSD). F302L does not affect KCNA2 subunit membrane trafficking. However, it does alter channel functional properties, accelerating channel opening at more hyperpolarized membrane potentials, indicating gain of function. F302L also caused loss of KV1.2 function via accelerated inactivation onset, decelerated recovery and shifted inactivation voltage dependence to more negative potentials. These effects, which are not fully rescued by coexpression of wild-type and mutant KCNA2 subunits, probably result from the enhancement of VSD function, as demonstrated by optically tracking VSD depolarization-evoked conformational rearrangements. In turn, molecular dynamics simulations suggest altered VSD exposure to membrane lipids. Compared to other encephalopathy patients with KCNA2 mutations, the proband exhibits mild neurological impairment, more characteristic of patients with KCNA2 loss of function. Based on this information, we propose a mechanism of epileptogenesis based on enhanced KV1.2 inactivation leading to increased synaptic release preferentially in excitatory neurons, and hence the perturbation of the excitatory/inhibitory balance of neuronal circuits.

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.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP₂ plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP₂ and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABA_A receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP₂ also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

HIV-1 Env gp41 Transmembrane Domain Dynamics Are Modulated by Lipid, Water, and Ion Interactions

The gp41 transmembrane domain (TMD) of the envelope glycoprotein of the human immunodeficiency virus modulates the conformation of the viral envelope spike, the only druggable target on the surface of the virion. Targeting the envelope glycoprotein with small-molecule and antibody therapies requires an understanding of gp41 TMD dynamics, which is often challenging given the difficulties in describing native membrane properties. Here, atomistic molecular dynamics simulations of a trimeric, prefusion gp41 TMD in a model, asymmetric viral membrane that mimics the native viral envelope were performed. Water and chloride ions were observed to permeate the membrane and interact with the highly conserved arginine bundle, (R696)₃, at the center of the membrane and influenced TMD stability by creating a network of hydrogen bonds and electrostatic interactions. We propose that this (R696)₃ - water - anion network plays an important role in viral fusion with the host cell by modulating protein conformational changes within the membrane. Additionally, R683 and R707 at the exofacial and cytofacial membrane-water interfaces, respectively, are anchored in the lipid headgroup region and serve as a junction point for stabilization of the termini. The membrane thins as a result of the tilting of the gp41 trimer with nearby lipids increasing in volume, leading to an entropic driving force for TMD conformational change. These results provide additional detail and perspective on the influence of certain lipid types on TMD dynamics and a rationale for targeting key residues of the TMD for therapeutic design. These insights into the molecular details of TMD membrane anchoring will build toward a greater understanding of the dynamics that lead to viral fusion with the host cell.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Towards understanding the molecular basis of ion channel modulation by lipids: Mechanistic models and current paradigms

Research on ion channel modulation has become a hot topic because of the key roles these membrane proteins play in both prokaryotic and eukaryotic organisms. In this respect, lipid modulation adds to the overall modulatory mechanisms as a potential via to find new pharmacological targets for drug design based on interfering with lipid/channel interactions. However, our knowledge in this field is scarce and often circumscribed to the sites where lipids bind and/or its final functional consequences. To fully understand this process it is necessary to improve our knowledge on its molecular basis, from the binding sites to the signalling pathways that derive in structural and functional effects on the ion channel. In this review, we have compiled information about such mechanisms and established a classification into four different modes of action. Afterwards, we have revised in more detail the lipid modulation of Cys-loop receptors and of the potassium channel KcsA, which were chosen as model channels modulated by specific lipids. This article is part of a Special Issue entitled: Membrane Lipid Therapy: Drugs Targeting Biomembranes edited by Pablo V. Escribá.

Proton currents constrain structural models of voltage sensor activation

The Hv1 proton channel is evidently unique among voltage sensor domain proteins in mediating an intrinsic 'aqueous' H⁺ conductance (G_AQ). Mutation of a highly conserved 'gating charge' residue in the S4 helix (R1H) confers a resting-state H⁺ 'shuttle' conductance (G_SH) in VGCs and Ci VSP, and we now report that R1H is sufficient to reconstitute G_SH in Hv1 without abrogating G_AQ. Second-site mutations in S3 (D185A/H) and S4 (N4R) experimentally separate G_SH and G_AQ gating, which report thermodynamically distinct initial and final steps, respectively, in the Hv1 activation pathway. The effects of Hv1 mutations on G_SH and G_AQ are used to constrain the positions of key side chains in resting- and activated-state VS model structures, providing new insights into the structural basis of VS activation and H⁺ transfer mechanisms in Hv1.

Capsaicin interaction with TRPV1 channels in a lipid bilayer: molecular dynamics simulation

Transient receptor potential vanilloid subtype 1 (TRPV1) is a heat-sensitive ion channel also involved in pain sensation, and is the receptor for capsaicin, the active ingredient of hot chili peppers. The recent structures of TRPV1 revealed putative ligand density within the S1 to S4 voltage-sensor-like domain of the protein. However, questions remain regarding the dynamic role of the lipid bilayer in ligand binding to TRPV1. Molecular dynamics simulations were used to explore behavior of capsaicin in a 1-palmitoyl-2-oleoyl phosphatidylcholine bilayer and with the target S1–S4 transmembrane helices of TRPV1. Equilibrium simulations reveal a preferred interfacial localization for capsaicin. We also observed a capsaicin molecule flipping from the extracellular to the intracellular leaflet, and subsequently able to access the intracellular TRPV1 binding site. Calculation of the potential of mean force (i.e., free energy profile) of capsaicin along the bilayer normal confirms that it prefers an interfacial localization. The free energy profile indicates that there is a nontrivial but surmountable barrier to the flipping of capsaicin between opposing leaflets of the bilayer. Molecular dynamics of the S1–S4 transmembrane helices of the TRPV1 in a lipid bilayer confirm that Y511, known to be crucial to capsaicin binding, has a distribution along the bilayer normal similar to that of the aromatic group of capsaicin. Simulations were conducted of the TRPV1 S1–S4 transmembrane helices in the presence of capsaicin placed in the aqueous phase, in the lipid, or docked to the protein. No stable interaction between ligand and protein was seen for simulations initiated with capsaicin in the bilayer. However, interactions were seen between TRPV1 and capsaicin starting from the cytosolic aqueous phase, and capsaicin remained stable in the majority of simulations from the docked pose. We discuss the significance of capsaicin flipping from the extracellular to the intracellular leaflet and mechanisms of binding site access by capsaicin.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Voltage Sensing in Membranes: From Macroscopic Currents to Molecular Motions

Voltage-sensing domains (VSDs) are integral membrane protein units that sense changes in membrane electric potential, and through the resulting conformational changes, regulate a specific function. VSDs confer voltage-sensitivity to a large superfamily of membrane proteins that includes voltage-gated Na[Formula: see text], K[Formula: see text], Ca[Formula: see text] ,and H[Formula: see text] selective channels, hyperpolarization-activated cyclic nucleotide-gated channels, and voltage-sensing phosphatases. VSDs consist of four transmembrane segments (termed S1 through S4). Their most salient structural feature is the highly conserved positions for charged residues in their sequences. S4 exhibits at least three conserved triplet repeats composed of one basic residue (mostly arginine) followed by two hydrophobic residues. These S4 basic side chains participate in a state-dependent internal salt-bridge network with at least four acidic residues in S1-S3. The signature of voltage-dependent activation in electrophysiology experiments is a transient current (termed gating or sensing current) upon a change in applied membrane potential as the basic side chains in S4 move across the membrane electric field. Thus, the unique structural features of the VSD architecture allow for competing requirements: maintaining a series of stable transmembrane conformations, while allowing charge motion, as briefly reviewed here.

Evolutionary imprint of activation: the design principles of VSDs

Voltage-sensor domains (VSDs) are modular biomolecular machines that transduce electrical signals in cells through a highly conserved activation mechanism. Here, we investigate sequence-function relationships in VSDs with approaches from information theory and probabilistic modeling. Specifically, we collect over 6,600 unique VSD sequences from diverse, long-diverged phylogenetic lineages and relate the statistical properties of this ensemble to functional constraints imposed by evolution. The VSD is a helical bundle with helices labeled S1-S4. Surrounding conserved VSD residues such as the countercharges and the S2 phenylalanine, we discover sparse networks of coevolving residues. Additional networks are found lining the VSD lumen, tuning the local hydrophilicity. Notably, state-dependent contacts and the absence of coevolution between S4 and the rest of the bundle are imprints of the activation mechanism on the VSD sequence ensemble. These design principles rationalize existing experimental results and generate testable hypotheses.

Mechanisms of activation of voltage-gated potassium channels

Voltage-gated potassium ion channels (Kv) play an important role in a variety of cellular processes, including the functioning of excitable cells, regulation of apoptosis, cell growth and differentiation, the release of neurotransmitters and hormones, maintenance of cardiac activity, etc. Failure in the functioning of Kv channels leads to severe genetic disorders and the development of tumors, including malignant ones. Understanding the mechanisms underlying Kv channels functioning is a key factor in determining the cause of the diseases associated with mutations in the channels, and in the search for new drugs. The mechanism of activation of the channels is a topic of ongoing debate, and a consensus on the issue has not yet been reached. This review discusses the key stages in studying the mechanisms of functioning of Kv channels and describes the basic models of their activation known to date.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Voltage-gated ion channel modulation by lipids: insights from molecular dynamics simulations

Cells commonly use lipids to modulate the function of ion channels. The lipid content influences the amplitude of the ionic current and changes the probability of voltage-gated ion channels being in the active or in the resting states. Experimental findings inferred from a variety of techniques and molecular dynamics studies have revealed a direct interaction between the lipid headgroups and the ion channel residues, suggesting an influence on the ion channel function. On the other hand the alteration of the lipids may in principle modify the overall electrostatic environment of the channel, and hence the transmembrane potential, leading to an indirect modulation, i.e. a global effect. Here we have investigated the structural and dynamical properties of the voltage-gated potassium channel Kv1.2 embedded in bilayers with modified upper or lower leaflet compositions corresponding to realistic biological scenarios: the first relates to the effects of sphingomyelinase, an enzyme that modifies the composition of lipids of the outer membrane leaflets, and the second to the effect of the presence of a small fraction of PIP2, a highly negatively charged lipid known to modulate voltage-gated channel function. Our molecular dynamics simulations do not enable to exclude the global effect mechanism in the former case. For the latter, however, it is shown that local interactions between the ion channel and the lipid headgroups are key-elements of the modulation.

Omega currents in voltage-gated ion channels: what can we learn from uncovering the voltage-sensing mechanism using MD simulations?

Ion channels conduct charged species through otherwise impermeable biological membranes. Their activity supports a number of physiological processes, and genetic mutations can disrupt their function dramatically. Among these channels, voltage gated cation channels (VGCCs) are ubiquitous transmembrane proteins involved in electrical signaling. In addition to their selectivity for ions, their function requires membrane-polarization-dependent gating. Triggered by changes in the transmembrane voltage, the activation and deactivation of VGCCs proceed through a sensing mechanism that prompts motion of conserved positively charged (basic) residues within the S4 helix of a four-helix bundle, the voltage sensor domain (VSD). Decades of experimental investigations, using electrophysiology, molecular biology, pharmacology, and spectroscopy, have revealed details about the function of VGCCs. However, in 2005, the resolution of the crystal structure of the activated state of one member of the mammalian voltage gated potassium (Kv) channels family (the Kv1.2) enabled researchers to make significant progress in understanding the structure-function relationship in these proteins on a molecular level. In this Account, we review the use of a complementary technique, molecular dynamics (MD) simulations, that has offered new insights on this timely issue. Starting from the "open-activated state" crystal structure, we have carried out large-scale all atom MD simulations of the Kv1.2 channel embedded in its lipidic environment and submitted to a hyperpolarizing (negative) transmembrane potential. We then used steered MD simulations to complete the full transition to the resting-closed state. Using these procedures, we have followed the operation of the VSDs and uncovered three intermediate states between their activated and deactivated conformations. Each conformational state is characterized by its network of salt bridges and by the occupation of the gating charge transfer center by a specific S4 basic residue. Overall, the global deactivation mechanism that we have uncovered agrees with proposed kinetic models and recent experimental results that point towards the presence of several intermediate states. The understanding of these conformations has allowed us to examine how mutations of the S4 basic residues analogous to those involved in genetic diseases affect the function of VGCCs. In agreement with electrophysiology experiments, mutations perturb the VSD structure and trigger the appearance of state-dependent "leak" currents. The simulation results unveil the key elementary molecular processes involved in these so-called "omega" currents. We generalize these observations to other members of the VGCC family, indicating which type of residues may generate such currents and which conditions might cause leaks that prevent proper function of the channel. Today, the understanding of the intermediate state conformations enables researchers to confidently tackle other key questions such as the mode of action of toxins or modulation of channel function by lipids.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1

The topological similarity of voltage-gated proton channels (H_V1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H_V1s have a similar structure. We present the construction and validation of a homology model of the human H_V1 (hH_V1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH_V1 by exploiting the homology of hH_V1 with VSDs of K⁺ and Na⁺ channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH_V1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp¹¹². Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH_V1. Each Arg residue in the S4 helix of hH_V1 was replaced by His to test accessibility using Zn²⁺ as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg²¹¹, the third Arg residue in the S4 helix in hH_V1, remains accessible to the internal solution and is located near the charge transfer center, Phe¹⁵⁰.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Challenges in the Development of Functional Assays of Membrane Proteins

Lipid bilayers are natural barriers of biological cells and cellular compartments. Membrane proteins integrated in biological membranes enable vital cell functions such as signal transduction and the transport of ions or small molecules. In order to determine the activity of a protein of interest at defined conditions, the membrane protein has to be integrated into artificial lipid bilayers immobilized on a surface. For the fabrication of such biosensors expertise is required in material science, surface and analytical chemistry, molecular biology and biotechnology. Specifically, techniques are needed for structuring surfaces in the micro- and nanometer scale, chemical modification and analysis, lipid bilayer formation, protein expression, purification and solubilization, and most importantly, protein integration into engineered lipid bilayers. Electrochemical and optical methods are suitable to detect membrane activity-related signals. The importance of structural knowledge to understand membrane protein function is obvious. Presently only a few structures of membrane proteins are solved at atomic resolution. Functional assays together with known structures of individual membrane proteins will contribute to a better understanding of vital biological processes occurring at biological membranes. Such assays will be utilized in the discovery of drugs, since membrane proteins are major drug targets.

Structural interactions between lipids, water and S1-S4 voltage-sensing domains

Membrane proteins serve crucial signaling and transport functions, yet relatively little is known about their structures in membrane environments or how lipids interact with these proteins. For voltage-activated ion channels, X-ray structures suggest that the mobile voltage-sensing S4 helix would be exposed to the membrane, and functional studies reveal that lipid modification can profoundly alter channel activity. Here, we use solid-state NMR to investigate structural interactions of lipids and water with S1-S4 voltage-sensing domains and to explore whether lipids influence the structure of the protein. Our results demonstrate that S1-S4 domains exhibit extensive interactions with lipids and that these domains are heavily hydrated when embedded in a membrane. We also find evidence for preferential interactions of anionic lipids with S1-S4 domains and that these interactions have lifetimes on the timescale of ≤ 10(-3)s. Arg residues within S1-S4 domains are well hydrated and are positioned in close proximity to lipids, exhibiting local interactions with both lipid headgroups and acyl chains. Comparative studies with a positively charged lipid lacking a phosphodiester group reveal that this lipid modification has only modest effects on the structure and hydration of S1-S4 domains. Taken together, our results demonstrate that Arg residues in S1-S4 voltage-sensing domains reside in close proximity to the hydrophobic interior of the membrane yet are well hydrated, a requirement for carrying charge and driving protein motions in response to changes in membrane voltage.

Molecular dynamics simulations of voltage-gated cation channels: insights on voltage-sensor domain function and modulation

Since their discovery in the 1950s, the structure and function of voltage-gated cation channels (VGCC) has been largely understood thanks to results stemming from electrophysiology, pharmacology, spectroscopy, and structural biology. Over the past decade, computational methods such as molecular dynamics (MD) simulations have also contributed, providing molecular level information that can be tested against experimental results, thereby allowing the validation of the models and protocols. Importantly, MD can shed light on elements of VGCC function that cannot be easily accessed through “classical” experiments. Here, we review the results of recent MD simulations addressing key questions that pertain to the function and modulation of the VGCC’s voltage-sensor domain (VSD) highlighting: (1) the movement of the S4-helix basic residues during channel activation, articulating how the electrical driving force acts upon them; (2) the nature of the VSD intermediate states on transitioning between open and closed states of the VGCC; and (3) the molecular level effects on the VSD arising from mutations of specific S4 positively charged residues involved in certain genetic diseases.

Membrane-perturbing properties of two Arg-rich paddle domains from voltage-gated sensors in the KvAP and HsapBK K(+) channels

Voltage-gated K(+) channels are gated by displacement of basic residues located in the S4 helix that together with a part of the S3 helix, S3b, forms a "paddle" domain, whose position is altered by changes in the membrane potential modulating the open probability of the channel. Here, interactions between two paddle domains, KvAPp from the K(v) channel from Aeropyrum pernix and HsapBKp from the BK channel from Homo sapiens, and membrane models have been studied by spectroscopy. We show that both paddle domains induce calcein leakage in large unilamellar vesicles, and we suggest that this leakage represents a general thinning of the bilayer, making movement of the whole paddle domain plausible. The fact that HsapBKp induces more leakage than KvAPp may be explained by the presence of a Trp residue in HsapBKp. Trp residues generally promote localization to the hydrophilic-hydrophobic interface and disturb tight packing. In magnetically aligned bicelles, KvAPp increases the level of order along the whole acyl chain, while HsapBKp affects the morphology, also indicating that KvAPp adapts more to the lipid environment. Nuclear magnetic resonance (NMR) relaxation measurements for HsapBKp show that overall the sequence has anisotropic motions. The S4 helix is well-structured with restricted local motion, while the turn between S4 and S3b is more flexible and undergoes slow local motion. Our results indicate that the calcein leakage is related to the flexibility in this turn region. A possibility by which HsapBKp can undergo structural transitions is also shown by relaxation NMR, which may be important for the gating mechanism.

The influence of lipids on voltage-gated ion channels

Voltage-gated ion channels are responsible for transmitting electrochemical signals in both excitable and non-excitable cells. Structural studies of voltage-gated potassium and sodium channels by X-ray crystallography have revealed atomic details on their voltage-sensor domains (VSDs) and pore domains, and were put in context of disparate mechanistic views on the voltage-driven conformational changes in these proteins. Functional investigation of voltage-gated channels in membranes, however, showcased a mechanism of lipid-dependent gating for voltage-gated channels, suggesting that the lipids play an indispensible and critical role in the proper gating of many of these channels. Structure determination of membrane-embedded voltage-gated ion channels appears to be the next frontier in fully addressing the mechanism by which the VSDs control channel opening. Currently electron crystallography is the only structural biology method in which a membrane protein of interest is crystallized within a complete lipid-bilayer mimicking the native environment of a biological membrane. At a sufficiently high resolution, an electron crystallographic structure could reveal lipids, the channel and their mutual interactions at the atomic level. Electron crystallography is therefore a promising avenue toward understanding how lipids modulate channel activation through close association with the VSDs.

In search of a consensus model of the resting state of a voltage-sensing domain

Voltage-sensing domains (VSDs) undergo conformational changes in response to the membrane potential and are the critical structural modules responsible for the activation of voltage-gated channels. Structural information about the key conformational states underlying voltage activation is currently incomplete. Through the use of experimentally determined residue-residue interactions as structural constraints, we determine and refine a model of the Kv channel VSD in the resting conformation. The resulting structural model is in broad agreement with results that originate from various labs using different techniques, indicating the emergence of a consensus for the structural basis of voltage sensing.

Functional extension of amino acid triads from the fourth transmembrane segment (S4) into its external linker in Shaker K(+) channels

The highly conserved fourth transmembrane segment (S4) is the primary voltage sensor of the voltage-dependent channel and would move outward upon membrane depolarization. S4 comprises repetitive amino acid triads, each containing one basic (presumably charged and voltage-sensing) followed by two hydrophobic residues. We showed that the triad organization is functionally extended into the S3-4 linker right external to S4 in Shaker K(+) channels. The arginine (and lysine) substitutes for the third and the sixth residues (Ala-359 and Met-356, respectively) external to the outmost basic residue (Arg-362) in S4 dramatically and additively stabilize S4 in the resting conformation. Also, Leu-361 and Leu-358 play a very similar role in stabilization of S4 in the resting position, presumably by their hydrophobic side chains. Moreover, the double mutation A359R/E283A leads to a partially extruded position of S4 and consequently prominent closed-state inactivation, suggesting that Glu-283 in S2 may coordinate with the arginines in the extruded S4 upon depolarization. We conclude that the triad organization extends into the S3-4 linker for about six amino acids in terms of their microenvironment. These approximately six residues should retain the same helical structure as S4, and their microenvironment serves as part of the "gating canal" accommodating the extruding S4. Upon depolarization, S4 most likely moves initially as a sliding helix and follows the path that is set by the approximately six residues in the S3-4 linker in the resting state, whereas further S4 translocation could be more like, for example, a paddle, without orderly coordination from the contiguous surroundings.

Water wires in atomistic models of the Hv1 proton channel

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.

Structural dynamics of the S4 voltage-sensor helix in lipid bilayers lacking phosphate groups

Voltage-dependent K(+) (Kv) channels require lipid phosphates for functioning. The S4 helix, which carries the gating charges in the voltage-sensing domain (VSD), inserts into membranes while being stabilized by a protein-lipid interface in which lipid phosphates play an essential role. To examine the physical basis of the protein-lipid interface in the absence of lipid phosphates, we performed molecular dynamics (MD) simulations of a KvAP S4 variant (S4mut) in bilayers with and without lipid phosphates. We find that, in dioleoyltrimethylammoniumpropane (DOTAP) bilayers lacking lipid phosphates, the gating charges are solvated by anionic counterions and, hence, lack the bilayer support provided by phosphate-containing palmitoyloleoylglycerophosphocholine (POPC) bilayers. The result is a water-permeable bilayer with significantly smaller deformations around the peptide. Together, these results provide an explanation for the nonfunctionality of VSDs in terms of a destabilizing protein-lipid interface.

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.

Interaction of diverse voltage sensor homologs with lipid bilayers revealed by self-assembly simulations

Voltage sensors (VS) domains couple the activation of ion channels/enzymes to changes in membrane voltage. We used molecular dynamics simulations to examine interactions with lipids of several VS homologs. VSs in intact channels in the activated state are exposed to phospholipids, leading to a characteristic local distortion of the lipid bilayer which decreases its thickness by ∼10 Å. This effect is mediated by a conserved hydrophilic stretch in the S4-S5 segment linking the VS and the pore domains, and may favor gating charges crossing the membrane. In cationic lipid bilayers lacking phosphate groups, VSs form fewer contacts with lipid headgroups. The S3-S4 paddle motifs show persistent interactions of individual lipid molecules, influenced by the hairpin loop. In conclusion, our results suggest common interactions with phospholipids for various VS homologs, providing insights into the molecular basis of their stabilization in the membrane and how they are altered by lipid modification.

3₁₀-helix conformation facilitates the transition of a voltage sensor S4 segment toward the down state

The activation of voltage-gated ion channels is controlled by the S4 helix, with arginines every third residue. The x-ray structures are believed to reflect an open-inactivated state, and models propose combinations of translation, rotation, and tilt to reach the resting state. Recently, experiments and simulations have independently observed occurrence of 3(10)-helix in S4. This suggests S4 might make a transition from α- to 3(10)-helix in the gating process. Here, we show 3(10)-helix structure between Q1 and R3 in the S4 segment of a voltage sensor appears to facilitate the early stage of the motion toward a down state. We use multiple microsecond-steered molecular simulations to calculate the work required for translating S4 both as α-helix and transformed to 3(10)-helix. The barrier appears to be caused by salt-bridge reformation simultaneous to R4 passing the F233 hydrophobic lock, and it is almost a factor-two lower with 3(10)-helix. The latter facilitates translation because R2/R3 line up to face E183/E226, which reduces the requirement to rotate S4. This is also reflected in a lower root mean-square deviation distortion of the rest of the voltage sensor. This supports the 3(10) hypothesis, and could explain some of the differences between the open-inactivated- versus activated-states.

Calculation of the gating charge for the Kv1.2 voltage-activated potassium channel

The atomic models of the Kv1.2 potassium channel in the active and resting state, originally presented elsewhere, are here refined using molecular dynamics simulations in an explicit membrane-solvent environment. With a minor adjustment of the orientation of the first arginine along the S4 segment, the total gating charge of the channel determined from >0.5 mus of molecular dynamics simulation is approximately 12-12.7 e, in good accord with experimental estimates for the Shaker potassium channel, indicating that the final models offer a realistic depiction of voltage-gating. In the resting state of Kv1.2, the S4 segment in the voltage-sensing domain (VSD) spontaneously converts into a 3(10) helix over a stretch of 10 residues. The 3(10) helical conformation orients the gating arginines on S4 toward a water-filled crevice within the VSD and allows salt-bridge interactions with negatively charged residues along S2 and S3. Free energy calculations of the fractional transmembrane potential, acting upon key charged residues of the VSD, reveals that the applied field varies rapidly over a narrow region of 10-15 A corresponding to the outer leaflet of the bilayer. The focused field allows the transfer of a large gating charge without translocation of S4 across the membrane.

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.

High-resolution orientation and depth of insertion of the voltage-sensing S4 helix of a potassium channel in lipid bilayers

The opening and closing of voltage-gated potassium (Kv) channels are controlled by several conserved Arg residues in the S4 helix of the voltage-sensing domain. The interaction of these positively charged Arg residues with the lipid membrane has been of intense interest for understanding how membrane proteins fold to allow charged residues to insert into lipid bilayers against free-energy barriers. Using solid-state NMR, we have now determined the orientation and insertion depth of the S4 peptide of the KvAP channel in lipid bilayers. Two-dimensional (15)N correlation experiments of macroscopically oriented S4 peptide in phospholipid bilayers revealed a tilt angle of 40 degrees and two possible rotation angles differing by 180 degrees around the helix axis. Remarkably, the tilt angle and one of the two rotation angles are identical to those of the S4 helix in the intact voltage-sensing domain, suggesting that interactions between the S4 segment and other helices of the voltage-sensing domain are not essential for the membrane topology of the S4 helix. (13)C-(31)P distances between the S4 backbone and the lipid (31)P indicate a approximately 9 A local thinning and 2 A average thinning of the DMPC (1,2-dimyristoyl-sn-glycero-3-phosphochloline)/DMPG (1,2-dimyristoyl-sn-glycero-3-phosphatidylglycerol) bilayer, consistent with neutron diffraction data. Moreover, a short distance of 4.6 A from the guanidinium C(zeta) of the second Arg to (31)P indicates the existence of guanidinium phosphate hydrogen bonding and salt bridges. These data suggest that the structure of the Kv gating helix is mainly determined by protein-lipid interactions instead of interhelical protein-protein interactions, and the S4 amino acid sequence encodes sufficient information for the membrane topology of this crucial gating helix.