- 1Hepatobiliary and Pancreatic Surgery, Cancer Hospital of Zhengzhou University, Zhengzhou City, China
- 2Pharmaceutical Medicinal Chemistry and Drug Design Department, Faculty of Pharmacy (Boys), Al-Azhar University, Cairo, Egypt
- 3Department of Pharmaceutical Chemistry, College of Pharmacy, Taif University, Taif, Saudi Arabia
- 4Zoology Department, Faculty of Science (Boys), Al-Azhar University, Cairo, Egypt
- 5Biophysics Department, Faculty of Women for Arts, Science and Education, Ain Shams University, Cairo, Egypt
- 6Department of Pharmaceutical Medicinal Chemistry, Faculty of Pharmacy, Horus University-Egypt, New Damietta, Egypt
Guided by the structural optimization principle and the promising anticancer effect of the quinoxaline nucleus, a new series of novel HDAC inhibitors were designed and synthesized. The synthesized compounds were designed to bear the reported pharmacophoric features of the HDAC inhibitors in addition to an extra moiety to occupy the non-used vacant deep pocket of the HDAC receptor. The newly prepared compounds were evaluated for their in vitro anti-proliferative activities against HepG-2 and HuH-7 liver cancer cell lines. The tested compounds showed promising anti-proliferative activities against both cell lines. The most active ten candidates (6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12) were further evaluated for their effect on the gene expression levels of Bax as an apoptotic marker and Bcl-2 as an anti-apoptotic one. Moreover, they were evaluated for their ability to inhibit histone deacetylase (HDAC1, HDAC4, and HDAC6) activities. Compound 6c achieved the best cytotoxic activities on both HepG-2 and HuH-7 cell lines with IC50 values of 1.53 and 3.06 µM, respectively, and also it showed the most inhibitory activities on HDAC1, HDAC4, and HDAC6 with IC50 values of 1.76, 1.39, and 3.46 µM, respectively, compared to suberoylanilide hydroxamic acid (SAHA) as a reference drug (IC50 = 0.86, 0.97, and 0.93 µM, respectively). Furthermore, it achieved a more characteristic arrest in the growth of cell population of HepG-2 at both G0/G1 and S phases with 1.23-, and 1.18-fold, respectively, compared to that of the control, as determined by cell cycle analysis. Also, compound 6c showed a marked elevation in the AnxV-FITC apoptotic HepG-2 cells percentage in both early and late phases increasing the total apoptosis percentage by 9.98-, and 10.81-fold, respectively, compared to the control. Furthermore, docking studies were carried out to identify the proposed binding mode of the synthesized compounds towards the prospective target (HDAC4). In silico ADMET and toxicity studies revealed that most of the synthesized compounds have accepted profiles of drug-likeness with low toxicity. Finally, an interesting SAR analysis was concluded to help the future design of more potent HDACIs in the future by medicinal chemists.
Highlights
• Twenty-seven quinoxaline derivatives were designed and synthesized.
• Cytotoxic activities were evaluated against two liver cancer cell lines (HepG-2 and HuH-7).
• In vitro histone deacetylase-4 inhibitory activities were evaluated.
• The effect on cell cycle analysis was studied.
• The effect on apoptosis was studied.
• Apoptotic and anti-apoptotic marker levels (Bcl-2 and Bax) were evaluated.
• Molecular docking studies were carried out against histone deacetylase-4.
• ADMET analysis was done for the newly synthesized derivatives.
• Computational toxicity studies were done.
• SAR analysis was concluded.
Introduction
Liver cancer, also known as hepatic cancer, may start in the liver (Yamashita and Wang, 2013) or spread from elsewhere to the liver, known as liver metastasis (McGuire, 2016). Symptoms include pain in the right side below the rib cage, swelling of the abdomen, yellowish skin, weight loss, and weakness (Liu, 2017). Liver cirrhosis from hepatitis B, hepatitis C, or alcohol is considered the main cause of liver cancer (Perz et al., 2006). The most common types of liver cancer are hepatocellular carcinoma (HCC), which contributes up to 80% of the cases, and cholangiocarcinoma (CCA), which are known as primary liver cancers (PLCs) (McGuire, 2016).
One of the most important challenges facing liver cancer types is their high resistance and poor response to chemotherapy. The produced resistance arises from synergistic interactions among diverse mechanisms of chemoresistance (MOC) in which about 100 genes are involved (Marin et al., 2018).
A very important pathway to fight cancer has appeared after the discovery of tumor angiogenesis by Judah Folkman approximately 50 years ago (Maj et al., 2016). Different effective antiangiogenic agents were approved, mostly targeting vascular endothelial growth factor (VEGF) (El-Helby et al., 2019a; El-Helby et al., 2019b). The antiangiogenic agents either increase the effectiveness of standard chemotherapy or even replace it completely. Now, there are novel strategies other than targeting the VEGF pathway, which are aimed at influencing the molecular factors involved in tumor angiogenesis (Maj et al., 2016). Subsequently, there is a great interest in the development of new antiangiogenic agents that could effectively inhibit tumor vascularization (Vasudev and Reynolds, 2014).
Hypoxia-inducible factor-α (HIF-1α), a central regulator of oxygen detection and adaptation at the cellular level, and its transcriptional activity are the key mediator of VEGF activity. Both HIF-1α and VEGF are crucial to angiogenesis and can be regulated by post-translational modifications (PTMs), including acetylation by histone acetyltransferases (HATs) and deacetylation by histone deacetylases (HDACs). Subsequently, many studies indicated HDAC inhibitors (HDACIs) as promising antiangiogenic compounds and recommended them as an effective class of anticancer therapeutics (Ellis et al., 2009; Deng et al., 2020).
The first HDACI approved by the FDA for the treatment of cutaneous T-cell lymphoma is vorinostat, followed by depsipeptide for the same purpose (Ververis et al., 2013). This was followed by the FDA approval of two other drugs, belinostat and panobinostat. Despite the success of HDACIs in the treatment of leukemias, they are still failing in the case of solid tumors (Sangwan et al., 2018). Many structurally diverse HDACIs are in different phases of clinical trials as a monotherapy and/or in combination with other anticancer agents (Marks and Xu, 2009).
Roquinimex I (linomide) is a quinoline derivative immunostimulant by increasing natural killer (NK) cell activity and macrophage cytotoxicity. It also inhibits angiogenesis and reduces the secretion of TNF-α. It was indicated for the treatment of some cancer types and autoimmune diseases, such as multiple sclerosis, and prevention of autoimmune diabetes mellitus (Gross et al., 2001; Banu et al., 2017). But several trials have been terminated due to observed cardiovascular toxicity (Tan et al., 2000). Tasquinimod II is a second-generation quinoline-3-carboxamide agent that is orally active in HDACI and is currently in phase III clinical trials for the treatment of castration-resistant prostate cancer (Gupta et al., 2014). It counteracts cancer development by inhibiting angiogenesis and metastasis and on the other hand by modulating the immune system (Isaacs et al., 2006; Jennbacken et al., 2012; Isaacs et al., 2013). Meanwhile, the mode of action for tasquinimod II is not fully understood, and several studies demonstrated its ability to interfere with tumor angiogenesis, cytokine production, macrophage infiltration, and autoimmune/inflammatory diseases. Furthermore, several studies including preclinical ones are required to investigate the real mechanisms of action for tasquinimod II (Gupta et al., 2014). Laquinimod III is an experimental immunomodulator and is being tested as an oral treatment for multiple sclerosis (MS). Phase III clinical trials for MS started in December 2007 and showed a slow progression of disability and reduction in the rate of relapse in patients with relapsing-remitting multiple sclerosis (Brück and Zamvil, 2012; Banu et al., 2017) (Figure 1).
The discovery of the histone deacetylases (HDACs) role has promised to be new hope for the treatment of various malignancies (Cappellacci et al., 2020). Thus, the intention for making new potent and safe anticancer therapeutic agents with minimal side effects is still a major concern for researchers nowadays. Besides, quinoxaline moiety is an important N-containing heterocycle in organic synthesis and drug discovery (Corona et al., 2009; Ibrahim, 2012; Ibrahim et al., 2013; Ibrahim et al., 2015; Eissa et al., 2018; Ibrahim et al., 2018; Abbass et al., 2020) due to its large scope of biological activity, especially antitumor activities (Corona et al., 2009). So, our goal is to design and synthesize new quinoxaline derivatives as antiangiogenic agents targeting HDAC with promising effects against liver cancer.
Rational of Molecular Design
Investigation of the common pharmacophoric features shared by various HDACIs revealed that most of them have three main features: i) a heterocyclic aromatic cap that occupies the narrow tubular pocket of HDAC and containing at least one H-bond acceptor. This H-bond acceptor constitutes the Zn2+-binding group, ii) an amide linker which occupies the linker region between the zinc-binding region and the hydrophobic tail region and contains one H-bond donor that forms an H-bond with the crucial amino acid His132, iii) a terminal hydrophobic group which protrudes to outside, similar to that of the co-crystallized inhibitor, and helps in the stabilization of the remaining part of the compound (Figure 2).
Interestingly, the newly synthesized quinoxaline candidates were found to form nearly the same binding mode compared to the co-crystallized inhibitor (TSA) with an additional extra positioning of the quinoxaline moiety inside the narrow tubular pocket of HDAC, which provides a more promising fitting as the N1 of quinoxaline, and it was found to be the zinc-binding group, as described in Figure 2.
In continuation of our previous works (Eissa et al., 2016; Eldehna et al., 2017; El-Naggar et al., 2017; Eissa et al., 2018; Gaber et al., 2018; Ibrahim et al., 2018; Elmetwally et al., 2019; Metwaly et al., 2019; Alesawy et al., 2020; Al-Karmalawy and Khattab, 2020; Eliaa et al., 2020; El-Zahabi et al., 2020; Ghanem et al., 2020; Al-Karmalawy and Elshal, 2021 and Elshal; El-Shershaby et al., 2021b; Khattab and Al-Karmalawy, 2021; Samra et al., 2021), that proved potential anticancer activities of novel chemical agents, a new series of quinoxaline derivatives were designed and synthesized to have the essential pharmacophoric features of the reported and clinically used HDACIs to get more potent antitumor molecules. The main core of our molecular design rational comprised bio-isosteric modification strategies of the lead compounds (roquinimex, tasquinimod, and laquinimod) at three different positions (Figure 3).
FIGURE 3. The basic structural requirements for roquinimex, tasquinimod, and laquinimod as reported HDAC inhibitors and their existence in our newly synthesized compounds.
The first position was the heterocyclic aromatic cap. Two different substituted quinoxaline moieties were used as follows: 2-chloroquinoxaline (compounds 6a-l, 7a,b, 8, and 9) and 2-methoxyquinoxaline (compounds 10a-i, 11, and 12). The choice of quinoxaline was based on some bio-isosteric considerations: i) the bicyclic structure of quinazoline core which occupies the narrow tubular pocket of HDAC superior to the co-crystallized inhibitor (Hou et al., 2003), ii) the nitrogen atoms serve as hydrogen-bond acceptors and one of them is acting as the zinc-binding region conferring excellent HDAC inhibitory activity. The second position was the linker (spacer) region. The amide linker of the lead compounds was modified to be methylenehydrazine (compounds 6a-l, 8, 9, 10, 11, and 12) or ethylidenehydrazine (compounds 7a,b). The third position was the terminal hydrophobic group. We used many different aromatic moieties to play the role of the terminal hydrophobic group. The wide variety of modifications enabled us to study the SAR of these compounds as effective anticancer agents with potential HDAC inhibitory activities, which is considered a crucial objective of our work. All modification pathways and molecular design rational were illustrated and summarized in Figure 3 and Supplementary Figure SI1.
Also, to explore and emphasize the mechanism of action of the synthesized compounds, molecular docking studies were conducted to understand the expected binding interactions of the target compounds with HDAC active sites. For the same purpose, other studies regarding the ability of the most active compounds to induce apoptosis and arrest cell cycle growth have been done.
Results and Discussion
Chemistry
For the synthesis of the target compounds, the reaction sequences were illustrated in Schemes 1, 2. At first, o-phenylenediamine was treated with oxalic acid in the presence of 4 N HCl to afford 2,3-(1H,4H)-quinoxalinedione 3. The latter was treated with thionyl chloride in the presence of a catalytic amount of DMF to produce 2,3- dichloroquinoxaline 4 (Romer, 2009). Stirring of the 2,3 dichloroquinoxaline with hydrazine hydrate in the presence of TEA at room temperature afforded 2-chloro-3-hydrazinylquinoxaline 5. Reflux of 2-chloro-3-hydrazinylquinoxaline with commercially available aldehydes and ketones, namely, benzaldehyde, 2-chlorobenzaldehyde, 4-chlorobenzaldehyde, 2,4-dichlorobenzaldehyde, 2,6-dichlorobenzaldehyde, 2-methoxybenzaldehyde, 3,4-dimethoxybenzaldehyde, 3,4,5-trimethoxybenzaldehyde, 3-nitrobenzaldehyde, 4-nitrobenzaldehyde, 2-hydroxybenzaldehyde, 4-N,N-dimethylaminobenzaldehyde, acetophenone, 2-hydroxyacetophenone, 2-hydroxy-1-naphthaldehyde, and isatin, in the presence of a catalytic amount of acetic acid produces target compounds 6a-l, 7a,b, 8, and 9, respectively (Scheme 1). The formation of benzylidene derivatives 6a-l, 7a,b, 8, and 9 was confirmed by 1H NMR spectra, which showed the appearance of a singlet signal for benzylidene methine proton in the range of δ 7.80–8.00 ppm. This methine proton was detected also in both 13C NMR spectra resonating around 127.00 ppm. Besides, 1H NMR and 13C NMR spectra exhibited increased aromatic signals indicating that the condensation reactions were completed.
The hydrazone derivatives 6c-h, 6j-l, 8, and 9 were treated with sodium methoxide to obtain final methoxy derivatives 10a-i, 11, and 12, respectively (Scheme 2). The formation of methoxy derivatives 6c-h, 6j-l, 8, and 9 was confirmed by 1H NMR spectra, which showed the appearance of a singlet signal for methoxy around δ 4.10 ppm. This methoxy group was detected also in both 13C NMR spectra resonating around 54.30 ppm.
Biological Evaluation
Anti-Proliferative Activities Towards Liver Cancer Cell Lines (HepG-2 and HuH-7)
Anti-proliferative of the synthesized compounds was examined towards liver cancer HepG-2 and HuH-7 cell lines, using MTT assay (Skehan et al., 1990). The IC50 values for the synthesized derivatives were compared with doxorubicin as a positive control (Table 1).
TABLE 1. Anti-proliferative activities for 6a-l, 7a,b, 8, 9, 10a-i, 11, and 12 compounds against HepG-2 and HuH-7 liver cancer cell lines.
To conclude a reasonable structure-activity relationship (SAR) of quinoxaline-based HDACIs, five series: 6, 7, 8, 9, and 10, was designed through hybridization between quinoxaline and different aldehydes, acetophenones, 2-hydroxy-1-naphthaldehyde, or isatin, and at the same time, substitution of 2-chloro group with 2-methoxy one to study their different effects.
Analyzing the IC50 values of the newly synthesized compounds on both HepG-2 and HuH-7 liver cancer cell lines revealed the following interesting results. Compounds of 6 series showed the best promising cytotoxic activities with IC50 ranging from 1.53 to 18.70 µM against the two cell lines except for compound 6i with a 3-NO2 side chain (IC50 = 63.87 and 65.72 µM, respectively). Compound 6c with a 4-Cl side chain showed the best cytotoxic activities against both cell lines among all the synthesized series (IC50 = 1.53 and 3.06 µM, respectively).
On the other hand, the derivatives of 10 series incorporating 2-OCH3 group at 2-position of quinoxaline nucleus showed a decrease in the cytotoxic activities ranging from 5.56 to 65.46 µM against the two cell lines, except for compound 10h incorporating 2-OH group in the side chain which elevated greatly its cytotoxic activity (IC50 = 5.56 and 6.72 µM) against both HepG-2 and HuH-7 cell lines, respectively.
The two compounds of the 7 series showed variable anticancer activities. On the one hand, compound 7b incorporates 2-OH (IC50 = 4.03 and 6.19 µM against HepG-2 and HuH-7, respectively). On the other hand, the unsubstituted member 7a showed decreased activities against the two cell lines (IC50 = 28.52 and 37.11 µM against HepG-2 and HuH-7, respectively).
The presence of 2-hydroxy-1-naphthyl moiety in the side chain of compound 8 made it superior to compound 9 incorporating isatin moiety against the two cell lines. Compound 8 (IC50 = 3.88 and 6.14 µM against HepG-2 and HuH-7, respectively) was 2.98-fold more active than compound 9 (IC50 = 11.58 and 16.43 µM against HepG-2 and HuH-7, respectively).
Insertion of the 2-Cl group in quinoxaline moiety of compound 8 (IC50 = 3.88 and 6.14 µM against HepG-2 and HuH-7, respectively) makes it more active than the corresponding member 11 incorporating 2-OCH3 (IC50 = 18.47 and 20.31 µM against HepG-2 and HuH-7, respectively). Compound 8 showed increase activities by 4.76-fold against the HepG-2 cell line and 3.30-fold against the HuH-7 cell line. Contrary to the above, the insertion of the 2-OCH3 group in quinoxaline moiety of compound 12 (IC50 = 1.99 and 3.08 µM against HepG-2 and HuH-7, respectively) make it more active than the corresponding member 9 incorporating 2-Cl (IC50 = 11.58 and 16.43 µM against HepG-2 and HuH-7, respectively) Compound 12 showed increase activities by 5.81-fold against the HepG-2 cell line and 5.33-fold against the HuH-7 cell line.
Histone Deacetylase Inhibitory Activities
The most active compounds (6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, 12) were further evaluated for their ability to inhibit histone deacetylase activities (HDAC1, HDAC4, and HDAC6) (Table 2). The HDAC4 enzyme was chosen since it is widely seen in early tumorigenesis through deacetylation and demethylation of the residues (especially the lysine residues) of the histone H4 as the most important feature in cancer prognosis (Montgomery et al., 2007). Moreover, both HDAC1 and HDAC6 enzymes were evaluated to further confirm the HDAC inhibitory activities of the newly designed and synthesized compounds. The tested compounds showed promising HDAC inhibitory activities compared to the reference drug, suberoylanilide hydroxamic acid (SAHA). Especially compounds 6c, 6g, and 12 achieved the highest inhibitory activities against HDAC1, HDAC4, and HDAC6 as well. Again, compound 6c showed better HDAC inhibitory activities compared to the reference standard that clarifies greatly its promising targeting as an HDAC inhibitor after further preclinical and clinical studies.
Cell Cycle Analysis
The most compound 6c incorporating 4-Cl benzylidene side chain was further evaluated through cell cycle analysis in HepG-2 cell line. Such a test was carried out to detect the exact phase at which cell cycle arrest takes place. The HepG-2 cells were treated with compound 6c at a concentration of 1.53 µM equal to its IC50, and its impact on the different phases of cell growth was recorded. Treatment of HepG-2 cells with 6c showed a significant decline in the cell population at both G0/G1 and S phases with 36.75% (1.23-fold) and 30.32% (1.18-fold), respectively, compared to that of the control which was 45.40 and 35.95% (Table 3 and Figure 4). Besides, there was an increase in the population of cells at the G2/M phase by 1.76-fold compared to the control. This indicates that the compound halted the cell cycle progression of HepG-2 cells in the G0/G1 phase.
Annexin V-FITC Apoptosis Assay
It is well-known that cell death may be due to programmed apoptosis or uncontrolled necrosis. Annexin V-based flow cytometry assay as a helpful strategy to determine the exact cause of death was carried out. Since compound 6c showed the highest anticancer activity toward the HepG-2 cell line, it was tested to investigate its apoptotic effect. The results revealed that the treatment of HepG-2 cells with 6c with a concentration of 1.53 µM showed a marked elevation in the AnxV-FITC apoptotic cells percentage in both early (from 0.62 to 6.19%, respectively) and late apoptosis (from 0.53 to 5.73%, respectively) phases (Figure 5). This refers to an increase in the total apoptosis percentage by 9.98-, and 10.81-fold, respectively, compared to the control. This confirms that the cytotoxic activity of compound 6c is due to programmed apoptosis, and not to nonspecific necrosis.
Apoptotic and Anti-Apoptotic Marker Levels (Bcl-2 and Bax)
To prove that the synthesized compounds exert their cytotoxic effects through driving cells to apoptosis, the effect of most active compounds against the level of Bax (as an apoptotic marker) and Bcl-2 (as an anti-apoptotic marker) was evaluated. Compounds 6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12 were utilized in this test.
Auspiciously, the tested members greatly elevated the level of proapoptotic Bax protein by a range from 164.3 to 274.3 pg/ml (from 6.74 to 11.25-fold increase), compared to the control (24.37 pg/ml). On the other hand, all the tested compounds were able to decrease the level of the anti-apoptotic Bcl-2 protein by a range from 2.91 to 1.27 pg/ml (from 1.73 to 3.96-fold decrease), compared to the control (5.04 pg/ml). Again, compound 6c among all of the tested derivatives achieved the highest increase in the level of Bax protein (274.3 pg/ml) and the lowest decrease in the level of Bcl-2 protein (1.27 pg/ml), compared to the control (24.37 and 5.04 pg/ml), respectively (Table 4).
TABLE 4. Effect of compounds 6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12 on the expression levels of Bcl-2 and Bax in HepG-2 cells.
In silico Studies Results
Docking Studies
Docking studies of the synthesized compounds were carried out to rationalize the obtained biological results and to understand the proposed binding mode of such compounds with the prospective target (HDAC). Trichostatin A (TSA) as an HDAC inhibitor was used as a reference drug in the docking studies.
At first, a validation process was performed for the target receptor by running a redocking process for only the co-crystallized inhibitor, and a low RMSD value indicated the valid performance (RMSD = 0.57) (Figure 6; Davis and Baker, 2009; Alnajjar et al., 2020; Abo Elmaaty et al., 2021).
FIGURE 6. 2D and 3D representations for the redocked co-crystallized TSA antagonist inside the HDAC receptor pocket.
The results of docking studies of TSA revealed that it occupied the tube-shaped pocket of the receptor with nearly a fingerprint binding mode compared to the co-crystallized one. It was found that the phenyl ring of TSA is located outside the receptor pocket and fits in the narrow portion of the pocket through its five-carbon-long branched aliphatic chain making multiple van der Waal interactions with the surrounding hydrophobic moieties lining the pocket. The hydroxamic acid group at the end of its aliphatic chain reaches the polar bottom part of the pocket, where it binds the zinc in a bidentate fashion and also contacts the crucial active site residues (His131 and His132). It recorded a binding score of −11.95 kcal/mol and an RMSD of 1.23 Å. It bound Zn2+ metal in a similar bidentate manner through its terminal charged and carbonyl oxygen atoms with 1.93 and 2.51 Å, respectively. Moreover, it formed one ionic bond through its charged oxygen atom with His131 with 3.17 Å, and two hydrogen bonds with Tyr297 and His132 with 2.67 and 3 Å through the carbonyl oxygen and amidic NH groups, respectively. Finally, it formed a pi-H interaction through its phenyl group with His170 with 4.74 Å. It was noted that the hydroxamic acid group at the end of the aliphatic chain of TSA fitted the polar bottom part of the pocket, leaving its hydrophobic moiety outside (Figure 7).
FIGURE 7. 3D representation and positioning for the docked TSA inside the tubular pocket of the HDAC receptor.
Interestingly, the synthesized derivatives were found to form nearly the same binding mode compared to the co-crystallized inhibitor (TSA) with an additional extra positioning of the quinoxaline moiety inside the narrow tubular pocket of HDAC, which provides a more promising fitting as the N1 of quinoxaline was found to be the zinc-binding group.
Compound 6K, as a representative example, was fitted inside the deep pocket through its quinoxaline ring with a binding score of − 6.46 kcal/mol and RMSD of 0.77 Å. It bound Zn2+ metal in a unidentate manner through its N1 atom of quinoxaline moiety with 2.63 Å. Furthermore, it formed a hydrogen bond with His132 at 2.96 Å through its linker NH group of the hydrazide moiety (Figure 8).
FIGURE 8. 3D representation and positioning for the docked 6k inside the tubular pocket of the HDAC receptor.
Besides, the quinoxaline ring of compound 10h was positioned inside the narrow tube-like pocket of HDAC forming a binding score of 6.45 kcal/mol and RMSD of 1.54 Å. As mentioned before as a general binding mode of our quinoxaline derivatives, it was stabilized in HDAC pocket through the formation of a unidentate bond with Zn2+ metal by N1 of quinoxaline ring with 2.58 Å. Moreover, it formed an H-bond with His132 through its NH group of hydrazide linker with 3.12 Å (Figure 9).
FIGURE 9. 3D representation and positioning for the docked 10h inside the tubular pocket of the HDAC receptor.
ADMET Analysis
ADMET studies were carried out for the synthesized compounds, including some descriptors. The predicted descriptors are listed in the Supplementary Table SI1.
ADMET-BBB penetration results revealed that compounds 6i, 6j, 10g, and 12 have low or very low levels; so that, these compounds were expected to be safe to CNS. The other compounds were predicted to have a very high, high, or medium level of BBB penetration. All the tested compounds showed low to very low range levels of ADMET aqueous solubility.
Intestinal absorption is defined as the percentage absorbed of a compound from the gut wall (Mannhold et al., 2012; Zaki et al., 2021). A well-absorbed compound can penetrate the bloodstream in humans by at least 90% (Klopman et al., 2002). According to ADMET studies, the absorption levels of all compounds appeared in the good range.
The cytochrome P450 2D6 (CYP2D6) model predicts CYP2D6 enzyme inhibition using a 2D chemical structure as input. CYP2D6 inhibition experiment is required as part of the regulatory procedures in the drug discovery and development process (Roy and Roy, 2009). All the tested compounds were predicted to be non-inhibitors of CYP2D6 except compounds 6b, 6c, 6d, 6e, and 10b. Consequently, a liver dysfunction side effect is not expected upon administration of these compounds. The plasma protein-binding model predicts whether a compound is likely to be highly bound (≥90% bound) to carrier proteins in the blood (Ghafourian and Amin, 2013). All compounds were expected to bind plasma protein over 90% (Figure 10).
Toxicity Studies
Toxicity prediction was carried out for the synthesized compounds based on the validated and constructed models in Discovery Studio software (Xia et al., 2004; Agoni et al., 2020).
It is worth mentioning that most compounds showed in silico low adverse effects and toxicity against the tested models. Concerning FDA rodent carcinogenicity, all the tested compounds were predicted to be non-carcinogenic except compounds 8, 10f, and 12 which were single carcinogens. For the carcinogenic potency TD50 rat model, the tested compounds showed TD50 values ranging from 1.797 to 52.581 mg/kg body weight/day. Regarding the rat maximum tolerated dose model, the compounds showed maximum tolerated dose with a range of 0.082–0.760 g/kg body weight. Additionally, all compounds were non-toxic against the developmental toxicity potential model except compound 6h. For the rat oral LD50 model, all compounds showed low oral LD50 values (from 0.102 to 1.109 mg/kg body weight/day). For the rat chronic LOAEL model, the compounds showed LOAEL values ranging from 0.055 to 0.413 g/kg body weight. Moreover, all compounds were predicted to be mild and non-irritant against ocular irritancy and skin irritancy models, respectively, as represented in Supplementary Table SI2.
Structure-Activity Relationship Studies
Studying the structure-activity relationship of our newly synthesized tested quinoxaline candidates according to their IC50 values towards liver cancer cell lines (HepG-2 and HuH-7) showed the following interesting results:
Generally, 2-chloro quinoxaline derivatives (6) with different substituted benzylidene side chains were found to exert the most promising cytotoxic activityranging from very strong to strong (except for 3-NO2 one). Besides, 2-methoxy quinoxaline compound (12) with isatin side chain maintains a very strong cytotoxic activity as well. Moreover, 2-chloro quinoxaline derivative with a 2-hydroxy-1-naphthaldehyde side chain (8) showed a very strong anticancer activity on both cell lines which exceeds that of some derivatives of (6). Furthermore, the introduction of a 2-hydroxy acetophenone side chain (7b) maintains the very strong cytotoxic activity of 2-chloro quinoxaline moiety in contrast to the plain acetophenone derivative (7a) which decreased the cytotoxic effect greatly. Also, the presence of an isatin side chain on 2-chloro quinoxaline (9) achieved a strong anti-proliferative effect. On the other hand, the analog of the previously mentioned compound (8) with a very strong cytotoxic activity, 2-methoxy quinoxaline derivative (11), showed only a strong cytotoxic activity. Finally, most of the 2-methoxy quinoxaline derivatives (10) compared to their 2-chloro analogs showed weak anti-proliferative activities except those with 2-OH, 4-NO2 and 2-OCH3 benzylidene side chains which showed very strong, strong, and moderate cytotoxic activities, respectively (Figure 11).
FIGURE 11. Structure-activity relationships of the newly synthesized quinoxaline candidates as HDACIs.
Conclusion
Twenty-seven new quinoxaline derivatives were synthesized with different substitutions to study their SAR as promising anticancer candidates targeting the HDAC enzymes (HDAC1, HDAC4, and HDAC6 subtypes). Twenty compounds showed cytotoxic effects ranging from very strong to moderate against two liver cancer cell lines (HepG-2 and HuH-7). Then, the most active ten compounds (6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12) were further evaluated as HDAC1, HDAC4, and HDAC6 inhibitors and revealed IC50 values ranging from 1.39 to 7.21 µM, compared to the reference drug, suberoylanilide hydroxamic acid (SAHA). Moreover, the most active compound 6c with a 4-Cl benzylidene side chain was further evaluated through cell cycle analysis in HepG-2 cell line and showed a significant decline in the cell population at both G0/G1 and S phases with 36.75% (1.23-fold) and 30.32% (1.18-fold), respectively, compared to that of the control which was 45.40 and 35.95%. Furthermore, it was subjected to annexin V-based flow cytometry assay and achieved an increase in the total apoptosis percentage by 9.98-, and 10.81-fold, respectively, compared to the control. Finally, the aforementioned most active ten compounds greatly elevated the level of proapoptotic Bax protein by a range from 6.74- to 11.25-fold increase, compared to the control. On the other hand, they were able to decrease the level of the anti-apoptotic Bcl-2 protein by a range from 1.73- to 3.96-fold decrease, compared to the control. Again, compound 6c among all of the tested derivatives achieved the highest increase in the level of Bax protein (274.3 pg/ml) and the lowest decrease in the level of Bcl-2 protein (1.27 pg/ml), compared to the control (24.37 and 5.04 pg/ml), respectively. Interestingly, docking studies revealed nearly the same binding mode compared to the co-crystallized inhibitor (TSA) with an additional extra positioning of the quinoxaline moiety inside the narrow tubular pocket of HDAC, which may explain the reason behind the previously discussed promising biological results, and confirm our proposed mechanism of action for them as HDACIs. Also, the SAR studies of our diverse synthesized derivatives based on their biological results may give a clear spot on the essential pharmacophoric features required for targeting HDAC as well in the future by medicinal chemists.
Materials and Methods
Chemistry
General
All melting points were carried out by the open capillary method on a Gallen Kamp apparatus. The infrared spectra were recorded on a Pye Unicam SP 1000 IR spectrophotometer using the potassium bromide disk technique. Proton and carbon magnetic resonance NMR spectra were recorded on a Bruker AVANCE-III 400 MHz-NMR spectrometer. TMS was used as an internal standard, and chemical shifts were measured on a d scale (ppm). The mass spectra were recorded on Varian MAT 311-A (70 e.v.) and Direct Inlet unit (DI-50) of SHIMADZU GC/MSeQP5050A. The reactions were monitored by thin-layer chromatography (TLC) using TLC sheets precoated with UV fluorescent silica gel Merck 60 F254 plates and were visualized using a UV lamp and different solvents as mobile phases. 2,3-(1H,4H)-Quinoxalinedione 3, 2,3- dichloroquinoxaline 4, 2-chloro-3-hydrazinylquinoxaline 5 were obtained according to the reported procedures(Romer, 2009).
General Procedure for the Synthesis of Compounds 6a-l, 7a,b, 8, and 9
Equimolar amounts of compound 5 (0.7 g, 0.002 mol) and appropriate carbonyl derivatives (0.002 mol), namely, benzaldehyde, 2-chlorobenzaldehyde, 4-chlorobenzaldehyde, 2,4-dichlorobenzaldehyde, 2,6-dichlorobenzaldehyde, 2-methoxybbenzaldehyde, 3,4-dimethoxybenzaldehyde, 3,4,5-trimethoxybenzaldehyde, 3-nitrobenzaldehyde, 4-nitrobenzaldehyde, 2-hydroxybenzaldehyde, acetophenone, 2-hydroxyacetophenone, 2-hydroxy-1-naphthaldehyde, and isatin were refluxed in absolute ethanol (25 ml) in the presence of a catalytic amount of glacial acetic acid for the appropriate time, and the reaction was followed up by TLC. The mixture was cooled, and the formed solid product was filtered, dried, and crystallized from ethanol to afford the corresponding compounds 6a-l, 7a,b, 8, and 9, respectively.
2-(2-Benzylidenehydrazinyl)-3-chloroquinoxaline 6a
Canary yellow crystal (yield, 89%); m. p. = 228–230°C; IR (KBr, cm1): 3,265 (NH), 3,056 (CH aromatic), 1,607 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 10.99 (s, 1H, NH), 8.62 (s, 1H, NCH), 8.32–7.63 (m, 3H, Ar-H, H2 and H6 of benzene, H8 of quinoxalne), 7.51 (dd, J = 14.9, 7.8 Hz, 5H, Ar-H, H5, H6 and H7 of quinoxalne, H3 and H5 of benzene), 7.16 (s, 1H, Ar-H, H4, of benzene); 13C NMR (101 MHz, DMSO-d6) δ 156.44 (Ar-C, C3 of quinoxaline), 146.82 (NCH), 140.42 (Ar-C2 of quinoxaline), 135.44 (Ar-C), 133.53 (Ar-C), 130.85(Ar-C), 129.13(Ar-C), 128.78 (Ar-C), 127.37 (Ar-C), 126.02 (Ar-C), 125.31(Ar-C), 123.05 (Ar-C), 115.36 (Ar-C); C15H11ClN4 (282.73).
2-Chloro-3-(2-(2-chlorobenzylidene)hydrazinyl)quinoxaline 6b
Bright yellow crystal (yield, 83%); m. p. = 224–226°C; IR (KBr, cm1): 3,361 (NH), 3,006 (CH aromatic), 1,610 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.36 (s, 1H, NH), 9.07 (s, 1H, NCH), 8.61 (d, J = 5.2 Hz, 1H, Ar-H, H5 of quinoxaline), 8.11 (d, J = 7.7 Hz, 1H, Ar-H, H8 of quinoxaline), 7.58 (t, J = 7.0 Hz, 1H, Ar-H, H6 of quinoxaline), 7.56–7.49 (m, 2H, Ar-H, H7 of quinoxaline and H6 of benzene), 7.46 (d, J = 7.8 Hz, 2H, Ar-H, H3 and H5 of benzene), 7.24–7.12 (m, 1H, H4 of benzene); C15H10Cl2N4 (317.17). 13C NMR (101 MHz, DMSO-d6) δ 152.16, 145.05, 142.73, 140.94, 134.08, 133.63, 132.84, 131.34, 130.38, 129.27, 127.83, 127.34, 125.60, 123.30, 115.48.
2-Chloro-3-(2-(4-chlorobenzylidene)hydrazinyl)quinoxaline 6c
Orang powder (yield, 84%); m. p. = 184–186°C; 1H NMR (400 MHz, DMSO-d6) δ 10.82 (s, 1H, NH), 8.61 (s, 1H, NCH), 7.87 (d, J = 8.2 Hz, 2H, Ar-H, H5 and H8 of quinoxaline), 7.75 (t, J = 9.1 Hz, 2H, Ar-H, H6 and H7of quinoxaline), 7.61 (d, J = 7.8 Hz, 1H, Ar-H, H2 of benzene), 7.52 (d, J = 8.1 Hz, 2H, 2H, Ar-H, H4 and H5 of benzene), 7.40 (d, J = 7.7 Hz, 1H, Ar-H, H6 of benzene); 13C NMR (101 MHz, DMSO-d6) 13C NMR (101 MHz, DMSO-d6) δ 145.05 (Ar-C, C3 of quinoxaline), 144.83, (NCH), 138.41 (Ar-C, C2 of quinoxaline), 136.33 (Ar-C), 134.96 (Ar-C), 134.08 (Ar-C), 131.06 (Ar-C), 129.55 (Ar-C), 129.29 (Ar-C), 128.06 (Ar-C), 125.40 (Ar-C).; C15H10Cl2N4 (317.17).
2-Chloro-3-(2-(2,4-dichlorobenzylidene)hydrazinyl)quinoxaline 6d
Yellow crystal (yield, 80%); m. p. = 182–184°C; 1H NMR (400 MHz, DMSO-d6) δ 11.50 (s, 1H, NH), 9.03 (s, 1H, NCH), 8.12 (d, J = 8.5 Hz, 1H, Ar-H, H6 of benzene), 7.89–7.83 (m, 1H, Ar-H, H3 of benzene), 7.75 (d, J = 7.3 Hz, 2H, Ar-H, H5 and H8 of quinoxaline), 7.58 (m, J = 8.3 Hz, 3H, Ar-H, H6 and H7 of quinoxaline, H5 of benzene); 13C NMR (101 MHz, DMSO-d6) δ 145.27 (Ar-C, C3 of quinoxaline), 142.25 (Ar-C, C2 of quinoxaline), 140.93 (NCH), 137.54 (Ar-C), 136.71 (Ar-C), 135.09 (Ar-C), 134.02 (Ar-C), 131.62 (Ar-C), 131.20 (Ar-C), 129.86 (Ar-C), 128.47 (Ar-C), 128.03 (Ar-C), 126.94 (Ar-C).; C15H9Cl3N4 (351.62).
2-Chloro-3-(2-(2,6-dichlorobenzylidene)hydrazinyl)quinoxaline 6e
Orange crystal (yield, 81%); m. p. = 220–222°C; 1H NMR (400 MHz, DMSO-d6) δ 11.50 (s, 1H, NH), 8.83 (s, 1H, NCH), 7.87 (d, J = 8.2 Hz, 1H, Ar-H, H5 of quinoxaline), 7.72 (t, J = 7.6 Hz, 1H, Ar-H, H8 of quinoxaline), 7.59 (t, J = 8.2 Hz, 3H, Ar-H, H6 of quinoxaline, H3 and H5 of benzene), 7.50–7.42 (m, 1H, Ar-H, H4 of benzene); 13C NMR (101 MHz, DMSO-d6) δ 13C NMR (101 MHz, DMSO-d6) δ 145.37 (Ar-C, C3 of quinoxaline), 142.69 (Ar-C, C3 of quinoxaline), 140.88 (Ar-C, C2 of quinoxaline), 137.58 (NCH), 136.59 (Ar-C), 134.44 (Ar-C), 131.60 (Ar-C), 131.15 (Ar-C), 129.52 (Ar-C), 128.01 (Ar-C), 127.02 (Ar-C); C15H9Cl3N4 (351.62).
2-Chloro-3-(2-(2-methoxybenzylidene)hydrazinyl)quinoxaline 6f
Canary yellow crystal (yield, 88%); m. p. = 210–212°C; IR (KBr, cm1): 3,369 (NH), 3,015 (CH aromatic), 2,952 (CH aliphatic), 1,627 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 10.95 (s, 1H, NH), 8.82 (s, 1H, NCH), 8.39 (d, J = 7.6 Hz, 1H, Ar-H, H5 of quinoxaline), 8.24 (s, 1H, Ar-H, H6 of benzene), 7.51 (d, J = 7.4 Hz, 1H, Ar-H, H8 of quinoxaline), 7.46 (t, J = 7.7 Hz, 1H, Ar-H, H6 of quinoxaline), 7.32 (d, J = 7.4 Hz, 1H, Ar-H, H3 of benzene), 7.13 (t, J = 8.5 Hz, 1H, Ar-H, H7 of quinoxaline), 7.07–7.03 (m, 2H, Ar-H, H4 and H5 of benzene), 3.89 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 158.51 (Ar-C, C3 of quinoxaline), 151.40 (Ar-C, C2 of benzene), 149.37 , 140.31 (NCH), 133.88 (Ar-C, C2 of quinoxaline), 132.30 (Ar-C), 128.58 (Ar-C), 127.73 (Ar-C), 124.77 (Ar-C), 123.75 (Ar-C), 123.36 (Ar-C), 122.73 (Ar-C), 120.93 (Ar-C), 115.21 (Ar-C), 112.18 (Ar-C), 56.17 (OCH3).; MS (m/z): 313 (M+ + 1, 8.18% %), 308.42 (100%, base beak); C16H13ClN4O (312.76).
2-Chloro-3-(2-(3,4-dimethoxybenzylidene)hydrazinyl)quinoxaline 6g
Orange powder (yield, 79%); m. p. = 164–166°C; IR (KBr, cm1): 3,292 (NH), 3,053 (CH aromatic), 2,967 (CH aliphatic); 1H NMR (400 MHz, DMSO-d6) δ 10.89 (s, 1H, NH), 8.56 (s, 1H, NCH), 7.86–7.66 (m, 2H, Ar-H, H5 and H8 of quinoxaline), 7.66–7.45 (m, 2H, Ar-H, H6 and H7 of quinoxaline), 7.39 (d, J = 11.9 Hz, 2H, Ar-H, H5 and H6 of benzene), 7.08 (d, J = 8.2 Hz, 1H, Ar-H, H2 of benzene), 3.89 (s, 3H, OCH3), 3.84 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 151.42 (Ar-C, C3 of quinoxaline), 149.66 (Ar-C, C4 and C3 of benzene), 145.24 (NCH), 144.73 (Ar-C, C2 of quinoxaline), 131.04 (Ar-C), 128.08 (Ar-C), 127.87 (Ar-C), 125.98 (Ar-C), 124.58 (Ar-C), 122.96 (Ar-C), 116.42 (Ar-C), 111.74 (Ar-C), 109.84 (Ar-C), 56.20 (OCH3), 56.07 (OCH3).; C16H13ClN4O (312.76).
2-Chloro-3-(2-(3,4,5-trimethoxybenzylidene)hydrazinyl)quinoxaline 6h
Orange powder (yield, 75%); m. p. = 210–212°C; IR (KBr, cm1): 3,275 (NH), 3,065 (CH aromatic), 2,998 (CH aliphatic), 1,607 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.2 (s, 1H, NH), 8.60 (s, 1H, NCH), 8.00–7.51 (m, 3H, Ar-H, H5, H8 and H6 of quinoxaline ), 7.52–6.94 (m, 3H, Ar-H, H7of quinoxaline, H2 and H6 of benzene), 3.89 (s, 6H, 2OCH3), 3.74 (s, 3H, OCH3); MS (m/z): 372 (M+, 100%); C18H17ClN4O3 (372.81); 13C NMR (101 MHz, DMSO-d6) δ 153.65 (Ar-C, C3 of quinoxaline), 145.03 (Ar-C, C3 and C5 of benzene), 144.82 (Ar-C, C2 of quinoxaline), 131.05 (Ar-C), 130.66 (Ar-C), 128.10 (Ar-C), 106.03 (Ar-C), 105.46 (Ar-C), 60.65 (OCH3), 56.50 (OCH3).
2-Chloro-3-(2-(3-nitrobenzylidene)hydrazinyl)quinoxaline 6i
Orange crystal (yield, 81%); m. p. = 214–216°C; IR (KBr, cm1): 3,260 (NH), 3,066 (CH aromatic), 1,615 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.47 (s, 1H, NH), 8.75 (s, 1H, Ar-H, H2 of benzene), 8.34 (d, J = 7.7 Hz, 3H, NCH, Ar-H, H5 and H6 of benzene) , 8.02 (d, J = 8.3 Hz, 1H, Ar-H, H5 of quinoxaline), 7.89 (d, J = 8.1 Hz, 1H, Ar-H, H8 of quinoxaline), 7.75 (t, J = 7.3 Hz, 1H, Ar-H, H5 of benzene), 7.68–7.20 (m, 2H, Ar-H, H6 and H7 of quinoxaline); C15H10ClN5O2 (327.73); 13C NMR (101 MHz, DMSO-d6) δ 155.63 (Ar-C, C3 of quinoxaline), 147.95 (Ar-C C3 of benzene), 145.20 (NCH), 141.46 (Ar-C, C2 of quinoxaline), 140.88 (Ar-C), 137.63 (Ar-C) , 136.73 (Ar-C), 131.22 (Ar-C), 129.70 (Ar-C), 128.15 (Ar-C), 127.06 (Ar-C), 124.60 (Ar-C), 115.95 (Ar-C).
2-Chloro-3-(2-(4-nitrobenzylidene)hydrazinyl)quinoxaline 6j
Yellow crystal (yield, 80%); m. p. = 231–232°C; IR (KBr, cm1): 3,308 )NH), 3,100 (CH aromatic), 1,615 (C=N1H NMR (400 MHz, DMSO-d6) δ 11.40 (s, 1H, NH), 8.77 (s, 1H, NCH), 8.58 (m, 1H, Ar-H, H4 of benzene), 8.37–8.26 (m, 1H, Ar-H, H5 of benzene), 8.19 (m, 1H, Ar-H, H2 of benzene), 7.96–7.85 (m, 1H, Ar-H, H6 of benzene), 7.84–7.70 (m, 2H, Ar-H, H5 and H8 of quinoxaline), 7.68–7.14 (m, 2H, Ar-H, H6 and H7 of quinoxaline); 13C NMR (101 MHz, DMSO-d6) δ 155.83 (Ar-C, C3 of quinoxaline), 148.78 (Ar-C C4 of benzene), 145.34 (NCH), 140.92 (Ar-C, C2 of quinoxaline), 137.54 (Ar-C), 136.96 (Ar-C), 134.95 (Ar-C), 133.67 (Ar-C), 131.22 (Ar-C), 128.03 (Ar-C), 127.05 (Ar-C), 124.41 (Ar-C), 121.17 (Ar-C).; C15H10ClN5O2 (327.73).
2-((2-(3-Chloroquinoxalin-2-yl)hydrazono)methyl)phenol 6k
Yellow crystal (yield, 79%); m. p. = 225–227°C; IR (KBr, cm1): 3,337 (NH), 3,055 (CH aromatic), 1,618 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.74 (s, 1H, OH), 11.51 (s, 1H, NH), 8.82 (s, 1H, NCH), 8.06–7.38 (m, 5H, Ar-H, H5, H6, H7 and H8 of quinoxaline, H6 of benzene), 7.33 (t, J = 7.7 Hz, 1H, H4 of benzene), 7.04–6.89 (m, 2H, H3 and H5 of benzene); 13C NMR (101 MHz, DMSO-d6) δ 157.92 (Ar-C, C3 of quinoxaline), 144.88 (Ar-C, C2 of benzene), 131.66 (Ar-C), 131.23 (Ar-C), 130.28 (Ar-C), 128.07 (Ar-C), 126.39 (Ar-C), 126.12 (Ar-C), 119.79 (Ar-C), 119.44 (Ar-C), 116.96 (Ar-C); C15H11ClN4O (298.73).
4-((2-(3-Chloroquinoxalin-2-yl)hydrazono)methyl)-N,N-dimethylaniline 6l
Yellow crystal (yield, 82%); m. p. = 202–204°C; IR (KBr, cm1): 3,386 (OH), 3,239 (NH), 3,050 (CH aromatic), 2,892 (CH aliphatic), 1,610 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.00 (s, 1H, NH), 8.51 (s, 1H, NCH), 7.87 (d, J = 8.6 Hz, 1H, Ar-H, H5 of quinoxaline), 7.81 (dd, J = 7.7, 5.1 Hz, 1H, Ar-H, H8 of quinoxaline), 7.60 (d, J = 8.7 Hz, 2H, Ar-H, H3 and H5 of benzene), 7.55–7.47 (m, 1H, Ar-H, H6 of quinoxaline), 7.14 (t, J = 7.6 Hz, 1H, Ar-H, H7 of quinoxaline), 6.81 (d, J = 3.0 Hz, 1H, Ar-H, H2 of benzene), 6.79 (d, J = 3.1 Hz, 1H, Ar-H, H6 of benzene), 3.03 (s, 3H, NCH3), 3.01 (s, 3H, NCH3).; C17H16ClN5 (325.80); 13C NMR (101 MHz, DMSO-d6) δ 158.56 (Ar-C, C3 of quinoxaline), 152.44 (Ar-C, C4 of benzene), 148.99 (NCH), 145.46 (Ar-C, C2 of quinoxaline), 141.38 (Ar-C), 136.73 (Ar-C), 133.09 (Ar-C), 130.66 (Ar-C), 128.87 (Ar-C), 128.08 (Ar-C), 126.61 (Ar-C), 122.37 (Ar-C), 115.44 (Ar-C), 111.97 (Ar-C), 41.30 (NCH3).
2-Chloro-3-(2-(1-phenylethylidene)hydrazinyl)quinoxaline 7a
Dark red crystal (yield, 74%); m. p. = 272–274°C; IR (KBr, cm1): 3,424 (NH), 3,020 (CH aromatic), 1,603 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.10 (s, 1H, NH), 8.29 (d, J = 29.0 Hz, 3H, Ar-H, H5 and H8 of quinoxaline, H2 of benzene), 8.02 (d, J = 8.2 Hz, 1H, Ar-H, H6 of benzene), 7.76–7.61 (m, 1H, Ar-H, H6 of quinoxaline), 7.56–7.42 (m, 1H, Ar-H, H7 of quinoxaline), 7.40–6.98 (m, 3H, Ar-H, H3 , H4 and H5 of benzene), 2.67 (s, 3H, CCH3); MS (m/z): 298 (M+ + 1, 28.08%), 297 (87.40%), 296 (100%, base beak); C16H13ClN4 (296.76); 13C NMR (101 MHz, DMSO-d6) δ 139.10 (Ar-C, C3 of quinoxaline), 132.64 (NHC), 131.50 (Ar-C, C2 of quinoxaline), 130.87 (Ar-C), 130.00 (Ar-C), 129.30 (Ar-C), 126.25 (Ar-C), 126.00 (Ar-C), 119.92 (Ar-C), 115.60 (Ar-C), 26.94 (CCH3).
2-(1-(2-(3-Chloroquinoxalin-2-yl)hydrazono)ethyl)phenol 7b
Yellow crystal (yield, 80%); m. p. = 248–250°C; 1H NMR (400 MHz, DMSO-d6) δ 13.33 (s, 1H, OH), 9.97 (s, 1H, NH), 7.69 (t, J = 7.0 Hz, 1H, Ar-H, H4 of benzene), 7.48 (d, J = 7.2 Hz, 1H, Ar-H, H5 of quinoxaline), 7.41 (d, J = 7.0 Hz, 1H, Ar-H, H8 of quinoxaline), 7.32 (dd, J = 14.3, 7.1 Hz, Ar-H, 1H, H4 of benzene), 7.23–7.12 (m, 1H, Ar-H, H6 of quinoxaline), 7.06–6.79 (m, 3H, Ar-H, H7 of quinoxaline, H6 of benzene), 2.69 (s, 3H, CCH3); 13C NMR (101 MHz, DMSO-d6) δ 166.14 (Ar-C, C3 of quinoxaline), 159.12 (Ar-C, C2 of benzene), 144.94 (NCH), 137.10 (Ar-C, C2 of quinoxaline), 133.04(Ar-C), 131.35(Ar-C), 129.52(Ar-C), 128.64 (Ar-C), 126.28 (Ar-C), 125.89 (Ar-C), 123.32 (Ar-C), 120.09 (Ar-C), 119.09 (Ar-C), 117.78 (Ar-C), 115.76 (Ar-C), 15.60 (CCH3). MS (m/z): 313 (M+ + 1, 10.2%), 312 (M+, 21.90%), 205 (100%, base beak); C16H13ClN4O (312.76).
1-((2-(3-Chloroquinoxalin-2-yl)hydrazono)methyl)naphthalen-2-ol 8
Yellow crystal (yield, 72%); m. p. = 218–220°C; IR (KBr, cm1): 3,342 (NH), 3,058 (CH aromatic), 1,656 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 13.18 (s, 1H, OH), 11.58 (s, 1H, NH), 9.76 (s, 1H, NCH), 8.20 (d, J = 8.6 Hz, 1H, Ar-H, H8 of naphthalene), 8.01–7.86 (m, 3H, Ar-H, H5 and H8 of quinoxaline, H4 of naphthalene), 7.78 (dt, J = 15.2, 8.2 Hz, 2H, Ar-H, H6 and H7 of quinoxaline), 7.65 (t, J = 7.7 Hz, 1H, Ar-H, H5 of naphthalene), 7.59 (t, J = 7.5 Hz, 1H, Ar-H, H7 of naphthalene), 7.44 (t, J = 7.5 Hz, 1H, Ar-H, H6 of naphthalene), 7.29 (d, J = 8.9 Hz, 1H, Ar-H, H3 of naphthalene); C19H13ClN4O (348.79); 13C NMR (101 MHz, DMSO-d6) δ 158.19 (Ar-C, C3 of quinoxaline), 146.46 (Ar-C, C2 of naphthalene), 144.82 (Ar-C, C2 of quinoxaline), 140.94 (NCH), 137.49 (Ar-C), 136.52 (Ar-C), 132.80 (Ar-C), 132.09 (Ar-C), 131.33(Ar-C), 129.47(Ar-C), 128.16 (Ar-C), 126.88 (Ar-C), 124.01 (Ar-C), 120.87 (Ar-C), 119.56 (Ar-C), 109.33 (Ar-C).
2-(2-(3-Chloroquinoxalin-2-yl)hydrazono)indolein-3-one 9
Red powder (yield, 89%); m. p. > 300°C; IR (KBr, cm1): 3,171 )NH), 3,062 (CH aromatic), 1715 (C=O), 1,616 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 12.35 (s, 1H, NH indole), 10.72 (s, 1H, NH), 8.52 (d, J = 7.6 Hz, 1H, Ar-H, H4 of indole), 7.86 (d, J = 8.2 Hz, 1H, Ar-H, H5 of quinoxaline), 7.71 (d, J = 7.9 Hz, 1H, Ar-H, H8 of quinoxaline), 7.58 (t, J = 7.7 Hz, 1H, Ar-H, H6 of quinoxaline), 7.37 (dq, J = 17.4, 8.0 Hz, 2H, Ar-H, H5 and H6 of indole), 7.07 (t, J = 7.5 Hz, 1H, Ar-H, H7 of quinoxaline), 6.91 (d, J = 7.7 Hz, 1H, Ar-H, H7 of indole); 13C NMR (101 MHz, DMSO-d6) δ 165.85 (Ar-C, C3 of indole), 148.07 (Ar-C, C3 of quinoxaline), 146.52 (Ar-C, C2 of indole), 144.26 (Ar-C, C2 of quinoxaline), 133.23 (Ar-C), 131.99 (Ar-C), 131.55 (Ar-C), 128.34 (Ar-C), 124.58 (Ar-C), 122.52 (Ar-C), 118.02 (Ar-C), 117.08 (Ar-C), 110.69 (Ar-C). C16H10ClN5O (323.74).
General Procedure for the Synthesis of Compounds 10a-i, 11, and 12
An equimolar of compounds 6c-h, 6j-l, 8, and 9 (0.01 mol) and sodium methoxide (0.02 mol) were refluxed in methanol (25 ml) for the appropriate time, and the reaction was followed up by TLC. The mixture was cooled then poured into water (50 ml). The formed precipitated was filtered, dried, and crystallized in ethanol to afford the corresponding compounds 10, 11, and 12, respectively.
2-(2-(4-Chlorobenzylidene)hydrazinyl)-3-methoxyquinoxaline 10a
Red orange powder (yield, 73%); m. p. = 184–186°C; 1H NMR (400 MHz, DMSO-d6) δ 11.16 (s, 1H, NH), 8.54 (s, 1H, NCH), 7.94–7.74 (m, 2H, Ar-H, H2 and H6 of benzene), 7.68 (m, Hz, 2H, 2H, Ar-H, H5 and H8 of quinoxaline), 7.53 (d, J = 8.1 Hz, 2H, 2H, Ar-H, H3 and H5 of benzene), 7.49–7.20 (m, 2H, 2H, Ar-H, H6 and H7 of quinoxaline), 4.11 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 134.33 (Ar-C), 129.34 (Ar-C), 127.19 (Ar-C), 126.49 (Ar-C), 54.44 (OCH3); C16H13ClN4O (312.76).
2-(2-(2,4-Dichlorobenzylidene)hydrazinyl)-3-methoxyquinoxaline 10b
Yellow crystal (yield, 85%); m. p. = 232–234°C 1H NMR (400 MHz, DMSO-d6) δ 11.42 (s, 1H, NH), 8.85 (s, 1H, NCH), 8.16 (d, J = 8.4 Hz, 1H, Ar-H, H5 of quinoxline), 7.71 (s, 1H, Ar-H, H3 of benzene), 7.64 (dd, J = 12.7, 8.0 Hz, 2H, Ar-H, H5 and H6 of benzene), 7.54 (d, J = 8.6 Hz, 1H, Ar-H, H8 of qunoxaline), 7.43 (t, J = 7.4 Hz, 1H, Ar-H, H6 of quinoxline), 7.34 (t, J = 6.9 Hz, 1H, Ar-H, H7 of quinoxline), 4.10 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 129.73 (Ar-C), 128.33 (Ar-C), 126.46 (Ar-C), 54.40 (OCH3); MS (m/z): 346 (M+ - 1, 5.21%), 331 (100%, base beak); C16H12Cl2N4O (347.20).
2-(2-(2,6-Dichlorobenzylidene)hydrazinyl)-3-methoxyquinoxaline 10c
Faint yellow crystal (yield, 88%); m. p. = 240–242°C; 1H NMR (400 MHz, DMSO-d6) δ 11.39 (s, 1H, NH), 8.70 (s, 1H, NCH), 7.65 (t, J = 7.9 Hz, 2H, Ar-H, H6, H7 of quinoxline), 7.59 (d, J = 8.0 Hz, 2H, Ar-H, H5 and H8 of quinoxline), 7.45 (t, J = 8.3 Hz, 2H, Ar-H, H3 and H5 of benzene), 7.39 (t, J = 7.3 Hz, 1H, Ar-H, H4 of benzene), 4.13 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 141.35 (Ar-C, C2 of quinoxline), 136.41 (NCH), 134.41 (Ar-C, C3 of quinoxline), 131.86 (Ar-C), 131.36 (Ar-C), 129.45 (Ar-C), 128.06 (Ar-C), 127.20 (Ar-C), 126.50 (Ar-C), 125.27 (Ar-C), 54.49 (OCH3); C16H12Cl2N4O (347.20).
2-Methoxy-3-(2-(2-methoxybenzylidene)hydrazinyl)quinoxaline 10d
Orange yellow powder (yield, 89%); m. p. = 225–227°C; IR (KBr, cm1): 3,449 (NH), 3,021 (CH aromatic), 2,935 (CH aliphatic), 1,656 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 9.04 (s, 1H, NH), 8.66 (s, 1H, NCH), 8.12 (d, J = 7.1 Hz, 1H, Ar-H, H5 of quinoxline), 7.95 (d, J = 8.1 Hz, 1H, Ar-H, H8 of quinoxline), 7.84 (d, J = 7.3 Hz, 1H, Ar-H, H3 of benzene), 7.73 (t, J = 7.5 Hz, 1H, Ar-H, H6 of quinoxline), 7.55 (d, J = 7.3 Hz, 1H, Ar-H, H6 of benzene), 7.45 (t, J = 7.6 Hz, 1H, Ar-H, H7 of quinoxline), 7.14 (d, J = 8.3 Hz, 1H, Ar-H, H5 of benzene), 7.08 (t, J = 7.4 Hz, 1H, Ar-H, H4 of benzene), 3.93 (s, 3H, OCH3 of quinoxaline), 3.79 (s, 3H, OCH3 of benzene); MS (m/z): 308 (M+, 26.23%), 91 (100%, base beak); C17H16N4O2 (308.34); 13C NMR (101 MHz, DMSO-d6) δ 158.05 (Ar-C, C2 of quinoxline), 137.69 (NCH), 131.39 (Ar-C, C3 of quinoxline), 129.43 (Ar-C), 126.32 (Ar-C), 125.97 (Ar-C), 123.03 (Ar-C), 122.62 (Ar-C), 121.24 (Ar-C), 112.32 (Ar-C), 56.22 (OCH3).
2-(2-(3,4-Dimethoxybenzylidene)hydrazinyl)-3-methoxyquinoxaline 10e
Faint yellow crystal (yield, 81%); m. p. = 235–237°C; IR (KBr, cm1): 3,419 (NH), 3,060 (CH aromatic), 2,941 (CH aliphatic), 1,661 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 12.95 (s, 1H, NH), 8.65 (s, 1H, NCH), 8.01 (t, J = 7.6 Hz, 1H, Ar-H, H5 of quinoxaline), 7.85–7.69 (m, 2H, Ar-H, H6 and H8 of quinoxaline), 7.54 (t, J = 7.5 Hz, 1H, Ar-H, H7 of quinoxaline), 7.46 (t, J = 8.0 Hz, 1H, Ar-H, H6 and H2 of benzene), 7.38 (d, J = 7.7 Hz, 1H, Ar-H, H6 of benzene), 7.10 (t, J = 8.8 Hz, 1H, Ar-H, H5 of benzene), 4.17 (s, 3H, OCH3), 3.92 (s, 3H, OCH3 of benzene), 3.85 (s, 3H, OCH3 of benzene); 13C NMR (101 MHz, DMSO-d6) δ 149.69 (Ar-C, C2 of quinoxline), 149.60 (NCH), 145.19 (Ar-C, C3 of quinoxline), 128.02 (Ar-C), 127.59 (Ar-C), 127.25 (Ar-C), 126.48 (Ar-C), 126.03 (Ar-C), 124.52 (Ar-C), 123.92 (Ar-C), 123.75 (Ar-C), 111.88 (Ar-C), 110.33 (Ar-C), 109.96 (Ar-C), 56.35 (OCH3), 56.26 (OCH3), 56.13 (OCH3); C18H18N4O3 (338.37).
2-Methoxy-3-(2-(3,4,5-trimethoxybenzylidene)hydrazinyl)quinoxaline 10f
Faint yellow crystal (yield, 88%); m. p. = 256–258°C; IR (KBr, cm1): 3,421 (NH), 3,025 (CH aromatic), 2,940(CH aliphatic), 1,657 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 10.96 (s, 1H, NH), 8.43 (s, 1H, NCH), 7.67–7.50 (m, 2H, Ar-H, H5 and H8 of quinoxaline), 7.38 (t, J = 6.4 Hz, 1H, Ar-H, H6 of quinoxaline), 7.31–7.20 (m, 1H, Ar-H, H7 of quinoxaline), 7.09 (s, 2H, Ar-H, H2 and H6 of benzene), 4.07 (s, 3H, OCH3 of quinoxaline), 3.88 (s, 6H, OCH3 of benzene), 3.72 (s, 3H. OCH3); 13C NMR (101 MHz, DMSO-d6) δ 153.62 (Ar-C, C2 of quinoxline), 126.7 (Ar-C), 126.33 (Ar-C), 106.45 (Ar-C), 104.72 (Ar-C), 60.62 (OCH3), 56.45 (OCH3). 54.14 (OCH3); C19H20N4O4 (368.39).
2-Methoxy-3-(2-(4-nitrobenzylidene)hydrazinyl)quinoxaline 10g
Yellow crystal (yield, 84%); m. p. = 228–230°C; IR (KBr, cm1): 3,310 )NH), 3,079 (CH aromatic), 2,991 (CH aliphatic), 1,623 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 11.38 (s, 1H, NH), 8.68 (s, 1H, NCH), 8.58 (d, J = 8.0 Hz, 2H, Ar-H, H3 and H5 of benzene), 8.27 (d, J = 8.0 Hz, 2H, Ar-H, H2 and H6 of benzene), 7.88–7.73 (m, 2H, Ar-H, H5 and H8 of quinoxaline), 7.45 (d, J = 24.2 Hz, 2H, Ar-H, H6 and H7 of quinoxaline), 4.14 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 148.77 (Ar-C, C4 of benzene), 143.70 (Ar-C, C2 of quinoxline), 137.24 (NCH), 133.47 (Ar-C), 130.92 (Ar-C), 127.28 (Ar-C), 126.54 (Ar-C), 124.66 (Ar-C), 123.63 (Ar-C), 123.33 (Ar-C), 54.53 (OCH3). C16H13N5O3 (323.31).
2-((2-(3-Methoxyquinoxalin-2-yl)hydrazono)methyl)phenol 10h
Yellow crystal (yield, 81%); m. p. = 247–249°C; 1H NMR (400 MHz, DMSO-d6) δ 12.91 (s, 1H, OH), 11.43 (s, 1H, NH), 8.89 (s, 1H, NCH), 7.83 (d, J = 7.8 Hz, 1H, Ar-H, H5 of quinoxaline), 7.76 (d, J = 6.7 Hz, 1H, Ar-H, H8 of quinoxaline), 7.71 (d, J = 7.8 Hz, 1H, Ar-H, H3 of benzene), 7.53 (t, J = 7.4 Hz, 1H, Ar-H, H6 of quinoxaline), 7.45 (t, J = 7.3 Hz, 1H, Ar-H, H7 of quinoxaline), 7.34 (q, J = 8.3, 7.2 Hz, 1H, Ar-H, H6 of benzene), 7.03 (t, J = 7.8 Hz, 1H, Ar-H, H4 of benzene), 6.94 (t, J = 7.0 Hz, 1H, Ar-H, H5 of benzene), 4.15 (s, 3H, OCH3); C16H14N4O2 (294.31).
4-((2-(3-Methoxyquinoxalin-2-yl)hydrazono)methyl)-N,N-dimethylaniline 10i
Yellow crystal (yield, 86%); m. p. = 214–216°C; IR (KBr, cm1): 3,185 (NH), 3,042 (CH aromatic), 2,941 (CH aliphatic), 1,616 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 10.77 (s, 1H, NH), 8.41 (s, 1H, NCH), 7.70–7.52 (m, 4H, Ar-H, H5 and H8 of quinoxaline, H2 and H6 of benzene), 7.40 (s, 1H. Ar-H, H6 of quinoxaline), 7.28 (s, 1H, Ar-H, H7 of quinoxaline), 6.79 (d, J = 8.4 Hz, 2H, H3 and H5 of benzene), 4.08 (s, 3H, OCH3), 3.00 (s, 6H, NCH3); 13C NMR (101 MHz, DMSO-d6) δ 152.14 (Ar-C, C2 of quinoxaline), 141.00 (Ar-C, C4 of benzene), 133.22 (Ar-C), 129.49 (Ar-C), 128.06 (Ar-C), 127.03 (Ar-C), 126.39 (Ar-C), 122.77 (Ar-C), 112.27 (Ar-C), 54.35 (OCH3), 41.57 (NH3); C18H19N5O (321.38).
1-((2-(3-Methoxyquinoxalin-2-yl)hydrazono)methyl)naphthalen-2-ol 11
Orange powder (yield, 89%); m. p. = 153–155°C; IR (KBr, cm1): 3,423 (NH), 3,053 (CH aromatic), 2,941 (CH aliphatic), 1,663 (C=N); 1H NMR (400 MHz, DMSO-d6) δ 12.71 (s, 1H, OH), 9.67 (s, 1H, NH), 8.29 (s, 1H, NCH), 8.05–7.77 (m, 3H, Ar-H, H5 and H8 of quinoxaline, H3 of naphthalene), 7.73 (d, J = 7.6 Hz, 1H, Ar-H, H3 of naphthalene), 7.64 (t, J = 7.1 Hz, 1H, Ar-H, H6 of quinoxaline), 7.53 (t, J = 7.4 Hz, 1H, Ar-H, H7 of quinoxaline), 7.45 (q, J = 7.6 Hz, 2H, Ar-H, H5 and H8 of quinoxaline), 7.29 (d, J = 11.3 Hz, 2H, Ar-H, H6 and H7 of quinoxaline), 4.19 (s, 3H, OCH3); 13C NMR (101 MHz, DMSO-d6) δ 158.18 (Ar-C, C3 of naphthalene), 140.39 (Ar-C, C2 of qunoxaline) , 132.03 (Ar-C, C3 of qunoxaline), 129.44, 128.44 (Ar-C), 128.19 (Ar-C), 127.54 (Ar-C), 126.68 (Ar-C), 125.87 (Ar-C), 124.34 (Ar-C), 124.00 (Ar-C), 121.24 (Ar-C), 119.47 (Ar-C), 116.00 (Ar-C), 109.52 (Ar-C), 54.74 (OCH3).; MS (m/z): 345 (M+ +1, 15.61%), 344 (M+, 21.63%), 332 (100%, base beak); C20H16N4O2 (344.37).
2-(2-(3-Methoxyquinoxalin-2-yl)hydrazono)indolein-3-one 12
Yellow crystal (yield, 87%); m. p. = 242–244°C; 1H NMR (400 MHz, DMSO-d6) δ 12.17 (s, 1H, NH of indole), 10.75 (s, 1H, NH of quinoxaline), 8.49 (d, J = 7.2 Hz, 1H, Ar-H, H5 of quinoxaline), 7.85 (d, J = 8.1 Hz, 1H, Ar-H, H8 of quinoxaline), 7.70 (d, J = 7.9 Hz, 1H, Ar-H, H4 of indole), 7.58 (t, J = 7.3 Hz, 1H, Ar-H, H6 of quinoxaline), 7.35 (q, J = 7.6 Hz, 2H, Ar-H, H7 of quinoxaline, H7 of indole), 7.06 (d, J = 7.4 Hz, 1H, Ar-H, H5 of indole), 6.91 (d, J = 7.7 Hz, 1H, Ar-H, H6 of indole), 4.13 (s, 3H, OCH3); C17H13N5O2 (319.32).
Biological Evaluation
Experimental protocols applied for our newly synthesized compounds in the different biological assays were provided in detail in the Supplementary Material.
Anti-proliferative Activities Against Human Liver Cancer Cell Lines
The new quinoxaline derivatives of the five series: 6, 7, 8, 9, and 10 were evaluated for their potential anti-proliferative activity against two liver cancer cell lines (HepG-2 and HuH-7) obtained from the American Type Culture Collection. Cytotoxicity was assessed following the SRB colorimetric assay protocol (Skehan et al., 1990), as reported earlier (Eldehna et al., 2016; Sabt et al., 2018; Ghanem et al., 2020).
Histone Deacetylase Inhibitory Activities
The new quinoxaline analogs (6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12) were further assayed for their deacetylase enzymes inhibitory activities (HDAC1, HDAC4, and HDAC6 subtypes) based on a homogeneous fluorescence release assay, as discussed before (Fournel et al., 2008).
Cell Cycle Analysis
The aforementioned most active compound 6c incorporating 4-Cl side chain was further evaluated through cell cycle analysis in HepG-2 cell line at IC50 = 1.53 µM, using BD FACS Calibur flow cytometer, as described previously (Eliaa et al., 2020; Sabt et al., 2020).
Annexin V-FITC Apoptosis Assay
Furthermore, our most active candidate 6c was assayed for apoptosis induction using the FITC Annexin-V/PI kit (Becton Dickenson, Franklin Lakes, NJ) following the manufacture’s protocol. The previous compound was analyzed by FACS as we previously described (Al-Rashood et al., 2020; Eliaa et al., 2020; Eldehna et al., 2021).
Apoptotic and Anti-Apoptotic Marker Levels (Bcl-2 and Bax)
Quantitative real-time PCR to evaluate the effects of compounds 6c, 6d, 6f, 6g, 6k, 6l, 7b, 8, 10h, and 12 on two important target genes (Bcl-2 and Bax) and the housekeeping gene (GAPDH) in HepG-2 cells was performed as well. The method was performed in detail as previously explained (Eliaa et al., 2020).
In silico Studies
Docking Studies
A molecular docking study of the newly synthesized quinoxaline derivatives at the histone deacetylase (HDAC) receptor was performed, and the co-crystallized inhibitor, trichostatin A (TSA), was used as a reference standard. Using MOE 2019.0102 drug design software (Inc, 2016), the binding mode of the compound against histone deacetylase (ID: 1C3R) was predicted (Finnin et al., 1999). The crystal structure of the target receptor (HDAC) was downloaded from Protein Data Bank (http://www.rcsb.org/, PDB code: 1C3R, resolution of 2.00 Å) (Finnin et al., 1999). The protein structure was prepared for docking studies by the default method (Alnajjar et al., 2020; Soltan et al., 2021; Soltane et al., 2021; Zaki et al., 2020). The deacetylase and deacetylase-TSA structures show an active site consisting of a tubular pocket, a zinc-binding site (which is the metal cofactor required for HDAC activity), and two Asp-His charge-relay systems, and explain the mechanism of HDAC inhibition (Finnin et al., 1999). Validation of the docking procedure was carried out by applying the docking process for the co-crystallized ligand (Elmaaty et al., 2021a; Elmaaty et al., 2021b; Kandeil et al., 2021). All of the newly synthesized quinoxaline derivatives were prepared and imported in the same database together with the co-crystallized inhibitor (TSA) and generally docked. After completion of the docking process, the obtained poses for each were carefully studied, and the ones having the best scores and binding modes with the protein pocket residues were selected.
ADMET Studies
ADMET descriptors (absorption, distribution, metabolism, excretion, and toxicity) of the synthesized compounds were determined using Discovery studio 4.0. i) Blood-brain barrier penetration predicts blood-brain barrier penetration of a molecule. ii) Intestinal absorption predicts human intestinal absorption (HIA) after oral administration. iii) Aqueous solubility predicts the solubility of each compound in the water at 25°C. iv) CYP2D6 binding predicts cytochrome P450 2D6 enzyme inhibition. v) Plasma protein binding predicts the fraction of drug bound to plasma proteins in the blood (Van De Waterbeemd and Gifford, 2003). Discovery studio 4.0 was used to predict ADMET descriptors for all compounds. At first, the CHARMM force field was applied, and then, the compounds were prepared and minimized according to the preparation of small molecule protocol (Al-Karmalawy et al., 2021; Al-Karmalawy and Eissa, 2021). Then, ADMET descriptors protocol was applied to carry out these studies (Ibrahim et al., 2017; El-Gamal et al., 2018; El-Zahabi et al., 2019; El-Shershaby et al., 2021a).
Toxicity Studies
The toxicity parameters of the synthesized compounds were calculated using Discovery studio 4.0. At first, the CHARMM force field was applied, and then, the compounds were prepared and minimized according to the preparation of small molecule protocol. Then, different parameters were calculated from the toxicity prediction (extensible) protocol as follows: i) FDA rodent carcinogenicity which computes the probability of a submitted chemical structure being a carcinogen, ii) carcinogenic potency TD50 which predicts the tumorigenic dose rate 50 (TD50) of a chemical in a rodent chronic exposure toxicity test of carcinogenic potency (Venkatapathy et al., 2009), iii) rat maximum tolerated dose which predicts the rat maximum tolerated dose (MTD) of a chemical (4,5) (Goodrnan and Wilson, 1992), iv) developmental toxicity potential which predicts whether a particular compound is likely to be toxic in a developmental toxicity potential assessment (Agency, 1991; Louisse et al., 2015), v) rat oral LD50 which predicts the rat oral acute median lethal dose (LD50) in the toxicity test of a chemical (Gonella Diaza et al., 2015), vi) rat chronic LOAEL which predicts the rat chronic lowest observed adverse effect level (LOAEL) value of a chemical (Venkatapathy et al., 2004; Benfenati, 2016), vii) ocular irritancy which predicts whether a particular compound is likely to be an ocular irritant and how severe the irritation is in the Draize test (Macfarlane et al., 2009), viii) skin irritancy predicts whether a particular compound is likely to be a skin irritant and how severe it is in a rabbit skin irritancy test (Macfarlane et al., 2009).
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding authors.
Author Contributions
CM: Conceptualization; MT and AB: Synthesis; AM and AN: Biology; NM: ADMET studies; IE: Supervision, interpretation, writing, and revising the manuscript. AA-K: Docking, interpretation, writing, and revising the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
MT, AM, AN, IE, and AA-K acknowledge Al-Azhar University, Cairo, Egypt, and Horus University-Egypt, New Damietta, Egypt for partial financial support. AB is thankful to Taif University Researchers Supporting Project number (TURSP-2020/35), Taif University, Taif, Saudi Arabia.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.725135/full#supplementary-material
References
Abbass, E. M., Khalil, A. K., Mohamed, M. M., Eissa, I. H., and El-Naggar, A. M. (2020). Design, Efficient Synthesis, Docking Studies, and Anticancer Evaluation of New Quinoxalines as Potential Intercalative Topo II Inhibitors and Apoptosis Inducers. Bioorg. Chem. 104, 104255. doi:10.1016/j.bioorg.2020.104255
Abo Elmaaty, A., Hamed, M. I. A., Ismail, M. I., B. Elkaeed, E., S. Abulkhair, H., Khattab, M., et al. (2021). Computational Insights on the Potential of Some NSAIDs for Treating COVID-19: Priority Set and Lead Optimization. Molecules 26 (12), 3772. doi:10.3390/molecules26123772
Agency, U. E. P. (1991). Guidelines for Developmental Toxicity Risk Assessment. Fed. Reg. 56 (234), 63798–63826.
Agoni, C., Olotu, F. A., Ramharack, P., and Soliman, M. E. (2020). Druggability and Drug-Likeness Concepts in Drug Design: Are Biomodelling and Predictive Tools Having Their Say? J. Mol. Model. 26 (6), 1–11. doi:10.1007/s00894-020-04385-6
Al-Karmalawy, A. A., and Elshal, M. F. (2021). Concanavalin-A Shows Synergistic Cytotoxicity with Tamoxifen via Inducing Apoptosis in Estrogen Receptor-Positive Breast Cancer: In Vitro and Molecular Docking Studies Pharm Sci. doi:10.34172/PS.2021.22
Al‐Karmalawy, A. A., and Eissa, I. H. J. P. S. (2021). Molecular Docking and Dynamics Simulations Reveal the Potential of Anti-HCV Drugs to Inhibit COVID-19 Main Protease Pharm Sci. doi:10.34172/PS.2021.3
Al-Karmalawy, A. A., Dahab, M. A., Metwaly, A. M., Elhady, S. S., Elkaeed, E. B., Eissa, I. H., et al. (2021). Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of ACEIs against SARS-CoV-2 Targeting the hACE2 Receptor. Front. Chem. 9 (227). doi:10.3389/fchem.2021.661230
Al-Karmalawy, A. A., and Khattab, M. (2020). Molecular Modelling of Mebendazole Polymorphs as a Potential Colchicine Binding Site Inhibitor. New J. Chem. 44 (33), 13990–13996. doi:10.1039/d0nj02844d
Al-Rashood, S. T., Hamed, A. R., Hassan, G. S., Alkahtani, H. M., Almehizia, A. A., Alharbi, A., et al. (2020). Antitumor Properties of Certain Spirooxindoles towards Hepatocellular Carcinoma Endowed with Antioxidant Activity. J. Enzyme Inhib. Med. Chem. 35 (1), 831–839. doi:10.1080/14756366.2020.1743281
Alesawy, M. S., Al‐Karmalawy, A. A., Elkaeed, E. B., Alswah, M., Belal, A., Taghour, M. S., et al. (2020). Design and Discovery of New 1, 2, 4‐triazolo [4, 3‐c] Quinazolines as Potential DNA Intercalators and Topoisomerase II Inhibitors, e2000237.
Alnajjar, R., Mostafa, A., Kandeil, A., and Al-Karmalawy, A. A. J. H. (2020). Molecular Docking, Molecular Dynamics, and In Vitro Studies Reveal the Potential of Angiotensin II Receptor Blockers to Inhibit the COVID-19 Main Protease. 6(12), e05641. doi:10.1016/j.heliyon.2020.e05641
Banu, S., Bollu, R., Bantu, R., Nagarapu, L., Polepalli, S., Jain, N., et al. (2017). Design, Synthesis and Docking Studies of Novel 1,2-Dihydro-4-Hydroxy-2-Oxoquinoline-3-Carboxamide Derivatives as a Potential Anti-proliferative Agents. Eur. J. Med. Chem. 125, 400–410. doi:10.1016/j.ejmech.2016.09.062
Brück, W., and Zamvil, S. S. (2012). Laquinimod, a Once-Daily Oral Drug in Development for the Treatment of Relapsing-Remitting Multiple Sclerosis. Expert Rev. Clin. Pharmacol. 5 (3), 245–256. doi:10.1586/ecp.12.12
Cappellacci, L., Perinelli, D. R., Maggi, F., Grifantini, M., and Petrelli, R. (2020). Recent Progress in Histone Deacetylase Inhibitors as Anticancer Agents. Cmc 27 (15), 2449–2493. doi:10.2174/0929867325666181016163110
Corona, P., Carta, A., Loriga, M., Vitale, G., and Paglietti, G. (2009). Synthesis and In Vitro Antitumor Activity of New Quinoxaline Derivatives. Eur. J. Med. Chem. 44 (4), 1579–1591. doi:10.1016/j.ejmech.2008.07.025
Davis, I. W., and Baker, D. (2009). RosettaLigand Docking with Full Ligand and Receptor Flexibility. J. Mol. Biol. 385 (2), 381–392. doi:10.1016/j.jmb.2008.11.010
Deng, B., Luo, Q., Halim, A., Liu, Q., Zhang, B., and Song, G. (2020). The Antiangiogenesis Role of Histone Deacetylase Inhibitors: Their Potential Application to Tumor Therapy and Tissue Repair. DNA Cel. Biol. 39 (2), 167–176. doi:10.1089/dna.2019.4877
Eissa, I. H., El-Naggar, A. M., and El-Hashash, M. A. (2016). Design, Synthesis, Molecular Modeling and Biological Evaluation of Novel 1H-Pyrazolo[3,4-B]pyridine Derivatives as Potential Anticancer Agents. Bioorg. Chem. 67, 43–56. doi:10.1016/j.bioorg.2016.05.006
Eissa, I. H., El-Naggar, A. M., El-Sattar, N. E. A. A., and Youssef, A. S. A. (2018). Design and Discovery of Novel Quinoxaline Derivatives as Dual DNA Intercalators and Topoisomerase II Inhibitors. Acamc 18 (2), 195–209. doi:10.2174/1871520617666170710182405
El-Helby, A. G. A., Sakr, H., Eissa, I. H., Al‐Karmalawy, A. A., and El‐Adl, K. (2019b). Benzoxazole/benzothiazole‐derived VEGFR‐2 Inhibitors: Design, Synthesis, Molecular Docking, and Anticancer Evaluations. Archiv Der Pharmazie 352 (12), 1900178.
El-Gamal, K. M., El-Morsy, A. M., Saad, A. M., Eissa, I. H., and Alswah, M. (2018). Synthesis, Docking, QSAR, ADMET and Antimicrobial Evaluation of New Quinoline-3-Carbonitrile Derivatives as Potential DNA-Gyrase Inhibitors. J. Mol. Struct. 1166, 15–33. doi:10.1016/j.molstruc.2018.04.010
El-Helby, A. A., Sakr, H., Eissa, I. H., Abulkhair, H., Al-Karmalawy, A. A., and El-Adl, K. (2019a). Design, Synthesis, Molecular Docking, and Anticancer Activity of Benzoxazole Derivatives as VEGFR-2 Inhibitors. Arch. Pharm. (Weinheim) 352 (10), e1900113. doi:10.1002/ardp.201900113
El-Naggar, A. M., Abou-El-Regal, M. M., El-Metwally, S. A., Sherbiny, F. F., and Eissa, I. H. (2017). Synthesis, Characterization and Molecular Docking Studies of Thiouracil Derivatives as Potent Thymidylate Synthase Inhibitors and Potential Anticancer Agents. Mol. Divers. 21 (4), 967–983. doi:10.1007/s11030-017-9776-1
El-Shershaby, M. H., El-Gamal, K. M., Bayoumi, A. H., El-Adl, K., Alswah, M., Ahmed, H. E. A., et al. (2021a). The Antimicrobial Potential and Pharmacokinetic Profiles of Novel Quinoline-Based Scaffolds: Synthesis and In Silico Mechanistic Studies as Dual DNA Gyrase and DHFR Inhibitors. New J. Chem. doi:10.1039/D1NJ02838C
El-Shershaby, M. H., Ghiaty, A., Bayoumi, A. H., Al-Karmalawy, A. A., Husseiny, E. M., El-Zoghbi, M. S., et al. (2021b). From Triazolophthalazines to Triazoloquinazolines: A Bioisosterism-Guided Approach toward the Identification of Novel PCAF Inhibitors with Potential Anticancer Activity. Bioorg. Med. Chem. 42, 116266. doi:10.1016/j.bmc.2021.116266
El-Zahabi, M. A., Elbendary, E. R., Bamanie, F. H., Radwan, M. F., Ghareib, S. A., and Eissa, I. H. (2019). Design, Synthesis, Molecular Modeling and Anti-hyperglycemic Evaluation of Phthalimide-Sulfonylurea Hybrids as PPARγ and SUR Agonists. Bioorg. Chem. 91, 103115. doi:10.1016/j.bioorg.2019.103115
El-Zahabi, M. A., Sakr, H., El-Adl, K., Zayed, M., Abdelraheem, A. S., Eissa, S. I., et al. (2020). Design, Synthesis, and Biological Evaluation of New Challenging Thalidomide Analogs as Potential Anticancer Immunomodulatory Agents. Bioorg. Chem. 104, 104218. doi:10.1016/j.bioorg.2020.104218
Eldehna, W., Fares, M., Ibrahim, H., Alsherbiny, M., Aly, M., Ghabbour, H., et al. (2016). Synthesis and Cytotoxic Activity of Biphenylurea Derivatives Containing Indolin-2-One Moieties. Molecules 21 (6), 762. doi:10.3390/molecules21060762
Eldehna, W. M., Abo-Ashour, M. F., Nocentini, A., Gratteri, P., Eissa, I. H., Fares, M., et al. (2017). Novel 4/3-((4-Oxo-5-(2-Oxoindolin-3-Ylidene)thiazolidin-2-Ylidene)amino) Benzenesulfonamides: Synthesis, Carbonic Anhydrase Inhibitory Activity, Anticancer Activity and Molecular Modelling Studies. Eur. J. Med. Chem. 139, 250–262. doi:10.1016/j.ejmech.2017.07.073
Eldehna, W. M., Al-Rashood, S. T., Al-Warhi, T., Eskandrani, R. O., Alharbi, A., and El Kerdawy, A. M. (2021). Novel Oxindole/benzofuran Hybrids as Potential Dual CDK2/GSK-3β Inhibitors Targeting Breast Cancer: Design, Synthesis, Biological Evaluation, and In Silico Studies. J. Enzyme Inhib. Med. Chem. 36 (1), 270–285. doi:10.1080/14756366.2020.1862101
Eliaa, S. G., Al-Karmalawy, A. A., Saleh, R. M., and Elshal, M. F. (2020). Empagliflozin and Doxorubicin Synergistically Inhibit the Survival of Triple-Negative Breast Cancer Cells via Interfering with the mTOR Pathway and Inhibition of Calmodulin: In Vitro and Molecular Docking Studies. ACS Pharmacol. Transl. Sci. 3 (6), 1330–1338. doi:10.1021/acsptsci.0c00144
Ellis, L., Hammers, H., and Pili, R. (2009). Targeting Tumor Angiogenesis with Histone Deacetylase Inhibitors. Cancer Lett. 280 (2), 145–153. doi:10.1016/j.canlet.2008.11.012
Elmaaty, A. A., Alnajjar, R., Hamed, M. I. A., Khattab, M., Khalifa, M. M., and Al-Karmalawy, A. A. (2021a). Revisiting Activity of Some Glucocorticoids as a Potential Inhibitor of SARS-CoV-2 Main Protease: Theoretical Study. RSC Adv. 11 (17), 10027–10042. doi:10.1039/d0ra10674g
Elmaaty, A. A., Darwish, K. M., Khattab, M., Elhady, S. S., Salah, M., Hamed, M. I. A., et al. (2021b). In a Search for Potential Drug Candidates for Combating COVID-19: Computational Study Revealed Salvianolic Acid B as a Potential Therapeutic Targeting 3CLpro and Spike Proteins. J. Biomol. Struct. Dyn., 1–28. doi:10.1080/07391102.2021.1918256
Elmetwally, S. A., Saied, K. F., Eissa, I. H., and Elkaeed, E. B. (2019). Design, Synthesis and Anticancer Evaluation of Thieno[2,3-D]pyrimidine Derivatives as Dual EGFR/HER2 Inhibitors and Apoptosis Inducers. Bioorg. Chem. 88, 102944. doi:10.1016/j.bioorg.2019.102944
Finnin, M. S., Donigian, J. R., Cohen, A., Richon, V. M., Rifkind, R. A., Marks, P. A., et al. (1999). Structures of a Histone Deacetylase Homologue Bound to the TSA and SAHA Inhibitors. Nature 401 (6749), 188–193. doi:10.1038/43710
Fournel, M., Bonfils, C., Hou, Y., Yan, P. T., Trachy-Bourget, M.-C., Kalita, A., et al. (2008). MGCD0103, a Novel Isotype-Selective Histone Deacetylase Inhibitor, Has Broad Spectrum Antitumor Activity In Vitro and In Vivo. Mol. Cancer Ther. 7 (4), 759–768. doi:10.1158/1535-7163.mct-07-2026
Gaber, A. A., Bayoumi, A. H., El-Morsy, A. M., Sherbiny, F. F., Mehany, A. B. M., and Eissa, I. H. (2018). Design, Synthesis and Anticancer Evaluation of 1H-Pyrazolo[3,4-D]pyrimidine Derivatives as Potent EGFRWT and EGFRT790M Inhibitors and Apoptosis Inducers. Bioorg. Chem. 80, 375–395. doi:10.1016/j.bioorg.2018.06.017
Ghafourian, T., and Amin, Z. (2013). QSAR Models for the Prediction of Plasma Protein Binding. Bioimpacts 3 (1), 21–27. doi:10.5681/bi.2013.011
Ghanem, A., Emara, H. A., Muawia, S., Abd El Maksoud, A. I., Al-Karmalawy, A. A., and Elshal, M. F. (2020). Tanshinone IIA Synergistically Enhances the Antitumor Activity of Doxorubicin by Interfering with the PI3K/AKT/mTOR Pathway and Inhibition of Topoisomerase II: In Vitro and Molecular Docking Studies. New J. Chem. 44 (40), 17374–17381. doi:10.1039/d0nj04088f
Gonella Diaza, R., Manganelli, S., Esposito, A., Roncaglioni, A., Manganaro, A., and Benfenati, E. (2015). Comparison Ofin Silicotools for Evaluating Rat Oral Acute Toxicity. SAR QSAR Environ. Res. 26 (1), 1–27. doi:10.1080/1062936x.2014.977819
Goodrnan, G., and Wilson, R. (1992). Comparison of the Dependence of the TD50 on Maximum Tolerated Dose for Mutagens and Nonmutagens. Risk Anal. 12 (4), 525–533. doi:10.1111/j.1539-6924.1992.tb00709.x
Gross, D. J., Weiss, L., Reibstein, I., Hedlund, G., Dahlén, E., Rapoport, M. J., et al. (2001). The Immunomodulator Linomide: Role in Treatment and Prevention of Autoimmune Diabetes Mellitus. Int. immunopharmacology 1 (6), 1131–1139. doi:10.1016/s1567-5769(01)00042-x
Gupta, N., Al Ustwani, O., Shen, L., and Pili, R. (2014). Mechanism of Action and Clinical Activity of Tasquinimod in Castrate-Resistant Prostate Cancer. Onco Targets Ther. 7, 223–234. doi:10.2147/OTT.S53524
Hou, T., Zhu, L., Chen, L., and Xu, X. (2003). Mapping the Binding Site of a Large Set of Quinazoline Type EGF-R Inhibitors Using Molecular Field Analyses and Molecular Docking Studies. J. Chem. Inf. Comput. Sci. 43 (1), 273–287. doi:10.1021/ci025552a
Ibrahim, M.-K., Abd-Elrahman, A. A., Ayyad, R. R. A., El-Adl, K., Mansour, A. M., and Eissa, I. H. (2013). Design and Synthesis of Some Novel 2-(3-Methyl-2-Oxoquinoxalin-1(2h)-Yl)-N-(4-(substituted)phenyl)acetamide Derivatives for Biological Evaluation as Anticonvulsant Agents. Bull. Fac. Pharm. Cairo Univ. 51 (1), 101–111. doi:10.1016/j.bfopcu.2012.11.003
Ibrahim, M.-K., El-Adl, K., and Al-Karmalawy, A. A. (2015). Design, Synthesis, Molecular Docking and Anticonvulsant Evaluation of Novel 6-Iodo-2-Phenyl-3-Substituted-Quinazolin-4(3h)-Ones. Bull. Fac. Pharm. Cairo Univ. 53 (2), 101–116. doi:10.1016/j.bfopcu.2015.05.001
Ibrahim, M. (2012). Design, Synthesis, Molecular Docking and Biological Evaluation of Some Novel Quinazolin-4(3h)-One Derivatives as Anti-inflammatory Agents. Al-Azhar J. Pharm. Sci. 46 (2), 185–203. doi:10.21608/ajps.2012.7145
Ibrahim, M. K., Eissa, I. H., Alesawy, M. S., Metwaly, A. M., Radwan, M. M., and ElSohly, M. A. (2017). Design, Synthesis, Molecular Modeling and Anti-hyperglycemic Evaluation of Quinazolin-4(3h)-One Derivatives as Potential PPARγ and SUR Agonists. Bioorg. Med. Chem. 25 (17), 4723–4744. doi:10.1016/j.bmc.2017.07.015
Ibrahim, M. K., Taghour, M. S., Metwaly, A. M., Belal, A., Mehany, A. B. M., Elhendawy, M. A., et al. (2018). Design, Synthesis, Molecular Modeling and Anti-proliferative Evaluation of Novel Quinoxaline Derivatives as Potential DNA Intercalators and Topoisomerase II Inhibitors. Eur. J. Med. Chem. 155, 117–134. doi:10.1016/j.ejmech.2018.06.004
Inc, C. C. G. (2016). Molecular Operating Environment (MOE). Chemical Computing Group Inc. 1010 Sherbooke St. West, Suite# 910, Montreal.
Isaacs, J. T., Antony, L., Dalrymple, S. L., Brennen, W. N., Gerber, S., Hammers, H., et al. (2013). Tasquinimod Is an Allosteric Modulator of HDAC4 Survival Signaling within the Compromised Cancer Microenvironment. Cancer Res. 73 (4), 1386–1399. doi:10.1158/0008-5472.can-12-2730
Isaacs, J. T., Pili, R., Qian, D. Z., Dalrymple, S. L., Garrison, J. B., Kyprianou, N., et al. (2006). Identification of ABR-215050 as lead Second Generation Quinoline-3-Carboxamide Anti-angiogenic Agent for the Treatment of Prostate Cancer. Prostate 66 (16), 1768–1778. doi:10.1002/pros.20509
Jennbacken, K., Welén, K., Olsson, A., Axelsson, B., Törngren, M., Damber, J.-E., et al. (2012). Inhibition of Metastasis in a Castration Resistant Prostate Cancer Model by the Quinoline-3-Carboxamide Tasquinimod (ABR-215050). Prostate 72 (8), 913–924. doi:10.1002/pros.21495
Kandeil, A., Mostafa, A., Kutkat, O., Moatasim, Y., Al-Karmalawy, A. A., Rashad, A. A., et al. (2021). Bioactive Polyphenolic Compounds Showing Strong Antiviral Activities against Severe Acute Respiratory Syndrome Coronavirus 2. Pathogens 10 (6), 758. doi:10.3390/pathogens10060758
Khattab, M., and Al‐Karmalawy, A. A. (2021). Revisiting Activity of Some Nocodazole Analogues as a Potential Anticancer Drugs Using Molecular Docking and DFT Calculations. Front. Chem. 9, 92. doi:10.3389/fchem.2021.628398
Klopman, G., Stefan, L. R., and Saiakhov, R. D. (2002). ADME Evaluation: 2. A Computer Model for the Prediction of Intestinal Absorption in Humans. Eur. J. Pharm. Sci. 17 (4-5), 253–263. doi:10.1016/s0928-0987(02)00219-1
Liu, D. (2017). Hepatic Cancer. Tumors and Cancers. Boca Raton: CRC Press, 159–164. doi:10.1201/b22275-29
Louisse, J., Bosgra, S., Blaauboer, B. J., Rietjens, I. M. C. M., and Verwei, M. (2015). Prediction of In Vivo Developmental Toxicity of All-Trans-Retinoic Acid Based on In Vitro Toxicity Data and In Silico Physiologically Based Kinetic Modeling. Arch. Toxicol. 89 (7), 1135–1148. doi:10.1007/s00204-014-1289-4
Macfarlane, M., Jones, P., Goebel, C., Dufour, E., Rowland, J., Araki, D., et al. (2009). A Tiered Approach to the Use of Alternatives to Animal Testing for the Safety Assessment of Cosmetics: Skin Irritation. Regul. Toxicol. Pharmacol. 54 (2), 188–196. doi:10.1016/j.yrtph.2009.04.003
Maj, E., Papiernik, D., and Wietrzyk, J. (2016). Antiangiogenic Cancer Treatment: The Great Discovery and Greater Complexity (Review). Int. J. Oncol. 49 (5), 1773–1784. doi:10.3892/ijo.2016.3709
Mannhold, R., Kubinyi, H., and Folkers, G. (2012). Pharmacokinetics and Metabolism in Drug Design, Vol. 51. John Wiley & Sons. doi:10.1002/9783527645763
Marin, J. J. G., Briz, O., Herraez, E., Lozano, E., Asensio, M., Di Giacomo, S., et al. (2018). Molecular Bases of the Poor Response of Liver Cancer to Chemotherapy. Clin. Res. Hepatol. Gastroenterol. 42 (3), 182–192. doi:10.1016/j.clinre.2017.12.006
Marks, P. A., and Xu, W.-S. (2009). Histone Deacetylase Inhibitors: Potential in Cancer Therapy. J. Cel. Biochem. 107 (4), 600–608. doi:10.1002/jcb.22185
McGuire, S. (2016). World Cancer Report 2014. Geneva, Switzerland: World Health Organization, International Agency for Research on Cancer, WHO Press, 2015. Adv. Nutr. 7 (2), 418–419. doi:10.3945/an.116.012211
Metwaly, A. M., Colotta, V., and Catarzi, D. (2019). Discovery and Antiproliferative Evaluation of New Quinoxalines as Potential DNA Intercalators and Topoisomerase. Arch. Pharm. Chem., e1900123.
Montgomery, R. L., Davis, C. A., Potthoff, M. J., Haberland, M., Fielitz, J., Qi, X., et al. (2007). Histone Deacetylases 1 and 2 Redundantly Regulate Cardiac Morphogenesis, Growth, and Contractility. Genes Dev. 21 (14), 1790–1802. doi:10.1101/gad.1563807
Perz, J. F., Armstrong, G. L., Farrington, L. A., Hutin, Y. J. F., and Bell, B. P. (2006). The Contributions of Hepatitis B Virus and Hepatitis C Virus Infections to Cirrhosis and Primary Liver Cancer Worldwide. J. Hepatol. 45 (4), 529–538. doi:10.1016/j.jhep.2006.05.013
Romer, D. R. (2009). Synthesis of 2,3-dichloroquinoxalinesviaVilsmeier Reagent Chlorination. J. Heterocyclic Chem. 46 (2), 317–319. doi:10.1002/jhet.56
Roy, P. P., and Roy, K. (2009). QSAR Studies of CYP2D6 Inhibitor Aryloxypropanolamines Using 2D and 3D Descriptors. Chem. Biol. Drug Des. 73 (4), 442–455. doi:10.1111/j.1747-0285.2009.00791.x
Sabt, A., Abdelhafez, O. M., El-Haggar, R. S., Madkour, H. M. F., Eldehna, W. M., El-Khrisy, E. E.-D. A. M., et al. (2018). Novel Coumarin-6-Sulfonamides as Apoptotic Anti-proliferative Agents: Synthesis, In Vitro Biological Evaluation, and QSAR Studies. J. Enzyme Inhib. Med. Chem. 33 (1), 1095–1107. doi:10.1080/14756366.2018.1477137
Sabt, A., Eldehna, W. M., Al-Warhi, T., Alotaibi, O. J., Elaasser, M. M., Suliman, H., et al. (2020). Discovery of 3,6-disubstituted Pyridazines as a Novel Class of Anticancer Agents Targeting Cyclin-dependent Kinase 2: Synthesis, Biological Evaluation and In Silico Insights. J. Enzyme Inhib. Med. Chem. 35 (1), 1616–1630. doi:10.1080/14756366.2020.1806259
Samra, R. M., Soliman, A. F., Zaki, A. A., Ashour, A., Al-Karmalawy, A. A., Hassan, M. A., et al. (2021). Bioassay-guided Isolation of a New Cytotoxic Ceramide from Cyperus Rotundus L. South Afr. J. Bot. 139, 210–216. doi:10.1016/j.sajb.2021.02.007
Sangwan, R., Rajan, R., and Mandal, P. K. (2018). HDAC as Onco Target: Reviewing the Synthetic Approaches with SAR Study of Their Inhibitors. Eur. J. Med. Chem. 158, 620–706. doi:10.1016/j.ejmech.2018.08.073
Skehan, P., Storeng, R., Scudiero, D., Monks, A., McMahon, J., Vistica, D., et al. (1990). New Colorimetric Cytotoxicity Assay for Anticancer-Drug Screening. JNCI J. Natl. Cancer Inst. 82 (13), 1107–1112. doi:10.1093/jnci/82.13.1107
Soltan, M. A., Elbassiouny, N., Gamal, H., Elkaeed, E. B., Eid, R. A., Eldeen, M. A., et al. (2021). In Silico Prediction of a Multitope Vaccine against Moraxella Catarrhalis: Reverse Vaccinology and Immunoinformatics. Vaccines 9 (6), 669. doi:10.3390/vaccines9060669
Soltane, R., Chrouda, A., Mostafa, A., Al-Karmalawy, A. A., Chouaïb, K., dhahri, A., et al. (2021). Strong Inhibitory Activity and Action Modes of Synthetic Maslinic Acid Derivative on Highly Pathogenic Coronaviruses: COVID-19 Drug Candidate. Pathogens 10 (5), 623. doi:10.3390/pathogens10050623
Tan, I. L., Lycklama à Nijeholt, G. J., Polman, C. H., Adèr, H. J., Barkhof, F., and Group, L. S. (2000). Linomide in the Treatment of Multiple Sclerosis: MRI Results from Prematurely Terminated Phase-III Trials. Mult. Scler. J. 6 (2), 99–104. doi:10.1191/135245800678827626
Van De Waterbeemd, H., and Gifford, E. (2003). ADMET In Silico Modelling: towards Prediction paradise? Nat. Rev. Drug Discov. 2 (3), 192–204. doi:10.1038/nrd1032
Vasudev, N. S., and Reynolds, A. R. (2014). Anti-angiogenic Therapy for Cancer: Current Progress, Unresolved Questions and Future Directions. Angiogenesis 17 (3), 471–494. doi:10.1007/s10456-014-9420-y
Venkatapathy, R., Moudgal, C. J., and Bruce, R. M. (2004). Assessment of the Oral Rat Chronic Lowest Observed Adverse Effect Level Model in TOPKAT, a QSAR Software Package for Toxicity Prediction. J. Chem. Inf. Comput. Sci. 44 (5), 1623–1629. doi:10.1021/ci049903s
Venkatapathy, R., Wang, N. C. Y., Martin, T. M., Harten, P. F., and Young, D. (2009). Structure–Activity Relationships for Carcinogenic Potential. Gen. Appl. Syst. Toxicol. doi:10.1002/9780470744307.gat079
Ververis, K., Hiong, A., Karagiannis, T. C., and Licciardi, P. V. (2013). Histone Deacetylase Inhibitors (HDACIs): Multitargeted Anticancer Agents. Biologics 7, 47–60. doi:10.2147/BTT.S29965
Xia, X., Maliski, E. G., Gallant, P., and Rogers, D. (2004). Classification of Kinase Inhibitors Using a Bayesian Model. J. Med. Chem. 47 (18), 4463–4470. doi:10.1021/jm0303195
Yamashita, T., and Wang, X. W. (2013). Cancer Stem Cells in the Development of Liver Cancer. J. Clin. Invest. 123 (5), 1911–1918. doi:10.1172/jci66024
Zaki, A. A., Al-Karmalawy, A. A., El-Amier, Y. A., and Ashour, A. (2020). Molecular Docking Reveals the Potential of Cleome Amblyocarpa Isolated Compounds to Inhibit COVID-19 Virus Main Protease. New J. Chem. 44 (39), 16752–16758. doi:10.1039/d0nj03611k
Keywords: quinoxaline, anti-proliferative, HDAC, apoptosis, molecular modeling, structure-activity relationship
Citation: Ma C, Taghour MS, Belal A, Mehany ABM, Mostafa N, Nabeeh A, Eissa IH and Al-Karmalawy AA (2021) Design and Synthesis of New Quinoxaline Derivatives as Potential Histone Deacetylase Inhibitors Targeting Hepatocellular Carcinoma: In Silico, In Vitro, and SAR Studies. Front. Chem. 9:725135. doi: 10.3389/fchem.2021.725135
Received: 15 June 2021; Accepted: 29 July 2021;
Published: 22 September 2021.
Edited by:
Simone Brogi, University of Pisa, ItalyReviewed by:
Sobhi M. Gomha, Cairo University, EgyptJose G Trujillo-Ferrara, Instituto Politécnico Nacional (IPN), Mexico
Copyright © 2021 Ma, Taghour, Belal, Mehany, Mostafa, Nabeeh, Eissa and Al-Karmalawy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mohammed S. Taghour, mohammad1533.el@azhar.edu.eg; Ahmed B. M. Mehany, abelal_81@azhar.edu.eg; Ahmed A. Al-Karmalawy, akarmalawy@horus.edu.eg