Ana M.
Ortuño
a,
Pablo
Reiné
a,
Sandra
Resa
a,
Luis
Álvarez de Cienfuegos
a,
Victor
Blanco
a,
José Manuel
Paredes
b,
Antonio J.
Mota
*c,
Giuseppe
Mazzeo
de,
Sergio
Abbate
de,
Jesus M.
Ugalde
f,
Vladimiro
Mujica
*fg,
Giovanna
Longhi
*de,
Delia
Miguel
*b and
Juan Manuel
Cuerva
*a
aDepartment of Organic Chemistry, Faculty of Science, Universidad de Granada (UGR), Unidad de Excelencia de Química (UEQ), E-18071 Granada, Spain. E-mail: jmcuerva@ugr.es
bDepartment of Physical Chemistry, Faculty of Pharmacy, UGR-UEQ. E-mail: dmalvarez@ugr.es
cDepartment of Inorganic Chemistry, Faculty of Science, UGR-UEQ. E-mail: mota@ugr.es
dDepartment of Molecular and Translational Medicine, Università di Brescia, Brescia, Italy
eIstituto Nazionale di Ottica – CNR, Brescia Research Unit, via Branze 45, 25123 Brescia, Italy. E-mail: giovanna.longhi@unibs.it
fKimika Fakultatea, Euskal Herriko Unibertsitatea and Donostia International Physics Center (DIPC), P. K. 1072, 20080 Donostia, Euskadi, Spain
gArizona State University, School of Molecular Sciences, Tempe, AZ 85287, USA. E-mail: vmujica@asu.edu
First published on 28th June 2021
A novel synthetic strategy based on a bifunctional stapled chiral nucleus from which segments of different lengths can be added to both ends of o-phenylene ethynylenes (o-OPEs) has been developed to obtain a new type of foldamer and a novel chiral Pd2L2 metallacycle. For the first time, an enantiopure fully conjugated helical foldamer having 14 phenyl rings and 13 alkynes is reported. The folded structure has four complete loops and is able to host three Ag(I) cations in their cavity with high binding constants. The complete photophysical and chiroptical (ECD, CPL and VCD) characterization of these foldamers has shown that these molecules show intense chiroptical responses with dissymmetry ratios in the range of 10−2. Theoretical modeling of these systems reveals the origin of these remarkable responses and points out a potential connection with the chiral induced spin selectivity (CISS) effect. The magnetic dipole moment is proposed as a key physical variable connecting the chiroptical properties and CISS-based spin filtering properties observed in chiral compounds.
Moreover, the remarkable geometry of [6]-helicenes has also been used in the context of chiral metallosupramolecular architectures12,13 as bidentate donors with an angle between the binding sites close to 0°.14 These systems have emerged as an unique research field owing to their electronic, optical and sensing capabilities. However, the main drawback in [n]-helicene chemistry is the challenging synthesis of members with increasing turns, even in their racemic form. In fact, the longest member of the series, [16]-helicene, composed of two and a half turns, represents an isolated example (Fig. 1, top).15
Although longer related helicene-like molecules (two and a half turns) have been described, easy access to fully conjugated helical systems continues to be a synthetic challenge.16,17 Therefore, easy-to-prepare extended and conjugated helical systems would be highly desirable for important applications like those mentioned above. Within this context, foldamers based on conjugated subunits constitute a privileged scaffold for a number of different reasons: (i) they can be easily prepared by iterative protocols using simple conjugated monomers with different geometries, (ii) the functionality present in such monomers can cooperatively act to strengthen their properties or even create new ones, and (iii) they can also interact dynamically with the environment. In fact, many aromatic foldamers can be switched between different states by different stimuli such as the presence of light,18–22 changes in temperature,23–26 solvent composition,27 anions,28–32 cations,33,34 or water35–38 binding, and redox processes.39,40 Moreover, foldamers are able to generate chirality in a dynamic way even using achiral monomers/repeating units. Therefore, the control of the P/M helicity41–46 yields systems that are able to act as chiroptical probes47,48 or even CPL emitters.49–53
In recent years, we have been working on conjugated helical structures based on o-OPEs54–58 as an alternative to [n]-helicenes exhibiting, in many cases, better chiroptical responses. Our approach relies on widening the inner cavity present in [n]-helicenes, which translates into an enhancement in the magnetic moment of optical transitions and the corresponding chiroptical response, especially at the longest wavelength. Thus, we have developed two strategies to obtain enantiopure folded fully conjugated helical systems: the first based on the use of chiral sulfoxides59,60 and, secondly, the use of chiral staples, demonstrating that enantiopure P/M helixes can be also prepared by efficient chirality transfer from a chiral staple to the o-OPE core, including even double staples to avoid racemization processes.61–64 However, the later systems could be synthesized with only one turn and a half.65
All the above approaches provide structures with remarkable chiroptical responses in ECD, CPL and VCD. It is worth noting that the response of ECD in terms of molar circular dichroism has been for some time correlated with the magnitude of the chiral-induced spin selectivity (CISS) effect.66 Actually, this phenomenon has been recently discovered, offering a promising change of paradigm for the discrimination and filtering of spin states. The spin of the electron, S = ½, yields two quantum states characterized by its projection on a given axis, ms = ±½. The two states are degenerate in the absence of either internal or external magnetic fields and, consequently, any electron current under such circumstances will consist of a half–half mixture of both states. The separation of the spin states of the electrons is a milestone of the greatest importance for the development of spintronics.
Although spin states of neutral particles67 can be selected passing through an in-homogeneous magnetic field, charged particles, like electrons, suffer from the Lorentz force arising from the interaction of the electric charge with the magnetic field, preventing a Zeeman's effect-based technique to be used for separating electron spin states. Although more complex architectures can be used,68 the (preferential) selection of the spin state of electrons travelling through chiral structures is highly appealing. The chemically relevant fact is that now it is the scattering of electrons through the molecular (helical) structure that induces the filtering of the electron spin states. Namely, depending on the helicity of the molecular structure, it preferentially allows the passage of one of the two spin states over the other, because helicity couples the spin of the electron with its own linear momentum, resulting in the selective transmission of electrons with the preferred spin state. This opens the way to great opportunities for “chemical control” in the selection and discrimination of the spin state of electric currents. The key concept resides in the achievement of precise chemical control on the factors affecting the molecular helical structure, i.e. length, diameter, pitch, chemical composition, chemical environment, either the solvent selection or its deposition on selected surfaces, including nanoparticles, etc. The absence of external electromagnetic fields avoids any interference with the spins in “uncontrollable” ways. The implications of the CISS effect for the control of the electronic spin state have already been documented.69 Nevertheless, a long-standing goal, not fully resolved at the time being, is to know beforehand which molecules are best for CISS spin-filtering. Within this context, remarkable effort has been made by Waldeck et al.,66 who showed conclusive experimental evidence for the direct correlation between CISS induced spin-filtering power and the magnitude of features of the molecule electronic circular dichroism (ECD) spectrum. A direct correlation of the dissymmetry of the electron-transfer rate for photoelectrons produced by circularly polarized light was found with both the sign and intensity of the ECD signal of electron acceptor chiral molecules. This suggests that the larger the ECD signal intensity the better the CISS spin-filtering a molecule can perform. The sign determines which spin is preferentially filtered.
Although the relationship between both phenomena seems plausible, it is not clear which factors are relevant for such a connection. One conceivable hypothesis is that the CISS effect may be associated with a transient magnetic moment in chiral molecules. Indeed, the intrinsic magnetic nature of the CISS effect suggests that induced magnetic transition moments (m) could be relevant to characterize the magnetic molecular response associated with the dynamic spin-polarized electron transport through a molecule. The optimal design of structures in which such a relationship could be validated is, therefore, essential.
Here we report on a versatile scaffold to explore new helical architectures using simple derivatizations. Chiral stapled nucleus (S,S,P)-1 can be used as a chiral seed for the synthesis of highly extended enantiopure conjugated helical systems (S,S,P)-2–(S,S,P)-6 (Fig. 1, bottom) containing one up to four complete turns. In the presence of Ag(I) we have obtained the most extended enantiopure systems with a well-defined structure described to date.70–73 This situation is optimal to gain insight into the importance of the length of helical systems, the magnitude of the magnetic dipole transition moment (m) associated with the less energetic transition, and their potential use as a predictor of CISS capabilities for an organic architecture. We must emphasize that a comprehensive theoretical model coupling optical dichroism with the CISS effect is still missing, but the work reported here may be a step forward in this direction.
Moreover, the usefulness of (S,S,P)-1 is also demonstrated with the preparation of ligand (S,S,P)-7. This compound increases the toolbox of the supramolecular research area showing that simple helical foldamers can be also used as a building block for the Pd(II)-directed self-assembly of a chiral Pd2L2 metallacycle, thus expanding the scarce number of examples of chiral metallosupramolecular architectures reported to date.13
Scheme 1 Synthesis of helical systems (S,S,P)-2–(S,S,P)-6 and ligand (S,S,P)-7. Reagents and conditions: (a) Ac2O (1.4 eq.), DMAP (2 eq.), CH2Cl2, rt, 5 min, 100%; (b) Pd(CH3CN)2Cl2 (0.03 eq.), PtBu3·HBF4 (0.06 eq.), CuI (0.03 eq.), iPr2NH, rt, 24 h, 70%; (c) K2CO3 (7 eq.), THF, MeOH, rt, 30 min, 88%; (d) CsCO3, CH3CN, reflux, 24 h, 56%; (e) Pd(CH3CN)2Cl2, CuI, PtBu3·HBF4, iPr2NH/THF, 60 °C, 24 h; (f) Pd(CH3CN)2Cl2, CuI, PtBu3·HBF4, iPr2NH/THF, 60 °C, 24 h; 81%. Yields of step e for each derivative are described in the ESI.† |
Circular dichroism of (S,S,P)-1 (CH2Cl2) showed a positive Cotton effect at the longest wavelength, Δε = +47 M−1 cm−1, which is characteristic of P helicity (Fig. 2b). We note that molar circular dichroism (Δε) depends on the intensity of absorption (that is, on electric transition dipole moment) and its magnitude cannot be always correlated with the intrinsic chirality of the sample. In that situation, a dimensionless gabs value, Δε/ε, is more representative and easy to compare with theoretical values. Nevertheless, extracting reliable experimental gabs values from ECD bands is tricky since calculations predict energetically close transitions with the opposite sign, such that the observed gabs value results from the algebraic sum of g-values for the individual transitions. Being interested in the first electronic transition, the best estimate can be obtained from the tail of the CD band. In that case, the experimental gabs value (+1 × 10−2) is comparable to previously described values for related systems.52,53
The origin of this helical chirality could be clarified using density functional theory (DFT) calculations. After a molecular mechanics (MM) conformational search, all structures within 5 kcal mol−1 have been optimized with both P and M helicities at the B3LYP/6-31G** level in CH2Cl2, treated with the polarizable continuum model (PCM). Only three conformers, corresponding to the relative spatial orientations of the bromine atoms (out–out, out–in and in–in), have been obtained with non-negligible populations (Fig. 2c and Table S10†). They all present P helicity and, considering the calculated energy, are nearly equally populated. As a consequence, the P epimer dominates the conformational equilibrium and makes possible a definite assignment of the (S,S,P)-1 absolute configuration in solution. This conclusion was strengthened by the similarity between the averaged simulated ECD and absorption spectra of the abovementioned P conformers (M06/6-31G** level) and the experimental ones (Fig. 2a and b). The calculated gabs value for the longest wavelength transition is in reasonable agreement with the experimental one.
Fortunately we could also obtain the single-crystal X-ray diffraction structure of (S,S,P)-1, which confirms its helical geometry with P preference in the solid state (Fig. 2c). This helical arrangement generates a dihedral angle between the central aromatic rings close to 40°. The outer aromatic rings adopt a parallel orientation, although displaced, with a distance between their main planes of 3.5 Å and a distance between centroids of 3.8 Å, suggesting the establishment of π-stacking interactions between them.
Being dynamic systems, foldamers are affected by the environment and, consequently, their properties are usually solvent dependent. In this case, we only observed a single set of signals in NMR, showing that any conformational equilibrium takes place faster than the NMR timescale. The absorbance and fluorescence of (S,S,P)-1 were analyzed in different solvents, including apolar, polar and protic ones. Although no significant differences were observed in the emission spectra, a slightly increment of the shoulder in the absorbance spectra was observed in some solvents (Fig. S109†). Fluorescence lifetimes were also measured at the emission maximum using an excitation source of 325 nm. In all tested solvents, decays were fit to a biexponential function where the shortest fluorescence lifetimes were around 0.5 ns in most cases, except in CH2Cl2 in which it presented a considerably lower value of 0.2 ns. However, we observed higher differences for the longest fluorescence lifetimes, which varied between 0.7 ns in CH2Cl2 to 4.1 ns in toluene (see Table S2† for more details). In view of these results we recorded the time-resolved emission spectra (TRES) of compound (S,S,P)-1 in CH2Cl2 covering all of the emission range, between 340 and 600 nm. This global analysis also afforded a biexponential decay showing, along the whole spectra, an exclusive contribution of the shortest lifetime species, the presence of the largest being almost negligible (see Fig. S106†). Regarding chiroptical properties and according to the abovementioned data, minor changes in the ECD of (S,S,P)-1 were observed. We noticed that the chiroptical response of (S,S,P)-1 is maximized in MeOH with an experimental gabs value of +1.3 × 10−2.
Since (S,S,P)-1 is fluorescent and chiral, CPL is a valuable tool to extract information of the dynamic system. Although the fluorescence was weak, reliable dissymmetry factors of the emission, glum values defined as 2(IL − IR)/(IL + IR), could be measured (CH2Cl2: glum = +1.1 × 10−2, MeOH: glum = +1.0 × 10−2, MeCN: glum = +1.0 × 10−2, n-hexane: glum = +0.9 × 10−2). These results suggest that the P helicity is preserved in the excited state. Consequently, the similarity of gabs and glum values also reveals that the structure is highly preserved in the excited state despite the apparent conformational flexibility.
Analogous studies were carried out with compounds (S,S)-2 to (S,S)-7 (Fig. 3). Again, chiroptical properties, especially ECD and CPL, were highly informative and solvent dependent. The compound presenting five alkynes, (S,S,P)-2, showed remarkable molar circular dichroism in MeCN (Δε = +41 M−1 cm−1, gabs = +1.0 × 10−2), being compatible with a relatively ordered structure. Again, we found a correlation between the polarity of the solvent and the gabs value for the longest wavelength transition, with similar values in MeOH (Δε = +36 M−1 cm−1, gabs = +1.4 × 10−2) to those obtained in MeCN, but being weaker in CH2Cl2 (Δε = +20 M−1 cm−1, gabs = +0.8 × 10−2) and hexane (Δε = +20 M−1 cm−1, gabs = +1 × 10−2). If we assume gabs values as the representative parameter of the system helicity, polar solvents are apparently favoring the folding process. Although this is a well-known effect for m-OPEs,27 this is the first time it is clearly observed for o-OPEs.
Theoretical calculations support the prevalence of the fully folded conformations both considering CH2Cl2 and acetonitrile in the polarizable continuum model using the integral equation formalism variant (iefpcm) (see the ESI†), showing again a clear preference for the P conformation.74 The rotational barrier for the unfolding process was theoretically estimated by turning one arm of the folded structure in 100 steps (Fig. S156†), achieving a maximum at 8.1 kcal mol−1. Therefore, the partially disordered situation favored by less polar solvents might explain the observed differences (see Fig. S119, in the ESI†).
Compound (S,S,P)-2 is also fluorescent, with quantum yield values varying from 19% to 29% depending on the solvent (see Table S3 in the ESI†). The fluorescence lifetime at the emission maximum showed again two values around 4 and 2 ns in all solvents. Interestingly, time-resolved fluorescence allows the deconvolution of the emission spectra, providing the ratio of the different emitting species by using the areas of the species-associated emission spectra (SAEMS). In the four solvents we could detect two main contributions coming from two different species with around 4.6 and 2.6 ns, which cannot interconvert in that timescale. Although relative proportion depends on the solvent nature (see Fig. S107†), two limit situations are present in hexane and MeOH, in which the corresponding populations are reversed. The coexistence of fully folded and partially folded structures in the excited state may be a reasonable explanation, fully folded structures being maximized in MeOH.
Since the chiroptical properties of the abovementioned species are expected to be different, CPL is once again a valuable tool to extract information on the dynamic system. Representative solvents as MeOH, MeCN, CH2Cl2 and hexane were evaluated using this technique, obtaining remarkably high glum values (MeOH: +1.3 × 10−2, MeCN: +1.3 × 10−2, CH2Cl2: +0.88 × 10−2, hexane: +1 × 10−2). Once more, a correlation with the polarity of the solvent was observed, the fully folded conformations with higher glum values being preferred in MeOH or MeCN. The positive sign of CPL is correlated with the positive sign of the ECD signal for the longest wavelength transition, P helicity being preserved in the excited state.4 The similarity of gabs and glum values also reveals that the geometrical structure of the ground state is preserved in the excited state. In fact, the optimized structure for the first excited state starting from the lowest-energy fully folded ground state structure resulted in a nearly superimposable helix in both ground and excited states, showing differences only in the bond length alteration (BLA) of the C–CC–C groups (see Fig. 4f and Table S12†).
Similar conclusions were extracted for (S,S,P)-3–(S,S,P)-5 and (S,S,P)-7. Molar circular dichroism and gabs values remain high, especially in MeCN ((S,S,P)-3: Δε = +21 M−1 cm−1, gabs = +0.9 × 10−2; (S,S,P)-4: Δε = +18 M−1 cm−1, gabs = +0.5 × 10−2; (S,S,P)-5: Δε = +47 M−1 cm−1, gabs = +0.9 × 10−2; (S,S,P)-7: Δε = +19 M−1 cm−1, gabs = +0.8 × 10−2), compared to CH2Cl2 ((S,S,P)-3: Δε = +6 M−1 cm−1, gabs = +0.3 × 10−2; (S,S,P)-4: Δε = +8 M−1 cm−1, gabs = +0.3 × 10−2; (S,S,P)-5: Δε = +13 M−1 cm−1, gabs = +0.3 × 10−2; (S,S,P)-7: Δε = +19 M−1 cm−1, gabs = +0.3 × 10−2), although exceptions can be also observed (see the ESI†). These higher gabs in MeCN suggest a relatively ordered structure in more polar solvents, in accordance with the stability of the helical structure. Nevertheless, the increasing number of potential conformations did not allow a complete modeling of the conformational space as in previous cases. Even so, it seems reasonable to assume that in the absence of solvophobic interactions many more disordered structures in solution could be present. For compound (S,S,P)-6, we found some solubility issues in hexane and MeOH and reliable CD spectra could be obtained only in CH2Cl2 (Δε = +13 M−1 cm−1, gabs = +0.3 × 10−2) and MeCN (Δε = +15 M−1 cm−1, gabs = +0.3 × 10−2).
Fluorescence lifetimes and quantum yields for such longer oligomers (S,S,P)-3-(S,S,P)-7 were evaluated in MeOH, MeCN, CH2Cl2 and hexane and are summarized in Table S4.† Quantum yields depend on both the solvent and number of alkynes present in the molecule, but no sensible relationships could be obtained in any of the cases. However, it is remarkable that, except for compound (S,S,P)-3, the intensity-averaged lifetime of longer oligomers is always lower in hexane (around 3 ns) than in MeOH and CH3CN (around 4 ns), showing again a higher contribution of the shortest lifetime species in apolar solvents, in accordance with the results obtained for TRES of the first member of the series (S,S,P)-2.
All helical systems were also very efficient CPL emitters, especially in MeOH ((S,S,P)-3, glum = +1.3 × 10−2, (S,S,P)-4, glum = +1.1 × 10−2, (S,S,P)-5, glum = +1.3 × 10−2, (S,S,P)-6, glum = +1.3 × 10−2, (S,S,P)-7, glum = +0.54 × 10−2). Modeling CPL properties has been be carried out for compound (S,S,P)-2 in the most stable folded conformer. However, the analysis of larger systems is computationally too demanding. Despite this limitation, the structures were optimized by DFT calculations (Fig. 4b–e). Although, the similarity between gabs and glum values, which is also observed in experimental trends, is expected, the glum value is controlled by the S1 → S0 transition, which cannot be assumed to be identical to the S0 → S1 transition because the structural changes in the excited state may affect the ground state transition parameters.
The first member of the series, (S,S,P)-1, presents evident structural changes after Ag(I) addition as we observed an increased resolution and downfield shielding of the aromatic protons (Fig. S49†). Regarding the 13C NMR spectra, alkyne carbons, which were initially located within 1 ppm to around 92 ppm, evolve to three different signals at 89, 92 and 94 ppm by coordination with the Ag(I) cation (Fig. S50†). The coordination process was also followed by ECD spectroscopy (Fig. S126†). In this particular case, we observed a decrease in the chiroptical response. This fact is in agreement with previous results in related o-OPEs with three alkynes owing to a planarization of the helical structure. ECD titration experiments resulted in a relatively strong single binding event, K1:Ag = 32480 ± 2% M−1. Theoretical calculations (M06/6-31G* plus LANL2DZ (valence + ECP) for silver in CH2Cl2) showed that the minimum energy structure corresponds to a distorted trigonal bipyramidal geometry (Fig. 6b) in which the bromine atoms participate in the coordination to the Ag(I) cation.
Fig. 5 Experimental (solid line) and calculated (dashed line) ECD spectra of: (a) (S,S,P)-2-(S,S,P)-3 and (b) (S,S,P)-4–(S,S,P)-6 in CH2Cl2 in the presence of an excess of Ag(I). |
Compound (S,S,P)-2 is expected to behave as a ligand for only one Ag(I) cation taking into account the previous studies of related o-OPE systems.61 The binding process can be again properly followed by 1H NMR spectroscopy. NMR titrations showed that the original signals of ortho hydrogen atoms of the inner phenyl rings experience the characteristic deshielding (i.e. from 6.45 ppm to 6.64 ppm, Fig. 6a) attributed to Ag(I) coordination. In addition, significant changes in the chemical shift of alkyne carbons were observed in the 13C NMR spectrum (Fig. S52†). The ECD response at the longest wavelength is intense (Δε = +85 M−1 cm−1, gabs = +2.5 × 10−2) in agreement with a well-defined helical structure upon the formation of the complex (Fig. 6). It is worth noting that such a dissymmetry ratio is huge compared with usual values for organic molecules at the longest wavelength. The corresponding ECD titration showed a high binding constant K2:Ag = 3.67 × 105 ± 5% M−1. As expected, the DFT-based conformational analysis (see the ESI†) of complex (S,S,P)-2:Ag is in agreement with experimental findings, showing a stabilization of the fully folded structure. The calculated ECD is also in good agreement with the experimental one, thus confirming the folding process with Ag(I) (Fig. 5a and 6c). In this case, the distorted trigonal bipyramidal geometry is generated by coordination with five alkynes. Ag–alkyne bond distances are predicted to lay within 2.51–2.69 Å, with the Ag–C (alkyne) distances ranging from 2.57 to 2.76 Å. The distance between the terminal aromatic rings of the side arms and those on the corners on the central core of the helix, which display a quite displaced parallel arrangement, is ca. 3.5–3.6 Å, which suggests some extent of π-interactions between them.
On the other hand, the increase to seven potential binding sites for Ag(I) cations in compound (S,S,P)-3 did not result in a significant increase of the chiroptical properties (Δε = +21 M−1 cm−1, gabs = +1 × 10−2). In principle, a single Ag(I) cation does not need all the alkynes in the binding event, making the remaining ones useless to stabilize a second Ag(I) cation. In this situation the presence of loosely or non-bound side arms can result in a more disordered structure compared to (S,S,P)-2, even in the presence of Ag(I). Unfortunately, the complex conformational space makes it difficult to perform any theoretical prediction. As a consequence of the previous assumption, ECD titrations showed that the compound is only bound to one Ag(I) cation with a binding constant K3:Ag = 3.7 × 106 ± 15% M−1. Minimum energy calculated structures showed that Ag(I) interacts with five alkynes. Only two alkynes remain free to bind a new cation. Our previous studies60 indicated that two alkynes are unable to provide an efficient coordination to another Ag(I) cation. The different spatial arrangement for these two additional unbound side chains and even in-cavity Ag(I) displacements could explain the weaker ECD response.
The presence of two additional alkynes in (S,S,P)-4 should be in principle beneficial for the coordination process. Nevertheless, the ECD response remains weak and similar to the previous case (Δε = +24 M−1 cm−1, gabs = +0.6 × 10−2). This fact suggest that at the concentrations used for the ECD experiments there is only a relevant binding constant, although subsequent Ag(I) coordination processes can take place with a modest binding constant. Accordingly, in the titration experiments a higher overall binding constant, K4:Ag2 = 2.4 × 1010 ± 12% M−2, was obtained (se Fig. S129–S131†). It is worth noting that the corresponding NMR studies were carried out at much higher concentrations. Under such experimental conditions more than two equivalents of the Ag(I) cation are required to fully resolve the signals, suggesting the existence of a symmetric dimetallated structure.
Compound (S,S,P)-5, presenting eleven alkynes, is able to accept two Ag(I) cations, thus yielding a longer helical structure. The first evidence came from the ECD spectra, showing a remarkable molar circular dichroism (Δε) at the longest wavelength (Δε = +85 M−1 cm−1, gabs = +1.7 × 10−2) in the presence of an excess of Ag(I). ECD titrations also showed strong coordination of two Ag(I) cations (K5:Ag2 = 18 × 109 ± 1.6% M−2). This is the first example of a pure and well-defined o-OPE system presenting an enantiopure helical structure with three and a half turns. This fact shows that the chiral nucleus can be useful with the aid of Ag(I) cations in the preparation of long and conjugated helical systems. Although different situations are possible, DFT calculations showed that the energetically favored case is that in which the Ag(I) cations are separated by a distance of approximately 5.2 Å (Fig. 6f).
The best chiroptical response was observed for compound (S,S,P)-6, which showed an even higher molar circular dichroism (Δε = +130 M−1 cm−1, gabs = +2.2 × 10−2). Again, this result strongly suggests a fully folded helical structure. ECD titrations also showed three strong binding events with three Ag(I) cations (K6:Ag3 = 2.7 × 1017 ± 19% M−3). As it can be seen in Fig. 6g, DFT calculations predict a highly symmetric structure with four complete loops. The interatomic distances between consecutive Ag(I) cations are 4.4 and 4.5 Å. Interestingly, on the basis of the Ag–alkyne distances, calculations predict that the metal centers are coordinated to five alkynes. As a result, there are two alkynes that coordinate simultaneously to two Ag(I) centers, which results in longer Ag–alkyne distances (2.74–2.89 Å) in comparison with those involving triple bonds coordinated just to one metal center (Ag–alkyne distances: 2.44–2.65 Å).
Fig. 7 Comparison of theoretical and experimental (a) IR and (b) VCD spectra of compound (S,S,P)-2 in CH2Cl2 solvent. |
To better understand the origin of the spectroscopic changes upon silver addition, we examined two representative cases with the aid of DFT calculations: the shortest oligomer (S,S,P)-2 and a longer case, (S,S,P)-5. For the first case, we considered the possible conformers and compared the calculated spectra for the most stable structures in the absence or presence of silver. The details of this analysis are given in the ESI (conformer characteristics and spectra).† Due to the staple, the structure of the short oligomer is rigid enough to present a VCD signal typical of the ordered OPE structures already observed.60,61 In particular, the configuration of the staple dictates the prevalent helicity sense and the most stable conformer presents already the correct shape to host a silver ion; however, also partially unfolded structures of (S,S,P)-2 maintain the spectroscopic pattern. In the presence of Ag(I), the calculated spectrum reproduces the wavenumber shift of the two higher energy components of the previously cited triplet (+,−,+) and the 1315 cm−1 feature intensification (Fig. 7). Also the IR spectra differences are correctly predicted.
In the case of the longer oligomer, Ag(I) cations give rise to the enhancement of the triplet at 1450–1500 cm−1. The phenomenon can be attributed to normal modes delocalized along the whole molecule, consisting of in-plane bending of the phenyl CH bonds. Structures with various Ag contents have been optimized: one central Ag ion, two Ag ions symmetrically or non-symmetrically disposed and three Ag ions (see Fig. S150† for the optimized structures and Fig. S151† for the corresponding calculated spectra). Overall, there is a good correspondence between calculated and experimental data. In fact, compound (S,S,P)-5 revealed to be quite efficient in incorporating silver ions and has sufficient length to show an enhancement of VCD signals when it is guided to form a regular helix in the presence of silver. In conclusion, in the presence of silver, the longer backbone assumes an ordered structure able to show intensification of the triplet; the increase in intensity is more regular considering VCD, while ECD shows a non-monotonic dependence on the alkyne number. A complete theoretical representation of all the situations and of the two spectroscopic responses (ECD and VCD) is not easy due to the complexity of the conformational landscape and of the evident solvent effect (see Fig. 2) which is not simple to model.
We have analyzed the relevant chiroptical parameters, such as electric and magnetic dipole transition moments |m| and |μ|, dipole (D) and rotational strengths (R), and angle between electric dipole transition moment and magnetic dipole transition moment, obtained by TD-DFT calculations (M06/6-31G**) of molecules (S,S,P)-2 to (S,S,P)-6 and their corresponding Ag(I) complexes. It is important to note that the studied systems are quite large and cannot be fully modelled at a high level of theory. Consequently, only tendencies and not precise values can be provided by the calculations. It is worth noting that theoretical calculations overestimate the predicted gabs values.77 In this case, we use CH2Cl2 as the solvent for a better comparison with the parent compounds and Ag(I) complexes derived.
Inspection of the data of Table 1 shows that the magnetic dipole transition moment for the S0 → S1 transition increases monotonically with the number of windings for the whole set of helical molecules (S,S,P)-2 to (S,S,P)-6. This is in agreement with the naïve view that holds that increasing the helicity of the compound should translate into an increase of the magnetic moment. It had been suggested that when the magnetic dipole transition moment stands for the power absorbed from the exciting light magnetic field in a conductive helix, a linear relationship between the transition magnetic dipole moment and the number of loops n exists.78
Compound | μ (10−20) | m (10−20) | D (10−40) length | R (10−40) length | E–M | g abs |
---|---|---|---|---|---|---|
(S,S,P)-2 | 204 | 3.74 | 41568 | 645 | 32 | 0.062 |
(S,S,P)-3 | 382 | 3.97 | 146279 | 1407 | 22 | 0.03 |
(S,S,P)-4 | 426 | 4.19 | 181829 | 1606 | 26 | 0.03 |
(S,S,P)-5 | 442 | 4.20 | 195208 | 1727 | 22 | 0.03 |
(S,S,P)-6 | 524 | 4.42 | 275017 | 2237 | 15 | 0.03 |
(S,S,P)-2:Ag | 215 | 3.73 | 46231 | 705 | 29 | 0.06 |
(S,S,P)-3:Ag | 249 | 3.45 | 62283 | 811 | 20 | 0.05 |
(S,S,P)-4:Ag2 | 216 | 3.84 | 46991 | 815 | 12 | 0.06 |
(S,S,P)-5:Ag2 | 185 | 3.94 | 34572 | 734 | 0 | 0.08 |
(S,S,P)-6:Ag3 | 289 | 4.28 | 83865 | 1190 | 17 | 0.05 |
In that case, the classical electromagnetic power dissipated by light absorption of intensity E02 and wavelength λ by the magnetic field alone in a helix of conductivity σ, cross section a, radius r, pitch length d, and pitch angle Φ, is dictated by P = aσ(E0rsinΦ/4λ)2(d2 + 4π2r2)1/2n. Such a linear relationship is nicely paralleled by the plot in Fig. 8 of the magnetic dipole transition moment |m| versus N, where N, the number of alkynes displayed on the x-axis, is related to the number of loops (n) through n = N/3. This is in agreement with a seminal paper about the CISS effect by Zacharias et al.79 in which the spin polarization increases as a function of the number of turns in helical polypeptide molecules adsorbed on a gold surface, a result that runs parallel with the observed dependence of the transition magnetic moment.
Fig. 8 (a) Linear fitting of |m| with the number of alkynes of compounds (S,S,P)-2–6 in CH2Cl2 solutions in the absence of Ag(I) (black squares, y = 3.38 + 0.8x, R2 = 0.955) and silver complexes (S,S,P)-2:Ag, (S,S,P)-3:Ag, (S,S,P)-4:Ag2, (S,S,P)-5:Ag2 and (S,S,P)-6:Ag3 (pink circles, y = 3.13 + 0.08x, R2 = 0.69). (b) Calculated structure including helix axis and magnetic dipole transition moment of (S,S,P)-6 (for the rest of compounds see Fig. S157 and 158†). |
This nice result is based on the alignment of the magnetic transition vector and the helix axis for the S0 → S1 transition. At this point we hypothesized that the projection of the magnetic transition vector on the helix axis could be even a better predictor owing to the fact that spin filtering implies some sort of directionality, implicit in the idea that the scattered electron should pass through the helix. In this way, transitions with magnetic transition vectors misaligned with respect to the direction of the electron transport (helix axis) would not be relevant for the desired relationship between CISS and chiroptical properties. Although the helix axis is not fully conserved in all the structures (see Fig. S157 and 158†) we could observe again an almost perfect alignment between m and the defined helix axis. Consequently, the linear relationship is preserved (Fig. S159†).
A similar analysis can be carried out with the Ag(I) complexes, although in this case the geometrical variations along the series are higher resulting in a worse linearity (see Fig. 8, pink circles). In any case, a clear increase in the magnetic moment of the transition (m) can be observed for the Ag(I)-saturated structures (S,S,P)-2:Ag, (S,S,P)-3:Ag, (S,S,P)-4:Ag2, (S,S,P)-5:Ag2, (S,S,P)-6:Ag3.
As suggested above we speculate that the magnetic transition vector or its projection in the electron transport direction of a selected transition could be the missing link between the CISS effect and chiroptical properties. This can be supported by the interpretation of the CISS effect as an additional contribution to molecular magnetic response that is triggered in the presence of electron transport through a chiral (helical) structure. A simple theoretical explanation for this dependence of the spin polarization with molecular length in helices can be given in terms of a scattering model including both coherent and incoherent contributions to electron transport.80 All these assumptions, although promising, require further theoretical inquiries.
In the case of the Ag(I) cations, two different binding sites are possible, the alkynes and the pyridines. ECD spectroscopy was highly informative to discriminate between the two alternatives. As in the case of related compound (S,S,P)-3, intensification of the ECD signal at the longest wavelength (Δε = +41 M−1 cm−1, gabs = +0.45 × 10−2) is observed upon the addition of a Ag(I) salt to a solution 2.5 × 10−5 M of (S,S,P)-7, compatible with an increase of the helicity of the system (Fig. S121†). Moreover, titration experiments also showed that this increase was related to the coordination of the ligand to one Ag(I) cation (K7:Ag = 8.21 × 105 ± 13% M−1) (Fig. S132†). Therefore, the interaction of Ag(I) cations with the alkynes seems to be the most favorable interaction for the first added equivalent. Nevertheless, in the presence of an excess of Ag(I) cation the ECD signals drastically decrease suggesting the formation of insoluble material in CH2Cl2, probably by subsequent coordination of the cation with the available pyridine groups. All these results allow us to conclude that carbophilic metals cannot easily direct the self-assembly of metallosupramolecular structures for our systems.
On the other hand, the addition of a solution of trans-PdCl2(MeCN)2 (10 mM) to solutions of (S,S,P)-7 in CD2Cl2 or DMSO-d6 (10 mM) resulted in the appearance of new signals in the corresponding 1H NMR spectra. This fact indicates the presence of new species in solution and suggests the formation of a metallacycle as a result of the coordination of the pyridine rings to the Pd(II) centers (see ESI, Fig. S92–102†). The DOSY NMR spectrum in DMSO-d6 of this mixture (Fig. S103†) agrees with the presence of two species as two diffusion coefficients are observed for the 1H signals. Thus, the new set of signals that appeared after Pd(II) addition shows a lower diffusion coefficient, supporting the coordination of (S,S,P)-7 to the metal center and the self-assembly of the larger metallacycle. By increasing the amount of the added Pd(II) complex, this new set of signals became the main one in the mixture, being practically the only one present when 1–1.2 equiv. of trans-PdCl2(MeCN)2 were added (Fig. S93 and S98†). In this case, the DOSY NMR spectrum confirms that all signals correspond to a single species that diffuses as a whole, as expected, as only one diffusion coefficient is observed (Fig. S65 and S70†).
The clear solution obtained in DMSO-d6 was monitored over time by 1H NMR spectroscopy, showing no significant changes in the spectrum over 6 days (see ESI, Fig. S102†). In CD2Cl2 the evolution of the system was more complicated since a precipitate was formed soon after the addition of the metal when a final concentration of 5 mM was used. After filtration, the 1H NMR spectrum of the remaining solution showed two sets of signals, one corresponding to the free ligand and the other one to the newly formed species, postulated as the ((S,S,P)-7)2Pd2 complex. However, the ratio between species was not maintained over time, as a decrease of the proportion of the coordination complex was observed by 1H NMR monitoring. Better results were obtained when the concentration was lowered to 2.5 mM, as no precipitation was observed. In this case, the set of signals attributed to the coordinated species was present almost exclusively after 24 h and good stability was observed over time (Fig. S96†).
Comparison of the 1H NMR of (S,S,P)-7 before and after the addition of Pd(II) both in DMSO-d6 and CD2Cl2 shows a downfield shift (ΔδHa = 0.23 ppm; ΔδHb = 0.32 ppm in CD2Cl2 and ΔδHa = 0.22 ppm; ΔδHb = 0.23 ppm in DMSO-d6) of the H atoms of the pyridine ring at positions 2 and 6 (Fig. 9a). This is the expected behavior for the coordination of the pyridine ring to the Pd(II) center. To lend further support to the fact that the observed shift corresponds to the binding event, 3-(phenylethynyl)pyridine was prepared as a model ligand. The similarities between the chemical shift for the signals of the pyridine nuclei and their variation in the two complexes (ΔδHa = 0.23 ppm; ΔδHb = 0.25 ppm in CD2Cl2 and ΔδHa = 0.23 ppm; ΔδHb = 0.19 ppm in DMSO-d6, see Fig. S93 and 104†) support the same coordination pattern in both cases.
Fig. 9 (a) Partial 1H NMR (500 MHz, CD2Cl2) spectra of ligand (S,S,P)-7 (top) and metallosupramolecular complex ((S,S,P)-7)2Pd2, (bottom). The lettering corresponds to that shown in Fig. 1; (b) observed (top) and calculated (bottom) HRMS (MALDI+) isotopic pattern for the [M − Cl]+ ion of metallacycle ((S,S,P)-7)2Pd2; (c) ECD spectra of complex ((S,S,P)-7)2Pd2 (black line) and its enantiomer ((R,R,M)-7)2Pd2 (red line). |
Dilution experiments down to 0.5 mM both in CD2Cl2 or DMSO-d6 show that the self-assembled complex is the only species at 2.5 mM in CD2Cl2 and in the 5–2.5 mM range in DMSO-d6 and the major one down to 1 mM in both solvents. This fact strongly suggests that the more entropically favored complex is being obtained. At 0.5 mM in CD2Cl2, the complex is still the major species in solution, although the free ligand is clearly observed. In contrast, a considerable amount of the bare (S,S,P)-7 ligand is observed in DMSO-d6 at the same concentration. This is not surprising taking into account the higher coordinating character of DMSO compared to CH2Cl2, which leads to a higher self-assembly concentration.33
The self-assembly of a ((S,S,P)-7)2Pd2 complex was further demonstrated by MALDI mass spectrometry of an evaporated sample of the complex. The high-resolution mass spectrum shows a signal at m/z = 2193.3, whose exact mass and isotopic distribution (Fig. 9b) nicely match those calculated for the [M − Cl]+ ion, resulting from the loss of one chloride ligand by one of the Pd centers.
Chiroptical properties were not particularly informative owing to just a moderate distortion of the electronic states after coordination. Therefore, very similar UV-Vis and almost identical ECD spectra to those of the free ligand were obtained (Fig. 9c and S122†). To confirm that we were really observing the signal associated with the complex, the 1H NMR spectra of the sample were recorded before and after the UV-CD measurements, performed at the NMR concentration (Fig. S123 and S124†).
Finally, the structures of metallosupramolecular complex ((S,S,P)-7)2Pd2 and its ((R,R,M)-7)2Pd2 enantiomer were investigated by means of DFT calculations at the ωB97XD/def2SVP level of theory both in DMSO and CD2Cl2 (Fig. 10 and Tables S25–28†). The optimized structures gave structural arrangements which were in good agreement with the experimental data. They showed the coordination of the pyridine rings to the Pd(II) metal centers and the generation of a chiral cavity defined by the helical cores.
Fig. 10 Optimized DFT (ωB97XD/def2SVP, DMSO) structure of ((S,S,P)-7)2Pd2. Color coding: C, grey; H, white; O, red, Pd, yellow; Cl, green. |
Footnote |
† Electronic supplementary information (ESI) available: Additional figures and tables referred to in the main text, the experimental details for NMR, IR, ECD, CPL and VCD measurements, the details of the conformational analysis, the computational aspects of the spectral calculations and data treatment. CCDC 2073558. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1qo00822f |
This journal is © the Partner Organisations 2021 |