skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


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
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. ; ; ; ; ; ; ; ; ; ; et al (, Membranes)
    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
  2. ; ; ; ; ; ; ; ; ; (, Nature Materials)
    Vincent Dusastre (Ed.)
    Alternative solid-electrolytes are the next key step in advancing lithium batteries with better thermal and chemical stability. A soft-solid electrolyte (Adpn)2LiPF6 (Adpn = adiponitrile) is synthesized and characterized, which exhibits high thermal and electrochemical stability and good ionic conductivity, overcoming several limitations of conventional organic and ceramic materials. The surface of the electrolyte possesses a liquid nano-layer of Adpn that links grains for a facile ionic conduction without high pressure/temperature treatments. Further, the material can quickly self-heal if fractured and provides liquid-like conduction paths via the grain boundaries. A significantly high ion conductivity (~ 10-4 S/cm) and lithium-ion transference number (0.54) are obtained due to weak interactions between “hard” (charge-dense) Li+ ions and “soft” (electronically polarizable) -C≡N group of Adpn. Molecular simulations predict that Li+ ions migrate at the co-crystal grain boundaries with a (preferentially) lower Ea and within the interstitial regions between the co-crystals with higher Ea, where the bulk conductivity comprises a smaller but extant contribution. These cocrystals establish a special concept of crystal design to increase the thermal stability of LiPF6 by separating ions in Adpn solvent matrix, and also exhibit a unique mechanism of ion-conduction via low-resistance grain-boundaries, which is contrasting to ceramics or gel-electrolytes. 
    more » « less
  3. ; ; ; (, Archives of Biochemistry and Biophysics)
    Antivirulence strategy has been explored as an alternative to traditional antibiotic development. The bacterial type IV pilus is a virulence factor involved in host invasion and colonization in many antibiotic resistant pathogens. The PilB ATPase hydrolyzes ATP to drive the assembly of the pilus filament from pilin subunits. We evaluated Chloracidobacterium thermophilum PilB (CtPilB) as a model for structure-based virtual screening by molecular docking and molecular dynamics (MD) simulations. A hexameric structure of CtPilB was generated through homology modeling based on an existing crystal structure of a PilB from Geobacter metallireducens. Four representative structures were obtained from molecular dynamics simulations to examine the conformational plasticity of PilB and improve docking analyses by ensemble docking. Structural analyses after 1 μs of simulation revealed conformational changes in individual PilB subunits are dependent on ligand presence. Further, ensemble virtual screening of a library of 4234 compounds retrieved from the ZINC15 database identified five promising PilB inhibitors. Molecular docking and binding analyses using the four representative structures from MD simulations revealed that top-ranked compounds interact with multiple Walker A residues, one Asp-box residue, and one arginine finger, indicating these are key residues in inhibitor binding within the ATP binding pocket. The use of multiple conformations in molecular screening can provide greater insight into compound flexibility within receptor sites and better inform future drug development for therapeutics targeting the type IV pilus assembly ATPase. 
    more » « less
  4. ; ; ; ; (, European Journal of Inorganic Chemistry)
    Abstract Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu3+ion complexed with EDTA in all its protonation states in aqueous solution. Eu‐EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X‐ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu3+aqua ion and Eu‐EDTA complexes at pH 3 and 11. This work details how Eu‐EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu3+ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation. 
    more » « less
  5. ; ; ; ; ; ; ; (, Molecules)
    Although the 3D structures of active and inactive cannabinoid receptors type 2 (CB2) are available, neither the X-ray crystal nor the cryo-EM structure of CB2-orthosteric ligand-modulator has been resolved, prohibiting the drug discovery and development of CB2 allosteric modulators (AMs). In the present work, we mainly focused on investigating the potential allosteric binding site(s) of CB2. We applied different algorithms or tools to predict the potential allosteric binding sites of CB2 with the existing agonists. Seven potential allosteric sites can be observed for either CB2-CP55940 or CB2-WIN 55,212-2 complex, among which sites B, C, G and K are supported by the reported 3D structures of Class A GPCRs coupled with AMs. Applying our novel algorithm toolset-MCCS, we docked three known AMs of CB2 including Ec2la (C-2), trans-β-caryophyllene (TBC) and cannabidiol (CBD) to each site for further comparisons and quantified the potential binding residues in each allosteric binding site. Sequentially, we selected the most promising binding pose of C-2 in five allosteric sites to conduct the molecular dynamics (MD) simulations. Based on the results of docking studies and MD simulations, we suggest that site H is the most promising allosteric binding site. We plan to conduct bio-assay validations in the future. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email