Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Novel HCN2 Mutation Contributes to Febrile Seizures by Shifting the Channel's Kinetics in a Temperature-Dependent Manner

  • Yuki Nakamura,

    Affiliations Department of Pediatrics, Faculty of Medicine, Fukuoka University, Fukuoka, Japan, The Research Institute for the Molecular Pathomechanisms of Epilepsy, Fukuoka University, Fukuoka, Japan

  • Xiuyu Shi,

    Affiliations The Research Institute for the Molecular Pathomechanisms of Epilepsy, Fukuoka University, Fukuoka, Japan, Department of Pediatrics, Chinese PLA General Hospital, Beijing, China

  • Tomohiro Numata,

    Affiliation Department of Synthetic Chemistry and Biological Chemistry, Graduate School of Engineering, Kyoto University, Kyoto, Japan

  • Yasuo Mori,

    Affiliation Department of Synthetic Chemistry and Biological Chemistry, Graduate School of Engineering, Kyoto University, Kyoto, Japan

  • Ryuji Inoue,

    Affiliation Department of Physiology, Faculty of Medicine, Fukuoka University, Fukuoka, Japan

  • Christoph Lossin,

    Affiliation Department of Neurology, School of Medicine University of California Davis, Sacramento, California, United States of America

  • Tallie Z. Baram,

    Affiliation Departments of Anatomy & Neurobiology, Pediatrics, and Neurology, University of California Irvine, Irvine, California, United States of America

  • Shinichi Hirose

    hirose@fukuoka-u.ac.jp

    Affiliations Department of Pediatrics, Faculty of Medicine, Fukuoka University, Fukuoka, Japan, The Research Institute for the Molecular Pathomechanisms of Epilepsy, Fukuoka University, Fukuoka, Japan

Abstract

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channel-mediated currents, known as Ih, are involved in the control of rhythmic activity in neuronal circuits and in determining neuronal properties including the resting membrane potential. Recent studies have shown that HCN channels play a role in seizure susceptibility and in absence and limbic epilepsy including temporal lobe epilepsy following long febrile seizures (FS). This study focused on the potential contributions of abnormalities in the HCN2 isoform and their role in FS. A novel heterozygous missense mutation in HCN2 exon 1 leading to p.S126L was identified in two unrelated patients with FS. The mutation was inherited from the mother who had suffered from FS in a pedigree. To determine the effect of this substitution we conducted whole-cell patch clamp electrophysiology. We found that mutant channels had elevated sensitivity to temperature. More specifically, they displayed faster kinetics at higher temperature. Kinetic shift by change of temperature sensitivity rather than the shift of voltage dependence led to increased availability of Ih in conditions promoting FS. Responses to cyclic AMP did not differ between wildtype and mutant channels. Thus, mutant HCN2 channels cause significant cAMP-independent enhanced availability of Ih during high temperatures, which may contribute to hyperthermia-induced neuronal hyperexcitability in some individuals with FS.

Introduction

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are nonselective cation channels in the heart and brain [1][3], and are encoded by four genes (HCN1, 2, 3, and 4) [4]. HCN2 protein is expressed in regions including thalamus, brainstem, cerebral cortex, hippocampus, and amygdala [5], [6]. In hippocampus and cortex, HCN2 may co-assemble with HCN1, at least in part, to form heterotetrameric channels [7][13], which contributes to native Ih, a hyperpolarization-activated cation current. Unlike in the thalamus, where Ih contributes to regulating rhythmic activity [1], [14], hippocampal and cortical Ih contributes to the resting membrane potential [1], [15], [16], and is involved in dendritic summation and other parameters of neuronal excitability [17][20].

A role for Ih in abnormal neuronal excitability has been documented [14], [17], [19], [21][25]. In addition to thalamocortical seizures [14], [26], HCN channels have been implicated in limbic seizures [25], [27][29], including those that arise after prolonged febrile seizures [27]. In support of a role of HCN2 channels in epileptic seizures, a mutation in HCN2 was recently causally linked to general epilepsy [30].

Febrile seizures (FS) are the most common seizures in young children, and mutations in ion channels have been found in variants of FS that include epilepsy such as Dravet syndrome and Genetic epilepsy with FS plus (GEFS+), a clinical subset of FS [31][34]. However it remains unknown if abnormalities in HCN channel function contribute to typical FS.

The aims of this study were to (1) identify novel HCN2 mutations in children with FS, (2) examine differential temperature sensitivity between wildtype and mutant HCN2, and (3) examine the putative role of HCN2 mutations in FS. We analyzed HCN2 for mutations in 160 children with FS, and identified one novel mutation in two individuals with FS. Characterization of the HCN2 mutant channels showed that the activation kinetics of mutant HCN2 were faster in hyperthermic environments, compared with the wildtype, increasing the availability of Ih. Because HCN2-mediated Ih promotes neuronal firing, the results suggest that, in the presence of mutant HCN2, fever would induce neuronal hyperexcitability promoting FS.

Materials and Methods

Ethics statement

The present study was approved by the Ethics Review Committees of Fukuoka University. Parents of each patient and the parents themselves provided signed informed consent before the study.

Sample collection and mutation screening

This study was performed in Japan. 160 Japanese children diagnosed with FS based on the guidelines of the American Academy of Pediatrics [35] and 125 healthy control subjects were examined. The eight exons and exon/intron boundaries of HCN2 were amplified by polymerase chain reaction (PCR). The sequences were analyzed using a direct sequencing method with a 3130xl Genetic Analyzer (Applied Biosystems).

HCN2 mutagenesis

Human HCN2 cDNA was a gift from Dr. A. Ludwig (Friedrich-Alexander-University, Erlangen-Nuremberg, Germany). This construct contains the HCN2 open reading frame described by GenBank reference sequence NM_001194.3. Site-directed mutagenesis employed double-overlap PCR. The mutagenic primers were (forward) 5′-GAGCTCGGATCCACTAGTCCAGTGTGGTGG-3′ and (reverse) 5′-GCCCCGCGGCACaAGAACGACACCT-3′, where the altered base (c.377C>T) encoding the p.Ser126Leu (abbreviated as S126L hereafter) exchange is underlined and shown in lowercase letters.

Expression of HCN channels and electrophysiological recordings

HEK293 cells were grown in Eagle's Minimum Essential Medium (MEM) containing 0.1 mM non-essential amino acids and 10% fetal bovine serum at 37°C in a humidified 5% CO2 environment. Transient expression of HCN2 was achieved by transfection of 1.0 µg plasmid DNA with Lipofectamine2000 (Invitrogen); pIRES2-EGFP (Clontech, Mountainview, CA) was co-transfected as a fluorescence marker in a 1∶10 mass ratio. In experiments where wildtype (W) and mutant DNA (M) were transfected simultaneously, equal amounts of both were transfected totaling 1.0 µg. Assuming random assembly and equal subunit supply, stochastically, this produced a mix of different channel types with the following distribution: MMMM and WWWW: 6.25% each; MWWW and WMMM in all 4 iterations 25% each, MMWW: 37.5% [36]. Three hours after transfection, the cells were dissociated with 0.25% trypsin/EDTA (Invitrogen) and seeded onto coverslips.

Electrophysiological analyses commenced 24–48 hrs after transfection on a fluorescence-enabled inverted microscope (Nikon, Japan) equipped with a temperature-controlled TC-324B recording chamber (Warner Instruments, Holliston, MA). The cells were continually superfused with extracellular solution containing (in mM): 110 NaCl, 30 KCl, 0.5 MgCl2, 1.8 CaCl2, and 5 HEPES, at pH 7.4 (NaOH). Recording pipettes were pulled from borosilicate glass on a Narishige PC-10 (Japan), and fire-polished on a Narishige MF-830 microforge. Filled with internal solution comprising (in mM) 130 KCl, 10 NaCl, 0.5 MgCl2, 1 EGTA, and 5 HEPES, at pH 7.4 (KOH), the pipettes had resistances of ∼3 MΩ. The liquid junction potential was negligible (∼3 mV, [37]) and corrected, along with the pipette capacitance, using the amplifier's circuitry (EPC-8, Heka Instruments, Bellmore, NY) prior to recording. Data sampled at 7 kHz were 7-pole Bessel filtered at 3 kHz and digitized with a PCI-6221 digitizer (National Instruments, Austin, TX). Series resistance was compensated by 40–70%. All experiments utilized a holding potential of −40 mV, which approximated the K+ Nernst potential (∼−38 mV) between 25 and 38°C. Cellular capacitance measurements relied on the amplifier's readout after cell capacitance compensation. Leak subtraction was not necessary because the holding current at EK+ (−40 mV) was commonly ≤10 pA.

Fitting and statistical analysis

Data were analyzed off-line using the Strathclyde Electrophysiological Software (John Dempster, University of Strathclyde, Glasgow, UK, downloadable at http://spider.science.strath.ac.uk/sipbs/software_ses.htm). Peak amplitudes from tail currents acquired immediately after stepping back to the holding potential (−40 mV) were normalized to the maximal current, plotted against the test potential, and fitted with a Boltzmann functionwhere y is the normalized tail current, V is the test potential, V1/2 is the half-activation potential and k is the slope factor. Activation kinetics were fitted with the double-exponential function y = Afast·exp(−t/τfast)+Aslow·exp(−t/τslow)+C. Data are presented as mean ± SEM. Statistical analysis employed one-way ANOVA, and significance was set at p<0.05.

Results

HCN2 gene analysis

An HCN2 mutation was identified as a heterozygous missense mutation (c. 377C>T) in two unrelated children with FS. The mother (I–2), who also suffered from FS in childhood, presented with the same HCN2 alteration as her child (II–1); the father (I–1), unaffected, did not, suggesting that the change was inherited (Fig. 1 A). The mutation is predicted to substitute a polar serine at position 126 with a non-polar leucine (S126L). Position 126 is within the intracellular N–terminus, close to the S1 trans-membrane segment (Fig. 1B). The amino acid in this position varies among the four different human HCN channels. However, an alignment of HCN2 protein sequences across species shows that p.Ser126 is evolutionarily conserved (Fig. 1C). This is reiterated by the results of a SIFT in silico analysis, which predicts S126L to be damaging to HCN2 function [38]. Mutations in other FS candidate genes such as SCN1A, SCN2A, SCN1B, SCN2B, GABRA1, GABRB2, and GABRG2 were not found in any of the patients.

thumbnail
Figure 1. Genetic analysis of S126L.

A left, Family pedigree; the proband is indicated by an arrow. A right, Double-strand partial sequencing data for HCN2 along with the amino acid translation for the unaffected father (I-1) and the mother and child with febrile seizure (I-2 and II-1). The missense mutation c.377C>T (S126L) is indicated by the arrow. B, Putative transmembrane topology of one HCN2 subunit showing the approximate position of S126 in the N-terminus. C, Sequence alignment of the N-terminus. The rectangle indicates the altered residue.

https://doi.org/10.1371/journal.pone.0080376.g001

Electrophysiological characterization at 25°C and effects of temperature on wildtype, S126L and heteromeric HCN2 channels

FS, by definition, occur during fever. If S126L causes or contributes to the phenotype, one would expect that evidence of HCN2 dysfunction presents itself preferentially at or above 38°C. To examine the effect of fever on wildtype, wildtype/S126L (heteromeric) channels, and S126L homomeric channels, we recorded whole cell currents at 25 and 38°C, using 3-s voltage steps to hyperpolarizing potentials (−140 to −40 mV) from a −40 mV holding potential (Fig. 2 A). S126L and heteromeric channels exhibited a voltage-dependent inward current that resembled wildtype channels at 25°C and 38°C. Cells stimulated in this fashion produced voltage-dependent inward currents followed by tail currents when stepped back to the −40 mV hold. We used the tail currents (Fig. 2 A, arrows) to deduce the channels' voltage dependence and voltage sensitivity of activation (expressed in Table 1), by normalizing to the maximal amplitude and fitting it to a Boltzmann function (Fig. 2B, C, and D).

thumbnail
Figure 2. Current traces and activation curves of HCN2 wildtype, S126L and heretomeric channels at 25, 35 and 38°C.

A, Whole-cell currents of HCN2 wildtype (top), S126L (middle), and heteromeric (bottom) channels at 25 and 38°C. The currents were recorded in response to voltage pulses from a holding potential of −40 mV to test potentials between −140 to −40 mV. B, C, D, Conductance-voltage curves for wildtype, S126L and heteromeric channels at 25, 35 and 38°C. Continuous lines show Boltzmann fits of the experimental data.

https://doi.org/10.1371/journal.pone.0080376.g002

thumbnail
Table 1. V1/2 for wildtype, S126L and heteromeric channels at 25°C and 38°C.

https://doi.org/10.1371/journal.pone.0080376.t001

The half maximal activation voltages (V1/2) at 25°C for wildtype (n = 12), S126L (n = 12) and heteromeric channels (n = 8) were −95.1±1.9, −99.0±2.1, and −94.9±1.3 mV, respectively, and the respective slope factors (k) were 10.3±0.9, 10.3±0.6, and 9.8±1.9. V1/2 at 38°C for wildtype (n = 8), S126L (n = 8) and heteromeric channels (n = 7) were −92.7±2.0, −90.1±2.6 and −91.5±2.9 mV, respectively, and k were 9.9±1.5, 10.0±1.0 and 8.1±0.7.

No significant differences were found when V1/2 of the wildtype, S126L and heteromeric channels were compared at same temperature. However, raising the temperature (from 25°C to 38°C) led to depolarizing shift (ΔV1/2) in S126L channels (ΔV1/2 = +8.9 mV) that was significantly larger than in wildtype channels (ΔV1/2: 2.4 mV, p<0.05), which speaks for increased temperature sensitivity in the mutant (Table 1). The current-voltage relationship was determined by the steady-state current in response to step pulse stimuli for 3 seconds (Fig. 3 A, B). At 25°C, the current density of S126L at −140 mV (135.9±16.2 pA/pF) was significant larger than that of wildtype (72.1±14.6 pA/pF, p<0.05). Although there was no significant difference between the channels at 38°C, the current density of S126L at −140 mV (150.2±27.1 pA/pF) and that of heteromeric channels (135.6±30.4 pA/pF) were larger than that of wildtype (93.3±18.0 pA/pF, p>0.05) (Fig. 3B). Current traces of these channels were fitted by two exponential functions. The fast and slow time constants (tau fast and tau slow, respectively) describing the activation kinetics at 25°C were voltage-dependent and smaller (faster) at more hyperpolarized voltages (Fig. 3 C, top). Between −90 and −100 mV, tau fast was significantly larger (slower) in S126L channels (at −100 mV: 252.6±13 ms, n = 10 wildtype, compared with 391.9±64.3 ms, n = 10 S126L, p<0.05, Table S1). In −90 mV, tau fast was significantly larger in heteromeric channels (295.1±29.0 ms, n = 10 wildtype, compared with 473.9±72.4 ms, n = 7 hetero, p<0.05)

thumbnail
Figure 3. Kinetic properties in HCN2 wildtype, S126L and heteromeric channels at 25 and 38°C.

A, B, Current density at 25 (A) and 38°C (B). A, The current density of S126L channels at 25°C is significantly larger than that of wildtype at −140 mV (*p<0.05). C, Fast activation time constants for wildtype, S126L and heteromeric channels at 25 and 38°C. denotes tau fast of S126L channels at 25°C is significantly larger than that of wildtype channels at −100 mV. * denotes tau fast of S126L channels at 25°C is significantly larger than that of wildtype channels at −90 mV. ** denotes tau fast of S126L channels at 25°C is significantly larger than that of heteromeric channels at −90 mV. †† denotes tau fast of S126L channels at 38°C is significantly smaller than that of wildtype channels at −90 and −100 mV (†, *, **, †† p<0.05). D, Alterations in tau fast by increasing temperature at −140, −100 and −90 mV. E, F, Relative amplitudes of fast exponential component as function of voltage at wildtype, S126L and heteromeric channels at 25 (E) and 38°C (F). E, Relative amplitude of S126L channels at 25°C was significantly larger than that of wildtype channels at −90 mV (* p<0.05). F, Relative amplitude of S126L and heteromeric channels at 38°C was significantly larger than that of wildtype channels at −90 mV (*, **p<0.05).

https://doi.org/10.1371/journal.pone.0080376.g003

Tau fast and slow at 38°C (Fig. 3 C, bottom) also exhibited voltage dependence and were faster (smaller) than those observed at 25°C. In contrast to the observation at 25°C, tau fast of S126L and heteromeric channels was faster than that of wild type from −90 to −140 mV. This phenomenon was most prominent at −90 mV. To better characterize this behavior, we subtracted tau fast at 38°C from that of 25°C in −140, −100 and −90 mV (Fig. 3 D). In −140 mV, “tau fast (25°C) – tau fast (38°C)” (Δtau) of wildtype was nearly equal to that of S126L and heteromeric channels. However at −100 and −90 mV, Δtau of S126L and heteromeric channels was larger than that of wildtype channels. These results indicate that the shift of the activation kinetics for S126L and heteromeric channels due to raising the temperature is larger in mutant channels, and that this is more evident in relatively depolarized (physiological) voltage ranges.

Because the relative contribution of tau fast and slow to the activation kinetics influence channel behavior [7], [39], we calculated the relative amplitude of the fast and slow exponential components of Ih (Fig. 3 E, F, Table S1). At 25°C, in S126L channels, the contribution of the fast component was significantly higher at −90 mV compared with wildtype channels (Fig. 3 E, p<0.05). These data indicate that the fast component (Afast) in S126L channels is increased at depolarized voltages. The contribution of the fast component to total Ih in each channel at 38°C was less voltage dependent than at 25°C (Fig. 3 F vs. 3 E). The contribution of the fast component of S126L and heteromeric channels was significantly increased compared with wildtype channels at −90 mV (Fig. 3 F, p<0.05). These results suggest that the fast component (Afast) in S126L and heteromeric channels is increased at more physiological potentials and at hyperthermic temperatures, and that, compared to wildtype, S126L and heteromeric channels exhibit altered temperature sensitivity. In summary, the shift of the activation kinetics by a raise in temperature was highly influenced rather than that of the conductance and the voltage dependence.

Comparison of temperature dependence in wildtype, S126L and heteromeric channels

To evaluate the temperature dependence in each channel, we measured the Ih amplitude at 25, 35 and 38°C. Figure 4 shows that the current density was increased by a rise in temperature in all channels. Between 25 to 35°C, increased current density was similar in S126L (slope: 1.2±3.0 pA·pF−1·K−1), heteromeric (2.0±3.1 pA·pF−1·K−1) and wildtype (0.8±3.0 pA·pF−1·K−1) channels (Table S2). However, raising the temperature from 35 to 38°C - a temperature commonly associated with fever - led to highly divergent effects: whereas a gradual increase continued in the current density of wildtype (0.6±10.1 pA·pF−1·K−1) channel, a drastic, steep increase in current density occurred in the S126L (4.7±12.6 pA·pF−1·K−1) and heteromeric (6.7±14.7 pA·pF−1·K−1) channels. These results suggest that specific modulation by this mutation was occurred between 35 and 38°C.

thumbnail
Figure 4. Temperature dependence of HCN2 wildtype, S126L and heteromeric channels.

The effect of temperature on wildtype (n = 7–11), S126L (n = 7–12) and heteromeric (n = 7–10) channels at −140 mV. Elevated temperature raises the current densities in all channels, but only S126L-containing channels show a distinctly steeper slope from 35 to 38°C.

https://doi.org/10.1371/journal.pone.0080376.g004

Role of cAMP in mutant HCN2 channel function

Because cAMP is a direct modulator of HCN channels [3], [40] and may influence the availability of Ih during hyperthermia, we investigated the effects of cAMP on S126L in comparison to wildtype channels. As expected, the activation curves of both channels shifted towards more depolarized potentials following application of cAMP at 25°C (Fig. 5 A, B). In addition, there was a tendency for higher sensitivity of the mutant channels to low doses (2 µM) of cAMP (Table S3, p<0.05).

thumbnail
Figure 5. Effects of cAMP concentrations and raising temperature.

A, B, Conductance-voltage curves of wildtype (A) and S126L (B) channels in the absence and presence of cAMP at concentrations of 1, 2, and 10 µM. C, D, Comparison of conductance-voltage curves between at 25°C and at 38°C in 2 µM cAMP of wildtype (C) and S126L (D) channels. E, F, Hill fits of the cAMP dose-response as measured by half-maximal activation at 25°C and 38°C.

https://doi.org/10.1371/journal.pone.0080376.g005

More specifically, focused on 2 µM cAMP – the lowest dose leading to consistent changes in the activation curves in all channels at 38°C – resulted in a +3.5 mV(p = 0.40) right-shift of the activation curve for S126L channels. Wildtype channel, on the other hand, showed almost no shift (−0.1 mV, p = 0.98, Table S4). Although no significant difference was noted between the wildtype and S126L channels, this change produces clear conductance differences that enables the mutant to mediate more current (Fig. 5 C, D).

The data produced in Figure 5A and B were re-plotted to produce cAMP dose-response curves based on the half-maximal activation (Fig. 5E, F). Hill fits of the data show that S126L channels were activated at more positive voltages than wildtype channels at 38°C and that both channels at 38°C were activated by lower cAMP concentrations compared with 25°C. However K1/2 and h values for S126L channels were close to these values for wildtype channels (Table S5). These results suggest that the altered profile of S126L channels is a result of changes in voltage dependence and temperature sensitivity rather than in cAMP dependence.

Discussion

Here we identified a novel heterozygous HCN2 missense mutation (S126L) in children with FS and characterized its functional properties. We found increased current densities and accelerated kinetics in HCN2 channel with the S126L mutant. These alterations produce a mutant that allows more positive charge to enter the cell more quickly compared to healthy cells, specifically at elevated temperature, resulting in depolarization and excitation. Our results suggest that current mediated by S126L-containing channels may augment neuronal excitability during hyperthermia.

One may ask whether our findings are truly representative of what is seen in the patients when our data show the strongest effects in setups using only S126L channels, when the patients are heterozygous for the alteration. To emulate the patient internal environment, we used setups with a mix of wildtype and mutant DNA. This commonly produced behaviors that lay in between wildtype and purely-mutant, which may relate to the stochastic distribution of the ensuing channel possibilities. Closer analysis of possible subunit compositions for 1∶1 transfection experiments reveals that only 37.5% of the channels were composed of 2 wildtype and 2 mutant subunits. Exactly half of the channels (50%) contained 3 subunits of one type and the remaining subunit was of the other type. Most importantly, however, even in WT/mutant transfections, 12.5% of the channels were in fact homomeric, either wildtype or mutant. A “dose-dependent” effect is therefore to be expected, and our data confirm this, because all heteromeric channels produced electrophysiological deviations similar (usually somewhat reduced) or equal to what was seen in the setups using only mutant DNA. We have no way of examining the stoichiometry of the mutant channels in the native environment. It is uncertain, for example, whether, in the patients, S126L protein assembles with wildtype and/or mutant protein as readily as it does with healthy protein. Conversely, we cannot exclude that S126L does not enhance subunit interaction between the mutants, which could explain the augmented current densities as we elaborate below.

In most cortical and hippocampal neurons, Ih channels are constitutively open near the resting membrane potential [41][43], and contribute to its maintenance [1], [16]. In addition, Ih controls the membrane potential via a depolarizing inward current and by facilitating hyperpolarization [44], [45]. The large shift of the kinetics in homomeric S126L and heteromeric channels near the resting membrane potential at high temperatures could lead to early activation of Ih and disrupt this Ih-based fine-tuning, promoting neuronal hyperexcitability.

In hippocampal CA1 neurons of immature rats, hyperthermia may reduce GABAA-receptor-mediated synaptic inhibition [46]. Such decrease, coupled with the shift in the activation kinetics of HCN2 S126L and heteromeric channels during hyperthermia could lead to an imbalance between neuronal excitation and inhibition, accelerating the development of FS.

Increased current densities at strongly hyperpolarized potentials have been reported before [24], [47], [48]. At this time, we cannot explain this phenomenon. Three factors contribute to macroscopic current (I), namely the number of contributing channels (N), single-channel current (i), and open probability (p). We deem it unlikely that either i or p are altered, simply based on the location of p.Ser126 within the HCN2 structure. We expect to find p.Ser126 in the periphery of the channel, although a definitive crystal structure is so far unavailable. This position would situate p.Ser126 far from the pore, which makes it hard to conceive an influence on the HCN2 single-channel conductance, which defines i in a voltage-dependent fashion. Similarly, p.Ser126 would be unlikely to impact the open probability of the channel. Unfortunately, enhanced surface expression is not a plausible explanation for the increase in Ih either, as this would imply a voltage-dependent change in the number of channels in the membrane. There is, however, the possibility of a change in functional-channel availability. How this could be mediated is currently unclear. Possible scenarios, given the residue's mapping to the N-terminus, include abnormal phosphorylation as well as alterations in protein-protein interactions. The latter is particular interesting, owing to a known importance of the N-terminus in inter-subunit interaction and formation of functional homo- and heteromeric channels [9]. The substitution of a polar serine with a nonpolar leucine might modify intermolecular forces or hydrogen bonds with other amino acids, leading to aberrant interaction with other HCN subunits required for tetramerization. Indeed, heteromerization of HCN2 with HCN1 [9], [10] is highly regulated in the brain [12], [13]. Augmented heteromerization after experimental FS [13] has been correlated with a shift in Ih properties that is similar to the one we found in S126L channels [24]. These data raise the possibility that S126L enhances HCN2∶HCN1 interaction, leading to aberrant channel function.

The mutation might interfere also with the interaction of HCN2 with auxiliary proteins. HCN channels interact with filamin [49], TRIP8b [50], [51] and enzymes such as p38 MAPK [52] and Src kinase [53]. Whereas the large majority of these proteins interact with HCN1 [54], tamalin, KCNE2, and others interact with HCN2. S126L might impair these interactions, enhancing a kinetic shift of the channel and promoting FS.

In conclusion, a mutation in HCN2 has been identified, with a potential to contribute to FS in a subgroup of children. It is now the 3rd HCN mutation found in seizures and epilepsy, all of which affect HCN2. This underlies the importance of HCN channels to normal and aberrant neuronal excitability, and the specific properties of HCN2 that make mutations in these channels (or their absence) likely to result in a hyperexcitable network.

Supporting Information

Table S1.

Kinetic parameters for HCN2 current activation at −90 mV based on double-exponential fits. For the sake of clarity, only fast kinetic parameters are listed. * indicates p<0.05 versus wildtype.

https://doi.org/10.1371/journal.pone.0080376.s001

(DOC)

Table S2.

Development of current density versus temperature. The difference in the averaged current densities between two temperature points at −140 mV was calculated as a slope value.

https://doi.org/10.1371/journal.pone.0080376.s002

(DOC)

Table S3.

cAMP sensitivity in wildtype and mutant channels: comparison of the shift in the voltage dependence of activation. The ΔV1/2 describes the cAMP-induced voltage shift in half-maximal activation compared to control levels. * indicates p<0.05 compared to the control.

https://doi.org/10.1371/journal.pone.0080376.s003

(DOC)

Table S4.

Temperature dependence of cAMP sensitivity in wildtype and mutant channels.

https://doi.org/10.1371/journal.pone.0080376.s004

(DOC)

Table S5.

cAMP dose-response curves based on half-maximal activation voltages.

https://doi.org/10.1371/journal.pone.0080376.s005

(DOC)

Acknowledgments

We thank Dr. A. Ludwig for providing human HCN2 cDNA. We are grateful to Ms. M. Yonetani and A. Hamachi for their technical assistance.

Author Contributions

Conceived and designed the experiments: YN RI TZB SH. Performed the experiments: YN XS TN. Analyzed the data: YN TN CL. Contributed reagents/materials/analysis tools: YN YM RI SH. Wrote the paper: YN CL TZB SH.

References

  1. 1. Biel M, Wahl-Schott C, Michalakis S, Zong X (2009) Hyperpolarization-activated cation channels: from genes to function. Physiol Rev 89: 847–885.
  2. 2. Accili EA, Proenza C, Baruscotti M, DiFrancesco D (2002) From funny current to HCN channels: 20 years of excitation. News Physiol Sci 17: 32–37.
  3. 3. Ludwig A, Zong X, Jeglitsch M, Hofmann F, Biel M (1998) A family of hyperpolarization-activated mammalian cation channels. Nature 393: 587–591.
  4. 4. Santoro B, Liu DT, Yao H, Bartsch D, Kandel ER, et al. (1998) Identification of a gene encoding a hyperpolarization-activated pacemaker channel of brain. Cell 93: 717–729.
  5. 5. Santoro B, Chen S, Luthi A, Pavlidis P, Shumyatsky GP, et al. (2000) Molecular and functional heterogeneity of hyperpolarization-activated pacemaker channels in the mouse CNS. J Neurosci 20: 5264–5275.
  6. 6. Notomi T, Shigemoto R (2004) Immunohistochemical localization of Ih channel subunits, HCN1-4, in the rat brain. J Comp Neurol 471: 241–276.
  7. 7. Chen S, Wang J, Siegelbaum SA (2001) Properties of hyperpolarization-activated pacemaker current defined by coassembly of HCN1 and HCN2 subunits and basal modulation by cyclic nucleotide. J Gen Physiol 117: 491–504.
  8. 8. Ulens C, Tytgat J (2001) Functional heteromerization of HCN1 and HCN2 pacemaker channels. J Biol Chem 276: 6069–6072.
  9. 9. Proenza C, Tran N, Angoli D, Zahynacz K, Balcar P, et al. (2002) Different roles for the cyclic nucleotide binding domain and amino terminus in assembly and expression of hyperpolarization-activated, cyclic nucleotide-gated channels. J Biol Chem 277: 29634–29642.
  10. 10. Much B, Wahl-Schott C, Zong X, Schneider A, Baumann L, et al. (2003) Role of subunit heteromerization and N-linked glycosylation in the formation of functional hyperpolarization-activated cyclic nucleotide-gated channels. J Biol Chem 278: 43781–43786.
  11. 11. Santoro B, Baram TZ (2003) The multiple personalities of h-channels. Trends Neurosci 26: 550–554.
  12. 12. Zha Q, Brewster AL, Richichi C, Bender RA, Baram TZ (2008) Activity-dependent heteromerization of the hyperpolarization-activated, cyclic-nucleotide gated (HCN) channels: role of N-linked glycosylation. J Neurochem 105: 68–77.
  13. 13. Brewster AL, Bernard JA, Gall CM, Baram TZ (2005) Formation of heteromeric hyperpolarization-activated cyclic nucleotide-gated (HCN) channels in the hippocampus is regulated by developmental seizures. Neurobiol Dis 19: 200–207.
  14. 14. Kuisle M, Wanaverbecq N, Brewster AL, Frere SG, Pinault D, et al. (2006) Functional stabilization of weakened thalamic pacemaker channel regulation in rat absence epilepsy. J Physiol 575: 83–100.
  15. 15. Lupica CR, Bell JA, Hoffman AF, Watson PL (2001) Contribution of the hyperpolarization-activated current (I(h)) to membrane potential and GABA release in hippocampal interneurons. J Neurophysiol 86: 261–268.
  16. 16. Maccaferri G, McBain CJ (1996) The hyperpolarization-activated current (Ih) and its contribution to pacemaker activity in rat CA1 hippocampal stratum oriens-alveus interneurones. J Physiol (London) 497: 119–130.
  17. 17. Noam Y, Bernard C, Baram TZ (2011) Towards an integrated view of HCN channel role in epilepsy. Curr Opin Neurobiol
  18. 18. Fan Y, Fricker D, Brager DH, Chen X, Lu HC, et al. (2005) Activity-dependent decrease of excitability in rat hippocampal neurons through increases in I(h). Nat Neurosci 8: 1542–1551.
  19. 19. Poolos NP, Migliore M, Johnston D (2002) Pharmacological upregulation of h-channels reduces the excitability of pyramidal neuron dendrites. Nat Neurosci 5: 767–774.
  20. 20. Magee JC (1998) Dendritic hyperpolarization-activated currents modify the integrative properties of hippocampal CA1 pyramidal neurons. J Neurosci 18: 7613–7624.
  21. 21. Lewis AS, Chetkovich DM (2011) HCN channels in behavior and neurological disease: too hyper or not active enough? Mol Cell Neurosci 46: 357–367.
  22. 22. Kanyshkova T, Meuth P, Bista P, Liu Z, Ehling P, et al. (2011) Differential regulation of HCN channel isoform expression in thalamic neurons of epileptic and non-epileptic rat strains. Neurobiol Dis
  23. 23. Santoro B, Lee JY, Englot DJ, Gildersleeve S, Piskorowski RA, et al. (2010) Increased seizure severity and seizure-related death in mice lacking HCN1 channels. Epilepsia 51: 1624–1627.
  24. 24. Chen K, Aradi I, Thon N, Eghbal-Ahmadi M, Baram TZ, et al. (2001) Persistently modified h-channels after complex febrile seizures convert the seizure-induced enhancement of inhibition to hyperexcitability. Nat Med 7: 331–337.
  25. 25. McClelland S, Flynn C, Dube C, Richichi C, Zha Q, et al. (2011) Neuron-restrictive silencer factor-mediated hyperpolarization-activated cyclic nucleotide gated channelopathy in experimental temporal lobe epilepsy. Ann Neurol 70: 454–464.
  26. 26. Coulon P, Herr D, Kanyshkova T, Meuth P, Budde T, et al. (2009) Burst discharges in neurons of the thalamic reticular nucleus are shaped by calcium-induced calcium release. Cell Calcium 46: 333–346.
  27. 27. McClelland S, Dube CM, Yang J, Baram TZ (2011) Epileptogenesis after prolonged febrile seizures: mechanisms, biomarkers and therapeutic opportunities. Neurosci Lett 497: 155–162.
  28. 28. Shin M, Brager D, Jaramillo TC, Johnston D, Chetkovich DM (2008) Mislocalization of h channel subunits underlies h channelopathy in temporal lobe epilepsy. Neurobiol Dis 32: 26–36.
  29. 29. Jung S, Jones TD, Lugo JN Jr, Sheerin AH, Miller JW, et al. (2007) Progressive dendritic HCN channelopathy during epileptogenesis in the rat pilocarpine model of epilepsy. J Neurosci 27: 13012–13021.
  30. 30. DiFrancesco JC, Barbuti A, Milanesi R, Coco S, Bucchi A, et al. (2011) Recessive loss-of-function mutation in the pacemaker HCN2 channel causing increased neuronal excitability in a patient with idiopathic generalized epilepsy. J Neurosci 31: 17327–17337.
  31. 31. Hawkins NA, Martin MS, Frankel WN, Kearney JA, Escayg A (2011) Neuronal voltage-gated ion channels are genetic modifiers of generalized epilepsy with febrile seizures plus. Neurobiol Dis 41: 655–660.
  32. 32. Wimmer VC, Reid CA, Mitchell S, Richards KL, Scaf BB, et al. (2010) Axon initial segment dysfunction in a mouse model of genetic epilepsy with febrile seizures plus. J Clin Invest 120: 2661–2671.
  33. 33. Escayg A, Goldin AL (2010) Sodium channel SCN1A and epilepsy: mutations and mechanisms. Epilepsia 51: 1650–1658.
  34. 34. Wallace RH, Hodgson BL, Grinton BE, Gardiner RM, Robinson R, et al. (2003) Sodium channel alpha1-subunit mutations in severe myoclonic epilepsy of infancy and infantile spasms. Neurology 61: 765–769.
  35. 35. Steering Committee on Quality Improvement and Management SoFS (2008) Febrile seizures: clinical practice guideline for the long-term management of the child with simple febrile seizures. Pediatrics 121: 1281–1286.
  36. 36. Hille B (2001) Ion channels of excitable membranes. Sunderland, Mass.: Sinauer. xviii, 814 p. p.
  37. 37. Barry PH (1994) JPCalc, a software package for calculating liquid junction potential corrections in patch-clamp, intracellular, epithelial and bilayer measurements and for correcting junction potential measurements. J Neurosci Methods 51: 107–116.
  38. 38. Kumar P, Henikoff S, Ng PC (2009) Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc 4: 1073–1081.
  39. 39. Surges R, Brewster AL, Bender RA, Beck H, Feuerstein TJ, et al. (2006) Regulated expression of HCN channels and cAMP levels shape the properties of the h current in developing rat hippocampus. Eur J Neurosci 24: 94–104.
  40. 40. Wainger BJ, DeGennaro M, Santoro B, Siegelbaum SA, Tibbs GR (2001) Molecular mechanism of cAMP modulation of HCN pacemaker channels. Nature 411: 805–810.
  41. 41. Aponte Y, Lien CC, Reisinger E, Jonas P (2006) Hyperpolarization-activated cation channels in fast-spiking interneurons of rat hippocampus. J Physiol 574: 229–243.
  42. 42. Ludwig A, Budde T, Stieber J, Moosmang S, Wahl C, et al. (2003) Absence epilepsy and sinus dysrhythmia in mice lacking the pacemaker channel HCN2. EMBO J 22: 216–224.
  43. 43. Rodrigues AR, Oertel D (2006) Hyperpolarization-activated currents regulate excitability in stellate cells of the mammalian ventral cochlear nucleus. J Neurophysiol 95: 76–87.
  44. 44. Nolan MF, Dudman JT, Dodson PD, Santoro B (2007) HCN1 channels control resting and active integrative properties of stellate cells from layer II of the entorhinal cortex. J Neurosci 27: 12440–12451.
  45. 45. Solomon JS, Nerbonne JM (1993) Hyperpolarization-activated currents in isolated superior colliculus-projecting neurons from rat visual cortex. J Physiol (London) 462: 393–420.
  46. 46. Qu L, Liu X, Wu C, Leung LS (2007) Hyperthermia decreases GABAergic synaptic transmission in hippocampal neurons of immature rats. Neurobiol Dis 27: 320–327.
  47. 47. Dibbens LM, Reid CA, Hodgson B, Thomas EA, Phillips AM, et al. (2010) Augmented currents of an HCN2 variant in patients with febrile seizure syndromes. Ann Neurol 67: 542–546.
  48. 48. Dyhrfjeld-Johnsen J, Morgan RJ, Foldy C, Soltesz I (2008) Upregulated H-current in hyperexcitable CA1 dendrites after febrile seizures. Front Cell Neurosci 2: 2.
  49. 49. Gravante B, Barbuti A, Milanesi R, Zappi I, Viscomi C, et al. (2004) Interaction of the pacemaker channel HCN1 with filamin A. J Biol Chem 279: 43847–43853.
  50. 50. Lewis AS, Schwartz E, Chan CS, Noam Y, Shin M, et al. (2009) Alternatively spliced isoforms of TRIP8b differentially control h channel trafficking and function. J Neurosci 29: 6250–6265.
  51. 51. Santoro B, Wainger BJ, Siegelbaum SA (2004) Regulation of HCN channel surface expression by a novel C-terminal protein-protein interaction. J Neurosci 24: 10750–10762.
  52. 52. Poolos NP, Bullis JB, Roth MK (2006) Modulation of h-channels in hippocampal pyramidal neurons by p38 mitogen-activated protein kinase. J Neurosci 26: 7995–8003.
  53. 53. Zong X, Eckert C, Yuan H, Wahl-Schott C, Abicht H, et al. (2005) A novel mechanism of modulation of hyperpolarization-activated cyclic nucleotide-gated channels by Src kinase. J Biol Chem 280: 34224–34232.
  54. 54. Lewis AS, Estep CM, Chetkovich DM (2010) The fast and slow ups and downs of HCN channel regulation. Channels (Austin) 4: 215–231.