Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA,PIIIA, and KIIIA in a NavAb-based model of Nav1.4

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.     

Sodium channels: ionic model of slow inactivation and state-dependent drug binding

Inactivation is a fundamental property of voltage-gated ion channels. Fast inactivation of Na(+) channels involves channel block by the III-IV cytoplasmic interdomain linker. The mechanisms of nonfast types of inactivation (intermediate, slow, and ultraslow) are unclear, although the ionic environment and P-loops rearrangement appear to be involved. In this study, we employed a TTX-based P-loop domain model of a sodium channel and the MCM method to investigate a possible role of P-loop rearrangement in the nonfast inactivation. Our modeling predicts that Na(+) ions can bind between neighboring domains in the outer-carboxylates ring EEDD, forming an ordered structure with interdomain contacts that stabilize the conducting conformation of the outer pore. In this model, the permeant ions can transit between the EEDD ring and the selectivity filter ring DEKA, retaining contacts with at least two carboxylates. In the absence of Na(+), the electrostatic repulsion between the EEDD carboxylates disrupts the permeable configuration. In this Na(+)-deficient model, the region between the EEDD and DEKA rings is inaccessible for Na(+) but is accessible for TMA. Taken together, these results suggest that Na(+)-saturated models are consistent with experimental characteristics of the open channels, whereas Na(+)-deficient models are consistent with experimentally defined properties of the slow-inactivated channels. Our calculations further predict that binding of LAs to the inner pore would depend on whether Na(+) occupies the DEKA ring. In the absence of Na(+) in the DEKA ring, the cationic group of lidocaine occurs in the focus of the pore helices' macrodipoles and would prevent occupation of the ring by Na(+). Loading the DEKA ring with Na(+) results in the electrostatic repulsion with lidocaine. Thus, there are antagonistic relations between a cationic ligand bound in the inner pore and Na(+) in the DEKA ring.     

Structural modeling of calcium binding in the selectivity filter of the L-type calcium channel

Calcium channels play crucial physiological roles. In the absence of high-resolution structures of the channels, the mechanism of ion permeation is unknown. Here we used a method proposed in an accompanying paper (Cheng and Zhorov in Eur Biophys J, 2009) to predict possible chelation patterns of calcium ions in a structural model of the L-type calcium channel. We compared three models in which two or three calcium ions interact with the four selectivity filter glutamates and a conserved aspartate adjacent to the glutamate in repeat II. Monte Carlo energy minimizations yielded many complexes with calcium ions bound to at least two selectivity filter carboxylates. In these complexes calcium-carboxylate attractions are counterbalanced by calcium-calcium and carboxylate-carboxylate repulsions. Superposition of the complexes suggests a high degree of mobility of calcium ions and carboxylate groups of the glutamates. We used the predicted complexes to propose a permeation mechanism that involves single-file movement of calcium ions. The key feature of this mechanism is the presence of bridging glutamates that coordinate two calcium ions and enable their transitions between different chelating patterns involving four to six oxygen atoms from the channel protein. The conserved aspartate is proposed to coordinate a calcium ion incoming to the selectivity filter from the extracellular side. Glutamates in repeats III and IV, which are most distant from the repeat II aspartate, are proposed to coordinate the calcium ion that leaves the selectivity filter to the inner pore. Published experimental data and earlier proposed permeation models are discussed in view of our model.     

Modeling P-loops domain of sodium channel: homology with potassium channels and interaction with ligands

A large body of experimental data on Na+ channels is available, but the interpretation of these data in structural terms is difficult in the absence of a high-resolution structure. Essentially different electrophysiological and pharmacological properties of Na+ and K+ channels and poor identity of their sequences obstruct homology modeling of Na+ channels. In this work, we built the P-loops model of the Na+ channel, in which the pore helices are arranged exactly as in the MthK bacterial K+ channel. The conformation of the selectivity-filter region, which includes residues in positions -2 through +4 from the DEKA locus, was shaped around rigid molecules of saxitoxin and tetrodotoxin that are known to form multiple contacts with this region. Intensive Monte Carlo minimization that started from the MthK-like conformation produced practically identical saxitoxin- and tetrodotoxin-based models. The latter was tested to explain a wide range of experimental data that were not used at the model building stage. The docking of tetrodotoxin analogs unambiguously predicted their optimal orientation and the interaction energy that correlates with the experimental activity. The docking of mu-conotoxin produced a binding model consistent with experimentally known toxin-channel contacts. Monte Carlo-minimized energy profiles of tetramethylammonium pulled through the selectivity-filter region explain the paradoxical experimental data that this organic cation permeates via the DEAA but not the AAAA mutant of the DEKA locus. The model is also consistent with earlier proposed concepts on the Na+ channel selectivity as well as Ca2+ selectivity of the EEEE mutant of the DEKA locus. Thus, the model integrates available experimental data on the Na+ channel P-loops domain, and suggests that it is more similar to K+ channels than was believed before.     

In silico activation of KcsA K+ channel by lateral forces applied to the C-termini of inner helices

Crystallographic studies of K(+) channels in the closed (KcsA) and open (MthK) states suggest that Gly(99) (KcsA numbering) in the inner helices serves as a gating hinge during channel activation. However, some P-loop channels have larger residues in the corresponding position. The comparison of x-ray structures of KcsA and MthK shows that channel activation alters backbone torsions and helical H-bonds in residues 95-105. Importantly, the changes in Gly(99) are not the largest ones. This raises questions about the mechanism of conformational changes upon channel gating. In this work, we have built a model of the open KcsA using MthK as a template and simulated opening and closing of KcsA by constraining C-ends of the inner helices at a gradually changing distance from the pore axis without restraining mobility of the helices along the axis. At each imposed distance, the energy was Monte Carlo-minimized. The channel-opening and channel-closing trajectories arrived to the structures in which the backbone geometry was close to that seen in MthK and KcsA, respectively. In the channel-opening trajectory, the constraints-induced lateral forces caused kinks at midpoints of the inner helices between Val(97) and Gly(104) but did not destroy interdomain contacts, the pore helices, and the selectivity filter. The simulated activation of the Gly(99)Ala mutant yielded essentially similar results. Analysis of interresidue energies shows that the N-terminal parts of the inner helices form strong attractive contacts with the pore helices and the outer helices. The lateral forces induce kinks at the position where the helix-breaking torque is maximal and the intersegment contacts vanish. This mechanism may be conserved in different P-loop channels.     

Use-dependent block of the voltage-gated Na+ channel by tetrodotoxin and saxitoxin: Effect of pore mutations that change ionic selectivity

Voltage-gated Na⁺ channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca²⁺ permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca²⁺ or Na⁺ ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca²⁺ permeability, suggesting that ion–toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.     

An electrostatic interaction between TEA and an introduced pore aromatic drives spring-in-the-door inactivation in Shaker potassium channels

Slow inactivation of Kv1 channels involves conformational changes near the selectivity filter. We examine such changes in Shaker channels lacking fast inactivation by considering the consequences of mutating two residues, T449 just external to the selectivity filter and V438 in the pore helix near the bottom of the selectivity filter. Single mutant T449F channels with the native V438 inactivate very slowly, and the canonical foot-in-the-door effect of extracellular tetraethylammonium (TEA) is not only absent, but the time course of slow inactivation is accelerated by TEA. The V438A mutation dramatically speeds inactivation in T449F channels, and TEA slows inactivation exactly as predicted by the foot-in-the-door model. We propose that TEA has this effect on V438A/T449F channels because the V438A mutation produces allosteric consequences within the selectivity filter and may reorient the aromatic ring at position 449. We investigated the possibility that the blocker promotes the collapse of the outer vestibule (spring-in-the-door) in single mutant T449F channels by an electrostatic attraction between a cationic TEA and the quadrupole moments of the four aromatic rings. To test this idea, we used in vivo nonsense suppression to serially fluorinate the introduced aromatic ring at the 449 position, a manipulation that withdraws electrons from the aromatic face with little effect on the shape, net charge, or hydrophobicity of the aromatic ring. Progressive fluorination causes monotonically enhanced rates of inactivation. In further agreement with our working hypothesis, increasing fluorination of the aromatic gradually transforms the TEA effect from spring-in-the-door to foot-in-the-door. We further substantiate our electrostatic hypothesis by quantum mechanical calculations.     

A molecular switch between the outer and the inner vestibules of the voltage-gated Na+ channel

Voltage-gated ion channels are transmembrane proteins that undergo complex conformational changes during their gating transitions. Both functional and structural data from K(+) channels suggest that extracellular and intracellular parts of the pore communicate with each other via a trajectory of interacting amino acids. No crystal structures are available for voltage-gated Na(+) channels, but functional data suggest a similar intramolecular communication involving the inner and outer vestibules. However, the mechanism of such communication is unknown. Here, we report that amino acid Ile-1575 in the middle of transmembrane segment 6 of domain IV (DIV-S6) in the adult rat skeletal muscle isoform of the voltage-gated sodium channel (rNa(V)1.4) may act as molecular switch allowing for interaction between outer and inner vestibules. Cysteine scanning mutagenesis of the internal part of DIV-S6 revealed that only mutations at site 1575 rescued the channel from a unique kinetic state ("ultra-slow inactivation," I(US)) produced by the mutation K1237E in the selectivity filter. A similar effect was seen with I1575A. Previously, we reported that conformational changes of both the internal and the external vestibule are involved in the generation of I(US). The fact that mutations at site 1575 modulate I(US) produced by K1237E strongly suggests an interaction between these sites. Our data confirm a previously published molecular model in which Ile-1575 of DIV-S6 is in close proximity to Lys-1237 of the selectivity filter. Furthermore, these functional data define the position of the selectivity filter relative to the adjacent DIV-S6 segment within the ionic permeation pathway.     

Mechanism of tetrodotoxin block and resistance in sodium channels

Tetrodotoxin (TTX) has been used for many decades to characterize the structure and function of biological ion channels. Yet, the precise mechanism by which TTX blocks voltage-gated sodium (Na_V) channels is not fully understood. Here molecular dynamics simulations are used to elucidate how TTX blocks mammalian voltage-gated sodium (Nav) channels and why it fails to be effective for the bacterial sodium channel, Na_VAb. We find that, in Na_VAb, a sodium ion competes with TTX for the binding site at the extracellular end of the filter, thus reducing the blocking efficacy of TTX. Using a model of the skeletal muscle channel, Na_V1.4, we show that the conduction properties of the channel observed experimentally are faithfully reproduced. We find that TTX occludes the entrance of Na_V1.4 by forming a network of hydrogen-bonds at the outer lumen of the selectivity filter. The guanidine group of TTX adopts a lateral orientation, rather than pointing into the filter as proposed previously. The acidic residues just above the selectivity filter are important in stabilizing the hydrogen-bond network between TTX and Na_V1.4. The effect of two single mutations of a critical tyrosine residue in the filter of Na_V1.4 on TTX binding observed experimentally is reproduced using computational mutagenesis.     

Overlapping Binding Sites of Structurally Different Antiarrhythmics Flecainide and Propafenone in the Subunit Interface of Potassium Channel Kv2.1

Kv2.1 channels, which are expressed in brain, heart, pancreas, and other organs and tissues, are important targets for drug design. Flecainide and propafenone are known to block Kv2.1 channels more potently than other Kv channels. Here, we sought to explore structural determinants of this selectivity. We demonstrated that flecainide reduced the K⁺ currents through Kv2.1 channels expressed in Xenopus laevis oocytes in a voltage- and time-dependent manner. By systematically exchanging various segments of Kv2.1 with those from Kv1.2, we determined flecainide-sensing residues in the P-helix and inner helix S6. These residues are not exposed to the inner pore, a conventional binding region of open channel blockers. The flecainide-sensing residues also contribute to propafenone binding, suggesting overlapping receptors for the drugs. Indeed, propafenone and flecainide compete for binding in Kv2.1. We further used Monte Carlo-energy minimizations to map the receptors of the drugs. Flecainide docking in the Kv1.2-based homology model of Kv2.1 predicts the ligand ammonium group in the central cavity and the benzamide moiety in a niche between S6 and the P-helix. Propafenone also binds in the niche. Its carbonyl group accepts an H-bond from the P-helix, the amino group donates an H-bond to the P-loop turn, whereas the propyl group protrudes in the pore and blocks the access to the selectivity filter. Thus, besides the binding region in the central cavity, certain K⁺ channel ligands can expand in the subunit interface whose residues are less conserved between K⁺ channels and hence may be targets for design of highly desirable subtype-specific K⁺ channel drugs.     

Gating pore currents in DIIS4 mutations of NaV1.4 associated with periodic paralysis: saturation of ion flux and implications for disease pathogenesis

S4 voltage–sensor mutations in CaV1.1 and NaV1.4 channels cause the human muscle disorder hypokalemic periodic paralysis (HypoPP). The mechanism whereby these mutations predispose affected sarcolemma to attacks of sustained depolarization and loss of excitability is poorly understood. Recently, three HypoPP mutations in the domain II S4 segment of NaV1.4 were shown to create accessory ionic permeation pathways, presumably extending through the aqueous gating pore in which the S4 segment resides. However, there are several disparities between reported gating pore currents from different investigators, including differences in ionic selectivity and estimates of current amplitude, which in turn have important implications for the pathological relevance of these aberrant currents. To clarify the features of gating pore currents arising from different DIIS4 mutants, we recorded gating pore currents created by HypoPP missense mutations at position R666 in the rat isoform of Nav1.4 (the second arginine from the outside, at R672 in human NaV1.4). Extensive measurements were made for the index mutation, R666G, which created a gating pore that was permeable to K⁺ and Na⁺. This current had a markedly shallow slope conductance at hyperpolarized voltages and robust inward rectification, even when the ionic gradient strongly favored outward ionic flow. These characteristics were accounted for by a barrier model incorporating a voltage-gated permeation pathway with a single cation binding site oriented near the external surface of the electrical field. The amplitude of the R666G gating pore current was similar to the amplitude of a previously described proton-selective current flowing through the gating pore in rNaV1.4-R663H mutant channels. Currents with similar amplitude and cation selectivity were also observed in R666S and R666C mutant channels, while a proton-selective current was observed in R666H mutant channels. These results add support to the notion that HypoPP mutations share a common biophysical profile comprised of a low-amplitude inward current at the resting potential that may contribute to the pathological depolarization during attacks of weakness.     

Altered KCNQ3 Potassium Channel Function Caused by the W309R Pore-Helix Mutation Found in Human Epilepsy

The second tryptophan (W) residue of the conserved WW motif in the pore helix of many K⁺ channel subunit is thought to interact with the tyrosine (Y) residues of the selectivity filter. A missense mutation causing the replacement of the corresponding residues with an arginine (W309R) occurs in KCNQ3 subunits forming part of M-channels. In this study, we examined the functional consequences of the W309R mutation in heterogously expressed KCNQ channels. Homomeric KCNQ3^W309R channels lacked KCNQ currents. Heteromeric KCNQ2/KCNQ3^W309R channels displayed a dominant-negative suppression of current and a significant modification in gating properties when compared with heteromeric KCNQ3/KCNQ2 channels mimicking the M-channels. A three-dimensional homology model in the W309R mutant indicated that the R side chain of pore helices is too far from the Y side chain of the selectivity filter to interact via hydrogen bonds with each other and stabilize the pore structure. Collectively, the present results suggest that the second W residues of pore helices and their chemical interaction with the Y residues of the selectivity filter are essential for normal K⁺ channel function. This pore-helix mutation, if occurs in the brain M channels, could thus lead to a channel dysfunction sufficient to trigger epileptic hyperexcitability.     

The tryptophan cluster: a hypothetical structure of the DNA-binding domain of the myb protooncogene product

In the DNA-binding domain of the c-myb protooncogene product (c-Myb) which consists of three repeats of 51-52 amino acids, there are 3 perfectly conserved tryptophans in each repeat. Site-directed mutagenesis of these tryptophans showed that any single or multiple mutations of tryptophan to hydrophilic residues or alanine abolished or greatly reduced the sequence-specific DNA-binding activity, but mutations to hydrophobic amino acids retained considerable activity. Raman spectroscopic study showed that these tryptophans were buried in the protein core. These 3 tryptophans are proposed to form a cluster in the hydrophobic core in each repeat. This hypothetical structure is referred to as the "tryptophan cluster," and it may represent a characteristic property of a group of DNA-binding proteins including the myb- and ets-related proteins.     

Identification of a novel pharmacophore for peptide toxins interacting with K+ channels

KappaM-conotoxin RIIIK blocks TSha1 K+ channels from trout with high affinity by interacting with the ion channel pore. As opposed to many other peptides targeting K+ channels, kappaM-RIIIK does not possess a functional dyad. In this study we combine thermodynamic mutant cycle analysis and docking calculations to derive the binding mode of kappaM-conotoxin RIIIK to the TSha1 channel. The final model reveals a novel pharmacophore, where no positively charged side chain occludes the channel pore. Instead the positive-charged residues of the toxin form a basic ring; kappaM-RIIIK is anchored to the K+ channel via electrostatic interactions of this basic ring with the loop and pore helix residues of the channel. The channel amino acid Glu-354 is likely to be a fundamental determinant of the selectivity of kappaM-RIIIK for the TSha1 channel. The Cgamma-OH of Hyp-15 is in contact with the carbonyls of the selectivity filter, disturbing the charge distribution pattern necessary for the coordination of K+ ions. This novel, experimentally based pharmacophore model proves the existence of diverse binding modes of peptidic toxins to K+ channels and underlines the role of intermolecular electrostatic interactions involving channel loop side chains in determining the selectivity of toxins for specific K+ channel types.     

Forced Dissociation of the Strand Dimer Interface between C-Cadherin Ectodomains

The force-induced dissociation of the strand dimer interface in C-cadherin has been studied using steered molecular dynamics simulations. The dissociation occurred, without domain unraveling, after the extraction of the conserved trypthophans (Trp2) from their respective hydrophobic pockets. The simulations revealed two stable positions for the Trp2 side chain inside the pocket. The most internal stable position involved a hydrogen bond between the ring Ne of Trp2 and the backbone carbonyl of Glu90. In the second stable position, the aromatic ring is located at the pocket entrance. After extracting the two tryptophans from their pockets, the complex exists in an intermediate bound state that involves a close packing of the tryptophans with residues Asp1 and Asp27 from both domains. Dissociation occurred after this residue association was broken. Simulations carried out with a complex formed between W2A mutants showed that the mutant complex dissociates more easily than the wild type complex does. These results correlate closely with the role of the conserved tryptophans suggested previously by site directed mutagenesis.     

Electrostatics studies and molecular dynamics simulations of a homology model of the Shaker K+ channel pore

A homology model of the pore domain of the Shaker K+ channel has been constructed using a bacterial K+ channel, KcsA, as a template structure. The model is in agreement with mutagenesis and sequence variability data. A number of structural features are conserved between the two channels, including a ring of tryptophan sidechains on the outer surface of the pore domain at the extracellular end of the helix bundle, and rings of acidic sidechains close to the extracellular mouth of the channel. One of these rings, that formed by four Asp447 sidechains at the mouth of the Shaker pore, is shown by pK(A) calculations to be incompletely ionized at neutral pH. The potential energy profile for a K+ ion moved along the central axis of the Shaker pore domain model selectivity filter reveals a shallow well, the depth of which is modulated by the ionization state of the Asp447 ring. This is more consistent with the high cation flux exhibited by the channel in its conductance value of 19 pS.     

T-type Channels Become Highly Permeable to Sodium Ions Using an Alternative Extracellular Turret Region (S5-P) Outside the Selectivity Filter

T-type (Ca_v3) channels are categorized as calcium channels, but invertebrate ones can be highly sodium-selective channels. We illustrate that the snail LCa_v3 T-type channel becomes highly sodium-permeable through exon splicing of an extracellular turret and descending helix in domain II of the four-domain Ca_v3 channel. Highly sodium-permeable T-type channels are generated without altering the invariant ring of charged residues in the selectivity filter that governs calcium selectivity in calcium channels. The highly sodium-permeant T-type channel expresses in the brain and is the only splice isoform expressed in the snail heart. This unique splicing of turret residues offers T-type channels a capacity to serve as a pacemaking sodium current in the primitive heart and brain in lieu of Na_v1-type sodium channels and to substitute for voltage-gated sodium channels lacking in many invertebrates. T-type channels would also contribute substantially to sodium leak conductances at rest in invertebrates because of their large window currents.     

Mutations in the K+ channel signature sequence.

Potassium channels share a highly conserved stretch of eight amino acids, a K+ channel signature sequence. The conserved sequence falls within the previously defined P-region of voltage-activated K+ channels. In this study we investigate the effect of mutations in the signature sequence of the Shaker channel on K+ selectivity determined under bi-ionic conditions. Nonconservative substitutions of two threonine residues and the tyrosine residue leave selectivity intact. In contrast, mutations at some positions render the channel nonselective among monovalent cations. These findings are consistent with a proposal that the signature sequence contributes to a selectivity filter. Furthermore, the results illustrate that the hydroxyl groups at the third and fourth positions, and the aromatic group at position seven, are not essential in determining K+ selectivity.     

Intersegment hydrogen bonds as possible structural determinants of the N/Q/R site in glutamate receptors.

Specific electrophysiological and pharmacological properties of ionic channels in NMDA, AMPA, and kainate subtypes of ionotropic glutamate receptors (GluRs) are determined by the Asn (N), Gln (Q), and Arg (R) residues located at homologous positions of the pore-lining M2 segments (the N/Q/R site). Presumably, the N/Q/R site is located at the apex of the reentrant membrane loop and forms the narrowest constriction of the pore. Although the shorter Asn residues are expected to protrude in the pore to a lesser extent than the longer Gln residues, the effective dimension of the NMDA channel (corresponding to the size of the largest permeant organic cation) is, surprisingly, smaller than that of the AMPA channel. To explain this paradox, we propose that the N/Q/R residues form macrocyclic structures (rings) stabilized by H-bonds between a NH(2) group in the side chain of a given M2 segment and a C==O group of the main chain in the adjacent M2 segment. Using Monte Carlo minimization, we have explored conformational properties of the rings. In the Asn, but not in the Gln ring, the side-chain oxygens protruding into the pore may facilitate ion permeation and accept H-bonds from the blocking drugs. In this way, the model explains different electrophysiological and pharmacological properties of NMDA and non-NMDA GluR channels. The ring of H-bonded polar residues at the pore narrowing resembles the ring of four Thr(75) residues observed in the crystallographic structure of the KcsA K(+) channel.