Ion-releasing state of a secondary membrane transporter

Sodium-Dependent Neutral Amino Acid Transporter 2 Can Serve as a Tertiary Carrier for l-Type Amino Acid Transporter 1-Utilizing Prodrugs

Membrane transporters are the key determinants of the homeostasis of endogenous compounds in the cells and their exposure to drugs. However, the substrate specificities of distinct transporters can overlap. In the present study, the interactions of l-type amino acid transporter 1 (LAT1)-utilizing prodrugs with sodium-coupled neutral amino acid transporter 2 (SNAT2) were explored. The results showed that the cellular uptake of LAT1-utilizing prodrugs into a human breast cancer cell line, MCF-7 cells, was mediated via SNATs as the uptake was increased at higher pH (8.5), decreased in the absence of sodium, and inhibited in the presence of unselective SNAT-inhibitor, (α-(methylamino)isobutyric acid, MeAIB). Moreover, docking the compounds to a SNAT2 homology model (inward-open conformation) and further molecular dynamics simulations and the subsequent trajectory and principal component analyses confirmed the chemical features supporting the interactions of the studied compounds with SNAT2, which was found to be the main SNAT expressed in MCF-7 cells.

Cryo-EM structure of the human NKCC1 transporter reveals mechanisms of ion coupling and specificity

The sodium-potassium-chloride transporter NKCC1 of the SLC12 family performs Na⁺ -dependent Cl^- - and K⁺ -ion uptake across plasma membranes. NKCC1 is important for regulating cell volume, hearing, blood pressure, and regulation of hyperpolarizing GABAergic and glycinergic signaling in the central nervous system. Here, we present a 2.6 Å resolution cryo-electron microscopy structure of human NKCC1 in the substrate-loaded (Na⁺ , K⁺ , and 2 Cl^- ) and occluded, inward-facing state that has also been observed for the SLC6-type transporters MhsT and LeuT. Cl^- binding at the Cl1 site together with the nearby K⁺ ion provides a crucial bridge between the LeuT-fold scaffold and bundle domains. Cl^- -ion binding at the Cl2 site seems to undertake a structural role similar to conserved glutamate of SLC6 transporters and may allow for Cl^- -sensitive regulation of transport. Supported by functional studies in mammalian cells and computational simulations, we describe a putative Na⁺ release pathway along transmembrane helix 5 coupled to the Cl2 site. The results provide insight into the structure-function relationship of NKCC1 with broader implications for other SLC12 family members.

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.

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

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

Electrogenic Binding of Intracellular Cations Defines a Kinetic Decision Point in the Transport Cycle of the Human Serotonin Transporter

The plasmalemmal monoamine transporters clear the extracellular space from their cognate substrates and sustain cellular monoamine stores even during neuronal activity. In some instances, however, the transporters enter a substrate-exchange mode, which results in release of intracellular substrate. Understanding what determines the switch between these two transport modes demands time-resolved measurements of intracellular (co-)substrate binding and release. Here, we report an electrophysiological investigation of intracellular solute-binding to the human serotonin transporter (SERT) expressed in HEK-293 cells. We measured currents induced by rapid application of serotonin employing varying intracellular (co-)substrate concentrations and interpreted the data using kinetic modeling. Our measurements revealed that the induction of the substrate-exchange mode depends on both voltage and intracellular Na⁺ concentrations because intracellular Na⁺ release occurs before serotonin release and is highly electrogenic. This voltage dependence was blunted by electrogenic binding of intracellular K⁺ and, notably, also H⁺. In addition, our data suggest that Cl⁻ is bound to SERT during the entire catalytic cycle. Our experiments, therefore, document an essential role of electrogenic binding of K⁺ or of H⁺ to the inward-facing conformation of SERT in (i) cancelling out the electrogenic nature of intracellular Na⁺ release and (ii) in selecting the forward-transport over the substrate-exchange mode. Finally, the kinetics of intracellular Na⁺ release and K⁺ (or H⁺) binding result in a voltage-independent rate-limiting step where SERT may return to the outward-facing state in a KCl- or HCl-bound form.

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

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

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

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

Stochastic steps in secondary active sugar transport

Secondary active transporters, such as those that adopt the leucine-transporter fold, are found in all domains of life, and they have the unique capability of harnessing the energy stored in ion gradients to accumulate small molecules essential for life as well as expel toxic and harmful compounds. How these proteins couple ion binding and transport to the concomitant flow of substrates is a fundamental structural and biophysical question that is beginning to be answered at the atomistic level with the advent of high-resolution structures of transporters in different structural states. Nonetheless, the dynamic character of the transporters, such as ion/substrate binding order and how binding triggers conformational change, is not revealed from static structures, yet it is critical to understanding their function. Here, we report a series of molecular simulations carried out on the sugar transporter vSGLT that lend insight into how substrate and ions are released from the inward-facing state of the transporter. Our simulations reveal that the order of release is stochastic. Functional experiments were designed to test this prediction on the human homolog, hSGLT1, and we also found that cytoplasmic release is not ordered, but we confirmed that substrate and ion binding from the extracellular space is ordered. Our findings unify conflicting published results concerning cytoplasmic release of ions and substrate and hint at the possibility that other transporters in the superfamily may lack coordination between ions and substrate in the inward-facing state.

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

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

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

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

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.

Conformational Changes in Two Inter-Helical Loops of Mhp1 Membrane Transporter

Mhp1 is a bacterial secondary transporter with high-resolution crystal structures available for both the outward- and inward-facing conformations. Through molecular dynamics simulations of the ligand-free Mhp1 as well as analysis of its crystal structures, here we show that two inter-helical loops, respectively located at the extra- and intracellular ends of the "hash motif" in the protein, play important roles in the conformational transition. In the outward- and inward-facing states of the protein, the loops adopt different secondary structures, either wrapped to the end of an alpha-helix, or unwrapped to extended conformations. In equilibrium simulations of 100 ns with Mhp1 in explicit lipids and water, the loop conformations remain largely stable. In targeted molecular dynamics simulations with the protein structure driven from one state to the other, the loops exhibit resistance and only undergo abrupt changes when other parts of the protein already approach the target conformation. Free energy calculations on the isolated loops further confirm that the wrapping/unwrapping transitions are associated with substantial energetic barriers, and consist of multiple sequential steps involving the rotation of certain backbone torsion angles. Furthermore, in simulations with the loops driven from one state to the other, a large part of the protein follows the loops to the target conformation. Taken together, our simulations suggest that changes of the loop secondary structures would be among the slow degrees of freedom in the conformational transition of the entire protein. Incorporation of detailed loop structures into the reaction coordinate, therefore, should improve the convergence and relevance of the resulting conformational free energy.

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

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

Metadynamics simulations reveal a Na+ independent exiting path of galactose for the inward-facing conformation of vSGLT

Sodium-Galactose Transporter (SGLT) is a secondary active symporter which accumulates sugars into cells by using the electrochemical gradient of Na⁺ across the membrane. Previous computational studies provided insights into the release process of the two ligands (galactose and sodium ion) into the cytoplasm from the inward-facing conformation of Vibrio parahaemolyticus sodium/galactose transporter (vSGLT). Several aspects of the transport mechanism of this symporter remain to be clarified: (i) a detailed kinetic and thermodynamic characterization of the exit path of the two ligands is still lacking; (ii) contradictory conclusions have been drawn concerning the gating role of Y263; (iii) the role of Na⁺ in modulating the release path of galactose is not clear. In this work, we use bias-exchange metadynamics simulations to characterize the free energy profile of the galactose and Na⁺ release processes toward the intracellular side. Surprisingly, we find that the exit of Na⁺ and galactose is non-concerted as the cooperativity between the two ligands is associated to a transition that is not rate limiting. The dissociation barriers are of the order of 11–12 kcal/mol for both the ion and the substrate, in line with kinetic information concerning this type of transporters. On the basis of these results we propose a branched six-state alternating access mechanism, which may be shared also by other members of the LeuT-fold transporters.

Membrane proteins are crucial for the communication of the cell with the environment. Among these, symporters are in charge of the transport of molecules (like sugars, amino acids, osmolytes) inside the cells, exploiting the concentration gradient of an ion to perform the task. Here we investigate by atomistic simulations the transport mechanism of the Sodium-Galactose symporter. Our results allow constructing a detailed and quantitative model of the release process of the two ligands. Surprisingly, we find that the galactose is released to the cytosol independently from the ion, unambiguously indicating that the coupling in their transport mechanism is associated to the steps preceding the release process. A large family of symporters shares the same fold and potentially the same transport mechanism. As such our results are important also because they can provide insights on common mechanistic features of these transporters.

Identification of a second substrate-binding site in solute-sodium symporters

The structure of the sodium/galactose transporter (vSGLT), a solute-sodium symporter (SSS) from Vibrio parahaemolyticus, shares a common structural fold with LeuT of the neurotransmitter-sodium symporter family. Structural alignments between LeuT and vSGLT reveal that the crystallographically identified galactose-binding site in vSGLT is located in a more extracellular location relative to the central substrate-binding site (S1) in LeuT. Our computational analyses suggest the existence of an additional galactose-binding site in vSGLT that aligns to the S1 site of LeuT. Radiolabeled galactose saturation binding experiments indicate that, like LeuT, vSGLT can simultaneously bind two substrate molecules under equilibrium conditions. Mutating key residues in the individual substrate-binding sites reduced the molar substrate-to-protein binding stoichiometry to ~1. In addition, the related and more experimentally tractable SSS member PutP (the Na(+)/proline transporter) also exhibits a binding stoichiometry of 2. Targeting residues in the proposed sites with mutations results in the reduction of the binding stoichiometry and is accompanied by severely impaired translocation of proline. Our data suggest that substrate transport by SSS members requires both substrate-binding sites, thereby implying that SSSs and neurotransmitter-sodium symporters share common mechanistic elements in substrate transport.

A mechanism for intracellular release of Na+ by neurotransmitter/sodium symporters

Neurotransmitter/sodium symporters (NSSs) terminate synaptic signal transmission by Na+-dependent reuptake of released neurotransmitters. Key conformational states have been reported for the bacterial homolog LeuT and an inhibitor-bound Drosophila dopamine transporter. However, a coherent mechanism of Na+-driven transport has not been described. Here, we present two crystal structures of MhsT, an NSS member from Bacillus halodurans, in occluded inward-facing states with bound Na+ ions and L-tryptophan, providing insight into the cytoplasmic release of Na+. The switch from outward- to inward-oriented states is centered on the partial unwinding of transmembrane helix 5, facilitated by a conserved GlyX9Pro motif that opens an intracellular pathway for water to access the Na2 site. We propose a mechanism, based on our structural and functional findings, in which solvation through the TM5 pathway facilitates Na+ release from Na2 and the transition to an inward-open state.

Substrate-bound outward-open state of the betaine transporter BetP provides insights into Na+ coupling

The Na(+)-coupled betaine symporter BetP shares a highly conserved fold with other sequence unrelated secondary transporters, for example, with neurotransmitter symporters. Recently, we obtained atomic structures of BetP in distinct conformational states, which elucidated parts of its alternating-access mechanism. Here, we report a structure of BetP in a new outward-open state in complex with an anomalous scattering substrate, adding a fundamental piece to an unprecedented set of structural snapshots for a secondary transporter. In combination with molecular dynamics simulations these structural data highlight important features of the sequential formation of the substrate and sodium-binding sites, in which coordinating water molecules play a crucial role. We observe a strictly interdependent binding of betaine and sodium ions during the coupling process. All three sites undergo progressive reshaping and dehydration during the alternating-access cycle, with the most optimal coordination of all substrates found in the closed state.

Functional identification and characterization of sodium binding sites in Na symporters

Sodium cotransporters from several different gene families belong to the leucine transporter (LeuT) structural family. Although the identification of Na(+) in binding sites is beyond the resolution of the structures, two Na(+) binding sites (Na1 and Na2) have been proposed in LeuT. Na2 is conserved in the LeuT family but Na1 is not. A biophysical method has been used to measure sodium dissociation constants (Kd) of wild-type and mutant human sodium glucose cotransport (hSGLT1) proteins to identify the Na(+) binding sites in hSGLT1. The Na1 site is formed by residues in the sugar binding pocket, and their mutation influences sodium binding to Na1 but not to Na2. For the canonical Na2 site formed by two -OH side chains, S392 and S393, and three backbone carbonyls, mutation of S392 to cysteine increased the sodium Kd by sixfold. This was accompanied by a dramatic reduction in the apparent sugar and phlorizin affinities. We suggest that mutation of S392 in the Na2 site produces a structural rearrangement of the sugar binding pocket to disrupt both the binding of the second Na(+) and the binding of sugar. In contrast, the S393 mutations produce no significant changes in sodium, sugar, and phlorizin affinities. We conclude that the Na2 site is conserved in hSGLT1, the side chain of S392 and the backbone carbonyl of S393 are important in the first Na(+) binding, and that Na(+) binding to Na2 promotes binding to Na1 and also sugar binding.

Arginine oscillation explains Na+ independence in the substrate/product antiporter CaiT

Most secondary-active transporters transport their substrates using an electrochemical ion gradient. In contrast, the carnitine transporter (CaiT) is an ion-independent, l-carnitine/γ-butyrobetaine antiporter belonging to the betaine/carnitine/choline transporter family of secondary transporters. Recently determined crystal structures of CaiT from Escherichia coli and Proteus mirabilis revealed an inverted five-transmembrane-helix repeat similar to that in the amino acid/Na(+) symporter LeuT. The ion independence of CaiT makes it unique in this family. Here we show that mutations of arginine 262 (R262) make CaiT Na(+)-dependent. The transport activity of R262 mutants increased by 30-40% in the presence of a membrane potential, indicating substrate/Na(+) cotransport. Structural and biochemical characterization revealed that R262 plays a crucial role in substrate binding by stabilizing the partly unwound TM1' helix. Modeling CaiT from P. mirabilis in the outward-open and closed states on the corresponding structures of the related symporter BetP reveals alternating orientations of the buried R262 sidechain, which mimic sodium binding and unbinding in the Na(+)-coupled substrate symporters. We propose that a similar mechanism is operative in other Na(+)/H(+)-independent transporters, in which a positively charged amino acid replaces the cotransported cation. The oscillation of the R262 sidechain in CaiT indicates how a positive charge triggers the change between outward-open and inward-open conformations as a unifying critical step in LeuT-type transporters.

The molecular mechanism of ion-dependent gating in secondary transporters

LeuT-like fold Na-dependent secondary active transporters form a large family of integral membrane proteins that transport various substrates against their concentration gradient across lipid membranes, using the free energy stored in the downhill concentration gradient of sodium ions. These transporters play an active role in synaptic transmission, the delivery of key nutrients, and the maintenance of osmotic pressure inside the cell. It is generally believed that binding of an ion and/or a substrate drives the conformational dynamics of the transporter. However, the exact mechanism for converting ion binding into useful work has yet to be established. Using a multi-dimensional path sampling (string-method) followed by all-atom free energy simulations, we established the principal thermodynamic and kinetic components governing the ion-dependent conformational dynamics of a LeuT-like fold transporter, the sodium/benzyl-hydantoin symporter Mhp1, for an entire conformational cycle. We found that inward-facing and outward-facing states of Mhp1 display nearly the same free energies with an ion absent from the Na2 site conserved across the LeuT-like fold transporters. The barrier separating an apo-state from inward-facing or outward-facing states of the transporter is very low, suggesting stochastic gating in the absence of ion/substrate bound. In contrast, the binding of a Na2 ion shifts the free energy stabilizing the outward-facing state and promoting substrate binding. Our results indicate that ion binding to the Na2 site may also play a key role in the intracellular thin gate dynamics modulation by altering its interactions with the transmembrane helix 5 (TM5). The Potential of Mean Force (PMF) computations for a substrate entrance displays two energy minima that correspond to the locations of the main binding site S1 and proposed allosteric S2 binding site. However, it was found that substrate's binds to the site S1 ∼5 kcal/mol more favorable than that to the site S2 for all studied bound combinations of ions and a substrate.

Insights on Na(+) binding and conformational dynamics in multidrug and toxic compound extrusion transporter NorM

MATE (multidrug and toxic compound extrusion) transporter proteins mediate metabolite transport in plants and multidrug resistance in bacteria and mammals. MATE transporter NorM from Vibrio cholerae is an antiporter that is driven by Na+ gradient to extrude the substrates. To understand the molecular mechanism of Na+-substrate exchange, molecular dynamics simulation was performed to study conformational changes of both wild-type and mutant NorM with and without cation bindings. Our results show that NorM is able to bind two Na(+) ions simultaneously, one to each of the carboxylic groups of E255 and D371 in the binding pocket. Furthermore, this di-Na(+) binding state is likely more efficient for conformational changes of NorM_VC toward the inward-facing conformation than single-Na(+) binding state. The observation of two Na(+) binding sites of NorM_VC is consistent with the previous study that two sites for ion binding (denoted as Na1/Na2 sites) are found in the transporter LeuT and BetP, another two secondary transporters. Taken together, our findings shed light on the structure rearrangements of NorM on Na(+) binding and enrich our knowledge of the transport mechanism of secondary transporters.

Visualizing functional motions of membrane transporters with molecular dynamics simulations

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

A Candidate Ion-Retaining State in the Inward-Facing Conformation of Sodium/Galactose Symporter: Clues from Atomistic Simulations

The recent Vibrio parahaemolyticus sodium/galactose (vSGLT) symporter crystal structure captures the protein in an inward-facing substrate-bound conformation, with the sodium ion placed, by structural alignment, in a site equivalent to the Na2 site of the leucine transporter (LeuT). A recent study, based on molecular dynamics simulations, showed that the sodium ion spontaneously leaves its initial position diffusing outside vSGLT, toward the intracellular space. This suggested that the crystal structure corresponds to an ion-releasing state of the transporter. Here, using metadynamics, we identified a more stable Na(+) binding site corresponding to a putative ion-retaining state of the transporter. In addition, our simulations, consistently with mutagenesis studies, highlight the importance of D189 that, without being one of the Na(+)-coordinating residues, regulates its binding/release.

Simulation studies of the mechanism of membrane transporters

Membrane transporters facilitate active transport of their specific substrates, often against their electrochemical gradients across the membrane, through coupling the process to various sources of cellular energy, for example, ATP binding and hydrolysis in primary transporters, and pre-established electrochemical gradient of molecular species other than the substrate in the case of secondary transporters. In order to provide efficient energy-coupling mechanisms, membrane transporters have evolved into molecular machines in which stepwise binding, translocation, and transformation of various molecular species are closely coupled to protein conformational changes that take the transporter from one functional state to another during the transport cycle. Furthermore, in order to prevent the formation of leaky states and to be able to pump the substrate against its electrochemical gradient, all membrane transporters use the widely-accepted "alternating access mechanism," which ensures that the substrate is only accessible from one side of the membrane at a given time, but relies on complex and usually global protein conformational changes that differ for each family of membrane transporters. Describing the protein conformational changes of different natures and magnitudes is therefore at the heart of mechanistic studies of membrane transporters. Here, using a number of membrane transporters from diverse families, we present common protocols used in setting up and performing molecular dynamics simulations of membrane transporters and in analyzing the results, in order to characterize relevant motions of the system. The emphasis will be on highlighting how optimal design of molecular dynamics simulations combined with mechanistically oriented analysis can shed light onto key functionally relevant protein conformational changes in this family of membrane proteins.

Investigation of the sodium-binding sites in the sodium-coupled betaine transporter BetP

Sodium-coupled substrate transport plays a central role in many biological processes. However, despite knowledge of the structures of several sodium-coupled transporters, the location of the sodium-binding site(s) often remains unclear. Several of these structures have the five transmembrane-helix inverted-topology repeat, LeuT-like (FIRL) fold, whose pseudosymmetry has been proposed to facilitate the alternating-access mechanism required for transport. Here, we provide biophysical, biochemical, and computational evidence for the location of the two cation-binding sites in the sodium-coupled betaine symporter BetP. A recent X-ray structure of BetP in a sodium-bound closed state revealed that one of these sites, equivalent to the Na2 site in related transporters, is located between transmembrane helices 1 and 8 of the FIRL-fold; here, we confirm the location of this site by other means. Based on the pseudosymmetry of this fold, we hypothesized that the second site is located between the equivalent helices 6 and 3. Molecular dynamics simulations of the closed-state structure suggest this second sodium site involves two threonine sidechains and a backbone carbonyl from helix 3, a phenylalanine from helix 6, and a water molecule. Mutating the residues proposed to form the two binding sites increased the apparent K(m) and K(d) for sodium, as measured by betaine uptake, tryptophan fluorescence, and (22)Na(+) binding, and also diminished the transient currents measured in proteoliposomes using solid supported membrane-based electrophysiology. Taken together, these results provide strong evidence for the identity of the residues forming the sodium-binding sites in BetP.

A conformational switch in a partially unwound helix selectively determines the pathway for substrate release from the carnitine/γ-butyrobetaine antiporter CaiT

CaiT is a homotrimeric antiporter that exchanges l-carnitine (CRN) with γ-butyrobetaine (GBB) across the bacterial membrane. Three structures have been resolved to date for CaiT, all in the inward-facing state: CRN-bound (with four CRNs per subunit), GBB-bound (two GBBs per subunit), and apo. One of the reported binding sites is the counterpart of the primary site observed in structurally similar transporters. However, the mechanism and pathway(s) of CRN/GBB unbinding and translocation, or even the ability of the substrates to dislodge from the reported binding sites, are yet to be determined. To shed light on these issues, we performed a total of 1.3 μs of molecular dynamics simulations and examined the dynamics of substrate-bound CaiT structures under different conditions. We find that both CRN and GBB are able to dissociate completely from their primary site into the cytoplasm. Substrate molecules initially located at the secondary sites dissociate even faster (within tens of nanoseconds) into the extra- or intracellular regions. Interestingly, the unbinding pathway from the primary site appears to be dictated by the geometry of the unwound part of the transmembrane (TM) helix 3, mostly around Thr¹⁰⁰ therein. Arg²⁶² on TM7, which apparently mimics the role of Na⁺ in CaiT structural homologues, plays a key role in triggering the dissociation of the substrate away from the primary site and guiding its release to the cytoplasm provided that the unwound part of TM3 switches from a shielding to a yielding pose.

Carbohydrate force fields

Carbohydrates present a special set of challenges to the generation of force fields. First, the tertiary structures of monosaccharides are complex merely by virtue of their exceptionally high number of chiral centers. In addition, their electronic characteristics lead to molecular geometries and electrostatic landscapes that can be challenging to predict and model. The monosaccharide units can also interconnect in many ways, resulting in a large number of possible oligosaccharides and polysaccharides, both linear and branched. These larger structures contain a number of rotatable bonds, meaning they potentially sample an enormous conformational space. This article briefly reviews the history of carbohydrate force fields, examining and comparing their challenges, forms, philosophies, and development strategies. Then it presents a survey of recent uses of these force fields, noting trends, strengths, deficiencies, and possible directions for future expansion.

A comparative study of structures and structural transitions of secondary transporters with the LeuT fold

Secondary active transporters from several protein families share a core of two five-helix inverted repeats that has become known as the LeuT fold. The known high-resolution protein structures with this fold were analyzed by structural superposition of the core transmembrane domains (TMDs). Three angle parameters derived from the mean TMD axes correlate with accessibility of the central binding site from the outside or inside. Structural transitions between distinct conformations were analyzed for four proteins in terms of changes in relative TMD arrangement and in internal conformation of TMDs. Collectively moving groups of TMDs were found to be correlated in the covariance matrix of elastic network models. The main features of the structural transitions can be reproduced with the 5 % slowest normal modes of anisotropic elastic network models. These results support the rocking bundle model for the major conformational change between the outward- and inward-facing states of the protein and point to an important role for the independently moving last TMDs of each repeat in occluding access to the central binding site. Occlusion is also supported by flexing of some individual TMDs in the collectively moving bundle and hash motifs.

The structural pathway for water permeation through sodium-glucose cotransporters

Although water permeation across cell membranes occurs through several types of membrane proteins, the only permeation mechanism resolved at atomic scale is that through aquaporins. Crystallization of the Vibrio parahaemolyticus sodium-galactose transporter (vSGLT) allows investigation of putative water permeation pathways through both vSGLT and the homologous human Na-glucose cotransporter (hSGLT1) using computational methods. Grand canonical Monte Carlo and molecular dynamics simulations were used to stably insert water molecules in both proteins, showing the presence of a water-filled pathway composed of ∼100 water molecules. This provides a structural basis for passive water permeation that is difficult to reconcile with the water cotransport hypothesis. Potential-of-mean-force calculations of water going through the crystal structure of vSGLT shows a single barrier of 7.7 kCal mol(-1), in agreement with previously published experimental data for cotransporters of the SGLT family. Electrophysiological and volumetric experiments performed on hSGLT1-expressing Xenopus oocytes showed that the passive permeation pathway exists in different conformational states. In particular, experimental conditions that aim to mimic the conformation of the crystal structure displayed passive water permeability. These results provide groundwork for understanding the structural basis of cotransporter water permeability.

A gate-free pathway for substrate release from the inward-facing state of the Na⁺-galactose transporter

Employing molecular dynamics (MD) simulations, the pathway and mechanism of substrate unbinding from the inward-facing state of the Na(+)-coupled galactose transporter, vSGLT, have been investigated. During a 200-ns equilibrium simulation, repeated spontaneous unbinding events of the substrate from its binding site have been observed. In contrast to the previously proposed gating role of a tyrosine residue (Y263), the unbinding mechanism captured in the present equilibrium simulation does not rely on the displacement and/or rotation of this side chain. Rather, the unbinding involves an initial lateral displacement of the substrate out of the binding site which allows the substrate to completely emerge from the region covered by the side chain of Y263 without any noticeable conformational changes of the latter. Starting with the snapshots taken from this equilibrium simulation with the substrate outside the binding site, steered MD (SMD) simulations were then used to probe the translocation of the substrate along the remaining of the release pathway within the protein's lumen and to characterize the nature of protein-substrate interactions involved in the process. Combining the results of the equilibrium and SMD simulations, we provide a description of the full translocation pathway for the substrate release from the binding site into the cytoplasm. Residues E68, N142, T431, and N267 facilitate the initial substrate's displacement out of the binding site, while the translocation of the substrate along the remainder of the exit pathway formed between TM6 and TM8 is facilitated by H-bond interactions between the substrate and a series of conserved, polar residues (Y138, N267, R273, S365, S368, N371, S372, and T375). The observed molecular events indicate that no gating is required for the release of the substrate from the crystallographically captured structure of the inward-facing state of SGLT, suggesting that this conformation might represent an open, rather than occluded, state of the transporter. This article is part of a Special Issue entitled: Membrane protein structure and function.

Biology of human sodium glucose transporters

There are two classes of glucose transporters involved in glucose homeostasis in the body, the facilitated transporters or uniporters (GLUTs) and the active transporters or symporters (SGLTs). The energy for active glucose transport is provided by the sodium gradient across the cell membrane, the Na(+) glucose cotransport hypothesis first proposed in 1960 by Crane. Since the cloning of SGLT1 in 1987, there have been advances in the genetics, molecular biology, biochemistry, biophysics, and structure of SGLTs. There are 12 members of the human SGLT (SLC5) gene family, including cotransporters for sugars, anions, vitamins, and short-chain fatty acids. Here we give a personal review of these advances. The SGLTs belong to a structural class of membrane proteins from unrelated gene families of antiporters and Na(+) and H(+) symporters. This class shares a common atomic architecture and a common transport mechanism. SGLTs also function as water and urea channels, glucose sensors, and coupled-water and urea transporters. We also discuss the physiology and pathophysiology of SGLTs, e.g., glucose galactose malabsorption and familial renal glycosuria, and briefly report on targeting of SGLTs for new therapies for diabetes.

Conserved tyrosine in the first transmembrane segment of solute:sodium symporters is involved in Na+-coupled substrate co-transport

Solute:sodium symporters (SSSs) transport vital molecules across the plasma membrane of all living organisms. vSGLT, the Na(+)/galactose transporter of Vibrio parahemeolyticus, is the only SSS for which high resolution structural information is available, revealing a LeuT-like fold and a Na(+)-binding site analogous to the Na2 site of LeuT. Whereas the core transmembrane segments (TMs) of SSSs share high structural similarity with other transporters of LeuT-like fold, TM1 does not correspond to any TM in those structural homologs and was only resolved for the backbone atoms in the initial vSGLT structure (Protein Data Bank code 3DH4). To assess the role of TM1 in Na(+)-coupled substrate symport by the SSSs, here we have studied the role of a conserved residue in TM1 by computational modeling in conjunction with radiotracer transport and binding studies. Based on our sequence alignment and much topological data for homologous PutP, the Na(+)/proline transporter, we have simulated a series of vSGLT models with shifted TM1 residue assignments. We show that in two converged vSGLT models that retained the original TM1 backbone conformation, a conserved residue, Tyr-19, is associated with the Na(+) binding interaction network. In silico and in vitro mutagenesis of homologous Tyr-14 in PutP revealed the involvement of this conserved residue in Na(+)-dependent substrate binding and transport. Thus, our combined computational and experimental data provide the first clues about the importance of a conserved residue in TM1, a unique TM in the proteins with LeuT-like fold, in the Na(+)-coupled symport mechanism of SSSs.

Role of Asp187 and Gln190 in the Na+/proline symporter (PutP) of Escherichia coli

Asp187 and Gln190 were predicted as conserved and closely located at the Na(+) binding site in a topology and homology model structure of Na(+)/proline symporter (PutP) of Escherichia coli. The replacement of Asp187 with Ala or Leu did not affect proline transport activity; whereas, change to Gln abolished the active transport. The binding affinity for Na(+) or proline of these mutants was similar to that of wild-type (WT) PutP. This result indicates Asp187 to be responsible for active transport of proline without affecting the binding. Replacement of Gln190 with Ala, Asn, Asp, Leu and Glu had no effect on transport or binding, suggesting that it may not have a role in the transport. However, in the negative D187Q mutant, a second mutation, of Gln190 to Glu or Leu, restored 46 or 7% of the transport activity of WT, respectively, while mutation to Ala, Asn or Asp had no effect. Thus, side chain at position 190 has a crucial role in suppressing the functional defect of the D187Q mutant. We conclude that Asp187 is responsible for transport activity instead of coupling-ion binding by constituting the translocation pathway of the ion and Gln190 provides a suppressing mutation site to regain PutP functional activity.

The mechanism of substrate release by the aspartate transporter GltPh: insights from simulations

Glutamate transporters regulate excitatory amino acid neurotransmission across neuronal and glial cell membranes by coupling the translocation of their substrate (aspartate or glutamate) into the intracellular (IC) medium to the energetically favorable transport of sodium ions or other cations. The first crystallographically resolved structure of this family, the archaeal aspartate transporter, Glt(Ph), has served as a structural paradigm for elucidating the mechanism of substrate translocation by these transporters. Two helical hairpins, HP2 and HP1, at the core domains of the three subunits that form this membrane protein have been proposed to act as the respective extracellular and IC gates for substrate intake and release. Molecular dynamics simulations using the outward-facing structure have confirmed that the HP2 loop acts as an EC gate. The mechanism of substrate release at atomic scale, however, remained unknown due to the lack of structural data until the recent determination of the inward-facing structure of Glt(Ph). In the present study, we use this recently resolved structure to simulate the release of substrate to the cytoplasm and the roles of HP1 and HP2 in this process. The highly flexible HP2 loop is observed to serve as an activator (or initiator) prompting the release of a gatekeeper Na(+) to the cytoplasm and promoting the influx of water molecules from the cytoplasm, which effectively disrupt substrate-protein interactions and drive the dislodging of the substrate from its binding site. The completion of substrate release and exit, however, entails the opening of the highly stable HP1 loop as well. Overall, the unique conformational flexibility of the HP2 loop, the dissociation of a Na(+), the hydration of binding pocket, and final yielding of the HP1 loop 3-Ser motif emerge as the successive events controlling the release of the bound substrate to the cell interior by glutamate transporters.

Modeling and simulation of ion-coupled and ATP-driven membrane proteins

The molecular mechanisms of membrane proteins that are activated either by ions or by ATP are just beginning to come into focus, as long-awaited structural data are revealed. This information is being leveraged and supplemented to great effect by molecular modeling and computer simulation studies. Important examples include the homology modeling of eukaryotic protein structures based on distantly related templates, as well as the use of internal structural symmetry for modeling different states in conformational cycles. Molecular simulation studies have elucidated the location and coordination structure of ion binding sites, and explained their selectivity, while also providing tantalizing insights into the mechanisms that couple conformational change to ion translocation or ATP hydrolysis.

Reconstructing a chloride-binding site in a bacterial neurotransmitter transporter homologue

In ion-coupled transport proteins, occupation of selective ion-binding sites is required to trigger conformational changes that lead to substrate translocation. Neurotransmitter transporters, targets of abused and therapeutic drugs, require Na(+) and Cl(-) for function. We recently proposed a chloride-binding site in these proteins not present in Cl(-)-independent prokaryotic homologues. Here we describe conversion of the Cl(-)-independent prokaryotic tryptophan transporter TnaT to a fully functional Cl(-)-dependent form by a single point mutation, D268S. Mutations in TnaT-D268S, in wild type TnaT and in serotonin transporter provide direct evidence for the involvement of each of the proposed residues in Cl(-) coordination. In both SERT and TnaT-D268S, Cl(-) and Na(+) mutually increased each other's potency, consistent with electrostatic interaction through adjacent binding sites. These studies establish the site where Cl(-) binds to trigger conformational change during neurotransmitter transport.

The substrate-driven transition to an inward-facing conformation in the functional mechanism of the dopamine transporter

Background

The dopamine transporter (DAT), a member of the neurotransmitter:Na⁺ symporter (NSS) family, terminates dopaminergic neurotransmission and is a major molecular target for psychostimulants such as cocaine and amphetamine, and for the treatment of attention deficit disorder and depression. The crystal structures of the prokaryotic NSS homolog of DAT, the leucine transporter LeuT, have provided critical structural insights about the occluded and outward-facing conformations visited during the substrate transport, but only limited clues regarding mechanism. To understand the transport mechanism in DAT we have used a homology model based on the LeuT structure in a computational protocol validated previously for LeuT, in which steered molecular dynamics (SMD) simulations guide the substrate along a pathway leading from the extracellular end to the intracellular (cytoplasmic) end.

Methodology/Principal Findings

Key findings are (1) a second substrate binding site in the extracellular vestibule, and (2) models of the conformational states identified as occluded, doubly occupied, and inward-facing. The transition between these states involve a spatially ordered sequence of interactions between the two substrate-binding sites, followed by rearrangements in structural elements located between the primary binding site and the cytoplasmic end. These rearrangements are facilitated by identified conserved hinge regions and a reorganization of interaction networks that had been identified as gates.

Conclusions/Significance

Computational simulations supported by information available from experiments in DAT and other NSS transporters have produced a detailed mechanistic proposal for the dynamic changes associated with substrate transport in DAT. This allosteric mechanism is triggered by the binding of substrate in the S2 site in the presence of the substrate in the S1 site. Specific structural elements involved in this mechanism, and their roles in the conformational transitions illuminated here describe, a specific substrate-driven allosteric mechanism that is directly amenable to experiment as shown previously for LeuT.

Sugar binding residue affects apparent Na+ affinity and transport stoichiometry in mouse sodium/glucose cotransporter type 3B

SGLT1 is a sodium/glucose cotransporter that moves two Na(+) ions with each glucose molecule per cycle. SGLT3 proteins belong to the same family and are described as glucose sensors rather than glucose transporters. Thus, human SGLT3 (hSGLT3) does not transport sugar, but extracellular glucose depolarizes the cell in which it is expressed. Mouse SGLT3b (mSGLT3b), although it transports sugar, has low apparent sugar affinity and partially uncoupled stoichiometry compared with SGLT1, suggesting that mSGLT3b is also a sugar sensor. The crystal structure of the Vibrio parahaemolyticus SGLT showed that residue Gln(428) interacts directly with the sugar. The corresponding amino acid in mammalian proteins, 457, is conserved in all SGLT1 proteins as glutamine. In SGLT3 proteins, glutamate is the most common residue at this position, although it is a glycine in mSGLT3b and a serine in rat SGLT3b. To test the contribution of this residue to the function of SGLT3 proteins, we constructed SGLT3b mutants that recapitulate residue 457 in SGLT1 and hSGLT3, glutamine and glutamate, respectively. The presence of glutamine at residue 457 increased the apparent Na(+) and sugar affinities, whereas glutamate decreased the apparent Na(+) affinity. Moreover, glutamate transported more cations per sugar molecule than the wild type protein. We propose a model where cations are released intracellularly without the release of sugar from an intermediate state. This model explains the uncoupled charge:sugar transport phenotype observed in wild type and G457E-mSGLT3b compared with SGLT1 and the sugar-activated cation transport without sugar transport that occurs in hSGLT3.

Mechanism of melibiose/cation symport of the melibiose permease of Salmonella typhimurium

The MelB permease of Salmonella typhimurium (MelB-ST) catalyzes the coupled symport of melibiose and Na(+), Li(+), or H(+). In right-side-out membrane vesicles, melibiose efflux is inhibited by an inwardly directed gradient of Na(+) or Li(+) and stimulated by equimolar concentrations of internal and external Na(+) or Li(+). Melibiose exchange is faster than efflux in the presence of H(+) or Na(+) and stimulated by an inwardly directed Na(+) gradient. Thus, sugar is released from MelB-ST externally prior to the release of cation in agreement with current models proposed for MelB of Escherichia coli (MelB-EC) and LacY. Although Li(+) stimulates efflux, and an outwardly directed Li(+) gradient increases exchange, it is striking that internal and external Li(+) with no gradient inhibits exchange. Furthermore, Trp → dansyl FRET measurements with a fluorescent sugar (2'-(N-dansyl)aminoalkyl-1-thio-β-D-galactopyranoside) demonstrate that MelB-ST, in the presence of Na(+) or Li(+), exhibits (app)K(d) values of ∼1 mM for melibiose. Na(+) and Li(+) compete for a common binding pocket with activation constants for FRET of ∼1 mM, whereas Rb(+) or Cs(+) exhibits little or no effect. Taken together, the findings indicate that MelB-ST utilizes H(+) in addition to Na(+) and Li(+). FRET studies also show symmetrical emission maximum at ∼500 nm with MelB-ST in the presence of 2'-(N-dansyl)aminoalkyl-1-thio-β-D-galactopyranoside and Na(+), Li(+), or H(+), which implies a relatively homogeneous distribution of conformers of MelB-ST ternary complexes in the membrane.

The mechanism of sodium and substrate release from the binding pocket of vSGLT

Membrane co-transport proteins that use a five-helix inverted repeat motif have recently emerged as one of the largest structural classes of secondary active transporters. However, despite many structural advances there is no clear evidence of how ion and substrate transport are coupled. Here we report a comprehensive study of the sodium/galactose transporter from Vibrio parahaemolyticus (vSGLT), consisting of molecular dynamics simulations, biochemical characterization and a new crystal structure of the inward-open conformation at a resolution of 2.7 Å. Our data show that sodium exit causes a reorientation of transmembrane helix 1 that opens an inner gate required for substrate exit, and also triggers minor rigid-body movements in two sets of transmembrane helical bundles. This cascade of events, initiated by sodium release, ensures proper timing of ion and substrate release. Once set in motion, these molecular changes weaken substrate binding to the transporter and allow galactose readily to enter the intracellular space. Additionally, we identify an allosteric pathway between the sodium-binding sites, the unwound portion of transmembrane helix 1 and the substrate-binding site that is essential in the coupling of co-transport.

Homology model of the Na+/proline transporter PutP of Escherichia coli and its functional implications

Na(+)/solute symporters are essential membrane integrated proteins that couple the flow of Na(+) ions driven by electrochemical Na(+) gradients to the transport of solutes across biological membranes. Here, we used a combination of molecular modeling techniques and evolutionary conservation analysis to construct and validate a first model of the Na(+)/proline symporter PutP of Escherichia coli based on the crystal structure of the bacterial Na(+)/galactose symporter vSGLT. Ligand docking experiments were employed to gain information about residues involved in proline binding. The proposed model is consistent with the available experimental data and was further validated by amino acid substitutions and kinetic and protein chemical analyses. Combination of the results of molecular modeling and functional studies predicts the location and organization of the Na(+) and proline binding sites. Remarkably, as proposed computationally and discovered here experimentally, residues Y140, W244, and Y248 of transmembrane segments 4 and 7 are found to be particularly important for PutP function and suggested to participate in proline binding and/or gating.

Capturing Functional Motions of Membrane Channels and Transporters with Molecular Dynamics Simulation

Conformational changes of proteins are involved in all aspects of protein function in biology. Almost all classes of proteins respond to changes in their environment, ligand binding, and interaction with other proteins and regulatory agents through undergoing conformational changes of various degrees and magnitudes. Membrane channels and transporters are the major classes of proteins that are responsible for mediating efficient and selective transport of materials across the cellular membrane. Similar to other proteins, they take advantage of conformational changes to make transitions between various functional states. In channels, large-scale conformational changes are mostly involved in the process of "gating", i.e., opening and closing of the pore of the channel protein in response to various signals. In transporters, conformational changes constitute various steps of the conduction process, and, thus, are more closely integrated in the transport process. Owing to significant progress in developing highly efficient parallel algorithms in molecular dynamics simulations and increased computational resources, and combined with the availability of high-resolution, atomic structures of membrane proteins, we are in an unprecedented position to use computer simulation and modeling methodologies to investigate the mechanism of function of membrane channels and transporters. While the entire transport cycle is still out of reach of current methodologies, many steps involved in the function of transport proteins have been characterized with molecular dynamics simulations. Here, we present several examples of such studies from our laboratory, in which functionally relevant conformational changes of membrane channels and transporters have been characterized using extended simulations.

The structural basis of secondary active transport mechanisms

Secondary active transporters couple the free energy of the electrochemical potential of one solute to the transmembrane movement of another. As a basic mechanistic explanation for their transport function the model of alternating access was put forward more than 40 years ago, and has been supported by numerous kinetic, biochemical and biophysical studies. According to this model, the transporter exposes its substrate binding site(s) to one side of the membrane or the other during transport catalysis, requiring a substantial conformational change of the carrier protein. In the light of recent structural data for a number of secondary transport proteins, we analyze the model of alternating access in more detail, and correlate it with specific structural and chemical properties of the transporters, such as their assignment to different functional states in the catalytic cycle of the respective transporter, the definition of substrate binding sites, the type of movement of the central part of the carrier harboring the substrate binding site, as well as the impact of symmetry on fold-specific conformational changes. Besides mediating the transmembrane movement of solutes, the mechanism of secondary carriers inherently involves a mechanistic coupling of substrate flux to the electrochemical potential of co-substrate ions or solutes. Mainly because of limitations in resolution of available transporter structures, this important aspect of secondary transport cannot yet be substantiated by structural data to the same extent as the conformational change aspect. We summarize the concepts of coupling in secondary transport and discuss them in the context of the available evidence for ion binding to specific sites and the impact of the ions on the conformational state of the carrier protein, which together lead to mechanistic models for coupling.

Water permeation through the sodium-dependent galactose cotransporter vSGLT

It is well accepted that cotransporters facilitate water movement by two independent mechanisms: osmotic flow through a water channel in the protein and flow driven by ion/substrate cotransport. However, the molecular mechanism of transport-linked water flow is controversial. Some researchers believe that it occurs via cotransport, in which water is pumped along with the transported cargo, while others believe that flow is osmotic in response to an increase in intracellular osmolarity. In this letter, we report the results of a 200-ns molecular dynamics simulation of the sodium-dependent galactose cotransporter vSGLT. Our simulation shows that a significant number of water molecules cross the protein through the sugar-binding site in the presence as well as the absence of galactose, and 70-80 water molecules accompany galactose as it moves from the binding site into the intracellular space. During this event, the majority of water molecules in the pathway are unable to diffuse around the galactose, resulting in water in the inner half of the transporter being pushed into the intracellular space and replaced by extracellular water. Thus, our simulation supports the notion that cotransporters act as both passive water channels and active water pumps with the transported substrate acting as a piston to rectify the motion of water.

Modeling and dynamics of the inward-facing state of a Na+/Cl- dependent neurotransmitter transporter homologue

The leucine transporter (LeuT) has recently commanded exceptional attention due mainly to two distinctions; it provides the only crystal structures available for a protein homologous to the pharmacologically relevant neurotransmitter: sodium symporters (NSS), and, it exhibits a hallmark 5-TM inverted repeat ("LeuT-fold"), a fold recently discovered to also exist in several secondary transporter families, underscoring its general role in transporter function. Constructing the transport cycle of "LeuT-fold" transporters requires detailed structural and dynamic descriptions of the outward-facing (OF) and inward-facing (IF) states, as well as the intermediate states. To this end, we have modeled the structurally unknown IF state of LeuT, based on the known crystal structures of the OF state of LeuT and the IF state of vSGLT, a "LeuT-fold" transporter. The detailed methodology developed for the study combines structure-based alignment, threading, targeted MD and equilibrium MD, and can be applied to other proteins. The resulting IF-state models maintain the secondary structural features of LeuT. Water penetration and solvent accessibility calculations show that TM1, TM3, TM6 and TM8 line the substrate binding/unbinding pathway with TM10 and its pseudosymmetric partner, TM5, participating in the extracellular and intracellular halves of the lumen, respectively. We report conformational hotspots where notable changes in interactions occur between the IF and OF states. We observe Na2 exiting the LeuT-substrate- complex in the IF state, mainly due to TM1 bending. Inducing a transition in only one of the two pseudosymmetric domains, while allowing the second to respond dynamically, is found to be sufficient to induce the formation of the IF state. We also propose that TM2 and TM7 may be facilitators of TM1 and TM6 motion. Thus, this study not only presents a novel modeling methodology applied to obtain the IF state of LeuT, but also describes structural elements involved in a possibly general transport mechanism in transporters adopting the "LeuT-fold".

Structural perspectives on secondary active transporters

Secondary active transporters catalyze the concentrative transport of substrates across lipid membranes by harnessing the energy of electrochemical ion gradients. These transporters bind their ligands on one side of the membrane, and undergo a global conformational change to release them on the other side of the membrane. Over the last few years, crystal structures have captured several bacterial secondary transporters in different states along their transport cycle, providing insight into possible molecular mechanisms. In this review, we summarize recent findings focusing on the emerging structural and mechanistic similarities between evolutionary diverse transporters. We also discuss the structural basis of substrate binding, ion coupling and inhibition viewed from the perspective of these similarities.