Diverse biophysical mechanisms in voltage-gated sodium channel Nav1.4 variants associated with myotonia.

Mutations in SCN4A gene encoding Nav1.4 channel α-subunit, are known to cause neuromuscular disorders such as myotonia or paralysis. Here, we study the effect of two amino acid replacements, K1302Q and G1306E, in the DIII-IV loop of the channel, corresponding to mutations found in patients with myotonia. We combine clinical, electrophysiological, and molecular modeling data to provide a holistic picture of the molecular mechanisms operating in mutant channels and eventually leading to pathology. We analyze the existing clinical data for patients with the K1302Q substitution, which was reported for adults with or without myotonia phenotypes, and report two new unrelated patients with the G1306E substitution, who presented with severe neonatal episodic laryngospasm and childhood-onset myotonia. We provide a functional analysis of the mutant channels by expressing Nav1.4 α-subunit in Xenopus oocytes in combination with ß1 subunit and recording sodium currents using two-electrode voltage clamp. The K1302Q variant exhibits abnormal voltage dependence of steady-state fast inactivation, being the likely cause of pathology. K1302Q does not lead to decelerated fast inactivation, unlike several other myotonic mutations such as G1306E. For both mutants, we observe increased window currents corresponding to a larger population of channels available for activation. To elaborate the structural rationale for our experimental data, we explore the contacts involving K/Q1302 and E1306 in the AlphaFold2 model of wild-type Nav1.4 and Monte Carlo-minimized models of mutant channels. Our data provide the missing evidence to support the classification of K1302Q variant as likely pathogenic and may be used by clinicians.

Biophysical mechanisms of myocardium sodium channelopathies.

Genetic variants of gene SCN5A encoding the alpha-subunit of cardiac voltage-gated sodium channel Nav1.5 are associated with various diseases, including long QT syndrome (LQT3), Brugada syndrome (BrS1), and progressive cardiac conduction disease (PCCD). In the last decades, the great progress in understanding molecular and biophysical mechanisms of these diseases has been achieved. The LQT3 syndrome is associated with gain-of-function of sodium channels Nav1.5 due to impaired inactivation, enhanced activation, accelerated recovery from inactivation or the late current appearance. In contrast, BrS1 and PCCD are associated with the Nav1.5 loss-of-function, which in electrophysiological experiments can be manifested as reduced current density, enhanced fast or slow inactivation, impaired activation, or decelerated recovery from inactivation. Genetic variants associated with congenital arrhythmias can also disturb interactions of the Nav1.5 channel with different proteins or drugs and cause unexpected reactions to drug administration. Furthermore, mutations can affect post-translational modifications of the channels and their sensitivity to pH and temperature. Here we briefly review the current knowledge on biophysical mechanisms of LQT3, BrS1 and PCCD. We focus on limitations of studies that use heterologous expression systems and induced pluripotent stem cells (iPSC) derived cardiac myocytes and summarize our understanding of genotype-phenotype relations of SCN5A mutations.

A lysine residue from an extracellular turret switches the ion preference in a Cav3 T-Type channel from calcium to sodium ions.

Cav3 T-type calcium channels from great pond snail Lymnaea stagnalis have a selectivity-filter ring of five acidic residues, EE(D)DD. Splice variants with exons 12b or 12a spanning the extracellular loop between the outer helix IIS5 and membrane-descending pore helix IIP1 (IIS5-P1) in Domain II of the pore module possess calcium selectivity or dominant sodium permeability, respectively. Here, we use AlphaFold2 neural network software to predict that a lysine residue in exon 12a is salt-bridged to the aspartate residue immediately C terminal to the second-domain glutamate in the selectivity filter. Exon 12b has a similar folding but with an alanine residue in place of lysine in exon 12a. We express LCav3 channels with mutated exons Ala-12b-Lys and Lys-12a-Ala and demonstrate that they switch the ion preference to high sodium permeability and calcium selectivity, respectively. We propose that in the calcium-selective variants, a calcium ion chelated between Domain II selectivity-filter glutamate and aspartate is knocked-out by the incoming calcium ion in the process of calcium permeation, whereas sodium ions are repelled. The aspartate is neutralized by the lysine residue in the sodium-permeant variants, allowing for sodium permeation through the selectivity-filter ring of four negatively charged residues akin to the prokaryotic sodium channels with four glutamates in the selectivity filter. The evolutionary adaptation in invertebrate LCav3 channels highlight the involvement of a key, ubiquitous aspartate, "a calcium beacon" of sorts in the outer pore of Domain II, as determinative for the calcium ion preference over sodium ions through eukaryotic Cav1, Cav2, and Cav3 channels.

Mutations in transmembrane proteins: diseases, evolutionary insights, prediction and comparison with globular proteins.

Membrane proteins are unique in that they interact with lipid bilayers, making them indispensable for transporting molecules and relaying signals between and across cells. Due to the significance of the protein's functions, mutations often have profound effects on the fitness of the host. This is apparent both from experimental studies, which implicated numerous missense variants in diseases, as well as from evolutionary signals that allow elucidating the physicochemical constraints that intermembrane and aqueous environments bring. In this review, we report on the current state of knowledge acquired on missense variants (referred to as to single amino acid variants) affecting membrane proteins as well as the insights that can be extrapolated from data already available. This includes an overview of the annotations for membrane protein variants that have been collated within databases dedicated to the topic, bioinformatics approaches that leverage evolutionary information in order to shed light on previously uncharacterized membrane protein structures or interaction interfaces, tools for predicting the effects of mutations tailored specifically towards the characteristics of membrane proteins as well as two clinically relevant case studies explaining the implications of mutated membrane proteins in cancer and cardiomyopathy.

Brevetoxin versus Brevenal Modulation of Human Nav1 Channels.

Brevetoxins (PbTx) and brevenal are marine ladder-frame polyethers. PbTx binds to and activates voltage-gated sodium (Nav) channels in native tissues, whereas brevenal antagonizes these actions. However, the effects of PbTx and brevenal on recombinant Nav channel function have not been systematically analyzed. In this study, the PbTx-3 and brevenal modulation of tissue-representative Nav channel subtypes Nav1.2, Nav1.4, Nav1.5, and Nav1.7 were examined using automated patch-clamp. While PbTx-3 and brevenal elicit concentration-dependent and subtype-specific modulatory effects, PbTx-3 is >1000-fold more potent than brevenal. Consistent with effects observed in native tissues, Nav1.2 and Nav1.4 channels were PbTx-3- and brevenal-sensitive, whereas Nav1.5 and Nav1.7 appeared resistant. Interestingly, the incorporation of brevenal in the intracellular solution caused Nav channels to become less sensitive to PbTx-3 actions. Furthermore, we generated a computational model of PbTx-2 bound to the lipid-exposed side of the interface between domains I and IV of Nav1.2. Our results are consistent with competitive antagonism between brevetoxins and brevenal, setting a basis for future mutational analyses of Nav channels' interaction with brevetoxins and brevenal. Our findings provide valuable insights into the functional modulation of Nav channels by brevetoxins and brevenal, which may have implications for the development of new Nav channel modulators with potential therapeutic applications.

Functional Analysis of SCN5A Genetic Variants Associated with Brugada Syndrome.

BACKGROUND: Brugada syndrome (BrS) is a rare inherited cardiac arrhythmia with increased risk of sudden cardiac death. Mutations in gene SCN5A, which encodes the α-subunit of cardiac voltage-gated sodium channel NaV1.5, have been identified in over 20% of patients with BrS. However, only a small fraction of NaV1.5 variants, which are associated with BrS, are characterized in electrophysiological experiments. RESULTS: Here we explored variants V281A and L1582P, which were found in our patients with BrS, and variants F543L and K1419E, which are reportedly associated with BrS. Heterologous expression of the variants in CHO-K1 cells and the Western blot analysis demonstrated that each variant appeared at the cell surface. We further measured sodium current in the whole-cell voltage clamp configuration. Variant F543L produced robust sodium current with a hyperpolarizing shift in the voltage dependence of steady-state fast inactivation. Other variants did not produce detectable sodium currents, indicating a complete loss of function. In a recent cryoEM structure of the hNaV1.5 channel, residues V281, K1419, and L1582 are in close contacts with residues whose mutations are reportedly associated with BrS, indicating functional importance of respective contacts. CONCLUSIONS: Our results support the notion that loss of function of NaV1.5 or decrease of the channel activity is involved in the pathogenesis of BrS.

Mapping the interaction surface of scorpion ß-toxins with an insect sodium channel.

The interaction of insect-selective scorpion depressant ß-toxins (LqhIT2 and Lqh-dprIT3 from Leiurus quinquestriatus hebraeus) with the Blattella germanica sodium channel, BgNav1-1a, was investigated using site-directed mutagenesis, electrophysiological analyses, and structural modeling. Focusing on the pharmacologically defined binding site-4 of scorpion ß-toxins at the voltage-sensing domain II (VSD-II), we found that charge neutralization of D802 in VSD-II greatly enhanced the channel sensitivity to Lqh-dprIT3. This was consistent with the high sensitivity of the splice variant BgNav2-1, bearing G802, to Lqh-dprIT3, and low sensitivity of BgNav2-1 mutant, G802D, to the toxin. Further mutational and electrophysiological analyses revealed that the sensitivity of the WT = D802E < D802G < D802A < D802K channel mutants to Lqh-dprIT3 correlated with the depolarizing shifts of activation in toxin-free channels. However, the sensitivity of single mutants involving IIS4 basic residues (K4E = WT << R1E < R2E < R3E) or double mutants (D802K = K4E/D802K = R3E/D802K > R2E/D802K > R1E/D802K > WT) did not correlate with the activation shifts. Using the cryo-EM structure of the Periplaneta americana channel, NavPaS, as a template and the crystal structure of LqhIT2, we constructed structural models of LqhIT2 and Lqh-dprIT3-c in complex with BgNav1-1a. These models along with the mutational analysis suggest that depressant toxins approach the salt-bridge between R1 and D802 at VSD-II to form contacts with linkers IIS1-S2, IIS3-S4, IIIP5-P1 and IIIP2-S6. Elimination of this salt-bridge enables deeper penetration of the toxin into a VSD-II gorge to form new contacts with the channel, leading to increased channel sensitivity to Lqh-dprIT3.

Computational Analysis of the Crystal and Cryo-EM Structures of P-Loop Channels with Drugs.

The superfamily of P-loop channels includes various potassium channels, voltage-gated sodium and calcium channels, transient receptor potential channels, and ionotropic glutamate receptors. Despite huge structural and functional diversity of the channels, their pore-forming domain has a conserved folding. In the past two decades, scores of atomic-scale structures of P-loop channels with medically important drugs in the inner pore have been published. High structural diversity of these complexes complicates the comparative analysis of these structures. Here we 3D-aligned structures of drug-bound P-loop channels, compared their geometric characteristics, and analyzed the energetics of ligand-channel interactions. In the superimposed structures drugs occupy most of the sterically available space in the inner pore and subunit/repeat interfaces. Cationic groups of some drugs occupy vacant binding sites of permeant ions in the inner pore and selectivity-filter region. Various electroneutral drugs, lipids, and detergent molecules are seen in the interfaces between subunits/repeats. In many structures the drugs strongly interact with lipid and detergent molecules, but physiological relevance of such interactions is unclear. Some eukaryotic sodium and calcium channels have state-dependent or drug-induced π-bulges in the inner helices, which would be difficult to predict. The drug-induced π-bulges may represent a novel mechanism of gating modulation.

Predicting novel disease mutations in the cardiac sodium channel.

BACKGROUND: Voltage-gated sodium channels Nav1.x mediate the rising phase of action potential in excitable cells. Variations in gene SCN5A, which encodes the hNav1.5 channel, are associated with arrhythmias and other heart diseases. About 1,400 SCN5A variants are listed in public databases, but for more than 30% of these the clinical significance is unknown and can currently only be derived by bioinformatics approaches. METHODS AND RESULTS: We used the ClinVar, SwissVar, Humsavar, gnomAD, and Ensembl databases to assemble a dataset of 1392 hNav1.5 variants (370 pathogenic variants, 602 benign variants and 420 variants of uncertain significance) as well as a dataset of 1766 damaging variants in 20 human sodium and calcium channel paralogs. Twelve in silico tools were tested for their ability to predict damaging mutations in hNav1.5. The best performing tool, MutPred, correctly predicted 93% of damaging variants in our hNav1.5 dataset. Among the 86 hNav1.5 variants for which electrophysiological data are also available, MutPred correctly predicted 82% of damaging variants. In the subset of 420 uncharacterized hNav1.5 variants MutPred predicted 196 new pathogenic variants. Among these, 74 variants are also annotated as damaging in at least one hNav1.5 paralog. CONCLUSIONS: Using a combination of sequence-based bioinformatics techniques and paralogous annotation we have substantially expanded the knowledge on disease variants in the cardiac sodium channel and assigned a pathogenic status to a number of mutations that so far have been described as variants of uncertain significance. A list of reclassified hNav1.5 variants and their properties is provided.

Optical control of L-type Ca2+ channels using a diltiazem photoswitch.

L-type Ca2+ channels (LTCCs) play a crucial role in excitation-contraction coupling and release of hormones from secretory cells. They are targets of antihypertensive and antiarrhythmic drugs such as diltiazem. Here, we present a photoswitchable diltiazem, FHU-779, which can be used to reversibly block endogenous LTCCs by light. FHU-779 is as potent as diltiazem and can be used to place pancreatic ß-cell function and cardiac activity under optical control.

Phenylalkylamines in calcium channels: computational analysis of experimental structures.

Experimental 3D structures of calcium channels with phenylalkylamines (PAAs) provide basis for further analysis of atomic mechanisms of these important cardiovascular drugs. In the crystal structure of the engineered calcium channel CavAb with Br-verapamil and in the cryo-EM structure of the Cav1.1 channel with verapamil, the ligands bind in the inner pore. However, there are significant differences between these structures. In the crystal structure the ligand ammonium group is much closer to the ion in the selectivity-filter region Site 3, which is most proximal to the inner pore, than in the cryo-EM structure. Here we used Monte Carlo energy minimizations to dock PAAs in calcium channels. Our computations suggest that in the crystal structure Site 3 is occupied by a water molecule rather than by a calcium ion. Analysis of the published electron density map does not rule out this possibility. In the cryo-EM structures the ammonium group of verapamil is shifted from the calcium ion in Site 3 either along the pore axis, towards the cytoplasm or away from the axis. Our unbiased docking reproduced these binding modes. However, in the cryo-EM structures detergent and lipid molecules interact with verapamil. When we removed these molecules, the nitrile group of verapamil bound to the calcium ion in Site 3. Models of Cav1.2 with different PAAs suggest similar binding modes and direct contacts of the ligands electronegative atoms with the calcium ion in Site 3. Such interactions explain paradoxes in structure-activity relationships of PAAs.

Chronology of sodium channel mutations associated with pyrethroid resistance in Aedes aegypti.

Aedes aegypti is the primary mosquito vector of dengue, yellow fever, Zika and chikungunya. Current strategies to control Ae. aegypti rely heavily on insecticide interventions. Pyrethroids are a major class of insecticides used for mosquito control because of their fast acting, highly insecticidal activities and low mammalian toxicity. However, Ae. aegypti populations around the world have begun to develop resistance to pyrethroids. So far, more than a dozen mutations in the sodium channel gene have been reported to be associated with pyrethroid resistance in Ae. aegypti. Co-occurrence of resistance-associated mutations is common in pyrethroid-resistant Ae. aegypti populations. As global use of pyrethroids in mosquito control continues, new pyrethroid-resistant mutations keep emerging. In this microreview, we compile pyrethroid resistance-associated mutations in Ae. aegypti in a chronological order, as they were reported, and summarize findings from functional evaluation of these mutations in an in vitro sodium channel expression system. We hope that the information will be useful for tracing possible evolution of pyrethroid resistance in this important human disease vector, in addition to the development of methods for global monitoring and management of pyrethroid resistance in Ae. aegypti.

Molecular basis of selective resistance of the bumblebee BiNav1 sodium channel to tau-fluvalinate.

Insecticides are widely used to control pests in agriculture and insect vectors that transmit human diseases. However, these chemicals can have a negative effect on nontarget, beneficial organisms including bees. Discovery and deployment of selective insecticides is a major mission of modern toxicology and pest management. Pyrethroids exert their toxic action by acting on insect voltage-gated sodium channels. Honeybees and bumblebees are highly sensitive to most pyrethroids, but are resistant to a particular pyrethroid, tau-fluvalinate (τ-FVL). Because of its unique selectivity, τ-FVL is widely used to control not only agricultural pests but also varroa mites, the principal ectoparasite of honeybees. However, the mechanism of bee resistance to τ-FVL largely remains elusive. In this study, we functionally characterized the sodium channel BiNav1-1 from the common eastern bumblebee (Bombus impatiens) in Xenopus oocytes and found that the BiNav1-1 channel is highly sensitive to six commonly used pyrethroids, but resistant to τ-FVL. Phylogenetic and mutational analyses revealed that three residues, which are conserved in sodium channels from 12 bee species, underlie resistance to τ-FVL or sensitivity to the other pyrethroids. Further computer modeling and mutagenesis uncovered four additional residues in the pyrethroid receptor sites that contribute to the unique selectivity of the bumblebee sodium channel to τ-FVL versus other pyrethroids. Our data contribute to understanding a long-standing enigma of selective pyrethroid toxicity in bees and may be used to guide future modification of pyrethroids to achieve highly selective control of pests with minimal effects on nontarget organisms.

Characterization of a novel SCN5A genetic variant A1294G associated with mixed clinical phenotype.

Mutations in gene SCN5A, which encodes cardiac voltage-gated sodium channel Nav1.5, are associated with multiple clinical phenotypes. Here we describe a novel A1294G genetic variant detected in a male patient with combined clinical phenotype including atrioventricular II block, Brugada-like ECG, septal fibrosis, right ventricular dilatation and decreased left ventricular contractility. Residue A1294 is located in the IIIS3-S4 extracellular loop, in proximity to several residues whose mutations are associated with sodium channelopathies. The wild-type channel Nav1.5 and mutant Nav1.5-A1294G were expressed in the CHO-K1 and HEK293T cells and whole-cell sodium currents were recorded using the patch-clamp method. The A1294G channels demonstrated a negative shift of steady-state inactivation, accelerated fast and slow inactivation and decelerated recovery from intermediate inactivation. Our study reveals biophysical mechanism of the Nav1.5-A1294G dysfunction, which may underlie the combined phenotypic manifestation observed in the patient.

Extremely Potent Block of Bacterial Voltage-Gated Sodium Channels by µ-Conotoxin PIIIA.

µ-Conotoxin PIIIA, in the sub-picomolar, range inhibits the archetypal bacterial sodium channel NaChBac (NavBh) in a voltage- and use-dependent manner. Peptide µ-conotoxins were first recognized as potent components of the venoms of fish-hunting cone snails that selectively inhibit voltage-gated skeletal muscle sodium channels, thus preventing muscle contraction. Intriguingly, computer simulations predicted that PIIIA binds to prokaryotic channel NavAb with much higher affinity than to fish (and other vertebrates) skeletal muscle sodium channel (Nav 1.4). Here, using whole-cell voltage clamp, we demonstrate that PIIIA inhibits NavBac mediated currents even more potently than predicted. From concentration-response data, with [PIIIA] varying more than 6 orders of magnitude (10-12 to 10-5 M), we estimated an IC50 = ~5 pM, maximal block of 0.95 and a Hill coefficient of 0.81 for the inhibition of peak currents. Inhibition was stronger at depolarized holding potentials and was modulated by the frequency and duration of the stimulation pulses. An important feature of the PIIIA action was acceleration of macroscopic inactivation. Docking of PIIIA in a NaChBac (NavBh) model revealed two interconvertible binding modes. In one mode, PIIIA sterically and electrostatically blocks the permeation pathway. In a second mode, apparent stabilization of the inactivated state was achieved by PIIIA binding between P2 helices and trans-membrane S5s from adjacent channel subunits, partially occluding the outer pore. Together, our experimental and computational results suggest that, besides blocking the channel-mediated currents by directly occluding the conducting pathway, PIIIA may also change the relative populations of conducting (activated) and non-conducting (inactivated) states.

Mutational analysis of state-dependent contacts in the pore module of eukaryotic sodium channels.

Voltage-gated sodium channels have residues that change or may change contacts upon gating. Contributions of individual contacts in stability of different states are incompletely understood. Pore-lining inner helices contain exceptionally conserved asparagines in positions i20. Here we explored how mutations in positions i20 and i29 affect electrophysiological properties of insect sodium channels. In repeat interfaces I/IV, III/II and IV/III, alanine substitutions caused positive activation shifts in positions i20 and i29, negative shifts of slow inactivation in positions i20 and positive shifts of slow inactivation in positions i29. The results support the hypothesis on open state inter-repeat H-bonding of residues i20 and i29. The shift magnitudes vary between interfaces, reflecting structural asymmetry of the channels. Mutations in positions i20 of repeats III and IV caused much longer recovery delay from the slow and fast inactivation than other mutations. In repeat IV, alanine substitution of tyrosine i30 rescued positive activation shift of mutation in position i29. Our data suggest that polar residues in positions i29 are involved in stabilization of both the open and slow-inactivated states. Transition between the states may involve switching of H-bonding partners of residues i29 from the conserved asparagines to currently unknown residues.

Structural Models of Ligand-Bound Sodium Channels.

X-ray and cryo-EM structures of tetrameric and pseudo-tetrameric P-loop channels are used to elaborate homology models of mammalian voltage-gated sodium channels with drugs and neurotoxins. Such models integrate experimental data, assist in planning new experiments, and may facilitate drug design. This chapter outlines sodium channel models with local anesthetics, anticonvulsants, and antiarrhythmics, which are used to manage pain and treat sodium channelopathies. Further summarized are sodium channel models with tetrodotoxin, mu-conotoxins, batrachotoxin, scorpion toxins, and insecticides. Possible involvement of sodium ions in the action of some ligands is discussed.

Action of six pyrethrins purified from the botanical insecticide pyrethrum on cockroach sodium channels expressed in Xenopus oocytes.

Pyrethrin I, pyrethrin II, cinerin I, cinerin II, jasmolin I and jasmolin II are six closely related insecticidal active esters, known as pyrethrins, found in the pyrethrum extract from the dry flowers of Tanacetum cinerariifolium. The chemical structures of the six compounds differ only in the terminal moieties at the acid and alcohol ends, but the compounds' in vivo toxicities are substantially different. Pyrethrins are lead compounds for pyrethroids, a large family of synthetic insecticides that alter nerve functions by prolonging the opening of voltage-gated sodium channels. However, data on the mechanism of action of natural pyrethrins are very limited. In this study, we examined the actions of all six pyrethrins on cockroach sodium channels expressed in Xenopus oocytes. Although the six compounds showed comparable potencies in inhibiting the inactivation of sodium channels, they had greatly variable potencies in inhibiting channel deactivation. Furthermore, unlike pyrethroids, the action of pyrethrins neither depend on nor were enhanced by repeated channel activation. We created a NavMs-based model of the cockroach sodium channel, in which pyrethrin II was docked at the pyrethroid receptor site 1 (PyR1), and proposed a rationale for the observed structure-activity relationship of the six pyrethrins. Our study sheds light on the molecular mechanism of pyrethrum action on sodium channels and reveled differences in the modes of action of the six bioactive constitutes of pyrethrum.

Evidence for Dual Binding Sites for 1,1,1-Trichloro-2,2-bis(p-chlorophenyl)ethane (DDT) in Insect Sodium Channels.

1,1,1-Trichloro-2,2-bis(p-chlorophenyl)ethane (DDT), the first organochlorine insecticide, and pyrethroid insecticides are sodium channel agonists. Although the use of DDT is banned in most of the world due to its detrimental impact on the ecosystem, indoor residual spraying of DDT is still recommended for malaria control in Africa. Development of resistance to DDT and pyrethroids is a serious global obstacle for managing disease vectors. Mapping DDT binding sites is necessary for understanding mechanisms of resistance and modulation of sodium channels by structurally different ligands. The pioneering model of the housefly sodium channel visualized the first receptor for pyrethroids, PyR1, in the II/III domain interface and suggested that DDT binds within PyR1. Previously, we proposed the second pyrethroid receptor, PyR2, at the I/II domain interface. However, whether DDT binds to both pyrethroid receptor sites remains unknown. Here, using computational docking of DDT into the Kv1.2-based mosquito sodium channel model, we predict that two DDT molecules can bind simultaneously within PyR1 and PyR2. The bulky trichloromethyl group of each DDT molecule fits snugly between four helices in the bent domain interface, whereas two p-chlorophenyl rings extend into two wings of the interface. Model-driven mutagenesis and electrophysiological analysis confirmed these propositions and revealed 10 previously unknown DDT-sensing residues within PyR1 and PyR2. Our study proposes a dual DDT-receptor model and provides a structural background for rational development of new insecticides.