Structural and energetic requirements for a second binding site at the dimeric β-lactoglobulin interface

Through experimental and theoretical approaches, it has been shown that bovine β-lactoglobulin (βlg) uses its hydrophobic cavity or calyx as the primary binding site for hydrophobic molecules, whereas the existence of a second ligand binding site at the dimeric interface has only been structurally identified for vitamin D3 (VD3). This binding exists even in the thermally denatured state, suggesting the prevalence of this secondary site. Although crystallographic experiments have suggested that VD3 can bind to both monomeric and dimeric states without significant structural differences, theoretical and experimental reports have proposed some structural requirements. Thus, in this study, based on known experimental data, the dynamic interaction of VD3 with the monomeric or dimeric forms of βlg was investigated through a protocol combining blind docking and 2 microsecond molecular dynamics simulations coupled with binding free energy and per-residue binding free energy decomposition analyses using the Molecular Mechanics Generalized Born Surface Area approach. Binding free energy calculations allowed us to estimate the energetic differences of coupling VD3 at the calyx and the dimeric interface for the monomeric or dimeric state, revealing that the dimeric structure is required to form a stable complex with VD3 at the dimeric interface. This also has an important impact on the dimerization process, whereas although the monomeric state also forms a stable complex with VD3 at the dimeric interface, the incorporation of the entropy component contributed to producing a marginally favorable binding free energy. Finally, the per-residue decomposition analysis provided energetic information about the most relevant residues in stabilizing the different systems.