Molecular simulations of β-lactoglobulin complexed with fatty acids reveal the structural basis of ligand affinity to internal and possible external binding sites.

The interaction of saturated fatty acids of different length (C8:0 to C18:0) with β-lactoglobulin (βLG) was investigated by molecular dynamics simulation and docking approaches. The results show that the presence of such ligands in the hydrophobic central cavity of βLG, known as the protein calyx, determines an enhancement of atomic fluctuations compared with the unliganded form, especially for loops at the entrance of the binding site. Concerted motions are evidenced for protein regions that could favor the binding of ligands. The mechanism of anchoring of fatty acids of different length is similar for the carboxylate head-group, through electrostatic interactions with the side chains of Lys60/Lys69. The key protein residues to secure the hydrocarbon chain are Phe105/Met107, which adapt their conformation upon ligand binding. In particular, Phe105 provides an additional hydrophobic clamp only for the tail of the two fatty acids with the longest chains, palmitic, and stearic acid, which are known to bind βLG with a high affinity. The search of additional external binding sites for fatty acids, distinct from the calyx, was also carried out for palmitic acid. Two external sites with a lower affinity were identified as secondary sites, one consisting in a hydrophobic cavity allowing two distinct binding modes for the fatty acid, and the other corresponding to a surface crevice close to the protein α-helix. The overall results provide a comprehensive picture of the dynamical behavior of βLG in complex with fatty acids, and elucidate the structural basis of the binding of these physiological ligands.