- 1Laboratory of Brain Modulation and Repair, Glenn Biggs Institute for Alzheimer’s’ and Neurodegenerative Diseases, University of Texas Health San Antonio, San Antonio, TX, United States
- 2Department of Computer Science and Artificial Intelligence, Andalusian Research Institute in Data Science and Computational Intelligence (DaSCI), University of Granada, Granada, Spain
- 3Instituto de Investigación Sanitaria ibs. GRANADA, Hospitales Universitarios de Granada-Universidad de Granada, Granada, Spain
- 4Department of Psychiatry, Washington University, St. Louis, MO, United States
The small-conductance calcium-activated potassium channel SK3, encoded by the KCNN3 gene, plays a critical role in regulating dopaminergic neuron (DN) firing patterns by modulating after hyperpolarization currents. SK3 dysfunction has been implicated in neuropsychiatric and neurodegenerative disorders. We analyzed structural and functional consequences of KCNN3 splicing and genetic variation. Alternative splicing variants of the KCNN3 gene were retrieved from the Ensembl database and aligned using T-Coffee, manually inspected and curated. Protein domains were identified with Pfam 35.0, SMART 9.0, and InterPro 98.0, and visualized. An AlphaFold2 model of SK3 full-length protein (UniProt: Q9UGI6) used as reference and structural models of its splicing variants were predicted with ColabFold. Functional domains (S1–S6 transmembrane helices, H5 pore loop, and calmodulin-binding) were defined and superimposed onto the AlphaFold2 reference. Domain integrity was assessed based on completeness of all expected residue indices within each functional region. SNPs and CNVs across all coding KCNN3 splicing variants were analyzed, classified, and filtered to isolate pathogenic variants prioritizing non-synonymous amino acid substitutions. Differential variant impacts across splicing isoforms were assessed by mapping variant positions to individual transcript protein sequences and used to predict functional consequences. Two long and two short splicing variants are known. Short variants lack the motif required for potassium channels. Pathogenic variants result from missense mutations resulting in amino acid substitutions. In all cases, the consequential effects depend on the specific location and role of the amino acid being changed.
Introduction
During embryonic development in humans, dopaminergic neurons DN are formed in the mesencephalon and differentiate into three major clusters corresponding to the Retrobural Field (RRF), Substantia Nigra pars compacta (SNc) and Ventral Tegmental Area (VTA) respectively (Hegarty et al., 2013). DN in vivo typically fire either regularly, in a pacemaker fashion, or in bursts (Grace and Bunney, 1984a, 1984b; Shepard & German, 1988a). Pacemaker activity produces a tonic release of dopamine (Grace, 1991; Overton and Clark, 1997), while burst firing results in phasic release of dopamine in terminal axonal fields (Bean and Roth, 1991; Gonon and Buda, 1985). When partially denervated in brain slices, DN only maintain pacemaker activity indicating that burst firing is driven by synaptic input (Kita et al., 1986; Sanghera et al., 1984). Bursts consists of 3–7 action potentials that progressively decrease in amplitude and increase in duration, followed by a pause. Non-burst firing patterns, whether pacemaking or not, consist of similarly shaped spikes (Grace and Bunney, 1984a, 1984b; Gu et al., 1991; Shepard & German, 1988b). The controlling switch between regular and burst firing patterns is the small conductance potassium channel sensitive to calcium (SK) (Gu et al., 1991) which is the subject of this work.
There are four subfamilies of SK channels: SK1, SK2, SK3 and the intermediate conductance channel, IKCa1. The latter is encoded by the gene KCNN4 which is not expressed in neurons and therefore will not be discussed further (González et al., 2012). On the other hand. SK1 and SK2 channels are highly expressed in the hippocampus and cortex, whereas SK3 channels are abundant in the thalamus and hypothalamus and have the highest expression levels in the midbrain (Sah and Faber, 2002). All SK channels are voltage independent but open in response to nanomolar changes in cytosolic free calcium (Park, 1994); the ensuing outward potassium current contributes to afterhyperpolarization after a single spike event as well as in the pacemaker properties of neurons (Herrik et al., 2010; Ji et al., 2009; Ji and Shepard, 2006; Wolfart and Roeper, 2002). Following an action potential, potassium currents cause an after-hyperpolarization event (AHP) consisting of fast, medium and slow components (Lancaster and Nicoll, 1987; Sah, 1996; Storm, 1987, 1990). In midbrain DNs, SK3 channels mediate the medium AHP component and control firing frequency (Seutin and Liégeois, 2007; Waroux et al., 2005; Wolfart and Roeper, 2002). Action potential-caused depolarization causes influx of ionized calcium into the cytosol via voltaje dependent calcium channels; further release from intracellular calcium stores is mediated by calcium-induced calcium release receptors (ryanodine receptors) leading to the activation of SK3 channels (Blatz and Magleby, 1986; Park, 1994) and Kd below 1 μmol/L (Halling et al., 2022). SK3 channels do not have a calcium-binding motif; instead, gating is mediated by calmodulin binding to the heteromeric complex of four pore-forming α subunits (Lee et al., 2003). Adult SNpc neurons richly express SK3 mRNAs, whereas DN in the VTA do so much less, and accordingly have smaller AHP currents (Wolfart et al., 2001), but expression is more ubiquitous during embryonic development (Sarpal et al., 2004).
SK channel protein subfamily members in humans are encoded by the genes KCNN1, KCNN2 and KCNN3 (Faber and Sah, 2007; Kshatri et al., 2018). Genetic mutations in KCNN3 are associated with risk of disease (Chandy et al., 1998; Gargus et al., 1998; Stöber et al., 1998) but their impact in channel conductance and neuronal function are poorly understood. Mutations in the KCNN3 gene are associated with risk of schizophrenia and bipolar disorder type 1 (Chandy et al., 1998; Gargus et al., 1998; Stöber et al., 1998), as well as of Huntington’s disease and anorexia nervosa (Dallérac et al., 2015; Koronyo-Hamaoui et al., 2002; Saleem et al., 2000).
The KCNN3 gene contains two regions with CAG trinucleotide repeats encoding polyglutamine segments in the protein (Chandy et al., 1998). Because trinucleotide expansions have known pathological impact, the triplet expansion found in the sequence of the KCNN3 gene is of interest. Somewhat surprisingly, an association was found between reduced number of CAG repeats in KCNN3 with schizophrenia (Chandy et al., 1998; Gargus et al., 1998; Ivković et al., 2006; Stöber et al., 1998) but several replication studies failed to confirm the finding (Glatt et al., 2003; Laurent et al., 2003; Li et al., 1998; Ritsner et al., 2003; Saleem et al., 2000; Tsai et al., 1999; Ujike et al., 2001; Wittekindt et al., 1998). Nevertheless, the number of CAG repeats and the differences in allele size may be related with specific disease features, such as time of onset or cognitive performance (Grube et al., 2011; Mitas, 1997).
Cells transfected with transgenic KCNN3 constructs with increasing number of CAG showed a significant reduction in overall conductance and stronger inward rectification (Grube et al., 2011). A 4-base pair deletion in this gene was identified in a patient with schizophrenia and this mutation generated a truncated protein (hSK3Δ) without channel function (Bowen et al., 2001; Miller et al., 2001). In DN, expression of hSK3Δ also suppressed the SK endogenous current and caused increased burst firing activity and increased intracellular calcium signaling (Soden et al., 2013). In this article we provide new data on the impact of known human genetic variations in KCNN3 on its probably channel properties and associated conductance, and discuss possible impact of SK3 channels modifications in DN function.
Materials and methods
Data acquisition and sequence alignments
Sequence retrieval and alignment
Alternative splicing variants of the potassium calcium-activated channel subfamily N member 3 (KCNN3) gene were retrieved from the Ensembl database (release 110) (Martin et al., 2023). Multiple sequence alignment of the variant transcripts was performed using T-Coffee version 11.00.8cbe486 with default parameters. The resulting alignments were manually inspected and curated to ensure accuracy, with particular attention to splice junction regions and protein-coding segments (Figure 1).
Figure 1. Comparative analysis of KCNN3 splicing variants. (A) Sequence alignment of four KCNN3 transcripts reveal two distinct length categories. Long variants (ENST00000271915.9 and ENST00000618040.4) differ only by a 16-amino acid insertion at position 484, while short variants (ENST00000358505.2 and ENST00000361147.8) exhibit a 5-amino acids difference in the region spanning positions 309–314. (B) Schematic representation of the coding regions for the main functional domains of the SK3 channel protein.
Functional domain analysis
Protein domain identification was performed using three complementary databases: Pfam 35.0, SMART 9.0, and InterPro 98.0. Results were cross-referenced to identify conserved and variant-specific domains. Domain architectures were visualized using the Ensembl genome browser protein schematic representation tool (Figure 2). Differences in domain content, organization, and integrity between splicing variants were systematically cataloged.
Figure 2. Structural characterization of KCNN3 protein isoforms. (A) Domain architecture of KCNN3 (743 amino acids) and its splicing variants. Functional domains and conserved motifs are mapped to their respective positions, with domain boundaries indicated by amino acid position. High-confidence domains (Pfam/SMART identified) are shown in darker colors, while lower-confidence predictions appear in lighter shades. This representation was generated using the Ensembl Protein Feature View tool, integrating data from multiple prediction algorithms and experimental validation. Only domain predictions with significant scores (E-value < 0.001) are included. Annotated protein structure of the human KCNN3 gene product (SK3 channel), highlighting key functional domains across its isoforms. Conserved regions include six transmembrane domains (S1–S6), a pore-forming loop between S5 and S6, and a C-terminal calmodulin-binding domain. Protein domains were mapped using Pfam, PRINTS, PROSITE, and SUPERFAMILY databases. Functional annotations include: PF07885: Two pore domain potassium channel domain; PF03530: Potassium channel, calcium-activated, SK subfamily; PF02888: Calmodulin-binding domain; PR01451/PS00322/PS00909: Signatures of calcium-activated potassium channels. These structural features are critical for the calcium sensitivity and electrophysiological function of SK3 channels. Variability across transcript isoforms may influence channel gating, localization, and regulatory interactions. (B–E) Three-dimensional protein structures of KCNN3 splicing variants predicted using the AlphaFold 3 (Abramson et al., 2024). (B) Functional variant ENST00000618040.4, (C) functional variant ENST00000271915.9, (D,E) shorter non-functional variants ENST00000361147.8 and ENST00000358505.2, respectively, revealing the absence of key transmembrane helices in the truncated isoforms (Martin et al., 2023). Renderings were generated in the Human Protein Atlas (https://www.proteinatlas.org/ENSG00000143603-KCNN3/structure+interaction).
Structural model and analysis
Reference structure acquisition
The high-confidence AlphaFold2 (Jumper et al., 2021) model for human SK3 full-length protein (AF-Q9UGI6-F1) was obtained from the AlphaFold Protein Structure Database (Varadi et al., 2024). The canonical KCNN3 sequence (UniProt: Q9UGI6) (Ahmad et al., 2025) served as the reference for comparative analysis. Two-dimensional (Figures 2A, 3) and three-dimensional (Figures 2B–E) structural models of SK3 splicing variants (Variant 202, Variant 203, Variant 205) were predicted using ColabFold (Mirdita et al., 2022), a high-efficiency AlphaFold2 (Jumper et al., 2021) implementation. Structural predictions were performed in batch mode using the following parameters: --model-type alphafold2_ptm, --num-models 5, --num-recycles 3, --rank-by confidence, --use-templates (Abramson et al., 2024), with the wild-type AlphaFold2 model supplied as a template. The top-ranked model based on predicted local distance difference test (pLDDT) scores was selected for each variant.
Figure 3. Structural and functional annotation of KCNN3 transcript isoforms. Schematic representation of four transcript isoforms of the KCNN3 gene (KCNN3-201, -202, -203, and -205), including exon-intron structures and corresponding protein features. Top panel: Genomic organization and alternative splicing patterns of the four KCNN3 transcripts. Black boxes represent coding exons; gray lines represent introns. Panels 201–205: Isoform-specific annotations including Antibody target regions (green bars) for antibodies HPA017990 and HPA057127. Sequence alignment reliability from Human Protein Atlas (HsID 10 and HsID 50), where the red-yellow line indicates alignment confidence across the amino acid sequence; Predicted functional and structural features, including: MDM domains (brown); Low complexity regions and intrinsically disordered regions (green and gray); Common and isoform-specific immunogenic regions (purple and orange). Predicted secondary structure elements (alpha helices in red, beta strands in blue, coils in green). Together, these annotations highlight key structural and functional differences among KCNN3 isoforms, which may influence protein stability, localization, and interaction profiles. Data adapted from the Human Protein Atlas (https://www.proteinatlas.org/ENSG00000143603-KCNN3/structure+interaction).
Structure processing and domain annotation
PDB models were parsed using Biopython’s PDBParser (Cock et al., 2009) to verify AlphaFold confidence scores in the B-factor field. Residues with confidence scores below 70 were flagged for cautious interpretation. Functional domains were defined based on literature-curated regions including S1–S6 transmembrane helices, the pore (H5) loop, and the C-terminal calmodulin-binding domain (Supplementary Table 1) (Ahmad et al., 2025).
Structural comparison analysis
Cα-based root-mean-square deviation (RMSD) values were calculated using Biopython’s Superimposer for both global (residues 1-731) and domain-specific alignments. Each variant was superimposed onto the AlphaFold2 reference using common residues within each domain. Domain RMSD calculations required ≥5 aligned Cα pairs; otherwise, values were assigned as not available. Per-residue confidence scores (pLDDT) were extracted and mapped to assess regional reliability, with mean pLDDT values reported per domain. Domain integrity was assessed based on completeness of all expected residue indices within each functional region.
Genetic variant analysis
Variant data collection
Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) across all coding KCNN3 gene splicing variants were analyzed using the Ensembl Variant Effect Predictor (VEP) tool (release 110) (Martin et al., 2023). Variant data were retrieved from gnomAD v3.1.2 (Chen et al., 2024) and ClinVar (Landrum et al., 2014) databases. VEP analysis was performed with parameters: “--species homo_sapiens --assembly GRCh38 --cache --offline --symbol --af --af_1kg --af_gnomad --pubmed --domains --regulatory --biotype --check_existing.” All KCNN3 transcript identifiers (ENST00000618040.4, ENST00000271915.9, ENST00000361147.8, ENST00000358505.2) were specified for targeted analysis.
Pathogenicity assessment
Variants were classified according to ACMG/AMP guidelines (Martin et al., 2023; Richards et al., 2015) and filtered to isolate pathogenic and likely pathogenic variants. Non-synonymous amino acid substitutions were prioritized and categorized into four primary types: (1) Alanine to Threonine, (2) Serine to Cysteine, (3) Glycine to Aspartic acid, and (4) Lysine to Glutamic acid conversions. Variant positions were analyzed relative to functional domains identified through domain analysis.
Transcript-specific impact analysis
Differential variant impacts across splicing variants were assessed by mapping variant positions to individual transcript protein sequences. Analysis focused on variants affecting conserved functional domains versus transcript-specific regions, with particular attention to differences between long and short SK3 isoforms. Functional consequences were predicted based on amino acid changes, structural location, and cross-species conservation patterns.
Results
Overview of variant structures and domains
We analyzed splicing variants and known structural mutations of the gene to predict their functional impact. Currently, there are five known splicing variants of the KCNN3 gene, only four of which are in coding regions. Variants ENST00000618040.4 and ENST00000271915.9 encompass 746 and 731 amino acids, respectively. The other two coding variants, ENST00000361147.8 and ENST00000358505.2 respectively, are much shorter encompassing 426 and 418 amino acids (Figure 1). Variants ENST00000618040.4 and ENST00000271915.9 have a nearly identical sequence, the only difference being a segment insertion in position 484. Variants ENST00000361147.8 and ENST00000358505.2 are also nearly identical, except for positions 309–314. Three dimensional structures obtained using AlphaFold reveal visible differences in the protein conformation as the result of each variant (Figures 2B–E).
Global structural divergence from the AlphaFold reference model
The global Cα RMSD values revealed a step-wise divergence among the splice isoforms (Table 1). Variant 202 and Variant 203 produced global RMSD values > 48 Å, whereas Variant_205 deviated by 17.6 Å. The wild type variant (Variant 201) aligned almost perfectly with the AlphaFold model (RMSD ~0 Å), validating the reference structure. We believe that the extreme RMSD of Variants 202/203 is driven primarily by sequence truncations rather than wholesale mis-folding; absent segments inflate RMSD because no equivalent C-α positions exist for superposition.
Analysis of the short variants and their coding regions with more detail using Ensembl data base (Martin et al., 2023) revealed absence of the motif required for potassium channels (PF03530 Potassium channel, calcium-activated SK) (Figures 2A, 3). In contrast, the two longer variants express all the segments for that motif and are both functional despite the segment insertion. This led us to do a more detailed investigation at the domain structure level.
Domain-level RMSD analysis and structural integrity
When dissected by domain, all transmembrane (TM) regions (S1–S6) and the pore loop (H5) were preserved in variant 201 and AlphaFold model, as expected (RMSD < 3e−14 Å, all domains present). Variant 205 showed very low RMSD values across all TM domains (≤0.38 Å), indicating a structurally intact pore-forming unit. In contrast, Variant 202 and Variant 203 retained the S1–S3 regions with minimal divergence (RMSD ~0.9–5.5 Å) but exhibited markedly higher RMSD values and missing data across S4, S5, H5, and S6 domains. Specifically; Variant_202 showed complete loss of domains from S4 onwards. Variant 203 retained S4 (RMSD = 9.00 Å) but was missing S5–S6 and the pore loop (Table 1).
Functional consequences of domain loss across SK3 variants
To contextualize the structural observations, we assessed the functional implications of domain-level disruptions observed in each SK3 isoform. Transmembrane (TM) integrity and the calmodulin-binding domain were evaluated for presence and AlphaFold-predicted confidence (Supplementary Table S1). Wild-Type and AlphaFold models presented all gating and pore motifs intact indicating that likely expected to conduct K+ and respond to Ca2+. In Variant 205, all TM and regulatory domains were also present, with minor structural deviations from AlphaFold but preserved architecture. Thus, this variant likely retains most functional properties of the wild-type channel, including potassium conductance and calcium sensitivity.
Variant 203 presents a deletion of S5, pore (H5), and S6 domains, which abolished the K + permeation pathway; the CaM-BD region loss removes the Ca2 sensor suggesting a disrupted gating and compromised ion selectivity, consistent with severe functional impairment. Lastly, Variant_202 exhibited an even greater truncation (S4-S6 absent) which predicts a complete loss-of-function.
Structural variants were studied in Ensembl in an effort to study the documented mutations that have a known pathogenic effect and further analyzed with SIFT, Polyphen, CADD, REVEL y MetaLR (Martin et al., 2023), eventually settling on SIFT. From the 2,394 reported variants (Supplementary Table S2) 13 variants were identified that were predicted as pathogenic with SIFT (Table 2) and were backed by supportive literature. As shown in Table 2, the 13 pathogenic variants result from missense mutations, in which the resulting protein is largely the same but will encode a different amino acid at a given position due to a nucleotide substitution. In all cases, the consequential effects depend on the specific location and role of the amino acid being changed.
Discussion
KCNN3 genomic variants have been associated with brain disease and SK channels are expressed widely in brain and play a key role in the control of neuronal activity; thus, drugs modulating SK3 activity may be useful to treat neuropsychiatric disorders (Blank et al., 2004).
SK3 conductance is specifically blocked by apamin and activated by 1-ethyl-1,3-dihydro-2H-benzimidazol-2-one (1-EBIO) and NS309. Another activator, cyclohexyl-[2-(3,5-dimethyl-pyrazol-1-yl)-6-methyl-pyrimidin-4-yl]-amine (CyPPA) is an SK2- and SK3-specific positive modulator (Hougaard et al., 2007). Generally, positive SK3 modulators decrease neuronal activity whereas negative modulation increases firing rate and bursting in neurons (Herrik et al., 2010). Blocking SK3 channels with apamin changes DN pacemaking into a burst firing pattern (Gu et al., 1991; Waroux et al., 2005).
In cell cultures of DN, the DN-specific death pathway produced by excitotoxicity (Figure 4), is prevented by SK3 activation with 1-EBIO or CyPPA and facilitated by the SK3 antagonist apamin (Benítez et al., 2011). Likewise, SK channel activation by NS309 confers neuroprotection against oxidative stress with rotenone (Dolga et al., 2014), and 1-EBIO protects DN in organotypic cultures from the neurotoxic effect of 6-hydroxydopamine (6-OHDA) (Wang et al., 2015). In intact animals, 6-OHDA reduces expression of SK3 channels in the substantia nigra pars compacta of the lesioned animals (Mourre et al., 2017).
Mutations in KCNN3 have been associated with schizophrenia (Chandy et al., 1998; Gargus et al., 1998; Stöber et al., 1998). KCNN3 contains two regions with CAG repeats (Chandy et al., 1998) with expansions reported in association with schizophrenia (Mitas, 1997), and with earlier onset of the disease (Ritsner et al., 2003). Other studies found no relation between the CAG trinucleotide repeat and schizophrenia disease or even the onset age (Glatt et al., 2003; Laurent et al., 2003; Li et al., 1998; Saleem et al., 2000; Tsai et al., 1999; Ujike et al., 2001; Wittekindt et al., 1998). A truncated variant of the SK3 with 283 amino acids instead of usual 731 (Bowen et al., 2001) results in a loss of function SK3 channel enhancing the excitability of DN in vivo (Soden et al., 2013).
Regarding the gene structure of KCNN3, splicing variants can result in lack of the potassium channel or the calmodulin-binding domain. A truncated SK3 variant similar to a schizophrenia-related mutation did not produce functional channels but selectively suppressed endogenous SK3 currents in a dominant-negative fashion (Tomita et al., 2003). Although the splicing variants ENST00000361147.8 and ENST00000358505.2 are not identical to this truncated variant, the results presented here strongly suggest a similar or identical effect.
The current data consolidate a dual structural/functional role for the Calmodulin-binding domain. Its deletion in Variants 202/203 coincides with wholesale loss of the distal TM bundle (S4-S6) and the pore helix, producing the highest global RMSD (> 48 Å) and predicting non-conductive, Ca2+-insensitive channels. Conversely, its retention in WT and Variant 205 correlates with an intact pore and minimal RMSD (<18 Å). These observations extend previous biochemical findings that Calmodulin binding both senses intracellular Ca2+ and stabilizes channel architecture. While structurally coherent, the conclusions remain predictive. Post-translational modifications, membrane lipid composition, and protein dynamics are not captured by static in-silico models. Rigorous functional assays are therefore essential before any assertion of clinical impact, particularly regarding neurodegeneration or psychiatric risk, can be made.
We studied in detail 13 documented variants reported in Zimmermann-Laband Syndrome (Bauer et al., 2019; Nam et al., 2023; Schwarz et al., 2022). In further review of the literature, no documented studies were found that focused solely on the possible link of the structural variants of KCNN3 with brain diseases. It is important to note, however, that mutations in KCNN3 that have been related with Zimmermann-Laband are a gain-of-function mutation or mutations that involve a loss of function of the channel have been linked with neurodevelopmental disorders (Nam et al., 2023). Figure 4 provides a model for the role of SK3 channels in dopaminergic neuronal death which may be used as a framework to understand the pathogenicity of channel variants (Maldonado et al., 2020; Benítez et al., 2011; de Erausquin et al., 2003; Dorsey et al., 2006). Further investigation is needed to elucidate the possible effects of the splicing variants in disease is likely to yield new therapeutic targets of possible relevance for Parkinson’s disease, schizophrenia, or both.
Figure 4. Proposed role of SK3 channels in excitotoxicity-mediated dopaminergic neuronal death: A toxic amount of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) contacts the GluRAMPA receptor in dopaminergic neurons (DNs), (2) allowing sodium (Na+) into the DNs, causing depolarization in the membrane. (3) This depolarization activates voltage-dependent L-type calcium channels (LTCC), triggering an influx of calcium (Ca2+) into the DNs, (4) which further releases intracellular Ca2+ stores through the IP3 pathway and ryanodine receptors (RyR). (5) At this point in a normal cell an afterhyperpolarization current (IAHP) would occur, with potassium (K+) exiting the cell through SK3. However, experimental data shows that under AMPA excitotoxicity, the IAHP is significantly reduced or does not take place at all. Meaning that, while the mechanism of this inhibition is not fully understood, it is safe to assume that calmodulin, the calcium-binding protein that activates SK3, is unable to bind normally during AMPA excitotoxicity, which is hypothesized in the above figure. Also, due to their unique sensitivity to AMPA-induced excitotoxicity, DNs are unable to regulate the intracellular levels of Ca2+. This creates a (6) destabilization of neuronal ionized Ca2+ ([Ca2+]i) homeostasis, (7) leading to oxidative stress, mitochondrial swelling, increased levels of reactive oxygen species (ROS) and the release of cytochrome C (cyt c). (8) The commitment to die is determined by the phosphatidylserine (PS) translocation to the cellular membrane, which is determined by its attachment to annexin V. (9) After phosphorylation, inhibitor of nuclear factor κ B (IκB) is separated from nuclear factor κ B (NFκB) and then degraded. NFκB then translocates to the nucleus, (10) where it activates transcription of protein p53, which later phosphorylates into phos-p53. (11) This is the last step in the dopaminergic neuron-specific programmed death pathway.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://www.google.com/url?sa=t&source=web&rct=j&opi=89978449&url=https://www.ensembl.org/index.html&ved=2ahUKEwi7iK6HqrCNAxUWLLkGHWNnGywQFnoECAkQAQ&usg=AOvVaw3m00q2JEfEz-KntGt1NVcQ.
Author contributions
LP: Data curation, Investigation, Visualization, Writing – original draft, Writing – review & editing. CD: Data curation, Formal analysis, Methodology, Resources, Software, Visualization, Writing – original draft, Writing – review & editing. DN: Formal analysis, Project administration, Supervision, Writing – original draft, Writing – review & editing. AK: Data curation, Visualization, Writing – original draft, Writing – review & editing. JM: Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. AT: Conceptualization, Supervision, Writing – original draft, Writing – review & editing. IZ: Conceptualization, Formal analysis, Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing. GE: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported in part by NIH grants U19AG076581 and P30AG066546 to GE and by the Spanish Ministry of Science and Technology project RTI2018-098983-B-100 and “Ethical, Responsible, and General-Purpose Artificial Intelligence: Applications in Risk Scenarios (IAFER)” Exp.: TSI-100927-2023-1.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2025.1631536/full#supplementary-material
References
Abramson, J., Adler, J., Dunger, J., Evans, R., Green, T., Pritzel, A., et al. (2024). Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature 630, 493–500. doi: 10.1038/s41586-024-07487-w
Ahmad, S., Jose da Costa Gonzales, L., Bowler-Barnett, E. H., Rice, D. L., Kim, M., Wijerathne, S., et al. (2025). The UniProt website API: facilitating programmatic access to protein knowledge. Nucleic Acids Res. 53, W547–W553. doi: 10.1093/NAR/GKAF394
Bauer, C. K., Schneeberger, P. E., Kortüm, F., Altmüller, J., Santos-Simarro, F., Baker, L., et al. (2019). Gain-of-function mutations in KCNN3 encoding the small-conductance Ca2+-activated K+ channel SK3 cause Zimmermann-Laband syndrome. Am. J. Hum. Genet. 104, 1139–1157. doi: 10.1016/J.AJHG.2019.04.012
Bean, A. J., and Roth, R. H. (1991). Extracellular dopamine and neurotensin in rat prefrontal cortex in vivo: effects of median forebrain bundle stimulation frequency, stimulation pattern, and dopamine autoreceptors. J. Neurosci. 11, 2694–2702. doi: 10.1523/JNEUROSCI.11-09-02694.1991
Benítez, B. A., Belálcazar, H. M., Anastasía, A., Mamah, D. T., Zorumski, C. F., Mascó, D. H., et al. (2011). Functional reduction of SK3-mediated currents precedes AMPA-receptor-mediated excitotoxicity in dopaminergic neurons. Neuropharmacology 60, 1176–1186. doi: 10.1016/J.NEUROPHARM.2010.10.024
Blank, T., Nijholt, I., Kye, M. J., and Spiess, J. (2004). Small conductance Ca2+−activated K+ channels as targets of CNS drug development. Curr. Drug Targets CNS Neurol. Disord. 3, 161–167. doi: 10.2174/1568007043337472
Blatz, A. L., and Magleby, K. L. (1986). Single apamin-blocked ca-activated K+ channels of small conductance in cultured rat skeletal muscle. Nature 323, 718–720. doi: 10.1038/323718A0
Bowen, T., Williams, N., Norton, N., Spurlock, G., Wittekindt, O. H., Morris-Rosendahl, D. J., et al. (2001). Mutation screening of the KCNN3 gene reveals a rare frameshift mutation. Mol. Psychiatry 6, 259–260. doi: 10.1038/SJ.MP.4000128
Chandy, K. G., Fantino, E., Wittekindt, O., Kalman, K., Tong, L. L., Ho, T. H., et al. (1998). Isolation of a novel potassium channel gene hSKCa3 containing a polymorphic CAG repeat: a candidate for schizophrenia and bipolar disorder? Mol. Psychiatry 3, 32–37. doi: 10.1038/SJ.MP.4000353
Chen, S., Francioli, L. C., Goodrich, J. K., Collins, R. L., Kanai, M., Wang, Q., et al. (2024). A genomic mutational constraint map using variation in 76,156 human genomes. Nature 625, 92–100. doi: 10.1038/S41586-023-06045-0
Cock, P. J. A., Antao, T., Chang, J. T., Chapman, B. A., Cox, C. J., Dalke, A., et al. (2009). Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25, 1422–1423. doi: 10.1093/BIOINFORMATICS/BTP163
de Erausquin, G. A., Hyrc, K., Dorsey, D. A., Mamah, D., Dokucu, M., Mascó, D. H., et al. (2003). Nuclear Translocation of Nuclear Transcription Factor-κB by α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic Acid Receptors Leads to Transcription of p53 and Cell Death in Dopaminergic Neurons. Molecular Pharmacology, 63, 784–790. doi: 10.1124/mol.63.4.784
Dallérac, G. M., Levasseur, G., Vatsavayai, S. C., Milnerwood, A. J., Cummings, D. M., Kraev, I., et al. (2015). Dysfunctional dopaminergic neurones in mouse models of Huntington’s disease: a role for SK3 channels. Neurodegener. Dis. 15, 93–108. doi: 10.1159/000375126
Dolga, A. M., De Andrade, A., Meissner, L., Knaus, H. G., Höllerhage, M., Christophersen, P., et al. (2014). Subcellular expression and neuroprotective effects of SK channels in human dopaminergic neurons. Cell Death Dis 5, e999–e911. doi: 10.1038/cddis.2013.530
Dorsey, D. A., Mascó, D. H., Dikranian, K., Hyrc, L., Masciotra, B., Faddis, M., et al. (2006). Ultrastructural characterization of alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid-induced cell death in embryonic dopaminergic neurons. Apoptosis 11, 535–544. doi: 10.1007/s10495-006-5268-y
Faber, E. S. L., and Sah, P. (2007). Functions of SK channels in central neurons. Clin. Exp. Pharmacol. Physiol. 34, 1077–1083. doi: 10.1111/J.1440-1681.2007.04725.X
Gargus, J. J., Fantino, E., and Gutman, G. A. (1998). A piece in the puzzle: an ion channel candidate gene for schizophrenia. Mol. Med. Today 4, 518–524. doi: 10.1016/S1357-4310(98)01358-6
Glatt, S. J., Faraone, S. V., and Tsuang, M. T. (2003). CAG-repeat length in exon 1 of KCNN3 does not influence risk for schizophrenia or bipolar disorder: a meta-analysis of association studies. Am J Med Genet Part B Neuropsychiatr Genet. 121B, 14–20. doi: 10.1002/AJMG.B.20048
Gonon, F. G., and Buda, M. J. (1985). Regulation of dopamine release by impulse flow and by autoreceptors as studied by in vivo voltammetry in the rat striatum. Neuroscience 14, 765–774. doi: 10.1016/0306-4522(85)90141-1
González, C., Baez-Nieto, D., Valencia, I., Oyarzún, I., Rojas, P., Naranjo, D., et al. (2012). K+ channels: function-structural overview. Compr. Physiol. 2, 2087–2149. doi: 10.1002/j.2040-4603.2012.tb00452.x
Grace, A. A. (1991). Phasic versus tonic dopamine release and the modulation of dopamine system responsivity: a hypothesis for the etiology of schizophrenia. Neuroscience 41, 1–24. doi: 10.1016/0306-4522(91)90196-U
Grace, A. A., and Bunney, B. S. (1984a). The control of firing pattern in nigral dopamine neurons: burst firing. J. Neurosci. 4, 2877–2890. doi: 10.1523/JNEUROSCI.04-11-02877.1984
Grace, A. A., and Bunney, B. S. (1984b). The control of firing pattern in nigral dopamine neurons: single spike firing. J. Neurosci. 4, 2866–2876. doi: 10.1523/jneurosci.04-11-02866.1984
Grube, S., Gerchen, M. F., Adamcio, B., Pardo, L. A., Martin, S., Malzahn, D., et al. (2011). A CAG repeat polymorphism of KCNN3 predicts SK3 channel function and cognitive performance in schizophrenia. EMBO Mol. Med. 3, 309–319. doi: 10.1002/emmm.201100135
Gu, X., Blatzt, A. L., and German, D. C. (1991). Subtypes of substantia nigra dopaminergic neurons revealed by apamin: autoradiogaphic and electrophysiological studies. Brain Res. Bull. 28, 435–440.
Halling, D. B., Philpo, A. E., and Aldrich, R. W. (2022). Calcium dependence of both lobes of calmodulin is involved in binding to a cytoplasmic domain of SK channels. eLife 11:303. doi: 10.7554/ELIFE.81303
Hegarty, S. V., Sullivan, A. M., and O’Keeffe, G. W. (2013). Midbrain dopaminergic neurons: a review of the molecular circuitry that regulates their development. Dev. Biol. 379, 123–138. doi: 10.1016/J.YDBIO.2013.04.014
Herrik, K. F., Christophersen, P., and Shepard, P. D. (2010). Pharmacological modulation of the gating properties of small conductance Ca2+−activated K+ channels alters the firing pattern of dopamine neurons in vivo. J. Neurophysiol. 104, 1726–1735. doi: 10.1152/JN.01126.2009
Hougaard, C., Eriksen, B. L., Jørgensen, S., Johansen, T. H., Dyhring, T., Madsen, L. S., et al. (2007). Selective positive modulation of the SK3 and SK2 subtypes of small conductance ca 2+−activated K + channels. Br. J. Pharmacol. 151, 655–665. doi: 10.1038/sj.bjp.0707281
Ivković, M., Ranković, V., Tarasjev, A., Orolicki, S., Damjanović, A., Paunović, V. R., et al. (2006). Schizophrenia and polymorphic CAG repeats array of calcium-activated potassium channel (KCNN3) gene in Serbian population. Int. J. Neurosci. 116, 157–164. doi: 10.1080/00207450341514
Ji, H., Hougaard, C., Herrik, K. F., Strøbæk, D., Christophersen, P., and Shepard, P. D. (2009). Tuning the excitability of midbrain dopamine neurons by modulating the Ca2+ sensitivity of SK channels. Eur. J. Neurosci. 29, 1883–1895. doi: 10.1111/j.1460-9568.2009.06735.x
Ji, H., and Shepard, P. D. (2006). SK Ca2+−activated K+ channel ligands alter the firing pattern of dopamine-containing neurons in vivo. Neuroscience 140, 623–633. doi: 10.1016/J.NEUROSCIENCE.2006.02.020
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589. doi: 10.1038/s41586-021-03819-2
Kita, T., Kita, H., and Kitai, S. T. (1986). Electrical membrane properties of rat substantia nigra compacta neurons in an in vitro slice preparation. Brain Res. 372, 21–30. doi: 10.1016/0006-8993(86)91454-X
Koronyo-Hamaoui, M., Danziger, Y., Frisch, A., Stein, D., Leor, S., Laufer, N., et al. (2002). Association between anorexia nervosa and the hsKCa3 gene: a family-based and case control study. Mol. Psychiatry 7, 82–85. doi: 10.1038/SJ.MP.4000931
Kshatri, A. S., Gonzalez-Hernandez, A., and Giraldez, T. (2018). Physiological roles and therapeutic potential of Ca2+ activated potassium channels in the nervous system. Front. Mol. Neurosci. 11:258. doi: 10.3389/FNMOL.2018.00258
Lancaster, B., and Nicoll, R. A. (1987). Properties of two calcium-activated hyperpolarizations in rat hippocampal neurones. J. Physiol. 389, 187–203. doi: 10.1113/JPHYSIOL.1987.SP016653
Landrum, M. J., Lee, J. M., Riley, G. R., Jang, W., Rubinstein, W. S., Church, D. M., et al. (2014). ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 42, D980–D985. doi: 10.1093/NAR/GKT1113
Laurent, C., Niehaus, D., Bauché, S., Levinson, D. F., Soubigou, S., Pimstone, S., et al. (2003). CAG repeat polymorphisms in KCNN3 (HSKCa3) and PPP2R2B show no association or linkage to schizophrenia. Am. J. Med. Genet. B Neuropsychiatr. Genet. 116B, 45–50. doi: 10.1002/AJMG.B.10797
Lee, W. S., Ngo-Anh, T. J., Bruening-Wright, A., Maylie, J., and Adelman, J. P. (2003). Small conductance Ca2+−activated K+ channels and calmodulin. Cell surface expression and gating. J. Biol. Chem. 278, 25940–25946. doi: 10.1074/jbc.M302091200
Li, T., Hu, X., Chandy, K. G., Fantino, E., Kalman, K., Gutman, G., et al. (1998). Transmission disequilibrium analysis of a triplet repeat within the hKCa3 gene using family trios with schizophrenia. Biochem. Biophys. Res. Commun. 251, 662–665. doi: 10.1006/bbrc.1998.9484
Maldonado, O., Jenkins, A., Belalcazar, H. M., Hernandez-Cuervo, H., Hyman, K. M., Ladaga, G., et al. (2020). Age-dependent neuroprotective effect of an SK3 channel agonist on excitotoxity to dopaminergic neurons in organotypic culture. PLoS One 15:e0223633. doi: 10.1371/journal.pone.0223633
Martin, F. J., Amode, M. R., Aneja, A., Austine-Orimoloye, O., Azov, A. G., Barnes, I., et al. (2023). Ensembl 2023. Nucleic Acids Res. 51, D933–D941. doi: 10.1093/NAR/GKAC958
Miller, M. J., Rauer, H., Tomita, H., Rauer, H., Gargus, J. J., Gutman, G. A., et al. (2001). Nuclear localization and dominant-negative suppression by a mutant SKCa3 N-terminal channel fragment identified in a patient with schizophrenia. J. Biol. Chem. 276, 27753–27756. doi: 10.1074/JBC.C100221200
Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S., and Steinegger, M. (2022). Colabfold: making protein folding accessible to all. Nat. Methods 19, 679–682. doi: 10.1038/S41592-022-01488-1
Mitas, M. (1997). Trinucleotide repeats associated with human disease. Nucleic Acids Res. 25, 2245–2253. doi: 10.1093/NAR/25.12.2245
Mourre, C., Manrique, C., Camon, J., Aidi-Knani, S., Deltheil, T., Turle-Lorenzo, N., et al. (2017). Changes in SK channel expression in the basal ganglia after partial nigrostriatal dopamine lesions in rats: functional consequences. Neuropharmacology 113, 519–532. doi: 10.1016/j.neuropharm.2016.11.003
Nam, Y. W., Downey, M., Rahman, M. A., Cui, M., and Zhang, M. (2023). Channelopathy of small- and intermediate-conductance Ca2+−activated K+ channels. Acta Pharmacol. Sin. 44, 259–267. doi: 10.1038/S41401-022-00935-1
Overton, P. G., and Clark, D. (1997). Burst firing in midbrain dopaminergic neurons. Brain Res. Rev. 25, 312–334. doi: 10.1016/S0165-0173(97)00039-8
Park, Y. B. (1994). Ion selectivity and gating of small conductance ca(2+)-activated K+ channels in cultured rat adrenal chromaffin cells. J. Physiol. 481, 555–570. doi: 10.1113/JPHYSIOL.1994.SP020463
Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 17, 405–424. doi: 10.1038/GIM.2015.30
Ritsner, M., Amir, S., Koronyo-Hamaoui, M., Gak, E., Ziv, H., Halperin, T., et al. (2003). Association study of CAG repeats in the KCNN3 gene in Israeli patients with major psychosis. Psychiatr. Genet. 13, 143–150. doi: 10.1097/00041444-200309000-00002
Sah, P. (1996). Ca(2+)-activated K+ currents in neurones: types, physiological roles and modulation. Trends Neurosci. 19, 150–154. doi: 10.1016/S0166-2236(96)80026-9
Sah, P., and Faber, E. S. L. (2002). Channels underlying neuronal calcium-activated potassium currents. Prog. Neurobiol. 66, 345–353. doi: 10.1016/S0301-0082(02)00004-7
Saleem, Q., Sreevidya, V., Sudhir, J., Vijaya Savithri, J., Gowda, Y., B-Rao, C., et al. (2000). Association analysis of CAG repeats at the KCNN3 locus in Indian patients with bipolar disorder and schizophrenia. Am J. Med. Genet 96, 744–748. doi: 10.1002/1096-8628(20001204)96:6<744::aid-ajmg9>3.0.co;2-z
Sanghera, M. K., Trulson, M. E., and German, D. C. (1984). Electrophysiological properties of mouse dopamine neurons: in vivo and in vitro studies. Neuroscience 12, 793–801. doi: 10.1016/0306-4522(84)90171-4
Sarpal, D., Koenig, J. I., Adelman, J. P., Brady, D., Prendeville, L. C., and Shepard, P. D. (2004). Regional distribution of SK3 mRNA-containing neurons in the adult and adolescent rat ventral midbrain and their relationship to dopamine-containing cells. Synapse 53, 104–113. doi: 10.1002/syn.20042
Schwarz, M., Ryba, L., Křepelová, A., Moslerová, V., Zelinová, M., Turnovec, M., et al. (2022). Zimmermann-Laband syndrome in monozygotic twins with a mild neurobehavioral phenotype lacking gingival overgrowth-a case report of a novel KCNN3 gene variant. Am. J. Med. Genet. A 188, 1083–1087. doi: 10.1002/AJMG.A.62616
Seutin, V., and Liégeois, J. F. (2007). SK channels are on the move. Br. J. Pharmacol. 151, 568–570. doi: 10.1038/sj.bjp.0707282
Shepard, P. D., and German, D. C. (1988). Electrophysiological and pharmacological evidence for the existence of distinct subpopulations of nigrostriatal dopaminergic neuron in the rat. Neuroscience 27, 537–546. doi: 10.1016/0306-4522(88)90287-4
Soden, M. E., Jones, G. L., Sanford, C. A., Chung, A. S., Güler, A. D., Chavkin, C., et al. (2013). Disruption of dopamine neuron activity pattern regulation through selective expression of a human KCNN3 mutation. Neuron 80, 997–1009. doi: 10.1016/j.neuron.2013.07.044
Stöber, G., Jatzke, S., Meyer, J., Okladnova, O., Knapp, M., Beckmann, H., et al. (1998). Short CAG repeats within the hSKCa3 gene associated with schizophrenia: results of a family-based study. Neuroreport 9, 3595–3599. doi: 10.1097/00001756-199811160-00010
Storm, J. F. (1987). Action potential repolarization and a fast after-hyperpolarization in rat hippocampal pyramidal cells. J. Physiol. 385, 733–759. doi: 10.1113/jphysiol.1987.sp016517
Storm, J. F. (1990). Potassium currents in hippocampal pyramidal cells. Progr Brain Res 83, 161–187. doi: 10.1016/S0079-6123(08)61248-0
Tomita, H., Shakkottai, V. G., Gutman, G. A., Sun, G., Bunney, W. E., Cahalan, M. D., et al. (2003). Novel truncated isoform of SK3 potassium channel is a potent dominant-negative regulator of SK currents: implications in schizophrenia. Mol. Psychiatry 8, 524–535. doi: 10.1038/sj.mp.4001271
Tsai, M. T., Shaw, C. K., Hsiao, K. J., and Chen, C. H. (1999). Genetic association study of a polymorphic CAG repeats array of calcium-activated potassium channel (KCNN3) gene and schizophrenia among the Chinese population from Taiwan. Mol. Psychiatry 4, 271–273. doi: 10.1038/SJ.MP.4000482
Ujike, H., Yamamoto, A., Tanaka, Y., Takehisa, Y., Takaki, M., Taked, T., et al. (2001). Association study of CAG repeats in the KCNN3 gene in Japanese patients with schizophrenia, schizoaffective disorder and bipolar disorder. Psychiatry Res. 101, 203–207. doi: 10.1016/S0165-1781(01)00229-3
Varadi, M., Bertoni, D., Magana, P., Paramval, U., Pidruchna, I., Radhakrishnan, M., et al. (2024). AlphaFold protein structure database in 2024: providing structure coverage for over 214 million protein sequences. Nucleic Acids Res. 52, D368–D375. doi: 10.1093/NAR/GKAD1011
Wang, Y., Qu, L., Wang, X. L., Gao, L., Li, Z. Z., Gao, G. D., et al. (2015). Firing pattern modulation through SK channel current increase underlies neuronal survival in an organotypic slice model of Parkinson’s disease. Mol. Neurobiol. 51, 424–436. doi: 10.1007/S12035-014-8728-3
Waroux, O., Massotte, L., Alleva, L., Graulich, A., Thomas, E., Liégeois, J. F., et al. (2005). SK channels control the firing pattern of midbrain dopaminergic neurons in vivo. Eur. J. Neurosci. 22, 3111–3121. doi: 10.1111/J.1460-9568.2005.04484.X
Wittekindt, O., Jauch, A., Burgert, E., Schärer, L., Holtgreve-Grez, H., Yvert, G., et al. (1998). The human small conductance calcium-regulated potassium channel gene (hSKCa3) contains two CAG repeats in exon 1, is on chromosome 1q21.3, and shows a possible association with schizophrenia. Neurogenetics 1, 259–265. doi: 10.1007/S100480050038
Wolfart, J., Neuhoff, H., Franz, O., and Roeper, J. (2001). Differential expression of the small-conductance, calcium-activated potassium channel SK3 is critical for pacemaker control in dopaminergic midbrain neurons. J. Neurosci. 21, 3443–3456. doi: 10.1523/jneurosci.21-10-03443.2001
Keywords: SK3 channels, neurons, dopaminergic neurons, brain development, neurodegeneration
Citation: Padilla L, Del Val C, Neidre DB, Kokenge AS, Martinez JE, Teixeira AL, Zwir I and de Erausquin GA (2025) In silico prediction of the impact of genomic variations in the small conductance calcium activated potassium channel SK3 structure and function. Front. Neurosci. 19:1631536. doi: 10.3389/fnins.2025.1631536
Edited by:
Hua Lou, Case Western Reserve University, United StatesReviewed by:
Luye Qin, University of South Dakota, United StatesGuenson Chevalier, National Scientific and Technical Research Council (CONICET), Argentina
Copyright © 2025 Padilla, Del Val, Neidre, Kokenge, Martinez, Teixeira, Zwir and de Erausquin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gabriel A. de Erausquin, ZGVlcmF1c3F1aW5nQHV0aHNjc2EuZWR1
†These authors have contributed equally to this work
Agustín S. Kokenge1