Abstract
Five homology models for honeybee (Apis mellifera) nicotinic acetylcholine receptor (nAChR) α1/β1, α3/β2, α4/β2, α6/β2 and α9/α9 subtypes were built from the Torpedo marmorata nAChR X-ray structure. Then, imidacloprid, fipronil and their metabolites were docked into the ligand binding domain (LBD) of these receptors and the corresponding scoring functions were calculated. The binding modes of the docked compounds were carefully analysed. Finally, multivariate analyses were used for deriving structure-activity relationships based on hydrogen bond number and scoring functions between the insecticides and the nAChR models.