Synthesis and Biological Evaluation of Imadazo[1,2-a]pyrazines as Anticancer and Antiviral Agents through Inhibition of CDK9 and Human Coronavirus

In this work, novel imadazo[1,2-a]pyrazine derivatives were synthesized and evaluated as CDK9 inhibitors. The results of CDK9 assay showed that the derivatives with pyridin-4-yl in position 2 and benzyl in position 3 of imadazo[1,2-a]pyrazine 3c displayed the most potent CDK9 inhibitory activity with IC50 of 0.16 µM. The anti-proliferative effect of the new compounds was examined against breast cancer (MCF7), colorectal cancer (HCT116), and chronic myelogenous leukaemia (K652) cell lines. The data of MTT assay showed that the cytotoxic effect of the inhibitors is correlated to their inhibitory activity against CDK9. Compound 3c exhibited the most potent cytotoxicity effect with average IC50s of three cell lines of 6.66 µM. The drug likeness properties of 3c were predicated in silico and demonstrated that 3c have reasonable physiochemical and pharmacokinetic properties. Selected derivatives were assessed in antiviral assay against human coronavirus 229E. The results of this assay showed that the derivative with pyridin-4-yl in position 2 and cyclohexyl in position 3 of imadazo[1,2-a]pyrazine 3b exhibited the most potent anti-coronaviral activity with IC50 of 56.96 µM and selectivity index of 7.14. The target predication result revealed that 3b showed high affinity to protease enzyme. Docking studies of 3b with COVID-19 main protease was conducted and showed good binding affinity, which confirmed the in vitro assay data.


Introduction
Cyclin-dependent kinase (CDK) is a family of enzymes that dimerize with specific proteins called cyclins. These CDK-cyclin heterodimers have crucial roles in gene transcription and cell cycle progression [1]. CDK9 is a member of the CDK family, which forms a heterodimer complex with cyclin T. This dimer is known as Positive Transaction Elongation Factor b (p-TEFB), which stimulates the transcriptional elongation by phosphorylation of RNA polymerase II at Ser2 and Ser5 [2,3].
There is sufficient evidence that supports the role of CDK9 as a molecular target for treatment of a number of malignancies including breast cancer, colorectal cancer, and acute myeloid leukaemia (AML) [4][5][6][7][8][9]. CDK9 conducts its oncogenic function through the regulation of anti-apoptotic agents such as Mc-1 and Myc and the promotion of cell proliferation [10,11].
Large pharmaceutical companies have developed a huge number of molecules as CDK9 inhibitors, and most of them are small molecules with good drug likeness properties and potent antitumor activity. Several of these were proceeded to clinical trials as anticancer drugs to treat various types of cancer [12][13][14][15][16][17][18][19][20]. The main chemical scaffolds of clinical and In 2019, the imidazo[1,2-a]pyrazines 9 and 10 were identified by virtual screening from natural resources as CDK9 inhibitors with IC50 values of 7.88 and 5.12 µM, respectively, and evaluated as anticancer agents against breast cancer [27]. In a separate patent, compounds with 2-phenylimidazo[1,2-a]pyrazin-3-amine scaffold were reported as antiviral agent against influenza A virus such as 11 and 12 [28]. Compounds with similar structures such as 13 and 14 were reported in another patent as inhibitors of the transcriptional factor SALL4 and evaluated as anticancer agents against AML cell lines ( Figure 2) [29].
In this study, imidazo [1,2-a]pyrazine was used as a scaffold for the preparation of CDK9 inhibitors and antiviral agents. Apart from the patent that reported the identification of two compounds having this scaffold by virtual screening from natural resources In 2019, the imidazo[1,2-a]pyrazines 9 and 10 were identified by virtual screening from natural resources as CDK9 inhibitors with IC 50 values of 7.88 and 5.12 µM, respectively, and evaluated as anticancer agents against breast cancer [27]. In a separate patent, compounds with 2-phenylimidazo[1,2-a]pyrazin-3-amine scaffold were reported as antiviral agent against influenza A virus such as 11 and 12 [28]. Compounds with similar structures such as 13 and 14 were reported in another patent as inhibitors of the transcriptional factor SALL4 and evaluated as anticancer agents against AML cell lines ( Figure 2) [29]. as inhibitors of CDK9 [27], no other published data, to the best of our knowledge, reports the development of imidazo [1,2-a]pyrazine as CDK9 inhibitors. Therefore, the rationale of choosing this scaffold is to develop a new class of CDK9 inhibitors bearing a different scaffold from those already extensively studied and currently available in the literature. The aim of this study is to develop novel 2-phenylimidazo[1,2-a]pyrazin-3-amine as CDK9 inhibitors and examine their cytotoxicity as anticancer agents against a panel of cancer cell lines, study their binding mode, and evaluate of their drug-likeness properties in silico. Pathogenic human coronaviruses are species that cause respiratory diseases, appearing early in 2002 by the identification of SARS-CoV and recently by SARS-CoV-2, which caused the COVID-19 pandemic in 2019 [30]. There are a few pathogenic species of human coronaviruses such as 229E, OC43, and NL63 that cause mild diseases but share a similar homology sequence [31]. In this study, the antiviral activity of the new compounds was also evaluated against human coronavirus.

Chemistry
The synthetic approach that was used to synthesize the novel compounds is described in Scheme 1. In this work, the new compounds were synthesized via one-pot multicomponent Groebke-Blackburn-Bienaymé reaction by reacting 2-aminopyrazine with various isocyanides and various aldehydes in the presence of scandium(III)-trifluoromethanesulfonate as catalyst and DCM:MeOH (3:1) as a solvent system. The reaction was carried out in a microwave under 150°C temperature for 30 min. The products were purified by column chromatography using ethyl acetate: hexane.
The first serial of compounds 1(a-d) have phenyl in position 2 of imadazo[1,2-a]pyrazine and various amines in position 3 including aliphatic, aromatic, cyclic, and acyclic amines. The yield of this series ranges from 84 to 96%. The second series 2(a-d) with phenyl-3,4 diol in position 2 of imadazo[1,2-a]pyrazine and the yield ranged from 75 to 91%. In the third series 3(a-d), pyridin-4-yl is placed on position 2 of imadazo[1,2-a]pyrazine, and this series obtained an excellent yield 88-97%. The fourth series 4(a-d) had thiophen- In this study, imidazo[1,2-a]pyrazine was used as a scaffold for the preparation of CDK9 inhibitors and antiviral agents. Apart from the patent that reported the identification of two compounds having this scaffold by virtual screening from natural resources as inhibitors of CDK9 [27], no other published data, to the best of our knowledge, reports the development of imidazo[1,2-a]pyrazine as CDK9 inhibitors. Therefore, the rationale of choosing this scaffold is to develop a new class of CDK9 inhibitors bearing a different scaffold from those already extensively studied and currently available in the literature.
The aim of this study is to develop novel 2-phenylimidazo[1,2-a]pyrazin-3-amine as CDK9 inhibitors and examine their cytotoxicity as anticancer agents against a panel of cancer cell lines, study their binding mode, and evaluate of their drug-likeness properties in silico. Pathogenic human coronaviruses are species that cause respiratory diseases, appearing early in 2002 by the identification of SARS-CoV and recently by SARS-CoV-2, which caused the COVID-19 pandemic in 2019 [30]. There are a few pathogenic species of human coronaviruses such as 229E, OC43, and NL63 that cause mild diseases but share a similar homology sequence [31]. In this study, the antiviral activity of the new compounds was also evaluated against human coronavirus.

Chemistry
The synthetic approach that was used to synthesize the novel compounds is described in Scheme 1. In this work, the new compounds were synthesized via one-pot multicomponent Groebke-Blackburn-Bienaymé reaction by reacting 2-aminopyrazine with various isocyanides and various aldehydes in the presence of scandium(III)-trifluoromethanesulfonate as catalyst and DCM:MeOH (3:1) as a solvent system. The reaction was carried out in a microwave under 150 • C temperature for 30 min. The products were purified by column chromatography using ethyl acetate: hexane. pound 2a showed the appearance of six aromatic signals at the region from 145 to 116 ppm, representing imadazopyrazine carbons, in addition to another six aromatic signals at the same region, representing the phenyl-3,4 diol ring at position 2 of imadazo[1,2-a]pyrazine. Two aliphatic carbon signals appeared at 56.63 and 30.54 ppm, representing the tert-butyl group attached to the amino group at position 3 of imadazo[1,2-a]pyrazine confirm the structure of 2a ( Figure 4). The first serial of compounds 1(a-d) have phenyl in position 2 of imadazo[1,2-a]pyrazine and various amines in position 3 including aliphatic, aromatic, cyclic, and acyclic amines. The yield of this series ranges from 84 to 96%. The second series 2(a-d) with phenyl-3,4 diol in position 2 of imadazo[1,2-a]pyrazine and the yield ranged from 75 to 91%. In the third series 3(a-d), pyridin-4-yl is placed on position 2 of imadazo[1,2-a]pyrazine, and this series obtained an excellent yield 88-97%. The fourth series 4(a-d) had thiophen-3-yl in position 2 of imadazo[1,2-a]pyrazine. These compounds were achieved in yield from 85% to 90%.
The formation of imadazo[1,2-a]pyrazine core was confirmed by 1 H and 13 C NMR spectra by appearances of three signals in aromatic region corresponding to pyrazine ring methylene protons in 1 H NMR spectrum and six signals in the aromatic region, corresponding to imadazopyrazine carbons in 13

CDK9 Inhibitory Activity
Cyclin-dependent kinase-9 (CDK9) enzyme forms a heterodimer with cyclin T. This complex is known as Positive Transaction Elongation Factor b (p-TEFB), which has a crucial role as elongation factor during RNA transcription [2,3]. CDK9 phosphorylates C-terminus of RNA polymerase II [9,26].
There is compelling evidence to support the role of CDK9 as a molecular target for treatment of a number of malignancies via the regulation of anti-apoptotic proteins and the promotion of cancer initiation and progression [4][5][6][32][33][34].
In this study, all the newly synthesized compounds were assessed against CDK9 in a biochemical assay, and dinaciclib was used as a standard compound. Dinaciclib is a well-established Pan-CDKs inhibitor with low nM potency against CDK9 (literature IC50 of 4 nM [35,36]). The results are presented in Table 1 as IC50 (µM).
The data indicate that compounds with pyridin-4-yl at position 2 of imadazo[1,2a]pyrazine 3(a-d) displayed the most potent CDK9 inhibitory activity with sub µM IC50. Among them, 3c, a compound with benzyl amine in position 3 of imadazo[1,2-a]pyrazine, was found as the most active derivative, with IC50 of 0.16 µM.

CDK9 Inhibitory Activity
Cyclin-dependent kinase-9 (CDK9) enzyme forms a heterodimer with cyclin T. This complex is known as Positive Transaction Elongation Factor b (p-TEFB), which has a crucial role as elongation factor during RNA transcription [2,3]. CDK9 phosphorylates C-terminus of RNA polymerase II [9,26].
There is compelling evidence to support the role of CDK9 as a molecular target for treatment of a number of malignancies via the regulation of anti-apoptotic proteins and the promotion of cancer initiation and progression [4][5][6][32][33][34].
In this study, all the newly synthesized compounds were assessed against CDK9 in a biochemical assay, and dinaciclib was used as a standard compound. Dinaciclib is a well-established Pan-CDKs inhibitor with low nM potency against CDK9 (literature IC50 of 4 nM [35,36]). The results are presented in Table 1 as IC50 (µM).
The data indicate that compounds with pyridin-4-yl at position 2 of imadazo[1,2a]pyrazine 3(a-d) displayed the most potent CDK9 inhibitory activity with sub µM IC50. Among them, 3c, a compound with benzyl amine in position 3 of imadazo[1,2-a]pyrazine, was found as the most active derivative, with IC50 of 0.16 µM.

CDK9 Inhibitory Activity
Cyclin-dependent kinase-9 (CDK9) enzyme forms a heterodimer with cyclin T. This complex is known as Positive Transaction Elongation Factor b (p-TEFB), which has a crucial role as elongation factor during RNA transcription [2,3]. CDK9 phosphorylates C-terminus of RNA polymerase II [9,26].
There is compelling evidence to support the role of CDK9 as a molecular target for treatment of a number of malignancies via the regulation of anti-apoptotic proteins and the promotion of cancer initiation and progression [4][5][6][32][33][34].
In this study, all the newly synthesized compounds were assessed against CDK9 in a biochemical assay, and dinaciclib was used as a standard compound. Dinaciclib is a well-established Pan-CDKs inhibitor with low nM potency against CDK9 (literature IC 50 of 4 nM [35,36]). The results are presented in Table 1 as IC 50 (µM).
The data indicate that compounds with pyridin-4-yl at position 2 of imadazo[1,2a]pyrazine 3(a-d) displayed the most potent CDK9 inhibitory activity with sub µM IC 50 . Among them, 3c, a compound with benzyl amine in position 3 of imadazo[1,2-a]pyrazine, was found as the most active derivative, with IC 50 of 0.16 µM. Furthermore, the other derivatives with benzyl amine in position 3 of imadazo[1,2a]pyrazine, 2c and 4c were also exhibited good CDK9 inhibitory activity with IC 50 of 0.31 and 0.71 µM, respectively.

Cytotoxicity Activity
All synthesized compounds were evaluated for their cytotoxicity against breast cancer (MCF7), colorectal cancer (HCT116), and chronic myelogenous leukaemia (K652) cell lines using staurosporine as positive control in MTT assay. The results are shown in Table 2 as IC 50 (µM).
The data of the MTT assay showed that the cytotoxic effect of the inhibitors is correlated to their inhibitory activity against CDK9. The most potent compounds in the biochemical assay against CDK9, for example, 3c, 3a, and 1b (IC 50 of 0.16, 0.26, and 0.25 µM, respectively), also exhibited the greatest cytotoxicity effect in the anti-proliferation assay (average IC 50 of 6.66, 17.66, and 18.24 µM, respectively). Moreover, the compounds with weak CDK9 inhibitory activity in the biochemical assay, for example, 1a, 2d, and 4d (IC 50 of 1.67, 2.03, and 8.5 µM, respectively), also exhibited a weak cytotoxic effect in the MTT assay (average IC 50 of 437.39, 126.11, and 18.24 µM, respectively). This could indicate that the anticancer mechanism of these compounds is through the inhibition of the CDK9.
The compounds that showed the greatest cytotoxic effect on cancer cell lines were further evaluated for their cytotoxicity on normal non-cancerous cells FHC. The results are shown in Table 3. Compounds 1b, 3a, 3c, and 4a showed IC 50 values of 122.51, 195.78, 97.19, and 76.07 µM on normal cells, respectively, and average selectivity over three cancer cells of 43, 41, 15, and 3.5 fold, respectively. This indicates that the most promising compounds displayed selective cytotoxicity effect on cancer cell lines.

Molecular Docking Study into CDK9 Active Site
A docking simulation of 2c, 3a, 3c and 4a into the ATP binding site of CDK9 was carried out with AutoDock Vina embedded in PyRx 0.8 software using CDK9 co-crystalized with CR8 inhibitor (PDB:3LQ5). The validation of the docking test was performed by re-docking of the co-crystalized ligand (CR8). It revealed a docking score of −8.7 kcal/mol with an RMSD value of 1.7 Å. The key hydrogen bond interactions between the Cys106 and aminopurine ring of CR8 were obtained.
Generally, all tested compounds were able to fit into ATP binding pocket and interact well with key amino acid residues in the CDK9 binding site. (Figure 5, Table 4). All docked compounds showed a hydrogen bond interaction between N1 of imadazo[1,2-a]pyrazine ring and Ile25. Additionally, in all compounds, the alkyl side chain attached to the amino group in position 3 of imadazo[1,2-a]pyrazine occupied the hydrophobic pocket in the CDK9 binding site and formed several hydrophobic interactions with Phe103, Val29, Ala166, Ala46, and Leu156.   Compounds 3a and 4a showed similar binding pose, were N6 of imadazo[1,2-a]pyrazine interacted with Gly28 of G-rich loop, and the pyridin-4-yl of 3a and thiophen-3-yl of 4a in position 2 of imadazo[1,2-a]pyrazine bind to the hinge region by forming a hydrogen bond with Cys106 ( Figure 5).
While compounds 2c and 4c showed a flipped orientation compared to compounds 3a and 4a, a pyrazine ring occupies the hinge region by forming several interactions with Glu107, Asp109, and Cys106, while the dihydroxy phenyl ring of 2c in position 2 of imadazo[1,2-a]pyrazine interacts with G-rich loop and forms two hydrogen bonds with Asp167.

Anti-Coronaviral Inhibitory Activity
The anti-coronaviral activity of the new compounds was assessed using a cytopathic effect (CPE) inhibition assay against human coronavirus 229E. Ribavirin (well-established antiviral agent) was used as standard in this assay.
The results showed that compound 3b demonstrated good anti-coronaviral activity with an IC 50 concentration of 56.96 µM; it also exhibited no cytotoxic effect on target cells at concentrations that produced an antiviral effect of 0.1 to 1000 µg/mL, with a 50% cellular cytotoxicity CC 50 concentration of 406.86 µM and a selectivity index of 7.14. Therefore, 3b showed better antiviral activity and selectivity than standard drug ribavirin (Table 5).

Anti-Coronavirus Target Identification and Docking Studies
The predication of target for the viral activity of our compound was carried out using SwissTarget [37]. The data of the target predication of compound 3b are presented in Figure 6. The result showed that 3b showed high affinity to protease enzyme as target. Therefore, the potent antiviral activity of 3b demonstrated in antiviral assay could be due to the inhibition of coronavirus protease.
± standard deviation of three individual experiments.

Anti-Coronavirus Target Identification and Docking Studie
The predication of target for the viral activity of our compoun SwissTarget [37]. The data of the target predication of compound 3b 6. The result showed that 3b showed high affinity to protease enzy the potent antiviral activity of 3b demonstrated in antiviral assay hibition of coronavirus protease. Human coronavirus (HCoV-229E) shares high homology seq navirus SARS-CoV-2 but is less pathogenic [38]. Therefore, SAR acute and lethal respiratory disease in 2019 (COVID-19) was selec ing studies in this study [39].
Docking studies of 3b into SARS-CoV-2 main protease enzym AutoDock Vina embedded in PyRx 0.8 software using COVID-19 talized with X77 inhibitor (PDB:6W63). The validation of the dock by re-docking the co-crystalized ligand (X77). It revealed a dockin with RMSD value of 1.7 Å. The key hydrogen bond interactions b Gly143 residues and X77 inhibitor were obtained.
The results of the docking studies displayed a hydrogen bo N-pyridine of 3b and Cys44. The imadazo[1,2-a]pyrazine of 3b Human coronavirus (HCoV-229E) shares high homology sequence similarity to coronavirus SARS-CoV-2 but is less pathogenic [38]. Therefore, SARS-CoV-2, which caused acute and lethal respiratory disease in 2019 (COVID-19) was selected for molecular docking studies in this study [39].
Docking studies of 3b into SARS-CoV-2 main protease enzyme were carried out with AutoDock Vina embedded in PyRx 0.8 software using COVID-19 main protease co-crystalized with X77 inhibitor (PDB:6W63). The validation of the docking test was performed by re-docking the co-crystalized ligand (X77). It revealed a docking score of −9.7 kcal/mol with RMSD value of 1.7 Å. The key hydrogen bond interactions between the Glu166 and Gly143 residues and X77 inhibitor were obtained.
The results of the docking studies displayed a hydrogen bond interaction between N-pyridine of 3b and Cys44. The imadazo[1,2-a]pyrazine of 3b interacted with Phe140 residue via a hydrogen bond. The pyridyl and cyclohexyl rings of 3b formed several hydrophobic interactions with Met165, Met49, His41, and Cys44. The imadazo[1,2-a]pyrazine also formed hydrophobic interactions with Cy145, Glu166, and His41, as shown in Figure 7 and Table 6.

Physiochemical and Drug Likeness Properties
Molecular and physiochemical parameters of compounds 1b, 3c, and 4a were calculated using SwissADMET [40] and presented in Table 7. All derivatives exhibited good drug-likeness properties and followed Lipinski role of five with zero violation. Their molecular weight is less than 500; their lipophilicity is within the optimal range (log P between 1 and 5); their polarity is within the optimal level (tPSA < 90 Å 2 ); and the number of hydrogen acceptors/donor < 5. Therefore, our compounds are expected to be orally bioavailable. Additionally, they demonstrated good aqueous solubility (Log S < −5).   Molecular and physiochemical parameters of compounds 1b, 3c, and 4a were calculated using SwissADMET [40] and presented in Table 7. All derivatives exhibited good drug-likeness properties and followed Lipinski role of five with zero violation. Their molecular weight is less than 500; their lipophilicity is within the optimal range (log P between 1 and 5); their polarity is within the optimal level (tPSA < 90 Å 2 ); and the number of hydrogen acceptors/donor < 5. Therefore, our compounds are expected to be orally bioavailable. Additionally, they demonstrated good aqueous solubility (Log S < −5).

Pharmacokinetic and Toxicity Studies
Pharmacokinetic and toxicity studies of compounds 1b, 3c, and 4a were estimated in silico using the ADMETLab platform [41], and the data of these studies are shown in Table 8. Our compounds demonstrated excellent oral absorption; they showed good caco-2 permeability and low papp efflux (<5.15), none of them are expected to be a substrate or inhibitor of p-glycoprotein, and they all tested positive for human intestine absorption and had a high score of oral bioavailability. They also showed acceptable distribution profiles, with a volume of distribution within the optimal range (0.04-20 L/kg) and reasonable binding to plasma protein within the accepted level (<90%), but all studied derivatives showed high probability to cross the blood-brain barrier. With respect to metabolism, they proved to be inhibitors of CYP3A4, but none of them are expected to be a substrate of the same CYP isoform. They showed excreted T 1/2 > 1.5 h and clearance lower than 2.5 min/kg/L, which is considered an accepted excretion profile during the drug development process. On the other hand, they showed a suboptimal toxicity profile by demonstrating high probability for the inhibition of hERG, high probability for acute toxicity, and hepatic injury, but they did not cause skin sensitivity at the maximum recommended dose.

Chemistry
Chemicals such as 2-aminopyrazine, aldehydes, and isocyanides were obtained from TCI and Alderich Chemical Companies. The measurement of melting points was conducted Pharmaceuticals 2022, 15, 859 13 of 18 using the electrothermal melting point apparatus (UK), and the results were not corrected. NMR spectra data were obtained using Bruker Ascend 700 NMR spectrometer (Fällanden, Switzerland) at 700.17 MHz for proton NMR and 176.08 MHz for carbon NMR, and DMSOd 6 was used as a solvent in all samples. Thin-layer chromatography (TLC) technique was used as a monitoring technique to check the reaction progression using silica-gel-precoated aluminum sheets (60 F254, Merck) and visualize it under UV at 365 and 254 nm wavelength. The reactions were carried out in a microwave using Biotage, Initiator Sweden AB using vial size 2.0-5.0 mL and normal absorption level.

In Vitro CDK9 Assay
The CDK9 assay kit was purchased from BPS Biosciences (San Diego, CA, USA) Catalogue# 79628, and the anti-CDK effect was measured following the manufacturer's instructions using kinase-GLO ® MAX as a detection component. The luminescence was measured using a microplate reader. The results were analyzed using GraphPad Prism 5.0, and the kinase inhibitory activity was expressed as IC 50 for three individual experiments. In the assay, dinaciclib was used as a positive control, and DMSO was used as negative control.

MTT Cytotoxicity Assay
The in vitro cytotoxicity of our compounds was determined using MTT assay. Cancer cells were exposed to different concentrations of the compounds, ranging from 0.1 to 10 µM for 48 h. The treated cells were then incubated with 10% v/v reconstituted MTT for 3 h before reading the plates on Wallac Victor2 1420 multilabel counter in fluorescence mode at a wavelength of 460/590 nm (ex/em). The results were analyzed using GraphPad Prism 5.0, and cytotoxicity was expressed as IC 50 for three individual experiments. In the assay, staurosporine was used as a positive control, and DMSO was used as a negative control.

Antiviral Assay
The crystal violet method was used to evaluate antiviral activity and cytotoxicity assays using cytopathic (CPE) inhibition effect as described by Choi et al., 2009 [42]. In this assay, the low pathogenic coronavirus (229E) was used. The percentage of antiviral activities of the test's compounds were calculated according to Pauwels et al. (1988) [43] using the following equation: antiviral activity = [(mean optical density of cell controls − mean optical density of virus controls)/(optical density of test − mean optical density of virus controls)] × 100%. Based on these results, the 50% CPE inhibitory dose (IC50) was calculated. The results were analyzed using GraphPad Prism 5.0, and cytotoxicity was expressed as IC 50 for three individual experiments. In the assay, ribavirin was used as a positive control, and DMSO was used as a negative control.

Docking Studies
The ligands used in these studies (compounds 2c, 3a, 3b, 3c, and 4a) were prepared using ChemBoiDraw Ultra 14.0, and each ligand file was saved as PDB format. The crystal structure of the proteins was obtained from Protein Data Bank (CDK9 co-crystalized with CR8 inhibitor; PDB:3LQ5 and COVID-19 co-crystalized with X77 inhibitor; PDB:6W63), and each protein file was downloaded as a PDB file. Then, the preparation of the protein was performed using Discovery studio by selecting one subunit of the protein and inverted the selection and remove other molecules. The clean protein was then saved as a PDB file. The co-crystalized ligand was also saved as a PDB file and used in the docking studies as a comparative reference in determining the binding pocket. Then, the AutoDock software was used to introduce polar hydrogens to the protein, and it was saved as PDBQT file. The docking studies was conducted using AutoDock vina and grid box embedded in PyRx software. The docking results were analyzed using Discovery Studio by the visualization of the ligand-receptor interactions.

Conclusions
In summary, novel imadazo[1,2-a]pyrazines were synthesized and evaluated for their anticancer and antiviral activities. All the newly synthesized compounds were evaluated in a CDK9 inhibitory assay as the molecular target for their anticancer activity. The data of CDK9 biochemical assay showed that the derivatives with pyridin-4-yl in position 2 of imadazo[1,2-a]pyrazine 3(a-d) displayed the most potent CDK9 inhibitory activity with IC 50 range 0.16-0.8 µM. The cytotoxicity effect of all synthesized compounds was evaluated in an MTT assay against breast cancer (MCF7), colorectal cancer (HCT116), and chronic myelogenous leukaemia (K652) cell lines. The results of this assay showed that the cytotoxic effect of the inhibitors is correlated to their inhibitory activity against CDK9. Compounds 3c, 3a, and 1b showed the most potent cytotoxicity effect, with average IC 50 s of three cell lines of 0.16, 0.26, and 0.25 µM, respectively. The new compounds demonstrated reasonable physiochemical and predicated pharmacokinetic properties including lipophilicity, aqueous solubility, absorption, distribution, and elimination. The most promising agents were evaluated in an antiviral assay against human coronavirus (HCoV-229E). In this assay, 3b displayed the most potent antiviral activity, with an IC 50 of 56.96 µM and selectivity index of 7.14, better than the standard drug ribavirin. The target predication result revealed that 3b showed high affinity to protease enzyme. Therefore, the potent anti-coronaviral in in vitro assay could be contributed to the inhibition of protease enzyme. Docking studies of 3b with COVID-19 main protease were conducted and showed good binding affinity, which supports the in vitro results.