Homology Modeling of the Serotonin Transporter: Insights into the Primary Escitalopram-binding Site

Substrate and inhibitor binding to the serotonin transporter: Insights from computational, crystallographic, and functional studies

The serotonin transporter (SERT) belongs to the monoamine transporter family, which also includes the dopamine and norepinephrine transporters. SERT is essential for regulating serotonergic signaling by the reuptake of serotonin from the synaptic cleft back into the presynaptic neuron. Dysregulation of SERT has been implicated in several major psychiatric disorders such as major depressive disorder (MDD). MDD was among the top five leading causes of years lived with disease in 2016 and is characterized as a major global burden. Several drugs have been developed to target SERT for use in the treatment of MDD, and their respective binding modes and locations within SERT have been studied. The elucidation of the first structure of a bacterial SERT homologue in 2005 has accelerated crystallographic, computational, and functional studies to further elucidate drug binding and method of action in SERT. Herein, we aim to highlight and compare these studies with an emphasis on what the different experimental methods conclude on substrate and inhibitor binding modes, and the potential caveats of using the different types of studies are discussed. We focus this review on the binding of cognate substrate and drugs belonging to the different families of antidepressants, including tricyclic antidepressants, selective serotonin reuptake inhibitors, serotonin-norepinephrine reuptake inhibitors, and multimodal drugs, as well as illicit drugs such as cocaine, amphetamines, and ibogaine. This article is part of the issue entitled 'Special Issue on Neurotransmitter Transporters'.

Characterization of a Novel Drosophila SERT Mutant: Insights on the Contribution of the Serotonin Neural System to Behaviors

A better comprehension on how different molecular components of the serotonergic system contribute to the adequate regulation of behaviors in animals is essential in the interpretation on how they are involved in neuropsychiatric and pathological disorders. It is possible to study these components in "simpler" animal models including the fly Drosophila melanogaster, given that most of the components of the serotonergic system are conserved between vertebrates and invertebrates. Here we decided to advance our understanding on how the serotonin plasma membrane transporter (SERT) contributes to serotonergic neurotransmission and behaviors in Drosophila. In doing this, we characterized for the first time a mutant for Drosophila SERT (dSERT) and additionally used a highly selective serotonin-releasing drug, 4-methylthioamphetamine (4-MTA), whose mechanism of action involves the SERT protein. Our results show that dSERT mutant animals exhibit an increased survival rate in stress conditions, increased basal motor behavior, and decreased levels in an anxiety-related parameter, centrophobism. We also show that 4-MTA increases the negative chemotaxis toward a strong aversive odorant, benzaldehyde. Our neurochemical data suggest that this effect is mediated by dSERT and depends on the 4-MTA-increased release of serotonin in the fly brain. Our in silico data support the idea that these effects are explained by specific interactions between 4-MTA and dSERT. In sum, our neurochemical, in silico, and behavioral analyses demonstrate the critical importance of the serotonergic system and particularly dSERT functioning in modulating several behaviors in Drosophila.

Molecular dynamics of conformation-specific dopamine transporter-inhibitor complexes

The recreational psychostimulant cocaine inhibits dopamine reuptake from the synapse, resulting in excessive stimulation of postsynaptic dopamine receptors in brain areas associated with reward and addiction. Cocaine binds to and stabilizes the outward- (extracellular-) facing conformation of the dopamine transporter (DAT) protein, while the low abuse potential DAT inhibitor benztropine prefers the inward- (cytoplasmic-) facing conformation. A correlation has been previously postulated between psychostimulant abuse potential and preference for the outward-facing DAT conformation. The 3β-aryltropane cocaine analogs LX10 and LX11, however, differ only in stereochemistry and share a preference for the outward-facing DAT, yet are reported to vary widely in abuse potential in an animal model. In search of the molecular basis for DAT conformation preference, complexes of cocaine, benztropine, LX10 or LX11 bound to each DAT conformation were subjected to 100ns of all-atom molecular dynamics simulation. Results were consistent with previous findings from cysteine accessibility assays used to assess an inhibitor's DAT conformation preference. The respective 2β- and 2α-substituted phenyltropanes of LX10 and LX11 interacted with hydrophobic regions of the DAT S1 binding site that were inaccessible to cocaine. Solvent accessibility measurements also revealed subtle differences in inhibitor positioning within a given DAT conformation. This work serves to advance our understanding of the conformational selectivity of DAT inhibitors and suggests that MD may be useful in antipsychostimulant therapeutic design.

Exploration of insights, opportunities and caveats provided by the X-ray structures of hSERT

The recently reported X-ray structures of the human serotonin (5-HT) transporter SERT with bound inhibitors open new opportunities for drug discovery at SERT, selectivity design with respect to other neurotransmitter sodium transporters, and enhanced understanding of the molecular events involved in SERT action. Through computational and structural analyses, we explore the binding and migration of 5-HT at SERT. Consistent with earlier studies of leucine migration at the bacterial homolog of SERT, LeuT, we find multiple potential 'stopover' sites for 5-HT binding at SERT including the two (transmembrane S1 and extracellular vestibule S2) seen in the binding of the SSRI (S)-citalopram (S-Cit) to SERT, as well as other sites. Docking studies reveal the possibility of both hetero- (S-Cit+5-HT) and homo-dimeric (5-HT+5-HT) co-binding at both these sites which may explain earlier published allosteric activity observations and provide novel design strategies. Comparisons with substrate bound X-ray structures of the dopamine transporter reveal a number of potential sources of selectivity, some of which may be 'artificial' including target based, species related, experimental design related, and ligand dependent examples including substrate versus inhibitor related features.

Insights to ligand binding to the monoamine transporters-from homology modeling to LeuBAT and dDAT

Understanding of drug binding to the human biogenic amine transporters (BATs) is essential to explain the mechanism of action of these pharmaceuticals but more importantly to be able to develop new and improved compounds to be used in the treatment of depression or drug addiction. Until recently no high resolution structure was available of the BATs and homology modeling was a necessity. Various studies have revealed experimentally validated binding modes of numerous ligands to the BATs using homology modeling. Here we examine and discuss the similarities between the binding models of substrates, antidepressants, psychostimulants, and mazindol in homology models of the human BATs and the recently published crystal structures of the Drosophila dopamine transporter and the engineered protein, LeuBAT. The comparison reveals that careful computational modeling combined with experimental data can be utilized to predict binding of molecules to proteins that agree very well with crystal structures.

Structural basis for action by diverse antidepressants on biogenic amine transporters

The biogenic amine transporters (BATs) regulate endogenous neurotransmitter concentrations and are targets for a broad range of therapeutic agents including selective serotonin reuptake inhibitors (SSRIs), serotonin-noradrenaline reuptake inhibitors (SNRIs) and tricyclic antidepressants (TCAs). Because eukaryotic BATs are recalcitrant to crystallographic analysis, our understanding of the mechanism of these inhibitors and antidepressants is limited. LeuT is a bacterial homologue of BATs and has proven to be a valuable paradigm for understanding relationships between their structure and function. However, because only approximately 25% of the amino acid sequence of LeuT is in common with that of BATs, and as LeuT is a promiscuous amino acid transporter, it does not recapitulate the pharmacological properties of BATs. Indeed, SSRIs and TCAs bind in the extracellular vestibule of LeuT and act as non-competitive inhibitors of transport. By contrast, multiple studies demonstrate that both TCAs and SSRIs are competitive inhibitors for eukaryotic BATs and bind to the primary binding pocket. Here we engineered LeuT to harbour human BAT-like pharmacology by mutating key residues around the primary binding pocket. The final LeuBAT mutant binds the SSRI sertraline with a binding constant of 18 nM and displays high-affinity binding to a range of SSRIs, SNRIs and a TCA. We determined 12 crystal structures of LeuBAT in complex with four classes of antidepressants. The chemically diverse inhibitors have a remarkably similar mode of binding in which they straddle transmembrane helix (TM) 3, wedge between TM3/TM8 and TM1/TM6, and lock the transporter in a sodium- and chloride-bound outward-facing open conformation. Together, these studies define common and simple principles for the action of SSRIs, SNRIs and TCAs on BATs.

New design strategies for antidepressant drugs

INTRODUCTION

In spite of research efforts spanning six decades, the most prominent antidepressant drugs to date still carry several adverse effects, often serious enough to warrant discontinuation of the drug. Molecular mechanisms of depression are now better understood such that some of the specific receptors responsible can be targeted for activation or inhibition. This advance, coupled with the recent availability of crystal structures of relevant drug targets or their homologs, has opened the door for new antidepressant therapeutic compounds.

AREAS COVERED

The authors review the evolution of monoamine-based antidepressant drugs, up to the selective serotonin reuptake inhibitors (SSRIs). The authors discuss classic and contemporary antidepressant drug design strategies, with a focus on virtual screening and fragment-based drug design methods. Furthermore, they discuss the recent advancements in the understanding of the serotonin transporter (SERT) structure/function relationship in the context of recognition of SSRIs and outline a strategy for the use of computational approaches in producing new SSRI lead compounds.

EXPERT OPINION

The authors suggest that given the long-awaited availability of credible three-dimensional structures for the SERT and related monoamine transporter proteins, cutting-edge computational methods should be the linchpin of future drug discovery efforts regarding monoamine-based antidepressant lead compounds. Because these transporter inhibitors cause a ubiquitous increase in extraneuronal neurotransmitter levels leading to side and adverse therapeutic effects, the drug discovery should extend to appropriate manipulation of the 'downstream' receptors affected by the neurotransmitter boost. Efficient use of new computational strategies will accelerate the drug discovery process and reduce its economic burden.

Ligand induced conformational changes of the human serotonin transporter revealed by molecular dynamics simulations

The competitive inhibitor cocaine and the non-competitive inhibitor ibogaine induce different conformational states of the human serotonin transporter. It has been shown from accessibility experiments that cocaine mainly induces an outward-facing conformation, while the non-competitive inhibitor ibogaine, and its active metabolite noribogaine, have been proposed to induce an inward-facing conformation of the human serotonin transporter similar to what has been observed for the endogenous substrate, serotonin. The ligand induced conformational changes within the human serotonin transporter caused by these three different types of ligands, substrate, non-competitive and competitive inhibitors, are studied from multiple atomistic molecular dynamics simulations initiated from a homology model of the human serotonin transporter. The results reveal that diverse conformations of the human serotonin transporter are captured from the molecular dynamics simulations depending on the type of the ligand bound. The inward-facing conformation of the human serotonin transporter is reached with noribogaine bound, and this state resembles a previously identified inward-facing conformation of the human serotonin transporter obtained from molecular dynamics simulation with bound substrate, but also a recently published inward-facing conformation of a bacterial homolog, the leucine transporter from Aquifex Aoelicus. The differences observed in ligand induced behavior are found to originate from different interaction patterns between the ligands and the protein. Such atomic-level understanding of how an inhibitor can dictate the conformational response of a transporter by ligand binding may be of great importance for future drug design.

Amino acid secondary transporters: toward a common transport mechanism

Solute carriers (SLC) that transport amino acids are key players in health and diseases in humans. Their prokaryotic relatives are often involved in essential physiological processes in microorganisms, e.g. in homeostasis and acidic/osmotic stress response. High-resolution X-ray structures of the sequence-unrelated amino acid transporters unraveled a striking structural similarity between carriers, which were formerly assigned to different families. The highly conserved fold is characterized by two inverted structural repeats of five transmembrane helices each and indicates common mechanistic transport concepts if not an evolutionary link among a large number of amino acid transporters. Therefore, these transporters are classified now into the structural amino acid-polyamine-organocation superfamily (APCS). The APCS includes among others the mammalian SLC6 transporters and the heterodimeric SLC7/SLC3 transporters. However, it has to be noted that the APCS is not limited entirely to amino acid transporters but contains also transporters for, e.g. amino acid derivatives and sugars. For instance, the betaine-choline-carnitine transporter family of bacterial activity-regulated Na(+)- and H(+)-coupled symporters for glycine betaine and choline is also part of this second largest structural superfamily. The APCS fold provides different possibilities to transport the same amino acid. Arginine can be transported by an H(+)-coupled symport or by antiport mechanism in exchange against agmatine for example. The convergence of the mechanistic concept of transport under comparable physiological conditions allows speculating if structurally unexplored amino acid transporters, e.g. the members of the SLC36 and SLC38 family, belong to the APCS, too. In the kidney, which is an organ that depends critically on the regulated amino acid transport, these different SLC transporters have to work together to account for proper function. Here, we will summarize the basic concepts of Na(+)- and H(+)-coupled amino acid symport and amino acid-product antiport in the light of the respective physiological requirements.

Comparative modeling of the human monoamine transporters: similarities in substrate binding

The amino acid compositions of the substrate binding pockets of the three human monoamine transporters are compared as is the orientation of the endogenous substrates, serotonin, dopamine, and norepinephrine, bound in these. Through a combination of homology modeling, induced fit dockings, molecular dynamics simulations, and uptake experiments in mutant transporters, we propose a common binding mode for the three substrates. The longitudinal axis of the substrates is similarly oriented with these, forming an ionic interaction between the ammonium group and a highly conserved aspartate, Asp98 (serotonin transporter, hSERT), Asp79 (dopamine transporter, hDAT), and Asp75 (norepinephrine transporter, hNET). The 6-position of serotonin and the para-hydroxyl groups of dopamine and norepinephrine were found to face Ala173 in hSERT, Gly153 in hDAT, and Gly149 in hNET. Three rotations of the substrates around the longitudinal axis were identified. In each mode, an aromatic hydroxyl group of the substrates occupied equivalent volumes of the three binding pockets, where small changes in amino acid composition explains the differences in selectivity. Uptake experiments support that the 5-hydroxyl group of serotonin and the meta-hydroxyl group norepinephrine and dopamine are placed in the hydrophilic pocket around Ala173, Ser438, and Thr439 in hSERT corresponding to Gly149, Ser419, Ser420 in hNET and Gly153 Ser422 and Ala423 in hDAT. Furthermore, hDAT was found to possess an additional hydrophilic pocket around Ser149 to accommodate the para-hydroxyl group. Understanding these subtle differences between the binding site compositions of the three transporters is imperative for understanding the substrate selectivity, which could eventually aid in developing future selective medicines.

Binding of serotonin to lipid membranes

Serotonin (5-hydroxytryptamine, 5-HT) is a prevalent neurotransmitter throughout the animal kingdom. It exerts its effect through the specific binding to the serotonin receptor, but recent research has suggested that neural transmission may also be affected by its nonspecific interactions with the lipid matrix of the synaptic membrane. However, membrane-5-HT interactions remain controversial and superficially investigated. Fundamental knowledge of this interaction appears vital in discussions of putative roles of 5-HT, and we have addressed this by thermodynamic measurements and molecular dynamics (MD) simulations. 5-HT was found to interact strongly with lipid bilayers (partitioning coefficient ~1200 in mole fraction units), and this is highly unusual for a hydrophilic solute like 5-HT which has a bulk, oil-water partitioning coefficient well below unity. It follows that membrane affinity must rely on specific interactions, and the MD simulations identified the salt-bridge between the primary amine of 5-HT and the lipid phosphate group as the most important interaction. This interaction anchored cationic 5-HT in the membrane interface with the aromatic ring system pointing inward and a prevailing residence between the phosphate and the carbonyl groups of the lipid. The unprotonated form of 5-HT shows the opposite orientation, with the primary amine pointing toward the membrane core. Partitioning of 5-HT was found to decrease lipid chain order. These distinctive interactions of 5-HT and model membranes could be related to nonspecific effects of this neurotransmitter.

Interaction of antidepressants with the serotonin and norepinephrine transporters: mutational studies of the S1 substrate binding pocket

The serotonin transporter (SERT) and the norepinephrine transporter (NET) are sodium-dependent neurotransmitter transporters responsible for reuptake of released serotonin and norepinephrine, respectively, into nerve terminals in the brain. A wide range of inhibitors of SERT and NET are used as treatment of depression and anxiety disorders or as psychostimulant drugs of abuse. Despite their clinical importance, the molecular mechanisms by which various types of antidepressant drugs bind and inhibit SERT and NET are still elusive for the majority of the inhibitors, including the molecular basis for SERT/NET selectivity. Mutational analyses have suggested that a central substrate binding site (denoted the S1 pocket) also harbors an inhibitor binding site. In this study, we determine the effect of mutating six key S1 residues in human SERT (hSERT) and NET (hNET) on the potency of 15 prototypical SERT/NET inhibitors belonging to different drug classes. Analysis of the resulting drug sensitivity profiles provides novel information on drug binding modes in hSERT and hNET and identifies specific S1 residues as important molecular determinants for inhibitor potency and hSERT/hNET selectivity.

Monoamine transporter structure, function, dynamics, and drug discovery: a computational perspective

With the breakthrough crystallization of the bacterial leucine transporter protein LeuT, the first available X-ray structure for the neurotransmitter/sodium symporter family, development of 3-D computational models is suddenly essential for structure-function studies on the plasmalemmal monoamine transporters (MATs). LeuT-based MAT models have been used to guide elucidation of substrate and inhibitor binding pockets, and molecular dynamics simulations using these models are providing insight into conformations involved in the substrate translocation cycle. With credible MAT models finally in hand, structure-based virtual screening for novel ligands is yielding lead compounds toward the development of new medications for psychostimulant dependence, attention deficit hyperactivity, depression, anxiety, schizophrenia, and other disorders associated with dopamine, norepinephrine, or serotonin dysregulation.

SLC6 neurotransmitter transporters: structure, function, and regulation

The neurotransmitter transporters (NTTs) belonging to the solute carrier 6 (SLC6) gene family (also referred to as the neurotransmitter-sodium-symporter family or Na(+)/Cl(-)-dependent transporters) comprise a group of nine sodium- and chloride-dependent plasma membrane transporters for the monoamine neurotransmitters serotonin (5-hydroxytryptamine), dopamine, and norepinephrine, and the amino acid neurotransmitters GABA and glycine. The SLC6 NTTs are widely expressed in the mammalian brain and play an essential role in regulating neurotransmitter signaling and homeostasis by mediating uptake of released neurotransmitters from the extracellular space into neurons and glial cells. The transporters are targets for a wide range of therapeutic drugs used in treatment of psychiatric diseases, including major depression, anxiety disorders, attention deficit hyperactivity disorder and epilepsy. Furthermore, psychostimulants such as cocaine and amphetamines have the SLC6 NTTs as primary targets. Beginning with the determination of a high-resolution structure of a prokaryotic homolog of the mammalian SLC6 transporters in 2005, the understanding of the molecular structure, function, and pharmacology of these proteins has advanced rapidly. Furthermore, intensive efforts have been directed toward understanding the molecular and cellular mechanisms involved in regulation of the activity of this important class of transporters, leading to new methodological developments and important insights. This review provides an update of these advances and their implications for the current understanding of the SLC6 NTTs.

Synthesis and serotonin transporter activity of 1,3-bis(aryl)-2-nitro-1-propenes as a new class of anticancer agents

Structural derivatives of 4-MTA, an illegal amphetamine analogue have been previously shown to have anticancer effects in vitro. In this study we report the synthesis of a series of novel 1,3-bis(aryl)-2-nitro-1-propene derivatives related in structure to 4-MTA. A number of these compounds containing a classic nitrostyrene structure are shown to have antiproliferative activities in vitro in a range of malignant cell lines, particularly against Burkitt's lymphoma derived cell lines, whilst having no effect on 'normal' peripheral blood mononuclear cells. Such effects appear to be independent of the serotonin transporter, a high affinity target for amphetamines and independent of protein tyrosine phosphatases and tubulin dynamics both of which have been previously associated with nitrostyrene-induced cell death. We demonstrate that a number of these compounds induce caspase activation, PARP cleavage, chromatin condensation and membrane blebbing in a Burkitt's lymphoma derived cell line, consistent with these compounds inducing apoptosis in vitro. Although no specific target has yet been identified for the action of these compounds, the cell death elicited is potent, selective and worthy of further investigation.

Structure-activity relationships for a novel series of citalopram (1-(3-(dimethylamino)propyl)-1-(4-fluorophenyl)-1,3-dihydroisobenzofuran-5-carbonitrile) analogues at monoamine transporters

(+/-)-Citalopram (1, 1-(3-(dimethylamino)propyl)-1-(4-fluorophenyl)-1,3-dihydroisobenzofuran-5-carbonitrile), and its eutomer, escitalopram (S-(+)-1) are selective serotonin reuptake inhibitors (SSRIs) that are used clinically to treat anxiety and depression. To further explore structure-activity relationships at the serotonin transporter (SERT), a series of (+/-)-4- and 5-substituted citalopram analogues were designed, synthesized, and evaluated for binding at the SERT, dopamine transporter (DAT) and norepinephrine transporter (NET) in native rodent tissue. Many of these analogues showed high SERT binding affinities (Ki=1-40 nM) and selectivities over both NET and DAT. Selected enantiomeric pairs of analogues were synthesized and both retained enantioselectivity as with S- and R-1, wherein S>R at the SERT. In addition, the enantiomeric pairs of 1 and 5 were tested for binding at the homologous bacterial leucine transporter (LeuT), wherein low affinities and the absence of enantioselectivity suggested distinctive binding sites for these compounds at SERT as compared to LeuT. These novel ligands will provide molecular tools to elucidate drug-protein interactions at the SERT and to relate those to behavioral actions in vivo.

Transmembrane domain 6 of the human serotonin transporter contributes to an aqueously accessible binding pocket for serotonin and the psychostimulant 3,4-methylene dioxymethamphetamine

The plasma membrane serotonin (5-HT) transporter (SERT, SLC6A4) clears 5-HT after release at nerve termini and is targeted by both antidepressant medications and psychostimulants (e.g. MDMA, cocaine). Homology modeling of human SERT (hSERT), based on high resolution structures of the microbial SLC6 family member LeuT(Aa), along with biochemical studies of wild type and mutant transporters, predicts transmembrane (TM) domains 1, 3, 6, and 8 comprise the 5-HT-binding pocket. We utilized the substituted cysteine accessibility method along with surface and site-specific biotinylation to probe TM6 for aqueous accessibility and differential interactions with 5-HT and psychostimulants. Our results are consistent with TM6 being composed of an aqueous-accessible, alpha-helical extracellular domain (TM6a) that is separated by a central, unwound section from a cytoplasmically localized domain (TM6b) with limited aqueous accessibility. The substitution G338C appears to lock hSERT in an outward-facing conformation that, although accessible to aminoethylmethanethiosulfonate-biotin, 5-HT, and citalopram, is incapable of inward 5-HT transport. Transport of 5-HT by G338C can be partially restored by the TM1 mutation Y95F. With regard to methanethiosulfonate (MTS) inactivation of uptake, TM6a Cys mutants demonstrate Na(+)-dependent [2-(trimethylammonium)ethyl]-MTS sensitivity. Studies with the centrally located substitution S336C reveal features of a common binding pocket for 5-HT and 3,4-methylenedioxymethamphetamine (MDMA). Interestingly, the substitution I333C reveals an MDMA-induced conformation not observed with 5-HT. In the context of prior studies on TM1, our findings document shared and unique features of TM6 contributing to hSERT aqueous accessibility, ligand recognition, and conformational dynamics.

Mutational mapping and modeling of the binding site for (S)-citalopram in the human serotonin transporter

The serotonin transporter (SERT) regulates extracellular levels of the neurotransmitter serotonin (5-hydroxytryptamine) in the brain by facilitating uptake of released 5-hydroxytryptamine into neuronal cells. SERT is the target for widely used antidepressant drugs, including imipramine, fluoxetine, and (S)-citalopram, which are competitive inhibitors of the transport function. Knowledge of the molecular details of the antidepressant binding sites in SERT has been limited due to lack of structural data on SERT. Here, we present a characterization of the (S)-citalopram binding pocket in human SERT (hSERT) using mutational and computational approaches. Comparative modeling and ligand docking reveal that (S)-citalopram fits into the hSERT substrate binding pocket, where (S)-citalopram can adopt a number of different binding orientations. We find, however, that only one of these binding modes is functionally relevant from studying the effects of 64 point mutations around the putative substrate binding site. The mutational mapping also identify novel hSERT residues that are crucial for (S)-citalopram binding. The model defines the molecular determinants for (S)-citalopram binding to hSERT and demonstrates that the antidepressant binding site overlaps with the substrate binding site.

Recent advances in the understanding of the interaction of antidepressant drugs with serotonin and norepinephrine transporters

The biogenic monoamine transporters are integral membrane proteins that perform active transport of extracellular dopamine, serotonin and norepinephrine into cells. These transporters are targets for therapeutic agents such as antidepressants, as well as addictive substances such as cocaine and amphetamine. Seminal advances in the understanding of the structure and function of this transporter family have recently been accomplished by structural studies of a bacterial transporter, as well as medicinal chemistry and pharmacological studies of mammalian transporters. This feature article focuses on antidepressant drugs that act on the serotonin and/or the norepinephrine transporters. Specifically, we focus on structure-activity relationships of these drugs with emphasis on relationships between their molecular properties and the current knowledge of transporter structure.

Computational molecular biology approaches to ligand-target interactions

Binding of small molecules to their targets triggers complex pathways. Computational approaches are keys for predictions of the molecular events involved in such cascades. Here we review current efforts at characterizing the molecular determinants in the largest membrane-bound receptor family, the G-protein-coupled receptors (GPCRs). We focus on odorant receptors, which constitute more than half GPCRs. The work presented in this review uncovers structural and energetic aspects of components of the cellular cascade. Finally, a computational approach in the context of radioactive boron-based antitumoral therapies is briefly described.

Structure and localisation of drug binding sites on neurotransmitter transporters

The dopamine (DAT), serotontin (SERT) and noradrenalin (NET) transporters are molecular targets for different classes of psychotropic drugs. The crystal structure of Aquifex aeolicus LeuT(Aa) was used as a template for molecular modeling of DAT, SERT and NET, and two putative drug binding sites (pocket 1 and 2) in each transporter were identified. Cocaine was docked into binding pocket 1 of DAT, corresponding to the leucine binding site in LeuT(Aa), which involved transmembrane helices (TMHs) 1, 3, 6 and 8. Clomipramine was docked into binding pocket 2 of DAT, involving TMHs 1, 3, 6, 10 and 11, and extracellular loops 4 and 6, corresponding to the clomipramine binding site in a crystal structure of a LeuT(Aa)-clomipramine complex. The structures of the proposed cocaine- and tricyclic antidepressant-binding sites may be of particular interest for the design of novel DAT interacting ligands.

Antidepressants targeting the serotonin reuptake transporter act via a competitive mechanism

Although several antidepressants (including fluoxetine, imipramine, citalopram, venlafaxine, and duloxetine) are known to inhibit the serotonin transporter (SERT), whether or not these molecules compete with 5-hydroxytryptamine (serotonin) (5-HT) for binding to SERT has remained controversial. We have performed radioligand competition binding experiments and found that all data can be fitted via a simple competitive interaction model, using Cheng-Prusoff analysis (Biochem Pharmacol 22:3099-3108, 1973). Two different SERT-selective radioligands, [(3)H]N,N-dimethyl-2-(2-amino-4-cyanophenyl thio)-benzylamine (DASB) and [(3)H]S-citalopram, were used to probe competitive binding to recombinantly expressed human SERT or native SERT in rat cortical membranes. All the SERT inhibitors that we tested were able to inhibit [(3)H]DASB and [(3)H]S-citalopram binding in a concentration-dependent manner, with unity Hill coefficient. In accordance with the Cheng-Prusoff relationship for a competitive interaction, we observed that test compound concentrations associated with 50% maximal inhibition of radiotracer binding (IC(50)) increased linearly with increasing radioligand concentration for all ligands: 5-HT, S-citalopram, R-citalopram, paroxetine, clomipramine, fluvoxamine, imipramine venlafaxine, duloxetine, indatraline, cocaine, and 2-beta-carboxy-3-beta-(4-iodophenyl)tropane. The equilibrium dissociation constant of 5-HT and SERT inhibitors were also derived using Scatchard analysis of the data set, and they were found to be comparable with the data obtained using the Cheng-Prusoff relationship. Our studies establish a reference framework that will contribute to ongoing efforts to understand ligand binding modes at SERT by demonstrating that 5-HT and the SERT inhibitors tested bind to the serotonin transporter in a competitive manner.

Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system

We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.

The molecular interactions of buspirone analogues with the serotonin transporter

A major problem with the selective serotonin reuptake inhibitors (SSRIs) is the delayed onset of action. A reason for that may be that the initial SSRI-induced increase in serotonin levels activates somatodendritic 5-HT(1A) autoreceptors, causing a decrease in serotonin release in major forebrain areas. It has been suggested that compounds combining inhibition of the serotonin transport protein with antagonistic effects on the 5-HT(1A) receptor will shorten the onset time. The anxiolytic drug buspirone is known as 5-HT(1A) partial agonist. In the present work, we are studying the inhibition of the serotonin transporter protein by a series of buspirone analogues by molecular modelling and by experimental affinity measurements. Models of the transporter protein were constructed using the crystal structure of the Escherichia coli major facilitator family transporter-LacY and the X-ray structure of the neurotransmitter symporter family (NSS) transporter-LeuT(Aa) as templates. The buspirone analogues were docked into both SERT models and the interactions with amino acids within the protein were analyzed. Two putative binding sites were identified on the LeuT(Aa) based model, one suggested to be a high-affinity site, and the other suggested to be a low-affinity binding site. Molecular dynamic simulations of the LacY based model in complex with ligands did not induce a helical architecture of the LacY based model into an arrangement more similar to that of the LeuT(Aa) based model.

Substrate binding and formation of an occluded state in the leucine transporter

Translocation through the extracellular vestibule and binding of leucine in the leucine transporter (LeuT) have been studied with molecular dynamics simulations. More than 0.1 mus of all-atom molecular dynamics simulations have been performed on different combinations of LeuT, bound substrate, and bound structural Na(+) ions to describe molecular events involved in substrate binding and in the formation of the occluded state and to investigate the dynamics of this state. Three structural features are found to be directly involved in the initial steps of leucine transport: a Na(+) ion directly coordinated to leucine (Na-1), two aromatic residues closing the binding site toward the extracellular vestibule (Tyr-108 and Phe-253), and a salt bridge in the extracellular vestibule (Arg-30 and Asp-404). These features account for observed differences between simulations of LeuT with and without bound substrate and for a possible pathway for leucine binding and thereby formation of the occluded LeuT binding site.