Synthesis and Rational Design of New Appended 1,2,3-Triazole-uracil Ensembles as Promising Anti-Tumor Agents via In Silico VEGFR-2 Transferase Inhibition

Angiogenesis inhibition is a key step towards the designing of new chemotherapeutic agents. In a view to preparing new molecular entities for cancer treatment, eighteen 1,2,3-triazole-uracil ensembles 5a–r were designed and synthesized via the click reaction. The ligands were well characterized using 1H-, 13C-NMR, elemental analysis and ESI-mass spectrometry. The in silico binding propinquities of the ligands were studied sequentially in the active region of VEGFR-2 using the Molegro virtual docker. All the compounds produced remarkable interactions and potentially inhibitory ligands against VEGFR-2 were obtained with high negative binding energies. Drug-likeness was assessed from the ADME properties. Cytotoxicity of the test compounds was measured against HeLa and HUH-7 tumor cells and NIH/3T3 normal cells by MTT assay. Compound 5h showed higher growth inhibition activity than the positive control, 5-fluorouracil (5-FU), against both HeLa and HUH-7 cells with IC50 values of 4.5 and 7.7 μM respectively. Interestingly, the compounds 5a–r did not show any cytotoxicity towards the normal cell lines. The results advance the position of substituted triazoles in the area of drug design with no ambiguity.


Introduction
Cancer is a disease with a high mortality rate that accounts for more than 9 million deaths every year [1,2]. The discovery of molecules to treat cancer with minimal or no side effects is of great importance in medicinal chemistry. To date, over 100 anticancer agents are approved for chemotherapy [3]. Due to the increasing resistance to existing drugs, there is a need for the synthesis of new agents which can effectively treat various cancers. The present investigation attempts to develop new therapeutic agents assembled using 1,2,3-triazoles and uracil moieties. 1,2,3-Triazoles contain three nitrogens in the ring system and are an important class of heterocyclic compounds with remarkable biological properties [4]. It is known that they possess promising antiviral, antifungal and anticancer, antileishmanial, antituberculosis, anti-inflammatory, antimicrobial and

Chemistry
The target compounds 5a-r (Scheme 1) were synthesized using previously reported literature methods [38,39]. The required aromatic azides 2a-e and benzyl or alkyl azides 2f-r were first synthesized by the reaction of aromatic anilines 1a-e with dil. HCl and NaNO 2 in water at 0-5 • C, followed by the addition of sodium azide in water at room temperature. The benzyl and other alkyl bromides 1f-r were made to react with sodium azide to give azide reactants 2f-r by the replacement of the terminal bromide. Finally, Cu(I) catalyzed coupling reactions between 6-chloro-3-methyl-1-(prop-2-yn-1-yl)pyrimidine-2,4(1H,3H)-dione (4) and the aromatic/benzyl/alkyl azides 2a-r was achieved by refluxing in DMF at 80 • C for 12 h to produce the target compounds 5a-r. As a representative case, the spectral analysis of compound 5e is discussed below.
The formation of the substituted triazole-uracil compound 5e was confirmed by analyzing its 400 MHz 1 H-NMR spectrum (Supplementary Materials, Figure S10) which was recorded in CDCl 3 . The protons of the -CH 2 -group bridging between the uracil and triazole rings appeared at δ (ppm) 5.44 (s, 2H). The uracil ring proton was present at δ 6.02 (s, 1H) and phenyl ring protons appeared at δ 7.65 (m, 2H) and 7.54 (m, 2H). The methoxy protons of the phenyl ring and N-CH 3 protons of uracil appeared at 3.94 (s, 3H) and 3.34 (s, 3H), respectively. The triazole ring proton appeared at δ 8.08 (s, 1H). The 100.6 MHz 13 C-NMR spectrum of the compound was also recorded in CDCl 3 , (Supplementary Materials, Figure S11). The uracil carbonyl carbons appeared at δ 164.3 (C 2 ) and δ 150.5 (C 5 ). The triazole ring carbons appeared at δ 134.0 and 126.9. The four phenyl carbons were assigned to the peaks at δ 160.5 (C4 ), 131.9 (C1 ), 131.3 (C2 ), 114.7 (C3 ). The carbon of the -CH 2group between the uracil and triazole rings was observed at δ 36.1. The methoxy and N-CH 3 carbons appeared at δ 66.6 and δ 28.4, respectively. The peaks that appeared at m/z 348 and 349 in the mass spectrum of 5e ( Figure S37) were assigned to the [M + H] + and [M + 2] + molecular ion peaks, respectively.

In Silico Docking Study
Molecular docking is an essential tool to screen the binding affinities and mode of binding of protein-ligand interactions [40]. The docking analysis of the synthesized ligands ( Figure 2) 5a-r was performed using the Molegro virtual docker. A well-established molecular target, VEGFR-2 (PDB ID: 4AGD), was selected to study the docking interactions. The active site pocket of the receptor protein consists of Leu840, Gly922, Leu1038, Phe1047, Asp1046, Glu885, Val916, Glu917, Phe918, Cys919, Val848, Val899, Ile392, Leu889, Val916, Lys868, Gly841, Cys1045, Asp1052, Lys1055, Tyr1082, Arg1027, Leu802, Ser803, Leu1053, Pro1068, Tyr1054 and Arg842 amino acid residues. The results (Table 1) indicated that all the compounds bound to VEGFR-2 through hydrogen bonds and Van der Waal interactions with high binding affinities ranging from −108 to −198.7 kcal/mol. Figures 3-10 depict the best binding poses for the reference drug 5-FU and compounds 5a, 5f, 5h, 5i, 5j, 5n and 5r.  The reference drug, 5-FU produced a binding affinity of −65.73 kcal/mol. The interactions of the reference drug with the target protein are depicted in Figure 3. The -NH and carbonyl groups of the uracil moiety were stabilized by conventional hydrogen bonds with amino acid residues Glu885 (2.37 Å) and Asp1046 (2.03 Å), respectively. The fluorine atom of the uracil core interacted with the amino acid Val899. Further, the uracil nucleus is anchored by π-alkyl interactions with the amino acids Ile888, Ile804 and Leu889.  The reference drug, 5-FU produced a binding affinity of −65.73 kcal/mol. The interactions of the reference drug with the target protein are depicted in Figure 3. The -NH and carbonyl groups of the uracil moiety were stabilized by conventional hydrogen bonds with amino acid residues Glu885 (2.37 Å) and Asp1046 (2.03 Å), respectively. The fluorine atom of the uracil core interacted with the amino acid Val899. Further, the uracil nucleus is anchored by π-alkyl interactions with the amino acids Ile888, Ile804 and Leu889. Molecules 2021, 26, x 6 of 21 The synthesized compound 5a showed a binding affinity of −108.00 kcal/mol and the corresponding interactions are presented in Figure 4. The uracil core is anchored by πalkyl interactions with the Ile804, Val898, Ile892 and Leu889 residues. The chlorine atom of the uracil ring interacted with the carbonyl backbone of the Ile1044 amino acid. The methyl group present on the nitrogen of the uracil ring interacted with the carbonyl oxygen of Ile888. However, the carbonyl oxygens of the uracil ring did not form any kind of interactions with the target protein. The triazole ring is located in the hydrophobic pocket formed by Val899, Cys1045, Lys868 and Leu889 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended phenyl ring of 5a showed π-π interactions with Phe1047 and π-alkyl interactions with the Val848, Val916, Val899, Leu1035 and Cys1045 amino acid residues. Compound 5f produced a binding affinity of −112.05 kcal/mol and the corresponding interactions are presented in Figure 5. The uracil moiety is anchored by π-alkyl interactions with the Leu889, Ile804, Ile888 and Ile892 residues. The chlorine atom of the uracil ring interacted with the carboxyl group of the Glu885 amino acid. The methyl group present on the nitrogen of the uracil ring settled in the hydrophobic pocket formed by the Ile892, Val898 and Leu1019 residues. The triazole ring is located in the hydrophobic The synthesized compound 5a showed a binding affinity of −108.00 kcal/mol and the corresponding interactions are presented in Figure 4. The uracil core is anchored by π-alkyl interactions with the Ile804, Val898, Ile892 and Leu889 residues. The chlorine atom of the uracil ring interacted with the carbonyl backbone of the Ile1044 amino acid. The methyl group present on the nitrogen of the uracil ring interacted with the carbonyl oxygen of Ile888. However, the carbonyl oxygens of the uracil ring did not form any kind of interactions with the target protein. The triazole ring is located in the hydrophobic pocket formed by Val899, Cys1045, Lys868 and Leu889 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended phenyl ring of 5a showed π-π interactions with Phe1047 and π-alkyl interactions with the Val848, Val916, Val899, Leu1035 and Cys1045 amino acid residues.  The synthesized compound 5a showed a binding affinity of −108.00 kcal/mol and the corresponding interactions are presented in Figure 4. The uracil core is anchored by πalkyl interactions with the Ile804, Val898, Ile892 and Leu889 residues. The chlorine atom of the uracil ring interacted with the carbonyl backbone of the Ile1044 amino acid. The methyl group present on the nitrogen of the uracil ring interacted with the carbonyl oxygen of Ile888. However, the carbonyl oxygens of the uracil ring did not form any kind of interactions with the target protein. The triazole ring is located in the hydrophobic pocket formed by Val899, Cys1045, Lys868 and Leu889 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended phenyl ring of 5a showed π-π interactions with Phe1047 and π-alkyl interactions with the Val848, Val916, Val899, Leu1035 and Cys1045 amino acid residues. Compound 5f produced a binding affinity of −112.05 kcal/mol and the corresponding interactions are presented in Figure 5. The uracil moiety is anchored by π-alkyl interactions with the Leu889, Ile804, Ile888 and Ile892 residues. The chlorine atom of the uracil ring interacted with the carboxyl group of the Glu885 amino acid. The methyl group present on the nitrogen of the uracil ring settled in the hydrophobic pocket formed by the Ile892, Val898 and Leu1019 residues. The triazole ring is located in the hydrophobic Compound 5f produced a binding affinity of −112.05 kcal/mol and the corresponding interactions are presented in Figure 5. The uracil moiety is anchored by π-alkyl interactions with the Leu889, Ile804, Ile888 and Ile892 residues. The chlorine atom of the uracil ring interacted with the carboxyl group of the Glu885 amino acid. The methyl group present on the nitrogen of the uracil ring settled in the hydrophobic pocket formed by the Ile892, Val898 and Leu1019 residues. The triazole ring is located in the hydrophobic pocket formed Molecules 2021, 26,1952 7 of 20 by Val916, Leu889, Val899 and Cys1045 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended benzyl ring of 5f is stabilized via π-π interactions with Phe1047 and π-alkyl interactions the with amino acid residues Val848, Val899, Leu1035 and Cys1045.
Molecules 2021, 26, x 7 of 21 pocket formed by Val916, Leu889, Val899 and Cys1045 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended benzyl ring of 5f is stabilized via π-π interactions with Phe1047 and π-alkyl interactions the with amino acid residues Val848, Val899, Leu1035 and Cys1045. Compound 5h showed a binding affinity of −124.79 kcal/mol and the corresponding interactions are presented in Figure 6. It produced three conventional hydrogen bonds in the active region of the target protein. The carbonyl oxygen of the uracil ring formed a hydrogen bond with Cys919 with a bond distance of 1.97 Å. The two nitrogens at the 2 and 3-positons of the triazole ring formed hydrogen bonds with Asp1046 (2.06 Å) and Lys868 (2.54 Å), respectively. The uracil core is anchored by π-alkyl interactions with the Ala868, Val848, Cys919 and Leu1035 residues. The methyl group present on the nitrogen of the uracil ring interacted with the alkyl groups of the amino acids Leu840 and Leu1035 and it also formed a π-sigma interaction with the Phe1047 residue. The triazole ring is anchored by π-alkyl interactions with Val916, Cys1045 and Val899. Further, the triazole core is stabilized byπ-π and π-cation interactions with the amino acids Phe1047 and Lys868, respectively. The appended benzyl ring of 5h showed π-alkyl interactions with the amino acids Ile804, Leu889 and Val899. The methyl group present on this phenyl ring is located in the hydrophobic pocket formed by the Ile888, Ile804 and Leu889 residues.  Compound 5h showed a binding affinity of −124.79 kcal/mol and the corresponding interactions are presented in Figure 6. It produced three conventional hydrogen bonds in the active region of the target protein. The carbonyl oxygen of the uracil ring formed a hydrogen bond with Cys919 with a bond distance of 1.97 Å. The two nitrogens at the 2 and 3-positons of the triazole ring formed hydrogen bonds with Asp1046 (2.06 Å) and Lys868 (2.54 Å), respectively. The uracil core is anchored by π-alkyl interactions with the Ala868, Val848, Cys919 and Leu1035 residues. The methyl group present on the nitrogen of the uracil ring interacted with the alkyl groups of the amino acids Leu840 and Leu1035 and it also formed a π-sigma interaction with the Phe1047 residue. The triazole ring is anchored by π-alkyl interactions with Val916, Cys1045 and Val899. Further, the triazole core is stabilized byπ-π and π-cation interactions with the amino acids Phe1047 and Lys868, respectively. The appended benzyl ring of 5h showed π-alkyl interactions with the amino acids Ile804, Leu889 and Val899. The methyl group present on this phenyl ring is located in the hydrophobic pocket formed by the Ile888, Ile804 and Leu889 residues. pocket formed by Val916, Leu889, Val899 and Cys1045 with π-alkyl interactions. The hydrogen atom of the triazole ring is bonded to the carbonyl oxygen of the Asp1046 residue. Further, the triazole moiety is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The appended benzyl ring of 5f is stabilized via π-π interactions with Phe1047 and π-alkyl interactions the with amino acid residues Val848, Val899, Leu1035 and Cys1045. Compound 5h showed a binding affinity of −124.79 kcal/mol and the corresponding interactions are presented in Figure 6. It produced three conventional hydrogen bonds in the active region of the target protein. The carbonyl oxygen of the uracil ring formed a hydrogen bond with Cys919 with a bond distance of 1.97 Å. The two nitrogens at the 2 and 3-positons of the triazole ring formed hydrogen bonds with Asp1046 (2.06 Å) and Lys868 (2.54 Å), respectively. The uracil core is anchored by π-alkyl interactions with the Ala868, Val848, Cys919 and Leu1035 residues. The methyl group present on the nitrogen of the uracil ring interacted with the alkyl groups of the amino acids Leu840 and Leu1035 and it also formed a π-sigma interaction with the Phe1047 residue. The triazole ring is anchored by π-alkyl interactions with Val916, Cys1045 and Val899. Further, the triazole core is stabilized byπ-π and π-cation interactions with the amino acids Phe1047 and Lys868, respectively. The appended benzyl ring of 5h showed π-alkyl interactions with the amino acids Ile804, Leu889 and Val899. The methyl group present on this phenyl ring is located in the hydrophobic pocket formed by the Ile888, Ile804 and Leu889 residues.  Compound 5i showed a binding affinity of −118.52 kcal/mol and the corresponding interactions are depicted in Figure 7. It showed three conventional hydrogen bonds in the active site of the receptor protein. The nitrogen at the 2-positon of the triazole ring formed a hydrogen bond with Asn923 with a bond distance of 2.48 Å. The oxygen of the appended benzyl group produced two hydrogen bonds with Arg1051 forming a stable six-membered ring with bond distances of 2.16 and 2.37 Å. The chlorine atom of the uracil moiety produced four hydrophobic interactions with the amino acids Leu840, Cys919, Phe918 and Ala866. The uracil core is anchored by π-alkyl interactions with Leu1035, Ala866 and Val848. N-methyl group of uracil ring interacted with Val848 and Phe1047 via alkyl and π-alkyl interactions. The triazole ring is stabilized by Leu1035 and Leu840 though π-alkyl interactions and Phe1047 with π-π interactions. The appended benzyl ring showed a π-π interaction with the amino acid Phe1047. Compound 5i showed a binding affinity of −118.52 kcal/mol and the corresponding interactions are depicted in Figure 7. It showed three conventional hydrogen bonds in the active site of the receptor protein. The nitrogen at the 2-positon of the triazole ring formed a hydrogen bond with Asn923 with a bond distance of 2.48 Å. The oxygen of the appended benzyl group produced two hydrogen bonds with Arg1051 forming a stable six-membered ring with bond distances of 2.16 and 2.37 Å. The chlorine atom of the uracil moiety produced four hydrophobic interactions with the amino acids Leu840, Cys919, Phe918 and Ala866. The uracil core is anchored by π-alkyl interactions with Leu1035, Ala866 and Val848. N-methyl group of uracil ring interacted with Val848 and Phe1047 via alkyl and π-alkyl interactions. The triazole ring is stabilized by Leu1035 and Leu840 though π-alkyl interactions and Phe1047 with π-π interactions. The appended benzyl ring showed a π-π interaction with the amino acid Phe1047. Compound 5j produced a binding affinity of −152.32 kcal/mol and the corresponding interactions are presented in Figure 8. It produced two hydrogen bonds in the active region of the target protein. The carbonyl oxygen of the uracil ring formed two hydrogen bonds with Cys919 and Phe918 with bond distances of 1.61 and 2.93 Å, respectively. The uracil core is anchored by π-alkyl interactions with Ala866, Val848, leu1035, Cys919 and Lue840. The N-methyl group of the uracil ring interacted with the Leu840 and Cys919 residues via alkyl interactions. The triazole ring is stabilized by π-alkyl interactions with Cys1045 and Val848 and π-π interactions with Phe1047. The appended benzyl ring showed a π-alkyl interaction with the amino acids Val916, Val899 and Leu889. Further, the benzyl ring is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The phenoxy ring produced π-alkyl interactions with the amino acids Ile804, Leu889, Ile888 and Ile892. Compound 5j produced a binding affinity of −152.32 kcal/mol and the corresponding interactions are presented in Figure 8. It produced two hydrogen bonds in the active region of the target protein. The carbonyl oxygen of the uracil ring formed two hydrogen bonds with Cys919 and Phe918 with bond distances of 1.61 and 2.93 Å, respectively. The uracil core is anchored by π-alkyl interactions with Ala866, Val848, leu1035, Cys919 and Lue840. The N-methyl group of the uracil ring interacted with the Leu840 and Cys919 residues via alkyl interactions. The triazole ring is stabilized by π-alkyl interactions with Cys1045 and Val848 and π-π interactions with Phe1047. The appended benzyl ring showed a π-alkyl interaction with the amino acids Val916, Val899 and Leu889. Further, the benzyl ring is stabilized by π-cation and π-anion interactions with the amino acids Lys868 and Glu885, respectively. The phenoxy ring produced π-alkyl interactions with the amino acids Ile804, Leu889, Ile888 and Ile892. Molecules 2021, 26, x 9 of 21 Compound 5n produced a binding affinity of −125 kcal/mol and the corresponding interactions are presented in Figure 9. It produced one conventional hydrogen bond with the target protein. The carbonyl oxygen of the isoquinoline ring formed one hydrogen bond with Leu802 with a bond distance of 1.96 Å. The uracil core is stabilized by a π-alkyl interaction with the amino acid Leu1067. The uracil chlorine atom produced two hydrophobic interactions with the Met1072 and Pro1068 amino acids and it also interacted with the carbonyl backbone of Arg1027. The triazole ring is stabilized by π-anionic interactions with the carboxylic acid group of Asp1028. The isoquinoline ring is anchored by π-alkyl interactions with the amino acids Ile1053, Leu1049 and Ala844.  Figure 10. In this molecule, the uracil and triazole rings are anchored by π-π stacking interactions between uracil-triazole, triazole-triazole rings and uracil-Phe1047 amino acid residue of the receptor protein. The carbonyl oxygen of two uracil rings interacted with the side chain hydrogens of the amino acids Gly922 and Gly841. Further, the uracil and triazole rings are stabilized by the amino acids Leu1035, Leu840, Val848, Lys868, Val916 and Val899 through π-alkyl interactions. Compound 5n produced a binding affinity of −125 kcal/mol and the corresponding interactions are presented in Figure 9. It produced one conventional hydrogen bond with the target protein. The carbonyl oxygen of the isoquinoline ring formed one hydrogen bond with Leu802 with a bond distance of 1.96 Å. The uracil core is stabilized by a π-alkyl interaction with the amino acid Leu1067. The uracil chlorine atom produced two hydrophobic interactions with the Met1072 and Pro1068 amino acids and it also interacted with the carbonyl backbone of Arg1027. The triazole ring is stabilized by π-anionic interactions with the carboxylic acid group of Asp1028. The isoquinoline ring is anchored by π-alkyl interactions with the amino acids Ile1053, Leu1049 and Ala844. Compound 5n produced a binding affinity of −125 kcal/mol and the corresponding interactions are presented in Figure 9. It produced one conventional hydrogen bond with the target protein. The carbonyl oxygen of the isoquinoline ring formed one hydrogen bond with Leu802 with a bond distance of 1.96 Å. The uracil core is stabilized by a π-alkyl interaction with the amino acid Leu1067. The uracil chlorine atom produced two hydrophobic interactions with the Met1072 and Pro1068 amino acids and it also interacted with the carbonyl backbone of Arg1027. The triazole ring is stabilized by π-anionic interactions with the carboxylic acid group of Asp1028. The isoquinoline ring is anchored by π-alkyl interactions with the amino acids Ile1053, Leu1049 and Ala844. Of all the docked compounds, 5r showed the highest binding affinity (−198 kcal/mol) with the receptor protein. The interactions of the molecule 5r with the target protein are depicted in Figure 10. In this molecule, the uracil and triazole rings are anchored by π-π stacking interactions between uracil-triazole, triazole-triazole rings and uracil-Phe1047 amino acid residue of the receptor protein. The carbonyl oxygen of two uracil rings interacted with the side chain hydrogens of the amino acids Gly922 and Gly841. Further, the uracil and triazole rings are stabilized by the amino acids Leu1035, Leu840, Val848, Lys868, Val916 and Val899 through π-alkyl interactions. Of all the docked compounds, 5r showed the highest binding affinity (−198 kcal/mol) with the receptor protein. The interactions of the molecule 5r with the target protein are depicted in Figure 10. In this molecule, the uracil and triazole rings are anchored by π-π stacking interactions between uracil-triazole, triazole-triazole rings and uracil-Phe1047 amino acid residue of the receptor protein. The carbonyl oxygen of two uracil rings interacted with the side chain hydrogens of the amino acids Gly922 and Gly841. Further, the uracil and triazole rings are stabilized by the amino acids Leu1035, Leu840, Val848, Lys868, Val916 and Val899 through π-alkyl interactions.

ADME Properties
A favorable absorption, distribution, metabolism and elimination profile is a prerequisite for any drug candidate to conclude its drug-likeness. The ADME properties for the reported ligands, were predicted using the Lipinski filters from the Molinspiration online tool kit. The properties like molecular weight, molecular volume, the logarithm of partition coefficient, hydrogen bond acceptors and donors, polar surface area, the number of rotatable bonds and Lipinski's rule of five were determined by using the Molinspiration online property calculation toolkit [41]. Absorption percentage was obtained by: %ABS = 109 − (0.345 × TPSA) [42].
The Molsoft software was utilized to calculate properties like physico-chemical parameters, pharmacokinetics, pharmacodynamics (drug-likeness model score) and was represented by a numerical value [43]. The pharmacokinetic parameters for the synthesized ligands are presented in Table 2.

ADME Properties
A favorable absorption, distribution, metabolism and elimination profile is a prerequisite for any drug candidate to conclude its drug-likeness. The ADME properties for the reported ligands, were predicted using the Lipinski filters from the Molinspiration online tool kit. The properties like molecular weight, molecular volume, the logarithm of partition coefficient, hydrogen bond acceptors and donors, polar surface area, the number of rotatable bonds and Lipinski's rule of five were determined by using the Molinspiration online property calculation toolkit [41]. Absorption percentage was obtained by: %ABS = 109 − (0.345 × TPSA) [42].
The Molsoft software was utilized to calculate properties like physico-chemical parameters, pharmacokinetics, pharmacodynamics (drug-likeness model score) and was represented by a numerical value [43]. The pharmacokinetic parameters for the synthesized ligands are presented in Table 2.
The results indicated good %Abs values for the triazole ligands ranging from 57.43% to 83.21%. Furthermore, compounds 5a-r obeyed Lipinski's rule of five (number of hydrogen bond acceptors (n-ON) ≤ 10) and followed the criteria for orally active drug and therefore, these compounds may have a good potential for eventual development as oral agents. A molecule likely to be developed as an orally active drug candidate should not produce more than one violation of the following four criteria: m i LogP (octanol-water partition coefficient) ≤ 5, molecular weight ≤ 500, number of hydrogen bond acceptors ≤ 10 and number of hydrogen bond donors ≤ 5 [44]. All the synthesized triazole-uracil analogues 5a-r obeyed the criteria and hence, can be considered for the development as orally active drugs.

In Vitro Cytotoxicity
The cytotoxic potency for 5a-r was assayed against HeLa and HUH-7 tumor cells by the standard MTT protocol. The ligands significantly inhibited the proliferation of the selected tumor cells and the corresponding viability profiles are presented in Figure 11. produce more than one violation of the following four criteria: miLogP (octanol-water partition coefficient) ≤ 5, molecular weight ≤ 500, number of hydrogen bond acceptors ≤ 10 and number of hydrogen bond donors ≤ 5 [44]. All the synthesized triazole-uracil analogues 5a-r obeyed the criteria and hence, can be considered for the development as orally active drugs.

In Vitro Cytotoxicity
The cytotoxic potency for 5a-r was assayed against HeLa and HUH-7 tumor cells by the standard MTT protocol. The ligands significantly inhibited the proliferation of the selected tumor cells and the corresponding viability profiles are presented in Figure 11. The half-maximal inhibitory concentrations (IC50) for compounds 5a-r are presented in Table 3. The ligands 5a, 5h, 5i, 5j, 5n and 5r produced higher inhibition effect against HeLa cells, with IC50 values ranging from 4.5 μM to 11.2 μM. HUH-7 cells were found to be more sensitive towards the treatment of 5f, 5h and 5i with IC50 values 10.8 μM, 7.7 μM and 9 μM, respectively, and the remaining compounds failed to hamper the proliferation of HUH-7 tumor cells. Compound 5h emerged as a hit compound and produced IC50 values of 4.5 μM and 7.7 μM, respectively, against both the cell lines. The obtained values for 5h are better than those of the standard drug, 5-FU (4.6 μM and 30 μM) and it has the The half-maximal inhibitory concentrations (IC 50 ) for compounds 5a-r are presented in Table 3. The ligands 5a, 5h, 5i, 5j, 5n and 5r produced higher inhibition effect against HeLa cells, with IC 50 values ranging from 4.5 µM to 11.2 µM. HUH-7 cells were found to be more sensitive towards the treatment of 5f, 5h and 5i with IC 50 values 10.8 µM, 7.7 µM and 9 µM, respectively, and the remaining compounds failed to hamper the proliferation of HUH-7 tumor cells. Compound 5h emerged as a hit compound and produced IC 50 values of 4.5 µM and 7.7 µM, respectively, against both the cell lines. The obtained values for 5h are better than those of the standard drug, 5-FU (4.6 µM and 30 µM) and it has the potential to be developed as a new chemotherapeutic agent for cervical and liver carcinoma. Table 3. IC 50 values (µM) for the new 1,2,3-triazole-uracil scaffolds 5a-r against the treatment of HeLa and HUH-7 tumor cells and NIH/3T3 normal cells after 72 h.
For all the compounds 5a-r, the 1,2,3-triazole and uracil rings are common, and hence some structure-activity relationship inferences can be drawn based on the different substituents. The compounds 5b-e bearing a phenyl ring with electron-donating or withdrawing groups at the para position showed lower activity than the unsubstituted compound 5a. The compounds 5f-j (except 5g) bearing a benzyl ring and with various electron-donating groups produced higher activity against HeLa or HUH-7 cells. The compounds 5k-m with chromene rings failed to show any remarkable activity against the selected tumor cell lines. The compound 5n with an isoquinoline ring showed good activity against HeLa cells. Compounds 5o and 5p with aliphatic hydrocarbon chains did not induce considerable activity against either of the cancer cell lines. The compound 5r produced a higher inhibition effect than its dimer analogue 5q. Interestingly, the compounds 5a-r did not show any cytotoxicity towards the normal NIH/3T3 cells and this offer a great future for the reported ligands in cancer research.
On the other hand, the approved drugs for the treatment of liver cancer are sparse and sorafenib (VEGFR inhibitor) is the only drug approved to treat patients with advanced hepatocellular carcinoma [45][46][47]. Even though it is effective in treating liver cancer, patients are reported to develop resistance within six months and no alternative treatment is currently available [48,49]. Thus, the current results, that are close to the in vitro cytotoxicity of sorafenib [22,28] against liver carcinoma, are exciting and the reported ligands 5f, 5h and 5i could be developed as new alternatives for treating liver cancer.

General Information
All the chemicals used in this study were obtained from Sigma Chemicals (St. Louis, MO, USA). Solvents of reagent grade were purchased from nearby commercial sources and distilled before the usage. DMEM media and the serum supplement FBS for the cytotoxicity assay were obtained from Sigma Chemicals and Gibco (Hyderabad, TS, India), respectively. The melting points of the ligands were measured in open capillaries using the electrothemal melting point apparatus. 1 H-and 13 C-NMR spectra were recorded using Avance II 400 and 100 MHz spectrometers (Bruker, Zurich, Switzerland) in 99.99% DMSOd 6 and 99.82% CDCl 3 (2 mg/mL) using TMS as the standard solvent. ESI-MS spectra were measured on a LCMS 2010 VG mass spectrometer (Thermo Finnigan, San Jose, CA, USA) in methanol/DCM solvents (10 µg/mL). Elemental analysis was carried out using a Carlo Erba 1108 analyzer (Heraeus, Hanau, Germany).

Syntheses
3.2.1. Synthesis of 6-Chloro-3-methyl-1-(prop-2-yn-1-yl)pyrimidine-2,4(1H,3H)-dione (4) A DMF (50 mL) solution of 3 (1 g, 0.0062 mol) was slowly added with propargyl bromide (1.33 mL, 0.015 mol) and anhydrous K 2 CO 3 (2.56 g, 0.0186 mol). The resulting mixture was allowed to reflux for 4 h and the course of the reaction was monitored by TLC until its completion. When the reaction was complete, the remaining solvent was evaporated under vacuum. Excess propargyl bromide was removed by washing the reaction mixture with hexane and then poured into ice-water (100 mL) and extracted with ethyl acetate. Then, the organic layer was added with the brine solution and residual water content after workup was removed by sodium sulphate.

. Preparation of Ligands and Target Protein
The structures of all the ligands 5a-r were drawn using the ChemDraw suite (Ver.14.0.0.117, PerkinElmer Inc., Waltham, MA, USA) [50] and energy minimization was accomplished by OPLS 2005 through the Ligprep module of the Schrodinger suite (Schrödinger, New York, NY, USA) [51]. Ionization states for ligands were retained at their original values and saved in .sdf format for docking simulations. The crystal structure of the VEGFR-2 protein was retrieved from the protein data bank (PDB ID: 4AGD). Molecular docking operations in the active site of receptor protein were executed by the Molegro virtual docking program (MVD 6.0, 2013, Molexus, Odder, Denmark). All the solvent molecules attached to the receptor protein were removed and structural parameters like hybridization state, bond order, explicit hydrogen atoms and relevant charges for ligands were assigned according to the docking software instructions. Five potential cavities for binding were obtained by using the detect cavities option from the preparation tools. The cavity, surrounding the anion binding site with a volume of 177 Å 3 , was further modified by side-chain minimization and used for virtual screening. Moldock score (grid-based) with a grid resolution of 0.20 Å was used as a scoring function for docking analysis [52]. Moldock and Rerank scores estimated the best binding poses for ligands. ALIGN program from COOT graphical program was used for molecular alignments [53].

Molecular Docking
MVD is a convenient and fast algorithm to estimate protein and ligand interactions. The virtual screening of synthesized ligands depends on Rerank score which is a mathematical expression for protein-ligand binding affinity and is based on the Moldock score derived from PLP (piecewise linear potential) [14]. Total energy for the protein-ligand system was minimized by Medler Mead simplex minimization with H bond directionality and non-grid force field [54]. Binding affinities of ligands with protein were computed based on internal electrostatic, hydrogen bond interactions and sp 2 -sp 2 tortions. All the docking interactions of the protein-ligand system were analyzed using Accelrys Discovery Studio v3.5 (Waltham, MA, USA) [55].

ADME Profiles
Drug-likeness of the derived ligands was assessed by Lipinski filters in the Molinspiration webserver ( © Molinspiration Cheminformatics 2018, Nova Ulica, Slovensky Grob, Slovak Republic).

In Vitro Cytotoxicity
The cell viability was measured by MTT colourimetric assay. The cell lines viz., HeLa (human cervical adenocarcinoma), HUH-7 (human hepatoma cancer cells) and NIH/3T3 (normal murine fibroblast cells) were supplemented with 10% fetal bovine serum and maintained in an environment of 5% CO 2 at 37 • C temperature. Cell lines at a density of 2 × 10 3 cells/well were seeded onto 96 well plates for 16 h and maintained in DMEM culture media. Five different concentrations of synthesized ligands and 5-fluorouracil (except for normal cells) i.e., 0.1, 0.5, 1, 5 and 10 µM were treated against each cell line and incubated for 3 days at 37 • C [56]. Then, the solutions were discarded and an MTT reagent with a final concentration of 500 µg/mL was added to each well and plate and incubated for another 4-6 h [57][58][59]. 200 µL of DMSO was added to solubilize purple precipitates (formazan product) and absorbance (optical density) was recorded at 570 nm on ELISA reader. IC 50 values were calculated using dose-response data by applying the linear regression method. The experiments were done in triplicate for each concentration of ligand.

Conclusions
Herein, eighteen new 1,2,3-triazole-uracil scaffolds have been designed, synthesized and tested as anti-cancer agents. The structures of the compounds were established by standard analytical techniques viz., 1 H-NMR, 13 C-NMR, elemental analysis and mass spectrometry. The analytical data were in good agreement with the molecular formulae of the compounds. A rational approach of drug design was adopted to present the anti-tumor efficacies of the synthesized ligands 5a-r. To realize the objective, the ligands 5a-r were sequentially screened in the active site of VEGFR-2 protein and the results depicted remarkable interactions with the receptor protein. The molecular docking simulations by Molegro virtual docker produced high binding affinities for all the compounds (−108 to −198.7 kcal/mol) with various hydrophobic, electrostatic and hydrogen bond interactions. The interactions observed for the compounds 5a, 5f, 5h, 5i, 5j, 5n and 5r determine proximity with receptor protein with different binding modes and justify their remarkable cytotoxic activity among all the synthesized compounds. The propinquity of ligands with receptor protein presumed the ability to inhibit in vitro cell proliferation and hence the ligands were screened for cell viability by the MTT method. Six compounds 5a, 5h, 5i, 5j, 5n and 5r were found to inhibit cell proliferation with low half-maximal inhibitory concentrations (IC 50 ) of 11.2 µM, 4.5 µM, 9.0 µM, 10 µM, 9.6 µM, 9.6 µM, respectively. Compound 5h emerged as a hit compound by producing higher inhibition activity than 5-FU against both the HeLa and HUH-7 tumor cell lines with IC 50 values 4.5 and 7.7 µM, respectively. The docking results were very well correlated to the experimental results. Interestingly, the compounds didn't show any cytotoxicity towards normal cells, paving the way for the reported ligands to be developed as new alternatives to existing drugs for cervical and liver cancers.  5b, 5c, 5d, 5e, 5g, 5k, 5l, 5m, 5o,