Microseconds simulations reveal a new sodium-binding site and the mechanism of sodium-coupled substrate uptake by LeuT

Molecular mechanisms of Na+-driven bile acid transport in human NTCP

Human Na+ taurocholate co-transporting protein (hNTCP) is a key bile salt transporter to maintain enterohepatic circulation and is responsible for the recognition of hepatitis B and D viruses. Despite landmark cryoelectron microscopy studies revealing open-pore and inward-facing states of hNTCP stabilized by antibodies, the transport mechanism remains largely unknown. To address this knowledge gap, we used molecular dynamics and enhanced sampling metadynamics simulations to elucidate the intrinsic mechanism of hNTCP-mediated taurocholate acid (TCA) transport driven by Na+ binding. We uncovered three TCA-binding modes, including one that closely matched the limited cryoelectron microscopy density observed in the open-pore hNTCP. We also captured several key hNTCP conformations in the substrate transport cycle, particularly including an outward-facing, substrate-bound state. Furthermore, we provided thermodynamic evidence supporting that changes in the Na+-binding state drive the TCA transport by exploiting the amphiphilic nature of the substrate and modulating the protein environment, thereby enabling the TCA molecule to flip through. Understanding these mechanistic details of Na+-driven bile acid transport may aid in the development of hNTCP-targeted therapies for liver diseases.

Sodium/Iodide Symporter Metastable Intermediates Provide Insights into Conformational Transition between Principal Thermodynamic States

The Sodium/Iodide Symporter (NIS), a 13-helix transmembrane protein found in the thyroid and other tissues, transports iodide, a required constituent of thyroid hormones T3 and T4. Despite extensive experimental information and clinical data, structural details of the intermediate microstates comprising the conformational transition of NIS between its inwardly and outwardly open states remain unresolved. We present data from a combination of enhanced sampling and transition path molecular dynamics (MD) simulations that elucidate the principal intermediate states comprising the inwardly to outwardly open transition of fully bound and apo NIS under an enforced ionic gradient. Our findings suggest that in both the absence and presence of bound physiological ions, NIS principally occupies a proximally inward to inwardly open state. When fully bound, NIS is also found to occupy a rare "inwardly occluded" state. The results of this work provide novel insight into the populations of NIS intermediates and the free energy landscape comprising the conformational transition, adding to a mechanistic understanding of NIS ion transport. Moreover, the knowledge gained from this approach can serve as a basis for studies of NIS mutants to target therapeutic interventions.

Structure-based virtual screening and molecular dynamics of potential inhibitors targeting sodium-bile acid co-transporter of carcinogenic liver fluke Clonorchis sinensis

Background

Clonorchis sinensis requires bile acid transporters as this fluke inhabits bile juice-filled biliary ducts, which provide an extreme environment. Clonorchis sinensis sodium-bile acid co-transporter (CsSBAT) is indispensable for the fluke’s survival in the final host, as it circulates taurocholate and prevents bile toxicity in the fluke; hence, it is recognized as a useful drug target.

Methodology and principal findings

In the present study, using structure-based virtual screening approach, we presented inhibitor candidates targeting a bile acid-binding pocket of CsSBAT. CsSBAT models were built using tertiary structure modeling based on a bile acid transporter template (PDB ID: 3zuy and 4n7x) and were applied into AutoDock Vina for competitive docking simulation. First, potential compounds were identified from PubChem (holding more than 100,000 compounds) by applying three criteria: i) interacting more favorably with CsSBAT than with a human homolog, ii) intimate interaction to the inward- and outward-facing conformational states, iii) binding with CsSBAT preferably to natural bile acids. Second, two compounds were identified following the Lipinski’s rule of five. Third, other two compounds of molecular weight higher than 500 Da (Mr > 500 Da) were presumed to efficiently block the transporter via a feasible rational screening strategy. Of these candidates, compound 9806452 exhibited the least hepatotoxicity that may enhance drug-likeness properties.

Conclusions

It is proposed that compound 9806452 act as a potential inhibitor toward CsSBAT and further studies are warranted for drug development process against clonorchiasis.

Sodium Binding Stabilizes the Outward-Open State of SERT by Limiting Bundle Domain Motions

The human serotonin transporter (hSERT) removes the neurotransmitter serotonin from the synaptic cleft by reuptake into the presynaptic nerve terminal. A number of neurologic diseases are associated with dysfunction of the hSERT, and several medications for their treatment are hSERT blockers, including citalopram, fluoxetine, and paroxetine. The substrate transport is energized by the high concentration of external NaCl. We showed through molecular dynamics simulations that the binding of NaCl stabilized the hSERT in the substrate-binding competent conformation, which was characterized by an open access path to the substrate-binding site through the outer vestibule. Importantly, the binding of NaCl reduced the dynamics of the hSERT by decreasing the internal fluctuations of the bundle domain as well as the movement of the bundle domain relative to the scaffold domain. In contrast, the presence of only the bound chloride ion did not reduce the high domain mobility of the apo state.

Investigating the Mechanism of Sodium Binding to SERT Using Direct Simulations

The serotonin transporter (SERT) terminates neurotransmission by transporting serotonin from the synapse into the pre-synaptic nerve terminal. Altered SERT function leads to several neurological diseases including depression, anxiety, mood disorders, and attention deficit hyperactivity disorders (ADHD). Accordingly SERT is the target for their pharmacological treatments, but also targeted by multiple drugs of abuse. Transport of serotonin by SERT is energized by the transmembrane electrochemical gradient of sodium. We used extensive molecular dynamics simulations to investigate the process of sodium binding to SERT, which is the first step in the transport cycle that leads to serotonin uptake. Comparing data from 51 independent simulations, we find a remarkably well-defined path for sodium entry and could identify two transient binding sites, while observing binding kinetics that are comparable to experimental data. Importantly, the structure and dynamics of the sodium binding sites indicate that sodium binding is accompanied by an induced-fit mechanism that leads to new conformations and reduces local dynamics.

Identification of multiple substrate binding sites in SLC4 transporters in the outward-facing conformation: Insights into the transport mechanism

Solute carrier family 4 (SLC4) transporters mediate the transmembrane transport of HCO₃^-, CO₃^2-, and Cl^- necessary for pH regulation, transepithelial H⁺/base transport, and ion homeostasis. Substrate transport with varying stoichiometry and specificity is achieved through an exchange mechanism and/or through coupling of the uptake of anionic substrates to typically co-transported Na⁺. Recently solved outward-facing structures of two SLC4 members (human anion exchanger 1 [hAE1] and human electrogenic sodium bicarbonate cotransporter 1 [hNBCe1]) with different transport modes (Cl^-/HCO₃^- exchange versus Na⁺-CO₃^2- symport) revealed highly conserved three-dimensional organization of their transmembrane domains. However, the exact location of the ion binding sites and their protein-ion coordination motifs are still unclear. In the present work, we combined site identification by ligand competitive saturation mapping and extensive molecular dynamics sampling with functional mutagenesis studies which led to the identification of two substrate binding sites (entry and central) in the outward-facing states of hAE1 and hNBCe1. Mutation of residues in the identified binding sites led to impaired transport in both proteins. We also showed that R730 in hAE1 is crucial for anion binding in both entry and central sites, whereas in hNBCe1, a Na⁺ acts as an anchor for CO₃^2- binding to the central site. Additionally, protonation of the central acidic residues (E681 in hAE1 and D754 in hNBCe1) alters the ion dynamics in the permeation cavity and may contribute to the transport mode differences in SLC4 proteins. These results provide a basis for understanding the functional differences between hAE1 and hNBCe1 and may facilitate potential drug development for diseases such as proximal and distal renal tubular acidosis.

Crystal structures of LeuT reveal conformational dynamics in the outward-facing states

The neurotransmitter:sodium symporter (NSS) homolog LeuT from Aquifex aeolicus has proven to be a valuable model for studying the transport mechanism of the NSS family. Crystal structures have captured LeuT in key conformations visited during the transport cycle, allowing for the construction of a nearly complete model of transport, with much of the conformational dynamics studied by computational simulations. Here, we report crystal structures of LeuT representing new intermediate conformations between the outward-facing open and occluded states. These structures, combined with binding and accessibility studies, reveal details of conformational dynamics that can follow substrate binding at the central substrate binding site (S1) of LeuT in outward-facing states, suggesting a potential competition for direction between the outward-open and outward-occluded states at this stage during substrate transport. Our structures further support an intimate interplay between the protonation state of Glu290 and binding of Na1 that may ultimately regulate the outward-open-to-occluded transition.

Interaction of Camptothecin with Model Cellular Membranes

Accurate and efficient prediction of drug partitioning in model membranes is of significant interest to the pharmaceutical industry. Herein, we utilize advanced sampling methods, specifically, the adaptive biasing force methodology to calculate the potential of mean force for a model hydrophobic anticancer drug, camptothecin (CPT), across three model interfaces. We consider an octanol bilayer, a thick octanol/water interface, and a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/water interface. We characterize the enthalpic and entropic contributions of the drug to the potential of mean force. We show that the rotational entropy of the drug is inversely related to the probability of hydrogen bond formation of the drug with the POPC membrane. In addition, in long-time microsecond simulations of a high concentration of CPT above the POPC membrane, we show that strong drug-drug aromatic interactions shift the spatial orientation of the drug with the membrane. Stacks of hydrophobic drugs form, allowing penetration of the drug just under the POPC head groups. These results imply that inhomogeneous membrane models need to take into account the effect of drug aggregation on the membrane environment.

Computational Dissection of Membrane Transport at a Microscopic Level

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Locking Two Rigid-body Bundles in an Outward-Facing Conformation: The Ion-coupling Mechanism in a LeuT-fold Transporter

Secondary active transporters use electrochemical gradient of ions to fuel the "uphill" translocation of the substrate following the alternating-access model. The coupling of ions to conformational dynamics of the protein remains one of the least characterized aspects of the transporter function. We employ extended molecular dynamics (MD) simulations to examine the Na+-binding effects on the structure and dynamics of a LeuT-fold, Na+-coupled secondary transporter (Mhp1) in its major conformational states, i.e., the outward-facing (OF) and inward-facing (IF) states, as well as on the OF ↔ IF state transition. Microsecond-long, unbiased MD simulations illustrate that Na+ stabilizes an OF conformation favorable for substrate association, by binding to a highly conserved site at the interface between the two helical bundles and restraining their relative position and motion. Furthermore, a special-protocol biased simulation for state transition suggests that Na+ binding hinders the OF ↔ IF transition. These synergistic Na+-binding effects allosterically couple the ion and substrate binding sites and modify the kinetics of state transition, collectively increasing the lifetime of an OF conformation with high substrate affinity, thereby facilitating substrate recruitment from a low-concentration environment. Based on the similarity between our findings for Mhp1 and experimental reports on LeuT, we propose that this model may represent a general Na+-coupling mechanism among LeuT-fold transporters.

Serotonin transport in the 21st century

Rudnick and Sandtner review the history of serotonin transporter research in light of structural and electrophysiological advances.

Serotonin (5-hydroxytryptamine [5-HT]) is accumulated within nerve endings by the serotonin transporter (SERT), which terminates its extracellular action and provides cytoplasmic 5-HT for refilling of synaptic vesicles. SERT is the target for many antidepressant medications as well as psychostimulants such as cocaine and ecstasy (3,4-methylenedioxymethamphetamine). SERT belongs to the SLC6 family of ion-coupled transporters and is structurally related to several other transporter families. SERT was studied in the 1970s and 1980s using membrane vesicles isolated from blood platelets. These studies led to a proposed stoichiometry of transport that has been challenged by high-resolution structures of SERT and its homologues and by studies of SERT electrophysiology. Here, we review the original evidence alongside more recent structural and electrophysiological evidence. A self-consistent picture emerges with surprising insights into the ion fluxes that accompany 5-HT transport.

Monoamine transporters: structure, intrinsic dynamics and allosteric regulation

Monoamine transporters (MATs) regulate neurotransmission via the reuptake of dopamine, serotonin and norepinephrine from extra-neuronal regions and thus maintain neurotransmitter homeostasis. As targets of a wide range of compounds, including antidepressants, substances of abuse and drugs for neuropsychiatric and neurodegenerative disorders, their mechanism of action and their modulation by small molecules have long been of broad interest. Recent advances in the structural characterization of dopamine and serotonin transporters have opened the way for structure-based modeling and simulations, which, together with experimental data, now provide mechanistic understanding of their transport function and interactions. Here we review recent progress in the elucidation of the structural dynamics of MATs and their conformational landscape and transitions, as well as allosteric regulation mechanisms.

Na+-binding modes involved in thrombin's allosteric response as revealed by molecular dynamics simulations, correlation networks and Markov modeling

The monovalent sodium ion (Na+) is a critical modulator of thrombin. However, the mechanism of thrombin's activation by Na+ has been widely debated for more than twenty years. Details of the linkage between thrombin and Na+ remain vague due to limited temporal and spatial resolution in experiments. In this work, we combine microsecond scale atomic-detailed molecular dynamics simulations with correlation network analyses and hidden Markov modeling to probe the detailed thermodynamic and kinetic picture of Na+-binding events and their resulting allosteric responses in thrombin. We reveal that ASP189 and ALA190 comprise a stable Na+-binding site (referred as "inner" Na+-binding site) along with the previously known one (referred as "outer" Na+-binding site). The corresponding newly identified Na+-binding mode introduces significant allosteric responses in thrombin's regulatory regions by stabilizing selected torsion angles of residues responsive to Na+-binding. Our Markov model indicates that the bound Na+ prefers to transfer between the two Na+-binding sites when an unbinding event takes place. These results suggest a testable hypothesis of a substrate-driven Na+ migration (ΔG ∼ 1.7 kcal mol-1) from the "inner" Na+-binding site to the "outer" one during thrombin's catalytic activities. The binding of a Na+ ion at the "inner" Na+-binding site should be inferred as a prerequisite for thrombin's efficient recognition to the substrate, which opens a new angle for our understanding of Na+-binding's allosteric activation on thrombin and sheds light on detailed processes in thrombin's activation.

Modification of a Putative Third Sodium Site in the Glycine Transporter GlyT2 Influences the Chloride Dependence of Substrate Transport

Neurotransmitter removal from glycine-mediated synapses relies on two sodium-driven high-affinity plasma membrane GlyTs that control neurotransmitter availability. Mostly glial GlyT1 is the main regulator of glycine synaptic levels, whereas neuronal GlyT2 promotes the recycling of synaptic glycine and supplies neurotransmitter for presynaptic vesicle refilling. The GlyTs differ in sodium:glycine symport stoichiometry, showing GlyT1 a 2:1 and GlyT2 a 3:1 sodium:glycine coupling. Sodium binds to the GlyTs at two conserved Na⁺ sites: Na1 and Na2. The location of GlyT2 Na3 site remains unknown, although Glu650 has been involved in the coordination. Here, we have used comparative MD simulations of a GlyT2 model constructed by homology to the crystalized DAT from Drosophila melanogaster by placing the Na3 ion at two different locations. By combination of in silico and experimental data obtained by biochemical and electrophysiological analysis of GlyTs mutants, we provide evidences suggesting the GlyT2 third sodium ion is held by Glu-250 and Glu-650, within a region with robust allosteric properties involved in cation-specific sensitivity. Substitution of Glu650 in GlyT2 by the corresponding methionine in GlyT1 reduced the charge-to-flux ratio to the level of GlyT1 without producing transport uncoupling. Chloride dependence of glycine transport was almost abolished in this GlyT2 mutant but simultaneous substitution of Glu250 and Glu650 by neutral amino acids rescued chloride sensitivity, suggesting that protonation/deprotonation of Glu250 substitutes chloride function. The differential behavior of equivalent GlyT1 mutations sustains a GlyT2-specific allosteric coupling between the putative Na3 site and the chloride site.

Hyperekplexia-associated mutations in the neuronal glycine transporter 2

Hyperekplexia or startle disease is a dysfunction of inhibitory glycinergic neurotransmission characterized by an exaggerated startle in response to trivial tactile or acoustic stimuli. Although rare, this disorder can have serious consequences, including sudden infant death. One of the most frequent causes of hyperekplexia are mutations in the SLC6A5 gene, encoding the neuronal glycine transporter 2 (GlyT2), a key component of inhibitory glycinergic presynapses involved in synaptic glycine recycling though sodium and chloride-dependent co-transport. Most GlyT2 mutations detected so far are recessive, but two dominant missense mutations have been described. The detailed analysis of these mutations has revealed structural cues on the quaternary structure of GlyT2, and opens the possibility that novel selective pharmacochaperones have potential therapeutic effects in hyperekplexia.

How structural elements evolving from bacterial to human SLC6 transporters enabled new functional properties

Background

Much of the structure-based mechanistic understandings of the function of SLC6A neurotransmitter transporters emerged from the study of their bacterial LeuT-fold homologs. It has become evident, however, that structural differences such as the long N- and C-termini of the eukaryotic neurotransmitter transporters are involved in an expanded set of functional properties to the eukaryotic transporters. These functional properties are not shared by the bacterial homologs, which lack the structural elements that appeared later in evolution. However, mechanistic insights into some of the measured functional properties of the eukaryotic transporters that have been suggested to involve these structural elements are sparse or merely descriptive.

Results

To learn how the structural elements added in evolution enable mechanisms of the eukaryotic transporters in ways not shared with their bacterial LeuT-like homologs, we focused on the human dopamine transporter (hDAT) as a prototype. We present the results of a study employing large-scale molecular dynamics simulations and comparative Markov state model analysis of experimentally determined properties of the wild-type and mutant hDAT constructs. These offer a quantitative outline of mechanisms in which a rich spectrum of interactions of the hDAT N-terminus and C-terminus contribute to the regulation of transporter function (e.g., by phosphorylation) and/or to entirely new phenotypes (e.g., reverse uptake (efflux)) that were added in evolution.

Conclusions

The findings are consistent with the proposal that the size of eukaryotic neurotransmitter transporter termini increased during evolution to enable more functions (e.g., efflux) not shared with the bacterial homologs. The mechanistic explanations for the experimental findings about the modulation of function in DAT, the serotonin transporter, and other eukaryotic transporters reveal separate roles for the distal and proximal segments of the much larger N-terminus in eukaryotic transporters compared to the bacterial ones. The involvement of the proximal and distal segments — such as the role of the proximal segment in sustaining transport in phosphatidylinositol 4,5-bisphosphate-depleted membranes and of the distal segment in modulating efflux — may represent an evolutionary adaptation required for the function of eukaryotic transporters expressed in various cell types of the same organism that differ in the lipid composition and protein complement of their membrane environment.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0495-6) contains supplementary material, which is available to authorized users.

Quantitative Assessment of the Energetics of Dopamine Translocation by Human Dopamine Transporter

Computational evaluation of the energetics of substrate binding, transport, and release events of neurotransmitter transporters at the molecular level is a challenge, as the structural transitions of these membrane proteins involve coupled global and local changes that span time scales of several orders of magnitude, from nanoseconds to seconds. Here, we provide a quantitative assessment of the energetics of dopamine (DA) translocation through the human DA transporter (hDAT), using a combination of molecular modeling, simulation, and analysis tools. DA-binding and -unbinding events, which generally involve local configurational changes, are evaluated using free-energy perturbation or adaptive biasing force methods. The global transitions between the outward-facing state and the inward-facing state, on the other hand, require a dual-boost accelerated molecular dynamics simulation. We present results on DA-binding/unbinding energetics under different conditions, as well as the conformational energy landscape of hDAT in both DA-bound and -unbound states. The study provides a tractable method of approach for quantitative evaluation of substrate-binding energetics and efficient estimation of conformational energy landscape, in general.

Allosteric modulation of human dopamine transporter activity under conditions promoting its dimerization

The human dopamine (DA) transporter (hDAT) is a key regulator of neurotransmission and a target for antidepressants and addictive drugs. Despite the recent resolution of dDAT structures from Drosophila melanogaster, complete understanding of its mechanism of function and even information on its biological assembly is lacking. The resolved dDAT structures are monomeric, but growing evidence suggests that hDAT might function as a multimer, and its oligomerization may be relevant to addictive drug effects. Here, using structure-based computations, we examined the possible mechanisms of hDAT dimerization and its dynamics in a lipid bilayer. Using a combination of site-directed mutagenesis, DA-uptake, and cross-linking experiments that exploited the capacity of Cys-306 to form intermonomeric disulfide bridges in the presence of an oxidizing agent, we tested the effects of mutations at transmembrane segment (TM) 6 and 12 helices in HEK293 cells. The most probable structural model for hDAT dimer suggested by computations and experiments differed from the dimeric structure resolved for the bacterial homolog, LeuT, presumably because of a kink at TM12 preventing favorable monomer packing. Instead, TM2, TM6, and TM11 line the dimer interface. Molecular dynamics simulations of the dimeric hDAT indicated that the two subunits tend to undergo cooperative structural changes, both on local (extracellular gate opening/closure) and global (transition between outward-facing and inward-facing states) scales. These observations suggest that hDAT transport properties may be allosterically modulated under conditions promoting dimerization. Our study provides critical insights into approaches for examining the oligomerization of neurotransmitter transporters and sheds light on their drug modulation.

Ligand Binding in the Extracellular Vestibule of the Neurotransmitter Transporter Homologue LeuT

The human monoamine transporters (MATs) facilitate the reuptake of monoamine neurotransmitters from the synaptic cleft. MATs are linked to a number of neurological diseases and are the targets of both therapeutic and illicit drugs. Until recently, no high-resolution structures of the human MATs existed, and therefore, studies of this transporter family have relied on investigations of the homologues bacterial transporters such as the leucine transporter LeuT, which has been crystallized in several conformational states. A two-substrate transport mechanism has been suggested for this transporter family, which entails that high-affinity binding of a second substrate in an extracellular site is necessary for the substrate in the central binding site to be transported. Compelling evidence for this mechanism has been presented, however, a number of equally compelling accounts suggest that the transporters function through a mechanism involving only a single substrate and a single high-affinity site. To shed light on this apparent contradiction, we have performed extensive molecular dynamics simulations of LeuT in the outward-occluded conformation with either one or two substrates bound to the transporter. We have also calculated the substrate binding affinity in each of the two proposed binding sites through rigorous free energy simulations. Results show that substrate binding is unstable in the extracellular vestibule and the substrate binding affinity within the suggested extracellular site is very low (0.2 and 3.3 M for the two dominant binding modes) compared to the central substrate binding site (14 nM). This suggests that for LeuT in the outward-occluded conformation only a single high-affinity substrate binding site exists.

Effect of Dimerization on the Dynamics of Neurotransmitter:Sodium Symporters

Dimerization is a common feature among the members of the neurotransmitter:sodium symporter (NSS) family of membrane proteins. Yet, the effect of dimerization on the mechanism of action of NSS members is not fully understood. In this study, we examined the collective dynamics of two members of the family, leucine transporter (LeuT) and dopamine transporter (DAT), to assess the significance of dimerization in modulating the functional motions of the monomers. We used to this aim the anisotropic network model (ANM), an efficient and robust method for modeling the intrinsic motions of proteins and their complexes. Transporters belonging to the NSS family are known to alternate between outward-facing (OF) and inward-facing (IF) states, which enables the uptake and release of their substrate (neurotransmitter) respectively, as the substrate is transported from the exterior to the interior of the cell. In both LeuT and DAT, dimerization is found to alter the collective motions intrinsically accessible to the individual monomers in favor of the functional transitions (OF ↔ IF), suggesting that dimerization may play a role in facilitating transport.

Glycine Transporters in Glia Cells: Structural Studies

Glycine, besides exerting essential metabolic functions, is an important inhibitory neurotransmitter in caudal areas of the central nervous system and also a positive neuromodulator at excitatory glutamate-mediated synapses. Glial cells provide metabolic support to neurons and modulate synaptic activity. Six transporters belonging to three solute carrier families (SLC6, SLC38, and SLC7) are capable of transporting glycine across the glial plasma membrane. The unique glial glycine-selective transporter GlyT1 (SLC6) is the main regulator of synaptic glycine concentrations, assisted by the neuronal GlyT2. The five additional glycine transporters ATB^0,+, SNAT1, SNAT2, SNAT5, and LAT2 display broad amino acid specificity and have differential contributions to glial glycine transport. Glial glycine transporters are divergent in sequence but share a similar architecture displaying the 5 + 5 inverted fold originally characterized in the leucine transporter LeuT. The availability of protein crystals solved at high resolution for prokaryotic and, more recently, eukaryotic homologues of this superfamily has advanced significantly our understanding of the mechanism of glycine transport.

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.

Energy landscape of LeuT from molecular simulations

The bacterial sodium-coupled leucine transporter (LeuT) has been broadly used as a structural model for understanding the structure-dynamics-function of mammalian neurotransmitter transporters as well as other solute carriers that share the same fold (LeuT fold), as the first member of the family crystallographically resolved in multiple states: outward-facing open, outward-facing occluded, and inward-facing open. Yet, a complete picture of the energy landscape of (sub)states visited along the LeuT transport cycle has been elusive. In an attempt to visualize the conformational spectrum of LeuT, we performed extensive simulations of LeuT dimer dynamics in the presence of substrate (Ala or Leu) and co-transported Na(+) ions, in explicit membrane and water. We used both conventional molecular dynamics (MD) simulations (with Anton supercomputing machine) and a recently introduced method, collective MD, that takes advantage of collective modes of motions predicted by the anisotropic network model. Free energy landscapes constructed based on ∼40 μs trajectories reveal multiple substates occluded to the extracellular (EC) and/or intracellular (IC) media, varying in the levels of exposure of LeuT to EC or IC vestibules. The IC-facing transmembrane (TM) helical segment TM1a shows an opening, albeit to a smaller extent and in a slightly different direction than that observed in the inward-facing open crystal structure. The study provides insights into the spectrum of conformational substates and paths accessible to LeuT and highlights the differences between Ala- and Leu-bound substates.

Na+ coordination at the Na2 site of the Na+/I- symporter

The sodium/iodide symporter (NIS) mediates active I(-) transport in the thyroid-the first step in thyroid hormone biosynthesis-with a 2 Na(+): 1 I(-) stoichiometry. The two Na(+) binding sites (Na1 and Na2) and the I(-) binding site interact allosterically: when Na(+) binds to a Na(+) site, the affinity of NIS for the other Na(+) and for I(-) increases significantly. In all Na(+)-dependent transporters with the same fold as NIS, the side chains of two residues, S353 and T354 (NIS numbering), were identified as the Na(+) ligands at Na2. To understand the cooperativity between the substrates, we investigated the coordination at the Na2 site. We determined that four other residues-S66, D191, Q194, and Q263-are also involved in Na(+) coordination at this site. Experiments in whole cells demonstrated that these four residues participate in transport by NIS: mutations at these positions result in proteins that, although expressed at the plasma membrane, transport little or no I(-) These residues are conserved throughout the entire SLC5 family, to which NIS belongs, suggesting that they serve a similar function in the other transporters. Our findings also suggest that the increase in affinity that each site displays when an ion binds to another site may result from changes in the dynamics of the transporter. These mechanistic insights deepen our understanding not only of NIS but also of other transporters, including many that, like NIS, are of great medical relevance.

Conformational Dynamics on the Extracellular Side of LeuT Controlled by Na+ and K+ Ions and the Protonation State of Glu290

Ions play key mechanistic roles in the gating dynamics of neurotransmitter:sodium symporters (NSSs). In recent microsecond scale molecular dynamics simulations of a complete model of the dopamine transporter, a NSS protein, we observed a partitioning of K(+) ions from the intracellular side toward the unoccupied Na2 site of dopamine transporter following the release of the Na2-bound Na(+) Here we evaluate with computational simulations and experimental measurements of ion affinities under corresponding conditions, the consequences of K(+) binding in the Na2 site of LeuT, a bacterial homolog of NSS, when both Na(+) ions and substrate have left, and the transporter prepares for a new cycle. We compare the results with the consequences of binding Na(+) in the same apo system. Analysis of >50-μs atomistic molecular dynamics and enhanced sampling trajectories of constructs with Glu(290), either charged or neutral, point to the Glu(290) protonation state as a main determinant in the structural reconfiguration of the extracellular vestibule of LeuT in which a "water gate" opens through coordinated motions of residues Leu(25), Tyr(108), and Phe(253) The resulting water channel enables the binding/dissociation of the Na(+) and K(+) ions that are prevalent, respectively, in the extracellular and intracellular environments.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

A conserved leucine occupies the empty substrate site of LeuT in the Na(+)-free return state

Bacterial members of the neurotransmitter:sodium symporter (NSS) family perform Na(+)-dependent amino-acid uptake and extrude H(+) in return. Previous NSS structures represent intermediates of Na(+)/substrate binding or intracellular release, but not the inward-to-outward return transition. Here we report crystal structures of Aquifex aeolicus LeuT in an outward-oriented, Na(+)- and substrate-free state likely to be H(+)-occluded. We find a remarkable rotation of the conserved Leu25 into the empty substrate-binding pocket and rearrangements of the empty Na(+) sites. Mutational studies of the equivalent Leu99 in the human serotonin transporter show a critical role of this residue on the transport rate. Molecular dynamics simulations show that extracellular Na(+) is blocked unless Leu25 is rotated out of the substrate-binding pocket. We propose that Leu25 facilitates the inward-to-outward transition by compensating a Na(+)- and substrate-free state and acts as the gatekeeper for Na(+) binding that prevents leak in inward-outward return transitions.

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.

The molecular mechanism of secondary sodium symporters elucidated through the lens of the computational microscope

Transport of molecules across cellular membranes is a key biological process for normal cell function. In this review we describe current state-of-the-art knowledge on molecular mechanism of secondary active transporters obtained by molecular simulations studies.

Two Na+ Sites Control Conformational Change in a Neurotransmitter Transporter Homolog

In LeuT, a prokaryotic homolog of neurotransmitter transporters, Na(+) stabilizes outward-open conformational states. We examined how each of the two LeuT Na(+) binding sites contributes to Na(+)-dependent closure of the cytoplasmic pathway using biochemical and biophysical assays of conformation. Mutating either of two residues that contribute to the Na2 site completely prevented cytoplasmic closure in response to Na(+), suggesting that Na2 is essential for this conformational change, whereas Na1 mutants retained Na(+) responsiveness. However, mutation of Na1 residues also influenced the Na(+)-dependent conformational change in ways that varied depending on the position mutated. Computational analyses suggest those mutants influence the ability of Na1 binding to hydrate the substrate pathway and perturb an interaction network leading to the extracellular gate. Overall, the results demonstrate that occupation of Na2 stabilizes outward-facing conformations presumably through a direct interaction between Na(+) and transmembrane helices 1 and 8, whereas Na(+) binding at Na1 influences conformational change through a network of intermediary interactions. The results also provide evidence that N-terminal release and helix motions represent distinct steps in cytoplasmic pathway opening.

Spontaneous inward opening of the dopamine transporter is triggered by PIP2-regulated dynamics of the N-terminus

We present the dynamic mechanism of concerted motions in a full-length molecular model of the human dopamine transporter (hDAT), a member of the neurotransmitter/sodium symporter (NSS) family, involved in state-to-state transitions underlying function. The findings result from an analysis of unbiased atomistic molecular dynamics simulation trajectories (totaling >14 μs) of the hDAT molecule immersed in lipid membrane environments with or without phosphatidylinositol 4,5-biphosphate (PIP₂) lipids. The N-terminal region of hDAT (N-term) is shown to have an essential mechanistic role in correlated rearrangements of specific structural motifs relevant to state-to-state transitions in the hDAT. The mechanism involves PIP₂-mediated electrostatic interactions between the N-term and the intracellular loops of the transporter molecule. Quantitative analyses of collective motions in the trajectories reveal that these interactions correlate with the inward-opening dynamics of hDAT and are allosterically coupled to the known functional sites of the transporter. The observed large-scale motions are enabled by specific reconfiguration of the network of ionic interactions at the intracellular end of the protein. The isomerization to the inward-facing state in hDAT is accompanied by concomitant movements in the extracellular vestibule and results in the release of an Na⁺ ion from the Na2 site and destabilization of the substrate dopamine in the primary substrate binding S1 site. The dynamic mechanism emerging from the findings highlights the involvement of the PIP₂-regulated interactions between the N-term and the intracellular loop 4 in the functionally relevant conformational transitions that are also similar to those found to underlie state-to-state transitions in the leucine transporter (LeuT), a prototypical bacterial homologue of the NSS.

Molecular Mechanism of Dopamine Transport by Human Dopamine Transporter

Dopamine transporters (DATs) control neurotransmitter dopamine (DA) homeostasis by reuptake of excess DA, assisted by sodium and chloride ions. The recent resolution of DAT structure (dDAT) from Drosophila permits us for the first time to directly view the sequence of events involved in DA reuptake in human DAT (hDAT) using homology modeling and full-atomic microseconds accelerated simulations. Major observations are spontaneous closure of extracellular gates prompted by DA binding; stabilization of a holo-occluded intermediate; disruption of N82-N353 hydrogen bond and exposure to intracellular (IC) water triggered by Na2 dislocation; redistribution of a network of salt bridges at the IC surface in the inward-facing state; concerted tilting of IC-exposed helices to enable the release of Na(+) and Cl(-) ions; and DA release after protonation of D79. The observed time-resolved interactions confirm the conserved dynamics of LeuT-fold family, while providing insights into the mechanistic role of specific residues in hDAT.

Designing modulators of monoamine transporters using virtual screening techniques

The plasma-membrane monoamine transporters (MATs), including the serotonin (SERT), norepinephrine (NET) and dopamine (DAT) transporters, serve a pivotal role in limiting monoamine-mediated neurotransmission through the reuptake of their respective monoamine neurotransmitters. The transporters are the main target of clinically used psychostimulants and antidepressants. Despite the availability of several potent and selective MAT substrates and inhibitors the continuing need for therapeutic drugs to treat brain disorders involving aberrant monoamine signaling provides a compelling reason to identify novel ways of targeting and modulating the MATs. Designing novel modulators of MAT function have been limited by the lack of three dimensional structure information of the individual MATs. However, crystal structures of LeuT, a bacterial homolog of MATs, in a substrate-bound occluded, substrate-free outward-open, and an apo inward-open state and also with competitive and non-competitive inhibitors have been determined. In addition, several structures of the Drosophila DAT have also been resolved. Together with computational modeling and experimental data gathered over the past decade, these structures have dramatically advanced our understanding of several aspects of SERT, NET, and DAT transporter function, including some of the molecular determinants of ligand interaction at orthosteric substrate and inhibitor binding pockets. In addition progress has been made in the understanding of how allosteric modulation of MAT function can be achieved. Here we will review all the efforts up to date that has been made through computational approaches employing structural models of MATs to design small molecule modulators to the orthosteric and allosteric sites using virtual screening techniques.

Insights into the Modulation of Dopamine Transporter Function by Amphetamine, Orphenadrine, and Cocaine Binding

Human dopamine (DA) transporter (hDAT) regulates dopaminergic signaling in the central nervous system by maintaining the synaptic concentration of DA at physiological levels, upon reuptake of DA into presynaptic terminals. DA translocation involves the co-transport of two sodium ions and the channeling of a chloride ion, and it is achieved via alternating access between outward-facing (OF) and inward-facing states of DAT. hDAT is a target for addictive drugs, such as cocaine, amphetamine (AMPH), and therapeutic antidepressants. Our recent quantitative systems pharmacology study suggested that orphenadrine (ORPH), an anticholinergic agent and anti-Parkinson drug, might be repurposable as a DAT drug. Previous studies have shown that DAT-substrates like AMPH or -blockers like cocaine modulate the function of DAT in different ways. However, the molecular mechanisms of modulation remained elusive due to the lack of structural data on DAT. The newly resolved DAT structure from Drosophila melanogaster opens the way to a deeper understanding of the mechanism and time evolution of DAT–drug/ligand interactions. Using a combination of homology modeling, docking analysis, molecular dynamics simulations, and molecular biology experiments, we performed a comparative study of the binding properties of DA, AMPH, ORPH, and cocaine and their modulation of hDAT function. Simulations demonstrate that binding DA or AMPH drives a structural transition toward a functional form predisposed to translocate the ligand. In contrast, ORPH appears to inhibit DAT function by arresting it in the OF open conformation. The analysis shows that cocaine and ORPH competitively bind DAT, with the binding pose and affinity dependent on the conformational state of DAT. Further assays show that the effect of ORPH on DAT uptake and endocytosis is comparable to that of cocaine.

Functional mechanisms of neurotransmitter transporters regulated by lipid-protein interactions of their terminal loops

The physiological functions of neurotransmitter:sodium symporters (NSS) in reuptake of neurotransmitters from the synapse into the presynaptic nerve have been shown to be complemented by their involvement, together with non-plasma membrane neurotransmitter transporters, in the reverse transport of substrate (efflux) in response to psychostimulants. Recent experimental evidence implicates highly anionic phosphatidylinositol 4,5-biphosphate (PIP(2)) lipids in such functions of the serotonin (SERT) and dopamine (DAT) transporters. Thus, for both SERT and DAT, neurotransmitter efflux has been shown to be strongly regulated by the presence of PIP(2) lipids in the plasma membrane, and the electrostatic interaction of the N-terminal region of DAT with the negatively charged PIP(2) lipids. We examine the experimentally established phenotypes in a structural context obtained from computational modeling based on recent crystallographic data. The results are shown to set the stage for a mechanistic understanding of physiological actions of neurotransmitter transporters in the NSS family of membrane proteins. This article is part of a Special Issue entitled: Lipid-protein interactions.

Computational characterization of structural dynamics underlying function in active membrane transporters

Active transport of materials across the cellular membrane is one the most fundamental processes in biology. In order to accomplish this task, membrane transporters rely on a wide range of conformational changes spanning multiple time and size scales. These molecular events govern key functional aspects in membrane transporters, namely, coordinated gating motions underlying the alternating access mode of operation, and coupling of uphill transport of substrate to various sources of energy, for example, transmembrane electrochemical gradients and ATP binding and hydrolysis. Computational techniques such as molecular dynamics simulations and free energy calculations have equipped us with a powerful repertoire of biophysical tools offering unparalleled spatial and temporal resolutions that can effectively complement experimental methodologies, and therefore help fill the gap of knowledge in understanding the molecular basis of function in membrane transporters.