The bacterial K+ channel structure and its implications for neuronal channels

Is Smaller Better? A Proposal to Use Bacteria For Neuroscientific Modeling

Bacteria are easily characterizable model organisms with an impressively complicated set of abilities. Among them is quorum sensing, a cell-cell signaling system that may have a common evolutionary origin with eukaryotic cell-cell signaling. The two systems are behaviorally similar, but quorum sensing in bacteria is more easily studied in depth than cell-cell signaling in eukaryotes. Because of this comparative ease of study, bacterial dynamics are also more suited to direct interpretation than eukaryotic dynamics, e.g., those of the neuron. Here we review literature on neuron-like qualities of bacterial colonies and biofilms, including ion-based and hormonal signaling, and a phenomenon similar to the graded action potential. This suggests that bacteria could be used to help create more accurate and detailed biological models in neuroscientific research. More speculatively, bacterial systems may be considered an analog for neurons in biologically based computational research, allowing models to better harness the tremendous ability of biological organisms to process information and make decisions.

A tyrosine substitution in the cavity wall of a k channel induces an inverted inactivation

Ion permeation and gating kinetics of voltage-gated K channels critically depend on the amino-acid composition of the cavity wall. Residue 470 in the Shaker K channel is an isoleucine, making the cavity volume in a closed channel insufficiently large for a hydrated K(+) ion. In the cardiac human ether-a-go-go-related gene channel, which exhibits slow activation and fast inactivation, the corresponding residue is tyrosine. To explore the role of a tyrosine at this position in the Shaker channel, we studied I470Y. The activation became slower, and the inactivation faster and more complex. At +60 mV the channel inactivated with two distinct rates (tau(1) = 20 ms, tau(2) = 400 ms). Experiments with tetraethylammonium and high K(+) concentrations suggest that the slower component was of the P/C-type. In addition, an inactivation component with inverted voltage dependence was introduced. A step to -40 mV inactivates the channel with a time constant of 500 ms. Negative voltage steps do not cause the channel to recover from this inactivated state (tau >> 10 min), whereas positive voltage steps quickly do (tau = 2 ms at +60 mV). The experimental findings can be explained by a simple branched kinetic model with two inactivation pathways from the open state.

Molecular dynamics simulations of potassium channels

Despite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.

Local sequence information-based support vector machine to classify voltage-gated potassium channels

In our previous work, we developed a computational tool, PreK-ClassK-ClassKv, to predict and classify potassium (K+) channels. For K+ channel prediction (PreK) and classification at family level (ClassK), this method performs well. However, it does not perform so well in classifying voltage-gated potassium (Kv) channels (ClassKv). In this paper, a new method based on the local sequence information of Kv channels is introduced to classify Kv channels. Six transmembrane domains of a Kv channel protein are used to define a protein, and the dipeptide composition technique is used to transform an amino acid sequence to a numerical sequence. A Kv channel protein is represented by a vector with 2000 elements, and a support vector machine algorithm is applied to classify Kv channels. This method shows good performance with averages of total accuracy (Acc), sensitivity (SE), specificity (SP), reliability (R) and Matthews correlation coefficient (MCC) of 98.0%, 89.9%, 100%, 0.95 and 0.94 respectively. The results indicate that the local sequence information-based method is better than the global sequence information-based method to classify Kv channels.

Long Distance Interactions within the Potassium Channel Pore Are Revealed by Molecular Diversity of Viral Proteins

Kcv is a 94-amino acid protein encoded by chlorella virus PBCV-1 that corresponds to the pore module of K(+) channels. Therefore, Kcv can be a model for studying the protein design of K(+) channel pores. We analyzed the molecular diversity generated by approximately 1 billion years of evolution on kcv genes isolated from 40 additional chlorella viruses. Because the channel is apparently required for virus replication, the Kcv variants are all functional and contain multiple and dispersed substitutions that represent a repertoire of allowed sets of amino acid substitutions (from 4 to 12 amino acids). Correlations between amino acid substitutions and the new properties displayed by these channels guided site-directed mutations that revealed synergistic amino acid interactions within the protein as well as previously unknown interactions between distant channel domains. The effects of these multiple changes were not predictable from a priori structural knowledge of the channel pore.

Localization of voltage-gated ion channels in mammalian brain

The intrinsic electrical properties of neurons are shaped in large part by the action of voltage-gated ion channels. Molecular cloning studies have revealed a large family of ion channel genes, many of which are expressed in mammalian brain. Much recent effort has focused on determining the contribution of the protein products of these genes to neuronal function. This requires knowledge of the abundance and distribution of the constituent subunits of the channels in specific mammalian central neurons. Here we review progress made in recent studies aimed at localizing specific ion channel subunits using in situ hybridization and immunohistochemistry. We then discuss the implications of these results in terms of neuronal physiology and neuronal mechanisms underlying the observed patterns of expression.

Potassium channels: structures, models, simulations

Potassium channels have been studied intensively in terms of the relationship between molecular structure and physiological function. They provide an opportunity to integrate structural and computational studies in order to arrive at an atomic resolution description of mechanism. We review recent progress in K channel structural studies, focussing on the bacterial channel KcsA. Structural studies can be extended via use of computational (i.e. molecular simulation) approaches in order to provide a perspective on aspects of channel function such as permeation, selectivity, block and gating. Results from molecular dynamics simulations are shown to be in good agreement with recent structural studies of KcsA in terms of the interactions of K(+) ions with binding sites within the selectivity filter of the channel, and in revealing the importance of filter flexibility in channel function. We discuss how the KcsA structure may be used as a template for developing structural models of other families of K channels. Progress in this area is explored via two examples: inward rectifier (Kir) and voltage-gated (Kv) potassium channels. A brief account of structural studies of ancillary domains and subunits of K channels is provided.

Inactivation determinants in segment IIIS6 of Ca(v)3.1

1. Low threshold, T-type, Ca(2+) channels of the Ca(v)3 family display the fastest inactivation kinetics among all voltage-gated Ca(2+) channels. The molecular inactivation determinants of this channel family are largely unknown. Here we investigate whether segment IIIS6 plays a role in Ca(v)3.1 inactivation as observed previously in high voltage-activated Ca(2+) channels. 2. Amino acids that are identical in IIIS6 segments of all Ca(2+) channel subtypes were mutated to alanine (F1505A, F1506A, N1509A, F1511A, V1512A, F1519A, FV1511/1512AA). Additionally M1510 was mutated to isoleucine and alanine. 3. The kinetic properties of the mutants were analysed with the two-microelectrode voltage-clamp technique after expression in Xenopus oocytes. The time constant for the barium current (I(Ba)) inactivation, tau(inact), of wild-type channels at -20 mV was 9.5 +/- 0.4 ms; the corresponding time constants of the mutants ranged from 9.2 +/- 0.4 ms in V1512A to 45.7 +/- 5.2 ms (4.8-fold slowing) in M1510I. Recovery at -80 mV was most significantly slowed by V1512A and accelerated by F1511A. 4. We conclude that amino acids M1510, F1511 and V1512 corresponding to previously identified inactivation determinants in IIIS6 of Ca(v)2.1 (Hering et al. 1998) have a significant role in Ca(v)3.1 inactivation. These data suggest common elements in the molecular architecture of the inactivation mechanism in high and low threshold Ca(2+) channels.

Controlling potassium channel activities: Interplay between the membrane and intracellular factors

Neural signaling is based on the regulated timing and extent of channel opening; therefore, it is important to understand how ion channels open and close in response to neurotransmitters and intracellular messengers. Here, we examine this question for potassium channels, an extraordinarily diverse group of ion channels. Voltage-gated potassium (Kv) channels control action-potential waveforms and neuronal firing patterns by opening and closing in response to membrane-potential changes. These effects can be strongly modulated by cytoplasmic factors such as kinases, phosphatases, and small GTPases. A Kv alpha subunit contains six transmembrane segments, including an intrinsic voltage sensor. In contrast, inwardly rectifying potassium (Kir) channels have just two transmembrane segments in each of its four pore-lining alpha subunits. A variety of intracellular second messengers mediate transmitter and metabolic regulation of Kir channels. For example, Kir3 (GIRK) channels open on binding to the G protein betagamma subunits, thereby mediating slow inhibitory postsynaptic potentials in the brain. Our structure-based functional analysis on the cytoplasmic N-terminal tetramerization domain T1 of the voltage-gated channel, Kv1.2, uncovered a new function for this domain, modulation of voltage gating, and suggested a possible means of communication between second messenger pathways and Kv channels. A yeast screen for active Kir3.2 channels subjected to random mutagenesis has identified residues in the transmembrane segments that are crucial for controlling the opening of Kir3.2 channels. The identification of structural elements involved in potassium channel gating in these systems highlights principles that may be important in the regulation of other types of channels.

Homology model of dihydropyridine receptor: implications for L-type Ca(2+) channel modulation by agonists and antagonists

L-type calcium channels (LCCs) are transmembrane (TM) proteins that respond to membrane depolarization by selectively permeating Ca(2+) ions. Dihydropyridine (DHP) agonists and antagonist modulate Ca(2+) permeation by stabilizing, respectively, the open and closed states of the channel. The mechanism of action of these drugs remains unclear. Using, as a template, the crystal structure of the KcsA K(+) channel (Doyle et al. (1998) Science 280, 69-77), we have built several homology models of LCC with alternative alignments of TM segments between the proteins. In each model, nifedipine was docked in the pore region and in the interface between repeats III and IV. Several starting structures were generated by constraining the ligand to residues whose mutations reportedly affect DHP binding (DHP-sensing residues). These structures were Monte Carlo-minimized with and without constraints. In the complex with the maximum number of contacts between the ligand and DHP-sensing residues and the lowest ligand-receptor energy, the drug fits snugly in the "water-lake" cavity between segments S6s, which were aligned with M2 segment of KcsA as proposed for Na(+) channel (Lipkind and Fozzard (2000) Biochemistry 39, 8161-8170). In the flattened-boat conformation of DHP ring, the NH group at the stern approaches the DHP-sensing tyrosines in segments IIIS6 and IVS6. Stacking interactions of IVS6 Tyr with the bowsprit aromatic ring stabilize the ligand's orientation in which the starboard COOMe group coordinates Ca(2+) ion chelated by two conserved glutamates in the selectivity filter. In the inverted teepee structure of LCC, the portside COOMe group approaches a bracelet of conserved hydrophobic residues at the helical-bundle crossing, which may function as the activation gate. The dimensions of the gate may readily change upon small rotation of the pore-forming TM segments. The end of the portside group is hydrophobic in nifedipine, (R)-Bay K 8644, and other antagonists. Favorable interactions of this group with the hydrophobic bracelet would stabilize its closed conformation. In contrast, (S)-Bay K 8644 and several other agonists have hydrophilic groups at the portside. Unfavorable interactions of the hydrophilic group with the hydrophobic bracelet would destabilize its closed conformation thereby stabilizing the open conformation. In the agonist-bound channel, Ca(2+) ions would permeate between the hydrophilic face of the ligand and conserved hydrophilic residues in segments IS6 and IIS6. Our model suggests mutational experiments that could further our understanding of the pharmacological modulation of voltage-gated ion channels.

Amino Acids in Segment IVS6 and β-Subunit Interaction Support Distinct Conformational Changes during Cav2.1 Inactivation

Ca(v)2.1 mediates voltage-gated Ca2+ entry into neurons and the release of neurotransmitters at synapses of the central nervous system. An inactivation process that is modulated by the auxiliary beta-subunits regulates Ca2+ entry through Ca(v)2.1. However, the molecular mechanism of this alpha1-beta-subunit interaction remains unknown. Herein we report the identification of new determinants within segment IVS6 of the alpha(1)2.1-subunit that markedly influence channel inactivation. Systematic substitution of residues within IVS6 with amino acids of different size, charge, and polarity resulted in mutant channels with rates of fast inactivation (k(inact)) ranging from a 1.5-fold slowing in V1818I (k(inact) = 0.98 +/- 0.09 s(-1) compared with wild type alpha(1)2.1/alpha2-delta/beta1a k(inact) = 1.35 +/- 0.25 s(-1) to a 75-fold acceleration in mutant M1811Q (k(inact) = 102 +/- 3 s(-1). Coexpression of mutant alpha(1)2.1-subunits with beta(2a) resulted in two different phenotypes of current inactivation: 1) a pronounced reduction in the rate of channel inactivation or 2) an attenuation of a slow component in I(Ba) inactivation. Simulations revealed that these two distinct inactivation phenotypes arise from a beta2a-subunit-induced destabilization of the fast-inactivated state. The IVS6- and beta2a-subunit-mediated effects on Ca(v)2.1 inactivation are likely to occur via independent mechanisms.

Structural similarities between glutamate receptor channels and K(+) channels examined by scanning mutagenesis

The pores of glutamate receptors and K(+) channels share sequence homology, suggesting a conserved secondary structure. Scanning mutagenesis with substitution of alanine and tryptophan in GluR6 channels was performed based on the structure of KcsA. Our assay used disruption of voltage-dependent polyamine block to test for changes in the packing of pore-forming regions. Alanine scanning from D567 to R603 revealed reduced rectification resulting from channel block in two regions. A periodic pattern from F575 to M589 aligned with the pore helix in KcsA, whereas a cluster of sensitive positions around Q590, a site regulated by RNA editing, mapped to the selectivity filter in KcsA. Tryptophan scanning from D567 to R603 revealed similar patterns, but with a complete disruption of spermine block for 7 out of the 37 positions and a pM dissociation constant for Q590W. Molecular modeling with KcsA coordinates showed that GluR6 pore helix mutants disrupting polyamine block pack against M1 and M2, and are not exposed in the ion channel pore. In the selectivity filter, tryptophan creates an aromatic cage consistent with the pM dissociation constant for Q590W. A scan with glutamate substitution was used to map the cytoplasmic entrance to the pore based on charge neutralization experiments, which established that E594 was uniquely required for high affinity polyamine block. In E594Q mutants, introduction of glutamate at positions S593-L600 restored polyamine block at positions corresponding to surface-exposed residues in KcsA. Our results reinforce proposals that the pore region of glutamate receptors contains a helix and pore loop analogous to that found in K(+) channels. At the cytoplasmic entrance of the channel, a negatively charged amino acid, located in an extended loop with solvent-exposed side chains, is required for high affinity polyamine block and probably attracts cations via a through space electrostatic mechanism.

Side-chain ionization states in a potassium channel

KcsA is a bacterial K+ channel that is gated by pH. Continuum dielectric calculations on the crystal structure of the channel protein embedded in a low dielectric slab suggest that side chains E71 and D80 of each subunit, which lie adjacent to the selectivity filter region of the channel, form a proton-sharing pair in which E71 is neutral (protonated) and D80 is negatively charged at pH 7. When K+ ions are introduced into the system at their crystallographic positions the pattern of proton sharing is altered. The largest perturbation is for a K+ ion at site S3, i.e., interacting with the carbonyls of T75 and V76. The presence of multiple K+ ions in the filter increases the probability of E71 being ionized and of D80 remaining neutral (i.e., protonated). The ionization states of the protein side chains influence the potential energy profile experienced by a K+ ion as it is translated along the pore axis. In particular, the ionization state of the E71-D80 proton-sharing pair modulates the shape of the potential profile in the vicinity of the selectivity filter. Such reciprocal effects of ion occupancy on side-chain ionization states, and of side-chain ionization states on ion potential energy profiles will complicate molecular dynamics simulations and related studies designed to calculate ion permeation energetics.

Signalling mechanisms: a decade of signalling

Knowledge of signaling mechanisms has increased dramatically during the past decade, particularly in the areas of development, biochemical signaling cascades, synaptic transmission and ion channel biophysics.

Effects of organic antagonists of Ca(2+), Na(+), and K(+) on chemotaxis and motility of escherichia coli

Various Ca(2+) antagonists used in animal research, many of them known to be Ca(2+) channel blockers, inhibited Escherichia coli chemotaxis (measured as entry of cells into a capillary containing attractant). The most effective of these, acting in the nanomolar range, was omega-conotoxin GVIA. The next most effective were gallopamil and verapamil. At concentrations around 100-fold higher than that needed for inhibition of chemotaxis, each of these antagonists inhibited motility (measured as entry of cells into a capillary lacking attractant). Various other Ca(2+) antagonists were less effective, though chemotaxis was almost always more sensitive to inhibition than was motility. Cells treated with each of these Ca(2+) antagonists swam with a running bias, i.e., tumbling was inhibited. Similarly, some Na(+) antagonists used in animal research inhibited bacterial chemotaxis. E. coli chemotaxis was inhibited by saxitoxin at concentrations above 10(-7) M, while more than 10(-4) M was needed to inhibit motility. Cells treated with saxitoxin swam with a tumbling bias. In the case of other Na(+) antagonists in animals, aconitine inhibited bacterial chemotaxis 10 times more effectively than it inhibited motility, and two others inhibited chemotaxis and motility at about the same concentration. In the case of K(+) antagonists used in animal research, 4-aminopyridine blocked E. coli chemotaxis between 10(-3) M and, totally, 10(-2) M, while motility was not affected at 10(-2) M; on the other hand, tetraethylammonium chloride failed to inhibit either chemotaxis or motility at 10(-2) M.

A potassium channel protein encoded by chlorella virus PBCV-1

The large chlorella virus PBCV-1, which contains double-stranded DNA (dsDNA), encodes a 94-codon open reading frame (ORF) that contains a motif resembling the signature sequence of the pore domain of potassium channel proteins. Phylogenetic analyses of the encoded protein, Kcv, indicate a previously unidentified type of potassium channel. The messenger RNA encoded by the ORF leads to functional expression of a potassium-selective conductance in Xenopus laevis oocytes. The channel blockers amantadine and barium, but not cesium, inhibit this conductance, in addition to virus plaque formation. Thus, PBCV-1 encodes the first known viral protein that functions as a potassium-selective channel and is essential in the virus life cycle.

Stable cation coordination at a single outer pore residue defines permeation properties in Kir channels

In epithelial Kir7.1 channels a non-conserved methionine in the outer pore region adjacent to the G-Y-G selectivity filter (position +2) was found to determine unique properties for permeant and blocking ions characteristic of a K(+) channel in a single-occupancy state. The monovalent cation permeability sequence of Kir7.1 channels expressed in Xenopus oocytes was Tl(+)>K(+)>Rb(+)NH(4)(+)>Cs(+)>Na(+)>Li(+), but the macroscopic conductance for Rb(+) was approximately 8-fold larger than for the smaller K(+) ions, and decreased approximately 40-fold with the conserved arginine at the +2 position (Kir7.1M125R). Moreover, in Kir7.1 Rb(+) restored the typical permeation properties of other multi-ion channels indicating that a stable coordination of permeant ions at the +2 position defines the initial step in the conduction pathway of Kir channels.

Functional characterization of a potassium-selective prokaryotic glutamate receptor

Ion channels are molecular pores that facilitate the passage of ions across cell membranes and participate in a range of biological processes, from excitatory signal transmission in the mammalian nervous system to the modulation of swimming behaviour in the protozoan Paramecium. Two particularly important families of ion channels are ionotropic glutamate receptors (GluRs) and potassium channels. GluRs are permeable to Na+, K+ and Ca2+, are gated by glutamate, and have previously been found only in eukaryotes. In contrast, potassium channels are selective for K+, are gated by a range of stimuli, and are found in both prokaryotes and eukaryotes. Here we report the discovery and functional characterization of GluR0 from Synechocystis PCC 6803, which is the first GluR found in a prokaryote. GluR0 binds glutamate, forms potassium-selective channels and is related in amino-acid sequence to both eukaryotic GluRs and potassium channels. On the basis of amino-acid sequence and functional relationships between GluR0 and eukaryotic GluRs, we propose that a prokaryotic GluR was the precursor to eukaryotic GluRs. GluR0 provides evidence for the missing link between potassium channels and GluRs, and we suggest that their ion channels have a similar architecture, that GluRs are tetramers and that the gating mechanisms of GluRs and potassium channels have some essential features in common.