Search

Small-molecule modulators of diverse voltage-gated K+ (Kv) channels may help treat severe neurological disorders. However, the development of selective modulators requires an understanding of their mechanism-of-action (MoA). We applied an orthogonal approach to elucidate the MoA of an imidazolidinedione derivative (AUT5), which is a highly specific positive allosteric modulator (PAM) of Kv3.1 and Kv3.2 channels. AUT5 modulation involves positive cooperativity and preferential stabilization of the open state. Critically, we found that the unique and highly conserved extracellular turret region of Kv3.1 and Kv3.2 essentially governs AUT5 modulation. Furthermore, leveraging on the cryo-EM structure of Kv3.1a, atomistic blind docking calculations revealed four equivalent AUT5 binding sites near the turrets and between the voltage-sensing and pore domains of the channel’s tetrameric assembly. Therefore, the unique Kv3 turret emerges as a novel structural correlate of the selective MoA of a new class of Kv3 channel PAMs with a therapeutic potential.

Tutorial project for Allosteric Modulation of Membrane Proteins by Small Low-Affinity Ligands. The project contains a complete set of scripts for analysis and molecular configurations for flooding MD simulations of the sevoflurane/Kv1.2 system., The work was supported by National Council for Scientific and Technological Development CNPq [grant number 302089/2019-5 and 200114/2020-4], Coordenação de Aperfeiçoamento de Pessoal de Nível Superior CAPES [grant number 23038.010052/2013-95], and Fundação de Apoio a Pesquisa do Distrito Federal FAPDF [grant number 193.001.202/2016].

The problem of finding the correct set of partners for a given pair of interacting protein families based on multi-sequence alignments (MSAs) has received great attention over the years. Recently, the native contacts of two interacting proteins were shown to store the strongest mutual information (MI) signal to discriminate MSA concatenations with the largest fraction of correct pairings. Although that signal might be of practical relevance in the search for an effective heuristic to solve the problem, the number of MSA concatenations with near-native MI is large, imposing severe limitations. Here, a Genetic Algorithm that explores possible MSA concatenations according to a MI maximization criteria is shown to find degenerate solutions with two error sources, arising from mismatches among (i) similar and (ii) non-similar sequences. If mistakes made among similar sequences are disregarded, type-(i) solutions are found to resolve correct pairings at best true positive (TP) rates of 70%—far above the very same estimates in type-(ii) solutions. A machine learning classification algorithm helps to show further that differences between optimized solutions based on TP rates are not artificial and may have biological meaning associated with the three-dimensional distribution of the MI signal. Type-(i) solutions may therefore correspond to reliable results for predictive purposes, found here to be more likely obtained via MI maximization across protein systems having a minimum critical number of amino acid contacts on their interaction surfaces (N > 200).

Graphical abstract, Here, we investigate the contributions of coevolutive, evolutive and stochastic information in determining protein-protein interactions (PPIs) based on primary sequences of two interacting protein families A and B. Specifically, under the assumption that coevolutive information is imprinted on the interacting amino acids of two proteins in contrast to other (evolutive and stochastic) sources spread over their sequences, we dissect those contributions in terms of compensatory mutations at physically-coupled and uncoupled amino acids of A and B. We find that physically-coupled amino-acids at short range distances store the largest per-contact mutual information content, with a significant fraction of that content resulting from coevolutive sources alone. The information stored in coupled amino acids is shown further to discriminate multi-sequence alignments (MSAs) with the largest expectation fraction of PPI matches – a conclusion that holds against various definitions of intermolecular contacts and binding modes. When compared to the informational content resulting from evolution at long-range interactions, the mutual information in physically-coupled amino-acids is the strongest signal to distinguish PPIs derived from cospeciation and likely, the unique indication in case of molecular coevolution in independent genomes as the evolutive information must vanish for uncorrelated proteins.

Donepezil is a second generation acetylcholinesterase (AChE) inhibitor for treatment of Alzheimer's disease (AD). AChE is important for neurotransmission at neuromuscular junctions and cholinergic brain synapses by hydrolyzing acetylcholine into acetate and choline. In vitro data support that donepezil is a reversible, mixed competitive and noncompetitive inhibitor of AChE. The experimental fact then suggests a more complex binding mechanism beyond the molecular view in X-ray models resolved at cryogenic temperatures that show a unique binding mode of donepezil in the active site of the enzyme. Aiming at clarifying the mechanism behind that mixed competitive and noncompetitive nature of the inhibitor, we have applied molecular dynamics (MD) simulations and docking and free-energy calculations to investigate microscopic details and energetics of donepezil association for conditions of substrate-free and -bound states of the enzyme. Liquid-phase MD simulation at room temperature shows AChE transits between "open" and "closed" conformations to control accessibility to the active site and ligand binding. As shown by docking and free-energy calculations, association of donepezil involves its reversible axial displacement and reorientation in the active site of the enzyme, assisted by water molecules. Donepezil binds equally well the main-door anionic binding site PAS, the acyl pocket, and the catalytic site CAS by respectively adopting outward-inward-inward orientations regardless of substrate occupancy-the overall stability of that reaction process depends however on co-occupancy of the enzyme being preferential for its substrate-free state. All together, our findings support a physiologically relevant mechanism of AChE inhibition by donepezil involving multistable interactions modes at the molecular origin of the inhibitor's activity.

Changes in chromatin structure regulate gene expression and genome maintenance. Molecules that bind to the nucleosome, the complex of DNA and histone proteins, are key modulators of chromatin structure. Previous work indicated that cholesterol, a ubiquitous cellular lipid, may bind to chromatin in vivo, suggesting a potential function for lipids in modulating chromatin architecture. However, the molecular mechanisms of cholesterol's action on chromatin structure have remained unclear. Here, we explored the biophysical impact of cholesterol on nucleosome and chromatin fibers reconstituted in vitro and characterized in silico the cholesterol binding to the nucleosome. Our findings support that cholesterol assists 10 and 30 nm chromatin formation and induces folding of long chromatin fibers as a result of direct interaction of the cholesterol to six nucleosomal binding sites.

The direct-site hypothesis assumes general anesthetics bind ion channels to impact protein equilibrium and function, inducing anesthesia. Despite advancements in the field, a first principle all-atom demonstration of this structure-function premise is still missing. We focus on the clinically used sevoflurane interaction to anesthetic-sensitive Kv1.2 mammalian channel to resolve if sevoflurane binds protein’s well-characterized open and closed structures in a conformation-dependent manner to shift channel equilibrium. We employ an innovative approach relying on extensive docking calculations and free-energy perturbation of all potential binding sites revealed by the latter, and find sevoflurane binds open and closed structures at multiple sites under complex saturation and concentration effects. Results point to a non-trivial interplay of site and conformation-dependent modes of action involving distinct binding sites that increase channel open-probability at diluted ligand concentrations. Given the challenge in exploring more complex processes potentially impacting channel-anesthetic interaction, the result is revealing as it demonstrates the process of multiple anesthetic binding events alone may account for open-probability shifts recorded in measurements., Author summary General anesthetics are central to modern medicine, yet their microscopic mechanism of action is still unknown. Here, we demonstrate that a clinically used anesthetic, sevoflurane, binds the mammalian voltage-gated potassium channel Kv1.2 effecting a shift in its open probability, even at low concentrations. The results, supported by recent experimental measurements, are promising as they demonstrate that the molecular process of direct binding of anesthetic to ion channels play a relevant role in anesthesia.

2,2,2-Trichloroethanol (TCE) is the active form of the sedative hypnotic drug chloral hydrate, one of the oldest sleep medications in the market. Understanding of TCE's action mechanisms to its many targets, particularly within the ion channel family, could benefit from the state-of-the-art computational molecular studies. In this direction, we employed de novo modeling aided by the force field toolkit to develop CHARMM36-compatible TCE parameters. The classical potential energy function was calibrated targeting molecular conformations, local interactions with water molecules, and liquid bulk properties. Reference data comes from both tabulated thermodynamic properties and ab initio calculations at the MP2 level. TCE solvation free energy calculations in water and oil reproduce a lipophilic, yet nonhydrophobic, behavior. Indeed, the potential mean force profile for TCE partition through the phospholipid bilayer reveals the sedative's preference for the interfacial region. The calculated partition coefficient also matches experimental measures. Further validation of the proposed parameters is supported by the model's ability to recapitulate quenching experiments demonstrating TCE binding to bovine serum albumin.

The direct-site hypothesis assumes general anesthetics bind ion channels to impact protein equilibrium and function, inducing anesthesia. Despite advancements in the field, a first-principle all-atom demonstration of this structure-function premise misses. We focus on clinically used sevoflurane interaction to anesthetic-sensitive Kv1.2 mammalian channel to resolve if sevoflurane binds the protein’s well-characterized open and closed structures in a conformation-dependent manner to shift channel equilibrium. We employ an innovative approach relying on extensive docking calculations and free-energy perturbation and find sevoflurane binds open and closed structures at multiple sites under complex saturation and concentration effects. Results point to a non-trivial interplay of conformation-dependent modes of action involving distinct binding sites that increase channel open-probability at diluted ligand concentrations. Given the challenge in exploring more complex processes potentially impacting channel-anesthetic interaction, the result is reassuring as demonstrates that the process of multiple binding events alone may account for open-probability shifts recorded in measurements.

Structural changes in chromatin regulate gene expression and define phenotypic outcomes. Molecules that bind to the nucleosome, the complex of DNA and histone proteins, are key modulators of chromatin structure. Most recently, the formation of condensed chromatin regions based on phase-separation in the cell, a basic physical mechanism, was proposed. Increased understanding of the mechanisms of interaction between chromatin and lipids suggest that small lipid molecules, such as cholesterol and short-chain fatty acids, can regulate important nuclear functions. New biophysical data has suggested that cholesterol interacts with nucleosome through multiple binding sites and affects chromatin structure in vitro. Regardless of the mechanism of how lipids bind to chromatin, there is currently little awareness that lipids may be stored in chromatin and influence its state. Focusing on lipids that bind to nuclear receptors, clinically relevant transcription factors, we discuss the potential interactions of the nucleosome with steroid hormones, bile acids and fatty acids, which suggest that other lipid chemotypes may also impact chromatin structure through binding to common sites on the nucleosome. Herein, we review the main impacts of lipids on the nuclear environment, emphasizing its role on chromatin architecture. We postulate that lipids that bind to nucleosomes and affect chromatin states are likely to be worth investigating as tools to modify disease phenotypes at a molecular level.

Membrane proteins are primary targets for most therapeutic indications in cancer and neurological diseases, binding over 50% of all known small molecule drugs. Understanding how such ligands impact membrane proteins requires knowledge on the molecular structure of ligand binding, a reasoning that has driven relentless efforts in drug discovery and translational research. Binding of small ligands appears however highly complex involving interaction to multiple transmembrane protein sites featuring single or multiple occupancy states. Within this scenario, looking for new developments in the field, we investigate the concentration-dependent binding of ligands to multiple saturable sites in membrane proteins. The study relying on docking and free-energy perturbation provides us with an extensive description of the probability density of protein-ligand states that allows for computation of thermodynamic properties of interest. It also provides one- and three-dimensional spatial descriptions for the ligand density across the protein-membrane system which can be of interest for structural purposes. Illustration and discussion of the results are shown for binding of the general anesthetic sevoflurane against Kv1.2, a mammalian ion channel for which experimental data are available.

A dynamic transmembrane voltage field has been suggested as an intrinsic element in voltage sensor (VS) domains. Here, the dynamic field contribution to the VS energetics was analyzed via electrostatic calculations applied to a number of atomistic structures made available recently. We find that the field is largely static along with the molecular motions of the domain, and more importantly, it is minimally modified across VS variants. This finding implies that sensor domains transfer approximately the same amount of gating charges when moving the electrically charged S4 helix between fixed microscopic configurations. Remarkably, the result means that the observed operational diversity of the domain, including the extension, rate, and voltage dependence of the S4 motion, as dictated by the free energy landscape theory, must be rationalized in terms of dominant variations of its chemical free energy.

The X-ray structure of the bacterial voltage-gated sodium channel NavAb has been reported in a conformation with a closed conduction pore. Comparison between this structure and the activated-open and resting-closed structures of the voltage-gated Kv1.2 potassium channel suggests that the voltage-sensor domains (VSDs) of the reported structure are not fully activated. Using the aforementioned structures of Kv1.2 as templates, molecular dynamics simulations are used to identify analogous functional conformations of NavAb. Specifically, starting from the NavAb crystal structure, conformations of the membrane-bound channel are sampled along likely pathways for activation of the VSD and opening of the pore domain. Gating charge computations suggest that a structural rearrangement comparable to that occurring between activated-open and resting-closed states is required to explain experimental values of the gating charge, thereby confirming that the reported VSD structure is likely an intermediate along the channel activation pathway. Our observation that the X-ray structure exhibits a low pore domain-opening propensity further supports this notion. The present molecular dynamics study also identifies conformations of NavAb that are seemingly related to the resting-closed and activated-open states. Our findings are consistent with recent structural and functional studies of the orthologous channels NavRh, NaChBac, and NavMs and offer possible structures for the functionally relevant conformations of NavAb.

Computational methods and experimental data are used to provide structural models for NaChBac, the homo-tetrameric voltage-gated sodium channel from the bacterium Bacillus halodurans, with a closed and partially open pore domain. Molecular dynamic (MD) simulations on membrane-bound homo-tetrameric NaChBac structures, each comprising six helical transmembrane segments (labeled S1 through S6), reveal that the shape of the lumen, which is defined by the bundle of four alpha-helical S6 segments, is modulated by hinge bending motions around the S6 glycine residues. Mutation of these glycine residues into proline and alanine affects, respectively, the structure and conformational flexibility of the S6 bundle. In the closed channel conformation, a cluster of stacked phenylalanine residues from the four S6 helices hinders diffusion of water molecules and Na+ ions. Activation of the voltage sensor domains causes destabilization of the aforementioned cluster of phenylalanines, leading to a more open structure. The conformational change involving the phenylalanine cluster promotes a kink in S6, suggesting that channel gating likely results from the combined action of hinge-bending motions of the S6 bundle and concerted reorientation of the aromatic phenylalanine side-chains.

Molecular dynamics (MD) calculations have been used to identify Na+ binding sites in the lumen of a bacterial voltage-gated ion channel. MD trajectories have been carried out starting from the rece...

International audience; Most studies of the interactions of neat and functionalized fullerenes with cells have focused so far on their ability to cross the cell membrane envelopes. Membranes are, however, also host to a large number of proteins responsible for various cellular functions. Among these, ion channels are prominent components of the nervous system. Recently, it was shown that fullerenes may act as blockers or modulators of a variety of K+ channels. Here we use computer simulations to investigate the propensity of such nanocompounds to bind to K+ channels. Our results based on extensive atomistic molecular dynamics simulations reveal a variety of specific binding sites depending on the structure and properties of the channel. The corresponding binding free energies and putative mechanisms suggest that C60 may indeed effectively hinder the function of K+ channels and hence induce toxicity.

Potassium channel subunits composed of two-pore domains arranged in tandem (K(2P)) are of paramount importance for neural function. A variety of stimuli, such as membrane depolarization and tension, acidification, and anesthetic action, activate K(2P) channels. Most of the channel sensitivity is attributed to its intracellular C-terminal moiety, which works as a sensor domain required for proper integration of the electrical, chemical, and mechanical signals into channel activity. Herein, the structure of K(2P) in a membrane environment has been studied using molecular dynamics (MD). Two distinct fully atomistic models for the most studied K(2P) channel, namely, the TWIK-related (TREK)-1 channel have been built. These constructs were then inserted into a fully hydrated zwitterionic lipid bilayer, and each relaxed by means of MD simulations spanning approximately 0.3 micros. Both simulated TREK-1 structures converged to a final conformation characterized by a closed pore and a C-terminal domain adsorbed onto the lipid bilayer surface. The C-terminus, which is physically linked to the pore and energetically coupled to the bilayer, is poised to gate the channel in response to membrane stimulation. The present study indicates the nature of the direct coupling between the C-terminal domain and the membrane, which is a key structural feature underlying K(2P) channel function.

International audience; Coarse grain modeling has recently emerged as an alternative to classical atomistic simulations in the study of spontaneous self assembly and structural organization of complex molecular systems. For surfactant and lipid systems, it was shown to allow, under appropriate conditions, in-silico self assembly of a variety of architectures. Recently, this approach has been extended to peptides for which force fields allowing self assembly of mixed peptides–lipid systems were proposed. Here we introduce elements of a coarse grained force field that accurately describe self assembly of hydrophobic cyclic peptides [Trp–Leu]4 and their reorganization within lipid membranes to form transmembrane channels in agreement with experiments. Extension to hydrophobic helical transmembrane, and amphipatic helical antimicrobial peptides show that the model is robust enough to constitute a building block for more complete and appropriate force field describing peptide interactions with membranes.

General anesthesia is a relatively safe medical procedure, which for nearly 170 years has allowed life saving surgical interventions in animals and people. However, the molecular mechanism of general anesthesia continues to be a matter of importance and debate. A favored hypothesis proposes that general anesthesia results from direct multisite interactions with multiple and diverse ion channels in the brain. Neurotransmitter-gated ion channels and two-pore K+ channels are key players in the mechanism of anesthesia; however, new studies have also implicated voltage-gated ion channels. Recent biophysical and structural studies of Na+ and K+ channels strongly suggest that halogenated inhalational general anesthetics interact with gates and pore regions of these ion channels to modulate function. Here, we review these studies and provide a perspective to stimulate further advances.

A high-resolution crystal structure of KvAP, an archeabacterial voltage-gated potassium (Kv) channel, complexed with a monoclonal Fab fragment has been recently determined. Based on this structure, a mechanism for the activation (opening) of Kv channels has been put forward. This mechanism has since been criticized, suggesting that the resolved structure is not representative of the family of voltage-gated potassium channels. Here, we propose a model of the transmembrane domain of Shaker B, a well-characterized Kv channel, built by homology modeling and docking calculations. In this model, the positively charged S4 helices are oriented perpendicular to the membrane and localized in the groove between segments S5 and S6 of adjacent subunits. The structure and the dynamics of the full atomistic model embedded in a hydrated lipid bilayer were investigated by means of two large-scale molecular dynamics simulations under transmembrane-voltage conditions known to induce, respectively, the resting state (closed) and the activation (opening) of voltage-gated channels. Upon activation, the model undergoes conformational changes that lead to an increase of the hydration of the charged S4 helices, correlated with an upward translation and a tilting of the latter, concurrently with movements of the S5 helices and the activation gate. Although small, these conformational changes ultimately result in an alteration of the ion-conduction pathway. Our findings support the transporter model devised by Bezanilla and collaborators, and further underline the crucial role played by internal hydration in the activation of the channel.

By Monte Carlo simulations, we explored the effect of single mutations on the thermodynamics and kinetics of the folding of a two-dimensional, energetically frustrated, hydrophobic protein model. Φ-Value analysis, corroborated by simulations beginning from given sets of judiciously chosen initial contacts, suggests that the transition state of the model consists of a limited region of the native structure, that is, a folding nucleus. It seems that the most important contacts in the transition state (large and positive Φ) are not the ones with the highest contact order, because in this case the entropic cost of their formation would be too high, but exactly the ones that decrease the entropic cost of difficult contacts, reducing their effective contact order. Mutations of internal monomers involved in high-order contacts were actually the ones resulting in the fastest kinetics (and Φ < 0), indicating they tend to make low order, non-native contacts of low entropic cost that stabilize the unfolded state with respect to the transition state. Folding acceleration by other non-native interactions was also observed and a simple general mechanism is proposed according to which non-native contacts can act indirectly over the folding nucleus, “chelating” out potentially harmful contacts. The polymer graph of our model, which facilitates the visualization of effective contact orders, successfully suggests the relative kinetic importance of different contacts and is reasonably consistent with analogous graphs for the well characterized family of SH3 domains. Proteins 2002;49:167–180. © 2002 Wiley-Liss, Inc.

The effects on the structural and functional properties of the Kv1.2 voltage-gated ion channel, caused by selective mutation of voltage sensor domain residues, have been investigated using classical molecular dynamics simulations. Following experiments that have identified mutations of voltage-gated ion channels involved in state-dependent omega currents, we observe for both the open and closed conformations of the Kv1.2 that specific mutations of S4 gating-charge residues destabilize the electrostatic network between helices of the voltage sensor domain, resulting in the formation of hydrophilic pathways linking the intra- and extracellular media. When such mutant channels are subject to transmembrane potentials, they conduct cations via these so-called “omega pores.” This study provides therefore further insight into the molecular mechanisms that lead to omega currents, which have been linked to certain channelopathies.

Ion channels assist and control the diffusion of ions through biological membranes. The conduction process depends on the structural characteristics of these nanopores, among which are the hydrophobicity and the afforded diameter of the conduction pathway. In this contribution, we use full atomistic free-energy molecular dynamics simulations to estimate the effect of such characteristics on the energetics of ion conduction through the activation gate of voltage-gated potassium (Kv) channels. We consider specifically the ionic translocation through three different permeation pathways, corresponding to the activation gate of an atomistic model of Shaker channels in closed and partially opened conformations, and that of the open conformation of the Kv1.2 channel. In agreement with experiments, we find that the region of Val478 constitutes the main gate. The conduction is unfavorable through this gate when the constriction is smaller than an estimated threshold of 4.5–5.0Å, mainly due to incomplete coordination-hydration of the ion. Above this critical size, e.g., for the Kv1.2, the valine gate is wide enough to allow fully coordination of the ion and therefore its diffusion on a flat energy surface. Similar to other ion channels, Kv channels appear therefore to regulate diffusion by constricting hydrophobic regions of the permeation pathway.

Recently, the structure of the Shaker channel Kv1.2 has been determined at a 2.9-angstroms resolution. This opens new possibilities in deciphering the mechanism underlying the function of voltage-gated potassium (Kv) channels. Molecular dynamics simulations of the channel, embedded in a membrane environment show that the channel is in its open state and that the gating charges carried by S4 are exposed to the solvent. The hydrated environment of S4 favors a local collapse of the electrostatic potential, which generates high electric-field gradients around the arginine gating charges. Comparison to experiments suggests furthermore that activation of the channel requires mainly a lateral displacement of S4. Overall, the results agree with the transporter model devised for Kv channels from electrophysiology experiments, and provide a possible pathway for the mechanistic response to membrane depolarization.

The recent crystal structure of Orai, the pore unit of a calcium release-activated calcium (CRAC) channel, is used as the starting point for molecular dynamics and free-energy calculations designed to probe this channel's conduction properties. In free molecular dynamics simulations, cations localize preferentially at the extracellular channel entrance near the ring of Glu residues identified in the crystal structure, whereas anions localize in the basic intracellular half of the pore. To begin to understand ion permeation, the potential of mean force (PMF) was calculated for displacing a single Na(+) ion along the pore of the CRAC channel. The computed PMF indicates that the central hydrophobic region provides the major hindrance for ion diffusion along the permeation pathway, thereby illustrating the nonconducting nature of the crystal structure conformation. Strikingly, further PMF calculations demonstrate that the mutation V174A decreases the free energy barrier for conduction, rendering the channel effectively open. This seemingly dramatic effect of mutating a nonpolar residue for a smaller nonpolar residue in the pore hydrophobic region suggests an important role for the latter in conduction. Indeed, our computations show that even without significant channel-gating motions, a subtle change in the number of pore waters is sufficient to reshape the local electrostatic field and modulate the energetics of conduction, a result that rationalizes recent experimental findings. The present work suggests the activation mechanism for the wild-type CRAC channel is likely regulated by the number of pore waters and hence pore hydration governs the conductance.

The crystal structure of NavAb, a bacterial voltage gated Na(+) channel, exhibits a selectivity filter (SF) wider than that of K(+) channels. This new structure provides the opportunity to explore the mechanism of conduction and help rationalize its selectivity for sodium. Recent molecular dynamics (MD) simulations of single- and two-ion permeation processes have revealed that a partially hydrated Na(+) permeates the channel by exploring three SF binding sites while being loosely coupled to other ions and/or water molecules; a finding that differs significantly from the behavior of K(+) selective channels. Herein, we present results derived from a combination of metadynamics and voltage-biased MD simulations that throws more light on the nature of the Na(+) conduction mechanism. Conduction under 0 mV bias explores several distinct pathways involving the binding of two ions to three possible SF sites. While these pathways are very similar to those observed in the presence of a negative potential (inward conduction), a completely different mechanism operates for outward conduction at positive potentials.

Because of their remarkable roles in electrical cell signaling, voltage-gated cation channels (VGCCs) have been the subject of intense investigations and debate for more than 50 years. Ultimately, the prospective implications of such studies have an impact on our understanding of the molecular properties of VGCCs involved in consciousness, anesthesia, and diseases, to mention a few. The following review aims to summarize our current knowledge of activation of VGCCs by highlighting major methodological innovations in the field and the breakthroughs they allowed. Focusing mainly on insights gained through computer simulations, while acknowledging important experimental findings, we hope to inspire experimentalists to benefit from these approaches in the generation of hypotheses and design of experiments. Also, we outline major future challenges for the field, such as channel modulation, lesser-known receptors, and molecular origins of channel dysfunctions.

Voltage-sensor domains (VSDs) are electrically-charged constructs controlling the voltage-dependent activity of ion channels in excitable cells. Four packed transmembrane (TM) helices, S1 through S4, form the domain in which S4 contains 4 to 7 positively charged basic amino acids, mostly arginines. VSDs operate essentially by transferring the S4 charges across the transmembrane electric field (E), giving rise to the observable Q, the so-called “gating charge”. Mostly supported by structure-function studies on voltage-gated potassium (Kv) channels, a focused E has been identified as one key electric property of the VSD machinery. The recent increasing availability of other VSD-containing ion channel structures, including the x-ray structures for the NavAb and NavRh voltage-gated Na+ channels, provides us with the opportunity to extend the structure-based investigation of the domain electrostatic properties over a larger set of distinct conformations and isoforms. Using all-atom MD simulations in combination with electrostatic calculations, founded on an energetic formalism, we show that, over the entire set of available VSD structures, a specific hydration of the voltage sensor focuses E over a narrow TM region across the domain, at the vicinity of the so-called catalytic center. Furthermore, its focalization and shape is largely preserved over distinct conformations of the construct. Our results support that a focused and conformation-independent TM field is a robust electric feature of the VSD machinery, despite sequence variations or local structural modifications of the domain. This electric fingerprint seems to favor a highly conserved sensing mechanism for VSDs over the large family of voltage-gated cation channels.

The relentless growth in computational power has seen increasing applications of molecular dynamics (MD) simulation to the study of membrane proteins in realistic membrane environments, which include explicit membrane lipids, water and ions. The concomitant increasing availability of membrane protein structures for ion channels, and transporters -- to name just two examples -- has stimulated many of these MD studies. In the case of voltage-gated cation channels (VGCCs) recent computational works have focused on ion-conduction and gating mechanisms, along with their regulation by agonist/antagonist ligands. The information garnered from these computational studies is largely inaccessible to experiment and is crucial for understanding the interplay between the structure and function as well as providing new directions for experiments. This article highlights recent advances in probing the structure and function of potassium channels and offers a perspective on the challenges likely to arise in making analogous progress in characterizing sodium channels.