Pharmacological Characterization of a Betaine/GABA Transporter 1 (BGT1) Inhibitor Displaying an Unusual Biphasic Inhibition Profile and Anti-seizure Effects

Focal epileptic seizures can in some patients be managed by inhibiting γ-aminobutyric acid (GABA) uptake via the GABA transporter 1 (GAT1) using tiagabine (Gabitril®). Synergistic anti-seizure effects achieved by inhibition of both GAT1 and the betaine/GABA transporter (BGT1) by tiagabine and EF1502, compared to tiagabine alone, suggest BGT1 as a target in epilepsy. Yet, selective BGT1 inhibitors are needed for validation of this hypothesis. In that search, a series of BGT1 inhibitors typified by (1R,2S)-2-((4,4-bis(3-methylthiophen-2-yl)but-3-en-yl)(methyl)amino)cyclohexanecarboxylic acid (SBV2-114) was developed. A thorough pharmacological characterization of SBV2-114 using a cell-based [3H]GABA uptake assay at heterologously expressed BGT1, revealed an elusive biphasic inhibition profile with two IC50 values (4.7 and 556 μM). The biphasic profile was common for this structural class of compounds, including EF1502, and was confirmed in the MDCK II cell line endogenously expressing BGT1. The possibility of two binding sites for SBV2-114 at BGT1 was assessed by computational docking studies and examined by mutational studies. These investigations confirmed that the conserved residue Q299 in BGT1 is involved in, but not solely responsible for the biphasic inhibition profile of SBV2-114. Animal studies revealed anti-seizure effects of SBV2-114 in two mouse models, supporting a function of BGT1 in epilepsy. However, as SBV2-114 is apparent to be rather non-selective for BGT1, the translational relevance of this observation is unknown. Nevertheless, SBV2-114 constitutes a valuable tool compound to study the molecular mechanism of an emerging biphasic profile of BGT1-mediated GABA transport and the putative involvement of two binding sites for this class of compounds. Electronic supplementary material The online version of this article (10.1007/s11064-020-03017-y) contains supplementary material, which is available to authorized users.


Animals
Male audiogenic-susceptible (AGS) Frings mice were obtained from an in-house breeding facility at the University of Utah (Salt Lake City, UT, USA). Male CF-1 mice for the maximal electroshock (MES) test and the subcutaneous pentylenetetrazol (s.c.PTZ) seizure threshold test were purchased from Charles River (Kingston, WA, USA). All mice were housed in a temperature, humidity, and light (on at 6:00 AM) controlled facility at the University of Utah with access to ad libitum food and water except during testing. All procedures were carried out at the University of Utah in concordance with the National Institutes of Health Guidelines for the Care and Use of Laboratory Animals and under approval from the University of Utah's Institutional Animal Care and Use Committee. The mice were euthanized in accordance with Public Health Service policies on the humane care of laboratory animals.

Audiogenic Seizures
The ability of SBV2-114 to prevent sound-induced seizures was studied in AGS Frings mice. The seizures were induced as previously described [20]. In brief, Frings mice were placed in a Plexiglas cylinder (diameter 15 cm, height 18 cm) connected to an audio transducer and exposed to a sound stimulus of 110 decibel (11 kHz) for a duration of 20 s. SBV2-114 (50 and 120 mg/kg, n = 4) was injected intraperitoneal (i.p.) at various time points to determine the time-to peak (TPE). The mice were considered protected if they failed to display full hind limb tonic extension seizures. All Frings mice were pre-screened the day prior to testing and only those mice that displayed a tonic extension seizure in response to the sound stimulus were included in the SBV2-114 studies.

Maximal Electroshock (MES) Test
The ability of SBV2-114 to prevent electroshock-induced seizures was assessed using the MES test. This is a model of generalised tonic-clonic seizures, and was conducted as described previously [21]. In brief, a drop of anaesthetic/ electrolyte solution (0.5% tetracaine hydrochloride in 0.9% saline) was applied to the cornea of each animal prior to placement of the silver-coated corneal electrodes. An alternating current of 50 mA (60 Hz) was delivered for 0.2 s to the cornea of the mouse using a Woodbury-Davenport stimulator. The supramaximal current applied to the cornea was sufficient to induce maximal tonic extension seizures in all non-treated mice. Separate groups of mice were treated with increasing doses of SBV2-114 (50-200 mg/kg, i.p., n = 3-7) and were exposed to the MES stimulus at the predetermined SBV2-114 TPE; i.e., 15 min. Those mice not displaying full hind limb tonic extension seizures were considered protected.

Subcutaneous Pentylenetetrazol Seizure Test
SBV2-114 was tested for its ability to prevent a minimal clonic seizure in the s.c.PTZ seizure test, which is a test of generalized myoclonic seizures. The s.c.PTZ seizure test was carried out as previously described [21], and SBV2-114 was injected i.p. at a dose of 150 mg/kg (n = 8). At the TPE, mice were challenged with a s.c. dose of PTZ (85 mg/kg) and observed for 30 min for the presence or absence of a minimal clonic seizure characterized by minimal forelimb clonus and/or vibrissae twitching. Mice not displaying this phenotypic behaviour were considered protected.

Behavioural Impairment
Minimal motor impairment of mice receiving SBV2-114 was evaluated using the rotarod procedure [43]. Impairment was noted if the mice failed to maintain their equilibrium for 1 min in three consecutive trials on a one-inch rotating rod (6 rpm). The rotarod test was performed just prior to the three seizure tests.

Quantification of Anti-seizure Efficacy and Behavioural Impairment
Probit analysis [21] was used to analyse the efficacy and tolerability results obtained from the dose-response studies and to calculate the median effective (ED 50 ) dose, the median toxic dose (TD 50 ), and the 95% confidence interval.

Screening for Selectivity
The selectivity of SBV2-114 for BGT1 for 42 different neuroreceptors and transporters was tested by examining the affinity of SBV2-114 in radioligand binding assays at the National Institute of Mental Health's Psychoactive Drug Screening Program (NIMH PDSP, Contract # HHSN-271-2008-00025-C). SBV2-114 was tested at a concentration of 10 μM followed by full curves if inhibition exceeded 50% at 10 μM.

Statistical Analysis
All statistical and pharmacological analyses were performed using GraphPad Prism (version 8.2.1, Yosemite, GraphPad Software, San Diego, CA, USA), and appropriate statistical tests are specified accordingly. Level of significance was set at P < 0.05.

Identification of a Biphasic Inhibition Profile of SBV2-114 at BGT1
Previous studies in our group have revealed that SBV2-114 has a potency in the mid μM range for mBGT1 and displays about 26, 43 and 13 times reduced inhibitory potency at mGAT1, mGAT2 and mGAT3, respectively [30]. In this study, a thorough pharmacological characterization of SBV2-114 using the competitive [ 3 H]GABA uptake assay was performed. This entailed 12 concentrations to generate a detailed CRC. This study revealed that SBV2-114 inhibited mBGT1 (stably expressed in HEK-293 cells) and hBGT1 (transiently expressed in tsA201 cells) in an unusual biphasic manner (confirmed by F tests) with IC 50 values in the low (m/hBGT1: 3.8/4.7 μM) and high μM range (m/hBGT1: 402.7/555.9 μM) corresponding to the high and low affinity components, respectively ( Fig. 2a and Table 1). The % inhibition (Frac value) at the low IC 50 was 45% at mBGT1 and 42% at hBGT1. This biphasic behaviour was not seen in previously performed 8-point mBGT1 inhibition CRCs [30]. In comparison, the endogenous ligand GABA, the BGT1 inhibitor BPDBA (Fig. 2a, Supplemental Fig. 2b), and the conventional GAT1 inhibitor tiagabine all inhibit [ 3 H]GABA uptake in a clearly monophasic manner [36]. Furthermore, the IC 50 values of SBV2-114 were in the same range regardless of the amount of plasmid DNA used for the transient expression of hBGT1 (Supplemental Table 1) indicating that the biphasic inhibition profile at BGT1 was not influenced by the expression level of BGT1.
To investigate if the biphasic inhibition profile of SBV2-114 was a consequence of recombinant overexpression in  Table 1 for details). Typical counts ranged from 100 to 8000 CPM at hBGT1 and mBGT1 and 500 to 2000 CPM at canine BGT1 the heterologous expression system, SBV2-114 was tested at MDCK II cells endogenously expressing canine BGT1 [34]. In this study with MDCK II cells, SBV2-114 retained its biphasic inhibition profile (P = 0.0067, F test) at BGT1 with IC 50 values of 0.1 μM and 250.6 μM ( Fig. 2b and Table 1) supporting the finding that the biphasic inhibition profile of SBV2-114 at BGT1 is not a phenomenon only seen in overexpressed cells and is independent of the species.
To ascertain whether the biphasic CRC of SBV2-114 at BGT1 was caused by kinetic factors, e.g. by a slower onkinetics, tsA201 cells expressing hBGT1 were pre-incubated for 10 min with SBV2-114 prior to the [ 3 H]GABA uptake assay. This had no impact on the curve fitting and the biphasic fit was still preferred (P < 0.0001, F test, data not shown) although pre-incubation reduced the % inhibition at the low IC 50 as the Frac value decreased from 0.45 to 0.34.
To test whether the biphasic inhibition profile was confined to BGT1, the detailed pharmacological characterization of SBV2-114 was extended to mGAT1-3 stably expressed in HEK-293 cells [32] and hGATs stably expressed in Flp-In CHO cells [36]. It is clear that SBV2-114 did not inhibit the m/hGAT1-3 in a typical monophasic manner (Supplemental Tables 2 and 3). Specifically, the fitting of the inhibition profile of SBV2-114 was predominantly ambiguous, meaning neither fitted to a monophasic or biphasic profile, at mGAT1-3-HEK-293 with IC 50 values > 1000 μM (Supplemental Table 2). SBV2-114 also produced ambiguous fitting results at the hGATs (Supplemental Table 3). In detail, SBV2-114 inhibited the [ 3 H]GABA uptake at hBGT1-CHO in a biphasic manner in two out of four experiments, and with a reduced % inhibition at the low IC 50 (Frac = 0.10). Additionally, SBV2-114 displayed a biphasic inhibition profile at hGAT3, but also with a reduced Frac value compared to hBGT1-tsA, and in one out of four experiments at GAT2 (Supplemental Table 3). These data collectively show that SBV2-114 generally inhibits the [ 3 H]GABA uptake by m/ hGATs in an atypical manner.

Several Structurally Related SBV2-114 Analogues Also Inhibit BGT1 in a Biphasic Manner
To investigate whether the biphasic inhibition profile at BGT1 is unique to SBV2-114, we also characterized the inhibition profile of a number of different compound classes in the [ 3 H]GABA uptake assay ( Table 2). We tested three compounds that all contain a β-amino acid moiety like SBV2-114, but with varying stereochemistry and different di-aromatic lipophilic side-chains including (1R,2R)-4a, (1S,2R)-5a, and (1S,2R)-5c [30]. (R)-EF1502 was also included in the detailed pharmacological characterization, with no prior notion of a biphasic inhibition profile, as it is structurally highly similar to SBV2-114. Finally, two BGT1 inhibitors belonging to other compound classes namely BPDBA [31] and NNC 05-2090 [44] were included (Fig. 1). Convincingly, all compounds containing a β-amino acid moiety with di-aromatic lipophilic side-chains including (R)-EF1502 showed a biphasic CRC at hBGT1 transiently expressed in tsA201 cells. The most potent component of the biphasic response ranged from 2.5-11.2 μM inhibiting a fraction of 0.19-0.42, where the largest inhibition was observed with SBV2-114. The high IC 50 values ranged from 190.5 to 818.5 μM (Table 2). (R)-EF1502 and SBV2-114 were the compounds with the highest F value, respectively, and hence statistically, the ones with the most distinct biphasic CRCs. BGT1 inhibitors, not comprising a β-amino acid moiety attached to di-aromatic lipophilic side-chains, such as BPDBA and NNC 05-2090, did not inhibit hBGT1 in a biphasic manner ( Table 2, Supplemental Fig. 2b).

Investigations into the Possible Mechanisms of Action
Having verified the biphasic response of SBV2-114 and structurally related compounds at BGT1 in the [ 3 H]GABA uptake assay, we wanted to address the molecular details. To this end, the possibility of SBV2-114 binding to two Table 1 Inhibitory activity of SBV2-114 at mouse, human and canine BGT1. SBV2-114 was tested for its ability to inhibit uptake of [ 3 H]GABA (see Fig. 1

for details)
All experiments were performed in triplicate in four-five independent experiments. The data fitting is based on the preferred model according to the extra-sum-of-squares F test, and the associated P, F, DFn, DFd, Frac values, and nH 1 , nH 2

Could SBV2-114 Bind to Two Non-equal Binding Sites?
The fact that SBV2-114 clearly inhibits BGT1 in a biphasic manner with two distinct IC 50 values could suggest that SBV2-114 possibly binds to two different sites in BGT1. This is in line with recent crystallographic evidence that the antidepressant escitalopram not only binds to the orthosteric pocket of the human SERT, but also to an allosteric pocket in the extracellular vestibule above the orthosteric site [15]. Since hSERT is highly related to hBGT1 (44% sequence identity), both belong to the SLC6A family as well as share the same fold, it is very likely that a similar allosteric binding site exist in BGT1 and the other GATs [13,15,37,46,47].
To address the possibility of two binding sites for SBV2-114 at BGT1, SBV2-114 was docked into the orthosteric and allosteric site of a hBGT1 homology model based on the hSERT crystal structure PDB ID 5I73 co-crystallized with escitalopram in the orthosteric and allosteric site [15,37]. After visually investigating the top scored poses of The compounds were tested for their ability to inhibit uptake of 30 nM [3H]GABA for 3 min. All experiments were performed in triplicates in three-five independent experiments. The curve fittings are based on the preferred model according to the extra-sum-of-squares F test with the associated P, F, DFn, DFd, Frac values, and nH 1 , nH 2 indicated accordingly. Non-significant (n.s.) SBV2-114 in the ortho and allosteric site, we conclude that it is sterically possible for SBV2-114 to bind to both sites simultaneously (Fig. 3a). In more detail, two potential binding modes for SBV2-11 in the orthosteric site as well as one possible binding mode in the allosteric site was identified (Fig. 3a). The first binding mode in the orthosteric site is characterized by the carboxyl group of SBV2-114, which is a key feature of GAT inhibitors, coordinating the sodium ion (Na 1) (Fig. 3b). The positively charged nitrogen of SBV2-114, which is the second key feature of GAT inhibitors, forms a hydrogen bond with the backbone of Phe293. This residue is part of the extracellular lid which is essential for the translocation process and therefore also highly conserved among the neurotransmitter transporters [13,48]. Finally, the lipophilic di-aromatic side chain of SBV2-114 expands into the primarily lipophilic extracellular vestibule. Since all interactions observed were either formed by functionally essential residues or ions, this binding mode was not accessible to validation by mutational experiments.
Interestingly, among the top scored poses we also observed a second binding mode in the orthosteric site where the carboxyl group of SBV2-114 forms a hydrogen bond with the side chain of glutamine 299 (Q299) and the protonated nitrogen undergoes hydrogen bonding with Ser395 (Fig. 3c). A similar alternative binding mode was already discussed by Vogensen et al. [30]. The lipophilic aromatic side chain extends analogously to the first binding mode into the extracellular vestibule. While Ser395 constitutes a highly conserved residue amongst the GATs and the monoamine transporter, Q299 is a unique residue for BGT1, corresponding to leucine in all other GATs (see alignment [46]). An interaction with this residue could possibly explain subtype selectivity and be important for the observed biphasic inhibition profile. This prompted us to investigate the activity of SBV2-114 in the hBGT1 mutant where we would predict a loss of activity. However, the [ 3 H] GABA uptake results for the hBGT1 Q299L mutant showed a clear biphasic inhibition profile and no significant altered activity for SBV2-114, which would contradict the second binding mode and the relevance of Q299 for binding (data not shown). Additionally, we also investigated the activity of SBV2-114 in the corresponding hGAT3 mutant L314Q mimicking the orthosteric Q299 site of hBGT1. Through this mutation we were able to strengthen the biphasic inhibition profile of SBV2-114 and introduce activity in hGAT3 L314Q similar to hBGT1 underlining the potential relevance of Q299 for hBGT1 activity. In detail, the profile of SBV2-114 at GAT3 L314Q had a higher F value compared to at hGAT3 (15.29 vs 8.14, respectively), and hence, statistically a more pronounced biphasic CRC, and a higher Frac value compared to both hGAT3 and hBGT1 (Fig. 4).
In terms of the allosteric site (third binding mode), the best poses showed that the carboxylic group undergoes hydrogen bonding with the side chains of tyrosine 453 (Tyr453), arginine 61 (Arg61) and glutamine 290 (Gln290), whereas the protonated nitrogen forms a hydrogen bond with aspartate 452 (Asp452) (Fig. 3d). Only Tyr453 has been proposed by others to be relevant for subtype selectivity [37,47]. Because this residue is also the only interacting residue not highly conserved among the GATs and monoamine transporters, we investigated the inhibition profile of SBV2-114 at the hBGT1 Y453A mutant. However, the [ 3 H] GABA uptake results for the hBGT Y453A mutant did not show the expected decrease in activity for SBV2-114 (data not shown).

SBV2-114 at High Concentrations Leads to an Intracellular Ca 2+ Increase
A link between changes in intracellular Ca 2+ levels and altered GAT1-and GAT3-mediated [ 3 H]GABA uptake has been reported [45], which prompted us to investigate if SBV2-114 triggers an intracellular Ca 2+ increase in tsA201 cells transiently expressing hBGT1. As depicted in Fig. 5a, SBV2-114 triggered an intracellular Ca 2+ increase at 200 μM and 500 μM within the 3 min timeframe used for [ 3 H]GABA uptake at hBGT1. As a positive control, ATP also induced a Ca 2+ increase. Both effects could be blocked by pre-incubation with the membrane permeable Ca 2+ chelator BAPTA/ AM (Fig. 5a). A similar Ca 2+ response was observed in  hGAT1 (b). SBV2-114 was tested for its ability to induce an intracellular Ca 2+ increase determined using the Fluo-4 Ca 2+ -assay with and without the presence of 25 μM BAPTA/ AM. The experiments were performed in triplicate in three-four independent experiments and depicted as mean ± S.E.M. One-Way ANOVA followed by Dunnett's multiple comparison test comparing to control levels, **P < 0.01, ***P < 0.001 tsA201 cells expressing hGAT1 (Fig. 5b), revealing that the effect was not BGT1 specific, although it is unknown if SBV2-114 triggers a Ca 2+ response in tsA201 cells from internal stores or influx and through which target(s).
Next, we wanted to investigate if the biphasic CRC of SBV2-114 was dependent on the Ca 2+ response, and where the Ca 2+ originates from. To this end, tsA201 cells expressing hBGT1 were pre-incubated with either the cell impermeant Ca 2+ chelator BAPTA to abolish any possible Ca 2+ influx or the cell permeant BAPTA/AM to chelate all intracellular Ca 2+ increases in response to SBV2-114. It is clear from these experiments that the pre-incubation itself (30 min at 37 °C), both in the assay buffer (control) and assay buffer supplemented with BAPTA/AM, affected the inhibition curve of SBV2-114 at hBGT1 in the [ 3 H] GABA uptake assay (Fig. 6), however without changing the IC 50 values. Specifically, lowered Frac values and Hill coefficients compared to Table 1 were observed both in the control (IC 50 = 6.1 μM and 507.0; μM, nH 1 = -0.76, nH 2 = -2.05; Frac = 0.21, F value = 62.75) and the BAPTA/ AM experiment (IC 50 = 5.2 μM and 512.9 μM; nH 1 = -2.02, nH 2 = -− 1.91; Frac = 0.15, F value = 5.61). The same scenario was observed in the BAPTA experiment (data not shown). Consequently, it is impossible to conclude on the effect of chelating Ca 2+ in response to SBV2-114 alone. Nevertheless, although the pre-incubation affected the Frac values and Hill coefficients, the preferred fit was still biphasic in both cases (F test), altogether indicating that the Ca 2+ response elicited by SBV2-114 is not responsible for the biphasic inhibition profile.
To address the involvement of Ca 2+ in the biphasic profile in another way, we tested if we could abolish the biphasic profile of SBV2-114 by inhibiting any potential release of Ca 2+ from the endoplasmic reticulum (ER), elicited by SBV2-114, by inhibiting phospholipase C (PLC), or by blocking any SBV2-114 mediated Ca 2+ influx through the P2X7 receptor [49]. For this, tsA201 cells were pre-incubated with either U-73122 or A 804595 for 30 min prior to and during the [ 3 H]GABA uptake. Inhibition of PLC with 5 μM U-73122 or P2X7 with 1 μM A 804595, however, did not convert the biphasic profile into a monophasic curve. Also, we still observed that the pre-incubation had an impact similar to pre-incubation with BAPTA/AM on the SBV2-114 inhibition curve at hBGT1 compared to Table 1 in terms of a lowered Frac value and Hill coefficients (Supplemental Fig. 1). Thus, it is not possible based on this approach either to deduce if the Ca 2+ increase is involved in the biphasic CRC of SBV2-114 at BGT1.
In order to avoid this pre-incubation, we tested if we could induce a biphasic inhibition profile of BPDBA and GABA by triggering Ca 2+ increases with ATP. Even though ATP modulated the BGT1 mediated [ 3 H]GABA uptake and decreased the total [ 3 H]GABA uptake by 9% (P < 0.001, One-way ANOVA, data not shown), the presence of ATP during the [ 3 H]GABA uptake had no effect on the type of inhibition profile of BPDBA and GABA at hBGT1 (Supplemental Fig. 2).

Anti-seizure Profile of SBV2-114
Several animal studies with the GAT1/BGT1 inhibitor EF1502 and the BGT1 inhibitor RPC-425 have reported anti-seizure effects of the compounds that cannot be ascribed solely to GAT1 [20][21][22][23]. This, together with the better selectivity of SBV2-114 (Supplemental Table 2) [46], impelled us to test the ability of SBV2-114 to prevent seizures in three mouse models. SBV2-114 was first tested in AGS Frings mice, which is a model of generalized seizures, to determine the ability of SBV2-114 to prevent sound-induced seizures and determine the approximate TPE for further testing. Frings mice were injected with 50 mg/kg SBV2-114 i.p. and tested after 15, 30, and 60 min. A subsequent time response study was conducted at a dose of 120 mg/kg and testing after 30, 60, and 120 min. The results from these studies show that SBV2-114 is brain-permeant and blocks sound-induced seizures in Frings mice in a dose-dependent manner at a TPE of 15 min (Table 3). Yet, a faster TPE cannot be excluded. In the subsequent MES test, SBV2-114 was administered at doses ranging from 50 to 200 mg/kg 15 min prior to the MES stimulus. SBV2-114 also prevented the seizures in the MES test and with an ED 50 of 139 mg/kg.  6 Inhibitory activity of SBV2-114 at hBGT1 transiently expressed in tsA201 cells and the effect of BAPTA/AM. Cells were pre-incubated 30 min with assay buffer ± 25 µM BAPTA/AM before SBV2-114 was tested for its ability to inhibit uptake of 30 nM [ 3 H] GABA for 3 min. The experiment shown is a representative example performed in triplicate and depicted as means ± S.D. According to the extra-sum-of-squares F test, the biphasic fit was preferred for both conditions P < 0.05 However, SBV2-114 produced rotarod impairment with a TD 50 of 163 mg/kg, and with 95% confidence intervals overlapping with the ED 50 . In the s.c.PTZ seizure test we only examined a single higher dose of 150 mg/kg SBV2-114. At this dose, SBV2-114 did not block s.c.PTZ-induced seizures or produce rotarod impairment ( Table 3). The reason for the mixed effect of SBV2-114 at rotarod impairment between the CF1 mice tested in the MES test and the s.c.PTZ seizure test is unclear. Further studies are needed to address if this could be due to e.g. sedative effects of SBV2-114.

Selectivity Studies of SBV2-114
Despite its dual activity at GAT1 and BGT1, EF1502 is completely devoid of any activity at a broad range of neurotransmitter systems including adrenergic, GABAergic and dopaminergic receptors in a radioligand binding study of 74 different neuroreceptors and ion channels [21]. Since SBV2-114 belongs to the same compound class as EF1502, and displayed anti-seizure effects, we wanted to address the overall selectivity of SBV2-114 in an effort to probe its potential as a tool compound to study BGT1 function in vivo. Hence, the selectivity of SBV2-114 for BGT1 over 42 neuroreceptors and transporters was tested in a similar manner as EF1502 through an affinity screening program supported by NIMH PDSP. From this evaluation, it is clear that SBV2-114 is fairly nonselective and competes with the radioligand binding to a number of other known drug targets including the serotonin receptors (5-HT 2B 5-HT 2C , 5-HT 5A , and 5-HT 6 ), the adrenergic receptor alpha 2A , and the dopamine transporter (DAT) (see Supplemental Table 4 for details) consequently hampering the usability of SBV2-114 to study BGT1 function in vivo. Additional studies are needed to investigate if interaction with any of these targets are implicated in the biphasic inhibition profile of SBV2-114 at hBGT1 in which case they would have to be endogenously expressed in our cell lines. Most likely, the involvement of any serotonin receptors and the adrenergic receptor alpha 2A can be disregarded, since they are not endogenously expressed in HEK-293 cells [50]. Recently, it was demonstrated that the highly related human serotonin transporter comprises an additional allosteric binding site in the extracellular vestibule above the orthosteric site [15]. This prompted us to investigate whether a similar allosteric site exists in BGT1 which could accommodate SBV2-114, and hence, could explain the observed biphasic behaviour. With the help of a homology model based on the SERT crystal structure along with docking studies, we identified hBGT1 Q299L in the orthosteric site and hBGT1 Tyr453 in the allosteric site to be possibly involved in SBV2-114 binding (Fig. 3). Furthermore, the top scored docking poses of SBV2-114 in the orthosteric and allosteric site suggested that it is sterically feasible for SBV2-114 to bind to both sites simultaneously (Fig. 3a). Nevertheless, we were not able to alter the biphasic inhibition profile of SBV2-114 by introducing a single mutation in the orthosteric site (hBGT1 Q299L) or the allosteric site (hBGT1 Y453A). However, by introducing a glutamine in hGAT3 at the position L314Q, which corresponds to Q299L in the orthosteric site of hBGT1, an increase in the activity of SBV2-114 and a strengthening of the biphasic profile was observed. Thus, residue Q299 in hBGT1 and residue L314 in hGAT3 appear to be important for the binding and activity of SBV2-114, but not solely responsible for the biphasic inhibition profile considering the unaltered inhibition profile of SBV2-114 at hBGT1 Q299L. The latter could be the result of possible compensation since the docking studies suggest that SBV2-114 can adopt two slightly different binding modes in the orthosteric pocket where the carboxyl group of SBV2-114 either interacts with the sodium ion Na1 or undergoes hydrogen bonding with Q299L ( Fig. 3b and  c). Lastly, the [ 3 H]GABA uptake data clearly shows that SBV2-114 inhibits GATs in an atypical manner, and hence, it is possible that other interacting residues not specific to BGT1 are accountable for the biphasic inhibition profile of SBV2-114.
To address the link between Ca 2+ , altered [ 3 H]GABA uptake [45], and the potential involvement in the biphasic inhibition profile of SBV2-114, we first demonstrated that SBV2-114 triggered an intracellular Ca 2+ increase at high concentrations independent of BGT1, via an unknown mechanism. Additionally, we showed that ATP, which triggered a large increase in intracellular Ca 2+ levels, lowered the BGT1-mediated [ 3 H]GABA uptake, supporting that a change in intracellular Ca 2+ levels modulate BGT1 activity as reported for GAT1 and GAT3 [45]. Nevertheless, the subsequent studies aimed at reducing the Ca 2+ involvement in the biphasic inhibition profile of SBV2-114 at hBGT1 were inconclusive. Thus, further studies are needed to investigate if the intracellular Ca 2+ increase triggered by the high concentrations of SBV2-114 is implicated in the observed biphasic behaviour, especially since the least potent part of the biphasic inhibition curve lies in the same range as the Ca 2+ releasing concentrations of SBV2-114.
Although a biphasic inhibition profile for the subfamily of GABA transporters has not been noted before, the phenomenon has been described within the family of SLC6A transporters [51][52][53][54]. A biphasic activity of a small molecule ATM7 has been reported at SERT, and was explained by an enhanced SERT uptake in response to nM concentrations of ATM7, while ATM7 inhibited SERT uptake at μM concentrations [53], with no notion of two binding sites for ATM7 at SERT. A mechanism like that for SBV2-114, however, is not in agreement with the [ 3 H]GABA uptake data presented here, since nM concentrations of SBV2-114 decreased the total [ 3 H]GABA uptake. Another example comes from a truncated version of the human DAT, which displayed a biphasic inhibition profile of [ 3 H]DA uptake presumably due to two DAT populations with different affinities for dopamine [51,55]. The fact that a biphasic inhibition profile of SBV2-114 was observed in cells transiently expressing BGT1 argues against the possibility that two BGT1 populations underlie the biphasic profile observed. Last, it seems unlikely that the presumed dimerization of BGT1 [56,57] accounts for the biphasic profile of SBV2-114 at BGT1, since other BGT1 inhibitors including NNC 05-2090, BPDBA as well as GABA inhibit BGT1 in a monophasic manner. If not, then a specific binding mode of SBV2-114 to the dimer should be in play.
BGT1 has attracted attention as a potential target for seizure management due to animal studies demonstrating antiseizure properties of the GAT1/BGT1 inhibitor EF1502 and the BGT1 inhibitor RPC-425 [20][21][22][23]. Still, the unaltered seizure susceptibility reported in BGT1 KO mice has questioned the importance of BGT1 as such [26], and BGT1 selective compounds are warranted to clarify the role of BGT1 in epilepsy. As reported here, SBV2-114 is indeed more selective for mBGT1 over mGAT1-3 (IC 50 > 1000 μM) as compared to both EF1502 and RPC-425 [46], which is an intriguing finding demonstrating how much a small chemical change from EF1502 to SBV2-114 can affect selectivity. The improved BGT1 preference of SBV2-114 encouraged us to test the anti-seizure properties of SBV2-114. We report that SBV2-114 dose-dependently prevented both audiogenic induced seizures in Frings mice and maximal electroshock induced seizures in CF-1 mice. In contrast, SBV2-114 did not prevent a minimal clonic seizure in the s.c.PTZ test at the dose (150 mg/kg, i.p) and pre-treatment time (30 min) tested. Though, this might be due to the apparent fast pharmacokinetics of SBV2-114 according to the testing in Frings mice illustrating that SBV2-114 was most effective at the earliest time point tested (15 min). Despite SBV2-114's preference for mBGT1 over mGAT1-3, results from the NIMH PDSP screen of neuro targets suggest that SBV2-114 also interacts with other molecular targets; a finding that renders it difficult to conclude that the anti-seizure effect of SBV2-114 is solely related to BGT1 inhibition. Nonetheless, the only known shared target for EF1502, RPC-425, and SBV2-114 is BGT. This, coupled with the fact that EF1502 has a good overall selectivity profile [21], collectively supports a role for BGT1 in epilepsy. This is in agreement with the synergistic effects of EF1502 and tiagabine in epilepsy models that cannot be linked to GAT1 exclusively [20][21][22][23]. Of note, the effect of combining RPC-425 and tiagabine did not have synergistic but only additive effects, which might be explained by differences in the pharmacokinetics between EF1502 and RPC-425 [20] rather than dismissing the importance of BGT1 in epilepsy.

Future Directions
While the molecular mechanism underlying the biphasic inhibition profile of SBV2-114 at BGT1 remains elusive, we show that it is a consistent phenomenon for this class of compounds all sharing a β-amino acid moiety attached to di-aromatic lipophilic side-chains. Although the in vivo suitability of SBV2-114 is poor due to its low selectivity and narrow therapeutic index, the data support previous findings suggesting an anti-seizure role for BGT1 inhibition. Further, SBV2-114 still constitutes a valuable in vitro tool compound to identify molecular mechanisms modulating GABA uptake, specifically the biphasic inhibition profile, and as a lead compound for the design of more selective BGT1 inhibitors. Ideally such studies would be aided by a BGT1 inhibitor radioligand for binding studies.