Toward Docking-Based Virtual Screening for Discovering Antitubulin Agents by Targeting Taxane and Colchicine Binding Sites

Discovery of microtubule stabilizers with novel scaffold structures based on virtual screening, biological evaluation, and molecular dynamics simulation

Disrupting the dynamics and structures of microtubules can perturb mitotic spindle formation, cause cell cycle arrest in G2/M phase, and subsequently lead to cellular death via apoptosis. In this investigation, the structure-based virtual screening methods, including molecular docking and rescoring, and similarity analysis of interaction molecular fingerprints, were developed to discover novel tubulin inhibitors from ChemDiv database with 1,601,806 compounds. The screened compounds were further filtered by PAINS, ADME/T, Toxscore, SAscore, and Drug-likeness analysis. Finally, 17 hit compounds were selected, and then submitted to the biologic evaluation. Among these hits, the P2 exhibited the strongest antiproliferative activity against four tumor cells including HeLa, HepG2, MCF-7, and A549. The in vitro tubulin polymerization assay revealed P2 could promote tubulin polymerization in a dose dependent manner. Finally, in order to analyze the interaction modes of complexes, the molecular dynamics simulation was performed to investigate the interactions between P2 and tubulin. The molecular dynamics simulation analysis showed that P2 could stably bind to taxane site, induced H6-H7, B9-B10, and M-loop regions changes. The ΔG_bind energies of tubulin-P2 and tubulin-paclitaxel were -68.25 ± 12.98 and -146.05 ± 16.17 kJ mol^-1, respectively, which were in line with the results of the experimental test. Therefore, P2 has been well characterized as lead compounds for developing new tubulin inhibitors with potential anticancer activity.

Identification of novel αβ-tubulin modulators with antiproliferative activity directed to cancer therapy using ligand and structure-based virtual screening

Among several strategies related to cancer therapy targeting the modulation of αβ-tubulin has shown encouraging findings, more specifically when this is achieved by inhibitors located at the colchicine binding site. In this work, we aim to fish new αβ-tubulin modulators through a diverse and rational VS study, and thus, exhibiting the development of two VS pipelines. This allowed us to identify two compounds 5 and 9 that showed IC₅₀ values of 19.69 and 21.97 μM, respectively, towards possible modulation of αβ-tubulin, such as assessed by in vitro assays in C6 glioma and HEPG2 cell lines. We also evaluated possible mechanisms of action of obtained hits towards the colchicine binding site of αβ-tubulin by using docking approaches. In addition, assessment of the stability of the active (5 and 9) and inactive compounds (3 and 13) within the colchicine binding site was carried out by molecular dynamics (MD) simulations, highlighting the solvent effect and revealing the compound 5 as the most stable in the complex. At last, deep analysis of these results provided some valuable insights on the importance of using mixed ligand- and structure-based strategies in VS campaigns, in order to achieve higher chemical diversity and biological effect as well.

Interdimeric Curvature in Tubulin-Tubulin Complexes Delineates the Microtubule-Destabilizing Properties of Plocabulin

Plocabulin is a novel microtubule (MT) destabilizer agent with potent antineoplastic activity. This compound binds to the maytansine site at the longitudinal interface between tubulin dimers and exerts a hinge-like effect that disrupts normal microtubule assembly. Plocabulin has emerged as a valuable model for the rational design of novel MT destabilizers because of its unique structural and mechanistic features. To make progress on this matter, detailed molecular-level understanding of the ligand-protein interactions responsible for plocabulin association and the conformation and energetic effects arising from plocabulin binding on the longitudinal interaction between tubulin dimers must be provided. In this work, fully atomistic MD simulations and MM/GBSA binding free-energy calculations were used to examine the association of plocabulin to one or two tubulin dimers in longitudinal arrangement. Our results revealed that plocabulin binding is favored by the addition of a second tubulin dimer and that this ligand promotes the assembly of curved tetrameric arrangements with strengthened longitudinal interdimeric interactions compared to ligand-free systems. The applicability of these findings to the rational discovery of novel MT destabilizers was tested using MD and MM/GBSA calculations as filtering tools to narrow the results of virtual screening among an FDA-approved drug database. Our results confirmed that tight-binding ligands do not necessarily exert the expected conformational and energetic effects on longitudinal tubulin-tubulin interactions, which is a matter to consider in future design strategies.

Docking and binding free energy calculations of sirtuin inhibitors

Sirtuins form a unique and highly conserved class of NAD(+)-dependent lysine deacylases. Among these the human subtypes Sirt1-3 has been implicated in the pathogenesis of numerous diseases such as cancer, metabolic syndromes, viral diseases and neurological disorders. Most of the sirtuin inhibitors that have been identified so far show limited potency and/or isoform selectivity. Here, we introduce a promising method to generate protein-inhibitor complexes of human Sirt1, Sirt2 and Sirt3 by means of ligand docking and molecular dynamics simulations. This method highly reduces the complexity of such applications and can be applied to other protein targets beside sirtuins. To the best of our knowledge, we present the first binding free energy method developed by using a validated data set of sirtuin inhibitors, where both a fair number of compounds (33 thieno[3,2-d]pyrimidine-6-carboxamide derivatives) was developed and tested in the same laboratory and also crystal structures in complex with the enzyme have been reported. A significant correlation between binding free energies derived from MM-GBSA calculations and in vitro data was found for all three sirtuin subtypes. The developed MM-GBSA protocol is computationally inexpensive and can be applied as a post-docking filter in virtual screening to find novel Sirt1-3 inhibitors as well as to prioritize compounds with similar chemical structures for further biological characterization.

Discovery of 7-hydroxy-6-methoxy-2-methyl-3-(3,4,5-trimethoxybenzoyl)benzo[b]furan (BNC105), a tubulin polymerization inhibitor with potent antiproliferative and tumor vascular disrupting properties

A structure-activity relationship (SAR) guided design of novel tubulin polymerization inhibitors has resulted in a series of benzo[b]furans with exceptional potency toward cancer cells and activated endothelial cells. The potency of early lead compounds has been substantially improved through the synergistic effect of introducing a conformational bias and additional hydrogen bond donor to the pharmacophore. Screening of a focused library of potent tubulin polymerization inhibitors for selectivity against cancer cells and activated endothelial cells over quiescent endothelial cells has afforded 7-hydroxy-6-methoxy-2-methyl-3-(3,4,5-trimethoxybenzoyl)benzo[b]furan (BNC105, 8) as a potent and selective antiproliferative. Because of poor solubility, 8 is administered as its disodium phosphate ester prodrug 9 (BNC105P), which is rapidly cleaved in vivo to return the active 8. 9 exhibits both superior vascular disrupting and tumor growth inhibitory properties compared with the benchmark agent combretastatin A-4 disodium phosphate 5 (CA4P).