Design, Synthesis, and Biological Evaluation of a Novel VEGFR-2 Inhibitor Based on a 1,2,5-Oxadiazole-2-Oxide Scaffold with MAPK Signaling Pathway Inhibition

Over the past few decades, the development of broad-spectrum anticancer agents with anti-angiogenic activity has witnessed considerable progress. In this study, a new series of pyrazolo[3,4-d]pyrimidines based on a phenylfuroxan scaffold were designed, synthesized, and evaluated, in terms of their anticancer activities. NCI-60 cell one-dose screening revealed that compounds 12a–c and 14a had the best MGI%, among the tested compounds. The target fluorinated compound 12b, as the most active one, showed better anticancer activity compared to the reference drug sorafenib, with IC50 values of 11.5, 11.6, and 13 µM against the HepG-2, A2780CP, and MDA-MB-231 cell lines, respectively. Furthermore, compound 12b (IC50 = 0.092 µM) had VEGFR-2-inhibitory activity comparable to that of the standard inhibitor sorafenib (IC50 = 0.049 µM). Furthermore, the ability of compound 12b in modulating MAPK signaling pathways was investigated. It was found to decrease the level of total ERK and its phosphorylated form, as well as leading to the down-regulation of metalloproteinase MMP-9 and the over-expression of p21 and p27, thus leading to subG1 cell-cycle arrest and, thus, the induction of apoptosis. Additionally, compound 12b decreased the rate of wound healing in the absence of serum, in comparison to DMSO-treated cells, providing a significant impact on metastasis inhibition. The quantitative RT-PCR results for E-cadherin and N-cadherin showed lower expression of the neuronal N-cadherin and increased expression of epithelial E-cadherin, indicating the ability of 12b to suppress metastasis. Furthermore, 12b-treated HepG2 cells expressed a low level of anti-apoptotic BCL-2 and over-expressed proapoptotic Bax genes, respectively. Using the DAF-FM DA fluorescence probe, compound 12b produced NO intracellularly as efficiently as the reference drug JS-K. In silico molecular docking studies showed a structural similarity through an overlay of 12b with sorafenib. Interestingly, the drug-likeness properties of compound 12b met the expectations of Pfizer’s rule for the design of new drug candidates. Therefore, this study presents a novel anticancer lead compound that is worthy of further investigation and activity improvement.


Introduction
Cancer is a large group of diseases, and is the result of the uncontrollable growth of abnormal cells [1]. In 2018, according to the WHO, the second leading cause of death worldwide was cancer, and it was estimated that 9.6 million people worldwide died from various cancers [1]. About 300,000 new cases are diagnosed every year among children  Furthermore, the pyrimidine-based core was selected as the ATP-binding domain isostere. It has recently been reported that numerous pyrazolo [3,4-d]pyrimidine derivatives exhibit potent anti-tumor activities [29][30][31][32]. In 2018, Wand et al. reported a potent VEGFR-2 inhibitor having a pyrazolo [3,4-d]pyrimidine scaffold [20]. Therefore, the privileged pyrazolo [3,4-d]pyrimidine nucleus was deemed to be a promising scaffold for developing a new VEGFR-2 inhibitor class. Furthermore, the pyrimidine-based core was selected as the ATP-binding domain isostere. It has recently been reported that numerous pyrazolo [3,4-d]pyrimidine derivatives exhibit potent anti-tumor activities [29][30][31][32]. In 2018, Wand et al. reported a potent VEGFR-2 inhibitor having a pyrazolo [3,4-d]pyrimidine scaffold [20]. Therefore, the privileged pyrazolo [3,4-d]pyrimidine nucleus was deemed to be a promising scaffold for developing a new VEGFR-2 inhibitor class. On the other hand, the nitric oxide (NO)-releasing furoxan moiety was shown to have potent anti-tumor activity in earlier reports [21,22]. It is worth mentioning that NO plays a pivotal role as a signaling free radical molecule, functioning as a vascular smooth muscle relaxant, platelet aggregation inhibitor, neurotransmitter, and immune regulator [33]. The anti-tumor activity of NO has been observed at high levels (more than 300 nM), through the initiation of DNA damage and nitrosative stress [34,35], the phosphorylation of P53, and the increasing of MKP-1 expression, leading to cellular respiration inhibition. High levels of NO also play a crucial role in the reversal of drug resistance through MRP efflux pump inhibition [36,37], hypoxia-induced drug-resistance attenuation [36], hypoxic radio-sensitization [38], transcription factor hypoxia-inducible factor (HIF-1) activation blocking [36,38], inhibition of epithelial-mesenchymal transition (EMT) [39], and metastatic inhibition [39]. Therefore, increasing the concentration of intracellular NO may have a significant impact on tumor growth and metastasis [35,36].
A convenient way to achieve high levels of intracellular NO is through the use of NO donors [40]. Among them, the most stable is the furoxan moiety [41], which releases NO in the presence of plasma, reduced glutathione (GSH), or albumin (i.e., through thiol mechanisms) and, therefore, may exert better biological activity and rarely leads to tolerance, compared to other nitric oxide-releasing motifs [40,41]. In summary, the synergistic effects discussed above have raised hopes for its therapeutic potential and, as such, endowing the pyrazolo [3,4-d]pyrimidine scaffold with a furoxan moiety may have a significant impact on tumor and metastasis suppression [20][21][22]41].
Enlightened by these findings, we were inspired to design and synthesize a new series of pyrazolo [3,4-d]pyrimidines tethered to a NO-releasing furoxan moiety through a central aryl bridge. In this context, as shown in Figure 3, our program began with the replacement of the pyridine ring with an isosteric pyrazolo [3,4-d]pyrimidine base-scaffold (in order to fit the ATP-binding domain), the replacement of the oxygen atom with the bioisosteric NH as a hydrogen bond donor (HBD), the keeping of the central aryl moiety without modification as a linker, the replacement of the hydrogen-bonding moiety urea with furoxan methoxy/amino to bind the DFG domain and, finally, the incorporating of the phenyl group as a hydrophobic tail to occupy the allosteric binding region. Furthermore, C3 and C6 functionalizations were studied. On the other hand, the nitric oxide (NO)-releasing furoxan moiety was shown to have potent anti-tumor activity in earlier reports [21,22]. It is worth mentioning that NO plays a pivotal role as a signaling free radical molecule, functioning as a vascular smooth muscle relaxant, platelet aggregation inhibitor, neurotransmitter, and immune regulator [33]. The anti-tumor activity of NO has been observed at high levels (more than 300 nM), through the initiation of DNA damage and nitrosative stress [34,35], the phosphorylation of P53, and the increasing of MKP-1 expression, leading to cellular respiration inhibition. High levels of NO also play a crucial role in the reversal of drug resistance through MRP efflux pump inhibition [36,37], hypoxia-induced drug-resistance attenuation [36], hypoxic radiosensitization [38], transcription factor hypoxia-inducible factor (HIF-1) activation blocking [36,38], inhibition of epithelial-mesenchymal transition (EMT) [39], and metastatic inhibition [39]. Therefore, increasing the concentration of intracellular NO may have a significant impact on tumor growth and metastasis [35,36].
A convenient way to achieve high levels of intracellular NO is through the use of NO donors [40]. Among them, the most stable is the furoxan moiety [41], which releases NO in the presence of plasma, reduced glutathione (GSH), or albumin (i.e., through thiol mechanisms) and, therefore, may exert better biological activity and rarely leads to tolerance, compared to other nitric oxide-releasing motifs [40,41]. In summary, the synergistic effects discussed above have raised hopes for its therapeutic potential and, as such, endowing the pyrazolo [3,4-d]pyrimidine scaffold with a furoxan moiety may have a significant impact on tumor and metastasis suppression [20][21][22]41].
Enlightened by these findings, we were inspired to design and synthesize a new series of pyrazolo [3,4-d]pyrimidines tethered to a NO-releasing furoxan moiety through a central aryl bridge. In this context, as shown in Figure 3, our program began with the replacement of the pyridine ring with an isosteric pyrazolo [3,4-d]pyrimidine base-scaffold (in order to fit the ATP-binding domain), the replacement of the oxygen atom with the bioisosteric NH as a hydrogen bond donor (HBD), the keeping of the central aryl moiety without modification as a linker, the replacement of the hydrogen-bonding moiety urea with furoxan methoxy/amino to bind the DFG domain and, finally, the incorporating of the phenyl group as a hydrophobic tail to occupy the allosteric binding region. Furthermore, C3 and C6 functionalizations were studied.

Chemistry
The synthetic process for the target compounds was initiated with a trial to prepare the furoxan moiety from allyl ether precursors 3a-d. As depicted in Scheme 1, the phenolic derivatives 2a-h were prepared from alkylidyne malononitrile, through a multi-step reaction, with reasonable yields [42]. Compounds 2a-d were alkylated, using allyl bromide, to obtain the allyloxy derivatives 3a-d [42][43][44][45]. Upon treating compounds 3a-d with sodium nitrite, an electrophilic aromatic substitution reaction took place at the ortho position to an allyloxy group, wherein the carbons were more shielded, thus providing the nitrated derivatives 5a-d, rather than the desired compound 4 [46]. The unaffected chemical shift of allyl protons through the reaction and disappearance of the doublet at 7.00 ppm corresponding to the protons ortho to allyloxy group evidenced the proposed chemical structures. On the other hand, the furoxan moiety was successfully installed through the alkylation of phenolic OH precursors 2a-h using phenyl furoxan methyl mesylate to afford the target compounds 6a-h [21].
Moreover, based on the above findings, compound 12b was identified as the most active analog among the tested compounds, displaying a better anticancer profile. As shown in Tables S1 (from the Supplementary Material File), compound 12b exhibited growth-inhibitory % activity against the tested panel of up to 77% in leukemia, 64% in non-small cell lung cancer, 62% in colon cancer, 45% in CNS cancer, 72% in melanoma, 40% in ovarian cancer, 69% in renal cancer, 52% in prostate cancer, and 61% in breast cancer cells. Hence, the anticancer profile of compound 12b against the NCI-59 cancer cell lines encouraged us to further explore its mode of action.

Preliminary SARs Study
The methyl substitution (EDG) at C3 and C6 (i.e., compounds 6a-h) had no significant effect on the anticancer activity. Notably, the replacement of EDG at C6 in the tested series (i.e., compounds 6a-h) with the electron-withdrawing group CF3, and the replacement of the oxygen of the central aryl ring with NH, caused compounds 12a-c to exhibit significant inhibitory effects in their anticancer activity, indicating that the electron-withdrawing substituents at C6 very likely led to enhanced potency. On the other hand, moving the substitution from C4 (compounds 12a-c) to C6, as in compounds 14a and 14b, reduced the growth-inhibitory activity while compound 14c enhanced the growth, indicating that appending the furoxan moiety at C6 and keeping both the pyrazole C3 and 1phenyl ring unsubstituted may have a good impact on the activity. Furthermore, to test the hybridization hypothesis, compounds 11 and 13a-c were also tested, and the results displayed a significant improvement in anticancer activity, especially for compound 14a, As shown in Figure 4, considering the mean growth inhibition % (MGI%) at 10 µM, compounds 6a and f demonstrated no mean growth inhibition percentage, while compounds 6b, d, and e enhanced the growth percentage by 1-5%. Compounds 6c, g, and h demonstrated low anticancer activity, with GI% = 1.58-4.57%. Additionally, compounds 12a-c exhibited significant inhibitory effects in their anticancer activity, with GI% in the range of 23.2-33.17%.
On the other hand, moving the substitution from C4 (compounds 12a-c) to C6, as in compounds 14a and 14b, reduced the growth-inhibitory activity (exhibiting GI% of 27.93 and 8.22, respectively), while compound 14c enhanced the growth by 3.35%.
Moreover, based on the above findings, compound 12b was identified as the most active analog among the tested compounds, displaying a better anticancer profile. As shown in Table S1 (from the Supplementary Material File), compound 12b exhibited growth-inhibitory % activity against the tested panel of up to 77% in leukemia, 64% in non-small cell lung cancer, 62% in colon cancer, 45% in CNS cancer, 72% in melanoma, 40% in ovarian cancer, 69% in renal cancer, 52% in prostate cancer, and 61% in breast cancer cells. Hence, the anticancer profile of compound 12b against the NCI-59 cancer cell lines encouraged us to further explore its mode of action.

Preliminary SARs Study
The methyl substitution (EDG) at C3 and C6 (i.e., compounds 6a-h) had no significant effect on the anticancer activity. Notably, the replacement of EDG at C6 in the tested series (i.e., compounds 6a-h) with the electron-withdrawing group CF 3 , and the replacement of the oxygen of the central aryl ring with NH, caused compounds 12a-c to exhibit significant inhibitory effects in their anticancer activity, indicating that the electron-withdrawing substituents at C6 very likely led to enhanced potency. On the other hand, moving the substitution from C4 (compounds 12a-c) to C6, as in compounds 14a and 14b, reduced the growth-inhibitory activity while compound 14c enhanced the growth, indicating that appending the furoxan moiety at C6 and keeping both the pyrazole C3 and 1-phenyl ring unsubstituted may have a good impact on the activity. Furthermore, to test the hybridization hypothesis, compounds 11 and 13a-c were also tested, and the results displayed a significant improvement in anticancer activity, especially for compound 14a, indicating synergistic efficiency. Consequently, appending the furoxan moiety at C4, utilizing an NH in the central aryl ring at C4, and substituting the C6 with an electron-withdrawing CF 3 group can be suggested as vital elements for the growth-inhibitory activity of the furoxan-based pyrazolo [3,4-d]pyrimidine hybrids reported herein. A summary, based on the above discussion, of the preliminary structure-activity relationships (SARs) observed in this study is provided in Figure 5. based on the above discussion, of the preliminary structure-activity relationships (SARs) observed in this study is provided in Figure 5. In Vitro Anti-Proliferative Activities The cell growth-inhibitory activities of the most active compound 12b were evaluated against human cancer HepG-2 [52][53][54], ovarian [55,56], breast [54], and colon [54] cell lines highly expressing the VEGFR-2 protein. The cells were treated with serial dilutions for 48 h, after which the cell viability was determined by standard MTT assay, using sorafenib as a reference drug. The initial screening of the target compounds showed that 12b was the most active compound, at a single 50 µM dose, against liver (HepG2), breast (MDA-MB-231 and MCF-7), Colon (HT-29), and ovarian (SKOV-3, A2780, and A2780CP) cancer cell lines. In order to determine the concentration that killed 50% of cells, serial dilutions were used in the treatment, and the determination of percentage of viability identified that 12b was the most active compound. The IC50 values are shown in Table 1 and Figure 6. Compared to sorafenib, the target compound 12b showed better anticancer activity, with IC50 values of 11.5, 11.6, and 13 µM against HepG-2, A2780CP, and MDA-MB-231 cell lines, respectively, while the IC50 for SKOV-3 was comparable to that of the reference drug sorafenib. In addition, the morphological changes in 12b-treated cancer cell lines at concentrations of 6.25 and 12.50 µM are shown in Figure 6. In contrast, the des-NO-releasing compounds 15a and b were tested, but did not inhibit the proliferation of cancer cells at concentrations as high as 50 µM, indicating the significance of the phenylfuroxan moiety in the anticancer activity of 12b.  In Vitro Anti-Proliferative Activities The cell growth-inhibitory activities of the most active compound 12b were evaluated against human cancer HepG-2 [52][53][54], ovarian [55,56], breast [54], and colon [54] cell lines highly expressing the VEGFR-2 protein. The cells were treated with serial dilutions for 48 h, after which the cell viability was determined by standard MTT assay, using sorafenib as a reference drug. The initial screening of the target compounds showed that 12b was the most active compound, at a single 50 µM dose, against liver (HepG2), breast (MDA-MB-231 and MCF-7), Colon (HT-29), and ovarian (SKOV-3, A2780, and A2780CP) cancer cell lines. In order to determine the concentration that killed 50% of cells, serial dilutions were used in the treatment, and the determination of percentage of viability identified that 12b was the most active compound. The IC 50 values are shown in Table 1 and Figure 6. Compared to sorafenib, the target compound 12b showed better anticancer activity, with IC 50 values of 11.5, 11.6, and 13 µM against HepG-2, A2780CP, and MDA-MB-231 cell lines, respectively, while the IC 50 for SKOV-3 was comparable to that of the reference drug sorafenib. In addition, the morphological changes in 12b-treated cancer cell lines at concentrations of 6.25 and 12.50 µM are shown in Figure 6. In contrast, the des-NO-releasing compounds 15a and b were tested, but did not inhibit the proliferation of cancer cells at concentrations as high as 50 µM, indicating the significance of the phenylfuroxan moiety in the anticancer activity of 12b.  Furthermore, the selectivity index (SI) of compound 12b towards the various cancer cell lines was evaluated by detecting the required concentration to kill 50% of normal  Furthermore, the selectivity index (SI) of compound 12b towards the various cancer cell lines was evaluated by detecting the required concentration to kill 50% of normal human skin fibroblast cells (HSF). The results revealed that compound 12b did not cause any morphological changes in HSF cells at a concentration as high as 50 µM, indicating its safety and selectivity towards cancer cells, rather than normal cells.

VEGFR Kinase Inhibitory Profile of 12b
The fluorinated pyrazolo [3,4-d]pyrimidine derivative based on the 1,2,5-oxadiazole-2-oxide scaffold (12b) was evaluated for its VEGFRx (VEGFR-1, VEGFR-2, and VEGFR-3) kinase inhibitory activities, using a VEGFR inhibitory assay with sorafenib as a reference drug. As shown in Figure 7, the results indicated that compound 12b possesses potency that is comparable to that of sorafenib, displaying noticeable selectivity against VEGFR-2 with an IC 50 of 0.09 µM. human skin fibroblast cells (HSF). The results revealed that compound 12b did not cause any morphological changes in HSF cells at a concentration as high as 50 µM, indicating its safety and selectivity towards cancer cells, rather than normal cells.

Flow Cytometric Studies
The mechanism of action of 12b was investigated by evaluating its effectiveness in modulating the cell cycle and induction of apoptosis. As seen in Figure 8A

Flow Cytometric Studies
The mechanism of action of 12b was investigated by evaluating its effectiveness in modulating the cell cycle and induction of apoptosis. As seen in Figure 8A,B, treating HepG2 cells with compound 12b for 48 h caused cell-cycle arrest at subG1 (24 vs. 5.6 in DMSO-treated cells), leading to the induction of apoptotic pathways ( Figure 8C,D).

Intracellular Measurement of NO
The NO-releasing properties of 12b were evaluated by measuring the level of intracellular NO by flow cytometry, using a cell-permeable and photo-stable NO fluorescent indicator, 3-Amino-4-aminomethyl-2',7'-difluorofluorescein diacetate (DAF-FM DA), with a detection limit of ~3 nM. DAF-FM diacetate penetrates the cells and is deacetylated by esterase enzyme to yield a non-fluorescent dye, DAF-FM, which, upon reaction with NO, produces a fluorescent benzotriazole product [57]. Compound 12b was efficiently able to increase the production of NO in the treated cells (9.9% vs. 0.1% in DMSO-treated cells; Figure 9). Furthermore, compound 12b increased NO production as efficiently as the positive control reference drug JS-k.

Intracellular Measurement of NO
The NO-releasing properties of 12b were evaluated by measuring the level of intracellular NO by flow cytometry, using a cell-permeable and photo-stable NO fluorescent indicator, 3-Amino-4-aminomethyl-2',7'-difluorofluorescein diacetate (DAF-FM DA), with a detection limit of~3 nM. DAF-FM diacetate penetrates the cells and is deacetylated by esterase enzyme to yield a non-fluorescent dye, DAF-FM, which, upon reaction with NO, produces a fluorescent benzotriazole product [57]. Compound 12b was efficiently able to increase the production of NO in the treated cells (9.9% vs. 0.1% in DMSO-treated cells; Figure 9). Furthermore, compound 12b increased NO production as efficiently as the positive control reference drug JS-k.

Wound-Healing Assay
A wound-healing assay is a simple and economical assay used to investigate the antimetastatic activity of compounds. Cells were treated with compound 12b in serum-free media for 24 h. After that, the cells were fixed with ice-cold methanol and air-dried, then stained with 0.5% crystal violet. As shown in Figure 10, compound 12b decreased the rate of wound healing in the absence of serum, in comparison to the DMSO-treated cells. These data revealed that 12b is an anticancer compound that possesses anti-metastatic activity.

Wound-Healing Assay
A wound-healing assay is a simple and economical assay used to investigate the antimetastatic activity of compounds. Cells were treated with compound 12b in serum-free media for 24 h. After that, the cells were fixed with ice-cold methanol and air-dried, then stained with 0.5% crystal violet. As shown in Figure 10, compound 12b decreased the rate of wound healing in the absence of serum, in comparison to the DMSO-treated cells. These data revealed that 12b is an anticancer compound that possesses anti-metastatic activity.

Apoptosis and Metastatic Proteins
Epithelial cells express specific genes involved in maintaining apical-basal polarity and adherence, such as epithelial E-cadherin. The epithelial cells migrate to distant locations by over-expressing mesenchymal genes, including neuronal N-cadherin and vimentin [58]. Cancer aggressiveness has been positively correlated with the over-expression of N-cadherin. Therefore, in this study, the partial mechanism of activity of 12b was investigated, by quantitative RT-PCR, with respect to E-cadherin and N-cadherin. The 12b-treated cells showed lower expression of neuronal N-cadherin and increased expres-

Apoptosis and Metastatic Proteins
Epithelial cells express specific genes involved in maintaining apical-basal polarity and adherence, such as epithelial E-cadherin. The epithelial cells migrate to distant locations by over-expressing mesenchymal genes, including neuronal N-cadherin and vimentin [58]. Cancer aggressiveness has been positively correlated with the over-expression of N-cadherin.
Therefore, in this study, the partial mechanism of activity of 12b was investigated, by quantitative RT-PCR, with respect to E-cadherin and N-cadherin. The 12b-treated cells showed lower expression of neuronal N-cadherin and increased expression of epithelial Ecadherin, indicating the ability of 12b to suppress metastasis ( Figure 11A). Additionally, the ability of 12b was confirmed by detecting the change in expression of both pro-apoptotic (Bax) and anti-apoptotic (BCL-2) genes. Compound 12b-treated HepG2 cells (for 12 h) expressed a low level of BCL-2 and over-expressed Bax ( Figure 11A).

Apoptosis and Metastatic Proteins
Epithelial cells express specific genes involved in maintaining apical-basal polarity and adherence, such as epithelial E-cadherin. The epithelial cells migrate to distant locations by over-expressing mesenchymal genes, including neuronal N-cadherin and vimentin [58]. Cancer aggressiveness has been positively correlated with the over-expression of N-cadherin. Therefore, in this study, the partial mechanism of activity of 12b was investigated, by quantitative RT-PCR, with respect to E-cadherin and N-cadherin. The 12b-treated cells showed lower expression of neuronal N-cadherin and increased expression of epithelial E-cadherin, indicating the ability of 12b to suppress metastasis ( Figure  11A). Additionally, the ability of 12b was confirmed by detecting the change in expression of both pro-apoptotic (Bax) and anti-apoptotic (BCL-2) genes. Compound 12b-treated HepG2 cells (for 12 h) expressed a low level of Bcl-2 and over-expressed Bax ( Figure 11A).  The mitogen-activated protein kinase (MAPK) pathway is considered to be a key player in the gene expression of survival, proliferative, and anti-apoptotic proteins. Aberrant dysregulation in the MAPK pathway activity has been linked to the overall poor prognosis in various cancers [59]. Compound 12b decreased the level of total ERK and its phosphorylated form ( Figure 11B), indicating its ability to modulate signaling pathways. Inhibition of the MAPK pathway led to the downregulation of metalloproteinase MMP-9 (metastatic protein), the over-expression of cell-cycle inhibitors (p21 and p27), and a slight increase in the expression of p53; see Figure 11B. Additionally, compound 12b was further tested against p38, MKK3, and JNK using ELISA assay protocols [60][61][62] at a concentration of 11 µM. The results indicated that 12b caused low expression in p38 and MKK3 while JNK was overexpressed; see Figure 12. The overexpression of JNK might have contributed to the NO released from 12b and reactive oxygen species (ROS). In summary, the novel anticancer lead compound 12b is worthy of further investigation and activity improvement. of 11 µM. The results indicated that 12b caused low expression in p38 and MKK3 wh JNK was overexpressed; see Figure 12. The overexpression of JNK might have contribut to the NO released from 12b and reactive oxygen species (ROS). In summary, the nov anticancer lead compound 12b is worthy of further investigation and activity improv ment.

EON Scaffold Hopping
The discovery of new molecular targets can be achieved by whole-molecule replac ments, through the use of physically realistic shapes and electrostatic potential similariti that permit the identification of biologically active molecules with significantly differe structures, with respect to existing active agents [63]. OpenEye's EON application [64] c contribute significantly to lead generation and library design, through the creation of ele trostatic Tanimoto (ET) grids using full Poisson-Boltzmann (PB) electrostatics and attac ing them to each output molecule [65]. EON is also influenced by the pKa state and form charges, which have major impacts on the electrostatics, as the query and database mol cule are adjusted to a neutral pH model [65]. The data manager in the VIDA module c easily read the EON report file [66]. Sorafenib was included in the library in order to va idate the EON results. The designed library was tested against sorafenib as a query refe ence, and the compounds were ranked according to their shape and electrostatic similari to sorafenib. As shown in Table 2, the top-scoring compounds were those of the C6-a pending series 14a-c, which showed the best alignment, followed by the fluorinated d rivatives 12a-c and, then, the rest of the designed compounds.

In Silico Molecular Modeling 2.3.1. EON Scaffold Hopping
The discovery of new molecular targets can be achieved by whole-molecule replacements, through the use of physically realistic shapes and electrostatic potential similarities that permit the identification of biologically active molecules with significantly different structures, with respect to existing active agents [63]. OpenEye's EON application [64] can contribute significantly to lead generation and library design, through the creation of electrostatic Tanimoto (ET) grids using full Poisson-Boltzmann (PB) electrostatics and attaching them to each output molecule [65]. EON is also influenced by the pKa state and formal charges, which have major impacts on the electrostatics, as the query and database molecule are adjusted to a neutral pH model [65]. The data manager in the VIDA module can easily read the EON report file [66]. Sorafenib was included in the library in order to validate the EON results. The designed library was tested against sorafenib as a query reference, and the compounds were ranked according to their shape and electrostatic similarity to sorafenib. As shown in Table 2, the top-scoring compounds were those of the C6-appending series 14a-c, which showed the best alignment, followed by the fluorinated derivatives 12a-c and, then, the rest of the designed compounds.

Docking Studies
In silico molecular docking studies were performed using the OpenEye Scientific Software, version 2021.spr, in order to explore the binding modes of the target compounds with the VEGFR-2 active site, adopting the crystal structure of VEGFR-2 in complex with sorafenib (PDB ID: 3WZE). The tested compounds were ranked according to their FRED Chemgauss4 scores-the sum of shape, HB, and desolvation energies-with the lowest score indicating the best binding interactions (see Table 3). To validate the docking results and rationalize the predictive protein-ligand interactions of the target compounds within the active site, the positive drug sorafenib was re-docked using the same procedure as the target compounds. The results demonstrated that the docking pose of the co-crystallized ligand (see Figure 13A), showing interactions similar to those existing in the crystal structure, and the hydrophobic interactions were roughly the same. Table 3. FRED Chemgauss4 scores of target compounds 6a-h, 12a-c, and 14a-   The most active compound against NCI-60 cell lines, 12b, was selected as the most plausible binding conformation and, thus, we explored its binding mode within the binding pocket. Compound 12b displayed a unique binding mode that could underlie its broadspectrum anticancer activities. As shown in Figure 13B, in spite of not displaying a strong H-bond, 12b stacked well with sorafenib and tightly fit into the protein pocket without any projection to the protein cavity. The depicted 2D interaction of the protein with 12b (see Figure 13C,D), generated by Discovery Studio Visualizer [67,68], showed two weak hydrogen bonds, coming from the NH group (HBD) on the furoxan methylamino linker and Glu:885:A (HBA, bond length 4.79Å), as confirmed by the 2D depiction of the protein-ligand interaction using the OpenEye docking report ( Figure S75; from the Supplementary Material), as well as between fluorine and Cys:919:A (HBD). Furthermore, there were hydrophobic interactions with the contacted residues. The central aryl ring, pyrazole, and 1-phenyl formed pi-cation electrostatic, pi-pi T-shaped, and pi-pi stacking hydrophobic interactions, respectively. Additionally, 12b formed Van der Waals and other mixed hydrophobic interactions. The docking report of OpenEye revealed that the furoxan-N-oxide formed a hydrogen bond (3%) with Lys:868:A (HBD), as shown in Figure S75 (from the Supplementary Material File). In addition, the FRED docking report showed an additional property relating to the guidance of the optimization of compound 12b, displaying shape, hydrogen bonding interaction, and protein and ligand desolvation energies, identified through the FRED Chemgauss4 score. Compound 12b had a FRED chemgauss4 score of −11.37 (where that of sorafenib was −18.58). Thus, either the hydrogen bonds or hydrophobic interactions of 12b were weaker than those of sorafenib. Simultaneously, as a whole molecule, 12b spatially stacked well with the co-crystallized ligand and tightly fit into the pocket of VEGFR-2; this provided a reasonable explanation as to why compound 12b had better anticancer activity among the tested analogs, as confirmed by the NCI-60 cell line screening and its molecular mechanisms.

ADME/Toxicity Analysis
Prediction of the physicochemical, pharmacokinetic, and drug-likeness properties of compound 12b was achieved using the web-based tool ADMETlab 2.0 [69] which offers a straightforward approach to the comprehensive, accurate, and efficient prediction of ADMET properties using a high-quality database of 250,000 entries covering 53 endpoints and a multi-task graph attention framework [70]. As shown in the radar graph in Figure S76 (from the Supplementary Material File), compound 12b displayed physico-chemical properties comparable to those of sorafenib [69]. In addition, as shown in Table  S3 (from the Supplementary Material File), the ADMET analysis revealed that the druglikeness properties of compound 12b met the expectations of Pfizer's rule for designing new drug candidates. The synthetic accessibility score (SAscore), which was designed to estimate the ease-of-synthesis of drug-like molecules, showed a value of 3.406, thus indicating the ease of 12b synthesis and its good economic cost. Furthermore, it displayed an excellent absorption profile. Notably, its distribution, metabolism, excretion, and toxicity profiles were comparable to those of sorafenib. Furthermore, according to the toxicophoric rules, the acute toxicity during oral administration for 12b did not show any alerts, in contrast with the reference drug sorafenib. On the other hand, the environmental toxicity rules uncovered that the bioconcentration factors, which are used for considering secondary poisoning potential and assessing risks to human health through the food chain, showed a value of 1. 38

Chemistry
1 H and 13 C NMR spectra were obtained using a Bruker AVANCE and Varian Unity INOVA-400 MHz NMR spectrometer in DMSO-d 6 using solvent residual peaks as an internal standard; CDCl 3 (δ = 7.27 ppm) and DMSO-d 6 (δ = 2.50 ppm) for 1 H NMR and CDCl 3 (δ = 77.23 ppm) and DMSO-d 6 (δ = 39.51 ppm) for 13 C NMR. The signal multiplicities were reported as s (singlet), d (doublet), t (triplet), q (quartet), dd (doublet of doublets), m (multiplet), or br (broad) and the J-coupling constants were in Hz. High-resolution mass spectra (HRESI-MS) were acquired on a Thermo Instruments MS system (LTQ XL/LTQ Orbitrap Discovery) linked to a Thermo Instruments HPLC system (Accela PDA detector, Accela PDA autosampler and pump). Thin-layer chromatography (TLC) was performed for reaction monitoring on pre-coated silica gel 60 F 254 (Merck) sheets and visualized by ultraviolet light (254 nm). Purification of synthesized compounds was achieved via either column chromatography using 150-250 µm silica gel or recrystallization from an appropriate solvent. Solvents and reagents were acquired from commercially available suppliers.

General Procedure for Preparation of Allyloxy Derivatives (3a-d)
To a stirred solution of phenolic derivatives 2a-d (0.002 mol) in absolute ethanol (15 mL), allyl bromide (1.209 g, 0.010 mol) and sodium hydroxide (0.080 g, 0.002 mol) were added. After 15 min of reflux, another equivalent of sodium hydroxide was added, and the reflux was continued for another 30 min (TLC monitoring). After cooling, the precipitated solid was filtered, washed with ethanol, dried, and recrystallized from acetonitrile to yield the allyloxy derivatives 3a-d [43][44][45].

General Procedure for Preparation of Nitrated Allyloxy Derivatives (5a-d)
To a flask charged with allyloxy derivatives 3a-d (0.015 mol) in glacial acetic acid (3 mL), saturated aqueous NaNO 2 (0.045 mol) solution was added dropwise so that the temperature did not exceed 70 • C [46]. After 1 h, the mixture was diluted with water and extracted with EtOAc. The organic layer was washed with water and brine, dried over anhydrous sodium sulfate, and concentrated in vacuo. The residue was then purified by column chromatography using hexane/ethyl acetate (8:2) as eluent to yield the nitrated derivatives 5a-d.

General Procedure for Preparation of Target Compounds 6a-h
To a flask charged with phenolic derivatives 2a-h (0.002 mol) in THF (15 mL), furoxan mesylate 10 (0.540 g, 0.002 mol) and NaOH (0.080 g, 0.002 mol) were added and the mixture was heated at reflux for 2 h (TLC monitoring) [21]. Then, the mixture was cooled to room temperature, diluted with water, and extracted with EtOAc. The organic layer was washed with brine, dried over anhydrous sodium sulfate, and concentrated in vacuo. The residue was then purified by column chromatography using EtOAc:Hex (1:4) as eluent to provide the target compounds 6a-h.

General Procedure for Preparation of Pyrazolopyrimidinone Derivatives 8b,c
To a suspension of 5-amino-4-cyanopyrazole intermediates 7b,c (0.002 mol) in TFA (6 mL), phosphorus oxychloride (0.4 mL) was added instantly and the mixture was heated at reflux for 2 h [47]. Then, the reaction mixture was cooled to room temperature and diluted with ice-water. The filtered precipitate was recrystallized from formic acid to afford compounds 8b,c in good yields.

General Procedure for Preparation of 4-Chloro Derivatives 9a-c
A flask charged with suspension of pyrazolopyrimidinone derivatives 8a-c (0.002 mol) in POCl 3 (16 mL) was heated at reflux for 10 h [48]. Then, the mixture was cooled to room temperature and added dropwise to ice-cooled water. The solid formed was filtered, washed with water, dried, and recrystallized from n-hexane to yield the 4-chlorinated derivatives 9a-c. To a flask charged with a solution of furoxan mesylate 10 (0.270 g, 0.001 mol) in THF (10 mL), p-phenylendiamine (0.642 g, 0.006 mol) and TEA (0.001 mol) were added, and the reaction mixture was heated at reflux for 2 h (TLC monitoring) [21]. Then, the mixture was cooled to room temperature and extracted with EtOAc, then washed with hot water 3 times or until the disappearance of unreacted p-phenylenediamine. The organic layer was dried over anhydrous Na 2 SO 4 and concentrated in vacuo. The residue was then purified by column chromatography using EtOAc:Hex (1:1) to give the phenylfuroxan derivative 11 in good yield.

General Procedure for Preparation of C6 Appending Target Compounds 14a-c
To a stirred solution of 6-chloromethyl derivatives 13a-c (0.001 mol) in THF (10 mL), phenylfuroxan derivative 11 (0.282 g, 0.001 mol), TEA (0.001 mol), and NaI (0.0005 mol) were added, and the reaction mixture was heated under reflux for 2 h (TLC monitoring) [50]. Then, the mixture was cooled to room temperature, diluted with water, and extracted with EtOAc. The organic layer was dried over anhydrous Na 2 SO 4 then concentrated in vacuo to afford the crude product, which was purified by column chromatography using EtOAc:Hex (3:7) as eluent to yield the C6 appended target compounds 14a-c. To a stirred solution of 6-trifluoromethyl-4-chloro derivatives 9b,c (0.001 mol) in isopropanol (10 mL), p-phenylenediamine (0.108 g, 0.001 mol), TEA (0.001 mol), and NaI (0.0005 mol) were added and the reaction mixture was heated under reflux for 2 h (TLC monitoring) [42]. Then, the mixture was cooled to room temperature and diluted with water. The formed precipitate was filtered and recrystallized form absolute ethanol to give the des-furoxan derivatives 15a,b in good yields.

NCI-60 Cell Line Screening
The in vitro NCI-60 cell line screening assay was conducted at the National Cancer Institute (NCI), Bethesda, USA against 59 cancer cell lines representing 9 types of cancers, as reported earlier [51,73]. Briefly, the human cancer cell lines were grown in RPMI 1640 medium containing 5% fetal bovine serum and 2 mM L-glutamine and inoculated into 96-well microtiter plates in 100 µL (5000 to 40,000 cells/well). Then, the plates were incubated at 37 • C, 5% CO 2 , 95% air, and 100% relative humidity for 24 h. After that, two plates of each cell line were fixed in situ with TCA, to represent cell population at time zero (Tz). The tested compounds were dissolved in DMSO at 400-fold and added so that the final concentration became 10 µM. The plates were then incubated for an additional 48 h followed by the addition of cold TCA to terminate the assay in case of adherent cells. Then, the cells were fixed in situ by cold 50% (w/v) TCA (50 µL; final concentration, 10% TCA) and incubated for 1 h at 4 • C. After discarding the supernatant, the plates were washed 5 times with tap water and air dried. Sulforhodamine B (SRB) solution (100 µL) at 0.4% (w/v) in 1% HOAc was added to each well, and plates were incubated for 10 min. at rt. Then, washing with 1% acetic acid was conducted 5 times to remove the unbound dye followed by air drying. After that, 10 mM trizma base was added to solubilize the bound stain, and the absorbance was read on an automated plate reader at a wavelength of 515 nm. On the other hand, the methodology for suspension cells was the same except that 50 µL of 80% TCA (final concentration, 16% TCA) was gently added to terminate the assay by fixing settled cells at the bottom of the wells. Percentage growth inhibition (GI%) was calculated as follows: Growth × 100 for Ti < Tz.
Data are expressed as a one-dose mean graph of growth percentages including the mean, delta, and range.

VEGFRx Kinase Assay
The kinase assay was performed according to BPS Bioscience ® kinase assay kit protocols [17,[74][75][76][77]. Briefly, 5x Kinase Buffer 1, ATP and Poly-(Glu,Tyr 4:1) (10 mg/mL) in case of FLT1 (Cat. No. 78019)/50x PTK substrate in case of KDR (Cat. No. 40325)/MBP (5 mg/mL) in case of FLT3 (Cat. No. 79797) were thawed and the master mixture (25 µL per well: N wells x (5x Kinase Buffer 1 + ATP (500 µM) + substrate + water)) was prepared. The master mixture (25 µL) was added to every well then inhibitor solution of each well, labeled as "Test Inhibitor", was added. For the "Positive Control" and "Blank", the same solution without inhibitor (Inhibitor buffer) was added. To the wells designated as "Blank", 1x Kinase Buffer 1 was added. FLT1/VEGFR2/FLT3 enzyme was thawed on ice and diluted with 1x Kinase Buffer 1 then added to the wells designated "Positive Control" and "Test Inhibitor Control". The plate was incubated at 30 • C for 45 min then ADP-Glo reagent/Kinase-Glo Max reagent was added to each well. The plate was covered with aluminum foil and incubated at room temperature for 15 min. The luminescence was measured using the microplate reader. In case of FLT1 and FLT3, the plate was incubated for 45 min, then Kinase Detection reagent was added to each well and the plate was covered with aluminum foil and incubated at room temperature for another 45 min. The luminescence was measured immediately using the Tecan-spark ® reader. A "Blank" value was subtracted from all readings. All samples and controls were tested in duplicate, and IC 50 s were represented as ±SD. Sorafenib was used as a reference drug.

Cell Lines and Reagents
All cell lines were purchased from Nawah Scientific Company (Almokattam Mall, Cairo, Egypt) except A2780 and A2780CP, which were a kind gift from Prof. Jan Brabek (BIOCEV, Vestec u Prahy, Czech Republic). Cells were grown in DMEM medium (BioWhit-taker™ [78,79], the cancer or normal cells were seeded in a 96-well plate (100 µL/well). After overnight incubation at 37 • C and 5% CO 2 , the cells were incubated with 50 µM of each tested compound or DMSO (0.5% v/v). After 48 h of incubation, MTT (3-(4,5-dimethylthiazoyl)-2,5-diphenyl-tetrazolium bromide) 5 mg/mL Phosphate-Buffered Saline (PBS) was added, and the plate was incubated for 4 h. After that, acidified sodium dodecyl sulfate (SDS) solution (10% SDS containing 0.01N HCl in 1× PBS) was used to solubilize formazan crystals. The absorbance was measured after 14 h of incubation at λ 570-630 nm by a Biotek plate reader (Gen5™). To determine the IC 50 for the most active compound 12b, the cells were treated with serial dilutions of the 12b compound for 48 h, then the viability was detected using MTT as mentioned above.

Selectivity of 12b towards Cancer Cells
According to the procedures of Bézivin et al. (with minor modifications) [80][81][82], the selectivity of compound 12b against cancer cell lines was evaluated by seeding human skin fibroblast cells (HSF) at 5 × 10 4 cells/mL in a 96-well plate and incubating them overnight at 37 • C and 5% CO 2 . On the second day, serial dilution of the compound 12b was added. After 48 h of incubation, the viability of HSF cells was determined by MTT assay as mentioned above and the selectivity index (SI) was calculated according to Equation (1)

Cell-Cycle Analysis
According to Gray et al. [83], HepG2 cells were seeded in a 6-well plate (2 mL/well) and incubated at 37 • C and 5% CO 2 for 12 h. Cell synchronization at the G1 phase was performed by incubating the cells with serum-free medium for 12 h, then cells were treated with 5.5 or 11 µM of 12b. After 48 h of incubation, cells were washed twice with ice-cold 1× PBS and subjected to trypsinization. The detached cells were centrifuged at 500× g for 5 min at 4 • C. Then, the cell pellet was resuspended in ice-cold 1× PBS and fixed with absolute ethanol. After 2 h of incubation at −20 • C, the fixed cells were washed with ice-cold 1× PBS. The cell pellet was resuspended in propidium (PI)/RNase (BD Biosciences, BDB550825, Franklin Lakes, NJ, USA) and incubated for 30 min at room temperature (RT) in the dark, then the DNA content in each phase was measured by flow cytometry using BD FACS Calibur [84].

Apoptosis Analysis
Briefly, HepG2 cells (10 6 cells/mL) were treated either with 12b (5.5 or 11 µM) or with DMSO for 24 h, and then were collected after trypsinization and washed twice with 1× PBS. The induction of apoptosis was analyzed by staining cells with PE-Annexin-V [85] in the presence of PI and analyzed using the BD FACS Calibur Flow cytometer.

Detection of Nitric Oxide Level by DAF-FM Diacetate
HepG2 cells were collected, and the number was adjusted at 10 6 cells/test, and then they were stained with 5 µM DAF-FM™ diacetate (Invitrogen, Cat. No. D-23844) for 60 min at 37 • C and 5% CO 2 . After that, the stained cells were washed to remove the excess probe, and then were re-incubated for another 30 min. Following that, the stained cells were treated either with 12b (5.5 or 11 µM) or JS-K as a positive nitric oxide donor or DMSO as a negative control for another 1 h. The fluorescence was detected by flow cytometry using the FITC detection system in a BD flow Cytometer [86][87][88].

RT-PCR for Gene Expression
HepG2 cells were treated either with DMSO, 5.5 or 11 µM of 12b for 12 h, and then the total RNA was extracted by TRIzol™ Reagent (Invitrogen, Cat. No. 15596018). The RNA concentration and purity were evaluated using a Thermo Scientific™ NanoDrop™ One Microvolume UV-Vis Spectrophotometer. Four µg RNA was reverse transcribed to cDNA using the RevertAid First Strand cDNA Synthesis kit. The level of N-cadherin, E-cadherin, BCL-2, Bax, and GAPDH was detected using the SensiFAST™ SYBR ® High-ROX Kit with the StepOne™ Real-Time PCR System (Applied Biosystems™). The expression level was firstly normalized to the level of GAPDH, then calculated as 2 −(∆∆Ct) , where ∆C t 12b = (C t − C t GAPDH ); ∆C t DMSO = (C t DMSO − C t GAPDH ); ∆∆C t = (∆C t 12b − ∆C t DMSO ). C t is the cycle threshold. Primers were synthesized by Eurofins Scientific and are listed in Table 4 [58,89]. The cells were seeded (2.5 × 10 5 cells/mL) and treated as described above. After 48 h of incubation with compound 12b, cells were lysed by RIPA buffer containing 1x protease and phosphatase inhibitor cocktail (Thermo Scientific™, Waltham, MA, USA). The supernatant was collected after centrifugation at 4 • C and the Pierce BCA Protein Assay Kit (Thermo Scientific™, 23225, Waltham, MA, USA) was used for protein concentration detection. Proteins were then separated using SDS-PAGE at 100 V for 45 min. After that, the proteins were transferred to a nitrocellulose membrane. The membrane was subjected to a suitable primary antibody overnight. Then, an appropriate HRP-conjugated secondary antibody was used. GAPDH was used as a loading control. The membrane was developed using WesternBright ® ECL (Advansta Inc., K-12045-D20, San Jose, CA, USA), and the signal was detected using ChemiDoc XRS + (1708265, Bio-Rad Laboratories, Hercules, CA, USA). All antibodies were purchased from Cell Signaling Technology (Danvers, MA, USA) are listed in Table 5 [59]. The p38, MKK3, and JNK assays were performed according to the manufacturer's ELISA assay protocols [60][61][62]. Briefly, a 96-well plate was seeded with a cell density of 10,000-40,000 cells/well. Then, 50 µL of compound 12b or control (DMSO) was added to appropriate wells in triplicate. Following this, 50 µL of the antibody cocktail (in case of p38; Cat. No. ab221012, Abcam (Waltham, MA, USA) and MKK3; Cat. No. PEL-MKK3-S189-T, RayBiotech, Peachtree Corners, GA, USA) was added to each well. The plate was then sealed and incubated for 1-2.5 h at room temperature with gentle shaking. Each well was washed with 1X Wash Buffer PT (4 × 350 µL), and then the plate was inverted and tapped gently against clean paper towels to remove excess liquid. In case of JNK (Cat. No. ab176662, Abcam, Waltham, MA, USA), this washing step was applied before and after adding 1X rabbit anti-phospho-MKK3 antibody, then 100 µL of HRP-conjugated anti-rabbit IgG solution was added and the plate was incubated for 1 h at room temperature with gentle shaking followed by washing (4 × 350 µL). After that, 100 µL of TMB Substrate was added to each well and incubated for 15-30 min in the dark with gentle shaking. Finally, 100 µL of Stop Solution was added to each well and the plate was shaken for 1 min to mix well. The OD was recorded at 450 nm.

ADME/Toxicity Analysis
The web-based tool ADMETlab 2.0 was used for prediction of the physicochemical, pharmacokinetic, and drug-likeness properties of compound 12b, with respect to the reference drug sorafenib. ADMETlab 2.0 offers a straightforward approach for comprehensive, accurate, and efficient prediction of ADMET properties, using a high-quality database of 250,000 entries covering 53 endpoints, and a multi-task graph attention framework [69,70].

Docking Studies
Molecular modeling studies were performed using the OpenEye ® scientific software (2021.spr) [19,42,[90][91][92][93]. A library of pyrazolopyrimidines based on 1,2,5-oxadiaole-2-oxide were minimized using the MMFF94 minimization force field. OMEGA ® application was used for the generation of multi-conformers and the whole library was aligned against the reference drug sorafenib, according to shape and electrostatic potential, using the EON module. The VEGFR-2 (PDB code: 3WEZ) receptor was prepared using the Spruce ® module, after which the FRED ® application was used for docking and generating FRED Chemgauss4 scores. Then, the visualization tool, Vida ® application, was employed to display the potential binding interactions of the target compounds to the receptor active site. The Discovery Studio Visualization software was used to generate a 2D depiction of the protein-ligand interaction [67,68]. The docking report was generated using the Docking-Report application through the OpenEye Applications 2021.1.1 command line.

Conclusions
In this paper, the design, synthesis, and biological evaluation of a new series of pyrazolo [3,4-d]pyrimidines based on a phenylfuroxan scaffold were detailed. One-dose screening of the NCI-60 cell revealed that compounds 12a-c and 14a had the best MGI% among the tested compounds. The target compound 12b, as the most active one, showed better anticancer activity, with IC 50 values of 11.5, 11.6, and 13 µM against the HepG-2, A2780CP, and MDA-MB-231 cell lines, respectively, than the reference anticancer drug sorafenib. Furthermore, compound 12b showed VEGFR-2-inhibitory activity comparable to that of sorafenib. The mechanistic study revealed that compound 12b decreased the level of total ERK and its phosphorylated form and led to the down-regulation of the metastatic protein metalloproteinase MMP-9 and the over-expression of the cell-cycle inhibitors p21 and p27. Cell-cycle analysis demonstrated that compound 12b can arrest cells in the subG1 phase and, thus, induce apoptosis. Compound 12b was found to inhibit wound healing in the absence of serum, in comparison to DMSO-treated cells, as shown by a wound-healing assay, while quantitative RT-PCR for E-cadherin and N-cadherin indicated lower expression of neuronal N-cadherin and increased expression of epithelial E-cadherin, demonstrating the ability of compound 12b to suppress metastasis. Furthermore, 12b-treated HepG2 cells expressed a low level of anti-apoptotic BCL-2 and over-expressed pro-apoptotic Bax genes. Intracellular NO measurement indicated that the NO released from 12b was similar to that from the diazeniumdiolate JS-K reference drug, using DAF-FM DA as a fluorescence probe. The docking results indicated that the most active derivative, 12b, overlaid well with sorafenib, while the anticancer activity of the less active analogs could be enhanced through further modification of the central aryl linker, which projected towards the protein cavity. Furthermore, the ADMET study revealed similarities between the potential compound 12b and sorafenib, as well as its passing of the Pfizer acceptance criteria. Therefore, in this study, we presented a novel anticancer lead compound that is worth further investigation and activity improvement.