Human serum albumin (HSA) has seven common fatty acid (FA) binding sites. In this study, we used the molecular mechanics Poisson-Boltzmann surface area method to identify high affinity FA binding sites on HSA in terms of binding free energy. Using multiple HSA-FA (myristate, palmitate) complex models constructed by molecular dynamics simulations, two methods were performed in molecular mechanics Poisson-Boltzmann surface area, the "three-trajectory method" and the "single-trajectory method". The former, which is less precise than the latter but may be more accurate as it includes the effects of conformational change upon binding, was used to classify high and low affinity sites. As a result, Sites 2, 4, and 5 were identified as high affinity sites for both FAs. The latter method, which is precise because energies are calculated from snapshots of the same trajectory for HSA-FA complex, was performed to compare the magnitude of binding free energy for these sites. The order of magnitude was 5 > 4 > 2, identical to that of a previous publication by others. In this way, a combination of the two methods was effectively used to identify high affinity sites. This study therefore provides an insight into the quantitative identification of high affinity FA binding sites on HSA.
Molecular simulation approaches have great potentials to provide detailed biophysical insights into HSA as well as the effects of the binding of FAs or other ligands to HSA. This article is part of a Special Issue entitled Serum Albumin.
Human serum albumin (HSA) binds with fatty acids under normal physiologic conditions. To date, there is little published information on the tertiary structure of HSA-fatty acid complex in aqueous solution. In the present study, we used molecular dynamics (MD) simulations to elucidate possible structural changes of HSA brought about by the binding of fatty acids. Both unliganded HSA and HSA-fatty acid complex models for MD calculations were constructed based on the X-ray crystal structures. Five myristates (MYRs) were bound in the HSA-fatty acid complex model. In the present MD study, the motion of domains I and III caused by the binding of MYR molecules increased the radius of gyration of HSA. Root-mean-square fluctuations from the MD simulations revealed that the atomic fluctuations of the specific amino acids at drug-binding site I that can regulate the drug-binding affinity were increased by the binding of MYR molecules. Primary internal motions, characterized by the first three principal components, were observed mainly at domains I and III in the principal component analysis for trajectory data. The directional motion projected on the first principal component of unliganded HSA was conserved in HSA-MYR complex as the third principal directional motion with higher frequency. However, the third principal directional motion in unliganded HSA turned into the first principal directional motion with lower frequency in the HSA-MYR complex. Thus, the present MD study provides insights into the possible conformational changes of HSA caused by the binding of fatty acids.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.