Jonathan
Álvarez-García
a,
Víctor
Rubio-Pisabarro
a,
Luis
García-Río
b and
María
Magdalena Cid
*a
aDepartamento de Química Orgánica, Edificio Ciencias Experimentais, Campus Lagoas-Marcosende, Vigo E-36310, Spain. E-mail: mcid@uvigo.gal
bDepartamento de Química Física, Facultade de Química, Avda das Ciencias, s/n, Santiago de Compostela, E-15782, Spain
First published on 18th September 2023
Proton transfer equilibria are of pivotal importance due to the role they play in a plethora of biological and phamaceutical processes. With the aim of explaining the relative position of the hydrogen to be transferred, we investigated the behavior of chiral pyrido-cyclophanes with different morphologies using circular dichroism in the presence of different acids in acetonitrile. The results showed that all three compounds underwent double protonation and formation of cascade ion-pairs, leading to the appearance of diagnostic signals in their ECD spectra. The presence of water fosters the crystallization of several intermediates that do not correspond to those in solution. By using Brønsted correlations, it was found that proton transfer from the acid to the pyridine occurred regardless of the acid's pKa.
Total or partial hydrogen transfer between acid and base partners results in the formation of salts or a hydrogen-bonding driven complex, respectively. The degree of proton transfer between two species in the solid state can be predicted using the pKa rule, that suggests the potential formation of a co-crystal when the difference in pKa values of the interacting partners, denoted as ΔpKa, is below 0, whereas if it exceeds 3, the formation of a molecular salt is likely.6,7 Nonetheless, there are always exceptions and for an intermediate value no accurate prediction can be made. On the other hand, in solution, salts are present as separated, contact or aggregated ion-pairs depending on the medium polarity. In particular, ion-pairing, in which understanding the degree of proton transfer is essential, explain unexpected observations in the behaviour of synthetic ion receptors.8–10
Shape-persistent pyridocyclophanes have been the subject of research in Supramolecular Chemistry and Material Science due, in part, to the properties of the pyridine motif.11–14 So, it is not surprising that proton transfer between acids and pyridine nitrogen atoms had drawn considerable attention.15–18
Our research group has exploited pyridines and allenes as construction motifs in the synthesis of chiral shape-persistent cyclophanes with applicability in the field of molecular recognition in non-aqueous solvents.19 We reported the preparation and full characterization of pyrido-cyclophane 1 (hereafter, pyrido-allenophane to evidence the presence of allenyl motifs, Fig. 1), which showed one of the strongest chiroptical responses described for small organic molecules with a g-factor‡ of 0.011.20 However, its rapid racemization under ambient conditions, driven by ring strain, did not allow us to exploit its sensing capabilities. With the aim of overcoming ring strain, we synthesized bipyrido-allenophane 2 and found that, in the presence of TfOH, it not only undergoes a stereoselective double protonation process, which is quite remarkable for a cis-locked bipyridine, but also captures one of the corresponding counteranions forming a contact ion-pair (also referred to as a cascade complex, in which the anion coordinates to the cation center).20 In this respect, pyridine-decorated chiral cyclic peptides have shown simultaneous complexation of cations and anions.21
![]() | ||
Fig. 1 Potential pyrido-allenophanes as ion-pair receptors: known 1, 2 and 4, and new [72]-pyridoallenophane 3. |
Our aim is to investigate the degree of the proton transfer in acid–base processes involving pyridine derivatives via chiroptical signals in non-aqueous media. For such endeavor, allenes are used to provide axial chirality to pyrido-cyclophanic systems (Fig. 1).
Having laid the foundation with bipyrido-allenophane 2, and taking into account the interest that proton transfer phenomena and pyridine H-bonding have, we decided to expand our study to a greater number of acids. Besides, we now envision the introduction of a pyridine motif into the butadiyne bond of 1 to afford the new pyridoallenophane 3, since, in addition to its role as a stability factor, the second pyridine moiety provides an extra recognition site. In addition, we want to sort out whether allenophane 4,22,23 whose conformational space in solution is defined mainly by two conformers, would undergo a conformation-selective molecular recognition process and thus susceptible to being monitored by circular dichroism (CD).
So, here we present the preparation and complete characterization of the new chiral macrocycle [72]-pyrido-allenophane 3, the behavior of allenophanes 2, 3 and 4 against different acids (oxoacids and HCl) in acetonitrile and their abilities to capture the counteranions and consequently reduce their conformational space via ion-pair recognition. These three chiral allenophanes each contain two pyridine units in three different spatial morphologies, as the aromatic spacers are separated by zero (2), one (3) or two (4) allene units, which could supply crucial information to correlate chirality signals with geometry rigidity and recognition properties.
![]() | ||
Scheme 1 Synthetic sequence followed to obtain macrocycle 3. Bottom: Crystal structure resolved by X-ray diffraction ((P2)-3) and conformation obtained by DFT computations (CAM-B3LYP/6-31g+(d,p)). |
Allenophane 3 exhibited a maximum dissymmetry factor‡g = 0.006, which is the same as that of the [142] analogue 4, even though the number of chiral units is reduced by 50%. This result supports the importance of the ratio between the number of chiral units and conformation stability to render outstanding chiral systems. Also, in contrast to compound 1, 3 is configurationally stable under ambient conditions of light and temperature.
Previous research by our group showed that both CD and NMR titration experiments gave consistent equilibrium constants (K) values for allenophanes in the presence of TFA in acetonitrile.20,31 In the current study, we therefore focused primarily on CD spectroscopy due to its sensitivity to conformational changes32,33 and the large chiroptical responses of this family of macrocycles. In addition, CD has proven to be a very useful tool for studying protonation of other chiral pyridocyclophanes.34 The details for each titration experiment may be found in the ESI.¶
For the five oxoacids tested, the CD spectrum (Fig. 2A, red line) evolved until a first saturation point was reached that is characterized by two negative bands centered at 290 and 370 nm (Fig. 2A, black line). Upon further increasing the acid concentration, the band at 370 nm moved to 388 nm and a new band appeared at 330 nm that is diagnostic for fixation of the bipyridine axis in the M configuration (Fig. 2A, blue line), as we previously showed for TFA.20
Given the results from allenophane 2, pyridoallenophane 3 was treated with TfOH, H2SO4 and TFA as representative acids. Upon increasing amounts of acid in acetonitrile, allenophane 3 showed important changes in its CD spectrum. Initial negative bands at 302 and 320 nm (Fig. 2B, red line) gradually transformed into a new positive band with a maximum centered at 302 nm and two new negative bands appeared with minima at 340 and 375 nm (Fig. 2B, blue line).
We found that the experimental data obtained from the different titrations for 2 and 3 could be fitted mathematically to the same set of equilibria (Scheme 2), although with different association constants. From a general point of view, at low acid concentrations, the free macrocycle begins to disappear in favor of the monoprotonated (AP·H)+ and AP·HA species (AP = allenophane 2 or 3). Later, when the concentration of acid in solution increases, both species disappear in favor of the doubly protonated macrocycle complexed with a counteranion (AP·H2A)+, which turns out to be the final species in all the titrations.
Finally, we evaluated the behavior of allenophane 4 in the presence of TfOH, H2SO4 and TFA. We observed the disappearance of the initial negative band centered at 330 nm (Fig. 2C, red line) followed by the gradual rise of three new diagnostic bands, one negative at 384 nm and two positives at 351 and 305 nm (Fig. 2C, blue line). In this case, the titration data suggest that 4·H2A2 is the final species (Scheme 2, orange), in which the free macrocycle captures two molecules of acid in contrast to macrocycles 2 and 3, whose endpoints correspond to the diprotonated species capturing one counteranion.
Inspired by the ability of macrocycle 4 to form 1:
2 complexes with several acids, we directed our efforts to achieve a 1
:
1 host–guest inclusion complex and obtain clues about the size and adaptability of the binding cavity. For this purpose, we thought that 1,2-ethanedisulfonic acid might be a suitable target for two reasons. First, it could function as a bidentate guest, and second, sulfonate counterions are key players in the pharmaceutical salt portfolio because of their unique properties, which include high melting points, fair aqueous solubility and stability, and the absence of hydrates. As a consequence, sulfonates often offer attractive advantages over the hydrochloride salt forms.35 Thus, the use of this bidentate acid as titrant led to the same experimental CD outcome as the other oxoacids but this time forming a 1
:
1 complex.
Our attempts to crystallize several of these complexes under anhydrous conditions were unsuccessful. However, using commercial TfOH, where traces of water are to be expected, we managed to crystallize allenophane 3 both in its monoprotonated (3·H)+ form (Fig. 3A) from acetonitrile and the diprotonated (3·H2)+2 (Fig. 3B) form from a 1:
1 acetonitrile
:
ethanol solution at 4 °C.§ Both have an encapsulated water molecule slightly displaced towards one of the pyridine rings. In (3·H)+, a water molecule acts as a bridge between both Npyr and the triflate anion. In (3·H2)+2, the water molecule acts as an H-bond acceptor for both NpyrHs and as a donor for the two counterions that are located on opposite sides of the macrocycle. In neither of the structures was the initial conformation of the starting macrocycle modified.
![]() | ||
Fig. 3 (A) Crystal structure of 3·H2O·OTf. (B) Crystal structure of 3·(H2O·OTf)2. (C) Crystal structure of 4·(H2O·OTf)2: asymmetric units I and II. Hydrogens were omitted for clarity. |
Under the same conditions, allenophane 4 crystallized in the form of (4·H2)+2 from a 8:
2 dichloromethane
:
ethanol solution.§ Here, the asymmetric unit contains two slightly different conformations of the allenophane framework: I- each of the two water molecules establishes H-bonding with NpyrHs and with a triflate anion (Fig. 3CI); II- both pyridine rings act independently due to the orientation of respective N–H bonds (Fig. 3CII). In the three crystals, the protons are located on the pyridine rings, in accord with predictions made by the pKa rule for triflic acid6 and the equilibrium constants previously obtained in our mathematical analysis (see Table S3†).
Next, we wanted to address the question “Where is the proton located?” in the final species of each titration. Different techniques have been used for this purpose, including NMR,36 IR,18,37 XRPD7,38 and polarography,39 although mainly in carboxylic acids. Here we apply a simple and intuitive Brønsted correlation, plotting the logK2 obtained from the mathematical fits (Table S3†) against the pKa of the tested oxoacids. As shown in Fig. 4, we observed linear correlations for all three macrocycles, indicating that regardless of the pKa the extent of proton transfer must be the same for all the acids tested.
![]() | ||
Fig. 4 Brønsted plots of log![]() |
For allenophanes 3 and 4, the lines are practically parallel to each other with a slope α ≈ |1|, which indicates complete transfer of the proton from the acid to the pyridine (Fig. 4a). In the case of 2, the slope turns out to be α ≈ |0.43|, which a priori could indicate that the proton is equally shared between the acid and the pyridine.|| However, the pyridines in 2 are not independent systems, so we are inclined to think that the proton has been transferred as in 3 and 4 and that it is shared by the bipyridine moiety (Fig. 4b).
Surprisingly, when bipyrido-allenophane 2 was titrated with HCl, only the first saturation point was reached (Fig. 2A, black line). The spectrum did not evolve further and diagnostic bands for the stereoselective discrimination of the bipyridine axis never appeared, pointing to a different final species. In contrast, treatment of a solution of pyrido-allenophane 3 in acetonitrile with HCl/dioxane led to the same spectrum changes as those seen for the oxoacids (Fig. 2B, blue line). Completely different behavior was observed in the case of pyrido-allenophane 4. At low acid concentrations, the experiment yielded the same outcome as the oxoacids (Fig. 2C, blue line), but as the concentration of acid increased, the bands at 384 and 351 nm intensified and the positive band at 305 nm became a weak negative band (Fig. 2C, black line).
These observations reveal a pivotal role of the counteranion in selecting the conformation of the final ion-pair complex. When fitting the experimental data to calculate the equilibrium constants, we have to keep in mind that the protonating species can be either H+, HCl or (H·dioxane)+. To make it simpler and since our main interest is to determine which species is the final one, we performed the mathematical analysis by just considering the pKa of HCl in acetonitrile (see the ESI† for details). The mathematical fitting suggested that the final species are 2·H2Cl2/(2·H2·Dioxane2)+2, (3·H2Cl)+/(3·H2·Dioxane)+2 and 4·H2Cl2/(4·H2·Dioxane2)+2, with calculated CD spectra that closely resemble the experimental ones (see below).
The conformational space of allenophane 2 was described by two possible conformers depending on the configuration adopted by the axis that joins both pyridine rings, namely P or M. When we computed and optimized the structure of the complexes with all five oxoacids (2·H2A)+, we found that the bipyridine nitrogens share the two protons with an oxygen atom of the corresponding oxoanion while the side chain of each of the acids is located towards the outside of the cavity, moving away from the tert-butyl groups (Fig. 5A). As a result of this interaction, the bipyridine axis is fixed in an M configuration. The calculated spectrum of this species (Fig. 5A, dash green line) coincides remarkably well with the experimental one, and when the axis adopts only the M configuration does the diagnostic band at 330 nm appear.
In the case of HCl, the computationally optimized structure of both 2·H2A2 complexes (A being Cl− or dioxane) shows the guest molecules located on each of the faces of the macrocycle (Fig. 5, bottom). Since bipyridine nitrogens do not coordinate the same molecule cooperatively, the bipyridine scaffold can continue to adopt both configurations, namely M and P, and therefore the diagnostic band at 330 nm never appeared. Interestingly, the optimized 2·H2A2 structure is very reminiscent of the structure obtained by X-ray diffraction of macrocycle 2 in chloroform20 (Fig. 5, bottom).
For allenophane 3, the calculated geometries of the complexes (3·H2A)+ showed the oxoanions acting as bidentate ligands and consequently forcing the macrocycle to adopt a V-shaped conformation (Fig. 5B). The calculated spectrum of these species (Fig. 5B, dash green line) coincides remarkably well with the final experimental spectrum, and only when the macrocycle acquires the new conformation, in addition to being diprotonated, does the positive diagnostic band at 302 nm appear, supporting the idea that one of the anions must be located inside the macrocycle. To ensure this V-like conformation in the case of HCl, we ascribed a bidentate character to the dioxane molecules. Indeed, the calculated and experimental CD spectra were coincident when a dioxane molecule was introduced as a bridge between the two N–Hs via the oxygen atoms, with an average NH⋯O distance of 1.86 Å.**
On the other hand, the conformational space in solution of allenophane 4 consists of two conformers: chair and twist.23 So, once again, the geometries of the 4·H2A2 complexes for the oxoacids were optimized and their CD spectra calculated. The outcomes showed a new half-chair conformation in which each anion is oriented towards the cavity of the allenophanic ring with an average distance of 10.50 Å between the two nitrogens of the pyridine moieties.** The calculated spectra of these complexes reproduced the experimental spectrum very well (Fig. 5C, dashed green line). In the case of 1,2-ethanedisulfonic acid, simultaneous interactions between the two acidic moieties and the two pyridine rings were observed (Fig. 5C), with the two sulfonic groups adopting an antiperiplanar disposition.
Modeling the obtained complex upon treatment of 4 with HCl proved to be more of a challenge. In order to arrive at a computed geometry whose CD spectrum resembled that seen experimentally, we devised a system in which the two protonated pyridines and two water molecules were bridged by one dioxane molecule (Fig. 5D, dash green line). The optimized geometry here showed the macrocycle adopting a new distorted twist conformation, with bond distances of 1.65 Å and 1.77 Å for the NH⋯O and HO⋯O(CH2)4O interactions, respectively. In addition, the distance between the two nitrogens of the pyridine moieties is larger (12.98 Å) than that in the complexes formed with the other acids (10.50 Å) (Fig. 5C).** The difference between the two new conformations in the final species explains the different outcomes in their ECD spectra.
It should be mentioned that the ECD spectra of the salt geometries obtained by X-ray diffraction were also computed. In none of the cases, the ECDs calculated coincided with the experimental end point of the titrations nor with those of the optimized crystallographic geometries (Fig. S73, S74 and S80†). Therefore, the complexes in solution at the end point must be different from those obtained in the solid state in the presence of water. This supports the idea that the geometries of the complexes shown in Fig. 5 should be responsible for the observed changes in the ECD during the titrations.
Furthermore, in all the optimized structures the protons were located on the pyridine moieties, matching their computed spectra with the experimental results. This agrees with what the Brønsted correlation called for.
Footnotes |
† Electronic supplementary information (ESI) available: Spectroscopic data and experimental details for the preparation of 3. Crystal data of 6, 3, 3·H2O·OTf, 3·(H2O·OTf)2 and 4·(H2O·OTf)2; geometries and spectra for the calculated species. CCDC 2280776, 2280778, 2280777, 2280775 and 2280779. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3qo01180a |
‡ Dissymmetry factor (g) quantifies the effective chiroptical activity of a system by correlating the circular dichroism signal intensity with the average absorption signal intensity at a given frequency. |
§ CCDC 2280776 contains the supplementary crystallographic data for (P)-6, 2280778 for (P2)-3, 2280777 for (rac)-3·TfOH, 2280775 for (rac)-3·(TfOH)2 and 2280779 for (P4)-4·(TfOH)2.† |
¶ The allenophanes 2, 3 and 4 are configurational stable under the tested conditions and upon addition of Et3N the original circular dichroism spectra were recovered. |
|| The hydrogen coordination to both nitrogens yields a larger N⋯H bond that is compatible with a partial hydrogen transfer. A plausible interpretation of this result implies that the counterion acts as a shuttle for the mobile H+, provided both pyridines behave independently. |
** The measurements were made using GaussView 6 software. |
This journal is © the Partner Organisations 2023 |