Lucía
González
a,
Iván
Marín
a,
Rosa M.
Tejedor
b,
Joaquín
Barberá
a,
Pilar
Romero
a,
Alberto
Concellón
a,
Santiago
Uriel
*a and
José L.
Serrano
*a
aDepartamento de Química Orgánica, Instituto de Nanociencia y Materiales de Aragón (INMA), CSIC-Universidad de Zaragoza, 50009 Zaragoza, Spain. E-mail: joseluis@unizar.es; suriel@unizar.es
bCentro Universitario de la Defensa, Academia General Militar, 50090 Zaragoza, Spain
First published on 4th December 2023
We have synthesized two geometric isomers of a cyclohexane-5-spirohydantoin derivative (1,3-diazaspiro[4.5]decane-2,4-dione) incorporating a hydrophobic phenyl 3,4,5-tris(dodecyloxy)benzoate unit at position 8. Separation of these diastereomers was accomplished through silica gel flash chromatography. The interplay of intermolecular hydrogen bonding and micro-segregation between the polar hydantoin unit and nonpolar aliphatic chains within the molecule endows them with remarkable self-assembly capabilities, both in solution and in the solid state. These hydantoin derivatives spontaneously form rosette-shaped structures composed of six molecules. In the solid state, these compounds display hexagonal columnar liquid crystal phases, with hydrogen-bonded disks as their fundamental building blocks. Similarly, when exposed to apolar solvents such as cyclohexane or dodecane, they adopt a columnar arrangement, resulting in gel formation comprising nanoscale fibers that intricately interlace to form a network. Remarkably, the two isomers exhibit markedly different properties. The major isomer behaves as a glassy liquid crystalline material, while the minor one exhibits liquid crystalline behavior with a high propensity to crystallize. Our experimental findings, in combination with theoretical studies, underscore the fundamentally distinct supramolecular organizations present in these isomers, shedding light on their unique self-assembling properties.
The novelty of this supramolecular approach lies in the simplicity of molecular geometry and the versatility of structural modifications. Indeed, a wide variety of LC organizations can be constructed by combining a small number of components, which typically include both rod-like (calamitic) and disc-like (discotic) mesogenic complexes. Many examples of supramolecular macrocycles have been reported based on self-complementary units such as melamine, cyanuric acid, folic acid, barbituric acid derivatives, or nucleobases.20–29 Specifically, “rosettes” (i.e., H-bonded macrocycles) have garnered significant attention recently as they enable the formation of complex columnar architectures with responsive behaviors.30–34 In this context, our research group has undertaken pioneering work on H-bonded rosette structures, consisting of six melamines and six carboxylic acids, which self-assemble into highly ordered columnar organizations.35,36
Nonetheless, hydantoin derivatives, or imidazoline-2,4-dione compounds, which feature a structure with two H-bond accepting carbonyl groups and two H-bond donating amino groups (two carbonyl and two N–H acceptors and donors), have been less frequently considered for building supramolecular architectures. These molecules are essential intermediates in amino acid synthesis and have been demonstrated to possess various biological properties, making them valuable compounds.37 This has led to intense research activity and the development of novel synthetic methods for accessing these valuable molecules.38 The unique geometry of hydantoins has also been exploited to prepare organometallic derivatives39–42 and to create different crystalline architectures linked by hydrogen bonds.43–46 In some of these derivatives, the introduction of a cyclopentane or cyclohexane unit at position 5 of the hydantoin heterocycle introduces additional steric factors that influence the supramolecular organization and properties of these molecules.47,48
Herein, we present two novel hydantoin derivatives with unprecedented self-assembly properties. Two geometric isomers of 1,3 diazaspiro[4.5]decane 2,4 dione ((Z)-4 and (E)-4 in Scheme 1) have been prepared bearing a hydrophobic phenyl 3,4,5-tris(dodecyloxy)benzoate unit at position 8. The intermolecular H-bonding and microsegregation between the polar (hydantoin unit) and nonpolar (aliphatic chains) regions of the molecule lead to an excellent self-assembly capability in solution and in the solid state. These hydantoin derivatives self-organize into rosette-shaped structures showing columnar hexagonal mesophases, in which six molecules form the hydrogen-bonded disk. In addition, in the presence of apolar solvents such as cyclohexane or dodecane, they form columnar-type organizations that give rise to gels made up of fibers of several tens of nanometers that entangle to form a network.
Fig. 1 Chemical structures of the diastereomers (Z)-4 and (E)-4 and the possible corresponding two conformers depending on the cyclohexane chair conformations. |
Fig. 2 (a) 1H-NMR spectrum (CDCl3, 298 K) of mixture 4 of diastereomers 4 diastereomers at a concentration of 40 mM. (b) Signals corresponding to the 1H-NMR spectrum (CDCl3, 298 K) of mixture 4 at different concentrations. The signals corresponding to the major isomer appear, linked by the red line in Fig. 2b and surrounded in red in Fig. 2a. In the case of the minor isomer, these lines and circles appear in blue. |
Fig. 2a shows the 1H NMR spectrum of the diastereomeric mixture of 4 in CDCl3 (40 mM), in which the two protons of the NH groups of the hydantoin ring appear as broad signals at different chemical shifts depending on the diastereomer. In particular, the NH protons of the main diastereomer appear at around 9.36 and 8.21 ppm, while in the minor diastereomer, the NH signals appear at around 9.09 and 6.71 ppm. The integration of these signals shows a 75 and 25% ratio. Moreover, we observed a concentration-dependent behaviour in the 1H-NMR spectra, wherein the NH signals originating from the hydantoin unit sharpen and shift to higher field values as the concentration decreases. In contrast, the other signals remain unaffected (see Fig. 2b). This phenomenon arises due to intermolecular hydrogen bonding.
The 1H-NMR spectra obtained at room temperature confirmed that, in both diastereomers, only the conformer that has the bulky aromatic group at the equatorial position appears. The configuration of H4 can be elucidated by studying the coupling constants between H4 and H3ax. The proton signal at approximately 2.6 ppm, corresponding to H4, exhibits a substantial coupling with H3ax (J ≈ 12 Hz), indicating that both H4 and H3ax are at the axial position. It is well-known that axial–axial couplings exhibit larger coupling constants (3Jax–ax = 10–12 Hz) compared to those of equatorial–axial (3Jax–eq = 2–6 Hz) or equatorial–equatorial (3Jeq–eq = 2–5 Hz) couplings. To further elucidate the relative configuration of the two hydantoins within the substituted cyclohexanes, ROESY experiments were conducted (Fig. 3). The 1H–1H ROESY experiment of diastereomer (Z)-4 reveals cross peaks between the hydantoin's NH and H3ax and H2ex, whereas the 1H–1H ROESY spectrum of (E)-4 indicates that hydantoin's NH exhibits cross peaks only with H2ax. Consequently, the NMR studies allow the assignment of (Z)-4 as the major diastereomer (75%) and (E)-4 as the minor one (25%) (Fig. S1–S10 in the ESI†). This observation aligns with the fact that the Bucherer–Bergs reaction of cyclohexanone derivatives predominantly yields hydantoins with the NH group at the axial position.51 In addition, a search of the Cambridge Structural Database (November CSD 2022 release) with both unsubstituted hydantoin NH groups and carbocyclic spiro fusion at the C5 position yielded 54 hits (excluding duplicate structure determinations), 19 of them being cyclohexanespiro 5 hydantoin derivatives.52 Importantly, all of these cyclohexane-5-spirohydantoin derivatives feature hydantoin's NHα, located close to cyclohexane, at the axial position. Besides this, in the ROESY studies, the NOE effect between the H6 proton of the benzene ring and the H3ax protons of the cyclohexane ring is observed in both diastereomers, confirming the substitution of the aromatic ring at the equatorial position. From these data, it is deduced that for the (Z)-4 diastereomer only the NHax–Φeq conformer appears, while for the (E)-4 diastereomer only the NHeq–Φeq conformer appears.
Fig. 4 1H DOSY NMR spectra at 298 K of (Z)-4 in CDCl3 (bottom trace in black) and in C6D12 (top trace in blue). The projection corresponds to the 1H-NMR spectrum in CDCl3. |
Diastereomer (E)-4 is not soluble in pure cyclohexane but it is soluble when adding 10% (v/v) chloroform (cyclohexane–chloroform, 9:1). The aggregation in these two solvents is different since in the latter case it tends to form dimers. There may be a competition between the formation of dimers in chloroform and rosettes in cyclohexane, so only a qualitative approximation of the aggregation was performed (see Fig. S14 in the ESI†).
Interestingly, during some NMR studies of both diastereomers in cyclohexane, we observed the formation of gels. Consequently, we conducted 1H–1H NOESY experiments using a HRMAS probe to investigate the supramolecular arrangement within these gels. A comparison of 1H–1H NOESY spectra in solution and in the gel phase confirmed that the organizational structure observed in solution is maintained in the gel phase (Fig. S15–17 in the ESI†).
This intriguing observation prompted us to examine the gelation capabilities of (Z)-4 and (E)-4 through H-bonding in various nonpolar or apolar solvents, as well as in aprotic or protic solvents. The solvents included dodecane, heptane, cyclohexane, toluene, dioxane, and octan-1-ol. In the gelation test, the solution was heated until the compound was completely dissolved, and the resulting clear solution was then cooled to room temperature. Gel formation was determined using the tube inversion method, which involves verifying whether the solution flows when the tube is turned upside down (Fig. 5a). Hydantoins 4 demonstrated the ability to form gels in a range of nonpolar solvents, such as dodecane, heptane, and cyclohexane, while they remained soluble in toluene, dioxane, and octan-1-ol. In all the cases, gel formation was immediate, resulting in transparent gels that remained stable at room temperature. To obtain information on the structures of these gels, the xerogels obtained in the three solvents at concentrations of 1 and 2% were studied by transmission electron microscopy (TEM) and scanning electron microscopy (SEM) (see Fig. 5b and c and Tables S2 and S3 in section 3.5 in the ESI†). Notably, the lowest critical gelling concentration was observed with cyclohexane, at just 0.5 wt%, indicating that it acted as a supergelator (Fig. 5b and c). TEM images revealed that the xerogels of (Z)-4 and (E)-4 consisted of a network of entangled fibers, of several tens of nanometers, capable of encapsulating the solvent. SEM observations yielded similar findings, confirming the presence of a fibrillar network, consistent with that observed in TEM (see Table S3 in the ESI†). It is worth noting that in both cases, the xerogel tended to collapse into a thin film due to solvent removal during the sample preparation.
To gain a more profound understanding of the molecular architecture within the fibers, we conducted X-ray diffraction (XRD) experiments (Fig. S18 in the ESI†). The XRD patterns obtained from the organogels displayed two sharp maxima in the small-angle region with d-spacings in the ratio of 1:1√3 arising from the (1 0 0) and (1 1 0) reflections of a two-dimensional hexagonal lattice in a columnar hexagonal arrangement. The calculated lattice parameter (a) of the organogels prepared in cyclohexane is 48.1 Å. These materials also displayed diffuse scattering in the high-angle region (around 4.4 Å) that corresponds to the aliphatic chains with conformational disorder and the solvent. Remarkably, these XRD patterns closely resemble those of the liquid crystal phase (vide infra), suggesting that the supramolecular arrangement along the fibers of the organogels corresponds to a columnar mesomorphic-like organization.
T 2% (°C) | Thermal transitionsb | XRD parametersc | |
---|---|---|---|
a Temperature at which 2% mass loss is detected in the thermogravimetric curve. b DSC data of the second heating and cooling scans at a rate of 10 °C min−1. Temperatures (°C) are read at the maximum of the corresponding peaks, and enthalpies (kJ mol−1) are given in brackets. I: isotropic liquid, Colh: hexagonal columnar mesophase, Cr: crystal. c a: lattice constant of the Colh phase (Å), h: estimated mean stacking distance (Å), Z: molecules per disc. | |||
(Z)-4 | 312 | Colh 112 (4.3) I | a = 48.4 Å |
I 77 (5.0) Colh | h = 4.1 Å | ||
Z = 6 | |||
(E)-4 | 285 | Cr1 51 (2.8) Cr2 122 (38.5) Colh 146 (0.9) I | a = 49.5 Å |
I 142 (1.0) Colh 74 (33.2) Cr2 15 (3.4) Cr1 | h = 3.9 Å | ||
Z = 6 |
The DSC traces of (Z)-4 exhibited a single peak corresponding to the isotropization temperature (Fig. 6a). This transition was represented by as a complex and broad peak, which appeared to result from the overlap of two distinct transitions and was likely associated with a mesophase-to-mesophase transformation. It is worth noting that this transition was not observed through POM analysis (further details below), and as a result, it was treated and integrated as a single transition. In contrast, (E)-4 displayed distinct thermal characteristics. It exhibited a notably higher isotropization temperature compared to (Z)-4. Furthermore, (E)-4 demonstrated a higher propensity to crystallize and exhibited a crystal-to-crystal transition at lower temperatures. The thermal properties were also studied by POM (Fig. 6b). (Z)-4 appeared as a waxy solid at room temperature. Following the first heating–cooling cycle, it exhibited a grained texture, characterized by a mosaic texture with small domains. On the other hand, (E)-4 initially appeared as a waxy substance containing microcrystals and displayed higher birefringence than (Z)-4. After the first heating–cooling cycle, it exhibited a well-defined schlieren mosaic texture.
The absolute assignment of the liquid crystal phases was performed by powder X ray diffraction (XRD). Each sample was heated up to the isotropic liquid and then cooled down to room temperature. The two compounds showed similar X ray patterns with a broad, diffuse halo in the high angle region that corresponds to an average distance of 4.5 Å (Fig. 6c). This halo is usually found in liquid crystal phases and is assigned to the short range correlations between the conformationally disordered hydrocarbon chains. In the low angle region, the patterns exhibited a set of four sharp maxima in a reciprocal ratio of 1:√3:√4:√7. This ratio is characteristic of a well-organized hexagonal columnar phase (Colh) and corresponds to the (1 0), (1 1), (2 0) and (2 1) reflections from the two dimensional hexagonal lattice. The distances measured by XRD are listed in Table S1 of the ESI,† and the calculated lattice constants (a) are displayed in Table 1.
These liquid crystal hydantoins do not possess the traditional discotic shape, typically required for the formation of columnar liquid crystal phases. However, they can engage in the formation of H-bonded aggregates that adopt a discotic shape, enabling further self-assembly into columns. To gain a deeper understanding of the molecular packing within these columns, we conducted density calculations.55–58
Taking into consideration the XRD parameters and assuming a density of approximately 1 g cm−3, a typical value for liquid crystals, we estimated that the number of molecules occupying the cross-sectional area of the column (denoted as Z) was six, giving mean stacking distances ranging from 3.9 to 4.1 Å (see Table 1). Based on these findings, we propose a model for the molecular arrangement within a columnar stratum, wherein six hydantoin molecules self-assemble through H-bonding to form a rosette-like supramolecular disc. This particular arrangement, favored by intermolecular H-bonding, facilitates efficient space filling within the column cross-section. Despite the notable differences in the thermal behavior between the two diastereomers, it is worth highlighting that they exhibited similar lattice parameters, indicative of a similar supramolecular organization.
In order to facilitate the theoretical calculation, the geometry (only of the most stable conformers NHa φe and NHe φe) of diastereomers (Z)-4 and (E)-4 was simplified by replacing the R group (Scheme 1) with a hydrogen atom. The resulting model compounds were named (Z)-4 M and (E)-4 M (Fig. 7). These models were optimized at the B3LYP D3/6 31G** level (see computational details in the ESI†). Hydrogen bonded dimers of optimized (Z)-4 M and (E)-4 M were optimized at the same level, and the calculated structures are shown in Fig. 7. The hydantoin rings of the optimized (E)-4 M dimer are coplanar, but, unexpectedly, those of (Z)-4 M show a dihedral angle of 28 degrees. The distances of hydrogen bonds in the coplanar dimer of (E)-4 M are 1.803 and 1.798 Å, while in the (Z)-4 M dimer, the H-bond distances are 1.881 and 1.782 Å (Fig. S22 in the ESI†). Despite the distance of hydrogen interactions being smaller on average in the (E)-4 M dimer, the interaction energy of the non-planar (Z)-4 M dimer is −77.0 kJ mol−1 and that of the coplanar (E)-4 M dimer is −66.3 kJ mol−1 (Table S4 in the ESI†). These results seem to suggest that the (Z)-4 M dimer is supported by two hydrogen bonds and, probably, by other non-covalent interactions that enhance the stability of this non-planar dimer.
To gain a deeper understanding of these results, the NCI analysis was performed.59Fig. 7 displays the three-dimensional image of this analysis for (E)-4 M, (Z)-4 M and their corresponding dimers (for NCI plots, see Fig. S23 in the ESI†). The first noteworthy result of the NCI analysis of the monomers is that the dispersive interactions (green isosurface in Fig. 7) are established between the hydrogen atoms of the cyclohexane ring and the hydantoin moiety. In fact, these interactions are significantly stronger in the (E)-4 M isomer, where they involved the oxygen atom. The NCI analysis also revealed dispersive interactions between the hydrogen atoms of the cyclohexane ring and the aromatic hydrogens, which are similar in both isomers. As expected, the NCI analysis of the dimers shows the hydrogen bonds as blue disc-shaped isosurfaces, characteristic of highly localized attractive interactions. The analysis of dispersive interactions in the dimers reveals the second remarkable result of the NCI analysis: in the coplanar dimer, (E)-4 M, the hydrogen bonds are the only interactions established between the two molecules, while in the non-planar dimer, (Z)-4 M, in addition to the hydrogen bonds, extended dispersive interactions are established throughout the length of the molecules. In summary, the (Z)-4 M dimer achieves stabilization by adopting a non-planar conformation that, while maintaining the hydrogen bonds, allows for the extension of dispersive interactions. However, the (E)-4 M dimer prefers to maintain the intramolecular dispersive interactions and relies on the hydrogen bonds alone to support its structure. These results justify the differences in the interaction energy between both dimers.
From the XRD data, it has been established that the liquid crystal phases of (Z)-4 and (E)-4 are hexagonal columnar. The columns of these phases are formed by the stacking of complexes consisting of six hydantoin molecules supported, mainly, by intermolecular H bonding. With this information and considering the results of the theoretical calculations, we built up hexamers of both diastereomers, (Z)-4 M and (E)-4 M, combining three optimized dimers. These complexes were optimized at the B3LYP D3/6 31G** level. Fig. 8 shows the structure of the optimized hexamers and Table S4 in the ESI† shows their interaction energies. The optimized (Z)-4 M hexamer exhibits a bowl-like geometry with C6 symmetry, while the optimized (E)-4 M hexamer is a disc-like structure with C3 symmetry. The NCI analysis of the hexamers confirmed that a bowl-like structure allows, in addition to hydrogen bonding, the dispersive interactions between the molecules of (Z)-4 M, while only hydrogen bonds support the coplanar disc-like structure of the (E)-4 M hexamer (Fig. 8). Table S4 in the ESI† shows the total interaction energies of optimized dimers and hexamers. As expected, the interaction energy of the (Z)-4-M hexamer (−444.9 kJ mol−1) is lower than that of the (E)-4-M hexamer (−317.5 kJ mol−1), which confirms more intense intermolecular interactions in the (Z)-4-M hexamer.
The application of the results from the theoretical studies of the (Z)-4-M and (E)-4-M models to the (Z)-4 and (E)-4 isomers provide very valuable information to explain the different properties of both diastereomers. These studies confirm a different structure of the rosette formed in the supramolecular organization of the two isomers. Thus, the bowl-like structure of the (Z)-4 hexamer, although favors the columnar stacking, seems to hinder the non-covalent interactions between the staked bowls. Consequently, a vitrification process is observed in the cooling process, affording a non-crystalline waxy material. In addition to this, after passing into an isotropic liquid this compound exhibits an important hysteresis in the appearance of the mesophase in the cooling process, probably due to the important reorganization of the molecules necessary to form the supramolecular structure. On the other hand, a disk-like planar structure appears as the most stable in the hexameric structure of the minority (E)-4 isomer. This structure favors the interaction between the discs, increasing the π–π interactions and favoring the formation of the crystalline state in the cooling process. Also, in this compound the transition temperatures among different phases are clearly defined in both, the heating and the cooling process, and do not present hysteresis in this last process. Finally, the bowl-like structure of the (Z)-4 M hexamer as opposed to the disc-like geometry of the (E)-4 M hexamer might well be the reason for the slightly-larger stacking distance estimated for the mesophase of (Z)-4 (h = 4.1 Å, Table 1), compared to the value estimated for the disc-like (E)-4 M hexamer (h = 3.9 Å, Table 1).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3qo01832f |
This journal is © the Partner Organisations 2024 |