A molecular dynamics study of the three-dimensional model of human synovial fluid phospholipase A2--transition state mimic complexes.

Abstract

Different modes of binding of transition state mimics: amide, phosphonate and difluoro ketone, to human synovial fluid phospholipase A2 (HSF PLA2) are studies by molecular dynamics simulations computed in solvent. The results are analysed in the light of primary binding sites. Hydrogen bonding interaction plays an important role for amino acids such as Gly32, Val30, and Glu55, apart from the well known active site residues viz Asp48, Gly25, Gly29, Gly31, His27, His47, Lys62, Phe23, Asn114 and Tyr112. In addition, the hydrogen bonding interaction between Sn-1 tetrahedral phosphonate group of amide and difluoro ketone inhibitors and crystallographic water molecules (H2O 523, H2O 524 and H2O 401) seems to have a significant role. Many of the active site charged residues display considerable movement upon ligand binding. The structural effects of ligand binding were analyzed from RMS deviations of C alpha in the resulting energy-minimized average structures of the receptor-ligand complexes. The values of the RMS deviations differ among the HSF PLA2s, in a pattern that is not the same for the three complexes. This suggests that ligands with different pharmacological efficacies induce different types of conformational changes of the receptor. Our active-orientation model is, at least qualitatively, consistent with experimental data and should be useful for the rational design of more potent inhibitors.