Yasmin M Syam1, Manal M Anwar1, Somaia S Abd El-Karim1, Khaled M Elokely2, Sameh H Abdelwahed3. 1. Department of Therapeutic Chemistry, National Research Center, Dokki, Cairo 12622, Egypt. 2. Institute for Computational Molecular Science, Department of Chemistry, Temple University, Philadelphia, PA 19122, USA. 3. Department of Chemistry, Prairie View A&M University, Prairie View, TX 77446, USA.
Abstract
Herein, 2,3-dioxo-1,2,3,4-tetrahydroquinoxaline was used as a bio-isosteric scaffold to the phthalazinone motif of the standard drug Olaparib to design and synthesize new derivatives of potential PARP-1 inhibitory activity using the 6-sulfonohydrazide analog 3 as the key intermediate. Although the new compounds represented the PARP-1 suppression impact of IC50 values in the nanomolar range, compounds 8a, 5 were the most promising suppressors, producing IC50 values of 2.31 and 3.05 nM compared to Olaparib with IC50 of 4.40 nM. Compounds 4, 10b, and 11b showed a mild decrease in the potency of the IC50 range of 6.35-8.73 nM. Furthermore, compounds 4, 5, 8a, 10b, and 11b were evaluated as in vitro antiproliferative agents against the mutant BRCA1 (MDA-MB-436, breast cancer) compared to Olaparib as a positive control. Compound 5 exhibited the most significant potency of IC50; 2.57 µM, whereas the IC50 value of Olaparib was 8.90 µM. In addition, the examined derivatives displayed a promising safety profile against the normal WI-38 cell line. Cell cycle, apoptosis, and autophagy analyses were carried out in the MDA-MB-436 cell line for compound 5, which exhibited cell growth arrest at the G2/M phase, in addition to induction of programmed apoptosis and an increase in the autophagic process. Molecular docking of the compounds 4, 5, 8a, 10b, and 11b into the active site of PARP-1 was carried out to determine their modes of interaction. In addition, an in silico ADMET study was performed. The results evidenced that compound 5 could serve as a new framework for discovering new potent anticancer agents targeting the PARP-1 enzyme.
Herein, 2,3-dioxo-1,2,3,4-tetrahydroquinoxaline was used as a bio-isosteric scaffold to the phthalazinone motif of the standard drug Olaparib to design and synthesize new derivatives of potential PARP-1 inhibitory activity using the 6-sulfonohydrazide analog 3 as the key intermediate. Although the new compounds represented the PARP-1 suppression impact of IC50 values in the nanomolar range, compounds 8a, 5 were the most promising suppressors, producing IC50 values of 2.31 and 3.05 nM compared to Olaparib with IC50 of 4.40 nM. Compounds 4, 10b, and 11b showed a mild decrease in the potency of the IC50 range of 6.35-8.73 nM. Furthermore, compounds 4, 5, 8a, 10b, and 11b were evaluated as in vitro antiproliferative agents against the mutant BRCA1 (MDA-MB-436, breast cancer) compared to Olaparib as a positive control. Compound 5 exhibited the most significant potency of IC50; 2.57 µM, whereas the IC50 value of Olaparib was 8.90 µM. In addition, the examined derivatives displayed a promising safety profile against the normal WI-38 cell line. Cell cycle, apoptosis, and autophagy analyses were carried out in the MDA-MB-436 cell line for compound 5, which exhibited cell growth arrest at the G2/M phase, in addition to induction of programmed apoptosis and an increase in the autophagic process. Molecular docking of the compounds 4, 5, 8a, 10b, and 11b into the active site of PARP-1 was carried out to determine their modes of interaction. In addition, an in silico ADMET study was performed. The results evidenced that compound 5 could serve as a new framework for discovering new potent anticancer agents targeting the PARP-1 enzyme.
Poly (ADP-ribose) polymerases (PARPs) constitute a group of at least 17 enzymes that are correlated to the DNA damage repair process. PARP-1is the most abundant member of this group and has emerged as one of the most auspicious molecular targets for cancer management in the past decade [1,2,3,4]. PARP-1 acts as a “molecular nick sensor” to DNA single-strand (ssDNA) breaks and catalyzes the transference of ADP-ribose units (utilizing nicotinamide adenine dinucleotide (NAD+) as a substrate) to acceptor proteins, facilitating the recruitment of the damaged DNA and promoting cell survival. It is an important stage in the base excision repair (BER) of single-strand DNA breaks [4], which is linked to the resistance that typically develops following traditional cancer treatments [5,6,7]. PARP-1 suppression enhances the damage of injured DNA resulting in synthetic lethality in DNA-repairing-deficient cancer cells, such as BRCA1/2-deficient cells. Thus, PARP-1 suppression synergizes the impact of various antiproliferative drugs such as topoisomerase-I inhibitors and DNA alkylating drugs in addition to ionizing radiation. Moreover, some PARP suppressors are effective as single agents against cancers bearing BRCA1- or BRCA2-mutations [8,9,10,11].The US FDA recently approved four PARP suppressors, Olaparib, Rucaparib, Niraparib, and Talazoparib, for curing BRCA-mutated, HER2-negative advanced, metastatic ovarian, or breast cancer. In addition, there are a number of PARP suppressors that are under study in various clinical phases such as Veliparib, Pamiparib, Simmiparib, and Fluzoparib [12,13] (Figure 1). Furthermore, recent studies investigated the therapeutic potential of various PARP-1 suppressors for other refractory diseases such as Alzheimer’s disease (AD) [14,15]. Accordingly, the development of effective PARP-1 inhibitors plays an important role in medicinal chemistry communities.
Figure 1
Structures of PARP-1 inhibitors approved by FDA and others under clinical studies.
It has been reported that the catalytic pocket of PARP1 is divided into three sub-pockets that are occupied by the substrate NAD+. The first sub-pocket is the nicotinamide-ribose binding site (NI site), the second is the phosphate-binding site (PH site), and the third is the adenine-ribose binding site (AD site) [16] (Figure 2).
Figure 2
The three catalytic sub-pockets of PARP-1enzyme.
It has been reported that most of the PARP-1 suppressors bind with the NI site via H-binding and π-π stacking interactions, and some of them produce further interactions in the adenine-ribose binding (AD) site, which is large enough to fit a variety of molecules, leading to enhancing their effectiveness and pharmacokinetic characteristics [17,18]. Many studies have determined that the design of PARP-1 inhibitors is based on the nicotinamide section of NAD+ to imitate the ligand-protein binding of NAD+ with PARP-1 [19]. Accordingly, PARP-1 suppressors shar common pharmacophoric features, which are an aromatic ring and a carboxamide core. The critical bindings between them are the H-bonding networks initiated between the carboxamide moiety and Gly863 (NH to Gly C=O and C=O to Gly NH) and Ser904 (C=O to Ser OH). Additionally, the phenyl ring of Tyr907 induces the π-π stacking interaction with the aryl ring. Additionally, an auxiliary appendage with a linking side chain is commonly conjugated with the polycyclic core as a solvent accessory region in the AD site [20,21,22,23] (Figure 3).
Figure 3
The design approach of the targeted quinoxaline-based derivatives 3–12 as PARP-1 inhibitors.
Plenty of research has shown that the improvement of PARP inhibitors’ binding affinity by restriction of the carboxamide’s free rotation greatly enhances the PARP1 inhibitory activity. The carboxamide moiety can be locked into the required confirmation by either inserting the aromatic ring heteroatoms or functionalities that can form an intramolecular hydrogen bond with the amide NH or conjugating the amide group in a bicyclic system [4,6,24].Quinoxaline is a privileged scaffold and one of the main blocks of different anticancer agents as it has been proven to be selective adenosine triphosphate (ATP) competitive as well as a bioisostere to benzimidazole, quinazolinones, isoquinolinones, phenanthridone, or phthalazinones, which are the basic scaffolds of the plurality of PARP-1 inhibitors [25,26,27]. In addition, sulfonyl and sulfonamide moieties conjugated to different heterocyclic ring systems have been reported as one of the most privileged scaffolds to inhibit the growth of various human cancer cell lines via different modes of action [28,29]. Cancer treatment is still a challenge due to the development of cancer cell resistance, toxicity, and the lack of selectivity of most commercialized anticancer medications. As a result, and in view of the continuation of our efforts in discovering new heterocyclic compounds of potential anticancer activity targeting the PARP-1 enzyme [30,31,32], the strategy of this study was focused on the design and synthesis of the new compounds based on the quinoxaline core to occupy the NI site of PARP-1 hybridized at its position-6 with different heterocycles, such as pyrrole, pyrazole, thiazole, imidazolidinone, and pyrimidine via sulfonyl, sulfonamide, and sulfonohydrazide linkers aiming to engage with the PARP-1 enzyme through different binding modes of action. The quinoxaline nucleus bears two carboxamide moieties that engage with the enzyme through additional hydrogen bonding (Figure 3). All new compounds were examined as PARP-1 inhibitors. Since PARP-1 inhibitors result in synthetic lethal effects, specifically in BRCA-mutated cells, MDA-MB-436 (BRCA-1-mutated breast cancer cell line) was selected to conduct a cell proliferation assay for the analogs that exhibited the most active inhibition effect against the PARP-1 enzyme. Thereafter, the safety margin of the most potent members was evaluated against WISH normal cells. A molecular docking study was also employed for the promising PARP-1 inhibiting candidates to rationalize and emphasize their mechanisms of binding with the active pocket of the target enzyme. Furthermore, in silico ADMET prediction was performed for the new compounds to explore their drug-likeness characteristics.
2. Results
2.1. Chemistry
This study was directed toward the design and construction of new quinoxaline compounds using various synthetic pathways illustrated in Scheme 1 and Scheme 2. The synthesis was initiated by reacting the starting material o-phenylenediamine with oxalic acid in the presence of HCl to provide quinoxaline-2,3(1H,4H)-dione (1), which was treated with chlorosulfonic acid to provide the corresponding key intermediate 6-sulfonyl chloride derivative 2 according to the reported methodology [33,34]. The latter derivative served as a facile intermediate for the nucleophilic substitution reaction with hydrazine hydrate to afford the 6- sulfonohydrazide derivative 3, which was utilized as a precursor for the ring closure reaction with different active methylene reagents, namely, ethyl-acetoacetate, acetylacetone, and diethyl malonate, to accomplish the corresponding pyrazole derivatives 4–6, respectively.
Scheme 1
Synthesis of different new quinoxaline-2,3-dione -based derivatives.
Scheme 2
Synthesis of different, new 2,3-dioxoquinoxaline-6-sulfonohydrazide-based derivatives.
IR spectra of compounds 3–6 demonstrated stretching bands at approximately 3441–3315 cm−1 corresponding to the NH2 and NH groups, other bands ranging from 1750 to 1678 cm−1 due to C=O groups, and at 1392–1138 corresponding to SO2 groups.1H NMR spectra of compounds 3–6 represented D2O exchangeable signals of NH2 and NH functionalities in the range of δ 11.78–11.83 ppm, alongside multiplet signals at the region of δ 7.01–8.03 ppm related to the aromatic protons. The sulfonohydrazide compound 3 showed an additional D2O exchangeable singlet at δ 7.85 ppm assignable to the NH2 group, while the target pyrazole derivatives 4, 5, and 6 represented new singlets at δ 2.12–2.74 ppm related to CH3, 2CH3, and CH2 functionalities, respectively, at δ 6.78–6.31 ppm due to the pyrazole-H4 of compounds 4 and 5, and at 12.35 ppm exchangeable with D2O related to the OH group of compound 4. Furthermore, 13C NMR spectra of compounds 4, 5, and 6 exhibited singlet signals at δ 12.53, 11.21, and 55.92 ppm assignable to CH3, 2CH3, and CH2 groups, respectively, at δ 102.07–155.73 related to the aromatic carbons, and at δ 154.82–167.89 ppm due to C=O groups.The further condensation reaction of 6- sulfonohydrazide derivative 3 with different acid anhydrides, namely succinic, maleic, and/or phthalic anhydride in glacial acetic acid, resulted in the achievement of the corresponding analogs 7a–c, respectively. 1H NMR spectra of the latter derivatives 7a–c revealed multiplet signals in the range of δ 7.32–7.93 ppm, contributing to the aromatic protons, as well as three D2O exchangeable signals at the region δ 9.61–12.01 ppm due to 3NH groups. Compound 7a represented a singlet signal at δ 2.73 ppm corresponding to the dioxopyrrolidine-2CH2 function, 7b exhibited a doublet signal at δ 7.20 ppm attributed to its vinylic protons, while 7c revealed an increase in the integration values in the aromatic region due to the phthalic protons. Moreover, 13C NMR spectra of compounds 7a–c exhibited singlet signals in the range of δ 115.13–138.42 ppm assigned to the aromatic carbons and in the range of δ 154.34–170.81 ppm representing C=O groups. A singlet signal appeared at δ 30.07 ppm due to the two methylene carbons of the pyrrolidine-2CH2 of compound 7a (Scheme 1).Furthermore, the reaction of 6-sulfonohydrazide derivative 3 with 4-methoxybenzene isothiocyanate and/or benzoyl isothiocyanate in refluxing DMF in the presence of a few drops of triethylamine resulted in the formation of the thiosemicarbazide derivatives 8a,b, respectively. 1H NMR spectra of the latter derivatives 8a,b represented singlet signals at δ 10.71–11.97 ppm exchangeable with D2O affordable to NH groups, multiplet signals at δ 6.90–8.14 ppm related to the aromatic protons, and a singlet signal at δ 3.81 ppm assignable to the methoxy protons (-OCH3) in the case of compound 8a. Furthermore, their 13C NMR spectra revealed the parent carbons of both derivatives in companion with a singlet signal at δ 55.44 ppm related to OCH3 of 8a and at δ 163 ppm due to ph-C=O of the benzoyl compound 8b.Thiosemicarbazide congeners are reported to be valuable intermediates in organic chemistry since they act as building blocks for the preparation of various heterocyclic compounds possessing biological importance [35,36]. Accordingly, the treatment of compounds 8a,b with diethyl malonic acid in refluxing ethanol furnished the corresponding 2-thioxo-3,4-dihydropyrimidine derivatives 9a,b, respectively. IR spectra of compounds 9a,b displayed different absorption bands at 3441–3160, 1710–1645, 1415, 1338, and 1196 cm−1 related to 3NH, 3C=O, C=S, and SO2 groups. 1H NMR spectra of compounds 9a,b represented an additional signal at δ 10.51 ppm exchangeable with D2O referring to the OH group, while the pyrimidine-H5 appeared as a singlet signal in the aromatic region at δ 7.12–7.35 ppm alongside the parent protons, which appeared in their expected regions. Furthermore, 13C NMR spectra of compounds 9a,b represented singlet signals at the region of δ 82.50–155.47 ppm related to the aromatic carbons, at the region of δ 155.40–155.50 due to C=O groups, and at δ 170.1, 189.68 ppm due to C=S groups. The methoxy carbon of compound 9a appeared as a singlet signal at δ 55.80 ppm.Moreover, the nucleophilic reaction of various α- halo ketones, namely chloroacetone and 3-chloroacetylacetone with compounds 8a,b, was carried out in the presence of sodium acetate to give the corresponding thiazolines, 10a,b, and 11a,b, respectively. 1H NMR spectra of compounds 10a,b exhibited a singlet signal in the region of δ 1.85–2.31 ppm contributing to the thiazoline-CH protons, while the thiazoline-H5 appeared as a singlet signal at δ 5.45 ppm, in addition to the precursor protons, which was presented in the correct regions. Furthermore, 1H NMR spectra of 11a,b exhibited two singlet signals at δ 2.33 and 2.51 ppm assigned to thiazoline-CH and COCH, respectively, alongside the signals of the parent protons. Similarly,13C NMR spectra of 10a,b, and 11a,b represented CH3 carbons as singlet signals in the region of δ 14.56–19.77 ppm, in addition to the parent carbons, which were presented in their correct regions. The acetyl carbon of compounds 11a,b was represented as a singlet signal at δ 25.13 and 25.62 ppm, respectively.Moreover, the treatment of 8a,b with ethyl bromoacetate in absolute ethanol containing a catalytic amount of sodium acetate accomplished the corresponding thiazolidine derivatives 12a,b, respectively. 1H NMR spectra of the target derivatives 12a,b showed an up-field signal in the region of δ 4.67 ppm corresponding to the thiazolidine-CH methylene protons alongside the parent protons, which appeared at their expected regions. 13C NMR of the latter derivatives showed a singlet signal at δ 25.13 assignable to the thiazolidine-CH2, and singlet signals in the regions of δ 114.66–155.68 ppm and 163.27–176.33 due to the aromatic and C=O carbons, respectively. The methoxy carbon of 12a appeared as a singlet signal at δ 63.09 ppm (Scheme 2). Mass spectra of the newly prepared compounds exhibited correct molecular ion peaks, which were in accordance with their molecular formulae.
2.2. Biological Evaluation
2.2.1. Assessment of PARP-1 Inhibitory Activity of the Target Quinoxaline Compounds
All the new target quinoxaline compounds 3–12a,b were evaluated as PARP-1 inhibitors to gain clear insight into the structure–activity relationship using the colorimetric 96-well PARP-1 assay kit [6,35]. The IC50 values of all the tested compounds against PARP-1 were expressed in nM concentrations utilizing Olaparib as a reference drug and are summarized in Table 1. Despite all the obtained IC50 values being in the nanomolar range, they showed a wide variation in PARP-1 inhibitory activity (IC50 range; 2.31–57.35 nM). Therefore, it can be supposed that the terminal position of the molecule can tolerate a wide variety of substituents and this point can be explained if these side chains are approaching the solvent surface and do not bind significantly with the enzyme. The key starting intermediate 6-sulfonohydrazide derivative 3 displayed 3-fold less inhibitory activity against PARP-1 than that of the reference drug Olaparib of IC50s; 12.86, 4.40 nM, respectively. We detected the direct conjugation of the parent 6-sulfonoquinoxaline core with a 3,5-dimethylpyrazole ring as compound 5 exhibited PARP-1 inhibitory activity higher than that of the control drug by 1.5 folds (IC50s = 3.05 nM). Conversely, the replacement of either one or both methyl groups of the pyrazole ring of compound 5 with OH or 2C=O groups as compounds 4 and 6, respectively, decreased the suppression impact by nearly 2- and 3-fold of IC50s = 8.73, 13.27 nM, respectively. This result indicated that the hydrophobic residues are favorable for PARP-1 inhibition activity.
Table 1
PARP-1 inhibitory activity and the cytotoxicity of the new quinoxaline-based compounds.
Compound No.
IC50 (mean ± SEM) (nM)
IC50 (mean ± SEM) (µM)
PARP-1
MDA-MB-436
WI-38
3
12.86 ± 4.73
4
8.73 ± 0.44
30.30 ± 1.78
5
3.05 ± 0.16
2.57 ± 0.15 SI = 31.77
81.67 ± 1.70
6
13.27 ± 0.68
7a
57.14 ± 2.91
7b
35.82 ± 1.80
7c
43.40 ± 0.10
8a
2.31 ± 0.30
10.70 ± 0.63 SI = 6.58
70.46 ± 0.43
8b
11.06 ± 0.56
9a
57.35 ± 2.90
9b
35.71 ± 1.82
10a
21.63 ± 1.10
10b
6.35 ± 0.32
9.62 ± 0.56 SI = 7.86
75.66 ± 0.51
11a
19.45 ± 0.99
11b
8.25 ± 0.42
11.50 ± 0.67 SI = 6.96
80.12 ± 0.82
12a
36.11 ± 1.84
12b
40.54 ± 2.06
Olaparib
4.40 ± 0.30
8.63 ± 1.25
Moreover, hybridization of the 6-sulfonoquinoxaline scaffold with the p-methoxyphenyl ring via the thiosemicarbazide linker as compound 8a represented a promising impact on the PARP-1 suppression effect, nearly 2-fold higher than that of Olaparib of IC50 = 2.31 nM. On the other hand, the inhibitory activity was 2.5 times less than the reference drug upon conjugation of the thiosemicarbazide side chain with the CO-unsubstituted phenyl ring as compound 8b of IC50 = 11.06 nM. The p-substitution of the phenyl ring with the electron-donating group OCH3 signified the inhibitory potency as depicted by compound 8a.On the other hand, the hybridization of the parent 6-sulfonoquinoxaline with 3-benzoyl-4-methylthiazoline and 5-acetyl- 3-benzoyl-4-methylthiazoline moieties via the hydrazide linker as congeners 10b, 11b produced a slight reduction in the inhibitory activity compared to the reference drug of IC50s = 6.35, 8.25 nM, respectively. A detectable decrease in the activity was further observed by the 4-methoxyphenyl analogs 10a and 11a of IC50s = 21.63 and 19.45, respectively. It could be noted that the decrease in the thiosemicarbazide length is not favorable for the potency of the desired activity. With respect to the series of pyrrole and isoindoline derivatives 7a–c, the 2-thioxopyrimidine derivatives 9a,b and the 4-oxothiazolidine derivatives 12a,b exhibited the lowest activity of IC50 values ranging from 35.82–57.14 nM. The SAR study of the most potent active congeners is depicted in Figure 4.
Figure 4
Summary of SAR study for PARP-1 suppression effect of the most potent derivatives.
2.2.2. Antiproliferative Activity
In order to find out the relationship between the anticancer potency and the PARP-1 suppression effect, the most effective compounds as PARP-1 inhibitors (4, 5, 8a, 10b, 11b) were further evaluated for their in vitro cytotoxicity against the mutant BRCA1 (MDA-MB-436, breast cancer) using an MTT assay [37]. Olaparib was used as a positive control. The IC50 values of all examined compounds are tabulated in Table 1. The resultant data exhibited that the dimethyl pyrazole compound 5 exhibited the best antiproliferative activity against the examined cancer cell line, being approximately 4 times more potent than the reference drug with IC50 values of 2.57, 8.90 µM, respectively. Furthermore, the tested compounds 8a, 10b, and 11b exhibited an approximately equal activity to that of Olaparib with IC50 values of 10.70, 9.62, and 11.50 µM, respectively. On the contrary, compound 4 displayed the weakest antitumor activity with an IC50 value of 30.30µM.This result represented an outstanding correlation between PARP-1 suppression activity and the anticancer activity of the tested compounds.It has been reported that the frequency and severity of the side effects on normal healthy cells at therapeutic levels are deemed to be critical factors that distinguish different anticancer drugs from each other. Accordingly, the cytotoxic activity of the potent members 5, 8a, 10b, and 11b was evaluated against the normal WI-38 cell line via an MTT assay to determine their safety profiles. It is worth mentioning that the IC50 values of all the representative compounds against the normal cells range from 70.46–81.67 µM, which are 7–8-fold higher than their IC50s values against the cancer cell line, confirming their promising safety profile (Table 1).
2.2.3. Cell Cycle Analysis in MDA-MB-436
Based on its well-balanced biological activity, i.e., promising PARP-1 inhibition and high antiproliferative activity, compound 5 was chosen as a representative example for further examining cellular mechanisms with respect to its impact on cell cycle progression and induction of apoptosis in MDA-MB-436 cells by using the flow cytometric technique [38,39]. In the present work, MDA-MB-436 cells were treated with compound 5 at its IC50 concentration of 2.57 µM and incubated for 48 h. The obtained results were compared to the results obtained by MDA-MB-436 cells incubated with dimethylsulfoxide for 48 h as a negative control. The obtained data are summarized in Figure 5 and Figure 6.
Figure 5
Cell cycle assessment of MDA-MB-436 before and after incubation with compound 5.
Figure 6
The percentage of cells in each phase was quantified using flow cytometry after PI staining of DNA. Data are presented as mean ± SD, n = 3.
Figure 5 represents a decrease in the cell distribution of MDA-MB-436 in the G1 phase from 51.96% (control) to 40.89% (5-treated cells) with a concomitant rise in the percentage of cells in the G2 stage from 25.12% (control) to 30.39% (5-treated cells), which proves that compound 5 arrested the MDA-MB-436 cell cycle at the G2/M stage. In addition, the cell number in the sub-G1 stage was 1.76% (control cells) and 1.66% in (5-treated cells).
2.2.4. Apoptosis Assay in MDA-MB-436 Cells
Annexin V-FITC and PI staining coupled with flow cytometry was utilized for further investigation of the apoptotic impact of compound 5 on MDA-MB-436. Treatment of MDA-MB-436 cells with the IC50 concentration of 5 induced both apoptosis and necrotic effects (Figure 7 and Figure 8). The early apoptotic cell population increased from 0.14% (control) to 0.27% (5-treated cells), the percentage of the late apoptotic cell population increased from 0.98% (control) to 1.02% (5-treated cells), and the necrotic cell population increased from 1.11% (control) to 1.72% (5-treated cells).
Figure 7
(a) Control, (b) Apoptosis/necrosis assessment of MDA-MB-436 after incubation with compound 5. The four quadrants are identified as the necrosis quadrant (Q2-1), late apoptosis quadrant (Q2-2), normal intact cells (Q2-3), and early apoptosis quadrant (Q2-4). Different cell populations were plotted as a percentage of total events. Data are presented as mean ± SD; n = 3.
Figure 8
Apoptosis induction analysis caused by compound 5.
2.2.5. Autophagy Assay
It has been reported that autophagy-induced programmed cell death is a hot topic in the scientific community. It was of interest to study the effect of compound 5 on the autophagy process within MDA-MB-436 cells utilizing Cyto-ID autophagy detection dye coupled with flow cytometry [40,41]. Treatment of the latter cells with 5 increased autophagic cell death by 68.65% (Figure 9).
Figure 9
Autophagic cell death assessment in MDA-MB-436 cells after exposure to compound 5.
2.3. Computational Studies
2.3.1. Molecular Docking
To find out the interaction modes of the most promising quinoxaline congeners 4, 5, 8a, 10b, and 11b, a standard docking protocol was used where Olaparib was utilized as the reference frame of the docking grid. The crustal structure of Olaparib in complex with the catalytic domain of PARP1 was downloaded from the protein databank (www.rcsb.org accessed 15 April 2021). The complex structure was processed with the Protein Preparation Wizard in Maestro to add missing atoms, sidechains, and residues, complete loops, add hydrogen atoms, and adjust bond orders for amino acids and ligands [42,43,44,45,46].The five compounds were docked with high affinity, and the prime MM-GBSA free energy of binding was computed as −93 kcal/mol for Olaparib, −79.3 kcal/mol for compound 8a, −63.9 for compound 10b, −60.9 kcal/mol for compound 5, −54.4 kcal/mol for compound 11b, and −54.3 kcal/mol for compound 4. The compounds fit well in the binding pocket and demonstrated several favorable interactions with the surrounding amino acids. Olaparib is complex with PARP1 in the crystal structure, and it showed the following interactions: Hydrogen bonds with Ser904, Gly863, Ser864, and Tyr896, water-bridged hydrogen bonds with Arg878, Ile879, and π-π contacts with Tyr896 and Tyr907. The compounds showed the following interactions: Compound 8a interacts with hydrogen bonds with Ser904, His862, Asp766, Ser864, π-π contacts with Tyr907 and His862, and cation-π contacts with Arg848; compound 5 showed hydrogen bonds with Gly863 and Ser904, and π-π contacts with Tyr907; compound 10b interacts through hydrogen bonds with Tyr896 and Ser894, a water-bridged hydrogen bond with Arg873, and π-π contacts with Tyr907 and Tyr889; compound 11b showed hydrogen bonds with Ser904, Asn906, and Lys903, water-bridged hydrogen bonds with Met890 and Glu988, π-π contacts with Tyr907, and cation-π contacts with Lys903; and compound 4 interacted through hydrogen bonds with Ser904, Gly863, and Asp766, and π-π contacts with Tyr907 and Tyr889 (Figure 10 and Figure 11 and Table 2).
Figure 10
The 2D interaction diagrams of Olaparib and the five potent compounds 4, 5, 8a, 10b, and 11b.
Figure 11
The binding modes of Olaparib (I), 4 (II), 5 (III), 10b (IV), 8a (V), and 11b (VI). The compounds are shown as sticks, proteins as a cartoon, and the interacting amino acids as lines.
Table 2
The docking results of the tested quinoxaline compounds 4, 5, 10b, 8a, and 11b.
Compd No.
Binding Free Energy (kcal/mol)
Hydrogen-Bond
Water-Bridged Hydrogen Bond
pi-pi/Cation-pi
4
−54.3
Ser904, Gly863, Asp766
Tyr907, Tyr889.
5
−60.9
Gly863, Ser904
Tyr907
8a
−79.3
Ser904, His862, Asp766, Ser864,
-
Tyr907, His862, Arg848
10b
−63.9
Tyr896, Ser894,
Arg873,
Tyr907, Tyr889
11b
−54.4
Ser904, Asn906, Lys903,
Met890, Glu988
Tyr907, Lys903.
Olaparib
−93
Ser904, Gly863, Ser864, Tyr896
Arg878, Ile879
Tyr896 and Tyr907
2.3.2. Prediction ADME parameters
The SwissADME tool was used to calculate the physicochemical properties of the tested compounds 4, 5, 8a, 10b, and 11b (Figure 12) [47,48,49,50,51,52,53]. The compounds showed low to moderate water solubility. Compound 4 showed the highest predicted solubility. Only compound 11b violated Lipinski’s role of five having more than 10 NH and OH groups. All the compounds were computed to have low GIT absorption with the exception of compound 5, which has a promising property. All the compounds are not expected to cross the BBB.
Figure 12
Left: Bioavailability radar from swissADME web tool for compounds 4, 5, 8a, 10b, and 11b. The pink area exhibits the range of the optimal property values for oral bioavailability and the red line is the predicted properties for each compound. Saturation (INSATU), size (SIZE), polarity (POLAR), solubility (INSOLU), lipophilicity (LIPO), and flexibility (FLEX) Right: Predicted Boiled-Egg plot for the compounds.
3. Materials and Methods
3.1. Chemistry
The instruments used for measuring the melting points, spectral data (IR, Mass, 1H NMR, and 13C NMR), and elemental analyses are provided in detail in Supplementary Material.
3.2. PARP-1 Inhibition Assay
PARP-1 enzyme inhibition activity was evaluated using a colorimetric 96-well PARP-1 assay kit (catalog no. 80580) (BPS Bioscience), according to the manufacturer’s protocol. More details are provided in Supplementary Material.
3.3. In Vitro Anticancer Screening
The in vitro cytotoxicity potency was screened against the MDA-MB-436 cancer cell line by MTT assay. The cytotoxicity was estimated as IC50 in µM for the tested compounds and the reference drug Olaparib. More details are provided in Supplementary Material.
3.4. Cell Cycle Analysis
The pre-calculated IC50 of compound 5 was applied to MDA-MB-436 breast cancer cells for 48 h. The cells were treated with trypsin, rinsed two times in PBS, fixed in ice-cold 60% ethanol at 40 °C, and washed again in PBS. More details are provided in Supplementary Material.
3.5. Apoptosis Analysis
MDA-MB-436 cells were treated with compound 5 for 48 h, then treated with trypsin and rinsed twice in PBS. Apoptosis assessment was performed via the “Annexin V-FITC/PI Apoptosis Detection Kit”, “BD Biosciences, San Diego, CA, USA”, as stated by the manufacturer. More details are provided in Supplementary Material.
3.6. Autophagy Analysis
To further confirm the cell death mechanism induced by the drugs, autophagic cell death was quantitatively analyzed using a Cyto-ID Autophagy Detection Kit (Abcam Inc., Cambridge Science Park, Cambridge, UK). More details are provided in Supplementary Material.
3.7. Docking Methodology
The crustal structure of Olaparib in complex with the catalytic domain of PARP1 [42] was downloaded from the protein databank (www.rcsb.org accessed on 15 April 2021). The complex structure was processed with the Protein Preparation Wizard [43,44] in Maestro [45] to add missing atoms, sidechains, and residues, complete loops, add hydrogen atoms, and adjust bond orders for amino acids and ligands. More details are provided in Supplementary Material.
3.8. Chemical Synthesis
3.8.1. Preparation of Quinoxaline-2,3(1H,4H)-dione (1)
A mixture of o-phenylenediamine (5.0 g, 46.3 mmol) in 100 mL 4N HCl and oxalic acid (4.17 g, 46.3 mmol) was heated under reflux for 4 h. Then the reaction mixture was cooled to room temperature and the resulting solid was filtered and washed with ethanol to give the required compound 1 as a white solid according to the reported method [33]. Yield 91.5%, m.p. > 300 °C
3.8.2. Preparation of 2,3-Dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonyl Chloride (2)
The starting compound 1,4-quinoxaline-2,3-dione (1) (1.90 g, 10 mmol) was added portion-wise to chlorosulfonic acid (2 mL, 3 mmol) at 65–90 °C over 3 h. The reaction mixture was cooled to room temperature and poured slowly onto the ice/water mixture. The formed precipitate was collected by filtration and washed with water and dried. The obtained product was crystallized from benzene/petroleum ether (40–60) to give the desired 6-sulfonyl chloride product as a yellowish-white solid according to the reported method [34]. Yield 75%, m.p. 280 (decomposed).
3.8.3. Preparation of 2,3-Dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonohydrazide (3)
A solution mixture of compound 2 (2.60 g, 10 mmol) and hydrazine hydrate 98% (40 mmol, 2 mL) in ethanol (30 mL) was stirred at room temperature for 7 h. The obtained ppt was filtered and crystallized from ethyl alcohol to give the target 6-sulfonohydrazide as a white powder.Yield (65%); mp. 285–287 °C; IR (KBr, cm−1): 3420 (NH2, forked), 3344–3320 (3NH), 3167 (CH, aromatic), 2990 (CH-alicyclic), 1750 (2C=O), 1332, 1138 (SO2); 1H NMR (DMSO-d, δ ppm): 7.01, 7.41 (2d, 2H, aromatic-H, J = 10.03 Hz), 7.85 (s, 2H, NH2, D2O exchangeable), 8.03 (s, 1H, aromatic-H), 8.63, 9.74, 10.01 (3s, 3H, 3NH, exchangeable with D2O); 13C NMR (DMSO-d, δ ppm): 120.29, 123.47, 126.14, 129.05, 133.60, 138.63 (aromatic-C), 155.66 (2C=O)MS, m/z (%): 258 [M+ + 2] (35.08), 256 [M+] (11.06); Analysis for C8H8N4O4S (256.24), Calcd.: %C, 37.50; H, 3.15; N, 21.78; S, 12.51; Found: %C, 37.62; H, 2.93; N, 21.93; S, 12.73.
3.8.4. Preparation of 6-(Substituted pyrazolyl)-sulfonyl-1,4-dihydroquinoxaline-2,3-dione Derivatives 4–6
A mixture of the sulfonohydrazide compound 3 (2.56 g, 10 mmol) with different active methylene reagents, namely ethyl acetoacetate, acetylacetone, and/or diethyl malonate, was refluxed in acetic acid (15 mL) for 8 h. The formed ppt was filtered after cooling, dried, and crystallized from the proper solvent to give the corresponding compounds 4, 5, and 6, respectively.
3.8.5. Preparation of N-Substituted-2,3-dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonamide Derivatives 7a–7c
A mixture of the sulfonohydrazide compound 3 (2.56 g, 10 mmol) and the appropriate acid anhydride derivatives, namely succinic anhydride, maleic anhydride, and/or phthalic anhydride (10 mmol) in acetic acid (15 mL), was heated under reflux for 8 h. The formed precipitate was filtered, dried, and recrystallized from dioxane to obtain the corresponding compounds 7a–c, respectively.
3.8.6. Preparation of 2-((2,3-Dioxo-1,2,3,4-tetrahydroquinoxalin-6-yl)sulfonyl)-N-substituted Hydrazine-1-carbothioamide 8a,b
A solution of the sulfonohydrazide derivative 3 (2.56 g, 10 mmol) and the appropriate isothiocyanate, namely 4-methoxybenzene isothiocyanate and/or benzoyl isothiocyanate (10 mmol) in DMF (20 mL) containing a few drops of triethylamine, was heated under reflux for 6 h. After reaction completion, the mixture was poured onto an ice/H2O mixture and neutralized with HCl. The formed ppt was collected by filtration, washed several times with water, and recrystallized from ethanol.
3.8.7. Preparation of 4-oxo-2-thioxo-3,4-dihydropyrimidin-1(2H)-yl)-2,3-dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonamide Derivatives 9a, 9b
A mixture of the thiosemicarbazide derivatives 8a,b (1 mmol) and diethyl malonate (1.6 mL, 1 mmol) in absolute ethanol (20 )mL was refluxed for 10 h. The reaction mixture was cooled and the formed precipitate was filtered, dried, and recrystallized with ethanol to give the target compounds 9a,b, respectively.
3.8.8. Preparation of N′-(3-aryl-4-methylthiazol-2(3H)-ylidene)-2,3-dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonohydrazide (10a,b) and N′-(5-acetyl-3-aryl-4-methylthiazol-2(3H)-ylidene)-2,3-dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonohydrazide (11a,b)
A mixture of compounds 8a,b (1 mmol) and the appropriate α-halo carbonyl compounds (1mmol), namely chloroacetone (0.82 mL) and/or 3-chloroacetylacetone (1.13 mL) in absolute ethanol (30 mL) containing sodium acetate (1.64 g, 20 mmol) was refluxed for 10–12 h. After cooling, the formed precipitate was filtered, washed with water, dried, and crystallized from ethanol to give the target derivatives 10a,b, and 11a,b, respectively.
3.8.9. Preparation of N’-(3-Aryl-4-oxothiazolidin-2-ylidene)-2,3-dioxo-1,2,3,4-tetrahydroquinoxaline-6-sulfonohydrazide derivatives 12a,b
A mixture of compounds 8a,b (1 mmol) and ethyl bromoacetate (1.11 mL, 1 mmol) in absolute ethanol (30 mL) containing sodium acetate (1.64 g, 20 mmol) was refluxed for 8–10 h. After cooling, the formed precipitate was filtered, washed with water, dried, and crystallized from isopropanol to afford the target derivatives 12a,b, respectively.
The current study deals with the design and synthesis of a novel set of derivatives 3–12a,b bearing the quinoxaline scaffold that is hybridized with various heterocyclic ring systems via a sulfonamide linkage. The new compounds were assessed for their suppression impact against the PARP-1 enzyme using Olaparib as a positive reference drug. Among the examined compounds, 4, 5, 8a, 10b, and 11b displayed the highest PARP-1 inhibitory suppression effect with IC50 values ranging from 2.31 to 8.25 nM, compared to IC50Olaparib of 4.40 nM. The latter compounds were further examined as antiproliferative agents in MDA-MB-436 in comparison with Olaparib as a reference drug. The compounds 5, 8a, 10b, and 11b exhibited promising inhibitory activity with IC50 values ranging from 2.57 to 11.50 µM, compared to IC50Olaparib of 4.40 µM, and confirmed a safety profile against the normal cells’ WI-38 cell lines. Due to the well-balanced activity of compound 5 as a promising PARP-1 inhibitor, as well as the antiproliferative agent, it was chosen as a representative example for further cellular mechanistic investigation regarding its impact on the cell cycle progression and induction of apoptosis in the MDA-MB-436 cell line. Treatment of the latter cells with compound 5 led to cell cycle arrest at the G2/M phase and demonstrated apoptotic and necrotic effects in comparison to the untreated control cells and increased the autophagic cell death (68.65%).Molecular docking of the newly synthesized hybrids 4, 5, 8a, 10b, and 11b in the PARP-1 active sites involved their good accommodation interacting with the various amino acid residues through hydrogen bonding and π-π contacts. The SwissADME tool represented the good GIT absorption of compound 5 and the inability of all the compounds to cross the BBB.
Authors: Richard A Friesner; Robert B Murphy; Matthew P Repasky; Leah L Frye; Jeremy R Greenwood; Thomas A Halgren; Paul C Sanschagrin; Daniel T Mainz Journal: J Med Chem Date: 2006-10-19 Impact factor: 7.446
Authors: Uday Kiran Velagapudi; Marie-France Langelier; Cristina Delgado-Martin; Morgan E Diolaiti; Sietske Bakker; Alan Ashworth; Bhargav A Patel; Xuwei Shao; John M Pascal; Tanaji T Talele Journal: J Med Chem Date: 2019-05-24 Impact factor: 7.446
Authors: Felipe Salech; Daniela P Ponce; Carol D SanMartín; Nicole K Rogers; Carlos Chacón; Mauricio Henríquez; Maria I Behrens Journal: Front Aging Neurosci Date: 2017-10-05 Impact factor: 5.750