Applying linear interaction energy method for binding affinity calculations of podophyllotoxin analogues with tubulin using continuum solvent model and prediction of cytotoxic activity

Podophyllotoxin and its analogues have important therapeutic value in the treatment of cancer, due to their ability to induce apoptosis in cancer cells in a proliferation-independent manner. These ligands bind to colchicine binding site of tubulin near the α- and β-tubulin interface and interfere with tubulin polymerization. The binding free energies of podophyllotoxin-based inhibitors of tubulin were computed using a linear interaction energy (LIE) method with a surface generalized Born (SGB) continuum solvation model. A training set of 76 podophyllotoxin analogues was used to build a binding affinity model for estimating the free energy of binding for 36 inhibitors (test set) with diverse structural modifications. The average root mean square error (RMSE) between the experimental and predicted binding free energy values was 0.56 kcal/mol which is comparable to the level of accuracy achieved by the most accurate methods, such as free energy perturbation (FEP) or thermodynamic integration (TI). The squared correlation coefficient between experimental and SGB–LIE estimates for the free energy for the test set compounds is also significant (R2 = 0.733). On the basis of the analysis of the binding energy, we propose that the three-dimensional conformation of the A, B, C and D rings is important for interaction with tubulin. On the basis of this insight, 12 analogues of varying ring modification were taken, tested with LIE methodology and then validated with their experimental potencies of tubulin polymerization inhibition. Low levels of RMSE for the majority of inhibitors establish the structure-based LIE method as an efficient tool for generating more potent and specific inhibitors of tubulin by testing rationally designed lead compounds based on podophyllotoxin derivatization.