|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
MEMBRANE TRANSPORTERS, ION CHANNELS, AND PUMPS
1Laboratory for Molecular Biophysics, Physiology, and Pharmacology, Department of Biomedical Sciences, University of Antwerp, Belgium; and 2Laboratory for Molecular Biophysics, Department of Biochemistry, University of Oxford, Oxford, United Kingdom
Submitted 11 October 2007 ; accepted in final form 26 March 2008
| ABSTRACT |
|---|
|
|
|---|
voltage-gated potassium channel
-subunits that surround a central ion-conducting pore. Each
-subunit has six membrane-spanning helices (S1–S6) with a pore loop (P) between S5 and S6. The detailed crystal structure of a voltage-independent bacterial K+ channel, KcsA, has resolved the molecular mechanism of K+ selectivity and permeation (14, 36, 51). The crystal structure of the mammalian voltage-gated Kv1.2 channel demonstrated that the structure of the ion-conducting pore (S5-pore loop-S6) is similar to that of KcsA (31, 32). Thus the molecular mechanism of K+ permeation is conserved between KcsA and Kv channels. The opening and closing of the ion-conducting pore (gating process) of voltage-gated channels is triggered by a change in the membrane potential. The cytoplasmic activation gate of Shaker-type channels [of which hKv1.5 is a human homolog (44)], which open or close the ion-conducting pore, is located in the inner section of the ion permeation pathway. It has been proposed that residues V478 and F481 in Shaker (V516 and F519 in hKv1.5, respectively) occlude the pore in the closed state (11, 16). However, it remains controversial which residues form the gating "hinge" or pivoting point. Extrapolation of the predicted gating mechanisms of the Kv channel from archeabacterium Aeropyrum pernix (KvAP), KcsA, and the Ca2+-gated K+ channel from Methanobacterium autotrophicum (MthK) to Shaker suggested that G466 in Shaker (G504 in hKv1.5; Fig. 1) would act as hinge or pivoting point (24–26, 49). However, the crystal structure of bacterial K+ inward-rectifying channel (KirBac1.1) (28) suggested that G143 (the equivalent of P475 in Shaker and P513 in hKv1.5) was the pivoting point for gating, whereas the former more central-located glycine in the transmembrane helix G134 (corresponding to G466 in Shaker) was more important for packing. In human ether-a-go-go-related gene K+ channels (hERG), both the central and lower-located glycine in S6 are also required for tight protein packing, since the S6 helix is already inherently flexible (17). A substitution scan in K+ inward-rectifying (Kir1.1) channels showed that both glycine residues facilitate gating in a similar manner but that they are not absolutely required for pH-dependent gating of this class of channels (40).
|
-helical configuration of S6 (10, 29, 48). The crystal structure of the voltage-dependent rat (r) Kv1.2 channel showed that hinging could occur in the vicinity of the central glycine, but the structure also showed that the S6 helix is indeed bent around the PXP motif (31) (Fig. 1). Apparently this
-helical disruption is crucial for voltage-dependent gating of this class of channels, but the proline-containing motif does not need to be located at the evolutionary conserved position (29, 48). Indeed, hKv1.5 channels with an altered PXP sequence, PVAV, AVPP, or AVAP (Fig. 1), were functional but displayed marked changes in channel gating. These channels activated very slowly, and the voltage dependence of activation did not reach saturation during 5-s voltage steps up to +110 mV (29). Here we report that channels with a AVPP sequence display two gating modes. | MATERIALS AND METHODS |
|---|
|
|
|---|
Starting from hKv1.5 cDNA cloned in a pBK-CMV vector, mutations were introduced in hKv1.5 with the QuikChange Site-Directed Mutagenesis kit (Stratagene, La Jolla, CA) using mutant primers. After the PCR-based mutagenesis, a Pst1-BamH1 fragment containing the mutation was cut out of the PCR-amplified vector and ligated in hKv1.5/pBK-CMV to replace the wild-type (WT) sequence. Double-stranded sequencing of the exchanged fragment and the adjacent sequence confirmed the presence of the desired modification and the absence of unwanted mutations.
Cell Culture and Transfection
Plasmid DNA for mammalian expression was obtained by amplification in XL2 blue script cells (Stratagene). The plasmid DNA was isolated from the bacterial cells with the GenElute HP plasmid maxiprep kit (Sigma, St. Louis, MO). The cDNA concentration was determined with ultraviolet absorption. Mouse Ltk– cells were cultured in DMEM medium supplemented with 10% horse serum and 1% penicillin/streptomycin. The culture was passed every 3–4 days. For experiments, the cells were transfected with 0.5–5 µg of WT or mutant subunit cDNA. Transfections were done according to the lipofection method using lipofectamine (Invitrogen, Carlsbad, CA). After transfection (12–24 h), the cells were trypsinized and used for analysis within 12 h. A stable cell line was created for the mutant with an AVPP sequence as described previously (43). For this purpose, the mutant AVPP was cloned in a pMSV-neo vector with an inducible promoter. For experiments, channel expression was induced by adding 100 µl dexamethasone (of a 100 µM stock solution) to the cell medium 16–24 h before patch-clamp analysis.
Electrophysiology
Current recordings were made with an Axopatch-200B amplifier (Axon instruments, Foster City, CA) in the whole cell configuration of the patch-clamp technique. Experiments were done at room temperature (20–23°C); current recordings were low pass filtered and sampled at 2–10 kHz with a digidata 1200A data acquisition system (Axon instruments). Command voltages and data storage were controlled with pClamp8 software (Axon Instruments). Patch pipettes were pulled from 1.2 mm kwik-fill borosilicate glass capillaries (World Precision Instruments, Sarasota, FL) with a P-2000 puller (Sutter Instruments, Novato, CA). After heat polishing, the resistance of the patch pipettes was <3 M
. The cells were perfused continuously with a bath solution containing (in mM) 130 NaCl, 4 KCl, 1.8 CaCl2, 1 MgCl2, 10 HEPES, and 10 glucose, adjusted to pH 7.35 with NaOH. The pipettes were filled with intracellular solution containing (in mM) 110 KCl, 5 K4BAPTA, 5 K2ATP, 1 MgCl2, and 10 HEPES and was adjusted to pH 7.2 using KOH. Junction potentials were zeroed with the filled pipette electrode in the bath solution. Experiments were excluded from analysis if the voltage errors originating from series resistance exceeded 5 mV.
The holding potential was –80 mV unless specified otherwise. The voltage protocols were adapted to adequately determine the different biophysical parameters of mutant channels. For WT hKv1.5, a 900-ms voltage ramp from –80 to +100 mV was optimal. For AVPP, a 1,200-ms ramp from –80 to +100 mV was used, since ramps longer than 1,200 ms were blurred by the "slow" activation state. Because of the Goldman-Hodgkin-Katz (GHK)-type rectification of hKv1.5 (43), the current recordings from ramp protocols were normalized with the GHK current equation;
![]() | (1) |
Molecular Dynamics Simulations
A series of simulations (each 10 ns long) of the isolated S6 helix of hKv1.5 K+ channel, modeled on the closed and open conformation, was carried out in the presence of a mixed 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphorylglycerol bilayer. The following systems have been set up: WT hKv1.5 and the PXP motif variants AVAP and AVPP. The initial protein coordinates for the closed conformation were modeled using Modeller and PDB entry 1K4C [PDB] as a template. The coordinates for the open conformation were extracted from the rKV1.2 X-ray structure (PDB entry: 1P7B [PDB] , www.rcsb.org). The transmembrane S6 helix was defined as extending from residue Val491 to Thr527, the PXP motif encompassing residue 511–513. The system (i.e., protein plus membrane) was solvated with SPC water molecules (4), retaining all the crystallographic waters.
Molecular dynamics (MD) simulations were performed with GROMACS 3.1.4 (30) (www.gromacs.org) with a modified version of the GROMOS-87 force field (47). Lipid parameters were based on those by Berger et al. (5) and Marrink et al. (35). For the lipid-protein interactions, the GROMOS parameters were used. Simulations were carried out at constant pressure and temperature, NPT statistical ensemble, with periodic boundary conditions. The initial velocities were taken randomly from a Maxwellian distribution at 300 K. The temperature was held constant by coupling to an external bath (22). Long-range electrostatic interactions were calculated using the Particle Mesh Ewald summation methods (9). Lennard-Jones interactions were calculated using a cutoff of 0.9 nm. The pair lists were updated every 10 steps. The LINCS algorithm (20) was used to constrain bond lengths. The time step was 2 fs, and coordinates were saved every 0.1 ps. Principal components analysis was performed as previously described to (1, 15). The computational Markov modeling is described in the supplemental data (Supplemental data for this article is available online at the American Journal of Physiology: Cell Physiology website.).
| RESULTS |
|---|
|
|
|---|
-helix by substituting alanine residues for both prolines (P511A + P513A) also resulted in a nonfunctional channel. However, both mutants could be rescued by the introduction of a proline at a position different from the conserved one. The double mutant P511A + V514P resulted in an AVPP sequence (Fig. 1) and rescued the P511A mutation, since voltage-dependent K+ currents were observed after patch disruption. Upon the first depolarization to +60 mV from a holding potential of –80 mV, the current activated quickly with time constants in the millisecond range and subsequently inactivated slowly (Fig. 2C). This time course of current activation was similar to WT hKv1.5 (Fig. 2A) and to that of the V514P mutation in WT, resulting in a PVPP sequence (Figs. 1 and 2B). However, during the second depolarization to +60 mV, the current activation was markedly slowed (s range) and displayed no time-dependent inactivation (Fig. 2C). This second type of activation (Fig. 2D) was similar to that of channels with an AVAP (triple mutant P511A + P513A + V514P; Fig. 2E) and a PVAV sequence (mutant P513A, Fig. 2F), all displaying markedly altered channel gating.
|
200 times slower than WT (PVPV) or the control mutant V514P (Fig. 3B). The time constant of activation of AVPP in the initial gating mode was 10.9 ± 1.0 ms (n = 8) at +60 mV, which is approximately three times slower than WT or the control mutant V514P (PVPP) but markedly (
55 times) faster than that of its second gating mode (Fig. 3B). Throughout we will use "fast" and "slow" mode as a shorthand for the initial gating mode that is characterized by fast WT-like activation or the second gating mode that shows slow current activation, respectively. The subsequent channel inactivation of the fast mode had kinetics that displayed no or a limited voltage dependency and were similar to those of WT hKv1.5 (Supplemental Fig. S1). Approximating the fast mode inactivation of AVPP with a single exponential function yielded a time constant of 702 ± 46 ms (n = 8) at +30 mV. The deactivation time constant of AVPP in the slow mode was comparable to that of the control mutant V514P with a PVPP sequence (Fig. 3B); both were approximately six times slower than WT (PVPV). All mutants had a reversal potential around –80 mV (data not shown), indicating that the K+ selectivity of the channels was not altered.
|
43% of the current density in this cell line. The expected values for the subsequent second and third depolarizing steps were 25 and 12% of the total value of 112 ± 11 pA/pF, being 28 ± 3 and 14 ± 1 pA/pF respectively. The measured values for this second and third step also compared very well to the expected ones (Fig. 4F). Furthermore, the sum of the current amplitudes obtained during all three depolarizing steps accounted for 90% of the activatable channels and therefore approached the current density value of this cell line (93 ± 15 pA/pF, n = 11; Fig. 4F). These results indicate that the activation (or deactivation) and not the subsequent inactivation trapped the channel in the slow mode; this trapping was observed at any depolarizing voltage above the apparent threshold (–10 mV; supplemental Fig. S2).
|
|
|
-helix, MD simulations were performed on the S6 transmembrane domain of WT hKv1.5 and the mutant channels AVPP and AVAP in a membrane-mimicking environment. For both WT and mutant channels, we ran three independent simulations starting from either a closed or open channel conformation, each time using a new set of initial velocities obtained at 300 K (see MATERIALS AND METHODS for details). Analysis of the stability and reliability of all simulations (based on calculation of root mean square deviation) showed that the isolated S6 helix of all systems was quite mobile and subjected to short time scale rearrangements (
1, 1.5 ns) (data not shown). Secondary structure content analysis was performed using the Dictionary of Secondary Structures of Proteins software (27 and Fig. 7). The analysis revealed that, in the simulations starting from the closed conformation (left part in Fig. 7), the
-helical content for WT was lost in the middle of the helix around residue 511, which corresponds to the first proline of the PXP motif. This is an indication for the presence of a hinge point around this position. For the AVAP variant, the loss of
-helical structure was more pronounced, and the COOH-terminal half of the helix adopted either a p-helix or turn conformation, which represents a different dynamic behavior for the AVAP mutant compared with WT. On the other hand, the AVPP mutant displayed an average
-helical content per residue similar to WT and displayed the marked loss of a structure and the gain of a turn-type conformation around residues 510–514 comparable to WT.
|
-helical content at position 511 that was clearly evident. This loss in
-helical content was absent in both the AVAP and the AVPP mutant. Interestingly, both showed a quite similar secondary structure content in the open helix simulations, suggesting that the open state of the S6 helix of these mutants is more stable than WT on this time scale. A further indication of this trend is shown in the analysis of fluctuations performed by measuring the mean square displacement from the average structure according to the first principal components of the covariance fluctuations matrix. Figure 8 clearly shows that, for WT, in both closed and open conformations, there was a clearly higher displacement located near residue 509 over the entire time scale of the simulation. Because these fluctuations were focused around the PXP motif, this region of the helix was less stable and likewise more flexible. Starting from a closed conformation, the AVPP variant displayed a higher flexibility around the residues 510–511, a behavior qualitatively comparable to the WT counterpart. However, in the open helix conformation, the flexibility pattern according to the principal components suggested that the intrinsic dynamical behavior of AVPP matched that of AVAP, indicating that the movements of the isolated S6 helix differ from that of WT.
|
| DISCUSSION |
|---|
|
|
|---|
-helix) resulted in a nonfunctional channel; even replacing the first proline by an alanine, creating the sequence AVPV, resulted in a nonfunctional channel. However, both channel mutants reached the plasma membrane, as assessed by confocal microscopy (29). This was in agreement with the results for the corresponding mutant in Shaker (P473A) that displayed clear gating currents but no ionic current (16). However, the introduction of an extra proline at position 514 (mutation V514P), creating a AVPP sequence, rescued the nonfunctionality of P511A. This AVPP sequence resulted in a channel that existed initially in a fast gating mode (Figs. 2C and 4), but upon the first activation or deactivation the channel became trapped in a slow one (Fig. 2, C and D, and 3). The gating kinetics of this slow mode were comparable to those of channels with a AVAP or a PVAV sequence (Fig. 2, E and F). However, these AVAP and PVAV channels did not display two gating modes, since the current activation upon the very first depolarization was the same as the second (data not shown). The biophysical properties of these channels and of AVPP in the slow mode were profoundly altered compared with WT hKv1.5 (PVPV). Because very short depolarizing steps were sufficient for AVPP channels to switch to their slow gating mode (Fig. 4, D–F), it was presumably not the inactivation that caused the switch. The fast gating mode of AVPP channels could be recovered, but this recovery was extremely slow (time constants in the order of 10 min) and voltage-independent for holding potentials between –50 and –80 mV (Fig. 6). Because of this very slow recovery process, we could not characterize the fast mode of AVPP completely. With the approximation of using a voltage ramp (Fig. 5C), the estimated voltage dependence of the fast mode was much closer to that of WT compared with that of the slow mode, which itself did not saturate (Fig. 3A). Furthermore, the activation time constant of the fast mode corresponded well to that of WT channels (PVPV) and the control mutant V514P with a PVPP sequence, especially if corrected for midpoint shifts. Remarkably, the deactivation time constants of AVPP fast and slow mode were similar and resembled those of the control mutant V514P, and both were markedly slower than WT (Fig. 3B). Such slowing of channel deactivation was also observed with alanine or isoleucine substitutions at this position (38). These V514A and V514I mutations displayed an isolated "deactivation" failure that was caused by stabilizing the open state while destabilizing the closed state. Such a slowing of deactivation is apparently common for mutations at this position given the results for the V514P mutant (PVPP) and the AVPP channel.
Previous results have suggested that the proline residues of the tandem proline motif (PVPV) disrupt the S6
-helix (10, 29, 48). The question remained whether this putative disruption of the
-helix is required solely for structural packing of the channel protein or whether it also serves a functional role by creating a gating hinge. If the PXP motif disrupts ("bends") the S6
-helix, this region should have greater degrees of freedom and could therefore generate flexibility. This idea is supported by molecular modeling that indicated that the PXP motif distorted the
-helix, creating a hinge region (6, 7, 33, 45, 46). If the PXP motif is important in WT to reorient the lower part of S6 (channel gate) during gate opening or closing, the same most likely accounts for the mutants. The biophysical data indicated that the Kv channel mutant with a AVPP sequence initially exists with the lower S6 region in a conformation that allows "normal" activation, but activation or subsequent deactivation locks S6 in a different conformation that precludes normal activation. MD simulations of the S6 helix of WT hKv1.5 revealed significant movements around the PVP motif in both the open and closed conformation. However, for the AVAP sequence, the simulations suggested that the helix is rigid in both the open and closed conformation, since the principal component analysis identified no specific movements (Fig. 8). The absence of clear movements may explain the slow kinetics of this AVAP mutant. Interestingly, the MD simulations of AVPP showed that this mutant displayed two different behaviors (as in the functional analysis). When starting from the open conformation, the S6 helix of the AVPP mutant was as rigid as AVAP, which could be correlated with the slow gating kinetics. In contrast, when starting from a closed conformation, AVPP displayed movements similar to WT, which suggests the presence of a flexible hinge and may explain the initial observed WT-like kinetics. This indicates that the AVPP sequence can either introduce flexibility in the S6 helix or be rigid, dependent on starting from a closed or open channel conformation, respectively. Because this switch in behavior occurs during channel gating, this indicates that this region (and presumably also the PVPV region in WT) undergoes conformational movements during channel opening or closing and could therefore act as a hinge.
Olcese et al. (37) described a double mutation in Shaker that converted the inactivated state into a second open state. One of these mutations, I470C (I508C in Kv1.5), is in the vicinity of the PXP motif. This raises the question whether the double-gating behavior of the channel with an AVPP sequence can be explained by a modified inactivated state. In this case, the effect of the AVPP mutation would not (only) be at the level of the S6 "bundle crossing" but also at the selectivity filter. Under this hypothesis, the reduced flexibility of the S6 helix with simulations starting from the open conformation would indicate that the cytoplasmic gate becomes locked in the open configuration after the first activation. The observed second slow gating mode would then be due to an alteration in the selectivity filter that converts from a new nonconducting state to the former "inactivated" state but which is now conducting. However, the inactivation kinetics of the fast mode are largely voltage independent and clearly slower than the activation time constants of the slow mode (supplemental Fig. S1). Furthermore, because the channel still deactivates, this would require that the normal-conducting (open) state of the selectivity filter becomes now a nonconducting one since the cytoplasmic gate is locked open and cannot seal off the channel pore anymore. Although we cannot completely rule out this possibility, it relies on several additional assumptions. Therefore, we prefer the proposal that the slow mode is triggered when the cytoplasmic S6 gate closes into a configuration that precludes fast activation.
Computational models to simulate the AVPP behavior indicated that a switch in gating mode at the level of an intermediate closed state matched the experimental data best (Supplemental Fig. S4).
|
| (2) |
If the switching is located closer to the open state, the behavior during depolarization can still be simulated, but not the relatively fast deactivation combined with slow activation kinetics (Supplemental Fig. S4). The observation that the deactivation after the first (fast) and subsequent (slow) activation proceeds with the same time course can be interpreted that the initial closing step(s) is the same in both gating modes but that a further reorientation upon deactivation shifts the S6 gate in a conformation from which subsequent activation is slow. Thus, if we approach the behavior of the AVPP mutant with a simplified version of the 15 or 16 state model (50), the state transitions of one subunit of the tetramer would schematically be as follow.
|
| (3) |
In equations 2 and 3, the A to O transition refers to the final concerted step of opening and presumably involves a movement(s) of the complete tetramer. Computational modeling suggested (Eq. 2, Supplemental Fig. S4) that the channel switches to a reluctant closed state (C*) at an intermediate closed state (Ci). This C* branch can thus be seen as a mirror of the normal closed C branch. The recovery process from this reluctant closed state C* to the normal C state is then extremely slow. If we combine this equation 2 with the results of the MD simulations, the normal "WT" closed C state is defined by a S6 segment that shows flexibility around the PVP motif and allows both fast channel activation and deactivation from and to this state or configuration. The reluctant C* branch on the other hand is characterized by a S6 segment that is rigid and results in slow gate opening. The results of the MD simulations are in good agreement with this interpretation, since the S6 helix revealed flexibility starting from the closed conformation but not when the channel was in the open state. Thus once the open or preopen activated conformation is reached, the S6 segment becomes less flexible, and the channel switches to the reluctant C* branch upon channel closure.
The cytoplasmic activation gate of Shaker-type channels is formed by the lower part of S6 (11, 12, 16, 38). In channels with a AVPP sequence, this lower part of S6 switches between a nearly normal fast conformation and an altered slow conformation. Although we cannot exclude the possibility that another associated channel segment locked the cytoplasmic activation gate in this new position, it is likely that the altered gating hinge reoriented and locked the channel in a slow mode. Presumably, channels with a PVAV or AVAP sequence instead of WT (PVPV) exist in a similar but permanent slow conformation, since MD simulations showed lack of flexibility in both conditions. One could speculate that the association with the KCNE1 subunit results in a similar permanent slow conformation in the KCNQ1 subunits, resulting in the known slow activation kinetics of the cardiac slow component of the delayed-rectifier current (IKs) (2, 42). Similarly, N-type and P/Q-type calcium channels switch from a rapidly activating "willing" (W) mode to a more-difficult-to-activate reluctant (R) mode upon binding of Gβ
subunits (3, 18, 23, 41). These channels can be shifted back from the R mode toward the W mode by stronger depolarizations or by protein kinase C phosphorylation. Apparently this R mode and its different modulation is an intrinsic property of the calcium channel itself (19). This indicates that, similar to our hKv1.5 channel with an AVPP sequence, these calcium channels reversibly switch between a fast activating W mode and a slow or R mode.
In conclusion, our results indicate that changing the evolutionary conserved PVPV motif in hKv1.5 into AVPP resulted in a fully folded channel that resides in the membrane and reversibly switches between a fast and a slow gating mode. This indicates that the region around the double proline (PXP) sequence reorients upon channel activation and/or deactivation. MD simulations of the isolated S6 helix show that the AVPP sequence can introduce flexibility in this region (fast mode) but at the same time being very rigid (slow mode). Because only a minor alteration around the crucial PXP motif may lead to faster or slower channel gating, the large variety in gating kinetics observed throughout the Kv family of channels may result from only subtle differences in the gating machinery that transduces voltage sensing to channel opening.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
Current address for A. Grottesi: Consorzio Interuniversitario per le Applicazioni del Supercalcolo per Universita e Ricerca, Rome, Italy.
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
2. Barhanin J, Lesage F, Guillemare E, Fink M, Lazdunski M, Romey G. K(v)LQT1 and IsK (minK) proteins associate to form the IKs cardiac potassium current. Nature 384: 78–80, 1996.[CrossRef][Medline]
3. Bean BP. Neurotransmitter inhibition of neuronal calcium currents by changes in channel voltage dependence. Nature 340: 153–156, 1989.[CrossRef][Medline]
4. Berendsen HJC, Postma JPM, van Gunsteren WF, Hermans J. Interaction models for water in relation to protein hydration. In: Intermolecular Forces. Dordrecht, The Netherlands: Reidel, 1981, p. 331–342.
5. Berger O, Edholm O, Jahnig F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys J 72: 2002–2013, 1997.[Web of Science][Medline]
6. Bright JN, Sansom MSP. Kv channel S6 helix as a molecular switch: simulation studies. IEE Proc Nanobiotechnol 151: 17–27, 2004.[CrossRef][Medline]
7. Bright JN, Shrivastava IH, Cordes FS, Sansom MS. Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies. Biopolymers 64: 303–313, 2002.[CrossRef][Web of Science][Medline]
8. Clay JR. Determining K+ channel activation curves from K+ channel currents. Eur Biophys J 29: 555–557, 2000.[CrossRef][Web of Science][Medline]
9. Darden T, York D, Pedersen L. Particle Mesh Ewald: an N. log(N) method for ewald sums in large systems. J Chem Phys 98: 10089–10092, 1993.[CrossRef]
10. del Camino D, Holmgren M, Liu Y, Yellen G. Blocker protection in the pore of a voltage-gated K+ channel and its structural implications. Nature 403: 321–325, 2000.[CrossRef][Medline]
11. del Camino D, Yellen G. Tight steric closure at the intracellular activation gate of a voltage-gated k(+) channel. Neuron 32: 649–656, 2001.[CrossRef][Web of Science][Medline]
12. Ding S, Horn R. Tail end of the s6 segment: role in permeation in shaker potassium channels. J Gen Physiol 120: 87–97, 2002.
13. Ding S, Ingleby L, Ahern CA, Horn R. Investigating the putative glycine hinge in shaker potassium channel. J Gen Physiol 126: 213–226, 2005.
14. Doyle DA, Cabral JM, Pfuetzner RA, Kuo A, Gulbis JM, Cohen SL, Chait BT, MacKinnon R. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science 280: 69–77, 1998.
15. Garcia AE. Large-Amplitude Nonlinear Motions in Proteins. Phys Rev Lett 68: 2696–2699, 1992.[CrossRef][Web of Science][Medline]
16. Hackos DH, Chang TH, Swartz KJ. Scanning the intracellular S6 activation gate in the Shaker k(+) channel. J Gen Physiol 119: 521–532, 2002.
17. Hardman RM, Stansfeld PJ, Dalibalta S, Sutcliffe MJ, Mitcheson JS. Activation gating of hERG potassium channels: S6 glycines are not required as gating hinges. J Biol Chem 282: 31972–31981, 2007.
18. Herlitze S, Garcia DE, Mackie K, Hille B, Scheuer T, Catterall WA. Modulation of Ca2+ channels by G-protein beta-gamma subunits. Nature 380: 258–262, 1996.[CrossRef][Medline]
19. Herlitze S, Zhong H, Scheuer T, Catterall WA. Allosteric modulation of Ca2+ channels by G proteins, voltage-dependent facilitation, protein kinase C, and Ca(v)beta subunits. Proc Natl Acad Sci USA 98: 4699–4704, 2001.
20. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM. LINCS: a linear constraint solver for molecular simulations. J Comp Chem 18: 1463–1472, 1997.[CrossRef]
21. Hille B. Ion Channels of Excitable Mmembranes. Sunderland, MA: Sinauer, 2001.
22. Hoover WG. Canonical dynamics: equilibrium phase-space distributions. Phys Rev A 31: 1695–1697, 1985.[CrossRef][Medline]
23. Ikeda SR. Voltage-dependent modulation of N-type calcium channels by G-protein beta gamma subunits. Nature 380: 255–258, 1996.[CrossRef][Medline]
24. Jiang Y, Lee A, Chen J, Cadene M, Chait BT, MacKinnon R. Crystal structure and mechanism of a calcium-gated potassium channel. Nature 417: 515–522, 2002.[CrossRef][Medline]
25. Jiang Y, Lee A, Chen J, Cadene M, Chait BT, MacKinnon R. The open pore conformation of potassium channels. Nature 417: 523–526, 2002.[CrossRef][Medline]
26. Jiang Y, Lee A, Chen J, Ruta V, Cadene M, Chait BT, MacKinnon R. X-ray structure of a voltage-dependent K+ channel. Nature 423: 33–41, 2003.[CrossRef][Medline]
27. Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 22: 2577–2637, 1983.[CrossRef][Web of Science][Medline]
28. Kuo A, Gulbis JM, Antcliff JF, Rahman T, Lowe ED, Zimmer J, Cuthbertson J, Ashcroft FM, Ezaki T, Doyle DA. Crystal structure of the potassium channel KirBac1.1 in the closed state. Science 300: 1922–1926, 2003.
29. Labro AJ, Raes AL, Bellens I, Ottschytsch N, Snyders DJ. Gating of Shaker-type channels requires the flexibility of S6 caused by prolines. J Biol Chem 278: 50724–50731, 2003.
30. Lindahl E, Hess B, van der Spoel O. 3 DGROMACS: a package for molecular simulation and trajectory analysis. J Mol Mod 7: 306–317, 2001.
31. Long SB, Campbell EB, MacKinnon R. Crystal structure of a mammalian voltage-dependent Shaker family K+ channel. Science 309: 897–903, 2005.
32. Long SB, Campbell EB, MacKinnon R. Voltage sensor of Kv1.2: structural basis of electromechanical coupling. Science 309: 903–908, 2005.
33. Luzhkov VB, Nilsson J, Arhem P, Aqvist J. Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine. Biochim Biophys Acta 1652: 35–51, 2003.[Medline]
34. Magidovich E, Yifrach O. Conserved gating hinge in ligand- and voltage-dependent K+ channels. Biochemistry 43: 13242–13247, 2004.[CrossRef][Web of Science][Medline]
35. Marrink SJ, Berger O, Tieleman P, Jahnig F. Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations. Biophys J 74: 931–943, 1998.[Web of Science][Medline]
36. Morais-Cabral JH, Zhou Y, MacKinnon R. Energetic optimization of ion conduction rate by the K+ selectivity filter. Nature 414: 37–42, 2001.[CrossRef][Medline]
37. Olcese R, Sigg D, Latorre R, Bezanilla F, Stefani E. A conducting state with properties of a slow inactivated state in a Shaker k(+) channel mutant. J Gen Physiol 117: 149–164, 2001.
38. Rich TC, Yeola SW, Tamkun MM, Snyders DJ. Mutations throughout the S6 region of the hKv1.5 channel alter the stability of the activation gate. Am J Physiol Cell Physiol 282: C161–C171, 2002.
39. Rosenhouse-Dantsker A, Logothetis DE. New roles for a key glycine and its neighboring residue in potassium channel gating. Biophys J 91: 2860–2873, 2006.[CrossRef][Web of Science][Medline]
40. Sackin H, Nanazashvili M, Palmer LG, Li H. Role of Conserved Glycines in pH Gating of Kir1.1 (ROMK). Biophys J 90: 3582–3589, 2006.[CrossRef][Web of Science][Medline]
41. Sandoz G, Lopez-Gonzalez I, Grunwald D, Bichet D, Altafaj X, Weiss N, Ronjat M, Dupuis A, De Waard M. Cavbeta-subunit displacement is a key step to induce the reluctant state of P/Q calcium channels by direct G protein regulation. Proc Natl Acad Sci USA 101: 6267–6272, 2004.
42. Sanguinetti MC, Curran ME, Zou A, Shen J, Spector PS, Atkinson DL, Keating MT. Coassembly of K(v)LQT1 and minK (IsK) proteins to form cardiac IKs potassium channel. Nature 384: 80–83, 1996.[CrossRef][Medline]
43. Snyders DJ, Tamkun MM, Bennett PB. A rapidly activating and slowly inactivating potassium channel cloned from human heart. Functional analysis after stable mammalian cell culture expression. J Gen Physiol 101: 513–543, 1993.
44. Tamkun MM, Knoth KM, Walbridge JA, Kroemer H, Roden DM, Glover DM. Molecular cloning and characterization of two voltage-gated K+ channel cDNAs from human ventricle. FASEB J 5: 331–337, 1991.[Abstract]
45. Tieleman DP, Shrivastava IH, Ulmschneider MR, Sansom MS. Proline-induced hinges in transmembrane helices: possible roles in ion channel gating. Proteins 44: 63–72, 2001.[CrossRef][Web of Science][Medline]
46. Treptow W, Maigret B, Chipot C, Tarek M. Coupled motions between pore and voltage-sensor domains: a model for Shaker B, a voltage-gated potassium channel. Biophys J 87: 2365–2379, 2004.[CrossRef][Web of Science][Medline]
47. van Gunsteren WF, Berendsen HJC. Gromos-87 Manual. Groningen, The Netherlands: Biomos, 1987.
48. Webster SM, del Camino D, Dekker JP, Yellen G. Intracellular gate opening in Shaker K+ channels defined by high-affinity metal bridges. Nature 428: 864–868, 2004.[CrossRef][Medline]
49. Yifrach O, MacKinnon R. Energetics of pore opening in a voltage-gated k(+) channel. Cell 111: 231–239, 2002.[CrossRef][Web of Science][Medline]
50. Zagotta WN, Hoshi T, Aldrich RW. Shaker potassium channel gating III: evaluation of kinetic models for activation. J Gen Physiol 103: 321–362, 1994.
51. Zhou Y, Morais-Cabral JH, Kaufman A, MacKinnon R. Chemistry of ion coordination and hydration revealed by a K+ channel-Fab complex at 2.0 Å resolution. Nature 414: 43–48, 2001.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |