Next Article in Journal
Acridine Based N-Acylhydrazone Derivatives as Potential Anticancer Agents: Synthesis, Characterization and ctDNA/HSA Spectroscopic Binding Properties
Next Article in Special Issue
Design, Synthesis, Docking, DFT, MD Simulation Studies of a New Nicotinamide-Based Derivative: In Vitro Anticancer and VEGFR-2 Inhibitory Effects
Previous Article in Journal
Insight into Spodium–π Bonding Characteristics of the MX2⋯π (M = Zn, Cd and Hg; X = Cl, Br and I) Complexes—A Theoretical Study
Previous Article in Special Issue
Arene Ru(II) Complexes with Difluorinated Ligands Act as Potential Inducers of S-Phase Arrest via the Stabilization of c-myc G-Quadruplex DNA
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 9
10.3390/molecules27092873
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin

by
Sergei Boichuk
1,2,3,*,
Kirill Syuzov
1,
Firuza Bikinieva
1,
Aigul Galembikova
1,
Svetlana Zykova
4,
Ksenia Gankova
4,
Sergei Igidov
4 and
Nazim Igidov
4
1
Department of Pathology, Kazan State Medical University, 420012 Kazan, Russia
2
Department of Radiotherapy and Radiology, Russian Medical Academy of Continuous Professional Education, 125993 Moscow, Russia
3
Biologically Active Terpenoids Laboratory, Kazan Federal University, 18 Kremlyovskaya St., 420008 Kazan, Russia
4
Department of Pharmacology, Perm State Academy of Pharmacy, 614990 Perm, Russia
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(9), 2873; https://doi.org/10.3390/molecules27092873
Submission received: 3 April 2022 / Revised: 27 April 2022 / Accepted: 29 April 2022 / Published: 30 April 2022
(This article belongs to the Special Issue Design, Synthesis and Applications of New Anti-cancer Agents II)
Browse Figures
Review Reports Versions Notes

Abstract

:
Despite the tubulin-binding agents (TBAs) that are widely used in the clinic for cancer therapy, tumor resistance to TBAs (both inherited and acquired) significantly impairs their effectiveness, thereby decreasing overall survival (OS) and progression-free survival (PFS) rates, especially for the patients with metastatic, recurrent, and unresectable forms of the disease. Therefore, the development of novel effective drugs interfering with the microtubules’ dynamic state remains a big challenge in current oncology. We report here about the novel ethyl 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) exhibiting potent anti-cancer activities against the breast and lung cancer cell lines in vitro. This was due to their ability to inhibit tubulin polymerization and induce cell cycle arrest in M-phase. As an outcome, the EAPC-treated cancer cells exhibited a significant increase in apoptosis, which was evidenced by the expression of cleaved forms of PARP, caspase-3, and increased numbers of Annexin-V-positive cells. By using the in silico molecular modeling methods (e.g., induced-fit docking, binding metadynamics, and unbiased molecular dynamics), we found that EAPC-67 and -70 preferentially bind to the colchicine-binding site of tubulin. Lastly, we have shown that the EAPCs indicated above and colchicine utilizes a similar molecular mechanism to inhibit tubulin polymerization via targeting the T7 loop in the β-chain of tubulin, thereby preventing the conformational changes in the tubulin dimers required for their polymerization. Collectively, we identified the novel and potent TBAs that bind to the colchicine-binding site and disrupt the microtubule network. As a result of these events, the compounds induced a robust cell cycle arrest in M-phase and exhibited potent pro-apoptotic activities against the epithelial cancer cell lines in vitro.

1. Introduction

The microtubules are known as important regulators of many aspects of cellular functions, including cell proliferation and migration, vesicle trafficking during endocytosis, chromosomal segregation during mitosis, etc. [1,2,3,4]. The last one has made them an attractive target for the development of anti-cancer drugs. The protein tubulin, a major compartment of the microtubules, contains several various binding sites for the small-molecule drugs. This includes the laulimalide, maytansine, taxane/epothilone, vinca alkaloid, and colchicine sites [5,6,7,8,9,10]. Based on their effects on the microtubule dynamic state, the tubulin-binding agents (TBAs) are categorized into 2 major groups: microtubule-stabilizing agents including the taxanes, epothilones, and laulimalide and microtubule-destabilizing agents that are composed of colchicine, the vinca alkaloids, and maytansine [11,12].
Despite the impressive response rate shortly after the initiation of taxanes- and/or vinca alkaloids-based therapies, the extended use of these chemotherapeutic agents is limited due to the acquired resistance of tumors to these drugs. This might be due to the broad spectrum of molecular mechanisms: e.g., overexpression of ABC-transporters, such as P-glycoprotein in tumor cells [13,14], mutations in β-tubulin, and the altered expression of specific β-tubulin isotypes [15,16], etc. In this view, the agents targeting the colchicine-binding site (CBS) exhibit several advantages, including overcoming multidrug resistance (MDR), which is supplemented by their well-known abilities to inhibit angiogenesis [17,18,19,20]. Moreover, CBS inhibitors (CBSi) have several significant advantages over other microtubule-targeting agents (MTAs) because this site of tubulin is more amenable to the molecules with favorable physicochemical properties that improve oral bioavailability over taxanes and vinca alkaloids binding sites, for example, have less drug-drug interactions, and are less prone to developing multi-drug resistance. Despite the clinical use of colchicine as an anti-cancer derivative, it is currently excluded based on its extreme toxicity; the broad number of CBSi is currently described [21] and can be used as structurally diverse scaffolds for the generation of potent anti-cancer derivatives exhibiting limited toxicities, improved solubility, etc. [22].
In general, the scaffolds to develop the novel TBAs targeting the CBS are based on the presence of diverse chemical groups, including the trimethoxyphenyl (TMP) group. This principle is raised from the structure of colchicine, which consists of a 3,4,5-trimethoxyphenyl ring (the A ring), a saturated seven-membered ring containing an acetamido group at position 7 (the B ring), and a tropolone ring (the C ring). It was shown that the TMP group of colchicine is oriented within β-tubulin close to Cys β241 and thereby allows colchicine to bind to β-tubulin at its interface with α-tubulin, subsequently inhibiting tubulin polymerization [23]. Thus, it was accepted that the most common conservative chemical groups in the compounds exhibiting tubulin depolymerizing properties are composed of trimethoxyphenyl (TMP) moiety (ring A) and the linker region connecting ring A and ring B [24,25,26]. This point of view was mainly based on the structure-activity relationship (SAR) analysis of colchicine and combretastatin A-4 and its analogs, the well-known earliest inhibitors of tubulin polymerization.
The present study aimed to examine the cytotoxic activities of the 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) synthesized in our lab and further assess their molecular mechanisms of anti-cancer activities as the novel chemical structures exhibiting tubulin-binding activities.
Here, we show for the first time that 2-aminopyrrole derivatives (EAPC-67 and -70) exhibit potent cytotoxic activities against a broad spectrum of cancer cell lines (e.g., triple-negative breast cancer HCC 1806, MDA-MB-231; non-small cell lung cancer H 1299, etc.). This was evidenced by the decreased viability of EAPC-treated cancer cells as measured by MTS-based assay, increased expression of apoptotic markers (e.g., cleaved forms of PARP, caspase-3), and the number of apoptotic (e.g., Annexin V-positive) cells. These activities of EAPCs were due to their abilities to halt a cell cycle progression and induce a robust cell cycle arrest in M-phase. This, in turn, was a consequence of the inhibition of the microtubule polymerization due to the ability of EAPCs to bind with the CBS of the tubulin. To reveal this possibility, we used different in silico molecular modeling methods. First, we used a blind docking strategy to determine the preferred drug-binding site on the tubulin dimer. We found the CBS of tubulin to be the most energetically favorable for EAPCs indicated above. Next, we used the induced docking of the ligands in the CBS. By utilizing this approach, we identified two of the most probable poses of EAPC-67 exhibiting the maximal induced-fit docking (IFD) scores. Afterward, by using the binding metadynamics approach, we determined the stability of the poses for these selected ligands by calculating the PoseScore and CompScore indexes. Finally, based on the metadynamics data, we performed 100 nanosecond unbiased molecular dynamics to determine the key interactions of the most stable ligand pose with the CBS. Data analysis of the induced fit docking and the trajectory of molecular dynamics revealed that the external orientation of the T7 loop of the β-chain of tubulin is required for the successful binding of the drug to the CBS. Otherwise, steric clashes are created between the ligand atoms and the amino acid atoms of the T7 loop that block access for the ligand to its binding site. This data indicates that the molecular mechanism of action of EAPCs might be similar to colchicine and is predominantly associated with the impossibility of internal turning of the T7 loop in the β-chain of tubulin. Based on our data, we also proposed that EAPCs block the rotation of the T7 loop and subsequently prevent the conformational changes in the tubulin dimer from the curved to the straight state which is required for tubulin polymerization.

2. Materials and Methods

2.1. Chemical Compounds

Paclitaxel (PTX), vinblastine (Vin), and colchicine (Col) were purchased from Sigma, St Louis, Missouri, USA, and dissolved in 100% dimethyl sulfoxide (DMSO) (Sigma Aldrich, St. Louis, MO, USA). 10 novel 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) as potential microtubule destabilizing agents were synthesized in our laboratory according to the standard protocols, as shown in Figure 1 and in previous studies [27].

2.2. Molecular Modeling Methods

Blind docking: To identify the potential binding sites of EAPCs on the tubulin, the molecular docking procedure was performed by using Schrodinger molecular modeling software (Schrödinger Inc., New York, NY, USA, 2021). The Glide Docking XP mode was used to find the tubulin-ligand supramolecular complex with the minimal energy scoring function. The tubulin structure (no. 4O2B) was downloaded from the Protein Data Bank (PDB) (4O2B. Available online: www.rcsb.org/structure/4O2B, accessed on 3 April 2022). The protein structure was prepared before docking using the Protein Preparation Wizard module. Ligand preparation was performed by the standard protocol using the module LigPrep [28]. To determine the possible interactions of the most active compounds with the tubulin molecule, “blind docking” was performed. Only tubulin inhibitor binding sites were used for analysis. After defining binding sites, a Grid Box for each of them was generated by using the receptor grid module.
Induced Fit Docking: The Standard Induced Fit Docking Protocol is used to generate possible induced CBS conformations [29]. Initially, the prepared ligands were docked using a softened potential with Glide XP [30]. The Coulomb-vdW scaling factor was set to 0.5 for the ligand and protein atoms, and a maximum of 20 poses were generated on one ligand tautomer. The generated poses were further processed by Prime for side-chain refinements within 7 Å of the binding site for better accommodation of the ligands. Next, the protein-ligand system was minimized using the OPLS4 force field [31]. In the final stage, the ligands were docked using the Glide XP module into the optimized binding site generated within 30 kcal mol−1 of the lowest energy structure obtained after Prime refinement to generate 20 poses per system. To rank the poses of the ligands, we used IFD score (IFDScore = 1.0 × GlideScore + 0.05 × PrimeEnergy). Based on the maximum IFD score the conformations of the ligand-protein complex were chosen for the next steps of molecular modeling.
Metadynamics: The metadynamics method was used to assess the stability of the EAPC-67 pose in the selected induced-fit conformations. 10 independent metadynamics simulations of 10 ns are performed using root-mean-square deviation (RMSD) of the ligand heavy atoms as the collective variable. The alignment before the RMSD calculation was done by selecting protein residues within 3 Å of the ligand. The hill height and width were set to 0.05 kcal/mol (about 1/10 of the characteristic thermal energy of the system, kBT) and 0.02 Å, respectively [32]. To remove bad contacts and/or strain in the initial starting structure and to slowly reach a temperature of 300 K for the simulated system, a 0.5 ns simulation was performed in the NVT ensemble before the start of the metadynamics itself. The EAPC-67-tubulin complex was solvated in a simple point charge (SPC) water model using an orthorhombic box with periodic boundary conditions. The total charge of the system was neutralized by the addition of sodium and potassium ions. The final snapshot of the short unbiased molecular dynamics (MD) simulation of 0.5 ns is then used as the reference structure for following the metadynamics production.
Molecular dynamics: After the metadynamics procedure, 100 ns unbiased molecular dynamics was performed to assess the position of the ligand in the CBS and determine key amino acid interactions over time. The Desmond program in Schrödinger suite 2021-2 was used to perform classical molecular dynamics of the EAPC-67-tubulin complex [33]. The EAPC-67-tubulin complex was solvated in a simple point charge (SPC) water model using an orthorhombic box with periodic boundary conditions. The total charge of the system was neutralized by the addition of sodium and potassium ions. To minimize and relax the system, a short simulation in the NVT ensemble of 100 picoseconds was performed. The main simulation was performed in the NPT ensemble. The simulation was run for a total of 100 ns with a recording interval of 100 ps. The temperature and pressure were kept constant at 300 K and 1.01325 bar, respectively, throughout the simulations. Calculation of the root means square deviations (RMSD) and analysis of the protein-ligand interaction were performed using the simulation interaction diagram (SID) program tool in the Schrödinger Maestro Suite.
Trajectory analysis and clusterization: To investigate conformational flexibility over MD simulation, the collective motions of the EAPC-67-tubulin complex were investigated for all snapshots in the trajectories. For the extraction of snapshots from 100 nanosecond molecular dynamics, VMD software was used [34]. After snapshot extraction using the pca.xyz() function in the Bio3D package in R, we performed a principal component analysis (PCA) analysis over cartesian coordinates of Cα atoms of all snapshots in the trajectories [35,36]. Each trajectory was projected onto the first eigenvector subspaces to reveal concerted atomic displacements, representing the “essential dynamics” of the protein. Using the PCA data, K-means clusterization analysis was performed over each trajectory described by the first three principal components. The optimal number of clusters was calculated using the Silhouette method [37]. After clustering, a medoid was selected from each cluster as a representative conformation.
The molecular mechanics/generalized Born surface area (MM-GBSA): The MM-GBSA method was used as a post-molecular dynamic validation tool [38]. Based on the medoids of clusters, we calculated the binding free energy (∆Gbind) of each EAPC-67-tubulin complex, including the solvation energy (∆Gsolv), using the implicit water model. A representative conformation of the ligand-receptor complex was chosen based on minimal ΔGbind Prime energy.

2.3. Cell Lines and Culture Conditions

The human basal-like TNBC cell lines HCC 1806, MDA-MB-231, and non-small cell lung cancer cell line H1299 were purchased from the American Type Culture Collection (Manassas, VA, USA). The cell lines indicated above were maintained in RPMI-1640 medium (Paneco, Moscow, Russia) supplemented with 15% fetal bovine serum (Gibco; Thermo Fisher Scientific Inc., Waltham, MA, USA), 50 U/mL penicillin, and 50 μg/mL streptomycin. The cell lines were cultured at 37 °C in a humidified atmosphere of 5% CO2 in an incubator (LamSystems, Miass, Russia).

2.4. Antibodies

The primary antibodies used for western blotting were anti-PARP (cat. no. 436400; Invitrogen; Thermo Fisher Scientific Inc., Waltham, MA, USA), cleaved caspase-3 (cat. no. 9661S), phospho-NuMA Ser395 (cat. no. 3429), phospho-Histone H3 Ser10 (cat. no. 53348) (Cell Signaling Technology Inc., Danvers, MA, USA), and beta-actin (cat. No. A00730-200, GenScript, Piscataway, NJ, USA); HRP-conjugated secondary antibodies, anti-mouse immunoglobulin (Ig)G (cat. no. sc-2005) and anti-rabbit IgG (cat. no. sc-2004), were purchased from Santa Cruz Biotechnology, Dallas, TX, USA.

2.5. Western Blotting Analysis

To examine the protein expression in parental and Tx-R cells, whole-cell lysates (WCL) were prepared by scraping the cells growing as monolayers into a radio-immunoprecipitation buffer (RIPA buffer) (25 mM Tris-HCl pH 7.6, 5 mM EDTA, 150 mM NaCl, 0.1% SDS, 1% NP-40, 1% sodium deoxycholate) supplemented with the cocktail of protease and phosphatase inhibitors. The cellular lysates were further incubated for 1 h at 4 °C and clarified by centrifugation for 30 min at 11,400 rpm at 2 °C. The protein concentrations in WCL were calculated by the Bradford assay. The protein samples (20 µg) were loaded on the 4–12% Bis-Tris or 3–8% Tris-acetate NuPAGE gels (Invitrogen, Carlsbad, CA, USA) and upon completion of electrophoresis transferred to a nitrocellulose membrane (Bio-Rad, Hercules, CA, USA), membranes were probed with primary (1:1000 and incubated overnight at 4 °C), and secondary antibodies (1:1000 and incubated for 1 h at room temperature) and visualized by enhanced chemiluminescence (Western Lightning Plus-ECL reagent, Perkin Elmer, Waltham, MA, USA).

2.6. Tubulin Polymerization Assay

The impact of the CAs on the dynamic state of the tubulin polymerization was assessed by using the Tubulin Polymerization kit (Cytoskeleton Inc., Denver, CO, USA) as specified by the manufacturer. Results were obtained on a SpectraFluor Plus microplate reader (Tecan GmbH, Grödig, Salzburg, Austria) and readings were taken every minute for 1 h (61 measurements in total).

2.7. Flow Cytometry

The number of apoptotic cells was counted by using Muse Annexin V Dead Cell Kit (Merck KGaA, Darmstadt, Germany) according to the manufacturer’s instructions. Unstained and single-stained untreated cells were used as controls. Cells were analyzed in a Muse Cell Analyzer (Merck KGaA, Darmstadt, Germany). For all the experiments indicated above at least 10,000 events were acquired for each sample. Results were presented as the percentage of desired cells relative to the total number of cells as mean ± standard deviation of four biological repeats.

2.8. Statistics

All the experiments have been repeated a minimum of 3 times. The results are presented as the mean ± standard error (SE) for each group. Statistical analyses (Student’s t-test) were performed using Statistical software program version 7.0 (S.A. Glantz, McGraw Hill Education, New York, NY, USA). p < 0.05 was considered to indicate a statistically significant difference.

3. Results

3.1. Synthesis of 2-Aminopyrrole Derivatives and Their Characteristics

In the current study we performed the synthesis of ethyl 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (shown in Figure 1 as 2a-j). These compounds were obtained by recycling furanones 1 under the action of cyanoacetic ester in the presence of a catalyst-triethylamine.
The chemical compounds indicated above were synthesized according to the standard protocols. Briefly, a solution of 0.002 mol of the corresponding 3-(4-methylbenzoyl) hydrazone 2,3-dihydro-2,3-furandione 1 in 20 mL of anhydrous dioxane was added to an equivalent amount of cyanoacetic acid ethyl ester and triethylamine. The resulting mixture was heated for 30–40 min and further cooled to 0 °C. The precipitate was filtered and recrystallized from ethanol.
The resulting compounds 2a-j are yellow-colored crystalline substances, soluble in dimethyl sulfoxide, dimethylformamide, poorly soluble in ethanol, and insoluble in benzene, water, and hexane. The infrared (IR) spectra of synthesized compounds 2a-j are characterized by the presence of absorption bands of stretching vibrations of NH groups in the region of 3109–3414 cm−1, carbonyl groups in the region of 1692–1712 and 1646–1676 cm−1, as well as stretching vibrations of the C = C double bond in the region of 1590–1628 cm−1. In the 1H nuclear magnetic resonance (NMR) spectra, in addition to the signals of substituents in the aromatic ring and the ester group, singlets of protons of the NH group of the amide fragment at 10.47–10.90 ppm, signals of protons of the amino group at 8.85–9.20 ppm and 8.27–8.48 ppm are recorded, respectively; the vinyl proton singlet was recorded at 6.32–6.64 ppm.
Infrared spectra were taken on an IR-Fourier spectrophotometer-spectrometer, model FSM-1201 (SPb Instruments, Saint-Petersburg, Russia) in disks with potassium bromide; 1H NMR spectra were taken on Bruker Avance III apparatus (Bruker Avance, Bremen, Germany) in DMSO-d6, operating frequencies 400 MHz (1H) and 100 MHz (13C), the internal standard is the residual signal from the deuterium solvent. Elemental analysis was performed on a Leco CHNS-932 (Leco Corporation, St. Joseph, MI, USA). The chemical purity of the compounds and the progress of the reactions were monitored by TLC on Sorbfil PTSKh P-A-UF-254 plates (Sorbpolymer, Krasnodar, Russia) in the ether-benzene-acetone solution (ratio 10:9:1); detection was carried out with iodine vapor. Melting points were determined on an SMP40 apparatus (Bibby Scientific Ltd., Stone, UK). 2D structure formulas and spectral characteristics of the synthesized compounds are shown in Supplementary Figures S1 and S2, respectively.

3.2. EAPCs Reduce the Viability of the Epithelial Cancer Cell Lines In Vitro

To examine whether EAPCs are active against epithelial cancer cell lines, we performed the MTS-based survival assay by using HCC1806 MDA-MB-231 breast cancer and H1299 non-small cell lung cancer (NSCLC) cancer cells. Each cancer cell line indicated above was treated with various concentrations of EAPCs (0.01–100 μM) for 48–72 h. PTX, Vin and Col were also included as the positive controls and exhibited cytotoxic activities in nanomolar concentrations for all cancer cell lines indicated above. The IC50 values for EAPCs are shown in Table 1, illustrating that EAPC-67, -70, and -71 exhibited the most potent cytotoxic activities for all cancer cell lines included in the present study.

3.3. EAPC-67 and -70 Arrest Cancer Cells in M-Phase and Inhibit Tubulin Polymerization

Given that EAPC-67, -70, and -71-treated cancer cells acquired a round-shape morphology (data is not shown), we proposed that the growth-inhibitory activities of these compounds were due to the ability to induce mitotic arrest and interfere with the microtubule dynamic state. To examine this possibility, we performed the western blotting analysis to assess the expression of the proteins specifically accumulated in the M-phase of the cell cycle. Indeed, the expression of pH3 Ser10 and pNuMA Ser395 (the well-known mitotic markers) was significantly increased in EAPC-treated cells, as shown in Figure 2. Importantly, in H1299 NSCLC cells of all three types of EAPCs exhibited much more potent activity to induce cell cycle arrest in M-phase when compared to the well-known TBAs, vinblastine, and colchicine.
To further reveal whether the accumulation of EAPC-treated cancer cells was due to their ability to interfere with the microtubule dynamic state, we performed the tubulin polymerization assay. PTX- and Vin-treated tubulin samples were used as positive and negative controls, respectively. We found that EAPC-67 and -70 effectively inhibited tubulin polymerization, as shown in Figure 3. As expected, PTX exhibited high polymerization activity with tubulin, whereas Vin effectively inhibited tubulin polymerization.

3.4. EAPCs Induce Apoptosis of Breast and Lung Cancer Cells

Next, we examined the pro-apoptotic activities of EAPC-67 and -70- in epithelial cancer cell lines. For this purpose, we initially assessed the expression of apoptotic markers (cleaved forms of caspase-3 and PARP) by western blotting. Indeed, we found a significant increase in the expression of apoptotic markers in breast and lung cell lines after EAPC treatment (Figure 4).
Similar to WB data, FACs analysis revealed a significant increase of apoptotic (i.e., Annexin V-positive) cells after EAPC treatment (Figure 5).
Given that EAPC-67 was found to be most active against both of the breast cancer cell lines and also exhibited potent anti-proliferative and cytotoxic activities against the NSCLC cell line (Table 1 and Figure 4 and Figure 5), we further performed computational-based analysis to determine the molecular mechanism of action of this compound.

3.5. Molecular Modeling Studies

The molecular modeling study was composed of four parts: (1) Blind docking to determine the most possible binding site in the tubulin for EAPC-67; (2) Binding metadynamics to select the most stable conformation for EAPC-67; (3) Unbiased molecular dynamics of the most stable pose of EAPC-67 within 100 nanoseconds; and (4) Analysis of the structural changes in tubulin upon EAPC-67 binding.
First, we generated EAPC-67 tautomers by using the LigPrep module and further used them for the blind docking to assess their binding with the well-known tubulin binding sites, including colchicine, and the vinca alkaloid, maytansine, and the pironetin binding site. The binding site with the lowest GlideXP score was considered the most possible for ligand binding. Based on the GlideXP score values, we found the CBS as the most energetically favorable for EAPC-67. It is well-known that the protein molecules are not rigid and might undergo conformational changes, in particular after binding with the ligand, thereby suggesting the obvious limitations associated with protein molecule rigidity. Therefore, to more accurately predict the ligand pose, we performed induced-fit docking of two EAPC-67 tautomers with the previously selected CBS.
For each of the EAPC-67 tautomers, 20 ligand poses were generated by the induced docking algorithm. According to the results of induced fit docking, the pose with the maximum IDF score was chosen for each of the two tautomers. To select an EAPC-67 tautomer that is more likely to bind to the tubulin molecule, we performed binding pose metadynamics (BPMD), which allows us to assess the ligand’s stability in solution. Ligand poses that are unstable under the bias of the metadynamics simulation are expected to be infrequently occupied in the energy landscape, thus making minimal contributions to the binding affinity. The results of BPMD were evaluated by the CompScore value, where CompScore is the linear combination of PoseScore and PersScore and is equal to PoseScore − 5 × PersScore. For each of the tautomers, the CompScore was 0.158 and 1.436, respectively (Figure 6). Lower values of the CompScore indicate that the protein-ligand complexes exhibit the highest stability.
Next, we subjected the EAPC-67 tautomer conformation with maximum CompScore and PoseScore to 100 nanosecond unbiased molecular dynamics under periodic conditions. Molecular dynamics was used for the refinement of an IFD Pose and to determine the key interactions between ligand atoms and the amino acids of the CBS. Also, MD allowed us to evaluate the ensemble of available conformations for the complex and generate the “best” single structure approximation (representative conformation) for the complex EAPC-67-tubulin by performing the clusterization assay. An analysis of the structure of the representative conformation provided a better understanding of the mechanism of inhibition of tubulin polymerization induced by the synthesized compounds. The trajectory analysis showed the low RMSD values of the ligands and proteins (Figure 7), which indicates the absence of large fluctuations during the entire time of molecular dynamics.
Over the 100 nanoseconds of molecular dynamics, EAPC-67 predominantly formed the non-covalent contacts with the following amino acids B: Leu 248, B: Asn 249, B: Ala 250, B: Asp 251 of the T7 loop, with B: Lys 254, B: Leu 255 of the H8 helix and with B: Lys 352 of the S9 sheet. Therefore, we found that the ligand interacted predominantly with the amino acids of the secondary elements of the T7 loop, H8 helix, and S9 sheet of CBS. The types of non-covalent interactions with these amino acids, including the formation of hydrogen, hydrophobic, ionic interactions, and the formation of water bridges, as well as the contribution of each interaction type to the overall interaction with a particular amino acid, are shown in the Figure 8A and Figure 9. The time-dependent total number of contacts between the ligand and protein and the number of contacts between the ligand and specific amino acids are shown in Figure 8B.
To select a representative conformation of the ligand-receptor complex, we performed clustering of the molecular dynamics snapshot by the K-means method based on PCA analysis. First, a PCA analysis of the EAPC-67-tubulin complex trajectory was performed, as shown in Figure 10A. This plot demonstrates the distribution of the EAPC-67-tubulin complex conformations in the PC 1–2 subspace, where the blue-white-red color scale reflects the time of the simulation. The contribution of each principal component to the overall variability of atomic positions during molecular dynamics is shown in Figure 10B. The cumulative contribution of the first three PCs cover 15.4%, 26.5%, and 34.5% of the total variance, respectively. Since the number of clusters is critical for the K-means method, we used the Silhouette method to determine their optimal numbers. Based on the Silhouette graph, the optimal number of clusters was determined as 3. Then, clustering of the EAPC-67-tubulin complex conformations by the K-means method was carried out in the PC 1–3 subspace, followed by the projection of the cluster members onto the 2D PC1-2 graph (Figure 10C). The cluster members were further projected onto the RMSD plot to visualize the relationship between conformations that belong to each particular cluster and time. Eventually, we found three large clusters of EAPC-67-tubulin complex conformations that were clearly separated in time with medoids (Snapshots 834, 1556, and 1987) corresponding to each cluster (Figure 10D).
To determine the representative conformation of three medoids, we utilized the MM-GBSA method to evaluate the protein-ligand complex exhibiting the minimal value of ΔGbind prime energy.
As shown in Table 2, the conformation of the EAPC-67-tubulin complex represented as a medoid of cluster 2 exhibited the lowest ΔGbind binding energy and ligand strain energy. The lowest value of ΔGbind indicates that the medoid of cluster 2 is the most probable conformation when compared to the other medoids. The lowest ligand strain energy values also illustrate that minimal energy is required for the adaptation of EAPC-67 to the receptor-bound conformation.
After selecting a representative conformation, we aligned the CBS site to the straight conformation of tubulin obtained from the PDB database (PDB ID: 1JFF). Structural differences were most prominent in the β-chain of tubulin, showing that upon ligand binding the rotation of the T7 loop is necessary. Otherwise, critical steric clashes are formed between the atoms of the ligand and the amino acids of the T7 loop, thereby making the interaction of the protein and the ligand impossible. We also found the displacements of other secondary elements, (e.g., translation of the H7 helix, upward movement, and rotation of S8 and S9 β-sheet) required for the enlargement of free space (Figure 11). Such structural changes increase the available space for the ligand but completely block the transition of the tubulin dimer from the curve to the straight conformation. As a consequence, an effect of the EAPC-67 was due to its ability to interfere with the concerted movements of these secondary structure elements required for tubulin to adopt its straight, microtubular conformation and, therefore, to assemble in the microtubules.

4. Discussion

Pyrrole-based compounds are known as one of the attractive scaffolds in medicinal chemistry [39]. Indeed, the privileged structures represent the molecular scaffolds with versatile binding properties, such that a single scaffold can provide potent and selective ligands for a range of different biological targets through modification of the functional groups [40,41]. Therefore, the use of pyrrole–based compounds is considered an attractive strategy to obtain potential drugs due to its ability to modulate several pharmacokinetics parameters, such as solubility, lipophilicity, polarity, and pharmacodynamics, including hydrogen-bonding capacity and the ability to form the complexes with coordinating metals [42,43]. This explains the existence of a broad spectrum of pyrrole-based drugs, including anti-fungal agents, anti-microbial agents, anti-inflammatory agents, HMG-CoA reductase inhibitors, antidepressants, antihypertensive agents, antimalarial agents, anticancer agents, anti-HIV-1 agents, etc. [44,45,46]. To our knowledge, pyrrole and pyrrole-fused heterocycles are suitable to develop novel and effective scaffolds for drugs exhibiting anti-cancer activities due to the targeting of the CBS. Indeed, these scaffolds exhibit high capacities for expanding the chemical space of tubulin inhibitors due to the existence of the various ways to synthesize the derivatives with various functional groups, as well as a wide modification of existing candidate molecules [47,48,49,50,51,52,53].
In addition to the well-established and informative in vitro assays to examine the biological activities of the newly synthesized compounds, computer-aided drug discovery (CADD) has emerged as a powerful and promising technology for cheaper, faster, and effective drug design. These CAAD approaches based on various molecular modeling methods, including standard docking, induced-fit docking, biased/unbiased molecular dynamics, MM-GBSA, etc. are widely used in the field of anti-cancer research [54,55,56]. All of them were shown as effective tools to develop novel inhibitors targeting the various binding sites of proteins, including the CBS of the tubulin. This is due to the availability of the multiple X-ray structures of tubulin-ligand complexes, which provide important information about the conformational changes in the tubulin molecule and interactions between various functional groups of ligands and amino acids of the CBS [57,58,59,60]. The use of the in silico approach significantly reduces the cost to develop novel drugs and allows the negative selection of non-effective molecules at the early stage of drug discovery. Moreover, SAR analysis of various functional chemical groups allows us to find out the novel patterns of interactions between ligand atoms and specific protein amino acids. In addition, such interaction analysis is effective for the prediction of the tubulin depolymerization rate, cytotoxicity against cancer cells, etc. For example, Lorenzo Pallante et al. modeled the β-III tubulin isotype protein structure based on the homology modeling method to discover novel compounds exhibiting high specificity to β-III tubulin isotype [61]. Laura Gallego-Yerga et al. performed a virtual screening of the Zinc database to find potential inhibitors of the colchicine site by using an ensemble docking approach, molecular dynamics, and pharmacophore modeling. As a result, the authors synthesized tetrazole derivatives and further tested their potential activity to interfere with tubulin polymerization in vitro. In good agreement with the design principles, this compound demonstrated potent in vitro activity against tubulin polymerization in a micromolar range and an anti-proliferative effect against human epithelioid carcinoma HeLa cell line in nanomolar concentrations [62]. Based on the 2D-quantitative SAR analysis of quinolines and using the docking method, Mohamed et al. synthesized 29 quinoline derivatives targeting the CBS and exhibiting potent anti-proliferative activities against various human cancer cell lines including colorectal carcinoma (HCT-116), hepatocellular carcinoma (HepG-2), and breast cancer (MCF-7) [63]. In addition, abundant examples illustrating the effectiveness of computer simulation modeling approaches for the design and synthesis of tubulin inhibitors targeting the CBS were shown and reviewed recently [64,65,66,67].
In agreement with current studies, we used a combination of computer modeling methods indicated above and biological activity assays to discover the anti-cancer activities of pyrrole-based compounds targeting the CBS of tubulin. As a result of our CAAD approaches, we report here about the novel pyrrole-based active compounds exhibiting potent cytotoxic activities against distinct cancer cell lines, including HCC1806, MDA-MB-231 triple-negative breast cancer, and H1299 non-small cell lung cancer cells. Among all synthesized compounds, EAPC-67 and EAPC-70 were the most active, demonstrating the highest cytotoxic activities against cancer cells in a dose-dependent manner. Next, by using the polymerization assay, we found that the cytotoxic activities of EAPCs are due to their ability to interfere with the microtubule dynamic state (Figure 3). Molecular modeling data revealed that the EAPCs interact directly with the CBS and induce the conformational changes of the tubulin dimer that are typical to the mode of action of the compounds targeting the CBS [68,69]. In particular, translation of the H7 helix, upward movement, rotation of the S9 β-sheet, and rotation of the T7 loop in the β-chain of tubulin was demonstrated. The transition of these secondary elements is impossible to the position characteristic of the straight tubulin conformation in the assembled microtubule because the EAPCs occupy the physical space between the α- and β-chains. Therefore, EAPCs-bound free tubulin molecules cannot assemble into organized microtubules (Figure 11). Our molecular modeling data was a proper fit with the biological activities observed for EAPCs. Indeed, we found that EAPC-67 and -70 effectively inhibited tubulin polymerization. As an outcome of EAPC-induced failure of assembly of microtubules, we observed a robust cell cycle arrest in the G2/M phase and the accumulation of cells in the M-phase. This was evidenced by the substantial increase of pH3 Ser10 and pNuMA Ser395 expression in cancer cells treated with EAPCs (Figure 2). As a result of the mitotic arrest, EAPCs-treated cells underwent apoptotic cell death, which was evidenced by a substantial increase of the cleaved forms of caspase-3 and PARP in the vast majority of cancer cell lines and also a significant increase of Annexin V-positive cells after EAPC treatment (Figure 4 and Figure 5, respectively). Collectively, new synthesized EAPCs compounds, in particular EAPC-67 and EAPC-70, can be used as a potential scaffold for the development of new and more promising MTAs targeting the colchicine-binding site.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27092873/s1, Figure S1: 2D structural formulas of synthesized EAPCs, Figure S2: Spectral characteristics of synthesized EAPCs.

Author Contributions

Conception and design, S.B.; development of methodology, S.B., A.G., F.B., K.S., K.G., S.Z., N.I. and S.I.; acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.), S.B., A.G., K.S., F.B., S.Z., N.I., K.G. and S.I.; analysis and interpretation of data (statistical analysis, biostatistics, computational; analysis), S.B., A.G., K.S., F.B., K.G. and S.Z.; writing—review and/or revision of the manuscript, S.B.; administrative—technical, or material support (reporting or organizing data, constructing databases), S.B.; study supervision, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by a grant from the Russian Science Foundation (RSF) (grant no. 20-015-0001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

The research was performed solely on the mammalian cell lines (informed consent and animals statement are not required).

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by Kazan Federal University Strategic Academic Leadership Program (PRIORITY-2030).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Avila, J. Microtubule functions. Life Sci. 1992, 50, 327–334. [Google Scholar] [CrossRef]
  2. Vukušić, K.; Buđa, R.; Tolić, I.M. Force-generating mechanisms of anaphase in human cells. J. Cell Sci. 2019, 132, jcs231985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. de Forges, H.; Bouissou, A.; Perez, F. Interplay between microtubule dynamics and intracellular organization. Int. J. Biochem. Cell Biol. 2012, 44, 266–274. [Google Scholar] [CrossRef] [PubMed]
  4. Bonifacino, J.S.; Neefjes, J. Moving and positioning the endolysosomal system. Curr. Opin. Cell Biol. 2017, 47, 266–274. [Google Scholar] [CrossRef] [PubMed]
  5. Mooberry, S.L.; Tien, G.; Hernandez, A.H.; Plubrukarn, A.; Davidson, B.S. Laulimalide and isolaulimalide, new paclitaxel-like microtubule-stabilizing agents. Cancer Res. 1999, 59, 653–660. [Google Scholar]
  6. Hamel, E. Natural products which interact with tubulin in the vinca domain: Maytansine, rhizoxin, phomopsin a, dolastatins 10 and 15 and halichondrin B. Pharmacol. Ther. 1992, 55, 31–51. [Google Scholar] [CrossRef]
  7. von Hoff, D.D. The taxoids: Same roots, different drugs. Semin. Oncol. 1997, 24 (Suppl. 13), S13-3–S13-10. [Google Scholar]
  8. Bollag, D.M.; McQueney, P.A.; Zhu, J.; Hensens, O.; Koupal, L.; Liesch, J.; Goetz, M.; Lazarides, E.; Woods, C.M. Epothilones, a new class of microtubule-stabilizing agents with a taxol-like mechanism of action. Cancer Res. 1995, 55, 2325–2333. [Google Scholar]
  9. Gigant, B.; Wang, C.; Ravelli, R.B.G.; Roussi, F.; Steinmetz, M.O.; Curmi, P.A.; Sobel, A.; Knossow, M. Structural basis for the regulation of tubulin by vinblastine. Nature 2005, 435, 519–522. [Google Scholar] [CrossRef]
  10. Hastie, S.B. Interactions of colchicine with tubulin. Pharmacol. Ther. 1991, 51, 377–401. [Google Scholar] [CrossRef]
  11. Jordan, M.A.; Wilson, L. Microtubules as a target for anticancer drugs. Nat. Cancer 2004, 4, 253–265. [Google Scholar] [CrossRef] [PubMed]
  12. Stanton, R.A.; Gernert, K.M.; Nettles, J.H.; Aneja, R. Drugs that target dynamic microtubules: A new molecular perspective. Med. Res. Rev. 2011, 31, 443–481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Krishna, R.; Mayer, L.D. Multidrug resistance (MDR) in cancer: Mechanisms, reversal using modulators of MDR and the role of MDR modulators in influencing the pharmacokinetics of anticancer drugs. Eur. J. Pharm. Sci. 2000, 11, 265–283. [Google Scholar] [CrossRef]
  14. Mechetner, E.; Kyshtoobayeva, A.; Zonis, S.; Kim, H.; Stroup, R.; Garcia, R.; Parker, R.J.; Fruehauf, J.P. Levels of multidrug resistance (MDR1) P-glycoprotein expression by human breast cancer correlate with in vitro resistance to taxol and doxorubicin. Clin. Cancer Res. 1998, 4, 389–398. [Google Scholar]
  15. Kavallaris, M.; Kuo, D.Y.; Burkhart, C.A.; Regl, D.L.; Norris, M.; Haber, M.; Horwitz, S.B. Taxol-resistant epithelial ovarian tumors are associated with altered expression of specific beta-tubulin isotypes. J. Clin. Investig. 1997, 100, 1282–1293. [Google Scholar] [CrossRef] [Green Version]
  16. Stengel, C.; Newman, S.P.; Leese, M.P.; Potter, B.V.L.; Reed, M.J.; Purohit, A. Class III β-tubulin expression and in vitro resistance to microtubule targeting agents. Br. J. Cancer 2009, 102, 316–324. [Google Scholar] [CrossRef]
  17. Chaplin, D.J.; Hill, S.A. The development of combretastatin A4 phosphate as a vascular targeting agent. Int. J. Radiat. Oncol. 2002, 54, 1491–1496. [Google Scholar] [CrossRef]
  18. Siemann, D.W.; Shi, W. Dual targeting of tumor vasculature: Combining Avastin and vascular disrupting agents (CA4P or OXi4503). Anticancer Res. 2008, 28, 2027–2031. [Google Scholar]
  19. Lindamulage, I.K.; Vu, H.-Y.; Karthikeyan, C.; Knockleby, J.; Lee, Y.-F.; Trivedi, P.; Lee, H. Novel quinolone chalcones targeting colchicine-binding pocket kill multidrug-resistant cancer cells by inhibiting tubulin activity and MRP1 function. Sci. Rep. 2017, 7, 10298. [Google Scholar] [CrossRef] [Green Version]
  20. Gupta, S.; Banerjee, M.; Poddar, A.; Banerjee, A.; Basu, G.; Roy, D.; Bhattacharyya, B. Biphasic Kinetics of the Colchicine−Tubulin Interaction: Role of Amino Acids Surrounding the A ring of Bound Colchicine Molecule. Biochemistry 2005, 44, 10181–10188. [Google Scholar] [CrossRef] [Green Version]
  21. McLoughlin, E.C.; O’Boyle, N.M. Colchicine-Binding Site Inhibitors from Chemistry to Clinic: A Review. Pharmaceuticals 2020, 13, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Arnst, K.E.; Banerjee, S.; Chen, H.; Deng, S.; Hwang, D.; Li, W.; Miller, D.D. Current advances of tubulin inhibitors as dual acting small molecules for cancer therapy. Med. Res. Rev. 2019, 39, 1398–1426. [Google Scholar] [CrossRef] [PubMed]
  23. Ravelli, R.B.G.; Gigant, B.; Curmi, P.A.; Jourdain, I.; Lachkar, S.; Sobel, A.; Knossow, M. Insight into tubulin regulation from a complex with colchicine and a stathmin-like domain. Nature 2004, 428, 198–202. [Google Scholar] [CrossRef] [PubMed]
  24. Quinn, F.R.; Neiman, Z.; Beisler, J.A. Toxicity and quantitative structure-activity relationships of colchicines. J. Med. Chem. 1981, 24, 636–639. [Google Scholar] [CrossRef] [PubMed]
  25. Ringel, I.; Jaffe, D.; Alerhand, S.; Boye, O.; Muzaffar, A.; Brossi, A. Fluorinated colchicinoids: Antitubulin and cytotoxic properties. J. Med. Chem. 1991, 34, 3334–3338. [Google Scholar] [CrossRef]
  26. Sapra, S.; Bhalla, Y.; Nandani; Sharma, S.; Singh, G.; Nepali, K.; Budhiraja, A.; Dhar, K.L. Colchicine and its various physicochemical and biological aspects. Med. Chem. Res. 2012, 22, 531–547. [Google Scholar] [CrossRef]
  27. Zykova, S.S.; Kizimova, I.A.; Syutkina, A.I.; Toksarova, Y.S.; Igidov, N.M.; Ibishov, D.F.; Boichuk, S.V.; Dunaev, P.D.; Galembikova, A.R.; Korochkina, R.R. Synthesis and Cytostatic Activity of (E)-Ethyl-2-Amino-5-(3,3-Dimethyl-4-Oxobutyliden)-4-Oxo-1- (2-Phenylaminobenzamido)-4,5-Dihydro-1Hpyrrol-3-Carboxylate. Pharm. Chem. J. 2020, 53, 895–898. [Google Scholar] [CrossRef]
  28. Shelley, J.C.; Cholleti, A.; Frye, L.L.; Greenwood, J.R.; Timlin, M.R.; Uchimaya, M. Epik: A software program for pK a prediction and protonation state generation for drug-like molecules. J. Comput. Mol. Des. 2007, 21, 681–691. [Google Scholar] [CrossRef]
  29. Sherman, W.; Day, T.; Jacobson, M.P.; Friesner, R.A.; Farid, R. Novel Procedure for Modeling Ligand/Receptor Induced Fit Effects. J. Med. Chem. 2006, 49, 534–553. [Google Scholar] [CrossRef]
  30. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [Green Version]
  31. Lu, C.; Wu, C.; Ghoreishi, D.; Chen, W.; Wang, L.; Damm, W.; Ross, G.A.; Dahlgren, M.K.; Russell, E.; Von Bargen, C.D.; et al. OPLS4: Improving Force Field Accuracy on Challenging Regimes of Chemical Space. J. Chem. Theory Comput. 2021, 17, 4291–4300. [Google Scholar] [CrossRef] [PubMed]
  32. Fusani, L.; Palmer, D.S.; Somers, D.O.; Wall, I.D. Exploring Ligand Stability in Protein Crystal Structures Using Binding Pose Metadynamics. J. Chem. Inf. Model. 2020, 60, 1528–1539. [Google Scholar] [CrossRef] [PubMed]
  33. Bowers, K.J.; Chow, D.E.; Xu, H.; Dror, R.O.; Eastwood, M.P.; Gregersen, B.A.; Klepeis, J.L.; Kolossvary, I.; Moraes, M.A.; Sacerdoti, F.D.; et al. Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters. In Proceedings of the ACM/IEEE SC 2006 Conference (SC’06), Tampa, FL, USA, 11–17 November 2006; IEEE: Piscataway, NJ, USA, 2006; p. 43. [Google Scholar] [CrossRef] [Green Version]
  34. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  35. Grant, B.J.; Rodrigues, A.P.; ElSawy, K.M.; McCammon, J.A.; Caves, L.S. Bio3d: An R package for the comparative analysis of protein structures. Bioinformatics 2006, 22, 2695–2696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: https://www.R-project.org/ (accessed on 28 April 2022).
  37. Chaimontree, S.; Atkinson, K.; Coenen, F. Best Clustering Configuration Metrics: Towards Multiagent Based Clustering. In International Conference on Advanced Data Mining and Applications, Chongqing, China, 19–21 November 2010; Springer: Berlin/Heidelberg, Germany, 2010; pp. 48–59. [Google Scholar] [CrossRef] [Green Version]
  38. Li, J.; Abel, R.; Zhu, K.; Cao, Y.; Zhao, S.; Friesner, R.A. The VSGB 2.0 model: A next generation energy model for high resolution protein structure modeling. Proteins: Struct. Funct. Bioinform. 2011, 79, 2794–2812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bhardwaj, V.; Gumber, D.; Abbot, V.; Dhiman, S.; Sharma, P. Pyrrole: A resourceful small molecule in key medicinal hetero-aromatics. RSC Adv. 2015, 5, 15233–15266. [Google Scholar] [CrossRef]
  40. DeSimone, R.W.; Currie, K.S.; Mitchell, S.A.; Darrow, J.W.; Pippin, D.A. Privileged Structures: Applications in Drug Discovery. Comb. Chem. High Throughput Screen. 2004, 7, 473–493. [Google Scholar] [CrossRef]
  41. Duarte, C.D.; Barreiro, E.J.; Fraga, C.A.M. Privileged structures: A useful concept for the rational design of new lead drug candidates. Mini-Rev. Med. Chem. 2007, 7, 1108–1119. [Google Scholar] [CrossRef]
  42. Petri, G.L.; Spanò, V.; Spatola, R.; Holl, R.; Raimondi, M.V.; Barraja, P.; Montalbano, A. Bioactive pyrrole-based compounds with target selectivity. Eur. J. Med. Chem. 2020, 208, 112783. [Google Scholar] [CrossRef]
  43. Walsh, C.T.; Garneau-Tsodikova, S.; Howard-Jones, A.R. Biological formation of pyrroles: Nature’s logic and enzymatic machinery. Nat. Prod. Rep. 2006, 23, 517–531. [Google Scholar] [CrossRef]
  44. Ahmad, S.; Alam, O.; Naim, M.J.; Shaquiquzzaman, M.; Alam, M.M.; Iqbal, M. Pyrrole: An insight into recent pharmacological advances with structure activity relationship. Eur. J. Med. Chem. 2018, 157, 527–561. [Google Scholar] [CrossRef] [PubMed]
  45. Bianco, M.d.C.A.D.; Marinho, D.I.L.F.; Hoelz, L.V.B.; Bastos, M.M.; Boechat, N. Pyrroles as Privileged Scaffolds in the Search for New Potential HIV Inhibitors. Pharmaceuticals 2021, 14, 893. [Google Scholar] [CrossRef] [PubMed]
  46. La Regina, G.; Bai, R.; Coluccia, A.; Famiglini, V.; Pelliccia, S.; Passacantilli, S.; Mazzoccoli, C.; Ruggieri, V.; Sisinni, L.; Bolognesi, A.; et al. New Pyrrole Derivatives with Potent Tubulin Polymerization Inhibiting Activity As Anticancer Agents Including Hedgehog-Dependent Cancer. J. Med. Chem. 2014, 57, 6531–6552. [Google Scholar] [CrossRef] [PubMed]
  47. Jadala, C.; Prasad, B.; Prasanthi, A.V.G.; Shankaraiah, N.; Kamal, A. Transition metal-free one-pot synthesis of substituted pyrroles by employing aza-Wittig reaction. RSC Adv. 2019, 9, 30659–30665. [Google Scholar] [CrossRef] [Green Version]
  48. Tang, S.; Zhou, Z.; Jiang, Z.; Zhu, W.; Qiao, D. Indole-Based Tubulin Inhibitors: Binding Modes and SARs Investigations. Molecules 2022, 27, 1587. [Google Scholar] [CrossRef]
  49. Romagnoli, R.; Oliva, P.; Salvador, M.K.; Manfredini, S.; Padroni, C.; Brancale, A.; Ferla, S.; Hamel, E.; Ronca, R.; Maccarinelli, F.; et al. A facile synthesis of diaryl pyrroles led to the discovery of potent colchicine site antimitotic agents. Eur. J. Med. Chem. 2021, 214, 113229. [Google Scholar] [CrossRef]
  50. Sun, J.; Chen, L.; Liu, C.; Wang, Z.; Zuo, D.; Pan, J.; Qi, H.; Bao, K.; Wu, Y.; Zhang, W. Synthesis and biological evaluations of 1,2-diaryl pyrroles as analogues of combretastatin A-4. Chem. Biol. Drug Des. 2015, 86, 1541–1547. [Google Scholar] [CrossRef]
  51. Ma, Z.; Ma, Z.; Zhang, D. Synthesis of Multi-Substituted Pyrrole Derivatives Through [3+2] Cycloaddition with Tosylmethyl Isocyanides (TosMICs) and Electron-Deficient Compounds. Molecules 2018, 23, 2666. [Google Scholar] [CrossRef] [Green Version]
  52. Mowery, P.; Mejia, F.B.; Franceschi, C.L.; Kean, M.H.; Kwansare, D.O.; Lafferty, M.M.; Neerukonda, N.D.; Rolph, C.E.; Truax, N.J.; Pelkey, E.T. Synthesis and evaluation of the anti-proliferative activity of diaryl-3-pyrrolin-2-ones and fused analogs. Bioorganic Med. Chem. Lett. 2017, 27, 191–195. [Google Scholar] [CrossRef] [Green Version]
  53. Boichuk, S.; Galembikova, A.; Syuzov, K.; Dunaev, P.; Bikinieva, F.; Aukhadieva, A.; Zykova, S.; Igidov, N.; Gankova, K.; Novikova, M.; et al. The Design, Synthesis, and Biological Activities of Pyrrole-Based Carboxamides: The Novel Tubulin Inhibitors Targeting the Colchicine-Binding Site. Molecules 2021, 26, 5780. [Google Scholar] [CrossRef]
  54. Mohammad, T.; Siddiqui, S.; Shamsi, A.; Alajmi, M.F.; Hussain, A.; Islam, A.; Ahmad, F.; Hassan, M.I. Virtual Screening Approach to Identify High-Affinity Inhibitors of Serum and Glucocorticoid-Regulated Kinase 1 among Bioactive Natural Products: Combined Molecular Docking and Simulation Studies. Molecules 2020, 25, 823. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Anwar, S.; Shamsi, A.; Shahbaaz, M.; Queen, A.; Khan, P.; Hasan, G.M.; Islam, A.; Alajmi, M.F.; Hussain, A.; Ahmad, F.; et al. Rosmarinic Acid Exhibits Anticancer Effects via MARK4 Inhibition. Sci. Rep. 2020, 10, 10300. [Google Scholar] [CrossRef] [PubMed]
  56. Anwar, S.; DasGupta, D.; Shafie, A.; Alhumaydhi, F.A.; Alsagaby, S.A.; Shahwan, M.; Anjum, F.; Al Abdulmonem, W.; Sharaf, S.E.; Hassan, I. Implications of tempol in pyruvate dehydrogenase kinase 3 targeted anticancer therapeutics: Computational, spectroscopic, and calorimetric studies. J. Mol. Liq. 2022, 350, 118581. [Google Scholar] [CrossRef]
  57. Prota, A.E.; Danel, F.; Bachmann, F.; Bargsten, K.; Buey, R.M.; Pohlmann, J.; Reinelt, S.; Lane, H.; Steinmetz, M.O. The Novel Microtubule-Destabilizing Drug BAL27862 Binds to the Colchicine Site of Tubulin with Distinct Effects on Microtubule Organization. J. Mol. Biol. 2014, 426, 1848–1860. [Google Scholar] [CrossRef] [PubMed]
  58. Zhao, W.; Zhou, C.; Guan, Z.-Y.; Yin, P.; Chen, F.; Tang, Y.-J. Structural Insights into the Inhibition of Tubulin by the Antitumor Agent 4β-(1,2,4-triazol-3-ylthio)-4-deoxypodophyllotoxin. ACS Chem. Biol. 2017, 12, 746–752. [Google Scholar] [CrossRef] [PubMed]
  59. Niu, L.; Wang, Y.; Wang, C.; Wang, Y.; Jiang, X.; Ma, L.; Wu, C.; Yu, Y.; Chen, Q. Structure of 4′-demethylepipodophyllotoxin in complex with tubulin provides a rationale for drug design. Biochem. Biophys. Res. Commun. 2017, 493, 718–722. [Google Scholar] [CrossRef] [PubMed]
  60. Arnst, K.E.; Wang, Y.; Hwang, D.-J.; Xue, Y.; Costello, T.; Hamilton, D.; Chen, Q.; Yang, J.; Park, F.; Dalton, J.T.; et al. A Potent, Metabolically Stable Tubulin Inhibitor Targets the Colchicine Binding Site and Overcomes Taxane Resistance. Cancer Res. 2017, 78, 265–277. [Google Scholar] [CrossRef] [Green Version]
  61. Pallante, L.; Rocca, A.; Klejborowska, G.; Huczynski, A.; Grasso, G.; Tuszynski, J.A.; Deriu, M.A. 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. 2020, 8, 108. [Google Scholar] [CrossRef]
  62. Gallego-Yerga, L.; Ochoa, R.; Lans, I.; Peña-Varas, C.; Alegría-Arcos, M.; Cossio, P.; Ramírez, D.; Peláez, R. Application of ensemble pharmacophore-based virtual screening to the discovery of novel antimitotic tubulin inhibitors. Comput. Struct. Biotechnol. J. 2021, 19, 4360–4372. [Google Scholar] [CrossRef]
  63. Hagras, M.; El Deeb, M.A.; Elzahabi, H.S.A.; Elkaeed, E.B.; Mehany, A.B.M.; Eissa, I.H. Discovery of new quinolines as potent colchicine binding site inhibitors: Design, synthesis, docking studies, and anti-proliferative evaluation. J. Enzym. Inhib. Med. Chem. 2021, 36, 640–658. [Google Scholar] [CrossRef]
  64. Wu, X.; Wang, Q.; Li, W. Recent advances in heterocyclic tubulin inhibitors targeting the colchicine binding site. Anticancer Agents Med. Chem. 2016, 16, 1325–1338. [Google Scholar] [CrossRef] [PubMed]
  65. Sun, K.; Sun, Z.; Zhao, F.; Shan, G.; Meng, Q. Recent advances in research of colchicine binding site inhibitors and their interaction modes with tubulin. Futur. Med. Chem. 2021, 13, 839–858. [Google Scholar] [CrossRef]
  66. Li, W.; Sun, H.; Xu, S.; Zhu, Z.; Xu, J. Tubulin inhibitors targeting the colchicine binding site: A perspective of privileged structures. Futur. Med. Chem. 2017, 9, 1765–1794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Lu, Y.; Chen, J.; Xiao, M.; Li, W.; Miller, D.D. An Overview of Tubulin Inhibitors That Interact with the Colchicine Binding Site. Pharm. Res. 2012, 29, 2943–2971. [Google Scholar] [CrossRef] [Green Version]
  68. Dorléans, A.; Gigant, B.; Ravelli, R.B.G.; Mailliet, P.; Mikol, V.; Knossow, M. Variations in the colchicine-binding domain provide insight into the structural switch of tubulin. Proc. Natl. Acad. Sci. USA 2009, 106, 13775–13779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Peng, L.X.; Hsu, M.T.; Bonomi, M.; Agard, D.A.; Jacobson, M.P. The Free Energy Profile of Tubulin Straight-Bent Conformational Changes, with Implications for Microtubule Assembly and Drug Discovery. PLoS Comput. Biol. 2014, 10, e1003464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Molecules 27 02873 g001 550
Figure 1. The common synthetic pathway of the derivatives of 2-amino pyrroles.
Figure 1. The common synthetic pathway of the derivatives of 2-amino pyrroles.
Molecules 27 02873 g001
Molecules 27 02873 g002 550
Figure 2. Immunoblot analysis for the expression of M-phase specific proteins (e.g., phospho-NuMA Ser395 and phospho-H3 Ser10) in H1299 lung cancer (A) and HCC1806 breast cancer (B) treated with DMSO (negative control), EAPC-67, EAPC-70, EAPC-71 (both 10 μM), Colchicine (0.05 μM) and Vinblastine (0.01 μM) for 24 h. Actin stain was used as a loading control.
Figure 2. Immunoblot analysis for the expression of M-phase specific proteins (e.g., phospho-NuMA Ser395 and phospho-H3 Ser10) in H1299 lung cancer (A) and HCC1806 breast cancer (B) treated with DMSO (negative control), EAPC-67, EAPC-70, EAPC-71 (both 10 μM), Colchicine (0.05 μM) and Vinblastine (0.01 μM) for 24 h. Actin stain was used as a loading control.
Molecules 27 02873 g002
Molecules 27 02873 g003 550
Figure 3. Dynamics of tubulin polymerization in samples treated with EAPC-67 (A) and EAPC-70 (B). Tubulin was also incubated with DMSO (control), Paclitaxel (PTX), Vinblastine (Vin), and CaCl2 at 37 °C, and absorbance was assessed every min for 1 h. A shift of the curve to the upper left of the control (DMSO) represents an increase in polymerized microtubules. A shift to the bottom right of the graph reflects the decrease in the rate of tubulin polymerization.
Figure 3. Dynamics of tubulin polymerization in samples treated with EAPC-67 (A) and EAPC-70 (B). Tubulin was also incubated with DMSO (control), Paclitaxel (PTX), Vinblastine (Vin), and CaCl2 at 37 °C, and absorbance was assessed every min for 1 h. A shift of the curve to the upper left of the control (DMSO) represents an increase in polymerized microtubules. A shift to the bottom right of the graph reflects the decrease in the rate of tubulin polymerization.
Molecules 27 02873 g003
Molecules 27 02873 g004 550
Figure 4. Immunoblot analysis for apoptosis markers (e.g., cleaved forms of PARP and caspase-3) in H1299 lung cancer (A) and HCC1806 breast cancer (B) after treatment with DMSO (negative control), EAPC-67, EAPC-70, EAPC-71 (10 μM), colchicine (0.05 μM) and vinblastine (0.01 μM) for 24 h. Actin stain is used as a loading control.
Figure 4. Immunoblot analysis for apoptosis markers (e.g., cleaved forms of PARP and caspase-3) in H1299 lung cancer (A) and HCC1806 breast cancer (B) after treatment with DMSO (negative control), EAPC-67, EAPC-70, EAPC-71 (10 μM), colchicine (0.05 μM) and vinblastine (0.01 μM) for 24 h. Actin stain is used as a loading control.
Molecules 27 02873 g004
Molecules 27 02873 g005 550
Figure 5. FACs analysis for apoptotic markers in HCC 1806 breast cancer cells treated with DMSO (control), paclitaxel, vinblastine, and colchicine (positive control), EAPC-67, and EAPC-70 for 24 h. (A) Representative dot plots are shown. (B) Quantitative analysis of the early-apoptotic cells after the treatment as indicated above. (C) Quantitative analysis of the total apoptotic cells after the treatment as indicated above. **** p < 0.0001.
Figure 5. FACs analysis for apoptotic markers in HCC 1806 breast cancer cells treated with DMSO (control), paclitaxel, vinblastine, and colchicine (positive control), EAPC-67, and EAPC-70 for 24 h. (A) Representative dot plots are shown. (B) Quantitative analysis of the early-apoptotic cells after the treatment as indicated above. (C) Quantitative analysis of the total apoptotic cells after the treatment as indicated above. **** p < 0.0001.
Molecules 27 02873 g005
Molecules 27 02873 g006 550
Figure 6. Average RMSD of two EAPC-67 (A,B) tautomers conformation over 10 × 10 ns metadynamics runs with calculated values of PoseScore and PersScore.
Figure 6. Average RMSD of two EAPC-67 (A,B) tautomers conformation over 10 × 10 ns metadynamics runs with calculated values of PoseScore and PersScore.
Molecules 27 02873 g006
Molecules 27 02873 g007 550
Figure 7. Change in RMSD (root means square deviation) for ligand (pink color) and for Cα atoms of tubulin dimer (blue color) during 100 ns of unbiased molecular dynamic simulation.
Figure 7. Change in RMSD (root means square deviation) for ligand (pink color) and for Cα atoms of tubulin dimer (blue color) during 100 ns of unbiased molecular dynamic simulation.
Molecules 27 02873 g007
Molecules 27 02873 g008 550
Figure 8. (A) Protein-ligand fraction interaction diagram between EAPC-67 and tubulin. H-bonds, ionic, hydrophobic interactions, and water bridges are shown in green, red, purple, and blue, respectively. (B) Contact points of EAPC-67 and amino acids of colchicine binding site during the whole simulation trajectory. The top panel shows the total number of protein-ligand contacts over 100 ns.
Figure 8. (A) Protein-ligand fraction interaction diagram between EAPC-67 and tubulin. H-bonds, ionic, hydrophobic interactions, and water bridges are shown in green, red, purple, and blue, respectively. (B) Contact points of EAPC-67 and amino acids of colchicine binding site during the whole simulation trajectory. The top panel shows the total number of protein-ligand contacts over 100 ns.
Molecules 27 02873 g008
Molecules 27 02873 g009 550
Figure 9. 2D plot of ligand-protein fraction interaction diagram between EAPC-67 and tubulin over 100 nanoseconds MD with labeled key amino acids. Each type of non-covalent interaction is highlighted in the corresponding color.
Figure 9. 2D plot of ligand-protein fraction interaction diagram between EAPC-67 and tubulin over 100 nanoseconds MD with labeled key amino acids. Each type of non-covalent interaction is highlighted in the corresponding color.
Molecules 27 02873 g009
Molecules 27 02873 g010 550
Figure 10. Principal component analysis (PCA) of EAPC-67-tubulin complex trajectory. (A) Plotting PC1 component (x-axis) versus PC2 component (y-axis). Each dot denotes the one conformation of the protein. The spread of blue and red color dots describes the degree of conformational changes in the simulation, where the color scale from blue to white to red is equivalent to the simulation time. The blue indicates the initial timestep, white is intermediate and the final timestep is represented by the red color. (B) The proportion of variance per each principal component. (C) Clustering of the EAPC-67-tubulin complex conformations in the PC1-PC3 subspace based on the K-means method and projection of cluster members onto the PC1-PC2 plot diagram. A member of each of the three clusters is indicated by a separate color. (D) Conformational clusters of the EAPC-67-tubulin complex are displayed on the RMSD plot. Each dot represents a complex conformer, and the yellow markers indicate the representative conformations (medoids–snapshots 834, 1556, and 1987) of each cluster.
Figure 10. Principal component analysis (PCA) of EAPC-67-tubulin complex trajectory. (A) Plotting PC1 component (x-axis) versus PC2 component (y-axis). Each dot denotes the one conformation of the protein. The spread of blue and red color dots describes the degree of conformational changes in the simulation, where the color scale from blue to white to red is equivalent to the simulation time. The blue indicates the initial timestep, white is intermediate and the final timestep is represented by the red color. (B) The proportion of variance per each principal component. (C) Clustering of the EAPC-67-tubulin complex conformations in the PC1-PC3 subspace based on the K-means method and projection of cluster members onto the PC1-PC2 plot diagram. A member of each of the three clusters is indicated by a separate color. (D) Conformational clusters of the EAPC-67-tubulin complex are displayed on the RMSD plot. Each dot represents a complex conformer, and the yellow markers indicate the representative conformations (medoids–snapshots 834, 1556, and 1987) of each cluster.
Molecules 27 02873 g010
Molecules 27 02873 g011 550
Figure 11. The structural difference between straight (PDB ID: 1JFF) and representative conformations of β-chain of tubulin was obtained from MD. The straight conformation (light brown color) of the β-chain of tubulin is aligned with the representative conformation (blue color). Arrows indicate the main secondary elements undergoing significant structural changes during the transition from one conformation to another. The straight confirmation is not suitable for interaction with the EAPC-67, since steric clusters are formed with the amino acids of the protein secondary elements (S8, S9, H7, T7) surrounding the colchicine-binding pocket.
Figure 11. The structural difference between straight (PDB ID: 1JFF) and representative conformations of β-chain of tubulin was obtained from MD. The straight conformation (light brown color) of the β-chain of tubulin is aligned with the representative conformation (blue color). Arrows indicate the main secondary elements undergoing significant structural changes during the transition from one conformation to another. The straight confirmation is not suitable for interaction with the EAPC-67, since steric clusters are formed with the amino acids of the protein secondary elements (S8, S9, H7, T7) surrounding the colchicine-binding pocket.
Molecules 27 02873 g011
Table
Table 1. IC50 values (in the micromolar range) for EAPCs in HCC1806, MDA-MB-231 triple-negative breast cancer, and H1299 non-small cell lung cancer cell lines.
Table 1. IC50 values (in the micromolar range) for EAPCs in HCC1806, MDA-MB-231 triple-negative breast cancer, and H1299 non-small cell lung cancer cell lines.
Number of EAPCCell Lines
MDA-MB 231H1299HCC 1806
6421.5 ± 2.0128.7 ± 1.5423.19 ± 3.72
6528.2 ± 1.1217.5 ± 2.3114.13 ± 2.91
6640.2 ± 3.4246.3 ± 4.2154.2 ± 4.87
672.9 ± 0.216.7 ± 0.539.4 ± 0.64
68>100>10056.16 ± 3.2
6926.1 ± 2.5621.5 ± 1.529.50 ± 1.43
706.7 ± 0.522.2 ± 0.2525.35 ± 2.12
7110.24 ± 1.362.5 ± 0.4229.73 ± 1.91
7285.1 ± 1.2579.6 ± 2.568.63 ± 4.17
7339.3 ± 2.97>100>100
Table
Table 2. ΔGbind energy and ligand strain energy values for each medoid of three clusters.
Table 2. ΔGbind energy and ligand strain energy values for each medoid of three clusters.
Medoids of Cluster MM-GBSA ΔG bind
(kcal/mol) 		Ligand Strain Energy
(kcal/mol)
Snapshot 834−24.211.4
Snapshot 1556−31.19.3
Snapshot 1987−25.610.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boichuk, S.; Syuzov, K.; Bikinieva, F.; Galembikova, A.; Zykova, S.; Gankova, K.; Igidov, S.; Igidov, N. Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin. Molecules 2022, 27, 2873. https://doi.org/10.3390/molecules27092873

AMA Style

Boichuk S, Syuzov K, Bikinieva F, Galembikova A, Zykova S, Gankova K, Igidov S, Igidov N. Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin. Molecules. 2022; 27(9):2873. https://doi.org/10.3390/molecules27092873

Chicago/Turabian Style

Boichuk, Sergei, Kirill Syuzov, Firuza Bikinieva, Aigul Galembikova, Svetlana Zykova, Ksenia Gankova, Sergei Igidov, and Nazim Igidov. 2022. "Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin" Molecules 27, no. 9: 2873. https://doi.org/10.3390/molecules27092873

APA Style

Boichuk, S., Syuzov, K., Bikinieva, F., Galembikova, A., Zykova, S., Gankova, K., Igidov, S., & Igidov, N. (2022). Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin. Molecules, 27(9), 2873. https://doi.org/10.3390/molecules27092873

Article Metrics

Back to TopTop