Molecular Dynamics of Ion Conduction through the Selectivity Filter of the Na_VAb Sodium Channel

The determination of the atomic structures of voltage-gated bacterial sodium channels using X-ray crystallography has provided a first view of this family of membrane proteins. Molecular dynamics simulations offer one approach to clarify the underlying mechanism of permeation and selectivity in these channels. However, it appears that the intracellular gate of the pore domain is either closed or only open partially in the available X-ray structures. The lack of structure with a fully open intracellular gate poses a special challenge to computational studies aimed at simulating ion conduction. To circumvent this problem, we simulated a model of the Na_VAb channel in which the transmembrane S5 and S6 helices of the pore domain have been truncated to provide direct open access to the intracellular entryway to the pore. Molecular dynamics simulations were carried out over a range of membrane potential and ion concentration of sodium and potassium. The simulations show that the Na_VAb selectivity filter is essentially a cationic pore supporting the conduction of ions at a rate comparable to aqueous diffusion with no significant selectivity for sodium. Conductance and selectivity vary as a function of ion concentration for both cations. Permeation occurs primarily via a knock-on mechanism for both sodium and potassium, although the ion ordering in single file along the pore is not strictly maintained. The character of the outward current appears quite different from the inward current, with a buildup on ions in the selectivity filter prior to escape toward the extracellular side, indicating the presence of a rectification effect that is overcome by nonphysiological applied voltages.