Abstract
The sugar transporter Lactose permease (LacY) of Escherichia coli has become a prototype to understand the underlying molecular details of membrane transport. Crystal structures have trapped the protein in sugar-bound states facing the periplasm, but with narrow openings unable to accommodate sugar. Therefore, the molecular details of sugar uptake remain elusive. In this work, we have used extended simulations and metadynamics sampling to explore a putative sugar-uptake pathway and associated free energy landscape. We found an entrance at helix-pair 2 and 11, which involved lipid head groups and residues Gln 241 and Gln 359. Furthermore, the protein displayed high flexibility on the periplasmic side of Phe 27, which is located at the narrowest section of the pathway. Interactions to Phe 27 enabled passage into the binding site, which was associated with a 24 ± 4 kJ/mol binding free energy in excellent agreement with an independent binding free energy calculation and experimental data. Two free energy minima corresponding to the two possible binding poses of the lactose analog β-D-galactopyranosyl-1-thio-β-D-galactopyranoside (TDG) were aligned with the crystal structure-binding pocket. This work outlines the chemical environment of a putative periplasmic sugar pathway and paves way for understanding substrate affinity and specificity in LacY.