Xiao-Jian Tan

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.

The interaction model between metal cation and tropylium: a quantum chemistry predication

Abstract The aromatic cation tropylium, C 7 H 7 + , predicted at the MP2/6-31G** level, is capable of binding with metal cations Be 2+ or Mg 2+ , forming M 2+ –C 7 H 7 + complexes. The obstacle for their binding is almost electrostatic repulsion, and the binding is from polarization and charge transfer. The orbital interaction between the M 2+ and C 7 H 7 + is mainly the s–π and p–π interactions. Interestingly, Be 2+ is possible to pass through the ring of C 7 H 7 + , while Mg 2+ is not. The intrinsic IR band of the M 2+ –C 7 H 7 + complex is below 600 cm −1 , which results from the vibration of the M 2+ along the normal axis of C 7 H 7 + .

Density functional theory (DFT) study on the interaction of ammonium (NH4+) and aromatic nitrogen heterocyclics

A DFT calculation was performed at the B3LYP/6-31G* level on the complexes formed by NH4+ and aromatic nitrogen heterocyclics, viz. pyrrole, imidazole, pyridine and indole, in order to investigate the mechanism and complexity of the interaction between the ammonium group and the aromatic heterocyclic in biomacromolecules. The optimized geometries suggested that there are two different types of complexes: one is a cation–π complex and the other is a hydrogen bond complex. A cation–π complex will be formed if the heteroatom has no localized lone-pair electrons. A hydrogen bond complex will be formed by proton transfer from NH4+ to the heteroatom if the heteroatom has localized lone-pair electrons. In the case of the cation–π complex, the predicted geometries, atomic charges and thermodynamic parameters revealed that ammonium binds more strongly to heterocyclics than it binds to benzene. The calculated orbital coefficient and the optimized structures implied that NH4+ interacts with the π electrons of the CC bond of heterocyclics to form a cation–π complex mainly through one hydrogen atom. Regarding the hydrogen bond complex, although the calculated binding strength is similar to that for the cation–π complex, the ΔH of the whole reaction process suggested that the formation of the hydrogen bond complex is favorable to the stability of the whole system. Calculated IR spectra showed that three groups of new bands appear when NH4+ binds to heterocyclics. Normal mode analysis showed that these new bands are all related to the relative motion of the two parts in the formed complexes. All these results suggest that the NH4+–heterocyclic system is a better model for studying the nature and complexity of the interaction between the ammonium group and the aromatic ring structure in biomolecules.

STRUCTURAL FEATURE OF AChE INHIBITOR HUPERZINE B IN NATURE AND IN THE BINDING SITE OF AChE: DENSITY FUNCTIONAL THEORY STUDY COMBINED WITH IR DETERMINATION

Quantum chemical DFT-B3LYP/6-31G method and IR spectrometry have been used to investigate the natural and binding structures of Huperzine B (HupB) in order to better understand the interaction nature between acetylcholinesterase (AChE) and its inhibitor, with the view of designing new AChE inhibitors. The predicted and experimental results reveal that both the natural state and binding form of HupB adopt the chair conformation. Furthermore, the B3LYP/6-31G results suggest that structure S1 should be the dominant form of the two possible chair structures (S1 and S2, Fig. 2). The calculated results also show that the condensed ring structure composing of rings A, B and C is very rigid. Therefore, its flexibility does not need to be considered when we try to dock this structure to its target. Indeed, this supposition is conrmed by the excellent alignment of the binding structure produced from our recent X-ray crystallographic structure of the HupB-AChE complex with the B3LYP/6-31G predicted geometry. Among all the 111 predicted vibrational bands, the mode 110, which is resulted from the stretching of the bond N2{H and having the second highest frequency, is essential for the geometrical identication. The dierence between our predicted strongest absorption band and experimental IR spectrum suggests

Quantum chemical HF/4-31G calculations on buckminsterfullerene intermediates

Quantum chemical ab initio (U)HF/4-31G investigation on buckminsterfullerene and some proposed intermediates in its formation is carried out in this study with a view to better understanding how small carbon species carry out self-assembly to form fullerenes. The calculations on 19 carefully designed fullerene intermediates reveal that the core of an intermediate, rather than the number of its dangling bonds or abutting pentagon rings, has an intrinsic effect on its energy. The computational results show that hexagonal-core structures have lower energies than pentagonal-core structures. In addition, the pentagonal core enclosed completely by hexagonal rings has the highest energy. The UHF/4-31G results also suggest that some intermediates such as C18, C21 and C30 with hexagonal cores have unusually low energies in comparison with their isomers or neighbours. Based on these calculated results, we outline the possible pathways from precursor to intermediates to fullerenes, subject to synthesis conditions and raw materials. These pathways support some existing proposals, such as medium monocyclic ring stacking and small ring polymerization mechanisms. However, our results do not suggest that the numbers of dangling bonds or abutting pentagonal rings have the highest impact on fullerene formation. The calculated thermodynamic parameters of the dimerization and addition reactions between two bowl-shaped intermediates suggest that these reactions are favorable to fullerene formation, and that the concentration of bowl-shaped fullerene intermediates should be very low in all detectable carbon species.