Potential colchicine binding site inhibitors unraveled by virtual screening, molecular dynamics and MM/PBSA

Microtubules have been widely studied in recent decades as an important pharmacological target for the treatment of cancer especially due to its key role in the mitosis process. Among the constituents of the microtubules, αβ-tubulin dimers stand out in view of their four distinct interaction sites, including the so-called colchicine binding site (CBS) - a promising target for the development of new tubulin modulators. When compared to other tubulin sites, targeting the CBS is advantageous because this site is able to host ligands with lower molecular volume and lipophilicity, thus reducing the chances of entailing the phenomenon of multiple drug resistance (MDR) - one of the main reasons of failure in chemotherapy. However, colchicine, the first ligand ever discovered with affinity towards the CBS, despite modulating the action of microtubules, has shown toxicity in clinical studies. Therefore, in order to expand the known chemical space of scaffolds capable of interacting with CBS and to design non-toxic colchicine binding site inhibitors, we conducted a robust virtual screening pipeline. This has been rigorously validated and consisted of ligand- and structure-based methodologies, which allowed us to select four promising CBS inhibitors called tubLCQF1-4. These four compounds were also evaluated with long trajectories molecular dynamics simulations and respective results were used for the theoretical determination of the free energy released in the formation of the complexes, using the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) methodology.