High-throughput characterization of photocrosslinker-bearing ion channel variants to map residues critical for function and pharmacology

Incorporation of noncanonical amino acids (ncAAs) can endow proteins with novel functionalities, such as crosslinking or fluorescence. In ion channels, the function of these variants can be studied with great precision using standard electrophysiology, but this approach is typically labor intensive and low throughput. Here, we establish a high-throughput protocol to conduct functional and pharmacological investigations of ncAA-containing human acid-sensing ion channel 1a (hASIC1a) variants in transiently transfected mammalian cells. We introduce 3 different photocrosslinking ncAAs into 103 positions and assess the function of the resulting 309 variants with automated patch clamp (APC). We demonstrate that the approach is efficient and versatile, as it is amenable to assessing even complex pharmacological modulation by peptides. The data show that the acidic pocket is a major determinant for current decay, and live-cell crosslinking provides insight into the hASIC1a–psalmotoxin 1 (PcTx1) interaction. Further, we provide evidence that the protocol can be applied to other ion channels, such as P2X2 and GluA2 receptors. We therefore anticipate the approach to enable future APC-based studies of ncAA-containing ion channels in mammalian cells.

This study describes a method to rapidly screen hundreds of ion channel variants containing non-canonical amino acids. A proof-of-principle introducing photocrosslinking non-canonical amino acids into the human ion channel hASIC1a shows how this approach can provide insights into function and pharmacology.     

Putative Binding Sites,and Pathways to Them,for Amidine and Guanidine Current Inhibitors on Acid‐Sensing Ion Channels (ASIC). A Theoretical Approach with hASIC1a Homology Model

Central inhibition of the acid‐sensing hASIC1a channel, acting upstream of the opiate system, might serve to treat any type of pain, avoiding the unwanted addiction problems of the opioid drugs. To this end, inhibition of hASIC1a channel by PcTx1, a peptide from the Trinidad chevron tarantula, is under development. New inhibitors of the hASIC1a channel are also being sought, in the hope of further modulating the activity, from which antiplasmodial amidine and guanidine phenyl drugs have emerged as promising candidates. However, how such current inhibition takes place remains obscure from the molecular point of view, hindering any further progress in developing drugs. Therefore, the nature of the binding sites, and how they are reached by the amidine‐guanidine drugs, was investigated here via automated docking and molecular dynamics with hASIC1a homology models. This study has revealed that this ion channel is rich in binding sites, and that flexible drugs, such as nafamostat, may penetrate it in a snake‐like elongated conformation. Then, crawling like a snake through temporary holes in the protein, nafamostat either simply flips, or changes to a high‐energy folded conformation to become adapted to the shape of the binding site.     

Functional Importance of the Carboxyl Tail Cysteine Residues in the Human D1 Dopamine Receptor

Abstract: To assess the importance of the cysteine residues Cys³⁴⁷ and Cys³⁵¹ in the carboxylic tail in the human D₁ dopamine receptor, seven mutant receptors were constructed by PCR. The pharmacological and functional properties of the wild-type and mutant receptors were assessed following transient expression in COS-7 cells. Affinities for [³H]SCH 23390 of mutant S347 (Cys³⁴⁷→ Gly), T348 (Tyr³⁴⁸→ stop), S351 (Cys³⁵¹→ Gly), T351 (Cys³⁵¹→ stop), T352 (Pro³⁵²→ stop), and S347/S351 (Cys³⁴⁷→ Gly and Cys³⁵¹→ Gly) were similar to that of wild-type receptor, whereas the expression levels were reduced up to 80%. The potency of dopaminergic antagonists for these mutant receptors was very similar to that of the wild-type receptor. However, mutant T347 (Cys³⁴⁷→ stop) showed a 15–25-fold reduced affinity for the antagonists SCH 23390, (+)-butaclamol, and cis-flupentixol, thus not allowing radioligand analysis. Wild-type and mutant receptors responded dose-dependently with similar potency to dopamine and SKF 38393 with an increased adenylyl cyclase activity. However, mutant receptors with the Cys³⁴⁷ residue changed or removed displayed a diminished ability to activate adenylyl cyclase. Dopamine preexposure desensitized wild-type and mutant S351 receptors. However, mutant receptors with Cys³⁴⁷ replaced or the distal part of the carboxyl tail removed were unable to desensitize. Thus, Cys³⁴⁷ in the cytoplasmic tail of the human D₁ dopamine receptor is important for the receptor in maintaining the conformation for antagonist binding, to play a crucial role in activation of adenylyl cyclase, and to be essential for agonist-induced desensitization.     

Putative binding sites, and pathways to them, for amidine and guanidine current inhibitors on acid-sensing ion channels (ASIC). A theoretical approach with hASIC1a homology model

Central inhibition of the acid-sensing hASIC1a channel, acting upstream of the opiate system, might serve to treat any type of pain, avoiding the unwanted addiction problems of the opioid drugs. To this end, inhibition of hASIC1a channel by PcTx1, a peptide from the Trinidad chevron tarantula, is under development. New inhibitors of the hASIC1a channel are also being sought, in the hope of further modulating the activity, from which antiplasmodial amidine and guanidine phenyl drugs have emerged as promising candidates. However, how such current inhibition takes place remains obscure from the molecular point of view, hindering any further progress in developing drugs. Therefore, the nature of the binding sites, and how they are reached by the amidine-guanidine drugs, was investigated here via automated docking and molecular dynamics with hASIC1a homology models. This study has revealed that this ion channel is rich in binding sites, and that flexible drugs, such as nafamostat, may penetrate it in a snake-like elongated conformation. Then, crawling like a snake through temporary holes in the protein, nafamostat either simply flips, or changes to a high-energy folded conformation to become adapted to the shape of the binding site.     

Identification of cystic fibrosis transmembrane conductance regulator channel-lining residues in and flanking the M6 membrane-spanning segment.

The cystic fibrosis transmembrane conductance regulator (CFTR) forms a chloride channel that is regulated by phosphorylation and ATP binding. Work by others suggested that some residues in the sixth transmembrane segment (M6) might be exposed in the channel and play a role in ion conduction and selectivity. To identify the residues in M6 that are exposed in the channel and the secondary structure of M6, we used the substituted cysteine accessibility method. We mutated to cysteine, one at a time, 24 consecutive residues in and flanking the M6 segment and expressed these mutants in Xenopus oocytes. We determined the accessibility of the engineered cysteines to charged, lipophobic, sulfhydryl-specific methanethiosulfonate (MTS) reagents applied extracellularly. The cysteines substituted for Ile331, Leu333, Arg334, Lys335, Phe337, Ser341, Ile344, Arg347, Thr351, Arg352, and Gln353 reacted with the MTS reagents, and we infer that they are exposed on the water-accessible surface of the protein. From the pattern of the exposed residues we infer that the secondary structure of the M6 segment includes both alpha-helical and extended regions. The diameter of the channel from the extracellular end to the level of Gln353 must be at least 6 A to allow the MTS reagents to reach these residues.     

Interaction of nucleotides with Asp(351) and the conserved phosphorylation loop of sarcoplasmic reticulum Ca(2+)-ATPase.

The nucleotide binding properties of mutants with alterations to Asp(351) and four of the other residues in the conserved phosphorylation loop, (351)DKTGTLT(357), of sarcoplasmic reticulum Ca(2+)-ATPase were investigated using an assay based on the 2', 3'-O-(2,4,6-trinitrophenyl)-8-azidoadenosine triphosphate (TNP-8N(3)-ATP) photolabeling of Lys(492) and competition with ATP. In selected cases where the competition assay showed extremely high affinity, ATP binding was also measured by a direct filtration assay. At pH 8.5 in the absence of Ca(2+), mutations removing the negative charge of Asp(351) (D351N, D351A, and D351T) produced pumps that bound MgTNP-8N(3)-ATP and MgATP with affinities 20-156-fold higher than wild type (K(D) as low as 0.006 microM), whereas the affinity of mutant D351E was comparable with wild type. Mutations K352R, K352Q, T355A, and T357A lowered the affinity for MgATP and MgTNP-8N(3)-ATP 2-1000- and 1-6-fold, respectively, and mutation L356T completely prevented photolabeling of Lys(492). In the absence of Ca(2+), mutants D351N and D351A exhibited the highest nucleotide affinities in the presence of Mg(2+) and at alkaline pH (E1 state). The affinity of mutant D351A for MgATP was extraordinarily high in the presence of Ca(2+) (K(D) = 0.001 microM), suggesting a transition state like configuration at the active site under these conditions. The mutants with reduced ATP affinity, as well as mutants D351N and D351A, exhibited reduced or zero CrATP-induced Ca(2+) occlusion due to defective CrATP binding.     

Recombinant production and solution structure of PcTx1, the specific peptide inhibitor of ASIC1a proton-gated cation channels

Acid-sensing ion channels (ASICs) are thought to be important ion channels, particularly for the perception of pain. Some of them may also contribute to synaptic plasticity, learning, and memory. Psalmotoxin 1 (PcTx1), the first potent and specific blocker of the ASIC1a proton-sensing channel, has been successfully expressed in the Drosophila melanogaster S2 cell recombinant expression system used here for the first time to produce a spider toxin. The recombinant toxin was identical in all respects to the native peptide, and its three-dimensional structure in solution was determined by means of (1)H 2D NMR spectroscopy. Surface characteristics of PcTx1 provide insights on key structural elements involved in the binding of PcTx1 to ASIC1a channels. They appear to be localized in the beta-sheet and the beta-turn linking the strands, as indicated by electrostatic anisotropy calculations, surface charge distribution, and the presence of residues known to be implicated in channel recognition by other inhibitor cystine knot (ICK) toxins.     

Identification of a glutamic acid and an aspartic acid residue essential for catalytic activity of aspergillopepsin II, a non-pepsin type acid proteinase

Aspergillopepsin II from Aspergillus niger var. macrosporus is a non-pepsin type or pepstatin-insensitive acid proteinase. To identify the catalytic residues of the enzyme, all acidic residues that are conserved in the homologous proteinases of family A4 were replaced with Asn, Gln, or Ala using site-directed mutagenesis. The wild-type and mutant pro-enzymes were heterologously expressed in Escherichia coli and refolded in vitro. The wild-type pro-enzyme was shown to be processed into a two-chain active enzyme under acidic conditions. Most of the recombinant mutant pro-enzymes showed significant activity under acidic conditions because of autocatalytic activation except for the D123N, D123A, E219Q, and E219A mutants. The D123A, E219Q, and E219A mutants showed neither enzymatic activity nor autoprocessing activity under acidic conditions. The circular dichroism spectra of the mutant pro- and mature enzymes were essentially the same as those of the wild-type pro- and mature enzyme, respectively, indicating that the mutant pro-enzymes were correctly folded. In addition, two single and one double mutant pro-enzyme, D123E, E219D, and D123E/E219D, did not show enzymatic activity under acidic conditions. Taken together, Glu-219 and Asp-123 are deduced to be the catalytic residues of aspergillopepsin II.     

Acidic residues and a predicted,highly conserved α‐helix are critical for the endonuclease/strand separation functions of bacteriophage λ’s TerL

Complementation, endonuclease, strand separation, and packaging assays using mutant TerL^λ’s, coupled with bioinformatic information and modeling of its endonuclease, identified five residues, D401, E408, D465, E563, and E586, as critical acidic residues of TerL^λ’s endonuclease. Studies of phage and viral TerL nucleases indicate acidic residues participate in metal ion‐binding, part of a two‐ion metal catalysis mechanism, where metal ion A activates a water for DNA backbone hydrolysis. Modeling places D401, D465, and E586 in locations analogous to those of the metal‐binding residues of many phage and viral TerLs. Our work leads to a model of TerL^λ’s endonuclease domain where at least three acidic residues from a ~185 residue segment (D401 to E586) are near each other in the structure, forming the endonuclease catalytic center at cosN, the nicking site. DNA interactions required to bring the rotationally symmetric cosN precisely to the catalytic center are proposed to rely on an ~60 residue region that includes a conserved α‐helix for dimerization. Metal ion A, positioned by TerL^λ’s acidic D401 and E586, would be placed at cosN for water activation, ensuring high accuracy for DNA backbone hydrolysis.     

Golgi Anti-apoptotic Proteins Are Highly Conserved Ion Channels That Affect Apoptosis and Cell Migration

Golgi anti-apoptotic proteins (GAAPs) are multitransmembrane proteins that are expressed in the Golgi apparatus and are able to homo-oligomerize. They are highly conserved throughout eukaryotes and are present in some prokaryotes and orthopoxviruses. Within eukaryotes, GAAPs regulate the Ca²⁺ content of intracellular stores, inhibit apoptosis, and promote cell adhesion and migration. Data presented here demonstrate that purified viral GAAPs (vGAAPs) and human Bax inhibitor 1 form ion channels and that vGAAP from camelpox virus is selective for cations. Mutagenesis of vGAAP, including some residues conserved in the recently solved structure of a related bacterial protein, BsYetJ, altered the conductance (E207Q and D219N) and ion selectivity (E207Q) of the channel. Mutation of residue Glu-207 or -178 reduced the effects of GAAP on cell migration and adhesion without affecting protection from apoptosis. In contrast, mutation of Asp-219 abrogated the anti-apoptotic activity of GAAP but not its effects on cell migration and adhesion. These results demonstrate that GAAPs are ion channels and define residues that contribute to the ion-conducting pore and affect apoptosis, cell adhesion, and migration independently.     

Two Salt Bridges Differentially Contribute to the Maintenance of Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) Channel Function

Previous studies have identified two salt bridges in human CFTR chloride ion channels, Arg³⁵²-Asp⁹⁹³ and Arg³⁴⁷-Asp⁹²⁴, that are required for normal channel function. In the present study, we determined how the two salt bridges cooperate to maintain the open pore architecture of CFTR. Our data suggest that Arg³⁴⁷ not only interacts with Asp⁹²⁴ but also interacts with Asp⁹⁹³. The tripartite interaction Arg³⁴⁷-Asp⁹²⁴-Asp⁹⁹³ mainly contributes to maintaining a stable s2 open subconductance state. The Arg³⁵²-Asp⁹⁹³ salt bridge, in contrast, is involved in stabilizing both the s2 and full (f) open conductance states, with the main contribution being to the f state. The s1 subconductance state does not require either salt bridge. In confirmation of the role of Arg³⁵² and Asp⁹⁹³, channels bearing cysteines at these sites could be latched into a full open state using the bifunctional cross-linker 1,2-ethanediyl bismethanethiosulfonate, but only when applied in the open state. Channels remained latched open even after washout of ATP. The results suggest that these interacting residues contribute differently to stabilizing the open pore in different phases of the gating cycle.     

Interaction of acid-sensing ion channel (ASIC) 1 with the tarantula toxin psalmotoxin 1 is state dependent

Acid-sensing ion channels (ASICs) are Na(+) channels gated by extracellular H(+). Six ASIC subunits that are expressed in neurons have been characterized. The tarantula toxin psalmotoxin 1 has been reported to potently and specifically inhibit homomeric ASIC1a and has been useful to characterize ASICs in neurons. Recently we have shown that psalmotoxin 1 inhibits ASIC1a by increasing its apparent affinity for H(+). However, the mechanism by which PcTx1 increases the apparent H(+) affinity remained unclear. Here we show that PcTx1 also interacts with ASIC1b, a splice variant of ASIC1a. However, PcTx1 does not inhibit ASIC1b but promotes its opening; under slightly acidic conditions, PcTx1 behaves like an agonist for ASIC1b. Our results are most easily explained by binding of PcTx1 with different affinities to different states (closed, open, and desensitized) of the channel. For ASIC1b, PcTx1 binds most tightly to the open state, promoting opening, whereas for ASIC1a, it binds most tightly to the open and the desensitized state, promoting desensitization.     

Mutations at Arginine 352 Alter the Pore Architecture of CFTR

Arginine 352 (R352) in the sixth transmembrane domain of the cystic fibrosis transmembrane conductance regulator (CFTR) previously was reported to form an anion/cation selectivity filter and to provide positive charge in the intracellular vestibule. However, mutations at this site have nonspecific effects, such as inducing susceptibility of endogenous cysteines to chemical modification. We hypothesized that R352 stabilizes channel structure and that charge-destroying mutations at this site disrupt pore architecture, with multiple consequences. We tested the effects of mutations at R352 on conductance, anion selectivity and block by the sulfonylurea drug glipizide, using recordings of wild-type and mutant channels. Charge-altering mutations at R352 destabilized the open state and altered both selectivity and block. In contrast, R352K-CFTR was similar to wild-type. Full conductance state amplitude was similar to that of wild-type CFTR in all mutants except R352E, suggesting that R352 does not itself form an anion coordination site. In an attempt to identify an acidic residue that may interact with R352, we found that permeation properties were similarly affected by charge-reversing mutations at D993. Wild-type-like properties were rescued in R352E/D993R-CFTR, suggesting that R352 and D993 in the wild-type channel may interact to stabilize pore architecture. Finally, R352A-CFTR was sensitive to modification by externally applied MTSEA⁺, while wild-type and R352E/D993R-CFTR were not. These data suggest that R352 plays an important structural role in CFTR, perhaps reflecting its involvement in forming a salt bridge with residue D993.     

Roles of Glu 349 and Asp 352 in membrane insertion and translocation by diphtheria toxin.

Acidic conditions within the endosomal lumen induce the T domain of receptor-bound diphtheria toxin (DT) to insert into the endosomal membrane and mediate translocation of the toxin's catalytic domain to the cytosol. A conformational rearrangement in the toxin occurring near pH5 allows a buried apolar helical hairpin of the native T domain (helices TH8 and TH9) to undergo membrane insertion. If the inserted hairpin spans the bilayer, as hypothesized, then the two acidic residues within the TL5 interhelical loop, Glu 349 and Asp 352, should become exposed at the neutral cytosolic face of the membrane and reionize. To investigate the roles of these residues in toxin action, we characterized mutant toxins in which one or both acidic residues had been replaced with nonionizable ones. Each of two double mutants examined showed a several-fold reduction in cytotoxicity in 24-h Vero cell assays (sixfold for E349A + D352A and fourfold for E349Q + D352N), whereas the individual E349Q and D352N mutations caused smaller reductions in toxicity. The single and double mutations also attenuated the toxin's ability to permeabilize Vero cells to Rb+ at low pH and decreased channel formation by the toxin in artificial planar bilayers. Neither of the double mutations affected the pH-dependence profile of the toxin's conformational rearrangement in solution, as measured by binding of the hydrophobic fluorophore, 2-p-toluidinyl-naphthalene 6-sulfonate. The results demonstrate that, although there is no absolute requirement for an acidic residue within the TL5 loop for toxicity, Glu 349 and Asp 352 do significantly enhance the biological activity of the protein. The data are consistent with a model in which ionization of these residues at the cytosolic face of the endosomal membrane stabilizes the TH8/TH9 hairpin in a transmembrane configuration, thereby facilitating channel formation and translocation of the toxin's catalytic chain.     

Orai1 mutations alter ion permeation and Ca2+-dependent fast inactivation of CRAC channels: evidence for coupling of permeation and gating

Ca²⁺ entry through store-operated Ca²⁺ release-activated Ca²⁺ (CRAC) channels is an essential trigger for lymphocyte activation and proliferation. The recent identification of Orai1 as a key CRAC channel pore subunit paves the way for understanding the molecular basis of Ca²⁺ selectivity, ion permeation, and regulation of CRAC channels. Previous Orai1 mutagenesis studies have indicated that a set of conserved acidic amino acids in trans membrane domains I and III and in the I–II loop (E106, E190, D110, D112, D114) are essential for the CRAC channel's high Ca²⁺ selectivity. To further dissect the contribution of Orai1 domains important for ion permeation and channel gating, we examined the role of these conserved acidic residues on pore geometry, properties of Ca²⁺ block, and channel regulation by Ca²⁺. We find that alteration of the acidic residues lowers Ca²⁺ selectivity and results in striking increases in Cs⁺ permeation. This is likely the result of enlargement of the unusually narrow pore of the CRAC channel, thus relieving steric hindrance for Cs⁺ permeation. Ca²⁺ binding to the selectivity filter appears to be primarily affected by changes in the apparent on-rate, consistent with a rate-limiting barrier for Ca²⁺ binding. Unexpectedly, the mutations diminish Ca²⁺-mediated fast inactivation, a key mode of CRAC channel regulation. The decrease in fast inactivation in the mutant channels correlates with the decrease in Ca²⁺ selectivity, increase in Cs⁺ permeability, and enlargement of the pore. We propose that the structural elements involved in ion permeation overlap with those involved in the gating of CRAC channels.     

Structural role for a conserved region in the CTP synthetase glutamine amide transfer domain.

Site-directed mutations were introduced into a conserved region of the Escherichia coli CTP synthetase glutamine amide transfer domain. The amino acid replacements, valine 349 to serine, glycine 351 to alanine, glycine 352 to proline, and glycine 352 to cysteine, all increased the lability of CTP synthetase. The proline 352 replacement abolished the capacity to form the covalent glutaminyl-cysteine 379 catalytic intermediate, thus preventing glutamine amide transfer function; NH3-dependent CTP synthetase activity was retained. In CTP synthetase (serine 349), both glutamine and NH3-dependent activities were increased approximately 30% relative to that of the wild type. CTP synthetase mutants alanine 351 and cysteine 352 were not overproduced because of apparent instability and proteolytic degradation. We conclude that the conserved region between residues 346 and 355 in the CTP synthetase glutamine amide transfer domain has an important structural role.     

Protons and Psalmotoxin-1 reveal nonproton ligand stimulatory sites in chicken acid-sensing ion channel: Implication for simultaneous modulation in ASICs

Acid-sensing ion channels (ASICs) are proton-sensitive, sodium-selective channels expressed in the nervous system that sense changes in extracellular pH. These ion channels are sensitive to an increasing number of nonproton ligands that include natural venom peptides and guanidine compounds. In the case of chicken ASIC1, the spider toxin Psalmotoxin-1 (PcTx1) activates the channel, resulting in an inward current. Furthermore, a growing class of ligands containing a guanidine group has been identified that stimulate peripheral ASICs (ASIC3), but exert subtle influence on other ASIC subtypes. The effects of the guanidine compounds on cASIC1 have not been the focus of previous study. Here, we investigated the interaction of the guanidine compound 2-guanidine-4-methylquinazoline (GMQ) on cASIC1 proton activation and PcTx1 stimulation. Exposure of expressed cASIC1 to PcTx1 resulted in biphasic currents consisting of a transient peak followed by an irreversible cASIC1 PcTx1 persistent current. This cASIC1 PcTx1 persistent current may be the result of locking the cASIC1 protein into a desensitized transition state. The guanidine compound GMQ increased the apparent affinity of protons on cASIC1 and decreased the half-maximal constant of the cASIC1 steady-state desensitization profile. Furthermore, GMQ stimulated the cASIC1 PcTx1 persistent current in a concentration-dependent manner, which resulted in a non-desensitizing inward current. Our data suggests that GMQ may have multiple sites within cASIC1 and may act as a “molecular wedge” that forces the PcTx1-desensitized ASIC into an open state. Our findings indicate that guanidine compounds, such as GMQ, may alter acid-sensing ion channel activity in combination with other stimuli, and that additional ASIC subtypes (along with ASIC3) may serve to sense and mediate signals from multiple stimuli.     

Molecular dynamics simulations investigate the mechanism of Psalmotoxin 1 regulating gating process of an acid-sensing ion channel 1a at pH 5.5

Acid-sensing ion channel 1a (ASIC1a) is a cation channel activated by protons and causes neuronal death through central nervous system. Psalmotoxin1 (PcTx1) is a gating modifier for ASIC1a. The process of PcTx1 regulating the channel gating from the extracellular domain to the transmembrane domain is unclear. Here we used molecular dynamics (MD) simulations method to investigate how PcTx1 regulates the gating of the ASIC1a. Our results indicated that PcTx1can mainly regulate ASIC1a gating process through hydrogen bonds, which can affect their relative positions of several key domains in ASIC1a, further, a long-range conformational changes path was determined, which is composed of β1, β2, β10, α6, α7, β11, and β12 in ASIC1a.     

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.