Investigation of the selectivity of one type of small-molecule inhibitor for three Na_v channel isoforms based on the method of computer simulation

Voltage-gated sodium (Na_v) channels play a pivotal role for the changes in membrane potential and belong to large membrane proteins that compose four voltage sensor domains (VSD1–4). In this study, we describe the binding mode and selectivity of one of the aryl sulfonamide sodium channel inhibitors, PF-04856264, for the VSD4s in Na_v1.4, Na_v1.5 and Na_v1.7, respectively, through molecular dynamics simulation and enhanced post-dynamics analyses. Our results show that there are three binding site regions (BSR1–3) in the combination of the ligand and receptors, of which BSR1 and BSR3 contribute to the selectivity and affinity of the ligand to the receptor. What’s more, the 39th residue (Y39 in VSD4^hNav1.4/ VSD4^hNav1.7 and A39 in VSD4^hNav1.5) and N42 in BSR1, the 84th residue (L84 in VSD4^hNav1.4, T84 in VSD4^hNav1.5, and M84 in VSD4^hNav1.7) in BSR2 and the conserved positive charged residues in BSR3 have major contributions to the interaction between the ligand and receptor. Further analysis reveals that if the 39th residue has a benzene ring structure, the connection of BSR1 and the ligand would be much stronger through π-stacking interaction. On the other hand, the strength and number of the hydrogen bonds formed by the ligand and the conserved arginines on S4 determine the contribution of BSR3 to the total free binding energy. We anticipate this study pave the way for the design of more effective and safe treatment for pain that selectively target Na_v1.7.