Hydropathicity-based prediction of pain-causing NaV1.7 variants

Anxiety and dysautonomia symptoms in patients with a NaV1.7 mutation and the potential benefits of low-dose short-acting guanfacine

PURPOSE

Guanfacine is an α2A-adrenergic receptor agonist, FDA-approved to treat attention-deficit hyperactivity disorder and high blood pressure, typically as an extended-release formulation up to 7 mg/day. In our dysautonomia clinic, we observed that off-label use of short-acting guanfacine at 1 mg/day facilitated symptom relief in two families with multiple members presenting with severe generalized anxiety. We also noted anecdotal improvements in associated dysautonomia symptoms such as hyperhidrosis, cognitive impairment, and palpitations. We postulated that a genetic deficit existed in these patients that might augment guanfacine susceptibility.

METHODS

We used whole-exome sequencing to identify mutations in patients with shared generalized anxiety and dysautonomia symptoms. Guanfacine-induced changes in the function of voltage-gated Na+ channels were investigated using voltage-clamp electrophysiology.

RESULTS

Whole-exome sequencing uncovered the p.I739V mutation in SCN9A in the proband of two nonrelated families. Moreover, guanfacine inhibited ionic currents evoked by wild-type and mutant NaV1.7 encoded by SCN9A, as well as other NaV channel subtypes to a varying degree.

CONCLUSION

Our study provides further evidence for a possible pathophysiological role of NaV1.7 in anxiety and dysautonomia. Combined with off-target effects on NaV channel function, daily administration of 1 mg short-acting guanfacine may be sufficient to normalize NaV channel mutation-induced changes in sympathetic activity, perhaps aided by partial inhibition of NaV1.7 or other channel subtypes. In a broader context, expanding genetic and functional data about ion channel aberrations may enable the prospect of stratifying patients in which mutation-induced increased sympathetic tone normalization by guanfacine can support treatment strategies for anxiety and dysautonomia symptoms.

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. As a consequence, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion-channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated. We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion-channel families. We collected and curated 3049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12 546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures. We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis centre and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain versus loss-of-channel function. In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional readouts and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.