A computational study of barium blockades in the KcsA potassium channel based on multi-ion potential of mean force calculations and free energy perturbation

Electrophysiological studies have established that the permeation of Ba²⁺ ions through the KcsA K⁺-channel is impeded by the presence of K⁺ ions in the external solution, while no effect is observed for external Na⁺ ions. This Ba²⁺ “lock-in” effect suggests that at least one of the external binding sites of the KcsA channel is thermodynamically selective for K⁺. We used molecular dynamics simulations to interpret these lock-in experiments in the context of the crystallographic structure of KcsA. Assuming that the Ba²⁺ is bound in site S₂ in the dominant blocked state, we examine the conditions that could impede its translocation and cause the observed “lock-in” effect. Although the binding of a K⁺ ion to site S₁ when site S₂ is occupied by Ba²⁺ is prohibitively high in energy (>10 kcal/mol), binding to site S₀ appears to be more plausible (ΔG > 4 kcal/mol). The 2D potential of mean force (PMF) for the simultaneous translocation of Ba²⁺ from site S₂ to site S₁ and of a K⁺ ion on the extracellular side shows a barrier that is consistent with the concept of external lock-in. The barrier opposing the movement of Ba²⁺ is very high when a cation is in site S₀, and considerably smaller when the site is unoccupied. Furthermore, free energy perturbation calculations show that site S₀ is selective for K⁺ by 1.8 kcal/mol when S₂ is occupied by Ba²⁺. However, the same site S₀ is nonselective when site S₂ is occupied by K⁺, which shows that the presence of Ba²⁺ affects the selectivity of the pore. A theoretical framework within classical rate theory is presented to incorporate the concentration dependence of the external ions on the lock-in effect.