Ion conductance vs. pore gating and selectivity in KcsA channel: Modeling achievements and perspectives

KcsA potassium channel belongs to a wide family of allosteric proteins that switch between closed and open states conformations in response to a stimulus, and act as a regulator of cation activity in living cells. The gating mechanism and cation selectivity of such channels have been extensively studied in the literature, with a revival emphasis these latter years, due to the publication of the crystallized structure of KcsA. Despite the increasing number of research and review papers on these topics, quantitative interpretation of these processes at the atomic scale is far from achieved. On the basis of available experimental and theoretical data, and by including our recent results, we review the progresses in this field of activity and discuss the weaknesses that should be corrected. In this spirit, we partition the channel into the filter, cavity, extra and intracellular media, in order to analyze separately the specificity of each region. Special emphasis is brought to the study of an open state for the channel and to the different properties generated by the opening. The influence of water as a structural and dynamical component of the channel properties in closed and open states, as well as in the sequential motions of the cations, is analyzed using molecular dynamics simulations and ab initio calculations. The polarization and charge transfer effects on the ions’ dynamics and kinetics are discussed in terms of partial charge models.     

Molecular dynamics estimates of ion diffusion in model hydrophobic and KcsA potassium channels

Molecular dynamics simulations are carried out to obtain estimates of diffusion coefficients of biologically important Na+, K+, Ca2+ and Cl- ions in hydrophobic cylindrical channels with varying radii and large reservoirs. Calculations for the cylindrical channels are compared to those for the KcsA potassium channel, for which the protein structure has recently been determined from X-ray diffraction experiments. Our results show that ion diffusion is maintained at reasonably high levels even within narrow channels, and does not support the very small diffusion coefficients used in some continuum models in order to fit experimental data. The present estimates of ion diffusion coefficients are useful in the calculation of channel conductance using the Poisson-Nernst-Planck theory, or Brownian dynamics.     

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).     

Electrostatic basis of valence selectivity in cationic channels

We examine how a variety of cationic channels discriminate between ions of differing charge. We construct models of the KcsA potassium channel, voltage gated sodium channel and L-type calcium channel, and show that they all conduct monovalent cations, but that only the calcium channel conducts divalent cations. In the KcsA and sodium channels divalent ions block the channel and prevent any further conduction. We demonstrate that in each case, this discrimination and some of the more complex conductance properties of the channels is a consequence of the electrostatic interaction of the ions with the charges in the channel protein. The KcsA and sodium channels bind divalent ions strongly enough that they cannot be displaced by other ions and thereby block the channel. On the other hand, the calcium channel binds them less strongly such that they can be destabilized by the repulsion of another incoming divalent ion, but not by the lesser repulsion from monovalent ions.     

Electrostatic basis of valence selectivity in cationic channels

We examine how a variety of cationic channels discriminate between ions of differing charge. We construct models of the KcsA potassium channel, voltage gated sodium channel and L-type calcium channel, and show that they all conduct monovalent cations, but that only the calcium channel conducts divalent cations. In the KcsA and sodium channels divalent ions block the channel and prevent any further conduction. We demonstrate that in each case, this discrimination and some of the more complex conductance properties of the channels is a consequence of the electrostatic interaction of the ions with the charges in the channel protein. The KcsA and sodium channels bind divalent ions strongly enough that they cannot be displaced by other ions and thereby block the channel. On the other hand, the calcium channel binds them less strongly such that they can be destabilized by the repulsion of another incoming divalent ion, but not by the lesser repulsion from monovalent ions.     

Permeation through an open channel: Poisson-Nernst-Planck theory of a synthetic ionic channel.

The synthetic channel [acetyl-(LeuSerSerLeuLeuSerLeu)3-CONH2]6 (pore diameter approximately 8 A, length approximately 30 A) is a bundle of six alpha-helices with blocked termini. This simple channel has complex properties, which are difficult to explain, even qualitatively, by traditional theories: its single-channel currents rectify in symmetrical solutions and its selectivity (defined by reversal potential) is a sensitive function of bathing solution. These complex properties can be fit quantitatively if the channel has fixed charge at its ends, forming a kind of macrodipole, bracketing a central charged region, and the shielding of the fixed charges is described by the Poisson-Nernst-Planck (PNP) equations. PNP fits current voltage relations measured in 15 solutions with an r.m.s. error of 3.6% using four adjustable parameters: the diffusion coefficients in the channel's pore DK = 2.1 x 10(-6) and DCl = 2.6 x 10(-7) cm2/s; and the fixed charge at the ends of the channel of +/- 0.12e (with unequal densities 0.71 M = 0.021e/A on the N-side and -1.9 M = -0.058e/A on the C-side). The fixed charge in the central region is 0.31e (with density P2 = 0.47 M = 0.014e/A). In contrast to traditional theories, PNP computes the electric field in the open channel from all of the charges in the system, by a rapid and accurate numerical procedure. In essence, PNP is a theory of the shielding of fixed (i.e., permanent) charge of the channel by mobile charge and by the ionic atmosphere in and near the channel's pore. The theory fits a wide range of data because the ionic contents and potential profile in the channel change significantly with experimental conditions, as they must, if the channel simultaneously satisfies the Poisson and Nernst-Planck equations and boundary conditions. Qualitatively speaking, the theory shows that small changes in the ionic atmosphere of the channel (i.e., shielding) make big changes in the potential profile and even bigger changes in flux, because potential is a sensitive function of charge and shielding, and flux is an exponential function of potential.     

A computational study of ion binding and protonation states in the KcsA potassium channel

We report results from microscopic molecular dynamics and free energy perturbation simulations of the KcsA potassium channel based on its experimental atomic structure. Conformational properties of selected amino acid residues as well as equilibrium positions of K(+) ions inside the selectivity filter and the internal water cavity are examined. Positions three and four (counting from the extracellular site) in the experimental structure correspond to distinctly separate binding sites for K(+) ions inside the selectivity filter. The protonation states of Glu71 and Asp80, which are close to each other and to the selectivity filter, as well as K(+) binding energies are determined using free energy perturbation calculations. The Glu71 residue which is buried inside a protein cavity is found to be most stable in the neutral form while the solvent exposed Asp80 is ionized. The channel altogether exothermically binds up to three ions, where two of them are located inside the selectivity filter and one in the internal water cavity. Ion permeation mechanisms are discussed in relation to these results.     

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K_v1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K_v1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) α-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K_v open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K_v1.5 channel is unable to close when bupivacaine is bound.     

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Voltage-dependent gating at the KcsA selectivity filter

The prokaryotic K(+) channel KcsA, although lacking a 'standard' voltage-sensing domain, shows voltage-dependent gating that leads to an increase in steady-state open probability of almost two orders of magnitude between +150 and -150 mV. Here we show that voltage-dependent gating in KcsA is associated with the movement of approximately 0.7 equivalent electronic charges. This charge movement produces an increase in the rate of entry into a long-lived inactivated state and seems to be independent of the proton-activation mechanism. Charge neutralization at position 71 renders the channel essentially voltage-independent by preventing entry into the inactivated state. A mechanism for voltage-dependent gating at the selectivity filter is proposed that is based on the reorientation of the carboxylic moiety of Glu71 and its influence in the conformational dynamics of the selectivity filter.     

Mechanism of intracellular block of the KcsA K+ channel by tetrabutylammonium: insights from X-ray crystallography, electrophysiology and replica-exchange molecular dynamics simulations

The mechanism of intracellular blockade of the KcsA potassium channel by tetrabutylammonium (TBA) is investigated through functional, structural and computational studies. Using planar-membrane electrophysiological recordings, we characterize the binding kinetics as well as the dependence on the transmembrane voltage and the concentration of the blocker. It is found that the apparent affinity of the complex is significantly greater than that of any of the eukaryotic K(+) channels studied previously, and that the off-rate increases with the applied transmembrane voltage. In addition, we report a crystal structure of the KcsA-TBA complex at 2.9 A resolution, with TBA bound inside the large hydrophobic cavity located at the center of the channel, consistent with the results of previous functional and structural studies. Of particular interest is the observation that the presence of TBA has a negligible effect on the channel structure and on the position of the potassium ions occupying the selectivity filter. Inspection of the electron density corresponding to TBA suggests that the ligand may adopt more than one conformation in the complex, though the moderate resolution of the data precludes a definitive interpretation on the basis of the crystallographic refinement methods alone. To provide a rationale for these observations, we carry out an extensive conformational sampling of an atomic model of TBA bound in the central cavity of KcsA, using the Hamiltonian replica-exchange molecular dynamics simulation method. Comparison of the simulated and experimental density maps indicates that the latter does reflect at least two distinct binding orientations of TBA. The simulations show also that the relative population of these binding modes is dependent on the ion configuration occupying the selectivity filter, thus providing a clue to the nature of the voltage-dependence of the binding kinetics.     

Probing the interaction of lipids with the non-annular binding sites of the potassium channel KcsA by magic-angle spinning NMR

The activity of the potassium channel KcsA is tightly regulated through the interactions of anionic lipids with high-affinity non-annular lipid binding sites located at the interface between the channel's subunits. Here we present solid-state phosphorous NMR studies that resolve the negatively charged lipid phosphatidylglycerol within the non-annular lipid-binding site. Perturbations in chemical shift observed upon the binding of phosphatidylglycerol are indicative of the interaction of positively charged sidechains within the non-annular binding site and the negatively charged lipid headgroup. Site directed mutagenesis studies have attributed these charge interactions to R64 and R89. Functionally the removal of the positive charges from R64 and R89 appears to act synergistically to reduce the probability of channel opening.     

Role of the KcsA channel cytoplasmic domain in pH-dependent gating

The KcsA channel is a representative potassium channel that is activated by changes in pH. Previous studies suggested that the region that senses pH is entirely within its transmembrane segments. However, we recently revealed that the cytoplasmic domain also has an important role, because its conformation was observed to change dramatically in response to pH changes. Here, to investigate the effects of the cytoplasmic domain on pH-dependent gating, we made a chimera mutant channel consisting of the cytoplasmic domain of the KcsA channel and the transmembrane region of the MthK channel. The chimera showed a pH dependency similar to that of KcsA, indicating that the cytoplasmic domain can act as a pH sensor. To identify how this region detects pH, we substituted certain cytoplasmic domain amino acids that are normally negatively charged at pH 7 for neutral ones in the KcsA channels. These mutants opened independently of pH, suggesting that electrostatic charges have a major role in the cytoplasmic domain's ability to sense and respond to pH.     

Voltage sensitivity and gating charge in Shaker and Shab family potassium channels.

The members of the voltage-dependent potassium channel family subserve a variety of functions and are expected to have voltage sensors with different sensitivities. The Shaker channel of Drosophila, which underlies a transient potassium current, has a high voltage sensitivity that is conferred by a large gating charge movement, approximately 13 elementary charges. A Shaker subunit's primary voltage-sensing (S4) region has seven positively charged residues. The Shab channel and its homologue Kv2.1 both carry a delayed-rectifier current, and their subunits have only five positively charged residues in S4; they would be expected to have smaller gating-charge movements and voltage sensitivities. We have characterized the gating currents and single-channel behavior of Shab channels and have estimated the charge movement in Shaker, Shab, and their rat homologues Kv1.1 and Kv2.1 by measuring the voltage dependence of open probability at very negative voltages and comparing this with the charge-voltage relationships. We find that Shab has a relatively small gating charge, approximately 7.5 e(o). Surprisingly, the corresponding mammalian delayed rectifier Kv2.1, which has the same complement of charged residues in the S2, S3, and S4 segments, has a gating charge of 12.5 e(o), essentially equal to that of Shaker and Kv1.1. Evidence for very strong coupling between charge movement and channel opening is seen in two channel types, with the probability of voltage-independent channel openings measured to be below 10(-9) in Shaker and below 4 x 10(-8) in Kv2.1.     

Influence of interionic interactions on functional state and blocker binding of voltage-gated potassium channels

A mechanism of ion conduction of a voltage-gated potassium channel KcsA was investigated in full-atomic approximation at a trajectory length of 100 ns using the Lomonosov supercomputer. Methods of molecular dynamics were employed. A structure of the KcsA channel in the open state obtained by X-ray structure analysis (PDB ID 3fb7) was used. Free energy profiles of the KcsA pore occupied with either one or three potassium ions were calculated. It was shown that, under physiological conditions, ions pass through the channel pore cooperatively and the mechanism most probably includes three ions permeating in concert. Interactions of the mammalian voltage-gated channel Kv1.2 with neurotoxin were investigated. It was demonstrated that the effect of interionic interactions on binding of a blocker is rather insufficient.     

On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps

Partial charges of atoms in a molecule and electrostatic potential (ESP) density for that molecule are known to bear a strong correlation. In order to generate a set of point‐field force field parameters for molecular dynamics, Kollman and coworkers have extracted atomic partial charges for each of all 20 amino acids using restrained partial charge‐fitting procedures from theoretical ESP density obtained from condensed‐state quantum mechanics. The magnitude of atomic partial charges for neutral peptide backbone they have obtained is similar to that of partial atomic charges for ionized carboxylate side chain atoms. In this study, the effect of these known atomic partial charges on ESP is examined using computer simulations and compared with the experimental ESP density recently obtained for proteins using electron microscopy. It is found that the observed ESP density maps are most consistent with the simulations that include atomic partial charges of protein backbone. Therefore, atomic partial charges are integral part of atomic properties in protein molecules and should be included in model refinement.     

Structural differences of bacterial and mammalian K+ channels

Using a peptide toxin, kaliotoxin (KTX), we gained new insight into the topology of the pore region of a voltage-gated potassium channel, mKv1.1. In order to find new interactions between mKv1.1 and KTX, we investigated the pH dependence of KTX block which was stronger at pH(o) 6.2 compared with pH(o) 7.4. Using site-directed mutagenesis on the channel and the toxin, we found that protonation of His(34) in KTX caused the pH(o) dependence of KTX block. Glu(350) and Glu(353) in mKv1.1, which interact with His(34) in KTX, were calculated to be 4 and 7 A away from His(34)/KTX, respectively. Docking of KTX into a homology model of mKv1.1 based on the KcsA crystal structure using this and other known interactions as constraints showed structural differences between mKv1.1 and KcsA within the turret (amino acids 348-357). To satisfy our data, we would have to modify the KcsA crystal structure for the mKv1.1 channel orienting Glu(350) 7 A and Glu(353) 4 A more toward the center of the pore compared with KcsA. This would place Glu(350) 15 A and Glu(353) 11 A away from the center of the pore instead of the distances for the equivalent KcsA residues with 22 A for Gly(53) and 15 A for Gly(56), respectively. Bacterial and mammalian potassium channels may have structural differences regarding the turret of the outer pore vestibule. This topological difference between both channel types may have substantial influence on structure-guided development of new drugs for mammalian potassium channels by rational drug design.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Functional evidence for a supramolecular structure for the Streptomyces lividans potassium channel KcsA

Here we present functional evidence for involvement of poly-(R)-3-hydroxybutyrate (PHB) and inorganic polyphosphate (polyP) in ion conduction and selection at the intracellular side of the Streptomyces lividans potassium channel, KcsA. At < or = 25 degrees C, KcsA forms channels in planar bilayers that display signal characteristics of PHB/polyP channels at the intracellular side; i.e., a preference for divalent Mg(2+) cations at pH 7.2, and a preference for monovalent K+ cations at pH 6.8. Between 25 and 26 degrees C, KcsA undergoes a transition to a new conformation in which the channel exhibits high selectivity for K+, regardless of solution pH. This suggests that basic residues of the C-terminal polypeptides have moved closer to the polyP end unit, reducing its negative charge. The data support a supramolecular structure for KcsA in which influx of ions is prevented by the selectivity pore, whereas efflux of K+ is governed by a conductive core of PHB/polyP in partnership with the C-terminal polypeptide strands.     

Local and global structure of the monomeric subunit of the potassium channel KcsA probed by NMR

KcsA is a homotetrameric 68-kDa membrane-associated potassium channel which selectively gates the flux of potassium ions across the membrane. The channel is known to undergo a pH-dependent open-to-closed transition. Here we describe an NMR study of the monomeric subunit of the channel (KcsAM), solubilized in SDS micelles. Chemical shift, solvent exchange, backbone 15N relaxation and residual dipolar coupling (RDC) data show the TM1 helix to remain intact, but the TM2 helix contains a distinct kink, which is subject to concentration-independent but pH-dependent conformational exchange on a microsecond time scale. The kink region, centered at G99, was previously implicated in the gating of the tetrameric KcsA channel. An RDC-based model of KcsAM at acidic pH orients TM1 and the two helical segments of the kinked TM2 in a configuration reminiscent of the open conformation of the channel. Thus, the transition between states appears to be an inherent capability of the monomer, with the tetrameric assembly exerting a modulatory effect upon the transition which gives the channel its physiological gating profile.