Discovery of novel-scaffold monoamine transporter ligands via in silico screening with the S1 pocket of the serotonin transporter.

Discovery of new inhibitors of the plasmalemmal monoamine transporters (MATs) continues to provide pharmacotherapeutic options for depression, addiction, attention deficit disorders, psychosis, narcolepsy, and Parkinson's disease. The windfall of high-resolution MAT structural information afforded by X-ray crystallography has enabled the construction of credible computational models. Elucidation of lead compounds, creation of compound structure-activity series, and pharmacologic testing are staggering expenses that could be reduced by using a MAT computational model for virtual screening (VS) of structural libraries containing millions of compounds. Here, VS of the PubChem small molecule structural database using the S1 (primary substrate) ligand pocket of a serotonin transporter homology model yielded 19 prominent "hit" compounds. In vitro pharmacology of these VS hits revealed four structurally unique MAT substrate uptake inhibitors with high nanomolar affinity at one or more of the three MATs. In vivo characterization of three of these hits revealed significant activity in a mouse model of acute depression at doses that did not elicit untoward locomotor effects. This constitutes the first report of MAT inhibitor discovery using exclusively the primary substrate pocket as a VS tool. Novel-scaffold MAT inhibitors offer hope of new medications that lack the many classic adverse effects of existing antidepressant drugs.

Current therapeutics targeting one or more of the three plasma membrane monoamine transporters (MATs), those for the substrates serotonin (SERT), norepinephrine (NET) or dopamine (DAT), carry numerous adverse effects such as gastrointestinal disorders, increased hunger, insomnia/hypersomnia, and impotence. Additionally, the therapeutic response is suboptimal for many patients.[1−6] Toward rectifying these issues, new selective serotonin reuptake inhibitors (SSRIs), serotonin–norepinephrine reuptake inhibitors (SNRIs), and serotonin–norepinephrine–dopamine triple reuptake inhibitors (TRIs) are in demand. TRIs, once viewed as useless in the process of searching for inhibitors for specific transporter proteins, are under renewed study as therapeutics not only for depression but also for obsessive-compulsive disorder, anhedonia, substance abuse, chronic pain, Parkinson’s disease, attention deficit hyperactivity disorder, autism, and obesity.[7−16] With the recent development of reliable 3D models of the plasmalemmal MAT proteins, rational discovery and design of novel MAT inhibitors should be achievable by focusing on the MAT ligand-binding pockets

At least two ligand-binding pockets, S1 and S2, are postulated to exist for the MATs (reviewed in refs (17 and 18)). Extrapolating from the position of leucine in the crystallized bacterial leucine transporter (LeuT) protein,[19] the MAT primary substrate pocket (S1) is at the approximate midpoint of the lipid bilayer and the center of the MAT transmembrane domains (Figure 1A). This positioning of the substrate is consistent with MAT mutagenesis reports.[20−22] The competitive inhibition of substrate binding by certain psychostimulants, SSRIs, and tricyclic antidepressant (TCA) drugs[22,23] suggested an inhibitor binding site that overlapped or coincided with the S1 site. Indeed, cocaine, benztropine, and various antidepressants including the SSRI fluoxetine (Prozac) and the TCA nortriptyline (Pamelor) have now been shown to reach the S1 pocket.[24−30]

Figure 1

Preparation of the SERT S1 pocket for VS. (A) Side view (lipid bilayer cross-section) of SERT computational model (gray tubes). The S1 (light gray spherical net) and S2 (approximated by the yellow spherical net) ligand-binding pockets are indicated. Bilayer interfaces with the extracellular space and cytoplasm are delineated (upper and lower dashed red lines, respectively). (B) Zoom view of the S1 pocket. Polypeptide backbones (gray tubes) of relevant TM domains are labeled. Donor projection (F1:Don2, orange spherical net), hydrophobic (F2:Hyd, blue spherical net) and volume constraint (ligand-encompassing gray spherical net) pharmacophore features are shown. Citalopram’s docking pose (carbon atoms, yellow; oxygen atoms, red; nitrogen atoms, blue; fluorine atom, green) was used to generate the pharmacophore. Distances between the Don and Hyd features, the Hyd feature and the Ile172 terminal methyl group, and the Hyd feature and the Ile172 backbone carbonyl oxygen are shown (orange dotted lines) in angstroms (orange text). Na+ ions (solid orange spheres) serving as transport cofactors are shown.

The longstanding observation that TCA drugs could also noncompetitively inhibit MAT substrate binding further implied a second, distinct MAT ligand-binding pocket.[31,32] TCA drugs cocrystallized with LeuT were bound in the protein’s extracellular vestibule, in a pocket just to the extracellular side of the four “external gate” residues that control entry into the S1 pocket.[33−35] Docking of substrates and inhibitors to the first LeuT-based MAT computational models indicated a similar position for this second binding pocket; inhibitors could progress to S1 only by application of external force, via molecular dynamics.[36,37] This vestibular pocket was suggested to be a “staging area” for substrates prior to their relocation to the S1 pocket via a MAT conformational change.[37] A similar scenario was described for binding of the MAT inhibitor cocaine.[36] Evidence for the vestibular, or S2, pocket serving as the initial substrate binding site was provided with the dopamine transporter, in which two dopamine molecules were seen to be necessary for substrate translocation through the cell membrane.[38] The S2 binding pocket has been equated with the MAT allosteric inhibitor binding site;[28] alternatively, these sites may be nonidentical, or multiple allosteric sites may exist.[39]

The necessity for, and existence of, the S2 pocket is in question. S2 as a substrate site does not appear to be necessary for LeuT,[40,41] although variability in experimental conditions has been proposed as a factor in detecting association of a second substrate molecule within the DAT.[42] Regarding the observed TCA ligand occupation of LeuT’s extracellular vestibule, the relevance of a TCA drug binding (weakly) to a bacterial leucine transporter has itself been questioned.[43] Replacement of residues in the vicinity of the S1 pocket of LeuT with their SERT counterparts (“LeuBAT”) rendered crystal structures in which the TCA clomipramine (Anafranil) occupied S1, not S2. LeuBAT crystals containing eight different SERT inhibitors covering different structural classes all localized the ligand to S1.[30] The first MAT crystal structure, a Drosophila DAT protein bound to the TCA nortriptyline, also positions the drug in S1. Similar to LeuBAT, the DAT-bound TCA drug cannot progress through the substrate pore as a substrate would because its binding extends into the region of the external gate, preventing gate closure.[27] Taken together, these findings suggest that MAT drug discovery efforts should include the S1 pocket.

Virtual screening (VS) has been successfully applied to a number of protein targets for the discovery of novel ligands.[44,45] VS employs a computational model of the drug receptor in question and involves a rapid in silico ligand docking survey of a structural library containing thousands to millions of chemical compounds. Herein, a VS hybrid approach that included both docking and structure-based pharmacophore filtering has been applied to the SERT S1 pocket, yielding SERT ligand chemotypes that one would be unlikely to find by conventional methods.

Results

Computational Model VS of a Small Molecule Structural Database for Novel SERT Ligands

Using induced-fit docking, citalopram (Celexa) was allowed to associate with the S1 pocket of the SERT model (Figure 1B). This SSRI drug, among the most SERT-preferring, has been localized to the S1 pocket[20,28,29,46] and was chosen as the template in building an S1 pocket pharmacophore. Features of the pharmacophore were based on the selected binding pose of citalopram and were added to refine the screening protocol prior to ligand VS (Figure 1C). The VS protocol was verified using an enrichment study in which 10 known non-TCA SERT ligands (Supporting Information, Figure S1) were used to seed a structural library of 1990 random compounds. (Because the evidence for TCA binding at S1 was equivocal at the time the model was optimized, TCA drugs were excluded in the 10 compound training set.) Seven of the 10 seeded compounds were among the 54 hit compounds retrieved by SERT S1 VS in screening the verification library. Following this verification step, a considerably larger structural library was screened for potential SERT ligands of novel structural scaffold.

SERT model S1 pocket screening of the PubChem database of almost half a million compounds yielded 13 378 VS hit compounds. From these, 49 were selected on the basis of visual inspection that focused on the presence of a protonatable amine, receptor placement, ligand conformation, and interactions with side chain residues. Nineteen of the 49 were found to be commercially available; these were purchased for in vitro pharmacologic characterization and labeled TN-01–TN-19 (Figure 2). All 19 VS hits contain a positively charged nitrogen atom and some aromaticity, consistent with the known SERT ligands; interestingly, only two of the 19 contain the indole ring shared with serotonin.

Figure 2

Structures of the final 19 VS hit compounds. The randomly numbered hit compounds TN-01, TN-05, TN-06, and TN-13 (boxed in red) were selected for additional pharmacologic characterization.

In Vitro Pharmacologic Characterization of VS Hit Compounds

Using the pan-specific MAT radioligand and cocaine analogue [125I]RTI-55, initial in vitro binding assays tested the ability of a single concentration (10 μM) of each nonradioactive VS hit compound in displacing the radioligand at the three plasma membrane MATs. A similar concentration of nonradioactive citalopram, mazindol, or nisoxetine served as a positive control for SERT-, DAT- or NET-selective [125I]RTI-55 binding inhibition, respectively. Depending on the transporter protein, one-quarter to one-half of the 19 VS hits displayed 50% or better inhibition of radioligand binding. Four VS hits, TN-01, TN-05, TN-06, and TN-13, displayed higher relative affinity for at least one MAT (Figure 3).

Figure 3

VS hit compound in vitro MAT binding screen. Compounds at 10 μM final concentration were tested for the ability to inhibit [125I]RTI-55 binding at hDAT N2A neuroblastoma cells (top panel), hNET N2A neuroblastoma cells (middle panel), or hSERT HEK293 cells (bottom panel). Nonspecific binding was assessed by the presence of 10 μM citalopram (CIT), mazindol (MAZ), and nisoxetine (NIS) for SERT, DAT, and NET, respectively. Data represent n = 3 independent experiments performed in duplicate. Data are presented as the mean ± SEM and were analyzed by one-way ANOVA with Dunnett’s multiple comparison posthoc test. *p < 0.01 vs total binding for that assay; ***p < 0.0001 vs total binding for that assay.

These four hit compounds were more rigorously characterized at each of the three MATs to determine binding affinities (Ki values) and substrate uptake inhibition potencies (IC50 values). Of the hit compounds, TN-05 displayed the strongest SERT affinity, with a Ki value of 668 nM. This compound had an even higher affinity at the NET (323 nM) and no detectable affinity at the DAT (Table 1). Of the four characterized compounds, TN-13 displayed the strongest NET (215 nM) and DAT (780 nM) affinities. Interestingly, this compound was a poor SERT ligand, registering a SERT affinity 50-fold weaker than that for the NET. TN-06, with an NRI profile, displayed the same rank order of affinities as that of TN-13, with modest NET affinity (841 nM) and selectivity over DAT and SERT. For TN-01 and TN-05, the NET was favored over the SERT by less than 2-fold. The TN-05 affinity for the NET was relatively high (323 nM) and no DAT affinity was detected, suggesting a SNRI classification for this hit compound. With only a slight bias toward the NET, TN-01 displayed high affinity for all three MATs, indicating potential as a lead TRI compound (Table 1).

Table 1

MAT Binding Affinities of Top VS Hit Compounds

	Ki (nM)
	SERT	DAT	NET	selectivity ratio SERT/DAT/NET
TN-01	1029 ± 81	3058 ± 403a,d	613 ± 162b,e	2:5:1
TN-05	668 ± 41	>20 000	323 ± 53c,d	2:>62:1
TN-06	>20 000	15 740 ± 2787	841 ± 225c,e	24:19:1
TN-13	12 600 ± 2122	780 ± 78a,e	215 ± 55c,e	58:4:1

SERT vs DAT.

DAT vs NET.

SERT vs NET.

p < 0.005.

p < 0.0005.

In terms of substrate uptake inhibition, the rank order of potencies for these compounds (Table 2) largely mirrored the selectivity ratios from the binding assays (Table 1). For the DAT, IC50 values for substrate uptake inhibition agreed well with binding affinity Ki values, in part because the two assays employed intact N2A cells and almost identical conditions. For the SERT and NET, binding and uptake inhibition constants did not correlate nearly as well, typically because these binding assays employed membrane preparations (a requirement for adequate radioligand binding signal-to-noise ratio).

Table 2

MAT Substrate Uptake Inhibition Potencies of Top VS Hit Compounds

	IC50 (nM)
	SERT	DAT	NET	potency ratio SERT/DAT/NET
TN-01	5025 ± 1894	2562 ± 450	2092 ± 291	2:1:1
TN-05	3845 ± 257	>20 000	635 ± 99b,d	6:>31:1
TN-06	>20 000	13 880 ± 4424	5574 ± 1607	>3:2:1
TN-13	18 870 ± 5566	720 ± 174a,c	615 ± 24b,c	30:1:1

SERT vs DAT.

SERT vs NET.

p < 0.05.

p < 0.005.

Antidepressant-Like Effects of VS Hit Compounds in Mice

Following VS identification and in vitro pharmacology, the most promising hit compounds were evaluated via the tail suspension test (TST) for the ability to induce antidepressant-like effects in mice.[47] The TST assay was first validated using citalopram and fluvoxamine (Luvox, another SSRI) as positive controls. Compared to saline-treated animals, both citalopram (10 mg/kg) and fluvoxamine (10 mg/kg) induced a significant decrease in immobility, indicative of a reduction in despair-like behavior (saline, 143.7 ± 3.3 s; citalopram, 29.7 ± 14.5 s; fluvoxamine, 59.0 ± 12.5 s; Dunnett’s multiple comparison, p < 0.001 for citalopram and fluvoxamine vs saline). Naïve mice were next acutely treated with TN-01 (0.5 mg/kg or 5 mg/kg), TN-06 (20 mg/kg), or TN-13 (10 mg/kg or 20 mg/kg). (TN-05 was not tested in vivo, as the compound was no longer commercially available at this point.) At one or more doses, all three VS hit compounds induced a statistically significant decrease in immobility compared to that of saline-treated mice (Figure 4A–C), suggesting antidepressant-like activity.

Figure 4

In vivo characterization of three VS hit compounds. (A) Naïve C57Bl/6J mice treated with TN-01 (0.5 mg/kg, n = 6 or 5 mg/kg, n = 5) showed significant decreases in immobility in the tail suspension test compared to saline-treated mice (n = 6). (p = 0.016; Dunnett’s test *p < 0.05 compared to saline.) (B) Naïve mice treated with TN-06 (1 mg/kg, n = 6; 10 mg/kg, n = 7; 20 mg/kg, n = 6) showed a significant decrease in immobility at the highest dose compared to saline (n = 6). (p = 0.018; Dunnett’s test *p < 0.05 compared to saline.) (C) Naïve mice treated with TN-13 (1 mg/kg, n = 6; 10 mg/kg, n = 6; 20 mg/kg, n = 6) showed a significant decrease in immobility with the two highest doses compared to saline (n = 6). (p = 0.0064; Dunnett’s test *p < 0.05, **p < 0.01 compared to saline.) The distance traveled in an open field compared to saline-treated mice (n = 5) is shown for mice treated with (D) TN-01 (5 mg/kg, n = 6), (E) TN-06 (20 mg/kg, n = 5), or (F) TN-13 (1 mg/kg, n = 3; 10 mg/kg, n = 3).

To test for possible nonspecific motor effects, an open field test was conducted following VS compound injection. The open field test allows for analysis of total locomotion as well as measures of anxiety (time in the center compared to the outside edge of the open field). Neither TN-01 (5 mg/kg), TN-06 (20 mg/kg), nor TN-13 (10 mg/kg) induced statistically significant changes in the total distance traveled compared to that of saline-treated mice (Figure 4D–F), nor was there a statistically significant change in center time between vehicle groups and any of the three VS hit compounds (Supporting Information, Table S2). These data indicate that nonspecific locomotor effects of the compounds cannot account for the TST results.

Discussion

Previously, this research group created a computational model of the human SERT based on the X-ray crystal structure of the LeuT “outward-occluded” conformation (PDB: 2a65),[19] within which the putative S2 substrate binding pocket in the SERT extracellular vestibule was used to find new MAT ligands via VS.[49] Here, the S1 primary substrate binding pocket of the same SERT homology model was employed to screen the PubChem small molecule structural library for SERT inhibitor candidates of atypical scaffold. Access of these inhibitors to the S1 site had not been permitted using MAT computational models unless the protein was afforded enough flexibility to accommodate induced-fit ligand docking.[50] Only recently have substrate- or inhibitor-free LeuT crystal structures become available in which the protein presents an S1-accessible conformation.

The VS approach was primarily employed to uncover potential MAT ligands of novel scaffold (i.e., those that are not clearly structure–activity relatives of established MAT ligands). All hits contained two conserved features: a protonated amine and a hydrophobic moiety ∼6.5 Å away. Three of the four characterized hit compounds, TN-01, TN-05, and TN-06, have scaffolds differing from other reported MAT ligands. While TN-13 is structurally related to a class of compounds recently reported as SERT imaging agents,[51−56] the structural differences confer a unique selectivity profile. Interestingly, TN-13 preferred NET and DAT, whereas the imaging agents were selective for SERT, suggesting that TN-13’s extended amine linker and lack of substituents are important for NET/DAT activity (Figure 5). The lack of selectivity for SERT was not unexpected given the high sequence similarity among the three transporters at the S1 binding site.

Figure 5

Docking placement of the top four VS hits (in yellow) in the S1 SERT pocket: (A) TN-01, (B) TN-05, (C) TN-06, and (D) TN-13. Polar hydrogen atoms are shown. Ligand atoms are color-coded (carbon atoms, yellow; oxygen atoms, red; nitrogen atoms, blue; fluoride atoms, light green; chloride atoms, dark green; sulfur atoms, pink). Relevant SERT side chains are annotated (carbon atoms, gray; oxygen atoms, red; nitrogen atoms, blue). Na+ ions serving as transport cofactors (solid orange spheres) are positioned consistent with the LeuT crystal structure, with the Na+ in the background hidden for clarity.

Ultimately, the goal is to use VS techniques to identify lead compounds with in vivo efficacy in treating depression and other monoamine-related disorders. Depression is associated with a variety of behavioral changes including loss of energy, sleep alterations, learning and memory impairments, anhedonia, and helplessness (despair). Of the four VS hits characterized in vitro for binding affinity and uptake inhibition potency, TN-01, TN-06, and TN-13 were used with naïve mice in the tail suspension test (TST). Reduced immobility (or increased movement) in the TST is associated with a decrease in despair-like behavior. Behavioral effects in the TST are often seen after acute treatment with SSRIs, even though chronic treatment is often necessary for efficacy in human patients.[57] This apparent disconnect may be resolved by evidence of acute[58] and subchronic[59] effects of SSRIs in human subjects. In fact, the TST has shown strong predictive validity in mice as a screen for efficacious human antidepressants.[57] The TST and locomotor data suggest that TN-01, TN-06, and TN-13 have antidepressant activity but not anxiolytic activity and that the antidepressant action is not confounded by nonspecific changes in locomotor output after treatment. The lack of an anxiolytic phenotype despite the presence of antidepressant-like activity is similar to that of established SSRIs such as fluoxetine.[60] Future analysis of these VS hit compounds will include tests for chronic antidepressant efficacy, effects on learning and memory, and additional assays for anxiety-like behavior.

Until recently, only the extracellular vestibule region of MAT computational models was employed for VS.[49,61] The first report of pharmacology-supported MAT structure-based VS exclusive to the S1 pocket was a drug-repurposing effort using a NET model. The VS yielded five novel NET inhibitors with moderate micromolar affinity; notably, all five compounds were structurally similar to that of the norepinephrine substrate.[65] More recently, the SERT S1 pocket was used as a VS tool to identify two novel compounds proposed to possess better SERT affinity than that of paroxetine (Paxil), a classic SSRI with subnanomolar SERT affinity. The affinity values, however, were extrapolated from computational modeling predictions as opposed to being pharmacologically verified.[62] Very recently, SERT structure-based VS has been employed utilizing an outward-facing (extracellular-facing) SERT model conformation[63] that allowed simultaneous access to the S1 and S2 pockets; thus, both pockets and the extracellular vestibule served as potential hit compound binding sites. Several VS hits were obtained that displayed nanomolar to low-micromolar Ki values and a degree of structural uniqueness; SERT selectivity was not addressed.[64]

Here, the discovery of novel MAT ligands through a hybrid VS approach has been described. Specifically, the screening of a large small molecule structural database using exclusively the S1 binding site of a SERT homology model has afforded four submicromolar affinity hits with varying MAT selectivity profiles. The hit compounds were confirmed to have true MAT activity, as measured by inhibition of substrate uptake, as opposed to merely having MAT binding affinity. Three of these compounds show antidepressant-like activity in a rodent model. This represents the first report in which antidepressant candidate compounds have been identified in silico based solely on the primary ligand-binding (S1) pocket of the SERT and validated using in vitro and in vivo pharmacology. These compounds may now serve both as leads for the discovery of new MAT therapeutics via SAR guided by the computational SERT model as well as new tools in investigating MAT mechanism of action. With such reliable computational models in hand, the VS approach to drug discovery is accessible to research universities as well as to the pharmaceutical industry, which lacks the cost-prohibitive high-throughput in vitro screening steps characteristic of classic pharmaceutical development.

Methods

Materials

Molecular modeling was performed using MOE v2010 and v2011.10 software from Chemical Computing Group (Montreal, Quebec, CA). The radioligands [3H]serotonin (∼28 Ci/mmol), [3H]dopamine (∼26 Ci/mmol), [3H]norepinephrine (∼26 Ci/mmol), and [125I]RTI-55 (∼2200 Ci/mmol) were obtained from PerkinElmer Life and Analytical Sciences (Foster City, CA). Nonradioactive citalopram, mazindol, nisoxetine, and fluvoxamine were obtained from Tocris Bioscience (Ellisville, MO). Virtual screening hit compounds were purchased from Ambinter (Orleans, FR). C57BL/6J mice were obtained from The Jackson Laboratory (Bar Harbor, ME). Data analysis was performed using GraphPad Prism 5.0 (GraphPad Software, San Diego, CA).

Molecular Modeling

Database Generation

The PubChem database, consisting of 473 965 compounds, was downloaded from www.ncbi.nlm.nih.gov. Three-dimensional structures were generated, and partial charges were determined using the MMFF94x force field in MOE. The database was then “washed” to remove salt fragments and metals, followed by generation of a maximum of 30 tautomers for each compound. Protonation states of each compound were considered, protonating strong bases and deprotonating strong acids. Using the stochastic search method, a conformational search was carried out on the modified database, now containing 1 091 982 entries, to generate a maximum of 30 conformations for each compound. During this stage, several filters were employed to eliminate undesirable compounds and to retain compounds with the following characteristics: MW ≤ 600, donor/acceptor atoms ≤ 12, chiral centers ≤ 4, rotatable bonds < 7, and a LogP range of −4 to 5.

Docking Site Preparation

A previously described SERT model was modified for use in this study.[49] An initial docking study of known SSRI and TCA compound SERT ligands revealed the inability of either class to comfortably dock into the S1 site of the occluded model (LeuT PDB: 2a65;[19]). On the basis of the literature, citalopram was assumed to bind in the S1 site and was therefore used as the reference compound for an induced-fit docking in an effort to expand the S1 pocket (Figure 1B).[20,46] Protonate 3D was used to prepare the protein for docking. Alpha spheres selected using SiteFinder defined the S1 site for VS docking. This site corresponds to the position of serotonin and citalopram binding. Structure-based pharmacophore features were introduced (Figure 1C) in order to retain only hits capable of matching three features: a volume constraint (radius = 6 Å, centered on citalopram docked into S1), a hydrogen bond donor projection (radius = 1 Å, centered on a carboxylate oxygen of Asp-98 in TM 1), and a hydrophobic element (radius = 1.4 Å, placed within 2 Å of Ile-172 in TM 3; measurements in angstroms shown in Figure 1C).

Enrichment Docking Study

An enrichment database containing 1990 unknown compounds seeded with 10 known SERT ligands including SSRIs and SNRIs (Table S1, Supporting Information) was generated to evaluate the ability of the employed screening parameters to retrieve known ligands from a database. Several trials of docking in MOE were carried out in order to fine tune the parameters to be used for the actual screening. Additionally, the volume constraint was adjusted in order to limit the number of hit compounds while still allowing known compounds to be selected. A final protocol consisted of the Proxy Triangle placement, the Affinity dG scoring function, and a volume constraint with a 6 Å radius, all features of MOE software. Compounds capable of docking into the S1 site and matching the above pharmacophore features were considered to be hits.

While the score proved to be useful for setting a cutoff limit for hit compounds, it was not capable of rank-ordering known ligands with respect to actual experimental binding affinities. The final enrichment docking retrieved 6 of the 10 known SERT ligands out of a total of 253 hits.

PubChem Database Screening

Using the above docking protocol, the 10 subsets of the PubChem database were screened, resulting in 13 378 hit compounds. As mentioned, the scoring function was used only as a cutoff limit to select compounds for visual inspection, not for ranking. Compounds with S < 0 and MW = 200–450 were inspected both for their fit with the pharmacophore as well as for the interactions formed with the protein. Chemical complexity included the number of stereocenters; synthetic feasibility was taken into consideration. On the basis of these criteria, 49 compounds were selected, of which 19 were commercially available.

In Vitro Pharmacology

In Vitro Pharmacologic Screening

The 19 hit compounds yielded by the S1 SERT virtual screening were purchased and initially tested at a single concentration (10 μM) with cells expressing one of the three MATs, toward detecting specific binding. hDAT or hNET N2A whole cells or hSERT HEK membranes were coincubated with VS hit compound and [125I]RTI-55 (∼0.1 nM concentration). Four hit compounds capable of inhibiting >50% of the radioligand binding were examined further.

SERT Membrane Binding

Membranes were prepared using SERT HEK stable cells grown in a 5% CO2 environment. Cell monolayers were washed twice with 10 mL cold phosphate-buffered saline (DPBS). An additional 10 mL cold DPBS was added, and cells were scraped from the plate, transferred to 15 mL tubes, and centrifuged at low speed (700g). Supernatant was removed followed by resuspension of the cell pellet in 500 μL cold TE buffer (50 mM Tris, pH 7.5, 1 mM EDTA). Homogenate was transferred to cold 1.5 mL microcentrifuge tubes and centrifuged at 100 000g at 4 °C for 30 min (Sorvall Discovery M150 centrifuge). Supernatant was discarded, and the pellet was frozen for later use or resuspended in cold binding buffer (50 mM Tris, pH 7.5, 100 mM NaCl) for immediate use in a membrane-binding assay. Each sample was analyzed for protein content using the micro-Bradford protein assay (Bio-Rad). For competition binding, membrane fractions were incubated with [125I]RTI-55 (∼0.1 nM concentration) radioligand and increasing concentrations of nonradioactive competitor (1 fM to 1 μM concentration) or 10 μM citalopram for nonspecific binding. Reactions were carried out in 12 × 75 mm borosilicate glass tubes at 22 °C for 1 h and terminated by rapid filtration through GF/B filters (Schleicher and Schuell, Keene, NH) presoaked in 0.5% polyethylenimine solution (v/v). Filters were washed twice with 5 mL of cold 50 mM Tris buffer and transferred to counting vials. Radioactivity was determined using a Beckman gamma counter.

DAT and NET Competition Binding

Competition binding assays were performed using hDAT- or hNET-expressing N2A cells grown on 24-well plates in a 5% CO2 environment. Cell monolayers were initially washed twice with 1 mL of KRH buffer (25 mM HEPES, pH 7.3, 125 mM NaCl, 4.8 mM KCl, 1.3 mM CaCl2, 1.2 mM MgSO4, 1.2 mM KH2PO4, 5.6 mM glucose) supplemented with 50 μM ascorbic acid (KRH/AA). Cells were then incubated for 15 min with [125I]RTI-55 (0.1 nM concentration) and either increasing concentrations of nonradioactive competitor (1 fM to 1 μM concentration) or 10 μM of mazindol or nisoxetine for assessing nonspecific binding of DAT or NET, respectively. Cell monolayers were washed twice with 1 mL of KRH/AA buffer and were then treated with 1 mL of 1% SDS with gentle shaking at room temperature for 1 h. Cell lysates were transferred into 5 mL scintillation fluid for radioactivity analysis using a liquid scintillation analyzer.

DAT, NET, and SERT Uptake Inhibition

[3H]substrate uptake assays were performed using hDAT- or hNET-expressing N2A cells or hSERT HEK cells grown on 24-well plates in a 5% CO2 environment. hSERT HEK cells were grown on poly l-lysine coated plates to enhance cell adhesion. Cell monolayers were washed twice with 1 mL of KRH/AA buffer, followed by preincubation for 10 min with either increasing concentrations of the drug of interest or 10 μM mazindol, nisoxetine, or citalopram for assessment of nonspecific uptake by DAT, NET or SERT, respectively. Cells were then treated with 10 nM [3H]substrate (dopamine, norepinephrine, or serotonin) supplemented with 99 μM tropolone (total volume of 500 μL) for 5 min. Reactions were terminated by washing twice with 1 mL of KRH/AA buffer to remove any remaining substrate radioligand from the extracellular milieu. Cells were then lysed with 1 mL of 1% SDS under gentle shaking at 22 °C for 1 h, after which cell lysates were transferred into 5 mL scintillation fluid tubes for radioactivity analysis using a liquid scintillation analyzer.

Data Analysis

Experimental data, expressed as counts per minute (cpm), were analyzed with GraphPad Prism. Nonspecific binding was subtracted, and data were transformed to percentages with respect to baseline levels (“no drug”). Nonlinear regression and one-site analysis were used to determine binding (Ki) and uptake (IC50) values.

Behavioral Testing

Animals

All mouse protocols were in accordance with National Institutes of Health guidelines and were approved by the Animal Care and Use Committee of Duquesne University (Pittsburgh, PA). Male C57Bl/6J mice were group-housed on a 12 h/12 h light/dark cycle with ad libitum access to rodent chow and water.

Behavioral Analysis

All behavioral analyses were performed by an observer blinded to treatment. Behavioral tests were conducted with adult male mice 8–20 weeks of age. Mice were treated with citalopram (10 mg/kg in 0.9% normal saline), fluvoxamine (10 mg/kg in 0.9% normal saline), TN-01 (0.5 or 5 mg/kg in 0.9% normal saline), TN-06 (1, 10, or 20 mg/kg in 0.9% normal saline), TN-13 (1, 10, or 20 mg/kg in 0.9% normal saline), or vehicle (0.9% normal saline, pH 7.5). Treatment was given 30 min prior to testing via intraperitoneal injection (0.3 mL). The highest safe dose of each compound showing efficacy in the tail suspension test was chosen for use in the open field test.

Tail Suspension Test (TST)

The TST apparatus consisted of a cubicle made of 1.2 cm Plexiglas with inside dimensions of 40 (w) × 40 (l) × 35 (d) cm3. Mice were suspended by the distal 1.5 cm of their tails with tape. Activity was continuously scored for immobility behavior during the entire 6 min trial. Immobility was defined as the lack of all motion except respiration.

Open Field Test (OFT)

Lighting was provided by a single 100 W incandescent light bulb placed 2 m above the Plexiglas box. Each mouse was placed in a corner of the box and subjected to one 10 min trial. Between sessions, the box was rinsed with 70% ethanol and dried with paper towels. Total distance traveled and time in center square (31 × 31 cm2) were analyzed using Any-Maze software (Stoelting).

Experimental data, expressed as mean ± SEM, were analyzed with GraphPad Prism. All behavior was analyzed with Student’s t test or one-way ANOVA followed by Dunnett’s multiple comparison tests. p < 0.05 was considered statistically significant.