Structural and Energetic Properties of Weak Noncovalent Interactions in Two Closely Related 3,6-Disubstituted-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazole Derivatives: In Vitro Cyclooxygenase Activity, Crystallography, and Computational Investigations

Two 3,6-disubstituted-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazole derivatives, namely, 3-(adamantan-1-yl)-6-(2-chloro-6-fluorophenyl)-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazole 1 and 6-(2-chloro-6-fluorophenyl)-3-phenyl-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazole 2, were prepared, and the detailed analysis of the weak intermolecular interactions responsible for the supramolecular self-assembly was performed using X-ray diffraction and theoretical tools. Analyses of Hirshfeld surface and 2D fingerprint plot demonstrated the effect of adamant-1-yl/phenyl moieties on intermolecular interactions in solid-state structures. The effect of these substituents on H···H/Cl/N contacts was more specific. The CLP-PIXEL and density functional theory methods provide information on the energetics of molecular dimers observed in these compounds. The crystal structure of compound 1 stabilizes with a variety of weak intermolecular interactions, including C–H···N, C–H···π, and C–H···Cl hydrogen bonds, a directional C–S···π chalcogen bond, and unconventional short F···C/N contacts. The crystal structure of compound 2 is stabilized by π-stacking interactions, C–H···N, C–H···π, and C–H···Cl hydrogen bonds, and highly directional attractive σ–hole interactions such as the C–Cl···N halogen bond and the C–S···N chalcogen bond. In addition, S(lp)···C(π) and short N···N contacts play a supportive role in the stabilization of certain molecular dimers. The final supramolecular architectures resulting from the combination of different intermolecular interactions are observed in both the crystal packing. The molecular electrostatic potential map reveals complementary electrostatic potentials of the interacting atoms. The quantum theory of atoms in molecules approach was used to delineate the nature and strength of different intermolecular interactions present in different dimers of compounds 1 and 2. The in vitro experiments suggest that both compounds showed selectivity against COX-2 targets rather than COX-1. Molecular docking analysis showed the binding pose of the compounds at the active sites of COX-1/2 enzymes.


RESULTS AND DISCUSSION
The purpose of this study is to investigate the effect of the adamantane and phenyl groups introduced in the core skeleton of the triazolothiadiazole ring on crystal packing and intermolecular interactions using various computational approaches. With the support of the in vitro data from the COX-1/2 assay, in silico molecular docking and free energy calculation, we describe the effect of the adamantane group on the observed biological activity.
2.1. Chemical Synthesis. The investigated compounds 1 and 2 were prepared using the following reaction sequences outlined in Scheme 1. The corresponding carboxylic acid hydrazide A was treated with carbon disulfide in an ethanolic potassium hydroxide solution to yield the corresponding potassium dithiocarbazate intermediates B, which were subsequently reacted with hydrazine to produce the corresponding 4-amino-5-substituted-4H-1,2,4-triazole-3-thiols C. 41,42 Compounds C were then reacted with 2-chloro-6fluorobenzoic acid and phosphorous oxychloride by heating under reflux for 2 hours to produce target compounds 1 and 2 in yields of 84 and 72%, respectively.

Description of Molecular and Crystal Structures.
For the two closely related triazolothiadiazoles, a lowtemperature (160 K) single-crystal structure determination has been performed to delineate the effect of adamantan-1-yl (compound 1) and phenyl (compound 2) moieties attached at position 3 (C9) carbon of the triazolothiadiazole ring on molecular conformation, crystal packing, and intermolecular interactions. Both compounds have a 2-chloro-6-fluorophenyl substituent at position 6 (C7) to study their role in crystal packing and intermolecular interactions. The X-ray analysis revealed that both compounds crystallize in the monoclinic system, and one molecule is present in the respective asymmetric unit. In compound 1, the 2-chloro-6-fluorophenyl ring was disordered over two orientations with the site occupation factors of 0.7764 (19) for the major disordered component. This major disordered component was used for all analyses reported in this work. Crystallographic data and refinement parameters for compounds 1 and 2 are presented in Table 1, and their thermal ellipsoid representations are shown in Figure 1.
The central triazolothiadiazole unit makes a right angle (86.6°) with the mean plane of the adamantane cage in compound 1, whereas the unsubstituted phenyl ring is nearly coplanar (3.4°) with the central triazolothiadiazole moiety in compound 2. We note that the orientation of the 2-chloro-6-Scheme 1. Synthetic Pathway for the Title Compounds 1 and 2 fluorophenyl ring is nearly perpendicular to the mean plane of the central unit in both structures (75.4°in 1 and 83.0°in 2). The Mogul geometry check suggests that both structures do not show unusual bond lengths or angles. However, it should be mentioned that the C7−S1 bond [1.761(1) Å in compound 1 and 1.757(1) Å in 2] exhibits bond lengthening compared to the C8−S1 bond [1.738(1) Å in 1 and 1.728(1) Å in 2]. A similar feature has been observed in related structures deposited in the CSD (version 5.42, November 2020). 43 This bond lengthening could be due to the formation of a σ− hole interaction (chalcogen bond). Furthermore, the CSD search detected 51 hits containing different substituents attached to atom C9 of the triazole ring and atom C7 of the thiadiazole ring. Five of these compounds have an adamantane cage attached to a thiadiazole ring (CSD refcodes: BOTYEP, UPAVIR, VUNLUM, VUNMAT, and WADHIU). 44−47 Unlike compound 1, which has an adamantane cage attached to a triazole ring, we identified eight structures containing monosubstituted (CSD refcodes: NITSIV, ILETOK, LEPQED, QEMMUR, and SAPRAD, UXIPOI), 36,48−52 disubstituted (CSD refcode: XOYWOZ), 53 and trisubstituted (csd refocde: GAJCUS) 38 halophenyl groups. The halogen atoms are either Cl or F or both present at different positions of the phenyl ring. Contrary to compounds 1 and 2, the substituted phenyl ring in these structures is coplanar with the central triazolothiadiazole ring. This planarity could be maintained due to the formation of intramolecular S···X (X = Cl or F) contact in all structures mentioned above except in 6-(4-chlorophenyl)-3-(thiophen-2-yl)- [1,2,4]triazolo [3,4-b]- [1,3,4]-thiadiazole (UXIPOI). 52 In this structure, Cl was substituted in the para position of the phenyl ring and does not participate in the S···Cl contact.
For a deeper understanding of why title compounds 1 and 2 do not exhibit either S···Cl or S···F intramolecular contact (or why disubstituted phenyl adopts non-coplanar with central unit?), as well as to uncover the energy barriers associated with it, we performed a relaxed potential surface scan around the C7−C1 bond in compound 2 at the B3LYP/6-31G+(d) level of theory ( Figure 2). The result suggests that the conformer "a" is the least energy conformer for the S1−C7−C1−C6 torsion angle of ±90°, which is very similar to that of the X-ray conformer (−97.08°). In the second least energy conformer "b", S and F atoms are displaced on the same side, whereas in the high energy conformer "c", S and Cl atoms are oriented on the same side. The energy difference between the conformers "a" and "b" is 2.8 kcal mol −1 , while the corresponding energy difference for the conformers "b" and "c" is only 1.3 kcal mol −1 (4.0 kcal mol −1 for the conformers "a" and "c"). This analysis demonstrates that the twisted orientation of the disubstituted phenyl ring has minimum energy on its potential energy surface, and a molecular conformation with S···F/Cl contact has a higher energy range from 2.8 to 4.0 kcal mol −1 concerning the minimum energy conformer.
As shown in Figure S1, the molecular conformation of compound 1 is stabilized by two intramolecular C−H···N interactions. Similarly, compound 2 has an intramolecular C− H···N interaction. The topological analysis confirms the existence of these interactions (Table S1). Some of the related structures are also stabilized by an intramolecular C−H···N interaction as observed in compounds 1 and 2 (CSD refcodes for these structures are given in Supporting Information: see Table S2). 36−38,44,48,49,52−63 Due to the existence of this intramolecular interaction, the central triazolothiadiazole moiety becomes coplanar with the substituent group attached to the triazole ring at atom C9 (corresponding atom in the related structures). In other structures, the substituent group on the triazole ring is slightly twisted to the mean plane of the central unit due to the lack of an intramolecular C−H···N interaction.

Hirshfeld Surface Analysis and 2D Fingerprint Plots.
Hirshfeld surface (HS) analysis has been widely used to qualitatively determine the intermolecular contacts within crystal structures. HS and decomposed 2D fingerprint (FP) plots reveal the effect of the phenyl and adamantane moieties on the intermolecular interactions. Figure 3a shows the HS over the d norm values for compounds 1 and 2. In compound 1, the intense red spots are shown for the intermolecular C−H··· N hydrogen bond and an attractive σ−hole interaction such as a chalcogen bond of type S···C(π) and an unorthodox bifurcated fluorine bond of type F···C/N contacts.
In compound 2, the S···N chalcogen bond shows intense and broad red surfaces, and a short N···N contact also shows an intense spot near the chalcogen interacting region ( Figure  3b). This feature suggests that the chalcogen bond along with a short N···N contact plays a significant role in the stabilization of the crystal structure. The relatively less intense spots were observed for C−H···N, C−H···Cl, and Cl···N interactions. Over the five-membered ring, we observed less intense red spots due to the π-stacking interaction. We also observed tiny red spots for the S(lp)···C(π) contact in the same π-stacking dimer, which might provide additional stabilization to this dimer. The 2D-FP plots show distinct distribution patterns for intercontacts in 1 and 2 (Figures 4, S2 and S3, Supporting Information). The intermolecular H···H interactions occupy most of the HS area toward the crystal packing of these structures. Due to the adamantane cage in compound 1, the H···H contacts increase (38.6%) and those in compound 2 decrease (19.7%) due to the phenyl ring. The next significant contribution comes from H···N contacts (14.3% in 1 and 13.1% in 2), and the contribution of these contacts is comparable. Furthermore, this contact shows double spikes on the 2D-FP and the tip distance is at 2.3 Å in compound 1, whereas the corresponding contact beyond 2.5 Å in compound 2 indicates a relatively weak nature. The contributions of other inter-contacts such as H···C (11% in 1 and 13.9% in 2) and H···S (6.3% in 1 and 7.2% in 2) interactions are comparable. We note that the shortest d i + d e distance for H···C contacts is very similar in these structures and located around 2.8 Å. The inter-H···S contacts are located beyond 3.0 Å, which is longer than the sum of the vdW radii, and these contacts may not have a significant role in stabilizing crystal structures. It is observed that the H···Cl contacts influence the crystal packing, and the relative contribution of this contact is higher in compound 1 (12.5%) than in compound 2 (7.2%). The shortest d i + d e distance is located for this contact at 2.8 Å in compound 2 and 3.0 Å in compound 1, suggesting that the short H···Cl contact plays a significant role in the stabilization of the crystal structure of compound 2.
In contrast to the H···Cl contact, the contribution of H···F contact is higher in compound 2 (9.2%) than in compound 1 (3.3%). The shortest d i + d e distance appears more extended than the sum of the vdW radii of H and F atoms suggesting that the H···F contact plays a minor role in stabilization. Additionally, the inter-C···C contacts, which represent πstacking interactions, contribute about 4.9% to compound 2, and the corresponding contact contributes only 0.3%, indicating its absence in the crystal structure of compound 1. The shape index plot shows the presence or absence of πstacking in compounds 2 and 1, respectively ( Figure 4).  In addition to the intermolecular interaction mentioned above, the contribution of some inter-contacts is relatively less compared to other contacts. For instance, F···C and F···N contacts contribute each about 1.9% and a highly directional S···C(π) chalcogen bond contributes only 2.0% to the crystal packing of compound 1 yet plays a significant role in stabilizing the crystal structure.
2.4. Crystal Packing. 2.4.1. Crystal Packing of Compound 1. In the solid state, the molecules of compound 1 are packed in a columnar way along the crystallographic ac plane (Figure 5a). CLP-PIXEL energy analysis revealed six dimers (D1−D6; Figure 5a−e) in the crystal structure of compound 1, and the intermolecular interaction energies for these dimers range from −12.7 to −3.8 kcal mol −1 . These energies for most dimers are comparable to those calculated by the density functional theory (DFT) method with the B97D3/def2-TZVP level of approximation. As shown in Figure 5a, the primary structural motif (dimer D3) is stabilized by an intermolecular C−H···C(π) interaction between the adamantyl moiety and the phenyl ring. This interaction links the neighboring molecules into a C(10) chain that runs parallel to the crystallographic b axis ( Figure S4). Furthermore, the D3 dimer is stabilized with an 80% dispersion energy component. The most stable dimer D1 stabilizes with the three centered F1A··· C8/N2 contacts. The stabilization of dimer D1 is further supported by an intermolecular C−H···N interaction, which makes a R 2 2 (18) motif. Furthermore, the dispersion and electrostatic energy components contribute 59 and 41%, respectively, toward the stabilization of dimer D1. As shown in Figure 5a, the primary structural motifs in the adjacent columns are interconnected by dimeric motif D1.
The second most stable dimer D2 is formed by two intermolecular C−H···π interactions involving H atoms of the adamantane and triazolothiadiazole core, with a dispersion energy contribution of approximately 74% for stabilization. These C−H···π interactions interconnect the neighboring molecules into a chain which runs parallel to the crystallographic a axis ( Figure S5). Dimer D4 is formed by a highly directional chalcogen bond of type C−S···C(π) interaction, which links the adjacent molecules into a C(6) chain that runs parallel to the crystallographic b axis ( Figure S6). The electrostatic (51%) and dispersion (49%) energy components are nearly equal to the stabilization of this dimer. Dimer D5 is stabilized by an intermolecular C−H···N interaction, and this interaction links neighboring molecules into a C(10) chain that runs parallel to the crystallographic b axis. Moreover, the electrostatic energy contributes about 65% toward the stabilization of this dimer. The least stable cyclic dimer [D6; R 2 2 (22) motif] is formed by intermolecular C−H··· Cl interactions. The H···Cl contact distance is slightly longer (by 0.06 Å) than the sum of the vdW radii of H and Cl atoms and Cl atoms. For stabilization, the dispersion energy contributes about 79%.
In addition, different intermolecular interactions observed in 1 are combined to form a supramolecular self-association in the solid state. The basic structural motif D3 and the dimeric motif D5 combined to assemble a molecular ribbon that runs parallel to the b axis. Further, three molecules of 1 generate a loop utilizing two D3 and a D5 motifs ( Figure 6a). As shown     in Figure 6b, two intermolecular C−H···N interactions (dimers D1 and D5) and a weak C−H···Cl (dimer D6) interaction collectively generate a supramolecular sheet. The adjacent loops formed by this weak C−H···Cl interaction interlinked by C−H···N interactions in both the directions. As mentioned above, the chalcogen bonding (dimer D4) interaction alone generates a molecular chain and forms a molecular ribbon when the chalcogen bond is combined to C−H···N interactions (dimer D5) and this molecular ribbon runs parallel to the b axis ( Figure S7). Furthermore, we calculated the MESP map for compound 1 to understand the molecule's charge distribution and gain further insights into the observed intermolecular interactions in the solid states (Figure 7). On the MESP map, we can see that the protons of the di-substituted phenyl ring had the most positive electrostatic potentials (V s,max ) compared to the H atoms of the adamantane moiety. The most negative electrostatic potentials (V s,min ) are observed for atoms N4 (−45.9 kcal mol −1 ) and N3 (−45.8 kcal mol −1 ) of the triazole ring. Atoms H4A (V s,max : 30.0 kcal mol −1 ) and H5A (V s,max : 29.7 kcal mol −1 ) have participated as donors for the intermolecular C−H···N interactions with triazole N4 atoms having the most negative electrostatic potential. The lone pair of F1A atom with a V s,min value of −10.4 kcal mol −1 is likely to make F···π contact (dimer D1) with the π-hole over the surface of C8−N2 bond with a V s,max value of 10.1 kcal mol −1 . The MESP map also reveals the σ−hole for the S atom with the positive electrostatic potential at the tip of the C7− S1 bond with a V s,max value of 14.3 kcal mol −1 and this atom forms a directional chalcogen bond (C−S···π type) with atom C4A(π). On the MESP map, we observed a negative electrostatic potential with a V s,min value of 3.8 kcal mol −1 near the C4A atom. From the calculated MESP map, we can also see the σ− hole for the Cl atom with a V s,max value of 15.2 kcal mol −1 at the tip of the C−Cl bond and the negative belt around Cl atom corresponds to lone pairs. It is interesting to note that there is no σ−hole interaction involving the Cl atom. However, the lone pair of the Cl atom makes a weak C−H···Cl hydrogen bond (dimer D6) with the H atom of the adamantane moiety.

Crystal Packing of Compound 2.
In the solid state, the molecules of compound 2 are also columnar packed in the crystallographic bc plane, but somewhat different from the crystal packing of molecule 1 (Figure 8). CLP-PIXEL energy analysis identifies four dimers (D1−D4) with significant intermolecular interaction energies (E tot ) observed in this structure (Table 2). These energies range from −11.1 to −6.3 kcal mol −1 . The energy values are comparable to those calculated by the DFT method with the B97D3/def2-TZVP level of theory. The basic structural motif (dimer D1) consists of inversion-related molecules stabilized by π-stacking interactions between phenyl and triazole rings (Figure 8a). The π-stacked dimer D1 is further strengthened by weak lp···π (involving atoms S1 and C14) and C−H···F interactions. The latter hydrogen bond is established slightly longer than the sum of the vdW radii of the H and F atoms and makes an R 2 2 (20) ring motif. The second strong dimer (D2) is also formed between inversion-related molecules and stabilized by C7−S1···N3 chalcogen bonds and further reinforced by a short N···N (triazole···triazole) contact with a distance of 2.89 Å ( Figure  8b). It is important to note that the electrostatic energy contributes 77% toward stabilizing dimer D2. The NCI (noncovalent interaction index) analysis was performed to visualize the nature of the interactions formed in this dimer based on electron density (ρ) and its reduced density gradient (s). The NCI plot shows a strong attractive nature of a chalcogen bond compared to the N···N contact formed in the same dimer. A similar feature has been observed in the related structure reported in the literature. 38 Dimer D3 stabilizes with intermolecular C−H···π interactions and three centered interactions. The Cl atom acts as an acceptor for C−H···Cl interactions and as a donor for a halogen bond of type C−Cl···N interactions. For stabilization of this dimer, the dispersion energy contributes about 76% (Figure 8c). The least stable dimer (D4) in this structure is stabilized by the three-centered C−H···N interaction (R 2 1 (5)ring motif) in which one of the triazole nitrogens (atom N4) acts as an acceptor (Figure 8d). This interaction links the molecules into a chain that runs parallel to the crystallographic a axis. The electrostatic energy contributes about 61% toward the stabilization of this dimer. It is observed that the adjacent dimers of D1 are interlinked by D3 dimer runs that alternately form a molecular ribbon (Figure 8e). We also observed that adjacent dimers of D2 are interconnected by dimer D3 that forms a herringbone-like supramolecular architecture ( Figure 8f). Furthermore, the deformation density map was calculated for dimers D2 and D3 to visualize σ−hole interactions (chalcogen and halogen bonds). Figure 9a shows the deformation density map displaying the charge concentration region (blue) near the N3 atom and charge depletion region (red) closer to the S1 atom. These features suggest the existence of the S···N chalcogen bond in dimer D2 (Figure 9a). Similarly, the charge concentration region is located near the N1 atom and the charge depletion region around the Cl1 atom, suggesting the presence of a highly directional Cl···N halogen bond (Figure 9b).
To further understand the formed interactions in compound 2, MESP map was computed. The MESP map of compound 2 shows that the most positive electrostatic potentials (V s,max ) are observed for protons of the di-substituted phenyl ring, which are higher than H atoms of the phenyl ring. This feature is very similar to that of compound 1. The H atoms of the disubstituted phenyl ring act as donors for intermolecular C− H···N interactions with atom N4 (V s,min : −41.7 kcal mol −1 ). The most positive electrostatic potential surface near atoms H4  Cl···N) with the lone pair of N1 atom having a negative electrostatic potential value of −9.4 kcal mol −1 . As observed in 1, the S1 atom shows σ−hole (17.6 kcal mol −1 ) and makes an attractive chalcogen bond (C−S··· N) with atom N3 having the negative electrostatic potential value of −43.0 kcal mol −1 . Atom H11 with a positive ESP value of 10.2 kcal mol −1 is likely to participate in a weak intermolecular C−H···Cl hydrogen bond with the negative belt of the Cl atom.

Topological Analysis of Intermolecular Interactions in Compounds 1 and 2.
The intermolecular interactions observed in various dimers of compounds 1 and 2 were characterized using the topological properties. The topological properties for the selected interactions in these dimers are summarized in Table 3. The molecular graphs showing the intermolecular interactions at the bond critical points in different dimers of compounds 1 and 2 are illustrated in Figures S8 and S9. The positive values of Laplacian of the electron density (∇ 2 ρ(r) > 0) and total electronic energy density (H(r) > 0) and |−V(r)/G(r)| < 1 indicate all the observed intermolecular interactions in compounds 1 and 2 are closed-shell in nature. In compound 1, the dissociation energy values (D e ) suggest that the C−H···N (dimer D4) interaction is found to be more assertive with a value of 2.9 kcal mol −1 , compared to other interactions. The next strongest interaction is the F···N contact with a D e value of 1.9 kcal mol −1 . The chalcogen bond of type S···π showed a similar strength with that of F···N contact. The other two interactions, C−H···N (dimer D1) and C−H···Cl (dimer D6) showed similar strength.
In compound 2, the chalcogen bond (S···N, D e = 4.2 kcal mol −1 ) and N···N contact (D e = 2.7 kcal mol −1 ) in dimer D1 showed significant strength compared to other interactions. The topological parameters for these contacts are comparable with the same contacts observed in the triazolothiadiazole derivatives reported earlier. 38 We also note that the strength of Cl···N halogen bond (dimer D3) is comparable with the strength of one of the C−H···N interactions observed in dimer D4. As expected, the intermolecular C−H···Cl bond is slightly weaker than C−H···N interactions.
2.6. In Vitro COX Inhibition Assay and Molecular Docking Analysis. The synthesized compounds in this study were subjected to evaluate the anti-inflammation potentials against two important enzymes (COX-1 and COX-2) involved in the inflammatory cell signaling pathway. The in vitro inhibitory activity of compounds 1 and 2 against these enzymes and the selectivity index (SI) are summarized in Table 4. The inhibition concentration (IC 50 ) values for control drugs, namely, celecoxib (selective COX-2 inhibitor) and diclofenac (non-selective COX inhibitor) are compared to assess the activity of compounds 1 and 2. The in vitro data reveals that the binding affinity of compound 2 is nearly in 3fold excess that of compound 1 toward the COX-1 enzyme. In contrast, the binding affinity of compound 1 is in 2-fold excess that of compound 2 toward the COX-2 enzyme. The assay results also indicate that molecule 1 with an adamantane moiety showed marked COX inhibitory activity with 10-fold selectivity toward COX-2 (SI = 9.82) compared to molecule 2 with a phenyl ring in the place of adamantane cage. The result reveals the effect of the substituents (adamantane and phenyl) on the observed COX inhibition activity. Molecules containing 1,2,4-triazole have been reported to exhibit anti-inflammatory potential, and in our earlier reports we described the selective inhibitory potential of 4-(4-chlorophenyl)-3-[(4-fluorobenzyl)sulfanyl]-5-(thiophen-2-yl)-4H-1,2,4-triazole with a COX-2 SI of 1.89. 35 To corroborate in vitro data, we performed an in silico molecular docking analysis for the title compounds and two known drugs (diclofenac and celecoxib) with COX-1/2 enzymes. The GOLD docking program offered a better   performing ChemPLP fitness score for pose prediction, 64 which can be used to assess the binding potential of the compound, and a higher fitness score can be considered to be a better docking pose at the active site. Table 5 summarizes the ChemPLP fitness score for compounds 1 and 2 and two known COX inhibitors. The scores obtained agree well with the in vitro data. The relative affinity (ΔG bind ) of the COX-1/ 2-ligand (title compounds) was obtained using the MM-GBSA calculation. From ΔG bind values, one can observe that compound 2 shows higher affinity toward COX-1 compared to COX-2, whereas both compounds display nearly equal affinity toward the COX-2 enzyme (Table 5). Overall, the molecule with an adamantane moiety shows relatively higher binding potential with the active site of the COX-2 compared to COX-1, which is in good agreement with the in vitro data.
To understand the nature of noncovalent interaction formed between title compounds and COX enzymes, we analyzed the docked complexes (refined models using the MM-GBSA approach) using a PLIP Web server, as mentioned in the experimental section. Compound 1 establishes only hydrophobic interactions with the active site residues (Arg 120, Leu 352, Leu 359, Trp 387, Tyr 385, and Leu 531) of COX-1 (Figure 10a). It is noted that the di-substituted phenyl and adamantane moieties are involved in the hydrophobic interactions. As can be seen in Figure 10b, there is a relatively higher number of interactions (residues involved are as follows: Arg 120, Tyr 355, Leu 359, Phe 381, leu 384, Tyr 385, Trp 387, Leu 531, and Ala 527) established between compound 1 and the COX-2 enzyme compared to the corresponding COX-1 complex. This could a possible reason for the higher selectivity of compound 2 toward COX-2 enzymes.
Compound 2 showed relatively higher affinity compared to compound 1 for COX-1 according to the calculation of the MM-GBSA and the in vitro assay. It can be seen in Figure 10c, in addition to the higher number of hydrophobic interactions, the T-shaped π-stacking interaction between the di-substituted phenyl ring of 1 and a conserved Tyr 385 residue of COX-1 and a short F···O contact (backbone O atoms of the backbone of Met 522 and Ile 523 are involved). The unsubstituted phenyl ring also establishes hydrophobic interactions with the key residues. In the COX-2-compound 2 complex, both disubstituted and unsubstituted phenyl rings participate in hydrophobic interactions with the active site residues of COX-2 (Figure 10d). In addition, a hydrogen bond (O−H··· N) formed between the central triazole ring and the Tyr 355 residue.

CONCLUSIONS
The present study demonstrated the interplay of hydrogen, halogen, and chalcogen bonding in addition to π-stacking interactions in two compounds containing a [1,2,4]triazolo- [3,4-b] [1,3,4]thiadiazole scaffold with adamantane and phenyl substituents. The X-ray structures and various theoretical tools were used to characterize the noncovalent interactions. Molecular conformation is locked by weak intramolecular C−H···N interactions in both compounds and one of the N atoms of the thiadiazole ring was involved as an acceptor. In related reported structures, the di-substituted ring was coplanar with the central triazolothiadiazole ring. The potential energy surface analysis revealed the preferred orientation of the disubstituted phenyl ring, which is in good agreement with the X-ray conformation of the investigated compounds. Through the crystal packing of compounds 1 and 2 are very similar (columnar packing mode), their basic structural motifs are different. The basic structural motif observed in 1 is stabilized by C−H···C(π), whereas the π-stacking interaction generates the basic structural motif in compound 2. The solid-state structure of 1 was also stabilized by C−H···N, C−H···π, and C−H···Cl hydrogen bonds and a directional C−S···π chalcogen bond and unconventional short F···C/N contacts. Similarly, the solid-state structure of 2 was also stabilized by C−H···N, C−H···π, and C−H···Cl hydrogen bonds and directional attractive σ−hole interactions such as the C−Cl··· N halogen bond and C−S···N chalcogen bond. In compound 2, S(lp) ···C(π) and a short N···N contacts play a supporting role in the stabilization of certain molecular dimers. The electrostatic potential surface map and the deformation density map revealed the characteristic features of σ−hole interactions. The contribution of H···H/Cl/N contacts was significantly varied due to the effect of phenyl/adamantane substituents, as evident from the HS analysis. The topological analysis suggested that the C−H···N interaction was found to be stronger in 1, whereas C−S···N chalcogen bond was stronger in compound 2. In vitro assay and in silico molecular coupling and subsequent calculation of MM-GBSA free energy supported that the compound with an adamantane scaffold was 10-fold more selective against the COX-2 enzyme compared to compound 1. This study identified that terminal groups (di-substituted and phenyl/adamantane moieties) are important in establishing noncovalent interactions with the active site residues of the COX-1/2 enzymes. In the future, these groups could be altered to improve the affinity and selectivity against these enzymes.

Chemicals and Instruments.
Melting points (°C, uncorrected) were measured in open glass capillaries using a Stuart SMP30 electro-thermal melting point apparatus. NMR spectra were obtained on a Bruker 400 Avance III at 400.20 MHz for 1 H and 100.63 MHz for 13 C using DMSO-d 6 as a solvent. Electrospray ionization mass spectroscopy (ESI-MS) was performed on an Agilent 6410 triple quad tandem mass spectrometer at 4.0 kV for positive ions. All chemical and solvents (HPLC grade) were purchased from commercial suppliers and were used without further purification. Monitoring the reactions and checking the purity of the final products were carried out by thin-layer chromatography using silica gel pre-coated aluminum sheets (60 F 254 , Merck) and visualization was done with ultraviolet light (UV) at 365 and 254 nm.

Single-Crystal X-ray Diffraction.
Single-crystal X-ray diffraction data were collected for crystals of compounds 1 and 2 at 160 (1) K on a Rigaku OD SuperNova/Atlas area-detector diffractometer using Cu Kα radiation (λ = 1.54184 Å) from a micro-focus X-ray source and an Oxford Instruments Cryojet XL cooler. The selected suitable single crystals were mounted using polybutene oil on a flexible loop fixed on a goniometer head and immediately transferred to the diffractometer. Preexperiment, data collection, data reduction, and analytical absorption correction 65 were performed with the program suite CrysAlisPro, version 1.171.40.68a (Rigaku Oxford Diffraction, England, 2019). Using Olex2, 66 the structure was solved with the SHELXT 67 small-molecule structure solution program, and it was refined with the SHELXL2018/3 program package 68 by full-matrix least-squares minimization on F 2 . In compound 1, the substituted benzene was disordered over two sets of positions (labeled A and B) with site-occupancy factors of 0.7764 (19) and 0.2236 (19). All H atoms were placed in calculated positions (C−H = 0.95−0.99 Å) and were constrained to ride on their parent atoms, with U iso (H) = 1.2U eq (C). The PLATON program was used to check the results of the X-ray analysis. 69 Crystal packing and molecular dimers were produced using the MERCURY program. 70

Theoretical Calculations.
For all calculations, we used the crystal structure geometry with the normalized H positions (C−H = 1.083 Å). Energy framework analysis 71,72 was performed on the CrystalExplorer-17.5 program 73 using the B3LYP/6-31G(d,p) level of approximation. The HS 74 and 2D-FP plots 75 were obtained from the CrystalExplorer-17.5 program. The deformation density was calculated for selected dimers at the HF/6-311G level using the CrystalExplorer program. We calculated the intermolecular interaction energies (E tot ) between molecular pairs using the CLP-PIXEL program. 76,77 For this computation, the electron density was obtained with the Gaussian 09 program 78 using MP2/6-31G** levels of theory. Furthermore, the accurate complexation energies of molecular pairs identified from the CLP-PIXEL energy calculation were calculated using the B97D3/def2-TZVP level of theory, and then, these complexation energies were corrected (ΔE cp ) for the basis set superposition error using the counterpoise method. 79 The topological analysis of the selected dimers was performed within the framework of Bader's QTAIM approach using the AIMALL program. 80 The wave functions were calculated at the M062X-D3/def2-TZVP level of theory for topological analysis. The NCIplot index was also used to characterize the nature of noncovalent interactions. 81 The MESP surfaces have been constructed using the 0.001 a.u. isosurface with the program WFA-SAS. 82 4.5. In Vitro COX Inhibition Assay. The in vitro inhibitory activity of compounds 1 and 2 against cyclooxygenases COX-1 and COX-2 was evaluated using an enzyme immunoassay kit of Cayman Chemical, Ann Arbor, MI, USA (catalog no. 560131). The preparation of the reagents and the testing procedure was performed according to the manufacturer recommendations using various concentrations (0.01− 100 μM) of compounds 1 and 2, celecoxib (selective COX-2 inhibitor), and diclofenac (nonselective COX inhibitor) in dimethylsulfoxide (DMSO). The concentration that causes 50% enzyme inhibition (IC 50 ) was calculated from the concentration inhibition response curve and the SI was calculated by dividing IC 50 COX-1 by IC 50 COX-2.
4.6. Molecular Docking and Protein−Ligand Interaction Analysis. Docking studies were performed for compounds 1 and 2 and two control inhibitors, namely, celecoxib and diclofenac, using GOLD software. 64 Crystallographic structures of COX-1 (PDB ID: 3KK6; Ovis aries cyclooxygenase-1 complexed with celecoxib) and COX-2 (PDB ID: 5IKR; human cyclooxygenase-2 complexed with mefenamic acid) were obtained from Protein Data Bank (www.rcsb.org). Prior to docking simulations, the water molecules and cocrystallized ligands except celecoxib and mefenamic acid inhibitors were removed. All protein atoms within 6.0 Å of celecoxib and mefenamic acid inhibitors were used for the binding site definition.
The highest scoring poses (ChemPLP) of compounds 1 and 2 were also subjected to the MM-GBSA method to calculate the free energy (ΔG bind ) of the binding of ligands to target proteins using the Schrodinger suite (Schrodinger Release 2022-2, Schrodinger, LLC, New York, NY, 2021). In this calculation, the title compounds and active residues within 6 Å from the ligands were treated to be flexible and the remaining protein residues were kept frozen. Protein−ligand interactions for the docked complexes were analyzed using the PLIP Web server (protein−ligand interaction profiler). 83 ■ ASSOCIATED CONTENT
Topological properties for compounds 1 and 2, molecular graphs for different dimers of 1 and 2, and 2D-fingerprint plots for various intercontacts observed in 1 and 2 (PDF)