Next Article in Journal
The Role of Nuclear Factor of Activated T Cells 5 in Hyperosmotic Stress-Exposed Human Lens Epithelial Cells
Next Article in Special Issue
Mood Stabilizers in Psychiatric Disorders and Mechanisms Learnt from In Vitro Model Systems
Previous Article in Journal
Effect of Exogenous Auxin Treatment on Cell Wall Polymers of Strawberry Fruit
Previous Article in Special Issue
Therapeutic Effect and Mechanisms of Essential Oils in Mood Disorders: Interaction between the Nervous and Respiratory Systems
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 12
10.3390/ijms22126293
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters

by
Magdalena Kowalska
,
Łukasz Fijałkowski
and
Alicja Nowaczyk
*
Department of Organic Chemistry, Faculty of Pharmacy, Collegium Medicum in Bydgoszcz, Nicolaus Copernicus University, 2 dr. A. Jurasza St., 85-094 Bydgoszcz, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(12), 6293; https://doi.org/10.3390/ijms22126293
Submission received: 4 May 2021 / Revised: 31 May 2021 / Accepted: 10 June 2021 / Published: 11 June 2021
(This article belongs to the Special Issue Molecular Mechanisms of Mood Stabilizers)
Browse Figures
Versions Notes

Abstract

:
Thus far, many hypotheses have been proposed explaining the cause of depression. Among the most popular of these are: monoamine, neurogenesis, neurobiology, inflammation and stress hypotheses. Many studies have proven that neurogenesis in the brains of adult mammals occurs throughout life. The generation of new neurons persists throughout adulthood in the mammalian brain due to the proliferation and differentiation of adult neural stem cells. For this reason, the search for drugs acting in this mechanism seems to be a priority for modern pharmacotherapy. Paroxetine is one of the most commonly used antidepressants. However, the exact mechanism of its action is not fully understood. The fact that the therapeutic effect after the administration of paroxetine occurs after a few weeks, even if the levels of monoamine are rapidly increased (within a few minutes), allows us to assume a neurogenic mechanism of action. Due to the confirmed dependence of depression on serotonin, norepinephrine, dopamine and γ-aminobutyric acid levels, studies have been undertaken into paroxetine interactions with these primary neurotransmitters using in silico and in vitro methods. We confirmed that paroxetine interacts most strongly with monoamine transporters and shows some interaction with γ-aminobutyric acid transporters. However, studies of the potency inhibitors and binding affinity values indicate that the neurogenic mechanism of paroxetine’s action may be determined mainly by its interactions with serotonin transporters.

1. Introduction

Thus far, modern psychiatry classifies many types of depression including episodic, chronic, post-schizophrenic, endogenous, exogenous and other [1]. To explain the cause of depression a variety of hypotheses have been proposed. These include the monoamine hypothesis, neurogenesis hypothesis, neurobiology hypothesis, inflammation hypothesis, and the stress hypothesis [2,3,4]. The monoamine hypothesis suggests a relationship between a patient’s depressed mood and a decrease in the level of selected neurotransmitters. Until now, over 100 different neurotransmitters have been discovered [5]. It is widely accepted that in depressive states the key neurotransmitters are serotonin (SER), noradrenaline (norepinephrine) (NE), dopamine (DA) [6,7,8] and γ-aminobutyric acid (GABA) [9,10,11]. SER, NE and DA contain a catechol moiety, and thus belong to catecholamines, known as excitatory neurotransmitters. GABA, on the other hand, belongs to the group of acidic neurotransmitters and is the main inhibitory neurotransmitter in the brain. The role of controlling neurotransmitter levels has been confirmed by the discovery of antidepressants such as tricyclic antidepressants (TCA’s), selective serotonin reuptake inhibitors (SSRI’s) and serotonin norepinephrine reuptake inhibitors (SNRI’s) that increase serotonin levels in the brain [4]. Many of the modern generations of antidepressants act as inhibitors of one or several of the monoamine transporters (MATs) [12,13] and/or GABA transporters (GAT) [14]. The neurotransmitter transporter family dependent on sodium and chloride (SLC6 transporters) and including norepinephrine transporter (NET), dopamine transporter (DAT), serotonin transporter (SERT) and GABA transporters (GATs), encoded by SLC6A1-4 genes in humans are specifically known to be important for efficient neuronal synaptic transmission, hence providing neurotransmitter homeostasis in the CNS [15]. Briefly, inhibition of MATs/GATs leads to reduced clearance of neurotransmitters after synaptic release, thus increasing the neuronal signal intensity [11,16,17]. The inhibitors act by increasing the amount of neurotransmitters in the synapse, and in consequence alleviate or subside symptoms of depression [18].
Currently, one of the most commonly used antidepressants is paroxetine (PRX, (3S,4R)-3-([benzo[d][1,3]diokso-5-iloksy]metylo)-4-(4-fluorofenylo)piperidina), which belongs to the SSRI group [19]. It effectively increases the concentration of endogenous SER in the synaptic cleft [20,21,22]. It is noteworthy that it usually takes several weeks to achieve the desired antidepressant therapeutic effect of PRX. This therapeutic delay suggests that slow adaptive changes in the neural circuits during long-term pharmacological treatment possibly involve changes in gene expression and protein translation [23,24]. This in turn supports the hypothesis that depression is caused by disruption of functional and structural connections of the neural circuits that underlie the regulation of mood. However, the detailed molecular mechanisms of PRX have not been fully recognized and explained so far.
The studies presented in this publication are intended to answer the question about the possibility of PRX acting in the neurogenic mechanism. The primary purpose of this study is to analyze thoroughly the interactions between PRX and human MAT (hMAT: hSERT, hDAT, hNET), human GAT (hGAT: hGAT1-3, hBGT1). According to an extensive literature survey we have established that, however, the literature is replete with different studies concerning the four key transporter proteins, and there is a lack of pharmacological data on the molecular level. Therefore, we have undertaken the effort to fill this gap. In order to achieve the goal of this study we have conducted the in silico research involving dual directional analyses of the issue. The first mode of the analysis was designed to compare the analyzed ligand and protein structures by means of superimposition procedures. The results of this part of study delivered information about the structural similarities, which are often used in the pharmacological preliminary studies. The second part of analysis comprised a detailed inquiry of the intermolecular interactions present in the complexes especially hydrogen bonding. To compare the binding strength of PRX with individual transporters, the following reference compounds were selected, SER, DA, NE and GABA. The interactions were investigated by means of docking studies of selected antidepressant compounds performed on the crystal structure of hSERT [25] and on homology modeled hNET [26], hDAT [27] and hGAT1 [28]. It is noteworthy that the sequence identity between the template (i.e., sequences with known 3D structure) and the modeled sequence is 63.3% for hNET, 52% for hDAT [29], and 47% for hGAT1 and to the best of our knowledge this is the highest available identity. Reaching a general aim of the study requires also pharmacological investigation to provide the reader with experimental context. Even though the literature study shows that PRX–MAT interactions are extensively investigated by means of pharmacology, there are still gaps in the assessment of PRX–hGAT interactions. The pharmacological part of the present study includes the presentation of PRX–hGAT data that helps to understand and assess the neurogenesis hypothesis. This part of the study covers the investigation of the inhibitory effect of PRX, SER, NE and DA at the four hGATs (hGAT1-3, BGT1) using the [3H]GABA uptake assay.

2. Results

2.1. Molecular Docking Studies

2.1.1. Structural Similarity Studies for Neurotransmitters and Their Transporters

The Tanimoto similarity coefficient (Tsc) for pairwise comparison of molecules is probably the most widely used estimator of molecular similarity. This parameter by definition ranges from 0.0 (completely different) to 1.0 (equal molecules) [30]. Results of 2D similarity tests can be arranged in descending order of the Tsc: TscDA/NE = 0.59, TscSER/DA = 0.46 TscSER/NE = 0.35, TscSER/GABA = 0.25, TscDA/GABA = 0.12, TscGABA/NE = 0.05. They show a large structural diversity of the studied neurotransmitters, among which the greatest similarity of DA was noted to NE and SER. To obtain additional information concerning the shape of the investigated neurotransmitters and transporters, the active conformations were chosen for superimposition. The atoms common to these molecules were selected for the fitting procedure. Their similarity was calculated as RMS fit. The RMS routine provided estimates of how closely molecules fit to each other. The lower the RMS value, the better the similarity. The docked structures of the studied molecules were compared by their superimposition using a least-squares algorithm that minimizes the distances between the corresponding non-hydrogen atoms as shown in Figure 1.
The analysis of similarity of obtained active structures shows that all four compounds adopt very similar conformations in the target site of the protein RMS ≈ 0.1–1.5. The only significant deviation of the structure concerns the alkyl regions.
Recent studies on the identification of homologous structures for the studied proteins indicate that the binding site of hSERT shares a sequence identity of 57% with hDAT and 68% with hNET. Additionally, the identity between hNET and hDAT is 86% [31]. The hDAT shows the greatest homology to hNET within the binding site with an amino acid sequence identity of 78% [32,33]. The homology for each human transporter related to the others was sequence identified is approximately hDAT/hNET = 67%, hDAT/hSERT = 50%, hNET/hSERT = 53%. The evidence from these studies presents the lower sequence similarity for hNET/hDAT = 58% [34] compared with hGAT1/dDAT = 66% [35] of identity. All these values are quite similar to the sequence of proteins in all the transporters analyzed. By definition the percentage of structurally equivalent positions was defined as the percentage of the two alpha carbon atoms in the shorter of the sequences that are less than 3.5 Å of the equivalent atoms in the superposed structure [36]. All these values indicate a significant similarity in the sequence of proteins in all transporters analyzed (Figure 1). However, the analysis of structural similarities counted as RMS (in Å) superimposition the structures of two proteins [37], with alpha carbons of corresponding amino acids being taken together [38,39], results for MAT as following data for hSERT: hSERT/hDAT = 5.70 Å, hSERT/hNET = 5.66 Å. These results indicate a comparable similarity of both studied pairs, with a slightly higher similarity between hSERT and hNET than hSERT and hDAT. On the contrary, the obtained RMS value for hSERT/hGAT1 = 19.62 Å pair clearly indicates a significantly different 3D structure of the protein analyses. Similar RMS values were observed for hGAT1 with the MAT structures studied: hGAT1/hSERT = 19.62 Å, hGAT1/hNET = 19.82 Å, hGAT1/hDAT = 19.6 Å. These data indicate a low probability of strong interactions between PRX and hGAT1, while it is likely that a better interaction of PRX is seen with hNET and hDAT although hNET seems to have a slightly better interaction compared to hDAT. Furthermore, despite sharing high sequence identity and a similar MAT-structural fold (Figure 1), the pharmacology of neurotransmitter transporters is diverse, primarily because of the amino acid variations in the central binding site [40].

2.1.2. hSERT

Based on many pharmacological studies previously presented in the literature, PRX has been classified as an SSRI due to its very high inhibitory potencies for SERT. On the basis of predicted binding affinity, the highest value was observed with PRX (calculated pKihSERT ≈ 10.20, Table 1) compared to all studied compounds and hSERT.
In various studies, pIC50SERT was found to be between 8.96 and 10.40 [16,21,41,42,43]. This is in agreement with our docking results that shows that the most stable complex is for hSERT-PRX (EB ≈ 13.92 kcal/mol). The second strongest interaction with hSERT was SER (calculated pKiSERT ≈ 7.11, EB ≈ 9.7 kcal/mol, Table 1). The results of the present study demonstrate that DA, R/S-NE and GABA have low stability (EB ϵ (3.86; 3.99 kcal/mol). Other studied hSERT-neurotransmitters complexed have calculated pKi value ≤ 4 indicating weak ligand/inhibitor binding affinity. Comparing these values to other currently used antidepressants (pIC50SERT ≈ 9.10; 9.54; 8.05; 8.92; for fluoxetine (Prozac) [41], sertraline (Zoloft) [16], citalopram (Celexa) and venlafaxine (Effexor) [16], respectively) confirm the fact that PRX is the strongest SERT inhibitor currently known. Additionally, it is stronger than SER (pIC50SERT ≈ 7.00 [21,43]). Given the extraordinary strength with which PRX inhibits hSERT, the study of molecular factors that influence the neurogenesis phenomenon seems warranted. It is worth highlighting that in the biological activity under consideration, hSERT is the only transporter studied whose crystallographic structure is known (PDB ID: 5i6x, 3.14 Å [25]). The localization of PRX binding sites in crystals is well understood and described [19]. It is a central substrate binding pocket, commonly located approximately halfway across the membrane bilayer and marketed as S1 pocket (Table 2, Figure 2).
From chemical point of view hSERT is composed of three different regions; A, B, C. Subsite A is defined by Tyr95, Ala96, Asp98, Gly100 form transmembrane segments helix 1 (TMH1), and Phe335, Ser336, Gly338, Phe341, Val343 from transmembrane segments helix 6 (TMH6) and Ser438, Thr439, Ala441, Gly442 from transmembrane segments helix 8 (TMH8) and is define as a polar region surrounding Asp98. In crystallographic data it was found that regions A accommodates the polar, amine moiety of the PRX which form one hydrogen bond between NH2 group and Ala96 (Figure 2). Subsite B is the groove delineated between Ser438, Thr439, Ala441, Gly442 (TMH8) and Ala 169, Ile172, Ala173, Tyr176 from transmembrane segments helix 3 (TMH3). It has been confirmed by X-ray research that PRX locates a benzodioxol group in this groove. Subsite C is defined as a region between TMH6 and TMH10 (Val489, Lys490, Glu493, Glu494, Thr497, Gly498, Pro499, Leu502) and is located in the extracellular vestibule. Subsite C can interact with bulky drugs. Structurally B region is located opposite to subsite C and both are largely hydrophobic regions [46]. Commonly the pose obtained from X-ray diffraction is often denoted by many authors as ABC [22]. The results obtained in our PRX docking and validation experiments also remain to indicate PRX orientations in the ABC pose. In our experiment, the amine group of piperidine ring donates one hydrogen bond to the oxygen of the carboxylic moieties of Asp98 (Figure 2) with EHB ≈ −1.68 kcal/mol and bond length of 2.76 Å. The docking experiment confirms that the PRX molecule has a hydrophobic effect with an identical set of amino acids compared to crystallographic data (Figure 2, red circles and ellipses identify equivalent residua (in 3D superposition of structures)).
All calculated molecular configurations indicate that the energy of the single hydrogen bond between studied compounds and the set of neurotransmitters is relatively low (Table 1). The highest hydrogen bond energy is found for the major hydrogen bond component of Gly65.

2.1.3. hDAT

hDAT is the major target of addictive psychostimulants such as cocaine, which bind to the active site and prevent the conformational transition of the transporter, thereby inhibiting the reuptake of dopamine [40]. DAT mediates reuptake of DA (pIC50SERT << 4, pIC50NET ≈ 5.28, pIC50DAT ≈ 5.06 [21] − 5.61 [32]) from the synaptic cleft and thereby controls the termination of dopaminergic signaling.
The DAT structure distinguishes binding places defined as hDAT/cocaine (pocket 1) and hDAT/clomipramine (pocket 2) [6]. Cocaine is a high-potencies inhibitor of DAT (pIC50SERT ≈ 6.54, pIC50NET ≈ 5.48, pIC50DAT ≈ 7.14 [47]) and it is thought that its binding to DAT causes a rapid increase in extracellular dopamine levels [48] that produce the reinforcing effects leading to cocaine abuse. In addition, cocaine, despite the increase in DA in the synapse, does not act as an antidepressant. As confirmed in many studies, cocaine abuse and addiction are associated with an increased risk of depression [49]. Binding site for DA and cocaine in DAT overlap and are named by many authors as the central binding site, surrounded by TMs 1, 3, 6, 8 [33,50]. Clomipramine is a tertiary amine belonging to a dibenzazepine TCA [51]. Clomipramine has a stronger potency for the serotonin transporter (pIC50SERT ≈ 9.52, pIC50NET ≈ 7.42, pIC50DAT ≈ 5.66 [16]), compared to other TCAs [51]. Clomipramine binding site in hDAT protein is surrounded by TMHs 1, 3, 6, 10 and 11 [6]. Thus, two binding pockets have particular interest for the design of novel DAT interacting ligands. As shown in Figure 3 in both pockets PRX forms a hydrogen bond with Asp476. Comparing hydrophobic interactions in the analyzed cases reveals partial interaction compatibility (Figure 3, red circles identify equivalent residua (in 3D superposition structures)) but it also points to other residua in both cases (Figure 3, red arcs identify residua in hDAT/cocaine and hDAT/clomipramine, respectively (in superposition 3D structures)). The resulting energy values (Table 1) indicate that PRX forms a more energy-stable complex in the cocaine pocket. It is noteworthy that all neurotransmitters analyzed interact with the hDAT transporter (Table 1). It seems somewhat astonishing that SER interacted more strongly with hDAT (pIC50DAT ≈ 5.57 [21]) than the DA itself (pIC50DAT ≈ 5.06 [21]). Racemate R/S-NE, on the other hand, has a potency of pIC50DAT ≈ order of 4.28 [21]. However, this is confirmed by literature data.

2.1.4. hNET

Norepinephrine (also called noradrenaline) is a neuromodulator that in multiple ways regulates the activity of neuronal and non-neuronal cells [52]. It was one of the first neurotransmitters to be discovered [53]. It was isolated in 1901 and is the first hormone obtained in the crystalline state [54]. The presence of norepinephrine in two main enantiomers (R/S-NE) was confirmed by structural (crystallographic and spectroscopic studies [54,55]). In crystallographic study it was proved that NE has a zwitterionic structure, formed by proton transfer from the meta phenolic hydroxyl group to the nitrogen atom (−O(OH)Ph-CH(OH)-CH2-NH3+) [56]. It was founded approx. 6.5–7.7 Å due to separation of charged centers [55], which determines the molecule of the interaction with neurons [28]. These conformers are stabilized by intramolecular hydrogen bonds between the nitrogen atom of the amino group and hydrogen of hydroxyl group attached to an ethyl amine carbon atom [55]. NE is the second major biogenic amine that has been proposed to be causally involved in the pathophysiology of major depression and in the mechanism of antidepressant drug action. The norepinephrine transporter (NET) is a transmembrane protein responsible for transporting norepinephrine into the synaptic terminals of the central and peripheral nervous systems as well as neuroendocrine adrenal chromaffin cells [57]. The selective NE reuptake potencies inhibitor, reboxetine (pIC50 ≈ 8.96 [16]), has demonstrated equivalent efficacy to the TCA in some studies and is approved as an antidepressant in Europe but not in the USA [58] (Figure 4).
The inhibitory potencies R/S-NE for its transporter lie in the range pIC50NET ϵ (6.17 [32]; 5.06 [21]). Comparison of the data obtained for our docking experiment (Table 1) clearly leads to the conclusion that S-NE shows stronger binding affinity than the other studied neurotransmitters. It was found in our experiment that pKiNET ≈ 5.59 for S-NE and pKiNET ≈ 4.12 for R-NE. Another strength of the interacting neurotransmitter is DA (calculated pKiNET ≈ 5.31), which is confirmed by the literature data (pIC50NET ≈ 5.28 [21]). Among the studied catecholamines, serotonin is characterized by the weakest affinity for NE transporters (calculated pKihNET ≈ 4.33 and pIC50NET ≈ 4.85 [21]). Calculated binding affinity pKihNET ≈ 3.31 for GABA suggests interaction with hNET. These data confirm that NET has some selectivity relative to NE. When it comes to the effects of PRX on NET, it is much stronger (calculated pKihNET ≈ 7.61 and pIC50NET ϵ 6.46 [21] − 7.40 [16]) compared to all the neurotransmitters in the studied set.

2.1.5. hGAT1

In regions such as the cerebral cortex, hippocampus, thalamus, basal ganglia, cerebellum, hypothalamus, and brainstem, GABA represents about one-third of the synapses. GAT1 is the most copiously expressed GAT in the CNS and is mainly localized into the presynaptic axon terminal and into a few astrocytic structures [59]. The GABA chain is made of three single carbon-carbon bonds, thus its structure is highly flexible. This determines the ability to adapt various conformations that determine binding, as ligand, in the proteins with which it interacts. Crystallographic studies have shown that both crystal and aqueous solution of GABA occur in the form of a zwitterion, which in the body makes it a carrier of electric current and thus creates the possibility of carrying stimuli. It is widely accepted that distance between cationic and anionic center approx. A value of 5 Å determines the neuronal and glial blocking effect in selective inhibitors. Computational methods have shown that in an aqueous solution more than 94% of GABA molecules adapt a folded conformation with separation of charged centers below 5 Å [28]. GABA is transported by GAT in a somewhat folded conformation (Figure 5). This means GABA interacts with neuronal and astrocytic structures.
In the preliminary analysis of the docking it was observed that GABA interacts with the active site of hGAT1 with calculated pKihGAT1 ≈ 4.49 and [3H]GABA pIC50hGAT1 ≈ 5.00 [60,61]. The studied catecholamines show significantly less interaction with hGAT1 as compared to GABA (Table 1). It is noteworthy that PRX shows the strongest interactions with hGAT1 (calculated pKihGAT1 ≈ 4.35) compared to SER, NE and DA. (Table 1). PRX creates a bifurcated hydrogen bond, similar to the GABA–hGAT1 complex, where NH2 from PRX is involved in a COO...NH2...CONH incorporating the COO from Asp451 and the NHCO from Asp451 (Figure 5). From the data in Table 1 it can be seen that this bifurcated system is almost geometrically (LHB ≈ 2.92 Å) and energetically (EHB ≈ −3.41 kcal/mol) symmetric, where the major hydrogen bond component is fixed into the side-NH-CO of Asp451 and the minor involves the side-COO of Asp451 residue (Figure 5).

2.1.6. Validation Experiment

The re-docking method was applied to confirm legitimacy of the docking procedure [62] and was performed using Vina software [63]. Results of redocking in Vina (optimal poses) were saved as pdbqt files. Subsequently, the AutoDock Tools (ADT) program was employed to calculate the superimposition of the ligand and protein complexes. The validation process indicates a high similarity of PRX corresponding poses obtained using AutoDock and Vina (Figure 6).

2.2. Inhibitory Activities of PRX, SER, NE and DA at hGATs

PRX concentration-dependently inhibits the [3H]GABA uptake mediated by all four hGATs with mid-high micromolar potency in the range of 85.6–256.1 µM with the rank order: hBGT1 < hGAT3 ≈ hGAT2 < hGAT1 (Figure 7a and Table 3).
SER, NE and DA display limited inhibitory activities at hGAT1-3 and hBGT1 with less than 10% inhibition observed at compound concentrations up to 100 μM (Figure 7b).
Structurally PRX can be thought of as a substituted derivative of piperidine where in position 3 there is a group of benzodioxol and in position 4 the fluorophenyl group. The closest to PRX active structure in hGAT is nipecotic acid, which belongs to the leading structures for the development of structural hGAT inhibitors and has the following inhibition potential: pIC50hGAT1 ≈ 4.7, pIC50hGAT2 ≈ 3.3, pIC50hGAT3 ≈ 4.0, pIC50hBGT1 ≈ 2.5 [61]. A comparison of the inhibition potential of PRX (Table 3) and nipecotic acid indicates that the replacement of the carboxyl group with benzodioxol and the addition of the fluorophenyl group change the activity profile from hGAT1 and hGAT3 to hGAT2, and hBGT1. Consequently, one can draw a preliminary conclusion that such a modification of the structure makes the molecule more susceptible to inhibition of hGAT located outside the CNS. Which in turn indicates that PRX has less influence on the neurogenesis process in the brain.

3. Discussion

Decreased amounts of SER, NE, DA and GABA have been observed in people with depression. Increasing SER in key CNS pathways and at desired serotonin receptor subtypes hypothetically mediates therapeutic actions in depression, and other diseases with similar symptoms such as OCD, PD and bulimia [64]. This is the result of the powerful influence of SER on emotionality that has been proven many times over [65]. SER and GABA are mainly, but not only, the basic inhibitory neurotransmitters in the CNS. They play an important role in processing neuronal information as well as regulating neurogenesis (proliferation, differentiation and migration of neural stem cells (NSC)) [66,67]. Numerous studies have shown that deficits in serotonin, norepinephrine and GABAergic neurotransmission and reduced neurogenesis are associated with the ethology of pathological anxiety and various mood disorders including depression [67,68,69]. Evidence of reduced neurogenesis has been observed in animal models of depression and in postmortem studies of individuals who had been diagnosed with major depressive disorder. However, the mechanisms driving these alterations in neurogenesis are not fully understood [70]. SER, NE and GABA have been shown to regulate hippocampal neurogenesis and neuronal development in both children and adults [71,72]. In the adult brain, as in the embryonic nervous system, SER, NE and GABA depolarizes neural progenitors and immature neurons [73]. GABA also regulates quiescence of NSCs, differentiation of NSCs into neural progenitor cells (NPCs), maturation of NPCs into granule cells and synaptic integration of adult-born granule cells into the existing circuitry of the hippocampus [72]. A common feature of depression and anxiety disorders are SER, NE and GABA deficits resulting in emotional control disorders [65]. Recently, a direct link between chronic antidepressant treatment and an enhancement of SER, NE and GABA transmission were found [65,67]. This is in fact that SER, NE and GABA concentration are reduced in cortical brain and CSF in major depression, but its deficit can be reversed by chronic SSRI and electroconvulsive therapy [74,75].
PRX is an antidepressant drug known by the commercial names Aropax, Paxil, Pexeva, Seroxat, Sereupin and Brisdelle. It exhibits the highest known inhibitory potency for the active site of the hSERT (Table 1). A plethora of publications devoted to PRX (e.g., PubMed database 1992–2021 the word “paroxetine” indicates in about 6200 publications, Science Direct about 22000, Google Scholar about 54,500 works) proves the great importance of this drug in current medical therapies. These works provide a lot of information on the different aspects of taking this medicine and its mechanisms of action and it is clear that PRX is currently clinically approved for the treatment of numerous neurological disorders and not only depression [22].
After administration of PRX the level of SER rapidly increases (after a few minutes), while the antidepressant effect does not appear before after several weeks of chronic treatment [18,53,67,76]. This may indicate a possible different mechanism of action of PRX, and the most probable seems to be the neurogenic mechanism of action [77]. For this reason, we undertook studies of PRX interactions with four basic transporters, such as three MAT and GAT, the inhibition of which is crucial for increasing the level of basic neurotransmitters in the brain. At the beginning using the results of research based on molecular modeling techniques, the structural similarity of studied neurotransmitters: SER, NE, DA and GABA, as well as their respective targets, i.e., MAT and GAT1 transporters were analyzed. Using various techniques, such as the Tanimoto coefficient [30] and superimposition tests [63], it was shown that some structural similarities were found within the tested neurotransmitters as well as their transporters. Nevertheless, the analysis of structural similarities calculated as RMS determined by superimposing the structures of two proteins seems to be of key importance in understanding the interaction of various compounds with the analyzed objectives [37]. These results, in contrast to the data routinely presented in the literature based on the amino acid sequence identity technique [29,32,33,34,35], better highlight the structural differences of the analyzed transporters. These differences are of particular importance within the ligand binding site. The information obtained in this study indicates some unsimilarities between MAT and GAT1 (Figure 1), which may explain the low inhibitors of potency and binding affinity values of PRX for GAT1-3 and BGT1. On the other hand, a more detailed analysis of the literature data on the interactions of PRX and nipecotinic acid (as one of the leaders for the hGAT inhibitor family) provides information about the direction of interaction of these molecules as substituted piperidine derivatives. The data we obtained based on the in silico and in vitro methods indicate that PRX interacts with each monoamine transporter and shows some interactions with the hGATs (Table 1 and Table 3, Figure 7). As shown in the Section 2, the obtained in silico data are confirmed by the relevant pharmacological studies for MAT previously published in the literature. It has recently been found that PRX may promote the proliferation of nerve cells in vitro, as well as hippocampal neurogenesis in human and animal models [78]. Nowadays it is well accepted that depletion of SER in the brain results in suppression of neurogenesis in the adult hippocampal. At the same time, raising levels of monoamine, serotonin and norepinephrine and GABA increases the rate of neurogenesis. Since PRX is the strongest known hSERT inhibitor, the study of the possible mechanism of antidepressant action clearly points to issues related to the phenomenon of neurogenesis. However, the neurogenesis process itself is not yet sufficiently explained at the molecular level. Since this phenomenon is known to be linked to an increase in serotonin and norepinephrine and GABA levels within neurons, it seems that the potential for blocking suitable transporters is at the heart of this phenomenon.
It is important to emphasize that the interaction of PRX with the hGATs seem to be a more complex issue due to the multitude of different subtypes of hGAT. Hitherto four subtypes of plasma membrane human transporters for GABA (hGAT1–3 and hBGT1) have been identified [79,80,81], of which hGAT1 and hGAT3 have the highest expression levels in the mammalian central nervous system (CNS). hBGT1 is only found in scarce amounts and hGAT2 is not found in the brain parenchyma at all [82]. hGAT1 is predominantly located on GABAergic nerve terminals [35,59,83], while hGAT3 and hBGT1 are commonly associated with perisynaptic and distal astrocytic sites [28,60]. The location and level of individual transporters in CNS seem to be crucial in the neurogenesis process. Our experiments show that PRX interacts with all hGAT (Table 1 and 3 and Figure 7) to varying degrees. Interaction of PRX with individual hGATs can be arranged in descending order: pIC50hBGT1 ≈ 4.10, pIC50hGAT3 ≈ 4.00, pIC50hGAT2 ≈ 3.90, pIC50hGAT1 ≈ 3.6 (Table 3) and pKihGAT1 ≈ 4.35 (Table 1). In vitro data clearly confirm the strongest interaction with hBGT1. However, due to the very low level of BGT1 expression in the brain (about 100–1000 times lower than GAT1 [84]), BGT1 can at most be responsible for 0.1–1.0% of the GABA transport. In addition, neurotransmitters diffuse rapidly out of the synaptic cleft on a low microsecond time scale until they bind to transporters and are removed [85,86,87], hence BGT1′s functional role in the neurogenesis process seems negligible. For these reasons and in light of the low inhibitory potency of PRX at hGAT1 and hGAT3, it appears that the interaction of PRX with the hGATs has little effect on the PRX induced neurogenesis process. Due to this analysis one can assume that if PRX induced neurogenesis it potentially connected with its inhibition to hSERT due to higher value inhibitory potency and binding affinity.

4. Materials and Methods

4.1. Pharmacological Studies

Flp-In Chinese hamster ovary (CHO) cell lines stably expressing hGATs (hGAT1-3 or hBGT1) used for the pharmacological studies have been described previously and were cultivated accordingly [88]. The [3H]GABA competition uptake assay was performed exactly as previously described [61], and the [3H]GABA uptake data were normalized to the percentage of total uptake in the individual experiments. Data presented are the pooled data of three independent experiments with three technical replicates. Concentration–response curves were fitted with GraphPad Prism (version 9.0.0, Yosemite, GraphPadSoftware, San Diego, CA, USA) as outlined previously [88].

4.2. In Silico Studies

4.2.1. Docking Study

The calculation procedures applied in the study are typical for processing of docking studies [63,89,90,91].

Ligand Preparation

For the 3D molecular structure calculations, the Gaussian 09 (version D.01. for Unix/Linux) package was used [92]. The initial accepTable 3D structures of 6 studied compounds were downloaded (as mol2 file) from ZINC [93]. Later, the GaussView [92,94] was applied for preparation of Gaussian input files. All the molecules were geometry-optimized in water described by the PCM (polarizable continuum model). DFT/B3LYP level of theory, 6-311 + G(d, p) basis set. After geometry optimization (the root-mean-square gradient value smaller than 10−6 a.u.) compounds were saved as mol2 files using the GaussView. Subsequently, torsionals and number of active torsions for ligands were defined and the Gasteiger charges were assigned to each compound via AutoDockTools (ADT) [95]. Finally, ligands prepared for docking were saved as pdbqt files.

Monoamine and γ-Aminobutyric Acid Transporters Preparation

Despite numerous protein sequences in RCSB Protein Data Bank, lack of three-dimensional structures of many important drug targets is still observed. As a consequence, homology modelling is the different option to construct an accepTable 3D model of the protein. Therefore, in this study docking experiments were mostly performed on the homology modeled proteins (hDAT, hNET and hGAT1), except hSERT, where the X-ray structure is available.
  • hSERT
The structure of the ts3 hSERT was gained from the RCSB Protein Data Bank (5I6X pdb access code) [25] First, the pdb structure was adapted for ADT software environment. Accordingly, one crystallographic water molecule and following ligands: dodecyl-beta-d-maltoside (LMT), cholesterol (CLR), Paroxetine (PRX), 2-acetamido-2-deoxy-beta-D-glucopyranose (NAG) were removed. Subsequently, missing hydrogen atoms were added and protein was saved in pdbqt format [96].
  • hDAT, hNET hGAT1
As mentioned before, for hDAT, hNET and hGAT1 homology human protein sequences were obtained from the Swiss-Prot database (accession numbers: Q01959, P23975 and P30531, respectively). In this study the most actual homology modeled channels were used [6,26,44,97]. After that, pdb files were opened in ADT, which read coordinates. Subsequently, charges were added, correct atom types were assigned and nonpolar hydrogens were merged. Finally, according to the procedure the prepared protein was saved as a pdbqt file.

Molecular Docking

Preparation of the inputs for the modeling was carried out in full compliance with the applicable procedures. A grid box with a dimension of 60 × 60 × 60 Å3 and a grid spacing of 0.375 Å was built in the middle of the pore forming sequences. In this in silico experiment we consider the interaction between ligands and intracellular pore gate formed from the transmembrane helices. For each model, drug-binding pockets in the cavity forming part of the substrate permeation pathway were identified (see Table 2). The rigid docking procedures were performed using the Lamarckian genetic algorithm of Autodock 4.2 software. The optimized docking parameters were set as default values, except the number of runs which was 100. Torsionals in the residuals of binding pockets were not rotatable. A cluster analysis was performed on the docked results using an RMS tolerance of 2 Å. In each case the best docking result was considered as complex with the lowest binding energy. Hydrogen bindings between docked potent agents and related macromolecules were analyzed using ADT, the AutoDockTools program (ADT. Version 1.5.4) [95].

4.3. Similarity and Superimposition Study

The Tanimoto similarity coefficient (Tsc) for pairwise was computed as available on the ChemMine tools server (http://chemminetools.ucr.edu/ (accessed on 1 June 2021)) [30]. To obtain additional information concerning the shape of the investigated drug, the active conformations were chosen for superimposition. The atoms (except hydrogen atoms) common to these molecules were selected for the fitting procedure using SCIGRESS software, version 3.4.4, www.scigress.com. Their similarity was calculated as RMS fit. The RMS routine provided estimates of how closely molecules fit with each other. The lower the RMS value, the better the similarity [37].

5. Conclusions and Perspectives

The mechanism of adult neurogenesis has not been fully understood and described so far. However, in light of current research, it has been found that it is initiated by an increase in the concentration of neurotransmitters in the brain. For this reason, it seems logical to suppose that drugs used in the pharmacotherapy of diseases with etiology resulting from a decreased level of neurotransmission with the participation of basic mediators (SER, NE, DA, GABA) in the nervous system may potentially act in the neurogenic mechanism. Undoubtedly, depression is one of these diseases, and PRX, as the most widely used antidepressant, meets the requirements of the research assumption. It should be emphasized that PRX belongs to the group of the most studied molecules. For this reason, we have a lot of well-documented data that can be used in the context of analyzing new mechanisms of action.
The study presented in this publication predominantly was meant to answer the question about the possibility of PRX acting in the neurogenic mechanism. Data from molecular studies indicate a potential neurogenic effect of PRX, which may be mainly due to the strong inhibition of hSERT. In order to identify the causes of different values of the inhibition parameters (pKi and pIC50), it was reasonable to compare the structural similarity of the analyzed structures of the analyzed objects. First, the structural similarity of four basic neurotransmitters: SER, NE, DA and GABA and the corresponding targets, i.e., MAT and GAT1 transporters were analyzed, using the results of research based on molecular modeling techniques.
Using various techniques such as the Tanimoto test and overlap tests, it was shown that some structural similarities were found in the tested carriers and their transporters. Nevertheless, the analysis of structural similarities calculated as RMS determined by the superposition of the structures of two proteins considering the alpha carbon atoms of the respective amino acids seems to be crucial for understanding the interactions of various compounds with the analyzed targets. These results, in contrast to the data routinely presented in the literature, based on the amino acid sequence identification technique, highlight the structural differences of the analyzed transporters. These differences are of particular importance within the ligand binding site. The information obtained from this study indicates significant differences between MAT and GAT1, which may explain the low GAT1-3 and BGT1 inhibition values (potency and affinity, respectively) by PRX.
An additional, but very important, prospective achievement of the research carried out is the confirmation of some activity of PRX towards hGAT. This enables the development of research based on the modification of the PRX structure in order to increase the affinity for particular types of hGAT.

Author Contributions

M.K., Ł.F., A.N. conceived and directed the project, designed the study, collected the data, analyzed the data, interpreted the results and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the research grant from the EIRU program at the NCU entitled EI DEBUTS “Finding and understanding mechanisms of action of selected compounds active on central nervous system monoamine transporters”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Petrine Wellendorph and Maria E.K. Lie from the Department of Drug Design and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Universitetsparken 2, 2100 Copenhagen, Denmark for conducting the in vitro part of this research. The authors would like to thank Wojciech Płonka from the Universität Hamburg, Center for Bioinformatics (ZBH), Hamburg, 20146, Germany and FQS Poland (Fujitsu Group), Parkowa 11, 30-538 Cracow, Poland for his technical support during superimposition part of this research.

Conflicts of Interest

The authors declare no conflict of interest.

Compliance with Ethical Standards

Not applicable.

Abbreviations

θhydrogen bond angle
#hydrogen bond components: from the protein
%hydrogen bond components: from the ligand
acchydrogen bond acceptor
ADTAutoDockTools
CHOChinese hamster ovary
CNScentral nervous system
CSFcerebrospinal fluid
DAdopamine
DATproteins of dopamine transporter
EBcomplex energy binding
EHBhydrogen bond energy
GABAγ-aminobutyric acid
GAT1GABA transporter isoform 1
GAT2GABA transporter isoform 2
GAT3GABA transporter isoform 3
GDAgeneralized anxiety disorder
hSERThuman proteins of serotonin transporter
hDAThuman proteins of dopamine transporter
hNEThuman proteins of norepinephrine transporter
hGAT1human GABA transporter isoform 1
hGAT2human GABA transporter isoform 2
hGAT3human GABA transporter isoform 3
hBGT1human betaine/GABA transporter
LHBhydrogen bond length
MATmonoamine transporters
MDDmajor depressive disorder
NEnorepinephrine
NETproteins of norepinephrine transporter
NPCneural progenitor cells
NSCneural stem cells
OCDobsessive-compulsive disorder
PDpanic disorder
PDB IDProtein Data Bank Identifier
pIC50xxxinhibitory potency of xxx transporters
pKixxxbinding affinity of xxx transporters
PMDDpremenstrual dysphoric disorder
PRXparoxetine
PTSDpost traumatic stress disorder
R/S-NER/S enantiomers norepinephrine
RMSroot mean square
SADsocial anxiety disorder
SERserotonin
SERTproteins of serotonin transporter
SNRIserotonin norepinephrine reuptake inhibitors
SSRIselective serotonin reuptake inhibitors
TCAtricyclic antidepressants
TscTanimoto similarity coefficient
TMtransmembrane segments
TMHtransmembrane segments helix

References

  1. Rantala, M.J.; Luoto, S.; Krams, I.; Karlsson, H. Depression subtyping based on evolutionary psychiatry: Proximate mechanisms and ultimate functions. Brain Behav. Immun. 2018, 69, 603–617. [Google Scholar] [CrossRef] [PubMed]
  2. Butterweck, V. Mechanism of action of St John’s wort in depression. CNS Drugs 2003, 17, 539–562. [Google Scholar] [CrossRef] [PubMed]
  3. Walker, F.R. A critical review of the mechanism of action for the selective serotonin reuptake inhibitors: Do these drugs possess anti-inflammatory properties and how relevant is this in the treatment of depression? Neuropharmacology 2013, 67, 304–317. [Google Scholar] [CrossRef] [PubMed]
  4. Dean, J.; Keshavan, M. The neurobiology of depression: An integrated view. Asian J. Psychiatry 2017, 27, 101–111. [Google Scholar] [CrossRef]
  5. Moon, J.-M.; Thapliyal, N.; Hussain, K.K.; Goyal, R.N.; Shim, Y.-B. Conducting polymer-based electrochemical biosensors for neurotransmitters: A review. Biosens. Bioelectron. 2018, 102, 540–552. [Google Scholar] [CrossRef]
  6. Ravna, A.W.; Sylte, I.; Dahl, S.G. Structure and localisation of drug binding sites on neurotransmitter transporters. J. Mol. Mod. 2009, 15, 1155–1164. [Google Scholar] [CrossRef]
  7. Zangrossi, H., Jr.; Del Ben, C.M.; Graeff, F.G.; Guimarães, F.S. Serotonin in panic and anxiety disorders. In Handbook of Behavioral Neuroscience; Elsevier: Amsterdam, The Netherlands, 2020; Volume 31, pp. 611–633. [Google Scholar]
  8. Whitaker-Azmitia, P.M. The discovery of serotonin and its role in neuroscience. Neuropsychopharmacology 1999, 21, 2–8. [Google Scholar] [CrossRef]
  9. Fijałkowski, Ł.; Sałat, K.; Podkowa, A.; Zaręba, P.; Nowaczyk, A. Potential role of selected antiepileptics used in neuropathic pain as human GABA transporter isoform 1 (GAT1) inhibitors—Molecular docking and pharmacodynamic studies. Eur. J. Pharm. Sci. 2017, 96, 362–372. [Google Scholar] [CrossRef]
  10. Monteleone, P.; Maj, M.; Iovino, M.; Steardo, L. GABA, depression and the mechanism of action of antidepressant drugs: A neuroendocrine approach. J. Affect. Disord. 1990, 20, 1–5. [Google Scholar] [CrossRef]
  11. Schwartz, T.L.; Nihalani, N. Tiagabine in anxiety disorders. Expert Opin. Pharmacother. 2006, 7, 1977–1987. [Google Scholar] [CrossRef]
  12. Iversen, L. Neurotransmitter transporters: Fruitful targets for CNS drug discovery. Mol. Psychiatry 2000, 5, 357–362. [Google Scholar] [CrossRef]
  13. Racagni, G.; Brunello, N. Physiology to functionality: The brain and neurotransmitter activity. Int. Clin. Psychopharmacol. 1999, 14, S3–S7. [Google Scholar] [CrossRef]
  14. Kalueff, A.V.; Nutt, D.J. Role of GABA in anxiety and depression. Depress. Anxiety 2007, 24, 495–517. [Google Scholar] [CrossRef]
  15. Kristensen, A.S.; Andersen, J.; Jørgensen, T.N.; Sørensen, L.; Eriksen, J.; Loland, C.J.; Strømgaard, K.; Gether, U. SLC6 neurotransmitter transporters: Structure, function, and regulation. Pharmacol. Rev. 2011, 63, 585–640. [Google Scholar] [CrossRef]
  16. Iversen, L. Neurotransmitter transporters and their impact on the development of psychopharmacology. Br. J. Pharmacol. 2006, 147, S82–S88. [Google Scholar] [CrossRef] [Green Version]
  17. Tuma, A.H.; Maser, J.D. Anxiety and the Anxiety Disorders; Routledge: New York, NY, USA, 2019. [Google Scholar]
  18. Boas, G.R.V.; de Lacerda, R.B.; Paes, M.M.; Gubert, P.; da Cruz Almeida, W.L.; Rescia, V.C.; de Carvalho, P.M.G.; de Carvalho, A.A.V.; Oesterreich, S.A. Molecular aspects of depression: A review from neurobiology to treatment. Eur. J. Pharmacol. 2019, 851, 99–121. [Google Scholar] [CrossRef]
  19. Coleman, J.A.; Navratna, V.; Antermite, D.; Yang, D.; Bull, J.A.; Gouaux, E. Chemical and structural investigation of the paroxetine-human serotonin transporter complex. eLife 2020, 9, e56427. [Google Scholar] [CrossRef]
  20. Davis, B.A.; Nagarajan, A.; Forrest, L.R.; Singh, S.K. Mechanism of paroxetine (paxil) inhibition of the serotonin transporter. Sci. Rep. 2016, 6, 23789. [Google Scholar] [CrossRef] [Green Version]
  21. Bourin, M.; Chue, P.; Guillon, Y. Paroxetine: A review. CNS Drug Rev. 2001, 7, 25–47. [Google Scholar] [CrossRef]
  22. Kowalska, M.; Nowaczyk, J.; Fijałkowski, Ł.; Nowaczyk, A. Paroxetine—Overview of the molecular mechanisms of action. Int. J. Mol. Sci. 2021, 22, 1662. [Google Scholar] [CrossRef]
  23. Oh, S.-J.; Cheng, J.; Jang, J.-H.; Arace, J.; Jeong, M.; Shin, C.-H.; Park, J.; Jin, J.; Greengard, P.; Oh, Y.-S. Hippocampal mossy cell involvement in behavioral and neurogenic responses to chronic antidepressant treatment. Mol. Psychiatry 2020, 25, 1215–1228. [Google Scholar] [CrossRef]
  24. Duman, R.S.; Aghajanian, G.K.; Sanacora, G.; Krystal, J.H. Synaptic plasticity and depression: New insights from stress and rapid-acting antidepressants. Nat. Med. 2016, 22, 238–249. [Google Scholar] [CrossRef] [Green Version]
  25. Coleman, J.A.; Green, E.M.; Gouaux, E. X-ray structures and mechanism of the human serotonin transporter. Nature 2016, 532, 334–339. [Google Scholar] [CrossRef] [Green Version]
  26. Selvaraj Manikandan, A.; Ramanathan, T.; Selvarajan, K.; Kesavanarayanan, T.L.; Salleh, M. Identification of Benzoxazolinone Derivatives Based Inhibitors for Depression and Pain Related Disorders Using Human Serotonin and Norepinephrine Transporter as Dual Therapeutic Target: A Computational Approach. Int. J. Drug Deliv. 2014, 6, 389–395. [Google Scholar]
  27. Penmatsa, A.; Wang, K.H.; Gouaux, E. X-ray structures of Drosophila dopamine transporter in complex with nisoxetine and reboxetine. Nat. Struct. Mol. Biol. 2015, 22, 506–508. [Google Scholar] [CrossRef] [Green Version]
  28. Nowaczyk, A.; Fijałkowski, Ł.; Kowalska, M.; Podkowa, A.K.S. Studies on the activity of selected highly lipophilic compounds toward hGAT1 inhibition: Part II. ACS Chem. Neurosci. 2019, 10, 337–347. [Google Scholar] [CrossRef]
  29. Haddad, Y.; Heger, Z.; Adam, V. Guidelines for homology modeling of dopamine, norepinephrine, and serotonin transporters. ACS Chem. Neurosci. 2016, 7, 1607–1613. [Google Scholar] [CrossRef]
  30. Bero, S.; Muda, A.; Choo, Y.; Muda, N.; Pratama, S. Similarity measure for molecular structure: A brief review. J. Phys. 2017, 2017, 012015. [Google Scholar] [CrossRef] [Green Version]
  31. Koldsø, H.; Christiansen, A.B.; Sinning, S.; Schiøtt, B. Comparative modeling of the human monoamine transporters: Similarities in substrate binding. ACS Chem. Neurosci. 2013, 4, 295–309. [Google Scholar] [CrossRef] [Green Version]
  32. Giros, B.; Wang, Y.-M.; Suter, S.; McLeskey, S.B.; Pifl, C.; Caron, M.G. Delineation of discrete domains for substrate, cocaine, and tricyclic antidepressant interactions using chimeric dopamine-norepinephrine transporters. J. Biol. Chem. 1994, 269, 15985–15988. [Google Scholar] [CrossRef]
  33. Wang, K.H.; Penmatsa, A.; Gouaux, E. Neurotransmitter and psychostimulant recognition by the dopamine transporter. Nature 2015, 521, 322–327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ortore, G.; Orlandini, E.; Betti, L.; Giannaccini, G.; Mazzoni, M.R.; Camodeca, C.; Nencetti, S. Focus on Human Monoamine Transporter Selectivity. New Human DAT and NET Models, Experimental Validation, and SERT Affinity Exploration. ACS Chem. Neurosci. 2020, 11, 3214–3232. [Google Scholar] [CrossRef] [PubMed]
  35. Zafar, S.; Nguyen, M.E.; Muthyala, R.; Jabeen, I.; Sham, Y.Y. Modeling and Simulation of hGAT1: A Mechanistic Investigation of the GABA Transport Process. Comput. Struct. Biotechnol. J. 2019, 17, 61–69. [Google Scholar] [CrossRef] [PubMed]
  36. Madhusudhan, M.; Marti-Renom, M.A.; Sanchez, R.; Sali, A. Variable gap penalty for protein sequence–structure alignment. Protein Eng. Des. Sel. 2006, 19, 129–133. [Google Scholar] [CrossRef] [Green Version]
  37. Wong, C.; Currie, J. Teaching with SCIGRES Molecular Modeling in Chemistry; Fujitsu Limited: Forest Grove, OR, USA, 2012. [Google Scholar]
  38. Chothia, C.; Lesk, A.M. The relation between the divergence of sequence and structure in proteins. EMBO J. 1986, 5, 823–826. [Google Scholar] [CrossRef]
  39. Fujitsu Limited. CAChe for Windows User Guide; Fujitsu Limited: Tokyo, Japan; Oxford Molecular Ltd.: Oxford, UK, 2003. [Google Scholar]
  40. Navratna, V.; Tosh, D.K.; Jacobson, K.A.; Gouaux, E. Thermostabilization and purification of the human dopamine transporter (hDAT) in an inhibitor and allosteric ligand bound conformation. PLoS ONE 2018, 13, e0200085. [Google Scholar] [CrossRef] [Green Version]
  41. Mattson, R.J.; Catt, J.D.; Denhart, D.J.; Deskus, J.A.; Ditta, J.L.; Higgins, M.A.; Marcin, L.R.; Sloan, C.P.; Beno, B.R.; Gao, Q. Conformationally restricted homotryptamines. 2. Indole cyclopropylmethylamines as selective serotonin reuptake inhibitors. J. Med. Chem. 2005, 48, 6023–6034. [Google Scholar] [CrossRef]
  42. Owens, M.J.; Krulewicz, S.; Simon, J.S.; Sheehan, D.V.; Thase, M.E.; Carpenter, D.J.; Plott, S.J.; Nemeroff, C.B. Estimates of serotonin and norepinephrine transporter inhibition in depressed patients treated with paroxetine or venlafaxine. Neuropsychopharmacology 2008, 33, 3201–3212. [Google Scholar] [CrossRef]
  43. Tatsumi, M.; Groshan, K.; Blakely, R.D.; Richelson, E. Pharmacological profile of antidepressants and related compounds at human monoamine transporters. Eur. J. Pharmacol. 1997, 340, 249–258. [Google Scholar] [CrossRef]
  44. Andersen, J.; Ringsted, K.B.; Bang-Andersen, B.; Strømgaard, K.; Kristensen, A.S. Binding site residues control inhibitor selectivity in the human norepinephrine transporter but not in the human dopamine transporter. Sci. Rep. 2015, 5, 15650. [Google Scholar] [CrossRef] [Green Version]
  45. Skovstrup, S.; Taboureau, O.; Bräuner-Osborne, H.; Jørgensen, F.S. Homology modelling of the GABA transporter and analysis of tiagabine binding. ChemMedChem 2010, 5, 986–1000. [Google Scholar] [CrossRef]
  46. Abramyan, A.M.; Slack, R.D.; Meena, S.; Davis, B.A.; Newman, A.H.; Singh, S.K.; Shi, L. Computation-guided analysis of paroxetine binding to hSERT reveals functionally important structural elements and dynamics. Neuropharmacology 2019, 161, 107411. [Google Scholar] [CrossRef]
  47. Kulkarni, S.S.; Grundt, P.; Kopajtic, T.; Katz, J.L.; Newman, A.H. Structure—activity relationships at monoamine transporters for a series of N-substituted 3α-(bis [4-fluorophenyl] methoxy) tropanes: Comparative molecular field analysis, synthesis, and pharmacological evaluation. J. Med. Chem. 2004, 47, 3388–3398. [Google Scholar] [CrossRef]
  48. Beuming, T.; Kniazeff, J.; Bergmann, M.L.; Shi, L.; Gracia, L.; Raniszewska, K.; Newman, A.H.; Javitch, J.A.; Weinstein, H.; Gether, U. The binding sites for cocaine and dopamine in the dopamine transporter overlap. Nat. Neurosci. 2008, 11, 780–789. [Google Scholar] [CrossRef] [Green Version]
  49. Rounsaville, B.J. Treatment of cocaine dependence and depression. Biol. Psychiatry 2004, 56, 803–809. [Google Scholar] [CrossRef]
  50. Xu, L.; Chen, L.Y. Identification of a new allosteric binding site for cocaine in dopamine transporter. J. Chem. Inf. Model. 2020, 60, 3958–3968. [Google Scholar] [CrossRef]
  51. Wilson, M.; Tripp, J. Clomipramine. Statpearls. 2020. Available online: https://www.ncbi.nlm.nih.gov/books/NBK541006/ (accessed on 1 June 2021).
  52. O’Donnell, J.; Zeppenfeld, D.; McConnell, E.; Pena, S.; Nedergaard, M. Norepinephrine: A neuromodulator that boosts the function of multiple cell types to optimize CNS performance. Neurochem. Res. 2012, 37, 2496–2512. [Google Scholar] [CrossRef] [Green Version]
  53. Eliwa, H.; Belzung, C.; Surget, A. Adult hippocampal neurogenesis: Is it the alpha and omega of antidepressant action? Biochem. Pharmacol. 2017, 141, 86–99. [Google Scholar] [CrossRef]
  54. Pidathala, S.; Mallela, A.K.; Joseph, D.; Penmatsa, A. Structural basis of norepinephrine recognition and transport inhibition in neurotransmitter transporters. bioRxiv 2020. [Google Scholar] [CrossRef]
  55. Altova, E.P.; Rykov, A.N.; Abaev, M.A.; Marochkin, I.I.; Vilkova, A.L.; Popik, M.V.; Shishkov, I.F. Molecular structure of noradrenaline according to gas-phase electron diffraction data and quantum chemical calculations. Mendeleev Commun. 2019, 29, 411–413. [Google Scholar] [CrossRef]
  56. Andersen, A. Structural studies of metabolic products of dopamine. IV. Crystal and molecular structure of (—)-noradrenaline. Acta Chem. Scand. 1975, 29, 871–876. [Google Scholar] [CrossRef]
  57. Pandit-Taskar, N.; Modak, S. Norepinephrine transporter as a target for imaging and therapy. J. Nucl. Med. 2017, 58, 39S–53S. [Google Scholar] [CrossRef] [Green Version]
  58. Moret, C.; Briley, M. The importance of norepinephrine in depression. Neuropsychiatry Dis. Treat. 2011, 7, 9. [Google Scholar]
  59. Brambilla, P.; Perez, J.; Barale, F.; Schettini, G.; Soares, J.C. GABAergic dysfunction in mood disorders. Mol. Psychiatry 2003, 8, 721–737. [Google Scholar] [CrossRef] [Green Version]
  60. Kickinger, S.; Hellsberg, E.; Schousboe, A.; Ecker, G.F.; Wellendorph, P. Structural and molecular aspects of betaine-GABA transporter 1 (BGT1) and its relation to brain function. Neuropharmacology 2019, 161, 107644. [Google Scholar] [CrossRef]
  61. Kvist, T.; Christiansen, B.; Jensen, A.A.; Brauner-Osborne, H. The four human γ-aminobutyric acid (GABA) transporters: Pharmacological characterization and validation of a highly efficient screening assay. Comb. Chem. High Throughput Screen. 2009, 12, 241–249. [Google Scholar] [CrossRef]
  62. Sliwoski, G.; Kothiwale, S.; Meiler, J.; Lowe, E.W. Computational methods in drug discovery. Pharmacol. Rev. 2014, 66, 334–395. [Google Scholar] [CrossRef] [Green Version]
  63. Shah, K.; Mujwar, S.; Gupta, J.K.; Shrivastava, S.K.; Mishra, P. Molecular Docking and In Silico Cogitation Validate Mefenamic Acid Prodrugs as Human Cyclooxygenase-2 Inhibitor. Assay Drug Dev. Technol. 2019, 17, 285–291. [Google Scholar] [CrossRef]
  64. Stahl, S.M. Serotonin: It's possible to have too much of a good thing. J. Clin. Psychiatry 1997, 58, 520–521. [Google Scholar] [CrossRef] [Green Version]
  65. Möhler, H. The GABA system in anxiety and depression and its therapeutic potential. Neuropharmacology 2012, 62, 42–53. [Google Scholar] [CrossRef]
  66. Wang, S.; Du, L.; Peng, G.; Li, W. GABA inhibits proliferation and self-renewal of mouse retinal progenitor cell. Cell Death Discov. 2019, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  67. Park, S.-C. Neurogenesis and antidepressant action. Cell Tissue Res. 2019, 377, 95–106. [Google Scholar] [CrossRef] [PubMed]
  68. Earnheart, J.C.; Schweizer, C.; Crestani, F.; Iwasato, T.; Itohara, S.; Mohler, H.; Lüscher, B. GABAergic control of adult hippocampal neurogenesis in relation to behavior indicative of trait anxiety and depression states. J. Neurosci. 2007, 27, 3845–3854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. De Oliveira, C.L.; Bolzan, J.A.; Surget, A.; Belzung, C. Do antidepressants promote neurogenesis in adult hippocampus? A systematic review and meta-analysis on naive rodents. Pharmacol. Ther. 2020, 210, 107515. [Google Scholar] [CrossRef]
  70. Huffman, J.; Taylor, G.T. Stress, Neurogenesis, and Mood: An Introduction to the Neurogenic Hypothesis of Depression; Wiley: New York, NY, USA, 2020; pp. 37–42. [Google Scholar]
  71. Kempermann, G.; Song, H.; Gage, F.H. Neurogenesis in the adult hippocampus. Cold Spring Harb. Perspect. Biol. 2015, 7, a018812. [Google Scholar] [CrossRef] [Green Version]
  72. Catavero, C.; Bao, H.; Song, J. Neural mechanisms underlying GABAergic regulation of adult hippocampal neurogenesis. Cell Tissue Res. 2018, 371, 33–46. [Google Scholar] [CrossRef]
  73. Zhao, C.; Deng, W.; Gage, F.H. Mechanisms and functional implications of adult neurogenesis. Cell 2008, 132, 645–660. [Google Scholar] [CrossRef] [Green Version]
  74. Pabba, M.; Sibille, E. GABA, depression and suicide. In Biological Aspects of Suicidal Behavior; Karger Publishers: Basel, Switzerland, 2016; Volume 30, pp. 37–50. [Google Scholar]
  75. Pytka, K.; Dziubina, A.; Młyniec, K.; Dziedziczak, A.; Żmudzka, E.; Furgała, A.; Olczyk, A.; Sapa, J.; Filipek, B. The role of glutamatergic, GABA-ergic, and cholinergic receptors in depression and antidepressant-like effect. Pharmacol. Rep. 2016, 68, 443–450. [Google Scholar] [CrossRef]
  76. Frazer, A.; Benmansour, S. Delayed pharmacological effects of antidepressants. Mol. Psychiatry 2002, 7, S23–S28. [Google Scholar] [CrossRef]
  77. Snyder, J.S. Recalibrating the relevance of adult neurogenesis. Trends Neurosci. 2019, 42, 164–178. [Google Scholar] [CrossRef]
  78. Jahromi, M.; Razavi, S.; Amirpour, N.; Khosravizadeh, Z. Paroxetine can enhance neurogenesis during neurogenic differentiation of human adipose-derived stem cells. Avicenna J. Med. Biotechnol. 2016, 8, 152. [Google Scholar]
  79. Schousboe, A.; Madsen, K.K. Delineation of the Role of Astroglial GABA Transporters in Seizure Control. Neurochem. Res. 2017, 42, 2019–2023. [Google Scholar] [CrossRef]
  80. Schousboe, A.; Wellendorph, P.; Frølund, B.; Clausen, R.P.; Krogsgaard-Larsen, P. Astrocytic GABA Transporters: Pharmacological Properties and Targets for Antiepileptic Drugs. In Glial Amino Acid Transporters; Springer: Cham, Switzerland, 2017; pp. 283–296. [Google Scholar]
  81. Annalisa, S. Structure, function, and plasticity of GABA transporters. Front. Cell. Neurosci. 2014, 8, 1–14. [Google Scholar]
  82. Lie, M.E.K.; Al-Khawaja, A.; Damgaard, M.; Haugaard, A.S.; Schousboe, A.; Clarkson, A.N.; Wellendorph, P. Glial GABA Transporters as Modulators of Inhibitory Signalling in Epilepsy and Stroke. Adv. Neurobiol. 2017, 16, 137–167. [Google Scholar]
  83. Forrest, L.R.; Tavoulari, S.; Zhang, Y.-W.; Rudnick, G.; Honig, B. Identification of a chloride ion binding site in Na+/Cl−-dependent transporters. Proc. Natl. Acad. Sci. USA 2007, 104, 12761–12766. [Google Scholar] [CrossRef] [Green Version]
  84. Kempson, S.A.; Zhou, Y.; Danbolt, N.C. The betaine/GABA transporter and betaine: Roles in brain, kidney, and liver. Front. Physiol. 2014, 5, 159. [Google Scholar] [CrossRef] [Green Version]
  85. Rusakov, D.A.; Saitow, F.; Lehre, K.P.; Konishi, S. Modulation of presynaptic Ca2+ entry by AMPA receptors at individual GABAergic synapses in the cerebellum. J. Neurosci. 2005, 25, 4930–4940. [Google Scholar] [CrossRef] [Green Version]
  86. Herculano-Houzel, S. Numbers of neurons as biological correlates of cognitive capability. Curr. Opin. Behav. Sci. 2017, 16, 1–7. [Google Scholar] [CrossRef]
  87. Drachman, D.A. Do We Have Brain to Spare? AAN Enterprises: Worcester, MA, USA, 2005. [Google Scholar]
  88. Al-Khawaja, A.; Petersen, J.G.; Damgaard, M.; Jensen, M.H.; Vogensen, S.B.; Lie, M.E.K.; Kragholm, B.; Bräuner-Osborne, H.; Clausen, R.P.; Frølund, B.; et al. Pharmacological identification of a guanidine-containing β-alanine analogue with low micromolar potency and selectivity for the betaine/GABA transporter 1 (BGT1). Neurochem. Res. 2014, 39, 1988–1996. [Google Scholar] [CrossRef]
  89. Furgała-Wojas, A.; Kowalska, M.; Nowaczyk, A.; Fijałkowski, Ł.; Sałat, K. Comparison of Bromhexine and its Active Metabolite-Ambroxol as Potential Analgesics Reducing Oxaliplatin-induced Neuropathic Pain-Pharmacodynamic and Molecular Docking Studies. Curr. Drug Metab. 2020, 21, 548–561. [Google Scholar] [CrossRef]
  90. Kowalska, M.; Nowaczyk, J.; Nowaczyk, A. KV11.1, NaV1.5, and CaV1.2 Transporter Proteins as Antitarget for Drug Cardiotoxicity. Int. J. Mol. Sci. 2020, 21, 8099. [Google Scholar] [CrossRef] [PubMed]
  91. Pagadala, N.S.; Syed, K.; Tuszynski, J. Software for molecular docking: A review. Biophys. Rev. 2017, 9, 91–102. [Google Scholar] [CrossRef] [PubMed]
  92. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Zakrzewski, V.G.; Montgomery, J.A.; Stratmann, R.E.; Burant, S.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  93. Irwin, J.J.; Sterling, T.; Mysinger, M.M.; Bolstad, E.S.; Coleman, R.G. ZINC: A free tool to discover chemistry for biology. J. Chem. Inf. Model. 2012, 52, 1757–1768. [Google Scholar] [CrossRef] [PubMed]
  94. Dennington, R.; Keith, T.; Millam, J. GaussView, Version 5; Semichem Inc.: Shawnee Mission, KS, USA, 2009. [Google Scholar]
  95. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 1998, 19, 1639–1662. [Google Scholar] [CrossRef] [Green Version]
  96. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  97. Boeckmann, B.; Bairoch, A.; Apweiler, R.; Blatter, M.-C.; Estreicher, A.; Gasteiger, E.; Martin, M.J.; Michoud, K.; ‘Donovan, C.; Phan, I. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 1, 365–370. [Google Scholar] [CrossRef]
Ijms 22 06293 g001 550
Figure 1. Superimposition of the investigated docked structures of neurotransmitters: serotonin (SER, brown), dopamine (DA; blue), noradrenaline (NE, green), γ-aminobutyric acid (GABA, yellow) (a), and human monoamine transporters: serotonin transporter (hSERT, brown), dopamine transporter (hDAT, blue), norepinephrine transporter (hNET, green) (b), human transporters: serotonin transporter (hSERT, brown), dopamine transporter (hDAT, blue), norepinephrine transporter (hNET, green), and γ-Aminobutyric acid transporters 1 (hGAT1, yellow) (c).
Figure 1. Superimposition of the investigated docked structures of neurotransmitters: serotonin (SER, brown), dopamine (DA; blue), noradrenaline (NE, green), γ-aminobutyric acid (GABA, yellow) (a), and human monoamine transporters: serotonin transporter (hSERT, brown), dopamine transporter (hDAT, blue), norepinephrine transporter (hNET, green) (b), human transporters: serotonin transporter (hSERT, brown), dopamine transporter (hDAT, blue), norepinephrine transporter (hNET, green), and γ-Aminobutyric acid transporters 1 (hGAT1, yellow) (c).
Ijms 22 06293 g001
Ijms 22 06293 g002 550
Figure 2. X-ray structure of the hSERT crystal (a) with the grid box showing the ligand-binding site (PDB ID: 5i6x, 3.14 Å). Comparison between calculated (b) and crystallographic data (c) for binding pocket in PRX–hSERT complex. The red circles and ellipses identify equivalent residua (in 3D superposition of structures).
Figure 2. X-ray structure of the hSERT crystal (a) with the grid box showing the ligand-binding site (PDB ID: 5i6x, 3.14 Å). Comparison between calculated (b) and crystallographic data (c) for binding pocket in PRX–hSERT complex. The red circles and ellipses identify equivalent residua (in 3D superposition of structures).
Ijms 22 06293 g002
Ijms 22 06293 g003 550
Figure 3. X-ray structure of the hDAT homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: Q01959) and calculated data for cocaine binding pocket (b) and clomipramine binding pocket in PRX–hNET complex (c). The red circles and ellipses identify equivalent residua (in 3D superposition of structures).
Figure 3. X-ray structure of the hDAT homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: Q01959) and calculated data for cocaine binding pocket (b) and clomipramine binding pocket in PRX–hNET complex (c). The red circles and ellipses identify equivalent residua (in 3D superposition of structures).
Ijms 22 06293 g003
Ijms 22 06293 g004 550
Figure 4. X-ray structure of the hNET homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: P23975) and calculated data (b) for binding pocket in PRX–hSERT complex.
Figure 4. X-ray structure of the hNET homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: P23975) and calculated data (b) for binding pocket in PRX–hSERT complex.
Ijms 22 06293 g004
Ijms 22 06293 g005 550
Figure 5. X-ray structure of the hGAT1 homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: P30531) and calculated data (b) for binding pocket in PRX–hGAT1 and (c) GABA–hGAT1 complex. The red circles and ellipses identify equivalent residues (in 3D superposition of structures).
Figure 5. X-ray structure of the hGAT1 homology model (a) with the grid box showing the ligand-binding site (Swiss-Prot ID: P30531) and calculated data (b) for binding pocket in PRX–hGAT1 and (c) GABA–hGAT1 complex. The red circles and ellipses identify equivalent residues (in 3D superposition of structures).
Ijms 22 06293 g005
Ijms 22 06293 g006 550
Figure 6. Validation redocking experiment. Docking results of the hSERT (a), hNET (b), hDAT/cocaine (c), hGAT1 (d) obtained using AutoDock (gray) and Vina (blue).
Figure 6. Validation redocking experiment. Docking results of the hSERT (a), hNET (b), hDAT/cocaine (c), hGAT1 (d) obtained using AutoDock (gray) and Vina (blue).
Ijms 22 06293 g006
Ijms 22 06293 g007 550
Figure 7. Inhibitory activity of PRX, SER, NE and DA at hGATs. The compounds were tested for their ability to inhibit the uptake of 30 nM [3H]GABA for 3 min at all hGATs stably expressed in CHO cells. Representative concentration–response curves of PRX at hGAT1-3 and hBGT1 (a) and inhibitory activity of SER, NE and DA at hGAT1-3 and hBGT1 (b). All experiments were performed with technical replicates in three independent experiments and depicted as means ± S.D in (a) or normalized means ± S.E.M (b).
Figure 7. Inhibitory activity of PRX, SER, NE and DA at hGATs. The compounds were tested for their ability to inhibit the uptake of 30 nM [3H]GABA for 3 min at all hGATs stably expressed in CHO cells. Representative concentration–response curves of PRX at hGAT1-3 and hBGT1 (a) and inhibitory activity of SER, NE and DA at hGAT1-3 and hBGT1 (b). All experiments were performed with technical replicates in three independent experiments and depicted as means ± S.D in (a) or normalized means ± S.E.M (b).
Ijms 22 06293 g007
Table
Table 1. The summary of hSERT, hDAT, hNET, hGAT1 and PRX, neurotransmitters docking experiment results.
Table 1. The summary of hSERT, hDAT, hNET, hGAT1 and PRX, neurotransmitters docking experiment results.
ComplexEBpKiAmino Acid ResiduesHBAngleLHBEHB
ProteinLigandkcal/molDonorAccθÅkcal/mol
hSERTPRX−13.9210.20ASP98%NH2+# COO139.952.76−1.68
SER−9.77.11ALA96%NH1+# CONH129.521.91−0.27
ASP98%NH2+# COO132.831.78−0.15
SER336%NH3+# OH164.351.86−2.54
SER438%IndOH# CONH138.771.92−0.27
R−NE−5.453.99ASP98%OH# COO134.191.81−0.12
PHE335%NH3+# CONH174.952.09−3.46
GLU493%m-Ph-OH# COO134.941.85−0.05
S−NE−5.313.89TYR95# Ph-OH%m-PH-OH156.361.77−5.21
ASP98%NH3+# COO147.641.81−1.97
ASP98%OH# COO141.191.64−6.26
DA−5.453.99PHE335%NH3+# CONH159.932.11−0.87
GLU493%p-Ph-OH# COO154.891.96−2.34
GLU493%m-PH-OH# COO155.262.01−5.80
GABA−6.773.86LYS490# NH1+%COO1150.891.68−0.62
LYS490# NH2+%COO2126.722.24−0.05
hDAT/cocainePRX−8.356.12ASP476%NH2+# COO144.542.64−0.02
SER−7.775.70ASP79%NH1+# COO142.741.99−0.12
ASN157%IndOH# CONH2120.702.17−0.23
SER422%NH2+# CONH135.292.03−2.57
R−NE−6.014.40PHE76%OH# CONH154.982.12−1.02
ALA77%NH3+# CONH139.521.85−1.81
ASP79%NH3+# COO149.091.81−0.38
ASP79%m-PH-OH# COO145.332.15−1.00
S−NE−6.014.41ASP79%p-PH-OH# COO168.081.89−3.85
ASP79%m-PH-OH# COO153.471.73−1.02
ASP476%OH# COO122.981.95−0.26
ASP476%NH3+# COO127.312.12−0.01
ASP476%NH3+# COO138.161.78−0.12
DA−7.225.29ASP79%NH1# COO131.462.12−2.93
SER149%p-Ph-OH# CONH158.762.19−3.40
TYR156%NH2# Ph-OH159.812.120.032
GABA−4.623.38ASP79%NH1+# COO144.111.74−1.10
hDAT/clomipraminePRX−8.025.88ASP476%NH2+# COO155.892.75−0.04
SER−7.055.17ASP476%NH1+# COO124.042.19−0.08
ASP476%NH2+# COO141.211.87−0.51
ASP476%IndOH# CONH165.741.83−4.42
R−NE−6.044.43ASP385%p-Ph-OH# COO159.511.77−0.46
ASP385%m-PH-OH# COO171.741.96−2.65
ASP476%OH# COO156.161.96−1.07
S−NE−6.074.45ASP385%p-Ph-OH# COO160.651.83−1.50
ASP385%m-PH-OH# COO158.362.11−1.50
ASP476%OH# COO168.501.84−1.29
DA−6.724.92ASP385%p-Ph-OH# COO156.931.87−1.38
ASP385%m-PH-OH# COO162.142.00−4.85
ASP476%NH3+# COO127.691.91−0.10
GABA−4.343.18ASN93%NH1+# CONH2126.491.99−0.85
SER309%NH2+# OH137.672.15−2.84
hNETPRX−10.387.61ALA73%NH1+# CONH125.522.73−1.06
ASP75%NH2+# COO162.172.74−0.12
SER−5.914.33PHE72%NH1+# CONH152.372.18−4.06
ASP75%NH2+# COO142.292.02−0.04
SER419%NH3+# OH135.901.96−2.85
SER420%IndOH# OH168.892.20−4.18
R−NE−5.624.12ASP75%p-Ph-OH# COO170.951.72−3.11
ASP75%m-PH-OH# COO153.611.74−1.64
ASP473%OH# COO143.252.01−0.01
S−NE−7.755.59PHE72%NH1+# CONH121.162.06−0.56
ASP75%NH2+# COO153.532.19−0.53
SER318%OH# CONH135.061.94−1.34
LEU319%m-PH-OH# CONH131.152.14−0.61
DA−7.255.31ASP75%NH1+# COO121.361.99−0.13
ALA145%p-Ph-OH# CONH134.72.11−0.93
TYR152%NH2+# Ph-OH169.392.18−3.64
SER419%NH3+# OH139.352.25−2.08
GABA−4.513.31ASP75%NH1+# COO151.091.81−1.60
hGAT1PRX−5.934.35ASP451%NH2+# CONH151.562.91−3.70
ASP451%NH2+# COO147.392.93−3.13
SER−4.863.56TYR139%NH1# Ph-OH172.421.99−5.78
ASP451%NH2# COO143.111.86−2.83
ASP451%NH3# CONH122.091.870.023
SER456%m-PH-OH# CONH139.312.06−2.94
MET458# CONH%m-PH-OH150.791.90−4.45
SER459# CONH%m-PH-OH164.482.08−5.04
R−NE−5.283.87TYR139%NH1# Ph-OH130.792.21−1.60
ASP451%NH2# CONH148.661.85−1.96
MET458# CONH%p-PH-OH163.771.99−5.67
SER459# CONH%p-PH-OH151.541.93−4.51
LEU460# CONH%m-PH-OH151.032.24−2.73
S−NE−5.283.86ASP451%NH1# COO135.121.76−0.47
ASP451%OH# CONH171.271.65−2.93
SER456%m-PH-OH# CONH152.492.06−3.09
MET458# CONH%m-PH-OH153.651.78−4.80
SER459# CONH%m-PH-OH162.422.005−5.5
DA−4.853.55ASP451%NH1# COO123.391.941−0.503
ASP451%NH2# CONH127.832.017−1.36
SER456%m-PH-OH# CONH172.051.937−4.819
MET458# CONH%m-PH-OH151.592.026−4.052
SER459# CONH%m-PH-OH162.662.124−4.493
LEU460# CONH%p-PH-OH155.461.797−5.329
GABA−6.244.49ASP451%NH1+# COO142.522.503−0.026
ASP451%NH2+# CONH156.842.451−3.32
HB—hydrogen bond, Acc—hydrogen bond acceptor, Hydrogen bond components: from the ABX % and from the protein #, EB—complex energy binding, θ—hydrogen bond angle, LHB—hydrogen bond length, EHB—hydrogen bond energy, pKi—calculated binding affinity.
Table
Table 2. The composition of the binding pocket of the analyzed channel models: hSERT, hDAT, hNET, hGAT1.
Table 2. The composition of the binding pocket of the analyzed channel models: hSERT, hDAT, hNET, hGAT1.
ProteinResidues
hSERTTyr95, Ala96, Asp98, Gly100, Ala169, Ile172, Ala173, Tyr176, Phe335, Ser336, Gly338, Phe341, Val343, Ser438, Thr439, Ala441, Gly442, Val489, Lys490, Glu493, Glu494, Thr497, Gly498, Pro499, Leu502 [6]
hNETAla145, Tyr151, Ile315, Phe316, Ser420, Ala426 [44].
hDATAla81, Tyr88, Asp385, Asp385, Asp476, Ala480, Phe472, Thr473, Asp476, His477, Ala480, Gly481, Thr482, Leu485 [6]
hGAT1Tyr60, Ala61, Gly63, Gly65, Leu136, Tyr140, Phe294, Ser295, Tyr296, Gly297, Leu300, Thr400 [45]
Table
Table 3. Inhibitory activity of PRX at hGAT1-3 and hBGT1 stably expressed in CHO cells using the [3H]GABA uptake assay.
Table 3. Inhibitory activity of PRX at hGAT1-3 and hBGT1 stably expressed in CHO cells using the [3H]GABA uptake assay.
IC50 (pIC50 ± S.E.M.) (µM)
hGAT1hGAT2hGAT3hBGT1
PRX256.1 (3.6 ± 0.03)122.4 (3.9 ± 0.08)110.6 (4.0 ± 0.04)85.6 (4.1 ± 0.05)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kowalska, M.; Fijałkowski, Ł.; Nowaczyk, A. Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters. Int. J. Mol. Sci. 2021, 22, 6293. https://doi.org/10.3390/ijms22126293

AMA Style

Kowalska M, Fijałkowski Ł, Nowaczyk A. Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters. International Journal of Molecular Sciences. 2021; 22(12):6293. https://doi.org/10.3390/ijms22126293

Chicago/Turabian Style

Kowalska, Magdalena, Łukasz Fijałkowski, and Alicja Nowaczyk. 2021. "Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters" International Journal of Molecular Sciences 22, no. 12: 6293. https://doi.org/10.3390/ijms22126293

APA Style

Kowalska, M., Fijałkowski, Ł., & Nowaczyk, A. (2021). Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters. International Journal of Molecular Sciences, 22(12), 6293. https://doi.org/10.3390/ijms22126293

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop