Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation

Plasticity of the selectivity filter is essential for permeation in lysosomal TPC2 channels

Two-pore channels are pathophysiologically important Na+- and Ca2+-permeable channels expressed in lysosomes and other acidic organelles. Unlike most other ion channels, their permeability is malleable and ligand-tuned such that when gated by the signaling lipid PI(3,5)P2, they are more Na+-selective than when gated by the Ca2+ mobilizing messenger nicotinic acid adenine dinucleotide phosphate. However, the structural basis that underlies such plasticity and single-channel behavior more generally remains poorly understood. A recent Cryo-electron microscopy (cryo-EM) structure of TPC2 bound to PI(3,5)P2 in a proposed open-channel conformation provided an opportunity to address this via molecular dynamics (MD) simulation. To our surprise, simulations designed to compute conductance through this structure revealed almost no Na+ permeation events even at very high transmembrane voltages. However further MD simulations identified a spontaneous transition to a dramatically different conformation of the selectivity filter that involved expansion and a flip in the orientation of two core asparagine residues. This alternative filter conformation was remarkably stable and allowed Na+ to flow through the channel leading to a conductance estimate that was in very good agreement with direct single-channel measurements. Furthermore, this conformation was more permeable for Na+ over Ca2+. Our results have important ramifications not just for understanding the control of ion selectivity in TPC2 channels but also more broadly in terms of how ion channels discriminate ions.

Multiscale Responsive Kinetic Modeling: Quantifying Biomolecular Reaction Flux under Varying Electrochemical Conditions

Attaining a complete thermodynamic and kinetic characterization for processes involving multiple interconnected rare-event transitions remains a central challenge in molecular biophysics. This challenge is amplified when the process must be understood under a range of reaction conditions. Herein, we present a condition-responsive kinetic modeling framework that can combine the strengths of bottom-up rate quantification from multiscale simulations with top-down solution refinement using experimental data. Although this framework can be applied to any process, we demonstrate its use for electrochemically driven transport through channels and transporters. Using the Cl- /H+ antiporter ClC-ec1 as a model system, we show how robust and predictive kinetic solutions can be obtained when the solution space is grounded by thermodynamic constraints, seeded through multiscale rate quantification, and further refined with experimental data, such as electrophysiology assays. Turning to the Shaker K+ channel, we demonstrate that robust solutions and biophysical insights can also be obtained with sufficient experimental data. This multi-pathway method proves capable of identifying single-pathway dominant mechanisms but also highlights that competing and off-pathway flux is still essential to replicate experimental findings and to describe concentration-dependent channel rectification.

Metabolic heterogeneity in tumor microenvironment - A novel landmark for immunotherapy

The surrounding non-cancer cells and tumor cells that make up the tumor microenvironment (TME) have various metabolic rhythms. TME metabolic heterogeneity is influenced by the intricate network of metabolic control within and between cells. DNA, protein, transport, and microbial levels are important regulators of TME metabolic homeostasis. The effectiveness of immunotherapy is also closely correlated with alterations in TME metabolism. The response of a tumor patient to immunotherapy is influenced by a variety of variables, including intracellular metabolic reprogramming, metabolic interaction between cells, ecological changes within and between tumors, and general dietary preferences. Although immunotherapy and targeted therapy have made great strides, their use in the accurate identification and treatment of tumors still has several limitations. The function of TME metabolic heterogeneity in tumor immunotherapy is summarized in this article. It focuses on how metabolic heterogeneity develops and is regulated as a tumor progresses, the precise molecular mechanisms and potential clinical significance of imbalances in intracellular metabolic homeostasis and intercellular metabolic coupling and interaction, as well as the benefits and drawbacks of targeted metabolism used in conjunction with immunotherapy. This offers insightful knowledge and important implications for individualized tumor patient diagnosis and treatment plans in the future.

Molecular dynamics simulations reveal differences in the conformational stability of FtsZs derived from Staphylococcus aureus and Bacillus subtilis

FtsZ is highly conserved among bacteria and plays an essential role in bacterial cell division. The tense conformation of FtsZ bound to GTP assembles into a straight filament via head-to-tail associations, and then the upper subunit of FtsZ hydrolyzes GTP bound to the lower FtsZ subunit. The subunit with GDP bound disassembles accompanied by a conformational change in the subunit from the tense to relaxed conformation. Although crystal structures of FtsZ derived from several bacterial species have been determined, the conformational change from the relaxed to tense conformation has only been observed in Staphylococcus aureus FtsZ (SaFtsZ). Recent cryo-electron microscopy analyses revealed the three-dimensional reconstruction of the protofilament, in which tense molecules assemble via head-to-tail associations. However, the lower resolution of the protofilament suggested that the flexibility of the FtsZ protomers between the relaxed and tense conformations caused them to form in less-strict alignments. Furthermore, this flexibility may also prevent FtsZs other than SaFtsZ from crystalizing in the tense conformation, suggesting that the flexibility of bacterial FtsZs differs. In this study, molecular dynamics simulations were performed using SaFtsZ and Bacillus subtilis FtsZ in several situations, which suggested that different features of the FtsZs affect their conformational stability.

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

Dynamics accelerate the kinetics of ion diffusion through channels: Continuous-time random walk models beyond the mean field approximation

Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels.

Evaluating sequential and allosteric activation models in IKs channels with mutated voltage sensors

The ion-conducting IKs channel complex, important in cardiac repolarization and arrhythmias, comprises tetramers of KCNQ1 α-subunits along with 1-4 KCNE1 accessory subunits and calmodulin regulatory molecules. The E160R mutation in individual KCNQ1 subunits was used to prevent activation of voltage sensors and allow direct determination of transition rate data from complexes opening with a fixed number of 1, 2, or 4 activatable voltage sensors. Markov models were used to test the suitability of sequential versus allosteric models of IKs activation by comparing simulations with experimental steady-state and transient activation kinetics, voltage-sensor fluorescence from channels with two or four activatable domains, and limiting slope currents at negative potentials. Sequential Hodgkin-Huxley-type models approximately describe IKs currents but cannot explain an activation delay in channels with only one activatable subunit or the hyperpolarizing shift in the conductance-voltage relationship with more activatable voltage sensors. Incorporating two voltage sensor activation steps in sequential models and a concerted step in opening via rates derived from fluorescence measurements improves models but does not resolve fundamental differences with experimental data. Limiting slope current data that show the opening of channels at negative potentials and very low open probability are better simulated using allosteric models of activation with one transition per voltage sensor, which implies that movement of all four sensors is not required for IKs conductance. Tiered allosteric models with two activating transitions per voltage sensor can fully account for IKs current and fluorescence activation kinetics in constructs with different numbers of activatable voltage sensors.

Insight into the K channel's selectivity from binding of K+, Na+ and water to N-methylacetamide

In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.

Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Nature of Overcharging and Charge Inversion in Electrical Double Layers

Understanding overcharging and charge inversion is one of the long-standing challenges in soft matter and biophysics. To study these phenomena, we employ the modified Gaussian renormalized fluctuation theory, which allows for the self-consistent accounting of spatially varying ionic strength as well as the spatial variations in dielectric permittivity and excluded volume effects. The underlying dependence of overcharging on the electrostatic coupling is elucidated by varying the surface charge, counterion valency, and dielectric contrast. Consistent with simulations, three characteristic regimes corresponding to weak, moderate, and strong coupling are identified. Important features like the inversion of zeta potential, crowding, and ionic layering at the surface are successfully captured. For weak coupling, there is no overcharging. In the moderate coupling regime, overcharging increases with the surface charge. Finally, in the strong coupling regime, ionic crowding and saturation in overcharging are observed. Our theory predicts a nonmonotonic dependence of charge inversion on multivalent salt concentration as well as the addition of monovalent salt, in quantitative agreement with experiments.

Light-Augmented Multi-ion Interaction in MXene Membrane for Simultaneous Water Treatment and Osmotic Power Generation

The mixing of wastewater and natural water releases abundant osmotic energy. Harvesting this energy could significantly reduce the energy and economic cost of water treatment, leading to sustainable wastewater treatment technology. Yet, such energy harvesting is highly challenging because it requires a material that is highly permeable to nontoxic ions while rejecting toxic ions in wastewater to reach high power density and prevent environmental pollution. In this work, we demonstrate that a light-augmented biomimetic multi-ion interaction in an MXene membrane can simultaneously realize high permeability of Na+ ions for enhanced osmotic power generation and high selectivity to heavy metal ions up to a ratio of 2050 for wastewater treatment. The Na+ permeability is enhanced by the photothermal effect of the MXene membrane. The transport of heavy metal ions, however, is suppressed because, under angstrom-confinement, heavy metal ions are strongly electrostatically repelled by the increased number of permeating Na+ ions. As a result, the membrane can stably generate osmotic power from simulated industrial wastewater, and the power density can be enhanced by 4 times under light illumination of approximate 1 sun intensity. This work highlights the importance of multi-ion interaction for the transport properties of ionic materials, which remains rarely investigated and poorly understood in previous studies.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Effects of membrane potentials on the electroporation of giant unilamellar vesicles

Living organisms maintain a resting membrane potential, which plays an important role in various biophysical and biological processes. In the context of medical applications, irreversible electroporation (IRE) is a non-thermal and minimally invasive technique that utilizes precisely controlled electric field pulses of micro- to millisecond durations to effectively ablate cancer and tumor cells. Previous studies on IRE-induced rupture of cell-mimetic giant unilamellar vesicles (GUVs) have primarily been conducted in the absence of membrane potentials. In this study, we investigated the electroporation of GUVs, including parameters such as the rate constant of rupture and the probability of rupture, in the presence of various negative membrane potentials. The membranes of GUVs were prepared using lipids and channel forming proteins. As the membrane potential increased from 0 to -90 mV, the rate constant of rupture showed a significant increase from (7.5 ± 1.6)×10-3 to (35.6 ± 5.5)×10-3 s-1. The corresponding probability of rupture also exhibited a notable increase from 0.40 ± 0.05 to 0.68 ± 0.05. To estimate the pore edge tension, the electric tension-dependent logarithm of the rate constant was fitted with the Arrhenius equation for different membrane potentials. The presence of membrane potential did not lead to any significant changes in the pore edge tension. The increase in electroporation is reasonably explained by the decrease in the prepore free energy barrier. The choice of buffer used in GUVs can significantly influence the kinetics of electroporation. This study provides valuable insights that can contribute to the application of electroporation techniques in the biomedical field.

Joule heating involving ion currents through channel proteins

Ion currents associated with channel proteins in the presence of membrane potential are ubiquitous in cellular and organelle membranes. When an ion current occurs through a channel protein, Joule heating should occur. However, this Joule heating seems to have been largely overlooked in biology. Here we show theoretical investigation of Joule heating involving channel proteins in biological processes. We used electrochemical potential to derive the Joule's law for an ion current through an ion transport protein in the presence of membrane potential, and we suggest that heat production and absorption can occur. Simulation of temperature distribution around a single channel protein with the Joule heating revealed that the temperature increase was as small as <10-3 K, although an ensemble of channel proteins was suggested to exhibit a noticeable temperature increase. Thereby, we theoretically investigated the Joule heating of systems containing ensembles of channel proteins. Nerve is known to undergo rapid heat production followed by heat absorption during the action potential, and our simulation of Joule heating for a squid giant axon combined with the Hodgkin-Huxley model successfully reproduced the feature of the heat. Furthermore, we extended the theory of Joule heating to uncoupling protein 1 (UCP1), a solute carrier family transporter, which is important to the non-shivering thermogenesis in brown adipose tissue mitochondria (BATM). Our calculations showed that the Joule heat involving UCP1 was comparable to the literature calorimetry data of BATM. Joule heating of ion transport proteins is likely to be one of important mechanisms of cellular thermogenesis.

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.

Elucidating Molecular Mechanisms of Protoxin-2 State-specific Binding to the Human NaV1.7 Channel

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Mechanism of external K+ sensitivity of KCNQ1 channels

KCNQ1 voltage-gated K+ channels are involved in a wide variety of fundamental physiological processes and exhibit the unique feature of being markedly inhibited by external K+. Despite the potential role of this regulatory mechanism in distinct physiological and pathological processes, its exact underpinnings are not well understood. In this study, using extensive mutagenesis, molecular dynamics simulations, and single-channel recordings, we delineate the molecular mechanism of KCNQ1 modulation by external K+. First, we demonstrate the involvement of the selectivity filter in the external K+ sensitivity of the channel. Then, we show that external K+ binds to the vacant outermost ion coordination site of the selectivity filter inducing a diminution in the unitary conductance of the channel. The larger reduction in the unitary conductance compared to whole-cell currents suggests an additional modulatory effect of external K+ on the channel. Further, we show that the external K+ sensitivity of the heteromeric KCNQ1/KCNE complexes depends on the type of associated KCNE subunits.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Multidisciplinary Advances Address the Challenges in Developing Drugs against Transient Receptor Potential Channels to Treat Metabolic Disorders

Transient receptor potential (TRP) channels are cation channels that regulate key physiological and pathological processes in response to a broad range of stimuli. Moreover, they systemically regulate the release of hormones, metabolic homeostasis, and complications of diabetes, which positions them as promising therapeutic targets to combat metabolic disorders. Nevertheless, there are significant challenges in the design of TRP ligands with high potency and durability. Herein we summarize the four challenges as hydrophobicity, selectivity, mono-target therapy, and interspecies discrepancy. We present 1134 TRP ligands with diversified modes of TRP-ligand interaction and provide a detailed discussion of the latest strategies, especially cryogenic electron microscopy (cryo-EM) and computational methods. We propose solutions to address the challenges with a critical analysis of advances in membrane partitioning, polypharmacology, biased agonism, and biochemical screening of transcriptional modulators. They are fueled by the breakthrough from cryo-EM, chemoinformatics and bioinformatics. The discussion is aimed to shed new light on designing next-generation drugs to treat obesity, diabetes and its complications, with optimal hydrophobicity, higher mode selectivity, multi-targeting and consistent activities between human and rodents.

Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca²⁺. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca²⁺ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca²⁺ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca²⁺ plays a pivotal role.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

Crown Ether-Based Ion Transporters in Bilayer Membranes

Bilayer membranes that enhance the stability of the cell are essential for cell survival, separating and protecting the interior of the cell from its external environment. Membrane-based channel proteins are crucial for sustaining cellular activities. However, dysfunction of these proteins would induce serial channelopathies, which could be substituted by artificial ion channel analogs. Crown ethers (CEs) are widely studied in the area of artificial ion channels owing to their intrinsic host-guest interaction with different kinds of organic and inorganic ions. Other advantages such as lower price, chemical stability, and easier modification also make CE a research hotspot in the field of synthetic transmembrane nanopores. And numerous CEs-based membrane-active synthetic ion channels were designed and fabricated in the past decades. Herein, the recent progress of CEs-based synthetic ion transporters has been comprehensively summarized in this review, including their design principles, functional mechanisms, controllable properties, and biomedical applications. Furthermore, this review has been concluded by discussing the future opportunities and challenges facing this research field. It is anticipated that this review could offer some inspiration for the future fabrication of novel CEs-derived ion transporters with more advanced structures, properties, and practical applications.

Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules

Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.

Structure of POPC Lipid Bilayers in OPLS3e Force Field

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters S_CH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured S_CH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca²⁺. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca²⁺ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel

Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K⁺ and Cl^- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K⁺, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl^- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Single ion free energy calculation in ASIC1: the importance of the HG loop

Acid Sensing Ion Channels (ASICs) are one of the most studied channels of the Epithelial Sodium Channel/Degenerin (ENaC/DEG) superfamily. They are responsible for excitatory responses following acidification of the extracellular medium and are involved in several important physiological roles. The ASIC1 subunit can form a functional homotrimeric channel and its structure is currently the most characterised of the whole ENaC/DEG family. Here we computed the free energy profiles for single ion permeation in two different structures of ASIC1 using both Na⁺ and Cl^- as permeating ions. The first structure is the open structure of the channel from the PDB entry 4NTW, and the second structure is the closed structure with the re-entrant loop which contains the highly conserved 'HG' motif form PDB entry 6VTK. Both structures show cation selective free energy profiles, however the profiles of the permeating Na⁺ differ significantly between the two structures. Indeed, whereas there is only a small energetically favorable (-0.5 kcal mol^-1) location for Na⁺ in the open channel (4NTW) near the end of the pore, we observed a clear ion binding site (-7.8 kcal mol^-1) located in between the 'GAS' belt and the 'HG' loop for the channel containing the re-entrant loop (6VTK). Knowing that the 'GAS' motif was determined as the selectivity filter, our results support previous observations while addressing the importance of the 'HG' motif for the interactions between the pore and the permeating cations.

Molecular Dynamics and Theratyping in Airway and Gut Organoids Reveal R352Q-CFTR Conductance Defect

A significant challenge to making targeted cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapies accessible to all individuals with cystic fibrosis (CF) are many mutations in the CFTR gene that can cause CF, most of which remain uncharacterized. Here, we characterized the structural and functional defects of the rare CFTR mutation R352Q, with a potential role contributing to intrapore chloride ion permeation, in patient-derived cell models of the airway and gut. CFTR function in differentiated nasal epithelial cultures and matched intestinal organoids was assessed using an ion transport assay and forskolin-induced swelling assay, respectively. CFTR potentiators (VX-770, GLPG1837, and VX-445) and correctors (VX-809, VX-445, with or without VX-661) were tested. Data from R352Q-CFTR were compared with data of 20 participants with mutations with known impact on CFTR function. R352Q-CFTR has residual CFTR function that was restored to functional CFTR activity by CFTR potentiators but not the corrector. Molecular dynamics simulations of R352Q-CFTR were carried out, which indicated the presence of a chloride conductance defect, with little evidence supporting a gating defect. The combination approach of in vitro patient-derived cell models and in silico molecular dynamics simulations to characterize rare CFTR mutations can improve the specificity and sensitivity of modulator response predictions and aid in their translational use for CF precision medicine.

Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond

Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.

Unifying Single-Channel Permeability From Rare-Event Sampling and Steady-State Flux

Various all-atom molecular dynamics (MD) simulation methods have been developed to compute free energies and crossing rates of ions and small molecules through ion channels. However, a systemic comparison across different methods is scarce. Using a carbon nanotube as a model of small conductance ion channel, we computed the single-channel permeability for potassium ion using umbrella sampling, Markovian milestoning, and steady-state flux under applied voltage. We show that a slightly modified inhomogeneous solubility-diffusion equation yields a single-channel permeability consistent with the mean first passage time (MFPT) based method. For milestoning, applying cylindrical and spherical bulk boundary conditions yield consistent MFPT if factoring in the effective bulk concentration. The sensitivity of the MFPT to the output frequency of collective variables is highlighted using the convergence and symmetricity of the inward and outward MFPT profiles. The consistent transport kinetic results from all three methods demonstrated the robustness of MD-based methods in computing ion channel permeation. The advantages and disadvantages of each technique are discussed, focusing on the future applications of milestoning in more complex systems.

Ion permeation, selectivity, and electronic polarization in fluoride channels

Fluoride channels (Flucs) export toxic F^- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein, we have applied constant-pH molecular dynamics and free-energy-sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F^-, by traversing through a non-bonded network. The extraordinary F^- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F^- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared with the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F^--Fluc interactions.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Predicting Ion Channel Conductance via Dissipation-Corrected Targeted Molecular Dynamics and Langevin Equation Simulations

Ion channels are important proteins for physiological information transfer and functional control. To predict the microscopic origins of their voltage-conductance characteristics, here we applied dissipation-corrected targeted molecular dynamics in combination with Langevin equation simulations to potassium diffusion through the gramicidin A channel as a test system. Performing a nonequilibrium principal component analysis on backbone dihedral angles, we find coupled protein-ion dynamics to occur during ion transfer. The dissipation-corrected free energy profiles correspond well to predictions from other biased simulation methods. The incorporation of an external electric field in Langevin simulations enables the prediction of macroscopic observables in the form of I-V characteristics.

Electropumping Phenomenon in Modified Carbon Nanotubes

Controlling the water transport in a given direction is essential to the design of novel nanofluidic devices, which is still a challenge because of thermal fluctuations on the nanoscale. In this work, we find an interesting electropumping phenomenon for charge-modified carbon nanotubes (CNTs) through a series of molecular dynamics simulations. In electric fields, the flowing counterions on the CNT inner surface provide a direct driving force for water conduction. Specifically, the dynamics of cations and anions exhibit distinct behaviors that lead to thoroughly different water dynamics in positively and negatively charged CNTs. Because of the competition between the increased ion number and ion-CNT interaction, the cation flux displays an interesting maximum behavior with the increase in surface charge density; however, the anion flux rises further at higher charge density because it is less attractive to the surface. Thus, the anion flux is always several times larger than cation flux that induces a higher water flux in positive CNTs with nearly 100% pumping efficiency, which highly exceeds the efficiency of pristine CNTs. With the change in charge density, the translocation time, occupancy number, and radial density profiles for water and ions also demonstrate a nontrivial difference for positive and negative CNTs. Furthermore, the ion flux exhibits an excellent linear relationship with the field strength, leading to the same water flux behavior. For the change in salt concentration, the pumping efficiency for positive CNTs is also nearly 100%. Our results provide significant new insight into the ionic transport through modified CNTs and should be helpful for the design of nanometer water pumps.

Ion Selectivity in the ENaC/DEG Family: A Systematic Review with Supporting Analysis

The Epithelial Sodium Channel/Degenerin (ENaC/DEG) family is a superfamily of sodium-selective channels that play diverse and important physiological roles in a wide variety of animal species. Despite their differences, they share a high homology in the pore region in which the ion discrimination takes place. Although ion selectivity has been studied for decades, the mechanisms underlying this selectivity for trimeric channels, and particularly for the ENaC/DEG family, are still poorly understood. This systematic review follows PRISMA guidelines and aims to determine the main components that govern ion selectivity in the ENaC/DEG family. In total, 27 papers from three online databases were included according to specific exclusion and inclusion criteria. It was found that the G/SxS selectivity filter (glycine/serine, non-conserved residue, serine) and other well conserved residues play a crucial role in ion selectivity. Depending on the ion type, residues with different properties are involved in ion permeability. For lithium against sodium, aromatic residues upstream of the selectivity filter seem to be important, whereas for sodium against potassium, negatively charged residues downstream of the selectivity filter seem to be important. This review provides new perspectives for further studies to unravel the mechanisms of ion selectivity.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Conductance selectivity of Na+ across the K+ channel via Na+ trapped in a tortuous trajectory

Ion selectivity of the potassium channel is crucial for regulating electrical activity in living cells; however, the mechanism underlying the potassium channel selectivity that favors large K⁺ over small Na⁺ remains unclear. Generally, Na⁺ is not completely excluded from permeation through potassium channels. Herein, the distinct nature of Na⁺ conduction through the prototypical KcsA potassium channel was examined. Single-channel current recordings revealed that, at a high Na⁺ concentration (200 mM), the channel was blocked by Na⁺, and this blocking was relieved at high membrane potentials, suggesting the passage of Na⁺ across the channel. At a 2,000 mM Na⁺ concentration, single-channel Na⁺ conductance was measured as one-eightieth of the K⁺ conductance, indicating that the selectivity filter allows substantial conduit of Na⁺ Molecular dynamics simulations revealed unprecedented atomic trajectories of Na⁺ permeation. In the selectivity filter having a series of carbonyl oxygen rings, a smaller Na⁺ was distributed off-center in eight carbonyl oxygen-coordinated sites as well as on-center in four carbonyl oxygen-coordinated sites. This amphipathic nature of Na⁺ coordination yielded a continuous but tortuous path along the filter. Trapping of Na⁺ in many deep free energy wells in the filter caused slow elution. Conversely, K⁺ is conducted via a straight path, and as the number of occupied K⁺ ions increased to three, the concerted conduction was accelerated dramatically, generating the conductance selectivity ratio of up to 80. The selectivity filter allows accommodation of different ion species, but the ion coordination and interactions between ions render contrast conduction rates, constituting the potassium channel conductance selectivity.

Independent Markov decomposition: Toward modeling kinetics of biomolecular complexes

To advance the mission of in silico cell biology, modeling the interactions of large and complex biological systems becomes increasingly relevant. The combination of molecular dynamics (MD) simulations and Markov state models (MSMs) has enabled the construction of simplified models of molecular kinetics on long timescales. Despite its success, this approach is inherently limited by the size of the molecular system. With increasing size of macromolecular complexes, the number of independent or weakly coupled subsystems increases, and the number of global system states increases exponentially, making the sampling of all distinct global states unfeasible. In this work, we present a technique called independent Markov decomposition (IMD) that leverages weak coupling between subsystems to compute a global kinetic model without requiring the sampling of all combinatorial states of subsystems. We give a theoretical basis for IMD and propose an approach for finding and validating such a decomposition. Using empirical few-state MSMs of ion channel models that are well established in electrophysiology, we demonstrate that IMD models can reproduce experimental conductance measurements with a major reduction in sampling compared with a standard MSM approach. We further show how to find the optimal partition of all-atom protein simulations into weakly coupled subunits.

Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential

Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C₄ term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M²⁺) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/E_b, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M²⁺ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.

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.

Polarization shapes the energy landscape of ion conduction in potassium channels.

Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore

A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.

Determinants of ion selectivity in ASIC1a- and ASIC2a-containing acid-sensing ion channels

Trimeric acid-sensing ion channels (ASICs) contribute to neuronal signaling by converting extracellular acidification into excitatory sodium currents. Previous work with homomeric ASIC1a implicates conserved leucine (L7') and consecutive glycine-alanine-serine (GAS belt) residues near the middle, and conserved negatively charged (E18') residues at the bottom of the pore in ion permeation and/or selectivity. However, a conserved mechanism of ion selectivity throughout the ASIC family has not been established. We therefore explored the molecular determinants of ion selectivity in heteromeric ASIC1a/ASIC2a and homomeric ASIC2a channels using site-directed mutagenesis, electrophysiology, and molecular dynamics free energy simulations. Similar to ASIC1a, E18' residues create an energetic preference for sodium ions at the lower end of the pore in ASIC2a-containing channels. However, and in contrast to ASIC1a homomers, ion permeation through ASIC2a-containing channels is not determined by L7' side chains in the upper part of the channel. This may be, in part, due to ASIC2a-specific negatively charged residues (E59 and E62) that lower the energy of ions in the upper pore, thus making the GAS belt more important for selectivity. This is confirmed by experiments showing that the L7'A mutation has no effect in ASIC2a, in contrast to ASIC1a, where it eliminated selectivity. ASIC2a triple mutants eliminating both L7' and upper charges did not lead to large changes in selectivity, suggesting a different role for L7' in ASIC2a compared with ASIC1a channels. In contrast, we observed measurable changes in ion selectivity in ASIC2a-containing channels with GAS belt mutations. Our results suggest that ion conduction and selectivity in the upper part of the ASIC pore may differ between subtypes, whereas the essential role of E18' in ion selectivity is conserved. Furthermore, we demonstrate that heteromeric channels containing mutations in only one of two ASIC subtypes provide a means of functionally testing mutations that render homomeric channels nonfunctional.

Molecular Insights Into the Gating Kinetics of the Cardiac hERG Channel, Illuminated by Structure and Molecular Dynamics

The rapidly activating delayed rectifier K⁺ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.

Potassium and sodium ion complexes with a partial peptide of the selectivity filter in K+ channels studied by cold ion trap infrared spectroscopy: the effect of hydration

Potassium channels allow K+ to rapidly diffuse, while the selectivity filter (SF) actively blocks Na+. The presence of water in the SF during ion translocation remains under debate due to the experimental and computational challenges in characterizing the interactions between water, ions, and the SF. Our bottom-up approach has been applied to a system composed of a partial peptide of the SF (Ac-tyrosine-NHMe) with a metal ion and a single water molecule to probe these interactions. The IR photodissociation spectra of M+Ac-tyrosine-NHMe(H2O) (M = Na, K) combined with quantum chemical calculations revealed that the water molecule binding sites are ion-dependent. In addition, the ion-peptide distances are elongated significantly for the K+ complex in comparison to the Na+ complex by the addition of a single water molecule. This striking structural difference with the water molecule is discussed in relation to ion selectivity and translocation within the K+ channel.

Electrophysiological Properties from Computations at a Single Voltage: Testing Theory with Stochastic Simulations

We use stochastic simulations to investigate the performance of two recently developed methods for calculating the free energy profiles of ion channels and their electrophysiological properties, such as current–voltage dependence and reversal potential, from molecular dynamics simulations at a single applied voltage. These methods require neither knowledge of the diffusivity nor simulations at multiple voltages, which greatly reduces the computational effort required to probe the electrophysiological properties of ion channels. They can be used to determine the free energy profiles from either forward or backward one-sided properties of ions in the channel, such as ion fluxes, density profiles, committor probabilities, or from their two-sided combination. By generating large sets of stochastic trajectories, which are individually designed to mimic the molecular dynamics crossing statistics of models of channels of trichotoxin, p7 from hepatitis C and a bacterial homolog of the pentameric ligand-gated ion channel, GLIC, we find that the free energy profiles obtained from stochastic simulations corresponding to molecular dynamics simulations of even a modest length are burdened with statistical errors of only 0.3 kcal/mol. Even with many crossing events, applying two-sided formulas substantially reduces statistical errors compared to one-sided formulas. With a properly chosen reference voltage, the current–voltage curves can be reproduced with good accuracy from simulations at a single voltage in a range extending for over 200 mV. If possible, the reference voltages should be chosen not simply to drive a large current in one direction, but to observe crossing events in both directions.

Asymmetric CorA Gating Mechanism as Observed by Molecular Dynamics Simulations

The CorA family of proteins plays a housekeeping role in the homeostasis of divalent metal ions in many bacteria and archaea as well as in mitochondria of eukaryotes, rendering it an important target to study the mechanisms of divalent transport and regulation across different life domains. Despite numerous studies, the mechanistic details of the channel gating and the transport of the metal ions are still not entirely understood. Here, we use all-atom and coarse-grained molecular dynamics simulations combined with in vitro experiments to investigate the influence of divalent cations on the function of CorA. Simulations reveal pronounced asymmetric movements of monomers that enable the rotation of the α7 helix and the cytoplasmic subdomain with the subsequent formation of new interactions and the opening of the channel. These computational results are functionally validated using site-directed mutagenesis of the intracellular cytoplasmic domain residues and biochemical assays. The obtained results infer a complex network of interactions altering the structure of CorA to allow gating. Furthermore, we attempt to reconcile the existing gating hypotheses for CorA to conclude the mechanism of transport of divalent cations via these proteins.