Next Article in Journal
Mitochondrial Oxidative Stress Induces Cardiac Fibrosis in Obese Rats through Modulation of Transthyretin
Next Article in Special Issue
Abnormal Expression of Synaptic and Extrasynaptic GABAA Receptor Subunits in the Dystrophin-Deficient mdx Mouse
Previous Article in Journal
Analysis of the Biological Properties of Blood Plasma Protein with GcMAF Functional Activity
Previous Article in Special Issue
Translocation of TMEM175 Lysosomal Potassium Channel to the Plasma Membrane by Dynasore Compounds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia

1
Department of Pharmacy-Drug Sciences, University of Bari “Aldo Moro”, 70125 Bari, Italy
2
Department of Physics, Université de Montréal, Montreal, QC H3C 3J7, Canada
3
CIRCA, Center for Interdisciplinary Research on Brain and Learning, Université de Montréal, Montreal, QC H3C 3J7, Canada
4
Department of Genetics, Hadassah Medical Center, Jerusalem 91120, Israel
5
Department of Medicine and Surgery, LUM University, 70010 Casamassima, Italy
6
Pediatric Neurology Clinic, Shamir Medical Center (Assaf Harofeh), Zerifin 7033001, Israel
7
Department of Pharmacology and Physiology, Université de Montréal, Montreal, QC H3C 3J7, Canada
8
Child Neurology Unit, IRCCS, Institute of Neurological Sciences of Bologna, 40139 Bologna, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(15), 8079; https://doi.org/10.3390/ijms23158079
Submission received: 13 June 2022 / Revised: 15 July 2022 / Accepted: 18 July 2022 / Published: 22 July 2022
(This article belongs to the Special Issue Implication of Ion Channels in Neurodevelopmental Disorders)

Abstract

:
Mutations in the KCNA1 gene, encoding the voltage-gated potassium channel Kv1.1, have been associated with a spectrum of neurological phenotypes, including episodic ataxia type 1 and developmental and epileptic encephalopathy. We have recently identified a de novo variant in KCNA1 in the highly conserved Pro-Val-Pro motif within the pore of the Kv1.1 channel in a girl affected by early onset epilepsy, ataxia and developmental delay. Other mutations causing severe epilepsy are located in Kv1.1 pore domain. The patient was initially treated with a combination of antiepileptic drugs with limited benefit. Finally, seizures and ataxia control were achieved with lacosamide and acetazolamide. The aim of this study was to functionally characterize Kv1.1 mutant channel to provide a genotype–phenotype correlation and discuss therapeutic options for KCNA1-related epilepsy. To this aim, we transfected HEK 293 cells with Kv1.1 or P403A cDNAs and recorded potassium currents through whole-cell patch-clamp. P403A channels showed smaller potassium currents, voltage-dependent activation shifted by +30 mV towards positive potentials and slower kinetics of activation compared with Kv1.1 wild-type. Heteromeric Kv1.1+P403A channels, resembling the condition of the heterozygous patient, confirmed a loss-of-function biophysical phenotype. Overall, the functional characterization of P403A channels correlates with the clinical symptoms of the patient and supports the observation that mutations associated with severe epileptic phenotype cluster in a highly conserved stretch of residues in Kv1.1 pore domain. This study also strengthens the beneficial effect of acetazolamide and sodium channel blockers in KCNA1 channelopathies.

1. Introduction

The voltage-gated potassium channel Kv1.1, encoded by the KCNA1 gene, is widely distributed in the central and peripheral nervous system, where it regulates nerve cell repolarization after an action potential and contributes to neuronal excitability, firing rate and neurotransmitter release [1,2]. The channel, belonging to the Shaker-like potassium channel family, is formed by the tetrameric assembly of four Kv1.1 alpha subunits, each comprising six transmembrane segments that include functionally critical voltage-sensing (S1–S4) and pore domains (S5–S6) [3,4]. The KCNA1 gene is well-known to scientists because several mutations in this gene have been associated with episodic ataxia type 1 (EA1), a rare disorder of cerebellar origin characterized by recurrent episodes of ataxia and myokymia since early childhood [5]. The improvements in clinical and genetic diagnosis have contributed to broaden the phenotypic spectrum associated with KCNA1 mutations, which now includes seizures, cognitive impairment, migraine, paroxysmal kinesigenic dyskinesia, hypomagnesemia, metabolic alterations and cataplexy [6,7]. Interestingly, EA1 patients carrying KCNA1 mutations show increased incidence of epilepsy and, in the last years, variants in KCNA1 have been associated with epilepsy phenotypes with or without ataxia [8,9,10,11,12,13] (Table S1). This is not surprising if one considers that kcna1 knock out mice have increased seizures susceptibility and undergo premature death (SUDEP) [14]. Since 2018, at least four mutations in the pore domain of Kv1.1 channel have been described in patients with early onset developmental and epileptic encephalopathy (DEE) [15,16,17,18]. Three of these mutations alter the two proline residues in the Pro-Val-Pro (PVP) motif at the C-terminal end of S6 (P403S, P405S, P405L), whereas the fourth is located close to the selectivity filter for potassium ions in the S5–S6 loop (V368L). The PVP region is highly conserved amongst the voltage-gated K+ channels and has been shown to be crucial for S6 flexibility and channel gating [19,20]. At the corresponding positions in KCNA2, the de novo P405L and P407A mutations have been associated with neurodevelopmental syndromes in which epilepsy, motor and cognitive delay co-exist [21].
Here, we report another de novo missense variant, P403A, in the PVP motif of the Kv1.1 channel in a 5-year-old girl who presents epilepsy, developmental delay and ataxia. Through patch-clamp electrophysiology and molecular dynamics simulations, we elucidated the impact of the amino acid change on channel gating and carried out a review of the literature to gain insight into the genotype–phenotype drug response correlation of nearby mutations. Our study supports previous findings, indicating that the functional impairment of residues within the pore domain result in more severe phenotypes in which epilepsy coexist with developmental delay. We also suggest the need for a careful revision of the clinical phenotype and of the therapeutic management of patients carrying KCNA1 mutations to personalize the treatment and improve the therapeutic outcome.

2. Results

2.1. Clinical Diagnosis and Follow Up

The patient is a 5-year-2-month-old female who was born full term via spontaneous vaginal delivery after an uncomplicated pregnancy. Since previous pregnancy was terminated due to fetal tracheo-esophageal fistula, a chromosomal microarray analysis obtained by amniocentesis was performed, resulting negative.
At 5 days of life, she presented with a sudden onset of many seizures per week, characterized by staring. EEG showed frequent and independent medium–high amplitude sharp waves over bilateral frontotemporal regions, activated during NREM sleep (Figure 1).
Brain MRI, lumbosacral and renal ultrasound, echocardiogram, ophthalmologic examination and metabolic screening (blood lactate, amino acids, ammonia, urine organic acids, cerebrospinal fluid analysis for amino acids, glucose and lactate, very long chain fatty acid, neurotransmitters, blood and cerebrospinal fluid cultures) were normal. No dysmorphic features were present. Neurodevelopmental delay was noticed in the months following the seizure onset. At 13–15 months of age she presented other neurological signs, including eye blinking, ataxia and intentional tremor. Seizures were never controlled by antiepileptic treatment, despite attempts with several antiepileptic drugs (AEDs), such as levetiracetam, sulthiame, topiramate, valproate, cannabidiol, oxcarbazepine and clobazam. Moreover, a ketogenic diet was tried without benefits. Follow-up EEGs showed normal background activity with independent medium–high amplitude sharp wave over bilateral centrotemporal regions, with variable activation of these abnormalities during NREM sleep, but without configuring a continuous spike and wave pattern.
Seizures freedom and improvement in motor function were obtained when acetazolamide was added to lacosamide.
At last follow-up (5 years of age), she was able to walk using a hard walker and presented with intellectual disability, with an intentional distal tremor at the upper limbs and ataxia.

2.2. Genetic Diagnosis

Trio exome analysis of the index patient yielded 43 million mapped reads with a mean coverage of 60×. Following alignment to the reference genome (hg19) and variant calling, we performed a series of filtering steps. These included removing variants that were called less than ×8, were off-target, synonymous or had MAF  >  0.1% at gnomAD v2.1.1. There were no homozygous variants or pairs of variants in the same gene, within the exons of CGD genes among them. Compared with the parents’ samples, a heterozygous non-inherited variant was found: Chr12: 5021751 C > G, NM_000217.2: c.1207C > G, P403A, in the KCNA1 gene. This variant does not appear in large exome databases (gnomAD), the position is preserved, and the change is predicted to be pathogenic. In addition, a variant in a homologous position of the KCNA2 gene (P405L) has been previously described as pathogenic [22].

2.3. Functional and Molecular Characterization of Kv1.1P403A Channels

The variant P403A is located in the PVP motif in the middle of the S6 segment of the Kv1.1 channel (Figure 2A). This residue is highly conserved within the members of the Kv1 family (Figure 2B).
As the patient is heterozygous for the mutation, Kv1.1 channels will likely be composed of Kv1.1 wild-type (WT) and P403A mutant subunits. To test the hypothesis that the P403A variant altered Kv1.1 function and caused epilepsy and ataxia in the affected patient, we transfected HEK 293 cells with equal amount of WT (5 μg) or P403A (5 μg) cDNAs alone or in 1:1 ratio (5 μg + 5 μg). As shown in Figure 3A–C, transfected HEK cells displayed typical fast activating non-inactivating potassium currents. The current amplitude and biophysical properties of potassium currents elicited by P403A and Kv1.1 WT+P403A channels were then compared with those obtained from WT currents. As shown in Figure 3D, P403A channels generate significantly smaller potassium currents at every tested potential compared with Kv1.1 WT channels (0.4 ± 0.04 vs. 2.1 ± 0.2 at + 20 mV for P403A and WT, respectively). The co-expression of Kv1.1 WT and P403A subunits gave rise to potassium currents that were less than half of the currents expected for WT channels alone (0.55 ± 0.08 at + 20 mV), suggesting a dominant-negative effect exerted by the mutant subunit on the Kv1.1 WT (Figure 3D; Table 1).
To determine whether the P403A variant may affect the voltage dependent activation, tail currents were recorded at −30 mV (for P403A and Kv1.1 WT+P403A) or −50 mV (for Kv 1.1 WT) following prepulse commands to several voltages (Figure 3E), and data points were fitted with a Boltzmann function. Mutant P403A channels displayed voltage-dependent activation significantly shifted by ~ 30 mV toward positive potentials compared to Kv1.1 WT, an effect that predicts loss-of-function (LoF) (Table 1). Potassium currents resulting from the co-transfection of WT and P403A showed voltage-dependent gating that falls between that of Kv1.1 WT and P403A homomeric channels, with a ~17 mV depolarizing shift of V1/2 compared to Kv1.1 WT. Conversely, the slope factor, k, calculated from the Boltzmann fit of tail currents was almost unaffected by the mutation (Figure 3E, Table 1).
To investigate if the P403A variant affected the kinetics of activation and deactivation of Kv1.1 channels, the activating and deactivating current traces of either Kv1.1 WT, P403A or Kv1.1 WT+P403A channels were fitted with a single-exponential function and the calculated time constants were plotted as a function of membrane potential. The kinetics of activation of P403A channels were slower than those of Kv1.1 WT at each tested potential, indicating a LoF effect. P403A channels had kinetics of deactivation comparable with those of Kv1.1 WT (Figure 4A,B, Table 1). No differences were observed in the C-type inactivation of the three channel types (Figure 4C, Table 1). However, P403A and Kv1.1 WT+P403A recovered slightly faster from the C-type inactivation with respect to Kv1.1 WT channels (Figure 4D, Table 1).
To understand the structural basis and molecular mechanism of the P403A mutation in Kv1.1, we created a Kv1.1 homology model based on the Kv1.2/2.1 crystal structure [4,6]. We introduced the P403A mutation in all four subunits (Figure 5A; see Section 4 for details). The sixth transmembrane helix (S6), which forms the central ion conducting pore together with the pore loop and the S5, has a pronounced kink in the WT channel caused by the PVP motif. P403A changes this motif to AVP, and, accordingly, the kink is less pronounced. The angle of the S6 kink and the angle of the C-terminal part with respect to the membrane were altered (Figure 5B). However, this effect was not consistent in all four subunits but more pronounced in one subunit versus the others (Figure 5C). In accordance with this differential conformational change, the extracellular ends of the S6 are no longer symmetrical (Figure 5B,D). Simulations of the WT channel did not show this asymmetry. We should note that in order to preserve the integrity of the selectivity filter, it had to be stabilized during the simulation. This need for stabilization and the induced asymmetry suggest that the open state of the P403A mutant is not a stable conformation. A destabilized open state would push the channel towards the closed state in accordance with the LoF phenotypes described above. A rotation of the S6 C-terminal to the PVP motif has been suggested to close the ion conduction pathway in Kv channels by turning a hydrophobic residue into the permeation pathway and dewetting the pore [23,24].
Since the simulations start with the folded structure of the WT channel and the open state structure of the P403A mutant seems to be strained, it remains a possibility that the P403A channel directly folds differently than the wildtype channel.

3. Discussion

3.1. The PVP Motif of Kv1.1 Channels as a Hotspot for Mutations Associated with Severe Epilepsy

As mentioned before, genetic variants in both KCNA1 and KCNA2, which disrupt the proline residues of the PVP motif, have been identified in children with severe epilepsy and comorbidities (P403S and P405S/L in Kv1.1 and P405L and P407A in Kv1.2; Table S1). In this study, we have characterized an additional mutation in the PVP motif of the Kv1.1 channel, the P403A, associated with an epileptic phenotype.
The contribution of the PVP motif to Kv1 channel gating and the role of proline-induced kinks in transmembrane helices are well known. As suggested by the 3D structure of the rat Kv1.2 channel and of the Kv1.2/2.1 (paddle) chimera [3,4] and shown in previous mutagenesis studies, pore opening is thought to involve bending of the inner S6 helices in the vicinity of glycine gating hinge and around the PVP motif [20,25,26]. The PVP motif thus acts as a helix-breaking sequence, which provides a flexible hinge that allows the lower half of the S6 to move [23,24]. Thus, replacing proline residues in the PVP region affects S6 flexibility and pore opening and leads to drastically reduced channel activity (LoF defect). Indeed, Kv1.1 P403S and P405L channels failed to generate measurable currents when expressed in CHO cells, whereas the P403A and P405S channels showed a marked depolarizing shift in voltage-dependent activation and a significant reduction in current density [18]. In addition, the P403A subunit exerts a clear dominant-negative effect on the wild-type Kv1.1 subunit, further lowering potassium flow. Similar results have been reported for the corresponding substitution P405L in Kv1.2, which has been found in several patients with severe drug resistant phenotypes. The P405L mutant shows a dramatic reduction in current density and exerts a dominant-negative effect on Kv1.2 WT channels [22]. Of note, a de novo copy number variation in KCNA1 extending from the PVP motif to the end of the S6 helix has been identified in a 3-year-old proband with serious myoclonic epilepsy in infancy who died of SUDEP [27].
As expected, our modelling study shows a destabilization of the open state and suggests a defective folding of the P403A channel that agrees with the reduced current amplitude and with the dominant-negative effect of the mutant subunit on the WT. Similarly, a modelling study performed on Kv1.2 channel proposed that the PVP motif may also interact with the S4–S5 linker. The substitutions of P405 and P407 in Kv1.2 by disrupting the symmetry of the PVP, might impair interactions with S4–S5 during gating-related conformational changes, leading to altered channel gating [21]. Additional evidence supports the critical role of the PVP motif for potassium channel activity and suggests that this motif can be a hotspot of mutations underlying severe epileptic phenotypes. The introduction of the P407A mutation in the Shaker Kv channel generates a non-conducting mutant channel [19]. In human Kv1.5, replacing the first proline of the PVP by an alanine results in a nonfunctional channel [28]. The mutation P403A in the PXP motif of Kv4.2 channel has been associated with early onset global developmental delay [29].
Trying to draw a correlation between epilepsy phenotype and mutation localization in KCNA1 diseases, it appears that KCNA1 pore mutations affecting the selectivity filter or the proline residues in the PVP motif are associated with epileptic phenotypes aggravated by developmental delay and sometimes drug-resistance, whereas aminoacidic substitutions outside the pore and the PVP domain (such as A261T, S3) seem to occur in patients with milder epileptic phenotypes, usually with a favorable seizure outcome and without intellectual disability [7,30,31] (Table S1).

3.2. Genotype–Phenotype–Drug Response Correlation

Regardless of the position on channel structure and clinical phenotype, all the KCNA1 mutations identified to date, except for the A261T, cause LoF of the Kv1.1 channel expressed in heterologous cells [1,18]. Therefore, a genotype–phenotype correlation is not simple for KCNA1 mutations and additional unclear factors, besides Kv1.1 impairment and perhaps secondary to Kv1.1 LoF, may contribute to explain different clinical phenotypes, such as those seen in twins (P403S) or among carriers of the same variant [6,21]. So far, epilepsy seems less likely to occur because of the KCNA1 mutation with respect to EA1, requiring perhaps more deleterious defects of Kv1.1 activity, and less frequent de novo modifications and recessive inheritance have been reported. The ratio of WT/mutant subunits in the tetramer, the presence of interacting proteins, the interplay of rare and common variant in the patient, the occurrence of neuron-selective disruption and the role of Kv1.1 in neurogenesis are among the factors that could impact Kv1.1 membrane expression and gating and hence contribute to neuronal network impairment and clinical heterogeneity [1,32,33].
The presentation of seizures alongside other neurological and mental health issues in KCNA1 carriers means that accurate diagnosis and AEDs prescribing is challenging. As in other rare neurologic channelopathies, standard therapy has not yet been established and different treatment attempts need to be evaluated case by case to reach a personalized therapeutic scheme able to reduce seizure frequency and improve patient’s quality of life. Acetazolamide, a carbonic anhydrase inhibitor, is often used to treat ataxia, with moderate benefit in most patients, including the girl described here. Sodium channel blockers proved effective to restore brain excitability in both EA1 and epilepsy [18,34,35] (Table S1). In severe epilepsy cases, however, seizures may remain intractable or poorly improve after trials with combinations of AEDs, including phenobarbital, clobazam, valproate, lamotrigine, carbamazepine and phenytoin. The P405S carrier experienced severe intractable seizures [16], whereas the patient bearing P405L slightly improved only with a combination of acetazolamide, lamotrigine, valproic acid and ACTH to specifically treat status epilepticus during sleep [17]. The P403S mutation was identified in twins showing epileptic phenotypes of different severity, with one boy presenting convulsive episodes controlled by lamotrigine and the other drug-resistant seizures [16]. The recessively inherited V368L mutation was found in a patient presenting with a combination of dyskinesia, motor and intellectual disability and neonatal epileptic encephalopathy, controlled with oxcarbazepine [15]. After having tried combinations of AEDs, seizures and ataxia were finally controlled in our patient (P403A mutation) with acetazolamide and lacosamide. Lacosamide, a drug used to treat focal epilepsy as monotherapy or as an add-on therapy, differs from other sodium channel blockers as it enhances the slow inactivation of the sodium channel without affecting the fast inactivation and also acts as a carbonic anhydrase inhibitor [36], a dual mechanism that may well motivate its efficacy. Cannabis-based medical products have recently been approved to treat seizures associated with intractable Lennox–Gastaut syndrome, Dravet syndrome or tuberous sclerosis complex in patients of 1 year of age and over [37]. Cannabidiol was unsuccessful in our patient and was rapidly discontinued, as well as ketogenic diet. Overall, the comparison of the treatment regimens of patients carrying a KCNA1 mutation confirms the need of a patient-centric approach and supports the use of sodium channel blockers, such as carbamazepine, lacosamide, lamotrigine, to improve the course of the disease [34,38]. Caution must be taken when using sodium valproate in children due to possible seizure aggravation and clonazepam and oxcarbazepine due to motor coordination disturbances [39].
In conclusion, in depth structure–function and pharmacological studies are essential to understand the consequences of pathogenic mutations on channel activity and clinical phenotype and can also be fundamental for the rational design of drugs targeting specific mutations defects. An improved knowledge of patients’ genetic background, including the presence of common variants, protein–protein interactions and brain networks functioning will be necessary to better understand the molecular mechanisms underlying KCNA1 channelopathies and to support therapeutic decisions.

4. Materials and Methods

4.1. Clinical and Genetic Analysis

Following informed consent, exonic sequences were enriched in the DNA sample of the proband and the parents using SureSelect Human All Exon v.5 50 Mb Kit (Agilent Technologies, Santa Clara, CA, USA). Sequences were determined by HiSeq2500 (Illumina, San Diego, CA, USA) as 100 bp paired-end runs. Data analysis, including read alignment and variant calling, was performed by DNAnexus software (Palo Alto, CA, USA) using the default parameters with the human genome assembly hg19 (GRCh37) as reference.

4.2. Electrophysiology

Whole-cell patch-clamp experiments were performed at room temperature (~20 °C) using an Axopatch 200B amplifier and Digidata 1550B AD/DA converter (Axon Instruments, Sunnyvale, CA, USA). Currents were low-pass filtered at 2 kHz and digitized with sampling rates of 50 kHz. The bath solution contained (mM) NaCl 142, KCl 2.8, MgCl2 1, CaCl2 1, HEPES 10, glucose 11, pH = 7.4, whereas the pipettes were filled with solution containing (mM) NaCl 10, K-glutamate 132, MgCl2 2, CaCl2 0.9, EGTA 1, HEPES 10, pH = 7.4 [40]. The pipettes were pulled from borosilicate glass and had ~4 MΩ resistance.
Potassium currents were evoked by 200 ms depolarizing commands from a holding potential of −80 to +60 mV (10 mV intervals). To measure tail currents, each prepulse was followed by a 150 ms voltage step at −50 or −30 mV for Kv1.1 wild-type and P403A channels, respectively. Normalized tail currents were plotted as a function of the prepulse potential, and data points fitted with the Boltzmann function I = 1/1 + exp{−(V − V1/2)/k} to determine the voltage dependence of channel activation. V1/2, the half-maximal activation potential, and k, the slope factor, were calculated from the fit.
To measure the kinetics of activation, currents were elicited by 200 ms depolarizing pulses from a holding potential of −80 to +60 mV (10 mV intervals). To measure the kinetics of deactivation, currents were elicited by 200 ms depolarizing pulse at +20 mV, followed by 200 ms depolarizing pulses from −80 to +20 mV (10 mV intervals). The kinetics of activation and deactivation were determined by fitting activating and deactivating current traces with a single exponential function. The calculated time constants were plotted as a function of voltage and fitted with the equation τ = τV1/2 exp(V−V1/2)/k, where τV1/2 is the time constant at the V1/2 of the channels, and k is the slope factor for the voltage dependence of the time constants.
To determine the slow C-type inactivation, a test pulse to +40 mV was delivered for 90 s to cells expressing wild-type and mutant channels. The kinetics of inactivation were estimated by fitting the time course of current decay with a double exponential function and calculating the fast (τfast) and slow (τslow) time constants and the relevant amplitudes (A%).
The recovery from the C-type inactivation was determined by using a double-pulse protocol to +20 mV, separated by interpulse intervals of increasing duration (from 0.1 to 17.1 s). The peak currents evoked by the second pulse (100 ms) were divided by the first pulse (20 s), normalized and plotted as a function of the interpulse duration. Data points were fitted with an exponential function, from which the time constants were calculated.
Data were analyzed by using pCLAMP 10.3 (Axon Instruments, Sunnyvale, CA, USA) and GraphPad Software. Results were reported as mean ± SEM from n cells, and statistical analysis was performed using Student’s t-test, with p < 0.05 (*) considered as significant.

4.3. Homology Modeling and Molecular Dynamics Simulations

A homology model of Kv1.1 was built from the crystal structure of the Kv1.2/2.1 chimera (PDB 2R9R) [4], as described previously [6]. Residues that differed or were missing in the Kv1.2/2.1 chimera as compared with the human Kv1.1 were replaced and modeled, respectively, using Modeller 9.25 [41,42]. We introduced the mutation P403A in the four subunits in the Kv1.1 tetramer. The N- (1–147) and C-termini (420–495) were truncated. The resulting homotetrameric transmembrane domain was introduced into a POPE: POPC: POPS membrane (3:2:1; 444 lipids total) and inserted in a hexagonal box (diagonal 170 Å, height 110 Å; filled with 150 mM KCl) using CHARMM-GUI [43,44] and was equilibrated using NAMD 2.14 with the CHARMM36 force field [45]. The simulation was run for 280 ns after equilibration.
As control condition, a homotetramer of the Kv1.1 WT channel was modelled and simulated identically. For some simulations, the selectivity filter was stabilized by imposing the geometry of one carbonyl group and the position of a single potassium ion in the selectivity filter. Analysis and graphical representation were created using in-house written python and Matlab routines using the open source Pymol library.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23158079/s1.

Author Contributions

Genetic analysis, H.M.S. and O.E.; clinical investigation, I.H. and A.R.; molecular biology, E.C. and G.D.; patch clamp experiments, M.C.D. and M.A.C.; molecular modeling, R.B. and M.M.; writing—original draft preparation, P.I., R.B. and A.R.; writing—review and editing, P.I., A.L. and A.D.L.; funding, R.B. and P.I. All authors have read and agreed to the published version of the manuscript.

Funding

We gratefully acknowledge the financial support by Fondo Ordinario Ricerca Scientifica Uniba 2017-2018 to P.I.; Canadian Institutes for Health Research (CIHR) Grant PJT-169160 and the Natural Sciences and Engineering Research Council (NSERC) Grant RGPIN-2017-06871 to R.B.; This research was enabled in part by support provided by Calcul Québec (www.calculquebec.ca) and Compute Canada (www.computecanada.ca).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Hadassah Medical Center (approval code 0306-10-HMO, 8 June 2021).

Informed Consent Statement

Informed consent for publication was obtained from the parents of the patient involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Acknowledgments

We thank the patient and her family for participation in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. D’Adamo, M.C.; Liantonio, A.; Rolland, J.-F.; Pessia, M.; Imbrici, P. Kv1.1 Channelopathies: Pathophysiological Mechanisms and Therapeutic Approaches. Int. J. Mol. Sci. 2020, 21, 2935. [Google Scholar] [CrossRef] [PubMed]
  2. Paulhus, K.; Ammerman, L.; Glasscock, E. Clinical Spectrum of KCNA1 Mutations: New Insights into Episodic Ataxia and Epilepsy Comorbidity. Int. J. Mol. Sci. 2020, 21, 2802. [Google Scholar] [CrossRef] [PubMed]
  3. Long, S.B.; Campbell, E.B.; Mackinnon, R. Crystal structure of a mammalian voltage-dependent Shaker family K+ channel. Science 2005, 309, 897–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Long, S.B.; Tao, X.; Campbell, E.B.; MacKinnon, R. Atomic structure of a voltage-dependent K+ channel in a lipid membrane like environment. Nature 2007, 450, 376–382. [Google Scholar] [CrossRef]
  5. Browne, D.L.; Gancher, S.T.; Nutt, J.G.; Brunt, E.R.P.; Smith, E.A.; Kramer, P.; Litt, M. Episodic ataxia/myokymia syndrome is associated with point mutations in the human potassium channel gene, KCNA1. Nat. Genet. 1994, 8, 136–140. [Google Scholar] [CrossRef]
  6. Imbrici, P.; Accogli, A.; Blunck, R.; Altamura, C.; Iacomino, M.; D’adamo, M.C.; Allegri, A.; Pedemonte, M.; Brolatti, N.; Vari, S.; et al. Musculoskeletal Features without Ataxia Associated with a Novel de novo Mutation in KCNA1 Impairing the Voltage Sensitivity of Kv1.1 Channel. Biomedicines 2021, 9, 75. [Google Scholar] [CrossRef]
  7. Zuberi, S.M.; Eunson, L.H.; Spauschus, A.; De Silva, R.; Tolmie, J.; Wood, N.W.; McWilliam, R.C.; Stephenson, J.B.; Stephenson, J.P.; Kullmann, D.M.; et al. A novel mutation in the human voltage-gated potassium channel gene (Kv1.1) associates with episodic ataxia type 1 and sometimes with partial epilepsy. Brain 1999, 122, 817–825. [Google Scholar] [CrossRef] [Green Version]
  8. Adelman, J.P.; Bond, C.T.; Pessia, M.; Maylie, J. Episodic ataxia results from voltage-dependent potassium channels with altered functions. Neuron 1995, 15, 1449–1454. [Google Scholar] [CrossRef] [Green Version]
  9. Yin, X.M.; Lin, J.H.; Cao, L.; Zhang, T.M.; Zeng, S.; Zhang, K.L.; Tian, W.T.; Hu, Z.M.; Li, N.; Wang, J.L.; et al. Familial paroxysmal kinesigenic dyskinesia is associated with mutations in the KCNA1 gene. Hum. Mol. Genet. 2018, 27, 757–758. [Google Scholar] [CrossRef] [Green Version]
  10. Tristán-Clavijo, E.; Scholl, F.G.; Macaya, A.; Iglesias, G.; Rojas, A.M.; Lucas, M.; Castellano, A.; Martinez-Mir, A. Dominant-negative mutation p.Arg324Thr in KCNA1 impairs Kv1.1 channel function in episodic ataxia. Mov. Disord. 2016, 31, 1743–1748. [Google Scholar] [CrossRef]
  11. Lee, H.; Wang, H.; Jen, J.C.; Sabatti, C.; Baloh, R.W.; Nelson, S.F. A novel mutation in KCNA1 causes episodic ataxia without myokymia. Hum. Mutat. 2004, 24, 536. [Google Scholar] [CrossRef]
  12. Yuan, H.; Yuan, H.; Wang, Q.; Ye, W.; Yao, R.; Xu, W.; Liu, Y. Two novel KCNA1 variants identified in two unrelated Chinese families affected by episodic ataxia type 1 and neurodevelopmental disorders. Mol. Genet. Genom. Med. 2020, 8, e1434. [Google Scholar] [CrossRef]
  13. Demos, M.K.; Macri, V.; Farrell, K.; Nelson, T.N.; Chapman, K.; Accili, E.; Armstrong, L. A novel KCNA1 mutation associated with global delay and persistent cerebellar dysfunction. Mov. Disord. 2009, 24, 778–782. [Google Scholar] [CrossRef]
  14. Smart, S.L.; Lopantsev, V.; Zhang, C.L.; Robbins, C.A.; Wang, H.; Chiu, S.Y.; Schwartzkroin, P.A.; Messing, A.; Tempel, B.L. Deletion of the Kv1.1 potassium channel causes epilepsy in mice. Neuron 1998, 20, 809–819. [Google Scholar] [CrossRef] [Green Version]
  15. Verdura, E.; Fons, C.; Schlüter, A.; Ruiz, M.; Fourcade, S.; Casasnovas, C.; Castellano, A.; Pujol, A. Complete loss of KCNA1 activity causes neonatal epileptic encephalopathy and dyskinesia. J. Med. Genet. 2020, 57, 132–137. [Google Scholar] [CrossRef]
  16. Rogers, A.; Golumbek, P.; Cellini, E.; Doccini, V.; Guerrini, R.; Wallgren-Pettersson, C.; Thuresson, A.C.; Gurnett, C.A. De novo KCNA1 variants in the PVP motif cause infantile epileptic encephalopathy and cognitive impairment similar to recurrent KCNA2 variants. Am. J. Med. Genet. A 2018, 176, 1748–1752. [Google Scholar] [CrossRef]
  17. Russo, A.; Gobbi, G.; Pini, A.; Møller, R.S.; Rubboli, G. Encephalopathy related to status epilepticus during sleep due to a de novo KCNA1 variant in the Kv-specific Pro-Val-Pro motif: Phenotypic description and remarkable electroclinical response to ACTH. Epileptic Disord. 2020, 22, 802–806. [Google Scholar] [CrossRef]
  18. Miceli, F.; Guerrini, R.; Nappi, M.; Soldovieri, M.V.; Cellini, E.; Gurnett, C.A.; Parmeggiani, L.; Mei, D.; Taglialatela, M. Distinct epilepsy phenotypes and response to drugs in KCNA1 gain- and loss-of function variants. Epilepsia 2022, 63, e7–e14. [Google Scholar] [CrossRef]
  19. Hackos, D.H.; Chang, T.H.; Swartz, K.J. Scanning the intracellular S6 activation gate in the shaker K+ channel. J. Gen. Physiol. 2002, 119, 521–532. [Google Scholar] [CrossRef]
  20. Imbrici, P.; Grottesi, A.; D’Adamo, M.C.; Mannucci, R.; Tucker, S.J.; Pessia, M. Contribution of the central hydrophobic residue in the PXP motif of voltage-dependent K+ channels to S6 flexibility and gating properties. Channels 2009, 3, 39–45. [Google Scholar] [CrossRef] [Green Version]
  21. Döring, J.H.; Schröter, J.; Jüngling, J.; Biskup, S.; Klotz, K.A.; Bast, T.; Dietel, T.; Korenke, G.C.; Christoph, S.; Brennenstuhl, H.; et al. Refining Genotypes and Phenotypes in KCNA2-Related Neurological Disorders. Int. J. Mol. Sci. 2021, 22, 2824. [Google Scholar] [CrossRef]
  22. Syrbe, S.; Hedrich, U.B.S.; Riesch, E.; Djémié, T.; Müller, S.; Møller, R.S.; Maher, B.; Hernandez-Hernandez, L.; Synofzik, M.; Caglayan, H.S.; et al. De novo loss- or gain-of-function mutations in KCNA2 cause epileptic encephalopathy. Nat. Genet. 2015, 47, 393–399. [Google Scholar] [CrossRef]
  23. Faure, É.; Starek, G.; McGuire, H.; Bernèche, S.; Blunck, R. A limited 4 Å radial displacement of the S4-S5 linker is sufficient for internal gate closing in Kv channels. J. Biol. Chem. 2012, 287, 40091–40098. [Google Scholar] [CrossRef] [Green Version]
  24. Jensen, M.Ø.; Jogini, V.; Borhani, D.W.; Leffler, A.E.; Dror, R.O.; Shaw, D.E. Mechanism of voltage gating in potassium channels. Science 2012, 336, 229–233. [Google Scholar] [CrossRef]
  25. Catacuzzeno, L.; Sforna, L.; Franciolini, F. Voltage-dependent gating in K channels: Experimental results and quantitative models. Pflug. Arch. 2020, 472, 27–47. [Google Scholar] [CrossRef]
  26. Blunck, R.; Batulan, Z. Mechanism of electromechanical coupling in voltage-gated potassium channels. Front. Pharmacol. 2012, 3, 166. [Google Scholar] [CrossRef] [Green Version]
  27. Klassen, T.L.; Bomben, V.C.; Patel, A.; Drabek, J.; Chen, T.T.; Gu, W.; Zhang, F.; Chapman, K.; Lupski, J.R.; Noebels, J.L.; et al. High-resolution molecular genomic autopsy reveals complex sudden unexpected death in epilepsy risk profile. Epilepsia 2014, 55, e6–e12. [Google Scholar] [CrossRef] [Green Version]
  28. Labro, A.J.; Grottesi, A.; Sansom, M.S.; Raes, A.L.; Snyders, D.J. A Kv channel with an altered activation gate sequence displays both “fast” and “slow” activation kinetics. Am. J. Physiol. Cell Physiol. 2008, 294, C1476–C1484. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, Y.; Tachtsidis, G.; Schob, C.; Koko, M.; Hedrich, U.B.S.; Lerche, H.; Lemke, J.R.; van Haeringen, A.; Ruivenkamp, C.; Prescott, T.; et al. KCND2 variants associated with global developmental delay differentially impair Kv4.2 channel gating. Hum. Mol. Genet. 2021, 30, 2300–2314. [Google Scholar] [CrossRef]
  30. Imbrici, P.; Gualandi, F.; D’Adamo, M.C.; Masieri, M.T.; Cudia, P.; De Grandis, D.; Mannucci, R.; Nicoletti, I.; Tucker, S.J.; Ferlini, A.; et al. A novel KCNA1 mutation identified in an Italian family affected by episodic ataxia type 1. Neuroscience 2008, 157, 577–587. [Google Scholar] [CrossRef]
  31. Eunson, L.H.; Rea, R.; Zuberi, S.M.; Youroukos, S.; Panayiotopoulos, C.P.; Liguori, R.; Avoni, P.; McWilliam, R.C.; Stephenson, J.B.; Hanna, M.G.; et al. Clinical, genetic, and expression studies of mutations in the potassium channel gene KCNA1 reveal new phenotypic variability. Ann. Neurol. 2000, 48, 647–656. [Google Scholar] [CrossRef]
  32. Chou, S.M.; Li, K.X.; Huang, M.Y.; Chen, C.; Lin King, Y.H.; Li, G.G.; Zhou, W.; Teo, C.F.; Jan, Y.N.; Jan, L.Y.; et al. Kv1.1 channels regulate early postnatal neurogenesis in mouse hippocampus via the TrkB signaling pathway. Elife 2021, 10, e58779. [Google Scholar] [CrossRef] [PubMed]
  33. Lin King, Y.H.; Chen, C.; Lin King, J.V.; Simms, J.; Glasscock, E.; Yang, S.B.; Jan, Y.N.; Jan, L.Y. Kv1.1 preserves the neural stem cell pool and facilitates neuron maturation during adult hippocampal neurogenesis. Proc. Natl. Acad. Sci. USA 2022, 119, e2118240119. [Google Scholar] [CrossRef] [PubMed]
  34. Lauxmann, S.; Sonnenberg, L.; Koch, N.A.; Bosselmann, C.; Winter, N.; Schwarz, N.; Wuttke, T.V.; Hedrich, U.B.S.; Liu, Y.; Lerche, H.; et al. Therapeutic Potential of Sodium Channel Blockers as a Targeted Therapy Approach in KCNA1-Associated Episodic Ataxia and a Comprehensive Review of the Literature. Front. Neurol. 2021, 12, 703970. [Google Scholar] [CrossRef]
  35. Du, J.; Vegh, V.; Reutens, D.C. Persistent sodium current blockers can suppress seizures caused by loss of low-threshold D-type potassium currents: Predictions from an in silico study of Kv1 channel disorders. Epilepsia Open 2020, 5, 86–96. [Google Scholar] [CrossRef] [Green Version]
  36. Temperini, C.; Innocenti, A.; Scozzafava, A.; Parkkila, S.; Supuran, C.T. The coumarin-binding site in carbonic anhydrase accommodates structurally diverse inhibitors: The antiepileptic lacosamide as an example and lead molecule for novel classes of carbonic anhydrase inhibitors. J. Med. Chem. 2010, 53, 850–854. [Google Scholar] [CrossRef]
  37. Caraballo, R.; Reyes, G.; Demirdjian, G.; Huaman, M.; Gutierrez, R. Long-term use of cannabidiol-enriched medical cannabis in a prospective cohort of children with drug-resistant developmental and epileptic encephalopathy. Seizure 2022, 95, 56–63. [Google Scholar] [CrossRef]
  38. Pozzi, M.; Zanotta, N.; Epifanio, R.; Baldelli, S.; Cattaneo, D.; Clementi, E.; Zucca, C. Lacosamide effectiveness and tolerability in patients with drug-resistant epilepsy and severe disability under polytherapy: Therapy optimization as emerging from an observational study. Epilepsy Behav. 2022, 128, 108598. [Google Scholar] [CrossRef]
  39. Franco, V.; Barbieri, M.A.; Cutroneo, P.M.; Arena, I.; Cicala, G.; Marchione, P.; Spina, E.; Perucca, E. Pediatric adverse reactions to antiseizure medications—An analysis of data from the Italian spontaneous reporting system (2001–2019). Epilepsy Behav. 2021, 119, 107989. [Google Scholar] [CrossRef]
  40. Imbrici, P.; Altamura, C.; Gualandi, F.; Mangiatordi, G.F.; Neri, M.; De Maria, G.; Ferlini, A.; Padovani, A.; D’Adamo, M.C.; Nicolotti, O.; et al. A novel KCNA1 mutation in a patient with paroxysmal ataxia, myokymia, painful contractures and metabolic dysfunctions. Mol. Cell. Neurosci. 2017, 83, 6–12. [Google Scholar] [CrossRef]
  41. Sali, A.; Blundell, T.L. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 1993, 234, 779–815. [Google Scholar] [CrossRef]
  42. Webb, B.; Sali, A. Comparative Protein Structure Modeling Using MODELLER. Curr. Protoc. Bioinform. 2016, 54, 5.6.1–5.6.37. [Google Scholar] [CrossRef] [Green Version]
  43. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  44. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef]
  45. Phillips, J.C.; Hardy, D.J.; Maia, J.D.C.; Stone, J.E.; Ribeiro, J.V.; Bernardi, R.C.; Buch, R.; Fiorin, G.; Hénin, J.; Jiang, W.; et al. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 2020, 153, 044130. [Google Scholar] [CrossRef]
Figure 1. EEG during NREM sleep showing frequent and independent epileptiform abnormalities over bilateral frontotemporal regions.
Figure 1. EEG during NREM sleep showing frequent and independent epileptiform abnormalities over bilateral frontotemporal regions.
Ijms 23 08079 g001
Figure 2. (A) Position of identified mutations in the Kv1.1 channel S5–S6 regions, including the P403A variant. (B) Amino acid alignment of Kv1 channels.
Figure 2. (A) Position of identified mutations in the Kv1.1 channel S5–S6 regions, including the P403A variant. (B) Amino acid alignment of Kv1 channels.
Ijms 23 08079 g002
Figure 3. (AC) Representative current traces evoked by 200 ms depolarizing steps from a holding potential of −80 to +60 mV from Kv1.1 WT (A), P403A (B), and Kv1.1 WT+P403A (C) channels expressed in HEK 293 cells. The voltage protocol is indicated in the upper panel in (A). (D) Current–voltage relationship for Kv1.1 WT (5 μg), P403A (5 μg) and Kv1.1 WT+P403A (5 μg + 5 μg) channels (n = 14–36). (E) Voltage-dependent activation curves for Kv1.1 WT, P403A, and Kv1.1 WT+P403A channels were obtained by plotting the normalized peak tail currents measured at −50 or −30 mV as a function of the prepulse potentials and fitting data points with a Boltzmann function (n = 10–18 cells).
Figure 3. (AC) Representative current traces evoked by 200 ms depolarizing steps from a holding potential of −80 to +60 mV from Kv1.1 WT (A), P403A (B), and Kv1.1 WT+P403A (C) channels expressed in HEK 293 cells. The voltage protocol is indicated in the upper panel in (A). (D) Current–voltage relationship for Kv1.1 WT (5 μg), P403A (5 μg) and Kv1.1 WT+P403A (5 μg + 5 μg) channels (n = 14–36). (E) Voltage-dependent activation curves for Kv1.1 WT, P403A, and Kv1.1 WT+P403A channels were obtained by plotting the normalized peak tail currents measured at −50 or −30 mV as a function of the prepulse potentials and fitting data points with a Boltzmann function (n = 10–18 cells).
Ijms 23 08079 g003
Figure 4. Kinetics of activation (A) and deactivation (B) measured for Kv1.1 WT, P403A and Kv1.1 WT+P403A channels. The time constants, resulting from the fit of the activating and deactivating current traces with a single exponential function, were plotted as a function of voltage (n = 10–19 cells). (C) Bar graphs showing the time constants of the C-type inactivation for the indicated channels calculated by fitting current decay with a single exponential function (n = 9–11 cells). (D) Recovery from inactivation for the indicated channels.
Figure 4. Kinetics of activation (A) and deactivation (B) measured for Kv1.1 WT, P403A and Kv1.1 WT+P403A channels. The time constants, resulting from the fit of the activating and deactivating current traces with a single exponential function, were plotted as a function of voltage (n = 10–19 cells). (C) Bar graphs showing the time constants of the C-type inactivation for the indicated channels calculated by fitting current decay with a single exponential function (n = 9–11 cells). (D) Recovery from inactivation for the indicated channels.
Ijms 23 08079 g004
Figure 5. (A) Side view of the homology model of the Kv1.1 channel built upon the crystal structure of Kv1.2/2.1 chimera showing the localization of the P403 residue and P403A mutation. P403A is highlighted in magenta. (B) Loss of symmetry of wildtype and P403A measured at the top of S6 (position G385). (C) Measure of the S6 kink angle (left) and symmetry angle (right) for the four subunits of WT and P403A channels. The kink angle of the S6 was measured as the angle between the helix axis of the S6 above and below the PVP/A-motif. The symmetry angles are formed by the vectors from G385 of one subunit to the G385 of the two neighboring subunits. Values are given as mean ± SD of the values measured every 0.1 ns over the last 50 ns.
Figure 5. (A) Side view of the homology model of the Kv1.1 channel built upon the crystal structure of Kv1.2/2.1 chimera showing the localization of the P403 residue and P403A mutation. P403A is highlighted in magenta. (B) Loss of symmetry of wildtype and P403A measured at the top of S6 (position G385). (C) Measure of the S6 kink angle (left) and symmetry angle (right) for the four subunits of WT and P403A channels. The kink angle of the S6 was measured as the angle between the helix axis of the S6 above and below the PVP/A-motif. The symmetry angles are formed by the vectors from G385 of one subunit to the G385 of the two neighboring subunits. Values are given as mean ± SD of the values measured every 0.1 ns over the last 50 ns.
Ijms 23 08079 g005
Table 1. Biophysical parameters of Kv1.1 WT, P403A and Kv1.1 WT+P403A channels expressed in HEK 293 cells.
Table 1. Biophysical parameters of Kv1.1 WT, P403A and Kv1.1 WT+P403A channels expressed in HEK 293 cells.
Voltage Dependence
of Activation
Kinetics of
Activation
Kinetics of
Deactivation
C-Type InactivationRecovery from
Inactivation
V1/2 (mV)k (mV)τV1/2 (ms)τV1/2 (ms)τfast (s)τslow (s)τ (s)
Kv1.1 WT−21.5 ± 0.4
(29)
10.0 ± 0.8
(29)
9.4 ± 0.9
(15)
41.5 ± 2.4
(16)
4.9 ± 1.4
(10)
37.4 ± 6.7
(10)
3.4 ± 0.6
(7)
P403A7.8 ± 0.4 *
(13)
9.0 ± 0.8
(13)
31.7 ± 6.0 *
(14)
47.2 ± 2.6
(14)
5.8 ± 1.7
(7)
42.1 ± 15.7
(7)
2.2 ± 0.3
(9)
Kv1.1 WT+P403A−4.3 ± 0.7 *
(10)
12.0 ± 0.9
(10)
14.7 ± 3.1
(10)
41.5 ± 1.6
(7)
5.8 ± 0.9
(9)
40.2 ± 9.5
(9)
2.8 ± 0.3
(9)
Data are mean ± SE. The number of cells is indicated in brackets. * p < 0.05, with respect to WT.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dinoi, G.; Morin, M.; Conte, E.; Mor Shaked, H.; Coppola, M.A.; D’Adamo, M.C.; Elpeleg, O.; Liantonio, A.; Hartmann, I.; De Luca, A.; et al. Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia. Int. J. Mol. Sci. 2022, 23, 8079. https://doi.org/10.3390/ijms23158079

AMA Style

Dinoi G, Morin M, Conte E, Mor Shaked H, Coppola MA, D’Adamo MC, Elpeleg O, Liantonio A, Hartmann I, De Luca A, et al. Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia. International Journal of Molecular Sciences. 2022; 23(15):8079. https://doi.org/10.3390/ijms23158079

Chicago/Turabian Style

Dinoi, Giorgia, Michael Morin, Elena Conte, Hagar Mor Shaked, Maria Antonietta Coppola, Maria Cristina D’Adamo, Orly Elpeleg, Antonella Liantonio, Inbar Hartmann, Annamaria De Luca, and et al. 2022. "Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia" International Journal of Molecular Sciences 23, no. 15: 8079. https://doi.org/10.3390/ijms23158079

APA Style

Dinoi, G., Morin, M., Conte, E., Mor Shaked, H., Coppola, M. A., D’Adamo, M. C., Elpeleg, O., Liantonio, A., Hartmann, I., De Luca, A., Blunck, R., Russo, A., & Imbrici, P. (2022). Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia. International Journal of Molecular Sciences, 23(15), 8079. https://doi.org/10.3390/ijms23158079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop