Brownian dynamics simulations of interaction between scorpion toxin Lq2 and potassium ion channel

The association of the scorpion toxin Lq2 and a potassium ion (K(+)) channel has been studied using the Brownian dynamics (BD) simulation method. All of the 22 available structures of Lq2 in the Brookhaven Protein Data Bank (PDB) determined by NMR were considered during the simulation, which indicated that the conformation of Lq2 affects the binding between the two proteins significantly. Among the 22 structures of Lq2, only 4 structures dock in the binding site of the K(+) channel with a high probability and favorable electrostatic interactions. From the 4 candidates of the Lq2-K(+) channel binding models, we identified a good three-dimensional model of Lq2-K(+) channel complex through triplet contact analysis, electrostatic interaction energy estimation by BD simulation and structural refinement by molecular mechanics. Lq2 locates around the extracellular mouth of the K(+) channel and contacts the K(+) channel using its beta-sheet rather than its alpha-helix. Lys27, a conserved amino acid in the scorpion toxins, plugs the pore of the K(+) channel and forms three hydrogen bonds with the conserved residues Tyr78(A-C) and two hydrophobic contacts with Gly79 of the K(+) channel. In addition, eight hydrogen-bonds are formed between residues Arg25, Cys28, Lys31, Arg34 and Tyr36 of Lq2 and residues Pro55, Tyr78, Gly79, Asp80, and Tyr82 of K(+) channel. Many of them are formed by side chains of residues of Lq2 and backbone atoms of the K(+) channel. Thirteen hydrophobic contacts exist between residues Met29, Asn30, Lys31 and Tyr36 of Lq2 and residues Pro55, Ala58, Gly79, Asp80 and Tyr82 of the K(+) channel. These favorable interactions stabilize the association between the two proteins. These observations are in good agreement with the experimental results and can explain the binding phenomena between scorpion toxins and K(+) channels at the level of molecular structure. The consistency between the BD simulation and the experimental data indicates that our three-dimensional model of Lq2-K(+) channel complex is reasonable and can be used in further biological studies such as rational design of blocking agents of K(+) channels and mutagenesis in both toxins and K(+) channels.