Outward Rectification of Voltage-Gated K+ Channels Evolved at Least Twice in Life History

Voltage-gated potassium (K+) channels are present in all living systems. Despite high structural similarities in the transmembrane domains (TMD), this K+ channel type segregates into at least two main functional categories—hyperpolarization-activated, inward-rectifying (Kin) and depolarization-activated, outward-rectifying (Kout) channels. Voltage-gated K+ channels sense the membrane voltage via a voltage-sensing domain that is connected to the conduction pathway of the channel. It has been shown that the voltage-sensing mechanism is the same in Kin and Kout channels, but its performance results in opposite pore conformations. It is not known how the different coupling of voltage-sensor and pore is implemented. Here, we studied sequence and structural data of voltage-gated K+ channels from animals and plants with emphasis on the property of opposite rectification. We identified structural hotspots that alone allow already the distinction between Kin and Kout channels. Among them is a loop between TMD S5 and the pore that is very short in animal Kout, longer in plant and animal Kin and the longest in plant Kout channels. In combination with further structural and phylogenetic analyses this finding suggests that outward-rectification evolved twice and independently in the animal and plant kingdom.


Introduction
Voltage-gated potassium (K + ) channels have been investigated in deep detail in various organisms ranging from pro-to eukaryotic species (reviewed in [1][2][3][4]).Having a common subunit structure, these channels consist of six transmembrane domains (TMD), S1 to S6, and a pore helix and selectivity filter between the last two TMD (S5-P-S6) (Fig 1A).Functional channels are tetramers, in which the four subunits are twisted with each other.The S5-P-S6 parts act together in the centre of the channel and form the ion conduction pathway.TMDs S1-S4 form the voltage sensors that are located in the periphery of the channel.They are connected to the pore via a short amino acid sequence, the S4-S5 linker.
Our current picture of channel gating assumes that the fourth TMD (S4) moves in a voltage-dependent manner within the membrane (Fig 1B,top).This movement is forwarded to the pore and results in conformational changes that cause the pore to open or close [5].Initially, the S4-S5 linker was supposed to transmit the S4 movement to the pore and to couple both channel modules [6,7].Recent evidence, however, indicates that this picture is probably too simple.Instead, the interaction of the S4-S5 linker with the C-terminal end of S6 has been described to be crucial for gating of various voltage-gated potassium channels [8][9][10].Additionally, a second connection between the voltage sensor and the pore has been discovered on the extracellular channel side with particular importance for channel functionality [11].In contrast to the intracellular interface (S4-S5 linker) that connects voltage sensor and pore on the primary protein structure level, the extracellular interface is established on the tertiary protein structure level by interacting amino acid side chains.The function of the extracellular interface is described as a fix point that is necessary to enable the intracellular interface to act as a mechanical lever on the pore's gate.During channel gating the S4 movement shifts the mechanical lever, while the extracellular interface provides the structural stability within the channel to allow an S4 movement relative to the pore.
Voltage-gated potassium channels appear in at least two different functional typeshyperpolarization-activated, inward-rectifying (K in ) and depolarization-activated, outwardrectifying (K out ) K + channels.For practical reasons, we use here in this study the unified nomenclature of K in and K out channels that is well established for voltage-gated K + channels of plants [12].We should emphasize, however, that inward-rectifying K in channels must not be confused with two transmembrane domain Kir channels (also called 'inward-rectifying K + channels'), which gain their rectification properties by a voltage-dependent block and not by the movement of a voltage sensor [13,14].
In the animal kingdom, we identified two channel families that can be classified as K in and K out channels, respectively.The Kv family represents animal K out channels and the HCN family Each subunit (grey representation) comprises six transmembrane domains (S1-S6) and a selectivity filter between S5 and S6.The first four TMDs form the voltage sensor that is located in the periphery of the functional channel.The last two TMDs contribute to the channel pore in the centre of the channel.(B) TMD S4 of the voltage sensor translocates in response to the transmembrane voltage within the membrane and affects the channel conformation.It moves towards the intracellular side of the membrane upon hyperpolarizing membrane voltages and towards the extracellular side upon membrane depolarization.The "down" movement of S4 causes K in channels to open, while the same S4 movement results in closure of K out channels.The opening of K out and closure of K in channels proceeds in the inverse manner in response to S4 "up" movements.doi:10.1371/journal.pone.0137600.g001might be considered as animal K in channels.Strictly speaking, HCN channels are cation channels that possess also a significant permeability to sodium (Na + ).Nevertheless, potassium permeates to a much larger extent, and in the absence of K + HCN channels hardly transport Na + [1,15].Besides, HCN channels show the same six TMD structure, including the voltage sensor, as animal Kv and plant voltage-gated K + channels, and contain the selectivity filter motif GYG that is present in highly selective potassium channels.Because HCN channels are furthermore voltage-gated and hyperpolarization activated as it has been described for plant K in channels it appears to be appropriate to include this channel class in our comparative analyses.For the sake of simplicity, we refer to HCN channels as animal K in channels.
In plants the family of voltage-gated K + channels combines both, K in and K out channels, together with two other functional types: weak-rectifying (K weak ) K + channels and functionally silent regulatory subunits [3].The first voltage-gated K + channels in plants had been described to be structurally similar to animal Kv channels [16][17][18].Therefore, in literature the term "plant Kv-like channels" has widely been established when referring to the group of plant voltage-gated K + channels.So we will use this term in the following as well.
The response of voltage-gated K + channels to a change in the transmembrane voltage begins with the transition of the fourth TMD.Interestingly, S4's movement is conserved in these K + channels.It has been shown experimentally that S4 moves downward upon hyperpolarizing voltages and upward upon depolarizing membrane voltages [19].Curiously, the same S4 movement results in opening of K in channels, while it results in the closure of K out channels and vice versa (Fig 1B, bottom) [20,21].Because the fundamental voltage sensing mechanism is not different, it is hypothesized that the coupling of voltage sensor and pore must be different in both channel types.Astonishingly, the main structural components that are important for channel gating, like voltage sensor, S4-S5 linker and pore, are quite conserved in voltage-gated K + channels and bear only slight differences [12].Thus, coupling of voltage sensor and pore as well as the opposite rectification of K in and K out channels seem to be subtly implemented within the channel structure.
In this study we analysed and compared sequence and structural data of voltage-gated K + channels with a focus on the four groups of animal K in (HCN), animal K out (Kv), plant K in and plant K out (plant Kv-like) channels.The results allowed us (i) to reason an independent evolution of outward-rectification in plants and animals and (ii) to develop a more sophisticated picture on voltage-gated K + channel gating that might be suited to account for differences in gating of these channels.

Sequence data analyses
Known sequences of voltage-gated K + channels were extracted from databases indicated in tables and figures.Putative sequences of K + channels from Klebsormidium flaccidum were extracted from http://www.plantmorphogenesis.bio.titech.ac.jp/~algae_genome_project/ klebsormidium/index.html [22].Sequences of the plant K in channel KAT1 and the plant K out channel SKOR were blasted against the database.Those sequences have been selected that contained features characteristic for voltage-gated K + channels, like the selectivity filter motif GYGD, positive charges in S4 and the S4-S5 linker as well as characteristic residues in the Cterminal end of S6.
The multiple sequence alignment (MSA) of K + channel sequences was performed using ProbCons [23].Subsequent phylogenetic analysis was carried out via MEGA6 by using the Maximum Likelihood method based on the amino acid replacement model developed by Whelan And Goldman (WAG).The WAG model combines attributes of the amino acid replacement models of Dayhoff and JTT and has been shown to provide a better overall fit to the evolutionary process [24,25].Initial trees for the heuristic search were obtained by applying the Neighbour-Joining method to a matrix of pairwise distances estimated using a JTT model.A resampling analysis with 1000 bootstrap replications was applied and all positions with less than 15% site coverage were eliminated during tree construction.
Clustering analyses were performed by calculating the average distance between sequences using the substitution matrix BLOSUM62.Analyses were executed through the function 'Average Distance Using BLOSUM62' that is implemented in JalView on the basis of the previous obtained MSA [26].Resulting trees were loaded in MEGA6 to prepare final tree representation.For average distance analyses on sequence fragments the full-length sequences were separated according to the proposed position of transmembrane domains provided by the multiple sequence alignment.New MSA were executed for fragmented sequences and average distance trees calculated as described for full-length sequences.
To calculate the conservation of sequence parts the conservation values given by JalView were used.For each column in the multiple sequence alignment a value between 0 (not conserved) and 11 (100% conserved) is assigned.The average of these values was calculated for the columns representing the N-and C-termini and the core channel structure, respectively.Average conservation of the cytosolic termini was 0.05, while the average conservation of the core channel part was 1.82.Thus, the core region shows a more than 30 times higher conservation than the termini.

Structural data analyses
The protein crystal structure of Kv1.2 in the open state (2A79) was downloaded from www. rcsb.org.Kv1.2 in the closed state [27] as well as KAT1 and SKOR in both conformational states [28] are only available as modelled protein structures.For the KAT1 and SKOR models 6ns of molecular dynamics simulation are available.
Analyses of protein structures were performed in VMD by executing tcl-scripts to calculate the data described in the following.Potential interactions between voltage sensor and pore on the extra-and intracellular side were identified by searching for residues that are within 3Å of residues in the C-terminal end of S1 or the S4-S5 linker, respectively.The diameter of the ring that is formed by the four S4-S5 linkers was calculated by measuring the distance of two opposite residues from the middle of the linker using the function 'bond' from VMD.With four linkers this calculation resulted in two distances, whose mean was calculated to determine the final S4-S5 linker ring diameter.The 'bond' function was also used to calculate the distance of neighbouring S6 ends.The mean was calculated of four distances to determine the average distance of S6 ends.To calculate the bend of S6 three points were defined within S6: One point in the bend around the PxP motive (animal Kv) or the conserved glycine (HCN, plant Kv-like), and two more at the beginning and the end of S6.To place the points in the centre of the αhelix three amino acids were selected and the centre between them measured.Subsequently, the angle at the point in the bend was measured by calculating the angle between the two lines that arise when the middle point is connected with the points at both ends of S6.

Results and Discussion
For our comparison of voltage-gated potassium (K + ) channels from plants and animals we used as prototypes the nine Kv-like channels from the model plant Arabidopsis thaliana, 12 hyperpolarization-activated cyclic nucleotide-gated cation channels (HCNs) and 10 Kv channels from Rattus norvegicus, Mus musculus and Homo sapiens, as well as the HCN from the sea urchin, Strongylocentrotus purpuratus, and the Shaker channel from Drosophila melanogaster.
According to their voltage-dependence, these channels were categorized into hyperpolarization-activated, inward-rectifying (K in ) and depolarization-activated, outward-rectifying (K out ) channels [12].Based on their voltage-dependence and on their highest permeability for K + , HCN channels are considered as animal K in channels although they are no strict K + channels but also permeable to other cations, like e.g.Na + [15].In total 33 protein sequences (Table 1) were investigated in phylogenetic as well as average distance analyses.Phylogenetic analysis allowed the classification of channels according to their evolutionary relationship, whereas the average distance analysis was performed on full-length sequences and sequence fragments to classify voltage-gated K + channels with respect to certain criteria, like their rectification, for instance.The latter analyses allowed the exclusion of structural channel parts that are not responsible for the variability in channel gating.
Furthermore, the inclusion of available crystal structures and homology protein models in the comparison enabled the investigation of channel attributes that are not accessible by the primary protein sequence, like distances between structural parts or intermolecular interactions.

K + channels of different function and origin vary most in the loop between S5 and the pore domain
The phylogenetic analysis grouped K + channels into three major clades that coincide with the separation according to their origin and function (Fig 2).Plant voltage-gated K + (plant Kvlike) channels form one clade, while the other two comprise animal K in (HCN) and animal K out (Kv) channels, respectively.Interestingly, plant Kv-like and animal K in channels are sister clades indicating a close evolutionary relationship of both channel groups.
The multiple sequence alignment (MSA) shows similar sequence features that support the close relationship of HCN and plant Kv-like channels.The accuracy of the MSA was tested by annotating the position of transmembrane domains (TMD).The positions of TMD S1 to S6 identified according to their presence in available crystal structures and protein models were well aligned in the MSA (S1 Fig).
Structural features of defined K + channel groups are illustrated in Fig 3 .In general, the different channel classes are very distinct in their cytosolic parts, while the core transmembrane structures are rather similar.In the MSA the average conservation of positions within the protein core structures is more than 30 times higher than the average conservation of positions in the cytosolic channel parts (for details see Materials and Methods).Plant Kv-like and animal K in channels had the longest C-termini with lengths of 295 to 680 amino acids, while among the considered animal K out channels, only Kv2.1 contains a long C-terminus with comparable length (440 amino acids).All other animal K out channels exhibit shorter C-termini ranging from less than 100 to 230 amino acids.In contrast, the N-termini are little conserved and variable even within the K + channel groups.
In the more similar transmembrane core region, plant K in and K out as well as animal K in channels have an extended loop of equal length between TMD S2 and S3 when compared to the other channels.The most striking difference, however, was the presence and length of a loop between TMD S5 and the pore domain.The S5-P loop is the only conserved structural difference that clearly distinguishes between animal K in and K out channels as well as between plant K in and K out channels, albeit not in a uniform manner.An extended S5-P loop is only present in animal K in and plant Kv-like channels, but not in animal K out channels.With about 16 amino acids, the S5-P loop has the same length in animal and plant K in channels.In plant K out channels, however, the loop is even longer with about 24 amino acids.
Further differences between main K + channel groups were identified within the distribution of charged residues in TMD S4 and the S4-S5 linker.Plant Kv-like channels show the lowest number of charged residues in S4 (5-6), but the highest number in the S4-S5 linker (5)(6).In contrast, animal K out channels exhibit a higher number of charged residues in S4 (7-8), but hardly contain charged residues in the S4-S5 linker (1-2).When both structural parts are considered as a unit, animal K in channels show the highest number of charged residues (~13).The variable number and distribution of charged residues in S4 and the linker may account for differences in the sensitivity towards the membrane voltage, but cannot explain the fundamental differences in coupling voltage-sensing to pore-opening; i.e. the opposite rectification behaviour.Ath-KAT2 AT4G18290 Ath-AKT1 AT2G26650 Ath-AKT5 AT4G32500 Ath-AKT6 AT2G25600 Ath-GORK plant K out AT5G37500

Rno-HCN3 Q9JKA8
Rno-HCN4 Q9JKA7 To better illustrate the classification, only the organization in K + channel groups is shown instead of displaying every K + channel separately.In the following step we calculated average distances to allow channel clustering; here, instead of using full-length sequences, only channel fragments were compared: (i) only the N-termini, (ii) the TMDs S1 to S4, (iii) the isolated S4 segment, (iv) the S4-S5 linker, (v) the pore module comprising the TMDs S5 to S6, (vi) the On the basis of the regions S1-S4, S5-S6, S6 and the C-termini multiple sequence alignments followed by average distance analyses clustered K + channels in a similar manner as observed for full-length sequences (Fig 4D , 4G, 4H and 4I).In all these cases, plant K in and K out channels were associated with animal K in channels in the same clade.A slightly different picture arose when considering the isolated S4 segment, the S4-S5 linker or the N-terminus (Fig 4C , 4E and 4F).On the basis of S4 sequences, animal K in and K out channels clustered in a clade separated from plant K + channels (Fig 4E The cluster analyses of sequence fragments of K + channels allowed three main conclusions: (i) Plant K in and K out channels clustered always together, while animal K in and K out channels showed independent classification patterns.(ii) The charged TMD S4 alone did not contain sequence information that discriminated between different voltage-dependencies of channels, whereas the whole voltage sensor, S1-S4, did allow to distinguish between animal K out channels on one side and animal K in and plant Kv-like channels on the other.This observation confirmed experimental data proving that the voltage-dependent S4 movement in plant K in , animal The N-termini of these K + channels are sparsely conserved.Even within channel groups there is a high variability in length (represented by the white to black gradient).Some animal K out channels show an extended sequence between the first two TMDs.In contrast, all plant Kv-like and animal K in channels have an extended loop between TMD S2 and S3.Animal K in channels have the highest number of charged residues in S4, whereas the lowest number is present in plant Kv-like channels.In contrast, plant Kv-like channels feature the highest number of charged residues in the S4-S5 linker, while animal K out channels hardly possess charged residues in this linker.The most striking structural difference that correlates with channel function is a loop between TMD S5 and the pore region.This loop is only present in plant Kv-like and animal K in channels.In addition, a small stretch in the loop discriminates between plant K in and K out channels.doi:10.1371/journal.pone.0137600.g003K in and animal K out channels is similar and does not explain the opposite rectification properties (Figs 4 and 5 and Table 2, [19][20][21]29].(iii) The S4-S5 linker did not comprise sequence information that discriminates between channel functions.The linker is a short stretch of a few amino acids that is tightly connected to S4.Therefore, the linker has to follow the universal S4 movement, and as a stand-alone structural element it could not explain the gating variability.Nevertheless, the S4-S5 linker is a structure with fundamental function in the gating process.Its peculiar importance derives from its interaction with the C-terminal part of S6 [8,10,30,31].

Spu
By all means, neither S4 nor the S4-S5 linker could explain the rectification behaviour of voltage-gated K + channels.Thus, rectification properties seem to derive from structures surrounding S4 and the S4-S5 linker, rather than from these structures itself.Previous experimental work supports that information on the rectification behaviour is likely spread over The simplified tree that represents the average distance analysis of full-length K + channel sequences corresponds to the tree in Fig 2 .(B) Arrows in the schematic channel representation mark separation points used for sequence fragmentation.Plant Kv-like and animal K in channels cluster in one clade if only those channel structures are considered that code for S1 to S4 (D), S5 to S6 (G), only S6 (H), and the C-terminus (I).If only the N-termini are considered plant Kv-like channels cluster together with only some of the animal K in channels (C).In contrast to the complete voltage sensor unit S1 to S4 (D), S4 alone does not differentiate between animal K in and K out channels (E).Here, animal K in and K out channels cluster together and plant Kv-like channels are present in a separated clade.The S4-S5 linker sequences are the only structural part that clusters plant Kv-like and animal K out channels in one clade, which is separated from animal K in channels (F).doi:10.1371/journal.pone.0137600.g004these structures.There, we obtained experimental results supporting the interplay of various channel parts in the rectification process of plant Kv-like channels [28,32,33].Domain swapping between plant K in and K out channels results in destruction of channel rectification or functionality, or does not affect the rectification at all.However, not a single K in -K out chimera with inverted rectification properties has been identified so far.This led to the conclusion that the interplay of various channel parts is essential to obtain rectification.

The shape of the sixth TMD correlates with conformational channel states
Next, we considered crystal structures of voltage-gated K + channels and available homology models built on the basis of these structures [27,28,34].In all these representations, the TMD S6 shows two structural characteristics that change during transition from the open to the close conformation and vice versa.The C-terminal ends of S6 bend around a glycine residue in plant Kvlike and animal K in channels or a PxP motif in animal K out channels [35][36][37][38].Due to this bend, the ends of S6 can change their position and approach or depart from each other during channel state transition.In all K + channel groups, the C-terminal S6 ends behaved in a conserved manner, they were closer in the closed conformation and more distant in the open conformation (Table 2 and Fig 5).The angle of the S6 bend also correlated with channel states.It was larger in the open conformation, which indicates that S6 is straightened during channel opening.

Two ubiquitous interfaces connect voltage sensor and pore but exhibit channel-class specific differences
In general, the voltage-sensing module and the pore module are structurally connected via a short amino acid sequence at the intracellular side of the channel-the S4-S5 linker.A second  2. Channel characteristics are illustrated for animal and plant Kv-like channels in the open and in the closed state.The representation of animal K in channels is grey shaded because no protein structures are available so far.In general, the voltage sensor is displayed in blue and the pore in orange.The pore opening at the intracellular side is according to the channel states wide or narrow.The width of this pore part correlates with the distance of the C-terminal S6 ends.S4 is marked within the voltage sensor in its up or down position, respectively.The S4-S5 linkers of all four subunits surround the inner part of the channel pore and form a kind of ring around it.Its diameter correlates with the S4 position and is larger when S4 is in the up position.A black dot at the extracellular side of Kv channels represents the occurrence of the connection of voltage sensor and pore.The loop between the C-terminal end of S5 and the pore TMD is displayed at the extracellular side of the pore in plant Kv-like and animal K in channels.The bending of S6 is illustrated within the pore.It is more bended in the closed states and straighter in the open states.Red stars in plant Kv-like and animal K out channels represent the interaction of the S4-S5 linker and S6.An interaction of these both structures has been described for animal K in channels as well.But it is not clear in which state it appears.Therefore, a yellow star is used to illustrate the interaction.doi:10.1371/journal.pone.0137600.g005interface has been described at the extracellular side between the C-terminal end of S1 and the junction between TMD S5 and the pore domain [11].This connection is formed by interacting amino acid side chains.The examination of structures of channel proteins in the open and closed conformation showed that in all K + channel groups certain amino acids in the C-terminal region of S1 were in close proximity (<3Å) to the S5-P region (Table 2 and S3 Fig).
Although these two links appeared to be ubiquitous features of voltage-gated K + channels, there were channel class-specific differences: (i) The extracellular S1-S5/P connection involved regions that were shown to discriminate between K + channel groups (Figs 3 and 5).It is thus likely that the functional features of this junction can be distinct in the different groups.However, with the available data it is not clear whether the S5-P loop is directly involved in this interaction or whether this interaction restricts rather to close-by sequence regions.For animal K out channels this connection has been identified by analysing Kv channel sequences and its importance has been experimentally confirmed by the examination of Shaker channel mutants [11].The plant Kv-like channel models do not contain the S5-P loop, because they have been modelled on the bases of the animal K out channel that lacks this loop.However, domain swaps of the loop between the plant Kv-like channels KAT1 and SKOR showed its critical role for channel functionality [32,33,39].(ii) At the intracellular side, the S4-S5 linker interacts with the C-terminal end of S6 as shown experimentally and on the basis of structural data (Table 2) [5,8,9,33,40].Molecular modelling data suggested that this interaction was not the same in all voltage-gated K + channels.Protein models of plant Kv-like channels proposed a permanent connection between S4/S5 and S6 in the case of the K out channel SKOR.In contrast, in the case of the K in channel KAT1 the connection was found in the open state only.Interestingly, also the crystal structure of the animal K out channel chimera Kv1.2/2.1 shows the S4/S5-S6 Data derived from: + protein models * crystal structures; E experiments; or S sequence analysis.'n.a.' not available data due to lacking structure.E' an interaction between the S4-S5 linker and the S6 end has been observed experimentally, but it is not clear in which state the interaction occurs.connection only in the open but not in the closed channel state [40].In context animal K out channels appear to be like plant K in channels.
For animal K in channels on the other hand it has been shown that the so called C-linker right after S6 is interacting with residues from the S4-S5 linker and has significant impact on the correct coupling between voltage-sensor and pore module [41,42].The C-linker of animal K in channels consists of six α-helices that connect S6 with the C-terminus.The helix closest to S6, the A' helix, is located underneath the channel and is able to interact with the S4-S5 linker.In animal K out channels the C-termini are rather short.Instead, the N-termini are pronounced and interact underneath the channel [34].In the case of plant Kv-like channels C-terminal structural data is missing.However, experimental data provide room for speculation that plant K in and K weak channels might share the same voltage-sensor and pore coupling like animal K in channels (reviewed in [43]).
Taken together, the behaviour of the voltage-sensing domain S4 and the pore-lining domain S6 were conserved among inward-and outward-rectifying K + channels in animals and plants.The two connections of voltage sensor and pore were conserved as well, but bear channel-class specific differences.These differences might provide variability in the flexibility of the overall K + channel structure, which in turn results in an unequal forwarding of the universal voltage sensor movement to the pore.
Plant K in and K out channels are both more closely related to animal K in channels than to animal K out channels Structural and cluster analyses clearly pointed to differences that allow discriminating between animal K in and K out channels.Such segregation into the two functional types was not possible for plant Kv-like channels.Phylogenetic and average distance analyses always grouped all plant Kv-like channels within the same clade, despite their functional differences (Figs 2 and 4).For the plant K weak channel AKT2 this observation is not surprising.As it has been described as specialised K in channel that switches between K in and K weak channel properties dependent on its phosphorylation status [44][45][46].Plant K out channels on the other hand exhibit a contrary physiological function.Surprisingly, they were included in the K in channel clade rather than in the K out channel clade indicating that plant K out channels are structurally closer to K in channels.This implies that plant and animal K out channels bear structural differences that do not affect the physiological function of K + efflux.To solve this paradoxical evolutionary and functionally classification we propose that animal K out and plant K out channels occurred independently of each other in life history.In other words this would mean that the feature of outward-rectification of K + channels has been invented twice in evolution.
To challenge this hypothesis in an extended phylogenetic analysis, we included sequences of voltage-gated K + channels from (i) the moss Physcomitrella patens, (ii) the ancient vascular plant species Selaginella moellendorffii, (iii) the freshwater alga Klebsormidium flaccidum, (iv) bacteria and (v) Caenorhabditis elegans (Fig 6).The phylogenetic analysis corroborated the close relationship of plant Kv-like channels to animal K in channels, while bacterial and animal K out channels form a distant group.Interestingly, the eleven K + channels found in the oldest of the considered plants, K. flaccidum, are not clustered in one clade.Instead they spread over the entire tree.Two of them, Kfl-21-370 and Kfl-119-90, are closely related to plant Kv-like channels sharing ~40% identity in the channel core region.Indeed, Kfl-21-370 and Kfl-119-90 show several features that are exclusive characteristics of plant Kv-like channels, like a comparable number of charged residues in S4 and the S4-S5 linker, additional amino acids between S2 and S3 as well as the longer loop between S5 and the pore domain ( Fig 7).Interestingly, from a structural point of view, the Klebsormidium K + channels are intermediates between plant K in Fig 6 .Two Klebsormidium voltage-gated K + channels are related to Kv-like channels of higher land plants.A maximum likelihood tree was calculated using the WAG model and a bootstrap value of 1000 with K + channel sequences from (i) bacteria (Aeropyrum pernix: Ape-KvAP (APE_0955); Listeria monocytogenes: Lmo-Kv (lmo2059); Methanocaldococcus jannaschii: Mja-Kv (MJ0139)), (ii) animals (see Table 1; Caenorhabditis elegans: Cel-Kv (CELEg_ZK1321.2)),(iii) early land plants (Selaginella moellendorffii: Smo-Kout (scaffold_121_2_ORF2559)); Physcomitrella patens: Ppa-Kin1 (Pp1s3_156V6.1),Ppa-Kin2 (Pp1s2_170V6.1),Ppa-Kin3 (Pp1s22_165V6.1),Ppa-Kin4 (Pp1s283_74V6.1),(iv) eudicots (see Table 1), and (vi) the early terrestrial alga Klebsormidium (Kfl-21-370 (kfl00021_0370), Kfl-21-430 (kfl00021_0430), Kfl-35-100 (kfl00035_0100), Kfl-43-330 (kfl00043_0330), Kfl-108-190 (kfl00108_0190), Kfl-119-90 (kfl00119_0090), Kfl-175-250 (kfl00175_0250), Kfl-373-10 (kfl00373_0010), Kfl-409-60 (kfl00409_0060), Kfl-437-30 (kfl00437_0030), Kfl-600-80 (kfl00600_0080)).The phylogenetic tree is mid-point rooted and displays the relationship of K + channels.The main K + channel families are grouped in separated clades.Plant K in and K out channels as well as animal K in and K out channels form individual clades.Thereby, bacterial and animal K out channels are separated from K + channels from green plants as well as animal K in channels.The K + channels from Klebsormidium are and K out channels: (i) In its core region Kfl-21-370 is 41% to the plant K in prototype KAT1 and 35% to the K out channel SKOR.Likewise, Kfl-119-90 is 43% identical to KAT1 and 42% to SKOR.(ii) The length of the S5-P loop is two amino acids longer than in plant K in channels and six amino acids shorter than in plant K out channels.(iii) A characteristic motif in TMD S6 (DMI in K out channels; NLG in K in channels), which is important for the K + sensing mechanism of plant K out channels [39], is not yet developed in Klebsormidium K + channels.(iv) In TMD S1 Kfl-21-370 and Kfl-119-90 share a motif with plant K in channels.(v) The expanded K + -selectivity filter (TxxTxGYGD) is that of plant K in channels and at the x-positions clearly different from the highly conserved motif of plant K out channels.Taken together, the two K + channels from Klebsormidium exhibit more plant specific structural channel features than Kv-like channels from bacteria or animals, and they are slightly closer to plant K in channels than to plant K out channels.Thus, it is likely that plant K in and K out channels evolved in the plant lineage from a common ancestor.This assumption is further encouraged by the studies performed by Gomez-Porras and colleagues [47].They demonstrated in comprehensive phylogenetic analyses that the transition of green plants from an aqueous to a terrestrial environment was attended by a dramatic reduction of the K + channel diversity.The fact that structural K in features can be found in K + channels from mosses and algae but typical structural K out features can first be found in K + channels from vascular plants suggests that plant K out channels evolved from K in channels.

Conclusion
The comparison of voltage-gated K + channels of different kingdoms revealed new insights into the structural and evolutionary background of K + channel rectification.The entire plant Kvlike channel family likely emerged from a common ancestor and is closer related to animal K in than to animal K out channels.Thus, plants probably developed outward-rectification of Kv-like channels anew and independently of the outward-rectification that is found in animals.Several structural differences and similarities supporting this finding have been identified throughout not grouped into one clade.Instead they spread over the entire tree indicating a broad diversity of different K + channels.Especially two, out of eleven, K + channels from Klebsormidium (marked by a square) are closely related to plant K in and plant K out channels.doi:10.1371/journal.pone.0137600.g006Fig 7 .Simplified representation of structural features of Kv-like channels from higher land plants and Klebsormidium.The core channel parts S1 to S6 are displayed including features that are characteristic for these K + channel groups, like the extended S2-S3 and S5-P loops, conserved positive charges in S4 and the S4-S5 linker, and the selectivity filter motif TxxTxGYGD.Furthermore, two motifs in the transmembrane domains S1 and S6 are displayed.Outward-rectifying K + channels contain an SFFT motif in S1, where inward-rectifying channels contain an AWxx motif.In S6 the characteristic DMI (K out channels) and NLG (K in channels) motifs are illustrated.In Klebsormidium the S6 motif is not yet specified, while the motif in S1 is that of K in channels.The extended loop between S5 and P is two amino acids longer than in plant K in channels and six amino acids shorter than the corresponding loop in plant K out channels.doi:10.1371/journal.pone.0137600.g007voltage-gated K + channel groups.Two known voltage sensor and interfaces have been observed as ubiquitous traits in voltage-gated K + channels, the S4-S5 linker on one hand and the connection of S1 with the S5-P region on the other hand; even though, both interfaces bear channel class-specific differences in their specifications.S1 forms constant interactions with or close to the variable S5-P loop, while the S4-S5 linker establishes variable interactions with the C-terminal end of S6.All in all, several smaller structural differences appear to play together to provide variable flexibilities to voltage-gated K + channel structures, which may result in variable transitions of the voltage sensor movement to the pore and therefore in an inverse coupling of voltage sensor and pore opening in K in and K out channels.

Fig 1 .
Fig 1. Structure and function of voltage-gated K + channels (A) Functional voltage-gated K + channels form tetramers (coloured representation).Each subunit (grey representation) comprises six transmembrane domains (S1-S6) and a selectivity filter between S5 and S6.The first four TMDs form the voltage sensor that is located in the periphery of the functional channel.The last two TMDs contribute to the channel pore in the centre of the channel.(B) TMD S4 of the voltage sensor translocates in response to the transmembrane voltage within the membrane and affects the channel conformation.It moves towards the intracellular side of the membrane upon hyperpolarizing membrane voltages and towards the extracellular side upon membrane depolarization.The "down" movement of S4 causes K in channels to open, while the same S4 movement results in closure of K out channels.The opening of K out and closure of K in channels proceeds in the inverse manner in response to S4 "up" movements.
.1371/journal.pone.0137600.t001K + channel rectification properties cannot be pinpointed to certain isolated channel parts K + channel sequence clustering by phylogenetic analyses reflected well a separation according to their origin and their function (Fig 2).Fig 4A illustrates the result of Fig 2 as a simplified tree.

Fig 2 .
Fig 2. Plant Kv-like channels are closer related to animal K in than to animal K out channels.Phylogenetic analysis classifies voltage-gated K + channels according to their function and appearance.The maximum likelihood tree was calculated using the WAG model and a bootstrap value of 1000.Plant Kv-like channels, animal K in and animal K out channels are each classified into one clade.Thereby, plant Kv-like and animal K in channels belong to the same superior clade, while animal K out channels are organized into a second superior clade.Dashed lines indicate the two mayor clades.Please note: The mid-point rooted tree allows conclusions about channel relationships, but not about assumptions of the ancestors.doi:10.1371/journal.pone.0137600.g002 ).When considering only the S4-S5 linker, plant Kv-like and animal K out channels were grouped together (Fig 4F).The classification based on N-terminal sequences was less strict, what is expected when taking into account the high variability in N-termini length (Fig 4C).

Fig 3 .
Fig 3. Simplified representation of structural characteristics of voltage-gated K + channel groups.Voltage-gated K + channel groups are defined according to the clades that arose from the phylogenetic analysis presented in Fig 2.The N-termini of these K + channels are sparsely conserved.Even within channel groups there is a high variability in length (represented by the white to black gradient).Some animal K out channels show an extended sequence between the first two TMDs.In contrast, all plant Kv-like and animal K in channels have an extended loop between TMD S2 and S3.Animal K in channels have the highest number of charged residues in S4, whereas the lowest number is present in plant Kv-like channels.In contrast, plant Kv-like channels feature the highest number of charged residues in the S4-S5 linker, while animal K out channels hardly possess charged residues in this linker.The most striking structural difference that correlates with channel function is a loop between TMD S5 and the pore region.This loop is only present in plant Kv-like and animal K in channels.In addition, a small stretch in the loop discriminates between plant K in and K out channels.

Fig 4 .
Fig 4. Simplified mid-point rooted trees based on average distance analyses of fragmented voltage-gated K + channel sequences.(A)The simplified tree that represents the average distance analysis of full-length K + channel sequences corresponds to the tree in Fig 2.(B) Arrows in the schematic channel representation mark separation points used for sequence fragmentation.Plant Kv-like and animal K in channels cluster in one clade if only those channel structures are considered that code for S1 to S4 (D), S5 to S6 (G), only S6 (H), and the C-terminus (I).If only the N-termini are considered plant Kv-like channels cluster together with only some of the animal K in channels (C).In contrast to the complete voltage sensor unit S1 to S4 (D), S4 alone does not differentiate between animal K in and K out channels (E).Here, animal K in and K out channels cluster together and plant Kv-like channels are present in a separated clade.The S4-S5 linker sequences are the only structural part that clusters plant Kv-like and animal K out channels in one clade, which is separated from animal K in channels (F).

Fig 5 .
Fig 5. Illustration of voltage-gated K + channel features described in Table2.Channel characteristics are illustrated for animal and plant Kv-like channels in the open and in the closed state.The representation of animal K in channels is grey shaded because no protein structures are available so far.In general, the voltage sensor is displayed in blue and the pore in orange.The pore opening at the intracellular side is according to the channel states wide or narrow.The width of this pore part correlates with the distance of the C-terminal S6 ends.S4 is marked within the voltage sensor in its up or down position, respectively.The S4-S5 linkers of all four subunits surround the inner part of the channel pore and form a kind of ring around it.Its diameter correlates with the S4 position and is larger when S4 is in the up position.A black dot at the extracellular side of Kv channels represents the occurrence of the connection of voltage sensor and pore.The loop between the C-terminal end of S5 and the pore TMD is displayed at the extracellular side of the pore in plant Kv-like and animal K in channels.The bending of S6 is illustrated within the pore.It is more bended in the closed states and straighter in the open states.Red stars in plant Kv-like and animal K out channels represent the interaction of the S4-S5 linker and S6.An interaction of these both structures has been described for animal K in channels as well.But it is not clear in which state it appears.Therefore, a yellow star is used to illustrate the interaction.
M and C in the column 'structure' indicate availability of protein models (M) or crystal structures (C).Number given in brackets following the 'C' indicates PDB ID of the respective structure.doi:10.1371/journal.pone.0137600.t002

Table 1 .
Overview of voltage-gated K + channel sequences used in this study.

Table 2 .
Overview of characteristic features found in voltage-gated K + channels of different taxa.