In this work, we present a study of the interaction between human serum albumin (HSA) and acetylsalicylic acid (ASA, C9H8O4) by molecular dynamics simulations (MD). Starting from an experimentally resolved structure of the complex, we performed the extraction of the ligand by means of the application of an external force. After stabilization of the system, we quantified the force used to remove the ASA from its specific site of binding to HSA and calculated the mechanical nonequilibrium external work done during this process. We obtain a reasonable value for the upper boundary of the Gibbs free energy difference (an equilibrium thermodynamic potential) between the complexed and noncomplexed states. To achieve this goal, we used the finite sampling estimator of the average work, calculated from the Jarzynski Equality. To evaluate the effect of the solvent, we calculated the so-called “viscous work,” that is, the work done to move the aspirin in the same trajectory through the solvent in absence of the protein, so as to assess the relevance of its contribution to the total work. The results are in good agreement with the available experimental data for the albumin affinity constant for aspirin, obtained through quenching fluorescence methods.
Conformational relaxations of human serum albumin studied by molecular dynamics simulations with pressure jumps