In this study, Quantitative Structure‐Activity/property relationships (QSAR/QSPR) by means of multiple linear regressions (MLR) was performed to investigate the relationship between the 48 compounds of Temephos (Tem) drivatives and their bioactivities against acetylcholinesterase (AChE) of Tribolium castaneum. Two compounds Propylene‐N, N′‐bis (O, O′‐di Ethyl Phosphorothioat (12) and Ethylene N, ′N‐ di Methyl bis (O, O′‐ di Phenyl Phosphoramidate (24) had the most mortality on the T. castaneum. The compound 24 (IC50=34.54 ppm) is a good alternative to Tem. Also, QSAR calculations indicated that the electrostatic characteristics of the most effective insecticide are applied. In docking data, Tem derivatives with the backbone of P (O)‐NH−P(O), P(O)‐NH‐NH−P(O) and P(O)‐X−P(O) are located in the active site gorge of both AChE and butyrylcholinesterase (BChE) so as to maximize the favorable contacts. These compounds relate to enzymes by non‐covalent interactions such as hydrogen bonding, electrostatic and hydrophobic (as Trp82 and Trp286). Also, these cannot reach the end of the hole because of compounds are locked in the middle between the peripheral and acyl pocket site in AChE and BChE enzymes. The results of MLR and GA‐QSAR/QSPR model of human ChE showed that topological parameters affect the inhibitory potencies of the compounds on both of logK and p(IC50). logK and p(IC50) results have a good linear relationship.