ABSTRACT: Based on a homology model of the Kv1.3 potassium channel, the recognitions of the six scorpion toxins, viz. agitoxin2, charybdotoxin, kaliotoxin, margatoxin, noxiustoxin, and Pandinus toxin, to the human Kv1.3 potassium channel have been investigated by using an approach of the Brownian dynamics (BD) simulation integrating molecular dynamics (MD) simulation. Reasonable three-dimensional structures of the toxin-channel complexes have been obtained employing BD simulations and triplet contact analyses. All of the available structures of the six scorpion toxins in the Research Collaboratory for Structural Bioinformatics Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformations of the toxin significantly affect both the molecular recognition and binding energy between the two proteins. BD simulations predicted that all the six scorpion toxins in this study use their beta-sheets to bind to the extracellular entryway of the Kv1.3 channel, which is in line with the primary clues from the electrostatic interaction calculations and mutagenesis results. Additionally, the electrostatic interaction energies between the toxins and Kv1.3 channel correlate well with the binding affinities (-logK(d)s), R(2) = 0.603, suggesting that the electrostatic interaction is a dominant component for toxin-channel binding specificity. Most importantly, recognition residues and interaction contacts for the binding were identified. Lys-27 or Lys-28, residues Arg-24 or Arg-25 in the separate six toxins, and residues Tyr-400, Asp-402, His-404, Asp-386, and Gly-380 in each subunit of the Kv1.3 potassium channel, are the key residues for the toxin-channel recognitions. This is in agreement with the mutation results. MD simulations lasting 5 ns for the individual proteins and the toxin-channel complexes in a solvated lipid bilayer environment confirmed that the toxins are flexible and the channel is not flexible in the binding. The consistency between the results of the simulations and the experimental data indicated that our three-dimensional models of the toxin-channel complex are reasonable and can be used as a guide for future biological studies, such as the rational design of the blocking agents of the Kv1.3 channel and mutagenesis in both toxins and the Kv1.3 channel. Moreover, the simulation result demonstrates that the electrostatic interaction energies combined with the distribution frequencies from BD simulations might be used as criteria in ranking the binding configuration of a scorpion toxin to the Kv1.3 channel.
