This paper describes a rhodopsin-based model of 5-HT(1A) serotonin receptor. The flexibility of the receptor was considered by using large number of models for ligand dockings. Rearrangements of the heptahelical bundle were introduced, which resulted in the improvement of correlation between computational results and experimental data. The model was validated by automated docking of conformationally restricted arylpiperazines. Specific interactions, responsible for the recognition of arylpiperazine derivatives, were identified. An ionic bond was formed between the protonated amine of ligands and Asp3.32. The aromatic moiety and its substituents specifically interacted with Phe6.52 and Ser5.42, respectively, while the carbonyl groups of imide part of ligands formed hydrogen bonds with Asn7.39 and Tyr7.43. The model reproduced the binding affinity of the test group of ligands (correlation r = 0.8 between pK(i) and docking score). It also gave the enrichment in virtual screening-like experiment (100 compounds), in which 34 high-affinity compounds were found among 50 top-scored ligands.

A set of 31 diversified 5-HT7 receptor antagonists was automatically docked to a conformational ensemble of rhodopsin-based 5-HT7R models (flexible docking). It was found that ergolines, aporphines, and tricyclic psychotropic agents were always docked in a pocket formed by TMHs 4-6, and besides the main ionic bond with Asp3.32, they had specific interactions with Phe6.52, Phe6.51, Trp6.48, and Ser5.42. The arylpiperidine, arylpiperazine, or beta-carboline fragment of the complex ligands occupied the same pocket, whereas the terminal amide/imide part of those compounds reached the second cavity formed by TMHs 7-3 and interacted with Phe3.28, Arg7.36, and Tyr7.43. A similar orientation of selective antagonists of the group of arylsulfonamidoalkylamines was observed, that is, the sulfonamide part was located in the second pocket. Coherent docking results allowed the generation of two receptor-based pharmacophores: first containing features necessary for high 5-HT7R affinity and the other defining selectivity for this receptor subtype. The latter model indicated the importance of specific interactions with the residues of the TMHs 7-3 pocket (especially nonconserved Arg7.36) for selectivity over other monoamine GPCRs.