Complete mapping of substrate translocation highlights the role of LeuT N-terminal segment in regulating transport cycle

Investigation of the interactions of HSA and SARS-CoV-2 papain-like protease against eugenol for novel COVID-19 drug discovery: spectroscopic and insilico study

Coronavirus family consist of a member known as SARS-CoV-2, spread drastically in 2019 (Covid-19), affecting millions of people worldwide. Till date there is no clear-clinical therapy or drug, targeted to cure this serious disease. Researches are going on to prevent this corona virus. Here, we tried to explore a novel SARS-CoV-2 papain-like protease as a potential inhibitor. Finally, eugenol was docked with this protease to find prime SARS-inhibitors. In silico studies revealed that eugenol binds to the active site of SARS-CoV-2 papain-like protease with appropriate binding. Moreover, the MD simulation for 100 ns and MMPBSA calculation reveals that eugenol possess potential phytotherapeutic properties against COVID-19. The interaction of eugenol with human serum albumin has been examined by using fluorescence, UV-vis spectroscopy, circular dichroism as well as computational techniques such as molecular docking, molecular dynamic simulation and MMPBSA calculation. Overall investigation shows eugenol having good affinity for HSA Ka 6.80 × 106 M-1.Communicated by Ramaswamy H. Sarma.

The Alkylamine Stimulant 1,3-Dimethylamylamine Exhibits Substrate-Like Regulation of Dopamine Transporter Function and Localization

The alkylamine stimulant 1,3-dimethylamylamine (DMAA) is used nonmedically as an appetite suppressant and exercise performance enhancer despite adverse cardiovascular effects that have limited its legal status. There is scant research describing the mechanism of action of DMAA, making it difficult to gauge risks or therapeutic potential. An important molecular target of structurally related phenethylamines, such as amphetamine, for regulating mood, cognition, movement, and the development of substance use disorder is the dopamine transporter, which limits the range and magnitude of dopamine signaling via reuptake from the extracellular space. The present studies were therefore initiated to characterize the effects of DMAA on dopamine transporter function. Specifically, we tested the hypothesis that DMAA exhibits substrate-like effects on dopamine transporter function and trafficking. In transport assays in human embryonic kidney cells, DMAA inhibited dopamine uptake by the human dopamine transporter in a competitive manner. Docking analysis and molecular dynamics simulations supported these findings, revealing that DMAA binds to the S1 substrate binding site and induces a conformational change from outward-facing open to outward-facing closed states, similar to the known substrates. Further supporting substrate-like effects of DMAA, the drug stimulated dopamine transporter endocytosis in a heterologous expression system via cocaine- and protein kinase A-sensitive mechanisms, mirroring findings with amphetamine. Together, these data indicate that DMAA elicits neurologic effects by binding to and regulating function of the dopamine transporter. Furthermore, pharmacologic distinctions from amphetamine reveal structural determinants for regulating transporter conformation and add mechanistic insight for the regulation of dopamine transporter endocytosis. SIGNIFICANCE STATEMENT: The alkylamine stimulant 1,3-dimethylamylamine (DMAA) is used as an appetite suppressant and athletic performance enhancer and is structurally similar to amphetamine, but there is scant research describing its mechanism of action. Characterizing the effects of DMAA on dopamine transporter function supports evaluation of potential risks and therapeutic potential while also revealing mechanistic details of dynamic transporter-substrate interactions.

Assessing the effect of a series of mutations on the dynamic behavior of phosphite dehydrogenase using molecular docking, molecular dynamics and quantum mechanics/molecular mechanics simulations

Discovered in Pseudomonas stutzeri, phosphite dehydrogenase (PTDH) is an enzyme that catalyzes the oxidation of phosphite to phosphate while simultaneously reducing NAD⁺ to NADH. Despite several investigations into the mechanism of reaction and cofactor regeneration, only a few studies have focused on improving the activity and stability of PTDH. In this study, we combine molecular docking, molecular dynamics (MD) simulation, and Quantum Mechanics/Molecular Mechanics (QM/MM) to identify the impact of 30 mutations on the activity and stability of PTDH. Molecular docking results suggest that E266Q, K76A, K76M, K76R, K76C, and R237K can act on the NAD⁺ binding site through relatively weak bond development due to their high free binding energy. Moreover, Mulliken population analysis and potential energy barrier indicate that T101A, E175A, E175A/A176R, A176R, and E266Q act on phosphite oxidation. The mutants M53N, M53A, K76R, D79N, D79A, T101A, W134A, W134F Y139F, A146S, E175A, F198I, F198M, E266Q, H292K, S295A, R301K, and R301A were found to act on the structural dynamic of PTDH. The remaining mutants cause the loss of the nitrogen atom of R237 and H292, respectively, inactivating the enzyme. This study provides specific explanations of how mutations affect weak interactions of PTDH. The results should allow researchers to conduct experimental studies to improve PTDH activity and stability.Communicated by Ramaswamy H. Sarma.

The substrate import mechanism of the human serotonin transporter

The serotonin transporter (SERT) initiates the reuptake of extracellular serotonin in the synapse to terminate neurotransmission. The cryogenic electron microscopy structures of SERT bound to ibogaine and the physiological substrate serotonin resolved in different states have provided a glimpse of the functional conformations at atomistic resolution. However, the conformational dynamics and structural transitions to intermediate states are not fully understood. Furthermore, the molecular basis of how serotonin is recognized and transported remains unclear. In this study, we performed unbiased microsecond-long simulations of the human SERT to investigate the structural dynamics to various intermediate states and elucidated the complete substrate import pathway. Using Markov state models, we characterized a sequential order of conformational-driven ion-coupled substrate binding and transport events and calculated the free energy barriers of conformation transitions associated with the import mechanism. We find that the transition from the occluded to inward-facing state is the rate-limiting step for substrate import and that the substrate decreases the free energy barriers to achieve the inward-facing state. Our study provides insights on the molecular basis of dynamics-driven ion-substrate recognition and transport of SERT that can serve as a model for other closely related neurotransmitter transporters.

Probing the Impact of Temperature and Substrates on the Conformational Dynamics of the Neurotransmitter:Sodium symporter LeuT

The crucial function of neurotransmitter:sodium symporters (NSS) in facilitating the reuptake of neurotransmitters into neuronal cells makes them attractive drug targets for treating multiple mental diseases. Due to the challenges in working with eukaryotic NSS proteins, LeuT, a prokaryotic amino acid transporter, has served as a model protein for studying structure-function relationships of NSS family proteins. With hydrogen-deuterium exchange mass spectrometry (HDX-MS), slow unfolding/refolding kinetics were identified in multiple regions of LeuT, suggesting that substrate translocation involves cooperative fluctuations of helical stretches. Earlier work has solely been performed at non-native temperatures (25 °C) for LeuT, which is evolutionarily adapted to function at high temperatures (85 - 95 °C). To address the effect of temperature on LeuT dynamics, we have performed HDX-MS experiments at elevated temperatures (45 °C and 60 °C). At these elevated temperatures, multiple regions in LeuT exhibited increased dynamics compared to 25 °C. Interestingly, coordinated slow unfolding/refolding of key regions could still be observed, though considerably faster. We have further investigated the conformational impact of binding the efficiently transported substrate alanine (Ala) relative to the much slower transported substrate leucine (Leu). Comparing the HDX of the Ala-bound versus Leu-bound state of LeuT, we observe distinct differences that could explain the faster transport rate (k_cat) of Ala relative to Leu. Importantly, slow unfolding/refolding dynamics could still be observed in regions of Ala-bound LeuT . Overall, our work brings new insights into the conformational dynamics of LeuT and provides a better understanding of the transport mechanism of LeuT and possibly other transporters bearing the LeuT fold.

Functional Characterization of the Dopaminergic Psychostimulant Sydnocarb as an Allosteric Modulator of the Human Dopamine Transporter

The dopamine transporter (DAT) serves a critical role in controlling dopamine (DA)-mediated neurotransmission by regulating the clearance of DA from the synapse and extrasynaptic regions and thereby modulating DA action at postsynaptic DA receptors. Major drugs of abuse such as amphetamine and cocaine interact with DATs to alter their actions resulting in an enhancement in extracellular DA concentrations. We previously identified a novel allosteric site in the DAT and the related human serotonin transporter that lies outside the central orthosteric substrate- and cocaine-binding pocket. Here, we demonstrate that the dopaminergic psychostimulant sydnocarb is a ligand of this novel allosteric site. We identified the molecular determinants of the interaction between sydnocarb and DAT at the allosteric site using molecular dynamics simulations. Biochemical-substituted cysteine scanning accessibility experiments have supported the computational predictions by demonstrating the occurrence of specific interactions between sydnocarb and amino acids within the allosteric site. Functional dopamine uptake studies have further shown that sydnocarb is a noncompetitive inhibitor of DAT in accord with the involvement of a site different from the orthosteric site in binding this psychostimulant. Finally, DA uptake studies also demonstrate that sydnocarb affects the interaction of DAT with both cocaine and amphetamine. In summary, these studies further strengthen the prospect that allosteric modulation of DAT activity could have therapeutic potential.

Psychomotor impairments and therapeutic implications revealed by a mutation associated with infantile Parkinsonism-Dystonia

Parkinson disease (PD) is a progressive, neurodegenerative disorder affecting over 6.1 million people worldwide. Although the cause of PD remains unclear, studies of highly penetrant mutations identified in early-onset familial parkinsonism have contributed to our understanding of the molecular mechanisms underlying disease pathology. Dopamine (DA) transporter (DAT) deficiency syndrome (DTDS) is a distinct type of infantile parkinsonism-dystonia that shares key clinical features with PD, including motor deficits (progressive bradykinesia, tremor, hypomimia) and altered DA neurotransmission. Here, we define structural, functional, and behavioral consequences of a Cys substitution at R445 in human DAT (hDAT R445C), identified in a patient with DTDS. We found that this R445 substitution disrupts a phylogenetically conserved intracellular (IC) network of interactions that compromise the hDAT IC gate. This is demonstrated by both Rosetta molecular modeling and fine-grained simulations using hDAT R445C, as well as EPR analysis and X-ray crystallography of the bacterial homolog leucine transporter. Notably, the disruption of this IC network of interactions supported a channel-like intermediate of hDAT and compromised hDAT function. We demonstrate that Drosophila melanogaster expressing hDAT R445C show impaired hDAT activity, which is associated with DA dysfunction in isolated brains and with abnormal behaviors monitored at high-speed time resolution. We show that hDAT R445C Drosophila exhibit motor deficits, lack of motor coordination (i.e. flight coordination) and phenotypic heterogeneity in these behaviors that is typically associated with DTDS and PD. These behaviors are linked with altered dopaminergic signaling stemming from loss of DA neurons and decreased DA availability. We rescued flight coordination with chloroquine, a lysosomal inhibitor that enhanced DAT expression in a heterologous expression system. Together, these studies shed some light on how a DTDS-linked DAT mutation underlies DA dysfunction and, possibly, clinical phenotypes shared by DTDS and PD.

Direct coupling of oligomerization and oligomerization-driven endocytosis of the dopamine transporter to its conformational mechanics and activity

Dopamine transporter (DAT) mediates the reuptake of synaptically released dopamine, and thus controls the duration and intensity of dopamine neurotransmission. Mammalian DAT has been observed to form oligomers, although the mechanisms of oligomerization and its role in DAT activity and trafficking remain largely unknown. We discovered a series of small molecule compounds that stabilize trimers and induce high-order oligomers of DAT and concomitantly promote its clathrin-independent endocytosis. Using a combination of chemical cross-linking, fluorescence resonance energy transfer microscopy, antibody-uptake endocytosis assay, live-cell lattice light sheet microscopy, ligand binding and substrate transport kinetics analyses, and molecular modeling and simulations, we investigated molecular basis of DAT oligomerization and endocytosis induced by these compounds. Our study showed that small molecule–induced DAT oligomerization and endocytosis are favored by the inward-facing DAT conformation and involve interactions of four hydrophobic residues at the interface between transmembrane (TM) helices TM4 and TM9. Surprisingly, a corresponding quadruple DAT mutant displays altered dopamine transport kinetics and increased cocaine-analog binding. The latter is shown to originate from an increased preference for outward-facing conformation and inward-to-outward transition. Taken together, our results demonstrate a direct coupling between conformational dynamics of DAT, functional activity of the transporter, and its oligomerization leading to endocytosis. The high specificity of such coupling for DAT makes the TM4-9 hub a new target for pharmacological modulation of DAT activity and subcellular localization.

The Amino Terminus of LeuT Changes Conformation in an Environment Sensitive Manner

Neurotransmitter:sodium symporters are highly expressed in the human brain and catalyze the uptake of substrate through the plasma membrane by using the electrochemical gradient of sodium as the energy source. The bacterial homolog LeuT, a small amino acid transporter isolated from the bacteria Aquifex aeolicus, is the founding member of the family and has been crystallized in three conformations. The N-terminus is structurally well defined and strongly interacts with the transporter core in the outward-facing conformations. However, it could not be resolved in the inward-facing conformation, which indicates enhanced mobility. Here we investigate conformations and dynamics of the N-terminus, by combining molecular dynamics simulations with experimental verification using distance measurements and accessibility studies. We found strongly increased dynamics of the N-terminus, but also that helix TM1A is subject to enhanced mobility. TM1A moves towards the transporter core in the membrane environment, reaching a conformation that is closer to the structure of LeuT with wild type sequence, indicating that the mutation introduced to create the inward-facing structure might have altered the position of helix TM1A. The mobile N-terminus avoids entering the open vestibule of the inward-facing state, as accessibility studies do not show any reduction of quenching by iodide of a fluorophore attached to the N-terminus.

Mapping of Ion and Substrate Binding Sites in Human Sodium Iodide Symporter (hNIS)

The human sodium iodide symporter (hNIS) is a theranostic reporter gene which concentrates several clinically approved SPECT and PET radiotracers and plays an essential role for the synthesis of thyroid hormones as an iodide transporter in the thyroid gland. Development of hNIS mutants which could enhance translocation of the desired imaging ions is currently underway. Unfortunately, it is hindered by lack of understanding of the 3D organization of hNIS and its relation to anion transport. There are no known crystal structures of hNIS in any of its conformational states. Homology modeling can be very effective in such situations; however, the low sequence identity between hNIS and relevant secondary transporters with available experimental structures makes the choice of a template and the generation of 3D models nontrivial. Here, we report a combined application of homology modeling and molecular dynamics refining of the hNIS structure in its semioccluded state. The modeling was based on templates from the LeuT-fold protein family and was done with emphasis on the refinement of the substrate-ion binding pocket. The consensus model developed in this work is compared to available biophysical and biochemical experimental data for a number of different LeuT-fold proteins. Some functionally important residues contributing to the formation of putative binding sites and permeation pathways for the cotransported Na⁺ ions and I^- substrate were identified. The model predictions were experimentally tested by generation of mutant versions of hNIS and measurement of relative (to WT hNIS) ¹²⁵I^- uptake of 35 hNIS variants.

X-ray structure of LeuT in an inward-facing occluded conformation reveals mechanism of substrate release

Neurotransmitter:sodium symporters (NSS) are conserved from bacteria to man and serve as targets for drugs, including antidepressants and psychostimulants. Here we report the X-ray structure of the prokaryotic NSS member, LeuT, in a Na⁺/substrate-bound, inward-facing occluded conformation. To obtain this structure, we were guided by findings from single-molecule fluorescence spectroscopy and molecular dynamics simulations indicating that L-Phe binding and mutation of the conserved N-terminal Trp8 to Ala both promote an inward-facing state. Compared to the outward-facing occluded conformation, our structure reveals a major tilting of the cytoplasmic end of transmembrane segment (TM) 5, which, together with release of the N-terminus but without coupled movement of TM1, opens a wide cavity towards the second Na⁺ binding site. The structure of this key intermediate in the LeuT transport cycle, in the context of other NSS structures, leads to the proposal of an intracellular release mechanism of substrate and ions in NSS proteins.

Shared Signature Dynamics Tempered by Local Fluctuations Enables Fold Adaptability and Specificity

Recent studies have drawn attention to the evolution of protein dynamics, in addition to sequence and structure, based on the premise structure-encodes-dynamics-encodes-function. Of interest is to understand how functional differentiation is accomplished while maintaining the fold, or how intrinsic dynamics plays out in the evolution of structural variations and functional specificity. We performed a systematic computational analysis of 26,899 proteins belonging to 116 CATH superfamilies. Characterizing cooperative mechanisms and convergent/divergent features that underlie the shared/differentiated dynamics of family members required a methodology that lends itself to efficient analyses of large ensembles of proteins. We therefore introduced, SignDy, an integrated pipeline for evaluating the signature dynamics of families based on elastic network models. Our analysis confirmed that family members share conserved, highly cooperative (global) modes of motion. Importantly, our analysis discloses a subset of motions that sharply distinguishes subfamilies, which lie in a low-to-intermediate frequency regime of the mode spectrum. This regime has maximal impact on functional differentiation of families into subfamilies, while being evolutionarily conserved among subfamily members. Notably, the high-frequency end of the spectrum also reveals evolutionary conserved features across and within subfamilies; but in sharp contrast to global motions, high-frequency modes are minimally collective. Modulation of robust/conserved global dynamics by low-to-intermediate frequency fluctuations thus emerges as a versatile mechanism ensuring the adaptability of selected folds and the specificity of their subfamilies. SignDy further allows for dynamics-based categorization as a new layer of information relevant to distinctive mechanisms of action of subfamilies, beyond sequence or structural classifications.

In silico characterization of residues essential for substrate binding of human cystine transporter, xCT

xCT is a sodium-independent amino acid antiporter that imports L-cystine and exports L-glutamate in a 1:1 ratio. It is a component of heterodimeric amino acid transporter system Xc- working at the cross-roads of maintaining neurological processes and regulating antioxidant defense. The transporter has 12 transmembrane domains with intracellular N- and C-termini, and like other transporter proteins can undergo various conformational changes while switching the ligand accessibilities from intracellular to extracellular site. In the present study, we generated two homology models of human xCT in two distinct conformations: inward-facing occluded state and outward-facing open state. Our results indicated the substrate translocation channel composed of transmembrane helices TMs 1, 3, 6, 8, and 10. We docked anionic L-cystine and L-glutamate within the cavities to assess the two distinct binding scenarios for xCT as antiporter. We also assessed the interactions between the ligands and transporter and observed that ligands bind to similar residues within the channel. Using MM-PBSA/MM-GBSA approach, we computed the binding energies of these ligands to different conformational states. Cystine and glutamate bind xCT with favorable binding energies, with more favorable binding observed in inward occluded state than in outward open state. We further computed the residue-wise decomposition of these binding energies and identified the residues as essential for substrate binding/permeation. Filtering the residues that form favorable energetic contributions to the ligand binding in both the states, our studies suggest T56, A60, R135, A138, V141, Y244, A247, F250, S330, L392, and R396 as critical residues for ligand binding as well as ligand transport for any conformational state adopted by xCT during its transport cycle. .Graphical Abstract.

Advances and Challenges in Rational Drug Design for SLCs

There are over 420 human solute carrier (SLC) transporters from 65 families that are expressed ubiquitously in the body. The SLCs mediate the movement of ions, drugs, and metabolites across membranes and their dysfunction has been associated with a variety of diseases, such as diabetes, cancer, and central nervous system (CNS) disorders. Thus, SLCs are emerging as important targets for therapeutic intervention. Recent technological advances in experimental and computational biology allow better characterization of SLC pharmacology. Here we describe recent approaches to modulate SLC transporter function, with an emphasis on the use of computational approaches and computer-aided drug design (CADD) to study nutrient transporters. Finally, we discuss future perspectives in the rational design of SLC drugs.

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.

Tales of tails in transporters

Cell nutrition, detoxification, signalling, homeostasis and response to drugs, processes related to cell growth, differentiation and survival are all mediated by plasma membrane (PM) proteins called transporters. Despite their distinct fine structures, mechanism of function, energetic requirements, kinetics and substrate specificities, all transporters are characterized by a main hydrophobic body embedded in the PM as a series of tightly packed, often intertwined, α-helices that traverse the lipid bilayer in a zigzag mode, connected with intracellular or extracellular loops and hydrophilic N- and C-termini. Whereas longstanding genetic, biochemical and biophysical evidence suggests that specific transmembrane segments, and also their connecting loops, are responsible for substrate recognition and transport dynamics, emerging evidence also reveals the functional importance of transporter N- and C-termini, in respect to transport catalysis, substrate specificity, subcellular expression, stability and signalling. This review highlights selected prototypic examples of transporters in which their termini play important roles in their functioning.

Trimerization of dopamine transporter triggered by AIM-100 binding: Molecular mechanism and effect of mutations

Recent work demonstrated the propensity of dopamine transporters (DATs) to form trimers or higher oligomers, enhanced upon binding a furopyrimidine, AIM-100. AIM-100 binding promotes DAT endocytosis and thereby moderates dopaminergic transmission. Despite the neurobiological significance of these events, the molecular mechanisms that underlie the stabilization of DAT trimer and the key interactions that modulate the trimerization of DAT, and not serotonin transporter SERT, remain unclear. In the present study, we determined three structural models, termed trimer-W238, -C306 and -Y303, for possible trimerization of DATs . To this aim, we used structural data resolved for DAT and its structural homologs that share the LeuT fold, advanced computational modeling and simulations, site-directed mutagenesis experiments and live-cell imaging assays. The models are in accord with the versatility of LeuT fold to stabilize dimeric or higher order constructs. Selected residues show a high propensity to occupy interfacial regions. Among them, D231-W238 in the extracellular loop EL2, including the intersubunit salt-bridge forming pair D231/D232-R237 (not present in SERT) (in trimer-W238), the loop EL3 (trimers-C306 and -Y303), and W497 on the intracellularly exposed IL5 loop (trimer-C306) and its spatial neighbors (e.g. K525) near the C-terminus are computationally predicted and experimentally confirmed to play important roles in enabling the correct folding and/or oligomerization of DATs in the presence of AIM-100. The study suggests the possibility of controlling the effective transport of dopamine by altering the oligomerization state of DAT upon small molecule binding, as a possible intervention strategy to modulate dopaminergic signaling. This article is part of the issue entitled 'Special Issue on Neurotransmitter Transporters'.

Feedback adaptation of synaptic excitability via Glu:Na+ symport driven astrocytic GABA and Gln release

Glutamatergic transmission composed of the arriving of action potential at the axon terminal, fast vesicular Glu release, postsynaptic Glu receptor activation, astrocytic Glu clearance and Glu→Gln shuttle is an abundantly investigated phenomenon. Despite its essential role, however, much less is known about the consequences of the mechanistic connotations of Glu:Na⁺ symport. Due to the coupled Na⁺ transport, Glu uptake results in significantly elevated intracellular astrocytic [Na⁺] that markedly alters the driving force of other Na⁺-coupled astrocytic transporters. The resulting GABA and Gln release by reverse transport through the respective GAT-3 and SNAT3 transporters help to re-establish the physiological Na⁺ homeostasis without ATP dissipation and consequently leads to enhanced tonic inhibition and replenishment of axonal glutamate pool. Here, we place this emerging astrocytic adjustment of synaptic excitability into the centre of future perspectives. This article is part of the issue entitled 'Special Issue on Neurotransmitter Transporters'.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Shared dynamics of LeuT superfamily members and allosteric differentiation by structural irregularities and multimerization

The LeuT-fold superfamily includes secondary active transporters from different functional families, which share a common tertiary structure, despite having a remarkably low sequence similarity. By identifying the common structural and dynamical features upon principal component analysis of a comprehensive ensemble of 90 experimentally resolved structures and anisotropic network model evaluation of collective motions, we provide a unified point of view for understanding the reasons why this particular fold has been selected by evolution to accomplish such a broad spectrum of functions. The parallel identification of conserved sequence features, localized at specific sites of transmembrane helices, sheds light on the role of broken helices (TM1 and TM6 in LeuT) in promoting ion/substrate binding and allosteric interconversion between the outward- and inward-facing conformations of transporters. Finally, the determination of the dynamics landscape for the structural ensemble provides a promising framework for the classification of transporters based on their dynamics, and the characterization of the collective movements that favour multimerization.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

In silico Description of LAT1 Transport Mechanism at an Atomistic Level

The molecular mechanism of transport mediated by LAT1, a sodium-independent antiporter of large neutral amino acids, was investigated through in silico procedures, specifically making reference to two transported substrates, tyrosine (Tyr) and leucine methyl ester (LME), and to 3,5-diiodo-L-tyrosine (DIT), a well-known LAT1 inhibitor. Two models of the transporter were built by comparative modeling, with LAT1 either in an outward-facing (OF) or in an inward-facing (IF) conformation, based, respectively, on the crystal structure of AdiC and of GadC. As frequently classic Molecular Dynamics (MD) fails to monitor large-scale conformational transitions within a reasonable simulated time, the OF structure was equilibrated for 150 ns then processed through targeted MD (tMD). During this procedure, an elastic force pulled the OF structure to the IF structure and induced, at the same time, substrates/inhibitor to move through the transport channel. This elastic force was modulated by a spring constant (k) value; by decreasing its value from 100 to 70, it was possible to comparatively account for the propensity for transport of the three tested molecules. In line with our expectations, during the tMD simulations, Tyr and LME behaved as substrates, moving down the transport channel, or most of it, for all k values. On the contrary, DIT behaved as an inhibitor, being (almost) transported across the channel only at the highest k value (100). During their transit through the channel, Tyr and LME interacted with specific amino acids (first with Phe252 then with Thr345, Arg348, Tyr259, and Phe262); this suggests that a primary as well as a putative secondary gate may contribute to the transport of substrates. Quite on the opposite, DIT appeared to establish only transient interactions with side chains lining the external part of the transport channel. Our tMD simulations could thus efficiently discriminate between two transported substrates and one inhibitor, and therefore can be proposed as a benchmark for developing novel LAT1 inhibitors of pharmacological interest.

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.

A partially-open inward-facing intermediate conformation of LeuT is associated with Na+ release and substrate transport

Neurotransmitter:sodium symporters (NSS), targets of antidepressants and psychostimulants, clear neurotransmitters from the synaptic cleft through sodium (Na+)-coupled transport. Substrate and Na+ are thought to be transported from the extracellular to intracellular space through an alternating access mechanism by coordinated conformational rearrangements in the symporter that alternately expose the binding sites to each side of the membrane. However, the mechanism by which the binding of ligands coordinates conformational changes occurring on opposite sides of the membrane is not well understood. Here, we report the use of single-molecule fluorescence resonance energy transfer (smFRET) techniques to image transitions between distinct conformational states on both the extracellular and intracellular sides of the prokaryotic NSS LeuT, including partially open intermediates associated with transport activity. The nature and functional context of these hitherto unidentified intermediate states shed new light on the allosteric mechanism that couples substrate and Na+ symport by the NSS family through conformational dynamics.

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.

Latch and trigger role for R445 in DAT transport explains molecular basis of DTDS

A recent study reports on five different mutations as sources of dopamine transporter (DAT) deficiency syndrome (DTDS). One of these mutations, R445C, is believed to be located on the intracellular side of DAT distal to the primary (S1) or secondary (S2) sites to which substrate binding is understood to occur. Thus, the molecular mechanism by which the R445C mutation results in DAT transport deficiency has eluded explanation. However, the recently reported X-ray structures of the endogenous amine transporters for dDAT and hSERT revealed the presence of a putative salt bridge between R445 and E428 suggesting a possible mechanism. To evaluate whether the R445C effect is a result of a salt bridge interaction, the mutants R445E, E428R, and the double mutant E428R/R445E were generated. The single mutants R445E and E428R displayed loss of binding and transport properties of the substrate [³H]DA and inhibitor [³H]CFT at the cell surface while the double mutant E428R/R445E, although nonfunctional, restored [³H]DA and [³H]CFT binding affinity to that of WT. Structure based analyses of these results led to a model wherein R445 plays a dual role in normal DAT function. R445 acts as a component of a latch in its formation of a salt bridge with E428 which holds the primary substrate binding site (S1) in place and helps enforce the inward closed protein state. When this salt bridge is broken, R445 acts as a trigger which disrupts a local polar network and leads to the release of the N-terminus from its position inducing the inward closed state to one allowing the inward open state. In this manner, both the loss of binding and transport properties of the R445C variant are explained.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Targeting of dopamine transporter to filopodia requires an outward-facing conformation of the transporter

Dopamine transporter (DAT) has been shown to accumulate in filopodia in neurons and non-neuronal cells. To examine the mechanisms of DAT filopodial targeting, we used quantitative live-cell fluorescence microscopy, and compared the effects of the DAT inhibitor cocaine and its fluorescent analog JHC1-64 on the plasma membrane distribution of wild-type DAT and two non-functional DAT mutants, R60A and W63A, that do not accumulate in filopodia. W63A did not bind JHC1-64, whereas R60A did, although less efficiently compared to the wild-type DAT. Molecular dynamics simulations predicted that R60A preferentially assumes an outward-facing (OF) conformation through compensatory intracellular salt bridge formation, which in turn favors binding of cocaine. Imaging analysis showed that JHC1-64-bound R60A mutant predominantly localized in filopodia, whereas free R60A molecules were evenly distributed within the plasma membrane. Cocaine binding significantly increased the density of R60A, but not that of W63A, in filopodia. Further, zinc binding, known to stabilize the OF state, also increased R60A concentration in filopodia. Finally, amphetamine, that is thought to disrupt DAT OF conformation, reduced the concentration of wild-type DAT in filopodia. Altogether, these data indicate that OF conformation is required for the efficient targeting of DAT to, and accumulation in, filopodia.

Mechanism of Substrate Translocation in an Alternating Access Transporter

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

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.

Substrate transport and anion permeation proceed through distinct pathways in glutamate transporters

Advances in structure-function analyses and computational biology have enabled a deeper understanding of how excitatory amino acid transporters (EAATs) mediate chloride permeation and substrate transport. However, the mechanism of structural coupling between these functions remains to be established. Using a combination of molecular modeling, substituted cysteine accessibility, electrophysiology and glutamate uptake assays, we identified a chloride-channeling conformer, iChS, transiently accessible as EAAT1 reconfigures from substrate/ion-loaded into a substrate-releasing conformer. Opening of the anion permeation path in this iChS is controlled by the elevator-like movement of the substrate-binding core, along with its wall that simultaneously lines the anion permeation path (global); and repacking of a cluster of hydrophobic residues near the extracellular vestibule (local). Moreover, our results demonstrate that stabilization of iChS by chemical modifications favors anion channeling at the expense of substrate transport, suggesting a mutually exclusive regulation mediated by the movement of the flexible wall lining the two regions.

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.

Conformational dynamics of a neurotransmitter:sodium symporter in a lipid bilayer

Neurotransmitter:sodium symporters (NSSs) are integral membrane proteins responsible for the sodium-dependent reuptake of small-molecule neurotransmitters from the synaptic cleft. The symporters for the biogenic amines serotonin (SERT), dopamine (DAT), and norepinephrine (NET) are targets of multiple psychoactive agents, and their dysfunction has been implicated in numerous neuropsychiatric ailments. LeuT, a thermostable eubacterial NSS homolog, has been exploited as a model protein for NSS members to canvass the conformational mechanism of transport with a combination of X-ray crystallography, cysteine accessibility, and solution spectroscopy. Despite yielding remarkable insights, these studies have primarily been conducted with protein in the detergent-solubilized state rather than embedded in a membrane mimic. In addition, solution spectroscopy has required site-specific labeling of nonnative cysteines, a labor-intensive process occasionally resulting in diminished transport and/or binding activity. Here, we overcome these limitations by reconstituting unlabeled LeuT in phospholipid bilayer nanodiscs, subjecting them to hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS), and facilitating interpretation of the data with molecular dynamics simulations. The data point to changes of accessibility and dynamics of structural elements previously implicated in the transport mechanism, in particular transmembrane helices (TMs) 1a and 7 as well as extracellular loops (ELs) 2 and 4. The results therefore illuminate the value of this strategy for interrogating the conformational mechanism of the more clinically significant mammalian membrane proteins including SERT and DAT, neither of which tolerates complete removal of endogenous cysteines, and whose activity is heavily influenced by neighboring lipids.

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.

A Markov State-based Quantitative Kinetic Model of Sodium Release from the Dopamine Transporter

The dopamine transporter (DAT) belongs to the neurotransmitter:sodium symporter (NSS) family of membrane proteins that are responsible for reuptake of neurotransmitters from the synaptic cleft to terminate a neuronal signal and enable subsequent neurotransmitter release from the presynaptic neuron. The release of one sodium ion from the crystallographically determined sodium binding site Na2 had been identified as an initial step in the transport cycle which prepares the transporter for substrate translocation by stabilizing an inward-open conformation. We have constructed Markov State Models (MSMs) from extensive molecular dynamics simulations of human DAT (hDAT) to explore the mechanism of this sodium release. Our results quantify the release process triggered by hydration of the Na2 site that occurs concomitantly with a conformational transition from an outward-facing to an inward-facing state of the transporter. The kinetics of the release process are computed from the MSM, and transition path theory is used to identify the most probable sodium release pathways. An intermediate state is discovered on the sodium release pathway, and the results reveal the importance of various modes of interaction of the N-terminus of hDAT in controlling the pathways of release.

What Can and Cannot Be Learned from Molecular Dynamics Simulations of Bacterial Proton-Coupled Oligopeptide Transporter GkPOT?

We have performed an extensive set of all-atom molecular dynamics (MD) simulations of a bacterial proton-coupled oligopeptide transporter (POT) in an explicit membrane environment. We have characterized both the local and global conformational dynamics of the transporter upon the proton and/or substrate binding, within a statistical framework. Our results reveal a clearly distinct behavior for local conformational dynamics in the absence and presence of the proton at the putative proton binding residue E310. Particularly, we find that the substrate binding conformation is drastically different in the two conditions, where the substrate binds to the protein in a lateral/vertical manner, in the presence/absence of the proton. We do not observe any statistically significant distinctive behavior in terms of the global conformational changes in different simulation conditions, within the time scales of our simulations. Our extensive simulations and analyses call into question the implicit assumption of many MD studies that local conformational changes observed in short simulations could provide clues to the global conformational changes that occur on much longer time scales. The linear regression analysis of quantities associated with the global conformational fluctuations, however, provides an indication of a mechanism involving the concerted motion of the transmembrane helices, consistent with the rocker-switch mechanism.

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.

Unveiling the Mechanism of Arginine Transport through AdiC with Molecular Dynamics Simulations: The Guiding Role of Aromatic Residues

Commensal and pathogenic enteric bacteria have developed several systems to adapt to proton leakage into the cytoplasm resulting from extreme acidic conditions. One such system involves arginine uptake followed by export of the decarboxylated product agmatine, carried out by the arginine/agmatine antiporter (AdiC), which thus works as a virtual proton pump. Here, using classical and targeted molecular dynamics, we investigated at the atomic level the mechanism of arginine transport through AdiC of E. coli. Overall, our MD simulation data clearly demonstrate that global rearrangements of several transmembrane segments are necessary but not sufficient for achieving transitions between structural states along the arginine translocation pathway. In particular, local structural changes, namely rotameric conversions of two aromatic residues, are needed to regulate access to both the outward- and inward-facing states. Our simulations have also enabled identification of a few residues, overwhelmingly aromatic, which are essential to guiding arginine in the course of its translocation. Most of them belong to gating elements whose coordinated motions contribute to the alternating access mechanism. Their conservation in all known E. coli acid resistance antiporters suggests that the transport mechanisms of these systems share common features. Last but not least, knowledge of the functional properties of AdiC can advance our understanding of the members of the amino acid-carbocation-polyamine superfamily, notably in eukaryotic cells.

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.

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.

Allosteric Mechanisms of Molecular Machines at the Membrane: Transport by Sodium-Coupled Symporters

Solute transport across cell membranes is ubiquitous in biology as an essential physiological process. Secondary active transporters couple the unfavorable process of solute transport against its concentration gradient to the energetically favorable transport of one or several ions. The study of such transporters over several decades indicates that their function involves complex allosteric mechanisms that are progressively being revealed in atomistic detail. We focus on two well-characterized sodium-coupled symporters: the bacterial amino acid transporter LeuT, which is the prototype for the "gated pore" mechanism in the mammalian synaptic monoamine transporters, and the archaeal GltPh, which is the prototype for the "elevator" mechanism in the mammalian excitatory amino acid transporters. We present the evidence for the role of allostery in the context of a quantitative formalism that can reconcile biochemical and biophysical data and thereby connects directly to recent insights into the molecular structure and dynamics of these proteins. We demonstrate that, while the structures and mechanisms of these transporters are very different, the available data suggest a common role of specific models of allostery in their functions. We argue that such allosteric mechanisms appear essential not only for sodium-coupled symport in general but also for the function of other types of molecular machines in the membrane.

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.

Atomic-level characterization of transport cycle thermodynamics in the glycerol-3-phosphate:phosphate antiporter

Membrane transporters actively translocate their substrate by undergoing large-scale structural transitions between inward- (IF) and outward-facing (OF) states ('alternating-access' mechanism). Despite extensive structural studies, atomic-level mechanistic details of such structural transitions, and as importantly, their coupling to chemical events supplying the energy, remain amongst the most elusive aspects of the function of these proteins. Here we present a quantitative, atomic-level description of the functional thermodynamic cycle for the glycerol-3-phosphate:phosphate antiporter GlpT by using a novel approach in reconstructing the free energy landscape governing the IF↔OF transition along a cyclic transition pathway involving both apo and substrate-bound states. Our results provide a fully atomic description of the complete transport process, offering a structural model for the alternating-access mechanism and substantiating the close coupling between global structural transitions and local chemical events.