RESUMEN
Development of a highly accurate prediction model for protein-ligand inhibition has been a major challenge in drug discovery. Herein, we describe a novel predictive model for the inhibition of HIV-1 integrase (IN)-LEDGF/p75 protein-protein interaction. The model was constructed using energy parameters approximated from molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. Chemometric analysis using partial least squares (PLS) regression revealed that solvent accessible surface area energy (ΔGSASA) is the major determinant parameter contributing greatly to the prediction accuracy. PLS prediction model on the ΔGSASA values collected from 41 complexes yielded a strong correlation between the predicted and the actual inhibitory activities (R2 = 0.9666, RMSEC of pIC50 values = 0.0890). Additionally, for the test set of 14 complexes, the model performed satisfactorily with very low pIC50 errors (Q2 = 0.5168, RMSEP = 0.3325). A strong correlation between the buried surface areas on the IN protein, when bound with IN-LEDGF/p75 inhibitors, and the respective ΔGSASA values was also obtained. Furthermore, the current method could identify 'hot spots'of amino acid residues highly influential to the inhibitory activity prediction. This could present fruitful implications in binding site determination and future inhibitor developments targeting protein-protein interactions.Communicated by Ramaswamy H. Sarma.