Abstract
The Major Facilitator Superfamily (MFS) transporters are an essential class of secondary active transporters involved in various physiological and pathological processes. The melibiose permease (MelB), which catalyzes the stoichiometric symport of the disaccharide melibiose and monovalent cations (e.g., Na+, H+, or Li+), is a key model for understanding the cation-coupled symport mechanisms. Extensive experimental data has established that positive cooperativity between the cargo melibiose and the coupling cation is central to the symport mechanism. However, the structural and energetic origins of this cooperativity remain unclear at the atomistic level for MelB and most other coupled transporters. Here, leveraging recently resolved structures in inward- and outward-facing conformations, we employed the string method and replica-exchange umbrella sampling simulation techniques to comprehensively map the all-atom free energy landscapes of the Na+-coupled melibiose translocation across the MelB in Salmonella enterica serovar Typhimurium (MelBSt), in comparison with the facilitated melibiose transport in a uniporter mutant. The simulation results unravel asymmetrical free energy profiles of melibiose translocation, which is tightly coupled to protein conformational changes in both the N- and C-terminal domains. Notably, the cytoplasmic release of the melibiose induces the simultaneous opening of an inner gate, resulting in a high-energy state of the system. Periplasmic sugar binding and cytoplasmic melibiose released are dynamically coupled with changes in the internal gating elements along the translocation pathway. The outward-facing sugar-bound state is thermodynamically most stable, while the occluded state is a transient state. The binding of Na+ facilitates melibiose translocation by increasing the melibiose-binding affinity and decreasing the overall free energy barrier and change. The cooperative binding of the two substrates results from the allosteric coupling between their binding sites instead of direct electrostatic interaction. These findings add substantial new atomic-level details into how Na+ binding facilitates melibiose translocation and deepen the fundamental understanding of the molecular basis underlying the symport mechanism of cation-coupled transporters.
Introduction
Major facilitator superfamily (MFS) membrane transporters (MFS) are ubiquitous across all kingdoms of life, accounting for more than 25% of all transmembrane proteins, and exhibit a wide range of functions and substrate specificities. They are critical in various physiological processes, including the uptake of nutrients and drugs, as well as the expulsion of xenobiotics1. Certain human homologs have become promising targets for drug development due to their vital roles in nutrient and drug transport2. For example, the MelB homolog MFSD2A is crucial for the uptake of the essential lipid lysophosphatidylcholine in the brain3. Melibiose permease of Salmonella enterica serovar Typhimurium (MelBSt) is a MFS sugar symporter that facilitates the simultaneous translocation of one galactoside-containing disaccharide (e.g., melibiose) and one cation (e.g., Na+, H+, or Li+) across the membrane in the same direction, with a strict 1:1 stoichiometric ratio 4 (Fig. 1). This symporter is a well-established and useful model for studying the cation-coupled transport mechanisms of MFS transporters1, 4–8. It has been experimentally well-characterized by various biochemical/biophysical techniques in combination with genetic modification9–15 and structural analysis4, 16–19. Additionally, extensive experimental data has been made available for MelB of E. coli20–27. Consistent results on the molecular recognitions for the galactoside and cation and their cooperative binding led to a hypothesis that the cooperative binding of melibiose and the coupling cation is critical for the symport mechanism10, 28. More recently, the correlations between substrate-binding affinities and protein conformational states have been identified through structural and binding analyses19, 29. While the cation-binding affinity appears largely independent of protein conformation, it was hypothesized that the sugar-binding affinity significantly decreases when MelBSt is in the inward-facing (IF) state19, 29.
To date, the structures for two major conformational states of MelBSt have been resolved: the apo- or sugar-bound outward-facing (OF) state16, 17, and the sugar-released, Na+-bound IF state18. However, the molecular mechanisms underlying the coupled symport of the two substrates remain poorly understood. Crucial information regarding the structures, energetics and dynamics of the intermediate states (Fig. 1) during the transitions between the two major conformational states remains scarce. Furthermore, previous experiments have established that the binding of the two substrates is cooperative28, but the structural and energetic basis of this cooperativity remains unclear. Although the regions of the mobile energy barriers for regulating the sugar translocation have been identified19, the global free energy landscape of the entire sugar translocation process— encompassing sugar-binding and releasing events along with conformational transitions between the OF and IF conformations—remains largely unknown. The lack of quantitative characterization of such free energy landscapes has impeded a fundamental understanding of the complex interplays between the shifting of the mobile barriers and the cooperative binding of the two substrates.
Molecular simulations can fully elucidate the structures, dynamics, and thermodynamics of the transport cycle at atomic-level detail. However, calculating the sugar-translocation free energy landscape is a highly challenging task, because it requires sampling the large structural changes in the transporter coupled with substrate binding/unbinding processes. Although numerous computational studies have mapped the free energy landscapes of secondary active transporters and uniporters30–42, none characterized the entire substrate translocation process coupled to global protein conformational transitions in cation-coupled MFS transporters. In addition, to the best of our knowledge, there has been no study on the structural and energetic origins of the cooperative transport of two substrates in cation-coupled MFS transporters.
To this end, we calculated the free energy profiles for the translocation of the melibiose across the WT MelBSt in the Na+-bound and -unbound states, as well as the uniport D59C mutant. Building upon the recently resolved structures for both OF and IF states17, 19, the string method43, 44 was employed to identify the minimum free energy pathway (MFEP) for the translocation of the melibiose molecule from the periplasmic to cytoplasmic sides, accompanied by the OF to IF conformational transition. Extensive replica-exchange umbrella sampling (REUS) simulations were performed along the MFEP to quantify the free energy surfaces of melibiose translocation with or without a bound Na+. The simulations correspond to the experimentally observed melibiose exchange4, 9, 16, 17, 45 driven by melibiose concentration gradients without cation transduction (Fig. 1). Explicit characterization of the thermodynamics underlying the entire melibiose translocation process reveals the structural and energetic underpinnings of how the sugar translocation across MelBSt is facilitated by the binding of the coupling cation Na+. To the best of our knowledge, this work is the first time that the free energy landscape dictating the functional cycle of a cation-coupled MFS symporter is characterized at full atomic-level detail, thus significantly deepening our understanding of cation-coupled transporters in general.
Results
Melibiose translocation coupled to protein conformational changes
The Fig. 2 depicts the free-energy profile (or potential of mean force, (PMF)) for the melibiose translocation from the periplasmic to the cytoplasmic sides of the membrane through the WT MelBSt with a Na+ bound at the cation-binding pocket. The free energy changes (ΔG′s) and barriers (ΔG‡′s) of key steps during the translocation are summarized in Table 1. The translocation process starts with a melibiose molecule in the periplasmic bulk and a MelBSt in its outward-facing sugar-free state (OFF) (image ID 0). Binding of melibiose from the periplasmic side leads to the outward-facing sugar-bound state (OFB) (image ID 13, Fig. S1), which exhibited a ΔG of −5.0 kcal/mol (Table 1, ΔGOF binding). This corresponds to a 5.0 kcal/mol melibiose binding affinity when the Na+-bound protein is in the OF conformation. Notably, this state has the lowest free energy in the entire translocation process and thus is thermodynamically most stable. Then, the protein conformational changes took the system from the OFB to the inward-facing sugar-bound state (IFB) (image ID 19), overcoming a barrier of ∼6.8 kcal/mol (Table 1, ΔG‡OF to IF) required for passing through the occluded transition state (OCB) during the conformational transition. The OCB state is a transient state, and the IFB state is higher in free energy than the OFB state by ∼6.4 kcal/mol (Fig. 2, Table 1). Then, the melibiose was released to the cytoplasmic side of the bulk solution, reaching the inward-facing sugar-free state (IFB, image ID 32), leaving behind a Na+-bound, sugar-free MelBSt. This cytoplasmic melibiose release process only needs to overcome a small ΔG‡of ∼2.8 kcal/mol (Table 1, ΔG‡OF to IF) and has a small ΔG of ∼2.3 kcal/mol, indicating a lower melibiose-binding affinity of the IF than the OF conformations (∼5.0 kcal/mol). The entire transport process has a net ΔG of ∼3.7 kcal/mol (Table 1 and Fig. 2) and an overall ΔG‡ of ∼8.4 kcal/mol (Table 1 and Fig. 2).
To characterize the coupling between the melibiose translocations and protein conformational transitions, the free energy landscape was projected onto 2D free energy surfaces (FES) spanned by different combinations of collective variables (CVs) using Eq S4 (Fig. 3). In Fig. 3A, the 2D FES is spanned by the first principle component (PC1) of the backbone Cα atoms of all 12 transmembrane helices and the melibiose transport coordinate Z (see Method in SI for definition). As Z increases from −30 Å to +30 Å, the melibiose molecule moves through MelBSt from the periplasmic to the cytoplasmic sides of the membrane. As the PC1 gradually increases from ∼-3.5 Å to +3 Å, the transmembrane helices transit from the OF to the IF conformational states.
As shown in Fig. 3A, the periplasmic sugar-binding process is correlated with dynamic fluctuations of the transmembrane helices, as exhibited by the shifts in the PC1 values across the multiple minima encountered in the OFF→OFB process. For example, from OFF (Z = −30 Å, PC1 = −3.2 Å) to OFB (Z = −3 Å, PC1 = −2.5 Å), the PC1 values fluctuate back and forth and has a net increase of ∼0.7 Å after reaching the OFB minimum. These fluctuations in PC1 correspond to non-negligible conformational fluctuations in the overall transmembrane protein backbones during the periplasmic-binding process. The binding-induced conformational change is more evident when the FES is projected to the 2D plane spanned by the translocation coordinate and the interhelical distance between helix I and VII near the periplasmic side of the membrane (Fig. 3B). There is a decrease in this interhelical distance from ∼22 Å at OFF to ∼19 Å at OFB as the melibiose approaches the binding site. Thus, the periplasmic sugar-binding slightly closes the periplasmic side of the sugar translocation pathway, preparing the protein for the subsequent conformational changes.
After the periplasmic sugar-binding event, the OFB→IFB process begins, largely increasing the PC1 from approximately −2.5 Å to +2 Å (Fig. 3A). This process features global protein conformational changes critical for the entire transport process. In this step, the sugar translocation pathway is fully closed on the periplasmic side and opened on the cytoplasmic side, preparing the protein for releasing sugar into the cytoplasm. A large decrease in multiple periplasmic interhelical distances is observed, such as between helices I and VII (from ∼19 Å to ∼14 Å) (Fig. 3B), as well as an increase in cytoplasmic interhelical distances, such as between the helices IV and X (from ∼13 Å to ∼20 Å) (Fig. 3C).
After the system reaches the IFB state, the melibiose release to the cytoplasmic bulk begins (IFB⟶IFF). The process goes through an intermediate state IFC, where the melibiose passes through the most constricted region of the entire pathway contributed by Val145, Val346 (helix X), and Tyr369 (helix XI) (Fig. 3A). From IFB to IFC, the cytoplasmic inter-helices distances further increase. For example, the distance between helices IV and X further increases to ∼23 Å near Z = + 11 Å to prepare the space for the melibiose to pass through (Fig. 3A). After this, the melibiose is eventually released to the cytoplasmic bulk and the interhelical distance between helices IV and X slightly decreases back to below 20 Å. The diagonal nature of the pathway linking IFB, IFC, and IFF (Fig. 3C) indicates that the motion of the substrate and the expansion/shrinkage of the cytoplasmic constricted region are tightly coupled. The finding is significant because it reveals that bound melibiose does not passively wait for the cytoplasmic path to widen, as previously speculated. Instead, it actively induces protein conformational changes to move out. This observation redefines the alternation-access model, demonstrating that both the periplasmic binding and cytoplasmic release of melibiose are tightly coupled with protein conformational changes.
Changes in the pore radius profile during melibiose translocation
The coupling between the conformational change and the melibiose translocation is also evident through analyzing the radius profiles of the internal cavities that form the sugar translocation path in different intermediate states. For example, we first monitored the radius of the cavity near the periplasmic bulk (Z = −25 Å) during the entire sugar translocation process. This radius slightly decreases during periplasmic binding (OFF→OFB, red vs. orange curves in Fig. 4A). It then largely decreases by more than 3 Å during conformational transition (OFB→IFB, orange vs. green curves in Fig. 4A). Following this, it remains almost the same during cytoplasmic release (IFB→IFF, green vs. blue curves in Fig. 4A). Monitoring the radius of the cavity near the cytoplasmic bulk (Z = + 25 Å), it remains small below 2 Å during the periplasmic binding, and then increases by approximately 3 Å following the OFB→IFB transition. Thus, in general, the OFF→IFF process is accompanied by the widening and shrinking of the periplasmic and cytoplasmic terminals of the path in a reciprocal manner.
The radius profile also revealed two constricted regions surrounding the bound sugar molecular, which functions as gates on both the periplasmic and cytoplasmic sides of the sugar translocation path. The periplasmic constriction region is located near Z = −15 to −10 Å and formed by the cavity-lining residues, particularly the Tyr26 and Met27 on the kink of the helix I, and Asn248 and Asn251 on the helix VII (Fig. 4A, right panel). Clearly, the helix I kink moves towards the helix VII positioning in the middle of the melibiose-accessing path upon melibiose binding. Tyr26 frequently interacts with the melibiose as it passes through and defines the periplasmic edge of the sugar-binding pocket, and its extended bulky sidechain partially occludes the bound sugar from the periplasmic side. This constricted region is further narrowed when MelBSt changes to the occluded transition state. Thus, the Tyr26 functions as a key gating residue. Experimentally, mutation of the Tyr26 to Cys26 residue resulted in the loss of active transport of the MelBSt, highlighting the importance of this residue in the functional cycle13. The cytoplasmic constriction region is located near Z = +11 Å, which is formed by the cavity-lining residues Val145, Val346 and Tyr369, as discussed above. When the periplasmic side of the transporter is open, the N- and C-terminal domains on the cytoplasmic side form a tightly connected salt-bridge network, as described previously13, 17. These ionic interactions majorly contribute to forming the inner barrier of cytoplasmic sugar release. As the melibiose passes through this region, the radius in this region expands to accommodate its motion (Fig. 4C), in line with the above-mentioned observation from the FES (Fig. 3C).
Another important observation is that periplasmic binding reduces the radius of the sugar-binding pocket (Fig. 4B), reaffirming the binding-induced conformational change discussed above (Fig. 3). The subsequent OF to IF conformational changes (OFB →IFB) increase the radius of the sugar-binding cavity near Z = −2.5 Å (Fig. 4B). The expansion of the binding pocket eventually facilitates the release of melibiose to the cytoplasmic side. This is consistent with the results from the free energy profile calculation. The simulation results thus corroborate the previous experimental result19 that the outward-facing conformation has a higher sugar-binding affinity than the inward-facing (Fig. 2-3) and explain its structural and energetic origin with atomic-level detail.
Cooperative motion of N- and C-terminal domains
Multiple helices in both the N- and C-terminal domains change their relative orientation with respect to the normal of the membrane plane during the melibiose translocation. To quantify this observation, the free energy of the system is projected to 2D surfaces spanned by two collective variables (CVs), which measure the directional tilt angle of a helix with respect to the membrane plane normal (see “Method” for definition) and melibiose translocation. Based on the topographies of the 2D FES’s (Fig. S4 B-G), it is obvious that in the N-terminal domain, all helices decrease their tilting angles by 10-15 degrees as Z moves from −35 to 35 Å. In contrast, in the C-terminal domain, all helices increase their tilt angles by 10-20 degrees as Z moves from −35 to 35 Å (Fig. S4 H-M). The simultaneous tilting of all helices suggests that both domains must reorient with respect to the membrane normal to complete the global protein conformational transition. This transition exhibits a symmetric-like movement of the two pseudo-symmetric helix bundles. Notably, this finding, which cannot be derived from experimental structures alone due to the lack of a membrane surface, helps resolve a longstanding debate regarding the alternating-access model.
Melibiose translocation energetically coupled to Na+ binding
Next, we characterize the PMF of melibiose translocation when the MelBSt is in the Na+-unbound state (Fig. 5, blue curve; Table 1). The comparison between the PMFs of Na+-unbound and the Na+-bound states is essential for understanding how the Na+ binding cooperates with the thermodynamics and kinetics of melibiose translocation. When the Na+ is unbound, the overall net ΔG and overall barrier ΔG‡are ∼11.8 kcal/mol and ∼12.2 kcal/mol, respectively (Table 1). These values are greater than the Na+-bound state by 8.1 kcal/mol and 3.8 kcal/mol, respectively, indicating that the melibiose translocation is thermodynamically and kinetically less favorable if the Na+ is unbound. Several factors contribute to the more endergonic and slower translocation process in the absence of bound Na+. First, periplasmic melibiose binding (image IDs from 0 to 15) is endergonic with ΔGOF binding ∼ 1.4 kcal/mol, in contrast to the exergonic process in the Na+ bound state (ΔGOF binding ∼ −5.0 kcal/mol). Second, the endergonic protein conformational transition (image IDs from 15 to 21) exhibits ΔGOF to IF ∼ 7.0 kcal/mol and ΔG‡OF to IF ∼ 8.0 kcal/mol, which are both greater than the Na+ bound state (6.4 kcal/mol and 6.8 kcal/mol, respectively). Third, the endergonic cytoplasmic unbinding event (images 21 to 32) exhibits ΔGIF binding ∼ 3.4 kcal/mol and ΔG‡IF unbinding ∼ 3.8 kcal/mol, again greater than the Na+ bound state (2.3 kcal/mol and 2.8 kcal/mol, respectively).
The well-characterized D59C mutation results in a melibiose uniporter due to the loss of the carboxyl group on the Asp59, which is the only residue critical for binding of all three types of coupled cations (Na+, H+, and Li+)12, 18, 28. This mutant thus loses cation-coupled active transport activity but can facilitate the transport of melibiose downhill its concentration gradient17 and the exchange of melibiose across the membrane 16. This unique feature makes the mutant a good model system for understanding the Na+/melibiose symport mechanism.
The melibiose translocation across the D59C mutant has an overall net ΔG of ∼ 8.3 kcal/mol and overall free energy barrier ΔG‡ of 10.6 kcal/mol (Fig. 5, green curve; Table 1). They are both greater than the WT Na+ bound system (3.7 kcal/mol and 8.4 kcal/mol, respectively), making the melibiose translocation thermodynamically and kinetically less favorable. Several factors contribute to the more endergonic and slower melibiose translocation in the mutant. The periplasmic sugar binding (image IDs from 0 to 13) has a ΔGOF binding of ∼ −2.0 kcal/mol, indicating that the mutation diminishes the binding affinity of the OF state by ∼3.0 kcal/mol as compared to the WT Na+-bound system. Unlike the WT Na+-unbound system, the periplasmic binding in the D59C mutant is exergonic despite its incapability of binding cations, which is consistent with previous experiential measurements17. After the OFB state is reached, the OFB⟶IFB transition together with the cytoplasmic sugar release (IFB⟶IFF) (image IDs from 13 to 32) results in a combined high barrier of ∼10.6 kcal/mol (ΔGOF to IF + ΔG‡IF unbinding). This rate-limiting barrier in the mutant is ∼2.2 kcal/mol higher than the WT Na+ bound system. Notably, the mutant features ∼3 kcal/mol higher barrier for the cytoplasmic sugar release, which is likely caused by the different orientations of the residues lining the cytoplasmic constricted region.
In both Na+ unbound systems (WT and D59C), as expected, melibiose translocation is still coupled with protein conformational changes, but to a lesser extent than in the WT, Na+ bound system (Figs. S5-S10). In both Na+ unbound systems, periplasmic binding and cytoplasmic release are coupled with changes in the overall conformation of the protein and specific interhelical distances (Figs. S5-S6), as well as pore radius of the sugar-binding pocket and cytoplasmic constricted region (Figs. S7-S8). These changes are smaller in scale compared to the WT though. The directional tilt angles of helices in the N- and C-terminal domains follow the same trend as the WT Na+ bound system (Figs. S9-S10).
Overall, the removal of Na+ from the cation-binding site, either due to low Na+ concentrations in bulk or loss of Na+ binding site, can largely reduce melibiose binding affinity, as has been well-documented by a variety of experimental tests in varied methods9, 12, 16–18, 28. This is one of the key factors contributing to the increased overall ΔG’s and ΔG‡’s for the melibiose translocation in the two Na+ unbound systems. These results explicitly demonstrate that the free energy profiles explicitly demonstrate that the co-transport of Na+ and melibiose is a direct result of binding cooperativity between the two substrates.
Allosteric coupling between Na+ and melibiose binding
To understand the different melibiose binding affinities in the WT Na+ bound and unbound systems, we performed a total of 500 ns unbiased MD simulations in the OFB state for each system. The interaction energies between the melibiose and MelBSt are compared between these two systems (Fig. S11). The interaction energy is decomposed into contributions from electrostatic and van der Waals components (Fig. S11 A&B, green and red bars, respectively). Importantly, the direct interaction energy between the bound Na+ and melibiose in the WT Na+ bound system was also evaluated.
Upon the unbinding of Na+ (Fig. S11A), the protein-melibiose interaction energy shifted from −113 ± 2 to −105.4 ± 0.5 kcal/mol in WT MelBSt, consistent with the reduction in the binding affinity of the melibiose. Such a reduction largely arises from the electrostatic component, which changes from −96 ± 2 to −88.0 ± 0.4 kcal/mol. Notably, the direct interaction energy between the bound Na+ and melibiose (∼ 0.2 ± 0.5 kcal/mol) in the WT Na+ bound state is negligible. This finding is important in that the decrease in the melibiose binding affinity upon Na+ unbinding is not due to losing the direct Na+-melibiose interaction but due to a weakened protein-melibiose interaction. The decreased melibiose-binding affinity in the absence of Na+ is mainly attributed to the alteration of local electrostatic interactions between melibiose and its nearby residues (Fig. S11B).
The WT Na+ bound system features a narrower distribution of the local melibiose-protein interaction energy that peaked at a lower energy range than the WT Na+ unbound system (Fig. S12A). This implies that the Na+ unbinding alters the binding cavity and disfavors melibiose binding. To further probe the structural origin of this effect, we analyzed the hydrogen-bonding interactions between the melibiose and protein (Table S1). Na+ unbinding from the WT MelBSt led to the loss of 0.9 ± 0.2 hydrogen bonds between the melibiose and protein. In addition, the hydrogen bonds with three charged residues Asp124, Asp19, and Arg149, which contribute most to the protein-melibiose hydrogen bonds in the OFB state, are affected by the Na+ unbinding the most. Thus, the reduction of hydrogen bonds between melibiose with the protein contributes to decreased binding affinity in the absence of a bound Na+ ion.
The change in the hydrogen-bonding interaction due to Na+ unbinding likely results from the allosteric coupling between the cation-binding and melibiose-binding sites. The Asp124 and Tyr120 residues in helix IV are essential in the coupling pathway since this helix also contributes the Thr121 residue to the Na+-binding site. Without a bound Na+ ion, the distance between the Asp55 and Asp124 is elongated (Fig. S12B). The same trend is also observed for the distances between Asp124 and Asp19 (Fig. S12C), between Tyr120 and Asp55 (Fig. S12D), and between Tyr120 and Asp19 (Fig. S12E). Furthermore, the distance between Asp55 and Asp19, i.e., the two key residues in the cation- and sugar-binding sites, respectively, is increased upon Na+ unbinding (Fig. S12F). The correlations between the changes in these pairwise distances suggest that both Asp124 and Tyr120 residues propagate structural changes from the Na+-binding site to the sugar-binding site upon Na+ unbinding, leading to a decrease in melibiose affinity. The unbinding of Na+ weakens the hydrogen bond between Asp19 and the 3-hydroxyl group of melibiose, as indicated by their elongated separation distance (Fig. S12G). This observation aligns with previous experimental studies13, 16, 24, 27, 46: the D124C or Y120C mutants can bind melibiose but lose Na+-coupled melibiose active transport16. Our simulations thus provide an atomistic-level explanation for these experimental observations in these decoupling mutants.
Discussion
Calculating the free energy landscape that describes the protein conformational transitions associated with substrate translocation through transporters is a challenging yet essential task for understanding their functional mechanisms. In this study, building upon two resolved ligand-bound structures representing the inward-19 and outward-facing conformational states17, 18, we uncovered a wealth of new mechanistic information about the functional cycle of MelB. This was achieved through extensive free energy calculations combined with state-of-the-art reaction path-finding techniques. Notably, the use of two experimental structures as initial guesses for the two endpoints of the minimum free energy pathway (outward-facing sugar-free and inward-facing sugar-free states) enhances the reliability of the string method (Fig. S3). This is an improvement over previous applications of this approach to other transporters, where only one major conformational state (outward-facing, inward-facing, or occluded) was experimentally available, and the other major conformations had to be generated by biased molecular dynamics simulations31, 33, 35. Our simulation results are consistent with a large body of experimental data collected in the past decades via varied biochemical and biophysical data and structural analyses9, 10, 12, 13, 18, 19, 28.
First, all simulation data consistently showed that Na+ binding increases the melibiose-binding affinity of WT MelBSt, in agreement with previous experimental measurements9, 28. Second, the free-energy landscape indicates that the inward-facing states are thermodynamically less stable than the outward-facing states, regardless of the binding of Na+ and melibiose, providing the energetic data supporting the conclusion drawn from structural analysis19. Third, the overall free energy landscape reveals that the outward-facing conformation of the WT has the highest sugar-binding affinity, and the OFB state is the thermodynamically most stable state for the entire translocation process (Table 1 & Fig. 2). This is a piece of direct evidence supporting that the IF conformational state has a lower sugar-binding affinity and the experimentally measured binding affinity is primarily related to the OF conformational state19.
During the melibiose translocation, only a modest rate-limiting barrier is encountered in the OFB→OFF process (∼5 kcal/mol, Fig. 2). The protein conformational transition and cytoplasmic release is a cooperative process that increases the free energy of the system. The opening of the cytoplasmic gate is narrow and is induced by the motion of the sugar, so the cytoplasmic sugar-release process can be appropriately described as sugar squeezing through the inner barrier. In addition, the lack of a thermodynamically stable occluded state (OCB) facilitates the OFB→IFB transition. All these features make MelBSt a highly effective transporter.
The free energy landscapes reveal the coupling mechanisms between the two substrates and the transporter, providing new insights into experimental measurements of MelBSt and its mutants10, 17, 19, 28. Without Na+, both WT and D59C mutants exhibit higher ΔG‡and ΔG than the WT Na+-bound system, making melibiose translocation less favorable. This suggests that the coupling between sugar and cation is primarily due to energetic coupling at the binding step. Na+binding increases melibiose affinity, lowers the free energy barrier and change, and facilitates the protein conformational transition, thus accelerating melibiose translocation. The molecular basis underlying the difference in the free energy landscapes involves strengthened hydrogen bonds between melibiose and its binding residues (Asp19, Asp124, Arg149) due to Na+ binding at the cation-binding pocket (Asp55, Asp59, Thr121, Asn58), with allosteric coupling through Asp124 and Tyr120 on helix IV, as proposed sed on experimental data16, 24, 27, 46. The interaction between the two binding sites underpins the energetic and kinetic coupling necessary for the cooperative transport of Na+ and melibiose.
The WT MelBSt mediates both the active transport and facilitated diffusion of melibiose 4. In contrast, the D59C mutation eliminates the active transport mode of WT but retains its melibiose uniport activity16, 17, similar to the E325A mutant of LacY7, 47. Although such a phenomenon is observed in the D59C MelB and has been an important clue for identifying the cation-binding site of cation-coupled secondary transporters (7, 47, 48), the molecular origin of why the mutant retains uniport activity remains elusive to date. Here, the comparison between the free energy landscapes of the WT and mutant offers a microscopic-level explanation. For the D59C mutant, the overall ΔG‡ of melibiose translocation is modestly higher than WT, Na+ bound MelBSt by only 2.2 kcal/mol (Table 1). Thus, this mutation does not kinetically block the inward sugar translocation entirely, and it can even facilitate the outward translocation process by reducing the periplasmic sugar-releasing barrier as a result of reduced binding affinity in the OF state (Fig. 5). Therefore, the bi-directional exchange and gradient-downhill transport of melibiose can still be observed so long as the melibiose concentration is high enough on at least one side of the membrane. Notably, the coupling between the melibiose translocation and conformational change in this uniporter mutant remains (Figs. S6, S8 & S10). From an evolutional point of view, it is likely that the introduction of a cation site into a uniporter can accelerate the substrate translocation and increase the substrate-binding affinity. This symport function evolved from a uniporter can better serve cellular needs by harvesting substrate from scarce conditions.
In summary, our all-atom free energy landscapes of MelBSt provided numerous critical, novel and fundamental insights into the structural and energetic origins of the coupling mechanism essential in cation-coupled MFS symporters: (1) the substrate translocation is tightly coupled to global protein conformational changes where all transmembrane helices are reoriented with respect to the membrane normal, (2) periplasmic binding is coupled with a partial closure of the periplasmic gate and shrinkage of sugar-binding site, (3) the outward-facing sugar-bound state is thermodynamically most stable during the entire translocation process, (4) the occluded state is a transient state, (5) the melibiose cytoplasmic release is coupled to the temporary expansion of the cytoplasmic gate, and (6) the allosteric coupling between the cation- and sugar-binding sites leads to cooperative binding of Na+ and melibiose, which facilitates sugar translocation thermodynamically and kinetically. For the first time, our free energy profiles explicitly demonstrate that the co-transport of Na+ and melibiose is a direct result of binding cooperativity between the two substrates. Notably, although such cooperativity has been observed experimentally on a macroscopic level, our free-energy simulations comprehensively characterized this core mechanism of symport with atomic-level detail, and provide explicit, crucial evidence for energetic coupling between the two co-transported solutes across all intermediate states visited during the entire melibiose translocation process (Figs. 1-2 & 5).
Methods
A brief summary of the methods is provided here, with full details included in the SI.
System setup
Three simulation systems were prepared to calculate the free energy landscapes of melibiose translocation coupled with protein conformational transition: (1) WT MelBSt with Na+ bound, (2) WT MelBSt without Na+, and (3) D59C mutant without Na+. The OF and IF conformational states were derived from crystal (PDB code 7L1617) and CryoEM (PDB code 8T6019) structures, respectively. Key steps in the structural preparation included mutating residues to match the desired states, positioning melibiose in the binding site, adding a lipid bilayer, and neutralizing with 0.15M NaCl. Each system consisted of ∼130,000 atoms in a simulation box of ∼100 ×100 ×130 ų. All systems were constructed using the CharmmGUI web interface48.
Preparation for string method simulations
Initial string images for string method simulations were prepared for all three systems. First, for melibiose binding/unbinding processes in OF and IF states, the system was relaxed with geometry optimization and then equilibrated with backbone restraints. The bound melibiose was gradually pulled out of the protein (cytoplasmic unbinding for the IF state and periplasmic unbinding for the OF state) using a series of harmonic potentials, creating intermediate snapshots for substrate release into both sides of the membrane. This provided initial string images for OFF→OFB and IFB→IFF transitions. Next, the geodesic interpolation algorithm49 generated 10 intermediate images of the transmembrane backbone for the OFB→IFB conformational transition. These backbone structures guided restrained simulations that dynamically relaxed all atoms in the images. All simulations used the CHARMM36m50–56 and TIP3P force fields57, and were performed using the NAMD software package 58.
String method simulations
The String Method with Swarms of Trajectories (SMwST)43, 44 was employed to identify the minimum free energy pathway (MFEP) for melibiose translocation in all three systems. This method has been successfully applied to other types of secondary active transporters31, 33, 35 and protein complexes59. The initial string consisted of 33 images, representing the entire OFF ↔ IFF transition. These images were projected into a 13-dimensional space defined by CVs related to melibiose translocation and protein conformational changes.
The SMwST simulations involved 500 iterations, where each image was equilibrated and propagated in each iteration. The string was updated iteratively until convergence was achieved, as monitored by RMSD in the 13-D space (Fig. S2). The endpoints of the final string corresponded closely to the experimental OF and IF structures (Fig. S3).
Replica-exchange umbrella sampling (REUS) simulations
The last iteration of the string from the SMwST simulation was used as the MFEP for performing REUS simulations60, 61. The 33 image centers were interpolated into equidistant 13-D window centers for 33 umbrella windows. A harmonic potential (0.5-1 kcal/mol) acted on each of the 13 CVs per window. Initial conditions for each window were taken from the images in the last SMwST iteration. Replica exchanges among 10 neighboring windows were attempted every 10 ps on a rotation basis59. The REUS simulations ran for 90 ns per window, with the first 10 ns discarded as equilibration, resulting in 7.9 μs of total sampling time for all three systems.
Analysis
The REUS simulations were unbiased using a generalized version33 of the weighted histogram analysis method62, 63 (Eqs. S1-S2). The free energies of each image were corrected using Eq S3. The 2D free energy surfaces were constructed using Eq. S4. The MFEP connecting the OFF and IFF states on the 2D free energy surfaces (FES) was traced using the zero-temperature string method. The HOLE program64 was used to analyze the pore radius profiles for all snapshots sampled by the REUS simulations.
Data Availability Statement
The data supporting this article have been included in the Supplementary Information.
Acknowledgements
This work was supported by the National Institutes of Health Grants R35GM150780 to R.L. and R35GM153222 to L.G. The researchers used GPU computing facilities provided by the High-Performance Computing Center at Texas Tech University. The authors also acknowledge the helpful insights and inputs provided by Prof. Abhishek Singharoy at Arizona State University, and the initial simulation system setup by Mr. Amirhossein Bakhtiiari at Texas Tech University.
Additional information
Author contributions
The manuscript was written through the contributions of all authors. Ruibin Liang designed the research, performed the simulations, analyzed the data, and wrote the manuscript. Lan Guan designed the research and revised the manuscript. All authors have given approval to the final version of the manuscript.
Conflict of interest
The authors declare no competing financial interests.
All-atom free energy landscapes of the entire melibiose translocation process in wild-type MelBSt and its uniport D59C mutant elucidate the structural and energetic basis of the positive cooperativity between melibiose and its driving cation. The new insights significantly deepen our understanding of the molecular basis underlying cation-coupled transport mechanisms in MFS symporters.
References
- 1.The major facilitator superfamilyJ Mol Microbiol Biotechnol 1:257–279
- 2.SLC transporters as therapeutic targets: emerging opportunitiesNat Rev Drug Discov 14:543–560
- 3.Mfsd2a is a transporter for the essential omega-3 fatty acid docosahexaenoic acidNature 509:503–506
- 4.Na(+)/Melibiose Membrane Transport ProteinMelB.Journal https://doi.org/10.1007/978-3-642-35943-9_10082-1
- 5.Cation coupling to melibiose transport in Salmonella typhimuriumJ. Bacteriol 144:192–199
- 6.Cation and sugar selectivity determinants in a novel family of transport proteinsMol. Microbiol 19:911–922
- 7.Lessons from lactose permeaseAnnu. Rev. Biophys. Biomol. Struct 35:67–91
- 8.Structural Biology of the Major Facilitator Superfamily TransportersAnnu Rev Biophys 44:257–283
- 9.Mechanism of melibiose/cation symport of the melibiose permease of Salmonella typhimuriumJ. Biol. Chem 286:6367–6374
- 10.Cooperative binding ensures the obligatory melibiose/Na+ cotransport in MelBJ. Gen. Physiol 153
- 11.Reduced Na(+) affinity increases turnover of Salmonella enterica serovar Typhimurium MelBJ. Bacteriol 194:5538–5544
- 12.Molecular Basis for the Cation Selectivity of Salmonella typhimurium Melibiose PermeaseJ. Mol. Biol 434
- 13.Complete cysteine-scanning mutagenesis of the Salmonella typhimurium melibiose permeaseJ. Biol. Chem 297
- 14.Structural and functional characterization of protein-lipid interactions of the Salmonella typhimurium melibiose transporter MelBBMC Biol 16
- 15.Substrate-binding guides individual melibiose permeases MelB to structurally soften and to destabilize cytoplasmic middle-loop C3Structure 31:58–67
- 16.Structure-based mechanism for Na(+)/melibiose symport by MelBNat. Commun 5
- 17.X-ray crystallography reveals molecular recognition mechanism for sugar binding in a melibiose transporter MelBCommun Biol 4
- 18.Distinct roles of the major binding residues in the cation-binding pocket of the melibiose transporter MelBJ. Biol. Chem https://doi.org/10.1016/j.jbc.2024.107427
- 19.Mobile barrier mechanisms for Na(+)-coupled symport in an MFS sugar transporterElife 12
- 20.Cation specificity for sugar substrates of the melibiose carrier in Escherichia coliBiochim Biophys Acta 904:191–200
- 21.Sodium-substrate cotransport in bacteriaBiochim Biophys Acta 1505:121–130
- 22.Sugar binding properties of the melibiose permease in Escherichia coli membrane vesicles. Effects of Na(+) and H(+) concentrationsJ. Biol. Chem 261:6893–6899
- 23.Mutagenesis of acidic residues in putative membrane-spanning segments of the melibiose permease of Escherichia coli. I. Effect on Na(+)-dependent transport and binding propertiesJ. Biol. Chem 268:3209–3215
- 24.Mutagenesis of acidic residues in putative membrane-spanning segments of the melibiose permease of Escherichia coli. II. Effect on cationic selectivity and coupling propertiesJ. Biol. Chem 268:3216–3221
- 25.The melibiose permease of Escherichia coli: importance of the NH2-terminal domains for cation recognition by the Na+/sugar cotransporterSoc. Gen. Physiol. Ser 48:213–227
- 26.Evidence for intraprotein charge transfer during the transport activity of the melibiose permease from Escherichia coliBiochemistry 40:13744–13752
- 27.Structural insights into the activation mechanism of melibiose permease by sodium bindingProc Natl Acad Sci U S A 107:22078–22083
- 28.Thermodynamic cooperativity of cosubstrate binding and cation selectivity of Salmonella typhimurium MelBJ. Gen. Physiol 149:1029–1039
- 29.In vivo and in vitro characterizations of melibiose permease (MelB) conformation-dependent nanobodies reveal sugar-binding mechanismsJ. Biol. Chem 299
- 30.How do antiporters exchange substrates across the cell membrane? An atomic-level description of the complete exchange cycle in NarKStructure 29:922–933
- 31.Free Energy Landscape for the Entire Transport Cycle of Triose-Phosphate/Phosphate TranslocatorStructure 26:1284–1296
- 32.Simulating the function of sodium/proton antiportersProc. Natl. Acad. Sci. U. S. A 112:12378–12383
- 33.Atomic-level characterization of transport cycle thermodynamics in the glycerol-3-phosphate:phosphate antiporterNat. Commun 6
- 34.Molecular Basis of the Glucose Transport Mechanism in PlantsACS Central Science 5:1085–1096
- 35.Molecular determinants for the thermodynamic and functional divergence of uniporter GLUT1 and proton symporter XylEPLoS Comput. Biol 13
- 36.Distinct Substrate Transport Mechanism Identified in Homologous Sugar TransportersJ. Phys. Chem. B 123:8411–8418
- 37.Free Energy Landscape of the Complete Transport Cycle in a Key Bacterial TransporterACS Central Science 4:1146–1154
- 38.Uptake dynamics in the Lactose permease (LacY) membrane protein transporterSci. Rep 8
- 39.The substrate import mechanism of the human serotonin transporterBiophys. J 121:715–730
- 40.Sugar Transport across Lactose Permease Probed by Steered Molecular DynamicsBiophys. J 93:92–102
- 41.Locking Two Rigid-body Bundles in an Outward-Facing Conformation: The Ion-coupling Mechanism in a LeuT-fold TransporterSci. Rep 9
- 42.Ligand coupling mechanism of the human serotonin transporter differentiates substrates from inhibitorsNat. Commun 15
- 43.A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and PitfallsJ. Chem. Theory Comput 18:1406–1422
- 44.Finding Transition Pathways Using the String Method with Swarms of TrajectoriesJ. Phys. Chem. B 112:3432–3440
- 45.Role of Gly117 in the cation/melibiose symport of MelB of Salmonella typhimuriumBiochemistry 51:2950–2957
- 46.Mutation of polar and charged residues in the hydrophobic NH2-terminal domains of the melibiose permease of Escherichia coliJ. Biol. Chem 269:24883–24889
- 47.lac permease of Escherichia coli: histidine-322 and glutamic acid-325 may be components of a charge-relay systemBiochemistry 25:4486–4488
- 48.CHARMM-GUI: A web-based graphical user interface for CHARMMJ. Comput. Chem 29:1859–1865
- 49.Geodesic interpolation for reaction pathwaysJ. Chem. Phys 150
- 50.CHARMM36m: an improved force field for folded and intrinsically disordered proteinsNat. Methods 14:71–73
- 51.Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone phi, psi and Side-Chain chi(1) and chi(2) Dihedral AnglesJ. Chem. Theory Comput 8:3257–3273
- 52.Improved Treatment of the Protein Backbone in Empirical Force FieldsJ. Am. Chem. Soc 126:698–699
- 53.All-atom empirical potential for molecular modeling and dynamics studies of proteinsJ. Phys. Chem. B 102:3586–3616
- 54.Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid TypesJ. Phys. Chem. B 114:7830–7843
- 55.Additive empirical force field for hexopyranose monosaccharidesJ. Comput. Chem 29:2543–2564
- 56.CHARMM Additive All-Atom Force Field for Glycosidic Linkages between HexopyranosesJ. Chem. Theory Comput 5:2353–2370
- 57.Comparison of simple potential functions for simulating liquid waterJ. Chem. Phys 79:926–935
- 58.Scalable molecular dynamics on CPU and GPU architectures with NAMDJ. Chem. Phys 153
- 59.Chemomechanical Coupling in Hexameric Protein–Protein Interfaces Harnesses Energy within V-Type ATPasesJ. Am. Chem. Soc 139:293–310
- 60.Replica-exchange molecular dynamics method for protein foldingChem. Phys. Lett 314:141–151
- 61.Multidimensional replica-exchange method for free-energy calculationsJ. Chem. Phys 113:6042–6051
- 62.The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules. 1. The MethodJ. Comput. Chem. 13:1011–1021
- 63.WHAM: the weighted histogram analysis methodGrossfield Lab
- 64.The pore dimensions of gramicidin ABiophys. J 65:2455–2460
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Copyright
© 2024, Ruibin Liang & Lan Guan
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 37
- downloads
- 0
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.