The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http: //www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its β-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.

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 b-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 (-logKds), R2 = 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.

to evaluate the second-order nonlinear susceptibilities for a series of fulleropyrrolidine-tetrathiafulvalene(C60PY-TTF) dyads, C60PY, and C60PY-(CH=CH)n-TTF (n = 0-6). The equilibrium geometries, electronic structures, and UV-vis spectra of these C60PY-TTFs D-B-A dyads were determined by using ZINDO methods. On the basis of electronic spectra, we calculated the β for the dyads using the ZINDO-SOS expression, and found that these dyads have relatively large β. The calculated values of p are 48.78, 62.34, 97.90, 123.24, 183.78, 189.43, 193.45, 195.82 × 10-30 esu for C60PY and C60PY-(CH=CH)n-TTF (n=0-6), respectively, which shows the enhancing effect when TTF is introduced. The β significantly increases as n increases from 0-3; and it tends to be saturated in the end.

The equilibrium geometries, electronic structures and UV-vis spectra of a series of spiroannelated quinone-type methanofullerenes have been determined by using Zerner’s intermediate neglect of differential overlap method. The results show that between fullerene and the addend there exists a special interaction, “periconjugation”, which results in through-space orbital interactions. The calculated UV-vis spectra are in good agreement with experiments. On the basis of the electronic spectra, the β values are calculated. The results show that spiroannelated quinone-type methanofullerenes have quite large β values. We attribute the large β values to both the charge transfer from C60 to benzoquinone and on the C60 three-dimensional conjugated sphere.

Using density functional theory and ZINDO methods, the structures, molecular orbitals and electronic spectra of a series of spirobifluorenes were investigated, in which two identical push-pull type chromophores are connected via an insulating carbon. Based on the correct electronic spectra, according to the sum-over-states formula, the calculation program was devised and nonlinear third-order optical susceptibilities were calculated. The results show a good agreement between the calculation and experiment. The calculated nonlinear third-order susceptibility γ of the spirobifluorenes are about six times as large as those of the corresponding mono-fluorenes without an increase of λmax. Since electronic coupling between monomer units in the spirobifluorenes is negligible owing to symmetry, we attribute the large enhancement to the spiroconjugation. Thus, we conclude that derived spirobifluorenes will be a hopeful type of third-order nonlinear optical material from the standpoint of the high transparency and relatively large γ value.