One-Pot Synthesis and Molecular Modeling Studies of New Bioactive Spiro-Oxindoles Based on Uracil Derivatives as SARS-CoV-2 Inhibitors Targeting RNA Polymerase and Spike Glycoprotein

The first outbreak in Wuhan, China, in December 2019 was reported about severe acute coronaviral syndrome 2 (SARS-CoV-2). The global coronavirus disease 2019 (COVID-19) pandemic in 2020 resulted in an extremely high potential for dissemination. No drugs are validated in large-scale studies for significant effectiveness in the clinical treatment of COVID-19 patients, despite the worsening trends of COVID-19. This study aims to design a simple and efficient cyclo-condensation reaction of 6-aminouracil derivatives 2a–e and isatin derivatives 1a–c to synthesize spiro-oxindoles 3a–d, 4a–e, and 5a–e. All compounds were tested in vitro against the SARS-CoV-2. Four spiro[indoline-3,5′-pyrido[2,3-d:6,5-d’]dipyrimidine derivatives 3a, 4b, 4d, and 4e showed high activities against the SARS-CoV-2 in plaque reduction assay and were subjected to further RNA-dependent-RNA-polymerase (RdRp) and spike glycoprotein inhibition assay investigations. The four compounds exhibited potent inhibitory activity ranging from 40.23 ± 0.09 to 44.90 ± 0.08 nM and 40.27 ± 0.17 to 44.83 ± 0.16 nM, respectively, when compared with chloroquine as a reference standard, which showed 45 ± 0.02 and 45 ± 0.06 nM against RdRp and spike glycoprotein, respectively. The computational study involving the docking studies of the binding mode inside two proteins ((RdRp) (PDB: 6m71), and (SGp) (PDB: 6VXX)) and geometrical optimization used to generate some molecular parameters were performed for the most active hybrids.


Introduction
Coronaviruses (CoV) are a type of virus that can cause mild to severe respiratory distress symptoms of cough, high fever, headache, rigor, myalgia, and dizziness [1,2]. Two coronavirus outbreaks, Serious Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS) have emerged as epidemics with high mortality in the last two decades. SARS CoV transmission from civet cats to humans occurred in China in 2002, and MERS CoV transmission from dromedary camels to humans occurred in Saudi Arabia in 2012 [3]. A cluster of pneumonia cases with an unknown cause surfaced in Wuhan, Hubei Province, China, on 31 December 2019, and China later reported that the outbreak was linked to a seafood market in Wuhan. China revealed the genomic code of the new coronavirus that caused the outbreak for diagnostic purposes on 12 January 2020 [4]. Among other coronaviruses, the virus causing COVID-19 has the advantage of the presence of a unique polybasic cleavage site leading to its increased pathogenicity [5]. In addition, it consists of some structural and some non-structural proteins [6]. The structural proteins include spike (S protein), envelope (E protein), membrane (M protein), nucleocapsid (N protein). In contrast, the non-structural proteins include: main protease (M pro ), papainlike protease (PL pro ), non-structural protein 13 (nsp13, helicase), non-structural protein 12 (nsp12, RNA-dependent RNA polymerase), N-terminal exoribonuclease and C-terminal guanine-N7 methyl transferase (nsp14), uridylate-specific endoribonuclease (nsp15), 2 -O-methyltransferase (nsp16) and nsp10 [7,8]. Because heterocyclic chemicals have been implicated in a variety of diseases, including viral infections, AIDS, and cancer, there is a significant potential for using these numerous nuclei to combat coronaviruses. Some antiviral drugs containing indolin-2-one, such as Arbidol and Marboran/Methisazone, are tested against COVID-19 infection. These hybrids, Arbidol I and Marboran/Methisazone II inhibit membrane fusion and mRNA and protein synthesis, respectively [9]. The use of variants of already recognized antiviral medicines is a useful option until more accurate treatment methodologies for COVID-19 become available [10]. In addition, the pyridine containing hits as 5-chloropyridin-3-yl-1H-indole-4-carboxylate (III) and chloroquine IV, show specificity against SARS CoV 3CL pro . All screened compounds showed potential electrophilic centers which may be capable of forming a covalent bond with the nucleophilic thiol of Cys 145 at the active site of 3CL pro [11] (Figure 1). China in 2002, and MERS CoV transmission from dromedary camels to humans occurred in Saudi Arabia in 2012 [3]. A cluster of pneumonia cases with an unknown cause surfaced in Wuhan, Hubei Province, China, on 31 December 2019, and China later reported that the outbreak was linked to a seafood market in Wuhan. China revealed the genomic code of the new coronavirus that caused the outbreak for diagnostic purposes on 12 January 2020 [4]. Among other coronaviruses, the virus causing COVID-19 has the advantage of the presence of a unique polybasic cleavage site leading to its increased pathogenicity [5]. In addition, it consists of some structural and some non-structural proteins [6]. The structural proteins include spike (S protein), envelope (E protein), membrane (M protein), nucleocapsid (N protein). In contrast, the non-structural proteins include: main protease (M pro ), papain-like protease (PL pro ), non-structural protein 13 (nsp13, helicase), non-structural protein 12 (nsp12, RNA-dependent RNA polymerase), N-terminal exoribonuclease and C-terminal guanine-N7 methyl transferase (nsp14), uridylate-specific endoribonuclease (nsp15), 2′ -O-methyltransferase (nsp16) and nsp10 [7,8]. Because heterocyclic chemicals have been implicated in a variety of diseases, including viral infections, AIDS, and cancer, there is a significant potential for using these numerous nuclei to combat coronaviruses. Some antiviral drugs containing indolin-2-one, such as Arbidol and Marboran/Methisazone, are tested against COVID-19 infection. These hybrids, Arbidol I and Marboran/Methisazone II inhibit membrane fusion and mRNA and protein synthesis, respectively [9]. The use of variants of already recognized antiviral medicines is a useful option until more accurate treatment methodologies for COVID-19 become available [10]. In addition, the pyridine containing hits as 5-chloropyridin-3-yl-1H-indole-4-carboxylate (III) and chloroquine IV, show specificity against SARS CoV 3CL pro . All screened compounds showed potential electrophilic centers which may be capable of forming a covalent bond with the nucleophilic thiol of Cys 145 at the active site of 3CL pro [11] (Figure 1). Furthermore, the purine and pyrimidine derivatives, such as Acyclovir V and Ganciclovir VI and Lopinavir VII, were tested against COVID-19 infection and were shown to be effective in the treatment of COVID-19 [12] involving more than two hydrogen bonds with M pro , which were further analyzed by SARS-CoV M pro inhibition assay. Additionally, purine derivatives have been shown to have significant antiviral action against a variety of viruses. Nucleoside analogs, based on purine moiety, were designed and analyzed for their activity against SARS-CoV using plaque reduction assay in SARS CoV Frankfurt-1 strain infected Vero E6 cells [13] (Figure 2). Furthermore, the purine and pyrimidine derivatives, such as Acyclovir V and Ganciclovir VI and Lopinavir VII, were tested against COVID-19 infection and were shown to be effective in the treatment of COVID-19 [12] involving more than two hydrogen bonds with M pro , which were further analyzed by SARS-CoV M pro inhibition assay. Additionally, purine derivatives have been shown to have significant antiviral action against a variety of viruses. Nucleoside analogs, based on purine moiety, were designed and analyzed for their activity against SARS-CoV using plaque reduction assay in SARS CoV Frankfurt-1 strain infected Vero E6 cells [13] (Figure 2). China in 2002, and MERS CoV transmission from dromedary camels to humans occurred in Saudi Arabia in 2012 [3]. A cluster of pneumonia cases with an unknown cause surfaced in Wuhan, Hubei Province, China, on 31 December 2019, and China later reported that the outbreak was linked to a seafood market in Wuhan. China revealed the genomic code of the new coronavirus that caused the outbreak for diagnostic purposes on 12 January 2020 [4]. Among other coronaviruses, the virus causing COVID-19 has the advantage of the presence of a unique polybasic cleavage site leading to its increased pathogenicity [5]. In addition, it consists of some structural and some non-structural proteins [6]. The structural proteins include spike (S protein), envelope (E protein), membrane (M protein), nucleocapsid (N protein). In contrast, the non-structural proteins include: main protease (M pro ), papain-like protease (PL pro ), non-structural protein 13 (nsp13, helicase), non-structural protein 12 (nsp12, RNA-dependent RNA polymerase), N-terminal exoribonuclease and C-terminal guanine-N7 methyl transferase (nsp14), uridylate-specific endoribonuclease (nsp15), 2′ -O-methyltransferase (nsp16) and nsp10 [7,8]. Because heterocyclic chemicals have been implicated in a variety of diseases, including viral infections, AIDS, and cancer, there is a significant potential for using these numerous nuclei to combat coronaviruses. Some antiviral drugs containing indolin-2-one, such as Arbidol and Marboran/Methisazone, are tested against COVID-19 infection. These hybrids, Arbidol I and Marboran/Methisazone II inhibit membrane fusion and mRNA and protein synthesis, respectively [9]. The use of variants of already recognized antiviral medicines is a useful option until more accurate treatment methodologies for COVID-19 become available [10]. In addition, the pyridine containing hits as 5-chloropyridin-3-yl-1H-indole-4-carboxylate (III) and chloroquine IV, show specificity against SARS CoV 3CL pro . All screened compounds showed potential electrophilic centers which may be capable of forming a covalent bond with the nucleophilic thiol of Cys 145 at the active site of 3CL pro [11] (Figure 1). Furthermore, the purine and pyrimidine derivatives, such as Acyclovir V and Ganciclovir VI and Lopinavir VII, were tested against COVID-19 infection and were shown to be effective in the treatment of COVID-19 [12] involving more than two hydrogen bonds with M pro , which were further analyzed by SARS-CoV M pro inhibition assay. Additionally, purine derivatives have been shown to have significant antiviral action against a variety of viruses. Nucleoside analogs, based on purine moiety, were designed and analyzed for their activity against SARS-CoV using plaque reduction assay in SARS CoV Frankfurt-1 strain infected Vero E6 cells [13] (Figure 2). Molecular hybridization/conjugation is a promising strategy involving combining two pharmacophoric groups via a covalent bond to create a single entity [14]. It was hypothesized that the hybridization of two or more pharmacologically active components in the molecular architecture of hybrid compound/conjugate would prove to be a promising chemotherapeutic agent [15][16][17][18]. As a result of these findings, and as part of our medicinal program aimed at the discovery of novel biologically important heterocyclic compounds with various biological activities [19,20], we integrated the structural features of isatin to design and synthesize a new class of isatin-pyridine-pyrimidine conjugates, hoping to identify novel functional molecules with potent antiviral effects. The results may provide useful information for the design of novel chemotherapeutic drugs ( Figure 3). Molecular hybridization/conjugation is a promising strategy involving combining two pharmacophoric groups via a covalent bond to create a single entity [14]. It was hypothesized that the hybridization of two or more pharmacologically active components in the molecular architecture of hybrid compound/conjugate would prove to be a promising chemotherapeutic agent [15][16][17][18]. As a result of these findings, and as part of our medicinal program aimed at the discovery of novel biologically important heterocyclic compounds with various biological activities [19,20], we integrated the structural features of isatin to design and synthesize a new class of isatin-pyridine-pyrimidine conjugates, hoping to identify novel functional molecules with potent antiviral effects. The results may provide useful information for the design of novel chemotherapeutic drugs ( Figure 3).

Chemistry
The synthetic strategy to design the target spiro-oxindoles based on uracil derivatives are outlined from Schemes 1-3. Both 5-(substituted-1-ylsulfonyl)indoline-2,3-dione 1a-c and 6-aminouracil 2a-e were chosen as the scaffold for annulations of the target congeners in our quest to create new spiro-oxindoles based on uracil derivatives with potential pharmacological significance via one-pot reaction.

Chemistry
The synthetic strategy to design the target spiro-oxindoles based on uracil derivatives are outlined from Schemes 1-3. Both 5-(substituted-1-ylsulfonyl)indoline-2,3-dione 1a-c and 6-aminouracil 2a-e were chosen as the scaffold for annulations of the target congeners in our quest to create new spiro-oxindoles based on uracil derivatives with potential pharmacological significance via one-pot reaction. The substituted 6-aminouracil derivatives 2a-e were obtained by adding urea, methyl urea, and/or methyl thiourea to ethyl cyanoacetate in absolute ethanol containing sodium (sodium ethoxide) by heating under reflux, according to the reported methods [21,22]. Additionally, 5-(substituted-sulfonyl)indoline-2,3-dione 1a-c were prepared according to the previously reported method [23,24].

Scheme 2.
Illustration of the proposed mechanism for synthesis of spiro-oxindole derivatives.
Additionally, the 1 H NMR spectra of 1′,9′-bis(4-chlorobenzyl)-1′H-spiro[indoline-3,5′-pyrido[2,3-d:6,5-d′]dipyrimidine 4d revealed four exchangeable signals at δ 9.84, 11.15, 11.38, 11.90 ppm attributed to four NH groups. Moreover, signals at δ 5.18, 5.28 ppm owned to four protons of two methylene groups (2CH2) and signals at δ 1.32, 1.45, 2.81 ppm related to piperidinyl protons (5CH2). In addition, the eleven aromatic protons ranged between δ 7.06-7.59 ppm. Furthermore, the 13   Additionally, the 1 H NMR spectra of 1 ,9 -bis(4-chlorobenzyl)-1 H-spiro[indoline-3,5pyrido[2,3-d:6,5-d ]dipyrimidine 4d revealed four exchangeable signals at δ 9.84, 11.15, 11.38, 11.90 ppm attributed to four NH groups. Moreover, signals at δ 5.18, 5.28 ppm owned to four protons of two methylene groups (2CH 2 ) and signals at δ 1.32, 1.45, 2.81 ppm related to piperidinyl protons (5CH 2 ). In addition, the eleven aromatic protons ranged between δ 7.06-7.59 ppm. Furthermore, the 13  Eventually, the development of an effective antiviral for SARS-CoV-2, if given to patients early in infection, could help to limit the viral load, prevent severe disease progression, and limit person-person transmission. Benchmark testing of those compounds against other potential antivirals for SARS-CoV-2 with alternative mechanisms of action would therefore be important as soon as practicable [25]. The synthetic compounds 3a-d, 4a-e, and 5a-e were tested in vitro against the SARS-CoV-2 strain isolated from Egyptian patients. The reference drug was chloroquine, a well-known chemotherapeutic agent. The results were presented as IC 50 values and described in the table below (Table 1 and Figure 4). The effect of different concentrations of the compounds on the cellular proliferation of the Vero E6 cell line following 24 h of treatment was determined using MTT assay [26]. According to the results, most of the compounds tested showed moderate to excellent cytotoxic

Plaque Reduction Assay (% of Inhibition SARS-CoV2)
As represented in Table  2 and Figure  5, four spiro[indoline-3,5′-pyrido[2,3-d:6,5-d′]dipyrimidine derivatives 3a, 4b, 4d, and 4e showed high activities against the SARS-CoV-2 that was isolated from Egyptian patients. Firstly, 5-(pyrrolidin-1-ylsulfonyl)-1′H-spiroindoline derivative 3a, which has R= CH3 and X= O, found the most active derivatives among these series 3a-d against replication of the virus with the percentage of inhibition = 84%. Additionally, replacing the methyl group with a benzyl or 4-Cl-benzyl moiety, as in compounds 3b and 3c, the activity decreased to 75% and 0%, respectively. This decrease in the activity because of the presence of benzylidene moiety in general, and the substitution of one hydrogen bond by chlorine atom at position four in the benzylidene moiety, leads to removing the activity.

Plaque Reduction Assay (% of Inhibition SARS-CoV2)
As represented in Table 2 and Figure 5, four spiro[indoline-3,5 -pyrido[2,3-d:6,5d ]dipyrimidine derivatives 3a, 4b, 4d, and 4e showed high activities against the SARS-CoV-2 that was isolated from Egyptian patients. Firstly, 5-(pyrrolidin-1-ylsulfonyl)-1 Hspiroindoline derivative 3a, which has R = CH 3 and X = O, found the most active derivatives among these series 3a-d against replication of the virus with the percentage of inhibition = 84%. Additionally, replacing the methyl group with a benzyl or 4-Cl-benzyl moiety, as in compounds 3b and 3c, the activity decreased to 75% and 0%, respectively. This decrease in the activity because of the presence of benzylidene moiety in general, and the substitution of one hydrogen bond by chlorine atom at position four in the benzylidene moiety, leads to removing the activity.  Three independent experiments were performed for each concentration.

In Vitro Enzymes Assay
The RNA-dependent RNA-polymerase (RdRp) and spike glycoprotein inhibition assay were investigated for the most potent anti-SARS-CoV-2 active hybrids 3a, 4b, 4d, and 4e using chloroquine as a reference drug. The assessment results were summarized in Table 3 as IC50 values in nM. All tested compounds exhibited potent inhibitory activity towards RdRp and spike glycoprotein ranging from 40.23 ± 0.09 to 44.90 ± 0.08 nM and 40.27 ± 0.17 to 44.83 ± 0.16 nM, respectively, when compared with chloroquine as a reference standard, which showed 45 ± 0.02 and 45 ± 0.06 nM against RdRp and spike glycoprotein, respectively. On the other hand, substituting the carbonyl by thiocarbonyl (X = S), as well as presence of the methyl group as described in 5-(pyrrolidin-1-ylsulfonyl)-2 ,8 -dithioxospiroindoline derivative 3d, displayed the percentage of inhibition as 55%. Furthermore, 5-(piperidin-1-ylsulfonyl)-1 H-spiroindoline derivatives 4b and 4e with a methyl moiety and carbonyl or thiocarbonyl (X = O or S) showed the highest activity among the tested derivatives against replication of SARS-CoV-2 with a percentage of inhibition reading 99% and 91%, respectively. Moreover, 1 H-spiroindoline derivative 4c demonstrated an inhibition percentage = 74%, and this decrease in activity may be related to the presence of benzyl moiety. Additionally, compound 4d exhibited good activity with the percentage of inhibition 80%, when the substitutions were R = 4-chloro benzyl and X = O.
Finally, 5-(morpholinosulfonyl)-1 H-spiroindoline derivatives 5a-e revealed weak antiviral activity against SARS-CoV-2, with a percentage of inhibition ranging from 20-70%. Furthermore, in pharmacological terms, the presence of the pyrrolidin-1-yl and piperidin-1yl moieties on the skeleton of 1 H-spiroindoline derivatives 3a-d and 4a-e may be directly responsible for their antiviral activities. In contrast, the presence of the morpholino moiety as in the series 5a-e exhibited weak antiviral agents.

In Vitro Enzymes Assay
The RNA-dependent RNA-polymerase (RdRp) and spike glycoprotein inhibition assay were investigated for the most potent anti-SARS-CoV-2 active hybrids 3a, 4b, 4d, and 4e using chloroquine as a reference drug. The assessment results were summarized in Table 3 as IC 50 values in nM. All tested compounds exhibited potent inhibitory activity towards RdRp and spike glycoprotein ranging from 40.23 ± 0.09 to 44.90 ± 0.08 nM and 40.27 ± 0.17 to 44.83 ± 0.16 nM, respectively, when compared with chloroquine as a reference standard, which showed 45 ± 0.02 and 45 ± 0.06 nM against RdRp and spike glycoprotein, respectively. It was found that 5-(pyrrolidin-1-ylsulfonyl)-1 H-spiroindoline derivative 3a strongly inhibited RdRp and spike glycoprotein (IC 50 = 40.23 ± 0.09 and 40.27 ± 0.17 nM, respectively) and had relatively higher potency than the standard (IC 50 = 45 ± 0.02 and 45 ± 0.06 nM). This increase of the reactivity of 1 ,9 -dimethyl-5-(pyrrolidin-1-ylsulfonyl)- Finally, it can be concluded that the existence of pyrrolidinyl/piperidinyl moiety, in addition to the N-substituted with methyl or 4-chlorobenzyl and X = O as 1 H-spiroindoline derivatives 4b and 4d is essential for antiviral activity (Table 3).

Computational Study of the Binding Mode
Recently, in regards to structure-based drug design, molecular docking has become an evermore essential medication discovery method and is considered the most common approach in pharmaceutical research [27,28]. Additionally, the docking study illustrates the interaction between newly synthesized compounds and active sites in selected proteins at the atomic level to determine the behavior of new promising compounds in the binding site of the target proteins [29]. Furthermore, docking simulation studies have been progressed to visualize, calculate, formulate, and hypothesize about the energy and orientation of new ligands in the active site in the pocket of a target protein [30,31]. Our study involved docking of the most active derivative inside the RNA polymerase (RdRp) (PDB: 6m71) and spike glycoprotein (SGp) (PDB: 6VXX), and the results displayed good binding energy with lower binding affinity values (S) (kcal/mol) with different types of interactions and the docking results represented in Table 4 and Figures 6-9. coprotein more effectively. The docking results were confirmed by the IC50 and inhibition percentage values in experimental activity and it can be said that these derivatives could resist the SARS-CoV-2.

Docking Study Inside RNA-Dependent RNA Polymerase (RdRp) (PDB: 6m71)
Docking of the most active derivatives 3a, 4b, 4d, and 4e, inside the active site of RNA polymerase (RdRp) (PDB: 6m71), exhibited lower binding energy ranging from −18.48 kcal/mol to −15.38 kcal/mol compared with chloroquine as positive control S = −14.94 kcal/mol and remdesivir S = −16.09 Kcal/mol. All the designed derivatives displayed lower binding energy than chloroquine with different types of interaction (H-bonds donor or H-bond acceptors and arene-cation interaction besides hydrophobic interactions). Furthermore, the four derivatives displayed hydrophobic interaction through an aliphatic group (pyrrolidinyl, piperidinyl, and morphinyl), oxygen of sulfonyl, phenyl of indolin-2-one, alongside carbonyl and methyl, and NH of pyridino-pyrimidine moiety with the adjacent residues.
Additionally, compound 3a revealed two hydrogen bonds donor between the amino acid residues, Asp623 and Lys621, with the carbonyl of indolin-2-one and the carbonyl of the pyridino-pyrimidine derivative with bond lengths of 2.90 and 2.70 Å, respectively. In addition, there were two hydrogen bond acceptors between the oxygen of sulfonyl derivatives with Arg553 and Arg555 with bond lengths of 2.83 Å and 2.74 Å (Figure 6), respectively. Furthermore, compound 4b displayed two hydrogen bond acceptors between Lys621 and carbonyl of indolin-2-one derivative, Arg553 and oxygen of sulfonyl, and one hydrogen bond donor between Asp623 with NH of pyridino-pyrimidine moiety with bond lengths of 2.88 Å (30%), 2.73 Å (32%) and 2.56 Å (70%), respectively. Moreover, compound 4d that excreted the highest binding energy S = −18.48 Kcal/mol, showed only one hydrogen bond side-chain acceptor between Lys551 and carbonyl of a pyridino-pyrimidine nucleus with 2.63 Å bond length and 44% strength. In addition, compound 4d can form two arenecation interactions between the two-aryl group of benzylidene derivatives with Lys708 and Arg553, alongside hydrophobic interactions (Figure 7).  On the other hand, compound 4e exhibited the lowest binding energy of the synthesized derivatives. It also exhibited less binding energy than chloroquine (with the binding energy S = −15.38 Kcal/mol) with three hydrogen bonds between Asp760 and Asp623 with NH of a pyridino-pyrimidine nucleus, and Arg553 with the oxygen of sulfonyl derivatives with bond lengths and strengths of 2.38 Å (48%), 2.81 Å (16%), and 2.53 Å (71%), respectively. (All docking figures are represented in the Supplementary Material File).
Finally, chloroquine established binding energy-14.94 Kcal/mol through a one sidechain hydrogen bond donor between the Asp760 and NH of quaternary salt in the side chain with a bond length of 2.42 Å. In addition, the arene-cation interaction between Arg553 and phenyl ring of quinoline ring occurred. Additionally, the remdesivir demonstrated binding energy S = −16.09 Kcal/mol, through three hydrogen bond side-chain acceptors with the residues Arg553, Arg553, alongside one hydrogen bond side-chain donor with Asp623 and one hydrogen bond backbone donor with Thr556 (See Table 4 and the Supplementary Material File for more details).

Docking Study Inside Spike Glycoprotein (SGp) (PDB: 6VXX)
The simulation study was extended to study the reactivity and interaction of the newly designed derivatives inside the active site of spike glycoprotein as a second target for these derivatives. Finally, it can be concluded that the combination of the indolin-2-one sulfonamide derivatives with 6-aminouracil derivatives to form a new hybrid target as 5-(substitutedsulfo nyl)-1 H-spiro[indoline-3,5 -pyrido[2,3-d:6,5-d ]dipyrimidine derivatives will probably hinder the binding of SARS-CoV-2 RNA polymerase and spike glycoprotein more effectively. The docking results were confirmed by the IC 50 and inhibition percentage values in experimental activity and it can be said that these derivatives could resist the SARS-CoV-2.

Geometrical Optimization and Molecular Parameters
The molecular modelling calculation using the DFT method was carried out to evaluate and specify the relationship between the structure of the most active derivatives and the experimental activity results. Additionally, the geometrical optimization of the most promising derivatives 3a, 4b, 4d, 4e and positive controls (Chloroquine and Remdesivir) were calculated by the DFT method that was used to determine the frontier energies of molecular orbital, and therefore deduce some quantum molecular pa rameters according to previously reported methods [32,33], as described in Tables 5 and 6. IP =ionization potential, EA= Electron affinity, X = Electronegativity, ï ï ï = Chemical hardness, S = Chemical softness, µ = Chemical potential, ω = Electrophilic index, CQ = Chloroquine, RDV = Remdesivir.
The results showed the highest occupied molecular orbitals E HOMO ranged between (−6.17 to −6.00 eV) and the lowest unoccupied molecular orbitals E LUMO between (−2.13 to −1.74 eV) compared with chloroquine (E HOMO = −5.66 eV and E LUMO = −1.17 eV) and remdesivir (E HOMO = −6.11 eV and E LUMO = −1.27 eV). The HOMO energy represented the ability of these derivatives to give electrons as electron donors and localized mainly on 2-oxoindoline derivatives. On the other hand, the LUMO energy displayed by a site that has the ability to act as an electron attractive, i.e., electron acceptors due to vacant orbitals and localized at pyrido[2,3-d:6,5-d ]dipyrimidine derivatives moiety.              (∆E) = 4.83 eV, i.e., the synthesized derivatives displayed the most stable conformer as well as the most polarizable. Generally, the molecules with lower energy gaps required less excitation energy and could offer electrons to neighboring biological receptors and were predicted to have high biological potency as confirmed by the experimental results [34]. Moreover, the values of softness S = 0.469-0.495 eV −1 and hardness ï = 2.01-2.13 eV displayed superior activity and lower hardness for the most active derivatives 3a, 4b, 4d, and 4e in comparison to chloroquine (S = 0.446 eV −1 and ï = 2.24 eV) and remdesivir (S = 0.41 eV −1 and ï = 2.41 eV). Additionally, all the tested derivatives exhibited higher electrophilic index (ω) = 3.52-4.27 eV rather than chloroquine (ω) = 2.59 eV and remdesivir (ω) = 2.81 eV.

Molecular Electrostatic Potential (MEP)
Designing new drugs via molecular electrostatic potential (MEP) is an important property in evaluating the shape, size, and charge distribution around the molecules. The distribution of electron density on the surface of molecules provides us information with the regions that have the ability to donate electrons (which act as nucleophiles) and accept electrons (which act as electrophiles) by appearing in different colors. The electrostatic potential increase is in order of blue > green > yellow > orange > red. The red region indicated an electron-rich region (negative charge sites), while the blue region designated an electron deficiency region (partially positive charged), and the green and yellow regions indicated the neutral sites [35,36].
The MEP map of the most promising derivatives 3a, 4b, 4d, 4e and positive controls (Chloroquine and Remdesivir) are presented in Figure 4. The electron-rich areas are represented at the oxygen of indolin-2-one derivatives, oxygen of sulfonyl group, and oxygen of carbonyl of pyrido[2,3-d:6,5-d ]dipyrimidine derivatives moiety, as well as nitrogen of quinoline in chloroquine, while an electron deficiency region is located at the nitrogen of pyrido[2,3-d:6,5-d ]dipyrimidine in spiro derivatives. Moreover, the yellow and green colors appear on carbon and hydrogen of all designed and chloroquine derivatives that characterized neutral sites as represented in Figure 10. The MEP of the most promising derivatives 3a, 4b, 4d, and 4e are abundant with positive and negative regions that are important in the interaction with biological targets, and these results supported a molecular docking study where the hydrogen bonds formed (donors or acceptors) between these regions are different on the active site than in the pocket, as represented in Table 4.
atives that characterized neutral sites as represented in Figure 10. The MEP of the most promising derivatives 3a, 4b, 4d, and 4e are abundant with positive and negative regions that are important in the interaction with biological targets, and these results supported a molecular docking study where the hydrogen bonds formed (donors or acceptors) between these regions are different on the active site than in the pocket, as represented in Table 4.

Chemistry
All melting points were measured with the Electrothermal LA9000 sequence, and the Digital Melting Point Apparatus was left uncorrected. At the Pharmaceutical Analytical Unit, Faculty of Pharmacy, Al-Azhar University, the IR Spectra were calculated using the KBr disc technique on a Nikolet IR 200 FT IR Spectrophotometer, and the values are expressed in (cm −1 ). The 1 H NMR and 13 C NMR Spectra were recorded on a Bruker 400 MHz Spectrometer and the 13 C-NMR spectra was ran at 125 MHz in dimethyl sulfoxide (DMSO-d 6 ) at Applied Nucleic Acid Research Center, Zagazig University, Egypt. Chemical changes were calculated in ppm relative to TMS as an internal norm using DMSO-d 6 as a solvent. At the Regional Center for Mycology and Biotechnology (RCMB) at Al-Azhar University, the mass spectrum was recorded at 70 ev on the DI-50 unit of the Schimadzu GC/ MS-QP5050A Spectrometer and represented as m/z (relative abundance percent). Moreover, the elemental analysis (C, H, N) was performed at Al-Azhar University's Regional Center for Mycology and Biotechnology, and the results were found to be within 0.4 percent of theoretical values, unless otherwise mentioned. The biological activities were performed at the Ministry of Defense's Chemical Warfare Department's Main Chemical Warfare Laboratories. TLC sheets precoated with UV fluorescent silica gel Merck 60 F254 plates were used to track the progress of the reaction, which was visualized using the UV lamp. A mixture of isatin derivatives 1a-c (0.01 mol) and 6-aminouracil derivatives 2a-e (0.02 mol) in acetic acid (7 mL) were heated under reflux for 6-8 hs. The product therefore formed was collected, filtrated and recrystallized from the proper solvent.   -1-ylsulfonyl)-1 H-spiro[indoline-3,5 -pyrido[2,3-d:6,5-d ]dipyrimidine]-2,  2 ,4 ,-6 ,8 (3 H,7 H,9 H,10 H)

Plaque Reduction Assay
Assay was carried out according to the method of [37] in a six-well plate where Vero E6 cells (10 5 cells/mL) were cultivated for 24 h at 37 • C. Cells were grouped into negative and positive control groups without treatment and treated groups. Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) virus (isolated from Egyptian patients with the help of Egy-Army with GenBank code (MT776904) and the link of this on-site represented as (https://www.ncbi.nlm.nih.gov/nuccore/MT776904, accessed on 5 February 2022) was diluted to give 10 3 PFU/well and mixed with the safe concentration of the tested compounds and incubated for 1 h at 37 • C before being added to the cells. Growth medium was removed from the cell culture plates and the cells were inoculated (100 µL/well) with the virus with the tested compounds. After 1 h contact time for virus adsorption, 3 mL of DMEM supplemented with 2% agarose was added and the tested compounds were then added onto the cell monolayer. Plates were left to solidify and were incubated at 37 • C until formation of viral plaques (3 to 4 days) appeared. Formalin (10%) was added for two hours before plates were stained with 0.1% crystal violet in distilled water. Control wells were included where the untreated virus was incubated with Vero E6 cells. Finally, the plaques were counted and the percentage reduction in plaques formation (% Reduction) in comparison to control wells was recorded according to the following Equation (1): Reduction (%) = Viral count (untreated) − Viral count(treated) Viral count(untreated) × 100 (1) where the viral count (untreated) is the viral count in wells where the virus was untreated with the compounds. Additionally, the viral count (treated) is the viral count in wells where the virus was treated with the compounds. Additionally, this performance is shown in Table 2, which presents the inhibition percentage with three independent experiments for each concentration. Enzyme assay was prepared according to reported method [38].

Molecular Docking Study
Firstly, the most active derivatives 3a, 4b, 4d, 4e and positive controls (Chloroquine and Remdesivir) were built using ChemBioDraw 2014, before being exported to Molecular Operating Environmental 10.2008 (MOE) [39,40]. Additionally, the structure was prepared by protonation, and the energy was minimized using the MMFF94x forcefield. The crystal structure of different proteins as a structure of RNA-dependent RNA polymerase (RdRp) (PDB: 6M71) and spike glycoprotein (SGp) (PDB: 6VXX) were obtained from the protein data bank [41,42]. The selected protein and docking process structure was prepared according to the default protocol and according to the previously reported method [43] using the dummy atoms to generate the active site. We selected only one chain for SARS-CoV-2 spike glycoprotein (PDB: 6VXX) (Chain A). Additionally, the SARS-CoV-2 RNA-dependent RNA polymerase (PDB: 6M71) chain A was selected for the docking proposed. The docking simulation inside the active site was performed using the trigonal matcher placement method using London DG as a scoring function. Its energy was represented by Kcal/mol and the top-scoring pose was inspected visually.

Computational Study
The geometrical optimization of the most active derivatives 3a, 4b, 4d, 4e and positive controls (Chloroquine and Remdesivir) were performed using DFT methods through the Gaussian 09 set of the program according to previously reported methods [32]. The global descriptors of chemical reactivity are related to frontier molecular orbital and are calculated simply depending on previously reported methods [44,45]. The global descriptors used in this study are: IP =ionization potential, EA = Electron affinity, X = Electronegativity, ï = Chemical hardness, S = Chemical softness, µ = Chemical potential, ω = Electrophilic index.

Conclusions
In summary, a series of newly synthesized compounds of spirooxindole based on uracil derivatives 3a-d, 4a-e, and 5a-e were prepared through a simple method by heating 5-(substituted-1-ylsulfonyl)indoline-2,3-dione 1a-c under reflux condition with 6-aminouracil derivatives 2a-e. Using the MTT assay, the new hybrids were investigated for anti-viral activity against SARS-CoV-2. As a result, the majority of the compounds tested had moderate to strong anti-SARS-CoV-2 action. In addition, all new compounds were tested for percentage of inhibition using the plaque reduction assay, which revealed that compounds 3a, 4b, 4d, and 4e had a high range of action. The mechanism of action on RNA-dependent RNA polymerase (RdRp) and spike glycoprotein was studied further with these four hybrids. The results were quite promising, as the new hits were found to be as equally effective as chloroquine, which was previously used to treat COVID-19. These findings motivate our team to conduct more advanced studies on the novel hybrids in order to learn more about the mechanism of action in future investigations.