Derivation of a three-dimensional pharmacophore model of substance P antagonists bound to the neurokinin-1 receptor

Abstract

Constrained systematic search was used in an exhaustive conformational analysis of a structurally diverse set of substance P (SP) antagonists to identify a unique hypothesis for their bound conformation at the neurokinin-1 receptor. In this conformation, two aromatic groups essential for high affinity adopt a perpendicular or edge-on arrangement. This pharmacophore hypothesis for the receptor-bound conformation was used in a comparative molecular field analysis (CoMFA) of an expanded set of SP antagonists, and the predictive ability of the resulting three-dimensional quantitative structure-activity relationship (3D-QSAR) was evaluated against a test set of SP antagonists different from those in the training set. This CoMFA model based on the Constrained Search alignment yielded significant cross-validated, conventional, and predictive r2 values equal to 0.70, 0.93, and 0.82, respectively. For comparison, the SP antagonists were forced into an alternative poorer alignment in which the two aromatic rings were parallel and then subjected to a CoMFA analysis. Both the parallel and perpendicular arrangements of the aromatic rings are seen in X-ray structures of SP antagonists and have been proposed as candidates for the receptor-bound conformation. The parallel (or stacked) conformation yielded a poorer correlation with a cross-validated r2 = 0.57, a conventional r2 = 0.90, and a predictive r2 = 0.78. Our results indicate that although both alignments could generate a reasonable CoMFA correlation, the stacked conformation is unlikely to be the receptor-bound conformation, as the covalent structure of the antagonists precludes a common geometry in which the aromatic rings are stacked.