Synthesis, In Silico Prediction and In Vitro Evaluation of Antimicrobial Activity, DFT Calculation and Theoretical Investigation of Novel Xanthines and Uracil Containing Imidazolone Derivatives

Novel xanthine and imidazolone derivatives were synthesized based on oxazolone derivatives 2a-c as a key intermediate. The corresponding xanthine 3-5 and imidazolone derivatives 6-13 were obtained via reaction of oxazolone derivative 2a-c with 5,6-diaminouracils 1a-e under various conditions. Xanthine compounds 3-5 were obtained by cyclocondensation of 5,6-diaminouracils 1a-c with different oxazolones in glacial acetic acid. Moreover, 5,6-diaminouracils 1a-e were reacted with oxazolones 2a-c in presence of drops of acetic acid under fused condition yielding the imidazolone derivatives 6-13. Furthermore, Schiff base of compounds 14-16 were obtained by condensing 5,6-diaminouracils 1a,b,e with 4-dimethylaminobenzaldehyde in acetic acid. The structural identity of the resulting compounds was resolved by IR, 1H-, 13C-NMR and Mass spectral analyses. The novel synthesized compounds were screened for their antifungal and antibacterial activities. Compounds 3, 6, 13 and 16 displayed the highest activity against Escherichia coli as revealed from the IC50 values (1.8–1.9 µg/mL). The compound 16 displayed a significant antifungal activity against Candia albicans (0.82 µg/mL), Aspergillus flavus (1.2 µg/mL) comparing to authentic antibiotics. From the TEM microgram, the compounds 3, 12, 13 and 16 exhibited a strong deformation to the cellular entities, by interfering with the cell membrane components, causing cytosol leakage, cellular shrinkage and irregularity to the cell shape. In addition, docking study for the most promising antimicrobial tested compounds depicted high binding affinity against acyl carrier protein domain from a fungal type I polyketide synthase (ACP), and Baumannii penicillin- binding protein (PBP). Moreover, compound 12 showed high drug- likeness, and excellent pharmacokinetics, which needs to be in focus for further antimicrobial drug development. The most promising antimicrobial compounds underwent theoretical investigation using DFT calculation.


Introduction
Increasing the incidence of microbial resistance to antibiotics due to the overuse and misuse of antimicrobial drugs is the major challenge for controlling the pathogenic bacteria. Recently, the World Health Organization reported that excessive use of antibiotics during the COVID-19 pandemic, as well as previous H 1 N 1 influenza pandemics, increased the incidence of antimicrobial resistance [1][2][3]. Also, the infections by antibiotic -resistant bacteria has been exceeded that of tuberculosis, HIV infections and influenza combined [4,5]. Scheme 1. Chemical structures of Fluconazole, Bacmethrin, and Prismsol, as authentic drugs, and the current synthesized compounds "imidazolopyrimidine-2,4-diones, 5-(imidazole-1-yl) uracils and Schiff base of uracils.

Chemistry
Here, we reported the synthesis of novel purine derivatives and the assessment of their antibacterial and antifungal activities. 5,6-Diaminouracils 1a-e, the key starting material for the obtained compounds, were synthesized by the reaction between ethyl cyanoacetate and different urea derivatives followed by nitrosation and reduction by conventional methods [25,38,39]. The Erlenmeyer azlactones (2a-c), also known as 4-arylideno-5(4H)-1,3-oxazolones, were synthesized according to literature [40][41][42][43] through an aldol-type condensation of a glycine-based azlactone (A) with different aldehydes (Scheme 2).  Moreover, N-Substituted-8-purinyl-2-substituted vinylacetamide 3-5 were prepared by heating 5,6-diaminouracils 1a-c in water bath with azlactone 2c in acetic acid for 1 h as illustrated in Scheme 4. The xanthine derivatives were crystallized from DMF/EtOH. Different analyses were applied to identify the resulting compounds. The most interesting observation in the 1 H-NMR spectra of compounds 3-5 showed the disappearance of both signals of NH 2 (5) and NH 2 (6) around δ 5.0-6.0 ppm in compounds 1a-c and the appearance of signals of NH (7) around δ 12.74-12.66 ppm, NH signals of acetamide group around δ 9.48-9.46 ppm, a singlet signal of CH benzylidene around δ 6.98-6.56 ppm and the upfield singlet signals of methyl group of acetamide around δ 1.91-1.90 ppm. While the IR spectra of the compounds 3-5 displayed a stretching band of NH within the range υ 3294-3278 cm −1 and the C=O groups (Amide I), red-shifted within the range υ 1705-1620 cm −1 . 13 C-NMR spectra indicates the formation of compounds 3-5 through the appearance of the CH benzylidene around δ 112 ppm and the upfield signals of methyl group of acetamide around δ 23 ppm.
Analogously, 5,6-diaminouracils 1a-e were heated under fusion with different 1,3-oxazolones in acetic acid affording 6-amino-1-substituted-5-arylimidazolyluracil 6-13 in a high yield (Scheme 4). Compounds 6-13 are characterized by the disappearance of the NH 2 (5) signals around δ 5.0-6.0 ppm, the presence of the NH 2 (6) signals around δ 7.41-7.23 ppm, in addition to the appearance of singlet signals of CH benzylidene around δ 7.15-6.95 ppm and the upfield singlet signals of methyl group at δ 2.22-2.09 ppm. The IR spectra of compounds 6-13 showed the stretching band of NH within the range υ 3347-3172 cm −1 and the C=O groups (Amide I) within the range υ 1735-1627 cm −1 , The strong asymmetric and symmetric NO 2 stretching bands of the nitro group were noted at υ 1519-1512 and 1388-1381 cm −1 . 13 C-NMR spectra prove the formation of compounds 6-13 through the appearance of the CH benzylidene around δ 124. .08 ppm and the upfield signals of methyl group at the range δ 16.41-15.04 ppm. Plausible mechanisms for these reactions are shown in Scheme 5 which takes place by two ways A or B. The acidic medium A first activates oxazolone carbonyl group via protonation which indicated by the blue arrows. Then, oxazolone ring is ruptured due to uracil amino group (5) attack. After that, pyrimidoimidazolone derivatives are obtained in the procedure of intramolecular cyclization B. Otherwise the amino group (6) attack the carbon of carbonyl group of oxazolone via protonation, then imidazolopyrimidine-2,4diones were formed. Scheme 5. The suggested mechanism for the synthesis of xanthine and imidazolone derivatives. Treatment of compounds 1a,b,e with 4-dimethylaminobenzaldehyde in acetic acid was heated under fusion for 5-7 min to yield the target Schiff base derivatives 14-16 in high yields (Scheme 6). The compound 14 was synthesized by heating 5,6-diamino-1-methyl-2thiouracil (1b) with 4-dimethylaminobenzaldehyde in dioxan at 80 • C with stirring for 2 h, as described by Nagamatsu [44]. Our improved method for synthesis of compound 14 has been conducted at 5-7 min, with 96% yield and affordable lower cost. The structures of the synthesized Schiff base derivatives were confirmed on the basis of their IR, 1 H-NMR, 13 C-NMR, mass spectral data and elemental analysis.

Biological Activity 2.2.1. Antimicrobial Activity
The antimicrobial efficacy of the synthesized uracil derivatives was evaluated by disc diffusion assay for the antibiotic resistant bacteria; E. coli, Pseudomonas sp., Salmonella sp., Staphylococcus aureus and Bacillus sp., in addition to pathogenic fungi (A. flavus and C. albicans). After desired incubation, the biological activity of the tested compounds was evaluated from the diameter of the inhibition zones. From the diameter of inhibition zones (Table 1), compounds 3 and 6 displayed the highest antifungal activity towards C. albicans (10 mm), A. flavus (2 mm), as well as antibacterial activity towards Gram positive bacteria S. aureus (1.0 mm) and Gram negative bacteria E. coli (4 mm). Practically, the compounds 3 and 6 exhibited the same physiological effect on the growth of bacteria and fungi revealing their chemical proximity. The compound 16 displayed the highest antifungal activity; A. flavus (10 mm) and C. albicans (12 mm), as well as highest activity against E. coli (7 mm), Pseudomonas sp. (1.9 mm) and Bacillus sp. (1.2 mm). The antibacterial and antifungal activities of the tested compounds were normalized to the authentic commercial antifungal (Fluconazole) and antibacterial (Bacimethrin, Primsol) drugs. As well as, the compounds # 11, 12 and 13 exhibited a strong antifungal activity against C. albicans, A. flavus and E. coli. However, the compounds 4, 5, 7, 8, 10, displayed a mild activity towards the tested microorganisms. Obviously, a fluctuation has been observed on the activity of the synthesized compounds towards the tested microorganisms, suggesting the structure-activity relationships of the synthesized compounds for each microorganism. DMSO at the 10% has been used as negative control, displaying no signs of antimicrobial activity towards the tested microorganisms.
For further authentication of the antimicrobial activity of the tested compounds, different concentration of each compound was used, and after the desired incubation the kinetics of microbial growth of each microorganism was determined ( Figure 1). The IC 50 values for each compound against the tested microorganisms were summarized in Table 2    The cytomorphological deformations of the microorganisms in response to the tested compounds were evaluated by TEM. The effects of the most potent uracil derivativecompounds on the cellular structure of Gram positive and Gram-negative bacteria were assessed. The influence of the most potent compounds on the cellular structures of E. coli was assessed ( Figure 2). An obvious cellular deformation on the cell wall and cell membrane was observed upon treatment with the compounds 3, 12, 13, 14 and 16. With treatment of E. coli by the compounds at their MIC values, a complete distortion and ruptures to the bacilli type, comparing to regular bacillus shape of control. The cellular structure of Pseudomonas sp., upon treatment with the compounds 3, 12, 13, 14, and 16 was evaluated by TEM. From the results of TEM analysis ( Figure 3), a complete morphological deformation to bacterial cells from short rod cells into rounded cells, was observed. A complete rupture on the cell walls and cell membranes was observed, with obvious leakage to the cellular organelles, caused by the partial ruptures on the cell walls. The rationality of destructive effect of the compound 16 may be due to the suppression of pivotal biological secondary metabolites especially the proper cellular morphogenesis as revealed from the binding affinity with ACP domain of polyketide synthase type I. The cellular deformation of the Klebsiella sp., in response to the compounds 3, 12, 13, 14 and 16 has been observed ( Figure 4), a complete distortion to the cellular identity of Klebsiella in response to treatment with these compounds. From the TEM micrograph, a dramatic deformation has been noticed to the bacterial cell walls and various organelles in response to the tested compounds. A remarkable alternations and deformation to the cells, irregularity on shape, loss of symmetry, equivocal distribution of the cytoplasmic organelles, especially the cell wall and plasma membrane that becoming a highly denser, have been observed. The powerful antibacterial activity of the compounds 12 and 16 elaborates from their higher binding affinity with the ACP domains of PKS type I [45] and DNA gyrase [46], blocking the synthesis of essential secondary metabolites and DNA replication. Increase on the plasma membrane permeability causes a strong cracking to the bacterial cell leading to a visual release of the cytoplasmic materials. With the treatment, almost all cellular organelles were strongly deforming, in contrary to the control cells which are usually short bacilli, with a multilayered cell surface consisting of an outer membrane, a peptidoglycan layer in the periplasmic space, the cytoplasmic membrane was relatively uniform [47].

Cytotoxicity Analysis Cytotoxicity Assay
The cytotoxicity of the experimental compounds was evaluated towards VERO-E6 cells using the (MTT) assay. As revealed from the half maximal cytotoxic concentration (CC 50 ), the tested compounds 3, S-Aza 1, 5, 6, 8, 11, 12, 13, 14, 9 displayed relatively no signs of toxicity in which the CC50 values were ranged from 160 to 676 µg/mL as shown from the Supplementary Material ( Figure S4). Negative control of 10% DMSO has been used. The DMSO (10%) did not display any signs of toxicity to the VERO-E6 cells.

Molecular Docking Study
To take one step further, and to determine the mode of action of the tested compounds as potential anti-microbial agents, molecular-docking study was employed to determine the binding modes against three important proteins implicated in fungal and bacterial targeting mechanism including APC, DNAg and PBP. These targets were selected based on structure activity relationship with the corresponding co-crystalized ligands with our tested compounds. The cocrystal ligand for the three proteins (PNS for APC protein, CBN for DNAg protein, and PNM for PBB protein) were redocked to assure the validity of the docking parameters and methods to represent the position and orientation of the ligand detected in the crystal structure. The difference of RMSD value between cocrystal ligands to the original cocrystal ligand was <2 Å which approved the accuracy of the docking protocols and parameters.
To elucidate the mode of action of the most promising antimicrobial compounds "3, 12, 13 and 16" docking studies were conducted. Compounds 3, 12, 13 showed an affinity to bind with APC protein active site, that was evidenced by low binding energy compared to reference ligand 4'-Phosphopantetheine (PNS) (−5.1, −5.1, −5.2 Kcal/mol, respectively), while compound 16 has the same binding free energy value of the reference ligand (−4.5 Kcal/mol) ( Table 3). All tested compounds formed one hydrogen bond with key amino acid residue Glu65 (Table 3). Docking results against DNAg protein, showed less binding affinity of all tested compounds against this protein, which was evident by high free of energy for all compounds compared to reference ligand Clorobiocin (CBN) ( Table 3). In addition, docking score against PBP protein, showed a favorable binding affinity to this protein, as revealed from the low values of energy of free binding compared to reference ligand penicillin G (PNG) ( Table 3). Compounds 3, and 12 showed the lowest values binging free energybinging free energy (−8.7 and −8.2 Kcal/mol), compared to PNG (−6.7 Kcal/mol), moreover compound 3 formed one hydrogen bond interaction with amino acid (ASN674: 2.006 Å), while compound 12 formed three hydrogen bonds with the key amino acid residues in the PBP pocket (SER470: 2.118; SER434: 3.269; THR672: 2.797) ( Table 3), and ( Figures 5 and 6).

Pharmacokinetics and ADME Activity
To predict the physicochemical and drug-likeness properties, bioinformatics study for the compound 12 was carried out (Table 4). Bioavailability Radar has been generated to check the suitable physicochemical properties for oral bioavailability, based on 6 characters including (lipophilicity: XLOGP3 between −0.7 and +5.0, size: M W between 150 and 500 g/mol, polarity: TPSA between 20 and 130 Å, solubility: log S not higher than 6, saturation: fraction of carbons in the sp 3 hybridization not less than 0.25, and flexibility: no more than 9 rotatable bonds) (Table 4, Figure 7). Compound 12 exhibited good bioavailability, and good water solubility, while it showed good pharmacokinetics with high gastrointestinal absorption (GI), and it is does not permit blood brain barrier (BBB) as it depicted in the boiled Egg model ( Figure 8). In addition, compound 12 has good drug likeness properties according to Lipinski rule of five for drug discovery without any violation.

Structural Parameters and Models
N-(1-(3-benzyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-8-yl)-2-(4-chlorophenyl)vinyl) acetamide, compound 3. The activity of hetero cyclic compounds is mainly determined by its structure, this compound has many characteristic structural features. The bond length N22-C21 is 1.395 Å, C21-N0 is 1.389 Å, there is a single bond characters C and N atoms [48], the bond lengths of N17-C10 and N11-C10 are 1.269 and 1.339 Å, respectively, whereas the N17-C10 bond has double bond characters, while the N11-C10 has single bond characters [48]. Detailed analysis of corresponding bond lengths in various hetero cyclic compounds were given elsewhere [49]. All distances and angles between the atoms of the ligand I was given in (Table 5). The C21-O24 bond length is 1.208 Å, while the C19-O23 bond length is 1.210 Å [50], C13-N12 is 1.384 Å [49] and the C13=O14 bond length of carbonyl group is 1.209 Å [50]. Table 5. Equilibrium geometric parameters bond lengths (Å), bond angles (°), dihedral angles (°), total energy (k cal/mol), heat of formation (k cal/mol), Mulliken charges over atoms and dipole moment of the N-(1-(3-benzyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-8-yl)-2-(4-chlorophenyl) vinyl)acetamide, compound 3 by using DFT calculations.   The molecule is a highly sterically-hindered, there are three aromatic rings. The molecule is non-planar, there are two plans one occupied by the two aromatic rings and other plane occupied by the purine ring and also, the two plans are perpendicular respect to each other. The purine ring is lying out of plane of the plane occupied by the two aromatic rings, the dihedral angles C26C25N22C18 is −84.71 • ≈ 90 • and the acetamide group is lying in perpendicular plane respect to purine ring system, the C10C9N12C13 is 82.58 • ≈ 90 • , where the values are neither zero nor 180 • . The C10N17C18N22, 178.05 • ≈ 180 • , also the dihedral angles: N22C25C26C27, −175.85 • ≈ 180 • and N22C25C26C20, −4.19 • ≈ 0.00 • , where the values are nearly zero or 180 • . Figure 9, the terminal benzyl group is lying in perpendicular plane respect to purine ring system, with a dihedral bond angle C8C9N10N17 is 102.31 • ≈ 90 • , Figure 9 shows the optimized geometrical structure of the compound (3), N-(1-(3-benzyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-8-yl)-2-(4-chlorophenyl)vinyl)acetamide, the dihedral angles O23C19N20C21 is 178.85 • ≈ 180 • and O24C21N20C19 is −179.53 • ≈ 180 • , which confirms that the O23 and O24 are lying in the same plane of the nitrogen atom N20 of purine ring. The oxygen atom O14 of acetamide group is lying in trans position respect to the C9 which bonded with purine ring with dihedral angle 168.98 • ≈ 180 • , while the others two oxygen atoms O23 and O24 are lying in the same plane away from the plane occupied by O14 of acetamide group, this confirms that one oxygen atoms of these three atoms can be chelated as donor with any acceptor according to the value of accumulated charge density lying on it. The value of bond angle C4C8C9 is 131.89 • , C26C25N22 is 113.98 • , C8C9C10 is 126.94 • and C9N12C13 is 124.18 • reflects on sp2 hybridization on the C8, C25, C9 and N12. There is building up of charge density on the hetero atoms, nitrogen and oxygen atoms, the highest charged atoms are O14 of carbonyl group in acetamide with −0.488, while the others two oxygen atoms O23 and O24 have lower charge density values, −0.446 and −0.432, so the active site in this compound is O14 as set in docking study in biological investigation, the compound 3 favors bonding with any receptor through the oxygen atom O14 of carbonyl group. The accumulated charge density on nitrogen atoms N11, N17, N20, N22 and N12 involved in this compound are −0.107, −0.334, −0.432, −0.198 and −0.284, respectively. The value of energy of this compound is −172,059.722 kcal/mol, heat of formation of this compound is −10,607.049 kcal/mol and slightly weak dipole 3.078D. Detailed analysis of corresponding bond lengths in various hetero cyclic compounds were given elsewhere [51,52]. All distances and angles between the atoms of the ligand are given in (Table 5).
The pyrimidine ring is lying out of plane of the plane occupied by the others two aromatic rings, the dihedral angles C17N13C1C3 and C17N13C1C2 are 109.33 • ≈ 90 • and −70.47 • ≈ 90 • , respectively, where the values are neither zero nor 180°. These values confirm that the pyrimidine ring is lying in perpendicular plane respect to the central aromatic imidazole ring of the compound. Also, the value of dihedral angle of C5N6C9C11 is 90.54 • ≈ 90 • , this value confirms that the terminal ethyl group is lying in perpendicular plane respect to the pyrimidine ring. The dihedral angle N16C15C20C21 is 171.89 • ≈ 180 • , also the dihedral angle, C15C20C21C22 is 176.18 • ≈ 18 • , these values confirm that the benzylidene ring is lying in the same plane of the imidazole ring, also the values of the dihedral angles C26C21C20C15 and C21C20C15C14 are 7.63 • ≈ 0.00 • and 4.55 • ≈ 0.00 • as shown in Figure 10. The main observation from the dihedral angles of the optimized geometrical structure of the compound 12 that the two oxygen atoms O7 and O8 are lying in the same plane, while the oxygen atom O18 is lying out of this plane, one of these atoms can be participate in the chelation as donor atom with any acceptor atom. The calculated mullkine charge density over these oxygen atoms helps us to detect which of these oxygen atoms can be behaved as active site in bonding with acceptor atom in docking calculations.
There is building up of charge density on the all hetero atoms nitrogen and oxygen atoms, the highest negative charged atoms are O7 and O8, −0.442 and −0.432, while the negative charges on nitrogen atoms of the pyrimidine ring, N4 and N6 are −0.328 and −0.208, also the negative charges on the nitrogen atoms of imidazole ring, N13 and N16 are −0.135 and −0.385. These values of charges confirmed that the oxygen atoms can be participate in bonding with any acceptor more than nitrogen atoms involved in the compound 12, so the oxygen atom, O7 can be considered as active site in this compound and can bond through it with any acceptor as shown in the docking calculations. The charge density on other atoms less negative than these atoms as listed in Table S1. The energy of the compound 12, 6-amino-5-(4-benzylidene-2-methyl-5-oxo-4,5-dihydro-1H-imidazol-1-yl)-1-ethylpyrimidine-2,4(1H,3H)-dione is−141,856.224 kcal/mol, heat of formation is −8661.756 kcal/mol and slightly high dipole 8.326D. From the value of energy we found that this compound has higher energy value greater than the compound 3, which has lower energy and lower dipole moment. This compound is similar to the compound (12) with replacement of the C=O of the pyrimidine ring with C=S, the bond length C3-N4 is and C5-N4 are 1.369 and 1.383 Å, there is a single bond characters C and N atoms in the pyrimidine ring [48], also the bond lengths of N12-C16, C13-N12 and C14-N15 of the imidazole ring are 1.352, 1.348 and 1.359 Å, these values reflect that these bonds have single bond characters, while the bond distance of C16-N15 of the imidazol ring is 1.339 Å, whereas this bond has double bond characters [48] as shown in the ( Figure 11). All calculated bond distances, bond angles and dihedral angles between the atoms of the compound 13 are given in (Table S2). The C5-S8 bond length is 1.692 Å [53] and the C3-O7 bond length is 1.209 Å [50], while the C13-O17 bond length is 1.208 Å [50].
There is no significant difference between the structure of this compound and compound (12), the benzylidene ring and dihydropyrimidine ring are lying in the same plane, the dihedral angles C13C14C19C20 and C14C19C20C25 are 6.42 • ≈ 0.00 • and 11.35 • ≈ 0.00 • , these values confirm that the benzylidene ring and dihydropyrimidine ring are lying in the same plane, while the dihedral angle C3C1N12C13 is −49.18, so the imidazol ring is lying out of the plane occupied by others two rings, as shown in the ( Figure 11). The dihedral angles C2N6C5S8 is −153.03 • and O7C3C1C2 is 169.59°, which confirms that the O7 and S8 are lying in the same plane of the nitrogen atom N4 of dihydropyrimidin ring, while the dihedral angle O17C13N12C1 of the imidazol ring is 79.9 • ≈ 90.0 • , so the oxygen atom O17 is lying away from the other plane occupied by O7 and S8. In this compound there are many hetero atoms can behave as active donating site, the active donating site can be detected from the value of accumulated charges spreading over hetero atoms.  Table S3, so that the compound 13 have a power full chelating sites oxygen atom, O7 and sulfur atom S8 more than others compounds, from the optimized structure the bonding of compound 13 favors bonding through the oxygen atom O7 not through sulfur atom S8, the dihedral angles C2N6C5S8 and C9N6C5S8 are 153.03 • and −43.49 • , these values reflected that sulfur atom is not completely planer, while the dihedral angles N12C1C3O7 and O7C3C1C2 are −13.93 • ≈ 00.0 • and 169.59 • ≈ 180.0 • , these values reflected that oxygen atom O7 is nearly planer, so oxygen atom, O7 is considered as active site in bonding with acceptor atom. The energy of the compound 13 is −132,805.565 kcal/mol, heat of formation is −8086.449 kcal/mol and highly dipole 8.586 D more than others previous compounds.
6-Amino-5-((4-(dimethylamino)benzylidene)amino)-1-ethylpyrimidine-2,4(1H,3H)dione, compound (16) This compound is more energetic than all studied compound, the energy of it is −123,990.259 kcal/mol, so this compound is less stable than others. The heat of formation of this compound is −8026.826 kcal/mol, the value of dipole moment is 6.051 D, this value is slightly high, this value is greater than dipole moment of the compound 3 and less than compounds 12 and 13, this result attributed to the presence of C25=N24 bond as bridge between the two aromatic rings, which facilitates the rotation of any aromatic ring around it. All geometrical parameters of the compound 16, 6-amino-5-((4-(dimethylamino) benzylidene)amino)-1-ethylpyrimidine -2,4(1H,3H)-dione, bond lengths, bond angles and dihedral angles are given in (Table S3). The bond length C16-N15 is 1.379 Å is longer slightly than bond length of C14-N15 is 1.376 Å, from the values of bond lengths, C13-N17 and N17-C16, there is a single bond characters C and N atoms [49], the bond length of C25-N24 is 1.347 Å, whereas this bond has a double bond characters [49]. Detailed analysis of corresponding bond lengths in various hetero cyclic compounds were given elsewhere [50]. The C14-O18 and C16-O19 bond lengths of the pyrimidine ring are 1.211 Å and 1.208 Å [50].
As seen in (Figure 12) of the optimized geometrical structure of the compound 16, The compound is non planar and there are two planes each of them ia perpendicular respect to each other, the dihedral angles C25N24C12C14 and C25N24C12C13 are 50.26°and −131.88 • , where the values are neither zero nor 180 • , these values reflected that the plane occupied by pyrimidine ring is perpendicular with the plane occupied by benzylidene ring, the dihedral angle C13N17C20C22 is −83.59 • ≈ 90 • , this value confirms that the terminal ethyl group is lying in plane perpendicular with the plane occupied by pyrimidine ring. Also, the value of dihedral angle of C5N6C9C11 is 90.54 • ≈ 90 • , this value confirms that the terminal ethyl group is lying in perpendicular plane respect to the pyrimidine ring. The dihedral angles C20N17C16O19 and N24C12C14O18 are 1.65 • ≈ 0.00 • and 0.10 • ≈ 0.00 • as shown in (Figure 12). The main observation from the dihedral angles of the optimized geometrical structure of the compound 16 that the two oxygen atoms O18 and O19 are lying in the same plane, one of them can be considered as an active site in this compound. The calculated mullkine charge density over these two oxygen atoms can be used to detect which of them can be used in bonding with acceptor atom in docking calculations. There is building up of charge density on the all hetero atoms nitrogen and oxygen atoms, the highest negative charged atoms are O18 and O19, −0.448 and −0.459, while the negative charges on nitrogen atoms of the pyrimidine ring, N15 and N17 are −0.429 and −0.203, also the negative charge on the terminal nitrogen atom of amino group, N9 is −0.437. These values of charges confirmed that the oxygen atom, O18 can behave as an active site in the compound 16 and can bond through it with any acceptor as shown in the docking calculations. The charge density on other atoms less negative than these atoms as listed in (Table S4).

Charge Distribution Analysis
The charge distribution analysis on the optimized geometry configuration of the all studied compounds was made on the basis of natural population analysis (NPA). The charge distribution on the hetero atoms, nitrogen, oxygen and sulfur of the different studied compounds indicates the presence of a net negative pole attributed to the presence of polar atoms beside the high degree of planarity of these molecules, as a result the molecules are a slightly high dipole in case of compound 3, 8.586D, while compound 2 has slightly weak dipole, 3.078D less than others compounds. There is building up of charge density on the donating oxygen, nitrogen and sulfur atoms involved in the studied compounds. All calculated parameters, dipole moment, mulliken charges and total energies are given in (Table 6). From the factors affecting on the degree of stability of these compounds is the value of the total energy and distribution charges over all atoms in compounds. From the values of energies, the compound 3, N-(1-(3-benzyl-2,6-dioxo-2,3,6,7-tetrahydro-1Hpurin-8-yl)-2-(4-chlorophenyl)vinyl)acetamide is considered as more stable than others, with more negative energy value less than second compound 12 by 30,203.498 kcal/mol, while the third compound in stability is 13, the difference in energy between it and the first compound 3 is 39,254.157 kcal/mol. The energy difference between the less stable compound 16 and the most stable compound 3 is 48,069.463 kcal/mol.

Molecular Orbitals and Frontier
Molecular orbitals play also an important role in the electric properties, as well as in UV-Vis [54]. An electronic system with smaller values of HOMO-LUMO gap should be more reactive than one having a greater energy gap [55]. The energy gap, ∆E of the studied compounds varied between 0.061 for Compound 13 which more reactive and 0.082 eV for compound 3 which less reactive, so electron movement between these orbitals could be easily occur by decreasing the value of energy gap, so that there is a peak around 250 nm in the UV-Vis spectra for all studied compounds as shown in UV-Vis spectra for all studied compounds. The energy gap is closely associated with the reactivity and stability of the executed compounds and shows the nature of the compound with low kinetic stability and slightly high chemical reactivity. On the other hand, the adjacent orbitals are often closely spaced on the frontier region.
The nodal properties of molecular orbitals of studied complexes in Figure 13 are illustrative and suggest orbital delocalization, strong orbital overlap, and low number of nodal planes. The energy difference between HOMO and LUMO (energy gap, ∆E) for the all studied compounds varied according to the composition of each compound or the type of substitutions as shown in (Table 7), and (Figure 13), shows the isodensity surface plots of HOMO and LUMO for the studied compounds. Hard molecules have high the HOMO-LUMO gap, and soft molecules have smaller HOMO-LUMO gap [56]. The values of η and ∆E (HOMO-LUMO) are given in (Table 7). It is obvious that the compound 3 is considered as most hard with η equals to 0.041, and compound 13 is considered as less hard molecule with η equals to 0.031, also the electronic transition within the compounds is easy as indicated from the ∆E. There are some quantum chemical parameters depending upon the energy values of HOMO and LUMO were calculated as global softness (S), electro negativity (χ), absolute softness (σ), chemical potential (Pi), global electrophilicity (ω) and additional electronic charge (∆N max ) of the all studied compounds. From these values, the compound 3 is absolute more soft according to the (σ = 24.390 eV), while the compound 13 is treated as less soft compound (σ = 32.258 eV). As seen in Figure 13, the electron density of HOMO, Ψ 113 of the compound 3 is delocalized over all atoms with highly percent on the two oxygen atoms of the C=O groups of the 2,6-dioxopurine ring with percent 55.5%, with spreading of the remaining percentage on the remaining components, with 20.2% on the acetamide group, 15.4 % on the chlorophenyl ring, 7.9% on the vinyl group and finally 1% on the benzyl group. The LUMO, Ψ 114, the electron density delocalized over all atoms of compound 3, the most percent (53.7%) is localized on the 2,6-dioxopurine ring with spreading of remaining percentage over the remaining components except benzyl group with 0.0%, as given in (Table 8). Also the HOMO, Ψ 89 of the compound 12, the electron density delocalized over all atoms, the most percent (61.2%) is localized on the 2-methyl-5-oxoimidazol ring and 25.6% on the 2,4-dioxo-6-aminopyrimidine ring with small percent 13.2% on the benzylidene ring, without any percent on the terminal ethyl group. The LUMO, Ψ 90 of the compound 12, the electron density delocalized over all atoms, the most percent (51.9%) is localized on the 2,4-dioxo-6-aminopyrimidine ring and 48.1% on the benzylidene ring only. The HOMO of the compound 13, Ψ 89, the electron density localized over 2-thio-4-oxo-6-aminopyrimidine ring with percent (56.3%) and 41.9% on the 2-methyl-5-oxoimidazol ring, with small percent on benzylidene ring (1.8%) and without any percent on the terminal methyl group. The LUMO, Ψ 90, the electron density delocalized over all atoms, the most percent (52.7%) is localized on the 2-methyl-5-oxoimidazol ring and 31.4% on the benzylidene ring with small percent on the terminal methyl group (15.9%), while 0.0% on the 2-thio-4-oxo-6-aminopyrimidine.
In compound 16, the HOMO, Ψ 80, the electron density delocalized over all atoms of this compound with percent (41.3%) on the 2,4-dioxo-6-aminopyrimidine group, 33.6% on the dimethylaminobenzylidene ring and 25.1% on the Amino group with 0.0% on the terminal ethyl group. The LUMO, Ψ 81, the electron density delocalized over all atoms in this compound with 58.7% localized on the dimethylaminobenzylidene ring, 21.8% on the amino group and 19.5% on the 2,4-dioxo-6-aminopyrimidine ring without any percent on the ethyl group 0.0%. The variation in the distribution of electron densities on the components of HOMO and LUMO for each compound play important role in the variation of electronic transitions occurred in compounds and effect on the degree of biological effect of them. Table 8. Computed excitation energies (eV), electronic transition configurations and wave lengths (nm) of the obtained stable compounds; 3 and 12 by using B3LYP/Cep-31G for the, (f ≥ 0.001) f = oscillator strengths compared with experimental data of UV-Visible spectra, λ max , nm. The TD-DFT at the B3LYP level by using G98W program proved to give accurate description of the UV-vis. spectra [57,58]. Time-dependent density functional response theory (TD-DFT) has been recently reformulated [59] to compute discrete transition energies and oscillator strengths and has been applied to a number of different atoms and molecule. Bauernschmitt and Ahlrichs [60] included hybrid functional proposed in the calculation of the excitation energies. These hybrid methods typically constitute a considerable improvement over conventional Hatree-Fock (HF) based methods. In this work, the optimized geometry was calculated and was used in all subsequent calculations; the wave functions of SCF MOs were explicitly analyzed. The calculated wave functions of the different MOs reflect and suggest the fraction of the different fragments of the compounds contributing to the total wave functions of different states. The results indicate that there is an extent of electron delocalization in the different molecular orbitals.
In compound 13, there are five peaks obtained in spectra, while theoretically there are nine transitions can be observed as given in (Table S4). The first peak obtained at 267 nm, this peak is composed of two theoretical transition states overlapped together at 260.77 nm and 263.14 nm, these transitions are assigned to π-π *. The second peak obtained at 290 nm, we found that this peak is consist of two overlapped excited states at 284.13 nm and 284.13 nm, the first one is assigned to π-π * and the second is assigned to n-π * transitions. Also, the third peak observed at 322 nm is composed of three overlapped transition states at 288.71 nm, 326.03 and 328.13 nm, the first and third transition states assigned to π-π * transitions, while the second transition states assigned to n-π* transition. The fourth peak observed experimentally at 339 nm and this peak is agree theoretically with these electronic transitions, H-10 → L (2.8%), H-9 → L (12.7%), H-9 → L + 4 (8.1%), H-6 → L (9.7%), H-6 → L + 5 (8.5%), H-5 → L (4.1%), H-5 → L + 5 (10.7%), H-4 → L (82.6%), H-2 → L + 2 (17.3%), which obtained at 339.43 nm and 3.7876 eV, this transition is assigned to n-π*. The last peak observed experimentally at 357nm, this peak can be represented by these transitions, H-2 → L + 1 (3.2%), H-3 → L (7.8%), this excited state is assigned to n-π * transition, the theoretical wave length corresponding with excited state is 359.12 nm with energy equals to 3.4525 eV.
In (Table S5), the number of excited states involved in compound 16 are five transition states, there are four of them are assigned to π-π *, which can be observed at 249.89, 313.44, 335.42 and 365.71nm, while fourth excited state is assigned to n-π *, this transition can be observed at 278.83 nm.

Materials and Instruments
All melting points were determined by an Electrothermal Mel.-Temp. II apparatus and were uncorrected. Element analyses were performed at Regional Center for Mycology and Biotechnology at Al-Azhar University. The infrared (IR) spectra were recorded using potassium bromide disc technique on Nikolet IR 200 FT IR (Perkin-Elmer, Akron, OH, USA). Mass spectra were recorded on DI-50 unit of Shimadzu GC/MS-QP 5050A at the Regional Center for Mycology and Biotechnology at Al-Azhar University. The proton nuclear magnetic resonance ( 1 H-NMR) spectra were recorded on Bruker 400 MHz Spectrometer and 13 C-NMR spectra were run at 125 MHz in dimethylsulfoxide (DMSO-d6) and TMS as an internal standard, Applied Nucleic Acid Research Center, Zagazig University, Egypt. All new compounds gave corresponding elemental analyses (C, H, N, typically ±0.3%). All reactions were monitored by TLC using precoated plastic sheets silica gel (Merck 60 F 254 ) and spots were visualized by irradiation with UV light (254 nm). The used solvent system was chloroform: methanol (9:1) and ethyl acetate: toluene (1:1).

Antibacterial and Antifungal Activity
The biological activity of the synthesized uracil derivative compounds was assessed against various antibiotic resistant bacteria such as E. coli, Pseudomonas sp., Salmonella sp., S. aureus and Bacillus sp., isolated from different human medical specimens (Al-Ahrar Hospital, Zagazig, Egypt). These bacterial isolates were identified based on their biochemical features according to Bergey's manual [61][62][63]. In addition, the activity of these compounds was evaluated against various pathogenic fungi; A. flavus and C. albicans. The antimicrobial activity of the tested compounds was determined by the disc-diffusion assay [64]. One mL of bacteria cell suspensions (10 7 CFU/mL) from 24 h old cultures, was seeded into the nutrient agar medium, shacked vigorously then poured into sterilized Petri-dish. After medium solidification, 20 µL of the different concentrations of the synthesized compounds were loaded into filter paper disc (8.0 mm) that placed on the plate surface, then, the cultures were incubated for 37 • C for 24 h. The experimental compounds were completely dissolved in 10% DMSO, as well as, the further dilutions were made using 10% DMSO. The solubility and stability of the synthesized compounds were assessed from the UV-Vis spectral analysis. The compounds were completely dissolved and with no signs of precipitation over 5 days at room temperature ( Figure S1). After incubation, observed for the zone of inhibition and diameter of these zones were measured in millimeters. The inhibition area was calculated comparing to the total bacterial growth area using Image J software portal [65].

Transmission Electron Microscopic Analysis
For evaluation of the cytomorphological alterations induced by the most active compounds, the most resistant bacteria were treated with the compounds at their MIC doses, and the cellular organelles were investigated. The IC 50 value was expressed by the concentration of compounds inhibiting the growth of 50% of initial microbial growth, as revealed from the regression curve. The Minimum Inhibitory Concentration (MIC) was expressed by the concentration of compound preventing the visible growth of microorganism. The IC 50 value and MIC value of the potent compounds 3, 12, 13, 14 and 16 were summarized in Table S6. The bacterial broth culture was incubated for 6 h, then amended with the tested compounds at their MIC dose, and continue for 24 h incubation at the same conditions. The bacterial cells were collected by centrifugation at 6000 rpm for 10 min, washed with sterile distilled water. The samples were prepared by immersing in primary fixative (2.5% glutaraldehyde buffered to pH 7.4 with 0.2 M phosphate buffer) for 3 h, and then post-fixed in 1 % osmium tetraoxide for 2 h, buffered with the phosphate buffer for 30 min. All steps of fixation were carried out at 4 • C, the samples were dehydrated in a series of ethanol (50% to 100%), then embedded in resin capsule, that sectioned into ultrathin section of 70 nm, loaded on copper grids and contrasted with uranium acetate and lead citrate prior examination by JEOL-1200 EX microscope (National Research Centre, Giza, Egypt).

Cytotoxicity Assay
The cytotoxicity of the experimental compounds was evaluated towards VERO-E6 cells using the 3-(4,5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT) assay with slight modifications [64]. The stock solutions of compounds were prepared by dissolving in 10% DMSO in ddH 2 O and diluted further to the working solutions with DMEM. Briefly, the cells were seeded in 96 well-plates (100 µL/well at a density of 3 × 10 5 cells/mL), incubated for 24 h at 37 • C in 5%CO 2 , the cells were treated with various concentrations of the tested compounds in triplicates, and further incubated for 24 h. The supernatant was discarded and cell monolayers were washed with sterile phosphate buffered saline and the MTT solution (20 µL of 5 mg/mL stock solution) was added to each well and incubated at 37 • C for 4 h followed by medium aspiration. The developed formazan crystals were dissolved in 200 µL of acidified isopropanol and measured at λ 540 nm using a multi-well plate reader. To assess the half maximal cytotoxic concentration (CC 50 ), the percentage of cytotoxicity compared to the untreated cells was determined with the following equation.
% cytotoxicity = ((absorbance of cells without treatment − absorbance of cells with treatment) × 100)/(absorbance of cells without treatment) The plot of % cytotoxicity versus sampleconcentration was used to calculate the concentration which exhibited 50% cytotoxicity (CC 50 ).

Molecular Docking Study
The structures of all tested compounds were modeled using the Chemsketch software (http://www.acdlabs.com/resources/freeware/ (accessed on 5 May 2021)). The structures were optimized and energy minimized using VEGAZZ software [66]. The optimized compounds were used to perform molecular docking to elucidate the anticancer activity against one fungal target protein which is acyl carrier protein domain from a fungal type I polyketide synthase (ACP), and two bacterial proteins including DNA gyrase (DNAg), and Baumannii penicillin-binding protein (PBP) to speculate the proposed antifungal and antibacterial activity for the most promising compounds. The three-dimensional structure of the molecular target was obtained from Protein Data Bank (PDB) (www.rcsb.org (accessed on 5 May 2021): (PBD: 2KR5, https://www.rcsb.org/structure/2KR5 (accessed on 5 May 2021), (PDB: 1KZN, https://www.rcsb.org/structure/1KZN (accessed on 5 May 2021), and (PDB: 3UDI, https://www.rcsb.org/structure/3UDI (accessed on 5 May 2021) respectively. The steps for receptor preparation included the removal of heteroatoms (water and ions), the addition of polar hydrogen, and the assignment of charge. The active sites were defined using grid boxes of appropriate sizes around the bound cocrystal ligands. The docking study was performed using Autodock vina [67] and Chimera for visualization [68]. All docking procedures and scoring were recorded according to our previous publications [69][70][71].

Pharmacokinetics and ADME Activity
To identify the biological targets for the most promising tested compound, we employed searching in database function integrated on Swiss institute Bioinformatics tools. Further on absorption, distribution, metabolism, and excretion "ADME" were calculated hypothetically.

Computational Method
The geometric parameters and energies were computed by density functional theory at the B3LYP/CEP-31G level of theory, using the GAUSSIAN 98W package of the programs [75], on geometries that were optimized at CEP-31G basis set. The high basis set was chosen to detect the energies at a highly accurate level. The atomic charges were computed using the natural atomic orbital populations. The B3LYP is the key word for the hybrid functional [76], which is a linear combination of the gradient functionals proposed by Becke [77] and Lee, Yang and Parr [78], together with the Hartree-Fock local exchange function [79]. UV spectra were recorded in Rigol, Ultra-3000 series in Enzymology and Fungal Biotechnology Lab, Faculty of Science, Zagazig University.

Conclusions
The corresponding xanthines 3-5, and imidazolone derivatives 6-13, were obtained via reaction of oxazolone derivative 2 with 5,6-diaminouracils 1 under various conditions in excellent yields. The condensation of 5,6-diaminouracils 1 with aldehyde afforded 5-Benzylideneaminouracils 14-16. The structural identities of the resulting compounds were resolved by IR, 1 H-, 13 C-NMR and Mass spectral analyses. The compounds 3, 6 and 16 displayed the highest activity against E. coli ((IC 50 value 1.8-1.9 µg/mL). The compound 16 displayed a significant antifungal activity against C. albicans (0.82 µg/mL), A. flavus (1.2 µg/mL) comparing to authentic antibiotics. This study paves the way for new generations of uracil and imidiazolone derivatives as potential antimicrobial and biologically active agents. Funding: The work has been funded by our own money.
Institutional Review Board Statement: Ethical review and approval were waived for this study, due to the absence of any direct work with animals or human.

Informed Consent Statement:
The study does not include any experiments related to human. Data Availability Statement: All the data are included on this manuscript as well as enclosed supplementary data.