Ion selectivity in channels and transporters

Evaluating Polarizable Biomembrane Simulations against Experiments

Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.

Bio-Inspired Subnanofluidics: Advanced Fabrication and Functionalization

Biological ion channels can realize high-speed and high-selective ion transport through the protein filter with the sub-1-nanometer channel. Inspired by biological ion channels, various kinds of artificial subnanopores, subnanochannels, and subnanoslits with improved ion selectivity and permeability are recently developed for efficient separation, energy conversion, and biosensing. This review article discusses the advanced fabrication and functionalization methods for constructing subnanofluidic pores, channels, tubes, and slits, which have shown great potential for various applications. Novel fabrication methods for producing subnanofluidics, including top-down techniques such as electron beam etching, ion irradiation, and electrochemical etching, as well as bottom-up approaches starting from advanced microporous frameworks, microporous polymers, lipid bilayer embedded subnanochannels, and stacked 2D materials are well summarized. Meanwhile, the functionalization methods of subnanochannels are discussed based on the introduction of functional groups, which are classified into direct synthesis, covalent bond modifications, and functional molecule fillings. These methods have enabled the construction of subnanochannels with precise control of structure, size, and functionality. The current progress, challenges, and future directions in the field of subnanofluidic are also discussed.

The Molecular Mechanism of Ion Selectivity in Nanopores

Ion channels exhibit strong selectivity for specific ions over others under electrochemical potentials, such as KcsA for K+ over Na+. Based on the thermodynamic analysis, this study is focused on exploring the mechanism of ion selectivity in nanopores. It is well known that ions must lose part of their hydration layer to enter the channel. Therefore, the ion selectivity of a channel is due to the rearrangement of water molecules when entering the nanopore, which may be related to the hydrophobic interactions between ions and channels. In our recent works on hydrophobic interactions, with reference to the critical radius of solute (Rc), it was divided into initial and hydrophobic solvation processes. Additionally, the different dissolved behaviors of solutes in water are expected in various processes, such as dispersed and accumulated distributions in water. Correspondingly, as the ion approaches the nanopore, there seems to exist the "repulsive" or "attractive" forces between them. In the initial process ( Collapse

Key Words

• hydrogen bondings
• hydrophobic interactions
• ion selectivity
• nanopore
• water

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Yan-Nan Chen
Yu-Zhen Liu
Qiang Sun Key Laboratory of Orogenic Belts and Crustal Evolution, The School of Earth and Space Sciences, Ministry of Education, Peking University, Beijing 100871, China; (Y.-N.C.); (Y.-Z.L.)

Collapse

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.

Ion Transport and Inhibitor Binding by Human NHE1: Insights from Molecular Dynamics Simulations and Free Energy Calculations

The human Na+/H+ exchanger (NHE1) plays a crucial role in maintaining intracellular pH by regulating the electroneutral exchange of a single intracellular H+ for one extracellular Na+ across the plasma membrane. Understanding the molecular mechanisms governing ion transport and the binding of inhibitors is of importance in the development of anticancer therapeutics targeting NHE1. In this context, we performed molecular dynamics (MD) simulations based on the recent cryo-electron microscopy (cryo-EM) structures of outward- and inward-facing conformations of NHE1. These simulations allowed us to explore the dynamics of the protein, examine the ion-translocation pore, and confirm that Asp267 is the ion-binding residue. Our free energy calculations did not show a significant difference between Na+ and K+ binding at the ion-binding site. Consequently, Na+ over K+ selectivity cannot be solely explained by differences in ion binding. Our MD simulations involving NHE1 inhibitors (cariporide and amiloride analogues) maintained stable interactions with Asp267 and Glu346. Our study highlights the importance of the salt bridge between the positively charged acylguanidine moiety and Asp267, which appears to play a role in the competitive inhibitory mechanism for this class of inhibitors. Our computational study provides a detailed mechanistic interpretation of experimental data and serves the basis of future structure-based inhibitor design.

Binding Behavior between Transforming-Growth-Factor-Beta1 and Its Receptor Reconstituted in Biomimetic Membranes

Transforming growth factor β1 (TGF-β1) is critical to cell differentiation, proliferation, and apoptosis. It is important to understand the binding affinity between TGF-β1 and its receptors. In this study, their binding force was measured using an atomic force microscope. Significant adhesion was induced by the interaction between the TGF-β1 immobilized on the tip and its receptor reconstituted in the bilayer. Rupture and adhesive failure occurred at a specific force around 0.4~0.5 nN. The relationship of the force to loading rate was used to estimate the displacement where the rupture occurred. The binding was also monitored in real time with surface plasmon resonance (SPR) and interpreted with kinetics to acquire the rate constant. Using the Langmuir adsorption, the SPR data were analyzed to estimate equilibrium and association constants to be approximately 10⁷ M^-1 and 10⁶ M^-1 s^-1. These results indicated that the natural release of the binding seldom occurred. Furthermore, the degree of binding dissociation, confirmed by the rupture interpretation, supported that the reverse of the binding hardly happened.

Structure-Based Function and Regulation of NCX Variants: Updates and Challenges

The plasma-membrane homeostasis Na⁺/Ca²⁺ exchangers (NCXs) mediate Ca²⁺ extrusion/entry to dynamically shape Ca²⁺ signaling/in biological systems ranging from bacteria to humans. The NCX gene orthologs, isoforms, and their splice variants are expressed in a tissue-specific manner and exhibit nearly 10⁴-fold differences in the transport rates and regulatory specificities to match the cell-specific requirements. Selective pharmacological targeting of NCX variants could benefit many clinical applications, although this intervention remains challenging, mainly because a full-size structure of eukaryotic NCX is unavailable. The crystal structure of the archaeal NCX_Mj, in conjunction with biophysical, computational, and functional analyses, provided a breakthrough in resolving the ion transport mechanisms. However, NCX_Mj (whose size is nearly three times smaller than that of mammalian NCXs) cannot serve as a structure-dynamic model for imitating high transport rates and regulatory modules possessed by eukaryotic NCXs. The crystal structures of isolated regulatory domains (obtained from eukaryotic NCXs) and their biophysical analyses by SAXS, NMR, FRET, and HDX-MS approaches revealed structure-based variances of regulatory modules. Despite these achievements, it remains unclear how multi-domain interactions can decode and integrate diverse allosteric signals, thereby yielding distinct regulatory outcomes in a given ortholog/isoform/splice variant. This article summarizes the relevant issues from the perspective of future developments.

Bioinspired Artificial Ion Pumps

Ion pumps are important membrane-spanning transporters that pump ions against the electrochemical gradient across the cell membrane. In biological systems, ion pumping is essential to maintain intracellular osmotic pressure, to respond to external stimuli, and to regulate physiological activities by consuming adenosine triphosphate. In recent decades, artificial ion pumping systems with diverse geometric structures and functions have been developing rapidly with the progress of advanced materials and nanotechnology. In this Review, bioinspired artificial ion pumps, including four categories: asymmetric structure-driven ion pumps, pH gradient-driven ion pumps, light-driven ion pumps, and electron-driven ion pumps, are summarized. The working mechanisms, functions, and applications of those artificial ion pumping systems are discussed. Finally, a brief conclusion of underpinning challenges and outlook for future research are tentatively discussed.

Angstrom-scale ion channels towards single-ion selectivity

Artificial ion channels with ion permeability and selectivity comparable to their biological counterparts are highly desired for efficient separation, biosensing, and energy conversion technologies. In the past two decades, both nanoscale and sub-nanoscale ion channels have been successfully fabricated to mimic biological ion channels. Although nanoscale ion channels have achieved intelligent gating and rectification properties, they cannot realize high ion selectivity, especially single-ion selectivity. Artificial angstrom-sized ion channels with narrow pore sizes <1 nm and well-defined pore structures mimicking biological channels have accomplished high ion conductivity and single-ion selectivity. This review comprehensively summarizes the research progress in the rational design and synthesis of artificial subnanometer-sized ion channels with zero-dimensional to three-dimensional pore structures. Then we discuss cation/anion, mono-/di-valent cation, mono-/di-valent anion, and single-ion selectivities of the synthetic ion channels and highlight their potential applications in high-efficiency ion separation, energy conversion, and biological therapeutics. The gaps of single-ion selectivity between artificial and natural channels and the connections between ion selectivity and permeability of synthetic ion channels are covered. Finally, the challenges that need to be addressed in this research field and the perspective of angstrom-scale ion channels are discussed.

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.

The Archaeal Na+/Ca2+ Exchanger (NCX_Mj) as a Model of Ion Transport for the Superfamily of Ca2+/CA Antiporters

The superfamily of Calcium/Cation (Ca²⁺/CA) antiporters extrude Ca²⁺ from the cytosol or subcellular compartments in exchange with Na⁺, K⁺, H⁺, Li⁺, or Mg²⁺ and thereby provide a key mechanism for Ca²⁺ signaling and ion homeostasis in biological systems ranging from bacteria to humans. The structure-dynamic determinants of ion selectivity and transport rates remain unclear, although this is of primary physiological significance. Despite wide variances in the ion selectivity and transport rates, the Ca²⁺/CA proteins share structural motifs, although it remains unclear how the ion recognition/binding is coupled to the ion translocation events. Here, the archaeal Na⁺/Ca²⁺ exchanger (NCX_Mj) is considered as a structure-based model that can help to resolve the ion transport mechanisms by using X-ray, HDX-MS, ATR-FTIR, and computational approaches in conjunction with functional analyses of mutants. Accumulating data reveal that the local backbone dynamics at ion-coordinating residues is characteristically constrained in apo NCX_Mj, which may predefine the affinity and stability of ion-bound species in the ground and transition states. The 3Na⁺ or 1Ca²⁺ binding to respective sites of NCX_Mj rigidify the backbone dynamics at specific segments, where the ion-dependent compression of the ion-permeating four-helix bundle (TM2, TM3, TM7, and TM8) induces the sliding of the two-helix cluster (TM1/TM6) on the protein surface to switch the OF (outward-facing) and IF (inward-facing) conformations. Taking into account the common structural elements shared by Ca²⁺/CAs, NCX_Mj may serve as a model for studying the structure-dynamic and functional determinants of ion-coupled alternating access, transport catalysis, and ion selectivity in Ca²⁺/CA proteins.

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.

Pore-forming transmembrane domains control ion selectivity and selectivity filter conformation in the KirBac1.1 potassium channel

Potassium (K⁺) channels are membrane proteins with the remarkable ability to very selectively conduct K⁺ ions across the membrane. High-resolution structures have revealed that dehydrated K⁺ ions permeate through the narrowest region of the pore, formed by the backbone carbonyls of the signature selectivity filter (SF) sequence TxGYG. However, the existence of nonselective channels with similar SF sequences, as well as effects of mutations in other regions on selectivity, suggest that the SF is not the sole determinant of selectivity. We changed the selectivity of the KirBac1.1 channel by introducing mutations at residue I131 in transmembrane helix 2 (TM2). These mutations increase Na⁺ flux in the absence of K⁺ and introduce significant proton conductance. Consistent with K⁺ channel crystal structures, single-molecule FRET experiments show that the SF is conformationally constrained and stable in high-K⁺ conditions but undergoes transitions to dilated low-FRET states in high-Na⁺/low-K⁺ conditions. Relative to wild-type channels, I131M mutants exhibit marked shifts in the K⁺ and Na⁺ dependence of SF dynamics to higher K⁺ and lower Na⁺ concentrations. These results illuminate the role of I131, and potentially other structural elements outside the SF, in controlling ion selectivity, by suggesting that the physical interaction of these elements with the SF contributes to the relative stability of the constrained K⁺-induced SF configuration versus nonselective dilated conformations.

Non-extensitivity and criticality of atomic hydropathicity around a voltage-gated sodium channel's pore: a modeling study

Voltage-gated sodium channels (NavChs) are pore-forming membrane proteins that regulate the transport of sodium ions through the cell membrane. Understanding the structure and function of NavChs is of major biophysical, as well as clinical, importance given their key role in cellular pathophysiology. In this work, we provide a computational framework for modeling system-size-dependent, i.e., cumulative, atomic properties around a NavCh's pore. We illustrate our methodologies on the bacterial NavAb channel captured in a closed-pore state where we demonstrate that the atomic environment around its pore exhibits a bi-phasic spatial organization dictated by the structural separation of the pore domains (PDs) from the voltage-sensing domains (VSDs). Accordingly, a mathematical model describing packing of atoms around NavAb's pore is constructed that allows-under certain conservation conditions-for a power-law approximation of the cumulative hydropathic dipole field effect acting along NavAb's pore. This verified the non-extensitivity hypothesis for the closed-pore NavAb channel and revealed a long-range hydropathic interactions law regulating atom-packing around the NavAb's selectivity filter. Our model predicts a PDs-VSDs coupling energy of [Formula: see text] kcal/mol corresponding to a global maximum of the atom-packing energy profile. Crucially, we demonstrate for the first time how critical phenomena can emerge in a single-channel structure as a consequence of the non-extensive character of its atomic porous environment.

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel

Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Zinc transporters and their functional integration in mammalian cells

Zinc is a ubiquitous biological metal in all living organisms. The spatiotemporal zinc dynamics in cells provide crucial cellular signaling opportunities, but also challenges for intracellular zinc homeostasis with broad disease implications. Zinc transporters play a central role in regulating cellular zinc balance and subcellular zinc distributions. The discoveries of two complementary families of mammalian zinc transporters (ZnTs and ZIPs) in the mid-1990s spurred much speculation on their metal selectivity and cellular functions. After two decades of research, we have arrived at a biochemical description of zinc transport. However, in vitro functions are fundamentally different from those in living cells, where mammalian zinc transporters are directed to specific subcellular locations, engaged in dedicated macromolecular machineries, and connected with diverse cellular processes. Hence, the molecular functions of individual zinc transporters are reshaped and deeply integrated in cells to promote the utilization of zinc chemistry to perform enzymatic reactions, tune cellular responsiveness to pathophysiologic signals, and safeguard cellular homeostasis. At present, the underlying mechanisms driving the functional integration of mammalian zinc transporters are largely unknown. This knowledge gap has motivated a shift of the research focus from in vitro studies of purified zinc transporters to in cell studies of mammalian zinc transporters in the context of their subcellular locations and protein interactions. In this review, we will outline how knowledge of zinc transporters has been accumulated from in-test-tube to in-cell studies, highlighting new insights and paradigm shifts in our understanding of the molecular and cellular basis of mammalian zinc transporter functions.

Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force Fields

In recent years, the K2P family of potassium channels has been the subject of intense research activity. Owing to the complex function and regulation of this family of ion channels, it is common practice to complement experimental findings with the atomistic description provided by computational approaches such as molecular dynamics (MD) simulations, especially, in light of the unprecedented timescales accessible at present. However, despite recent substantial improvements, the accuracy of MD simulations is still undermined by the intrinsic limitations of force fields. Here, we systematically assessed the performance of the most popular force fields employed to study ion channels at timescales that are orders of magnitude greater than the ones accessible when these energy functions were first developed. Using 32 μs of trajectories, we investigated the dynamics of a member of the K2P ion channel family, the TRAAK channel, using two established force fields in simulations of biological systems: AMBER and CHARMM. We found that while results are comparable on the nanosecond timescales, significant inconsistencies arise at microsecond timescales.

A Series of Discrete Dichloride Dihydrates: Characterisation and Symmetry Effects

A series of three discrete dichloride dihydrates [Cl₂ (H₂ O)₂ ]^2- have been isolated with different triaminocyclopropenium (TAC) cations and with different crystallographic symmetries. The cluster exhibits D_2h symmetry with the tris(dimethylamino)cyclopropenium cation [C₃ (NMe₂ )₃ ]⁺ , C_2h symmetry with the fluorinated cation [C₃ (N(CH₂ CF₃ )₂ )(NBu₂ )₂ ]⁺ (containing two 2,2,2-trifluoroethyl substituents) and C_2v symmetry with the more fluorinated [C₃ (N(CH₂ CF₃ )₂ )₂ (NBu₂ )]⁺ cation. The effect of symmetry on the infrared spectra of the dichloride ion-pair clusters, as well as deuterated analogues, has been investigated. The D_2h - and C_2h -symmetric clusters each exhibit two stretching bands in the infrared at 3427 and 3368 cm^-1 for D_2h symmetry and 3444 and 3392 cm^-1 for C_2h symmetry, whereas the C_2v -symmetric cluster exhibits three bands at 3475, 3426 and 3373 cm^-1 . Computational studies were carried out on a [Cl₂ (H₂ O)₂ ]^2- cluster with C_2v symmetry to aid the infrared band assignments.

One water to bind a chloride-chloride ion pair: isolation of discrete [Cl2(H2O)]2- in the solid state

A discrete dichloride ion pair in the form of a monohydrate, [Cl₂(H₂O)]^2-, was isolated using the triaminocyclopropenium cation [C₃(NHex₂)(N(CH₂CF₃)₂)₂]⁺. Although this ion pair is calculated to be unstable in the gas phase, the ionic lattice and weak CH-Cl hydrogen bonds assist the stabilization of the cluster. The D₂O and HDO isotopomers were also prepared and characterized.

Bicarbonate permeation through anion channels: its role in health and disease

Many anion channels, frequently referred as Cl^- channels, are permeable to different anions in addition to Cl^-. As the second-most abundant anion in the human body, HCO₃^- permeation via anion channels has many important physiological roles. In addition to its classical role as an intracellular pH regulator, HCO₃^- also controls the activity and stability of dissolved proteins in bodily fluids such as saliva, pancreatic juice, intestinal fluid, and airway surface liquid. Moreover, HCO₃^- permeation through these channels affects membrane potentials that are the driving forces for transmembrane transport of solutes and water in epithelia and affect neuronal excitability in nervous tissue. Consequently, aberrant HCO₃^- transport via anion channels causes a number of human diseases in respiratory, gastrointestinal, genitourinary, and neuronal systems. Notably, recent studies have shown that the HCO₃^- permeabilities of several anion channels are not fixed and can be altered by cellular stimuli, findings which may have both physiological and pathophysiological significance. In this review, we summarize recent progress in understanding the molecular mechanisms and the physiological roles of HCO₃^- permeation through anion channels. We hope that the present discussions can stimulate further research into this very important topic, which will provide the basis for human disorders associated with aberrant HCO₃^- transport.

Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor Channel

Ion channel proteins form water-filled nanoscale pores within lipid bilayers, and their properties are dependent on the complex behavior of water in a nanoconfined environment. Using a simplified model of the pore of the 5-HT3 receptor (5HT3R) which restrains the backbone structure to that of the parent channel protein from which it is derived, we compare additive with polarizable models in describing the behavior of water in nanopores. Molecular dynamics simulations were performed with four conformations of the channel: two closed state structures, an intermediate state, and an open state, each embedded in a phosphatidylcholine bilayer. Water density profiles revealed that for all water models, the closed and intermediate states exhibited strong dewetting within the central hydrophobic gate region of the pore. However, the open state conformation exhibited varying degrees of hydration, ranging from partial wetting for the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. Water dipole moments calculated using polarizable force fields also revealed that water molecules remaining within dewetted sections of the pore resemble gas phase water. Free energy profiles for Na+ and for Cl- ions within the open state pore revealed more rugged energy landscapes using polarizable force fields, and the hydration number profiles of these ions were also sensitive to induced polarization resulting in a substantive reduction of the number of waters within the first hydration shell of Cl- while it permeates the pore. These results demonstrate that induced polarization can influence the complex behavior of water and ions within nanoscale pores and provides important new insights into their chemical properties.

A Discrete Dichloride Tetrahydrate Trapped by a Cyclopropenium Cation: Structure and Spectroscopic Properties

A discrete dichloride tetrahydrate cluster, [Cl₂ (H₂ O)₄ ]^2- , was obtained as a salt of the bis(diphenylamino)diethylamino cyclopropenium cation [C₃ (NPh₂ )₂ (NEt₂ )]⁺ and characterized by single-crystal X-ray diffraction and infrared spectroscopy. This chloride-chloride ion-pair cluster consists of a [Cl₂ (H₂ O)₂ ]^2- square with opposite edges bridged by water molecules to give a chair-like structure of the non-hydrogen atoms. The solid-state structure is essentially the same as the calculated gas-phase structure. Infrared spectra were also collected on the deuterium analogue [Cl₂ (D₂ O)₄ ]^2- . Computational studies were carried out on gas-phase [Cl₂ (H₂ O)₄ ]^2- to confirm the infrared band assignments in the solid state. The structure and infrared spectrum are consistent with the discrete nature of the cluster.

Hydrogen-Deuterium Exchange Mass-Spectrometry of Secondary Active Transporters: From Structural Dynamics to Molecular Mechanisms

Membrane transporters allow the selective transport of otherwise poorly permeable solutes across the cell membrane and thus, play a key role in maintaining cellular homeostasis in all kingdoms of life. Importantly, these proteins also serve as important drug targets. Over the last decades, major progress in structural biology methods has elucidated important structure-function relationships in membrane transporters. However, structures obtained using methods such as X-ray crystallography and high-resolution cryogenic electron microscopy merely provide static snapshots of an intrinsically dynamic, multi-step transport process. Therefore, there is a growing need for developing new experimental approaches capable of exploiting the data obtained from the high-resolution snapshots in order to investigate the dynamic features of membrane proteins. Here, we present the basic principles of hydrogen-deuterium exchange mass-spectrometry (HDX-MS) and recent advancements in its use to study membrane transporters. In HDX-MS experiments, minute amounts of a protein sample can be used to investigate its structural dynamics under native conditions, without the need for chemical labelling and with practically no limit on the protein size. Thus, HDX-MS is instrumental for resolving the structure-dynamic landscapes of membrane proteins in their apo (ligand-free) and ligand-bound forms, shedding light on the molecular mechanism underlying the transport process and drug binding.

Selectivity of the KcsA potassium channel: Analysis and computation

Ion channels regulate the flux of ions through cell membranes and play significant roles in many physiological functions. Most of the existing literature focuses on computational approaches based on molecular dynamics simulation or numerical solution of the modified Poisson-Nernst-Planck (PNP) system. In this paper, we present an analytical and computational study of a mathematical model of the KcsA potassium channel, including the effects of ion size (Bikerman model) and solvation energy (Born model). Under equilibrium conditions, we obtain an analytical solution of our modified PNP system, which is used to explain selectivity of KcsA of various ions (K^{+}, Na^{+}, Cl^{-}, Ca^{2+}, and Ba^{2+}) due to negative permanent charges inside the filter region and the effect of ion sizes. Our results show that K^{+} is always selected over Na^{+}, as smaller Na^{+} ions have larger solvation energy. As the amount of negative charges in the filter exceeds a critical value, divalent ions (Ca^{2+} and Ba^{2+}) can enter the filter region and block the KcsA channel. For the nonequilibrium cases, due to difficulties associated with a pure analytical or numerical approach, we use a hybrid analytical-numerical method to solve the modified PNP system. Our predictions of selectivity of KcsA channels and saturation phenomenon of the current-voltage (I-V) curve agree with experimental observations.

Pore-modulating toxins exploit inherent slow inactivation to block K+ channels

Voltage-dependent potassium channels (K_vs) gate in response to changes in electrical membrane potential by coupling a voltage-sensing module with a K⁺-selective pore. Animal toxins targeting K_vs are classified as pore blockers, which physically plug the ion conduction pathway, or as gating modifiers, which disrupt voltage sensor movements. A third group of toxins blocks K⁺ conduction by an unknown mechanism via binding to the channel turrets. Here, we show that Conkunitzin-S1 (Cs1), a peptide toxin isolated from cone snail venom, binds at the turrets of K_v1.2 and targets a network of hydrogen bonds that govern water access to the peripheral cavities that surround the central pore. The resulting ectopic water flow triggers an asymmetric collapse of the pore by a process resembling that of inherent slow inactivation. Pore modulation by animal toxins exposes the peripheral cavity of K⁺ channels as a novel pharmacological target and provides a rational framework for drug design.

Water molecules mediate zinc mobility in the bacterial zinc diffusion channel ZIPB

Regulated ion diffusion across biological membranes is vital for cell function. In a nanoscale ion channel, the active role of discrete water molecules in modulating hydrodynamic behaviors of individual ions is poorly understood because of the technical challenge of tracking water molecules through the channel. Here we report the results of a hydroxyl radical footprinting analysis of the zinc-selective channel ZIPB from the Gram-negative bacterium, Bordetella bronchiseptica Irradiating ZIPB by microsecond X-ray pulses activated water molecules to form covalent hydroxyl radical adducts at nearby residues, which were identified by bottom-up proteomics to detect residues that interact either with zinc or water in response to zinc binding. We found a series of residues exhibiting reciprocal changes in water accessibility attributed to alternating zinc and water binding. Mapping these residues to the previously reported crystal structure of ZIPB, we identified a water-reactive pathway that superimposed on a zinc translocation pathway consisting of two binuclear metal centers and an interim zinc-binding site. The cotranslocation of zinc and water suggested that pore-lining residues undergo a mode switch between zinc coordination and water binding to confer zinc mobility. The unprecedented details of water-mediated zinc transport identified here highlight an essential role of solvated waters in driving zinc coordination dynamics and transmembrane crossing.

Potassium channel selectivity filter dynamics revealed by single-molecule FRET

Potassium (K) channels exhibit exquisite selectivity for conduction of K⁺ ions over other cations, particularly Na⁺. High-resolution structures reveal an archetypal selectivity filter (SF) conformation in which dehydrated K⁺ ions, but not Na⁺ ions, are perfectly coordinated. Using single-molecule FRET (smFRET), we show that the SF-forming loop (SF-loop) in KirBac1.1 transitions between constrained and dilated conformations as a function of ion concentration. The constrained conformation, essential for selective K⁺ permeability, is stabilized by K⁺ but not Na⁺ ions. Mutations that render channels nonselective result in dilated and dynamically unstable conformations, independent of the permeant ion. Further, while wild-type KirBac1.1 channels are K⁺ selective in physiological conditions, Na⁺ permeates in the absence of K⁺. Moreover, whereas K⁺ gradients preferentially support ⁸⁶Rb⁺ fluxes, Na⁺ gradients preferentially support ²²Na⁺ fluxes. This suggests differential ion selectivity in constrained versus dilated states, potentially providing a structural basis for this anomalous mole fraction effect.

Allosteric Modalities for Membrane-Bound Receptors: Insights from Drug Hunting for Brain Diseases

Medicinal chemists are accountable for embedding the appropriate drug target profile into the molecular architecture of a clinical candidate. An accurate characterization of the functional effects following binding of a drug to its biological target is a fundamental step in the discovery of new medicines, informing the translation of preclinical efficacy and safety observations into human trials. Membrane-bound proteins, particularly ion channels and G protein-coupled receptors (GPCRs), are biological targets prone to allosteric modulation. Investigations using allosteric drug candidates and chemical tools suggest that their functional effects may be tailored with a high degree of translational alignment, making them molecular tools to correct pathophysiological functional tone and enable personalized medicine when a causative target-to-disease link is known. We present select examples of functional molecular fine-tuning of allosterism and discuss consequences relevant to drug design.

Mechanistic insights into the Orai channel by molecular dynamics simulations

Highly Ca²⁺ selective channels trigger a large variety of cellular signaling processes in both excitable and non-excitable cells. Among these channels, the Orai channel is unique in its activation mechanism and its structure. It mediates Ca²⁺ influx into the cytosol with an extremely small unitary conductance over longer time-scales, ranging from minutes up to several hours. Its activation is regulated by the Ca²⁺ content of the endoplasmic reticulum (ER). Depletion of luminal [Ca²⁺]_ER is sensed by the STIM1 single transmembrane protein that directly binds and gates the Orai1 channel. Orai mediated Ca²⁺ influx increases cytosolic Ca²⁺ from 100 nM up to low micromolar range close to the pore and thereby forms Ca²⁺ microdomains. Hence, these features of the Orai channel can trigger long-term signaling processes without affecting the overall Ca²⁺ content of a single living cell. Here we focus on the architecture and dynamic conformational changes within the Orai channel. This review summarizes current achievements of molecular dynamics simulations in combination with live cell recordings to address gating and permeation of the Orai channel with molecular precision.

Ion-peptide interactions between alkali metal ions and a termini-protected dipeptide: modeling a portion of the selectivity filter in K+ channels

Potassium channels have the unique ability to allow the selective passage of potassium ions at near diffusion-free rates while inhibiting the passage of more abundant sodium ions. Local interactions between chemical functional groups and the ions are responsible for both selectivity and transport. As an initial step in characterizing these interactions, the structures of Na+ and K+ complexed to the Ac-Tyr-NHMe peptide have been determined from infrared laser spectroscopy and supporting ab initio calculations. Ac-Tyr-NHMe, a termini-protected peptide sequence, replicates the GYG portion of one of the four peptide chains comprising the selectivity filter of a K+ channel. This peptide contains two carbonyl groups, among the eight C[double bond, length as m-dash]O groups forming the S1 binding site of the selectivity filter. Three conformations have been identified for both ions by laser IR-IR double resonance methods. Two conformations have the ion bound to the two C[double bond, length as m-dash]O groups. The third conformation has, in addition, a cation-π interaction with the aromatic ring of tyrosine, i.e. tridentate binding. The relative contributions of the three conformers are approximately the same for K+Ac-Tyr-NHMe, while the tridentate conformer is preferred for Na+Ac-Tyr-NHMe. These differences will be discussed in the context of ion mobility and selectivity.

Patch-Clamp Analysis of the Mitochondrial Calcium Uniporter

Mitochondria accumulate significant amounts of calcium when cytosolic calcium is elevated above the resting levels of 50-100 nM during signaling events. This calcium uptake is primarily mediated by a macromolecular protein assembly called mitochondrial calcium uniporter (MCU) that resides in the mitochondrial inner membrane. In 2004, we applied patch-clamp technique for the first time to record calcium currents from the mitochondrial inner membrane and proved unequivocally that MCU is a highly selective calcium channel. This chapter describes how patch-clamp technique can be applied to record the Ca²⁺ uniporter currents from the mitochondrial inner membrane, isolation of mitochondria from the heart tissue, and preparation of mitoplast using French Press. We also discuss advantages of patch-clamp technique as compared to other methods of determining mitochondrial uniporter activity and important considerations in applying patch-clamp technique to such a small subcellular organelle. With small variations in the bath and pipette solution composition, the same methodology can be applied to study any other currents (e.g., H⁺ or Cl^-) from the mitochondrial inner membrane.

Cation-specific interactions of protein surface charges in dilute aqueous salt solutions: a combined study using dielectric relaxation spectroscopy and Raman spectroscopy

We exploited glycine as a zwitterionic model system to experimentally probe the cation specific interaction of protein surface charges in dilute (≤0.25 mol L-1) aqueous solutions of four biologically relevant inorganic salts, NaCl, KCl, MgCl2 and CaCl2, via dielectric relaxation spectroscopy (DRS) and Raman spectroscopy. Glycine is the simplest building block of proteins and it exposes the same charged groups (carboxylate and ammonium) to the solvent that dominate the protein-water interface. As a counter ion, we selected Cl- due to its biological importance. For all systems, we performed simultaneous fitting of the real (ε') and imaginary (ε″) parts of the dielectric functions, assuming a multimodal relaxation model, obtained from concentration dependent dielectric measurements at ∼293 K. We observe a reduction of the dielectric amplitude for the glycine relaxation while the corresponding time constant shows only small (<7%) deviations compared to aqueous glycine solutions. We propose that the observed reduction in dielectric amplitude is due to a reduction of the effective dipole moment (µeff) of zwitterionic glycine caused by the interaction of glycine with the ion even at very low (0.05 M) salt concentrations. The interaction between divalent metal ions and zwitterionic glycine is increased compared to the monovalent cation-zwitterion interaction; a finding that is also supported by Raman spectroscopy. Our combined dielectric relaxation and Raman spectroscopic study indicates that ion-glycine interactions are weak and mediated by the solvent. Cation-specificity of protein surface charges is also observed in dilute salt solutions (≤0.25 mol L-1), where electrostatic interactions dominate.

Intramembrane ionic protein-lipid interaction regulates integrin structure and function

Protein transmembrane domains (TMDs) are generally hydrophobic, but our bioinformatics analysis shows that many TMDs contain basic residues at terminal regions. Physiological functions of these membrane-snorkeling basic residues are largely unclear. Here, we show that a membrane-snorkeling Lys residue in integrin αLβ2 (also known as lymphocyte function-associated antigen 1 [LFA-1]) regulates transmembrane heterodimer formation and integrin adhesion through ionic interplay with acidic phospholipids and calcium ions (Ca²⁺) in T cells. The amino group of the conserved Lys ionically interacts with the phosphate group of acidic phospholipids to stabilize αLβ2 transmembrane association, thus keeping the integrin at low-affinity conformation. Intracellular Ca²⁺ uses its charge to directly disrupt this ionic interaction, leading to the transmembrane separation and the subsequent extracellular domain extension to increase adhesion activity. This Ca²⁺-mediated regulation is independent on the canonical Ca²⁺ signaling or integrin inside-out signaling. Our work therefore showcases the importance of intramembrane ionic protein–lipid interaction and provides a new mechanism of integrin activation.

Integrin αLβ2 is the major integrin in T cells and plays a vital role in regulating T-cell activation, adhesion, and migration. The transmembrane association of αL and β2 is crucial for maintaining the integrin at low-affinity conformation. Here, we find that the conserved basic residue (K702) in the transmembrane domain of β2 contributes to transmembrane association through ternary ionic interaction with acidic phospholipid and αL cytoplasmic residue. Upon T-cell activation, influxed calcium ions (Ca²⁺) can directly disrupt the ionic K702–lipid interaction through its positive charges, which leads to transmembrane separation and subsequent extracellular domain extension to switch αLβ2 to high-affinity conformation. This Ca²⁺-mediated regulation is through the modulation of the ionic Lys–lipid interaction but not through the canonical Ca²⁺ signaling or integrin inside-out signaling. Our study thus reports a new regulatory mechanism of integrin activation and showcases the importance of intramembrane ionic protein–lipid interaction. This finding might have general relevance, as bioinformatics analysis shows the presence of membrane-snorkeling basic residue is a common feature of transmembrane proteins.

Comparison of permeation mechanisms in sodium-selective ion channels

Voltage-gated sodium channels are the molecular components of electrical signaling in the body, yet the molecular origins of Na⁺-selective transport remain obscured by diverse protein chemistries within this family of ion channels. In particular, bacterial and mammalian sodium channels are known to exhibit similar relative ion permeabilities for Na⁺ over K⁺ ions, despite their distinct signature EEEE and DEKA sequences. Atomic-level molecular dynamics simulations using high-resolution bacterial channel structures and mammalian channel models have begun to describe how these sequences lead to analogous high field strength ion binding sites that drive Na⁺ conduction. Similar complexes have also been identified in unrelated acid sensing ion channels involving glutamate and aspartate side chains that control their selectivity. These studies suggest the possibility of a common origin for Na⁺ selective binding and transport.

Ion Permeation Mechanism in Epithelial Calcium Channel TRVP6

Calcium is the most abundant metal in the human body that plays vital roles as a cellular electrolyte as well as the smallest and most frequently used signaling molecule. Calcium uptake in epithelial tissues is mediated by tetrameric calcium-selective transient receptor potential (TRP) channels TRPV6 that are implicated in a variety of human diseases, including numerous forms of cancer. We used TRPV6 crystal structures as templates for molecular dynamics simulations to identify ion binding sites and to study the permeation mechanism of calcium and other ions through TRPV6 channels. We found that at low Ca²⁺ concentrations, a single calcium ion binds at the selectivity filter narrow constriction formed by aspartates D541 and allows Na⁺ permeation. In the presence of ions, no water binds to or crosses the pore constriction. At high Ca²⁺ concentrations, calcium permeates the pore according to the knock-off mechanism that includes formation of a short-lived transition state with three calcium ions bound near D541. For Ba²⁺, the transition state lives longer and the knock-off permeation occurs slower. Gd³⁺ binds at D541 tightly, blocks the channel and prevents Na⁺ from permeating the pore. Our results provide structural foundations for understanding permeation and block in tetrameric calcium-selective ion channels.

Chloride Ion Transport by the E. coli CLC Cl-/H+ Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study

We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl^- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl^-/H⁺ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl^-. Moving the anion into the external binding site (S_ext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl^- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl^- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl^- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl^- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl^- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl^- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl^- transport.

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.

Identification of Intracellular β-Barrel Residues Involved in Ion Selectivity in the Mechanosensitive Channel of Thermoanaerobacter tengcongensis

The mechanosensitive channel of small conductance (MscS) is a bacterial membrane pore that senses membrane tension and protects cells from lysis by releasing osmolytes. MscS is a homoheptameric channel with a cytoplasmic domain with seven portals and a β-barrel opening to the cytoplasm. TtMscS, an MscS channel from Thermoanaerobacter tengcongensis, is an anion-selective channel. A previous study from our laboratory has defined the crucial role of β-barrel in the anion selectivity of TtMscS (Zhang et al., 2012). However, the mechanistic details by which the β-barrel determines anion selectivity remain unclear. Here, using mutagenesis and patch-clamp recordings, we investigated the function and structural correlations between β-barrels and the anion selectivity of TtMscS at the atomic level. Our results indicated that mutation of V274, a residue at the center of the inner wall of the β-barrel in TtMscS, caused the anion selectivity of TtMscS reverse to cation selectivity. Moreover, the electrostatic potential (T272) and physical size (L276) of residues in the inner wall of β-barrel also determine the anion selectivity of TtMscS. In summary, the present study confirmed that the β-barrel region of TtMscS acts as a “selective filter” that renders TtMscS anion selectivity.