A Molecular Docking Study Reveals That Short Peptides Induce Conformational Changes in the Structure of Human Tubulin Isotypes αβI, αβII, αβIII and αβIV

Abstract

Microtubules are cylindrical protein polymers assembled in the cytoplasm of all eukaryotic cells by polymerization of aβ tubulin dimers, which are involved in cell division, migration, signaling, and intracellular traffic. These functions make them essential in the proliferation of cancerous cells and metastases. Tubulin has been the molecular target of many anticancer drugs because of its crucial role in the cell proliferation process. By developing drug resistance, tumor cells severely limit the successful outcomes of cancer chemotherapy. Hence, overcoming drug resistance motivates the design of new anticancer therapeutics. Here, we retrieve short peptides obtained from the data repository of antimicrobial peptides (DRAMP) and report on the computational screening of their predicted tertiary structures for the ability to inhibit tubulin polymerization using multiple combinatorial docking programs, namely PATCHDOCK, FIREDOCK, and ClusPro. The interaction visualizations show that all the best peptides from the docking analysis bind to the interface residues of the tubulin isoforms αβl, αβll, αβlll, and αβlV, respectively. The docking studies were further confirmed by a molecular dynamics simulation, in which the computed root-mean-square deviation (RMSD), and root-mean-square fluctuation (RMSF), verified the stable nature of the peptide–tubulin complexes. Physiochemical toxicity and allergenicity studies were also performed. This present study suggests that these identified anticancer peptide molecules might destabilize the tubulin polymerization process and hence can be suitable candidates for novel drug development. It is concluded that wet-lab experiments are needed to validate these findings.

Development and Identification of Novel α-Glucosidase Inhibitory Peptides from Mulberry Leaves

“…Root mean square deviation (RMSD) represents the average displacement deviation in a specific simulation duration of the complex to reference, and a lower value indicates better stability of the system [39]. As demonstrated in Figure 6, the receptor protein in every docking system was relatively stable with an RMSD value less than 0.25 nm.…”