Homology modeling of the serotonin transporter: insights into the primary escitalopram-binding site.

The serotonin transporter (SERT) is one of the neurotransmitter transporters that plays a critical role in the regulation of endogenous amine concentrations and therefore is an important target for therapeutic agents affecting the central nervous system. The recently published, high resolution X-ray structure of the closely related amino acid transporter, Aquifex aeolicus leucine transporter (LeuT), provides an opportunity to develop a three-dimensional model of the structure of SERT. We present herein a homology model of SERT using LeuT as the template and containing escitalopram as a bound ligand. Our model explains selectivities known from mutational studies and varying ligand data, which are discussed and illustrated in the paper.

Probing the pharmacophore of ginkgolides as glycine receptor antagonists.

Ginkgolides are antagonists of the inhibitory ligand-gated ion channels for the neurotransmitters glycine and gamma-aminobutyric acid (GABA). In this study the ginkgolide structure was modified in order to investigate the minimum structural requirements for glycine receptor antagonism. The five native ginkgolides and a series of 29 ginkgolide derivatives were characterized at the three glycine receptor subtypes alpha1, alpha1beta, and alpha2, which revealed that only minor changes in the ginkgolide skeleton were allowed for maintaining glycine receptor antagonism. A pharmacophore model was generated and applied in a virtual screening of a compound database (300000 compounds), resulting in the identification of 31 hits. Twenty-seven of these hits were screened for biological activity, but none displayed antagonist activity at the glycine receptors. This strongly suggests the importance of other pharmacophore components in the binding of ginkgolides to glycine receptors, and we propose that the structural rigidity of the ginkgolide molecule may be crucial for its glycine receptor activity.

A 3D QSAR study on a set of dopamine D4 receptor antagonists.

The molecular alignments obtained from a previously reported pharmacophore model have been employed in a three-dimensional quantitative structure-activity relationship (3D QSAR) study, to obtain a more detailed insight into the structure-activity relationships for D(2) and D(4) receptor antagonists. The frequently applied CoMFA method and the related CoMSIA method were used. Statistically significant models have been derived with these two methods, based on a set of 32 structurally diverse D(2) and D(4) receptor antagonists. The CoMSIA and the CoMFA methods produced equally good models expressed in terms of q(2) values. The predictive power of the derived models were demonstrated to be high. Graphical interpretation of the results, provided by the CoMSIA method, brings to light important structural features of the compounds related to either low- or high-affinity D(2) or D(4) antagonism. The results of the 3D QSAR studies indicate that bulky N-substituents decrease D(2) binding, whereas D(4) binding is enhanced. Electrostatically favorable and unfavorable regions exclusive to D(2) receptor binding were identified. Likewise, certain hydrogen-bond acceptors can be used to lower D(2) affinity. These observations may be exploited for the design of novel dopamine D(4) selective antagonists.

Pharmacophore and receptor models for neurokinin receptors.

Three neurokinin (NK) antagonist pharmacophore models (Models 1-3) accounting for hydrogen bonding groups in the 'head' and 'tail' of NK receptor ligands have been developed by use of a new procedure for treatment of hydrogen bonds during superimposition. Instead of modelling the hydrogen bond acceptor vector in the strict direction of the lone pair, an angle is allowed between the hydrogen bond acceptor direction and the ideal lone pair direction. This approach adds flexibility to hydrogen bond directions and produces more realistic RMS values. By using this approach, two novel pharmacophore models were derived (Models 2 and 3) and a hydrogen bond acceptor was added to a previously published NK2 pharmacophore model [Poulsen et al., J. Comput.-Aided Mol. Design, 16 (2002) 273] (Model 1). Model 2 as well as Model 3 are described by seven pharmacophore elements: three hydrophobic groups, three hydrogen bond acceptors and a hydrogen bond donor. Model 1 contains the same hydrophobic groups and hydrogen bond donor as Models 2 and 3, but only one hydrogen bond acceptor. The hydrogen bond acceptors and donor are represented as vectors. Two of the hydrophobic groups are always aromatic rings whereas the other hydrophobic group can be either aromatic or aliphatic. In Model 1 the antagonists bind in an extended conformation with two aromatic rings in a parallel displaced and tilted conformation. Model 2 has the same two aromatic rings in a parallel displaced conformation whereas Model 3 has the rings in an edge to face conformation. The pharmacophore models were evaluated using both a structure (NK receptor homology models) and a ligand based approach. By use of exhaustive conformational analysis (MMFFs force field and the GB/SA hydration model) and least-squares molecular superimposition studies, 21 non-peptide antagonists from several structurally diverse classes were fitted to the pharmacophore models. More antagonists could be fitted to Model 2 with a low RMS and a low conformational energy penalty than to Models 1 and 3. Pharmacophore Model 2 was also able to explain the NK1, NK2 and NK3 subtype selectivity of the compounds fitted to the model. Three NK 7TM receptor models were constructed, one for each receptor subtype. The location of the antagonist binding site in the three NK receptor models is identical. Compounds fitted to pharmacophore Model 2 could be docked into the NK1, NK2 and NK3 receptor models after adjustment of the conformation of the flexible linker connecting the head and tail. Models I and 3 are not compatible with the receptor models.

A pharmacophore model for NK2 antagonist comprising compounds from several structurally diverse classes.

A neurokinin 2 (NK2) antagonist pharmacophore model has been developed on the basis of five non-peptide antagonists from several structurally diverse classes. To evaluate the pharmacophore model, another 20 antagonists were fitted to the model. By use of exhaustive conformational analysis (MMFFs force field and the GB/SA hydration model) and least-squares molecular superimposition studies, 23 of the 25 antagonists were fitted to the model in a low energy conformation with a low RMS value. The pharmacophore model is described by four pharmacophore elements: Three hydrophobic groups and a hydrogen bond donor represented as a vector. The hydrophobic groups are generally aromatic rings, but this is not a requirement. The antagonists bind in an extended conformation with two aromatic rings in a parallel displaced and tilted conformation. The model was able to explain the enantioselectivity of SR48968 and GR159897.