Tyrosine Residues from the S4-S5 Linker of Kv11.1 Channels Are Critical for Slow Deactivation*

Slow deactivation of Kv11.1 channels is critical for its function in the heart. The S4-S5 linker, which joins the voltage sensor and pore domains, plays a critical role in this slow deactivation gating. Here, we use NMR spectroscopy to identify the membrane-bound surface of the S4S5 linker, and we show that two highly conserved tyrosine residues within the KCNH subfamily of channels are membrane-associated. Site-directed mutagenesis and electrophysiological analysis indicates that Tyr-542 interacts with both the pore domain and voltage sensor residues to stabilize activated conformations of the channel, whereas Tyr-545 contributes to the slow kinetics of deactivation by primarily stabilizing the transition state between the activated and closed states. Thus, the two tyrosine residues in the Kv11.1 S4S5 linker play critical but distinct roles in the slow deactivation phenotype, which is a hallmark of Kv11.1 channels.

The rhythm of the heartbeat is primarily controlled by the coordinated activity of a range of voltage-gated ion channels (1,2). Kv11.1, which is encoded by the human ether-à-go-go-related gene (official name, KCNH2) (3), is responsible for the delayed rectifier current (I Kr ), which plays a critical role in cardiac repolarization (4,5). Kv11.1 channels have unusual kinetics of gating, namely slow activation and deactivation but fast inactivation and recovery from inactivation (6). Consequently, these channels pass relatively little current during the plateau phase of the action potential but deliver a large current during repolarization (7). The slow kinetics of deactivation during the early diastolic period also enable Kv11.1 channels to play a crit-ical role in terminating ectopic beats that may otherwise lead to cardiac arrhythmias and sudden cardiac arrest (6,8).
Kv11.1 channels have the typical topology of voltage-gated K ϩ channels, which includes cytoplasmic N-and C-terminal domains, a voltage sensor domain (VSD) 5 composed of the first four transmembrane segments (S1-S4), and a pore domain composed of the 5th and 6th transmembrane segments (S5 and S6). Pore domains from each of the four subunits co-assemble to form the ion conduction pathway, as well as the gates that control activation and inactivation. The slow activation and deactivation kinetics of Kv11.1 channels are due in part to the slow movement of the VSD (9,10). Motion of the VSD, in response to changes in the membrane electric field, facilitates opening and closing of the activation gate at the intracellular end of the pore domain by pulling on the S4S5 linker (residues 539 -550), which connects the VSD to the activation gate at the cytoplasmic end of the pore (11). Numerous studies have emphasized the importance of the S4S5 linker for Kv11.1 gating (11)(12)(13)(14)(15)(16)(17). In addition to providing a direct link between the voltage sensor and pore domains, the S4S5 linker interacts with the cytoplasmic domains and is thought to mediate the influence of cytoplasmic domains on the kinetics of deactivation gating (12,13,15). Recently, Lörinczi et al. (18) showed that when the S4S5 linker in Kv11.1 channels was disrupted (i.e. co-expressing a construct composed of residues 1-545 with another composed of residues 546 -1159), the channels can still assemble and retain their slow activation kinetics, as well as rapid inactivation and recovery from inactivation kinetics, but they have much faster deactivation kinetics compared with full-length wildtype (WT) channels (18). These data indicate that an intact S4S5 linker is not required for normal activation gating, but once the channels have opened, an intact S4S5 linker is required to stabilize the activated state and thereby contributes to slow deactivation. Another important implication of this study is that activation and deactivation gating are not simply reciprocal processes, but rather the transition from the closed state to the open state likely follows a different pathway to the transition from the open state back to the closed state.
In the crystal structure of the Kv1.2/2.1 chimeric channel, which is thought to have captured the channel in the activated state (19), the S4S5 linker forms an ␣-helix that is located at the cytoplasmic membrane interface, with one surface interacting with the membrane and/or transmembrane regions of the channel when the activation gates are in the open conformation. In this study, we sought to identify key S4S5 linker residues that interact with the membrane and/or intramembrane protein domains to stabilize the open state of the channel. A combination of nuclear magnetic resonance (NMR) and electrophysiology experiments suggest that the two tyrosine side chains (Tyr-542 and Tyr-545) point toward the membrane interior. Thermodynamic mutant cycle analysis experiments indicate that Tyr-542 interacts with both pore domain and voltage sensor residues to stabilize the open state of the channel. Conversely, mutations to Tyr-545 cause marked changes to the kinetics of deactivation but have much smaller effects on the voltage dependence of the distribution between open and closed states, which suggests that Tyr-545 stabilizes an intermediate transition state in the pathway between the activated and closed states.

Results
An alignment of the S4S5 linkers of different voltage-gated potassium channels is shown in Fig. 1. The most highly conserved residues are the charged arginine residues from the distal S4 and a glycine residue in the middle of the S4S5 linker. There are also two leucine residues in the S4S5 linker that are highly conserved in the Kv1-9 channels but are replaced by tyrosines in Kv10 -12 channels. In the crystal structure of the Kv1.2/2.1 chimera, the leucine residues in the S4S5 linker are oriented toward the lipid membrane. It has been shown previously that the S4S5 linker is critical for channel activation and deactivation in Kv11.1 (12,14,16,17) and that many of the S4S5 residues contribute to the slow gating kinetics. Previous studies, however, have largely focused on how the S4S5 linker interacts with cytoplasmic domains. For example, Asp-540 interacts with the bottom of the S6 transmembrane domain to stabilize the closed state of Kv11.1 channel (16), and cysteine scanning studies have suggested that the S4S5 domain interacts with the N-terminal tail (13). The tyrosine residues, however, have not been extensively studied. Our sequence homology analysis, coupled to the location of the S4S5 linker in crystal structures of the activated state, suggest that the two tyrosine residues are likely to contribute to interactions with the membrane and/or intra-membrane domains of the channel.
To determine whether the equivalent residues in Kv11.1 (Tyr-542 and Tyr-545) are likely to be oriented toward the lipid membrane, we studied the interaction of a peptide (Leu-532-Phe-551) that contains the S4S5 linker of the Kv11.1 channel with dimyristoylphosphatidylcholine (DMPC) or DMPC/dodecylphosphocholine (DPC) model membranes (Fig. 2). First, a  solution NMR was used to gain insight into whether the S4S5 linker interacts with the lipid membrane. We recorded onedimensional 1 H NMR spectra of the S4S5 peptides in water ( Fig.  2A, panel i) and in the presence of DMPC/DPC bicelles ( Fig. 2A,  panel ii). Comparison of the spectra shows changes in the chemical shifts and significant broadening of the peptide resonances, consistent with an interaction between the peptide and the model membranes (20). Phosphorus ( 31 P) and deuterium ( 2 H) solid-state NMR experiments were then employed to further investigate the interaction of the S4S5 peptide with the membrane using DMPC multilamellar vesicles (MLVs). Because phospholipids harbor a phosphate in their polar headgroup, 31 P NMR spectra were acquired to monitor change at the membrane surface. Broad power spectra are characteristic of non-oriented MLVs and represent the superposition of all chemical shifts from the distribution of orientations in the magnetic field (Fig. 2B, panel i), with the 90°edge on the right and 0°e dge on the left. To probe the bilayer's hydrophobic core, phospholipids with perdeuterated acyl chains (DMPC-d 54 ) were used. The spacing between each doublet on the 2 H NMR spectrum is known as the quadrupolar splitting (⌬ Q ) and gives insight on the degree of organization of the acyl chains (Fig. 2B, panel ii). The largest splitting corresponds to the CD 2 bonds closer to the headgroup, whereas the smaller splitting corresponds to the most mobile CD 3 group at the end of the lipid chain. Pure lipids 31 P (Fig. 2B, panel i, gray) and 2 H (Fig. 2B, panel ii, gray) spectra were recorded as a reference. The peptide disrupts model DMPC membranes suggesting the formation of fast-tumbling lipid assemblies, as revealed by the isotropic resonance on the 31 P and 2 H spectra (Fig. 2B, black). Moreover, the 31 P spectra for the peptide (Fig. 2B, panel i, black) shows a change in intensity between the 0 and 90°edges of model membranes (Fig. 2B, panel i, gray), with less lipids on the 90°edge, which can be ascribed to a change in vesicle shape or the transverse relaxation (T 2 ) effect caused by an increase in slow motions. The peptide mainly interacts with the surface of the membrane because no changes in ⌬ Q of 2 H spectra are observed, indicative of no deep insertion in the hydrophobic core of the membrane (21). Although these data suggest that the S4S5 linker is directly involved in membrane interaction, the orientation of the tyrosine residues was further probed by adding soluble paramagnetic gadolinium lanthanide ions (Gd 3ϩ ) to the samples. The Gd 3ϩ ions attenuate the proton signals of the peptide residues that are solvent-exposed (22). In agreement with solidstate NMR, the 1 H total correlation spectroscopy (TOCSY) spectrum of the linker (Fig. 2C) shows that the tyrosine ring protons are protected from the effect of Gd 3ϩ . This suggests that Tyr-542 and Tyr-545 are located facing the membrane.
We first performed alanine-scanning mutagenesis on the S4S5 residues (Asp-540 -Val-549) to investigate how individual residues in the S4S5 linker contribute to the voltage dependence of the distribution between closed and open states (activation) as well as the rate of activation of Kv11.1 channels.  Table 1. B, envelope of tail protocol to elicit the tail current at Ϫ70 mV after activation at 20 mV. One of the traces and its respective protocol are highlighted (panel i). The time constants of activation at 20 mV for WT and S4S5 mutants are summarized in panel ii. All summary data are presented as mean Ϯ S.E. and statistical comparisons made using one-way ANOVA with Bonferroni post-test. Significant differences (p Ͻ 0.05) are highlighted with asterisks. The values for all mutants are summarized in Table 1.
(102.5 Ϯ 6.1 ms, n ϭ 8) and mutants are summarized in Fig. 3B, panel ii, and in Table 1. Notably, both Y542A and Y545A did not significantly perturb the rate of activation at 20 mV compared with WT.
Previous studies have suggested that activation and deactivation are not simply reciprocal processes but rather occur via different pathways (18,23). We therefore investigated whether any of the S4S5 mutants preferentially altered the deactivation phenotype. Fig. 4A shows typical examples of 3-s tail current traces recorded at Ϫ70 mV (black) and Ϫ120 mV (gray) for WT (i), Y545A (ii), and Y542A (iii) following a 1-s step to 40 mV to fully activate the channels. The time constants for the fast component of deactivation in the voltage range from Ϫ70 to Ϫ150 mV for the three channels shown in Fig. 4B. The rates of deactivation of Y545A (in Fig. 4B, triangle) are 2-3 times faster than WT (circle) at all voltages measured. For Y542A (in Fig. 4B, square), the rate of deactivation is more than 10 times faster than WT at Ϫ70 mV, but at Ϫ150 mV there is almost no difference in the rates of deactivation for WT and Y542A channels. A summary of the time constant for the fast component of deactivation at Ϫ70 mV (black) and Ϫ120 mV (gray) for all S4S5 alanine mutants is shown in Fig. 4C and in Table 2. Fig. 4D shows a plot of the voltage dependence of the rates of deactivation for all the S4S5 alanine mutants. Mutation to the three charged residues (D540A, R541A, and E544A) all reduced the voltage dependence of the rates of deactivation. Conversely, mutations to the uncharged residues caused only modest changes to the voltage dependence of the rates of deactivation, with the notable exception of Y542A.
We next investigated how S4S5 mutants affected the voltage dependence of the distribution between open and closed states when starting from the open state, i.e. we measured 3-s isochronal deactivation curves to derive the V 0.5 of steady-state deactivation (24). Fig. 5A shows WT Kv11.1 channels recorded from a 3-s isochronal deactivation protocol. The midpoints of the voltage dependence of deactivation of V 0.5 derived by fitting the normalized amplitude of the Ϫ70 mV tail currents for WT (circle), Y542A (square), and Y545A (triangle) using Boltzmann function (Equation 1) were Ϫ61.3 Ϯ 0.5 mV (n ϭ 5), Ϫ1.6 Ϯ 1.1 mV (n ϭ 5), and Ϫ46.4 Ϯ 1.1 mV (n ϭ 5), respectively (Fig. 5B). Y542A caused the largest depolarizing shift in the V 0.5 of deactivation of all the S4S5 alanine mutants, whereas Y545A caused a more modest but still statistically significant depolarizing shift in the V 0.5 of deactivation (Fig. 5C).
Perturbations to the rates of deactivation in the S4S5 alanine mutants could be caused by either a change in the electrochemical driving force for deactivation and/or due to a change in the transition state energy barrier separating the open and closed states, i.e. there can be thermodynamic and/or kinetic components. To explore these two components, we have plotted the time constants for the fast component of deactivation at Ϫ120 and Ϫ70 mV versus the V 0.5 for 3-s isochronal deactivation for all S4S5 alanine mutants (Fig. 5, D and E). When kinetics were measured at Ϫ120 mV, most of the alanine mutants, with the notable exception of Y545A, reside within the 95% confidence interval obtained from a regression line fitted to the rate of deactivation versus isochronal deactivation V 0.5 data. This suggests that much of the change in the rates of deactivation at Ϫ120 mV for these S4S5 alanine mutants can be explained by changes in the thermodynamics of deactivation gating. Conversely, for Y545A, the rate of deactivation is much faster than would have expected from the simple change in the V 0.5 of deactivation. Consistent with the data at Ϫ120 mV, the change in the rates of deactivation at Ϫ70 mV for the S4S5 alanine mutants can also be explained by changes in the thermodynamics of deactivation gating, except for D540A and Y545A, whose mean values lie outside of the 5-95% confidence intervals. Note that S543A and G546A were excluded from analysis as these mutant channels did not deactivate at Ϫ70 mV. Our data show that whereas Y542A clearly affected the thermodynamics of deactivation gating, Y545A appeared to have a significant kinetic component on top of the alteration to the thermodynamics of deactivation gating.
As alanine mutagenesis of the two tyrosine residues appear to have quite distinct phenotypes in terms of effects on activation/ deactivation, both resulted in accelerated rates of deactivation with minimal perturbation on the rate of activation. We thus focus our investigation on the deactivation phenotypes of Tyr-545 and Tyr-542. Each tyrosine residue was mutated to Ala, Leu, Ile, Phe, or Trp. All of the Tyr-542 and Tyr-545 mutants resulted in accelerated rates of deactivation, with the exception of Y545F (Fig. 6A). Mutations at Tyr-542 tended to cause greater perturbations to the V 0.5 of deactivation than mutations at Tyr-545 (Fig. 6B). A plot of the perturbation to the rate of deactivation (measured at Ϫ70 mV) versus the change in V 0.5 of deactivation, for the groups of mutants at Tyr-542 and Tyr-545, highlights that mutations to Tyr-542 predominantly affect the thermodynamics of channel deactivation, whereas mutations to Tyr-545 cause relatively more perturbation to the kinetics of deactivation gating (Fig. 6C).

Discussion
Kv11.1 channels have unusual kinetics, which include rapid recovery from inactivation and slow deactivation. In the heart, the slow deactivation kinetics play a particularly important role in the suppression of propagation of ectopic beats that occur during the early diastolic interval (6). The S4S5 linker clearly plays a critical role in deactivation gating of voltage-gated K ϩ channels, including Kv11.1 channels (12,15,16). Here, we have shown that two tyrosine residues in the S4S5 linker, Tyr-542 and Tyr-545, play important, but distinct, roles in stabilizing the open state of the Kv11.1 channel and thus in mediating the slow deactivation gating of Kv11.1 channels.
In this study, mutations to Tyr-542 resulted in large shifts in the voltage dependence of the equilibrium distribution between the open and closed states, with all mutants shifting the equilibrium in favor of the closed state (Fig. 6B, panel i). All Tyr-542 mutants also resulted in an acceleration of the rates of deactivation (Fig. 6A, panel i). These results are consistent with a role for Tyr-542 in stabilizing the open state of the channel. Mutating Tyr-542 to alanine caused a marked acceleration of the rates of deactivation with a shallow slope (Fig. 4D). The fact that the slope of the voltage dependence of rates of deactivation for Y542A is much shallower than that of WT or Y545A suggests that the native Tyr-542 interacts with, and so influences, the motion of the VSD during deactivation. This conclusion is consistent with our mutant cycle analysis experiments indicating an energetic interaction between Tyr-542 and Arg-537, as well as with Val-535 and Ile-560 (Fig. 7), although only the interaction with Ile-560 exceeds the 1 kcal mol Ϫ1 threshold that is often used to indicate biologically relevant interactions (25)(26)(27)(28)(29)(30).
In contrast to Tyr-542, mutations to Tyr-545 caused less perturbation to the voltage dependence of the equilibrium distribution between the open and closed states (Fig. 6B, panel ii) but still caused a marked acceleration of the rates of deactivation (Fig. 6A, panel ii). This suggests that Tyr-545 might form  Table 2.
important interactions that stabilize the transition state complex; however, these interactions are destabilized by mutations, i.e. lower the energy barrier for channel deactivation. This suggestion is also consistent with the thermodynamic mutant cycle analysis experiments, which showed that there were no significant energetic interactions involving Tyr-545 and other residues in the stable end states. However, at this stage we cannot determine the nature of such interactions. They could be with transmembrane domain residues and/or with the membrane lipids. Given that the NMR experiments suggest that the peptide mainly interacts with the surface of the membrane (because there were no changes in ⌬ Q of 2 H spectra), if there are any specific interactions between Tyr-545 and the membrane, these would most likely involve interactions with the lipid headgroups.
The S4S5 linker is thought to be essential to transmit conformational changes of the voltage sensor to the pore domain, as evidenced by its proximity to the S6 helix in the crystal structure of the Kv1.2 channel, as well as its direct interaction with the S6 domain (16). Previous studies of the S4S5 linker in Kv11.1 channels have focused on its role in linking the voltage sensor domain to the activation gate at the intracellular end of the pore domain (11,16), as well as its role in mediating interactions with cytoplasmic domains that can also modulate deactivation gating (13). Our NMR studies of a peptide that contains the Kv11.1 S4S5 linker in a model membrane indicate that both Tyr-542 and Tyr-545 are likely to point toward the membrane, which is consistent with the homology model based on the structure of the activated Kv1.2/2.1 channel (see Fig. 7A). This suggests that the S4S5 linker could also have important interactions with the membrane directly and/or with residues in the transmembrane domains of the channel. We have clear evidence that Tyr-542 interacts with other transmembrane domain residues. However, the relatively small magnitude of this interaction (e.g. Ϫ4.6 kJ mol Ϫ1 for Y542A/I560A) suggests that other factors might be involved, which could include interactions between Tyr-542 and the membrane. Tyr-545 interactions with the membrane contributing to the stabilization of the transition state in the pathway connecting the activated state with the closed states are also plausible, but we do not yet have definitive proof of such interactions.
Recently, Lorinczi et al. (18) showed that cutting the S4S5 linker of Kv11.1 channels in half could still produce functional channels with normal activation but accelerated deactivation. This indicates that the major role of S4S5 linker is to stabilize the open state, and it is conceivable that the split S4S5 linker in the Kv11.1 prevents the interaction promoted by the S4S5 linker with other parts of the channel in the activated state. It is also possible that the split S4S5 may disrupt its interaction with the lipid membrane. It is unlikely that a split S4S5 linker would maintain a continuous helix that lies parallel to the membrane in the activated state of the channel (31)(32)(33)(34)(35)(36). Our functional studies provide complementary information that is consistent with recent findings (18) by showing that the two conserved tyrosine residues play important but distinct roles in stabilizing the open state of Kv11.1 channels.

NMR Spectroscopy
Materials-Peptide that contains the S4S5 linker (LVR-VARKLDRYSEYGAVLF) was synthesized by GL Biochem Ltd. (Shanghai, China) with Ͼ98% purity. Protonated and deuterated dimyristoylphosphatidylcholine (DMPC, DMPC-d 54 ) were purchased from Avanti Polar Lipids (Alabaster, AL), and deuterated dodecylphosphocholine (DPC-d 38 ), deuterium-depleted water, and diethylenetriaminepentaacetic acid gadolinium (III) dihydrogen salt hydrate were purchased from Sigma. Deuterium oxide (D 2 O) was obtained from CDN Isotopes (Pointe-Claire, Quebec, Canada). Phosphatidylcholines are abundantly found in biological membranes; consequently, DMPC is frequently used in model membranes, although DPC has been previously used to determine the Kv11.1 S4S5 linker structure by NMR (12). DMPC/DPC bicelles are a more biologically relevant alternative and have recently been introduced as a new membrane mimicking system for the NMR analysis of membrane proteins (37).
NMR Sample Preparation-Bicelle and MLV samples used for solution-and solid-state NMR, respectively, were prepared by mixing freeze-dried DMPC and/or DPC in water. The mixture was submitted to a series of at least three freeze (liquid N 2 )/thaw (50°C)/vortex shaking cycles. In all experiments, the peptide concentration was kept between 2 and 5 mM, and the lipid concentration was maintained well above the critical micelle concentration of the lipids. To ensure the stability of the model membranes, the hydration percentages were kept at 75% for solid-state NMR and 90% for high resolution NMR using deuterium-depleted water and 10% D 2 O, respectively. Solution

TABLE 2
Summary of the values for the 3-s isochronal deactivation V 0.5 , time constants of deactivation at ؊70 and ؊120 mV for S4S5 alanine mutants Data are presented as mean Ϯ S.E., and n denotes number of individual recordings. Asterisk indicates p Ͻ 0.05 versus WT using one-way ANOVA with Bonferroni post-test.

Mutant (n)
Deactivation V 0.5 (mV) NMR experiments used lipid molar ratios (q ϭ DMPC/DPC) of 0.25 and 1, as well as a lipid/peptide (L/P) molar ratio ranging from of 50:1 to 58:1. Solid-state NMR samples were prepared at a L/P molar ratio of 100:1. NMR Experiments-All experiments, except high resolution solution experiments, were recorded on a hybrid solution/ solid-state Varian Inova Unity 600 (Agilent Technologies, Mississauga, Ontario, Canada) spectrometer operating at frequencies of 599.95 MHz for 1 H, 246.86 MHz for 31 P, and 92.125 MHz for 2 H. Solid-state NMR experiments used a 4-mm broadband/ 1 H dual-frequency magic-angle spinning probehead. The 1 H chemical shifts were internally referenced by adding 0.5 mM of 2,2-dimethyl-2-silapentane-5-sulfonic acid set to 0.0 ppm, and 31 P NMR spectra were externally referenced with respect to the signal of 85% phosphoric acid set to 0 ppm. All data were processed using MestReNova (Mestrelab Research, Santiago de Compostela, Spain). 1 H total correlation spectroscopy (TOCSY) spectra were recorded on a hybrid solution/solid-state Bruker Avance III-HD narrow bore spectrometer (Milton, Ontario, Canada) operating at a frequency of 599.95 MHz for 1 H with a broad-  Table 2. D and E, relationship between time constants of deactivation for S4S5 alanine mutants at Ϫ120 and Ϫ70 mV, respectively, and their corresponding deacti-vation V 0.5 values from C. Linear regression was applied (straight line) with 95% confidence interval (CI) of the fit. Residues outside the 95% confidence interval are highlighted as filled circles. band double resonance 5-mm probe. A 90°pulse of 13 s and spin-lock time of 80 ms with MLEV17 spin-lock field of 10 kHz were employed (38). A spectral width of 6 kHz was used in both dimensions with 2048 and 256 complex data points in the direct and indirect dimensions, respectively. A total of 40 transients were accumulated with a repetition delay of 2 s. 31 P NMR spectra were recorded using a phase-cycled Hahn echo pulse sequence with gated broadband proton continuous wave decoupling at a field strength of 20 kHz (39). A 90°pulse length of 13 s was used with interpulse delays of 33 s. Typically 1024 scans were acquired with a recycle delay of 5 s. The acquisition time was set at 10 ms with a 5-s dwell time. 2 H NMR spectra were obtained using a solid-echo pulse sequence (40) with a 90°pulse length of 2.8 s, interpulse and repetition delays of 20 s and 0.5 s, respectively. At least 5000 data points were obtained, and typically 2048 scans were acquired with an acquisition time of 5 ms and 1-s dwell time.
Electrophysiology Methods-We chose to use Xenopus oocytes for our electrophysiological recordings as this system enables high level expression of almost all channels, including mutant channels that result in expression defects in mammalian systems (41).
Molecular Biology-The cDNA of Kv11.1 (a gift from Dr. Gail Robertson, University of Wisconsin) was subcloned into a pBluescript vector, which contains the 5Ј-untranslated region (UTR) and 3Ј UTR of the Xenopus laevis ␤-globin gene (a gift from Dr. Robert Vandenberg, University of Sydney). Site-directed mutagenesis of Kv11.1 cDNA was performed using the QuikChange method (Agilent Technologies, Mulgrave, Victoria, Australia), and mutations were confirmed by Sanger DNA sequencing. WT and mutant channel cDNAs were linearized with BamHI-HF (New England Biolabs, Ipswich, MA), and cRNA was transcribed with T7 RNA polymerase using the mMessage mMachine kit (Ambion, Austin, TX).  Oocyte Preparation-Female X. laevis frogs were purchased from Nasco (Fort Atkinson, WI). The Garvan/St Vincent's Animal Ethics Committee granted approval (ID 14/30) for all animal experiments described in this study. Frogs were anesthetized by immersion in 0.17% w/v Tricaine. One ovarian lobe was mobilized, and ϳ5 ml were removed. The follicular cell layer was digested by using 1 mg/ml collagenase A (Roche Applied Science) in Ca 2ϩfree ND96 solution containing (in mM): NaCl 96, KCl 2, MgCl 2 1.0 and Hepes 5 (pH adjusted to 7.5 with 5 M NaOH) for ϳ2 h. After rinsing with ND96 (as above, plus 1.8 mM CaCl 2 ) to remove collagenase A, stage V and VI oocytes were isolated and stored at 18°C in ND96 supplemented with 2.5 mM pyruvic acid sodium salt, 0.5 mM theophylline, and 50 g/ml gentamicin. X. laevis oocytes were injected with cRNA and incubated at 18°C for 24 -48 h prior to electrophysiological recordings.
Electrophysiology-Two-electrode voltage-clamp experiments were performed at room temperature (ϳ21°C) using a Geneclamp 500B amplifier (Molecular Devices Corp., Sunnyvale, CA). We used glass microelectrodes with tip resistances of 0.3-1.0 megohms, when filled with 3 M KCl. Oocytes were perfused with ND96 solution (see above) during all experiments. A 50-ms 20-mV depolarization step from the holding potential of Ϫ90 mV was applied at the start of each sweep to enable off-line leak-current subtraction. We assumed that the current leakage was linear in the voltage range Ϫ150 to ϩ40 mV. Data acquisition and analysis were performed using pCLAMP software (Version 10.2, Molecular Devices, Sunnyvale, CA), Excel software (Microsoft, Seattle), and Prism 6 (GraphPad Software Inc., La Jolla, CA). All parameter values are reported as means Ϯ S.E. of the mean (S.E.) for n experiments, where n denotes the number of different oocytes studied for each construct.
Data Analysis-To measure the voltage dependence of activation, cells were depolarized from Ϫ70 to ϩ60 mV for 3 s (exact voltage range depending on the mutant tested) followed by a step to Ϫ70 mV to measure tail current amplitude. Tail current amplitudes were normalized to the maximum tail current value and fitted with a Boltzmann function as shown in Equation 1, I/I max ϭ ͑1 ϩ e ͑v0.5Ϫvt͒/k ͒ Ϫ 1 (Eq. 1) where I/I max is the relative tail current amplitude; V 0.5 is the voltage at which 50% of channels are deactivated; V t is the test potential, and k is the slope factor. To measure rates of activation, an envelope of tail protocol was used to obtain the rates of activation at 20 mV (42). The peak tail current amplitudes recorded at Ϫ70 mV using different activation duration at 20 mV were fitted with a single exponential function to obtain a time constant for activation.
Rates of deactivation were measured from the hooked tail current traces recorded at voltages in the range Ϫ70 to Ϫ150 mV, following a 1-s depolarization step to 40 mV from a holding potential of Ϫ90 mV. The decaying phase of the current traces was fitted with a double exponential component to derive both components of channel deactivation. However, only the fast component is reported in this study. To measure the voltage dependence of deactivation, cells were depolarized to ϩ40 mV for 1 s to ensure channels were fully activated and then stepped to voltages between ϩ60 and Ϫ100 mV (exact voltage range depending on the mutant tested) for 3 s, followed by a step to Ϫ70 mV to measure tail current amplitude (23,24,43). Tail current amplitudes were normalized to the maximum tail current value and fitted with a Boltzmann function using Equation 1. In experiments in which the free energies associated with voltage-dependent deactivation were examined, the same datasets were fitted with the thermodynamic form of the Boltzmann distribution as shown in Equation 2, where ⌬G 0 is the work done at 0 mV; z g is the effective number of gating charges moving across the membrane electric field (E); F is the Faraday constant; R is the universal gas constant; and T is the absolute temperature. From Equation 2 we can also calculate the effect of a mutation on ⌬G 0 as shown in Equation 3, Double Mutant Cycle Analysis-To investigate whether two residues interact to stabilize the open and/or closed states we used thermodynamic mutant cycle analysis. The interaction energy, ⌬⌬G 0 int , between two mutants was calculated as the difference in free energy of the double mutant compared with the sum of the free energy differences of the individual mutants shown in Equation 4, Absolute ⌬⌬G 0 int values of greater than 4.2 kJ mol Ϫ1 (i.e. 1 kcal mol Ϫ1 ) were considered significant (25)(26)(27)(28)(29)(30).
Statistical Analysis-We used one-way ANOVA with Bonferroni post-test analysis to compare rates of deactivation and V 0.5 of deactivation or an unpaired t test when comparing the ⌬⌬G 0 values of double mutants (⌬⌬G 0 Mut1 ϩ Mut2 ) compared with the sum of ⌬⌬G 0 values of individual mutants (⌬⌬G 0 Mut1 ϩ ⌬⌬G 0 Mut2 ) in double mutant cycle analysis experiments. p values Ͻ0.05 were considered significant.
Homology Modeling-The crystal structure of a Kv1.2/2.1 channel chimera (19) was used to generate the homology model of Kv11.1 using Swiss PdbViewer (44), and the final structure was optimized by using SWISS-MODEL Workspace (45,46).