{'fraction_non_alphanumeric': 0.03732966077123804, 'fraction_numerical': 0.015149318643084952, 'mean_word_length': 5.083933558724092, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Design principles for K+ selectivity in membrane transport

Selection among competing alternatives is always interesting, but when organisms select between K and Na for transport across biological membranes, it is especially intriguing for a couple of reasons (Hille, 2001). First, the results are physiologically significant. Second, K and Na are nearly as similar as they could be while not being the same things. Perhaps the separation of isotopes of the same chemical species is more demanding. Despite decades of research, the question remains: how do K-selective ion channels catalyze K movement but recognize the detailed molecular-scale differences of Na and discriminate against it? Suggestions based on channel size and coordination chemistry (Bezanilla and Armstrong, 1972; Eisenman and Horn, 1983) have been available for a long time, but the determination of a KcsA K channel crystal structure (Doyle et al., 1998) enabled molecularly specific modeling studies of this K/Na selectivity. In the subsequent flood of computational studies, finding consistency in results and interpretations has proven challenging. Here, we describe our perspective on how molecular modeling has advanced our understanding of the specific chemical and structural design elements of biological molecules that enable selective ion transport.

P e r s p e c t i v e Selection among competing alternatives is always interesting, but when organisms select between K + and Na + for transport across biological membranes, it is especially intriguing for a couple of reasons (Hille, 2001). First, the results are physiologically significant. Second, K + and Na + are nearly as similar as they could be while not being the same things. Perhaps the separation of isotopes of the same chemical species is more demanding. Despite decades of research, the question remains: how do K + -selective ion channels catalyze K + movement but recognize the detailed molecular-scale differences of Na + and discriminate against it? Suggestions based on channel size and coordination chemistry (Bezanilla and Armstrong, 1972;Eisenman and Horn, 1983) have been available for a long time, but the determination of a KcsA K + channel crystal structure (Doyle et al., 1998) enabled molecularly specific modeling studies of this K + /Na + selectivity. In the subsequent flood of computational studies, finding consistency in results and interpretations has proven challenging. Here, we describe our perspective on how molecular modeling has advanced our understanding of the specific chemical and structural design elements of biological molecules that enable selective ion transport.

Background
The simplest idea for K + /Na + discrimination is just that these ions differ in size. Because permeation of the K + ion is preferred over the slightly smaller Na + ion, the possibility that K + /Na + selectivity occurs because the channel is too narrow is quickly discarded. Instead, we must consider the energetic barriers that discourage Na + ions from crossing the membrane (Bezanilla and Armstrong, 1972). Clearly, nature has capitalized on the size difference between the two ions to create energetic barriers for Na + , leaving K + with a less resistive conduction pathway, but how?
Complications for explaining selectivity arise because the family of K + channels exhibits a range of selectivities, from weakly selective HCN pacemaker K + channels to strongly selective "maxi"-type K + channels or their bacterial homologue, KcsA. Fig. 1 A illustrates the structure of a representative K + channel , the family of proteins primarily responsible for passive K + -selective transport. The narrowest region of the pore, referred to as the selectivity filter, is understood to impart K + selectivity, although other portions of the channel also play a role. Diverse K + channels typically share this filter architecture (Hille, 2001) according to sequence alignment (Shealy et al., 2003), x-ray structures of a variety of K + channels, and a vast amount of physiological data. In light of this structural similarity, how does the variability in selectivity come about?
The length of the selectivity filter with several evident ion-binding locations, and the known multiple-ion occupancy of the selectivity filter (Doyle et al., 1998;Åqvist and Luzhkov, 2000;Bernèche and Roux, 2001;Morais-Cabral et al., 2001;Zhou et al., 2001) provide yet further complications for explaining selectivity. Both kinetic modeling and free energy calculations show that multi-ion interactions within the filter affect the selective movement of K + (Åqvist and Luzhkov, 2000;Morais-Cabral et al., 2001). In addition, K + and Na + have multiple binding modes within the dense dipole array of the carbonyl and hydroxyl-lined selectivity filter (Fig. 1, A-C). These underlie the complex behavior of K + channels. For example, in the rare event that Na + overcomes a barrier to entry into the filter, it obstructs K + conductance (Nimigean and Miller, 2002;Thompson et al., 2009). This suggests that, in addition to barriers, the filter also presents Na + with some binding modes that, in fact, compete with K + .
The goal of understanding ion selectivity extends beyond explaining the activities of highly selective K + channels.
In light of these considerations, we present here progress that has resulted from studies that consider simplified model systems and how these studies, alongside traditional whole molecule observations, can advance our understanding of selective K + movement across membranes. We begin our discussion by defining a thermodynamic measure of K + selectivity over Na + .
Free energy of ion selectivity The free energy difference, ∆G where G x is the free energy for ion X, quantifies the thermodynamic preference for K + versus Na + in a medium, e.g., a solvent, a membrane protein, or an ionophore. Given additionally the free energy cost of replacing K + with Na + in bulk liquid water, ∆G , the selectivity of the medium relative to liquid water for K + over Na + can be given in terms of the selectivity free energy, This is simply a comparison of transfer free energies. If G > 0, then the medium favors K + , if G < 0, the medium favors Na + , and if G = 0, then the medium exhibits no preference for either ion. In a bulk liquid setting, the experimental free energy differences ∆∆G K Na → are unambiguous. For example, Fig. 2 displays the selectivity free energy for a variety of organic solvents relative to water as determined from experiments. As noted previously (Varma and Rempe, 2007;Asthagiri et al., 2010), some solvents display positive (K + ) selectivity and some display negative (Na + ) selectivity. The organic solvents composed of molecular analogues of the protein backbone, formamide and N-methylacetamide (NMA), are CNG channels share the same basic architecture of the K + channel selectivity filter (Fig. 1,D and E). Despite their similarity, the CNG channels exhibit no K + /Na + selectivity (Alam and Jiang, 2009). How do the subtle differences between the CNG and K + channels lead to the inability of CNG channels to select for K + ? In contrast, the x-ray structure of a bacterial toxin (Dobler, 1981), valinomycin, bears little resemblance to the selectivity filters of K + and CNG channels (Fig. 1 F). Nevertheless, valinomycin binds K + ions selectively over Na + ions and transports them, one by one, across the membrane. Furthermore, experiments show that the selectivity of valinomycin is sensitive to solvent properties. What explains the behavior of this ionophore?
Given the complexity in structure and the variability in selectivity of these biological molecules, it helps to focus on model systems where specific design principles may be isolated and their ramifications for selectivity directly tested. By determining ion-binding free energies from carefully specified structural and chemical elements, we can seek clear physical ideas to guide our understanding of design principles that give rise to macroscopic behaviors of ion transport molecules. The study of simplified binding site models, with implicit or explicit inclusion of environmental effects, provides a natural route for learning about the underpinnings of selective ion transport, even without proceeding to the selectivity of any specific transport case.
Whole molecule approaches such as electrophysiological or thermodynamic measurements, or computer simulations of entire channels, characterize in situ selectivity. These can also provide tests of mechanistic hypotheses, as illustrated with the experiments by Valiyaveetil et al. (2006) and molecular simulations by Fowler et al. (2008). Representative binding modes of Na + and K + ions in structural motifs of K + -selective membrane transport molecules. (Note that only two units from the tetrameric selectivity filters are shown for clarity.) (A) The selectivity filter of KcsA  adopts different configurations under conditions of high and (B) low K + concentrations, presenting K + with different sets of binding modes. (C) Under rare conditions when Na + binds to the KcsA filter, Na + prefers a binding site different from K + (Nimigean and Miller, 2002;Shrivastava et al., 2002;Lockless et al., 2007;Thompson et al., 2009). (D) The bacterial NaK channel, which belongs to the family of CNG channels, has a selectivity filter architecture similar to KcsA, but is only weakly selective for K + . Initially, low temperature x-ray data suggested binding modes for K + that are identical to Na + . (E) Newer higher resolution crystallographic studies show more variety in Na + binding, attributing electron density at the S3 site to competitive binding of a contaminant (orange) with Na + , with other Na + -binding sites between planes of carbonyl or hydroxyl oxygens (Alam and Jiang, 2009). (F) In comparison to KcsA and the NaK channel, the K + -selective bacterial toxin molecule (Dobler, 1981), valinomycin, binds K + differently, using six (or fewer) carbonyl oxygens, instead of eight (or fewer) as in KcsA and NaK. When reconciling information from approaches based on simplified models or whole channels, it is important to keep in mind that the optimal binding sites of K + and Na + could be different. For example, in the K + and NaK channels (Fig. 1), the K + ion finds a position to organize eight ligands in the selectivity filters, whereas a Na + ion optimally uses fewer ligands. Being distinct from each other, the optimal binding configuration for K + can be a barrier configuration for Na + . The K + -binding regionsbarriers for Na + -can be part of the larger explanation for the low permeability of Na + (Bezanilla and Armstrong, 1972). Overlooking this important point can lead to discrepancies (Bucher and Rothlisberger, 2010) in correlating coordination structure and selectivity.
Here, we will analyze the design principles that have emerged to explain K + /Na + selectivity. We start with a host-guest steric model (also known as "snug fit"), proposed decades before a K + channel was crystallized, that has been used to explain the selectivity of both K + -selective ionophores and K + channels.
A host-guest steric model The first x-ray structure of a strongly selective K + channel, KcsA (Doyle et al., 1998), produced a tangible picture of the molecular bases of K + channel activity. Initial inspection of the structure supported an earlier hypothesis particularly interesting. These coordinate K + /Na + with carbonyl oxygen atoms and provide Na + selectivity. Also providing Na + selectivity are the liquids that coordinate these ions with hydroxyl oxygen atoms, including methanol, ethanol, and propanol.
These simple observations alone tell us something important about K + channels and valinomycin, which bind K + selectively. Despite coordinating the cation (K + or Na + ) with carbonyl oxygen atoms ( Fig. 1), the binding sites in these biological molecules, and the nonselective CNG channels, do not behave like liquid formamide or NMA. If they did, they would be Na + selective. A model of K + /Na + selectivity must, therefore, explain the design differences between the K + -selective sites of these molecules and the liquid states of their chemical analogues.
In contrast to liquid solutions, an ion-binding site in a channel is specified relative to the surrounding protein matrix and defined by the experimental crystal structure. Of course, the ion and the protein matrix, or ionophore matrix, jointly establish the binding configuration. Thermal motion, apparent in molecular dynamics simulation, blurs ideal geometrical consideration of a binding site. The inherent structure concept (Stillinger and Weber, 1982) eliminates that thermal blurring by cataloging basins of optimized configurations that are visited by the thermal motion. Inherent structure analysis has been widely implemented over recent decades (Rao and Karplus, 2010). Quasi-chemical theory, which is relevant to the ion selectivities discussed here (Asthagiri et al., 2010), follows directly from inherent structure concepts (Hummer et al., 1997). If the system presents more than one binding mode, perhaps including multiion binding, these additional complications can be treated by calculating a set of structurally constrained free energies (Rogers and Rempe, 2011).
Free energy differences, determined experimentally, are not as clear for binding sites as in liquids. For example, some experimental methods (Neyton and Miller, 1988;Lockless et al., 2007) probe the equilibrium thermodynamics of ion binding in K + channels. Nonequilibrium electrophysiological techniques (Latorre and Miller, 1983;LeMasurier et al., 2001) might also be used to infer effective K + -and Na + -binding constants in K + channel filters. These inferences rely on mechanistic assumptions about K + /Na + binding in the filter-assumptions about the conformation of the filter, the number of binding sites and their characteristics, and about what other components such as water and additional ions might be doing. The determination of driving forces for selective ion binding, G, is simplest in theoretical and computational settings where it is clear exactly which binding site, conformation, and ion/water occupancy is being considered. In these cases, simple theoretical models and hypotheses pertaining to the design of K + -selective complexes may be unambiguously evaluated in terms of G.
bulk , for different organic solvents with respect to bulk liquid water (Cox and Parker, 1973;Marcus, 1983;Schmid et al., 2000;Yu et al., 2010a). For these estimates, , is the partial molar Gibbs free energy for ion X in the specified medium relative to an ideal standard state, so that G vanishes when ion-medium interactions vanish.
on September 21, 2017 jgp.rupress.org Downloaded from of valinomycin's intramolecular hydrogen bonds, resulting in a selectivity that depends on solvent hydrogenbonding strength. Although this solvent-dependent selectivity of valinomycin was observed experimentally decades ago (Eisenman et al., 1991), an explanation for the role of the surrounding environment arose only recently .
Other selectivity concepts matured alongside the steric mechanism for K + channel selectivity (Bezanilla and Armstrong, 1972). We focus our discussion on the ideas developed by Eisenman and coworkers in the context of K + /Na + selectivity.
The conventional field strength model Studies by Eisenman on the effect of ligand chemistry, or field strength, on the ion selectivity of functionalized surfaces (Eisenman and Horn, 1983) led to two important conclusions: (1) "the primary physical variable controlling equilibrium cationic specificity is the field strength of the anion (i.e., the electronegative coordinating ligand)," and (2) "the primary factor controlling the magnitude of selectivity among ions at a given field strength is the amount of water admitted into the vicinity of the site." If this conventional ligand field strength model is applied to a K + /Na + -selective binding site, then decreasing the field strength of the ligand should increase selectivity for the larger K + ion. Quantum mechanical calculations confirmed this trend for a hypothetical ion-ligand system in which substitution of ion-coordinating formamide by water (a carbonyl-to-water substitution) increases K + /Na + selectivity (Varma and Rempe, 2007). The same behavior also occurred upon varying the ligand field strength in simulations (Fig. 3 A) of a simplified binding site of eight nonpolarizable linear ligands confined within a sphere of radius 3.5 Å around a central ion (Noskov et al., 2004;Thomas et al., 2007;. Although these results illustrate the promise of ligand field strength arguments in explaining K + selectivity, it is important to remember that a successful model should explain the design differences between K + -selective binding sites of molecules and the liquid states of their chemical analogues (Fig. 2).

A revised field strength model
In computational investigations of KcsA, one study concluded that "the carbonyl groups coordinating the ion in the narrow pore are indeed very dynamic ('liquid-like') and that their intrinsic electrostatic properties control ion selectivity" (Noskov et al., 2004). By bringing ligand flexibility into consideration with a ligand field strength model, extended to include not only ion-ligand but also ligand-ligand interactions, this work demonstrated that a thermally fluctuating binding site could be selective.
Although not a complete solution, the focus on simplified binding site models suggested a path for understanding the elusive selectivity problem. that the K + /Na + selectivity in K + channels might be attributable to a steric mechanism (Bezanilla and Armstrong, 1972). The specific positions of the K + ions in the selectivity filter indicated that K + ions preferred to occupy positions coordinated by eight carbonyl oxygen atoms ( Fig. 1 A). Because of the relatively stronger van der Waals interactions of K + compared with Na + at this distance, a selective environment might result if the oxygen atoms of these sites maintained a cavity size that snugly fit the radius of K + but did not collapse onto the smaller Na + ion, as the central assumption to the steric model demands. Thus, a binding site designed to maintain a K + -specific cavity size would hypothetically provide the necessary energetic barrier to Na + passage across the membrane.
Additional crystal diffraction studies produced KcsA structures that contained a K + ion trapped in the waterfilled central cavity of the structure, which also coordinated a cage of exactly eight water oxygen atoms. This finding reinforced the notion that the eightfold carbonyllined sites in the filter served as surrogate hydration shells for K + , mimicking the coordination geometry of K + in the aqueous phase, and thus favoring its fast selective movement Zhou et al., 2001). This view of ion partitioning, however, eventually conflicted with new theoretical and experimental data (Varma and Rempe, 2006) showing that the most probable coordination number of both Na + and K + in liquid water is fewer than eight, and that coordination of Na + or K + with eight water molecules is unlikely in the aqueous phase (Bostick and Brooks, 2007;Varma and Rempe, 2007).
Molecular dynamics simulations of the whole channel embedded in a surrounding membrane (Guidoni et al., 1999;Shrivastava and Sansom, 2000;Bernèche and Roux, 2001) also demonstrated that fluctuations of the selectivity filter backbone atoms are larger than the size difference in Pauling radii between Na + and K + ions, which appeared to conflict with the host-guest steric mechanism of selectivity. Although this cavity-size argument does not satisfactorily explain K + /Na + selectivity in the eightfold binding sites of K + channels, quantum mechanical calculations subsequently demonstrated how this mechanism is operative in the bacterial ionophore molecule, valinomycin .
Valinomycin provides only six carbonyl ligands for ion coordination (Fig. 1 F). Experimental and simulation studies concur that the ionophore achieves K + /Na + selectivity by physically resisting collapse onto the smaller Na + ion. If the ligands are given complete freedom in simplified models, binding sites composed of six carbonyl ligands do not bind K + selectively. Quantum mechanical studies of the whole molecule showed that geometric constraints on the carbonyl ligands in valinomycin are enforced through a combination of intramolecular hydrogen bonds and other structural features, including a specific ring size and the spacing between connected ligands. Competition with solvent can decrease the stability on September 21, 2017 jgp.rupress.org Downloaded from with water molecules modifies the magnitude of K + selectivity but does not abolish or reverse it. Reductions observed in the selectivity free energy were accounted for by differences in molecular geometry (e.g., size, shape, and charge distribution), which are controlled by the full balance of electrostatic, dispersion, and core repulsions in a force field model, rather than simply molecular dipole moment from electrostatic charge-charge interactions . Although the disagreement in Fig. 3 B has not been resolved, more recent simulations showed that K + selectivity persisted in water-based binding site models, even those composed purely of water Brooks, 2007, 2010;Varma and Rempe, 2007;.
Ligand composition clearly plays an important role in modulating selectivity (see Figs. 2 and 3), but a focus on chemical forces as a primary design principle does not account for differences in K + /Na + selectivity between systems possessing comparable ligand composition and/ or field strength. Although the conventional field strength model anticipates the effect of ligand type on ion binding, and the revised field strength model predicts K + selectivity in a fluctuating binding site of strong dipolar ligands, several questions remain. Specifically, why do binding sites in K + channels composed of carbonyl ligands select for K + , but liquids of carbonyl analogues, such as NMA and formamide, inherently select for Na + (see Fig. 2 and Table I)? Furthermore, why do the chemically identical binding sites (i.e., sites S1, S2, and S3) of the K + channel selectivity filter exhibit varying degrees of selectivity noted in whole channel simulations Separating out electrostatic interactions in computations can be problematic (Onsager, 1939), however, because classical electrostatic interactions between charges are balanced non-uniquely by dispersion and core repulsion in force fields (MacKerell, 2005). Nevertheless, observing a correlation between strong ligand-ligand electrostatic interactions and high selectivity might suggest an enhancement of the K + selectivity in binding sites composed of strongly dipolar carbonyls compared with the case of lower-dipole water ligands available in solution. Because this trend would contrast with the simple conventional ligand field strength model (see Figs. 2 and 3 A), we refer to this description as a revised field strength model.
A test of this revised field strength explanation of selectivity for weakly selective CNG (NaK) channels used free energy perturbation simulations Roux, 2006, 2007) in which eight linear carbonyl moieties around an ion were progressively replaced with molecules of a weaker ligand field strength, water. Those substitutions resulted in a loss and reversal of K + /Na + selectivity ( Fig. 3 B), which supported the view that increasing the ligand dipole moment or field strength (water to carbonyl) should enhance the free energy of K + selectivity, G.
Although attractive, this interpretation is difficult to reconcile with previously observed trends in simulations (Fig. 3 A) and experiments (Eisenman and Horn, 1983). Additionally, calculations (Fig. 3 B) using an identical methodology and force field to the initial works found that progressive substitution of linear carbonyl ligands In contrast to the trend in A, recent work predicted a systematic loss of selectivity for each water molecule that replaces a carbonyl group (1.8 kcal/mol per water using the CHARMM force field) Roux, 2006, 2007). The red line illustrates this trend toward Na + selectivity, which supports the revised field strength model. Data from subsequent calculations ), using the same force field and either a Lennard-Jones (LJ, black solid lines/circles) or a half-harmonic (black dashed lines and open circles) confining potential, do not eliminate the K + selectivity. The constrained cavity size of valinomycin effected through intra-molecular hydrogen bonds, illustrates one particular form of architectural control . Constraints imposing over-coordination with respect to liquid solutions provide yet another example of architectural control (Bostick and Brooks, 2007;Varma and Rempe, 2007).
Simulations of binding site models illustrate the consequences of constraints on coordination (Fig. 4). Coordination by five or six carbonyl oxygen atoms typically provides the weak selectivity observed in bulk liquids (Fig. 2), whereas strong K + selectivity is associated with seven and eight coordinate complexes that resemble K + channel architectures (Figs. 1 A and 4 B). In the absence of external forces imposing over-coordination, carbonyl ligands naturally assume a lower coordination number around each ion and K + /Na + selectivity is lost (Fig. 4 A). Quantum mechanical studies of a model binding site composed of diglycine molecules show the same trends with Na + binding (Varma and Rempe, 2007). These results support the argument that an environment that achieves high coordination numbers (Noskov et al., 2004), and why do NaK channels show reduced selectivity despite sharing two chemically identical binding sites with strongly selective K + channels? The parsimonious conclusion is that the focus on ligand type is incomplete; the role of the surrounding matrix must be incorporated into a model that explains selectivity.

A primary role for the surrounding matrix
Recent studies shed light on issues concerning variable selectivity in systems possessing comparable ligand composition or field strength. Using physiologically motivated ion-binding site models, several investigations suggested that external forces derived from the surrounding environment and exerted on the ion-ligated complex play a primary role in determining selectivity Brooks, 2007, 2009;Thomas et al., 2007;Rempe, 2007, 2008;. The most unambiguous demonstrations of these principles explicitly include electronic polarization (Guidoni and Carloni, 2002;Bostick and Brooks, 2010;Varma and Rempe, 2010) because the non-transferability of pairwise-additive force fields can lead to discrepancies (see Table I).
Returning to Fig. 2, we notice that ionic selectivities in bulk liquids, particularly for the carbonyl ligands, provide a control that illustrates K + /Na + selectivity in the absence of macromolecular structure. Because selectivity differs between organic solvents and K + channels, it must be concluded that the environment surrounding an ion-binding site is also necessary for explaining the selectivity of K + channels. The challenge then centers on understanding physiological mechanisms that modulate binding site properties, and thus modulate selective ion binding.
Interactions between the ion-ligand complex and the remainder of the protein, membrane surfaces, other ions, or the solution can be modeled as external fields acting on the ion and its coordinating ligands. From the quasi-chemical development (Asthagiri et al., 2010) and other statistical mechanical theories , it is clear that these interactions need not be pairdecomposable. The effects of these interactions have been referred to in the literature as "topological" or "architectural" control (Bostick and Brooks, 2007, 2010Thomas et al., 2007;Varma and Rempe, 2007;). Those interactions typically affect both the selectivity free energy and the structural attributes apparent in Fig. 1, such as the coordination number the ion-ligand complex adopts, the orientations of its ligands, or its deformability. In fact, recent crystallographic results confirm that the structures of channel-binding sites are not fixed; instead, the KcsA selectivity filter can take on radically different configurations depending on the surroundings (Cordero-Morales et al., 2006). 0.3 Na + K + /Na + selectivity of bulk liquid formamide and NMA as determined by experimental work, pairwise-additive (nonpolarizable) force field calculations, polarizable force field (AMOEBA and Drude) calculations, and QM/MM (hybrid quantum mechanical and molecular mechanical) calculations. Note that the pairwise-additive force field calculations display positive (K + ) selectivity, whereas all other determinations display negative (Na + ) selectivity. a Marcus, 1983. b Cox andParker, 1973, andSchmid et al., 2000. c Grossfield et al., 2003. d Yu et al., 2010a. e Noskov et al., 2004   on September 21, 2017 jgp.rupress.org

Downloaded from
Here, ∆U K Na field → is the internal energy contribution from an imposed external field, ∆U K Na IL → , ∆U K Na LL → , and ∆U K Na intra → are contributions from ion-ligand interactions, ligand-ligand interactions, and intramolecular interactions, respectively, and − → T S ∆ K Na is the entropic contribution. Decomposition of the selectivity free energy in this form shows several large contributions (Fig. 4 C), with the unfavorable ligand-ligand interactions ( ∆U K Na LL → ) providing the largest contribution favoring K + selectivity.
These results might be interpreted as suggesting that the ligand identity controls K + /Na + selectivity through local interactions, as proposed by the revised field strength model. It is, however, important to note that each contribution in Eq. 2, including ligand-ligand interactions ∆U K Na LL → , depends on the constraints imposed on the ion-ligated complex-in the case of Fig. 4, a constraint on coordination number Bostick and Brooks, 2010). In fact, without requiring the larger results in K + /Na + selectivity, whereas an environment that permits relaxed coordination numbers produces lower selectivity. Note that this does not suggest that selectivity, in general, is solely controlled by coordination number: other types of constraints on an ion-bound complex can also play a role in selectivity, and ligand chemistry remains an important modulatory factor.
To probe the energetic details giving rise to selectivity, the K + selectivity free energy can be parsed in several ways (see, for example, Asthagiri et al., 2006;Bostick and Brooks, 2007;Varma and Rempe, 2007;Rogers and Rempe, 2011). One way is to decompose the selectivity free energy into contributions from different types of interactions in the force field. Using cases considered in Fig. 4 as an example, this decomposition can be written as (Bostick and Brooks, 2010), Figure 4. Results from binding site models demonstrating the effect of an external field, in the form of an ion coordination number constraint, on K + /Na + selectivity. Calculations were performed using the AMOEBA polarizable force field (Bostick and Brooks, 2010).
(A) G versus the number of included molecules, N I , in gas-phase clusters around K + and Na + in the absence of a constraint on coordination. As N I increases, the observed selectivity approaches values expected for bulk liquids (for water, G ≈ 0; for formamide/ NMA, G < 0). Note that the number of included molecules, N I , is not necessarily equal to the number of coordinating molecules, N C , because of the absence of a coordination constraint. (B) G versus the number of molecules, N C , directly coordinating K + and Na + . In agreement with quantum mechanical calculations Rempe, 2007, 2008), G is larger in the water-based models (blue) than in the carbonyl-based models (black and red). Because of the presence of an external field (half-harmonic confinement) that imposes a specific coordination number, K + selectivity is observed for seven or more ligands. In the models that coordinate K + or Na + with carbonyl ligands, K + selectivity is determined by the external field rather than the ligand identity (NMA, formamide, and water) because the binding site models are Na + selective in the absence of the constraint on coordination number, N C . (C) Contributions from specified individual components of the selectivity free energy (Eq. 2) in models composed of eight formamide molecules: ligand-ligand interactions, ∆U K Na LL → , ion-ligand interactions, ∆U K Na IL → , intramolecular interactions, ∆U K Na intra → , and entropy − → T S ∆ K Na . In the case considered here, the contribution from the external field, ∆U K Na field → ≈ 0, is negligible. In the absence of the field (left), the components yield net Na + selectivity (G < 0). When a field enforces eightfold coordination, the distribution of these individual components changes, producing net K + selectivity (G > 0). Thus, the redistribution of the individual components, and therefore the net K + selectivity, is an effect of the applied external field.
on September 21, 2017 jgp.rupress.org Downloaded from K + channels. Hydrogen bond-donating groups in weakly selective K + channels can reduce the availability of ligands to bind ions, resulting in reduced selectivities. In contrast, the presence of proximal hydrogen bondaccepting groups can enhance K + selectivity by constraining optimal filter orientation, as predicted earlier Rempe, 2007, 2008) and shown recently via electrophysiology and crystallography (Cheng et al., 2011). This provides one mechanism for the potassium channel family to exhibit a range of selectivities, despite sharing a common filter architecture (Varma and Rempe, 2007).
Alternatively, the observed diversity in selectivity exhibited by K + channels could also arise from differences in the flexibilities of the ion-ligand complexes, which are influenced by the external environment as well (Bostick and Brooks, 2007;Varma and Rempe, 2007). Observations of large shifts in the crystal structures of KcsA selectivity filters with changes to the surrounding protein matrix lend support to this proposition (Cordero-Morales et al., 2006;Cheng et al., 2011).
Yet another environmental factor that can modulate a site's selectivity is the availability of water molecules to a binding site. Because water molecules can serve as both hydrogen bond donors or acceptors, their presence provides external interactions that can reduce selectivity by either competing with ion-coordinating ligands for the ion, or by competing with the ion for ion-coordinating ligands (Bostick and Brooks, 2007;Varma and Rempe, 2007). This provides an explanation for the historical observations of the modulating effect of water exposure on ion selectivity (Eisenman and Horn, 1983).
In alignment with earlier predictions (Bostick and Brooks, 2007;Varma and Rempe, 2007), computer simulations of membrane-embedded whole channels (Fowler et al., 2008) indicated that the lack of K + selectivity in CNG channels is accounted for by the combined effect of the enhanced exposure of binding sites to water and the enhanced flexibilities of the ion-binding sites in their selectivity filters, both of which foster lower coordination numbers for K + and Na + .

Summary and conclusions
The advent of atomic resolution structures of ion-selective channels made possible the transition from black box models to molecular descriptions for considering the design principles underlying selectivity in a flexible selectivity filter. Simple theoretical models applied to interpret the complex behaviors observed in whole molecule experiments and molecular simulations led to suggestions that ligand type is of primary importance in determining selectivity. As discussed in this Perspective, this view is incomplete. Models that focus on chemical forces as a primary design principle do not explain why carbonyl ligands produce K + selectivity in eightfold K + channel-binding sites, but yield Na + selectivity in liquid analogues (Fig. 2). Additionally, these models have not coordination number, these various contributions should combine to yield net Na + selectivity, consistent with the experimental liquid solution results of Fig. 2. Fig. 4 C explicitly demonstrates the shift in contributions caused by imposing a constraint on coordination number for a test case of eight formamide molecules interacting with K + and Na + .
Acknowledging the influence of the environment in ion binding, the question remains: how are those interactions involved in K + channels? Protein composition and water occupancy help determine properties of the environment surrounding channel-binding sites. Significant environmental characteristics that could modulate these binding sites include, for example, steric and electrostatic interactions that define the flexibility of the local ion-ligand complex, and the distribution of functional groups that can hydrogen bond with the ioncoordinating ligands Rempe, 2007, 2008).
During formation of an ionic binding site complex, ligands are extracted from the exterior environment and rearranged. Depending on characteristics of the environment, this ligand extraction step can be associated with high free energy penalties, even when the overall process of ion solvation is favorable. Computational studies Rempe, 2007, 2008) demonstrate that alteration of these penalties affects the distribution of ligands in the ion complex. For example, reducing extraction penalties increases the stabilities of high coordination numbers, thus bringing down the energetic cost of transferring ions from water into overcoordinated binding sites. A significant reduction in this penalty for binding site formation can even drive up preferred coordination numbers. In contrast, an increased penalty can lower ion coordination preferences, as observed in highly concentrated salt solutions.
In biomolecules that bind ions in over-coordinated states, the free energy penalties for extracting ligands can be lowered by reducing direct favorable interactions of the ligands with atoms other than the ions themselves. In the case of highly selective K + channels, a bioinformatics and molecular simulation analysis  suggests that lowered ligand-extraction penalties that stabilize high coordinations are accomplished by a lowered density of free hydrogen bond donors near the carbonyl ligands of the binding sites, and by a constraint on the motions of selectivity filters and their neighboring environments.
Furthermore, a sequence alignment of K + channels (Shealy et al., 2003) indicates specific chemical differences in hydrogen-bonding ability between the ionbinding neighborhoods of strongly and weakly selective channels. In weakly selective channels, the immediate neighborhoods of the ion-coordinating carbonyl ligands contain specific amino acids (hydrogen bond donors) that can form hydrogen bonds with those ligands. These donors are typically absent in strongly selective on September 21, 2017 jgp.rupress.org Downloaded from explained why chemically identical binding sites in the strongly selective KcsA channel show significantly different selectivities in simulations of whole channels. The same question applies to NaK channels, which lack selectivity despite sharing two chemically identical binding sites with strongly selective K + channels. We are left with the conclusion that ion selectivity requires consideration of both ligand characteristics and the forces influencing binding site composition and structure.
Recent simulations support the view that interactions with the more distant environment of the membrane transport molecule can modify properties of the binding site and influence selective binding. In the case of K + channels, restricted carbonyl motion or a decrease in their availability to the protein environment can drive up ion coordination numbers in the selectivity filter, leading to K + selectivity. A snug fit as in valinomycin can be enforced by an environment that stabilizes intramolecular hydrogen bonding, achieving selectivity without over-coordination through a constraint on cavity size. The importance of environmental controls on binding sites has since been recognized in several recent works (see, for example, Fowler et al., 2008;Miloshevsky and Jordan, 2008;Vora et al., 2008;Dixit et al., 2009;Dudev and Lim, 2009;Yu et al., , 2010bRoux, 2010;Rogers and Rempe, 2011). This Perspective provides a foundation necessary for understanding the more complex behavior of selective ion transport through membranes.
This Perspectives series includes articles by Andersen, Alam and Jiang, Nimigean and Allen, Roux et al., and Dixit and Asthagiri. D.M. Rogers and S.B. Rempe acknowledge funding from the National Institutes of Health (NIH) through the NIH Road Map for Medical Research and the Sandia Laboratory-Directed Research and Development program. Sandia National Laboratories is a multi-program laboratory operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.