Role of protein dynamics in ion selectivity and allosteric coupling in the NaK channel

Non-ionotropic voltage-gated calcium channel signaling

Voltage-gated calcium channels (VGCCs) are the major conduits for calcium ions (Ca2+) within excitable cells. Recent studies have highlighted the non-ionotropic functionality of VGCCs, revealing their capacity to activate intracellular pathways independently of ion flow. This non-ionotropic signaling mode plays a pivotal role in excitation-coupling processes, including gene transcription through excitation-transcription (ET), synaptic transmission via excitation-secretion (ES), and cardiac contraction through excitation-contraction (EC). However, it is noteworthy that these excitation-coupling processes require extracellular calcium (Ca2+) and Ca2+ occupancy of the channel ion pore. Analogous to the "non-canonical" characterization of the non-ionotropic signaling exhibited by the N-methyl-D-aspartate receptor (NMDA), which requires extracellular Ca2+ without the influx of ions, VGCC activation requires depolarization-triggered conformational change(s) concomitant with Ca2+ binding to the open channel. Here, we discuss the contributions of VGCCs to ES, ET, and EC coupling as Ca2+ binding macromolecules that transduces external stimuli to intracellular input prior to elevating intracellular Ca2+. We emphasize the recognition of calcium ion occupancy within the open ion-pore and its contribution to the excitation coupling processes that precede the influx of calcium. The non-ionotropic activation of VGCCs, triggered by the upstroke of an action potential, provides a conceptual framework to elucidate the mechanistic aspects underlying the microseconds nature of synaptic transmission, cardiac contractility, and the rapid induction of first-wave genes.

Full opening of helix bundle crossing does not lead to NaK channel activation

A critical part of ion channel function is the ability to open and close in response to stimuli and thus conduct ions in a regulated fashion. While x-ray diffraction studies of ion channels suggested a general steric gating mechanism located at the helix bundle crossing (HBC), recent functional studies on several channels indicate that the helix bundle crossing is wide-open even in functionally nonconductive channels. Two NaK channel variants were crystallized in very different open and closed conformations, which served as important models of the HBC gating hypothesis. However, neither of these NaK variants is conductive in liposomes unless phenylalanine 92 is mutated to alanine (F92A). Here, we use NMR to probe distances at near-atomic resolution of the two NaK variants in lipid bicelles. We demonstrate that in contrast to the crystal structures, both NaK variants are in a fully open conformation, akin to Ca2+-bound MthK channel structure where the HBC is widely open. While we were not able to determine what a conductive NaK structure is like, our further inquiry into the gating mechanism suggests that the selectivity filter and pore helix are coupled to the M2 helix below and undergo changes in the structure when F92 is mutated. Overall, our data show that NaK exhibits coupling between the selectivity filter and HBC, similar to K+ channels, and has a more complex gating mechanism than previously thought, where the full opening of HBC does not lead to channel activation.

A selectivity filter mutation provides insights into gating regulation of a K+ channel

G-protein coupled inwardly rectifying potassium (GIRK) channels are key players in inhibitory neurotransmission in heart and brain. We conducted molecular dynamics simulations to investigate the effect of a selectivity filter (SF) mutation, G154S, on GIRK2 structure and function. We observe mutation-induced loss of selectivity, changes in ion occupancy and altered filter geometry. Unexpectedly, we reveal aberrant SF dynamics in the mutant to be correlated with motions in the binding site of the channel activator Gβγ. This coupling is corroborated by electrophysiological experiments, revealing that GIRK2_wt activation by Gβγ reduces the affinity of Ba²⁺ block. We further present a functional characterization of the human GIRK2_G154S mutant validating our computational findings. This study identifies an allosteric connection between the SF and a crucial activator binding site. This allosteric gating mechanism may also apply to other potassium channels that are modulated by accessory proteins.

Gly selectivity filter (TIGYGYR) mutant of the GIRK2 channel causes rare but severe neurological disorder called the Keppen-Lubinsky syndrome. Here, the authors explore the molecular mechanism of action of this glycine to serine mutant causing disease and identify an allosteric connection between the selectivity filter and a crucial activator binding site.

A biomimetic metal–organic framework with cuboid inner cavities for enantioselective separation

A biomimetic metal–organic framework with cuboid inner cavities and multiple recognition sites was constructed from a phenylalanine-derived ligand. It can enantioselectively separate various racemic alcohols, diols and epoxides with ee up to 99.5%.

Ion-dependent structure, dynamics, and allosteric coupling in a non-selective cation channel

The selectivity filter (SF) determines which ions are efficiently conducted through ion channel pores. NaK is a non-selective cation channel that conducts Na⁺ and K⁺ with equal efficiency. Crystal structures of NaK suggested a rigid SF structure, but later solid-state NMR and MD simulations questioned this interpretation. Here, we use solution NMR to characterize how bound Na⁺ vs. K⁺ affects NaK SF structure and dynamics. We find that the extracellular end of the SF is flexible on the ps-ns timescale regardless of bound ion. On a slower timescale, we observe a structural change between the Na⁺ and K⁺-bound states, accompanied by increased structural heterogeneity in Na⁺. We also show direct evidence that the SF structure is communicated to the pore via I88 on the M2 helix. These results support a dynamic SF with multiple conformations involved in non-selective conduction. Our data also demonstrate allosteric coupling between the SF and pore-lining helices in a non-selective cation channel that is analogous to the allosteric coupling previously demonstrated for K⁺-selective channels, supporting the generality of this model.

NaK is a bacterial non-selective cation channel. Here, the authors use solution NMR to show that selectivity filter (SF) in NaK is dynamic, with structural differences between the Na+ and K + -bound states. The conformation of the SF is communicated to the pore-lining helices similarly as in the K + -selective channels.

Function-Related Dynamics in Multi-Spanning Helical Membrane Proteins Revealed by Solution NMR

A primary biological function of multi-spanning membrane proteins is to transfer information and/or materials through a membrane by changing their conformations. Therefore, particular dynamics of the membrane proteins are tightly associated with their function. The semi-atomic resolution dynamics information revealed by NMR is able to discriminate function-related dynamics from random fluctuations. This review will discuss several studies in which quantitative dynamics information by solution NMR has contributed to revealing the structural basis of the function of multi-spanning membrane proteins, such as ion channels, GPCRs, and transporters.

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.

TOK channels use the two gates in classical K+ channels to achieve outward rectification

TOKs are outwardly rectifying K⁺ channels in fungi with two pore-loops and eight transmembrane spans. Here, we describe the TOKs from four pathogens that cause the majority of life-threatening fungal infections in humans. These TOKs pass large currents only in the outward direction like the canonical isolate from Saccharomyces cerevisiae (ScTOK), and distinct from other K⁺ channels. ScTOK, AfTOK1 (Aspergillus fumigatus), and H99TOK (Cryptococcus neoformans grubii) are K⁺ -selective and pass current above the K⁺ reversal potential. CaTOK (Candida albicans) and CnTOK (Cryptococcus neoformans neoformans) pass both K⁺ and Na⁺ and conduct above a reversal potential reflecting the mixed permeability of their selectivity filter. Mutations in CaTOK and ScTOK at sites homologous to those that open the internal gates in classical K⁺ channels are shown to produce inward TOK currents. A favored model for outward rectification is proposed whereby the reversal potential determines ion occupancy, and thus, conductivity, of the selectivity filter gate that is coupled to an imperfectly restrictive internal gate, permitting the filter to sample ion concentrations on both sides of the membrane.

Spotlight on the Ballet of Proteins: The Structural Dynamic Properties of Proteins Illuminated by Solution NMR

Solution NMR spectroscopy is a unique and powerful technique that has the ability to directly connect the structural dynamics of proteins in physiological conditions to their activity and function. Here, we summarize recent studies in which solution NMR contributed to the discovery of relationships between key dynamic properties of proteins and functional mechanisms in important biological systems. The capacity of NMR to quantify the dynamics of proteins over a range of time scales and to detect lowly populated protein conformations plays a critical role in its power to unveil functional protein dynamics. This analysis of dynamics is not only important for the understanding of biological function, but also in the design of specific ligands for pharmacologically important proteins. Thus, the dynamic view of structure provided by NMR is of importance in both basic and applied biology.

Conformational changes upon gating of KirBac1.1 into an open-activated state revealed by solid-state NMR and functional assays

Inward rectifier K⁺ (Kir) channels play an important role in reestablishing the resting membrane state of the action potential of excitable cells in humans. KirBac1.1 is a prokaryotic Kir channel with a high degree of homology to human Kir channels and can be isotopically labeled in NMR quantities for structural studies. Functional assays and NMR assignments confirm that KirBac1.1 is in a constitutively conductive state. Solid-state NMR assignments further reveal alternate conformations at key sites in the protein that are well conserved through human Kir channels, hinting at a possible allosteric network between channels. These underlying sequential and structural motifs could explain abnormal conductive properties of these channels fundamental to their native gating processes.

The conformational changes required for activation and K⁺ conduction in inward-rectifier K⁺ (Kir) channels are still debated. These structural changes are brought about by lipid binding. It is unclear how this process relates to fast gating or if the intracellular and extracellular regions of the protein are coupled. Here, we examine the structural details of KirBac1.1 reconstituted into both POPC and an activating lipid mixture of 3:2 POPC:POPG (wt/wt). KirBac1.1 is a prokaryotic Kir channel that shares homology with human Kir channels. We establish that KirBac1.1 is in a constitutively active state in POPC:POPG bilayers through the use of real-time fluorescence quenching assays and Förster resonance energy transfer (FRET) distance measurements. Multidimensional solid-state NMR (SSNMR) spectroscopy experiments reveal two different conformers within the transmembrane regions of the protein in this activating lipid environment, which are distinct from the conformation of the channel in POPC bilayers. The differences between these three distinct channel states highlight conformational changes associated with an open activation gate and suggest a unique allosteric pathway that ties the selectivity filter to the activation gate through interactions between both transmembrane helices, the turret, selectivity filter loop, and the pore helix. We also identify specific residues involved in this conformational exchange that are highly conserved among human Kir channels.

NMR Structural Analysis of Isolated Shaker Voltage-Sensing Domain in LPPG Micelles

The voltage-sensing domain (VSD) is a conserved structural module that regulates the gating of voltage-dependent ion channels in response to a change in membrane potential. Although the structures of many VSD-containing ion channels are now available, our understanding of the structural dynamics associated with gating transitions remains limited. To probe dynamics with site-specific resolution, we utilized NMR spectroscopy to characterize the VSD derived from Shaker potassium channel in 1-palmitoyl-2-hydroxy-sn-glycero-3-phospho-(1'-rac-glycerol) (LPPG) micelles. The backbone dihedral angles predicted based on secondary chemical shifts using torsion angle likeliness obtained from shift (TALOS+) showed that the Shaker-VSD shares many structural features with the homologous Kv1.2/2.1 chimera, including a transition from α-helix to 3₁₀ helix in the C-terminal portion of the fourth transmembrane helix. Nevertheless, there are clear differences between the Shaker-VSD and Kv1.2/2.1 chimera in the S2-S3 linker and S3 transmembrane region, where the organization of secondary structure elements in Shaker-VSD appears to more closely resemble the KvAP-VSD. Comparison of microsecond-long molecular dynamics simulations of Kv 1.2-VSD in LPPG micelles and a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayer showed that LPPG micelles do not induce significant structural distortion in the isolated voltage sensor. To assess the integrity of the tertiary fold, we directly probed the binding of BrMT analog 2-[2-({[3-(2-amino-ethyl)-6-bromo-1H-indol-2-yl]methoxy}k7methyl)-6-bromo-1H-indol-3-yl]ethan-1-amine (BrET), a gating modifier toxin, and identified the location of the putative binding site. Our results suggest that the Shaker-VSD in LPPG micelles is in a native-like fold and is likely to provide valuable insights into the dynamics of voltage-gating and its regulation.

A single NaK channel conformation is not enough for non-selective ion conduction

NaK and other non-selective channels are able to conduct both sodium (Na⁺) and potassium (K⁺) with equally high efficiency. In contrast to previous crystallographic results, we show that the selectivity filter (SF) of NaK in native-like lipid membranes adopts two distinct conformations that are stabilized by either Na⁺ or K⁺ ions. The atomic differences of these conformations are resolved by solid-state NMR (ssNMR) spectroscopy and molecular dynamics (MD) simulations. Besides the canonical K⁺ permeation pathway, we identify a side entry ion-conduction pathway for Na⁺ permeation unique to NaK. Moreover, under otherwise identical conditions ssNMR spectra of the K⁺ selective NaK mutant (NaK2K) reveal only a single conformational state. Therefore, we propose that structural plasticity within the SF and the selection of these conformations by different ions are key molecular determinants for highly efficient conduction of different ions in non-selective cation channels.

How to Connect Cardiac Excitation to the Atomic Interactions of Ion Channels

Many have worked to create cardiac action potential models that explicitly represent atomic-level details of ion channel structure. Such models have the potential to define new therapeutic directions and to show how nanoscale perturbations to channel function predispose patients to deadly cardiac arrhythmia. However, there have been significant experimental and theoretical barriers that have limited model usefulness. Recently, many of these barriers have come down, suggesting that considerable progress toward creating these long-sought models may be possible in the near term.

Mechanisms and models of cardiac sodium channel inactivation

Shortly after cardiac Na+ channels activate and initiate the action potential, inactivation ensues within milliseconds, attenuating the peak Na+ current, INa, and allowing the cell membrane to repolarize. A very limited number of Na+ channels that do not inactivate carry a persistent INa, or late INa. While late INa is only a small fraction of peak magnitude, it significantly prolongs ventricular action potential duration, which predisposes patients to arrhythmia. Here, we review our current understanding of inactivation mechanisms, their regulation, and how they have been modeled computationally. Based on this body of work, we conclude that inactivation and its connection to late INa would be best modeled with a "feet-on-the-door" approach where multiple channel components participate in determining inactivation and late INa. This model reflects experimental findings showing that perturbation of many channel locations can destabilize inactivation and cause pathological late INa.

A computationally efficient algorithm for fitting ion channel parameters

Continuous time Markov models have been widely used to describe ion channel kinetics, providing explicit representation of channel states and transitions. Fitting models to experimental data remains a computationally demanding task largely due to the high cost of model evaluation. Here, we propose a method to efficiently optimize model parameters and structure. Voltage clamp channel protocols can be decomposed into a series of fixed steps of constant voltage resulting in a set of linear systems of differential equations. Given the linear systems, ODE integration can be swapped for the faster matrix exponential routine. With our parallelized implementation, optimized models are able to reproduce a wide range of experimentally collected data within one minute, a 50 times speedup over ODE integration. •The cost of the objective function is reduced by employing the matrix exponential•The likelihood of convergence is improved by applying synchronous start simulated annealing•The approach was tested by optimizing parameters for a model of the cardiac voltage-gated Na⁺ channel, Na_V1.5, and the KCNQ1 K⁺ channel.