A Kinetic Map of the Homomeric Voltage-Gated Potassium Channel (Kv) Family

The voltage-gated potassium (Kv) channels, encoded by 40 genes, repolarize all electrically excitable cells, including plant, cardiac, and neuronal cells. Although these genes were fully sequenced decades ago, a comprehensive kinetic characterization of all Kv channels is still missing, especially near physiological temperature. Here, we present a standardized kinetic map of the 40 homomeric Kv channels systematically characterized at 15, 25, and 35°C. Importantly, the Kv kinetics at 35°C differ significantly from commonly reported kinetics, usually performed at room temperature. We observed voltage-dependent Q10 for all active Kv channels and inherent heterogeneity in kinetics for some of them. Kinetic properties are consistent across different host cell lines and conserved across mouse, rat, and human. All electrophysiology data from all Kv channels are made available through a public website (Channelpedia). This dataset provides a solid foundation for exploring kinetics of heteromeric channels, roles of auxiliary subunits, kinetic modulation, and for building accurate Kv models.


INTRODUCTION
Ion channels (ICs) are proteins that selectively allow ions to diffuse through the cell membrane, creating an electrical potential across the membrane. They are classified in terms of which ion passes through the channel (sodium, potassium, chloride, calcium etc.) and by gating activity (voltage gating, ligand gating, other forms of gating; Hille, 2001;Alberts et al., 2008).
The ∼350 IC types expressed in the mammalian brain include 145 voltage-gated channels, of which 40 are voltage-gated potassium (Kv) channels, divided into 12 sub-families, Kv1-Kv12 ( Figure 1A). The Kv channels are expressed in many different tissues/organs, including muscle, heart, and brain (Post et al., 1992;Lai and Jan, 2006;Li and Dong, 2010;Southan et al., 2016), and also exhibit specific expression patterns at the subcellular level (Rasband, 2010;Jensen et al., 2011;Trimmer, 2015). In all these tissues, Kv channels are activated by variations in the voltage across the cell membrane. These changes regulate the return of the membrane to a resting state (hyperpolarization) after depolarization, thus controlling the excitability of different cell types (neurons, cardiomyocytes, skeletal muscle cells, etc.).
Functionally, Kv channels are known to regulate the threshold potential for firing, the duration of action potentials and, the firing rates (Storm, 1988;Gabel and Nisenbaum, 1998;Glazebrook et al., 2002;Begum et al., 2016). They are also involved in cell proliferation, and play a significant role in neurotoxicity, neuroprotection, and neuroregulation (Kaczmarek, 2006;Shah and Aizenman, 2014). Hence, malfunction in Kv channels are implicated in a range of neurological diseases known as channelopathies (Graves and Hanna, 2005;Kullmann and Waxman, 2010).
Structurally, mammalian Kv channels as homotetramers, consist of four identical alpha (α)-subunits arranged around a central axis that forms a pore (Coetzee et al., 1999). Each αsubunit is made up of six α-helical transmembrane spanning segments (S1-S6), five loops connecting successive segments, and cytoplasmic regions at the N and C terminal ends. The first four segments (S1-S4) form the voltage sensor domain (VSD). Within VSD, the fourth segment, S4, which contains a positively charged amino acid (arginine or lysine) at every third position (Bezanilla, 2000), is considered the main voltage sensor. The pore is formed by the S5-S6 segments and the corresponding connecting loop, which contains a conserved sequence (glycine-tyrosine-glycine or glycine-phenylalanine-glycine) that acts as a K + ion selectivity filter (Heginbotham et al., 1994).
In mammalian genomes, the Kv channels are encoded by 40 genes. Each Kv gene encodes a corresponding protein (αsubunit) with a unique structure ( Figure 1B). While all Kv channels share a similar core structure, differences in their connecting loops, and in the length and sequence of their N and C terminals, make each one of them unique. The genes coding for the Kv channels have been cloned and studied in cell lines for many decades. However, a comprehensive and standardized kinetic characterization of all homomeric Kv channels is still missing. The literature focuses on just a few Kv channels (e.g., Kv7.1,Kv11.1,Kv1.5,Kv1.1,Kv1.3,Kv1.2,Kv2.1,etc.), and neglects many others (e.g., Kv1.7,Kv1.8,Kv6.x,Kv9.x,Kv12.x; Table 1). Moreover, different studies have used different species, cloning procedures, host cell lines, cell culture conditions, intracellular and extracellular solutions, and stimulation protocols, making it difficult to compare and integrate their data and sometimes producing discrepant results. In addition, the raw electrophysiology data needed to build realistic IC models are not publicly available; available data are only in the form of figures or extracted features. Eventually, the majority of previous studies on Kv channels has been Abbreviations: AP, Action potential; APC, Automated patch clamp; CHO, Chinese hamster ovary; CV1, Cercopithecus aethiops kidney; FT, Flp-In TM T-Rex TM ; HEK, Human embryonic kidney; HSP, Heat shock proteins; IC, ion channel; Kv, Voltage-gated potassium channel; RT, Room temperature. conducted at room temperature (RT), which can vary from 18 • C (Heinemann et al., 1996) to 28 • C (Hatton et al., 2001) that could produce inconsistencies in results.
The thermal sensitivity of any biological process can be described by its temperature coefficient (Q 10 ). Classically, it is defined as the ratio of a reaction rate (α) measured at two temperatures 10 degrees apart (Bělehrádek, 1935). In ion channel research, instead of reaction rates, current amplitudes or time constants are often used to calculate Q 10 value. A single value of Q 10 is typically reported to indicate the temperature dependence of a channel. However, it is important to note that the Q 10 for an IC can be different for different temperature ranges (Beam and Donaldson, 1983). Moreover, the kinetic properties (activation, inactivation, deactivation, recovery from inactivation) of an IC can have different temperature dependence (Q 10 value; Lee and Deutsch, 1990). These additional complexities have often been overlooked especially in IC modeling, where a single approximate Q 10 value between 1 and 5 is used to account for temperature dependence. This approximation is used mainly because of the limited number of studies conducted near physiological temperature (last column of Table 2).
Taken together, the non-standardized experimental conditions, the lack of studies near physiological temperature and the unavailability of the electrophysiological traces, prevent reaching a consensus on Kv channel's kinetic properties. We have therefore developed a standardized method to systematically characterize the electrophysiology of the homomeric Kv channels at different temperatures: 15 and 25 • C to compare against the literature, and 35 • C to provide data on the behavior of these channels near physiological temperature. Briefly, we cloned the Kv genes from the rat brain and generated a library of isogenic cell lines over-expressing single homomeric Kv channels, in an inducible manner. We then used automated patch clamping to characterize the biophysics of each channel. This allowed us to construct comprehensive maps of the kinetics of the homomeric Kv channels at 15, 25, and at 35 • C-a temperature close to physiological conditions. The large amount of data produced using automated patch clamping, especially near physiological temperature, allowed us to observe unexpected qualitative changes for specific Kv channels, beside quantitative changes in kinetics associated with changes in temperature. We found that the Q 10 is non-linear not only with temperature but also with voltage. Furthermore, the systematic kinetic characterization revealed that some ICs exhibit inherent kinetic variability (heterogeneity), which might have been overlooked as artifacts. We illustrate how these data on temperature and voltage dependencies in IC kinetics can also be used to build more accurate temperature-dependent models of Kv channels, using the data for the Kv1.1 channel as an example. Contrarily to previous reports (Petersen and Nerbonne, 1999;Fernandez et al., 2003) that show that host cell lines affect Kv kinetic properties, we show that, under standardized conditions and with large sample sizes, the kinetic properties of the homomeric Kv channels are largely consistent across host cell lines (CHO, HEK, CV1). We also show that kinetics is well-conserved across mouse, rat, and human species. The raw and processed data from our study have been made publicly FIGURE 1 | Voltage-gated potassium channels, kinetic characterization workflow. (A) Radial phylogenetic tree of the 145 known voltage-gated ion channel genes, where half are potassium channels (in green). Our study focuses on the 40 voltage-gated potassium (Kv) channels, divided into 12 subfamilies (expanded area). (B) Phylogenetic tree obtained from pairwise alignment of the coding sequences of isoform1 of all rat Kv channel genes; scale represents normalized phylogenetic distance; the corresponding Kv protein names are indicated (prefix r stands for rat). The schematic and scaled depiction of the protein structure of each Kv channel is based on the amino acid length of each domain. (C) Standardized workflow for Kv channels kinetic characterization. Each Kv gene is amplified from a rat brain cDNA library and cloned in a mammalian-expression vector (Table S6), building a library of expression vectors. A library of stable cell lines is generated after transfection of each expression vector in host cell lines (CHO, HEK, or CV1; See also Figures S1, S2 and Table S7). The kinetics of each cell line is characterized using an automated patch clamp setup (APC) in a dedicated temperature-controlled room (14-36 • C; Table S8). Each recorded cell is assigned a unique ID (Cell ID). Electrophysiology data are analyzed, stored and shared through the "Channelpedia" website.

Standardized Kinetic Characterization of Kv Channels
To establish a standardized kinetic map of the Kv family, we developed a screening workflow for the kinetic characterization of Kv channels. The workflow consists of four main steps: cloning, cell line generation, automated electrophysiology, and data sharing ( Figure 1C). Heterologous over-expression of ICs in host cell lines is commonly used to characterize IC kinetics. Previous studies have used many different host cell lines (e.g., CHO, COS7, CV1, HEK, LTK, ND7-23, NG108, NIH3T3; Lalik et al., 1993;Stephens et al., 1997;John et al., 2004). We chose CHO cells as they meet a broad range of criteria for successful IC experiments (Gamper et al., 2005), in particular: an efficient exogenous expression of recombinant proteins, very low endogenous ionic current and a good compatibility with planar automated patch clamp method. To ensure highly standardized conditions, we produced isogenic cell lines using the Flp-In TM system (O'Gorman et al., 1991), which provides reproducible conditions within and across different cell lines. We combined the Flp-In system with the tetracycline induction system T-Rex TM , that enables expression of IC on demand, avoiding possible side effects of constitutive expression (Yao et al., 1998).
We generated a cDNA library from rat brain tissue, and attempted to amplify the 40 known Kv coding genes using the corresponding RefSeqs from the GenBank database (NCBI, Resource Coordinators, 2013; Figure S1). Thirty-five out of 40 Kv genes were successfully amplified from rat brain tissue. The remaining five channels (Kv1.8, Kv6.2,Kv7.4,Kv11.3,and Kv12.1), which we did not succeed to amplify, were obtained by commercial synthesis (see Methods "IC gene cloning"). All the amplified genes were fully sequenced and verified against the original RefSeqs. With minor exceptions, all cloned sequences were identical to the RefSeqs or contained synonymous variations ( Table S1). The amplification process also revealed several new splicing variants ( Table S2). As we only focused on the main isoforms (commonly known as isoform 1), these new variants are not included in this study.
We then cloned each amplified Kv gene in a customized mammalian expression vector and transfected them in CHO Flp-In TM T-Rex TM cells to construct a library of isogenic tetracyclineinducible CHO-Kv stable cell lines ( Figure S1). Each Kv cell line was validated before kinetic characterization, as shown for the CHO rKv1.1 cell line as an example ( Figure S2). Gene expression screening of the entire library confirmed that each cell line over-expressed a single target Kv gene ( Figure S3 and Table S4).

Kinetic Characterization of the Rat Kv1.1 Channel
To illustrate our method, we describe the detailed kinetic characterization of Kv1.1 at 25 • C (Figure 2), a well-studied ion channel that has been implicated in several diseases including episodic ataxia, malignant hyperthermia, and hypomagnesaemia (Rajakulendran et al., 2007;van der Wijst et al., 2010;D'Adamo et al., 2014). Activation properties were analyzed in terms of their I-V relation, time constant for activation (Act τ ), and activation voltage (Act_Volt) (Figure 2A). Normalized peak currents against command voltage (I-V) curves were fitted to a Boltzmann function, yielding V 1/2 = 4.54 ± 5.87 mV, and slope k = 19.84 ± 1.72 mV. Act τ was measured by fitting a single exponential curve to the recorded current trace from start of stimulus to peak current. For voltages from −10 to +80 mV, the value of Act τ decreased from 4.68 ± 1.76 to 0.71 ± 0.15 ms. The median activation voltage (Act_Volt)-defined as the voltage where the channel current exceeds 10% of the peak current-was −31.0 mV (IQR = −34.2 to −27.38).
Deactivation properties were characterized by measuring the tail currents evoked by 200 ms hyperpolarizing stimuli at voltages increasing from −80 to +50 mV in 10 mV steps ( Figure 2B). For voltages between −60 and +10 mV, the deactivation time constant (Deact τ ) increased from 14.96 ± 7.32 ms to 120.68 ± 64.9 ms.
Inactivation properties were analyzed in terms of their I-V relation, time constant (Inact τ ) and inactivation factor (Figure 2C), yielding a V 1/2 value of -35.26 ± 4.0 mV and slope k = 5.62 ± 2.01 mV. For voltages from −20 mV to +70 mV, Inact τ from the peak to the end of the first pulse decreased from 313.25 ± 81.7 ms to 166 ± 48.38 ms. Kv1.1 has been reported to be a non-inactivating or slowly inactivating IC which becomes rapidly inactivating in presence of Kvβ1 subunits (Heinemann et al., 1996;Jow et al., 2004). Our characterization study shows that at +70 mV, this channel inactivates as much as 57.4 ± 10%. However, targeted PCR and full transcriptome screening showed no significant Kvβ1 expression in CHO cells ( Figure S4). This led us to conclude that Kv1.1 is indeed an inactivating IC, and that Kvβ1 is not essential for its inactivation.
Recovery from inactivation was analyzed at a recovery potential of −80 mV using 1.5 s conditioning pulse (+50 mV) to induce inactivation, followed by 150 ms test pulses at varying intervals ( Figure 2D). The time constant for recovery from inactivation (Rec τ ) was measured by fitting a single exponential curve to the maximum current values during the test pulse, yielding median Rec τ = 0.85 s (IQR = 0.72 to 0.95).
To gain more insight into the behavior of Kv1.1, we used two in vivo-like stimuli: a slow voltage ramp and a train of 27 action potentials (APs). With the slow ramp we observed maximum conductance at +12.03 mV (IQR = 1.40 to 18.71; Figure 2E). The response to the train of APs displayed median inactivation of 39% (IQR = 29 to 47), providing further evidence that Kv1.1 is indeed inactivating ( Figure 2F).

Kinetic Characterization of All Rat Kv Channels at 25 and 15 • C
Each of the Kv cell lines was characterized using the method described above for Kv1.1. Figure 3A shows typical behavior of all the 40 Kv channels in response to the activation protocol. This map provides the first comparative overview of all rat Kv channels kinetics obtained in standardized condition.
As a first analysis, in order to differentiate between active and silent channels we measured voltage-dependent activity with activity index (AI) for each cell. AI is calculated by combining signal to noise ratio (SNR) and non-linearity factor (NLF) ( Figure S5). An AI value close to 0 for a cell indicates low voltage dependence denoting an electrically silent channel. Whereas, an AI value close to 1 represents a highly active channel. We produced box plot of AI values for all cells of a Kv cell line and assigned a cut-off value of 0.3 for third quartile (i.e., 75% of the population) to distinguish between active and silent channel. Below this cut-off channels are considered as silent. The remaining channels are considered as active, with a wide range of activity levels. We further used the median values of AI to divide active channels in highly active (median value above 0.3) and low active channels (median value below 0.3).
According to this criterion, at 25 • C all members of the Kv5, Kv6, Kv8, and Kv9 families were classified as silent ( Figure 3B), which is consistent with the literature; these channels are known to have regulatory effects on active channels (Salinas et al., 1997;Kramer et al., 1998;Bocksteins et al., 2014). Several other channels were silent according to the AI criterion: Kv1.8, Kv7.2, Kv7.3, Kv7.4, Kv7.5, Kv11.2, and Kv12.2. Previous studies reported Kv7.2, Kv7.3, Kv7.4, and Kv7.5 channels as active, showing however that the homomeric channels conduct relatively small current, becoming substantially more active as FIGURE 2 | Automated kinetic characterization of CHO rKv1.1 at 25 • C. (A-F) The kinetic properties of rat Kv1.1 channel in CHO cells are extracted from the current response of six different voltage protocols: Activation, Deactivation, Inactivation, Inactivation recovery, and two in vivo-like stimuli, Ramp, and Action potential. Each panel shows the applied voltage stimulus (black traces), representative traces of the current responses (gray traces), region of interest for analysis (red box), and extracted features. For all analysis, current traces are first normalized to overall maximum current (I max ). All data are presented as means ± S.D. (A) I-V curve, activation voltage (Act_Volt), and activation time constant (Act τ ) are extracted for all cells (n = 71). Peak current (I peak ) for each trace is identified and plotted against command voltage to get I-V curve (top right panel). For each cell the voltage where normalized current (I peak ) exceeds 0.1 in I-V curve is considered as the activation voltage (red arrow). Act τ is calculated by fitting single exponential curve from 0 to I peak for each current trace. I peak and Act τ from each cell are plotted and fitted with Boltzmann and single exponential function, respectively to get mean I-V curve and mean Act τ . (B) Deactivation time constant (Deact τ ) for each trace is calculated by fitting a single exponential to current response during the second stimulus pulse (400-600 ms) and then plotted against command voltage (n = 70). (C) Inactivation curve, time constant (Inact τ ), and inactivation factor are features extracted for all cells (n = 67). I peak from the second stimulus pulse (red box) are plotted against command voltage to get Inactivation curve. Inact τ is calculated by fitting a single exponential to each current trace from I peak to the end of the first stimulus pulse. Inactivation factor (x) is the difference from I peak to the end of the first stimulus pulse. Inactivation I-V, Inact τ , and Inactivation factors from all cells are averaged and plotted against command voltage and fitted with Boltzmann, single exponential, and Boltzmann function, respectively. (D) Inactivation recovery time constant (Rec τ ) is obtained by fitting a single exponential function to the peak current values of the responses to recovery pulses (red box) (n = 59). (E) Maximum conductance (V max_Cond ) is calculated on the rising phase of the first Ramp (n = 69). (F) AP-Inactivation is measured by subtracting last AP (AP27) amplitude from normalized maximum value (1) (n = 70). Act_Volt, Rec τ , V max_Cond, and AP-inactivation values for cell population are reported with histograms and box plots.
heteromers (Wang et al., 1998;Nakajo and Kubo, 2008;Gilling et al., 2013), or in association with their accessory subunits KCNE (Barhanin et al., 1996). According to transcriptomic data, KCNE subunits are not expressed in CHO cells (Table S3), that could explain why these channels are silent in this system. Kv11.2 and Kv12.2 have not been formally considered as silent Representative traces of the typical response to activation stimulus for each Kv channel, recorded at 25 • C, is shown. Current traces are sorted by Kv subfamily. ⋆ indicates a channel with inherent kinetic heterogeneity (see also Figure 6); for these channels, only one from a range of recorded responses is shown. (B) Box plots of AI values for each Kv channel at 25 • C; ion channels are ordered by their median AI values and categorized as active or silent based on the 0.3 cut-off value for the 3rd quartile of the box plot (N = 3,409 cells). Active channels are further divided into highly active or low active based on 0.3 cut-off value for the median value of the box plot. channels but have been shown to have very low conductance (Engeland et al., 1998;Wimmers et al., 2002;Sturm et al., 2005), that is consistent with our data. The remaining 23 Kv channels 3) are all electrically active at 25 • C, with wide range of activity levels ( Figure 3B). As control, screening for mRNA expression confirmed that all cell lines, for both active and silent channels, correctly expressed the target gene ( Figure S3 and Table S4). In addition, analysis of membrane fraction by western-blot on a sample of two active and two silent channels further showed that they had been correctly translocated to the membrane ( Figure S6).
The kinetics of all active channels have been already reported ( Table 2), but mainly at room temperature which can vary between 18 and 28 • C from one study to another. These reports contain many inconsistencies, with agreement only for a few channels (Kv1.4, Kv2.2, Kv4.1-4.3, Kv10.1, Kv10.2). The largest inconsistencies concern Kv1.3, Kv1.5, and Kv3.3 (sometimes described as non-inactivating and sometimes as fast inactivating channels), that could be due to differences in experimental conditions across studies (species, host cell lines, method of expression, temperature, intracellular/extracellular solutions, patch clamp technique etc.). The standardized workflow used in our study makes it possible to compare the kinetic properties for all Kv channels, potentially reconciling these inconsistencies.
In our study at 25 • C, most of the active ICs shows more inactivation than previously reported; for example, Kv1.1, Kv1.6, Kv2.1, and Kv2.2, mainly reported as non-inactivating (Table 2), show a consistent inactivating pattern at 25 • C ( Figure 3A). One possible reason could be that the data reported in literature were acquired at temperatures below 25 • C ( Table 2). To assess this possible effect of lower temperature, we recharacterized the complete Kv cell line library at 15 • C. Indeed, all cell lines show less inactivation at 15 • C than at 25 • C (Figures S7A, S8A); for example, Kv1.1, Kv1.6, Kv2.1, and Kv2.2 are almost non-inactivating at 15 • C providing a closer match to the results in literature ( Table 2). This substantial temperature sensitivity points out that measurements of IC kinetics around room temperature do not necessarily reflect their behavior at physiological temperature.
Kinetic Characterization of All Kv Channels at 35 • C Like other biological processes, Kv channel kinetics is known to become faster with increase in temperature (Table 3). However, there are very few studies on Kv kinetics near physiological temperature ( Table 2, right column), and moreover temperature varies between 30 and 39 • C depending on the study. To obtain detailed kinetic behavior for all Kv channels near physiological temperature, we screened the complete Kv cell line library at 35 • C (Figure 4), close to the maximum limit of the temperature controller in our setup. From a technical perspective, it is more difficult to patch cells at 35 • C than at lower temperatures: cell membrane becomes less stable, making it more difficult to achieve and maintain a good seal-likely the reason why so few studies are conducted at higher temperature. In our case, for example, the success rate at 35 • C could be as low as ∼15% compared to ∼80% at 25 and 15 • C. However, the automated patch clamp robot used in our study produced a sufficient number of recordings to compensate for this high failure rate. The map in Figure 4A, provides the first comparative overview of the Kv channel kinetics near physiological temperature (35 • C). The kinetics is qualitatively very different from 15, to 25 to 35 • C (especially for inactivation; see Figure S7 and Figures 3, 4). Quantitatively, Kv channels activate and inactivate faster at 35 • C than at 25 • C (Figures S8B,C). In addition, the comparison of AI values of all channels at 15, 25, and 35 • C ( Figure S9) shows a surprising temperature-dependent effect on the overall level of activity of specific Kv channels. With increases in temperature, Kv10.2, Kv11.1, and Kv11.3 become significantly more active, while Kv4.1 becomes significantly less active (Figure 4C). At 35 • C, Kv7.1 becomes even completely silent; all these specific temperature-dependent effects would require further investigations. The voltage-dependent activity of the other Kv channels is not affected by the temperature; as shown in Figure S9, the channels with high activity stay highly active at all three temperatures. In addition, Kv5, Kv6, Kv8, and Kv9 families stay silent at all temperatures. This observation is consistent with the literature where these channels are described to have regulatory effects on active channels (Salinas et al., 1997;Kramer et al., 1998;Bocksteins et al., 2014). As mentioned before, Kv7.2-Kv7.5, Kv11.2, and Kv12.2 homomeric channels are reported in literature to have low conductance at room temperature; our study shows that they remain silent also at 35 • C in CHO cells. These channels might require heteromerization or modulatory subunits to be active at physiological temperature.

Comparative Kinetic Properties of Kv Channels at 35 • C
Our standardized kinetic characterization of active ICs at 35 • C enables us to provide the first comparative analysis of the kinetic properties of the Kv channels near physiological temperature (Figure 5).
The activation protocol used in our characterization captures the opening of an IC in response to a given voltage stimulus. As illustrated in Figure 5A, at 35 • C, all active Kv channels show a range of current responses from 0.67 to 18.35 nA (median values). The ICs with the largest current responses (mainly Kv3 family) activate at the highest voltages ( Figure 5B) and have the fastest activation rate (Figure 5C, see also Figure S10 in blue). This is consistent with reports of neuronal AP repolarization by currents mediated by Kv3 family (Rudy and McBain, 2001;Labro et al., 2015). In contrast, ICs of the Kv11 and Kv12 families, that show low current responses, activate at lower voltages and have the slowest activation rate (Figures 5A-C and Figure S10 in red). These ICs may play a role in repolarizing the cell membrane at sub-threshold voltages and during sub-threshold oscillations. The behavior of the other Kv ICs lies between these extremes. In general, we observe that channels with high current responses tend to be activated by high voltages and vice-versa (Figure S10 in black).
The Kv4 family is an exception, with all members activating at high voltage but having a low current response (Figure S10 in green).
Deactivation refers to the closing of a channel when a stimulus is removed. The lowest deactivation time constants at 35 • C are found for the Kv3 and Kv10 families. The ICs with the highest deactivation time constants (up to 60 ms) belong to Kv11 and Kv12 families (Figure 5D). Inactivation refers to the closing of a channel during stimulation. Inactivation is an important property of most of the Kv channels, determining how long the channel affects the cell, under suitable activation conditions. There are several conflicting reports about the inactivation properties of Kv channels ( Table 2). Our study shows that their kinetic behavior span from highly inactivating (e.g., Kv1.4, Kv1.3, Kv3.4, Kv1.7) to non-inactivating (e.g., Kv12.1, Kv10.1, Kv12.3; Figure 5E). In some cases, we confirm previous reports. For example, we confirm that Kv1.3, Kv1.4, and Kv.4.3 inactivates more than 80% at 35 • C. We also confirm that Kv10.1, Kv10.2, Kv12.1 are truly non-inactivating ICs, maintaining this characteristic at 35 • C ( Figure 5E). However, in other cases, the extrapolation from 25 to 35 • C does not hold. For example, at 35 • C, Kv1.6, Kv1.1, and Kv1.2, which are non-inactivating at 15 • C, and slowly inactivating at 25 • C, actually inactivate by 83, 74, and 70%, respectively at 35 • C. Other Kv channels show similar trends. For example, Kv2.1 and Kv2.2, which have been reported in literature to be less inactivating, actually inactivate by 65% and 75% respectively, at 35 • C.
We also observed inactivation in response to AP-like physiological stimuli. When applying a train of APs, inactivating channels like Kv1.4 and Kv1.7 show up to 80% inactivation, whereas non-inactivating channels like Kv10.1 and Kv12.1 show no inactivation at all (Figure 5G).

Inherent Kinetic Heterogeneity and Delayed Inactivation
The effect of temperature on inactivation explains some inconsistencies in ion channel literature. For example, Kv1.1, Kv1.6, and Kv2.1 are reported as non-inactivating or slowinactivating depending on the study (Table 2); that we show is a temperature-dependent effect as these channels are all non-inactivating at 15 • C but become inactivating at 25 • C ( Figure S7 and Figure 3). This systematic characterization led to an additional observation: in our study, the majority of the Kv cell lines shows responses that are consistent across different cells from the same cell line. However, in some cases, different cells from the same cell line systematically show different responses (kinetic heterogeneity), despite the fact that they are isogenic (Figures 6A,B). To quantify this kinetic heterogeneity, we have measured the variance in current response for all cells of each cell line, in response to +80 mV command stimulus from activation protocol ( Figure 6C). This analysis shows that at 25 • C, four ICs, namely Kv1.3, Kv1.5, Kv3.3, Kv3.4, exhibit particularly high kinetic heterogeneity. This phenomenon was observed independently of the temperature or the host cell line (data not shown).
We have verified that the heterogeneity observed is not due to gene mutations (Table S1), cell contamination (Figure S3), recording artifacts (series resistance, intracellular/extracellular composition, extracellular K + accumulation, etc.), cell handling or cell environment. Since the methods used were highly consistent (standardized workflow, isogenic cell lines, recording of ∼100 cells for each of these channels), we conclude that the kinetic heterogeneity is an inherent property of specific Kv channels. The underlying mechanism may include transcriptional, post-translational and/or translocational factors and requires further investigation. Except for one previous report about two distinct phenotypes for Kv1.2 cells (Rezazadeh et al., 2007), this is the first report of inherent kinetic heterogeneity in Kv channels. The inherent heterogeneity of Kv1.3, Kv1.5, and Kv3.3 could also explain inconsistencies in previously reported kinetics for these channels. However, the high heterogeneity we observed in the kinetics of Kv3.4, contradicts all previous reports, which describe this channel as a fast inactivating channel ( Table 2).
Three of the abovementioned channels (Kv1.3, Kv3.3, Kv3.4) also show striking delayed inactivation, with inactivation starting after a delay that varies from cell to cell-a behavior that differs from classical inactivation patterns (slow, intermediate, or fast inactivation; Figures S11A,B). With longer stimuli, delayed inactivation also appears in two further channels (Kv3.1, Kv3.2) (data not shown). The delayed inactivation was observed at all three temperatures ( Figure S11C). To assess if the delayed inactivation is linked with extracellular potassium accumulation, we performed control experiment on Kv3.4 cells. A single voltage pulse of +70 mV was applied to evoke delayed inactivation and then extracellular space was washed with extracellular FIGURE 4 | with inherent kinetic heterogeneity (see also Figure 6); for these channels only one response from the range of recorded responses is shown. (B) Box plots of AI values for each ion channel at 35 • C; ion channels are ordered by their median AI values and categorized as active or silent based on the 0.3 cut-off value for the 3rd quartile of the box plot (N = 2,350 cells). Active channels are further divided into highly active or low active based on 0.3 cut-off value for the median value of the box plot. (C) AI values at 15 • C (blue), 25 • C (black), and 35 • C (red) are plotted for Kv channels that show significant change in activity over temperature. AI values for the non-transfected CHO-FT cell line is plotted as a control (see also Figure S9). ***p < 0.001, Student's t-test.
solution for 90 s to remove potassium accumulation. The delayed inactivation was observed even after multiple washes with extracellular solution (Figure S12A). Moreover, for Kv3.4 cells, Person's linear correlation coefficient between maximum current amplitude and delay in inactivation resulted in a value of 0.12, indicating that the delayed inactivation is not correlated with outward potassium current amplitude (Figures S12B1-B3).
Full transcriptome analysis showed no significant expression of Ih, Ca or Na voltage-gated channels, indicating that this delay was not caused by inward current (Figure S12C and Table S3). The mechanism behind this delayed inactivation requires further investigations.

Voltage-Dependent Q 10
The Q 10 has been widely used in biology as a convenient measure of temperature effects. For ion channel modeling, Hodgkin and Huxley (1952) used a constant value of Q 10 = 3 to model Na + and K + conductances in the squid axon. Since then, other authors have used a fixed Q 10 value between 1 and 5 in models to estimate the kinetic rates for Kv channels at different temperatures. Our experiments at 15, 25, and 35 • C, allowed us to measure empirical Q 10 values for all kinetic parameters, using Van't Hoff 's equation: where T and T1 are the temperatures (in • C) at which experiments were carried out, and τ (T) and τ (T1) are the time constants.
As an example, Figure 7 shows Q 10 for Act τ for selected Kv channels with very different activation profiles. Act τ was obtained by fitting the current traces to a single exponential and plotted against command voltage ( Figure 7A). Q 10 for 15 and 35 • C were calculated considering 25 • C as the base temperature, the usual temperature used for ion channel studies: Applying Equation (2) to our empirical data for all active Kv channels yields values that vary not only as a function of temperature but also of voltage ( Figure 7B). For example, Kv1.1 has a Q 10 15c value around 5 and a Q 10 35c value between 3 and 4, varying with voltages. Data for other IC show that Q 10 can be lower than 2 (ex: Kv7.1 at 15 • C) or as high as 8 (e.g., Kv2.1 at 15 • C) depending on voltages and temperature range. Our detailed kinetic characterization of all Kv channels provides data to calculate accurate temperature and voltage-dependent Q 10 for all kinetic parameters.

Voltage-Dependent Q 10 Hodgkin-Huxley Model
One case where extrapolating across temperatures may lead to inaccurate results is the implementation of IC kinetics in computational models. Neuronal models often use the Hodgkin-Huxley (H-H) formalism to describe IC conductances. In this section, we describe a revised model-fitting procedure, which accurately incorporates voltage dependence in temperature sensitivity. We use Kv1.1 as an example. Briefly, current traces from individual cells were fitted to single activation and inactivation gates in the H-H formulation, and model parameters were extracted for each cell. The medians of these parameters were fitted with Q 10 functions.
The net transmembrane current flow, I Kv1.1 , and the conductance g Kv1.1 were computed using Equations (3) and (4): Where V m is the membrane potential, and E k denotes the Nernst potential. Variable g scales for maximum channel conductance. m and h represent the proportion of open activation and inactivation gates, respectively; p and q are the numbers of independent gates required to account for the observed time course of activation and inactivation. Since single gates are used both for activation and for inactivation, p = q = 1. Gating variables are modeled as a first-order kinetic process Where m ∞ , h ∞ represent the voltage-dependent steady state and m τ , h τ are the voltage-dependent time constants for activation and inactivation gates. The free parameters in Equations (3)-(6) were fitted to the normalized conductance for each cell (Figure 8A for a sample fit). The fitted parameter values for all cells were plotted and median parameter values for each temperature were calculated ( Figure 8B).
The m ∞ function was considered temperature independent, thus a single Boltzmann function was used for all three temperatures. The remaining three parameters, h ∞ , m τ and h τ , were fitted with three different Q 10 functions. For the steady state inactivation function h ∞ , we used a temperature-dependent linear function h ∞ Q 10 . m τ curves were fitted using a double exponential function of voltage and temperature, mτ Q 10 . The hτ Q 10 function was represented by a constant value of  Figure 2, are plotted for the 22 active Kv channels at 35 • C (see Figure 4B). Features are represented as box plots for each ion channel and sorted by their median values. N is the total number of cells used for the analysis of each feature (see also    Figure S11). (A) The evoked current trace corresponding to the activation stimulus at +80 mV is normalized to maximum current for each cell. The overlay plot presents the normalized currents from 25 cells that statistically represent the whole group (same variance). (B) Overlay plots as described in panel A are shown for all active channels at 25 • C (Kv11.1 and Kv11.3 are not included due to low current and high noise). The overlay plots visually represent the kinetic heterogeneity. (C) Kinetic heterogeneity is quantified by calculating Var score , obtained by summing five variances of normalized currents evaluated at five time points for all cells of a given cell line, as shown for Kv12.1 as an example. The channels that show the highest heterogeneity (highest Var score ) are Kv1.3,Kv1.5,Kv3.3,and Kv3.4 Table S5. Figure 8D). Figure 8C shows the final plots for each gating variable, and the normalized conductance plot for each temperature.

(
The main difference between our approach and regular H-H modeling was that, instead of using a constant value for Q 10 (usually between 2 and 3), we used a voltage and temperaturedependent function for m τ and a temperature-dependent linear function for h ∞ . With these changes, we were able to fit the temperature sensitivity of Kv1.1 to the experimental data, which was not possible with the standard model. The resulting revised H-H model can thus be used to extrapolate kinetics from the temperature at which the channel was studied to physiological temperatures.

Kv Channel Kinetics Across Host Cell Lines and Species
Previous studies of IC kinetics used many different host cell lines, with different origins, morphologies, intracellular environments, cell handling etc. These differences may also have contributed to all those conflicting reports on the kinetics of specific ICs. To verify possible dependencies on the specific cell line (i.e., CHO) used in our study, we used our standardized workflow to characterize the kinetics of five significantly different active Kv channels (Kv1.1, Kv1.4, Kv1.5, Kv1.6, and Kv2.1) in two additional cell lines (CV1 and HEK). The three cell lines are all adherents and have different morphologies ( Figure 9A). In contrast with CHO cells, CV1 and HEK cells both generate small outward currents, possibly due to higher levels of endogenous IC expression (Table S3). Our analysis of three kinetic features (I-V, Act τ , Inactivation factor) shows that, except for Kv1.5's inactivation factor, the kinetic behaviors of the ICs are largely consistent across the three cell lines at 25 • C (Figures 9B,C). Furthermore, the kinetics are comparable also at 35 • C for Kv1.1, Kv1.4, and Kv2.1 (Figures S13A,B). IC genes are highly conserved across species (Siepel et al., 2005); 36 out of 40 rat Kv channels have >97% homology with mouse orthologs ( Table 4) and even the least homologous channels (Kv6.4, Kv8.2) share >95% homology in their protein sequence. Similarly, 23 out of 40 rat Kv channels show >95% homology with human orthologs, and 36 out of 40 share >90% homology. However, some channels do show differences in the length of their N-terminus and/or C-terminus regions ( Table 4). What is not known is whether and how these small differences in sequences translate into a difference in Kv channels kinetics. To address this question, we cloned the mouse and human versions of the same five active channels previously used for host cells comparison, and we characterized them in CHO cells at 25 • C (Figures 9D,E) and at 35 • C for Kv1.1, Kv1.4, and Kv2.1  also Table 4). (E) Median values for three kinetic features (I-V curve, activation time constant, and inactivation factor calculated as explained in Figure 2) obtained from rat, mouse, and human genes of each Kv channel, overlaid for comparison. Error bars are ± S.D. The amplitude in nA and time in ms are indicated with scale bars. ⋆: For Kv1.5, only one representative trace from a range of responses is shown. The number of cells (N) for each group is listed in Table S5.
( Figures S13C,D). Among these five characterized channels, Kv1.1, Kv1.4, Kv1.6, and Kv2.1 have high homology between species, whereas Kv1.5 was chosen as an example of an active channel with one of the lowest homology across species. We detected only subtle quantitative differences in some kinetic parameters across species both at 25 and 35 • C (Figures 9D,E and  Figures S13C,D). For example at 25 • C, we measured a slightly slower activation time constant for human Kv2.1 compared to rat Kv2.1, as previously reported (Ju et al., 2003; Figure 9E); however, the difference no longer remains significant at 35 • C ( Figure S13D). We conclude that the behavior of these five Kv channels is relatively conserved across the three species.

A Data Resource of Ion Channel Kinetics
To facilitate sharing and further validation of our results, we have published all raw and processed data from our study as a wiki-like web resource (https://channelpedia.net or https:// channelpedia.epfl.ch), freely available for academic use. Each cell has been assigned a unique cell ID, which is appended to the data file name. The raw data file provides the recorded currents for each repetition of the voltage protocols and a description of the stimuli applied. The analyzed data file contains features extracted from recordings of each repetition of each protocol. The experimental and analyzed data for each Kv IC are made available in Neurodata without Borders (NWB) format under "Experimental Data" section on the Channelpedia website. NWB is a new format designed to promote data standardization and sharing (Teeters et al., 2015). NWB is based on the HDF5 file format, thus can be viewed with any HDF5 viewer. The raw data, stimulus and metadata are stored under "acquisition/timeseries, " "stimulus/presentation, " and "general, " respectively. The data can be accessed via a web interface or a REST API. Over 1 million recording traces from over 18,700 cells are available and can be downloaded for academic use.

DISCUSSION
In this study, we produced a standardized map of the kinetics of the homomeric form of all members of the Kv family, at 15, 25, and 35 • C. The map includes all the channels, also those not previously studied or underrepresented in the literature ( Table 1). The data we provide near physiological temperature were not available for most of the Kv family members. We show that under controlled conditions, these channels kinetics is reproducible and mostly conserved across host cell lines and species. We find that for many Kv channels, the kinetics near physiological temperature not only differ quantitatively as expected from temperature-dependent effects, but also qualitatively. We find that for all Kv channels, Q 10 is not only temperature-dependent, but also voltagedependent. To allow extrapolation from the large number of low temperature kinetics available in the literature to kinetics at physiological temperatures, we provide the data on the voltage-dependent Q 10 for each channel and a revised Hodgkin-Huxley model to capture the channel kinetics accurately. We found a novel form of delayed inactivation and inherent heterogeneity in kinetics for some Kv members, behavior that was likely interpreted as artifacts in previous studies. We suggest that different experimental conditions, non-trivial temperature dependencies and inherent heterogeneity can account for most of the inconsistencies in kinetics reported in the ion channel literature. The full dataset of over one million traces is available as a publicly accessible resource (https://channelpedia.net or https://channelpedia.epfl.ch).

Reference Kinetic Dataset for Kv Channels
A vast amount of data on Kv channel kinetics is available in the literature. The lack of consistency limits the value of these data and prevents a consensus on ion channel kinetics. The reference dataset we provide was generated using a standardized heterologous system with isogenic cell lines and an inducibleexpression system to achieve a high level of reproducibility and reliability.
Our system, however, does have limitations. In particular, levels of expression are probably much higher than in vivo, and it is not possible to control levels of expression across cells. We also cannot account for possible effects of transcriptional regulation in different cells of the body at different stages of development.
The need of special recording solutions for APC experiments is another limitation of this study. Specific fluoride-based intracellular solution (ICS) (50 mM KCl, 10 mM NaCl, 20 mM EGTA, 60 mM KF, 10 mM Hepes), extracellular solution (ECS) (140 mM NaCl, 4 mM KCl, 1 mM MgCl 2 , 2 mM CaCl 2 , 5 mM Dglucose monohydrate, 10 mM Hepes) and seal enhancer solution (SES) (80 mM NaCl, 3 mM KCl, 10 mM MgCl 2 , 35 mM CaCl 2 , 10 mM Hepes) were used for APC experiments. In ICS, instead of usual 2-10 mM, 20 mM EGTA was used as per specifications of APC robot to keep cells longer in whole cell configuration and to maintain constant access resistance. However, we verified that reducing EGTA to 10 or 0 mM did not affect the kinetics of the selected Kv ion channels (data not shown). Inclusion of fluoride in ICS has been shown to improve patch clamp seal quality and to stabilize the cell membrane, resulting in longer and more stable patch clamp recordings (Kostyuk et al., 1975). While the effect of internal fluoride on Kv channel kinetics has been reported, directly (Adams and Oxford, 1983) or through PIP2 (Rodriguez-Menchaca et al., 2012), there are other reports claiming a lack of side effects on Kv channel kinetics (Clay,  1988; López-López et al., 1993;Zeng et al., 2013). Along with fluoride in ICS, SES with high Ca 2+ and Mg 2+ was also applied extracellularly to achieve the giga seal. In our experience, without fluoride in ICS and high Ca 2+ extracellularly, the APC robot could not form the giga seal. The SES was washed out with ECS after reaching whole cell configuration as well as between each repetition of voltage protocols. We observed that high Ca 2+ and Mg 2+ in SES did cause membrane potential oscillation artifacts for Kv11.1 and Kv11.3 cell lines, but it disappeared after washing with ECS solution (Figure S14). In conclusion, fluoride in ICS and high Ca 2+ in SES were necessary to achieve giga seal in our APC experiments. Previous studies have found no significant differences in channel kinetics between manual and automated patch clamp system (Li et al., 2017). In addition, we did not find any difference in Kv2.1 kinetics when comparing with data obtained with another automated patch clamp robot (Fejtl et al., 2007) that worked with physiological ICS, ECS, and without the need of SES solution (data not shown). Despite these limitations, the solutions used in this study still contain physiological levels of sodium and potassium concentration to better reflect in vivo ion composition, and not bi-ionic high potassium (K + ) solution as used in many previous studies on potassium channels. Currently, the reference dataset we provide, contains the main isoform of each Kv channel and does not include other known splicing variants. It is known that Kv channel genes, such as those involved in the immune system, are highly prone to alternative splicing (Lipscombe, 2005) that substantially increases the repertoire of Kv channels proteins. For example, the 40 known human Kv genes have 73 validated, and 113 predicted splicing isoforms. Our study additionally identified several unreported splicing variants ( Table S2), showing that the full repertoire is yet to be mapped. The investigation of the kinetic activity of the different variants using our standardized system would also be of high value.

Kv Kinetics Near Physiological Temperature
Data on Kv kinetics near physiological temperature were not available for most of the Kv members ( Table 2). This is the first complete map of Kv channels kinetics at 35 • C. The qualitative differences we found at 35 • C clearly suggests the need of a reassessment of previous functional characterizations of ICs. For example, Kv1.1 is generally considered to be a non-inactivating channel, which only inactivates in presence of Kvß1 or Kvß3 subunits ( Table 2). We show that Kv1.1 actually inactivates strongly at 35 • C, even in absence of these subunits in CHO cells (Figure S4). Similarly, Kv1.2 and Kv1.6, reported to become fast inactivating in presence of Kvß1 or Kvß3 subunits (Heinemann et al., 1996;Bähring et al., 2004), actually are strongly inactivating at 35 • C even in absence of these subunits. As another example, we show that Kv7.1, one of the most studied ion channel, reported as active at room temperature and also at 35 • C (Loussouarn et al., 1997), is actually only active at low temperatures and becomes silent at 35 • C in CHO cells. Furthermore, it is generally accepted that KCNE1 subunit (mink) is required for proper functioning of Kv7.1 (Sanguinetti et al., 1996;Aromolaran et al., 2014); the KCNE1 subunit might make Kv7.1 more active, but it is not required to form a functional channel, at least at low temperature. These examples illustrate how understanding the homomeric channel kinetics across temperatures is essential to correctly interpreting the effects of the several levels of regulation, such as heteromerization, accessory subunit associations or effects of post translational modifications (e.g., phosphorylation, glycosylation).

Implications for Ion Channel Modeling
The biophysical map and the reference dataset provide data for both improving the accuracy of genetically specified ion channel models and for cellular models (e.g., neurons and heart cells) which often use generic conductances (e.g., A-type channels, delayed rectifier channels) that do not capture the diversity of kinetics present within these generic classes. For example, Kv1.1 and Kv1.7 both belong to the Kv1 family but have completely different time constants for activation and inactivation. Generic IC models are unable to exploit the growing volume of data from single cell transcriptomics, which makes it possible to identify the full spectrum of IC genes expressed by specific cell types.
The data we publish will enable researchers to develop improved models for all the Kv channels. Particularly important is the availability of temperature and voltage-dependent Q 10 for all kinetic parameters of Kv channels. Patch clamp experiments on cell lines at physiological temperature are difficult and hence majority of the experiments in the future might continue to be performed at lower temperatures. The availability of Q 10 values from our data makes it possible to interpret those data in the context of physiological temperature. This prospect is very important for ongoing brain modeling programs such as the Blue Brain Project and the Allen Institute for Brain Sciences as well as for initiatives which develop on their work, such as the Human Brain Project. In addition, since specific receptors affect different subsets of ICs in different cells, the new data open the prospect of modeling the effects of neuromodulators on specific neuron types (Nicoll, 1988;Khorkova and Golowasch, 2007).

Inherent Heterogeneity and Delayed Inactivation
Our standardized kinetic characterization has revealed inherent kinetic heterogeneity in four Kv channels, namely Kv1.3, Kv1.5, Kv3.3, and Kv3.4, three of which also show a novel form of delayed inactivation. This inherent heterogeneity is present in isogenic cell lines that have not been manipulated in any way. The large number of recordings in standardized condition and screening of all Kv channels, allowed us to conclude that it is a truly inherent property. To our knowledge, the only previous study that mentions inherent heterogeneity is a report of two distinct phenotypes for rat Kv1.2 expressed in CHO cells after transient transfection (Rezazadeh et al., 2007). Specifically, the authors describe two gating phenotypes (fast and slow) influenced by a cytoplasmic regulator. Modulation in IC kinetics have been reported to occur upon glycosylation (Brooks et al., 2006;Noma et al., 2009), phosphorylation (Desai et al., 2008;Ritter et al., 2012), oxidation (Hoshi and Heinemann, 2001), or by co-expression of beta-subunits (Heinemann et al., 1996;Bähring et al., 2004) or other binding proteins (Pongs and Schwarz, 2010). We suggest that especially these four channels (Kv1.3, Kv1.5,Kv3.3,and Kv3.4) are susceptible to abovementioned modifications. The inherent heterogeneity we described may also explain some of the conflicting reports present in the literature ( Table 2).
The delayed inactivation has been observed only for human Kv3.1, and only at high (>35 • C) temperatures (Oliver et al., 2017). The mechanism underlying this novel form of delayed inactivation is unknown and would require further investigation.

Conserved Kinetics Across Species
Using a sample set of five Kv channels with different kinetics and different degree of homology, we found only subtle differences in some kinetics parameters between mouse, rat and human channel kinetics. Compared to the strong temperaturedependent effect on ion channel kinetics, the inter-species differences we observed are minor. According to these data on a sample set and the high degree of homology of ion channels (Table 4), we postulate that data on ion channels kinetics could be generalized over species. However, this has to be done with caution as even subtle differences observed in kinetics for a single ion channel type could have a strong impact over a neuronal network.
Our observation is consistent with the fact that transmembrane regions S4 to S6, believed to be primarily responsible for channel kinetics, are structurally almost identical across the three species. Nearly all the interspecies differences in the Kv channels sequences concern the N-terminus and C-terminus regions and the S1-S2 loop. The N and C-terminus contain protein binding motifs that are involved in the distribution of IC through the cell (Lai and Jan, 2006;Vacher et al., 2008;Duménieu et al., 2017). The N-terminus region is also involved in tetramerization and N-type inactivation via the "ball and chain" mechanism (Bezanilla and Armstrong, 1977). Hence significant differences in the length of N-terminus or C-terminus regions could affect channels kinetics between different species. For example, human Kv1.7, that has a 32 amino acids shorter N-terminus, has been reported as slowinactivating while the longer versions in mouse and rat are reported as fast inactivating (Finol-Urdaneta et al., 2006. Such channels with significant difference in the length of their N-terminus and/or C-terminus regions ( Table 4) would require separate studies to clarify the effect of such possible interspecies differences on their kinetics.

Implications for Drug Discovery
The methods we have established to create the kinetic map of Kv channels can be used to systematically screen drug candidates for potentially positive or deleterious effects on their kinetics. They can also be used to test the effects of targeted mutations and the ability of drugs to reverse their effects. Such studies have the potential to improve our understanding of disease mechanisms, and have obvious applications in drug discovery. In a longer-term perspective, detailed neuron models incorporating genetically specified ion channels could allow simulation of channelopathies at the cellular and network level.

Future Outlook
For the first time in ion channel research, electrophysiology data from a complete family of ICs have been obtained and made publicly available. We have developed Channelpedia as an online tool for data sharing and we believe that the availability and the ease of data access on this website may encourage other researchers to make their electrophysiology data available through Channelpedia or similar platforms. The availability of the new dataset allows us to investigate links between IC genetics, structure, and biophysics. For example, preliminary results from pairwise alignment of the sequences for the S1-S2 loops, reveal a genetic basis for at least two big groups-inactivating and non-inactivating ICs (data not shown). Large-scale bioinformatics studies aimed at establishing genotype-phenotype correlations are now possible. A key issue for future works is to unravel the mechanism for inherent heterogeneity of IC kinetics.
Homotetrameric Kv channels consist of four identical alpha (α)-subunits, but additional diversity is introduced when different α-subunits combine to form heterotetramers (Dodson et al., 2002;Hadley et al., 2003;Zhang et al., 2016). The function of homotetramers and heterotetramers is further modulated by a broad range of auxiliary (signaling/scaffolding) proteins (Levitan, 2006;Li et al., 2006), including Kvβ-subunits (Pongs and Schwarz, 2010), chaperone proteins (HSP90, HSP70) (Ficker et al., 2003), and other regulatory protein families like K + channel-Interacting proteins (KChIPs), K + channel-associated proteins (KChAP), 14-3-3 proteins, A-kinase-anchoring proteins (AKAP) (Zhang et al., 2016), Dipeptidyl aminopeptidase-like protein (DPP), and KCNE-encoded proteins (Coetzee et al., 1999). Among all possible combinations for heterotetramers and auxiliary subunits, some have been shown to have modulatory effects on kinetic properties, trafficking, permeation, and stability of Kv channels. However, there are also contradictory results on the role of auxiliary subunits; for example on the role of Kvβ3 on Kv1.1 activity ((Bähring et al., 2004) shows that coexpression of Kvβ3 modulates the inactivation of Kv1.1 whereas (Majumder et al., 1995) shows no effect of Kvβ3 on Kv1.1 kinetics) or on the effect of KCNE1(Mink) on Kv7.1 activity ((Bett et al., 2006) shows that Mink modulates KCNQ1 kinetic whereas (Sanguinetti et al., 1996) shows that Mink is necessary for KCNQ1 to be active). The next layers of the kinetic map for the Kv family are the heterotetramers map (kinetics of viable combinations of α-subunits), the β-subunits map (kinetics when bound to β-subunits), and the neuromodulatory map (kinetics for all forms of modulation). The biggest challenge in building these maps is to identify principles for biologically viable combinations of subunits and modulation since all theoretical combinations cannot be mapped. Knowledge of the detailed kinetics of the homotetrameric channels with their structural and genetic correlates, combined with transcriptomics data on the expression of channels in single cells or cell types, will provide a starting point for investigation of this vast, as yet unexplored, combinatorial space.

ETHICS STATEMENT
All procedures were conducted in conformity with the Swiss Welfare Act and the Swiss National Institutional Guidelines on Animal Experimentation for the ethical use of animals.

AUTHOR CONTRIBUTIONS
RR, EL, and HM: conceptualization, writing, and supervision. RR, EL, and MM: methodology and data curation. RR: data analysis. RR and ES: coding and channelpedia. RR, EL, MM, MH, VT and VB: investigation.

FUNDING
This work was funded by the EPFL Blue Brain Project Fund, ETH board funding to the Blue Brain Project and the funding from the EPFL to the Laboratory of Neural Microcircuitry (LNMC).