β-Lactoglobulin's Conformational Requirements for Ligand Binding at the Calyx and the Dimer Interphase: a Flexible Docking Study

β-lactoglobulin (BLG) is an abundant milk protein relevant for industry and biotechnology, due significantly to its ability to bind a wide range of polar and apolar ligands. While hydrophobic ligand sites are known, sites for hydrophilic ligands such as the prevalent milk sugar, lactose, remain undetermined. Through the use of molecular docking we first, analyzed the known fatty acid binding sites in order to dissect their atomistic determinants and second, predicted the interaction sites for lactose with monomeric and dimeric BLG. We validated our approach against BLG structures co-crystallized with ligands and report a computational setup with a reduced number of flexible residues that is able to reproduce experimental results with high precision. Blind dockings with and without flexible side chains on BLG showed that: i) 13 experimentally-determined ligands fit the calyx requiring minimal movement of up to 7 residues out of the 23 that constitute this binding site. ii) Lactose does not bind the calyx despite conformational flexibility, but binds the dimer interface and an alternate Site C. iii) Results point to a probable lactolation site in the BLG dimer interface, at K141, consistent with previous biochemical findings. In contrast, no accessible lysines are found near Site C. iv) lactose forms hydrogen bonds with residues from both monomers stabilizing the dimer through a claw-like structure. Overall, these results improve our understanding of BLG's binding sites, importantly narrowing down the calyx residues that control ligand binding. Moreover, our results emphasize th