The resulting noncovalent bonding of the salicylic acid to ovine COX-1 after bromoaspirin and aspirin acetylation by Ser530 is investigated within the scope of density functional theory considering a 6.5 Å radius binding pocket. We have not only took full advantage of published X-ray structural data for the ovine COX-1 cocrystallized with bromoaspirin, but we also have improved that data through computation, finding good estimates for the hydrogen atom positions at the residues of the binding pocket, and repositioning the Ser530Ac[Br;H] lateral chain and salicylic acid by total energy minimization procedures employing LDA and GGA+D exchange-correlation functionals. Using bromoaspirin as a template, we have simulated the positioning of aspirin in the binding pocket, estimating its interaction energy with each of its neighbor COX-1 residues. We demonstrate that the binding energies of bromoaspirin and aspirin to COX-1 are very close when second-order quantum refinements of the structural data are performed, which points to an explanation on why the IC(50) values for the 126 μM COX-1 activity of both bromoaspirin and aspirin are practically the same. Attracting and repelling residues were identified, being shown that Arg120 is the most effective residue attracting the salicylic acid, followed by Ala527, Leu531, Leu359, and Ser353. On the other hand, Glu524 was found the most effective repulsive residue (strength interaction comparable to Arg120).
© 2012 American Chemical Society