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

Effect of drug amlodipine on the charged lipid bilayer cell membranes DMPS and DMPS + DMPC: a molecular dynamics simulation study

In this work, the effects of the anti-hypertensive drug amlodipine in native and PEGylated forms on the malfunctioning of negatively charged lipid bilayer cell membranes constructed from DMPS or DMPS + DMPC were studied by molecular dynamics simulation. The obtained results indicate that amlodipine alone aggregates and as a result its diffusion into the membrane is retarded. In addition, due to their large size aggregates of the drug can damage the cell, rupturing the cell membrane. It is shown that PEGylation of amlodipine prevents this aggregation and facilitates its diffusion into the lipid membrane. The interaction of the drug with negatively charged membranes in the presence of an aqueous solution of NaCl, as the medium, is investigated and its effects on the membrane are considered by evaluating the structural properties of the membrane such as area per lipid, thickness, lipid chain order and electrostatic potential difference between bulk solution and lipid bilayer surface. The effect of these parameters on the diffusion of the drug into the cell is critically examined and discussed.

A direct interaction of cholesterol with the dopamine transporter prevents its out-to-inward transition

Monoamine transporters (MATs) carry out neurotransmitter reuptake from the synaptic cleft, a key step in neurotransmission, which is targeted in the treatment of neurological disorders. Cholesterol (CHOL), a major component of the synaptic plasma membrane, has been shown to exhibit a modulatory effect on MATs. Recent crystal structures of the dopamine transporter (DAT) revealed the presence of two conserved CHOL-like molecules, suggesting a functional protein-CHOL direct interaction. Here, we present extensive atomistic molecular dynamics (MD) simulations of DAT in an outward-facing conformation. In the absence of bound CHOL, DAT undergoes structural changes reflecting early events of dopamine transport: transition to an inward-facing conformation. In contrast, in the presence of bound CHOL, these conformational changes are inhibited, seemingly by an immobilization of the intracellular interface of transmembrane helix 1a and 5 by CHOL. We also provide evidence, from coarse grain MD simulations that the CHOL sites observed in the DAT crystal structures are preserved in all human monoamine transporters (dopamine, serotonin and norepinephrine), suggesting that our findings might extend to the entire family.

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.

The Environment Shapes the Inner Vestibule of LeuT

Human neurotransmitter transporters are found in the nervous system terminating synaptic signals by rapid removal of neurotransmitter molecules from the synaptic cleft. The homologous transporter LeuT, found in Aquifex aeolicus, was crystallized in different conformations. Here, we investigated the inward-open state of LeuT. We compared LeuT in membranes and micelles using molecular dynamics simulations and lanthanide-based resonance energy transfer (LRET). Simulations of micelle-solubilized LeuT revealed a stable and widely open inward-facing conformation. However, this conformation was unstable in a membrane environment. The helix dipole and the charged amino acid of the first transmembrane helix (TM1A) partitioned out of the hydrophobic membrane core. Free energy calculations showed that movement of TM1A by 0.30 nm was driven by a free energy difference of ~15 kJ/mol. Distance measurements by LRET showed TM1A movements, consistent with the simulations, confirming a substantially different inward-open conformation in lipid bilayer from that inferred from the crystal structure.

Crystal structures of the bacterial small amino acid transporter LeuT provided structural evidence for the alternating access model. Thereby, these structures shaped our understanding of the mechanisms underlying substrate translocation by neurotransmitter transporters. However, it has been questioned, if the crystallized inward-open conformation of LeuT can exist in the membrane environment. Here we show that, while stable in detergent micelles, the inward-open conformation of LeuT is of high energy and undergoes structural readjustments. We use a multi-faceted approach including molecular dynamics simulations, scintillation proximity assays, free energy calculations and apply for the first time lanthanide resonance energy transfer measurements to verify the in silico predictions. In silico and in vitro approaches using the same conditions allowed us to combine the macroscopic experimental data with microscopic all atom results from simulations to identify the underlying driving forces: partitioning of charged and polar groups from the hydrophobic membrane interior to the hydrophilic environment. We propose that the inward-facing state shows a much smaller movement of TM1A, but large enough to create an access path to the S1 substrate binding site from the vestibule.

Exploring the Inhibitory Mechanism of Approved Selective Norepinephrine Reuptake Inhibitors and Reboxetine Enantiomers by Molecular Dynamics Study

Selective norepinephrine reuptake inhibitors (sNRIs) provide an effective class of approved antipsychotics, whose inhibitory mechanism could facilitate the discovery of privileged scaffolds with enhanced drug efficacy. However, the crystal structure of human norepinephrine transporter (hNET) has not been determined yet and the inhibitory mechanism of sNRIs remains elusive. In this work, multiple computational methods were integrated to explore the inhibitory mechanism of approved sNRIs (atomoxetine, maprotiline, reboxetine and viloxazine), and 3 lines of evidences were provided to verify the calculation results. Consequently, a binding mode defined by interactions between three chemical moieties in sNRIs and eleven residues in hNET was identified as shared by approved sNRIs. In the meantime, binding modes of reboxetine's enantiomers with hNET were compared. 6 key residues favoring the binding of (S, S)-reboxetine over that of (R, R)-reboxetine were discovered. This is the first study reporting that those 11 residues are the common determinants for the binding of approved sNRIs. The identified binding mode shed light on the inhibitory mechanism of approved sNRIs, which could help identify novel scaffolds with improved drug efficacy.

Identification of the inhibitory mechanism of FDA approved selective serotonin reuptake inhibitors: an insight from molecular dynamics simulation study

The binding mode shared by 4 FDA approved SSRIs treating major depression was identified by integrating multiple computational methods.

Monoamine transporters: insights from molecular dynamics simulations

The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson's disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na(+) symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors.

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 dynamics of the monoamine transporter homolog LeuT from accelerated conformational sampling and channel analysis

The bacterial leucine transporter LeuT retains significant secondary structure similarities to the human monoamine transporters (MAT) such as the dopamine and serotonin reuptake proteins. The primary method of computational study of the MATs has been through the use of the crystallized LeuT structure. Different conformations of LeuT can give insight into mechanistic details of the MAT family. A conformational sampling performed through accelerated molecular dynamics simulations testing different combinations of the leucine substrate and bound sodium ions revealed seven distinct conformational clusters. Further analysis has been performed to target salt-bridge residues R30-D404, Y108-F253, and R5-D369 and transmembrane domains on both the seven isolated structures and the total trajectories. In addition, solvent accessibility of LeuT and its substrate binding pockets has been analyzed using a program for calculating channel radii. Occupation of the Na2 site stabilizes the outward conformation and should bind to the open outward conformation before the leucine and Na1 sodium while two possible pathways were found to be available for intracellular transport.

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.

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.

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.

Insights from molecular dynamics: the binding site of cocaine in the dopamine transporter and permeation pathways of substrates in the leucine and dopamine transporters

The dopamine transporter (DAT) facilitates the regulation of synaptic neurotransmitter levels. As a target for therapeutic and illicit psycho-stimulant drugs like antidepressants and cocaine, DAT has been studied intensively. Despite a wealth of mutational and physiological data regarding DAT, the structure remains unsolved and details of the transport mechanism, binding sites and conformational changes remain debated. A bacterial homolog of DAT, the leucine transporter (LeuT(Aa)) has been used as a template and framework for modeling and understanding DAT. Free energy profiles obtained from Multi-Configuration Thermodynamic Integration simulations allowed us to correctly identify the primary and secondary binding pockets of LeuT(Aa). A comparison of free energy profiles for dopamine and cocaine in DAT suggests that the binding site of cocaine is located in a secondary pocket, not the primary substrate site. Two recurring primary pathways for intracellular substrate release from the primary pocket are identified in both transporters using the Random Acceleration Molecular Dynamics method. One pathway appears to follow transmembranes (TMs) 1a and 6b while the other pathway follows along TMs 6b and 8. Interestingly, we observe that a single sodium ion is co-transported with leucine during both simulation types.

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.

Substrate binding and translocation of the serotonin transporter studied by docking and molecular dynamics simulations

The serotonin (5-HT) transporter (SERT) plays an important role in the termination of 5-HT-mediated neurotransmission by transporting 5-HT away from the synaptic cleft and into the presynaptic neuron. In addition, SERT is the main target for antidepressant drugs, including the selective serotonin reuptake inhibitors (SSRIs). The three-dimensional (3D) structure of SERT has not yet been determined, and little is known about the molecular mechanisms of substrate binding and transport, though such information is very important for the development of new antidepressant drugs. In this study, a homology model of SERT was constructed based on the 3D structure of a prokaryotic homologous leucine transporter (LeuT) (PDB id: 2A65). Eleven tryptamine derivates (including 5-HT) and the SSRI (S)-citalopram were docked into the putative substrate binding site, and two possible binding modes of the ligands were found. To study the conformational effect that ligand binding may have on SERT, two SERT–5-HT and two SERT–(S)-citalopram complexes, as well as the SERT apo structure, were embedded in POPC lipid bilayers and comparative molecular dynamics (MD) simulations were performed. Our results show that 5-HT in the SERT–5-HT^B complex induced larger conformational changes in the cytoplasmic parts of the transmembrane helices of SERT than any of the other ligands. Based on these results, we suggest that the formation and breakage of ionic interactions with amino acids in transmembrane helices 6 and 8 and intracellular loop 1 may be of importance for substrate translocation.

Graphical abstract

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.

Single-molecule dynamics of gating in a neurotransmitter transporter homologue

Neurotransmitter:Na⁺ symporters (NSS) remove neurotransmitters from the synapse in a reuptake process driven by the Na⁺ gradient. Drugs that interfere with this reuptake mechanism, such as cocaine and antidepressants, profoundly influence behavior and mood. In order to probe the nature of conformational changes associated with substrate binding and transport, we have developed a single-molecule fluorescence imaging assay, in combination with functional and computational studies, using the prokaryotic NSS homolog LeuT. Here we show molecular details of the modulation of intracellular gating of LeuT by substrates and inhibitors, as well as by mutations that alter binding and/or transport. Our direct observations of single-molecule transitions, reflecting structural dynamics of the intracellular region of the transporter that may be masked by ensemble averaging or suppressed under crystallographic conditions, are interpreted in the context of an allosteric mechanism coupling ion and substrate binding to transport.

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.

Molecular dynamics of leucine and dopamine transporter proteins in a model cell membrane lipid bilayer

The dopamine transporter (DAT) operates via facilitated diffusion, harnessing an inward Na(+) gradient to drive dopamine from the extracellular synaptic cleft to the neuron interior. The DAT is relevant to central nervous system disorders such as Parkinson disease and attention-deficit hyperactivity disorder and is the primary site of action for the abused psychostimulants cocaine and amphetamines. Crystallization of a DAT homolog, the bacterial leucine transporter LeuT, provided the first reliable 3-D DAT template. Here, the LeuT crystal structure and the DAT molecular model have been combined with their respective substrates, leucine and dopamine, in lipid bilayer molecular dynamics simulations toward tracking substrate movement along the protein's substrate/ion permeation pathway. Specifically, movement of residue pairs that comprise the "external gate" was followed as a function of substrate presence. The transmembrane (TM) 1 arginine-TM 10 aspartate strut formed less readily in DAT compared with LeuT, with or without substrate present. For LeuT but not DAT, the addition of substrate enhanced the chances of forming the TM 1-10 bridge. Also, movement of the fourth extracellular loop EL-4 in the presence of substrate was more pronounced for DAT, the EL-4 unwinding to a degree. The overall similarity between the LeuT and DAT molecular dynamics simulations indicated that LeuT was a legitimate model to guide DAT structure-function predictions. There were, nevertheless, differences significant enough to allow for DAT-unique insights, which may include how cocaine, methylphenidate (Ritalin, NIDA Drug Supply, Rockville, MD), and other DAT blockers are not recognized as substrates even though they can access the primary substrate binding pocket. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

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.

Structural determinants of species-selective substrate recognition in human and Drosophila serotonin transporters revealed through computational docking studies

To identify potential determinants of substrate selectivity in serotonin (5-HT) transporters (SERT), models of human and Drosophila serotonin transporters (hSERT, dSERT) were built based on the leucine transporter (LeuT(Aa)) structure reported by Yamashita et al. (Nature 2005;437:215-223), PBDID 2A65. Although the overall amino acid identity between SERTs and the LeuT(Aa) is only 17%, it increases to above 50% in the first shell of the putative 5-HT binding site, allowing de novo computational docking of tryptamine derivatives in atomic detail. Comparison of hSERT and dSERT complexed with substrates pinpoints likely structural determinants for substrate binding. Forgoing the use of experimental transport and binding data of tryptamine derivatives for construction of these models enables us to critically assess and validate their predictive power: A single 5-HT binding mode was identified that retains the amine placement observed in the LeuT(Aa) structure, matches site-directed mutagenesis and substituted cysteine accessibility method (SCAM) data, complies with support vector machine derived relations activity relations, and predicts computational binding energies for 5-HT analogs with a significant correlation coefficient (R = 0.72). This binding mode places 5-HT deep in the binding pocket of the SERT with the 5-position near residue hSERT A169/dSERT D164 in transmembrane helix 3, the indole nitrogen next to residue Y176/Y171, and the ethylamine tail under residues F335/F327 and S336/S328 within 4 A of residue D98. Our studies identify a number of potential contacts whose contribution to substrate binding and transport was previously unsuspected.

Driving forces for ligand migration in the leucine transporter

The potential mechanisms for driving a ligand to and through a transporter are examined for the leucine transporter using the recently published X-ray structure. Through analyses with computational methods, including investigation of electrostatic properties, site and channel identification, and docking studies, a picture emerges whereby dipolar patches which characterize the electrostatic field serve as the primary driving force to attract the ligand to the protein and begin propagation into it. The electrostatic forces are then augmented by hydrophobic forces in the transport stages, with conformational changes in the protein helping to accommodate the migration. We identify 12 sites that might be involved in ligand recognition and migration. One of these sites corresponds to the tricyclic antidepressant binding site observed in the recently published X-ray structures.

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.

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.

Molecular basis of phospholipase A2 activity toward phospholipids with sn-1 substitutions

We studied secretory phospholipase A(2) type IIA (sPLA(2)) activity toward phospholipids that are derivatized in the sn-1 position of the glycerol backbone. We explored what type of side group (small versus bulky groups, hydrophobic versus polar groups) can be introduced at the sn-1 position of the glycerol backbone of glycerophospholipids and at the same time be hydrolyzed by sPLA(2). The biophysical characterization revealed that the modified phospholipids can form multilamellar vesicles, and several of the synthesized sn-1 functionalized phospholipids were hydrolyzed by sPLA(2). Molecular dynamics simulations provided detailed insight on an atomic level that can explain the observed sPLA(2) activity toward the different phospholipid analogs. The simulations revealed that, depending on the nature of the side chain located at the sn-1 position, the group may interfere with an incoming water molecule that acts as the nucleophile in the enzymatic reaction. The simulation results are in agreement with the experimentally observed sPLA(2) activity toward the different phospholipid analogs.

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.