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

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Analysis of the Methylation Level of the DAT1 Dopamine Transporter Gene in Patients Addicted to Stimulants, Taking into Account an Analysis of Personality Traits

Drug addiction is a chronic biochemical drug use disorder that affects the human brain and behavior and leads to the uncontrolled use of legal or illicit drugs. It has been shown that three factors are involved in the development of addiction: genetic factors, a diverse environment, and the effect of medication on gene expression. The comprehensive approach and holistic analysis of the problem are due to the multigenic and multifactorial nature of addiction. Dopamine, one of the major neurotransmitters in the brain, is believed to be the "culprit" that leads to a drug abuse-induced "high". That is why, in our research, we focused mainly on the genes related to dopaminergic reuptake. In the present study, we chose methylation of the DAT1 dopamine transporter gene based on molecular reasons related to the dopaminergic theory of addiction. This study included two groups: 226 stimulant-dependent and 290 non-stimulant-dependent subjects. The analysis consisted of a case-control comparison of people addicted to psychostimulants compared to a control group of healthy and non-addicted people. There were differences in the levels of statistical significance between the groups. Our research shows lower methylation of islands 1, 9, and 14 in addicted people and greater methylation of islands 32 and 33. The difference in individual CpG methylation islands of the gene under study provides valuable information about the DNA methylation process in patients addicted to psychostimulants. Pearson's linear correlation analysis in stimulant dependence showed a negative correlation between total methylation island levels and the NEO-FFI Neuroticism scale. In subjects with neuroticism, the methylation level was statistically significantly lower. Pearson's linear correlation analysis of stimulant-dependent subjects showed a positive correlation between total methylation island levels and the NEO-FFI Openness scale and the NEO-FFI Conscientiousness scale.

Molecular mechanisms of action of stimulant novel psychoactive substances that target the high-affinity transporter for dopamine

Drug misuse is a significant social and public health problem worldwide. Misused substances exert their neurobehavioural effects through changing neural signalling within the brain, many of them leading to substance dependence and addiction in the longer term. Among drugs with addictive liability, there are illicit classical stimulants such as cocaine and amphetamine, and their more recently available counterparts known as novel psychoactive substances (NPS). Stimulants normally increase dopamine availability in the brain, including the pathway implicated in reward-related behaviour. This pattern is observed in both animal and human brain. The main biological target of stimulants, both classical and NPS, is the dopamine transporter (DAT) implicated in the dopamine-enhancing effects of these drugs. This article aims at reviewing research on the molecular mechanisms underpinning the interactions between stimulant NPS, such as benzofurans, cathinones or piperidine derivatives and DAT, to achieve a greater understanding of the core phenomena that decide about the addictive potential of stimulant NPS. As the methodology is essential in the process of experimental research in this area, we review the applications of in vitro, in vivo and in silico approaches. The latter, including molecular dynamics, attracts the focus of the present review as the method of choice in molecular and atomistic investigations of the mechanisms of addiction of stimulant NPS. Research of this kind is of interest to not only scientists but also health professionals as updated knowledge of NPS, their modes of action and health risks, is needed to tackle the challenges posed by NPS misuse.

Recent Advances and Challenges of the Drugs Acting on Monoamine Transporters

Background:

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

Objective:

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

Methods:

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

Results:

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

Conclusion:

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

Feasible Classified Models for Parkinson Disease from 99mTc-TRODAT-1 SPECT Imaging

The neuroimaging techniques such as dopaminergic imaging using Single Photon Emission Computed Tomography (SPECT) with ^99mTc-TRODAT-1 have been employed to detect the stages of Parkinson’s disease (PD). In this retrospective study, a total of 202 ^99mTc-TRODAT-1 SPECT imaging were collected. All of the PD patient cases were separated into mild (HYS Stage 1 to Stage 3) and severe (HYS Stage 4 and Stage 5) PD, according to the Hoehn and Yahr Scale (HYS) standard. A three-dimensional method was used to estimate six features of activity distribution and striatal activity volume in the images. These features were skewness, kurtosis, Cyhelsky’s skewness coefficient, Pearson’s median skewness, dopamine transporter activity volume, and dopamine transporter activity maximum. Finally, the data were modeled using logistic regression (LR) and support vector machine (SVM) for PD classification. The results showed that SVM classifier method produced a higher accuracy than LR. The sensitivity, specificity, PPV, NPV, accuracy, and AUC with SVM method were 0.82, 1.00, 0.84, 0.67, 0.83, and 0.85, respectively. Additionally, the Kappa value was shown to reach 0.68. This claimed that the SVM-based model could provide further reference for PD stage classification in medical diagnosis. In the future, more healthy cases will be expected to clarify the false positive rate in this classification model.

Modulation of Hippocampal Antioxidant Defense System in Chronically Stressed Rats by Lithium

This study examined the effects of lithium on gene expression and activity of the antioxidant enzymes copper zinc superoxide dismutase (SOD1), manganese superoxide dismutase (SOD2), catalase (CAT), glutathione peroxidase (GPx), and glutathione reductase (GR) in the hippocampus of chronically stressed rats. In addition, we examined the effects of lithium on anxiety behaviors, hippocampal concentrations of dopamine (DA) and malondialdehyde (MDA), protein levels of brain-derived neurotrophic factor (BDNF), tyrosine hydroxylase (TH), dopamine transporter (DAT), and catechol-O-methyltransferase (COMT), as well as activity of monoamine oxidase (MAO) in chronically stressed rats. The investigated parameters were quantified by real-time RT-PCR, Western blot analyses, and assays of enzyme activities. We found that lithium did not change gene expression of SOD1, CAT, GPx, and GR but decreased gene expression of SOD2 in chronically stressed rats. A very important result in this study was that lithium treatment decreased the enzyme activities of SOD1 and SOD2 but increased the enzyme activities of GPx and GR in stress condition, which indicates the control of redox balance. The reduced concentration of MDA confirms this. In addition, we found that lithium treatment decreased high protein levels of BDNF and DAT in chronically stressed rats to the level found in unstressed animals. Also, lithium treatment increased the expression of TH but decreased the enzyme activity of MAO B, which contributed to the increase of hippocampal concentration of DA in chronically stressed rats to the level of unstressed animals. Finally, lithium treatment in animals exposed to chronic stress increased the time spent in open arms. Lithium-induced modulation of hippocampal antioxidant status and attenuation of oxidative stress stabilized behavior in animals with high anxiety index. In addition, reduced oxidative stress was followed by the changes of both turnover of DA and levels of BDNF protein in chronically stressed rats treated with lithium. These findings may be important in preclinical research of the effects of lithium on oxidative stress level in pathological conditions.

Inhibitor mechanisms in the S1 binding site of the dopamine transporter defined by multi-site molecular tethering of photoactive cocaine analogs

Dopamine transporter (DAT) blockers like cocaine and many other abused and therapeutic drugs bind and stabilize an inactive form of the transporter inhibiting reuptake of extracellular dopamine (DA). The resulting increases in DA lead to the ability of these drugs to induce psychomotor alterations and addiction, but paradoxical findings in animal models indicate that not all DAT antagonists induce cocaine-like behavioral outcomes. How this occurs is not known, but one possibility is that uptake inhibitors may bind at multiple locations or in different poses to stabilize distinct conformational transporter states associated with differential neurochemical endpoints. Understanding the molecular mechanisms governing the pharmacological inhibition of DAT is therefore key for understanding the requisite interactions for behavioral modulation and addiction. Previously, we leveraged complementary computational docking, mutagenesis, peptide mapping, and substituted cysteine accessibility strategies to identify the specific adduction site and binding pose for the crosslinkable, photoactive cocaine analog, RTI 82, which contains a photoactive azide attached at the 2β position of the tropane pharmacophore. Here, we utilize similar methodology with a different cocaine analog N-[4-(4-azido-3-I-iodophenyl)-butyl]-2-carbomethoxy-3-(4-chlorophenyl)tropane, MFZ 2-24, where the photoactive azide is attached to the tropane nitrogen. In contrast to RTI 82, which crosslinked into residue Phe319 of transmembrane domain (TM) 6, our findings show that MFZ 2-24 adducts to Leu80 in TM1 with modeling and biochemical data indicating that MFZ 2-24, like RTI 82, occupies the central S1 binding pocket with the (+)-charged tropane ring nitrogen coordinating with the (-)-charged carboxyl side chain of Asp79. The superimposition of the tropane ring in the three-dimensional binding poses of these two distinct ligands provides strong experimental evidence for cocaine binding to DAT in the S1 site and the importance of the tropane moiety in competitive mechanisms of DA uptake inhibition. These findings set a structure-function baseline for comparison of typical and atypical DAT inhibitors and how their interactions with DAT could lead to the loss of cocaine-like behaviors.

The Environment Shapes the Inner Vestibule of LeuT

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

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

Properties of an inward-facing state of LeuT: conformational stability and substrate release

The leucine transporter (LeuT) is a bacterial homolog of the human monoamine transporters, which are important pharmaceutical targets. There are no high-resolution structures of the human transporters available; however, LeuT has been crystallized in several different conformational states. Recently, an inward-facing conformation of LeuT was solved revealing an unexpectedly large movement of transmembrane helix 1a (TM1a). We have performed molecular dynamics simulations of the mutated and wild-type transporter, with and without the cocrystallized Fab antibody fragment, to investigate the properties of this inward-facing conformation in relation to transport by LeuT within the membrane environment. In all of the simulations, local conformational changes with respect to the crystal structure are consistently observed, especially in TM1a. Umbrella sampling revealed a soft potential for TM1a tilting. Furthermore, simulations of inward-facing LeuT with Na(+) ions and substrate bound suggest that one of the Na(+) ion binding sites is fully disrupted. Release of alanine and the second Na(+) ion is also observed, giving insight into the final stage of the translocation process in atomistic detail.

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

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

Monoamine transporters: insights from molecular dynamics simulations

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

Insights to ligand binding to the monoamine transporters-from homology modeling to LeuBAT and dDAT

Understanding of drug binding to the human biogenic amine transporters (BATs) is essential to explain the mechanism of action of these pharmaceuticals but more importantly to be able to develop new and improved compounds to be used in the treatment of depression or drug addiction. Until recently no high resolution structure was available of the BATs and homology modeling was a necessity. Various studies have revealed experimentally validated binding modes of numerous ligands to the BATs using homology modeling. Here we examine and discuss the similarities between the binding models of substrates, antidepressants, psychostimulants, and mazindol in homology models of the human BATs and the recently published crystal structures of the Drosophila dopamine transporter and the engineered protein, LeuBAT. The comparison reveals that careful computational modeling combined with experimental data can be utilized to predict binding of molecules to proteins that agree very well with crystal structures.

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

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

Changes in dopamine transporter binding in nucleus accumbens following chronic self-administration cocaine: heroin combinations

Concurrent use of cocaine and heroin (speedball) has been shown to exert synergistic effects on dopamine neurotransmission in the nucleus accumbens (NAc), as observed by significant increases in extracellular dopamine levels and compensatory elevations in the maximal reuptake rate of dopamine. The present studies were undertaken to determine whether chronic self-administration of cocaine, heroin or a combination of cocaine:heroin led to compensatory changes in the abundance and/or affinity of high- and low-affinity DAT binding sites. Saturation binding of the cocaine analog [(125) I] 3β-(4-iodophenyl)tropan-2β-carboxylic acid methyl ester ([(125) I]RTI-55) in rat NAc membranes resulted in binding curves that were best fit to two-site binding models, allowing calculation of dissociation constant (Kd ) and binding density (Bmax ) values corresponding to high- and low-affinity DAT binding sites. Scatchard analysis of the saturation binding curves clearly demonstrate the presence of high- and low- affinity binding sites in the NAc, with low-affinity sites comprising 85 to 94% of the binding sites. DAT binding analyses revealed that self-administration of cocaine and a cocaine:heroin combination increased the affinity of the low-affinity site for the cocaine congener RTI-55 compared to saline. These results indicate that the alterations observed following chronic speedball self-administration are likely due to the cocaine component alone; thus further studies are necessary to elaborate upon the synergistic effect of cocaine:heroin combinations on the dopamine system in the NAc.

Structural dynamics of the monoamine transporter homolog LeuT from accelerated conformational sampling and channel analysis

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