High-Throughput Computational Evaluation of Low Symmetry Pd2 L4 Cages to Aid in System Design*

Photoswitchable Catalysis by a Self-Assembled Molecular Cage

A heteroleptic [Pd2L2L'2]4+ coordination cage containing a photoswitchable azobenzene-derived ligand catalyzes the Michael addition reaction between methyl vinyl ketone and benzoyl nitromethane within its cavity. The corresponding homoleptic cages are catalytically inactive. The heteroleptic cage can be reversibly disassembled and reassembled using 530 and 405 nm light, respectively, allowing catalysis within the cage to be switched OFF and ON at will.

Non-Templated Assembly of D5h-Symmetric Pd5L10 Rings by Precise Ligand Angle Adjustment

We report a series of Pd(II)nL2n coordination rings for which nuclearity is controlled by the binding angle of the corresponding bis-monodentate bridging ligands. Judicious choice of the angle within a family of rather rigid ligands allowed for the first-time to synthesize a homoleptic five-membered Pd5L10 ring that does not require any template to form. We demonstrate that control over the ring size is maintained both in the solid-, solution-, and gas-phase. Two X-ray structures of five-membered rings from ligands with ideal angles (yielding a perfect pentagonal ring) vs. suboptimal angles (resulting in a highly distorted structure) illustrate the importance of the correct ligand geometry. A mathematical model for estimating the expected ring size based on the ligand angle was derived and DFT computations show that ring-strain is the major factor determining the assembly outcome.

Orientational Self-Sorting in Octahedral Palladium Cages: Scope and Limitations of the "cis Rule"

Octahedral coordination cages of the general formula [Pd6L12](BF4)12 were obtained by combining [Pd(CH3CN)4](BF4)2 with heteroditopic N-donor ligands. Four different ligands were employed. These ligands have 3-pyridyl donor groups at one end and 4-pyridyl, imidazolyl, or triazolyl donor groups at the other end. According to a geometric analysis, cages with a cis configuration at the six metal centers should be preferred ("cis rule"). This prediction was corroborated by spectroscopic data and crystallographic analyses. Limitations of the "cis rule" were also encountered, and possible explanations are discussed.

Self-assembly of conformation-adaptive dihydrophenazine-based coordination cages

In this study, we designed and synthesized three conformation-adaptive Pd2L4- and Pd3L6-type coordination cages based on three dihydrophenazine-based ligands with different lengths. Interestingly, the shorter ligands L1 and L2 self-assembled into Pd2L4-type coordination cages while the longer ligand L3 formed Pd3L6-type one, mainly driven by the anion template effect. All coordination cages were confirmed through single-crystal X-ray diffraction, and their structural conformations underwent great changes compared with those of their corresponding ligands. Moreover, the conformational changes also significantly affected their photophysical and electrochemical properties which were distinct from their parent ligands.

Heterometallic cages: synthesis and applications

High symmetry metallosupramolecular architectures (MSAs) have been exploited for a range of applications including molecular recognition, catalysis and drug delivery. Recently there have been increasing efforts to enhance those applications by generating reduced symmetry MSAs. While there are several emerging methods for generating lower symmetry MSAs, this tutorial review examines the general methods used for synthesizing heterometallic MSAs with a particular focus on heterometallic cages. Additionally, the intrinsic properties of the cages and their potential emerging applications as host-guest systems and reaction catalysts are described.

Diastereoselective Self-Assembly of Low-Symmetry Pdn L2n Nanocages through Coordination-Sphere Engineering

Metal-organic cages (MOCs) are popular host architectures assembled from ligands and metal ions/nodes. Assembling structurally complex, low-symmetry MOCs with anisotropic cavities can be limited by the formation of statistical isomer libraries. We set out to investigate the use of primary coordination-sphere engineering (CSE) to bias isomer selectivity within homo- and heteroleptic Pdn L2n cages. Unexpected differences in selectivities between alternative donor groups led us to recognise the significant impact of the second coordination sphere on isomer stabilities. From this, molecular-level insight into the origins of selectivity between cis and trans diastereoisomers was gained, highlighting the importance of both host-guest and host-solvent interactions, in addition to ligand design. This detailed understanding allows precision engineering of low-symmetry MOC assemblies without wholesale redesign of the ligand framework, and fundamentally provides a theoretical scaffold for the development of stimuli-responsive, shape-shifting MOCs.

Systematic exploration of accessible topologies of cage molecules via minimalistic models

Cages are macrocyclic structures with an intrinsic internal cavity that support applications in separations, sensing and catalysis. These materials can be synthesised via self-assembly of organic or metal-organic building blocks. Their bottom-up synthesis and the diversity in building block chemistry allows for fine-tuning of their shape and properties towards a target property. However, it is not straightforward to predict the outcome of self-assembly, and, thus, the structures that are practically accessible during synthesis. Indeed, such a prediction becomes more difficult as problems related to the flexibility of the building blocks or increased combinatorics lead to a higher level of complexity and increased computational costs. Molecular models, and their coarse-graining into simplified representations, may be very useful to this end. Here, we develop a minimalistic toy model of cage-like molecules to explore the stable space of different cage topologies based on a few fundamental geometric building block parameters. Our results capture, despite the simplifications of the model, known geometrical design rules in synthetic cage molecules and uncover the role of building block coordination number and flexibility on the stability of cage topologies. This leads to a large-scale and systematic exploration of design principles, generating data that we expect could be analysed through expandable approaches towards the rational design of self-assembled porous architectures.

Exploiting reduced-symmetry ligands with pyridyl and imidazole donors to construct a second-generation stimuli-responsive heterobimetallic [PdPtL4]4+ cage

A new sequential metalation strategy that enables the assembly of a new more robust reduced symmetry heterobimetallic [PdPtL4]4+ cage C is reported. By exploiting a low-symmetry ditopic ligand (L) that features imidazole and pyridine donor units we were able to selectively form a [Pt(L)4]2+ "open-cage" complex. When this was treated with Pd(ii) ions the cage C assembled. 1H and DOSY nuclear magnetic resonance (NMR) spectroscopy and electrospray ionisation mass spectrometry (ESIMS) data were consistent with the quantitative formation of the cage and the heterobimetallic structure was confirmed by single crystal X-ray crystallography. The cage C was shown to bind anionic guest molecules. NMR studies suggested that these guests interacted with the cavity of the cage in a specific orientation and this was confirmed for the mesylate ion (MsO-) : C host-guest adduct using X-ray crystallography. In addition, the system was shown to be stimulus-responsive and could be opened and closed on demand when treated with appropriate stimuli. If a guest molecule was bound within the cage, the opening and closing was accompanied by the release and re-uptake of the guest molecule.

Maximized axial helicity in a Pd2L4 cage: inverse guest size-dependent compression and mesocate isomerism

Helicity is an archetypal structural motif of many biological systems and provides a basis for molecular recognition in DNA. Whilst artificial supramolecular hosts are often helical, the relationship between helicity and guest encapsulation is not well understood. We report a detailed study on a significantly coiled-up Pd₂L₄ metallohelicate with an unusually wide azimuthal angle (∼176°). Through a combination of NMR spectroscopy, single-crystal X-ray diffraction, trapped ion mobility mass spectrometry and isothermal titration calorimetry we show that the coiled-up cage exhibits extremely tight anion binding (K of up to 10⁶ M^-1) by virtue of a pronounced oblate/prolate cavity expansion, whereby the Pd-Pd separation decreases for mono-anionic guests of increasing size. Electronic structure calculations point toward strong dispersion forces contributing to these host-guest interactions. In the absence of a suitable guest, the helical cage exists in equilibrium with a well-defined mesocate isomer that possesses a distinct cavity environment afforded by a doubled Pd-Pd separation distance.

Exploiting Supramolecular Interactions to Control Isomer Distributions in Reduced-Symmetry [Pd2L4]4+ Cages

High-symmetry metallosupramolecular architectures (MSAs) have been exploited for a range of applications including molecular recognition, catalysis, and drug delivery. Recently, there have been increasing efforts to enhance those applications by generating reduced-symmetry MSAs. Here we report our attempts to use supramolecular (dispersion and hydrogen-bonding) forces and solvophobic effects to generate isomerically pure [Pd₂(L)₄]⁴⁺ cage architectures from a family of new reduced-symmetry ditopic tripyridyl ligands. The reduced-symmetry tripyridyl ligands featured either solvophilic polyether chains, solvophobic alkyl chains, or amino substituents. We show using NMR spectroscopy, high-performance liquid chromatography, X-ray diffraction data, and density functional theory calculations that the combination of dispersion forces and solvophobic effects does not provide any control of the [Pd₂(L)₄]⁴⁺ isomer distribution with mixtures of all four cage isomers (HHHH, HHHT, cis-HHTT, or trans-HTHT, where H = head and T = tail) obtained in each case. More control was obtained by exploiting hydrogen-bonding interactions between amino units. While the cage assembly with a 3-amino-substituted tripyridyl ligand leads to a mixture of all four possible isomers, the related 2-amino-substituted tripyridyl ligand generated a cis-HHTT cage architecture. Formation of the cis-HHTT [Pd₂(L)₄]⁴⁺ cage was confirmed using NMR studies and X-ray crystallography.

Computational modeling to assist in the discovery of supramolecular materials

Computational modeling is increasingly used to assist in the discovery of supramolecular materials. Supramolecular materials are typically primarily built from organic components that are self-assembled through noncovalent bonding and have potential applications, including in selective binding, sorption, molecular separations, catalysis, optoelectronics, sensing, and as molecular machines. In this review, the key areas where computational prediction can assist in the discovery of supramolecular materials, including in structure prediction, property prediction, and the prediction of how to synthesize a hypothetical material are discussed, before exploring the potential impact of artificial intelligence techniques on the field. Throughout, the importance of close integration with experimental materials discovery programs will be highlighted. A series of case studies from the author's work across some different supramolecular material classes will be discussed, before finishing with a discussion of the outlook for the field.

Pseudo-heterolepticity in Low-Symmetry Metal-Organic Cages

Heteroleptic metal-organic cages, formed through integrative self-assembly of ligand mixtures, are highly attractive as reduced symmetry supramolecular hosts. Ensuring high-fidelity, non-statistical self-assembly, however, presents a significant challenge in molecular engineering due to the inherent difficulty in predicting thermodynamic energy landscapes. In this work, two conceptual strategies are described that circumvent this issue, using ligand design strategies to access structurally sophisticated metal-organic hosts. Using these approaches, it was possible to realise cavity environments described by two inequivalent, unsymmetrical ligand frameworks, representing a significant step forward in the construction of highly anisotropic confined spaces.

Orientational self-sorting in cuboctahedral Pd cages

Cuboctahedral coordination cages of the general formula [Pd12L24]24+ (L = low-symmetry ligand) were analyzed theoretically and experimentally. With 350 696 potential isomers, the structural space of these assemblies is vast. Orientational self-sorting refers to the preferential formation of particular isomers within the pool of potential structures. Geometric and computational analyses predict the preferred formation of cages with a cis arrangement at the metal centers. This prediction was corroborated experimentally by synthesizing a [Pd12L24]24+ cage with a bridging 3-(4-(pyridin-4-yl)phenyl)pyridine ligand. A crystallographic analysis of this assembly showed exclusive cis coordination of the 3- and the 4-pyridyl donor groups at the Pd2+ ions.

Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities

Self-assembled metallo-organic cages have emerged as promising biomimetic platforms that can encapsulate whole substrates akin to an enzyme active site. Extensive experimental work has enabled access to a variety of structures, with a few notable examples showing catalytic behavior. However, computational investigations of metallo-organic cages are scarce, not least due to the challenges associated with their modeling and the lack of accurate and efficient protocols to evaluate these systems. In this review, we discuss key molecular principles governing the design of functional metallo-organic cages, from the assembly of building blocks through binding and catalysis. For each of these processes, computational protocols will be reviewed, considering their inherent strengths and weaknesses. We will demonstrate that while each approach may have its own specific pitfalls, they can be a powerful tool for rationalizing experimental observables and to guide synthetic efforts. To illustrate this point, we present several examples where modeling has helped to elucidate fundamental principles behind molecular recognition and reactivity. We highlight the importance of combining computational and experimental efforts to speed up supramolecular catalyst design while reducing time and resources.

Beyond Platonic: How to Build Metal-Organic Polyhedra Capable of Binding Low-Symmetry, Information-Rich Molecular Cargoes

The field of metallosupramolecular chemistry has advanced rapidly in recent years. Much work in this area has focused on the formation of hollow self-assembled metal-organic architectures and exploration of the applications of their confined nanospaces. These discrete, soluble structures incorporate metal ions as ‘glue’ to link organic ligands together into polyhedra.Most of the architectures employed thus far have been highly symmetrical, as these have been the easiest to prepare. Such high-symmetry structures contain pseudospherical cavities, and so typically bind roughly spherical guests. Biomolecules and high-value synthetic compounds are rarely isotropic, highly-symmetrical species. To bind, sense, separate, and transform such substrates, new, lower-symmetry, metal-organic cages are needed. Herein we summarize recent approaches, which taken together form the first draft of a handbook for the design of higher-complexity, lower-symmetry, self-assembled metal-organic architectures.

Metallosupramolecular cages: from design principles and characterisation techniques to applications

Although metallosupramolecular cages are self-assembled from seemingly simple building blocks, metal ions and organic ligands, architectures of increasingly large size and complexity are accessible and exploited in applications from catalysis to the stabilisation of reactive species. This Tutorial Review gives an introduction to the principles for designing metallosupramolecular cages and highlights advances in the design of large and lower symmetry cages. The characterisation and identification of cages relies on a number of complementary techniques with NMR spectroscopy, mass spectrometry, X-ray crystallography and computational methods being the focus of this review. Finally, examples of cages are discussed where these design principles and characterisation techniques are put into practice for an application or function of the cage.

Unlocking the computational design of metal-organic cages

Metal-organic cages are macrocyclic structures that can possess an intrinsic void that can hold molecules for encapsulation, adsorption, sensing, and catalysis applications. As metal-organic cages may be comprised from nearly any combination of organic and metal-containing components, cages can form with diverse shapes and sizes, allowing for tuning toward targeted properties. Therefore, their near-infinite design space is almost impossible to explore through experimentation alone and computational design can play a crucial role in exploring new systems. Although high-throughput computational design and screening workflows have long been known as powerful tools in drug and materials discovery, their application in exploring metal-organic cages is more recent. We show examples of structure prediction and host-guest/catalytic property evaluation of metal-organic cages. These examples are facilitated by advances in methods that handle metal-containing systems with improved accuracy and are the beginning of the development of automated cage design workflows. We finally outline a scope for how high-throughput computational methods can assist and drive experimental decisions as the field pushes toward functional and complex metal-organic cages. In particular, we highlight the importance of considering realistic, flexible systems.

Heteroleptic tripalladium(II) cages

There is a concerted attempt to develop self-assembled metallo-cages of greater structural complexity, and heteroleptic Pd II cages are emerging as prime candidates in these efforts. Most of these are dinuclear: few examples of higher nuclearity have been reported. We demonstrate here a robust method for the formation of tripalladium(II) cages from the 2:3:3 combination of a tritopic ligand, Pd II , and a selection of ditopic ligands of the correct size and geometry.

PoreMatMod.jl: Julia Package for in Silico Postsynthetic Modification of Crystal Structure Models

PoreMatMod.jl is a free, open-source, user-friendly, and documented Julia package for modifying crystal structure models of porous materials such as metal-organic frameworks (MOFs). PoreMatMod.jl functions as a find-and-replace algorithm on crystal structures by leveraging (i) Ullmann's algorithm to search for subgraphs of the crystal structure graph that are isomorphic to the graph of a query fragment and (ii) the orthogonal Procrustes algorithm to align a replacement fragment with a targeted substructure of the crystal structure for installation. The prominent application of PoreMatMod.jl is to generate libraries of hypothetical structures for virtual screenings. For example, one can install functional groups on the linkers of a parent MOF, mimicking postsynthetic modification. Other applications of PoreMatMod.jl to modify crystal structure models include introducing defects with precision and correcting artifacts of X-ray structure determination (adding missing hydrogen atoms, resolving disorder, and removing guest molecules). The find-and-replace operations implemented by PoreMatMod.jl can be applied broadly to diverse atomistic systems for various in silico structural modification tasks.

Molecular engineering of confined space in metal–organic cages

The host–guest chemistry of metal–organic cages can be modified through tailoring of structural aspects such as size, shape and functionality. In this review, strategies, opportunities and challenges of such molecular engineering are discussed.

Entropy directs the self-assembly of supramolecular palladium coordination macrocycles and cages

The self-assembly of palladium-based cages is frequently rationalized via the cumulative enthalpy (ΔH) of bonds between coordination nodes (M, i.e., Pd) and ligand (L) components. This focus on enthalpic rationale limits the complete understanding of the Gibbs free energy (ΔG) for self-assembly, as entropic (ΔS) contributions are overlooked. Here, we present a study of the M₂^linL₃ intermediate species (M = dinitrato(N,N,N′,N′-tetramethylethylenediamine)palladium(ii), ^linL = 4,4′-bipyridine), formed during the synthesis of triangle-shaped (M₃^linL₃) and square-shaped (M₄^linL₄) coordination macrocycles. Thermochemical analyses by variable temperature (VT) ¹H-NMR revealed that the M₂^linL₃ intermediate exhibited an unfavorable (relative) ΔS compared to M₃^linL₃ (triangle, ΔTΔS = +5.22 kcal mol⁻¹) or M₄^linL₄ (square, ΔTΔS = +2.37 kcal mol⁻¹) macrocycles. Further analysis of these constructs with molecular dynamics (MD) identified that the self-assembly process is driven by ΔG losses facilitated by increases in solvation entropy (ΔS_solv, i.e., depletion of solvent accessible surface area) that drives the self-assembly from “open” intermediates toward “closed” macrocyclic products. Expansion of our computational approach to the analysis of self-assembly in Pd_n^benL_2n cages (^benL = 4,4'-(5-ethoxy-1,3-phenylene)dipyridine), demonstrated that ΔS_solv contributions drive the self-assembly of both thermodynamic cage products (i.e., Pd₁₂^benL₂₄) and kinetically-trapped intermediates (i.e., Pd₈^cL₁₆).

These studies demonstrate that ΔS drives the self-assembly of supramolecular palladium-based coordination macrocycles and cages. As this ΔS contribution arises from solvation, these findings broadly reflect the thermodynamic drive of self-assembly to form compact structures.

Conformational Control of a Metallo-Supramolecular Cage via the Dissymmetrical Modulation of Ligands

In nature as well as life systems, the presence of asymmetrical and dissymmetrical structures with specific functions is extremely common. However, the construction of metallo-supramolecular assemblies based on dissymmetrical ligands still remains a considerable challenge for avoiding the generation of unexpected isomers with similar thermodynamic stabilities, especially for three-dimensional supramolecular structures. In this study, a strategy for the conformational control of metallo-supramolecular cages via the enhancement of ligand dissymmetry was proposed. Four dissymmetrical ditopic ligands were designed and synthesized. By increasing the dissymmetry of length or angle, conformations of assemblies were precisely controlled to form discrete cis-Pd_n L_2n molecular cages.

Orientational self-sorting: formation of structurally defined Pd4L8 and Pd6L12 cages from low-symmetry dipyridyl ligands

Tetra- and hexanuclear coordination cages were obtained in reactions of [Pd(CH₃CN)₄](BF₄)₂ with low-symmetry dipyridyl ligands. In both cases, only one structurally defined complex was formed out of a vast pool of potential isomers.