Skip to main content

An extracellular cation coordination site influences ion conduction of OsHKT2;2



HKT channels mediate sodium uniport or sodium and potassium symport in plants. Monocotyledons express a higher number of HKT proteins than dicotyledons, and it is only within this clade of HKT channels that cation symport mechanisms are found. The prevailing ion composition in the extracellular medium affects the transport abilities of various HKT channels by changing their selectivity or ion transport rates. How this mutual effect is achieved at the molecular level is still unknown. Here, we built a homology model of the monocotyledonous OsHKT2;2, which shows sodium and potassium symport activity. We performed molecular dynamics simulations in the presence of sodium and potassium ions to investigate the mutual effect of cation species.


By analyzing ion-protein interactions, we identified a cation coordination site on the extracellular protein surface, which is formed by residues P71, D75, D501 and K504. Proline and the two aspartate residues coordinate cations, while K504 forms salt bridges with D75 and D501 and may be involved in the forwarding of cations towards the pore entrance. Functional validation via electrophysiological experiments confirmed the biological relevance of the predicted ion coordination site and identified K504 as a central key residue. Mutation of the cation coordinating residues affected the functionality of HKT only slightly. Additional in silico mutants and simulations of K504 supported experimental results.


We identified an extracellular cation coordination site, which is involved in ion coordination and influences the conduction of OsHKT2;2. This finding proposes a new viewpoint in the discussion of how the mutual effect of variable ion species may be achieved in HKT channels.


Proper sodium (Na+)/potassium (K+) homeostasis is a crucial requirement for high yielding plant growth. Although Na+ ions can promote plant growth at low concentrations, they turn hazardous at increasing levels [1,2,3]. Due to similar physicochemical properties, Na+ ions can mimic functions of K+ ions, bridging periods of K+ shortage [4, 5]. However, at high concentrations, Na+ ions compete with K+ and cause K+ deficiency symptoms in plants, since Na+ mimics the K+ function only incompletely [6, 7].

Plants have developed sophisticated mechanisms to cope with salt stress. These systems aim to avoid high cytosolic Na+ concentrations in cells of plant shoots by compartmentalization and retrieval of Na+ from the xylem sap [8, 9]. HKT channels form one comprehensive family that is involved in Na+ usage and detoxification [10,11,12,13]. In monocots, like rice, two types of HKT channels exist: (i) class I-type channels, which are mainly Na+ selective, and (ii) channels of class II acting as Na+/K+ co-transporters [14, 15]. The transmembrane structure of HKT channels is composed of one subunit containing four transmembrane-pore-transmembrane (MPM) units, which group together and form the conduction pathway in the center of the protein. Among HKTs, there is a wide variety of transport kinetics, cation selectivity and rectification properties even within the two classes [16,17,18,19,20,21].

Several studies have demonstrated the interactive effect of different ion species for various HKT channels. In particular, the mutual effects of Na+ and K+ have been documented [18, 21, 22]. However, the mechanism by which ion species affect each other and modulate the ion transport process is so far not well understood. Initially, it was discussed that HKT channels form multi-ion pores [17, 21]. In this scenario, multiple ion binding sites exist in the pore allowing a coupled ion transport without additional conformational changes [23]. However, recently Böhm and colleagues experimentally demonstrated that in the Venus flytrap HKT1 a maximum of one ion at a time occupies the selectivity filter [24]. Therefore, the mechanism behind mutual ion species effects remains controversial.

In this study, we provide new functional-structural perspectives for the understanding of ion species-specific effects on ion conduction in HKT channels. We built a homology model of the class II-type Oryza sativa HKT channel OsHKT2;2 embedded in a lipid membrane, performed molecular dynamics simulations in the presence of Na+ and K+ ions and evaluated protein-ion interactions. Through systematic analysis, we identified a potential cation coordination site at the extracellular surface of the protein – a region that ions have to pass before entering the pore. The cation coordination site was characterized by computational studies and its functional importance for OsHKT2;2 was experimentally validated using electrophysiological methods. Additionally, analyses of in silico mutants further underpinned the experimental insights.


Electrophysiological characterization of OsHKT2;2

OsHKT2;2 is one out of eight to nine HKT channels (depending on the cultivar) expressed in rice plants [14]. Members of the HKT family show remarkable functional diversity regarding ion selectivity, rectification properties and effect of external cation compositions [20]. To provide a solid basis for this study, we first characterized the functional properties of OsHKT2;2 expressed in Xenopus leavis oocytes using the Two-Electrode Voltage-Clamp (TEVC) technique.

Under voltage clamp conditions, cells expressing OsHKT2;2 conducted mainly Na+ inward currents in the absence of K+. Thereby, ion transport rates raised with increasing Na+ concentrations (Fig. 1b, d, i). At low extracellular Na+ concentrations (0.3 mM NaCl), ion currents of several hundred nanoamperes were detected, while at high Na+ concentrations (30 mM NaCl) currents increased to several microamperes consistent with Na+ being the ion carrying the current. A 100-fold increment of Na+ concentration (from 0.3 mM to 30 mM) resulted in an almost 5-fold increase of ion conduction.

Fig. 1
figure 1

Electrophysiological characteristics of OsHKT2;2 in Xenopus laevis oocytes. a Representative currents elicited in control cells (not injected with cRNA). Only currents for Na30K1 solution are shown since highest currents are expected under these conditions. Control currents for all conditions are shown in Additional file 2: Figure S2. b-h Representative OsHKT2;2 mediated whole cell currents recorded two days after cRNA injection. Currents were recorded at indicated Na+ and K+ concentrations: Na0K1 - no NaCl and 1 mM KCl, Na0K30 - no NaCl and 30 mM KCl, Na03K0 - 0.3 mM NaCl without KCl, Na03K1 - 0.3 mM NaCl and 1 mM KCl, Na30K0 - 30 mM NaCl without KCl, Na30K1 - 30 mM NaCl and 1 mM KCl, Na30K30 - 30 mM NaCl and 30 mM KCl. The holding potential was set to the zero current level (approximately -45 mV with solution Na30K1) and 1 s voltage pulses were stepped from +20 to -160 mV in -15 mV decrements. A 1.5 s resting interval was allowed between successive voltage steps. The current (μA) and time (ms) scale for all traces shown in A through H is shown at the bottom of panel A. i Mean current-voltage (IV) curves. Standard deviation calculated from data obtained from five to eight measurements obtained on five different days and from four different oocyte batches. Currents were normalized to currents measured at -145 mV in Na30K1 solution. Inset: Reversal potentials (Vrev) were determined based on the data plotted in I and plotted against the extracellular Na+ concentration. Symbols for Vrev and IV-curve data are identical and correspond to the same solutions in both plots. The external Na+ dependency of the Vrev values obtained with 1 mM KCl (filled triangle and diamond) or without KCl (empty triangle and circle) in extracellular solutions was fitted with a semi-logarithmic function (line K1 and K0). The dashed line represents the curve for an ideally Na+ permeable channel. j Currents elicited by a -145 mV voltage pulse, measured at the beginning (dark grey) or end (white) of the voltage pulse, and as average currents over the whole voltage pulse (light grey). Currents were normalized to currents measured in Na30K1 solution. Numbers in the plot indicate the degree of current reduction over time (from beginning to end of the voltage pulse). The magnitude of deactivation is stated in percentage. Means and standard deviation were obtained from data of five to eight measurements

The addition of low K+ concentrations (1 mM KCl) had a stimulating effect on the ion transport and led to a significant increase in ion conduction (Fig. 1d, e). One millimolar of KCl almost doubled the ion conduction at high and low Na+ concentrations. We are assuming a stimulating effect of K+ ions for two reasons: (1) In the absence of NaCl in the bath solution, no inward currents could be measured in the presence of 1 or 30 mM KCl (Fig. 1g, h). Currents at negative voltages were comparable to those recorded in control oocytes (compare Fig. 1g, h, i and Additional file 2: Figure S2). Therefore, potassium ions, at least in the absence of sodium ions, do not seem to pass the conduction pathway, or, they do pass with low velocity in a way the resulting currents cannot be clearly distinguished from background currents. (2) The stimulating effect was revoked under rising K+ concentration (30 mM NaCl/30 mM KCl; Fig. 1f, i, j). Interestingly, the ionic current measured with 30 mM NaCl/30 mM KCl was comparable to that observed without K+ ions in the extracellular solution indicating that K+ stimulates ion conduction only at low concentrations. Not only the increase of the Na+ concentration but also raising K+ concentrations provoked a shift of the reversal potential (Vrev) towards less negative membrane voltages (inset of Fig. 1i, for the K+-dependency compare Vrev values in 0 mM (K0), 1 mM (K1) and 30 mM KCl (K30) at constant Na+), indicating that -besides Na+- also K+ ions may pass OsHKT2;2. The observation that despite this shift of Vrev the current amplitude only rises at low KCl concentrations in the media but not at high concentrations might be explained by different ion permeation velocities depending on the ion composition. While the presence of a few K+ ions could break electrostatic interactions between Na+ and the channel and thus might increase the ion mobility, elevated K+ concentrations could induce new hindering interactions and reduce ion mobility again.

Additionally, under high Na+ concentrations, inward currents partially deactivated at negative voltages. This effect was most pronounced in the absence of K+ and reduced with increasing K+ concentrations (Fig. 1d-f, j). At 30 mM NaCl (no KCl), the current amplitude was reduced by 12.0% (±2.8%) over one-second voltage pulses. In the presence of 1 mM KCl, the current reduction was 10.7% (±4.6%), and at 30 mM KCl current deactivation was almost abolished (3.2% ± 2.2%).

Overall, we identified and characterized three functional properties for wild-type OsHKT2;2 which we used later in the study to compare functional changes in mutated OsHKT2;2 channels (Table 1). The functional properties of wild-type OsHKT2;2 include: (1) ionic currents enhance with increasing Na+ concentration, (2) low K+ concentrations (1 mM KCl) stimulate Na+ ion transport, and (3) high K+ concentrations (30 mM KCl) lack the stimulating effect on Na+ ion transport.

Table 1 Electrophysiological characterisation of OsHKT2;2 and its mutants. The presence of three current characteristics are specified according to their degree of occurrence: very strong (+++), medium (++) or weak (+)

Identification of an extracellular cation coordination site

In order to identify amino acids that are approached by ions before entering the pore, we built a homology model of the OsHKT2;2 wild-type channel on the basis of the structurally and functionally related bacterial channels KtrB (PDB ID 4J7C) and TrkH (PDB ID 3PJZ). KtrB conducts Na+-dependent K+ transport for which reason its selectivity filter provides a good template for OsHKT2;2 that conducts both cations as well [25, 26]. Furthermore, it has been shown that KtrB alone (in absence of the regulatory KtrA unit) conducts K+ and Na+ ions, which further justifies the use of KtrB as template for homology modelling [27]. For structural comparison of template and model structure see Additional file 1: Figure S1. Since the N-terminus of KtrB is not included in the crystal structure, N-terminal amino acids of OsHKT2;2 were modeled on the basis of the TrkH (see Methods for details). Subsequently, we performed molecular dynamics (MD) simulation in the presence of 10 mM NaCl and 10 mM KCl for 100 ns. This allowed evaluating the frequency at which ions approached residues over time. To demonstrate reproducibility, the MD simulation was performed threefold and a total of 300 ns simulation time was reached, which formed the basis of the following results.

Analyses of the contact frequency between ions and amino acids during MD simulations led to the identification of a potential extracellular cation coordination site (Table 2). This coordination site is located in the outer extracellular protein region approximately 20 Å away from the pore entrance (Fig. 2a). We identified three polar residues and one proline to be involved in the formation of a putative cation coordination site – P71, D75, D501, and K504. Within these, three are highly conserved among plant HKT channels (P71, D75, and K504), while at position D501 a negative charge is preserved in 75% of HKT channel sequences. Figure 2b illustrates the conservation of these four residues within the rice HKT family and additionally states their conservation within a multiple sequence alignment based on 20 experimentally characterized HKT channels.

Table 2 Contact frequencies of cations and selected amino acids. Absolute and relative frequencies of cations approaching (being within 4 Å) the four amino acids forming the potential extracellular cation coordination site (bold entries) and other amino acids in their proximity. Frequencies are based on 300 ns MD simulations and are displayed for sodium and potassium ions separately, respectively, as well as accumulated for both cations
Fig. 2
figure 2

Localization and conservation of extracellular cation coordination site. a OsHKT2;2 is displayed as transparent surface representation with the selectivity filter and surrounding residues marked in the center of the channel as red surface representation. The cation coordination site is marked on the outer extracellular edge of the protein also as red surface representation, and the position of the four contributing residues is shown. The approximate distance between coordination site and pore is indicated in the top view representation. Detailed position of coordination site forming residues is illustrated on the right in the cartoon representation of OsHKT2;2. The four MPM units with their respective transmembrane (M) and pore domains (P) as well as P71, D75, D501 and K504 are marked. b An extract of the multiple sequence alignment of all eight HKT members present in Oryza sativa. Conservation of the residues P71, D75, D501 and K504 forming the coordination site is indicated. ConSurf grades were calculated on the basis of a multiple sequence alignment of 20 experimentally characterized HKT sequences, and are shown below the sequence alignment. The grades range from 1 (variable) to 9 (conserved) [29]. D75 and K504 are 100% conserved (ConSurf grade 9). P71 is in 90% of HKT sequences conserved (ConSurf grade 8), and contains in 5% of the sequences a lysine or serine. At position 501 a negative charge is somewhat conserved as aspartic in 60% of the sequences (ConSurf grade 4) or glutamic acid (in 15% of the sequences). In the remaining 25% of the sequences, this position holds an asparagine (10%), proline (5%), lysine (5%), or tryptophan (5%). Amino acid sequences used for this analysis are the following: AtHKT1 (Q84TI7), EcHKT1;1 (Q9LLM4), EcHKT1;2 (Q9FVK0), HvHKT2;1 (Q4VWZ3), McHKT1;1 (Q84UZ6), OsHKT1;1 (Q7XPF8), OsHKT1;3 (Q6H501), OsHKT1;4 (Q7XPF7), OsHKT1;5 (Q0JNB6), OsHKT2;1 (Q0D9S3), OsHKT2;2 (Q93XI5), OsHKT2;3 (Q8L481), OsHKT2;4 (Q8L4K5), PhaHKT2;1 (AB234305), PutHKT2;1 (C7B119), SlHKT1;1 (CCJ09641), SlHKT1;2 (CCJ09642), TaHKT1;4 (A0N0D7), TaHKT1;5 (A6XCE1), TaHKT2;1 (Q41515)

Structural attributes of the extracellular cation coordination site

Both sodium and potassium ions approached a defined area on the protein surface during MD simulations more often than residues in close proximity. K+ ions were found in this region between 0.5 to 1% of the simulation time, while Na+ stayed there about three times longer (1.2 to 5.2% of the time; Fig. 3a, Table 2). At this point, it is important to clarify that we were not searching for a region in the sense of an ion binding site that binds ions for extended time periods with potential allosteric impacts on protein function. Instead, we were searching for areas, which ions approach and pass on their way into the pore – in other words, coordination sites that attract cations and guide them to the pore entrance. This implicates that the total time of residence of an ion within the identified region is comparatively short as observed in our MD simulations. Nevertheless, cations approach the site frequently and are coordinated by the identified residues, in contrast to randomly approached residues where no ion coordination occurs. Throughout the study, we will refer to this region as extracellular cation coordination site.

Fig. 3
figure 3

Oxygen atoms of P71, D75 and D501 coordinate cations. a Total time that K+ and Na+ stay in the coordination site during each of the three simulations expressed in percentage. b Coordination of cations happens via the backbone oxygen atom of P71, one side chain oxygen atom of D75 and both side chain oxygen atoms of D501 (red spheres). Average ion-atom distances for sodium and potassium are indicated in the figure. Sodium: 1) 2.38 ± 0.14 Å, 2) 2.29 ± 0.09 Å, 3) 4.02 ± 0.21 Å, 4) 2.90 ± 0.65 Å, 5) 2.58 ± 0.31 Å and 6) 3.74 ± 0.25 Å. Potassium: 1) 2.74 ± 0.19 Å, 2) 2.63 ± 0.30 Å, 3) 4.15 ± 0.40 Å, 4) 2.92 ± 0.47 Å, 5) 2.90 ± 0.38 Å and 6) 3.96 ± 0.28 Å. c Exemplary plot of ion-atom distances for sodium. Distances of one sodium ion and all six charged atoms facing the potential coordination site are represented for 5 ns of MD system 1 while Na+ is inside the pocket. The close proximity (< 2.5 Å) of Na+ to the backbone oxygen of P71, one side chain oxygen of D75 and both side chain oxygen atoms of D501 is illustrated, indicating their ability to coordinate sodium

Despite their short residence in the coordination site, the cations are coordinated by surrounding amino acids. A closer examination revealed four oxygen atoms coordinating the cations in this region: the backbone oxygen of P71 (Fig. 3b, c.1), one side chain oxygen atom of D75 (Fig. 3b, c.2) and the two side chain oxygen atoms of D501 (Fig. 3b, c.4 and c.5). The average distances between the mentioned oxygen atoms and Na+ were 2.38 Å, 2.29 Å, 2.90 Å and 2.58 Å, respectively, and, 2.74 Å, 2.63 Å, 2.92 Å and 2.90 Å for K+ (Fig. 3b). The coordination of one Na+ ion over time is exemplarily shown in Fig. 3c. Besides the proximity (< 2.5 Å) of the Na+ ion to the four mentioned oxygen atoms (plots c.1, c.2, c.4, c.5), the proximity to other side chain atoms facing towards the predicted coordination site is also shown. Namely, the proximity to the second side chain oxygen atom of D75 (plot c.3) and the positively charged nitrogen atom of K504 (plot c.6). The second side chain oxygen atom of D75 is oriented away from the coordination site most of the time and cannot contribute to the coordination of cations. As expected, the nitrogen of lysine 504 is also not involved in the coordination, since the positive charges of cations and the nitrogen atom repel each other.

Nevertheless, lysine 504 is located between the two aspartic acids D75 and D501, which suggests an interaction between these residues. Salt bridges may be formed between oppositely charged residues when the charged atoms are closer than 4 Å [28]. Average distances between the positively charged nitrogen atom of lysine 504 and the negatively charged oxygen atoms of both aspartates fulfill this requirement (Fig. 4a, b). Indeed, both oxygen atoms of D75 are in constant proximity to the nitrogen of lysine 504 (Fig. 4b, upper two traces) with an average distance of 2.83 Å and 2.82 Å, respectively (Fig. 4c). The two oxygen atoms of D501 alternate in their proximity to this nitrogen (Fig. 4b, lower two traces). If in proximity, the distance is 2.68 Å on average (Fig. 4c). Thus, lysine 504 may interact with both aspartates and acts as a structural component to obtain the geometry of the coordination site by keeping the aspartates in place rather than coordinating ions.

Fig. 4
figure 4

K504 can form salt bridges with D75 and D501. a Representation of potential salt bridges between K504 and D75 and D501. The median oxygen-nitrogen distance is stated. b Exemplary plot of oxygen-nitrogen distances throughout one complete simulation. The formation of salt bridges is considered possible when oxygen-nitrogen distances are below the threshold of 3.2 Å. The red line in each plot indicates this threshold. c Distance distribution of indicated oxygen-nitrogen distances over all three simulations. The median is stated in each box plot. Outliers are not drawn. For each of the two D501 oxygen atoms only those time spans where the respective oxygen atom was orientated towards K504’s nitrogen atom were used

Overall, our computational studies suggest that the extracellular cation coordination site is constituted by a positively charged lysine (K504) that holds two negatively charged aspartates (D75 and D501) in place, which form a negative environment that attracts cations. The side chain oxygen atoms of both aspartates, in conjunction with the backbone oxygen of P71 coordinate Na+ and K+ ions.

Experimental validation of the predicted cation coordination site

To study the relevance of the cation coordination site for the ion conduction of OsHKT2;2, we mutated P71, D75, D501, and K504 and experimentally characterized the changes in functionality in these mutants. Alanine mutations of all four residues were performed to assess the effect of side chain removal, and asparagine mutants of D75 and D501 allowed the evaluation of negative charge neutralization. Besides, K504 was mutated to arginine to examine the relevance of the side chain size, to glutamine to analyze the effect of charge neutralization, and to glutamic acid to evaluate the inversion from a positive to a negative charge. Overall, mutation of residue K504 had the most significant impact on OsHKT2;2’s functionality. Modifications of D75 and D501 affected the kinetic features slightly, whereas mutation of P71 did not cause changes in the transport mechanism as compared to the wild-type channel (Table 1).

Mutations of P71, D75, and D501 show no or minor functional alterations

To evaluate the effect of mutations of the cation coordinating residues we modified P71, D75 and D501 by the removal of side chains (alanine mutants) and charges (asparagine mutants). The mutant P71A behaved similarly to the wild-type channel regarding the three functional properties that have been established before on the basis of wild-type OsHKT2;2: Na+ conductivity (1), the mutual ion-species effect at low (2) and high (3) K+ concentrations (Table 1). P71A transported Na+ in the absence of K+, and the ion transport increased at high Na+ concentration. As observed in the wild-type channel, ion currents increased in the presence of 1 mM K+ but not 30 mM K+, when compared to currents recorded in the absence of K+ (Fig. 5).

Fig. 5
figure 5

Mutation of P71 does not affect channel functionality. a Mean current-voltage (IV) curves and standard deviation calculated from data obtained from at least three different oocytes. b Representative currents measured at the indicated NaCl and KCl concentrations in Xenopus laevis oocytes two days after cRNA injection. A pulse at holding potential (zero current level) was followed by 1 s voltage pulses from +20 to -160 mV in -15 mV decrements and continued with a final pulse at holding potential for 1.5 s

Mutations of residues D75 and D501 affected the ion conduction slightly. As observed in the wild-type channel, D75A/N and D501A/N transported Na+ in the absence of K+ in a concentration-dependent manner (0.3 mM NaCl and 30 mM NaCl). However, the increase of the current amplitude was slightly lower for D75A and D501A (4-fold), and only 3-fold in the D501N mutant, compared to a 5-fold increase for the wild-type channel (Fig. 6). The strong stimulating effect of low K+ concentrations on the ion conduction was observed in all mutants (30 mM NaCl/1 mM KCl). However, in contrast to wild-type OsHKT2;2, low magnitudes of current activation in the presence of high K+ concentrations (30 mM NaCl/30 mM KCl) were still observed in D75N and D501N mutants (Fig. 6c, f, g and Additional file 3: Figure S3) indicating that the mutual effect of K+ on Na+ currents had slightly shifted.

Fig. 6
figure 6

Mutations of D75 and D501 maintain channel functionality with only small changes in the kinetics. a, c, d, f Mean current-voltage (IV) curves and standard deviation calculated from data obtained from at least three different oocytes. b, e, g Currents measured at -145 mV at the beginning (dark grey) or end (white) of voltage pulses and of average currents over the whole voltage pulse (light grey) of D75A and D75N (b), D501A (e) and D501N (g). Currents were normalized to currents measured in Na30K1 solution. Means and standard deviations were obtained from data of three to five different oocytes

In summary, residues D75 and D501 were weakly sensitive towards modifications, while P71 was not being affected by the mutation. The latter was expected since proline coordinates cations via its backbone oxygen atom, which remains unaltered by mutations.

K504 mutants showed substantial differences in their functionality compared to wild-type OsHKT2;2

Lysine 504 was identified as a structural component of the cation coordination site that seems to take a crucial position in the overall geometry. Indeed, K504 was very sensitive towards residue substitutions. Removal of the side chain (mutant K504A) and inversion of the charge (mutant K504E) rendered OsHKT2;2 in a manner that no ion conduction could be detected with the established measuring protocol. K504A and K504E injected oocytes showed no currents two days after cRNA injection (Fig. 7e, f). Each, fifteen injected oocytes have been analyzed on three different measuring days. However, oocytes injected with K504A cRNA sporadically conducted low currents in the nanoampere range, which differed significantly from those recorded in control oocytes. These observations were occasionally made three days after cRNA injection. We speculated that the K504A mutant may be functional but conducts ions with extremely low efficiency. In this case, the increased protein expression over time may lead to increased whole-cell currents through OsHKT2;2-K504A that are only detectable after long incubation times.

Fig. 7
figure 7

Mutations of K504 affect conduction behaviour of Na+ substantially. a-c, e-i Representative currents recorded in Xenopus laevis oocytes two days after cRNA injection at indicated Na+ and K+ concentrations: Na03K0 - 0.3 mM NaCl without KCl, Na30K0 - 30 mM NaCl without KCl, Na30K1 - 30 mM NaCl and 1 mM KCl. A pulse at holding potential (zero current level) was followed by 1 s voltage pulses from +20 to -160 mV in -15 mV decrements and continued with a final pulse at holding potential for 1.5 s. d, j Mean current-voltage (IV) curves and standard deviation calculated from data obtained from three to six measurements

In contrast, lysine mutants K504R and K504Q were conducting currents although with substantial differences in their transport characteristics as compared to wild-type OsHKT2;2 and all other mutants. Na+ currents in the absence of K+ were low even at high Na+ concentrations (30 mM NaCl in comparison to 0.3 mM NaCl). In fact, the current amplitude rose not more than 2-fold with increasing Na+ concentration (compared to the 5-fold increase in wild-type channels) and remained in the nanoampere range. Therefore, the increase in the external Na+ concentration was not able to enhance the ion conduction to the same extent as in wild-type OsHKT2;2 (Fig. 7b, d, h, j). Regardless, low K+ concentrations did not only maintain its strong stimulating effect but also was this effect even more pronounced in the lysine mutants (30 mM NaCl and 30 mM NaCl/1 mM KCl, Fig. 7c, i). In the background of 30 mM NaCl, the addition of 1 mM KCl resulted in a more than 4- (K504Q) to 6-fold (K504R) increase of the current amplitude, while the wild-type was stimulated just about 2-fold. A further increase in the external K+ concentration was less inhibitory in the K504Q and K504R mutants than in the wild-type. While the wild-type showed almost identical current amplitudes at -160 mV measured in 30 mM Na+ and at 30 mM Na+ /30 mM K+, the amplitudes of the two mutants were about twice as large in 30 mM Na+ /30 mM K+ compared to those in 30 mM Na+ only (Fig. 7d, j and Additional file 4: Figure S4). This observation indicates that the Na+ transportability of OsHKT2;2 is affected in the K504R and K504Q mutants, while K+ ions are still able to promote and, in these mutants, also to rescue proper ion conduction.

In silico mutations support experimental results

To gain further understanding of the nature of the experimental results, the four lysine mutants were studied in silico by performing 100 ns MD simulations. In each simulation, we examined the time that Na+ and K+ ions approached one of the four coordination site residues (< 4 Å) and the total number of approaches throughout the MD simulation (Table 3). The combined assessment of the duration of stay and number of approaches allowed the evaluation of the average time an ion spent in the cation coordination site. It is striking that Na+ ions approached the coordination site of the slow- or non-conducting lysine mutants K504A and K504E more frequently than in the conducting mutants K504Q and K504R (35 and 85 times versus 15 and 16 times; Table 3). The approaches observed in the mutants K504Q/R were comparable to the number of Na+ ion approaches in wild-type OsHKT2;2.

Table 3 Ion approaches of external coordination site in K504 in silico mutants compared to OsHKT2;2-wt. The proportion of time which cations stay within 4 Å of coordination site forming residues. Time is expressed as a percentage during 100 ns MD simulations. In addition, the number of approaches per ion species is indicated

In OsHKT2;2-K504R, where the positive charge was kept but the side chain slightly increased, Na+, as well as K+, approached the coordination site only shortly (see the short duration of stay and the low number of approaches in Table 3). Both cations were most of the time close to residue D501 (Na+: 2.5% and K+: 3.2% of simulation time), while in the wild-type the duration of stay was more balanced between the four coordination site residues. It is plausible that a more voluminous arginine residue at position 504 is preventing cations from entering into the coordination site due to the repulsion of positive charges. Although, the surface area around the coordination site is still negatively charged and able to attract cations, the positive charge at position 504 is bigger in the case of the mutant and further exposed to the extracellular space than in the wild-type channel, which may cause repulsion of cations (Fig. 8e, j compared to a, f).

Fig. 8
figure 8

Coordination site constitution of in silico mutants K504A, K504E, K504Q and K504R. a-e Electrostatic potential mapped to the protein surface. The cation coordination site is highlighted and the position of P71, D75, D501 and K504 are indicated. f-j Representative illustration of the coordination site constitution of all four K504 mutants and the wild-type with ions approaching from the extracellular site or being inside the coordination site

In mutant K504Q, where the charge of the amino acid was neutralized, cations stayed at least five times longer in the coordination site (12.9 and 26.8% for Na+ and K+, respectively) than in mutant K504R, although both mutants showed comparable ion transport abilities in electrophysiological experiments. In contrast to mutant K504R, cations entered into the coordination site in mutant K504Q, which may account for the more extended stay in the pocket (compare Fig. 8i, j). Additionally, K+ ions approached the coordination site six times more frequently than in K504R, which may account at least in part for the longer duration of stay of K+ in the coordination site. Compared to the wild-type, cations stayed significantly longer close to the coordination site. The absence of the positive charge may allow cations to stay longer in the negatively charged coordination site before continuing to the pore entrance of the channel.

In the slow- or non-conducting mutants K504A and K504E, Na+ approached the coordination site with high frequency (35 and 85 times) and stayed up to 9.6% (K504A) and 27.4% (K594E) of the simulation time. K+ ions stayed shorter time than Na+ ions in the coordination site in mutant K504A (up to 7.3%) and longer in mutant K504E (up to 38.7%). The insertion of a negatively charged glutamate in mutant K504E generated a large negative area on the protein surface compared to OsHKT2;2-wt (Fig. 7c compared to a). A negative electrostatic potential highly attracts cations, and more than one cation approached the ion coordination site at a time, which may trap cations in the coordination site (Fig. 7g).

In summary, in silico data provide insights into potential reasons for functional alterations and the non-conduction of K504 mutants, respectively, and support experimental data.


HKT channels have been described as crucial players in salt tolerance. They transport Na+ as uniporter (class I) or in symport with K+ ions (class II). Class II-type HKT channels function as Na+ uniporter in the absence of K+. However, they become Na+/K+ symporter in the presence of K+. Analyses of the reversal potential suggest that K+ may be transported as well in the absence of Na+. However, it may be transported at such low velocity that resulting conduction cannot be distinguished from background currents in electrophysiological measurements. However, it has been shown in yeast complementation experiments that a K+ uptake deficient yeast strain could be complemented by expressing OsHKT2;2, and that K+ was depleted from the medium cultivating OsHKT2;2-expressing yeast [22]. Besides, K+ uptake by OsHKT2;2 was assessed in a plant expression system by analyzing rubidium uptake in tobacco protoplasts. Also here, a very low rubidium uptake could be detected in the absence of Na+ [21]. The mutual effect of ion-species on the ion conduction has been observed for several HKT channels. Dependent on the channel member the influence on ion transport varies [20]. Our electrophysiological studies confirmed the stimulating effect that K+ ions exert on Na+ transport in OsHKT2;2. Interestingly, this effect was only seen at low K+ concentrations and decreased as the extracellular K+ concentration increased. This phenomenon has been previously reported for OsHKT2;2 [21, 22]. From the physiological point of view, this is a consistent behavior. During K+ starvation (low K+ concentrations), ion uptake increases and with it the uptake of Na+ ions. Physico-chemically similar Na+ ions may complement functions of K+ up to a certain degree and may enable the survival of the plant while K+ availability is low. With increasing concentration of K+, the ion conduction is reduced and with it the uptake of Na+ ions. In the presence of sufficient K+ ions, high Na+ influx is not desired since Na+ ions compete with K+ and may cause symptoms of K+ starvation. However, while the physiological response of HKT is well investigated, the ion-specific effect at the molecular level is still unknown.

To get deeper insights into the molecular basis behind the mutual ion-species effect, we modeled the structure of OsHKT2;2 and analyzed cation-protein contacts. Our computational studies show that four residues P71, D75, D501 and K504 form an extracellular cation coordination site. Thereby, P71, D75 and D501 coordinate cations and the positively charged lysine (K504) holds the two negatively charged aspartates, D75 and D501, in place. Thus, a negative environment that attracts cations is formed.

We confirmed the relevance of the coordination site for the ion conduction of OsHKT2;2 in functional electrophysiological studies. Out of the three residues involved in ion coordination, residue D501 is slightly more sensitive towards modifications (see Table 1 for summary). D75 and P71 are located nearby on one extracellular end of a transmembrane helix and contribute to the cation coordination with each one oxygen atom. D501 is placed on the opposite side on the extracellular end of another helix and coordinates with two oxygen atoms (Fig. 2a, right). Therefore, mutating D501 removes all coordinating atoms from this side of the coordination site, while mutating D75 is only eliminating one of two coordinating atoms from the opposite side, which may account for the slightly stronger effect seen in D501 mutants. In either case, these mutations maintain the overall functionality of OsHKT2;2 and affect its kinetics only slightly. Considering this coordination site as one of many regions that ions pass before entering the pore, it is plausible that mutations affect the channel kinetics instead of rendering the channel non-functional. The identified coordination site would be involved in the attraction of cations and their forwarding towards the pore entrance, rather than in the tight and constant binding of ions. Therefore, modifying the coordinating oxygen atoms would first affect transport dynamics.

Lysine 504, on the other hand, is a crucial residue in the coordination site. Inverting the charge from positive to negative by a glutamate substitution results in a non-conductive mutant. With a third negative charge in the coordination site, the two aspartates repulse and drift apart, which leads to an expansion of the coordination site. Additionally, the three negative charges form a large negative region on the protein surface that attracts cations and traps them in the coordination site. If the identified coordination site is understood as one stopover for ions on their way into the pore, the trapped ions may explain the inhibited ion conduction, since they might block the path for ions leading into the pore.

Overall, the positively charged lysine 504 may be the necessary impulse that cations need to be forwarded towards the pore. The positive charge and moreover its correct size and position are crucial for channel functionality. The sensitivity towards changes of one of the parameters size or charge indicates how fine-tuned the constitution of the coordination site is. Removing the side chain by mutations to alanine renders OsHKT2;2 seemingly slow-conducting. Neutralizing the positive charge (K504Q) or increasing the size of the amino acid side chain with simultaneous maintenance of its charge (K504R), led to substantial changes in the conductance of OsHKT2;2. In the latter, the Na+ uniport is strongly reduced at high Na+ concentrations. A recent study of Xu and co-workers [30] showed that mutating G490 to arginine in HKT1;5 from Triticum monoccocum abolished its Na+ transport ability completely. G490 from TmHKT1;5 is a highly conserved glycine among class I and class II HKT channels which corresponds to G503 in OsHKT2;2 and is the neighboring residue of K504 (analyzed in this study). The introduction of a voluminous positive charge in TmHKT1;5-G490R is abolishing the Na+ transport, which supports our finding of a delicate constitution of the coordination site. Despite the reduced Na+ transport ability in K504 mutants, at low K+ concentrations, these mutants still mediate large ion currents indicating that either the coordination site still works correctly for K+ ions or that the positive effect of K+ ions has another origin. One possibility may be, e.g. the existence of another ion coordination site that is more relevant for the K+ forwarding than the one described here. In fact, during simulations, we could identify additional regions that are approached by ions on their way into the pore. These regions are currently under investigation.

The importance of K504 has been shown earlier for the wheat TaHKT1 channel [31]. In this study, positively charged residues of the second transmembrane domain of the last MPM unit were mutated to investigate their contribution to ion conduction. Mutant K508Q, which corresponds to K504Q in our present study, showed reduced functionality in comparison to the wild-type channel in electrophysiological experiments and yeast complementation studies. Compared to wild-type TaHKT1, TaHKT1-K508Q complemented a K+ uptake deficient yeast strain although less effective than the wild-type. In electrophysiology experiments, the Na+ transport ability of TaHKT1-K508Q was affected, a result that we also observed in OsHKT2;2-K504Q. The authors also proposed a potential salt bridge formation between K508 and D78 based on a homology model of TaHKT1. TaHKT1-D78 corresponds to D75 in OsHKT2;2, which we found to form salt bridges with K504, which is in agreement with the prediction made by Kato and colleagues. Those results and the ones presented in our study point to the functional relevance of K504 in combination with D75. Our study additionally provides a biological and structural context by identifying these residues as part of a cation coordination site.


Overall, using computational methods, we identified an extracellular cation coordination site and validated its functional relevance experimentally. Although being about 20 Å far from the pore entrance, the identified coordination site is essential for proper ion conduction. Crucial ion coordination sites on the extracellular protein surface provide an example of how different ion species may affect the channel conduction behavior even before entering the pore. We are still far from understanding how K+ ions influence the ion transport of OsHKT2;2 on the molecular level. This study, however, provides new evidence, which opens the discussion for new explanations for the mutual effect of ion species.


Homology modelling and molecular dynamics simulation

The structural model of OsHKT2;2 (accession number: Q93XI5) was generated using the I-Tasser server [32]. Full-length structural models were built from multiple threading alignments. For threading, the bacterial channels KtrB (PDB ID 4J7C) and TrkH (PDB ID 3PJZ) were used. The mayor part of the channel was modeled on the basis of KtrB (residues 39–530), while the N-terminus of OsHKT2;2 (residues 1–38) was modeled on the basis of TrkH, since the crystal structure of KtrB does not contain the N-terminal region. I-TASSER calculates two values, which are recommended for evaluating model quality and correct folding. Model quality is estimated by the C-score, where a value higher than −1.5 indicates a generally correctly folded protein. The second value, the TM-value, gives a measure of structural similarity compared to the native structure. A TM-value higher than 0.5 indicates proteins of similar folds. C-score and TM-value of the best model generated by I-TASSER were − 0.4 and 0.66, respectively, which is above the threshold, and indicates that the generated model represents a correctly folded protein.

OsHKT2;2 contains a 29 amino acid long extracellular loop that is not present in bacterial channels. Due to a missing template, the loop spanning residues 473–501 were modeled ab initio by I-TASSER due to the missing template. To improve the model quality, we modeled that part with the FALC-Loop server, a program that is specialized in loop modeling, and which combines statistical and knowledge-based methods [33]. The improved full-length model was refined using i3Drefine [34] and minimized using CHARMM via its web portal CHARMMing [35, 36]. Throughout the whole process, model quality was monitored independently by quality assessment tools like Procheck [37], QMEAN [38] and ModFold [39, 40]. The refined model was embedded in a pre-equilibrated POPC bilayer solvated with pre-equilibrated TIP3P water molecules [41] in a periodic boundary condition box using VMD [42]. The system was neutralized by the addition of sixteen chloride ions. Each, six Na+ and K+ ions as well as twelve further chloride ions were added to mimic the concentration of 10 mM NaCl and 10 mM KCl. The entire system underwent several minimization steps and was equilibrated with successive restraint reduction until restraints were removed entirely to ensure step-by-step equilibration and avoid degeneration of the protein. Each reduction of restrictions was monitored by checking system temperature, pressure and energy, as well as, the RMSD value of the protein. The consequent reduction of restriction was not initiated unless the RMSD value was not stable. This procedure was continued until all restrictions were removed entirely and the RMSD converged to a stable value. Not until then, were the production runs of the MD simulations initiated. Three replicates of 100 ns molecular dynamics (MD) simulations were generated using NAMD [43, 44].

Coordination site residues were identified by counting the frames in which an ion was within 4 Å of a given residue during MD simulations. Ion coordination was examined by measuring distances between charged atoms of P71, D75, D501, K504 and the ion present in the coordination site. Salt bridges were evaluated using VMD’s Salt Bridge plugin. Electrostatic potentials were calculated using the APBS web server with the ionic strength set to simulate 10 mM NaCl and 10 mM KCl [45]. APBS input files were generated using the PDB2PQR Server version 2.1.1 [46, 47].

In silico mutations of OsHKT2;2 were generated via VMD’s Mutate plugin. Mutants were minimized and equilibrated followed by 15 ns of MD simulation in NAMD.


Wild-type and mutant OsHKT2;2 constructs were cloned into the pNB1u vector for expression in Xenopus laevis oocytes using the USER-cloning technique [48]. Utilized primers were: GGCTTAAUatgacgagcatttaccaagaa (forward) and GGTTTAAUctaccatagcctccaatatt (reverse). Vector-specific and OsHKT2;2-specific sequences are given in uppercase and lowercase format, respectively. The position of uracil is underlined. Standard fusion PCR technique was used for site-directed mutagenesis. After DNA linearization with NotI, cRNA was synthesized using the mMessage mMachine in vitro transcription kit following the manufacturer’s guidelines.

For expression in oocytes, stage V and VI oocytes were harvested from Xenopus laevis and kept in ND96 solution (96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 1 mM MgCl2, 2.5 mM Na Pyruvate, 5 mM Hepes – pH 7.5, 50 mg mL− 1 gentamycin and 0.4 g L− 1 BSA). Oocytes were defolliculated by collagenase treatment in ND96 without CaCl2, gentamycin and BSA and kept overnight at 18 °C in complete ND96 solution as described initially. All animal procedures, including husbandry, oocyte harvesting and post-treatment, were performed in accordance with Cornell University IACUC Protocol number 2017–0139. Animals were purchased from Xenopus Express Inc. (Brooksville, FL, USA). 50 nl of cRNA (500 ng/μl) were microinjected into oocytes using an oil-driven injection system (Nanoject II Auto-Nanoliter Injector, Drummond Scientific Company, US). Cells were then incubated for two days at 18 °C in complete ND96 solution. Whole-cell currents were recorded using conventional Two-Electrode Voltage-Clamp technique (GeneClamp 500 amplifier and Digidata 1320A-PClamp 10 data acquisition system, Axon Instruments). Recordings were carried out under constant perfusion of bath solutions containing 2 mM MgCl2, 1.8 mM CaCl2, 10 mM MES, pH 5.5 with Tris-Base with the addition of: 1 mM KCl and 165 mM Sorbitol (Na0K1); 30 mM KCl and 130 mM Sorbitol (Na0K30); 0.3 mM NaCl and 165 mM Sorbitol (Na03K0); 0.3 mM NaCl, 1 mM KCl and 165 mM Sorbitol (Na03K1); 30 mM NaCl and 130 mM Sorbitol (Na30K0); 30 mM NaCl, 1 mM KCl, and 130 mM Sorbitol (Na30K1); 30 mM NaCl, 30 mM KCl and 100 mM Sorbitol (Na30K30).

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.


  1. Leigh RA, Wyn Jones RG. A hypothesis relating critical potassium concentrations for growth to the distribution and functions of this ion in the plant cell. New Phytol. 1984;97(1):1–13.

    Article  CAS  Google Scholar 

  2. Maathuis FJM. Physiological functions of mineral macronutrients. Curr Opin Plant Biol. 2009;12(3):250–8.

    Article  CAS  PubMed  Google Scholar 

  3. Marschner P. Marschner's mineral nutrition of higher plants. London: Academic; 2012.

    Google Scholar 

  4. Benito B, Haro R, Amtmann A, Cuin TA, Dreyer I. The twins K+ and Na+ in plants. J Plant Physiol. 2014;171(9):723–31.

    Article  CAS  PubMed  Google Scholar 

  5. Horie T, Costa A, Kim TH, Han MJ, Horie R, Leung H-Y, et al. Rice OsHKT2;1 transporter mediates large Na+ influx component into K+-starved roots for growth. EMBO J. 2007;26(12):3003–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Hasegawa PM, Bressan RA, Zhu J-K, Bohnert HJ. Plant cellular and molecular responses to high salinity. Annu Rev Plant Physiol Plant Mol Biol. 2000;51:463–99.

    Article  CAS  PubMed  Google Scholar 

  7. Pardo JM, Quintero FJ. Plants and sodium ions: keeping company with the enemy. Genome Biol. 2002;3(6):reviews1017.1–4.

    Article  Google Scholar 

  8. Flowers T, Läuchli A. Sodium versus potassium: substitution and compartmentation. In: Inorganic plant nutrition, vol. 15B; 1983. p. 651–81.

    Google Scholar 

  9. Adams E, Shin R. Transport, signaling, and homeostasis of potassium and sodium in plants. J Integr Plant Biol. 2014 Mar;56(3):231–49.

    Article  CAS  PubMed  Google Scholar 

  10. Maser P, Eckelman B, Vaidyanathan R, Horie T. Altered shoot/root Na+ distribution and bifurcating salt sensitivity in Arabidopsis by genetic disruption of the Na+ transporter AtHKT1. FEBS Lett. 2002;531(2002):157–61.

    Article  CAS  PubMed  Google Scholar 

  11. Horie T, Hauser F, Schroeder JI. HKT transporter-mediated salinity resistance mechanisms in Arabidopsis and monocot crop plants. Trends Plant Sci. 2009;14(12):660–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Møller IS, Gilliham M, Jha D, Mayo GM, Roy SJ, Coates JC, et al. Shoot Na+ exclusion and increased salinity tolerance engineered by cell type-specific alteration of Na+ transport in Arabidopsis. Plant Cell. 2009;21(7):2163–78.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Munns R, James RA, Xu B, Athman A, Conn SJ, Jordans C, et al. Wheat grain yield on saline soils is improved by an ancestral Na+ transporter gene. Nat Biotechnol. 2012;30(4):360–4.

    Article  CAS  PubMed  Google Scholar 

  14. Garciadeblás B, Senn ME, Bañuelos MA, Rodríguez-Navarro A. Sodium transport and HKT transporters: the rice model. Plant J. 2003;34(6):788–801.

    Article  PubMed  Google Scholar 

  15. Almeida P, Katschnig D, de Boer AH. HKT transporters--state of the art. Int J Mol Sci. 2013;14(10):20359–85.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Rubio F, Gassmann W, Schroeder JI. Sodium-driven potassium uptake by the plant potassium transporter HKT1 and mutations conferring salt tolerance. Science. 1995;270(5242):1660–3.

    Article  CAS  PubMed  Google Scholar 

  17. Gassmann W, Rubio F, Schroeder JI. Alkali cation selectivity of the wheat root high-affinity potassium transporter HKT1. Plant J. 1996;10(5):869–82.

    Article  CAS  Google Scholar 

  18. Horie T, Yoshida K, Nakayama H, Yamada K, Oiki S, Shinmyo A. Two types of HKT transporters with different properties of Na+ and K+ transport in Oryza sativa. Plant J. 2001;27(2):129–38.

    Article  CAS  PubMed  Google Scholar 

  19. Haro R, Bañuelos MA, Senn ME, Barrero-Gil J, Rodríguez-Navarro A. HKT1 mediates sodium uniport in roots. Pitfalls in the expression of HKT1 in yeast. Plant Physiol. 2005;139(3):1495–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Jabnoune M, Espeout S, Mieulet D, Fizames C, Verdeil J-L, Conéjéro G, et al. Diversity in expression patterns and functional properties in the rice HKT transporter family. Plant Physiol. 2009;150(4):1955–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Yao X, Horie T, Xue S, Leung H-Y, Katsuhara M, Brodsky DE, et al. Differential sodium and potassium transport selectivities of the rice OsHKT2;1 and OsHKT2;2 transporters in plant cells. Plant Physiol. 2010;152(1):341–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Oomen RJFJ, Benito B, Sentenac H, Rodríguez-Navarro A, Talón M, Véry A-A, et al. HKT2;2/1, a K+-permeable transporter identified in a salt-tolerant rice cultivar through surveys of natural genetic polymorphism. Plant J. 2012;71(5):750–62.

    Article  CAS  PubMed  Google Scholar 

  23. Su A, Mager S, Mayo SL, Lester HA. A multi-substrate single-file model for ion-coupled transporters. Biophys J. 1996;70(2):762–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Böhm J, Scherzer S, Shabala S, Krol E, Neher E, Mueller TD, et al. Venus flytrap HKT1-type channel provides for prey sodium uptake into carnivorous plant without conflicting with electrical excitability. Mol Plant. 2015;9(3):428–36.

    Article  PubMed  Google Scholar 

  25. Tholema N, Bakker EP, Suzuki A, Nakamura T. Change to alanine of one out of four selectivity filter glycines in KtrB causes a two orders of magnitude decrease in the affinities for both K+ and Na+ of the Na+ dependent K+ uptake system KtrAB from Vibrio alginolyticus. FEBS Lett. 1999;450(3):217–20.

    Article  CAS  PubMed  Google Scholar 

  26. Matsuda N, Kobayashi H, Katoh H, Ogawa T, Futatsugi L, Nakamura T, et al. Na+−dependent K+ uptake Ktr system from the cyanobacterium Synechocystis sp. PCC 6803 and its role in the early phases of cell adaptation to hyperosmotic shock. J Biol Chem. 2004;279(52):54952–62.

    Article  CAS  PubMed  Google Scholar 

  27. Tholema N, Vor der Brüggen M, Mäser P, Nakamura T, Schroeder JI, Kobayashi H, et al. All four putative selectivity filter glycine residues in KtrB are essential for high affinity and selective K+ uptake by the KtrAB system from Vibrio alginolyticus. J Biol Chem. 2005;280(50):41146–54.

    Article  CAS  PubMed  Google Scholar 

  28. Kumar S, Nussinov R. Close-range electrostatic interactions in proteins. Chembiochem. 2002;3(7):604–17.

    Article  CAS  PubMed  Google Scholar 

  29. Landau M, Mayrose I, Rosenberg Y, Glaser F, Martz E, Pupko T, et al. ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures. Nucleic Acids Res. 2005;33(Web Server):W299–302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Xu B, Waters S, Byrt CS, Plett D, Tyerman SD, Tester M, et al. Structural variations in wheat HKT1;5 underpin differences in Na+ transport capacity. Cell Mol Life Sci. 2018;75(6):1133–44.

    Article  CAS  PubMed  Google Scholar 

  31. Kato N, Akai M, Zulkifli L, Matsuda N, Kato Y, Goshima S, et al. Role of positively charged amino acids in the M2D transmembrane helix of Ktr/Trk/HKT type cation transporters. Channels (Austin). 2007;1(3):161–71.

    Article  Google Scholar 

  32. Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5(4):725–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Ko J, Lee D, Park H, Coutsias EA, Lee J, Seok C. The FALC-loop web server for protein loop modeling. Nucleic Acids Res. 2011;39(Web Server issue:W210–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Bhattacharya D, Cheng J. i3Drefine software for protein 3D structure refinement and its assessment in CASP10. PLoS One. 2013;8(7):e69648.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Miller BT, Singh RP, Klauda JB, Hodoscek M, Brooks BR, Woodcock HL. CHARMMing: a new, flexible web portal for CHARMM. J Chem Inf Model. 2008;48(9):1920–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B, et al. CHARMM: the biomolecular simulation program. J Comput Chem. 2009;30(10):1545–614.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr. 1993;26(2):283–91.

    Article  CAS  Google Scholar 

  38. Benkert P, Biasini M, Schwede T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011;27(3):343–50.

    Article  CAS  PubMed  Google Scholar 

  39. Maghrabi AHA, McGuffin LJ. ModFOLD6: an accurate web server for the global and local quality estimation of 3D protein models. Nucleic Acids Res. 2017;45(W1):W416–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. McGuffin LJ, Shuid AN, Kempster R, Maghrabi AHA, Nealon JO, Salehe BR, et al. Accurate template-based modeling in CASP12 using the IntFOLD4-TS, ModFOLD6, and ReFOLD methods. Proteins. 2018;86(Suppl 1(Suppl 10)):335–44.

    Article  CAS  PubMed  Google Scholar 

  41. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79(2):926–35.

    Article  CAS  Google Scholar 

  42. Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14(1):33–8 27–8.

    Article  CAS  PubMed  Google Scholar 

  43. Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, et al. Scalable molecular dynamics with NAMD. J Comput Chem. 2005;26(16):1781–802.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kalé L, Skeel R, Bhandarkar M, Brunner R, Gursoy A, Krawetz N, et al. NAMD2: greater scalability for parallel molecular dynamics. J Comput Phys. 1999;151(1):283–312.

    Article  Google Scholar 

  45. Unni S, Huang Y, Hanson RM, Tobias M, Krishnan S, Li WW, et al. Web servers and services for electrostatics calculations with APBS and PDB2PQR. J Comput Chem. 2011;32(7):1488–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004;32(Web Server issue):W665–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Dolinsky TJ, Czodrowski P, Li H, Nielsen JE, Jensen JH, Klebe G, et al. PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res. 2007;35(Web Server):W522–5.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Nour-Eldin HH, Hansen BG, Nørholm MHH, Jensen JK, Halkier BA. Advancing uracil-excision based cloning towards an ideal technique for cloning PCR fragments. Nucleic Acids Res. 2006;34(18):e122.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33(7):2302–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank Dr. Julian Schroeder (University of California, San Diego) for providing the clone of OsHKT2;2 wild-type. We also thank Alison Coluccio for technical assistance.


This work was supported by grants from the Chilean Fondo Nacional de Desarrollo Científico y Tecnológico ( to JR (No. 3150173), AV-J (No. 11170223), ID (No. 1150054) and WG (No. 1191133), and from Conicyt, programme PAI, Convocatoria nacional subvención a la instalación en la academia 2017 to JR (No. PAI77170035), and Conicyt, Programma de Cooperación Internacional, Concurso de Apoyo a la Formación de Redes Internacionales to JR, MP, ID and WG (REDES180117). Data collection was funded by Fondecyt No. 3150173 and Conicyt No. PAI77170035. Study design, analysis, data interpretation and writing was supported by all previously mentioned grants.

Author information

Authors and Affiliations



JR conceived and coordinated the study and wrote the paper. JR and AVJ designed, performed and analyzed the computational experiments. JR and MP designed, performed and analyzed the electrophysiological experiments. JR, ID and WG interpreted the data. All authors reviewed the results, revised the manuscript and approved its final version.

Corresponding author

Correspondence to Janin Riedelsberger.

Ethics declarations

Ethics approval and consent to participate

All animal procedures were performed in accordance with Cornell University IACUC Protocol number 2017–0139.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Structural comparison of KtrB template and OsHKT2;2 model. Template and model structure were aligned using the structural alignment algorithm TM-align [49]. A TM-score of 0.90536 was calculated indicating that both structures have the same fold. (A) Overlay of KtrB (PDB ID 4J7C) (purple) and modeled OsHKT2;2 (green) in side and top view. Alpha helices are displayed as tubes. Residues forming the selectivity filter and cation coordination site are represented in licorice. Selectivity filter forming residues in (1) KtrB: G62, G177, G280 and G392, (2) OsHKT2;2: G88, G244, G368 and G469. Cation coordination site forming residues in OsHKT2;2: P71, D75, D501 and K504. Corresponding residues in KtrB: L45, D49, T403 and K406. (B) Closer top view on cation coordination site (left) and selectivity filter (right). Positions of respective residues in the structure are illustrated by schematic representations to the left and right of the zoom. (TIF 7748 kb)

Additional file 2:

Figure S2. Currents elicited in control cells. (A) Representative currents recorded in Xenopus laevis oocytes at indicated Na+ and K+ concentrations: Na03K0 - 0.3 mM NaCl without KCl, Na03K1 - 0.3 mM NaCl and 1 mM KCl, Na30K0 - 30 mM NaCl without KCl, Na30K1 - 30 mM NaCl and 1 mM KCl, Na30K30 - 30 mM NaCl and 30 mM KCl. Control oocytes underwent the same handling and incubation procedure as injected oocytes and were measured on the same day as cRNA injected oocytes. A pulse at holding potential (zero current level) was followed by 1 s voltage pulses from +20 to -160 mV in -15 mV decrements and continued with a final pulse at holding potential for 1.5 s. (B) Representative current-voltage (IV) curves extracted form current traces shown in (A). For comparison the mean IV curve of OsHKT2;2-wt from Fig. 1g is presented (x). (TIF 574 kb)

Additional file 3:

Figure S3. D75A, D75N, D501A and D501N mutants behave comparable to OsHKT2;2-wt. Representative currents recorded in Xenopus laevis oocytes two days after cRNA injection at indicated Na+ and K+ concentrations: Na03K0 - 0.3 mM NaCl without KCl, Na03K1 - 0.3 mM NaCl and 1 mM KCl, Na30K0 - 30 mM NaCl without KCl, Na30K1 - 30 mM NaCl and 1 mM KCl, Na30K30 - 30 mM NaCl and 30 mM KCl. A pulse at holding potential (zero current level) was followed by 1 s voltage pulses from +20 to -160 mV in -15 mV decrements and continued with a final pulse at holding potential for 1.5 s. (TIF 1307 kb)

Additional file 4:

Figure S4. Mutants K504R and K504Q show altered kinetic characteristics in comparison to OsHKT2;2-wt. Representative currents recorded in Xenopus laevis oocytes two days after cRNA injection at indicated Na+ and K+ concentrations: Na03K0 - 0.3 mM NaCl without KCl, Na03K1 - 0.3 mM NaCl and 1 mM KCl, Na30K0 - 30 mM NaCl without KCl, Na30K1 - 30 mM NaCl and 1 mM KCl, Na30K30 - 30 mM NaCl and 30 mM KCl. A pulse at holding potential (zero current level) was followed by 1 s voltage pulses from +20 to -160 mV in -15 mV decrements and continued with a final pulse at holding potential for 1.5 s. (TIF 638 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Riedelsberger, J., Vergara-Jaque, A., Piñeros, M. et al. An extracellular cation coordination site influences ion conduction of OsHKT2;2. BMC Plant Biol 19, 316 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: