Supramolecular self-sorting predicted by a simple harmonic oscillator model

This journal is © The Royal Society of Chemistry 2024 Chem. Commun., 2024, 60, 14109–14112 | 14109 Cite this: Chem. Commun., 2024, 60, 14109 Supramolecular self-sorting predicted by a simple harmonic oscillator model† Adriana Sacrista ´n-Martı ´n, ‡ a Nerea A ´lvarez-Llorente, ‡ b Eric Masson, * c Alberto Diez-Varga * b and He ´ctor Barbero * b Macrocycles that encapsulate two guests can self-sort those into homo- and heterodimers. We report here a family of self-sorting homobimetallic Pt(II) terpyridyl acetylide dimers secured together with a pair of Cucurbit[8]uril macrocycles (CB[8]). The rigid bridging unit between both Pt centers introduces varying ‘‘hinge’’ angles, resulting in disparities in Pt–Pt distances in the heterodimers, and leads to recognition motif mismatch. We found that the self-sorting process can be quantified using a simple model, in which each complex behaves as a simple harmonic oscillator, whose heteroassembly tends to minimize geometry distorsions through C(aryl)–CR R RC–Pt axis deformation. Cucurbit[n]urils (CB[n]s) are barrel-shaped macrocycles consisting of nglycoluril units connected via methylene bridges. Their capacity to form strong inclusion complexes in aqueous media, primarily due to the release of energetically constrained water molecules, has established them as key components for selfsorting in a wide variety of chemical systems. 1–7 CB[8] is capable of bringing together up to two hydrophobic guests, allowing homo/ hetero pair formation. 7–10 In those tertiary assemblies, cationic units on the guest tend to minimize electrostatic repulsion with the other positive guest, while maximizing attractive interaction with both carbonyl rims, resulting in head-to-tail (HT) assemblies. Nonetheless, these two driving forces might be overcome by employing square planar Pt(II)-terpyridine (tpy) complexes, wherein both extended dispersive ligand attractive forces and metallophilic interactions cooperate to exclusively afford head-to-head (HH) assemblies. 11 This new arrangement has led to several applications in peptide engineering 12,13 and hydrogen photoreduction. 14 Homobimetallic Pt(II) acetylide complexes form quaternary assemblies upon encapsulation with two CB[8] hosts (Scheme 1). The self-sorting process can be quantified using equilibrium (1) 15 and the corresponding free energy term DG (eqn (2)),§where A and B are guests, AA and BB homoquaternary complexes, AB (or BA) hetero-quaternary complexes and [AB] tot is the total concentration of hetero-complex. 1 2AA þ1 2BB$1 2AB þ1 2BA (1) DG¼RTln ½ABtot 2½AA1 2½BB1 2 (2) The equilibrium is considered fully social when the homodimers cannot be detected by NMR spectroscopy (i.e. less than 1%); Scheme 1 Previously reported equilibrium between two homoquaternary inclusion complexes (HH) with tolyl and 3,5-difluorophenyl heads as recognition motifs. 16 a Institute of Organic Chemistry, Ulm University, 89081 Ulm, Germany b GIR MIOMeT, IU CINQUIMA/Quı ´mica Inorga ´nica, Facultad de Ciencias, Universidad de Valladolid, Valladolid, E47011, Spain. E-mail: [email protected], [email protected] c Department of Chemistry and Biochemistry, Ohio University, Athens, Ohio 45701, USA. E-mail: mass[email protected] †Electronic supplementary information (ESI) available: Detailed synthetic methods and characterization. Self-sorting process study. Computational data. See DOI: https://doi.org/10.1039/d4cc05336b ‡These two authors contributed equally to this work. Received 9th October 2024, Accepted 25th October 2024 DOI: 10.1039/d4cc05336b rsc.li/chemcomm ChemComm COMMUNICATION Open Access Article. Published on 28 October 2024. Downloaded on 2/17/2025 8:23:37 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online View Journal | View Issue 14110 | Chem. Commun., 2024, 60, 14109–14112 This journal is © The Royal Society of Chemistry 2024 this corresponds to a 2.3 kcal mol 1 free energy term. Similarly, self-sorting is fully narcissistic when DG4+2.3 kcal mol 1 .Inthe case of bimetallic Pt(II) acetylides, self-sorting properties primarily depend on (1) the relative binding free energies of the aryl head groups; and (2) the relative orientation of the bridging group. 16 In order to isolate the second contribution, we use here rigid spacers to assess the importance of structural flexibility on self-sorting. A family of five homobimetallic Pt(II)-tpy complexes 1–5 bearing distinct aromatic bridging ligands have been synthesized (Fig. 1). These species are expected to form well-defined homoquaternary assemblies whose self-sorting will be evaluated. To avoid the contribution of the aryl heads, a uniform tpy head group (3,5-difluophenyl) was utilized (see highlighted R group in Scheme1);self-sortingisthussolelygovernedbyinteractionsand angle mismatches between the central bridging units. The angle of bridging ligands (that we call ‘‘hinge’’ angle) is defined by two vectors extending from the quaternary carbons bearing the alkynyl moiety and the opposite quaternary (para) carbons. The range spans 64 degrees (see Fig. 1).¶ The distance between both substituted carbons in the spacer group is similar so that a good overlap is expected, with the exception of the 2,7-substituted naphthyl unit 16 (complex 5,Fig.1),whoselengthissignificantly shorter than any other ligand. This spacer group has been studied to account for ligand length differences when the angle disparity is small (61)incombinationwithcomplex4.Inanymixture,atpy head group mismatch (i.e. aryl units from distinct complexes that are misaligned within the CB[8] cavity) is therefore expected. This results in the deformation of both complexes in the heteroquaternary assembly to align both pairs of head groups; the deformation would occur at the most flexible part of the molecule, i.e. the acetylide moiety. 7,8,17,18 Hence, experimental free energy terms on our self-sorting scale 15 provide a measurement of Ptalkynyl flexibility. Pt(II)-tpy complexes 1–4 were prepared from the corresponding diacetylide ligands, chloro-[40-(3,5-difluoro phenyl)- 2,20:60,200-terpyridine]platinum(II) chloride, copper(I)iodideand isopropyl amine.§These complexes readily formed HH homoquaternary assemblies 1 2 CB[8] 2 –4 2 CB[8] 2 upon addition of 1.0 equivalent of CB[8] to their suspension in a 2: 1 mixture of D 2 O and CD 3 CN (Fig. S91–S114, ESI†).§Pairs of homoquaternary assemblies 1 2 CB[8] 2 –5 2 CB[8] 2 were then combined in equimolar amounts, and kept at 40 1C for up to 12 days until equilibrium was reached, as shown by 1 H- and 19 F-NMR spectroscopy.§All mixtures are in slow exchange regime, allowing the observation of new (integrable) signals corresponding to the heteroquaternary assemblies formed at the expense of their homoquaternary counterparts (Fig. 2 and Fig. S122–S151, ESI†).§All heteroquaternary assemblies could also be fully characterized by high resolution mass spectrometry (HRMS).§The observed decline in social behaviour with increasing calculated hinge angle mismatch (or d Pt–Pt difference) was expected. Most mixtures self-sort narcissistically (from +0.90 to +0.02 kcal mol 1 , corresponding to 18% and 49% of the heteroassembly, respectively, Table 1). A notable exception is heteroassembly 23CB[8] 2 , that bears dibenzothiophene and fluorenone spacer groups, which is favoured over the corresponding pair of homoassemblies (free energy term 0.45 (0.05) kcal mol 1 ,i.e. 68% of heterodimer, Table 1). This might suggest some favourable Coulombic interactions in the heteropair of central units; the low hinge angle mismatch (61)andvery short Pt–Pt distance difference (0.22 Å) are not significant enough to cause strain in the heteroassembly. On the other hand, the mixture between assemblies 4 2 CB[8] 2 and 5 2 CB[8] 2 ,whichwas expected to account for distinct spacer group length (d Pt–Pt difference of 2.08 Å) while keeping a nearly negligible hinge angle mismatch (61), returns a small free energy term of +0.16 (0.01) kcal mol 1 , or 43% of heteroquaternary assembly. To rationalize these findings, we questioned whether a set of simple descriptors for each Pt complex 1–5 could be used to predict the self-sorting behaviour of any pair. To that aim, we carried out geometry optimizations with the GFN2-xTB 19–21 semiempirical method§at various constrained Pt–Pt distances.8 Remarkably, they were found to behave as simple harmonic oscillators, as stabilities as a function of Pt–Pt distances can be fit with quadratic eqn (3), where k a isthespringconstantin Fig. 1 Novel complexes (1–5) synthesized in this work along with their ‘‘hinge’’ angle definition.¶ Oligo(ethylene glycol) chains have been grafted to compound 1to increase its solubility. Previously reported homobimetallic Pt-tpy complexes. 16 GFN2-xTB optimized structure of assembly 31 CB[8] 2 .8See Scheme 1 for the structure of group R. Fig. 2 Stacked 1 H-NMR (500 MHz, D 2 O/CD 3 CN 2:1) of (a) assembly 3 2  CB[8] 2 , (b) assembly 1 2 CB[8] 2 , (c) their mixture at t= 0, and (d) their mixture after 12 days. Communication ChemComm Open Access Article. Published on 28 October 2024. Downloaded on 2/17/2025 8:23:37 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2024 Chem. Commun., 2024, 60, 14109–14112 | 14111 kcal mol 1 Å 2 ,a 0 the optimal Pt–Pt distance, and athe Pt–Pt distance in the constrained systems (Fig. 3). The spring constants k0were also expressed in the more commonly used mDyn Å 1 unit using eqn (4), where N A istheAvogadroconstant(Table2). DE¼1 2kaa0a ðÞ 2(3) k0¼4:184 1021 k NA ¼6:948 103k(4) As complexes 1–5 behave like simple harmonic oscillators along a Pt–Pt axis, we then questioned whether self-sorting was simply governed by the deformation of the pair of oscillators when forced to stack on top of one another upon CB[8] encapsulation of the tpy head groups. In that case, a harmonic correlation should be observed between the free energy terms of the self-sorting process and the differences in Pt–Pt distances between the complexes. Remarkably,thisisindeedthecase (Fig. 4, R 2 = 0.997), with the exception of complex 23CB[8] that self-sorts socially (Table 1, reddotinFig.4).Onecanreadily notice that major mismatches in Pt–Pt distances (up to 5.0 Å) or hinge angles (up to 641) still allow the formation of heteroassemblies; free energy terms of the self-sorting processes do not exceed +0.90 kcal mol 1 .Yet,Fig.3showsthatcompressingor expanding the Pt–Pt distance in complexes 1–5 does result in significant destabilization (0.6–2.3 kcal mol 1 for 1 Å distortions). Distortions imposed to Pt complexes to reach identical Pt–Pt distances in hetero-quaternary assemblies can be calculated using eqn (5)–(7) and the following hypothesis: the distortion energy DE hetero is modelled as the sum of the harmonic oscillator of both Pt complexes, where k a and k b are the spring constants (eqn (5)), aand bthe Pt–Pt distances in the distorted hetero-quaternary assemblies, and a 0 and b 0 the distances in the optimized (relaxed) Pt complexes. With both Pt–Pt distances forced to be equal in the hetero-quaternary complexes (i.e. a =b, Fig. 5), eqn (5) and (6) afford the Pt–Pt distances that minimize the energetic penalty DE hetero asshownineqn(7)(Tables3and4,tophalf). DEhetero ¼1 2kaaa0 ðÞ 2þ1 2kbbb0 ðÞ 2(5) dDEhetero da¼0(6) a¼b¼kaa0þkbb0 kaþkb (7) The comparison of experimental self-sorting energy terms (Table 1) with the much higher calculated penalties imposed to match Pt–Pt distances in the hetero-complexes (Table 3) clearly shows that Pt–Pt distances do not need to be equal to allow hetero-quaternary assembly formation. In other terms, CB[8] encapsulation of the tpy heads remains highly favourable Table 1 Free energy term DG[kcal mol 1 ] for the self-sorting of Pt complexes 1–5 in the presence of CB[8]. Experiments repeated in triplicate 12 34 2+0.02 (0.01) 3+0.07 (0.01) 0.45 (0.05) 4+0.90 (0.01) +0.62 (0.01) +0.56 (0.02) 5+0.16 (0.01) Fig. 3 Energy profiles [kcal mol 1 ] of Pt complexes 1–5 as a function of Pt–Pt distance [Å] (electronic component at 0 K relative to the optimized relaxed structure). Table 2 Physical parameters of complexes 1–5 when treated as harmonic oscillators Hinge angle [1]d Pt–Pt [Å] k[kcal mol 1 Å 2 ]k0[mDyn Å 1 ] 162 10.21 1.22 (0.07) 8.5 (0.5) 10 3 277 11.23 1.26 (0.05) 8.7 (0.4) 10 3 383 11.45 1.16 (0.04) 8.0 (0.3) 10 3 4126 15.24 2.3 (0.2) 1.6 (0.1) 10 2 5120 13.16 4.6 (0.5) 3.2 (0.3) 10 2 Fig. 4 Free energy terms for self-sorting processes between pairs of Pt complexes 1–5, as a function of their Pt–Pt distance mismatches. Harmonic correlation (blue line); outlier in red; ‘‘nn’’ refers to all pairs of identical Pt complexes (returning DG=0). Fig. 5 Extension and contraction imposed to a pair of Pt–Pt harmonic oscillators to reach even lengths. ChemComm Communication Open Access Article. Published on 28 October 2024. Downloaded on 2/17/2025 8:23:37 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online 14112 | Chem. Commun., 2024, 60, 14109–14112 This journal is © The Royal Society of Chemistry 2024 despite significant geometrical mismatches between the complexes. We therefore propose the following model to estimate the Pt–Pt distortions Daand Dbin the hetero assemblies: from eqn (8) (equivalent to eqn (5), with DG exp being the experimental self-sorting energy term), the difference between both distortions is minimized to extract both Daand Dbterms (eqn (9)–(11)).§Table 4 (bottom half) shows that required Pt– Pt distance distortions in the hetero assemblies are only 30 (1)% of what they would otherwise have to be to match Pt–Pt distances (Table 4, top half, and Fig. S154 and S155, ESI†). DGexp ¼1 2kaDa2þ1 2kbDb2(8) D=DaDb(9) dD dDa¼0(10) Da¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2kbDGexp kakaþkb ðÞ s;Db¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2kaDGexp kbkaþkb ðÞ s(11) In summary, we showed that two CB[8] units can readily secure hetero-pairs of Pt complexes 1–5, despite major structural differences between the Pt complexes. While Pt arylacetylides are flexible along their C(aryl)–CRC–Pt axis, it is the tolerance of CB[8] for the orientation and position of the encapsulated aryl tpy head groups that drives the assembly. Self-sorting can be quantified using a harmonic oscillator model based on Pt–Pt distance mismatches, predicting most self-sorting patterns. However, the model is limited, as short Pt– Pt distance mismatches may introduce favourable coulombic interactions between compatible spacer groups that promote social behaviour. Ongoing studies are addressing this limitation. We thank the Spanish Ministry of Science and Innovation (PID2021-124691NB-I00, funded by MCIN/AEI/10.13039/ 501100011033/FEDER, UE). EM is grateful to the National Science Foundation (grants CHE-1507321 and CHE-1905238), the American Chemical Society Petroleum Research Fund (grant 56375- ND4), the Roenigk Family Foundation and Ohio University for their continuing financial support. A. S.-M. acknowledges support from Fundacio ´nRamo ´n Areces for a Postdoctoral fellowship. N. A.-L. acknowledges the University of Valladolid and Santander Bank for a predoctoral contract. A. D.-V. acknowledges the Ministerio de Universidades and the European Union – NextGenerationEU for the Marı ´a Zambrano contract (ADV, CONVREC-2021-264). Data availability All synthetic and analytical data is provided in the ESI.†All details about mathematical and computational methods, including associated coordinates are available in the ESI.† Conflicts of interest There are no conflicts to declare. Notes and references §See ESI†for further discussions and details. ¶ According to crystal structures obtained from the Cambridge Structural Database (CSD). 8We must note that the Pt–Pt distances in the hetero-quaternary complexes presented herein should be seen as average distances from a large ensemble of conformers with similar energies but varying geometries; for example, the mere rotation of the CB[8]-bound aryl substituents at position 40of the tpy ligand clockwise or counterclockwise relative to the tpy planes yields at least 4 different hetero-quaternary assemblies with similar stabilities but different Pt–Pt distances. Reported optimized structures only represent one of these arrangements. 1 S. Liu, et al.,J. Am. Chem. Soc., 2005, 127, 15959–15967. 2 J. Lagona, et al.,Angew. Chem., Int. Ed., 2005, 44, 4844–4870. 3 E. Masson, et al.,RSC Adv., 2012, 2, 1213–1247. 4 S. J. Barrow, et al.,Chem. Rev., 2015, 115, 12320–12406. 5 K. I. Assaf and W. M. Nau, Chem. Soc. Rev., 2015, 44, 394–418. 6 J. Murray, et al.,Chem. Soc. Rev., 2017, 46, 2479–2496. 7 X. Yang, et al.,Angew. Chem., Int. Ed., 2020, 59, 21280–21292. 8 G. Wu, et al.,J. Am. Chem. Soc., 2022, 144, 14962–14975. 9 G. Wu, et al.,Chem. Sci., 2020, 11, 812–825. 10 G. Wu, et al.,Angew. Chem., Int. Ed., 2020, 59, 15963–15967. 11 K. Kotturi and E. Masson, Chem. – Eur. J., 2018, 24, 8670–8678. 12 H. Barbero, et al.,J. Am. Chem. Soc., 2020, 142, 867–873. 13 H. Barbero and E. Masson, Chem. Sci., 2021, 12, 9962–9968. 14 R. Rabbani, et al.,Chem. Sci., 2021, 12, 15347–15352. 15 M. Raeisi, et al.,J. Am. Chem. Soc., 2018, 140, 3371–3377. 16 N. A. Thompson, et al.,Chem. Commun., 2019, 3, 12160–12163. 17 P. S. Bols and H. L. Anderson, Acc. Chem. Res., 2018, 51, 2083–2092. 18 K. Miki and K. Ohe, Chem. – Eur. J., 2020, 26, 2529–2575. 19 S. Grimme, et al.,J. Chem. Theory Comput., 2017, 13, 1989–2009. 20 C. Bannwarth, et al.,J. Chem. Theory Comput., 2019, 15, 1652–1671. 21 C. Bannwarth, et al.,Wiley Interdiscip. Rev. Comput. Mol. Sci., 2021, 11, e1493. Table 3 Energy penalty in kcal mol 1 required to match Pt–Pt distances in hetero-quaternary assemblies of Pt complexes 1–5 in the presence of CB[8] 1234 2+0.31 3+0.48 +0.02 4+5.31 +3.57 +2.77 5+1.67 Table 4 Top half: proposed average Pt–Pt distortion [Å] of Pt complexes Aupon quaternary CB[8]-secured complex formation with Pt complex B with identical Pt–Pt distances in the hetero-complexes, based on eqn (5)– (7). Bottom half: Proposed average Pt–Pt distortion [Å] of Pt complexes A upon quaternary CB[8]-secured complex formation with Pt complex B, based on eqn (8)–(11) B12345 A 10 +0.51 +0.60 +3.29 20.49 0 a +2.60 30.63 a 0 +2.53 41.73 1.41 1.26 0 1.38 5+0.70 0 10 +0.13 +0.24 +0.98 20.13 0 a +0.80 30.25 a 0 +0.80 40.52 0.43 0.40 0 0.30 5+0.15 0 a Social self-sorting observed experimentally; model not applicable. Communication ChemComm Open Access Article. Published on 28 October 2024. Downloaded on 2/17/2025 8:23:37 AM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. View Article Online