Essential oils (EOs) and their volatile secondary metabolites have been proved to be effective on storage pests control, while restricted on the application due to unclear mechanism. Molecular dynamics (MD) simulations and binding free energies analysis provided an effective approach to reveal mechanism on conformational calculation. In this work, the insecticidal and repellent capacities of Praxelis clematidea and Ageratum houstonianum oils and their main components identified by gas chromatography–mass spectrometry (GC–MS) were scientifically measured. Interestingly, P. clematidea oil exhibited strong fumigant toxicity against Tribolium castaneum (LC₅₀?=?7.07?mg/L air). Moreover, two EOs exhibited over 80% repellent rate against T. castaneum at the highest concentration of 78.63?nL/cm². Furthermore, hundreds of enzymes related to the regulation of biological processes of T. castaneum were screened to explore the underlying molecular mechanism and develop promising insecticides. Besides, top hits were subjected to MD simulations and binding free energies analysis to elucidate complex inter-molecular stability and affinity over simulated time. The results demonstrated that isolongifolene, δ-cadinene, β-caryophyllene and caryophyllene oxide were prioritized as they were establishing conserved and stable interactions with residues of nuclear hormone receptor 3 (TcHR3) of T. castaneum, which suggested that the four sesquiterpenes have potential to be promising insecticides on storage pests control.