Novel 1,2,3-Triazole-Based Benzothiazole Derivatives: Efficient Synthesis, DFT, Molecular Docking, and ADMET Studies

A new series of 1,2,3-triazole derivatives 5a–f based on benzothiazole were synthesized by the 1,3-dipolar cycloaddition reaction of S-propargyl mercaptobenzothiazole and α-halo ester/amide in moderate to good yields (47–75%). The structure of all products was characterized by 1H NMR, 13C NMR, and CHN elemental data. This protocol is easy and green and proceeds under mild and green reaction conditions with available starting materials. The structural and electronic analysis and 1H and 13C chemical shifts of the characterized structure of 5e were also calculated by applying the B3LYP/6-31 + G(d, p) level of density functional theory (DFT) method. In the final section, all the synthesized compounds were evaluated for their anti-inflammatory activity by biochemical COX-2 inhibition, antifungal inhibition with CYP51, anti-tuberculosis target protein ENR, DPRE1, pks13, and Thymidylate kinase by molecular docking studies. The ADMET analysis of the molecules 5a–f revealed that 5d and 5a are the most-promising drug-like molecules out of the six synthesized molecules.

Due to the importance and application of 1,2,3-triazoles bearing the benzothiazole group [27,28], the synthesis of these derivatives is essential for organic chemists and biologists, and efforts are underway to design drugs based on these compounds with high activity and lower toxicity. As a result of our group's experience synthesizing 1,4disubstituted 1,2,3-triazole derivatives [29,30,31,32,33,34], the new 1,2,3-triazoles were prepared using benzothiazole groups 5a-f via a 1,3-dipolar cycloaddition of Spropargylated mercaptobenzothiazole 3, sodium azide, and α-halo ester/amide 4 in the presence of CuSO4/sodium ascorbate in a mixture of H2O/t-BuOH (Scheme 1). The density functional theory (DFT) calculations, with the B3LYP/6-31 + G(d, p) method, were applied to investigate the physicochemical properties of the characterized structure of 5e, including structural, electronic, and spectral analysis. Moreover, the study concludes by examining each of the synthesized compounds to specific targets using molecular docking to predict the possible biological activity of the proposed molecules. As mentioned above, triazoles have anti-inflammatory, anti-tuberculosis, and anti-fungal properties. For this purpose, the COX-2 protein was chosen for the anti-inflammatory target. COX-2 has been an important target for pain reduction drugs through the years and even shows some evidence relating this protein to cancer [35,36,37]. Numerous studies have been conducted on this protein with triazole structures in silico, in vitro, and in vivo that have shown promising results [38,39]. For anti-tuberculosis activity, four different targets were selected. DPRE1 has been one of the potential targets for resistant tuberculosis. Antituberculosis target protein ENR was another protein target on our list. The thymidylate kinase of Mycobacterium tuberculosis and Pks13 are the other anti-tuberculosis targets we selected for this research due to their critical role in this disease. For anti-fungal activity, we analyzed our synthesized molecule's activity against Saccharomyces cerevisiae CYP51 [40,41,42,43,44,45,46,47,48,49,50,51]. Absorption, distribution, metabolism, excretion, and toxicities (ADMET) prediction is one of the essential steps for drug discovery and development. After the docking study, we analyzed the molecules 5a-f for ADMET properties to analyze and predict these molecules against the control molecules for each target to better understand the chances of introducing them to the next phase of drug discovery and development. Our molecular docking experiments for the molecule 5e, which has a chiral center, were conducted to check how the enantiomers affect its potential activity and interactions with proteins. Due to the limitations of ADMET servers, which predict properties based on the 2D structure of molecules, we could not generate separate properties for each of the two enantiomers.

Results and Discussion
2.1. Synthesis and Spectroscopic Characterization of 1,2,3-Triazoles 5a-f 2-Mercaptobenzothiazole 2 was obtained according to the literature [52] and then was propargylated by propargyl bromide at 45 • C in dioxane in the presence of Et 3 N (Scheme 1). The azide intermediates were obtained by the in situ reaction of the corresponding α-halo ester/amide 4 with sodium azide. The desired 1,2,3-triazoles 5a-f were prepared by the azide-alkyne [3 + 2] cycloaddition reaction at room temperature. To optimize the reaction conditions, the solvent, catalyst, and temperature effect were investigated on the model reaction of 3 with 4a in the presence of sodium azide (Table 1). It should be mentioned that increasing the temperature had no good effect on the yield of the reaction (Entries 1-3 and 5, Table 1). With the inclusion of a low amount of copper catalyst (30 or 10 mol%), the reaction took longer, as well as resulted in a low yield of the product. CuSO 4 (50 mol%) as a catalyst, sodium ascorbate as a reductant agent (50 mol%), and H 2 O/t-BuOH (1:1) as solvent at room temperature were applied as the best reaction conditions (Entry 4, Table 1).

Synthesis and Spectroscopic Characterization of 1,2,3-Triazoles 5a-f
2-Mercaptobenzothiazole 2 was obtained according to the literature [52] and then was propargylated by propargyl bromide at 45 °C in dioxane in the presence of Et3N (Scheme 1). The azide intermediates were obtained by the in situ reaction of the corresponding α-halo ester/amide 4 with sodium azide. The desired 1,2,3-triazoles 5a-f were prepared by the azide-alkyne [3 + 2] cycloaddition reaction at room temperature. To optimize the reaction conditions, the solvent, catalyst, and temperature effect were investigated on the model reaction of 3 with 4a in the presence of sodium azide (Table 1). It should be mentioned that increasing the temperature had no good effect on the yield of the reaction (Entries 1-3 and 5, Table 1). With the inclusion of a low amount of copper catalyst (30 or 10 mol%), the reaction took longer, as well as resulted in a low yield of the product. CuSO4 (50 mol%) as a catalyst, sodium ascorbate as a reductant agent (50 mol%), and H2O/t-BuOH (1:1) as solvent at room temperature were applied as the best reaction conditions (Entry 4, Table 1). The generality of the protocol was then studied. As shown in Scheme 2, 1,2,3-triazole derivatives (5a-f) were synthesized in moderate to good yields (47-75%) by the reaction of a series of α-halo ester/amide 4 and S-propargyl mercaptobenzothiazole 3. The structure of triazoles 5a-f was established based on 1 H NMR, 13 C NMR spectral data, and CHN analysis. As a representative example, the 1 H NMR and 13 C NMR data of the product 5b are discussed. The 1 H NMR spectrum of 5b consisted of a singlet peak at δ = 3.72 ppm for the hydrogens of the methoxy group, a broad line at δ = 4.75 ppm for the methylene of the -CH2Sgroup, a broad line at δ = 5.27 ppm for -CH2CO-, eight aromatic hydrogens at δ = 6.91-8.04 ppm, a sharp characteristic signal at δ = 8.16 ppm for hydrogen on the triazole ring, and a broad line at δ = 10.34 ppm correlating with -CONH (amidic hydrogen). The amidic hydrogen (acceptor) (N-H...N) forms a strong intramolecular hydrogen bond with the nitrogen atom of the triazole ring (donor), resulting in a deshielded proton for the amide moiety in 1 H NMR. The 1 H-decoupled 13 CNMR spectrum of 5b showed 17 distinct resonances, a resonance at δ = 166.30 ppm for the carbon of carbonyl group, 13 distinct resonances for aromatic carbons between δ = 114.44-164.02 ppm, two peaks at δ = 55.60 and 52.61 ppm for the carbons of the methylene groups, and a resonance at δ = 27.85 ppm for the carbon of the methoxy group, which is in agreement with the proposed structure. The generality of the protocol was then studied. As shown in Scheme 2, 1,2,3-triazole derivatives (5a-f) were synthesized in moderate to good yields (47-75%) by the reaction of a series of α-halo ester/amide 4 and S-propargyl mercaptobenzothiazole 3. The structure of triazoles 5a-f was established based on 1 H NMR, 13 C NMR spectral data, and CHN analysis. As a representative example, the 1 H NMR and 13 C NMR data of the product 5b are discussed. The 1 H NMR spectrum of 5b consisted of a singlet peak at δ = 3.72 ppm for the hydrogens of the methoxy group, a broad line at δ = 4.75 ppm for the methylene of the -CH 2 S-group, a broad line at δ = 5.27 ppm for -CH 2 CO-, eight aromatic hydrogens at δ = 6.91-8.04 ppm, a sharp characteristic signal at δ = 8.16 ppm for hydrogen on the triazole ring, and a broad line at δ = 10.34 ppm correlating with -CONH (amidic hydrogen). The amidic hydrogen (acceptor) (N-H . . . N) forms a strong intramolecular hydrogen bond with the nitrogen atom of the triazole ring (donor), resulting in a deshielded proton for the amide moiety in 1 H NMR. The 1 H-decoupled 13 CNMR spectrum of 5b showed 17 distinct resonances, a resonance at δ = 166.30 ppm for the carbon of carbonyl group, 13 distinct resonances for aromatic carbons between δ = 114.44-164.02 ppm, two peaks at δ = 55.60 and 52.61 ppm for the carbons of the methylene groups, and a resonance at δ = 27.85 ppm for the carbon of the methoxy group, which is in agreement with the proposed structure.

Theoretical Studies on the 1,2,3-Triazole 5e
In continuation, the physicochemical properties of 5e, such as structural, electronic, and spectral data, were investigated by the DFT-B3LYP computational method in conjunction with the 6-31 + G(d, p) basis set.

Benchmarking of the Computational Method
DFT calculations have successfully provided the theoretical background of popular qualitative chemical concepts [53]. Furthermore, organic chemists use the DFT-B3LYP method for various chemical calculations [54][55][56][57]. This method has been utilized in the present study for performing calculations. The experimental 1 H NMR and 13 C NMR data for 5e could pave the way to detect the accurate computational method.
The experimental 1 H NMR spectrum of 5e consisted of a triplet at δ = 1.10 for methyl, a doublet for methyl at 1.72, a multiplet at 4.41 for methylene (-OCH2-), and a singlet at 4.71 for methylene (-SCH2−), a quartet at 5.62 for -NCH-(1H), four protons for phenyl rings at δ = 7.35-7.99 ppm, and a characteristic signal at δ = 8.24 ppm for hydrogen on the triazole ring. The 1 H-decoupled 13 C NMR spectrum of 5e showed 15 resonances, which is in agreement with the proposed structure, while the carbonyl appearing at δ = 169.21 ppm, nine distinct resonances for the aromatic carbons of phenyl and triazole rings between δ = 121.30-165.74 ppm, and five resonances at δ = 13.80-61.66 ppm for aliphatic carbons. The simulated NMR spectra for 5e were obtained at the B3LYP/6-31 + G(d, p) level with the aid of the Gauss View software [58]. Table 2 summarizes the calculated and experimental 1 H and 13 C chemical shifts for 5e. The 1 H and 13 C chemical shifts are reported in parts per million (ppm) relative to tetramethylsilane as the reference, and its corresponding shifts were calculated at the same theoretical level. For the 6-31 + G(d, p)-optimized geometry of 5e in DMSO, the NMR chemical shifts were also calculated using the same basis set and solvent as the experimental NMR spectra.

Theoretical Studies on the 1,2,3-Triazole 5e
In continuation, the physicochemical properties of 5e, such as structural, electronic, and spectral data, were investigated by the DFT-B3LYP computational method in conjunction with the 6-31 + G(d, p) basis set.

Benchmarking of the Computational Method
DFT calculations have successfully provided the theoretical background of popular qualitative chemical concepts [53]. Furthermore, organic chemists use the DFT-B3LYP method for various chemical calculations [54][55][56][57]. This method has been utilized in the present study for performing calculations. The experimental 1 H NMR and 13 C NMR data for 5e could pave the way to detect the accurate computational method.
The experimental 1 H NMR spectrum of 5e consisted of a triplet at δ = 1.10 for methyl, a doublet for methyl at 1.72, a multiplet at 4.41 for methylene (-OCH 2 -), and a singlet at 4.71 for methylene (-SCH 2 −), a quartet at 5.62 for -NCH-(1H), four protons for phenyl rings at δ = 7.35-7.99 ppm, and a characteristic signal at δ = 8.24 ppm for hydrogen on the triazole ring. The 1 H-decoupled 13 C NMR spectrum of 5e showed 15 resonances, which is in agreement with the proposed structure, while the carbonyl appearing at δ = 169.21 ppm, nine distinct resonances for the aromatic carbons of phenyl and triazole rings between δ = 121.30-165.74 ppm, and five resonances at δ = 13.80-61.66 ppm for aliphatic carbons. The simulated NMR spectra for 5e were obtained at the B3LYP/6-31 + G(d, p) level with the aid of the Gauss View software [58]. Table 2 summarizes the calculated and experimental 1 H and 13 C chemical shifts for 5e. The 1 H and 13 C chemical shifts are reported in parts per million (ppm) relative to tetramethylsilane as the reference, and its corresponding shifts were calculated at the same theoretical level. For the 6-31 + G(d, p)-optimized geometry of 5e in DMSO, the NMR chemical shifts were also calculated using the same basis set and solvent as the experimental NMR spectra.  The calculated 1 H and 13 C chemical shifts are in excellent agreement with the experimental data; regression was obtained R 2 = 0.996 for 1 H NMR and R 2 = 0.995 for 13 C NMR data ( Figure 2). 1 H NMR comparison showed a better correlation rather than 13 C NMR. Based on the results, it is possible to continue the calculations and obtain the physicochemical properties of 5e using the B3LYP/6-31 + G(d, p) method.  The calculated 1 H and 13 C chemical shifts are in excellent agreement with the experimental data; regression was obtained R 2 = 0.996 for 1 H NMR and R 2 = 0.995 for 13 C NMR data ( Figure 2). 1 H NMR comparison showed a better correlation rather than 13 C NMR. Based on the results, it is possible to continue the calculations and obtain the physicochemical properties of 5e using the B3LYP/6-31 + G(d, p) method. a For the numbering of atoms, refer to Figure 1. The calculated 1 H and 13 C chemical shifts are in excellent agreement with the experimental data; regression was obtained R 2 = 0.996 for 1 H NMR and R 2 = 0.995 for 13 C NMR data (Figure 2). 1 H NMR comparison showed a better correlation rather than 13 C NMR. Based on the results, it is possible to continue the calculations and obtain the physicochemical properties of 5e using the B3LYP/6-31 + G(d, p) method.   Table 3. The calculated structural parameters of 5e by the B3LYP/6-31 + G(d, p) method.

Geometrical Parameters a Theoretical Value
Bond lengths (Å) C 1 -C 2 1.397 C 2 -S 12 1.757 N 13 -C 11 1.298 C 11 -S 14 1.761 S 14 -C 15 1.856 N 21 -N 22 1.306 C 30 -O 39 1.215 C 38 -C 31 1.463 C 20 -C 18 1.382 Bond angles ( • ) C 1 -C 2 -C 3 121.7 H 10 -C 6 -C 5 119.7 C 11 -S 14  As mentioned above, there is strong intramolecular hydrogen bonding (IMHB) between the amidic hydrogen and nitrogen atom of the triazole ring (N-H . . . N) in the chemical structures of 5a-d. This interaction was confirmed from the 1 H NMR spectra of 5a-d, resulting in a deshielded proton for the amide group. As can be seen from the optimized structures of 5a (Figure 3), there is a very strong IMHB with the hydrogen atom of the amide group as the acceptor center and the nitrogen atom of the triazole ring as the donor site. By analyzing the natural bond orbital (NBO) of 5a, we can gain valuable insight into the intermolecular/intramolecular interactions, such as hydrogen bonding, numerically according to two-order perturbation energy (E (2) ) [59]. The LP (1) N 22 → BD*(1) N 28 -H 41 as a remarkable interaction with energy of 4.85 kcal mol −1 in 5a is attributed to the IMHB interaction. The structures of 5b-d exhibit also an analogous trend.  Table 4 shows the electronic data calculated by the B3LYP/6-31 + G(d, p) meth 5e in DMSO. The dipole moment in DMSO for 5e is 5.81 Debye, indicating the polarity. The presence of polar groups, such as the amide group and unsymme structure, are the reasons for this high dipolar moment data.   Table 4 shows the electronic data calculated by the B3LYP/6-31 + G(d, p) method of 5e in DMSO. The dipole moment in DMSO for 5e is 5.81 Debye, indicating the high polarity. The presence of polar groups, such as the amide group and unsymmetrical structure, are the reasons for this high dipolar moment data. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are the frontier orbitals that reveal the ability to donate or accept an electron, respectively. It is worth mentioning that removing electrons from HOMO is easier, and this orbital can donate electrons. As shown in Figure 4, the HOMO and LUMO orbitals are localized mainly on the benzothiazole ring. A useful parameter for determining the organic compound's reactivity is the gap between the HOMO and LUMO energy levels (E g ). The difference between E HOMO and E LUMO (E g ) is 4.98 eV (Table 4).  Table 4 shows the electronic data calculated by the B3LYP/6-31 + G(d, p) method of 5e in DMSO. The dipole moment in DMSO for 5e is 5.81 Debye, indicating the high polarity. The presence of polar groups, such as the amide group and unsymmetrical structure, are the reasons for this high dipolar moment data. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are the frontier orbitals that reveal the ability to donate or accept an electron, respectively. It is worth mentioning that removing electrons from HOMO is easier, and this orbital can donate electrons. As shown in Figure 4, the HOMO and LUMO orbitals are localized mainly on the benzothiazole ring. A useful parameter for determining the organic compound's reactivity is the gap between the HOMO and LUMO energy levels (Eg). The difference between EHOMO and ELUMO (Eg) is 4.98 eV (Table 4). The global electrophilicity index (GEI) is a criterion of the molecule's ability to accept electrons, which has been used to determine some adequate information on the structure, reactivity, dynamics, aromaticity, Lewis acidity, and toxicity [60,61]. Parr defined the GEI (ω) as [62]: The global electrophilicity index (GEI) is a criterion of the molecule's ability to accept electrons, which has been used to determine some adequate information on the structure, reactivity, dynamics, aromaticity, Lewis acidity, and toxicity [60,61]. Parr defined the GEI (ω) as [62]:

Electronic Characterization of 5e
where µ and η are the chemical potential and hardness, respectively, given by Equation (2); 5e has a GEI of 2.94 eV. The GEI can be categorized based on Domingo's scale [63]. Electrophiles are divided into three groups based on the GEI value: strong, mild, and weak. When the electrophilic index ω is more than 1.5 eV, the electrophile is strong; in the range of 0.8 to 1.5 eV, it is mild, and when it is less than 0.8 eV, it is a weak electrophile. Based on Domingo's scale, the compound 5e is a strong electrophile and can easily interact with a typical nucleophile. The NBO charge distribution for the structure of 5e is shown in Figure 5. According to the NBO analysis, there is a remarkable excess of positive charge (1.466 e) on the carbon atom of the carbonyl group. The study also indicates that the carbon atom of the imine unit (N = C) in the benzothiazole ring, methylene (-SCH 2 -), and -NCH-bear a positive charge of 1.106, 0.338, and 0.421e, respectively. NBO analysis and the electrophilicity index of 5e show that a nucleophile can attack these positive carbon centers.
Possible nucleophilic centers in 5e, based on NBO analysis and the electrophilicity index, are shown in Figure 5.
range of 0.8 to 1.5 eV, it is mild, and when it is less than 0.8 eV, it is a weak electrophile. Based on Domingo's scale, the compound 5e is a strong electrophile and can easily interact with a typical nucleophile. The NBO charge distribution for the structure of 5e is shown in Figure 5. According to the NBO analysis, there is a remarkable excess of positive charge (1.466 e) on the carbon atom of the carbonyl group. The study also indicates that the carbon atom of the imine unit (N = C) in the benzothiazole ring, methylene (-SCH2-), and -NCH-bear a positive charge of 1.106, 0.338, and 0.421e, respectively. NBO analysis and the electrophilicity index of 5e show that a nucleophile can attack these positive carbon centers. Possible nucleophilic centers in 5e, based on NBO analysis and the electrophilicity index, are shown in Figure 5.

Molecular Docking and ADMET Studies
Molecular docking results show that molecule 5d has the highest binding energy for 1G3U (−9.7 kcal mol −1 ), 1QSG (−9.1 kcal mol −1 ), and 5F19 (−9.6 kcal mol −1 ), even more than the control molecules. Based on this information, 5d is likely to be a promising inhibitor of the proposed target. Molecule 5b has the same binding energy as 5d for 1QSG, giving the idea of the potential anti-tuberculosis activity of this molecule. In 4ZE2, the binding energy of the control molecule is similar to 5a, equal to −9.7 kcal mol −1 . The synthesized molecules 5a-f show higher binding energy and more minor interactions with two other proteins (4P8N, 5V3X), indicating that these molecules have a lower anti-tuberculosis activity than their reference molecules. In 4P8N, the control molecule BTZ043 shows a binding energy of −10.6 kcal mol −1 . In addition, 5d has the lowest binding energy of all proposed molecules, with a binding energy of −9.6 kcal mol −1 . As shown by the binding energy of 5V3X, the reference molecule has a binding energy of −10.1 kcal mol −1 , which indicates a strong interaction between this molecule and the protein target. Given that the six synthesized molecules do not bypass the reference molecule and have lower inhibitory activity against the target, 5d has the best binding energy of −9.7 kcal mol −1 . A summary of the docking results can be found in Table 5.

Molecular Docking and ADMET Studies
Molecular docking results show that molecule 5d has the highest binding energy for 1G3U (−9.7 kcal mol −1 ), 1QSG (−9.1 kcal mol −1 ), and 5F19 (−9.6 kcal mol −1 ), even more than the control molecules. Based on this information, 5d is likely to be a promising inhibitor of the proposed target. Molecule 5b has the same binding energy as 5d for 1QSG, giving the idea of the potential anti-tuberculosis activity of this molecule. In 4ZE2, the binding energy of the control molecule is similar to 5a, equal to −9.7 kcal mol −1 . The synthesized molecules 5a-f show higher binding energy and more minor interactions with two other proteins (4P8N, 5V3X), indicating that these molecules have a lower antituberculosis activity than their reference molecules. In 4P8N, the control molecule BTZ043 shows a binding energy of −10.6 kcal mol −1 . In addition, 5d has the lowest binding energy of all proposed molecules, with a binding energy of −9.6 kcal mol −1 . As shown by the binding energy of 5V3X, the reference molecule has a binding energy of −10.1 kcal mol −1 , which indicates a strong interaction between this molecule and the protein target. Given that the six synthesized molecules do not bypass the reference molecule and have lower inhibitory activity against the target, 5d has the best binding energy of −9.7 kcal mol −1 . A summary of the docking results can be found in Table 5. Table 5. Binding energy data for the triazoles 5a-f and control molecules (in kcal mol −1 ). All compounds 5a-f can bind to the active site of the protein targets. The synthesized compounds can form hydrogen bonds (H-bond) with the amino acids in the active site of the protein targets. The number of H-bonds of 5a-f with the related protein targets is collected in Table 6. H-bond interactions were found between nitrogen atoms on the triazole ring of the compounds 5a-f and amino acids in the active site of the protein targets, It is interesting to note that the compound 5e has two enantiomers (5e-R and 5e-S). Molecular docking studies of these enantiomers can be interesting. Except for the protein 1QSG, the binding energy of 5e-S with other proteins is higher than 5e-R (Table 6). For example, we considered only the interaction of 5e-R and 5e-S with 1G3U and 4P8N. For 1G3U, the 5e-R interacts with ARG-160 through a hydrogen bond, one carbon-hydrogen interaction with PRO-37, six van der Waals interactions with SER-99, TYR-103, ASP-9, MG-300, ARG-14, and ALA-35, and several other interactions such as P-cation, Pi-anion, Pi-sulfur, Pi-Pi stacked, and Pi-alkyl. 5e-S shows more interactions, specifically four hydrogen bonds with LYS-13, ARG-95, TYR-39, and ARG-160. Moreover, it also has a carbon-hydrogen bond with ASP-9. The 5e-S molecule also has eight van der Waals interactions with ARG-153, ARG-14, ASP-94, ARG-74, SER-99, ASN-100, LEU-52, and ASP-162. Like 5e-R, 5e-S also interacts with the active site, interactions such as Pi-cation, P-anion, Pi-sulfur, Pi-Pi stacked, Pi-Pi T-shaped, and Pi-alkyl. These results could It is interesting to note that the compound 5e has two enantiomers (5e-R and 5e-S). Molecular docking studies of these enantiomers can be interesting. Except for the protein 1QSG, the binding energy of 5e-S with other proteins is higher than 5e-R (Table 6). For example, we considered only the interaction of 5e-R and 5e-S with 1G3U and 4P8N. For 1G3U, the 5e-R interacts with ARG-160 through a hydrogen bond, one carbon-hydrogen interaction with PRO-37, six van der Waals interactions with SER-99, TYR-103, ASP-9, MG-300, ARG-14, and ALA-35, and several other interactions such as P-cation, Pi-anion, Pi-sulfur, Pi-Pi stacked, and Pi-alkyl. 5e-S shows more interactions, specifically four hydrogen bonds with LYS-13, ARG-95, TYR-39, and ARG-160. Moreover, it also has a carbon-hydrogen bond with ASP-9. The 5e-S molecule also has eight van der Waals interactions with ARG-153, ARG-14, ASP-94, ARG-74, SER-99, ASN-100, LEU-52, and ASP-162. Like 5e-R, 5e-S also interacts with the active site, interactions such as Pi-cation, P-anion, Pi-sulfur, Pi-Pi stacked, Pi-Pi T-shaped, and Pi-alkyl. These results could indicate the higher binding affinity of 5e-S in comparison to 5e-R (Figure 7). Based on the interaction results for these enantiomers for 4P8N, 5e-R has two hydrogen bonds with The Lipinski, Pfizer, and GSK rules are the focus of the ADMET study. While the molecules 5a, 5b, and 5d showed the best docking results, the ADMET analysis reveals that all three mentioned molecules passed two out of the three rules above. The control molecule for 4ZE2, Posaconazole, passes only the Pfizer rule, suggesting that the molecules 5a, 5b, and 5d have a better chance than others, mainly molecule 5d, owing to its robust docking against anti-inflammatory and anti-tuberculosis targets. The ADMET The Lipinski, Pfizer, and GSK rules are the focus of the ADMET study. While the molecules 5a, 5b, and 5d showed the best docking results, the ADMET analysis reveals that all three mentioned molecules passed two out of the three rules above. The control molecule for 4ZE2, Posaconazole, passes only the Pfizer rule, suggesting that the molecules 5a, 5b, and 5d have a better chance than others, mainly molecule 5d, owing to its robust docking against anti-inflammatory and anti-tuberculosis targets. The ADMET results are summarized in Table 6.

Compound ID Protein PDB ID
Lipophilicity (permeability of a drug to reach the target tissue in the body) and solubility are two important molecular properties for the absorption of a drug. The calculated solubility in water (Log S) of the synthesized compounds 5a-f followed the order: 5f > 5e > 5a > 5b > 5d > 5c. The estimated lipophilicity (Log P) also followed the order: 5c > 5d > 5b > 5a > 5e > 5f. After analyzing the chemical structures of 5a-f (Scheme 2), it is easy to understand that, as the number of hydrocarbon chains increases, the solubility in water decreases and lipophilicity increases.

General Information
The chemicals and solvents were purchased from Sigma-Aldrich (St. Louis, MO. USA), Fluka (Neu-Ulm, Germany), and Merck (Darmstadt, Germany). The structure of the synthesized triazole derivatives 5a-f was confirmed by 1 H NMR, 13 C NMR, and CHN elemental data. A Bruker (DRX-500 Avance) NMR was used to record the 1 H and 13 C NMR spectra in CDCl 3 and DMSO at room temperature (The Supplementary Material). The multiplicity of the proton signals is abbreviated with s for singlet, d for doublet, m for multiplet, and dd for the doublet of doublet peaks. Quantitative microanalyses (CHN) of the products 5a-f were carried out on a Thermo Finnigan Flash-1112EA microanalyzer.

Computational Details
Geometric optimization of the characterized structure 5e was performed at the DFT/B3LYP level using the 6-31 + G(d, p) basis set utilizing the Gaussian 09 software [64]. The vibrational frequencies were calculated simultaneously to ensure that each stationary point was a real minimum. NMR data calculations were performed using the gaugeindependent atomic orbital (GIAO) method [65]. Relative chemical shifts were calculated using the corresponding tetramethylsilane for 1 H NMR and 13 C NMR shielding, calculated at the same level of theory as the reference. Further, the solvent effects were taken into account by using the conductor-like polarizable continuum model (CPCM) [66] through DFT calculations performed at the B3LYP/6-31 + G(d, p) level of theory in DMSO as the solvent.

Molecular Docking and ADMET Details
Molecular docking is a crucial component of the in silico study of drug discovery and development. We performed in silico docking studies with Pyrx to demonstrate the synthesized compounds' anti-tuberculosis, anti-fungal, and anti-inflammatory properties. Pyrx was used to determine the interactions between the studied compounds and the specific protein targets [67]. All the molecules' smiles were extracted by Cheminfo [68]. Afterward, all the structures were minimized using chimera and saved as PDB files. The next step involved unifying the six PDB structures using OpenBabel and converting them to a single SDF file [69]. RSCB PDB was used to extract all protein structures. Chimera was applied in all protein structures to remove solvents, water molecules, unnecessary ions, and ligands [70]. We selected a single control molecule per protein target for molecular docking to compare our synthesized molecules 5a-f with the reference molecules. All the procedures to prepare the reference molecules for docking were the same as for 5a-f, as shown above. Each target protein's PDB ID and the control molecules are listed below in Table 7. Pyrx was then used to dock the preprocessed data. For the ADMET approach, we used the ADMETlab2 web server to analyze and predict the synthesized molecules' properties [71]. In this section, all synthesized molecules 5a-f and control molecules are converted into their smile strings and fed to ADMETlab2. To a solution of 2-aminothiophenol 1 (1 mmol) in 30 mL of KOH solution (10%), CS 2 (6 mmol) was added and refluxed for three hours. The mixture was cooled, then 25 mL water and 25 mL ethanol were added and heated until 40 • C. At this temperature, 4.5 mL acetic acid was added dropwise to the reaction mixture to adjust the pH = 6. The mixture was cooled to room temperature, filtrated, and dried overnight to obtain the desired 2-mercaptobenzothiazole 2 in 90% yield.

General Synthesis Procedure for S-Propargylated Mercaptobenzothiazole 3
To a solution of 2-mercaptobenzothiazole 2 (5 mmol) in 10 mL dioxane, Et 3 N (5 mmol) was added and stirred at room temperature for 45 min. Propargyl bromide (5 mmol, 80% in toluene) was added dropwise to the reaction mixture and refluxed at 45 • C for one hour. After completion of the reaction, the mixture was cooled, and then, 20 mL of water was added. The reaction mixture was extracted by DCM (3 × 5 mL). The organic layer was dried over Na 2 SO 4 , filtrated, and evaporated to give the desired S-propargyl mercaptobenzothiazole 3 in 95% yield.
3.6. General Synthesis Procedure for 1,2,3-Triazole Derivatives 5a-f To a solution of t-BuOH/H 2 O (1:1, 5 mL), α-halo ester/amide 4 (1 mmol) and NaN 3 (1.3 mmol) were added and stirred at room temperature for one hour, then CuSO 4 and sodium ascorbate (50% mol) and S-propargyl mercaptobenzothiazole 3 (1 mmol) were added. The reaction mixture was stirred at room temperature for 12 h. After completion of the reaction, as monitored by TLC, 10 mL of water was added, and the mixture was extracted with DCM (3 × 5 mL). The organic layer was dried over Na 2 SO 4 , filtrated, and evaporated under a vacuum. The desired products 5a-f were purified by column chromatography using ethyl acetate/hexane (1:1). The structure of the synthesized triazole derivatives 5a-f was confirmed by 1 H and 13 C NMR spectra data and CHN analysis.

Conclusions
We successfully reported synthesizing a series of new 1,2,3-triazoles based on benzothiazole by the 1,3-dipolar cycloaddition reaction of S-propargyl mercaptobenzothiazole and α-halo ester/amide in moderate to good yields. This protocol is easy and green and proceeds under mild reaction conditions with available starting materials. We also studied the full structural, electronic, and NMR data of the synthesized compound 5e using the B3LYP/6-31 + G(d, p) computational method. Molecular docking and ADMET studies indicate that the 5d scaffold could be the key to binding with the protein targets more than the other five synthesized molecules in terms of their anti-inflammatory and anti-tuberculosis activity. All synthesized molecules did not achieve a better docking result in terms of antifungal activity than the reference molecules; however, 5a has the same binding affinity for 4ZE2 as the control molecule, suggesting that this molecule may possess antifungal properties. Further in vitro analysis is needed to prove the in silico results.