Structural determinants of water permeation through the sodium-galactose transporter vSGLT

A Futile Cycle?: Tissue Homeostatic Trans-Membrane Water Co-Transport: Kinetics, Thermodynamics, Metabolic Consequences

The phenomenon of active trans-membrane water cycling (AWC) has emerged in little over a decade. Here, we consider H2O transport across cell membranes from the origins of its study. Historically, trans-membrane water transport processes were classified into: A) compensating bidirectional fluxes ("exchange"), and B) unidirectional flux ("net flow") categories. Recent literature molecular structure determinations and molecular dynamic (MD) simulations indicate probably all the many different hydrophilic substrate membrane co-transporters have membrane-spanning hydrophilic pathways and co-transport water along with their substrates, and that they individually catalyze category A and/or B water flux processes, although usually not simultaneously. The AWC name signifies that, integrated over the all the cell's co-transporters, the rate of homeostatic, bidirectional trans-cytolemmal water exchange (category A) is synchronized with the metabolic rate of the crucial Na+,K+-ATPase (NKA) enzyme. A literature survey indicates the stoichiometric (category B) water/substrate ratios of individual co-transporters are often very large. The MD simulations also suggest how different co-transporter reactions can be kinetically coupled molecularly. Is this (Na+,K+-ATPase rate-synchronized) cycling futile, or is it consequential? Conservatively representative literature metabolomic and proteinomic results enable comprehensive free energy analyses of the many transport reactions with known water stoichiometries. Free energy calculations, using literature intracellular pressure (Pi) values reveals there is an outward trans-membrane H2O barochemical gradient of magnitude comparable to that of the well-known inward Na+ electrochemical gradient. For most co-influxers, these gradients are finely balanced to maintain intracellular metabolite concentration values near their consuming enzyme Michaelis constants. The thermodynamic analyses include glucose, glutamate-, gamma-aminobutyric acid (GABA), and lactate- transporters. 2%-4% Pi alterations can lead to disastrous concentration levels. For the neurotransmitters glutamate- and GABA, very small astrocytic Pi changes can allow/disallow synaptic transmission. Unlike the Na+ and K+ electrochemical steady-states, the H2O barochemical steady-state is in (or near) chemical equilibrium. The analyses show why the presence of aquaporins (AQPs) does not dissipate the trans-membrane pressure gradient. A feedback loop inherent in the opposing Na+ electrochemical and H2O barochemical gradients regulates AQP-catalyzed water flux as an integral AWC aspect. These results also require a re-consideration of the underlying nature of Pi. Active trans-membrane water cycling is not futile, but is inherent to the cell's "NKA system" - a new, fundamental aspect of biology.

Collective Domain Motion Facilitates Water Transport in SGLT1

The human sodium-glucose cotransporter protein (SGLT1) is an important representative of the sodium solute symporters belonging to the secondary active transporters that are critical to the homeostasis of sugar, sodium, and water in the cell. The underlying transport mechanism of SGLT1 is based on switching between inward- and outward-facing conformations, known as the alternating access model, which is crucial for substrate transport, and has also been postulated for water permeation. However, the nature of water transport remains unclear and is disputed along the passive and active transport, with the latter postulating the presence of the pumping effect. To better examine the water transport in SGLT1, we performed a series of equilibrium all-atom molecular dynamics simulations, totaling over 6 μs of sample representative conformational states of SGLT1 and its complexes, with the natural substrates, ions, and inhibitors. In addition to elucidating the basic physical factors influencing water permeation, such as channel openings and energetics, we focus on dynamic flexibility and its relationship with domain motion. Our results clearly demonstrate a dependence of instantaneous water flux on the channel opening and local water diffusion in the channel, strongly supporting the existence of a passive water transport in SGLT1. In addition, a strong correlation found between the local water diffusion and protein domain motion, resembling the "rocking-bundle" motion, reveals its facilitating role in the water transport.

Arginine transportation mechanism through cationic amino acid transporter 1: insights from molecular dynamics studies

Metabolic and signaling mechanisms in mammalian cells are facilitated by the transportation of L-arginine (Arg) across the plasma membrane through cationic amino acid transporter (CAT) proteins. Due to a lack of argininosuccinate synthase (ASS) activity in various tumor cells such as acute myeloid leukemia, acute lymphocytic leukemia, and chronic lymphocytic leukemia, these tumor entities are arginine-auxotrophic and therefore depend on the uptake of the amino acid arginine. Cationic amino acid transporter-1 (CAT-1) is the leading arginine importer expressed in the aforementioned tumor entities. Hence, in the present study, to investigate the transportation mechanism of arginine in CAT-1, we performed molecular dynamics (MD) simulation methods on the modeled human CAT-1. The MM-PBSA approach was conducted to determine the critical residues interacting with arginine within the corresponding binding site of CAT-1. In addition, we found out that the water molecules have the leading role in forming the transportation channel within CAT-1. The conductive structure of CAT-1 was formed only when the water molecules were continuously distributed across the channel. Steered molecular dynamics (SMD) simulation approach showed various energy barriers against arginine transportation through CAT-1, especially while crossing the bottlenecks of the related channel. These findings at the molecular level might shed light on identifying the crucial amino acids in the binding of arginine to eukaryotic CATs and also provide fundamental insights into the arginine transportation mechanisms through CAT-1. Understanding the transportation mechanism of arginine is essential to developing CAT-1 blockers, which can be potential medications for some types of cancers.Communicated by Ramaswamy H. Sarma.

Modeling of SGLT1 in Reconstituted Systems Reveals Apparent Ion-Dependencies of Glucose Uptake and Strengthens the Notion of Water-Permeable Apo States

The reconstitution of secondary active transporters into liposomes shed light on their molecular transport mechanism. The latter are either symporters, antiporters or exchangers, which use the energy contained in the electrochemical gradient of ions to fuel concentrative uptake of their cognate substrate. In liposomal preparations, these gradients can be set by the experimenter. However, due to passive diffusion of the ions and solutes through the membrane, the gradients are not stable and little is known on the time course by which they dissipate and how the presence of a transporter affects this process. Gradient dissipation can also generate a transmembrane potential (VM). Because it is the effective ion gradient, which together with VM fuels concentrative uptake, knowledge on how these parameters change within the time frame of the conducted experiment is key to understanding experimental outcomes. Here, we addressed this problem by resorting to a modelling approach. To this end, we mathematically modeled the liposome in the assumed presence and absence of the sodium glucose transporter 1 (SGLT1). We show that 1) the model can prevent us from reaching erroneous conclusions on the driving forces of substrate uptake and we 2) demonstrate utility of the model in the assignment of the states of SGLT1, which harbor a water channel.

Cerebrospinal fluid production by the choroid plexus: a century of barrier research revisited

Cerebrospinal fluid (CSF) envelops the brain and fills the central ventricles. This fluid is continuously replenished by net fluid extraction from the vasculature by the secretory action of the choroid plexus epithelium residing in each of the four ventricles. We have known about these processes for more than a century, and yet the molecular mechanisms supporting this fluid secretion remain unresolved. The choroid plexus epithelium secretes its fluid in the absence of a trans-epithelial osmotic gradient, and, in addition, has an inherent ability to secrete CSF against an osmotic gradient. This paradoxical feature is shared with other 'leaky' epithelia. The assumptions underlying the classical standing gradient hypothesis await experimental support and appear to not suffice as an explanation of CSF secretion. Here, we suggest that the elusive local hyperosmotic compartment resides within the membrane transport proteins themselves. In this manner, the battery of plasma membrane transporters expressed in choroid plexus are proposed to sustain the choroidal CSF secretion independently of the prevailing bulk osmotic gradient.

Paracellular Filtration Secretion Driven by Mechanical Force Contributes to Small Intestinal Fluid Dynamics

Studies of fluid secretion by the small intestine are dominated by the coupling with ATP-dependent generation of ion gradients, whereas the contribution of filtration secretion has been overlooked, possibly by the lack of a known mechanistic basis. We measured apical fluid flow and generation of hydrostatic pressure gradients by epithelia of cultured mouse enterocytes, Caco-2 and T-84 cells, and fibroblasts exposed to mechanical force provided by vigorous aeration and in response to ion gradients, inhibitors of ion channels and transporters and in vitro using intact mouse and rat small intestine. We describe herein a paracellular pathway for unidirectional filtration secretion that is driven by mechanical force, requires tight junctions, is independent of ionic and osmotic gradients, generates persistent hydrostatic pressure gradients, and would contribute to the fluid shifts that occur during digestion and diarrhea. Zinc inhibits the flow of fluid and the paracellular marker fluorescein isothyocyanate conjugated dextran (MW = 4 kD) across epithelia of cultured enterocytes (>95%; p < 0.001) and intact small intestine (>40%; p = 0.03). We propose that mechanical force drives fluid secretion through the tight junction complex via a “one-way check valve” that can be regulated. This pathway of filtration secretion complements chloride-coupled fluid secretion during high-volume fluid flow. The role of filtration secretion in the genesis of diarrhea in intact animals needs further study. Our findings may explain a potential linkage between intestinal motility and intestinal fluid dynamics.

A kinetic mechanism for enhanced selectivity of membrane transport

Membrane transport is generally thought to occur via an alternating access mechanism in which the transporter adopts at least two states, accessible from two different sides of the membrane to exchange substrates from the extracellular environment and the cytoplasm or from the cytoplasm and the intracellular matrix of the organelles (only in eukaryotes). In recent years, a number of high resolution structures have supported this general framework for a wide class of transport molecules, although additional states along the transport pathway are emerging as critically important. Given that substrate binding is often weak in order to enhance overall transport rates, there exists the distinct possibility that transporters may transport the incorrect substrate. This is certainly the case for many pharmaceutical compounds that are absorbed in the gut or cross the blood brain barrier through endogenous transporters. Docking studies on the bacterial sugar transporter vSGLT reveal that many highly toxic compounds are compatible with binding to the orthosteric site, further motivating the selective pressure for additional modes of selectivity. Motivated by recent work in which we observed failed substrate delivery in a molecular dynamics simulation where the energized ion still goes down its concentration gradient, we hypothesize that some transporters evolved to harness this 'slip' mechanism to increase substrate selectivity and reduce the uptake of toxic molecules. Here, we test this idea by constructing and exploring a kinetic transport model that includes a slip pathway. While slip reduces the overall productive flux, when coupled with a second toxic molecule that is more prone to slippage, the overall substrate selectivity dramatically increases, suppressing the accumulation of the incorrect compound. We show that the mathematical framework for increased substrate selectivity in our model is analogous to the classic proofreading mechanism originally proposed for tRNA synthase; however, because the transport cycle is reversible we identified conditions in which the selectivity is essentially infinite and incorrect substrates are exported from the cell in a 'detoxification' mode. The cellular consequences of proofreading and membrane slippage are discussed as well as the impact on future drug development.

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.

Computational Estimation of Microsecond to Second Atomistic Folding Times

Despite the development of massively parallel computing hardware including inexpensive graphics processing units (GPUs), it has remained infeasible to simulate the folding of atomistic proteins at room temperature using conventional molecular dynamics (MD) beyond the microsecond scale. Here, we report the folding of atomistic, implicitly solvated protein systems with folding times τ ranging from ∼10 μs to ∼100 ms using the weighted ensemble (WE) strategy in combination with GPU computing. Starting from an initial structure or set of structures, WE organizes an ensemble of GPU-accelerated MD trajectory segments via intermittent pruning and replication events to generate statistically unbiased estimates of rate constants for rare events such as folding; no biasing forces are used. Although the variance among atomistic WE folding runs is significant, multiple independent runs are used to reduce and quantify statistical uncertainty. Folding times are estimated directly from WE probability flux and from history-augmented Markov analysis of the WE data. Three systems were examined: NTL9 at low solvent viscosity (yielding τf = 0.8-9 μs), NTL9 at water-like viscosity (τf = 0.2-2 ms), and Protein G at low viscosity (τf = 3-200 ms). In all cases, the folding time, uncertainty, and ensemble properties could be estimated from WE simulation; for Protein G, this characterization required significantly less overall computing than would be required to observe a single folding event with conventional MD simulations. Our results suggest that the use and calibration of force fields and solvent models for precise estimation of kinetic quantities is becoming feasible.

Novel and Unexpected Functions of SGLTs

It has been 30 years since the intestinal sodium glucose cotransporter SGLT1 was cloned, and, in the intervening years, there have been many advances that have influenced physiology and medicine. Among the first was that SGLT1 is the founding member of the human gene family SLC5, containing 11 diverse transporters and a glucose sensor. Equally surprising was that SGLTs are members of a structural family of cotransporters and exchangers in different gene families. This led to the conclusion that these proteins operate by a mechanism where transport involves the opening and closing of external and internal gates. The mechanism is shared by a wide variety of transporters in different structural families, e.g., the human facilitated glucose transporters (SLC2) in the huge major facilitator superfamily (MFS). Not surprising is the finding that mutations in Sglt genes cause the rare diseases glucose-galactose-malabsorption (GGM) and familial renal glucosuria (FRG). However, it was not envisaged that SGLT inhibitors would be used to treat diabetes mellitus, and these drugs may be able to treat cancer. Finally, in 2017, we have just learned that SGLT1 may be required to resist infection and to avoid recurrent pregnancy loss.

Directional Fluid Transport across Organ-Blood Barriers: Physiology and Cell Biology

Directional fluid flow is an essential process for embryo development as well as for organ and organism homeostasis. Here, we review the diverse structure of various organ-blood barriers, the driving forces, transporters, and polarity mechanisms that regulate fluid transport across them, focusing on kidney-, eye-, and brain-blood barriers. We end by discussing how cross talk between barrier epithelial and endothelial cells, perivascular cells, and basement membrane signaling contribute to generate and maintain organ-blood barriers.

Structural and functional significance of water permeation through cotransporters

Membrane transporters, in addition to their major role as specific carriers for ions and small molecules, can also behave as water channels. However, neither the location of the water pathway in the protein nor their functional importance is known. Here, we map the pathway for water and urea through the intestinal sodium/glucose cotransporter SGLT1. Molecular dynamics simulations using the atomic structure of the bacterial transporter vSGLT suggest that water permeates the same path as Na+ and sugar. On a structural model of SGLT1, based on the homology structure of vSGLT, we identified and mutated residues lining the sugar transport pathway to cysteine. The mutants were expressed in Xenopus oocytes, and the unitary water and urea permeabilities were determined before and after modifying the cysteine side chain with reversible methanethiosulfonate reagents. The results demonstrate that water and urea follow the sugar transport pathway through SGLT1. The changes in permeability, increases or decreases, with side-chain modifications depend on the location of the mutation in the region of external or internal gates, or the sugar binding site. These changes in permeability are hypothesized to be due to alterations in steric hindrance to water and urea, and/or changes in protein folding caused by mismatching of side chains in the water pathway. Water permeation through SGLT1 and other transporters bears directly on the structural mechanism for the transport of polar solutes through these proteins. Finally, in vitro experiments on mouse small intestine show that SGLT1 accounts for two-thirds of the passive water flow across the gut.

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.

Compressing molecular dynamics trajectories: Breaking the one-bit-per-sample barrier

Molecular dynamics simulations yield large amounts of trajectory data. For their durable storage and accessibility an efficient compression algorithm is paramount. State of the art domain-specific algorithms combine quantization, Huffman encoding and occasionally domain knowledge. We propose the high resolution trajectory compression scheme (HRTC) that relies on piecewise linear functions to approximate quantized trajectories. By splitting the error budget between quantization and approximation, our approach beats the current state of the art by several orders of magnitude given the same error tolerance. It allows storing samples at far less than one bit per sample. It is simple and fast enough to be integrated into the inner simulation loop, store every time step, and become the primary representation of trajectory data. © 2016 Wiley Periodicals, Inc.

The Sodium Glucose Cotransporter SGLT1 Is an Extremely Efficient Facilitator of Passive Water Transport

The small intestine is void of aquaporins adept at facilitating vectorial water transport, and yet it reabsorbs ∼8 liters of fluid daily. Implications of the sodium glucose cotransporter SGLT1 in either pumping water or passively channeling water contrast with its reported water transporting capacity, which lags behind that of aquaporin-1 by 3 orders of magnitude. Here we overexpressed SGLT1 in MDCK cell monolayers and reconstituted the purified transporter into proteoliposomes. We observed the rate of osmotic proteoliposome deflation by light scattering. Fluorescence correlation spectroscopy served to assess (i) SGLT1 abundance in both vesicles and plasma membranes and (ii) flow-mediated dilution of an aqueous dye adjacent to the cell monolayer. Calculation of the unitary water channel permeability, pf, yielded similar values for cell and proteoliposome experiments. Neither the absence of glucose or Na(+), nor the lack of membrane voltage in vesicles, nor the directionality of water flow grossly altered pf Such weak dependence on protein conformation indicates that a water-impermeable occluded state (glucose and Na(+) in their binding pockets) lasts for only a minor fraction of the transport cycle or, alternatively, that occlusion of the substrate does not render the transporter water-impermeable as was suggested by computational studies of the bacterial homologue vSGLT. Although the similarity between the pf values of SGLT1 and aquaporin-1 makes a transcellular pathway plausible, it renders water pumping physiologically negligible because the passive flux would be orders of magnitude larger.

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.

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.

Membrane Protein Properties Revealed through Data-Rich Electrostatics Calculations

The electrostatic properties of membrane proteins often reveal many of their key biophysical characteristics, such as ion channel selectivity and the stability of charged membrane-spanning segments. The Poisson-Boltzmann (PB) equation is the gold standard for calculating protein electrostatics, and the software APBSmem enables the solution of the PB equation in the presence of a membrane. Here, we describe significant advances to APBSmem, including full automation of system setup, per-residue energy decomposition, incorporation of PDB2PQR, calculation of membrane-induced pKa shifts, calculation of non-polar energies, and command-line scripting for large-scale calculations. We highlight these new features with calculations carried out on a number of membrane proteins, including the recently solved structure of the ion channel TRPV1 and a large survey of 1,614 membrane proteins of known structure. This survey provides a comprehensive list of residues with large electrostatic penalties for being embedded in the membrane, potentially revealing interesting functional information.

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.

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.