|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Laboratory of Chemical Physics, NIDDK, National Institutes of Health, Bethesda, Maryland 20892, USA
2 Department of Biochemistry, Howard Hughes Medical Institute, Brandeis University, Waltham, Massachusetts 02454, USA
(RECEIVED November 6, 2005; FINAL REVISION January 5, 2005; ACCEPTED January 5, 2005)
| Abstract |
|---|
|
|
|---|
Keywords: backbone assignment; detergent; global fold; membrane protein; NMR; NOE; potassium channel; secondary chemical shifts
| Introduction |
|---|
|
|
|---|
Several structures of potassium channels have been elucidated in recent years by X-ray crystallography (Doyle et al. 1998; Jiang et al. 2003; Kuo et al. 2003) and electron microscopy (Jiang et al. 2004), revealing the molecular basis of channel function and regulation. The first of these to be solved, the channel from Streptomyces lividans (KcsA), is an assembly of four identical membrane subunits characterized by an "inverted teepee" structure exhibiting fourfold symmetry. Each subunit consists of two membrane-spanning helices (TM1 and TM2) connected by a pore domain that includes the pore helix and the well-conserved TVGYG selectivity motif (Heginbotham et al. 1994). The main-chain carbonyl oxygen atoms of this segment form a K+-selective "coordination cage" capable of stabilizing K+, but not Na+, as these lose their hydration shell upon entering the pore. The aqueous pore narrows toward the cytoplasmic side of the membrane, and is lined by the selectivity filter and the TM2 helix of each subunit (Doyle et al. 1998).
The rapidly increasing structural information for K+ channels has been invaluable in understanding the molecular mechanisms underlying channel selectivity and high throughput (Gulbis and Doyle 2004). These crystal structures provide high-resolution "snapshots" of the channel in multiple states, but important questions remain regarding the precise dynamic processes underlying the gating mechanism. To date, most studies have focused on the transmembrane and pore regions, whereas the N- and C-terminal domains, which typically exhibit insufficient electron density, appear unstructured or were excised for obtaining high-quality crystals. These cytoplasmic domains have been implicated in mediating the formation of channel tetramers and pH- dependent gating (Cortes et al. 2001; Molina et al. 2004). For these reasons, structural studies of full-length KcsA that address dynamic issues at the molecular level should enhance our understanding of potassium channel function.
Nuclear magnetic resonance (NMR) has evolved in recent years into an important tool for studying the structure and dynamics of proteins in solution (Ferentz and Wagner 2000; Wider 2005). Several key methodological advances have allowed NMR to address soluble proteins of increasing size and complexity (Bax and Grzesiek 1993; Pervushin et al. 1997; Salzmann et al. 1998; Riek et al. 2002; Tugarinov et al. 2002; Frueh et al. 2005). As membrane-associated proteins play pivotal roles in a wide range of biological processes, application of NMR methods to the study of such systems has been a prime objective. Solid-state NMR methods have been used to study membrane peptides and small proteins, particularly transmembrane helical domains (Fu and Cross 1999; Marassi and Opella 1999; Mesleh et al. 2003). High-resolution NMR studies have afforded several structures of proteins in the 1530 kDa range (Arora et al. 2001; Hwang et al. 2002; Fernández and Wüthrich 2003; Oxenoid and Chou 2005). However, large (>30 kDa) membrane-associated proteins continue to present a formidable challenge in NMR studies. These highly hydrophobic proteins must be sequestered from the aqueous environment in relatively fluid macromolecular lipid or detergent assemblies. As a result, samples are prone both to line broadening due to exchange processes as well as aggregation. Furthermore, such proteindetergent systems rapidly reach the limit at which long tumbling times preclude the acquisition of meaningful NMR data (Arora and Tamm 2001). This is particularly true for helical proteins, for which spectral dispersion is typically poor. Clearly, NMR methods suitable for the study of large membrane proteins must be established if it is to remain relevant for structural investigation of membrane- associated systems.
These considerations motivated us to examine the KcsA tetramer as a model system for developing NMR methods that bring relatively large helical transmembrane systems within the reach of solution NMR. Considering that it can be expressed and purified in sufficient quantities, is suitable for isotopic labeling, and possesses high thermostability, KcsA represents an attractive system for this purpose. In addition, NMR is well suited for probing dynamic processes and also provides access to the terminal regions of the KcsA chain, which may play an important role in tetramerization and regulation of ion conductance. By surveying samples of KcsA solubilized in various detergent micelles, we established that sodium dodecyl sulfate (SDS) micelles containing the homotetrameric 68-kDa channel were amenable to backbone assignment. The secondary structure of KcsA determined by secondary chemical shifts and hydrogen-exchange protection factors closely agrees with the KcsA crystal structure, and indicates the presence of a relatively long a-helix in the C-terminal cytoplasmic domain. Quantitative cross-relaxation attenuation experiments have been used to map interactions between detergent and protein, and results indicate that the detergent forms a hemitoroid around the membrane-spanning hydrophobic domains of the channel. The KcsA channel in SDS micelles exhibits K+-induced chemical shift changes in selectivity filter, indicating that the channel maintains its biological activity under these conditions. Correlated pH-induced changes observed for both TM domains suggest a possible role for an intersubunit interaction in the pH-gating of the channel.
| Results |
|---|
|
|
|---|
As has been reported previously for KcsA in various detergent micelles (Heginbotham et al. 1997; Molina et al. 2004), KcsA in SDS micelles behaves as a highly stable molecular ensemble. It migrates as a tetramer in conventional SDS-PAGE even after exposure to elevated temperatures (75°C for 60 min), pH 10 or 4 M urea. Dynamic light scattering (DLS) measurements indicate a molecular weight of 115 ± 10 kDa for the proteindetergent ensemble, corresponding to an aggregation number of 160 ± 35 SDS molecules per KcsA tetramer (Fig. 1A). Circular dichroism (CD) measurements of KcsA in SDS indicate that it is ~50% helical, and highly similar in helical content to KcsA in dodecylphosphocholine (DPC) and DM micelles (Fig. 1B). DPC has been acknowledged as a suitable membrane surrogate in structural studies (Oxenoid et al. 2004), and DM is a mild and often-used nonionic detergent for stabilization of membrane proteins. This evidence, together with specific chemical shift changes in the selectivity filter upon addition of K+ (vide infra), strongly suggests that KcsA in SDS micelles embodies a valid representation of the cellular conformation of the channel.
General approach to backbone assignment of KcsA
As a large proteindetergent assembly, exceeding 100 kDa, backbone assignment of KcsA requires deuteration of all 13C nuclei, while maintaining full protonation in backbone amide sites, without which triple resonance experiments cannot be efficiently acquired (Grzesiek et al. 1993; Tugarinov et al. 2004). This is typically achieved by performing expression in D2O-based medium, followed by base-catalyzed back-exchange of amide protons during purification. In cases where this is insufficient, denaturation and refolding of the protein may be used to obtain full protonation of backbone amide sites (Tugarinov et al. 2004). However, the highly stable KcsA resisted all attempts at back-exchange of amide protons at its detergent-embedded core. As a result, KcsA expressed in D2O-based medium exhibited only 68% of the expected amide proton cross-peaks. Extended exposure to pH values as high as 10 and temperatures as high as 75°C, just under its melting point, failed to introduce protons into these well-protected sites. While successful refolding of denatured KcsA into its tetrameric form has been reported (Valiyaveetil et al. 2002), in our hands, the overall yield of this process remains insufficient for preparation of NMR samples.
We resolved this difficulty by treating KcsA as a two-domain protein, extending the approach posited by Lohr and coworkers (Lohr et al. 2003). At the basis of this approach lies the assumption that nonexchangeable protons will be heavily clustered along the sequence of KcsA. Accordingly, backbone assignment was performed concurrently for two samples. The first contained 2H,13C,15N-labeled KcsA expressed in D2O-based M9 medium, allowing the assignment of the solvent-accessible (exchangeable) domain of KcsA (KcsAE). The second contained KcsA expressed in glucose-deficient H2O-based medium supplemented with 2H,13C,15N-labeled nutrients. This afforded a 2H(1H
),U-13C,15N-labeled KcsA sample, for which all amide sites were fully protonated, side-chain 13C sites were highly deuterated, and 13C
sites were fractionally deuterated at a level that varied with amino acid type. This sample was lyophilized and reconstituted in D2O, selectively quenching signals arising from the solvent-accessible domain, and allowing the assignment of the transmembrane, or core domain of KcsA (KcsATM). The 1H,15N-TROSY-HSQC (tr-HSQC) "fingerprint" regions for the two domains are shown in Figure 2.
Although this approach requires the acquisition of NMR data for both samples, doubling the overall measurement time, these two complementary samples facilitated the assignment process by dividing the data into two self-contained subsets.
Backbone assignment of KcsA
The differences in resonance line width between the KcsAE and KcsATM domains (Fig. 2) are striking, and indicate that internal mobility is strongly decreased for protein segments enclosed in the micelle. An average tC value of 40 nsec is estimated for the KcsATM domain on the basis of 15N relaxation measurements, whereas KcsAE transverse relaxation rates were quite heterogeneous and, on average, about twofold lower, indicative of extensive internal mobility for the KcsAE domain. These differences in mobility determined the appropriate strategy for backbone assignment of the two domains. In the case of KcsAE, line widths were sufficiently narrow to obtain intra- and inter-residual connectivities for both 13C' and 13C
nuclei. Thus, the assignment of its backbone resonances was predominantly based upon correlations of 13C' and 13C
nuclei to their adjacent amide protons. In contrast, the broader line widths of KcsATM rendered the out-and-back three-step transfer experiments tr-HN(CO)CA and tr-HN(CA)CO (Salzmann et al. 1998) inefficient, precluding detection of the inter- residue connectivities. The intrinsically helical character of this domain and its favorable dispersion of peaks in the tr-HSQC spectrum enabled us to complete the assignment utilizing a set of 15N-separated NOESY experiments. This strategy is demonstrated in Figure 3. By independent optimization of the assignment strategy for each domain, an overall assignment level in excess of 97% of native non-proline residues was obtained.
Figure 4 describes the distribution of KcsA residues between these two domains. The KcsATM domain is characterized by amide protons that remain nonexchangeable, or very slowly exchangeable, over a period of several weeks in a 99% D2O sample at pH 6 and 323 K. Protein segments included in KcsATM were residues 3550, 6579, and 87111. The first and third segments roughly correspond to the two transmembrane helices (TM1 and TM2, respectively) in the KcsA structure, and the second overlaps with the pore helix and selectivity filter. All other residues were accessible to solvent, with backbone amide solvent exchange rates in the 0.00110 sec1 range, and comprise the KcsAE domain. Several residues could be observed in both samples, including L35, L49-A50, G79, W87-G88, and T107-A111, corresponding to residues located at the boundaries of the two domains.
The secondary structure of KcsA
Using the well-known correlation between backbone13 C chemical shifts and secondary structure, the local backbone structure of KcsA in SDS micelles could be established (Fig. 5). NMR data define four helical segments in KcsA, connected by loops with extended backbone conformation. Strong helical indicators are found for residues 3152, 6273, and 86115, corresponding to the positions of TM1, the pore helix, and TM2, respectively. A fourth helical segment, spanning residues 142155, represents an additional helix in the C-terminal cytoplasmic domain. The 13C chemical shift helical indicators in this helix are structurally and statistically significant, although not as strong as observed for the TM domains. The four KcsA helices are similarly identified and located by an analysis of NOE cross-peaks between adjacent backbone amide protons. The vast majority of residues in the proposed helical sequences exhibit a strong 1HN(i)-1HN(i + 1) interaction, which is characteristic in segments with helical structure. The distinction between structured and disordered segments in the C-terminal domain (residues 121160) is particularly striking. The total helical content of KcsA as determined by NMR (78 of 152 residues) is consistent with the results of CD measurements (Fig. 1B).
Rates of exchange of amide protons with the solvent (HX rates) have been extensively used as a measure of secondary structure in proteins (Wand and Englander 1996; Maity et al. 2003). These were quantitatively derived for the KcsAE domain by comparing the amide proton intensities observed with and without selective inversion of the H2O resonance prior to acquisition of a 3D-TROSY-HNCO (tr-HNCO) spectrum. At pH 8.0, residues 121135 and 142155 have average exchange rates of 10.5 (maximum of 36) and 1.5 (maximum of 5) sec1, suggesting the presence of a helix in the latter segment. Similarly, the four N-terminal residues of KcsA (1619) and the next six residues (2025) average an exchange rate of 2.4 (maximum of 3.9) and 19.5 (maximum of 36) sec1, respectively. In these segments the large increase in hydrogen exchange rate of G21 over that of G18 is particularly striking (Fig. 5). Consistent with secondary shift data, this suggests the existence of an additional helix in the N-terminal domain, the majority of which is excluded from the construct used for this work.
Proteindetergent interactions
Having analyzed the secondary chemical shifts in KcsA, we turn our focus to the assembly of detergent molecules, which surround it and stabilize it in aqueous solution. By measuring dipoledipole cross-relaxation between detergent and protein backbone amide protons, NMR can characterize the spatial organization of micelle-solubilized KcsA. Considering that identification of individual backbone amides required three-dimensional dispersion of the signals in a tr-HNCO spectrum, we opted for a set of selective SDS inversion experiments to probe the NOE interaction between SDS and KcsA, rather than extending the three-dimensional experiment to four dimensions, or resorting to reduced-dimensionality experiments. The selective NOE measurements are carried out by inverting individual SDS resonances prior to acquisition of the 3D-tr-HNCO data set. Inversion pulses afforded a 70%80% inversion of the selected detergent proton, whereas its effect upon nonexcited SDS protons was minimal. The NOE-related attenuation during a mixing time of 250 msec was measured as the ratio between the intensities associated with a given tr-HNCO signal in spectra recorded in interleaved fashion with and without the inversion pulse.
The 1H SDS spectrum exhibits four signals, corresponding to the C1H2, C2H2, C311H2, and C12H3 protons, resonating at 3.8, 1.5, 1.25 (broad), and 0.85 ppm, respectively. Results from the experiment involving inversion of the downfield protons are difficult to interpret due to overlap with non-negligible levels of C
protonation in KcsATM, thereby giving rise to strong intraprotein NOEs. The results for selective inversion of each of the three upfield SDS protons are summarized in Figure 6. Overall, only moderate attenuation of the HN magnetization, even after 250 msec, is observed. As expected, the proteindetergent interaction is strongest for membrane-buried domains, exhibiting maximum NOE attenuations of 15%, 30%, and 15%, for inversion of the C2H2, C311H2, and C12H3 protons, respectively. The proteindetergent interactions observed for the "outer" and "inner" transmembrane helices (TM1 and TM2, respectively) are comparable. In contrast, the entire C-terminal domain (residues 121160) is hardly affected by inversion of detergent protons, indicating that few, if any, micellar interactions exist for this domain. Notably, NOE attenuation increases at the N-terminal residues 1619, particularly when compared with the C-terminal segment. Interactions between these residues and the detergent hint once again at the existence of a membrane- embedded N-terminal helix in full-length KcsA.
A comparison between attenuation maps for different types of SDS protons is revealing. Generally speaking, highest attenuation levels are observed for the C311H2 protons, making them the most sensitive probe of proteindetergent interactions. However, since the overlapping C311H2 protons span the length of the SDS monomer, the results for the C12H3 and C2H2 protons are used to interpret the relative orientation (in the head-to-tail sense) of the detergent in the proteinmicelle assembly. NOE attenuations for the C311H2 and C12H3 protons are highly correlated for the TM1 and TM2 helices, whereas for the pore helix it is the attenuations of the C311H2 and C2H2 protons that are correlated. These results are summarized in Figure 7, displaying the results of SDS-inversion experiments juxtaposed upon the structure of KcsA. The attenuation maps for C2H2 (Fig. 7A) and C12H3 (Fig. 7C) protons are highly complementary. Whereas inversion of the SDS-charged headgroups affects KcsA residues at the membranesolvent interface, inversion of the hydrophobic tails affects residues in the membrane-spanning domain. The selectivity filter appears relatively sheltered from detergent NOEs, as expected for this highly buried sequence. With the exception of N-terminal residues 1619, KcsA segments in the extracellular and cytoplasmic domains are characterized by low attenuation levels in all experiments, thereby clearly demarcating the membrane interface.
K+-induced chemical shift changes in the KcsA selectivity filter
The KcsA selectivity filter uses a coordinated array of backbone carbonyl oxygens and side-chain hydroxyl groups of residues 75TVGY79G to selectively stabilize incoming K+ ions as they traverse the membrane (Doyle et al. 1998). Titration of tetrameric KcsA in SDS micelles with KCl affects the chemical shifts of residues 73ATTV77G, a segment that forms the two "inner" coordination sites (of four) for K+. The largest difference is observed for the backbone amide of V76, whose 15N chemical shift changes by about 1.4 ppm (Fig. 8A). Additional residues affected by the presence of K+ include L41 and A42 of the TM1 helix and A98 of TM2. Overall, the results indicate an equilibrium between K+-bound and K+-deficient states, which, for conditions at which NMR data were acquired, is governed by a dissociation constant of KD ~3 mM (Fig. 8B). For V76, an additional 13 Hz 15N line broadening is observed at the midpoint of the titration relative to the endpoints (Fig. 8A), which, using standard equations for two-site exchange, indicates an exchange rate kex
104 s1.
Recently solved channel structures have established that the filter is electrostatically driven to adopt an alternative conformation in the absence of K+, and that residue V76 is most affected by this change (Zhou et al. 2001; Zhou and MacKinnon 2003). Furthermore, comparable KD values for K+ high-affinity sites have been reported in previous physiological (LeMasurier et al. 2001) and structural (Zhou and MacKinnon 2003) titration studies. Therefore, our NMR results strongly suggest that the functionality of the channel is preserved in SDS micelles.
pH-induced changes in KcsA
KcsA exhibits pH-dependent gating, shifting to the open state in bilayers when its cytoplasmic face is exposed to pH values below ~5 (Cuello et al. 1998; Heginbotham et al. 1999). It is therefore interesting to follow pH-induced changes in KcsA chemical shifts in an attempt to account for this gating mechanism on the molecular level. The KcsA backbone chemical shifts were reassigned for samples in 20 mM Tris (pH 8.0) buffer and compared with values observed in 25 mM MES (pH 6.0). Weighted chemical shift deviations (
=sqrt((
H)2+(
N/5)2)) in the tr-HSQC spectra are summarized in Figure 9. Most changes in chemical shift found in the C-terminal cytoplasmic domain can be attributed to local electrostatic effects of charged residues affected by the pH change. An important exception in this regard is the change observed for residues 114117 at the cytoplasmic end of the TM2 helix. Equally interesting are the smaller, yet significant chemical shift changes in the KcsATM domain, despite its being sheltered from the macroscopic environment by the SDS micelle. Notably, these changes cluster in the TM1 helix; of the eight residues exhibiting the largest chemical shift changes (>0.07 ppm), and which are unaffected by local electrostatic effects, six are located in this domain. Residue T75 of the selectivity filter shows considerable pH-induced chemical shift change as well. In both KcsAE and KcsATM, the changes in backbone 13C chemical shifts are relatively small, exhibiting RMS values of 0.2 and 0.1 ppm, respectively. The small nature of the chemical shift changes indicates that the KcsA time-averaged structure undergoes minimal change in the range of pH values studied.
| Discussion |
|---|
|
|
|---|
Comparison with KcsA structures
A comparison of our secondary chemical shift, amide proton NOE connectivities, and hydrogen-exchange NMR data with the crystal structure reported for residues 23119 of KcsA (Doyle et al. 1998) reveals excellent agreement. The three previously identified structural elements of KcsA, the TM1, pore, and TM2 helices are located in the crystal structure at residues 2751, 6274, and 86112. NMR C
secondary chemical shifts identify these structural elements at residues 3152, 6273, and 86115, respectively. The R27AAGA31 segment preceding the TM1 helix shows slow HX rates and strong sequential amide NOEs, suggestive of a helical conformation for this region too. Ala and Gly residues within the transmembrane helices tend to have smaller 13C secondary chemical shifts than other residues, possibly accounting for the modest secondary chemical shifts seen in this region. The unusually elevated apparent HX rate of residue T112 is likely to be an artifact resulting from OH-mediated intraresidue NOE attenuation.
The cytoplasmic N- and C-terminal segments of KcsA represent obvious putative sites for binding cellular factors that regulate channel opening. Therefore, structural elements in these regions are of great interest, particularly since previously they were removed in order to obtain sufficiently diffracting crystals. In a site-directed spin-labeling (SDSL) study of full-length KcsA (residues 1160) in liposomes (Cortes et al. 2001), N-terminal residues (120) were found to form a helix that enters the membrane at residues 1516 and is linked by a disordered loop of 56 residues to the TM1 domain. In our study, both secondary shifts and HX rates suggest that residues 1619 form a helical turn, whereas the following residues are disordered. We conclude, therefore, that the helical character of native KcsA residues is preserved in our shorter construct. In contrast, SDSL and NMR data provide conflicting views regarding secondary structure in the C-terminal domain. The spin-labeling study predicts the C-terminal region to contain two helical stretches, spanning residues 129148 and 154160, whereas our NMR data locate a single C-terminal helix at residues 142155, and residues 124137 and 156160 appear to lack any significant structure. Notably, the SDSL study was performed at low pH (~4) and in liposomes, rather than a strongly anionic detergent. A recent study established structural differences between KcsA solubilized in native-like membranes and detergents (Encinar et al. 2005). Nevertheless, it is noteworthy that the sequence of the helical segment proposed in our NMR results is amphiphilic, making it conducive to the formation of a tetrameric helical bundle. In this assembly, residues L144, F148, L151, M154, and L155 would create the hydrophobic surface of the helix facing inward to the fourfold symmetry axis. R147 is the only residue opposing such an arrangement, but it is conceivable that with moderate distortions in the helices, its relatively long side chain allows the guanidino group to reach the solvent. The opposite face of the helix is rich in charged Asp and Arg residues, which could stabilize this region in the aqueous environment. Although we cannot unambiguously establish the presence of a tetrameric arrangement for the C-terminal helix from the data at hand, we would expect lower transverse 15N relaxation rates than observed for this helix and the unstructured region preceding the helix if it were monomeric. Therefore, our results appear more compatible with at least partial occupation of a parallel four-helical bundle arrangement.
Residue-specific modulation of secondary chemical shifts in the membrane environment
With the increasing available volume of backbone 13C chemical shift data for proteins, secondary chemical shifts are a primary source of structural information used in numerous structure determination algorithms (Kuszewski et al. 1995; Cornilescu et al. 1999; Wishart 2001). However, data on the behavior of chemical shifts in membrane proteins lags behind the abundant information available for their soluble counterparts. The backbone assignment of KcsA helps to address this paucity of data. Although based upon limited data (a total of 60 residues), a characteristic pattern of residue- specific secondary shift magnitudes appears evident. This is best appreciated by comparing 13C
chemical shifts in the three buried KcsA helices to average helical values (as opposed to random coil values) (Zhang et al. 2003). Analysis shows that the largest deviations are observed for the
-OH-substituted residues Ser (3) and Thr (4), averaging downfield shifts of 2.4 and 1.2 ppm, respectively, and the
-branched residues Ile (2) and Val (12), averaging downfield shifts of 1.6 and 0.9 ppm (the number of averaged residues appears in parentheses). Similar effects are observed when analyzing the chemical shifts of the transmembrane helical residues in diacylglycerol kinase (DAGK) (Oxenoid et al. 2004).
Notably, the four amino acids whose chemical shifts are most affected in membrane helices also experience the largest downfield shift when moving from random coil to helical conformation (Zhang et al. 2003). This suggests that the lipophilic environment is amplifying effects exerted upon helices in soluble proteins as well. Formation of hydrogen bonds between
-OH groups and the carbonyl in the preceding turn of the helix is favored in both aqueous (Gray and Matthews 1984) and membrane-spanning (Chamberlain and Bowie 2004)
-helices. We propose that as this interaction becomes more dominant in the lipophilic environment, the resulting equilibrium shift accentuates the "helicity effect" and accounts for the larger downfield deviation observed. This consideration also correctly predicts a more pronounced effect upon Ser compared with Thr residues. Since the Ser side chain, lacking the
-substituent, would have more conformational freedom prior to the formation of the
-OH hydrogen bond, the relative effect of the hydrophobic environment upon its chemical shift would be greater.
Cross-relaxation effects and structural characterization of SDSKcsA mixed micelles
Detergents and lipids are an inseparable part of membrane protein structure. Previous structural studies have focused upon this key participant in the membrane protein assembly by monitoring the appearance of NOESY cross-peaks between various detergent and protein protons (Fernández et al. 2002; Rooslid et al. 2005), the relaxation-enhancing effects of spin-labeled detergent molecules (Hilty et al. 2004), or estimating detergentprotein cross-relaxation (Ulmer and Bax 2005). Here, we study the assembly of SDS molecules around the KcsA tetramer by following the attenuation of amide proton signals resulting from an inversion of specific detergent proton populations. The slow NOE buildup between SDS and KcsA implies that while effective correlation times
int for the SDSKcsA interaction are longer than the 
int = 1 limit (
int >200 psec), they are much shorter than the overall tumbling time of the KcsA-containing micelle. Notably, the NMR spectrum shows narrow SDS resonances, with relaxation properties not measurably affected by the KcsA, indicative of rapid exchange between monomeric SDS, free SDS micelles, and KcsA-containing micelles. The detergent must therefore be treated as a highly fluid medium, and NOE attenuations are a function of both the occupancy and residence time of detergent at the protein interface. In contrast, intraprotein NOEs between amide protons are governed by the long protein
C and build up far more rapidly, particularly in helical domains, and result in spin diffusion. The slow SDSKcsA NOE buildup makes the use of short mixing times (
20 msec) that would avoid such spin diffusion impractical. Intraprotein NOEs account for the comparable attenuation observed for the two TM domains. Although the "inner" helix (TM2) is less exposed to the environment, one in every three to four residues does interact with detergent protons, and the attenuation then rapidly propagates by spin diffusion in the slow tumbling KcsATM domain. Since the relation
int <<
C generally holds for large membrane proteins, proteindetergent NOES may be better suited for mapping the exposure of protein segments, rather than particular residues, to detergent.
Despite these limitations, a clear structural picture emerges from our results. The data in Figure 7 can be explained by a population of SDS molecules oriented parallel (or antiparallel) to the TM helices, "hugging" the surface of the membrane-spanning domains. This would account for the observed distribution of C2H2 -andC12H3 -induced attenuations along the TM domains. In fast exchange with this population must be additional SDS molecules oriented perpendicularly to the channel fourfold axis, with their hydrophobic tails contacting KcsA, since two SDS molecules fail to span the entire membrane (Fig. 10). In this hemitoroidal arrangement, the micelle optimally covers the surface of KcsA, while the interface between the hydrophobic detergent molecules and the aqueous environment is minimized. To relate this macromolecular assembly with the physical parameters of KcsA, we model the channel as an open-ended cone, with a height of 34 Å and base diameters of 47 Å and 25 Å for the extra-cellular and cytoplasmic sides, respectively. These dimensions include an additional 2 Å on either side of the backbone to account for the spatial requirements of the hydrophobic side chains. Assuming the SDS hydrocarbon chains occupy a cylinder 16 Å long and with a 4.3 Å diameter, we calculate that ~190 SDS molecules are required to fulfill the hydrophobic requirements of KcsA, and estimate the total micellar molecular weight at ~122 kDa. This number agrees well with the results of DLS measurements, which predict a molecular weight of 115 ± 10 kDa, as well as the KcsATM correlation time of 40 nsec at 323 K, which corresponds to a molecular weight of ~120 kDa.
Conformational changes associated with gating in KcsA
The pH dependence of ion conductance through the potassium channel has been a heatedly debated subject in recent years (Cuello et al. 1998; Heginbotham et al. 1999; Zakharian and Reusch 2004). The biological significance of a pH-gating mechanism is not clear, in view of the fact that intracellular pH must be closely controlled (Corvini et al. 2000). It is assumed that the pH change mimics the action of a cellular effector, which binds to the cytoplasmic domain of the channel. Since KcsATM is devoid of residues with pH-dependent electrostatic characteristics, and presumably insulated by the hydrophobic interior of the micelle from the macroscopic pH, changes in chemical shift clustered in the TM1 helix may reflect structural differences between high- and low-pH states of the channel. Notably, somebut not allof the changes observed when moving from pH 8 to pH 6 are echoed in the KCl-titration. This suggests the K+- and pH-induced changes in KcsATM and the selectivity filter may be structurally related.
Because of their unique side-chain pKa, histidine residues are often involved in pH-dependent conformational changes. The crystal structure of KcsA23119 (Doyle et al. 1998) suggests a single possible long-range histidine interaction, residue H25, which precedes the TM1 domain packing against residues 114118 at the base of the TM2 helix of the adjacent subunit. This segment has been postulated to represent a binding site for cellular modulators of channel activity (Cortes et al. 2001). Indeed, chemical shift changes in this segment represent the only effect observed in the C-terminal region unaccounted for by local electrostatic changes. The NMR data suggest that conformational changes at the convergence of the four subunits could be propagated to the TM1 helix and influence the flux of K+ ions through the channel by a mechanism yet to be uncovered. Although chemical shift changes in TM1 are not consistent with a large-scale conformational change in KcsA, it should be noted that we have not monitored these changes at lower pH values, since SDS-solubilized tetrameric KcsA loses its thermostability at such conditions. Chemical shift changes at pH 6 may therefore reflect only a small alteration in the equilibrium between two channel states, merely hinting at the magnitude of the effect expected when the channel moves to a fully open state.
| Materials and Methods |
|---|
|
|
|---|
Expression of KcsA followed general guidelines previously proposed (Gardner and Kay 1998). OverExpress C41(DE3) Escherichia coli cells (Avidis) (Miroux and Walker 1996) containing the KcsA plasmid were cultured in a modified M9 medium (Cai et al. 1998) at 37°C until an OD600 of 0.6 was reached. Cells were then transferred into a D2O-based M9 medium containing 4 g/L 13C-glucose and supplemented with 1 g/L 2H,13C,15N-labeled Isotec medium (Sigma-Aldrich) and grown to OD600 of 0.8. Isopropylthiogalactoside (IPTG) was added to a final concentration of 1 mM and shaking was continued for 16 h at 27°C. For expression of KcsA with amide protonation in the transmembrane domains (KcsATM), the D2O-based medium was replaced by a glucose-deficient H2O-based M9 medium containing 47 g/L 2H,13C,15N-labeled Isogro medium (Isotec).
Purification of tetrameric KcsA followed the general protocol previously outlined by Shimada and coworkers (Takeuchi et al. 2003) with minor modifications. After harvesting, cells were lysed by sonication in buffer containing 20 mM Tris (pH8), 150 mM KCl, 150 mM NaCl, 20 mM imidazole, and 5 mM benzamidine. KcsA was extracted for 6090 min at 4°C using 20 mM DM, and the suspension was clarified by centrifugation. The supernatant was loaded on a Ni++-charged HiTrap Chelate column (Amersham Biosciences) and eluted using 400 mM imidazole in Tris (pH 8) buffer containing 5 mM DM. KcsA was concentrated to 11.5 mL in a Centriprep YM10 device (Millipore) and fractionated on a size-exclusion column (Superdex 200 16/60, Amersham Biosciences) equilibrated with 25 mM N-morpholino-ethane-sulfonate (MES) buffer at pH
6.0 (20 mM Tris at pH 8.0 for samples at pH 8) and 10 mM sodium dodecyl sulfate (SDS). For preparation of samples in other detergents, KcsA in SDS buffer was pooled and concentrated before reloading again on the same size-exclusion column, equilibrated with appropriate buffer and appropriate detergent. KcsA in SDS micelles was concentrated using Centriprep and Amicon Ultra devices (Millipore) to volumes of 300700 µL for NMR measurements.
Sample preparation
NMR measurements were recorded at concentrations of 0.250.35 mM tetrameric KcsA (1.01.4 mM monomeric concentration) in 25 mM MES (pH 6.0) or 20 mM Tris (pH 8.0) and a 100200:1 excess of detergent over monomer. Samples with amide protonation in the solvent-accessible domains (KcsAE) were prepared in 93% H2O, 7% D2O. Samples with amide protonation in the transmembrane domains (KcsATM) were lyophilized once and dissolved in 99.9% D2O.
NMR measurements
All NMR measurements were conducted at 323 K. Sample stability at elevated temperatures was verified using SDSpolyacyrlamide gels and by NMR. Data were collected on DRX600 and DRX800 Bruker spectrometers using cryogenic triple-resonance probeheads equipped with z-axis pulsed field gradients. Triple resonance, TROSY versions (Salzmann et al. 1998) of the HNCO, HN(CO)CA, HN(CA)CO, and HNCA experiments for purposes of backbone assignment were carried out for 2H,13C,15N-labeled KcsAE at 600 MHz field using a TROSY (tr) block similar to that suggested by Pervushin and coworkers (Pervushin et al. 1998) optimized for water suppression on cryoprobes (A. Bax, unpubl.). Triple-resonance experiments were typically acquired with 3236 complex points and an acquisition time of 20.222.8 msec in the 15N dimension, and 512 complex points and an acquisition time of 63.2 msec in the observed dimension. Experiments with 13CO (13Ca) evolution were acquired with 3240 complex points and 21.226.5 (77.5) msec acquisition time. 3D 1H-15N NOESY-HMQC and 15N HMQC-NOESY-HMQC spectra were acquired at 800 MHz field with mixing times of 7080 (250) msec for KcsATM (KcsAE) samples, respectively. These were acquired with 4050 complex points and an acquisition time of 1923.7 msec in the 15N {F2} dimension, and 512 complex points and an acquisition time of 45.9 msec in the observed dimension. Experiments with 1HN (15N) evolution were acquired with 72 (50) complex points and an acquisition time of 13.6 (23.6) msec. NOESY experiments utilized a selective E-BURP first 1H pulse to minimize excitation of solvent and detergent protons.
To overcome the poor dispersion of 1H-15N signals for KcsAE, proteindetergent NOE attenuations and hydrogen-exchange rates were measured using three-dimensional tr-HNCO-based experiments. In the case of proteindetergent dipolar interactions, SDS protons were selectively inverted using an IBURP2 sequence (Geen and Freeman 1991) and a 40-msec selective pulse, and allowed a 50250 msec mixing time preceding the tr-HNCO sequence. This was repeated for each of the four distinct proton chemical shifts in the 1D spectrum of SDS. The SDS-to-protein NOE effect increases steadily with mixing time, and the highest signal-to-noise data, collected at 250 msec, were used for analysis. For estimating hydrogen-exchange rates of SDS-KcsA in 20 mM Tris (pH 8.0), the H2O protons were selectively inverted prior to the tr-HNCO sequence, with radiation damping suppressed using a weak 1 G/cm gradient throughout a mixing time of 25 msec. In all experiments, a reference spectrum, in which the inversion pulse was replaced by an appropriate delay, was acquired in interleaved fashion, and care was taken to purge steady-state 15N magnetization.
K+-titration of KcsA was performed by successively adding aliquots of a 50 mM KCl solution in 100 mM MES (pH 6.0) in 99.9% D2O to a NMR sample of 0.18 mM tetrameric KcsA (determined by absorbance at 280 nm) in 200 mM SDS and 25 mM MES (pH 6.0) in 99% D2O. A 1H-15N-tr-HSQC (108 complex points and 50-msec acquisition time in the 15N dimension, total measurement time of 100 min.) was acquired at K+ concentrations of 09 mM on a DRX800 Bruker spectrometer at 323 K. Peak position for the amide 15N of residue V76 was used to estimate the K+-binding constant.
Data analysis
In all NMR experiments, FIDs were weighted by sine square and sine functions (both shifted 72° and truncated at 168°), in acquisition and evolution dimensions, respectively, and Fourier transformed. All spectra were processed and examined using the NMRPipe/NMRDraw software package (Delaglio et al. 1995). Backbone assignment was performed with the aid of in-house scripts. The dipolar interaction between protein and detergent protons was qualitatively assessed using the ratio of tr-HNCO peaks in the interleaved inversion and reference spectra. Solvent exchange rates for amide protons were estimated in a similar fashion, and this intensity ratio was used to calculate exchange rates assuming a mono-exponential recovery curve.
CD measurements
CD spectra of KcsA (at 25 µM monomeric concentration) were recorded in various buffer/detergent mixtures on a JASCO J-700 spectropolarimeter using a 0.05-cm path length cell at 25°C. Quantitative evaluation of secondary structure from the CD spectrum was carried out using the program CDNN (www.bioinformatik.bio-chemtech.uni-halle.de/cd_spect/index.html).
Dynamic light scattering measurements
Light scattering data were obtained using an analytical Super-dex-200 column 10/30 (Amersham Biosciences) with in-line multiangle light scattering (DAWN EOS, Wyatt Technology, Inc.) and refractive index detectors (OPTILAB DSP, Wyatt Technology, Inc.) and processed using the Astra V 5.1.3.0
[EC]
software. Samples of 150300 µg of protein in 125 µL of either 20 mM Tris-HCl (pH 8) or 25 mM sodium phosphate (pH 6) containing 10 mM SDS were applied to the pre-equilibrated S200 column running at a flow rate of 0.5 mL/min at room temperature and eluted with the same buffer.
|
|
|
|
|
|
|
|
|
|
| Footnotes |
|---|
Article published online ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.051954706.
Abbreviations: CD, circular dichroism; DLS, dynamic light scattering; DM, dodecyl maltopyranoside; DPC, dodecylphosphocholine; HMQC, heteronuclear multiple-quantum coherence; HX, hydrogen solvent- exchange rate; KcsAE, exchangeable domain of KcsA; KcsATM, trans- membrane domain of KcsA; MES, 2-[4-morpholino]-ethanesulfonic acid; NOE, nuclear Overhauser effect; NOESY, nuclear Overhauser effect spectroscopy; NMR, nuclear magnetic resonance; SDS, sodium dodecyl sulfate; SDSL, site-directed spin-labeling; TM1(2), helical trans- membrane domain 1(2); tr-HSQC, 1H-15N transverse relaxation-optimized heteronuclear single-quantum coherence.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Arora A., Abildgaard F., Bushweller J.H., Tamm L.K. 2001. Structure of outer membrane protein A transmembrane domain by NMR spectroscopy Nat. Struct. Biol. 8: 334338.[CrossRef][Medline]
Bax A. and Grzesiek S. 1993. Methodological advances in protein NMR Acc. Chem. Res. 26: 131138.[CrossRef]
Cai M., Huang Y., Sakaguchi K., Clore G.M., Gronenborn A.M., Craigie R. 1998. An efficient and cost-effective isotope labeling protocol for proteins expressed in Escherichia coli J. Biomol. NMR 11: 97102.[CrossRef][Medline]
Chamberlain A.K. and Bowie J.U. 2004. Analysis of side-chain rotamers in transmembrane proteins Biophys. J. 87: 34603469.[CrossRef][Medline]
Cornilescu G., Delaglio F., Bax A. 1999. Protein backbone angle constraints from searching a database for chemical shift and sequence homology J. Biomol. NMR 13: 289302.[CrossRef][Medline]
Cortes D.M., Cuello L.G., Perozo E. 2001. Molecular architecture of full-length KcsA: Role of cytoplasmic domains in ion permeation and activation gating J. Gen. Physiol. 117: 165180.
Corvini P.F., Gautier H., Rondags E., Vivier H., Goergen J.L., Germain P. 2000. Intracellular pH determination of pristinamycin-producing Streptomyces pristinaespiralis by image analysis Microbiol. 146: 26712678.
Cuello L.G., Romero J.G., Cortes D.M., Perozo E. 1998. pH-dependent gating in the Streptomyces lividansK+channel Biochemistry 37: 32293236.[CrossRef][Medline]
Delaglio F., Grzesiek S., Vuister G.W., Zhu G., Pfeifer J., Bax A. 1995. NMRPipe: A multidimensional spectral processing system based on UNIX pipes J. Biomol. NMR 6: 277293.[Medline]
Doyle D.A., Cabral J.M., Pfuetzner R.A., Kuo A., Gulbis J.M., Cohen S.L., Chait B.T., MacKinnon R. 1998. The structure of the potassium channel: Molecular basis of K+conduction and selectivity Science 280: 6977.
Encinar J.A., Molina M.L., Poveda J.A., Barrera F.N., Renart M.L., Fernandez A.M., Gonzalez-Ros J.M. 2005. The influence of a membrane environment on the structure and stability of a prokaryotic potassium channel KcsA FEBS Lett. 579: 51995204.[CrossRef][Medline]
Ferentz A.E. and Wagner G. 2000. NMR spectroscopy: A multifaceted approach to macromolecular structure Q. Rev. Biophys. 33: 2965.[CrossRef][Medline]
Fernández C. and Wüthrich K. 2003. NMR solution structure determination of membrane proteins reconstituted in detergent micelles FEBS Lett. 555: 144150.[CrossRef][Medline]
Fernández C., Hilty C., Wider G., Wüthrich K. 2002. Lipidprotein interactions in DHPC micelles containing the integral membrane protein OmpX investigated by NMR spectroscopy Proc. Natl. Acad. Sci. 99: 1353313537.
Frueh D.P., Ito T., Li T.S., Wagner G., Glaser S.J., Khaneja N. 2005. Sensitivity enhancement in NMR of macromolecules by application of optimal control theory J. Biomol. NMR 32: 2330.[Medline]
Fu R. and Cross T.A. 1999. Solid-state nuclear magnetic resonance investigation of protein and polypeptide structure. Annu. Rev. Biophys. Biomol Struct. 28: 235268.
Gardner K.H. and Kay L.E. 1998. The use of 2H, 13C, 15N multidimensional NMR to study the structure and dynamics of proteins. Annu. Rev. Biophys. Biomol Struct. 27: 357406.
Geen H. and Freeman R. 1991. Band-selective radiofrequency pulses J. Magn. Reson. 93: 93141.
Gray T.M. and Matthews B.W. 1984. Intrahelical hydrogen bonding of serine, threonine and cysteine residues within a-helices and its relevance to membrane-bound proteins J. Mol. Biol. 175: 7581.[CrossRef][Medline]
Grzesiek S., Anglister J., Ren H., Bax A. 1993. 13C line-narrowing by 2H decoupling in 2H/13C/15N enriched proteins. Application to triple resonance 4D J-connectivity of sequential amides J. Am. Chem. Soc. 115: 43694370.
Gulbis J.M. and Doyle D.A. 2004. Potassium channel structures: Do they conform? Curr. Opin. Struct. Biol. 14: 440446.[CrossRef][Medline]
Heginbotham L., Lu Z., Abramson T., MacKinnon R. 1994. Mutations in the K+channel signature sequence Biophys. J. 66: 10611067.[Medline]
Heginbotham L., Odessey E., Miller C. 1997. Tetrameric stoichiometry of a prokaryotic K+ channel Biochemistry 36: 1033510342.[CrossRef][Medline]
Heginbotham L., LeMasurier M., Kolmakova-Partensky L., Miller C. 1999. Single streptomyces lividans K+channels: Functional asymmetries sidedness of proton activation J. Gen. Physiol. 114: 551559.
Hille B. In Ion channels of excitable membranes . 2001. 3rd ed. Sinauer Associates, Sunderland, MA.
Hilty C., Wider G., Fernández C., Wüthrich K. 2004. Membrane proteinlipid interactions in mixed micelles studies by NMR spectroscopy with the use of paramagnetic reagents Chembiochem. 5: 467473.[CrossRef][Medline]
Ho K., Nichols C.G., Lederer W.J., Lytton J., Vassiliev P.M., Kanazirska M.V., Herbert S.C. 1993. Cloning and expression of an inwardly rectifying ATP-regulated potassium channel Nature 362: 3138.[CrossRef][Medline]
Hwang P.M., Choy W.-Y., Lo E.I., Forman-Kay J.D., Raetz C.R.H., Prive G.G., Bishop R.E., Kay L.E. 2002. Solution structure and dynamics of the outer m