An estimation method of binding free energy in terms of ABEEMσπ/MM and continuum electrostatics fused into LIE method

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. The linear interaction energy method is combined with atom-bond electronegativity equalization method at σπ level Force field (fused into molecular mechanics) and generalized Born continuum model calculation of electrostatic solvation for the estimation of the absolute free energy of binding. The parameters of this method are calibrated by using a training set of 24 HIV-1 protease-inhibitor complexes (PDB entry 1AAQ). A correlation coefficient of 0.93 was obtained with a root mean square deviation of 0.70 kcal mol(-1) . This approach is further tested on seven inhibitor and protease complexes, and it provides small root mean square deviation between the calculated binding free energy and experimental binding free energy without reparametrization. By comparing the radii of gyration and the hydrogen bond distances between ligand and protein of three training model molecules, the consistent comparison result of binding free energy is obtained. It proves that this method of calculating the binding free energy with appropriate structural analysis can be applied to quickly assess new inhibitors of HIV-1 proteases. To test whether the parameters of this method can apply to other drug targets, we have validated this method for the drug target cyclooxygenase-2.