More environment-friendly solvents in analytical chemistry are almost inevitable in both sample preparation and analysis to meet the high standards of sustainable chemistry and human welfare.

Commercially available pesticides were examined as Mus musculus and Homo sapiens acetylcholinesterase (mAChE and hAChE) inhibitors by means of ligand-based (LB) and structure-based (SB) in silico approaches. Initially, the crystal structures of simazine, monocrotophos, dimethoate, and acetamiprid were reproduced using various force fields. Subsequently, LB alignment rules were assessed and applied to determine the inter synaptic conformations of atrazine, propazine, carbofuran, carbaryl, tebufenozide, imidacloprid, diuron, monuron, and linuron. Afterwards, molecular docking and dynamics SB studies were performed on either mAChE or hAChE, to predict the listed pesticides’ binding modes. Calculated energies of global minima (E_{glob_min}) and free energies of binding (∆G_binding) were correlated with the pesticides’ acute toxicities (i.e., the LD₅₀ values) against mice, as well to generate the model that could predict the LD₅₀s against humans. Although for most of the pesticides the low E_{glob_min} correlates with the high acute toxicity, it is the ∆G_binding that conditions the LD₅₀ values for all the evaluated pesticides. Derived pLD₅₀ = f(∆G_binding) mAChE model may predict the pLD₅₀ against hAChE, too. The hAChE inhibition by atrazine, propazine, and simazine (the most toxic pesticides) was elucidated by SB quantum mechanics (QM) DFT mechanistic and concentration-dependent kinetic studies, enriching the knowledge for design of less toxic pesticides.