skip to main content


Title: Thermodynamics of ion binding and occupancy in potassium channels
Potassium channels modulate various cellular functions through efficient and selective conduction of K + ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K + ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.  more » « less
Award ID(s):
1856173
NSF-PAR ID:
10327045
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Chemical Science
Volume:
12
Issue:
25
ISSN:
2041-6520
Page Range / eLocation ID:
8920 to 8930
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Voltage-gated sodium, potassium, and calcium channels consist of four voltage-sensing domains (VSDs) that surround a central pore domain and transition from a down state to an up state in response to membrane depolarization. While many types of drugs bind pore domains, the number of organic molecules known to bind VSDs is limited. The Hv1 voltage-gated proton channel is made of two VSDs and does not contain a pore domain, providing a simplified model for studying how small ligands interact with VSDs. Here, we describe a ligand, named HIF, that interacts with the Hv1 VSD in the up and down states. We find that HIF rapidly inhibits proton conduction in the up state by blocking the open channel, as previously described for 2-guanidinobenzimidazole and its derivatives. HIF, however, interacts with a site slowly accessible in the down state. Functional studies and MD simulations suggest that this interaction traps the compound in a narrow pocket lined with charged residues within the VSD intracellular vestibule, which results in slow recovery from inhibition. Our findings point to a “wrench in gears” mechanism whereby side chains within the binding pocket trap the compound as the teeth of interlocking gears. We propose that the use of screening strategies designed to target binding sites with slow accessibility, similar to the one identified here, could lead to the discovery of new ligands capable of interacting with VSDs of other voltage-gated ion channels in the down state.

     
    more » « less
  2. KCNE3 is a potassium channel accessory transmembrane protein that regulates the function of various voltage-gated potassium channels such as KCNQ1. KCNE3 plays an important role in the recycling of potassium ion by binding with KCNQ1. KCNE3 can be found in the small intestine, colon, and in the human heart. Despite its biological significance, there is little information on the structural dynamics of KCNE3 in native-like membrane environments. Molecular dynamics (MD) simulations are a widely used as a tool to study the conformational dynamics and interactions of proteins with lipid membranes. In this study, we have utilized all-atom molecular dynamics simulations to characterize the molecular motions and the interactions of KCNE3 in a bilayer composed of: a mixture of POPC and POPG lipids (3:1), POPC alone, and DMPC alone. Our MD simulation results suggested that the transmembrane domain (TMD) of KCNE3 is less flexible and more stable when compared to the N- and C-termini of KCNE3 in all three membrane environments. The conformational flexibility of N- and C-termini varies across these three lipid environments. The MD simulation results further suggested that the TMD of KCNE3 spans the membrane width, having residue A69 close to the center of the lipid bilayers and residues S57 and S82 close to the lipid bilayer membrane surfaces. These results are consistent with previous biophysical studies of KCNE3. The outcomes of these MD simulations will help design biophysical experiments and complement the experimental data obtained on KCNE3 to obtain a more detailed understanding of its structural dynamics in the native membrane environment. 
    more » « less
  3. C-type inactivation is a time-dependent process of great physiological significance that is observed in a large class of K+ channels. Experimental and computational studies of the pH-activated KcsA channel show that the functional C-type inactivated state, for this channel, is associated with a structural constriction of the selectivity filter at the level of the central glycine residue in the signature sequence, TTV(G)YGD. The structural constriction is allosterically promoted by the wide opening of the intracellular activation gate. However, whether this is a universal mechanism for C-type inactivation has not been established with certainty because similar constricted structures have not been observed for other K+ channels. Seeking to ascertain the general plausibility of the constricted filter conformation, molecular dynamics simulations of a homology model of the pore domain of the voltage-gated potassium channel Shaker were performed. Simulations performed with an open intracellular gate spontaneously resulted in a stable constricted-like filter conformation, providing a plausible nonconductive state responsible for C-type inactivation in the Shaker channel. While there are broad similarities with the constricted structure of KcsA, the hypothetical constricted-like conformation of Shaker also displays some subtle differences. Interestingly, those are recapitulated by the Shaker-like E71V KcsA mutant, suggesting that the residue at this position along the pore helix plays a pivotal role in determining the C-type inactivation behavior. Free energy landscape calculations show that the conductive-to-constricted transition in Shaker is allosterically controlled by the degree of opening of the intracellular activation gate, as observed with the KcsA channel. The behavior of the classic inactivating W434F Shaker mutant is also characterized from a 10-μs MD simulation, revealing that the selectivity filter spontaneously adopts a nonconductive conformation that is constricted at the level of the second glycine in the signature sequence, TTVGY(G)D.

     
    more » « less
  4. Metal-ion batteries (e.g., lithium and sodium ion batteries) are the promising power sources for portable electronics, electric vehicles, and smart grids. Recent metal-ion batteries with organic liquid electrolytes still suffer from safety issues regarding inflammability and insufficient lifetime.1 As the next generation energy storage devices, all-solid-state batteries (ASSBs) have promising potentials for the improved safety, higher energy density, and longer cycle life than conventional Li-ion batteries.2 The nonflammable solid electrolytes (SEs), where only Li ions are mobile, could prevent battery combustion and explosion since the side reactions that cause safety issues as well as degradation of the battery performance are largely suppressed. However, their practical application is hampered by the high resistance arising at the solid–solid electrode–electrolyte interface (including cathode-electrolyte interface and anode-electrolyte interface).3 Several methods have been introduced to optimize the contact capability as well as the electrochemical/chemical stability between the metal anodes (i.e.: Li and Na) and the SEs, which exhibited decent results in decreasing the charge transfer resistance and broadening the range of the stable energy window (i.e., lowing the chemical potential of metal anode below the highest occupied molecular orbital of the SEs).4 Nevertheless, mitigation for the cathode in ASSB is tardily developed because: (1) the porous structure of the cathode is hard to be infiltrated by SEs;5 (2) SEs would be oxidized and decomposed by the high valence state elements at the surface of the cathode at high state of charge.5 Herein, we demonstrate a universal cathode design strategy to achieve superior contact capability and high electrochemical/chemical stability with SEs. Stereolithography is adopted as a manufacturing technique to realize a hierarchical three-dimensional (HTD) electrode architecture with micro-size channels, which is expected to provide larger contact areas with SEs. Then, the manufactured cathode is sintered at 700 °C in a reducing atmosphere (e.g.: H2) to accomplish the carbonization of the resin, delivering sufficiently high electronic conductivity for the cathode. To avoid the direct exposure of the cathode active materials to the SEs, oxidative chemical vapor deposition technique (oCVD) is leveraged to build conformal and highly conducting poly(3,4-ethylenedioxythiophene) (PEDOT) on the surface of the HTD cathode.6 To demonstrate our design strategy, both NCM811 and Na3V2(PO4)3 is selected as active materials in the HTD cathode, then each cathode is paired with organic (polyacrylonitrile-based) and inorganic (sulfur-based) SEs assembled into two batteries (total four batteries). SEM and TEM reveal the micro-size HTD structure with built-in channels. Featured by the HTD architecture, the intrinsic kinetic and thermodynamic conditions will be enhanced by larger surface contact areas, more active sites, improved infusion and electrolyte ion accessibility, and larger volume expansion capability. Disclosed by X-ray computed tomography, the interface between cathode and SEs in the four modified samples demonstrates higher homogeneity at the interface between the cathode and SEs than that of all other pristine samples. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than that of pristine samples, which confirms a weakened space charge layer by the enhanced contact capability. In addition, through Electron Energy Loss Spectroscopy coupled with Scanning Transmission Electron Microscopy, the preserved interface between HTD cathode and SE is identified; however, the decomposing of the pristine cathode is clearly observed. In addition, Finite element method simulations validate that the diffusion dynamics of lithium ions is favored by HTD structure. Such a demonstrated universal strategy provides a new guideline to engineer cathode electrolyte interface by reconstructing electrode structures that can be applicable to all solid-state batteries in a wide range of chemical conditions. 
    more » « less
  5. Hybrid capacitive deionization (HCDI), which combines a capacitive carbon electrode and a redox active electrode in a single device, has emerged as a promising method for water desalination, enabling higher ion removal capacity than devices containing two carbon electrodes. However, to date, the desalination performance of few redox active materials has been reported. For the first time, we present the electrochemical behavior of manganese oxide nanowires with four different tunnel crystal structures as faradaic electrodes in HCDI cells. Two of these phases are square tunnel structured manganese oxides, α-MnO2 and todorokite-MnO2. The other two phases have novel structures that cross-sectional scanning transmission electron microscopy analysis revealed to have ordered and disordered combinations of structural tunnels with different dimensions. The ion removal performance of the nanowires was evaluated not only in NaCl solution, which is traditionally used in laboratory experiments, but also in KCl and MgCl2 solutions, providing better understanding of the behavior of these materials for desalination of brackish water that contains multiple cation species. High ion removal capacities (as large as 27.8 mg g−1, 44.4 mg g−1, and 43.1 mg g−1 in NaCl, KCl, and MgCl2 solutions, respectively) and high ion removal rates (as large as 0.112 mg g−1 s−1, 0.165 mg g−1 s−1, and 0.164 mg g−1 s−1 in NaCl, KCl, and MgCl2 solutions, respectively) were achieved. By comparing ion removal capacity to structural tunnel size, it was found that smaller tunnels do not favor the removal of cations with larger hydrated radii, and more efficient removal of larger hydrated cations can be achieved by utilizing manganese oxides with larger structural tunnels. Extended HCDI cycling and ex situ X-ray diffraction analysis revealed the excellent stability of the manganese oxide electrodes in repeated ion removal/ion release cycles, and compositional analysis of the electrodes indicated that ion removal is achieved through both surface redox reactions and intercalation of ions into the structural tunnels. This work contributes to the understanding of the behavior of faradaic materials in electrochemical water desalination and elucidates the relationship between the electrode material crystal structure and the ion removal capacity/ion removal rate in various salt solutions. 
    more » « less