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

Applications and Potential of In Silico Approaches for Psychedelic Chemistry

Molecular-level investigations of the Central Nervous System have been revolutionized by the development of computational methods, computing power, and capacity advances. These techniques have enabled researchers to analyze large amounts of data from various sources, including genomics, in vivo, and in vitro drug tests. In this review, we explore how computational methods and informatics have contributed to our understanding of mental health disorders and the development of novel drugs for neurological diseases, with a special focus on the emerging field of psychedelics. In addition, the use of state-of-the-art computational methods to predict the potential of drug compounds and bioinformatic tools to integrate disparate data sources to create predictive models is also discussed. Furthermore, the challenges associated with these methods, such as the need for large datasets and the diversity of in vitro data, are explored. Overall, this review highlights the immense potential of computational methods and informatics in Central Nervous System research and underscores the need for continued development and refinement of these techniques and more inclusion of Quantitative Structure-Activity Relationships (QSARs).

Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 μM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.

Single amino acid residue mediates reciprocal specificity in two mosquito odorant receptors

The southern house mosquito, Culex quinquefasciatus, utilizes two odorant receptors, CquiOR10 and CquiOR2, narrowly tuned to oviposition attractants and well conserved among mosquito species. They detect skatole and indole, respectively, with reciprocal specificity. We swapped the transmembrane (TM) domains of CquiOR10 and CquiOR2 and identified TM2 as a specificity determinant. With additional mutations, we showed that CquiOR10A73L behaved like CquiOR2. Conversely, CquiOR2L74A recapitulated CquiOR10 specificity. Next, we generated structural models of CquiOR10 and CquiOR10A73L using RoseTTAFold and AlphaFold and docked skatole and indole using RosettaLigand. These modeling studies suggested space-filling constraints around A73. Consistent with this hypothesis, CquiOR10 mutants with a bulkier residue (Ile, Val) were insensitive to skatole and indole, whereas CquiOR10A73G retained the specificity to skatole and showed a more robust response than the wildtype receptor CquiOR10. On the other hand, Leu to Gly mutation of the indole receptor CquiOR2 reverted the specificity to skatole. Lastly, CquiOR10A73L, CquiOR2, and CquiOR2L74I were insensitive to 3-ethylindole, whereas CquiOR2L74A and CquiOR2L74G gained activity. Additionally, CquiOR10A73G gave more robust responses to 3-ethylindole than CquiOR10. Thus, we suggest the specificity of these receptors is mediated by a single amino acid substitution, leading to finely tuned volumetric space to accommodate specific oviposition attractants.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate–protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Glycine Transporter 2: Mechanism and Allosteric Modulation

Neurotransmitter sodium symporters (NSS) are a subfamily of SLC6 transporters responsible for regulating neurotransmitter signalling. They are a major target for psychoactive substances including antidepressants and drugs of abuse, prompting substantial research into their modulation and structure-function dynamics. Recently, a series of allosteric transport inhibitors have been identified, which may reduce side effect profiles, compared to orthosteric inhibitors. Allosteric inhibitors are also likely to provide different clearance kinetics compared to competitive inhibitors and potentially better clinical outcomes. Crystal structures and homology models have identified several allosteric modulatory sites on NSS including the vestibule allosteric site (VAS), lipid allosteric site (LAS) and cholesterol binding site (CHOL1). Whilst the architecture of eukaryotic NSS is generally well conserved there are differences in regions that form the VAS, LAS, and CHOL1. Here, we describe ligand-protein interactions that stabilize binding in each allosteric site and explore how differences between transporters could be exploited to generate NSS specific compounds with an emphasis on GlyT2 modulation.

Recent Advances and Challenges of the Drugs Acting on Monoamine Transporters

Background:

The human Monoamine Transporters (hMATs), primarily including hSERT, hNET and hDAT, are important targets for the treatment of depression and other behavioral disorders with more than the availability of 30 approved drugs.

Objective:

This paper is to review the recent progress in the binding mode and inhibitory mechanism of hMATs inhibitors with the central or allosteric binding sites, for the benefit of future hMATs inhibitor design and discovery. The Structure-Activity Relationship (SAR) and the selectivity for hit/lead compounds to hMATs that are evaluated by in vitro and in vivo experiments will be highlighted.

Methods:

PubMed and Web of Science databases were searched for protein-ligand interaction, novel inhibitors design and synthesis studies related to hMATs.

Results:

Literature data indicate that since the first crystal structure determinations of the homologous bacterial Leucine Transporter (LeuT) complexed with clomipramine, a sizable database of over 100 experimental structures or computational models has been accumulated that now defines a substantial degree of structural variability hMATs-ligands recognition. In the meanwhile, a number of novel hMATs inhibitors have been discovered by medicinal chemistry with significant help from computational models.

Conclusion:

The reported new compounds act on hMATs as well as the structures of the transporters complexed with diverse ligands by either experiment or computational modeling have shed light on the poly-pharmacology, multimodal and allosteric regulation of the drugs to transporters. All of the studies will greatly promote the Structure-Based Drug Design (SBDD) of structurally novel scaffolds with high activity and selectivity for hMATs.

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'.

Mapping the Binding Site for Escitalopram and Paroxetine in the Human Serotonin Transporter Using Genetically Encoded Photo-Cross-Linkers

In spite of the important role of the human serotonin transporter (hSERT) in depression treatment, the molecular details of how antidepressant drugs bind are still not completely understood, in particular those related to potential high- and low-affinity binding sites in hSERT. Here, we utilize amber codon suppression in hSERT to encode the photo-cross-linking unnatural amino acid p-azido-l-phenylalanine into the suggested high- and low-affinity binding sites. We then employ UV-induced cross-linking with azF to map the binding site of escitalopram and paroxetine, two prototypical selective serotonin reuptake inhibitors (SSRIs). We find that the two antidepressant drugs exclusively cross-link to azF incorporated at the high-affinity binding site of hSERT, while cross-linking is not observed at the low-affinity binding site. Combined with previous homology models and recent structural data on hSERT, our results provide important information to understand the molecular details of these clinical relevant binding sites.

Investigation into the Mode of Phosphate Activation in the 4-Hydroxy-4-Methyl-2-Oxoglutarate/4-Carboxy-4-Hydroxy-2-Oxoadipate Aldolase from Pseudomonas putida F1

The 4-hydroxy-4-methyl-2-oxoglutarate (HMG)/4-carboxy-4-hydroxy-2-oxoadipate (CHA) aldolase is the last enzyme of both the gallate and protocatechuate 4,5-cleavage pathways which links aromatic catabolism to central cellular metabolism. The enzyme is a class II, divalent metal dependent, aldolase which is activated in the presence of inorganic phosphate (Pi), increasing its turnover rate >10-fold. This phosphate activation is unique for a class II aldolase. The aldolase pyruvate methyl proton exchange rate, a probe of the general acid half reaction, was increased 300-fold in the presence of 1 mM Pi and the rate enhancement followed saturation kinetics giving rise to a KM of 397 ± 30 μM. Docking studies revealed a potential Pi binding site close to, or overlapping with, the proposed general acid water site. Putative Pi binding residues were substituted by site-directed mutagenesis which resulted in reductions of Pi activation. Significantly, the active site residue Arg-123, known to be critical for the catalytic mechanism of the enzyme, was also implicated in supporting Pi mediated activation.

Mechanism of Paroxetine (Paxil) Inhibition of the Serotonin Transporter

The serotonin transporter (SERT) is an integral membrane protein that exploits preexisting sodium-, chloride-, and potassium ion gradients to catalyze the thermodynamically unfavorable movement of synaptic serotonin into the presynaptic neuron. SERT has garnered significant clinical attention partly because it is the target of multiple psychoactive agents, including the antidepressant paroxetine (Paxil), the most potent selective serotonin reuptake inhibitor known. However, the binding site and orientation of paroxetine in SERT remain controversial. To provide molecular insight, we constructed SERT homology models based on the Drosophila melanogaster dopamine transporter and docked paroxetine to these models. We tested the predicted binding configurations with a combination of radioligand binding and flux assays on wild-type and mutant SERTs. Our data suggest that the orientation of paroxetine, specifically its fluorophenyl ring, in SERT’s substrate binding site directly depends on this pocket’s charge distribution, and thereby provide an avenue toward understanding and enhancing high-affinity antidepressant activity.

Computational and biochemical docking of the irreversible cocaine analog RTI 82 directly demonstrates ligand positioning in the dopamine transporter central substrate-binding site

The dopamine transporter (DAT) functions as a key regulator of dopaminergic neurotransmission via re-uptake of synaptic dopamine (DA). Cocaine binding to DAT blocks this activity and elevates extracellular DA, leading to psychomotor stimulation and addiction, but the mechanisms by which cocaine interacts with DAT and inhibits transport remain incompletely understood. Here, we addressed these questions using computational and biochemical methodologies to localize the binding and adduction sites of the photoactivatable irreversible cocaine analog 3β-(p-chlorophenyl)tropane-2β-carboxylic acid, 4'-azido-3'-iodophenylethyl ester ([(125)I]RTI 82). Comparative modeling and small molecule docking indicated that the tropane pharmacophore of RTI 82 was positioned in the central DA active site with an orientation that juxtaposed the aryliodoazide group for cross-linking to rat DAT Phe-319. This prediction was verified by focused methionine substitution of residues flanking this site followed by cyanogen bromide mapping of the [(125)I]RTI 82-labeled mutants and by the substituted cysteine accessibility method protection analyses. These findings provide positive functional evidence linking tropane pharmacophore interaction with the core substrate-binding site and support a competitive mechanism for transport inhibition. This synergistic application of computational and biochemical methodologies overcomes many uncertainties inherent in other approaches and furnishes a schematic framework for elucidating the ligand-protein interactions of other classes of DA transport inhibitors.

Conformational changes in dopamine transporter intracellular regions upon cocaine binding and dopamine translocation

The dopamine transporter (DAT), a member of the neurotransmitter:sodium symporter family, mediates the reuptake of dopamine at the synaptic cleft. DAT is the primary target for psychostimulants such as cocaine and amphetamine. We previously demonstrated that cocaine binding and dopamine transport alter the accessibility of Cys342 in the third intracellular loop (IL3). To study the conformational changes associated with the functional mechanism of the transporter, we made cysteine substitution mutants, one at a time, from Phe332 to Ser351 in IL3 of the background DAT construct, X7C, in which 7 endogenous cysteines were mutated. The accessibility of the 20 engineered cysteines to polar charged sulfhydryl reagents was studied in the absence and presence of cocaine or dopamine. Of the 11 positions that reacted with methanethiosulfonate ethyl ammonium, as evidenced by inhibition of ligand binding, 5 were protected against this inhibition by cocaine and dopamine (S333C, S334C, N336C, M342C and T349C), indicating that reagent accessibility is affected by conformational changes associated with inhibitor and substrate binding. In some of the cysteine mutants, transport activity is disrupted, but can be rescued by the presence of zinc, most likely because the distribution between inward- and outward-facing conformations is restored by zinc binding. The experimental data were interpreted in the context of molecular models of DAT in both the inward- and outward-facing conformations. Differences in the solvent accessible surface area for individual IL3 residues calculated for these states correlate well with the experimental accessibility data, and suggest that protection by ligand binding results from the stabilization of the outward-facing configuration. Changes in the residue interaction networks observed from the molecular dynamics simulations also revealed the critical roles of several positions during the conformational transitions. We conclude that the IL3 region of DAT undergoes significant conformational changes in transitions necessary for both cocaine binding and substrate transport.

Molecular basis for selective serotonin reuptake inhibition by the antidepressant agent fluoxetine (Prozac)

Inhibitors of the serotonin transporter (SERT) are widely used antidepressant agents, but the structural mechanism for inhibitory activity and selectivity over the closely related norepinephrine transporter (NET) is not well understood. Here we use a combination of chemical, biological, and computational methods to decipher the molecular basis for high-affinity recognition in SERT and selectivity over NET for the prototypical antidepressant drug fluoxetine (Prozac; Eli Lilly, Indianapolis, IN). We show that fluoxetine binds within the central substrate site of human SERT, in agreement with recent X-ray crystal structures of LeuBAT, an engineered monoamine-like version of the bacterial amino acid transporter LeuT. However, the binding orientation of fluoxetine is reversed in our experimentally supported model compared with the LeuBAT structures, emphasizing the need for careful experimental verification when extrapolating findings from crystal structures of bacterial transporters to human relatives. We find that the selectivity of fluoxetine and nisoxetine, a NET selective structural congener of fluoxetine, is controlled by residues in different regions of the transporters, indicating a complex mechanism for selective recognition of structurally similar compounds in SERT and NET. Our findings add important new information on the molecular basis for SERT/NET selectivity of antidepressants, and provide the first assessment of the potential of LeuBAT as a model system for antidepressant binding in human transporters, which is essential for future structure-based drug development of antidepressant drugs with fine-tuned transporter selectivity.

Computational methods in drug discovery

Computer-aided drug discovery/design methods have played a major role in the development of therapeutically important small molecules for over three decades. These methods are broadly classified as either structure-based or ligand-based methods. Structure-based methods are in principle analogous to high-throughput screening in that both target and ligand structure information is imperative. Structure-based approaches include ligand docking, pharmacophore, and ligand design methods. The article discusses theory behind the most important methods and recent successful applications. Ligand-based methods use only ligand information for predicting activity depending on its similarity/dissimilarity to previously known active ligands. We review widely used ligand-based methods such as ligand-based pharmacophores, molecular descriptors, and quantitative structure-activity relationships. In addition, important tools such as target/ligand data bases, homology modeling, ligand fingerprint methods, etc., necessary for successful implementation of various computer-aided drug discovery/design methods in a drug discovery campaign are discussed. Finally, computational methods for toxicity prediction and optimization for favorable physiologic properties are discussed with successful examples from literature.

The two Na+ sites in the human serotonin transporter play distinct roles in the ion coupling and electrogenicity of transport

Neurotransmitter transporters of the SLC6 family of proteins, including the human serotonin transporter (hSERT), utilize Na(+), Cl(-), and K(+) gradients to induce conformational changes necessary for substrate translocation. Dysregulation of ion movement through monoamine transporters has been shown to impact neuronal firing potentials and could play a role in pathophysiologies, such as depression and anxiety. Despite multiple crystal structures of prokaryotic and eukaryotic SLC transporters indicating the location of both (or one) conserved Na(+)-binding sites (termed Na1 and Na2), much remains uncertain in regard to the movements and contributions of these cation-binding sites in the transport process. In this study, we utilize the unique properties of a mutation of hSERT at a single, highly conserved asparagine on TM1 (Asn-101) to provide several lines of evidence demonstrating mechanistically distinct roles for Na1 and Na2. Mutations at Asn-101 alter the cation dependence of the transporter, allowing Ca(2+) (but not other cations) to functionally replace Na(+) for driving transport and promoting 5-hydroxytryptamine (5-HT)-dependent conformational changes. Furthermore, in two-electrode voltage clamp studies in Xenopus oocytes, both Ca(2+) and Na(+) illicit 5-HT-induced currents in the Asn-101 mutants and reveal that, although Ca(2+) promotes substrate-induced current, it does not appear to be the charge carrier during 5-HT transport. These findings, in addition to functional evaluation of Na1 and Na2 site mutants, reveal separate roles for Na1 and Na2 and provide insight into initiation of the translocation process as well as a mechanism whereby the reported SERT stoichiometry can be obtained despite the presence of two putative Na(+)-binding sites.

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.

Stereoselective inhibition of serotonin transporters by antimalarial compounds

The serotonin (5-HT) transporter (SERT) is an integral membrane protein that functions to reuptake 5-HT released into the synapse following neurotransmission. This role serves an important regulatory mechanism in neuronal homeostasis. Previous studies have demonstrated that several clinically important antimalarial compounds inhibit serotonin (5-hydroxytryptamine, 5-HT) reuptake. In this study, we examined the details of antimalarial inhibition of 5-HT transport in both Drosophila (dSERT) and human SERT (hSERT) using electrophysiologic, biochemical and computational approaches. We found that the cinchona alkaloids quinidine and cinchonine, which have identical stereochemistry about carbons 8 and 9, exhibited the greatest inhibition of dSERT and hSERT transporter function whereas quinine and cinchonidine, enantiomers of quinidine and cinchonine, respectively, were weaker inhibitors of dSERT and hSERT. Furthermore, SERT mutations known to decrease the binding affinity of many antidepressants affected the cinchona alkaloids in a stereo-specific manner where the similar inhibitory profiles for quinine and cinchonidine (8S,9R) were distinct from quinidine and cinchonine (8R,9S). Small molecule docking studies with hSERT homology models predict that quinine and cinchonidine bind to the central 5-HT binding site (S1) whereas quinidine and cinchonine bind to the S2 site. Taken together, the data presented here support binding of cinchona alkaloids to two different sites on SERT defined by their stereochemistry which implies separate modes of transporter inhibition. Notably, the most potent antimalarial inhibitors of SERT appear to preferentially bind to the S2 site. Our findings provide important insight related to how this class of drugs can modulate the serotonergic system as well as identify compounds that may discriminate between the S1 and S2 binding sites and serve as lead compounds for novel SERT inhibitors.

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.

Computational design of protein-small molecule interfaces

The computational design of proteins that bind small molecule ligands is one of the unsolved challenges in protein engineering. It is complicated by the relatively small size of the ligand which limits the number of intermolecular interactions. Furthermore, near-perfect geometries between interacting partners are required to achieve high binding affinities. For apolar, rigid small molecules the interactions are dominated by short-range van der Waals forces. As the number of polar groups in the ligand increases, hydrogen bonds, salt bridges, cation-π, and π-π interactions gain importance. These partial covalent interactions are longer ranged, and additionally, their strength depends on the environment (e.g. solvent exposure). To assess the current state of protein-small molecule interface design, we benchmark the popular computer algorithm Rosetta on a diverse set of 43 protein-ligand complexes. On average, we achieve sequence recoveries in the binding site of 59% when the ligand is allowed limited reorientation, and 48% when the ligand is allowed full reorientation. When simulating the redesign of a protein binding site, sequence recovery among residues that contribute most to binding was 52% when slight ligand reorientation was allowed, and 27% when full ligand reorientation was allowed. As expected, sequence recovery correlates with ligand displacement.

Small-molecule ligand docking into comparative models with Rosetta

Structure-based drug design is frequently used to accelerate the development of small-molecule therapeutics. Although substantial progress has been made in X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy, the availability of high-resolution structures is limited owing to the frequent inability to crystallize or obtain sufficient NMR restraints for large or flexible proteins. Computational methods can be used to both predict unknown protein structures and model ligand interactions when experimental data are unavailable. This paper describes a comprehensive and detailed protocol using the Rosetta modeling suite to dock small-molecule ligands into comparative models. In the protocol presented here, we review the comparative modeling process, including sequence alignment, threading and loop building. Next, we cover docking a small-molecule ligand into the protein comparative model. In addition, we discuss criteria that can improve ligand docking into comparative models. Finally, and importantly, we present a strategy for assessing model quality. The entire protocol is presented on a single example selected solely for didactic purposes. The results are therefore not representative and do not replace benchmarks published elsewhere. We also provide an additional tutorial so that the user can gain hands-on experience in using Rosetta. The protocol should take 5-7 h, with additional time allocated for computer generation of models.

Synthesis, antidepressant evaluation and docking studies of long-chain alkylnitroquipazines as serotonin transporter inhibitors

Twelve alkyl analogues (1-12) of the high-affinity serotonin transporter (SERT) inhibitor 6-nitroquipazine (6-NQ) were synthesized and studied using in vitro radioligand competition binding assays to determine their binding affinity (Ki ). The putative antidepressant activity of five of the binders with the highest SERT binding affinities was studied by the forced swim and locomotor activity mouse tests. The three-dimensional (3D) structures of 8 and 9 were determined using NOE NMR technique. Flexible docking of the compounds was undertaken to illustrate the binding of the compounds in the SERT model. Our results showed that several of the 6-NQ analogues are high-affinity SERT inhibitors and indicated that the octyl (8), decyl (10) and dodecyl (12) 6-NQ analogues exhibit moderate antidepressant activity.

Mutational analysis of the high-affinity zinc binding site validates a refined human dopamine transporter homology model

The high-resolution crystal structure of the leucine transporter (LeuT) is frequently used as a template for homology models of the dopamine transporter (DAT). Although similar in structure, DAT differs considerably from LeuT in a number of ways: (i) when compared to LeuT, DAT has very long intracellular amino and carboxyl termini; (ii) LeuT and DAT share a rather low overall sequence identity (22%) and (iii) the extracellular loop 2 (EL2) of DAT is substantially longer than that of LeuT. Extracellular zinc binds to DAT and restricts the transporter‚s movement through the conformational cycle, thereby resulting in a decrease in substrate uptake. Residue H293 in EL2 praticipates in zinc binding and must be modelled correctly to allow for a full understanding of its effects. We exploited the high-affinity zinc binding site endogenously present in DAT to create a model of the complete transmemberane domain of DAT. The zinc binding site provided a DAT-specific molecular ruler for calibration of the model. Our DAT model places EL2 at the transporter lipid interface in the vicinity of the zinc binding site. Based on the model, D206 was predicted to represent a fourth co-ordinating residue, in addition to the three previously described zinc binding residues H193, H375 and E396. This prediction was confirmed by mutagenesis: substitution of D206 by lysine and cysteine affected the inhibitory potency of zinc and the maximum inhibition exerted by zinc, respectively. Conversely, the structural changes observed in the model allowed for rationalizing the zinc-dependent regulation of DAT: upon binding, zinc stabilizes the outward-facing state, because its first coordination shell can only be completed in this conformation. Thus, the model provides a validated solution to the long extracellular loop and may be useful to address other aspects of the transport cycle.

The dopamine transporter (DAT) regulates dopaminergic neurotransmission in the brain and is implicated in numerous human disease states. DAT is unique among the monoamine neurotransmitter transporter family because its substrate transport is inhibited by extracellular zinc. DAT homology models rely upon the crystal structure of LeuT solved in 2005. LeuT and DAT share a relatively low overall sequence identity of 22%. In addition, the length of the second extracellular loop of DAT exceeds that of LeuT by 21 residues. The zinc binding site cannot be directly modeled from the LeuT template alone because of these differences. Current available homology models of DAT focused on substrate or inhibitor binding rather than on the second extracellular loop. We exploited the specificity of the zinc binding site to build and calibrate a DAT homology model of the complete transmembrane domain. Our model predicted that the zinc binding site in DAT consists of four zinc co-ordinating residues rather than three that had been previously identified. We verified this hypothesis by site-directed mutagenesis and uptake inhibition studies.

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.

Using RosettaLigand for small molecule docking into comparative models

Computational small molecule docking into comparative models of proteins is widely used to query protein function and in the development of small molecule therapeutics. We benchmark RosettaLigand docking into comparative models for nine proteins built during CASP8 that contain ligands. We supplement the study with 21 additional protein/ligand complexes to cover a wider space of chemotypes. During a full docking run in 21 of the 30 cases, RosettaLigand successfully found a native-like binding mode among the top ten scoring binding modes. From the benchmark cases we find that careful template selection based on ligand occupancy provides the best chance of success while overall sequence identity between template and target do not appear to improve results. We also find that binding energy normalized by atom number is often less than -0.4 in native-like binding modes.

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.

Single molecule analysis of serotonin transporter regulation using antagonist-conjugated quantum dots reveals restricted, p38 MAPK-dependent mobilization underlying uptake activation

The presynaptic serotonin (5-HT) transporter (SERT) is targeted by widely prescribed antidepressant medications. Altered SERT expression or regulation has been implicated in multiple neuropsychiatric disorders, including anxiety, depression and autism. Here, we implement a generalizable strategy that exploits antagonist-conjugated quantum dots (Qdots) to monitor, for the first time, single SERT proteins on the surface of serotonergic cells. We document two pools of SERT proteins defined by lateral mobility, one that exhibits relatively free diffusion, and a second, localized to cholesterol and GM1 ganglioside-enriched microdomains, that displays restricted mobility. Receptor-linked signaling pathways that enhance SERT activity mobilize transporters that, nonetheless, remain confined to membrane microdomains. Mobilization of transporters arises from a p38 MAPK-dependent untethering of the SERT C terminus from the juxtamembrane actin cytoskeleton. Our studies establish the utility of ligand-conjugated Qdots for analysis of the behavior of single membrane proteins and reveal a physical basis for signaling-mediated SERT regulation.

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.

Bioinformatics and variability in drug response: a protein structural perspective

Marketed drugs frequently perform worse in clinical practice than in the clinical trials on which their approval is based. Many therapeutic compounds are ineffective for a large subpopulation of patients to whom they are prescribed; worse, a significant fraction of patients experience adverse effects more severe than anticipated. The unacceptable risk-benefit profile for many drugs mandates a paradigm shift towards personalized medicine. However, prior to adoption of patient-specific approaches, it is useful to understand the molecular details underlying variable drug response among diverse patient populations. Over the past decade, progress in structural genomics led to an explosion of available three-dimensional structures of drug target proteins while efforts in pharmacogenetics offered insights into polymorphisms correlated with differential therapeutic outcomes. Together these advances provide the opportunity to examine how altered protein structures arising from genetic differences affect protein-drug interactions and, ultimately, drug response. In this review, we first summarize structural characteristics of protein targets and common mechanisms of drug interactions. Next, we describe the impact of coding mutations on protein structures and drug response. Finally, we highlight tools for analysing protein structures and protein-drug interactions and discuss their application for understanding altered drug responses associated with protein structural variants.

Rosetta Ligand docking with flexible XML protocols

RosettaLigand is premiere software for predicting how a protein and a small molecule interact. Benchmark studies demonstrate that 70% of the top scoring RosettaLigand predicted interfaces are within 2Å RMSD from the crystal structure [1]. The latest release of Rosetta ligand software includes many new features, such as (1) docking of multiple ligands simultaneously, (2) representing ligands as fragments for greater flexibility, (3) redesign of the interface during docking, and (4) an XML script based interface that gives the user full control of the ligand docking protocol.

Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

The recent determination of high-resolution crystal structures of several transporters offers unprecedented insights into the structural mechanisms behind secondary transport. These proteins utilize the facilitated diffusion of the ions down their electrochemical gradients to transport the substrate against its concentration gradient. The structural studies revealed striking similarities in the structural organization of ion and solute binding sites and a well-conserved inverted-repeat topology between proteins from several gene families. In this paper we will overview recent atomistic simulations applied to study the mechanisms of selective binding of ion and substrate in LeuT, Glt, vSGLT and hSERT as well as its consequences for the transporter conformational dynamics. This article is part of a Special Issue entitled: Membrane protein structure and function.

Probing binding pocket of serotonin transporter by single molecular force spectroscopy on living cells

The serotonin transporter (SERT) terminates neurotransmission by removing serotonin from the synaptic cleft. In addition, it is the site of action of antidepressants (which block the transporter) and of amphetamines (which induce substrate efflux). The interaction energies involved in binding of such compounds to the transporter are unknown. Here, we used atomic force microscopy (AFM) to probe single molecular interactions between the serotonin transporter and MFZ2-12 (a potent cocaine analog) in living CHOK1 cells. For the AFM measurements, MFZ2-12 was immobilized on AFM tips by using a heterobifunctional cross-linker. By varying the pulling velocity in force distance cycles drug-transporter complexes were ruptured at different force loadings allowing for mapping of the interaction energy landscape. We derived chemical rate constants from these recordings and compared them with those inferred from inhibition of transport and ligand binding: k_off values were in good agreement with those derived from uptake experiments; in contrast, the k_on values were scaled down when determined by AFM. Our observations generated new insights into the energy landscape of the interaction between SERT and inhibitors. They thus provide a useful framework for molecular dynamics simulations by exploring the range of forces and energies that operate during the binding reaction.

Molecular mechanism of serotonin transporter inhibition elucidated by a new flexible docking protocol

The two main groups of antidepressant drugs, the tricyclic antidepressants (TCAs) and the selective serotonin reuptake inhibitors (SSRIs), as well as several other compounds, act by inhibiting the serotonin transporter (SERT). However, the binding mode and molecular mechanism of inhibition in SERT are not fully understood. In this study, five classes of SERT inhibitors were docked into an outward-facing SERT homology model using a new 4D ensemble docking protocol. Unlike other docking protocols, where protein flexibility is not considered or is highly dependent on the ligand structure, flexibility was here obtained by side chain sampling of the amino acids of the binding pocket using biased probability Monte Carlo (BPMC) prior to docking. This resulted in the generation of multiple binding pocket conformations that the ligands were docked into. The docking results showed that the inhibitors were stacked between the aromatic amino acids of the extracellular gate (Y176, F335) presumably preventing its closure. The inhibitors interacted with amino acids in both the putative substrate binding site and more extracellular regions of the protein. A general structure-docking-based pharmacophore model was generated to explain binding of all studied classes of SERT inhibitors. Docking of a test set of actives and decoys furthermore showed that the outward-facing ensemble SERT homology model consistently and selectively scored the majority of active compounds above decoys, which indicates its usefulness in virtual screening.

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.

A conserved asparagine residue in transmembrane segment 1 (TM1) of serotonin transporter dictates chloride-coupled neurotransmitter transport

Na(+)- and Cl(-)-dependent uptake of neurotransmitters via transporters of the SLC6 family, including the human serotonin transporter (SLC6A4), is critical for efficient synaptic transmission. Although residues in the human serotonin transporter involved in direct Cl(-) coordination of human serotonin transport have been identified, the role of Cl(-) in the transport mechanism remains unclear. Through a combination of mutagenesis, chemical modification, substrate and charge flux measurements, and molecular modeling studies, we reveal an unexpected role for the highly conserved transmembrane segment 1 residue Asn-101 in coupling Cl(-) binding to concentrative neurotransmitter uptake.

The natural resistance-associated macrophage protein from the protozoan parasite Perkinsus marinus mediates iron uptake

Microbial pathogens succeed in acquiring essential metals such as iron and manganese despite their limited availability because of the host's immune response. The eukaryotic natural resistance-associated macrophage proteins mediate uptake of divalent metals and, during infection, may compete directly for metal acquisition with the pathogens' transporters. In this study, we characterize the Nramp gene family of Perkinsus marinus, an intracellular parasite of the eastern oyster, and through yeast complementation, we demonstrate for the first time for a protozoan parasite that Nramp imports environmental Fe. Three PmNramp isogenes differ in their exon-intron structures and encode transcripts that display a trans splicing leader at the 5' end. The protein sequences share conserved properties predicted for the Nramp/Solute carrier 11 (Slc11) family, such as 12-transmembrane segment (TMS) topology (N- and C-termini cytoplasmic) and preferential conservation of four TMS predicted to form a pseudosymmetric proton/metal symport pathway. Yeast fet3fet4 mutant complementation assays showed iron transport activity for PmNramp1 and a fusion chimera of the PmNramp3 hydrophobic core and PmNramp1 N- and C-termini. PmNramp1 site-directed mutagenesis demonstrated that Slc11 invariant and predicted pseudosymmetric motifs (TMS1 Asp-Pro-Gly and TMS6 Met-Pro-His) are key for transport function. PmNramp1 TMS1 mutants D76E, G78A, and D76E/G78A prevented membrane protein expression, while TMS6 M250A, H252Y, and M250A/H252Y specifically abrogated Fe uptake; the TMS6 H252Y mutation also correlates with divergence from Nramp specificity for divalent metals.

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

Transgenic elimination of high-affinity antidepressant and cocaine sensitivity in the presynaptic serotonin transporter

Serotonin [i.e., 5-hydroxytryptamine (5-HT)]-targeted antidepressants are in wide use for the treatment of mood disorders, although many patients do not show a response or experience unpleasant side effects. Psychostimulants, such as cocaine and 3,4-methylenedioxymethamphetamine (i.e., "ecstasy"), also impact 5-HT signaling. To help dissect the contribution of 5-HT signaling to the actions of these and other agents, we developed transgenic mice in which high-affinity recognition of multiple antidepressants and cocaine is eliminated. Our animals possess a modified copy of the 5-HT transporter (i.e., SERT, slc6a4) that bears a single amino acid substitution, I172M, proximal to the 5-HT binding site. Although the M172 substitution does not impact the recognition of 5-HT, this mutation disrupts high-affinity binding of many competitive antagonists in transfected cells. Here, we demonstrate that, in M172 knock-in mice, basal SERT protein levels, 5-HT transport rates, and 5-HT levels are normal. However, SERT M172 mice display a substantial loss of sensitivity to the selective 5-HT reuptake inhibitors fluoxetine and citalopram, as well as to cocaine. Through a series of biochemical, electrophysiological, and behavioral assays, we demonstrate the unique properties of this model and establish directly that SERT is the sole protein responsible for selective 5-HT reuptake inhibitor-mediated alterations in 5-HT clearance, in 5-HT1A autoreceptor modulation of raphe neuron firing, and in behaviors used to predict the utility of antidepressants.

Y95 and E444 interaction required for high-affinity S-citalopram binding in the human serotonin transporter

The human serotonin (5-hydroxytryptamine, 5-HT) transporter (hSERT) is responsible for the reuptake of 5-HT following synaptic release, as well as for import of the biogenic amine into several non-5-HT synthesizing cells including platelets. The antidepressant citalopram blocks SERT and thereby inhibits the transport of 5-HT. To identify key residues establishing high-affinity citalopram binding, we have built a comparative model of hSERT and Drosophila melanogaster SERT (dSERT) based on the Aquifex aeolicus leucine transporter (LeuT(Aa)) crystal structure. In this study, citalopram has been docked into the homology model of hSERT and dSERT using RosettaLigand. Our models reproduce the differential binding affinities for the R- and S-isomers of citalopram in hSERT and the impact of several hSERT mutants. Species-selective binding affinities for hSERT and dSERT also can be reproduced. Interestingly, the model predicts a hydrogen bond between E444 in transmembrane domain 8 (TM8) and Y95 in TM1 that places Y95 in a downward position, thereby removing Y95 from a direct interaction with S-citalopram. Mutation of E444D results in a 10-fold reduced binding affinity for S-citalopram, supporting the hypothesis that Y95 and E444 form a stabilizing interaction in the S-citalopram/hSERT complex.

Monoamine transporters: vulnerable and vital doorkeepers

Transporters of dopamine, serotonin, and norepinephrine have been empirically used as medication targets for several mental illnesses in the last decades. These protein-targeted medications are effective only for subpopulations of patients with transporter-related brain disorders. Since the cDNA clonings in early 1990s, molecular studies of these transporters have revealed a wealth of information about the transporters' structure-activity relationship (SAR), neuropharmacology, cell biology, biochemistry, pharmacogenetics, and the diseases related to the human genes encoding these transporters among related regulators. Such new information creates a unique opportunity to develop transporter-specific medications based on SAR, mRNA, DNA, and perhaps transporter trafficking regulation for a number of highly relevant diseases including substance abuse, depression, schizophrenia, and Parkinson's disease.

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.

Design and directed evolution of a dideoxy purine nucleoside phosphorylase

Purine nucleoside phosphorylase (PNP) catalyzes the synthesis and phosphorolysis of purine nucleosides, interconverting nucleosides with their corresponding purine base and ribose-1-phosphate. While PNP plays significant roles in human and pathogen physiology, we are interested in developing PNP as a catalyst for the formation of nucleoside analog drugs of clinical relevance. Towards this aim, we describe the engineering of human PNP to accept 2',3'-dideoxyinosine (ddI, Videx((R))) as a substrate for phosphorolysis using a combination of site-directed mutagenesis and directed evolution. In human PNP, we identified a single amino acid, Tyr-88, as a likely modulator of ribose selectivity. RosettaLigand was employed to calculate binding energies for substrate and substrate analog transition state complexes for single mutants of PNP where Tyr-88 was replaced with another amino acid. In parallel, these mutants were generated by site-directed mutagenesis, expressed and purified. A tyrosine to phenylalanine mutant (Y88F) was predicted by Rosetta to improve PNP catalytic activity with respect to ddI. Kinetic characterization of this mutant determined a 9-fold improvement in k(cat) and greater than 2-fold reduction in K(M). Subsequently, we used directed evolution to select for improved variants of PNP-Y88F in Escherichia coli cell extracts resulting in an additional 3-fold improvement over the progenitor strain. The engineered PNP may form the basis for catalysts and pathways for the biosynthesis of ddI.

Practically useful: what the Rosetta protein modeling suite can do for you

The objective of this review is to enable researchers to use the software package Rosetta for biochemical and biomedicinal studies. We provide a brief review of the six most frequent research problems tackled with Rosetta. For each of these six tasks, we provide a tutorial that illustrates a basic Rosetta protocol. The Rosetta method was originally developed for de novo protein structure prediction and is regularly one of the best performers in the community-wide biennial Critical Assessment of Structure Prediction. Predictions for protein domains with fewer than 125 amino acids regularly have a backbone root-mean-square deviation of better than 5.0 Å. More impressively, there are several cases in which Rosetta has been used to predict structures with atomic level accuracy better than 2.5 Å. In addition to de novo structure prediction, Rosetta also has methods for molecular docking, homology modeling, determining protein structures from sparse experimental NMR or EPR data, and protein design. Rosetta has been used to accurately design a novel protein structure, predict the structure of protein−protein complexes, design altered specificity protein−protein and protein−DNA interactions, and stabilize proteins and protein complexes. Most recently, Rosetta has been used to solve the X-ray crystallographic phase problem.

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.

Structural analysis of the extracellular entrance to the serotonin transporter permeation pathway

Neurotransmitter transporters are responsible for removal of biogenic amine neurotransmitters after release into the synapse. These transporters are the targets for many clinically relevant drugs, such as antidepressants and psychostimulants. A high resolution crystal structure for the monoamine transporters has yet to be solved. We have developed a homology model for the serotonin transporter (SERT) based on the crystal structure of the leucine transporter (LeuT(Aa)) from Aquifex aeolicus. The objective of the present studies is to identify the structural determinants forming the entrance to the substrate permeation pathway based on predictions from the SERT homology model. Using the substituted cysteine accessibility method, we identified residues predicted to reside at the entrance to the substrate permeation pathway that were reactive with methanethiosulfonate (MTS) reagents. Of these residues, Gln(332) in transmembrane helix (TMH) VI was protected against MTS inactivation in the presence of serotonin. Surprisingly, the reactivity of Gln(332) to MTS reagents was enhanced in the presence of cocaine. Bifunctional MTS cross-linkers also were used to examine the distances between helices predicted to form the entrance into the substrate and ion permeation pathway. Our studies suggest that substrate and ligand binding may induce conformational shifts in TMH I and/or VI, providing new opportunities to refine existing homology models of SERT and related monoamine transporters.

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.

Fluoxetine (Prozac) binding to serotonin transporter is modulated by chloride and conformational changes

Serotonin transporter (SERT) is the main target for widely used antidepressant agents. Several of these drugs, including imipramine, citalopram, sertraline, and fluoxetine (Prozac), bound more avidly to SERT in the presence of Cl(-). In contrast, Cl(-) did not enhance cocaine or paroxetine binding. A Cl(-) binding site recently identified in SERT, and shown to be important for Cl(-) dependent transport, was also critical for the Cl(-) dependence of antidepressant affinity. Mutation of the residues contributing to this site eliminated the Cl(-)-mediated affinity increase for imipramine and fluoxetine. Analysis of ligand docking to a single state of SERT indicated only small differences in the energy of interaction between bound ligands and Cl(-). These differences in interaction energy cannot account for the affinity differences observed for Cl(-) dependence. However, fluoxetine binding led to a conformational change, detected by cysteine accessibility experiments, that was qualitatively different from that induced by cocaine or other ligands. Given the known Cl(-) requirement for serotonin-induced conformational changes, we propose that Cl(-) binding facilitates conformational changes required for optimal binding of fluoxetine and other antidepressant drugs.

cGMP-dependent protein kinase Ialpha associates with the antidepressant-sensitive serotonin transporter and dictates rapid modulation of serotonin uptake

Background

The Na⁺/Cl^--dependent serotonin (5-hydroxytryptamine, 5-HT) transporter (SERT) is a critical element in neuronal 5-HT signaling, being responsible for the efficient elimination of 5-HT after release. SERTs are not only targets for exogenous addictive and therapeutic agents but also can be modulated by endogenous, receptor-linked signaling pathways. We have shown that neuronal A3 adenosine receptor activation leads to enhanced presynaptic 5-HT transport in vitro and an increased rate of SERT-mediated 5-HT clearance in vivo. SERT stimulation by A3 adenosine receptors derives from an elevation of cGMP and subsequent activation of both cGMP-dependent protein kinase (PKG) and p38 mitogen-activated protein kinase. PKG activators such as 8-Br-cGMP are known to lead to transporter phosphorylation, though how this modification supports SERT regulation is unclear.

Results

In this report, we explore the kinase isoform specificity underlying the rapid stimulation of SERT activity by PKG activators. Using immortalized, rat serotonergic raphe neurons (RN46A) previously shown to support 8-Br-cGMP stimulation of SERT surface trafficking, we document expression of PKGI, and to a lower extent, PKGII. Quantitative analysis of staining profiles using permeabilized or nonpermeabilized conditions reveals that SERT colocalizes with PKGI in both intracellular and cell surface domains of RN46A cell bodies, and exhibits a more restricted, intracellular pattern of colocalization in neuritic processes. In the same cells, SERT demonstrates a lack of colocalization with PKGII in either intracellular or surface membranes. In keeping with the ability of the membrane permeant kinase inhibitor DT-2 to block 8-Br-cGMP stimulation of SERT, we found that DT-2 treatment eliminated cGMP-dependent kinase activity in PKGI-immunoreactive extracts resolved by liquid chromatography. Similarly, treatment of SERT-transfected HeLa cells with small interfering RNAs targeting endogenous PKGI eliminated 8-Br-cGMP-induced regulation of SERT activity. Co-immunoprecipitation studies show that, in transporter/kinase co-transfected cells, PKGIα specifically associates with hSERT.

Conclusion

Our findings provide evidence of a physical and compartmentalized association between SERT and PKGIα that supports rapid, 8-Br-cGMP-induced regulation of SERT. We discuss a model wherein SERT-associated PKGIα supports sequentially the mobilization of intracellular transporter-containing vesicles, leading to enhanced surface expression, and the production of catalytic-modulatory SERT phosphorylation, leading to a maximal enhancement of 5-HT clearance capacity.

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.