Investigation of Paclitaxel Resistant R306C Mutation in β-Tubulin—A Computational Approach

Paclitaxel is the most effective chemotherapeutic agent used for the treatment of a broad spectrum of solid tumors. However, observed paclitaxel resistance in clinical trials presents one of the major obstacles for cancer chemotherapy. Most importantly, resistance due to β-tubulin mutations (R306C) has been intensely debated in recent years. Despite all efforts, mechanism of resistance is still not well understood. In this study, computational techniques were employed to uncover the effect of R306C mutation in the β-tubulin structure and its function. The tools such as I-Mutant, CUPSAT and Fold-X were employed to address the consequence of R306C mutation in the structural stability of β-tubulin. Further, molecular docking and molecular dynamics study was employed to understand the functional impact of β-tubulin mutation. Our results suggest that the R306C mutation causes a significant reduction in the binding affinity between β-tubulin and paclitaxel. Further, docked complex analysis indicates that destruction of conservative hydrogen bond maintained by the residues Arg282 and Gly360 should be responsible for the large conformation changes of the binding pocket in R306C mutant. Finally, molecular dynamics simulations study confirms the stable binding of paclitaxel with native type β-tubulin structure rather than mutant (R306C) type. We certainly believe that this study will provide useful guidance for the development of novel inhibitors that are less susceptible to drug resistance.