Determination of the charge profile in the KcsA selectivity filter using ab initio calculations and molecular dynamics simulations.

The charge profile of K+ and Na+ ions moving in a single file across the filter of the KcsA channel is determined using both molecular dynamics simulations and ab initio calculations. We show a strong correlation between the charge variation and the ion location resulting in a saw-tooth profile, which provides additional information on the influence of charge transfer on the permeation and selectivity of the channel.

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.

Incidence of partial charges on ion selectivity in potassium channels.

Potassium channels are membrane proteins known to select potassium over sodium ions at a high diffusion rate. We conducted ab initio calculations on a filter model of KcsA of about 300 atoms at the Hartree-Fock level of theory. Partial charges were derived from the quantum mechanically determined electrostatic potential either with Merz-Kollman or Hinsen-Roux schemes. Large polarization and/or charge transfer occur on potassium ions located in the filter, while the charges on sodium ions remain closer to unity. As a result, a weaker binding is obtained for K(+) ions. Using a simplified version of a permeation model based on the concerted-motion mechanism for ion translocation within the single-file ion channel [P. H. Nelson, J. Chem. Phys. 117, 11396 (2002)], we discuss how differences in polarization effects in the adducts with K(+) and Na(+) can play a role as for ionic selectivity and conductance.

Targeted molecular dynamics of an open-state KcsA channel.

Pore opening of KcsA channel is studied using targeted molecular dynamics simulations. Conformational changes of the protein are determined, starting from the crystallized refined 2.0 A structure (pdb 1K4C) determined in x-ray experiments and arriving to the open-state structure constructed on the basis of electron paramagnetic resonance spectroscopic data (pdb 1JQ1). Our results corroborate the essential role played by the terminal residues located on the transmembrane helices M2 which were not taken into account at that time. The aperture mechanism of the channel appears to be ziplike. A small constraint (approximately equal to 5 x 10(-2) kcal mol(-1) A(-2) per C(alpha)) applied to the terminal residues located on the intracellular side is sufficient to initialize the pore opening at the innermost part of the gate, but additional constraint must be applied to definitely complete the pore aperture. The open structure is proved to be a metastable state since releasing the constraint leads to another relaxed open conformation which seems to reach stability.

Molecular dynamics study of the KcsA channel at 2.0-A resolution: stability and concerted motions within the pore.

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-A resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K+ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K+ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K+ ions and water molecules displacements shows concerted file motions on the simulated time scale (approximately 1 ns), which could act as precursor to the diffusion of K+ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.