A nonlinear QSAR study was conducted on a series of 4-phenylpiperidine derivatives (4PPs) acting as mu opioid agonists by three-layer back-propagation neural network (NN) method. At first a variety of molecular descriptors were calculated and then selected with two-stage least squares combining partial least squares (PLS) method. The selected four molecular descriptors, out of 292 ones, were correlated with the known analgesic activities of 38 4PPs by NN method. The established QSAR model was further validated by five additional 4PPs, as an external testing set. Moreover, a pharmacophore model was hypothesized based on the results, which would be helpful for structural optimization of 4PPs.