Synthesis, Molecular Docking Study, and Cytotoxicity Evaluation of Some Novel 1,3,4-Thiadiazole as Well as 1,3-Thiazole Derivatives Bearing a Pyridine Moiety

Pyridine, 1,3,4-thiadiazole, and 1,3-thiazole derivatives have various biological activities, such as antimicrobial, analgesic, anticonvulsant, and antitubercular, as well as other anticipated biological properties, including anticancer activity. The starting 1-(3-cyano-4,6-dimethyl-2-oxopyridin-1(2H)-yl)-3-phenylthiourea (2) was prepared and reacted with various hydrazonoyl halides 3a–h, α-haloketones 5a–d, 3-chloropentane-2,4-dione 7a and ethyl 2-chloro-3-oxobutanoate 7b, which afforded the 3-aryl-5-substituted 1,3,4-thiadiazoles 4a–h, 3-phenyl-4-arylthiazoles 6a–d and the 4-methyl-3- phenyl-5-substituted thiazoles 8a,b, respectively. The structures of the synthesized products were confirmed by spectral data. All of the compounds also showed remarkable anticancer activity against the cell line of human colon carcinoma (HTC-116) as well as hepatocellular carcinoma (HepG-2) compared with the Harmine as a reference under in vitro condition. 1,3,4-Thiadiazole 4h was found to be most promising and an excellent performer against both cancer cell lines (IC50 = 2.03 ± 0.72 and 2.17 ± 0.83 µM, respectively), better than the reference drug (IC50 = 2.40 ± 0.12 and 2.54 ± 0.82 µM, respectively). In order to check the binding modes of the above thiadiazole derivatives, molecular docking studies were performed that established a binding site with EGFR TK.


Introduction
Designing new, effective, selective, highly potent, although more tolerant, anticancer drugs through the identification of novel structures remains a considerable challenge for the researchers in the field of medicinal chemistry. Hybrid drug design has emerged during the past few years as a leading technique for the creation of innovative anticancer medicines that, in theory, can address many of the pharmacokinetic drawbacks of conventional anticancer medications [1]. Medical researchers have focused on pyridines, 1,3,4-thiadiazole and 1,3-thiazole systems that have led to somewhat more effective and promising results in recent years. To name just a few, the following pyridine-based small compounds have Under the influence of the findings mentioned above, continuous efforts were per formed to synthesize innovative anticancer compounds [13,[35][36][37][38][39][40][41][42][43][44][45][46]. The present repor aims to elaborate on the new series of thiadiazole-pyridines as well as thiazole-pyridines that might have cytotoxic effects via the inhibition of protein Epidermal Growth Factor Tyrosine Kinase receptor (EGFR TK), which plays an essential mediating role in cell pro liferation, angiogenesis, apoptosis, and metastatic spread compared with reference drugs. Under the influence of the findings mentioned above, continuous efforts were performed to synthesize innovative anticancer compounds [13,[35][36][37][38][39][40][41][42][43][44][45][46]. The present report aims to elaborate on the new series of thiadiazole-pyridines as well as thiazole-pyridines that might have cytotoxic effects via the inhibition of protein Epidermal Growth Factor Tyrosine Kinase receptor (EGFR TK), which plays an essential mediating role in cell proliferation, angiogenesis, apoptosis, and metastatic spread compared with reference drugs.
All the structures of products 4a-h were confirmed by elemental analyses followed by spectral data. In general, the 1 H-NMR spectra of 4a (see Supplementary Materials), taken as an example, showed a singlet (1H) at δ 6.39 ppm corresponding to the pyridine proton, three singlets at δ 2.01, 2.33 and 2.46 ppm corresponding to the three CH 3 groups and a multiplet δ 7.09-7.46 ppm corresponding to the five aromatic protons. The 13 C-NMR spectrum of 4a revealed two signals at δ = 25.6, 194.7 ppm, which is characteristic for the acetyl group (CH 3 C=O). The disappearance of the two NH absorption bands was also observed in the IR spectra because of the elimination of amine groups from the starting material 2. Moreover, the mass spectrum of the products 4a-h revealed a molecular ion peaks at the expected m/z values.
All the structures of products 4a-h were confirmed by elemental analyses followed by spectral data. In general, the 1 H-NMR spectra of 4a (see supplementary material), taken as an example, showed a singlet (1H) at δ 6.39 ppm corresponding to the pyridine proton, three singlets at δ 2.01, 2.33 and 2.46 ppm corresponding to the three CH3 groups and a multiplet δ 7.09-7.46 ppm corresponding to the five aromatic protons. The 13 C-NMR spectrum of 4a revealed two signals at δ = 25.6, 194.7 ppm, which is characteristic for the acetyl group (CH3C=O). The disappearance of the two NH absorption bands was also observed in the IR spectra because of the elimination of amine groups from the starting material 2. Moreover, the mass spectrum of the products 4a-h revealed a molecular ion peaks at the expected m/z values.

Scheme 1. Synthesis of thiadiazoles 4a-h.
In addition to this, the chemical reactivity of compound 2 with several of α-haloketones was also investigated to synthesize a series of novel thiazole derivatives. Accordingly, compound 2 reacted with α-haloketones 5a-d in the presence of TEA as a catalyst under refluxing condition using EtOH as a solvent that resulted a corresponding thiazoles 6a-d series, as shown in Scheme 2 (see Experimental).

Scheme 1. Synthesis of thiadiazoles 4a-h.
In addition to this, the chemical reactivity of compound 2 with several of α-haloketones was also investigated to synthesize a series of novel thiazole derivatives. Accordingly, compound 2 reacted with α-haloketones 5a-d in the presence of TEA as a catalyst under refluxing condition using EtOH as a solvent that resulted a corresponding thiazoles 6a-d series, as shown in Scheme 2 (see Experimental).
All the structures of the series of products 6a-e were also confirmed through the analytical followed by spectral data analysis (see Experimental). Compound 6c showed a typical singlet signal that appeared at δ 3.82, 6.55, and 6.80 ppm due to the OCH 3 group, pyridine-H5, and thiazole-H5, respectively, in the 1 H-NMR spectra. In addition to this, a multiplet is observed in the region: 7.01-7.47 ppm assignable to the nine aromatic hydrogens. On the other hand, the 13 C-NMR spectrum of 6c showed four signals at δ = 17.8, 21.2, 56.7 and 162.4 ppm characteristic for 2 Ar-CH 3 , Ar-OCH 3 and C=O groups, respectively, in addition to sixteen aromatic carbon signals in the range of 107.3-149.4 ppm. All the structures of the series of products 6a-e were also confirmed through the analytical followed by spectral data analysis (see Experimental). Compound 6c showed a typical singlet signal that appeared at δ 3.82, 6.55, and 6.80 ppm due to the OCH3 group, pyridine-H5, and thiazole-H5, respectively, in the 1 H-NMR spectra. In addition to this, a multiplet is observed in the region: 7.01-7.47 ppm assignable to the nine aromatic hydrogens. On the other hand, the 13 C-NMR spectrum of 6c showed four signals at δ = 17.8, 21.2, 56.7 and 162.4 ppm characteristic for 2 Ar-CH3, Ar-OCH3 and C=O groups, respectively, in addition to sixteen aromatic carbon signals in the range of 107.3-149.4ppm.
The structure of the isolated product 8 was inferred from its IR and 1 H-NMR spectral data and elemental analysis (see Experimental).

Anti-cancer Activity
The series of prepared compounds 4a-h and 6a-d were investigated against human colon carcinoma (HCT-116) followed by the hepatocellular carcinoma (HepG2) cell lines to obtain pharmacological activities using Harmine as a reference drug through colorimetric MTT assay under in vitro conditions. The survival curve was obtained by plotting the relation between the concentrations of the drugs against the surviving cells, resulting in the 50% inhibitory concentration (IC50). The anti-proliferative activity is also achieved through the expression of the mean IC50 by three independent experiments (µM) ± standard deviation calculated from three replicates. Table 1 and Figure 2 summarize the structure-and concentration-dependent anticancer activities of the series of compounds against HTC-116 cell lines. An in vitro inhibition activity shows a positive trend along with all tested compounds. Compounds such as 4c, 4d, 4f and 4g show comparable activity to that of Harmine (IC50 = 2.40 ± 0.12 µM) as a reference, whereas compound 4h demonstrates even better results compared with the same reference. A similar trend of results was also observed for the hepatocellular carcinoma (HepG2) cell line assay, where 4c, 4d, 4f, 4g and 4h show either comparable or improved inhibitory activity, with 6h (IC50 = 2.17 ± 0.83 µM) showing the maximum effect in comparison with the reference Harmine (IC50 = 2.54 ± 0.82 µM). Finally, compound 2 reacted with α-chloro compounds 7a,b under refluxing condition in the presence of EtOH as a solvent and TEA as a catalyst that resulted a single product, identified as 4,6-dimethyl-1-((4-methyl-3-phenyl-5-substitutedthiazol-2(3H)-ylidene) amino)4yridine-2(1H)-ones 8a,b, as outlined in Scheme 2.
The structure of the isolated product 8 was inferred from its IR and 1 H-NMR spectral data and elemental analysis (see Experimental).

Anti-Cancer Activity
The series of prepared compounds 4a-h and 6a-d were investigated against human colon carcinoma (HCT-116) followed by the hepatocellular carcinoma (HepG2) cell lines to obtain pharmacological activities using Harmine as a reference drug through colorimetric MTT assay under in vitro conditions. The survival curve was obtained by plotting the relation between the concentrations of the drugs against the surviving cells, resulting in the 50% inhibitory concentration (IC 50 ). The anti-proliferative activity is also achieved through the expression of the mean IC 50 by three independent experiments (µM) ± standard deviation calculated from three replicates. Table 1 and Figure 2 summarize the structure-and concentration-dependent anticancer activities of the series of compounds against HTC-116 cell lines. An in vitro inhibition activity shows a positive trend along with all tested compounds. Compounds such as 4c, 4d, 4f and 4g show comparable activity to that of Harmine (IC 50 = 2.40 ± 0.12 µM) as a reference, whereas compound 4h demonstrates even better results compared with the same reference. A similar trend of results was also observed for the hepatocellular carcinoma (HepG2) cell line assay, where 4c, 4d, 4f, 4g and 4h show either comparable or improved inhibitory activity, with 6h (IC 50

Docking Study for Cytotoxicity
The protein Epidermal Growth Factor Receptor Tyrosine Kinase Domain (EGFR TK) was selected for this study, where the ability to inhibit this receptor ultimately leads to the blockade of the growth pathway, giving a promising anti-cancer agent [48]. The lower binding energy resulting from the association of the compound with the targeted protein is an indication of a higher binding efficiency. The results of the docking protocol were validated by the re-docking of the co-crystallized ligand (W19) inside the active site of EGFR TK (Figures 4 and 5). Harmine was used in this study as an EGFR TK inhibitor. By comparing the binding affinity of different screened synthesized compounds with Harmine (ΔG of -7.1), it was found that compound 4h showed the best binding affinity with ΔG of −10.8, and 4b, 4c, 4d, 4e, 4f, 6a, 6b, and 6c showed binding activity ΔG −8.1 to −9.2. The screened compounds showed a possible interaction with EGFR TK active sites as depicted in Table 2     For thiadiazoles 4a-h: 4h (amidophenyl, has a phenyl ring along with electron withdrawing amido group resulting strongest activity) > 4g (strong electron withdrawing nitro group, increases activity) > 4f (with ester group along with one electron with drawing Cl atom) > 4d (acetyl group with electron withdrawing 2 Cl atom) > 4c (acety group with mild electron withdrawing one Cl atom) > 4a (with acetyl group with un-substituted phenyl group) > 4e (ester with methyl group) > 4b (acetyl group with methyl group, electron donating group decreases activity). Overall, electron releasing groups decrease the activity, whereas strong electron withdrawing groups increase th For thiadiazoles 4a-h: 4h (amidophenyl, has a phenyl ring along with electron withdrawing amido group resulting strongest activity) > 4g (strong electron withdrawing nitro group, increases activity) > 4f (with ester group along with one electron withdrawing Cl atom) > 4d (acetyl group with electron withdrawing 2 Cl atom) > 4c (acetyl group with mild electron withdrawing one Cl atom) > 4a (with acetyl group with un-substituted phenyl group) > 4e (ester with methyl group) > 4b (acetyl group with methyl group, electron donating group decreases activity). Overall, electron releasing groups decrease the activity, whereas strong electron withdrawing groups increase the activity. A selective high activity is observed, particularly with 4h, possibly due to the fact that 4h possess one extra phenyl ring connected with the pyridine group in the amido side, which significantly enhances its aromatic π-π interaction with the Phe, Tyr and Trp residues. This is in coherence with the harmine where one phenyl and one pyridine is fused with an indole group, resulting in the same kind of interaction with amino acid residues containing an aromatic group. In addition to this, a noteworthy mention would be, in 4h, the nitrogen in the amido group is engaged in a tautomeric structure, thereby restricting the electron releasing power of nitrogen in the phenyl ring. Such interactions are absent for the rest of thiadiazole derivatives 4a-g, where only electron withdrawing group is present in the lone available phenyl ring. This is in analogous with the thiazoles derivatives 6a-d, where only electron releasing groups are present with the single phenyl group available, resulting in less activity compared to thiadiazoles 4a-h with the only exception of 6d (an electron withdrawing Cl atom is present), where moderate anticancer activity is detected (Figure 3).

Docking Study for Cytotoxicity
The protein Epidermal Growth Factor Receptor Tyrosine Kinase Domain (EGFR TK) was selected for this study, where the ability to inhibit this receptor ultimately leads to the blockade of the growth pathway, giving a promising anti-cancer agent [48]. The lower binding energy resulting from the association of the compound with the targeted protein is an indication of a higher binding efficiency. The results of the docking protocol were validated by the re-docking of the co-crystallized ligand (W19) inside the active site of EGFR TK (Figures 4 and 5). Harmine was used in this study as an EGFR TK inhibitor. By comparing the binding affinity of different screened synthesized compounds with Harmine (ΔG of -7.1), it was found that compound 4h showed the best binding affinity with ΔG of −10. 8, and 4b, 4c, 4d, 4e, 4f, 6a, 6b, and 6c showed binding activity ΔG −8.1 to −9.2. The screened compounds showed a possible interaction with EGFR TK active sites as depicted in Table 2 and Figures 6-13.

Docking Study for Cytotoxicity
The protein Epidermal Growth Factor Receptor Tyrosine Kinase Domain (EGFR TK) was selected for this study, where the ability to inhibit this receptor ultimately leads to the blockade of the growth pathway, giving a promising anti-cancer agent [48]. The lower binding energy resulting from the association of the compound with the targeted protein is an indication of a higher binding efficiency. The results of the docking protocol were validated by the re-docking of the co-crystallized ligand (W19) inside the active site of EGFR TK (Figures 4 and 5). Harmine was used in this study as an EGFR TK inhibitor. By comparing the binding affinity of different screened synthesized compounds with Harmine (∆G of -7.1), it was found that compound 4h showed the best binding affinity with ∆G of −10. 8, and 4b, 4c, 4d, 4e, 4f, 6a, 6b, and 6c showed binding activity ∆G −8.1 to −9.2. The screened compounds showed a possible interaction with EGFR TK active sites as depicted in Table 2

Experimental
Elementar vario LIII CHNS analyzer (Elementar Analysensysteme GmbH, Langenselbold, Germany) is used to measure all elemental analysis. Electrothermal IA 9000 series Digital Melting Point Apparatus (Shanghai Jiahang Instruments Co., Jiading District, shanghai, China) was used to obtain melting points data. Shimadzu FTIR 8101 PC infrared spectrophotometers (Shimadzu Co., Kyoto, Japan) were used to record IR spectra data in KBr discs on Pye Unicam SP 3300. Varian Mercury VX-300 NMR spectrometer (Bruker Biospin, Karlsruhe, Germany) was used with the operating frequency of 300 MHz ( 1 H-NMR) in deuterated dimethylsulfoxide (DMSO-d6) solvent to record NMR spectra, where chemical shifts were related to the solvent used. Shimadzu GCeMS-QP1000 EX mass spectrometer(Shimadzu Co., Kyoto, Japan) was used to record

Experimental
Elementar vario LIII CHNS analyzer (Elementar Analysensysteme GmbH, Langenselbold, Germany) is used to measure all elemental analysis. Electrothermal IA 9000 series Digital Melting Point Apparatus (Shanghai Jiahang Instruments Co., Jiading District, Shanghai, China) was used to obtain melting points data. Shimadzu FTIR 8101 PC infrared spectrophotometers (Shimadzu Co., Kyoto, Japan) were used to record IR spectra data in KBr discs on Pye Unicam SP 3300. Varian Mercury VX-300 NMR spectrometer (Bruker Biospin, Karlsruhe, Germany) was used with the operating frequency of 300 MHz ( 1 H-NMR) in deuterated dimethylsulfoxide (DMSO-d6) solvent to record NMR spectra, where chemical shifts were related to the solvent used. Shimadzu GCeMS-QP1000 EX mass spectrometer (Shimadzu Co., Kyoto, Japan) was used to record mass spectra at 70 eV. The cytotoxicity of the prepared compounds was measured by the Regional Center for Mycology and Biotechnology in Al-Azhar University, Cairo, Egypt.

Synthesis of thiazoles 6a-d or 8a,b.
A mixture of 2 (0.298 g, 1 mmol) and α-haloketones 5a-d or 3-chloropentane-2,4-dione (7a) or ethyl 2-chloro-3-oxobutanoate (7b) (1 mmol) in DMF (15 mL) was refluxed for 4-6 h and was continuously monitored in TLC. The separation of the product was clearly observed during the course of the reaction. The resultant solid product was then filtered, washed several times with water, dried and recrystallized in the proper solvent to give the corresponding thiazoles 6a-d or 8a,b, respectively.

Cytotoxic Activity
The cytotoxicity of the newly synthesized series of compounds was studied against HCT-116 and HepG2 cells using the MTT assay through an incubation period of 24 h [49,50].
Mammalian cell line: HCT-116 and HepG2 cells were collected from VACSERA Tissue Culture Unit, Cairo, Egypt.

Docking Method
The MOE 2019.012 suite (Chemical Computing Group ULC, Montreal, Canada) [51] was applied in order to carry all docking studies for the newly synthesized compounds to suggest their plausible mechanism of action as the protein Epidermal Growth Factor Receptor Tyrosine Kinase Domain (EGFR TK) inhibitors by evaluating their binding grooves and modes to compare with Harmine as a reference standard.
The newly prepared derivatives were placed into the MOE window, where they were treated to partial charge addition and energy minimization [52,53]. The produced compounds also were placed into a single database with the Harmine and saved as an MDB file, which was then uploaded to the ligand icon during the docking process. The Protein Data Bank was used to generate an X-ray of the targeted Epidermal Growth Factor Receptor Tyrosine Kinase Domain (EGFR TK) 3W33. Available online: https://www.rcsb.org/ structure/3W33 (accessed on 17 July 2022). [48]. Furthermore, it was readied for the docking process by following the previously detailed stages [54,55]. Additionally, the downloaded protein was error-corrected, 3D hydrogen-loaded and energy-minimized [56,57].
In a general docking procedure, the newly prepared derivatives were substituted for the ligand site. After modifying the default program requirements previously stated, the co-crystallized ligand site was selected as the docking site, and the docking process was started [58]. In a nutshell, the docking site was chosen using the dummy atoms method [59]. The placement and scoring procedures, respectively, Triangle matcher and London dG, were chosen. The stiff receptor was used as the scoring method, and the GBVI/WSA dG was used as the refining method, respectively, to select the top 10 poses out of a total of 100 poses for each docked molecule [60,61]. For further research, the optimal pose for each ligand with the highest score, binding mode and RMSD value was chosen. It is important to note that the applied MOE program underwent the first step of program verification by docking Harmine to its ligand binding of the prepared Target [62,63]. By obtaining a low RMSD value (1.43) between the newly created compounds with docked Harmine, a valid performance was demonstrated.

Conclusions
In this paper, two new series of aryl substituted novel pyridine-1,3,4-thiadiazoles, and pyridine-thiazoles were synthesized starting with 1-(3-cyano-4,6-dimethyl-2-oxopyridin-1(2H)-yl)-3-phenylthiourea and several available reagents. All the structures were confirmed through elemental and spectral analysis, where the plausible mechanistic approach for their formation was discussed. All the prepared derivatives showed effectiveness towards the inhibition of human colon carcinoma (HCT-116) as well as hepatocellular carcinoma (HepG2) cell lines through in vitro evaluation and an in silico docking study. From the obtained results, compound 4h (amidophenyl has a phenyl ring along with the electron withdrawing amido group resulting in the strongest activity) was found to be the strongest and most effective with 2.03 ± 0.72 µM against HCT-116, contributing its activity through a variety of interactions, such as hydrophobic interaction, hydrogen bonding in addition to aromatic stacking interactions with selected target (EGFR TK) pockets, compared with Harmine as a reference drug. A detailed analysis of all the series confirmed the electron withdrawing group present in the aryl substitution, resulting in the enhancement of anticancer activity, which could be promising for the future generation of new efficient anticancer drugs based on 1,3,5-thiadiazole and 1,3-thiazole derivatives.

Data Availability Statement:
The data presented in this study are available on request from corresponding author.