Two mechanisms of ion selectivity in protein binding sites

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

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

Alkali metal cations modulate the geometry of different binding sites in HCN4 selectivity filter for permeation or block

Hyperpolarization-activated cyclic-nucleotide gated (HCN) channels are important for timing biological processes like heartbeat and neuronal firing. Their weak cation selectivity is determined by a filter domain with only two binding sites for K+ and one for Na+. The latter acts as a weak blocker, which is released in combination with a dynamic widening of the filter by K+ ions, giving rise to a mixed K+/Na+ current. Here, we apply molecular dynamics simulations to systematically investigate the interactions of five alkali metal cations with the filter of the open HCN4 pore. Simulations recapitulate experimental data like a low Li+ permeability, considerable Rb+ conductance, a block by Cs+ as well as a punch through of Cs+ ions at high negative voltages. Differential binding of the cation species in specific filter sites is associated with structural adaptations of filter residues. This gives rise to ion coordination by a cation-characteristic number of oxygen atoms from the filter backbone and solvent. This ion/protein interplay prevents Li+, but not Na+, from entry into and further passage through the filter. The site equivalent to S3 in K+ channels emerges as a preferential binding and presumably blocking site for Cs+. Collectively, the data suggest that the weak cation selectivity of HCN channels and their block by Cs+ are determined by restrained cation-generated rearrangements of flexible filter residues.

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.

Stable potassium isotopes (41K/39K) track transcellular and paracellular potassium transport in biological systems

As the most abundant cation in archaeal, bacterial, and eukaryotic cells, potassium (K⁺) is an essential element for life. While much is known about the machinery of transcellular and paracellular K transport–channels, pumps, co-transporters, and tight-junction proteins—many quantitative aspects of K homeostasis in biological systems remain poorly constrained. Here we present measurements of the stable isotope ratios of potassium (⁴¹K/³⁹K) in three biological systems (algae, fish, and mammals). When considered in the context of our current understanding of plausible mechanisms of K isotope fractionation and K⁺ transport in these biological systems, our results provide evidence that the fractionation of K isotopes depends on transport pathway and transmembrane transport machinery. Specifically, we find that passive transport of K⁺ down its electrochemical potential through channels and pores in tight-junctions at favors ³⁹K, a result which we attribute to a kinetic isotope effect associated with dehydration and/or size selectivity at the channel/pore entrance. In contrast, we find that transport of K⁺ against its electrochemical gradient via pumps and co-transporters is associated with less/no isotopic fractionation, a result that we attribute to small equilibrium isotope effects that are expressed in pumps/co-transporters due to their slower turnover rate and the relatively long residence time of K⁺ in the ion pocket. These results indicate that stable K isotopes may be able to provide quantitative constraints on transporter-specific K⁺ fluxes (e.g., the fraction of K efflux from a tissue by channels vs. co-transporters) and how these fluxes change in different physiological states. In addition, precise determination of K isotope effects associated with K⁺ transport via channels, pumps, and co-transporters may provide unique constraints on the mechanisms of K transport that could be tested with steered molecular dynamic simulations.

ION BEHAVIOUR IN THE SELECTIVITY FILTER OF HCN1 CHANNELS

Hyperpolarization-activated cyclic nucleotide-gated channels (HCNs) are responsible for the generation of pacemaker currents (I_f or I_h) in cardiac and neuronal cells. Despite the overall structural similarity to voltage-gated potassium (Kv) channels, HCNs show much lower selectivity for K⁺ over Na⁺ ions. This increased permeability to Na⁺ is critical to their role in membrane depolarization. HCNs can also select between Na⁺ and Li⁺ ions. Here we investigate the unique ion selectivity properties of HCNs using molecular dynamics simulations. Our simulations suggest that the HCN1 pore is flexible and dilated compared to Kv channels with only one stable ion binding site within the selectivity filter. We also observe that ion coordination and hydration differ within the HCN1 selectivity filter compared to those in Kv and CNG channels. Additionally, the C358T mutation further stabilizes the symmetry of the binding site and provides a more fit space for ion coordination, particularly for Li⁺. STATEMENT OF SIGNIFICANCE: Hyperpolarization-activated cyclic-nucleotide gated (HCN) channels represent the molecular correlate of the currents I_f or I_h in cardiomyocytes and neurons. Here we study the unique low conductance and semi-selective properties of HCNs. The conductance and selectivity mechanisms of ion channels are tightly associated with their physiological role and contribute to the specific properties of the excitable cells in which they are expressed.

Cooperative Metal Ion Coordination to the Short Self-Assembling Peptide Promotes Hydrogelation and Cellular Proliferation

Non-covalent interactions among short peptides and proteins led to their molecular self-assembly into supramolecular packaging, which provides the fundamental basis of life. These biomolecular assemblies are highly susceptible to the environmental conditions, including temperature, light, pH, and ionic concentration, thus inspiring the fabrication of a new class of stimuli-responsive biomaterials. Here, we report for the first time the cooperative effect of the divalent metal ions to promote hydrogelation in the short collagen inspired self-assembling peptide for developing advanced biomaterials. Introduction of the biologically relevant metal ions (Ca²⁺ /Mg²⁺ ) to the peptide surpasses its limitation to self-assemble into a multi-scale structure at physiological pH. In particular, in presence of metal ions, the negatively charged peptide showed a distinct shift in its equilibrium point of gelation and demonstrated conversion from sol to gel and thus enabling the scope of fabricating an advanced biomaterial for controlling cellular behaviour. Interestingly, tunable mechanical strength and improved cellular response were observed within ion-coordinated peptide hydrogels compared to the peptide gelator. Microscopic analyses, rheological assessment, and biological studies established the importance of utilizing a novel strategy by simply using metal ions to modulate the physical and biological attributes of CIPs to construct next-generation biomaterials. This article is protected by copyright. All rights reserved.

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.

Water-Mediated Interactions Enhance Alkaline Earth Cation Chelation in Neighboring Cavities of a Cytosine Quartet in the DNA Quadruplex

Larger Coulombic repulsion between divalent cations compared to the monovalent counterparts dictates the cation-cation distance in the central ion channel of quadruplexes. In this work, density functional theory and a continuum solvation model were employed to study bond energies of alkaline earth cations in adjacent cavities of the central ion channel. Four crystallized tetramolecular quadruplexes with various geometric constraints and structural motifs available in the Protein Data Bank were examined in order to understand how the cation binding affinities could be increased in aqueous solution. A cytosine quartet sandwiched between guanine quartets has a larger bond energy of the second alkaline earth cation in comparison with guanine and uracil quartets. Four highly conserved hydrogen-bonded water molecules in the center of the cytosine quartet are responsible for a higher electrostatic interaction with the cations in comparison with guanines' carbonyl groups. The reported findings are valuable for the design of synthetic quadruplexes templated with divalent cations for optoelectronic applications.

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.

Ion Selectivity and Permeation Mechanism in a Cyclodextrin-Based Channel

Synthetic ion channels are a promising technology in the medical and materials sciences because of their ability to conduct ions. Channels based on cyclodextrin, a cyclic oligomer of glucose, are of particular interest because of their nontoxicity and biocompatibility. Using molecular dynamics-based free energy calculations, this study identifies cyclodextrin channel types that are best suited to serve as synthetic ion channels. Free energy profiles show that the connectivity in the channel determines whether the channel is cation-selective or anion-selective. Furthermore, the energy barrier for ion transport is governed by the number of glucose molecules making up the cyclodextrin units of the channel. A detailed mechanism is proposed for ion transport through these channels. Findings from this study will aid in designing cyclodextrin-based channels that could be either cation-selective or anion-selective, by modifying the linkages of the channel or the number of glucose molecules in the cyclodextrin rings.

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.

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.

Bases of Bacterial Sodium Channel Selectivity Among Organic Cations

Hille's (1971) seminal study of organic cation selectivity of eukaryotic voltage-gated sodium channels showed a sharp size cut-off for ion permeation, such that no ion possessing a methyl group was permeant. Using the prokaryotic channel, NaChBac, we found some similarity and two peculiar differences in the selectivity profiles for small polyatomic cations. First, we identified a diverse group of minimally permeant cations for wildtype NaChBac, ranging in sizes from ammonium to guanidinium and tetramethylammonium; and second, for both ammonium and hydrazinium, the charge-conserving selectivity filter mutation (E191D) yielded substantial increases in relative permeability (PX/PNa). The relative permeabilities varied inversely with relative Kd calculated from 1D Potential of Mean Force profiles (PMFs) for the single cations traversing the channel. Several of the cations bound more strongly than Na+, and hence appear to act as blockers, as well as charge carriers. Consistent with experimental observations, the E191D mutation had little impact on Na+ binding to the selectivity filter, but disrupted the binding of ammonium and hydrazinium, consequently facilitating ion permeation across the NaChBac-like filter. We concluded that for prokaryotic sodium channels, a fine balance among filter size, binding affinity, occupancy, and flexibility seems to contribute to observed functional differences.

Thermodynamics of Cation Binding to the Sarcoendoplasmic Reticulum Calcium ATPase Pump and Impacts on Enzyme Function

Sarcoendoplasmic reticulum Ca²⁺-ATPase (SERCA) is a transmembrane pump that plays an important role in transporting calcium into the sarcoplasmic reticulum (SR). While calcium (Ca²⁺) binds SERCA with micromolar affinity, magnesium (Mg²⁺) and potassium (K⁺) also compete with Ca²⁺ binding. However, the molecular bases for these competing ions' influence on the SERCA function and the selectivity of the pump for Ca²⁺ are not well-established. We therefore used in silico methods to resolve molecular determinants of cation binding in the canonical site I and II Ca²⁺ binding sites via (1) triplicate molecular dynamics (MD) simulations of Mg²⁺, Ca²⁺, and K⁺-bound SERCA, (2) mean spherical approximation (MSA) theory to score the affinity and selectivity of cation binding to the MD-resolved structures, and (3) state models of SERCA turnover informed from MSA-derived affinity data. Our key findings are that (a) coordination at sites I and II is optimized for Ca²⁺ and to a lesser extent for Mg²⁺ and K⁺, as determined by MD-derived cation-amino acid oxygen and bound water configurations, (b) the impaired coordination and high desolvation cost for Mg²⁺ precludes favorable Mg²⁺ binding relative to Ca²⁺, while K⁺ has limited capacity to bind site I, and (c) Mg²⁺ most likely acts as inhibitor and K⁺ as intermediate in SERCA's reaction cycle, based on a best-fit state model of SERCA turnover. These findings provide a quantitative basis for SERCA function that leverages molecular-scale thermodynamic data and rationalizes enzyme activity across broad ranges of K⁺, Ca²⁺, and Mg²⁺ concentrations.

Different roles for aspartates and glutamates for cation permeation in bacterial sodium channels

A key driving force for ion channel selectivity is represented by the negative charge of the Selectivity Filter carried by aspartate (D) and glutamate (E) residues. However, the structural effects and specific properties of D and E residues have not been extensively studied. In order to investigate this issue we studied the mutants of NaChBac channel with all possible combinations of D and E in the charged rings in position 191 and 192. Electrophysiological measurements showed significant Ca²⁺ currents only when position 191 was occupied by E. Equilibrium Molecular Dynamics simulations revealed the existence of two binding sites, corresponding to the charged rings and another one, more internal, at the level of L190. The simulations showed that the ion in the innermost site can interact with the residue in position 191 only when this is glutamate. Based on the MD simulations, we suggest that a D in position 191 leads to a high affinity Ca²⁺ block site resulting from a significant drop in the free energy of binding for an ion moving between the binding sites; in contrast, the free energy change is more gradual when an E residue occupies position 191, resulting in Ca²⁺ permeability. This scenario is consistent with the model of ion channel selectivity through stepwise changes in binding affinity proposed by Dang and McCleskey. Our study also highlights the importance of the structure of the selectivity filter which should contribute to the development of more detailed physical models for ion channel selectivity.

Many-body effect determines the selectivity for Ca2+ and Mg2+ in proteins

Calcium ion is a versatile messenger in many cell-signaling processes. To achieve their functions, calcium-binding proteins selectively bind Ca2+ against a background of competing ions such as Mg2+ The high specificity of calcium-binding proteins has been intriguing since Mg2+ has a higher charge density than Ca2+ and is expected to bind more tightly to the carboxylate groups in calcium-binding pockets. Here, we showed that the specificity for Ca2+ is dictated by the many-body polarization effect, which is an energetic cost arising from the dense packing of multiple residues around the metal ion. Since polarization has stronger distance dependence compared with permanent electrostatics, the cost associated with the smaller Mg2+ is much higher than that with Ca2+ and outweighs the electrostatic attraction favorable for Mg2+ With the AMOEBA (atomic multipole optimized energetics for biomolecular simulation) polarizable force field, our simulations captured the relative binding free energy between Ca2+ and Mg2+ for proteins with various types of binding pockets and explained the nonmonotonic size dependence of the binding free energy in EF-hand proteins. Without electronic polarization, the smaller ions are always favored over larger ions and the relative binding free energy is roughly proportional to the net charge of the pocket. The many-body effect depends on both the number and the arrangement of charged residues. Fine-tuning of the ion selectivity could be achieved by combining the many-body effect and geometric constraint.

H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel

Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K⁺ ion through KcsA.

Determinants of the host-guest interactions between α-, β- and γ-cyclodextrins and group IA, IIA and IIIA metal cations: a DFT/PCM study

The most widely used native cyclodextrins are α-, β- and γ-cyclodextrins containing six, seven or eight α-d-glucopyranoside units in the ring, respectively. Although the ligation properties of these host molecules have been extensively studied, a number of questions regarding their metal binding and selectivity remain unaddressed: to what extent do the size and flexibility of the host α-, β- and γ-cyclodextrins influence their metal affinity/selectivity? Which metal is the most preferred binding partner of α-, β- and γ-cyclodextrins? How do the charge, size and preferred coordination number of the metal cation shape its interactions with the host cyclodextrin? Can the guest metal cation inflict structural alterations in the host molecule and, if so, how do these changes correlate with the metal's properties? In the present study, by employing density functional theory (DFT) calculations combined with polarizable continuum model (PCM) computations, we try to answer these questions by evaluating the thermodynamic parameters of the IA, IIA and IIIA group metal binding to α, β- and γ-cyclodextrins. We assess how the interaction between the two binding partners depends on (1) the size, valence state and preferred coordination number of the guest metal cations, (2) the size and flexibility of the host molecule, and (3) the dielectric properties of the environment. The series of group IA (Na⁺ and Rb⁺), IIA (Mg²⁺ and Sr²⁺) and IIIA (Al³+ and In³⁺) metal cations have been chosen for the task as they allow us to study the effect of various metal parameters (variable charge, ionic radius and coordination number) on the strength and form of the interactions with the host cyclodextrins.

Ion channels and ion selectivity

Specific macromolecular transport systems, ion channels and pumps, provide the pathways to facilitate and control the passage of ions across the lipid membrane. Ion channels provide energetically favourable passage for ions to diffuse rapidly and passively according to their electrochemical potential. Selective ion channels are essential for the excitability of biological membranes: the action potential is a transient phenomenon that reflects the rapid opening and closing of voltage-dependent Na⁺-selective and K⁺-selective channels. One of the most critical functional aspects of K⁺ channels is their ability to remain highly selective for K⁺ over Na⁺ while allowing high-throughput ion conduction at a rate close to the diffusion limit. Permeation through the K⁺ channel selectivity filter is believed to proceed as a 'knockon' mechanism, in which 2-3 K⁺ ions interspersed by water molecules move in a single file. Permeation through the comparatively wider and less selective Na⁺ channels also proceeds via a loosely coupled knockon mechanism, although the ions do not need to be fully dehydrated. While simple structural concepts are often invoked to rationalize the mechanism of ion selectivity, a deeper analysis shows that subtle effects play an important role in these flexible dynamical structures.

The Environment Shapes the Inner Vestibule of LeuT

Human neurotransmitter transporters are found in the nervous system terminating synaptic signals by rapid removal of neurotransmitter molecules from the synaptic cleft. The homologous transporter LeuT, found in Aquifex aeolicus, was crystallized in different conformations. Here, we investigated the inward-open state of LeuT. We compared LeuT in membranes and micelles using molecular dynamics simulations and lanthanide-based resonance energy transfer (LRET). Simulations of micelle-solubilized LeuT revealed a stable and widely open inward-facing conformation. However, this conformation was unstable in a membrane environment. The helix dipole and the charged amino acid of the first transmembrane helix (TM1A) partitioned out of the hydrophobic membrane core. Free energy calculations showed that movement of TM1A by 0.30 nm was driven by a free energy difference of ~15 kJ/mol. Distance measurements by LRET showed TM1A movements, consistent with the simulations, confirming a substantially different inward-open conformation in lipid bilayer from that inferred from the crystal structure.

Crystal structures of the bacterial small amino acid transporter LeuT provided structural evidence for the alternating access model. Thereby, these structures shaped our understanding of the mechanisms underlying substrate translocation by neurotransmitter transporters. However, it has been questioned, if the crystallized inward-open conformation of LeuT can exist in the membrane environment. Here we show that, while stable in detergent micelles, the inward-open conformation of LeuT is of high energy and undergoes structural readjustments. We use a multi-faceted approach including molecular dynamics simulations, scintillation proximity assays, free energy calculations and apply for the first time lanthanide resonance energy transfer measurements to verify the in silico predictions. In silico and in vitro approaches using the same conditions allowed us to combine the macroscopic experimental data with microscopic all atom results from simulations to identify the underlying driving forces: partitioning of charged and polar groups from the hydrophobic membrane interior to the hydrophilic environment. We propose that the inward-facing state shows a much smaller movement of TM1A, but large enough to create an access path to the S1 substrate binding site from the vestibule.

The selectivity of the Na(+)/K(+)-pump is controlled by binding site protonation and self-correcting occlusion

The Na⁺/K⁺-pump maintains the physiological K⁺ and Na⁺ electrochemical gradients across the cell membrane. It operates via an 'alternating-access' mechanism, making iterative transitions between inward-facing (E₁) and outward-facing (E₂) conformations. Although the general features of the transport cycle are known, the detailed physicochemical factors governing the binding site selectivity remain mysterious. Free energy molecular dynamics simulations show that the ion binding sites switch their binding specificity in E₁ and E₂. This is accompanied by small structural arrangements and changes in protonation states of the coordinating residues. Additional computations on structural models of the intermediate states along the conformational transition pathway reveal that the free energy barrier toward the occlusion step is considerably increased when the wrong type of ion is loaded into the binding pocket, prohibiting the pump cycle from proceeding forward. This self-correcting mechanism strengthens the overall transport selectivity and protects the stoichiometry of the pump cycle.

DOI:http://dx.doi.org/10.7554/eLife.16616.001

A protein called the sodium-potassium pump resides in the membrane that surrounds living cells. The role of this protein is to 'pump' sodium and potassium ions across the membrane to help restore their concentration inside and outside of the cell. About 25% of the body's energy is used to keep the pump going, rising to nearly 70% in the brain. Problems that affect the pump have been linked to several disorders, including heart, kidney and metabolic diseases, as well as severe neurological conditions.

The sodium-potassium pump must be able to effectively pick out the correct ions to transport from a mixture of many different types of ions. However, it was not clear how the pump succeeds in doing this efficiently. Rui et al. have now used a computational method called molecular dynamics simulations to model how the sodium-potassium pump transports the desired ions across the cell membrane.

The pump works via a so-called 'alternating-access' mechanism, repeatedly transitioning between inward-facing and outward-facing conformations. In each cycle, it binds three sodium ions from the cell’s interior and exports them to the outside. Then, the pump binds to two potassium ions from outside the cell and imports them inside. Although the bound sodium and potassium ions interact with similar binding sites in the pump, the pump sometimes preferentially binds sodium, and sometimes potassium. The study performed by Rui et al. shows that this preference is driven by how protons (hydrogen ions) bind to the amino acids that make up the binding site.

The simulations also suggest that the pump uses a ‘self-correcting’ mechanism to prevent the pump from transporting the wrong types of ions. When incorrect ions are present at the binding sites, the pump cycle pauses temporarily until the ions detach from the pump. Only when the correct ions are bound will the pump cycle continue again.

In the future, Rui et al. hope to use long time-scale molecular dynamics simulations to show the conformational transition in action. In addition, the 'self-correcting' mechanism will be directly tested by letting the wrong and correct ions compete for the binding sites to see whether the pump will transport only the correct ions.

DOI:http://dx.doi.org/10.7554/eLife.16616.002

LiCl solvation in N-methyl-acetamide (NMA) as a model for understanding Li(+) binding to an amide plane

The thermodynamics of ion solvation in non-aqueous solvents remains of great significance for understanding cellular transport and ion homeostasis for the design of novel ion-selective materials and applications in molecular pharmacology. Molecular simulations play pivotal roles in connecting experimental measurements to the microscopic structures of liquids. One of the most useful and versatile mimetic systems for understanding biological ion transport is N-methyl-acetamide (NMA). A plethora of theoretical studies for ion solvation in NMA have appeared recently, but further progress is limited by two factors. One is an apparent lack of experimental data on solubility and thermodynamics of solvation for a broad panel of 1 : 1 salts over an appropriate temperature and concentration range. The second concern is more substantial and has to do with the limitations hardwired in the additive (fixed charge) approximations used for most of the existing force-fields. In this submission, we report on the experimental evaluation of LiCl solvation in NMA over a broad range of concentrations and temperatures and compare the results with those of MD simulations with several additive and one polarizable force-field (Drude). By comparing our simulations and experimental results to density functional theory computations, we discuss the limiting factors in existing potential functions. To evaluate the possible implications of explicit and implicit polarizability treatments on ion permeation across biological channels, we performed potential of mean force (PMF) computations for Li(+) transport through a model narrow ion channel with additive and polarizable force-fields.

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled-stranded RNA to deactivate interferon. The interface of VP35-dsRNA would be a feasible target for structure-based antiviral agent design. To explore the essence of VP35-dsRNA interaction, molecular dynamics simulation combined with MM-GBSA calculations were performed on Marburg virus VP35-dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35-dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35-dsRNA recognition mechanism will provide some insights for development of antiviral drug.

Pore dimensions and the role of occupancy in unitary conductance of Shaker K channels

The resistance of the inner vestibule limits Shaker’s conductance.

K channels mediate the selective passage of K⁺ across the plasma membrane by means of intimate interactions with ions at the pore selectivity filter located near the external face. Despite high conservation of the selectivity filter, the K⁺ transport properties of different K channels vary widely, with the unitary conductance spanning a range of over two orders of magnitude. Mutation of Pro475, a residue located at the cytoplasmic entrance of the pore of the small-intermediate conductance K channel Shaker (Pro475Asp (P475D) or Pro475Gln (P475Q)), increases Shaker’s reported ∼20-pS conductance by approximately six- and approximately threefold, respectively, without any detectable effect on its selectivity. These findings suggest that the structural determinants underlying the diversity of K channel conductance are distinct from the selectivity filter, making P475D and P475Q excellent probes to identify key determinants of the K channel unitary conductance. By measuring diffusion-limited unitary outward currents after unilateral addition of 2 M sucrose to the internal solution to increase its viscosity, we estimated a pore internal radius of capture of ∼0.82 Å for all three Shaker variants (wild type, P475D, and P475Q). This estimate is consistent with the internal entrance of the Kv1.2/2.1 structure if the effective radius of hydrated K⁺ is set to ∼4 Å. Unilateral exposure to sucrose allowed us to estimate the internal and external access resistances together with that of the inner pore. We determined that Shaker resistance resides mainly in the inner cavity, whereas only ∼8% resides in the selectivity filter. To reduce the inner resistance, we introduced additional aspartate residues into the internal vestibule to favor ion occupancy. No aspartate addition raised the maximum unitary conductance, measured at saturating [K⁺], beyond that of P475D, suggesting an ∼200-pS conductance ceiling for Shaker. This value is approximately one third of the maximum conductance of the large conductance K (BK) channel (the K channel of highest conductance), reducing the energy gap between their K⁺ transport rates to ∼1 kT. Thus, although Shaker’s pore sustains ion translocation as the BK channel’s does, higher energetic costs of ion stabilization or higher friction with the ion’s rigid hydration cage in its narrower aqueous cavity may entail higher resistance.

Elucidating factors important for monovalent cation selectivity in enzymes: E. coli β-galactosidase as a model

Many enzymes require a specific monovalent cation (M(+)), that is either Na(+) or K(+), for optimal activity. While high selectivity M(+) sites in transport proteins have been extensively studied, enzyme M(+) binding sites generally have lower selectivity and are less characterized. Here we study the M(+) binding site of the model enzyme E. coli β-galactosidase, which is about 10 fold selective for Na(+) over K(+). Combining data from X-ray crystallography and computational models, we find the electrostatic environment predominates in defining the Na(+) selectivity. In this lower selectivity site rather subtle influences on the electrostatic environment become significant, including the induced polarization effects of the M(+) on the coordinating ligands and the effect of second coordination shell residues on the charge distribution of the primary ligands. This work expands the knowledge of ion selectivity in proteins to denote novel mechanisms important for the selectivity of M(+) sites in enzymes.

Two Na+ Sites Control Conformational Change in a Neurotransmitter Transporter Homolog

In LeuT, a prokaryotic homolog of neurotransmitter transporters, Na(+) stabilizes outward-open conformational states. We examined how each of the two LeuT Na(+) binding sites contributes to Na(+)-dependent closure of the cytoplasmic pathway using biochemical and biophysical assays of conformation. Mutating either of two residues that contribute to the Na2 site completely prevented cytoplasmic closure in response to Na(+), suggesting that Na2 is essential for this conformational change, whereas Na1 mutants retained Na(+) responsiveness. However, mutation of Na1 residues also influenced the Na(+)-dependent conformational change in ways that varied depending on the position mutated. Computational analyses suggest those mutants influence the ability of Na1 binding to hydrate the substrate pathway and perturb an interaction network leading to the extracellular gate. Overall, the results demonstrate that occupation of Na2 stabilizes outward-facing conformations presumably through a direct interaction between Na(+) and transmembrane helices 1 and 8, whereas Na(+) binding at Na1 influences conformational change through a network of intermediary interactions. The results also provide evidence that N-terminal release and helix motions represent distinct steps in cytoplasmic pathway opening.

Mechanism of the Association between Na+ Binding and Conformations at the Intracellular Gate in Neurotransmitter:Sodium Symporters

Neurotransmitter:sodium symporters (NSSs) terminate neurotransmission by Na⁺-dependent reuptake of released neurotransmitters. Previous studies suggested that Na⁺-binding reconfigures dynamically coupled structural elements in an allosteric interaction network (AIN) responsible for function-related conformational changes, but the intramolecular pathway of this mechanism has remained uncharted. We describe a new approach for the modeling and analysis of intramolecular dynamics in the bacterial NSS homolog LeuT. From microsecond-scale molecular dynamics simulations and cognate experimental verifications in both LeuT and human dopamine transporter (hDAT), we apply the novel method to identify the composition and the dynamic properties of their conserved AIN. In LeuT, two different perturbations disrupting Na⁺ binding and transport (i.e. replacing Na⁺ with Li⁺ or the Y268A mutation at the intracellular gate) affect the AIN in strikingly similar ways. In contrast, other mutations that affect the intracellular gate (i.e. R5A and D369A) do not significantly impair Na⁺ cooperativity and transport. Our analysis shows these perturbations to have much lesser effects on the AIN, underscoring the sensitivity of this novel method to the mechanistic nature of the perturbation. Notably, this set of observations holds as well for hDAT, where the aligned Y335A, R60A, and D436A mutations also produce different impacts on Na⁺ dependence. Thus, the detailed AIN generated from our method is shown to connect Na⁺ binding with global conformational changes that are critical for the transport mechanism. That the AIN between the Na⁺ binding sites and the intracellular gate in bacterial LeuT resembles that in eukaryotic hDAT highlights the conservation of allosteric pathways underlying NSS function.

Realizing synchronous energy harvesting and ion separation with graphene oxide membranes

A synchronous ion separation and electricity generation process has been developed using G-O membranes. In addition to the size effect proposed prevsiouly, the separation of ions can be attributed to the different interactions between ions and G-O membranes; the generation of electricity is due to the confinement of G-O membranes, and the mobility difference of ions. Efficient energy transduction has been achieved with G-O membranes, converting magnetic, thermal and osmotic energy to electricity, distinguishing this material from other commercial semi-permeable membranes. Our study indicated that G-O membranes could find potential applications in the purification of wastewater, while producing electricity simultaneously. With G-O membranes, industrial magnetic leakage and waste heat could also be used to produce electricity, affording a superior approach for energy recovery.

Molecular strategies to achieve selective conductance in NaK channel variants

A recent crystallization of several ion channels has provided strong impetus for efforts aimed at understanding the different strategies employed by nature for selective ion transport. In this work, we used two variants of the selectivity filter of NaK channel to explore molecular mechanisms that give rise to K(+)-selectivity. We computed one-dimensional (1D) and two-dimensional (2D) potentials of mean force (PMFs) for ion permeation across the channel. The results indicate that the energies for Na(+) and K(+) permeation across the selectivity filter display significant differences in positions of the binding sites and barriers. One characteristic signature of a K(+)-selective channel is the apparent preservation of the site analogous to that of S2 in KcsA. The S2-bound ion can be almost ideally dehydrated and coordinated by 6 to 8 carbonyls. In a striking contrast, the PMFs controlling transport of ions in a nonselective variant show almost identical profiles for either K(+) or Na(+) and significant involvement of water molecules in ion coordination across the entire selectivity filter. An analysis of differences in 1D PMFs for Na(+) and K(+) suggests that coordination number alone is an insufficient predictor of site selectivity, while chemical composition (ratio of carbonyls and water molecules) correlates well with preference for K(+). Multi-ion effects such as dependence of the barriers and wells for permeant ion on the type of copermeant ion were found to play a significant role in the selectivity signature of the channel as well.

A ribokinase family conserved monovalent cation binding site enhances the MgATP-induced inhibition in E. coli phosphofructokinase-2

The presence of a regulatory site for monovalent cations that affects the conformation of the MgATP-binding pocket leading to enzyme activation has been demonstrated for ribokinases. This site is selective toward the ionic radius of the monovalent cation, accepting those larger than Na(+). Phosphofructokinase-2 (Pfk-2) from Escherichia coli is homologous to ribokinase, but unlike other ribokinase family members, presents an additional site for the nucleotide that negatively regulates its enzymatic activity. In this work, we show the effect of monovalent cations on the kinetic parameters of Pfk-2 together with its three-dimensional structure determined by x-ray diffraction in the presence of K(+) or Cs(+). Kinetic characterization of the enzyme shows that K(+) and Na(+) alter neither the kcat nor the KM values for fructose-6-P or MgATP. However, the presence of K(+) (but not Na(+)) enhances the allosteric inhibition induced by MgATP. Moreover, binding experiments show that K(+) (but not Na(+)) increases the affinity of MgATP in a saturable fashion. In agreement with the biochemical data, the crystal structure of Pfk-2 obtained in the presence of MgATP shows a cation-binding site at the conserved position predicted for the ribokinase family of proteins. This site is adjacent to the MgATP allosteric binding site and is only observed in the presence of Cs(+) or K(+). These results indicate that binding of the monovalent metal ions indirectly influences the allosteric site of Pfk-2 by increasing its affinity for MgATP with no alteration in the conformation of residues present at the catalytic site.

Sodium channel selectivity and conduction: prokaryotes have devised their own molecular strategy

The molecular strategy for alkali cation selectivity by a bacterial sodium channel resembles those of eukaryotic calcium and potassium channels, rather than those of eukaryotic sodium channels.

Striking structural differences between voltage-gated sodium (Nav) channels from prokaryotes (homotetramers) and eukaryotes (asymmetric, four-domain proteins) suggest the likelihood of different molecular mechanisms for common functions. For these two channel families, our data show similar selectivity sequences among alkali cations (relative permeability, Pion/PNa) and asymmetric, bi-ionic reversal potentials when the Na/K gradient is reversed. We performed coordinated experimental and computational studies, respectively, on the prokaryotic Nav channels NaChBac and NavAb. NaChBac shows an “anomalous,” nonmonotonic mole-fraction dependence in the presence of certain sodium–potassium mixtures; to our knowledge, no comparable observation has been reported for eukaryotic Nav channels. NaChBac’s preferential selectivity for sodium is reduced either by partial titration of its highly charged selectivity filter, when extracellular pH is lowered from 7.4 to 5.8, or by perturbation—likely steric—associated with a nominally electro-neutral substitution in the selectivity filter (E191D). Although no single molecular feature or energetic parameter appears to dominate, our atomistic simulations, based on the published NavAb crystal structure, revealed factors that may contribute to the normally observed selectivity for Na over K. These include: (a) a thermodynamic penalty to exchange one K⁺ for one Na⁺ in the wild-type (WT) channel, increasing the relative likelihood of Na⁺ occupying the binding site; (b) a small tendency toward weaker ion binding to the selectivity filter in Na–K mixtures, consistent with the higher conductance observed with both sodium and potassium present; and (c) integrated 1-D potentials of mean force for sodium or potassium movement that show less separation for the less selective E/D mutant than for WT. Overall, tight binding of a single favored ion to the selectivity filter, together with crucial inter-ion interactions within the pore, suggests that prokaryotic Nav channels use a selective strategy more akin to those of eukaryotic calcium and potassium channels than that of eukaryotic Nav channels.

The two Na+ sites in the human serotonin transporter play distinct roles in the ion coupling and electrogenicity of transport

Neurotransmitter transporters of the SLC6 family of proteins, including the human serotonin transporter (hSERT), utilize Na(+), Cl(-), and K(+) gradients to induce conformational changes necessary for substrate translocation. Dysregulation of ion movement through monoamine transporters has been shown to impact neuronal firing potentials and could play a role in pathophysiologies, such as depression and anxiety. Despite multiple crystal structures of prokaryotic and eukaryotic SLC transporters indicating the location of both (or one) conserved Na(+)-binding sites (termed Na1 and Na2), much remains uncertain in regard to the movements and contributions of these cation-binding sites in the transport process. In this study, we utilize the unique properties of a mutation of hSERT at a single, highly conserved asparagine on TM1 (Asn-101) to provide several lines of evidence demonstrating mechanistically distinct roles for Na1 and Na2. Mutations at Asn-101 alter the cation dependence of the transporter, allowing Ca(2+) (but not other cations) to functionally replace Na(+) for driving transport and promoting 5-hydroxytryptamine (5-HT)-dependent conformational changes. Furthermore, in two-electrode voltage clamp studies in Xenopus oocytes, both Ca(2+) and Na(+) illicit 5-HT-induced currents in the Asn-101 mutants and reveal that, although Ca(2+) promotes substrate-induced current, it does not appear to be the charge carrier during 5-HT transport. These findings, in addition to functional evaluation of Na1 and Na2 site mutants, reveal separate roles for Na1 and Na2 and provide insight into initiation of the translocation process as well as a mechanism whereby the reported SERT stoichiometry can be obtained despite the presence of two putative Na(+)-binding sites.

Stereoselective inhibition of serotonin transporters by antimalarial compounds

The serotonin (5-HT) transporter (SERT) is an integral membrane protein that functions to reuptake 5-HT released into the synapse following neurotransmission. This role serves an important regulatory mechanism in neuronal homeostasis. Previous studies have demonstrated that several clinically important antimalarial compounds inhibit serotonin (5-hydroxytryptamine, 5-HT) reuptake. In this study, we examined the details of antimalarial inhibition of 5-HT transport in both Drosophila (dSERT) and human SERT (hSERT) using electrophysiologic, biochemical and computational approaches. We found that the cinchona alkaloids quinidine and cinchonine, which have identical stereochemistry about carbons 8 and 9, exhibited the greatest inhibition of dSERT and hSERT transporter function whereas quinine and cinchonidine, enantiomers of quinidine and cinchonine, respectively, were weaker inhibitors of dSERT and hSERT. Furthermore, SERT mutations known to decrease the binding affinity of many antidepressants affected the cinchona alkaloids in a stereo-specific manner where the similar inhibitory profiles for quinine and cinchonidine (8S,9R) were distinct from quinidine and cinchonine (8R,9S). Small molecule docking studies with hSERT homology models predict that quinine and cinchonidine bind to the central 5-HT binding site (S1) whereas quinidine and cinchonine bind to the S2 site. Taken together, the data presented here support binding of cinchona alkaloids to two different sites on SERT defined by their stereochemistry which implies separate modes of transporter inhibition. Notably, the most potent antimalarial inhibitors of SERT appear to preferentially bind to the S2 site. Our findings provide important insight related to how this class of drugs can modulate the serotonergic system as well as identify compounds that may discriminate between the S1 and S2 binding sites and serve as lead compounds for novel SERT inhibitors.

Ion permeation in the NanC porin from Escherichia coli: free energy calculations along pathways identified by coarse-grain simulations

Using the X-ray structure of a recently discovered bacterial protein, the N-acetylneuraminic acid-inducible channel (NanC), we investigate computationally K(+) and Cl(-) ions' permeation. We identify ion permeation pathways that are likely to be populated using coarse-grain Monte Carlo simulations. Next, we use these pathways as reaction coordinates for umbrella sampling-based free energy simulations. We find distinct tubelike pathways connecting specific binding sites for K(+) and, more pronounced, for Cl(-) ions. Both ions permeate the porin preserving almost all of their first hydration shell. The calculated free energy barriers are G(#) ≈ 4 kJ/mol and G(#) ≈ 8 kJ/mol for Cl(-) and K(+), respectively. Within the approximations associated with these values, discussed in detail in this work, we suggest that the porin is slightly selective for Cl(-) versus K(+). Our suggestion is consistent with the experimentally observed weak Cl(-) over K(+) selectivity. A rationale for the latter is suggested by a comparison with previous calculations on strongly anion selective porins.

Absolute free energies of biomolecules from unperturbed ensembles

Computing the absolute free energy of a macromolecule's structural state, F, is a challenging problem of high relevance. This study presents a method that computes F using only information from an unperturbed simulation of the macromolecule in the relevant conformational state, ensemble, and environment. Absolute free energies produced by this method, dubbed Valuation of Local Configuration Integral with Dynamics (VALOCIDY), enable comparison of alternative states. For example, comparing explicitly solvated and vaporous states of amino acid side-chain analogs produces solvation free energies in good agreement with experiments. Also, comparisons between alternative conformational states of model heptapeptides (including the unfolded state) produce free energy differences in agreement with data from μs molecular-dynamics simulations and experimental propensities. The potential of using VALOCIDY in computational protein design is explored via a small design problem of stabilizing a β-turn structure. When VALOCIDY-based estimation of folding free energy is used as the design metric, the resulting sequence folds into the desired structure within the atomistic force field used in design. The VALOCIDY-based approach also recognizes the distinct status of the native sequence regardless of minor details of the starting template structure, in stark contrast with a traditional fixed-backbone approach.

SLC6 transporters: structure, function, regulation, disease association and therapeutics

The SLC6 family of secondary active transporters are integral membrane solute carrier proteins characterized by the Na(+)-dependent translocation of small amino acid or amino acid-like substrates. SLC6 transporters, which include the serotonin, dopamine, norepinephrine, GABA, taurine, creatine, as well as amino acid transporters, are associated with a number of human diseases and disorders making this family a critical target for therapeutic development. In addition, several members of this family are directly involved in the action of drugs of abuse such as cocaine, amphetamines, and ecstasy. Recent advances providing structural insight into this family have vastly accelerated our ability to study these proteins and their involvement in complex biological processes.

Relative Free Energies for Hydration of Monovalent Ions from QM and QM/MM Simulations

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.

Pore helices play a dynamic role as integrators of domain motion during Kv11.1 channel inactivation gating

Proteins that form ion-selective pores in the membrane of cells are integral to many rapid signaling processes, including regulating the rhythm of the heartbeat. In potassium channels, the selectivity filter is critical for both endowing an exquisite selectivity for potassium ions, as well as for controlling the flow of ions through the pore. Subtle rearrangements in the complex hydrogen-bond network that link the selectivity filter to the surrounding pore helices differentiate conducting (open) from nonconducting (inactivated) conformations of the channel. Recent studies suggest that beyond the selectivity filter, inactivation involves widespread rearrangements of the channel protein. Here, we use rate equilibrium free energy relationship analysis to probe the structural changes that occur during selectivity filter gating in Kv11.1 channels, at near atomic resolution. We show that the pore helix plays a crucial dynamic role as a bidirectional interface during selectivity filter gating. We also define the molecular bases of the energetic coupling between the pore helix and outer helix of the pore domain that occurs early in the transition from open to inactivated states, as well as the coupling between the pore helix and inner helix late in the transition. Our data demonstrate that the pore helices are more than just static structural elements supporting the integrity of the selectivity filter; instead they play a crucial dynamic role during selectivity filter gating.

An entropic mechanism of generating selective ion binding in macromolecules

Several mechanisms have been proposed to explain how ion channels and transporters distinguish between similar ions, a process crucial for maintaining proper cell function. Of these, three can be broadly classed as mechanisms involving specific positional constraints on the ion coordinating ligands which arise through: a “rigid cavity”, a ‘strained cavity’ and ‘reduced ligand fluctuations’. Each operates in subtly different ways yet can produce markedly different influences on ion selectivity. Here we expand upon preliminary investigations into the reduced ligand fluctuation mechanism of ion selectivity by simulating how a series of model systems respond to a decrease in ligand thermal fluctuations while simultaneously maintaining optimal ion-ligand binding distances. Simple abstract-ligand models, as well as simple models based upon the ion binding sites in two amino acid transporters, show that limiting ligand fluctuations can create ion selectivity between Li⁺, Na⁺ and K⁺ even when there is no strain associated with the molecular framework accommodating the different ions. Reducing the fluctuations in the position of the coordinating ligands contributes to selectivity toward the smaller of two ions as a consequence of entropic differences.

Differentiating between Na⁺ and K⁺ ions is important for many cellular processes, such as nerve conduction and the regulation of membrane potentials. Different biological molecules utilise different methods to discriminate between ions. In this work, the reduced ligand fluctuation mechanism of ion selectivity is described. This entropy-driven mechanism is due to the limited thermal fluctuations of the atoms in some macromolecular ion binding sites. The elucidation of this mechanism offers a more complete picture of the ways in which the fundamental process of ion selectivity can be achieved.