Zhuoer
Wang
,
Zhaozhen
Cao
,
Aiyou
Hao
and
Pengyao
Xing
*
Key Laboratory of Colloid and Interface Chemistry of Ministry of Education and School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, People's Republic of China. E-mail: xingpengyao@sdu.edu.cn
First published on 5th April 2024
Pnictogen bonding (PnB) is an attraction interaction that originates from the anisotropic distribution of electron density of pnictogen elements, which however has been rarely found in nitrogen atoms. In this work, for the first time, we unveil the general presence of N-involved PnB in aromatic or aliphatic imide groups and reveal its implications in chiral self-assembly of folding. This long-neglected interaction was consolidated by Cambridge structural database (CSD) searching as well as subsequent computational studies. Though the presence of PnB has limited effects on spectroscopic properties in the solution phase, conformation locking effects are sufficiently expressed in the chiral folding and self-assembly behavior. PnB anchors the chiral conformation to control the emergence and inversion of chiroptical signals, while intramolecular PnB induces the formation of supramolecular tilt chirality. It also enables the chiral folding of imide-containing amino acid or peptide derivatives, which induces the formation of unique secondary structural sequences such as β-sheets. Finally, the effects of PnB in directing folded helical structures were revealed. Examples of cysteine and cystine derivatives containing multiple N⋯O and N⋯S PnBs constitute an α-helix like secondary structure with characteristic circular dichroism. This work discloses the comprehensive existence of imide-involved PnB, illustrates its important role in folding and self-assembly, and sheds light on the rational fabrication of conformation-locked compounds and polymers with controllable chiroptical activities.
σ-hole (or π-hole) based interactions have implications on the self-assembly, which have shown great impacts on the properties of self-assembled arrays, phase behaviors and supramolecular chirality.19–22 Inspired by σ-holes in some important biomacromolecules, exploration of weak force driven chiral self-assembly provides a novel methodology for functional supramolecular aggregated or folded chiral molecules.23–25 Developing new types of weak forces and expanding rarely explored elements as donors would initiate improved functionality, which however remains considerably challenging. Nitrogen widely exists in amino acids, peptides and proteins, which are in the forms of amides and amines.26–28 Though in the amide form, the appended carbonyl groups would withdraw electrons from nitrogen, it could hardly induce the emergence of practical σ-holes, and thus no amide-based PnBs are found. Bickelhaupt et al. obtained the interaction energy (ΔEint) values of halogenated nitrogen derivatives with halides, suggesting the ΔEint of X3N⋯X- (X = F, Cl, Br) ranging from −3.5 to −86.1 kcal mol−1.29 This work provides theoretical hints that by appropriately controlling electron-withdrawing groups, intermolecular N-involved PnB might be realized in the solution or solid phase.
On top of the above developments, we considered that imides consisting of two acyl groups bound to nitrogen would exert a great electron-withdrawing effect, which might induce the emergence of positive hole regions. Imides were widely used in biochemical agents, pharmacy, industrial chemical intermediates and some high strength or conductive polymers such as Kapton .30–34 Imide containing compounds such as succinimide, maleimide, glutarimide, phthalimide and naphthalimide greatly facilitate the fabrication of functional compounds such as the famous perylene diimide (PDI)-based fluorophores.35 Arylene imides and diimides are also widely employed as electron-rich and electron-deficient scaffolds.36 Nevertheless, the potential imide N-involved PnB has seldom been referred to, which although may have great impacts on the folding and self-assembly of imide-containing compounds and biomacromolecules (Scheme 1). The imide-based PnB interaction can be found in some proteins (Scheme 1a), which stabilize structural formation and stabilization, as well as their biofunctions.37 Through this inspiration, in this work, we for the first time unveil the general presence of intra- and intermolecular PnB in imide-based compounds where imide N atoms act as donors. CSD searching indicates that there are thousands of X-ray structures containing imide-based PnB, most of which occur intramolecularly. The occurrence of intramolecular PnB is in most cases aided by the short tether (less than 3 atoms) between PnB acceptors (N, O, and S) and imide N. And the N⋯O distances (for instance) are within 2.5∼3.0 Å, shorter than the sum of the van der Waals radii (3.07 Å). The PnB locks the conformation by anchoring the side substitutes, which, in the chiral compounds, enables the fixation of chiral conformation and subsequently allows for chirality transfer with controlled chiroptical activities (Scheme 1b and c). The implications for chiral properties at the molecular and supramolecular level were illustrated in intra- and intermolecular PnB by highlighting several representative examples, especially in the simple folding of amino acid or peptide imide derivatives and the formation of supramolecular tilt chirality initiated by intermolecular PnB (Scheme 1d). This work systematically unveils the occurrence, properties and applications of a new type of PnB in chiral folding and self-assembly and sheds light on the rational design of imide-based functional small organic compounds or polymers.
Representatively, we screened some typical compounds with imide groups and potential PnB shown in Chart 1. Except for 3, other compounds possess imides conjugated with aromatic groups. Most of these compounds show short contacts between O and N either in an intramolecular or intermolecular manner (Fig. 2a and S10†). Phthalimide derivatives 1 and 2 whose alcohol moieties are conjugated to imide N via two carbon atom tethers feature N⋯O short contacts. The distances are determined to be 2.871 and 2.793 Å respectively. Structural extension to naphthalimide groups leads to such a short contact as well (compound 6, d = 2.946 Å). It also applies for other types of aromatic imide structures of compounds 5 and 8, giving N⋯O contacts of 2.709 and 2.750 Å respectively. Nonaromatic succinimide derivative 3 contains a similar N⋯O short contact (2.772 Å), indicating that the aromatic imide structure is not a necessity. Electrostatic potential (ESP) maps were calculated based on the X-ray structures of 6 and 7 (Fig. 2b, S11 and Table S5†).44 On both sides of imide N, two most positive (VS, max) regions were found with values ranging from 7.3 to 13.9 kcal mol−1, which might be related to the strong electron-withdrawing effect of the carbonyl group that induces the emergence of π-holes. It is quite similar to the π-holes on nitro groups. The uneven distribution of VS, max on both sides of N is possibly attributed to the neighbouring groups. And in both cases, O points to the regions with the highest VS, max values (12.7 and 13.9 kcal mol−1 respectively for 6 and 7). Noncovalent interaction (NCI) analysis was performed on compounds 4, 7 and 9 respectively. Compound 4 features intermolecular stacking supported by the N⋯O short contact, which corresponds to the van der Waals forces (vdWs, green region). In an intramolecular manner, 7 and 9 comprise attraction force between O and N regions. NCI analysis verifies the crucial role of N⋯O interaction in directing the self-assembly and folding of imide structures. By rotating the dihedral angle θ and freezing other parts, the changes in electronic energy are summarized in Fig. 2d. The lowest energy (θ ≈ 90°) corresponds to the formation of a N⋯O short contact, while the highest energy (θ ≈ 0°) is contributed by the repulsion between amide O and hydroxyl groups. The energy gap is determined to be 35.3 kcal mol−1, which effectively hampers the free rotation of the N–C bond. This energy is contributed by the O/O repulsion and N⋯O attraction, which in a synergistic manner anchor the very conformation. The rotary transition state (TS) of compound 7 was further calculated at the B3LYP-D3/def2SVP level of theory, which afforded a rotary energy barrier of 15.66 kcal mol−1, lower than the extremum energy gap. NBO calculations on the compounds 1, 9 and 10 were performed. However, the values are too small to be observed, which means that the orbital–orbital components in the N–O type PnB are negligible. We also checked the possibility of n to π* type of interactions and the Burgi–Dunitz angle from such interactions, as shown in Table S6.† It shows that most of the studied compounds do not contain n–π* interaction, while only compounds 4, 5 and 11 show n–π* interaction. These interactions synergistically contribute to the stability of PnB.
Will the weak N⋯O short contacts influence the properties in the solution and solid states? To figure out this question, we compared compounds 9 and 10 with respect to their folding and self-assembly arrangements as well as the spectroscopic data in the solution and solid phases (Fig. 3). Compounds 9 and 10 are isomerides whose O atoms are at different positions (Fig. 3a). This results in the formation of N⋯O PnB (2.957 Å) and H⋯O HB (2.046 Å) respectively for 9 and 10. In compound 10, the N⋯O distance of 3.632 Å is out of the sum of vdWs distances, and thus no PnB is expected. The different types of folding influence the self-assembled arrangements. The presence of Br shall introduce HaB that occurs between Br⋯O. Complementary Br⋯O HaB (3.184 Å) aids in the formation of box-like dimers which further stack through π–π stacking (P space group). The complementary HaB is apparently facilitated by the presence of PnB that draws the distance between Br and O. In comparison, compound 10 with HB would not favor complementary HaB, yet adopts a parallelly packed arrangement within a P21/c space group. The distinct packing in the solid state reflects the crucial role of PnB in directing self-assembly arrays. In the solution phase (CDCl3), chemical shifts of aromatic protons of 9 and 10 show differences (Fig. 3b) that protons of 9 shift to lower fields. As the tether between O and N contains two or more CH2, the inductive effect of O is negligible. The shifts might be associated with the N⋯O PnB and H⋯O HB. The nature of PnB is that O donates electrons to electrophilic regions of imide N, while for HB, imide O donates electrons to the electrophilic H. The above difference will slightly alter the electronic properties of the naphthalimide core both in solution and solid states. In the Fourier transform infrared spectroscopy (FT-IR) spectra, two characteristic absorption peaks at 1656 and 1342 cm−1 respectively are assigned to the CO and C–N stretching vibrational bands, which shift to 1646 and 1331 cm−1 ascribed to the formation of PnB (Fig. 3c). In the diluted solution, we also observed slight emission and absorption spectroscopic changes (Fig. 3d). Compound 10 shows bathochromic shifts with respect to the absorption and emission spectra ascribed to the distinct properties of PnB and HB. The effect was more pronounced in the solid-state diffusion reflectance spectra (Fig. 3e). Compound 10 exhibits a more profound bathochromic shift with a newly emerged peak at around 450 nm. This means that the design of PnB may realize chromic properties. The heavy atomic effects of Br enable room temperature phosphorescence in the solid state (Fig. 3f). Compounds 9 and 10 show different decay rates that afford a lifetime of 6.898 ms and 5.398 ms respectively. And in the corresponding delayed emission spectra, 9 also gives rise to a slight hypsochromic shift compared to 10 as well. Through the above comparison, the differences are less pronounced for 9 and 10 due to the weak nature of the PnB. Nevertheless, such a difference is potentially boosted in future design such as introducing PnB in polymeric materials, which may significantly influence self-assembly, photophysical properties and performances.
Then the chiroptical properties were probed using circular dichroism (CD) spectra (Fig. 4a, S13–S18†). In the solution state, compounds 1 and 2 with similar chromophores and absolute chirality yet with different substitutes yielded active Cotton effects (Fig. 4a and b). The active Cotton effects range from 230 to 300 nm and fall into the absorption of phthalimides, which are not intrinsically chiral. The conformation-locked structure by N⋯O PnB would facilitate chirality transfer to aromatic entities from amino alcohols. L- and D-enantiomers afforded positive and negative CD signs respectively, and this principle is also observed in an analogue shown in Fig. 4c. We also explored the chiroptical properties of 6, 7 and 8 in the solution and solid state (Fig. 4d–f). All three compounds share a similar naphthalimide aromatic core conjugated to the chiral entities. Using the nanoprecipitation method, after pouring the DMSO stock solutions into poor solvent water, self-assemblies were obtained. As shown in Fig. S19,† most CD spectra show great similarity compared to the solid-state CD spectra, indicating that the molecular arrangements in the solution-based self-assemblies follow a similar modality compared to the solid or crystalline states. Nevertheless, compound 7 exhibits an opposite sign to 6 though 6 and 7 possess identical bromo-naphthalimide chromophores (Fig. 4d and e). Their single crystal structures displayed in Fig. 4g suggest that the Br group has different relative orientation to the chiral parts. And the two relative orientations could be switched by rotating the N–C bond, which however is anchored by PnB and CO/O electrostatic repulsion forces (Fig. 2d). Calculated ECD spectra also verified the opposite mirror CD signal (Fig. S20†). The difference caused by the conformation is covered by the solid-state assembly where supramolecular packing leads to pronounced supramolecular chirality (7 and 8 exhibited a similar sign in the solid state). The above results illustrate the fact that N⋯O PnB could be used to anchor or manipulate the chiral conformation with active chiroptical properties. The high rotational barrier illustrated the integrity of PnB-locked conformation that might be resistant to temperature. Temperature-variable CD spectra of 6 and 7 in solution were collected (Fig. 4h and i). Heating from 0 to 70 °C produces the slightly declined CD signal. Subsequent cooling from 70 °C to 0 °C fully recovers the Cotton effect intensity and the corresponding absorption (Fig. S21†), suggesting that the chiral conformation has good resistance to temperature variations.
Then we explored the impact of PnB occurring intermolecularly in self-assemblies. Compounds 4 and 5 feature N⋯O PnB between imide CO and N respectively with distances of 2.979 and 3.040 Å (Fig. 5). The self-assembly is also aided by other supramolecular forces such as arene-perfluoroarene (AP) force between benzene and perfluorinated phthalimide moieties, which finally afforded one-dimensional packing arrays with specific homochirality (Fig. 5a and d). Within the one-dimensional growth, adjacent perfluorinated phthalimide moieties demonstrated a discontinuous tilt chirality with a designated 21 M-handedness. As the compound 4 or 5 has no inherent chirality, the self-resolution during crystallization can be regarded as the conglomeration behaviour (Fig. S22†).45 Conglomerate crystals tend to form chiral assemblies with unbiased M- or P-handedness. Consequently, the random measurement of different crystals or solids produced opposite handedness with mirror CD signs (Fig. 5b and e). Solid B and solid A respectively possessed positive and negative Cotton effects. To correlate the CD signal and handedness, ECD calculation was performed based on the tilt chirality found in the crystals. The fragmented M-helices of 4 and 5 gave exciton-type positive signals centered at around 350 nm (Fig. 5c and f). This strongly evidences the correlation between handedness and CD signs that solid A and solid B from 4 are P- and M-helices respectively. The observations imply the rational usage of N⋯O PnB in designing supramolecular chiral helices with predicted chiroptical properties.
The intramolecular or intermolecular N⋯O PnB favors the folding or self-assembly of amino acids or peptides into secondary structures. Some examples are shown in Fig. 6a, which are folded amino acid imide derivatives including Ala, Val, His, Gln, Phe and Leu residues. The N⋯O distances range from 2.599 to 2.774 Å. Due to the presence of carboxylic acid ester, there are two PnB acceptor O atoms, which did not show apparent selectivity. Adaptative to the structures of substituents, two O atoms form PnB in a nonpreferred selective manner. The conformation-locked geometries generally exist in the amino acid-containing imide derivatives, which shall facilitate the fabrication of chiroptical materials.
Intermolecular N⋯O PnBs between Ala–Ala dipeptides (compound 11) conjugated to aromatic imide structures were found to assist in the formation of supramolecular secondary structures. The N⋯O interaction between carbonyl and imide N crosslinks adjacent structures to afford supramolecular helices (Fig. 6b). By highlighting aromatic imide moieties, the 21 tilt helix was observed within a P212121 space group, and importantly, the dipeptide sequences adopt β-sheet arrangements. Interestingly, compound 11 does not have intramolecular N⋯O PnB in the X-ray structure, which correspondingly barely showed active Cotton effects (Fig. 6c) in the solution phase. This phenomenon further establishes the relationship between the active Cotton effects and PnB-locked conformations displayed in Fig. 4. In the solid state, an active CD signal of 11 was observed thanks to the presence of tilt chirality induced by PnB and other noncovalent forces.
Further implications of PnB on the intramolecular conformation locking or folding are given by the cystine and cysteine-based compounds 12 and 13 (Fig. 7a). Single and double functionalities of sulfur atoms perform as PnB acceptors, in addition to the copresence of ester oxygen atoms. Single crystal analysis indicated the presence of four PnB interactions that occurred intramolecularly for 12 and 13 respectively (Fig. 7b and c). The PnBs are continuous and adjacent. Compound 12 gives rise to CO‧‧‧N‧‧‧S‧‧‧N‧‧‧O arrays with distances of 2.816, 3.211, 3.211 and 2.846 Å respectively. Compound 13 gives rise to CO‧‧‧N‧‧‧S and S‧‧‧N‧‧‧OC arrays with distances of 2.764, 3.168, 3.062 and 2.649 Å respectively. Clearly, the S‧‧‧N distances are shorter than the sum of van der Waals radii of 332 pm. The conformation locking of cysteine and cystine derivatives through multiple PnBs is interesting and provides hints for the design of peptides or peptoids with special folding geometries. Highlighting the skeleton of compounds 12 and 13 witnessed the formation of P- and M-helical geometry respectively (Fig. 7b and c). It is not surprising because the absolute configuration of 12 and 13 is R,R and S,S respectively. Folding shall transfer the absolute chirality of cystine to the skeletons. We further synthesized compound 13 and measured the CD spectra (Fig. 7d). In diluted solution, RR13 and SS13 respectively gave negative and positive CD bands at around 250 nm, indicating the transfer of chirality to imide chromophores. This characteristic CD signal is verified using the calculated ECD spectra, showing similar curves at around 250 nm. Indeed, the derivation of either aromatic or aliphatic imide structures shall indiscriminately introduce PnB to enable new functionalities. Normally the formation of helical structures for peptides or peptoids requires at least four amino acid sequences with sufficient skeleton length and hydrogen bonds. The present cases shown by representatives 12 and 13 illustrate the feasible realization of helically folded structures using dimers.
The formation of intramolecular hydrogen bonding involving imide would show the shift of CO bands to lower energies, which shall appear at 1656 cm−1.46 As shown in Fig. S28,† the CO stretching bands of compounds 1, 2, 5 and 8 could be found at above 1700 cm−1, indicative of the failed formation of H-bonded pseudocycles. In contrast, other compounds 6, 7, 9, 10 and 11 demonstrate CO stretching bands at 1650, 1653, 1659, 1652 and 1668 cm−1 respectively, which verify the formation of intramolecular hydrogen bonding, accompanied by the potential formation of PnB. We also carried out temperature-variable 1H NMR tests of compounds 1 and 10 for comparison in dichloroethane-d4. As shown in 1H NMR spectra (Fig. S29†), increasing temperature from 298 K to 333 K led to the shifts of hydroxyl protons as well as the aromatic protons, which shift to high fields and lower fields respectively. Compound 1 mainly forms PnB without an intramolecular hydrogen bond, while compound 10 mainly forms a hydrogen bond. Increasing temperature would disassociate both the hydrogen bond and PnB, causing a 0.183 ppm shift of the OH group for compound 1. This value is larger (0.269 ppm) in compound 10 due to the breakdown of hydrogen bonds, and thus the shifts of aromatic protons (0.03 ppm) are more prominent than that of compound 1 (0.013 ppm). The different shift values verify the dominant formation of PnB and hydrogen bonds in 1 and 10 respectively. Based on the single crystal X-ray results, though the PnB and hydrogen bonds might coexist in solutions, they could hardly coexist in solid states ascribed to the highly directional angle of PnB or hydrogen bonds.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2309600, 2309601, 2309599, 2309598, 2309602. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc00554f |
This journal is © The Royal Society of Chemistry 2024 |