Water and potassium dynamics inside the KcsA K(+) channel

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.     

Molecular dynamics of the KcsA K(+) channel in a bilayer membrane

Molecular dynamics (MD) simulations of an atomic model of the KcsA K(+) channel embedded in an explicit dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayer solvated by a 150 mM KCl aqueous salt solution are performed and analyzed. The model includes the KcsA K(+) channel, based on the recent crystallographic structure of, Science. 280:69-77), 112 DPPC, K(+) and Cl(-) ions, as well as over 6500 water molecules for a total of more than 40,000 atoms. Three K(+) ions are explicitly included in the pore. Two are positioned in the selectivity filter on the extracellular side and one in the large water-filled cavity. Different starting configurations of the ions and water molecules in the selectivity filter are considered, and MD trajectories are generated for more than 4 ns. The conformation of KcsA is very stable in all of the trajectories, with a global backbone root mean square (RMS) deviation of less than 1.9 A with respect to the crystallographic structure. The RMS atomic fluctuations of the residues surrounding the selectivity filter on the extracellular side of the channel are significantly lower than those on the intracellular side. The motion of the residues with aromatic side chains surrounding the selectivity filter (Trp(67), Trp(68), Tyr(78), and Tyr(82)) is anisotropic with the smallest RMS fluctuations in the direction parallel to the membrane plane. A concerted dynamic transition of the three K(+) ions in the pore is observed, during which the K(+) ion located initially in the cavity moves into the narrow part of the selectivity filter, while the other two K(+) ions move toward the extracellular side. A single water molecule is stabilized between each pair of ions during the transition, suggesting that each K(+) cation translocating through the narrow pore is accompanied by exactly one water molecule, in accord with streaming potential measurements (, Biophys. J. 55:367-371). The displacement of the ions is coupled with the structural fluctuations of Val(76) and Gly(77), in the selectivity filter, as well as the side chains of Glu(71), Asp(80), and Arg(89), near the extracellular side. Thus the mechanical response of the channel structure at distances as large as 10-20 A from the ions in the selectivity filter appears to play an important role in the concerted transition.     

The ionization state and the conformation of Glu-71 in the KcsA K(+) channel.

The side chain of Glu-71 of the KcsA K(+) channel, an important residue in the vicinity of the selectivity filter, was not resolved in the crystallographic structure of Doyle et al. (Doyle, D. A., J. M. Cabral, R. A. Pfuetzner, A. Kuo, J. M. Gulbis, S. L. Cohen, B. T. Chait, and R. MacKinnon. 1998. Science. 280:69-77). Its atomic coordinates are undetermined and its ionization state is unknown. For meaningful theoretical and computational studies of the KcsA K(+) channel, it is essential to address questions about the conformation and the ionization state of this residue in detail. In previous MD simulations in which the side chain of Glu-71 is protonated and forming a strong hydrogen bond with Asp-80 it was observed that the channel did not deviate significantly from the crystallographic structure (Bernèche, S., and B. Roux. 2000. Biophys. J. 78:2900-2917). In contrast, we show here that the structure of the selectivity filter of the KcsA channel is significantly disrupted when these side chains are fully ionized on each of the four monomers. To further resolve questions about the ionization state of Glu-71 we calculated the pK(a) value of this residue using molecular dynamics free energy simulations (MD/FES) with a fully flexible system including explicit solvent and membrane and finite-difference Poisson-Boltzmann (PB) continuum electrostatics. It is found that the pK(a) of Glu-71 is shifted by approximately +10 pK(a) units. These results strongly suggest that Glu-71 is protonated under normal conditions.     

Extracellular blockade of K(+) channels by TEA: results from molecular dynamics simulations of the KcsA channel

TEA is a classical blocker of K(+) channels. From mutagenesis studies, it has been shown that external blockade by TEA is strongly dependent upon the presence of aromatic residue at Shaker position 449 which is located near the extracellular entrance to the pore (Heginbotham, L., and R. MacKinnon. 1992. Neuron. 8:483-491). The data suggest that TEA interacts simultaneously with the aromatic residues of the four monomers. The determination of the 3-D structure of the KcsA channel using X-ray crystallography (Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. Science. 280:69-77) has raised some issues that remain currently unresolved concerning the interpretation of these observations. In particular, the center of the Tyr82 side chains in KcsA (corresponding to position 449 in Shaker) forms a square of 11.8-A side, a distance which is too large to allow simultaneous interactions of a TEA molecule with the four aromatic side chains. In this paper, the external blockade by TEA is explored by molecular dynamics simulations of an atomic model of KcsA in an explicit phospholipid bilayer with aqueous salt solution. It is observed, in qualitative accord with the experimental results, that TEA is stable when bound to the external side of the wild-type KcsA channel (with Tyr82), but is unstable when bound to a mutant channel in which the tyrosine residue has been substituted by a threonine. The free energy profile of TEA relative to the pore is calculated using umbrella sampling simulations to characterize quantitatively the extracellular blockade. It is found, in remarkable agreement with the experiment, that the TEA is more stably bound by 2.3 kcal/mol to the channel with four tyrosine residues. In the case of the wild-type KcsA channel, TEA (which has the shape of a flattened oblate spheroid) acts as an ideal plug blocking the pore. In contrast, it is considerably more off-centered and tilted in the case of the mutant channel. The enhanced stability conferred by the tyrosine residues does not arise from Pi-cation interactions, but appears to be due to differences in the hydration structure of the TEA. Finally, it is shown that the experimentally observed voltage dependence of TEA block, which is traditionally interpreted in terms of the physical position of the TEA along the axis of the pore, must arise indirectly via coupling with the ions in the pore.     

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.     

A mutant KcsA K(+) channel with altered conduction properties and selectivity filter ion distribution

The selectivity filter of K(+) channels is comprised of a linear queue of four equal-spaced ion-binding sites spanning a distance of 12A. Each site is formed of eight oxygen atoms from the protein. The first three sites, numbered 1-3 from the extracellular side, are made of exclusively main-chain carbonyl oxygen atoms. The fourth site, closest to the intracellular side, is made of four main-chain carbonyl oxygen atoms and four threonine side-chain hydroxyl oxygen atoms. Here we characterize the effects of mutating the threonine to cysteine on the distribution of ions in the selectivity filter and on the conduction of ions through the filter. The mutation influences the occupancy of K(+) at sites 2 and 4 and it reduces the maximum rate of conduction in the limit of high K(+) concentration. The mutation does not affect the conduction of Rb(+). These results can be understood in the context of a conduction mechanism in which a pair of K ions switch between energetically balanced 1,3 and 2,4 configurations.     

K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations

Microscopic molecular dynamics free energy perturbation calculations of the K(+)/Na(+) selectivity in the KcsA potassium channel, based on its experimental three-dimensional structure, are reported. The relative binding free energies for K(+) and Na(+) in the most relevant ion occupancy states of the four-site selectivity filter are calculated. The previously proposed mechanism for ion permeation through the KcsA channel is predicted, in agreement with available experimental data, to have a significant selectivity for K(+) over Na(+). The calculations also show that the individual 'binding site' selectivities are generally not additive and the doubly loaded states of the filter thus display cooperative effects. The only site that is not K(+) selective is that which is located at the entrance to the internal water cavity, suggesting the possibility that internal Na(+) could block outward currents.     

Permeation of water through the KcsA K+ channel

Previous studies have reported that the KcsA potassium channel has an osmotic permeability coefficient of 4.8 x 10(-12) cm3/s, giving it a significantly higher osmotic permeability coefficient than that of some membrane channels specialized in water transport. This high osmotic permeability is proposed to occur when the channel is depleted of potassium ions, the presence of which slow down the water permeation process. The atomic structure of the potassium-depleted KcsA channel and the mechanisms of water permeation have not been well characterized so far. Here, all-atom molecular dynamics simulations, in conjunction with an umbrella sampling strategy and a nonequilibrium approach to simulate pressure gradients are employed to illustrate the permeation of water in the absence of ions through the KcsA K+ channel. Equilibrium molecular dynamics simulations (95 ns combined total length) identified a possible structure of the potassium-depleted KcsA channel, and umbrella sampling calculations (160 ns combined total length) revealed that this structure is not permeable by water molecules moving along the channel axis. The simulation of a pressure gradient across the channel (30 ns combined total length) identified an alternative permeation pathway with a computed osmotic permeability of approximately (2.7 +/- 0.9) x 10(-13) cm3/s. Water fluxes along this pathway did not proceed through collective water motions or transitions to vapor state. All of the major results of this study were robust against variations in a wide set of simulation parameters (force field, water model, membrane model, and channel conformation).     

Physical and functional interaction of the Arabidopsis K(+) channel AKT2 and phosphatase AtPP2CA

The AKT2 K(+) channel is endowed with unique functional properties, being the only weak inward rectifier characterized to date in Arabidopsis. The gene is expressed widely, mainly in the phloem but also at lower levels in leaf epiderm, mesophyll, and guard cells. The AKT2 mRNA level is upregulated by abscisic acid. By screening a two-hybrid cDNA library, we isolated a protein phosphatase 2C (AtPP2CA) involved in abscisic acid signaling as a putative partner of AKT2. We further confirmed the interaction by in vitro binding studies. The expression of AtPP2CA (beta-glucuronidase reporter gene) displayed a pattern largely overlapping that of AKT2 and was upregulated by abscisic acid. Coexpression of AtPP2CA with AKT2 in COS cells and Xenopus laevis oocytes was found to induce both an inhibition of the AKT2 current and an increase of the channel inward rectification. Site-directed mutagenesis and pharmacological analysis revealed that this functional interaction involves AtPP2CA phosphatase activity. Regulation of AKT2 activity by AtPP2CA in planta could allow the control of K(+) transport and membrane polarization during stress situations.     

Structural compatibility between the putative voltage sensor of voltage-gated K+ channels and the prokaryotic KcsA channel

Sequence similarity among and electrophysiological studies of known potassium channels, along with the three-dimensional structure of the Streptomyces lividans K(+) channel (KcsA), support the tenet that voltage-gated K(+) channels (Kv channels) consist of two distinct modules: the "voltage sensor" module comprising the N-terminal portion of the channel up to and including the S4 transmembrane segment and the "pore" module encompassing the C-terminal portion from the S5 transmembrane segment onward. To substantiate this modular design, we investigated whether the pore module of Kv channels may be replaced with the pore module of the prokaryotic KcsA channel. Biochemical and immunocytochemical studies showed that chimeric channels were expressed on the cell surface of Xenopus oocytes, demonstrating that they were properly synthesized, glycosylated, folded, assembled, and delivered to the plasma membrane. Unexpectedly, surface-expressed homomeric chimeras did not exhibit detectable voltage-dependent channel activity upon both hyperpolarization and depolarization regardless of the expression system used. Chimeras were, however, strongly dominant-negative when coexpressed with wild-type Kv channels, as evidenced by the complete suppression of wild-type channel activity. Notably, the dominant-negative phenotype correlated well with the formation of stable, glycosylated, nonfunctional, heteromeric channels. Collectively, these findings imply a structural compatibility between the prokaryotic pore module and the eukaryotic voltage sensor domain that leads to the biogenesis of non-responsive channels. Our results lend support to the notion that voltage-dependent channel gating depends on the precise coupling between both protein domains, probably through a localized interaction surface.     

Determination of the interface of a large protein complex by transferred cross-saturation measurements

In an earlier paper, it was shown that the cross-saturation method enables us to identify the contact residues of large protein complexes in a more rigorous manner than is possible using chemical shift perturbation and hydrogen-deuterium exchange experiments. However, there are limitations within the determination of the contact residues by the cross-saturation method, in that the method is difficult to apply to protein complexes with a molecular mass over 150 kDa and/or with weak binding, since the resonances originating from the complexes should be observed directly in the method. In the present work, to overcome these limitations, we carried out the cross-saturation measurements under conditions of a fast exchange between free and bound states on the NMR time-scale, and determined the contact residues of the complex of the B domain of protein A and intact IgG, which has a molecular mass of 164 kDa and shows weak binding.     

Structural basis for the interaction of Bordetella pertussis adenylyl cyclase toxin with calmodulin

CyaA is crucial for colonization by Bordetella pertussis, the etiologic agent of whooping cough. Here we report crystal structures of the adenylyl cyclase domain (ACD) of CyaA with the C-terminal domain of calmodulin. Four discrete regions of CyaA bind calcium-loaded calmodulin with a large buried contact surface. Of those, a tryptophan residue (W242) at an alpha-helix of CyaA makes extensive contacts with the calcium-induced, hydrophobic pocket of calmodulin. Mutagenic analyses show that all four regions of CyaA contribute to calmodulin binding and the calmodulin-induced conformational change of CyaA is crucial for catalytic activation. A crystal structure of CyaA-calmodulin with adefovir diphosphate, the metabolite of an approved antiviral drug, reveals the location of catalytic site of CyaA and how adefovir diphosphate tightly binds CyaA. The ACD of CyaA shares a similar structure and mechanism of activation with anthrax edema factor (EF). However, the interactions of CyaA with calmodulin completely diverge from those of EF. This provides molecular details of how two structurally homologous bacterial toxins evolved divergently to bind calmodulin, an evolutionarily conserved calcium sensor.     

Structural basis of Na(+)/K(+)-ATPase adaptation to marine environments

Throughout evolution, enzymes have adapted to perform in different environments. The Na(+)/K(+) pump, an enzyme crucial for maintaining ionic gradients across cell membranes, is strongly influenced by the ionic environment. In vertebrates, the pump sees much less external Na(+) (100-160 mM) than it does in osmoconformers such as squid (450 mM), which live in seawater. If the extracellular architecture of the squid pump were identical to that of vertebrates, then at the resting potential, the pump's function would be severely compromised because the negative voltage would drive Na(+) ions back to their binding sites, practically abolishing forward transport. Here we show that four amino acids that ring the external mouth of the ion translocation pathway are more positive in squid, thereby reducing the pump's sensitivity to external Na(+) and explaining how it can perform optimally in the marine environment.     

Crystallographic study of the tetrabutylammonium block to the KcsA K+ channel

K(+) channels play essential roles in regulating membrane excitability of many diverse cell types by selectively conducting K(+) ions through their pores. Many diverse molecules can plug the pore and modulate the K(+) current. Quaternary ammonium (QA) ions are a class of pore blockers that have been used for decades by biophysicists to probe the pore, leading to important insights into the structure-function relation of K(+) channels. However, many key aspects of the QA-blocking mechanisms remain unclear to date, and understanding these questions requires high resolution structural information. Here, we address the question of whether intracellular QA blockade causes conformational changes of the K(+) channel selectivity filter. We have solved the structures of the KcsA K(+) channel in complex with tetrabutylammonium (TBA) and tetrabutylantimony (TBSb) under various ionic conditions. Our results demonstrate that binding of TBA or TBSb causes no significant change in the KcsA structure at high concentrations of permeant ions. We did observe the expected conformational change of the filter at low concentration of K(+), but this change appears to be independent of TBA or TBSb blockade.     

Structural basis of pore formation by the bacterial toxin pneumolysin

The bacterial toxin pneumolysin is released as a soluble monomer that kills target cells by assembling into large oligomeric rings and forming pores in cholesterol-containing membranes. Using cryo-EM and image processing, we have determined the structures of membrane-surface bound (prepore) and inserted-pore oligomer forms, providing a direct observation of the conformational transition into the pore form of a cholesterol-dependent cytolysin. In the pore structure, the domains of the monomer separate and double over into an arch, forming a wall sealing the bilayer around the pore. This transformation is accomplished by substantial refolding of two of the four protein domains along with deformation of the membrane. Extension of protein density into the bilayer supports earlier predictions that the protein inserts beta hairpins into the membrane. With an oligomer size of up to 44 subunits in the pore, this assembly creates a transmembrane channel 260 A in diameter lined by 176 beta strands.     

Molecular basis of the voltage-dependent gating of TREK-1, a mechano-sensitive K(+) channel

TREK-1 is a member of the mammalian two P domain K(+) channel family. Mouse TREK-1 activity, in transiently transfected COS cells, is reduced at negative resting membrane potentials by both an external Mg(2+) block and an intrinsic voltage-dependent gating mechanism leading to a strong outward rectification. Deletional and chimeric analysis demonstrates that the carboxy terminal domain of TREK-1, but not the PKA phosphorylation site S333, is responsible for voltage-dependent gating. Since the same region is also critically required for TREK-1 mechano-gating, both mechanisms might be functionally linked. Preferential opening of TREK-1 at depolarized potentials will greatly affect action potential duration, recovery from inactivation and neuronal repetitive firing activity.     

Solution structure of the potassium channel inhibitor agitoxin 2: caliper for probing channel geometry.

The structure of the potassium channel blocker agitoxin 2 was solved by solution NMR methods. The structure consists of a triple-stranded antiparallel beta-sheet and a single helix covering one face of the beta-sheet. The cysteine side chains connecting the beta-sheet and the helix form the core of the molecule. One edge of the beta-sheet and the adjacent face of the helix form the interface with the Shaker K+ channel. The fold of agitoxin is homologous to the previously determined folds of scorpion venom toxins. However, agitoxin 2 differs significantly from the other channel blockers in the specificity of its interactions. This study was thus focused on a precise characterization of the surface residues at the face of the protein interacting with the Shaker K+ channel. The rigid toxin molecule can be used to estimate dimensions of the potassium channel. Surface-exposed residues, Arg24, Lys27, and Arg31 of the beta-sheet, have been identified from mutagenesis studies as functionally important for blocking the Shaker K+ channel. The sequential and spatial locations of Arg24 and Arg31 are not conserved among the homologous toxins. Knowledge on the details of the channel-binding sites of agitoxin 2 formed a basis for site-directed mutagenesis studies of the toxin and the K+ channel sequences. Observed interactions between mutated toxin and channel are being used to elucidate the channel structure and mechanisms of channel-toxin interactions.     

Functional influence of the pore helix glutamate in the KcsA K+ channel

The E71 residue is buried near the selectivity filter in the KcsA K+ channel and forms a carboxyl-carboxylate bridge with D80. We have investigated the importance of E71 by examining neutralization mutants at this position using biochemical and electrophysiological methods. E71 mutations differentially destabilize the detergent-solubilized tetramer; among them, the E71V neutralization mutant has a relatively subtle effect. The E71V channel displays electrical activity when reconstituted into planar lipid bilayers. In single channel recordings, the mutant retains K+/Na+ selectivity, and its conductance in the outward direction is unaltered. Some conduction properties are changed: inward conductance is increased. Our results show that that the E71 side chain is not a primary determinant of ion selectivity or conduction in the wild-type channel, either directly or through the E71:D80 carboxyl-carboxylate bridge.