Steered molecular dynamics simulations of ligand–receptor interaction in lipocalins

Retinol binding protein (RBP) and an engineered lipocalin, DigA16, have been studied using molecular dynamics simulations. Special emphasis has been placed on explaining the ligand–receptor interaction in RBP–retinol and DigA16–digoxigenin complexes, and steered molecular dynamics simulations of 10–20 ns have been carried out for the ligand expulsion process. Digoxigenin is bound deep inside the cavity of DigA16 and forms several stable hydrogen bonds in addition to the hydrophobic van der Waals interaction with the aromatic side-chains. Four crystalline water molecules inside the ligand-binding cavity remain trapped during the simulations. The strongly hydrophobic receptor site of RBP differs considerably from DigA16, and the main source of ligand attraction comes from the phenyl side-chains. The hydrogen bonds between digoxigenin and DigA16 cause the rupture forces on ligand removal in DigA16 and RBP to differ. The mutated DigA16 residues contribute approximately one-half of the digoxigenin interaction energy with DigA16 and, of these, the energetically most important are residues His35, Arg58, Ser87, Tyr88, and Phe114. Potential “sensor loops” were found for both receptors. These are the outlier loops between residues 114–121 and 63–67 for DigA16 and RBP, respectively, and they are located near the entrance of the ligand-binding cavity. Especially, the residues Glu119 (DigA16) and Leu64 (RBP) are critical for sensing. The ligand binding energies have been estimated based on the linear response approximation of binding affinity by using a previous parametrization for retinoids and RBP.