Intracellular gating in an inward-facing state of aspartate transporter Glt(Ph) is regulated by the movements of the helical hairpin HP2

Conformational free energy landscape of a glutamate transporter and microscopic details of its transport mechanism

Removing glutamate from the synaptic cleft is vital for proper function of the brain. Excitatory amino acid transporters mediate this process by uptaking the neurotransmitter from the synaptic cleft back to the cell after its release. The archaeal homolog, GltPh, an aspartate transporter from Pyrococcus horikoshii, presents the best structurally characterized model for this family of transporters. In order to transport, GltPh undergoes elevator-like conformational changes between inward-facing (IF) and outward-facing (OF) states. Here, we characterize, at an atomic level, the OF⇌IF transition of GltPh in different apo/bound states using a combination of ensemble-based enhanced sampling techniques, employing more than two thousand of coupled simulation replicas of membrane-embedded GltPh. The resulting free-energy profiles portray the transition of apo/bound states as a complex four-stage process, while sodium binding alone locks the structure in one of its states. Along the transition, the transport domain (TD) disengages from the scaffold domain (SD), allowing it to move as a piston sliding vertically with respect to the membrane during the elevator-like motion of TD. Lipid interactions with residues comprising the SD-TD interface directly influence the large-scale conformational changes and, consequently, the energetics of transport. Structural intermediates formed during the transition leak water molecules and may correlate to the uncoupled Cl- ion conductance observed experimentally in both prokaryotic and mammalian glutamate transporters. Mechanistic insights obtained from our study provide a structural framework for better development of therapeutic for neurological disorders.

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.

Linking function to global and local dynamics in an elevator-type transporter

Transporters cycle through large structural changes to translocate molecules across biological membranes. The temporal relationships between these changes and function, and the molecular properties setting their rates, determine transport efficiency-yet remain mostly unknown. Using single-molecule fluorescence microscopy, we compare the timing of conformational transitions and substrate uptake in the elevator-type transporter GltPh We show that the elevator-like movements of the substrate-loaded transport domain across membranes and substrate release are kinetically heterogeneous, with rates varying by orders of magnitude between individual molecules. Mutations increasing the frequency of elevator transitions and reducing substrate affinity diminish transport rate heterogeneities and boost transport efficiency. Hydrogen deuterium exchange coupled to mass spectrometry reveals destabilization of secondary structure around the substrate-binding site, suggesting that increased local dynamics leads to faster rates of global conformational changes and confers gain-of-function properties that set transport rates.

Mechanism and potential sites of potassium interaction with glutamate transporters

In the mammalian glutamate transporters, countertransported intracellular K⁺ is essential for relocating the glutamate binding site to the extracellular side of the membrane. This K⁺-dependent process is believed to be rate limiting for the transport cycle. In contrast, extracellular K⁺ induces glutamate release upon transporter reversal. Here, we analyzed potential K⁺ binding sites using molecular dynamics (MD) simulations and site-directed mutagenesis. Two candidate sites were identified by spontaneous K⁺ binding in MD simulations, one site (K1 site) overlapping with the Na1 Na⁺ binding site and the K2 site being localized under hairpin loop 2 (HP2). Mutations to conserved amino acid residues in these sites resulted in several transporters that were defective in K⁺-induced reverse transport and which bound K⁺ with reduced apparent affinity compared with the wild-type transporter. However, external K⁺ interaction was abolished in only one mutant transporter EAAC1_D454A in the K1 site. Our results, for the first time, directly demonstrate effects of K1-site mutations on K⁺ binding, in contrast to previous reports on K⁺ binding sites based on indirect evidence. We propose that K⁺ binding to the K1 site is responsible for catalyzing the relocation step, whereas binding to the K2 site may have an as-of-yet unidentified regulatory function.

ASCT1 and ASCT2: Brother and Sister?

The SLC1 family includes seven members divided into two groups, namely, EAATs and ASCTs, that share similar 3D architecture; the first one includes high-affinity glutamate transporters, and the second one includes SLC1A4 and SLC1A5, known as ASCT1 and ASCT2, respectively, responsible for the traffic of neutral amino acids across the cell plasma membrane. The physiological role of ASCT1 and ASCT2 has been investigated over the years, revealing different properties in terms of substrate specificities, affinities, and regulation by physiological effectors and posttranslational modifications. Furthermore, ASCT1 and ASCT2 are involved in pathological conditions, such as neurodegenerative disorders and cancer. This has driven research in the pharmaceutical field aimed to find drugs able to target the two proteins.This review focuses on structural, functional, and regulatory aspects of ASCT1 and ASCT2, highlighting similarities and differences.

Observing spontaneous, accelerated substrate binding in molecular dynamics simulations of glutamate transporters

Glutamate transporters are essential for removing the neurotransmitter glutamate from the synaptic cleft. Glutamate transport across the membrane is associated with elevator-like structural changes of the transport domain. These structural changes require initial binding of the organic substrate to the transporter. Studying the binding pathway of ligands to their protein binding sites using molecular dynamics (MD) simulations requires micro-second level simulation times. Here, we used three methods to accelerate aspartate binding to the glutamate transporter homologue Gltph and to investigate the binding pathway. 1) Two methods using user-defined forces to prevent the substrate from diffusing too far from the binding site. 2) Conventional MD simulations using very high substrate concentrations in the 0.1 M range. The final, substrate bound states from these methods are comparable to the binding pose observed in crystallographic studies, although they show more flexibility in the side chain carboxylate function. We also captured an intermediate on the binding pathway, where conserved residues D390 and D394 stabilize the aspartate molecule. Finally, we investigated glutamate binding to the mammalian glutamate transporter, excitatory amino acid transporter 1 (EAAT1), for which a crystal structure is known, but not in the glutamate-bound state. Overall, the results obtained in this study reveal new insights into the pathway of substrate binding to glutamate transporters, highlighting intermediates on the binding pathway and flexible conformational states of the side chain, which most likely become locked in once the hairpin loop 2 closes to occlude the substrate.

On the Role of a Conserved Methionine in the Na+-Coupling Mechanism of a Neurotransmitter Transporter Homolog

Excitatory amino acid transporters (EAAT) play a key role in glutamatergic synaptic communication. Driven by transmembrane cation gradients, these transporters catalyze the reuptake of glutamate from the synaptic cleft once this neurotransmitter has been utilized for signaling. Two decades ago, pioneering studies in the Kanner lab identified a conserved methionine within the transmembrane domain as key for substrate turnover rate and specificity; later structural work, particularly for the prokaryotic homologs Glt_Ph and Glt_Tk, revealed that this methionine is involved in the coordination of one of the three Na⁺ ions that are co-transported with the substrate. Albeit extremely atypical, the existence of this interaction is consistent with biophysical analyses of Glt_Ph showing that mutations of this methionine diminish the binding cooperativity between substrates and Na⁺. It has been unclear, however, whether this intriguing methionine influences the thermodynamics of the transport reaction, i.e., its substrate:ion stoichiometry, or whether it simply fosters a specific kinetics in the binding reaction, which, while influential for the turnover rate, do not fundamentally explain the ion-coupling mechanism of this class of transporters. Here, studies of Glt_Tk using experimental and computational methods independently arrive at the conclusion that the latter hypothesis is the most plausible, and lay the groundwork for future efforts to uncover the underlying mechanism.

The high-energy transition state of the glutamate transporter homologue GltPh

Membrane transporters mediate cellular uptake of nutrients, signaling molecules, and drugs. Their overall mechanisms are often well understood, but the structural features setting their rates are mostly unknown. Earlier single-molecule fluorescence imaging of the archaeal model glutamate transporter homologue GltPh from Pyrococcus horikoshii suggested that the slow conformational transition from the outward- to the inward-facing state, when the bound substrate is translocated from the extracellular to the cytoplasmic side of the membrane, is rate limiting to transport. Here, we provide insight into the structure of the high-energy transition state of GltPh that limits the rate of the substrate translocation process. Using bioinformatics, we identified GltPh gain-of-function mutations in the flexible helical hairpin domain HP2 and applied linear free energy relationship analysis to infer that the transition state structurally resembles the inward-facing conformation. Based on these analyses, we propose an approach to search for allosteric modulators for transporters.

Na+-dependent gate dynamics and electrostatic attraction ensure substrate coupling in glutamate transporters

Excitatory amino acid transporters (EAATs) harness [Na⁺], [K⁺], and [H⁺] gradients for fast and efficient glutamate removal from the synaptic cleft. Since each glutamate is cotransported with three Na⁺ ions, [Na⁺] gradients are the predominant driving force for glutamate uptake. We combined all-atom molecular dynamics simulations, fluorescence spectroscopy, and x-ray crystallography to study Na⁺:substrate coupling in the EAAT homolog Glt_Ph A lipidic cubic phase x-ray crystal structure of wild-type, Na⁺-only bound Glt_Ph at 2.5-Å resolution revealed the fully open, outward-facing state primed for subsequent substrate binding. Simulations and kinetic experiments established that only the binding of two Na⁺ ions to the Na1 and Na3 sites ensures complete HP2 gate opening via a conformational selection-like mechanism and enables high-affinity substrate binding via electrostatic attraction. The combination of Na⁺-stabilized gate opening and electrostatic coupling of aspartate to Na⁺ binding provides a constant Na⁺:substrate transport stoichiometry over a broad range of neurotransmitter concentrations.

Millisecond dynamics of an unlabeled amino acid transporter

Excitatory amino acid transporters (EAATs) are important in many physiological processes and crucial for the removal of excitatory amino acids from the synaptic cleft. Here, we develop and apply high-speed atomic force microscopy line-scanning (HS-AFM-LS) combined with automated state assignment and transition analysis for the determination of transport dynamics of unlabeled membrane-reconstituted Glt_Ph, a prokaryotic EAAT homologue, with millisecond temporal resolution. We find that Glt_Ph transporters can operate much faster than previously reported, with state dwell-times in the 50 ms range, and report the kinetics of an intermediate transport state with height between the outward- and inward-facing states. Transport domains stochastically probe transmembrane motion, and reversible unsuccessful excursions to the intermediate state occur. The presented approach and analysis methodology are generally applicable to study transporter kinetics at system-relevant temporal resolution.

Excitatory amino acid transporters (EAATs) are crucial for the removal of excitatory amino acids from the synaptic cleft. Here authors combined high-speed atomic force microscopy line-scanning with automated state assignment for the determination of transport dynamics of Glt_Ph, a prokaryotic EAAT homologue, with millisecond temporal resolution.

Allosteric gate modulation confers K+ coupling in glutamate transporters

Excitatory amino acid transporters (EAATs) mediate glial and neuronal glutamate uptake to terminate synaptic transmission and to ensure low resting glutamate concentrations. Effective glutamate uptake is achieved by cotransport with 3 Na⁺ and 1 H⁺, in exchange with 1 K⁺. The underlying principles of this complex transport stoichiometry remain poorly understood. We use molecular dynamics simulations and electrophysiological experiments to elucidate how mammalian EAATs harness K⁺ gradients, unlike their K⁺‐independent prokaryotic homologues. Glutamate transport is achieved via elevator‐like translocation of the transport domain. In EAATs, glutamate‐free re‐translocation is prevented by an external gate remaining open until K⁺ binding closes and locks the gate. Prokaryotic Glt_Ph contains the same K⁺‐binding site, but the gate can close without K⁺. Our study provides a comprehensive description of K⁺‐dependent glutamate transport and reveals a hitherto unknown allosteric coupling mechanism that permits adaptions of the transport stoichiometry without affecting ion or substrate binding.

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.

Glutamate transporters: a broad review of the most recent archaeal and human structures

Glutamate transporters play important roles in bacteria, archaea and eukaryotes. Their function in the mammalian central nervous system is essential for preventing excitotoxicity, and their dysregulation is implicated in many diseases, such as epilepsy and Alzheimer's. Elucidating their transport mechanism would further the understanding of these transporters and promote drug design as they provide compelling targets for understanding the pathophysiology of diseases and may have a direct role in the treatment of conditions involving glutamate excitotoxicity. This review outlines the insights into the transport cycle, uncoupled chloride conductance and modulation, as well as identifying areas that require further investigation.

Investigation of the allosteric coupling mechanism in a glutamate transporter homolog via unnatural amino acid mutagenesis

Glutamate transporters harness the ionic gradients across cell membranes for the concentrative uptake of glutamate. The sodium-coupled Asp symporter, Glt_Ph is an archaeal homolog of glutamate transporters and has been extensively used to understand the transport mechanism. A critical aspect of the transport cycle in Glt_Ph is the coupled binding of sodium and aspartate. Previous studies have suggested a major role for hairpin-2 (HP2), which functions as the extracellular gate for the aspartate binding site, in the coupled binding of sodium and aspartate to Glt_Ph In this study, we develop a fluorescence assay for monitoring HP2 movement by incorporating tryptophan and the unnatural amino acid, p-cyanophenylalanine into Glt_Ph We use the HP2 assays to show that HP2 opening with Na⁺ follows an induced-fit mechanism. We also determine how residues in the substrate binding site affect the opening and closing of HP2. Our data, combined with previous studies, provide the molecular sequence of events in the coupled binding of sodium and aspartate to Glt_Ph.

Free-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPh

Glutamate transporters clear up excess extracellular glutamate by cotransporting three Na⁺ and one H⁺ with the countertransport of one K⁺. The archaeal homologs are selective to aspartate and only cotransport three Na⁺. The crystal structures of Glt_Ph from archaea have been used in computational studies to understand the transport mechanism. Although some progress has been made with regard to the ligand-binding sites, a consistent picture of transport still eludes us. A major concern is the discrepancy between the computed binding free energies, which predict high-affinity Na⁺-low-affinity aspartate binding, and the experimental results in which the opposite is observed. Here, we show that the binding of the first two Na⁺ ions involves an intermediate state near the Na1 site, where two Na⁺ ions coexist and couple to aspartate with similar strengths, boosting its affinity. Binding free energies for Na⁺ and aspartate obtained using this intermediate state are in good agreement with the experimental values. Thus, the paradox in binding affinities arises from the assumption that the ligands bind to the sites observed in the crystal structure following the order dictated by their binding free energies with no intermediate states. In fact, the presence of an intermediate state eliminates such a correlation between the binding free energies and the binding order. The intermediate state also facilitates transition of the first Na⁺ ion to its final binding site via a knock-on mechanism, which induces substantial conformational changes in the protein consistent with experimental observations.

Thermodynamic first law efficiency of membrane proteins

Proteins are nature's biomolecular machines. Proteins, such as transporters, pumps and motors, have complex function/operating-machinery/mechanisms, comparable to the macro-scaled machines that we encounter in our daily life. These proteins, as it is for their macro-scaled counterparts, convert (part of) other/various forms of energy into work. In this study, we are performing the first law analysis on a set of proteins, including the dopamine transporter, glycine transporters I and II, glutamate transporter, sodium-potassium pump and Ca²⁺ ATPase. Each of these proteins operates on a thermodynamic/mechanic cycle to perform their function. In each of these cycles, they receive energy from a source, convert part of this energy into work and reject the remaining part of the energy to the environment. Conservation of energy principle was applied to the thermodynamic/mechanic cycle of each protein, and thermodynamic first law efficiency was evaluated for each cycle, which shows how much of the energy input per cycle was converted into useful work. Interestingly, calculations based on experimental data indicate that proteins can operate under a range of efficiencies, which vary based on the extracellular and intracellular ion and substrate concentrations. The lowest observed first law efficiency was 50%, which is a very high value if compared to the efficiency of the macro-scaled heat engines we encounter in our daily lives.Communicated by Ramaswamy H. Sarma.

Kinetic mechanism of coupled binding in sodium-aspartate symporter GltPh

Many secondary active membrane transporters pump substrates against concentration gradients by coupling their uptake to symport of sodium ions. Symport requires the substrate and ions to be always transported together. Cooperative binding of the solutes is a key mechanism contributing to coupled transport in the sodium and aspartate symporter from Pyrococcus horikoshii Glt_Ph. Here, we describe the kinetic mechanism of coupled binding for Glt_Ph in the inward facing state. The first of the three coupled sodium ions, binds weakly and slowly, enabling the protein to accept the rest of the ions and the substrate. The last ion binds tightly, but is in rapid equilibrium with solution. Its release is required for the complex disassembly. Thus, the first ion serves to ‘open the door’ for the substrate, the last ion ‘locks the door’ once the substrate is in, and one ion contributes to both events.

Transport mechanism of a glutamate transporter homologue GltPh

Glutamate transporters are responsible for uptake of the neurotransmitter glutamate in mammalian central nervous systems. Their archaeal homologue GltPh, an aspartate transporter isolated from Pyrococcus horikoshii, has been the focus of extensive studies through crystallography, MD simulations and single-molecule FRET (smFRET). Here, we summarize the recent research progress on GltPh, in the hope of gaining some insights into the transport mechanism of this aspartate transporter.

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.

Structural insights into the elevator-like mechanism of the sodium/citrate symporter CitS

The sodium-dependent citrate transporter of Klebsiella pneumoniae (KpCitS) belongs to the 2-hydroxycarboxylate transporter (2-HCT) family and allows the cell to use citrate as sole carbon and energy source in anaerobic conditions. Here we present crystal structures of KpCitS in citrate-bound outward-facing, citrate-bound asymmetric, and citrate-free inward-facing state. The structures reveal that the KpCitS dimerization domain remains stationary throughout the transport cycle due to a hydrogen bond network as well as extensive hydrophobic interactions. In contrast, its transport domain undergoes a ~35° rigid-body rotation and a ~17 Å translocation perpendicular to the membrane to expose the substrate-binding site alternately to either side of the membrane. Furthermore, homology models of two other 2-HCT proteins based on the KpCitS structure offer structural insights into their differences in substrate specificity at a molecular level. On the basis of our results and previous biochemical data, we propose that the activity of the 2-HCT CitS involves an elevator-like movement in which the transport domain itself traverses the lipid bilayer, carrying the substrate into the cell in a sodium-dependent manner.

Characterization of the Inward- and Outward-Facing Substrate Binding Sites of the Prokaryotic Aspartate Transporter, GltPh

Crystal structures of the prokaryotic aspartate transporter, Glt_Ph, have provided important insights into the mechanism of amino acid transport by Glt_Ph and related eukaryotic members of the glutamate transporter family (SLC1A family). Identification of inhibitors of Glt_Ph can provide valuable tools for understanding the molecular basis for substrate and inhibitor specificity and selectivity of SLC1A members, but at present, few inhibitors of Glt_Ph have been identified. We have screened a collection of commercially available aspartate analogues and identified new transportable and nontransportable Glt_Ph inhibitors. We have explored the inhibition profile of Glt_Ph by utilizing a thiol modification assay that isolates sided populations of the transporters reconstituted in liposomes to determine if any aspartate analogues display a preference for either the inwardly or outwardly directed binding sites. Here, we have characterized several new inhibitors of Glt_Ph and identified three β-carbon-substituted molecules that display a strong preference for the outwardly directed binding site of Glt_Ph.

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.

Mechanism of Na(+)-dependent citrate transport from the structure of an asymmetrical CitS dimer

The common human pathogen Salmonella enterica takes up citrate as a nutrient via the sodium symporter SeCitS. Uniquely, our 2.5 Å x-ray structure of the SeCitS dimer shows three different conformations of the active protomer. One protomer is in the outside-facing state. Two are in different inside-facing states. All three states resolve the substrates in their respective binding environments. Together with comprehensive functional studies on reconstituted proteoliposomes, the structures explain the transport mechanism in detail. Our results indicate a six-step process, with a rigid-body 31° rotation of a helix bundle that translocates the bound substrates by 16 Å across the membrane. Similar transport mechanisms may apply to a wide variety of related and unrelated secondary transporters, including important drug targets.

DOI:http://dx.doi.org/10.7554/eLife.09375.001

Cells have specialized proteins known as transporters in their surface membranes that move molecules into or out of the cell. Transporters pass their cargo through the membrane by changing shape. This process requires energy and is sometimes driven by simultaneously transporting a charged ion such as sodium. There are different classes of transporters and researchers have described a range of structural changes, and therefore transport mechanisms, that different transporters use.

Citrate transporters are found in a wide range of organisms. In bacteria, they bring the citrate substrate molecule into the cell to be used as a nutrient. In humans, citrate transporters are important in metabolism, and are of interest as targets for drugs that could potentially treat obesity and diabetes. This requires an understanding of the atomic structure and the transport mechanisms used by citrate transporters, which were not known.

Wöhlert et al. now use a technique called X-ray crystallography to uncover the structure of a citrate transporter called SeCitS in high detail. This transporter is found in a bacterium called Salmonella enterica, a well-known human pathogen that causes typhoid. The crystallized protein simultaneously showed three different structural states – one where the citrate binding site faces the outside of the cell, and two where the binding site faces the inside of the cell. The simultaneous occurrence of different functional states in one and the same crystal structure of a membrane transporter is so far unique.

Combining the detailed structures of SeCitS with biochemical studies allowed Wöhlert et al. to deduce that citrate is transported in a six-step process. Sodium ions attach to SeCitS, and then citrate binds to the transporter from outside the cell. This binding causes part of the protein to undergo a substantial rotation, shifting it to an inward-facing state and moving the citrate and sodium ions inside the cell. The release of the citrate and sodium ions then triggers the reverse rotation of the transporter, bringing the empty binding site back to the outside of the cell for a repeat of the cycle.

After working out the mechanisms of a bacterial citrate transporter, the next challenge is to extend the analysis to the structure of similar transporters in more complex organisms, including human cells. This could provide an accurate basis for drug development.

DOI:http://dx.doi.org/10.7554/eLife.09375.002

Computational Studies of Glutamate Transporters

Glutamate is the major excitatory neurotransmitter in the human brain whose binding to receptors on neurons excites them while excess glutamate are removed from synapses via transporter proteins. Determination of the crystal structures of bacterial aspartate transporters has paved the way for computational investigation of their function and dynamics at the molecular level. Here, we review molecular dynamics and free energy calculation methods used in these computational studies and discuss the recent applications to glutamate transporters. The focus of the review is on the insights gained on the transport mechanism through computational methods, which otherwise is not directly accessible by experimental probes. Recent efforts to model the mammalian glutamate and other amino acid transporters, whose crystal structures have not been solved yet, are included in the review.

Refinement of the Central Steps of Substrate Transport by the Aspartate Transporter GltPh: Elucidating the Role of the Na2 Sodium Binding Site

Glutamate homeostasis in the brain is maintained by glutamate transporter mediated accumulation. Impaired transport is associated with several neurological disorders, including stroke and amyotrophic lateral sclerosis. Crystal structures of the homolog transporter GltPh from Pyrococcus horikoshii revealed large structural changes. Substrate uptake at the atomic level and the mechanism of ion gradient conversion into directional transport remained enigmatic. We observed in repeated simulations that two local structural changes regulated transport. The first change led to formation of the transient Na2 sodium binding site, triggered by side chain rotation of T308. The second change destabilized cytoplasmic ionic interactions. We found that sodium binding to the transiently formed Na2 site energized substrate uptake through reshaping of the energy hypersurface. Uptake experiments in reconstituted proteoliposomes confirmed the proposed mechanism. We reproduced the results in the human glutamate transporter EAAT3 indicating a conserved mechanics from archaea to humans.

We used the archaeal homolog GltPh of the human glutamate transporters to refine our understanding how large scale conformational changes are translated into substrate translocation. We identified the structural changes that accompany substrate transport and convert the energy stored in the ion gradient into a directional transport. Insights into the mechanics of these transporters are likely to increase our understanding of how they contribute to excitotoxicity and to develop drugs, which preclude the underlying accumulation of glutamate in the synaptic cleft.

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.

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.

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

Neurotransmitter: sodium symporters (NSSs) regulate neuronal signal transmission by clearing excess neurotransmitters from the synapse, assisted by the co-transport of sodium ions. Extensive structural data have been collected in recent years for several members of the NSS family, which opened the way to structure-based studies for a mechanistic understanding of substrate transport. Leucine transporter (LeuT), a bacterial orthologue, has been broadly adopted as a prototype in these studies. This goal has been elusive, however, due to the complex interplay of global and local events as well as missing structural data on LeuT N-terminal segment. We provide here for the first time a comprehensive description of the molecular events leading to substrate/Na+ release to the postsynaptic cell, including the structure and dynamics of the N-terminal segment using a combination of molecular simulations. Substrate and Na+-release follows an influx of water molecules into the substrate/Na+-binding pocket accompanied by concerted rearrangements of transmembrane helices. A redistribution of salt bridges and cation-π interactions at the N-terminal segment prompts substrate release. Significantly, substrate release is followed by the closure of the intracellular gate and a global reconfiguration back to outward-facing state to resume the transport cycle. Two minimally hydrated intermediates, not structurally resolved to date, are identified: one, substrate-bound, stabilized during the passage from outward- to inward-facing state (holo-occluded), and another, substrate-free, along the reverse transition (apo-occluded).

Coupled ion binding and structural transitions along the transport cycle of glutamate transporters

Membrane transporters that clear the neurotransmitter glutamate from synapses are driven by symport of sodium ions and counter-transport of a potassium ion. Previous crystal structures of a homologous archaeal sodium and aspartate symporter showed that a dedicated transport domain carries the substrate and ions across the membrane. Here, we report new crystal structures of this homologue in ligand-free and ions-only bound outward- and inward-facing conformations. We show that after ligand release, the apo transport domain adopts a compact and occluded conformation that can traverse the membrane, completing the transport cycle. Sodium binding primes the transport domain to accept its substrate and triggers extracellular gate opening, which prevents inward domain translocation until substrate binding takes place. Furthermore, we describe a new cation-binding site ideally suited to bind a counter-transported ion. We suggest that potassium binding at this site stabilizes the translocation-competent conformation of the unloaded transport domain in mammalian homologues.

DOI:http://dx.doi.org/10.7554/eLife.02283.001

Molecules of glutamate can carry messages between cells in the brain, and these signals are essential for thought and memory. Glutamate molecules can also act as signals to build new connections between brain cells and to prune away unnecessary ones. However, too much glutamate outside of the cells kills the brain tissue and can lead to devastating brain diseases.

In a healthy brain, special pumps called glutamate transporters move these molecules back into the brain cells, where they can be stored safely. However, when brain cells are damaged—by, for example, a stroke or an injury,—the glutamate stored inside spills out, killing the surrounding cells. This leads to a cascade of dying cells and leaking glutamate, which causes even more damage and slows the recovery.

Glutamate transporters ensure that there are more glutamate molecules inside cells than outside. However, it requires energy to maintain this gradient in the concentration of glutamate molecules. The transporters get this energy by moving three sodium ions into the cell with each glutamate molecule, and moving one potassium ion out of the cell. However, it is not clear how these transporters ensure that they move the glutamate molecules and the sodium ions at the same time.

Now, Verdon, Oh et al. have uncovered the 3D structure of a glutamate transporter homologue at each step of the transport process. These structures reveal that, on the outside of the cell membrane, sodium ions attach to the so-called ‘transporter domain’ and make it better able to bind glutamate. The transporter domain then carries the sodium ions and glutamate through the cell membrane and releases them into the cell. Verdon, Oh et al. suggest that a potassium ion then binds to the empty transport domain, stabilizing it into a more compact shape that easily makes the return trip to the outside of the cell.

Most experiments on glutamate transporters, including the work of Verdon, Oh et al., are carried out on model proteins taken from bacteria. An important challenge for the future will be to obtain structural information on human glutamate transporters, as these could be therapeutic targets for the treatment of various neurological conditions.

DOI:http://dx.doi.org/10.7554/eLife.02283.002

Global motions exhibited by proteins in micro- to milliseconds simulations concur with anisotropic network model predictions

The Anton supercomputing technology recently developed for efficient molecular dynamics simulations permits us to examine micro- to milli-second events at full atomic resolution for proteins in explicit water and lipid bilayer. It also permits us to investigate to what extent the collective motions predicted by network models (that have found broad use in molecular biophysics) agree with those exhibited by full-atomic long simulations. The present study focuses on Anton trajectories generated for two systems: the bovine pancreatic trypsin inhibitor, and an archaeal aspartate transporter, GltPh. The former, a thoroughly studied system, helps benchmark the method of comparative analysis, and the latter provides new insights into the mechanism of function of glutamate transporters. The principal modes of motion derived from both simulations closely overlap with those predicted for each system by the anisotropic network model (ANM). Notably, the ANM modes define the collective mechanisms, or the pathways on conformational energy landscape, that underlie the passage between the crystal structure and substates visited in simulations. In particular, the lowest frequency ANM modes facilitate the conversion between the most probable substates, lending support to the view that easy access to functional substates is a robust determinant of evolutionarily selected native contact topology.

Exploring the conformational transitions of biomolecular systems using a simple two-state anisotropic network model

Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results.

Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.

Molecular dynamics simulations of the mammalian glutamate transporter EAAT3

Excitatory amino acid transporters (EAATs) are membrane proteins that enable sodium-coupled uptake of glutamate and other amino acids into neurons. Crystal structures of the archaeal homolog GltPh have been recently determined both in the inward- and outward-facing conformations. Here we construct homology models for the mammalian glutamate transporter EAAT3 in both conformations and perform molecular dynamics simulations to investigate its similarities and differences from GltPh. In particular, we study the coordination of the different ligands, the gating mechanism and the location of the proton and potassium binding sites in EAAT3. We show that the protonation of the E374 residue is essential for binding of glutamate to EAAT3, otherwise glutamate becomes unstable in the binding site. The gating mechanism in the inward-facing state of EAAT3 is found to be different from that of GltPh, which is traced to the relocation of an arginine residue from the HP1 segment in GltPh to the TM8 segment in EAAT3. Finally, we perform free energy calculations to locate the potassium binding site in EAAT3, and find a high-affinity site that overlaps with the Na1 and Na3 sites in GltPh.