Open Access Article
Benjámin
Kovács
a,
Tamás
Földes
b,
Márk
Szabó
a,
Éva
Dorkó
b,
Bianka
Kótai
b,
Gergely
Laczkó
bc,
Tamás
Holczbauer
b,
Attila
Domján
*a,
Imre
Pápai
*b and
Tibor
Soós
*b
aCentre for Structural Science, HUN-REN Research Centre for Natural Sciences, Magyar tudósok körútja 2, H-1117 Budapest, Hungary. E-mail: domjan.attila@ttk.hu
bInstitute of Organic Chemistry, HUN-REN Research Centre for Natural Sciences, Magyar tudósok körútja 2, H-1117 Budapest, Hungary. E-mail: papai.imre@ttk.hu; soos.tibor@ttk.hu
cHevesy György PhD School of Chemistry, Eötvös Loránd University, P. O. Box 32, Budapest, H-1518, Hungary
First published on 26th August 2024
Using the principle that constrained conformational spaces can generate novel and hidden molecular properties, we challenged the commonly held perception that a single-centered Lewis acid reacting with a single-centered Lewis base always forms a single Lewis adduct. Accordingly, the emergence of single-centered but multiple Lewis acidity among sterically hindered and non-symmetric triaryl-boranes is reported. These Lewis acids feature several diastereotopic faces providing multiple binding sites at the same Lewis acid center in the interaction with Lewis bases giving rise to adducts with diastereomeric structures. We demonstrate that with a proper choice of the base, atropisomeric adduct species can be formed that interconvert via the dissociative mechanism rather than conformational isomerism. The existence of this exotic and peculiar molecular phenomenon was experimentally confirmed by the formation of atropisomeric piperidine-borane adducts using state-of-the-art NMR techniques in combination with computational methods.
As molecular geometry and function are interrelated, there has been growing interest in remodeling the conformational space and designing molecules that can adopt the desired shape with minimal or no energetic penalty. The common strategy that guided these efforts was to impose structural constraints via steric repulsion or ring strains. As a result of these challenging undertakings, peculiar reactivities,3 well-designed emissive systems,4 rarely observed metastable conformational and redox states5 and enhanced ligand bindings6 have been engineered.
Along the above lines, we took up the challenge of demonstrating the existence of a peculiar molecular phenomenon that enables multiple Lewis acidity (MLA) to be assigned to a single-centered Lewis acid (LA). IUPAC defines Lewis acidity as ‘the thermodynamic tendency of a substrate to act as a Lewis acid’7 which is routinely described by a single equilibrium constant associated with a single dative adduct. In contrast, we envisioned that sterically crowded and non-symmetric triaryl-boranes, i.e. those involving an LA center with diastereotopic faces and confinement,8 are able to form distinct dative LA–LB adducts9 with the same single-centered achiral Lewis base (LB). Accordingly, these Lewis adducts represent different stereoisomeric states formed in competitive association pathways. Due to hindered rotations around single bonds, these LA–LB stereoisomers are diastereomeric atropisomers, which can be distinguished at ambient temperature, and thus, the Lewis acidity of these boranes can be described by not a single, but multiple equilibrium constants. Atropisomerism is a unique case of stereoisomerism caused by restricted rotation around a single covalent bond that is not constrained in a ring. While this phenomenon commonly refers to C(sp2)–C(sp2) single bonds, there are some less common examples for atropisomerism as a result of restricted rotation around C(sp3)–C(sp2) or even C(sp3)–C(sp3) single bonds (see Fig. 1 for selected examples).10 In this study we demonstrate that atropisomeric LA–LB adducts with dual restricted rotation around B(sp3)–N(sp3) and B(sp3)–Ar(sp2) bonds are formed between fluoro-chloro triaryl-boranes and amines within the same reaction.
![]() | ||
| Fig. 1 Selected examples of molecules featuring (A) C(sp2)–C(sp2) and (B) C(sp3)–C(sp3) atropisomeric states10 and a (C) dual B(sp3)–N(sp3) and B(sp3)–Ar(sp2) type atropisomeric adducts reported in our present work. Bonds corresponding to restricted rotations are highlighted in red. | ||
The challenge outlined inevitably poses an additional fundamental question: when do we consider configurational, conformational, and stereochemical isomers as distinct, isomeric molecules with different properties, including reactivity and symmetry? This question cannot be answered without specifying important experimental conditions such as the time scale of observation, temperature, pressure, etc. That is, these molecular properties are definable only within limits. As an example, owing to the sterical congestion, R- and S-binaphthol (Fig. 1A) are atropisomers at ambient temperature and fundamental chiral building blocks for chiral catalysis. So, these conformational isomers are distinct molecules at room temperature with different reactivities in asymmetric reactions. Nevertheless, over 200 °C reaction temperature or million years of observation period, R- and S-binaphthols should be treated as non-distinct, equivalent molecules.
In virtue of above, we had to specify the external parameters under which the existence of MLA was investigated. Considering practical aspects and the potential of future use, we have confined ourselves to the most relevant parameters of routine organic chemistry (e.g. ambient temperature, 1 atm pressure, and ms-days timescale). Thus, to be distinguished as distinct molecules, the formed isomeric dative LA–LB adducts must be detectable by NMR at ambient temperature. Additionally, the structure of the LA should be such that its conformational isomers can be handled as non-distinct, equivalent molecules at ambient temperature. So, the sterical congestion around the Lewis acidic center must be large enough to form atropisomeric dative adducts, but not too large to halt the conformational isomerism of the LA itself at ambient temperature.
It is then worthwhile to further specify the phenomenon of single-centered, but multiple Lewis acidity and to discuss the challenges associated with its observation. First and foremost, the presence of non-equivalent, diastereotopic faces as distinct ‘binding sites’ of the Lewis acid is a necessary but not sufficient condition for MLA. It is possible that, for steric and/or electronic reasons, only one dative adduct is produced, i.e. a diastereoselective adduct formation occurs and single Lewis acidity can be determined (Fig. 2A).9a,e If the formation of multiple adducts is kinetically and thermodynamically possible, we might expect that all observable atropisomeric dative adducts are formed in different, mutually competitive chemical reactions between the same LA and LB. However, if the thermodynamic equilibrium between the detected LA–LB adducts is also established by rapid conformational isomerisation (via rotation around single bonds) as well as slow dissociation–association, then the situation should be considered a case of single Lewis acidity, which, again, can be described by a single equilibrium constant (Fig. 2B). Thus, multiple Lewis acidity can only emerge, if either stable, non-dissociable atropisomeric dative adducts are formed between the diastereotopic LA and the same LB (in a kinetically controlled manner) or the interconversion of the atropisomeric dative adducts requires the breaking and reforming of the dative bond (Fig. 2C). Indeed, in the latter case each LA–LB atropisomer represents a reversible association pathway, and thus, the Lewis acidity can be described by as many equilibrium constants as the number of LA–LB atropisomer detected. Finally, it should be noted that not only the strength of Lewis acidity (hydricity, Gutmann–Beckett, etc.), but also its multiplicity is a situation-dependent phenomenon, so it can happen that the very same LA can exhibit single or multiple Lewis acidity with different LBs.
In this paper, we report chemical evidence that the hypothesized phenomenon that enables multiple Lewis acidity to be assigned exists, even at ambient temperature. Indeed, combined computational and NMR studies revealed that the crowded LA–LB atropisomeric adducts (Fig. 1C) display such high rotational barriers that their isomerization pathway is dissociative.
O or the hydride anion (H−) through the counteracting change in steric effects (front- and back-strains).
![]() | ||
| Fig. 3 (A) Illustration of the diastereotopicity of non-symmetrically halogenated triaryl-borane catalysts I, II and III (highlighted with blue, orange and cream backgrounds). Due to the presence of only one non-symmetrically substituted aryl ring, borane III has only one rotamer. The rotameric states of each borane are shown in van der Waals representations from top and bottom views. The corresponding rotameric states are denoted with Greek letters which refer to the order of their calculated stability, with the α-state being the respective lowest energy state (Table S2†). (B) Illustration of the torsional anisotropy upon dative bond formation with appropriate Lewis bases (having no C3 symmetry) using the top and bottom views of the I(α) structure. The colored arrows indicate distinct molecular confinements for the incoming Lewis base. The F atoms are colored blue, and the Cl atoms are colored green. | ||
This interesting and somewhat unusual behavior upon dative adduct formation prompted us to conduct further and more detailed structural studies with boranes I–III. After a thorough analysis of their structures, we postulated that each of these special boranes individually features diastereotopic faces (Fig. 3), possibly resulting in diastereomeric adducts upon interaction with an achiral LB and thus might possess MLA. In these BAr3 structures, the aryl rings adopt propeller shaped conformations (see the computed I–III rotameric states in ESI Section 2 in the ESI†); however, the particular o-F and o-Cl substitution patterns decrease the molecular symmetry and these boranes can display both directional and torsional structural anisotropy upon dative adduct formation. In contrast to III, a more complex situation is expected to emerge for boranes I and II, as the rotation of the non-symmetrically substituted aryl rings around the B(sp2)–Ar(sp2) single bonds results in multiple rotamers (2 and 3 for I and II, respectively), which, in principle, would further diversify the global and intrinsic Lewis acidities and their diastereotopicity. Additionally, we also hypothesized that the distinct conformations of the dative borane adducts might be detected as the enhanced rotational barriers imposed upon the dative adduct formation would hinder the rotational isomerism of the thermally accessible conformational states. Accordingly, these isomeric LA–LB adducts would be atropisomers, stemming from directional and/or torsional anisotropy of the LA upon dative bond formation and the hindered internal rotations around the B(sp3)–Ar(sp2) and B(sp3)–N(sp3) single bonds (using N-based LB).
![]() | ||
| Fig. 4 The schematic representation of the triaryl-borane structures and their hydride adducts (F: blue, Cl: green, and B: pink spheres). The ‘top’ and ‘bottom’ approaches refer to the diastereotopic views presented in Fig. 3. In the sp3 states the spatial density of halogen ligands differs on the two sides of the adducts leading to increased back-strain (congested ligands) and decreased front-strain (sparser ligands). The combination of rotameric states and diastereotopic faces engenders multiple, atropisomeric hydride adduct conformers for each borane. The respective ΔGh values are shown in Table S3.† Here, the ΔΔGh values with respect to the most stable conformational state (α) are indicated for the higher energy hydride adduct conformers in kcal mol−1 at 298.15 K. | ||
![]() | ||
Fig. 5 The 1D 19F NMR spectrum (564 MHz, 298.2 K) of a mixture of p and II in toluene-d8 using a p to II ratio of 3 : 1 without (top) and with (bottom) 1H-decoupling during acquisition. The two F atoms in p–II combined with the three distinct atropisomeric states (α, β, or γ) generate six resolved p–II19F signals. The represented conformation and specific F atom position (aV, aT, hV, hT, sV or sT, see Fig. 6A) are indicated for each p–II signal. For the p–II(α) signals the effect of 19F–19F through a space J-coupling (JFF) of 39.5 Hz between the aT and hT positioned 19F spins is apparent as in the 1H-decoupled spectrum they remain doublets as opposed to the p–II(β) and p–II(γ) signals which simplify to a singlet. The broad signal marked as ‘w–II’ that remains close to the sp2-hybridized borane peak position represents slowly diffusing species that comprise interacting II, p, and residual water molecules encompassed in a single solvent cage (the characterization of these aqueous species is detailed in ESI Section 9†). | ||
To aid the structural assignment of the expected atropisomeric states of p–I, p–II and p–III, computational studies were conducted (ESI Section 4†) following the structure calculation protocol used for the hydride adduct calculations (all computed p-borane structures along with the corresponding ΔGa values are reported in Tables S4–S6†). Importantly, the calculated dative p-adduct structures display clear differences from the hydride adducts as a result of enhanced steric effects and the lower symmetry of the p LB. Nevertheless, they share a common molecular topology regarding the relative orientation of p and the borane's aryl rings. This is illustrated in Fig. 6A with the energetically lowest lying (α) conformer of p–II. The structurally common features are (i) the N atom of p donates the equatorial lone electron pair to the boron to establish the dative bond, (ii) one borane’ aryl ring is oriented vertically thereby pointing one of its halogen atoms towards the amine proton of p (HN) to form a hydrogen bond, and (iii) the second borane’ aryl ring, which also interacts with HN, is tilted about 45° from the B–N axis (iv) while the third aromatic ring is in horizontal orientation, that is, nearly perpendicular to the B–N axis. We designated the emerging halogen atom positions on the vertical and tilted aryl rings as being in the ‘syn’ or ‘anti’ position with respect to the incoming p. Thus, given the observed ring topology preference, any p-adduct atropisomer can be defined in a straightforward way through the specific F/Cl spatial distribution pattern between the horizontal (h), syn (s) and anti (a) position pairs. Notably, due to the emerging hydrogen bond with the syn positioned aryl halides and back-strain interactions, in all three p-borane systems, the vertical and tilted fluorinated ring(s) prefer the F atom in the syn and anti positions, respectively. Also, the horizontal fluorinated ring adopts a more stable orientation with the F atom pointing towards the tilted ring. Based on the computed Gibbs free energies of adduct formation (ΔGa) the expected affinity of the boranes for dative bonding with p follows the general order of III < II < I which indicates the impact of the growing back-strain to destabilize the dative adducts (Tables S4–S6†). The importance of back-strain effects in adduct formation is supported by activation-strain analysis presented in ESI Section 11.† We also note that the presence of hydrogen bonds in these piperidine-borane adducts may contribute considerably to the thermodynamic stability of the various atropisomeric adducts (although the most stable isomers feature the weaker Cl⋯H–N bond), but it is unlikely that the hydrogen bonding interactions (acting as a ‘structural pin’17) are the prerequisite of the formation of distinguishable atropisomeric states. This is confirmed by the fact that borane-LB adducts formed with bases having no N–H protons (for instance LB = quinuclidine, see Section 12 of the ESI†) also show multiple resonances in the 19F NMR spectra, pointing to multiple atropisomeric adduct states.
![]() | ||
| Fig. 6 (A) The common ring topology of the computed p-adduct structures illustrated with the lowest lying (α) conformation of p–II. The emerging halogen atom positions are denoted as horizontal (h), syn (s) or anti (a) whereas the subscripts specify whether the halogen atom is located on (s, a) or points towards (h) the tilted (T) or vertical (V) ring. The hydrogen bonding interaction between the HN and the syn-positioned ligands is indicated with black dashed lines. The p–II(α) features the two F atoms (blue) in the aT and hT positions through which a through-space JFF coupling of 39.5 Hz emerges (see α(aT) and α(hT) signals in Fig. 5). (B) The 19F–1H HOESY correlations of the p–II atropisomers (α, β and γ) are highlighted with circles. The s-fluorines (of the β and γ atropisomers) are found sufficiently close in space to most of the p ring to be correlated with the HN and most of the alkyl p hydrogens; the h-fluorines (α, γ) are correlated with a limited number of alkyl p hydrogens; in contrast, the a-fluorines (α and β) do not display HOESY contacts with p. | ||
To validate the assignment of the experimentally observed p-borane atropisomeric states to the computed structures, 19F chemical shift calculations (ESI Section 5†) at the B3LYP-D3/6-311+G(2d,p) level of DFT and 19F–1H HOESY18 measurements (Fig. 6B, S28 and S29†) were performed. These revealed correlations in the frequency, linewidth, and fine structure of the 19F signals with the fluorine atom allocated to the proposed horizontal, syn or anti position (for the full 19F NMR signal assignment process see ESI Section 6†). That is, the hydrogen bonding interaction of syn positioned fluorines with the p N–H group is reflected by the upfield shift of their 19F signals and the unusually high 1H–19F J-coupling of 20–25 Hz that exists through the H⋯F hydrogen bond.19 By distinguishing the downfield-shifted horizontal and anti-positioned 19F signals as well, we concluded that the respective 4, 3 and 3 computed atropisomeric states of p–I, p–II and p–III indeed form under experimental conditions (Table 1). Our assignment was further supported by the good agreement between the experimental and computed 19F chemical shifts (Fig. S30†) and conformational probability distributions (Table S10†). The spatial vicinity of fluorine atoms in the p–I(α), p–I(β), p–I(δ) and p–II(α) structures could be confirmed by 19F–19F through space J-couplings20 (JFF) ranging between 15 Hz and 80 Hz (all the extracted JFF values are specified in Table S9†).
| p-Borane adduct | Detected atropisomer | |||
|---|---|---|---|---|
| α | β | γ | δ | |
| p–I | a T s V h T (85%) | a T a V h T (10%) | a T s V h V (3%) | s V s T h T (2%) |
| p–II | a T h T (51%) | a T s V (30%) | h T s V (19%) | |
| p–III | a T (61%) | h T (30%) | s V (9%) | |
:
borane = 1
:
1 ratio or excess of piperidine, the sp2 borane signal remains unobservable due to the high conversion rate of piperidine and the boranes into the corresponding dative adducts. Thus, to detect the sp2 signal, excess of the borane is required.
Further insights into the kinetic stability of the observed p-adducts were also gained through the temperature dependence of the linewidths in the corresponding 19F-NMR spectra. The NMR peaks of p–III resolve at low temperatures (<253 K) only; however, they display severe line broadenings throughout (Fig. S15†). At ∼273 K, the signals start to coalesce, and at higher temperatures a weighted average signal is visible that represents the fast-exchanging ensemble of the p–III atropisomers and the borane III (which notedly incorporates the ‘frustrated’ state). The signals of p–II display moderate line broadenings starting at 293 K and they remain in the slow exchange regime even at 333 K (Fig. S16†); while the signals of p–I do not display any broadenings unless the temperature is raised to 323 K (Fig. S17†). These results imply that the kinetic stability of the p-adducts follows the order of III < II < I which is in agreement with the thermodynamic trend predicted by computation.
:
borane ratio was inverted to 1
:
2 in the mixtures of p with I or II. Under optimal experimental conditions (temperature, mixing time) the corresponding 2D 19F EXSY spectra show qualitatively the same behavior for all the p-adducts, see that of the p + II system as illustrated in Fig. 7 (p + III: ESI Fig. S31, p + I: Fig. S33†). That is, each fluorine resonance of any chosen p-adduct atropisomer shows correlation peaks with those of its other atropisomers and with the corresponding sp2 form as well. These respectively indicate that on the timescale of the experiments the simultaneously existing p-adduct atropisomeric states are in slow exchange with each other; meanwhile, they undergo rapid dissociation and reassociation of the B–N dative bond through the dissociated states of the atropisomeric adduct. The exchange processes and the related barriers are illustrated in Fig. 8 in terms of schematic free energy surface projections. Based on the corresponding conformational exchange rates and dissociation rates (p–I: Table S16, p–II: Table S13, and p–III: Table S11†) extracted from the EXSY peak integrals (ESI Section 1.1†), for any given atropisomeric state of p–I or p–III its exchange with another atropisomer is generally 1–2 orders of magnitude slower than its dissociation. Unfortunately, signal overlaps prevented us from unequivocally quantifying the exchange rates for p–II with the free borane II (Fig. S10†). Nevertheless, a 2D 19F EXSY spectrum was recorded on the p–II system with minimal mixing time (2 ms) that displays the exchanges of p–II with the free borane II while the atropisomeric exchanges remain unobservable (Fig. S32†). Furthermore, for p–I and p–II the atropisomeric exchanges appear to be parallel in terms of fluorine NMR exchanges, that is, any F atom of a detected atropisomer undergoes chemical exchanges with all the F atoms of each detected atropisomer. The exchange rates of the different F-pairs that can be associated with the same atropisomeric exchange fall within an order of magnitude (Table S12†). Consequentially, upon an atropisomeric exchange process there is no correlation between the initial and final position of the F atoms, that is, the orientation of the aryl rings. All these imply that, for the atropisomeric exchanges of these p-borane adducts, neither free nor correlated rotation22 around B(sp3)–Ar(sp2) or B(sp3)–N(sp3) single bonds of the dative adduct is operative. Thus, the atropisomeric exchanges for the investigated p-adducts, i.e. their conformational isomerism rather take place through the sp2 borane state where the rotations of the aryl rings are free and essentially, they occur as the side-processes of the dissociation–reassociation equilibrium between p and the borane (Fig. 8). Accordingly, all p-borane atropisomers are formed in different chemical reactions between the same Lewis acid (I, II or III) and base (p) which directly proves the concept of MLA for these triaryl-boranes. To further corroborate this finding, a series of DFT potential energy surface scans were performed using the p–II(α) and II structures (ESI Section 8†). These showed that in contrast with the sp2 borane state (12–14 kcal mol−1), rotations around the central B(sp3) atom (i.e. B–C and B–N bonds) are hampered in the dative bonded p-adduct state with the associated barrier heights ranging between 25 and 53 kcal mol−1. These values are also significantly higher than the experimentally measured barrier heights (17–20 kcal mol−1) for the atropisomeric exchanges of the p–II adduct (Table S14†) which demote the interconversion of the p-adduct atropisomeric states via internal rotations in comparison to the breaking of the dative B–N bond. Hence, in the case of these dative adducts, the non-stable dative bonding underlies the conformational dynamic equilibrium between their respective atropisomers via breaking and reforming the dative bond and it seems that these atropisomers have a rather rigid structure, ultimately adhering to the nomenclature ‘atropos’, i.e. ‘without turn’.
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 2004799 and 2004800. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc00925h |
| This journal is © The Royal Society of Chemistry 2024 |