Ramya Tokala1, Sravani Sana1, Uppu Jaya Lakshmi1, Prasanthi Sankarana1, Dilep Kumar Sigalapalli1, Nikhil Gadewal2, Jyoti Kode3, Nagula Shankaraiah4. 1. Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), Hyderabad 500 037, India. 2. Bioinformatics Centre, Advanced Centre for Treatment, Research & Education in Cancer, Tata Memorial Centre, Kharghar, Navi Mumbai 410210, India. 3. Tumor Immunology & Immunotherapy Group, Advanced Centre for Treatment, Research & Education in Cancer, Tata Memorial Centre, Kharghar, Navi Mumbai 410210, India; Homi-Bhabha National Institute (HBNI), Training School Complex, Anushakti Nagar, Mumbai 400094, India. 4. Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), Hyderabad 500 037, India. Electronic address: shankar@niperhyd.ac.in.
Abstract
The conjoining of salient pharmacophoric properties directing the development of prominent cytotoxic agents was executed by constructing thiadiazolo-carboxamide bridged β-carboline-indole hybrids. On the evaluation of in vitro cytotoxic potential, 12c exhibited prodigious cytotoxicity among the synthesized new molecules 12a-k, with an IC50 < 5 μM in all the tested cancer cell lines (A549, MDA-MB-231, BT-474, HCT-116, THP-1) and the best cytotoxic potential was expressed in lung cancer cell line (A549) with an IC50 value of 2.82 ± 0.10 μM. Besides, another compound 12a also displayed impressive cytotoxicity against A549 cell line (IC50: 3.00 ± 1.40 μM). Further target-based assay of these two compounds 12c and 12a revealed their potential as DNA intercalative topoisomerase-IIα inhibitors. Additionally, the antiproliferative activity of compound 12c was measured in A549 cells by traditional apoptosis assays revealing the nuclear, morphological alterations, and depolarization of membrane potential in mitochondria and externalization of phosphatidylserine in a concentration-dependent manner. Cell cycle analysis unveiled the G0/G1 phase inhibition and wound healing assay inferred the inhibition of in vitro cell migration by compound 12c in lung cancer cells. Remarkably, the safety profile of compound 12c was disclosed by screening against normal human lung epithelial cell line (BEAS-2B: IC50: 71.2 ± 7.95 μM) with a selectivity index range of 14.9-25.26. Moreover, Molecular modeling studies affirm the intercalative binding of compound 12c and 12a in the active pocket of topo-IIα. Furthermore, in silico prediction of physico-chemical parameters divulged the propitious drug-like properties of the synthesized derivatives.
The conjoining of salient pharmacophoric properties directing the development of prominent cytotoxic agents was executed by constructing thiadiazolo-carboxamide bridged β-carboline-indole hybrids. On the evaluation of in vitro cytotoxic potential, 12c exhibited prodigious cytotoxicity among the synthesized new molecules 12a-k, with an IC50 < 5 μM in all the tested cancer cell lines (A549, MDA-MB-231, BT-474, HCT-116, THP-1) and the best cytotoxic potential was expressed in lung cancer cell line (A549) with an IC50 value of 2.82 ± 0.10 μM. Besides, another compound 12a also displayed impressive cytotoxicity against A549 cell line (IC50: 3.00 ± 1.40 μM). Further target-based assay of these two compounds 12c and 12a revealed their potential as DNA intercalative topoisomerase-IIα inhibitors. Additionally, the antiproliferative activity of compound 12c was measured in A549 cells by traditional apoptosis assays revealing the nuclear, morphological alterations, and depolarization of membrane potential in mitochondria and externalization of phosphatidylserine in a concentration-dependent manner. Cell cycle analysis unveiled the G0/G1 phase inhibition and wound healing assay inferred the inhibition of in vitro cell migration by compound 12c in lung cancer cells. Remarkably, the safety profile of compound 12c was disclosed by screening against normal human lung epithelial cell line (BEAS-2B: IC50: 71.2 ± 7.95 μM) with a selectivity index range of 14.9-25.26. Moreover, Molecular modeling studies affirm the intercalative binding of compound 12c and 12a in the active pocket of topo-IIα. Furthermore, in silico prediction of physico-chemical parameters divulged the propitious drug-like properties of the synthesized derivatives.
Promising anticancer agents should be necessarily developed as cancer remains as public health concern worldwide [1] along with the additional origin of a newer pandemic like COVID-19 [2]. The drug repurposing strategy has been employed to reduce COVID-19 induced inflammation, coagulopathy and dysfunction by using anticancer agents like interleukin inhibitors, corticosteroids, kinase inhibitors, etc [2c]. The types, mutations [3], metastasis [4], resistance [5], and after-effects of clinical treatment, make it a gruesome disease [6] without ideal treatment. Consequently, there lies a lacuna in the progression of novel agents to combat cancer. Inhibition of more than one target with a key role in cell replication and proliferation set forth magnificent directions in the generation of new chemical entities (NCEs) towards chemical research.The discovery and comprehension of the structure and function of DNA are paramount of all the discoveries in the last century [7]. Hitherto, DNA is a major pharmacological target for small molecules which are effective in various ailments which originated from virus [8], bacteria [9], (a), (b), (c) and also in cancer [10]. Based on the interaction of the small molecules with DNA, they have been broadly divided into (i) Covalent inhibitors, which are irreversible binders that form covalent adducts, thereby cause serious damage to the DNA process and induce cell death. e.g., cisplatin (ii) Non-covalent inhibitors are reversible binders; further classified based on modes of binding with duplex DNA, (a) Surface binders (exhibits electrostatic interactions) (b) Intercalators (between DNA base pairs) (c) Groove binders (major and minor groove binding) [11], (a), (b), (c). Usually, the small molecules bind with DNA with more than one binding mode. Also, enzymes acting on DNA play a promising role in replication and repair. One such enzyme is DNA Topoisomerase II (Topo II), with a potential role in cell division and demarcation by causing nicking and relegation of DNA and is an imperative target in drug discovery [12]. There are two isoforms of Topo II with gene codes Topo IIα and Topo IIβ which are analogous in structure and are uniform in catalytic properties. However, topo IIα has positive coordination on the mitotic process of cells such as rate of cell proliferation, expression of mRNA, and also on different phases of the cell cycle nevertheless, cellular functions of topo IIβ are yet unclear [13], (a), (b), (c), (d), (e), (f). Topo II, in general, modifies DNA topology and induces double-strand breaks through multiple steps that require ATP hydrolysis and Mg2+ for activity [14a] and agents inhibiting it are categorized into interfacial poisons (IFPs) and catalytic inhibitory compounds (CICs). Anticancer agents like doxorubicin and etoposide are renowned Topo II poisons that cause accretion of deplorable linear DNA by inhibition of the relegation of cleaved DNA strands through the stabilization of the Topo-DNA complex. Catalytic inhibitors such as novobiocin and merbarone, block binding and nicking of DNA respectively. On the contrary, agents such as suramin, a potent catalytic inhibitor block topo II’s ATP binding site [14], (b), (c). In addition to the above features, agents that induce apoptosis are considered as remarkable anticancer agents as they target an important hallmark of cancer [15]. Thus the discovery of new chemical entities targeting topo II through DNA intercalation is an important area of cancer research. Based on the literature, the rationale to exhibit intercalative topo II inhibition by NCE’s is to possess a chromophore (head part) linked to a side chain (tail part) through a cationic species. The mechanistic details indicate that the chromophore is a polyaromatic ring system which intercalates between DNA base pairs. The linker and extending heterocycles impart stabilization by interacting with phosphate groups and by orienting into the groove of DNA respectively [16].The structural features of naturally occurring β-carboline alkaloids make it a potential chromophore to intercalate with DNA in between the base pairs. The tricyclic pyrido [3,4-b] ring system (norharmane A) possesses an affinity to stack between GC residues of DNA [17]. Furthermore, these Peganum harmala derived agents (harmine B) followed up by their chemical functionalization were reported as DNA damaging and photocleavage agents [18], (a), (b), (c). In synergism, inhibitions of other enzymes like topoisomerases [12], [19], (a), (b), kinases [20], (a), (b), (c), (d) make β-carboline a promising entity in cancer research along with indications in several diseases [21], (a), (b), (c).1,3,4-Thiadiazole is another privileged heterocyclic scaffold with excellent therapeutic potential against diverse health problems. The potential of this bioactive scaffold was far recognized as an antiepileptic (acetazolamide C), antibacterial, antiprotozoal (megazol D), and to treat glaucoma [22]. Farooqi et al., revealed the affinity of 1,3,4-thiadiazole towards DNA interaction by hybridization of 1,3,4-thiadiazole to the clinical drugs (Ibuprofen and Ciprofloxacin)[23]. Plech et. al., proved the affinity of 1, 3, 4-thiadiazole towards topo II by disubstitution on the heterocycles [24]. Indeed, this entity reveals wide pharmacological activities such as antimicrobial, antitubercular, anti-inflammatory, antiviral, anticonvulsant, and activity against central nervous system (CNS) disorders, cardiovascular system (CVS) diseases and cancer [22].Indole is a versatile and privileged scaffold with an extensive range of biological actions principally in cancer by targeting diverse enzymes (Topoisomerases, Tubulin polymerases, Aromatase) pathways (NF-kB/mTOR/PI3K/AKT), and mechanisms like DNA binding and apoptosis [25], (a), (b), (c), (d), (e). Naturally occurring antitumor antibiotic like duocarmycin has a strong affinity for DNA interaction, nonetheless unable to reach the clinic. However, synthetic analogs like adozelesin with substitution on 2nd position (E) reached clinical trials in the treatment of cancer. Highly effective D-64131(F) and D 68144 (G) were recognized as potent anticancer agents even in paclitaxel-resistant tumors [26] (Fig. 1
).
Fig. 1
Representative structures bearing pharmacologically active β-carbolines, thiadiazole, and indole nucleus.
Representative structures bearing pharmacologically active β-carbolines, thiadiazole, and indole nucleus.
Rationale
In maintaining our persistence to achieve the goal of reaching successful clinical candidates towards cancer [27], by considering the advantage of structural features of the chemically diverse scaffolds like β-carboline, 1,3,4-thiadiazole, indole and exploiting pharmacophore hybridization to carboxamide formation, we have methodologically designed a new series (12a–k) of DNA intercalative topo II inhibitors with classical apoptosis-inducing features. The rationale for this design has been established based on literature and structural functionalities as illustrated in Fig. 2
. Foremost requirements for intercalative topoisomerase II inhibition are the presence of a chromophore (intercalation into DNA), adjoining cationic center (electrostatic interaction with base pairs), and a groove extending side-chain/heterocycles (stabilization). DNA intercalative topoisomerase II inhibitors such as doxorubicin and PD-115934 also possess this kind of arrangement (chromophore + cationic center + groove extending side chain) [16] and support the rationale to this design. The chemical synthesis of β-carboline acid from amino acid -tryptophan satisfies the utmost feature as a chromophore for dual activity. Further, the cationic center of 1,3,4-thiadiazole on the 2nd position of indole (extending heterocycles) has been constructed by utilizing chemical functionalities of indole-2-acid and thiosemicarbazide followed by cyclization resulting in the free amino moiety. The chromophore has been adjoined with the rest of the scaffold by pharmacophore hybridization [28], (a), (b)
via carboxamide formation through acid amine coupling.
Fig. 2
Rationale in the design of target molecules (thiadiazolo-carboxamide bridged β-carboline-indole derivatives 12a–k).
Rationale in the design of target molecules (thiadiazolo-carboxamide bridged β-carboline-indole derivatives 12a–k).
Chemistry
The desired hybrids 1,3,4-thiadiazolo carboxamide bridged β-carboline-indole hybrids 12a–k were synthesized via amide formation between β-carboline acid and indol-2-yl-1,3,4-thiadiazol-2-amine as depicted in Scheme 2. 1-Aryl-9H-pyrido[3,4-b] indole-3-carboxylic acids 6a–f were synthesized by earlier reported procedures, [12] from commercially available -tryptophan (Spectrochem; CAS NO: 73-22-3) 1 which was esterified 2 and subjected to Pictet-Spengler condensation using acid catalyst trifluoroacetic acid with substituted aromatic benzaldehydes 3a–f to obtain tetrahydro-β-carboline esters 4a–f. The β-carboline esters 5a–f were obtained by further aromatization of 4a–f using KMnO4 and were subjected to basic hydrolysis to yield the respective β-carboline acids 6a–f (Scheme 1
).
Scheme 2
Synthesis of thiadiazolo-carboxamide bridged β-carboline-indole hybrids 12a–k.
Scheme 1
Synthesis of β-carboline component 6a–f (1-aryl-9H-pyrido[3,4-b]indole-3-carboxylic acids).
Synthesis of β-carboline component 6a–f (1-aryl-9H-pyrido[3,4-b]indole-3-carboxylic acids).Synthesis of thiadiazolo-carboxamide bridged β-carboline-indole hybrids 12a–k.The amine components 11a–c were attained by esterification of indole-2-acid 7, which (8) was further alkylated using different alkyl halides to yield alkylated indole-2-esters 9a–c. Thus obtained esters, on basic hydrolysis result in respective alkylated indole-2-acids 10a–c. This acid component undergoes cyclization with thiosemicarbazide in POCl3 to succumb amine counterpart. Further, the title compounds 12a–k was achieved by acid–amine coupling between 6a–f and 11a–c in good yields (Scheme 2). All these compounds were purified by using column chromatography.The representative compound 12a of synthesized derivatives was thus characterized by 1H NMR spectroscopy. The appearance of singlet at δ = 12.14 corresponds to the amidic N—H proton and confirms the formation of titled compounds. The two singlets at δ = 12.58 and 12.19 ppm account for the NH proton of indole and 9H-indolic proton of β-carboline ring respectively. The proton at the C4 position of β-carboline appeared as a singlet at δ = 9.13 ppm. All the other 14 aromatic protons of both β-carboline and indole are accounted for in the range of δ = 8.55 to 7.13 ppm. Further analysis by 13C NMR spectroscopy disclosed that amidic carbonyl carbon resonated at δ = 164.1 ppm. The signals at δ = 158.2 and 158.0 ppm account for amide attached C2 carbon and indole attached C5 carbon of 1,3,4-thiadiazole respectively. All the remaining carbons of 12a appeared in the region of δ = 142.1–104.5 ppm. Next, IR spectra show the band at 1673 cm−1 which indicates the CO stretching of the amide of the compound 12a. The HRMS exhibited specific [M+H]+ peaks that matched with their molecular formula. Similarly, spectral data (1H, 13C, FT-IR) of all the remaining compounds 12b–k matched accurately with their respective structures.
Biological evaluation
In vitro cytotoxicity assay
Thus, the synthesized compounds 12a–k were evaluated for their cytotoxic potential against a panel of human cancer cell lines and compared with the standard harmine and doxorubicin, VP-16 (etoposide) and represented in Table 1
. The MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] assay [29], (a) was performed for all the derivatives on different cancer cell lines such as (i) adherent cell lines viz lung cancer (A549), breast cancer (MDA-MB-231, BT-474), colon cancer (HCT-116), and (ii) suspension cell lines viz leukemia (HL-60) and acute monocytic leukemia (THP-1) (Figs. S1 and S2 of ESI) and compared to normal human cell line (BEAS-2B). Notably, all the tested compounds exhibited potent cytotoxicity against all the tested cell lines, particularly 12c with the highest cytotoxic potential of IC50 of 2.82 ± 0.10 µM in the A549 cell line. Further, it has leading cytotoxicity in all the cell lines tested with an IC50 of 3.18 ± 0.49 µM and 3.83 ± 0.94 µM in THP-1 and HL-60 respectively. Also, in both breast cancer cell lines, 12c expressed an IC50 of 3.28 ± 0.43 µM and 4.09 ± 0.59 µM in MDA-MB-231 (triple-negative) and BT-474. Also, the compound 12c has exhibited in vitro cytotoxicity in HCT-116 with an IC50 value of 4.78 ± 1.29 µM and the order of cytotoxicity for compound 12c in different cell lines is A549 > THP-1 > MDA-MB-231 > HL-60 > BT-474 > HCT-116.
Table 1
Cytotoxic potential (IC50 values in μM) a of thiadiazolo-carboxamide-β-carboline-indole derivatives 12a–k against human cell lines.
Entry
A549 b
MDA-MB-231c
BT-474 d
HCT-116 e
HL-60f
THP-1g
BEAS-2B h
12a
3.00 ± 1.40
3.94 ± 0.31
5.68 ± 0.58
5.12 ± 1.37
8.25 ± 1.46
3.57 ± 0.55
ND
12b
9.14 ± 2.87
7.93 ± 0.33
17.55 ± 1.06
10.33 ± 1.48
33.25 ± 2.54
14.23 ± 1.25
ND
12c
2.82 ± 0.10
3.28 ± 0.43
4.09 ± 0.59
4.78 ± 1.29
3.83 ± 0.94
3.18 ± 0.49
71.25 ± 7.95
12d
8.20 ± 1.74
7.85 ± 1.77
14.08 ± 1.88
8.07 ± 0.55
30.20 ± 0.54
14.45 ± 1.17
ND
12e
14.21 ± 1.98
9.39 ± 0.52
20.13 ± 1.19
19.89 ± 1.93
46.21 ± 1.68
48.59 ± 2.11
ND
12f
7.85 ± 1.77
7.23 ± 1.43
9.15 ± 2.68
6.11 ± 2.04
10.01 ± 0.50
8.78 ± 0.27
ND
12g
10.88 ± 1.54
8.66 ± 0.30
19.47 ± 0.74
18.72 ± 2.39
43.21 ± 1.33
40.22 ± 1.21
ND
12h
3.64 ± 0.98
6.79 ± 0.57
9.20 ± 2.47
6.22 ± 0.99
9.85 ± 1.90
6.11 ± 0.43
ND
12i
9.45 ± 0.27
8.43 ± 1.30
18.36 ± 1.22
14.34 ± 2.62
37.25 ± 1.01
32.84 ± 0.97
ND
12j
3.33 ± 1.60
4.45 ± 0.91
8.68 ± 1.49
5.98 ± 0.34
9.31 ± 1.08
4.38 ± 1.28
ND
12k
9.51 ± 0.33
8.57 ± 0.40
19.53 ± 1.12
15.09 ± 3.79
39.62 ± 0.23
38.89 ± 1.57
ND
Hari
19.41 ± 0.66
16.74 ± 0.79
20.22 ± 0.98
25.12 ± 1.18
ND
ND
ND
Doxj
1.92 ± 0.86
1.48 ± 0.56
1.82 ± 0.98
ND
1.08 ± 2.02
3.25 ± 2.20
ND
VP-16k
2.00 ± 1.47
ND
ND
ND
ND
ND
ND
ND: Not determined.
50% Inhibitory concentration after 72 h of drug treatment; all values are expressed as mean ± SEM, and each treatment was performed in triplicate wells.
Lung cancer.
Breast cancer.
Colon cancer.
Leukemia.
Acute monocytic leukemia.
Human bronchial lung epithelial cells. ijkReference compound.
Harmine.
Doxorubicin.
VP-16 = Etoposide
Cytotoxic potential (IC50 values in μM) a of thiadiazolo-carboxamide-β-carboline-indole derivatives 12a–k against human cell lines.ND: Not determined.50% Inhibitory concentration after 72 h of drug treatment; all values are expressed as mean ± SEM, and each treatment was performed in triplicate wells.Lung cancer.Breast cancer.Colon cancer.Leukemia.Acute monocytic leukemia.Human bronchial lung epithelial cells. ijkReference compound.Harmine.Doxorubicin.VP-16 = EtoposideEncouragingly, the compound 12c displayed the lowest in vitro cytotoxicity with an IC50 value of 71.25 ± 7.95 µM in normal human lung epithelial cells (BEAS-2B). Remarkably, the compounds 12a, 12f, 12h, and 12j exhibited notable cytotoxicity potential in all the cell lines tested with overall IC50 ≤ 10 µM. Moreover, compounds 12b, 12d, 12i, and 12k also expressed a reasonable in vitro cytotoxicity on A549 and MDA-MB-231 cell lines with an IC50 value < 10 µM. From the detailed analysis of the in vitro cytotoxicity table, it was comprehended that 12c is a potent compound and possesses the least inhibition in BEAS-2B and the cytospecificity profile on all the human cancer cell lines along with normal human cell line was represented in graphical form in Fig. 3a
.
Fig. 3a
Cytospecificity exhibited by compound 12c towards human cancer cell lines in comparison to normal human lung epithelial cells (BEAS-2B).
Cytospecificity exhibited by compound 12c towards human cancer cell lines in comparison to normal human lung epithelial cells (BEAS-2B).
Selectivity index
After contemplating the specificity exhibited by potent compound 12c in different cancer cell lines and normal human cell line (Fig. 3a), the selectivity index (SI) was calculated and portrayed in Fig. 3b
. The selectivity index is the ratio of IC50 in the non-tumor cell line to IC50 in the tumor cell line [29b] and a compound is said to be selective only if the value of SI ≥ 10 [29c]. It was observed that compound 12c exhibited selectivity in all the cell lines tested with selectivity range of 25.26–14.90 in the following order A549 > THP-1 > MDA-MB-231 > HL-60 > BT-474 > HCT-116.
Fig. 3b
Selectivity index exhibited by compound 12c in different human cancer cell lines.
Selectivity index exhibited by compound 12c in different human cancer cell lines.
Structure-Activity Relationship (SAR)
Based on the established in vitro cytotoxicity data against human cancer cell lines, the structure–activity relationship (SAR) was constructed to provide insights regarding the pattern of efficiency of the synthesized NCEs as illustrated in Fig. 4
. Two-point diverse modifications have been performed on the C1-phenyl of the β-carboline (R1) and nitrogen of indole (R2). Noticeably, the compounds 12a, 12c, 12h, and 12j without substitution on indole (R2 = H) were found to be potent than substituted derivatives 12b, 12d–g, 12i, and 12k. Additionally, increase in length of alky chain (R2 = H > Me > Et) on indole resulted in derivatives with lower cytotoxicity potential (12c > 12d > 12e & 12a > 12b). The better cytotoxicity was observed when R1 is substituted by electron-donating group (EDG) (12c–g, 12k: R1 = 4-Me, 4-OMe & 3,4,5-triOMe) than no substitution (12a, b: R1 = H) and followed by electron-withdrawing group (EWG) (12h–j: R1 = 4-F & 4-Cl). The in vitro cytotoxic potential was compared with EDG and EWG in the following order: (a) EDG: (i) In derivatives where R2 = methyl; and R1 is 3,4,5-trimethoxy substituted, such derivative exhibited superior inhibition than 4-methoxy containing corresponding derivative (12f > 12d); (ii) In compounds with R2 = ethyl; the compound with 4-methyl substitution on R1 displayed better inhibition than derivative with 3,4,5-trimethoxy followed by 4-methoxy substituted compound (12k>12g > 12e); (b) EWG: In compounds with R2 = H; the derivative with 4-chloro substitution at R1 position exhibited higher inhibition than 4-fluoro substituted derivative (12j > 12h). Interestingly, a similar pattern of inhibition was observed for all the compounds in the tested cell lines.
Fig. 4
Structure Activity Relationship (SAR) developed on target molecules 12a–k.
Structure Activity Relationship (SAR) developed on target molecules 12a–k.
Target-based assays
Topoisomerase IIα inhibition
Topoisomerase IIα is located in the nucleus and functions as a DNA nick inducing enzyme with established importance in DNA replication. It induces nicks on both the DNA strands, religates the phosphodiester bonds, and thereby alters the DNA topology. Thus, inhibition of topoisomerase IIα results in either blockage of the enzyme to induce nicks which thereby interferes the cell replication process (CICs) or prevent the relegation of nicks caused by the enzyme by stabilization of DNA-ternary complex which results in the initiation of recombination pathways (IFPs) and ultimately results in cell death [30], (a), (b), (c). Topo-IIα inhibitory assay has been performed using Topo-IIα drug screening kit (TG1019-1, TopoGEN, USA) with enzyme human topoisomerase IIα (TG2000H-2, TopoGEN, USA). In this assay, the two most potent compounds 12c and 12a from IC50 data were evaluated for their topo-IIα inhibitory activity. The catenated kDNA and topo II (5 units) along with compound 12c (10 μM) and 12a (10 μM) (lane 1 and 2 respectively) didn’t result in the formation of any linear DNA (a & b). Further, decatenated kDNA and linear DNA in lines 3 and 4 were the reference markers provided with the screening kit. From Fig. 5
, it is evident that the compound 12c and 12a inhibited topo IIα to induce nicks and linearize kDNA and fall into the category of catalytic inhibitory compounds.
Fig. 5
Inhibition of enzyme topo IIα by compound 12c and 12a using horizontal gel electrophoresis. (a) Lane 1: kDNA + Topo IIα + 12c (b) Lane 2: kDNA + Topo IIα + 12a (c) Lane 3: DNA marker = Decatenated kDNA (d) Lane 4: DNA marker = Linear DNA. The DNA markers were obtained directly along with the kit (TG1019-1, TopoGEN, USA).
Inhibition of enzyme topo IIα by compound 12c and 12a using horizontal gel electrophoresis. (a) Lane 1: kDNA + Topo IIα + 12c (b) Lane 2: kDNA + Topo IIα + 12a (c) Lane 3: DNA marker = Decatenated kDNA (d) Lane 4: DNA marker = Linear DNA. The DNA markers were obtained directly along with the kit (TG1019-1, TopoGEN, USA).
DNA intercalation assays
Intercalation of the ligand into DNA doesn’t always cause alterations in hydrogen bonds between base pairs, but the ligand-DNA complex is stabilized by van der Waals, charge transfer, π-π stacking with the planar scaffold, hydrophobic and hydrogen-bonding interactions [31], (a), (b). Intercalation stiffens, lengthens, reduces helical twist, and unwinds the DNA double helix [32], (a), (b), thereby resulting in inhibition of cell proliferation and triggers cell death. [33] Among various classes of DNA binding agents, intercalating agents were widely studied by considering different biophysical assays. By relating absorbance and viscosity measurements using CT-DNA with a concentration of the tested compound, these assays suggest the mode of binding of the novel synthesized compounds.
Relative viscosity experiment
To understand the binding mode of newly synthesized series, a relative viscosity experiment has been performed for potent compounds 12c and 12a using CT-DNA as per the reported protocol [34]. The different modes of binding of ligand result in alterations in the viscosity of DNA. The increase in viscosity is a result of an increase in the axial length of DNA to accommodate the ligand in between the base pairs, which indicates the intercalative mode of binding. However, no or minimal alteration in viscosity is observed with compounds possessing minor groove binding affinity. Conversely, covalent binders such as cisplatin result in a decrease of viscosity due to absolute functional inhibition of DNA through irreversible binding [35], (a), (b). Ethidium bromide (EtBr) and Hoechst 33258 generally cause an increase and comparative alteration in relative viscosity of DNA indicating their mode of binding as intercalator and minor groove binder respectively [34].In the current study, different concentrations of compound 12c and 12a along with harmine, Hoechst 33258, and EtBr are incubated with CT-DNA, and changes in viscosity were recorded. The controls EtBr and harmine displayed an increase in viscosity whereas; Hoechst 33248 exhibited a minimal change in viscosity with the increase in concentration. The tested compounds exhibited an increase in viscosity in a concentration-dependent manner indicating intercalation as their mode of DNA binding. Therefore, a graph is plotted with the concentration of compound [compound/CT-DNA] on X-axis and viscosity [η/ηo]1/3 on Y-axis and is represented in Fig. 6
.
Fig. 6
Graphical representation of relative viscosity measurements of different compounds 12c and 12a along with different standards with CT-DNA. Hoechst33258 is negative control and EtBr: Ethidium Bromide is a positive intercalative standard.
Graphical representation of relative viscosity measurements of different compounds 12c and 12a along with different standards with CT-DNA. Hoechst33258 is negative control and EtBr: Ethidium Bromide is a positive intercalative standard.
Absorbance spectroscopy
Intercalation of ligands with DNA results in a slight bathochromic and significant hypochromic shift in absorbance at the specific wavelength [36]. This principle has been applied to find the efficacy of compounds 12c and 12a using CT-DNA.
Nanodrop method
Nanodrop is a kind of UV-spectrophotometer which uses a small quantity (nano quantity) to measure the quality of the nucleic acids [37]. The specific absorbance was measured at 260/230 nm wavelength with the different concentrations of compound 12c and 12a along with control CT-DNA. Harmine, Dox, Hoechst-33258, and EtBr were used as standards. It is evident from Fig. 7
A and 7B, that compound 12c and 12a caused hypochromic shift at specific wavelength respectively, and are potent intercalators.
Fig. 7
Spectrophotometric analysis was performed by nanodrop spectrophotometer using CT-DNA with (A) Test compound: 12c (B) Test compound: 12a. EtBr and Dox were standard intercalators and control DNA is CT-DNA without any ligand. Harmine is a representative β-carboline standard.
Spectrophotometric analysis was performed by nanodrop spectrophotometer using CT-DNA with (A) Test compound: 12c (B) Test compound: 12a. EtBr and Dox were standard intercalators and control DNA is CT-DNA without any ligand. Harmine is a representative β-carboline standard.
UV–Visible spectroscopy
In addition to the nanodrop method, the absorbance alterations induced by compounds 12a and 12c on CT-DNA were evaluated at a specific concentration along with potent intercalator EtBr as a positive control. The plain CT-DNA without the addition of any ligand was used as a control, and from Fig. 8
, it is evident that the order of absorption is EtBr > 12c > 12a. The results of the absorbance spectroscopic methods are consistent with relative viscosity measurements and establish the intercalative mode of binding of compounds tested with CT-DNA. [38]
Fig. 8
Absorbance spectra of compound 12c and 12a with CT-DNA. EtBr is standard, whereas control is CT-DNA alone in the absence of ligand.
Absorbance spectra of compound 12c and 12a with CT-DNA. EtBr is standard, whereas control is CT-DNA alone in the absence of ligand.
Determination of apoptosis
Programmed cell death is evaded by cancer cells and prevails in the different organs; however, new cancer agents that induce apoptosis are always the prime choice of cancer research [39]. To understand the ability of compound 12c in apoptosis induction, various qualitative (morphological) and quantitative assays were performed on the A549 lung cancer cell line.
Determination of morphological changes
The morphological assays reveal the formation of characteristic apoptotic features in the cells with the aid of different nuclear affinity dyes. Henceforth, these are the elemental line of apoptosis identification in diverse cancer cells.
Phase-contrast microscopy
The changes in the morphological features of A549 cells were observed after the treatment with compound 12c at different concentrations along with the untreated control cells. Further, images reported in Fig. 9
A were captured using phase-contrast microscopy after 72 h, reveal the characteristic apoptotic features like changes in morphology (shape, shrinkage) of the cell, reduction in the number of live cells.
Fig. 9
Microscopic observation of the effect of different concentrations (1.25, 2.5, and 5 μM) of compound 12c in comparison to control (untreated) on A549 cells after 72 h of treatment. The red scale bar denotes 50 μm in all the images. (A) Morphological changes were observed through phase-contrast microscopy. (B) Morphological changes such as membrane blebs were visualized through a fluorescent microscope by performing acridine orange staining. (C) Nuclear changes were observed by DAPI staining using a fluorescent microscope. Insets at different concentrations of treatment indicate the changes induced by treatment and red-colored arrows specify the area of effect.
Microscopic observation of the effect of different concentrations (1.25, 2.5, and 5 μM) of compound 12c in comparison to control (untreated) on A549 cells after 72 h of treatment. The red scale bar denotes 50 μm in all the images. (A) Morphological changes were observed through phase-contrast microscopy. (B) Morphological changes such as membrane blebs were visualized through a fluorescent microscope by performing acridine orange staining. (C) Nuclear changes were observed by DAPI staining using a fluorescent microscope. Insets at different concentrations of treatment indicate the changes induced by treatment and red-colored arrows specify the area of effect.
Acridine orange staining
Acridine orange is selective toward nucleic acid and penetrates live cells and stains the nuclei green [40]. In the present study, A549 cells after treatment with compound 12c for 72 h exhibited the formation of apoptotic features such as the appearance of membrane blebs and inverse proportion in the number of cells with a concentration of compound tested. In addition, the cationic dye has the potential to bind lysosomes and thereby also exhibiting the vacuoles as represented in Fig. 9
B.
DAPI staining
The initial discovery of 4,6-diamidino-2-phenylindole (DAPI) in the treatment of protozoal diseases retarded due to a strong binding affinity towards DNA base pairs. Later, this ability of passage through cell membranes made it an efficient stain to detect nuclear changes in the apoptosis cells [41]. The DAPI staining was executed according to the reported procedure [40] to detect apoptosis induction. Compound 12c treated A549 cells after 72 h, on staining with DAPI visualized the chromatin condensation, pyknotic (inset of 1.25 μM), and condensed (bright colored: inset of 2.5 μM) nuclei formation as depicted in Fig. 9
C.
Flow cytometry analysis
Quantification of apoptosis by dual staining assay
To quantify apoptotic cell death and determine the phase of apoptosis, dual staining assay with annexin V-FITC/PI was employed by using flow cytometry. The effect of different concentrations of compound 12c on A549 cells was studied using the principle of phosphatidylserine externalization [42], (a), (b) and compared with control (untreated) cells. The percentage of apoptotic cells at different phases are represented on a four quadrant graph. From Fig. 10
, it is clear that induction of slight early (+annexin V/– PI) and significant late apoptosis (+annexin V/+ PI) after 72 h of incubation with compound 12c.
Fig. 10
Annexin V-FITC/propodium staining dual staining assay on A549 cells to determine the apoptosis induced at different concentrations (1.25, 2.5, and 5 μM) of 12c and compared with untreated (control) cells after 72 h. (A) Flow cytometric analysis of each sample was analyzed and the percentage of Annexin V-FITC/PI positive cells was represented in four quadrants. (B) Quantification of cells in each quadrant (%) with respect to the concentration.
Annexin V-FITC/propodium staining dual staining assay on A549 cells to determine the apoptosis induced at different concentrations (1.25, 2.5, and 5 μM) of 12c and compared with untreated (control) cells after 72 h. (A) Flow cytometric analysis of each sample was analyzed and the percentage of Annexin V-FITC/PI positive cells was represented in four quadrants. (B) Quantification of cells in each quadrant (%) with respect to the concentration.
Depolarization of mitochondrial membrane potential (Δψm)
Mitochondria are the source of energy and also the crucial targets for oxidative stress in a cell [43a-b]. JC-1 staining was performed by standard protocol [43], (a), (b), (c) to assess the impact of 12c on the mitochondrial membrane potential of A549 lung cancer cells. After the treatment period of 72 h, JC-1 dye was added and normal polarized healthy mitochondria labeled as J aggregates (red color) in the control remained intact, whereas 12c treated resulted in the formation of depolarized mitochondria categorized as J monomers. However, the concentration of compound 12c and the formation of J monomers increased proportionately (Fig. 11
).
Fig. 11
Effect of different concentrations (1.25, 2.5, and 5 μM) of compound 12c on the mitochondrial membrane potential of A549 cells after the treatment for 72 h relative to control (untreated). (A) Flow cytometric analysis determines the percentage of J aggregates and J monomers (B) Bar graph represents the quantification of J aggregates and monomers at different concentrations represented in red and blue color respectively.
Effect of different concentrations (1.25, 2.5, and 5 μM) of compound 12c on the mitochondrial membrane potential of A549 cells after the treatment for 72 h relative to control (untreated). (A) Flow cytometric analysis determines the percentage of J aggregates and J monomers (B) Bar graph represents the quantification of J aggregates and monomers at different concentrations represented in red and blue color respectively.
Cell cycle analysis
Cell cycle analysis was conducted on A549 cells using different concentrations (1.25, 2.5, 5 μM) of compound 12c to figure out the phase where cells are arrested. Further, after the treatment for 48 h, cells were interpreted through the flow analyzer after staining with propidium iodide (PI). The number of cells viz., 55.68 in the G0/G1 phase in control increased exponentially to 73.27 in compound 12c (5 μM) treated cells. It is apparent from Fig. 12
A & B that the treatment arrested the G0/G1 phase of the cell cycle in a concentration-dependent manner.
Fig. 12
Effect of varying concentrations (1.25, 2.5, and 5 μM) of potent compound 12c on cell cycle progression in A549 cells for 48 h and compared with control (untreated). The propidium staining method was used to analyze cell cycle distribution. (A) Flow cytometric analysis of cell cycle distribution. (B) Quantification of cells at different phases of the cell cycle.
Effect of varying concentrations (1.25, 2.5, and 5 μM) of potent compound 12c on cell cycle progression in A549 cells for 48 h and compared with control (untreated). The propidium staining method was used to analyze cell cycle distribution. (A) Flow cytometric analysis of cell cycle distribution. (B) Quantification of cells at different phases of the cell cycle.
Metastatic assay
Cancer cells detach from the site of origin, migrate, and develop the disease at the secondary site is known as metastasis, which makes the disease perilous and abating the survival rate of the patient [4], [44].
Cell migration assay
To understand the capability of 12c on cell migration, wound healing assay has been performed on A549 cells. This assay is based on the ability of the healing of an artificial wound called “scratch” created using a sterile pipette tip (200 μL).[45], (a), (b) The cells were treated with a range of concentrations varying from 1.25 to 5 μM of 12c, along with control, and images were captured at 0 and 72 h and compared. It is inferred from Fig. 13
, that the wound in control was practically healed, whereas, the cell migration is inhibited in 12c treated cells after the period of treatment.
Fig. 13
Effect of compound 12c on A549 cell migration. Cells were treated with various concentrations (1.25, 2.5, and 5 μM) of 12c, and scratch was made using a sterile pipette tip and compared with control (untreated cells). Images were captured using a phase-contrast microscope at different periods. (A) Images captured at 0 h (B) Images captured at 72 h.
Effect of compound 12c on A549 cell migration. Cells were treated with various concentrations (1.25, 2.5, and 5 μM) of 12c, and scratch was made using a sterile pipette tip and compared with control (untreated cells). Images were captured using a phase-contrast microscope at different periods. (A) Images captured at 0 h (B) Images captured at 72 h.
Computational analysis
Molecular docking
Further affirmation of intercalative topo-IIα inhibition by synthesized derivatives, molecular docking [46] was performed against active site of human topoisomerase II alpha complexed with DNA [47] using Glide (Schrödinger 2017–1) [48] and the best-scored compounds were detailed along with their interactive pattern. Structural accommodation of the ligand 12c and 12a revealed that β-carboline moiety exhibited intercalation into the DNA in the active site whereas, indole chain extended into the hydrophobic region via thiadiazolo carboxamide as illustrated in A, C of Fig. 14
.
Fig. 14
Molecular docking studies. The ligands were represented in ball and stick model, carbons in ligand are colored pink and DNA in light red and cyan ribbon. (A) Compound 12c in the active pocket of 5GWK. (B) Intercalative binding of compound 12c along with residual interactions. (C) Compound 12a in the active site of 5GWK. (D) Key interactions of compound 12a with amino acids and DNA base pairs in the active site.
Molecular docking studies. The ligands were represented in ball and stick model, carbons in ligand are colored pink and DNA in light red and cyan ribbon. (A) Compound 12c in the active pocket of 5GWK. (B) Intercalative binding of compound 12c along with residual interactions. (C) Compound 12a in the active site of 5GWK. (D) Key interactions of compound 12a with amino acids and DNA base pairs in the active site.The analysis of docking interactions indicated that all the structural fragments (β-carboline, thiadiazole, and indole) of ligand 12c established key interactions with residues in the active pocket such as the formation of two hydrogen bonds between amide attached thiadiazole’s nitrogen and indolic nitrogen with dC8 (d = 3.65 Å) and dT9 (d = 3.51 Å) residues respectively as represented in B of Fig. 14. Further, π-cation interaction was exhibited by tricyclic pyridine NH with dA12 residue and π- π interaction was observed with phenyl ring of β-carboline with dA12 and dG13 residues. Notably, several hydrophobic interactions were observed with different amino acid residues such as Leu486, Met762, Met766, Ile769, Ala801, and Pro803. Additionally, compound 12a exhibited another crucial interaction on the β-carboline viz hydrogen bond between indolic N—H and Arg487 (d = 2.87 Å), along with other interactions as displayed in D of Fig. 14. The docking results reveal the topo II inhibition is a collective contribution by all the structural fragments (Fig. S3
, ESI) of the compounds and are compatible with the in vitro biological evaluation and encourage the structural diversification on the ligands in the near future.
Binding energy and Lipinski’s rules
Binding energies for the molecules were calculated to identify the affinity and complex forming energy of the ligands in the active sites of the human topo-IIα complex with DNA and topo I using the prime MM/GBSA method. The lead molecules exhibited superlative binding energy values of −58.91(12c) and −53.67 (12a) kcal/mol in the active domain of 5GWK (topoisomerase IIα) in comparison to binding energy values of −46.87 (12c) and −40.52 (12a) in the active pocket of 1SC7 (topoisomerase I) (Table S1, ESI). The binding energy values of harmine and co-crystal (etoposide) are −40.94 and −72.21 kcal/mol respectively.The drug likeliness of the potent molecules and standards (Harmine, Etoposide) was analyzed by the evaluation of the Lipinski rule of five and is incorporated in Table 2
. Although the compounds other than 12a violate the Lipinski rule in terms of molecular weight, the optimal range [49] obtained by evaluating the 95% of known drugs recommends the acceptable range up to 725 Daltons reveals the compounds tested satisfies the entire parameters and possess good drug-like properties [50], [51].
Table 2
Lipinski rule of five for potent hybrids and standards.
Entry
R1
R2
Molecular weight
H-bond donors
H-bond acceptors
QPlog P (o/w)a
12a
phenyl
H
486.55
3
5.5
5.63
12c
4-methoxy phenyl
H
516.58
3
6.25
5.717
12f
3,4,5-trimethoxyphenyl
Me
590.66
2
7.75
6.453
12h
4-fluorophenyl
H
504.54
3
5.5
5.863
12j
4-chlorophenyl
H
521
3
5.5
6.12
Harmine
–
–
212.25
1
1.75
3.05
VP-16b
–
–
588.56
3
16.95
0.60
Range or recommended values [49]
130 – 725
0 – 6.0
2 – 20
−2.0 – 6.5
Predicted octanol/water partition co-efficient log P.
VP-16 = Etoposide.
Lipinski rule of five for potent hybrids and standards.Predicted octanol/water partition co-efficient log P.VP-16 = Etoposide.
In silico ADME/T profile
The prediction of pharmacokinetic parameters provides an idea regarding the acceptable range for human use which was obtained by analyzing 95% of known drugs [51], [52]. The Physico-chemical parameters of the five best performers in the cytotoxicity table and standards (Harmine, Etoposide) were predicted by using the QikProp module [53] of Schrödinger and some key pharmacokinetic descriptors along with the recommended range are tabulated. The evaluation of computed parameters such as QPlogKhsa predicts the binding of a ligand to serum albumin which determines the availability of free drug and assessment of QPPCaco, where Caco cells act as a model for gut-blood barrier and QPPMDCK, where MDCK (Madin-Darby Canine Kidney) cells are a mimic for blood–brain barrier assist in the estimation of bioavailability. From the in silico data delineated in Table 3
, all the newly synthesized compounds assessed through QikProp fall in the optimal range and exhibit pharmaceutical significance.
Table 3
Physico-chemical parameters of some of the potent thiadiazolo-carboxamide bridged β-carboline-indole hybrids and standards predicted by QikProp [53]
Descriptors
Ligand entry
Standards
Recommended range
12a
12c
12f
12h
12j
Harmine
VP-16
EAa (ev)
1.24
1.25
1.25
1.28
1.27
0.48
0.22
−0.9 to 1.7
PSAb
101.26
106.32
108.42
101.26
101.27
33.4
166.59
70–200
SAamideOc
0
0
0
0
0
0
0
0–35
QPlogKhsad
1.17
1.16
1.43
1.21
1.29
0.203
−0.65
−1.5 to 1.5
QPpolrze
56.17
57.47
63.73
56.45
57.48
24.38
52.15
13–70
CNSf
−2
−2
−1
−2
−2
1
−2
−2 to +2
QPPCacog
566.66
562.31
1229.128
571.09
571.56
4338.87
212.72
<25 poor, >500 great
QPlogBBh
−1.185
−1.146
−0.881
−1.082
−1.039
0.18
−1.59
−3 to 1.2
QPPMDCKi
430.3
431.40
996.531
632.96
864.13
2417.01
92.85
<25 poor, >500 great
QPlogKpj
−1.534
−1.483
−1.006
−1.667
−1.701
−1.34
−3.76
−8 to −1
IPk (ev)
8.433
8.437
8.383
8.458
8.45
8.35
8.68
7.9 to 10.5
Percent human oral absorption
94.67
83.724
95.867
83.078
84.581
100
46.24
>80% high, <25% poor
PM3 calculated electron affinity.
Van der Waals surface area of polar nitrogen and oxygen atoms.
Solvent-accessible surface area of amide oxygen atoms.
Prediction of binding to human serum albumin.
Predicted polarizability in cubic angstorms.
Predicted central nervous system activity.
Predicted apparent Caco-2 cell permeability in nm/sec. Caco-2 cells are a model for the gut-blood barrier.
Predicted brain/blood partition coefficient.
Predicted apparent MDCK cell permeability in nm/sec. MDCk cells are considered to be a good mimic for the blood–brain barrier.
Predicted skin permeability, log Kp.
PM3 calculated ionization potential.
Physico-chemical parameters of some of the potent thiadiazolo-carboxamide bridged β-carboline-indole hybrids and standards predicted by QikProp [53]PM3 calculated electron affinity.Van der Waals surface area of polar nitrogen and oxygen atoms.Solvent-accessible surface area of amide oxygen atoms.Prediction of binding to human serum albumin.Predicted polarizability in cubic angstorms.Predicted central nervous system activity.Predicted apparent Caco-2 cell permeability in nm/sec. Caco-2 cells are a model for the gut-blood barrier.Predicted brain/blood partition coefficient.Predicted apparent MDCK cell permeability in nm/sec. MDCk cells are considered to be a good mimic for the blood–brain barrier.Predicted skin permeability, log Kp.PM3 calculated ionization potential.
In silico toxicity profiling
In order to understand the toxicity profile of the synthesized compounds and harmine, TOPKAT (TOxicity Prediction by Komputer Assisted Technology) tool from Discovery Studio v19.1.0.18287 was utilized. TOPKAT employs robust and cross-validated Quantitative Structure Toxicity Relationship (QSTR) models to predict toxicological endpoints of given molecules chemical structure. Toxicity model development involves the use of 2D molecular, spatial and electronic descriptors. The optimal Predictive Space validation method of TOPKAT gives assurance to the prediction [54], [55].In silico toxicity profiling results for some of the synthesized compounds and harmine was depicted in Table 4
. The toxicity profile assessment has been predicted using different animal models. The toxicity of synthesized compounds was assessed and compared with harmine in different animal models. Broadly, the probability values of the synthesized compounds showed a lower toxicity profile than harmine. All the analyzed compounds including harmine have shown higher toxicity probabilities values for Ames mutagenicity. Compounds 12a, 12c, 12f, 12h and 12j did not show carcinogenicity call in National Toxicology Program (NTP) carcinogenicity (male and female mouse) model and found to be safe, while harmine has shown positive carcinogenicity signs in the NTP model. All the analyzed molecules had shown a positive sign for carcinogenicity in NTP (male rat) and weight of evidence rodent carcinogenicity (except 12j) model. Except for 12f and harmine, all the molecules had shown a negative sign for carcinogenicity in NTP (female rat) model. Moreover, all the compounds are non-degradable and negative for skin sensitization, skin irritancy and developmental toxicity potential. These compounds had shown higher probability for various animal models which include Rat Male FDA None vs Carcinogen, Rat Male FDA Single vs Multiple, Rat Female FDA None vs Carcinogen and Rat Female FDA Single vs Multiple. However, in silico toxicity studies indicated that the synthesized compounds were found to be safer in many models with a minimal toxicological profile.
Table 4
Probability values of different toxicity models for lead compounds and harmine by TOPKAT analysis.
TOPKAT MODEL
12a
12c
12f
12h
12j
Harmine
Aerobic_Biodegradability
−12.1 (Non-Degradable)
−10.9 (Non-Degradable)
−8.81 (Non-Degradable)
−13.5 (Non-Degradable)
−20.6 (Non-Degradable)
−4.31 (Non-Degradable)
Ames_Mutagenicity
8.02 (Mutagen)
6.11 (Mutagen)
5.98 (Mutagen)
5.99 (Mutagen)
5.7 (Mutagen)
7.15 (Mutagen)
Developmental_Toxicity_Potential
−3.72 (Non-Toxic)
−0.859 (Non-Toxic)
−1.13 (Non-Toxic)
−1.35 (Non-Toxic)
−2.51 (Non-Toxic)
−0.421 (Non-Toxic)
Mouse_Female_FDA_None_vs_Carcinogen
−0.797 (Carcinogen)
−4.72 (Non-Carcinogen)
−5.39 (Non-Carcinogen)
−1.46 (Non-Carcinogen)
−3.87 (Non-Carcinogen)
1.83 (Carcinogen)
Mouse_Female_FDA_Single_vs_Multiple
−4.48 (Single-Carcinogen)
–
–
–
–
0.879 (Multiple-Carcinogen)
Mouse_Female_NTP
−0.856 (Non-Carcinogen)
−3.07 (Non-Carcinogen)
0.818 (Non-Carcinogen)
−1.43 (Non-Carcinogen)
−2.11 (Non-Carcinogen)
2.58 (Carcinogen)
Mouse_Male_FDA_None_vs_Carcinogen
1.49 (Carcinogen)
−0.662 (Carcinogen)
−2.13 (Non-Carcinogen)
−1.72 (Carcinogen)
−2.31 (Non-Carcinogen)
2.8 (Carcinogen)
Mouse_Male_FDA_Single_vs_Multiple
−10.8 (Single-Carcinogen)
−10 (Single-Carcinogen)
–
−13.7 (Single-Carcinogen)
–
−4.03 (Multiple-Carcinogen)
Mouse_Male_NTP
−5.02 (Non-Carcinogen)
−5.35 (Non-Carcinogen)
−4.97 (Non-Carcinogen)
−5.26 (Non-Carcinogen)
−3.75 (Non-Carcinogen)
1.84 (Carcinogen)
Ocular_Irritancy_Mild_vs_Moderate_Severe
−3.2 (Mild)
−4.54 (Mild)
−4.62 (Mild)
−2.61 (Mild)
−2.61 (Mild)
1.22 (Moderate_Severe)
Rat_Female_FDA_None_vs_Carcinogen
−1.2 (Carcinogen)
3.25 (Carcinogen)
3.49 (Carcinogen)
−0.979 (Carcinogen)
−2.83 (Non-Carcinogen)
4.66 (Carcinogen)
Rat_Female_FDA_Single_vs_Multiple
−0.0137 (Single-Carcinogen)
1.12 (Multiple-Carcinogen)
−1.49 (Single-Carcinogen)
1.99 (Multiple-Carcinogen)
–
−0.139 (Single-Carcinogen)
Rat_Female_NTP
−1.36 (Non-Carcinogen)
0.0595 (Non-Carcinogen)
2.93 (Carcinogen)
−2.7 (Non-Carcinogen)
−2.7 (Non-Carcinogen)
1.48 (Carcinogen)
Rat_Male_FDA_None_vs_Carcinogen
3.45 (Carcinogen)
4.55 (Carcinogen)
7.61 (Carcinogen)
5.36 (Carcinogen)
3.06 (Carcinogen)
4.21 (Carcinogen)
Rat_Male_FDA_Single_vs_Multiple
−5.08 (Single-Carcinogen)
−4.07 (Single-Carcinogen)
−4.74 (Single-Carcinogen)
−5.85 (Single-Carcinogen)
−6.54 (Single-Carcinogen)
−1.15 (Single-Carcinogen)
Rat_Male_NTP
3.21 (Carcinogen)
2.1 (Carcinogen)
7.24 (Carcinogen)
3.3 (Carcinogen)
1.25 (Carcinogen)
0.968 (Carcinogen)
Skin_Irritancy_None_vs_Irritant
−1.82 (Non-Irritant)
−1.9 (Non-Irritant)
−1.56 (Non-Irritant)
−2.28 (Non-Irritant)
−2.28 (Non-Irritant)
−4.66 (Mild)
Skin_Sensitization_None_vs_Sensitizer
−4.14 (Non-Sensitizer)
−3.54 (Non-Sensitizer)
−5.06 (Non-Sensitizer)
−3.53 (Non-Sensitizer)
−3.53 (Non-Sensitizer)
−2.95 (Non-Sensitizer)
Weight_of_Evidence_Rodent_Carcinogenicity
1.46 (Carcinogen)
0.889 (Carcinogen)
1.21 (Carcinogen)
2.05 (Carcinogen)
−0.0222 (Non-Carcinogen)
1.08 (Carcinogen)
Carcinogenic_Potency_TD50_Mouse
0.506 mg/kg
0.622 mg/kg
1.08 mg/kg
0.477 mg/kg
0.364 mg/kg
32 mg/kg
Carcinogenic_Potency_TD50_Rat
0.228 mg/kg
0.0373 mg/kg
0.0266 mg/kg
0.0649 mg/kg
0.0601 mg/kg
3.01 mg/kg
Chronic_LOAEL
0.036 g/kg
0.0211 g/kg
0.0177 g/kg
0.0288 g/kg
0.0264 g/kg
0.0139 g/kg
Daphnia_EC50
0.0539 mg/l
0.0257 mg/l
0.0382 mg/l
0.053 mg/l
0.0375 mg/l
1.6 mg/l
Fathead_Minnow_LC50
2.07e-05 g/l
9.2e-06 g/l
2.19e-06 g/l
2.21e-05 g/l
1.12e-05 g/l
0.0392 g/l
Rat_Inhalational_LC50
7.65e + 03 mg/m3/h
1.03e + 04 mg/m3/h
6.33e + 03 mg/m3/h
1.25e + 04 mg/m3/h
5.14e + 03 mg/m3/h
2.43e + 04 mg/m3/h
Rat_Maximum_Tolerated_Dose_Feed
0.165 g/kg
0.0771 g/kg
0.039 g/kg
0.216 g/kg
0.199 g/kg
0.0441 g/kg
Rat_Maximum_Tolerated_Dose_Gavage
0.00422 g/kg
0.000328 g/kg
4.86e-05 g/kg
0.000435 g/kg
0.000336 g/kg
0.00477 g/kg
Rat_Oral_LD50
2.89 g/kg
2.69 g/kg
1.88 g/kg
1.75 g/kg
3.37 g/kg
0.202 g/kg
Probability values of different toxicity models for lead compounds and harmine by TOPKAT analysis.
Conclusion
In summary, new molecules 12a–k were rationally designed by taking into consideration of structural fragments and pharmacophoric hybridization approach to build intercalative topoisomerase IIα inhibitors targeting cancer. Based on the earlier presented pharmacophoric profile of these heterocycles (β-carboline, indolo-thiadiazole), the synthesis was performed and adjoined by carboxamide formation which resulted in the target molecules with good yields. From the preliminary in vitro screening, it can be observed that the compounds 12c and 12a remarkably expressed its cytotoxicity among all the derivatives tested in all the cell lines (solid and liquid tumors) with magnificent IC50 values of 2.82 ± 0.10 and 3.00 ± 1.40 μM, respectively in A549 cell line. Also, cytospecificity and selectivity index of potent compound 12c towards cancer cells was revealed by calculating cytotoxicity (IC50: 71.25 ± 7.95) against normal human lung epithelial cell line (BEAS-2B). Catalytic inhibition of topo-IIα by compounds 12c and 12a was evinced through gel-based electrophoresis and intercalation into DNA was revealed by evaluating alterations in parameters such as absorbance, viscosity, and correlation with standard intercalators. Additionally, deeper insights into mechanistic features by different flow cytometric assays (Annexin V/FITC, JC-1) indicated that compound 12c is a potent inducer of apoptosis. Cell cycle progression in A549 cells was arrested at the G0/G1 phase by compound 12c in a dose-dependent manner. Furthermore, cellular and nucleo-morphological changes caused by 12c were also visualized using diverse staining assays (Acridine-orange and DAPI), revealed the induction of abundant apoptotic features in A549 cells. Scratch wound assay by compound 12c indicated the inhibition of in vitro cell migration and deterrence of wound closure in A549 cells. Further, computational studies supported the in vitro target-based assays as compounds 12c and 12a accommodated well in the catalytic region of topo-IIα with residual amino acid and DNA base pairs interactions (PDB code: 5GWK). Further, in silico evaluation of ADME/T profile and drug-likeliness indicated that all the molecules fall under the specified range deduced from known drugs. As a concluding remark, the current study primarily emphasizes on assembling varied pharmacological aspects of designed compounds and can be suitably derivatized to mark their progression as promising anticancer agents.
Experimental section
Chemistry
Materials and methods: All the reagents and solvents were obtained from commercial suppliers and were used without further purification. Analytical thin-layer chromatography (TLC) was performed on MERCK precoated silica gel 60-F254 (0.5 mm) aluminium plates. Visualization of the spots on TLC plates was achieved by UV light. Wherever required, column chromatography was performed using silica gel (60–120). The reactions wherever anhydrous conditions are required were carried under nitrogen positive pressure using freshly distilled solvents. All evaporation of solvents was carried out under reduced pressure using a rotary evaporator below 45 °C. 1H and 13C NMR spectra were recorded on Bruker 500 MHz by making a solution of samples in the DMSO using tetramethylsilane (TMS) as the internal standard. Chemical shifts for 1H and 13C are reported in parts per million (ppm) downfield from tetramethylsilane. Spin multiplicities are described as s (singlet), d (doublet), t (triplet), q (quartet), and m (multiplet). Coupling constant (J) values are reported in hertz (Hz). HRMS was determined with Agilent QTOF mass spectrometer 6540 series instrument. Melting points were determined with an electrothermal digital melting point apparatus IA9100 and are uncorrected. The names of all the compounds given in the experimental section were taken from Chem Bio Draw Ultra, Version 12.0.
General procedure for N-(5-(1-alkyl-1H-indol-2-yl)-1,3,4-thiadiazol-2-yl)-1-aryl-9H-pyrido[3,4-b]indole-3-carboxamide (12a–k)
To a mixture of 1-aryl-9H-pyrido[3,4-b]indole-3-carboxylic acid (6a–f, 1.1 equiv.), EDC (1.5 equiv.), HOBt (1.5 equiv.) under nitrogen conditions, dichloromethane was added and allowed to stir at room temperature for 15 min. Later, the amino component (5-(1-alkyl-indol-2-yl)-1,3,4-thiadiazol-2-amine 11a–c, 1 equiv.) dissolved in dichloromethane, was added and the reaction was allowed to stir at room temperature till complete consumption of the starting materials as determined by TLC. The reaction mixture was extracted with dichloromethane and water (3 × 50 mL) and purified by column chromatography to yield N-(5-(1-alkyl-1H-indol-2-yl)-1,3,4-thiadiazol-2-yl)-1-phenyl-9H-pyrido[3,4-b]indole-3-carboxamide 12a–k in moderate to good yields (72–85%).
Normal human bronchial epithelial cell line (BEAS-2B) and cancer cell lines like lung (A549), breast (MDA-MB-231, BT-474), colon (HCT-116), Leukemia (HL-60, THP-1), were maintained in appropriate media supplemented with 10% fetal bovine serum (FBS) stabilized with 1% antibiotic–antimycotic solution (Sigma Aldrich) and cells were maintained at 37 °C with 5% CO2 and 98% relative humidity in the incubator. When the cells reached up to 80–90% of confluency, sub-culturing was performed using a 0.25% trypsin/1 mM EDTA solution for further passage. The stock solution of 10 mM was prepared by dissolving test compounds and standards in DMSO. Required concentrations were further obtained by suitable dilutions with respective media.
MTT assay
MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) is a dye that converts into insoluble formazan by mitochondrial succinate dehydrogenase enzyme, where MTT assay is a colorimetric assay, thereby measures the reduction of MTT. Metabolically active cells can reduce the MTT and the level of activity indicates the viability of cells. Briefly, cells were seeded in 96-well plates at a density of 1000 to 4000 cells per well in 100 µL of complete medium and allowed to grow overnight for attachment onto the wells. After the media was replaced with fresh media, the cells were treated with various concentrations of the compounds for a period of 72 h. After incubation, the media was aspirated and 100 µL of MTT (0.5 mg/mL) was added and incubated at 37 °C for 4 h. Then MTT reagent was aspirated and the formazan crystals formed were dissolved by the addition of 200 µL of DMSO for 20 min. at 37 °C. The formazan product quantity was measured by using a spectrophotometric microtiter plate reader (Spectra Max, M4 Molecular Devices, USA) at 570 nm wavelength.
Target-based assays
Topoisomerase IIα inhibition
Topo IIα inhibition of potent compounds 12c and 12a at a concentration of 10 μM was evaluated by performing ATP dependent decatenation of kDNA by the protocol described in the Drug Screening Kit of Topo II (TG1019-1, TopoGEN, USA). Reactions were carried out in 20 μL and contained two different types of buffer in equal ratio and 200–300 ng of kDNA, and topoisomerase II (TG2000H-2, TopoGEN, USA). The quantity of topoisomerase II (5 units) was attuned in pilot experiments to decatenate around 100% of the kDNA under our assay conditions. Authentic decatenated DNA and linear DNA was used as a control to identify the changes. The samples were made with buffer, kDNA, test compounds, topo-IIα and incubated at 37 °C for 30 min. and terminated by the addition of 2 μL of a stop buffer containing 10% (w/v) SDS and 2 μL of 0.5 mg/mL proteinase-K and incubated for 10 min at 37 °C. After completion of the reaction, the products in the reaction mixture were separated by 1% agarose gel. The products in the agarose gel were visualized after staining with ethidium bromide (0.2 mg/mL). The gels were run at 100 V for about 30 min. and visualized under UV transillumination (BIO-RAD gel doc XR+, USA).
DNA intercalation assays
Relative viscosity studies
The viscosities of the DNA-ligand complexes were determined by the Lovis 2000 M/ME Rolling-ball viscometer (Anton Paar GmbH, Graz, Austria), based on the falling ball principle. The temperature was controlled at ±0.005 K through an internal Peltier thermostat. A calibrated 1.59 mm glass capillary containing a steel ball was filled with the sample for measuring the ball falling time at angles in the range from 20° to 70°. The kinematic, as well as dynamic viscosities at 25 °C, were estimated based on the ball falling time and densities DNA solution was prepared in 100 mM Tris-HCl (pH 7.4) and viscosity was measured while each derivative 12c and 12a (5 µM) was added to CT-DNA solution (50 µM). Ethidium bromide, Harmine and Hoechst 33258 at different concentrations were used as controls. Data was represented graphically as (ŋ/ŋ0)1/3 vs. the ratio of the concentration of the hybrid to CT-DNA, where ŋ is the viscosity of CT-DNA in the presence of the various concentrations of derivatives and ŋ0 is the viscosity of CT-DNA solution.
Absorbance spectroscopy
DNA nanodrop method
DNA intercalation is determined by nanodrop spectrophotometric analysis. The intercalating agents decrease the absorption and increase the wavelengths. In our experiment, we incubated the 50 μM of calf thymus DNA (Sigma-Aldrich, USA) with compound 12c and 12a, Ethidium bromide (EtBr), Doxorubicin (DOX) and Harmine for 10 min at 1 μM concentration. Later, the TAE buffer used as blank and absorbance and concentrations were determined by NanoDrop™ 2000/2000c Spectrophotometer (Thermo Fisher Scientific, USA).
UV–Visible spectroscopy
DNA intercalation is determined by UV–Visible spectroscopic analysis. In our experiment, we incubated the 50 μM of calf thymus DNA (Sigma-Aldrich, USA) with compound 12c and 12a, and Ethidium bromide (EtBr), for 10 min at 1 μM concentration and compared with control (i.e. CT-DNA alone). Later, the TAE buffer used as blank and absorbance and concentrations were determined by UV–Visible Spectrophotometer (Jasco; V-650).
Apoptosis detection studies
Morphological observations
Phase-contrast microscopy
A549 cells were plated in 12 well culture plates with a cell density of 1 × 106 cells/mL and allowed to adhere overnight. The cells were incubated with different (1.25, 2.5, and 5 µM) concentrations of the compound 12c for 72 h treatment. Later cells were checked for the morphological changes, and images were taken using a phase-contrast microscope (Nikon, Inc. Japan).
Acridine orange (AO) staining
A549 cells were plated at a concentration of 1 × 106 cells/mL and treated with various concentrations of compound 12c and the plates were incubated for 72 h. 10 µL of fluorescent dye, Acridine Orange (AO) was added into each well in equal volumes (10 µg/ml) respectively then the cells were visualized immediately under a fluorescence microscope (Nikon, Inc. Japan) with excitation (488 nm) and emission (550 nm) at 200X magnification.
DAPI nucleic acid staining
DAPI staining aids in the observation of morphological changes in the nucleus. After treatment with compound 12c for 72 h, lung cancer A549 cells were washed with PBS and permeabilized with 0.1% Triton X for 10 min followed by staining with 1 µM DAPI. Control and treated cells were observed with a fluorescence microscope with excitation at 359 nm and emission at 461 nm using DAPI filter at 200X magnification.
Flow cytometric analysis
Annexin V assay
This assay was performed by the method [56] given by Rieger et al. with slight modifications. Briefly, 1 × 105 cells were seeded in a 12-well plate and treated with different concentrations of compound 12c for 72 h. The collected cells were washed twice with ice-cold PBS, then incubated with 200 μL of 1× binding buffer containing 1 μL propidium iodide (PI) for 15 min. at room temperature in the dark. After incubation, cells were analyzed for apoptosis using a flow cytometer (BD FACSVerseTM, USA). Apoptosis and necrosis were analyzed with quadrant statistics on propidium iodide-negative cells, fluorescein positive cells, and propidium iodide (PI)-positive cells, respectively.
Measurement of mitochondrial membrane potential
A549 cells (1 × 106 cells/ml) were seeded in 12 well plates and allowed to adhere overnight. The cells were incubated with compound 12c at 1.25, 2.5, and 5 µM concentrations for 72 h. Cells were collected and washed with PBS and resuspended in a solution of JC-1 (1 µM) and incubated for 30 min. in an incubator at 37 °C. The cells were washed twice with PBS and analyzed by flow cytometer (BD FACSVerse™, USA).
Cell cycle analysis
In general, the novel compounds exert their cytotoxic or growth inhibitory effect by arresting the specific checkpoint in the cell cycle. Flow cytometric analysis (FACS) was performed to calculate the distribution of the cell population in various cell cycle phases. Here A549 cancer cells were incubated with compound 12c at various concentrations from 1.25 to 5 µM for 48 h. Untreated and treated cells were harvested, washed, and fixed overnight in 70% ethanol in PBS at −20 °C. Fixed cells were pelleted and stained with cell cycle analysis reagent propidium iodide (50 μg/ml) with RNase A for 20 min at 37 °C in dark according to the manual instructions and about 10,000 events were acquired and analyzed on a flow cytometer BD FACSVerseTM (BD Biosciences, USA).
Metastatic assays
Cell migration assay/wound healing assay
A549 cells were plated at a cell density of 3.5 × 105 cells/well into 12-well plates in RPMI medium and allowed to form a confluent monolayer, and then the scratch is made with 200 μL sterile pipette tip and then treated with 1.25, 2.5 and 5 μM of the compound 12c. The migration of lung cancer cells were captured by microscopic observations at 0 h and 72 h.
Computational studies
The DNA crystal structure has been retrieved from Protein Data Bank (PDB ID: 5GWK) [46] and 1SC7. The protein preparation tool was used for the preparation of the human topoisomerase II alpha with DNA and topoisomerase I. The missing atoms will be added up and peripheral water molecules will be removed with a distance of less than 5 Å from the pocket through the protein preparation wizard. The grid is generated by picking the active site where the co-crystal is located with a grid box of 10 × 10 × 10 Å (Schrödinger 2017–1) [47]. The potent hybrids were sketched (12c & a) by using 2D sketcher and energy minimized, and prepared using Ligprep for the generation of different conformers (Schrödinger 2017–1). The different conformers thus obtained were subjected to molecular docking with SP Glide (Schrödinger 2017–1). The poses generated were evaluated and the best one was described. Five potent compounds and standards were subjected to ligand preparation and the diverse conformers thus generated were envisaged for their Physico-chemical parameters using QikProp [52]. The different descriptors were thus reported.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Authors: X Shui; M E Peek; L A Lipscomb; A P Wilkinson; L D Williams; M Gao; C Ogata; B P Roques; C Garbay-Jaureguiberry; A P Wilkinson; L D Williams Journal: Curr Med Chem Date: 2000-01 Impact factor: 4.530
Authors: Kamal S Saini; Carlo Lanza; Marco Romano; Evandro de Azambuja; Javier Cortes; Begoña de Las Heras; Javier de Castro; Monika Lamba Saini; Sibylle Loibl; Giuseppe Curigliano; Chris Twelves; Manuela Leone; Mrinal M Patnaik Journal: Br J Cancer Date: 2020-06-22 Impact factor: 7.640