Understanding subtype specific ion channel pore blockage by natural peptide-based toxins is crucial for developing such compounds into promising drug candidates. Herein, docking and molecular dynamics simulations were employed in order to understand the dynamics and binding states of the µ-conotoxins, PIIIA, SIIIA, and GIIIA, at the voltage-gated potassium channels of the KV1 family, and they were correlated with their experimental activities recently reported by Leipold et al. Their different activities can only adequately be understood when dynamic information about the toxin-channel systems is available. For all of the channel-bound toxins investigated herein, a certain conformational flexibility was observed during the molecular dynamic simulations, which corresponds to their bioactivity. Our data suggest a similar binding mode of µ-PIIIA at KV1.6 and KV1.1, in which a plethora of hydrogen bonds are formed by the Arg and Lys residues within the α-helical core region of µ-PIIIA, with the central pore residues of the channel. Furthermore, the contribution of the K+ channel's outer and inner pore loops with respect to the toxin binding. and how the subtype specificity is induced, were proposed.
Keywords: ion channel; potassium channel; toxin interactions; µ-GIIIA; µ-PIIIA; µ-SIIIA; µ-conotoxins.