Molecular modeling of sigma receptor ligands: A model of binding based on conformational and electrostatic considerations

We have performed molecular modeling studies on four representative sigma receptor specific ligands, (+)haloperidol, (+)3-PPP, (+)pentazocine and progesterone, to develop a model for sigma receptor-ligand binding. The modeling studies have investigated the conformational and electrostatic properties of the ligands. Based on the complementarity of the conformational and electrostatic properties of the ligands, a model of binding has been proposed which shows that the four ligands can fit a common receptor sit. Unlike the binding model for haloperidol that was previously proposed by Manallack and Andrews, our model binds haloperidol in the gauche conformation. The first site binds the fluorophenyl group and the second site the lone pair of the piperidine nitrogen. This pharmacophore can be presented by (+)3-PPP and (+)pentazocine, but for progesterone the binding model requires the ring junction of the cyclohexenyl ring A and ring B to fit the fluorophenyl region, while the lone pair of the acetylcarbonyl oxygen at ring D emulates the nitrogen lone pair of the piperidine ring. Calculations were performed using RCG5 for generating conformations, molecular mechanics for calculating steric energies, quantum mechanical methods for generating charges, and ARCHEM for calculating electrostatic potentials on the Van der Waals surface.