Open Access Article
Ho-Yin
Lai
a,
Yao-Chih
Lu
a,
Hsiang-Jen
Cheng
b,
Hsiu-Fu
Hsu
*b and
Long-Li
Lai
*a
aDepartment of Applied Chemistry, National Chi Nan University, Puli, Nantou 545, Taiwan. E-mail: lilai@ncnu.edu.tw
bDepartment of Chemistry, Tamkang University, Tamsui, 251, Taiwan
First published on 21st November 2025
Four macrocycles with the penta-triazine and penta-piperazine scaffold were efficiently synthesized via self-cyclization in 15–23% yields; the synthesized macrocycles were assumed to have crown-shaped conformations by NMR spectroscopy and computational simulations and contain void spaces in the bulky state by CO2 isotherms.
On the other hand, metal–organic frameworks (MOFs),10 covalent–organic frameworks (COFs),11 and porous organic polymers (POPs)12 have been extensively investigated as representative porous materials. However, their recrystallization or purification after use is still of concern due to their poor solubility in common solvents.13 Although hydrogen-bonded organic frameworks (HOFs) do not suffer from this solubility issue, they are generally unstable, and their cavities often collapse or shrink significantly upon solvent evaporation.9f,14
Thus, the development of macrocycles that exhibit good solubility in organic solvents and also function as porous materials is highly desirable. Such macrocycles not only enable efficient purification after synthesis but also facilitate reprocessing after usage. Moreover, the intrinsic cavities within these macrocycles can serve as hosts for guest molecules in the solid and solution states.15 In this work, we first synthesized compounds 2a–2d, which underwent self-cyclization to afford the corresponding macrocycles 3a–3d with the penta-triazine and penta-piperazine [5T+5P] scaffold in 15–23% yields. Macrocycle 3b, as a representative member of this series, was further investigated in the solid state by CO2 adsorption isotherms. We now report these results.
Compounds 2a–2d, serving as precursors to the target macrocycles, were efficiently synthesized following the reported procedure (Scheme 1).4b,c,9h Macrocycles 3a–3c were subsequently obtained from 2a–2c, respectively, by self-cyclization at 70 °C. In contrast, a higher temperature of 110 °C was required to completely convert 2d into 3d. It is noteworthy that compounds 2a–2d, which contain both electron-donating and electron-withdrawing substituents, were found to undergo slight decomposition either in solution or in the solid state and were therefore employed directly in the subsequent cyclization immediately after preparation.
As mentioned previously, macrocycles featuring a tri-triazine/tri-piperazine [3T+3P] scaffold have been prepared via combinatorial synthesis as potential molecular receptors.9h This observation prompted us to target the synthesis of triazine- and piperazine-based macrocycles with larger [5T+5P] scaffolds, which possess significantly larger cavities and are expected to exhibit more diverse host–guest recognition capabilities. However, the stepwise synthetic route employing protection–deprotection strategies reported in the literature for the preparation of [3T+3P] macrocycles9h is not well suited for the synthesis of macrocycles 3a–3d. The approach to construct the larger macrocycles via self-cyclization from precursors 2a–2d should greatly simplify the synthesis as shown in Scheme 1.
The multistep procedures containing protection and deprotection are rather laborious and may give the targeted macrocycles in very low yields. Additionally, they may be costly because lots of solvents and silica may be used in purification. In addition to the stepwise [3T+3P] systems,9h another macrocyclic framework containing triazine and piperazine units with a [6T+6P] scaffold has been reported as an intermediate in the one-pot synthesis of covalent organic polymers (COPs).16 However, this hexagonal macrocyclic scaffold, either presented as an intermediate or within COPs, was not spectroscopically confirmed but was merely proposed as the thermodynamically stable structure. Instead of the [6T+6P] scaffold,16 the [5T+5P] scaffold was obtained for 3a–3d in our synthesis, as directly evidenced by high-resolution mass spectrometry (HRMS), and further indirectly confirmed by nuclear magnetic resonance (NMR) analyses. HRMS analysis of 3a revealed a molecular ion peak at m/z 1452.1137 with an isotope distribution matching the theoretically calculated pattern for the cyclic [5T+5P] framework (Fig. S1). The HRMS spectra of 3b–3d are also provided in Fig. S1.
The 1H NMR spectrum of 3a (Fig. 1), which displays a single resonance for each type of proton, indicates a highly symmetrical molecular structure. Similarly, the 13C NMR spectrum (Fig. S2) exhibits signals consistent with a symmetric framework, in agreement with the 1H NMR observations. These spectroscopic results collectively support the cyclic structure of 3a–3d, which should have a [5T+5P] scaffold further on the analysis of HRMS.
The conformation of the macrocycle was further investigated by 1H NMR spectroscopy. The axial–equatorial interconversion of the CH2 units and chair–boat interconversion of the piperazine in planar macrocycles have been reported to be rapid even at temperatures as low as 203 K.17 At room temperature, the two distinct resonances observed at 3.79 and 3.75 ppm for the piperazine CH2 protons (HA and HB) do not arise from the axial–equatorial and chair–boat interconversion, but indicate that these protons reside in non-equivalent chemical environments, implying the nonplanar conformation of macrocycle 3a. A high energy barrier for inversion of the crown conformation18 is expected to result in two inequivalent piperazine CH2 environments, with HA located inside the crown and HB positioned outside (Fig. 1).
Variable-temperature 1H NMR spectra of 3a recorded between 20 and –60 °C in the 0.6–4.4 ppm region (Fig. S3) also showed that the axial–equatorial and chair–boat interconversion of the piperazine CH2 units in a non-planar macrocycle is also rapid, as that observed in a planar macrocycle.17 This further confirms that the resonances at 3.79 and 3.75 ppm actually arise from non-equivalent chemical environments of Hs. In the 13C NMR spectrum collected at room temperature (Fig. S2), two distinct signals were also observed for the piperazine CH2 carbons at 43.18 and 43.34 ppm (CA and CB), further supporting the presence of two chemically non-equivalent CH2 environments. These findings are consistent with the optimized crown-shaped conformation and corroborate the conclusions drawn from 1H NMR analysis.
To gain structural insights, macrocycle 3a was subjected to theoretical calculations for geometry optimization and heat of formation. Calculations were performed using the CaChe program with the MM2 force field19 and density functional theory (DFT) in the gas phase. In addition to 3a with a [5T+5P] scaffold, the analogous macrocycle with a [6T+6P] scaffold was also simulated (Fig. 2). Both MM2 and DFT calculations converged to a crown-shaped conformation for the [5T+5P] macrocycle and a nearly planar conformation for the [6T+6P] macrocycle (Fig. S4 and S5).
For macrocycles bearing identical peripheral groups, R = N(C4H9)2, the calculated formation energies of the [5T+5P] and [6T+6P] scaffolds were 0.108 and 0.107 kcal·mmol−1·atom−1, respectively. The corresponding total energies were −12
112 and −12
218 kcal·mmol−1·atom−1 using the B3LYP functional with the 6-311+G** basis set (Table S1).20 Both the MM2 and DFT results are in good agreement, indicating that the [6T+6P] macrocycle is only slightly more thermodynamically stable than the [5T+5P] macrocycle. These results suggest that the preferential formation of [5T+5P] macrocycles is governed by kinetic rather than thermodynamic control.
Furthermore, conformers of 3a with varying ratios of chair/boat forms in the piperazine units were also calculated, with results summarized in Tables S2, S3 and Fig. S6. The calculations revealed that each conversion from a chair to a boat conformation required approximately 2.7 kcal·mol−1. Among all eight possible chair/boat conformers, the fully chair conformation 3a with crown shape was found to be the most thermodynamically stable.
The 1H and 13C NMR spectra of 3b–3d (Fig. S7 and S8) further support the nonplanar conformation of these macrocycles. For 3c and 3d, two distinct resonances were observed at approximately 3.79 and 3.75 ppm, consistent with the presence of two chemically non-equivalent piperazine CH2 environments. In the case of 3b, a broadened signal for the piperazine CH2 protons was detected near 3.79 ppm. This broadening can be attributed to a more dynamic conformation, likely arising from the reduced steric hindrance of its peripheral substituents.
The thermal stability of macrocycle 3a was evaluated by thermogravimetric analysis (TGA). Compound 3a exhibited thermal stability up to 280 °C (Fig. S9). The slight weight loss observed is attributed to solvent evaporation during heating, further supporting the presence of a porous framework in the bulky state. Similar TGA profiles were obtained for compounds 3b–3d (Fig. S9).
The porous properties of 3b, as a representative compound of this series, were further investigated by CO2 adsorption isotherms measured at 195, 273, and 298 K (Fig. 3a). The Brunauer–Emmett–Teller (BET) surface area21 of 3b was determined to be 154.2 m2·g−1 (Langmuir surface area: 193.0 m2·g−1) based on the CO2 isotherm at 195 K, which is comparable with that of a porous triazine-based macrocycle.9f The pore size distribution of 3b in the bulk state was calculated from the CO2 isotherm of 273 K using the DFT model (density functional theory model) according to the literature,22 which falls within the range of 4.5–8.5 Å (Fig. 3b). Based on the conformation of 3b without considering the peripheral substituents, the rim-to-rim distances in larger and smaller openings of crown-shape molecules are between 5.6 and 11.7 Å (Fig. 4), which indicates the rationality of pore size distribution from the CO2 isotherm of 273 K.
![]() | ||
| Fig. 4 The height and rim-to-rim distances in the upper and lower openings of the crown-shape molecule. | ||
From the CO2 isotherms collected at 273 and 298 K, the isosteric heat of adsorption (Qst) of 3b was calculated to be approximately 28.2 kJ·mol−1 at zero coverage and was observed to decrease progressively with increasing CO2 loading, as determined by the virial method (Fig. 3c).23 This trend can be rationalized by the preferential occupation of the strongest adsorption sites at low coverage, likely located near the triaminotriazine cores of 3b, where strong polar–polar interactions occur between CO2 molecules and triaminotriazine moieties.24 Consequently, Qst is initially high but decreases significantly as the nitrogen-containing sites become saturated with CO2.25
In summary, four triazine- and piperazine-based macrocycles with [5T+5P] scaffolds were successfully synthesized via a self-cyclization strategy, thereby eliminating the need for repetitive protection–deprotection steps. HRMS analysis confirmed the molecular compositions, while NMR spectroscopy and computational simulations indicated that macrocycles 3a–3d adopt crown-shaped conformations. The porosity of 3b was further demonstrated through surface area measurements and pore-size distribution analysis by CO2 adsorption. Collectively, these results establish a straightforward and efficient synthetic approach to access triazine–piperazine macrocycles with large internal cavities, which hold promise for future host–guest and adsorption-related applications.
| This journal is © The Royal Society of Chemistry 2026 |