Gene variant effects across sodium channelopathies predict function and guide precision therapy

Pathogenic variants in the voltage-gated sodium channel gene family lead to early onset epilepsies, neurodevelopmental disorders, skeletal muscle channelopathies, peripheral neuropathies and cardiac arrhythmias. Disease-associated variants have diverse functional effects ranging from complete loss-of-function to marked gain-of-function. Therapeutic strategy is likely to depend on functional effect. Experimental studies offer important insights into channel function but are resource intensive and only performed in a minority of cases. Given the evolutionarily conserved nature of the sodium channel genes, we investigated whether similarities in biophysical properties between different voltage-gated sodium channels can predict function and inform precision treatment across sodium channelopathies. We performed a systematic literature search identifying functionally assessed variants in any of the nine voltage-gated sodium channel genes until 28 April 2021. We included missense variants that had been electrophysiologically characterized in mammalian cells in whole-cell patch-clamp recordings. We performed an alignment of linear protein sequences of all sodium channel genes and correlated variants by their overall functional effect on biophysical properties. Of 951 identified records, 437 sodium channel-variants met our inclusion criteria and were reviewed for functional properties. Of these, 141 variants were epilepsy-associated (SCN1/2/3/8A), 79 had a neuromuscular phenotype (SCN4/9/10/11A), 149 were associated with a cardiac phenotype (SCN5/10A) and 68 (16%) were considered benign. We detected 38 missense variant pairs with an identical disease-associated variant in a different sodium channel gene. Thirty-five out of 38 of those pairs resulted in similar functional consequences, indicating up to 92% biophysical agreement between corresponding sodium channel variants (odds ratio = 11.3; 95% confidence interval = 2.8 to 66.9; P < 0.001). Pathogenic missense variants were clustered in specific functional domains, whereas population variants were significantly more frequent across non-conserved domains (odds ratio = 18.6; 95% confidence interval = 10.9-34.4; P < 0.001). Pore-loop regions were frequently associated with loss-of-function variants, whereas inactivation sites were associated with gain-of-function (odds ratio = 42.1, 95% confidence interval = 14.5-122.4; P < 0.001), whilst variants occurring in voltage-sensing regions comprised a range of gain- and loss-of-function effects. Our findings suggest that biophysical characterisation of variants in one SCN-gene can predict channel function across different SCN-genes where experimental data are not available. The collected data represent the first gain- versus loss-of-function topological map of SCN proteins indicating shared patterns of biophysical effects aiding variant analysis and guiding precision therapy. We integrated our findings into a free online webtool to facilitate functional sodium channel gene variant interpretation (http://SCN-viewer.broadinstitute.org).Peer reviewe

Fast automated placement of polar hydrogen atoms in protein-ligand complexes

<p>Abstract</p> <p>Background</p> <p>Hydrogen bonds play a major role in the stabilization of protein-ligand complexes. The ability of a functional group to form them depends on the position of its hydrogen atoms. An accurate knowledge of the positions of hydrogen atoms in proteins is therefore important to correctly identify hydrogen bonds and their properties. The high mobility of hydrogen atoms introduces several degrees of freedom: Tautomeric states, where a hydrogen atom alters its binding partner, torsional changes where the position of the hydrogen atom is rotated around the last heavy-atom bond in a residue, and protonation states, where the number of hydrogen atoms at a functional group may change. Also, side-chain flips in glutamine and asparagine and histidine residues, which are common crystallographic ambiguities must be identified before structure-based calculations can be conducted.</p> <p>Results</p> <p>We have implemented a method to determine the most probable hydrogen atom positions in a given protein-ligand complex. Optimality of hydrogen bond geometries is determined by an empirical scoring function which is used in molecular docking. This allows to evaluate protein-ligand interactions with an established model. Also, our method allows to resolve common crystallographic ambiguities such as as flipped amide groups and histidine residues. To ensure high speed, we make use of a dynamic programming approach.</p> <p>Conclusion</p> <p>Our results were checked against selected high-resolution structures from an external dataset, for which the positions of the hydrogen atoms have been validated manually. The quality of our results is comparable to that of other programs, with the advantage of being fast enough to be applied on-the-fly for interactive usage or during score evaluation.</p

Synthesis of 5-Hydroxyectoine from Ectoine: Crystal Structure of the Non-Heme Iron(II) and 2-Oxoglutarate-Dependent Dioxygenase EctD

As a response to high osmolality, many microorganisms synthesize various types of compatible solutes. These organic osmolytes aid in offsetting the detrimental effects of low water activity on cell physiology. One of these compatible solutes is ectoine. A sub-group of the ectoine producer's enzymatically convert this tetrahydropyrimidine into a hydroxylated derivative, 5-hydroxyectoine. This compound also functions as an effective osmostress protectant and compatible solute but it possesses properties that differ in several aspects from those of ectoine. The enzyme responsible for ectoine hydroxylation (EctD) is a member of the non-heme iron(II)-containing and 2-oxoglutarate-dependent dioxygenases (EC 1.14.11). These enzymes couple the decarboxylation of 2-oxoglutarate with the formation of a high-energy ferryl-oxo intermediate to catalyze the oxidation of the bound organic substrate. We report here the crystal structure of the ectoine hydroxylase EctD from the moderate halophile Virgibacillus salexigens in complex with Fe3+ at a resolution of 1.85 Å. Like other non-heme iron(II) and 2-oxoglutarate dependent dioxygenases, the core of the EctD structure consists of a double-stranded β-helix forming the main portion of the active-site of the enzyme. The positioning of the iron ligand in the active-site of EctD is mediated by an evolutionarily conserved 2-His-1-carboxylate iron-binding motif. The side chains of the three residues forming this iron-binding site protrude into a deep cavity in the EctD structure that also harbours the 2-oxoglutarate co-substrate-binding site. Database searches revealed a widespread occurrence of EctD-type proteins in members of the Bacteria but only in a single representative of the Archaea, the marine crenarchaeon Nitrosopumilus maritimus. The EctD crystal structure reported here can serve as a template to guide further biochemical and structural studies of this biotechnologically interesting enzyme family

Precision and accuracy of single-molecule FRET measurements - a multi-laboratory benchmark study

Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards to ensure the reproducibility and accuracy of measurements of FRET efficiencies are currently lacking. Here we report the results of a comparative blind study in which 20 labs determined the FRET efficiencies (E) of several dye-labeled DNA duplexes. Using a unified, straightforward method, we obtained FRET efficiencies with s.d. between ±0.02 and ±0.05. We suggest experimental and computational procedures for converting FRET efficiencies into accurate distances, and discuss potential uncertainties in the experiment and the modeling. Our quantitative assessment of the reproducibility of intensity-based smFRET measurements and a unified correction procedure represents an important step toward the validation of distance networks, with the ultimate aim of achieving reliable structural models of biomolecular systems by smFRET-based hybrid methods

Gene variant effects across sodium channelopathies predict function and guide precision therapy

Pathogenic variants in the voltage-gated sodium channel gene family lead to early onset epilepsies, neurodevelopmental disorders, skeletal muscle channelopathies, peripheral neuropathies and cardiac arrhythmias. Disease-associated variants have diverse functional effects ranging from complete loss-of-function to marked gain-of-function. Therapeutic strategy is likely to depend on functional effect. Experimental studies offer important insights into channel function but are resource intensive and only performed in a minority of cases. Given the evolutionarily conserved nature of the sodium channel genes, we investigated whether similarities in biophysical properties between different voltage-gated sodium channels can predict function and inform precision treatment across sodium channelopathies. We performed a systematic literature search identifying functionally assessed variants in any of the nine voltage-gated sodium channel genes until 28 April 2021. We included missense variants that had been electrophysiologically characterized in mammalian cells in whole-cell patch-clamp recordings. We performed an alignment of linear protein sequences of all sodium channel genes and correlated variants by their overall functional effect on biophysical properties. Of 951 identified records, 437 sodium channel-variants met our inclusion criteria and were reviewed for functional properties. Of these, 141 variants were epilepsy-associated (SCN1/2/3/8A), 79 had a neuromuscular phenotype (SCN4/9/10/11A), 149 were associated with a cardiac phenotype (SCN5/10A) and 68 (16%) were considered benign. We detected 38 missense variant pairs with an identical disease-associated variant in a different sodium channel gene. Thirty-five out of 38 of those pairs resulted in similar functional consequences, indicating up to 92% biophysical agreement between corresponding sodium channel variants (odds ratio = 11.3; 95% confidence interval = 2.8 to 66.9; P < 0.001). Pathogenic missense variants were clustered in specific functional domains, whereas population variants were significantly more frequent across non-conserved domains (odds ratio = 18.6; 95% confidence interval = 10.9-34.4; P < 0.001). Pore-loop regions were frequently associated with loss-of-function variants, whereas inactivation sites were associated with gain-of-function (odds ratio = 42.1, 95% confidence interval = 14.5-122.4; P < 0.001), whilst variants occurring in voltage-sensing regions comprised a range of gain- and loss-of-function effects. Our findings suggest that biophysical characterisation of variants in one SCN-gene can predict channel function across different SCN-genes where experimental data are not available. The collected data represent the first gain- versus loss-of-function topological map of SCN proteins indicating shared patterns of biophysical effects aiding variant analysis and guiding precision therapy. We integrated our findings into a free online webtool to facilitate functional sodium channel gene variant interpretation (http://SCN-viewer.broadinstitute.org).Peer reviewe

Phosphorylation of residues inside the SNARE complex suppresses secretory vesicle fusion

Membrane fusion is essential for eukaryotic life, requiring SNARE proteins to zipper up in an α-helical bundle to pull two membranes together. Here, we show that vesicle fusion can be suppressed by phosphorylation of core conserved residues inside the SNARE domain. We took a proteomics approach using a PKCB knockout mast cell model and found that the key mast cell secretory protein VAMP8 becomes phosphorylated by PKC at multiple residues in the SNARE domain. Our data suggest that VAMP8 phosphorylation reduces vesicle fusion in vitro and suppresses secretion in living cells, allowing vesicles to dock but preventing fusion with the plasma membrane. Markedly, we show that the phosphorylation motif is absent in all eukaryotic neuronal VAMPs, but present in all other VAMPs. Thus, phosphorylation of SNARE domains is a general mechanism to restrict how much cells secrete, opening the door for new therapeutic strategies for suppression of secretion