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

CageCavityCalc (C3): A Computational Tool for Calculating and Visualizing Cavities in Molecular Cages

Organic(porous) and metal-organic cages are promising biomimetic platforms with diverse applications spanning recognition, sensing, and catalysis. The key to the emergence of these functions is the presence of well-defined inner cavities capable of binding a wide range of guest molecules and modulating their properties. However, despite the myriad cage architectures currently available, the rational design of structurally diverse and functional cages with specific host-guest properties remains challenging. Efficiently predicting such properties is critical for accelerating the discovery of novel functional cages. Herein, we introduce CageCavityCalc (C3), a Python-based tool for calculating the cavity size of molecular cages. The code is available on GitHub at https://github.com/VicenteMartiCentelles/CageCavityCalc. C3 utilizes a novel algorithm that enables the rapid calculation of cavity sizes for a wide range of molecular structures and porous systems. Moreover, C3 facilitates easy visualization of the computed cavity size alongside hydrophobic and electrostatic potentials, providing insights into host-guest interactions within the cage. Furthermore, the calculated cavity can be visualized using widely available visualization software, such as PyMol, VMD, or ChimeraX. To enhance user accessibility, a PyMol plugin has been created, allowing nonspecialists to use this tool without requiring computer programming expertise. We anticipate that the deployment of this computational tool will significantly streamline cage cavity calculations, thereby accelerating the discovery of functional cages.

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.

Water-Soluble Molecular Cages for Biological Applications

The field of molecular cages has attracted increasing interest in relation to the development of biological applications, as evidenced by the remarkable examples published in recent years. Two key factors have contributed to this achievement: First, the remarkable and adjustable host-guest chemical properties of molecular cages make them highly suitable for biological applications. This allows encapsulating therapeutic molecules to improve their properties. Second, significant advances have been made in synthetic methods to create water-soluble molecular cages. Achieving the necessary water solubility is a significant challenge, which in most cases requires specific chemical groups to overcome the inherent hydrophobic nature of the molecular cages which feature the organic components of the cage. This can be achieved by either incorporating water-solubilizing groups with negative/positive charges, polyethylene glycol chains, etc.; or by introducing charges directly into the cage structure itself. These synthetic strategies allow preparing water-soluble molecular cages for diverse biological applications, including cages' anticancer activity, anticancer drug delivery, photodynamic therapy, and molecular recognition of biological molecules. In the review we describe selected examples that show the main concepts to achieve water solubility in molecular cages and some selected recent biological applications.

Modeling Amine Methylation in Methyl Ester Cavitand

Methylation of amines inside an introverted resorcinarene-based deep methyl ester cavitand is investigated by means of molecular dynamics simulations and quantum chemical calculations. Experimentally, the cavitand has been shown to bind a number of amines and accelerate the methylation reaction by more than four orders of magnitude for some of them. Eight different amines are considered in the present study, and the geometries and energies of their binding to the cavitand are first characterized and analyzed. Next, the methyl transfer reactions are investigated and the calculated barriers are found to be in generally good agreement with experimental results. In particular, the experimentally-observed rate acceleration in the cavitand as compared to the solution reaction is well reproduced by the calculations. The origins of this rate acceleration are analyzed by computational modifications made to the structure of the cavitand, and the role of the solvent is discussed.

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.

Dynamics of organophosphate guest encapsulation in heteroleptic coordination cages

Heteroleptic coordination cages allow the design of different host structures that can bind guest molecules within their cavities. In a previous work, the energetics of organophosphate encapsulation in palladium(II)-based heteroleptic coordination cages that differ in terms of their ability to form hydrogen bonds have been investigated [Platzek et al., Endohedrally Functionalized Heteroleptic Coordination Cages for Phosphate Ester Binding, Angew. Chem., Int. Ed. 2022, 61, e2022093]. The present work focuses on the dynamics of this system. Dynamic information is obtained through the application of a Markov state model (MSM) to unbiased multi-microsecond atomistic molecular dynamics simulations of guest binding and release. The MSM reveals that both the bound state and the binding/unbinding pathways are highly dynamic, with different types of interactions mediating the binding of the diphenylphosphate guest. Thus, the simulations highlight the dynamic nature of the nanoconfinement in the host-guest systems, with possible implications for the use of such coordination cages as catalysts.

Guest Encapsulation Alters the Thermodynamic Landscape of a Coordination Host

The architecture of self-assembled host molecules can profoundly affect the properties of the encapsulated guests. For example, a rigid cage with small windows can efficiently protect its contents from the environment; in contrast, tube-shaped, flexible hosts with large openings and an easily accessible cavity are ideally suited for catalysis. Here, we report a "Janus" nature of a Pd6L4 coordination host previously reported to exist exclusively as a tube isomer (T). We show that upon encapsulating various tetrahedrally shaped guests, T can reconfigure into a cage-shaped host (C) in quantitative yield. Extracting the guest affords empty C, which is metastable and spontaneously relaxes to T, and the T⇄C interconversion can be repeated for multiple cycles. Reversible toggling between two vastly different isomers paves the way toward controlling functional properties of coordination hosts "on demand".

A porphyrin-based molecular cage guided by designed local-electric field is highly selective and efficient

The present work outlines a general methodology for designing efficient catalytic machineries that can easily be tweaked to meet the demands of the target reactions. This work utilizes a principle of the designed local electric field (LEF) as the driver for an efficient catalyst. It is demonstrated that by tweaking the LEF, we can catalyze the desired hydroxylation products with enantioselectivity that can be changed at will. Using computation tools, we caged a synthetic analog of heme porphyrin (HM1) and investigated the pharmaceutically relevant conversion of tetralin to tetralol, inside the modified supramolecular cage. The QM/MM calculations demonstrate a resulting catalytic efficiency with virtually absolute R-selectivity for the tetralin hydroxylation. Our calculations show that the LEF of the supramolecular cage and HM1 exert a strong electric field along the Fe-O reaction axis, which is the main driving force for enhanced reactivity. At the same time, the supramolecular cage applies a lateral LEF that regulates the enantioselectivity. We further demonstrate that swapping the charged/polar substitution in the supramolecular cage switches the lateral LEF which changes the enantioselectivity of hydroxylation from R to S.

PyConSolv: A Python Package for Conformer Generation of (Metal-Containing) Systems in Explicit Solvent

We introduce PyConSolv, a freely available Python package that automates the generation of conformers of metal- and nonmetal-containing complexes in explicit solvent, through classical molecular dynamics simulations. Using a streamlined workflow and interfacing with widely used computational chemistry software, PyConSolv is an all-in-one tool for the generation of conformers in any solvent. Input requirements are minimal; only the geometry of the structure and the desired solvent in xyz (XMOL) format are needed. The package can also account for charged systems, by including arbitrary counterions in the simulation. A bonded model parametrization is performed automatically, utilizing AmberTools, ORCA, and Multiwfn software packages. PyConSolv provides a selection of preparametrized solvents and counterions for use in classical molecular dynamics simulations. We show the applicability of our package on a number of (transition-metal-containing) systems. The software is provided open source and free of charge.

Constructing Ultrastable Metallo-Cages via In Situ Deprotonation/Oxidation of Dynamic Supramolecular Assemblies

Coordination-driven self-assembly enables the spontaneous construction of metallo-supramolecules with high precision, facilitated by dynamic and reversible metal-ligand interactions. The dynamic nature of coordination, however, results in structural lability in many metallo-supramolecular assembly systems. Consequently, it remains a formidable challenge to achieve self-assembly reversibility and structural stability simultaneously in metallo-supramolecular systems. To tackle this issue, herein, we incorporate an acid-/base-responsive tridentate ligand into multitopic building blocks to precisely construct a series of metallo-supramolecular cages through coordination-driven self-assembly. These dynamic cagelike assemblies can be transformed to their static states through mild in situ deprotonation/oxidation, leading to ultrastable skeletons that can withstand high temperatures, metal ion chelators, and strong acid/base conditions. This in situ transformation provides a reliable and powerful approach to manipulate the kinetic features and stability of metallo-supramolecules and allows for modulation of encapsulation and release behaviors of metallo-cages when utilizing nanoscale quantum dots (QDs) as guest molecules.

Cavity Characterization in Supramolecular Cages

Confining molecular guests within artificial hosts has provided a major driving force in the rational design of supramolecular cages with tailored properties. Over the last 30 years, a set of design strategies have been developed that enabled the controlled synthesis of a myriad of cages. Recently, there has been a growing interest in involving in silico methods in this toolbox. Cavity shape and size are important parameters that can be easily accessed by inexpensive geometric algorithms. Although these algorithms are well developed for the detection of nonartificial cavities (e.g., enzymes), they are not routinely used for the rational design of supramolecular cages. In order to test the capabilities of this tool, we have evaluated the performance and characteristics of seven different cavity characterization software in the context of 22 analogues of well-known supramolecular cages. Among the tested software, KVFinder project and Fpocket proved to be the most software to characterize supramolecular cavities. With the results of this work, we aim to popularize this underused technique within the supramolecular community.

Theoretical and computational methodologies for understanding coordination self-assembly complexes

This perspective highlights three theoretical and computational methods to capture the coordination self-assembly processes at the molecular level: quantum chemical modeling, molecular dynamics, and reaction network analysis. These methods cover the different scales from the metal-ligand bond to a more global aspect, and approaches that are best suited to understand the coordination self-assembly from different perspectives are introduced. Theoretical and numerical researches based on these methods are not merely ways of interpreting the experimental studies but complementary to them.

Solvent Induced Conversion of a Self-Assembled Gyrobifastigium to a Barrel and Encapsulation of Zinc-Phthalocyanine within the Barrel for Enhanced Photodynamic Therapy

A rare gyrobifastigium architecture (GB) was constructed by self-assembly of a tetradentate donor (L) with Pd^II acceptor in DMSO. The GB was converted to its isomeric tetragonal barrel (MB) upon treatment with water. The hydrophobic cavity of MB has been explored for the encapsulation of zinc-phthalocyanine (ZnPc), which is an excellent photosensitizer for photodynamic therapy (PDT). However, the poor water-solubility and aggregation tendency are the main reasons for the suboptimal PDT performance of free ZnPc in the aqueous medium. Effective solubilization of ZnPc in an aqueous medium was achieved by encapsulating it in the cavity of MB. The inclusion complex (ZnPc⊂MB) showed enhanced singlet oxygen generation in water. Higher cellular uptake and anticancer activity of the ZnPc⊂MB compared to free ZnPc on HeLa cells indicate that encapsulation of ZnPc in an aqueous host is a potential strategy for enhancement of its PDT activity in water.

Effector Regulated Catalytic Cyclization of Alkynoic Acids Using Pt2 L4 Cages

Metabolic pathways are highly regulated by effector molecules that influences the rate of enzymatic reactions. Inspired by the catalytic regulation found in living cells, we report a Pt₂ L₄ cage of which the activity can be controlled by effectors that bind inside the cage. The cage shows catalytic activity in the lactonization of alkynoic acids, with the reaction rates dependent on the effector guest bound in the cage. Some effector guests enhance the rate of the lactonization by up to 19-fold, whereas one decreases it by 5-fold. When mixtures of specific substrates are used, both starting materials and products act as guests for the Pt₂ L₄ cage, enhancing its catalytic activity for one substrate while reducing its activity for the other. The reported regulatory behavior obtained by the addition of effector molecules paves the way to the development of more complex, metabolic-like catalyst systems.

Supramolecular Coordination Cages for Artificial Photosynthesis and Synthetic Photocatalysis

Because sunlight is the most abundant energy source on earth, it has huge potential for practical applications ranging from sustainable energy supply to light driven chemistry. From a chemical perspective, excited states generated by light make thermodynamically uphill reactions possible, which forms the basis for energy storage into fuels. In addition, with light, open-shell species can be generated which open up new reaction pathways in organic synthesis. Crucial are photosensitizers, which absorb light and transfer energy to substrates by various mechanisms, processes that highly depend on the distance between the molecules involved. Supramolecular coordination cages are well studied and synthetically accessible reaction vessels with single cavities for guest binding, ensuring close proximity of different components. Due to high modularity of their size, shape, and the nature of metal centers and ligands, cages are ideal platforms to exploit preorganization in photocatalysis. Herein we focus on the application of supramolecular cages for photocatalysis in artificial photosynthesis and in organic photo(redox) catalysis. Finally, a brief overview of immobilization strategies for supramolecular cages provides tools for implementing cages into devices. This review provides inspiration for future design of photocatalytic supramolecular host-guest systems and their application in producing solar fuels and complex organic molecules.

Efficient Computation of the Interaction Energies of Very Large Non-covalently Bound Complexes

We present a new benchmark set consisting of 16 large non-covalently bound systems (LNCI16) ranging from 380 up to 1988 atoms and featuring diverse interaction motives. Gas-phase interaction energies are calculated with various composite DFT, semi-empirical quantum mechanical (SQM), and force field (FF) methods and are evaluated using accurate DFT reference values. Of the employed QM methods, PBEh-3c proves to be the most robust for large systems with a relative mean absolute deviation (relMAD) of 8.5% with respect to the reference interaction energies. r2SCAN-3c yields an even smaller relMAD, at least for the subset of complexes for which the calculation could be converged, but is less robust for systems with smaller HOMO–LUMO gaps. The inclusion of Fock-exchange is therefore important for the description of very large non-covalent interaction (NCI) complexes in the gas phase. GFN2-xTB was found to be the best performer of the SQM methods with an excellent result of only 11.1% deviation. From the assessed force fields, GFN-FF and GAFF achieve the best accuracy. Considering their low computational costs, both can be recommended for routine calculations of very large NCI complexes, with GFN-FF being clearly superior in terms of general applicability. Hence, GFN-FF may be routinely applied in supramolecular synthesis planning.

1 Introduction

2 The LNCI16 Benchmark Set

3 Computational Details

4 Generation of Reference Values

5 Results and Discussion

6 Conclusions

Encapsulation Enhances the Catalytic Activity of C-N Coupling: Reaction Mechanism of a Cu(I)/Calix[8]arene Supramolecular Catalyst

Development of C-N coupling methodologies based on Earth-abundant metals is a promising strategy in homogeneous catalysis for sustainable processes. However, such systems suffer from deactivation and low catalytic activity. We here report that encapsulation of Cu(I) within the phenanthroyl-containing calix[8]arene derivative 1,5-(2,9-dimethyl-1,10-phenanthroyl)-2,3,4,6,7,8-hexamethyl-p-tert-butylcalix[8]arene (C8PhenMe6 ) significantly enhances C-N coupling activity up to 92 % yield in the reaction of aryl halides and aryl amines, with low catalyst loadings (2.5 % mol). A tailored multiscale computational protocol based on Molecular Dynamics simulations and DFT investigations revealed an oxidative addition/reductive elimination process of the supramolecular catalyst [Cu(C8PhenMe6)I]. The computational investigations uncovered the origins of the enhanced catalytic activity over its molecular analogues: Catalyst deactivation through dimerization is prevented, and product release facilitated. Capturing the dynamic profile of the macrocycle and the impact of non-covalent interactions on reactivity allows for the rationalization of the behavior of the flexible supramolecular catalysts employed.

Helicate versus Mesocate in Quadruple-Stranded Lanthanide Cages: A Computational Insight

To drive the synthesis of metallo-supramolecular assemblies (MSAs) and to fully exploit their functional properties, robust computational tools are crucial. The capability to model and to rationalize different parameters that can influence the outcome is mandatory. Here, we report a computational insight on the factors that can determine the relative stability of the supramolecular isomers helicate and mesocate in lanthanide-based quadruple-stranded assemblies. The considered MSAs have the general formula [Ln₂L₄]²⁻ and possess a cavity suitable to allocate guests. The analysis was focused on three different factors: the ligand rigidity and the steric hindrance, the presence of a guest inside the cavity, and the guest dimension. Three different quantum mechanical calculation set-ups (in vacuum, with the solvent, and with the solvent and the dispersion correction) were considered. Comparison between theoretical and experimental outcomes suggests that all calculations correctly estimated the most stable isomer, while the inclusion of the dispersion correction is mandatory to reproduce the geometrical parameters. General guidelines can be drawn: less rigid and less bulky is the ligand and less stable is the helicate, and the presence of a guest can strongly affect the isomerism leading to an inversion of the stability by increasing the guest size when the ligand is flexible.

Automated Rational Design of Metal-Organic Polyhedra

Metal-organic polyhedra (MOPs) are hybrid organic-inorganic nanomolecules, whose rational design depends on harmonious consideration of chemical complementarity and spatial compatibility between two or more types of chemical building units (CBUs). In this work, we apply knowledge engineering technology to automate the derivation of MOP formulations based on existing knowledge. For this purpose we have (i) curated relevant MOP and CBU data; (ii) developed an assembly model concept that embeds rules in the MOP construction; (iii) developed an OntoMOPs ontology that defines MOPs and their key properties; (iv) input agents that populate The World Avatar (TWA) knowledge graph; and (v) input agents that, using information from TWA, derive a list of new constructible MOPs. Our result provides rapid and automated instantiation of MOPs in TWA and unveils the immediate chemical space of known MOPs, thus shedding light on new MOP targets for future investigations.