
==== Front
Chem Sci
Chem Sci
SC
CSHCBM
Chemical Science
2041-6520
2041-6539
The Royal Society of Chemistry

39239481
d4sc03699a
10.1039/d4sc03699a
Chemistry
Correlated shapeshifting and configurational isomerization†
† Electronic supplementary information (ESI) available. CCDC 2294194 (E,γ-1a). For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc03699a

Hussein Burhan A. a
Maturi William ab
Rylands Mary Kate a
https://orcid.org/0000-0002-0588-4939
Bismillah Aisha N. ab
Wen Yuzhen ab
https://orcid.org/0000-0001-9181-8892
B-5315-2015
Aguilar Juan A. a
https://orcid.org/0000-0003-2128-6733
Ayub Rabia b
https://orcid.org/0000-0002-7104-847X
Rankine Conor D. b
https://orcid.org/0000-0002-8538-7579
J-8291-2012
McGonigal Paul R. ab
a Department of Chemistry, Durham University, Lower Mountjoy Stockton Road Durham DH1 3LE UK paul.mcgonigal@york.ac.uk

b Department of Chemistry, University of York Heslington York YO10 5DD UK
23 8 2024
18 9 2024
23 8 2024
15 36 1461814624
6 6 2024
23 8 2024
This journal is © The Royal Society of Chemistry
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by/3.0/ This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given.
Herein we demonstrate that the rapid ‘shapeshifting’ constitutional isomerization of a substituted bullvalene is influenced by the E-to-Z configurational isomerization of a remote carbamate group, giving rise to correlated motion. We find that, while the E-configurational isomer of a bulky carbamate favors the β-bullvalene constitutional isomer, a noncovalent bonding interaction within the Z-carbamate tips the equilibrium toward the γ-bullvalene form. Using DFT modelling and NMR spectroscopy, this long-range interaction is identified as being between the bullvalene core and a pendant phenyl group connected to the carbamate. Coupling the constitutional changes of a bullvalene to a reciprocal configurational isomerization through a long-range interaction in this way will allow shapeshifting rearrangements to be exploited as part of collective motion in extended structures.

Rearrangements of a bullvalene cage between constitutional isomers are directed by configurational isomerization of a pendant group through noncovalent interactions.

Engineering and Physical Sciences Research Council 10.13039/501100000266 EP/L015536/1 EP/T518001/1 EP/V040049/2 EP/V047817/2 EP/X021564/1 Leverhulme Trust 10.13039/501100000275 RPG-2020-218 pubstatusPaginated Article
==== Body
pmcIntroduction

The rapid, successive Cope rearrangements of bullvalenes1 produce a series (Fig. 1a) of constitutional isomers. Recent advances in the synthesis of substituted derivatives,2 alongside the development of related fluxional cages,3 have led to renewed interest in exploiting these ‘shapeshifting’ structures as part of functional molecules and materials. Part of the appeal of using bullvalene derivatives is that they juxtapose the rigidity and well-defined bond angles1e that are typical of cage-like structures with the rapid structural dynamics that are commonly only found in flexible molecules. They exhibit ‘rigid dynamics’ at the single-molecule level.4 Accordingly, investigations have been reported into bullvalene-containing fluorophores,5 antibiotics,6 small-molecule receptors,7 transition metal complexes,8 rigid-rod polymers,9 and single-molecule junctions.10

Fig. 1 The combination of (a) bullvalene shapeshifting rearrangements between constitutional isomers with (b) conformational and configurational changes by single bond rotations and E-to-Z isomerization are (c) correlated through noncovalent interactions (red) to bias the two bullvalene rearrangement equilibria in opposite directions (blue arrows). Black and purple circles represent appended functional groups.

By contrast, most organic structures, particularly those containing rotatable single bonds, rapidly sample stereoisomeric structures through changes in bond angles and dihedral angles.11 The resulting isomerization (Fig. 1b) is not only influenced by local electronic and steric factors, but also by long-range interactions. Bond rotations of groups that are separated by several covalent bonds but are close together in space can, therefore, become correlated with one another,12i.e., the isomerization or motion of one part of the molecule is coupled to a reciprocal conformational or configurational change in another part.

To exploit the shapeshifting rearrangements of bullvalenes in larger collective motions of extended structures,12d–f it is necessary to understand how their fluxional Cope rearrangements are influenced (Fig. 1c) by isomerization of surrounding groups. Here, we report the isomeric distribution of carbamate-functionalized bullvalenes 1 (Scheme 1) and demonstrate that, even in this structurally simple case, conformational changes can become correlated to shapeshifting isomerization. For one of the compounds (benzhydryl derivative 1b), we find that while the E-carbamate is most stable as its β-substituted bullvalene isomer, bond rotation to the Z-carbamate biases the bullvalene unit towards its γ-substituted form instead, i.e., the energetics of bullvalene isomerization and carbamate rotation are coupled to one another. We show that through-space interactions subtly alter the energetics of the dynamic system. As part of this investigation, we have also assessed the extent to which this isomer distribution can be accurately modelled by comparing the calculated energies of the lowest energy conformers, as opposed to analyzing the full conformational landscape.

Scheme 1 (a) Synthesis of sterically crowded bullvalene carbamates. Reagents and conditions: (i) 2, Boc2O, NEt3, CH2Cl2, rt, 24 h; (ii) 1, NaH, DMF, 0 °C, 30 min, 2. propargyl bromide, rt, 5 h, 70% 3a from 2a, 64% 3b from 2b; (iii) 3, cyclooctatetraene, CoBr2(dppe), ZnI2, Zn, TFE, 55 °C, 18 h, 55% 4a, 50% 4b; (iv) 4, hν (365 nm), thioxanthone, THF, 25 °C, 3 h, 50% 1a, 54% 1b. DMF = dimethylformamide, dppe = 1,2-bis(diphenylphosphino)ethane, TFE = 2,2,2-trifluoroethanol, and THF = tetrahydrofuran. (b) Eight isomers interchange by Cope rearrangements (horizontal equilibria) and by rotation around the carbamate N–C bond (vertical equilibria).

Results and discussion

We targeted (Scheme 1) tertiary carbamates 1a and 1b derived from benzylamine 2a and benzhydrylamine 2b, respectively. The partial double bond character of the carbamate C–N bond (c.f., tertiary amide C–N bonds13) gives rise to E-to-Z configurational isomerization that can be slowed down and observed by dynamic NMR (DNMR) spectroscopy at low temperatures, which we reasoned would allow us to experimentally measure distributions of the configurational and constitutional isomers of 1a and 1b.

Starting from 2, we performed (Scheme 1a) a sequence of carbamate formation, alkylation with propargyl bromide, and cobalt(i)-catalyzed [6 + 2] cycloaddition with cyclooctatetraene2c,14 to produce bicyclo[4.2.2]deca-2,4,7,9-tetraene intermediates 4. A final photochemical di-π-methane rearrangement step15 using thioxanthone2e as a photosensitizer gave rise to the target bullvalenes 1.

There are four possible constitutional isomers of 1 (Scheme 1b), which we label as α–δ according to the attachment point of the carbamate substituent to the bullvalene cage.1e Consecutive Cope rearrangements from the α-isomer lead to the δ-, β-, then γ-isomers in a reversible linear sequence. Alternatively, rotation around the carbamate C–N bond interconverts the E- and Z-configurational isomers. In combination, these two pathways give rise to a set of eight isomers (Scheme 1b).

1H NMR Spectroscopic analysis of the benzyl derivative 1a in CDCl3 at 313 K (Fig. 2a) shows two broad resonances corresponding to the two sets of methylene protons, H1 and H2. Rapid chemical exchange at this temperature averages out the contributions from the different isomers to these two methylene resonances and causes the signals arising from the bullvalene methine groups, Hα–Hδ, to broaden and partially merge. At 219 K, the reduced rates of Cope rearrangements and carbamate rotation allow signals from the different species to be resolved.

Fig. 2 1H DNMR (499 MHz, CDCl3) and crystallographic analyses of 1a. (a) Partial 1H NMR spectra at temperatures ranging from 219 K to 313 K. Signals are assigned using the labels in Scheme 1. Regions in orange boxes correspond to those shown in panels b, d, and e. (b) Line fitting was used to deconvolute the H1 signals at 219 K. An overlay of the experimentally measured (black) and simulated (orange) traces is shown, including the integrals of the four simulated peaks (blue). (c) X-ray crystal structure of E,γ-1a. (d) Partial EXSY NMR spectrum (499 MHz, CDCl3, tm = 500 ms). Blue boxes indicate cross peaks arising from Cope rearrangements. (e) Line-shape analysis for the H2 signal at 308 K, comparing the experimentally measured (exp.) and simulated (sim.) line shapes. See Table S1† for simulation parameters.

Although there are eight possible isomers, the α- and δ-forms of mono-substituted bullvalenes are usually several kJ mol−1 higher in energy than the olefin-substituted β- and γ-forms,16 so they are not present in sufficient concentrations to be observed. β-Substitution is usually preferred over γ-substitution by ∼1–2 kJ mol−1.2c In keeping with this expectation, line fitting of the H1 peaks indicates (Fig. 1b) that four isomers are present in a 45 : 23 : 21 : 11 ratio, which can be assigned to Z,β-1a, E,β-1a, Z,γ-1a, and E,γ-1a, respectively, using 2D NMR (Fig. S15–S22†). Shape-selective crystallization‡2e,17 from this mixture yields E,γ-1a (the least populated of the four isomers in solution) as the sole isomer in the solid state (Fig. 2c). In energetic terms (Table 1), the solution-phase population corresponds to small Gibbs energy biases ΔGexp of ∼1.3 kJ mol−1 toward both β-isomers over their γ-substituted counterparts, and ΔGexp of ∼1 kJ mol−1 for each Z-carbamate over its corresponding E-carbamate. Therefore, for 1a, there is no evidence that the E/Z-state of the molecule substantially changes the energetics of β/γ-isomerization, or vice versa. The two isomerization processes do not appear to be correlated.

Populations and relative Gibbs energies of the isomers present in the dynamic mixtures and their transition states

State	1a, p/%	1a, ΔG/kJ mol−1	1b, p/%	1b, ΔG/kJ mol−1	
p exp a	p calc1 c	p calcG d	ΔGexpa	ΔGcalc1f	ΔGcalcGh	p exp a	p calc1 c	p calcG d	ΔGexpa	ΔGcalc1f	ΔGcalcGh	
E,α	—b	<0.1	<0.1	—b	25.5	25.4	—b	<0.1	<0.1	—b	34.5	34.0	
Z,α	—b	<0.1	<0.1	—b	23.1	23.2	—b	<0.1	<0.1	—b	26.5	26.8	
E,β	23	18	22	1.2	1.7	1.4	16	24	23	1.9	1.4	1.2	
Z,β	45	47	46	0.0	0.0	1.0	27	18	28	1.0	1.5	0.8	
E,γ	11	10	13	2.5	2.8	2.3	8	6	5	3.2	3.8	3.8	
Z,γ	21	14	18	1.4	2.2	1.7	48	51	44	0.0	0.0	0.0	
E,δ	—b	<1	<1	—b	8.7	8.2	—b	<0.1	<0.1	—b	14.8	14.1	
Z,δ	—b	<1	1	—b	7.2	6.8	—b	<1	<1	—b	9.5	8.8	
TS1	—	—	—	55.6	—g	—g	—	—	—	54.5	—g	—g	
TS2	—	—	—	57.2	—g	—g	—	—	—	57.4	—g	—g	
TS3	—	—	—	64.9e	—g	—g	—	—	—	54.8	—g	—g	
TS4	—	—	—	—	—g	—g	—	—	—	55.3	—g	—g	
a Measured at 219 K in CDCl3.

b The population of α- and δ-isomers is below the detection limit of the equipment used to acquire 1H NMR spectra.

c Population calculated as a Boltzmann distribution based on the DFT calculated energies of the most stable conformer for each constitutional and configurational isomer (ΔGcalc1).

d Population calculated as a global Boltzmann distribution by summing the populations of all accessible conformers calculated for each constitutional and configurational isomer (ΔGcalcG).

e Apparent Gibbs energy barrier for carbamate rotation (weighted average of TS3 and TS4) at 308 K.

f Relative Gibbs energies of only the most stable conformer for each isomer calculated at 219 K using PBE0-D3/def2-SV(P).20i–k,p

g Transition states were not modelled.

h Effective relative Gibbs energies based on the DFT calculated global Boltzmann distribution (pcalcG) for the purpose of comparison with ΔGexp.

This system does, however, provide an unusual case of reversible covalent changes occurring in a molecule more rapidly than its bond rotation, i.e., the carbamate isomerization of 1a proceeds at a slower rate than its Cope rearrangement. Using 1H–1H Exchange NMR Spectroscopy18 (EXSY) at 219 K (Fig. 2d), we found that the Cope rearrangements pass through transition states (TS1 and TS2, Scheme 1b) that are ∼55–57 kJ mol−1 above the lowest energy isomer, Z,β-1a, which is typical19 for bullvalene constitutional isomerization. There are no cross peaks visible between E- and Z-isomers at this temperature because of the slow carbamate rotation rate. At the higher temperature of 308 K, 1H NMR line-shape analysis (Fig. 2e) gives an averaged rate of carbamate rotation, kEZ, of 62 s−1, indicating that TS3 and TS4 lie almost 10 kJ mol−1 higher than TS1 at ΔGexp of ∼65 kJ mol−1 relative to Z,β-1a (Table 1).

Unlike 1a, the increased steric crowding present in the benzhydryl derivative 1b leads to correlated isomerization. The additional phenyl ring of 1b experiences long-range interactions with the bullvalene that influence the solution-phase isomerization equilibria. There are some key clues to this behavior in the variable-temperature NMR spectra (Fig. 3a). Crucially, the distribution of H1 methylene resonances (Fig. 3b) has been altered substantially relative to 1a. The Z,β-1b, E,β-1b, Z,γ-1b, and E,γ-1b isomers are present in a 27 : 16 : 48 : 8 ratio, as assigned by 2D NMR (Fig. S15–S22†). Therefore, the Z,γ-isomer is the most favored form of 1b, despite being only the third most populated isomer of 1a. The changes in relative Gibbs energy (Grel) between the 1a and 1b isomers is summarized graphically as a potential energy diagram in Fig. 4. A structural change from a hydrogen substituent to a phenyl group at a remote site has overridden the inherent preference2c,16 for the bullvalene β-isomer in the reaction network, but only when the carbamate is in its Z-form (Fig. 5a). The benzhydryl carbamate of 1b also rotates at a faster rate than the benzyl carbamate of 1a, exhibiting near-identical energy barriers of ∼54–58 kJ mol−1 for the four isomerization processes (Table 1).

Fig. 3 1H DNMR analysis (499 MHz, CDCl3) of 1b. (a) Partial 1H NMR spectra at temperatures ranging from 219 K to 313 K. Signals are assigned using the labels in Scheme 1. The region in an orange box corresponds to that shown in panels b and c. (b) Line fitting was used to deconvolute the H1 signals at 219 K. An overlay of the experimentally measured (black) and simulated (orange) traces is shown, including the integrals of the four simulated peaks (blue). (c) Partial EXSY NMR spectrum (499 MHz, CDCl3, tm = 500 ms). Blue boxes indicate cross peaks arising from Cope rearrangements. Dashed black boxes indicate cross peaks arising from carbamate rotation.

Fig. 4 Potential energy diagrams for the observed isomers of (a) 1a and (b) 1b, showing the relative stabilization of Z,γ-1b. Transition states in blue correspond to Cope rearrangement steps, those in purple correspond to E/Z-isomerisation.

Fig. 5 (a) Network diagrams for the isomerization of 1a and 1b showing the experimentally measured populations (pexp) of the β- and γ-isomers at 219 K as filled grey circles and the calculated populations (pcalG) as hollow orange circles. According to DFT calculations, the α- and δ-isomers are present in <1%. Blue arrows point to the major species of each Cope rearrangement process, highlighting that the preference for the β- and γ-bullvalene isomers is correlated to the E-to-Z isomerization of 1b. (b) NCI plots [LC-PBE/def2-TZVP]22 of Z,γ-1b and Z,β-1b showing the increased Ph–BV interaction in Z,γ-1b.

To understand the unusual constitutional isomer distribution of 1b, we carried out DFT modelling. We generated all the constitutional isomers of 1a and 1b and carried out distance-geometry-based generation of their configurational and conformational isomers before optimizing the geometries at the PBE0-D3/def2-SV(P) level of theory.1e,20 As there are several rotatable single bonds in the structure of 1, each of the isomers observed by low-temperature NMR itself represents a rapidly interconverting population of conformers. The relative energies of all the isomers are tabulated in Tables S3 and S4.† Cartesian coordinates are available in the ESI.† To take one representative example, there are six conformers of Z,β-1a spanning a range of ∼16 kJ mol−1 in energy, of which the two lowest energy structures are within just 2.5 kJ mol−1 of one another.

Despite there being numerous different structures contributing to the overall energetics of E-to-Z and β-to-γ isomerization, the experimentally observed populations, pexp, can be approximated reasonably well by analyzing just the lowest energy conformers (Table 1). For 1a, the calculated energies, ΔGcalc1, of the lowest energy conformer for each of the four observable isomers gives a Boltzmann distribution, pcalc1, that matches closely to pexp. For 1b, however, this simplified approach incorrectly predicts that more of E,β-1b is present than Z,β-1b. Instead, it is necessary to perform a global Boltzmann population analysis (giving pcalcG) that takes into account all the accessible conformers in order to reproduce the experimentally observed hierarchy of isomers. This observation indicates that for bullvalene derivatives with several rotatable single bonds, accurate modelling of the shapeshifting rearrangements requires consideration of the full energetic landscape of accessible conformers.

Noncovalent interaction (NCI) plots21 generated using the DFT-optimized geometries reveal (Fig. 5b and S26–S29†) the interactions that are responsible for the correlated isomerization of 1b. In addition to several close contacts between the crowded tert-butyl, carbonyl, and phenyl groups, which are present in all of the isomers, the lowest energy conformer of Z,γ-1b shows evidence of substantial interaction between an ortho-hydrogen of one phenyl ring with the bullvalene cyclopropyl ring (Ph–BV interaction, Fig. 5b). This long-range interaction is less significant in the NCI analyses of the other 1b isomers (Fig. S26–S29†), such as Z,β-1b (Fig. 5b). Experimental evidence for the interaction in Z,γ-1b is apparent in its low-temperature 1H NMR spectrum (Fig. 3a) as its Hδ resonance is shifted downfield relative to the Hδ signals of the other isomers of 1b, which is indicative of deshielding caused by the aromatic ring current. The NOESY NMR spectrum (Fig. S23†) also confirms Hδ undergoes through-space interaction with a phenyl group. Although this Ph–BV interaction appears to be a relatively weak van der Waals contact, the combination of it together with other subtle differences in the attractive and repulsive noncovalent bonding interactions present in 1b is evidently sufficient to provide enough of an energetic bias to overcome the ∼1–2 kJ mol−1 preference for the β-isomer.

Conclusions

In summary, the shapeshifting constitutional rearrangements of even structurally simple bullvalenes can be influenced by a complementary, remote isomerization process. Given that the energetic differences between bullvalene isomers are often as little as a few kJ mol−1, relatively subtle noncovalent bonding interactions, including weak van der Waals contacts, are sufficient to reshuffle the proportions of each isomer. In the system presented here, the introduction of an additional phenyl ring to the structure of 1a causes the Z,γ-isomer to change from being the third most populous to being the major species of 1b as a result of a Ph–BV interaction. Programming in switchable long-range interactions may allow for control to be exercised over the large number of constitutional isomers that arise in multi-substituted bullvalene derivatives.

Data availability

Experimental synthetic procedures, characterization data and theoretical calculation results are available in the ESI.†

Author contributions

Conceptualization: BAH, PRM. Investigation: BAH, WM, MKR, ANB, JAA, CR. Writing: BAH, WM, ANB, YW, RA, PRM.

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

SC-015-D4SC03699A-s001

SC-015-D4SC03699A-s002

SC-015-D4SC03699A-s003

B. A. H. acknowledges an EPSRC SOFI CDT Studentship (EP/L015536/1) and an NSERC Scholarship. W. M. (EP/T518001/1), R. A (EP/X021564/1) and P. R. M (EP/V047817/2, EP/V040049/2) thank the EPSRC for funding. A. N. B. and P. R. M. acknowledge a Leverhulme Trust Research Project Grant (RPG-2020-218). Y. W. thanks the Chinese Scholarship Council. We thank Dr Dmitry Yufit for performing X-ray crystallography. We acknowledge the University of York High Performance Computing service, Viking.

‡ The energetics of crystal packing can outweigh the small solution-phase energetic bias towards certain isomers within fluxional mixtures.2e,17
==== Refs
Notes and references

(a) Doering W. v. E. Roth W. R. Tetrahedron 1963 19 715 737 10.1016/S0040-4020(01)99207-5
(b) Schröder G. Angew. Chem., Int. Ed. Engl. 1963 2 481 482 10.1002/anie.196304814
(c) Ferrer S. Echavarren A. Synthesis 2019 51 1037 1048 10.1055/s-0037-1611637
(d) Bismillah A. N. Chapin B. M. Hussein B. A. McGonigal P. R. Chem. Sci. 2020 11 324 332 10.1039/C9SC05482K 32206269
(e) Ives R. A. Maturi W. Gill M. T. Rankine C. D. McGonigal P. R. Chem. Sci. 2024 10.1039/D4SC03700F
(a) Lippert A. R. Kaeobamrung J. Bode J. W. J. Am. Chem. Soc. 2006 128 14738 14739 10.1021/ja063900+ 17105247
(b) Ferrer S. Echavarren A. M. Angew. Chem., Int. Ed. 2016 55 11178 11182 10.1002/anie.201606101 27538835
(c) Yahiaoui O. Pašteka L. F. Judeel B. Fallon T. Angew. Chem., Int. Ed. 2018 57 2570 2574 10.1002/anie.201712157 29314602
(d) Yahiaoui O. Pašteka L. F. Blake C. J. Newton C. G. Fallon T. Org. Lett. 2019 21 9574 9578 10.1021/acs.orglett.9b03737 31746207
(e) Patel H. D. Tran T. H. Sumby C. J. Pašteka L. F. Fallon T. J. Am. Chem. Soc. 2020 142 3680 3685 10.1021/jacs.9b12930 32041409
(f) Sanchez A. Gonzalez V. M. Sakamoto J. Gurajapu A. Maimone T. J. J. Am. Chem. Soc. 2024 146 17573 17579 10.1021/jacs.4c03323 38901002
(a) McGonigal P. R. de León C. Wang Y. Homs A. Solorio-Alvarado C. R. Echavarren A. M. Angew. Chem., Int. Ed. 2012 51 13093 13096 10.1002/anie.201207682
(b) Wang X. Shyshov O. Hanževački M. Jäger C. M. von Delius M. J. Am. Chem. Soc. 2019 141 8868 8876 10.1021/jacs.9b01350 31117548
(c) Sanz-Novo M. Mato M. León Í. Echavarren A. M. Alonso J. L. Angew. Chem., Int. Ed. 2022 61 e202117045 10.1002/anie.202117045
(d) Sanchez A. Maimone T. J. J. Am. Chem. Soc. 2022 144 7594 7599 10.1021/jacs.2c02627 35420799
(e) Sanchez A. Gurajapu A. Guo W. Kong W.-Y. Laconsay C. J. Settineri N. S. Tantillo D. J. Maimone T. J. J. Am. Chem. Soc. 2023 145 13452 13461 10.1021/jacs.3c03960 37279177
(f) Bismillah A. N. Johnson T. G. Hussein B. A. Turley A. T. Saha P. K. Wong H. C. Aguilar J. A. Yufit D. S. McGonigal P. R. Nat. Chem. 2023 15 615 624 10.1038/s41557-023-01156-7 36914791
(a) Deng H. Olson M. A. Stoddart J. F. Yaghi O. M. Nat. Chem. 2010 2 439 443 10.1038/nchem.654 20489710
(b) Krause S. Milić J. V. Commun. Chem. 2023 6 151 10.1038/s42004-023-00945-y 37452112
(c) Saha P. K. Tran Ngoc T. McGonigal P. R. Teichert J. F. Nat. Synth. 2024 3 684 697 10.1038/s44160-024-00526-4
Dohmen C. Ihmels H. Paululat T. Eur. J. Org Chem. 2022 2022 e202201172 10.1002/ejoc.202201172
Ottonello A. Wyllie J. A. Yahiaoui O. Sun S. Koelln R. A. Homer J. A. Johnson R. M. Murray E. Williams P. Bolla J. R. Robinson C. V. Fallon T. Soares da Costa T. P. Moses J. E. Proc. Natl. Acad. Sci. U. S. A. 2023 120 e2208737120 10.1073/pnas.2208737120 37011186
(a) Larson K. K. He M. Teichert J. F. Naganawa A. Bode J. W. Chem. Sci. 2012 3 1825 10.1039/C2SC20238G
(b) Teichert J. F. Mazunin D. Bode J. W. J. Am. Chem. Soc. 2013 135 11314 11321 10.1021/ja404981q 23802800
(a) Birve A. Patel H. D. Price J. R. Bloch W. M. Fallon T. Angew. Chem., Int. Ed. 2021 61 e202115468 10.1002/anie.202115468
(b) Dohmen C. Paululat T. Ihmels H. Chem.–Eur. J. 2024 30 e202304311 10.1002/chem.202304311 38275100
(a) Pomfret M. N. Sun P. B. Huang Z. Freund A. C. Miyoshi T. Golder M. R. Angew. Chem., Int. Ed. 2023 62 e202301695 10.1002/anie.202301695
(b) Sun P. B. Pomfret M. N. Elardo M. J. Suresh A. Rentería-Gómez Á. Lalisse R. F. Keating S. Chen C. Shayna Hilburg L. Chakma P. Wu Y. Bell R. C. Rowan S. J. Golder M. R. J. Am. Chem. Soc. 2024 146 19229 19238 10.1021/jacs.4c04401 38961828
Reimers J. R. Li T. Birvé A. P. Yang L. Aragonès A. C. Fallon T. Kosov D. S. Darwish N. Nat. Commun. 2023 14 6089 10.1038/s41467-023-41674-z 37789027
Mazzanti A. Casarini D. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012 2 613 641
(a) Mislow K. Acc. Chem. Res. 1976 9 26 33 10.1021/ar50097a005
(b) Oki M. Angew. Chem., Int. Ed. Engl. 1976 15 87 93 10.1002/anie.197600871
(c) Rappoport Z. Biali S. E. Acc. Chem. Res. 1997 30 307 314 10.1021/ar960216z
(d) Erbas-Cakmak S. Leigh D. A. McTernan C. T. Nussbaumer A. L. Chem. Rev. 2015 115 10081 10206 10.1021/acs.chemrev.5b00146 26346838
(e) Dattler D. Fuks G. Heiser J. Moulin E. Perrot A. Yao X. Giuseppone N. Chem. Rev. 2020 120 310 433 10.1021/acs.chemrev.9b00288 31869214
(f) Corra S. de Vet C. Baroncini M. Credi A. Silvi S. Chem 2021 7 2137 2150 10.1016/j.chempr.2021.04.010 34435161
(g) Van Beek C. L. F. Feringa B. L. J. Am. Chem. Soc. 2024 146 5634 5642 10.1021/jacs.3c14430 38350104
(a) Johnston E. R. Fortt R. Barborak J. C. Magn. Reson. Chem. 2000 38 932 936 10.1002/1097-458X(200011)38:11<932::AID-MRC752>3.0.CO;2-6
(b) Clayden J. Pink J. H. Angew. Chem., Int. Ed. 1998 37 1937 1939 10.1002/(SICI)1521-3773(19980803)37:13/14<1937::AID-ANIE1937>3.0.CO;2-4
(c) Bragg R. A. Clayden J. Org. Lett. 2000 2 3351 3354 10.1021/ol0064462 11029208
(d) Bragg R. A. Clayden J. Morris G. A. Pink J. H. Chem.–Eur. J. 2002 8 1279 1289 10.1002/1521-3765(20020315)8:6<1279::AID-CHEM1279>3.0.CO;2-7 11921211
Achard M. Mosrin M. Tenaglia A. Buono G. J. Org. Chem. 2006 71 2907 2910 10.1021/jo052630v 16555854
(a) Jones M. Scott L. T. J. Am. Chem. Soc. 1967 89 150 151 10.1021/ja00977a031
(b) Jones M. Reich S. D. Scott L. T. J. Am. Chem. Soc. 1970 92 3118 3126 10.1021/ja00713a033
Oth J. F. M. Merényi R. Nielsen J. Schröder G. Chem. Ber. 1965 98 3385 3400 10.1002/cber.19650981040
(a) Bismillah A. N. Sturala J. Chapin B. M. Yufit D. S. Hodgkinson P. McGonigal P. R. Chem. Sci. 2018 9 8631 8636 10.1039/C8SC04303E 30647882
(b) Nyman J. Day G. M. CrystEngComm 2015 17 5154 5165 10.1039/C5CE00045A
Nikitin K. O'Gara R. Chem.–Eur. J. 2019 25 4551 4589 10.1002/chem.201804123 30421834
Merényi R. Oth J. F. M. Schröder G. Chem. Ber. 1964 97 3150 3161 10.1002/cber.19640971126
(a) bullviso, 2024, https://gitlab.com/conorrankine/bullviso
(b) Riniker S. Landrum G. A. J. Chem. Inf. Model. 2015 55 2562 2574 10.1021/acs.jcim.5b00654 26575315
(c) Wang S. Witek J. Landrum G. A. Riniker S. J. Chem. Inf. Model. 2020 60 2044 2058 10.1021/acs.jcim.0c00025 32155061
(d) RDKit: Open-Source Cheminformatics. https://rdkit.org
(e) RDKit, 2023, https://github.com/rdkit/rdkit
(f) Grimme S. Bannwarth C. Shushkov P. J. Chem. Theory Comput. 2017 13 1989 2009 10.1021/acs.jctc.7b00118 28418654
(g) Bannwarth C. Ehlert S. Grimme S. J. Chem. Theory Comput. 2019 15 1652 1671 10.1021/acs.jctc.8b01176 30741547
(h) Bannwarth C. Caldeweyher E. Ehlert S. Hansen A. Pracht P. Seibert J. Spicher S. Grimme S. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2020 11 e01493
(i) Perdew J. P. Burke K. Ernzerhof M. Phys. Rev. Lett. 1996 77 3865 3868 10.1103/PhysRevLett.77.3865 10062328
(j) Perdew J. P. Burke K. Ernzerhof M. Phys. Rev. Lett. 1997 78 1396 10.1103/PhysRevLett.78.1396
(k) Adamo C. Barone V. J. Chem. Phys. 1999 110 6158 6169 10.1063/1.478522
(l) Grimme S. Antony J. Ehrlich S. Krieg H. J. Chem. Phys. 2010 132 154104 10.1063/1.3382344 20423165
(m) Neese F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012 2 73 78
(n) Neese F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2017 8 e1327
(o) Neese F. Wennmohs F. Becker U. Riplinger C. J. Chem. Phys. 2020 152 224108 10.1063/5.0004608 32534543
(p) Weigend F. Ahlrichs R. Phys. Chem. Chem. Phys. 2005 7 3297 3305 10.1039/B508541A 16240044
(q) Weigend F. Phys. Chem. Chem. Phys. 2006 8 1057 1065 10.1039/B515623H 16633586
(a) Iikura H. Tsuneda T. Yanai T. Hirao K. J. Chem. Phys. 2001 115 3540 3544 10.1063/1.1383587
(b) Najibi A. Casanova-Paez M. Goerigk L. J. Phys. Chem. A 2021 125 4026 4035 10.1021/acs.jpca.1c02549 33938224
(a) Johnson E. R. Keinan S. Mori-Sánchez P. Contreras-García J. Cohen A. J. Yang W. J. Am. Chem. Soc. 2010 132 6498 6506 10.1021/ja100936w 20394428
(b) Contreras-García J. Johnson E. R. Keinan S. Chaudret R. Piquemal J.-P. Beratan D. N. Yang W. J. Chem. Theory Comput. 2011 7 625 632 10.1021/ct100641a 21516178
