Analyses of activity for factor Xa inhibitors based on Monte Carlo simulations

Monte Carlo/Extended Linear Response (MC/ELR) simulations have been conducted on 60 inhibitors of human factor Xa to determine the important interactions associated with their activity. A variety of physicochemical descriptors were configurationally averaged during the course of the simulations of each inhibitor bound to factor Xa and free in water. A regression equation was then derived; it reproduces the experimental inhibition data with a correlation coefficient, r(2), of 0.74, an rms error of 0.67 kcal/mol, and an average unsigned error of 0.60 kcal/mol using only two physically reasonable descriptors. The two factors that emerged as important in determining inhibitory potential are (1) favorable van der Waals interactions between protein and ligand and (2) direct hydrogen bonding between the inhibitor and protein. The conclusions were supported with structural analyses and results of MC/free energy perturbation (FEP) calculations.