Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 21 February 2020
Sec. Theoretical and Computational Chemistry
This article is part of the Research Topic In Silico Methods for Drug Design and Discovery View all 35 articles

In silico Investigations of the Mode of Action of Novel Colchicine Derivatives Targeting β-Tubulin Isotypes: A Search for a Selective and Specific β-III Tubulin Ligand

  • 1PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
  • 2Department of Chemistry, Adam Mickiewicz University, Poznań, Poland
  • 3Dalle Molle Institute for Artificial Intelligence (IDSIA), University of Applied Sciences of Southern Switzerland (SUPSI), University of Italian Switzerland (USI), Manno, Switzerland
  • 4Department of Oncology, University of Alberta, Edmonton, AB, Canada

The cardinal role of microtubules in cell mitosis makes them interesting drug targets for many pharmacological treatments, including those against cancer. Moreover, different expression patterns between cell types for several tubulin isotypes represent a great opportunity to improve the selectivity and specificity of the employed drugs and to design novel compounds with higher activity only on cells of interest. In this context, tubulin isotype βIII represents an excellent target for anti-tumoral therapies since it is overexpressed in most cancer cells and correlated with drug resistance. Colchicine is a well-known antimitotic agent, which is able to bind the tubulin dimer and to halt the mitotic process. However, it shows high toxicity also on normal cells and it is not specific for isotype βIII. In this context, the search for colchicine derivatives is a matter of great importance in cancer research. In this study, homology modeling techniques, molecular docking, and molecular dynamics simulations have been employed to characterize the interaction between 55 new promising colchicine derivatives and tubulin isotype βIII. These compounds were screened and ranked based on their binding affinity and conformational stability in the colchicine binding site of tubulin βIII. Results from this study point the attention on an amide of 4-chlorine thiocolchicine. This colchicine-derivative is characterized by a unique mode of interaction with tubulin, compared to all other compounds considered, which is primarily characterized by the involvement of the α-T5 loop, a key player in the colchicine binding site. Information provided by the present study may be particularly important in the rational design of colchicine-derivatives targeting drug resistant cancer phenotypes.

Introduction

The pivotal role of microtubules (MTs) in the mitotic process make them important targets for anticancer therapies since cancerous cells proliferate by unregulated cell division (Gajewski et al., 2013). By either stabilizing MTs or enhancing their depolymerization, it is possible to halt the mitotic process and eventually lead cells to apoptosis (Nettles et al., 2002). Among antimitotic agents, colchicine is able to block cell division (Bhattacharyya et al., 2008) by destabilizing MT assembly kinetics and dynamics. In particular, when colchicine binds in its specific binding site (located at the interface between tubulin α and β monomers) the structural conformation of the tubulin dimer is affected in such a way that tubulin integration into the MT lattice is hampered.

However, one of the main drawbacks of colchicine is its general toxicity (Wallace, 1974; Finkelstein et al., 2010). Several studies in the past have proposed less toxic colchicine derivatives as an alternative to colchicine (Lu et al., 2012; Wang et al., 2016; Johnson et al., 2017; Majcher et al., 2018a,b; Klejborowska et al., 2019). Moreover, these novel colchicine derivatives may be designed to show high specificity only for tubulin isotypes, which are over-expressed in cancer, in order to maximize their effect only on tumor cells and reduce side effects of the drug due to its toxicity on normal cells (Lu and Luduena, 1994; Luduena et al., 1995).

Differing in point or restricted sequence variations, several tubulin isotypes (Leandro-García et al., 2010) are differently expressed by cells under both physiological and pathological conditions. For example, tubulin isotype αβIII is considered as an excellent target for anti-tumoral therapies because it is over-expressed in tumoral cells and it is less widespread than other isotypes, such as αβI, αβII and αβIV, in normal cells (Ferlini et al., 2007; Tseng et al., 2010). Moreover, an over-expression of tubulin isotype αβIII by cancer cells is considered as one among several known drug resistance mechanisms (Derry et al., 1997; Ludueña, 1998; Katsetos et al., 2003; Kamath et al., 2005; Seve, 2005; Ferlini et al., 2007; Sève and Dumontet, 2008; Tseng et al., 2010). Thereby, it is of primary importance to identify specific compounds, which selectively target isotype αβIII.

In this context, computational molecular modeling techniques, such as molecular dynamics (MD) and molecular docking, represent powerful tools to shed light on the molecular mechanisms concerning protein functions and their interaction between different ligands and a specific receptor (Lepre et al., 2017; Omar et al., 2018; Brogi, 2019; Sirous et al., 2019). These computational methods can be applied to investigate the action of different ligands on tubulin dimers (Mitra and Sept, 2008; Natarajan and Senapati, 2012; Gajewski et al., 2013; Kumbhar et al., 2016). Computational drug discovery may help to accelerate and economize the drug discovery process as a complementary tool for experimental research of novel inhibitors.

In this work, ensemble molecular docking, molecular dynamics simulations, and binding energy estimation methods have been employed to characterize the binding of 55 novel colchicine derivatives to the βIII tubulin isotype. We have identified an interesting 4-chlorine thiocolchicine derivative characterized by similar affinity but a different mode of binding to tubulin with respect to its parent compound, colchicine. The main findings of our study indicate this ligand as a promising candidate to overcome colchicine drawbacks and provide information for further developments in designing more selected and specific colchicine derivatives with an intended use as cancer chemotherapy agents.

Materials and Methods

Atomic Models of Investigated Compounds

Several series of novel colchicine derivatives (Majcher et al., 2018a,b; Klejborowska et al., 2019) were considered in this work. All 55 compounds have shown in vitro anti-proliferative effects on normal and cancer cells. In particular, they were tested on human lung adenocarcinoma, human breast adenocarcinoma, human colon adenocarcinoma cell lines and a doxorubicin-resistant subline (Majcher et al., 2018a,b; Klejborowska et al., 2019).

These compounds can be divided into five classes: 4-Br-Amides (10 compounds), 4-Cl-Amides (10 compounds), DT-and-4I-Amides (19 compounds), 4-Cl-Carbamates (8 compounds) and 4-I-Carbamates (8 compounds). The chemical structures of colchicine (C01) and its derivatives (C02-C56) are summarized in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Colchicine and its derivatives considered in this work.

The 2D structures of the colchicine derivatives have been drawn using ChemDraw 12.0, whereas their 3D structure was designed by AVOGADRO (Hanwell et al., 2012).

Human αβIII Tubulin Modeling and Conformational Dynamics

The atomic structures of human βIII tubulin isotype were obtained by homology modeling, starting from the Protein Data Bank (PDB) entry 4O2B model (Prota et al., 2014) as a template. This structure was chosen due to its high resolution (2.3 Å) and a low number of missing residues (Aryapour et al., 2017). First, from the starting template the information concerning tubulin αβ, GTP, GDP, Mg2+ ion and colchicine was extracted. Missing residues in β tubulin (from 276 to 281) were added by MODELER 9.20 (Šali and Blundell, 1993) where the best model was selected on the basis of the obtained DOPE (Discrete optimized protein energy) score. Then, the Fasta sequences Q71U36 and Q13509 were selected from the Uniprot website, respectively, for the α and β subunits. The above-mentioned amino acid sequences pertain to the isotype αβIII (Gajewski et al., 2013; Kumbhar et al., 2016). Homology modeling was then employed by MODELER 9.20 to generate a 3D structure of the αβIII sequence using the 4O2B model. The quality and the reliability of the generated model were evaluated using PROCHECK (Laskowski et al., 1993), VERIFY3D (Colovos and Yeates, 1993) and ERRAT (Bowie et al., 1991), as reported in previous literature in this area (Huzil et al., 2006; Deriu et al., 2007; Mane et al., 2008; Kumbhar et al., 2016).

Two systems were subsequently considered: (I) tubulin, GTP, GDP, and Mg2+ ion and (II) tubulin, GTP, GDP, Mg2+ ion and colchicine bound to tubulin. Information on colchicine binding was taken from the 4O2B model. The AMBER ff99SB-ILDN forcefield (Lindorff-Larsen et al., 2010) was used to describe protein, water and ion topology. GTP, GDP, and ligands were described by the General Amber Force Field (GAFF) (Wang et al., 2004) and AM1-BCC charge method (Jakalian et al., 2002), as applied in many previous studies (Gajewski et al., 2013; Kumbhar et al., 2016; Klejborowska et al., 2019; Sahakyan et al., 2019).

Each of the above mentioned protein systems (I and II) was then inserted in a dodecahedron box filled with TIP3P explicit water molecules (Mark and Nilsson, 2001), sodium and chlorine ions (150 mM). Particle-mesh Ewald (PME) method (Darden et al., 1993) was used to treat electrostatics (cut-off = 1.0 nm) whereas Van der Waals (VdW) interactions were treated by a plain cut-off at 1 nm (Natarajan and Senapati, 2012; Natarajan et al., 2013; Bueno et al., 2018). Each system was then energy minimized by the steepest descent algorithm for 1,000 steps with a maximum force of 100 kJmol−1nm−1. All systems were simulated in an NVT and NPT ensemble with position restraints applied on protein and ligand atoms. In detail, a 100 ps position restrained MD simulation in the NVT ensemble (Bussi et al., 2007), was followed by a 300 ps position restrained MD in the NPT (T = 300 K, P = 1 bar) ensemble (Berendsen et al., 1984; Bussi et al., 2007).

Temperature and pressure were controlled by weak coupling algorithms (Berendsen et al., 1984; Bussi et al., 2007).

Finally, production MD simulations (without restraints) were carried out for 100 ns in presence and in absence of colchicine, respectively. Ten configurations of each system were extracted as representative of structural equilibrium. The above-mentioned system configurations were then used for ensemble docking procedure.

Ensemble Docking and Binding Energy Refinement

Ensemble docking was performed using AUTODOCK VINA 1.1.2 (Trott and Olson, 2010). The center of the search space was defined by taking, from the 4O2B model, the relative position of the colchicine in its binding site. The docking was performed using a grid space 2 × 2 × 2 nm around the center of the search space and an exhaustiveness equal to 64 was set. Each compound was docked to the ten different isotype configurations extracted from the production MD, as explained above.

Then, for each VINA pose, the binding energy refinement was performed by running short 1 ns MD simulations on the ligand-protein complex starting from the VINA best pose for each considered ligand. Each complex was followed by solvation, neutralization, energy minimization, position restrained MD, and short production MD. Simulation set up was the same as described in the previous section. On the last 100 ps of MD production the ligand-protein binding was evaluated using two criteria. Firstly, the binding energy was quantified by the Molecular Mechanics Generalized Born Surface Area (MMGBSA) method (Genheden and Ryde, 2015). The parameters were set according to the previous literature (Nguyen et al., 2013, 2015; Su et al., 2015). Secondly, the ligand conformational displacement in the binding site was quantified by calculating the Root Mean Square Deviation (RMSD) of ligand carbon rings (a common feature of all considered compounds with colchicine). In particular, for each ligand, the MD protein trajectory was fitted on a reference structure (the starting configuration of the colchicine-protein complex). In this way, the RMSD quantifies the relative deviation of each ligand with respect to the colchicine starting position throughout the overall MD trajectory. Based on the above-mentioned criteria, the best colchicine derivative and colchicine, both bound to the βIII isotype, were simulated for 100 ns in order to highlight binding conformational differences at equilibrium.

All MD simulations were carried out using GROMACS 2018.3 (Abraham et al., 2015). The Visual Molecular Dynamics (VMD) package was employed for the visual inspection of the simulated systems (Humphrey et al., 1996). Dedicated GROMACS tools were used for a quantitative analysis in terms of Root-Mean-Square Deviation (RMSD), Root-Mean-Square Fluctuation (RMSF), and clustering, while analysis of the secondary structure was performed by applying the STRuctural IDEntification (STRIDE) algorithm (Heinig and Frishman, 2004).

Results

Human αβIII Tubulin Model Development and Conformational Dynamics

The Ramachandran plot (see also Figure S1) obtained by PROCHECK highlighted the 95.6% of residues in most favored regions, 4.2% in additional allowed regions, and 0.1% in generously allowed regions. No residues were found in disallowed regions. Since a good quality model is expected to have at least 90% of the residues in the most favored regions (Santoshi and Naik, 2014), the built model was considered reliable. Moreover, the Overall Quality Factor obtained by the ERRAT tool for the isotype αβIII was 80.29 for the α and 84.73 for the β tubulin monomer model. It is worth mentioning that the generally accepted range is higher than 50 for a high quality model (Messaoudi et al., 2013). Finally, the VERIFY3D confirmed that 98.15% of residues showed an averaged 3D-1D score higher than 0.2 (Messaoudi et al., 2013).

First, the backbone RMSD was calculated for isotype βIII both in presence and in absence of colchicine during the overall MD simulation (100 ns): all the simulated structures reached structural equilibrium, with values under 0.3 nm (see also Figure S2). Moreover, the cluster analysis on the last 50 ns of the simulations highlighted only one cluster using an RMSD cut-off of 0.15 nm, indicating a strong stability of the simulated systems. Moreover, the cluster analysis indicated that the colchicine presence did not modify significantly the conformation of the interaction site.

Ensemble Docking and Binding Energy Calculation

The 55 colchicine derivatives were docked to ten different configurations of βIII tubulin, extracted from the last 50 ns of the MD simulation described above. Only the best ligand pose in terms of binding affinity was considered (see also Figure S3). In order to take into account also the dynamic nature of the binding process, we have performed a MD simulation of 1 ns for each ligand-receptor complex. Throughout the quick MD run, the binding energy was quantified by means of the MM-GBSA method (Huzil et al., 2010; Gajewski et al., 2013; Kumbhar et al., 2016). Moreover, the ligand displacement in the binding site was quantified by the RMSD calculated as described in Materials and Methods. It is worth mentioning that low RMSD values indicate a compound which is stable in a spot close to the starting colchicine position, whereas high RMSD values identify a compound moving further apart (Figure 2). Most compounds showed RMSD lower than 0.2 nm, suggesting that the derivatives investigated here behaves similarly to colchicine (highly stable in its binding site during the short MD run). The only exception found is represented by compound C19 which displays high variation from the colchicine starting position (RMSD = 0.47 nm).

FIGURE 2
www.frontiersin.org

Figure 2. Representative snapshots of ligand conformational displacements in the colchicine binding site. Colchicine is represented in green, whereas two different derivatives with low (A) and high (B) RMSD with respect to the colchicine starting pose are depicted orange and red, respectively.

Our binding energy analysis highlights four specific compounds (C19, C20, C29, and C48) as possible hits. In fact, they exhibit similar values of their binding energy for βIII tubulin compared to colchicine. All binding energy values are reported in Supporting Information text (see also Figure S4). In order to better describe differences between investigated compounds, we have merged RMSD and binding energy information in a single plot (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Polar scatter plot representing ligands' RMSD from the colchicine starting pose (radial coordinate) in nm and their binding energy (angular coordinate) in kJ/mol; the diamond marker represents the C19 compound.

Interestingly, compound C19 features a peculiar behavior, i.e., it exhibits a large deviation from the colchicine starting position in the binding site (high RMSD value), with a significant difference to all the other compounds. It is also characterized by binding energy values comparable to colchicine and higher than most other derivatives. This result points the attention on the compound C19 as a promising candidate able to bind strongly to βIII human tubulin with a different mode of action with respect to colchicine.

Conformational Dynamics of Colchicine and C19 Bound to βIII Human Tubulin

Conformational dynamics of the C19-tubulin complex has been investigated by a 100 ns long MD simulation. For comparison, a 100 ns-long MD was also carried out on the colchicine-tubulin complex. Systems were replicated to confirm the consistency of the data (Figures S5, S7).

Structural modifications of the colchicine binding site were first analyzed by computing the RMSD of the tubulin binding cleft, i.e., residues within 1 nm from the ligand, from its starting position and the secondary structure probability during the last 50 ns of the simulation (Figure 4). A second replica was performed to ensure the repeatability of the results (see also Figure S5). The binding site was characterized by a structural stability throughout the overall MD, exhibiting low RMSD values (lower than 0.18 nm) and highly conserved secondary structures. The only noteworthy difference is represented by the α-T5 loop, which exhibits tendency to rearrange in a more structured shape only in the presence of C19.

FIGURE 4
www.frontiersin.org

Figure 4. (A) RMSD of the colchicine binding site from its starting position when colchicine (black) or compound C19 (gray) are bound to the tubulin dimer. (B) Secondary structure probability of residues in the colchicine binding site when colchicine (B1) or compound C19 (B2) are bound to the tubulin dimer.

The binding energy was estimated by the MM-GBSA approach in order to compare the binding affinities of the analyzed compounds at the structural equilibrium: again, compound C19 and colchicine showed similar binding energy for isotype αβIII, respectively, −229.98 ± 22.26 kJ/mol and −223.70 ± 22.31 kJ/mol (see also Figure S6). Nevertheless, the energy decomposition over the tubulin binding cleft residues reveals that the compound C19 shows a higher binding energy compared to the colchicine for residues 178–180 of the α tubulin, which belong to the αT5 loop (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Binding Energy decomposition over the residues of the colchicine binding site (only residues with the highest energies are reported). Compound C19 has a significantly higher affinity than colchicine for the αT5 loop.

In light of these results, the ligands' behavior in the binding site and their interaction with the αT5 loop were investigated in more depth (Figure 6). First, ligand RMSD (used to quantify the ligand movement in the binding site throughout the simulation) showed that compound C19 has a more marked tendency than colchicine to move apart, reaching a more favorable pose for the interaction with the αT5 loop (Figure 6A and see also Figure S7). Second, the interaction surface between each ligand and the αT5 loop, which quantifies the available area for their binding, is higher for C19 than colchicine (Figure 6B). Figures 6C1,C2 and 6D represent ligand structures and their relative position in the tubulin binding cleft (see also Movies S1, S2).

FIGURE 6
www.frontiersin.org

Figure 6. (A) Ligands' RMSD from their starting position (colchicine in black, C19 in red). (B) Probability density function of the buried surface between the ligands and the αT5 loop (colchicine in black, C19 in red), averaged between two replicas during the last 20 ns of simulation. (C) Chemical structures of colchicine (C1) and compound C19 (C2). (D) representative snapshot of the simulation, which shows that compound C19 (red) is closer to αT5 loop (yellow) than colchicine (green).

In conclusion, compound C19 was shown to be stable in the tubulin binding site with a relative position differing from the colchicine site. Specifically, C19 is predicted to be mostly stabilized by its interaction with the αT5 loop.

Discussion

In this study, 55 colchicine derivatives were screened for their binding properties to tubulin isotype βIII. The research work was aimed at identifying alternative compounds able to overcome colchicine's well-known limitations. After the docking of all compounds to the target isotype of tubulin, a molecular dynamics simulation of 1 ns was performed on each generated receptor-ligand complex. The obtained trajectories were analyzed considering the deviations of the compounds from the colchicine's starting pose, using the RMSD, and the binding energy evaluated with the MM-GBSA method. All compounds were characterized by low RMSD values, except for compound C19, which showed high deviations (RMSD = 0.47 nm). This evidence suggests a different particular pose for this derivative. From the affinity analysis we found out that the binding energies for compounds C19, C20, C29, and C48 are similar to that for colchicine and higher than those found for most other ligands. These results indicate that C19 is a promising compound to be further investigated and experimentally validated. Its specific binding to tubulin is characterized by a different conformational organization and dynamics in the tubulin binding site with high affinity. RMSD analysis indicates that C19 is able to be accommodated in the binding site by moving toward more favorable poses for interaction with the αT5 loop. This feature is less pronounced by colchicine. Moreover, the buried surface between C19 and the tubulin isotype βIII, which measures the available area for the binding, is greater than the one exhibited with colchicine, confirming a higher stability of C19 in the binding site. Finally, the ligand binding to the αT5 loop may affect its secondary structure toward a more structured arrangement. Therefore, a compound able to influence the αT5 loop structure could affect the dynamics of the entire microtubule.

The above mentioned evidences might be of a significant interest given that the αT5 loop is a key player region in the colchicine binding site and for intra-dimer contacts (Ravelli et al., 2004). Nonetheless, previous literature (Bueno et al., 2018) already highlighted the importance of the αT5 loop, identified as relevant for the binding of a promising anti-proliferative compound (Bueno et al., 2018).

In conclusion, our study clarifies some features characterizing the βIII tubulin binding mode of a promising novel 4-chlorine thiocolchicine derivative, which differs profoundly from that known for colchicine. The specific interaction of compound C19 with the αT5 loop is a promising feature that could be related to an increased destabilizing activity of the ligand with respect to microtubule dynamics. Moreover, this unique behavior exhibited in complex with the βIII tubulin isotype is of primary importance since this isotype is overexpressed in cancer cells, while very insignificantly represented in most normal cells and also implicated in drug resistance (Katsetos et al., 2003; Kamath et al., 2005; Seve, 2005; Sève and Dumontet, 2008; Leandro-García et al., 2010). In light of these results, C19 or similar compounds, as promising candidates able to possibly overcome some colchicine's drawbacks, deserve further investigations, including biological toxicity assessment and cancer cell cytotoxicity experiments to prove its specificity and selectivity for βIII isotype of tubulin.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

MD, JT, and AH conceived the research. LP, AR, and GG did the molecular dynamics simulations. LP, AR, GG, and GK analyzed and rationalized the data. All authors wrote the paper and critically commented to the manuscript, read, and approved the final manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Financial support by grant of the Polish National Science Centre (NCN)–No. 2016/21/B/ST5/00111 is gratefully acknowledged. GK gratefully acknowledges the doctoral scholarship ETIUDA–No. 2018/28/T/ST5/00041 financed by the Polish National Science Centre (NCN).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00108/full#supplementary-material

References

Abraham, M. J., Murtola, T., Schulz, R., Páll, S., Smith, J. C., Hess, B., et al. (2015). GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25. doi: 10.1016/j.softx.2015.06.001

CrossRef Full Text | Google Scholar

Aryapour, H., Dehdab, M., Sohraby, F., and Bargahi, A. (2017). Prediction of new chromene-based inhibitors of tubulin using structure-based virtual screening and molecular dynamics simulation methods. Comput. Biol. Chem. 71, 89–97. doi: 10.1016/j.compbiolchem.2017.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A., and Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690. doi: 10.1063/1.448118

CrossRef Full Text | Google Scholar

Bhattacharyya, B., Panda, D., Gupta, S., and Banerjee, M. (2008). Anti-mitotic activity of colchicine and the structural basis for its interaction with tubulin. Med. Res. Rev. 28, 155–183. doi: 10.1002/med.20097

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowie, J., Luthy, R., and Eisenberg, D. (1991). A method to identify protein sequences that fold into a known three-dimensional structure. Science 253, 164–170. doi: 10.1126/science.1853201

PubMed Abstract | CrossRef Full Text | Google Scholar

Brogi, S. (2019). Computational approaches for drug discovery. Molecules 24:3061. doi: 10.3390/molecules24173061

PubMed Abstract | CrossRef Full Text | Google Scholar

Bueno, O., Estévez Gallego, J., Martins, S., Prota, A. E., Gago, F., Gómez-SanJuan, A., et al. (2018). High-affinity ligands of the colchicine domain in tubulin based on a structure-guided design. Sci. Rep. 8:4242. doi: 10.1038/s41598-018-22382-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bussi, G., Donadio, D., and Parrinello, M. (2007). Canonical sampling through velocity rescaling. J. Chem. Phys. 126:014101. doi: 10.1063/1.2408420

PubMed Abstract | CrossRef Full Text | Google Scholar

Colovos, C., and Yeates, T. O. (1993). Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 2, 1511–1519. doi: 10.1002/pro.5560020916

PubMed Abstract | CrossRef Full Text | Google Scholar

Darden, T., York, D., and Pedersen, L. (1993). Particle mesh ewald: an N log(N) method for ewald sums in large systems. J. Chem. Phys. 98:10089. doi: 10.1063/1.464397

CrossRef Full Text | Google Scholar

Deriu, M. A., Enemark, S., Soncini, M., Montevecchi, F. M., and Redaelli, A. (2007). Tubulin: from atomistic structure to supramolecular mechanical properties. J. Mater. Sci. 42, 8864–8872. doi: 10.1007/s10853-007-1784-6

CrossRef Full Text | Google Scholar

Derry, W. B., Wilson, L., Khan, I. A., Ludueña, R. F., and Jordan, M. A. (1997). Taxol differentially modulates the dynamics of microtubules assembled from unfractionated and purified β-tubulin isotypes. Biochemistry 36, 3554–3562. doi: 10.1021/bi962724m

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferlini, C., Raspaglio, G., Cicchillitti, L., Mozzetti, S., Prislei, S., Bartollino, S., et al. (2007). Looking at drug resistance mechanisms for microtubule interacting drugs: does TUBB3 work? Curr. Cancer Drug Targets 7, 704–712. doi: 10.2174/156800907783220453

PubMed Abstract | CrossRef Full Text | Google Scholar

Finkelstein, Y., Aks, S. E., Hutson, J. R., Juurlink, D. N., Nguyen, P., Dubnov-Raz, G., et al. (2010). Colchicine poisoning: the dark side of an ancient drug. Clin. Toxicol. 48, 407–414. doi: 10.3109/15563650.2010.495348

PubMed Abstract | CrossRef Full Text | Google Scholar

Gajewski, M. M., Tuszynski, J. A., Barakat, K., Torin Huzil, J., and Klobukowski, M. (2013). Interactions of laulimalide, peloruside, and their derivatives with the isoforms of β-tubulin. Can. J. Chem. 91, 511–517. doi: 10.1139/cjc-2012-0360

CrossRef Full Text | Google Scholar

Genheden, S., and Ryde, U. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Exp. Opin. Drug Discov. 10, 449–461. doi: 10.1517/17460441.2015.1032936

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanwell, M. D., Curtis, D. E., Lonie, D. C., Vandermeersch, T., Zurek, E., and Hutchison, G. R. (2012). Avogadro: an advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 4:17. doi: 10.1186/1758-2946-4-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Heinig, M., and Frishman, D. (2004). STRIDE: a web server for secondary structure assignment from known atomic coordinates of proteins. Nucleic Acids Res. 32, W500–W502. doi: 10.1093/nar/gkh429

PubMed Abstract | CrossRef Full Text | Google Scholar

Humphrey, W., Dalke, A., and Schulten, K. (1996). VMD: visual molecular dynamics. J. Mol. Graph. 14, 33–38. doi: 10.1016/0263-7855(96)00018-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Huzil, J. T., Ludueña, R. F., and Tuszynski, J. (2006). Comparative modelling of human β tubulin isotypes and implications for drug binding. Nanotechnology 17, S90–S100. doi: 10.1088/0957-4484/17/4/014

CrossRef Full Text | Google Scholar

Huzil, J. T., Mane, J., and Tuszynski, J. A. (2010). Computer assisted design of second-generation colchicine derivatives. Interdiscip. Sci. Comput. Life Sci. 2, 169–174. doi: 10.1007/s12539-010-0076-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakalian, A., Jack, D. B., and Bayly, C. I. (2002). Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. parameterization and validation. J. Comput. Chem. 23, 1623–1641. doi: 10.1002/jcc.10128

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, L., Goping, I. S., Rieger, A., Mane, J. Y., Huzil, T., Banerjee, A., et al. (2017). Novel colchicine derivatives and their anti-cancer activity. Curr. Top. Med. Chem. 17, 2538–2558. doi: 10.2174/1568026617666170104143618

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamath, K., Wilson, L., Cabral, F., and Jordan, M. A. (2005). βIII-tubulin induces paclitaxel resistance in association with reduced effects on microtubule dynamic instability. J. Biol. Chem. 280, 12902–12907. doi: 10.1074/jbc.M414477200

PubMed Abstract | CrossRef Full Text | Google Scholar

Katsetos, C. D., Herman, M. M., and Mörk, S. J. (2003). Class III β-tubulin in human development and cancer. Cell Motil. Cytoskelet. 55, 77–96. doi: 10.1002/cm.10116

CrossRef Full Text | Google Scholar

Klejborowska, G., Urbaniak, A., Preto, J., Maj, E., Moshari, M., Wietrzyk, J., et al. (2019). Synthesis, biological evaluation and molecular docking studies of new amides of 4-bromothiocolchicine as anticancer agents. Bioorg. Med. Chem. 76:115144. doi: 10.1016/j.bmc.2019.115144

CrossRef Full Text | Google Scholar

Kumbhar, B. V., Borogaon, A., Panda, D., and Kunwar, A. (2016). Exploring the origin of differential binding affinities of human tubulin isotypes αβII, αβIII and αβIV for DAMA-colchicine using homology modelling, molecular docking and molecular dynamics simulations. PLoS ONE 11:e0156048. doi: 10.1371/journal.pone.0156048

PubMed Abstract | CrossRef Full Text | Google Scholar

Laskowski, R. A., MacArthur, M. W., Moss, D. S., and Thornton, J. M. (1993). PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26, 283–291. doi: 10.1107/S0021889892009944

CrossRef Full Text | Google Scholar

Leandro-García, L. J., Leskelä, S., Landa, I., Montero-Conde, C., López-Jiménez, E., Letón, R., et al. (2010). Tumoral and tissue-specific expression of the major human β-tubulin isotypes. Cytoskeleton 67, 214–223. doi: 10.1002/cm.20436

PubMed Abstract | CrossRef Full Text | Google Scholar

Lepre, M., Omar, S., Grasso, G., Morbiducci, U., Deriu, M., and Tuszynski, J. (2017). Insights into the effect of the G245S single point mutation on the structure of p53 and the binding of the protein to DNA. Molecules 22:E1358. doi: 10.3390/molecules22081358

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindorff-Larsen, K., Piana, S., Palmo, K., Maragakis, P., Klepeis, J. L., Dror, R. O., et al. (2010). Improved side-chain torsion potentials for the amber ff99SB protein force field. Proteins 78, 1950–1958. doi: 10.1002/prot.22711

PubMed Abstract | CrossRef Full Text

Lu, Q., and Luduena, R. F. (1994). In vitro analysis of microtubule assembly of isotypically pure tubulin dimers: intrinsic differences in the assembly properties of alpha beta II, alpha beta III, and alpha beta IV tubulin dimers in the absence of microtubule-associated proteins. J. Biol. Chem. 269, 2041–2047.

PubMed Abstract | Google Scholar

Lu, Y., Chen, J., Xiao, M., Li, W., and Miller, D. D. (2012). An overview of tubulin inhibitors that interact with the colchicine binding site. Pharm. Res. 29, 2943–2971. doi: 10.1007/s11095-012-0828-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludueña, R. F. (1998). Multiple forms of tubulin: different gene products and covalent modifications. Int. Rev. Cytol.178, 207–275. doi: 10.1016/S0074-7696(08)62138-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Luduena, R. F., Roach, M. C., Prasad, V., Chaudhuri, A. R., Tomita, I., Mizuhashi, F., et al. (1995). Interaction of bovine brain tubulin with the 4(1H)-pyrizinone derivative IKP104, an antimitotic drug with a complex set of effects on the conformational stability of the tubulin molecule. Biochemistry 34, 15751–15759. doi: 10.1021/bi00048a020

PubMed Abstract | CrossRef Full Text | Google Scholar

Majcher, U., Klejborowska, G., Moshari, M., Maj, E., Wietrzyk, J., Bartl, F., et al. (2018a). Antiproliferative activity and molecular docking of novel double-modified colchicine derivatives. Cells 7:192. doi: 10.3390/cells7110192

PubMed Abstract | CrossRef Full Text | Google Scholar

Majcher, U., Urbaniak, A., Maj, E., Moshari, M., Delgado, M., Wietrzyk, J., et al. (2018b). Synthesis, antiproliferative activity and molecular docking of thiocolchicine urethanes. Bioorg. Chem. 81, 553–566. doi: 10.1016/j.bioorg.2018.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Mane, J. Y., Klobukowski, M., Huzil, J. T., and Tuszynski, J. (2008). Free energy calculations on the binding of colchicine and its derivatives with the α/β-tubulin isoforms. J. Chem. Inf. Model. 48, 1824–1832. doi: 10.1021/ci800054n

PubMed Abstract | CrossRef Full Text | Google Scholar

Mark, P., and Nilsson, L. (2001). Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J. Phys. Chem. A 105, 9954–9960. doi: 10.1021/jp003020w

CrossRef Full Text | Google Scholar

Messaoudi, A., Belguith, H., and Ben Hamida, J. (2013). Homology modeling and virtual screening approaches to identify potent inhibitors of VEB-1 β-lactamase. Theor. Biol. Med. Model. 10, 1–10. doi: 10.1186/1742-4682-10-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitra, A., and Sept, D. (2008). Taxol allosterically alters the dynamics of the tubulin dimer and increases the flexibility of microtubules. Biophys. J. 95, 3252–3258. doi: 10.1529/biophysj.108.133884

PubMed Abstract | CrossRef Full Text | Google Scholar

Natarajan, K., Mohan, J., and Senapati, S. (2013). Relating nucleotide-dependent conformational changes in free tubulin dimers to tubulin assembly. Biopolymers 99, 282–291. doi: 10.1002/bip.22153

PubMed Abstract | CrossRef Full Text | Google Scholar

Natarajan, K., and Senapati, S. (2012). Understanding the basis of drug resistance of the mutants of αβ-tubulin dimer via molecular dynamics simulations. PLoS ONE 7:e42351. doi: 10.1371/journal.pone.0042351

PubMed Abstract | CrossRef Full Text | Google Scholar

Nettles, J. H., Snyder, J. P., Nogales, E., Cornett, B., and Downing, K. H. (2002). The binding conformation of taxol in -tubulin: a model based on electron crystallographic density. Proc. Natl. Acad. Sci. U.S.A. 98, 5312–5316. doi: 10.1073/pnas.051309398

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, H., Pérez, A., Bermeo, S., and Simmerling, C. (2015). Refinement of generalized born implicit solvation parameters for nucleic acids and their complexes with proteins. J. Chem. Theory Comput. 11, 3714–3728. doi: 10.1021/acs.jctc.5b00271

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, H., Roe, D. R., and Simmerling, C. (2013). Improved generalized born solvent model parameters for protein simulations. J. Chem. Theory Comput. 9, 2020–2034. doi: 10.1021/ct3010485

PubMed Abstract | CrossRef Full Text | Google Scholar

Omar, S. I., Lepre, M. G., Morbiducci, U., Deriu, M. A., and Tuszynski, J. A. (2018). Virtual screening using covalent docking to find activators for G245S mutant p53. PLoS ONE 13:e0200769. doi: 10.1371/journal.pone.0200769

PubMed Abstract | CrossRef Full Text | Google Scholar

Prota, A. E., Danel, F., Bachmann, F., Bargsten, K., Buey, R. M., Pohlmann, J., et al. (2014). The novel microtubule-destabilizing drug BAL27862 binds to the colchicine site of tubulin with distinct effects on microtubule organization. J. Mol. Biol. 426, 1848–1860. doi: 10.1016/j.jmb.2014.02.005

CrossRef Full Text | Google Scholar

Ravelli, R. B. G., Gigant, B., Curmi, P. A., Jourdain, I., Lachkar, S., Sobel, A., et al. (2004). Insight into tubulin regulation from a complex with colchicine and a stathmin-like domain. Nature 428, 198–202. doi: 10.1038/nature02393

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahakyan, H., Abelyan, N., Arakelov, V., Arakelov, G., and Nazaryan, K. (2019). In silico study of colchicine resistance molecular mechanisms caused by tubulin structural polymorphism. PLoS ONE 14:e0221532. doi: 10.1371/journal.pone.0221532

PubMed Abstract | CrossRef Full Text | Google Scholar

Šali, A., and Blundell, T. L. (1993). Comparative protein modelling by satisfaction of Spatial restraints. J. Mol. Biol. 234, 779–815. doi: 10.1006/jmbi.1993.1626

PubMed Abstract | CrossRef Full Text | Google Scholar

Santoshi, S., and Naik, P. K. (2014). Molecular insight of isotypes specific β-tubulin interaction of tubulin heterodimer with noscapinoids. J. Comput. Aided Mol. Des. 28, 751–763. doi: 10.1007/s10822-014-9756-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Seve, P. (2005). Class III -tubulin expression in tumor cells predicts response and outcome in patients with non-small cell lung cancer receiving paclitaxel. Mol. Cancer Ther. 4, 2001–2007. doi: 10.1158/1535-7163.MCT-05-0244

PubMed Abstract | CrossRef Full Text | Google Scholar

Sève, P., and Dumontet, C. (2008). Is class III β-tubulin a predictive factor in patients receiving tubulin-binding agents? Lancet Oncol. 9, 168–175. doi: 10.1016/S1470-2045(08)70029-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sirous, H., Chemi, G., Campiani, G., and Brogi, S. (2019). An integrated in silico screening strategy for identifying promising disruptors of p53-MDM2 interaction. Comput. Biol. Chem. 83:107105. doi: 10.1016/j.compbiolchem.2019.107105

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, P. C., Tsai, C. C., Mehboob, S., Hevener, K. E., and Johnson, M. E. (2015). Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F. tularensis enoyl-ACP reductase (FabI). J. Comput. Chem. 36, 1859–1873. doi: 10.1002/jcc.24011

PubMed Abstract | CrossRef Full Text | Google Scholar

Trott, O., and Olson, A. J. (2010). AutoDock vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455–461. doi: 10.1002/jcc.21334

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, C.-Y., Mane, J. Y., Winter, P., Johnson, L., Huzil, T., Izbicka, E., et al. (2010). Quantitative analysis of the effect of tubulin isotype expression on sensitivity of cancer cell lines to a set of novel colchicine derivatives. Mol. Cancer 9:131. doi: 10.1186/1476-4598-9-131

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallace, S. L. (1974). Colchicine. Semin. Arthritis Rheum. 3, 369–381. doi: 10.1016/0049-0172(74)90006-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2004). Development and testing of a general amber force field. J. Comput. Chem. 25, 1157–1174. doi: 10.1002/jcc.20035

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhang, H., Gigant, B., Yu, Y., Wu, Y., Chen, X., et al. (2016). Structures of a diverse set of colchicine binding site inhibitors in complex with tubulin provide a rationale for drug discovery. FEBS J. 283, 102–111. doi: 10.1111/febs.13555

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: molecular modeling, drug discovery, microtubule, cancer, drug resistance, tubulin, colchicine, colchicine derivatives

Citation: Pallante L, Rocca A, Klejborowska G, Huczynski A, Grasso G, Tuszynski JA and Deriu MA (2020) In silico Investigations of the Mode of Action of Novel Colchicine Derivatives Targeting β-Tubulin Isotypes: A Search for a Selective and Specific β-III Tubulin Ligand. Front. Chem. 8:108. doi: 10.3389/fchem.2020.00108

Received: 29 October 2019; Accepted: 04 February 2020;
Published: 21 February 2020.

Edited by:

Kamil Kuca, University of Hradec Králové, Czechia

Reviewed by:

John Holmes Miller, Victoria University of Wellington, New Zealand
Ambarish Kunwar, Indian Institute of Technology Bombay, India

Copyright © 2020 Pallante, Rocca, Klejborowska, Huczynski, Grasso, Tuszynski and Deriu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marco A. Deriu, marco.deriu@polito.it

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.