Phenoloxidase (PO), also known as tyrosinase, is a key enzyme in insect development, responsible for catalyzing the hydroxylation of tyrosine into o-diphenols and the oxidation of o-diphenols into o-quinones. Inhibition of PO may provide a basis for novel environmentally friendly insecticides. In the present study, we determined the inhibitory activities and IC50 values of 57 compounds belonging to the benzaldehyde thiosemicarbazone, benzaldehyde, and benzoic acid families against phenoloxidase from Pieris rapae (Lepidoptera) larvae. In addition, the inhibitory kinetics of 4-butylbenzaldehyde thiosemicarbazone against PO was measured in air-saturated solutions for the oxidation of L-3,4-dihydroxyphenylalanine (L-DOPA). The results indicated that the compound is a reversible noncompetitive inhibitor. The bioactivity results were used to construct three-dimensional quantitative structure-activity relationship (3D-QSAR) models using two molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). After carrying out superimposition using common substructure-based alignment, robust and predictive 3D-QSAR models were obtained from CoMFA (q2=0.926, r2=0.986) and CoMSIA (q2=0.933, r2=0.984) with six optimum components. The 3D-QSAR model built here will provide hints for the design of novel PO inhibitors. The molecular interactions between the ligands and the target were studied using a flexible docking method (FlexX). The best scored candidates were docked flexibly, and the interaction between the representative compound 4-butylbenzaldehyde thiosemicarbazone and the active site was elucidated in detail.