Regioselective Synthesis and Molecular Docking Studies of 1,5-Disubstituted 1,2,3-Triazole Derivatives of Pyrimidine Nucleobases

1,2,3-triazoles are versatile building blocks with growing interest in medicinal chemistry. For this reason, organic chemistry focuses on the development of new synthetic pathways to obtain 1,2,3-triazole derivatives, especially with pyridine moieties. In this work, a novel series of 1,5-disubstituted-1,2,3-triazoles functionalized with pyrimidine nucleobases were prepared via 1,3-dipolar cycloaddition reaction in a regioselective manner for the first time. The N1-propargyl nucleobases, used as an alkyne intermediate, were obtained in high yields (87–92%) with a new two-step procedure that selectively led to the monoalkylated compounds. Then, FeCl3 was employed as an efficient Lewis acid catalyst for 1,3-dipolar cycloaddition between different aryl and benzyl azides and the N1-propargyl nucleobases previously synthesized. This new protocol allows the synthesis of a series of new 1,2,3-triazole derivatives with good to excellent yields (82–92%). The ADME (Absorption, Distribution, Metabolism, and Excretion) analysis showed good pharmacokinetic properties and no violations of Lipinsky’s rules, suggesting an appropriate drug likeness for these new compounds. Molecular docking simulations, conducted on different targets, revealed that two of these new hybrids could be potential ligands for viral and bacterial protein receptors such as human norovirus capsid protein, SARS-CoV-2 NSP13 helicase, and metallo-β-lactamase.

Other transition metal catalysts (Au, Ir, Ni, Ag) were tested [15], but despite the huge efforts, the employment of heavy metals in the synthesis restricts their application due to their hazardous nature, toxicity, and high cost. In recent years, great attention was paid to Other transition metal catalysts (Au, Ir, Ni, Ag) were tested [15], but despite the huge efforts, the employment of heavy metals in the synthesis restricts their application due to their hazardous nature, toxicity, and high cost. In recent years, great attention was paid to the development of copper-and ruthenium-based alternative routes, and different protocols were developed [16]. Besides 1,2,3-triazoles, natural and modified nucleosides are another class of heterocyclic compounds with widespread applications as pharmaceutically active compounds [17,18] and diagnostic probes [19]. Considering the numerous advantages of using both triazoles and nucleosides as building blocks in drug discovery, many researchers synthesized nucleosides containing the triazole moiety [20,21]. Starting from the discovery of the broad-spectrum antiviral compound Ribavirin I (Figure 1), in which a 1,2,4-triazole was used as a nucleobase, different new hybrids were designed and tested, showing antiviral and/or antitumor activity. These nucleoside analogues were characterized by the introduction of a triazole group in place of the nucleobase I [22], or the sugar moiety II [23], as a linkage between them III [24], and as a modification group of the sugar IV [25] (Figure 1). To date, only a few examples of 1,4-disubstituted 1,2,3-triazole nucleoside derivatives were synthesized. These compounds were tested as corrosion inhibitors for steel [26], against HCV (Hepatitis C virus) [23], and against influenza virus A (H3N2) [27]. However, to the best of our knowledge, there have been no reported synthetic efforts toward 1,5disubstituted analogues.
In this work, considering our experience in the development of alternative synthetic routes catalyzed by Lewis acids [28,29] and cycloaddition reactions [30,31], we report the synthesis of new 1,5-disubstituted 1,2,3-triazole derivatives containing pyrimidine nucleobases in high yields and a regioselective way. All of the triazole derivatives were synthesized starting from different aryl azides and non-commercial N1-propargyl pyrimidine nucleobases using a common Lewis acid catalyst. In detail, iron(III) chloride was selected as a suitable catalyst because of its low price, easy availability, sustainability, nontoxicity, and environmentally friendly characteristic [32,33]. Moreover, the use of N1-propargyl pyrimidine nucleobases allows the introduction of a methylene bridge at C-5 of the triazole, which could mimic the behavior of "fleximers", a kind of nucleoside analogue characterized by enhanced conformational freedom [34]. The pharmacokinetic properties of the synthesized derivatives were predicted through ADME (Absorption, Distribution, Metabolism, and Excretion) analysis to evaluate the medicinal chemistry friendliness. Furthermore, an in silico screening was performed towards selected receptors from the Protein Data Bank, suggesting biological potential for our products. To date, only a few examples of 1,4-disubstituted 1,2,3-triazole nucleoside derivatives were synthesized. These compounds were tested as corrosion inhibitors for steel [26], against HCV (Hepatitis C virus) [23], and against influenza virus A (H3N2) [27]. However, to the best of our knowledge, there have been no reported synthetic efforts toward 1,5-disubstituted analogues.
In this work, considering our experience in the development of alternative synthetic routes catalyzed by Lewis acids [28,29] and cycloaddition reactions [30,31], we report the synthesis of new 1,5-disubstituted 1,2,3-triazole derivatives containing pyrimidine nucleobases in high yields and a regioselective way. All of the triazole derivatives were synthesized starting from different aryl azides and non-commercial N1-propargyl pyrimidine nucleobases using a common Lewis acid catalyst. In detail, iron(III) chloride was selected as a suitable catalyst because of its low price, easy availability, sustainability, nontoxicity, and environmentally friendly characteristic [32,33]. Moreover, the use of N1-propargyl pyrimidine nucleobases allows the introduction of a methylene bridge at C-5 of the triazole, which could mimic the behavior of "fleximers", a kind of nucleoside analogue characterized by enhanced conformational freedom [34]. The pharmacokinetic properties of the synthesized derivatives were predicted through ADME (Absorption, Distribution, Metabolism, and Excretion) analysis to evaluate the medicinal chemistry friendliness. Furthermore, an in silico screening was performed towards selected receptors from the Protein Data Bank, suggesting biological potential for our products.

Chemistry
N1-propargyl nucleobases are important starting materials for the synthesis of nucleoside analogues with biological activity [35]. Usually, propargyl nucleobases are synthesized in a one-step reaction between the appropriate nucleobase and propargyl bromide under basic conditions [36] or employing an intermediate bis(trimethylsilyl)pyrimidine nucleobase using N,O-bis(trimethylsilyl)-acetamide (BSA) [37]. Unfortunately, the first methodology yields the products with low selectivity [38] because a mixture composed of N1-monoalkylated and N1,N3-dialkylated pyrimidines is obtained, while the second approach results in long reaction times and low yield [26]. With this in mind, selective nucleobase propargylation is necessary, and, due to the feasibility of performing selective alkylation at the N-1 position of pyrimidine nucleobases, the propargylation is performed in two-reaction steps via O-protection by a transient group, as shown in Scheme 1.

Chemistry
N1-propargyl nucleobases are important starting materials for the synthesis of nucleoside analogues with biological activity [35]. Usually, propargyl nucleobases are synthesized in a one-step reaction between the appropriate nucleobase and propargyl bromide under basic conditions [36] or employing an intermediate bis(trimethylsilyl)pyrimidine nucleobase using N,O-bis(trimethylsilyl)-acetamide (BSA) [37]. Unfortunately, the first methodology yields the products with low selectivity [38] because a mixture composed of N1-monoalkylated and N1,N3-dialkylated pyrimidines is obtained, while the second approach results in long reaction times and low yield [26]. With this in mind, selective nucleobase propargylation is necessary, and, due to the feasibility of performing selective alkylation at the N-1 position of pyrimidine nucleobases, the propargylation is performed in two-reaction steps via O-protection by a transient group, as shown in Scheme 1. Here, we developed a methodology that involves the N1-selective propargylation of bis(trimethylsilyl)pyrimidine intermediates (4-6) to prepare compounds 7-9. Thus, pyrimidine nucleobases 1-3 were treated under an inert atmosphere with hexamethyldisilazane (HMDS), trimethylsilyl chloride (TMS-Cl), and (NH4)2SO4 to give the silylated nucleobases 4-6, which were used without further purification and isolated by under vacuum evaporation of HMDS. In situ propargylation of 4-6, conducted in dry acetonitrile at reflux under an inert atmosphere, furnished the desired products 7-9 with excellent reaction yields and in a regioselective way because only N-1 monosubstituted nucleobases were observed in all cases.
To begin, we chose the cycloaddition reaction between 1-propargylthymine 7 and phenyl azide 10 in the presence of Lewis acid catalysts as the model system to optimize the reaction conditions (Scheme 2). The results are reported in Table 1.
To begin, we chose the cycloaddition reaction between 1-propargylthymine 7 and phenyl azide 10 in the presence of Lewis acid catalysts as the model system to optimize the reaction conditions (Scheme 2). The results are reported in Table 1.

Chemistry
N1-propargyl nucleobases are important starting materials for the synthesis of nucleoside analogues with biological activity [35]. Usually, propargyl nucleobases are synthesized in a one-step reaction between the appropriate nucleobase and propargyl bromide under basic conditions [36] or employing an intermediate bis(trimethylsilyl)pyrimidine nucleobase using N,O-bis(trimethylsilyl)-acetamide (BSA) [37]. Unfortunately, the first methodology yields the products with low selectivity [38] because a mixture composed of N1-monoalkylated and N1,N3-dialkylated pyrimidines is obtained, while the second approach results in long reaction times and low yield [26]. With this in mind, selective nucleobase propargylation is necessary, and, due to the feasibility of performing selective alkylation at the N-1 position of pyrimidine nucleobases, the propargylation is performed in two-reaction steps via O-protection by a transient group, as shown in Scheme 1.

Scheme 1. Synthesis of propargyl nucleobases with two reaction steps.
Here, we developed a methodology that involves the N1-selective propargylation of bis(trimethylsilyl)pyrimidine intermediates (4-6) to prepare compounds 7-9. Thus, pyrimidine nucleobases 1-3 were treated under an inert atmosphere with hexamethyldisilazane (HMDS), trimethylsilyl chloride (TMS-Cl), and (NH4)2SO4 to give the silylated nucleobases 4-6, which were used without further purification and isolated by under vacuum evaporation of HMDS. In situ propargylation of 4-6, conducted in dry acetonitrile at reflux under an inert atmosphere, furnished the desired products 7-9 with excellent reaction yields and in a regioselective way because only N-1 monosubstituted nucleobases were observed in all cases.
To begin, we chose the cycloaddition reaction between 1-propargylthymine 7 and phenyl azide 10 in the presence of Lewis acid catalysts as the model system to optimize the reaction conditions (Scheme 2). The results are reported in Table 1.
Firstly, the reaction between 1-propargylthymine 7 and phenyl azide 10, in a 1:2 ratio, was performed with 20 mol% Er(OTf) 3 as the catalyst and heated at 60 • C both in CH 2 Cl 2 ( Table 1, entry 1) and THF (Table 1, entry 2). In both cases, no product formation was observed after 24 h of time reaction. An increase in temperature to 100 • C favored a slight formation of the product with a yield of 15% (  [8][9][10][11] and the best reaction conditions were found when FeCl 3 was used as the catalyst (Table 1, entry 11), leading to 88% yield in 8 h. The role of the catalyst was remarkable since it selectively led to one regioisomer, corresponding to the 1,5-disubstituted 1,2,3-triazole 14a. However, the regioisomer 14b was also isolated and characterized by 1 H NMR and 13 C NMR, and an 88:12 ratio was calculated by 1 H NMR analysis of the reaction crude (see Supplementary Materials). In fact, without any catalyst, reagents 7 and 10 under the same reaction conditions gave the 1,5-disubstituted triazoles 14 in low yield after a long reaction time (Table 1, entry 12) and a mixture of the 1,5-disubstituted/1,4-disubstituted products was observed in a 63:37 ratio. This last result suggests that the Lewis acid catalyst accelerates the reaction by increasing the electrophilicity of the alkyne group through coordination. Although the iron(III) chloride efficacy as a catalyst in an eliminative azide-olefin cycloaddition was already demonstrated for the synthesis of 1,5-disubstituted triazoles [39], we were delighted to observe such regioselectivity. Moreover, a comparison with classical cycloaddition reaction conditions was made in toluene without any catalyst. A regioisomer mixture in a 60:40 ratio was obtained in low yield after a long reaction time (Table 1, entry 13). The addition of a catalyst slightly increased the yield without significantly improving regioselectivity (Table 1, entry  14). Finally, the 1,3-dipolar cycloaddition reaction was carried out in [mPy](OTf) ionic liquid ( Table 1, entry 15), prepared as reported in the literature [40], with the purpose of hypothetically recycling the solvent. Unfortunately, due to the high solubility of the pyrimidine in the ionic liquid, it was not possible to use the latter as a solvent because it was impossible to purify the product by simple liquid-liquid extraction. Once the reaction conditions were optimized, we extended the protocol to various N1-propargyl nucleobases 7-9 and different azides 10-13 (Scheme 3, Table 2).
Aromatic azides 10, 12, and 13 exhibited different reactivities depending on the substitution of the aromatic ring. In fact, the presence of a strong electron-withdrawing group (-NO 2 ) at the para position reduced the nucleophilicity of the azide 12 (Table 2, entries 7-9), while the presence of the electron-donor group (CH 3 O − ) at the same position increased its nucleophilicity (azide 13, Table 2, entries [10][11][12]. Hence, the reaction yield for the re-gioisomer a, as an isolated compound, was excellent in all cases. In addition, only a single regioisomer was obtained for all synthesized compounds, proving the high regioselectivity of the iron(III)-catalyzed reaction.  As shown in Table 2, the N1-propargyl nucleobases 7-9 showed a high reactivity towards the iron(III)-catalyzed 1,3-dipolar cycloaddition, and high reaction yields were obtained. Conversely, the reactivity of the azides depended on the nature of the aryl or alkyl group. In fact, benzyl azide 11 provided products 17a-19a (Table 2, entries 4-6) with slightly lower yields than phenyl azide 10 ( Table 2, entries 1-3).
Aromatic azides 10, 12, and 13 exhibited different reactivities depending on the substitution of the aromatic ring. In fact, the presence of a strong electron-withdrawing group (-NO2) at the para position reduced the nucleophilicity of the azide 12 (Table 2, entries 7-9), while the presence of the electron-donor group (CH3O − ) at the same position increased its nucleophilicity (azide 13, Table 2, entries 10-12). Hence, the reaction yield for the regioisomer a, as an isolated compound, was excellent in all cases. In addition, only a single regioisomer was obtained for all synthesized compounds, proving the high regioselectivity of the iron(III)-catalyzed reaction.
Finally, we proposed a possible catalytic reaction mechanism for the reaction between 7 and 10 as illustrated in Scheme 4. The first reaction step is the coordination of iron(III) chloride to the alkyne group of N-propargylthymine 7, generating intermediate a with increased electrophilicity. Activated dipolarophile a reacts with phenyl azide 10 through a concerted 1,3-dipolar cycloaddition reaction to give intermediate c. Subsequent heterocycle aromatization of c gives product 14 and allows catalyst turnover.  Finally, we proposed a possible catalytic reaction mechanism for the reaction between 7 and 10 as illustrated in Scheme 4. The first reaction step is the coordination of iron(III) chloride to the alkyne group of N-propargylthymine 7, generating intermediate a with increased electrophilicity. Activated dipolarophile a reacts with phenyl azide 10 through a concerted 1,3-dipolar cycloaddition reaction to give intermediate c. Subsequent heterocycle aromatization of c gives product 14 and allows catalyst turnover.

Molecular Docking
Finally, to evaluate the potential biological activity of our products, molecular docking simulations were carried out using GOLD (CCDC Discovery) [41] and the ChemScore scoring function [42]. An in silico study was performed over 16 targets for the full set of the compounds synthesized in this work (see Supplementary Materials for further details). Compounds 20a and 21a were the best ranked ones (Figure 2), matching or exceeding in some cases the score of the co-crystallized ligand in the original structure, suggesting that these compounds might be able to bind the selected targets. Such targets have been at the center of great attention in the last years due to their possible implication in different pathologies. Human norovirus is one of the major causes of nonbacterial gastroenteritis in humans, and targeting the protruding P domain dimer (P-dimer) of a GII.10 Norovirus strain (Figure 2A) could be a successful strategy in drug discovery [43]. Eosinophil-derived neurotoxin (EDN) ( Figure 2B) is a member of the Ribonuclease A (RNase A) superfamily involved in inflammatory disorders and the immune response system [44]. Metallo βlactamases ( Figure 2C) are a family of enzymes employed by bacteria to hydrolyze β-lactam drugs as carbapenems, determining the resistance to antibiotics [45]. Hence, the discovery of new inhibitors capable of blocking such receptors could be of interest in combating bacterial infective diseases. As for the target reported in Figure 2D, SARS-CoV-2 NSP13 helicase was described by Newman et al. [46] in 2021 as a potential target for new antivirals due to its essential role in viral replication and its high sequence conservation.

Molecular Docking
Finally, to evaluate the potential biological activity of our products, molecular docking simulations were carried out using GOLD (CCDC Discovery) [41] and the ChemScore scoring function [42]. An in silico study was performed over 16 targets for the full set of the compounds synthesized in this work (see Supplementary Materials for further details). Compounds 20a and 21a were the best ranked ones (Figure 2), matching or exceeding in some cases the score of the co-crystallized ligand in the original structure, suggesting that these compounds might be able to bind the selected targets. Such targets have been at the center of great attention in the last years due to their possible implication in different pathologies. Human norovirus is one of the major causes of nonbacterial gastroenteritis in humans, and targeting the protruding P domain dimer (P-dimer) of a GII.10 Binding interactions mostly involve hydrogen bonds with the protein backbone or with charged residues, which are in general very much favored from the energetic point of view and capable of dictating the local structure [47].

Pharmacokinetics and ADMET Study
The drug-likeness of the newly synthesized hybrids was further evaluated. The analysis of absorption, distribution, metabolism, and excretion (ADME) was determined in silico by using the online database ADMETlab 2.0 [48]. As reported in Table 3, all hybrids expressed good ADME properties. They showed good water solubility and no violation of the Lipinsky rules of 5 [49]. Although the Caco-2 permeability was not excellent in some cases, effective intestinal absorption was found for all compounds, not only with the ADMETlab 2.0 software but also using the Brain Or IntestinaL EstimateD permeation method (BOILED-Egg) that evaluates the accessibility of compounds to the gastrointestinal (GI) tract and blood-brain barrier [50] (Figure 3). The boiled egg model revealed that all molecules had satisfactory GI absorption and no sufficient permeability across the blood-brain barrier (BBB), thus indicating good safety for the Central Nervous System (CNS). Moreover, all compounds had a good clearance rate (excretion rate) from the body. The encouraging results reported in Table 3 suggest that these compounds could be eligible as drug candidates, confirming the use of 20a and 21a as potential lead compounds for a new class of active pharmaceutical ingredients. Binding interactions mostly involve hydrogen bonds with the protein backbone or with charged residues, which are in general very much favored from the energetic point of view and capable of dictating the local structure [47]. Ligands and binding residues are shown as sticks. Non-polar hydrogens have been omitted for clarity. In the table, the docking scores of compounds 14a-25a and literature ligands were compared. The green color refers to scores higher than the literature ligands. The orange color refers to scores lower than the literature ligands. For those cases in which the docking score was higher than literature ligands the receptors are reported in bold.

General Procedure for Nucleobase Propargylation 7-9
In a three-necked round-bottomed flask, equipped with a bubble condenser and

General Procedure for Nucleobase Propargylation 7-9
In a three-necked round-bottomed flask, equipped with a bubble condenser and magnetic stir bar, the opportune nucleobase 1-3 (39.6 mmol, 1 eq) in dry hexamethyldisilazane (HMDS, 139 mmol, 3.5 eq) was suspended under nitrogen atmosphere. Subsequently, trimethylsilyl chloride (8.71 mmol, 0.22 eq) and (NH 4 ) 2 SO 4 (1.98 mmol, 0.05 eq) were added, and the mixture was stirred at 145 • C for 2 h. After completion, the solution was cooled to room temperature, and the HMDS was evaporated under vacuum. Then, the obtained silylated nucleobase 4-6 was dissolved in dry acetonitrile (150 mL) without any further purification. Propargyl bromide (39.6 mmol, 1 eq) was added dropwise at 80 • C for 30 min, and the reaction was stirred under reflux for 12 h. After cooling to room temperature, acetonitrile was removed under vacuum, and the crude was purified by silica flash chromatography (eluent mixture CHCl 3 /CH 3 OH 8:2) to give a solid product 7-9. For characterization data, see Supplementary Materials.

General Procedure for Nucleobase-Containing 1,5-Disubstituted 1,2,3-Triazoles 14a-25a
In a 50 mL two-necked round-bottomed flask equipped with a bubble condenser and magnetic stir bar, propargyl nucleobase 7-9 (1.52 mmol, 1 eq) was dissolved in DMF (8 mL). Subsequently, FeCl 3 (0.304 mmol, 0.2 eq) and opportune azide 10-13 (3.05 mmol, 2 eq) were added, and the mixture was stirred at 120 • C for 8 h. DMF was removed under vacuum by generating an azeotrope with toluene, and the obtained crude solid was purified on a flash silica gel column (eluent mixture: CHCl 3 /acetone/CH 3 OH 8:1:1 v/v/v) to obtain the desired solid product 14a-25a. The configuration of regioisomers was determined by spectroscopic data reported in Supplementary Materials. Furthermore, to calculate the regioisomeric ratio, the C6H proton on the nucleobases was chosen for both regioisomers (see Supplementary Materials).

Docking Studies
For each receptor, the docking cavity was centered on the binding site of the crystallographic ligand and allowed to extend in a spherical surrounding volume with a radius of 15 Å. In cases where a metal ion was present at the binding site, the docking cavity was centered on it, and metal parameters were set to maintain the same coordination number as that in the crystallographic structure. In the absence of metals, the XYZ coordinates that defined the center of the cavity were obtained from the position of the co-crystalized ligand, choosing an atom that was reasonably at the center of the ligand. The number of genetic algorithm runs was set to 20 for each analyzed ligand. Protein structures were prepared using UCSF Chimera [51], by reverting selenomethionine to methionine, eliminating alternate locations of side chains, adding hydrogen atoms, assigning appropriate protein atom types, and removing the co-crystalized ligand and solvent molecules. Crystallographic ligands were docked after adding hydrogen atoms with UCSF Chimera and without optimizing their geometries. Conversely, the geometries of the screened ligands 14a-25a were optimized quantum mechanically. Geometry optimizations and frequency calculations for stationary point characterization were carried out with Gaussian16 [52] using the M06-2X hybrid functional [53], the 6-31G(d,p) basis set, and ultrafine integration grids. Bulk solvent effects in water were considered implicitly through the IEF-PCM polarizable continuum model [54]. As for the potential receptors, 26 targets were initially selected from the Protein Data Bank ( Figure S1). The main selection criterion was the presence of a triazole (i.e., 1,2,3and 1,2,4-triazoles) scaffold or structurally similar heterocycles (i.e., imidazoles, thiazoles) in the crystallographic structure of the ligand-receptor complex. Docking simulations were performed, keeping the coordinates of the protein fixed while allowing flexibilization of the ligands around their rotatable bonds.

Prediction of Pharmacokinetic Properties
The absorption, distribution, metabolism, and excretion (ADME) analysis and assessment of Lipinski parameters were carried out using the free software "ADMETlab 2.0" supported by the Xiangya School of Pharmaceutical Sciences, Central South University ADMETlab 2.0 (2021) https://admetmesh.scbdd.com (accessed on 12 September 2022) [48]. The main properties considered were physical and chemical properties such as molecular weight, oil/water partition coefficient, and water solubility. Medicinal chemistry parameters were evaluated: the synthetic accessibility score (designed to estimate the ease of synthesis of drug-like molecules) and the violation of Lipinsky's rules. Absorption was studied mainly considering the permeability after oral administration (Caco-2 cell permeability) and human intestinal absorption (HIA). Among the distribution parameters, Plasma Protein Binding and blood brain-barrier (BBB) penetration were selected. CYP2D6 and CYP3A4 were selected as key enzymes in the metabolism of the compounds. Furthermore, clearance was evaluated as an excretion parameter. Finally, the Brain Or IntestinaL EstimateD permeation method (BOILED-Egg) was used to clarify the main distribution tissue for compounds 14a-25a [50].

Conclusions
We have developed a new synthetic approach towards the highly regioselective synthesis of 1,5-disubstituted 1,2,3-triazole derivatives containing pyrimidine nucleobases. N1-propargyl nucleobases were first prepared in excellent yields through an alternative route that led only to the N1-monoalkylated derivatives. These alkyne compounds were characterized and used in a 1,3-dipolar cycloaddition reaction to obtain the 1,2,3 triazole heterocycles. The 1,5-disubstituted-1,2,3-triazol moiety was obtained in a regioselective way by using FeCl 3 as an inexpensive and non-toxic catalyst. It was demonstrated that this procedure worked well for different aryl and benzyl azides. Although some differences were found with aryl azides bearing strong electron-withdrawing and electron-donating groups, reaction yields ranged from good to excellent in all cases. Finally, compounds 20a and 21a exhibited the best docking scores against four of the selected protein targets, suggesting potential biologic activity for these scaffolds. In addition, ADMET analysis was also performed using an in silico online database, which predicted that all hybrids might have qualified pharmacokinetic parameters and Lipinski's rule of five properties to claim their intestinal absorption and good clearance rate.