Mitochondrial uncouplers induce proton leak by activating AAC and UCP1.

Mitochondria generate heat due to H+ leak (IH) across their inner membrane1. IH results from the action of long-chain fatty acids on uncoupling protein 1 (UCP1) in brown fat2-6 and ADP/ATP carrier (AAC) in other tissues1,7-9, but the underlying mechanism is poorly understood. As evidence of pharmacological activators of IH through UCP1 and AAC is lacking, IH is induced by protonophores such as 2,4-dinitrophenol (DNP) and cyanide-4-(trifluoromethoxy) phenylhydrazone (FCCP)10,11. Although protonophores show potential in combating obesity, diabetes and fatty liver in animal models12-14, their clinical potential for treating human disease is limited due to indiscriminately increasing H+ conductance across all biological membranes10,11 and adverse side effects15. Here we report the direct measurement of IH induced by DNP, FCCP and other common protonophores and find that it is dependent on AAC and UCP1. Using molecular structures of AAC, we perform a computational analysis to determine the binding sites for protonophores and long-chain fatty acids, and find that they overlap with the putative ADP/ATP-binding site. We also develop a mathematical model that proposes a mechanism of uncoupler-dependent IH through AAC. Thus, common protonophoric uncouplers are synthetic activators of IH through AAC and UCP1, paving the way for the development of new and more specific activators of these two central mediators of mitochondrial bioenergetics.

Regulation of the mammalian-brain V-ATPase through ultraslow mode-switching.

Vacuolar-type adenosine triphosphatases (V-ATPases)1-3 are electrogenic rotary mechanoenzymes structurally related to F-type ATP synthases4,5. They hydrolyse ATP to establish electrochemical proton gradients for a plethora of cellular processes1,3. In neurons, the loading of all neurotransmitters into synaptic vesicles is energized by about one V-ATPase molecule per synaptic vesicle6,7. To shed light on this bona fide single-molecule biological process, we investigated electrogenic proton-pumping by single mammalian-brain V-ATPases in single synaptic vesicles. Here we show that V-ATPases do not pump continuously in time, as suggested by observing the rotation of bacterial homologues8 and assuming strict ATP-proton coupling. Instead, they stochastically switch between three ultralong-lived modes: proton-pumping, inactive and proton-leaky. Notably, direct observation of pumping revealed that physiologically relevant concentrations of ATP do not regulate the intrinsic pumping rate. ATP regulates V-ATPase activity through the switching probability of the proton-pumping mode. By contrast, electrochemical proton gradients regulate the pumping rate and the switching of the pumping and inactive modes. A direct consequence of mode-switching is all-or-none stochastic fluctuations in the electrochemical gradient of synaptic vesicles that would be expected to introduce stochasticity in proton-driven secondary active loading of neurotransmitters and may thus have important implications for neurotransmission. This work reveals and emphasizes the mechanistic and biological importance of ultraslow mode-switching.

Cilia-Associated Oxysterols Activate Smoothened.

Primary cilia are required for Smoothened to transduce vertebrate Hedgehog signals, but how Smoothened accumulates in cilia and is activated is incompletely understood. Here, we identify cilia-associated oxysterols that promote Smoothened accumulation in cilia and activate the Hedgehog pathway. Our data reveal that cilia-associated oxysterols bind to two distinct Smoothened domains to modulate Smoothened accumulation in cilia and tune the intensity of Hedgehog pathway activation. We find that the oxysterol synthase HSD11ß2 participates in the production of Smoothened-activating oxysterols and promotes Hedgehog pathway activity. Inhibiting oxysterol biosynthesis impedes oncogenic Hedgehog pathway activation and attenuates the growth of Hedgehog pathway-associated medulloblastoma, suggesting that targeted inhibition of Smoothened-activating oxysterol production may be therapeutically useful for patients with Hedgehog-associated cancers.

Elucidation of the molecular interactions that enable stable assembly and structural diversity in multicomponent immune receptors.

Multicomponent immune receptors are essential complexes in which distinct ligand-recognition and signaling subunits are held together by interactions between acidic and basic residues of their transmembrane helices. A 2:1 acidic-to-basic motif in the transmembrane domains of the subunits is necessary and sufficient to assemble these receptor complexes. Here, we study a prototype for these receptors, a DAP12-NKG2C 2:1 heterotrimeric complex, in which the two DAP12 subunits each contribute a single transmembrane Asp residue, and the NKG2C subunit contributes a Lys to form the complex. DAP12 can also associate with 20 other subunits using a similar motif. Here, we use molecular-dynamics simulations to understand the basis for the high affinity and diversity of interactions in this group of receptors. Simulations of the transmembrane helices with differing protonation states of the Asp-Asp-Lys triad identified a structurally stable interaction in which a singly-protonated Asp-Asp pair forms a hydrogen-bonded carboxyl-carboxylate clamp that clasps onto a charged Lys side chain. This polar motif was also supported by density functional theory and a Protein Data Bank-wide search. In contrast, the helices are dynamic at sites distal to the stable carboxyl-carboxylate clamp motif. Such a locally stable but globally dynamic structure is well suited to accommodate the sequence and structural variations in the transmembrane helices of multicomponent receptors, which mix and match subunits to create combinatorial functional diversity from a limited number of subunits. It also supports a signaling mechanism based on multisubunit clustering rather than propagation of rigid conformational changes through the membrane.

Allosteric cooperation in a de novo-designed two-domain protein.

We describe the de novo design of an allosterically regulated protein, which comprises two tightly coupled domains. One domain is based on the DF (Due Ferri in Italian or two-iron in English) family of de novo proteins, which have a diiron cofactor that catalyzes a phenol oxidase reaction, while the second domain is based on PS1 (Porphyrin-binding Sequence), which binds a synthetic Zn-porphyrin (ZnP). The binding of ZnP to the original PS1 protein induces changes in structure and dynamics, which we expected to influence the catalytic rate of a fused DF domain when appropriately coupled. Both DF and PS1 are four-helix bundles, but they have distinct bundle architectures. To achieve tight coupling between the domains, they were connected by four helical linkers using a computational method to discover the most designable connections capable of spanning the two architectures. The resulting protein, DFP1 (Due Ferri Porphyrin), bound the two cofactors in the expected manner. The crystal structure of fully reconstituted DFP1 was also in excellent agreement with the design, and it showed the ZnP cofactor bound over 12 Å from the dimetal center. Next, a substrate-binding cleft leading to the diiron center was introduced into DFP1. The resulting protein acts as an allosterically modulated phenol oxidase. Its Michaelis-Menten parameters were strongly affected by the binding of ZnP, resulting in a fourfold tighter Km and a 7-fold decrease in kcat These studies establish the feasibility of designing allosterically regulated catalytic proteins, entirely from scratch.

Taking the Monte-Carlo gamble: How not to buckle under the pressure!

Consistent buckling distortions of a large membrane patch (200 × 200 Å) are observed during molecular dynamics (MD) simulations using the Monte-Carlo (MC) barostat in combination with a hard Lennard-Jones (LJ) cutoff. The buckling behavior is independent of both the simulation engine and the force field but requires the MC barostat-hard LJ cutoff combination. Similar simulations of a smaller patch (90 × 90 Å) do not show buckling, but do show a small, systematic reduction in the surface area accompanied by ~1 Å thickening suggestive of compression. We show that a mismatch in the way potentials and forces are handled in the dynamical equations versus the MC barostat results in a compressive load on the membrane. Moreover, a straightforward application of elasticity theory reveals that a minimal compression of the linear dimensions of the membrane, inversely proportional to the edge length, is required for buckling, explaining this differential behavior. We recommend always using LJ force or potential-switching when the MC barostat is employed to avoid undesirable membrane deformations.

Chemically induced vesiculation as a platform for studying TMEM16F activity.

Calcium-activated phospholipid scramblase mediates the energy-independent bidirectional translocation of lipids across the bilayer, leading to transient or, in the case of apoptotic scrambling, sustained collapse of membrane asymmetry. Cells lacking TMEM16F-dependent lipid scrambling activity are deficient in generation of extracellular vesicles (EVs) that shed from the plasma membrane in a Ca2+-dependent manner, namely microvesicles. We have adapted chemical induction of giant plasma membrane vesicles (GPMVs), which require both TMEM16F-dependent phospholipid scrambling and calcium influx, as a kinetic assay to investigate the mechanism of TMEM16F activity. Using the GPMV assay, we identify and characterize both inactivating and activating mutants that elucidate the mechanism for TMEM16F activation and facilitate further investigation of TMEM16F-mediated lipid translocation and its role in extracellular vesiculation.

A systems-biology approach to molecular machines: Exploration of alternative transporter mechanisms.

Motivated by growing evidence for pathway heterogeneity and alternative functions of molecular machines, we demonstrate a computational approach for investigating two questions: (1) Are there multiple mechanisms (state-space pathways) by which a machine can perform a given function, such as cotransport across a membrane? (2) How can additional functionality, such as proofreading/error-correction, be built into machine function using standard biochemical processes? Answers to these questions will aid both the understanding of molecular-scale cell biology and the design of synthetic machines. Focusing on transport in this initial study, we sample a variety of mechanisms by employing Metropolis Markov chain Monte Carlo. Trial moves adjust transition rates among an automatically generated set of conformational and binding states while maintaining fidelity to thermodynamic principles and a user-supplied fitness/functionality goal. Each accepted move generates a new model. The simulations yield both single and mixed reaction pathways for cotransport in a simple environment with a single substrate along with a driving ion. In a "competitive" environment including an additional decoy substrate, several qualitatively distinct reaction pathways are found which are capable of extremely high discrimination coupled to a leak of the driving ion, akin to proofreading. The array of functional models would be difficult to find by intuition alone in the complex state-spaces of interest.

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.

Conformational transitions of the sodium-dependent sugar transporter, vSGLT.

Sodium-dependent transporters couple the flow of Na+ ions down their electrochemical potential gradient to the uphill transport of various ligands. Many of these transporters share a common core structure composed of a five-helix inverted repeat and deliver their cargo utilizing an alternating-access mechanism. A detailed characterization of inward-facing conformations of the Na+-dependent sugar transporter from Vibrio parahaemolyticus (vSGLT) has previously been reported, but structural details on additional conformations and on how Na+ and ligand influence the equilibrium between other states remains unknown. Here, double electron-electron resonance spectroscopy, structural modeling, and molecular dynamics are utilized to deduce ligand-dependent equilibria shifts of vSGLT in micelles. In the absence and presence of saturating amounts of Na+, vSGLT favors an inward-facing conformation. Upon binding both Na+ and sugar, the equilibrium shifts toward either an outward-facing or occluded conformation. While Na+ alone does not stabilize the outward-facing state, gating charge calculations together with a kinetic model of transport suggest that the resting negative membrane potential of the cell, absent in detergent-solubilized samples, may stabilize vSGLT in an outward-open conformation where it is poised for binding external sugars. In total, these findings provide insights into ligand-induced conformational selection and delineate the transport cycle of vSGLT.

Protonation state of glutamate 73 regulates the formation of a specific dimeric association of mVDAC1.

The voltage-dependent anion channel (VDAC) is the most abundant protein in the outer mitochondrial membrane and constitutes the primary pathway for the exchange of ions and metabolites between the cytosol and the mitochondria. There is accumulating evidence supporting VDAC's role in mitochondrial metabolic regulation and apoptosis, where VDAC oligomerization has been implicated with these processes. Herein, we report a specific pH-dependent dimerization of murine VDAC1 (mVDAC1) identified by double electron-electron resonance and native mass spectrometry. Intermolecular distances on four singly spin-labeled mVDAC1 mutants were used to generate a model of the low-pH dimer, establishing the presence of residue E73 at the interface. This dimer arrangement is different from any oligomeric state previously described, and it forms as a steep function of pH with an apparent pKa of 7.4. Moreover, the monomer-dimer equilibrium affinity constant was determined using native MS, revealing a nearly eightfold enhancement in dimerization affinity at low pH. Mutation of E73 to either alanine or glutamine severely reduces oligomerization, demonstrating the role of protonated E73 in enhancing dimer formation. Based on these results, and the known importance of E73 in VDAC physiology, VDAC dimerization likely plays a significant role in mitochondrial metabolic regulation and apoptosis in response to cytosolic acidification during cellular stress.

Protons in small spaces: Discrete simulations of vesicle acidification.

The lumenal pH of an organelle is one of its defining characteristics and central to its biological function. Experiments have elucidated many of the key pH regulatory elements and how they vary from compartment-to-compartment, and continuum mathematical models have played an important role in understanding how these elements (proton pumps, counter-ion fluxes, membrane potential, buffering capacity, etc.) work together to achieve specific pH setpoints. While continuum models have proven successful in describing ion regulation at the cellular length scale, it is unknown if they are valid at the subcellular level where volumes are small, ion numbers may fluctuate wildly, and biochemical heterogeneity is large. Here, we create a discrete, stochastic (DS) model of vesicular acidification to answer this question. We used this simplified model to analyze pH measurements of isolated vesicles containing single proton pumps and compared these results to solutions from a continuum, ordinary differential equations (ODE)-based model. Both models predict similar parameter estimates for the mean proton pumping rate, membrane permeability, etc., but, as expected, the ODE model fails to report on the fluctuations in the system. The stochastic model predicts that pH fluctuations decrease during acidification, but noise analysis of single-vesicle data confirms our finding that the experimental noise is dominated by the fluorescent dye, and it reveals no insight into the true noise in the proton fluctuations. Finally, we again use the reduced DS model explore the acidification of large, lysosome-like vesicles to determine how stochastic elements, such as variations in proton-pump copy number and cycling between on and off states, impact the pH setpoint and fluctuations around this setpoint.

Antibacterial membrane attack by a pore-forming intestinal C-type lectin.

Human body-surface epithelia coexist in close association with complex bacterial communities and are protected by a variety of antibacterial proteins. C-type lectins of the RegIII family are bactericidal proteins that limit direct contact between bacteria and the intestinal epithelium and thus promote tolerance to the intestinal microbiota. RegIII lectins recognize their bacterial targets by binding peptidoglycan carbohydrate, but the mechanism by which they kill bacteria is unknown. Here we elucidate the mechanistic basis for RegIII bactericidal activity. We show that human RegIIIα (also known as HIP/PAP) binds membrane phospholipids and kills bacteria by forming a hexameric membrane-permeabilizing oligomeric pore. We derive a three-dimensional model of the RegIIIα pore by docking the RegIIIα crystal structure into a cryo-electron microscopic map of the pore complex, and show that the model accords with experimentally determined properties of the pore. Lipopolysaccharide inhibits RegIIIα pore-forming activity, explaining why RegIIIα is bactericidal for Gram-positive but not Gram-negative bacteria. Our findings identify C-type lectins as mediators of membrane attack in the mucosal immune system, and provide detailed insight into an antibacterial mechanism that promotes mutualism with the resident microbiota.

Harmonizing Experimental Data with Modeling to Predict Membrane Protein Insertion in Yeast.

Membrane proteins must adopt their proper topologies within biological membranes, but achieving the correct topology is compromised by the presence of marginally hydrophobic transmembrane helices (TMHs). In this study, we report on a new model membrane protein in yeast that harbors two TMHs fused to an unstable nucleotide-binding domain. Because the second helix (TMH2) in this reporter has an unfavorable predicted free energy of insertion, we employed established methods to generate variants that alter TMH2 insertion free energy. We first found that altering TMH2 did not significantly affect the extent of protein degradation by the cellular quality control machinery. Next, we correlated predicted insertion free energies from a knowledge-based energy scale with the measured apparent free energies of TMH2 insertion. Although the predicted and apparent insertion energies showed a similar trend, the predicted free-energy changes spanned an unanticipated narrow range. By instead using a physics-based model, we obtained a broader range of free energies that agreed considerably better with the magnitude of the experimentally derived values. Nevertheless, some variants still inserted better in yeast than predicted from energy-based scales. Therefore, molecular dynamics simulations were performed and indicated that the corresponding mutations induced conformational changes within TMH2, which altered the number of stabilizing hydrogen bonds. Together, our results offer insight into the ability of the cellular quality control machinery to recognize conformationally distinct misfolded topomers, provide a model to assess TMH insertion in vivo, and indicate that TMH insertion energy scales may be limited depending on the specific protein and the mutation present.

Atomistic insight into lipid translocation by a TMEM16 scramblase.

The transmembrane protein 16 (TMEM16) family of membrane proteins includes both lipid scramblases and ion channels involved in olfaction, nociception, and blood coagulation. The crystal structure of the fungal Nectria haematococca TMEM16 (nhTMEM16) scramblase suggested a putative mechanism of lipid transport, whereby polar and charged lipid headgroups move through the low-dielectric environment of the membrane by traversing a hydrophilic groove on the membrane-spanning surface of the protein. Here, we use computational methods to explore the membrane-protein interactions involved in lipid scrambling. Fast, continuum membrane-bending calculations reveal a global pattern of charged and hydrophobic surface residues that bends the membrane in a large-amplitude sinusoidal wave, resulting in bilayer thinning across the hydrophilic groove. Atomic simulations uncover two lipid headgroup-interaction sites flanking the groove. The cytoplasmic site nucleates headgroup-dipole stacking interactions that form a chain of lipid molecules that penetrate into the groove. In two instances, a cytoplasmic lipid interdigitates into this chain, crosses the bilayer, and enters the extracellular leaflet, and the reverse process happens twice as well. Continuum membrane-bending analysis carried out on homology models of mammalian homologs shows that these family members also bend the membrane-even those that lack scramblase activity. Sequence alignments show that the lipid-interaction sites are conserved in many family members but less so in those with reduced scrambling ability. Our analysis provides insight into how large-scale membrane bending and protein chemistry facilitate lipid permeation in the TMEM16 family, and we hypothesize that membrane interactions also affect ion permeation.

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.

Entry from the Lipid Bilayer: A Possible Pathway for Inhibition of a Peptide G Protein-Coupled Receptor by a Lipophilic Small Molecule.

The pathways that G protein-coupled receptor (GPCR) ligands follow as they bind to or dissociate from their receptors are largely unknown. Protease-activated receptor-1 (PAR1) is a GPCR activated by intramolecular binding of a tethered agonist peptide that is exposed by thrombin cleavage. By contrast, the PAR1 antagonist vorapaxar is a lipophilic drug that binds in a pocket almost entirely occluded from the extracellular solvent. The binding and dissociation pathway of vorapaxar is unknown. Starting with the crystal structure of vorapaxar bound to PAR1, we performed temperature-accelerated molecular dynamics simulations of ligand dissociation. In the majority of simulations, vorapaxar exited the receptor laterally into the lipid bilayer through openings in the transmembrane helix (TM) bundle. Prior to full dissociation, vorapaxar paused in metastable intermediates stabilized by interactions with the receptor and lipid headgroups. Derivatives of vorapaxar with alkyl chains predicted to extend between TM6 and TM7 into the lipid bilayer inhibited PAR1 with apparent on rates similar to that of the parent compound in cell signaling assays. These data are consistent with vorapaxar binding to PAR1 via a pathway that passes between TM6 and TM7 from the lipid bilayer, in agreement with the most consistent pathway observed by molecular dynamics. While there is some evidence of entry of the ligand into rhodopsin and lipid-activated GPCRs from the cell membrane, our study provides the first such evidence for a peptide-activated GPCR and suggests that metastable intermediates along drug binding and dissociation pathways can be stabilized by specific interactions between lipids and the ligand.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations.

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models.

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

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.