A theoretical elucidation of bilirubin interaction with HSA's lysines: First electrostatic binding site in IIA subdomain

The electrostatic interaction of amino acid lysines 190, 195 and 199 of human serum albumin (HSA) with bilirubin have been investigated using molecular dynamic simulations, QM and QM/MM minimization methods. In this study two methodological approaches have been employed. In the first approach X-ray structure and the structure obtained from the molecular dynamic simulation of subdomain IIA of HSA in vacuum have been utilized. Interactions have been evaluated with the segment 186-200 of the cited subdomain. Calculations on the X-ray structure of above segment indicate an effective interaction of the lysine 195 with bilirubin, although that of the lysine 190 is also found considerable in this structure. Performing simulation in vacuum, it has been revealed that except for the lysine 195, the other two lysine residues (190 and 199) could not be considered as centers of interaction. Such finding, which is in accord with experimental data, lends support to the procedure employed in this study. NBO analyses suggest that tasks to achieve a structure indicating bilirubin interaction with the lysine 195 from the 186-200 segment extracted from X-ray structure, results in a structure that lacks any electrostatic interaction. In fact, it has been found that the stability of the latter species can be attributed to the H-bonding interaction of the glutamate 188 with both bilirubin and the lysine 195. Further NBO analysis on the structure of the same species, while achieved after molecular dynamic simulation on subdomain IIA in vacuum has revealed that a favorable electrostatic interaction between the lysine 195 and bilirubin has occurred. Besides, H-bonding interaction of the glutamate 188 with bilirubin has been evident in the same species. For the second approach, presence of water molecules and ions has been considered to simulate condensed medium. Applying docking, conformational sampling, and QM/MM minimization steps in sequence, a structure has been achieved which presents a specific interaction between ε-NH3+ group of the lysine 195 residue and the lactam oxygen atom of bilirubin. NBO analyses suggest that above electrostatic interaction is combined with hydrogen bonding interaction between same two groups. Moreover, a hydrogen bond between oxygen atom of bilirubin's acetate group and α-NH group of lysine 195 has been observed. Molecular orbital calculations have been presented which support the NBO analyses.