Search

Many voltage-gated potassium (Kv) channels display a time-dependent phenomenon called C-type inactivation, whereby prolonged activation by voltage leads to the inhibition of ionic conduction, a process that involves a conformational change at the selectivity filter toward a non-conductive state. Recently, a high-resolution structure of a strongly inactivated triple-mutant channel kv1.2-kv2.1-3m revealed a novel conformation of the selectivity filter that is dilated at its outer end, distinct from the well-characterized conductive state. While the experimental structure was interpreted as the elusive non-conductive state, our molecular dynamics simulations and electrophysiological measurements show that the dilated filter of kv1.2-kv2.1-3m is conductive and, as such, cannot completely account for the inactivation of the channel observed in the structural experiments. The simulation shows that an additional conformational change, implicating isoleucine residues at position 398 along the pore lining segment S6, is required to effectively block ion conduction. The I398 residues from the four subunits act as a state-dependent hydrophobic gate located immediately beneath the selectivity filter. These observations are corroborated by electrophysiological experiments showing that ion permeation can be resumed in the kv1.2-kv2.1-3m channel when I398 is mutated to an asparagine-a mutation that does not abolish C-type inactivation since digitoxin (AgTxII) fails to block the ionic permeation of kv1.2-kv2.1-3m_I398N. As a critical piece of the C-type inactivation machinery, this structural feature is the potential target of a broad class of quaternary ammonium (QA) blockers and negatively charged activators thus opening new research directions toward the development of drugs that specifically modulate gating states of Kv channels., Competing Interests: WT, YL, CB, BP, JA, RU, CC, FB, BR No competing interests declared, (© 2024, Treptow, Liu, Bassetto et al.)

Maurotoxin (MTX) is a 34-residue peptide from Scorpio maurus venom. It is reticulated by four disulfide bridges with a unique arrangement compared to other scorpion toxins that target potassium (K + ) channels. Structure-activity relationship studies have not been well performed for this toxin family. The screening of Scorpio maurus venom was performed by different steps of fractionation, followed by the ELISA test, using MTX antibodies, to isolate an MTX-like peptide. In vitro, in vivo and computational studies were performed to study the structure-activity relationship of the new isolated peptide. We isolated a new peptide designated MTX1, structurally related to MTX. It demonstrated toxicity on mice eight times more effectively than MTX. MTX1 blocks the Kv1.2 and Kv1.3 channels, expressed in Xenopus oocytes, with IC 50 values of 0.26 and 180 nM, respectively. Moreover, MTX1 competitively interacts with both 125 I-apamin (IC 50 = 1.7 nM) and 125 I-charybdotoxin (IC 50 = 5 nM) for binding to rat brain synaptosomes. Despite its high sequence similarity (85%) to MTX, MTX1 exhibits a higher binding affinity towards the Kv1.2 and SKCa channels. Computational analysis highlights the significance of specific residues in the β-sheet region, particularly the R27, in enhancing the binding affinity of MTX1 towards the Kv1.2 and SKCa channels.

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.

Shaker potassium channels have been an essential model for studying inactivation of ion channels and shaped our earliest understanding of N-type vs. C-type mechanisms. In early work describing C-type inactivation, López-Barneo and colleagues systematically characterized numerous mutations of Shaker residue T449, demonstrating that this position was a key determinant of C-type inactivation rate. In most of the closely related mammalian Kv1 channels, however, a persistent enigma has been that residue identity at this position has relatively modest effects on the rate of inactivation in response to long depolarizations. In this study, we report alternative ways to measure or elicit conformational changes in the outer pore associated with C-type inactivation. Using a strategically substituted cysteine in the outer pore, we demonstrate that mutation of Kv1.2 V381 (equivalent to Shaker T449) or W366 (Shaker W434) markedly increases susceptibility to modification by extracellularly applied MTSET. Moreover, due to the cooperative nature of C-type inactivation, Kv1.2 assembly in heteromeric channels markedly inhibits MTSET modification of this substituted cysteine in neighboring subunits. The identity of Kv1.2 residue V381 also markedly influences function in conditions that bias channels toward C-type inactivation, namely when Na + is substituted for K + as the permeant ion or when channels are blocked by an N-type inactivation particle (such as Kvβ1.2). Overall, our findings illustrate that in mammalian Kv1 channels, the identity of the T449-equivalent residue can strongly influence function in certain experimental conditions, even while having modest effects on apparent inactivation during long depolarizations. These findings contribute to reconciling differences in experimental outcomes in many Kv1 channels vs. Shaker., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Thermally stable full-length scorpion toxin peptides and partially degraded peptides with complete disulfide bond pairing are valuable natural peptide resources in traditional Chinese scorpion medicinal material. However, their pharmacological activities are largely unknown. This study discovered BmKcug1a-P1, a novel N-terminal degraded peptide, in this medicinal material. BmKcug1a-P1 inhibited hKv1.2 and hKv1.3 potassium channels with IC 50 values of 2.12 ± 0.27 μM and 1.54 ± 0.28 μM, respectively. To investigate the influence of N-terminal amino acid loss on the potassium channel inhibiting activities, three analogs (i.e., full-length BmKcug1a, BmKcug1a-P1-D2 and BmKcug1a-P1-D4) of BmKcug1a-P1 were prepared, and their potassium channel inhibiting activities on hKv1.3 channel were verified by whole-cell patch clamp technique. Interestingly, the potassium channel inhibiting activity of full-length BmKcug1a on the hKv1.3 channel was significantly improved compared to its N-terminal degraded form (BmKcug1a-P1), while the activities of two truncated analogs (i.e., BmKcug1a-P1-D2 and BmKcug1a-P1-D4) were similar to that of BmKcug1a-P1. Extensive alanine-scanning experiments identified the bonding interface (including two key functional residues, Asn30 and Arg34) of BmKcug1a-P1. Structural and functional dissection further elucidated whether N-terminal residues of the peptide are located at the bonding interface is important in determining whether the N-terminus significantly influences the potassium channel inhibiting activity of the peptide. Altogether, this research identified a novel N-terminal degraded active peptide, BmKcug1a-P1, from traditional Chinese scorpion medicinal material and elucidated how the N-terminus of peptides influences their potassium channel inhibiting activity, contributing to the functional identification and molecular truncation optimization of full-length and degraded peptides from traditional Chinese scorpion medicinal material Buthus martensii Karsch., (© 2024. The Author(s).)

Membrane proteins play critical physiological roles as receptors, channels, pumps, and transporters. Despite their importance, however, low expression levels often hamper the experimental characterization of membrane proteins. We present an automated and web-accessible design algorithm called mPROSS (https://mPROSS.weizmann.ac.il), which uses phylogenetic analysis and an atomistic potential, including an empirical lipophilicity scale, to improve native-state energy. As a stringent test, we apply mPROSS to the Kv1.2-Kv2.1 paddle chimera voltage-gated potassium channel. Four designs, encoding 9-26 mutations relative to the parental channel, were functional and maintained potassium-selective permeation and voltage dependence in Xenopus oocytes with up to 14-fold increase in whole-cell current densities. Additionally, single-channel recordings reveal no significant change in the channel-opening probability nor in unitary conductance, indicating that functional expression levels increase without impacting the activity profile of individual channels. Our results suggest that the expression levels of other dynamic channels and receptors may be enhanced through one-shot design calculations., (© 2024 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

In this paper, the influence of external terahertz electromagnetic fields with different frequencies of 4 THz, 10 THz, 15 THz, and 20 THz on the permeability of the Kv1.2 voltage-gated potassium ion channel on the nerve cell membrane was studied using the combined model of the "Constant Electric Field-Ion Imbalance" method by molecular dynamics. We found that although the applied terahertz electric field does not produce strong resonance with the -C=O groups of the conservative sequence T-V-G-Y-G amino acid residue of the selective filter (SF) of the channel, it would affect the stability of the electrostatic bond between potassium ions and the carbonyl group of T-V-G-Y-G of SF, and it would affect the stability of the hydrogen bond between water molecules and oxygen atoms of the hydroxyl group of the 374THR side chain at the SF entrance, changing the potential and occupied states of ions in the SF and the occurrence probability of the permeation mode of ions and resulting in the change in the permeability of the channel. Compared with no external electric field, when the external electric field with 15 THz frequency is applied, the lifetime of the hydrogen bond is reduced by 29%, the probability of the "soft knock on" mode is decreased by 46.9%, and the ion flux of the channel is activated by 67.7%. Our research results support the view that compared to "direct knock-on", "soft knock-on" is a slower permeation mode.

Kv1.2 channels, encoded by the KCNA2 gene, are localized in the central and peripheral nervous system, where they regulate neuronal excitability. Recently, heterozygous mutations in KCNA2 have been associated with a spectrum of symptoms extending from epileptic encephalopathy, intellectual disability, and cerebellar ataxia. Patients are treated with a combination of antiepileptic drugs and 4-aminopyridine (4-AP) has been recently trialed in specific cases. We identified a novel variant in KCNA2 , E236K, in a Serbian proband with non-progressive congenital ataxia and early onset epilepsy, treated with sodium valproate. To ascertain the pathogenicity of E236K mutation and to verify its sensitivity to 4-AP, we transfected HEK 293 cells with Kv1.2 WT or E236K cDNAs and recorded potassium currents through the whole-cell patch-clamp. In silico analysis supported the electrophysiological data. E236K channels showed voltage-dependent activation shifted towards negative potentials and slower kinetics of deactivation and activation compared with Kv1.2 WT. Heteromeric Kv1.2 WT+E236K channels, resembling the condition of the heterozygous patient, confirmed a mixed gain- and loss-of-function (GoF/LoF) biophysical phenotype. 4-AP inhibited both Kv1.2 and E236K channels with similar potency. Homology modeling studies of mutant channels suggested a reduced interaction between the residue K236 in the S2 segment and the gating charges at S4. Overall, the biophysical phenotype of E236K channels correlates with the mild end of the clinical spectrum reported in patients with GoF/LoF defects. The response to 4-AP corroborates existing evidence that KCNA2 -disorders could benefit from variant-tailored therapeutic approaches, based on functional studies.

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τ fold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization., Competing Interests: The authors declare no competing interest.

Graphene quantum dots (GQDs) are emerging as a versatile nanomaterial with numerous proposed biomedical applications. Despite the explosion in potential applications, the molecular interactions between GQDs and complex biomolecular systems, including potassium-ion (K + ) channels, remain largely unknown. Here, we use molecular dynamics (MD) simulations and electrophysiology to study the interactions between GQDs and three representative K + channels, which participate in a variety of physiological processes and are closely related to many disease states. Using MD simulations, we observed that GQDs adopt distinct contact poses with each of the three structurally distinct K + channels. Our electrophysiological characterization of the effects of GQDs on channel currents revealed that GQDs interact with the extracellular voltage-sensing domain (VSD) of a Kv1.2 channel, augmenting current by left-shifting the voltage dependence of channel activation. In contrast, GQDs form a "lid" cluster over the extracellular mouth of inward rectifier Kir3.2, blocking the channel pore and decreasing the current in a concentration-dependent manner. Meanwhile, GQDs accumulate on the extracellular "cap domain" of K2P2 channels and have no apparent impact on the K + flux through the channel. These results reveal a surprising multifaceted regulation of K + channels by GQDs, which might help de novo design of nanomaterial-based channel probe openers/inhibitors that can be used to further discern channel function.

Understanding the gating mechanism of ion channel proteins is key to understanding the regulation of cell signaling through these channels. Channel opening and closing are regulated by diverse environmental factors that include temperature, electrical voltage across the channel, and proton concentration. Low permeability in voltage-gated potassium ion channels (Kv) is intimately correlated with the prolonged action potential duration observed in many acidosis diseases. The Kv channels consist of voltage-sensing domains (S1-S4 helices) and central pore domains (S5-S6 helices) that include a selectivity filter and water-filled cavity. The voltage-sensing domain is responsible for the voltage-gating of Kv channels. While the low permeability of Kv channels to potassium ion is highly correlated with the cellular proton concentration, it is unclear how an intracellular acidic condition drives their closure, which may indicate an additional pH-dependent gating mechanism of the Kv family. Here, we show that two residues E327 and H418 in the proximity of the water cavity of Kv1.2 play crucial roles as a pH switch. In addition, we present a structural and molecular concept of the pH-dependent gating of Kv1.2 in atomic detail, showing that the protonation of E327 and H418 disrupts the electrostatic balance around the S6 helices, which leads to a straightening transition in the shape of their axes and causes dewetting of the water-filled cavity and closure of the channel. Our work offers a conceptual advancement to the regulation of the pH-dependent gating of various voltage-gated ion channels and their related biological functions., Competing Interests: The authors declare no competing interests.

Scorpion venom contains diverse bioactive peptides that can recognize and interact with membrane proteins such as ion channels. These natural toxins are believed to be useful tools for exploring the structure and function of ion channels. In this study, we characterized a K + -channel toxin gene, ImKTx96, from the venom gland cDNA library of the scorpion Isometrus maculates. The peptide deduced from the ImKTx96 precursor nucleotide sequence contains a signal peptide of 27 amino acid residues and a mature peptide of 29 residues with three disulfide bridges. Multiple sequence alignment indicated that ImKTx96 is similar with the scorpion toxins that typically target K + -channels. The recombined ImKTx96 peptide (rImKTx96) was expressed in the Escherichia coli system, and purified by GST-affinity chromatography and RP-HPLC. Results from whole-cell patch-clamp experiments revealed that rImKTx96 can inhibit the current of the Kv1.2 ion channel expressed in HEK293 cells. The 3D structure of ImKTx96 was constructed by molecular modeling, and the complex formed by ImKTx96 interacting with the Kv1.2 ion channel was obtained by molecular docking. Based on its structural features and pharmacological functions, ImKTx96 was identified as one member of K + -channel scorpion toxin α-KTx10 group and may be useful as a molecular probe for investigating the structure and function of the Kv1.2 ion channel., (Copyright © 2019. Published by Elsevier Inc.)

Ion channels in cell membranes control entry and exit of ions; their gating (opening and closing) is key to their functioning. It is known that protons can pass through the voltage-sensing domain (VSD) of channels such as Kv1.2. Quantum calculations for a section of the VSD show the steps protons take in responding to voltage and show no major displacement of the protein backbone with voltage change; 70 amino acids are included, 42 with side chains (9 directly in the proton path), 28 as backbone only, and 24 water molecules. Protons provide much of the gating current, the capacitative current immediately preceding channel opening with significant additional contributions from charge transfer to other groups. Most gating models, in contrast, require major protein displacement during gating. Energy terms without classical analogues (exchange plus correlation energy, which are greater than thermal energy) show that quantum calculations are required. Energy as a function of voltage for a key proton transfer leads to, approximately, the correct voltage for channel opening. Calculated total charge transfer (not only protons) for gating is reasonable compared to experimental values. We are also able to account, at least qualitatively, for two mutations, one with the gating current curve left-shifted and one right-shifted, and show the alternate proton paths that are required to account for these.

The Tityus stigmurus scorpion is widely distributed in the Northeast of Brazil and is the main causal agent of human envenoming. The venom produced by this scorpion includes neurotoxins, which are peptides belonging to Family 2 toxins and are able to interact with ion channels. The KTx subfamily displays selectivity and affinity for Kv channel subtypes and the result of this interaction is the blockade of potassium channels, impairing vital functions. We report the optimized structural model of a transcript encoding a potassium channel blocker toxin from T. stigmurus. LC-MS analysis confirmed the presence of the toxin in the venom and the three-dimensional structure was obtained by computational homology modeling and refined by molecular dynamic simulations. Furthermore, docking simulations were performed using a Shaker kV-1.2 potassium channel from rats as receptor model and the contacts were identified revealing which amino acid residues and interactions could be involved in its blockade. These residues were mapped and their contact and electrostatic interactions were evaluated revealing the influence of positive lysine residues and the additional contribution of an asparagine to the stabilization of the complex, bringing new insights into the mechanism of action of this toxin., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Background & Objective: Kunitz-type venoms are bioactive proteins isolated from a wide variety of venomous animals. These venoms are involved in protease inhibitory activity or potassium channel blocking activity. Therefore, they are reported as an important source for lead drug candidates towards protease or channel associated diseases like neurological, metabolic and cardiovascular disorders., Methods: This study aimed to check the inhibitory action of Kunitz-type venoms against potassium channels using computational tools., Results: Among potassium channels, Human Voltage-Gated Potassium Channel 1.2 (hKv1.2) was used as a receptor whereas Kunitz-type peptides from the venoms of various species were selected as ligand dataset., Conclusion: This study helped in finding the binding interface between the receptor and ligand dataset for their potential therapeutic use in treating potassium channelopathies., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Voltage-gated potassium channels play crucial roles in regulating membrane potential. They are activated by membrane depolarization, allowing the selective permeation of K + ions across the plasma membrane, and enter a nonconducting state after lasting depolarization, a process known as inactivation. Inactivation in voltage-activated potassium channels occurs through two distinct mechanisms, N-type and C-type inactivation. C-type inactivation is caused by conformational changes in the extracellular mouth of the channel, whereas N-type inactivation is elicited by changes in the cytoplasmic mouth of the protein. The W434F-mutated Shaker channel is known as a nonconducting mutant and is in a C-type inactivation state at a depolarizing membrane potential. To clarify the structural properties of C-type inactivated protein, we performed molecular dynamics simulations of the wild-type and W366F (corresponding to W434F in Shaker) mutant of the Kv1.2-2.1 chimera channel. The W366F mutant was in a nearly nonconducting state with a depolarizing voltage and recovered from inactivation with a reverse voltage. Our simulations and three-dimensional reference interaction site model analysis suggested that structural changes in the selectivity filter upon membrane depolarization trap K + ions around the inner mouth of the selectivity filter and prevent ion permeation. This pore restriction is involved in the molecular mechanism of C-type inactivation.

Voltage-activated potassium (Kv) channels open to conduct K + ions in response to membrane depolarization, and subsequently enter non-conducting states through distinct mechanisms of inactivation. X-ray structures of detergent-solubilized Kv channels appear to have captured an open state even though a non-conducting C-type inactivated state would predominate in membranes in the absence of a transmembrane voltage. However, structures for a voltage-activated ion channel in a lipid bilayer environment have not yet been reported. Here we report the structure of the Kv1.2-2.1 paddle chimera channel reconstituted into lipid nanodiscs using single-particle cryo-electron microscopy. At a resolution of ~3 Å for the cytosolic domain and ~4 Å for the transmembrane domain, the structure determined in nanodiscs is similar to the previously determined X-ray structure. Our findings show that large differences in structure between detergent and lipid bilayer environments are unlikely, and enable us to propose possible structural mechanisms for C-type inactivation., Competing Interests: DM, CB, GT, TF, AB No competing interests declared, SS, KS Reviewing editor, eLife

We introduce a two-dimensional version of the method called on-the-fly free energy parametrization (OTFP) to reconstruct free-energy surfaces using Molecular Dynamics simulations, which we name OTFP-2D. We first test the new method by reconstructing the well-known dihedral angles free energy surface of solvated alanine dipeptide. Then, we use it to investigate the process of K + ions translocation inside the Kv1.2 channel. By comparing a series of two-dimensional free energy surfaces for ion movement calculated with different conditions on the intercalated water molecules, we first recapitulate the widely accepted knock-on mechanism for ion translocation and then confirm that permeation occurs with water molecules alternated among the ions, in accordance with the latest experimental findings. From a methodological standpoint, our new OTFP-2D algorithm demonstrates the excellent sampling acceleration of temperature-accelerated molecular dynamics and the ability to efficiently compute 2D free-energy surfaces. It will therefore be useful in large variety complex biomacromolecular simulations.

Membrane potential regulates the activity of voltage-dependent ion channels via specialized voltage-sensing modules, but the mechanisms involved in coupling voltage-sensor movement to pore opening remain unclear owing to a lack of resting state structures and robust methods to identify allosteric pathways. Here, using a newly developed interaction-energy analysis, we probe the interfaces of the voltage-sensing and pore modules in the Drosophila Shaker K + channel. Our measurements reveal unexpectedly strong equilibrium gating interactions between contacts at the S4 and S5 helices in addition to those between S6 and the S4-S5 linker. Network analysis of MD trajectories shows that the voltage-sensor and pore motions are linked by two distinct pathways: a canonical pathway through the S4-S5 linker and a hitherto unknown pathway akin to rack-and-pinion coupling involving the S4 and S5 helices. Our findings highlight the central role of the S5 helix in electromechanical transduction in the voltage-gated ion channel (VGIC) superfamily.

We analyze the entropic effects of inner pore geometry changes of Kv 1.2 channel during membrane depolarization and their implications for the rate of transmembrane transport of potassium ions. We base this on the idea that spatial confinements within the channel pore give rise to entropic barriers which can both effectively affect the stability of open macroconformation and influence channel's ability to conduct the potassium ions through the membrane. First, we calculate the differences in entropy between voltage-activated and resting states of the channel. As a template, we take a set of structures of channel pore in an open state at different membrane potentials generated in our previous research. The obtained results indicate that tendency to occupy open states at membrane depolarization is entropy facilitated. Second, we describe the differences in rates of K + transport through the channel pore at different voltages based on the results of appropriate random walk simulations in entropic and electric potentials. The simulated single channel currents (I) suggest that the geometry changes during membrane depolarization are an important factor contributing to the observed flow of potassium ions through the channel. Nevertheless, the charge distribution within the channel pore (especially at the extracellular entrance) seems most prominent for the observed I/I max relation at a qualitative level at analyzed voltages.

The K v 1.2 channel plays an important role in the maintenance of resting membrane potential and the regulation of the cellular excitability of neurons, whose silencing or mutations can elicit neuropathic pain or neurological diseases (e.g., epilepsy and ataxia). Scorpion venom contains a variety of peptide toxins targeting the pore region of this channel. Despite a large amount of structural and functional data currently available, their detailed interaction modes are poorly understood. In this work, we choose four K v 1.2-targeted scorpion toxins (Margatoxin, Agitoxin-2, OsK-1, and Mesomartoxin) to construct their complexes with K v 1.2 based on the experimental structure of ChTx-K v 1.2. Molecular dynamics simulation of these complexes lead to the identification of hydrophobic patches, hydrogen-bonds, and salt bridges as three essential forces mediating the interactions between this channel and the toxins, in which four K v 1.2-specific interacting amino acids (D353, Q358, V381, and T383) are identified for the first time. This discovery might help design highly selective K v 1.2-channel inhibitors by altering amino acids of these toxins binding to the four channel residues. Finally, our results provide new evidence in favor of an induced fit model between scorpion toxins and K⁺ channel interactions., Competing Interests: The authors declare no conflict of interest.

Kv1.2 is a prominent potassium channel subtype in the nervous system and serves as an important structural template for investigation of ion channel function. However, Kv1.2 voltage-dependence exhibits dramatic cell-to-cell variability due to a gating mode shift that is regulated by an unknown mechanism. We report that this variable behavior is regulated by the extracellular redox environment. Exposure to reducing agents promotes a shift in gating properties towards an 'inhibited' gating mode that resists opening, and causes channels to exhibit pronounced use-dependent activation during trains of repetitive depolarizations. This sensitivity to extracellular redox potential is absent in other Kv1 channels, but is apparent in heteromeric channels containing Kv1.2 subunits, and overlaps with the reported physiological range of extracellular redox couples. Mutagenesis of candidate cysteine residues fails to abolish redox sensitivity. Therefore, we suggest that an extrinsic, redox-sensitive binding partner imparts these properties.

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.

We report isolation, sequencing, and electrophysiological characterization of OSK3 (α-KTx 8.8 in Kalium and Uniprot databases), a potassium channel blocker from the scorpion Orthochirus scrobiculosus venom. Using the voltage clamp technique, OSK3 was tested on a wide panel of 11 voltage-gated potassium channels expressed in Xenopus oocytes, and was found to potently inhibit Kv1.2 and Kv1.3 with IC 50 values of ~331nM and ~503nM, respectively. OdK1 produced by the scorpion Odontobuthus doriae differs by just two C-terminal residues from OSK3, but shows marked preference to Kv1.2. Based on the charybdotoxin-potassium channel complex crystal structure, a model was built to explain the role of the variable residues in OdK1 and OSK3 selectivity., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Current through the σ-pore was first detected in hKv1.3_V388C channels, where the V388C mutation in hKv1.3 channels opened a new pathway (σ-pore) behind the central α-pore. Typical for this mutant channel was inward current at potentials more negative than -100 mV when the central α-pore was closed. The α-pore blockers such as TEA+ and peptide toxins (CTX, MTX) could not reduce current through the σ-pore of hKv1.3_V388C channels. This new pathway would proceed in parallel to the α-pore in the S6-S6 interface gap. To see whether this phenomenon is restricted to hKv1.3 channels we mutated hKv1.2 at the homologue position (hKv1.2_V370C). By overexpression of hKv1.2_V370C mutant channels in COS-7 cells we could show typical σ-currents. The electrophysiological properties of the σ-pore in hKv1.3_V388C and hKv1.2_V370C mutant channels were similar. The σ-pore of hKv1.2_V370C channels was most permeable to Na+ and Li+ whereas Cl- and protons did not influence current through the σ-pore. Tetraethylammonium (TEA+), charybdotoxin (CTX) and maurotoxin (MTX), known α-pore blockers, could not reduce current through the σ-pore of hKv1.2_V370C channels. Taken together we conclude that the observation of σ-pore currents is not restricted to Kv1.3 potassium channels but can also be observed in a closely related potassium channel. This finding could have implications in the treatment of different ion channel diseases linked to mutations of the respective channels in regions close to homologue position investigated by us.

Electrical cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the molecular details of voltage sensing remain insufficiently understood. Here, we present a molecular dynamics simulation-based method to determine the underlying charge movement across the membrane-the gating charge-by measuring electrical capacitor properties of membrane-embedded proteins. We illustrate the approach by calculating the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K + channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against experimental data. We then use the gating charge analysis to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na + ions in Na + -coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with experiments, and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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

The mechanisms of ion permeation through potassium channels have been extensively examined. Molecular dynamics (MD) simulations have demonstrated that rapidly permeating ions collide near the selectivity filter (SF) ("knock-on" mechanism), but this oversimplified mechanism is insufficient to account for the experimentally observed single-channel current amplitudes. Here, we analyzed the MD-simulated ion trajectories through a Kv1.2 potassium channel using an event-oriented analysis method, and surprisingly, we found that the nanocavity (NC) governs ion permeation in a digital fashion. The NC has a maximal diameter of 10 Å and stands between the intracellular bulk solution and the SF, which holds only up to one K(+) during permeation. Accordingly, the K(+) concentration in the intracellular solution is translated as a digitalized zero or one K(+) in the NC. When the ion number in the NC is zero, the multiple ions in the SF are mostly immobilized. By contrast, when the number of ions in the NC is one, the structured water in the NC mediates the ion-occupied status to the queueing ions in the SF, and the ions then initiate a collective outward motion. Accordingly, the one ion in the NC serves as a catalytic intermediate for permeation, which quantitatively accounts for the experimentally obtained conductance-concentration relationships. We conclude that the ion movements are coherent across the entire pore.

The mechanism of ion conduction by potassium channels is one of the central issues in physiology. In particular, it is still unclear how the ion concentration and the membrane voltage drive ion conduction. We have investigated the dynamics of the ion conduction processes in the Kv1.2 pore domain, by molecular dynamics (MD) simulations with several different voltages and ion concentrations. By focusing on the detailed ion movements through the pore including selectivity filter (SF) and cavity, we found two major conduction mechanisms, called the III-IV-III and III-II-III mechanisms, and the balance between the ion concentration and the voltage determines the mechanism preference. In the III-IV-III mechanism, the outermost ion in the pore is pushed out by a new ion coming from the intracellular fluid, and four-ion states were transiently observed. In the III-II-III mechanism, the outermost ion is pulled out first, without pushing by incoming ions. Increases in the ion concentration and voltage accelerated ion conductions, but their mechanisms were different. The increase in the ion concentrations facilitated the III-IV-III conductions, while the higher voltages increased the III-II-III conductions, indicating that the pore domain of potassium channels permeates ions by using two different driving forces: a push by intracellular ions and a pull by voltage.

1,1,1-Trichloro-2,2-bis(p-chlorophenyl)ethane (DDT), the first organochlorine insecticide, and pyrethroid insecticides are sodium channel agonists. Although the use of DDT is banned in most of the world due to its detrimental impact on the ecosystem, indoor residual spraying of DDT is still recommended for malaria control in Africa. Development of resistance to DDT and pyrethroids is a serious global obstacle for managing disease vectors. Mapping DDT binding sites is necessary for understanding mechanisms of resistance and modulation of sodium channels by structurally different ligands. The pioneering model of the housefly sodium channel visualized the first receptor for pyrethroids, PyR1, in the II/III domain interface and suggested that DDT binds within PyR1. Previously, we proposed the second pyrethroid receptor, PyR2, at the I/II domain interface. However, whether DDT binds to both pyrethroid receptor sites remains unknown. Here, using computational docking of DDT into the Kv1.2-based mosquito sodium channel model, we predict that two DDT molecules can bind simultaneously within PyR1 and PyR2. The bulky trichloromethyl group of each DDT molecule fits snugly between four helices in the bent domain interface, whereas two p-chlorophenyl rings extend into two wings of the interface. Model-driven mutagenesis and electrophysiological analysis confirmed these propositions and revealed 10 previously unknown DDT-sensing residues within PyR1 and PyR2. Our study proposes a dual DDT-receptor model and provides a structural background for rational development of new insecticides., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

In this study, we report on a series of molecular dynamics simulations that were used to examine the effects of shock waves on a membrane-bound ion channel. A planar shock wave was found to compress the ion channel upon impact, but the protein geometry resembles the crystal structure as soon as the solvent density begins to dissipate. When a void was placed in close proximity to the membrane, the shock wave proved to be more destructive to the protein due to formation of a nanojet that results from the asymmetric collapse of the void. The nanojet was able to cause significant structural changes to the protein even at low piston velocities that are not able to directly cause poration of the membrane., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Inhalational general anesthesia results from the poorly understood interactions of haloethers with multiple protein targets, which prominently includes ion channels in the nervous system. Previously, we reported that the commonly used inhaled anesthetic sevoflurane potentiates the activity of voltage-gated K+ (Kv) channels, specifically, several mammalian Kv1 channels and the Drosophila K-Shaw2 channel. Also, previous work suggested that the S4-S5 linker of K-Shaw2 plays a role in the inhibition of this Kv channel by n-alcohols and inhaled anesthetics. Here, we hypothesized that the S4-S5 linker is also a determinant of the potentiation of Kv1.2 and K-Shaw2 by sevoflurane. Following functional expression of these Kv channels in Xenopus oocytes, we found that converse mutations in Kv1.2 (G329T) and K-Shaw2 (T330G) dramatically enhance and inhibit the potentiation of the corresponding conductances by sevoflurane, respectively. Additionally, Kv1.2-G329T impairs voltage-dependent gating, which suggests that Kv1.2 modulation by sevoflurane is tied to gating in a state-dependent manner. Toward creating a minimal Kv1.2 structural model displaying the putative sevoflurane binding sites, we also found that the positive modulations of Kv1.2 and Kv1.2-G329T by sevoflurane and other general anesthetics are T1-independent. In contrast, the positive sevoflurane modulation of K-Shaw2 is T1-dependent. In silico docking and molecular dynamics-based free-energy calculations suggest that sevoflurane occupies distinct sites near the S4-S5 linker, the pore domain and around the external selectivity filter. We conclude that the positive allosteric modulation of the Kv channels by sevoflurane involves separable processes and multiple sites within regions intimately involved in channel gating.

Voltage-gated potassium channels of the Kv1 family play a crucial role in the generation and transmission of electrical signals in excitable cells affecting neuronal and cardiac activities. Small-molecule blockage of these channels has been proposed to occur via a cooperative mechanism involving two main blocking sites: the inner-pore site located below the selectivity filter, and a side-pocket cavity located between the pore and the voltage sensor. Using 0.5 μs molecular dynamics simulation trajectories complemented by docking calculations, the potential binding sites of the PAP-1 (5-(4-phenoxybutoxy)psoralen) blocker to the crystal structure of Kv1.2 channel have been studied. The presence of both mentioned blocking sites at Kv1.2 is confirmed, adding evidence in favor of a cooperative channel blockage mechanism. These observations provide insight into drug modulation that will guide further developments of Kv inhibitors.

Much has been learned about the voltage sensors of ion channels since the x-ray structure of the mammalian voltage-gated potassium channel Kv1.2 was published in 2005. High resolution structural data of a Kv channel enabled the structural interpretation of numerous electrophysiological findings collected in various ion channels, most notably Shaker, and permitted the development of meticulous computational simulations of the activation mechanism. The fundamental premise for the structural interpretation of functional measurements from Shaker is that this channel and Kv1.2 have the same characteristics, such that correlation of data from both channels would be a trivial task. We tested these assumptions by measuring Kv1.2 voltage-dependent gating and charge per channel. We found that the Kv1.2 gating charge is near 10 elementary charges (eo), ∼25% less than the well-established 13-14 eo in Shaker. Next, we neutralized positive residues in the Kv1.2 S4 transmembrane segment to investigate the cause of the reduction of the gating charge and found that, whereas replacing R1 with glutamine decreased voltage sensitivity to ∼50% of the wild-type channel value, mutation of the subsequent arginines had a much smaller effect. These data are in marked contrast to the effects of charge neutralization in Shaker, where removal of the first four basic residues reduces the gating charge by roughly the same amount. In light of these differences, we propose that the voltage-sensing domains (VSDs) of Kv1.2 and Shaker might undergo the same physical movement, but the septum that separates the aqueous crevices in the VSD of Kv1.2 might be thicker than Shaker's, accounting for the smaller Kv1.2 gating charge., (© 2015 Ishida et al.)

Venom-derived neurotoxins are ideal probes for the investigation of structure-function relationship of ion channels and promising scaffolds for the design of ion channel-targeted drug leads as well. The discovery of highly selective toxins against a specific channel subtype facilitates the development of drugs with reduced side effects. Here, we describe the systemic characterization of a new scorpion short-chain K(+) channel blocker from Mesobuthus martensii, termed mesomartoxin (MMTX). MMTX is synthesized as a precursor comprising a signal peptide and a mature peptide of 29 residues. Nuclear magnetic resonance analysis confirmed that recombinant MMTX adopts a typical cysteine-stabilized α-helical and β-sheet fold. Electrophysiological experiments showed that MMTX exhibits high affinity for the Drosophila Shaker K(+) channel but differential selectivity on different members of the rat voltage-gated K(+) channel (Kv) family, with nanomolar affinity (IC50=15.6 nM) for rKv1.2, micromolar affinity for rKv1.3 (IC50=12.5 μM) and no activity on rKv1.1 at >50 μM. Site-directed mutagenesis of the channel pore identified a key site located on the selectivity filter of the pore, which is directly implicated in toxin binding and controls target's selectivity of the toxin. Given a key role of Kv1.2 in epilepsy, MMTX might serve as a potential drug lead for the disease., (Copyright © 2014 Elsevier Inc. All rights reserved.)

The gating mechanism of voltage sensitive ion channels is generally considered to be the motion of the S4 transmembrane segment of the voltage sensing domains (VSD). The primary supporting evidence came from R → C mutations on the S4 transmembrane segment of the VSD, followed by reaction with a methanethiosulfonate (MTS) reagent. The cys side chain is -SH (reactive form -S-); the arginine side chain is much larger, leaving space big enough to accommodate the MTS sulfonate head group. The cavity created by the mutation has space for up to seven more water molecules than were present in wild type, which could be displaced irreversibly by the MTS reagent. Our quantum calculations show there is major reorientation of three aromatic residues that face into the cavity in response to proton displacement within the VSD. Two phenylalanines reorient sufficiently to shield/unshield the cysteine from the intracellular and extracellular ends, depending on the proton positions, and a tyrosine forms a hydrogen bond to the cysteine sulfur with its side chain -OH. These could produce the results of the experiments that have been interpreted as evidence for physical motion of the S4 segment, without physical motion of the S4 backbone. The computations strongly suggest that the interpretation of cysteine substitution reaction experiments be re-examined in the light of these considerations.

Voltage sensor domains (VSDs) are membrane-bound protein modules that confer voltage sensitivity to membrane proteins. VSDs sense changes in the transmembrane voltage and convert the electrical signal into a conformational change called activation. Activation involves a reorganization of the membrane protein charges that is detected experimentally as transient currents. These so-called gating currents have been investigated extensively within the theoretical framework of so-called discrete-state Markov models (DMMs), whereby activation is conceptualized as a series of transitions across a discrete set of states. Historically, the interpretation of DMM transition rates in terms of transition state theory has been instrumental in shaping our view of the activation process, whose free-energy profile is currently envisioned as composed of a few local minima separated by steep barriers. Here we use atomistic level modeling and well-tempered metadynamics to calculate the configurational free energy along a single transition from first principles. We show that this transition is intrinsically multidimensional and described by a rough free-energy landscape. Remarkably, a coarse-grained description of the system, based on the use of the gating charge as reaction coordinate, reveals a smooth profile with a single barrier, consistent with phenomenological models. Our results bridge the gap between microscopic and macroscopic descriptions of activation dynamics and show that choosing the gating charge as reaction coordinate masks the topological complexity of the network of microstates participating in the transition. Importantly, full characterization of the latter is a prerequisite to rationalize modulation of this process by lipids, toxins, drugs, and genetic mutations.

Many drug molecules inhibit the conduction of several families of cation channels by binding to a small cavity just below the selectivity filter of the channel protein. The exact mechanisms governing drug-channel binding and the subsequent inhibition of conduction are not well understood. Here the inhibition of two K(+) channel isoforms, Kv1.2 and KCa3.1, by two drug molecules, lidocaine and TRAM-34, is examined in atomic detail using molecular dynamics simulations. A conserved valine-alanine-valine motif in the inner cavity is found to be crucial for drug binding in both channels, consistent with previous studies of similar systems. Potential of mean force calculations show that lidocaine in its charged form creates an energy barrier of ∼6 kT for a permeating K(+) ion when the ion is crossing over the drug, while the neutral form of lidocaine has no significant effect on the energetics of ion permeation. On the other hand, TRAM-34 in the neutral form is able to create a large energy barrier of ∼10 kT by causing the permeating ion to dehydrate. Our results suggest that TRAM-34 analogues that remain neutral and permeable to membranes under acidic conditions common to inflammation may act as possible drug scaffolds for combating local anesthetic failure in inflammation.

Signaling proteins such as ion channels largely exist in two functional forms, corresponding to the active and resting states, connected by multiple intermediates. Multiparametric kinetic models based on sophisticated electrophysiological experiments have been devised to identify molecular interactions of these conformational transitions. However, this approach is arduous and is not suitable for large-scale perturbation analysis of interaction pathways. Recently, we described a model-free method to obtain the net free energy of activation in voltage- and ligand-activated ion channels. Here we extend this approach to estimate pairwise interaction energies of side chains that contribute to gating transitions. Our approach, which we call generalized interaction-energy analysis (GIA), combines median voltage estimates obtained from charge-voltage curves with mutant cycle analysis to ascertain the strengths of pairwise interactions. We show that, for a system with an arbitrary gating scheme, the nonadditive contributions of amino acid pairs to the net free energy of activation can be computed in a self-consistent manner. Numerical analyses of sequential and allosteric models of channel activation also show that this approach can measure energetic nonadditivities even when perturbations affect multiple transitions. To demonstrate the experimental application of this method, we reevaluated the interaction energies of six previously described long-range interactors in the Shaker potassium channel. Our approach offers the ability to generate detailed interaction energy maps in voltage- and ligand-activated ion channels and can be extended to any force-driven system as long as associated "displacement" can be measured., (© 2014 Chowdhury et al.)

Today, drug discovery routinely uses experimental assays to determine very early if a lead compound can yield certain types of off-target activity. Among such off targets is hERG. The ion channel plays a primordial role in membrane repolarization and altering its activity can cause severe heart arrhythmia and sudden death. Despite routine tests for hERG activity, rather little information is available for helping medicinal chemists and molecular modelers to rationally circumvent hERG activity. In this article novel insights into the dynamics of hERG channel closure are described. Notably, helical pairwise closure movements have been observed. Implications and relations to hERG inactivation are presented. Based on these dynamics novel insights on hERG blocker placement are presented, compared to literature, and discussed. Last, new evidence for horizontal ligand positioning is shown in light of former studies on hERG blockers.

Voltage sensor domains (VSDs) regulate ion channels and enzymes by transporting electrically charged residues across a hydrophobic VSD constriction called the gating pore or hydrophobic plug. How the gating pore controls the gating charge movement presently remains debated. Here, using saturation mutagenesis and detailed analysis of gating currents from gating pore mutations in the Shaker Kv channel, we identified statistically highly significant correlations between VSD function and physicochemical properties of gating pore residues. A necessary small residue at position S240 in S1 creates a "steric gap" that enables an intracellular access pathway for the transport of the S4 Arg residues. In addition, the stabilization of the depolarized VSD conformation, a hallmark for most Kv channels, requires large side chains at positions F290 in S2 and F244 in S1 acting as "molecular clamps," and a hydrophobic side chain at position I237 in S1 acting as a local intracellular hydrophobic barrier. Finally, both size and hydrophobicity of I287 are important to control the main VSD energy barrier underlying transitions between resting and active states. Taken together, our study emphasizes the contribution of several gating pore residues to catalyze the gating charge transfer. This work paves the way toward understanding physicochemical principles underlying conformational dynamics in voltage sensors.

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Members of the voltage-gated ion channel superfamily (VGIC) regulate ion flux and generate electrical signals in excitable cells by opening and closing pore gates. The location of the gate in voltage-gated sodium channels, a founding member of this superfamily, remains unresolved. Here we explore the chemical modification rates of introduced cysteines along the S6 helix of domain IV in an inactivation-removed background. We find that state-dependent accessibility is demarcated by an S6 hydrophobic residue; substituted cysteines above this site are not modified by charged thiol reagents when the channel is closed. These accessibilities are consistent with those inferred from open- and closed-state structures of prokaryotic sodium channels. Our findings suggest that an intracellular gate composed of a ring of hydrophobic residues is not only responsible for regulating access to the pore of sodium channels, but is also a conserved feature within canonical members of the VGIC superfamily.

Members of the six-transmembrane segment family of ion channels share a common structural design. However, there are sequence differences between the members that confer distinct biophysical properties on individual channels. Currently, we do not have 3D structures for all members of the family to help explain the molecular basis for the differences in their biophysical properties and pharmacology. This is due to low-level expression of many members in native or heterologous systems. One exception is rat Kv1.2 which has been overexpressed in Pichia pastoris and crystallised. Here, we tested chimaeras of rat Kv1.2 with the hERG channel for function in Xenopus oocytes and for overexpression in Pichia. Chimaera containing the S1-S6 transmembrane region of HERG showed functional and pharmacological properties similar to hERG and could be overexpressed and purified from Pichia. Our results demonstrate that rat Kv1.2 could serve as a surrogate to express difficult-to-overexpress members of the six-transmembrane segment channel family.

Quantitative structure-based modeling of voltage activation of ion channels is very challenging. For example, it is very hard to reach converging results, by microscopic simulations while macroscopic treatments involve major uncertainties regarding key features. The current work overcomes some of the above challenges by using our recently developed coarse-grained (CG) model in simulating the activation of the Kv1.2 channel. The CG model has allowed us to explore problems that cannot be fully addressed at present by microscopic simulations, while providing insights on some features that are not usually considered in continuum models, including the distribution of the electrolytes between the membrane and the electrodes during the activation process and thus the physical nature of the gating current. Here, we demonstrate that the CG model yields realistic gating charges and free energy landscapes that allow us to simulate the fluctuating gating current in the activation processes. Our ability to simulate the time dependence of the fast gating current allows us to reproduce the observed trend and provides a clear description of its relationship to the landscape involved in the activation process.

The open gate of the Kv1.2 voltage-gated potassium channel can just hold a hydrated K(+) ion. Quantum calculations starting from the x-ray coordinates of the channel confirm this, showing little change from the x-ray coordinates for the protein. Water molecules not in the x-ray coordinates, and the ion itself, are placed by the calculation. The water molecules, including their orientation and hydrogen bonding, with and without an ion, are critical for the path of the ion, from the solution to the gate. A sequence of steps is postulated in which the potential experienced by the ion in the pore is influenced by the position of the ion. The gate structure, with and without the ion, has been optimized. The charges on the atoms and bond lengths have been calculated using natural bond orbital calculations, giving K(+) ~0.77 charges, rather than 1.0. The PVPV hinge sequence has been mutated in silico to PVVV (P407V in the 2A79 numbering). The water structure around the ion becomes discontinuous, separated into two sections, above and below the ion. PVPV conservation closely relates to maintaining the water structure. Finally, these results have implications concerning gating., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Molecular dynamics simulations is an important application in theoretical chemistry, and with the large high-performance computing resources available today the programs also generate huge amounts of output data. In particular in life sciences, with complex biomolecules such as proteins, simulation projects regularly deal with several terabytes of data. Apart from the need for more cost-efficient storage, it is increasingly important to be able to archive data, secure the integrity against disk or file transfer errors, to provide rapid access, and facilitate exchange of data through open interfaces. There is already a whole range of different formats used, but few if any of them (including our previous ones) fulfill all these goals. To address these shortcomings, we present "Trajectory Next Generation" (TNG)--a flexible but highly optimized and efficient file format designed with interoperability in mind. TNG both provides state-of-the-art multiframe compression as well as a container framework that will make it possible to extend it with new compression algorithms without modifications in programs using it. TNG will be the new file format in the next major release of the GROMACS package, but it has been implemented as a separate library and API with liberal licensing to enable wide adoption both in academic and commercial codes., (Copyright © 2013 Wiley Periodicals, Inc.)

The voltage-gated potassium channel Kv1.3 is an established target for treatment of autoimmune diseases. Hence, there are intense efforts to develop immunosuppressant drugs from Kv1.3 blockers. ShK toxin from sea anemone is the most advanced peptide in this regard, but its lack of selectivity for Kv1.3 over Kv1.1 is an ongoing concern. The scorpion toxin HsTx1 is an equally potent blocker of Kv1.3, which is also selective for Kv1.3. It is of interest to understand the molecular basis of this selectivity as the lessons learned may suggest new avenues for enhancing the selectivity of other Kv1.3 blockers. Here we construct accurate models of Kv1.x-HsTx1 complexes using docking and molecular dynamics simulations. For each complex, the binding free energy of HsTx1 is determined from the potential of mean force calculations. Good agreement is found between the computed and experimental binding free energies, which increases confidence in the complex models. Comparison of the binding modes of HsTx1 with Kv1.1 and Kv1.3 reveals that the lower affinity of HsTx1 for Kv1.1 is due to its inability to come close to the pore domain, which prevents the pore inserting lysine from making proper contacts with the tyrosine carbonyls in the selectivity filter.

A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific β-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.