Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter

Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter View ORCID Profile Ruibin Liang , View ORCID Profile Lan Guan doi: https://doi.org/10.1101/2024.08.21.608993 Ruibin Liang 1 Department of Chemistry and Biochemistry, Texas Tech University , Lubbock, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Ruibin Liang For correspondence: rliang{at}ttu.edu lan.guan{at}ttuhsc.edu Lan Guan 2 Department of Cell Physiology and Molecular Biophysics, Center for Membrane Protein Research, Texas Tech University Health Sciences Center, School of Medicine , Lubbock, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Lan Guan For correspondence: rliang{at}ttu.edu lan.guan{at}ttuhsc.edu Abstract Full Text Info/History Metrics Supplementary material Preview PDF 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 (MelB St ), 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 xenobiotics 1 . Certain human homologs have become promising targets for drug development due to their vital roles in nutrient and drug transport 2 . For example, the MelB homolog MFSD2A is crucial for the uptake of the essential lipid lysophosphatidylcholine in the brain 3 . Melibiose permease of Salmonella enterica serovar Typhimurium (MelB St ) 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 transporters 1 , 4 – 8 . It has been experimentally well-characterized by various biochemical/biophysical techniques in combination with genetic modification 9 – 15 and structural analysis 4 , 16 – 19 . Additionally, extensive experimental data has been made available for MelB of E. coli 20 – 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 mechanism 10 , 28 . More recently, the correlations between substrate-binding affinities and protein conformational states have been identified through structural and binding analyses 19 , 29 . While the cation-binding affinity appears largely independent of protein conformation, it was hypothesized that the sugar-binding affinity significantly decreases when MelB St is in the inward-facing (IF) state 19 , 29 . Download figure Open in new tab Figure 1. The transport cycle of the MelB symporter. The current study describes the melibiose (Mel, black double square) translocation between the periplasm and cytoplasm in the presence or absence of a bound coupling cation (Na + , red dot). The translocation of melibiose into cytoplasm starts from the OF F state ([2]) and proceeds to the IF F state ([6]) through several intermediate states such as the OF B ([3]), OC B ([4]), and IF B ([5]). Reversal of this process results in melibiose translocation into periplasm. These key intermediate states are highlighted in bold text, with the protein depicted in solid colors. See the main text for a detailed definition of the states. To date, the structures for two major conformational states of MelB St have been resolved: the apo- or sugar-bound outward-facing (OF) state 16 , 17 , and the sugar-released, Na + -bound IF state 18 . 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 cooperative 28 , 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 identified 19 , 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 uniporters 30 – 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 MelB St 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 states 17 , 19 , the string method 43 , 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 exchange 4 , 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 MelB St 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 MelB St 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 MelB St in its outward-facing sugar-free state (OF F ) (image ID 0). Binding of melibiose from the periplasmic side leads to the outward-facing sugar-bound state (OF B ) (image ID 13, Fig. S1 ), which exhibited a ΔG of −5.0 kcal/mol ( Table 1 , ΔG OF 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 OF B to the inward-facing sugar-bound state (IF B ) (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 (OC B ) during the conformational transition. The OC B state is a transient state, and the IF B state is higher in free energy than the OF B 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 (IF B , image ID 32), leaving behind a Na + -bound, sugar-free MelB St . 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 ). Download figure Open in new tab Figure 2. The potential of mean force (PMF) for the translocation of the melibiose across the WT MelB St in the Na + -bound state along the minimum free energy pathway (represented as a string of images 0 to 32). Top insets: representative snapshots of the MelB St (blue ribbons) and melibiose (spheres) in varied intermediate states. The bottom and top sides of each protein structure indicate the periplasmic and cytoplasmic sides of the membrane, respectively. The Asp19 and Asp124 residues, which are the major melibiose binding residues, are depicted as orange balls and sticks. The atoms in melibiose molecule are depicted as spheres. View this table: View inline View popup Download powerpoint Table 1. Free energy barriers (ΔG′s) and net free energy changes (ΔG ‡ ′s). Free energy barriers and net free energy changes (in kcal/mol) for melibiose translocation in the WT Na + bound, WT Na + unbound, and the D59C mutant Na + unbound systems were analyzed for the overall translocation process and several key steps in this process. 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 MelB St 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. Download figure Open in new tab Figure 3. Free energy surfaces (FES) for the transport of the melibiose through the WT MelB St in the Na + bound state. (A) FES spanned by the first principle component of the backbone (PC1) vs. melibiose transport coordinate Z. An illustrative structure indicating the scale of the Z coordinate is shown on the right, and key residues in the periplasmic gate, binding site and cytoplasmic gate are highlighted in red, orange and green, respectively (see main text). (B) FES spanned by interhelical distance between helices I and VII on the periplasmic side. An illustrative structure indicating the helices I (red) and VII (green) is shown on the right. (C) FES spanned by interhelical distance between helices IV and X on the cytoplasmic side. An illustrative structure indicating the helices IV (red) and X (green) is shown on the right. The approximate minimum free energy pathways (MFEP) tracking the major basins on the FES are indicated as black lines, and the regions corresponding to the intermediate states are labeled by arrows. All FES plots share the same color bar as (A). 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 OF F →OF B process. For example, from OF F (Z = −30 Å, PC1 = −3.2 Å) to OF B (Z = −3 Å, PC1 = −2.5 Å), the PC1 values fluctuate back and forth and has a net increase of ∼0.7 Å after reaching the OF B 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 OF F to ∼19 Å at OF B 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 OF B →IF B 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 IF B state, the melibiose release to the cytoplasmic bulk begins (IF B ⟶IF F ). The process goes through an intermediate state IF C , 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 IF B to IF C , 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 IF B , IF C , and IF F ( 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 (OF F →OF B , red vs. orange curves in Fig. 4A ). It then largely decreases by more than 3 Å during conformational transition (OF B →IF B , orange vs. green curves in Fig. 4A ). Following this, it remains almost the same during cytoplasmic release (IF B →IF F , 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 OF B →IF B transition. Thus, in general, the OF F →IF F process is accompanied by the widening and shrinking of the periplasmic and cytoplasmic terminals of the path in a reciprocal manner. Download figure Open in new tab Figure 4. Change of pore radius profile coupled to melibiose transport through WT MelB St in the Na + bound state. (A): Pore radii profile across the MelB as a function of the Z coordinate (relative to the COM of Asp124 and Asp19) along the entire protein (x-axis) for each different state during the melibiose translocation. As the melibiose is translocated from the periplasm to cytoplasm, the system transitions from the OF F (red) to IF F (blue) states through the OF B (orange), OC B (black), IF B (green), and IF C (brown) intermediate states, each of which features different pore radius profiles. The IF C state corresponds to the melibiose passing through the cytoplasmic constricted region near Z=11 Å. An illustrative structure indicating the scale of the Z coordinate is shown on the right. (B) and (C): Translocation of melibiose from the periplasmic to cytoplasmic sides of MelB St induces the change of pore radii measured at Z=-2.5 Å (periplasmic sugar-binding site) and Z=11 Å (cytoplasmic constricted region). The melibiose translocation process is represented as the transitioning of the system from the OF F to IF F states through multiple intermediate states along the x-axes. The protein structures below the plots serve as visual guides for the location of the melibiose in different states. 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 MelB St 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 MelB St , highlighting the importance of this residue in the functional cycle 13 . 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 previously 13 , 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 (OF B →IF B ) 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 result 19 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 MelB St 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 ΔG OF binding ∼ 1.4 kcal/mol, in contrast to the exergonic process in the Na + bound state (ΔG OF binding ∼ −5.0 kcal/mol). Second, the endergonic protein conformational transition (image IDs from 15 to 21) exhibits ΔG OF 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 ΔG IF 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). Download figure Open in new tab Figure 5. PMFs for the transport of the melibiose across MelB St along the MFEPs represented as strings of images (0 to 33). The PMFs are calculated in the WT, Na + bound (red) and unbound (blue) states, and the D59C mutant in the Na + unbound state (green). 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 gradient 17 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 ΔG OF 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 measurements 17 . After the OF B state is reached, the OF B ⟶IF B transition together with the cytoplasmic sugar release (IF B ⟶IF F ) (image IDs from 13 to 32) results in a combined high barrier of ∼10.6 kcal/mol (ΔG OF 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 methods 9 , 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 OF B state for each system. The interaction energies between the melibiose and MelB St 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 MelB St , 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 MelB St 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 OF B 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 studies 13 , 16 , 24 , 27 , 46 : the D124C or Y120C mutants can bind melibiose but lose Na + -coupled melibiose active transport 16 . 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 states 17 , 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 simulations 31 , 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 analyses 9 , 10 , 12 , 13 , 18 , 19 , 28 . First, all simulation data consistently showed that Na + binding increases the melibiose-binding affinity of WT MelB St , in agreement with previous experimental measurements 9 , 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 analysis 19 . Third, the overall free energy landscape reveals that the outward-facing conformation of the WT has the highest sugar-binding affinity, and the OF B 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 state 19 . During the melibiose translocation, only a modest rate-limiting barrier is encountered in the OF B →OF F 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 (OC B ) facilitates the OF B →IF B transition. All these features make MelB St 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 MelB St and its mutants 10 , 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 data 16 , 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 MelB St 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 activity 16 , 17 , similar to the E325A mutant of LacY 7 , 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 MelB St 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 MelB St 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 MelB St with Na + bound, (2) WT MelB St without Na + , and (3) D59C mutant without Na + . The OF and IF conformational states were derived from crystal (PDB code 7L16 17 ) and CryoEM (PDB code 8T60 19 ) 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 interface 48 . 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 OF F →OF B and IF B →IF F transitions. Next, the geodesic interpolation algorithm 49 generated 10 intermediate images of the transmembrane backbone for the OF B →IF B conformational transition. These backbone structures guided restrained simulations that dynamically relaxed all atoms in the images. All simulations used the CHARMM36m 50 – 56 and TIP3P force fields 57 , 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 transporters 31 , 33 , 35 and protein complexes 59 . The initial string consisted of 33 images, representing the entire OF F ↔ IF F 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 simulations 60 , 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 basis 59 . 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 version 33 of the weighted histogram analysis method 62 , 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 OF F and IF F states on the 2D free energy surfaces (FES) was traced using the zero-temperature string method. The HOLE program 64 was used to analyze the pore radius profiles for all snapshots sampled by the REUS simulations. 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. Data Availability Statement The data supporting this article have been included in the Supplementary Information. Supplemental information The supplemental information contains Figure S1-S12 . It includes a snapshot of the melibiose binding site in the WT, Na + bound MelB St , the evolution of RMSD of the string in the 13-dimensional CV space as a function of iteration number during the string method simulation for all three systems, the backbone RMSD of the string images with respect to the experimental structures, the free energy surfaces and pore radius profiles during the entire melibiose translocation across the WT and D59C mutant of MelB St in the Na + unbound state, the energy decomposition analysis of Na + -protein interaction, and the distances between the coupling residues (D124 and Y120) to the cation-binding and sugar-binding sites. The supplemental information also contains the MFEPs of all three systems as separate text files. This information is available online, free of charge. Download figure Open in new tab All-atom free energy landscapes of the entire melibiose translocation process in wild-type MelB St 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. Acknowledgments 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. Footnotes ↵ 3 Leading Contact Minor modifications to the abstract and discussions to highlight the major findings and significance of the results. References 1. ↵ M. H. Saier , Jr. , J. T. Beatty , A. Goffeau , K. T. Harley , W. H. Heijne , S. C. Huang , D. L. Jack , P. S. Jahn , K. Lew , J. Liu , S. S. Pao , I. T. Paulsen , T. T. Tseng and P. S. Virk , The major facilitator superfamily , J Mol Microbiol Biotechnol , 1999 , 1 , 257 – 279 . OpenUrl PubMed 2. ↵ L. Lin , S. W. Yee , R. B. Kim and K. M. Giacomini , SLC transporters as therapeutic targets: emerging opportunities , Nat Rev Drug Discov , 2015 , 14 , 543 – 560 . OpenUrl 3. ↵ L. N. Nguyen , D. Ma , G. Shui , P. Wong , A. Cazenave-Gassiot , X. Zhang , M. R. Wenk , E. L. Goh and D. L. Silver , Mfsd2a is a transporter for the essential omega-3 fatty acid docosahexaenoic acid , Nature , 2014 , 509 , 503 – 506 . OpenUrl CrossRef PubMed Web of Science 4. ↵ L. Guan , Na(+)/Melibiose Membrane Transport Protein , MelB.Journal , 2018 , DOI: 10.1007/978-3-642-35943-9_10082-1 . OpenUrl CrossRef 5. S. Niiya , Y. Moriyama , M. Futai and T. Tsuchiya , Cation coupling to melibiose transport in Salmonella typhimurium , J. Bacteriol ., 1980 , 144 , 192 – 199 . OpenUrl Abstract / FREE Full Text 6. B. Poolman , J. Knol , C. van der Does , P. J. Henderson , W. J. Liang , G. Leblanc , T. Pourcher and I. Mus-Veteau , Cation and sugar selectivity determinants in a novel family of transport proteins , Mol. Microbiol ., 1996 , 19 , 911 – 922 . OpenUrl CrossRef PubMed Web of Science 7. ↵ L. Guan and H. R. Kaback , Lessons from lactose permease , Annu. Rev. Biophys. Biomol. Struct ., 2006 , 35 , 67 – 91 . OpenUrl CrossRef PubMed Web of Science 8. ↵ N. Yan , Structural Biology of the Major Facilitator Superfamily Transporters , Annu Rev Biophys , 2015 , 44 , 257 – 283 . OpenUrl CrossRef PubMed 9. ↵ L. Guan , S. Nurva and S. P. Ankeshwarapu , Mechanism of melibiose/cation symport of the melibiose permease of Salmonella typhimurium , J. Biol. Chem ., 2011 , 286 , 6367 – 6374 . OpenUrl Abstract / FREE Full Text 10. ↵ P. Hariharan and L. Guan , Cooperative binding ensures the obligatory melibiose/Na+ cotransport in MelB , J. Gen. Physiol ., 2021 , 153 . 11. S. V. Jakkula and L. Guan , Reduced Na(+) affinity increases turnover of Salmonella enterica serovar Typhimurium MelB , J. Bacteriol ., 2012 , 194 , 5538 – 5544 . OpenUrl Abstract / FREE Full Text 12. ↵ S. Katsube , R. Liang , A. Amin , P. Hariharan and L. Guan , Molecular Basis for the Cation Selectivity of Salmonella typhimurium Melibiose Permease , J. Mol. Biol ., 2022 , 434 , 167598 . 13. ↵ K. J. Markham , E. B. Tikhonova , A. C. Scarpa , P. Hariharan , S. Katsube and L. Guan , Complete cysteine-scanning mutagenesis of the Salmonella typhimurium melibiose permease , J. Biol. Chem ., 2021 , 297 , 101090 . 14. P. Hariharan , E. Tikhonova , J. Medeiros-Silva , A. Jeucken , M. V. Bogdanov , W. Dowhan , J. F. Brouwers , M. Weingarth and L. Guan , Structural and functional characterization of protein-lipid interactions of the Salmonella typhimurium melibiose transporter MelB , BMC Biol ., 2018 , 16 , 85 . 15. ↵ N. Blaimschein , P. Hariharan , S. Manioglu , L. Guan and D. J. Muller , Substrate-binding guides individual melibiose permeases MelB to structurally soften and to destabilize cytoplasmic middle-loop C3 , Structure , 2023 , 31 , 58 – 67 e54 . OpenUrl CrossRef 16. ↵ A. S. Ethayathulla , M. S. Yousef , A. Amin , G. Leblanc , H. R. Kaback and L. Guan , Structure-based mechanism for Na(+)/melibiose symport by MelB , Nat. Commun ., 2014 , 5 , 3009 . OpenUrl CrossRef PubMed 17. ↵ L. Guan and P. Hariharan , X-ray crystallography reveals molecular recognition mechanism for sugar binding in a melibiose transporter MelB , Commun Biol , 2021 , 4 , 931 . OpenUrl 18. ↵ P. Hariharan , A. Bakhtiiari , R. Liang and L. Guan , Distinct roles of the major binding residues in the cation-binding pocket of the melibiose transporter MelB , J. Biol. Chem ., 2024 , DOI: 10.1016/j.jbc.2024.107427 . , 107427 . OpenUrl CrossRef 19. ↵ P. Hariharan , Y. Shi , S. Katsube , K. Willibal , N. D. Burrows , P. Mitchell , A. Bakhtiiari , S. Stanfield , E. Pardon , H. R. Kaback , R. Liang , J. Steyaert , R. Viner and L. Guan , Mobile barrier mechanisms for Na(+)-coupled symport in an MFS sugar transporter , Elife , 2024 , 12 . 20. ↵ D. M. Wilson and T. H. Wilson , Cation specificity for sugar substrates of the melibiose carrier in Escherichia coli , Biochim Biophys Acta , 1987 , 904 , 191 – 200 . OpenUrl CrossRef PubMed Web of Science 21. T. H. Wilson and P. Z. Ding , Sodium-substrate cotransport in bacteria , Biochim Biophys Acta , 2001 , 1505 , 121 – 130 . OpenUrl CrossRef PubMed Web of Science 22. E. Damiano-Forano , M. Bassilana and G. Leblanc , Sugar binding properties of the melibiose permease in Escherichia coli membrane vesicles. Effects of Na(+) and H(+) concentrations , J. Biol. Chem ., 1986 , 261 , 6893 – 6899 . OpenUrl Abstract / FREE Full Text 23. T. Pourcher , M. L. Zani and G. Leblanc , Mutagenesis of acidic residues in putative membrane-spanning segments of the melibiose permease of Escherichia coli . I. Effect on Na(+)-dependent transport and binding properties , J. Biol. Chem ., 1993 , 268 , 3209 – 3215 . OpenUrl Abstract / FREE Full Text 24. ↵ M. L. Zani , T. Pourcher and G. Leblanc , Mutagenesis of acidic residues in putative membrane-spanning segments of the melibiose permease of Escherichia coli . II. Effect on cationic selectivity and coupling properties , J. Biol. Chem ., 1993 , 268 , 3216 – 3221 . OpenUrl Abstract / FREE Full Text 25. G. Leblanc , T. Pourcher and M. L. Zani , The melibiose permease of Escherichia coli: importance of the NH2-terminal domains for cation recognition by the Na+/sugar cotransporter , Soc. Gen. Physiol. Ser ., 1993 , 48 , 213 – 227 . OpenUrl PubMed 26. C. Ganea , T. Pourcher , G. Leblanc and K. Fendler , Evidence for intraprotein charge transfer during the transport activity of the melibiose permease from Escherichia coli , Biochemistry , 2001 , 40 , 13744 – 13752 . OpenUrl CrossRef PubMed Web of Science 27. ↵ M. Granell , X. Leon , G. Leblanc , E. Padros and V. A. Lorenz-Fonfria , Structural insights into the activation mechanism of melibiose permease by sodium binding , Proc Natl Acad Sci U S A , 2010 , 107 , 22078 – 22083 . OpenUrl Abstract / FREE Full Text 28. ↵ P. Hariharan and L. Guan , Thermodynamic cooperativity of cosubstrate binding and cation selectivity of Salmonella typhimurium MelB , J. Gen. Physiol ., 2017 , 149 , 1029 – 1039 . OpenUrl Abstract / FREE Full Text 29. ↵ S. Katsube , K. Willibal , S. Vemulapally , P. Hariharan , E. Tikhonova , E. Pardon , H. R. Kaback , J. Steyaert and L. Guan , In vivo and in vitro characterizations of melibiose permease (MelB) conformation-dependent nanobodies reveal sugar-binding mechanisms , J. Biol. Chem ., 2023 , 299 , 104967 . 30. ↵ J. Feng , B. Selvam and D. Shukla , How do antiporters exchange substrates across the cell membrane? An atomic-level description of the complete exchange cycle in NarK , Structure , 2021 , 29 , 922 – 933.e923 . OpenUrl 31. ↵ M. Takemoto , Y. Lee , R. Ishitani and O. Nureki , Free Energy Landscape for the Entire Transport Cycle of Triose-Phosphate/Phosphate Translocator , Structure , 2018 , 26 , 1284 – 1296.e1284 . OpenUrl CrossRef PubMed 32. R. Alhadeff and A. Warshel , Simulating the function of sodium/proton antiporters , Proc. Natl. Acad. Sci. U. S. A ., 2015 , 112 , 12378 – 12383 . OpenUrl Abstract / FREE Full Text 33. ↵ M. Moradi , G. Enkavi and E. Tajkhorshid , Atomic-level characterization of transport cycle thermodynamics in the glycerol-3-phosphate:phosphate antiporter , Nat. Commun ., 2015 , 6 , 8393 . OpenUrl CrossRef PubMed 34. B. Selvam , Y.-C. Yu , L.-Q. Chen and D. Shukla , Molecular Basis of the Glucose Transport Mechanism in Plants , ACS Central Science , 2019 , 5 , 1085 – 1096 . OpenUrl CrossRef 35. ↵ M. Ke , Y. Yuan , X. Jiang , N. Yan and H. Gong , Molecular determinants for the thermodynamic and functional divergence of uniporter GLUT1 and proton symporter XylE , PLoS Comput. Biol ., 2017 , 13 , e1005603 . OpenUrl CrossRef 36. K. J. Cheng , B. Selvam , L.-Q. Chen and D. Shukla , Distinct Substrate Transport Mechanism Identified in Homologous Sugar Transporters , J. Phys. Chem. B , 2019 , 123 , 8411 – 8418 . OpenUrl CrossRef 37. B. Selvam , S. Mittal and D. Shukla , Free Energy Landscape of the Complete Transport Cycle in a Key Bacterial Transporter , ACS Central Science , 2018 , 4 , 1146 – 1154 . OpenUrl 38. D. Kimanius , E. Lindahl and M. Andersson , Uptake dynamics in the Lactose permease (LacY) membrane protein transporter , Sci. Rep ., 2018 , 8 , 14324 . 39. M. C. Chan , B. Selvam , H. J. Young , E. Procko and D. Shukla , The substrate import mechanism of the human serotonin transporter , Biophys. J ., 2022 , 121 , 715 – 730 . OpenUrl CrossRef 40. M. Ø. Jensen , Y. Yin , E. Tajkhorshid and K. Schulten , Sugar Transport across Lactose Permease Probed by Steered Molecular Dynamics , Biophys. J ., 2007 , 93 , 92 – 102 . OpenUrl CrossRef PubMed 41. J. Li , Z. Zhao and E. Tajkhorshid , Locking Two Rigid-body Bundles in an Outward-Facing Conformation: The Ion-coupling Mechanism in a LeuT-fold Transporter , Sci. Rep ., 2019 , 9 , 19479 . 42. ↵ R. Gradisch , K. Schlögl , E. Lazzarin , M. Niello , J. Maier , F. P. Mayer , L. Alves da Silva , S. M. C. Skopec , R. D. Blakely , H. H. Sitte , M. D. Mihovilovic and T. Stockner , Ligand coupling mechanism of the human serotonin transporter differentiates substrates from inhibitors , Nat. Commun ., 2024 , 15 , 417 . 43. ↵ H. Chen , D. Ogden , S. Pant , W. Cai , E. Tajkhorshid , M. Moradi , B. Roux and C. Chipot , A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls , J. Chem. Theory Comput ., 2022 , 18 , 1406 – 1422 . OpenUrl CrossRef 44. ↵ A. C. Pan , D. Sezer and B. Roux , Finding Transition Pathways Using the String Method with Swarms of Trajectories , J. Phys. Chem. B , 2008 , 112 , 3432 – 3440 . OpenUrl CrossRef PubMed 45. ↵ L. Guan , S. V. Jakkula , A. A. Hodkoff and Y. Su , Role of Gly117 in the cation/melibiose symport of MelB of Salmonella typhimurium , Biochemistry , 2012 , 51 , 2950 – 2957 . OpenUrl CrossRef PubMed Web of Science 46. ↵ M. L. Zani , T. Pourcher and G. Leblanc , Mutation of polar and charged residues in the hydrophobic NH2-terminal domains of the melibiose permease of Escherichia coli , J. Biol. Chem ., 1994 , 269 , 24883 – 24889 . OpenUrl Abstract / FREE Full Text 47. ↵ N. Carrasco , L. M. Antes , M. S. Poonian and H. R. Kaback , lac permease of Escherichia coli : histidine-322 and glutamic acid-325 may be components of a charge-relay system , Biochemistry , 1986 , 25 , 4486 – 4488 . OpenUrl CrossRef PubMed Web of Science 48. ↵ S. Jo , T. Kim , V. G. Iyer and W. Im , CHARMM-GUI: A web-based graphical user interface for CHARMM , J. Comput. Chem ., 2008 , 29 , 1859 – 1865 . OpenUrl CrossRef PubMed 49. ↵ X. Zhu , K. C. Thompson and T. J. Martínez , Geodesic interpolation for reaction pathways , J. Chem. Phys ., 2019 , 150 , 164103 . 50. ↵ J. Huang , S. Rauscher , G. Nawrocki , T. Ran , M. Feig , B. L. de Groot , H. Grubmüller and A. D. MacKerell , CHARMM36m: an improved force field for folded and intrinsically disordered proteins , Nat. Methods , 2017 , 14 , 71 – 73 . OpenUrl CrossRef PubMed 51. R. B. Best , X. Zhu , J. Shim , P. E. M. Lopes , J. Mittal , M. Feig and A. D. MacKerell , 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 Angles , J. Chem. Theory Comput ., 2012 , 8 , 3257 – 3273 . OpenUrl CrossRef PubMed 52. A. D. MacKerell , Jr. , M. Feig and C. L. Brooks , Improved Treatment of the Protein Backbone in Empirical Force Fields , J. Am. Chem. Soc ., 2004 , 126 , 698 – 699 . OpenUrl CrossRef PubMed Web of Science 53. A. D. MacKerell , D. Bashford , M. Bellott , R. L. Dunbrack , J. D. Evanseck , M. J. Field , S. Fischer , J. Gao , H. Guo , S. Ha , D. Joseph-McCarthy , L. Kuchnir , K. Kuczera , F. T. K. Lau , C. Mattos , S. Michnick , T. Ngo , D. T. Nguyen , B. Prodhom , W. E. Reiher , B. Roux , M. Schlenkrich , J. C. Smith , R. Stote , J. Straub , M. Watanabe , J. Wiorkiewicz-Kuczera , D. Yin and M. Karplus , All-atom empirical potential for molecular modeling and dynamics studies of proteins , J. Phys. Chem. B , 1998 , 102 , 3586 – 3616 . OpenUrl CrossRef PubMed 54. J. B. Klauda , R. M. Venable , J. A. Freites , J. W. O’Connor , D. J. Tobias , C. Mondragon-Ramirez , I. Vorobyov , A. D. MacKerell and R. W. Pastor , Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types , J. Phys. Chem. B , 2010 , 114 , 7830 – 7843 . OpenUrl CrossRef PubMed 55. O. Guvench , S. N. Greene , G. Kamath , J. W. Brady , R. M. Venable , R. W. Pastor and A. D. Mackerell Jr , Additive empirical force field for hexopyranose monosaccharides , J. Comput. Chem ., 2008 , 29 , 2543 – 2564 . OpenUrl CrossRef PubMed 56. ↵ O. Guvench , E. Hatcher , R. M. Venable , R. W. Pastor and A. D. MacKerell , Jr . ., CHARMM Additive All-Atom Force Field for Glycosidic Linkages between Hexopyranoses , J. Chem. Theory Comput ., 2009 , 5 , 2353 – 2370 . OpenUrl CrossRef PubMed 57. ↵ W. Jorgensen , J. Chandrasekhar , J. Madura , R. Impey and M. Klein , Comparison of simple potential functions for simulating liquid water , J. Chem. Phys ., 1983 , 79 , 926 – 935 . OpenUrl CrossRef PubMed Web of Science 58. ↵ J. C. Phillips , D. J. Hardy , J. D. C. Maia , J. E. Stone , J. V. Ribeiro , R. C. Bernardi , R. Buch , G. Fiorin , J. Hénin , W. Jiang , R. McGreevy , M. C. R. Melo , B. K. Radak , R. D. Skeel , A. Singharoy , Y. Wang , B. Roux , A. Aksimentiev , Z. Luthey-Schulten , L. V. Kalé , K. Schulten , C. Chipot and E. Tajkhorshid , Scalable molecular dynamics on CPU and GPU architectures with NAMD , J. Chem. Phys ., 2020 , 153 , 044130 . 59. ↵ A. Singharoy , C. Chipot , M. Moradi and K. Schulten , Chemomechanical Coupling in Hexameric Protein–Protein Interfaces Harnesses Energy within V-Type ATPases , J. Am. Chem. Soc ., 2017 , 139 , 293 – 310 . OpenUrl CrossRef 60. ↵ Y. Sugita and Y. Okamoto , Replica-exchange molecular dynamics method for protein folding , Chem. Phys. Lett ., 1999 , 314 , 141 – 151 . OpenUrl CrossRef PubMed Web of Science 61. ↵ Y. Sugita , A. Kitao and Y. Okamoto , Multidimensional replica-exchange method for free-energy calculations , J. Chem. Phys ., 2000 , 113 , 6042 – 6051 . OpenUrl CrossRef Web of Science 62. ↵ S. Kumar , D. Bouzida , R. H. Swendsen , P. A. Kollman and J. M. Rosenberg , The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules. 1. The Method , J. Comput. Chem. , 1992 , 13 , 1011 - 1021 . OpenUrl CrossRef PubMed Web of Science 63. ↵ A. Grossfield , WHAM: the weighted histogram analysis method . Journal . 64. ↵ O. S. Smart , J. M. Goodfellow and B. A. Wallace , The pore dimensions of gramicidin A , Biophys. J ., 1993 , 65 , 2455 – 2460 . OpenUrl CrossRef PubMed Web of Science View the discussion thread. Back to top Previous Next Posted September 14, 2024. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter Ruibin Liang , Lan Guan bioRxiv 2024.08.21.608993; doi: https://doi.org/10.1101/2024.08.21.608993 Share This Article: Copy Citation Tools Atomic-Level Free Energy Landscape Reveals Cooperative Symport Mechanism of Melibiose Transporter Ruibin Liang , Lan Guan bioRxiv 2024.08.21.608993; doi: https://doi.org/10.1101/2024.08.21.608993 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Biophysics Subject Areas All Articles Animal Behavior and Cognition (7644) Biochemistry (17728) Bioengineering (13916) Bioinformatics (42037) Biophysics (21489) Cancer Biology (18637) Cell Biology (25553) Clinical Trials (138) Developmental Biology (13401) Ecology (19941) Epidemiology (2067) Evolutionary Biology (24367) Genetics (15622) Genomics (22547) Immunology (17764) Microbiology (40475) Molecular Biology (17208) Neuroscience (88747) Paleontology (667) Pathology (2842) Pharmacology and Toxicology (4834) Physiology (7659) Plant Biology (15175) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9835) Zoology (2272)