Interactions between charged residues in the transmembrane segments of the voltage-sensing domain in the hERG channel.

Studies on voltage-gated K channels such as Shaker have shown that positive charges in the voltage-sensor (S4) can form salt bridges with negative charges in the surrounding transmembrane segments in a state-dependent manner, and different charge pairings can stabilize the channels in closed or open states. The goal of this study is to identify such charge interactions in the hERG channel. This knowledge can provide constraints on the spatial relationship among transmembrane segments in the channel's voltage-sensing domain, which are necessary for modeling its structure. We first study the effects of reversing S4's positive charges on channel activation. Reversing positive charges at the outer (K525D) and inner (K538D) ends of S4 markedly accelerates hERG activation, whereas reversing the 4 positive charges in between either has no effect or slows activation. We then use the 'mutant cycle analysis' to test whether D456 (outer end of S2) and D411 (inner end of S1) can pair with K525 and K538, respectively. Other positive charges predicted to be able, or unable, to interact with D456 or D411 are also included in the analysis. The results are consistent with predictions based on the distribution of these charged residues, and confirm that there is functional coupling between D456 and K525 and between D411 and K538.

Structural models of the MscL gating mechanism.

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

A putative prokaryote voltage-gated Ca(2+) channel with only one 6TM motif per subunit.

Until now, voltage-gated Ca(2+) channel proteins have been found only in eukaryotes. Here we report that a gene recently discovered in the eubacterium Bacillus halodurans codes for a protein closely related to eukaryotic Ca(2+) channels, but that has only one 6-transmembrane-segement (6TM) motif, instead of four, in its pore-forming subunit. This is supported by the comparison of consensus sequences, which, along with the patterns of residue conservation, indicates a similar structure in the membrane to voltage-gated K(+) channels. From this we hypothesize that Ca(2+) channels originally evolved in bacteria, and that the specific eubacteria protein highlighted here is an ideal candidate for structure determination efforts.

The gating mechanism of the large mechanosensitive channel MscL.

The mechanosensitive channel of large conductance, MscL, is a ubiquitous membrane-embedded valve involved in turgor regulation in bacteria. The crystal structure of MscL from Mycobacterium tuberculosis provides a starting point for analysing molecular mechanisms of tension-dependent channel gating. Here we develop structural models in which a cytoplasmic gate is formed by a bundle of five amino-terminal helices (S1), previously unresolved in the crystal structure. When membrane tension is applied, the transmembrane barrel expands and pulls the gate apart through the S1-M1 linker. We tested these models by substituting cysteines for residues predicted to be near each other only in either the closed or open conformation. Our results demonstrate that S1 segments form the bundle when the channel is closed, and crosslinking between S1 segments prevents opening. S1 segments interact with M2 when the channel is open, and crosslinking of S1 to M2 impedes channel closing. Gating is affected by the length of the S1-M1 linker in a manner consistent with the model, revealing critical spatial relationships between the domains that transmit force from the lipid bilayer to the channel gate.

A family of putative Kir potassium channels in prokaryotes.

BACKGROUND: Prior to this report, members of the inward rectifier family, or Kir, have been found only in eukaryotes. Like most K+ channels, the pore-forming part of the protein is formed by four identical, or closely related, subunits. Each subunit contains a transmembrane M1-P-M2 motif that is followed by a relatively large C-terminus region unique to Kir's. RESULTS: In searching unfinished microbial genomes for K+ channels, we identified five sequences in the prokaryote Burkholderia pseudomallei, Burkholderia cepacia, Burkholderia fungorum LB400, Magentospirillum magnetotacticum, and Nostoc Punctiforme genomes that code for proteins whose closest relatives in current sequence databases are eukaryote Kir's. The sequence similarity includes the C-terminus portion of Kir's, for which there are no other close homologs in current prokaryote sequences. Sequences of the pore-forming P and M2 segments of these proteins, which we call KirBac, is intermediate between those of eukaryotic Kir's and several other K+ channel families. CONCLUSIONS: Although KirBac's are more closely related to Kir's than to other families of K+ channels, the intermediate nature of their pore-forming P and M2 segments suggests that they resemble an ancestral precursor to the eukaryotic Kir's. The similarity of KirBac to the bacterial KcsA channel, whose transmembrane structure has been solved, helps align Kir's with KcsA. KirBac's may assist in solving the three-dimensional structure of a member of the Kir family since bacterial membrane proteins are more easily expressed in the quantities necessary for crystallography.

Two critical cysteine residues implicated in disulfide bond formation and proper folding of Kir2.1.

Inwardly rectifying potassium channels are important in cellular repolarization of many excitable tissues. Amino acid sequence alignment of different mammalian inward rectifier K(+) channels revealed two absolutely conserved cysteine residues in the putative extracellular face, suggesting a possible disulfide bond. Replacement of these cysteine residues in the Kir2.1 channel (i.e., C122 and C154) with either alanine or serine abolished current in Xenopus laevis oocytes although Western blotting established that the channels were fully expressed. The digestion pattern of channels treated with V8 protease combined with Western blotting under reducing and nonreducing conditions confirmed intrasubunit cross-linking of C122 and C154. Whole-cell and single channel current recordings of oocytes expressing tandem tetrameric constructs with one or two of the mutant subunits suggested that insertion of one mutant subunit is sufficient to eliminate channel function. Coexpression studies confirmed that the cysteine mutant channels eliminate wild-type Kir2.1 currents in a dominant-negative manner. Despite these results, sulfhydryl reduction did not alter the functional properties of Kir2.1 currents. Molecular modeling of Kir2.1 with the two cysteines cross-linked predicted that the extracellular loop between the first transmembrane domain and the pore helix contains a beta-hairpin structure. Distinct from the KcsA structure, the disulfide bond together with the beta-hairpin structure is expected to constrain and stabilize the P-loop and selectivity filter. Taken together, these results suggest that intramolecular disulfide bond exists between C122 and C154 of Kir2.1 channel and this cross-link might be required for proper channel folding.

Does the KdpA subunit from the high affinity K(+)-translocating P-type KDP-ATPase have a structure similar to that of K(+) channels?

Evidence is presented that the transmembrane KdpA subunit of the high affinity K(+)-translocating P-type Kdp-ATPase is evolutionarily derived from the superfamily of 2TM-type K(+) channels in bacteria. This extends a previous study relating the K(+) channels to the KtrAB, Trk, Trk1,2, and HKT1 K(+) symporter superfamily of both prokaryotes and eukaryotes. Although the channels are formed by four single-MPM motif subunits, the transmembrane KdpA subunit and the transmembrane subunit of the symporter proteins are postulated to have four corresponding MPM motifs within a single sequence. Analysis of 17 KdpA sequences reveals a pattern of residue conservation similar to that of the symporters and channels, and consistent with the crystal structure of the KcsA K(+) channel. In addition, the most highly conserved residues between the families, specifically the central glycines of the P2 segments, are those previously identified as crucial for the property of K(+)-selectivity that is common to each protein. This hypothesis is consistent with an experimental study of mutations that alter K(+) binding affinity of the Kdp transporter. Although most of the results of a previous study of the transmembrane topology of KdpA are consistent with the 4-MPM model, the one deviation can be explained by a plausible change in the structure due to the experimental method.

Evolutionary relationship between K(+) channels and symporters.

The hypothesis is presented that at least four families of putative K(+) symporter proteins, Trk and KtrAB from prokaryotes, Trk1,2 from fungi, and HKT1 from wheat, evolved from bacterial K(+) channel proteins. Details of this hypothesis are organized around the recently determined crystal structure of a bacterial K(+) channel: i. e., KcsA from Streptomyces lividans. Each of the four identical subunits of this channel has two fully transmembrane helices (designated M1 and M2), plus an intervening hairpin segment that determines the ion selectivity (designated P). The symporter sequences appear to contain four sequential M1-P-M2 motifs (MPM), which are likely to have arisen from gene duplication and fusion of the single MPM motif of a bacterial K(+) channel subunit. The homology of MPM motifs is supported by a statistical comparison of the numerical profiles derived from multiple sequence alignments formed for each protein family. Furthermore, these quantitative results indicate that the KtrAB family of symporters has remained closest to the single-MPM ancestor protein. Strong sequence evidence is also found for homology between the cytoplasmic C-terminus of numerous bacterial K(+) channels and the cytoplasm-resident TrkA and KtrA subunits of the Trk and KtrAB symporters, which in turn are homologous to known dinucleotide-binding domains of other proteins. The case for homology between bacterial K(+) channels and the four families of K(+) symporters is further supported by the accompanying manuscript, in which the patterns of residue conservation are demonstrated to be similar to each other and consistent with the known 3D structure of the KcsA K(+) channel.

Structural models of the KtrB, TrkH, and Trk1,2 symporters based on the structure of the KcsA K(+) channel.

Three-dimensional computer modeling is used to further investigate the hypothesis forwarded in the accompanying paper of an evolutionary relationship between four related families of K(+) sympoter proteins and the superfamily of K(+) channel proteins. Atomic-scale models are developed for the transmembrane regions of one member from each of the three more distinct symporter families, i.e., a TrkH protein from Escherichia coli, a KtrB protein from Aquifex aeolicus, and a Trk1,2 protein from Schizosaccharomyces pombe. The portions of the four consecutive M1-P-M2 motifs in the symporters that can be aligned with K(+) channel sequences are modeled directly from the recently determined crystal structure of the KcsA K(+) channel from Streptomyces lividans. The remaining portions are developed using our previously accumulated theoretical modeling criteria and principles. Concurrently, the use of these criteria and principles is further supported by the now verified predictions of our previous K(+) channel modeling efforts and the degree to which they are satisfied by the known structure of the KcsA protein. Thus the observed ability of the portions of the symporter models derived from the KcsA crystal structure to also satisfy the theoretical modeling criteria provides additional support for an evolutionary link with K(+) channel proteins. Efforts to further satisfy the criteria and principles suggest that the symporter proteins from fungi and plants (i.e., Trk1,2 and HKT1) form dimeric and/or tetrameric complexes in the membrane. Furthermore, analysis of the atomic-scale models in relation to the sequence conservation within and between the protein families suggests structural details for previously proposed mechanisms for the linked symport of K(+) with Na(+) and H(+). Suggestions are also given for experiments to test these structures and hypotheses.

Epitope mapping of CCR5 reveals multiple conformational states and distinct but overlapping structures involved in chemokine and coreceptor function.

The chemokine receptor CCR5 is the major coreceptor for R5 human immunodeficiency virus type-1 strains. We mapped the epitope specificities of 18 CCR5 monoclonal antibodies (mAbs) to identify domains of CCR5 required for chemokine binding, gp120 binding, and for inducing conformational changes in Env that lead to membrane fusion. We identified mAbs that bound to N-terminal epitopes, extracellular loop 2 (ECL2) epitopes, and multidomain (MD) epitopes composed of more than one single extracellular domain. N-terminal mAbs recognized specific residues that span the first 13 amino acids of CCR5, while nearly all ECL2 mAbs recognized residues Tyr-184 to Phe-189. In addition, all MD epitopes involved ECL2, including at least residues Lys-171 and Glu-172. We found that ECL2-specific mAbs were more efficient than NH2- or MD-antibodies in blocking RANTES or MIP-1beta binding. By contrast, N-terminal mAbs blocked gp120-CCR5 binding more effectively than ECL2 mAbs. Surprisingly, ECL2 mAbs were more potent inhibitors of viral infection than N-terminal mAbs. Thus, the ability to block virus infection did not correlate with the ability to block gp120 binding. Together, these results imply that chemokines and Env bind to distinct but overlapping sites in CCR5, and suggest that the N-terminal domain of CCR5 is more important for gp120 binding while the extracellular loops are more important for inducing conformational changes in Env that lead to membrane fusion and virus infection. Measurements of individual antibody affinities coupled with kinetic analysis of equilibrium binding states also suggested that there are multiple conformational states of CCR5. A previously described mAb, 2D7, was unique in its ability to effectively block both chemokine and Env binding as well as coreceptor activity. 2D7 bound to a unique antigenic determinant in the first half of ECL2 and recognized a far greater proportion of cell surface CCR5 molecules than the other mAbs examined. Thus, the epitope recognized by 2D7 may represent a particularly attractive target for CCR5 antagonists.

Structural models of the transmembrane region of voltage-gated and other K+ channels in open, closed, and inactivated conformations.

A large collaborative, multidisciplinary effort involving many research laboratories continues which uses indirect methods of molecular biology and membrane biophysics to analyze the three-dimensional structures and functional mechanisms of K+ channels. This work also extends to the distant relatives of these channels, including the voltage-gated Na+ and Ca2+ channels. The role that our group plays in this process is to combine the information gained from experimental studies with molecular modeling techniques to generate atomic-scale structural models of these proteins. The modeling process involves three stages which are summarized as: (I) prediction of the channel sequence transmembrane topology, including the functionality and secondary structure of the segments; (II) prediction of the relative positions of the transmembrane segments, and (III) filling in all atoms of the amino acid residues, with conformations for energetically stabilized interactions. Both physiochemical and evolutionary principles (including sequence homology analysis) are used to guide the development. In addition to testing the steric and energetic feasibilities of different structural hypotheses, the models provide guidance for the design of new experiments. Structural modeling also serves to "fill in the gaps" of experimental data, such as predicting additional residue interactions and conformational changes responsible for functional processes. The modeling process is currently at the stage that experimental studies have definitely confirmed most of our earlier predictions about the transmembrane topology and functionality of different segments. Additionally, this report describes the detailed, three-dimensional models we have developed for the entire transmembrane region and important functional sites of the voltage-gated Shaker K+ channel in the open, closed, and inactivated conformations (including the ion-selective pore and voltage-sensor regions). As part of this effort, we also describe how our development of structural models for many of the other major K+ channel families aids in determining common structural motifs. As an example, we also present a detailed model of the smaller, bacterial K+ channel from Streptomyces lividans. Finally, we discuss strategies for using newly developed experimental methods for determining the structures and analyzing the functions of these channel proteins.

NMR structure of inactivation gates from mammalian voltage-dependent potassium channels.

The electrical signalling properties of neurons originate largely from the gating properties of their ion channels. N-type inactivation of voltage-gated potassium (Kv) channels is the best-understood gating transition in ion channels, and occurs by a 'ball-and-chain' type mechanism. In this mechanism an N-terminal domain (inactivation gate), which is tethered to the cytoplasmic side of the channel protein by a protease-cleavable chain, binds to its receptor at the inner vestibule of the channel, thereby physically blocking the pore. Even when synthesized as a peptide, ball domains restore inactivation in Kv channels whose inactivation domains have been deleted. Using high-resolution nuclear magnetic resonance (NMR) spectroscopy, we analysed the three-dimensional structure of the ball peptides from two rapidly inactivating mammalian K. channels (Raw3 (Kv3.4) and RCK4 (Kv1.4)). The inactivation peptide of Raw3 (Raw3-IP) has a compact structure that exposes two phosphorylation sites and allows the formation of an intramolecular disulphide bridge between two spatially close cysteine residues. Raw3-IP exhibits a characteristic surface charge pattern with a positively charged, a hydrophobic, and a negatively charged region. The RCK4 inactivation peptide (RCK4-IP) shows a similar spatial distribution of charged and uncharged regions, but is more flexible and less ordered in its amino-terminal part.

Membrane topology and multimeric structure of a mechanosensitive channel protein of Escherichia coli.

We have studied the membrane topology and multimeric structure of a mechanosensitive channel, MscL, which we previously isolated and cloned from Escherichia coli. We have localized this 15-kDa protein to the inner membrane and, by PhoA fusion, have shown that it contains two transmembrane domains with both the amino and carboxyl termini on the cytoplasmic side. Mutation of the glutamate at position 56 to histidine led to changes in channel kinetics which were dependent upon the pH on the periplasmic, but not cytoplasmic side of the membrane, providing additional evidence for the periplasmic positioning of this part of the molecule. Tandems of two MscL subunits expressed as a single polypeptide formed functional channels, suggesting an even number of transmembrane domains per subunit (amino and carboxyl termini on the same side of the membrane), and an even number of subunits per functional complex. Finally, cross-linking studies suggest that the functional MscL complex is a homohexamer. In summary, these data are all consistent with a protein domain assignment and topological model which we propose and discuss.

Structural model of the outer vestibule and selectivity filter of the Shaker voltage-gated K+ channel.

A new generation of structural models were developed of the outer vestibule and ion-selective portion of the voltage-gated Shaker K+ channel. Some features of these models are similar to those that we have developed previously [Durrel S. R. and Guy H. R. (1992) Biophys. J. 62, 238-250; Guy H. R. (1990) In Monovalent Cations in Biological Systems (Pasternak C. A., Ed.), pp. 31-58, CRC Press, Boca Raton, FL; Guy H. R. and Durell S. R. (1994) In Molecular Evolution of Physiological processes (Fambrough D., Ed.), pp. 197-212, The Rockefeller University Press, NY; Guy H. R. and Durell S. R. (1995) In Ion Channels and Genetic Diseases (Dawson D., Ed.), pp. 1-16, The Rockefeller University Press, NY] and other features were modified to make the models more consistent with recent experimental findings. The first part of the P segment is postulated, as always, to form a short alpha helix that spans only the outer portion of the membrane. The helix is tilted so that its C-terminal is nearer the pore than its N-terminal. The latter part of the P segment, P2, is postulated to have a relatively elongated conformation that is positioned approximately parallel to the axis of the pore. Four of the P2 segments assemble to form an ion-selective region that has two narrow regions; one formed by the Y445 side-chains at the outer entrance of the pore and one formed by the backbone of the T442 residues near the innermost part of the P segments. The S6 segment is postulated to form two alpha helices. The first S6 helix packs next to the P segments in our models. The NMR structures of two scorpion toxins, charybdotoxin and agitoxin 2, have been docked into the models of the outer vestibules. The shape of the outer vestibule has been modeled so that specific toxin-channel residue-residue interactions correspond to those that have been identified experimentally.

The S. cerevisiae outwardly-rectifying potassium channel (DUK1) identifies a new family of channels with duplicated pore domains.

Potassium channel subunits have six or two transmembrane segments in addition to a conserved pore-forming (P) domain; four subunits come together to form a channel. A gene was identified in S. cerevisiae (J0911) encoding a protein with eight probable membrane-spanning segments and two such P regions. This protein (Duk1p) is a potassium channel because Xenopus oocytes injected with the corresponding RNA express potassium currents activated by depolarization that are not seen in control oocytes. Similar potassium currents were recorded from wildtype S. cerevisiae spheroplasts, but not from those in which the DUK1 locus had been disrupted. Cells carrying the duk1 delta 1::HIS disruption in addition to a chimeric gene comprising DUK1 behind the GAL1 promoter showed outward currents when grown in galactose, but not when grown in glucose. Additional sequences with the duplicate pore motif were found in C. elegans, suggesting that these proteins represent a novel structural family of potassium channel proteins.

min K channels form by assembly of at least 14 subunits.

Injection of min K mRNA into Xenopus oocytes results in expression of slowly activating voltage-dependent potassium channels, distinct from those induced by expression of other cloned potassium channels. The min K protein also differs in structure, containing only a single predicted transmembrane domain. While it has been demonstrated that all other cloned potassium channels form by association of four independent subunits, the number of min K monomers which constitute a functional channel is unknown. In rat min K, replacement of Ser-69 by Ala (S69A) causes a shift in the current-voltage (I-V) relationship to more depolarized potentials; currents are not observed at potentials negative to 0 mV. To determine the subunit stoichiometry of min K channels, wild-type and S69A subunits were coexpressed. Injections of a constant amount of wild-type mRNA with increasing amounts of S69A mRNA led to potassium currents of decreasing amplitude upon voltage commands to -20 mV. Applying a binomial distribution to the reduction of current amplitudes as a function of the different coinjection mixtures yielded a subunit stoichiometry of at least 14 monomers for each functional min K channel. A model is presented for how min K subunits may form a channel.

Structural models of Na+, Ca2+, and K+ channels.

The structure determination of voltage-gated channels by the combination of molecular modeling and mutagenesis experiments is a long term, iterative process. As such, the models should be considered as a work in progress, with changes expected as more data becomes available. The primary role of the models is that they assimilate the known data and provide ideas for further experiments to elucidate the real structures. Although the models presented here have already gone through two or three cycles of development and testing, many aspects remain tentative. Perhaps the most significant result so far is that the P segment was experimentally confirmed to form the ion-selective part of the channel. In a subsequent cycle of testing and modeling, the specific residues responsibility for Na+ and Ca2+ selectivity have been identified and the selectivity filter of K+ channels is now predicted to be formed by the side chains, but rather by the carbonyl oxygens of the conserved Gly-Tyr-Gly sequence backbone. As another example, the 9p residue of the P segment of K+ channels was originally modeled as either being buried in the protein or accessible channels was originally modeled as either being buried in the protein or accessible from inside the cell only. However, once mutation of this residue to histidine was found to affect blockade by extracellular TEA, protons, Zn2+ and histidine reagents (DeBiasi et al., 1993), the models were updated to have this and the hydrophilic residues in the first part of P form a helix that comprises part of the extracellular, outer vestibular of the pore. While this motif was used also for Na+ and Ca2+ pore models (see Fig. 2) where the putative helices are amphipathic, it remains to be verified. Modeling of the size and shape of the outer vestibule of K+ channels was also aided by the data for the binding of CTX in the extracellular entrance to the pore. Similarly, experiments with peptide toxins such as mu and omega conotoxins may prove useful in modeling the outer vestibules of the Na+ and/or Ca2+ channels. While important advances have been made, it is important to realize that these approaches are still very new. In the future we are likely to see improvements on both the theoretical and experimental sides which will greatly advance the process.(ABSTRACT TRUNCATED AT 400 WORDS)

Theoretical models of the ion channel structure of amyloid beta-protein.

Theoretical methods are used to develop models for the ion channel structure of the membrane-bound amyloid beta-protein. This follows recent observations that the beta-protein forms cation-selective channels in lipid bilayers in vitro. Amyloid beta-protein is the main component of the extracellular plaques in the brain that are characteristic of Alzheimer's disease. Based on the amino acid sequence and the unique environment of the membrane, the secondary structure of the 40-residue beta-protein is predicted to form a beta-hairpin followed by a helix-turn-helix motif. The channel structures were-designed as aggregates of peptide subunits in identical conformations. Three types of models were developed that are distinguished by whether the pore is formed by the beta-hairpins, the middle helices, or by the more hydrophobic C-terminal helices. The latter two types can be converted back and forth by a simple conformational change, which would explain the variable conduction states observed for a single channel. It is also demonstrated how lipid headgroups could be incorporated into the pore lining, and thus affect the ion selectivity. The atomic-scale detail of the models make them useful for designing experiments to determine the real structure of the channel, and thus further the understanding of peptide channels in general. In addition, if beta-protein-induced channel activity is found to be the cause of cell death in Alzheimer's disease, then the models may be helpful in designing counteracting drugs.