in commonly accepted opinion every searching for new , more effective drugs should be rationalized i.e. , determined by the low cost and non time - consuming procedures . these procedures are especially useful on the preliminary stage of searching for new chemical structures of potential biological activity ( jorgensen , 2004 ; leeson and springthorpe , 2007 ; ou - yang et al . , 2012 ) . in general , on this purpose there are employed various correlation qsar methods ( dudek et al . , 2006 ; yang and huang , 2006 ; shailesh et al . , 2012 ) . however , in particular cases it is more convenient to develop the procedure of selection of the appropriate structures based on more direct and easier interpretatively criteria . it seems that just such a case is a search for effective ligands of 5ht1a , 5ht2a , and d2 receptors since many structural data on their agonist and antagonist as well as the models of these receptors are well - known ( klabunde and hessler , 2002 ; bissantz et al . , 2003 ; teeter et al . , 1994 ; chambers and nichols , 2002 ; homan et al . , 1999 ) . in addition , wide availability of various bases containing a lot of structural data on very active ligands allows to generate pretty accurate pharmacophore patterns ( nelson , 1991 ; bojarski , 2006 ) . thanks to these all literature data it is possible to estimate the affinity of potential ligand for receptor of interest . the chemical structure of pharmacophore of being selected potential ligand and its affinity to the receptor seem to be sufficiently unambiguous discriminators , on a preliminary stage , in the search for new effective antipsychotics . to verify this hypothesis , the first step includes determination of pharmacophores for two tested compounds of well - known affinity ( previously in vitro determined ) to the same receptors as well as pharmacophore pertinent to well - known d2 receptor agonists or antagonists and finally comparison of their properties to in vitro binding data . the pharmacophore model of d2 receptor ligands was found on the basis of 15 compounds of high affinity to d2 receptor reported in literature ( sowiski et al . , 2011 ) . these two tested compounds were 3-acylamine derivatives of tropane : n-(8-furan-2-ylmethyl-8-azabicyclo[3.2.1]oct-3-yl)-2-methoxybenzamide ( compounds i ) and n-(8-furan-2-ylmethyl-8-azabicyclo[3.2.1]oct-3-yl)-2.3-dimethoxybenzamide ( compound ii ) ( fig .  1 ) . their synthesis have been developed and described in the previously published paper on tropane derivatives ( sowiski et al . .  1the chemical formulas of compound i and compound ii the chemical formulas of compound i and compound ii the pharmacophores of compounds i and ii were found on the basis of their structures determined by x - ray diffraction method . the ccdc ( cambridge crystallographic data centre ) numbers of compounds i and ii are : 905689 and 905690 , respectively ( figs .  2 , 3).fig . 3the x - ray diffraction structure of compound ii the x - ray diffraction structure of compound i the x - ray diffraction structure of compound ii the molecular structure of compound i shows an intramolecular hydrogen bond between the o atom of the methoxy group and the nh of the amide function leads to a six - membered ring . the dihedral angle between the least - squares planes of the phenyl and this virtual ring is only 2.50(7). the piperidine moiety adopts a chair conformation . the best plane of the furan ring and the c1/c2/c4/c5 plane make an angle 69.42(9) and the dihedral angle between the planes of the furan and benzene rings is 72.50(8). the compound ii molecule adopts a folded conformation with an angle between the furan and benzene rings of 63.29(8) and between the best plane of the furan ring and the c1/c2/c4/c5 plane of 87.56(9). this conformation is stabilized by an intramolecular n15h15ao25 and c26h26co27 hydrogen bonds . as a result of n15h15ao25 interaction a six - membered ring the piperidine moiety assumes a chair conformation and the substituent at n8 is in an equatorial position . conformations of both methoxy groups are different . the disposition of these groups with respect to the phenyl ring can be described by the torsion angles c18c19o25c26 of 107.8(2) and c21c20o27c28 of 11.1(3). in consequence , the methyl carbon atom c26 is found to be 1.107(4 )   out of the phenyl plane , and c28 atom is almost coplanar with this ring . the pharmacophore structure is a reflection template of the geometrical distribution of property centers localized in molecule and determines to large extent its biological activity . it means that even subtle differences in the geometry of structurally similar molecules can significantly impact on their affinity to receptor binding site . the comparative analysis of the studied pharmacophores was intended to find the specific properties and geometrical parameters which are crucial for the strength of binding of potential ligands to the receptors of interest . the second step of the applied procedure devoted to the selection of the potential agonists or antagonists of the studied receptors relies on docking of the reference compounds i and ii to the models of the d2 receptor ( sakhteman et al . , 2011 ) . from analysis of in vitro results ( table  1 ) follows that the both studied compounds ( i , ii ) are very poorly being bounded to 5-ht1a and 5-ht2a receptors . indeed , the model docking of compounds i and ii to these receptors also showed that such binding can not take place . the both molecules of compounds i and ii were placed outside the receptor binding pockets . thus , only docking of compounds i and ii to d2 receptor is detailed analyzed . the most discriminative parameters which distinctly classify the quality of docking are number and strength ( equivalently length and geometry ) of the hydrogen bonds formed between ligand and specific amino acids not only inside the receptor binding pocket but also , although to a less degree , intermolecular interactions of other types e.g. , hydrophobic and edge-to-face.table  15ht1a , 5ht2a , and d2 receptor affinitiesligandreceptor [ k(nm)]5ht1a5ht2ad2compound i6,1006,0001,000compound ii3,000744.526.3 5ht1a , 5ht2a , and d2 receptor affinities the used 3d homology model of d2 receptor has been revealed by comparative modeling using the crystal structure of the human 2-adrenergic receptor and the bovine rhodopsin as the templates ( sakhteman et al . , 2011 ; strzelczyk et al . , 2004 ; the quite recently reported x - ray structure of the human 2-adrenergic receptor opens new possibilities for modeling of the correct structures of the dopamine ones . currently , the human 2-adrenergic receptor is considered to be more homologous to the dopamine receptors than bovine rhodopsin ( cherezov et al . , 2007 ) . all modeling of the pharmacophores as well as docking of the compounds i and ii to the d2 receptor model were done by discovery studio software ( accelrys software inc . , discovery studio modeling environment , 2005 ) . crystals of compounds i and ii suitable for x - ray analysis were grown by slow evaporation from acetate / diisopropyl ether ( compound i ) and hexane / ethanol ( compound ii ) solutions . the data were collected on an oxford diffraction km4ccd diffractometer at 293  k , using graphite - monochromated mo k radiation . the unit cell parameters were determined by least - squares treatment of setting angles of highest - intensity reflections chosen from the whole experiment . the structure was solved by direct methods using the shelxs97 program ( sheldric , 1990 ) and refined by the full - matrix least - squares method with the shelxl97 program ( sheldric , 1997 ) . the function w(|fo| |fc| ) was minimized with w=  [ (fo )  +  ( 0.0688p ) ] , where p  =  ( fo2  +  2fc2)/3 . =  kfc[1  +  ( 0.001fc2/sin2 ) ] ( sheldric , 1997 ) and the extinction coefficient  was equal to 0.014(2 ) . the coordinates of the hydrogen atoms were calculated in idealized positions and refined as a riding model with their thermal parameters calculated as 1.2 ( 1.5 for methyl group ) times ueq of the respective carrier carbon atom . crystals of compounds i and ii suitable for x - ray analysis were grown by slow evaporation from acetate / diisopropyl ether ( compound i ) and hexane / ethanol ( compound ii ) solutions . the data were collected on an oxford diffraction km4ccd diffractometer at 293  k , using graphite - monochromated mo k radiation . the unit cell parameters were determined by least - squares treatment of setting angles of highest - intensity reflections chosen from the whole experiment . the structure was solved by direct methods using the shelxs97 program ( sheldric , 1990 ) and refined by the full - matrix least - squares method with the shelxl97 program ( sheldric , 1997 ) . the function w(|fo| |fc| ) was minimized with w=  [ (fo )  +  ( 0.0688p ) ] , where p  =  ( fo2  +  2fc2)/3 . =  kfc[1  +  ( 0.001fc2/sin2 ) ] ( sheldric , 1997 ) and the extinction coefficient  was equal to 0.014(2 ) . the coordinates of the hydrogen atoms were calculated in idealized positions and refined as a riding model with their thermal parameters calculated as 1.2 ( 1.5 for methyl group ) times ueq of the respective carrier carbon atom . the in vitro binding data for compounds i , ii as ligands of 5ht1a , 5ht2a , and d2 receptors are given in table  1 ( sowiski et al . , 2011 ) . these experimental binding data unambiguously points at very low affinity of compound i to 5ht1a and 5ht2a receptors and somewhat better to d2 one , yet , compound ii displayed very weak binding activity to 5ht1a , moderate to 5ht2a and very high to d2 receptors . the differences between parameters ( geometrical and property types ) of the reference pharmacophores and the pharmacophores pertinent to compounds i and ii are expected to reflect the differences in affinity of tested compounds to the receptors of interest . the found structures of pharmacophores described by their specific properties are given on  figs . the particular colors denote the following properties : red  positive ionization ( nitrogen atom ) , green  hydrogen bond acceptor , magenta  hydrogen bond donor , pale blue  hydrophobic , aromatic , ultramarine  hydrophobic , aliphatic , orange  aromatic ring ( table  2).fig . 4the spatial distribution of pharmacophore properties on a background of compound i x - ray diffraction structure . a green square depicts the plane of a phenyl ring ( color figure online)fig . 5the spatial distribution of pharmacophore properties on a background of compound ii x - ray diffraction structure . a green square depicts the plane of a phenyl ring ( color figure online)fig . aliphatic property ( color figure online)table  2pharmacophore properties of compound i and iipharmacophore feature / propertycompound icompound iipositive ionization ( red)nitrogen atomnitrogen atomhydrogen bond acceptor ( hba , green)carbonyl group of amide bondcarbonyl group of amide bondaromatic ring ( orange)benzene ring substituted with methoxy groupbenzene ring substituted with two methoxy groupshydrophobic , aromatic ( pale blue)furane ringfurane ringhydrophobic , aliphatic ( ultramarine)one methyl group in methoxy moiety attached to the benzene ringtwo methyl groups in methoxy moieties attached to the benzene ring the spatial distribution of pharmacophore properties on a background of compound i x - ray diffraction structure . a green square depicts the plane of a phenyl ring ( color figure online ) the spatial distribution of pharmacophore properties on a background of compound ii x - ray diffraction structure . a green square depicts the plane of a phenyl ring ( color figure online ) the spatial distribution of pharmacophore properties of d2 receptor ligands . aliphatic property ( color figure online ) pharmacophore properties of compound i and ii the geometry of a spatial distribution of pharmacophore properties in obtained models is an exact reflection of the x - ray diffraction structure of compounds i and ii ( table  3 ) . it is worthy to note that in spite of the high similarity of chemical structures of these compounds , that their conformations significantly differ each from other . obviously , it should be taken into account some flexibility of the spatial pharmacophore geometry and possibility of its change during docking of studied compounds to particular receptors . however , such changes are often possible only to small degree or impossible at all on account of the high energetic rotation barriers . in this context , the presence of two separate aliphatic  hydrophobic centers in pharmacophore of compound ii takes on a special importance for explanation of very high affinity of this compound , in contrast to compound i , for d2 receptor . it is likely that just second methoxy group in compound ii molecule underlies its high binding to d2 receptor while the same group do not affect the affinity of compound ii to 5-ht1a and 5-ht2a receptors . the comparative analysis of the d2 receptor ligand pharmacophore ( fig .  6 ) and pharmacophores of compounds i and ii also leads to the same conclusion ( figs .  4 and 5 ) . the pharmacophore of d2 ligand quite well matches the pharmacophore of compound ii but does not the pharmacophore of compound i ( c.f . in addition , specificity of the structural relation between these pharmacophores results from the identical spatial localization of the aliphatic property of d2 ligand pharmacophore and its analog present in pharmacophore of compound ii but absent in pharmacophore of compound i ( table  2 ) . 
 fig . 7superposition of the d2 receptor ligand pharmacophore and pharmacophore of compound iitable  3pharmacophore geometry parameterspharmacophore geometry parameterscompound icompound iidistance between piperidine nitrogen atom and center of the benzene ring7.85  7.76  dihedral angle between benzene ring plane and furane ring plane72.5063.29dihedral angle between piperidine ring ( c1/c2/c4/c5 ) plane and benzene ring plane65.7950.97dihedral angle between piperidine ring ( c1/c2/c4/c5 ) plane and furane ring plane69.4287.56dihedral angle between carbonyl group plane and piperidine ring plane73.5086.72 superposition of the d2 receptor ligand pharmacophore and pharmacophore of compound ii pharmacophore geometry parameters docking of both tested compounds to d2 receptor model turned out to be non discriminative investigation not giving criteria for explanation of difference in ability to the binding of compounds i and ii with d2 receptor . both compounds docked to d2 receptor interact with its amino acids via the same hydrogen bonds . in case of compound i the hydrogen bonds are : ligand  thyrosine 379 ( length 2.198   ) , ligand  alanine 185 ( length 2.315   ) , and compound ii ligand thyrosine 379 ( length 2.310   ) , ligand  alanine 185 ( length 2.139   ) . in addition , both compounds interact similarly with d2 receptor with hydrophobic forces ( fig . 8the molecules of compounds i and ii ( green ) inside binding pocket of d2 receptor . yellow dashed lines denote hydrogen bonds ( color figure online ) the molecules of compounds i and ii ( green ) inside binding pocket of d2 receptor . yellow dashed lines denote hydrogen bonds ( color figure online ) the obtained docking results are not unexpected since , purposely , the structurally similar compounds were investigated to point out that even very subtle differences in the chemical structure of compounds , to which docking procedure is  insensitive  , may impact crucially on their therapeutic activity . thus , it should be stated that two stages  pharmacophore  and  docking  investigations are necessary to estimate properly an affinity of newly designed receptor ligands . on the whole , these studies were intended to prove that postulated two - stages procedure can be applied to verification of the properties of even very similar structurally potential and being designed antipsychotics .
the two - stages studies of structure  activity relationship for model ligands of 5ht1a , 5ht2a , and d2 receptors were performed . on the first stage , the pharmacophores of two potential ligands of known in vitro binding to 5ht1a , 5ht2a , d2 receptors and model pharmacophore of strongly interacting d2 receptor ligands were found and their parameters were related to affinity data . 
 the analyzed parameters were hydrophobic , hydrophilic , aromatic , donor and acceptor of proton centers . 
 the geometry of spatial distribution of these properties was also investigated in comparative analysis . 
 the studied , model compounds were two 3-acylamine derivatives of tropane . 
 the second stage includes docking of studied compounds to d2 receptor model and the comparison of its quality with in vivo binding data . 
 the obtained results are consistent with in vitro binding data and applied procedure accurate estimates the affinity of potential ligands to d2 receptors .