Three-dimensional (3-D) models of the human serotonin 5-HT1A and 5-HT2A receptors were constructed, energy refined, and used to study the interactions with a series of buspirone analogues. For both receptors, the calculations showed that the main interactions of the ligand imide moieties were with amino acids in transmembrane helix (TMH) 2 and 7, while the main interactions of the ligand aromatic moieties were with amino acids in TMH5, 6 and 7. Differences in binding site architecture in the region of highly conserved serine and tyrosine residues in TMH7 gave slightly different binding modes of the buspirone analogues at the 5-HT1A and 5-HT2A receptors. Molecular dynamics simulations of receptor-ligand interactions indicated that the buspirone analogues did not alter the interhelical hydrogen bonding patterns upon binding to the 5-HT2A receptor, while interhelical hydrogen bonds were broken and others were formed upon ligand binding to the 5-HT1A receptor. The ligand-induced changes in interhelical hydrogen bonding patterns of the 5-HT1A receptor were followed by rigid body movements of TMH2, 4 and 6 relative to each other and to the other TMHs, which may reflect the structural conversion into an active receptor structure.