Affinity and selectivity of ShK toxin for the Kv1 potassium channels from free energy simulations

The voltage-gated potassium channel Kv1.3 is an attractive target for treatment of autoimmune diseases. ShK toxin from sea anemone is one of the most potent blockers of Kv1.3, and therefore ShK and its analogues have been proposed as therapeutic leads for such diseases. Increasing the selectivity of the proposed leads for Kv1.3 over other Kv1 channels is a major issue in this endeavor. Here we study binding of ShK toxin to Kv1 channels using free energy simulation methods. Homology models for Kv1.1 and Kv1.3 channels are constructed using the crystal structure of Kv1.2. The initial poses for the Kv1.x-ShK complexes are obtained using HADDOCK, which are then refined via molecular dynamics simulations. The binding mode in each complex is characterized by identifying the strongly interacting residues, which compare well with available mutagenesis studies. For each complex, the potential of mean force is calculated from umbrella sampling simulations, and the corresponding absolute binding free energy is determined. The computed binding free energies are in good agreement with the experimental data, which increases the confidence on the model complexes. The insights gained on Kv1.x-ShK binding modes will be valuable in the development of new ShK analogues with better selectivity properties.