Incidence of partial charges on ion selectivity in potassium channels

Chloride Ion Transport by the E. coli CLC Cl-/H+ Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study

We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl^- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl^-/H⁺ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl^-. Moving the anion into the external binding site (S_ext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl^- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl^- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl^- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl^- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl^- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl^- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl^- transport.

K+/Na+ selectivity in toy cation binding site models is determined by the 'host'

The macroscopic ion-selective behavior of K(+) channels is mediated by a multitude of physiological factors. However, considering the carbonyl-lined binding site of a conductive K(+) channel as a canonical eightfold coordinated construct can be useful in understanding the principles that correlate the channel's structure with its function. We probe the effects of structure and chemical composition on the K(+)/Na(+) selectivity provided by a variety of simplified droplet-like ion binding site models. We find that when carbonyl- and water-based models capture the qualitative structural features of the K(+) channel binding site, a selective preference for K(+) emerges. Thus our findings suggest that the preference for K(+) over Na(+) exhibited by such models is principally built-in, and is not due to a unique K(+)-selective property of carbonyl functional groups. This suggestion is confirmed by a general thermodynamic assessment, which provides a basis for using simplified models to study the design principles underlying the molecular evolution of K(+) channels.

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.

Many-body effects and simulations of potassium channels

The electronic polarizability of an ion or a molecule is a measure of the relative tendency of its electron cloud to be distorted from its normal shape by an electric field. On the molecular scale, in a condensed phase, any species sits in an electric field due to its neighbours, and the resulting polarization is an important contribution to the total interaction energy. Electrostatic interactions are crucial for determining the majority of chemical–physical properties of the system and electronic polarization is a fundamental component of these interactions. Thus, polarization effects should be taken into account if accurate descriptions are desired. In classical computer simulations, the forces required to drive the system are typically based on interatomic interaction potentials derived in part from electronic structure calculations or from experimental data. Owing to the difficulties in including polarization effects in classical force fields, most of them are based just on pairwise additive interaction potentials. At present, major efforts are underway to develop polarizable interaction potentials for biomolecular simulations. In this review, various ways of introducing explicit polarizability into biomolecular models and force fields are reviewed, and the progress that might be achieved in applying such methods to study potassium channels is described.

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

Quantum mechanical calculations on selectivity in the KcsA channel: the role of the aqueous cavity

We have carried out quantum calculations on selected residues at the intracellular side of the selectivity filter of the KcsA potassium channel, using the published X-ray coordinates as starting points. The calculations involved primarily the side chains of residues lining the aqueous cavity on the intracellular side of the selectivity filter, in addition to water molecules, plus a K+ or Na+ ion. The results showed unambiguously that Na+ significantly distorts the symmetry of the channel at the entrance to the selectivity filter (at the residue T75), while K+ does so to a much smaller extent. In all, three ion positions have been calculated: the S4 (lowest) position at the bottom of the selectivity filter, the top of the cavity, and the midpoint of the cavity; Na+ is trapped at the cavity top, while K+ is cosolvated by the selectivity filter carbonyl groups plus threonine hydroxyl groups so that it can traverse the filter. Only one water molecule remains in the K+ solvation shell at the upper position in the cavity; this solvation shell also contains four threonine (T75) hydroxyl oxygens and two backbone carbonyls, while Na+ is solvated by five molecules of water and one oxygen from threonine hydroxyls. T75 at the entrance to the selectivity filter has a key role in recognition of the alkali ion, and T74 has secondary importance. The energetic basis for the preferential bonding of potassium by these residues is briefly discussed, based on additional calculations. Taken together, the results suggest that Na+ would have difficulty entering the cavity, and if it did, it would not be able to enter the selectivity filter.

Selectivity in K+ channels is due to topological control of the permeant ion's coordinated state

The selectivity filter of K+ channels provides specific coordinative interactions between dipolar carbonyl ligands, water, and the permeant cation, which allow for selective flow of K+ over (most importantly) Na+ across the cell membrane. Although a structural viewpoint attributes K+ selectivity to coordination geometry provided by the filter, recent molecular dynamics simulation studies attribute it to dynamic and unique chemical/electrostatic properties of the filter's carbonyl ligands. Here we provide a simple theoretical analysis of K+ and Na+ complexation with water in the context of simplified binding site models and bulk solution. Our analysis reveals that water molecules and carbonyl groups can both provide K+ selective environments if equivalent constraints are imposed on the coordination number of the complex. Absence of such constraints annihilates selectivity, demonstrating that whether a coordinating ligand is a water molecule or a carbonyl group, "external" or "topological" constraints/forces must be imposed on an ion-coordinated complex to elicit selective binding. These forces must inevitably originate from the channel protein, because in bulk water, which, by definition, presents a nonselective medium, the coordination number is allowed to relax to suit the ion. We show that the coordination geometry of K+ channel binding sites is replicated by 8-fold complexation of K+ in both water and simplified binding site models due to dominance of local interactions within a complex and is thus a requirement for topologically constraining the coordination number to a specific value.

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.

Insight into the origins of the barrier-less knock-on conduction in the KcsA channel: molecular dynamics simulations and ab initio calculations

Since the pioneering work of Zhou et al. (Y. Zhou, J. H. Morais-Cabral, A. Kaufman and R. MacKinnon, Nature, 2001, 414, 43-48) it is now well established that the streptomyces lividans potassium channel (KcsA) can accommodate more than one ion, namely between 2 and 3. As a result, it is usually assumed that the conduction of ions proceeds through a barrier-less knock-on mechanism. This one is an alternation of two sequences containing either 2 or 3 ions which have nearly the same energies. However, the origin of such knock-on mechanism is not clearly known. The entry and the exit of ion in or out of the selectivity filter are suspected to be due to the repulsive interactions between ions. In this work, molecular dynamics simulations running over nanoseconds have been done in order to identify such events. Two specific situations, namely (S(1), S(3)) containing 2 ions and (S(2), S(4)) containing 3 ions, have been investigated regarding the different locations that ions can occupy during their diffusion through the selectivity filter of KcsA. We show that contractions of the (S(1), S(3)) file and dilation of the (S(2), S(4)) file are at the origin of the passage from one sequence to the other. The comparison between the experimentally observed diffusion rate and the occurrence's frequency of such contractions or dilation confirm the importance of such events. Ab initio calculations have also been conducted in order to examine the effect of ion polarization in the filter of KcsA. During the contraction of the ion/water file, one charge at the extra-cellular mouth of the channel strongly deviates from the others. This behavior could guide the diffusion direction to a certain extent since the contraction of the (S(1), S(3)) is favored.

Ion selectivity in potassium channels

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.

Mechanism of Ion Permeation in a Model Channel: Free Energy Surface and Dynamics of K+Ion Transport in an Anion-Doped Carbon Nanotube

The mechanism of the ion permeation is investigated for an anion-doped carbon nanotube, as a model of the K+ channel, by analyzing the free energy surface and the dynamics of the ion permeation through the model channel. It is found that the main rate-determining step is how an ion enters the channel. The entrance of the ion is mostly blocked by a water molecule located at this entrance. Only about 10% of K+ ions which reach the mouth of the channel can really enter the channel. The rejection rate sensitively depends on the location of this water molecule, which is easily controlled by the charge of the carbon nanotube; for example, the maximum permeation is obtained when the anion charge is at a certain value, -5.4e in the present model. At this charge, the facile translocation of the ion inside the channel is also induced due to the number of fluctuations of the ions inside the channel. Therefore, the so-called "Newton's balls", a toy model, combined with a simple ion diffusion model for explaining the fast ion permeation should be modified. The present analysis thus suggests that there exists an optimum combination of the length and the charge of the carbon nanotube for the most efficient ion permeation.