Modeling P-loops domain of sodium channel: Homology with potassium channels and interaction with ligands

A large body of experimental data on Na+ channels is available, but the interpretation of these data in structural terms is difficult in the absence of a high-resolution structure. Essentially different electrophysiological and pharmacological properties of Na+ and K+ channels and poor identity of their sequences obstruct homology modeling of Na+ channels. In this work, we built the P-loops model of the Na+ channel, in which the pore helices are arranged exactly as in the MthK bacterial K+ channel. The conformation of the selectivity-filter region, which includes residues in positions - 2 through +4 from the DEKA locus, was shaped around rigid molecules of saxitoxin and tetrodotoxin that are known to form multiple contacts with this region. Intensive Monte Carlo minimization that started from the MthK-like conformation produced practically identical saxitoxin- and tetrodotoxin-based models. The latter was tested to explain a wide range of experimental data that were not used at the model building stage. The docking of tetrodotoxin analogs unambiguously predicted their optimal orientation and the interaction energy that correlates with the experimental activity. The docking of mu-conotoxin produced a binding model consistent with experimentally known toxin-channel contacts. Monte Carlo-minimized energy profiles of tetramethylammonium pulled through the selectivity-filter region explain the paradoxical experimental data that this organic cation permeates via the DEAA but not the AAAA mutant of the DEKA locus. The model is also consistent with earlier proposed concepts on the Na+ channel selectivity as well as Ca2+ selectivity of the EEEE mutant of the DEKA locus. Thus, the model integrates available experimental data on the Na+ channel P-loops domain, and suggests that it is more similar to K+ channels than was believed before.