The Split Personality of Glutamate Transporters: A Chloride Channel and a Transporter

Transporters and ion channels are conventionally categorised into distinct classes of membrane proteins. However, some membrane proteins have a split personality and can function as both transporters and ion channels. The excitatory amino acid transporters (EAATs) in particular, function as both glutamate transporters and chloride (Cl^?) channels. The EAATs couple the transport of glutamate to the co-transport of three Na⁺ ions and one H⁺ ion into the cell, and the counter-transport of one K⁺ ion out of the cell. The EAAT Cl^? channel is activated by the binding of glutamate and Na⁺, but is thermodynamically uncoupled from glutamate transport and involves molecular determinants distinct from those responsible for glutamate transport. Several crystal structures of an EAAT archaeal homologue, Glt_Ph, at different stages of the transport cycle, alongside numerous functional studies and molecular dynamics simulations, have provided extensive insights into the mechanism of substrate transport via these transporters. However, the molecular determinants involved in Cl^? permeation, and the mechanism by which this channel is activated are not entirely understood. Here we will discuss what is currently known about the molecular determinants involved in EAAT-mediated Cl^? permeation and the mechanisms that underlie their split personality.     

Molecular mechanism of ion-ion and ion-substrate coupling in the Na+-dependent leucine transporter LeuT

Ion-coupled transport of neurotransmitter molecules by neurotransmitter:sodium symporters (NSS) play an important role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling and the corresponding ion-substrate stoichiometry in NSS transporters has yet to be understood. The recent determination of a high-resolution structure for a bacterial homolog of Na⁺/Cl⁻-dependent neurotransmitter transporters, LeuT, offers a unique opportunity to analyze the functional roles of the multi-ion binding sites within the binding pocket. The binding pocket of LeuT contains two metal binding sites. The first ion in site NA1 is directly coupled to the bound substrate (Leu) with the second ion in the neighboring site (NA2) only ∼7 Å away. Extensive, fully atomistic, molecular dynamics, and free energy simulations of LeuT in an explicit lipid bilayer are performed to evaluate substrate-binding affinity as a function of the ion load (single versus double occupancy) and occupancy by specific monovalent cations. It was shown that double ion occupancy of the binding pocket is required to ensure substrate coupling to Na⁺ and not to Li⁺ or K⁺ cations. Furthermore, it was found that presence of the ion in site NA2 is required for structural stability of the binding pocket as well as amplified selectivity for Na⁺ in the case of double ion occupancy.     

Identification of a 3rd Na+ Binding Site of the Glycine Transporter,GlyT2

The Na⁺/Cl^- dependent glycine transporters GlyT1 and GlyT2 regulate synaptic glycine concentrations. Glycine transport by GlyT2 is coupled to the co-transport of three Na⁺ ions, whereas transport by GlyT1 is coupled to the co-transport of only two Na⁺ ions. These differences in ion-flux coupling determine their respective concentrating capacities and have a direct bearing on their functional roles in synaptic transmission. The crystal structures of the closely related bacterial Na⁺-dependent leucine transporter, LeuT_Aa, and the Drosophila dopamine transporter, dDAT, have allowed prediction of two Na⁺ binding sites in GlyT2, but the physical location of the third Na⁺ site in GlyT2 is unknown. A bacterial betaine transporter, BetP, has also been crystallized and shows structural similarity to LeuT_Aa. Although betaine transport by BetP is coupled to the co-transport of two Na⁺ ions, the first Na⁺ site is not conserved between BetP and LeuT_Aa, the so called Na1'' site. We hypothesized that the third Na⁺ binding site (Na3 site) of GlyT2 corresponds to the BetP Na1'' binding site. To identify the Na3 binding site of GlyT2, we performed molecular dynamics (MD) simulations. Surprisingly, a Na⁺ placed at the location consistent with the Na1'' site of BetP spontaneously dissociated from its initial location and bound instead to a novel Na3 site. Using a combination of MD simulations of a comparative model of GlyT2 together with an analysis of the functional properties of wild type and mutant GlyTs we have identified an electrostatically favorable novel third Na⁺ binding site in GlyT2 formed by Trp263 and Met276 in TM3, Ala481 in TM6 and Glu648 in TM10.     

Mechanism of Cation Binding to the Glutamate Transporter EAAC1 Probed with Mutation of the Conserved Amino Acid Residue Thr101

The glutamate transporter excitatory amino acid carrier 1 (EAAC1) catalyzes the co-transport of three Na⁺ ions, one H⁺ ion, and one glutamate molecule into the cell, in exchange for one K⁺ ion. Na⁺ binding to the glutamate-free form of the transporter generates a high affinity binding site for glutamate and is thus required for transport. Moreover, sodium binding to the transporters induces a basal anion conductance, which is further activated by glutamate. Here, we used the [Na⁺] dependence of this conductance as a read-out of Na⁺ binding to the substrate-free transporter to study the impact of a highly conserved amino acid residue, Thr¹⁰¹, in transmembrane domain 3. The apparent affinity of substrate-free EAAC1 for Na⁺ was dramatically decreased by the T101A but not by the T101S mutation. Interestingly, in further contrast to EAAC1_WT, in the T101A mutant this [Na⁺] dependence was biphasic. This behavior can be explained by assuming that the binding of two Na⁺ ions prior to glutamate binding is required to generate a high affinity substrate binding site. In contrast to the dramatic effect of the T101A mutation on Na⁺ binding, other properties of the transporter, such as its ability to transport glutamate, were impaired but not eliminated. Our results are consistent with the existence of a cation binding site deeply buried in the membrane and involving interactions with the side chain oxygens of Thr¹⁰¹ and Asp³⁶⁷. A theoretical valence screening approach confirms that the predicted site of cation interaction has the potential to be a novel, so far undetected sodium binding site.     

High-Resolution Views and Transport Mechanisms of the NKCC1 and KCC Transporters

Cation-chloride cotransporters (CCCs) are responsible for the coupled co-transport of Cl^- with K⁺ and/or Na⁺ in an electroneutral manner. They play important roles in myriad fundamental physiological processes––from cell volume regulation to transepithelial solute transport and intracellular ion homeostasis––and are targeted by medicines commonly prescribed to treat hypertension and edema. After several decades of studies into the functions and pharmacology of these transporters, there have been several breakthroughs in the structural determination of CCC transporters. The insights provided by these new structures for the Na⁺/K⁺/Cl^- cotransporter NKCC1 and the K⁺/Cl^- cotransporters KCC1, KCC2, KCC3 and KCC4 have deepened our understanding of their molecular basis and transport function. This focused review discusses recent advances in the structural and mechanistic understanding of CCC transporters, including architecture, dimerization, functional roles of regulatory domains, ion binding sites, and coupled ion transport.     

Free energy simulations of ligand binding to the aspartate transporter Glt(Ph)

Glutamate/Aspartate transporters cotransport three Na⁺ and one H⁺ ions with the substrate and countertransport one K⁺ ion. The binding sites for the substrate and two Na⁺ ions have been observed in the crystal structure of the archeal homolog Glt_Ph, while the binding site for the third Na⁺ ion has been proposed from computational studies and confirmed by experiments. Here we perform detailed free energy simulations of Glt_Ph, giving a comprehensive characterization of the substrate and ion binding sites, and calculating their binding free energies in various configurations. Our results show unequivocally that the substrate binds after the binding of two Na⁺ ions. They also shed light into Asp/Glu selectivity of Glt_Ph, which is not observed in eukaryotic glutamate transporters.     

Insights on Na+ binding and conformational dynamics in multidrug and toxic compound extrusion transporter NorM

MATE (multidrug and toxic compound extrusion) transporter proteins mediate metabolite transport in plants and multidrug resistance in bacteria and mammals. MATE transporter NorM from Vibrio cholerae is an antiporter that is driven by Na+ gradient to extrude the substrates. To understand the molecular mechanism of Na+‐substrate exchange, molecular dynamics simulation was performed to study conformational changes of both wild‐type and mutant NorM with and without cation bindings. Our results show that NorM is able to bind two Na⁺ ions simultaneously, one to each of the carboxylic groups of E255 and D371 in the binding pocket. Furthermore, this di‐Na⁺ binding state is likely more efficient for conformational changes of NorM_VC toward the inward‐facing conformation than single‐Na⁺ binding state. The observation of two Na⁺ binding sites of NorM_VC is consistent with the previous study that two sites for ion binding (denoted as Na1/Na2 sites) are found in the transporter LeuT and BetP, another two secondary transporters. Taken together, our findings shed light on the structure rearrangements of NorM on Na⁺ binding and enrich our knowledge of the transport mechanism of secondary transporters. Proteins 2014; 82:240–249. © 2013 Wiley Periodicals, Inc.     

Sodium translocation by the iminoglycinuria associated imino transporter (SLC6A20)

The system IMINO transporter plays an essential role in the transport of proline and hydroxyproline in the intestine and kidney. Its molecular correlate has been identified and named SIT1 or IMINO (SLC6A20). Initial characterization of the transporter showed it to be Na⁺ and Cl^?-dependent, but the stoichiometry remained unresolved. Using homology modeling along the structure of the bacterial leucine transporter LeuT, we identified two highly conserved Na⁺-binding sites and a putative Cl^?-binding site. Mutation of all residues in the two proposed Na⁺-binding sites revealed that most of them were essential for uptake and completely inactivated the transporter. However, mutants A22V (Na⁺-binding site 1) and mutants S20A, S20G, S20G/G405S (Na⁺-binding site 2) were partially active and characterized further. Flux studies suggested that mutations of Na⁺-binding site 1 caused a decrease of the Na⁺-K_0.5, whereas mutations of site 2 increased the K_0.5. Mutation of Na⁺-binding site 1 also changed the ion selectivity of the IMINO transporter. IMINO actively translocates ³⁶Cl^? demonstrating that the proposed chloride binding site is used in the transporter. Accumulation experiments and flux measurements at different holding potentials showed that the transporter can work as a 2Na⁺/1Cl^?-proline cotransporter. The proposed homology model allows to study mutations in IMINO associated with iminoglycinuria.     

Control of ion selectivity in LeuT: two Na+ binding sites with two different mechanisms

The x-ray structure of LeuT, a bacterial homologue of Na⁺/Cl⁻-dependent neurotransmitter transporters, provides a great opportunity to better understand the molecular basis of monovalent cation selectivity in ion-coupled transporters. LeuT possesses two ion binding sites, NA1 and NA2, which are highly selective for Na⁺. Extensive all-atom free-energy molecular dynamics simulations of LeuT embedded in an explicit membrane are performed at different temperatures and various occupancy states of the binding sites to dissect the molecular mechanism of ion selectivity. The results show that the two binding sites display robust selectivity for Na⁺ over K⁺ or Li⁺, the competing ions of most similar radii. Of particular interest, the mechanism primarily responsible for selectivity for each of the two binding sites appears to be different. In NA1, selectivity for Na⁺ over K⁺ arises predominantly from the strong electrostatic field arising from the negatively charged carboxylate group of the leucine substrate coordinating the ion directly. In NA2, which comprises only neutral ligands, selectivity for Na⁺ is enforced by the local structural restraints arising from the hydrogen-bonding network and the covalent connectivity of the polypeptide chain surrounding the ion according to a “snug-fit” mechanism.     

Adaptation of plants to salt stress: the role of the ion transporters

Adaptation to high salinity is achieved by cellular ion homeostasis which involves regulation of toxic sodium ion (Na⁺) and Chloride ion (Cl⁻) uptake, preventing the transport of these ions to the aerial parts of the plants and vacuolar sequestration of these toxic ions. Ion transporters have long been known to play roles in maintaining ion homeostasis. Na⁺ enters the cell through various voltage dependent selective and non-selective ion channels. High Na⁺ concentration in the plasma membrane is balanced either by uptake of potassium ion (K⁺) by various potassium importing channels, by salt exclusion mechanism or by sequestration of Na⁺ in the vacuoles. Therefore, the role of high-affinity potassium transporter, the salt overly sensitive pathway, the most well-defined Na⁺ exclusion pathway that exports Na⁺ from cell into xylem and tonoplast localized cation transporters that compartmentalizes Na⁺ in vacuoles need to be studied in detail and applied to make the plant adaptable to saline soil. Knowledge on the regulation of expression of these transporters by the hormones, microRNAs and other non-coding RNAs can be utilized to manipulate the ion transport. Here, we reviewed paradigm of the ion transporters in salt stress signalling pathways from the recent and past studies aiding transformation of basic knowledge into biotechnological applications to generate engineered salt stress tolerant crops.

     

Properties of an Inward-Facing State of LeuT: Conformational Stability and?Substrate Release

The leucine transporter (LeuT) is a bacterial homolog of the human monoamine transporters, which are important pharmaceutical targets. There are no high-resolution structures of the human transporters available; however, LeuT has been crystallized in several different conformational states. Recently, an inward-facing conformation of LeuT was solved revealing an unexpectedly large movement of transmembrane helix 1a (TM1a). We have performed molecular dynamics simulations of the mutated and wild-type transporter, with and without the cocrystallized Fab antibody fragment, to investigate the properties of this inward-facing conformation in relation to transport by LeuT within the membrane environment. In all of the simulations, local conformational changes with respect to the crystal structure are consistently observed, especially in TM1a. Umbrella sampling revealed a soft potential for TM1a tilting. Furthermore, simulations of inward-facing LeuT with Na⁺ ions and substrate bound suggest that one of the Na⁺ ion binding sites is fully disrupted. Release of alanine and the second Na⁺ ion is also observed, giving insight into the final stage of the translocation process in atomistic detail.     

Ion-Releasing State of a Secondary Membrane Transporter

The crystal structure of Na⁺-coupled galactose symporter (vSGLT) reports the transporter in its substrate-bound state, with a Na⁺ ion modeled in a binding site corresponding to that of a homologous protein, leucine transporter (LeuT). In repeated molecular dynamics simulations, however, we find the Na⁺ ion instable, invariably and spontaneously diffusing out of the transporter through a pathway lined by D189, which appears to facilitate the diffusion of the ion toward the cytoplasm. Further analysis of the trajectories and close structural examination, in particular, comparison of the Na⁺-binding sites of vSGLT and LeuT, strongly indicates that the crystal structure of vSGLT actually represents an ion-releasing state of the transporter. The observed dynamics of the Na⁺ ion, in contrast to the substrate, also suggests that the cytoplasmic release of the Na⁺ ion precedes that of the substrate, thus shedding light on a key step in the transport cycle of this secondary transporter.     

Atomistic Detailed Mechanism and Weak Cation-Conducting Activity of HIV-1 Vpu Revealed by Free Energy Calculations

The viral protein U (Vpu) encoded by HIV-1 has been shown to assist in the detachment of virion particles from infected cells. Vpu forms cation-specific ion channels in host cells, and has been proposed as a potential drug target. An understanding of the mechanism of ion transport through Vpu is desirable, but remains limited because of the unavailability of an experimental structure of the channel. Using a structure of the pentameric form of Vpu – modeled and validated based on available experimental data – umbrella sampling molecular dynamics simulations (cumulative simulation time of more than 0.4 µs) were employed to elucidate the energetics and the molecular mechanism of ion transport in Vpu. Free energy profiles corresponding to the permeation of Na⁺ and K⁺ were found to be similar to each other indicating lack of ion selection, consistent with previous experimental studies. The Ser23 residue is shown to enhance ion transport via two mechanisms: creating a weak binding site, and increasing the effective hydrophilic length of the channel, both of which have previously been hypothesized in experiments. A two-dimensional free energy landscape has been computed to model multiple ion permeation, based on which a mechanism for ion conduction is proposed. It is shown that only one ion can pass through the channel at a time. This, along with a stretch of hydrophobic residues in the transmembrane domain of Vpu, explains the slow kinetics of ion conduction. The results are consistent with previous conductance studies that showed Vpu to be a weakly conducting ion channel.     

Molecular mechanisms of potassium and sodium uptake in plants

Potassium (K⁺) is an essential nutrient and the most abundant cation in plants, whereas the closely related ion sodium (Na⁺) is toxic to most plants at high millimolar concentrations. K⁺ deficiency and Na⁺ toxicity are both major constraints to crop production worldwide. K⁺ counteracts Na⁺ stress, while Na⁺, in turn, can to a certain degree alleviate K⁺ deficiency. Elucidation of the molecular mechanisms of K⁺ and Na⁺ transport is pivotal to the understanding – and eventually engineering – of plant K⁺ nutrition and Na⁺ sensitivity. Here we provide an overview on plant K⁺ transporters with particular emphasis on root K⁺ and Na⁺ uptake. Plant K⁺-permeable cation transporters comprise seven families: Shaker-type K⁺ channels, `two-pore' K⁺ channels, cyclic-nucleotide-gated channels, putative K⁺/H⁺ antiporters, KUP/HAK/KT transporters, HKT transporters, and LCT1. Candidate genes for Na⁺ transport are the KUP/HAK/KTs, HKTs, CNGCs, and LCT1. Expression in heterologous systems, localization in plants, and genetic disruption in plants will provide insight into the roles of transporter genes in K⁺ nutrition and Na⁺ toxicity.     

Identification of the Third Na Site and the Sequence of Extracellular Binding Events in the Glutamate Transporter

The transport cycle in the glutamate transporter (GlT) is catalyzed by the cotransport of three Na⁺ ions. However, the positions of only two of these ions (Na1 and Na2 sites) along with the substrate have been captured in the crystal structures reported for both the outward-facing and the inward-facing states of Glt_ph. Characterizing the third ion binding site (Na3) is necessary for structure-function studies attempting to investigate the mechanism of transport in GlTs at an atomic level, particularly for the determination of the sequence of the binding events during the transport cycle. In this study, we report a series of molecular dynamics simulations performed on various bound states of Glt_ph (the apo state, as well as in the presence of Na⁺, the substrate, or both), which have been used to identify a putative Na3 site. The calculated trajectories have been used to determine the water accessibility of potential ion-binding residues in the protein, as a prerequisite for their ion binding. Combined with conformational analysis of the key regions in the protein in different bound states and several additional independent simulations in which a Na⁺ ion was randomly introduced to the interior of the transporter, we have been able to characterize a putative Na3 site and propose a plausible binding sequence for the substrate and the three Na⁺ ions to the transporter during the extracellular half of the transport cycle. The proposed Na3 site is formed by a set of highly conserved residues, namely, Asp³¹², Thr⁹², and Asn³¹⁰, along with a water molecule. Simulation of a fully bound state, including the substrate and the three Na⁺ ions, reveals a stable structure—showing closer agreement to the crystal structure when compared to previous models lacking an ion in the putative Na3 site. The proposed sequence of binding events is in agreement with recent experimental models suggesting that two Na⁺ ions bind before the substrate, and one after that. Our results, however, provide additional information about the sites involved in these binding events.     

Dynamics of the extracellular gate and ion-substrate coupling in the glutamate transporter

Glutamate transporters (GluTs) are the primary regulators of extracellular concentration of the neurotransmitter glutamate in the central nervous system. In this study, we have investigated the dynamics and coupling of the substrate and Na⁺ binding sites, and the mechanism of cotransport of Na⁺ ions, using molecular dynamics simulations of a membrane-embedded model of GluT in its apo (empty form) and various Na⁺- and/or substrate-bound states. The results shed light on the mechanism of the extracellular gate and on the sequence of binding of the substrate and Na⁺ ions to GluT during the transport cycle. The results suggest that the helical hairpin HP2 plays the key role of the extracellular gate for the substrate binding site, and that the opening and closure of the gate is controlled by substrate binding. GluT adopts an open conformation in the absence of the substrate exposing the binding sites of the substrate and Na⁺ ions to the extracellular solution. Based on the calculated trajectories, we propose that Na1 is the first element to bind GluT, as it is found to be important for the completion of the substrate binding site. The subsequent binding of the substrate, in turn, is shown to result in an almost complete closure of the extracellular gate and the formation of the Na2 binding site. Finally, binding of Na2 locks the extracellular gate and completes the formation of the occluded state of GluT.     

Microseconds Simulations Reveal a New Sodium-binding Site and the Mechanism of Sodium-coupled Substrate Uptake by LeuT

The bacterial sodium-coupled leucine/alanine transporter LeuT is broadly used as a model system for studying the transport mechanism of neurotransmitters because of its structural and functional homology to mammalian transporters such as serotonin, dopamine, or norepinephrine transporters, and because of the resolution of its structure in different states. Although the binding sites (S1 for substrate, and Na1 and Na2 for two co-transported sodium ions) have been resolved, we still lack a mechanistic understanding of coupled Na⁺- and substrate-binding events. We present here results from extensive (>20 μs) unbiased molecular dynamics simulations generated using the latest computing technology. Simulations show that sodium binds initially the Na1 site, but not Na2, and, consistently, sodium unbinding/escape to the extracellular (EC) region first takes place at Na2, succeeded by Na1. Na2 diffusion back to the EC medium requires prior dissociation of substrate from S1. Significantly, Na⁺ binding (and unbinding) consistently involves a transient binding to a newly discovered site, Na1″, near S1, as an intermediate state. A robust sequence of substrate uptake events coupled to sodium bindings and translocations between those sites assisted by hydration emerges from the simulations: (i) bindings of a first Na⁺ to Na1″, translocation to Na1, a second Na⁺ to vacated Na1″ and then to Na2, and substrate to S1; (ii) rotation of Phe²⁵³ aromatic group to seclude the substrate from the EC region; and (iii) concerted tilting of TM1b and TM6a toward TM3 and TM8 to close the EC vestibule.     

Ion-Controlled Conformational Dynamics in the Outward-Open Transition from an Occluded State of LeuT

Neurotransmitter:sodium symporter (NSS) proteins are secondary Na⁺-driven active transporters that terminate neurotransmission by substrate uptake. Despite the availability of high-resolution crystal structures of a bacterial homolog of NSSs—Leucine Transporter (LeuT)—and extensive computational and experimental structure-function studies, unanswered questions remain regarding the transport mechanisms. We used microsecond atomistic molecular-dynamics (MD) simulations and free-energy computations to reveal ion-controlled conformational dynamics of LeuT in relation to binding affinity and selectivity of the more extracellularly positioned Na⁺ binding site (Na1 site). In the course of MD simulations starting from the occluded state with bound Na⁺, but in the absence of substrate, we find a spontaneous transition of the extracellular vestibule of LeuT into an outward-open conformation. The outward opening is enhanced by the absence of Na1 and modulated by the protonation state of the Na1-associated Glu-290. Consistently, the Na⁺ affinity for the Na1 site is inversely correlated with the extent of outward-open character and is lower than in the occluded state with bound substrate; however, the Na1 site retains its selectivity for Na⁺ over K⁺ in such conformational transitions. To the best of our knowledge, our findings shed new light on the Na⁺-driven transport cycle and on the symmetry in structural rearrangements for outward- and inward-open transitions.     

Contribution of residues in second transmembrane domain of ASIC1a protein to ion selectivity

Acid-sensing ion channels (ASICs) are proton-gated cation-selective channels expressed in the peripheral and central nervous systems. The ion permeation pathway of ASIC1a is defined by residues 426–450 in the second transmembrane (TM2) segment. The gate, formed by the intersection of the TM2 segments, localizes near the extracellular boundary of the plasma membrane. We explored the contribution to ion permeation and selectivity of residues in the TM2 segment of ASIC1a. Studies of accessibility with positively charged methanethiosulfonate reagents suggest that the permeation pathway in the open state constricts below the gate, restricting the passage to large ions. Substitution of residues in the intracellular vestibule at positions 437, 438, 443, or 446 significantly increased the permeability to K⁺ versus Na⁺. ASIC1a shows a selectivity sequence for alkali metals of Na⁺>Li⁺>K⁺≫Rb⁺>Cs⁺. Alanine and cysteine substitutions at position 438 increased, to different extents, the relative permeability to Li⁺, K⁺, Rb⁺, and Cs⁺. For these mutants, ion permeation was not a function of the diameter of the nonhydrated ion, suggesting that Gly-438 encompasses an ion coordination site that is essential for ion selectivity. M437C and A443C mutants showed slightly increased permeability to K⁺, Rb⁺, and Cs⁺, suggesting that substitutions at these positions influence ion discrimination by altering molecular sieving. Our results indicate that ion selectivity is accomplished by the contribution of multiple sites in the pore of ASIC1a.     

Exploring the Ion Selectivity Properties of a Large Number of Simplified Binding Site Models

The ability to discriminate between different cations efficiently is essential for the proper physiological functioning of many membrane transport proteins. One obvious mechanism of ion selectivity is when a binding site is structurally constrained by the protein architecture and its geometry is precisely adapted to fit an ion of a given size. This mechanism is not effective in the case of flexible protein binding sites that are able to deform structurally or to adapt to a bound ion. In this study, the concept of nontrivial ion selectivity arising in a highly flexible protein binding site conceptually represented as a microdroplet of ligands confined to a small volume is explored. The environment imposed by the spatial confinement is a critical feature of the reduced models. A large number of reduced binding site models (1077) comprising typical ion-coordinating ligands (carbonyl, hydroxyl, carboxylate, water) are constructed and characterized for Na⁺/K⁺ and Ca²⁺/Ba²⁺ size selectivity using free energy perturbation molecular dynamics simulations. Free energies are highly correlated with the sum of ion-ligand and ligand-ligand mean interactions, but the relative balance of those two contributions is different for K⁺-selective and Na⁺-selective binding sites. The analysis indicates that both the number and the type of ligands are important factors in ion selectivity.