Structural, Biochemical, and Functional Characterization of the Cyclic Nucleotide Binding Homology Domain from the Mouse EAG1 Potassium Channel

KCNH channels are voltage-gated potassium channels with important physiological functions. In these channels, a C-terminal cytoplasmic region, known as the cyclic nucleotide binding homology (CNB-homology) domain displays strong sequence similarity to cyclic nucleotide binding (CNB) domains. However, the isolated domain does not bind cyclic nucleotides. Here, we report the X-ray structure of the CNB-homology domain from the mouse EAG1 channel. Through comparison with the recently determined structure of the CNB-homology domain from the zebrafish ELK (eag-like K(+)) channel and the CNB domains from the MlotiK1 and HCN (hyperpolarization-activated cyclic nucleotide-gated) potassium channels, we establish the structural features of CNB-homology domains that explain the low affinity for cyclic nucleotides. Our structure establishes that the "self-liganded" conformation, where two residues of the C-terminus of the domain are bound in an equivalent position to cyclic nucleotides in CNB domains, is a conserved feature of CNB-homology domains. Importantly, we provide biochemical evidence that suggests that there is also an unliganded conformation where the C-terminus of the domain peels away from its bound position. A functional characterization of this unliganded conformation reveals a role of the CNB-homology domain in channel gating.     

Predicted structure of the extracellular region of ligand-gated ion-channel receptors shows SH2-like and SH3-like domains forming the ligand-binding site.

Fast synaptic neurotransmission is mediated by ligand-gated ion-channel (LGIC) receptors, which include receptors for acetylcholine, serotonin, GABA, glycine, and glutamate. LGICs are pentamers with extracellular ligand-binding domains and form integral membrane ion channels that are selective for cations (acetylcholine and serotonin 5HT3 receptors) or anions (GABAA and glycine receptors and the invertebrate glutamate-binding chloride channel). They form a protein superfamily with no sequence similarity to any protein of known structure. Using a 1D-3D structure mapping approach, we have modeled the extracellular ligand-binding domain based on a significant match with the SH2 and SH3 domains of the biotin repressor structure. Refinement of the model based on knowledge of the large family of SH2 and SH3 structures, sequence alignments, and use of structure templates for loop building, allows the prediction of both monomer and pentamer models. These are consistent with medium-resolution electron microscopy structures and with experimental structure/function data from ligand-binding, antibody-binding, mutagenesis, protein-labeling and subunit-linking studies, and glycosylation sites. Also, the predicted polarity of the channel pore calculated from electrostatic potential maps of pentamer models of superfamily members is consistent with known ion selectivities. Using the glycine receptor alpha 1 subunit, which forms homopentamers, the monomeric and pentameric models define the agonist and antagonist (strychnine) binding sites to a deep crevice formed by an extended loop, which includes the invariant disulfide bridge, between the SH2 and SH3 domains. A detailed binding site for strychnine is reported that is in strong agreement with known structure/function data. A site for interaction of the extracellular ligand-binding domain with the activation of the M2 transmembrane helix is also suggested.     

Short-range molecular rearrangements in ion channels detected by tryptophan quenching of bimane fluorescence

Ion channels are allosteric membrane proteins that open and close an ion-permeable pore in response to various stimuli. This gating process provides the regulation that underlies electrical signaling events such as action potentials, postsynaptic potentials, and sensory receptor potentials. Recently, the molecular structures of a number of ion channels and channel domains have been solved by x-ray crystallography. These structures have highlighted a gap in our understanding of the relationship between a channel's function and its structure. Here we introduce a new technique to fill this gap by simultaneously measuring the channel function with the inside-out patch-clamp technique and the channel structure with fluorescence spectroscopy. The structure and dynamics of short-range interactions in the channel can be measured by the presence of quenching of a covalently attached bimane fluorophore by a nearby tryptophan residue in the channel. This approach was applied to study the gating rearrangements in the bovine rod cyclic nucleotide-gated ion channel CNGA1 where it was found that C481 moves towards A461 during the opening allosteric transition induced by cyclic nucleotide. The approach offers new hope for elucidating the gating rearrangements in channels of known structure.     

Crystal structure of the atypical protein kinase domain of a TRP channel with phosphotransferase activity

Transient receptor potential (TRP) channels modulate calcium levels in eukaryotic cells in response to external signals. A novel transient receptor potential channel has the ability to phosphorylate itself and other proteins on serine and threonine residues. The catalytic domain of this channel kinase has no detectable sequence similarity to classical eukaryotic protein kinases and is essential for channel function. The structure of the kinase domain, reported here, reveals unexpected similarity to eukaryotic protein kinases in the catalytic core as well as to metabolic enzymes with ATP-grasp domains. The inclusion of the channel kinase catalytic domain within the eukaryotic protein kinase superfamily indicates a significantly wider distribution for this group of signaling proteins than suggested previously by sequence comparisons alone.     

A voltage-sensor water pore

Voltage-sensor (VS) domains cause the pore of voltage-gated ion channels to open and close in response to changes in transmembrane potential. Recent experimental studies suggest that VS domains are independent structural units. This independence is revealed dramatically by a voltage-dependent proton-selective channel (Hv), which has a sequence homologous to the VS domains of voltage-gated potassium channels (Kv). Here we show by means of molecular dynamics simulations that the isolated open-state VS domain of the KvAP channel in a lipid membrane has a configuration consistent with a water channel, which we propose as a common feature underlying the conductance of protons, and perhaps other cations, through VS domains.     

Three-dimensional structure of a voltage-gated potassium channel at 2.5 nm resolution

BACKGROUND: The voltage-gated potassium channel Shaker from Drosophila consists of a tetramer of identical subunits, each containing six transmembrane segments. The atomic structure of a bacterial homolog, the potassium channel KcsA, is much smaller than Shaker. It does not have a voltage sensor and other important domains like the N-terminal tetramerization (T1) domain. The structure of these additional elements has to be studied in the more complex voltage-gated channels. RESULTS: We determined the three-dimensional structure of the entire Shaker channel at 2.5 nm resolution using electron microscopy. The four-fold symmetric structure shows a large and a small domain linked by thin 2 nm long connectors. To interpret the structure, we used the crystal structures of the isolated T1 domain and the KcsA channel. A unique density assignment was made based on the symmetry and dimensions of the crystal structures and domains, identifying the smaller domain as the cytoplasmic mass of Shaker containing T1 and the larger domain as embedded in the membrane. CONCLUSIONS: The two-domain architecture of the Shaker channel is consistent with the recently proposed "hanging gondola" model for the T1 domain, putting the T1 domain at a distance from the membrane domain but attached to it by thin connectors. The space between the two domains is sufficient to permit cytoplasmic access of ions and the N-terminal inactivation domain to the pore region. A hanging gondola architecture has also been observed in the nicotinic acetylcholine receptor and the KcsA structure, suggesting that it is a common element of ion channels.     

Comparative models of P2X2 receptor support inter-subunit ATP-binding sites

ATP-gated P2X receptors (P2XRs) are ligand-gated ion channels (LGICs) presumably trimeric. To date, no experimental high-resolution structures are available. Recent X-ray structure of the acid-sensing ion channel 1 (ASIC1) revealed an unexpected trimeric ion channel. Beside their quaternary structure, P2XR and ASIC1 share common membrane topologies, but no significant sequence similarity. In order to overcome this low sequence resemblance, we have developed comparative models of P2X₂R based on secondary structure predictions using the crystal structure of ASIC1 as template. These models were constrained to be consistent with known arrangement of disulfide bridges. They agreed with cross-linking experiments and supported inter-subunit ATP-binding sites. One of our models reconciled most existing data and provides new structural insights for a plausible mechanism of gating, thus encouraging new experiments.     

Domain architecture comparison for multidomain homology identification.

Homology identification is the first step for many genomic studies. Current methods, based on sequence comparison, can result in a substantial number of mis-assignments due to the similarity of homologous domains in otherwise unrelated sequences. Here we propose methods to detect homologs through explicit comparison of protein domain content. We developed several schemes for scoring the homology of a pair of protein sequences based on methods used in the field of information retrieval. We evaluate the proposed methods and methods used in the literature using a benchmark of fifteen sequence families of known evolutionary history. The results of these studies demonstrate the effectiveness of comparing domain architectures using these similarity measures. We also demonstrate the importance of both weighting promiscuous domains and of compensating for the statistical effect of having a large number of domains in a protein. Using logistic regression, we demonstrate the benefit of combining similarity measures based on domain content with sequence similarity measures.     

Structure,Dynamics and Implied Gating Mechanism of a Human Cyclic Nucleotide-Gated Channel

Cyclic nucleotide-gated (CNG) ion channels are nonselective cation channels, essential for visual and olfactory sensory transduction. Although the channels include voltage-sensor domains (VSDs), their conductance is thought to be independent of the membrane potential, and their gating regulated by cytosolic cyclic nucleotide–binding domains. Mutations in these channels result in severe, degenerative retinal diseases, which remain untreatable. The lack of structural information on CNG channels has prevented mechanistic understanding of disease-causing mutations, precluded structure-based drug design, and hampered in silico investigation of the gating mechanism. To address this, we built a 3D model of the cone tetrameric CNG channel, based on homology to two distinct templates with known structures: the transmembrane (TM) domain of a bacterial channel, and the cyclic nucleotide-binding domain of the mouse HCN2 channel. Since the TM-domain template had low sequence-similarity to the TM domains of the CNG channels, and to reconcile conflicts between the two templates, we developed a novel, hybrid approach, combining homology modeling with evolutionary coupling constraints. Next, we used elastic network analysis of the model structure to investigate global motions of the channel and to elucidate its gating mechanism. We found the following: (i) In the main mode of motion, the TM and cytosolic domains counter-rotated around the membrane normal. We related this motion to gating, a proposition that is supported by previous experimental data, and by comparison to the known gating mechanism of the bacterial KirBac channel. (ii) The VSDs could facilitate gating (supplementing the pore gate), explaining their presence in such ‘voltage-insensitive’ channels. (iii) Our elastic network model analysis of the CNGA3 channel supports a modular model of allosteric gating, according to which protein domains are quasi-independent: they can move independently, but are coupled to each other allosterically.     

mu-conotoxin GIIIA interactions with the voltage-gated Na(+) channel predict a clockwise arrangement of the domains

Voltage-gated Na(+) channels underlie the electrical activity of most excitable cells, and these channels are the targets of many antiarrhythmic, anticonvulsant, and local anesthetic drugs. The channel pore is formed by a single polypeptide chain, containing four different, but homologous domains that are thought to arrange themselves circumferentially to form the ion permeation pathway. Although several structural models have been proposed, there has been no agreement concerning whether the four domains are arranged in a clockwise or a counterclockwise pattern around the pore, which is a fundamental question about the tertiary structure of the channel. We have probed the local architecture of the rat adult skeletal muscle Na(+) channel (mu1) outer vestibule and selectivity filter using mu-conotoxin GIIIA (mu-CTX), a neurotoxin of known structure that binds in this region. Interactions between the pore-forming loops from three different domains and four toxin residues were distinguished by mutant cycle analysis. Three of these residues, Gln-14, Hydroxyproline-17 (Hyp-17), and Lys-16 are arranged approximately at right angles to each other in a plane above the critical Arg-13 that binds directly in the ion permeation pathway. Interaction points were identified between Hyp-17 and channel residue Met-1240 of domain III and between Lys-16 and Glu-403 of domain I and Asp-1532 of domain IV. These interactions were estimated to contribute -1.0+/-0.1, -0.9+/-0.3, and -1.4+/-0.1 kcal/mol of coupling energy to the native toxin-channel complex, respectively. mu-CTX residues Gln-14 and Arg-1, both on the same side of the toxin molecule, interacted with Thr-759 of domain II. Three analytical approaches to the pattern of interactions predict that the channel domains most probably are arranged in a clockwise configuration around the pore as viewed from the extracellular surface.     

Domain size distributions can predict domain boundaries

MOTIVATION: The sizes of protein domains observed in the 3D-structure database follow a surprisingly narrow distribution. Structural domains are furthermore formed from a single-chain continuous segment in over 80% of instances. These observations imply that some choices of domain boundaries on an otherwise uncharacterized sequence are more likely than others, based solely on the size and segment number of predicted domains. This property might be used to guess the locations of protein domain boundaries. RESULTS: To test this possibility we enumerate putative domain boundaries and calculate their relative likelihood under a probability model that considers only the size and segment number of predicted domains. We ask, in a cross-validated test using sequences with known 3D structure, whether the most likely guesses agree with the observed domain structure. We find that domain boundary predictions are surprisingly successful for sequences up to 400 residues long and that guessing domain boundaries in this way can improve the sensitivity of threading analysis.     

Structural order in Pannexin 1 cytoplasmic domains

Pannexin 1 forms ion and metabolite permeable hexameric channels with abundant expression in the central nervous system and elsewhere. Although pannexin 1 does not form intercellular channels, a common channel topology and oligomerization state, as well as involvement of the intracellular carboxyl terminal (CT) domain in channel gating, is shared with connexins. In this study, we characterized the secondary structure of the mouse pannexin 1 cytoplasmic domains to complement structural studies of the transmembrane segments and compare with similar domains from connexins. A combination of structural prediction tools and circular dichroism revealed that, unlike connexins (predominately intrinsically disordered), cytosolic regions of pannexin 1 contain approximately 50% secondary structure, a majority being α-helical. Moreover, prediction of transmembrane domains uncovered a potential membrane interacting region (I360-G370) located upstream of the caspase cleavage site (D375-D378) within the pannexin 1 CT domain. The α-helical content of a peptide containing these domains (G357-S384) increased in the presence of detergent micelles providing evidence of membrane association. We also purified a pannexin 1 CT construct containing the caspase cleavage site (M374-C426), assigned the resonances by NMR, and confirmed cleavage by Caspase-3 in vitro. On the basis of these structural studies of the cytoplasmic domains of pannexin 1, we propose a mechanism for the opening of pannexin 1 channels upon apoptosis, involving structural changes within the CT domain.     

A model for the nucleotide-binding domains of ABC transporters based on the large domain of aspartate aminotransferase

ABC transporters are a large superfamily of integral membrane proteins involved in ATP-dependent transport across biological membranes. Members of this superfamily play roles in a number of phenomena of biomedical interest, including cystic fibrosis (CFTR) and multidrug resistance (P-glycoprotein, MRP). Most ABC transporters are predicted to consist of four domains, two membrane-spanning domains and two cytoplasmic domains. The latter contain conserved nucleotide-binding motifs. Attempts to determine the structure of ABC transporters and of their separate domains are in progress but have not yet been successful. To aid structure determination and possibly learn more about the domain boundaries, we set out to model nucleotide-binding domains (NBDs) of ABC transporters based on a known structure. Previous attempts to predict the 3D structure of NBDs were based solely on sequence similarity with known nucleotide-binding folds. We have analyzed the sequences of a number of nucleotide-binding domains with the algorithm THREADER, developed by D.T. Jones, and a possible fold was found in the structure of aspartate aminotransferase. We present a model for the N-terminal NBD of CFTR, based on the large domain of the A chain of aspartate aminotransferase. The model is refined using multiple sequence alignment, secondary structure prediction, and 3D-1D profiles. Our model seems to be in good agreement with known properties of nucleotide-binding domains and has some appealing characteristics compared with the previous models. Proteins 30:275–286, 1998. © 1998 Wiley-Liss, Inc.     

Amino acid residues contributing to function of the heteromeric insect olfactory receptor complex

Olfactory receptors (Ors) convert chemical signals--the binding of odors and pheromones--to electrical signals through the depolarization of olfactory sensory neurons. Vertebrates Ors are G-protein-coupled receptors, stimulated by odors to produce intracellular second messengers that gate ion channels. Insect Ors are a heteromultimeric complex of unknown stoichiometry of two seven transmembrane domain proteins with no sequence similarity to and the opposite membrane topology of G-protein-coupled receptors. The functional insect Or comprises an odor- or pheromone-specific Or subunit and the Orco co-receptor, which is highly conserved in all insect species. The insect Or-Orco complex has been proposed to function as a novel type of ligand-gated nonselective cation channel possibly modulated by G-proteins. However, the Or-Orco proteins lack homology to any known family of ion channel and lack known functional domains. Therefore, the mechanisms by which odors activate the Or-Orco complex and how ions permeate this complex remain unknown. To begin to address the relationship between Or-Orco structure and function, we performed site-directed mutagenesis of all 83 conserved Glu, Asp, or Tyr residues in the silkmoth BmOr-1-Orco pheromone receptor complex and measured functional properties of mutant channels expressed in Xenopus oocytes. 13 of 83 mutations in BmOr-1 and BmOrco altered the reversal potential and rectification index of the BmOr-1-Orco complex. Three of the 13 amino acids (D299 and E356 in BmOr-1 and Y464 in BmOrco) altered both current-voltage relationships and K(+) selectivity. We introduced the homologous Orco Y464 residue into Drosophila Orco in vivo, and observed variable effects on spontaneous and evoked action potentials in olfactory neurons that depended on the particular Or-Orco complex examined. Our results provide evidence that a subset of conserved Glu, Asp and Tyr residues in both subunits are essential for channel activity of the heteromeric insect Or-Orco complex.     

Structural order in Pannexin 1 cytoplasmic domains

Pannexin 1 forms ion and metabolite permeable hexameric channels with abundant expression in the central nervous system and elsewhere. Although pannexin 1 does not form intercellular channels, a common channel topology and oligomerization state, as well as involvement of the intracellular carboxyl terminal (CT) domain in channel gating, is shared with connexins. In this study, we characterized the secondary structure of the mouse pannexin 1 cytoplasmic domains to complement structural studies of the transmembrane segments and compare with similar domains from connexins. A combination of structural prediction tools and circular dichroism revealed that, unlike connexins (predominately intrinsically disordered), cytosolic regions of pannexin 1 contain approximately 50% secondary structure, a majority being α-helical. Moreover, prediction of transmembrane domains uncovered a potential membrane interacting region (I360-G370) located upstream of the caspase cleavage site (D375-D378) within the pannexin 1 CT domain. The α-helical content of a peptide containing these domains (G357-S384) increased in the presence of detergent micelles providing evidence of membrane association. We also purified a pannexin 1 CT construct containing the caspase cleavage site (M374-C426), assigned the resonances by NMR, and confirmed cleavage by Caspase-3 in vitro. On the basis of these structural studies of the cytoplasmic domains of pannexin 1, we propose a mechanism for the opening of pannexin 1 channels upon apoptosis, involving structural changes within the CT domain.     

Potassium channel structure: domain by domain

Since the determination of the structure of a bacterial potassium channel, the ion channel community has managed to gain momentum in the quest for a complete picture. The information is coming at a steady flow, on a domain by domain basis. Recent discoveries are starting to reveal clues to the complex manner in which potassium channels show enormous diversity of function and also to their methods of regulation. Currently, the structures of four domains are known, with the most recent addition being the Kvbeta structure. As efforts continue in the study of the transmembrane domains, especially the voltage-sensing apparatus, there has been a new realization with respect to the identification and role of the cytoplasmic domains in protein-protein interactions in particular. An additional discovery, considerably aided by recent genomic analysis, is that potassium channels comprising subunits with two pore regions and four transmembrane helices combined in a dimeric fashion are abundant and are probable targets for local anesthetics.     

Identification of a novel tetramerization domain in large conductance K(ca) channels

More than 50 genes are known to encode K(+) channel monomers and can coassemble to form hetero-tetrameric K(+) channels. However, only a subset of possible monomer combinations come together to form functional ion channels. The assembly and tetramerization of appropriate channel monomers is mediated by association domains (ADs). To identify such domains in human large-conductance Ca(2+)-activated K(+) channels (hSlo1), we screened hSlo1 domains for self-association using yeast two-hybrid assays. Putative ADs were subjected to functional assays in Xenopus oocytes and further characterized by coprecipitation, native gel electrophoresis, and sucrose density gradient centrifugation assays. This led to the identification of a single intracellular association domain localized near the channel pore and required for channel function. We conclude that this novel tetramerization domain, referred to as BK-T1, promotes the assembly of hSlo1 monomers into functional K(Ca) channels.     

A localized interaction surface for voltage-sensing domains on the pore domain of a K+ channel

Voltage-gated K+ channels contain a central pore domain and four surrounding voltage-sensing domains. How and where changes in the structure of the voltage-sensing domains couple to the pore domain so as to gate ion conduction is not understood. The crystal structure of KcsA, a bacterial K+ channel homologous to the pore domain of voltage-gated K+ channels, provides a starting point for addressing this question. Guided by this structure, we used tryptophan-scanning mutagenesis on the transmembrane shell of the pore domain in the Shaker voltage-gated K+ channel to localize potential protein-protein and protein-lipid interfaces. Some mutants cause only minor changes in gating and when mapped onto the KcsA structure cluster away from the interface between pore domain subunits. In contrast, mutants producing large changes in gating tend to cluster near this interface. These results imply that voltage-sensing domains interact with localized regions near the interface between adjacent pore domain subunits.     

Base of the thumb domain modulates epithelial sodium channel gating

The activity of the epithelial sodium channel (ENaC) is modulated by multiple external factors, including proteases, cations, anions and shear stress. The resolved crystal structure of acid-sensing ion channel 1 (ASIC1), a structurally related ion channel, and mutagenesis studies suggest that the large extracellular region is involved in recognizing external signals that regulate channel gating. The thumb domain in the extracellular region of ASIC1 has a cylinder-like structure with a loop at its base that is in proximity to the tract connecting the extracellular region to the transmembrane domains. This loop has been proposed to have a role in transmitting proton-induced conformational changes within the extracellular region to the gate. We examined whether loops at the base of the thumb domains within ENaC subunits have a similar role in transmitting conformational changes induced by external Na(+) and shear stress. Mutations at selected sites within this loop in each of the subunits altered channel responses to both external Na(+) and shear stress. The most robust changes were observed at the site adjacent to a conserved Tyr residue. In the context of channels that have a low open probability due to retention of an inhibitory tract, mutations in the loop activated channels in a subunit-specific manner. Our data suggest that this loop has a role in modulating channel gating in response to external stimuli, and are consistent with the hypothesis that external signals trigger movements within the extracellular regions of ENaC subunits that are transmitted to the channel gate.     

Isoflurane alters the structure and dynamics of GLIC

Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.