Solvation rules: aromatic interactions outcompete cation-π interactions in synthetic host-guest complexes in water

Chemical double mutant cycles were used to measure the interaction of a N-methyl pyridinium cation with a π-box in a calix[4]pyrrole receptor. Although the cation-π interaction is attractive (-11 kJ mol-1), it is 7 kJ mol-1 less favourable than the corresponding aromatic interaction with the isosteric but uncharged tolyl group.

Aromatic interactions with heterocycles in water

Conformationally well-defined supramolecular complexes that can be studied in different solvents provide a platform for separating and quantifying free energy contributions due to functional group interactions and desolvation. Here 1:1 complexes formed between four different calix[4]pyrrole receptors and eleven different pyridine N-oxide guests have been used to dissect the factors that govern aromatic interactions with heterocycles in water and in chloroform solution. 1H NMR spectroscopy shows that the three-dimensional structures of the complexes are fixed by four H-bonding interactions between the pyrrole donors at the bottom of the receptor and the N-oxide acceptor on the guest, locking the geometrical arrangement of interacting functional groups in the binding pocket at the other end of the receptor. An aromatic heterocycle on the guest makes two stacking interactions and two edge-to-face interactions with the side walls of the receptor. Chemical double mutant cycles were used to measure the free energy contribution of these four aromatic interactions to the overall stability of the complex. In chloroform, the aromatic interactions measured with pyridine, pyrimidine, furan, thiophene and thiazole are similar to the interactions with a phenyl group, but the effect of introducing a heteroatom depends on where it sits with respect to the aromatic side-walls of the cavity. A nitrogen lone pair directed into a π-face of the side-walls of the binding site leads to repulsive interactions of up to 8 kJ mol-1. In water, the heterocycle aromatic interactions are all significantly more favourable (by up to 12 kJ mol-1). For the non-polar heterocycles, furan and thiophene, the increase in interaction energy correlates directly with hydrophobicity, as measured by the free energy of transfer of the heterocycle from n-hexadecane into water (ΔG°(water-hex)). For the heterocycles with polar nitrogen H-bond acceptors, water can access cracks in the walls of the receptor binding site to solvate the edges of the heterocycles without significantly affecting the geometry of the aromatic interactions, and these nitrogen-water H-bonds stabilise the complexes by about 15 kJ mol-1. The results highlight the complexity of the solvation processes that govern molecular recognition in water.

Substituent effects on aromatic interactions in water

Molecular recognition in water involves contributions due to polar functional group interactions, partial desolvation of polar and non-polar surfaces and changes in conformational flexibility, presenting a challenge for rational design and interpretation of supramolecular behaviour. Conformationally well-defined supramolecular complexes that can be studied in both water and non-polar solvents provide a platform for disentangling these contributions. Here 1 : 1 complexes formed between four different calix[4]pyrrole receptors and thirteen different pyridine N-oxide guests have been used to dissect the factors that govern substituent effects on aromatic interactions in water. H-bonding interactions between the receptor pyrrole donors and the guest N-oxide acceptor at one end of the complex lock the geometrical arrangement of a cluster of aromatic interactions at the other end of the complex, so that a phenyl group on the guest makes two edge-to-face and two stacking interactions with the four aromatic side-walls of the receptor. The thermodynamic contribution of these aromatic interactions to the overall stability of the complex was quantified by chemical double mutant cycles using isothermal titration calorimetry and ¹H NMR competition experiments. Aromatic interactions between the receptor and a phenyl group on the guest stabilise the complex by a factor of 1000, and addition of substituents to the guest phenyl group further stabilises the complex by an additional factor of up to 1000. When a nitro substituent is present on the guest phenyl group, the complex has a sub-picomolar dissociation constant (370 fM). The remarkable substituent effects observed in water for these complexes can be rationalised by comparison with the magnitude of the corresponding substituent effects measured in chloroform. In chloroform, the double mutant cycle free energy measurements of the aromatic interactions correlate well with the substituent Hammett parameters. Electron-withdrawing substituents increase the strength of the interactions by a factor of up to 20, highlighting the role of electrostatics in stabilising both the edge-to-face and stacking interactions. The enhanced substituent effects observed in water are due to entropic contributions associated with the desolvation of hydrophobic surfaces on the substituents. The flexible alkyl chains that line the open end of the binding site assist the desolvation of the non-polar π-surfaces of polar substituents, like nitro, but at the same time allow water to interact with the polar H-bond acceptor sites on the substituent. This flexibility allows polar substituents to maximise non-polar interactions with the receptor and polar interactions with the solvent, leading to remarkably high binding affinities.

Using Solid Catalysts in Disulfide-Based Dynamic Combinatorial Solution- and Mechanochemistry

It was shown for the first time that solid amines can act as catalysts for disulfide-based dynamic combinatorial chemistry (DCC) by ball mill grinding. The mechanochemical equilibrium for the two disulfide reactions studied was reached within 1-3 h using ten different amine catalysts. This contrasts with the weeks to months to achieve solution equilibrium for most solid amine catalysts at 2 %mol mol^-1 concentration in a 2 mMolar disulfide dynamic combinatorial library in a suitable solvent. The final mechanochemical equilibrium was independent of the catalyst used but varied with other ball mill grinding factors such as the presence of traces of solvent. The different efficiencies of the amines tested were discussed.

Quantitative reversible one pot interconversion of three crystalline polymorphs by ball mill grinding

We demonstrate here using a disulfide system the first example of reversible, selective, and quantitative transformation between three crystalline polymorphs by ball mill grinding.

What has Royal Society Open Science achieved in its first few years?

It has been a pleasure and a privilege to serve as the first Editor-in-Chief of Royal Society Open Science for the past 6 years. I step down at the end of December 2021, having completed two 3-year terms, and am taking the opportunity here to reflect on some of the successes and challenges that the journal has experienced and the innovations that we have introduced. When I was first approached back in 2015, the breadth of the journal, covering the whole of science, resonated with my own interests: my research career has ranged across the entire landscape of chemistry, while my leadership roles have embraced all of science, technology and medicine. The open access ethos, the objective refereeing policy that rejects the idea of only publishing what is in fashion, and the opportunities offered by a new venture that could transcend traditional disciplinary boundaries also all appealed to me. Among our successful innovations are Registered Reports, Replication Studies and the new 'Science, Society and Policy' section. The challenges have included the transition to paid article processing charges (APCs), whether to resist pressure to retract a controversial paper, and bullying of young female authors by established senior males in the same field. I explore all of these below, provide some statistics on the journal's performance, also cover some of the notable papers we have published, and provide some concluding thoughts.

Understanding the unexpected effect of frequency on the kinetics of a covalent reaction under ball-milling conditions

We here explore how ball-mill-grinding frequency affects the kinetics of a disulfide exchange reaction. Our kinetic data show that the reaction progress is similar at all the frequencies studied (15–30 Hz), including a significant induction time before the nucleation and growth process starts. This indicates that to start the reaction an initial energy accumulation is necessary. Other than mixing, the energy supplied by the mechanical treatment has two effects: (i) reducing the crystal size and (ii) creating defects in the structure. The crystal-breaking process is likely to be dominant at first becoming less important later in the process when the energy supplied is stored at the molecular level as local crystal defects. This accumulation is taken here to be the rate-determining step. We suggest that the local defects accumulate preferentially at or near the crystal surface. Since the total area increases exponentially when the crystal size is reduced by the crystal-breaking process, this can further explain the exponential dependence of the onset time on the milling frequency.

On the prevalence of smooth polymorphs at the nanoscale: implications for pharmaceuticals

Ball mill neat grinding leads to smoother whereas liquid assisted grinding leads to rougher pharmaceutical forms.

Understanding the Influence of Surface Solvation and Structure on Polymorph Stability: A Combined Mechanochemical and Theoretical Approach

We explore the effect of solvent concentration on the thermodynamic stability of two polymorphs of a 1:1 cocrystal of theophylline and benzamide subjected to ball-mill liquid assisted grinding (LAG) and we investigate how this can be related to surface solvent solvation phenomena. In this system, most stable bulk polymorph form II converts to metastable bulk polymorph form I upon neat grinding (NG), while form I can fully or partially transform into form II under LAG conditions, depending on the amount of solvent used. Careful and strict experimental procedures were designed to achieve polymorph equilibrium under ball-mill LAG conditions for 16 different solvents. This allowed us to determine 16 equilibrium polymorph concentration curves as a function of solvent concentration. Ex-situ powder X-ray diffraction (PXRD) was used to monitor the polymorph concentration and crystallite size. The surface site interactions point (SSIP) description of noncovalent interactions was used in conjunction with the SSIMPLE method for calculating solvation energies to determine which functional groups are more or less exposed on the polymorph crystal surfaces. Our results demonstrate that (i) ball-mill LAG equilibrium curves can be successfully achieved experimentally for a cocrystal system; (ii) the equilibrium curves vary from solvent to solvent in onset values and slopes, thus confirming the generality of the interconversion phenomenon that we interpret here in terms of cooperativity; (iii) the concentration required for a switch in polymorphic outcome is dependent on the nature of the solvent; (iv) the SSIP results indicate that the theophylline π-system face is more exposed on the surface of form I while the theophylline N-methyl groups are more exposed in form II; and (v) for some solvents, form II has a significantly smaller crystal size at equilibrium than form I in the investigated solvent concentration range. Therefore, the free energy of the 1:1 cocrystal of theophylline and benzamide polymorphs studied here must be affected by surface solvation under ball-mill LAG conditions.

Editorial for 'New talent' collection: introducing New Talent in chemistry and material science

We introduce 14 articles published as part of the 'New talent' special collection of invited articles to showcase some of the exciting work being funded by the Royal Society. As Royal Society University Research Fellows or Dorothy Hodgkin Fellowship holders, the contributors to this collection are rising stars in their areas of research. This collection also illustrates the close collaboration between Royal Society Open Science and the Royal Society of Chemistry. The collection spans the range of the chemical and material sciences, reflecting the breadth of research areas receiving Royal Society grant support.

Reliable Mechanochemistry: Protocols for Reproducible Outcomes of Neat and Liquid Assisted Ball-mill Grinding Experiments

The equilibrium outcomes of ball mill grinding can dramatically change as a function of even tiny variations in the experimental conditions such as the presence of very small amounts of added solvent. To reproducibly and accurately capture this sensitivity, the experimentalist needs to carefully consider every single factor that can affect the ball mill grinding reaction under investigation, from ensuring the grinding jars are clean and dry before use, to accurately adding the stoichiometry of the starting materials, to validating that the delivery of solvent volume is accurate, to ensuring that the interaction between the solvent and the powder is well understood and, if necessary, a specific soaking time is added to the procedure. Preliminary kinetic studies are essential to determine the necessary milling time to achieve equilibrium. Only then can exquisite phase composition curves be obtained as a function of the solvent concentration under ball mill liquid assisted grinding (LAG). By using strict and careful procedures analogous to the ones here presented, such milling equilibrium curves can be obtained for virtually all milling systems. The system we use to demonstrate these procedures is a disulfide exchange reaction starting from the equimolar mixture of two homodimers to obtain at equilibrium quantitative heterodimer. The latter is formed by ball mill grinding as two different polymorphs, Form A and Form B. The ratio R = [Form B] / ([Form A] + [Form B]) at milling equilibrium depends on the nature and concentration of the solvent in the milling jar.

Selective C70 encapsulation by a robust octameric nanospheroid held together by 48 cooperative hydrogen bonds

Self-assembly of multiple building blocks via hydrogen bonds into well-defined nanoconstructs with selective binding function remains one of the foremost challenges in supramolecular chemistry. Here, we report the discovery of a enantiopure nanocapsule that is formed through the self-assembly of eight amino acid functionalised molecules in nonpolar solvents through 48 hydrogen bonds. The nanocapsule is remarkably robust, being stable at low and high temperatures, and in the presence of base, presumably due to the co-operative geometry of the hydrogen bonding motif. Thanks to small pore sizes, large internal cavity and sufficient dynamicity, the nanocapsule is able to recognize and encapsulate large aromatic guests such as fullerenes C₆₀ and C₇₀. The structural and electronic complementary between the host and C₇₀ leads to its preferential and selective binding from a mixture of C₆₀ and C₇₀.

Individual hydrogen bonds are weak, so self-assembling multiple components via hydrogen bonding is a significant challenge. Here the authors report a robust, enantiopure nanocapsule held together by 48 cooperative hydrogen bonds, and use it for the selective binding of C₇₀.

Solvation and surface effects on polymorph stabilities at the nanoscale

We explore the effects of particle size and solvent environment on the thermodynamic stability of two pairs of polymorphs subjected to ball-mill neat grinding (NG) and liquid assisted grinding (LAG).

We explore the effects of particle size and solvent environment on the thermodynamic stability of two pairs of polymorphs subjected to ball-mill neat grinding (NG) and liquid assisted grinding (LAG). Two systems were studied: (i) forms I and II of a 1 : 1 theophylline : benzamide cocrystal and (ii) forms A and B of an aromatic disulfide compound. For both systems, the most stable-bulk polymorph converted to the metastable-bulk polymorph upon NG. LAG experiments yielded different outcomes depending on the amount of solvent used. This was further investigated by performing carefully controlled LAG experiments with increasing μL amounts of solvents of different nature. With these experiments, we were able to monitor form A to B and form I to II conversions as a function of solvent concentration and derive polymorph equilibrium curves. The concentration required for a switch in polymorphic outcome was found to be dependent on solvent nature. We propose that these experiments demonstrate a switch in thermodynamic stability of the polymorphs in the milling jar. Form B, the stable-bulk polymorph, has less stable surfaces than form A, thus becoming metastable at the nanoscale when surface effects become important. Ex situ diffraction and electron microscopy data confirm crystal sizes in the order of tens of nanometers after the ball mill grinding experiments reach equilibrium. DFT-d computations of the polymorph particles stabilities support these findings and were used to calculate cross-over sizes of forms A and B as a function of solvent. Attachment energies and surface stabilities of the various crystalline faces exposed were found to be very sensitive to the solvent environment. Our findings suggest that surface effects are significant in polymorphism at the nanoscale and that the outcomes of equilibrium ball-mill NG and LAG experiments are in general controlled by thermodynamics.

Catenation and encapsulation induce distinct reconstitutions within a dynamic library of mixed-ligand Zn4L6 cages

Two new Zn4L6 cages composed of diamine subcomponents containing either naphthalene diimide (NDI) or porphyrin moieties are described. Their structural differences allow these cages to exhibit distinct interactions with different chemical stimuli, yielding different supramolecular products. The electron-poor NDI subunits of the first cage were observed to thread through electron-rich aromatic crown-ether macrocycles, forming mechanically-interlocked species up to a [3]catenane, whereas the porphyrin ligands of the second cage interacted favourably with C70, causing it to be bound as a guest. When mixed, the two cages were observed to form a dynamic combinatorial library (DCL) of seven constitutionally distinct mixed-ligand Zn4L6 cages. The DCL was observed to reconstitute in opposing ways when treated with either the crown ether or C70: the electron-rich macrocycle templated the formation of heteroleptic catenanes, whereas C70 caused the DCL to self-sort into homoleptic structures.

Guest-induced transformation of a porphyrin-edged Fe(II)4L6 capsule into a Cu(I)Fe(II)2L4 fullerene receptor

The combination of a bent diamino(nickel(II) porphyrin) with 2-formylpyridine and Fe(II) yielded an Fe(II) 4 L6 cage. Upon treatment with the fullerenes C60 or C70 , this cage was found to transform into a new host-guest complex incorporating three Fe(II) centers and four porphyrin ligands, in an arrangement that is hypothesized to maximize π interactions between the porphyrin units of the host and the fullerene guest bound within its central cavity. The new complex shows coordinative unsaturation at one of the Fe(II) centers as the result of the incommensurate metal-to-ligand ratio, which enabled the preparation of a heterometallic cone-shaped Cu(I) Fe(II) 2 L4 adduct of C60 or C70 .

Molecular motion of donor–acceptor catenanes in water

In this article, we use ¹H NMR spectroscopy to study the spontaneous molecular motion of donor–acceptor [2]catenanes in water.

Direct observation of intermediates in a thermodynamically controlled solid-state dynamic covalent reaction

We present the first polymorph interconversion study that uses solid-state dynamic covalent chemistry (DCC). This system exhibits unexpected and rich behavior, including the observation that under appropriate conditions the polymorph interconversion of a heterodimer proceeds through reversible covalent chemistry intermediates, and this route is facilitated by one of the two disulfide homodimers involved in the reaction. Furthermore, we demonstrate experimentally that in all cases a dynamic equilibrium is reached, meaning that changing the milling conditions affects the free energy difference between the two polymorphs and thus their relative stability. We suggest that this effect is due to the surface solvation energy combined with the high surface to volume ratio of the nanocrystalline powder.

Disulfide exchange: exposing supramolecular reactivity through dynamic covalent chemistry

A Tutorial Review of the subtle supramolecular interactions influencing the outcomes of equilibrating systems, focusing on the dynamic covalent chemistry (DCC) of disulfide exchange reactions, is presented. We discuss the topics of cation-π interactions (2.1), hydrophobic effects (2.2), hydrogen bonding interactions (2.3) aromatic donor-acceptor interactions (2.4), and metal-ligand interactions (2.5) in the context of dynamic disulfide chemistry.

Generation of a dynamic system of three-dimensional tetrahedral polycatenanes

Seven of the best: A dynamic combinatorial library of polycatenated tetrahedra was prepared by complexation between a dynamic Fe₄L₆ tetrahedral cage, constructed from ligands containing an electron‐deficient naphthalenediimide core, and an electron‐rich aromatic crown ether, 1,5‐dinaphtho[38]crown‐10. The highest order species in the library is the tetrahedral [7]catenane.WILEY-VCH

Diverse topologies in dynamic combinatorial libraries from tri- and mono-thiols in water: sensitivity to weak supramolecular interactions

We describe here the use of disulfide exchange reactions to generate distinct supramolecular topologies within a dynamic combinatorial library (DCL) derived from tri- and mono-functionalised thiol components in water. This approach allows for the effective expansion of structural diversity within a DCL, and also provides a sensitive probe for weak supramolecular interactions.

Evolution of dynamic combinatorial chemistry

Since its inception in the mid-1990s, dynamic combinatorial chemistry (DCC), the chemistry of complex systems under thermodynamic control, has proved valuable in identifying unexpected molecules with remarkable binding properties and in providing effective synthetic routes to complex species. Essentially, in this approach, one designs the experiment rather than the molecule. DCC has also provided us with insights into how some chemical systems respond to external stimuli. Using examples from the work of our laboratory and others, this Account shows how the concept of DCC, inspired by the evolution of living systems, has found an increasing range of applications in diverse areas and has evolved conceptually and experimentally. A dynamic combinatorial library (DCL) is a thermodynamically controlled mixture of interconverting species that can respond to various stimuli. The Cambridge version of dynamic combinatorial chemistry was initially inspired by the mammalian immune system and was conceived as a way to create and identify new unpredictable receptors. For example, an added template can select and stabilize a strongly binding member of the library which is then amplified at the expense of the unsuccessful library members, minimizing the free energy of the system. But researchers have exploited DCC in a variety of other ways: over the past two decades, this technique has contributed to the evolution of chemistry and to applications in the diverse fields of catalysis, fragrance release, and responsive materials. Among these applications, researchers have built intricate and well-defined architectures such as catenanes or hydrogen-bonded nanotubes, using the ability of complex chemical systems to reach a high level of organization. In addition, DCC has proved a powerful tool for the study of complex molecular networks and systems. The use of DCC is improving our understanding of chemical and biological systems. The study of folding or self-replicating macrocycles in DCLs has served as a model for appreciating how complex organisations such as life can emerge from a pool of simple chemicals. Today, DCC is no longer restricted to thermodynamic control, and new systems have recently appeared in which kinetic and thermodynamic control coexist. Expanding the realm of DCC to unexplored and promising new territories, these hybrid systems show that the concept of dynamic combinatorial chemistry continues to evolve.

Metal induced folding: synthesis and conformational analysis of the lanthanide complexes of two 44-membered hydrazone macrocycles

Six new lanthanide complexes of two 44-membered macrocycles have been prepared and characterised in solution. An analysis of the conformations of the free macrocycles and their lanthanide complexes both in solution (2D NMR) and in solid state (X-ray crystallography) demonstrate that the complexation induces changes in folding of the macrocycles.

Dynamic combinatorial libraries for the recognition of heavy metal ions

We present the use of hydrazone dynamic combinatorial libraries (DCLs) to identify macrocyclic receptors that are selective for alkaline earth metal ions over alkali metal ions. In particular, the toxic heavy metal ions Sr(2+) and Ba(2+) induce characteristic changes in the DCLs. Four macrocycles were isolated and characterised by LCMS, HRMS, NMR and X-ray crystallography; binding studies by UV-Vis spectroscopy confirm the selectivity observed in the DCLs.

From static to dynamic: escaping kinetic traps in hydrazone-based dynamic combinatorial libraries

Thermodynamic control over kinetically-trapped mixtures of hydrazone-based macrocycles is achieved by addition of an aromatic monohydrazide to generate dynamic combinatorial libraries (DCLs) of linear and macrocyclic oligomers.

Discovery of linear receptors for multiple dihydrogen phosphate ions using dynamic combinatorial chemistry

We describe the use of dynamic combinatorial chemistry to discover a new series of linear hydrazone-based receptors that bind multiple dihydrogen phosphate ions. Through the use of a template-driven, selection-based approach to receptor synthesis, dynamic combinatorial chemistry allows for the identification of unexpected host structures and binding motifs. Notably, we observed the unprecedented selection of these linear receptors in preference to competing macrocyclic hosts. Furthermore, linear receptors containing up to nine building blocks and three different building blocks were amplified in the dynamic combinatorial library. The receptors were formed using a dihydrazide building block based on an amino acid-disubstituted ferrocene scaffold. A detailed study of the linear pentamer revealed that it forms a helical ditopic receptor that employs four acylhydrazone hydrogen-bond donor motifs to cooperatively bind two dihydrogen phosphate ions.