Next Article in Journal
Cannabis, the Endocannabinoid System and Immunity—the Journey from the Bedside to the Bench and Back
Next Article in Special Issue
Computational SNP Analysis and Molecular Simulation Revealed the Most Deleterious Missense Variants in the NBD1 Domain of Human ABCA1 Transporter
Previous Article in Journal
The Strategic Alliance between Clinical and Molecular Science in the War against SARS-CoV-2, with the Rapid-Diagnostics Test as an Indispensable Weapon for Front Line Doctors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Case Report

Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome

by
Pedro A. Lazo
1,2,*,
Juan L. García
1,2,
Paulino Gómez-Puertas
3,
Íñigo Marcos-Alcalde
3,4,
Cesar Arjona
5,6,
Alvaro Villarroel
7,
Rogelio González-Sarmiento
1,2,8 and
Carmen Fons
5,9
1
Molecular Mechanisms of Cancer Program, Instituto de Biología Molecular y Celular del Cáncer, Consejo Superior de Investigaciones Científicas (CSIC), Universidad de Salamanca, 30007 Salamanca, Spain
2
Instituto de Investigación Biomédica de Salamanca (IBSAL), Hospital Universitario de Salamanca, 30007 Salamanca, Spain
3
Centro de Biología Molecular Severo Ochoa, CSIC-Universidad Autónoma de Madrid, Cantoblanco, E-28049 Madrid, Spain
4
Biosciences Research Institute, School of Experimental Sciences, Universidad Francisco de Vitoria, 28223 Madrid, Spain
5
Institut de Recerca Sant Joan de Déu, Esplugues de Llobregat, 08950 Barcelona, Spain
6
Instituto Pediátrico de Enfermedades Raras (IPER), Hospital Sant Joan de Déu, 08950 Barcelona, Spain
7
Instituto de Biofísica, Consejo Superior de Investigaciones Científicas (CSIC), Universidad del País Vasco, 48940 Bilbao, Spain
8
Unidad de Genética Molecular, Departamento de Medicina, Universidad de Salamanca, 37008 Salamanca, Spain
9
Neurology Department, Hospital Sant Joan de Déu, Sant Joan de Déu Research Institute and CIBERER, Instituto de Salud Carlos III, 08950 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(12), 4447; https://doi.org/10.3390/ijms21124447
Submission received: 20 May 2020 / Revised: 17 June 2020 / Accepted: 18 June 2020 / Published: 23 June 2020
(This article belongs to the Special Issue Macromolecular Modeling to Understand Genetic Disorders)

Abstract

:
Complex neurodevelopmental syndromes frequently have an unknown etiology, in which genetic factors play a pathogenic role. This study utilizes whole-exome sequencing (WES) to examine four members of a family with a son presenting, since birth, with epileptic-like crises, combined with cerebral palsy, severe neuromotor and developmental delay, dystonic tetraparexia, axonal motor affectation, and hyper-excitability of unknown origin. The WES study detected within the patient a de novo heterozygous in-frame duplication of thirty-six nucleotides within exon 7 of the human KCNQ2 gene. This insertion duplicates the first twelve amino acids of the calmodulin binding site I. Molecular dynamics simulations of this KCNQ2 peptide duplication, modelled on the 3D structure of the KCNQ2 protein, suggest that the duplication may lead to the dysregulation of calcium inhibition of this protein function.

1. Introduction

Newborns with severe non-hereditary neuromotor and neurodevelopmental delays, cerebral palsy, and epileptic crises often have an unknown etiology, which are likely to have a heterogeneous origin, in which genetic alterations can play a relevant pathogenic role.
Voltage-gated sodium channels generate neuronal action potentials [1]. Voltage-gated potassium channels re-establish these potentials back to normal [2], and the functional coordination of these channels regulate the neuronal excitability, controlling neuronal networks [3]. These two types of ion channels have specific subcellular locations on neuron membranes [4]. The pathogenesis of severe early neuromotor developmental delays and epilepsy is associated with mutations or variants in several genes. Among them, KCNQ2 forms dimers with KCNQ3 in order to assemble the potassium M channel [5], which limits repetitive firing [6,7], controlling the excitability of neurons [8] and their responsiveness to synaptic inputs [9]. Human KCNQ2 pathogenic variants have been associated to benign familial neonatal epilepsy/infantile epileptic encephalopathy type 7 (BFNS, OMIM: 121200) and to early onset epileptic encephalopathy (EOEE; OMIM: 613720), but as the number of cases increases, the associated clinical phenotypes are becoming more heterogeneous. Known KCNQ2 pathogenic variants are distributed throughout the length of the protein, with two-thirds located within trans-membrane domains. However, there is no correlation between the location of the variant amino acid on the protein sequence, the type of pathogenic variant, and the severity of the neurological phenotype. This variability suggests that KCNQ2 alterations are necessary for epileptic syndromes [10]. However, the severity and variability of the associated full neurological phenotype is very likely determined by the contribution of additional genetic alterations or variants not yet identified [10,11,12]. The clinical heterogeneity of these neurological syndromes reflects a very complex situation in which an unknown combination of alterations in gene dosage and pathogenic variants may contribute to the pathogenesis of complex neurodevelopmental phenotypes, including epileptic disorders [11,13,14,15], cerebral palsy [16,17], and autism [18,19], among others.

2. Results

2.1. Epileptic, Cerebral Palsy, and Neuromotor Developmental Delay Phenotype of Patient

The case is a 6-year-old boy presenting with epileptic encephalopathy of neonatal onset, cerebral palsy, and a severe neuromotor delay of unknown origin. Pregnancy with gestational diabetes was controlled with diet. Delivery was uneventful. After 48 h of life, the boy presented with episodes of cyanosis, generalized hypertonia, and tonic asymmetric postures followed by apnea. Video-EEG at 5 days of life showed bilateral and asynchronous spike-and-wave. Seizures were refractory to phenobarbital, but were controlled with phenytoin. At 3 months of age, the boy presented with startle episodes without EEG correlation related to the wake and sleep transition, together with dystonic postures. One month later, epileptic spasms without hypsarrhythmia were observed. There was no response to levetiracetam and valproic acid, but spasms stopped after lacosamide (LCS) treatment was initiated (11 months of age). Metabolic tests (including cerebro-spinal fluid(CSF) studies), karyotype, and brain MRI were normal. The progressive increase of Thyroid stimulating hormone (TSH) was treated with oral L-T4. At 20 months, the boy presented with a recurrence of spasms that were controlled with LCS dose adjustment. After several years of good seizure control, LCS was discontinued. The boy suffered from the decompensation of seizures with a febrile viral infection, and at 5 years of age he presented with startle episodes associated to auditory stimulus without EEG correlation, whereas during sleep he presented with epileptic spasms with EEG correlation. Treatment with vigabatrin was initiated with a clinical and EEG improvement. Video-EEG showed abnormal background activity (Figure 1A) and paroxysmal multifocal discharges that showed activation and generalization during sleep (Figure 1B). The final video-EEG monitoring showed bilateral polyspike discharges and clinical spasms in the polygraph recording, which did not show any clear EEG correlation.
Many KCNQ2 encephalopathy patients present a very good response to sodium channel blockers, as our patient did with phenytoin and LCS during the neonatal period, and later when treatment with LCS was reintroduced.
During a physical exam, the five-year old child presented with a severe axial hypotonia without head control, hypertonic movements of four limbs with dystonic movements of the upper limbs, no hand use, and was not able to sit or crawl. The child had a dystonic tetraparesia and peripheral neuro-axonal motor affectation. The child had no development of expressive language.

2.2. Partial de Novo KCNQ2 Exon 7 Duplication

This complex and severe neurodevelopmental syndrome could be due to the phenotypic convergence of several genetic alterations whose biological functions contribute to different aspects of this complex clinical phenotype. The initial approach used to identify them was a whole-exome sequencing (WES) study.
WES detected in heterozygosis a de novo alteration, not present in the parents, in the KCNQ2 gene that was confirmed by Sanger sequencing. This genetic alteration presented as an in-frame sequence duplication (DUP) of thirty-six nucleotides within exon 7 (GRCh37/hg19: chr20:g.62070007_62070008dupCTTCTCAAAGTGCTTCTGCCTGTGCTGCTCCTGAAC; c.958_993dupGTTCAGGAGCAGCACAGGCAGAAGCACTTTGAGAAG), (GRCh38/hg38: chr20:g.63438654_63438655dupCTTCTCAAAGTGCTTCTGCCTGTGCTGCTCCTGAAC), which resulted in the addition of twelve extra amino acids in the protein (NC_000020.11:g.63438659_63438694dup; p.Val320_Lys331 (ValGlnGluGlnHisArgGlnLysHisPheGluLys). This peptide duplication generated a tandem repeat of the first twelve residues (320–331) of the calmodulin binding site 1 (321–358) [20] and was located at the start of the cytosolic regions of KCNQ2, immediately after the end of the last transmembrane domain of the KCNQ2 protein (Figure 2). This peptide duplication was considered pathogenic by Polyphen-2 [21] and VarElect [22] program analysis; it was also considered pathogenic according to the ACMG criteria [23].
The location of the thirty-six nucleotide duplication was confirmed by Sanger sequencing (Supplementary Figure S1). Thus, this case represents the first association of this very rare type of variant to a neurological phenotype. All the genes identified by WES in the patient are listed in Table S1. The WES analysis also detected a RARS2 mutation in the heterozygosis that codes for the mutant p.S443P, which was inherited from the father (Supplementary Figure S2). RARS2 codes for an arginyl tRNA synthetase, and RARS2 mutations have been associated to pontocerebellar hypoplasia [24].

2.3. Dynamic Molecular Modelling of KCNQ2 Partial Duplication

This peptide insertion can affect the ion accessibility to the entrance of the channel by altering its properties, or alternatively, altering its regulation by calmodulin (CaM) and calcium. An increase in calcium resulted in its binding to CaM, which was responsible for the allosteric regulation of the transmembrane region of KCNQ2, which modulated the closing of the K+ channel [25,26,27]. The effect of the VQEQHRQKHFEK peptide duplication on the activity of the KCNQ2 channel was studied. For this, the structure modeling and molecular dynamics simulation was performed using 6FEH and 6FEG crystal three-dimensional structures [27]. The model of the duplicated peptide sequence suggests that the inserted amino acids form a loop in the helix hA of the CaM binding domain of KCNQ2, located in the vicinity of the C-lobe of CaM. Although the inserted sequence is located in the CaM binding domain, its presence does not represent a steric impediment to the interaction between the cytosolic domain of the channel and the CaM monomer. It has recently been described that a small change in the orientation of the helices hA and hB of KCNQ2 caused by a structural movement of CaM associated with the binding of two Ca++ ions in its C-lobe could be responsible for the allosteric regulation of the channel in response to an increase in the concentration of cytosolic calcium [20]. Dynamic modelling of the wild-type and variant KCNQ2 (DUP-KCNQ2) proteins was performed in response to increasing concentrations of calcium. The modelling was applied to the structural changes that the KCNQ2 protein underwent, mediated by CaM, when the calcium concentration was raised from basal cytosolic values (10–100 nM), in which Ca++ ions were bound only to the N-lobe of CaM, to high Ca++ values (1–10 μM) that also triggered the association of Ca++ to the CaM C-lobe. The position of hA and hB helices in WT-KCNQ2 modelled at high Ca2++ levels (Figure 3A) remained unchanged after 40 ns of free MD simulation (Figure 3B). In the case of the DUP-KCNQ2 model, at basal low Ca++ (Figure 3C), the organization was essentially identical to the WT- KCNQ2 model (Figure 3A). However, the presence of the duplicated sequence (green loop in Figure 3C) promoted a significant variation in the position of helix hA after 10 ns of free MD simulation (Figure 3E, black line) that was stably maintained after 40 ns of MD (Figure 3D, red). This new arrangement was equivalent to the structure of the same hA helix when Ca++ ions are bound only to the N-lobe of CaM [27]. In fact, the value of the root mean square deviation (rmsd), measured for helix hA under both conditions was 1.72 Å, indicating a virtually complete superposition. Since it has been proposed that the position of hA is responsible for the transmission of the high Ca++ signal from CaM to the K+ channel [27], this behavior predicts that this pathogenic KCNQ2 with the peptide insertion cannot respond to the increase in Ca++, and thus interferes with its regulation by CaM. This arrangement suggests that the variant protein, even in the presence of high concentrations of cytosolic calcium, can maintain a conformation equivalent to what it would have in the presence of basal calcium concentrations, and therefore would be functionally equivalent to a lack of response to calcium changes. Because this domain allosterically regulates the closure of the channel at high calcium concentrations, it suggests that within the variant protein, the channel could remain constitutively open, and would be unable to respond to the increase in calcium concentration in the cytoplasm and prevent its inhibitory effects, which would in turn lead to neuronal hyper-excitability.

3. Discussion

Known KCNQ2 pathogenic variants were very heterogeneous, and were distributed throughout the length of the protein (Figure 2). Pathogenic variants located within the transmembrane-pore domain altered the in- or out-flux of potassium ions, or their timing, resulting in a defective potassium transport. Other KCNQ2 pathogenic variants affected its interaction or regulation with calmodulin, or alternatively, facilitated the retention of the KCNQ2 protein in the endoplasmic reticulum, affecting its subcellular localization and leading to its reduction on the neuron surface [28]. The heterogeneity of KCNQ2 pathogenic variants suggests that what is important is the impairment of the normal KCNQ2 channel function via different mechanisms. KCNQ2 pathogenic variants alter, either quantitatively or qualitatively, the correct opening and closing of this potassium channel, resulting in a loss of function [27]. Despite their common loss of function, the associated neurological phenotype is clinically very heterogeneous. Therefore, whether the channel is locked in an open conformation or alters the balance between the in- and out-fluxes of potassium, the result is an alteration of its dampening signal on the sodium channel excitatory effects and a gain of neuronal excitability via the hyper-activation of sodium channels [28]. A reduction in the surface expression of KCNQ2 fails to suppress neuronal excitability [29]. Even KCNQ2 pathogenic variants having opposite functional effects on M potassium channels can produce similar EOEE neurological phenotypes [30]. Recently, a deletion in KCNQ2 was reported (c.913_915del [p.Phe305del)]) [31], which partially resembles the epileptic phenotype of this case, but without an associated cerebral palsy and distal neuromotor phenotype [31]. In murine heterozygous KCNQ2-knockout mice there was a lower seizure threshold [32]. KCNQ2-null mice have an increased excitability of primary sensory neurons, but do not have neuromotor alterations [33], and the inhibition of the Nav1.6 sodium channel counteracts the neuronal hyperexcitability [34]. Moreover, the replacement in murine embryos of KCNQ2 with the I250V mutant resulted in an increased neuronal excitability of pyramidal neurons [34].
There was a large KCNQ2 pathogenic variant heterogeneity in patients: 73 different mutated residues, including insertions and deletions, were detected in 172 mutant cases (Figure 2). However, there was no correlation between the locations of the pathogenic variants, despite their common loss of function, and the severity of the neurological phenotype. Moreover, there was no specific individual pathogenic variant associated to any specific phenotype, and the same pathogenic variant could be associated to different phenotypes, benign or severe, depending on the individual. The need to have additional genetic alterations to cause complex neurological phenotypes could explain why in cases in which the (c.628C>T; p.Arg210Cys) KCNQ2 pathogenic variant is inherited, the carrier parent has no, or mild, pathology and the affected offspring presents a very severe epileptic neurodevelopmental delay syndrome [35]. In a more recent study of novel mutations, there was no correlation between the location of mutation and the severity of phenotype [36]. This suggests that KCNQ2 alterations are necessary, but the severity and complexity of the neurological phenotype in each patient are a likely consequence of the contribution of additional gene variations in the patient genotype. The KCNQ2 pathogenic variant needs an environment in order to manifest [10]. A recent review on KCNQ2 [10] indicates the need for additional genetic events to explain the variability in the phenotype, which are currently unknown [10]. This means that both inherited and de novo KCNQ2 pathogenic variants require the cooperation of additional genetic alterations in order to be able to cause the full spectrum of the clinical phenotype [37].
We conclude that the genetic heterogeneity of early epileptic encephalopathy, cerebral palsy, and severe neurodevelopmental delays will have to be characterized in the context of the initiating pathogenic variant that is modulated by several additional genetic changes associated to neuronal functions, which remains to be identified. When unraveled, they might share common pathogenic pathways, although they might involve different genes. The identification of gene/protein networks associated to the clinical neurological phenotypes could set the basis for designing novel therapeutic approaches for the management of these patients, and minimize, or compensate, the neurological consequences of these pathogenic gene combinations.

4. Methods

4.1. Standard Protocol Approvals, Registrations, and Patient Consents

The genomic studies were performed for the diagnosis of a complex neurological syndrome of unknown origin. Total DNA was obtained from peripheral blood samples and used to perform the genomic diagnostic study. Both parents gave their informed consent for the genomic study.

4.2. Whole-Exome Sequencing

Whole-exome analysis (WES) was performed in the four members of the family with the affected boy as described in the Supplementary Methods. Briefly, 50 Mbp of captured exome content was selected for sequencing using Agilent SureSelect v5 and following the supplier’s Agilent SureSelect enrichment Protocol Guide. The generated DNA fragments (DNA libraries) were sequenced in the Illumina Hiseq 2000 platform with a sequence read length of 100 bp (Paired-End), obtaining an average coverage depth of 100 x in the target region. The Genome Reference Consortium Human Builds 37 and 38 (GRCh37 and GRCh38) were used as reference sequences for the identification of changes in the four family members. The exome analysis confirmed the paternity of the children. Sanger sequencing was used to confirm the mutations.

4.3. Analysis of Gene Alterations and Clinical Phenotype

All gene alterations or variants detected by WES were analyzed using the VarElect program [22] (LifeMap Sciences Inc., Tel Aviv, Israel) to search for correlations between the gene/protein biological function and the clinical phenotypes of the case, in order to detect a gene’s contribution to clinical symptoms.

4.4. Molecular Dynamics Simulation and Structural Modeling of KCNQ2

The 3D structures of the calmodulin (CaM) binding domain from the wild-type KCNQ2 channel (isoform 1, UniProtKB ID: O43526-1), and from the variant protein containing a duplicated VQEQHRQKHFEK sequence (residues 320–331) were modeled. As a template, the structure of the CaM binding domain of the human KCNQ2 channel was solved in the presence of CaM in the saturated calcium state (holoCaM/KCNQ2-hAB complex) or in the intermediate calcium state (intCaM/KCNQ2-hAB complex). The three-dimensional X-ray crystal structures of KCNQ2, 6FEH, and 6FEG are available in the Protein Data Bank [27].
Models of normal KCNQ2 and its peptide-duplication mutant were built using the SWISS-MODEL server (http://swissmodel.expasy.org), and their structural quality was within the range of those accepted for homology-based structures (Anolea/Gromos/QMEAN4). Modeled structures were subjected to MD simulation using the AMBER14 molecular dynamics package (University of California, San Francisco), as previously reported [38].
The 3D structures were solved with periodic octahedral pre-equilibrated solvent boxes using the LEaP module of AMBER, with 12 Å as the shortest distance between any atom in the protein and the periodic box boundaries. Free MD simulations were performed essentially as previously described [38], using the PMEMD program of AMBER and the ff14SB force field (University of California San Francisco). The SHAKE algorithm was used, allowing for a time step of 2 fs. Systems were initially relaxed over 15,000 steps of energy minimization with a cut-off of 12 Å. Simulations were then started with a 20 ps heating phase, raising the temperature from 0 to 300 K in 10 temperature-change steps, after each of which velocities were reassigned. During minimization and heating, the Cα trace dihedrals were restrained with a force constant of 500 kcal mol−1 rad−2 and gradually released in an equilibration phase in which the force constant was progressively reduced to 0 over 200 ps. After the equilibration phase, 40 ns of unrestricted MD simulations were obtained for wild-type and variant structures. Root mean square deviation (rmsd) of Cα trace was monitored along the trajectories, and was measured every 20 ps. MD trajectories were analyzed using VMD software [39]. Figures were generated using the Pymol Molecular Graphics System (Schrödinger, LLC).

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/12/4447/s1.

Author Contributions

P.A.L. planned and coordinated the study, interpreted data, and wrote the manuscript. J.L.G. performed and analyzed the aCGH study. P.G.-P. and C.A. analyzed the WES study. P.G.-P., Í.M.-A., and A.V. analyzed the structural data molecular modelling. R.G.-S. interpreted the data and performed validation. C.F. coordinated and interpreted the neurological studies. All authors have read and agreed to the published version of the manuscript.

Funding

Grants from Agencia Estatal de Investigación-Ministerio de Economía y Competitividad-FEDER-Fondo Social Europeo (SAF2016-75744-R, PID2019-105610RB-I00) and Consejería de Educación-Junta de Castilla y León (CLC-2017-01 and UIC-258) to P.A.L. Ministerio de Economía, Industria y Competitividad (CTQ2015-68756-R, BFU2015-66910-R) to A.V. Ministerio de Economía y Competitividad-Fondo Social Europeo Grant from Instituto de Salud Carlos III-Fondo Social Europeo (PI16/01920) to R. G-S. Agencia Estatal de Investigación-Ministerio de Economía y Competitividad-FEDER-Fondo Social Europeo (RTC-2017-6494-1 and RTI2018-094434-B-I00) to P. G-P. Financial support of the Departament de Salut de la Generalitat de Catalunya-URDCat Project (SLT002/16/00174) to C.A. and C.F.

Acknowledgments

The computational support of the “Centro de Computación Científica CCC-UAM” is gratefully recognized.

Conflicts of Interest

The authors declare no conflict of interest.

Data availability statement

The KCNQ2 partial duplication has been deposited in the ClinVar database. (https://www.ncbi.nlm.nih.gov/clinvar/variation/617505). Whole Exome Sequence (WES) raw data files are available in: https://www.ncbi.nlm.nih.gov/sra/PRJNA629061 https://www.ncbi.nlm.nih.gov/bioproject/PRJNA629061 https://digital.csic.es/handle/10261/170405.

References

  1. Wang, J.; Ou, S.W.; Wang, Y.J. Distribution and function of voltage-gated sodium channels in the nervous system. Channels (Austin) 2017, 11, 534–554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Smith, R.S.; Walsh, C.A. Ion Channel Functions in Early Brain Development. Trends Neurosci. 2020, 43, 103–114. [Google Scholar] [CrossRef] [PubMed]
  3. Battefeld, A.; Tran, B.T.; Gavrilis, J.; Cooper, E.C.; Kole, M.H. Heteromeric Kv7.2/7.3 channels differentially regulate action potential initiation and conduction in neocortical myelinated axons. J. Neurosci. 2014, 34, 3719–3732. [Google Scholar] [CrossRef] [PubMed]
  4. Misonou, H. Precise localizations of voltage-gated sodium and potassium channels in neurons. Dev. Neurobiol. 2018, 78, 271–282. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, H.S.; Pan, Z.; Shi, W.; Brown, B.S.; Wymore, R.S.; Cohen, I.S.; Dixon, J.E.; McKinnon, D. KCNQ2 and KCNQ3 potassium channel subunits: Molecular correlates of the M-channel. Science 1998, 282, 1890–1893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Weckhuysen, S.; Mandelstam, S.; Suls, A.; Audenaert, D.; Deconinck, T.; Claes, L.R.; Deprez, L.; Smets, K.; Hristova, D.; Yordanova, I.; et al. KCNQ2 encephalopathy: Emerging phenotype of a neonatal epileptic encephalopathy. Ann. Neurol. 2012, 71, 15–25. [Google Scholar] [CrossRef]
  7. Cooper, E.C.; Harrington, E.; Jan, Y.N.; Jan, L.Y. M channel KCNQ2 subunits are localized to key sites for control of neuronal network oscillations and synchronization in mouse brain. J. Neurosci. 2001, 21, 9529–9540. [Google Scholar] [CrossRef] [Green Version]
  8. Lombardo, J.; Sun, J.; Harrington, M.A. Rapid activity-dependent modulation of the intrinsic excitability through up-regulation of KCNQ/Kv7 channel function in neonatal spinal motoneurons. PLoS ONE 2018, 13, e0193948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Perez-Ramirez, M.B.; Laville, A.; Tapia, D.; Duhne, M.; Lara-Gonzalez, E.; Bargas, J.; Galarraga, E. KV7 Channels Regulate Firing during Synaptic Integration in GABAergic Striatal Neurons. Neural. Plast. 2015, 2015, 472676. [Google Scholar] [CrossRef]
  10. Lee, I.C.; Yang, J.J.; Li, S.Y. KCNQ2-Associated Epilepsy: A Review of Variable Phenotypes and Neurodevelopmental Outcomes. Neuropsychiatry 2018, 8, 318–323. [Google Scholar] [CrossRef]
  11. Mefford, H.C.; Muhle, H.; Ostertag, P.; von Spiczak, S.; Buysse, K.; Baker, C.; Franke, A.; Malafosse, A.; Genton, P.; Thomas, P.; et al. Genome-wide copy number variation in epilepsy: Novel susceptibility loci in idiopathic generalized and focal epilepsies. PLoS Genet. 2010, 6, e1000962. [Google Scholar] [CrossRef] [PubMed]
  12. Fahey, M.C.; Maclennan, A.H.; Kretzschmar, D.; Gecz, J.; Kruer, M.C. The genetic basis of cerebral palsy. Dev. Med. Child. Neurol. 2017, 59, 462–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rosenfeld, J.A.; Patel, A. Chromosomal Microarrays: Understanding Genetics of Neurodevelopmental Disorders and Congenital Anomalies. J. Pediatric Genet. 2017, 6, 42–50. [Google Scholar] [CrossRef]
  14. Lindy, A.S.; Stosser, M.B.; Butler, E.; Downtain-Pickersgill, C.; Shanmugham, A.; Retterer, K.; Brandt, T.; Richard, G.; McKnight, D.A. Diagnostic outcomes for genetic testing of 70 genes in 8565 patients with epilepsy and neurodevelopmental disorders. Epilepsia 2018, 59, 1062–1071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Borlot, F.; Regan, B.M.; Bassett, A.S.; Stavropoulos, D.J.; Andrade, D.M. Prevalence of Pathogenic Copy Number Variation in Adults With Pediatric-Onset Epilepsy and Intellectual Disability. JAMA Neurol. 2017, 74, 1301–1311. [Google Scholar] [CrossRef]
  16. McMichael, G.; Bainbridge, M.N.; Haan, E.; Corbett, M.; Gardner, A.; Thompson, S.; van Bon, B.W.; van Eyk, C.L.; Broadbent, J.; Reynolds, C.; et al. Whole-exome sequencing points to considerable genetic heterogeneity of cerebral palsy. Mol. Psychiatry 2015, 20, 176–182. [Google Scholar] [CrossRef]
  17. Emrick, L.T.; DiCarlo, S.M. The Expanding Role of Genetics in Cerebral Palsy. Phys. Med. Rehabil. Clin. N. Am. 2020, 31, 15–24. [Google Scholar] [CrossRef]
  18. Risch, N.; Spiker, D.; Lotspeich, L.; Nouri, N.; Hinds, D.; Hallmayer, J.; Kalaydjieva, L.; McCague, P.; Dimiceli, S.; Pitts, T.; et al. A genomic screen of autism: Evidence for a multilocus etiology. Am. J. Hum. Genet. 1999, 65, 493–507. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, X.; Shimada, T.; Otowa, T.; Wu, Y.Y.; Kawamura, Y.; Tochigi, M.; Iwata, Y.; Umekage, T.; Toyota, T.; Maekawa, M.; et al. Genome-wide Association Study of Autism Spectrum Disorder in the East Asian Populations. Autism. Res. 2015, 9, 340–349. [Google Scholar] [CrossRef]
  20. Yus-Najera, E.; Santana-Castro, I.; Villarroel, A. The identification and characterization of a noncontinuous calmodulin-binding site in noninactivating voltage-dependent KCNQ potassium channels. J. Biol. Chem. 2002, 277, 28545–28553. [Google Scholar] [CrossRef] [Green Version]
  21. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Stelzer, G.; Plaschkes, I.; Oz-Levi, D.; Alkelai, A.; Olender, T.; Zimmerman, S.; Twik, M.; Belinky, F.; Fishilevich, S.; Nudel, R.; et al. VarElect: The phenotype-based variation prioritizer of the GeneCards Suite. BMC Genom. 2016, 17 (Suppl. 2), 444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Green, R.C.; Berg, J.S.; Grody, W.W.; Kalia, S.S.; Korf, B.R.; Martin, C.L.; McGuire, A.L.; Nussbaum, R.L.; O’Daniel, J.M.; Ormond, K.E.; et al. ACMG recommendations for reporting of incidental findings in clinical exome and genome sequencing. Genet. Med. 2013, 15, 565–574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Paciorkowski, A.R.; Darras, B.T. Making sense of genetic heterogeneity: Emergence of pathways in developmental brain disorders. Neurology 2013. [Google Scholar] [CrossRef]
  25. Hernandez, C.C.; Zaika, O.; Tolstykh, G.P.; Shapiro, M.S. Regulation of neural KCNQ channels: Signalling pathways, structural motifs and functional implications. J. Physiol. 2008, 586, 1811–1821. [Google Scholar] [CrossRef]
  26. Alaimo, A.; Alberdi, A.; Gomis-Perez, C.; Fernandez-Orth, J.; Bernardo-Seisdedos, G.; Malo, C.; Millet, O.; Areso, P.; Villarroel, A. Pivoting between calmodulin lobes triggered by calcium in the Kv7.2/calmodulin complex. PLoS ONE 2014, 9, e86711. [Google Scholar] [CrossRef] [Green Version]
  27. Bernardo-Seisdedos, G.; Nunez, E.; Gomis-Perez, C.; Malo, C.; Villarroel, A.; Millet, O. Structural basis and energy landscape for the Ca(2+) gating and calmodulation of the Kv7.2 K(+) channel. Proc. Natl. Acad. Sci. USA 2018, 115, 2395–2400. [Google Scholar] [CrossRef] [Green Version]
  28. Abidi, A.; Devaux, J.J.; Molinari, F.; Alcaraz, G.; Michon, F.X.; Sutera-Sardo, J.; Becq, H.; Lacoste, C.; Altuzarra, C.; Afenjar, A.; et al. A recurrent KCNQ2 pore mutation causing early onset epileptic encephalopathy has a moderate effect on M current but alters subcellular localization of Kv7 channels. Neurobiol. Dis. 2015, 80, 80–92. [Google Scholar] [CrossRef]
  29. Kim, E.C.; Zhang, J.; Pang, W.; Wang, S.; Lee, K.Y.; Cavaretta, J.P.; Walters, J.; Procko, E.; Tsai, N.P.; Chung, H.J. Reduced axonal surface expression and phosphoinositide sensitivity in Kv7 channels disrupts their function to inhibit neuronal excitability in Kcnq2 epileptic encephalopathy. Neurobiol. Dis. 2018, 118, 76–93. [Google Scholar] [CrossRef]
  30. Devaux, J.; Abidi, A.; Roubertie, A.; Molinari, F.; Becq, H.; Lacoste, C.; Villard, L.; Milh, M.; Aniksztejn, L. A Kv7.2 mutation associated with early onset epileptic encephalopathy with suppression-burst enhances Kv7/M channel activity. Epilepsia 2016, 57, e87–e93. [Google Scholar] [CrossRef] [Green Version]
  31. Spagnoli, C.; Salerno, G.G.; Iodice, A.; Frattini, D.; Pisani, F.; Fusco, C. KCNQ2 encephalopathy: A case due to a de novo deletion. Brain Dev. 2018, 40, 65–68. [Google Scholar] [CrossRef] [PubMed]
  32. Otto, J.F.; Yang, Y.; Frankel, W.N.; White, H.S.; Wilcox, K.S. A spontaneous mutation involving Kcnq2 (Kv7.2) reduces M-current density and spike frequency adaptation in mouse CA1 neurons. J. Neurosci. 2006, 26, 2053–2059. [Google Scholar] [CrossRef] [PubMed]
  33. King, C.H.; Lancaster, E.; Salomon, D.; Peles, E.; Scherer, S.S. Kv7.2 regulates the function of peripheral sensory neurons. J. Comp. Neurol. 2014, 522, 3262–3280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Niday, Z.; Hawkins, V.E.; Soh, H.; Mulkey, D.K.; Tzingounis, A.V. Epilepsy-Associated KCNQ2 Channels Regulate Multiple Intrinsic Properties of Layer 2/3 Pyramidal Neurons. J. Neurosci. 2017, 37, 576–586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hewson, S.; Puka, K.; Mercimek-Mahmutoglu, S. Variable expressivity of a likely pathogenic variant in KCNQ2 in a three-generation pedigree presenting with intellectual disability with childhood onset seizures. Am. J. Med. Genet. A 2017, 173, 2226–2230. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, I.C.; Chang, T.M.; Liang, J.S.; Li, S.Y. KCNQ2 mutations in childhood nonlesional epilepsy: Variable phenotypes and a novel mutation in a case series. Mol. Genet. Genom. Med. 2019, 7, e00816. [Google Scholar] [CrossRef] [Green Version]
  37. Perucca, P.; Bahlo, M.; Berkovic, S.F. The Genetics of Epilepsy. Annu Rev. Genom. Hum. Genet. 2020, 20. [Google Scholar] [CrossRef]
  38. Marcos-Alcalde, I.; Mendieta-Moreno, J.I.; Puisac, B.; Gil-Rodriguez, M.C.; Hernandez-Marcos, M.; Soler-Polo, D.; Ramos, F.J.; Ortega, J.; Pie, J.; Mendieta, J.; et al. Two-step ATP-driven opening of cohesin head. Sci. Rep. 2017, 7, 3266. [Google Scholar] [CrossRef] [Green Version]
  39. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38, 27–38. [Google Scholar] [CrossRef]
Figure 1. Electroencephalogram (EEG) of the patient at six years of age. (A) Interictal EEG. Bilateral polyspike-and-wave discharges with frontal predominance. (B) Multifocal epileptiform abnormalities.
Figure 1. Electroencephalogram (EEG) of the patient at six years of age. (A) Interictal EEG. Bilateral polyspike-and-wave discharges with frontal predominance. (B) Multifocal epileptiform abnormalities.
Ijms 21 04447 g001
Figure 2. Partial duplication of exon 7 KCNQ2 in the calmodulin interaction site. The twelve duplicated amino acids are underlined, and the inserted peptide is indicated in red. In the sequence, the amino acids in yellow indicate the 73 residues affected by the pathogenic variants reported in 159 different patients obtained from the ClinVar database and publications referenced in the text. The human KCNQ2 protein reference is UniProtKB (O43526-1). The triangles indicate short deletions or insertions that caused a change in reading frame and a truncated protein. Information in blue boxes indicates the transmembrane domain (TM). The pore region is indicated with the black line. Information in purple indicates the two regions that interacted with calmodulin (CaM) that have two alpha helixes (α-hA in green and α-hB in dark purple).
Figure 2. Partial duplication of exon 7 KCNQ2 in the calmodulin interaction site. The twelve duplicated amino acids are underlined, and the inserted peptide is indicated in red. In the sequence, the amino acids in yellow indicate the 73 residues affected by the pathogenic variants reported in 159 different patients obtained from the ClinVar database and publications referenced in the text. The human KCNQ2 protein reference is UniProtKB (O43526-1). The triangles indicate short deletions or insertions that caused a change in reading frame and a truncated protein. Information in blue boxes indicates the transmembrane domain (TM). The pore region is indicated with the black line. Information in purple indicates the two regions that interacted with calmodulin (CaM) that have two alpha helixes (α-hA in green and α-hB in dark purple).
Ijms 21 04447 g002
Figure 3. Analysis of unrestricted molecular dynamics simulation of the wild-type and variant calmodulin domain of the KCNQ2 proteins. (A) Structural model of the CaM binding domain of WT-KCNQ2 (helices hA in red, and hB in blue) in the presence of CaM in the calcium-saturated state (holoCaM, orange). Calcium ions are depicted as yellow spheres. (B) Same structure as in (A) after 40 ns of unrestricted molecular dynamics simulation. (C) Structural model of the CaM binding domain of the variant DUP-KCNQ2, containing the duplicated sequence (VQEQHRQKHFEK sequence, residues 320–331, green loop). (D) Same structure as in (C) after 40 ns of unrestricted molecular dynamics simulation. (E) Root mean square deviation (rmsd) values measured over the 40 ns molecular dynamics trajectories of the wild-type (WT-KCNQ2) and the variant (DUP-KCNQ2) structures. Note the higher variation of the rmsd values in the case of the DUP-KCNQ2 structure, caused by the presence of the duplicated sequence, while the WT-KCNQ2 structure does not change significantly.
Figure 3. Analysis of unrestricted molecular dynamics simulation of the wild-type and variant calmodulin domain of the KCNQ2 proteins. (A) Structural model of the CaM binding domain of WT-KCNQ2 (helices hA in red, and hB in blue) in the presence of CaM in the calcium-saturated state (holoCaM, orange). Calcium ions are depicted as yellow spheres. (B) Same structure as in (A) after 40 ns of unrestricted molecular dynamics simulation. (C) Structural model of the CaM binding domain of the variant DUP-KCNQ2, containing the duplicated sequence (VQEQHRQKHFEK sequence, residues 320–331, green loop). (D) Same structure as in (C) after 40 ns of unrestricted molecular dynamics simulation. (E) Root mean square deviation (rmsd) values measured over the 40 ns molecular dynamics trajectories of the wild-type (WT-KCNQ2) and the variant (DUP-KCNQ2) structures. Note the higher variation of the rmsd values in the case of the DUP-KCNQ2 structure, caused by the presence of the duplicated sequence, while the WT-KCNQ2 structure does not change significantly.
Ijms 21 04447 g003

Share and Cite

MDPI and ACS Style

Lazo, P.A.; García, J.L.; Gómez-Puertas, P.; Marcos-Alcalde, Í.; Arjona, C.; Villarroel, A.; González-Sarmiento, R.; Fons, C. Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome. Int. J. Mol. Sci. 2020, 21, 4447. https://doi.org/10.3390/ijms21124447

AMA Style

Lazo PA, García JL, Gómez-Puertas P, Marcos-Alcalde Í, Arjona C, Villarroel A, González-Sarmiento R, Fons C. Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome. International Journal of Molecular Sciences. 2020; 21(12):4447. https://doi.org/10.3390/ijms21124447

Chicago/Turabian Style

Lazo, Pedro A., Juan L. García, Paulino Gómez-Puertas, Íñigo Marcos-Alcalde, Cesar Arjona, Alvaro Villarroel, Rogelio González-Sarmiento, and Carmen Fons. 2020. "Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome" International Journal of Molecular Sciences 21, no. 12: 4447. https://doi.org/10.3390/ijms21124447

APA Style

Lazo, P. A., García, J. L., Gómez-Puertas, P., Marcos-Alcalde, Í., Arjona, C., Villarroel, A., González-Sarmiento, R., & Fons, C. (2020). Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome. International Journal of Molecular Sciences, 21(12), 4447. https://doi.org/10.3390/ijms21124447

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