Structural and Thermodynamic Properties of Selective Ion Binding in a K+ Channel

Thermodynamic measurements of ion binding to the Streptomyces lividans K⁺ channel were carried out using isothermal titration calorimetry, whereas atomic structures of ion-bound and ion-free conformations of the channel were characterized by x-ray crystallography. Here we use these assays to show that the ion radius dependence of selectivity stems from the channel's recognition of ion size (i.e., volume) rather than charge density. Ion size recognition is a function of the channel's ability to adopt a very specific conductive structure with larger ions (K⁺, Rb⁺, Cs⁺, and Ba²⁺) bound and not with smaller ions (Na⁺, Mg²⁺, and Ca²⁺). The formation of the conductive structure involves selectivity filter atoms that are in direct contact with bound ions as well as protein atoms surrounding the selectivity filter up to a distance of 15 Å from the ions. We conclude that ion selectivity in a K⁺ channel is a property of size-matched ion binding sites created by the protein structure.     

Determinants of cation transport selectivity: Equilibrium binding and transport kinetics

ConclusionThe equilibrium ion-binding properties of ion channels and transporters can be difficult to discern from crystal structures alone, as proteins often adopt different lowest energy states depending on the ions bound. In cases where transport is slow, their inherent ion-binding preferences can be used to infer their transport preferences. However, in cases where transport is fast, the transport selectivity can hide their equilibrium preferences by accentuating the kinetics of ions hopping through a channel over its inherent ion-binding preferences. Thus, depending on the arrangement of ion-binding sites in a channel’s selectivity filter, one can achieve either selective or nonselective ion transport.The equilibrium K⁺ selectivity of some nonselective channels suggests a potential mechanism whereby they could evolve into a fast K⁺-selective channel. K⁺ channels and nonselective channels like CNG and HCN are related to one another in both sequence and structure, suggesting an evolutionary link between them. Swap experiments show that only a few mutations separate a nonselective channel from a K⁺-selective channel. One might imagine an evolutionary path between these channels in which the equilibrium preference for a K⁺ ion in a nonselective channel evolves into a K⁺-selective channel through these few mutations to create the selective ion queue. Alternatively, a slow single-ion channel with an equilibrium and transport preference for K⁺ ions could be transformed into a fast multi-ion channel through mutations that create a queue of K⁺-selective ion-binding sites, as is seen in most K⁺ channels studied to date.In the case of multi-ion selectivity filters, such as those found in K⁺ channels, the selectivity filter can be viewed as the active site that interacts with different queues of ions and water molecules. At least three properties emerge from multi-ion queues: (1) high conductance by reducing the affinity of multiple bound ions versus single ions; (2) high selectivity by allowing disfavored ions time to dissociate back into solution; and, consequently, (3) robust selectivity in an environment where ion concentrations can change. For transporters and carriers, the equilibrium preference and slow transport naturally create robust selectivity. In all these cases, equilibrium-based ion selectivity is achieved by slowing transport enough so that the disfavored ion is able to dissociate back into solution before transport takes place.     

The selectivity of K+ ion channels: testing the hypotheses

How K⁺ channels are able to conduct certain cations yet not others remains an important but unresolved question. The recent elucidation of the structure of NaK, an ion channel that conducts both Na⁺ and K⁺ ions, offers an opportunity to test the various hypotheses that have been put forward to explain the selectivity of K⁺ ion channels. We test the snug-fit, field-strength, and over-coordination hypotheses by comparing their predictions to the results of classical molecular dynamics simulations of the K⁺ selective channel KcsA and the less selective channel NaK embedded in lipid bilayers. Our results are incompatible with the so-called strong variant of the snug-fit hypothesis but are consistent with the over-coordination hypothesis and neither confirm nor refute the field-strength hypothesis. We also find that the ions and waters in the NaK selectivity filter unexpectedly move to a new conformation in seven K⁺ simulations: the two K⁺ ions rapidly move from site S4 to S2 and from the cavity to S4. At the same time, the selectivity filter narrows around sites S1 and S2 and the carbonyl oxygen atoms rotate 20°−40° inwards toward the ion. These motions diminish the large structural differences between the crystallographic structures of the selectivity filters of NaK and KcsA and appear to allow the binding of ions to S2 of NaK at physiological temperature.     

Molecular basis of ion permeability in a voltage‐gated sodium channel

Voltage‐gated sodium channels are essential for electrical signalling across cell membranes. They exhibit strong selectivities for sodium ions over other cations, enabling the finely tuned cascade of events associated with action potentials. This paper describes the ion permeability characteristics and the crystal structure of a prokaryotic sodium channel, showing for the first time the detailed locations of sodium ions in the selectivity filter of a sodium channel. Electrostatic calculations based on the structure are consistent with the relative cation permeability ratios (Na⁺ ≈ Li⁺ ≫ K⁺, Ca²⁺, Mg²⁺) measured for these channels. In an E178D selectivity filter mutant constructed to have altered ion selectivities, the sodium ion binding site nearest the extracellular side is missing. Unlike potassium ions in potassium channels, the sodium ions in these channels appear to be hydrated and are associated with side chains of the selectivity filter residues, rather than polypeptide backbones.     

Dynamics, Energetics, and Selectivity of the Low-K KcsA Channel Structure

Potassium channels are a diverse family of integral membrane proteins through which K⁺ can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K⁺ channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K⁺ concentration is thought to be representative of a K⁺ channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K⁺ KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K⁺ KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K⁺ structure of KcsA and KirBac. The low-K⁺ KcsA structure is not permeable by Na⁺, K⁺, or Rb⁺, and the selectivity of its binding sites is different from that of the high-K⁺ structure.     

Energy distribution and ion selectivity of the bacterial potassium channel

The ion selectivity of the bacterial potassium channel K_CSA is explained upon comparing the energy characteristics of the interaction of cations (Li⁺, Na⁺, K⁺) with atoms of the selectivity filter of the protein pore. Quantum-chemical calculations reveal a deeper potential well for potassium ions, which accounts for preferred K⁺ permeation. It is shown that the conventional methods with AMBER, CHARMM, OPLS force fields in standard parametrization as well as partial re-parametrization give incorrect estimates of ion energy distribution in the channel.     

Ion binding properties and structure stability of the NaK channel

Ion distribution in the selectivity filter and ion-water and ion-protein interactions of NaK channel are systematically investigated by all-atom molecular dynamics simulations, with the tetramer channel protein being embedded in a solvated phospholipid bilayer. Analysis of the simulation results indicates that K⁺ ions prefer to bind within the sites formed by two adjacent planes of oxygen atoms from the selectivity filter, while Na⁺ ions are inclined to bind to a single plane of four oxygen atoms. At the same time, both K⁺ and Na⁺ ions can diffuse in the vestibule, accompanying with movements of the water molecules confined in a complex formed by the vestibule together with four small grottos connecting to it. As a result, K⁺ ions show a wide range of coordination numbers (6-8), while Na⁺ ions display a constant coordination number of ∼ 6 in the selectivity filter, which may result in the loss of selectivity of NaK. It is also found that a Ca²⁺ can bind at the extracellular site as reported in the crystal structure in a partially hydrated state, or at a higher site in a full hydration state. Furthermore, the carbonyl group of Asp66 can reorient to point towards the center pore when an ion exists in the vestibule, while that of Gly65 always aligns tangentially to the channel axis, as in the crystallographic structures.     

Nonselective Conduction in a Mutated NaK Channel with Three Cation-Binding Sites

The NaK channel is a cation-selective protein with similar permeability for K⁺ and Na⁺ ions. Crystallographic structures are available for the wild-type and mutated NaK channels with different numbers of cation-binding sites. We have performed a comparison between the potentials of mean force governing the translocation of K⁺ ions and mixtures of one Na⁺ and three K⁺ ions in a mutated NaK channel with only three cation-binding sites (NaK-CNG). Since NaK-CNG is not selective for K⁺ over Na⁺, analysis of its multi-ion potential energy surfaces can provide clues about how selectivity originates. Comparison of the potentials of mean force of NaK-CNG and K⁺-selective channels yields observations that strongly suggest that the number of contiguous ion binding sites in a single-file mechanism is the key determinant of the channel’s selectivity properties, as already proposed by experimental studies. We conclude that the presence of four binding sites in K⁺-selective channels is essential for highly selective and efficient permeation of K⁺ ions, and that a key difference between K⁺-selective and nonselective channels is the absence/presence of a binding site for Na⁺ ions at the boundary between S2 and S3 in the context of multi-ion permeation events.     

Preferential binding of K+ ions in the selectivity filter at equilibrium explains high selectivity of K+ channels

K⁺ channels exhibit strong selectivity for K⁺ ions over Na⁺ ions based on electrophysiology experiments that measure ions competing for passage through the channel. During this conduction process, multiple ions interact within the region of the channel called the selectivity filter. Ion selectivity may arise from an equilibrium preference for K⁺ ions within the selectivity filter or from a kinetic mechanism whereby Na⁺ ions are precluded from entering the selectivity filter. Here, we measure the equilibrium affinity and selectivity of K⁺ and Na⁺ ions binding to two different K⁺ channels, KcsA and MthK, using isothermal titration calorimetry. Both channels exhibit a large preference for K⁺ over Na⁺ ions at equilibrium, in line with electrophysiology recordings of reversal potentials and Ba²⁺ block experiments used to measure the selectivity of the external-most ion-binding sites. These results suggest that the high selectivity observed during ion conduction can originate from a strong equilibrium preference for K⁺ ions in the selectivity filter, and that K⁺ selectivity is an intrinsic property of the filter. We hypothesize that the equilibrium preference for K⁺ ions originates in part through the optimal spacing between sites to accommodate multiple K⁺ ions within the selectivity filter.     

Non-Equilibrium Dynamics Contribute to Ion Selectivity in the KcsA Channel

The ability of biological ion channels to conduct selected ions across cell membranes is critical for the survival of both animal and bacterial cells. Numerous investigations of ion selectivity have been conducted over more than 50 years, yet the mechanisms whereby the channels select certain ions and reject others are not well understood. Here we report a new application of Jarzynski’s Equality to investigate the mechanism of ion selectivity using non-equilibrium molecular dynamics simulations of Na⁺ and K⁺ ions moving through the KcsA channel. The simulations show that the selectivity filter of KcsA adapts and responds to the presence of the ions with structural rearrangements that are different for Na⁺ and K⁺. These structural rearrangements facilitate entry of K⁺ ions into the selectivity filter and permeation through the channel, and rejection of Na⁺ ions. A mechanistic model of ion selectivity by this channel based on the results of the simulations relates the structural rearrangement of the selectivity filter to the differential dehydration of ions and multiple-ion occupancy and describes a mechanism to efficiently select and conduct K⁺. Estimates of the K⁺/Na⁺ selectivity ratio and steady state ion conductance for KcsA from the simulations are in good quantitative agreement with experimental measurements. This model also accurately describes experimental observations of channel block by cytoplasmic Na⁺ ions, the “punch through” relief of channel block by cytoplasmic positive voltages, and is consistent with the knock-on mechanism of ion permeation.     

KcsA 通道对Na+、K+及Rb+离子选择性的统计热力学研究

钾离子的通透率至少比钠离子的通透率大10000倍,这个问题至今没有很好地解决.为了在分子水平阐释钾离子通道的选择性机制,以KcsA钾通道X射线衍射结构为基础,采用密度泛函理论计算了不同离子在离子通道中的位能.计算结果表明,Rb+离子具有与K+离子相类似的位能曲线,但是其在通透过程遇到的位垒要比K+离子的位垒高,因而所对应的通透率也就小于钾离子的通透率,而钠离子的的通透率仅仅是钾离子通透率的0.0067%.文中所涉及的系统仅仅包含269个原子,而用分子动力学虽然也可以得到相近的结果,但是它的系统大小为41 000个原子.     

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

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

Ion–protein interactions of a potassium ion channel studied by attenuated total reflection Fourier transform infrared spectroscopy

An understanding of ion–protein interactions is key to a better understanding of the molecular mechanisms of proteins, such as enzymes, ion channels, and ion pumps. A potassium ion channel, KcsA, has been extensively studied in terms of ion selectivity. Alkali metal cations in the selectivity filter were visualized by X-ray crystallography. Infrared spectroscopy has an intrinsically higher structural sensitivity due to frequency changes in molecular vibrations interacting with different ions. In this review article, I attempt to summarize ion-exchange-induced differences in Fourier transform infrared spectroscopy, as applied to KcsA, to explain how this method can be utilized to study ion–protein interactions in the KcsA selectivity filter. A band at 1680 cm^?1 in the amide I region would be a marker band for the ion occupancy of K⁺, Rb⁺, and Cs⁺ in the filter. The band at 1627 cm^?1 observed in both Na⁺ and Li⁺ conditions suggests that the selectivity filter similarly interacts with these ions. In addition to the structural information, the results show that the titration of K⁺ ions provides quantitative information on the ion affinity of the selectivity filter.     

Calcium-potassium selectivity: kinetic analysis of current-voltage relationships of the open, slowly activating channel in the vacuolar membrane of Vicia faba guard-cells

Single-channel current-voltage (IV ) relationships of the open, slowly activating vacuolar (SV) channel of Vicia faba L. were recorded in solutions with different activities of Ca²⁺ and K⁺, and have been analyzed for Ca²⁺/K⁺ selectivity. Two models with one binding site have been examined. A rigid-pore model with a main binding site between two energy barriers (nine free parameters) provides fair fits. Slightly better fits are obtained with an alternative, dynamic-pore model, where the selectivity filter is located between two Mitchellian ion wells of the cytoplasmic and luminal pore sections, and where the selectivity filter alternates the orientation of the binding site between the two faces of the pore (ten free parameters). Using sets of IV-relationships with only Ca²⁺ or only K⁺ as transportable substrates, both models consistently predict open-channel IV-relationships in the presence of both substrates. Fits of both models to the entire ensemble of␣data yield very similar flux-voltage characteristics for␣Ca²⁺ and for K⁺ in experimental conditions, and consistently predict such flux-voltage characteristics over physiologically relevant ranges of voltage and substrate concentrations. In a very general sense, physiological Ca⁺ fluxes through the open SV channel are predominantly inward and about 50 times smaller than K⁺ fluxes. The ions Cl⁻, OH⁻, and H⁺, do not pass the SV channel at significant rates. Kinetic details of the SV channel with respect to binding and passage of Ca²⁺ and K⁺ are discussed on the basis of the consistent results of the reaction-kinetic analysis of the experimental data by the two models. Received: 14 July 1997 / Accepted: 26 September 1997     

Selectivity sequences in a model calcium channel: role of electrostatic field strength

The energetics that give rise to selectivity sequences of ionic binding selectivity of Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺ in a model of a calcium channel are considered. This work generalizes Eisenman’s classic treatment (Biophys J 2(Suppl. 2):259, 1962) by including multiple, mobile binding site oxygens that coordinate many permeating ions (all modeled as charged, hard spheres). The selectivity filter of the model calcium channel allows the carboxyl terminal groups of glutamate and aspartate side chains to directly interact with and coordinate the permeating ions. Ion dehydration effects are represented with a Born energy between the dielectric coefficients of the selectivity filter and the bath. High oxygen concentration creates a high field strength site that prefers small ions, as in Eisenman’s model. On the other hand, a low filter dielectric constant also creates a high field strength site, but this site prefers large ions, contrary to Eisenman’s model. These results indicate that field strength does not have a unique effect on ionic binding selectivity sequences once entropic, electrostatic, and dehydration forces are included in the model. Thus, Eisenman’s classic relationship between field strength and selectivity sequences must be supplemented with additional information about selectivity filters such as the calcium channel that has amino acid side chains mixing with ions to make a crowded permeation pathway.     

Sodium Ions Do Not Stabilize the Selectivity Filter of a Potassium Channel

Ion conduction is an essential function for electrical activity in all organisms. The non-selective ion channel NaK was previously shown to adopt two stable conformations of the selectivity filter. Here, we present solid-state NMR measurements of NaK demonstrating a population shift between these conformations induced by changing the ions in the sample while the overall structure of NaK is not affected. We show that two K⁺-selective mutants (NaK2K and NaK2K-Y66F) suffer a complete loss of selectivity filter stability under Na⁺ conditions, but do not collapse into a defined structure. Widespread chemical shift perturbations are seen between the Na⁺ and K⁺ states of the K⁺-selective mutants in the region of the pore helix indicating structural changes. We conclude that the stronger link between the selectivity filter and the pore helix in the K⁺-selective mutants, compared to the non-selective wild-type NaK channel, reduces the ion-dependent conformational flexibility of the selectivity filter.     

Widening the scope of constriction

JGP modeling study suggests that selectivity filter constriction is a plausible mechanism for C-type inactivation of the Shaker voltage-gated potassium channel.

In response to prolonged activation, many K⁺ channels spontaneously reduce the membrane conductance by undergoing C-type inactivation, a kinetic process crucial for the pacing of cardiac action potentials and the modulation of neuronal firing patterns. In the pH-activated bacterial channel KcsA, C-type inactivation appears to involve constriction of the channel’s selectivity filer that prohibits ion conduction, but whether voltage-gated channels like Drosophila Shaker use a similar mechanism is controversial (1). In this issue of JGP, a computational study by Li et al. suggests that filter constriction is indeed a plausible mechanism for the C-type inactivation of Shaker (2).(Left to right) Jing Li, Benoît Roux, and colleagues use computational modeling to show that selectivity filter constriction, allosterically promoted by opening of the intracellular activation gate, is a plausible mechanism for the C-type inactivation of voltage-gated K⁺ channels such as Drosophila Shaker. The selectivity filter is conductive (left) when the intracellular gate is partially open, but adopts a constricted conformation (right) when the gate is open wide.Various structural approaches have shown that C-type inactivation of KcsA channels is associated with the symmetrical constriction of all four channel subunits at the level of the central glycine residue in the selectivity filter. Benoît Roux and colleagues at The University of Chicago used MD simulations to show that the KcsA pore can transition from the conductive to the constricted conformation on an appropriate timescale, and that this transition is allosterically promoted by the wide opening of the pore’s intracellular gate (3). Modeling by Roux and colleagues suggests that C-type inactivation of cardiac hERG channels could also involve selectivity filter constriction, though in this case it appears to be an asymmetric process in which only two of the channel’s subunits move closer together (4).“In view of the high similarity between the pore domains of Shaker and KcsA (almost 40% sequence identity), we wanted to examine if it’s possible for the Shaker selectivity filter to constrict and, if so, how similar it is to KcsA,” Roux explains. Led by first author Jing Li—now an assistant professor at the University of Mississippi—Roux and colleagues developed several homology models of the Shaker pore domain with the intracellular gate open to various degrees (2).MD simulations and free energy calculations revealed that the Shaker selectivity filter can dynamically transition from a conductive to a constricted conformation, and that this transition is allosterically coupled to the intracellular gate; the constricted conformation is stable when the gate is wide open. “Our computations strongly suggest that constriction is a plausible mechanism for the C-type inactivation of Shaker,” Roux says. “There’s no reason based on the currently available information to reject the existence of a constricted state in Shaker channels.”As with KcsA, Shaker channels appear to constrict symmetrically at the level of the selectivity filter’s central glycine. But Li et al.’s simulations revealed some small variations between the two channels, including differences in the number of water molecules bound to each channel subunit and the arrangement of the hydrogen-bond network they form to stabilize the constricted state.Li et al. also modeled the pore domain of the Shaker W434F mutant, which is widely assumed to be trapped in a C-type inactivated state. The simulation suggests that the mutant channel’s filter adopts a stable constricted conformation even when the intracellular gate is only partially open, although the constriction is asymmetric and occurs at the level of a different filter residue (2).Constriction may therefore be a universal mechanism of C-type inactivation, even if the exact conformation varies from channel to channel. But, says Roux, confirming this will require more experimental work using the right conditions and mutations to capture the structure of inactivated channels.     

The conserved potassium channel filter can have distinct ion binding profiles: Structural analysis of rubidium,cesium, and barium binding in NaK2K

Potassium channels are highly selective for K⁺ over the smaller Na⁺. Intriguingly, they are permeable to larger monovalent cations such as Rb⁺ and Cs⁺ but are specifically blocked by the similarly sized Ba²⁺. In this study, we used structural analysis to determine the binding profiles for these permeant and blocking ions in the selectivity filter of the potassium-selective NaK channel mutant NaK2K and also performed permeation experiments using single-channel recordings. Our data revealed that some ion binding properties of NaK2K are distinct from those of the canonical K⁺ channels KcsA and MthK. Rb⁺ bound at sites 1, 3, and 4 in NaK2K, as it does in KcsA. Cs⁺, however, bound predominantly at sites 1 and 3 in NaK2K, whereas it binds at sites 1, 3, and 4 in KcsA. Moreover, Ba²⁺ binding in NaK2K was distinct from that which has been observed in KcsA and MthK, even though all of these channels show similar Ba²⁺ block. In the presence of K⁺, Ba²⁺ bound to the NaK2K channel at site 3 in conjunction with a K⁺ at site 1; this led to a prolonged block of the channel (the external K⁺-dependent Ba²⁺ lock-in state). In the absence of K⁺, however, Ba²⁺ acts as a permeating blocker. We found that, under these conditions, Ba²⁺ bound at sites 1 or 0 as well as site 3, allowing it to enter the filter from the intracellular side and exit from the extracellular side. The difference in the Ba²⁺ binding profile in the presence and absence of K⁺ thus provides a structural explanation for the short and prolonged Ba²⁺ block observed in NaK2K.     

A Structural Model of the Pore-Forming Region of the Skeletal Muscle Ryanodine Receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca²⁺ over K⁺ in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca²⁺ in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K⁺ but not Ca²⁺. Our simulations with G4898R likewise show a decrease in the preference of Ca²⁺ over K⁺ in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level.     

Ionic interactions of Ba2+ blockades in the MthK K+ channel

The movement and interaction of multiple ions passing through in single file underlie various fundamental K⁺ channel properties, from the effective conduction of K⁺ ions to channel blockade by Ba²⁺ ions. In this study, we used single-channel electrophysiology and x-ray crystallography to probe the interactions of Ba²⁺ with permeant ions within the ion conduction pathway of the MthK K⁺ channel. We found that, as typical of K⁺ channels, the MthK channel was blocked by Ba²⁺ at the internal side, and the Ba²⁺-blocking effect was enhanced by external K⁺. We also obtained crystal structures of the MthK K⁺ channel pore in both Ba²⁺–Na⁺ and Ba²⁺–K⁺ environments. In the Ba²⁺–Na⁺ environment, we found that a single Ba²⁺ ion remained bound in the selectivity filter, preferably at site 2, whereas in the Ba²⁺–K⁺ environment, Ba²⁺ ions were predominantly distributed between sites 3 and 4. These ionic configurations are remarkably consistent with the functional studies and identify a molecular basis for Ba²⁺ blockade of K⁺ channels.