**3. Results**

## *3.1. Structural Description*

According to the single crystal X-ray diffraction data, molecules **1**–**4** comprise the 1,3,5-triazacyclohexane ring bearing three substituents at the nitrogen atoms (see Figure 1).

**Figure 1.** Single crystal X-ray diffraction structures of triazinanes **1**–**4**.

The asymmetric unit contains one molecule of each triazine. The general geometrical features of these systems are similar, which are the slightly distorted chair conformation of the six-membered ring, with the *N*-alkyl substituents (*n*-Bu or *t*-Bu) occupying the axial position and the two *N*-sulfamide fragments occupying the sterically favorable pseudo-equatorial orientation (see Table S2 for torsion angles of *N*-substituents). The atoms C2, C6, N3, and N5 of the triazinane cycle in all structures lie nearly in one plane (the deviation of one atom from the plane of the other three is less than 0.03 Å), while the deviation of N1 and C4 atoms from this plane range from −0.633(2) Å to −0.661(2) Å and from 0.643(14) Å to 0.712(2) Å, correspondingly, therefore the molecules have a classical chair-conformation of the central heterocycle (see Scheme 2 for atom numbering scheme). All CH2-N bond lengths and CH2–N–CH2 bond angles are typical for 1,3,5-triazinanes and are listed in the corresponding tables in (Supplementary Materials Table S2–S15). The torsion angles CNCN in the triazacyclohexane ring are close to 60◦ (see Table S2).

The most intriguing and distinguishing feature of the triazinanes under discussion is the unprecedented geometry of the N3 and N5 nitrogen atoms of the sulfonamide fragments in the *N*-butyl substituted heterocycles **1** and **2**. As is generally known, the nitrogen atom in a sulfonamide group can adopt both *sp*<sup>2</sup> and *sp*<sup>3</sup> hybridization depending on substituents at the nitrogen atom [22–24]. However, according to the data of the CCDC, there are no known examples of 1,3,5-triazinanes or other saturated six-membered azaheterocycles simultaneously possessing two chemically identical sulfonamide nitrogen atoms in di fferent hybridization. Analysis of the values of the sums of valence angles at nitrogen atoms in positions 3 and 5 allows to clearly identify atoms in *sp*<sup>2</sup> or in *sp*<sup>3</sup> hybridization (Figure 1 and Table 2). The N3-atoms in compounds **1**–**3** are *sp*2-hybridizated and, as a result, they assume the flat trigonal configuration (the sum of the angles is close to 360◦). N5-Atoms, chemically equivalent to N3-atoms, in the same molecules adopt the tetrahedral configuration (the sum of the angles lies in the diapason of 344–351◦) and, therefore, are *sp*3-hybridizated. This is most clearly seen in the examples of compounds **1** and **2**. The presence of the *tert*-butyl group at the N1 position in compounds **3** and **4**, probably due to its high steric volume, symmetrizes the molecules, leveling the di fference between both sulfonamide nitrogen atoms in a crystal. This is also observed in the equalization of the S–N distances of the sulfonamide groups in compounds **3** and **4** compared to **1** and **2**, see Table 2. The largest di fference between both S–N distances is observed in compound **2**, i.e., 0.016 Å. In fact, the short S2–N5 distance is an indication of a partial double bond character, in agreemen<sup>t</sup> with the *sp*2-hybridization.

**Table 2.** Sums of angles at the sulfonamide nitrogen atoms N3 and N5 and S–N distances in 1,3,5-triazinanes **1**–**4**.


Compound **1** crystalizes in the monoclinic crystal system in the space group *P*21/*<sup>c</sup>*. Selected bond lengths, angles and torsional angles are shown in Tables S2–S5. The S2–C21 and S1–N3 bond lengths of 1.758(1) and 1.642(1) Å respectively are in agreemen<sup>t</sup> with a single bond character of these bonds [24]. The S–O bond distances of the sulfamide moiety are in the range 1.432(1)–1.436(1) Å, which is typical for *N*-sulfamides. The SO2 group has typical angles if compare with CSD data for N-sulfamides (Figure S1 in Supplementary Materials): O–S–O angle is around 120◦ [O1–S1–O2 is 119.95(6)◦], while N–S–O angles are between 105◦ and 115◦ [N3–S1–O2 is 106.07(6)◦, N3–S1–O1 is 106.63(5)◦].

The crystal structure of this compound exhibits interesting assemblies in the solid state (see Table S3 for H-bonds). For instance, Figure 2a shows a self-assembled dimer dominated by C–H··· O interactions where the methyl group in para acts as a H-bond donor. The acidity of these protons is higher than usual for a methyl group due to the presence of the electron withdrawing sulfamide group. Quite remarkable is the ternary assembly shown in Figure 2b, where the aromatic π-cloud interacts simultaneously with the methyl group at one side and an aromatic C–H bond at the opposite side, thus forming a C–H··· <sup>π</sup>···H–C assembly. It is worth noting that the C–H··· π distances are very short (2.60 and 2.70 Å) thus, confirming their relevance in the solid state of this compound.

Compound **2** crystalizes in the monoclinic crystal system in the space group *P*21/*n* and the main di fference with compound **1** is the absence of methyl groups. In addition, the sulfonamide groups are attached to the 1,3,5-triazacyclohexane ring in di fferent orientations as reflected by the C6–N5–S2–C21 and C2–N3–S1–C11 torsional angles of −105.0(1)◦ and <sup>−</sup>65.5(1)◦, respectively. It also forms self-assembled dimers in the solid state, where both C–H··· π and C–H···O (Table S3) interactions are established, as shown in Figure 3a. Moreover, it also forms infinite 1D supramolecular chains in the solid state promoted by C–H···O interactions involving the butyl chain and the sulfonamide group (see Figure 3b).

**Figure 2.** (**a**) Self assembled dimer in compound **1**. (**b**) C–H···π···H–C assembly in the solid state of structure **1**. Distances in Å. For the C–H···<sup>π</sup> interactions, the distances are measured from the H-atom to the ring centroid.

**Figure 3.** (**a**) Self-assembled dimer of compound **2**. (**b**) 1D supramolecular chain in **2**. Distances in Å. The C–H···<sup>π</sup> interaction distances are measured from the H-atom to the closest C-atom of the ring.

Compound **3** crystalizes in the monoclinic crystal system in the space group *P*21/*n* and, similarly to compound **1**, also forms ternary assemblies where the same aromatic ring establishes C–H···<sup>π</sup> interactions at both sides of the ring, thus forming a C–H···π···H–C assembly (see Figure 4). This compound also forms 1D supramolecular chains in the solid state where the H-atoms of the triazinane ring interact with the O-atoms of the sulfonamide groups (Table S3), as shown in Figure 4b. The formation of this assembly is facilitated by the relative orientation of the *p*-methyl-benzene-sulfonamide groups.

**Figure 4.** (**a**) C–H···π···H–C assembly in the solid state of structure **3**. Distances in Å. (**b**) 1D supramolecular chain in **3**. Distances in Å. For the C–H···<sup>π</sup> interactions, the distances are measured from the H-atom to the ring centroid.

Compound **4** crystallizes in the orthorhombic crystal system in the space group *P*212121 and also exhibits several motifs in the solid state that are mainly dominated by C–H···O interactions (Table S3). As examples, two motifs are given in Figure 5, one corresponds to a discrete dimer where three H-bonds are formed and the other one to a 1D supramolecular polymer also governed by the formation of C–H···O bonds involving the aromatic H-atom located in para to the sulfonamide group (Figure 5b).

**Figure 5.** (**a**) Self-assembled dimer of compound **4**. (**b**) 1D supramolecular chain in **4**. Distances in Å.

## *3.2. Hirshfeld Surfaces*

The Hirshfeld surface analysis is a very convenient tool for analyzing intermolecular interactions. The HS surfaces mapped over *d*norm property are displayed in Figure 6 highlighting the main intermolecular interactions and scheme of labels. The patterns of intermolecular interactions are similar in all structures, which prompted us to evaluate the contributions of the weak non-covalent contacts in the supramolecular assembly, as well as the importance of C–H···<sup>π</sup> interactions in stabilization of the crystal packing. The 2D fingerprint plots (Figure 7) of the molecules illustrate significant differences between the intermolecular interaction patterns. The surfaces are shown as transparent to allow the visualization of the molecules. Contacts with distances equal to the sum of van der Waals (vdW) radii are represented as white regions and contacts with distances shorter than and longer than the vdW radii are shown as red and blue colors, respectively.

**Figure 6.** Hirshfeld surfaces mapped over *d*norm function for compounds **1**–**4**. The labels are discussed in the text. For compound (**2**), the second molecule is rotated by 180◦ around the vertical axis of the plot.

**Figure 7.** Full and decomposed 2D fingerprint plots for compounds: (**a**) **1**; (**b**) **2**; (**c**) **3** and (**d**) **4**.

The vdW forces (H···H contacts) have the largest contribution to the HS, and they are highlighted in the scattered middle points in the fingerprint plots with a minimum value of (de + d1) ~ 2.2 Å (Figure 7), which is the sum of the vdW radii. All red areas that are visible on the surfaces mapped over *d*norm function correspond to C–H···O contacts. For (**1**), the largest bright-red spot labeled 1 on the HS shows O···H/H···O contact associated with C17–H17C···O1 interaction, which constitutes the strongest among all interactions present in this compound. Two medium sized red spots labeled as 2 and 3 are associated with C15–H15A···O1 and C27–H27C···O4, respectively. These interactions are also visible as symmetrical sharp spikes centered at (de + di) - 2.4 Å in the fingerprint plots (Figure 7a) with 25.2% contribution to the Hirshfeld contact surface. The intermolecular C25–H25A···N1 contact is visible in the HS as a red spot labelled 4, which comprises 2.3% of the total HS area. The HS of (**2**) mapped over *d*norm function (see Figure 6) shows four red spots, indicating the presence of C–H···O hydrogen bonds

[C13–H13A···O1 (**1**), C26–H16A···O4 (**2**), C10–H10B···O1 (**3**), C14–H14A···O2 (**4**)]. The decomposed fingerprint plots (Figure 7b) show that intermolecular O··· H/H···O contacts contribute 27.7% to the total HS area. The O··· H/H···O contacts appeared as sharp spikes with (de + di) - 2.35 Å. In the HS of **(3)**, the O··· H/H···O contacts (Figure 6) are visible as six red spots attributed to C6–H6A···O3 (labeled 1), C10–H10B···O3 (labeled 1), C15–H15A···O4 (labeled 2) and C17–H17C···O4 (labeled 2). These interactions comprise the 24.1% of the HS area. A similar behavior was observed in the HS mapped over *d*norm function for (**4**), which the six red spots observed (Figure 6) are attributed to C26–H26···O4 (1), C16–H16···O4 (2), C24,H24···O2 (3), C12–12···O3 (4) and C6–H6A···O3 (5). Intermolecular interactions O··· H/H···O are observed around 2.3 Å which is slightly shorter than those of other compounds with 30% contribution to the Hirshfeld contact surface.

As was described previously, the structure of (**1**) is also stabilized by C–H··· π interactions. The red area labeled 5 in the HS mapped over *d*norm is attributed to C9–H9B··· π. These C–H··· π interactions are also evident from a pair of "wings" in the top left and bottom right region of the fingerprint plots for compounds **1**–**4** (Figure 7). The shape of the wings and the sum of *d*e and *d*i show the importance of this interaction. The decomposition of the fingerprint plots shows that the C··· H/H···C contributions comprising 17.4%, 17.2%, 18.8% and 18.7% of the total HS for each molecule of **1**, **2**, **3** and **4**, respectively.

## *3.3. DFT Calculations*

The DFT study is focused to analyze the supramolecular assemblies commented above in Figures 2–5, where combinations of C–H··· π and C–H···O H-bonding networks are commonly formed in compounds **1**–**4**. First of all, the molecular electrostatic potential (MEP) surfaces of compounds **1** and **2** have been computed in order to analyze the electron rich and electron poor regions of the molecules. The surfaces are represented in Figure 8 and it can be observed that the most negative regions correspond to the O-atoms of the sulfonamide group. The N-atoms of the triazinane ring are not good H-bond acceptors, likely because either the lone pair is delocalized into the SO2-group, in accordance with the S2–N5 and S1-N3 bond lengths of 1.632(1) and 1.642(1) Å, respectively indicating a double bond character of these bonds. The most positive region corresponds to the middle of the three axial H-atoms of the triazinane ring (+27 kcal/mol). The aromatic H-atoms and the CH3 substituents also present positive MEP values (+23 and +20 kcal/mol, respectively). Finally, the MEP value over the aromatic rings is negative (−8 kcal/mol), thus adequate for establishing C–H··· π interactions. The MEP analysis evidences that the most favored interactions from an electrostatic point of view are those involving the O-atoms as electron donors and either aromatic or triazinane protons as electron acceptors.

**Figure 8.** MEP surfaces for compounds **1** (**a**) and **2** (**b**). The energies at selected points of the surfaces are given in kcal/mol.

We have selected two supramolecular assemblies commented above in Figure 2 to analyze the energetic features of the H-bonds and C–H···<sup>π</sup> interactions in **1**. The QTAIM distribution of critical points and bond paths are also given in Figure 9. The existence of a bond CP and bond path connecting two atoms is a good indicator of interaction [21]. For the self-assembled dimer (Figure 9a), in addition to the symmetrically related H-bonds (characterized by a bond CP (critical point) and bond path interconnecting the H and O-atoms), the QTAIM analysis reveals the existence of a π···π stacking interaction that further stabilizes the formation of the dimer. The dimerization energy is moderately strong (ΔE1 = −9.1 kcal/mol) due to the contribution of both H-bonds and the π-stacking. We have also analyzed the other motif, where C–H···<sup>π</sup> interactions are established. The interaction energy is very strong (ΔE2 = −20.4 kcal/mol), because in addition to the C–H···<sup>π</sup> contacts (two bond CPs and bond paths connect two aromatic H-atoms to two carbon atoms of the adjacent ring) an intricate network of H-bonds is established where six C–H···O and one C–H···N contacts are formed, which are highlighted in Figure 9b by yellow circles.

**Figure 9.** (**<sup>a</sup>**,**b**) Distribution of bond, ring and cage critical points (green, yellow and blue spheres, respectively) and bond paths in two dimers of complex **1**. The C–H···O bonds are highlighted in (**b**).

Figure 10a shows the self-assembled dimer of compound 2 where up to eight C–H···O contacts are established between either aromatic or aliphatic H-atoms and the O-atoms of sulfonamide (each one characterized by a bond CP and bond path, see yellow circles in Figure 10a). Moreover, two symmetrically distributed C–H···<sup>π</sup> interactions are also present and characterized by a bond CP and bond path connecting the aliphatic H-atom to one C-atom of the aromatic ring. As a consequence of this combination of interactions, the dimerization energy is very large ΔE3 = −16.7 kcal/mol, thus confirming the importance of this motif in the solid state of compound **2**. Figure 10b shows a dimer extracted from the infinite 1D chain represented in Figure 3b. In this case, the interaction energy is modest (ΔE4 = −3.4 kcal/mol) because only one H-bond is established. The distribution of bond CPs and bond path also reveals the existence of van der Waals interactions between the alkyl chain and the aromatic ring.

Figure 11 shows two dimers of compound **3**, one corresponds to the C–H···<sup>π</sup> assembly commented above in Figure 4, where in addition to the C–H···<sup>π</sup> interaction (characterized by a bond CP and bond path) the assembly is further characterized by a C–H···O bond involving the methyl group. This assembly presents a modest interaction energy of ΔE4 = −5.9 kcal/mol. In contrast, the dimer shown in Figure 11b, extracted from the infinite 1D assembly, exhibits a strong interaction energy (ΔE5 = −17.3 kcal/mol) due to the formation of four C–H···O contacts, which are characterized by a bond CP and bond path (see yellow circles in Figure 11b). The strong interaction energy agrees well with the MEP surface analysis commented above, since the H-bond donors belong to the triazinane

ring that exhibit the most positive MEP values. Moreover, the H-bond acceptors are the O-atoms of the sulfonamide groups that present the most negative MEP values (see Figure 8).

**Figure 10.** (**<sup>a</sup>**,**b**) Distribution of bond, ring and cage critical points (green, yellow and blue spheres, respectively) and bond paths in two dimers of complex **2**. The C–H···O bonds are highlighted.

**Figure 11.** (**<sup>a</sup>**,**b**) Distribution of bond, ring and cage critical points (green, yellow and blue spheres, respectively) and bond paths in two dimers of complex **3**. The C–H···O bonds are highlighted.

Finally, Figure 12 shows the dimeric motifs analyzed in compound **4**. The dimer of Figure 12a presents an intricate combination of C–H···O bonds in addition to two C–H···<sup>π</sup> interactions involving the t-butyl group. As a consequence of the formation of six concurrent H-bonds, the dimerization energy is very large, ΔE7 = −20.2 kcal/mol. Figure 12b shows the dimer extracted from the 1D supramolecular polymer (see Figure 5b), which presents a modest interaction energy due to the formation of a single H-bond along with van der Waals contacts between the aromatic and aliphatic C–H bonds.

**Figure 12.** (**<sup>a</sup>**,**b**) Distribution of bond, ring and cage critical points (green, yellow and blue spheres, respectively) and bond paths in two dimers of complex **4**. The C–H···O bonds are highlighted.

## **4. Concluding Remarks**

In summary, the synthesis and single crystal X-ray diffraction characterization of four N-substituted 1,3,5-triazinanes are reported along with a detailed analysis of the noncovalent interactions observed in the solid state. All complexes have in common the formation of several motifs characterized by a network of C–H···O interactions that exhibits very strong binding energies as a consequence of these cooperative H-bonds. Moreover, several structures also form interesting C–H···π···H–C ternary assemblies that have been described in detail. Besides, the MEP surfaces have been used to rationalize the noncovalent interactions and the QTAIM method to confirm the existence of the intricate combinations of H-bonds. Finally, the Hirshfeld surface analysis provides further evidence for the importance of C–H···<sup>π</sup> and C–H···O in the crystal packing of compounds **1**–**4**. We assume, that observed case of intramolecular *sp2-sp3* disequalization makes the corresponding family of *<sup>N</sup>*-alkyl-*N*,*N"*-substituted triazinanes an interesting object for research in the domain of local molecular disorder in organic crystals [25].

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4352/10/5/369/s1, Figure S1: Histogram of O–S–O angle in N-sulfamides (a) N–S–O angle in N-sulfamides (b) according CSD analysis. Table S1: Data and structure refinement for **1**–**4**. Table S2: Selected distances (Å), angles (◦) and torsion angles (◦) for compounds **1**–**4**. Table S3: Hydrogen bonds for compounds **1**–**4**. Table S4: Bond Lengths for **1**. Table S5: Bond Angles for **1**. Table S6: Torsion Angles for **1**. Table S7: Bond Lengths for **2**. Table S8: Bond Angles for **2**. Table S9: Torsion Angles for **2**. Table S10: Bond Lengths for **3**. Table S11: Bond Angles for **3**. Table S12: Torsion Angles for **3**. Table S13: Bond Lengths for **4**. Table S14: Bond Angles for **4**. Table S15. Torsion Angles for **4**.

**Author Contributions:** Conceptualization, D.M.G., V.P.Z., N.L.M. and A.F.; formal analysis, D.M.G., K.R.B. and A.A.S. investigation, D.M.G., A.F., F.I.Z., M.V.G. and A.V.K.; Methodology, M.S.G.; writing—original draft preparation, D.M.G. and A.F.; writing—review and editing, D.M.G., A.F., F.I.Z. and A.V.K.; supervision, A.F. and F.I.Z.; funding acquisition, D.M.G., A.F. and F.I.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partly funded by the MICIU/AEI of Spain (project CTQ2017-85821-R FEDER funds). The publication has been prepared with the support of the "RUDN University Program 5-100" and with the support of the Russian Foundation for Basic Research (project No 19-03-00807 A).

**Acknowledgments:** We thank the "Centre de Tecnologies de la Informació" (CTI), Universitat de les Illes Balears for computational facilities.

**Conflicts of Interest:** The authors declare no conflict of interest.
