Methods

To further elucidate the genotype-phenotype correlation of SCN2A variants we investigated the functional effects of six variants representing the phenotypic spectrum by whole-cell patch-clamp studies in transfected HEK293T cells and in-silico structural modeling.

Results

The two variants p.L1342P and p.E1803G detected in patients with early onset epileptic encephalopathy (EE) showed profound and complex changes in channel gating, whereas the BFNIE variant p.L1563V exhibited only a small gain of channel function. The three variants identified in ID patients without seizures, p.R937C, p.L611Vfs*35 and p.W1716*, did not produce measurable currents. Homology modeling of the missense variants predicted structural impairments consistent with the electrophysiological findings.

Conclusions

Our findings support the hypothesis that complete loss-of-function variants lead to ID without seizures, small gain-of-function variants cause BFNIE and EE variants exhibit variable but profound Nav1.2 gating changes. Moreover, structural modeling was able to predict the severity of the variant impact, supporting a potential role of structural modeling as a prognostic tool. Our study on the functional consequences of SCN2A variants causing the distinct phenotypes of EE, BFNIE and ID contributes to the elucidation of mechanisms underlying the broad phenotypic variability reported for SCN2A variants.

Background

SCN2A (OMIM 182390) encodes the pore-forming α-subunit of the voltage gated sodium channel type 2 (Nav1.2). It is predominantly expressed at the axon initial segment and nodes of Ranvier of excitatory neurons of the central nervous system and plays a crucial role for axon potential initiation and propagation in early development (Kaplan et al. 2001, Kole and Stuart 2012). At age 1–2 years, it is partially replaced by Nav1.6 (encoded by SCN8A) and its main function is then thought to be related to the support of backpropagating action potentials (Hu et al. 2009, Sanders et al. 2018). Variants in SCN2A were initially identified as disease cause in patients with generalized epilepsy with febrile seizures plus (GEFS+) (Sugawara et al. 2001) and benign familial neonatal-infantile epilepsy (BFNIE) (Heron et al. 2002). In recent years, SCN2A has been recognized as one of the most frequently mutated genes in neurodevelopmental disorders including autism spectrum disorders (ASD), intellectual disability (ID), developmental and epileptic encephalopathy (DEE) and schizophrenia (Li et al. 2016).

Here, we investigated the functional consequences of six SCN2A variants identified in patients with ID without seizures, early-onset EE or BFNIE using whole-cell patch-clamp experiments in HEK293T cells transiently expressing wild-type or mutant Nav1.2. Furthermore, the description of homologous structures in 2017 (Yan et al. 2017, Shen et al. 2017) allowed us to generate a 3D-model for SCN2A and assess the impact of the different variants on protein structure. We further re-evaluated the first family reported with a pathogenic SCN2A variant causing DEE (p.R102*) (Kamiya et al. 2004) which stands out as the only functionally studied truncating variant leading to DEE.

Methods

Patients and variants

Five previously reported patients carrying a pathogenic SCN2A variant (Hackenberg et al. 2014, Rauch et al. 2012, Wolff et al. 2017, Papuc et al. 2018) with a phenotype of either early-onset EE (n = 2) or ID without seizures (n = 3) were included in this study for functional and clinical evaluation. One BFNIE variant previously studied extensively by others was also investigated as a reference.

Whole exome sequencing

To re-evaluate the DEE case published by Kamiya et al. in 2004 carrying the p.R102* variant detected by targeted sequencing of SCN2A, we now performed trio whole exome sequencing (WES) of this family (see Additional file 1).

Electrophysiology

Whole-cell voltage-clamp recordings were performed at room temperature 18–30 h after transfection using 2–4 MΩ borosilicate glass pipettes, an EPC7 amplifier and Patchmaster v2.11 software (HEKA Elektronik Dr. Schulze GmbH, Lambrecht, Germany). The pipette solution contained (in mM) 110 CsF, 10 NaF, 20 CsCl, 2 EGTA and 10 HEPES ~ 310 mOsm (pH 7.35 with CsOH). Bath solution consisted of (in mM) 145 NaCl, 4 KCl, 1.8 CaCl2, 1 MgCl2 and 10 HEPES ~ 310 mOsm (pH 7.35 with NaOH). Cells with peak current amplitudes < 0.6 nA were not used for analysis of biophysical parameters due to possible contamination with endogenous currents of HEK cells. Whole-cell capacitance was determined by integrating the area under the capacitive currents in response to a hyperpolarizing voltage step. Capacitive transients were compensated electronically. Leak currents were subtracted by use of a P/N procedure. Sodium currents were evoked by 10 ms voltage steps from − 100 mV to + 50 mV from a holding potential of − 120 mV. Current densities were obtained by dividing the peak currents by the capacitance. Charge transfer was calculated as the area under the current trace and normalized to the capacitance. Sodium conductance (GNa) was calculated as GNa = INa / (Vm − Vrev), where INa is the measured peak current at the test potential Vm, and Vrev is the calculated sodium reversal potential. To obtain activation curves (conductance-voltage relationship), the normalized conductance was plotted against the test potentials and fitted with a Boltzmann sigmoidal function for quantitative analysis. The 10–90% rise time was obtained from the I-V current traces. Fast and slow inactivation time constants were obtained by fitting a second order exponential function to the current decay of the I-V current traces. Steady-state inactivation was determined using 500 ms conditioning prepulses to potentials varying from − 130 to + 10 mV followed by a depolarizing test pulse to + 5 mV. For quantitative analysis, normalized peak currents were plotted against prepulse potential and fitted with a Boltzmann sigmoidal function. Recovery from fast inactivation was recorded by a two-pulse protocol from a holding potential of − 90 mV to + 5 mV with increasing interpulse intervals. For quantitative analyses, the curves were fitted with a first order exponential function where the tau represents the recovery time constant. To analyze use-dependent inactivation, depolarizing pulse trains from − 90 mV to + 10 mV at different frequencies were applied and the residual peak current of the last pulse was normalized to the peak current of the first pulse.

Protein modeling

The structure of SCN2A was modeled with Modeller 9.16 (Webb and Sali 2017) using the recently published structure of a homologous voltage-gated sodium channel (PDB: 5X0M (Shen et al. 2017)) as template. Residues 440–741 and 990–1184 were excluded from the modeling procedure because they were either absent in the template or they did not exhibit a sufficient sequence similarity between template and target sequence. RasMol (Sayle and Milner-White 1995) was used for structure analysis and visualization.

Results

Patient phenotypes and variants

Comprehensive phenotypes of our patients and bioinformatic annotations of their respective variants are presented in Table 1 and their location in the 2D-scheme of the Nav1.2 channel is depicted in Fig. 1. All missense variants affect a highly conserved amino acid and have damaging in-silico predictions.

Schematic drawing of SCN2A indicating the location of studied pathogenic variants. The three patients with ID without seizures carried de novo heterozygous variants in SCN2A (shown in blue), two of which are nonsense variants (p.L611Vfs*35 and p.W1716*) and one is a missense variant in the pore-forming loop of domain 2 of the channel (p.R937C). The two patients with early onset EE harbored de novo heterozygous missense variants (shown in red) in the fifth transmembrane segment of the third channel domain (p.L1342P) or the C-terminal domain (p.E1803G), respectively. The heterozygous inherited missense variant previously reported by others in a family with BFNIE (Heron et al. 2002) (shown in green) is located in the second transmembrane segment of domain 4 (p.L1563V). Adapted from Meisler and Kearney (2005)

Patch-clamp recordings in HEK293T cells expressing recombinant Nav1.2 wild-type or mutant channels. Variants E1803G (blue) and L1342P (green) have been detected in patients with EE, L1563V (red) in a family with BFNIE, and R937C (purple), L611Vfs*35 (orange) and W1716* (light blue) cause ID. a Family of whole-cell sodium currents of wild-type and conducting mutant Nav1.2. b Peak current amplitudes normalized to cell capacitance (bars are mean ± SEM) were not significantly different. c Voltage dependence of activation obtained by plotting the normalized conductance against test potentials and fitted with a Boltzmann function. L1342P showed a significant negative shift (p < 0.001). d Voltage dependence of steady-state inactivation obtained by displaying the normalized current amplitude against the prepulse potential fitted with a Boltzmann function. A negative shift was observed for both L1342P (p < 0.001) and L1563 V (p = 0.017), while E1803G had a shallower slope (p < 0.001). e Activation and inactivation curves of L1342P in comparison to wild-type demonstrating a hyperpolarized shift with unchanged activation-inactivation coupling. f Activation and inactivation curves of E1803G in comparison to wild-type illustrating the increase of the window current. g Time-dependent recovery from inactivation fitted with a single exponential function and plotted against a logarithmic time scale. A significantly slower recovery rate was observed in L1342P (p = 0.002). h Use-dependent inactivation. The frequencies of depolarizing pulse trains are plotted against the residual peak current amplitude of the last pulse normalized to the peak current amplitude of the first pulse. i 10–90% rise time depicted against test potentials. L1342P exhibited a significantly faster rise time at the voltage step to − 30 mV (p = 0.02). k Fast inactivation time constants plotted against test potentials

The EE mutant E1803G exhibits a shallower slope of the inactivation curve that leads to an increased window current

While no changes in voltage dependence of activation were observed for the mutant E1803G compared to wild-type, the slope of the voltage dependence of inactivation was significantly less steep (Fig. 2c-d, Table 2), leading to a larger overlap of activation and inactivation curves. This increased window current suggests the presence of a persistent current at membrane voltages between − 45 to − 25 mV (Fig. 2f). Additionally, the E1803G mutant tended to have slower activation and inactivation kinetics (Fig. 2i-k, Additional file 1: Table S1) and hence tended to let more charge pass compared to wild-type channels (mean charge transfer normalized to conductance at a test potential of − 10 mV: wild-type = 0.09 ± 0.02 C/F, n = 8; E1803G = 0.16 ± 0.03 C/F, n = 8; p = 0.092). The use-dependent inactivation was similar to wild-type (Fig. 2h, Additional file 1: Table S2). The larger window current may increase neuronal excitability and thereby promote the presence of epileptic seizures in affected patients.

The BFNIE mutant L1563V exhibits a small depolarizing shift of the inactivation curves

Activation, recovery from inactivation, kinetics and use-dependent inactivation did not show relevant differences between L1563V and wild-type channels (Fig. 2, Table 2 and Additional file 1: Tables S1–2). However, the inactivation curve was slightly but significantly shifted by about 4 mV in the depolarizing direction (Fig. 2d, Table 2), which suggests a small gain of channel function.

Structural modeling of the ID variant R937C predicts a destabilizing effect on the selectivity filter

The α-subunits of voltage-gated sodium channels contain four homologous repeats, which contribute to the formation of one central pore region (Fig. 3a-b). The SCN2A pore region contains the typical DEKA motif, which is formed by four conserved Asp/Glu/Lys/Ala residues (Fig. 3c). These residues represent the sodium selectivity filter and are therefore crucial for channel function (Shen et al. 2017). The position of one of these residues (E942) is stabilized by electrostatic interactions with the adjacent R937 (Fig. 3d). In the R937C variant these interactions cannot be formed by the shorter and uncharged cysteine sidechain resulting in a less defined position of the E942 sidechain. These structural changes offer an explanation for the electrophysiologically observed loss of conductance.

Fig. 3

Effect of amino acid exchanges on the SCN2A structure. a Top view on a model of the SCN2A structure. The four homologous repeats are shown in different colors (blue, cyan, orange, red) and the III-IV domain linker is shown in purple. The site of the pore is marked by a black diamond and the positions of the missense variants investigated in the present study are indicated. The extracellular sequence stretch spanning residues 275–359 has been omitted for clarity. b Side-view on the SCN2A structure. Color coding as in (a). c Model of the SCN2A pore region. The four residues that are critical for Na+ selectivity are shown in stick presentation. d In the wild-type SCN2A (left panel) the sidechain orientation of E942, which is part of the selectivity filter, is stabilized by polar interactions with R937 (see encircled region). These interactions cannot be formed in the R937C variant (right panel) by the shorter and uncharged cysteine sidechain thus impeding fixation of the E942 sidechain (flexibility is indicated by the magenta arrow). e E1803 forms interactions with residues Q1510 and P1512 of the III-IV domain linker, which cannot be formed in the (f) E1803G variant. The lacking interactions are highlighted by a blue dotted circle. g L1342 forms sidechain interactions with W864 of the adjacent repeat. h In the L1342P variant, these interactions cannot be formed by the less-extended proline sidechain. i L1563 forms hydrophobic interactions within the fourth domain. k These interactions are partially lost in the L1563V variant by the shorter valine sidechain (indicated by a blue dotted circle)

Residue E1803 is located in the globular C-terminal domain of SCN2A and forms sidechain interactions with two residues of the III-IV domain linker, which cannot be formed in the E1803G variant (Fig. 3e-f). The III-IV linker plays a key role for fast inactivation of voltage-gated sodium channels (Shen et al. 2017, Yan et al. 2017) suggesting that an altered interaction with adjacent domains, as observed for the E1803G variant, might affect inactivation properties.

The EE variant L1342P is predicted to disturb interactions with the adjacent S4-S5 linker

Residue L1342 is located in helix S4 and the effect of a L1342P exchange has been assessed previously based on the isolated helix (Hackenberg et al. 2014) due to the lack of a homologous 3D-structure at that time. This previous study revealed the lack of a backbone hydrogen bond in the mutant resulting in helix destabilization (Hackenberg et al. 2014). The availability of homologous structures in 2017 (Yan et al. 2017, Shen et al. 2017) allowed the generation of a complete 3D-model for SCN2A and a more comprehensive evaluation of the effects of the L1342P variant. An inspection of the SCN2A model reveals that L1342 forms sidechain interactions with W864, which is located in the S4-S5 linker of the adjacent repeat (Fig. 3g). These interactions cannot be formed in the L1342P variant because the cyclic proline sidechain cannot adopt an extended conformation required for W864 interaction (Fig. 3h). Based on this significant structural disturbance and the observation that a number of disease-related variants of Nav channels mapped to the interface between S4 and S5 segments in the neighboring repeats (Shen et al. 2017, Huang et al. 2017), L1342P is also likely to affect the biophysical properties of SCN2A.

Hydrophobic interactions within the fourth domain are expected to be diminished in the L1563V BFNIE variant

Residue L1563 is located in helix S2 of the fourth voltage-sensing domain and forms hydrophobic sidechain interactions in the core of the domain (Fig. 3i). In the L1563V variant some of these hydrophobic interactions are lost due to the shorter valine sidechain, but the overall packing of helices S1 and S2 remains unaffected (Fig. 3k). Thus, this amino acid change is expected to have only a moderate effect on the local structure of the fourth domain.

Trio WES in the published DEE patient carrying the R102* variant did not reveal a second hit

The p.R102* variant is to date the only functionally evaluated protein truncating SCN2A variant that leads to an encephalopathic phenotype with seizures (Kamiya et al. 2004). The reported patch-clamp studies showed a complete loss of channel function for this variant. This patient was reported to exhibit intractable seizures with an onset at age 1y 7m, ASD and severe ID. However, the mother was also reported to have suffered from recurrent febrile seizures until the age of 6y, suggesting the possibility of a second, maternally inherited variant in an epilepsy gene to cause the seizures in the patient and her mother while the de novo truncating SCN2A variant could explain the ASD/ID phenotype. This hypothesis is further supported by a recent report of a second patient carrying this variant (Monies et al. 2017) with a phenotype of severe ID and ASD but no seizures at the age of 6.5 years (personal communication Fowzan S. Alkuraya). To re-evaluate this case, we performed WES in the patient and her parents and filtered for dominant maternally inherited and recessive rare variants in known epilepsy and ID genes, as well as de novo variants. However, we did not detect a likely pathogenic sequence or copy number variant in any known epilepsy gene other than the already known SCN2A de novo variant.

Discussion

All three variants from patients with ID without seizures studied here showed a complete loss of Nav1.2-mediated currents. While the truncating variant L611Vfs*35 is likely to result in nonsense mediated mRNA decay (NMD) in vivo, the stop variant W1716* is located in the last exon and the mRNA should therefore not be subject to NMD (Sanders et al. 2018), but is expected to direct the synthesis of a C-terminally truncated protein. Thus, the complete loss of function in this mutant must be due to other mechanisms such as protein instability or disturbed transport. In line with this concept, Ben-Shalom et al. (2017) showed that in all three investigated truncating SCN2A variants no channel could be detected in the membrane of overexpressing HEK293 cells by immunocytochemistry. Of note, multiple variants truncating the protein at later positions (E1777, K1863, I1877) have also been reported to result in complete loss of conductance (Mantegazza et al. 2001).

It was less apparent how the missense variant R937C caused a complete loss of function, as the R937C mutant protein has been previously shown to be correctly located in the cell membrane (Ben-Shalom et al. 2017). Its critical location in the pore loop let us speculate that the pore could be completely blocked by the amino acid substitution. Our homology model supports this idea, predicting that R937 stabilizes the glutamate residue which is an essential part of the selectivity filter of Nav1.2. The same structural mechanism also applies to the previously described R937H variant causing ASD, which also exhibited no sodium flux (Ben-Shalom et al. 2017). Interestingly, two other variants for which no or reduced sodium conductance was reported (Ben-Shalom et al. 2017) are also located in the immediate vicinity of the DEKA selectivity filter (Fig. 3c). R379H is located near D384 and T1420M near K1422 of the selectivity filter. Thus, it is likely that these variants interfere with sodium permeability by a mechanism similar to that deduced above for the R937C variant.

To date, only one other study has functionally evaluated SCN2A variants with a non-epileptic phenotype (Additional file 1: Table S3) (Ben-Shalom et al. 2017). In their analysis of variants associated with ASD, all truncating variants and most missense variants were non-conducting, including R937C also investigated in our study. Notably, all non-conducting variants were detected in patients without additional seizures.

Both variants causing a phenotype of EE led to profound changes of Nav1.2 channel function. While L1342P affected several gating properties, E1803G mainly showed an increased window current caused by a shallower slope of the inactivation curve. So far, 10 SCN2A variants detected in EE patients have been studied electrophysiologically (Additional file 1: Table S3) (Kamiya et al. 2004, Ogiwara et al. 2009, Liao, Anttonen, et al. 2010, Lossin et al. 2012, Wolff et al. 2017, Berecki et al. 2018, Lauxmann et al. 2018). Two of them, E1211K and R1312T, exhibited gating changes similar to the ones observed in our mutant L1342P with a negative shift in both activation and inactivation curves and slowed recovery from inactivation (Ogiwara et al. 2009, Lossin et al. 2012). Interpretation based on these profoundly altered biophysical properties alone is challenging. However, a computational model of a cortical pyramidal neuron carrying the E1211K variant showed a marked neuronal hyperexcitability with a hyperpolarized spike threshold and an increased spike rate (Additional file 1: Table S4) (Ben-Shalom et al. 2017). We therefore also assume an increase in excitability for neurons carrying the L1342P variant. To date, three other patients carrying the de novo variant L1342P have been reported (Matalon et al. 2014, Dimassi et al. 2016, Wolff et al. 2017). All four affected patients share a common phenotype with an early seizure onset between age 3 to 6 months, severe ID, muscle tone abnormalities, (progressive) brain atrophy in MRI, and no eye contact (3/4) or cortical visual impairment (1/4). The seizures were intractable in 3/4 patients, and the EEG showed hypsarrhythmia in all three patients where EEG was reported and multifocal epileptic activity in 2/3. Acquired microcephaly and choreiform movements were noted in 2/4 patients.

Several missense variants located in the cytoplasmic C-terminal domain (CTD) leading to EE have been reported (Ben-Shalom et al. 2017). Interestingly, the only functionally studied variant, R1882Q, was found to have altered activation and inactivation curves leading to an increased window current, similar to our findings for E1803G. Our structural model predicts that the contacts between the CTD and the inactivation gate are disturbed (Fig. 3e-f), providing a mechanistic explanation for the altered inactivation properties observed in the E1803G mutant. Interestingly, gating changes leading to a persistent current seem to be a recurrent pathomechanism for SCN2A variants detected in EE patients, as this dysfunction has also been described previously as the mechanism of disease for the EE variants A263V (Liao, Anttonen, et al. 2010a) and V423L (Wolff et al. 2017). This pathomechanism was also observed in hippocampal pyramidal neurons of a Q54 transgenic mouse model carrying a gain-of-function Scn2a variant and showing a progressive epilepsy phenotype (Kearney et al. 2001). Persistent currents have been found to substantially influence action potential thresholds in axons and are considered to be important for generation of high-frequency action potential bursts in axons and nodes of Ranvier, the main expression sites of Nav1.2 (Kole and Stuart 2012).

Our study of six SCN2A variants representing the phenotypic variability known for this gene supports the association of their respective functional impact with the observed patient phenotype. Our data corroborates the assumption that 1) a loss of channel function caused by de novo truncating or missense variants is the underlying mechanism that leads to ID 2) benign epilepsy is due to usually inherited missense variants that lead to a small but relevant gain of channel function and 3) epileptic encephalopathies are the result of profound and lasting gating changes due to de novo missense variants (Ben-Shalom et al. 2017, Sanders et al. 2018, Berecki et al. 2018). Additionally, we suggest that a complete loss of Nav1.2 conductance is associated with ID/ASD without seizures, as supported by our own data and the results reported by Ben-Shalom et al. (2017), although like in our study long-term follow-up is required since seizures may appear later in childhood. This assumption is in line with data from a heterozygous knock-out mouse model which did not show a seizure phenotype (Planells-Cases et al. 2000, Middleton et al. 2018). However, another study on knock-out mouse models despite the absence of behavioral seizures showed very short absence seizures-like patterns in electrocorticography (Ogiwara et al. 2018), and multiple truncating SCN2A variants have been reported in patients with EE that pose exceptions to this proposition (Kamiya et al. 2004, Carvill et al. 2013, Horvath et al. 2016, Moller et al. 2016, Wolff et al. 2017). These truncating variants, though functionally studied in only one case (Kamiya et al. 2004), are expected to result in a complete loss of function and would therefore well explain the developmental delay and mental impairment but not the development of seizures. A closer look at the respective reports revealed that a patient reported with moderate ID and seizures between the age of 4y 7m and 9y by Wolff et al. (2017) and carrying the variant W1716* is the same patient reported in this study, who was suspected but never proven to have seizures of absence type that were later interpreted as rather being mutistic episodes. The patient carrying the de novo variant R102* reported with intractable epilepsy, severe ID and ASD had a mother with recurrent febrile seizures (Kamiya et al. 2004) suggesting an inherited second hit contributing to the seizure phenotype, albeit not detected by our WES analysis. Accordingly, the same variant in an unrelated patient (Monies et al. 2017) with follow-up until age 6.5 years did not cause any seizures (personal communication Fowzan S. Alkuraya). The variant R1435* was reported twice as well, once in a patient with ASD, ID and seizures (Trump et al. 2016) and once in a patient with ID and ASD without seizures at the last follow-up at age 7y 8m (Wolff et al. 2017), also suggesting a possible second hit or other modifier in the first patient. The assumption of a modifying effect of the genetic background is supported by two studies where different strains of mice carrying the same Scn2a variant exhibited differences in seizure severity (Kearney et al. 2001, Thompson et al. 2017). In addition, acquired factors such as hypoxia may also impact Nav1.2 function (Plant et al. 2016, Winquist and Cohen 2018).

It has recently been shown that voltage-gated sodium channels assemble and gate as dimers (Clatot et al. 2017) and that some pathogenic Nav1.5 variants exert a dominant-negative effect through impairing trafficking and gating of the wild-type channel (Clatot et al. 2012, Clatot et al. 2018). A possible dominant-negative effect was also discussed for pathogenic Nav1.2 variants since Kamiya et al. (2004) described a reduced whole-cell sodium current when co-expressing the truncating variant R102* together with wild-type. However, this effect disappeared when co-expressing β-subunits. Ben-Shalom et al. (2017) did not find a dominant-negative effect either when testing a non-conducting missense variant. Currently it can only be speculated if some of the variants might exert a dominant-negative effect and thereby could explain some outliers of the assumed genotype-phenotype correlation. Hence, this potential pathomechanism should be kept in mind and addressed in future studies using suitable models resembling the heterozygous constellation in the patient.

Another point to consider regarding the pathomechanism of SCN2A variants is the isoform-dependent effect of Nav1.2 variants. Next to the ‘adult’ Nav1.2 isoform, which was used by all previous electrophysiological investigations utilizing overexpression of human SCN2A as well as in our study, a ‘neonatal’ isoform exists differing only in a single amino acid at position 209 of the protein (Kasai et al. 2001, Gazina et al. 2015). The ‘neonatal’ isoform is expressed at different levels at birth and then decreases until levelling at 10–20% of the ‘adult’ isoform at 20 weeks postnatally in mice (Gazina et al. 2010). Electrophysiological studies indicate that pathogenic variants can exert different effects in this ‘neonatal’ isoform compared to the predominantly studied ‘adult’ isoform (Xu et al. 2007, Liao, Deprez, et al. 2010b).

Conclusions

In conclusion, our work corroborates a strong genotype-phenotype correlation of deleterious SCN2A variants and demonstrates that structural modeling can be a useful tool to predict the severity of the variant impact, suggesting its potential as a useful and fast tool to support variant interpretation in diagnostics. We further illustrate the mechanism of how missense variants located in the pore-forming loop can lead to a loss of channel function and how C-terminally located missense variants can exert a severe disturbance on Nav1.2 gating properties.

Funding

Supported by radiz—Rare Disease Initiative Zürich, Clinical Research Priority Program for Rare Diseases of the University of Zurich to AR, BP and HUZ, by the Swiss National Science Foundation (SNSF) grant 320030_179547 to AR, and by the Forschungskredit UZH and the Josef Huwyler Ruth Bernet-Engeli Stiftung to ABe.

Availability of data and materials

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Search for Reza Asadollahi in:

Search for Heinrich Sticht in:

Search for Hanns Ulrich Zeilhofer in:

Search for Anita Rauch in:

Contributions

ABe and MA performed the experiments, analyzed and interpreted data and drafted the manuscript. MZ interpreted the data and revised the manuscript. MV, KS, RBG, AH, LA, BP, JKS, ABa, KY, and YI contributed clinical and/or genetic data. RA contributed to the design of the study. HS performed and analyzed the structural modeling. HUZ and AR designed the study, interpreted the data and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Ethics declarations

Ethics approval and consent to participate

This study was approved by the ethical committee of the Kanton of Zurich (reference StV 11/09) and informed consent was obtained from the patient or their legal guardians.

Consent for publication

Consent for publication of molecular and clinical data has been obtained from all reported patients or their legal guardians.

Competing interests

BP has received travel reimbursement from Nutricia in 2018 and honorarium for a lecture on vitamin B6 dependent seizures by Alexion in 2017 (Bad Nauheim). The remaining authors have no conflicts of interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.